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 depletion in California's Central Valley accelerates during megadrought.

Groundwater provides nearly half of irrigation water supply, and it enables resilience during drought, but in many regions of the world, it remains poorly, if at all managed. In heavily agricultural regions like California's Central Valley, where groundwater management is being slowly implemented over a 27-year period that began in 2015, groundwater provides two-thirds or more of irrigation water during drought, which has led to falling water tables, drying wells, subsiding land, and its long-term disappearance. Here we use nearly two decades of observations from NASA's GRACE satellite missions and show that the rate of groundwater depletion in the Central Valley has been accelerating since 2003 (1.86 km3/yr, 1961-2021; 2.41 km3/yr, 2003-2021; 8.58 km3/yr, 2019-2021), a period of megadrought in southwestern North America. Results suggest the need for expedited implementation of groundwater management in the Central Valley to ensure its availability during the increasingly intense droughts of the future.

A 10 per cent increase in global land evapotranspiration from 2003 to 2019.

Accurate quantification of global land evapotranspiration is necessary for understanding variability in the global water cycle, which is expected to intensify under climate change1-3. Current global evapotranspiration products are derived from a variety of sources, including models4,5, remote sensing6,7 and in situ observations8-10. However, existing approaches contain extensive uncertainties; for example, relating to model structure or the upscaling of observations to a global level11. As a result, variability and trends in global evapotranspiration remain unclear12. Here we show that global land evapotranspiration increased by 10 ± 2 per cent between 2003 and 2019, and that land precipitation is increasingly partitioned into evapotranspiration rather than runoff. Our results are based on an independent water-balance ensemble time series of global land evapotranspiration and the corresponding uncertainty distribution, using data from the Gravity Recovery and Climate Experiment (GRACE) and GRACE-Follow On (GRACE-FO) satellites13. Variability in global land evapotranspiration is positively correlated with El Niño-Southern Oscillation. The main driver of the trend, however, is increasing land temperature. Our findings provide an observational constraint on global land evapotranspiration, and are consistent with the hypothesis that global evapotranspiration should increase in a warming climate.

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 km3/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.

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.

Long-term, non-anthropogenic groundwater storage changes simulated by three global-scale hydrological models.

This study examined long-term, natural (i.e., excluding anthropogenic impacts) variability of groundwater storage worldwide. Groundwater storage changes were estimated by forcing three global-scale hydrological models with three 50+ year meteorological datasets. Evaluation using in situ groundwater observations from the U.S. and terrestrial water storage derived from the Gravity Recovery and Climate Experiment (GRACE) satellites showed that these models reasonably represented inter-annual variability of water storage, as indicated by correlations greater than 0.5 in most regions. Empirical orthogonal function analysis revealed influences of the El Niño Southern Oscillation (ENSO) on global groundwater storage. Simulated groundwater storage, including its global average, exhibited trends generally consistent with that of precipitation. Global total (natural) groundwater storage decreased over the past 5-7 decades with modeled rates ranging from 0.01 to 2.18 mm year-1. This large range can be attributed in part to groundwater's low frequency (inter-decadal) variability, which complicates identification of real long-term trends even within a 50+ year time series. Results indicate that non-anthropogenic variability in groundwater storage is substantial, making knowledge of it fundamental to quantifying direct human impacts on groundwater storage.

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.

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.

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.

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).

Groundwater rejuvenation in parts of India influenced by water-policy change implementation.

The dwindling groundwater resource of India, supporting almost one fifth of the global population and also the largest groundwater user, has been of great concern in recent years. However, in contrary to the well documented Indian groundwater depletion due to rapid and unmanaged groundwater withdrawal, here for the first time, we report regional-scale groundwater storage (GWS) replenishment through long-term (1996-2014, using more than 19000 observation locations) in situ and decadal (2003-2014) satellite-based groundwater storage measurements in western and southern parts of India. In parts of western and southern India, in situ GWS (GWSobs) has been decreasing at the rate of -5.81 ± 0.38 km3/year (in 1996-2001) and -0.92 ± 0.12 km3/year (in 1996-2002), and reversed to replenish at the rate of 2.04 ± 0.20 km3/year (in 2002-2014) and 0.76 ± 0.08 km3/year (in 2003-2014), respectively. Here, using statistical analyses and simulation results of groundwater management policy change effect on groundwater storage in western and southern India, we show that paradigm shift in Indian groundwater withdrawal and management policies for sustainable water utilization appear to have started replenishing the aquifers in western and southern parts of India.

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.

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.

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.

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.

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.