Evaluation and validation of landslide-susceptibility maps obtained by a GIS matrix method: examples from the Betic Cordillera (southern Spain)

This work presents the results of applying the matrix method in a Geographic Information System (GIS) to the drawing of maps of susceptibility to slope movements in different sectors of the Betic Cordillera (southern Spain). In addition, the susceptibility models built by the matrix method were compared with a multivariate statistical method, and the first method gave the best results. The susceptibility maps drawn by the GIS matrix method were validated by calculating the coefficients of association with the degree of fit between recent slope movements registered in 1997 and the different levels of susceptibility of previously drawn maps (1995–1996) in different representative zones of the Betic Cordillera (southern Spain). The first sector studied showed excellent degrees of fit, with an error of less than 10% for all the slope failures and 3% when considering only failures of natural origin. In the second sector, the relative errors were less than 5%. In the third sector, the error hardly exceeded 6%. The results are discussed in the different zones and for each type of slope movement. In any case, these results evidence the predictive capacity of susceptibility maps drawn in GIS by the matrix method, for a great number of slope movements.     

Landslide susceptibility mapping of vicinity of Yaka Landslide (Gelendost,Turkey) using conditional probability approach in GIS

On 19 February 2007, a landslide occurred on the Alaard?ç Slope, located 1.6 km south of the town of Yaka (Gelendost, Turkey.) Subsequently, the displaced materials transformed into a mud flow in E?lence Creek and continued 750 m downstream towards the town of Yaka. The mass poised for motion in the Yaka Landslide source area and its vicinity, which would be triggered to a kinetic state by trigger factors such as heavy or sustained rainfall and/or snowmelt, poise a danger in the form of loss of life and property to Yaka with its population of 3,000. This study was undertaken to construct a susceptibility mapping of the vicinity of the Yaka Landslide’s source area and to relate it to movement of the landslide mass with the goal of prevention or mitigation of loss of life and property. The landslide susceptibility map was formulated by designating the relationship of the effecting factors that cause landslides such as lithology, gradient, slope aspect, elevation, topographical moisture index, and stream power index to the landslide map, as determined by analysis of the terrain, through the implementation of the conditional probability method. It was determined that the surface area of the Goksogut formation, which has attained lithological characteristics of clayey limestone with a broken and separated base and where area landslides occur, possesses an elevation of 1,100–1,300 m, a slope gradient of 15 °–35 ° and a slope aspect between 0 °–67.5 ° and 157 °–247 °. Loss of life and property may be avoided by the construction of structures to check the debris mass in E?lence Creek, the cleaning of the canal which passes through Yaka, the broadening of the canal’s base area, elevating the protective edges along the canal and the establishment of a protective zone at least 10-m wide on each side of the canal to deter against damage from probable landslide occurrence and mud flow.     

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.     

The slope aspect: A predisposing factor for landsliding?

The influence of slope aspect on the distribution of landslides was studied in the Milia and Roglio basins in Tuscany, Italy. For each basin, the new Tuscany region landslide inventory that was initiated in 2010 was used. The landslides were split into separate datasets based on their prevailing movement typology. To assess the results that were obtained from the different slope aspect values, maps of the lithology, slope angle, distances to streams, and distances to tectonic lineaments were included in the bivariate statistical analysis as comparison terms. For each basin, all of the geo-environmental factor maps were compared with the different landslide typologies with GIS software. Pearson's Chi² (χ²) coefficient was used to test the degree of spatial association between each predictor variable and landslide type. In addition, Cramer's V test was used to quantify the strength of the degree of association. Next, a conditional analysis was applied to all of the possible combinations that occurred between the slope aspect and other landslide-predisposing factors. Overall, the slope aspect significantly affected the distribution of superficial landslide types, but apparently not that of other landslide types.     

Landslide susceptibility zonation mapping and its validation in part of Garhwal Lesser Himalaya, India, using binary logistic regression analysis and receiver operating characteristic curve method

A landslide susceptibility zonation (LSZ) map helps to understand the spatial distribution of slope failure probability in an area and hence it is useful for effective landslide hazard mitigation measures. Such maps can be generated using qualitative or quantitative approaches. The present study is an attempt to utilise a multivariate statistical method called binary logistic regression (BLR) analysis for LSZ mapping in part of the Garhwal Lesser Himalaya, India, lying close to the Main Boundary Thrust (MBT). This method gives the freedom to use categorical and continuous predictor variables together in a regression analysis. Geographic Information System has been used for preparing the database on causal factors of slope instability and landslide locations as well as for carrying out the spatial modelling of landslide susceptibility. A forward stepwise logistic regression analysis using maximum likelihood estimation method has been used in the regression. The constant and the coefficients of the predictor variables retained by the regression model have been used to calculate the probability of slope failure for the entire study area. The predictive logistic regression model has been validated by receiver operating characteristic curve analysis, which has given 91.7% accuracy for the developed BLR model.     

Use of satellite remote sensing data in the mapping of global landslide susceptibility

Satellite remote sensing data has significant potential use in analysis of natural hazards such as landslides. Relying on the recent advances in satellite remote sensing and geographic information system (GIS) techniques, this paper aims to map landslide susceptibility over most of the globe using a GIS-based weighted linear combination method. First, six relevant landslide-controlling factors are derived from geospatial remote sensing data and coded into a GIS system. Next, continuous susceptibility values from low to high are assigned to each of the six factors. Second, a continuous scale of a global landslide susceptibility index is derived using GIS weighted linear combination based on each factor’s relative significance to the process of landslide occurrence (e.g., slope is the most important factor, soil types and soil texture are also primary-level parameters, while elevation, land cover types, and drainage density are secondary in importance). Finally, the continuous index map is further classified into six susceptibility categories. Results show the hot spots of landslide-prone regions include the Pacific Rim, the Himalayas and South Asia, Rocky Mountains, Appalachian Mountains, Alps, and parts of the Middle East and Africa. India, China, Nepal, Japan, the USA, and Peru are shown to have landslide-prone areas. This first-cut global landslide susceptibility map forms a starting point to provide a global view of landslide risks and may be used in conjunction with satellite-based precipitation information to potentially detect areas with significant landslide potential due to heavy rainfall.     

Determination and application of the weights for landslide susceptibility mapping using an artificial neural network

The purpose of this study is the development, application, and assessment of probability and artificial neural network methods for assessing landslide susceptibility in a chosen study area. As the basic analysis tool, a Geographic Information System (GIS) was used for spatial data management and manipulation. Landslide locations and landslide-related factors such as slope, curvature, soil texture, soil drainage, effective thickness, wood type, and wood diameter were used for analyzing landslide susceptibility. A probability method was used for calculating the rating of the relative importance of each factor class to landslide occurrence. For calculating the weight of the relative importance of each factor to landslide occurrence, an artificial neural network method was developed. Using these methods, the landslide susceptibility index (LSI) was calculated using the rating and weight, and a landslide susceptibility map was produced using the index. The results of the landslide susceptibility analysis, with and without weights, were confirmed from comparison with the landslide location data. The comparison result with weighting was better than the results without weighting. The calculated weight and rating can be used to landslide susceptibility mapping.     

基于ArcGIS Engine的工程勘察柱状图绘制系统的设计与实现

柱状图的自动绘制作为工程勘察信息系统的一个重要组成部分,在GIS环境下的设计和实现是文章探讨的主题。这里将基于刚刚推出的ArcGIS 9.0全新开发环境ArcGIS Engine,先论述柱状图绘制系统中“模板”的设计以及以“模板”为基础的柱状图绘制系统的总体结构设计,再详细阐述它们的实现方法,最后总结ArcGIS Engine的一些实用技巧。     

Landslide susceptibility mapping in the Damrei Romel area, Cambodia using frequency ratio and logistic regression models

This study applied, tested and compared a probability model, a frequency ratio and statistical model, a logistic regression to Damre Romel area, Cambodia, using a geographic information system. For landslide susceptibility mapping, landslide locations were identified in the study area from interpretation of aerial photographs and field surveys, and a spatial database was constructed from topographic maps, geology and land cover. The factors that influence landslide occurrence, such as slope, aspect, curvature and distance from drainage were calculated from the topographic database. Lithology and distance from lineament were extracted and calculated from the geology database. Land cover was classified from Landsat TM satellite imagery. The relationship between the factors and the landslides was calculated using frequency ratio and logistic regression models. The relationships, frequency ratio and logistic regression coefficient were overlaid to make landslide susceptibility map. Then the landslide susceptibility map was compared with known landslide locations and tested. As the result, the frequency ratio model (86.97%) and the logistic regression (86.37%) had high and similar prediction accuracy. The landslide susceptibility map can be used to reduce hazards associated with landslides and to land cover planning.     

Landslide hazard mapping at Selangor, Malaysia using frequency ratio and logistic regression models

The aim of this study is to evaluate the landslide hazards at Selangor area, Malaysia, using Geographic Information System (GIS) and Remote Sensing. Landslide locations of the study area were identified from aerial photograph interpretation and field survey. Topographical maps, geological data, and satellite images were collected, processed, and constructed into a spatial database in a GIS platform. The factors chosen that influence landslide occurrence were: slope, aspect, curvature, distance from drainage, lithology, distance from lineaments, land cover, vegetation index, and precipitation distribution. Landslide hazardous areas were analyzed and mapped using the landslide-occurrence factors by frequency ratio and logistic regression models. The results of the analysis were verified using the landslide location data and compared with probability model. The comparison results showed that the frequency ratio model (accuracy is 93.04%) is better in prediction than logistic regression (accuracy is 90.34%) model.     

Landslide susceptibility mapping for a landslide-prone area (Findikli,NE of Turkey) by likelihood-frequency ratio and weighted linear combination models

Landslides are very common natural problems in the Black Sea Region of Turkey due to the steep topography, improper use of land cover and adverse climatic conditions for landslides. In the western part of region, many studies have been carried out especially in the last decade for landslide susceptibility mapping using different evaluation methods such as deterministic approach, landslide distribution, qualitative, statistical and distribution-free analyses. The purpose of this study is to produce landslide susceptibility maps of a landslide-prone area (Findikli district, Rize) located at the eastern part of the Black Sea Region of Turkey by likelihood frequency ratio (LRM) model and weighted linear combination (WLC) model and to compare the results obtained. For this purpose, landslide inventory map of the area were prepared for the years of 1983 and 1995 by detailed field surveys and aerial-photography studies. Slope angle, slope aspect, lithology, distance from drainage lines, distance from roads and the land-cover of the study area are considered as the landslide-conditioning parameters. The differences between the susceptibility maps derived by the LRM and the WLC models are relatively minor when broad-based classifications are taken into account. However, the WLC map showed more details but the other map produced by LRM model produced weak results. The reason for this result is considered to be the fact that the majority of pixels in the LRM map have high values than the WLC-derived susceptibility map. In order to validate the two susceptibility maps, both of them were compared with the landslide inventory map. Although the landslides do not exist in the very high susceptibility class of the both maps, 79% of the landslides fall into the high and very high susceptibility zones of the WLC map while this is 49% for the LRM map. This shows that the WLC model exhibited higher performance than the LRM model.     

Landslide Susceptibility Models Utilising Spatial Data Analysis Techniques. A Case Study from the Lower Deba Valley,Guipuzcoa (Spain)

A quantitative methodology for landslide susceptibility zonationis described and its application to a study area in the lower part of the Deba Valley (Guipúzcoa,Spain) presented. Susceptibility models were obtained on the basis of statisticalrelationships between known mass movements and conditioning factors. A landslide rupturehypothesis was set and a digital database consisting of seventeen causal factors layers constructed.The modelling procedure was implemented utilising a GIS. The susceptibility analysis methodis based on the Favourability Functions approach, and two different mathematicalframeworks: probability theory and Zadehïs fuzzy set theory. Several landslidesusceptibility models were produced and validated using different sets of independent landslide data.The predictive capability of models was determined.     

An approach for GIS-based statistical landslide susceptibility zonation—with a case study in the Himalayas

Landslide susceptibility zonation (LSZ) is necessary for disaster management and planning development activities in mountainous regions. A number of methods, viz. landslide distribution, qualitative, statistical and distribution-free analyses have been used for the LSZ studies and they are again briefly reviewed here. In this work, two methods, the Information Value (InfoVal) and the Landslide Nominal Susceptibility Factor (LNSF) methods that are based on bivariate statistical analysis have been applied for LSZ mapping in a part of the Himalayas. Relevant thematic maps representing various factors (e.g., slope, aspect, relative relief, lithology, buffer zones along thrusts, faults and lineaments, drainage density and landcover) that are related to landslide activity, have been generated using remote sensing and GIS techniques. The LSZ derived from the LNSF method, has been compared with that produced from the InfoVal method and the result shows a more realistic LSZ map from the LNSF method which appears to conform to the heterogeneity of the terrain.     

Landslide susceptibility analysis and its verification using likelihood ratio, logistic regression, and artificial neural network models: case study of Youngin, Korea

The likelihood ratio, logistic regression, and artificial neural networks models are applied and verified for analysis of landslide susceptibility in Youngin, Korea, using the geographic information system. From a spatial database containing such data as landslide location, topography, soil, forest, geology, and land use, the 14 landslide-related factors were calculated or extracted. Using these factors, landslide susceptibility indexes were calculated by likelihood ratio, logistic regression, and artificial neural network models. Before the calculation, the study area was divided into two sides (west and east) of equal area, for verification of the models. Thus, the west side was used to assess the landslide susceptibility, and the east side was used to verify the derived susceptibility. The results of the landslide susceptibility analysis were verified using success and prediction rates. The verification results showed satisfactory agreement between the susceptibility map and the existing data on landslide locations.     

Landslide susceptibility mapping using GIS-based weighted linear combination,the case in Tsugawa area of Agano River,Niigata Prefecture,Japan

A spatial database of 791 landslides is analyzed using GIS to map landslide susceptibility in Tsugawa area of Agano River. Data from six landslide-controlling parameters namely lithology, slope gradient, aspect, elevation, and plan and profile curvatures are coded and inserted into the GIS. Later, an index-based approach is adopted both to put the various classes of the six parameters in order of their significance to the process of landsliding and weigh the impact of one parameter against another. Applying primary and secondary-level weights, a continuous scale of numerical indices is obtained with which the study area is divided into five classes of landslide susceptibility. Slope gradient and elevation are found to be important to delineate flatlands that will in no way be subjected to slope failure. The area which is at high scale of susceptibility lies on mid-slope mountains where relatively weak rocks such as sandstone, mudstone and tuff are outcropping as one unit.     

GIS-based weights-of-evidence modelling of rainfall-induced landslides in small catchments for landslide susceptibility mapping

Landslide susceptibility mapping is a vital tool for disaster management and planning development activities in mountainous terrains of tropical and subtropical environments. In this paper, the weights-of-evidence modelling was applied, within a geographical information system (GIS), to derive landslide susceptibility map of two small catchments of Shikoku, Japan. The objective of this paper is to evaluate the importance of weights-of-evidence modelling in the generation of landslide susceptibility maps in relatively small catchments having an area less than 4 sq km. For the study area in Moriyuki and Monnyu catchments, northeast Shikoku Island in west Japan, a data set was generated at scale 1:5,000. Relevant thematic maps representing various factors (e.g. slope, aspect, relief, flow accumulation, soil depth, soil type, land use and distance to road) that are related to landslide activity were generated using field data and GIS techniques. Both catchments have homogeneous geology and only consist of Cretaceous granitic rock. Thus, bedrock geology was not considered in data layering during GIS analysis. Success rates were also estimated to evaluate the accuracy of landslide susceptibility maps and the weights-of-evidence modelling was found useful in landslide susceptibility mapping of small catchments.     

Landslide susceptibility mapping using machine learning algorithms and comparison of their performance at Abha Basin,Asir Region,Saudi Arabia

The current study aimed at evaluating the capabilities of seven advanced machine learning techniques(MLTs),including,Support Vector Machine(SVM),Random Forest(RF),Multivariate Adaptive Regression Spline(MARS),Artificial Neural Network(ANN),Quadratic Discriminant Analysis(QDA),Linear Discriminant Analysis(LDA),and Naive Bayes(NB),for landslide susceptibility modeling and comparison of their performances.Coupling machine learning algorithms with spatial data types for landslide susceptibility mapping is a vitally important issue.This study was carried out using GIS and R open source software at Abha Basin,Asir Region,Saudi Arabia.First,a total of 243 landslide locations were identified at Abha Basin to prepare the landslide inventory map using different data sources.All the landslide areas were randomly separated into two groups with a ratio of 70%for training and 30%for validating purposes.Twelve landslide-variables were generated for landslide susceptibility modeling,which include altitude,lithology,distance to faults,normalized difference vegetation index(NDVI),landuse/landcover(LULC),distance to roads,slope angle,distance to streams,profile curvature,plan curvature,slope length(LS),and slope-aspect.The area under curve(AUC-ROC)approach has been applied to evaluate,validate,and compare the MLTs performance.The results indicated that AUC values for seven MLTs range from 89.0%for QDA to 95.1%for RF.Our findings showed that the RF(AUC=95.1%)and LDA(AUC=941.7%)have produced the best performances in comparison to other MLTs.The outcome of this study and the landslide susceptibility maps would be useful for environmental protection.     

Implementation of reconstructed geomorphologic units in landslide susceptibility mapping: the Melen Gorge (NW Turkey)

In the international literature, although considerable amount of publications on the landslide susceptibility mapping exist, geomorphology as a conditioning factor is still used in limited number of studies. Considering this factor, the purpose of this article paper is to implement the geomorphologic parameters derived by reconstructed topography in landslide susceptibility mapping. According to the method employed in this study, terrain is generalized by the contours passed through the convex slopes of the valleys that were formed by fluvial erosion. Therefore, slope conditions before landsliding can be obtained. The reconstructed morphometric and geomorphologic units are taken into account as a conditioning parameter when assessing landslide susceptibility. Two different data, one of which is obtained from the reconstructed DEM, have been employed to produce two landslide susceptibility maps. The binary logistic regression is used to develop landslide susceptibility maps for the Melen Gorge in the Northwestern part of Turkey. Due to the high correct classification percentages and spatial effectiveness of the maps, the landslide susceptibility map comprised the reconstructed morphometric parameters exhibits a better performance than the other. Five different datasets are selected randomly to apply proper sampling strategy for training. As a consequence of the analyses, the most proper outcomes are obtained from the dataset of the reconstructed topographical parameters and geomorphologic units, and lithological variables that are implemented together. Correct classification percentage and root mean square error (RMSE) values of the validation dataset are calculated as 86.28% and 0.35, respectively. Prediction capacity of the different datasets reveal that the landslide susceptibility map obtained from the reconstructed parameters has a higher prediction capacity than the other. Moreover, the landslide susceptibility map obtained from the reconstructed parameters produces logical results.