Estimating groundwater recharge and its associated uncertainty: Use of regression kriging and the chloride mass balance method

The chloride mass balance method of estimating groundwater recharge is generally applied at a point scale but water resources management generally requires information at a regional scale. To estimate groundwater recharge regionally previously used upscaling methods have ranged from simple averaging to empirical relationships and geostatistical interpolation. This study combined the best components of these methods by using regression kriging: in data-sparse areas the recharge is upscaled using global regression equations with gridded rainfall and surface geology as covariates, while in data-dense areas the kriging of the regression equation residuals ensures that the upscaled recharge estimates respect the point estimates of recharge. The uncertainty in the recharge estimates was quantified using 1000 stochastic replicates of the chloride deposition of rainfall, the chloride exported in runoff, the chloride concentration of the groundwater and the regression equations used to perform the upscaling. This study focused on the coal bearing Sydney, Gunnedah, Gloucester and Surat Basins of Eastern New South Wales (Australia). Historically, groundwater recharge to the Permian units of these basins has received little attention due to their low yields of poor quality water. The increased potential extraction of regional groundwater due to coal development has demanded a greater understanding of the water balance. We found that recharge is highest in the younger productive aquifers (up to 20% of long-term mean annual rainfall) and lowest in the Permian units (∼1% of long-term mean annual rainfall). The magnitude of the uncertainty is often close to the magnitude of the median recharge estimate. The method developed here for upscaling the point estimates of recharge using the chloride mass balance provides robust estimates of recharge and its associated uncertainty. The method is applicable to any regional study with variable density input data.

[1]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[2]  Lina Schelin,et al.  Kriging Prediction Intervals Based on Semiparametric Bootstrap , 2010 .

[3]  K. Ivkovic,et al.  A top–down approach to characterise aquifer–river interaction processes , 2009 .

[4]  G. B. Allison,et al.  Recharge in karst and dune elements of a semi-arid landscape as indicated by natural isotopes and chloride , 1985 .

[5]  G. Heuvelink,et al.  A generic framework for spatial prediction of soil variables based on regression-kriging , 2004 .

[6]  O. Batelaan,et al.  Spatial variability of chloride deposition in a vegetated coastal area: Implications for groundwater recharge estimation , 2014 .

[7]  S. Schiff,et al.  A validation of the 3H/3He method for determining groundwater recharge , 1993 .

[8]  B. Scanlon,et al.  Estimating Groundwater Recharge , 2010 .

[9]  R. Crosbie,et al.  Deep drainage estimates using multiple linear regression with percent clay content and rainfall , 2012 .

[10]  E. Custodio,et al.  Estimating natural recharge by means of chloride mass balance in a volcanic aquifer: northeastern Gran Canaria (Canary Islands, Spain) , 2015 .

[11]  S. T. Buckland,et al.  An Introduction to the Bootstrap. , 1994 .

[12]  T. McVicar,et al.  Observational evidence from two mountainous regions that near‐surface wind speeds are declining more rapidly at higher elevations than lower elevations: 1960–2006 , 2010 .

[13]  Arthur C. Petersen,et al.  A method for the analysis of assumptions in model-based environmental assessments , 2011, Environ. Model. Softw..

[14]  J. Williams,et al.  Groundwater residence time in a dissected and weathered sandstone plateau: Kulnura–Mangrove Mountain aquifer, NSW, Australia , 2014 .

[15]  W. Edmunds,et al.  Published online in Wiley InterScience (www.interscience.wiley.com) DOI: 10.1002/hyp.6335 Global synthesis of groundwater recharge in semiarid andaridregions , 2022 .

[16]  M. Bakker,et al.  "Groundwater recharge: processes and quantification": Preface , 2013 .

[17]  E. Eriksson.  Principles and Applications of Hydrochemistry , 1985 .

[18]  Jin Teng,et al.  The Australian Water Resource Assessment Modelling System (AWRA) , 2013 .

[19]  W. W. Wood,et al.  Chemical and Isotopic Methods for Quantifying Ground‐Water Recharge in a Regional, Semiarid Environment , 1995 .

[20]  Dominique Fasbender,et al.  Bayesian data fusion in a spatial prediction context: a general formulation , 2007 .

[21]  C. Nichol,et al.  Hydrogeochemical and isotopic constraints on the origins of dryland salinity, Murray Basin, Victoria, Australia , 2004 .

[22]  D. Jones,et al.  High-quality spatial climate data-sets for Australia , 2009 .

[23]  Glen R. Walker,et al.  Factors affecting groundwater recharge following clearing in the south western Murray Basin , 1994 .

[24]  D. Mcphail,et al.  Chloride and bromide sources in water: Quantitative model use and uncertainty , 2017 .

[25]  A. Herczeg,et al.  Origin of dissolved salts in a large, semi-arid groundwater system: Murray Basin, Australia , 2001 .

[26]  E. Custodio,et al.  Spatial average aquifer recharge through atmospheric chloride mass balance and its uncertainty in continental Spain , 2014 .

[27]  D. McWhorter,et al.  Groundwater recharge assessment in an upland sandstone aquifer of southern California , 2016 .

[28]  C. Simmons,et al.  Groundwater recharge to a sedimentary aquifer in the topographically closed Uley South Basin, South Australia , 2012, Hydrogeology Journal.

[29]  Trevor Pickett,et al.  An assessment of the climate change impacts on groundwater recharge at a continental scale using a probabilistic approach with an ensemble of GCMs , 2013, Climatic Change.

[30]  E. Custodio,et al.  Natural uncertainty of spatial average aquifer recharge through atmospheric chloride mass balance in continental Spain , 2015 .

[31]  I. Cartwright,et al.  Constraining modern and historical recharge from bore hydrographs, 3H, 14C, and chloride concentrations: Applications to dual-porosity aquifers in dryland salinity areas, Murray Basin, Australia , 2007 .

[32]  B. Smerdon,et al.  Hydrostratigraphy, hydrogeology and system conceptualisation of the Great Artesian Basin. A technical report to the Australian Government from the CSIRO Great Artesian Basin Water Resource Assessment , 2012 .

[33]  J. Ball,et al.  Australian Rainfall and Runoff: A Guide to Flood Estimation , 2016 .

[34]  Brett A. Bryan,et al.  Mapping agriculture's impact by combining farm management handbooks, life-cycle assessment and search engine science , 2016, Environ. Model. Softw..

[35]  E. Eriksson,et al.  Chloride concentration in groundwater, recharge rate and rate of deposition of chloride in the Israel Coastal Plain , 1969 .

[36]  R. Crosbie,et al.  Mapping the spatial distribution of chloride deposition across Australia , 2018 .

[37]  Warren W. Wood,et al.  Chloride mass-balance method for estimating ground water recharge in arid areas: examples from western Saudi Arabia , 1996 .

[38]  Rodger B. Grayson,et al.  Towards a framework for predicting impacts of land-use on recharge: 1. A review of recharge studies in Australia , 2002 .

[39]  S. Landsberger,et al.  Characterization and Identification of Na‐Cl Sources in Ground Water , 2006, Ground water.

[40]  P. Cook,et al.  Spatial variability of groundwater recharge in a semiarid region , 1989 .

[41]  Jack P. C. Kleijnen,et al.  The correct Kriging variance estimated by bootstrapping , 2006, J. Oper. Res. Soc..

[42]  E. Custodio,et al.  Integrating soil water and tracer balances, numerical modelling and GIS tools to estimate regional groundwater recharge: Application to the Alcadozo Aquifer System (SE Spain). , 2016, The Science of the total environment.

[43]  R. Reedy,et al.  Groundwater depletion and sustainability of irrigation in the US High Plains and Central Valley , 2012, Proceedings of the National Academy of Sciences.

[44]  M. Andersen,et al.  Stream-aquifer interactions in the Maules Creek catchment, Namoi Valley, New South Wales, Australia , 2009 .

[45]  P. Binning,et al.  A time series approach to inferring groundwater recharge using the water table fluctuation method , 2005 .

[46]  H. Hofmann,et al.  Using hydrogeochemistry to understand inter-aquifer mixing in the on-shore part of the Gippsland Basin, southeast Australia , 2013 .

[47]  N. Laftouhi,et al.  Estimate of regional groundwater recharge rate in the Central Haouz Plain, Morocco, using the chloride mass balance method and a geographical information system , 2017, Applied Water Science.

[48]  J. Böhlke,et al.  Groundwater recharge and agricultural contamination , 2002 .

[49]  R. Doble,et al.  Review: Current and emerging methods for catchment-scale modelling of recharge and evapotranspiration from shallow groundwater , 2017, Hydrogeology Journal.

[50]  W. Wood.  Use and Misuse of the Chloride‐Mass Balance Method in Estimating Ground Water Recharge , 1999 .

[51]  Cuan Petheram,et al.  Can the dataset of field based recharge estimates in Australia be used to predict recharge in data-poor areas? , 2010 .

[52]  E. Custodio,et al.  Using the Cl/Br ratio as a tracer to identify the origin of salinity in aquifers in Spain and Portugal , 2008 .

[53]  M. Norman,et al.  Salinisation processes in a sub-catchment of Wybong Creek, Hunter Valley, Australia , 2011 .

[54]  Fernando Bação,et al.  Exploratory data analysis and clustering of multivariate spatial hydrogeological data by means of GEO3DSOM, a variant of Kohonen's Self-Organizing Map , 2006 .

[55]  Luk J M Peeters,et al.  Assumption Hunting in Groundwater Modeling: Find Assumptions Before They Find You , 2017, Ground water.

[56]  Edzer J. Pebesma,et al.  Multivariable geostatistics in S: the gstat package , 2004, Comput. Geosci..

[57]  R. Reedy,et al.  Impact of land use and land cover change on groundwater recharge and quality in the southwestern US , 2005 .

[58]  P. Cook,et al.  Spatial and Temporal Variability of Ground Water Recharge in Central Australia: A Tracer Approach , 2002, Ground water.

[59]  Alain Dassargues,et al.  Bayesian data fusion for water table interpolation: Incorporating a hydrogeological conceptual model in kriging , 2010 .

[60]  Mike Grundy,et al.  Soil and landscape grid of Australia. , 2015 .

[61]  L. Fifield,et al.  Differences in groundwater and chloride residence times in saline groundwater: The Barwon River Catchment of Southeast Australia , 2017 .

[62]  M. Keywood,et al.  The accession of chloride to the western half of the Australian continent , 1997 .

[63]  M. Hughes,et al.  The use of environmental chloride and tritium to estimate total recharge to an unconfined aquifer , 1978 .

[64]  B. Scanlon,et al.  Choosing appropriate techniques for quantifying groundwater recharge , 2002 .

[65]  Brian Smerdon,et al.  Water resource assessment for the Surat region. A report to the Australian Government from the CSIRO Great Artesian Basin Water Resource Assessment , 2012 .

[66]  Robert B. Jackson,et al.  A Global Analysis of Groundwater Recharge for Vegetation, Climate, and Soils , 2012 .