Next Article in Journal
Estimation of Natural Radionuclides’ Concentration of the Plutonic Rocks in the Sakarya Zone, Turkey Using Multivariate Statistical Methods
Previous Article in Journal
A Feature Selection Model for Network Intrusion Detection System Based on PSO, GWO, FFA and GA Algorithms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Landslide Susceptibility Mapping Using the Slope Unit for Southeastern Helong City, Jilin Province, China: A Comparison of ANN and SVM

1
College of Construction Engineering, Jilin University, Changchun 130026, China
2
Jilin Team of Geological Survey Center of China Building Materials Industrial, Changchun 130026, China
*
Author to whom correspondence should be addressed.
Submission received: 1 June 2020 / Revised: 19 June 2020 / Accepted: 21 June 2020 / Published: 23 June 2020

Abstract

:
The purpose of this study is to produce a landslide susceptibility map of Southeastern Helong City, Jilin Province, Northeastern China. According to the geological hazard survey (1:50,000) project of Helong city, a total of 83 landslides were mapped in the study area. The slope unit, which is classified based on the curvature watershed method, is selected as the mapping unit. Based on field investigations and previous studies, three groups of influencing Factors—Lithological factors, topographic factors, and geological environment factors (including ten influencing factors)—are selected as the influencing factors. Artificial neural networks (ANN’s) and support vector machines (SVM’s) are introduced to build the landslide susceptibility model. Five-fold cross-validation, the receiver operating characteristic curve, and statistical parameters are used to optimize model. The results show that the SVM model is the optimal model. The landslide susceptibility maps produced using the SVM model are classified into five grades—very high, high, moderate, low, and very low—and the areas of the five grades were 127.43, 151.60, 198.77, 491.19, and 506.91 km2, respectively. The very high and high susceptibility areas included 79.52% of the total landslides, demonstrating that the landslide susceptibility map produced in this paper is reasonable. Consequently, this study can serve as a guide for landslide prevention and for future land planning in the southeast of Helong city.

1. Introduction

Landslides are one of the most common geological disasters in the mountainous areas of China [1,2,3]. In recent years, the Chinese economy has developed steadily and rapidly, and has become the second largest economy in the world. With such rapid development of the economy, many previously inaccessible areas have carried out corresponding infrastructure construction and various engineering activities in succession. Therefore, it is necessary to evaluate the risk of natural disasters such as landslides in these areas. Landslide susceptibility mapping is the basis of landslide hazard and risk assessment [4]. The purpose of landslide susceptibility mapping is to answer the question “what is the geological background of landslides and where are the areas which are most prone to them?” [2]. Based on the review of relevant papers, it was found that landslide susceptibility mapping mainly includes the following [5,6]: (a) landslide inventory data acquisition; (b) mapping unit selection; (c) influencing factor selection; (d) establishment of the evaluation model; and (e) production of the landslide susceptibility map.
Landslide inventory data is the basis of landslide susceptibility mapping. The early acquisition of landslide inventory data mainly depends on field geological disaster surveys. With the development of remote sensing technology, landslide interpretation based on remote sensing technology has become one of the most popular tools to quickly acquire landslide inventory data [7,8]. Remote sensing techniques commonly used in landslide interpretation mainly include (a) visible optical remote sensing and (b) InSAR (Interferometric Synthetic Aperture Radar) remote sensing. The key to remote sensing interpretation of landslides based on visible optics is to interpret the landslide based on the special topographic and geomorphic features of the landslide in the remote sensing image. However, this method is greatly affected by the spatial resolution of remote sensing data and is unable to identify the surface deformation characteristics (i.e., the potential and slow deformation of the landslide); thus, it has certain limitations. However, landslide interpretation based on InSAR technology can effectively identify the surface deformation and identify the potential landslides that are under deformation [9], which makes up for the limitations of visible optical remote sensing interpretation. The combination of these two methods has become a common method for landslide interpretation. As for mapping units, they are the basic units in landslide susceptibility mapping [10]. At present, the commonly used mapping units are the grid unit, watershed unit, slope unit, region unit, and uniform condition unit [4]. Among these, the grid unit is the most widely used but, due to its lack of connection with the topography and other geological information, the slope unit [11,12] has become more and more respected by relevant scholars. The selection of influencing factors is one of the key issues in the study of landslide susceptibility mapping [13]. The choice of influencing factors used in the evaluation will affect the final mapping results. The selection of influencing factors mostly relies on field investigation, cause analysis, expert experience, and related research of landslides; the factors with the highest correlation with the occurrence of landslides in the study area are typically selected as the influencing factors to evaluate the susceptibility to landslides [2]. At present, many scholars also use methods of automatic forward selection of parameters that weight and evaluate the effectiveness of each influencing factor and make use only of the “best” ones [14,15]. The evaluation models applied to landslide susceptibility mapping are mainly divided into qualitative and quantitative evaluation [6,16]. Qualitative evaluation is mainly based on the knowledge and experience of relevant experts, where the susceptibility to geological hazards such as landslides is determined by means of grading and weighting [17,18]. Quantitative evaluation methods include statistics, machine learning, deep mining, and other methods [6]. With the rapid improvement of calculation ability, an increasing number of quantitative evaluation methods have been applied to landslide susceptibility mapping. In addition, evaluation methods of the performance of landslide susceptibility models have also been developed rapidly [16,19]. Finally, the landslide susceptibility map can be produced according to the established model.
According to the basic steps of landslide susceptibility mapping, this paper carries out a relative evaluation of the southeast of Helong City, Jilin Province, northeastern China. Firstly, according to the geological hazard survey (1:50,000) project, 83 landslides were mapped in the study. Then, the slope unit classified by the curvature watershed method was applied as the mapping unit, which can identify the horizontal surface and the inclined surface. Moreover, the area of the slope unit is more concentrated, the shapes mostly range between isosceles triangles and squares, and there are few elongated units. Next, ten influencing factors were selected. Artificial neural networks (ANN) and support vector machines (SVM) were introduced to build the landslide susceptibility model. Five-fold cross-validation, the receiver operating characteristic curve (ROC), and statistical parameters were used to optimize the model. Finally, the landslide susceptibility map was divided into five classes. This study can serve as a guide for landslide prevention and for future land planning in the southeast of Helong city.

2. Study Area

The study area is in the Chongshan and Nanping counties of the southeast Helong city, Jinlin province, China, along the Tumen river. The study area is marked by valleys and rugged topography, and covers an area of 1476 km2 (Figure 1a–c). More than 90% of the total area of the study area consists of mountainous areas. In the study area, the maximum elevation is 1450 m, and the minimum elevation is 350 m. The southern part of the study area is the Chinese and Korean quasi-platform, and the northern part is the Jihei fold system in the Tianshan-Xingan geosyncline fold area, which are bounded by the deep and large fault of Gudong River. The study area belongs to the temperate monsoon sub-humid climate zone. Based on rainfall data collected from 1960 to 2012, the maximum daily rainfall of the study area is 164.8 mm. The perennial average temperature is 5.6 °C. The vegetation coverage in the study area is relatively high. The seismic intensity of the study area has a degree of VI on the modified Mercalli index. At present, no earthquake-induced landslides have been detected in the study area. According to the field investigation, the landslides in the study area are mainly distributed along the Tumen River. According to a geological map (downloaded using the 91 Weitu software, with a scale of 1:500,000) (Figure 1d), it can be seen that the mainly exposed strata in the study area are Quaternary (Q), Neogene (N), Cretaceous (K), Jurassic (J), Middle Proterozoic (Pt), and New Archean (Ar). The lithological information of the study area is listed in Table 1. There are several reverse faults in the study area, which affect the regional stability. The Tumen river and its tributaries flow through the study area, which has a significant impact on landslide risk.

3. Methodology

Figure 2 shows the three main steps of this paper, namely (a) data preparation, (b) landslide susceptibility modelling, and (c) validation and selection of the optimal models.

3.1. The Mapping Unit

Reasonable mapping units should be determined before landslide susceptibility mapping is carried out [20]. Whether the mapping unit is reasonable or not directly affects the accuracy of the final evaluation result. Slope units are often used to study landslide susceptibility mapping due to their close relationship with the topography. The dividing principle of slope units is to divide the research area into many small areas with different sizes by cutting along the ridge and valley lines [21]. At present, hydrologic analysis is the most used method for dividing slope units. However, this method cannot identify the boundary between the horizontal and the inclined surface and generates a large number of parallel river networks, which increases the difficulty of manual modification [21]. By comparing the slope unit division results of the curvature watershed method and the hydrologic analysis method, it was found that the slope unit divided based on the curvature watershed method has a uniform size, a regular shape, and small terrain variation inside, which is obviously better than the result based on the hydrological method. Therefore, we chose the curvature watershed method to divide the slope units in this study. Slope units are divided by valley and ridge lines, where there are abrupt changes in slope angle and slope aspect. Therefore, slope units can be divided according to the change of slope angle and slope aspect. Profile curvature and plan curvature are the derivatives of slope angle and slope aspect, respectively. Their maximum and minimum values can reflect abrupt changes of slope angle and slope aspect. Therefore, the curvature can be used to divide the slope units. The specific process of dividing slope units based on curvature is shown in Figure 3. It is mainly used to identify the boundary of concave terrain and convex terrain by the curvature and reverse curvature, respectively, to divide the slope units. The classification steps of the slope units can be divided into two parts: positive relief extraction and negative relief extraction. The original curvature was used to extract the positive relief, and the flip curvature was used to extract the negative relief, which can be obtained by multiplying the curvature by −1. Before calculating the curvature, to remove the influence of the roughness of the original Digital Elevation Model (DEM) on the curvature calculation results, focal statistics were made on the DEM with three pixels as the radius. After the flow direction and sink calculations, the positive and negative relief boundary can be obtained. The results of dividing slope units can be obtained by merging the positive and negative relief boundaries and manually modifying the unreasonable units.

3.2. Landslide Inventory

Landslide inventory data is the basis of landslide susceptibility mapping. In order to grasp the spatial distribution characteristic information and develop a mechanism of landslide hazards in the study area, the Jilin team of the Geological Survey Center of China Industrial Building Materials carried out the geological hazard survey (1:50,000) project of Helong city. The characteristic geomorphological features of the landslides, such as the chair-like landform of the landslide, a color difference between the landslide and the surrounding environment, and a change of vegetation before and after the landslide, can be clearly identified in a remote sensing image. Thus, the project team preliminarily identified 112 landslides in the study area using the visible optical remote sensing technology and Google images. Then, during the field geological survey, the remote sensing interpretation results were reviewed, and 42 interpretation results were excluded (such as exposed steep cliffs of bedrock and a waste rock dump in a quarry) and the uninterpreted landslides were supplemented. Finally, a total of a total of 83 landslides were mapped in the study area (Figure 1c), including soil and rock slope deformation, soil slide, and collapse. Figure 1e,f are examples of the landslides. Table A1 shows the basic information for the landslides. According to Table A1, the landslides in the study area are mostly rock landslides, which are mainly developed in crystalline rocks. These landslides are mainly controlled by the rock mass structural plane, and a small number of landslides are controlled by the contact surface between the overburden layer and the rock mass. The scale of landslide is mainly small and medium size. The failure mode of the landslides is mainly pull-type and toppling-type.

3.3. Influencing Factors

The occurrence of a landslide is a complicated process. Therefore, the factors influencing the occurrence of a landslide are also varied. Pourghasemi and Rossi [13] reviewed a total of 220 related papers and found that the factors with the highest application frequency in landslide susceptibility mapping were slope angle, lithology, slope aspect, land use, distance to river, elevation, distance to faults, curvature, distance to road, and soil type, among others. The selection of influencing factors for the landslide susceptibility mapping in the study area should be based on the understanding of the characteristics of the landslide and the geological environment. Therefore, the relationship between the landslide and the geological environment in the study area was analyzed as follows.

3.3.1. Relationship between Geological Environment and Landslides

Relationship between Topographic Features and Landslides

The topography is the basic condition of a landslide and determines the formation and development of landslide to a great extent. For example, convex and bedding slopes are prone to landslides. Based on field investigation and analysis, it was found that most landslides in the study area occurred in hills and low mountains with an elevation range of 200–1000 m and a slope angle more than 30° (Figure 4). Furthermore, the number of landslides increased with the increase of the slope angle. The number of landslides in a convex slope was much larger than that in concave slope. The micro-landform of the landslide site was a steep slope or steep coast (Table A1).

Relationship between Lithological Features and Landslides

The lithology of landslides in the study area was mainly crystalline rocks, such as granite, diorite, and other magmatic rocks (Table A1). The number of landslides in the crystalline rock in the study area was the largest and was obviously higher than that in other geotechnical rock type distribution areas. This is because the crystalline rock is usually strongly influenced by joints, cracks, and faults. The presence of these structural planes greatly reduces the shear strength of the rock mass, resulting in the sliding of the upper rock mass or overlying loose rock.

Relationship between Geologic Features and Landslides

The number of landslides in fault geological tectonic units was the largest in the study area. This is mainly because the in-situ stress of the fracture is stronger than that of other parts, and its influence range is larger and wider. The rock in these parts has deformations such as bending, squeezing, and tearing, which reduces the structural mechanical properties of the rock mass.

Relationship between Rainfall Features and Landslides

According to statistics, the landslides that occurred in the study area were mostly caused by heavy rainfall, and the two are positively correlated. Continuous heavy rainfall, particularly during the flood season, results in a high incidence of landslides.

Relationship between Other Features and Landslides

Landslides occurred more frequently in areas with low vegetation coverage in the study area. Landslides occurred extensively along the Tumen river, which is the main river in the study area. Earthquakes are an important trigger of landslides. They damage the integrity of rock and soil mass and make slopes prone to landslide. However, there is no historical record of landslides triggered by earthquakes in the study area. Thus, seismic activity has had little effect on the occurrence of landslides in the study area.

3.3.2. Selection of Influencing Factors

Based on the above analysis, and the statistical results of Pourghasemi and Rossi [13], three groups of influencing factors were selected for this study: geologic factors, topographic factors, and environment factors. The geologic factors only consist of geology, while there were five topographic factors: elevation, slope angle, slope aspect, topographic relief, and curvature. Finally, there were four environment factors: land use, rainfall, distance to river, and distance to faults. In this study, all the influencing factor maps were extracted using the slope unit. Furthermore, the dominant category of lithology, land use, was estimated as the influencing factor value of each slope unit, with the average value used for the other factors.

3.3.3. Extraction of Influencing Factors

Geologic Factor

The study area is covered with six types of geologic formations, composed of different lithologies pertaining to different geologic ages. The combination of these two characteristics is one of the most important factors influencing the occurrence of landslides [22]. The shear strength, weathering resistance, and crushing degree of different lithologies are greatly different [23]. In addition, the same lithology with different structural planes will have different effects on the stability of the slope. A rock mass with joint and fissure development is more prone to landslides than a complete rock mass. A highly weathered rock mass is also less stable than a fresh rock mass. In this study, the geology map was obtained based on a geological map with a scale of 1:500,000 (downloaded using 91 Weitu software); the geology map of the study area is shown in Figure 5a.

Topographic Factors

Elevation usually affects the depth of the water table aquifer [1,24]. In addition, the higher the elevation, the less human activity takes place, which is also more conducive to the stability of the slope. Regarding the slope angle, a change of slope angle changes the original stress state of the slope [25]. Within a certain range of slope, the self-gravity stress and shear stress in the slope generally increase with an increase of the slope angle, and the probability of slope instability increases accordingly. The intensity of light exposure, the type and extent of vegetation cover, and the supply of surface water vary greatly with differing slope aspect [2,26]. For example, the illumination time on a sunny slope is much longer than that on a shady slope; therefore, the temperature difference between day and night on a sunny slope is also larger than that on a shady slope and the dry–wet cycle is also faster. In this case, the weathering strength of the rock mass on a sunny slope is larger than that on a shady slope, which reduces the strength and stability of rock and soil mass on the sunny slope, in turn increasing the probability of landslide [4]. The topographic Relief—The difference between the highest and lowest Points—Can reflect the degree of relief in a specific area [4]. Curvature can reflect the shape of a slope body. According to the shape characteristics, a slope can be divided into three forms: convex slope, straight slope, and concave slope [27]. A DEM (Digital Elevation Model) with a resolution of 10 × 10 m was used to extract the five influencing factor maps through the slope unit using ArcGis software. The five influencing factor maps are shown in Figure 5b–d and Figure 6a,b.

Environment Factors

Vegetation has a positive effect on the stability of a slope, which can enhance the resistance of the slope surface to water erosion and influence rainfall infiltration. Rainfall has a significant impact on slope risk. Rainwater infiltration will increase the gravity of rock and soil mass, reducing their shear strength parameters, which will have an adverse effect on slope stability. The distance to river factor has an important influence on the occurrence of landslides [28]. The downcut of a river increases the slope angle along the river. To adapt to the rapid downcut, disasters such as landslides often occur along the river, thus reducing the slope angle. The existence of faults leads to the development of joints and fissures in the surrounding rock mass, which leads to fragmentation of the rock mass and a reduction of weathering resistance [29]. Therefore, many geological hazards such as landslides often develop around deep and large fissures. These four influencing factors are extracted through the slope unit using ArcGis software. The five influencing factor maps are shown in Figure 6c,d and Figure 7a,b.

3.4. Landslide Susceptibility Modeling

At present, an increasing number of models are being applied to the study of landslide susceptibility mapping [30,31,32,33]. Each of these models has its advantages and disadvantages; combined with the complexity of the factors affecting slope stability and the different geological environments in different regions, it is impossible to apply one model to the study of landslide susceptibility in all areas. Therefore, in the relevant evaluation of the study area, the evaluation model adopted should be optimized to determine the most suitable landslide susceptibility evaluation model. Recently, some machine learning approaches have been developed with the aim of evaluating the landslide susceptibility [21]. Due to its strong learning ability, fast calculation speed, and strong fault tolerance, the artificial neural network (ANN) model has become one of the most used machine learning methods in landslide susceptibility mapping. By comparison, the support vector machine (SVM) model requires fewer modeling data and is suitable for binary classification problems. Thus, the SVM model may be more suitable for landslide susceptibility prediction with fewer data, such as in the case of this study. Therefore, we selected the artificial neural network (ANN) and support vector machine (SVM) models to optimize the evaluation model of landslide susceptibility in the study area.

3.4.1. Artificial Neural Network (ANN)

As a machine learning method, the artificial neural network model has been widely used in the study of landslide susceptibility mapping [26,34,35]. The ANN model has the following advantages [36]: (a) better non-linear mapping capability; (b) highly self-learning and adaptive; (c) strong generalization ability; and (d) strong fault tolerance. Thus, the ANN model can simulate the complex non-linear interactions between influencing factors and landslides through the interaction between neurons. Moreover, in an ANN model it is not necessary to describe the interactions between factors by complex mathematical formulae [37]. The fitting effect is relatively good, which is especially suitable for the simulation of complex geological phenomena influenced by many factors and presents advantages for the simulation of phenomena such as landslides, in which various factors interact and have complex relationships. A complete artificial neural network model is usually composed of an input layer, an output layer, and one or more hidden layers. Figure 8 shows a schematic diagram of a simple artificial neural network model. Each layer contains several of the model’s basic units: neurons. The data processing and storage of an artificial neural network are represented as the mutual relationships and connections between neurons. Through training using samples, the ANN changes the weights of its internal connections to minimize the error between the output value and the target value, to achieve the purpose of accurate modeling. Each node in the input layer corresponds to a predictive variable, while the node in the output layer corresponds to the target variable. A hidden layer is a regular layer connecting the input layer and output layer, and the number of the hidden layers and the number of nodes in each layer determine the complexity of the network. The process of forecasting with an ANN can be divided into a learning process and a prediction process [38]. The learning process, in which a large number of learning samples and the iterative function of the network are used to train the network by optimizing the process via minimizing the network error, includes the forward transmission of input information and the reverse transmission of error [34]. The prediction process is the process of substituting the unknown sample into the model after training. According to the rules of the learning sample, the process finally obtains the output result of the sample through the forward transmission of the input information [39].

3.4.2. Support Vector Machine (SVM)

The support vector machine (SVM) model is also a machine learning model. Its theoretical basis is statistical learning theory [40,41,42,43]. The SVM model has the following advantages [40,44]: (a) low data volume requirement; (b) strong generalization ability; (c) strong optimization ability; (d) adaptability to high-dimensional samples; and (e) strong learning ability and fast convergence. It can realize the linear segmentation of data by transforming each evaluation index from low dimensional space to high dimensional space. Thus, it can analyze and evaluate non-linear problems in low dimensional space. Due to its low requirement in terms of data volume, it has also been widely used in the study of landslide susceptibility mapping. The process of modeling for landslide susceptibility evaluation by SVM is summarized as follows [40,44]:
Consider some linearly separable data points xi (i = 1, 2,…,n) that fall into two different classes yi = ±1. The goal of SVM is to find a hyperplane in n-dimensional data space which can separate the two classes of data based on the maximum interval. The hyperplane can be expressed mathematically as follows:
L = 1 2 | | w | | 2 ,
which should satisfy the following constraint conditions:
y i ( ( w · x i ) + b ) 1 ,
where ||w|| is the norm of the normal vector of the hyperplane, b is a scalar, and represents the scalar product. Based on the Lagrange multiplier, the cost function can be expressed as follows:
L = 1 2 | | w | | 2 i = 1 n λ i ( y i ( ( w · x i ) + b ) 1 ) ,
where λi is the Lagrange multiplier.
In the case of linear indivisibility, the constraint condition can introduce a slack variable ξi, which can be expressed as follows:
y i ( ( w · x i ) + b ) 1 ξ i .
Then, Equation (1) can be converted into the following form:
L L = 1 2 | | w | | 2 1 v n i = 1 n ξ i .
In Equation (5), v in [0,1] is introduced to consider the case of misclassification. In addition, the kernel function K (xi, yi) is introduced to explain the non-linear decision boundary problem in SVM.

3.5. Data for Landslide Susceptibility Modeling

When using an artificial neural network or support vector machine to establish a landslide susceptibility model, the same amount of data is needed for landslide units and non-landslide units. In this study, we determined the number of slope units with landslides according to the division results of slope units in the study area and the landslide inventory map. An equal number of units were randomly selected at a minimum distance of 800 m from these units, in order to avoid the effects of landslides [37]. The five-fold cross-validation method [39,40] was used to validate the models and to overcome the shortage of landslide data and the problem of model overfitting. All the ten influencing factors were involved in landslide susceptibility model building.

3.6. Validation Method

3.6.1. Receiver Operating Characteristic Curve (ROC)

The receiver operating characteristic curve (ROC) [18,45] is a quantitative analysis method to evaluate the prediction accuracy of the landslide susceptibility model. This method evaluates the prediction accuracy of the model using the area under the curve (AUC). The AUC value lies between 0 and 1, and the greater its value, the higher the prediction accuracy of the model.

3.6.2. Statistical Analysis Method

Statistical indices are also widely used for evaluating the prediction ability of landslide susceptibility mapping models [21]. The most used statistical indices are the following:
A c = ( T P + T N ) / ( T P + T N + F P + F N ) ,
S e n = T P / ( T P + F N ) ,
S p = T N / ( T N + F P ) ,
P P V = T P / ( T P + F P ) ,
N P V = T N / ( T N + F N )
where Ac is the accuracy; Sen is the sensitivity; Sp is the specificity; PPV is the positive predictive value; NPV is the negative positive value; TP is the true positive; TN is the true negative; FP is the false positive; and FN is the false negative.

4. Results

4.1. Division Result of the Slope Units

In this paper, a DEM with a resolution of 10 × 10 m was adopted to the divided slope units of the study area. In the process of division, we found that the size of slope units divided by the curvature watershed method was related to the DEM resolution. Therefore, the DEM resolution was converted to 10 × 10 m, 30 × 30 m, 50 × 50 m, 80 × 80 m, 100 × 100 m, and 120 × 120 m for slope classification. By comparing the slope unit classification results with Google images of the study area, it was found that the slope unit classification result was the most consistent with the actual terrain with the DEM with a resolution of 80 × 80 m. A total of 2956 slope units were obtained (Figure 9). The maximum area of a slope unit was 18.45 × 105 m2 and the minimum area was 0.11 × 105 m2.

4.2. Model Fitting Results

According to the landslide inventory map, there are 83 slope units containing the entire known landslide body. Therefore, the 83 slope units that experienced landslides were used as the modeling data. To meet the modeling requirements, the same number of non-landslide units (83) were randomly selected at least 800 m away from the landslide units (Figure 9). To establish the ANN model, the numbers of input, hidden, and output layers should be determined first. In this study, each of the input, hidden, and output layers consisted of a single layer. Secondly, the number of neurons in each layer was determined, with the number of neurons in the input layer the same as the number of influencing factors. There were ten influencing factors used for modeling in this study, thus, the number of neurons in the input layer was ten. The output layer was used to determine whether a landslide occurs, so the number of neurons was two. The number of hidden layer neurons can be determined by the following empirical formula [46]:
N = A + B + k
where N is the recommended value for the number of neurons in the hidden layer, A is the number of neurons in the input layer, B is the number of neurons in the output layer, and k is an empirical coefficient with value between 0 and 10. According to the empirical formula, the ideal number of hidden layer neurons in the artificial neural network in this study ranged from 4 to 14. By using all the data to establish the ANN model, the number of hidden layer neurons was optimized. It can be seen from Figure 10 that, when the number of hidden layer neurons is 6, the ANN model has the highest prediction accuracy. Thus, the number of neurons in the hidden layer was finally selected as 6.
For the ANN model, the learning rate, momentum, and training time were set as 0.3, 0.3, and 500, respectively [34,37]. The choice of kernel function affects the prediction accuracy of SVM. In this study, the radial basis function (RBF) was selected as the kernel function, which was influenced by the regularization parameter (C) and the kernel parameter (g). C and g were set as 0.8 and 0.5 [34,47,48], respectively. All the parameters were obtained based on the previous research experience and experiments in the calculation process. The model fitting results are shown in Table 2.

4.3. Landslide Susceptibility Mapping Results

By comparing the five statistical parameters and AUC values of the two models, the SVM model was determined to be the optimal model. Therefore, the SVM model was used in this study to produce the landslide susceptibility map of the study area. The final model was established using the model with high accuracy and AUC value in the process of five-fold cross-validation. The natural breaks method was used to divide the landslide susceptibility of the study area into five grades: Very low, low, moderate, high, and very high. The landslide susceptibility map of the study area is shown in Figure 11.
Figure 11 and Table 3 show that the areas of the five susceptibility classes for the ANN model (very high, high, moderate, low, and very low) were 146.24, 297.95, 423.33, 310.44, and 297.95 km2, respectively. For landslide occurrence, the number of landslides in the five susceptibility classes were 43, 18, 10, 6, and 6, respectively. For the SVM model, the areas of the five susceptibility classes were 127.43, 151.60, 198.77, 491.19, and 506.91 km2, respectively. For landslide occurrence, the number of landslides in the five susceptibility classes were 52, 14, 8, 4, and 5, respectively.

5. Discussion

5.1. Slope Unit Classification Results

To evaluate the effect of dividing slope units, the area and shape indices of slope units were statistically analyzed. The ramp unit was used to extract the various influences, so the area it covers should be approximately the same size overall. The uniformity of slope unit area can be reflected by the distribution of slope unit area. It can be seen, from Figure 12a, that the slope unit area was concentrated between 4 and 8 × 105 m2, accounting for 44.2% for the total number of slope units. If the slope unit is too flat, or there is an elongated unit, the uniformity inside the unit will be destroyed to a great extent. The shape index can be used to evaluate the shape of the slope unit. The shape index of the slope units can be calculated using the following equation:
S = L ^ 2 / 4 π A ,
where S is the shape index, L is the perimeter of the slope units, and A is the area of the slope units. According to the definition of the shape index, the shape index of a circle is 1, that of a square is 1.27, and that of an equilateral triangle is 1.59. Based on Figure 12b, the slope units with shape index below 1.59 accounted for 81.0%, which means that most of the slope units were between circle and equilateral triangle shapes; however, there were a few elongated slope units. Thus, the effect of dividing slope units was reasonable.
In addition, in the process of dividing slope units, it was found that compared with the hydrologic subdivision method, there is no large amount of parallel river network in the process of dividing slope units by the curvature watershed method. Moreover, the number of unreasonable units generated is much less than the hydrologic analysis method, so the later manual modification work is much less.

5.2. Comparison between ANN and SVM Model

Landslide susceptibility mapping is a popular research issue due to its non-linear characteristics. Although many mathematical methods have been applied to landslide susceptibility mapping, the prediction accuracy of these models is not very stable. Therefore, the same model cannot be applied to all studies. As a result, one of the key problems in landslide susceptibility mapping is to find a landslide susceptibility model which is suitable for a specific study area. In this study, the ANN and SVM models were introduced to establish a landslide susceptibility model for the area southeast of Helong city.
Through the training and testing of the two selected models, a confusion matrix was obtained, and the corresponding statistical parameters of each model were calculated to evaluate their respective advantages and disadvantages. From Table 2, the mean AUC values of ANN and SVM models differ greatly: in the training stage they were 88.60% and 93.22%, respectively; in the testing stage they were 84.52% and 89.74%, respectively. The AUC values of the two models decreased in the testing stage (by 4.08% for the ANN model and 3.48% for the SVM model), and that of the ANN model decreased significantly. According to the mean AUC value, the SVM model was slightly better than the ANN model (4.62% for the training stage, 5.22% for the testing stage). For the standard deviation of the AUC value in the training stage, that of the ANN model was obviously larger than that of the SVM model (2.48% for the ANN model, 0.83% for the SVM model). This indicated that the stability of the SVM model was better than that of ANN model in the training stage. In the testing stage, for the ANN model, the standard deviation of the AUC value (2.88%) increased slightly, and similarly for the SVM model (1.50%); however, the SVM model was more stable. From the mean AUC value alone, we believe that the SVM model could be the best.
For the statistical parameters, in the training stage, the mean accuracies were 83.27% and 88.56%, respectively, for the ANN model and the SVM model; in the testing stage, these were 71.23% and 84.30%. The standard deviation of the accuracy, in the training stage, was 6.11% and 1.71%, respectively, for the ANN model and the SVM model; in the testing stage, these were 1.75% and 2.72%. In terms of accuracy, the SVM model performed better than the ANN model (5.29% for training stage, 13.07% for testing stage). In terms of stability (according to accuracy) in the training stage, the SVM model was better than the ANN model (4.40%); however, in the testing stage, the ANN model was better than the SVM model (0.97%). From the aspect of mean accuracy, the SVM model was superior to the ANN model in both prediction accuracy and stability.
For the other four statistical parameters, the stability of the two models declined in the testing stage. In particular, the stability of the positive predictive value and negative positive value of the ANN model decreased greatly, which means the accuracy of the ANN model in predicting landslide units and non-landslide units is very unstable. The mean values of the four statistical parameters of the two models declined, but the ANN model showed a large drop (more than 10%). From the aspect of these four statistical parameters, the SVM model could also be considered as the optimal model.

5.3. Comparison with Other Models

In the evaluation of related geological disasters in Helong city, we used the information content method (ICM) and analytic hierarchy process (AHP) to evaluate the landslide susceptibility of the whole of Helong city based on both the grid unit and slope unit; the results are shown in Table 4 First, from the perspective of different mapping units, it can be seen that when ICM and AHP use grid units and slope units, respectively, for the landslide susceptibility evaluation of Helong city, their prediction accuracy is quite different. The prediction accuracy of the slope unit is obviously higher than that of the grid unit. This shows that it is reasonable to use the slope unit as the mapping unit of landslide susceptibility in this paper. Furthermore, the prediction accuracy of the two models established in this paper is higher than that of ICM and AHP models. In particular, the prediction accuracy of SVM model is the highest. This is because landslide susceptibility mapping is a dichotomous problem, and the SVM model is better for making predictions in such problems.

5.4. Landslide Suceptibility Map analysis

The landslide susceptibility map should meet the following two requirements [37]:
(1) The landslide points should be distributed in the areas with high susceptibility as much as possible. The purpose of this is to evaluate the accuracy of grading the sensitivity of landslides.
(2) In the landslide susceptibility map, the points predicted to be of high susceptibility should account for as low a proportion as possible. The purpose of this is to reduce the redundancy of high landslide susceptibility prediction and improve the hit ratio of susceptibility assessment.
Table 3 shows that the very high and high susceptibility areas for SVM model had a combined area of 279.03 km2, accounting for 18.90% of the total study area. Regarding the landslide occurrence, the very high and high susceptibility areas had 66 landslides, accounting for 79.52% of the total landslides. For the ANN model, the very high and high susceptibility areas had a combined area of 444.19 km2, accounting for 30.01% of the total study area. Regarding the landslide occurrence, the very high and high susceptibility areas had 61 landslides, accounting for 70.49% of the total landslides. This shows that the landslide susceptibility map produced by the SVM model in this paper is more reasonable than that of the ANN model.
According to Figure 11, the very high, high, and moderate susceptibility areas were mainly distributed along rivers. The reason for this is that the erosion caused by a river increases the slope angle along the river. To adapt to the erosion of the river, landslides and other disasters often occur along rivers, thus reducing the slope angle adjacent to the river. Furthermore, for slopes located close to a river, the slope foot is soaked by the river water, which reduces the strength of the rock and soil mass inside, thus leading to landslides. The low susceptibility area was mainly distributed in the northwest of the study area. In this area, the vegetation coverage is high and human engineering activities are weak; thus, the geological environment of this area is relatively stable, and landslides are less common.

6. Conclusions

In this study, we produced a landslide susceptibility map of the area to the southeast of Helong city. A total of 83 landslides were mapped in the study area through remote sensing interpretation and field investigation. The slope unit divided by the curvature watershed method was selected as the mapping unit. Based on field investigations and previous studies, three groups of influencing factors—lithological factors, topographic factors, and geological environment factors (a total of ten influencing factors: lithology, elevation, slope angle, slope aspect, topographic relief, curvature, land use, rainfall, distance to river, and distance to faults)—were selected to establish landslide susceptibility models. The ANN and SVM methods were used to build the models. Five-fold cross-validation, the ROC curve, and statistical parameters were used to optimize the models. Finally, the landslide susceptibility map was classified into five classes—very high, high, moderate, low, and very low.
According to the model fitting results, the SVM model was the optimal model for our purposes. For the five classes—very high, high, moderate, low, and very low—the areas of the grades from the SVM model were 127.43, 151.60, 198.77, 491.19, and 506.91 km2, respectively. For landslide occurrence, the number of landslides in the five susceptibility classes were 52, 14, 8, 4, and 5, respectively; therefore, the very high and high susceptibility areas included 79.52% of the total landslides. This indicates that the landslide susceptibility map produced in this paper is reasonable.
In conclusion, the following inferences were obtained:
(1) Because the slope units are more closely related to the real topographic and geomorphic features, the obtained landslide susceptibility map is more reasonable. Thus, it is suggested that the slope unit should be used as the mapping unit in related work.
(2) The SVM model needs less data and is more suitable for binary classification, so it can be given priority in the study of landslide susceptibility mapping.

Author Contributions

C.Y. contributed to data analysis and manuscript writing. J.C. proposed the main structure of this study. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundations of China grant number [U1702241].

Acknowledgments

Thanks to anonymous reviewers for their valuable feedback on the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Description of the Landslides.
Table A1. Description of the Landslides.
No.LithologySlope
Angle
Slope
Aspect
Slope
Height
Slope
Shape
MicroreliefLandslide
Scale
Failure
Mode
1Granodiorite601986ConvexSteep SlopeSmallPull-Type
2Granodiorite7916217ConvexSteep CoastMiddlePull-Type
3Granodiorite6316224ConvexSteep CoastSmallToppling
4Monzonitic Granite655020ConcaveSteep CoastMiddlePull-Type
5Monzonitic Granite714511ConvexSteep CoastSmallPull-Type
6Granodiorite5221524ConvexSteep SlopeSmallSliding
7Granodiorite5711015ConcaveSteep SlopeMiddlePull-Type
8Monzonitic Granite6914817ConvexSteep CoastSmallPull-Type
9Monzonitic Granite6222117ConvexSteep CoastSmallSliding
10Dioritic Porphyrite5227512ConvexSteep SlopeSmallPull-Type
11Basalt8117850ConvexSteep CoastMiddleToppling
12Basalt8620242ConvexSteep CoastMiddlePull-Type
13Basalt7510517ConvexSteep CoastSmallSliding
14Basalt5812823ConcaveSteep SlopeMiddlePull-Type
15Basalt8717832ConvexSteep CoastMiddlePull-Type
16Basalt8213020ConcaveSteep CoastMiddlePull-Type
17Basalt6817818ConvexSteep CoastSmallSliding
18Basalt4027025StraightSteep SlopeSmallPull-Type
19Basalt8425312ConvexSteep CoastSmallPull-Type
20Basalt8628025ConvexSteep CoastMiddlePull-Type
21Basalt8418320ConvexSteep CoastMiddlePull-Type
22Basalt8121012ConcaveSteep CoastSmallPull-Type
23Granodiorite7119323ConvexSteep CoastSmallPull-Type
24Granodiorite3719549ConvexSteep SlopeSmallPull-Type
25Granodiorite4019643ConvexSteep SlopeSmallToppling
26Granodiorite3815425ConcaveSteep SlopeSmallPull-Type
27Granodiorite422689StraightSteep SlopeMiddlePull-Type
28Granodiorite552548ConcaveSteep SlopeSmallPull-Type
29Monzonitic Granite4319012StraightSteep SlopeSmallPull-Type
30Granite4123011ConvexSteep SlopeMiddlePull-Type
31Monzonitic Granite6213010ConvexSteep CoastSmallPull-Type
32Basalt8725012ConcaveSteep CoastSmallPull-Type
33Diorite7115510ConvexSteep CoastSmallPull-Type
34Granodiorite421368ConvexSteep SlopeSmallPull-Type
35Granodiorite551525ConvexSteep SlopeSmallPull-Type
36Basalt5013514StraightSteep SlopeSmallPull-type
37Granodiorite71909ConvexSteep CoastSmallToppling
38Monzonitic Granite3518615ConvexSteep SlopeSmallPull-Type
39Monzonitic Granite3613530StraightSteep SlopeSmallPull-Type
40Diorite3717450StraightSteep SlopeSmallToppling
41Granodiorite4814816ConvexSteep SlopeSmallPull-Type
42Granodiorite8512812ConvexSteep CoastSmallToppling
43Granodiorite7722415ConvexSteep CoastSmallToppling
44Granodiorite6718877ConvexSteep CoastMiddleSliding
45Monzonitic Granite5720461ConcaveSteep SlopeMiddleStaggered Breaking
46Monzonitic Granite8122825ConvexSteep CoastMiddleToppling
47Monzonitic Granite7225052ConvexSteep CoastMiddleToppling
48Monzonitic Granite7218531ConcaveSteep CoastMiddleToppling
49Monzonitic Granite76190142ConvexSteep CoastMiddleToppling
50Monzonitic Granite7412013ConvexSteep CoastSmallToppling
51Monzonitic Granite69204108ConvexSteep CoastLargePull-Splitting
52Diorite78172160ConvexSteep CoastLargeToppling
53Diorite62134224ConvexSteep CoastLargeToppling
54Granodiorite682423ConvexSteep CoastMiddleToppling
55Granodiorite65150136ConcaveSteep CoastLargeSliding
56Granodiorite76162228ConvexSteep CoastLargePull-Splitting
57Monzonitic Granite70155151ConcaveSteep CoastLargeStaggered Breaking
58Monzonitic Granite5115532ConvexSteep SlopeSmallToppling
59Granodiorite707232StraightSteep CoastSmallPull-Splitting
60Basalt7212546ConvexSteep CoastMiddleToppling
61Basalt5820045ConvexSteep SlopeMiddleToppling
62Basalt8321828ConvexSteep CoastSmallToppling
63Basalt6015289ConcaveSteep SlopeMiddleSliding
64Basalt521057ConvexSteep SlopeSmallPull-Splitting
65Granodiorite6013023ConvexSteep SlopeSmallPull-Splitting
66Andesite3822065ConvexSteep SlopeMiddlePull-Splitting
67Andesite5024555ConcaveSteep SlopeMiddlePull-Splitting
68Andesite7013515StraightSteep CoastSmallSliding
69Monzonitic Granite4421512ConcaveSteep SlopeSmallToppling
70Granodiorite7019024StraightSteep CoastSmallToppling
71Granite5317014ConvexSteep CoastSmallPull-Splitting
72Quartzite7313040ConvexSteep CoastSmallPull-Splitting
73Quartzite7218442ConcaveSteep CoastSmallStaggered Breaking
74Quartzite5215450ConvexSteep SlopeSmallStaggered Breaking
75Quartzite7019245ConvexSteep CoastMiddleStaggered Breaking
76Monzonitic Granite5819555ConvexSteep SlopeMiddleStaggered Breaking
77Granodiorite80568StraightSteep CoastSmallPull-Splitting
78Monzonitic Granite7117075ConvexSteep CoastSmallToppling
79Monzonitic Granite7513840ConvexSteep CoastSmallToppling
80Monzonitic Granite518515StraightSteep SlopeSmallToppling
81Granodiorite681757StraightSteep CoastSmallToppling
82Basalt621847ConvexSteep CoastSmallToppling
83Granodiorite4016019ConvexSteep SlopeSmallPull-Splitting
Note: Landslide scale: (a) Small; Landslide volume <1 × 104 m3; (b) Middle: 100 × 104 m3 < Landslide volume <104 m3.

References

  1. Cao, C.; Wang, Q.; Chen, J.; Ruan, Y.; Zheng, L.; Song, S.; Niu, C. Landslide susceptibility mapping in vertical distribution law of precipitation area: Case of the xulong hydropower station reservoir, southwestern China. Water 2016, 8, 270. [Google Scholar] [CrossRef] [Green Version]
  2. Sun, X.; Chen, J.; Bao, Y.; Han, X.; Zhan, J.; Peng, W. Landslide susceptibility mapping using logistic regression analysis along the jinsha river and its tributaries close to derong and deqin county, southwestern China. ISPRS Int. J. Geo-Inf. 2018, 7, 438. [Google Scholar] [CrossRef] [Green Version]
  3. Zhan, J.; Chen, J.; Zhang, W.; Han, X.; Sun, X.; Bao, Y. Mass movements along a rapidly uplifting river valley: An example from the upper Jinsha River, southeast margin of the Tibetan Plateau. Environ. Earth Sci. 2018, 77, 634. [Google Scholar] [CrossRef]
  4. Wang, F.; Xu, P.; Wang, C.; Wang, N.; Jiang, N. Application of a GIS-Based slope Unit Method for landslide susceptibility mapping along the Longzi River, Southeastern Tibetan Plateau, China. ISPRS Int. J. Geo-Inf. 2017, 6, 172. [Google Scholar] [CrossRef] [Green Version]
  5. Hearn, G.J.; Hart, A.B. Landslide susceptibility mapping: A practitioner’s view. Bull. Eng. Geol. Environ. 2019, 78, 5811–5826. [Google Scholar] [CrossRef]
  6. 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]
  7. Guzzetti, F.; Mondini, A.C.; Cardinali, M.; Fiorucci, F.; Santangelo, M.; Chang, K.T. Landslide inventory maps: New tools for an old problem. Earth-Sci. Rev. 2012, 112, 42–66. [Google Scholar] [CrossRef] [Green Version]
  8. Yang, X.J.; Chen, L.D. Geoinformation. Using multi-temporal remote sensor imagery to detect earthquake-triggered landslides. Int. J. Appl. Earth Obs. Geoinf. 2010, 12, 487–495. [Google Scholar] [CrossRef]
  9. Mulas, M.; Corsini, A.; Cuozzo, G.; Callegari, M.; Thiebes, B.; Mair, V. Quantitative Monitoring of Surface Movements on Active Landslides by Multi-Temporal, High-Resolution X-Band SAR Amplitude Information: Preliminary Results; CRC Press: Boca Raton, FL, USA, 2016; pp. 1511–1516. [Google Scholar]
  10. Schlogel, R.; Marchesini, I.; Alvioli, M.; Reichenbach, P.; Rossi, M.; Malet, J.P. Optimizing landslide susceptibility zonation: Effects of DEM spatial resolution and slope unit delineation on logistic regression models. Geomorphology 2018, 301, 10–20. [Google Scholar] [CrossRef]
  11. Rotigliano, E.; Cappadonia, C.; Conoscenti, C.; Costanzo, D.; Agnesi, V. Slope units-based flow susceptibility model: Using validation tests to select controlling factors. Nat. Hazards 2012, 61, 143–153. [Google Scholar] [CrossRef]
  12. Alvioli, M.; Marchesini, I.; Reichenbach, P.; Rossi, M.; Ardizzone, F.; Fiorucci, F.; Guzzetti, F. Automatic delineation of geomorphological slope units with r.slopeunits v1.0 and their optimization for landslide susceptibility modeling. Geosci. Model Dev. 2016, 9, 3975–3991. [Google Scholar] [CrossRef] [Green Version]
  13. Pourghasemi, H.R.; Rossi, M. Landslide susceptibility modeling in a landslide prone area in Mazandarn Province, north of Iran: A comparison between GLM, GAM, MARS, and M-AHP methods. Theor. Appl. Climatol. 2017, 130, 609–633. [Google Scholar] [CrossRef]
  14. Lagomarsino, D.; Tofani, V.; Segoni, S.; Catani, F.; Casagli, N. A Tool for Classification and Regression Using Random Forest Methodology: Applications to Landslide Susceptibility Mapping and Soil Thickness Modeling. Environ. Model. Assess. 2017, 22, 201–214. [Google Scholar] [CrossRef]
  15. Dou, J.; Bui, D.T.; Yunus, A.P.; Jia, K.; Song, X.; Revhaug, I.; Xia, H.; Zhu, Z. Optimization of Causative Factors for Landslide Susceptibility Evaluation Using Remote Sensing and GIS Data in Parts of Niigata, Japan. PLoS ONE 2015, 10, e0133262. [Google Scholar] [CrossRef] [Green Version]
  16. Frattini, P.; Crosta, G.; Carrara, A. Techniques for evaluating the performance of landslide susceptibility models. Eng. Geol. 2010, 111, 62–72. [Google Scholar] [CrossRef]
  17. Chen, W.; Li, W.; Chai, H.; Hou, E.; Li, X.; Ding, X. GIS-based landslide susceptibility mapping using analytical hierarchy process (AHP) and certainty factor (CF) models for the Baozhong region of Baoji City, China. Environ. Earth Sci. 2016, 75, 63. [Google Scholar] [CrossRef]
  18. Cao, C.; Xu, P.; Wang, Y.; Chen, J.; Zheng, L.; Niu, C. Flash Flood Hazard Susceptibility Mapping Using Frequency Ratio and Statistical Index Methods in Coalmine Subsidence Areas. Sustainability 2016, 8, 948. [Google Scholar] [CrossRef] [Green Version]
  19. Xiao, T.; Segoni, S.; Chen, L.; Yin, K.; Casagli, N. A step beyond landslide susceptibility maps: A simple method to investigate and explain the different outcomes obtained by different approaches. Landslides 2020, 17, 627–640. [Google Scholar] [CrossRef] [Green Version]
  20. Bai, S.-B.; Wang, J.; Lue, G.-N.; Zhou, P.-G.; Hou, S.-S.; Xu, S.-N. GIS-based logistic regression for landslide susceptibility mapping of the Zhongxian segment in the Three Gorges area, China. Geomorphology 2010, 115, 23–31. [Google Scholar] [CrossRef]
  21. Sun, X.; Chen, J.; Han, X.; Bao, Y.; Zhan, J.; Peng, W. Application of a GIS-based slope unit method for landslide susceptibility mapping along the rapidly uplifting section of the upper Jinsha River, South-Western China. Bull. Eng. Geol. Environ. 2019. [Google Scholar] [CrossRef]
  22. Segoni, S.; Pappafico, G.; Luti, T.; Catani, F. Landslide susceptibility assessment in complex geological settings: Sensitivity to geological information and insights on its parameterization. Landslides 2020. [Google Scholar] [CrossRef] [Green Version]
  23. Ren, D.; Zhou, D.; Liu, D.; Dong, F.; Ma, S.; Huang, H. Formation mechanism of the Upper Triassic Yanchang Formation tight sandstone reservoir in Ordos Basin-Take Chang 6 reservoir in Jiyuan oil field as an example. J. Pet. Sci. Eng. 2019, 178, 497–505. [Google Scholar] [CrossRef]
  24. Ayalew, L.; Yamagishi, H. The application of GIS-based logistic regression for landslide susceptibility mapping in the Kakuda-Yahiko Mountains, Central Japan. Geomorphology 2005, 65, 15–31. [Google Scholar] [CrossRef]
  25. Lee, S.; Min, K. Statistical analysis of landslide susceptibility at Yongin, Korea. Environ. Geol. 2001, 40, 1095–1113. [Google Scholar] [CrossRef]
  26. Conforti, M.; Pascale, S.; Robustelli, G.; Sdao, F. Evaluation of prediction capability of the artificial neural networks for mapping landslide susceptibility in the Turbolo River catchment (northern Calabria, Italy). Catena 2014, 113, 236–250. [Google Scholar] [CrossRef]
  27. Oh, H.-J.; Pradhan, B. Application of a neuro-fuzzy model to landslide-susceptibility mapping for shallow landslides in a tropical hilly area. Comput. Geosci. 2011, 37, 1264–1276. [Google Scholar] [CrossRef]
  28. Park, S.; Choi, C.; Kim, B.; Kim, J. Landslide susceptibility mapping using frequency ratio, analytic hierarchy process, logistic regression, and artificial neural network methods at the Inje area, Korea. Environ. Earth Sci. 2013, 68, 1443–1464. [Google Scholar] [CrossRef]
  29. Hong, H.; Pradhan, B.; Xu, C.; Tien Bui, D. 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]
  30. Binh Thai, P.; Dieu Tien, B.; Dholakia, M.B.; Prakash, I.; Ha Viet, P.; Mehmood, K.; Hung Quoc, L. A novel ensemble classifier of rotation forest and Naive Bayer for landslide susceptibility assessment at the Luc Yen district, Yen Bai Province (Viet Nam) using GIS. Geomat. Nat. Hazards Risk 2017, 8, 649–671. [Google Scholar] [CrossRef] [Green Version]
  31. Chen, W.; Xie, X.; Wang, J.; Pradhan, B.; Hong, H.; Bui, D.T.; Duan, Z.; Ma, J. A comparative study of logistic model tree, random forest, and classification and regression tree models for spatial prediction of landslide susceptibility. Catena 2017, 151, 147–160. [Google Scholar] [CrossRef] [Green Version]
  32. Dieu Tien, B.; Tran Anh, T.; Klempe, H.; Pradhan, B.; Revhaug, I. Spatial prediction models for shallow landslide hazards: A comparative assessment of the efficacy of support vector machines, artificial neural networks, kernel logistic regression, and logistic model tree. Landslides 2016, 13, 361–378. [Google Scholar] [CrossRef]
  33. Kumar, D.; Thakur, M.; Dubey, C.S.; Shukla, D.P. Landslide Susceptibility Mapping & Prediction using Support Vector Machine for Mandakini River Basin, Garhwal Himalaya, India. Geomorphology 2017, 295, 115–125. [Google Scholar]
  34. Chen, W.; Pourghasemi, H.R.; Kornejady, A.; Zhang, N. Landslide spatial modeling: Introducing new ensembles of ANN, MaxEnt, and SVM machine learning techniques. Geoderma 2017, 305, 314–327. [Google Scholar] [CrossRef]
  35. Yesilnacar, E.; Topal, T. Landslide susceptibility mapping: A comparison of logistic regression and neural networks methods in a medium scale study, Hendek region (Turkey). Eng. Geol. 2005, 79, 251–266. [Google Scholar] [CrossRef]
  36. Binh Thai, P.; Dieu Tien, B.; Prakash, I.; Dholakia, M.B. Hybrid integration of Multilayer Perceptron Neural Networks and machine learning ensembles for landslide susceptibility assessment at Himalayan area (India) using GIS. Catena 2017, 149, 52–63. [Google Scholar] [CrossRef]
  37. Su, Q.; Zhang, J.; Zhao, S.; Wang, L.; Liu, J.; Guo, J. Comparative Assessment of Three Nonlinear Approaches for Landslide Susceptibility Mapping in a Coal Mine Area. ISPRS Int. J. Geo-Inf. 2017, 6, 228. [Google Scholar] [CrossRef] [Green Version]
  38. Aditian, A.; Kubota, T.; Shinohara, Y. Comparison of GIS-based landslide susceptibility models using frequency ratio, logistic regression, and artificial neural network in a tertiary region of Ambon, Indonesia. Geomorphology 2018, 318, 101–111. [Google Scholar] [CrossRef]
  39. Jiang, P.; Zeng, Z.; Chen, J.; Huang, T. Generalized Regression Neural Networks with K-Fold Cross-Validation for Displacement of Landslide Forecasting. In Advances in Neural Networks—Isnn 2014; Zeng, Z., Li, Y., King, I., Eds.; Springer: Cham, Switzerland, 2014; Volume 8866, pp. 533–541. [Google Scholar]
  40. Jiang, T.; Lei, P.; Qin, Q. An Application of SVM-Based Classification in Landslide Stability. Intell. Autom. Soft Comput. 2016, 22, 267–271. [Google Scholar] [CrossRef]
  41. Li, X.Z.; Kong, J.M.J.N.H.; Sciences, E.S. Application of GA-SVM method with parameter optimization for landslide development prediction. Nat. Hazards Earth Syst. Sci. 2014, 14, 5295–5322. [Google Scholar] [CrossRef] [Green Version]
  42. Pradhan, B.; Sameen, M.I. Manifestation of SVM-Based Rectified Linear Unit (ReLU) Kernel Function in Landslide Modelling; Springer: Singapore, 2018. [Google Scholar]
  43. San, B.T. An evaluation of SVM using polygon-based random sampling in landslide susceptibility mapping: The Candir catchment area (western Antalya, Turkey). Int. J. Appl. Earth Obs. Geoinf. 2014, 26, 399–412. [Google Scholar] [CrossRef]
  44. Ballabio, C.; Sterlacchini, S. Support Vector Machines for Landslide Susceptibility Mapping: The Staffora River Basin Case Study, Italy. Math. Geosci. 2012, 44, 47–70. [Google Scholar] [CrossRef]
  45. Chen, W.; Pourghasemi, H.R.; Panahi, M.; Kornejady, A.; Wang, J.; Xie, X.; Cao, S. Spatial prediction of landslide susceptibility using an adaptive neuro-fuzzy inference system combined with frequency ratio, generalized additive model, and support vector machine techniques. Geomorphology 2017, 297, 69–85. [Google Scholar] [CrossRef]
  46. Zhang, X.; Wang, Q.; Huo, Z.; Yu, T.; Wang, G.; Liu, T.; Wang, W. Prediction of Frost-Heaving Behavior of Saline Soil in Western Jilin Province, China, by Neural Network Methods. Math. Probl. Eng. 2017. [Google Scholar] [CrossRef] [Green Version]
  47. Hong, H.; Pradhan, B.; Jebur, M.N.; Bui, D.T.; Xu, C.; Akgun, A. Spatial prediction of landslide hazard at the Luxi area (China) using support vector machines. Environ. Earth Sci. 2016, 75. [Google Scholar] [CrossRef]
  48. Huang, F.; Yin, K.; Huang, J.; Gui, L.; Wang, P. Landslide susceptibility mapping based on self-organizing-map network and extreme learning machine. Eng. Geol. 2017, 223, 11–22. [Google Scholar] [CrossRef]
Figure 1. Geographical position and landslide inventory of the study area. (a,b): geographical position of the study area; (c) landslide inventory of the study area; (d) geological map of the study area; (e,f) typical landslides of the study area.
Figure 1. Geographical position and landslide inventory of the study area. (a,b): geographical position of the study area; (c) landslide inventory of the study area; (d) geological map of the study area; (e,f) typical landslides of the study area.
Symmetry 12 01047 g001
Figure 2. Flowchart of this study.
Figure 2. Flowchart of this study.
Symmetry 12 01047 g002
Figure 3. Construction process of the slope units based on the curvature watershed method.
Figure 3. Construction process of the slope units based on the curvature watershed method.
Symmetry 12 01047 g003
Figure 4. (a) Elevation map; (b) slope angle map.
Figure 4. (a) Elevation map; (b) slope angle map.
Symmetry 12 01047 g004
Figure 5. Influencing factor maps of the study area: (a) lithology (downloaded using 91 Weitu software), (b) elevation, (c) slope angle, and (d) slope aspect.
Figure 5. Influencing factor maps of the study area: (a) lithology (downloaded using 91 Weitu software), (b) elevation, (c) slope angle, and (d) slope aspect.
Symmetry 12 01047 g005
Figure 6. Influencing factors maps of the study area: (a) topographic relief, (b) curvature, (c) land-use, and (d) rainfall.
Figure 6. Influencing factors maps of the study area: (a) topographic relief, (b) curvature, (c) land-use, and (d) rainfall.
Symmetry 12 01047 g006
Figure 7. Influencing factors maps of the study area: (a) distance to river and (b) distance to fault.
Figure 7. Influencing factors maps of the study area: (a) distance to river and (b) distance to fault.
Symmetry 12 01047 g007
Figure 8. Schematic diagram of a simple artificial neural network (ANN) model.
Figure 8. Schematic diagram of a simple artificial neural network (ANN) model.
Symmetry 12 01047 g008
Figure 9. Division of slope units of the study area.
Figure 9. Division of slope units of the study area.
Symmetry 12 01047 g009
Figure 10. Area under the curve (AUC) values of model evaluation parameters when the hidden neurons of ANN model varied.
Figure 10. Area under the curve (AUC) values of model evaluation parameters when the hidden neurons of ANN model varied.
Symmetry 12 01047 g010
Figure 11. Landslide susceptibility map. (a) ANN model; and (b) support vector machine (SVM) model.
Figure 11. Landslide susceptibility map. (a) ANN model; and (b) support vector machine (SVM) model.
Symmetry 12 01047 g011
Figure 12. Statistics of morphological characteristics of slope units. (a) Slope unit area distribution diagram; (b) slope unit shape index distribution diagram.
Figure 12. Statistics of morphological characteristics of slope units. (a) Slope unit area distribution diagram; (b) slope unit shape index distribution diagram.
Symmetry 12 01047 g012
Table 1. Description of the Lithology.
Table 1. Description of the Lithology.
SystemCodeLithology
QuaternaryQAlluvial–Diluvial, Gravel, Sub-Sandy Soil, Sub-Clay, and Basalt
NeogeneNSandstone, Conglomerate, and Siltstone with Basalt
CretaceousKSandstone, Conglomerate, Siltstone with Limestone, and Oil Limestone
JurassicJAndesite and Tuff
Middle ProterozoicPtMarble
New ArcheanArBlack Cloud Amphibolic Granulite and Granulite
Table 2. Performance of the Two Slope-Unit-Fitted Landslide Susceptibility Models (%).
Table 2. Performance of the Two Slope-Unit-Fitted Landslide Susceptibility Models (%).
StageMethodStatistical IndexK = 1K = 2K = 3K = 4K = 5MeanStandard Deviation
TrainingANNAUC88.2091.1090.3084.7088.7088.602.48
AC86.5789.3985.0773.4881.8283.276.11
SE91.5390.6383.1072.4682.8184.117.68
SP82.6788.2487.3074.6080.8882.745.49
PP80.6087.8888.0675.7680.3082.525.33
NP92.5490.9182.0971.2183.3384.028.49
SVMAUC92.7093.4092.3093.2094.5093.220.83
AC87.3189.3986.5788.6490.9188.561.71
SE89.0691.9487.6990.4893.5590.542.31
SP85.7187.1485.5186.9688.5786.781.24
PP85.0786.3685.0786.3687.8886.151.16
NP89.5592.4288.0690.9193.9490.982.31
TestingANNAUC83.2087.0082.1088.2082.1084.522.88
AC71.8870.5968.7571.4373.5371.231.75
SE73.3381.8265.0066.6772.2271.816.62
SP70.5965.2275.0078.5775.0072.885.13
PP68.7552.9481.2582.3576.4772.3512.10
NP75.0088.2456.2564.7170.5970.9611.94
SVMAUC88.7089.6091.2091.3087.9089.741.50
AC81.2588.2484.3885.2982.3584.302.72
SE81.2593.3382.3583.3382.3584.524.98
SP81.2584.2186.6787.5082.3584.402.69
PP81.2582.3587.5088.2482.3584.343.26
NP81.2594.1281.2582.3582.3584.265.54
Note: K is the number of cross validations.
Table 3. Statistical Results of the Landslide Susceptibility Map.
Table 3. Statistical Results of the Landslide Susceptibility Map.
ModelSusceptibilityLandslide OccurredTotal Study Area
CountRatioArea (km2)Ratio
ANNVery Low67.23%297.9520.19%
Low67.23%310.4421.03%
Moderate1012.05%423.3328.68%
High1821.69%297.9520.19%
Very High4351.81%146.249.91%
SVMVery Low56.02%506.9134.35%
Low44.82%491.1933.28%
Moderate89.64%198.7713.47%
High1416.87%151.6010.27%
Very High5262.65%127.438.63%
Table 4. Modeling Fitting Results of the Study Area.
Table 4. Modeling Fitting Results of the Study Area.
Mapping UnitsMethodPrediction Accuracy (Mean)
Slope UnitsANNTraining89.72%
Validating88.08%
Slope UnitsSVMTraining90.72%
Validating88.96%
Grid UnitsICM-83.42%
Grid UnitsAHP-70.93%
Slope UnitsICM-87.11%
Slope UnitsAHP-80.54%

Share and Cite

MDPI and ACS Style

Yu, C.; Chen, J. Landslide Susceptibility Mapping Using the Slope Unit for Southeastern Helong City, Jilin Province, China: A Comparison of ANN and SVM. Symmetry 2020, 12, 1047. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12061047

AMA Style

Yu C, Chen J. Landslide Susceptibility Mapping Using the Slope Unit for Southeastern Helong City, Jilin Province, China: A Comparison of ANN and SVM. Symmetry. 2020; 12(6):1047. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12061047

Chicago/Turabian Style

Yu, Chenglong, and Jianping Chen. 2020. "Landslide Susceptibility Mapping Using the Slope Unit for Southeastern Helong City, Jilin Province, China: A Comparison of ANN and SVM" Symmetry 12, no. 6: 1047. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12061047

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