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.     

GIS-based landslide susceptibility for Arsin-Yomra (Trabzon,North Turkey) region

The purpose of this study is to assess the susceptibility of landslides around Yomra and Arsin towns near Trabzon, in northeast of Turkey, using a geographical information system (GIS). Landslide inventory of the area was made by detailed field surveys and the analyses of the topographical map. The landslide triggering factors are considered to be slope angle, slope aspect, distance from drainage, distance from roads and the weathered lithological units, which were called as “geotechnical units” in the study. Idrisi and ArcGIS packages manipulated all the collected data. Logistic regression (LR) and weighted linear combination (WLC) statistical methods were used to create a landslide susceptibility map for the study area. The results were assessed within the scope of two different points: (a) effectiveness of the methods used and (b) effectiveness of the environmental casual parameters influencing the landslides. The results showed that the WLC model is more suitable than the LR model. Regarding the casual parameters, geotechnical units and slopes were found to be the most important variables for estimating the landslide susceptibility in the study area.     

Application of logistic regression model for slope instability prediction in Cuyahoga River Watershed,Ohio, USA

High incidences of slope movement are observed throughout Cuyahoga River watershed in northeast Ohio, USA. The major type of slope failure involves rotational movement in steep stream walls where erosion of the banks creates over-steepened slopes. The occurrence of landslides in the area depends on a complex interaction of natural as well as human induced factors, including: rock and soil strength, slope geometry, permeability, precipitation, presence of old landslides, proximity to streams and flood-prone areas, land use patterns, excavation of lower slopes and/or increasing the load on upper slopes, alteration of surface and subsurface drainage. These factors were used to evaluate the landslide-induced hazard in Cuyahoga River watershed using logistic regression analysis, and a landslide susceptibility map was produced in ArcGIS. The map classified land into four categories of landslide susceptibility: low, moderate, high, and very high. The susceptibility map was validated using known landslide locations within the watershed area. The landslide susceptibility map produced by the logistic regression model can be efficiently used to monitor potential landslide-related problems, and, in turn, can help to reduce hazards associated with landslides.     

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.     

Assessment of geo-environmental problems of the Zonguldak province (NW Turkey)

The Zonguldak province is a coastal settlement area that has been suffering from serious natural and human-induced environmental problems sourced by its geology and geomorphology. Since the province locates at the heart of a coal-producing basin, the geo-environmental problems related to mining activities such as esthetic degradation, disposal of mining wastes and subsidence of the abandoned coal galleries are badly affecting every day life in Zonguldak province. Disposal of municipal wastes is also a big problem since only one municipality out of 32 has a sanitary disposal area. The rest of the municipalities dispose their solid wastes to rivers or to the sea. The province has also some health problems, which are pointed out in the literature, related to coal mining and geologic environment. These are cytogenetic damage in peripheral lymphocytes and pheumoconiosis (most commonly seen at coal workers), goiter and cancer. Landslides are the most important hazards in the area since 13% of the total surface of the Zonguldak is affected by landslides. In this study, considering the hazard potential special attention is given to deep landslides and using the stepwise forward conditional logistic regression technique, the landslide susceptibility map for the Zonguldak province is produced. The results showed that the most important independent variables governing the landslides are slope gradient, volcanic, and sedimentary rocks of Eocene and clastic and carbonate units of Cretaceous. The landslide map is used as a base map for the production of geo-hazard reconnaissance map on which areas subjected to other important geo-hazards (flood, earthquake and subsidence) are also shown to provide guidance for both existing settlement areas to take the necessary preventive measures and for new developing settlement areas to avoid the problematic areas.     

Landslide susceptibility assessment in the İzmir (West Anatolia,Turkey) city center and its near vicinity by the logistic regression method

A landslide susceptibility assessment for İzmir city (Western Turkey), which is the third biggest city of Turkey, was performed by a logistic regression method. A database of landslide characteristics was prepared using detailed field surveys. The major landslides in the study area are generally observed in the field, dominated by weathered volcanics, and 39.63% of the total landslide area is in this unit. The parameters of lithology, slope gradient, slope aspect, distance to drainage, distance to roads and distance to fault lines were used as variables in the logistic regression analysis. The effect of each parameter on landslide occurrence was assessed from the corresponding coefficients that appear in the logistic regression function. On the basis of the obtained coefficients, lithology plays the most important role in determining landslide occurrence and distribution. Slope gradient has a more significant effect than the other geomorphological parameters, such as slope aspect and distance to drainage. Using a predicted map of probability, the study area was classified into five categories of landslide susceptibility: very low, low, moderate, high and very high. Whereas 49.65% of the total study area has very low susceptibility, very high susceptibility zones make up 11.69% of the area.     

A comparison of landslide susceptibility maps produced by logistic regression, multi-criteria decision, and likelihood ratio methods: a case study at ?zmir, Turkey

The main purpose of this study is to compare the use of logistic regression, multi-criteria decision analysis, and a likelihood ratio model to map landslide susceptibility in and around the city of İzmir in western Turkey. Parameters, such as lithology, slope gradient, slope aspect, faults, drainage lines, and roads, were considered. Landslide susceptibility maps were produced using each of the three methods and then compared and validated. Before the modeling and validation, the observed landslides were separated into two groups. The first group was for training, and the other group was for validation steps. The accuracy of models was measured by fitting them to a validation set of observed landslides. For validation process, the area under curvature (AUC) approach was applied. According to the AUC values of 0.810, 0.764, and 0.710 for logistic regression, likelihood ratio, and multi-criteria decision analysis, respectively, logistic regression was determined to be the most accurate method among the other used landslide susceptibility mapping methods. Based on these results, logistic regression and likelihood ratio models can be used to mitigate hazards related to landslides and to aid in land-use planning.     

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

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

Landslide susceptibility mapping of the slopes in the residual soils of the Mengen region (Turkey) by deterministic stability analyses and image processing techniques

The aim of the present study is to prepare a landslide susceptibility map of a region of about 120 km², between Gökcesu and Pazarköy (around Mengen, NW Turkey) at approximately 10 km north of the North Anatolian Fault Zone, where frequent landslides occur. For this purpose, mechanisms of the landslides were studied by two-dimensional stability analyses together with field observations, and the parameters controlling the development of such slides were identified. Field observations indicated that the failures generally developed within the unconsolidated and/or semiconsolidated soil units in forms of rotational, successive shallow landslides within the weathered zone in Mengen, Cukurca and Sazlar formations. Although consisting of residual soils, Capak and Gökdag formations do not exhibit landslides as the natural slopes formed on these, do not exceed the critical slope angles. Statistical evaluations and distribution of the landslides on the topographical map showed that such parameters as cohesion, angle of internal friction, slope, relative height, orientation of slopes, proximity to drainage pattern, vegetation cover and proximity to major faults were the common features on the landslides. Digital images were obtained to represent all these parameters on gray scale on the SPOT image and on the digital elevation model (DEM) of the area using image processing techniques. Soil mechanics tests were carried out on 36 representative samples collected from different units, and parameters were determined for two-dimensional stability analyses basing on “sensitivity approach” and for the preparation of digital shear strength map. In order to determine the critical slope angle values for the residual soils, a series of sensitivity analyses were realized by using two-dimensional deterministic slope stability analyses techniques for varying values of cohesion, angle of internal friction and slope height along with varying saturation conditions. According to the results of the sensitivity analyses, the Mengen formation was found to be most susceptible unit to landslides, covering about 33.5% of the region studied in terms of surface area. The distribution of the critical slopes were determined by superimposing the critical slope values from sensitivity analyses on slope map of the study area. On the other hand, iso-cohesion and iso-friction maps were produced by locating the values of cohesion and internal friction angles in a geographic coordinate system such that they coincide with sample locations on the DEM and by further interpolation of the values concerned. The pixel values were evaluated in gray scale from 0 to 255, 0 representing the lowest pixel value and 255 representing the highest. Sensitivity analyses on cohesion and angle of internal friction to investigate the effects of these parameters only on stability, revealed that cohesion was effective at a rate of 70% by itself, while angle of internal friction alone controlled the stability by a rate of 30%. The iso-cohesion and iso-friction maps previously obtained were digitally combined in these rates and a “shear strength map” was prepared. The geographic setting of the study area is such that northern slopes usually receive dense precipitation. In relation to this fact, about 42% of the landslides are due north. Thus, a slope orientation map was prepared using the DEM, and slopes facing north were evaluated as being more susceptible to sliding. Proximity to the drainage pattern was another important factor in the evaluation, as streams could adversely affect the stability by either eroding the toe or saturating the slope, or both. When considered together, in conjunction with the field observations, faults and landslides showed a close association. In the area, about 88% of the landslides were detected within an area closer than 250 m to major faults, therefore, a main discontinuity map was produced using the SPOT image of the region, and “proximity to major faults” was evaluated as a parameter as most of the landslides developed in areas where the vegetation was rather sparse. A vegetation cover map was therefore obtained from the SPOT image, and the areas with denser vegetation were considered to be less susceptible to sliding with respect to the areas with less or no vegetation. Having prepared the maps accounting for the distribution of critical slopes, shear strength properties, relative height, slope angle, orientation of the slopes, vegetation cover, proximity to the drainage pattern, geographic corrections were carried on each of these, and a potential failure map was obtained for the residual soils by superimposing all these maps. Next, a classification was performed on the final map and five relative zones of susceptibility were defined. When compared with this map, all of the landslides identified in the field were found to be located in the most susceptible zone. The performance of the method used in processing the images appears to be quite high, the zones determined on the map being the zones of relative susceptibility.     

Rainfall-triggered large landslides on 15 December 2005 in Van Canh District, Binh Dinh Province, Vietnam

Landslides are one of the most dangerous hazards in Vietnam. Most landslides occur at excavated slopes, and natural slope failures are rare in the country. However, the volume of natural slope failures can be very significant and can badly affect large areas. After a long period of heavy rainfall in the fourth quarter of 2005 in Van Canh district, a series of landslides with volumes of 20,000–195,000 m³ occurred on 15 December 2005. The travel distances for the landslides reached over 300–400 m, and the landslides caused some remarkable loud booming noises. The failures took place on natural slopes with unfavorable geological settings and slope angles of 28–31°. The rainfall in the fourth quarter of 2005 is estimated to have a return period of 100 years and was the main triggering factor. Because of the large affected area and low population density, resettling people from the dangerous landslide-prone residential areas to safer sites was the most appropriate solution. In order to do so, a map of landslide susceptibility was produced that took into account slope angle, distance to faults, and slope aspect. The map includes four levels from low to very high susceptibility to landslides.     

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.     

遗传算法优化BP网络在滑坡灾害预测中的应用研究

在陕西省宝鸡市附近长寿沟地区滑坡详细调查和遥感解译的基础上,完成了1∶10000滑坡编目图。通过使用GIS的水文分析功能,运用正反DEM技术,将长寿沟地区划分为216个自然斜坡单元,其中包括123个滑坡单元和93个未发生滑坡单元,分析滑坡发生与坡高、坡度、坡向、坡形、人类工程活动和水文地质条件影响因子之间的统计规律。利用经遗传算法优化后的BP神经网络对80个滑坡样本和40个未滑坡样本进行训练学习,然后再利用训练好的网络对预测样本进行评价分析。结果表明:43个已滑坡单元中只有3个被误判为无滑坡,正确率为9302%,53个未滑坡单元中有10个被预测为滑坡,正确率为8113%,总体正确率为8646%。通过对被预测为滑坡的10个斜坡单元进行分析,发现这些单元在坡形、坡高等影响因素的组合上已经具备了发生滑坡的条件,虽然目前没有发生滑坡,但作为潜在的滑坡危险区,可以为滑坡灾害预测预报和防灾减灾工作提供参考。     

A GIS method for landslide inventory and susceptibility mapping in the Río El Estado watershed,Pico de Orizaba volcano,México

In volcanic terrains, dormant stratovolcanoes are very common and can trigger landslides and debris flows continually along stream systems, thereby affecting human settlements and economic activities. It is important to assess their potential impact and damage through the use of landslide inventory maps and landslide models. In Mexico, numerous geographic information systems (GIS)-based applications have been used to represent and assess slope stability. However, there is no practical and standardized landslide mapping methodology under a GIS. This work provides an overview of the ongoing research project from the Institute of Geography at the National Autonomous University of Mexico that seeks to conduct a multi-temporal landslide inventory and produce a landslide susceptibility map by using GIS. The Río El Estado watershed on the southwestern flank of Pico de Orizaba volcano, the highest mountain in Mexico, is selected as a study area. The geologic and geomorphologic factors in combination with high seasonal precipitation, high degree of weathering, and steep slopes predispose the study area to landslides. The method encompasses two main levels of analysis to assess landslide susceptibility. First, the project aims to derive a landslide inventory map from a representative sample of landslides using aerial orthophotographs and field work. Next, the landslide susceptibility is modelled by using multiple logistic regression implemented in a GIS platform. The technique and its implementation of each level in a GISs-based technology is presented and discussed.     

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.     

Landslides in Sado Island of Japan: Part II. GIS-based susceptibility mapping with comparisons of results from two methods and verifications

Slope instability research and susceptibility mapping is a fundamental component of hazard management and an important basis for provision of measures aimed at decreasing the risk of living with landslides. On this basis, this paper presents the result of a comprehensive study on slope stability analyses and landslide susceptibility mapping carried out in part of Sado Island of Japan. Various types of landslides occurred in the island throughout history. Little is known about the triggering factors and severity of old landslides, but for many of the recent slope failures, the slope characteristics and stratigraphy are such that ground surfaces retain water perennially and landslides occur when additional moisture is induced during rainfall and snowmelt. A range of methods are available in literature for preparation of landslide susceptibility maps. In this study we used two methods namely, the analytical hierarchy process (AHP) and logistic regression, to produce and later compare two susceptibility maps. AHP is a semi-qualitative method, which involves a matrix-based pair-wise comparison of the contribution of different factors for landsliding. Logistic regression on the other hand promotes a multivariate statistical analysis with an objective to find the best-fitting model that describes the relationship between the presence or absence of landslides (dependent variable) and a set of causal factors (independent parameters). Elevation, lithology and slope gradient were casual factors in this study. The determinations of factor weights by AHP and logistic regression were preceded by the calculation of class weights (landslide densities) based on bivariate statistical analyses (BSA). The differences between the AHP derived susceptibility map and the logistic regression counterpart are relatively minor when broad-based classifications are considered. However, with an increase in the number of susceptibility classes, the logistic regression map gave more details but the one derived by AHP failed to do so. The reason is that the majority of pixels in the AHP map have high values, and an increase in the number of classes gives little change in the spatial distribution of susceptibility zones in the middle. To verify the practicality of the two susceptibility maps, both of them were compared with a landslide activity map containing 18 active landslide zones. The outcome was that the active landslide zones do not completely fit into the very high susceptibility class of both maps for various reasons. But 70% of these landslide zones fall into the high and very high susceptibility zones of the AHP map while this is 63% in the case of logistic regression. This indicates that despite the skewed distribution of susceptibility indices, the AHP map was better to capture the reality on the ground than the logistic regression equivalent.     

Shallow-landslide susceptibility in the Costa Viola mountain ridge (southern Calabria,Italy) with considerations on the role of causal factors

The “Costa Viola” mountain ridge (southern Calabria), in the sector between Bagnara Calabra and Scilla, is particularly exposed to geo-hydrological risk conditions. The study area has repeatedly been affected by slope instability events in the last decades, mainly related to debris slides, rock falls and debris flows. These types of slope movements are among the most destructive and dangerous for people and infrastructures, and are characterized by abrupt onset and extremely rapid movements. Susceptibility evaluations to shallow landslides have been performed by only focusing on source activation. A logistic regression approach has been applied to estimating the presence/absence of sources in terms of probability, on the basis of linear statistical relationships with a set of territorial variables. An inventory map of 181 sources, obtained from interpretation of air photographs taken in 1954–1955, has been used as training set, and another map of 81 sources, extracted from 1990 to 1991 photographs, has been adopted for validation purposes. An initial set of 12 territorial variables (i.e. lithology, land use, soil sand percentage, elevation, slope angle, aspect, across-slope and down-slope curvatures, topographic wetness index, distance to road, distance to fault and index of daily rainfall) has been considered. The adopted regression procedure consists of the following steps: (1) parameterization of the independent variables, (2) sampling, (3) calibration, (4) application and (5) evaluation of the forecasting capability. The “best set” of variables could be identified by iteratively excluding one variable at a time, and comparing the ROC results. Through a sensitivity analysis, the role of the considered factors in predisposing shallow slope failures in the study area has been evaluated. The results obtained for the Costa Viola mountain ridge can be considered acceptable, as 98.1 % of the cells are correctly classified. According to the susceptibility map, the village of Scilla and its surroundings fall in the highest susceptibility class.     

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.     

Slope units-based flow susceptibility model: using validation tests to select controlling factors

A susceptibility map for an area, which is representative in terms of both geologic setting and slope instability phenomena of large sectors of the Sicilian Apennines, was produced using slope units and a multiparametric univariate model. The study area, extending for approximately 90 km², was partitioned into 774 slope units, whose expected landslide occurrence was estimated by averaging seven susceptibility values, determined for the selected controlling factors: lithology, mean slope gradient, stream power index at the foot, mean topographic wetness index and profile curvature, slope unit length, and altitude range. Each of the recognized 490 landslides was represented by its centroid point. On the basis of conditional analysis, the susceptibility function here adopted is the density of landslides, computed for each class. Univariate susceptibility models were prepared for each of the controlling factors, and their predictive performance was estimated by prediction rate curves and effectiveness ratio applied to the susceptibility classes. This procedure allowed us to discriminate between effective and non-effective factors, so that only the former was subsequently combined in a multiparametric model, which was used to produce the final susceptibility map. The validation of this map latter enabled us to verify the reliability and predictive performance of the model. Slope unit altitude range and length, lithology and, subordinately, stream power index at the foot of the slope unit demonstrated to be the main controlling factors of landslides, while mean slope gradient, profile curvature, and topographic wetness index gave unsatisfactory results.     

Susceptibility and risk assessment of earthquake-induced landslides based on Landslide Response Units in the Subao River basin,China

Generally, pixels are the basic unit for assessment of landslide susceptibility. However, even if the results facilitate the comparison, a pixel-based analysis does not clearly illustrate the distribution relationships. To eliminate this deficiency, the concept of the Landslide Response Unit (LRU) is proposed in this study, for which adjacent pixels that have similar properties are combined as a basic unit for susceptibility assessment. The Subao River basin, seriously impacted by the Wenchuan Earthquake, was selected as the study area, and three factors including slope gradient, slope aspect, and slope shape, which have a significant impact on landslides, were chosen to divide the basin into 25,984 LRUs. Then topographic, geologic, and distance factors were applied for the landslide susceptibility evaluation. The logistic regression method was used to establish the susceptibility assessing model by analyzing 2,000 susceptible LRUs and 2,000 un-susceptible LRUs. The model accuracy was defined in terms of the ROC curve value and the κ value, 0.531 and 0.84, respectively. The susceptibility of landslides was divided into low, moderate, high, and very high in Subao River basin, and 73% of historical landslides and all four new landslides are in the highly susceptible zone and very highly susceptible zones. Finally, the LRUs with houses, farmlands, and roads prone to sliding and burial hazard were assessed separately. On the basis of considering the potential movement directions of the LRUs, the result found that 1,001 and 835 LRUs probably would be destroyed by slope sliding and landslide burial, respectively.     

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.