Gravity Recovery and Climate Experiment (GRACE) satellites can effectively monitor terrestrial water storage (TWS) changes in large-scale areas. However, due to the coarse resolution of GRACE products, there is still a large number of deficiencies that need to be considered when investigating TWS changes in small-scale areas. Hence, it is necessary to downscale the GRACE products with a coarse resolution. First, in order to solve this problem, the present study employs modeling windows of different sizes (Window Size, WS) combined with multiple machine learning algorithms to develop a new machine learning spatial downscaling method (MLSDM) in the spatial dimension. Second, The MLSDM is used to improve the spatial resolution of GRACE observations from 0.5° to 0.25°, which is applied to Guantao County. The present study has verified the downscaling accuracy of the model developed through the combination of WS3, WS5, WS7, and WS9 and jointed with Random Forest (RF), Extra Tree Regressor (ETR), Adaptive Boosting Regressor (ABR), and Gradient Boosting Regressor (GBR) algorithms. The analysis shows that the accuracy of each combined model is improved after adding the residuals to the high-resolution downscaled results. In each modeling window, the accuracy of RF is better than that of ETR, ABR, and GBR. Additionally, compared to the changes in the TWS time series that are derived by the model before and after downscaling, the results indicate that the downscaling accuracy of WS5 is slightly more superior compared to that of WS3, WS7, and WS9. Third, the spatial resolution of the GRACE data was increased from 0.5° to 0.05° by integrating the WS5 and RF algorithm. The results are as follows: (1) The TWS (GWS) changes before and after downscaling are consistent, decreasing at −20.86 mm/yr and −21.79 mm/yr (−14.53 mm/yr and −15.46 mm/yr), respectively, and the Nash–Sutcliffe efficiency coefficient (NSE) and correlation coefficient (CC) values of both are above 0.99 (0.98). (2) The CC between the 80% deep groundwater well data and the downscaled GWS changes are above 0.70. Overall, the MLSDM can not only effectively improve the spatial resolution of GRACE products but also can preserve the spatial distribution of the original signal, which can provide a reference scheme for research focusing on the downscaling of GRACE products.

The launch of GRACE satellites has provided a new avenue for studying the terrestrial water storage anomalies (TWSA) with unprecedented accuracy. However, the coarse spatial resolution greatly limits its application in hydrology researches on local scales. To overcome this limitation, this study develops a machine learning-based fusion model to obtain high-resolution (0.25°) groundwater level anomalies (GWLA) by integrating GRACE observations in the North China Plain. Specifically, the fusion model consists of three modules, namely the downscaling module, the data fusion module, and the prediction module, respectively. In terms of the downscaling module, the GRACE-Noah model outperforms traditional data-driven models (multiple linear regression and gradient boosting decision tree (GBDT)) with the correlation coefficient (CC) values from 0.24 to 0.78. With respect to the data fusion module, the groundwater level from 12 monitoring wells is incorporated with climate variables (precipitation, runoff, and evapotranspiration) using the GBDT algorithm, achieving satisfactory performance (mean values: CC: 0.97, RMSE: 1.10 m, and MAE: 0.87 m). By merging the downscaled TWSA and fused groundwater level based on the GBDT algorithm, the prediction module can predict the water level in specified pixels. The predicted groundwater level is validated against 6 in-situ groundwater level data sets in the study area. Compare to the downscaling module, there is a significant improvement in terms of CC metrics, on average, from 0.43 to 0.71. This study provides a feasible and accurate fusion model for downscaling GRACE observations and predicting groundwater level with improved accuracy.

Groundwater is crucial for economic development in arid and semiarid areas. The Shiyang River Basin (SRB) has the most prominent water use issues in northwestern China, and overexploited groundwater resources have led to continuous groundwater-level decline. The key governance planning project of the SRB was issued in 2007. This paper synthetically combines remote-sensing data from Gravity Recovery and Climate Experiment (GRACE) data and precipitation, actual evapotranspiration, land use, and in situ groundwater-level data to evaluate groundwater storage variations on a regional scale. Terrestrial water storage anomalies (TWSA) and groundwater storage anomalies (GWSA), in addition to their influencing factors in the SRB since the implementation of the key governance project, are analyzed in order to evaluate the effect of governance. The results show that GRACE-derived GWS variations are consistent with in situ observation data in the basin, with a correlation coefficient of 0.68. The GWS in the SRB had a slow downward trend from 2003 to 2016, and this increased by 0.38 billion m³/year after 2018. As the meteorological data did not change significantly, the changes in water storage are mainly caused by human activities, which are estimated by using the principle of water balance. The decline in GWS in the middle and lower reaches of the SRB has been curbed since 2009 and has gradually rebounded since 2014. GWS decreased by 2.2 mm EWH (equivalent water height) from 2011 to 2016, which was 91% lower than that from 2007 to 2010. The cropland area in the middle and lower reaches of the SRB also stopped increasing after 2011 and gradually decreased after 2014, while the area of natural vegetation gradually increased, indicating that the groundwater level and associated ecology significantly recovered after the implementation of the project.

Drought monitoring is useful to minimize the impact of drought on human production and the natural environment. Gravity Recovery and Climate Experiment (GRACE) satellites can directly capture terrestrial water storage anomalies (TWSA) in the large basin, which represents a new source of hydrological information. In this study, the GRACE-based total storage deficit index (TSDI) is employed to investigate the temporal evolution and spatial distribution of drought in Southwest China from 2003 to 2016. The comparison results of TSDI with the standardized precipitation index (SPI), the standardized precipitation evapotranspiration index (SPEI), and the self-calibrating Palmer drought severity index (SC-PDSI) show that TSDI has significant consistency with them, which verifies the reliability of TSDI. The spatial distribution of TSDI was more consistent with the governmental drought reports than SC-PDSI in the most severe drought event from September 2009 to April 2010. Finally, the links between drought and climate indicators are investigated using the partial least square regression (PLSR) model. The results show that insufficient precipitation has the most significant impact on drought in Southwest China, followed by excessive evaporation. Although Southwest China is selected as a case study in this paper, the method can be applied in other regions as well.

Dense Global Position System (GPS) arrays can be used to invert the terrestrial water-storage anomaly (TWSA) with higher accuracy. However, the uneven distribution of GPS stations greatly limits the application of GPS to derive the TWSA. Aiming to solve this problem, we grid the GPS array using regression to raise the reliability of TWSA inversion. First, the study uses the random forest (RF) model to simulate crustal deformation in unobserved grids. Meanwhile, the new Machine-Learning Loading-Inverted Method (MLLIM) is constructed based on the traditional GPS derived method to raise the truthfulness of TWSA inversion. Second, this research selects southwest China as the study region, the MLLIM and traditional GPS inversion methods are used to derive the TWSA, and the inverted results are contrasted with datasets of the Gravity Recovery and Climate Experiment (GRACE) Mascon and the Global Land Data Assimilation System (GLDAS) model. The comparison shows that values of Pearson Correlation Coefficient (PCC) between the MLLIM and GRACE and GRACE Follow-On (GRACE-FO) are equal to 0.91 and 0.88, respectively; and the values of R-squared (R²) are equal to 0.76 and 0.65, respectively; the values of PCC and R² between MLLIM and GLDAS solutions are equal to 0.79 and 0.65. Compared with the traditional GPS inversion, the MLLIM improves PCC and R² by 8.85% and 7.99% on average, which indicates that the MLLIM can improve the accuracy of TWSA inversion more than the traditional GPS method. Third, this study applies the MLLIM to invert the TWSA in each province of southwest China and combines the precipitation to analyze the change of TWSA in each province. The results are as follows: (1) The spatial distribution of TWSA and precipitation is coincident, which is highlighted in southwest Yunnan and southeast Guangxi; (2) this study compares TWSA of MLLIM with GRACE and GLDAS solutions in each province, which indicates that the maximum value of PCC is as high as 0.86 and 0.94, respectively, which indicates the MLLIM can be used to invert the TWSA in the regions with sparse GPS stations. The TWSA based on the MLLIM can be used to fill the vacancy between GRACE and GRACE-FO.