Integrating expert knowledge with statistical analysis for landslide susceptibility assessment at regional scale

Abstract: In this paper, an integration landslide susceptibility model by combining expert-based and bivariate statistical analysis (Landslide Susceptibility Index—LSI) approaches is presented. Factors related with the occurrence of landslides—such as elevation, slope angle, slope aspect, lithology, land cover, Mean Annual Precipitation (MAP) and Peak Ground Acceleration (PGA)—were analyzed within a GIS environment. This integrated model produced a landslide susceptibility map which categorized the study area according to the probability level of landslide occurrence. The accuracy of the final map was evaluated by Receiver Operating Characteristics (ROC) analysis depending on an independent (validation) dataset of landslide events. The prediction ability was found to be 76% revealing that the integration of statistical analysis with human expertise can provide an acceptable landslide susceptibility assessment at regional scale

GIS-Based landslide susceptibility modeling: a comparison between best-first decision tree and its two ensembles (BagBFT and RFBFT)

This study aimed to explore and compare the application of current state-of-the-art machine learning techniques, including bagging (Bag) and rotation forest (RF), to assess landslide susceptibility with the base classifier best-first decision tree (BFT). The proposed two novel ensemble frameworks, BagBFT and RFBFT, and the base model BFT, were used to model landslide susceptibility in Zhashui County (China), which suffers from landslides. Firstly, we identified 169 landslides through field surveys and image interpretation. Then, a landslide inventory map was built. These 169 historical landslides were randomly classified into two groups: 70% for training data and 30% for validation data. Then, 15 landslide conditioning factors were considered for mapping landslide susceptibility. The three ensemble outputs were estimated with a receiver operating characteristic (ROC) curve and statistical tests, as well as a new approach, the improved frequency ratio accuracy. The areas under the ROC curve (AUCs) for the training data (success rate) of the three algorithms were 0.722 for BFT, 0.869 for BagBFT, and 0.895 for RFBFT. The AUCs for the validating groups (prediction rates) were 0.718, 0.834, and 0.872, respectively. The frequency ratio accuracy of the three models was 0.76163 for the BFT model, 0.92220 for the BagBFT model, and 0.92224 for the RFBFT model. Both BagBFT and RFBFT ensembles can improve the accuracy of the BFT base model, and RFBFT was relatively better. Therefore, the RFBFT model is the most effective approach for the accurate modeling of landslide susceptibility mapping (LSM). All three models can improve the identification of landslide-prone areas, enhance risk management ability, and afford more detailed information for land-use planning and policy setting.National Natural Science Foundation of China | Ref. 41977228Key Research Program of Shaanxi | Ref. 2022SF-33

GIS Modeling of the Prominent Geohazards in Arkansas

The State of Arkansas is prone to numerous geohazards. This thesis is a twofold study of prominent geohazards in Arkansas: the first fold includes a novel triggerless approach for mass wasting susceptibility modeling applied to the Boston Mountains in NW Arkansas, and the second fold is a GIS-based regression modeling of the extreme weather patterns at the state level. Each study fold is presented in this thesis as a separate chapter embracing a published peer-reviewed paper. In the first paper, I have used the analytical hierarchy process to assign preliminary statistical weights to the most cogent variables influencing mass wasting in the central Boston Mountains. These most significant variables are then incorporated in Fuzzy modeling of mass wasting susceptibility within the 1200 km2 study area. For comparison and accuracy assessment, a second model has been established using a conventional weighted overlay (WO) approach. Results indicate that the developed novel approach is superior, with approximately 83% accuracy, to the traditional WO approach that has a marginal success of about 28% accuracy. Road related mass wasting events recorded by the Arkansas Department of Transportation have been used to validate both models. In the second paper, I have conducted a systematically gridded analysis of severe weather events, including tornadoes, derechos, and hail, during 1955-2015. The study examines and statistically determines the most significant explanatory variables contributing to the spatial patterns of severe weather events between 1955 and 2015, consequently it identifies severity indices for the entire state. These weather-related hazards and their associated risk will always abide; therefore, the best defense is employ geospatial technologies to plan for hazard mitigation. The mass wasting model developed in this study contributes pivotal information for identifying zones of high risk along roadways in NW Arkansas, which definitely can be adapted to avoid disastrous road failures. In addition, the weather-related severity indices determined at the state level can profoundly benefit state and federal agencies focused on increasing the availability of public and private storm shelters in previously under-represented zones of high risk. This undoubtedly will save lives from unavoidable catastrophic events across the entire state

Landslide susceptibility mapping using certainty factor, index of entropy and logistic regression models in GIS and their comparison at Mugling-Narayanghat road section in Nepal Himalaya

Landslide susceptibility maps are vital for disaster management and for planning development activities in the mountainous country like Nepal. In the present study, landslide susceptibility assessment of Mugling-Narayanghat road and its surrounding area is made using bivariate (certainty factor and index of entropy) and multivariate (logistic regression) models. At first, a landslide inventory map was prepared using earlier reports and aerial photographs as well as by carrying out field survey. As a result, 321 landslides were mapped and out of which 241 (75 %) were randomly selected for building landslide susceptibility models, while the remaining 80 (25 %) were used for validating the models. The effectiveness of landslide susceptibility assessment using GIS and statistics is based on appropriate selection of the factors which play a dominant role in slope stability. In this case study, the following landslide conditioning factors were evaluated: slope gradient; slope aspect; altitude; plan curvature; lithology; land use; distance from faults, rivers and roads; topographic wetness index; stream power index; and sediment transport index. These factors were prepared from topographic map, drainage map, road map, and the geological map. Finally, the validation of landslide susceptibility map was carried out using receiver operating characteristic (ROC) curves. The ROC plot estimation results showed that the susceptibility map using index of entropy model with AUC value of 0.9016 has highest prediction accuracy of 90.16 %. Similarly, the susceptibility maps produced using logistic regression model and certainty factor model showed 86.29 and 83.57 % of prediction accuracy, respectively. Furthermore, the ROC plot showed that the success rate of all the three models performed more than 80 % accuracy (i.e. 89.15 % for IOE model, 89.10 % for LR model and 87.21 % for CF model). Hence, it is concluded that all the models employed in this study showed reasonably good accuracy in predicting the landslide susceptibility of Mugling-Narayanghat road section. These landslide susceptibility maps can be used for preliminary land use planning and hazard mitigation purpose.ArticleNATURAL HAZARDS. 65(1):135-165 (2013)journal articl

Evaluation of the landslide susceptibility and its spatial difference in the whole Qinghai-Tibetan Plateau region by five learning algorithms

AbstractLandslides are considered as major natural hazards that cause enormous property damages and fatalities in Qinghai-Tibetan Plateau (QTP). In this article, we evaluated the landslide susceptibility, and its spatial differencing in the whole Qinghai-Tibetan Plateau region using five state-of-the-art learning algorithms; deep neural network (DNN), logistic regression (LR), Naïve Bayes (NB), random forest (RF), and support vector machine (SVM), differing from previous studies only in local areas of QTP. The 671 landslide events were considered, and thirteen landslide conditioning factors (LCFs) were derived for database generation, including annual rainfall, distance to drainage $${(\mathrm{Ds}}_{\mathrm{d}})$$ ( Ds d ) , distance to faults $${(\mathrm{Ds}}_{\mathrm{f}})$$ ( Ds f ) , drainage density ($${D}_{d})$$ D d ) , elevation (Elev), fault density $$({F}_{d})$$ ( F d ) , lithology, normalized difference vegetation index (NDVI), plan curvature $${(\mathrm{Pl}}_{\mathrm{c}})$$ ( Pl c ) , profile curvature $${(\mathrm{Pr}}_{\mathrm{c}})$$ ( Pr c ) , slope $${(S}^{^\circ })$$ ( S ∘ ) , stream power index (SPI), and topographic wetness index (TWI). The multi-collinearity analysis and mean decrease Gini (MDG) were used to assess the suitability and predictability of these factors. Consequently, five landslide susceptibility prediction (LSP) maps were generated and validated using accuracy, area under the receiver operatic characteristic curve, sensitivity, and specificity. The MDG results demonstrated that the rainfall, elevation, and lithology were the most significant landslide conditioning factors ruling the occurrence of landslides in Qinghai-Tibetan Plateau. The LSP maps depicted that the north-northwestern and south-southeastern regions ( 45% of total area). Moreover, among the five models with a high goodness-of-fit, RF model was highlighted as the superior one, by which higher accuracy of landslide susceptibility assessment and better prone areas management in QTP can be achieved compared to previous results. Graphical Abstrac

Spatial prediction of landslide susceptibility/intensity through advanced statistical approaches implementation: applications to the Cinque Terre (Eastern Liguria, Italy)

Landslides are frequently responsible for considerable huge economic losses and casualties in mountainous regions especially nowadays as development expands into unstable hillslope areas under the pressures of increasing population size and urbanization (Di Martire et al. 2012). People are not the only vulnerable targets of landslides. Indeed, mass movements can easily lay waste to everything in their path, threatening human properties, infrastructures and natural environments. Italy is severely affected by landslide phenomena and it is one of the most European countries affected by this kind of phenomena. In this framework, Italy is particularly concerned with forecasting landslide effects (Calcaterra et al. 2003b), in compliance with the National Law n. 267/98, enforced after the devastating landslide event of Sarno (Campania, Southern Italy). According to the latest Superior Institute for the Environmental Protection and Research (ISPRA, 2018) report on "hydrogeological instability" of 2018, it emerges that the population exposed to landslides risk is more than 5 million and in particular almost half-million falls into very high hazard zones. The slope stability can be compromised by both natural and human-caused changes in the environment. The main reasons can be summarised into heavy rainfalls, earthquakes, rapid snow-melts, slope cut due to erosions, and variation in groundwater levels for the natural cases whilst slopes steepening through construction, quarrying, building of houses, and farming along the foot of mountainous zone correspond to the human component. This Ph.D. thesis was carried out in the Liguria region, inside the Cinque Terre National Park. This area was chosen due to its abundance of different types of landslides and its geological, geomorphological and urban characteristics. The Cinque Terre area can be considered as one of the most representative examples of human-modified landscape. Starting from the early centuries of the Middle Ages, local farmers have almost completely modified the original slope topography through the construction of dry-stone walls, creating an outstanding terraced coastal landscape (Terranova 1984, 1989; Terranova et al. 2006; Brandolini 2017). This territory is extremely dynamic since it is characterized by a complex geological and geomorphological setting, where many surficial geomorphic processes coexist, along with peculiar weather conditions (Cevasco et al. 2015). For this reason, part of this research focused on analyzing the disaster that hit the Cinque Terre on October, 25th, 2011. Multiple landslides took place in this occasion, triggering almost simultaneously hundreds of shallow landslides in the time-lapse of 5-6 hours, causing 13 victims, and severe structural and economic damage (Cevasco et al. 2012; D\u2019Amato Avanzi et al. 2013). Moreover, this artificial landscape experienced important land-use changes over the last century (Cevasco et al. 2014; Brandolini 2017), mostly related to the abandonment of agricultural activity. It is known that terraced landscapes, when no longer properly maintained, become more prone to erosion processes and mass movements (Lesschen et al. 2008; Brandolini et al. 2018a; Moreno-de-las-Heras et al. 2019; Seeger et al. 2019). Within the context of slope instability, the international community has been focusing for the last decade on recognising the landslide susceptibility/hazard of a given area of interest. Landslide susceptibility predicts "where" landslides are likely to occur, whereas, landslide hazard evaluates future spatial and temporal mass movement occurrence (Guzzetti et al., 1999). Although both definitions are incorrectly used as interchangeable. Such a recognition phase becomes crucial for land use planning activities aimed at the protection of people and infrastructures. In fact, only with proper risk assessment governments, regional institutions, and municipalities can prepare the appropriate countermeasures at different scales. Thus, landslide susceptibility is the keystone of a long chain of procedures that are actively implemented to manage landslide risk at all levels, especially in vulnerable areas such as Liguria. The methods implemented in this dissertation have the overall objective of evaluating advanced algorithms for modeling landslide susceptibility. The thesis has been structured in six chapters. The first chapter introduces and motivates the work conducted in the three years of the project by including information about the research objectives. The second chapter gives the basic concepts related to landslides, definition, classification and causes, landslide inventory, along with the derived products: susceptibility, hazard and risk zoning, with particular attention to the evaluation of landslide susceptibility. The objective of the third chapter is to define the different methodologies, algorithms and procedures applied during the research activity. The fourth chapter deals with the geographical, geological and geomorphological features of the study area. The fifth chapter provides information about the results of the applied methodologies to the study area: Machine Learning algorithms, runout method and Bayesian approach. Furthermore, critical discussions on the outcomes obtained are also described. The sixth chapter deals with the discussions and the conclusions of this research, critically analysing the role of such work in the general panorama of the scientific community and illustrating the possible future perspectives

Landslide Susceptibility Assessment Using Bivariate Method: A Case Study from River Neelum and Jehlum Catchment Area

Landslide is a frequently occurring natural calamity in the northern areas of Pakistan. The current study is aimed to assess the susceptibility of landslide hazard to highlight the vulnerable areas for the purpose of risk reduction along Neelum and Jehlum rivers in district Muzaffarabad. A data-driven predictive approach was adopted to conduct this study by using Weight of Evidence (WOE) model along with eleven conditioning factors. A spatial distribution map of landslides was prepared using orthophoto, previous records, and derivatives (hill shad, topographic openness, slope, aspect, curvature) of Digital Elevation Model (DEM). The results show that the roads, lithology, and rivers are the most important triggering factors for landslides in both valleys. Approximately 30% of the area is under low susceptibility zone in Jhelum valley while only 13% of the area falls under low susceptibility zone in Neelum river valley. In Jhelum river valley the medium susceptibility zone covers 35% of the total area whereas, Neelam river valley has 26% of the total area under medium susceptibility zone. Around 61% of the land in the Neelam river valley and 35% of the land in the Jhelum river valley are under high susceptibility zone. The area under high hazard lies in the north-east of the district due to multiple conducive factors to trigger landslides including weak lithology (mudstone, sandstone, shales, and clays), high altitude along steep slopes and excessive precipitation (1800 mm/ year). Furthermore, the high hazard zone in study area is not suitable for construction purpose but was suitable for plantation. The validation result (89.41%) is justifying the performance of this model

Evaluation of multi-hazard map produced using MaxEnt machine learning technique

Natural hazards are diverse and uneven in time and space, therefore, understanding its complexity is key to save human lives and conserve natural ecosystems. Reducing the outputs obtained after each modelling analysis is key to present the results for stakeholders, land managers and policymakers. So, the main goal of this survey was to present a method to synthesize three natural hazards in one multi-hazard map and its evaluation for hazard management and land use planning. To test this methodology, we took as study area the Gorganrood Watershed, located in the Golestan Province (Iran). First, an inventory map of three different types of hazards including flood, landslides, and gullies was prepared using field surveys and different official reports. To generate the susceptibility maps, a total of 17 geo-environmental factors were selected as predictors using the MaxEnt (Maximum Entropy) machine learning technique. The accuracy of the predictive models was evaluated by drawing receiver operating characteristic-ROC curves and calculating the area under the ROC curve-AUCROC. The MaxEnt model not only implemented superbly in the degree of fitting, but also obtained significant results in predictive performance. Variables importance of the three studied types of hazards showed that river density, distance from streams, and elevation were the most important factors for flood, respectively. Lithological units, elevation, and annual mean rainfall were relevant for detecting landslides. On the other hand, annual mean rainfall, elevation, and lithological units were used for gully erosion mapping in this study area. Finally, by combining the flood, landslides, and gully erosion susceptibility maps, an integrated multi-hazard map was created. The results demonstrated that 60% of the area is subjected to hazards, reaching a proportion of landslides up to 21.2% in the whole territory. We conclude that using this type of multi-hazard map may be a useful tool for local administrators to identify areas susceptible to hazards at large scales as we demonstrated in this research