Climate change and land use change are the two main factors affecting the regional water cycle and water resources management. However, runoff studies in the karst basin based on future scenario projections are still lacking. To fill this gap, this study proposes a framework consisting of a future land use simulation model (FLUS), an automated statistical downscaling model (ASD), a soil and water assessment tool (SWAT) and a multi-point calibration strategy. This frameword was used to investigate runoff changes under future climate and land use changes in karst watersheds. The Chengbi River basin, a typical karst region in southwest China, was selected as the study area. The ASD method was developed for climate change projections based on the CanESM5 climate model. Future land use scenarios were projected using the FLUS model and historical land use data. Finally, the SWAT model was calibrated using a multi-site calibration strategy and was used to predict future runoff from 2025–2100. The results show that: (1) the developed SWAT model obtained a Nash efficiency coefficient of 0.83, which can adequately capture the spatial heterogeneity characteristics of karst hydro-climate; (2) land use changes significantly in all three future scenarios, with the main phenomena being the interconversion of farmland and grassland in SSPs1-2.6, the interconversion of grassland, farmland and artificial surfaces in SSPs2-4.5 and the interconversion of woodland, grassland and artificial surfaces in SSPs5-8.5; (3) the average annual temperature will show an upward trend in the future, and the average annual precipitation will increase by 11.53–14.43% and (4) the future annual runoff will show a significant upward trend, with monthly runoff mainly concentrated in July–September. The variability and uncertainty of future runoff during the main-flood period may increase compared to the historical situation. The findings will benefit future water resources management and water security in the karst basin.

Runoff prediction plays an extremely important role in flood prevention, mitigation, and the efficient use of water resources. Machine learning runoff prediction models have become popular due to their high computational efficiency. To select a model with a better runoff simulation and to validate the stability of the model, the following studies were done. Firstly, the support vector machine Model (SVM), the Elman Neural Network Model (ENN), and the multi-model mean model (MMM) were used for the runoff prediction, with the monthly runoff data from 1963–2007 recorded by the Pingtang hydrological station in the Chengbi River Karst Basin, China. Secondly, the comprehensive rating index method was applied to select the best model. Thirdly, the indicators of the hydrologic alteration–range of variability approach (IHA-RVA) was introduced to measure the model stability with different data structure inputs. According to the comprehensive rating index method, the SVM model outperformed the other models and was the best runoff prediction model with a score of 0.53. The overall change of the optimal model was 10.52%, which was in high stability.

Accurate hydrological simulation plays an important role in the research of hydrological problems; the accuracy of the watershed hydrological model is seriously affected by model-parameter uncertainty and model-input uncertainty. Thus, in this study, different calibration methods and rainfall inputs were introduced into the SWAT (Soil and Water Assessment Tool) model for watershed hydrological simulation. The Chengbi River basin, a typical karst basin in Southwest China, was selected as the target basin. The indicators of the NSE (Nash efficiency coefficient), Re (relative error) and R² (coefficient of determination) were adopted to evaluate the model performance. The results showed that: on the monthly and daily scales, the simulated runoff with the single-site method calibrated model had the lowest NSE value of 0.681 and highest NSE value of 0.900, the simulated runoff with the multi-site method calibrated model had the lowest NSE value of 0.743 and highest NSE value of 0.953, increased correspondingly, indicating that adopting the multi-site method could reduce the parameter uncertainty and improve the simulation accuracy. Moreover, the NSE values with IMERG (Integrated Multisatellite Retrievals for Global Rainfall Measurement) satellite rainfall data were the lowest, 0.660 on the monthly scale and 0.534 on the daily scale, whereas the NSE values with fusion rainfall data processed by the GWR (geographical weighted regression) method greatly increased to 0.854 and 0.717, respectively, and the NSE values with the measured rainfall data were the highest, 0.933 and 0.740, respectively, demonstrating that the latter two rainfall inputs were more suitable sources for hydrological simulation.

Drought poses a significant constraint on economic development. Drought assessment using the standardized precipitation index (S_PI) uses only precipitation data, eliminating other redundant and complex calculation processes. However, the sparse stations in southwest China and the lack of information on actual precipitation measurements make drought assessment highly dependent on satellite precipitation data whose accuracy cannot be guaranteed. Fortunately, the Chengbi River Basin in Baise City is rich in station precipitation data. In this paper, based on the evaluation of the accuracy of IMERG precipitation data, geographically weighted regression (GWR), geographic difference analysis (GDA), and cumulative distribution function (CDF) are used to fuse station precipitation data and IMERG precipitation data, and finally, the fused precipitation data with the highest accuracy are selected to evaluate the drought situation. The results indicate that the accuracy of IMERG precipitation data needs to be improved, and the quality of CDF-fused precipitation data is higher than the other two. The drought analysis indicated that the Chengbi River Basin is in a cyclical drought and flood situation, and from October to December 2014, the S_PI was basically between +1 and −1, showing a spatial pattern of slight flooding, normal conditions, and slight drought.

Karst basins have a relatively low capacity for water retention, rendering them very vulnerable to drought hazards. However, karst geo-climatic features are highly spatially heterogeneous, making reliable drought assessment challenging. To account for geo-climatic heterogeneous features and to enhance the reliability of drought assessment, a framework methodology is proposed. Firstly, based on the history of climate (1963–2019) from the Global Climate Model (GCM) and station observations within the Chengbi River karst basin, a multi-station calibration-based automated statistical downscaling (ASD) model is developed, and the Kling–Gupta efficiency (KGE) and Nash–Sutcliffe efficiency (NSE) are selected as performance metrics. After that, future climate (2023–2100) under three GCM scenarios (SSP1-2.6, SSP2-4.5, and SSP5-8.5) are obtained by using the ASD model. Finally, the Standardized Precipitation Evapotranspiration Index (SPEI), calculated by future climate is applied to assess drought conditions. The results indicate that the multi-station calibration-based ASD model has good performance and thus can be used for climate data downscaling in karst areas. Precipitation mainly shows a significant upward trend under all scenarios with the maximum variation (128.22%), while the temperature shows a slow upward trend with the maximum variation (3.44%). The drought condition in the 2040s is still relatively severe. In the 2060s and 2080s, the basin is wetter compared with the historical period. The percentage of drought duration decreases in most areas from the 2040s to the 2080s, demonstrating that the future drought condition is alleviated. From the SSP1-2.6 scenario to the SSP5-8.5 scenario, the trend of drought may also increase.

The complex and unique topography of karst regions highlights the weaknesses of traditional hydrological models which fail to fully generalize them. The successive proposals of standard artificial intelligence (AI) models, pre-processing techniques, and post-processing methods have provided new opportunities to enhance the accuracy of runoff prediction in karst areas. In this study, first, the BP neural network model and the Elman neural network model were used for runoff prediction. Then, the performance of four coupled models—formed by combining two AI pre-processing techniques, Empirical Modal Decomposition (EMD) and Ensemble Empirical Modal Decomposition (EEMD), with the previously mentioned AI models—was investigated. Finally, the accuracy of triple-coupled models, formed by applying the post-processing method of quantile mapping (QM) to the previous coupled models, was estimated. The Nash–Sutcliffe efficiency (NSE), the mean absolute percentage error (MAPE), the root mean square error (RMSE), and the peak percentage of threshold statistics (PPTS) were selected to evaluate and analyze the forecasting results of the above models. The results demonstrated that the BP model had the best prediction effect of the standard AI models, the coupled forecasting models had better prediction accuracy than the standard AI models, and the triple-coupled QM–EMD–Elman model had the best forecasting effect with an NSE value of 0.73, MAPE value of 0.75, RMSE value of 34.60, and PPTS value of 2.36.