The robustness of landslide prediction models has become a major focus of researchers worldwide. We developed two novel hybrid predictive models that combine the self-organizing, deep-learning group method of data handling (GMDH) with two swarm intelligence optimization algorithms, i.e., cuckoo search algorithm (CSA) and whale optimization algorithm (WOA) for spatially explicit prediction of landslide susceptibility. Eleven landslide-causing factors and 334 historic landslides in a 31,340 km2 landslide-prone area in Iran were used to produce geospatial training and validation datasets. The GMDH model was employed to develop a basic predictive model that was then restructured and its parameters were optimized using the CSA and WOA algorithms, yielding the novel hybrid GMDH-CSA and GMDH-WOA models. The hybrid models were validated and compared to the standalone GMDH model by calculating the area under the receiver operating characteristic (AUC) curve and root mean square error (RMSE). The results demonstrated that the hybrid models overcame the computational shortcomings of the basic GMDH model and significantly improved landslide susceptibility prediction (GMDH-CSA, AUC = 0.909 and RMSE = 0.089; GMDH-WOA, AUC = 0.902 and RMSE = 0.129; standalone GMDH, AUC = 0.791 and RMSE = 0.226). Further, the hybrid models were more robust than the standalone GMDH model, showing consistently excellent performance when the training and validation datasets were changed. Overall, the swarm intelligence-optimized models, but not the standalone model, identified the best trade-offs among objectives, accuracy, and robustness.
