The increasingly severe phenomenon of groundwater drought poses a dual threat to the development and construction of a region, as well as its ecological environment. Traditional groundwater drought monitoring methods rely on observation wells, which makes it difficult to obtain dynamic drought information in areas with limited measurement data. Based on Gravity Recovery and Climate Experiment (GRACE) satellite technology and data, the suitability of the standardized groundwater index (GRACE_SGI) was explored for drought characterization in the Mu Us Sandy Land. Multiscale and seasonal trend changes in groundwater drought in the study area from 2002 to 2021 were comprehensively identified. Subsequently, the characteristics of hysteresis time between the GRACE_SGI and the standardized precipitation index (SPI) were clarified. The results show that (1) different fitting functions impact the parameterized GRACE_SGI fitting results. The Anderson–Darling method was used to find the best-fitting function for groundwater data in the study area: the Pearson III distribution. (2) The gain and loss characteristics of the GRACE_SGI are similar, showing downward trends at different time scales, including seasonal scales. (3) The absolute values based on the maximum correlation coefficients between the SPI and the GRACE_SGI at different time scales were 0.1296, 0.2483, 0.2427, and 0.5224, with time lags of 0, 0, 12, and 11 months, respectively. The vulnerability of semiarid ecosystems to hydroclimatic changes is highlighted by these findings, and a satellite-based framework for monitoring groundwater drought in data-scarce regions is provided.

Taking the Yellow River Basin (YRB) as an example, this study explores the non-stationary drought evolution features in large river basins under climate change. This study utilized precipitation and multiple climate factor data to establish the non-stationary standardized precipitation index (NSPI) through the GAMLSS model. Combined with the run theory, Copula function and a cascaded RF-LSTM machine learning model, the drought characteristics and retrospective predictive patterns were systematically assessed. The results show that: (1) The Arctic Oscillation, the Pacific Decadal Oscillation, the Southern Oscillation and the North Pacific Index are the primary climate drivers of non-stationary precipitation variation in the YRB, with the former three being selected most frequently and NPI additionally influencing April–June and September, and their effects are both different and lagging. Compared with the traditional SPI, the NSPI assigned higher drought grades and greater severity to typical drought years (e.g., the 1974 event was rated D3 with a severity of 17.935 by NSPI versus D2 with 11.733 by SPI), and thus better captured non-stationary drought evolution. (2) The duration of droughts exhibited a decreasing trend that was not statistically significant (p > 0.05), whereas drought intensity and severity decreased significantly (p < 0.05); the peak severity showed a significant upward trend (p = 0.0078). Spatially, the northwest of the Loess Plateau was a compound core area with high severity, high frequency and long duration of droughts, while the upper reaches were mainly characterized by low severity, short duration and sudden droughts. (3) The drought risk in the YRB shows a higher frequency in the lower reaches and a lower frequency in the upper reaches. The middle and lower reaches were high-risk areas, with shorter AND-type joint exceedance return periods for moderate drought (2.46–5.83 years) and severe drought (3.77–9.15 years). The upper reaches were low-risk areas, with longer return periods reaching up to 5.83 years for moderate drought and 9.15 years for severe drought. The study shows that the NSPI, considering the driving of multiple climate factors, can more effectively identify and assess non-stationary drought risks, providing a scientific basis for drought prevention and control in river basins.