GIS-based landslide susceptibility mapping using bivariate statistical analysis in Devrek (Zonguldak-Turkey)

Devrek town with increasing population is located in a hillslope area where some landslides exist. Therefore, landslide susceptibility map of the area is required. The purpose of this study was to generate a landslide susceptibility map using a bivariate statistical index and evaluate and compare the results of the statistical analysis conducted with three different approaches in seed cell concept resulting in different data sets in Geographical Information Systems (GIS) based landslide susceptibility mapping applied to the Devrek region. The data sets are created from the seed cells of (a) crowns and flanks, (b) only crowns, and (c) only flanks of the landslides by using ten different causative parameters of the study area. To increase the data dependency of the analysis, all parameter maps are classified into equal frequency classes based directly on the percentile divisions of each corresponding seed cell data set. The resultant maps of the landslide susceptibility analysis indicate that all data sets produce fairly acceptable results. In each data set analysis, elevation, lithology, slope, aspect, and drainage density parameters are found to be the most contributing factors in landslide occurrences. The results of the three data sets are compared using Seed Cell Area Indexes (SCAI). This comparison shows that the crown data set produces the most accurate and successful landslide susceptibility map of the study area.     

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.     

Application and verification of a fractal approach to landslide susceptibility mapping

Landslide susceptibility mapping is essential for land-use activities and management decision making in hilly or mountainous regions. The existing approaches to landslide susceptibility zoning and mapping require many different types of data. In this study, we propose a fractal method to map landslide susceptibility using historical landslide inventories only. The spatial distribution of landslides is generally not uniform, but instead clustered at many different scales. In the method, we measure the degree of spatial clustering of existing landslides in a region using a box-counting method and apply the derived fractal clustering relation to produce a landslide susceptibility map by means of GIS-supported spatial analysis. The method is illustrated by two examples at different regional scales using the landslides inventory data from Zhejiang Province, China, where the landslides are mainly triggered by rainfall. In the illustrative examples, the landslides from the inventory are divided into two time periods: The landslides in the first period are used to produce a landslide susceptibility map, and those in the late period are taken as validation samples for examining the predictive capability of the landslide susceptibility maps. These examples demonstrate that the landslide susceptibility map created by the proposed technique is reliable.     

Use of fuzzy relations to produce landslide susceptibility map of a landslide prone area (West Black Sea Region, Turkey)

Preparation of landslide susceptibility maps is important for engineering geologists and geomorphologists. However, due to complex nature of landslides, producing a reliable susceptibility map is not easy. For this reason, many procedures have been used to produce such maps. In this study, a new attempt is tried to produce landslide susceptibility map of a part of West Black Sea Region of Turkey. To obtain the fuzzy relations for producing the susceptibility map, a landslide inventory database is compiled by both field surveys and airphoto studies. A total of 266 landslides are identified in the study area, and dominant mode of failure is rotational slide while the other mode of failures are soil flow and shallow translational slide. The landslide inventory and the parameter maps are analyzed together using a computer program (FULLSA) developed in this study. The computer program utilizes the fuzzy relations and produces the landslide susceptibility map automatically. According to this map, 9.6% of the study area is classified as very high susceptibility, 10.3% as high susceptibility, 8.9% as moderate susceptibility, 27.5% as low susceptibility and 43.8% as very low susceptibility or nonsusceptible areas. The prediction performance of the susceptibility map is checked by considering actual landslides in the study area. For this purpose, strength of the relation (r_ij) and the root mean square error (RMSE) values are calculated as 0.867 and 0.284, respectively. These values show that the produced landslide susceptibility map in the present study has a sufficient reliability. It is believed that the approach employed in this study mainly prevents the subjectivity sourced from the parameter selection and provides a support to improve the landslide susceptibility mapping studies.     

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.     

Assessment of landslide susceptibility for a landslide-prone area (north of Yenice, NW Turkey) by fuzzy approach

. Regional landslide susceptibility assessments pose complex problems. To solve these problems, numerous approaches, such as statistical analysis, geotechnical engineering approach, geomorphologic approach and fuzzy logic, have been employed. However, all the available methods for regional landslide susceptibility assessments have some uncertainties due to a lack of knowledge and variability. Minimizing these uncertainties provides realistic approaches. Use of the fuzzy logic approach to produce a landslide susceptibility map of a landslide-prone area in NW Turkey is the main purpose of the present study. For this purpose, the study includes five main stages, these being the preparation of a landslide inventory of the study area, the application of factor analysis, the extraction of fuzzy if-then rules, the use of a geographical information system, and the control of the reliability of the resulting landslide susceptibility map. Slope angle, slope aspect, land use, weathering depth, water conditions and topographical elevation were considered as landslide conditioning factors for the study area. A total of 23 if-then rules was extracted from the field data. Employing these rules, fuzzified index maps representing each parameter were obtained. Finally, combining these maps, the landslide susceptibility map of the area was prepared. When compared with the landslide susceptibility map, the landslides identified in the area were found to be located in the very high- and high-susceptibility zones. As far as the performance of the fuzzy approach for processing is concerned, the images appear to be quite satisfactory, the zones determined on the map being zones of relative susceptibility.     

基于分形的滑坡易发程度区划方法

本文发展了一种基于分形统计的滑坡易发程度评价方法,该方法仅使用已有的滑坡数据,首先通过分形统计获得滑坡分布的分形丛集关系,再通过GIS的空间操作与分析生成滑坡易发程度区划图。提出一种对滑坡易发程度区划图的可信度和预测效果进行评价的方法。本文介绍了这些方法及其在浙江地区应用的结果。     

Application of an analytical hierarchical process approach for landslide susceptibility mapping in A Luoi district,Thua Thien Hue Province,Vietnam

A Luoi is a Vietnamese–Laotian border district situated in the western part of Thua Thien Hue province, central Vietnam, where landslides occur frequently and seriously affect local living conditions. This study focuses on the spatial analysis of landslide susceptibility in this 263-km² area. To analyze landslide manifestation in the study area, causative factor maps are derived of slope angle, weathering, land use, geomorphology, fault density, geology, drainage distance, elevation, and precipitation. The analytical hierarchical process approach is used to combine these maps for landslide susceptibility mapping. A landslide susceptibility zonation map with four landslide susceptibility classes, i.e. low, moderate, high, and very high susceptibility for landsliding, is derived based on the correspondence with an inventory of observed landslides. The final map indicates that about 37% of the area is very highly susceptible for landsliding and about 22% is highly susceptible, which means that more than half of the area should be considered prone to landsliding.     

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.     

Influence of seismic acceleration on landslide susceptibility maps: a case study from NE Turkey (the Kelkit Valley)

Particularly in the last decade, landslide susceptibility and hazard maps have been used for urban planning and site selection of infrastructures. Most of the procedures for preparing of landslide susceptibility maps need high-quality landslide inventory map. Although the rainfall and seismic activities are accepted as triggering factor for landslides, designation of the triggering factor for each landslide in the inventory is almost impossible when well-documented records are unavailable. Therefore, during preparation of landslide susceptibility map, whole landslide records in the inventory map are used together without classifying based on the triggering factors. Although seismic activity is accepted as a triggering factor, possible effect of the use of seismic activity on production of landslide susceptibility map was investigated in this study, and the subject is open to discussion. For this purpose, a series of stability analyses based on circular failure and infinite slope model were performed considering different pseudostatic conditions. The results of analyses show that gentle slopes have higher susceptibility to failure than steeper ones, even if their stability conditions (susceptibilities) are similar for static condition. The seismic forces acting on failure surfaces may not be sufficiently taken into consideration in the conventionally prepared landslide susceptibility maps. Employing the general decreasing trend in stability condition based on slope face angle and the seismic acceleration, a new procedure was introduced for preparing of the landslide susceptibility map for a scenario earthquake. The prediction performance of occurring landslides increased after the procedure was applied to the conventionally prepared landslide susceptibility map. According to the threshold independent spatial performance analyses of the proposed methodology and the produced landslide susceptibility maps, the area under ROC curve values were calculated as 0.801, 0.933, and 0.947 for the maps prepared by considering conventional method and scenario earthquakes having M _w values of 5.5 and 7.5, respectively.     

Evaluation and comparison of GIS based landslide susceptibility mapping procedures in Kulekhani watershed, Nepal

The goal of this paper is to evaluate and compare the consistency of GIS-based heuristic and bivariate landslide susceptibility mapping techniques in the Himalayan region, taking the Kulekhani watershed of central Nepal as an example. For this purpose, a heuristic and two statistical bivariate landslide susceptibility mapping methods are applied, and three separate landslide susceptibility zonation maps are produced. The maps are compared using three approaches: landslide density analysis, success rate analysis, and agreed area analysis. A comparison of the values obtained from landslide density analysis and the curves of success rate analysis indicate that the two bivariate methods produce almost identical results, whereas the map produced with the heuristic method differs significantly from the others. On the other hand, the agreed area analysis highlights significant spatial differences in the maps obtained from the three methods. Although the three approaches evaluate the consistency of susceptibility maps, only the agreed area analysis is capable of spatially comparing them. Hence, this approach proves to be more suitable for spatially and quantitatively evaluating the consistency of various landslide susceptibility zonation maps.     

Comparison of the implementation of rock engineering system and analytic hierarchy process methods,upon landslide susceptibility mapping,using GIS: a case study from the Eastern Achaia County of Peloponnesus,Greece

As landslides are very common in Greece, causing serious problems to the social and economic welfare of many communities, the implementation of a proper hazard analysis system will help the creation of a reliable susceptibility map. Τhis will help local communities to define a safe land use and urban development. The purpose of this study is to compare the implementation of two semi-quantitative landslide assessment approaches, using landslide susceptibility maps compiled in a GIS environment. The compared methods are rock engineering system (RES) and the analytic hierarchy process (AHP). For the landslide susceptibility analysis, the Northeastern part of the Achaia County was examined. This area suffers from many landslides, because of its neighborhood with the tectonically active Corinthian Gulf and its geological setting (Neogene sediments, flysch and other bedrock formations, with local overthrusts). Ten parameters were used in both methodologies, and each one was separated into five categories ranging from 0 to 4, representing their specific conditions derived from the investigation of the landslides in the western part of the study area (ranking area). A layer map was generated for each parameter, using GIS, while the weighting coefficients of each methodology were used for the compilation of RES and AHP final maps of the eastern part of the study area (validating area). By examining these two maps, it is revealed that even though both correctly show the landslide status of the second site, the RES map reveals a better behavior in the spatial distribution of the various landslide susceptibility zones.     

Landslide susceptibility analysis in the Hoa Binh province of Vietnam using statistical index and logistic regression

The purpose of this study is to evaluate and compare the results of applying the statistical index and the logistic regression methods for estimating landslide susceptibility in the Hoa Binh province of Vietnam. In order to do this, first, a landslide inventory map was constructed mainly based on investigated landslide locations from three projects conducted over the last 10 years. In addition, some recent landslide locations were identified from SPOT satellite images, fieldwork, and literature. Secondly, ten influencing factors for landslide occurrence were utilized. The slope gradient map, the slope curvature map, and the slope aspect map were derived from a digital elevation model (DEM) with resolution 20 × 20 m. The DEM was generated from topographic maps at a scale of 1:25,000. The lithology map and the distance to faults map were extracted from Geological and Mineral Resources maps. The soil type and the land use maps were extracted from National Pedology maps and National Land Use Status maps, respectively. Distance to rivers and distance to roads were computed based on river and road networks from topographic maps. In addition, a rainfall map was included in the models. Actual landslide locations were used to verify and to compare the results of landslide susceptibility maps. The accuracy of the results was evaluated by ROC analysis. The area under the curve (AUC) for the statistical index model was 0.946 and for the logistic regression model, 0.950, indicating an almost equal predicting capacity.     

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.     

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.     

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.     

Different ways of landslide geometry interpretation in a process of statistical landslide susceptibility and hazard assessment: Horná Súča (western Slovakia) case study

The article draws a comparison between different ways of landslide geometry interpretation in the scope of the statistical landslide hazard and risk assessment processing. The landslides are included as a major input variable, which are compared with all of the input parametric factors. Based on the above comparison the input data are classified and the final map of landslide susceptibility is constructed. Methodology of multivariate conditional analysis has been used for the construction of final maps. Unique condition units was developed by combination of geological map (lithological units) and slope angle map. Lithological units were derived from geological map and subsequently reclassified into 22 classes. Slope angle map was calculated from digital elevation model (contour map at a scale 1:10,000) and reclassified into nine classes. As a case study, a wide area of Horná Súča (western Slovakia) strongly affected by landsliding (predominantly made of Flysch) has been chosen. Spatial data in the form of parametric maps, as well as final statistical data set were processed in GIS GRASS environment. Four different approaches are used for landslides interpretation: (1) area of landslide body including accumulation zone, (2) area of depletion zone, (3) lines of elongated main scarps, (4) lines of main scarp upper edge. For each approach, a zoning map of landslide susceptibility was compiled and these were compared with each other. Depending on the interpretation approach, the final susceptibility zones are markedly different (in tens of percent).     

Macro-zonation of landslide susceptibility in Garamaura-Swarghat-Gambhar section of national highway 21, Bilaspur District, Himachal Pradesh (India)

Landslides the most common geo-hazard in hilly terrain are short lived phenomena but cause extraordinary landscape changes and destruction of life and property. The frequency and intensity of landslides occurrences along NH-21 during the rainy season not only disrupts traffic movement but also misbalance the agro-economic and developmental activities of the region frittering away thousand crores of rupees from the exchequer. An assessment of landslide susceptibility is, therefore, a prerequisite for sustainable development of the region. The present study deals with the preparation of macro-zonation maps of landslide susceptibility in an area of about 100 sq km on 1:50,000 scale across Garamaura-Swarghat section of National Highway-21. The map has been prepared by superimposing the terrain evaluation maps in a particular zone such as lithological map, structural map, slope morphometry map, relative relief map, land use and land cover map and hydrological condition map using landslide susceptibility evaluation factor rating scheme and calculating the total estimated susceptibility as per the guidelines of IS: 14496 (Part-2) 1998). Numerical weightages are assigned to the prime causative factors of slope instability such as lithology, structure, slope morphometery, relative relief, land use and groundwater conditions as per the scheme approved by Bureau of Indian Standard for the purpose of landslide susceptibility zonation. The area depicts zones of different instability. The identified susceptibility zones compared with landslide intensity in the area show some congruence with the weightages of the inputs. The incongruence in intensity and frequency of landslide occurrences and the inferred susceptibility zones of BIS scheme allow other geotechnical considerations and causative factors to be incorporated for the landslide susceptibility zonation.     

Use of Geomorphological Information in Indirect Landslide Susceptibility Assessment

The objective of this paper is to evaluate the importance of geomorphological expert knowledge in the generation of landslide susceptibility maps, using GIS supported indirect bivariate statistical analysis. For a test area in the Alpago region in Italy a dataset was generated at scale 1:5,000. Detailed geomorphological maps were generated, with legends at different levels of complexity. Other factor maps, that were considered relevant for the assessment of landslide susceptibility, were also collected, such as lithology, structural geology, surficial materials, slope classes, land use, distance from streams, roads and houses. The weights of evidence method was used to generate statistically derived weights for all classes of the factor maps. On the basis of these weights, the most relevant maps were selected for the combination into landslide susceptibility maps. Six different combinations of factor maps were evaluated, with varying geomorphological input. Success rates were used to classify the weight maps into three qualitative landslide susceptibility classes. The resulting six maps were compared with a direct susceptibility map, which was made by direct assignment of susceptibility classes in the field. The analysis indicated that the use of detailed geomorphological information in the bivariate statistical analysis raised the overall accuracy of the final susceptibility map considerably. However, even with the use of a detailed geomorphological factor map, the difference with the separately prepared direct susceptibility map is still significant, due to the generalisations that are inherent to the bivariate statistical analysis technique.     

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.