A comparison of the GIS based landslide susceptibility assessment methods: multivariate versus bivariate

The purpose of this study is to evaluate and to compare the results of multivariate (logical regression) and bivariate (landslide susceptibility) methods in Geographical Information System (GIS) based landslide susceptibility assessment procedures. In order to achieve this goal the Asarsuyu catchment in NW Turkey was selected as a test zone because of its well-known landslide occurrences interfering with the E-5 highway mountain pass.Two methods were applied to the test zone and two separate susceptibility maps were produced. Following this a two-fold comparison scheme was implemented. Both methods were compared by the Seed Cell Area Indexes (SCAI) and by the spatial locations of the resultant susceptibility pixels.It was found that both of the methods converge in 80% of the area; however, the weighting algorithm in the bivariate technique (landslide susceptibility method) had some severe deficiencies, as the resultant hazard classes in overweighed areas did not converge with the factual landslide inventory map. The result of the multivariate technique (logical regression) was more sensitive to the different local features of the test zone and it resulted in more accurate and homogeneous susceptibility maps.     

GIS-based landslide susceptibility mapping using bivariate statistical analysis in Devrek (Zonguldak-Turkey)

Devrek town with increasing population is located in a hillslope area where some landslides exist. Therefore, landslide susceptibility map of the area is required. The purpose of this study was to generate a landslide susceptibility map using a bivariate statistical index and evaluate and compare the results of the statistical analysis conducted with three different approaches in seed cell concept resulting in different data sets in Geographical Information Systems (GIS) based landslide susceptibility mapping applied to the Devrek region. The data sets are created from the seed cells of (a) crowns and flanks, (b) only crowns, and (c) only flanks of the landslides by using ten different causative parameters of the study area. To increase the data dependency of the analysis, all parameter maps are classified into equal frequency classes based directly on the percentile divisions of each corresponding seed cell data set. The resultant maps of the landslide susceptibility analysis indicate that all data sets produce fairly acceptable results. In each data set analysis, elevation, lithology, slope, aspect, and drainage density parameters are found to be the most contributing factors in landslide occurrences. The results of the three data sets are compared using Seed Cell Area Indexes (SCAI). This comparison shows that the crown data set produces the most accurate and successful landslide susceptibility map of the study area.     

A GIS-based landslide susceptibility evaluation using bivariate and multivariate statistical analyses

Bivariate and multivariate statistical analyses were used to predict the spatial distribution of landslides in the Cuyahoga River watershed, northeastern Ohio, U.S.A. The relationship between landslides and various instability factors contributing to their occurrence was evaluated using a Geographic Information System (GIS) based investigation. A landslide inventory map was prepared using landslide locations identified from aerial photographs, field checks, and existing literature. Instability factors such as slope angle, soil type, soil erodibility, soil liquidity index, landcover pattern, precipitation, and proximity to stream, responsible for the occurrence of landslides, were imported as raster data layers in ArcGIS, and ranked using a numerical scale corresponding to the physical conditions of the region. In order to investigate the role of each instability factor in controlling the spatial distribution of landslides, both bivariate and multivariate models were used to analyze the digital dataset. The logistic regression approach was used in the multivariate model analysis. Both models helped produce landslide susceptibility maps and the suitability of each model was evaluated by the area under the curve method, and by comparing the maps with the known landslide locations. The multivariate logistic regression model was found to be the better model in predicting landslide susceptibility of this area. The logistic regression model produced a landslide susceptibility map at a scale of 1:24,000 that classified susceptibility into four categories: low, moderate, high, and very high. The results also indicated that slope angle, proximity to stream, soil erodibility, and soil type were statistically significant in controlling the slope movement.     

Flood susceptibility mapping using integrated bivariate and multivariate statistical models

Flooding can have catastrophic effects on human lives and livelihoods and thus comprehensive flood management is needed. Such management requires information on the hydrologic, geotechnical, environmental, social, and economic aspects of flooding. The number of flood events that took place in Busan, South Korea, in 2009 exceeded the normal situation for that city. Mapping the susceptible areas helps us to understand flood trends and can aid in appropriate planning and flood prevention. In this study, a combination of bivariate probability analysis and multivariate logistic regression was used to produce flood susceptibility maps of Busan City. The main aim of this research was to overcome the weakness of logistic regression regarding bivariate probability capabilities. A flood inventory map with a total of 160 flood locations was extracted from various sources. Then, the flood inventory was randomly split into a testing dataset 70 % for training the models and the remaining 30 %, which was used for validation. Independent variables datasets included the rainfall, digital elevation model, slope, curvature, geology, green farmland, rivers, slope, soil drainage, soil effect, soil texture, stream power index, timber age, timber density, timber diameter, and timber type. The impact of each independent variable on flooding was evaluated by analyzing each independent variable with the dependent flood layer. The validation dataset, which was not used for model generation, was used to evaluate the flood susceptibility map using the prediction rate method. The results of the accuracy assessment showed a success rate of 92.7 % and a prediction rate of 82.3 %.     

Comparing predictive capability of statistical and deterministic methods for landslide susceptibility mapping: a case study in the northern Apennines (Reggio Emilia Province, Italy)

Statistical and deterministic methods are widely used in geographic information system based landslide susceptibility mapping. This paper compares the predictive capability of three different models, namely the Weight of Evidence, the Fuzzy Logic and SHALSTAB, for producing shallow earth slide susceptibility maps, to be included as informative layers in land use planning at a local level. The test site is an area of about 450 km² in the northern Apennines of Italy where, in April 2004, rainfall combined with snowmelt triggered hundreds of shallow earth slides that damaged roads and other infrastructure. An inventory of the landslides triggered by the event was obtained from interpretation of aerial photos dating back to May 2004. The pre-existence of mapped landslides was then checked using earlier aerial photo coverage. All the predictive models were run on the same set of geo-environmental causal factors: soil type, soil thickness, land cover, possibility of deep drainage through the bedrock, slope angle, and upslope contributing area. Model performance was assessed using a threshold-independent approach (the ROC plot). Results show that global accuracy is as high as 0.77 for both statistical models, while it is only 0.56 for SHALSTAB. Besides the limited quality of input data over large areas, the relatively poorer performance of the deterministic model maybe also due to the simplified assumptions behind the hydrological component (steady-state slope parallel flow), which can be considered unsuitable for describing the hydrologic behavior of clay slopes, that are widespread in the study area.     

Validation and evaluation of two multivariate statistical models for predictive shallow landslide susceptibility mapping of the Eastern Pyrenees (Spain)

This paper deals with the quality of two multivariate statistical models based on the Geographical Information System for shallow landslide susceptibility assessment in a test area at La Pobla de Lillet (Eastern Pyrenees, Spain). The quality, which was guaranteed by a rigorous methodology based on a suitable diagnosis, validation, and evaluation of the models, ensured a reliable contrast of the final susceptibility maps. This enables us to transfer the best results to the end user. Landslide susceptibility models were carried out by logistic regression and discriminant analysis of the significant conditioning factors related to the characteristics of the slope and the upslope contributing area captured from the digital elevation model and landslide distribution. The explanatory variables were tested (KS test, principal components and one-way and T-test) to select the most statistically significant ones before being introduced into the logistic and discriminant analyses. Accuracy statistics and the receiver operating characteristic curve used for diagnosis and validation showed similar prediction skills and a good fit to the data with more than 85% of unfailed cells properly classified for the two models. The evaluation of the study area and the correlation function (R ² = 0.83) between the models revealed that the discriminant model overestimated the susceptibility of the most stable zones with respect to the logistic model. Different methods of producing susceptibility maps showed marked differences in matching the models. Substantial spatial agreement (Kappa = 0.741) between binary maps produced by the standard cut-off value descended moderately (Kappa = 0.540) as a result of superimposing maps with five susceptibility levels defined by landslide percentage. Despite the fact that the two statistical models are similar in assessing susceptibility in the study area, the implications for hazard and risk management can be different because of the conservative nature of the discriminant model.     

Improvement of statistical landslide susceptibility mapping by using spatial and global regression methods in the case of More and Romsdal (Norway)

Statistical models are one of the most preferred methods among many landslide susceptibility assessment methods. As landslide occurrences and influencing factors have spatial variations, global models like neural network or logistic regression (LR) ignore spatial dependence or autocorrelation characteristics of data between the observations in susceptibility assessment. However, to assess the probability of landslide within a specified period of time and within a given area, it is important to understand the spatial correlation between landslide occurrences and influencing factors. By including these relations, the predictive ability of the developed model increases. In this respect, spatial regression (SR) and geographically weighted regression (GWR) techniques, which consider spatial variability in the parameters, are proposed in this study for landslide hazard assessment to provide better realistic representations of landslide susceptibility. The proposed model was implemented to a case study area from More and Romsdal region of Norway. Topographic (morphometric) parameters (slope angle, slope aspect, curvature, plan, and profile curvatures), geological parameters (geological formations, tectonic uplift, and lineaments), land cover parameter (vegetation coverage), and triggering factor (precipitation) were considered as landslide influencing factors. These influencing factors together with past rock avalanche inventory in the study region were considered to obtain landslide susceptibility maps by using SR and LR models. The comparisons of susceptibility maps obtained from SR and LR show that SR models have higher predictive performance. In addition, the performances of SR and LR models at the local scale were investigated by finding the differences between GWR and SR and GWR and LR maps. These maps which can be named as comparison maps help to understand how the models estimate the coefficients at local scale. In this way, the regions where SR and LR models over or under estimate the landslide hazard potential were identified.     

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.     

A comparative study on the landslide susceptibility mapping using logistic regression and statistical index models

The logistic regression and statistical index models are applied and verified for landslide susceptibility mapping in Daguan County, Yunnan Province, China, by means of the geographic information system (GIS). A detailed landslide inventory map was prepared by literatures, aerial photographs, and supported by field works. Fifteen landslide-conditioning factors were considered: slope angle, slope aspect, curvature, plan curvature, profile curvature, altitude, STI, SPI, and TWI were derived from digital elevation model; NDVI was extracted from Landsat ETM7; rainfall was obtained from local rainfall data; distance to faults, distance to roads, and distance to rivers were created from a 1:25,000 scale topographic map; the lithology was extracted from geological map. Using these factors, the landslide susceptibility maps were prepared by LR and SI models. The accuracy of the results was verified by using existing landslide locations. The statistical index model had a predictive rate of 81.02%, which is more accurate prediction in comparison with logistic regression model (80.29%). The models can be used to land-use planning in the study area.     

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.     

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.     

A statistical approach to landslide risk modelling at basin scale: from landslide susceptibility to quantitative risk assessment

A procedure for landslide risk assessment is presented. The underlying hypothesis is that statistical relationships between past landslide occurrences and conditioning variables can be used to develop landslide susceptibility, hazard and risk models. The latter require also data on past damages. Landslides occurred during the last 50 years and subsequent damages were analysed. Landslide susceptibility models were obtained by means of Spatial Data Analysis techniques and independently validated. Scenarios defined on the basis of past landslide frequency and magnitude were used to transform susceptibility into quantitative hazard models. To assess vulnerability, a detailed inventory of exposed elements (infrastructures, buildings, land resources) was carried out. Vulnerability values (0–1) were obtained by comparing damages experienced in the past by each type of element with its actual value. Quantitative risk models, with a monetary meaning, were obtained for each element by integrating landslide hazard and vulnerability models. Landslide risk models showing the expected losses for the next 50 years were thus obtained for the different scenarios. Risk values obtained are not precise predictions of future losses but rather a means to identify areas where damages are likely to be greater and require priority for mitigation actions.     

Development of hybrid wavelet packet-statistical models (WP-SM) for landslide susceptibility mapping

Landslides are recognized as one of the most important natural hazards in many areas throughout the world. Producing landslide susceptibility maps have received particular attention from a wide range of scientists. The main objective of this study was to produce landslide susceptibility maps using hybrid wavelet packet-statistical models (WP-SM). In the first step, landslide susceptibility maps were produced using single artificial neural network (ANN), support vector machine (SVM), maximum entropy (MaxEnt), and generalized linear model (GLM). In the next step, the input maps were preprocessed using different mother wavelets in different levels. Then, the hybrid models were developed using the wavelet-based preprocessed maps. Results showed that the wavelet packet transform can be effectively used to produce precise landslide susceptibility maps. It was shown that wavelet packet transform significantly enhanced the ability of the single statistical models. The kappa coefficients were increased from 0.829 to 0.941, 0.846 to 0.978, 0.744 to 0.829, and 0.735 to 0.817 in hybrid ANN, SVM, MaxEnt, and GLM, respectively. The best wavelet transform was performed using bior1.5 with a three-level decomposition. It was also recognized that MaxEnt and GLM produced approximately poor results. However, SVM performed better than the other three models both in single and hybrid forms. ANN also outperformed MaxEnt and GLM models. Spatial distribution of the susceptible area is consistent with the observed landslide distribution pattern particularly in maps obtained from the hybrid models. The produced maps showed that the general pattern of susceptible area intensively followed the pattern of roads and sensitive geological formations.     

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.     

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.     

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.