Next Article in Journal
Investigating Roundabout Properties and Bicycle Accident Occurrence at Swiss Roundabouts: A Logistic Regression Approach
Next Article in Special Issue
Applicability of Remote Sensing-Based Vegetation Water Content in Modeling Lightning-Caused Forest Fire Occurrences
Previous Article in Journal
Delimitating Urban Commercial Central Districts by Combining Kernel Density Estimation and Road Intersections: A Case Study in Nanjing City, China
Previous Article in Special Issue
Scenario-Based Risk Assessment of Earthquake Disaster Using Slope Displacement, PGA, and Population Density in the Guyuan Region, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

Department of Geoinformatics—Z_GIS, University of Salzburg, 5020 Salzburg, Austria
*
Author to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2019, 8(2), 94; https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi8020094
Submission received: 23 December 2018 / Revised: 2 February 2019 / Accepted: 12 February 2019 / Published: 18 February 2019
(This article belongs to the Special Issue Natural Hazards and Geospatial Information)

Abstract

:
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.

Graphical Abstract

1. Introduction

Mass movements are one of the major geological hazards in Nepal; almost 80% of the total area of the country is prone to landslides [1]. Nepal is in the central part of the Himalayan mountain range, located in the intercontinental collision zone where the Indian plate is subducting under the Eurasian plate [2]. On 25 April 2015, an earthquake of magnitude (M) 7.8 struck Nepal and was followed by several aftershocks. The earthquake shook most of the country and neighboring countries, and it was one of the most powerful earthquakes since the Bihar–Nepal earthquake of 1934. The main event and aftershocks triggered mass movements and resulted in damage to infrastructure and human life. The Gorkha earthquake resulted in the deaths of more than 9000 people and is estimated to affect about eighty-five million people in Nepal and neighboring countries [3]. The estimates of previous studies as to the number of landslides triggered by the Gorkha earthquake range between 3000 and 25,000 [4,5,6,7,8,9]. A landslide inventory is one of the essential aspects for landslide susceptibility evaluation [10,11]. In most of the studies related to the Gorkha earthquake, the emphasis was on the development of a landslide inventory while using this information to develop a landslide susceptibility zonation for the earthquake-affected region was largely neglected. There has been extensive work done on the preparation of coseismic landslide inventories for susceptibility evaluations worldwide. The first coseismic landslide inventory was prepared after the 1957 San Francisco earthquake [12,13]. Since then, several researchers have attempted to prepare coseismic landslide inventories over the years. Due to advancements in geographic information systems (GIS) and Earth observation (EO) technologies, including semi-automatic and machine learning techniques, the preparation of a coseismic landslide inventory is becoming more attainable and less time-consuming.
Landslide susceptibility mapping (LSM) is done by assessing the probability of a landslide occurrence in a given region [14]. Over the years, landslide susceptibility mapping has become a practical approach with the aim to obtain better insights into the potential slope failures. Remote sensing and GIS techniques are the foundation of landslide susceptibility mapping. In this study, we used landslide inventory data obtained from very high-resolution Digital Globe Worldview 2–3 imagery. Satellite imagery of high spatial resolution helps in identifying even the smallest of the slope failures. Many previous studies have used very high-resolution satellite imagery to obtain landslide inventory data, for instance [15] used 2 m Worldview-2 imagery to detect landslides caused by the 2010 Haiti earthquake. Another study used IKONOS-1 imagery for landslide detection in the Tangjiaoa, Wanzhou district of China [16]. Reference [17] used different CNN models and three different machine learning models, namely artificial neural network (ANN), support vector machines (SVM), and random forest (RF) for landslide detection and inventory dataset production. High-resolution satellite imagery data is useful to generate input datasets for LSM. Integration of different models for the identification of probable areas of landslide occurrence through LSM is a recognized approach for sustainable risk mitigation [18]. Access to accurate inventory data along with hazard conditioning factors using both data-driven and knowledge-based models have yielded successful results for LSM in different studies [19,20,21]. Several studies have used expert knowledge-based models for conducting LSM. A range of models is used in these studies, including the analytical hierarchy process (AHP) [22,23,24,25]. Although the AHP yields accurate results in natural hazard susceptibility modelling and mapping, it has been proved that there is a degree of uncertainty associated with the resulting maps. Therefore, the AHP model is not always a satisfactory model for susceptibility mapping based on decision makers’ preferences. A number of attempts have been made to improve the accuracy of the AHP as it is potentially a valuable tool for susceptibility mapping [19,26].
Currently, data-based models show fairly good results in the natural hazard susceptibility assessments [27,28,29]. We used a data-based frequency ratio (FR) model for LSM in the Gorkha region. FR is based on data from previous landslides and causal factors and their spatial coverage [30]. Several studies used FR to generate LSMs at regional or small scales, such as [24,30,31,32,33,34,35,36]. A study carried out in Selangor, Malaysia by [37] concluded that the FR model shows a correlation between landslide causal factors and the location of past landslides in a given study area. As the weightings of the FR model are directly based on the density of the past landslides in the study area, it is sensitive to the use of different training datasets.
Several researchers have integrated different models to improve the susceptibility assessment of landslides. These studies used a range of models. For example, the entropy method, logistic regression (LR) and support vector machine (SVM) were combined for LSM by [34], and the ensemble weights of evidence, LR and random forest models were integrated into a hybrid approach for LSM in Shaanxi Province, China by [38]. The integration of the AHP model with a sensitivity and uncertainty analysis led to better results of LSM in a study that have done with [39,40].
In this study, the integration of the AHP and FR models was used as the theoretical background of the spatial multi-criteria evaluation (SMCE) model. The multi-criteria evaluation (MCE) is a decision-based approach in which alternatives are evaluated in a tree-like hierarchy of criteria and objectives [41]. There are several stages in carrying out a hybrid SMCE. Firstly, a problem analysis is carried out by making a tree-like hierarchy; this is followed by the standardization of the weights of the factors and the generation of the final composite maps [41,42,43,44]. The hybrid SMCE model yielded a higher prediction accuracy, resulting in higher reliability of the generated landslide maps, compared to FR, SVM, index of entropy (IOE) and the weighted linear combination (WLC) models implemented in a study done in the Wunning area in China by [45].
We used the AHP, FR, and hybrid SMCE models to carry out LSM for the study area. First, the AHP model was created based on expert knowledge, which was used to assign weights to factors. Then, the FR model was used to generate landslide susceptibility maps based on a data-driven approach. LSM was carried out based on different class frequency ratios of landslide inventory data. Finally, the third landslide susceptibility map was generated using the hybrid SMCE model, which is the integration of the AHP and FR models. In this study, we selected the region most affected by landslides after the earthquake; the districts are Gorkha, Rasuwa, Sindhupalchok, Dhading, Nuwakot, and Dholaka. GIS-based landslide susceptibility maps were produced using the three statistical models AHP, FR and hybrid SMCE. The validation of the three models was carried out using two different accuracy assessment methods and the validation inventory dataset of the past landslide occurrences in our case study area.

2. Study Area and Inventory Data Set

The epicenter of the earthquake was in the Gorkha district and its aftershock was about 140 km away in the Dholaka district. But landslides triggered by the main event and several aftershocks were scattered all over central Nepal. We have selected seven districts that were severely affected by the earthquake. The study area covers a 14,502 km² region spanning most of central Nepal. The area has two major drainage systems, namely Narayani and Saptakoshi [6]. A total of 23,439 landslides have been selected for the landslide susceptibility evaluation (see Figure 1).
Geologically, the study area lies within the Himalayan fold and thrust zone of central Nepal, which resulted from the collision of the Indian plate with the Eurasian plate [2]. The collision led to extensive crustal shortening and upheaval, resulting in the formation of the Himalaya, the quintessential collision orogen. Himalayan seismicity is attributed to the movement of the Indian plate relative to the Eurasian plate at a rate of about 5 cm per year. When the convergence is locked in some areas, the energy is ultimately released in the form of tremors along the mountain range and in its surroundings [46].

Landslide Inventory Data Set

The landslide inventory map illustrates the active landslide sites and their properties, such as landslide type, structural attributes, and distance to roads. These slope deformations are linked to the morphological, geological and climatic conditions of the locations; thus, these attributes can predict the conditions that could cause future landslides in the area. The first step was to identify the landslide locations in the satellite imagery and evaluate landslide-prone areas. Then active landslide locations were mapped and an inventory was prepared using satellite image interpretation, an extensive field survey, and literature search for historical landslide records [47]. The landslide inventory dataset was generated by comparing high-resolution imagery of pre- and post-event (Figure 2) data taken from the USGS data catalogue [48]. The satellite imagery was obtained from Digital Globe Worldview 2 and 3, and the spatial resolution of the imagery used for mapping the landslides was very high in the study area. The images were collected between 26th April and 15th June 2015, thus covering the main event of the Gorkha earthquake and its aftershocks. For our study area, a total of 23,439 landslide locations were identified.
Our landslide inventory dataset was separated into a training and a testing dataset. Dividing the inventory dataset randomly for training and validation purposes is a general approach that has been used in several natural hazard studies [49,50]. Samples are chosen based on the size of the study area, completeness of the inventory, and the methodology used. Currently, there is no standard methodology for the selection of testing and training samples [51]. Reference [52] have given different ratios for various methods, and [53] also used the same ratio for prediction of landslides using computational methods. To make our models computationally robust, the inventory dataset was randomly divided into two sections, with 70% (16,407) used for training the models and 30% (7032) used for validating the results (see Figure 2).

3. Materials and Methods

3.1. Landslide Conditioning Factors

The nine factors that most affect landslide occurrence in our study area are categorized into the four main groups of topographical, geological, hydrological and anthropogenic factors. Topographical factors include slope, elevation and aspect. The SRTM digital elevation model (DEM) with a 30 m spatial resolution was used as the basis of our slope and aspect layers. The slope angle is considered as the most critical factor of the slope stability assessment in any LSM study [54]. The elevation is another crucial factor of LSM as most geomorphological and geological processes are controlled by this topographical factor [55]. It was also used to describe the local relief of the study area. Relief classes denote the elevation range from the lowest to the highest point in the study area [56]. Five altitude classes were defined to determine the density of landslides in different relief classes. The aspect factor affects the landslide occurrence as it is related to meteorological factors such as the direction of the precipitation and the average amount of sunshine.
The geological factors in this study are defined as lithology and distance to faults. These geological factors were extracted from the available 1:250,000 geological map of central Nepal. Different lithologic units and proximity to faults play a significant role in controlling the landslides. More than 40% of the lithology of the study area is categorized as gneiss and migmatite. The distance to faults is essential, especially in our case study, since areas closer to the faults were more affected by the earthquake [6].
Precipitation and the distance to drainage were our hydrological factors. Drainage is also one of the causal factors of the landslides because they cause erosion and saturation of the materials in the lower parts of valleys [57,58].
The anthropogenic factors are the distance to roads and land use. Human activities, and road construction in particular, cause a diminution in the load on either the toe or the topography of slopes. Finally, to apply the three models, the layers of all conditioning factors were included in the GIS environment as raster layers with a 30 m resolution (see Figure 3). The scale of 30 m was selected to match the spatial resolution of the SRTM 30 m data on which the topography data set factors were based.

3.2. Methods

Landslide susceptibility analysis was carried out using the AHP, FR and hybrid SMCE models in the region affected by the Gorkha earthquake.

3.2.1. Analytical Hierarchy Process (AHP)

The AHP developed by [59] has been used to weight related factors of any spatial problem in GIS environments [60]. It is a common tool for analysis to support site selection, urban planning, and natural hazard susceptibility analysis [61]. The AHP is a decision-making process based on multi-criteria and multiple objectives and involves the participation of experts [62]. A hierarchical order of factors and numerical values is established based on the importance of each factor [63]. Subsequently, these factors are integrated, and each factor is weighted according to its importance [64]. The AHP is used to establish the correlative pairwise comparison matrix. This matrix is constructed using values that represent experts’ judgments by comparing the importance of each factor in relation to all the other related factors [65]. Each layer is based on a nine-point rating scale and added to the matrix as developed by [66] (See Table 1). The decision maker specifies the value of each factor. Each factor weight from the matrix class was multiplied by the weight class. Local representation of factors in the study area determined the results of the susceptibility map. The factors were weighted using the pairwise comparison matrices of the AHP based on expert knowledge. The principle of transitivity is important in AHP for any given three factors (such as f1, f2 and f3), and is defined as follows; if f1 > f2 and f2 > f3, then f1 > f3. The principle of transitivity is a basis for conditioning factors weighing in AHP. Due to this principle, a consistent pairwise comparison matrix would require that if 2f1 > f2 (i.e., f1 is two times more preferable than f2) and 4f2 > f3, then 8f1 > f3 to account for the transitivity principle [19,67]. Therefore, it was necessary to compute the consistency of expert comparisons in matrices in each stage [20]. Inconsistency can be defined based on the observation that λmax > n for comparison matrices and λmax = n, if C is a consistent comparison. The consistency ratio (CR) is defined by Equation (1):
CR = ( λ max n ) / ( RI ( n 1 ) ) ,
where RI is the random index of a randomly created pairwise comparison matrix for n = 2, 3, 4, 5, 6, 7, 8, and 9. If the consistency ratio is <0.10 it indicates an acceptable level of consistency, whereas a CR > 0.10 points to a degree of inconsistency [68]. Our case study in this research is to map areas susceptible to landslides. Accordingly, to calculate the criteria weights using the AHP model, questionnaire forms were prepared to gather expert knowledge for layer weighting. Therefore, nine local academic experts and twelve researchers and geologists from the international center for integrated mountain development (ICIMOD) and national organizations involved in post-earthquake reconstruction planning were asked to fill in the questionnaire, and a total of twenty-one forms were obtained for layer weighting. The assessed weights of the nine layers were calculated using the AHP model based on the resulting pairwise comparison matrices.

3.2.2. Frequency Ratio (FR)

FR is a useful geospatial assessment tool that provides the probabilities of the distribution of the occurrence and non-occurrence of landslides for each class of related factors [36]. Classes are weighted based on the ratio of observed landslides to the whole study area [58]. FR is one of the best geospatial assessment tools to determine the spatial correlation between inventory data and the class of related factors [69]. The number or percentage of the inventory data in each class indicates the significance of the correlation with the landslide [70]. For computing the FR weights, the frequency ratios of the factor classes were calculated based on our landslide inventory dataset. The percentage ratio for different classes of each factor to the whole study area was calculated. The final LSM was generated by a linear combination model (see Equation (2)) of the FR weights for all the factors.
L S M F R = i = 1 n J = 1 k F R i j [ f a c t o r i ( c l a s s j ) ] ,
where n is the number of factors and k is the number of classes in each factor that contributed to the landslide susceptibility map generation. Thus, FRij is the FR weight of the j-th class in the i-th factor. To prepare qualitative landslide susceptibility maps, the FR model was implemented using GIS tools. The frequency ratio of each factor was calculated from the relationship of the area of the corresponding class with the landslide events. The FR was calculated for all classes of each conditioning factor. According to the FR model, the ratio is defined as the area where the landslide occurred to the total area. This means that an FR value of 1 is an average value for the occurrence of a landslide in a particular area compared to the whole area. If the value is >1, it means the percentage of landslides that occurred in the area is above the regional average and shows a high correlation, whereas values lower than 1 indicate a lower correlation [71,72]. The generated landslide susceptibility map based on the FR weights is shown in Figure 4.

3.2.3. Hybrid Spatial Multi-Criteria Evaluation (SMCE)

The hybrid SMCE model enables users to solve the spatial problems of natural hazard susceptibility mapping [70]. Spatial features are defined as lines, points, and areas in this approach and the generated final maps result from the consideration of landslide conditioning factors [73]. SMCE is an approach that makes it possible to incorporate spatial analysis and GIS to use both spatial and non-spatial input data to produce final maps [74]. In hybrid SMCE, input layers are spatially represented as factors. Based on the criteria tree, input layers are grouped weighted and standardized. The input layers are standardized from their original values to the 0–1 value range.
The outputs of the hybrid SMCE were the index maps, which represent the match and mismatch of criteria in different areas. The multi-criteria evaluation (MCE) of the AHP was used as the theory of the hybrid SMCE model. The steps involved in the operation of the hybrid SMCE were the problem analysis, weighting the factors, standardization, ranking the factors using priority ordering, and finally producing the output map [41]. The conditioning factors were grouped into the four main groups of topographical, geological, hydrological, and anthropogenic. The significance of each group on the landslide susceptibility was calculated based on the AHP pairwise comparison matrix completed by the same experts (see Table 2). The same methodology was used for weighting the factors categorized into different hybrid SMCE groups. However, as there were some differences in the expert’s reports regarding the pairwise comparison of the hybrid SMCE groups, we used the weighted average of reported values for this matrix. Using the average of the comparison values of Saaty’s rating scale (see Table 1) from different expert knowledge is a common approach for dealing with the uncertainty sources in the multi-criteria decision analyses [60,63]. The hybrid SMCE resulting weight for each conditioning factor is represented in Table 3. These weights are the same with those of used for the AHP model resulting from the pairwise comparison matrix. However, in case of the hybrid SMCE, we used the resulting weights from the FR model for the classes of each conditioning factor. Therefore, the weightings of the factors were produced based on the AHP pairwise comparison matrices, whereas the FR values were used for the classes. Hence, the final impact of each hybrid SMCE factor on the resulting landslide susceptibility map was an integration of the weighting processes based on both the AHP and FR models, and it was considered the result of both models.

4. Results

To generate the landslide susceptibility maps to identify the areas that are highly susceptible to landslides, three different models were used in this study. The output values of the AHP and FR models were represented in Table 4 and Table 5, respectively. The resulting weights from each model were applied to weight the conditioning factors and classes in a GIS environment. The derived FR value of over unity shows a strong relationship between the landslide data of the training inventory dataset high landslide susceptibility. FR values of less than unity show a low susceptibility. This is also true for the resulting factor weights from the AHP. Therefore, higher AHP weights indicated a high susceptibility. The weightings of the FR model were derived based on our produced inventory data set. However, those of the AHP model were generated according to the expert preferences in the pairwise comparison matrices. The weightings for sub-criteria classes from both the FR and AHP models followed the same pattern, there were some big differences in the first class of land use, distance to faults, distance to roads, and slope factors. According to the results based on expert knowledge, the factor of lithology containing the quartzite rock formation had the most significant impact on landslide occurrence, and this class got the lowest FR value. However, the rest of the classes were approximately the same in terms of the FR value. This is also the case for the weightings of the classes of the land use factor.
For the hybrid SMCE approach, the main groups and their corresponding factors were weighted based on the AHP model, and the FR values were considered for the classes of the factors. Thus, the resulting landslide susceptibility map based on the hybrid SMCE approach was a combination of the results of both the AHP and the FR models. Figure 4 shows the resulting landslide susceptibility maps obtained from each model. The natural breaks classification technique was used to classify similar values separated by breakpoints. This is a common and useful technique for classification of any hazard susceptibility map. The results were interpreted as belonging to the same class when they were values close to each class boundary, e.g., values between “low” and “moderate” susceptibility. Thus, the resulting landslide susceptibility maps were classified into five classes of susceptibility using the natural breaks classification technique. These results were validated using the 30% validation inventory dataset. We explain the validation approach and results in the next section.

5. Validation

In this section, we outline two different validation methods, which are widely used in the accuracy assessment of natural hazard susceptibility mapping. Validation is an important part of preparing natural hazard susceptibility maps to predict possible future events [76]. To ascertain the success of applied landslide susceptibility models, we compared the resulting maps with the landslide inventory data from the study area. Analyzing the conformity between the inventory data and the resulting maps of the applied models can give a clear indication as to the effectiveness of each model for LSM and indicate whether the applied models could correctly predict areas that are susceptible to landslides. The locations of the historical landslides support the modelling and pattern recognition. Addressing more accurate areas that were exposed and susceptible to landslides becomes an advantageous driver to inform spatial managers in order to inform sustainable risk management and landslide mitigation. It would be useful for minimizing the adverse impacts, mainly by saving lives of the local people who live in such susceptible areas and protecting public properties in the area under investigation [44]. In this study, 30% (7032 points) of the landslide points of the inventory data set was reserved for the validation. The receiver operating characteristics (ROC) curves and the relative landslide density method (R-index) were used for the accuracy assessment of the resulting landslide susceptibility maps based on the AHP, FR and hybrid SMCE models. After preparing susceptibility maps using the three models, we compared the results with the ground truth at a randomly selected sub-area as shown in Figure 5. The presented graphical verification indicates that the hybrid SMCE model gave a more accurate representation of the landslides that fell into the very high susceptibility class. Validation of the results showed good agreement between the observed and the predicted values for the hybrid SMCE model.

5.1. Receiver Operating Characteristics (ROC)

The ROC curve [77] was used to validate different resulting landslide susceptibility maps using our validation data. The ROC approach displays a full scene of trade-off among the true-positive rate (TPR) and the false-positive rate (FPR) in the landslide susceptibility maps [78]. ROC curves were calculated for all three resulting landslide susceptibility maps. The vertical axis indicates the TPR, while the horizontal axis shows the FPR [79] (see Figure 6). TPRs are the pixels that correctly referred to the landslide areas, while FPRs are the pixels wrongly labelled as landslides. To generating ROC curves, we obtain the correctly and incorrectly labelled pixels in the inventory dataset and then plot TPRs versus FPRs across the values that they took. The area under the curve (AUC) is the measure that indicates the accuracy of the landslide susceptibility maps. The resulting AUCs indicate the probability that more pixels were correctly labelled than incorrectly labelled [27]. Therefore, greater AUC values indicate a higher accuracy of the resulting susceptibility map. The AUC values close to unity indicate a comprehensive susceptibility map; a value of 0.5 shows a worthless map because it means the map was generated by chance [80]. The AUC values obtained from the ROC approach for AHP, FR and hybrid SMCE were 0.902, 0.905 and 0.910, respectively. Therefore, according to the results, the hybrid SMCE model indicates a more accurate landslide susceptibility map compared to the other two methods. Moreover, the resulting landslide susceptibility map based on the FR model was more accurate than that based on the AHP, which means that the data-based model performed better than the knowledge-based model. However, the result of the knowledge-based models may vary if we use other experts for criteria weightings or other models (e.g., technique for order of preference by similarity to ideal solution (TOPSIS) and conventional AHP). The resulting maps based on the data-based models are also sensitive to the training data, and they will change even when using another section of the inventory data set. The resulting landslide susceptibility map based on the FR contained 6701 landslides in its very high susceptibility class, while those of the AHP and hybrid SMCE contained 7948 and 4454, respectively.

5.2. Relative Landslide Density (R-index)

The accuracy of the landslide susceptibility prediction was evaluated using an index of relative landslide density (R-index). We used a total of 23,439 landslides to validate the results. The R-index is given as follows, by [81] Equation (3):
R = (ni/Ni)/Σ(ni/Ni)) ×100
where ni is the percentage of the area that is susceptible to landslides in each susceptibility class, and Ni is the percentage of landslides in each susceptibility class. Results show that the hybrid SMCE and FR models have higher R-index values in the very high susceptibility classes. For low and moderate susceptibility classes, the R-index values are very similar for all three models. The R-index values for the very high susceptibility class in the hybrid SMCE, AHP and FR maps are 59%, 53%, and 58%, respectively (Table 6).

6. Discussion

The area covered by the very high susceptibility class varies from 648,798,300 m² in the hybrid SMCE model to almost three times that area in the AHP model (1,894,469,000 m²), and in the FR model (1,090,157,000 m²) it is more than the area covered by the hybrid SMCE model. In the resulting AHP map, 13.17% of the area falls into the very high susceptibility class, and the number of landslides located in this class was 7948. In the resulting hybrid SMCE map, the corresponding values were 4.51% and 4454, respectively, which resulted in the highest R-index of all three models. As shown in Figure 7, the R- index of the FR model was a value between that of the AHP and the hybrid SMCE. In the map resulting from the FR model, 7.55% of the area falls into the very high susceptibility class with a total of 6701 landslides in this area. The distribution of the locations of landslides in the different classes of the resulting landslide susceptibility maps of a sub-area is shown in Figure 5. For the high susceptibility class, the resulting R-index values and area covered by the corresponding class for the implemented models were almost similar. But the number of landslides that fall into that class of the hybrid SMCE model is higher than in the other two models. When considering the differences in areas of the susceptibility classes in the resulting landslide susceptibility maps and the subsequent R-index for each class, the hybrid SMCE model achieved better results overall.
The overall results of the AUC values obtained from the ROC approach for the three applied models showed that the hybrid SMCE generated a more accurate landslide susceptibility map than the other two models. As the hybrid SMCE is an integration of the AHP and FR models, we can say that the integration of knowledge-based models with the data-based ones achieves a more accurate result for our case study. The R-index value for the very high susceptibility category of the resulting susceptibility map based on the FR model is higher than that of the AHP. The higher accuracy of the resulting landslide susceptibility map based on the FR model compared to the AHP model shows that the data-based model performs better than the knowledge-based model. However, the results of knowledge-based models may vary if we use other experts for criteria weightings or use other models. Thus, the limitation of our work is that we used only the AHP model as our main knowledge-based model. Different multi-criteria decision-making models have been used in several spatial problem studies based on expert knowledge, and in some cases, the resulting accuracies of the models were compared to each other. The AHP model was compared to the TOPSIS by [82] within a fuzzy calculus. Both models were found to be suitable for their study. However, the TOPSIS model yielded better results than the AHP. In another study, the AHP model was compared to the interval pairwise comparison matrix (IPCM) by [83] and the second model was slightly more accurate in LSM. Generally, the AHP model was optimized through integration with other methodologies, such as fuzzy theory and the Dempster Shafer, to improve the accuracy of the final resulting landslide susceptibility maps [84]. However, the conventional AHP model was used for LSM in our study. Consequently, as it was not the intent of our study to analyze the data-based and knowledge-based models separately, but rather to determine the value in combining them, the results of our study should not be considered as the best in terms of the individual models. There are also several studies in the field of landslide susceptibility that implement and compare data-based models [85]. The resulting maps based on the data-based models are sensitive to the training data, and they will change by using a different selection of the inventory dataset. On the one hand, the results from the knowledge-based models may vary if using other experts for the criteria weightings. On the other hand, the resulting natural hazard susceptibility maps based on the data-based models can easily be affected by the training data, and the results will change if using training data from another source. Based on the data presented in Figure 8, although most of the weightings of the sub-criteria based on the AHP and FR models are the same, there are also some significant differences in the first class of the precipitation and distance to roads factors. According to the results based on the AHP, the class of 2500–3275 (mm/yr) of precipitation is the most important in terms of the landslide susceptibility. However, based on the weighting of the FR model, the class of 950–1275 (mm/yr) is the most significant. The weightings of the AHP and FR models for the criteria of elevation and distance to drainage follow a very similar pattern. The weightings of both models show that the impact of the drainage on the landslide susceptibility continuously decreases with increasing distance.

7. Conclusions

The production of landslide susceptibility maps is the first step in managing a sustainable risk mitigation program in any landslide-prone area. For the entire area affected by the Gorkha earthquake, no previous attempts have been made to produce susceptibility maps using multiple statistical approaches. In this research, we produced landslide susceptibility maps for the entire Gorkha region. As the purpose of LSM is to predict the probability of future landslide occurrence in an area, our aim was to identify which of the statistical techniques (AHP, FR, or hybrid SMCE) is more successful for LSM. For this purpose, we used a total of 23,439 landslides and randomly divided them into two datasets, 70% for training and 30% for validation. The degree of fit and accuracy of the resulting landslide susceptibility maps was validated using two approaches, namely the ROC and R-index. The results of these validation approaches show that although all models were applicable for landslide susceptibility mapping, the hybrid SMCE model (0.910) gives slightly better results for ROC than the AHP (0.902) and FR (0.905) models. The hybrid SMCE resulted in better prediction accuracy and reliability of output landslide susceptibility maps. In summary, these investigations have resulted in the understanding that the integration of knowledge-based and data-based models outperforms the use of these models individually for the generation of landslide susceptibility maps. The FR model slightly outperforms the knowledge-driven AHP model. Reliable and precise susceptibility maps can minimize the costs and damage from natural disasters, such as landslides, if prepared in advance. The output maps can assist decision makers and planners to identify areas that are susceptible to future hazards. It will lead to the mitigation of future damage to infrastructure and human life.

Author Contributions

Conceptualization, O.G., S.R.M., and T.B.; Data curation, S.R.M.; Funding acquisition, T.B.; Investigation, S.R.M. and O.G..; Methodology, O.G. and S.R.M.; Supervision, T.B.; Validation, S.R.M. and O.G.; Visualization, S.R.M.; Writing—original draft, O.G. and S.R.M.; Writing—review & editing, T.B. All authors read and approved the final manuscript.

Funding

This research is partly funded by the Austrian Science Fund (FWF) through the GIScience Doctoral College (DK W 1237-N23).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Manandhar, P.K. Disaster Review 2013; Department of Water Induced Disaster Prevention: Kathmandu, Nepal, 2015.
  2. Ni, J.; Barazangi, M. Seismotectonics of the himalayan collision zone: Geometry of the underthrusting indian plate beneath the himalaya. J. Geophys. Res. Solid Earth 1984, 89, 1147–1163. [Google Scholar] [CrossRef]
  3. Dahal, R.K. Engineering geological issues after Gorkha earthquake 2015 in Nepal—A preliminary understanding. In Proceedings of the 10th Asian Regional Conference of IAEG, Kyoto, Japan, 26–29 September 2015. [Google Scholar]
  4. Kargel, J.; Leonard, G.; Shugar, D.H.; Haritashya, U.; Bevington, A.; Fielding, E.; Fujita, K.; Geertsema, M.; Miles, E.; Steiner, J. Geomorphic and geologic controls of geohazards induced by Nepal’s 2015 Gorkha earthquake. Science 2016, 351, aac8353. [Google Scholar] [CrossRef] [PubMed]
  5. Martha, T.R.; Roy, P.; Mazumdar, R.; Govindharaj, K.B.; Kumar, K.V. Spatial characteristics of landslides triggered by the 2015 m w 7.8 (Gorkha) and m w 7.3 (Dolakha) earthquakes in Nepal. Landslides 2017, 14, 697–704. [Google Scholar] [CrossRef]
  6. Regmi, A.D.; Dhital, M.R.; Zhang, J.-Q.; Su, L.-J.; Chen, X.-Q. Landslide susceptibility assessment of the region affected by the 25 April 2015 Gorkha earthquake of Nepal. J. Mt. Sci. 2016, 13, 1941–1957. [Google Scholar] [CrossRef]
  7. Roback, K.; Clark, M.K.; West, A.J.; Zekkos, D.; Li, G.; Gallen, S.F.; Chamlagain, D.; Godt, J.W. The size, distribution, and mobility of landslides caused by the 2015 mw 7.8 Gorkha earthquake, Nepal. Geomorphology 2018, 301, 121–138. [Google Scholar] [CrossRef]
  8. Chigira, M. Inventory Mapping of Landslides Induced by the Gorkha Earthquake. Available online: https://www.jst.go.jp/inter/sicp/country/j-rapid/nepal/pdf/j-rapid_nepal_chigira.pdf (accessed on 5 February 2019).
  9. Gnyawali, K.R.; Maka, S.; Adhikari, B.R.; Chamlagain, D.; Duwal, S.; Dhungana, A.R. Spatial implications of earthquake induced landslides triggered by the April 25 Gorkha earthquake mw 7.8: Preliminary analysis and findings. In Proceedings of the International Conference on Earthquake Engineering and Post Disastor Reconstruction Planning, Bhaktapur, Nepal, 24–26 April 2016; pp. 50–58. [Google Scholar]
  10. Guzzetti, F.; Reichenbach, P.; Ardizzone, F.; Cardinali, M.; Galli, M. Estimating the quality of landslide susceptibility models. Geomorphology 2006, 81, 166–184. [Google Scholar] [CrossRef]
  11. Hervás, J.; Van Den Eeckhaut, M.; Legorreta, G.; Trigila, A. Landslide inventory and susceptibility and hazard zoning. Landslide Sci. Pract. 2013, 1, 1–2. [Google Scholar]
  12. Keefer, D.K. Landslides caused by earthquakes. Geol. Soc. Am. Bull. 1984, 95, 406–421. [Google Scholar] [CrossRef]
  13. Keefer, D.K. Investigating landslides caused by earthquakes—A historical review. Surv. Geophys. 2002, 23, 473–510. [Google Scholar] [CrossRef]
  14. Gorsevski, P.V.; Brown, M.K.; Panter, K.; Onasch, C.M.; Simic, A.; Snyder, J. Landslide detection and susceptibility mapping using lidar and an artificial neural network approach: A case study in the Cuyahoga Valley National Park, Ohio. Landslides 2016, 13, 467–484. [Google Scholar] [CrossRef]
  15. Saba, S.B.; Ali, M.; van der Meijde, M.; van der Werff, H. Co-seismic landslides automatic detection on regional scale with sub-pixel analysis of multi temporal high resolution optical images: Application to southwest of Port Au Prince, Haiti. J. Himal. Earth Sci. 2017, 50, 74–92. [Google Scholar]
  16. Huang, F.; Chen, L.; Yin, K.; Huang, J.; Gui, L. Object-oriented change detection and damage assessment using high-resolution remote sensing images, Tangjiao Landslide, Three Gorges Reservoir, China. Environ. Earth Sci. 2018, 77, 183. [Google Scholar] [CrossRef]
  17. Ghorbanzadeh, O.; Blaschke, T.; Gholamnia, K.; Meena, S.R.; Tiede, D.; Aryal, J. Evaluation of different machine learning methods and deep-learning convolutional neural networks for landslide detection. Remote Sens. 2019, 11, 196. [Google Scholar] [CrossRef]
  18. Feizizadeh, B.; Jankowski, P.; Blaschke, T. A gis based spatially-explicit sensitivity and uncertainty analysis approach for multi-criteria decision analysis. Comput. Geosci. 2014, 64, 81–95. [Google Scholar] [CrossRef] [PubMed]
  19. Ghorbanzadeh, O.; Feizizadeh, B.; Blaschke, T. Multi-criteria risk evaluation by integrating an analytical network process approach into GIS-based sensitivity and uncertainty analyses. Geomat. Nat. Hazards Risk 2018, 9, 127–151. [Google Scholar] [CrossRef]
  20. Hasekioğulları, G.D.; Ercanoglu, M. A new approach to use ahp in landslide susceptibility mapping: A case study at Yenice (Karabuk, Nw Turkey). Nat. Hazards 2012, 63, 1157–1179. [Google Scholar] [CrossRef]
  21. Kanwal, S.; Atif, S.; Shafiq, M. GIS based landslide susceptibility mapping of northern areas of Pakistan, a case study of Shigar and Shyok Basins. Geomat. Nat. Hazards Risk 2017, 8, 348–366. [Google Scholar] [CrossRef]
  22. Mondal, S.; Maiti, R. Integrating the analytical hierarchy process (AHP) and the frequency ratio (FR) model in landslide susceptibility mapping of Shiv-Khola Watershed, Darjeeling Himalaya. Int. J. Disaster Risk Sci. 2013, 4, 200–212. [Google Scholar] [CrossRef]
  23. Myronidis, D.; Papageorgiou, C.; Theophanous, S. Landslide susceptibility mapping based on landslide history and analytic hierarchy process (AHP). Nat. Hazards 2016, 81, 245–263. [Google Scholar] [CrossRef]
  24. Park, S.; Choi, C.; Kim, B.; Kim, J. Landslide susceptibility mapping using frequency ratio, analytic hierarchy process, logistic regression, and artificial neural network methods at the Inje Area, Korea. Environ. Earth Sci. 2013, 68, 1443–1464. [Google Scholar] [CrossRef]
  25. Pourghasemi, H.R.; Pradhan, B.; Gokceoglu, C. Application of fuzzy logic and analytical hierarchy process (AHP) to landslide susceptibility mapping at Haraz Watershed, iran. Nat. Hazards 2012, 63, 965–996. [Google Scholar] [CrossRef]
  26. Malczewski, J.; Rinner, C. Multicriteria Decision Analysis in Geographic Information Science; Springer: New York, NY, USA, 2016. [Google Scholar]
  27. Ghorbanzadeh, O.; Blaschke, T.; Aryal, J.; Gholaminia, K. A new gis-based technique using an adaptive neuro-fuzzy inference system for land subsidence susceptibility mapping. J. Spat. Sci. 2018. [Google Scholar] [CrossRef]
  28. Ghorbanzadeh, O.; Rostamzadeh, H.; Blaschke, T.; Gholaminia, K.; Aryal, J. A new gis-based data mining technique using an adaptive neuro-fuzzy inference system (ANFIS) and k-fold cross-validation approach for land subsidence susceptibility mapping. Nat. Hazards 2018, 94, 497–517. [Google Scholar] [CrossRef]
  29. Pourghasemi, H.R.; Kerle, N. Random forests and evidential belief function-based landslide susceptibility assessment in western Mazandaran province, Iran. Environ. Earth Sci. 2016, 75, 185. [Google Scholar] [CrossRef]
  30. Wu, Y.; Li, W.; Wang, Q.; Liu, Q.; Yang, D. Landslide Susceptibility Assessment Using Frequency Ratio, Statistical Index and Certainty Factor Models for the Gangu County, China. Arab. J. Geosci. 2016. [Google Scholar] [CrossRef]
  31. Ozdemir, A.; Altural, T. A comparative study of frequency ratio, weights of evidence and logistic regression methods for landslide susceptibility mapping: Sultan Mountains, SW Turkey. J. Asian Earth Sci. 2013, 64, 180–197. [Google Scholar] [CrossRef]
  32. Pradhan, B. Landslide susceptibility mapping of a catchment area using frequency ratio, fuzzy logic and multivariate logistic regression approaches. J. Indian Soc. Remote Sens. 2010, 38, 301–320. [Google Scholar] [CrossRef]
  33. Razandi, Y.; Pourghasemi, H.R.; Neisani, N.S.; Rahmati, O. Application of analytical hierarchy process, frequency ratio, and certainty factor models for groundwater potential mapping using gis. Earth Sci. Inform. 2015, 8, 867–883. [Google Scholar] [CrossRef]
  34. Umar, Z.; Pradhan, B.; Ahmad, A.; Jebur, M.N.; Tehrany, M.S. Earthquake induced landslide susceptibility mapping using an integrated ensemble frequency ratio and logistic regression models in west sumatera province, indonesia. CATENA 2014, 118, 124–135. [Google Scholar] [CrossRef]
  35. Wang, L.-J.; Guo, M.; Sawada, K.; Lin, J.; Zhang, J. A comparative study of landslide susceptibility maps using logistic regression, frequency ratio, decision tree, weights of evidence and artificial neural network. Geosci. J. 2016, 20, 117–136. [Google Scholar] [CrossRef]
  36. Wang, Q.; Li, W. A gis-based comparative evaluation of analytical hierarchy process and frequency ratio models for landslide susceptibility mapping. Phys. Geogr. 2017, 38, 318–337. [Google Scholar] [CrossRef]
  37. Lee, S.; Pradhan, B. Landslide hazard mapping at selangor, malaysia using frequency ratio and logistic regression models. Landslides 2007, 4, 33–41. [Google Scholar] [CrossRef]
  38. Chen, W.; Sun, Z.; Han, J. Landslide susceptibility modeling using integrated ensemble weights of evidence with logistic regression and random forest models. Appl. Sci. 2019, 9, 171. [Google Scholar] [CrossRef]
  39. Feizizadeh, B.; Blaschke, T. An uncertainty and sensitivity analysis approach for gis-based multicriteria landslide susceptibility mapping. Int. J. Geogr. Inf. Sci. 2014, 28, 610–638. [Google Scholar] [CrossRef] [PubMed]
  40. Omid Ghorbanzadeh, B.F.; Blaschke, T.; Khosravi, R. Spatially explicit sensitivity and uncertainty analysis for the landslide risk assessment of the gas pipeline networks. In Proceedings of the 21st AGILE Conference on Geo-Information Science, Lund, Sweden, 12–15 June 2018. [Google Scholar]
  41. Boerboom, L.; Flacke, J.; Sharifi, A.; Alan, O. Web-Based Spatial Multi-Criteria Evaluation (SMCE) Software, ITC Working Paper 1, for the ForestClim Project 25. Available online: https://www.researchgate.net/publication/267999302_Web-based_spatial_multi-criteria_evaluation_SMCE_software (accessed on 5 February 2019).
  42. Alcaraz Tarragüel, A.; Krol, B.; van Westen, C. Analysing the possible impact of landslides and avalanches on cultural heritage in Upper Svaneti, Georgia. J. Cult. Herit. 2012, 13, 453–461. [Google Scholar] [CrossRef]
  43. Pourghasemi, H.R. Landslide susceptibility mapping using a spatial multi criteria evaluation model at Haraz Watershed, Iran. In Terrigenous Mass Movements; Springer: Berlin/Heidelberg, Germany, 2012; pp. 1–398. [Google Scholar]
  44. Nsengiyumva, J.B.; Luo, G.; Nahayo, L.; Huang, X.; Cai, P. Landslide susceptibility assessment using spatial multi-criteria evaluation model in rwanda. Int. J. Environ. Res. Public Health 2018, 15, 243. [Google Scholar] [CrossRef] [PubMed]
  45. Hong, H.; Shahabi, H.; Shirzadi, A.; Chen, W.; Chapi, K.; Ahmad, B.B.; Roodposhti, M.S.; Hesar, A.Y.; Tian, Y.; Bui, D.T. Landslide susceptibility assessment at the Wuning Area, China: A comparison between multi-criteria decision making, bivariate statistical and machine learning methods. Nat. Hazards 2018, 1–40. [Google Scholar] [CrossRef]
  46. Bilham, R.; Larson, K.; Freymueller, J. Gps measurements of present-day convergence across the nepal himalaya. Nature 1997, 386, 61. [Google Scholar] [CrossRef]
  47. Meena, S.R. Web based landslide management system for nepal. In Proceedings of the 33rd Himalaya-Karakorum-Tibet Workshop (HKT), Lausanne, Switzerland, 10–12 September 2018; pp. 109–110. [Google Scholar]
  48. Roback, K.; Clark, M.; West, A.; Zekkos, D.; Li, G.; Gallen, S.; Godt, J. Map data of landslides triggered by the 25 April 2015 mw 7.8 Gorkha, Nepal earthquake. US Geol. Surv. Data Release 2017. [Google Scholar] [CrossRef]
  49. Termeh, S.V.R.; Kornejady, A.; Pourghasemi, H.R.; Keesstra, S. Flood susceptibility mapping using novel ensembles of adaptive neuro fuzzy inference system and metaheuristic algorithms. Sci. Total Environ. 2018, 615, 438–451. [Google Scholar] [CrossRef]
  50. Pourghasemi, H.R. Gis-based forest fire susceptibility mapping in Iran: A comparison between evidential belief function and binary logistic regression models. Scand. J. For. Res. 2016, 31, 80–98. [Google Scholar] [CrossRef]
  51. Guzzetti, F.; Mondini, A.C.; Cardinali, M.; Fiorucci, F.; Santangelo, M.; Chang, K.-T. Landslide inventory maps: New tools for an old problem. Earth-Sci. Rev. 2012, 112, 42–66. [Google Scholar] [CrossRef] [Green Version]
  52. Basheer, I.A.; Hajmeer, M. Artificial neural networks: Fundamentals, computing, design, and application. J. Microbiol. Methods 2000, 43, 3–31. [Google Scholar] [CrossRef]
  53. Bui, D.T.; Pradhan, B.; Lofman, O.; Revhaug, I.; Dick, O.B. Landslide susceptibility mapping at hoa binh province (vietnam) using an adaptive neuro-fuzzy inference system and GIS. Comput. Geosci. 2012, 45, 199–211. [Google Scholar]
  54. Feizizadeh, B.; Roodposhti, M.S.; Jankowski, P.; Blaschke, T. A gis-based extended fuzzy multi-criteria evaluation for landslide susceptibility mapping. Comput. Geosci. 2014, 73, 208–221. [Google Scholar] [CrossRef] [PubMed]
  55. Mohammadi, M. Mass Movement Hazard Analysis and Presentation of Suitable Regional Model Using GIS (Case Study: A Part of Haraz Watershed). Master’s Thesis, Tarbiat Modarres University International Campus, Tehran, Iran, 2008; p. 80. [Google Scholar]
  56. Raja, N.B.; Įiįek, I.; Türkoğlu, N.; Aydin, O.; Kawasaki, A. Correction to: Landslide susceptibility mapping of the sera river basin using logistic regression model. Nat. Hazards 2018, 91, 1423. [Google Scholar] [CrossRef]
  57. Gökceoglu, C.; Aksoy, H. Landslide susceptibility mapping of the slopes in the residual soils of the Mengen Region (turkey) by deterministic stability analyses and image processing techniques. Eng. Geol. 1996, 44, 147–161. [Google Scholar] [CrossRef]
  58. Pradhan, A.M.S.; Kim, Y.-T. Relative effect method of landslide susceptibility zonation in weathered granite soil: A case study in Deokjeok-ri Creek, South Korea. Nat. Hazards 2014, 72, 1189–1217. [Google Scholar] [CrossRef]
  59. Saaty, T.L. Analytic heirarchy process. Wiley Statsref Stat. Ref. Online 1980. [Google Scholar]
  60. Cabrera-Barona, P.; Ghorbanzadeh, O. Comparing classic and interval analytical hierarchy process methodologies for measuring area-level deprivation to analyze health inequalities. Int. J. Environ. Res. Public Health 2018, 15, 140. [Google Scholar] [CrossRef]
  61. Ayalew, L.; Yamagishi, H.; Ugawa, N. Landslide susceptibility mapping using gis-based weighted linear combination, the case in Tsugawa Area of Agano River, Niigata Prefecture, Japan. Landslides 2004, 1, 73–81. [Google Scholar] [CrossRef]
  62. Saaty, T.L.; Vargas, L.G. Inconsistency and rank preservation. J. Math. Psychol. 1984, 28, 205–214. [Google Scholar] [CrossRef]
  63. Ghorbanzadeh, O.; Moslem, S.; Blaschke, T.; Duleba, S. Sustainable urban transport planning considering different stakeholder groups by an interval-ahp decision support model. Sustainability 2018, 11, 9. [Google Scholar] [CrossRef]
  64. Sahnoun, H.; Serbaji, M.M.; Karray, B.; Medhioub, K. GIS and multi-criteria analysis to select potential sites of agro-industrial complex. Environ. Earth Sci. 2012, 66, 2477–2489. [Google Scholar] [CrossRef]
  65. Pirnazar, M.; Karimi, A.Z.; Feizizadeh, B.; Ostad-Ali-Askari, K.; Eslamian, S.; Hasheminasab, H.; Ghorbanzadeh, O.; Hamedani, M.H. Assessing flood hazard using gis based multi-criteria decision making approach; study area: East-Azerbaijan Province (Kaleybar Chay Basin). J. Flood Eng. 2017, 8, 203–223. [Google Scholar]
  66. Saaty, T.L. A scaling method for priorities in hierarchical structures. J. Math. Psychol. 1977, 15, 234–281. [Google Scholar] [CrossRef]
  67. Malczewski, J.; Rinner, C. Introduction to GIS-mcda. In Multicriteria Decision Analysis in Geographic Information Science; Springer: New York, NY, USA, 2015; pp. 23–54. [Google Scholar]
  68. Saaty, T.L. The Analytic Hierarchy Process: Planning, Priority Setting, Resource Allocation; McGraw-Hill International Book Co: New York, NY, USA, 1980; p. 287. [Google Scholar]
  69. Pal, S.; Talukdar, S. Application of frequency ratio and logistic regression models for assessing physical wetland vulnerability in punarbhaba river basin of indo-bangladesh. Hum. Ecol. Risk Assess. Int. J. 2018, 24, 1291–1311. [Google Scholar] [CrossRef]
  70. Shahabi, H.; Hashim, M. Landslide susceptibility mapping using gis-based statistical models and remote sensing data in tropical environment. Sci. Rep. 2015, 5, 9899. [Google Scholar] [CrossRef] [PubMed]
  71. Akgün, A.; Bulut, F. GIS-based landslide susceptibility for Arsin-Yomra (Trabzon, North Turkey) region. Environ. Geol. 2007, 51, 1377–1387. [Google Scholar] [CrossRef]
  72. Shahabi, H.; Hashim, M.; Ahmad, B.B. Remote sensing and gis-based landslide susceptibility mapping using frequency ratio, logistic regression, and fuzzy logic methods at the Central Zab Basin, Iran. Environ. Earth Sci. 2015, 73, 8647–8668. [Google Scholar] [CrossRef]
  73. Günther, A.; Van Den Eeckhaut, M.; Malet, J.-P.; Reichenbach, P.; Hervás, J. Climate-physiographically differentiated pan-european landslide susceptibility assessment using spatial multi-criteria evaluation and transnational landslide information. Geomorphology 2014, 224, 69–85. [Google Scholar] [CrossRef]
  74. Malczewski, J. GIS and Multicriteria Decision Analysis; John Wiley & Sons: New York, NY, USA, 1999. [Google Scholar]
  75. Rahman, R.; Saha, S. Remote sensing, spatial multi criteria evaluation (SMCE) and analytical hierarchy process (AHP) in optimal cropping pattern planning for a flood prone area. J. Spat. Sci. 2008, 53, 161–177. [Google Scholar] [CrossRef]
  76. Abdollahi, S.; Pourghasemi, H.R.; Ghanbarian, G.A.; Safaeian, R. Prioritization of effective factors in the occurrence of land subsidence and its susceptibility mapping using an SVM model and their different kernel functions. Bull. Eng. Geol. Environ. 2018, 1–18. [Google Scholar] [CrossRef]
  77. Linden, A. Measuring diagnostic and predictive accuracy in disease management: An introduction to receiver operating characteristic (ROC) analysis. J. Eval. Clin. Pract. 2006, 12, 132–139. [Google Scholar] [CrossRef] [PubMed]
  78. Ghorbanzadeh, O.; Feizizadeh, B.; Blaschke, T. An interval matrix method used to optimize the decision matrix in AHP technique for land subsidence susceptibility mapping. Environ. Earth Sci. 2018, 77, 584. [Google Scholar] [CrossRef]
  79. Pourghasemi, H.; Gayen, A.; Park, S.; Lee, C.-W.; Lee, S. Assessment of landslide-prone areas and their zonation using logistic regression, logitboost, and naïvebayes machine-learning algorithms. Sustainability 2018, 10, 3697. [Google Scholar] [CrossRef]
  80. Baird, C.; Healy, T.; Johnson, K.; Bogie, A.; Dankert, E.W.; Scharenbroch, C. A Comparison of Risk Assessment Instruments in Juvenile Justice; National Council on Crime and Delinquency: Madison, WI, USA, 2013. [Google Scholar]
  81. Baeza, C.; Corominas, J. Assessment of shallow landslide susceptibility by means of multivariate statistical techniques. Earth Surf. Process. Landf. 2001, 26, 1251–1263. [Google Scholar] [CrossRef]
  82. Junior, F.R.L.; Osiro, L.; Carpinetti, L.C.R. A comparison between fuzzy ahp and fuzzy topsis methods to supplier selection. Appl. Soft Comput. 2014, 21, 194–209. [Google Scholar] [CrossRef]
  83. Feizizadeh, B.; Ghorbanzadeh, O. GIS-based interval pairwise comparison matrices as a novel approach for optimizing an analytical hierarchy process and multiple criteria weighting. GI_Forum 2017 2017, 1, 27–35. [Google Scholar] [CrossRef]
  84. Feizizadeh, B.; Blaschke, T. GIS-multicriteria decision analysis for landslide susceptibility mapping: Comparing three methods for the urmia lake basin, iran. Nat. Hazards 2013, 65, 2105–2128. [Google Scholar] [CrossRef]
  85. Shahabi, H.; Khezri, S.; Ahmad, B.B.; Hashim, M. Landslide susceptibility mapping at Central Zab Basin, Iran: A comparison BETWEEN analytical hierarchy process, frequency ratio and logistic regression models. Catena 2014, 115, 55–70. [Google Scholar] [CrossRef]
Figure 1. Map showing the location of the study area and field photographs: (a) Mailung Khola and (b) camp near Mailung Khola hydropower plant (c) road section near Ramche (d) road section near Syaprubesi.
Figure 1. Map showing the location of the study area and field photographs: (a) Mailung Khola and (b) camp near Mailung Khola hydropower plant (c) road section near Ramche (d) road section near Syaprubesi.
Ijgi 08 00094 g001
Figure 2. Landslide inventory map of the study area with the distribution of training and validation datasets.
Figure 2. Landslide inventory map of the study area with the distribution of training and validation datasets.
Ijgi 08 00094 g002
Figure 3. Thematic maps used in this study (a) land cover, (b) distance to roads, (c) elevation, (d) lithology, (e) faults, (f) drainage, (g) precipitation, (h) slope angle, and (i) slope aspect.
Figure 3. Thematic maps used in this study (a) land cover, (b) distance to roads, (c) elevation, (d) lithology, (e) faults, (f) drainage, (g) precipitation, (h) slope angle, and (i) slope aspect.
Ijgi 08 00094 g003aIjgi 08 00094 g003b
Figure 4. Output landslide susceptibility maps for each model; (a) analytical hierarchy process (AHP), (b) frequency ratio (FR), and (c) hybrid spatial multi-criteria evaluation (SMCE).
Figure 4. Output landslide susceptibility maps for each model; (a) analytical hierarchy process (AHP), (b) frequency ratio (FR), and (c) hybrid spatial multi-criteria evaluation (SMCE).
Ijgi 08 00094 g004
Figure 5. An example the performance of the three landslide susceptibility models in a sub-area of the case study area. (a) Frequency ratio; (b) analytical hierarchy process; (c) spatial multi-criteria evaluation.
Figure 5. An example the performance of the three landslide susceptibility models in a sub-area of the case study area. (a) Frequency ratio; (b) analytical hierarchy process; (c) spatial multi-criteria evaluation.
Ijgi 08 00094 g005
Figure 6. Receiver operating characteristics (ROC) representing quality method success rate curves for the three methods.
Figure 6. Receiver operating characteristics (ROC) representing quality method success rate curves for the three methods.
Ijgi 08 00094 g006
Figure 7. R-index in susceptibility classes for hybrid SMCE, AHP and FR.
Figure 7. R-index in susceptibility classes for hybrid SMCE, AHP and FR.
Ijgi 08 00094 g007
Figure 8. The comparison of the resulting weights for the classes of each factor based on the FR and AHP model.
Figure 8. The comparison of the resulting weights for the classes of each factor based on the FR and AHP model.
Ijgi 08 00094 g008
Table 1. Pairwise comparison point based rating scale [75].
Table 1. Pairwise comparison point based rating scale [75].
ImportanceDefinitionExplanation
1Equal importanceContribution to objective is equal
3Moderate importanceThe attribute is slightly favored over another
5Strong importanceThe attribute is strongly favored over another
7Very strong importanceThe attribute is very strongly favored over another
9Extreme importanceEvidence favoring one attribute is of the highest possible order of affirmation
2, 4, 6, 8Intermediate valuesWhen compromise is needed
Table 2. The pairwise comparison matrix of the hybrid spatial multi-criteria evaluation (SMCE) groups and the resulting weights.
Table 2. The pairwise comparison matrix of the hybrid spatial multi-criteria evaluation (SMCE) groups and the resulting weights.
GroupsTopographyGeologicalHydrologicalManmadeWeights
Topography1---0.30
Geological11--0.30
Hydrological0.830.831-0.25
Manmade0.50.50.610.15
Table 3. Weights for groups and factors based on the hybrid SMCE model.
Table 3. Weights for groups and factors based on the hybrid SMCE model.
GroupsHybrid SMCE Group WeightsFactorsHybrid SMCE Factor Weights
Topography30DEM0.3
Slope0.36
Aspect0.34
Geological30Faults0.4
Lithology0.6
Hydrological25Drainage0.5
Precipitation0.5
Manmade15Land use0.8
Distance to roads0.2
Table 4. Pairwise comparison matrix for the analytical hierarchy process (AHP) model.
Table 4. Pairwise comparison matrix for the analytical hierarchy process (AHP) model.
FactorsClassesPairwise Comparison MatricesEigenvaluesCRNumber of Landslides
Land useForest1 0.37 11,854
Shrub land1/51 0.166 782
Grassland1/31/31 0.126 4803
Agriculture1/31/311 0.153 5555
Barren1/51/21/31/51 0.063 316
Waterbody1/51/21/31/31/21 0.053 31
Snow cover1/511/21/31110.068 97
0.092
Precipitation (mm)950–17251 0.067 1349
1725–250021 0.147 5296
2500–3275731 0.493 9456
3275–4050521/21 0.291 7275
0.006
LithologyGlacier1 0.034 8
Fluvial calcareous21 0.062 79
Fluvial non-calcareous321 0.098 260
Slate/ phyllite5321 0.161 1151
Quartzite84321 0.293 2429
Gneiss954311 0.35 19,512
0.009
Distance to fault (m)(1) 0–20001 0.641 4161
(2) 2000–40001/21 0.221 3471
(3) 4000–60001/31/31 0.086 2690
(4) 6000 <1/21/211 0.050 13,117
0.03
Distance to drainage (m)<2001 0.41 8942
200–4001/21 0.254 4708
400–6001/31/21 0.152 3047
600–8001/41/31/21 0.104 2038
>8001/41/31/21/21 0.078 4703
0.032
Slope (%)0–101 0.053 733
10–2031 0.067 2380
20–30871 0.235 5600
30–409831 0.325 8922
>4098331 0.320 5773
0.054
Elevation (m)(1) <10001 0.133 1704
(2) 1000–300051 0.566 18,520
(3) 3000–450021/31 0.206 3001
(4) > 450011/51/21 0.093 188
0.024
Aspect(1) Flat1 0.195 1303
(2) North1/21 0.172 4419
(3) East211 0.291 9219
(4) West111/221 0.193 5399
(5) South111/21/211 0.147 3068
0.058
Distance to roads (m)(1) 0–2001 0.527 11,316
(2) 200–4001/31 0.315 5403
(3) 400–6001/51/51 0.095 2634
(4) 600<1/61/61/21 0.061 4085
0.069
Table 5. The normalized weights for the classes of the factors based on the frequency ratio (FR) model.
Table 5. The normalized weights for the classes of the factors based on the frequency ratio (FR) model.
ClassesPixels of Each Class% of PixelsLandslide Pixels% of PixelsFR
Land useForest6,769,54642.227,357,50050.630.207
Shrub land436,3822.72495,9003.410.216
Grassland2,428,46115.143,038,40020.900.238
Agriculture3,691,18323.023,378,60023.250.174
Barren951,4855.93182,7001.250.037
Waterbody34,9930.2122,5000.150.122
Snow cover1,718,69010.7255,8000.380.006
Precipitation (mm)950–1725200,43913.874,824,73638.0300.641
1725–2500409,97228.372,546,43220.0720.166
2500–3275576,68539.911,629,93612.8480.075
3275–4050258,00817.851,135,2328.9480.117
LithologyGlacier94,8290.5145000.030.021
Fluvial calcareous341,5481.8543,2000.300.055
Fluvial non-calcareous427,6762.31157,5001.080.160
Slate/phyllite3,334,36618.04717,3004.940.093
Quartzite3,671,63719.861,473,30010.150.174
Gneiss10,615,03857.4212,123,00083.500.496
Distance to fault (m)(1) 0–20001,378,8397.952,583,90017.770.353
(2) 2000–40001,361,0417.852,114,10014.540.293
(3) 4000–60001,329,3247.671,678,50011.550.238
(4) 6000<13,273,15076.548,160,30056.140.116
Distance to drainage (m)<2003,285,08317.834,824,73638.030.362
200–4002,625,57014.252,546,43220.070.239
400–6002,280,17712.381,629,93612.850.176
600–8001,907,24510.351,135,2328.950.147
>8008,323,42445.182,550,35220.100.076
Slope (%)0–101,158,7246.76425,0553.050.076
10–203,806,96722.211,432,08310.270.078
20–306,322,93936.893,342,67823.980.110
30–404,377,81725.545,345,52538.350.254
>401,474,3048.603,393,54624.350.480
Elevation (m)(1) <10002,976,27517.331,030,3077.390.152
(2) 1000–30007,962,47846.3710,985,91678.800.607
(3) 3000–45003,586,95620.891,816,61613.030.223
(4) >45002,646,44015.41109,4970.790.018
Aspect(1) Flat1,899,66011.081,480,36610.620.113
(2) North3,304,52219.283,542,70425.420.279
(3) East4,238,42324.734,625,60433.180.298
(4) West4,004,39123.362,661,55419.090.190
(5) South3,693,75521.551,628,66111.680.119
Distance to roads (m)(1) 0–2009,064,14049.0986,143,42448.390.234
(2) 200–4003,973,07821.5214,304,16033.910.374
(3) 400–6001,237,4346.7031,048,9928.260.293
(4) 600<4,186,51622.6771,197,9529.440.099
Table 6. Resulting R-indexes for the landslide susceptibility mappings (LSMs) based on AHP, FR and hybrid SMCE models.
Table 6. Resulting R-indexes for the landslide susceptibility mappings (LSMs) based on AHP, FR and hybrid SMCE models.
ModelsSusceptibility ClassNumber of PixelsArea (m²)Area Percentage (ni)Number of LandslidesLandslide Percentage (Ni)R- Index
Hybrid SMCEVery Low53,100481,282,2003.351300.553
Low656,1002,926,548,00020.357553.222
Moderate4,664,7006,508,665,00045.26532022.707
High11,235,6003,814,145,00026.5212,78054.5229
Very high3,942,000648,798,3004.51445419.0059
AHPVery Low6300167,919,3001.17710.305
Low64,8001,554,764,00010.81740.321
Moderate3,978,9006,563,949,00045.65452819.329
High9,473,4004,198,336,00029.2010,81646.1532
Very high7,028,1001,894,469,00013.17794833.9153
FRVery Low61,200154,024,2001.07690.294
Low998,1003,768,698,00026.1111504.913
Moderate3,970,8005,621,287,00038.94450419.228
High9,657,9003,801,150,00026.3311,01546.9927
Very high5,905,8001,090,157,0007.55670128.5958

Share and Cite

MDPI and ACS Style

Meena, S.R.; Ghorbanzadeh, O.; Blaschke, T. A Comparative Study of Statistics-Based Landslide Susceptibility Models: A Case Study of the Region Affected by the Gorkha Earthquake in Nepal. ISPRS Int. J. Geo-Inf. 2019, 8, 94. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi8020094

AMA Style

Meena SR, Ghorbanzadeh O, Blaschke T. A Comparative Study of Statistics-Based Landslide Susceptibility Models: A Case Study of the Region Affected by the Gorkha Earthquake in Nepal. ISPRS International Journal of Geo-Information. 2019; 8(2):94. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi8020094

Chicago/Turabian Style

Meena, Sansar Raj, Omid Ghorbanzadeh, and Thomas Blaschke. 2019. "A Comparative Study of Statistics-Based Landslide Susceptibility Models: A Case Study of the Region Affected by the Gorkha Earthquake in Nepal" ISPRS International Journal of Geo-Information 8, no. 2: 94. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi8020094

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop