Estimation of Land Deformation and Groundwater Storage Dynamics in Shijiazhuang–Baoding–Cangzhou–Hengshui Using Multi-Temporal Interferometric Synthetic Aperture Radar

Groundwater resources are crucial to socio-economic development and the ecosystem, and over-extraction can cause the groundwater level to drop, deplete reserves, and trigger geological hazards like land subsidence. The North China Plain (NCP) has experienced both subsidence and groundwater depletion due to over-extraction in the past 70 years. In this study, we used MT-InSAR technology and ascending C-band Sentinel-1 SAR data from 2017 to 2023 to study land deformation in the junction area of Shijiazhuang–Baoding–Cangzhou–Hengshui. We identified multiple subsidence funnels with a maximum rate exceeding −150 mm/year and a total deformation surpassing 600 mm. Seasonal decomposition methods accurately separated seasonal signals in the time-series deformation and groundwater level data. An exponential function model applied to long-term deformation showed no significant decrease in subsidence in severely affected areas. By modeling seasonal deformation and seasonal groundwater levels, we determined the elastic skeletal storage coefficients (Ske) to be in the range of 1.02 × 10−3~6.53 × 10−3 in subsidence areas. We obtained the spatiotemporal evolution of the total groundwater storage (TGWS), irreversible ground storage (IGWS), and recoverable ground storage (RGWS). The TGWS and IGWS decreased annually while the RGWS increased, which is attributable to the implementation of the South-to-North Water Diversion Project (SNWDP) and the issuance of groundwater withdrawal policies in the NCP.