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

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

Application of logistic regression and fuzzy operators to landslide susceptibility assessment in Azdavay (Kastamonu,Turkey)

Landslides and their assessments are of great importance since they damage properties, infrastructures, environment, lives and so on. Particularly, landslide inventory, susceptibility, and hazard or risk mapping have become important issues in the last few decades. Such maps provide useful information and can be produced by qualitative or quantitative methods. The work presented in this paper aimed to assess landslide susceptibility in a selected area, covering 570.625 km² in the Western Black Sea region of Turkey, by two quantitative methods. For this purpose, in the first stage, a detailed landslide inventory map was prepared by extensive field studies. A total of 96 landslides were mapped during these studies. To perform landslide susceptibility analyses, six input parameters such as topographical elevation, lithology, land use, slope, aspect and distance to streams were considered. Two quantitative methods, logistic regression and fuzzy approach, were used to assess landslide susceptibility in the selected area. For the fuzzy approach, the fuzzy and, or, algebraic product, algebraic sum and gamma operators were considered. At the final stage, 18 landslide susceptibility maps were produced by the logistic regression and fuzzy operators in a GIS (Geographic Information System) environment. Two performance indicators such as ROC (relative operating characteristics) and cosine amplitude method (r _ij ) were used to validate the final susceptibility maps. Based on the analyses, the landslide susceptibility map produced by the fuzzy gamma operator with a level of 0.975 showed the best performance. In addition, the maps produced by the logistic regression, fuzzy algebraic product and the higher levels of gamma operators showed more satisfactory results, while the fuzzy and, or, algebraic sum maps were not sufficient to provide reliable outputs.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

GIS-based landslide susceptibility evaluation using fuzzy-AHP multi-criteria decision-making techniques in the Abha Watershed,Saudi Arabia

Landslides are natural geological disasters causing massive destructions and loss of lives, as well as severe damage to natural resources, so it is essential to delineate the area that probably will be affected by landslides. Landslide susceptibility mapping (LSM) is making increasing implications for GIS-based spatial analysis in combination with multi-criteria evaluation (MCE) methods. It is considered to be an effective tool to understand natural disasters related to mass movements and carry out an appropriate risk assessment. This study is based on an integrated approach of GIS and statistical modelling including fuzzy analytical hierarchy process (FAHP), weighted linear combination and MCE models. In the modelling process, eleven causative factors include slope aspect, slope, rainfall, geology, geomorphology, distance from lineament, distance from drainage networks, distance from the road, land use/land cover, soil erodibility and vegetation proportion were identified for landslide susceptibility mapping. These factors were identified based on the (1) literature review, (2) the expert knowledge, (3) field observation, (4) geophysical investigation, and (5) multivariate techniques. Initially, analytical hierarchy process linked with the fuzzy set theory is used in pairwise comparisons of LSM criteria for ranking purposes. Thereafter, fuzzy membership functions were carried out to determine the criteria weights used in the development of a landslide susceptibility map. These selected thematic maps were integrated using a weighted linear combination method to create the final landslide susceptibility map. Finally, a validation of the results was carried out using a sensitivity analysis based on receiver operator curves and an overlay method using the landslide inventory map. The study results show that the weighted overlay analysis method using the FAHP and eigenvector method is a reliable technique to map landslide susceptibility areas. The landslide susceptibility areas were classified into five categories, viz. very low susceptibility, low susceptibility, moderate susceptibility, high susceptibility, and very high susceptibility. The very high and high susceptibility zones account for 15.11% area coverage. The results are useful to get an impression of the sustainability of the watershed in terms of landsliding and therefore may help decision makers in future planning and mitigation of landslide impacts.     

Landslide susceptibility mapping and factor effect analysis using frequency ratio in a catchment scale: a case study from Garuwa sub-basin,East Nepal

For assessing landslide susceptibility, the spatial distribution of landslides in the field is essential. The landslide inventory map is prepared on the basis of historical information of individual landslide events from different sources such as previously published reports, satellite imageries, aerial photographs and interview with local inhabitants. Then, the distribution of landslides in the study area is verified with field surveys. However, the selection of contributing factors for modelling landslide susceptibility is an inhibit task. The previous studies show that the factors are chosen as per availability of data. This paper documents the landslide susceptibility mapping in the Garuwa sub-basin, East Nepal using frequency ratio method. Nine different contributing factors are considered: slope aspect, slope angle, slope shape, relative relief, geology, distance from faults, land use, distance from drainage and annual rainfall. To analyse the effect of contributing factors, the landslide susceptibility index maps are generated four times using (a) topographical factors and geological factors, (b) topographical factors, geological factors and land use, (c) topographical factors, geological factors, land use and drainage and (d) all nine causative factors. By comparing with the pre-existing landslides, the fourth case (considering all nine causative factors) yields the best success rate accuracy, i.e. 81.19 %, which is then used to produce the final landslide susceptibility zonation map. Then, the final landslide susceptibility map is validated through chi-square test. The standard chi-square value with 3 degrees of freedom at the 0.001 significance level is 16.3, whereas the calculated chi-square value is 7,125.79. Since the calculated chi-square value is greater than the standard chi-square value, it can be concluded that the landslide susceptibility map is considered as statistically significant. Moreover, the results show that the predicted susceptibility levels are found to be in good agreement with the past landslide occurrences.     

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.     

Use of GIS-based fuzzy logic relations and its cross application to produce landslide susceptibility maps in three test areas in Malaysia

Landslides are one of the most frequent and common natural hazards in Malaysia. Preparation of landslide susceptibility maps is one of the first and most important steps in the landslide hazard mitigation. However, due to complex nature of landslides, producing a reliable susceptibility map is not easy. For this reason, a number of different approaches have been used, including direct and indirect heuristic approaches, deterministic, probabilistic, statistical, and data mining approaches. Moreover, these landslides can be systematically assessed and mapped through a traditional mapping framework using geoinformation technologies. Since the early 1990s, several mathematical models have been developed and applied to landslide hazard mapping using geographic information system (GIS). Among various approaches, fuzzy logic relation for mapping landslide susceptibility is one of the techniques that allows to describe the role of each predisposing factor (landslide-conditioning parameters) and their optimal combination. This paper presents a new attempt at landslide susceptibility mapping using fuzzy logic relations and their cross application of membership values to three study areas in Malaysia using a GIS. The possibility of capturing the judgment and the modeling of conditioning factors are the main advantages of using fuzzy logic. These models are capable to capture the conditioning factors directly affecting the landslides and also the inter-relationship among them. In the first stage of the study, a landslide inventory was complied for each of the three study areas using both field surveys and airphoto studies. Using total 12 topographic and lithological variables, landslide susceptibility models were developed using the fuzzy logic approach. Then the landslide inventory and the parameter maps were analyzed together using the fuzzy relations and the landslide susceptibility maps produced. Finally, the prediction performance of the susceptibility maps was checked by considering field-verified landslide locations in the studied areas. Further, the susceptibility maps were validated using the receiver-operating characteristics (ROC) success rate curves. The ROC curve technique is based on plotting model sensitivity—true positive fraction values calculated for different threshold values versus model specificity—true negative fraction values on a graph. The ROC curves were calculated for the landslide susceptibility maps obtained from the application and cross application of fuzzy logic relations. Qualitatively, the produced landslide susceptibility maps showed greater than 82% landslide susceptibility in all nine cases. The results indicated that, when compared with the landslide susceptibility maps, the landslides identified in the study areas were found to be located in the very high and high susceptibility zones. This shows that as far as the performance of the fuzzy logic relation approach is concerned, the results appeared to be quite satisfactory, the zones determined on the map being zones of relative susceptibility.     

Landslide hazard zonation using analytical hierarchy process along National Highway-3 in mid Himalayas of Himachal Pradesh,India

Landslide zonation studies emphasize on preparation of landslide hazard zonation maps considering major instability factors contributing to occurrence of landslides. This paper deals with geographic information system-based landslide hazard zonation in mid Himalayas of Himachal Pradesh from Mandi to Kullu by considering nine relevant instability factors to develop the hazard zonation map. Analytical hierarchy process was applied to assign relative weightages over all ranges of instability factors of the slopes in study area. To generate landslide hazard zonation map, layers in geographic information system were created corresponding to each instability factor. An inventory of existing major landslides in the study area was prepared and combined with the landslide hazard zonation map for validation purpose. The validation of the model was made using area under curve technique and reveals good agreement between the produced hazard map and previous landslide inventory with prediction accuracy of 79.08%. The landslide hazard zonation map was classified by natural break classifier into very low hazard, low hazard, moderate hazard, high hazard and very high landslide hazard classes in geographic information system depending upon the frequency of occurrence of landslides in each class. The resultant hazard zonation map shows that 14.30% of the area lies in very high hazard zone followed by 15.97% in high hazard zone. The proposed model provides the best-fit classification using hierarchical approach for the causative factors of landslides having complex structure. The developed hazard zonation map is useful for landslide preparedness, land-use planning, and social-economic and sustainable development of the region.     

Landslide susceptibility analysis with a bivariate approach and GIS in Northern Iran

Globally, landslides cause hundreds of billions of dollars in damage and hundreds of thousands of deaths and injuries each year. A landslide susceptibility map describes areas where landslides are likely to occur in the future by correlating some of the principal factors that contribute to landslides with the past distribution of landslides. A case study is conducted in the mountainous northern Iran. In this study, a landslide susceptibility map of the study area was prepared using bivariate method with the help of the geographic information system. Area density (bivariate) method was used to weight landslide-influencing data layers. 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 hazard classes: very low, low, moderate, high, and very high. The percentage distribution of landslide susceptibility degrees was calculated. It was found that about 26% of the study area is classified as very high and high hazard classes.     

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.     

The 2010 Yushu earthquake triggered landslide hazard mapping using GIS and weight of evidence modeling

The Yushu County, Qinghai Province, China, April 14, 2010, earthquake triggered thousands of landslides in a zone between 96°20′32.9″E and 97°10′8.9″E, and 32°52′6.7″N and 33°19′47.9″N. This study examines the use of geographic information system (GIS) technology and Bayesian statistics in creating a suitable landslide hazard-zone map of good predictive power. A total of 2,036 landslides were interpreted from high-resolution aerial photographs and multi-source satellite images pre- and post-earthquake, and verified by selected field checking before a final landslide-inventory map of the study area could be established using GIS software. The 2,036 landslides were randomly partitioned into two subsets: a training dataset, which contains 80 % (1,628 landslides), for training the model; and a testing dataset 20 % (408 landslides). Twelve earthquake triggered landslide associated controlling parameters, such as elevation, slope gradient, slope aspect, slope curvature, topographic position, distance from main surface ruptures, peak ground acceleration, distance from roads, normalized difference vegetation index, distance from drainages, lithology, and distance from all faults were obtained from variety of data sources. Landslide hazard indices were calculated using the weight of evidence model. The landslide hazard map was compared with training data and testing data to obtain the success rate and predictive rate of the model, respectively. The validation results showed satisfactory agreement between the hazard map and the existing landslide distribution data. The success rate is 80.607 %, and the predictive rate is 78.855 %. The resulting landslide hazard map showed five classes of landslide hazard, i.e., very high, high, moderate, low and very low. The landslide hazard evaluation map should be useful for environmental recovery planning and reconstruction work.     

Integration of statistical and heuristic approaches for landslide risk analysis: a case of volcanic mountains in West Java Province,Indonesia

Owing to fragile geo-morphology, extreme climatic conditions, and densely populated settlements and rapid development activities, West Java Province is the most landslide hazardous area in Indonesia. So, a landslide risk map for this province bears a great importance such as for land-use planning. It is however widely accepted that landslide risk analysis is often difficult because of the difficulties involved in landslide hazard assessment and estimation of consequences of future landslide events. For instance, lack of multi-temporal inventory map or records of triggering events is often a major problem in landslide hazard mapping. In this study, we propose a simple technique for converting a landslide susceptibility map into a landslide hazard map, which we have employed for landslide risk analysis in one ideally hazardous part of volcanic mountains in West Java Province. The susceptibility analysis was carried out through correlation between past landslides and eight spatial parameters related to instability, i.e. slope, aspect, relative relief, distance to river, geological units, soil type, land use and distance to road. The obtained susceptibility map was validated using cross-time technique, and was collaborated with the frequency-area statistics to respond to ‘when landslide will occur’ and ‘how large it will be’. As for the judgment of the consequences of future landslides, expert opinion was used considering available literature and characteristic of the study area. We have only considered economic loss in terms of physical damage of buildings, roads and agricultural lands for the landslide risk analysis. From this study, we understand the following: (1) the hazard map obtained from conversion of the susceptibility map gives spatial probability and the area of an expected landslide will be greater than 500m² in the next 2 years, (2) the landslide risk map shows that 24% of the total area is in high risk; 30% in moderate risk; 45% in low risk and no risk covers only 1% of the total area, and (3) the loss will be high in agricultural lands, while it will be low in the road structures and buildings.     

Landslide risk assessment using concepts of danger pixels and fuzzy set theory in Darjeeling Himalayas

Landslide risk assessment (LRA) is a key component of landslide studies. The landslide risk can be defined as the potential for adverse consequences or loss to human population and property due to the occurrence of landslides. The LRA can be regional or site-specific in nature and is an important information for planning various developmental activities in the area. LRA is considered as a function of landslide potential (LP) and resource damage potential (RDP). The LP and RDP are typically characterized by the landslide susceptibility zonation map and the resource map (i.e., land use land cover map) of the area, respectively. Development of approaches for LRA has always been a challenge. In the present study, two approaches for LRA, one based on the concept of danger pixels and the other based on fuzzy set theory, have been developed and implemented to generate LRA maps of Darjeeling Himalayas, India. The LRA map based on the first approach indicates that 1,015 pixels of habitation and 921 pixels of road section are under risk due to landslides. The LRA map derived from fuzzy set theory based approach shows that a part of habitat area (2,496 pixels) is under very high risk due to landslides. Also, another part of habitat area and a portion of road network (7,204 pixels) are under high risk due to landslides. Thus, LRA map based on the concept of danger pixels gives the pixels under different resource categories at risk due to landslides whereas the LRA map based on the concept of fuzzy set theory further refines this result by defining the degree of severity of risk to these categories by putting these into high and low risk zones. Hence, the landslide risk assessment study carried out using two approaches in this paper can be considered in cohesion for assessing the risks due to landslides in a region.     

Inclusion of earthquake strong ground motion in a geographic information system-based landslide susceptibility zonation in Garhwal Himalayas

Garhwal Himalayas are seismically very active and simultaneously suffering from landslide hazards. Landslides are one of the most frequent natural hazards in Himalayas causing damages worth more than one billion US$ and around 200 deaths every year. Thus, it is of paramount importance to identify the landslide causative factors to study them carefully and rank them as per their influence on the occurrence of landslides. The difference image of GIS-derived landslide susceptibility zonation maps prepared for pre- and post-Chamoli earthquake shows the effect of seismic shaking on the occurrence of landslides in the Garhwal Himalaya. An attempt has been made to incorporate seismic shaking parameters in terms of peak ground acceleration with other static landslide causative factors to produce landslide susceptibility zonation map in geographic information system environment. In this paper, probabilistic seismic hazard analysis has been carried out to calculate peak ground acceleration values at different time periods for estimating seismic shaking conditions in the study area. Further, these values are used as one of the causative factors of landslides in the study area and it is observed that it refines the preparation of landslide susceptibility zonation map in seismically active areas like Garhwal Himalayas.