Application of the group method of data handling (GMDH) approach for landslide susceptibility zonation using readily available spatial covariates

Abstract Landslide susceptibility (LS) mapping is an essential tool for landslide risk assessment. This study aimed to provide a new approach with better performance for landslide mapping and adopting readily available variables. In addition, it investigates the capability of a state-of-the-art model developed using the group method of data handling (GMDH) to spatially model LS. Furthermore, hybridized models of GMDH were developed using different metaheuristic algorithms. The study area was the Bonghwa region of South Korea, for which an accurate landslide inventory dataset is available. We considered a total of 13 spatial covariates (altitude, slope, aspect, topographic wetness index, valley depth, plan curvature, profile curvature, distance from fault, distance from river, distance from road, land use, density of forest, and lithology were chosen as independent variables). Two benchmark models—random forest and boosted regression trees—were used to compare their results with the standalone GMDH and hybridized models. We compared model accuracy using the two most robust evaluation metrics, root mean square error (RMSE) and area under the receiver operating characteristic curve (AUROC). The validation results showed that hybridized models outperformed the standalone GMDH model. Moreover, the hybridized GMDH-PSO (AUC = 0.83, RMSE = 0.108), GMDH-IWO (AUC = 0.81, RMSE = 0.111), GMDH-BBO (AUC = 0.8; RMSE = 0.12), and GMDH-ICA (AUC = 0.8; RMSE = 0.117) had a better predictive performance than both RF and BRT. Therefore, the proposed approach could successfully produce landslide susceptibility maps using relatively few readily available variables and can be repeated in data-scarce regions.