Groundwater is the world’s largest freshwater resource and is critically important for irrigation, and hence for global food security1,2,3. Already, unsustainable groundwater pumping exceeds recharge from precipitation and rivers4, leading to substantial drops in the levels of groundwater and losses of groundwater from its storage, especially in intensively irrigated regions5,6,7. When groundwater levels drop, discharges from groundwater to streams decline, reverse in direction or even stop completely, thereby decreasing streamflow, with potentially devastating effects on aquatic ecosystems. Here we link declines in the levels of groundwater that result from groundwater pumping to decreases in streamflow globally, and estimate where and when environmentally critical streamflows—which are required to maintain healthy ecosystems—will no longer be sustained. We estimate that, by 2050, environmental flow limits will be reached for approximately 42 to 79 per cent of the watersheds in which there is groundwater pumping worldwide, and that this will generally occur before substantial losses in groundwater storage are experienced. Only a small decline in groundwater level is needed to affect streamflow, making our estimates uncertain for streams near a transition to reversed groundwater discharge. However, for many areas, groundwater pumping rates are high and environmental flow limits are known to be severely exceeded. Compared to surface-water use, the effects of groundwater pumping are markedly delayed. Our results thus reveal the current and future environmental legacy of groundwater use.
This is a preview of subscription content
Subscription info for Chinese customers
We have a dedicated website for our Chinese customers. Please go to naturechina.com to subscribe to this journal.
Rent or Buy article
Get time limited or full article access on ReadCube.
All prices are NET prices.
All data needed to evaluate the conclusions in the paper are presented in the paper (Figs. 1, 2) and are available through the University of Victoria, https://doi.org/10.5683/SP2/D7I7CC. Additional model outputs (as part of the sensitivity analysis and model evaluation presented in the Extended Data) are prohibitively large to store in a repository but are available from the corresponding author on reasonable request.
The model code used to run the global-scale surface water—groundwater model is provided through a GitHub repository, https://github.com/UU-Hydro/PCR-GLOBWB_model/tree/develop/modflow.
Giordano, M. Global groundwater? Issues and solutions. Annu. Rev. Environ. Resour. 34, 153–178 (2009).
Döll, P. & Siebert, S. Global modelling of irrigation water requirements. Wat. Resour. Res. 38, 8-1–8-10 (2002).
Siebert, S. & Döll, P. Quantifying blue and green virtual water contents in global crop production as well as potential production losses without irrigation. J. Hydrol. 384, 198–217 (2010).
Gleeson, T. & Richter, B. How much groundwater can we pump and protect environmental flow through time? Presumptive standards for conjunctive management of aquifers and rivers. River Res. Appl. 34, 83–92 (2018).
Rodell, M., Velicogna, I. & Famigletti, J. S. Satellite-based estimates of groundwater depletion in India. Nature 460, 999–1002 (2009).
Aeschbach-Hertig, W. & Gleeson, T. Regional strategies for the accelerating global problem of groundwater depletion. Nat. Geosci. 5, 853–861 (2012).
Wada, Y., Wisser, D. & Bierkens, M. F. P. Global modelling of withdrawal, allocation and consumptive use of surface water and groundwater resources. Earth Syst. Dynam. 5, 15–40 (2014).
Taylor, R. G. et al. Groundwater and climate change. Nat. Clim. Chang. 3, 322–329 (2013).
Konikow, L. F. & Kendy, E. Groundwater depletion: a global problem. Hydrogeol. J. 13, 317–320 (2005).
Lawford, R. et al. Basin perspectives on the water–energy–food security nexus. Curr. Opin. Environ. Sustain. 5, 607–616 (2013).
Postel, S. L. Entering an era of water scarcity: the challenges ahead. Ecol. Appl. 10, 941–948 (2000).
de Graaf, I. E. M., van Beek, L. P. H., Wada, Y. & Bierkens, M. F. P. Dynamic attribution of global water demand to surface water and groundwater resources: effects of abstractions and return flows on river discharges. Adv. Water Resour. 64, 21–33 (2014).
van Beek, L. P. H., Wada, Y. & Bierkens, M. F. P. Global monthly water stress: 1. Water balance and water availability. Wat. Resour. Res. 47, W07517 (2011).
de Graaf, I. E. M. et al. A global-scale two-layer transient groundwater model: development and application to groundwater depletion. Adv. Water Resour. 102, 53–67 (2017).
Sutanudjaja, E. H., et al. PCR-GLOBWB 2: a 5 arc-minute global hydrological and water resources model. Geosci. Model Dev. 11, 2429–2453 (2018).
Ficklin, D. L., Robeson, S. M. & Knouft, J. H. Impacts of recent climate change on trends in baseflow and stormflow in United States watersheds. Geophys. Res. Lett. 43, 5079–5088 (2016).
Weedon, G. P. et al. Creation of the WATCH forcing data and its use to assess global and regional reference crop evaporation over land during the twentieth century. J. Hydrometeorol. 12, 823–848 (2011).
Hempel, S., Frieler, K., Warszawski, L., Schewe, J. & Piontek, F. A trend-reserving bias correction—the ISI-MIP approach. Earth Syst. Dynam. 4, 219–236 (2013).
World Water Assessment Programme The United Nations World Water Development Report 4: Managing Water under Uncertainty and Risk. Report No. 978-92-3-104235-5, 407 (UNESCO, 2012).
Döll, P., Fiedler, K. & Zhang, J. Global-scale analysis of river flow alterations due to water withdrawals and reservoirs. Hydrol. Earth Syst. Sci. 13, 2413–2432 (2009).
India Central Ground Water Board Ground Water Yearbook 2013–14 http://cgwb.gov.in/Documents/Ground%20Water%20Year%20Book%202013-14.pdf (2014).
Perrone, D. & Jasechko, S. Dry groundwater wells in the western United States. Environ. Res. Lett. 12, 104002 (2017).
Lehner, B., Verdin, K. & Jarvis, A. New global hydrography derived from spaceborne elevation data. Eos 89, 93–94 (2008).
Winter, T. C., Harvey, J. W., Franke, O. L. & Alley, W. M. Ground Water and Surface Water: A Single Resource. Circular 1139 (United States Geological Survey, 1998).
Alley, W. M., Reilly, T. E. & Franke, O. L. Sustainability of Groundwater Resources. Circular 1186 (United States Geological Survey, 1999).
Konikow, L. F. & Leake, S. A. Depletion and capture: revisiting “the source of water derived from wells”. Ground Water 52, 100–111 (2014).
Wada, Y. et al. Global depletion of groundwater resources. Geophys. Res. Lett. 37, L20402 (2010).
Hartmann, J. & Moosdorf, N. The new global lithological map database glim: a representation of rock properties at the earth surface. Geochem. Geophy. Geosy. 13, Q12004 (2012).
Gleeson, T., Moosdorf, N., Hartmann, J. & van Beek, L. P. H. A glimpse beneath earth’s surface: GLobal HYdrogeological MaPS (GLHYMPS) of permeability and porosity. Geophys. Res. Lett. 41, 3891–3898 (2014).
de Graaf, I. E. M., van Beek, L. P. H., Sutanudjaja, E. H. & Bierkens, M. F. P. A high-resolution global-scale groundwater model. Hydrol. Earth Syst. Sci. 19, 823–837 (2015).
Sutanudjaja, E. H., van Beek, L. P. H., De Jong, S. M., Van Geer, F. C. & Bierkens, M. F. P. Calibrating a large-extent high-resolution coupled groundwater–land surface model using soil moisture and discharge data. Wat. Resour. Res. 50, 687–705 (2014).
Verdin, K. L. & Greenlee, S. K. Development of continental scale digital elevation models and extraction of hydrographic features. In Proc. Third Intl Conf., Workshop on Integrating GIS and Environmental Modeling, 1996 (National Center for Geographic Information and Analysis, 1996).
Fan, Y., Li, H. & Miguez-Macho, G. Global patterns of groundwater table depth. Science 339, 940–943 (2013).
McGuire, V. L. Water Level Changes in the High Plain Aquifers, Predevelopment to 2007, 2005–06, and 2006–2007. Scientific Investigations report 2009-5019 (United States Geological Survey, 2009).
Williamson, A. K., Prudic, D. E. & Swain, L. A. Ground-Water Flow in the Central Valley, California. Professional paper 1401-D (United States Geological Survey, 1989).
Konikow, L. F. Contribution of global groundwater depletion since 1900 to sea-level rise. Geophys. Res. Lett. 38, L17401 (2011).
Richey, A. S. et al. Uncertainty in global groundwater storage estimates in a total groundwater stress framework. Wat. Resour. Res. 51, 5198–5216 (2015).
MacDonald, A. M. et al. Groundwater quality and depletion in the Indo–Gangetic basin mapped from in situ observations. Nat. Geosci. 9, 762–766 (2016).
Republican River Compact Administration, Groundwater Model Update 2001–2003 http://www.republicanrivercompact.org/ (2005).
Van Loon, A. F. & Van Lanen, H. A. J. Making the distinction between water scarcity and drought using an observation-modelling framework. Wat. Resour. Res. 49, 1483–1502 (2013).
Ahner, B. W. Towards a Representative Scale of the Groundwater Footprint Based on the Ecohydrological Impacts of Groundwater Abstraction, MSc thesis, Georg-August-Universität Göttingen, Göttingen, Germany. (2018)
We thank S. Rangecroft, Y. Wada, E. Luijendijk, B. Ahner and S. de Kock for generously sharing data that supported the analysis done in this study. This study’s calculations were computed on the Dutch national supercomputer Cartesius with the support of SURFsara. This research was funded by the Netherlands Organization for Scientific Research (NWO) under the programme ‘Planetary boundaries of the global fresh water cycle’.
The authors declare no competing interests.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Peer review information Nature thanks Laura Foglia, Mary C. Hill and Timothy R. Green for their contribution to the peer review of this work.
Extended data figures and tables
a, Cumulative groundwater depletion trend (in km3) since 1960 for different climate scenarios. After 2010, the estimate assumes a ‘business-as-usual’ scenario for water demands and uses three GCMs—HadGEM2-HS, GFDL-ESM2M and MIROC-ESM-CHEM—using RCP 8.5. b, c, The first time the environmental flow limit is or will be reached (by year) averaged over watersheds using the sub-watershed level of hydroBASINS28 for the climate scenarios GFDL-ESM2M (a) and MIROC-ESM-CHEM (b) using RCP 8.5.
Extended Data Fig. 2 Gridded estimates of cumulative groundwater water depletion (in m3 per m2) for 1960–2099.
Four major heavily pumped aquifers are magnified. Aquifer magnifications are from WHYMAP, BGS/UNESCO.
Extended Data Fig. 3 Model evaluation of simulated groundwater head changes owing to pumping compared to observations.
Three of the largest, intensively pumped, and best monitored alluvial aquifer systems of the world are shown: the High Plains aquifer (left) and Central Valley aquifer, USA, (middle) and the Upper Ganges and Indus basin, India (right). a, Observed data and published maps. b, This study’s estimates. A comparison between a and b shows that the model results matches the observations well. Nonlinear colour scales are used. The Ganges basin map is from a previous work38.
Extended Data Fig. 4 Model evaluation of the estimated first time that the environmental flow limits are reached.
a, Observed versus simulated first time that the environmental flow limits are reached (x and y axes in years) for several groundwater-pumping-impacted catchments. b, Locations of the studied catchments are indicated on the map; the three dots in Kansas represent nine sub-catchments used in the analysis.
Extended Data Fig. 5 Distribution of surface water–groundwater interaction classes for north America as simulated using the physically based global-scale GSGM.
The figure distinguishes four classes: gaining streams, losing streams, intermittently disconnected streams and continuously disconnected streams. Results are shown for the month of July (generally the driest month of the year, on the basis of monthly discharge) and represent five-year moving averages of river drainage and groundwater levels. A stream is classified as a continuously disconnected stream if the stream is disconnected for at least two years in a row over the moving average window of five years. a, Spatial maps of surface water–groundwater interaction for July 1964 and 2004 under natural conditions and including human water withdrawal. b, Temporal variation of the total area covered by each surface water–groundwater interaction class. The thinner continuous lines are yearly values, thicker dashed and dotted lines are the five-year averages.
Extended Data Fig. 6 Evaluation of observed versus simulated water table averaged for sub-watersheds.
a, Scatter plot; the red line shows the 1:1 slope. b, Histogram of relative residuals, calculated as (observed − simulated)/observed. c, Global map of relative residuals. All watersheds with no available data are mapped in grey.
Extended Data Fig. 7 Model evaluation of simulated groundwater table depth to well observations for two well monitored aquifer systems in the USA.
a, b, Location of the wells used for this analysis within the Central Valley aquifer system (a) the High Plains aquifer system (b). c–h, Averaged water table depths (‘wtd’) (c, d); standard deviation (‘std’) of monthly wtd (e, f) and head drops (g, h) were estimated (‘obs avg’) and compared to simulated results (‘sim avg’). In each plot, the values per well in the watershed are given in grey (well data) and the statistics show the wide spread in observations. The value n indicates the number of wells within the watershed. The watersheds are numbered from north to south over both aquifers, indicated by the watershed number (the exact location of the watersheds is not relevant).
a, b, Comparison of groundwater head declines (a) and groundwater discharge (GD; b) at different spatial levels, simulated by a calibrated regional-scale model (‘calb.’); results modified from previous work41 using the Republican River Groundwater Model39 and the global-scale surface water–groundwater model (‘glob.’). c, The full basin covers the entire Republican River basin, which is situated in the central-north of the High Plains aquifer, USA. The larger sub-basins are the level 2 and 3 regional basins that consist of more level 1 basins. In b, the solid lines present the simulated groundwater discharge, the dashed lines present the trends of the groundwater discharge. Groundwater discharge trends simulated by both models are comparable.
a, Table showing the different parameter settings used; varying the sub-surface conductivity, the river depth and the river conductance by decreasing or increasing the settings of the baseline (‘bl’) run. The baseline run uses the average parameter settings. b, Frequency plot of the first time the environmental flow limits are reached under different parameter settings for 1960–2010. The fractional increase or decrease of estimated environmental limits compared to the baseline run is given in the fourth column of the table in a. The runs with the smallest and largest limits are indicated in bold in the table and presented in blue and yellow, respectively, in the graph.
a, Table giving the different criteria of the Q90 windows and consecutive years used to estimate the environmental flow limits. b, c, Histograms of the limits reached, estimated using the Q90 over five years (b) and the Q90 over ten years (c). The fractional increase or decrease of the estimated environmental limits compared to the baseline run (Q90 over five years, for two consecutive years) is given in the third column of the table in a. Difference in the estimated environmental flow limits are only limited when using different criteria.
About this article
Cite this article
de Graaf, I.E.M., Gleeson, T., (Rens) van Beek, L.P.H. et al. Environmental flow limits to global groundwater pumping. Nature 574, 90–94 (2019). https://doi.org/10.1038/s41586-019-1594-4
Nature Communications (2021)
Nature Climate Change (2021)
Scientific Reports (2021)
Nature Communications (2021)