Understanding the spatiotemporal characteristics of hydrological components and their impacts on vegetation are critical for comprehending hydrological, climatological, and ecological processes under environmental change and solving future water management challenges. Innovative methods need to be developed in semiarid areas to analyze the special hydrological factors in the water resource systems of these areas. Gravity Recovery and Climate Experiment (GRACE) and Global Land Data Assimilation System (GLDAS) were applied with the normalized difference vegetation index (NDVI) data in this paper to analyze spatiotemporal changes of hydrological factors in the Xiliaohe River Basin (XRB). The results showed that precipitation (P), evapotranspiration (ET) and temperature (T) had similar seasonal change patterns at rates of 0.05 cm/yr., 0.01 cm/yr. and −0.05 °C/yr., respectively. Total water storage change (TWSC) was consistent with the change trend of soil moisture change (SMC) and showed a fluctuating trend. Groundwater change (GWC) showed a decreasing trend at a rate of −0.43 cm/yr. P and ET had a greater impact on GLDAS data (R = 0.634, P < 0.05 and R = 0.686, P < 0.01, respectively) than on other factors. GWC was more sensitive to changes in T (R = 0.570, P < 0.05). Furthermore, a lag period of 0 to 1 months was observed for the effects of P and ET on TWSC and GLDAS. NDVI showed an upward trend at a rate of 0.001 yr⁻¹ between 2002 and 2014. A spatial distribution of NDVI was heterogeneous in the study area. ET, GLDAS and GWC in growing season limited vegetation growth and were more important than other factors in XRB. The results may contribute to an understanding of the relationships between the hydrological cycle and climate change and provide scientific support for local environmental management.

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.

Climate change remains among the most pressing environmental challenges confronting the world, exerting profound pressure on both ecological systems and socio-economic development. To advance understanding of the evolution patterns and driving mechanisms governing hydroclimatic systems in arid and semi-arid regions, this study employed an integrated framework encompassing trend testing, change-point detection, periodicity and persistence analysis, and machine learning-based attribution. Focusing on the Mu Us Sandy Land from 1982 to 2023, we systematically investigated the spatiotemporal evolution, periodic characteristics, and driving mechanisms of hydroclimatic factors. Furthermore, future climate risks were assessed using CMIP6 multi-model data. The results showed that: (1) All four variables exhibited positive slopes, but only soil moisture showed a statistically significant long-term wetting trend (β = 0.025 × 10⁻³, p = 0.0008) and a clear global abrupt change in 2011; the upward tendencies of precipitation (p = 0.3946), potential evapotranspiration (p = 0.4970), and surface runoff (p = 0.1097) did not reach the 0.05 significance level. (2) Meteorological elements showed weak periodicity and strong anti-persistence (mean Hurst index = 0.379 for precipitation and 0.222 for PET), whereas hydrological elements exhibited clear seasonal–interannual periods and more random future variability with greater spatial heterogeneity (mean Hurst index = 0.436 for runoff and 0.414 for soil moisture). (3) Monthly changes were mainly associated with local surface processes. Vegetation dynamics were key predictors of precipitation, runoff, and soil moisture, while potential evapotranspiration was dominated by atmospheric demand, with limited influence from large-scale climate indices. (4) Under high-emission scenarios, imbalanced water–heat increases may lead to a higher likelihood of drought conditions.

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.