Strong ground motions usually trigger lots of slope failures in the affected area. In this work, we analyse the occurrence likelihood of earthquake-triggered landslide by employing the ensembles of adaptive neuro-fuzzy inference systems (ANFIS) with four well-known metaheuristics techniques, namely particle swarm optimization (PSO), genetic algorithm (GA), ant colony optimization (ACO), and differential evolution (DE) algorithms. Twelve landslide conditioning factors namely, elevation, slope degree, lithology, peak ground acceleration (PGA), stream power index (SPI), topographic wetness index (TWI), distance to road, distance to river, distance to fault, normalized difference vegetation index (NDVI), slope aspect, and plan curvature are considered within the geographic information system (GIS) to produce the required spatial database. In this paper, frequency ratio (FR) model is used to evaluate the spatial interaction between the landslides and conditioning factors. Meantime, among a total of 458 marked earthquake-induced landslides, 366 (80%) are specified to the learning process, and the remaining 92 (20%) landslides are used to evaluate the accuracy of applied models. The landslide susceptibility maps are generated in the GIS environment. Three accuracy criteria of mean square error (MSE), root mean square error (RMSE), and area under the receiving operating characteristic curve (AUROC) are used to develop a ranking system for comparing the integrity of the designed models. The total ranking scores (TRSs) of 15, 8, 10, and 18, respectively, obtained for PSO-ANFIS, GA-ANFIS, ACO-ANFIS, and DE-ANFIS revealed the superiority of the DE algorithm compared to other metaheuristics techniques. Also, the DE-ANFIS emerged as the fastest ensemble, due to the highest convergence speed obtained for this model.
