^{1}

^{2}

This article takes Pingshan County of Yibin City as an example. After collecting various basic data and reviewing the literature, the slope, aspect, Normalized Difference Vegetation Index (NDVI), distance from water system, distance from road, stratum, distance from fault, rainfall 8 kinds of influencing factors are used as evaluation factors. The index of entropy model (index of entropy, IOE) and information values (IV) were used to evaluate the susceptibility of landslide disasters in the study area. On this basis, a weighted information values model (WIV) was proposed. Using the AUC area under the ROC curve to test and compare the three models, the results show that the success rates of the three models of IOE, IV and WIV are 76.2%, 75.6%, and 79.1%, respectively. Finally, according to the natural discontinuity method, the landslide susceptibility index calculated by ArcGIS software is divided into extremely high, high, medium, and low susceptible areas, and the landslide susceptibility zoning map of Pingshan County is obtained, which can provide reference for relevant departments for disaster prevention and mitigation.

Landslide disasters are very serious all over the world, and our country is one of the countries with the most landslide disasters. Every year, the casualties and economic losses caused by landslide disasters are very large. In recent years, with the acceleration of my country’s economic development, the development of mountainous areas has also increased, and the number and scale of landslides have also increased sharply. Pingshan County, Yibin City is a mountainous area. The increasing number of landslide disasters has seriously affected social stability and economic development, and threatened people’s lives. Therefore, the evaluation of landslide susceptibility in this area can provide a scientific basis for local relevant departments to carry out disaster prevention and mitigation and carry out various engineering activities.

Landslide susceptibility evaluation research lays the foundation for subsequent prediction, prevention and management of landslide disasters, and it is also an important scientific basis for the orderly progress of various engineering activities. At present, Chinese and foreign scholars have done a lot of research on the evaluation of landslide susceptibility, and put forward many different evaluation methods. These methods can be roughly distinguished based on qualitative and quantitative, and quantitatively divided into semi-quantitative evaluation and quantitative evaluation [

Pingshan County is located in the west of Yibin City, Yibin County to the east, Suijiang County of Zhaotong City, Yunnan Province to the south across the Jinsha River, Leibo County and Mabian Yi Autonomous Prefecture to the west, and Muchuan County to the north. Pingshan County has jurisdiction over 8 towns and 3 townships (including 2 Yi townships). The geographic location is the southern margin of the Sichuan Basin, on the north bank of the lower reaches of the Jinsha River, at 103˚36'31'' - 104˚23'36'' east longitude and 28˚28'26'' - 28˚53'51'' north latitude [^{2}. The county’s annual average temperature is 14.9˚C, and the annual average rainfall is 1066 mm. It belongs to the mid-subtropical humid monsoon climate zone, with large vertical climate differences. There are many rivers in the county area, mainly including the Jinsha River and the Minjiang River, as well as larger rivers such as Xining River, Zhongdu River and Longxi River.

In terms of geomorphology, the eastern part of Pingshan County is the hills of southern Sichuan, the west and north are connected with the mountains of southwest Sichuan, and the south is connected with the eastern Yunnan plateau [

The information model is based on the actual occurrence of historical landslide disasters in the study area and the information provided by various geological surveying and mapping data, and transforms the factors that have a greater impact on the stability of the landslide into the information value reflecting the stability of the landslide in the study area. The greater the information value, the greater the degree of correlation between the influencing factor and the occurrence of landslide disasters [

I ( X i , Y ) = ln N i / N S i / S (2)

In formula (2): N_{i}/N is the factor corresponding to the proportion of graded landslides, and S_{i}/S is the factor corresponding to the proportion of graded area.

Generally speaking, landslide disasters occur under the combined action of multiple influencing factors. Therefore, the sum of the information values of various influencing factors in the same study area is the total information value of the study area.

I i = ∑ i = 1 p ln ( N i / N S i / S ) (3)

The entropy index model is a classification prediction model based on probability theory and mathematical statistical analysis. It can also be used as the input data of a mixed model. The model can calculate the weight of the variable in the event through the attribute value of each input variable, and the weight of each variable is represented by the entropy index. The value of the entropy index is positively correlated with the effect of the influence factor on the landslide. The weight value is calculated by the following formula:

P i j = b a (4)

( P i j ) = P i j ∑ j = 1 s P i j (5)

H i = − ∑ j = 1 s ( P i j ) log 2 ( P i j ) , i = 1 , 2 , ⋯ , n (6)

H i max = log 2 s (7)

I i = H i max − H i H i max , i = 1 , 2 , ⋯ , n (8)

P i = 1 S j ∑ j = 1 s P i j (9)

W i = I i P i (10)

In the above formula: i is the evaluation factor; j is the grade of each factor; a and b are the area and the percentage of the number of landslides; P_{ij} is the frequency ratio of the factor grade; (P_{ij}) is the probability density; H_{i} is the entropy value; H_{i}_{max} is the maximum of the corresponding factor Entropy; s is the number of factor levels; I_{i} is the information rate; W_{i} is the entropy index value.

The entropy index model only considers the first-level index weight of the impact factor, and the information model only considers the relative weight of the second-level index of the impact factor, which cannot reflect the true weight of each factor. Based on this, the model multiplies the weight calculated by the entropy index model with the weight calculated by the information volume model to obtain the comprehensive weight of the factor, which is closer to its true weight. Then the weight values of each factor are superimposed and summed to obtain the comprehensive weight value of the entire study area, that is, the landslide hazard susceptibility index LSI. The specific calculation formula is as shown in formula (11):

LSI WIV = ∑ i = 1 F W i I V i j (11)

F is the number of impact factors.

All data mainly include the Pingshan County landslide catalog database, 30 m × 30 m Digital Elevation Model (DEM), road network map, water system map, geological map with a scale of 1:50,000, remote sensing image data of Sentinel No. 2 with a resolution of 10 m, and rainfall. Among them, DEM is downloaded from the cloud of geospatial data, geological maps are downloaded from 91 Satellite Map Assistant, remote sensing images are downloaded from the ESA website, and rainfall maps are obtained from rainfall site data downloaded by China Meteorological Data Network through Kriging interpolation. And randomly select non-landslide points equal to the number of landslide points.

The selection of evaluation factors is a very critical step in the evaluation of landslide hazard susceptibility. The incorrect selection of evaluation factors will affect the accuracy of the evaluation results. Too few evaluation factors will cause the evaluation model to be under-fitted and affect the evaluation accuracy. Too many numbers will cause over-fitting of the model, reduce the calculation speed of the entire model, and also affect the accuracy of the model.

Based on existing data and field surveys, combined with predecessors’ experience in selecting factors affecting landslide hazards [

With reference to previous experience [

There are many factors that affect the occurrence of landslides. Some of the influencing factors may have a high correlation between the two, and the relatively

high correlation factors may reduce the operating speed of the model and make the model more complicated. Reduce the accuracy of the entire model. Therefore, when evaluating the susceptibility of landslide hazards, the selected evaluation factors must have low correlation or irrelevance. This paper uses Kendall’s correlation evaluation method, and the correlation coefficient is represented by R. When the absolute value of R is less than |R| < 0.3, it indicates that the correlation between factors is weak, and the factor can be used for landslide susceptibility evaluation. In addition, the correlation between factors is required to have significant reference value, that is, the significance level of correlation between factors sig < 0.05. The data of each factor obtained after ArcGIS software processing is imported into SPSS software for correlation analysis to obtain the correlation coefficient between factors (

DEM | slope | aspect | NDVI | road | stratum | fault | water system | rainfall | |
---|---|---|---|---|---|---|---|---|---|

DEM | 1 | ||||||||

slope | 0.161 | 1 | |||||||

aspect | 0.064 | 0.070 | 1 | ||||||

NDVI | 0.364 | 0.191 | 0.049 | 1 | |||||

road | 0.368 | 0.102 | 0.046 | 0.261 | 1 | ||||

stratum | −0.037 | −0.051 | 0.103 | −0.050 | 0.094 | 1 | |||

fault | 0.110 | 0.047 | −0.011 | 0.059 | −0.061 | −0.103 | 1 | ||

water system | 0.282 | −0.017 | 0.091 | 0.199 | 0.273 | 0.057 | 0.121 | 1 | |

rainfall | −0.058 | −0.118 | 0.062 | 0.260 | 0.056 | 0.040 | −0.100 | 0.039 | 1 |

Before evaluating the landslide susceptibility of the study area, in order to make subsequent calculations simple and convenient and ensure the scientificity of the results, it is necessary to re-divide the entire study area and divide all evaluation factors into a unified measurement. Because the grid unit has the advantages of being simple and easy to superimpose, it is widely used when dividing the research area. Use the empirical formula obtained from the DEM test to calculate the appropriate grid size [

G s = 7.49 + 0.0006 S − 2.0 × 10 − 9 S 2 + 2.9 × 10 − 15 C 3 (12)

In formula (12): C is the scale of the topographic map.

Through calculation, a grid unit with a size of 30 m × 30 m is selected as the grid in this study area, with a total of 2,166,304 grid units.

The information value IV_{ij} of each factor calculated by the above formula is shown in

In ArcGIS software, the information value of different evaluation factors in the study area is superimposed to obtain the total information value of the entire study area, that is, the landslide hazard susceptibility index LSI, which ranges from −4.85 to 4.22. Use the natural discontinuity method to classify the LSI, which are extremely high, high, medium, and low prone areas, and draw a zoning map of the proneness level, as shown in

The relative weight (W_{i})of the first-level index of each evaluation factor obtained by using the above-mentioned entropy index model calculation formula is shown in

L S I I O E = ( N D V I × 0.11 ) + ( s l o p e × 0.12 ) + ( r a i n f a l l × 0.04 ) + ( a s p e c t × 0.06 ) + ( w a t e r s y s t e m × 0.07 ) + ( f a u l t × 0.01 ) + ( s t r a t u m × 0.08 ) + ( r o a d × 0.06 ) (13)

The value range of Pingshan County LSI calculated by the entropy index model is 0.72 - 3.22. The LSI is classified using the natural discontinuity method, which are extremely high, high, medium, and low prone areas, and the landslide proneness level is drawn. Figure, as shown in

The weighted information value under the classification state of each evaluation factor calculated by formula (11) is shown in

This article uses the ROC curve (success rate curve) method to evaluate the model [

Factor | slope | NDVI | rainfall | aspect | Water system | fault | stratum | road |
---|---|---|---|---|---|---|---|---|

IOE W_{i} | 0.12 | 0.11 | 0.04 | 0.06 | 0.07 | 0.01 | 0.08 | 0.06 |

Factor | Grade | IV_{ij} | W_{i}IV_{ij} | Factor | Grade | IV_{ij} | W_{i}IV_{ij} |
---|---|---|---|---|---|---|---|

Slope | 0 - 11.92 | 0.15 | 0.02 | Water system | 200 | 0.47 | 0.03 |

11.92 - 21.00 | 0.35 | 0.04 | 400 | 0.41 | 0.03 | ||

21.00 - 30.40 | −0.21 | −0.03 | 600 | −0.57 | −0.04 | ||

30.40 - 41.69 | −0.47 | −0.06 | 800 | 0.69 | 0.05 | ||

41.69 - 80.25 | −1.02 | −0.12 | >800 | −0.21 | −0.01 | ||

NDVI | <−0.10 | 0.82 | 0.09 | Fault | 2000 | 0.02 | 0.00023 |

−0.10 - 0.41 | 0.02 | 0.002 | 4000 | 0.03 | 0.00034 | ||

0.41 - 0.65 | 1.00 | 0.11 | 6000 | 0.25 | 0.0026 | ||

0.65 - 0.79 | 0.70 | 0.08 | 8000 | −0.33 | −0.003 | ||

0.79 - 0.98 | −0.59 | −0.07 | >8000 | −0.01 | 0.0001 | ||

Rainfall | <935.02 | −0.72 | −0.03 | Stratum | T_{3} | 0.32 | 0.03 |

935.02 - 957.80 | −0.30 | −0.01 | T_{1} | −1.4 | −0.12 | ||

957.80 - 978.30 | 0.12 | 0.004 | K_{2} | −0.56 | −0.05 | ||

978.30 - 993.33 | 0.26 | 0.009 | T_{2} | −1.42 | −0.12 | ||

>993.33 | −0.57 | −0.02 | K_{1} | −0.73 | −0.06 | ||

Aspect | North | 0.02 | 0.001 | J_{2} | 0.03 | 0.003 | |

North east | 0.34 | 0.02 | J_{3} | 0.24 | 0.02 | ||

East | −0.17 | −0.01 | J_{1} | −0.15 | −0.01 | ||

Southeast | 0.18 | 0.01 | P | 0 | 0 | ||

South | 0.02 | 0.001 | Road | 200 | 0.87 | 0.05 | |

South West | −0.19 | −0.01 | 400 | 0.19 | 0.01 | ||

West | 0.73 | 0.04 | 600 | 0.27 | 0.02 | ||

Northwest | −0.04 | −0.0002 | 800 | 0.50 | 0.03 | ||

>800 | −0.39 | −0.02 |

AUC (area under the curve) is used to evaluate the accuracy of the model. The larger the AUC, the closer the simulated value is to the true value of the sample, and the more accurate the model evaluation result is. It can be seen from

1) This paper uses three models of IOE, IV and WIV to evaluate the susceptibility of landslides in Pingshan County. The susceptibility zoning map of the corresponding model is obtained by the ArcGIS raster overlay, and the accuracy of the three models is 76.2%, 75.6%, and 79.1% after ROC curve inspection. Among them, the susceptibility zoning map generated by the WIV model has the highest accuracy, indicating that this model is more suitable for landslide hazard susceptibility evaluation research in this type of study area.

2) Based on the topography and geomorphology characteristics of the study area, this paper selects 9 evaluation factors, and after correlation analysis, the slope, aspect, NDVI, rainfall, distance to road, distance to fault, distance to water system, and strata are selected. Evaluation factors establish a landslide susceptibility evaluation system. According to the calculation of the model, the factors that account for the higher weight of the evaluation factors are NDVI, distance from the road, and slope. It can be concluded that landslide hazards are mainly concentrated in areas with medium vegetation coverage, short distance from the road, and gentle slope.

3) According to the zoning map of landslide hazard susceptibility, most areas of Taiping Township, Pingshan Town, Fuyan Town, Jinping Town, and Xin’an Town are located in areas with high landslide hazard and extremely high proneness, which is roughly consistent with the actual situation. The proportion of high-density landslide disaster areas is almost the same as that of high-risk areas.

The authors declare no conflicts of interest.

Yang, Y.J. and Zhang, W.J. (2021) Assessment of Landslide Susceptibility Based on Weighted Information Value Model in Pingshan County. Open Access Library Journal, 8: e8199. https://doi.org/10.4236/oalib.1108199