Landslide susceptibility zonation mapping using logistic regression and its validation in Hashtchin Region,northwest of Iran

Landslide susceptibility zonation mapping assists researchers greatly to understand the spatial distribution of slope failure probability in a region. Being extremely useful in reducing landslide hazards, such maps could simply be produced using both qualitative and quantitative methods. In the present study, a multivariate statistical method called ‘logistic regression’ was used to assess landslide susceptibility in Hashtchin region, situated in west of Alborz Mountainsnorthwest of Iran. In this study, two independent variables, categorical (predictor) and continuous, were drawn on together in the model. To identify the region’s landslides use was made of aerial photographs, field studies and topographic maps. To prepare the database of factors affecting the region’s landslides and to determine landslide zones, geographic information system (GIS) was used. Using such information, landslide susceptibility modeling was accomplished. The data related to factors causing landslides were extracted as independent variables in each cell (in 50 m×50 m cells). Then, the whole data were input into the SPSS, Version 18. The prepared database was later analyzed using logistic regression, the forward stepwise method and based on maximum likelihood estimation. Regression equation was determined using obtained constants and coefficients and the landslide susceptibility of the area in grid-cells (pixels) was computed between 0 and 0.9954. The Receiver Operating Characteristic (ROC) curve was used to assess the accuracy of the logistic regression model. The predicting ability of the model was 84.1% given the area under ROC curve. Finally, the degree of success of landslide susceptibility zonation mapping was estimated to be 79%.     

Application of logistic regression for landslide susceptibility zoning of Cekmece Area, Istanbul, Turkey

As a result of industrialization, throughout the world, cities have been growing rapidly for the last century. One typical example of these growing cities is Istanbul, the population of which is over 10 million. Due to rapid urbanization, new areas suitable for settlement and engineering structures are necessary. The Cekmece area located west of the Istanbul metropolitan area is studied, because the landslide activity is extensive in this area. The purpose of this study is to develop a model that can be used to characterize landslide susceptibility in map form using logistic regression analysis of an extensive landslide database. A database of landslide activity was constructed using both aerial-photography and field studies. About 19.2% of the selected study area is covered by deep-seated landslides. The landslides that occur in the area are primarily located in sandstones with interbedded permeable and impermeable layers such as claystone, siltstone and mudstone. About 31.95% of the total landslide area is located at this unit. To apply logistic regression analyses, a data matrix including 37 variables was constructed. The variables used in the forwards stepwise analyses are different measures of slope, aspect, elevation, stream power index (SPI), plan curvature, profile curvature, geology, geomorphology and relative permeability of lithological units. A total of 25 variables were identified as exerting strong influence on landslide occurrence, and included by the logistic regression equation. Wald statistics values indicate that lithology, SPI and slope are more important than the other parameters in the equation. Beta coefficients of the 25 variables included the logistic regression equation provide a model for landslide susceptibility in the Cekmece area. This model is used to generate a landslide susceptibility map that correctly classified 83.8% of the landslide-prone areas.     

A comparison of logistic regression-based models of susceptibility to landslides in western Colorado,USA

The Paonia-McClure Pass area of Colorado has been recognized as a region highly susceptible to mass movement. Because of the dynamic nature of this landscape, accurate methods are needed to predict susceptibility to movement of these slopes. The area was evaluated by coupling a geographic information system (GIS) with logistic regression methods to assess susceptibility to landslides. We mapped 735 shallow landslides in the area. Seventeen factors, as predictor variables of landslides, were mapped from aerial photographs, available public data archives, ETM + satellite data, published literature, and frequent field surveys. A logistic regression model was run using landslides as the dependent factor and landslide-causing factors as independent factors (covariates). Landslide data were sampled from the landslide masses, landslide scarps, center of mass of the landslides, and center of scarp of the landslides, and an equal amount of data were collected from areas void of discernible mass movement. Models of susceptibility to landslides for each sampling technique were developed first. Second, landslides were classified as debris flows, debris slides, rock slides, and soil slides and then models of susceptibility to landslides were created for each type of landslide. The prediction accuracies of each model were compared using the Receiver Operating Characteristic (ROC) curve technique. The model, using samples from landslide scarps, has the highest prediction accuracy (85 %), and the model, using samples from landslide mass centers, has the lowest prediction accuracy (83 %) among the models developed from the four techniques of data sampling. Likewise, the model developed for debris slides has the highest prediction accuracy (92 %), and the model developed for soil slides has the lowest prediction accuracy (83 %) among the four types of landslides. Furthermore, prediction from a model developed by combining the four models of the four types of landslides (86 %) is better than the prediction from a model developed by using all landslides together (85 %).     

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.     

A distributed hydrological–geotechnical model using satellite-derived rainfall estimates for shallow landslide prediction system at a catchment scale

This paper describes the potential applicability of a hydrological–geotechnical modeling system using satellite-based rainfall estimates for a shallow landslide prediction system. The physically based distributed model has been developed by integrating a grid-based distributed kinematic wave rainfall-runoff model with an infinite slope stability approach. The model was forced by the satellite-based near real-time half-hourly CMORPH global rainfall product prepared by NOAA-CPC. The method combines the following two model outputs necessary for identifying where and when shallow landslides may potentially occur in the catchment: (1) the time-invariant spatial distribution of areas susceptible to slope instability map, for which the river catchment is divided into stability classes according to the critical relative soil saturation; this output is designed to portray the effect of quasi-static land surface variables and soil strength properties on slope instability and (2) a produced map linked with spatiotemporally varying hydrologic properties to provide a time-varying estimate of susceptibility to slope movement in response to rainfall. The proposed hydrological model predicts the dynamic of soil saturation in each grid element. The stored water in each grid element is then used for updating the relative soil saturation and analyzing the slope stability. A grid of slope is defined to be unstable when the relative soil saturation becomes higher than the critical level and is the basis for issuing a shallow landslide warning. The method was applied to past landslides in the upper Citarum River catchment (2,310 km²), Indonesia; the resulting time-invariant landslide susceptibility map shows good agreement with the spatial patterns of documented historical landslides (1985–2008). Application of the model to two recent shallow landslides shows that the model can successfully predict the effect of rainfall movement and intensity on the spatiotemporal dynamic of hydrological variables that trigger shallow landslides. Several hours before the landslides, the model predicted unstable conditions in some grids over and near the grids at which the actual shallow landslides occurred. Overall, the results demonstrate the potential applicability of the modeling system for shallow landslide disaster predictions and warnings.     

Landslide hazard mapping using logistic regression model in Mackenzie Valley,Canada

A logistic regression model is developed within the framework of a Geographic Information System (GIS) to map landslide hazards in a mountainous environment. A case study is conducted in the mountainous southern Mackenzie Valley, Northwest Territories, Canada. To determine the factors influencing landslides, data layers of geology, surface materials, land cover, and topography were analyzed by logistic regression analysis, and the results are used for landslide hazard mapping. In this study, bedrock, surface materials, slope, and difference between surface aspect and dip direction of the sedimentary rock were found to be the most important factors affecting landslide occurrence. The influence on landslides by interactions among geologic and geomorphic conditions is also analyzed, and used to develop a logistic regression model for landslide hazard mapping. The comparison of the results from the model including the interaction terms and the model not including the interaction terms indicate that interactions among the variables were found to be significant for predicting future landslide probability and locating high hazard areas. The results from this study demonstrate that the use of a logistic regression model within a GIS framework is useful and suitable for landslide hazard mapping in large mountainous geographic areas such as the southern Mackenzie Valley.     

Predictive landslide susceptibility mapping using spatial information in the Pechabun area of Thailand

For predictive landslide susceptibility mapping, this study applied and verified probability model, the frequency ratio and statistical model, logistic regression at Pechabun, Thailand, using a geographic information system (GIS) and remote sensing. Landslide locations were identified in the study area from interpretation of aerial photographs and field surveys, and maps of the topography, geology and land cover were constructed to spatial database. The factors that influence landslide occurrence, such as slope gradient, slope aspect and curvature of topography and distance from drainage were calculated from the topographic database. Lithology and distance from fault were extracted and calculated from the geology database. Land cover was classified from Landsat TM satellite image. The frequency ratio and logistic regression coefficient were overlaid for landslide susceptibility mapping as each factor’s ratings. Then the landslide susceptibility map was verified and compared using the existing landslide location. As the verification results, the frequency ratio model showed 76.39% and logistic regression model showed 70.42% in prediction accuracy. The method can be used to reduce hazards associated with landslides and to plan land cover.     

Application of logistic regression and fuzzy operators to landslide susceptibility assessment in Azdavay (Kastamonu,Turkey)

Landslides and their assessments are of great importance since they damage properties, infrastructures, environment, lives and so on. Particularly, landslide inventory, susceptibility, and hazard or risk mapping have become important issues in the last few decades. Such maps provide useful information and can be produced by qualitative or quantitative methods. The work presented in this paper aimed to assess landslide susceptibility in a selected area, covering 570.625 km² in the Western Black Sea region of Turkey, by two quantitative methods. For this purpose, in the first stage, a detailed landslide inventory map was prepared by extensive field studies. A total of 96 landslides were mapped during these studies. To perform landslide susceptibility analyses, six input parameters such as topographical elevation, lithology, land use, slope, aspect and distance to streams were considered. Two quantitative methods, logistic regression and fuzzy approach, were used to assess landslide susceptibility in the selected area. For the fuzzy approach, the fuzzy and, or, algebraic product, algebraic sum and gamma operators were considered. At the final stage, 18 landslide susceptibility maps were produced by the logistic regression and fuzzy operators in a GIS (Geographic Information System) environment. Two performance indicators such as ROC (relative operating characteristics) and cosine amplitude method (r _ij ) were used to validate the final susceptibility maps. Based on the analyses, the landslide susceptibility map produced by the fuzzy gamma operator with a level of 0.975 showed the best performance. In addition, the maps produced by the logistic regression, fuzzy algebraic product and the higher levels of gamma operators showed more satisfactory results, while the fuzzy and, or, algebraic sum maps were not sufficient to provide reliable outputs.     

Using multiple logistic regression and GIS technology to predict landslide hazard in northeast Kansas, USA

Landslides in the hilly terrain along the Kansas and Missouri rivers in northeastern Kansas have caused millions of dollars in property damage during the last decade. To address this problem, a statistical method called multiple logistic regression has been used to create a landslide-hazard map for Atchison, Kansas, and surrounding areas. Data included digitized geology, slopes, and landslides, manipulated using ArcView GIS. Logistic regression relates predictor variables to the occurrence or nonoccurrence of landslides within geographic cells and uses the relationship to produce a map showing the probability of future landslides, given local slopes and geologic units. Results indicated that slope is the most important variable for estimating landslide hazard in the study area. Geologic units consisting mostly of shale, siltstone, and sandstone were most susceptible to landslides. Soil type and aspect ratio were considered but excluded from the final analysis because these variables did not significantly add to the predictive power of the logistic regression. Soil types were highly correlated with the geologic units, and no significant relationships existed between landslides and slope aspect.     

Geo-diversity and its hydrological response in relation to landslide susceptibility in the Himalaya: a GIS-based case study

Assessment and inventory of landslide susceptibility are essential for the formulation of successful disaster mitigation plans. The objective of this study was to assess landslide susceptibility in relation to geo-diversity and its hydrological response in the Lesser Himalaya with a case study using Geographic Information System (GIS) technology. The Dabka watershed, which constitutes a part of the Kosi Basin in the Lesser Himalaya, India, in the district of Nainital, has been selected for the case illustration. The study constitutes three GIS modules: geo-diversity informatics, hydro informatics and landslide informatics. Through the integration and superimposing of spatial data and attribute data of all three GIS modules, Landslide Susceptibility Index (LSI) has been prepared to identify the level of susceptibility for landslide hazards. This resonance study, carried out over a period of five years (2007–2011), found that areas of most stressed geo-diversity (comprising very steep slopes above 30°, geology of Lower Krol and Lariakanta formation, geomorphology of moist areas and debris sites, land use of barren land with a very high drainage frequency and spring density) have a high landslide susceptibility because of high rate of average runoff (33 l/s/km²), flood magnitude (307.28 l/s/km²), erosion (398 tons/km²) and landslide density (5–10 landslides/km²). The areas of least stressed geo-diversity (comprising gentle slopes below 10°, geology of Kailakhan and Siwalik formation, geomorphology of depositional terraces, land use of dense forest with low drainage frequency and spring density) have the lowest landslide susceptibility because of the low rate of average runoff (6.27 l/s/km²), flood magnitude (20.49 l/s/km²), erosion (65.80 tons/km²) and landslide density (1–2 landslides/km²).     

Quantitative landslide susceptibility mapping at Pemalang area,Indonesia

For quantitative landslide susceptibility mapping, this study applied and verified a frequency ratio, logistic regression, and artificial neural network models to Pemalang area, Indonesia, using a Geographic Information System (GIS). Landslide locations were identified in the study area from interpretation of aerial photographs, satellite imagery, and field surveys; a spatial database was constructed from topographic and geological maps. The factors that influence landslide occurrence, such as slope gradient, slope aspect, curvature of topography, and distance from stream, were calculated from the topographic database. Lithology was extracted and calculated from geologic database. Using these factors, landslide susceptibility indexes were calculated by frequency ratio, logistic regression, and artificial neural network models. Then the landslide susceptibility maps were verified and compared with known landslide locations. The logistic regression model (accuracy 87.36%) had higher prediction accuracy than the frequency ratio (85.60%) and artificial neural network (81.70%) models. The models can be used to reduce hazards associated with landslides and to land-use planning.     

Implementing landslide path dependency in landslide susceptibility modelling

Landslide susceptibility modelling—a crucial step towards the assessment of landslide hazard and risk—has hitherto not included the local, transient effects of previous landslides on susceptibility. In this contribution, we implement such transient effects, which we term “landslide path dependency”, for the first time. Two landslide path dependency variables are used to characterise transient effects: a variable reflecting how likely it is that an earlier landslide will have a follow-up landslide and a variable reflecting the decay of transient effects over time. These two landslide path dependency variables are considered in addition to a large set of conditioning attributes conventionally used in landslide susceptibility. Three logistic regression models were trained and tested fitted to landslide occurrence data from a multi-temporal landslide inventory: (1) a model with only conventional variables, (2) a model with conventional plus landslide path dependency variables, and (3) a model with only landslide path dependency variables. We compare the model performances, differences in the number, coefficient and significance of the selected variables, and the differences in the resulting susceptibility maps. Although the landslide path dependency variables are highly significant and have impacts on the importance of other variables, the performance of the models and the susceptibility maps do not substantially differ between conventional and conventional plus path dependent models. The path dependent landslide susceptibility model, with only two explanatory variables, has lower model performance, and differently patterned susceptibility map than the two other models. A simple landslide susceptibility model using only DEM-derived variables and landslide path dependency variables performs better than the path dependent landslide susceptibility model, and almost as well as the model with conventional plus landslide path dependency variables—while avoiding the need for hard-to-measure variables such as land use or lithology. Although the predictive power of landslide path dependency variables is lower than those of the most important conventional variables, our findings provide a clear incentive to further explore landslide path dependency effects and their potential role in landslide susceptibility modelling.     

Controlling parameter analyses and hazard mapping for earthquake-triggered landslides: an example from a square region in Beichuan County, Sichuan Province, China

On May 12, 2008, at 1428 hours (Beijing time), a catastrophic earthquake, with a magnitude of Ms 8.0, struck the Sichuan Province, China. About 200,000 landslides, as a secondary geological hazard associated with the earthquake, were triggered over a broad area. These landslides were of almost all types such as shallow, disrupted landslides, rock falls, deep-seated landslides, and rock avalanches. Some of these landslides damaged and destroyed large part of some towns, blocked roads, dammed rivers, and caused other serious damages. The purpose of this study is to detect correlations between landslide occurrence and the surface rupture plane, ground shaking conditions (measured by peak ground acceleration, PGA), lithology, slope gradient, slope aspect, topographic position, and distance from drainages by using two indices, landslide area percentage (LAP) and the landslide number density (LND), based on geographic information system (GIS) technology and statistical analysis method in a square region (study area) of Beichuan County, Sichuan Province, China. There were 5,096 landslides related with the earthquake which were delineated by visual interpretation and selected field checking throughout the study area. The total area (horizontal projection) of the 5,096 landslides is about 41.103 km². The LAP, which is defined as the percentage of the plane area affected by landslides, was 10.276 %, and the LND, means the number of landslides per square kilometers, was 12.74 landslides/km². Statistical analysis results show that both LAP and LND have a positive correlation with slope gradient and a negative correlation with distance from the surface rupture. However, the correlation between the occurrence of landslides with PGA, topographic position, and distance from drainages are uncertain, or has just a little positive correlation. The correlation between landslide and slope aspect also shows the effect of the directivity of the seismic wave. The Zbq formation had the most concentrated landslide activity with the LND value of 21.78 landslides/km _, ² and the ∈₁ q Gr. geological units had the highest LAP value. Furthermore, weight index (W _i) model is performed with a GIS platform to derive landslide hazard index map. The success rate of the model was 71.615 % and, thus, it was valid. In addition, comparison of five landslide controlling parameters’ influence on landslide occurrences was also carried out.     

逻辑回归模型在玉树地震滑坡危险性评价中的应用与检验

2010年4月14日07时49分(北京时间),青海省玉树县发生了Ms7.1级大地震。作者基于高分辨率遥感影像解译与现场调查验证的方法,圈定了2036处本次地震诱发滑坡。这些滑坡受地震地表破裂控制强烈,规模相对较小,常常密集成片分布。滑坡类型多样,以崩塌型滑坡为主,还包括滑动型、流滑型、碎屑流型、复合型等类型的滑坡。本文基于地理信息系统(GIS)与遥感(RS)技术,应用逻辑回归模型开展玉树地震滑坡危险性评价,并对结果合理性进行检验。应用GIS技术建立玉树地震滑坡灾害及相关滑坡影响因子空间数据库,选择高程、斜坡坡度、斜坡坡向、斜坡曲率、与水系距离、坡位、断裂、地层岩性、归一化植被指数(NDVI)、公路、同震地表破裂、地震动峰值加速度(PGA)共12个因子作为玉树地震滑坡影响因子,在GIS平台下将这些因子专题图层栅格化。应用逻辑回归模型得到每个因子分级的回归系数,然后建立滑坡危险性指数分布图。利用玉树地震滑坡空间分布图对滑坡危险性指数图进行检验,正确率达到83.21%。滑坡危险性分级结果表明,在占研究区总面积4.97%的"很高危险度"的较小范围内,实际发育滑坡数量为766个,占总滑坡面积的比例高达37.62%,表明地震滑坡危险性评价结果良好。不同危险性级别的滑坡点密度统计结果表明,滑坡点密度随着危险性级别的升高而非常迅速的升高。     

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.     

A test of transferability for landslides susceptibility models under extreme climatic events: application to the Messina 2009 disaster

A model building strategy is tested to assess the susceptibility for extreme climatic events driven shallow landslides. In fact, extreme climatic inputs such as storms typically are very local phenomena in the Mediterranean areas, so that with the exception of recently stricken areas, the landslide inventories which are required to train any stochastic model are actually unavailable. A solution is here proposed, consisting in training a susceptibility model in a source catchment, which was implemented by applying the binary logistic regression technique, and exporting its predicting function (selected predictors regressed coefficients) in a target catchment to predict its landslide distribution. To test the method, we exploit the disaster that occurred in the Messina area (southern Italy) on 1 October 2009 where, following a 250-mm/8-h storm, approximately two thousand debris flow/debris avalanches landslides in an area of 21 km² triggered, killing 37 people and injuring more than 100, and causing 0.5 M € worth of structural damage. The debris flows and debris avalanches phenomena involved the thin weathered mantle of the Varisican low to high-grade metamorphic rocks that outcrop in the eastern slopes of the Peloritani Mounts. Two 10-km²-wide stream catchments, which are located inside the storm core area, were exploited: susceptibility models trained in the Briga catchment were tested when exported to predict the landslides distribution in the Giampilieri catchment. The prediction performance (based on goodness of fit, prediction skill, accuracy and precision assessment) of the exported model was then compared with that of a model prepared in the Giampilieri catchment exploiting its landslide inventory. The results demonstrate that the landslide scenario observed in the Giampilieri catchment can be predicted with the same high performance without knowing its landslide distribution: we obtained, in fact, a very poor decrease in predictive performance when comparing the exported model to the native random partition-based model.     

Numerical Modeling of the Seasonal Dynamic Characteristics of the Koxkar Glacier,in West Tianshan,China

Landslides commonly occurs in hilly areas and causes an enormous loss iof life and property every year. National highway-1D (NH-1D) is the only road link between the two districts (Kargil and Leh) of Ladakh region that connects these districts with Kashmir valley. The landslide failure record of the recent past along this sector of the highway is not available. The present study documents landslide susceptible zones and records occurrence of 60 landslides during the last 4 years showing an increasing trend in the occurrence of landslides over these years in this sector. The landslide susceptibility zonation map has been prepared based on the numerical rating of ten major factors viz. slope morphometry, lithology, structure, relative relief, land cover, landuse, rainfall, hydrological conditions, landslide incidences and Slope Erosion, categorised the area in different zones of instability based on the intensity of susceptibility. The landslide susceptibility map of the area encompassing 73.03 km² is divided into 150 facets. Out of the total of 150 facets, 85 facets fall in low susceptibility zone covering 43.56 km2 which constitute about 59.65% of the total area under investigation with a record of 5 landslides; 40 facets fall in the moderate susceptibility zone covering 16.94km² which constitutes about 23.19% of the study area with a record of 20 landslides; and 25 facets fall in the high susceptibility zone covering 12.53 km² which constitute about 17.15% of the study area with a record of 35 landslides. Most of the facets which fall in HSZ are attributed to slope modification for road widening.     

Landslide susceptibility mapping using GIS-based multi-criteria decision analysis,support vector machines,and logistic regression

Identification of landslides and production of landslide susceptibility maps are crucial steps that can help planners, local administrations, and decision makers in disaster planning. Accuracy of the landslide susceptibility maps is important for reducing the losses of life and property. Models used for landslide susceptibility mapping require a combination of various factors describing features of the terrain and meteorological conditions. Many algorithms have been developed and applied in the literature to increase the accuracy of landslide susceptibility maps. In recent years, geographic information system-based multi-criteria decision analyses (MCDA) and support vector regression (SVR) have been successfully applied in the production of landslide susceptibility maps. In this study, the MCDA and SVR methods were employed to assess the shallow landslide susceptibility of Trabzon province (NE Turkey) using lithology, slope, land cover, aspect, topographic wetness index, drainage density, slope length, elevation, and distance to road as input data. Performances of the methods were compared with that of widely used logistic regression model using ROC and success rate curves. Results showed that the MCDA and SVR outperformed the conventional logistic regression method in the mapping of shallow landslides. Therefore, multi-criteria decision method and support vector regression were employed to determine potential landslide zones in the study area.     

Scaling and multivariate analysis of medium to large landslide events: Haida Gwaii, British Columbia

Gimbarzevsky (1988) collected an exceptional landsliding inventory for Haida Gwaii, British Columbia that included over 8,000 landsliding vectors covering an area of approximately 10,000 km². This database was never published in the referred literature, despite its regional significance. It was collected prior to widespread application of GIS technologies in landsliding studies, limiting the analyses undertaken at the time. Gimbarzevsky identified landslides using 1:50,000 aerial photographs, and transferred the information to NTS map sheets. In our study, we digitized the landslide vectors from these original map sheets and connected each landslide to a digital elevation model. Lengths of landslide vectors are compared to the landsliding inventory for Haida Gwaii analyzed in Rood (1984), Martin Y et al. BC Can J Earth Sci 39:289–305 (2002); the latter inventory is based on larger-scale aerial photographs (~1:12,000). Rood’s database contains a more complete record of smaller landslides, while the inventory of Gimbarzevsky provides improved statistical representation of less frequent, medium to large landslides. It is suggested that combined landslide delineation at different scales could provide a more complete landslide record. Discriminant analysis was undertaken to assess which of nine predictor variables, chosen on the basis of mechanical theory, best predict failed versus unfailed locations. Seven of the nine variables were found to be statistically significant in discriminating amongst failed and unfailed locations. Results show that 81.7% of original grouped cases were correctly classified.