Multidimensional Copula-Based Assessment, Propagation, and Prediction of Drought in the Lower Songhua River Basin

As global climate change intensifies, understanding drought mechanisms is crucial for managing water resources and agriculture. This study employs the Standardized Precipitation–Actual Evapotranspiration Index (SPAEI), Standardized Runoff Index (SRI), and Standardized Soil Moisture Index (SSMI) to analyze meteorological, hydrological, and agricultural droughts in the lower Songhua River basin. The PLUS model was used to predict future land types, with model accuracy validated using four evaluation metrics. The projected land cover was integrated with CMIP6 data into the SWAT model to simulate future runoff, which was used to calculate future SRI. Drought events were extracted using run theory, while drought occurrence probability and return period were calculated via a Copula-based joint distribution model. Bayesian conditional probability was employed to explore propagation mechanisms. The results indicate a significant increase in multidimensional drought risk, particularly when the cumulative frequency of univariate droughts reaches 25%, 50%, or 75%. Although increased duration and intensity enhance the likelihood of combined droughts, extremely high values cause a decline in joint probability under “OR” and “AND” conditions. Under different climate scenarios, the recurrence intervals of meteorological, hydrological, and agricultural droughts in the lower reaches of the Songhua River exhibit increased sensitivity with severity, demonstrating consistent propagation patterns across the meteorological–hydrological–agricultural system. Meteorological drought was found to propagate to hydrological and agricultural drought within ~6.00 months and ~3.67 months, respectively, with severity amplifying this effect. Propagation thresholds between drought types decreased with increasing intensity. This study combined SWAT and CMIP6 models with PLUS-based land-use scenarios, highlighting that land-use changes significantly influence spatiotemporal drought patterns. Model validation (Kappa = 0.83, OA = 0.92) confirmed robust predictive accuracy. Overall, this study proposes a multidimensional drought risk model integrating Copula and Bayesian networks, offering valuable insights for drought management under climate change.