Landslide susceptibility assessment: what are the effects of mapping unit and mapping method?

Landslide susceptibility assessment forms the basis of any hazard mapping, which is one of the essential parts of quantitative risk mapping. For the same study area, different susceptibility maps can be achieved depending on the type of susceptibility mapping methods, mapping unit, and scale. Although there are various methods of obtaining susceptibility maps, the efficiency and performance of each method should be evaluated. In this study the effect of mapping unit and susceptibility mapping method on landslide susceptibility assessment is investigated. When analyzing the effect of susceptibility mapping method, logistic regression (LR) which is widely used in landslide susceptibility mapping and, spatial regression (SR), which have not been used for landslide susceptibility mapping, are selected. The susceptibility maps with logistic and spatial regression models are obtained using two different mapping units namely slope unit-based and grid-based mapping units. The procedure for investigation of effect of mapping unit on different susceptibility mapping methods is applied to Kumluca watershed, in Bartin Province of Western Black Sea Region, Turkey. 18 factor maps are prepared for landslide susceptibility assessment in the study region. Geographic information systems and remote sensing techniques are used to create the landslide factor maps, to obtain susceptibility maps and to compare the results. The relative operating characteristics (ROC) curve is used to compare the predictive abilities of each model and mapping unit and also the accuracy is evaluated depending on the observations made during field surveys. By analyzing the area under the ROC curve for grid-based and slope unit-based mapping units, it can be concluded that SR model provide better predictive performance (0.774 in grids and 0.898 in slope units) as compared to the LR model (0.744 in grids and 0.820 in slope units). This result is also supported by the accuracy analysis. For both mapping units, the SR model provides more accurate result (0.55 for grids and 0.57 for slope units) than the LR model (0.50 for grids and 0.48 for slopes). The main reason for this better performance is that the spatial correlations between the mapping units are incorporated into the model in SR while this fact is not considered in LR model.     

Landslide susceptibility mapping using certainty factor, index of entropy and logistic regression models in GIS and their comparison at Mugling�CNarayanghat road section in Nepal Himalaya

Landslide susceptibility maps are vital for disaster management and for planning development activities in the mountainous country like Nepal. In the present study, landslide susceptibility assessment of Mugling?CNarayanghat road and its surrounding area is made using bivariate (certainty factor and index of entropy) and multivariate (logistic regression) models. At first, a landslide inventory map was prepared using earlier reports and aerial photographs as well as by carrying out field survey. As a result, 321 landslides were mapped and out of which 241 (75?%) were randomly selected for building landslide susceptibility models, while the remaining 80 (25?%) were used for validating the models. The effectiveness of landslide susceptibility assessment using GIS and statistics is based on appropriate selection of the factors which play a dominant role in slope stability. In this case study, the following landslide conditioning factors were evaluated: slope gradient; slope aspect; altitude; plan curvature; lithology; land use; distance from faults, rivers and roads; topographic wetness index; stream power index; and sediment transport index. These factors were prepared from topographic map, drainage map, road map, and the geological map. Finally, the validation of landslide susceptibility map was carried out using receiver operating characteristic (ROC) curves. The ROC plot estimation results showed that the susceptibility map using index of entropy model with AUC value of 0.9016 has highest prediction accuracy of 90.16?%. Similarly, the susceptibility maps produced using logistic regression model and certainty factor model showed 86.29 and 83.57?% of prediction accuracy, respectively. Furthermore, the ROC plot showed that the success rate of all the three models performed more than 80?% accuracy (i.e. 89.15?% for IOE model, 89.10?% for LR model and 87.21?% for CF model). Hence, it is concluded that all the models employed in this study showed reasonably good accuracy in predicting the landslide susceptibility of Mugling?CNarayanghat road section. These landslide susceptibility maps can be used for preliminary land use planning and hazard mitigation purpose.     

A comparative study of sequential minimal optimization-based support vector machines,vote feature intervals,and logistic regression in landslide susceptibility assessment using GIS

Landslide susceptibility assessment using GIS has been done for part of Uttarakhand region of Himalaya (India) with the objective of comparing the predictive capability of three different machine learning methods, namely sequential minimal optimization-based support vector machines (SMOSVM), vote feature intervals (VFI), and logistic regression (LR) for spatial prediction of landslide occurrence. Out of these three methods, the SMOSVM and VFI are state-of-the-art methods for binary classification problems but have not been applied for landslide prediction, whereas the LR is known as a popular method for landslide susceptibility assessment. In the study, a total of 430 historical landslide polygons and 11 landslide affecting factors such as slope angle, slope aspect, elevation, curvature, lithology, soil, land cover, distance to roads, distance to rivers, distance to lineaments, and rainfall were selected for landslide analysis. For validation and comparison, statistical index-based methods and the receiver operating characteristic curve have been used. Analysis results show that all these models have good performance for landslide spatial prediction but the SMOSVM model has the highest predictive capability, followed by the VFI model, and the LR model, respectively. Thus, SMOSVM is a better model for landslide prediction and can be used for landslide susceptibility mapping of landslide-prone areas.     

GIS-based spatial prediction of debris flows using logistic regression and frequency ratio models for Zêzere River basin and its surrounding area,Northwest Covilhã, Portugal

Landslide susceptibility mapping (LSM) is important for catastrophe management in the mountainous regions. They focus on generating susceptibility maps beginning from landslide inventories and considering the main predisposing parameters. The aim of this study was to assess the susceptibility of the occurrence of debris flows in the Zêzere River basin and its surrounding area using logistic regression (LR) and frequency ratio (FR) models. To achieve this, a landslide inventory map was created using historical information, satellite imagery, and extensive field works. One hundred landslides were mapped, of which 75% were randomly selected as training data, while the remaining 25% were used for validating the models. The landslide influence factors considered for this study were lithology, elevation, slope gradient, slope aspect, plan curvature, profile curvature, normalized difference vegetation index (NDVI), distance to roads, topographic wetness index (TWI), and stream power index (SPI). The relationships between landslide occurrence and these factors were established, and the results were then evaluated and validated. Validation results show that both methods give acceptable results [the area under curve (AUC) of success rates is 83.71 and 76.38 for LR and FR, respectively]. Furthermore, the AUC results for prediction accuracy revealed that LR model has the highest predictive performance (AUC of predicted rate?=?80.26). Hence, it is concluded that the two models showed reasonably good accuracy in predicting the landslide susceptibility in the study area. These two models have the potential to aid planners in development and land-use planning and to offer tools for hazard mitigation measures.     

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.     

基于GIS和层次分析法的沙溪流域滑坡地质灾害易发性评价

滑坡是沙溪流域主要地质灾害类型之一,开展滑坡灾害易发性评价可为区域地质灾害防治提供数据基础和决策依据。通过沙溪流域生态地质调查,分析了滑坡灾害分布规律和影响因素之间的关系,选取岩性建造、地貌、坡度、坡向、降雨量、距河流距离和距断层距离7项指标,利用层次分析法及地理信息系统空间分析技术,开展沙溪流域滑坡地质灾害易发性评价。结果显示: 沙溪流域滑坡易发性影响因子依次为岩性建造、多年年均降水量、地形地貌、坡度、距河流距离、距断层距离和坡向; 沙溪流域滑坡灾害易发性与坡度、岩性建造、年均降水量表现出明显正相关,即坡度越大、岩性建造性质越软弱、越易风化,年均降水量越多,越易引发滑坡灾害; 滑坡灾害易发性与断裂构造、河流距离与滑坡灾害易发性呈负相关,即距离越近越容易诱发地质灾害; 流域整体以低易发区和极低易发区为主,高易发区主要分布在沙溪流域中南部、东部及东北部地区。这为沙溪流域地质灾害防治提供了基础数据和决策依据。     

Combining landslide susceptibility maps obtained from frequency ratio,logistic regression,and artificial neural network models using ASTER images and GIS

Landslide-related factors were extracted from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) images, and integrated techniques were developed, applied, and verified for the analysis of landslide susceptibility in Boun, Korea, using a geographic information system (GIS). Digital elevation model (DEM), lineament, normalized difference vegetation index (NDVI), and land-cover factors were extracted from the ASTER images for analysis. Slope, aspect, and curvature were calculated from a DEM topographic database. Using the constructed spatial database, the relationships between the detected landslide locations and six related factors were identified and quantified using frequency ratio (FR), logistic regression (LR), and artificial neural network (ANN) models. These relationships were used as factor ratings in an overlay analysis to create landslide susceptibility indices and maps. Three landslide susceptibility maps were then combined and applied as new input factors in the FR, LR, and ANN models to make improved susceptibility maps. All of the susceptibility maps were verified by comparison with known landslide locations not used for training the models. The combined landslide susceptibility maps created using three landslide-related input factors showed improved accuracy (87.00% in FR, 88.21% in LR, and 86.51% in ANN models) compared to the individual landslide susceptibility maps (84.34% in FR, 85.40% in LR, and 74.29% in ANN models) generated using the six factors from the ASTER images.     

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.     

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

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

Landslide susceptibility mapping using analytic hierarchy process and information value methods along a highway road section in Constantine,Algeria

This research work deals with the landslide susceptibility assessment using Analytic hierarchy process (AHP) and information value (IV) methods along a highway road section in Constantine region, NE Algeria. The landslide inventory map which has a total of 29 single landslide locations was created based on historical information, aerial photo interpretation, remote sensing images, and extensive field surveys. The different landslide influencing geoenvironmental factors considered for this study are lithology, slope gradient, slope aspect, distance from faults, land use, distance from streams, and geotechnical parameters. A thematic layer map is generated for every geoenvironmental factor using Geographic Information System (GIS); the lithological units and the distance from faults maps were extracted from the geological database of the region. The slope gradient, slope aspect, and distance from streams were calculated from the Digital Elevation Model (DEM). Contemporary land use map was derived from satellite images and field study. Concerning the geotechnical parameters maps, they were determined making use of the geotechnical data from laboratory tests. The analysis of the relationships between the landslide-related factors and the landslide events was then carried out in GIS environment. The AUC plot showed that the susceptibility maps had a success rate of 77 and 66% for IV and AHP models, respectively. For that purpose, the IV model is better in predicting the occurrence of landslides than AHP one. Therefore, the information value method could be used as a landslide susceptibility mapping zonation method along other sections of the A1 highway.     

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

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

基于GIS与WOE-BP模型的滑坡易发性评价

区域滑坡易发性研究对地质灾害风险管理具有重要意义.以往研究中,将多元统计模型与机器学习方法相结合用于滑坡易发性评价的研究较少.以三峡库区万州区为例,首先选取9种指标因子（坡度、坡向、剖面曲率、地表纹理、地层岩性、斜坡结构、地质构造、水系分布及土地利用类型）作为滑坡易发性评价指标.基于证据权模型（weights of evidence,WOE）计算得到的对比度和滑坡面积比与分级面积比的相对大小,对各指标因子进行状态分级;再利用粒子群法优化的BP神经网络模型（PSO-BP）得到各指标因子权重.综合两种模型确定的状态分级权重和指标因子权重（WOE-BP）计算滑坡易发性指数（landslide susceptibility index,LSI）,基于GIS平台得到全区滑坡易发性分区图.结果表明:水系、地层岩性和地质构造是影响万州区滑坡发育的主要指标因子;WOE-BP模型的预测精度为80.8%,优于WOE模型的73.1%和BP神经网络模型的71.6%,可为定量计算指标因子权重和优化滑坡易发性评价提供有效途径.      

A Comparative Study of Least Square Support Vector Machines and Multiclass Alternating Decision Trees for Spatial Prediction of Rainfall-Induced Landslides in a Tropical Cyclones Area

The objective of this study is to explore and compare the least square support vector machine (LSSVM) and multiclass alternating decision tree (MADT) techniques for the spatial prediction of landslides. The Luc Yen district in Yen Bai province (Vietnam) has been selected as a case study. LSSVM and MADT are effective machine learning techniques of classification applied in other fields but not in the field of landslide hazard assessment. For this, Landslide inventory map was first constructed with 95 landslide locations identified from aerial photos and verified from field investigations. These landslide locations were then divided randomly into two parts for training (70 % locations) and validation (30 % locations) processes. Secondly, landslide affecting factors such as slope, aspect, elevation, curvature, lithology, land use, distance to roads, distance to faults, distance to rivers, and rainfall were selected and applied for landslide susceptibility assessment. Subsequently, the LSSVM and MADT models were built to assess the landslide susceptibility in the study area using training dataset. Finally, receiver operating characteristic curve and statistical index-based evaluations techniques were employed to validate the predictive capability of these models. As a result, both the LSSVM and MADT models have high performance for spatial prediction of landslides in the study area. Out of these, the MADT model (AUC = 0.853) outperforms the LSSVM model (AUC = 0.803). From the landslide study of Luc Yen district in Yen Bai province (Vietnam), it can be conclude that the LSSVM and MADT models can be applied in other areas of world also for and spatial prediction. Landslide susceptibility maps obtained from this study may be helpful in planning, decision making for natural hazard management of the areas susceptible to landslide hazards.     

GIS based landslide susceptibility mapping using a fuzzy logic approach: A case study from Ghurmi-Dhad Khola area, Eastern Nepal

Landslides cause extensive loss of life and property in the Nepal Himalaya. Since the late 1980s, different mathematical models have been developed and applied for landslide susceptibility mapping and hazard assessment in Nepal. The main goal of this paper is to apply fuzzy logic to landslide susceptibility mapping in the Ghurmi-Dhad Khola area, Eastern Nepal. Seven causative factors are considered: slope angle, slope aspect, distance from drainage, land use, geology, distance from faults and folds, soil and rock type. Likelihood ratios are obtained for each class of causative factors by comparison with past landslide occurrences. The ratios are normalized between zero and one to obtain fuzzy membership values. Further, different fuzzy operators are applied to generate landslide susceptibility maps. Comparison with the landslide inventory map reveals that the fuzzy gamma operator with a γ-value of 0.60 yields the best prediction accuracy. Consequently, this operator is used to produce the final landslide susceptibility zonation map.     

RS & GIS-based landslide susceptibility mapping of the Balason River basin,Darjeeling Himalaya,using logistic regression (LR) model

The present study deals with the preparation of a landslide susceptibility map of the Balason River basin, Darjeeling Himalaya, using a logistic regression model based on Geographic Information System and Remote Sensing. The landslide inventory map was prepared with a total of 295 landslide locations extracted from various satellite images and intensive field survey. Topographical maps, satellite images, geological, geomorphological, soil, rainfall and seismic data were collected, processed and constructed into a spatial database in a GIS environment. The chosen landslide-conditioning factors were altitude, slope aspect, slope angle, slope curvature, geology, geomorphology, soil, land use/land cover, normalised differential vegetation index, drainage density, lineament number density, distance from lineament, distance to drainage, stream power index, topographic wetted index, rainfall and peak ground acceleration. The produced landslide susceptibility map satisfied the decision rules and ?2 Log likelihood, Cox &; Snell R-Square and Nagelkerke R-Square values proved that all the independent variables were statistically significant. The receiver operating characteristic curve showed that the prediction accuracy of the landslide probability map was 96.10%. The proposed LR method can be used in other hazard/disaster studies and decision-making.     

Comparison of landslide susceptibility mapping methodologies for Koyulhisar,Turkey: conditional probability,logistic regression,artificial neural networks,and support vector machine

This case study presented herein compares the GIS-based landslide susceptibility mapping methods such as conditional probability (CP), logistic regression (LR), artificial neural networks (ANNs) and support vector machine (SVM) applied in Koyulhisar (Sivas, Turkey). Digital elevation model was first constructed using GIS software. Landslide-related factors such as geology, faults, drainage system, topographical elevation, slope angle, slope aspect, topographic wetness index, stream power index, normalized difference vegetation index, distance from settlements and roads were used in the landslide susceptibility analyses. In the last stage of the analyses, landslide susceptibility maps were produced from ANN, CP, LR, SVM models, and they were then compared by means of their validations. However, area under curve values obtained from all four methodologies showed that the map obtained from ANN model looks like more accurate than the other models, accuracies of all models can be evaluated relatively similar. The results also showed that the CP is a simple method in landslide susceptibility mapping and highly compatible with GIS operating features. Susceptibility maps can be easily produced using CP, because input process, calculation and output processes are very simple in CP model when compared with the other methods considered in this study.     

Effects of soil heterogeneity on susceptibility of shallow landslides

Uncertainties in parameters of landslide susceptibility models often hinder them from providing accurate spatial and temporal predictions of landslide occurrences. Substantial contribution to the uncertainties in landslide assessment originates from spatially variable geotechnical and hydrological parameters. These input parameters may often vary significantly through space, even within the same geological deposit, and there is a need to quantify the effects of the uncertainties in these parameters. This study addresses this issue with a new three-dimensional probabilistic landslide susceptibility model. The spatial variability of the model parameters is modeled with the random field approach and coupled with the Monte Carlo method to propagate uncertainties from the model parameters to landslide predictions (i.e., factor of safety). The resulting uncertainties in landslide predictions allow the effects of spatial variability in the input parameters to be quantified. The performance of the proposed model in capturing the effect of spatial variability and predicting landslide occurrence has been compared with a conventional physical-based landslide susceptibility model that does not account for three-dimensional effects on slope stability. The results indicate that the proposed model has better performance in landslide prediction with higher accuracy and precision than the conventional model. The novelty of this study is illustrating the effects of the soil heterogeneity on the susceptibility of shallow landslides, which was made possible by the development of a three-dimensional slope stability model that was coupled with random field model and the Monte Carlo method.

     

A comparison of slope units and grid cells as mapping units for landslide susceptibility assessment

Landslides lead to a great threat to human life and property safety. The delineation of landslide-prone areas achieved by landslide susceptibility assessment plays an important role in landslide management strategy. Selecting an appropriate mapping unit is vital for landslide susceptibility assessment. This paper compares the slope unit and grid cell as mapping unit for landslide susceptibility assessment. Grid cells can be easily obtained and their matrix format is convenient for calculation. A slope unit is considered as the watershed defined by ridge lines and valley lines based on hydrological theory and slope units are more associated with the actual geological environment. Using 70% landslide events as the training data and the remaining landslide events for verification, landslide susceptibility maps based on slope units and grid cells were obtained respectively using a modified information value model. ROC curve was utilized to evaluate the landslide susceptibility maps by calculating the training accuracy and predictive accuracy. The training accuracies of the grid cell-based susceptibility assessment result and slope unit-based susceptibility assessment result were 80.9 and 83.2%, and the prediction accuracies were 80.3 and 82.6%, respectively. Therefore, landslide susceptibility mapping based on slope units performed better than grid cell-based method.     

信息量模型、确定性系数模型与逻辑回归模型组合评价地质灾害敏感性的对比研究

为探索区域地质灾害敏感性评价方法,以宁夏盐池县为研究区域,选取坡度、坡向、坡高、高程、地层、距河流距离、距道路距离、植被覆盖度等8个影响地质灾害发生的评价因子,分别采用信息量模型+逻辑回归模型(I+LR)和确定性系数模型+逻辑回归模型(CF+LR)2种组合模型对盐池县地质灾害敏感性进行评价,将该区域地质灾害划分为极低、低、中和高敏感区4类,并完成结果检验。结果表明：(1)2种组合模型得到的低、中敏感区面积基本相当,而高敏感区面积相差较大,CF+LR模型较I+LR模型高敏感区面积增加约5336 km2,而极低敏感区面积减少约6%;(2)2种组合模型的合理性均符合检验要求,且ROC精度检验AUC值分别为0868和0829,渐进Sigb均小于005,表明2种组合评价模型都能较为客观准确地评价盐池县地质灾害敏感性;(3)ROC检验精度与盐池县地质灾害发育情况均表明I+LR模型精度更高。     

GIS-based landslide susceptibility mapping with probabilistic likelihood ratio and spatial multi-criteria evaluation models (North of Tehran,Iran)

The aim of this study is to produce landslide susceptibility mapping by probabilistic likelihood ratio (PLR) and spatial multi-criteria evaluation (SMCE) models based on geographic information system (GIS) in the north of Tehran metropolitan, Iran. The landslide locations in the study area were identified by interpretation of aerial photographs, satellite images, and field surveys. In order to generate the necessary factors for the SMCE approach, remote sensing and GIS integrated techniques were applied in the study area. Conditioning factors such as slope degree, slope aspect, altitude, plan curvature, profile curvature, surface area ratio, topographic position index, topographic wetness index, stream power index, slope length, lithology, land use, normalized difference vegetation index, distance from faults, distance from rivers, distance from roads, and drainage density are used for landslide susceptibility mapping. Of 528 landslide locations, 70 % were used in landslide susceptibility mapping, and the remaining 30 % were used for validation of the maps. Using the above conditioning factors, landslide susceptibility was calculated using SMCE and PLR models, and the results were plotted in ILWIS-GIS. Finally, the two landslide susceptibility maps were validated using receiver operating characteristic curves and seed cell area index methods. The validation results showed that area under the curve for SMCE and PLR models is 76.16 and 80.98 %, respectively. The results obtained in this study also showed that the probabilistic likelihood ratio model performed slightly better than the spatial multi-criteria evaluation. These landslide susceptibility maps can be used for preliminary land use planning and hazard mitigation purpose.