We developed a new boundary-included inversion model to improve the terrestrial water storage (TWS) inverted from regional GPS vertical deformation data. Through defining a new disc load empirical function (DLEF) and considering the mass change effect from the near but outside region, the result shows the TWS is more reasonable than the one inverted directly. Six simulation tests further confirmed the effectiveness of the boundary-included model. Finally, our new boundary-included model was used to derive the TWS in the Pacific Rim of the western United States based on the GPS-observed vertical deformation information. The inversion results show that our boundary-included inversion model can effectively improve the inversion results by 10–20% in terms of variance reduction in the boundary regions.

This paper aims to address the limitations of the distribution number and uniformity of Continuously Operating Reference Stations (CORS) and their impact on the reliability of inverting regional groundwater storage (GWS) based on Green’s function method and using global navigation satellite system (GNSS) data. A fusion method on the inversion of regional GWS changes from GNSS and the Gravity Recovery and Climate Experiment (GRACE) was proposed in this paper. Taking the Shaanxi–Gansu–Ningxia (SGN) region as an example, the in situ groundwater level data from ten CORS stations and eight wells were used for test analyses. In this paper, an atmospheric pressure model from the European Centre for Medium-Range Weather Forecasts (ECMWF), a global land data assimilation system (GLDAS), a WaterGAP global hydrology model (WGHM), and mean sea level anomaly (MSLA) data were used to quantitatively monitor the influence of vertical deformation caused by non-tidal environmental load. After deducing these loading deformations from the filtered time series of non-linear monthly geodetic height from the GNSS, the GWS changes in the SGN region from 2011 to 2014 were inverted. Meanwhile, the change in surface water storage from the GLDAS and WGHM models were removed from the terrestrial water storage (TWS) changes derived from GRACE. On this basis, the remove–restore theory in the Earth’s gravity field was introduced to both fuse the inversion results and obtain the regional GWS changes based on the fusion method. The results showed the following: (1) The local characteristics from the fusion results were more prominent than those of GRACE on the spatial scale, such as in the southwest and northeast in the study area. In addition, the fusion results were more uniform than those from GNSS, especially for the sparse and missing areas in which CORS stations were located, and the local effect was weakened. (2) On the time scale, compared with GRACE, the trends in GWS changes obtained from the fusion method and from GNSS inversion were roughly the same as the in situ groundwater level changes. (3) For the in situ groundwater wells “6105010031” and “6101260010”, the correlation coefficients of the fusion result were 0.53 and 0.56, respectively. The accuracy of the fusion method was slightly higher than that from GNSS, which indicates that the fusion method may be more effective for areas where CORS stations are missing or sparsely distributed. The methods in this paper can provide significant reference material for hydrodynamic research, sustainable management of water resources, and the dynamic maintenance of height data.

A global gravity field model (GGM) is essential to be validated with ground-based or airborne observational data for the accurate application of the GGM at a regional scale. Furthermore, accurately understanding the commission errors between the GGM and observational data are crucial for improving regional gravity fields. Taking the North China region as an example, to circumvent the omission errors, it is necessary to unify the spatial resolutions of the EIGEN-6C4 model and terrestrial gravity observational data to 110 km (determined by the distribution of gravity stations) by employing the spherical harmonic function for the EIGEN-6C4 model and the Slepian basis function for the gravity data, respectively. However, the application of spherical harmonic function expansions in the gravity model results in the Gibbs phenomenon, which may be a primary factor contributing to commission errors and impedes the accurate validation of the EIGEN-6C4 model with terrestrial gravity data. To effectively mitigate this issue, this study proposes a combination approach of window function filtering and regional eigenvalue constraint (based on the Slepian basis). Utilizing the EIGEN-6C4 gravity model to derive the gravity disturbance field at a resolution of 110 km (with spherical harmonic expansion up to the 180th degree and order), the combination approach effectively suppresses over 90% of high-degree (above the 120th degree) Gibbs phenomena. This approach also reduces signal leakage outside the region, thus enhancing the spatial accuracy of the regional gravity disturbance field. A subsequent comparison of the regional gravity disturbance field derived from the true model and terrestrial gravity data in North China indicates excellent consistency, with a root mean squared error (RMSE) of 0.80 mGal. This validation confirms that the combined approach of window function filtering and regional eigenvalue constraints effectively mitigates the Gibbs phenomenon and yields precise regional gravity fields. This approach is anticipated to significantly benefit scientific applications such as improving the accuracy of regional elevation benchmarks and accurately inverting the Earth’s internal structure.