Machine learning (ML) is increasingly utilized in Landslide Susceptibility Mapping (LSM), though challenges remain in interpreting the predictions of ML models. To reveal the response relationship between landslide susceptibility and evaluation factors, an interpretability model was constructed to analyze how the results of the ML model are realized. This study focuses on Zhenba County in Shaanxi Province, China, employing both Random Forest (RF) and Support Vector Machine (SVM) to develop LSM models optimized through Random Search (RS). To enhance interpretability, the study incorporates techniques such as Partial Dependence Plot (PDP), Local Interpretable Model-Agnostic Explanations (LIMEs), and Shapley Additive Explanations (SHAP). The RS-optimized RF model demonstrated superior performance, achieving an Area Under the Curve (AUC) of 0.965. The interpretability model identified the NDVI and distance from road as important factors influencing landslides occurrence. NDVI plays a positive role in the occurrence of landslides in this region, and the landslide-prone areas are within 500 m from the road. These analyses indicate the importance of improved hyperparameter selection in enhancing model accuracy and performance. The interpretability model provides valuable insights into LSM, facilitating a deeper understanding of landslide formation mechanisms and guiding the formulation of effective prevention and control strategies.

The reliability of data-driven approaches in generating landslide susceptibility maps depends on data quality, analytical method selection, and sampling techniques. Selecting optimal datasets and determining the most effective analytical methods pose significant challenges. This study assesses the performance of seven machine learning classifiers in the Himalayan region of the China–Pakistan Economic Corridor, utilizing statistical techniques and validation metrics. Thirteen geo-environmental variables were analyzed, including topographic (8), land cover (1), hydrological (1), geological (2), and meteorological (1) factors. These variables were evaluated for multicollinearity, feature importance, and their influence on landslide incidences. Our findings indicate that Support Vector Machines and Logistic Regression were highly effective, particularly near fault zones and roads, due to their effectiveness in handling complex, non-linear terrain interactions. Conversely, Random Forest and Logistic Regression demonstrated variability in their results. Each model distinctly identified landslide susceptibility zones ranging from very low to very high risk. Significant conditioning variables such as elevation, rainfall, lithology, slope, and land use were identified, reflecting the unique geomorphological conditions of the Himalayas. Further analysis using the Variance Inflation Factor and Pearson correlation coefficient showed minimal multicollinearity among the variables. Moreover, evaluations of Area Under the Receiver Operating Characteristic Curve (AUC-ROC) values confirmed the strong predictive capabilities of the models, with the Random Forest Classifier performing exceptionally well, achieving an AUC of 0.96 and an F-Score of 0.86. This study shows the importance of model selection based on dataset characteristics to enhance decision-making and strategy effectiveness.

The effectiveness of data-driven landslide susceptibility mapping relies on data integrity and advanced geospatial analysis; however, selecting the most suitable method and identifying key regional factors remains a challenging task. To address this, this study assessed the performance of six machine learning models, including Convolutional Neural Networks (CNNs), Random Forest (RF), Categorical Boosting (CatBoost), their CNN-based hybrid models (CNN+RF and CNN+CatBoost), and a Stacking Ensemble (SE) combining CNN, RF, and CatBoost in mapping landslide susceptibility along the Karakoram Highway in northern Pakistan. Twelve geospatial factors were examined, categorized into Topography/Geomorphology, Land Cover/Vegetation, Geology, Hydrology, and Anthropogenic Influence. A detailed landslide inventory of 272 occurrences was compiled to train the models. The proposed stacking ensemble and hybrid models improve landslide susceptibility modeling, with the stacking ensemble achieving an AUC of 0.91. Hybrid modeling enhances accuracy, with CNN–RF boosting RF’s AUC from 0.85 to 0.89 and CNN–CatBoost increasing CatBoost’s AUC from 0.87 to 0.90. Chi-square (χ²) values (9.8–21.2) and p-values (<0.005) confirm statistical significance across models. This study identifies approximately 20.70% of the area as from high to very high risk, with the SE model excelling in detecting high-risk zones. Key factors influencing landslide susceptibility showed slight variations across the models, while multicollinearity among variables remained minimal. The proposed modeling approach reduces uncertainties, enhances prediction accuracy, and supports decision-makers in implementing effective landslide mitigation strategies.

A rational landslide susceptibility mapping (LSM) can minimize the losses caused by landslides and enhance the efficiency of disaster prevention and reduction. At present, frequency ratio (FR), information value (IV), and certainty factor (CF) are widely used to quantify the relationships between landslides and their causative factors; however, it remains unclear which method is the most effective. Moreover, existing landslide susceptibility zoning methods lack full automation; thus, the results are full of uncertainties. To address this, the FR, IV, and CF were used to analyze the relationship between landslides and causative factors. Subsequently, three distinct sets of models were developed, namely random forest models (RF_FR, RF_IV, and RF_CF), support vector machine models (SVM_FR, SVM_IV, and SVM_CF), and logistic regression models (LR_FR, LR_IV, and LR_CF) using the analysis results as inputs. A Gaussian mixture model (GMM) was introduced as a new method for landslide susceptibility zoning, classifying the LSM into five distinct levels. An accuracy evaluation of the models and a rationality analysis of the LSM indicated that the FR is superior to the IV and CF in quantifying the relationship between landslides and causative factors. Additionally, the quantile method was employed as a comparative approach to the GMM, further validating the effectiveness of the GMM. This research contributes to more effective and efficient LSM, ultimately enhancing landslide prevention measures.