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.     

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

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

Landslide susceptibility analysis and its verification using likelihood ratio, logistic regression, and artificial neural network models: case study of Youngin, Korea

The likelihood ratio, logistic regression, and artificial neural networks models are applied and verified for analysis of landslide susceptibility in Youngin, Korea, using the geographic information system. From a spatial database containing such data as landslide location, topography, soil, forest, geology, and land use, the 14 landslide-related factors were calculated or extracted. Using these factors, landslide susceptibility indexes were calculated by likelihood ratio, logistic regression, and artificial neural network models. Before the calculation, the study area was divided into two sides (west and east) of equal area, for verification of the models. Thus, the west side was used to assess the landslide susceptibility, and the east side was used to verify the derived susceptibility. The results of the landslide susceptibility analysis were verified using success and prediction rates. The verification results showed satisfactory agreement between the susceptibility map and the existing data on landslide locations.     

Landslide detection and susceptibility mapping using LiDAR and an artificial neural network approach: a case study in the Cuyahoga Valley National Park,Ohio

The purpose of this study was to detect shallow landslides using hillshade maps derived from light detection and ranging (LiDAR)-based digital elevation model (DEM) derivatives. The landslide susceptibility mapping used an artificial neural network (ANN) approach and backpropagation method that was tested in the northern portion of the Cuyahoga Valley National Park (CVNP) located in northeast Ohio. The relationship between landslides and predictor attributes, which describe landform classes using slope, profile and plan curvatures, upslope drainage area, annual solar radiation, and wetness index, was extracted from LiDAR-based DEM using geographic information system (GIS). The approach presented in this paper required a training study area for the development of the susceptibility model and a validation study area to test the model. The results from the validation showed that within the very high susceptibility class, a total of 42.6 % of known landslides that were associated with 1.56 % of total area were correctly predicted. In contrast, the very low susceptibility class that represented 82.68 % of the total area was associated with 1.20 % of known landslides. The results suggest that the majority of the known landslides occur within a small portion of the study area, consistent with field investigation and other studies. Sample probabilistic maps of landslide susceptibility potential and other products from this approach are summarized and presented for visualization to help park officials in effective management and planning.     

Landslide susceptibility mapping using frequency ratio, analytic hierarchy process, logistic regression, and artificial neural network methods at the Inje area, Korea

Every year, the Republic of Korea experiences numerous landslides, resulting in property damage and casualties. This study compared the abilities of frequency ratio (FR), analytic hierarchy process (AHP), logistic regression (LR), and artificial neural network (ANN) models to produce landslide susceptibility index (LSI) maps for use in predicting possible landslide occurrence and limiting damage. The areas under the relative operating characteristic (ROC) curves for the FR, AHP, LR, and ANN LSI maps were 0.794, 0.789, 0.794, and 0.806, respectively. Thus, the LSI maps developed by all the models had similar accuracy. A cross-tabulation analysis of landslide occurrence against non-occurrence areas showed generally similar overall accuracies of 65.27, 64.35, 65.51, and 68.47 % for the FR, AHP, LR, and ANN models, respectively. A correlation analysis between the models demonstrated that the LR and ANN models had the highest correlation (0.829), whereas the FR and AHP models had the lowest correlation (0.619).     

Landslide susceptibility mapping based on frequency ratio and logistic regression models

The aim of this study is to apply and compare a probability model, frequency ratio and statistical model, and a logistic regression to Sajaroud area, Northern Iran using geographic information system. Landslide locations of the study area were detected from interpretation of aerial photographs and field surveys. Landslide-related factors such as elevation, slope gradient, slope aspect, slope curvature, rainfall, distance to fault, distance to drainage, distance to road, land use, and geology were calculated from the topographic and geology map and LANDSAT ETM satellite imagery. The spatial relationships between the landslide location and each landslide-related factor were analyzed and then landslide susceptibility maps were produced using the frequency ratio and forward stepwise logistic regression methods. Finally, the maps were tested and compared using known landslide locations, and success rates were calculated. Predicted accuracy values for frequency ratio (79.48%) and logistic regression models showed that the map obtained from frequency ratio model is more accurate than the logistic regression (77.4%) model. The models used in this study have shown a great deal of importance for watershed management and land use planning.     

Application of frequency ratio,statistical index,and weights-of-evidence models and their comparison in landslide susceptibility mapping in Central Nepal Himalaya

The Mugling–Narayanghat road section falls within the Lesser Himalaya and Siwalik zones of Central Nepal Himalaya and is highly deformed by the presence of numerous faults and folds. Over the years, this road section and its surrounding area have experienced repeated landslide activities. For that reason, landslide susceptibility zonation is essential for roadside slope disaster management and for planning further development activities. The main goal of this study was to investigate the application of the frequency ratio (FR), statistical index (SI), and weights-of-evidence (WoE) approaches for landslide susceptibility mapping of this road section and its surrounding area. For this purpose, the input layers of the landslide conditioning factors were prepared in the first stage. A landslide inventory map was prepared using earlier reports, aerial photographs interpretation, and multiple field surveys. A total of 438 landslide locations were detected. Out these, 295 (67 %) landslides were randomly selected as training data for the modeling using FR, SI, and WoE models and the remaining 143 (33 %) were used for the validation purposes. The landslide conditioning factors considered for the study area are slope gradient, slope aspect, plan curvature, altitude, stream power index, topographic wetness index, lithology, land use, distance from faults, distance from rivers, and distance from highway. The results were validated using area under the curve (AUC) analysis. From the analysis, it is seen that the FR model with a success rate of 76.8 % and predictive accuracy of 75.4 % performs better than WoE (success rate, 75.6 %; predictive accuracy, 74.9 %) and SI (success rate, 75.5 %; predictive accuracy, 74.6 %) models. Overall, all the models showed almost similar results. The resultant susceptibility maps can be useful for general land use planning.     

Landslide susceptibility mapping at Zonouz Plain,Iran using genetic programming and comparison with frequency ratio,logistic regression,and artificial neural network models

Without a doubt, landslide is one of the most disastrous natural hazards and landslide susceptibility maps (LSMs) in regional scale are the useful guide to future development planning. Therefore, the importance of generating LSMs through different methods is popular in the international literature. The goal of this study was to evaluate the susceptibility of the occurrence of landslides in Zonouz Plain, located in North-West of Iran. For this purpose, a landslide inventory map was constructed using field survey, air photo/satellite image interpretation, and literature search for historical landslide records. Then, seven landslide-conditioning factors such as lithology, slope, aspect, elevation, land cover, distance to stream, and distance to road were utilized for generation LSMs by various models: frequency ratio (FR), logistic regression (LR), artificial neural network (ANN), and genetic programming (GP) methods in geographic information system (GIS). Finally, total four LSMs were obtained by using these four methods. For verification, the results of LSM analyses were confirmed using the landslide inventory map containing 190 active landslide zones. The validation process showed that the prediction accuracy of LSMs, produced by the FR, LR, ANN, and GP, was 87.57, 89.42, 92.37, and 93.27 %, respectively. The obtained results indicated that the use of GP for generating LSMs provides more accurate prediction in comparison with FR, LR, and ANN. Furthermore; GP model is superior to the ANN model because it can present an explicit formulation instead of weights and biases matrices.     

Regional-scale landslide activity and landslide susceptibility zonation in the Nepal Himalaya

Landslide susceptibility zonation mapping is a fundamental procedure for geo-disaster management in tropical and sub-tropical regions. Recently, various landslide susceptibility zonation models have been introduced in Nepal with diverse approaches of assessment. However, validation is still a problem. Additionally, the role of various predisposing causative parameters for landslide activity is still not well understood in the Nepal Himalaya. To address these issues of susceptibility zonation and landslide activity, about 4,000 km² area of central Nepal was selected for regional-scale assessment of landslide activity and susceptibility zonation mapping. In total, 655 new landslides and 9,229 old landslides were identified with the study area with the help of satellite images, aerial photographs, field data and available reports. The old landslide inventory was “blind landslide database” and could not explain the particular rainfall event responsible for the particular landslide. But considering size of the landslide, blind landslide inventory was reclassified into two databases: short-duration high-intensity rainfall-induced landslide inventory and long-duration low-intensity rainfall-induced landslide inventory. These landslide inventory maps were considered as proxy maps of multiple rainfall event-based landslide inventories. Similarly, all 9,884 landslides were considered for the activity assessment of predisposing causative parameters. For the Nepal Himalaya, slope, slope aspect, geology and road construction activity (anthropogenic cause) were identified as most affective predisposing causative parameters for landslide activity. For susceptibility zonation, multivariate approach was considered and two proxy rainfall event-based landslide databases were used for the logistic regression modelling, while a relatively recent landslide database was used in validation. Two event-based susceptibility zonation maps were merged and rectified to prepare the final susceptibility zonation map and its prediction rate was found to be more than 82 %. From this work, it is concluded that rectification of susceptibility zonation map is very appropriate and reliable. The results of this research contribute to a significant improvement in landslide inventory preparation procedure, susceptibility zonation mapping approaches as well as role of various predisposing causative parameters for the landslide activity.     

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.     

Prediction of rainfall using artificial neural networks for synoptic station of Mashhad: a case study

In this paper, we have utilized ANN (artificial neural network) modeling for the prediction of monthly rainfall in Mashhad synoptic station which is located in Iran. To achieve this black-box model, we have used monthly rainfall data from 1953 to 2003 for this synoptic station. First, the Hurst rescaled range statistical (R/S) analysis is used to evaluate the predictability of the collected data. Then, to extract the rainfall dynamic of this station using ANN modeling, a three-layer feed-forward perceptron network with back propagation algorithm is utilized. Using this ANN structure as a black-box model, we have realized the complex dynamics of rainfall through the past information of the system. The approach employs the gradient decent algorithm to train the network. Trying different parameters, two structures, M₅₃₁ and M₇₄₁, have been selected which give the best estimation performance. The performance statistical analysis of the obtained models shows with the best tuning of the developed monthly prediction model the correlation coefficient (R), root mean square error (RMSE), and mean absolute error (MAE) are 0.93, 0.99, and 6.02 mm, respectively, which confirms the effectiveness of the developed models.     

Landslide hazard and risk assessment mapping of mountainous terrains - a case study from Kumaun Himalaya, India

Landslides are studied systematically in order to evaluate the nature of hazard and the damages to the human life, land, roads, buildings and other properties. This can be expressed in terms of risk, which is a function of hazard probability and damage potential. A risk map will indicate the priorities for landslide hazard management. A new approach to risk assessment mapping using a risk assessment matrix (RAM) is presented.     

Evaluation of the consistency of landslide susceptibility mapping: a case study from the Kankai watershed in east Nepal

GIS-based landslide susceptibility maps for the Kankai watershed in east Nepal are developed using the frequency ratio method and the multiple linear regression technique. The maps are derived from comparing observed landslides with possible causative factors: slope angle, slope aspect, slope curvature, relative relief, distance from drainage, land use, geology, distance from faults and mean annual rainfall. The consistency of the maps is evaluated using landslide density analysis, success rate analysis and spatially agreed area approach. The first two analyses produce almost identical quantitative results, whereas the last approach is able to reveal spatial differences between the maps and also to improve predictions in the agreed high landslide-susceptible area.     

Seismotectonics of the Nepal Himalaya from a local seismic network

The National Seismological Network of Nepal consists of 17 short period seismic stations operated since 1994. It provides an exceptional view of the microseismic activity over nearly one third of the Himalayan arc, including the only segment, between longitudes 78°E and 85°E, that has not produced any M>8 earthquakes over the last century. It shows a belt of seismicity that follows approximately the front of the Higher Himalaya with most of the seismic moment being released at depths between 10 and 20 km. This belt of seismicity is interpreted to reflect interseismic stress accumulation in the upper crust associated with creep in the lower crust beneath the Higher Himalaya. The seismic activity is more intense around 82°E in Far-Western Nepal and around 87°E in Eastern Nepal. Western Nepal, between 82.5 and 85°E, is characterized by a particularly low level of seismic activity. We propose that these lateral variations are related to segmentation of the Main Himalayan Thrust Fault. The major junctions between the different segments would thus lie at about 87°E and 82°E with possibly an intermediate one at about 85°E. These junctions seem to coincide with some of the active normal faults in Southern Tibet. Lateral variation of seismic activity is also found to correlate with lateral variations of geological structures suggesting that segmentation is a long-lived feature. We infer four 250–400 km long segments that could produce earthquakes comparable to the M=8.4 Bihar–Nepal earthquake that struck eastern Nepal in 1934. Assuming the model of the characteristic earthquake, the recurrence interval between two such earthquakes on a given segment is between 130 and 260 years.     

Landslide occurrence and its relation with terrain factors in the Siwalik Hills,Nepal: case study of susceptibility assessment in three basins

Landslides in the three studied basins of the Siwalik Hills are not random in distribution; they tend to cluster in certain areas implying the control of certain in situ factors or their combination. Landslide controlling in situ factors were reviewed and analyzed from maps, aerial photos and imageries using GIS. Chi square analysis was carried out to test the significance of landslide distribution vis-à-vis in situ factors. Slope gradient and relative relief were consistently significant in landslide distribution. Geology, dip-topography relation, land use and land cover, and vegetation conditions appeared important in landslide occurrence in all three basins either in terms of area or count in any two basins. Slope aspect and altitude tested significant for landslide occurrence in at least two basins. However, upslope flow contributing area, drainage density and distance to lineament were found insignificant in all three basins. In situ factors that tested significant in any two basins were used for landslide susceptibility analysis using a bivariate statistical approach. The distribution of landslides strongly correlates with susceptibility indices. With in situ factors, landslide susceptibility had good correlation with slope gradient and relative relief. Incorporating calculated factor weight values from one basin to the other two basins, proxy susceptibility index maps were also prepared. A moderate to good positive correlation appeared between them implying certain range of confidence for replicating the result for whole of the Siwalik Hills. Slope gradient and relative relief can be used as proxy indicators of landslide susceptible areas in the Siwalik Hills.     

Probability-weighted hazard maps for comparing different flood risk management strategies: a case study

The study proposes an original methodology for producing probability-weighted hazard maps based on an ensemble of numerical simulations. These maps enable one to compare different strategies for flood risk management. The methodology was applied over a 270-km² flood-prone area close to the left levee system of a 28-km reach of the river Reno (Northern Central Italy). This reach is characterised by the presence of a weir that allows controlled flooding of a large flood-prone area during major events. The proposed probability-weighted hazard maps can be used to evaluate how a structural measure such as the mentioned weir alters the spatial variability of flood hazard in the study area. This article shows an application by constructing two different flood hazard maps: a first one which neglects the presence of the weir using a regular levee system instead, and a second one that reflects the actual geometry with the weir. Flood hazard maps were generated by combining the results of several inundation scenarios, simulated by coupling 1D- and 2D-hydrodynamic models.     

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.     

An assessment on the use of logistic regression and artificial neural networks with different sampling strategies for the preparation of landslide susceptibility maps

The main purpose of this study is to highlight the conceptual differences of produced susceptibility models by applying different sampling strategies: from all landslide area with depletion and accumulation zones and from a zone which almost represents pre-failure conditions. Variations on accuracy and precision values of the models constructed considering different algorithms were also investigated. For this purpose, two most popular techniques, logistic regression analysis and back-propagation artificial neural networks were taken into account. The town Ispir and its close vicinity (Northeastern part of Turkey), suffered from landsliding for many years was selected as the application site of this study. As a result, it is revealed that the back-propagation artificial neural network algorithms overreact to the samplings in which the presence (1) data were taken from the landslide masses. When the generalization capacities of the models are taken into consideration, these reactions cause imprecise results, even though the area under curve (AUC) values are very high (0.915 < AUC < 0.949). On the other hand, the susceptibility maps, based on the samplings in which the presence (1) data were taken from a zone which almost represents pre-failure conditions constitute more realistic susceptibility evaluations. However, considering the spatial texture of the final susceptibility values, the maps produced using the outputs of the back-propagation artificial neural networks could be interpreted as highly optimistic, while of those generated using the resultant probabilities of the logistic regression equations might be evaluated as pessimistic. Consequently, it is evident that, there are still some needs for further investigations with more realistic validations and data to find out the appropriate accuracy and precision levels in such kind of landslide susceptibility studies.     

Use of backpropagation neural network for landslide monitoring: a case study in the higher Himalaya

Static and dynamic rock slope stability analyses were performed using a numerical discontinuum modelling technique for a 700-m high rock slope in western Norway. The rock slope has been investigated by the Geological Survey of Norway (NGU), which has been carrying out rock slide studies for the county Møre and Romsdal in western Norway. The purpose of numerical modelling was to estimate the volume of the rock mass that could potentially slide under static and dynamic forces. This estimation was required to assess the run-up heights (tsunami) in a fjord that could potentially be caused by the rockslide. Three cases have been simulated for predicting the behaviour of the rock slope. First, an initial static loading is applied in the numerical model to simulate the prevailing rock mass conditions at the site. Second, saturated and weathered joint conditions are modelled by reducing the residual friction angle along the discontinuities of the rock mass. In doing so, the model simulates the effect of degradation of discontinuities in the rock slope. Third, a dynamic loading, based on peak ground accelerations expected in the area, is applied to simulate dynamic earthquake conditions.

These numerical studies have provided some useful insights into the deformation mechanisms in the rock slope. Both sliding and rotation of blocks start to occur once the residual friction angle along the discontinuities is reduced and when the region is shaken by a strong earthquake. The results indicate that, due to variations in the inclination of discontinuities, the entire slope does not become unstable and that down-slope sliding and rotation of blocks occur mainly on the top layers of the slope. Within the range of parameter values considered for this study, it is unlikely that the whole rock slope can be destabilised. The study provides an illustration of how the geo-mechanical properties of a rock mass can be integrated in a discontinuum rock slope model, which is used for predicting the behaviour of the slope under existing environmental and earthquake conditions. This model has helped not only to better understand the dynamics of the rockslide but also to estimate the potential rock volume that can become unstable when subjected to static and dynamic loads.