Multi-Model Machine Learning Mapping of Gully Erosion Susceptibility in the Heihe Region of the Xiaoxingán Mountains, China

Gully erosion is a major driver of irreversible soil loss in Northeast China’s Mollisol belt, a region that supplies roughly one-quarter of the national grain output. Existing susceptibility assessments in this region have rarely combined multi-model comparison with spatially explicit cross-validation, and the predictive contribution of composite anthropogenic indicators such as the Human Footprint Index (HFI) has not been quantitatively benchmarked against conventional topographic variables. This study addresses these gaps for the Heihe region by combining an inventory of 4020 gully polygons supported by field checks in Xunke County, 16 VIF-screened environmental factors, three tree-based ensemble models and a logistic regression baseline. Under stratified random splitting, XGBoost achieved the highest discrimination (AUC = 0.95, κ = 0.74); under leave-one-district-out spatial cross-validation all tree-based models retained AUC above 0.83, confirming that random-split metrics overestimate discrimination by approximately 0.11 AUC units due to spatial autocorrelation and inter-district covariate shift. SHAP analysis identified LULC and HFI as the dominant predictors, exceeding all topographic variables, while slope gradient contributed least—consistent with the low-relief, intensively cultivated character of the study area. Susceptibility was highest in the southwestern agricultural lowlands. A one-factor sensitivity test in which only NDVI was increased by 20% suggested a reduction in modelled high-susceptibility area of approximately 12%, although co-occurring land-cover and hydrological changes were not simulated. The multi-model framework, integrating spatial cross-validation and post hoc interpretability, provides an explicit estimate of conventional evaluation optimism and supports spatially differentiated erosion management.