Terrestrial water storage (TWS) is a pivotal component of the global water cycle, profoundly impacting water resource management, hazard monitoring, and agriculture production. The Gravity Recovery and Climate Experiment (GRACE) and its successor, the GRACE Follow-On (GFO), have furnished comprehensive monthly TWS data since April 2002. However, there are 35 months of missing data over the entire GRACE/GFO observational period. To address this gap, we developed an operational approach utilizing singular spectrum analysis and principal component analysis (SSA-PCA) to fill these missing data over mainland China. The algorithm was demonstrated with good performance in the Southwestern River Basin (SWB, correlation coefficient, CC: 0.71, RMSE: 6.27 cm), Yangtze River Basin (YTB, CC: 0.67, RMSE: 3.52 cm), and Songhua River Basin (SRB, CC: 0.66, RMSE: 7.63 cm). Leveraging two decades of continuous time-variable gravity data, we investigated the spatiotemporal variations in TWS across ten major Chinese basins. According to the results of GRACE/GFO, mainland China experienced an average annual TWS decline of 0.32 ± 0.06 cm, with the groundwater storage (GWS) decreasing by 0.54 ± 0.10 cm/yr. The most significant GWS depletion occurred in the Haihe River Basin (HRB) at −2.07 ± 0.10 cm/yr, significantly substantial (~1 cm/yr) depletions occurred in the Yellow River Basin (YRB), SRB, Huaihe River Basin (HHB), Liao-Luan River Basin (LRB), and Southwest River Basin (SWB), and moderate losses were recorded in the Northwest Basin (NWB, −0.34 ± 0.03 cm/yr) and Southeast River Basin (SEB, −0.24 ± 0.10 cm/yr). Furthermore, we identified that interannual TWS variations in ten basins of China were primarily driven by soil moisture water storage (SMS) anomalies, exhibiting consistently and relatively high correlations (CC > 0.60) and low root-mean-square errors (RMSE < 5 cm). Lastly, through the integration of GRACE/GFO and Global Land Data Assimilation System (GLDAS) data, we unraveled the contrasting water storage patterns between northern and southern China. Southern China experienced drought conditions, while northern China faced flooding during the 2020–2023 La Niña event, with the inverse pattern observed during the 2014–2016 El Niño event. This study fills in the missing data and quantifies water storage variations within mainland China, contributing to a deeper insight into climate change and its consequences on water resource management.

Accurately monitoring spatio-temporal changes in lake water levels is important for studying the impacts of climate change on freshwater resources, and for predicting natural hazards. In this study, we applied multi-mission radar satellite altimetry data from the Laurentian Great Lakes, North America to optimally reconstruct multi-decadal lake-wide spatio-temporal changes of water level. We used the results to study physical processes such as teleconnections of El Niño and southern oscillation (ENSO) episodes over approximately the past three-and-a-half decades (1985–2018). First, we assessed three reconstruction methods, namely the standard empirical orthogonal function (EOF), complex EOF (CEOF), and complex independent component analysis (CICA), to model the lake-wide changes of water level. The performance of these techniques was evaluated using in-situ gauge data, after correcting the Glacial Isostatic Adjustment (GIA) process using a contemporary GIA forward model. While altimeter-measured water level was much less affected by GIA, the averaged gauge-measured water level was found to have increased up to 14 cm over the three decades. Our results indicate that the CICA-reconstructed 35-year lake level was more accurate than the other two techniques. The correlation coefficients between the CICA reconstruction and the in situ water-level data were 0.96, 0.99, 0.97, 0.97, and 0.95, for Lake Superior, Lake Michigan, Lake Huron, Lake Erie, and Lake Ontario, respectively; ~7% higher than the original altimetry data. The root mean squares of errors (RMSE) were 6.07 cm, 4.89 cm, 9.27 cm, 7.71 cm, and 9.88 cm, respectively, for each of the lakes, and ~44% less than differencing with the original altimetry data. Furthermore, the CICA results indicated that the water-level changes in the Great Lakes were significantly correlated with ENSO, with correlation coefficients of 0.5–0.8. The lake levels were ~25 cm higher (~30 cm lower) than normal during EI Niño (La Niña) events.