Droughts caused by runoff are an important environmental issue in the context of global climate change, with profound impacts on ecosystems, agriculture and water resource management. To assess the impact of future climate change on the hydrological response of watersheds, this study combines the SWAT (Soil and Water Assessment Tool) and MODFLOW (MODular groundwater FLOW model) models to predict future changes in runoff and hydrological drought in watersheds using data from two scenarios under 15 CMIP6 climate models. The results show that: (1) The R² and NSE values of monthly runoff at the Caizuzi station in the Naoli River basin are greater than 0.60 in different periods; (2) the ensemble of climate models after screening can effectively improve the accuracy of runoff simulation and reduce the prediction uncertainty of a single climate model; (3) under different scenarios, the temperature generally increases, the precipitation increases and evapotranspiration increased under the SSP2-4.5 scenario and decreased under the SSP5-8.5 scenario; (4) runoff showed an increasing trend under the SSP2-4.5 scenario and the opposite trend under the SSP5-8.5 scenario; (5) the frequency of winter runoff droughts decreased in the future period, while the frequency of spring and summer droughts increased, with the change trend being more pronounced under the SSP5-8.5 scenario; (6) compared with the baseline period (1965–2014), under the SSP2-4.5 and SSP5-8.5 scenarios, the average annual temperature in the watershed increased by 1.89 °C and 3.22 °C, respectively, and the annual precipitation increased by 32% and 36.19%, respectively, but the summer and autumn runoff decreased; and (7) The SRI-3 model analysis indicates that hydrological droughts will significantly intensify under both future emission scenarios. Under the SSP5-8.5 scenario, droughts will worsen earlier and the abrupt change will occur earlier, while under the SSP2-4.5 scenario, although the abrupt change will occur later, the drought intensity will be higher. The critical drought transition periods are 2030–2047 (SSP5-8.5) and 2045–2055 (SSP2-4.5). This study provides important scientific basis for adaptive water resources management and drought mitigation strategies in cold-region watersheds under future climate scenarios.

Climate and land use changes significantly affect runoff and hydrological drought, presenting challenges for water resource management. This study focuses on the Naoli River Basin, utilizing the SWAT model integrated with PLUS land use projections under the CMIP6 SSP245 and SSP585 scenarios to assess trends in runoff and drought characteristics from 2025 to 2100. The Standardized Runoff Index (SRI) and run theory are applied to analyze drought frequency and duration. Key findings include the following: (1) Under the SSP585 scenario (2061–2100), land use changes—specifically, a reduction in cropland and an increase in forest cover—resulted in a 12.59% decrease in runoff compared to the baseline period (1970–2014), with notable differences when considering climate-only scenarios. (2) The SSP585 scenario exhibits a significant rise in drought frequency and duration, particularly during summer, whereas SSP245 shows milder trends. (3) Based on the Taylor plot evaluation, the ensemble average MMM-Best (r = 0.80, RMSE = 26.15) has been identified as the optimal prediction model for the 2025–2100 period. Deviation analysis revealed that NorESM2-MM and IPSL-CM6A-LR demonstrated the greatest stability, while EC-Earth3 exhibited the largest deviation and highest uncertainty. (4) Land use changes under the SSP245 scenario help mitigate drought by enhancing water retention, although their effectiveness diminishes under SSP585 due to the dominant influence of climate factors, including increased temperature and precipitation variability. And (5) SRI-3 mutation analysis indicated that the mutation point occurred in July 2074 under the SSP245 scenario and in April 2060 under the SSP585 scenario (p < 0.05). The trend for SSP245 revealed significant fluctuations, with the number of crossover points rising to 40 following land use changes; conversely, the SSP585 trend remained stable with only seven crossover points, as high-emission scenarios predominantly influenced early mutations. These findings illuminate the interactive effects of land use and climate change, providing a scientific foundation for optimizing water resource management and developing effective drought mitigation strategies.

This study assesses future agricultural drought risk in the Ganjiang River Basin under climate change and land use change. A coupled analysis framework was established using the SWAT hydrological model, the CMIP6 climate models (SSP1-2.6, SSP2-4.5, SSP5-8.5), and the PLUS land use simulation model. Key methods included the Standardized Soil Moisture Index (SSMI), travel time theory for drought event identification and duration analysis, Mann–Kendall trend test, and the Pettitt change-point test to examine soil moisture dynamics from 2027 to 2100. The results indicate that the CMIP6 ensemble performs excellently in temperature simulations, with a correlation coefficient of R² = 0.89 and a root mean square error of RMSE = 1.2 °C, compared to the observational data. The MMM-Best model also performs well in precipitation simulations, with R² = 0.82 and RMSE = 15.3 mm, compared to observational data. Land use changes between 2000 and 2020 showed a decrease in forestland (−3.2%), grassland (−2.8%), and construction land (−1.5%), with an increase in water (4.8%) and unused land (2.7%). Under all emission scenarios, the SSMI values fluctuate with standard deviations of 0.85 (SSP1-2.6), 1.12 (SSP2-4.5), and 1.34 (SSP5-8.5), with the strongest drought intensity observed under SSP5-8.5 (minimum SSMI = −2.8). Drought events exhibited spatial and temporal heterogeneity across scenarios, with drought-affected areas ranging from 25% (SSP1-2.6) to 45% (SSP5-8.5) of the basin. Notably, abrupt changes in soil moisture under SSP5-8.5 occurred earlier (2045–2050) due to intensified land use change, indicating strong human influence on hydrological cycles. This study integrated the CMIP6 climate projections with high-resolution human activity data to advance drought risk assessment methods. It established a framework for assessing agricultural drought risk at the regional scale that comprehensively considers climate and human influences, providing targeted guidance for the formulation of adaptive water resource and land management strategies.

With the intensification of global climate change, extreme weather events increasingly threaten water resources and agricultural systems. This study focuses on the Naoli River Basin, employing the Standardized Precipitation Actual Evapotranspiration Index (SPAEI), the Standardized Runoff Index (SRI), and the Standardized Surface Moisture Index (SSMI) to assess the spatiotemporal variability of meteorological, hydrological, and agricultural droughts. Drought events are identified based on travel time theory, and joint distributions of drought characteristics are modeled using optimized two- and three-dimensional copula functions. Lagged correlation and Bayesian conditional probability analyses are used to explore drought propagation processes. Key findings include (1) the SWAT model showed strong runoff simulation performance (R² > 0.75, NSE > 0.97), while the PLUS model achieved high land use simulation accuracy (overall accuracy > 0.93, Kappa > 0.85); (2) future projections suggest continued forest expansion and farmland decline, with water areas increasing under SSP245 and urban areas expanding under SSP585; (3) five CMIP6 models with high skill (r = 0.80, RMSE = 26.15) were selected via a Taylor diagram for scenario simulation; (4) copula-based joint drought probabilities vary temporally, with meteorological drought risks increasing under long-term moderate-emission scenarios, while hydrological and agricultural droughts show contrasting trends; (5) and under extreme meteorological drought, the conditional probability of extreme agricultural drought doubles from 0.12 (SSP245) to 0.24 (SSP585), indicating heightened vulnerability under high-emission pathways. These results offer critical insights for regional drought risk assessment and adaptive management under future climate scenarios.