The previous comparative studies on watersheds were mostly based on the comparison of dispersive characteristics, which lacked systemicity and causality. We proposed a causal structure-based framework for basin comparison based on the Bayesian network (BN) and focus on the basin-scale water–energy–food–ecology (WEFE) nexus. We applied it to the Syr Darya River basin (SDB) and the Amu Darya River basin (ADB), of which poor water management caused the Aral Sea disaster. The causality of the nexus was effectively compared and universality of this framework was discussed. In terms of changes in the nexus, the sensitive factor for the water supplied to the Aral Sea changed from the agricultural development during the Soviet Union period to the disputes in the WEFE nexus after the disintegration. The water–energy contradiction of the SDB is more severe than that of the ADB, partly due to the higher upstream reservoir interception capacity. It further made management of the winter surplus water downstream of the SDB more controversial. Due to this, the water–food–ecology conflict between downstream countries may escalate and turn into a long-term chronic problem. Reducing water inflow to depressions and improving the planting structure prove beneficial to the Aral Sea ecology, and this effect of the SDB is more significant. The construction of reservoirs on the Panj River of the upstream ADB should be cautious to avoid an intense water–energy conflict such as the SDB's. It is also necessary to promote the water-saving drip irrigation and to strengthen the cooperation.

Climate change-induced precipitation variability is the leading cause of rainfall erosivity that leads to excessive soil losses in most countries of the world. In this paper, four global climate models (GCMs) were used to characterize the spatiotemporal prediction of rainfall erosivity and assess the effect of variations of rainfall erosivity in Central Asia. The GCMs (BCCCSM1-1, IPSLCM5BLR, MIROC5, and MPIESMLR) were statistically downscaled using the delta method under Representative Concentration Pathways (RCPs) 2.6 and 8.5 for two time periods: “Near” and “Far” future (2030s and 2070s). These GCMs data were used to estimate rainfall erosivity and its projected changes over Central Asia. WorldClim data was used as the present baseline precipitation scenario for the study area. The rainfall erosivity (R) factor of the Revised Universal Soil Loss Equation (RUSLE) was used to determine rainfall erosivity. The results show an increase in the future periods of the annual rainfall erosivity compared to the baseline. For all GCMs, with an average change in rainfall erosivity of about 5.6% (424.49 MJ mm ha⁻¹ h⁻¹ year⁻¹) in 2030s and 9.6% (440.57 MJ mm ha⁻¹ h⁻¹ year⁻¹) in 2070s as compared to the baseline of 402 MJ mm ha⁻¹ h⁻¹ year⁻¹. The magnitude of the change varies with the GCMs, with the largest change being 26.6% (508.85 MJ mm ha⁻¹ h⁻¹ year⁻¹), occurring in the MIROC-5 RCP8.5 scenario in the 2070s. Although annual rainfall erosivity shows a steady increase, IPSLCM5ALR (both RCPs and periods) shows a decrease in the average erosivity. Higher rainfall amounts were the prime causes of increasing spatial-temporal rainfall erosivity.