Next Article in Journal
Review of Physical Layer Security in 5G Wireless Networks
Next Article in Special Issue
Dynamic Monitoring and Analysis of Mining Land Subsidence in Multiple Coal Seams in the Ehuobulake Coal Mine Based on FLAC3D and SBAS-InSAR Technology
Previous Article in Journal
Static Analysis of Mobile Pump Truck Frame for Four Typical Working Conditions
Previous Article in Special Issue
Characterizing Crustal Deformation of the Weihe Fault, Weihe Basin (Central China), Using InSAR and GNSS Observations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploration and Comparison of the Effect of Conventional and Advanced Modeling Algorithms on Landslide Susceptibility Prediction: A Case Study from Yadong Country, Tibet

1
Guangzhou Urban Planning & Design Survey Research Institute, Guangzhou 510060, China
2
Guangdong Enterprise Key Laboratory for Urban Sensing, Monitoring and Early Warning, Guangzhou 510060, China
3
Guangzhou Collaborative Innovation Center of Natural Resources Planning and Marine Technology, Guangzhou 510060, China
4
College of Construction Engineering, Jilin University, Changchun 130012, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(12), 7276; https://doi.org/10.3390/app13127276
Submission received: 16 May 2023 / Revised: 13 June 2023 / Accepted: 13 June 2023 / Published: 19 June 2023
(This article belongs to the Special Issue Remote Sensing Technology in Landslide and Land Subsidence)

Abstract

:
Shallow landslides pose serious threats to human existence and economic development, especially in the Himalayan areas. Landslide susceptibility mapping (LSM) is a proven way for minimizing the hazard and risk of landslides. Modeling as an essential step, various algorithms have been applied to LSM, but no consensus exists on which model is most suitable or best. In this study, information value (IV) and logistic regression (LR) were selected as representatives of the conventional algorithms, categorical boosting (CatBoost), and conventional neural networks (CNN) as the advanced algorithms, for LSM in Yadong County, and their performance was compared. To begin with, 496 historical landslide events were compiled into a landslide inventory map, followed by a list of 11 conditioning factors, forming a data set. Secondly, the data set was randomly divided into two parts, 80% of which was used for modeling and 20% for validation. Finally, the area under the curve (AUC) and statistical metrics were applied to validate and compare the performance of the models. The results showed that the CNN model performed the best (sensitivity = 79.38%, specificity = 91.00%, accuracy = 85.28%, and AUC = 0.908), while the LR model performed the worst (sensitivity = 79.38%, specificity = 76.00%, accuracy = 77.66%, and AUC = 0.838) and the CatBoost model performed better (sensitivity = 76.28%, specificity = 85.00%, accuracy = 80.81%, and AUC = 0.893). Moreover, the LSM constructed by the CNN model did a more reasonable prediction of the distribution of susceptible areas. As for feature selection, a more detailed analysis of conditioning factors was conducted, but the results were uncertain. The result analyzed by GI may be more reliable but fluctuates with the amount of data. The conclusion reveals that the accuracy of LSM can be further improved with the advancement of algorithms, by determining more representative features, which serve as a more effective guide for land use planning in the study area or other highlands where landslides are frequent.

1. Introduction

Landslides are a geological phenomenon that can cause extensive property damage and pose a threat to the safety of residents worldwide [1,2,3]. Notably, China experiences a higher frequency and larger scale of landslides than any other country [4,5]. Effective prevention measures require the identification and mapping of existing landslides, with a focus on predicting potential disaster areas [6,7]. To minimize the hazard and risk of landslides, landslide susceptibility mapping (LSM) is widely recognized as a proven approach.
The precision of LSM is contingent upon dependable data and modeling algorithms. Throughout time, a range of algorithms have been utilized in LSM, allowing for incremental improvements in accuracy. Both quantitative (data-driven) and qualitative (knowledge-based or physically based) algorithms are viable options for landslide susceptibility modeling [8,9]. Data-driven algorithms are usually classified as bivariate methods (such as certainty factor (CF) and information value (IV)), multivariate methods (such as logistic regression (LR) and cluster analysis), conventional machine learning [10,11,12,13], and deep learning [14]. Notwithstanding the fact that knowledge-based methods tend to be subjective and confined to small-scale domains, physically based methods, which rely on physical experiments, often entail significant time and financial resources and exhibit limited scalability. In recent years, the advent of advanced remote sensing has facilitated the acquisition of landslide samples and critical factors, rendering data-driven approaches more accessible and precise. Statistical methods are sensitive to the distribution and collinearity of data, and certain assumptions must be satisfied before application. There have been a number of conventional machine learning methods (such as random forest and categorical boosting (CatBoost) observed due to their ability to solve non-linear geo-environmental problems without making unnecessary assumptions [15,16]. The structure of conventional machine learning methods is relatively simple, usually involving only 1~2 hidden layers. While certain studies have undertaken a thorough examination and comparison of conventional statistical and machine learning models, a consensus has yet to be reached regarding the optimal or most appropriate model [17,18].
Deep learning represents a progression beyond machine learning, exhibiting superior performance in identifying representative features and potentially enhancing classification accuracy [19,20,21]. Convolutional neural networks (CNN), a prominent example of deep learning, have been extensively validated in domains such as image classification and facial recognition, yet their application to LSM remains limited [22,23,24]. A multi-level structure of CNN is constructed to explore the complex non-linear relationships between variables, which is accorded with the characteristics of LSM.
The superiority of the study is to assess the impact of algorithms on prediction accuracy and feature recognition by comparing the performance of conventional and advanced algorithms, as well as four representative models (IV, LR, CatBoost, and CNN). Yadong County in Southeastern Tibet was selected as the study area because of its topographic and geological conditions, resulting in frequent shallow landslides. Four models, IV, LR, CatBoost, and CNN, were explored and compared for the effect on landslide susceptibility prediction.

2. Materials

2.1. Study Area

Yadong County, situated at the southern foot of the middle section of the Himalayas, falls under the jurisdiction of Shigatse City, Tibet Autonomous Region (Figure 1). The study area, spanning over 4240 km2, is home to a population of over 120,000 individuals. The region is characterized by a high mountain landform, with an average elevation of 3500 m (ranging from 1747~7057 m). The study area is geographically divided into two distinct parts, namely north and south, based on the line from Pali Town to Kangbu Town, owing to the significant difference in altitude. Consequently, two distinct climates have emerged in these regions [5]. The study area exhibits distinct geographical and geological characteristics. The northern region is situated at an elevation of 4300 m above sea level and experiences a cold and dry climate, with an average annual rainfall of approximately 410 mm. Conversely, the southern region has a lower average altitude of 2800 m, a humid climate, and an annual precipitation of 873 mm. The geological composition of the study area is primarily composed of shale, limestone, and dolomite, with the southern region being characterized by faults and folds. There is a degree of VIII on the modified Mercalli index, indicating a high seismic intensity. The southern area also exhibits a higher density of geological disasters, such as collapses and landslides, while the northern region is primarily affected by debris flows. Figure 1 also shows the location of landslides that occurred over the years and roads distributed in the study areas.
On 18 September 2011, a magnitude 6.8 earthquake occurred in Sikkim, India. The study area was affected by the earthquake, which induced a series of secondary geological disasters, causing great harm to local residents and roads (Figure 2 and Figure 3). Therefore, it is of great significance to compile an accurate landslide susceptibility map for the study area.

2.2. Data Preparation

2.2.1. Landslide Inventory

Typically, data-driven methods for LSM assume that landslides have a greater chance of recurring under the same conditions as they did before [25]. Therefore, a comprehensive and exact landslide inventory that shows the locations and numbers of landslides is essential [7]. It is therefore essential to develop a complete and accurate landslide inventory (Figure 4) that shows the locations and quantities of landslides. Landslides were obtained by aerial photographs, collecting literature and historical reports, and conducting extensive fieldwork. Landslides are bounded by polygons containing the entire perimeter, and 496 polygons representing the landslide perimeter were identified. The landslide locations are the centroids of the polygons (Figure 1). Landslides are mainly small- and medium-sized shallow landslides with volumes ranging from 100 m3 to 1000 m3. Using a 1:1 sampling strategy [26], the model was trained and validated with 992 samples, including 496 landslides with a sign of “1” and 496 non-landslides with a sign of “0”. Non-landslide samples were randomly selected from non-landslide dense areas.

2.2.2. Conditioning Factors

A total of 11 condition factors were selected based on the study region’s characteristics, the data’s availability, reliability, and practicality, including elevation, slope, maximum elevation difference (MED), plan curvature, profile curvature, topographic wetness index, distance to faults, distance to roads, distance to streams, annual rainfall, and NDVI [27].
Elevation, which was divided into five sub-classes by 1000 m intervals, has an influence on both rainfall and vegetation [28,29] (Figure 4a). Slope, which controls shear strength at potential slide surfaces as well as subsurface flow [30], was reclassified into six classes by 10° intervals (Figure 4b). It was calculated in ArcGIS 10.2 that the MED reflects potential kinetic energy of slope units. By 100 m intervals, the thematic map was reclassified into six classes (Figure 4c). Erosion and deposition can be determined from slope curvature since it is essential to the geometry of slopes [31]. Six classes were established for profile and plan curvatures (Figure 4d,e). The TWI represents basic terrain, which was divided into six categories (Figure 4f). With the spatial resolution of 30 m, six topographic factors were extracted from the digital elevation model (DEM). The DEM was originally sourced from the Geospatial Data Cloud, which was accessed on 6 January 2023.
Bulk-rock strength could be reduced as a result of faults acting as potential weak planes in slopes, and distance to fault was constructed with six classes by 2000 m intervals (Figure 4g). In a similar manner, distances to the road and river were divided into six categories with an interval of 2000 m (Figure 4h,i).
An existing 1:50,000 geological map was used to extract fault information, and Landsat 8 LOI images were used to determine road and river network information. Using the Euclidean Distance ArcGIS Tool, the distances from each raster unit of the area to the nearest fault, road, stream, and highway were calculated.
The unique trigger factor considered in the study was rainfall, and it has been applied numerous times. Based on the data from 14 precipitation stations near the study area, an ordinary kriging interpolation was used to generate the thematic map in ArcGIS. The thematic map was reclassified into 4 classes (Figure 4j).
The consolidation of vegetation roots can stabilize soil, which alleviates the effect of rainfall on the stability of the slope [32]. The vegetation normalization index (NDVI) was applied to evaluate the vegetation, and it is calculated by:
N D V I = N I R R E D N I R + R E D
where RED is the reflection value of the red wavelength and NIR is the near-infrared wavelength.
NDVI ranges from −1~1 with a positive value indicating coverage, and the greater the value, the greater the coverage. The thematic map was reclassified into 6 classes (Figure 4k). The NDVI was derived from Landsat 4–5 TM satellite images.

2.2.3. Mapping Units

LSM commonly uses grid cells, slope units, and unique-condition units as mapping units. The choice of mapping unit is controversial; however, grid cells are the most common [33,34]. Landslides can be represented better by a slope unit, which represents their source, transport, and accumulation areas [35,36,37]. Thus, the study area is divided into 25,483 slope units with the hydrologic analysis tool in ArcGIS and necessary artificial corrections are accompanied according to remote sensing images.

3. Methods

The main aims of the study are to explore and compare the effect of conventional and advanced modeling algorithms on landslide susceptibility prediction. The methodology followed in this study mainly contains four steps (Figure 5). Firstly, data sets including landslide samples, conditioning factors, and mapping units are prepared for modeling. The second step is to divide the data set into training and validation parts. After that, four representative algorithms, IV, LR, CatBoost, and CNN, are applied to LSM. Finally, the performance of LR, CatBoost, and CNN models are analyzed and compared based on some key measures.

3.1. IV

By using the frequency or density of landslides, the IV method reflects the magnitude of the hazards associated with different influencing factors and their sub-intervals. IV was first proposed by Yin and Yan (1988) and later modified by Van Westen (1993) [38]. Equation (1) shows the method for the calculation of the information values:
I ( A i j ) = ln s i j / S n i j / N
where i = 1, 2, 3, …, n; j = 1, 2, 3, …, m; si−j represents the area of the landslide of the i-th conditioning factor in j-th interval; ni−j represents the area of the i-th conditioning factor in j-th interval; S represents the total area of the landslide; N represents the total area.
Depending on the IV method, the value might be positive or negative. When IV values are positive, they indicate that the factor is conducive to the occurrence of landslides in a particular interval; a higher IV indicates a higher likelihood of landslides, and vice versa.
The total information value I can be determined by:
I i = i = 1 n l n s i / S n i / N
Finally, the information value calculated by each unit is processed by linear normalization.

3.2. LR

The LR model is used for statistical analysis of binary dependent variables (the dependent variable y has two values: 0 and 1) [39]. The LR model is advantageous since the data distribution can either be nominal or continuous [40]. It is computed using the following equation:
P = 1 1 + e y
where P represents the occurrence probability of an event (such as a landslide) on the s-curve in the range of 0 to 1; y represents a linear combination function and is calculated by using Equation (5).
y = b 0 + b 1 x 1 + b 2 x 2 + b 3 x 3 + b n x n
where b0 is the constant value or intercept of the equation, and b1, b2, ..., bn are the regression coefficients of the explanatory variables x1, x2, ..., xn.
In the current study, the LR model was performed in SPSS software, and the forward stepwise method was adopted to screen valuation indexes. The conditioning factors were calculated as independent variables, whereas dependent variables represent the occurrence of landslides. During the last step of the analysis, all variables were significant at less than 0.05, so no additional variables were included.

3.3. CatBoost

CatBoost is an improved implementation based on the Gradient Boosted Decision Tree (GBDT) framework, with fewer parameters. It is an open-source machine-learning library of Yandex and a member of the Boosting family. The gradient deviation and prediction offset can be effectively solved by the CatBoost algorithm, so as to reduce the occurrence of over-fitting and improve the accuracy and generalization ability of the algorithm [41]. It is worth noting that the selection of different random numbers can have an impact on the predictive outcome of the model. Detailed information about CatBoost can be found in other literature [42].

3.4. CNN

CNN provides an end-to-end learning model that can select model parameters through traditional gradient descent methods. The trained CNN can learn image features and complete the extraction and classification of image features [43]. A CNN consists of the network structures as convolutional layer, down-sampling layer, and fully connected layer, each of which contains a number of independent neurons [44]. Its notable feature is that the weight sharing and local connection of the convolution kernel in the hidden layer can greatly reduce the number of weights, thereby reducing the complexity of the convolutional network model.
The unique convolutional layer and pooling layer structure of CNN can be combined arbitrarily to obtain an infinite variety of network models. This study applied the AlexNet model (including 4 convolutional layers and 3 fully connected layers) to LSM [45]. Before modeling, conditioning factors are requested to be superimposed to obtain a multi-band two-dimensional image as the input raw data in LSM. Since LSM belongs to a two-classification problem of landslide disasters, two neurons are placed in the output layer classifier to represent landslide and non-landslide. All parameters have been optimized through a trial-and-error approach. The number of training cycles is 10, the initial learning rate is 0.01, the loss function uses standard cross-entropy, the optimizer is Adam, and the activation function is Relu.
The machine learning and deep learning-based algorithms are implemented in Python3.7 based on the Package of Numpy, Scikit-Learn, and Tensorflow. The statistics-based algorithms are implemented in SPSS.

3.5. Models Evaluation

A predictive model will not be convincing without scientific validation, so existing data will need to be split into training and validation sets. Previous research derived performance measures for single-hold-out models by employing a single training set and an independent test set. In the study, 80% of the data sets were randomly chosen for training and 20% for validation in LSM, and the split was repeated 5 times. Finally, these models were evaluated using the average of 5 groups.
Four evaluation measures—accuracy, sensitivity, specificity, and AUC—were combined to analyze the performance of the models. Accuracy, sensitivity, and specificity are calculated from the Confusion Matrix, which is an N × N matrix (Table 1). The TP represents the number of landslides that have been correctly predicted as unstable, and TN represents the number of non-landslides that have been correctly predicted as stable, while FP represents the number of non-landslides that have been predicted incorrectly as unstable and FN represents the number of landslide units that have been predicted incorrectly as stable.
An additional indicator of model validity is the area under the receiver operating characteristic curve (AUROC). The value of AUROC ranges from 0.5~1; the larger the value, the better the generalization ability of the model and prediction performance. Finally, the natural breaks method was applied to reclassify landslide susceptibility into five classes: very low, low, moderate, high, and very high.

4. Results

4.1. Performance and Comparison of Conventional and Advanced Algorithms

Before modeling, the data were normalized with Z-scores to eliminate the effects of different dimensions. Additionally, a correlation analysis was conducted to test the collinearity among the independent variables using the variance inflation factor (VIF) [44]. There is severe collinearity between the selected variables if the VIF value exceeds 10. Table 2 showed the VIF values of the chosen independent variables, and the factor elevation has the highest VIF value (5.711). The result indicates that no severe collinearity problems exist among the chosen variables, and thus 11 conditioning factors were applied to the modeling.
The performance of the three models using the confusion matrix is shown in Table 3 As for training, the CNN model achieved the highest value of sensitivity (84.21%), followed by the CatBoost model (sensitivity = 82.96%) and the LR model (sensitivity = 81.45%). As for specificity, the CNN model performed best with 93.52%, followed by the CatBoost model at 86.63% and the LR model at 84.79%. The CNN model had the best accuracy and ROC values of 88.88% and 0.944, while the LR model had the worst values of 83.13% and 0.897. As well, CatBoost performed well, scoring 86.63% and 0.930, respectively.
The verification data set is more useful and important for evaluating the ability of these models to generalize. In Figure 6, we found that the CNN model had the highest sensitivity, specificity, accuracy, and AUC values, namely 79.38%, 91.00%, 85.28%, and 0.908. Additionally, CatBoost performed well with 76.28%, 85.00%, 80.71%, and 0.893, respectively. The LR model remained the worst with values of 79.38%, 76.00%, 77.66%, and 0.838 (Table 3).
Model performance declined in verification, particularly for the LR and CatBoost models (the accuracy value decreased by about 6%), which indicated that the models were overfitting and generalizability was suspicious. Meanwhile, the performance of the CNN model was stable as the value of AUC reached 0.908, and the accuracy value was close to the training data set.

4.2. Evaluation of Conditioning Factors

4.2.1. Application of Conventional Algorithms

Bivariate methods were used to establish relationships between the conditioning factors and the occurrence of the landslides. From Table 4, as for elevation, the percentages of landslide area for 3700~4700 m and 4700~5700 m were 43.80% and 28.63%, respectively, which means that over 70% of landslide areas were distributed among the two classes but the IV value of these two classes were both negative, while the highest IV of elevation is 1.098 for the class of >5700. For slope, the IV values increased as the slope increased, as well as the factor MED, which indicated that landslides are more likely to occur in steep areas. The IV of slope ranged from −1.106~1.709, in which the classes of >20 were all positive. In terms of TWI, the class of 6~7 had the highest IV of 0.371 and the IV were all negative at the class of >8.
For curvature, the highest IV of the plan curvature was 0.538, located at the class of <−0.4, and 1.106 of profile curvature at the class of >0.3. It was found that landslides were concentrated in convex terrain. On the other hand, landslides are highly concentrated near the faults, streams, and roads as the IV reached the maximum value at first class. It indicated that the development of faults, streams, and human engineering activities was conducive to the occurrence of landslides.
As for rainfall, the highest probability of landslides occurring appeared in the class of >780 mm (IV reached 1.107). Meanwhile, the IV did not increase with increasing rainfall and it indicated that landslide occurring is complex. In the case of NDVI, the class of 0.6~0.75 has the highest IV of 0.804. The IV changed erratically for NDVI although high vegetation cover helped enhance the stability of the landslide.
As for the LR model, the final regression equation is as follows:
y = ( 1.756   M E D ) + ( 0.395   p r o f i l e . c u r v a t u r e ) + ( 1.799 rainfall ) + ( 1.641   D T R ) + ( 0.684   D T S ) + 0.28
It can be found that MED, rainfall, and DTR were essential to landslide occurring as the coefficients were relatively large. While profile curvature and DTS were considered as secondary factors. The coefficient of DTR was negative, which indicated a negative effect on landslide occurring.
Conventional algorithms such as IV and LR methods were used to establish the relationships between the conditioning factors and the occurrence of the landslides in this study. In summary, the conditioning factors have an effect on landslides occurring and the impact of different factors on landslides varies in different intervals.

4.2.2. Application of Advanced Algorithms

Actually, conditioning factors have different influences on the occurrence of landslides while the bivariate methods fail to recognize the difference. Gini importance (GI) is defined as the total reduction in average nodal impurities across all trees. GI was also applied to analyze the relative importance of different factors, exploring the factor’s contribution to landslide occurring.
The bigger the GI express greater the importance of factors of landslides occurring [44]. Table 5 showed the rank of the main conditioning factors of this study. The result showed that DTF was the most important factor responsible for landslides as the GI reached 5.36. Moreover, the factor MED, DTS, elevation, plan curvature, and slope were also pivotal since the GI were all greater than zero. While the other factors as TWI, rainfall, profile curvature, NDVI, and DTR showed little contribution to landslides in the study area (Figure 7).

4.3. Landslide Susceptibility Mapping Results

In the case of LR model, LSI ranges from 0.11 to 0.96, and the corresponding area percentages were 40.17% (very low), 12.24% (low), 5.54% (moderate), 17.75% (high), and 24.30% (very high), respectively. Similarly, five reclassified classes of the CatBoost model accounted for 30.16%, 23.81%, 15.21%, 26.82%, and 14.94%, respectively, of the entire area. LSM constructed by the CNN model was also divided into five classes: very low (<0.18), low (<0.42), moderate (<0.653), high (<0.82), and very high (>0.82), accounting for 28.99%, 19.65%, 14.21%, 9.19%, and 27.96% of the whole area. It is noticed that LR, CatBoost, and CNN models predicted the largest proportion of very low susceptibility while IV predicted the largest proportion of very high susceptibility.
A logical landslide susceptibility map should meet two rules: (1) the density of landslide samples should increase with the increase of susceptibility class and be mainly located in the highest susceptibility area; (2) the landslide susceptibility map should be spatially continuous and smooth, and the very high-susceptibility class area should occupy a small proportion. Concomitant with these maps (Figure 8), the landslide samples (dark spots) were mainly located in the red areas, and the non-landslide samples (blue spots) were in the green areas. The high or very high susceptibility areas were concentrated in the south of the areas, which was consistent with the distribution of landslides (Figure 9). Thus, the maps predicted by these models were logical on the whole.
There are distinct differences among the landslide susceptibility maps derived from the models. Table 6 showed that the IV increased as the susceptibility class increased for the model and the IVs were all negative for the very low, low, and moderate susceptibility class while positive for the very high susceptibility class. For high susceptibility, the IVs were positive for LR, CatBoost, and CNN models while negative for the IV model. On the other hand, the maps predicted by IV and LR models were spatially discontinuous while CatBoost and CNN models produced smoother patterns. Moreover, the percentage of very high susceptibility class was similar for LR, CatBoost, and CNN models, while the IV model predicted the highest percentage, reaching 37.27%.

5. Discussion

With the progress of related technologies and the maturity of theory, various approaches have been developed and applied to LSM to improve prediction accuracy and reliability [45,46,47,48,49]. Recently, deep learning is becoming popular and beginning to be used for LSM. In this study, CNN was applied to LSM, and its performance was compared to the conventional statistical methods and machine learning approaches.
In comparison to the other three models, the CNN model shows superior fitting and generalization capabilities in predicting landslide susceptibility. It has been found that very high landslide susceptibility areas are mainly associated with the Yarlung Zangbo River and its tributaries. A large part of the eroded slopes is scoured by the river network. The areas near streams are densely populated, and the occurrence of landslides can threaten lives.
Different algorithms have different emphases and, generally, their performance varies with different study areas [50,51]. However, advanced algorithms usually perform better in terms of accuracy compared to conventional statistical methods [52,53]. The result of our study also found that CNN and CatBoost performed better in terms of accuracy and AUC, and CNN did the best in generalization. There was a certain gap between the three models. The improvement benefits from the characteristics of the algorithms themselves for decreasing the bias, discrepancy, and over-fitting problems. It is easy to implement and acceptable in conventional statistical methods to establish a mathematical equation for investigating the relationship between factors related to landslides and landslides occurring. In machine learning methods, optimization is stressed, so the multiple parameters involved need to be tuned before application, which is unattainable for non-experts although some algorithms for optimizing parameters have been applied [54,55]. Deep learning algorithms use a more complex modeling architecture consisting of convolutional, activation, pooling, and fully connected layers, taking images as input parameters. Feature selection and information filtration are performed in the pooling layer, which is a robust step in CNN. Importantly, the dimensionality reduction is performed without changing the depth of the maps. Thus, deep learning algorithms have the ability to process data efficiently, feature extraction of high-dimensional data, and keep high prediction accuracy.
The accuracy of LSM should not be the only priority. Identifying the major conditions that lead to landslides is also critical, which assists in furthering the process [56,57]. Identifying subjective weights and objective weights allows for separating their contributions. Analytic hierarchy processes (AHP) and factor analyses (FA) are commonly used methods [58,59]. Landslides and their underlying conditioning factors can be directly correlated using LR models [60,61]. The relative importance of factors can be determined by the magnitude, plus, or minus of the coefficients. In addition, bivariate methods are capable of distinguishing factors with different susceptibilities to landslides across interval ranges. Accordingly, LR and IV are recommended to be combined to analyze the factor’s contribution to a landslide occurring. GI describes the contribution of the conditioning factors by calculating the total decrease, and the final coefficients reflect the relative importance of different factors. Some studies have applied GI for feature selection, which helps to decrease redundant information [62,63].
As for bivariate methods such as IV, FR, and CF, the performance is mainly up to whether the division of the conditioning factors interval is reasonable. Some factors such as DTF and DTR in the study may show a logical distribution while the NDVI showed a chaotic distribution. However, there is no consensus on the size or number of intervals, although some methods such as natural breaks and equal intervals have been applied. Moreover, the total IV is a linear addition of the IV of all conditioning factors, which further amplifies the uncertainty of the final result. Thus, the susceptibility analysis of conditioning factors on different intervals by bivariate methods may unreliable, and the landslide susceptibility maps are also difficult to be verified quantitatively. In terms of LR, the establishment of the equation derives from the distribution characteristics of data and is sensitive to linear correlation. Thus, the performance is up to the data partition and will fluctuate during validation. The results of relative importance for conditioning factors may be confused or even contrary to our experience on landslides. Machine learning or deep learning emphasizes iterative operations and repeated verification and requires a large amount of data. Landslide samples are limited in a restricted area and difficult to collect. The performance will also decline as the data decrease.
Moreover, lithology as a categorical variable is a crucial factor used by practically most of the researchers of LSM. The geological composition of the study area is primarily composed of shale, limestone, and dolomite. The lithology distribution in the study area is simple. Thus, we select the distance to faults as a more important factor.

6. Conclusions

In the current study, IV and LR were selected as representatives of the conventional algorithms, with CatBoost and CNN as the advanced algorithms, for LSM in Yadong county, and their performance was compared. The following conclusions can be drawn:
  • There was a certain gap between the models. Compared to conventional algorithms, advanced algorithms performed better in terms of prediction accuracy and CNN performed the best in generalization, thus it is regarded as the best model in this study.
  • The landslide susceptibility map predicted by CNN was more reasonable and the very high susceptibility areas were mainly distributed along the Yarlung Zangbo River.
  • As for feature selection, IV and LR performed a more detailed analysis of conditioning factors, but the results were uncertain. The result analyzed by GI may be more reliable but fluctuates with the amount of data.
  • The conventional algorithms are inferior to the advanced algorithms in accuracy and feature selection, but conventional algorithms have better resolvability and operability.
However, there are also some limitations of the present study:
  • There are possibilities for the combination of conventional and advanced algorithms, and further exploration is needed to improve prediction accuracy obviously.
  • Models need to be validated more reliably.

Author Contributions

Z.L., writing—original draft, methodology, and software; W.P., W.L. and H.H., review and validation; J.H. and K.L., reviewing and editing. K.J. and G.L. are responsible for the review and validation. All authors have read and agreed to the published version of the manuscript.

Funding

The completion of this work was supported by Guangzhou Collaborative Innovation Center of Natural Resources Planning and Marine Technology (No. 2023B04J0301), Key-Area Research and Development Program of Guangdong Province (No. 2020B0101130009), Guangdong Enterprise Key Laboratory for Urban Sensing, Monitoring and Early Warning (No. 2020B121202019), The Science and Technology Foundation of Guangzhou Urban Planning and Design Survey Research Institute (Grant No. RDI2220204128, RDI2220204031, RDI2220204037), and Postdoctoral Research Project of Guangzhou (20220402).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Das, I.; Sahoo, S.; van Westen, C.; Stein, A.; Hack, R. Landslide susceptibility assessment using logistic regression and its comparison with a rock mass classification system, along a road section in the northern Himalayas (India). Geomorphology 2010, 114, 627–637. [Google Scholar] [CrossRef]
  2. Achour, Y.; Boumezbeur, A.; Hadji, R.; Chouabbi, A.; Cavaleiro, V.; Bendaoud, E.A. Landslide susceptibility mapping using analytic hierarchy process and information value methods along a highway road section in Constantine, Algeria. Arab. J. Geosci. 2017, 10, 194. [Google Scholar] [CrossRef]
  3. Arabameri, A.; Pradhan, B.; Rezaei, K.; Lee, S.; Sohrabi, M. An ensemble model for landslide susceptibility mapping in a forested area. Geocarto Int. 2020, 35, 1680–1705. [Google Scholar] [CrossRef]
  4. Ni, H.Y.; Zheng, W.M.; Li, Z.L.; Ba, R.J. Recent catastrophic debris flows in Luding county, SW China: Geological hazards, rainfall analysis and dynamic characteristics. Nat. Hazards 2016, 55, 523–542. [Google Scholar] [CrossRef]
  5. Tong, L.; Qi, W.; An, G.; Liu, C. Remote sensing survey of major geological disasters in the Himalayas. J. Eng. Geol. 2019, 27, 496. [Google Scholar]
  6. Rai, P.K.; Mohan, K.; Kumra, V.K. Landslide hazard and its mapping using remote sensing and GIS. J. Sci. Res. 2014, 58, 1–13. [Google Scholar]
  7. Zhang, Q.; Liang, Z.; Liu, W.; Peng, W.; Huang, H.; Zhang, S.; Chen, L.; Jiang, K.; Liu, L. Landslide Susceptibility Prediction: Improving the Quality of Landslide Samples by Isolation Forests. Sustainability 2022, 14, 16692. [Google Scholar] [CrossRef]
  8. Corominas, J.; van Westen, C.; Frattini, P.; Cascini, L.; Malet, J.-P.; Fotopoulou, S.; Catani, F.; Van Den Eeckhaut, M.; Mavrouli, O.; Agliardi, F.; et al. Recommendations for the quantitative analysis of landslide risk. Bull. Eng. Geol. Environ. 2014, 73, 209–263. [Google Scholar] [CrossRef] [Green Version]
  9. Reichenbach, P.; Rossi, M.; Malamud, B.D.; Mihir, M.; Guzzetti, F. A Review of Statistically-Based Landslide Susceptibility Models. Earth Sci. Rev. 2018, 180, 60–91. [Google Scholar] [CrossRef]
  10. Shortliffe, E.H.; Buchanan, B.G. A model of inexact reasoning in medicine. Math. Biosci. 1975, 23, 351–379. [Google Scholar] [CrossRef]
  11. Korup, O.; Stolle, A. Landslide prediction from machine learning. Geol. Today 2014, 30, 26–33. [Google Scholar] [CrossRef]
  12. Chang, K.-T.; Merghadi, A.; Yunus, A.P.; Pham, B.T.; Dou, J. Evaluating scale effects of topographic variables in landslide susceptibility models using GIS-based machine learning techniques. Sci. Rep. 2019, 9, 12296. [Google Scholar] [CrossRef] [Green Version]
  13. Liang, Z.; Wang, C.; Zhang, Z.-M.; Khan, K.-U.-J. A comparison of statistical and machine learning methods for debris flow susceptibility mapping. Stoch. Environ. Res. Risk Assess. 2020, 34, 1887–1907. [Google Scholar] [CrossRef]
  14. 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] [Green Version]
  15. Liang, Z.; Wang, C.; Khan, K.U.J. Application and comparison of different ensemble learning machines combining with a novel sampling strategy for shallow landslide susceptibility mapping. Stoch. Environ. Res. Risk Assess. 2021, 35, 1243–1256. [Google Scholar] [CrossRef]
  16. Merghadi, A.; Yunus, A.P.; Dou, J.; Whiteley, J.; ThaiPham, B.; Bui, D.T.; Avtar, R.; Abderrahmane, B. Machine learning methods for landslide susceptibility studies: A comparative overview of algorithm performance. Earth Sci. Rev. 2020, 207, 103225. [Google Scholar] [CrossRef]
  17. Youssef, A.M. Landslide susceptibility delineation in the Ar-Rayth area, Jizan, Kingdom of Saudi Arabia, using analytical hierarchy process, frequency ratio, and logistic regression models. Environ. Earth Sci. 2015, 73, 8499–8518. [Google Scholar] [CrossRef]
  18. Li, L.; Lan, H.; Guo, C.; Zhang, Y.; Li, Q.; Wu, Y. A modified frequency ratio method for landslide susceptibility assessment. Landslides 2017, 14, 727–741. [Google Scholar] [CrossRef]
  19. Chen, Z.; Liang, S.; Ke, Y.; Yang, Z.; Zhao, H. Landslide susceptibility assessment using evidential belief function, certainty factor and frequency ratio model at Baxie River basin, NW China. Geocarto Int. 2019, 34, 348–367. [Google Scholar] [CrossRef]
  20. Gu, J.; Wang, Z.; Kuen, J.; Ma, L.; Shahroudy, A.; Shuai, B.; Liu, T.; Wang, X.; Wang, G.; Cai, J.; et al. Recent Advances in Convolutional Neural Networks. Pattern Recognit. 2018, 77, 354–377. [Google Scholar] [CrossRef] [Green Version]
  21. Ji, S.; Yu, D.; Shen, C.; Li, W.; Xu, Q. Landslide detection from an open satellite imagery and digital elevation model dataset using attention boosted convolutional neural networks. Landslides 2020, 17, 1337–1352. [Google Scholar] [CrossRef]
  22. Park, Y.; Yang, H.S. Convolutional neural network based on an extreme learning machine for image classification. Neurocomputing 2019, 339, 66–76. [Google Scholar] [CrossRef]
  23. Sameen, M.I.; Pradhan, B.; Lee, S. Application of convolutional neural networks featuring Bayesian optimization for landslide susceptibility assessment. Catena 2020, 186, 104249. [Google Scholar] [CrossRef]
  24. Fang, Z.; Wang, Y.; Peng, L.; Hong, H. Integration of convolutional neural network and conventional machine learning classifiers for landslide susceptibility mapping. Comput. Geosci. 2020, 139, 104470. [Google Scholar] [CrossRef]
  25. Varnes, D.J. Landslide Hazard Zonation: A Review of Principles and Practice; United Nations: New York, NY, USA, 1984.
  26. Kornejady, A.; Ownegh, M.; Bahremand, A. Landslide susceptibility assessment using maximum entropy model with two different data sampling methods. Catena 2017, 152, 144–162. [Google Scholar] [CrossRef]
  27. Van Westen, C.J.; Castellanos, E.; Kuriakose, S.L. Spatial data for landslide susceptibility, hazard, and vulnerability assessment: An overview. Eng. Geol. 2008, 102, 112–131. [Google Scholar] [CrossRef]
  28. Hong, H.; Pradhan, B.; Xu, C.; Bui, D.T. Spatial prediction of landslide hazard at the Yihuang area (China) using two-class kernel logistic regression, alternating decision tree and support vector machines. Catena 2015, 133, 266–281. [Google Scholar] [CrossRef]
  29. Sameen, M.I.; Pradhan, B.; Bui, D.T.; Alamri, A.M. Systematic sample subdividing strategy for training landslide susceptibility models. Catena 2019, 187, 104358. [Google Scholar] [CrossRef]
  30. Magliulo, P.; DiLisio, A.; Russo, F.; Zelano, A. Geomorphology and landslide susceptibility assessment using GIS and bivariate statistics: A case study in southern Italy. Nat. Hazards 2008, 47, 411–435. [Google Scholar] [CrossRef]
  31. Zevenbergen, L.W.; Thorne, C.R. Quantitative analysis of land surface topography. Earth Surf. Process. Landf. 1987, 12, 47–56. [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. Evans, I.S. An integrated system of terrain analysis and slope mapping. Z. Geomorphol. Suppl. Stuttg. 1980, 36, 274–295. [Google Scholar]
  34. Liang, Z.; Wang, C.; Duan, Z.; Liu, H.; Liu, X.; Jan Khan, K.U. A hybrid model consisting of supervised and unsupervised learning for landslide susceptibility mapping. Remote Sens. 2021, 13, 1464. [Google Scholar] [CrossRef]
  35. Carrara, A.; Cardinali, M.; Detti, R.; Guzzetti, F.; Pasqui, V.; Reichenbach, P. GIS techniques and statistical models in evaluating landslide hazard. Earth Surf. Process. Landf. 1991, 16, 427–445. [Google Scholar] [CrossRef]
  36. Carrara, A.; Cardinali, M.; Guzzetti, F.; Reichenbach, P. GIS technology in mapping landslide hazard. In Geographical Information Systems in Assessing Natural Hazards; Springer: Berlin/Heidelberg, Germany, 1995; pp. 135–175. [Google Scholar]
  37. Carrara, A.; Guzzetti, F. (Eds.) Geographical Information Systems in Assessing Natural Hazards; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013; Volume 5. [Google Scholar]
  38. Yin, K.L. Statistical prediction model for slope instability of metamorphosed rocks. In Proceedings of the 5th International Symposium on Landslides, Lausanne, Switzerland, 10–15 July 1988; Volume 2, pp. 1269–1272. [Google Scholar]
  39. Chen, W.; Yang, Z. Landslide susceptibility modeling using bivariate statistical-based logistic regression, naïve Bayes, and alternating decision tree models. Bull. Eng. Geol. Environ. 2023, 82, 190. [Google Scholar] [CrossRef]
  40. Ganga, A.; Elia, M.; D’Ambrosio, E.; Tripaldi, S.; Capra, G.F.; Gentile, F.; Sanesi, G. Assessing landslide susceptibility by coupling spatial data analysis and logistic model. Sustainability 2022, 14, 8426. [Google Scholar] [CrossRef]
  41. Ye, P.; Yu, B.; Chen, W.; Liu, K.; Ye, L. Rainfall-induced landslide susceptibility mapping using machine learning algorithms and comparison of their performance in Hilly area of Fujian Province, China. Nat. Hazards 2022, 113, 965–995. [Google Scholar] [CrossRef]
  42. Sahin, E.K. Comparative analysis of gradient boosting algorithms for landslide susceptibility mapping. Geocarto Int. 2022, 37, 2441–2465. [Google Scholar] [CrossRef]
  43. Krizhevsky, A.; Sutskever, I.; Hinton, G.E. ImageNet Classification with Deep Convolutional Neural Networks. Commun. ACM 2017, 60, 84–90. [Google Scholar] [CrossRef] [Green Version]
  44. Renza, D.; Cárdenas, E.A.; Martinez, E.; Weber, S.S. CNN-Based Model for Landslide Susceptibility Assessment from Multispectral Data. Appl. Sci. 2022, 12, 8483. [Google Scholar] [CrossRef]
  45. Saha, S.; Saha, A.; Hembram, T.K.; Mandal, K.; Sarkar, R.; Bhardwaj, D. Prediction of spatial landslide susceptibility applying the novel ensembles of CNN, GLM and random forest in the Indian Himalayan region. Stoch. Environ. Res. Risk Assess. 2022, 36, 3597–3616. [Google Scholar] [CrossRef]
  46. Liang, Z.; Liu, W.; Peng, W.; Chen, L.; Wang, C. Improved shallow landslide susceptibility prediction based on statistics and ensemble learning. Sustainability 2022, 14, 6110. [Google Scholar] [CrossRef]
  47. Lv, L.; Chen, T.; Dou, J.; Plaza, A. A hybrid ensemble-based deep-learning framework for landslide susceptibility mapping. Int. J. Appl. Earth Obs. Geoinf. 2022, 108, 102713. [Google Scholar] [CrossRef]
  48. Ling, X.; Zhu, Y.; Ming, D.; Chen, Y.; Zhang, L.; Du, T. Feature Engineering of Geohazard Susceptibility Analysis Based on the Random Forest Algorithm: Taking Tianshui City, Gansu Province, as an Example. Remote Sens. 2022, 14, 5658. [Google Scholar] [CrossRef]
  49. Huang, F.; Tao, S.; Chang, Z.; Huang, J.; Fan, X.; Jiang, S.-H.; Li, W. Efficient and automatic extraction of slope units based on multi-scale segmentation method for landslide assessments. Landslides 2021, 18, 3715–3731. [Google Scholar] [CrossRef]
  50. Chang, Z.; Catani, F.; Huang, F.; Liu, G.; Meena, S.R.; Huang, J.; Zhou, C. Landslide susceptibility prediction using slope unit-based machine learning models considering the heterogeneity of conditioning factors. J. Rock Mech. Geotech. Eng. 2023, 15, 1127–1143. [Google Scholar] [CrossRef]
  51. Huang, F.; Pan, L.; Fan, X.; Jiang, S.H.; Huang, J.; Zhou, C. The uncertainty of landslide susceptibility prediction modeling: Suitability of linear conditioning factors. Bull. Eng. Geol. Environ. 2022, 81, 182. [Google Scholar] [CrossRef]
  52. Bravo-López, E.; Del Castillo, T.F.; Sellers, C.; Delgado-García, J. Landslide susceptibility mapping of landslides with artificial neural networks: Multi-approach analysis of back propagation algorithm applying the neuralnet package in Cuenca, Ecuador. Remote Sens. 2022, 14, 3495. [Google Scholar] [CrossRef]
  53. Chen, W.; Li, Y. GIS-based evaluation of landslide susceptibility using hybrid computational intelligence models. Catena 2020, 195, 104777. [Google Scholar] [CrossRef]
  54. Dou, J.; Yunus, A.P.; Bui, D.T.; Merghadi, A.; Sahana, M.; Zhu, Z.; Chen, C.W.; Han, Z.; Pham, B.T. Improved landslide assessment using support vector machine with bagging, boosting, and stacking ensemble machine learning framework in a mountainous watershed, Japan. Landslides 2020, 17, 641–658. [Google Scholar] [CrossRef]
  55. Pourghasemi, H.R.; Kariminejad, N.; Amiri, M.; Edalat, M.; Zarafshar, M.; Blaschke, T.; Cerda, A. Assessing and mapping multi-hazard risk susceptibility using a machine learning technique. Sci. Rep. 2020, 10, 3203. [Google Scholar] [CrossRef] [Green Version]
  56. Murillo-García, F.G.; Steger, S.; Alcántara-Ayala, I. Landslide susceptibility: A statistically-based assessment on a depositional pyroclastic ramp. J. Mt. Sci. 2019, 16, 561–580. [Google Scholar] [CrossRef]
  57. Zhang, Y.-X.; Lan, H.-X.; Li, L.-P.; Wu, Y.-M.; Chen, J.-H.; Tian, N.-M. Optimizing the frequency ratio method for landslide susceptibility assessment: A case study of the Caiyuan Basin in the southeast mountainous area of China. J. Mt. Sci. 2020, 17, 340–357. [Google Scholar] [CrossRef]
  58. Saygin, F.; Şişman, Y.; Dengiz, O.; Şişman, A. Spatial assessment of landslide susceptibility mapping generated by fuzzy-AHP and decision tree approaches. Adv. Space Res. 2023, 71, 5218–5235. [Google Scholar] [CrossRef]
  59. Liang, Z.; Wang, C.; Han, S.; Khan, K.U.J.; Liu, Y. Classification and susceptibility assessment of debris flow based on a semiquantitative method combination of the fuzzy C-means algorithm, factor analysis and efficacy coefficient. Nat. Hazards Earth Syst. Sci. 2020, 20, 1287–1304. [Google Scholar] [CrossRef]
  60. Pourghasemi, H.R.; Moradi, H.R.; Fatemi Aghda, S.M. Landslide susceptibility mapping by binary logistic regression, analytical hierarchy process, and statistical index models and assessment of their performances. Nat. Hazards 2013, 69, 749–779. [Google Scholar] [CrossRef]
  61. Rai, D.K.; Xiong, D.; Zhao, W.; Zhao, D.; Zhang, B.; Dahal, N.M.; Wu, Y.; Baig, M.A. An Investigation of Landslide Susceptibility Using Logistic Regression and Statistical Index Methods in Dailekh District, Nepal. Chin. Geogr. Sci. 2022, 32, 834–851. [Google Scholar] [CrossRef]
  62. Bukhari, M.H.; da Silva, P.F.; Pilz, J.; Istanbulluoglu, E.; Görüm, T.; Lee, J.; Karamehic-Muratovic, A.; Urmi, T.; Soltani, A.; Wilopo, W.; et al. Community perceptions of landslide risk and susceptibility: A multi-country study. Landslides 2023, 20, 1321–1334. [Google Scholar] [CrossRef]
  63. Sajadi, P.; Sang, Y.F.; Gholamnia, M.; Bonafoni, S.; Mukherjee, S. Evaluation of the landslide susceptibility and its spatial difference in the whole Qinghai-Tibetan Plateau region by five learning algorithms. Geosci. Lett. 2022, 9, 9. [Google Scholar] [CrossRef]
Figure 1. Location of the study area showing elevation and landslide samples.
Figure 1. Location of the study area showing elevation and landslide samples.
Applsci 13 07276 g001
Figure 2. The Sikkim earthquake triggered landslides in Yadong County: (a) S204 road was interrupted; (b) houses collapsed.
Figure 2. The Sikkim earthquake triggered landslides in Yadong County: (a) S204 road was interrupted; (b) houses collapsed.
Applsci 13 07276 g002
Figure 3. Remote sensing image interpretation: (a) debris flow in Kangbu township; (b) debris flow in Duina township.
Figure 3. Remote sensing image interpretation: (a) debris flow in Kangbu township; (b) debris flow in Duina township.
Applsci 13 07276 g003
Figure 4. Study area thematic maps: (a) elevation; (b) slope; (c) MED; (d) plan curvature; (e) profile curvature; (f) TWI; (g) DTF; (h) DTR; (i) DTS; (j) rainfall; (k) NDVI.
Figure 4. Study area thematic maps: (a) elevation; (b) slope; (c) MED; (d) plan curvature; (e) profile curvature; (f) TWI; (g) DTF; (h) DTR; (i) DTS; (j) rainfall; (k) NDVI.
Applsci 13 07276 g004aApplsci 13 07276 g004b
Figure 5. Flow chart of this study.
Figure 5. Flow chart of this study.
Applsci 13 07276 g005
Figure 6. Analysis of ROC curve for the landslide susceptibility map: (a) success rate curve of landslide using the training data set; (b) prediction rate curve of landslide using the validation dataset.
Figure 6. Analysis of ROC curve for the landslide susceptibility map: (a) success rate curve of landslide using the training data set; (b) prediction rate curve of landslide using the validation dataset.
Applsci 13 07276 g006
Figure 7. Parametric importance graphics designed by Gini index.
Figure 7. Parametric importance graphics designed by Gini index.
Applsci 13 07276 g007
Figure 8. Landslide susceptibility maps: (a) IV; (b) LR model; (c) CatBoost model; (d) CNN model.
Figure 8. Landslide susceptibility maps: (a) IV; (b) LR model; (c) CatBoost model; (d) CNN model.
Applsci 13 07276 g008
Figure 9. Percentages of areas in different susceptibility classes for landslide.
Figure 9. Percentages of areas in different susceptibility classes for landslide.
Applsci 13 07276 g009
Table 1. Confusion matrix analysis with evaluation measures.
Table 1. Confusion matrix analysis with evaluation measures.
Actual ValuesAccuracySensitivitySpecificity
Positive (1)Negative (0)(TP + TN)/(TP + TN + FP + FN)TN/(TN + FP)TN/(TN + FP)
Predicted
Values
Positive (1)True Positive (TP)False Negative (FN)
Negative (0)False Positive (FP)True Negative (TN)
Table 2. Multicollinearity diagnosis indexes for variables.
Table 2. Multicollinearity diagnosis indexes for variables.
VariablesVIF
Elevation5.117
Slope3.426
MED5.726
Plan curvature1.499
Profile curvature1.291
TWI6.071
Distance to fault2.641
Distance to stream4.492
Distance to road4.302
Average annual precipitation1.763
NDVI2.697
Table 3. Models’ performance using related indices.
Table 3. Models’ performance using related indices.
TrainingValidation
ParameterLRCatBootsCNNLRCatBootsCNN
Sensitivity (%)81.4582.9684.2179.3876.2879.38
Specificity (%)84.7990.2793.5276.0085.0091.00
Accuracy (%)83.1386.6388.8877.6680.7185.28
AUC0.8970.9300.9440.8380.8930.908
Table 4. Landslide conditioning factors in this study.
Table 4. Landslide conditioning factors in this study.
Conditioning FactorsZoneNi/NSi/SIV
Elevation (m)<27001.21%1.09%0.100
2700~370021.98%8.26%0.978
3700~470043.80%49.88%−0.131
4700~570028.63%39.29%−0.317
>57004.44%1.48%1.098
Slope (°)<1011.90%35.95%−1.106
10~2025.60%29.88%−0.154
20~3038.31%25.67%0.400
30~4022.18%8.13%1.003
>402.02%0.36%1.709
MED (m)<1007.46%36.41%−1.585
100~20013.91%21.80%−0.449
200~30011.90%15.82%−0.285
300~40020.56%10.76%0.648
>40046.17%15.21%1.110
Plan curvature<−0.40.20%0.12%0.538
−0.4~−0.14.44%4.89%−0.097
−0.1~0.292.74%93.11%−0.004
>0.22.62%1.89%0.328
Profile curvature<−0.30.20%0.59%−1.078
−0.3~031.20%36.57%−0.157
0~0.367.54%62.50%0.077
>0.31.01%0.33%1.106
TWI<67.26%5.85%0.215
6~746.57%32.13%0.371
7~833.27%27.54%0.189
8~98.87%14.48%−0.490
9~102.42%9.11%−1.326
>101.61%10.89%−1.910
Distance to faults (km)<431.05%20.27%0.426
4~816.33%14.71%0.105
8~127.46%14.57%−0.670
12~1622.18%16.00%0.327
16~2014.31%14.99%−0.046
20~248.67%12.22%−0.343
24~280.60%6.10%−2.311
>280.20%1.14%−1.734
Distance to streams (km)<147.58%14.51%1.187
1~319.35%25.62%−0.280
3~57.06%19.97%−1.040
5~76.65%14.66%−0.790
7~99.07%11.07%−0.199
9~115.44%7.00%−0.252
11~133.63%4.38%−0.189
>131.61%2.78%−0.544
Distance to roads (km)<357.66%32.20%0.583
3~611.69%25.23%−0.769
6~914.11%18.23%−0.256
9~126.20%9.08%−0.374
12~155.44%6.37%−0.158
15~181.41%4.46%−1.151
18~212.62%2.07%0.237
>210.81%2.35%−1.070
Average annual precipitation (mm)<48032.26%50.99%−0.458
480~58029.44%24.55%0.182
580~68024.60%11.55%0.756
680~7803.02%9.38%−1.132
>78010.69%3.53%1.107
NDVI<0.1516.53%12.68%0.265
0.15~0.325.20%36.62%−0.374
0.3~0.4515.12%19.74%−0.266
0.45~0.615.93%17.94%−0.119
0.6~0.7522.18%9.93%0.804
>0.755.04%3.09%0.490
Table 5. Conditioning factors assigned by the Gini index.
Table 5. Conditioning factors assigned by the Gini index.
MethodDTFMEDDTSElevationPlan CurvatureSlope
Gini index5.365.004.804.103.692.80
Table 6. The IV of different landslide susceptibility levels.
Table 6. The IV of different landslide susceptibility levels.
ModelClassPercentage of Area (%)Percentage of Landslide Area (%)IV
IVVery low32.05%8.06%−1.74
Low8.64%7.66%−1.19
Moderate7.96%5.04%−0.50
High14.08%20.36%−0.14
Very high37.27%58.87%0.70
LRVery low40.17%5.24%−1.60
Low12.24%4.44%−0.47
Moderate5.54%6.05%−0.10
High17.75%21.17%0.14
Very high24.30%63.10%0.88
CatBoostVery low30.16%5.24%−1.74
Low14.42%4.44%−1.18
Moderate7.82%6.05%−0.25
High20.28%21.17%0.04
Very high27.31%63.10%0.84
CNNVery low28.99%4.44%−1.88
Low19.65%8.67%−0.82
Moderate14.21%12.30%−0.14
High9.19%11.09%0.19
Very high27.96%63.51%0.69
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liang, Z.; Peng, W.; Liu, W.; Huang, H.; Huang, J.; Lou, K.; Liu, G.; Jiang, K. Exploration and Comparison of the Effect of Conventional and Advanced Modeling Algorithms on Landslide Susceptibility Prediction: A Case Study from Yadong Country, Tibet. Appl. Sci. 2023, 13, 7276. https://doi.org/10.3390/app13127276

AMA Style

Liang Z, Peng W, Liu W, Huang H, Huang J, Lou K, Liu G, Jiang K. Exploration and Comparison of the Effect of Conventional and Advanced Modeling Algorithms on Landslide Susceptibility Prediction: A Case Study from Yadong Country, Tibet. Applied Sciences. 2023; 13(12):7276. https://doi.org/10.3390/app13127276

Chicago/Turabian Style

Liang, Zhu, Weiping Peng, Wei Liu, Houzan Huang, Jiaming Huang, Kangming Lou, Guochao Liu, and Kaihua Jiang. 2023. "Exploration and Comparison of the Effect of Conventional and Advanced Modeling Algorithms on Landslide Susceptibility Prediction: A Case Study from Yadong Country, Tibet" Applied Sciences 13, no. 12: 7276. https://doi.org/10.3390/app13127276

APA Style

Liang, Z., Peng, W., Liu, W., Huang, H., Huang, J., Lou, K., Liu, G., & Jiang, K. (2023). Exploration and Comparison of the Effect of Conventional and Advanced Modeling Algorithms on Landslide Susceptibility Prediction: A Case Study from Yadong Country, Tibet. Applied Sciences, 13(12), 7276. https://doi.org/10.3390/app13127276

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