Landslide Susceptibility Assessment Using Spatial Multi-Criteria Evaluation Model in Rwanda

Slide type landslide susceptibility assessment of the Büyük Menderes watershed using artificial neural network method

The Büyük Menderes watershed is the largest drainage watershed in Western Anatolia with an area of approximately 26,000 km². In the study area, almost 863 landslides occurred, extending over 222 km² with a mean landslide area of 0.21 km². In this study, landslide susceptibility assessments were carried out using artificial neural network method, which is one of the data-driven methods. In this study, that will contribute to the mitigation or control of the landslides caused by the reasons controlling the spatial and temporal distribution of landslides created in the GIS and MATLAB environment by using scientific and technological approaches within the framework. Since derivative activation function is also used in back-propagation artificial neural networks, its derivative is easily calculated in order not to slow down the calculation. Levenberg-Marquardt back-propagation (LM), resilient back propagation back-propagation (trainrp), scaled conjugate gradient back-propagation (trainscg), conjugate gradient with Powell/Beale restarts back-propagation (traincgb), and Fletcher-Powell conjugate gradient back-propagation (traincgf) algorithms are used, which constantly interrogate the link between the input parameter and the result output, and at least one cell's output is given as an input to any other cell. Geology, digital elevation model, slope, topographic wetness index, roughness index, plan, profile curvatures, and proximity to active faults and rivers were used as landslide conditioning factors. In susceptibility assessments, landslides were separated by 70% analysis, 15% test, and 15% validation datasets by random selection method. The performances of the landslide susceptibility maps were assessed by the area under the ROC curve (AUC), accuracy (ACC), precision, recall, F1 score, Kappa test error histogram, and confusion matrix, respectively. The area under the receiver operating characteristic curves, analysis, testing, validation, landslides, and study areas were found between 0.873 and 0.911. The susceptibility map had a high prediction rate in which high and very high susceptible zones corresponded to 26% of the study area including 82% of the recorded landslides.

Impact and a Novel Representation of Spatial Data Uncertainty in Debris Flow Susceptibility Analysis

In a study of debris flow susceptibility on the European continent, an analysis of the impact between known location and a location accuracy offset for 99 debris flows demonstrates the impact of uncertainty in defining appropriate predisposing factors and consequent analysis for areas of susceptibility. The dominant predisposing environmental factors, as determined through Maximum Entropy modeling, are presented and analyzed with respect to the values found at debris flow event points versus a buffered distance of locational uncertainty around each point. Maximum Entropy susceptibility models are developed utilizing the original debris flow inventory of points, randomly generated points, and two models utilizing a subset of points with an uncertainty of 5 km, 1 km, and a model utilizing only points with a known location of “exact”. The AUCs are 0.891, 0.893, 0.896, 0.921, and 0.93, respectively. The “exact” model, with the highest AUC, is ignored in final analyses due to the small number of points and localized distribution, and hence susceptibility results are likely non-representational of the continent. Each model is analyzed with respect to the AUC, highest contributing factors, factor classes, susceptibility impact, and comparisons of the susceptibility distributions and susceptibility value differences. Based on model comparisons, geographic extent, and the context of this study, the models utilizing points with a location uncertainty of less than or equal to 5 km best represent debris flow susceptibility for the continent of Europe. A novel representation of the uncertainty is expressed and included in a final susceptibility map, as an overlay of standard deviation and mean of susceptibility values for the two best models, providing additional insight for subsequent action.

A Scientometric Analysis of Predicting Methods for Identifying the Environmental Risks Caused by Landslides

The effects of a landslide can represent a very big problem, including the death of people, damage to the land, environmental pollution and the loss of natural resources. Landslides are the most important medium for transferring sediments and polluting waterways by earth and organic materials. An excess of sediments reduces the quality of fish habitat and the potability of water. In order to understand landslides in depth, a thorough study was conducted using a scientometric analysis, as well as a thorough practical examination of landslide analysis and monitoring techniques. This review focused on methods used for landslide analysis, including physical models requiring easily prepared event-based landslide inventory, probabilistic methods which are useful for both shallow and earthquake-based landslides, and landslide monitoring performed by remote sensing techniques, which provide data helpful for prediction, monitoring and mapping. The fundamental principles of each method are described in terms of the method used, and its advantages, and limits. People and infrastructure are at danger from landslides caused by heavy rain, so this report highlights landslide-prone regions and considers the analysis methods for landslides used in these countries, with a view to identifying mitigation measures for coping with landslide risks in hilly areas. Furthermore, future landslide research possibilities, as well as possible modeling methods, are addressed. The report summarizes some landslide prediction and monitoring techniques used in landslide-prone countries which can help inform researchers seeking to protect the public from danger in landslide areas.

Nationwide Susceptibility Mapping of Landslides in Kenya Using the Fuzzy Analytic Hierarchy Process Model

Landslide susceptibility mapping (LSM) is a cost-effective tool for landslide hazard mitigation. To date, no nationwide landslide susceptibility maps have been produced for the entire Kenyan territory. Hence, this work aimed to develop a landslide susceptibility map at the national level in Kenya using the fuzzy analytic hierarchy process method. First, a hierarchical evaluation index system containing 10 landslide contributing factors and their subclasses was established to produce a susceptibility map. Then, the weights of these indexes were determined through pairwise comparisons, in which triangular fuzzy numbers (TFNs) were employed to scale the relative importance based on the opinions of experts. Ultimately, these weights were merged in a hierarchical order to obtain the final landslide susceptibility map. The entire Kenyan territory was divided into five susceptibility levels. Areas with very low susceptibility covered 5.53% of the Kenyan territory, areas with low susceptibility covered 20.58%, areas with the moderate susceptibility covered 29.29%, areas with high susceptibility covered 29.16%, and areas with extremely high susceptibility covered 15.44% of Kenya. The resulting map was validated using an inventory of 425 historical landslides in Kenya. The results indicated that the TFN-AHP model showed a significantly improved performance (AUC = 0.86) compared with the conventional AHP (AUC = 0.72) in LSM for the study area. In total, 31.53% and 29.88% of known landslides occurred within the “extremely high” and “high” susceptibility zones, respectively. Only 8.24% and 1.65% of known landslides fell within the “low” and “very low” susceptibility zones, respectively. The map obtained as a result of this study is beneficial to inform planning and land resource management in Kenya.

Quantifying Land Use Land Cover Changes in the Lake Victoria Basin Using Satellite Remote Sensing: The Trends and Drivers between 1985 and 2014

The Lake Victoria Basin (LVB) is a significant resource for five states within East Africa, which faces major land use land cover changes that threaten ecosystem integrity and ecosystem services derived from the basin’s resources. To assess land use land cover changes between 1985 and 2014, and subsequently determine the trends and drivers of these changes, we used a series of Landsat images and field data obtained from the LVB. Landsat image pre-processing and band combinations were done in ENVI 5.1. A supervised classification was applied on 118 Landsat scenes using the maximum likelihood classifier in ENVI 5.1. The overall accuracy of classified images was computed for the 2014 images using 124 reference data points collected through stratified random sampling. Computations of area under various land cover classes were calculated between the 1985 and 2014 images. We also correlated the area from natural vegetation classes to farmlands and settlements (urban areas) to explore relationships between land use land cover conversions among these classes. Based on our land cover classifications, we obtained overall accuracy of 71% and a moderate Kappa statistic of 0.56. Our results indicate that the LVB has undergone drastic changes in land use land cover, mainly driven by human activities that led to the conversion of forests, woodlands, grasslands, and wetlands to either farmlands or settlements. We conclude that information from this work is useful not only for basin-scale assessments and monitoring of land cover changes but also for targeting, prioritizing, and monitoring of small scale, community led efforts to restore degraded and fragmented areas in the basin. Such efforts could mitigate the loss of ecosystem services previously derived from large contiguous land covers which are no longer tenable to restore. We recommend adoption of a basin scale, operational, Earth observation-based, land use change monitoring framework. Such a framework can facilitate rapid and frequent assessments of gains and losses in specific land cover classes and thus focus strategic interventions in areas experiencing major losses, through mitigation and compensatory approaches.

Rainfall-Induced Landslide Prediction Using Machine Learning Models: The Case of Ngororero District, Rwanda

Landslides fall under natural, unpredictable and most distractive disasters. Hence, early warning systems of such disasters can alert people and save lives. Some of the recent early warning models make use of Internet of Things to monitor the environmental parameters to predict the disasters. Some other models use machine learning techniques (MLT) to analyse rainfall data along with some internal parameters to predict these hazards. The prediction capability of the existing models and systems are limited in terms of their accuracy. In this research paper, two prediction modelling approaches, namely random forest (RF) and logistic regression (LR), are proposed. These approaches use rainfall datasets as well as various other internal and external parameters for landslide prediction and hence improve the accuracy. Moreover, the prediction performance of these approaches is further improved using antecedent cumulative rainfall data. These models are evaluated using the receiver operating characteristics, area under the curve (ROC-AUC) and false negative rate (FNR) to measure the landslide cases that were not reported. When antecedent rainfall data is included in the prediction, both models (RF and LR) performed better with an AUC of 0.995 and 0.997, respectively. The results proved that there is a good correlation between antecedent precipitation and landslide occurrence rather than between one-day rainfall and landslide occurrence. In terms of incorrect predictions, RF and LR improved FNR to 10.58% and 5.77% respectively. It is also noted that among the various internal factors used for prediction, slope angle has the highest impact than other factors. Comparing both the models, LR model’s performance is better in terms of FNR and it could be preferred for landslide prediction and early warning. LR model’s incorrect prediction rate FNR = 9.61% without including antecedent precipitation data and 3.84% including antecedent precipitation data.

Landslide Susceptibility Prediction Modeling Based on Remote Sensing and a Novel Deep Learning Algorithm of a Cascade-Parallel Recurrent Neural Network

Landslide susceptibility prediction (LSP) modeling is an important and challenging problem. Landslide features are generally uncorrelated or nonlinearly correlated, resulting in limited LSP performance when leveraging conventional machine learning models. In this study, a deep-learning-based model using the long short-term memory (LSTM) recurrent neural network and conditional random field (CRF) in cascade-parallel form was proposed for making LSPs based on remote sensing (RS) images and a geographic information system (GIS). The RS images are the main data sources of landslide-related environmental factors, and a GIS is used to analyze, store, and display spatial big data. The cascade-parallel LSTM-CRF consists of frequency ratio values of environmental factors in the input layers, cascade-parallel LSTM for feature extraction in the hidden layers, and cascade-parallel full connection for classification and CRF for landslide/non-landslide state modeling in the output layers. The cascade-parallel form of LSTM can extract features from different layers and merge them into concrete features. The CRF is used to calculate the energy relationship between two grid points, and the extracted features are further smoothed and optimized. As a case study, the cascade-parallel LSTM-CRF was applied to Shicheng County of Jiangxi Province in China. A total of 2709 landslide grid cells were recorded and 2709 non-landslide grid cells were randomly selected from the study area. The results show that, compared with existing main traditional machine learning algorithms, such as multilayer perception, logistic regression, and decision tree, the proposed cascade-parallel LSTM-CRF had a higher landslide prediction rate (positive predictive rate: 72.44%, negative predictive rate: 80%, total predictive rate: 75.67%). In conclusion, the proposed cascade-parallel LSTM-CRF is a novel data-driven deep learning model that overcomes the limitations of traditional machine learning algorithms and achieves promising results for making LSPs.

Landslide Susceptibility Evaluation and Management Using Different Machine Learning Methods in The Gallicash River Watershed, Iran

This analysis aims to generate landslide susceptibility maps (LSMs) using various machine learning methods, namely random forest (RF), alternative decision tree (ADTree) and Fisher’s Linear Discriminant Function (FLDA). The results of the FLDA, RF and ADTree models were compared with regard to their applicability for creating an LSM of the Gallicash river watershed in the northern part of Iran close to the Caspian Sea. A landslide inventory map was created using GPS points obtained in a field analysis, high-resolution satellite images, topographic maps and historical records. A total of 249 landslide sites have been identified to date and were used in this study to model and validate the LSMs of the study region. Of the 249 landslide locations, 70% were used as training data and 30% for the validation of the resulting LSMs. Sixteen factors related to topographical, hydrological, soil type, geological and environmental conditions were used and a multi-collinearity test of the landslide conditioning factors (LCFs) was performed. Using the natural break method (NBM) in a geographic information system (GIS), the LSMs generated by the RF, FLDA, and ADTree models were categorized into five classes, namely very low, low, medium, high and very high landslide susceptibility (LS) zones. The very high susceptibility zones cover 15.37% (ADTree), 16.10% (FLDA) and 11.36% (RF) of the total catchment area. The results of the different models (FLDA, RF, and ADTree) were explained and compared using the area under receiver operating characteristics (AUROC) curve, seed cell area index (SCAI), efficiency and true skill statistic (TSS). The accuracy of models was calculated considering both the training and validation data. The results revealed that the AUROC success rates are 0.89 (ADTree), 0.92 (FLDA) and 0.97 (RF) and predication rates are 0.82 (ADTree), 0.79 (FLDA) and 0.98 (RF), which justifies the approach and indicates a reasonably good landslide prediction. The results of the SCAI, efficiency and TSS methods showed that all models have an excellent modeling capability. In a comparison of the models, the RF model outperforms the boosted regression tree (BRT) and ADTree models. The results of the landslide susceptibility modeling could be useful for land-use planning and decision-makers, for managing and controlling the current and future landslides, as well as for the protection of society and the ecosystem.

Stability Analysis of Plant-Root-Reinforced Shallow Slopes along Mountainous Road Corridors Based on Numerical Modeling

Engineering methods such as soil nails, geosynthetic reinforcement, retaining structures, gabions, and shotcrete are implemented to stabilize road cut slopes along mountainous areas. However, these structures are not environmentally friendly and, particularly in Ethiopia, it is impossible to address all road problems due to financial limitations. Nowadays, soil reinforcement with plant roots is recognized as an environmentally sustainable alternative to improve shallow slope failure along mountainous transportation corridors. The aims of this study was, therefore, to conduct slope stability analysis along a road corridor by incorporating the effect of plant roots. Five plant species were selected for the analysis based on their mechanical characteristics. Namely, Eucalyptus globules (tree), Psidium guajava (shrub), Salix subserrata (shrub), Chrysopogon zizanioides, and Pennisetum macrourum (grasses). The roots’ tensile strength and soil parameters were determined through tensile strength testing and triaxial compression tests, respectively. The factor of safety of the slope was calculated by the PLAXIS-2D software. The study showed that when the slope was reinforced with plant roots, the factor of safety (FOS) improved from 22–34%. The decreasing effect of vegetation on slope stability was observed when soil moisture increased. The sensitivity analysis also indicated that: (1) as the spacing between plants decreased, the effect of vegetation on the slope increased. (2) Slope angle modification with a combination of plant roots had a significant impact on slope stabilization. Of the five-selected plant species, Salix subserrata was the promising plant species for slope stabilization as it exhibited better root mechanical properties among selected plant species.

Land Use Change Impacts on Water Erosion in Rwanda

Rwanda has experienced accelerated soil erosion as a result of unsustainable human activities and changes in land use. Therefore, this study aimed at applying the RUSLE (Revised Universal Soil Loss Equation) model using GIS (Geographical Information System) and remote sensing to assess water erosion in Rwanda, focusing on the erosion-prone lands for the time span 2000 to 2015. The estimated mean annual soil losses were 48.6 t ha−1 y−1 and 39.2 t ha−1 y−1 in 2000 and 2015, respectively, resulting in total nationwide losses of approximately 110 and 89 million tons. Over the 15 years, 34.6% of the total area of evaluated LULC (land use/land cover) types have undergone changes. The highest mean soil loss of 91.6 t ha−1 y−1 occurred in the area changing from grassland to forestland (0.5%) while a mean soil loss of 10.0 t ha−1 y−1 was observed for grassland converting to cropland (4.4%). An attempt has been made to identify the embedded driving forces of soil erosion in Rwanda. As a result, we found that mean soil loss for Rwanda’s districts in 2015 was significantly correlated with poverty (r = 0.45, p = 0.013), increased use of chemical fertilizers (r = 0.77, p = 0.005), and especially was related to extreme poverty (r = 0.77, p = 0.000). The soil conservation scenario analysis for Rwanda’s cropland in 2015 revealed that terracing could reduce the soil loss by 24.8% (from 14.6 t ha−1 y−1 to 11.7 t ha−1 y−1). Most importantly, the study suggests that (1) terracing integrated with mulching and cover crops could effectively control water erosion while ameliorating soil quality and fertility, and (2) reforestation schemes targeting the rapid-growing tree species are therefore recommended as an important feature for erosion control in the study area.

Comparative Analysis of Deterministic and Semiquantitative Approaches for Shallow Landslide Risk Modeling in Rwanda

The use of appropriate approaches to produce risk maps is critical in landslide disaster management. The aim of this study was to investigate and compare the stability index mapping (SINMAP) and the spatial multicriteria evaluation (SMCE) models for landslide risk modeling in Rwanda. The SINMAP used the digital elevation model in conjunction with physical soil parameters to determine the factor of safety. The SMCE method used six layers of landslide conditioning factors. In total, 155 past landslide locations were used for training and model validation. The results showed that the SMCE performed better than the SINMAP model. Thus, the receiver operating characteristic and three statistical estimators-accuracy, precision, and the root mean square error (RMSE)-were used to validate and compare the predictive capabilities of the two models. Therefore, the area under the curve (AUC) values were 0.883 and 0.798, respectively, for the SMCE and SINMAP. In addition, the SMCE model produced the highest accuracy and precision values of 0.770 and 0.734, respectively. For the RMSE values, the SMCE produced better prediction than SINMAP (0.332 and 0.398, respectively). The overall comparison of results confirmed that both SINMAP and SMCE models are promising approaches for landslide risk prediction in central-east Africa.

Landslide Susceptibility Prediction Using Particle-Swarm-Optimized Multilayer Perceptron: Comparisons with Multilayer-Perceptron-Only, BP Neural Network, and Information Value Models

Landslides are one type of serious geological hazard which cause immense losses of local life and property. Landslide susceptibility prediction (LSP) can be used to determine the spatial probability of landslide occurrence in a certain area. It is important to implement LSP for landslide hazard prevention and reduction. This study developed a particle-swarm-optimized multilayer perceptron (PSO-MLP) model for LSP implementation to overcome the drawbacks of the conventional gradient descent algorithm and to determine the optimal structural parameters of MLP. Shicheng County in Jiangxi Province of China was used as the study area. In total, 369 landslides, randomly selected non-landslides, and 14 landslide-related predisposing factors were used to train and test the present PSO-MLP model and three other comparative models (an MLP-only model with the gradient descent algorithm, a back-propagation neural network (BPNN), and an information value (IV) model). The results showed that the PSO-MLP model had the most accurate prediction performance (area under the receiver operating characteristic curve (AUC) of 0.822 and frequency ratio (FR) accuracy of 0.856) compared with the MLP-only (0.791 and 0.829), BPNN (0.800 and 0.840), and IV (0.788 and 0.824) models. It can be concluded that the proposed PSO-MLP model addresses the drawbacks of the MLP-only model well and performs better than conventional artificial neural networks (ANNs) and statistical models. The spatial probability distribution law of landslide occurrence in Shicheng County was well revealed by the landslide susceptibility map produced using the PSO-MLP model. Furthermore, the present PSO-MLP model may have higher prediction and classification performances in some other fields compared with conventional ANNs and statistical models.

Impact of Spatial Resolution of Digital Elevation Model on Landslide Susceptibility Mapping: A case Study in Kullu Valley, Himalayas

Landslides are one of the most damaging geological hazards in mountainous regions such as the Himalayas. The Himalayan region is, tectonically, the most active region in the world that is highly vulnerable to landslides and associated hazards. Landslide susceptibility mapping (LSM) is a useful tool for understanding the probability of the spatial distribution of future landslide regions. In this research, the landslide inventory datasets were collected during the field study of the Kullu valley in July 2018, and 149 landslide locations were collected as global positioning system (GPS) points. The present study evaluates the LSM using three different spatial resolution of the digital elevation model (DEM) derived from three different sources. The data-driven traditional frequency ratio (FR) model was used for this study. The FR model was used for this research to assess the impact of the different spatial resolution of DEMs on the LSM. DEM data was derived from Advanced Land Observing Satellite-1 (ALOS) Phased Array type L-band Synthetic Aperture Radar (PALSAR) ALOS-PALSAR for 12.5 m, the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Global for 30 m, and the Shuttle Radar Topography Mission (SRTM) for 90 m. As an input, we used eight landslide conditioning factors based on the study area and topographic features of the Kullu valley in the Himalayas. The ASTER-Global 30m DEM showed higher accuracy of 0.910 compared to 0.839 for 12.5 m and 0.824 for 90 m DEM resolution. This study shows that that 30 m resolution is better suited for LSM for the Kullu valley region in the Himalayas. The LSM can be used for mitigation and future planning for spatial planners and developmental authorities in the region.

Estimating landslides vulnerability in Rwanda using analytic hierarchy process and geographic information system

Landslides are among hazards that undermine the social, economic, and environmental well-being of the vulnerable community. Assessment of landslides vulnerability reveals the damages that could be recorded, estimates the severity of the impact, and increases the preparedness, response, recovery, and mitigation as well. This study aims to estimate landslides vulnerability for the western province of Rwanda. Field survey and secondary data sources identified 96 landslides used to prepare a landslides inventory map. Ten factors-altitude, slope angles, normalized difference vegetation index (NVDI), land use, distance to roads, soil texture, rainfall, lithology, population density, and possession rate of communication tools-were analyzed. The Analytical Hierarchy Process (AHP) model was used to weight and rank the vulnerability conditioning factors. Then the Weighted Linear Combination (WLC) in geographic information system (GIS) spatially estimated landslides vulnerability over the study area. The results indicated the altitude (19.7%), slope angles (16.1%), soil texture (14.3%), lithology (13.5%), and rainfall (12.2%) as the major vulnerability conditioning parameters. The produced landslides vulnerability map is divided into 5 classes: very low, low, moderate, high and very high. The proposed method is validated by using the relative landslides density index (R-index) method, which revealed that 35.4%, 25%, and 23.9% of past landslides are observed within moderate, high, and very high vulnerability zones, respectively. The consistency of validation indicates good performance of the methodology used and the vulnerability map prepared. The results can be used by policy makers to recognize hazard vulnerability lessening and future planning needs. Integr Environ Assess Manag 2019;00:000-000. © 2019 SETAC.

Comparing probabilistic and statistical methods in landslide susceptibility modeling in Rwanda/Centre-Eastern Africa

Application of suitable methods to generate landslide susceptibility maps (LSM) can play a key role in risk management. Rwanda, located in centre-eastern Africa experiences frequent and intense landslides which cause substantial impacts. The main aim of the current study was to effectively generate susceptibility maps through exploring and comparing different statistical and probabilistic models. These included weights of evidence (WoE), logistic regression (LR), frequency ratio (FR) and statistical index (SI). Experiments were conducted in Rwanda as a study area. Past landslide locations have been identified through extensive field surveys and historical records. Totally, 692 landslide points were collected and prepared to produce the inventory map. This was applied to calibrate and validate the models. Fourteen maps of conditioning factors were produced for landslide susceptibility modeling, namely: elevation, slope degree, topographic wetness index (TWI), curvature, aspect, distance from rivers and streams, distance to main roads, lithology, soil texture, soil depth, topographic factor (LS), land use/land cover (LULC), precipitation and normalized difference vegetation index (NDVI). Thus, the produced susceptibility maps were validated using the receiver operating characteristic curves (ROC/AUC). The findings from this study disclosed that prediction rates were 92.7%, 86.9%, 81.2% and 79.5% respectively for WoE, FR, LR and SI models. The WoE achieved the highest AUC value (92.7%) while the SI produced a lowest AUC value (79.5%). Additionally, 20.42% of Rwanda (5048.07 km²) was modeled as highly susceptible to landslides with the western part the highly susceptible comparing to other parts of the country. Conclusively, the comparison of produced maps revealed that all applied models are promising approaches for landslide susceptibility studying in Rwanda. The results of the present study may be useful for landslide risk mitigation in the study area and in other areas with similar terrain and geomorphological conditions. More studies should be performed to include other important conditioning factors that exacerbate increases in susceptibility especially anthropogenic factors.

A Comparative Study of Statistics-Based Landslide Susceptibility Models: A Case Study of the Region Affected by the Gorkha Earthquake in Nepal

As a result of the Gorkha earthquake in 2015, about 9000 people lost their lives and many more were injured. Most of these losses were caused by earthquake-induced landslides. Sustainable planning and decision-making are required to reduce the losses caused by earthquakes and related hazards. The use of remote sensing and geographic information systems (GIS) for landslide susceptibility mapping can help planning authorities to prepare for and mitigate the consequences of future hazards. In this study, we developed landslide susceptibility maps using GIS-based statistical models at the regional level in central Nepal. Our study area included the districts affected by landslides after the Gorkha earthquake and its aftershocks. We used the 23,439 landslide locations obtained from high-resolution satellite imagery to evaluate the differences in landslide susceptibility using analytical hierarchy process (AHP), frequency ratio (FR) and hybrid spatial multi-criteria evaluation (SMCE) models. The nine landslide conditioning factors of lithology, land cover, precipitation, slope, aspect, elevation, distance to roads, distance to drainage and distance to faults were used as the input data for the applied landslide susceptibility mapping (LSM) models. The spatial correlation of landslides and these factors were identified using GIS-based statistical models. We divided the inventory into data used for training the statistical models (70%) and data used for validation (30%). Receiver operating characteristics (ROC) and the relative landslide density index (R-index) were used to validate the results. The area under the curve (AUC) values obtained from the ROC approach for AHP, FR and hybrid SMCE were 0.902, 0.905 and 0.91, respectively. The index of relative landslide density, R-index, values in sample datasets of AHP, FR and hybrid SMCE maps were 53%, 58% and 59% for the very high hazard classes. The final susceptibility results will be beneficial for regional planning and sustainable hazard mitigation.