The terrestrial water resources in Shaanxi–Gansu–Ningxia (SGN) region are relatively scarce, and its climate change is unstable. Research on the deformation driven by terrestrial water load is of great significance to the dynamic maintenance of reference station networks. In this paper, data derived from Gravity Recovery and Climate Experiment (GRACE) and Global Navigation Satellite System (GNSS) from 2010 to 2014 were combined to monitor the spatiotemporal characteristics of surface vertical deformation caused by terrestrial water load change. The single scale factor was calculated by comparing CPC, WGHM, and GLDAS hydrological model to restore filtering leakage signal. The singular spectrum analysis (SSA) method was used to extract the principal component of temporal vertical deformation, and its spatial distribution was analyzed. At the same time, in order to study the relationship between the terrestrial water load deformation from GRACE and that from GNSS, the first-order term correction, the Atmosphere and Ocean De-aliasing Level-1B product (GAC) correction, and the first-order load LOVE number correction for GRACE were adopted in this paper. In addition, a quantitative comparative analysis of both the monitoring results was carried out. The results show that the time-variable characteristics of surface vertical deformation characterized by the filtered three hydrological models were consistent with those of GRACE. The correlation coefficient and Nash–Sutcliffe efficiency coefficient (NSE) values were the highest in the GLDAS model and the GRACE model, respectively; the former index is 0.93, while the latter is 0.85. The crustal vertical deformation from terrestrial water load showed a declining rate from 2010 to 2014. Its spatial change rate showed an obvious ladder distribution, with the surface subsidence rate gradually decreasing from south to north. In addition, weighted root mean square (WRMS) contribution rate of the crustal vertical deformation resulting from GRACE with GAC correction between the different GNSS stations ranged from 18.52% to 54.82%. The correlation coefficient between them was close to 0.70. After deducting the mass load impact of GRACE only, the WRMS contribution rate of the corresponding stations decreased from −8.42% to 21.18%. The correlation coefficient between them reduced noticeably. Adding GAC back can increase the comparability with GRACE and GNSS in terms of monitoring the crustal vertical deformation. The annual amplitude and phase of surface vertical deformation resulting from GRACE with GAC correction were close to those of GNSS. The research results can help to explore the motion mechanism between water migration and surface deformation, which is of benefit in the protection of the water ecological environment in the region.

Monitoring regional terrestrial water load deformation is of great significance to the dynamic maintenance and hydrodynamic study of the regional benchmark framework. In view of the lack of a spatial interpolation method based on the GNSS (Global Navigation Satellite System) elevation time series for obtaining terrestrial water load deformation information, this paper proposes to employ a CORS (Continuously Operating Reference Stations) network combined with environmental loading data, such as ECMWF (European Centre for Medium-Range Weather Forecasts) atmospheric data, the GLDAS (Global Land Data Assimilation System) hydrological model, and MSLA (Mean Sea Level Anomaly) data. Based on the load deformation theory and spherical harmonic analysis method, we took 38 CORS stations in southeast Zhejiang province as an example and comprehensively determined the vertical deformation of the crust as caused by regional terrestrial water load changes from January 2015 to December 2017, and then compared these data with the GRACE (Gravity Recovery and Climate Experiment) satellite. The results show that the vertical deformation value of the terrestrial water load in southeast Zhejiang, as monitored by the CORS network, can reach a centimeter, and the amplitude changes from −1.8 cm to 2.4 cm. The seasonal change is obvious, and the spatial distribution takes a ladder form from inland to coastal regions. The surface vertical deformation caused by groundwater load changes in the east–west–south–north–central sub-regions show obvious fluctuations from 2015 to 2017, and the trends of the five sub-regions are consistent. The amplitude of surface vertical deformation caused by groundwater load change in the west is higher than that in the east. We tested the use of GRACE for the verification of CORS network monitoring results and found a relatively consistent temporal distribution between both data sets after phase delay correction on GRACE, except for in three months—November in 2015, and January and February in 2016. The results show that the comprehensive solution based on the CORS network can effectively improve the monitoring of crustal vertical deformation during regional terrestrial water load change.

In this study, the basic concept of height nonlinear velocity field modeling in the CORS station is described. The noise results in a large deviation between the observation and predicted height. An ARCH testing method for heteroscedasticity of CORS height residual square series was proposed and the non-stationary characteristic of CORS height residual square time series was proved. A CORS height nonlinear velocity field reconstruction method based on the GARCH model was proposed. First, a nonlinear LS periodic fitting model was established for CORS height series data. Then, a GARCH model was established for the fitted non-stationary residual series. Finally, the signal term, linear trend term, and GARCH model noise term of nonlinear LS modeling were combined to reconstruct the nonlinear velocity field of the CORS height. The RMSE of nonlinear LS cycle modeling for 25 CORS stations worldwide ranged from 5 to 10 mm. The differences between the velocity, approximate annual and semi-annual amplitudes, and SOPAC results were 0.73 mm/a, 0.94 mm, and 0.51 mm, respectively. Compared with the centimeter amplitude of the CORS station height, the accuracy of the nonlinear model established in this study met the requirements. The results of height nonlinear velocity field reconstruction at 25 CORS stations worldwide showed that the mean square error of prediction of the one-year height movement reached 9 mm, and the average prediction accuracy of the semi-annual was 7 mm. Compared with the calculation accuracy of the current global CORS elevation component of 3–5 mm, the prediction error in this study was about 3 mm. The expected goal was achieved regarding the accuracy of the CORS station height nonlinear velocity field model.