Identifying groundwater recharge and discharge zones using geostatistical simulation of hydraulic head and its derivatives

Identifying groundwater flow directions and the locations of recharge and discharge areas is critical for effective groundwater management. Groundwater flow directions, the concavity and the locations of extrema (i.e., minima and maxima) can be assessed using the first and second derivatives of the hydraulic head surface. We developed a geostatistical method to jointly simulate hydraulic head and its first and second derivatives using sequential Gaussian simulation. The derivative values were used to identify regional groundwater flow directions, and the second derivative test was used to probabilistically map the concavity and the locations of extrema in the hydraulic head surface. By comparing the mapped concavity and extrema to known features, it was possible to attribute areas of recharge and discharge to physical features of the system, such as rivers, lakes and geological outcrops. This was applied to Triassic aquifers in the Galilee Basin (Queensland, Australia) to delineate the likely recharge and discharge areas. This provided an objective assessment of likely recharge and discharge zones and their uncertainty, which is an important addition to a region where the hydrogeology has been the subject of much conjecture.