Landslide susceptibility mapping by using fuzzy logic: a case study of Cham-gardalan catchment,Ilam, Iran

Landslides are introduced as regional movements, which influence different engineering structures such as roads, railways, and dams and cause the person’s death. Identification of landslide zones may decrease the financial losses and human injuries or deaths. This study tries to achieve a landslide susceptibility mapping in Cham-gardalan catchment by weighting the main criteria and the membership functions of fuzzy logic. For this, we applied the best relationship function between the presence and absence of landslides as well as a collection of the elements. At first, the landslide points were identified by the means of some components those of satellite images, topographical (1:50,000) and geographical (1:100,000) maps, field visits, and Google Earth software followed by the preparation of landslide distribution maps. Then, all effective landslide factors such as percentage of slope, slope aspect, height, geology, land uses, distance from roads, distance from drainages, distance from breakage, and precipitation map have been utilized in order to conduct the fuzzy analyses. Landslide susceptibility map was performed by fuzzy operators (Gamma, Product, Sum, Or, And) in the study area. After fuzzificating and weighting, the effective criteria of landslides were determined through fuzzy Gamma operators with the landaus of 0.2, 0.5, 0.8, and 0.9 and by comparing final maps for making an appropriate model of landslide susceptibility mapping. The regional susceptibility map represents the landslide-prone areas in five categories those of very low, low, moderate, high, and very high. Our results indicated that among the applied operators, Gamma with landau of 0.9 can be used as an appropriate method for mapping the landslide susceptibility due to the suitable fuzzification of given criteria based on landslide distribution maps. In addition, the elements of road, percentage of slope, distance from drainage, and geology were recognized as the most important factors for occurring the landslides.     

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.     

Evaluation of the consistency of landslide susceptibility mapping: a case study from the Kankai watershed in east Nepal

GIS-based landslide susceptibility maps for the Kankai watershed in east Nepal are developed using the frequency ratio method and the multiple linear regression technique. The maps are derived from comparing observed landslides with possible causative factors: slope angle, slope aspect, slope curvature, relative relief, distance from drainage, land use, geology, distance from faults and mean annual rainfall. The consistency of the maps is evaluated using landslide density analysis, success rate analysis and spatially agreed area approach. The first two analyses produce almost identical quantitative results, whereas the last approach is able to reveal spatial differences between the maps and also to improve predictions in the agreed high landslide-susceptible area.     

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.     

Comparison of the implementation of rock engineering system and analytic hierarchy process methods,upon landslide susceptibility mapping,using GIS: a case study from the Eastern Achaia County of Peloponnesus,Greece

As landslides are very common in Greece, causing serious problems to the social and economic welfare of many communities, the implementation of a proper hazard analysis system will help the creation of a reliable susceptibility map. Τhis will help local communities to define a safe land use and urban development. The purpose of this study is to compare the implementation of two semi-quantitative landslide assessment approaches, using landslide susceptibility maps compiled in a GIS environment. The compared methods are rock engineering system (RES) and the analytic hierarchy process (AHP). For the landslide susceptibility analysis, the Northeastern part of the Achaia County was examined. This area suffers from many landslides, because of its neighborhood with the tectonically active Corinthian Gulf and its geological setting (Neogene sediments, flysch and other bedrock formations, with local overthrusts). Ten parameters were used in both methodologies, and each one was separated into five categories ranging from 0 to 4, representing their specific conditions derived from the investigation of the landslides in the western part of the study area (ranking area). A layer map was generated for each parameter, using GIS, while the weighting coefficients of each methodology were used for the compilation of RES and AHP final maps of the eastern part of the study area (validating area). By examining these two maps, it is revealed that even though both correctly show the landslide status of the second site, the RES map reveals a better behavior in the spatial distribution of the various landslide susceptibility zones.     

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.     

基于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在区域滑坡易发性评价与制图中得到了可靠的结果,为滑坡灾害空间预测提供了新的技术支撑。     

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.     

Landslide susceptibility mapping using GIS and digital photogrammetric techniques: a case study from Ardesen (NE-Turkey)

Ardesen is a settlement area which has been significantly damaged by frequent landslides which are caused by severe rainfalls and result in many casualties. In this study a landslide susceptibility map of Ardesen was prepared using the Analytical Hierarchy Process (AHP) with the help of Geographical Information Systems (GIS) and Digital Photogrametry Techniques (DPT). A landslide inventory, lithology–weathering, slope, aspect, land cover, shear strength, distance to the river, stream density and distance to the road thematics data layers were used to create the map. These layer maps are produced using field, laboratory and office studies, and by the use of GIS and DPT. The landslide inventory map is also required to determine the relationship between these maps and landslides using DPT. In the study field in the Hemsindere Formation there are units that have different weathering classes, and this significantly affects the shear strength of the soil. In this study, shear strength values are calculated in great detail with field and laboratory studies and an additional layer is evaluated with the help of the stability studies used to produce the landslide susceptibility map. Finally, an overlay analysis is carried out by evaluating the layers obtained according to their weight, and the landslide susceptibility map is produced. The study area was classified into five classes of relative landslide susceptibility, namely, very low, low, moderate, high, and very high. Based on this analysis, the area and percentage distribution of landslide susceptibility degrees were calculated and it was found that 28% of the region is under the threat of landslides. Furthermore, the landslide susceptibility map and the landslide inventory map were compared to determine whether the models produced are compatible with the real situation resulting in compatibility rate of 84%. The total numbers of dwellings in the study area were determined one by one using aerial photos and it was found that 30% of the houses, with a total occupancy of approximately 2,300 people, have a high or very high risk of being affected by landslides.     

Slope stability analysis on a regional scale using GIS: a case study from Dhading,Nepal

A spatially distributed physically based slope stability model combined with a hydrological model is presented and applied to a 350-km² area located in Dhading district, Nepal. Land slide safety factor maps are generated for five cases, including three steady state conditions assuming either completely dry soils, half saturated soils, or fully saturated soils, and two quasi-dynamic conditions, i.e. soil wetness resulting from storm events with, respectively a 2 or 25-year return period. For the quasi-dynamic cases, two methods are used, one based on accumulation of groundwater flow from upstream areas, and the other on accumulation of soil water from direct infiltration. The methodology delineates areas most prone to shallow land sliding in function of readily available data as topography, land-use and soil types. For the study area only 29% of the soils are unconditionally stable, while 25% of the soils are found to be unstable under fully saturated conditions. The comparison between the methods based on contributing area or on infiltration for quasi-dynamic conditions show that the approach based on infiltration is more reliable for the study area. The proposed methodology for predicting landslide susceptibility on a regional scale, based on basic data in GIS form, may be useful for other remote regions where detailed information is not available.     

Geomorphology and landslide susceptibility assessment using GIS and bivariate statistics: a case study in southern Italy

In this article, the results of a study aimed to assess the landslide susceptibility in the Calaggio Torrent basin (Campanian Apennines, southern Italy) are presented. The landslide susceptibility has been assessed using two bivariate-statistics-based methods in a GIS environment. In the first method, widely used in the existing literature, weighting values (Wi) have been calculated for each class of the selected causal factors (lithology, land-use, slope angle and aspect) taking into account the landslide density (detachment zones + landslide body) within each class. In the second method, which is a modification of the first method, only the landslide detachment zone (LDZ) density has been taken into account to calculate the weighting values. This latter method is probably characterized by a major geomorphological coherence. In fact, differently from the landslide bodies, LDZ must necessarily occur in geoenvironmental classes prone to failure. Thus, the calculated Wi seem to be more reliable in estimating the propensity of a given class to generate failure. The thematic maps have been reclassified on the basis of the calculated Wi and then overlaid, with the purpose to produce landslide susceptibility maps. The used methods converge both in indicating that most part of the study area is characterized by a high–very high landslide susceptibility and in the location and extent of the low-susceptible areas. However, an increase of both the high–very high and moderate–high susceptible areas occurs in using the second method. Both the produced susceptibility maps have been compared with the geomorphological map, highlighting an excellent coherence which is higher using method-2. In both methods, the percentage of each susceptibility class affected by landslides increases with the degree of susceptibility, as expected. However, the percentage at issue in the lowest susceptibility class obtained using method-2, even if low, is higher than that obtained using method-1. This suggests that method-2, notwithstanding its major geomorphological coherence, probably still needs further refinements.     

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.     

Comparing data-driven landslide susceptibility models based on participatory landslide inventory mapping in Purwosari area,Yogyakarta, Java

There are different approaches and techniques for landslide susceptibility mapping. However, no agreement has been reached in both the procedure and the use of specific controlling factors employed in the landslide susceptibility mapping. Each model has its own assumption, and the result may differ from place to place. Different landslide controlling factors and the completeness of landslide inventory may also affect the different result. Incomplete landslide inventory may produce significance error in the interpretation of the relationship between landslide and controlling factor. Comparing landslide susceptibility models using complete inventory is essential in order to identify the most realistic landslide susceptibility approach applied typically in the tropical region Indonesia. Purwosari area, Java, which has total 182 landslides occurred from 1979 to 2011, was selected as study area to evaluate three data-driven landslide susceptibility models, i.e., weight of evidence, logistic regression, and artificial neural network. Landslide in the study area is usually affected by rainfall and anthropogenic activities. The landslide typology consists of shallow translational and rotational slide. The elevation, slope, aspect, plan curvature, profile curvature, stream power index, topographic wetness index, distance to river, land use, and distance to road were selected as landslide controlling factors for the analysis. Considering the accuracy and the precision evaluations, the weight of evidence represents considerably the most realistic prediction capacities (79%) when comparing with the logistic regression (72%) and artificial neural network (71%). The linear model shows more powerful result than the nonlinear models because it fits to the area where complete landslide inventory is available, the landscape is not varied, and the occurence of landslide is evenly distributed to the class of controlling factor.     

A fuzzy set based approach for integration of thematic maps for landslide susceptibility zonation

Spatial prediction of landslides is termed landslide susceptibility zonation (LSZ). In this study, an objective weighting approach based on fuzzy concepts is used for LSZ in a part of the Darjeeling Himalayas. Relevant thematic layers pertaining to landslide causative factors have been generated using remote sensing and geographic information system (GIS) techniques. The membership values for each category of thematic layers have been determined using the cosine amplitude fuzzy similarity method and are used as ratings. The integration of these ratings led to the generation of LSZ map. The integration of different ratings to generate an LSZ map has been performed using a fuzzy gamma operator apart from the arithmetic overlay approach. The process is based on determination of combined rating known as the landslide susceptibility index (LSI) for all the pixels using the fuzzy gamma operator and classification using the success rate curve method to prepare the LSZ map. The results indicate that as the gamma value increases, the accuracy of the LSZ map also increases. It is observed that the LSZ map produced by the fuzzy algebraic sum has reflected a more real situation in terms of landslides in the study area.     

A multi-method approach to study the stability of natural slopes and landslide susceptibility mapping

In this paper, a multi-method approach for the assessment of the stability of natural slopes and landslide hazard mapping applied to the Dakar coastal region is presented. This approach is based on the effective combination of geotechnical field and laboratory works, of GIS, and of mechanical (deterministic and numerical) stability analysis. By using this approach, valuable results were gained regarding instability factors, landslide kinematics, simulation of slope failure and coastal erosion. This led to a thorough assessment and strong reduction in the subjectivity of the slope stability and hazard assessment and to the development of an objective landslide danger map of the SW coast of Dakar. Analysis of the results shows that the slides were influenced by the geotechnical properties of the soil, the weathering, the hydrogeological situation, and the erosion by waves. The landslide susceptibility assessment based on this methodological approach has allowed for an appropriate and adequate consideration of the multiple factors affecting the stability and the optimization of planning and investment for land development in the city.     

Three-dimensional mapping of a landslide using a multi-geophysical approach: the Quesnel Forks landslide

A landslide located on the Quesnel River in British Columbia, Canada is used as a case study to demonstrate the utility of a multi-geophysical approach to subsurface mapping of unstable slopes. Ground penetrating radar (GPR), direct current (DC) resistivity and seismic reflection and refraction surveys were conducted over the landslide and adjacent terrain. Geophysical data were interpreted based on stratigraphic and geomorphologic observations, including the use of digital terrain models (DTMs), and then integrated into a 3-dimensional model. GPR surveys yielded high-resolution data that were correlated with stratigraphic units to a maximum depth of 25 m. DC electrical resistivity offered limited data on specific units but was effective for resolving stratigraphic relationships between units to a maximum depth of 40 m. Seismic surveys were primarily used to obtain unit boundaries up to a depth of >80 m. Surfaces of rupture and separation were successfully identified by GPR and DC electrical resistivity techniques.     

Application of fuzzy logic and analytical hierarchy process (AHP) to landslide susceptibility mapping at Haraz watershed, Iran

The main goal of this study is to produce landslide susceptibility maps of a landslide-prone area (Haraz) in Iran by using both fuzzy logic and analytical hierarchy process (AHP) models. At first, landslide locations were identified by aerial photographs and field surveys, and a total of 78 landslides were mapped from various sources. Then, the landslide inventory was randomly split into a training dataset 70?% (55 landslides) for training the models and the remaining 30?% (23 landslides) was used for validation purpose. Twelve data layers, as the landslide conditioning factors, are exploited to detect the most susceptible areas. These factors are slope degree, aspect, plan curvature, altitude, lithology, land use, distance from rivers, distance from roads, distance from faults, stream power index, slope length, and topographic wetness index. Subsequently, landslide susceptibility maps were produced using fuzzy logic and AHP models. For verification, receiver operating characteristics curve and area under the curve approaches were used. The verification results showed that the fuzzy logic model (89.7?%) performed better than AHP (81.1?%) model for the study area. The produced susceptibility maps can be used for general land use planning and hazard mitigation purpose.     

Landslide susceptibility mapping using LiDAR and DMC data: a case study in the Three Gorges area, China

The objective of this study is to map landslide susceptibility in Zigui segment of the Yangtze Three Gorges area that is known as one of the most landslide-prone areas in China by using data from light detection and ranging (LiDAR) and digital mapping camera (DMC). The likelihood ratio (LR) and logistic regression model (LRM) were used in this study. The work is divided into three phases. The first phase consists of data processing and analysis. In this phase, LiDAR and DMC data and geological maps were processed, and the landslide-controlling factors were derived such as landslide density, digital elevation model (DEM), slope angle, aspect, lithology, land use and distance from drainage. Among these, the landslide inventories, land use and drainage were constructed with both LiDAR and DMC data; DEM, slope angle and aspect were constructed with LiDAR data; lithology was taken from the 1:250,000 scale geological maps. The second phase is the logistic regression analysis. In this phase, the LR was applied to find the correlation between the landslide locations and the landslide-controlling factors, whereas the LRM was used to predict the occurrence of landslides based on six factors. To calculate the coefficients of LRM, 13,290,553 pixels was used, 29.5 % of the total pixels. The logical regression coefficients of landslide-controlling factors were obtained by logical regression analysis with SPSS 17.0 software. The accuracy of the LRM was 88.8 % on the whole. The third phase is landslide susceptibility mapping and verification. The mapping result was verified using the landslide location data, and 64.4 % landslide pixels distributed in “extremely high” zone and “high” zone; in addition, verification was performed using a success rate curve. The verification result show clearly that landslide susceptibility zones were in close agreement with actual landslide areas in the field. It is also shown that the factors that were applied in this study are appropriate; lithology, elevation and distance from drainage are primary factors for the landslide susceptibility mapping in the area, while slope angle, aspect and land use are secondary.