PROMETHEE II and fuzzy AHP: an enhanced GIS-based landslide susceptibility mapping

This paper describes the application of a well-known multi-criteria decision-making technique, called preference ranking organization method for enrichment evaluation (PROMETHEE II), in combination with fuzzy analytical hierarchy process (FAHP), as a weighting technique to explore landslide susceptibility mapping (LSM). To this end, eight landslide-related geodata layers of the Minoo Dasht located in the Gorgan province of Iran, involving slope, aspect, distance to river, drainage density, distance to fault, mean annual rainfall, distance to road and lithology have been integrated using the PROMETHEE II enhanced by FAHP technique. Afterward, the receiver operating characteristics (ROC) curves for the proposed LSM were drawn using an inventory of landslides containing 83 recent and historic landslide points, and the area under curve = 0.752 value was calculated accordingly. Additionally, to further verify the practicality of such susceptibility map, it was also evaluated against the landslide inventory using simple overlay. The outcome was that about 11 % of the occurred landslide points fall into the very high susceptibility class of the LSM, but approximately 52 % of them indeed fall into the high and very high susceptibility zones together. Also, it resulted that no recorded landslide occurred in the zone of very low susceptibility. According to the results of the ROC curves analysis and simple overlay evaluation, the produced map has exhibited good performance.     

Landslide susceptibility mapping using analytical hierarchy process (AHP) in Tehri reservoir rim region,Uttarakhand

A comprehensive use of analytical hierarchy process (AHP) method in landslide susceptibility mapping (LSM) has been presented for rim region of Tehri reservoir. Using remote sensing data, various landslide causative factors responsible for inducing instability in the area were derived. Ancillary data such as geological map, soil map, and topographic map were also considered along with remote sensing data. Exhaustive field checks were performed to define the credibility of the random landslide conditioning factors considered in this study. Apart from universally acceptable inherent causative factors used in the susceptibility mapping, others such as impact of reservoir impoundment on terrain, topographic wetness index and stream power index were found to be important causative factors in rim region of the Tehri reservoir. The AHP method was used to acquire weights of factors and their classes respectively. Weights achieved from AHP method matched with the existing field conditions. Acceptable consistency ratio (CR) value was achieved for each AHP matrix. Weights of each factor were integrated with weighted sum technique and a landslide susceptibility index map was generated. Jenk’s natural break classifier was used to classify LSI map into very low, low, moderate, high and very high landslide susceptible classes. Validation of the susceptibility map was performed using cumulative percentage/success rate curve technique. Area under curve value of the success rate curve was converted to percentage validation accuracy and a reasonable 78.7% validation accuracy was achieved.     

Landslide Susceptibility Mapping Using Fuzzy Logic System and Its Influences on Mainlines in Lashgarak Region,Tehran, Iran

Landslide susceptibility mapping is among the useful tools applied in disaster management and planning development activities in mountainous areas. The susceptibility maps prepared in this research provide valuable information for landslide hazard management in Lashgarak region of Tehran. This study was conducted to, first, prepare landslide susceptibility maps for Lashgarak region and evaluate landslide effect on mainlines and, second, to analyze the main factors affecting landslide hazard increase in the study area in order to propose efficient strategies for landslide hazard mitigation. A GIS-based multi-criteria decision analysis model (fuzzy logic) is used in the present work for scientific evaluation of landslide susceptible areas in Lashgarak region. To this end, ArcGIS, PCIGeomatica, and IDIRISI software packages were used. Eight information layers were selected for information analysis: ground strength class, slope angle, terrain roughness, normalized difference moisture index, normalized difference vegetation index, distance from fault, distance from the river, and distance from the road. Next, eight different scenarios were created to determine landslide susceptibility of the study area using different operators (intersection (AND), union (OR), algebraic sum (SUM), multiplication (PRODUCT), and different fuzzy gamma values) of fuzzy overlay approach. After that, the performance of various fuzzy operators in landslide susceptibility mapping was empirically compared. The results revealed the excellent consistency of landslide susceptibility map prepared using the fuzzy union (OR) operator with landslide distribution map in the study area. Eventually, the accuracy of landslide susceptibility map prepared using the fuzzy union (OR) operator was evaluated using the frequency ratio diagram. The results showed that frequency values of the landslides gradually increase from “low susceptibility” to high “susceptibility” as 88.34% of the landslides are categorized into two “high” and “very high” susceptibility classes, implying the satisfactory consistency between the landslide susceptibility map prepared using fuzzy union (OR) operator and landslide distribution map.     

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.

     

A fuzzy set based approach for integration of thematic maps for landslide susceptibility zonation

Spatial prediction of landslides is termed landslide susceptibility zonation (LSZ). In this study, an objective weighting approach based on fuzzy concepts is used for LSZ in a part of the Darjeeling Himalayas. Relevant thematic layers pertaining to landslide causative factors have been generated using remote sensing and geographic information system (GIS) techniques. The membership values for each category of thematic layers have been determined using the cosine amplitude fuzzy similarity method and are used as ratings. The integration of these ratings led to the generation of LSZ map. The integration of different ratings to generate an LSZ map has been performed using a fuzzy gamma operator apart from the arithmetic overlay approach. The process is based on determination of combined rating known as the landslide susceptibility index (LSI) for all the pixels using the fuzzy gamma operator and classification using the success rate curve method to prepare the LSZ map. The results indicate that as the gamma value increases, the accuracy of the LSZ map also increases. It is observed that the LSZ map produced by the fuzzy algebraic sum has reflected a more real situation in terms of landslides in the study area.     

Landslide susceptibility mapping for a landslide-prone area (Findikli,NE of Turkey) by likelihood-frequency ratio and weighted linear combination models

Landslides are very common natural problems in the Black Sea Region of Turkey due to the steep topography, improper use of land cover and adverse climatic conditions for landslides. In the western part of region, many studies have been carried out especially in the last decade for landslide susceptibility mapping using different evaluation methods such as deterministic approach, landslide distribution, qualitative, statistical and distribution-free analyses. The purpose of this study is to produce landslide susceptibility maps of a landslide-prone area (Findikli district, Rize) located at the eastern part of the Black Sea Region of Turkey by likelihood frequency ratio (LRM) model and weighted linear combination (WLC) model and to compare the results obtained. For this purpose, landslide inventory map of the area were prepared for the years of 1983 and 1995 by detailed field surveys and aerial-photography studies. Slope angle, slope aspect, lithology, distance from drainage lines, distance from roads and the land-cover of the study area are considered as the landslide-conditioning parameters. The differences between the susceptibility maps derived by the LRM and the WLC models are relatively minor when broad-based classifications are taken into account. However, the WLC map showed more details but the other map produced by LRM model produced weak results. The reason for this result is considered to be the fact that the majority of pixels in the LRM map have high values than the WLC-derived susceptibility map. In order to validate the two susceptibility maps, both of them were compared with the landslide inventory map. Although the landslides do not exist in the very high susceptibility class of the both maps, 79% of the landslides fall into the high and very high susceptibility zones of the WLC map while this is 49% for the LRM map. This shows that the WLC model exhibited higher performance than the LRM model.     

Land suitability assessment for sugarcane cultivation in Bijnor district,India using geographic information system and fuzzy analytical hierarchy process

Land suitability evaluation is prerequisite for assessing the limitations for sustainable land use planning. We used ten site specific criteria (rainfall, texture, drainage, soil depth, slope, distance to major road, distance to nearest sugar mill, erosion hazard, risk of flooding and pH) and applied weighted multi-criteria evaluation (MCE) technique in a geographic information system (GIS) environment to evaluate land suitability for sugarcane cultivation in Bijnor district, India. The weightage of all the parameters was calculated through fuzzy analytical hierarchy process. Sugarcane suitability map was prepared integrating various parameters through weighted overlay analysis. The map was categorized as highly suitable (S1), moderately suitable (S2), marginally suitable (S3) and unsuitable (N). The analysis revealed that of the total cultivable land of the district, largest area (61%) was highly suitable followed by moderately suitable (24%), marginally suitable (7%) and unsuitable (8%) for sugarcane cultivation. Nagina, Najibabad and Bijnor sub-districts need attention of land managers and policy makers to remove the limitations and increase the suitability of sugarcane in such areas. Only 7% area was unsuitable for sugarcane cultivation. Slope, soil depth and erosion hazard were the major limiting factors making the land unsuitable for sugarcane cultivation. Therefore, these areas should be given priority for land and soil restoration efforts. The study showed effectiveness of integrated GIS and MCE approach for land suitability analysis of sugarcane.     

基于要素分析和二元统计模型的区域滑坡危险等级制图——以国道212线陇南段为例

国道212线陇南段是我国地质灾害最发育的地区之一,绘制该区的滑坡危险等级地图对灾害管理和发展规划是极其必要的。基于滑坡的野外调查、机理研究和室内试验等工作,分析了滑坡与各种要素的相关性,选择控制滑坡的9个重要要素作为评价要素,利用GIS和二元统计的信息值模型和滑坡先验风险要素模型绘制了研究区的滑坡危险等级地图。最后,选用区内11个具有明显滑动位移的活动滑坡与滑坡危险等级地图比较,检验其可靠度。结果表明,活动的滑坡绝大部分都位于危险等级很高和高的范围内,说明两种模型的评价结果与研究区实际情况相吻合,同时也反映出信息值模型与实际情况更加相符。     

Landslide susceptibility mapping by using fuzzy logic: a case study of Cham-gardalan catchment,Ilam, Iran

Landslides are introduced as regional movements, which influence different engineering structures such as roads, railways, and dams and cause the person’s death. Identification of landslide zones may decrease the financial losses and human injuries or deaths. This study tries to achieve a landslide susceptibility mapping in Cham-gardalan catchment by weighting the main criteria and the membership functions of fuzzy logic. For this, we applied the best relationship function between the presence and absence of landslides as well as a collection of the elements. At first, the landslide points were identified by the means of some components those of satellite images, topographical (1:50,000) and geographical (1:100,000) maps, field visits, and Google Earth software followed by the preparation of landslide distribution maps. Then, all effective landslide factors such as percentage of slope, slope aspect, height, geology, land uses, distance from roads, distance from drainages, distance from breakage, and precipitation map have been utilized in order to conduct the fuzzy analyses. Landslide susceptibility map was performed by fuzzy operators (Gamma, Product, Sum, Or, And) in the study area. After fuzzificating and weighting, the effective criteria of landslides were determined through fuzzy Gamma operators with the landaus of 0.2, 0.5, 0.8, and 0.9 and by comparing final maps for making an appropriate model of landslide susceptibility mapping. The regional susceptibility map represents the landslide-prone areas in five categories those of very low, low, moderate, high, and very high. Our results indicated that among the applied operators, Gamma with landau of 0.9 can be used as an appropriate method for mapping the landslide susceptibility due to the suitable fuzzification of given criteria based on landslide distribution maps. In addition, the elements of road, percentage of slope, distance from drainage, and geology were recognized as the most important factors for occurring the landslides.     

A subjective and objective integrated weighting method for landslides susceptibility mapping based on GIS

The purpose of this study is to present a weighting method, integrating subjective weight with objective weight, for landslides susceptibility mapping based on geographical information system (GIS). First, the landslide inventory, aspect, slope, proximity to streams of drainage network, proximity to railway, proximity to road, topography, elevation, lithology, tectonic activity and annual precipitation, including their subclasses, were taken as independent landslide causal factors. Second, objective weights of the causal factors were calculated according to the landslide area density based on entropy weighting method, and key factors were selected according to the rank of the objective weights. Third, trapezoidal fuzzy number weighting approach was used to assess the sub-classes of each key factor. Finally, a case study was carried out in Guizhou province, China. A landslide susceptibility map was created using weighted linear combination model based on GIS. Using a predicted map of probability, the study area was classified into four categories of landslide susceptibility: low, moderate, moderate-high, and high.     

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.

     

Use of satellite remote sensing data in the mapping of global landslide susceptibility

Satellite remote sensing data has significant potential use in analysis of natural hazards such as landslides. Relying on the recent advances in satellite remote sensing and geographic information system (GIS) techniques, this paper aims to map landslide susceptibility over most of the globe using a GIS-based weighted linear combination method. First, six relevant landslide-controlling factors are derived from geospatial remote sensing data and coded into a GIS system. Next, continuous susceptibility values from low to high are assigned to each of the six factors. Second, a continuous scale of a global landslide susceptibility index is derived using GIS weighted linear combination based on each factor’s relative significance to the process of landslide occurrence (e.g., slope is the most important factor, soil types and soil texture are also primary-level parameters, while elevation, land cover types, and drainage density are secondary in importance). Finally, the continuous index map is further classified into six susceptibility categories. Results show the hot spots of landslide-prone regions include the Pacific Rim, the Himalayas and South Asia, Rocky Mountains, Appalachian Mountains, Alps, and parts of the Middle East and Africa. India, China, Nepal, Japan, the USA, and Peru are shown to have landslide-prone areas. This first-cut global landslide susceptibility map forms a starting point to provide a global view of landslide risks and may be used in conjunction with satellite-based precipitation information to potentially detect areas with significant landslide potential due to heavy rainfall.     

Pathways for adaptive and integrated disaster resilience

The GIS-multicriteria decision analysis (GIS-MCDA) technique is increasingly used for landslide hazard mapping and zonation. It enables the integration of different data layers with different levels of uncertainty. In this study, three different GIS-MCDA methods were applied to landslide susceptibility mapping for the Urmia lake basin in northwest Iran. Nine landslide causal factors were used, whereby parameters were extracted from an associated spatial database. These factors were evaluated, and then, the respective factor weight and class weight were assigned to each of the associated factors. The landslide susceptibility maps were produced based on weighted overly techniques including analytic hierarchy process (AHP), weighted linear combination (WLC) and ordered weighted average (OWA). An existing inventory of known landslides within the case study area was compared with the resulting susceptibility maps. Respectively, Dempster-Shafer Theory was used to carry out uncertainty analysis of GIS-MCDA results. Result of research indicated the AHP performed best in the landslide susceptibility mapping closely followed by the OWA method while the WLC method delivered significantly poorer results. The resulting figures are generally very high for this area, but it could be proved that the choice of method significantly influences the results.     

A comparative evaluation of heuristic and bivariate statistical modelling for landslide susceptibility mappings in Ghurmi–Dhad Khola, east Nepal

Landslides are one of the major natural disasters that occur in the Himalayan range with recurring frequency, causing enormous loss of life and property every year. Preparation of landslide inventory maps and landslide susceptibility zonation maps are the important tasks to be taken into account initially for safe mitigation measures. The present paper focuses on landslide susceptibility maps of the Ghurmi–Dhad Khola area, east Nepal, using Geographic Information System. For this purpose, the landslide susceptibility maps are prepared by using the heuristic and bivariate statistical methods. The parameters considered for the study are slope angle, slope aspect, elevation, distance from drainage, geology, land cover, rock and soil type, and distance from faults and folds. The landslide susceptibility zonation map produced from the heuristic method shows that 42.59 % of the observed landslide falls under the very high susceptible zone and 33.00 % under the high susceptible zone. Likewise, the landslide susceptibility zonation map produced from the bivariate method depicts that 44.19 % of the observed landslide falls under the very high susceptible zone and 31.59 % under the high susceptible zone. Both the landslide susceptibility zonation maps are identical, and success rates of both the maps are above 80 %. While comparing the landslide susceptibility maps obtained from two different methods, about 78 % of the study area falls in the identical susceptible zones. Special attention should be taken into consideration for the construction works in the areas which have been spatially agreed as very high and high susceptible zones from both techniques. Moreover, these maps can be used for slope management, land use planning, disaster management planning, etc., by the concerned authorities.     

Application of a RS- and GIS-based semi-quantitative approach (analytical hierarchy process – AHP) in landslide hazard risk assessment of the Shivkhola Watershed,Darjiling Himalaya

The present study deals with the application of analytical hierarchy process to prepare landslide hazard risk map of the Shivkhola Watershed applying remote sensing and geographic information system (GIS). Firstly, to integrate all the required thematic data layers and to prepare landslide susceptibility map, prioritised class rating value and prioritised factor rating value were obtained by developing couple-comparing matrix with a reasonable consistency and with the help of MATLAB software after Saaty. Three important risk factor/element maps, that is, weighted land use/land cover map, road contributing area map and settlement density map, were developed and their weighted linear combination was performed to prepare landslide risk exposure map. Then by integrating landslide susceptibility map and landslide risk exposure map, a classification was incorporated on ARC GIS Platform to prepare landslide hazard risk map. To evaluate the validity of the landslide hazard risk map, probability/chance of landslide hazard risk event has been estimated by means of frequency ratio between landslide hazard risk area (%) and number of risk events (%) for each landslide hazard risk class. Finally, an accuracy assessment was also made on ERDAS Imagine (8.5) which depicts that the classification accuracy of the landslide hazard risk map was 92.89 with overall Kappa statistics of 0.8929.     

Comparing GIS-based support vector machine kernel functions for landslide susceptibility mapping

This study compares the predictive performance of GIS-based landslide susceptibility mapping (LSM) using four different kernel functions in support vector machines (SVMs). Nine possible causal criteria were considered based on earlier similar studies for an area in the eastern part of the Khuzestan province of southern Iran. Different models and the resulting landslide susceptibility maps were created using information on known landslide events from a landslide inventory dataset. The models were trained using landslide inventory dataset. A two-step accuracy assessment was implemented to validate the results and to compare the capability of each function. The radial basis function was identified as the most efficient kernel function for LSM with the resulting landslide susceptibility map showing the highest predictive accuracy, followed by the polynomial kernel function. According to the obtained results, it concluded that using SVMs can generally be considered to be an effective method for LSM while it demands careful consideration of kernel function. The results of the present research will also assist other researchers to select the best SVM kernel function to use for LSM.     

Landslide susceptibility mapping using GIS and weighted overlay method: a case study from North Jordan

Road instability along the Jerash–Amman highway was assessed using the weighted overlay method in Geographic Information System environment. The landslide susceptibility map was developed from nine contributing parameters. The map of landslide susceptibility was classified into five zones: very low (very stable), low (stable), moderate (moderately stable), high (unstable), and very high (highly unstable). The very high susceptibility and high susceptibility zones covered 15.14% and 31.81% of the study area, respectively. The main factors that made most parts of study area prone to landslides include excessive drainage channels, road cuts, and unfavorable rock strata such as marl and friable sandstone intercalated with clay and highly fractured limestone. Fracture zones are a major player in land instability. The moderate and high susceptibility zones are the most common in urban (e.g., Salhoub and Gaza camp) and agricultural areas. About 34% of the urban areas and 28.82% of the agricultural areas are characterized by the high susceptibility zone. Twenty percent of the Jerash–Amman highway length and 58% of the overall highway length are located in the very high susceptibility zone. The landslide susceptibility map was validated by the recorded landslides. More than 80 of the inventoried landslides are in unstable zones, which indicate that the selected causative factors are relevant and the model performs properly.     

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.     

Landslides susceptibility mapping based on geographical information system,GuiZhou, south-west China

The purpose of this study is to assess the susceptibility of landslides around the area of Guizhou province, in south-west of China, using a geographical information system (GIS). The base map is prepared by visiting the field area and mapping individual landslide at a scale of 1:500,000 topographic maps. In the study, slope, lithology, landslide inventory, tectonic activity, drainage distribution and annual precipitation were taken as independent causal factors. Therefore, six causal factors maps are prepared by collecting information from various authorized sources and converting them in to GIS maps. The susceptibility assessment is based on the qualitative map combination model and trapezoidal fuzzy number weighting (TFNW) approach. Using a predicted map of probability, the study area was classified into four categories of landslide susceptibility: low, moderate, high and very high. In addition, the weighting procedure showed that the TFNW is an efficient method for landslide causal factors weighting.     

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.