Landslide Susceptibility Mapping Using GIS-based Fuzzy Logic and the Analytical Hierarchical Processes Approach: A Case Study in Constantine (North-East Algeria)

The main purpose of this study was to compare and evaluate the performance of two multicriteria models for landslide susceptibility assessment in Constantine, north-east of Algeria. The landslide susceptibility maps were produced using the analytic hierarchy process (AHP) and Fuzzy AHP (FAHP) via twelve landslides conditioning factors, including the slope gradient, lithology, land cover, distance from drainage network, distance from the roads, distance from faults, topographic wetness index, stream power index, slope curvature, Normalized Difference Vegetation Index, slope aspect and elevation. In this study, the mentioned models were used to derive the weighting value of the conditioning factors. For the validation process of these models, the receiver operating characteristic analysis, and the area under the curve (AUC) were applied by comparing the obtained results to The landslide inventory map which prepared using the archives of scientific publications, reports of local authorities, and field survey as well as analyzing satellite imagery. According to the AUC values, the FAHP model had the highest value (0.908) followed by the AHP model (0.777). As a result, the FAHP model is more consistent and accurate than the AHP in this case study. The outcome of this paper may be useful for landslide susceptibility assessment and land use management.

     

Application of fuzzy logic and analytical hierarchy process (AHP) to landslide susceptibility mapping at Haraz watershed, Iran

The main goal of this study is to produce landslide susceptibility maps of a landslide-prone area (Haraz) in Iran by using both fuzzy logic and analytical hierarchy process (AHP) models. At first, landslide locations were identified by aerial photographs and field surveys, and a total of 78 landslides were mapped from various sources. Then, the landslide inventory was randomly split into a training dataset 70?% (55 landslides) for training the models and the remaining 30?% (23 landslides) was used for validation purpose. Twelve data layers, as the landslide conditioning factors, are exploited to detect the most susceptible areas. These factors are slope degree, aspect, plan curvature, altitude, lithology, land use, distance from rivers, distance from roads, distance from faults, stream power index, slope length, and topographic wetness index. Subsequently, landslide susceptibility maps were produced using fuzzy logic and AHP models. For verification, receiver operating characteristics curve and area under the curve approaches were used. The verification results showed that the fuzzy logic model (89.7?%) performed better than AHP (81.1?%) model for the study area. The produced susceptibility maps can be used for general land use planning and hazard mitigation purpose.     

Estimating potential landslide sites of an upland sub-watershed in Western Ghat’s of Kerala (India) through frequency ratio and GIS

The purpose of this study is to assess the susceptibility of landslides in parts of Western Ghats, Kerala, India, using a geographical information system (GIS). Landslide inventory of the area was made by detailed field surveys and the analysis of the topographical maps. The landslide triggering factors are considered to be slope angle, slope aspect, slope curvature, slope length, distance from drainage, distance from lineaments, lithology, land use and geomorphology. ArcGIS version 8.3 was used to manipulate and analyse all the collected data. Probabilistic-likelihood ratio was used to create a landslide susceptibility map for the study area. The result was validated using the Area under Curve (AUC) method and temporal data of landslide occurrences. The validation results showed satisfactory agreement between the susceptibility map and the existing data on landslide locations. As the result, the success rate of the model was (84.46%) and the prediction rate of the model was (82.38%) shows high prediction accuracy. In the reclassified final landslide susceptibility zone map, 5.68% of the total area is classified as critical in nature. The landslide susceptibility map thus produced can be used to reduce hazards associated with landslides and to land cover planning.     

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

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

Landslide 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.     

Frequency ratio model based landslide susceptibility mapping in lower Mae Chaem watershed, Northern Thailand

The purpose of this study is to produce a landslide susceptibility map for the lower Mae Chaem watershed, northern Thailand using a Geographic Information System (GIS) and remotely sensed images. For this purpose, past landslide locations were identified from satellite images and aerial photographs accompanied by the field surveys to create a landslide inventory map. Ten landslide-inducing factors were used in the susceptibility analysis: elevation, slope angle, slope aspect, lithology, distance from lineament, distance from drainage, precipitation, soil texture, land use/land cover (LULC), and NDVI. The first eight factors were prepared from their associated database while LULC and NDVI maps were generated from Landsat-5 TM images. Landslide susceptibility was analyzed and mapped using the frequency ratio (FR) model that determines the level of correlation between locations of past landslides and the chosen factors and describes it in terms of frequency ratio index. Finally, the output map was validated using the area under the curve (AUC) method where the success rate of 80.06% and the prediction rate of 84.82% were achieved. The obtained map can be used to reduce landslide hazard and assist with proper planning of LULC in the future.     

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.     

Landslide susceptibility mapping using rough sets and back-propagation neural networks in the Three Gorges, China

In the Three Gorges of China, there are frequent landslides, and the potential risk of landslides is tremendous. An efficient and accurate method of generating landslide susceptibility maps is very important to mitigate the loss of lives and properties caused by these landslides. This paper presents landslide susceptibility mapping on the Zigui-Badong of the Three Gorges, using rough sets and back-propagation neural networks (BPNNs). Landslide locations were obtained from a landslide inventory map, supported by field surveys. Twenty-two landslide-related factors were extracted from the 1:10,000-scale topographic maps, 1:50,000-scale geological maps, Landsat ETM + satellite images with a spatial resolution of 28.5 m, and HJ-A satellite images with a spatial resolution of 30 m. Twelve key environmental factors were selected as independent variables using the rough set and correlation coefficient analysis, including elevation, slope, profile curvature, catchment aspect, catchment height, distance from drainage, engineering rock group, distance from faults, slope structure, land cover, topographic wetness index, and normalized difference vegetation index. The initial, three-layered, and four-layered BPNN were trained and then used to map landslide susceptibility, respectively. To evaluate the models, the susceptibility maps were validated by comparing with the existing landslide locations according to the area under the curve. The four-layered BPNN outperforms the other two models with the best accuracy of 91.53 %. Approximately 91.37 % of landslides were classified as high and very high landslide-prone areas. The validation results show sufficient agreement between the obtained susceptibility maps and the existing landslide locations.     

Applying different scenarios for landslide spatial modeling using computational intelligence methods

Landslides every year impose extensive damages to human beings in various parts of the world; therefore, identifying prone areas to landslides for preventive measures is essential. The main purpose of this research is applying different scenarios for landslide susceptibility mapping by means of combination of bivariate statistical (frequency ratio) and computational intelligence methods (random forest and support vector machine) in landslide polygon and point formats. For this purpose, in the first step, a total of 294 landslide locations were determined from various sources such as aerial photographs, satellite images, and field surveys. Landslide inventory was randomly split into a testing dataset 70% (206 landslide locations) for training the different scenarios, and the remaining 30% (88 landslides locations) was used for validation purposes. To providing landslide susceptibility maps, 13 conditioning factors including altitude, slope angle, plan curvature, slope aspect, topographic wetness index, lithology, land use/land cover, distance from rivers, drainage density, distance from fault, distance from roads, convergence index, and annual rainfall are used. Tolerance and the variance inflation factor indices were used for considering multi-collinearity of conditioning factors. Results indicated that the smallest tolerance and highest variance inflation factor were 0.31 and 3.20, respectively. Subsequently, spatial relationship between classes of each landslide conditioning factor and landslides was obtained by frequency ratio (FR) model. Also, importance of the mentioned factors was obtained by random forest (RF) as a machine learning technique. The results showed that according to mean decrease accuracy, factors of altitude, aspect, drainage density, and distance from rivers had the greatest effect on the occurrence of landslide in the study area. Finally, the landslide susceptibility maps were produced by ten scenarios according to different ensembles. The receiver operating characteristics, including the area under the curve (AUC), were used to assess the accuracy of the models. Results of validation of scenarios showed that AUC was varying from 0.668 to 0.749. Also, FR and seed cell area index indicators show a high correlation between the susceptibility classes with the landslide pixels and field observations in all scenarios except scenarios 10_RF and 10_SVM. The results of this study can be used for landslides management and mitigation and development activities such as construction of settlements and infrastructure in the future.     

A new approach to use AHP in landslide susceptibility mapping: a case study at Yenice (Karabuk, NW Turkey)

This study aimed to investigate the parameter effects in preparing landslide susceptibility maps with a data-driven approach and to adapt this approach to analytical hierarchy process (AHP). For this purpose, at the first stage, landslide inventory of an area located in the Western Black Sea region of Turkey covering approximately 567?km² was prepared, and a total of 101 landslides were mapped. In order to assess the landslide susceptibility, a total of 13 parameters were considered as the input parameters: slope, aspect, plan curvature, topographical elevation, vegetation cover index, land use, distance to drainage, distance to roads, distance to structural elements, distance to ridges, stream power index, sediment transport capacity index, and wetness index. AHP was selected as the major assessment methodology since the adapted approach and AHP work in data pairs. Adapted to AHP, a similarity relation?Cbased approach, namely landslide relation indicator (LRI) for parameter selection method, was also proposed. AHP and parametric effect analyses were performed by the proposed approach, and seven landslide susceptibility maps were produced. Among these maps, the best performance was gathered from the landslide susceptibility map produced by 9 parameter combinations using area under curve (AUC) approach. For this map, the AUC value was calculated as 0.797, while the others ranged between 0.686 and 0.771. According to this map, 38.3?% of the study area was classified as having very low, 8.5?% as low, 15.0?% as moderate, 20.3?% as high, and 17.9?% as very high landslide susceptibility, respectively. Based on the overall assessments, the proposed approach in this study was concluded as objective and applicable and yielded reasonable results.     

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.     

Spatial distribution analysis and susceptibility mapping of landslides triggered before and after Mw7.8 Gorkha earthquake along Upper Bhote Koshi,Nepal

The 2015 Mw7.8 Gorkha earthquake triggered thousands of landslides of various types scattered over a large area. In the current study, we utilized pre- and post-earthquake high-resolution satellite imagery to compile two landslide inventories before and after earthquake and prepared three landslide susceptibility maps within 404 km² area using frequency ratio (FR) model. From the study, we could map about 519 landslides including 178 pre-earthquake slides and 341 coseismic slides were identified. This study investigated the relationship between landslide occurrence and landslide causative factors, i.e., slope, aspect, altitude, plan curvature, lithology, land use, distance from streams, distance from road, distance from faults, and peak ground acceleration. The analysis showed that the majority of landslides both pre-earthquake and coseismic occurred at slope >30°, preferably in S, SE, and SW directions and within altitude ranging from 1000 to 1500 m and 1500 to 3500 m. Scatter plots between number of landslides per km^?2 (LN) and percentage of landslide area (LA) and causative factors indicate that slope is the most influencing factor followed by lithology and PGA for the landslide formation. Higher landslide susceptibility before earthquake is observed along the road and rivers, whereas landslides after earthquake are triggered at steeper slopes and at higher altitudes. Combined susceptibility map indicates the effect of topography, geology, and land cover in the triggering of landslides in the entire basin. The resultant landslide susceptibility maps are verified through AUC showing success rates of 78, 81, and 77%, respectively. These susceptibility maps are helpful for engineers and planners for future development work in the landslide prone area.     

Relative effect method of landslide susceptibility zonation in weathered granite soil: a case study in Deokjeok-ri Creek,South Korea

The objective of this study was to produce and evaluate a landslide susceptibility map for weathered granite soils in Deokjeok-ri Creek, South Korea. The relative effect (RE) method was used to determine the relationship between landslide causative factors (CFs) and landslide occurrence. To determine the effect of CFs on landslides, data layers of aspect, elevation, slope, internal relief, curvature, distance to drainage, drainage density, stream power index, sediment transport index, topographic wetness index, soil drainage character, soil type, soil depth, forest type, timber age, and geology were analyzed in a geographical information system (GIS) environment. A GIS-based landslide inventory map of 748 landslide locations was prepared using data from previous reports, aerial photographic interpretation, and extensive field work. A RE model was generated from a training set consisting of 673 randomly selected landslides in the inventory map, with the remaining 75 landslides used for validation of the susceptibility map. The results of the analysis were verified using the landslide location data. According to the analysis, the RE model had a success rate of 86.3 % and a predictive accuracy of 88.6 %. The validation results showed satisfactory agreement between the susceptibility map and the existing data on landslide locations. The results of this study can therefore be used to mitigate landslide-induced hazards and to plan land use.     

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-based landslide susceptibility mapping for a problematic segment of the natural gas pipeline,Hendek (Turkey)

A segment of natural gas pipeline was damaged due to landsliding near Hendek. Re-routing of the pipeline is planned, but it requires the preparation of a landslide susceptibility map. In this study, the statistical index (Wi) and weighting factor (Wf) methods have been used with GIS to prepare a landslide susceptibility map of the problematic segment of the pipeline. For this purpose, thematic layers including landslide inventory, lithology, slope, aspect, elevation, land use/land cover, distance to stream, and drainage density were used. In the study area, landslides occur in the unconsolidated to semi-consolidated clayey unit and regolith. The Wf method gives better results than the Wi method. Lithology is found to be the most important aspect in the study area. Based on the findings obtained in this study, the unconsolidated to semi-consolidated clayey unit and alluvium should be avoided during re-routing. Agricultural activities should not be allowed in the close vicinity of the pipeline.     

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.     

Landslide susceptibility evaluation in the Chemoga watershed,upper Blue Nile,Ethiopia

Landslide susceptibility consists of an essential component in the day-to-day activity of human beings. Landslide incidents are typically happening at a low rate of recurrence when compared and in contrast to other events. This might be generated into main natural catastrophes relating to widespread and undesirable sound effects. Landslide hotspot area identification and mapping are used for the regional community to secure from this disaster. Therefore, this research aims to identify the hotspot areas of landslide and to generate maps using GIS, AHP, and multi-criteria decision analysis (MCDA). MCDA techniques are applied under such circumstances to categorize and class decisions for successive comprehensive estimation or else to state possible from impossible potentiality with various landslides. Analytical hierarchy process (AHP) constructively applies for conveying influence to different criteria within multi-criteria decision analysis. The causative landslide identifying factors utilized in this research were elevation, slope, aspect, soil type, lithology, distance to stream, land use/land cover, rainfall, and drainage density achieved from various sources. Subsequently, to explain the significance of each constraint into landslide susceptibility, all factors were found using the AHP technique. Generally, landslide susceptibility map factors were multiplied by their weights to acquire with the AHP technique. The result showed that the AHP methods are comparatively good quality estimators of landslide susceptibility identification in the Chemoga watershed. As the result, the Chemoga watershed landslide susceptibility map classes were classified as 46.52%, 13.83%.18.71%, 15.39%, and 5.55% of the occurred landslide fall to very low, low, moderate, high, and very high susceptibility zones, respectively. Performance and accuracy of modeled maps have been established using GPS field data and Google earth data landslide map and area under curve (AUC) of the receiver operating characteristic curve (ROC). As the result, validation depends on the ROC specifies the accuracy of the map formed with the AHP merged through weighted overly method illustrated very good accuracy of AUC value 81.45%. In general, the research outcomes inveterate the very good test consistency of the generated maps.

     

A comparative study on the landslide susceptibility mapping using evidential belief function and weights of evidence models

The purpose of this study is to produce landslide susceptibility map of a landslide-prone area (Daguan County, China) by evidential belief function (EBF) model and weights of evidence (WoE) model to compare the results obtained. For this purpose, a landslide inventory map was constructed mainly based on earlier reports and aerial photographs, as well as, by carrying out field surveys. A total of 194 landslides were mapped. Then, the landslide inventory was randomly split into a training dataset; 70% (136 landslides) for training the models and the remaining 30% (58 landslides) was used for validation purpose. Then, a total number of 14 conditioning factors, such as slope angle, slope aspect, general curvature, plan curvature, profile curvature, altitude, distance from rivers, distance from roads, distance from faults, lithology, normalized difference vegetation index (NDVI), sediment transport index (STI), stream power index (SPI), and topographic wetness index (TWI) were used in the analysis. Subsequently, landslide susceptibility maps were produced using the EBF and WoE models. Finally, the validation of landslide susceptibility map was accomplished with the area under the curve (AUC) method. The success rate curve showed that the area under the curve for EBF and WoE models were of 80.19% and 80.75% accuracy, respectively. Similarly, the validation result showed that the susceptibility map using EBF model has the prediction accuracy of 80.09%, while for WoE model, it was 79.79%. The results of this study showed that both landslide susceptibility maps obtained were successful and would be useful for regional spatial planning as well as for land cover planning.     

Landslide susceptibility mapping using GIS and digital photogrammetric techniques: a case study from Ardesen (NE-Turkey)

Ardesen is a settlement area which has been significantly damaged by frequent landslides which are caused by severe rainfalls and result in many casualties. In this study a landslide susceptibility map of Ardesen was prepared using the Analytical Hierarchy Process (AHP) with the help of Geographical Information Systems (GIS) and Digital Photogrametry Techniques (DPT). A landslide inventory, lithology–weathering, slope, aspect, land cover, shear strength, distance to the river, stream density and distance to the road thematics data layers were used to create the map. These layer maps are produced using field, laboratory and office studies, and by the use of GIS and DPT. The landslide inventory map is also required to determine the relationship between these maps and landslides using DPT. In the study field in the Hemsindere Formation there are units that have different weathering classes, and this significantly affects the shear strength of the soil. In this study, shear strength values are calculated in great detail with field and laboratory studies and an additional layer is evaluated with the help of the stability studies used to produce the landslide susceptibility map. Finally, an overlay analysis is carried out by evaluating the layers obtained according to their weight, and the landslide susceptibility map is produced. The study area was classified into five classes of relative landslide susceptibility, namely, very low, low, moderate, high, and very high. Based on this analysis, the area and percentage distribution of landslide susceptibility degrees were calculated and it was found that 28% of the region is under the threat of landslides. Furthermore, the landslide susceptibility map and the landslide inventory map were compared to determine whether the models produced are compatible with the real situation resulting in compatibility rate of 84%. The total numbers of dwellings in the study area were determined one by one using aerial photos and it was found that 30% of the houses, with a total occupancy of approximately 2,300 people, have a high or very high risk of being affected by landslides.