Evaluation and comparison of GIS based landslide susceptibility mapping procedures in Kulekhani watershed, Nepal

The goal of this paper is to evaluate and compare the consistency of GIS-based heuristic and bivariate landslide susceptibility mapping techniques in the Himalayan region, taking the Kulekhani watershed of central Nepal as an example. For this purpose, a heuristic and two statistical bivariate landslide susceptibility mapping methods are applied, and three separate landslide susceptibility zonation maps are produced. The maps are compared using three approaches: landslide density analysis, success rate analysis, and agreed area analysis. A comparison of the values obtained from landslide density analysis and the curves of success rate analysis indicate that the two bivariate methods produce almost identical results, whereas the map produced with the heuristic method differs significantly from the others. On the other hand, the agreed area analysis highlights significant spatial differences in the maps obtained from the three methods. Although the three approaches evaluate the consistency of susceptibility maps, only the agreed area analysis is capable of spatially comparing them. Hence, this approach proves to be more suitable for spatially and quantitatively evaluating the consistency of various landslide susceptibility zonation maps.     

GIS based landslide susceptibility mapping using a fuzzy logic approach: A case study from Ghurmi-Dhad Khola area, Eastern Nepal

Landslides cause extensive loss of life and property in the Nepal Himalaya. Since the late 1980s, different mathematical models have been developed and applied for landslide susceptibility mapping and hazard assessment in Nepal. The main goal of this paper is to apply fuzzy logic to landslide susceptibility mapping in the Ghurmi-Dhad Khola area, Eastern Nepal. Seven causative factors are considered: slope angle, slope aspect, distance from drainage, land use, geology, distance from faults and folds, soil and rock type. Likelihood ratios are obtained for each class of causative factors by comparison with past landslide occurrences. The ratios are normalized between zero and one to obtain fuzzy membership values. Further, different fuzzy operators are applied to generate landslide susceptibility maps. Comparison with the landslide inventory map reveals that the fuzzy gamma operator with a γ-value of 0.60 yields the best prediction accuracy. Consequently, this operator is used to produce the final landslide susceptibility zonation map.     

PFR model and GiT for landslide susceptibility mapping: a case study from Central Alborz, Iran

In northern parts of Iran such as the Alborz Mountain belt, frequent landslides occur due to a combination of climate and geologic conditions with high tectonic activities. This results in millions of dollars of financial damages annually excluding casualties and unrecoverable resources. This paper evaluates the landslide susceptible areas in Central Alborz using the probabilistic frequency ratio (PFR) model and Geo-information Technology (GiT). The landslide location map in this study has been generated based on image elements interpreted from IRS satellite data and field observations. The display, manipulation and analysis have been carried out to evaluate layers such as geology, geomorphology, soil, slope, aspect, land use, distance from faults, lineaments, roads and drainages. The validation group of actual landslides and relative operation curve method has been used to increase the accuracy of the final landslide susceptibility map. The area under the curve evaluates how well the method predicts landslides. The results showed a satisfactory agreement of 91% between prepared susceptibility map and existing data on landslide locations.     

Landslide susceptibility mapping using the weight of evidence method in the Tinau watershed, Nepal

Mountainous areas in Nepal are prone to landslides, resulting in an enormous loss of life and property every year. As a first step towards mitigating or controlling such problems, it is necessary to prepare landslide susceptibility maps. Various methodologies have been proposed for landslide susceptibility mapping. This study applies the weight of evidence method to the Tinau watershed in west Nepal. A landslide susceptibility map is prepared on the basis of field observations and available data of geology, land use, topography and hydrology. Predicted susceptibility levels are found to be in good agreement with the locations of past landslides. The results show that about 30?% of the area is highly susceptible to landsliding. The present results provide useful information to the authorities concerning the landslide susceptibility zones and possible improvements for disaster management activities and sustainable development.     

Regional-scale landslide activity and landslide susceptibility zonation in the Nepal Himalaya

Landslide susceptibility zonation mapping is a fundamental procedure for geo-disaster management in tropical and sub-tropical regions. Recently, various landslide susceptibility zonation models have been introduced in Nepal with diverse approaches of assessment. However, validation is still a problem. Additionally, the role of various predisposing causative parameters for landslide activity is still not well understood in the Nepal Himalaya. To address these issues of susceptibility zonation and landslide activity, about 4,000 km² area of central Nepal was selected for regional-scale assessment of landslide activity and susceptibility zonation mapping. In total, 655 new landslides and 9,229 old landslides were identified with the study area with the help of satellite images, aerial photographs, field data and available reports. The old landslide inventory was “blind landslide database” and could not explain the particular rainfall event responsible for the particular landslide. But considering size of the landslide, blind landslide inventory was reclassified into two databases: short-duration high-intensity rainfall-induced landslide inventory and long-duration low-intensity rainfall-induced landslide inventory. These landslide inventory maps were considered as proxy maps of multiple rainfall event-based landslide inventories. Similarly, all 9,884 landslides were considered for the activity assessment of predisposing causative parameters. For the Nepal Himalaya, slope, slope aspect, geology and road construction activity (anthropogenic cause) were identified as most affective predisposing causative parameters for landslide activity. For susceptibility zonation, multivariate approach was considered and two proxy rainfall event-based landslide databases were used for the logistic regression modelling, while a relatively recent landslide database was used in validation. Two event-based susceptibility zonation maps were merged and rectified to prepare the final susceptibility zonation map and its prediction rate was found to be more than 82 %. From this work, it is concluded that rectification of susceptibility zonation map is very appropriate and reliable. The results of this research contribute to a significant improvement in landslide inventory preparation procedure, susceptibility zonation mapping approaches as well as role of various predisposing causative parameters for the landslide activity.     

基于XGBoost模型的三峡库区燕山乡滑坡易发性评价与区划

滑坡易发性评价是精细化滑坡灾害风险评价的基础。为了提升滑坡易发性评价模型的精度和稳健性,以三峡库区万州区燕山乡为例,选取工程地质岩组、堆积层厚度等九个影响因子构建滑坡易发性评价指标体系,应用信息量模型定量分析滑坡发育与指标之间的关系。在此基础上,随机选取70%/30%的滑坡样本作为训练/验证数据集,应用极致梯度提升模型(extreme gradient boosting, XGBoost)开展易发性评价。随后从模型预测精度和模型稳定性两方面将其与决策树模型(decision tree, DT)和梯度提升树模型(gradient boosting decision tree, GBDT)进行对比。结果表明：研究区堆积层滑坡主要受长江水系、堆积层厚度和工程地质岩组影响。XGBoost模型具有最高的准确率(94.3%)和预测精度(97.3%)。在模型稳定性验证中,平均预测精度最高(97.3%),优于DT(91.3%)和GBDT(95.7%),模型标准差和变异系数均为0.01,低于其余两种模型。XGBoost在区域滑坡易发性评价与制图中得到了可靠的结果,为滑坡灾害空间预测提供了新的技术支撑。     

Landslide susceptibility mapping and factor effect analysis using frequency ratio in a catchment scale: a case study from Garuwa sub-basin,East Nepal

For assessing landslide susceptibility, the spatial distribution of landslides in the field is essential. The landslide inventory map is prepared on the basis of historical information of individual landslide events from different sources such as previously published reports, satellite imageries, aerial photographs and interview with local inhabitants. Then, the distribution of landslides in the study area is verified with field surveys. However, the selection of contributing factors for modelling landslide susceptibility is an inhibit task. The previous studies show that the factors are chosen as per availability of data. This paper documents the landslide susceptibility mapping in the Garuwa sub-basin, East Nepal using frequency ratio method. Nine different contributing factors are considered: slope aspect, slope angle, slope shape, relative relief, geology, distance from faults, land use, distance from drainage and annual rainfall. To analyse the effect of contributing factors, the landslide susceptibility index maps are generated four times using (a) topographical factors and geological factors, (b) topographical factors, geological factors and land use, (c) topographical factors, geological factors, land use and drainage and (d) all nine causative factors. By comparing with the pre-existing landslides, the fourth case (considering all nine causative factors) yields the best success rate accuracy, i.e. 81.19 %, which is then used to produce the final landslide susceptibility zonation map. Then, the final landslide susceptibility map is validated through chi-square test. The standard chi-square value with 3 degrees of freedom at the 0.001 significance level is 16.3, whereas the calculated chi-square value is 7,125.79. Since the calculated chi-square value is greater than the standard chi-square value, it can be concluded that the landslide susceptibility map is considered as statistically significant. Moreover, the results show that the predicted susceptibility levels are found to be in good agreement with the past landslide occurrences.     

Frequency ratio model based landslide susceptibility mapping in lower Mae Chaem watershed, Northern Thailand

The purpose of this study is to produce a landslide susceptibility map for the lower Mae Chaem watershed, northern Thailand using a Geographic Information System (GIS) and remotely sensed images. For this purpose, past landslide locations were identified from satellite images and aerial photographs accompanied by the field surveys to create a landslide inventory map. Ten landslide-inducing factors were used in the susceptibility analysis: elevation, slope angle, slope aspect, lithology, distance from lineament, distance from drainage, precipitation, soil texture, land use/land cover (LULC), and NDVI. The first eight factors were prepared from their associated database while LULC and NDVI maps were generated from Landsat-5 TM images. Landslide susceptibility was analyzed and mapped using the frequency ratio (FR) model that determines the level of correlation between locations of past landslides and the chosen factors and describes it in terms of frequency ratio index. Finally, the output map was validated using the area under the curve (AUC) method where the success rate of 80.06% and the prediction rate of 84.82% were achieved. The obtained map can be used to reduce landslide hazard and assist with proper planning of LULC in the future.     

Assessing LNRF,FR, and AHP models in landslide susceptibility mapping index: a comparative study of Nojian watershed in Lorestan province,Iran

Landslides and slope instabilities are major risks for human activities which often lead to economic losses and human fatalities all over the world. The main purpose of this study is to evaluate and compare the results of Landslide Nominal Risk Factor (LNRF), Frequency Ratio (FR), and Analytical Hierarchy Process (AHP) models in mapping Landslide Susceptibility Index (LSI). The study case, Nojian watershed with an area of 344.91 km², is located in Lorestan province of Iran. The procedure was as follows: first, the effective factors of the landslide basin were prepared for each layer in the GIS software. Then, the layers and the landslides of the basin were also prepared using aerial photographs, satellite images, and fieldwork. Next, the effective factors of the layers were overlapped with the map of landslide distribution to specify the role of units in such distribution. Finally, nine factors including lithology, slope, aspect, altitude, distance from the fault, distance from river, fault land use, rainfall, and altitude were found to be effective elements in landslide occurrence of the basin. The final maps of LSI were prepared based on seven factors using LNRF, FR, and AHP models in GIS. The index of the quality sum (Qs) was also used to assess the accuracy of the LSI maps. The results of the three models with LNRF (40%), FR (39%), and AHP (44%) indicated that the whole study area was located in the classes of high to very high hazard. The Qs values for the three models above were also found to be 0.51, 0.70 and 0.70, respectively. In comparison, according to the amount of Qs, the results of AHP and FR models have slightly better performed than the LNRF model in determining the LSI maps in the study area. Finally, the study watershed was classified into five classes based on LSI as very low, low, moderate, high, and very high. The landslide susceptibility maps can be helpful to select sites and mitigate landslide hazards in the study area and the regions with similar conditions.     

Landslide susceptibility maps comparing frequency ratio and artificial neural networks: a case study from the Nepal Himalaya

This study considers landslide susceptibility mapping by means of frequency ratio and artificial neural network approaches using geographic information system (GIS) techniques as a basic analysis tool. The selected study area was that of the Panchthar district, Nepal. GIS was used for the management and manipulation of spatial data. Landslide locations were identified from field survey and aerial photographic interpretation was used for location of lineaments. Ten factors in total are related to the occurrence of landslides. Based on the same set of factors, landslide susceptibility maps were produced from frequency ratio and neural network models, and were then compared and evaluated. The weights of each factor were determined using the back-propagation training method. Landslide susceptibility maps were produced from frequency ratio and neural network models, and they were then compared by means of their checking. The landslide location data were used for checking the results with the landslide susceptibility maps. The accuracy of the landslide susceptibility maps produced by the frequency ratio and neural networks is 82.21 and 78.25%, respectively.     

Neural networks and landslide susceptibility: a case study of the urban area of Potenza

For those working in the field of landslide prevention, the estimation of hazard levels and the consequent production of thematic maps are principal objectives. They are achieved through careful analytical studies of the characteristics of landslide prone areas, thus, providing useful information regarding possible future phenomena. Such maps represent a fundamental step in the drawing up of adequate measures of landslide hazard mitigation. However, for a complete estimation of landslide hazard, meant as the degree of probability that a landslide occurs in a given area, within a given space of time, detailed and uniformly distributed data regarding their incidence and causes are required. This information, while obtainable through laborious historical research, is usually partial, incomplete and uneven, and hence, unsatisfactory for zoning on a regional scale. In order to carry this out effectively, the utilization of spatial estimation of the relative levels of landslide hazard in the various areas was considered opportune. These areas were classified according to their levels of proneness to landslide activity without taking recurrence periods into account. Various techniques were developed in order to obtain upheaval numerical estimates. The method used in this study, which was applied in the area of Potenza, is based on techniques derived from artificial intelligence (Artificial Neural Network—ANN). This method requires the definition of appropriate thematic layers, which parameterize the area under study. These are recognized by means of specific analyses in a functional relationship to the event itself. The parameters adopted are: slope gradient, slope aspect, topographical index, topographical shape, elevation, land use and lithology.     

A new approach to use AHP in landslide susceptibility mapping: a case study at Yenice (Karabuk, NW Turkey)

This study aimed to investigate the parameter effects in preparing landslide susceptibility maps with a data-driven approach and to adapt this approach to analytical hierarchy process (AHP). For this purpose, at the first stage, landslide inventory of an area located in the Western Black Sea region of Turkey covering approximately 567?km² was prepared, and a total of 101 landslides were mapped. In order to assess the landslide susceptibility, a total of 13 parameters were considered as the input parameters: slope, aspect, plan curvature, topographical elevation, vegetation cover index, land use, distance to drainage, distance to roads, distance to structural elements, distance to ridges, stream power index, sediment transport capacity index, and wetness index. AHP was selected as the major assessment methodology since the adapted approach and AHP work in data pairs. Adapted to AHP, a similarity relation?Cbased approach, namely landslide relation indicator (LRI) for parameter selection method, was also proposed. AHP and parametric effect analyses were performed by the proposed approach, and seven landslide susceptibility maps were produced. Among these maps, the best performance was gathered from the landslide susceptibility map produced by 9 parameter combinations using area under curve (AUC) approach. For this map, the AUC value was calculated as 0.797, while the others ranged between 0.686 and 0.771. According to this map, 38.3?% of the study area was classified as having very low, 8.5?% as low, 15.0?% as moderate, 20.3?% as high, and 17.9?% as very high landslide susceptibility, respectively. Based on the overall assessments, the proposed approach in this study was concluded as objective and applicable and yielded reasonable results.     

Accuracy assessment of InSAR derived input maps for landslide susceptibility analysis: a case study from the Swiss Alps

In recent years SAR interferometry has become a widely used technique for measuring altitude and displacement of the surface of the earth. Both these capabilities are highly relevant for landslide susceptibility studies. Although there are many problems that make the use of SAR interferometry less suitable for landslide inventory mapping, it’s use in landslide monitoring and in the generation of input maps for landslide susceptibility assessment looks very promising. The present work attempts to evaluate the usefulness and limitations of this technique based on a case study in the Swiss Alps. Input maps were generated from ERS repeat pass data using SAR interferometry. A land cover map has been generated by image classification of multi-temporal SAR intensity images. An InSAR DEM was generated and a number of maps were derived from it, such as slope-, aspect, altitude- and slope form classes. These maps were used to generate landslide and rockfall susceptibility maps, which give fairly well acceptable results. However, a comparison of the InSAR DEM with the conventional Swisstopo DEM, indicated significant errors in the absolute height and slope angles derived from InSAR, especially along the ridges and in the valleys. These errors are caused by low coherence mostly due to layover and shadow effects. Visual comparison of stereo images created from hillshading maps and corresponding DEMs demonstrate that a considerable amount of topographic details have been lost in the InSAR-derived DEM. It is concluded that InSAR derived input maps are not ideal for landslide susceptibility assessment, but could be used if more accurate data is lacking.     

Evaluation of deep learning algorithms for national scale landslide susceptibility mapping of Iran

The identification of landslide-prone areas is an essential step in landslide hazard assessment and mitigation of landslide-related losses.In this study,we applied two novel deep learning algorithms,the recurrent neural network(RNN)and convolutional neural network(CNN),for national-scale landslide susceptibility mapping of Iran.We prepared a dataset comprising 4069 historical landslide locations and 11 conditioning factors(altitude,slope degree,profile curvature,distance to river,aspect,plan curvature,distance to road,distance to fault,rainfall,geology and land-sue)to construct a geospatial database and divided the data into the training and the testing dataset.We then developed RNN and CNN algorithms to generate landslide susceptibility maps of Iran using the training dataset.We calculated the receiver operating characteristic(ROC)curve and used the area under the curve(AUC)for the quantitative evaluation of the landslide susceptibility maps using the testing dataset.Better performance in both the training and testing phases was provided by the RNN algorithm(AUC=0.88)than by the CNN algorithm(AUC=0.85).Finally,we calculated areas of susceptibility for each province and found that 6%and 14%of the land area of Iran is very highly and highly susceptible to future landslide events,respectively,with the highest susceptibility in Chaharmahal and Bakhtiari Province(33.8%).About 31%of cities of Iran are located in areas with high and very high landslide susceptibility.The results of the present study will be useful for the development of landslide hazard mitigation strategies.     

Statistical evaluation of the effect of earthquake with other related factors on landslide susceptibility: using the watershed area of Shihmen reservoir in Taiwan as a case study

Landslide susceptibility evaluation is one of the most important issues in watershed management. After an earthquake, the landslide susceptibility decreases functionally with increases in the distance from the epicenter. Under the same rainfall intensity, landslides are more likely to occur in an area where earthquakes occur more frequently. However, the questions of how much an earthquake should be weighted and how to evaluate the effects of an earthquake still need to be studied. To understand how earthquakes affect rainfall-triggered landslides, the horizontal peak ground acceleration (PGA) data from the Central Weather Bureau Seismic Network is used as the earthquake factor and combined with other factors to determine the weight of earthquakes in landslide susceptibility using logistic regression. The results indicate that the ability of landslide prediction is better when considering the earthquake factor. This study also proved that although there are no co-seismic landslides (after earthquakes) in the study area, the earthquake factor is still required to increase the model accuracy. PGA has been described as a usable factor. In areas with frequent earthquakes and high geological activity, when using historical data to evaluate landslide susceptibility, the earthquake factor should be taken into consideration to prevent errors.     

Landslide susceptibility mapping using GIS and weighted overlay method: a case study from North Jordan

Road instability along the Jerash–Amman highway was assessed using the weighted overlay method in Geographic Information System environment. The landslide susceptibility map was developed from nine contributing parameters. The map of landslide susceptibility was classified into five zones: very low (very stable), low (stable), moderate (moderately stable), high (unstable), and very high (highly unstable). The very high susceptibility and high susceptibility zones covered 15.14% and 31.81% of the study area, respectively. The main factors that made most parts of study area prone to landslides include excessive drainage channels, road cuts, and unfavorable rock strata such as marl and friable sandstone intercalated with clay and highly fractured limestone. Fracture zones are a major player in land instability. The moderate and high susceptibility zones are the most common in urban (e.g., Salhoub and Gaza camp) and agricultural areas. About 34% of the urban areas and 28.82% of the agricultural areas are characterized by the high susceptibility zone. Twenty percent of the Jerash–Amman highway length and 58% of the overall highway length are located in the very high susceptibility zone. The landslide susceptibility map was validated by the recorded landslides. More than 80 of the inventoried landslides are in unstable zones, which indicate that the selected causative factors are relevant and the model performs properly.