Temporal and Spatial Dynamics of Groundwater Drought Based on GRACE Satellite and Its Relationship with Agricultural Drought

Terrestrial water storage includes soil water storage, groundwater storage, surface water storage, snow water equivalent, plant canopy water storage, biological water storage, etc., which can comprehensively reflect the total change in water volume during processes such as precipitation, evapotranspiration, runoff, and human water use in the basin hydrological cycle. The Gravity Recovery and Climate Experiment (GRACE) satellite provides a powerful tool and a new approach for observing changes in terrestrial water storage and groundwater storage. The North China Plain (NCP) is a major agricultural region in the northern arid area of China, and long-term overexploitation of groundwater has led to increasingly prominent ecological vulnerability issues. This study uses GRACE and Global Land Data Assimilation System (GLDAS) hydrological model data to assess the spatiotemporal patterns of groundwater drought in the NCP and its various sub-regions from 2003 to 2022, identify the locations, occurrence probabilities, and confidence intervals of seasonal and trend mutation points, quantify the complex interactive effects of multiple climate factors on groundwater drought, and reveal the propagation time from groundwater drought to agricultural drought. The results show that: (1) from 2003 to 2022, the linear tendency rate of groundwater drought index (GDI) was −0.035 per 10 years, indicating that groundwater drought showed a gradually worsening trend during the study period; (2) on an annual scale, the most severe groundwater drought occurred in 2021 (GDI = −1.59). In that year, the monthly average GDI in the NCP ranged from −0.58 to −2.78, and the groundwater drought was most severe in July (GDI = −2.02); (3) based on partial wavelet coherence, the best univariate, bivariate for groundwater drought were soil moisture (PASC = 19.13%); and (4) in Beijing, Tianjin and Hebei, the propagation time was mainly concentrated in 1–5 months, with average lag times of 2.87, 3.20, and 2.92 months, respectively. This study can not only reduce and mitigate the harm of groundwater drought to agricultural production, social life, and ecosystems by monitoring changes in groundwater storage, but also provide a reference for the quantitative identification of the dominant factors of groundwater drought.