The West Liaohe River Basin (WLRB) is one of the most sensitive areas to climate change in China and an important grain production base in the Inner Mongolia Autonomous Region of China. Groundwater depletion in this region is becoming a critical issue. Here, we used the Gravity Recovery and Climate Experiment (GRACE) satellite data and in situ well observations to estimate groundwater storage (GWS) variations and discussed the driving factors of GWS changes in the WLRB. GRACE detects a GWS decline rate of −0.92 ± 0.49 km³/yr in the WLRB during 2005–2011, consistent with the estimate from in situ observations (−0.96 ± 0.19 km³/yr). This long-term GWS depletion is attributed to reduced precipitation and extensive groundwater overexploitation in the 2000s. Long-term groundwater level observations and reconstructed total water storage variations since 1980 show favorable agreement with precipitation anomalies at interannual timescales, both of which are significantly influenced by the El Niño-Southern Oscillation (ENSO). Generally, the WLRB receives more/less precipitation during the El Niño/La Niña periods. One of the strongest El Niño events on record in 1997–1998 and a subsequent strong La Niña drastically transform the climate of WLRB into a decade-long drought period, and accelerate the groundwater depletion in the WLRB after 1998. This study demonstrates the significance of integrating satellite observations, ground-based measurements, and climatological data for interpreting regional GWS changes from a long-term perspective.

Evapotranspiration (ET) is usually difficult to estimate at the regional scale due to scarce direct measurements. This study uses the water balance equation to calculate the regional ET with observations of precipitation, runoff, and terrestrial water storage changes (TWSC) in nine exorheic catchments of China. We compared the regional ET estimates from a water balance perspective with and without considering TWSC (ET_WB: ET estimates with considering TWSC, and ET_PQ: ET estimates from precipitation minus runoff without considering TWSC). Results show that the regional annual ET ranges from 417.7 mm/yr to 831.5 mm/yr in the nine exorheic catchments based on the water balance equation. The impact of ignoring TWSC on calculating ET is notable, as the root mean square errors (RMSEs) of annual ET between ET_WB and ET_PQ range from 12.0–105.8 mm/yr (2.6–12.7% in corresponding annual ET) among the exorheic catchments. We also compared the estimated regional ET with other ET products. Different precipitation products are assessed to explain the inconsistency between different ET products and regional ET from a water balance perspective. The RMSEs between ET estimates from Gravity Recovery and Climate Experiment (GRACE) and ET from land surface models can be reduced if the deviation of precipitation forcing data is considered. ET estimates from Global Land Evaporation Amsterdam Model (GLEAM) can be improved by reducing the uncertainty of precipitation forcing data in three semiarid catchments. This study emphasizes the importance of considering TWSC when calculating the regional ET using a water balance equation and provides more accurate ET estimates to help improve modeled ET results.

Terrestrial water storage (TWS) can be influenced by both climate change and anthropogenic activities. While the Gravity Recovery and Climate Experiment (GRACE) satellites have provided a global view on long-term trends in TWS, our ability to disentangle human impacts from natural climate variability remains limited. Here we present a quantitative method to isolate these two contributions with reconstructed climate-driven TWS anomalies (TWSA) based on long-term precipitation data. Using the Haihe River Basin (HRB) as a case study, we find a higher human-induced water depletion rate (−12.87 ± 1.07 mm/yr) compared to the original negative trend observed by GRACE alone for the period of 2003–2013, accounting for a positive climate-driven TWSA trend (+4.31 ± 0.72 mm/yr). We show that previous approaches (e.g., relying on land surface models) provide lower estimates of the climate-driven trend, and thus likely underestimated the human-induced trend. The isolation method presented in this study will help to interpret observed long-term TWS changes and assess regional anthropogenic impacts on water resources.

The terrestrial water storage anomaly (TWSA) gap between the Gravity Recovery and Climate Experiment (GRACE) and its follow-on mission (GRACE-FO) is now a significant issue for scientific research in high-resolution time-variable gravity fields. This paper proposes the use of singular spectrum analysis (SSA) to predict the TWSA derived from GRACE. We designed a case study in six regions in China (North China Plain (NCP), Southwest China (SWC), Three-River Headwaters Region (TRHR), Tianshan Mountains Region (TSMR), Heihe River Basin (HRB), and Lishui and Wenzhou area (LSWZ)) using GRACE RL06 data from January 2003 to August 2016 for inversion, which were compared with Center for Space Research (CSR), Helmholtz-Centre Potsdam-German Research Centre for Geosciences (GFZ), Jet Propulsion Laboratory (JPL)’s Mascon (Mass Concentration) RL05, and JPL’s Mascon RL06. We evaluated the accuracy of SSA prediction on different temporal scales based on the correlation coefficient (R), Nash–Sutcliffe efficiency (NSE), and root mean square error (RMSE), which were compared with that of an auto-regressive and moving average (ARMA) model. The TWSA from September 2016 to May 2019 were predicted using SSA, which was verified using Mascon RL06, the Global Land Data Assimilation System model, and GRACE-FO results. The results show that: (1) TWSA derived from GRACE agreed well with Mascon in most regions, with the highest consistency with Mascon RL06 and (2) prediction accuracy of GRACE in TRHR and SWC was higher. SSA reconstruction improved R, NSE, and RMSE compared with those of ARMA. The R values for predicting TWS in the six regions using the SSA method were 0.34–0.98, which was better than those for ARMA (0.26–0.97), and the RMSE values were 0.03–5.55 cm, which were better than the 2.29–5.11 cm RMSE for ARMA as a whole. (3) The SSA method produced better predictions for obvious periodic and trending characteristics in the TWSA in most regions, whereas the detailed signal could not be effectively predicted. (4) The predicted TWSA from September 2016 to May 2019 were basically consistent with Global Land Data Assimilation System (GLDAS) results, and the predicted TWSA during June 2018 to May 2019 agreed well with GRACE-FO results. The research method in this paper provides a reference for bridging the gap in the TWSA between GRACE and GRACE-FO.