Next Article in Journal
Concept and Implementation of Solutions Improving Water Relations in the Area of the Flooded Opencast Lignite Mine Kazimierz Północ in the East Wielkopolska Region (Central-West Poland)
Next Article in Special Issue
Mechanism Study of Differential Permeability Evolution and Microscopic Pore Characteristics of Soft Clay under Saturated Seepage: A Case Study in Chongming East Shoal
Previous Article in Journal
Design and Experimental Studies on a Single Slope Solar Still for Water Desalination
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of Machine Learning and Traditional Statistical Methods in Debris Flow Susceptibility Assessment: A Case Study of Changping District, Beijing

1
College of Construction Engineering, Jilin University, Changchun 130026, China
2
Department of Earth Sciences and Engineering, Taiyuan University of Technology, Taiyuan 030024, China
3
Key Laboratory of Shale Gas and Geoengineering, Institute of Geology and Geophysics, University of Chinese Academy of Sciences, Beijing 100029, China
*
Author to whom correspondence should be addressed.
Water 2023, 15(4), 705; https://doi.org/10.3390/w15040705
Submission received: 12 January 2023 / Revised: 5 February 2023 / Accepted: 7 February 2023 / Published: 10 February 2023
(This article belongs to the Special Issue Effects of Groundwater and Surface Water on the Natural Geo-Hazards)

Abstract

:
As a common geological hazard, debris flow is widely distributed around the world. Meanwhile, due to the influence of many factors such as geology, geomorphology and climate, the occurrence frequency and main inducing factors are different in different places. Therefore, the evaluation of debris flow sensitivity can provide a very important theoretical basis for disaster prevention and control. In this research, 43 debris flow gullies in Changping District, Beijing were cataloged and studied through field surveys and the 3S technology (GIS (Geography Information Systems), GPS (Global Positioning Systems), RS (Remote Sensing)). Eleven factors, including elevation, slope, plane curvature, profile curvature, roundness, geomorphic information entropy, TWI, SPI, TCI, NDVI and rainfall, were selected to establish a comprehensive evaluation index system. The watershed unit is directly related to the development and activities of debris flow, which can fully reflect the geomorphic and geological environment of debris flow. Therefore, the watershed unit was selected as the basic mapping unit to establish four evaluation models, namely ACA–PCA–FR (Analytic Hierarchy Process–Principal Component Analysis–Frequency Ratio), FR (Frequency Ratio), SVM (Support Vector Machines) and LR (Logistic Regression). In other words, this research evaluates debris flow susceptibility by comparingit with two traditional weight methods (ACA–PCA–FR and FR) and two machine learning methods (SVM and LR). The results show that the SVM evaluation model is superior to the other three models, and thevalueofthe area under the receiver-operating characteristic curve (AUC) is 0.889 from the receiver operating characteristic curve (ROC). It verifies that the SVM model has strong adaptability to small sample data. The study was divided into five regions, which were very low, low, moderate, high and very high, accounting for 22.31%, 25.04%, 17.66%, 18.85% and 16.14% of the total study area, respectively, by SVM model. The results obtained in this researchagree with the actual survey results, and can provide theoretical help for disaster prevention and reduction projects.

1. Introduction

Debris flow is a special flood that flows along slopes or gullies with a large amount of sediment, debris and other solid materials mixed with water. Debris flow is also a sudden geological disaster that often occurs in mountainous areas, bringing huge losses to the people. The debris flow in Zhouqu County, Gansu Province, China in 2010 resulted in 1487 deaths, and the direct economic loss alone was up to 400 million Yuan, which was the most disastrous geological disaster in China in the past 20 years [1]. In recent years, the frequency of debris flow in China is high, which seriously threatens the safety of life and property. Zhang et al. [2] counted all 10,927 debris flow disasters in China from 2005 to 2015 and found that the peak occurrence period of debris flow disasters in China was from May to September each year, which was the same time as the rainy season in most parts of China.
The research on debris flow susceptibility has mainly gone through three stages [3]: (1) In 1976, the United Nations entrusted the International Federation of Engineering Geology to research debris flow susceptibility [4], and since then it has gradually become an important part of the evaluation and prevention of debris flow disaster. KOVACS et al. [5] used the qualitative evaluation method to analyze debris flow susceptibility in the 1980s, which provided ideas for debris flow susceptibility assessment. (2) With the successful application of mathematical statistical methods and the 3S technology [6] in the susceptibility analysis of debris flow, debris flow susceptibility assessment gradually developed towarda comprehensive quantitative evaluation. These methods transition from qualitative to semi-quantitative and quantitative, and have been applied for a long time with a wide range of applications. The commonly used methods include the certainty factor (CF) [7], the analytic hierarchy process (AHP) [8], the information content model (ICM) [9], the principal component analysis (PCA) [10] and the frequency ratio method (FR) [11]. Li et al. [12], based on the analytic hierarchy process and principal component analysis method for elevation, combined eight influencing factors, such as slope, to give weight, and then the results of debris flow susceptibility assessment in Mentougou District of Beijing are obtained by the information content model. Li et al. [13] compared debris flow susceptibility in Pinggu District of Beijing by using the certainty factor (CF), the Information value (IV) model and the logistic regression (LR) method, and the results show that the sensitivity result obtained by the logistic regression (LR) model is the most accurate. The calculation process of these methods is complex and the calculation time is long [14]. Moreover, in the case of many damage points and evaluation factors, the above methods cannot gain advantages in the calculation of complex and big data. (3) With the rapid development of computer technology, more intelligent algorithms were put forward and widely used [15]. Among them, the machine learning algorithms were very widely used, and many scholars had applied them to the study and prediction of debris flow susceptibility. The commonly used machine learning algorithms were: support vector machine (SVM) [16], logistic regression (LR) [17], random forests (RF) [18], decision tree (DT) [19] and artificial neural network (ANN) [20]. Qiu et al. [21] improved and obtained more accurate GA-SVM and CF-GA-SVM basedon SVM. By comparing with SVM and LR models, debris flow susceptibility in the Jalong corridor area was researched. The results showed that the SVM model has higher accuracy than the LR model in the case of small sample data, and the improved SVM model has the most accurate prediction result. Ke et al. [22] used LR, SVM, RF and enhanced regression tree (RT) models to research debris flow susceptibility in Sichuan Province, China. The results showed that the enhanced RT model was more accurate than other models in prediction. Compared with mathematical statistical methods, machine learning methods were faster to calculate, capable of self-organizing learning, and more accurate when predicting and evaluating debris flow susceptibility [23]. Ahmad et al. [24] used Logistic Regression (LR), Shannon Entropy (SE), Weights-of-Evidence (WoE), and Frequency Ratio (FR) models to assess the geohazards susceptibility along the Upper Indus Basin; the results showed that the Logistic Regression model performed the best of all the models. Pal et al. [25] selected an ensemble of Bayesian generalized linear model (BGLM), sparse partial least squares (SPLS), boosted tree (BT), and random forest (RF) algorithms to evaluate debris flow and landslide hazards in Iran. Ciccarese et al. [26] put forward Frequency Ratio [FR], Weight of Evidence [WOE] and Logistic Regression [LR] models to map the debris flow hazard in Italy. Vianello et al. [27] used the Rock Engineering System (RES) method for debris flow susceptibility mapping in the Western Italian Alps.
Although there are many debris flow susceptibility assessment models, how to choose the model suitable for the study area is the most important [28]. Based on field investigation, remote sensing interpretation and historical data, the catchment unit was selected as the basic mapping unit in this research. After reading a large amount of scholarly literature on debris flow susceptibility assessment, rainfall, elevation, slope and another 11 factors were chosen to establish the FR, ACA–PCA–FR, SVM and LR models, and finally the debris flow susceptibility mapping results were obtained. The results show that the SVM model has the best performance among the four models, and the AUC value of the SVM model is the largest, which is 0.889. The result verifies that the SVM model has strong adaptability to small simple data and is very suitable for the evaluation of geological disasters. Finally, according to the results of the SVM model, the study area is divided into five regions, which are very low, low, moderate, high and very high zones, and accounted for 22.31%, 25.04%, 17.66%, 18.85% and 16.14%, respectively, of the total study area. Combined with the actual investigation and relevant data, the susceptibility assessment results obtained in this paper are reliable and can provide better help for disaster prevention and reduction in this region.

2. Study Area

The study area is located in the mountainous area of Changping District, which is located in the northwest of Beijing. The southernmost part of Changping district is only 40 km away from downtown Beijing. Changping District gradually forms a gentle slope zone from northwest to southeast, with mountains and half mountains in the west and north, and plains in the southeast, as presented in Figure 1 [29]. Changping District has a warm-temperate, semi-humid continental monsoon climate. Due to climate influence, the distribution of rainfall time is uneven in the study area. The average annual rainfall is 584.2 mm and is concentrated in summer. The rainfall from June to August is 443.3 mm, accounting for 76% of the rainfall. The temperature and rainfall show similar change laws. The average temperature of Changping District for many years has been 11.5–11.8 °C. The coldest month is January, with an average temperature of −4.1 °C. The hottest month is July, with an average temperature of 25.7 °C. Dolomite, granite, andesite, sandstone and alluvial diluvium are mainly exposed in the study area, as presented in Figure 2. The strength of rock mass is affected by fault and weathering, and broken rock blocks are produced. According to relevant data records, a large amount of rainfall in a short periodis the main cause of debris flow outbreak in the study area, and geological and geomorphic conditions are the basic conditions for the debris flow outbreak.
According to historical debris flow gully data, 43 debris flow gullies were investigated in this study, with a maximum area of 5.79 km2 and a minimum area of 0.10 km2. In the field survey, the vegetation in the study area is relatively lush. There are loose rocks in the gully, mainly from the broken rocks on both sides of the debris flow gully. There are several collapses on the slopes on both sides. There are some villages scattered in the study area, and the population density is relatively low. However, a quarry was found in the survey, which needs attention. The basic situation of debris flow gullies during the field investigation is presented in Figure 3. The basic data of 43 debris flow gullies can be found in Table A1 of Appendix A.
According to incomplete statistics, from 1950 to 1999, there were more than 200 debris flows in 29 times in Beijing, which killed 515 people and destroyed more than 8200 houses, as presented in Figure 1. On average, a disastrous debris flow occurred every 1.8 years. From 21 to 22 July 2012, most regions in China suffered heavy rain, including Beijing and its surrounding areas, which suffered the strongest rainstorm and flood disaster in 61 years. According to survey data, the rainstorm caused 79 deaths, 10,660 houses collapsed, 1.602 million people were affected, and the economic loss was 11.64 billion yuan [30].

3. Materials and Methods

3.1. Modeling Flow Chart

The debris flow susceptibility assessment mainly includes four processes: debris flow extraction and cataloging, model building, susceptibility mapping and result verification [31]. Among them, the establishment of the debris flow susceptibility assessment model is very important, which is directly related to the scientific reliability of debris flow assessment results. According to the historical debris flow disaster point data, this paper investigated 43 known debris flow gullies in detail, based on field surveys and the 3S technique, and also used the Borderline-SMOTE algorithm to generate 43 non-debris flow trenches [32], and established four models, FR, ACA–PCA–FR, SVM and LR, respectively. Then, the susceptibility assessment result chart was obtained according to the established model. Finally, the model was compared and verified by receiver-operating characteristic curve (ROC) and area under curve (AUC) values. The modeling flow chart of this research is presented in Figure 4.

3.2. Mapping Unit

The previous research on debris flow susceptibility mostly takes the valley watershed as the natural catchment basin [33]. Considering that the formation area and circulation area of debris flow in the study area are not distinguished, and the distribution of material sources is relatively wide, the use of the catchment unit as the evaluation unit can flexibly reflect the regional hydrological characteristics and surface topography by adjusting the water collection, and the evaluation results cover the whole area, with a wider application range [34]. At present, in the study of debris flow, relevant scholars have selected catchment units and achieved good results [35]. Thus, this research selects the catchment unit as the basic mapping unit. The main steps to extract the catchment units in the study area are as follows: (1) Carrying out hydrological analysis on DEM dataand dividing the catchment units according to different catchment volumes; (2) Comparing the terrain conditions and combining with the field investigation, finding the most practical extraction threshold; (3) Through manual comparison and identification, adjacent units are merged to obtain the modified division results. In this study, 5000, 10,000, 15,000 and 20,000 were, respectively, set as the thresholds for water-catchment extraction. After comparing the results, it was found that the watershed units divided when the threshold was 10,000 were the most reasonable. Therefore, 10,000 was selected as the extraction threshold, and the study area was, finally, divided into 273 catchment units.

3.3. Determination of Causative Factors

The number of factors that affect the development of debris flow is large and complex, and the factors are interrelated and interact with each other [36]. The causative factors of debris flow in this paper are mainly based on the distribution and development status of debris flow disasters in the study area. After referring to a large number of susceptibility assessment documents, 11 impact factors are selected, including elevation (H), slope, plane curvature (Pl.Cv), profile curvature (Pr.Cv), roundness (Rd), geomorphic information entropy (GIE), TWI, SPI, TCI, NDVI and rainfall (Rf). The elevation data is from the digital elevation data (DEM) of 91 Weitu software; with a resolution of 7.3 m × 7.3 m, the elevation directly reflects the topographic relief changes in the study area, and to a certain extent, reflects the changes in valleys, vegetation and the state of deposits [37]. The slope is one of the important factors affecting the occurrence of debris flow. In general, the probability of debris flow occurring in the area with a steep terrain slope is greater than that in the area with a gentle terrain slope. Secondly, the terrain slope also affects the distribution, migration and accumulation of material sources on the slope, as well as the erosion of runoff on the slope [38]. The slope can be extracted from DEM. Pl.Cv and Pr.Cv represent the concavity and convexity of the terrain surface, which indirectly affect the development range of debris flow [39]. When the Pl.Cv value is positive, it indicates that the terrain surface is convex upwards, while negative values indicate that the terrain is concave downwards. When the value is 0, the surface is horizontal. Pr.Cv is opposite to the plane curvature. Pl.Cv and Pl.Cv can be extracted from DEM. Roundness can be used to reflect the shape characteristics of debris flow basins [40]. It refers to the ratio of the area of a basin to the area of a circle with the same circumference, which can be calculated by the following equation:
Rd = 4 π A L 2  
where Rd is the roundness of the watershed; A is the area of each catchment unit; and L is the week of each catchment unit. In 1899, Davis et al. [41] put forward the theory of the geomorphic cycle and believed that the evolution process of the surface could be divided into three stages: initial uplift stage, weathering erosion stage and stable stage. These three stages were defined as the young, middle-aged and old age of watershed evolution, respectively. Ai et al. [42] combined analogical information entropy and the Straler integral theory, and proposed the calculation method of geomorphic information entropy GIE:
S = G mean G min G max G min
GIE = SlnS − 1
where S is the Straler integral value. Gmean is the average elevation of the small watershed; and Gmax and Gmin are the highest and lowest elevation values of the small watershed, respectively. The value of GIE is the same as the Straler integral value, which can reflect the development stage of the watershed. The larger the value is, the more stable the basin is; otherwise, the more intense the erosion is. When GIE < 0.111, the watershed development stage is juvenile; when 0.111 < GIE < 0.400, the development stage of the watershed is the prime period; when GIE > 0.400, the watershed development stage is the old age stage. The topographic wetness index (TWI) is defined as the theoretical measurement of flow accumulation and soil moisture at a point in the watershed [43]. Topographic moisture indices are commonly used to quantify topographic control over hydrological processes and as indicators of soil conditions and sediment and material accumulation. Therefore, debris flow susceptibility can be estimated as a function of the response of topography to hydrology in a region. Beven and Kirkby et al. [44] (1979) proposed the following formula to calculate TWI:
TWI = ln ( A s tan s )
SPI = A s   ×   tan s
TCI = k   ×   lnA s
where As is the specific catchment area; S is slope Angle; k is curvature; and TWI, SPI and TCI are dimensionless exponents.The stream power index (SPI) reflects the erosion ability of water flow to soil [45]. The higher the SPI value, the more the runoff concentration may lead to soil erosion, and the greater the probability of debris flow. The topographic characteristic index (TCI) mainly represents the migration ability of the river, to reflect the possibility of debris flow [46]. Rainfall is the most important factor for the formation of debris flow in the study area [47]. According to the historical debris flow disaster data, the occurrence of debris flow in the study area is mainly caused by continuous heavy rainfall in a short period. Therefore, this paper obtained the 24 h maximum rainfall data of the study area through Beijing Changping hydrological stations. Then, the rainfall factor map was obtained by the Kriging interpolation method in the Arcgis10.5. The Normalized Difference Vegetation Index (NDVI) reflected the size of the vegetation cover density in the study area [48]. The larger the NDVI value, the better the vegetation on the surface, which can effectively preserve the soil and water, so the more unfavorable it is to the occurrence of debris flow. The NDVI data used in this paper were acquired from Landsat8 OLI images with a resolution of 30 m. The spatial distribution of impact factors in this study area is presented in Figure 5.

3.4. Methods

3.4.1. Analytic Hierarchy Process (AHP)

The analytic hierarchy process (AHP) carries out the multi-criteria decision analysis method of relative importance calculation by decomposing all relevant elements in the whole process of weighting [49]. At present, the AHP method has been widely used in geological hazard mapping. When the AHP method is used to evaluate the risk of geological disasters, the corresponding hierarchical model should be established first. Then, the importance of the two influencing factors is compared subjectively and a judgment matrix is established. The judgment matrix can be expressed as:
A = a i j = a 11 a 12 a 1 n a 21 a 22 a 2 n a n 1 a n 2 a n n
where A is the judgment matrix; aij is the result of the comparison between the importance of factor i and factor j, which has the following characteristics:
a i j = 1 a j i
The relative importance of each factor can be rated on a scale of 1 to 9, from less important to more important. Finally, the following formula is used to check the consistency:
C I = λ m a x n n 1
CR = CI/RI
where CR is the random consistency ratio, when CR < 0.1, the consistency test passes; CI is the consistency index; n is the order of the judgment matrix; λmax is the largest eigenvalue of the judgment matrix. RI is the average random consistency index, as shown in Table 1.

3.4.2. Frequency Ratio (FR)

Frequency ratio (FR) is a binary statistical method used to calculate the probability relationship between the dependent variable and the independent variable, and it is a simple and effective zoning model for the debris flow hazard. Secondly, the frequency ratio method can also be used for multi-classification spatial data prediction [50]. This method determines the correlation between the hazard degree of debris flow and different evaluation indices by calculating the FR value. When the FR value is greater than 1, it indicates that the correlation between the debris flow risk degree and the evaluation index is high; when the FR value is less than 1, it indicates that the correlation between the debris flow risk degree and the evaluation index is small. The FR value can be calculated according to the following formula:
F R = n / N m / M
where n is the number of debris flow in the subinterval of the evaluation index; N is the total number of debris flow in the study area; m is the number of grids occupied by each evaluation index subinterval; and M is the total number of grids in the study area.

3.4.3. Principal Component Analysis (PCA)

In the research of practical problems, people always hope to use the factors involved in the problem as much as possible, in order to avoid the problem of unreliable analysis results due to insufficient information. However, the more complicated the consideration problem becomes, the higher the requirement for calculation and modeling analysis. This can also lead to a large amount of redundant information leading to unreliable results, so before establishing a model for data preprocessing, people remove the correlation between sample characteristic dimensions. PCA is a common method of data preprocessing, and it is used to process high-dimensional data in many fields [51]. When processing high-dimensional data, it can be represented by fewer principal components, and the dimensionality of sample data can be reduced under the premise of little data loss, to greatly simplify the cost of data calculation and research.

3.4.4. Support Vector Machine (SVM)

Support Vector Machine (SVM) is based on statistical learning theory [52]. It has the following advantages: (1) low requirements of data volume; (2) strong generalization ability; (3) strong adaptability to high-dimensional samples; and (4) strong learning ability and fast convergence. By transforming each evaluation index from low-dimensional space to high-dimensional space, it realizes the linear segmentation of data. By this method, nonlinear problems in low-dimensional space can be analyzed and evaluated. SVM is widely used in the research of debris flow hazard zoning because of its low requirement for data volume [53,54,55,56,57,58,59].

3.4.5. Logistic Regression (LR)

The main idea of the logistic regression (LR) model is to determine the possibility of debris flow in the future after converting each factor into a logical variable [60,61,62,63,64,65]. Logistic regression uses the maximum likelihood to find the best-fit function. The simplified logistic regression method is shown as follows:
P = 1 ( 1 + e y )
where P is the sensitivity index of debris flow, whose value is between 0 and 1. The larger the p-value, the greater the possibility of debris flow. y is the weighted sum of each evaluation index, which can be expressed as follows:
y = c0 + c1x1 + c2x2…+ cjxj
where c0 is a constant term, c1, c2, …, cj is the logistic regression coefficient; x1, x2, …, xj is the value of each evaluation index.

4. Result Analysis

4.1. Calculation Results ofWeights

The main step of the ACA–PCA–FR model is to obtain the weight of 11 influencing factors through the combination of AHP and PCA. FR is used to obtain the weight of each subclass in each influencing factor. Finally, the susceptibility map is obtained by stacking all factor maps. Based on the historical data and field investigation results of the research area, this paper ranked the 11 influencing factors in order of importance, Rf > H > slope > Pl.Cv > Pr.Cv > Rd > GIE > TWI > TCI > SPI > NDVI. The established judgment matrix is shown in Table 2.
Through calculations, the consistency index CI = 0.1208 and the consistency ratio CR = 0.0795 were obtained. Since CR < 0.10, the consistency of the judgment matrix A was acceptable.
The objective weights can be obtained by AHP. The data of 11 influencing factors of 273 water catchment units in the study area were imported into the SPSS software for principal component analysis, and the weights of 11 influencing factors were calculated by the software, as shown in Table 3. Finally, the combined weight value is obtained by averaging the weights obtained by the two methods. It can be seen from the results that the combined weight is more reasonable than the weight value obtained by AHP and PCA, as shown in Figure 6 and Figure 7.

4.2. Distribution of Debris Flow in Different Classes of Factors

In the process of FR and ACA–PCA–FR modeling, it is necessary to count the ratio of the debris flow grid number in each category and the ratio of each category to the total grid number in each influencing factor. In this paper, each influencing factor is divided into six categories according to the natural discontinuous method, and the frequency ratio of each subcategory of each influencing factor is obtained through the distribution of 43 debris flow trenches investigated in the field, as shown in Figure 8 and Table 4.The FR model is used to make a statistical calculation of the frequency ratio method for 11 evaluation factors by using the probability statistics method.Then, it makes a superposition of the evaluation factors to get the debris flow sensitivity susceptibility map.

4.3. Correlation Analysis

Before the susceptibility assessment of debris flow, many studies carried out independent tests for the evaluation indices selected. However, many susceptibility assessment models are sensitive to the multicollinearity of evaluation indices, such as the LR model. The linear correlation of each evaluation index will increase the prediction error of debris flow susceptibility assessment. Some studies use independent tests to test the mutual independence of each evaluation index [66]. For example, by using the variance inflation factor (VIF) and the conditional independence test [67], evaluation indices with high linear correlation can be proposed. However, compared with these methods, the PCA method can also eliminate the multicollinearity problem among evaluation indices, and has less information loss on the original evaluation indices, that is, the principal components obtained by this method are obtained through all the pre-selected price indices instead of reducing the dimension by directly eliminating some evaluation indicators. Therefore, the PCA method was used to reduce the dimension of the pre-selected evaluation indices, and the re-selected evaluation indices were transformed to make them independent of each other. Then, the sensitivity zoning of debris flow is carried out by using the evaluation index after the dimensionality reduction to eliminate the influence of linear correlation among factors on the prediction results.The correlation matrix of impact factors is shown in Table 5.

5. Discussion

In this paper, the natural break point method was used to classify the susceptibility assessment results of debris flow obtained by the four models into five susceptibility areas, which were very low, low, moderate, high and very high susceptibility areas. Finally, the susceptibility assessment map was obtained, as shown in Figure 9. The percentage of the five susceptibility areas in the total area can be obtained by extraction, as shown in Figure 10. As can be seen from the figure, the susceptibility results of FR and ACA–PCA–FR are very similar. On the whole, the areas of high and very high susceptibility areas in the two methods are large, accounting for about 50% of the total area. The 43 debris flow gullies in the field survey are all located in the highly sensitive area, which is reasonable.
SVM and LR are two machine learning languages widely used in the evaluation of geological hazards. In this study, the results obtained by SVM and LR models are similar. On the whole, the areas of high and very high susceptibility areas in these two methods are significantly smaller than the traditional weight assignment methods, accounting for 34.99% and 29.95% of the total area, respectively. However, the 43 debris flow gullies in the field survey are located in the high susceptibility area, so it is also reasonable.
However, in the process of field investigation, the vegetation in the study area is very dense, and the sources in the gullies are relatively few. There are relatively more sources in the gullies with a large area. Therefore, from the perspective of field investigation, the results obtained by SVM and LR models are more reasonable than the sensitivity results obtained by FR and ACA–PCA–FR models. In order to further verify the accuracy of the four sensitivity models, ROC curves were drawn, and AUC values were calculated for the susceptibility results obtained from the four models.
The ROC curve, also known as the receiver characteristic curve, is widely used in medical laboratories and disease prediction [68]. The individual of each disaster point and its evaluation factor subgroup is equivalent to the subject in the medical laboratory. The occurrence and non-occurrence of debris flow are taken as two categories (positive category and negative category) for analysis [69]. The ROC curve and the area under the curve (AUC) corresponding to each mode were obtained. AUC is the standard to judge the advantages and disadvantages of the four models used. When AUC = 0.5, it means that the model results have no reference value, while when it is less than 0.5, it means that the model does not conform to the real situation. When AUC is greater than 0.5, the closer its value is to 1, the more accurate the model effect is.
The ROC curve and AUC value drawn by the four models are shown in Figure 11. It can be seen from Figure 11 that the prediction accuracy of the SVM model is better than the other three models, with the largest AUC value of 0.889, followed by LR (AUC = 0.842), ACA–PCA–FR (AUC = 0.829), FR (AUC = 0.797). The results show that SVM has low requirements on the amount of data and can still maintain very high accuracy in the case of a small sample data, which is very suitable for the evaluation of geological disasters. Meanwhile, the AUC values of the SVM and LR models are larger than those of LR and ACA–PCA–FR, indicating that the performance of the machine learning algorithm is significantly better than that of the traditional weight method. In addition, the main difference between FR and ACA–PCA–FRis that when determining the weight of each influence factor, the FR model adopts the frequency ratio method, while ACA–PCA–FR determines the weight through the combination of the subjective weight determined by AHP and objective weight determined by PCA. According to the AUC value, it can be concluded that the method of subjective and objective combination weighting is better than the method of objective weighting.
According to the above results, the susceptibility results obtained from the SVM model were taken as the final evaluation results, among which the very low, low, moderate, high and very high susceptibility zones accounted for 22.31%, 25.04%, 17.66%, 18.85% and 16.14%, respectively, of the total area of the study area.
The main advantages of logical regression are simple implementation, wide application, small amount of calculation, fast speed, low storage resources, and the ability to solve multicollinearity problems. The main disadvantage is that it is easy to be under-fitted, and the accuracy is generally not high. It can not handle a large number of multi-class features or variables well. For nonlinear features, it needs to be converted.
The main advantages of the support vector machine are high classification accuracy, accurate classification in small sample size, good generalization ability, and easy to solve nonlinear problems, classification and regression problems of high-dimensional features. The main disadvantage is that there is no universal and effective standard for the selection of kernel functions of support vector machines. At the same time, because there is no reliable sample selection mode, the sample selection method also has a great impact on the evaluation results.
Compared with LR model, SVM model has better accuracy and reliability, and also has better classification ability: SVM tends to divide more watershed units into unstable units, and the over-prediction characteristics of SVM should be fully considered in practical work; however, LR is relatively conservative, and tends to predict more basin units as stable units, and it is easier to divide dangerous areas into stable areas. Therefore, the characteristics of the two models should be fully considered in practical work. Both models have certain limitations, and the single model has defects in varying degrees. In the study area, for the advantages, disadvantages and application conditions of different models, using multiple models to solve problems will also be a problem that needs further study in the future.
Comparing this paper with the research of Ahmad et al. [24], the theme of the two articles is to compare with machine learning and mathematical statistics methods to assess debris flow susceptibility. Among them, there are FR model and LR model in both articles, and seven of the influence factors used are the same. Both articles have concluded that the machine learning model is better than the mathematical statistics model, which also verifies that the conclusion of this paper is indeed correct. The difference is that this paper also makes a deeper comparative study between different machine learning models and different mathematical statistical models, in order to better apply them to the evaluation of geological hazards.

6. Conclusions

Based on field investigations and the 3S technology, this research selects 11 influence factors, such as elevation and slope, and selects the catchment unit that is suitable as the basic mapping unit. By setting different flow thresholds for comparative research, the study area is finally divided into 273 catchment units. Then, four models, R, ACA–PCA–FR, LR and SVM, are established, respectively. In the LR and SVM models, 43 debris flow gullies surveyed in the field are taken as known debris flows. The same number of 43 non-debris flow gullies were obtained based on the Borderline-SMOTE algorithm. These two kinds of data were substituted into the model as test data. Finally, the sensitivity evaluation results of the study area were obtained, and the results were verified by the ROC curve and AUC value, and the following conclusions were drawn:
1.
Among the four models, the SVMmodel has the best performance and the highest prediction accuracy, with AUC = 0.889, followed by LR (AUC = 0.842), ACA–PCA–FR (AUC = 0.829) and FR (AUC = 0.797). The results show that SVM can still maintain very high prediction accuracy in the case of small sample data, learning can be strong and have a fast convergence, and has strong adaptability to high-dimensional samples, which is very suitable for the evaluation and analysis of geological disasters.
2.
Among the four models, the results of the FR and ACA–PCA–FR models are relatively similar. These two methods are traditional weight evaluation methods. According to the field survey results and AUC values, the accuracy of these two methods is relatively low. The results of LR and SVM, as two widely-used machine learning algorithms, are similar, more consistent with the field survey results, and the AUC value is relatively high, so in this study, the machine learning algorithm is more accurate and reasonable.
At the same time, this paper also has some shortcomings. Because it is difficult to obtain historical data and the number of known debris flow gullies is small, this directly affects the reliability of machine learning algorithms. However, there is no doubt that more reasonable and efficient algorithms can be put forward to assess debris flow susceptibility in the future.

Author Contributions

Formal analysis, F.G.; funding acquisition, J.C.; methodology, F.G.; project administration, J.C.; software, Y.L.; supervision, Q.W.; validation, F.G.; visualization, Y.Z.; writing—original draft, F.G.; writing—review and editing, J.C., X.S. and Q.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the key project of the National Natural Science Foundation of China (Grant no. U1702241).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank anonymous reviewers for their comments and suggestions, which helped to improve this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The basic statistics of 43 debris flow gullies have been given in Table A1.
Table A1. The basic statistics of 43 debris flow gullies.
Table A1. The basic statistics of 43 debris flow gullies.
No.H (m)Slope (°)Pl.CvPr.CvTWISPITCIRdGIERfNDVI
1243.6520.130.040.185.0451.76−1.550.420.6685.360.41
2505.3825.920.050.044.63219.43−1.780.370.4289.650.38
31038.4325.570.040.064.67206.57−1.840.430.2175.510.39
4602.1526.640.030.064.70381.37−1.900.460.3990.600.38
5931.0126.220.040.054.71135.54−1.740.430.2177.370.40
6749.6730.090.050.004.56219.89−1.580.420.4780.980.44
71008.4025.660.010.014.84273.88−1.460.430.2177.230.41
8967.5727.500.030.064.77286.41−1.650.410.3180.010.41
9869.6726.780.030.064.72156.73−1.740.450.5280.430.41
10784.3326.030.000.044.4880.38−2.180.450.5380.580.41
11437.4424.66−0.020.094.81119.62−2.060.560.4787.570.46
12269.9428.29−0.010.084.52148.31−2.300.680.4679.670.48
13263.1827.44−0.060.044.54153.34−2.460.430.4579.760.48
14418.2827.990.06−0.014.54309.73−1.560.380.3678.210.44
15738.4626.25−0.020.094.70160.45−2.170.630.1975.400.46
16820.8627.170.040.024.64187.79−1.680.630.1976.350.41
17546.4927.140.000.004.54123.96−1.850.390.3674.290.46
18392.3027.960.010.074.72295.73−1.970.460.3578.190.43
19307.4428.72−0.020.074.49171.67−2.530.680.4679.390.46
20630.3327.450.030.064.70213.02−1.870.500.2588.020.42
21388.8626.870.020.094.68174.03−1.960.430.3983.890.43
22565.6229.210.020.014.64568.68−2.010.550.2570.200.35
23473.6127.670.040.064.74265.07−1.730.560.2585.870.38
24532.1131.210.02−0.014.64196.51−1.520.410.1682.790.42
25389.1223.220.050.064.73163.41−1.840.590.2583.820.42
26432.2923.280.050.054.88249.88−1.750.570.3382.150.40
27598.2126.25−0.040.014.6997.54−1.890.450.2992.200.36
28532.1927.010.040.084.72246.06−1.840.360.3380.850.41
29397.5829.74−0.130.144.47137.77−3.180.450.3691.340.40
30564.0529.76−0.110.054.42108.60−2.720.450.2992.240.34
31411.3525.690.030.074.77320.24−1.930.550.3294.100.38
32525.7325.930.070.084.70140.00−1.750.450.2989.900.41
33410.4315.580.070.045.2181.75−1.080.540.34106.720.41
34404.5813.640.080.125.39138.10−1.260.480.34106.590.41
35461.0825.490.020.114.70123.40−2.140.520.30102.130.40
36508.7626.660.030.064.66261.97−1.910.600.19102.290.39
37453.5523.810.020.044.97223.28−1.410.720.38109.660.37
38655.9828.490.00−0.024.67222.83−1.620.550.2570.360.36
39379.8623.350.060.134.82148.32−1.740.550.34108.490.40
40387.3420.130.060.054.99117.33−1.160.580.38105.980.41
41470.4724.170.010.074.80162.16−1.930.690.36102.310.40
42411.8024.660.020.035.03319.65−1.320.620.3090.490.39
43413.6620.880.010.114.9168.87−1.870.700.3298.690.43

References

  1. Chong, Y.; Chen, G.; Meng, X.; Yang, Y.; Shi, W.; Bian, S.; Zhang, Y.; Yue, D. Quantitative analysis of artificial dam failure effects ondebrisflows-A case study of theZhouqu′8.8′debrisflowin northwestern China. Sci. Total Environ. 2021, 792, 148439. [Google Scholar] [CrossRef] [PubMed]
  2. Zhang, S.; Zhang, L.M.; Chen, H.X. Relationships among three repeated large-scale debris flows at Pubugou Ravine in the Wenchuan earthquake zone. Can. Geotech. J. 2014, 51, 951–965. [Google Scholar] [CrossRef]
  3. Ma, C.; Deng, J.Y.; Wang, R. Analysis of the triggering conditions and erosion of a runoff-triggered debris flow in Miyun County, Beijing, China. Landslides 2018, 15, 2475–2485. [Google Scholar] [CrossRef]
  4. Dormann, C.F.; Elith, J.; Bacher, S.; Buchmann, C.; Carl, G.; Carré, G.; Marquéz, J.R.G.; Gruber, B.; Lafourcade, B.; Leitão, P.J.; et al. Collinearity: A review of methods to deal with it and a simulation study evaluating their performance. Ecography 2013, 36, 27–46. [Google Scholar] [CrossRef]
  5. Robert, H. Soil slumps and debris flows: Prediction and Protection. Bull. Assoc. Eng. Geol. 1981, 18, 17–28. [Google Scholar]
  6. Feizizadeh, B.; Roodposhti, M.S.; Blaschke, T.; Aryal, J. Comparing GIS-based support vector machine kernel functions for landslide susceptibility mapping. Arab. J. Geosci. 2017, 10, 1–13. [Google Scholar] [CrossRef]
  7. Devkota, K.C.; Regmi, A.D.; Pourghasemi, H.R.; Yoshida, K.; Pradhan, B.; Ryu, I.C.; Dhital, M.R.; Althuwaynee, O.F. Landslide susceptibility mapping using certainty factor, index of entropy and logistic regression models in GIS and their comparison at Mugling-Narayanghat road section in Nepal Himalaya. Nat. Hazards 2013, 65, 135–165. [Google Scholar] [CrossRef]
  8. Chen, X.; Chen, H.; You, Y.; Liu, J. Susceptibility assessment of debris flows using the analytic hierarchy process method-a case study in Subao river valley, China. J. Rock Mech. Geotech. Eng. 2015, 7, 404–410. [Google Scholar] [CrossRef]
  9. Xu, W.; Yu, W.; Jing, S.; Zhang, G.; Huang, J. Debris flow susceptibility assessment by GIS and information value model in a large-scale region, Sichuan Province (China). Nat. Hazards 2013, 65, 1379–1392. [Google Scholar] [CrossRef]
  10. Shi, M.Y.; Chen, J.P.; Song, Y.; Zhang, W.; Song, S.Y.; Zhang, X.D. Assessing debris flow susceptibility in Heshigten Banner, Inner Mongolia, China, using principal component analysis and an improved fuzzy C-means algorithm. Bull. Eng. Geol. Environ. 2016, 75, 909–922. [Google Scholar] [CrossRef]
  11. Lee, S.; Sambath, T. Landslide susceptibility mapping in the Damrei Romel area, Cambodia using frequency ratio and logistic regression models. Environ. Geol. 2006, 50, 847–855. [Google Scholar] [CrossRef]
  12. Li, Z.H.; Chen, J.P.; Tan, C.; Zhou, X.; Li, Y.C.; Han, M.X. Debris flow susceptibility assessment based on topo-hydrological factors at different unit scales: A case study of Mentougou district, Beijing. Environ. Earth Sci. 2021, 80, 365. [Google Scholar] [CrossRef]
  13. Li, Y.C.; Chen, J.P.; Tan, C.; Li, Y.; Gu, F.F.; Zhang, Y.W.; Mehmood, Q. Application of the borderline-SMOTE method in susceptibility assessments of debris flows in Pinggu District, Beijing, China. Nat. Hazards 2021, 105, 2499–2522. [Google Scholar] [CrossRef]
  14. Kang, S.H.; Lee, S.R. Debris flow susceptibility assessment based on an empirical approach in the central region of South Korea. Geomorphology 2018, 308, 1–12. [Google Scholar] [CrossRef]
  15. Liang, Z.; Wang, C.M.; Zhang, Z.M. A comparison of statistical and machine learning methods for debris flow susceptibility mapping. Stoch. Environ. Res. Risk A 2020, 34, 1887–1907. [Google Scholar] [CrossRef]
  16. Lin, G.F.; Chang, M.J.; Huang, Y.C.; Ho, J.Y. Assessment of susceptibility to rainfall induced landslides using improved self-organizing linear output map, support vector machine, and logistic regression. Eng. Geol. 2017, 224, 62–74. [Google Scholar] [CrossRef]
  17. Cao, C.; Zhang, W.; Chen, J.P.; Shan, B.; Song, S.Y.; Zhan, J.W. Quantitative estimation of debris flow source materials by integrating multi-source data: A case study. Eng. Geol. 2021, 291, 106222. [Google Scholar] [CrossRef]
  18. Sun, X.H.; Chen, J.P.; Bao, Y.D.; Han, X.D.; Zhan, J.W.; 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]
  19. McSherry, D. Strategic induction of decision trees. Knowl. Based Syst. 1988, 12, 269–275. [Google Scholar] [CrossRef]
  20. Sun, X.H.; Chen, J.P.; Li, Y.R.; Rene, N.N. Landslide Susceptibility Mapping along a Rapidly Uplifting River Valley of the Upper Jinsha River, Southeastern Tibetan Plateau, China. Remote Sens. 2022, 14, 1730. [Google Scholar] [CrossRef]
  21. Qiu, C.C.; Su, L.J.; Zou, Q.; Geng, X.Y. A hybrid machine-learning model to map glacier-related debris flow susceptibility along Gyirong Zangbo watershed under the changing climate. Sci. Total Environ. 2022, 818, 151752. [Google Scholar] [CrossRef]
  22. Ke, X.; Basanta, R. Comparison of Different Machine Learning Methods for Debris Flow Susceptibility Mapping: A Case Study in the Sichuan Province, China. Remote Sens. 2020, 12, 295. [Google Scholar]
  23. Elkadiri, R.; Sultan, M.; Youssef, A.M.; Elbayoumi, T.; Chase, R.; Bulkhi, A.B.; Al-Katheeri, M.M. A remote sensing-based approach for debris-flow susceptibility assessment using artificial neural networks and logistic regression modeling. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 4818–4835. [Google Scholar] [CrossRef]
  24. Ahmad, H.; Chen, N.S.; Rahman, M.; Islam, M.M.; Pourghasemi, H.R.; Hussain, S.F.; Habumugisha, J.M.; Liu, E.L.; Zheng, H.; Ni, H.Y.; et al. Geohazards Susceptibility Assessment along the Upper Indus Basin Using Four Machine Learning and Statistical Models. ISPRS Int.Geo-Inf. 2021, 10, 315. [Google Scholar] [CrossRef]
  25. Pal, S.C.; Chakrabortty, R.; Saha, A.; Bozchaloei, S.K.; Pham, Q.B.; Linh, N.T.T.; Anh, D.T.; Janizadeh, S.; Ahmadi, K. Evaluation of debris flow and landslide hazards using ensemble framework of Bayesian- and tree-based models. Bull. Eng. Geol. Environ. 2022, 81, 5. [Google Scholar] [CrossRef]
  26. Ciccarese, G.; Mulas, M.; Corsini, A. Combining spatial modelling andregionalization ofrainfall thresholds fordebris flows hazard map- ping in the Emilia-Romagna Apennines (Italy). Landslides 2021, 18, 3513–3529. [Google Scholar]
  27. Vianello, D.; Vagnon, F.; Bonetto, S.; Mosca, P. Debris flowsusceptibilitymapping using the Rock Engineering System (RES) method: A case study. Landslides 2022, 1–22. [Google Scholar] [CrossRef]
  28. Cao, J.; Zhang, Z.; Du, J.; Zhang, L.; Song, Y.; Sun, G. Multi-geohazards susceptibility mapping based on machine learning—A case study in Jiuzhaigou, China. Nat. Hazards 2020, 102, 851–871. [Google Scholar] [CrossRef]
  29. Ma, C.; Wang, Y.J.; Du, C.; Wang, Y.Q.; Li, Y.P. Variation in initiation condition of debris flow in the mountain regions surrounding Beijing. Geomorphology 2016, 273, 323–334. [Google Scholar] [CrossRef]
  30. Li, Y.; Chen, J.; Zhang, Y.; Song, S.; Han, X.; Ammar, M. Debris flow susceptibility assessment and runout prediction: A case study in Shiyang Gully, Beijing, China. Int J. Environ. Res. 2020, 14, 365–383. [Google Scholar] [CrossRef]
  31. Cheng, W.M.; Wang, N.; Zhao, M.; Zhao, S.M. Relative tectonics and debris flow hazards in the Beijing mountain area from DEM-derived geomorphic indices and drainage analysis. Geomorphology 2016, 257, 134–142. [Google Scholar] [CrossRef]
  32. Han, H.; Wang, W.Y.; Mao, B.H. Borderline-SMOTE: A new over-sampling method in imbalanced data sets learning. Lect. Notes Comput. Sci. 2005, 3644, 878–887. [Google Scholar]
  33. Qin, S.; Lv, J.; Cao, C.; Ma, Z.; Hu, X.; Liu, F.; Qiao, S.; Dou, Q. Mapping debris flow susceptibility based on watershed unit and grid cell unit: A comparison study. Geomat. Nat. Hazards Risk 2019, 10, 1648–1666. [Google Scholar] [CrossRef]
  34. 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]
  35. Tiranti, D.; Cavalli, M.; Crema, S.; Zerbato, M.; Graziadei, M.; Barbero, S.; Cremonini, R.; Silvestro, C.; Bodrato, G.; Tresso, F. Semi-quantitative method for the assessment of debris supply from slopes to river in ungauged catchments. Sci. Total Environ. 2016, 554, 337–348. [Google Scholar] [CrossRef]
  36. Zhang, Y.; Ge, T.; Tian, W.; Liou, Y.-A. Debris flow susceptibility mapping using machine-learning techniques in Shigatse Area, China. Remote Sens. 2019, 11, 2801. [Google Scholar] [CrossRef]
  37. Yilmaz, I. Comparison of landslide susceptibility mapping methodologies for Koyulhisar, Turkey: Conditional probability, logistic regression, artificial neural networks, and support vector machine. Environ. Earth Sci. 2010, 61, 821–836. [Google Scholar] [CrossRef]
  38. Suthaharan, S. Support Vector Machine. In Machine Learning Models and Algorithms for Big Data Classification; Springer: Boston, MA, USA, 2016; pp. 207–235. [Google Scholar]
  39. Xu, C.; Dai, F.; Xu, X.; Lee, Y.H. GIS-based support vector machine modeling of earthquake-triggered landslide susceptibility in the Jianjiang River watershed, China. Geomorphology 2012, 145–146, 70–80. [Google Scholar] [CrossRef]
  40. Li, T.; Qiu, S.; Mao, S.X.; Bao, R.; Deng, H.B. Evaluating water resource accessibility in Southwest China. Water 2019, 11, 1708. [Google Scholar] [CrossRef]
  41. Wang, Q.; Guo, Y.; Li, W.; He, J.; Wu, Z. Predictive modeling of landslide hazards in Wen County, northwestern China based on information value, weights-of-evidence, and certainty factor. Geomat. Nat. Hazards Risk 2019, 10, 820–835. [Google Scholar] [CrossRef]
  42. 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]
  43. Moore, I.D.; Grayson, R.; Ladson, A. Digital terrain modeling: A review of hydrological, geomorphological, and biological applications. Hydrol Process 1991, 5, 3–30. [Google Scholar] [CrossRef]
  44. Beven, K.J.; Kirkby, M.J. A physically based variable contributing area model of basin hydrology. Hydrol. Sci. Bull. 1979, 24, 43–69. [Google Scholar] [CrossRef]
  45. Wilson, J.P.; Gallant, J.C. Digital Terrain Analysis. In Terrain Analysis: Principles and Applications; Wilson, J.P., Gallant, J.C., Eds.; Wiley: New York, NY, USA, 2000; pp. 1–27. [Google Scholar]
  46. Park, S.J.; McSweeney, K.; Lowery, B. Identifcation of the spatial distribution of soils using a process-based terrain characterization. Geoderma 2001, 103, 249–272. [Google Scholar] [CrossRef]
  47. Tang, C.; Zhu, J.; Li, W.L. Rainfall-triggered debris flows following the Wenchuan earthquake. Bull. Eng. Geol. Environ. 2009, 68, 187–194. [Google Scholar] [CrossRef]
  48. Tehrany, M.S.; Pradhan, B.; Mansor, S.; Ahmad, N. Flood susceptibility assessment using GIS based support vector machine model with different kernel types. Catena 2015, 125, 91–101. [Google Scholar] [CrossRef]
  49. Saaty, T.L. Modeling unstructured decision problems—The theory of analytical hierarchies. Math. Comput. Simul. 1978, 20, 147–158. [Google Scholar] [CrossRef]
  50. Chen, W.; Wang, J.; Xie, X.; Hong, H.; Van Trung, N.; Bui, D.T.; Wang, G.; Li, X. Spatial prediction of landslide susceptibility using integrated frequency ratio with entropy and support vector machines by different kernel functions. Environ. Earth Sci. 2016, 75, 1344. [Google Scholar] [CrossRef]
  51. Sabokbar, H.F.; Roodposhti, M.S.; Tazik, E. Landslide susceptibility mapping using geographically-weighted principal component analysis. Geomorphology 2014, 226, 15–24. [Google Scholar] [CrossRef]
  52. Zhou, C.; Yin, K.; Cao, Y.; Ahmed, B. Application of time series analysis and PSO–SVM model in predicting the Bazimen landslide in the Three Gorges Reservoir, China. Eng. Geol. 2016, 204, 108–120. [Google Scholar] [CrossRef]
  53. Yilmaz, I. Landslide susceptibility mapping using frequency ratio, logistic regression, artificial neural networks and their comparison: A case study from Kat landslides (Tokat-Turkey). Comput. Geosci. 2009, 35, 1125–1138. [Google Scholar] [CrossRef]
  54. Zhang, X.; Wu, Y.; Zhai, E.; Ye, P. Coupling analysis of the heat-water dynamics and frozen depth in a seasonally frozen zone. J. Hydrol. 2021, 593, 125603. [Google Scholar] [CrossRef]
  55. Wu, X.; Ren, F.; Niu, R. Landslide susceptibility assessment using object mapping units, decision tree, and support vector machine models in the Three Gorges of China. Environ. Earth Sci. 2013, 71, 4725–4738. [Google Scholar] [CrossRef]
  56. Bui, D.T.; Tuan, T.A.; Hoang, N.D.; Thanh, N.Q.; Nguyen, D.B.; Van Liem, N.; Pradhan, B. Spatial prediction of rainfall-induced landslides for the Lao Cai area (Vietnam) using a hybrid intelligent approach of least squares support vector machines inference model and artificial bee colony optimization. Landslides 2016, 14, 447–458. [Google Scholar]
  57. Bui, D.T.; Pradhan, B.; Lofman, O.; Revhaug, I. Landslide susceptibility assessment in Vietnam using support vector machines, decision tree, and Naïve Bayes models. Math. Probl. Eng. 2012, 2012, 1–26. [Google Scholar]
  58. Zhang, X.; Zhai, E.; Wu, Y.; Sun, D.; Lu, Y. Theoretical and Numerical Analyses on Hydro–Thermal–Salt–Mechanical Interaction of Unsaturated Salinized Soil Subjected to Typical Unidirectional Freezing Process. Int. J. Geomech. 2021, 21, 04021104. [Google Scholar] [CrossRef]
  59. Pradhan, B. A comparative study on the predictive ability of the decision tree, support vector machine and neuro-fuzzy models in landslide susceptibility mapping using GIS. Comput. Geosci. 2013, 51, 350–365. [Google Scholar] [CrossRef]
  60. Umar, Z.; Pradhan, B.; Ahmad, A.; Jebur, M.N.; Tehrany, M.S. Earthquake induced landslide susceptibility mapping using an integrated ensemble frequency ratio and logistic regression models in West Sumatera Province, Indonesia. Catena 2014, 118, 124–135. [Google Scholar] [CrossRef]
  61. Bui, D.T.; Tuan, T.A.; 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 2015, 13, 361–378. [Google Scholar]
  62. Pradhan, B.; Lee, S. Landslide susceptibility assessment and factor effect analysis: Backpropagation artificial neural networks and their comparison with frequency ratio and bivariate logistic regression modelling. Environ. Model. Softw. 2010, 25, 747–759. [Google Scholar] [CrossRef]
  63. Pham, B.T.; Bui, D.T.; Prakash, I.; Nguyen, L.H.; Dholakia, M.B. A comparative study of sequential minimal optimization-based support vector machines, vote feature intervals, and logistic regression in landslide susceptibility assessment using GIS. Environ. Earth Sci. 2017, 76, 371. [Google Scholar] [CrossRef]
  64. Zhang, X.; Ye, P.; Wu, Y. Enhanced technology for sewage sludge advanced dewatering from an engineering practice perspective: A review. J. Environ. Manag. 2022, 321, 115938. [Google Scholar] [CrossRef] [PubMed]
  65. Meng, Q.; Miao, F.; Zhen, J.; Wang, X.; Wang, A.; Peng, Y.; Fan, Q. GIS-based landslide susceptibility mapping with logistic regression, analytical hierarchy process, and combined fuzzy and support vector machine methods: A case study from Wolong Giant Panda Natural Reserve, China. Bull. Eng. Geol. Environ. 2015, 75, 923–944. [Google Scholar] [CrossRef]
  66. Marjanović, M.; Kovačević, M.; Bajat, B.; Voženílek, V. Landslide susceptibility assessment using SVM machine learning algorithm. Eng. Geol. 2011, 123, 225–234. [Google Scholar] [CrossRef]
  67. Mather, P.M. The use of backpropagating artificial neural networks in land cover classification. Int. J. Remote Sens. 2003, 24, 4907–4938. [Google Scholar]
  68. Li, L.; Lan, H.; Guo, C.; Zhang, Y.; Li, Q.; Wu, Y. A modified frequency ratio method for landslide susceptibility assessment. Landslides 2017, 14, 727–741. [Google Scholar] [CrossRef]
  69. Kavzoglu, T.; Sahin, E.K.; Colkesen, I. Landslide susceptibility mapping using GIS-based multi-criteria decision analysis, support vector machines, and logistic regression. Landslides 2013, 11, 425–439. [Google Scholar] [CrossRef]
Figure 1. Debris flow distribution in Beijing and geographical location of study area. (a) The location of the research area in the whole of Beijing; (b) The location of the research area in the whole of China.
Figure 1. Debris flow distribution in Beijing and geographical location of study area. (a) The location of the research area in the whole of Beijing; (b) The location of the research area in the whole of China.
Water 15 00705 g001
Figure 2. The geological and tectonic map of the study area.
Figure 2. The geological and tectonic map of the study area.
Water 15 00705 g002
Figure 3. Field investigation in the study area. (a) Waste broken slag; (b) loose material source in the channel; (c) water erosion in the channel; (d) small collapse on both sides of the channel.
Figure 3. Field investigation in the study area. (a) Waste broken slag; (b) loose material source in the channel; (c) water erosion in the channel; (d) small collapse on both sides of the channel.
Water 15 00705 g003
Figure 4. The modeling flow chart of this research.
Figure 4. The modeling flow chart of this research.
Water 15 00705 g004
Figure 5. Maps showing causative factors in this study area with GIS.
Figure 5. Maps showing causative factors in this study area with GIS.
Water 15 00705 g005
Figure 6. Results of combined weights of each impact factor.
Figure 6. Results of combined weights of each impact factor.
Water 15 00705 g006
Figure 7. Weights of each impact factor.
Figure 7. Weights of each impact factor.
Water 15 00705 g007
Figure 8. Distribution of debris flow in different classes of factors.
Figure 8. Distribution of debris flow in different classes of factors.
Water 15 00705 g008aWater 15 00705 g008b
Figure 9. Debris flow susceptibility maps produced by four models.
Figure 9. Debris flow susceptibility maps produced by four models.
Water 15 00705 g009
Figure 10. Distribution of the different debris flow susceptibility classes from the four models.
Figure 10. Distribution of the different debris flow susceptibility classes from the four models.
Water 15 00705 g010
Figure 11. The ROC curves of four debris flow susceptibility models.
Figure 11. The ROC curves of four debris flow susceptibility models.
Water 15 00705 g011
Table 1. Average random consistency indicator value table.
Table 1. Average random consistency indicator value table.
n1234567891011
RI000.580.901.121.241.321.411.451.491.51
Table 2. Judgment matrix of impact factors.
Table 2. Judgment matrix of impact factors.
FactorsX1X2X3X4X5X6X7X8X9X10X11Weight
X1123445656730.260
X21/212334223330.154
X31/31/21223233420.118
X41/41/31/2122323230.095
X51/41/31/21/212332320.084
X61/51/41/31/21/21223230.066
X71/61/21/21/31/31/2123320.060
X81/51/21/31/21/31/21/212340.055
X91/61/31/31/31/21/31/31/21430.044
X101/71/31/41/21/31/21/31/31/4120.031
X111/31/31/21/31/21/31/21/41/31/210.033
X1 = Rf, X2 = H, X3 = slope, X4 = Pl.Cv, X5 = Pr.Cv, X6 = Rd, X7 = GIE, X8 = TWI, X9 = TCI, X10 = SPI, X11 = NDVI.
Table 3. Results of combined weights of each impact factor.
Table 3. Results of combined weights of each impact factor.
FactorsWeight
AHPPCAAverage
Rf0.2600.0830.172
H0.1540.1010.127
Slope0.1180.0910.104
Pl.Cv0.0950.0820.088
Pr.Cv0.0840.0790.082
Rd0.0660.0920.079
GIE0.0600.0860.073
TWI0.0550.1050.080
TCI0.0440.0990.072
SPI0.0310.0880.059
NDVI0.0330.0930.063
Table 4. The frequency ratio of different classes of influence factors.
Table 4. The frequency ratio of different classes of influence factors.
FactorClassStudy AreaDebris Flows AreaFR
CountRatio (%)CountRatio (%)
Elevation
(m)
<2461,763,58212.7732,3733.540.278
246–3522,943,22121.3032,8573.600.169
352–4472,749,38319.90232,60625.471.280
447–5712,615,98718.94314,92734.481.821
571–7732,147,02215.54164,59718.021.160
>10661,596,18011.55136,00714.891.289
Slope (°)<10355,1832.573280.040.014
10–181,706,12412.3555030.600.049
18–222,289,90516.58113,06912.380.747
22–253,142,47722.75232,49425.451.119
25–283,564,60425.80441,17648.301.872
>282,757,08219.96120,79713.230.663
Pl.Cv<0.020989,0537.1634,7873.810.532
0.020–0.0342,834,11020.51278,74430.521.488
0.034–0.0463,350,59224.25357,43739.131.614
0.046–0.0574,254,67030.80204,98922.440.729
0.057–0.0782,036,61014.7435,3193.870.262
>0.078350,3402.5420910.230.090
Pr.Cv<0.0018776,0905.6234,5463.780.673
0.0018–0.0332,850,84920.64135,74014.860.720
0.033–0.0464,174,34730.22454,51049.761.647
0.046–0.0573,751,04727.15235,96225.830.951
0.057–0.0781,892,81213.7050,1865.490.401
>0.078370,2302.6824230.270.099
TWI<4.73,530,25825.55334,72136.651.434
4.7–4.93,486,82025.24381,47941.771.655
4.9–5.13,477,97125.17161,00517.630.700
5.1–5.62,304,52216.6835,3933.880.232
5.6–6.5748,5185.424410.050.009
>6.5267,2861.933260.040.018
SPI<3217,044,59450.99681,34174.61.463
321–6523,294,69623.85197,75021.650.908
652–11871,661,89412.0323,4402.570.213
1187–2079626,9984.5483920.920.202
2079–32,2591,063,2587.703530.040.005
>32,259123,9350.9020910.230.255
TCI<−1.761,830,82613.25201,63422.081.666
−1.76–1.62,956,30921.40291,42131.911.491
−1.6–1.443,267,55223.65228,03524.971.056
−1.44–1.242,652,37919.20170,63418.680.973
−1.24–0.872,401,89617.3921,3152.330.134
>−0.87706,4135.113280.040.007
Roundness<0.231,149,6358.325760.060.008
0.23–0.322,869,93520.77249,89227.361.317
0.32–0.422,807,72720.32174,58319.110.941
0.42–0.552,830,49920.49110,03712.050.588
0.55–0.812,530,30018.32188,65020.651.128
>0.811,627,27911.78189,62920.761.763
GIE<0.331,518,87110.99151,47316.581.508
0.33–0.434,122,93429.84328,98136.021.207
0.43–0.53,229,95223.38354,52538.821.660
0.5–0.562,467,23717.8674,4208.150.456
0.56–0.641,953,78114.1431990.350.025
>0.64522,6003.787690.080.022
Rainfall
(mm)
<74.52,053,68514.8712,6871.390.093
74.5–80.42,559,07518.52293,25932.111.733
80.4–86.12,344,33516.97126,49313.850.816
86.1–94.22,629,60519.03255,50127.971.470
94.2–103.41,974,84314.29104,32311.420.799
>103.42,253,82216.31121,10413.260.813
NDVI<0.35454,5003.29000.000
0.35–0.382,669,16219.32122,15613.370.692
0.38–0.393,380,91024.47315,81934.581.413
0.39–0.413,899,55528.23284,58731.161.104
0.41–0.432,474,05817.91150,11616.440.918
>0.43937,1906.7840,6894.450.657
Table 5. Parameter correlation matrix.
Table 5. Parameter correlation matrix.
FactorsGIEHNDVIPl.CvPr.CvSlopeRfSPITCITWIRd
GIE1.000−0.550−0.3170.1990.188−0.7270.1910.1660.6580.696−0.257
H−0.5501.0000.155−0.170−0.1390.603−0.535−0.115−0.488−0.5240.210
NDVI−0.3170.1551.000−0.029−0.0410.378−0.277−0.116−0.383−0.4320.116
Pl.Cv0.199−0.170−0.0291.0000.980−0.1120.1330.277−0.019−0.040−0.032
Pr.Cv0.188−0.139−0.0410.9801.000−0.0890.1050.311−0.046−0.056−0.024
Slope−0.7270.6030.378−0.112−0.0891.000−0.299−0.005−0.821−0.9040.388
Rf0.191−0.535−0.2770.1330.105−0.2991.000−0.1430.1830.130−0.041
SPI0.166−0.115−0.1160.2770.311−0.005−0.1431.0000.0780.132−0.021
TCI0.658−0.488−0.383−0.019−0.046−0.8210.1830.0781.0000.923−0.300
TWI0.696−0.524−0.432−0.040−0.056−0.9040.1300.1320.9231.000−0.362
Rd−0.2570.2100.116−0.032−0.0240.388−0.041−0.021−0.300−0.3621.000
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gu, F.; Chen, J.; Sun, X.; Li, Y.; Zhang, Y.; Wang, Q. Comparison of Machine Learning and Traditional Statistical Methods in Debris Flow Susceptibility Assessment: A Case Study of Changping District, Beijing. Water 2023, 15, 705. https://doi.org/10.3390/w15040705

AMA Style

Gu F, Chen J, Sun X, Li Y, Zhang Y, Wang Q. Comparison of Machine Learning and Traditional Statistical Methods in Debris Flow Susceptibility Assessment: A Case Study of Changping District, Beijing. Water. 2023; 15(4):705. https://doi.org/10.3390/w15040705

Chicago/Turabian Style

Gu, Feifan, Jianping Chen, Xiaohui Sun, Yongchao Li, Yiwei Zhang, and Qing Wang. 2023. "Comparison of Machine Learning and Traditional Statistical Methods in Debris Flow Susceptibility Assessment: A Case Study of Changping District, Beijing" Water 15, no. 4: 705. https://doi.org/10.3390/w15040705

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