Abstract
This study develops an interpretable and validated XGBoost–SHAP framework integrating multisource geospatial data and historical flood observations, with model performance evaluated using independent validation and future projections driven by bias-corrected CMIP6 climate scenarios to characterize the spatiotemporal variations and contributions of flood driving factors across three regions of Kazakhstan (2000–2025). The results demonstrate pronounced spatial differences in flood-driving factors: delayed snowmelt coupled with orographic rainfall dominates flood variability in the mountainous Almaty Region; hydrological memory effects regulate flood responses in the Akmola plains; and socioeconomic exposure shows an increasing contribution to flood risk evolution in Turkestan. Future multi-scenario simulations indicate that the flood-affected area in the Almaty Region is projected to increase by 10.8% under SSP2-4.5, which is associated with enhanced snowmelt processes, whereas the Akmola Region may experience a 23.2% reduction under SSP5-8.5, which is associated with changes in evaporation–soil moisture interactions. The interaction between socioeconomic development and natural hazards results in divergent risk trajectories: urban expansion in Akmola and Turkestan may offset declining hydroclimatic hazards, creating a potential risk paradox, whereas mountainous regions remain sensitive to concurrent increases in hazard intensity and exposure. These findings indicate that flood risk evolution in the studied regions of arid Central Asia is being increasingly influenced by socioeconomic dynamics in addition to natural hazards, highlighting the importance of differentiated adaptive planning strategies.