Abstract
Landslides are one of the most widespread disasters and threaten people’s lives and properties in many areas worldwide. Landslide susceptibility mapping (LSM) plays a crucial role in the evaluation and extenuation of risk. To date, a large number of machine learning approaches have been applied to LSM. Of late, a high-level convolutional neural network (CNN) has been applied with the intention of raising the forecast precision of LSM. The primary contribution of the research was to present a model which was based on the CNN for LSM and methodically compare its capability with the traditional machine learning approaches, namely, support vector machine (SVM), logistic regression (LR), and random forest (RF). Subsequently, we used this model in the Wenchuan region, where a catastrophic earthquake happened on 12 May 2008 in China. There were 405 valuable landslides in the landslide inventory, which were divided into a training set (283 landslides) and validation set (122 landslides). Furthermore, 11 landslide causative factors were selected as the model’s input, and each model’s output was reclassified into five intervals according to the sensitivity. We also evaluated the model’s performance by the receiver operating characteristic (ROC) curve and several statistical metrics, such as precision, recall, F1-score, and other measures. The results indicated that the CNN-based methods achieved the best performance, with the success-rate curve (SRC) and prediction-rate curve (PRC) approaches reaching 93.14% and 91.81%, respectively. The current research indicated that the approach based on the CNN for LSM had both outstanding goodness-of-fit and excellent prediction capability. Generally, the LSM in our research is capable of advancing the ability to assess landslide susceptibility.
1 Introduction
Landslides are one of the main disasters which exist in mountainous zones worldwide, threatening people, properties, and the natural environment (Qiang and Xiu-June 2011; Dou et al., 2019; Mengistu et al., 2019; Shano et al., 2020). These events represent one of the most severe geological disasters and are caused by landforms, geology, hydrology, and human activities (Jaafari et al., 2019), possibly causing massive casualties and damaging many infrastructures when large-scale landslides occur. Especially in Wenchuan China, from the magnitude 8 earthquake in 2008 (Chen M. et al., 2020) to now, this region maintains a relatively high activity level, seriously threatening the rebuilding of the county and the local people (Chong et al., 2013; Zhou et al., 2021). Hence, reliability prediction and timely management of landslides are of vital importance for preventing and reducing the cost caused by landslides (Wu et al., 2015).
LSM describes the space division of landslides in a specific zone on the foundation of past and present landslide incidents. In LSM, the most important step in determining the evaluation accuracy involves defining which modeling approach is to be used (Akinci and Zeybek, 2021). In recent years, people have applied qualitative methods and quantitative methods (Depicker et al., 2020) in LSM. These qualitative methods are based on a knowledge-driven approach. The statistical models used in LSM include the analytical hierarchy process (AHP) (Hepdeniz, 2020), weight of evidence (Saha and Saha, 2020), frequency ratio (Zhang et al., 2020), information gain (Chen T. et al., 2020), etc. However, these methods have been gradually replaced by machine learning methods, which have achieved better performance in LSM, and a data-driven model is regarded as a more useful landslide spatiotemporal mapping tool in the quantitative method, including logistic regression (LR) (Akgun, 2012; Pourghasemi et al., 2013), support vector machine (SVM) (Pham et al., 2019), naïve Bayes (Xie M. et al., 2021), decision tree (DT) (Nhu et al., 2020a), and the artificial neural network (ANN) (Nhu et al., 2020b). While some progresses have been made in LSM, the performance of the methods in LSM still needs to be improved. For example, for some general methods such as DT and the ANN, their simple model structure makes it hard to fully capture the complex linear relationship or nonlinear relations between landslides and their surroundings (Kong et al., 2021).
Apparently, neither knowledge-driven approaches nor traditional machine learning approaches such as the SVM can make full use of the environmental information around landslides. More recently, except for the methods given above, convolutional neural networks (CNNs) have been successfully used in LSM (Wang et al., 2019; Fang et al., 2020; Pham et al., 2020; Wei et al., 2021). What is worth mentioning is that Wang et al., (2019) thoroughly compared one-dimensional, two-dimensional, and three-dimensional CNNs according to information from Yanshan County. A multiresolution sampling method (Yi et al., 2020) that can further take advantage of the surroundings of the landslides was proposed. Despite these developments, there is no consensus on the most appropriate method for LSM in the literature. Although it is novel to use CNN to analyze landslide susceptibility compared with other approaches, previous studies have mainly focused on the comparison between machine learning and statistical methods (Goetz et al., 2015; Pham et al., 2016; Merghadi et al., 2020), with few attempts to compare these approaches with conventional machine learning. To fill this gap, our research aims to compare the CNN with conventional machine learning approaches (LR, RF, and SVM), selecting 11 related factors such as terrain, environmental conditions, and the engineering behaviors of human beings to construct a landslide susceptibility assessment mechanism. We adopted the ROC curve, precision, recall, F1-score, and MCC coefficient to carry out assessment and verify our models, facilitating the prediction of landslide susceptibility.
As for this research, the mapping of landslide susceptibility in the Wenchuan region involves four steps:
1. For assembling the data required for the research, a causative factor determination procedure was performed by using multicollinearity assessment and Pearson’s analysis to choose appropriate factors, predisposing factors to establish datasets for CNN-based models.
2. Modeling landslide susceptibility in terms of CNN, logistic regression, random forest, and support vector machine.
3. Generating the landslide susceptibility mapping of the study area.
4. Evaluating the goodness-of-fit and prediction capability in training datasets and validation datasets.
2 Description of the Study Area
As illustrated in Figure 1, the study area is on the northwest edge of the Sichuan Basin, located in Sichuan Province, Southwest China. This area covers an area of 4,084 square kilometers and spans 102° 39′ E and 103° 44′ E longitudes and 30° 45′ N and 31° 43′ N latitudes. The pie chart in Figure 1 shows the percentage of landslides that occurred in the area for each year from 2013 to 2018. In this region, the variation in rock characteristics is complex and dominated by soft rock and soil layers, with a fourth-century loose soil layer and strong decaying magmatic layer as its symbols. The terrain in Wenchuan is mostly dominated by two mountain ranges: the Longmen mountain system and Qionglai mountain system. This area also contains the Maowen fault and Beichuan–Yinxiu fault (Xie W. et al., 2021). With prosperous water systems and highly fractured mounds, the main water systems in Wenchuan County are the Minjiang River, Zagunao River, and Yuzixi River.
FIGURE 1
The magnitude 8 earthquake that occurred on 12 May 2008 caused a large-scale loosened rock mass and loose sediments (Cui et al., 2021). Under the action of rainfall and other natural factors, it leads to the saturation of the soil and rock layer on the slope and even water accumulation in the water barrier at the lower part of the slope, which increases the weight of the landslide and reduces the shear strength of the soil and rock layer, so the frequency of landslide increases. In this region, the main forms of landslide collapse, landslides, and debris flow severely disrupt the local traffic and economy. Thus, analyzing the factors of landslides and evaluating landslide susceptibility will contribute greatly to distinguish and figure out where the landslide-prone areas are.
3 Data Preparation
3.1 Landslide Inventory
The basis of landslide susceptibility studies is landslide inventory. Generally, a complete landslide inventory should keep track of the sites, size, outline, and time of the landslide, which can be developed into a basic database for further study. The data of this research mainly originated from the detailed landslide data of Wenchuan County at the GDEMV2 30 m resolution digital elevation. In total, 405 landslides were selected from the detailed geological hazard survey data of Wenchuan (2013F02D2018). Taking the image resolution and aims into consideration, we converted all the data into 30 m × 30 m raster data. The landslide inventory was divided into two parts at random for the intention of evaluating the landslide susceptibility: 70% of landslides were utilized for the purpose of model training, and the remaining 30% were utilized for the purpose of validation. In addition to this, the same quantity of non-landslide sites was chosen at random outside the buffer area to perform the training and validation. Taking the scale of landslides into consideration, this research sets the buffer distance as 100 m.
3.2 Landslide Causative Factors
At present, there are no general rules existing for choosing the landslide causative factors (Yi et al., 2020). It is generally believed that not a single factor can contribute to the appearance of landslides, requiring a close combination of many environmental factors such as terrain and geology. Therefore, choosing appropriate landslide factors can better explain the occurrence of landslides. Because of the complexity and diversity in the trigger source of landslides, we chose the factors according to the variation in the specific region (Van Westen et al., 2008). Considering the conditions of the nature and availability of data in the research area, 11 causative factors were selected: elevation, slope, aspect, plan curvature, profile curvature, relief amplitude, distance to rivers, distance to roads, distance to faults, NDVI, and lithology (Figure 2). The data that are utilized in this research are mainly derived from Wenchuan landslide survey data (2013–2018), as shown in Table 1.
FIGURE 2
TABLE 1
| No | Factors | Data Sources and Scale |
|---|---|---|
| 1 | Elevation | Digital elevation model (DEM) data of 90 m |
| 2 | Slope | |
| 3 | Aspect | |
| 4 | Relief amplitude | |
| 5 | Plan curvature | |
| 6 | Profile curvature | |
| 7 | NDVI | Landsat 8 OLI image |
| 8 | Distance to faults | |
| 9 | Lithology | National Basic Geographic Database 1: 250,000 |
| 10 | Distance to rivers | |
| 11 | Distance to roads |
Data sources of landslide-impacting elements.
Among these influencing factors, elevation occupies a decisive role in landslides. As for the study area, landslides are concentrated between 875 and 2613 m. The aspect in the research zone is associated with the location of the landslides as slope deposits or quaternary cover are usually distinct in different directions (Nandi and Shakoor, 2010). The slope reflects the sharpness of the research zone, and generally, steep slopes are more prone to landslides than gentle slopes (Dai et al., 2001). Plan and profile curvature is efficient to mirror the complication of the terrain (Oh and Pradhan, 2011). In mountain areas, the normalized vegetation index (NDVI) is an indispensable index of slope stability, and we classified the NDVI into six ranges: 0, 0–0.08, 0.08–0.13, 0.13–0.18, 0.18–0.24, and 0.24. Relief amplitude is efficient in mirroring the gravitational underlying power of terrain, and it is often taken into consideration in the research studies related to the landslide. Geological factors are another key aspect causing landslides to occur; we divided the geological factors into two categories: engineering rock groups and fault distance groups. Engineering rock groups are one of the causative factors of LSM, and some geological formations will be better off for landslide generation (Wang et al., 2019). In the study area, engineering rock groups are divided into three groups: soft and hard interbedded rock groups (A), soft rock groups (B), and hard rock groups (C). Distances to faults were obtained through the buffer classification of ArcGIS 10.8 software, which is the key factor influencing the size and subsequent spread of landslides (Fang et al., 2020), and generally, landslides occur mainly near active faults. It is particularly significant to notice that landslides in the study area have occurred mainly near active faults since the large earthquake that occurred in 2008. Human engineering activities (distance to roads) significantly impact the scale and expansion of landslides, and the fluctuations in river flow significantly affect the surrounding side slopes and slopes.
4 Methodology
4.1 Convolutional Neural Network
With the landslide inventory and causative factors, the next step is to format these data into a uniform input dataset to construct the CNN model. The layers of the 11 landslide causative factors are superimposed together to form a tensor of size 11 × w × h, where w and h are the width and height of the research area, respectively. The established landslide inventory is then superimposed onto the abovementioned layers to obtain the raster location corresponding to the landslide point, which is expanded into square data of size n × n centered on that raster pixel so that more environmental information around the landslide will be used in model construction. It is to be noted that the scale of the raster influences the CNN model’s performance. The raster used for training should be set up on the resolution of the raster (Yi et al., 2020), and the final raster size used for learning in this study was determined to be 17 × 17 by using the trial-and-error approach.
As a characteristic deep learning algorithm, convolutional neural network (Krizhevsky et al., 2012) is widely used for image classification and recognition (Szegedy et al., 2015; Gu et al., 2018). The classical CNN structure consists of an input layer, multiple convolutional layers, multiple pooling layers, and fully connected layers (He et al., 2016). In this study, the input layer imports a normalized (uniform size) 11-dimensional training tensor, and the neurons extract some basic pixel features, including edges, corner points, etc. Next, the convolutional layer obtains the feature maps by convolutional operations, and at each position, the units that are derived from various feature maps gain their own kinds of features. Usually, a convolutional layer consists of numerous feature maps with various weight vectors in order to retain the richer features. The convolutional layer is later connected to a pooling layer for downsampling operation. For LSM, it can reduce the resolution of the image and the number of parameters and obtains robust translation and deformation. Finally, the chance of incidence of landslides is calculated by the fully connected layer to output the classification result. As for this research, the organization of the designed convolutional neural network is illustrated in Figure 3.
FIGURE 3
4.2 Conventional Machine Learning Methods
1) Logistic regression: It is a comprehensive linear model that is commonly utilized in dichotomous problems and widely used for landslide susceptibility modeling problems (Kavzoglu et al., 2014; Colkesen et al., 2016). The cost function is first set up, after which the optimum model parameters are addressed iteratively by optimization methods. Logistic regression commonly uses the sigmoid functional equation as the cost function.
In the landslide susceptibility partition rating, the probability of landslide occurrence can be formulated aswhere βi is the regression coefficient, xi is the individual influence factor, P is the possibility of the incidence of landslide, and Z is the sum of all the influence factors after simultaneous loading.
2) Random forest: It is an integrated learning approach which is utilized widely in order to perform the regression (Nhu et al., 2020a; Wang et al., 2020). The integrated learning approach produces numerous independent training sets and a quantity of classification and regression trees (CARTs) by combined bagging, and it is expressed aswhere h (X, θk) is the unprocessed categorical regression tree produced by the CART algorithm, x is the input vector, and θk is an independently distributed random vector.
A singular decision tree is a poor classifier because of high-level bias or variance among data. Therefore, random forest creates multiple decision trees and then calculates the classification by aggregating these decision trees. The integrated model is found to perform better than the individual model (Nhu et al., 2020a). Therefore, creating decision trees individually to form a random forest is considered a better approach. In this study, a spatial training set was created by using the bagging method, and CARTs were created by 11 different landslide causative factors. Finally, all the decision trees were combined, and the class with the most votes was used as the result of the model.
3) Support vector machine: The core idea of the support vector machine is to classify samples in higher-dimensional space by converting the initial space to a new higher-dimensional one by the definition of the inner product function. By means of kernel functions, the linear classification boundary can be transformed into a nonlinear classification boundary to obtain better generalization performance (Sevgen et al., 2019; Sun et al., 2020). As for this research, the radial basis kernel function (RBF) (Nhu et al., 2020c) is expressed as follows:in which z represents the center of the kernel function and γ represents the width parameter, which mainly controls the radial range of action of the function. The best parameters are found cyclically by using the GridSearchCV method (Vapnik, 1999) in the sklearn package.
5 Results and Discussion
5.1 Multicollinearity Analysis of Causative Factors
Feature selection is particularly important during spatial prediction of landslides. There are many features to choose, among which there may be features that are not related to landslides, and the features may depend on each other. The more features, the easier it will cause over-fitting and dimensional disasters, and the more complex the model will be. In this study, tolerance (TOL), variance inflation factor (VIF), and Pearson’s correlations (Al-Najjar et al., 2019) are utilized for the purpose of performing multicollinearity analysis.
1) VIF and TOL: To evaluate the correlation between the landslide impact factors, this study calculated the TOL and VIF so as to confirm the multicollinearity between the different impact factors. It is statistically demonstrated that the variance of the parameter estimator of the explanatory variables is expressed aswhere is the decidable coefficient about xj of the auxiliary regression model on the other explanatory variables, and the second factor is called the variance inflation factor and is denoted as VIF, taking 10 as a judgment boundary.
When the VIF is no more than 10, no multicollinearity exists; when 10 ≤ VIF ≤100, strong multicollinearity exists; and when VIF exceeds 100, there is severe multicollinearity. The tolerance level (TOL) is the reciprocal of the VIF and takes a value between 0 and 1. It is usually considered that when the tolerance level is less than 0.1, there is severe multicollinearity. If the VIF exceeds 10 or TOL is smaller than 0.1, the factor should be detached from the landslide forecasting model.
2) Pearson’s coefficient: The correlation of factors can be analyzed using Pearson’s coefficient. The correlation coefficient R represents the correlation between different factors, and its value generally ranges between [−1,1]in which cov (x, y) is the covariance of x and y, and Var [x] and Var [y] are the variances.
Deriving the formula further, x and y are the two landslide causative factors in this study, and the following equation can be derivedin which are their average, and n represents the number of samples.
The correlation increases as the absolute value of the correlation coefficient increases. Typically, when the absolute value of Pearson’s coefficient exceeds 0.8, it indicates that the correlation is obvious. As for this research, Pearson’s coefficient is used for the purpose of evaluating the magnitude of correlation between two causative factors.
3) Results: The final outcomes of multicollinearity analysis for the 11 landslide causative factors are presented in Table 2. It was found that the largest VIF value for lithology was 8.1569, but still less than 10. The VIF and TOL values for the remaining 10 influences were within the normal interval, and the results indicated that all of the above 11 causative factors could be included in the susceptibility modeling.
TABLE 2
| Landslide causative factors | Tolerance (TOL) | VIF |
|---|---|---|
| Elevation | 0.4801 | 2.0828 |
| Aspect | 0.3015 | 3.3158 |
| Slope | 0.1766 | 5.6616 |
| Plan curvature | 0.4199 | 2.3814 |
| Profile curvature | 0.4567 | 2.1896 |
| NDVI | 0.4898 | 2.0416 |
| Lithology | 0.1225 | 8.1569 |
| Relief amplitude | 0.5416 | 1.8462 |
| Distance to rivers | 0.1328 | 7.5288 |
| Distance to roads | 0.2191 | 5.5630 |
| Distance to faults | 0.2792 | 3.5812 |
Multicollinearity analysis of landslide causative factors.
Figure 4 shows the magnitude of correlation between different causative factors. The correlation to the fault distance and topographic relief shows a maximum positive correlation of 0.25, and the maximum negative correlation is NDVI with a plan curvature of −0.29. Typically, when Pearson’s coefficient exceeds 0.8, it indicates high covariance. The results show that there is no covariance problem with the influencing factors. The findings are in line with the validation results of VIF and TOL, and all the causative factors selected are able to be utilized for landslide susceptibility modeling, which ensures the accuracy and precision of the model.
FIGURE 4
5.2 Landslide Susceptibility Mapping
This subsection conducts experiments on the above-introduced CNN model and three machine learning models for landslide susceptibility analysis. Data processing, model building, and training are implemented in PyTorch by using Python programing. For the parameters in the CNN model, this study uses the trial-and-error approach to optimize the main parameters that may be optimized during the training process, as detailed in Table 3.
TABLE 3
| Parameters | Value |
|---|---|
| Conventional kernel size | 3 × 3 |
| Max pooling size | 2 × 2 |
| Loss function | Cross-entropy |
| Activation function | RELU |
| Optimizer | SGD |
| Epoch | 200 |
| Batch size | 64 |
| Learning rate/momentum | 0.0001/0.2 |
Parameter settings of the CNN model.
Finally, the trained CNN model was applied for the purpose of producing a landslide susceptibility map. The LSM is described in Figure 5. The susceptibility was treated with reclassification into five levels using the Jenks natural fracture approach with ArcGIS 10.8 software, namely, very low, low, moderate, high, and very high. The LSM of the CNN model assumed that the susceptibility classes of the above-occupied areas of 14%, 9%, 6%, 10%, and 61% of the research zone, respectively (Figure 6). The very high susceptibility zones are mainly located in the fault concentration area (southeast area of Wenchuan County) and the water system around the region, accounting for 9.4% of the total area of the research zone (Figure 6). In addition, it can be noticed that the majority of the landslides are situated in the susceptibility zones which are very high, and no landslide occurs in the susceptibility zones which are low, pointing out that the model which is based on the CNN approach can precisely predict the majority of the landslides that are witnessed in the high-susceptibility area.
FIGURE 5
FIGURE 6
The landslide susceptibility maps are produced by using three machine learning models which are LR, RF, and SVM to perform better comparison (Figure 7). Visually, the very high susceptibility areas for all three methods are similar to the CNN model distribution pattern, mainly concentrated in fault concentration areas and around water systems. In terms of the susceptibility which is very high, the outcomes of the CNN and the three machine learning algorithms are basically the same, but the moderate susceptibility based on the LR method is relatively scattered. Notably, the CNN-based method generates significantly lower regions of low susceptibility than the other three machine learning methods.
FIGURE 7
5.3 Validation and Comparison of the Results
The final step in LSM research is to validate the results of the model (Akinci and Zeybek, 2021). Chung and Fabbri, (2003) states that a central point in measuring the importance of a prediction is that it allows for interpretation of the results. If the validation is not carried out, it is impossible to obtain the explanation and the support for the approach (Chung and Fabbri, 2003; He and Kusiak, 2017). Numerous precision evaluation methods were used for the purpose of assessing the property of the CNN model (Saha et al., 2021). As shown in Figure 8, the property of the CNN model was determined by introducing five metrics such as ROC curve and precision, which were based on the foundation of confusion matrices with true positive (TP), false positive (FP), true negative (TN), and false negative (FN).
FIGURE 8
As far as the ROC curves are concerned, the area under the curve (AUC) (Jiao et al., 2019) is the most significant metric during the course of assessing the precision of the quality of the susceptibility model. As previously reported, AUC values between 0.5 and 0.6 show weak performance, 0.6–0.7 indicate average performance, 0.7–0.8 indicate good performance, 0.8–0.9 indicate very good performance, and 0.9–1.0 indicate excellent performance (Cao et al., 2020). As shown in Figure 9, numerically, the CNN structure outperforms the other three machine learning approaches on the training set and ranks first, with SRC-AUC and PRC-AUC values of 93.14% and 91.81%, respectively. This model is followed by the SVM (88.36%), RF (86.30%), and LR (84.89%). It is to be noted that the order of the AUC values on the validation set is not the same. The CNN model has the highest PRC-AUC value (91.81%), and SVM (89.75%) has the highest value among the three machine learning methods, followed by LR (87.87%), and the worst is RF (86.14%). However, their overall degrees of difference are not significant, approximately five to six percentage points.
FIGURE 9
In Table 4, the precision is the proportion of successful classification of landslide locations (Wang et al., 2019). The recall is the proportion of landslide points successfully classified by the CNN model in all landslide locations (Liu et al., 2020). The F1-score combines the output of precision and recall, which is defined as the harmonic average of precision and recall (Fukuda et al., 2013; Li et al., 2022). The higher the F1-score is, the better the performance of the model will be. The Matthews correlation coefficient (MCC) is useful in comparing the binary classification of unbalanced datasets. In addition, its value range is −1 to 1; the closer to 1, the better of the classification result will be (Fukuda et al., 2013).
TABLE 4
| Using the validation dataset | ||||
|---|---|---|---|---|
| Method | Precision | Recall | F1-score | MCC |
| CNN | 0.99* | 0.91* | 0.94* | 0.77* |
| LR | 0.86 | 0.88 | 0.86 | 0.63 |
| RF | 0.78 | 0.89 | 0.83 | 0.65 |
| SVM | 0.97 | 0.91* | 0.93 | 0.76 |
Accuracy statistics of the models.
Note: * denotes the best performance.
In summary, the CNN model is feasible and effective in LSM, with excellent goodness-of-fit and predictive power compared to those of the other three machine learning models.
6 Conclusion
This study is one of the few works that analyzes the performance of CNN models in deep learning compared with common machine learning approaches, conducting a case study in Wenchuan County, China. We summarize the method of constructing deep learning datasets in detail, namely, using images as input, which is rarely found in present studies in the literature. To prepare the training samples, 11 landslide causative factors, namely, elevation, slope, aspect, relief amplitude, plan curvature, profile curvature, NDVI, lithology, distance to roads, distance to rivers, and distance to faults were selected. A total of 405 landslides were adopted in this research, and 70% of landslides (283 landslides) were chosen at random to perform the training models, with the remaining (122 landslides) to perform the validation. The SRC and PRC were used for the purpose of validating the results of the approaches based on the CNN, and logistic regression, random forest, and support vector machine were applied for comparison. The results indicated that the CNN-based approach has excellent goodness-of-fit and strong prediction capability compared with those of LR, RF, and SVM, exhibiting a superior SRC-AUC (93.14%) and PRC-AUC (91.81%), with successful application to landslide susceptibility mapping. In addition, compared with machine learning methods, CNN’s precision, recall, F1-score, and MCC also achieve the best results.
We only considered the region of Wenchuan for our study. It is necessary to conduct the investigation on the incidence of landslides in other areas of Sichuan for the purpose of obtaining information that are more relevant. In addition, the limitation of this study is that the terrain condition elements used here is all derived from the 30-m resolution DEM. In future research, higher resolution input data should be used to deal with the LSM, which may lead to better performance.
Statements
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
SZ: investigation, data curation, methodology, writing—original draft, algorithm design, and program development. LB: methodology, funding acquisition and supervision. YL: conceptualization, investigation, data curation, program development, and writing—review and editing. WL: funding acquisition and methodology. MX: data curation.
Funding
This study was supported by the National Key Research and Development Program of China (2021YFC3000401), the Nature Resource Science and Technology Program of Sichuan Province (No. Kj-2022-7), the China Geological Survey Project (No. DD20211365), and the Project of Remote Sensing Recognition and Monitoring of Potential Geological Hazard in Sichuan Province (No. 510201202110324).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AkgunA. (2012). A Comparison of Landslide Susceptibility Maps Produced by Logistic Regression, Multi-Criteria Decision, and Likelihood Ratio Methods: a Case Study at İzmir, Turkey. Landslides9, 93–106. 10.1007/s10346-011-0283-7
2
AkinciH.ZeybekM. (2021). Comparing Classical Statistic and Machine Learning Models in Landslide Susceptibility Mapping in Ardanuc (Artvin), turkey. Nat. Hazards1–29. 10.1007/s11069-021-04743-4
3
Al-NajjarH. A.KalantarB.PradhanB.SaeidiV. (2019). “Conditioning Factor Determination for Mapping and Prediction of Landslide Susceptibility Using Machine Learning Algorithms,” in Earth Resources and Environmental Remote Sensing/GIS Applications X (International Society for Optics and Photonics), 11156, 111560K. 10.1117/12.2532687
4
CaoJ.ZhangZ.DuJ.ZhangL.SongY.SunG. (2020). Multi-geohazards Susceptibility Mapping Based on Machine Learning-A Case Study in Jiuzhaigou, China. Nat. Hazards102, 851–871. 10.1007/s11069-020-03927-8
5
ChenM.TangC.XiongJ.ShiQ. Y.LiN.GongL. F. (2020a). The Long-Term Evolution of Landslide Activity Near the Epicentral Area of the 2008 Wenchuan Earthquake in china. Geomorphology367, 107317. 10.1016/j.geomorph.2020.107317
6
ChenT.ZhuL.NiuR.-q.TrinderC. J.PengL.LeiT. (2020b). Mapping Landslide Susceptibility at the Three Gorges Reservoir, china, Using Gradient Boosting Decision Tree, Random forest and Information Value Models. J. Mountain Sci.17. 10.1007/s11629-019-5839-3
7
ChongX.XiweiX.XiyanW.FuchuD.XinY.QiY. (2013). Detailed Catalog of Landslides Triggered by the 2008 Wenchuan Earthquake and Statistical Analyses of Their Spatial Distribution. J. Eng. Geology.21, 25–44.
8
ChungC.-J. F.FabbriA. G. (2003). Validation of Spatial Prediction Models for Landslide hazard Mapping. Nat. Hazards30, 451–472. 10.1023/b:nhaz.0000007172.62651.2b
9
ColkesenI.SahinE. K.KavzogluT. (2016). Susceptibility Mapping of Shallow Landslides Using Kernel-Based Gaussian Process, Support Vector Machines and Logistic Regression. J. Afr. Earth Sci.118, 53–64. 10.1016/j.jafrearsci.2016.02.019
10
CuiS.PeiX.JiangY.WangG.FanX.YangQ.et al (2021). Liquefaction within a Bedding Fault: Understanding the Initiation and Movement of the Daguangbao Landslide Triggered by the 2008 Wenchuan Earthquake (Ms= 8.0). Eng. Geology.295, 106455. 10.1016/j.enggeo.2021.106455
11
DaiF.LeeC.LiJ.XuZ. (2001). Assessment of Landslide Susceptibility on the Natural Terrain of Lantau Island, hong kong. Environ. Geology.40, 381–391.
12
DepickerA.JacobsL.DelvauxD.HavenithH.-B.MatesoJ.-C. M.GoversG.et al (2020). The Added Value of a Regional Landslide Susceptibility Assessment: The Western branch of the East African Rift. Geomorphology353, 106886. 10.1016/j.geomorph.2019.106886
13
DouJ.YunusA. P.BuiD. T.MerghadiA.SahanaM.ZhuZ.et al (2019). Assessment of Advanced Random forest and Decision Tree Algorithms for Modeling Rainfall-Induced Landslide Susceptibility in the Izu-Oshima Volcanic Island, japan. Sci. total Environ.662, 332–346. 10.1016/j.scitotenv.2019.01.221
14
FangZ.WangY.PengL.HongH. (2020). Integration of Convolutional Neural Network and Conventional Machine Learning Classifiers for Landslide Susceptibility Mapping. Comput. Geosciences139, 104470. 10.1016/j.cageo.2020.104470
15
FukudaS.De BaetsB.WaegemanW.VerwaerenJ.MoutonA. M. (2013). Habitat Prediction and Knowledge Extraction for Spawning European grayling (thymallus thymallus l.) Using a Broad Range of Species Distribution Models. Environ. Model. Softw.47, 1–6. 10.1016/j.envsoft.2013.04.005
16
GoetzJ.BrenningA.PetschkoH.LeopoldP. (2015). Evaluating Machine Learning and Statistical Prediction Techniques for Landslide Susceptibility Modeling. Comput. geosciences81, 1–11. 10.1016/j.cageo.2015.04.007
17
GuJ.WangZ.KuenJ.MaL.ShahroudyA.ShuaiB.et al (2018). Recent Advances in Convolutional Neural Networks. Pattern Recognition77, 354–377. 10.1016/j.patcog.2017.10.013
18
HeK.ZhangX.RenS.SunJ. (2016). “Deep Residual Learning for Image Recognition,” in Proceedings of the IEEE conference on computer vision and pattern recognition, 770–778. 10.1109/cvpr.2016.90
19
HeY.KusiakA. (2017). Performance Assessment of Wind Turbines: Data-Derived Quantitative Metrics. IEEE Trans. Sust. Energ.9, 65–73.
20
HepdenizK. (2020). Using the Analytic Hierarchy Process and Frequency Ratio Methods for Landslide Susceptibility Mapping in Isparta-Antalya Highway (D-685), turkey. Arabian J. Geosciences13, 1–16. 10.1007/s12517-020-05764-2
21
JaafariA.PanahiM.PhamB. T.ShahabiH.BuiD. T.RezaieF.et al (2019). Meta Optimization of an Adaptive Neuro-Fuzzy Inference System with Grey Wolf Optimizer and Biogeography-Based Optimization Algorithms for Spatial Prediction of Landslide Susceptibility. Catena175, 430–445. 10.1016/j.catena.2018.12.033
22
JiaoY.ZhaoD.DingY.LiuY.XuQ.QiuY.et al (2019). Performance Evaluation for Four Gis-Based Models Purposed to Predict and Map Landslide Susceptibility: A Case Study at a World Heritage Site in Southwest china. Catena183, 104221. 10.1016/j.catena.2019.104221
23
KavzogluT.SahinE. K.ColkesenI. (2014). Landslide Susceptibility Mapping Using Gis-Based Multi-Criteria Decision Analysis, Support Vector Machines, and Logistic Regression. Landslides11, 425–439. 10.1007/s10346-013-0391-7
24
KongC.WangJ.MaX.TianY.ZhangZ.XuK. (2021). Optimization of Random forest Model for Assessing and Predicting Geological Hazards Susceptibility in Lingyun County.
25
KrizhevskyA.SutskeverI.HintonG. E. (2012). Imagenet Classification with Deep Convolutional Neural Networks. Adv. Neural Inf. Process. Syst.25, 1097–1105.
26
LiH.HeY.XuQ.DengJ.LiW.WeiY. (2022). Detection and Segmentation of Loess Landslides via Satellite Images: a Two-phase Framework. Landslides, 1–14. 10.1007/s10346-021-01789-0
27
LiuP.WeiY.WangQ.ChenY.XieJ. (2020). Research on post-earthquake Landslide Extraction Algorithm Based on Improved U-Net Model. Remote Sensing12, 894. 10.3390/rs12050894
28
MengistuF.SuryabhagavanK.RaghuvanshiT. K.LewiE. (2019). Landslide hazard Zonation and Slope Instability Assessment Using Optical and Insar Data: a Case Study from Gidole Town and its Surrounding Areas, Southern ethiopia. Remote Sensing of Land3, 1–14. 10.21523/gcj1.19030101
29
MerghadiA.YunusA. P.DouJ.WhiteleyJ.ThaiPhamB.BuiD. T.et al (2020). Machine Learning Methods for Landslide Susceptibility Studies: A Comparative Overview of Algorithm Performance. Earth-Science Rev.207, 103225. 10.1016/j.earscirev.2020.103225
30
NandiA.ShakoorA. (2010). A Gis-Based Landslide Susceptibility Evaluation Using Bivariate and Multivariate Statistical Analyses. Eng. Geology.110, 11–20. 10.1016/j.enggeo.2009.10.001
31
NhuV.-H.MohammadiA.ShahabiH.AhmadB. B.Al-AnsariN.ShirzadiA.et al (2020a). Landslide Susceptibility Mapping Using Machine Learning Algorithms and Remote Sensing Data in a Tropical Environment. Int. J. Environ. Res. Public Health17, 4933. 10.3390/ijerph17144933
32
NhuV.-H.ShirzadiA.ShahabiH.SinghS. K.Al-AnsariN.ClagueJ. J.et al (2020b). Shallow Landslide Susceptibility Mapping: A Comparison between Logistic Model Tree, Logistic Regression, Naïve Bayes Tree, Artificial Neural Network, and Support Vector Machine Algorithms. Int. J. Environ. Res. Public Health17, 2749. 10.3390/ijerph17082749
33
NhuV.-H.ZandiD.ShahabiH.ChapiK.ShirzadiA.Al-AnsariN.et al (2020c). Comparison of Support Vector Machine, Bayesian Logistic Regression, and Alternating Decision Tree Algorithms for Shallow Landslide Susceptibility Mapping along a Mountainous Road in the West of iran. Appl. Sci.10, 5047. 10.3390/app10155047
34
OhH.-J.PradhanB. (2011). Application of a Neuro-Fuzzy Model to Landslide-Susceptibility Mapping for Shallow Landslides in a Tropical Hilly Area. Comput. Geosciences37, 1264–1276. 10.1016/j.cageo.2010.10.012
35
PhamB. T.PradhanB.BuiD. T.PrakashI.DholakiaM. (2016). A Comparative Study of Different Machine Learning Methods for Landslide Susceptibility Assessment: A Case Study of Uttarakhand Area (india). Environ. Model. Softw.84, 240–250. 10.1016/j.envsoft.2016.07.005
36
PhamB. T.PrakashI.ChenW.LyH.-B.HoL. S.OmidvarE.et al (2019). A Novel Intelligence Approach of a Sequential Minimal Optimization-Based Support Vector Machine for Landslide Susceptibility Mapping. Sustainability11, 6323. 10.3390/su11226323
37
PhamV. D.NguyenQ.-H.NguyenH.-D.PhamV.-M.BuiQ.-T.et al (2020). Convolutional Neural Network—Optimized Moth Flame Algorithm for Shallow Landslide Susceptible Analysis. IEEE Access8, 32727–32736. 10.1109/access.2020.2973415
38
PourghasemiH.MoradiH.AghdaS. F. (2013). Landslide Susceptibility Mapping by Binary Logistic Regression, Analytical Hierarchy Process, and Statistical index Models and Assessment of Their Performances. Nat. hazards69, 749–779. 10.1007/s11069-013-0728-5
39
QiangX.Xiu-JunD. (2011). Genetic Types of Large-Scale Landslides Induced by Wenchuan Earthquake. Earth Science-Journal China Univ. Geosciences36, 1134–1142.
40
SahaA.SahaS. (2020). Comparing the Efficiency of Weight of Evidence, Support Vector Machine and Their Ensemble Approaches in Landslide Susceptibility Modelling: A Study on Kurseong Region of Darjeeling Himalaya, india. Remote Sensing Appl. Soc. Environ.19, 100323. 10.1016/j.rsase.2020.100323
41
SahaS.RoyJ.HembramT. K.PradhanB.DikshitA.Abdul MauludK. N.et al (2021). Comparison between Deep Learning and Tree-Based Machine Learning Approaches for Landslide Susceptibility Mapping. Water13, 2664. 10.3390/w13192664
42
SevgenE.KocamanS.NefesliogluH. A.GokceogluC. (2019). A Novel Performance Assessment Approach Using Photogrammetric Techniques for Landslide Susceptibility Mapping with Logistic Regression, Ann and Random forest. Sensors19, 3940. 10.3390/s19183940
43
ShanoL.RaghuvanshiT. K.MetenM. (2020). Landslide Susceptibility Evaluation and hazard Zonation Techniques–A Review. Geoenvironmental Disasters7, 1–19. 10.1186/s40677-020-00152-0
44
SunD.XuJ.WenH.WangY. (2020). An Optimized Random forest Model and its Generalization Ability in Landslide Susceptibility Mapping: Application in Two Areas of Three Gorges Reservoir, china. J. Earth Sci.31, 1068–1086. 10.1007/s12583-020-1072-9
45
SzegedyC.LiuW.JiaY.SermanetP.ReedS.AnguelovD.et al (2015). “Going Deeper with Convolutions,” in Proceedings of the IEEE conference on computer vision and pattern recognition, 1–9. 10.1109/cvpr.2015.7298594
46
Van WestenC. J.CastellanosE.KuriakoseS. L. (2008). Spatial Data for Landslide Susceptibility, hazard, and Vulnerability Assessment: An Overview. Eng. Geology.102, 112–131. 10.1016/j.enggeo.2008.03.010
47
VapnikV. N. (1999). An Overview of Statistical Learning Theory. IEEE Trans. Neural networks10, 988–999. 10.1109/72.788640
48
WangY.FangZ.HongH. (2019). Comparison of Convolutional Neural Networks for Landslide Susceptibility Mapping in Yanshan County, china. Sci. total Environ.666, 975–993. 10.1016/j.scitotenv.2019.02.263
49
WangY.SunD.WenH.ZhangH.ZhangF. (2020). Comparison of Random forest Model and Frequency Ratio Model for Landslide Susceptibility Mapping (Lsm) in Yunyang County (Chongqing, china). Int. J. Environ. Res. Public Health17, 4206. 10.3390/ijerph17124206
50
WeiX.ZhangL.LuoJ.LiuD. (2021). A Hybrid Framework Integrating Physical Model and Convolutional Neural Network for Regional Landslide Susceptibility Mapping. Nat. Hazards109, 471–497. 10.1007/s11069-021-04844-0
51
WuX.ChenX.ZhanF. B.HongS. (2015). Global Research Trends in Landslides during 1991–2014: a Bibliometric Analysis. Landslides12, 1215–1226. 10.1007/s10346-015-0624-z
52
XieM.JuN.ZhaoJ.FanQ.HeC. (2021a). Comparative Analysis on Classification Methods of Geological Disaster Susceptibility Assessment. Geomatics Inf. Sci. Wuhan Univ.46, 1003–1014.
53
XieW.NieW.SaffariP.RobledoL. F.DescoteP.-Y.JianW. (2021b). Landslide hazard Assessment Based on Bayesian Optimization–Support Vector Machine in Nanping City, china. Nat. Hazards, 1–18.
54
YiY.ZhangZ.ZhangW.JiaH.ZhangJ. (2020). Landslide Susceptibility Mapping Using Multiscale Sampling Strategy and Convolutional Neural Network: A Case Study in Jiuzhaigou Region. Catena195, 104851. 10.1016/j.catena.2020.104851
55
ZhangY.-x.LanH.-x.LiL.-p.WuY.-m.ChenJ.-h.TianN.-m. (2020). Optimizing the Frequency Ratio Method for Landslide Susceptibility Assessment: A Case Study of the Caiyuan basin in the Southeast Mountainous Area of china. J. Mountain Sci.17, 340–357. 10.1007/s11629-019-5702-6
56
ZhouJ.WeiJ.YangT.ZhangP.LiuF.ChenJ. (2021). Seepage Channel Development in the crown Pillar: Insights from Induced Microseismicity. Int. J. Rock Mech. Mining Sci.145, 104851. 10.1016/j.ijrmms.2021.104851
Summary
Keywords
convolutional neural network, landslide susceptibility, machine learning, GIS, Wenchuan County
Citation
Zhang S, Bai L, Li Y, Li W and Xie M (2022) Comparing Convolutional Neural Network and Machine Learning Models in Landslide Susceptibility Mapping: A Case Study in Wenchuan County. Front. Environ. Sci. 10:886841. doi: 10.3389/fenvs.2022.886841
Received
01 March 2022
Accepted
11 April 2022
Published
16 May 2022
Volume
10 - 2022
Edited by
Jingren Zhou, Sichuan University, China
Reviewed by
Zhongping Yang, Chongqing University, China
Yueqin Zhu, Ministry of Emergency Management, China
Updates
Copyright
© 2022 Zhang, Bai, Li, Li and Xie.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lin Bai, bailin@cdut.edu.cn
This article was submitted to Environmental Informatics and Remote Sensing, a section of the journal Frontiers in Environmental Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.