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.

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.

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.

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.

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.

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.

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.

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.

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 km3 of the total 64.8 km3 of freshwater loss. The rapid rate of depletion of groundwater storage (-5.6 ± 0.4 km3 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 km3 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 km3 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.

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.

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.

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.

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.

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.