Seasonal Water Storage Variations as Impacted by Water Abstractions: Comparing the Output of a Global Hydrological Model with GRACE and GPS Observations

Better quantification of continental water storage variations is expected to improve our understanding of water flows, including evapotranspiration, runoff and river discharge as well as human water abstractions. For the first time, total water storage (TWS) on the land area of the globe as computed by the global water model WaterGAP (Water Global Assessment and Prognosis) was compared to both gravity recovery and climate experiment (GRACE) and global positioning system (GPS) observations. The GRACE satellites sense the effect of TWS on the dynamic gravity field of the Earth. GPS reference points are displaced due to crustal deformation caused by time-varying TWS. Unfortunately, the worldwide coverage of the GPS tracking network is irregular, while GRACE provides global coverage albeit with low spatial resolution. Detrended TWS time series were analyzed by determining scaling factors for mean annual amplitude (f _GRACE) and time series of monthly TWS (f _GPS). Both GRACE and GPS indicate that WaterGAP underestimates seasonal variations of TWS on most of the land area of the globe. In addition, seasonal maximum TWS occurs 1 month earlier according to WaterGAP than according to GRACE on most land areas. While WaterGAP TWS is sensitive to the applied climate input data, none of the two data sets result in a clearly better fit to the observations. Due to the low number of GPS sites, GPS observations are less useful for validating global hydrological models than GRACE observations, but they serve to support the validity of GRACE TWS as observational target for hydrological modeling. For unknown reasons, WaterGAP appears to fit better to GPS than to GRACE. Both GPS and GRACE data, however, are rather uncertain due to a number of reasons, in particular in dry regions. It is not possible to benefit from either GPS or GRACE observations to monitor and quantify human water abstractions if only detrended (seasonal) TWS variations are considered. Regarding GRACE, this is mainly caused by the attenuation of the TWS differences between water abstraction variants due to the filtering required for GRACE TWS. Regarding GPS, station density is too low. Only if water abstractions lead to long-term changes in TWS by depletion or restoration of water storage in groundwater or large surface water bodies, GRACE may be used to support the quantification of human water abstractions.     

Groundwater Storage Changes: Present Status from GRACE Observations

Satellite gravity measurements from the Gravity Recovery and Climate Experiment (GRACE) provide quantitative measurement of terrestrial water storage (TWS) changes with unprecedented accuracy. Combining GRACE-observed TWS changes and independent estimates of water change in soil and snow and surface reservoirs offers a means for estimating groundwater storage change. Since its launch in March 2002, GRACE time-variable gravity data have been successfully used to quantify long-term groundwater storage changes in different regions over the world, including northwest India, the High Plains Aquifer and the Central Valley in the USA, the North China Plain, Middle East, and southern Murray–Darling Basin in Australia, where groundwater storage has been significantly depleted in recent years (or decades). It is difficult to rely on in situ groundwater measurements for accurate quantification of large, regional-scale groundwater storage changes, especially at long timescales due to inadequate spatial and temporal coverage of in situ data and uncertainties in storage coefficients. The now nearly 13 years of GRACE gravity data provide a successful and unique complementary tool for monitoring and measuring groundwater changes on a global and regional basis. Despite the successful applications of GRACE in studying global groundwater storage change, there are still some major challenges limiting the application and interpretation of GRACE data. In this paper, we present an overview of GRACE applications in groundwater studies and discuss if and how the main challenges to using GRACE data can be addressed.     

GRACE估算陆地水储量季节和年际变化

利用最新公布的GRACE GFZ RL04数据,分析了2003年1月～2007年12月全球27条流域和陆地水储量的季节性和年际变化.结果表明,相近流域季节性变化相位接近.2003年1月～2007年12月陆地水储量季节性变化为1572.4 km3,其中变化最大流域为亚马逊河,其次分别为鄂毕河、尼罗河和尼日尔河等流域.5年来 GRACE陆地水储量的年际变化为-75.4±40.3 km3/a,其中亚马逊河、勒拿河和马更些河等流域的年际变化呈现正增长,而刚果河、密西西比河、恒河、育空河和雅鲁藏布江等流域则相反.GRACE与GLDAS数据均表明2006年后陆地水储量年际变化存在明显增加.     

Terrestrial Water Storage Anomalies Associated with Drought in Southwestern USA from GPS Observations

The Earth’s surface fluid mass redistribution, e.g., groundwater depletion and severe drought, causes the elastic surface deformation, which can be measured by global positioning system (GPS). In this paper, the continuous GPS observations are used to estimate the terrestrial water storage (TWS) changes in southwestern USA, which have a good agreement with TWS changes derived from Gravity Recovery And Climate Experiment (GRACE) and hydrological models. The seasonal variation is mostly located in the Rocky mountain range and Mississippi river watershed. The largest amplitude of the seasonal variation is between 12 and 15 cm in equivalent water thickness. The timing and duration of TWS anomalies caused by the severe drought in 2012 are observed by the GPS-derived TWS, which are confirmed by the GRACE results. Different hydrological models are further used for comparison with GPS and GRACE results. The magnitude of TWS depletion from GRACE and GPS observations during the drought is larger than that from hydrological models, which indicates that the drought was caused by comparable groundwater and surface water depletion. The interannual TWS changes from GPS are also consistent with the precipitation pattern over the past 6 years, which further confirms the severe drought in 2012. This study demonstrates that continuous GPS observations have the potential as real-time drought indicator.     

Large-Scale Total Water Storage and Water Flux Changes over the Arid and Semiarid Parts of the Middle East from GRACE and Reanalysis Products

Previous studies indicate that water storage over a large part of the Middle East has been decreased over the last decade. Variability in the total (hydrological) water flux (TWF, i.e., precipitation minus evapotranspiration minus runoff) and water storage changes of the Tigris–Euphrates river basin and Iran’s six major basins (Khazar, Persian, Urmia, Markazi, Hamun, and Sarakhs) over 2003–2013 is assessed in this study. Our investigation is performed based on the TWF that are estimated as temporal derivatives of terrestrial water storage (TWS) changes from the Gravity Recovery and Climate Experiment (GRACE) products and those from the reanalysis products of ERA-Interim and MERRA-Land. An inversion approach is applied to consistently estimate the spatio-temporal changes of soil moisture and groundwater storage compartments of the seven basins during the study period from GRACE TWS, altimetry, and land surface model products. The influence of TWF trends on separated water storage compartments is then explored. Our results, estimated as basin averages, indicate negative trends in the maximums of TWF peaks that reach up to ?5.2 and ?2.6 (mm/month/year) over 2003–2013, respectively, for the Urmia and Tigris–Euphrates basins, which are most likely due to the reported meteorological drought. Maximum amplitudes of the soil moisture compartment exhibit negative trends of ?11.1, ?6.6, ?6.1, ?4.8, ?4.7, ?3.8, and ?1.2 (mm/year) for Urmia, Tigris–Euphrates, Khazar, Persian, Markazi, Sarakhs, and Hamun basins, respectively. Strong groundwater storage decrease is found, respectively, within the Khazar ?8.6 (mm/year) and Sarakhs ?7.0 (mm/year) basins. The magnitude of water storage decline in the Urmia and Tigris–Euphrates basins is found to be bigger than the decrease in the monthly accumulated TWF indicating a contribution of human water use, as well as surface and groundwater flow to the storage decline over the study area.     

Global surface mass from a new combination of GRACE,modelled OBP and reprocessed GPS data

Weekly surface loading variations are estimated from a joint least squares inversion of load-induced GPS site displacements, GRACE gravimetry and simulated ocean bottom pressure (OBP) from the finite element sea-ice ocean model (FESOM).In this study, we directly use normal equations derived from reprocessed GPS observations, where station and satellite positions are estimated simultaneously. The OBP weight of the model in the inversion is based on a new error model, obtained from 2 FESOM runs forced with different atmospheric data sets.Our findings indicate that the geocenter motion derived from the inversion is smooth, with non-seasonal RMS values of 1.4, 0.9 and 1.9 mm for the X, Y and Z directions, respectively. The absolute magnitude of the seasonal geocenter motion varies annually between 2 and 4.5 mm. Important hydrological regions such as the Amazon, Australia, South-East Asia and Europe are mostly affected by the geocenter motion, with magnitudes of up to 2 cm, when expressed in equivalent water height.The chosen solar radiation pressure model, used in the GPS processing, has only a marginal effect on the joint inversion results. Using the empirical CODE model slightly increases the annual amplitude of the Z component of the geocenter by 0.8 mm. However, in case of a GPS-only inversion, notable larger differences are found for the annual amplitude and phase estimates when applying the older physical ROCK models. Regardless of the used radiation pressure model the GPS network still exhibits maximum radial expansions in the order of 3 mm (0.45 ppb in terms of scale), which are most likely caused by remaining GPS technique errors.In an additional experiment, we have used the joint inversion solution as a background loading model in the GPS normal equations. The reduced time series, compared to those without a priori loading model, show a consistent decrease in RMS. In terms of the annual height component, 151 of the 189 stations show a reduction of at least 10% in seasonal amplitude.On the ocean floor, we find a positive overall correlation (0.51) of the inversion solution with time series from globally distributed independent bottom pressure recorders.Even after removing a seasonal fit we still find a correlation of 0.45. Furthermore, the geocenter motion has a significant effect on ocean bottom pressure as neglecting it causes the correlation to drop to 0.42.     

Groundwater Storage Change in the Jinsha River Basin from GRACE,Hydrologic Models,and In Situ Data

Groundwater plays a major role in the hydrological processes driven by climate change and human activities, particularly in upper mountainous basins. The Jinsha River Basin (JRB) is the uppermost region of the Yangtze River and the largest hydropower production region in China. With the construction of artificial cascade reservoirs increasing in this region, the annual and seasonal flows are changing and affecting the water cycles. Here, we first infer the groundwater storage changes (GWSC), accounting for sediment transport in JRB, by combining the Gravity Recovery and Climate Experiment mission, hydrologic models and in situ data. The results indicate: (1) the average estimation of the GWSC trend, accounting for sediment transport in JRB, is 0.76 ± 0.10 cm/year during the period 2003 to 2015, and the contribution of sediment transport accounts for 15%; (2) precipitation (P), evapotranspiration (ET), soil moisture change, GWSC, and land water storage changes (LWSC) show clear seasonal cycles; the interannual trends of LWSC and GWSC increase, but P, runoff (R), surface water storage change and SMC decrease, and ET remains basically unchanged; (3) the main contributor to the increase in LWSC in JRB is GWSC, and the increased GWSC may be dominated by human activities, such as cascade damming and climate variations (such as snow and glacier melt due to increased temperatures). This study can provide valuable information regarding JRB in China for understanding GWSC patterns and exploring their implications for regional water management.     

Global Terrestrial Water Storage Changes and Connections to ENSO Events

Improved data quality of extended record of the Gravity Recovery and Climate Experiment (GRACE) satellite gravity solutions enables better understanding of terrestrial water storage (TWS) variations. Connections of TWS and climate change are critical to investigate regional and global water cycles. In this study, we provide a comprehensive analysis of global connections between interannual TWS changes and El Niño Southern Oscillation (ENSO) events, using multiple sources of data, including GRACE measurements, land surface model (LSM) predictions and precipitation observations. We use cross-correlation and coherence spectrum analysis to examine global connections between interannual TWS changes and the Niño 3.4 index, and select four river basins (Amazon, Orinoco, Colorado, and Lena) for more detailed analysis. The results indicate that interannual TWS changes are strongly correlated with ENSO over much of the globe, with maximum cross-correlation coefficients up to ~0.70, well above the 95% significance level (~0.29) derived by the Monte Carlo experiments. The strongest correlations are found in tropical and subtropical regions, especially in the Amazon, Orinoco, and La Plata basins. While both GRACE and LSM TWS estimates show reasonably good correlations with ENSO and generally consistent spatial correlation patterns, notably higher correlations are found between GRACE TWS and ENSO. The existence of significant correlations in middle–high latitudes shows the large-scale impact of ENSO on the global water cycle.     

Assessing Groundwater Depletion and Dynamics Using GRACE and InSAR: Potential and Limitations

In the last decade, remote sensing of the temporal variation of ground level and gravity has improved our understanding of groundwater dynamics and storage. Mass changes are measured by GRACE (Gravity Recovery and Climate Experiment) satellites, whereas ground deformation is measured by processing synthetic aperture radar satellites data using the InSAR (Interferometry of Synthetic Aperture Radar) techniques. Both methods are complementary and offer different sensitivities to aquifer system processes. GRACE is sensitive to mass changes over large spatial scales (more than 100,000 km²). As such, it fails in providing groundwater storage change estimates at local or regional scales relevant to most aquifer systems, and at which most groundwater management schemes are applied. However, InSAR measures ground displacement due to aquifer response to fluid‐pressure changes. InSAR applications to groundwater depletion assessments are limited to aquifer systems susceptible to measurable deformation. Furthermore, the inversion of InSAR‐derived displacement maps into volume of depleted groundwater storage (both reversible and largely irreversible) is confounded by vertical and horizontal variability of sediment compressibility. During the last decade, both techniques have shown increasing interest in the scientific community to complement available in situ observations where they are insufficient. In this review, we present the theoretical and conceptual bases of each method, and present idealized scenarios to highlight the potential benefits and challenges of combining these techniques to remotely assess groundwater storage changes and other aspects of the dynamics of aquifer systems.     

Local hydrology, the Global Geodynamics Project and CHAMP/GRACE perspective: some case studies

We first quantify the influence of aquifers on gravity variations by considering local, regional and continental scales. We show that locally only the direct attraction of the underlying aquifer has to be taken into account. At continental (or global scales), the underground water masses act by direct attraction (due to the earth curvature), loading flexure and potential redistribution. We show that at the intermediate regional scale (saying a few kilometres to a few hundreds of kilometres), groundwater contributions can be neglected in practice. Afterwards, we illustrate the difficulties in tackling the local hydrological context by studying comparatively the geological and hydrogeological surroundings of three European Global Geodynamics Project (GGP) superconducting gravimeter stations (Strasbourg, Moxa, and Vienna). Finally, it appears clearly that hydrological variability and cycle characterisations constitute the up-to-date challenge while studying gravity variations in a large spectral range. That is why, gravity is used to quantify hydrological transfers, and overall when seeking for small signals from the Earth's deep interior or other environmental signals (atmosphere, oceans) where groundwater influence can be seen as a disturbance.     

Comparison of Daily GRACE Gravity Field and Numerical Water Storage Models for De-aliasing of Satellite Gravimetry Observations

Reducing aliasing effects of insufficiently modelled high-frequent, non-tidal mass variations of the atmosphere, the oceans and the hydrosphere in gravity field models derived from the Gravity Recovery and Climate Experiment (GRACE) satellite mission is the topic of this study. The signal content of the daily GRACE gravity field model series (ITG-Kalman) is compared to high-frequency bottom pressure variability and terrestrially stored water variations obtained from recent numerical simulations from an ocean circulation model (OMCT) and two hydrological models (WaterGAP Global Hydrology Model, Land Surface Discharge Model). Our results show that daily estimates of ocean bottom pressure from the most recent OMCT simulations and the daily ITG-Kalman solutions are able to explain up to 40 % of extra-tropical sea-level variability in the Southern Ocean. In contrast to this, the daily ITG-Kalman series and simulated continental total water storage variability largely disagree at periods below 30 days. Therefore, as long as no adequate hydrological model will become available, the daily ITG-Kalman series can be regarded as a good initial proxy for high-frequency mass variations at a global scale. As a second result of this study, based on monthly solutions as well as daily observation residuals, it is shown that applying this GRACE-derived de-aliasing model supports the determination of the time-variable gravity field from GRACE data and the subsequent geophysical interpretation. This leads us to the recommendation that future satellite concepts for determining mass variations in the Earth system should be capable of observing higher frequeny signals with sufficient spatial resolution.     

Thermal Anomalies and Earthquakes： Evidence from Wenchuan, China

Earthquake prediction is a difficult problem in Earth sciences.Unsuccessful predictions one after another urged people to explore more synthetic and comprehensive methods for earthquake prediction.The Lithosphere-Atmosphere-Ionosphere(LAI)coupling theory pays great attention to the processes taking place within the near ground layer of atmosphere.It has achieved great results recently,and can enlighten us about the nature of an earthquake's precursor.Based on the NCEP reanalysis dataset,this paper attempts to track the anomalies of the surface's upward long wave radiation flux(ULWRF),the temperature at the depth of 10cm～20cm below ground surface layer(BGL)and the air temperature at 2 meters above ground surface(AIR)around the time of the strong Wenchuan earthquake.Thermal anomalies were observed before and after May 12,2008,the time of the Wenchuan earthquake.Perhaps the thermal anomaly that occurred prior to the earthquake can be taken as indicators of the earthquake,but in view of the complexity of the earthquake phenomena,using thermal anomaly as a precursor should be done with caution.     

Assessing and Improving Land Surface Model Outputs Over Africa Using GRACE,Field, and Remote Sensing Data

The Gravity Recovery and Climate Experiment (GRACE), along with other relevant field and remote sensing datasets, was used to assess the performance of two land surface models (LSMs: CLM4.5-SP and GLDAS-Noah) over the African continent and improve the outputs of the CLM4.5-SP model. Spatial and temporal analysis of monthly (January 2003–December 2010) Terrestrial Water Storage (TWS) estimates extracted from GRACE (TWS_GRACE), CLM4.5-SP (TWS_CLM4.5), and GLDAS-Noah (TWS_GLDAS) indicates the following: (1) compared to GRACE, LSMs overestimate TWS in winter months and underestimate them in summer months; (2) the amplitude of annual cycle (AAC) of TWS_GRACE is higher than that of TWS_LSM (AAC: TWS_GRACE > TWS_GLDAS > TWS_CLM4.5); (3) higher, and statistically significant correlations were observed between TWS_GRACE and TWS_GLDAS compared to those between TWS_GRACE and TWS_CLM4.5; (4) differences in forcing precipitation and temperature datasets for GLDAS-Noah and CLM4.5-SP models are unlikely to be the main cause for the observed discrepancies between TWS_GRACE and TWS_LSM; and (5) the CLM4.5-SP model overestimates evapotranspiration (ET) values in summer months and underestimates them in winter months compared to ET estimates extracted from field-based (FLUXNET-MTE) and satellite-based (MOD16 and GLEAM) ET measurements. A first-order correction was developed and applied to correct the CLM4.5-derived ET, soil moisture, groundwater, and TWS. The corrections improved the correspondence (i.e., higher correlation and comparable AAC) between TWS_CLM4.5 and TWS_GRACE over various climatic settings. Our findings suggest that similar straightforward correction approaches could potentially be developed and used to assess and improve the performance of a wide range of LSMs.     

Toward a High-Resolution Monitoring of Continental Surface Water Extent and Dynamics,at Global Scale: from GIEMS (Global Inundation Extent from Multi-Satellites) to SWOT (Surface Water Ocean Topography)

Up to now, high-resolution mapping of surface water extent from satellites has only been available for a few regions, over limited time periods. The extension of the temporal and spatial coverage was difficult, due to the limitation of the remote sensing technique [e.g., the interaction of the radiation with vegetation or cloud for visible observations or the temporal sampling with the synthetic aperture radar (SAR)]. The advantages and the limitations of the various satellite techniques are reviewed. The need to have a global and consistent estimate of the water surfaces over long time periods triggered the development of a multi-satellite methodology to obtain consistent surface water all over the globe, regardless of the environments. The Global Inundation Extent from Multi-satellites (GIEMS) combines the complementary strengths of satellite observations from the visible to the microwave, to produce a low-resolution monthly dataset (\(0.25^\circ \,\times \,0.25^\circ\)) of surface water extent and dynamics. Downscaling algorithms are now developed and applied to GIEMS, using high-spatial-resolution information from visible, near-infrared, and synthetic aperture radar (SAR) satellite images, or from digital elevation models. Preliminary products are available down to 500-m spatial resolution. This work bridges the gaps and prepares for the future NASA/CNES Surface Water Ocean Topography (SWOT) mission to be launched in 2020. SWOT will delineate surface water extent estimates and their water storage with an unprecedented spatial resolution and accuracy, thanks to a SAR in an interferometry mode. When available, the SWOT data will be adopted to downscale GIEMS, to produce a long time series of water surfaces at global scale, consistent with the SWOT observations.     

Challenges and Opportunities in Water Cycle Research: WCRP Contributions

The state of knowledge and outstanding challenges and opportunities in global water cycle observations, research and modeling are briefly reviewed to set the stage for the reasons behind the new thrusts promoted by the World Climate Research Programme (WCRP) as Grand Challenges to be addressed on a 5- to 10-year time frame. Those focused on water are led by the GEWEX (Global Energy and Water Exchanges) project. A number of GEWEX science questions are being brought forward within GEWEX and the WCRP under guidance of the Joint Scientific Committee. Here, we describe what are some imperatives and opportunities for major advancements in observations, understanding, modeling and product development for water resources and climate that will enable a wide range of climate services and inform decisions on water resources management and practices.     

Detection of Continental Hydrology and Glaciology Signals from GRACE: A Review

Since its launch in March 2002, the Gravity Recovery and Climate Experiment (GRACE) has provided a global mapping of the time-variations of the Earth’s gravity field. Tiny variations of gravity from monthly to decadal time scales are mainly due to redistributions of water mass inside the surface fluid envelops of our planet (i.e., atmosphere, ocean and water storage on continents). In this article, we present a review of the major contributions of GRACE satellite gravimetry in global and regional hydrology. To date, many studies have focused on the ability of GRACE to detect, for the very first time, the time-variations of continental water storage (including surface waters, soil moisture, groundwater, as well as snow pack at high latitudes) at the unprecedented resolution of ~400–500 km. As no global complete network of surface hydrological observations exists, the advances of satellite gravimetry to monitor terrestrial water storage are significant and unique for determining changes in total water storage and water balance closure at regional and continental scales.     

Assessing Arkansas Ground Water for Pesticides: Methodology and Findings

During 1985 to 1987, 119 wells, Springs and municipal drinking water supplies throughout Arkansas were monitored for the presence of pesticides. Pesticides selected for analysis included acifluorfen, alachlor, aldicarb, atrazine, benomyl, cyanazine, cypermethrin, 2,4-D, dichlorprop, diuron, fenvalerate, fluometuron, hexazinone, linuron, metolachlor, permethrin, picloram, and propanil. Not every sample was analyzed for every pesticide. Overall, results indicated that the 18 herbicides, fungicides, and insecticides were not present in the ground water samples studied. (Note: Detectable concentrations of three herbicides – alachlor, atrazine, and metolachlor – were found in one irrigation well, at 5.5,5.8, and 6.9 μg/L, respectively. However, since previous and subsequent sampling failed to detect these compounds, their presence is attributed to a localized spill or handling error rather than agricultural application.)     

山西夏县东郭井水位水温巨幅异常变化调查与分析

山西夏县东郭井是山西省南部地区映震性能较好的一口流体观测井,自2007年观测以来,水位、水温的几次阶变异常与山西及周边地区的中等地震都有较好的对应关系。2011年4月16日,东郭井水位、水温再次出现阶变变化,且变化幅度较大,异常出现后,对观测系统、地下水开采情况、井孔结构三方面进行调查与分析,最终确定此次异常变化是由井孔井管错位造成,非地震前兆异常。