Satellite Gravimetry for Detecting Water Scarcity
Explore the source record for details and available documents.
Bailing Li↗
Engineering topics
Bailing Li
Publications and source records attributed to Bailing Li.
Terrestrial Water Storage in 2023
Global terrestrial water storage (TWS) anomalies reached a record low of –9.94 cm in 2023, decreasing 0.80 cm from 2022. These reductions largely reflect ongoing TWS losses from glacial melt and groundwater use for irrigation, offset by gains in central and eastern Antarctica and La Niña-related tropical wetting.
Bailing Li↗
How Have Hydrological Extremes Changed Over the Past 20 Years?
Severe floods and droughts, including their back-to-back occurrences (weather whiplash), have been increasing in frequency and severity around the world. Improved understanding of systematic changes in hydrological extremes is essential for preparation and adaptation. In this study, we identified and quantified extreme wet and dry events globally by applying a clustering algorithm to terrestrial water storage (TWS) data from the Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow-On (FO). The most intense events, ranked using an intensity metric, often reflect impacts of large-scale oceanic oscillations such as El Niño–Southern Oscillation and consequences of climate change. The severity of both wet and dry events, represented by standardized TWS anomalies, increased significantly in most cases, likely associated with intensification of wet and dry weather regimes in a warmer world, and consequently, exhibited strongest correlation with global temperature. In the Dry climate, the number of wet events decreased while the number of dry events increased significantly, suggesting a drying trend that may be attributed to climate variability and possible increases in irrigation and reliance on groundwater. In the Continental climate where temperature has risen faster than global average, dry events increased significantly. Characteristics of extreme events often showed strong correlations with global temperature, especially when averaged over all climates. These results suggest changes in hydrological extremes and underscore the importance of quantifying total water storage changes when studying hydrological extremes. Extending the GRACE/FO record, which spans 2002 to the present, is essential to continuously tracking changes in TWS and hydrological extremes.
Bailing Li↗
Changing Intensity of Hydroclimatic Extreme Events Revealed by GRACE and GRACE-FO
Distortion of the water cycle, particularly of its extremes (droughts and pluvials), will be among the most conspicuous consequences of climate change. We applied a novel approach with terrestrial water storage observations from the GRACE and GRACE-FO satellites to delineate and characterize 1,056 extreme events during 2002-2021. Dwarfing all other events was an ongoing pluvial that began in 2019 and engulfed central Africa. Total intensity of extreme events was strongly correlated with global mean temperature, more so than with the El Nino Southern Oscillation or other climate indicators, suggesting that continued warming of the planet will cause more frequent, more severe, longer, and/or larger droughts and pluvials. In three regions, including a vast swath extending from southern Europe to southwestern China, the ratio of wet to dry extreme events decreased substantially over the study period, while the opposite was true in two regions, including sub-Saharan Africa from 5°N to 20°N.
Climate change↗
A High-Resolution Land Data Assimilation System Optimized for the Western United States
The Western Land Data Assimilation System (WLDAS) is a custom instance of the NASA Land Information System that combines land surface parameters, meteorological forcing data, and satellite products within a land surface model to produce daily estimates of the water and energy budget variables for the western United States. WLDAS was configured through discussions with partners, with the goal of groundwater sustainability planning for the state of California in mind. The publicly available output dataset has a 1-km grid resolution and spans 1979–present, which makes it suitable for water resources assessments. The data are also able to contextualize the recent drought events in California. Assimilation of Leaf Area Index, which is demonstrated herein to improve simulation over agricultural areas in California, specifically in terms of evapotranspiration in irrigation regions, will be included along with other data assimilation in subsequent releases of WLDAS.
land surface modeling↗
Groundwater Recharge Estimated by Land Surface Models: An Evaluation in the Conterminous U.S.
Estimating diffuse recharge of precipitation is fundamental to assessing groundwater sustainability. Diffuse recharge is also the process through which climate and climate change directly affect groundwater. In this study, we evaluated diffuse recharge over the conterminous U.S. simulated by a suite of land surface models (LSMs) that were forced using a common set of meteorological input data. Simulated annual recharge exhibited spatial patterns that were similar among the LSMs, with the highest values in the eastern U.S. and Pacific Northwest. However, the magnitudes of annual recharge varied significantly among the models and were associated with differences in simulated ET, runoff and snow. Evaluation against two independent datasets did not answer the question of whether the ensemble mean performs the best, due to inconsistency between those datasets. The amplitude and timing of seasonal maximum recharge differed among the models, influenced strongly by model physics governing deep soil moisture drainage rates and, in cold regions, snowmelt. Evaluation using in situ soil moisture observations suggested that true recharge peaks 1-3 months later than simulated recharge, indicating systematic biases in simulating deep soil moisture. However, recharge from lateral flows and through preferential flows cannot be inferred from soil moisture data, and the seasonal cycle of simulated groundwater storage actually compared well with in situ groundwater observations. Long-term trends in recharge were not consistently correlated with either precipitation trends or temperature trends. This study highlights the need to employ dynamic flow models in LSMs, among other improvements, to enable more accurate simulation of recharge.
snowmelt↗
Groundwater Drought: Environmental Controls and Monitoring
Arising from a long-term deficit of precipitation, groundwater drought often lags meteorological drought by months to years, depending on climate and subsurface hydrogeological conditions. Due to the paucity of in situ groundwater data in many parts of the world, global monitoring of groundwater storage variations and drought may be best achieved using satellite observations and/or groundwater time series simulated by hydrological models. The Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow On (GRACE-FO) missions have greatly benefitted the modeling and monitoring of groundwater storage changes and drought at the global scale. In this chapter, we first review environmental controls on the temporal variability of groundwater using in situ data. We then describe an approach that infuses GRACE and GRACE-FO observations into a land surface model for assessing groundwater storage changes and drought globally. We also discuss characteristics of simulated groundwater drought and the limitations of current groundwater drought monitoring approaches.
Bailing Li↗
Groundwater and Terrestrial Water Storage
Groundwater, soil moisture, surface water, snow, and ice are the components of terrestrial water storage (TWS). On multi-annual timescales groundwater typically controls TWS variability, except in permanently frozen regions (Li et al. 2015). Even on an annual basis TWS changes are a reasonable proxy for groundwater storage changes, the latter being insufficiently observed in most of the world. From 2002 to 2017 the Gravity Recovery and Climate Experiment (GRACE; Tapley et al. 2004) and since 2018 the GRACE Follow On (GRACE-FO) satellite missions have enabled estimation of TWS anomalies (departures from the long-term mean) based on precise observations of variations in Earth’s gravity field. To bridge the 11-month gap between GRACE and GRACE-FO we make use of output from a land surface model that assimilates data from both missions (Li et al. 2019)
Matthew Rodell↗
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-47 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 degreesx0.125 degrees 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 in the Upper Great Plains and Northeast to -1 to -9 in the West and South; net radiation flux trends range from +0.05 to +0.20 W/m(exp2)/yr in the East to -0.05 to -0.20 in the West; U.S.-wide temperature trends average about +0.03K/yr. 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 while trend detection increased 11%; Trends in latent heat flux were hardly affected, RMSE decreasing only 0.17 to 0.16 W/m(exp 2)/yr, while trend detection increased 2%; NCA-LDAS runoff trends degraded significantly from 2.6 to 16.1 mm/yr 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 gaged regions. Overall, NCA-LDAS demonstrates capability for quantifying physically consistent, U.S. hydrologic climate trends over the satellite era.
Snow depth↗
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° 3 0.125° 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 in the upper Great Plains and Northeast to -1 to -9 mm/yr in the West and South, net radiation flux trends range from +0.05 to +0.20 W/sq. m yr in the East to -0.05 to -0.20 W/sq. m yr in the West, and U.S.-wide temperature trends average about +0.03K/yr. 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 while trend detection increased 11%. Trends in latent heat flux were hardly affected, with RMSE decreasing only from 0.17 to 0.16 W/sq. m yr, while trend detection increased 2%. NCA-LDAS runoff trends degraded significantly from 2.6 to 16.1 mm/yr 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.
Michael F. Jasinski↗
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.
Groundwater storage↗
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.
Bailing Li↗
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. Ground water 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.
Groundwater monitoring↗
Groundwater Variability Across Temporal and Spatial Scales in the Central and Northeastern U.S.
Depth-to-water measurements from 181 monitoring wells in unconfined or semi-confined aquifers in nine regions of the central and northeastern U.S. were analyzed. Groundwater storage exhibited strong seasonal variations in all regions, with peaks in spring and lows in autumn, and its interannual variability was nearly unbounded, such that the impacts of droughts, floods, and excessive pumping could persist for many years. We found that the spatial variability of groundwater storage anomalies (deviations from the long term mean) increases as a power function of extent scale (square root of area). That relationship, which is linear on a log-log graph, is common to other hydrological variables but had never before been shown with groundwater data. We describe how the derived power function can be used to determine the number of wells needed to estimate regional mean groundwater storage anomalies with a desired level of accuracy, or to assess uncertainty in regional mean estimates from a set number of observations. We found that the spatial variability of groundwater storage anomalies within a region often increases with the absolute value of the regional mean anomaly, the opposite of the relationship between soil moisture spatial variability and mean. Recharge (drainage from the lowest model soil layer) simulated by the Variable Infiltration Capacity (VIC) model was compatible with observed monthly groundwater storage anomalies and month-to-month changes in groundwater storage.
Groundwater↗
Groundwater Variability Across Temporal and Spatial Scales in the Central and Northeastern U.S.
Depth-to-water measurements from 181 monitoring wells in unconfined or semi-confined aquifers in nine regions of the central and northeastern U.S. were analyzed. Groundwater storage exhibited strong seasonal variations in all regions, with peaks in spring and lows in autumn, and its interannual variability was nearly unbounded, such that the impacts of droughts, floods, and excessive pumping could persist for many years. We found that the spatial variability of groundwater storage anomalies (deviations from the long term mean) increases as a power function of extent scale (square root of area). That relationship, which is linear on a log-log graph, is common to other hydrological variables but had never before been shown with groundwater data. We describe how the derived power function can be used to determine the number of wells needed to estimate regional mean groundwater storage anomalies with a desired level of accuracy, or to assess uncertainty in regional mean estimates from a set number of observations. We found that the spatial variability of groundwater storage anomalies within a region often increases with the absolute value of the regional mean anomaly, the opposite of the relationship between soil moisture spatial variability and mean. Recharge (drainage from the lowest model soil layer) simulated by the Variable Infiltration Capacity (VIC) model was compatible with observed monthly groundwater storage anomalies and month-to-month changes in groundwater storage.
water cycle↗
Evaluation of a Model-Based Groundwater Drought Indicator in the Conterminous U.S.
Monitoring groundwater drought using land surface models is a valuable alternative given the current lack of systematic in situ measurements at continental and global scales and the low resolution of current remote sensing based groundwater data. However, uncertainties inherent to land surface models may impede drought detection, and thus should be assessed using independent data sources. In this study, we evaluated a groundwater drought index (GWI) derived from monthly groundwater storage output from the Catchment Land Surface Model (CLSM) using a GWI similarly derived from in situ groundwater observations. Groundwater observations were obtained from unconfined or semi-confined aquifers in eight regions of the central and northeastern U.S. Regional average GWI derived from CLSM exhibited strong correlation with that from observation wells, with correlation coefficients between 0.43 and 0.92. GWI from both in situ data and CLSM was generally better correlated with the Standard Precipitation Index (SPI) at 12 and 24 month timescales than at shorter timescales, but it varied depending on climate conditions. The correlation between CLSM derived GWI and SPI generally decreases with increasing depth to the water table, which in turn depends on both bedrock depth (a CLSM parameter) and mean annual precipitation. The persistence of CLSM derived GWI is spatially varied and again shows a strong influence of depth to groundwater. CLSM derived GWI generally persists longer than GWI derived from in situ data, due at least in part to the inability of coarse model inputs to capture high frequency meteorological variability at local scales. The study also showed that groundwater can have a significant impact on soil moisture persistence where the water table is shallow. Soil moisture persistence was estimated to be longer in the eastern U.S. than in the west, in contrast to previous findings that were based on models that did not represent groundwater. Assimilation of terrestrial water storage data from the Gravity Recovery and Climate Experiment (GRACE) satellite mission improved the correlation between CLSM based regional average GWI and that based on in situ data in six of the eight regions. Practical issues regarding the application of GRACE assimilated groundwater storage for drought detection are discussed. An important conclusion of this study is that model parameters that control the depth to the water table, including bedrock depth, strongly influence the evolution and persistence of simulated groundwater and require careful configuration for drought monitoring.
hydrology↗
Spatial and Temporal Variability of Recharge in the U.S.
Explore the source record for details and available documents.
subsurface runoff↗
A Visual Guide to the GRACE Data Assimilation for Drought Monitoring (GRACE-DA-DM) Data Products
Explore the source record for details and available documents.
Stephanie Stettz↗