Using Satellite-Based Terrestrial Water Storage Data: A Review

Land water storage plays a key role for the Earth's climate, natural ecosystems, and human activities. Since the launch of the first Gravity Recovery and Climate Experiment (GRACE) mission in 2002, spaceborne observations of changes in terrestrial water storage (TWS) have provided a unique, global perspective on natural and human-induced changes in freshwater resources. Even though they have become much used within the broader Earth system science community, space-based TWS datasets still incorporate important and case-specific limitations which may not always be clear to users not familiar with the underlying processing algorithms. Here, we provide an accessible and illustrated overview of the measurement concept, of the main available data products, and of some frequently encountered technical terms and concepts. We summarize concrete recommendations on how to use TWS data in combination with other hydrological or climatological datasets, and guidance on how to avoid possible pitfalls. Finally, we provide an overview of some of the main applications of GRACE TWS data in the fields of hydrology and climate science. This review is written with the intention of supporting future research and facilitating the use of satellite-based terrestrial water storage datasets in interdisciplinary contexts.

Groundwater storage change detection from in situ and GRACE-based estimates in major river basins across India

India has been the subject of many recent groundwater studies due to the rapid depletion of groundwater in large parts of the country. However, few if any of these studies have examined groundwater storage conditions in all of India's river basins individually. Herein we assess groundwater storage changes in all 22 of India's major river basins using in situ data from 3420 observation locations for the period 2003-2014. One-month and 12-month standardized precipitation index measures (SPI-1 and SPI-12) indicate fluctuations in the long-term pattern. The Ganges and Brahmaputra basins experienced long-term decreasing trends in precipitation in both 1961-2014 and the study period, 2003-2014. Indeterminate or increasing precipitation trends occurred in other basins. Satellite-based and in situ groundwater storage time series exhibited similar patterns, with increases in most of the basins. However, diminishing groundwater storage (at rates of >0.4 km³/year) was revealed in the Ganges-Brahmaputra river basin based on in situ observations, which is particularly important due to its agricultural productivity.

GRACE improves seasonal groundwater forecast initialization over the U.S

We evaluate the impact of Gravity Recovery and Climate Experiment data assimilation (GRACE-DA) on seasonal hydrological forecast initialization over the U.S., focusing on groundwater storage. GRACE-based terrestrial water storage (TWS) estimates are assimilated into a land surface model for the 2003-2016 period. Three-month hindcast (i.e., forecast of past events) simulations are initialized using states from the reference (no data assimilation) and GRACE-DA runs. Differences between the two initial hydrological condition (IHC) sets are evaluated for two forecast techniques at 305 wells where depth-to-water-table measurements are available. Results show that using GRACE-DA-based IHC improves seasonal groundwater forecast performance in terms of both RMSE and correlation. While most regions show improvement, degradation is common in the High Plains, where withdrawals for irrigation practices affect groundwater variability more strongly than the weather variability, which demonstrates the need for simulating such activities. These findings contribute to recent efforts towards an improved U.S. drought monitor and forecast system.

NCA-LDAS: Overview and Analysis of Hydrologic Trends for the National Climate Assessment

Terrestrial hydrologic trends over the conterminous United States are estimated for 1980-2015 using the National Climate Assessment Land Data Assimilation System (NCA-LDAS) reanalysis. NCA-LDAS employs the uncoupled Noah version 3.3 land surface model at 0.125°× 1258° forced with NLDAS-2 meteorology, rescaled Climate Prediction Center precipitation, and assimilated satellite-based soil moisture, snow depth, and irrigation products. Mean annual trends are reported using the nonparametric Mann-Kendall test at p < 0.1 significance. Results illustrate the interrelationship between regional gradients in forcing trends and trends in other land energy and water stores and fluxes. Mean precipitation trends range from +3 to +9 mm yr-1 in the upper Great Plains and Northeast to -1 to -9 mm yr-1 in the West and South, net radiation flux trends range from 10.05 to 10.20 W m-2 yr-1 in the East to -0.05 to -0.20 W m-2 yr-1 in the West, and U.S.-wide temperature trends average about +0.03 K yr-1. Trends in soil moisture, snow cover, latent and sensible heat fluxes, and runoff are consistent with forcings, contributing to increasing evaporative fraction trends from west to east. Evaluation of NCA-LDAS trends compared to independent data indicates mixed results. The RMSE of U.S.-wide trends in number of snow cover days improved from 3.13 to 2.89 days yr-1 while trend detection increased 11%. Trends in latent heat flux were hardly affected, with RMSE decreasing only from 0.17 to 0.16 W m-2 yr-1, while trend detection increased 2%. NCA-LDAS runoff trends degraded significantly from 2.6 to 16.1 mm yr-1 while trend detection was unaffected. Analysis also indicated that NCA-LDAS exhibits relatively more skill in low precipitation station density areas, suggesting there are limits to the effectiveness of satellite data assimilation in densely gauged regions. Overall, NCA-LDAS demonstrates capability for quantifying physically consistent, U.S. hydrologic climate trends over the satellite era.

Contributions of GRACE to understanding climate change

Time-resolved satellite gravimetry has revolutionized understanding of mass transport in the Earth system. Since 2002, the Gravity Recovery and Climate Experiment (GRACE) has enabled monitoring of the terrestrial water cycle, ice sheet and glacier mass balance, sea level change and ocean bottom pressure variations and understanding responses to changes in the global climate system. Initially a pioneering experiment of geodesy, the time-variable observations have matured into reliable mass transport products, allowing assessment and forecast of a number of important climate trends and improve service applications such as the U.S. Drought Monitor. With the successful launch of the GRACE Follow-On mission, a multi decadal record of mass variability in the Earth system is within reach.

Earth observations and integrative models in support of food and water security

Global food production depends upon many factors that Earth observing satellites routinely measure about water, energy, weather, and ecosystems. Increasingly sophisticated, publicly-available satellite data products can improve efficiencies in resource management and provide earlier indication of environmental disruption. Satellite remote sensing provides a consistent, long-term record that can be used effectively to detect large-scale features over time, such as a developing drought. Accuracy and capabilities have increased along with the range of Earth observations and derived products that can support food security decisions with actionable information. This paper highlights major capabilities facilitated by satellite observations and physical models that have been developed and validated using remotely-sensed observations. Although we primarily focus on variables relevant to agriculture, we also include a brief description of the growing use of Earth observations in support of aquaculture and fisheries.

Author Correction: Emerging trends in global freshwater availability

In Fig. 2 of this Analysis, the tick-mark labels on the colour bars in the second and third images from the top were inadvertently swapped. In addition, the citation at the end of the sentence, "On a monthly basis GRACE can resolve TWS changes with sufficient accuracy over scales that range from approximately 200,000 km² at low latitudes to about 90,000 km² near the poles" should be to ref. ⁴ not ref. ¹. These errors have been corrected online.

Evaluation of Simulated Snow and Snowmelt Timing in the Community Land Model Using Satellite-based Products and Streamflow Observations

The purpose of this study was to evaluate snow and snowmelt simulated by version 4 of the Community Land Model (CLM4). We performed uncoupled CLM4 simulations, forced by Modem-Era Retrospective Analysis for Research and Applications Land-only (MERRA-Land) meteorological fields. GlobSnow snow cover fraction (SCF), snow water equivalent (SWE) and satellite-based passive microwave (PMW) snowmelt-off day of year (MoD) data were used to evaluate SCF, SWE, and snowmelt simulations. Simulated runoff was then fed into a river routing scheme and evaluation was performed at 408 snow-dominated catchments using gauge observations. CLM4 and GlobSnow snow cover extent showed a strong agreement, especially during the peak snow cover months. Overall there was a good correlation between simulated and observed SWE (correlation coefficient, R = 0.6). Simulated and observed SWE were similar over areas with relatively flat terrain and moderate forest density. The simulated MoD agreed (MoD differences (CLM4-PMW) = +/-7 days) with observations over 39.4% of the study domain. Snowmelt-off occurred earlier in the model compared to the observations over 39.5 % of the domain and later over 21.1% of the domain. Large differences of MoD were seen in the areas with complex terrain and dense forest cover. We also found that, although streamflow seasonal phase was accurately modeled (R=0.9), the peaks controlled by snowmelt were underestimated. Routed CLM4 streamflow tended to occur early (by 10 days on average).

Emerging trends in global freshwater availability

Freshwater availability is changing worldwide. Here we quantify 34 trends in terrestrial water storage observed by the Gravity Recovery and Climate Experiment (GRACE) satellites during 2002-2016 and categorize their drivers as natural interannual variability, unsustainable groundwater consumption, climate change or combinations thereof. Several of these trends had been lacking thorough investigation and attribution, including massive changes in northwestern China and the Okavango Delta. Others are consistent with climate model predictions. This observation-based assessment of how the world's water landscape is responding to human impacts and climate variations provides a blueprint for evaluating and predicting emerging threats to water and food security.

Emerging trends in global freshwater availability

Freshwater availability is changing worldwide. Here we quantify 34 trends in terrestrial water storage observed by the Gravity Recovery and Climate Experiment (GRACE) satellites during 2002-2016 and categorize their drivers as natural interannual variability, unsustainable groundwater consumption, climate change or combinations thereof. Several of these trends had been lacking thorough investigation and attribution, including massive changes in northwestern China and the Okavango Delta. Others are consistent with climate model predictions. This observation-based assessment of how the world's water landscape is responding to human impacts and climate variations provides a blueprint for evaluating and predicting emerging threats to water and food security.

The Value of Remotely Sensed Information: The Case of GRACE-Enhanced Drought Severity Index

A decision framework is developed for quantifying the economic value of information (VOI) from the Gravity Recovery and Climate Experiment (GRACE) satellite mission for drought monitoring, with a focus on the potential contributions of groundwater storage and soil moisture measurements from the GRACE Data Assimilation (GRACE-DA) System. The study consists of: (a) the development of a conceptual framework to evaluate the socioeconomic value of GRACE-DA as a contributing source of information to drought monitoring; (b) structured listening sessions to understand the needs of stakeholders who are affected by drought monitoring; (c) econometric analysis based on the conceptual framework that characterizes the contribution of GRACE-DA to the US Drought Monitor (USDM) in capturing the effects of drought on the agricultural sector; and (d) a demonstration of how the improved characterization of drought conditions may influence decisions made in a real-world drought disaster assistance program. Results show that GRACE-DA has the potential to lower the uncertainty associated with our understanding of drought, and that this improved understanding has the potential to change policy decisions that lead to tangible societal benefits.

The sensitivity of US wildfire occurrence to pre-season soil moisture conditions across ecosystems

It is generally accepted that year-to-year variability in moisture conditions and drought are linked with increased wildfire occurrence. However, quantifying the sensitivity of wildfire to surface moisture state at seasonal lead-times has been challenging due to the absence of a long soil moisture record with the appropriate coverage and spatial resolution for continental-scale analysis. Here we apply model simulations of surface soil moisture that numerically assimilate observations from NASA's Gravity Recovery and Climate Experiment (GRACE) mission with the US Forest Service's historical Fire-Occurrence Database over the contiguous United States. We quantify the relationships between pre-fire-season soil moisture and subsequent-year wildfire occurrence by land-cover type and produce annual probable wildfire occurrence and burned area maps at 0.25-degree resolution. Cross-validated results generally indicate a higher occurrence of smaller fires when months preceding fire season are wet, while larger fires are more frequent when soils are dry. This result is consistent with the concept of increased fuel accumulation under wet conditions in the pre-season. These results demonstrate the fundamental strength of the relationship between soil moisture and fire activity at long lead-times and are indicative of that relationship's utility for the future development of national-scale predictive capability.

How might recharge change under projected climate change in western US?

Although groundwater is a major resource of water in the western US, little research has been done on the impacts of climate change on groundwater storage and recharge in the West. Here we assess the impact of projected changes in climate on groundwater recharge in the near (2021-2050) and far (2071-2100) future across the western US. Recharge is expected to decrease slightly (highly certain) in the West (-1.6%) and Southwest (-2.9%) regions in the near future and decrease considerably (highly certain) in the South region (-10.6%) in the far future. The Northern Rockies region is expected to get more recharge (highly certain) in both the near (+5.0%) and far (+9.0%) future. In general, southern portions of the western US are expected to get less recharge in the future and northern portions will get more. This study also shows that climate change interacts with land surface properties to affect the amount of recharge that occurs in the future.

The Future of Earth Observation in Hydrology

In just the past five years, the field of Earth observation has progressed beyond the offerings of conventional space agency based platforms to include a plethora of sensing opportunities afforded by CubeSats, Unmanned Aerial Vehicles (UAVs), and smartphone technologies that are being embraced by both for-profit companies and individual researchers. Over the previous decades, space agency efforts have brought forth well-known and immensely useful satellites such as the Landsat series and the Gravity Research and Climate Experiment (GRACE) system, with costs typically on the order of one billion dollars per satellite and with concept-to-launch timelines on the order of two decades (for new missions). More recently, the proliferation of smartphones has helped to miniaturise sensors and energy requirements, facilitating advances in the use of CubeSats that can be launched by the dozens, while providing ultra-high (3-5 m) resolution sensing of the Earth on a daily basis. Start-up companies that did not exist five years ago now operate more satellites in orbit than any space agency, and at costs that are a mere fraction of the cost of traditional satellite missions. With these advances come new space-borne measurements, such as real-time high-definition video for tracking air pollution, storm-cell development, flood propagation, precipitation monitoring, or even for constructing digital surfaces using structure-from-motion techniques. Closer to the surface, measurements from small unmanned drones and tethered balloons have mapped snow depths, floods, and estimated evaporation at sub-meter resolutions, pushing back on spatio-temporal constraints and delivering new process insights. At ground level, precipitation has been measured using signal attenuation between antennae mounted on cell phone towers, while the proliferation of mobile devices has enabled citizen-scientists to catalogue photos of environmental conditions, estimate daily average temperatures from battery state, and sense other hydrologically important variables such as channel depths using commercially available wireless devices. Global internet access is being pursued via high altitude balloons, solar planes, and hundreds of planned satellite launches, providing a means to exploit the Internet of Things as an entirely new measurement domain. Such global access will enable real-time collection of data from billions of smartphones or from remote research platforms. This future will produce petabytes of data that can only be accessed via cloud storage and will require new analytical approaches to interpret. The extent to which today's hydrologic models can usefully ingest such massive data volumes is unclear. Nor is it clear whether this deluge of data will be usefully exploited, either because the measurements are superfluous, inconsistent, not accurate enough, or simply because we lack the capacity to process and analyse them. What is apparent is that the tools and techniques afforded by this array of novel and game-changing sensing platforms present our community with a unique opportunity to develop new insights that advance fundamental aspects of the hydrological sciences. To accomplish this will require more than just an application of the technology: in some cases, it will demand a radical rethink on how we utilise and exploit these new observing systems to enhance our understanding of the Earth and its linked processes.

Assessment of Irrigation Physics in a Land Surface Modeling Framework using Non-Traditional and Human-Practice Datasets

Irrigation increases soil moisture, which in turn controls water and energy fluxes from the land surface to the planetary boundary layer and determines plant stress and productivity. Therefore, developing a realistic representation of irrigation is critical to understanding land-atmosphere interactions in agricultural areas. Irrigation parameterizations are becoming more common in land surface models and are growing in sophistication, but there is difficulty in assessing the realism of these schemes, due to limited observations (e.g., soil moisture, evapotranspiration) and scant reporting of irrigation timing and quantity. This study uses the Noah land surface model run at high resolution within NASA's Land Information System to assess the physics of a sprinkler irrigation simulation scheme and model sensitivity to choice of irrigation intensity and greenness fraction datasets over a small, high resolution domain in Nebraska. Differences between experiments are small at the interannual scale but become more apparent at seasonal and daily time scales. In addition, this study uses point and gridded soil moisture observations from fixed and roving Cosmic Ray Neutron Probes and co-located human practice data to evaluate the realism of irrigation amounts and soil moisture impacts simulated by the model. Results show that field-scale heterogeneity resulting from the individual actions of farmers is not captured by the model and the amount of irrigation applied by the model exceeds that applied at the two irrigated fields. However, the seasonal timing of irrigation and soil moisture contrasts between irrigated and non-irrigated areas are simulated well by the model. Overall, the results underscore the necessity of both high-quality meteorological forcing data and proper representation of irrigation for accurate simulation of water and energy states and fluxes over cropland.

Tradeoff between cost and accuracy in large-scale surface water dynamic modeling

Recent efforts have led to the development of the local inertia formulation (INER) for an accurate but still cost-efficient representation of surface water dynamics, compared to the widely used kinematic wave equation (KINE). In this study, both formulations are evaluated over the Amazon basin in terms of computational costs and accuracy in simulating streamflows and water levels through synthetic experiments and comparisons against ground-based observations. Varying time steps are considered as part of the evaluation and INER at 60-second time step is adopted as the reference for synthetic experiments. Five hybrid (HYBR) realizations are performed based on maps representing the spatial distribution of the two formulations that physically represent river reach flow dynamics within the domain. Maps have fractions of KINE varying from 35.6% to 82.8%. KINE runs show clear deterioration along the Amazon river and main tributaries, with maximum RMSE values for streamflow and water level reaching 7827m3.s-1 and 1379cm near the basin's outlet. However, KINE is at least 25% more efficient than INER with low model sensitivity to longer time steps. A significant improvement is achieved with HYBR, resulting in maximum RMSE values of 3.9-292m3.s-1 for streamflows and 1.1-28.5cm for water levels, and cost reduction of 6-16%, depending on the map used. Optimal results using HYBR are obtained when the local inertia formulation is used in about one third of the Amazon basin, reducing computational costs in simulations while preserving accuracy. However, that threshold may vary when applied to different regions, according to their hydrodynamics and geomorphological characteristics.

Benefits and Pitfalls of GRACE Data Assimilation: a Case Study of Terrestrial Water Storage Depletion in India

This study investigates some of the benefits and drawbacks of assimilating Terrestrial Water Storage (TWS) observations from the Gravity Recovery and Climate Experiment (GRACE) into a land surface model over India. GRACE observes TWS depletion associated with anthropogenic groundwater extraction in northwest India. The model, however, does not represent anthropogenic groundwater withdrawals and is not skillful in reproducing the interannual variability of groundwater. Assimilation of GRACE TWS introduces long-term trends and improves the interannual variability in groundwater. But the assimilation also introduces a negative trend in simulated evapotranspiration whereas in reality evapotranspiration is likely enhanced by irrigation, which is also unmodeled. Moreover, in situ measurements of shallow groundwater show no trend, suggesting that the trends are erroneously introduced by the assimilation into the modeled shallow groundwater, when in reality the groundwater is depleted in deeper aquifers. The results emphasize the importance of representing anthropogenic processes in land surface modeling and data assimilation systems.

Comparing potential recharge estimates from three Land Surface Models across the Western US

Groundwater is a major source of water in the western US. However, there are limited recharge estimates available in this region due to the complexity of recharge processes and the challenge of direct observations. Land surface Models (LSMs) could be a valuable tool for estimating current recharge and projecting changes due to future climate change. In this study, simulations of three LSMs (Noah, Mosaic and VIC) obtained from the North American Land Data Assimilation System (NLDAS-2) are used to estimate potential recharge in the western US. Modeled recharge was compared with published recharge estimates for several aquifers in the region. Annual recharge to precipitation ratios across the study basins varied from 0.01-15% for Mosaic, 3.2-42% for Noah, and 6.7-31.8% for VIC simulations. Mosaic consistently underestimates recharge across all basins. Noah captures recharge reasonably well in wetter basins, but overestimates it in drier basins. VIC slightly overestimates recharge in drier basins and slightly underestimates it for wetter basins. While the average annual recharge values vary among the models, the models were consistent in identifying high and low recharge areas in the region. Models agree in seasonality of recharge occurring dominantly during the spring across the region. Overall, our results highlight that LSMs have the potential to capture the spatial and temporal patterns as well as seasonality of recharge at large scales. Therefore, LSMs (specifically VIC and Noah) can be used as a tool for estimating future recharge rates in data limited regions.

Spatio-temporal variability of groundwater storage in India

Groundwater level measurements from 3907 monitoring wells, distributed within 22 major river basins of India, are assessed to characterize their spatial and temporal variability. Groundwater storage (GWS) anomalies (relative to the long-term mean) exhibit strong seasonality, with annual maxima observed during the monsoon season and minima during pre-monsoon season. Spatial variability of GWS anomalies increases with the extent of measurements, following the power law relationship, i.e., log-(spatial variability) is linearly dependent on log-(spatial extent). In addition, the impact of well spacing on spatial variability and the power law relationship is investigated. We found that the mean GWS anomaly sampled at a 0.25 degree grid scale closes to unweighted average over all wells. The absolute error corresponding to each basin grows with increasing scale, i.e., from 0.25 degree to 1 degree. It was observed that small changes in extent could create very large changes in spatial variability at large grid scales. Spatial variability of GWS anomaly has been found to vary with climatic conditions. To our knowledge, this is the first study of the effects of well spacing on groundwater spatial variability. The results may be useful for interpreting large scale groundwater variations from unevenly spaced or sparse groundwater well observations or for siting and prioritizing wells in a network for groundwater management. The output of this study could be used to maintain a cost effective groundwater monitoring network in the study region and the approach can also be used in other parts of the globe.

Quantifying renewable groundwater stress with GRACE

Renewable groundwater stress is quantified in the world's largest aquifersCharacteristic stress regimes are defined to determine the severity of stressOverstressed aquifers are mainly in rangeland biomes with some croplands.

Uncertainty in global groundwater storage estimates in a Total Groundwater Stress framework

Groundwater resilience is defined and quantified with remote sensing from GRACETimescales of aquifer depletion are assessed as a Total Groundwater Stress ratioThe volume of usable global groundwater storage is found to be largely unknown.

Groundwater depletion during drought threatens future water security of the Colorado River Basin

Streamflow of the Colorado River Basin is the most overallocated in the world. Recent assessment indicates that demand for this renewable resource will soon outstrip supply, suggesting that limited groundwater reserves will play an increasingly important role in meeting future water needs. Here we analyze 9 years (December 2004 to November 2013) of observations from the NASA Gravity Recovery and Climate Experiment mission and find that during this period of sustained drought, groundwater accounted for 50.1 km³ of the total 64.8 km³ of freshwater loss. The rapid rate of depletion of groundwater storage (-5.6 ± 0.4 km³ yr^-1) far exceeded the rate of depletion of Lake Powell and Lake Mead. Results indicate that groundwater may comprise a far greater fraction of Basin water use than previously recognized, in particular during drought, and that its disappearance may threaten the long-term ability to meet future allocations to the seven Basin states.

Groundwater depletion in the Middle East from GRACE with implications for transboundary water management in the Tigris-Euphrates-Western Iran region

In this study, we use observations from the Gravity Recovery and Climate Experiment (GRACE) satellite mission to evaluate freshwater storage trends in the north-central Middle East, including portions of the Tigris and Euphrates River Basins and western Iran, from January 2003 to December 2009. GRACE data show an alarming rate of decrease in total water storage of approximately -27.2±0.6 mm yr^-1 equivalent water height, equal to a volume of 143.6 km³ during the course of the study period. Additional remote-sensing information and output from land surface models were used to identify that groundwater losses are the major source of this trend. The approach used in this study provides an example of "best current capabilities" in regions like the Middle East, where data access can be severely limited. Results indicate that the region lost 17.3±2.1 mm yr^-1 equivalent water height of groundwater during the study period, or 91.3±10.9 km³ in volume. Furthermore, results raise important issues regarding water use in transboundary river basins and aquifers, including the necessity of international water use treaties and resolving discrepancies in international water law, while amplifying the need for increased monitoring for core components of the water budget.

Satellite-based estimates of groundwater depletion in India

Groundwater is a primary source of fresh water in many parts of the world. Some regions are becoming overly dependent on it, consuming groundwater faster than it is naturally replenished and causing water tables to decline unremittingly. Indirect evidence suggests that this is the case in northwest India, but there has been no regional assessment of the rate of groundwater depletion. Here we use terrestrial water storage-change observations from the NASA Gravity Recovery and Climate Experiment satellites and simulated soil-water variations from a data-integrating hydrological modelling system to show that groundwater is being depleted at a mean rate of 4.0 +/- 1.0 cm yr(-1) equivalent height of water (17.7 +/- 4.5 km(3) yr(-1)) over the Indian states of Rajasthan, Punjab and Haryana (including Delhi). During our study period of August 2002 to October 2008, groundwater depletion was equivalent to a net loss of 109 km(3) of water, which is double the capacity of India's largest surface-water reservoir. Annual rainfall was close to normal throughout the period and we demonstrate that the other terrestrial water storage components (soil moisture, surface waters, snow, glaciers and biomass) did not contribute significantly to the observed decline in total water levels. Although our observational record is brief, the available evidence suggests that unsustainable consumption of groundwater for irrigation and other anthropogenic uses is likely to be the cause. If measures are not taken soon to ensure sustainable groundwater usage, the consequences for the 114,000,000 residents of the region may include a reduction of agricultural output and shortages of potable water, leading to extensive socioeconomic stresses.