Improving the Resolution of GRACE/InSAR Groundwater Storage Estimations Using a New Subsidence Feature Weighted Combination Scheme

GRACE observations and land subsidence data derived from InSAR both assess groundwater storage changes. However, GRACE data at local scales are restricted by the coarser spatial resolution of satellite systems, and inversion of Groundwater Storage Anomalies (GWSA) by InSAR requires extensive and unavailable lithological data. Here, we propose a New Subsidence Feature Weighted Combination (NSFWC) scheme to enhance the spatial resolution of GRACE-derived GWSA from 0.5° to 0.05°. This method can not only retain the spatial distribution of groundwater changes but also reflect local details related to surface subsidence. A case study was executed to evaluate the performance of the NSFWC scheme in the Beijing Plain, which has seriously overexploited groundwater. Results showed that the simulated GWSA were consistent with in situ measurements in most regions, with a correlation coefficient of 0.85 and an RMSE of 4.41 mm/year. Additionally, there were 22 overexploited wells in the Beijing Plain, although groundwater levels generally recovered after the South to North Water Diversion Project. Simultaneously, four cones of depression were detected by the InSAR technology, where the maximum cumulative subsidence and subsidence rate achieved −198.52 mm and −53.09 mm/year, respectively. This paper provides data support and technical guarantees for small-scale groundwater resources management.