GIS-Based Landslide Susceptibility Mapping with a Blended Ensemble Model and Key Influencing Factors in Sentani, Papua, Indonesia

Landslides represent a recurrent hazard in tropical mountain environments, where rapid urbanization and extreme rainfall amplify disaster risk. The Sentani region of Papua, Indonesia, is highly vulnerable, as demonstrated by the catastrophic debris flows of March 2019 that caused fatalities and widespread losses. This study developed high-resolution landslide susceptibility maps for Sentani using an ensemble machine learning framework. Three base learners—Random Forest, eXtreme Gradient Boosting (XGBoost), and CatBoost—were combined through a logistic regression meta-learner. Predictor redundancy was controlled using Pearson correlation and Variance Inflation Factor/Tolerance (VIF/TOL). The landslide inventory was constructed from multitemporal satellite imagery, integrating geological, topographic, hydrological, environmental, and seismic factors. Results showed that lithology, Slope Length and Steepness Factor (LS Factor), and earthquake density consistently dominated model predictions. The ensemble achieved the most balanced predictive performance, Area Under the Curve (AUC) > 0.96, and generated susceptibility maps that aligned closely with observed landslide occurrences. SHapley Additive Explanations (SHAP) analyses provided transparent, case-specific insights into the directional influence of key factors. Collectively, the findings highlight both the robustness and interpretability of ensemble learning for landslide susceptibility mapping, offering actionable evidence to support disaster preparedness, land-use planning, and sustainable development in Papua.