Informative landslide hazard estimates are needed to support landslide mitigation strategies to reduce landslide risk across the United States (U.S.). While existing national-scale landslide susceptibility products assess where landslides are likely to occur, they do not address how often, which is a critical element of landslide hazard and risk assessments. In particular, the U.S. Federal Emergency Management Agency’s National Risk Index (NRI) requires landslide frequency estimates by county, which are U.S. administrative regions ranging from 120 km2 to 377,055 km2 in size, to inform expected annual loss estimates. In this study, we present county-level landslide frequency (landslides area-1 y-1) estimates for the 50 U.S. states. We applied Bayesian negative binomial regression to estimate both the expected (average) landslide frequency and full distribution of annual landslide counts for each county as a function of landslide susceptible area, frequency of potentially triggering precipitation, and propensity for triggering earthquakes. We trained our model with 62,720 reported landslides from 316 counties with the most comprehensive records available nationwide and used zero-inflated negative binomial distributions as an incompleteness model to correct for temporal reporting gaps. We found that average annual landslide frequencies vary by nearly three orders of magnitude across U.S. counties, ranging from 0.05 (0.04–0.07) landslides 1000 km-2y-1 in Midland County, Texas to 31 (21–43) landslides 1000 km-2y-1 in Lake County, California and reflecting the country’s strong variations in landslide susceptibility, earthquake probability, and precipitation frequency. Counties with estimated frequencies in the top 20% of all counties are predominately along the West Coast of the continental United States, in mountainous regions of the Pacific Northwest and Intermountain West, in locally steep or earthquake prone regions of the Midwest and South, along the Appalachians, in southern Alaska, and on the big island of Hawaii. By examining the number of landslides predicted in 99th percentile years for each county, we identified that 31% of U.S. counties have potential for widespread landsliding, even when such large events have not been reported in the training data for that county. Overall, our results better represent the range of possible landslide frequencies and spatial variations across the entire United States than previous national-scale estimates reported in the NRI and can inform other risk reduction and loss mitigation efforts across the United States.
How to cite: Luna, L., Woodard, J., Bytheway, J., Belair, G., and Mirus, B.: Constraining landslide frequency across the United States to inform county-level risk reduction, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-1724, https://doi.org/10.5194/egusphere-egu25-1724, 2025.
In May 2024, a Mega Disaster hit 96% of the municipalities in Rio Grande do Sul (RS) state, southern Brazil, causing 182 casualties and impacting approximately 2.4 million people. In addition to the floods that hit the capital Porto Alegre, more than 15,000 landslides were recorded due to the extreme rainfall event (> 600 mm in some regions), severely impacting an area of nearly 18,000 km². Although other landslide events have been recorded in RS in the past, none of them have matched the magnitude of this one. In this sense, we aimed to generate a landslide susceptibility model, based on historical data and evaluate its capacity to forecast the areas affected by landslides in 2024. This retrospective assessment was performed using an inventory of four past events between 1995 and 2017, totalling 1,211 landslides, represented by 15,580 points. We randomly selected the same number of points (15,580) over the RS to represent non-landslide areas and split the entire sample set into training (70%) and validation (30%). A Random Forest model, leveraging seven morphometric parameters, was employed to generate the map, which was evaluated with the validation sample set. A second validation was carried out considering the landslides in 2024, represented by 324,500 points. This validation was based on the relationship closeness between 2024 landslides and each susceptibility class using frequency ratio. The last evaluation consisted of analysing landslide areas (rupture, propagation, and deposition) and their distribution in each susceptibility class. To achieve this, we automatically divided the 2024 landslide points into three sets, according to the altitude difference found in each polygon. Our landslide susceptibility map presented a high performance, with an overall accuracy of 0.9, being capable of correctly classifying 64% of 2024 landslides into susceptible areas (very high, high, and moderate susceptibility classes). The very high susceptibility class accounted for 31% of the 2024 landslides and had a frequency ratio of 13.04, showing a high correlation between landslide locations and the analysed class. Further analysis revealed that the model successfully predicted 79% of rupture zones, highlighting its robustness in identifying key prone areas. While the model performed well in identifying rupture and propagation areas as susceptible, its predictions for deposition zones were less accurate, likely due to limitations in the historical inventory, which was carried out after 2017, when most landslide deposition areas were no longer visible in remote sensing imagery. Furthermore, even though the 2024 Mega Disaster was responsible for 12.5 times more landslides than all the previous inventory, our model based on 1,211 landslides correctly classified around 9,600 landslides (64%) in susceptible areas. Therefore, although the 2024 extreme rainfall event was much larger than any previously recorded in the region, many areas could have already been identified as susceptible. Finally, the existence of a more complete landslide inventory (including rupture, propagation, and deposition areas) provides more accurate susceptibility maps, which can support territorial planning, contributing to disaster risk management, mitigation strategies, and land use policies.
How to cite: Quevedo, R. P., Maciel, D. A., Andrades Filho, C. O., Mexias, L. F. S., Oliveira, G. G., Herrmann, P. B., Alves, F. C., and Glade, T.: Could such a large landslide event be expected in Rio Grande do Sul, southern Brazil? Using past events to predict the area impacted by the 2024 Mega Disaster, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-15304, https://doi.org/10.5194/egusphere-egu25-15304, 2025.
Estimating the potential runout distance of landslides and their associated impacted areas is a critical component of landslide hazard and risk analysis. Traditionally, back-analysis of past landslides has been employed to predict the runout behaviour of potential future events. To refine landslide runout models and characterize co-seismic landslide dynamics, we conducted an in-depth analysis of a subset of landslides triggered by the Mw 7.8 Kaikōura earthquake in New Zealand (14 November 2016), focusing on the Kowhai Valley in Kaikōura.
First, we mapped polylines connecting landslide sources to their corresponding deposits. Given that all landslides were triggered during the same seismic event within steep upland catchments, source areas did not consistently correspond directly to mapped debris trails. Second, we attributed these polylines with information on confinement, substrate type, connectivity, geometry, and physiographic attributes, analysing their relationships with travel length and fall height to identify controls on runout distance. Third, we applied three regional-scale runout modelling approaches—1) a Fahrböschung angle method, 2) the Gravitational Path Process Model, and 3) Flow-R—to evaluate their effectiveness in predicting travel distances and patterns of co-seismic landslide runout.
Our mapping identified 3,535 landslide polylines linking 3,105 source areas to 2,652 debris trails. Approximately two-thirds of the landslides exhibited a one-to-one relationship between source and deposit, while the remainder displayed more complex linkages, including multiple deposits from a single source, single deposits from multiple sources, or interactions involving multiple sources and deposits. Statistical analysis revealed significant relationships between runout distance and factors such as substrate type, confinement, coupling, and geometry, although no significant relationship was observed with landslide volume.
Model accuracy assessments, using goodness of fit metrics, showed that most approaches either displayed weak accuracy or overestimated landslide runout areas. The best fit models indicated that the landslides triggered in the Kaikōura earthquake travelled a shorter distance than expected from the international literature. Further analysis revealed considerable variability in model accuracy for individual landslides, with larger landslides showing better goodness-of-fit metrics than smaller ones. Landslides located in the lower reaches of the Kowhai Valley also demonstrated higher model accuracy, potentially as a function of landscape relief. These findings underscore the complex controls influencing co-seismic landslide runout and highlight the importance of accounting for uncertainties in regional-scale landslide runout models.
How to cite: de Vilder, S., Wolter, A., Lukovic, B., Leith, K., Hill, S., and Cox, S.: Runout characteristics of landslides triggered by the 2016 Kaikoura Earthquake, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-14509, https://doi.org/10.5194/egusphere-egu25-14509, 2025.
Rapid, giant landslides, or sturzstroms, are among the most powerful natural hazards on Earth. They are produced by catastrophic, deep-seated slope collapses with minimum volumes on the order of 10⁶–10⁷ m³. Such collapses are often the final stage of accelerating slope movement that may have continued for years. Measurements made over 60 years ago before the failure of Mt. Toc into the Vajont reservoir in the Italian Alps remain one of the best records of pre-collapse slope movement. Numerous studies have recognized that the rate of movement increased hyperbolically during at least two months of heavy rainfall before the mountainside collapsed on 9 October 1963. Two hundred million m³ of rock sent a wave of water over the Vajont dam, killing approximately 2,500 people in the downstream communities of Longarone, Pirago, Villanova, Rivalta, and Fae. Analysis of the extended record shows that the hyperbolic trend was preceded by an exponential acceleration during 1962. The earlier trend was interrupted in December 1962 when the reservoir was temporarily drained to install engineering safety measures. The acceleration resumed in July–August 1963 after the reservoir was refilled to its pre-drainage level. This combined exponential-hyperbolic acceleration trend is consistent with the activation and eventual linkage of cracks along the future failure plane.This suggests that the surface movements were a consequence of fracturing as deep as 200 m underground, rather than cracking being a result of slope movement. This interpretation points to the weakening of deep rock as the primary driver of failure, caused by factors such as increases in pore water pressure and water-induced corrosion, rather than the destabilization from the weight of a water-saturated slope.Since rock cracking occurs within a restricted range of physical conditions, this case study demonstrates that medium-term forecasts of catastrophic slope failure are a feasible goal. By identifying and quantifying these conditions, we can advance predictive capabilities and mitigate the devastating impacts of rapid landslides, such as tsunamis, seismic shocks, and downstream flooding.
How to cite: Pan, F. and Kilburn, C. R. J.: Controls on rates of slope movement before catatsrophic collapse, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-5374, https://doi.org/10.5194/egusphere-egu25-5374, 2025.
Landslide failures pose a severe threat to society, especially when valley bottoms become blocked, ponding rivers and burying critical infrastructure. The erratic and spatially distributed occurrence of those rapid mass wasting processes makes it eminent to understand major drivers and find reliable predictors that can help early warning.
Here, we present results of a systematic study on a progressively developing landslide near the town of Müsch, in one of the narrowest sections of the Ahr Valley, Germany. The slope instability had been reactivated by the 2021 summer flood and shows accelerated toppling and rotational movement at the 100 m wide front, as well as surface evidence of distributed movement in the 200 m long hinterland. Partial failure of the frontal sector had been modelled, indicating the damming of the 30 m wide valley bottom, causing rapid inundation of upstream settlements.
We analyse 2.5 years of continuous seismic data from a small geophone network. Seismic coda wave interferometry and resonance frequency analysis yields insights to cyclic and progressive rock stress evolution as well as the effect of water content at and below the surface. More than 3000 discrete crack emissions due to brittle rock mass failure were detected, located and quantified. The precise timing of the crack signals reveals a strong control of working time hours, suggesting an external anthropogenic forcing of the slope instability. We discuss the generic applicability of the multi-proxy seismic approach in light of further, post-flood reactivations of slope instabilities in the Ahr Valley and elsewhere.
How to cite: Dietze, M., Fracica Gonzalez, L., Bell, R., Schrott, L., and Hovius, N.: Strong daily landslide cracking activity – Does traffic drive slope failure?, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-8438, https://doi.org/10.5194/egusphere-egu25-8438, 2025.
Landslide Susceptibility Modeling (LSM) is an important method for mitigating regional landslide risks. However, the scarcity of landslide inventories and the prevalence of low-quality non-landslide samples significantly limit the further development of traditional LSM frameworks. To address this issue, this paper develops a next generation of LSM framework that redefines landslide and non-landslide samples from the perspective of deformation. By integrating deformation time-series data from Global Navigation Satellite System (GNSS) and InSAR, the framework introduces deformation samples defined by deformation rates and obtains a greater number of landslide samples and high-quality non-landslide samples through the establishment of appropriate deformation thresholds. A series of ablation experiments were conducted in Wanzhou District, Chongqing, China. The results indicate that when the deformation threshold is set to 0.6, the proposed LSM framework achieves an AUC value of 0.94, a TPR of 0.92, and a TNR of 0.94, representing a significant improvement compared to the traditional LSM framework (AUC = 0.85, TPR = 0.74, TNR = 0.58). Additionally, the experimental results demonstrate that when using InSAR data to obtain deformation samples, either a large number of low-quality InSAR deformation samples or a small number of high-quality but spatially uneven InSAR deformation samples can result in the proposed LSM framework performing worse than the traditional LSM framework. Therefore, special attention must be paid to balancing the quality and quantity of InSAR data.
How to cite: Feng, X., Du, J., Chai, B., and Bogaard, T.: Combining landslide inventories with deformation time series: A methodology, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-8480, https://doi.org/10.5194/egusphere-egu25-8480, 2025.
Landslide inventories are fundamental for susceptibility mapping, hazard modeling, and risk management. For decades, the geoscientific community has relied on manual visual interpretation of satellite and aerial imagery for landslide inventory generation. However, manual methods pose significant challenges, including subjectivity in landslide boundary delineation, limited data sharing within the scientific community, and the substantial time and expertise required for accurate mapping. Recent advancements in artificial intelligence (AI) have spurred a surge in research on semi-automated and fully automated landslide inventory mapping. Despite this progress, AI-generated inventories remain in their developmental phase, with no existing models capable of consistently producing ground-truth representations of landslide events following a triggering event. Current studies utilizing AI-based models report F1-scores ranging between 50% and 80%, with only a few achieving over 80%, often limited to the same study areas used for model training. This highlights a significant research gap in the reliability and generalizability of AI-generated inventories for hazard and risk assessments. The geoscientific community must critically assess the accuracy and transferability of AI-generated landslide data to ensure their applicability in subsequent phases of landslide response and mitigation. Further collaborative efforts and benchmark datasets are needed to establish standardized protocols for validating AI-generated landslide inventories across diverse geomorphological settings.
How to cite: Meena, S. R., Singh, S., Bhookya, R., and Floris, M.: Evaluating the Potential of AI-Generated Landslide Inventories for Hazard and Risk Management: Advancements and Limitations, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-17609, https://doi.org/10.5194/egusphere-egu25-17609, 2025.
Forecasting catastrophic slope failures is one of the most challenging tasks in landslide hazard analysis. Reliable landslide forecast is essential for civil authorities to effectively inform the public about potential mountain collapses and their timing, facilitating timely evacuations and the implementation of other safety measures. Over the past decades, great efforts have been devoted to develop and deploy high-precision monitoring technologies to observe unstable slope movements. Various empirical or physical approaches have also been proposed to forecast imminent slope collapses, the predictability of which, however, still remains elusive. One major uncertainty arises from the intermittency of geomaterial rupture behaviour, which is typically characterised by a series of progressively shorter quiescent phases interrupted by sudden accelerations, rather than a smooth continuous progression of deformation and damage. This seemingly erratic pattern complicates landslide prediction. Here, we propose a generalised failure law based on the log-periodic power law singularity model for more reliable time-to-failure forecast of catastrophic landslides. Incorporating a discrete hierarchy of time scales and rooted in the fundamental principles of statistical physics, this novel failure law accurately captures the intermittent rupture dynamics of heterogeneous geomaterials at the site scale. It ensures robustness while maintaining a strong connection to the underlying physical processes. By "locking" into the oscillatory structure of rupture dynamics, this parsimonious model transforms intermittency from traditionally perceived noise into essential information to constrain its prediction. We extensively validate this new failure law on a large dataset of 49 historical landslide events, across a wide range of contexts including rockfalls, rockslides, clayslides, and embankment slopes. The results indicate that our method is general and robust, with significant potential to mitigate landslide hazards and enhance existing early warning systems.
How to cite: Lei, Q. and Sornette, D.: A physics-based generalised failure law for forecasting catastrophic landslides, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-18378, https://doi.org/10.5194/egusphere-egu25-18378, 2025.
Forecasting landslides induced by rainfall is a challenging task that involves the interaction of multiple factors, such as soil conditions, topography, and rainfall intensity. The complex nature of these events, combined with the lack of complete data on landslide occurrences, makes it difficult to produce accurate predictions. Traditional deterministic models struggle to account for the variability and uncertainty inherent in the processes leading to landslides. On the contrary, probabilistic approaches can incorporate uncertainty and provide more reliable description of this phenomenon. In this work we develop a probabilistic framework for forecasting rainfall-induced landslide occurrence addressing the challenges of data sampling, uncertainty, and variability. For the study, we collected a dataset of shallow rainfall-induced landslides in an area in southern Italy, spanning 22 years of rainfall records. The dataset includes the locations and date of occurrences of the landslides, and daily rainfall measurements. Using a Bayesian approach, we calculate the posterior probability of landslide occurrence given specific daily cumulated rainfall thresholds. To account for the uncertainty in the landslide and rainfall data, we employed probabilistic distributions i.e., uniform and beta distributions, to model the uncertainty in the prior and likelihood functions. The uncertainty was further addressed through random sampling techniques, allowing for the integration of data variability and the dependencies between landslides and rainfall, obtaining posterior probability distributions of landslide occurrence for each rainfall threshold. The results offer a probabilistic approach to landslide forecasting that can be used for better-informed decision-making in risk management and early warning systems. By accounting for the uncertainties in the data and model parameters, our approach provides a more robust method for landslide prediction under varying rainfall conditions.
How to cite: Ferriero, F., Guzzetti, F., Urciuoli, G., and Marzocchi, W.: Bayesian probabilistic forecasting of rainfall-induced landslides, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-10443, https://doi.org/10.5194/egusphere-egu25-10443, 2025.
Landslides triggered by earthquakes evolve over time, leading to repeated damage in the affected areas. These slope movements are influenced by a range of factors, including climatic, seismic, and terrain conditions, which vary both temporally and spatially [1]. To predict the likelihood of landslides occurring across different times and locations, statistical models must account for these spatial and temporal dependencies. In this study, we employ the Markov Switching Spatiotemporal Generalized Additive Model (MSST-GAM), as introduced by Sridharan et al. [2]. Their research highlighted how this model effectively captures the spatial and temporal influences of various landslide-related factors, offering accurate susceptibility estimates for the Wenchuan area in China.
In this work, we further extend the model for hazard prediction. The model is used on a multitemporal dataset of landslides that occurred in New Zealand during and following the 2016 Kaikoura earthquake [3]. The years in which the landslides were mapped were used to separate the temporal units. Twelve covariates were used, including terrain (slope, aspect, curvature, distance from features like faults, etc.), climatic (rainfall and soil moisture), and seismic (when the year coincided with a major seismic event). We employ zero-inflated Poisson and Gaussian emission probabilities [4] for the dependent variables, which are the areas and counts of landslides in slope units. A Markov-switching GAM is used to predict the dependent variables from the covariables based on two hidden risk states (high risk and low risk). We introduce soil moisture as an additional dynamic variable to parametrize the transition probabilities between the hidden states.
We tested the model using a five-fold spatiotemporal cross-validation. The results compare favourably to a number of cross-sectional models [5]. The model predictions indicate that MSST-GAM can capture the spatial and temporal dependence of the landslide occurrences in slope units when compared with other cross-sectional and spatiotemporal models in literature.
References
[1] Keefer, D., “Investigating landslides caused by earthquakes - A historical review,” Surv. Geophys., vol. 23, no. 6, pp. 473–510, 2002
[2] Sridharan, A., Gutjahr, G., and Gopalan, S., “Markov–Switching Spatio–Temporal Generalized Additive Model for Landslide Susceptibility,” Environ. Model. Softw., vol. 173, no. August, p. 105892, Feb. 2024
[3] Bhuyan, K., Tanyaş, H., Nava, L. et al. “Generating multi-temporal landslide inventories through a general deep transfer learning strategy using HR EO data”. Sci Rep 13, 162, 2023
[4] Wagh, Y.S. and Kamalja, K.K., 2018. “Zero-inflated models and estimation in zero-inflated Poisson distribution”. Communications in Statistics-Simulation and Computation, 47(8), pp.2248-2265.
[5] Reichenbach, P., Rossi, M., Malamud, B., Mihir, M., Guzzetti, F. “A review of statistically-based landslide susceptibility models”. Earth-science reviews. 2018 May 1;180:60-91.
How to cite: Gutjahr, G., Sridharan, A., and Gopalan, S.: A Markov Switching Spatiotemporal GAM for Landslide Hazards in New Zealand, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-12309, https://doi.org/10.5194/egusphere-egu25-12309, 2025.
Rock avalanches are large volume landslides composed of flowing fragments of rock that can reach velocities in excess of 50 m/s, impact large areas, and can seriously threaten the safety of people and infrastructure. Numerical models play a crucial role in forecasting the hazard and risk associated with rock avalanches. The Orin3D model, based on the equivalent fluid concept, can be used to simulate rock avalanche motion, however it is unknown what the best model parameterization is for forecasting. However, Orin3D is implemented to run on a graphical processing unit (GPU), which improves simulation times by two orders of magnitude, making large-scale calibration feasible, as is investigated herein.
In the present work, we use a posterior analysis based on Bayesian statistics to calibrate Orin3D for three different parameterizations: 1) Frictional rheology, 2) Voellmy rheology, and 3) the combination of Frictional and Voellmy rheology, using a data set containing 22 historical rock avalanche cases, and requiring over 450,000 model runs. Based on the calibration results, a probabilistic prediction framework is then tested that generates pseudo-predictions for the cases in the database, incorporating key features of rock avalanches, such as path materials and topographic constraints. We find that, among these three rheological settings, the best prediction results for most cases are obtained with the combination of Frictional and Voellmy rheology. We further use these results to suggest a prediction procedure that considers the volume, path material and topographic confinements of rock avalanches, which provide guidance for the rheological setting in the model and important basis for the prediction and mitigation of rock avalanche hazards in practice.
How to cite: Deng, G. and Aaron, J.: Calibration and prediction procedure of rock avalanche through advancing numerical simulation, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-8845, https://doi.org/10.5194/egusphere-egu25-8845, 2025.
The field of landslide susceptibility modelling has seen the adoption of many different data-driven approaches, spanning from linear models to the most recent deep-learning solutions. In short, simpler models offer greater interpretability, while predictions derived from complex architectures are more difficult to explain. For this reason, complex algorithms are often referred to as black-box models. However, in the context of landslide susceptibility mapping, the ability to provide highly accurate results along with interpretable predictions is highly valuable. In light of these considerations, this study presents a landslide susceptibility mapping by exploring the capabilities of a new generation of interpretable models, namely Explainable Boosting Machines (EBMs). Unlike the majority of explainable approaches that unveil the decisions of a complex model in a post-processing phase, EBMs offer direct interpretability and full transparency. As a consequence, EBMs fall into the category of glass-box models. Notably, the incorporation of these models within studies focusing on the relationship between landslide occurrence and extreme rainfall events raises considerable interest and represents the aim of this work. Therefore, this contribution focuses on landslides triggered by a heavy rainfall event on September 15, 2022, in Central Italy. To analyze the interaction between landslide occurrence and the event, a novel rainfall variable is introduced among the set of predictors, capturing the event’s intensity relative to historical rainfall patterns. Specifically, this rainfall variable is computed as the percentage of precipitation attributed to the event compared to the mean annual rainfall. The rainfall variable also introduces a dynamic component to the proposed modelling, since it may vary at every future rainfall event. As a consequence, by combining the dynamic nature of the rainfall variable with the exact intelligibility of EBMs, the study also presents a landslide susceptibility mapping under potentially different rainfall scenarios with respect to the September 15, 2022 event.
How to cite: Caleca, F., Confuorto, P., Raspini, F., Segoni, S., Tofani, V., Casagli, N., and Moretti, S.: Modelling landslide susceptibility through a glass-box machine learning , EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-5730, https://doi.org/10.5194/egusphere-egu25-5730, 2025.
Shallow landslides are frequently observed at natural slopes and often lead to more destruction through flow-like disasters. Traditionally, physically-based landslide susceptibility models utilized infinite slope stability analysis to determine slope stability in terms of factor of safety (FS) over regional scales. Although the infinite slope model is computationally less demanding, it cannot account for the spatial variability of soil properties and the three-dimensional (3D) effects arising from complex topography. However, using 3D slope stability models is computationally demanding and suffers from discontinuity introduced by abrupt changes in soil thickness. Therefore, this research proposes a new Three-Dimensional Translational Shallow (3DTS) slope stability model to overcome these drawbacks of the existing models with complex 3D sliding surfaces.
The developed 3DTS model utilizes the Green-Ampt (GA) infiltration model and the 3D extension of the Janbu simplified method of slope stability. The 3DTS utilizes a generalized GA model to account for non-uniform infiltration history and compute the surface runoff. In 3D limit equilibrium slope models, the failing soil mass must be subdivided into rigid soil columns; however, the developed 3DTS uses the cells from a digital elevation model (DEM) as the rigid soil columns. The shear strength, modeled with the Mohr-Coulomb criterion, is provided by the soil frictional resistance on the base and the side regions of the outermost soil columns. Additional strength from the vegetation roots at the shallow surfaces is modeled.
The method used in the developed 3DTS model for generating slip surfaces from DEM cells was verified by comparing computed FS with the 3-Dimensional Probabilistic Landslide Susceptibility (3DPLS) model, which uses ellipsoidal slip surfaces. A parametric study analyzed the sensitivity of the slip surface's shape, the side soil resistance, and the vegetation resistance to shallow translational failures. The applicability and computational efficiency of the developed 3DTS for large-scale landslide susceptibility assessment were demonstrated by analyzing landslide case studies in Norway and South Korea.
This work is the result of collaboration between Norwegian Geotechnical Institute (NGI) and Korea Advanced Institute of Science and Technology (KAIST) through the project GEOMME (2021-2026; Pnr. 322469), “Climate-induced geohazards mitigation, management, and education in Japan, South Korea, and Norway”, supported by the Research Council of Norway.
How to cite: Cheon, E., Ahmet Oguz, E., DiBiagio, A., Piciullo, L., Kwon, T. H., and Lee, S. R.: A Physically-based 3D Landslide Susceptibility Model for Shallow Translational Landslides using DEM, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-17059, https://doi.org/10.5194/egusphere-egu25-17059, 2025.
Landslides pose significant threats to human lives and economies, with their frequency and intensity increasingly exacerbated by climate change. This was demonstrated in May 2023 in the Emilia-Romagna region (northern Italy), where 80,000 landslides were triggered by two rare, high-intensity rainfall events with a return period of 300 years, occurring just 14 days apart. The landslides exhibited diverse types and materials, necessitating tailored risk management approaches due to differences in volume, velocity, and post-event behaviour. To address these complexities, susceptibility maps must integrate both static predictors and dynamic triggering factors to better understand the relationships between rainfall and landslide types.
Using a detailed landslide inventory developed through collaboration between the Emilia-Romagna Geological Service and the universities of Modena and Bologna, we analysed the relationship between rainfall and five distinct mapped landslide types: debris slide, debris flow, earth slide, earth flow, and rock slide. This study introduces a Transformer Neural Network (TNN) to integrate static predictors (e.g., slope, aspect, geology, land cover) with dynamic rainfall data from the 30 days preceding the second rainfall event (16th May), capturing the influence of antecedent wet/dry conditions. The TNN processes rainfall time series data similarly to speech recognition algorithms, allowing it to model temporal dependencies effectively.
We evaluated the TNN with rainfall data at different temporal resolutions (daily and hourly intervals) and compared its performance against models using only static predictors or cumulative rainfall. The TNN was trained on 70% of the dataset, targeting specific landslide types to generate susceptibility maps tailored for each type. Model performance was assessed using a comprehensive set of metrics, including Area Under the Curve (AUC), Accuracy, Recall, F1 and F2 scores, Matthew’s Correlation Coefficient, and Kappa Coefficient. Additionally, we applied the SHAP (Gradient Explained) method to analyse the influence of rainfall on susceptibility values, revealing the model's internal decision-making processes.
The results demonstrate that integrating rainfall time series significantly enhances susceptibility mapping accuracy. The TNN using daily rainfall data produced the most reliable maps for all landslide types, except debris flows, where hourly intervals yielded slightly better results. SHAP analysis further illuminated the role of rainfall in susceptibility variations, providing valuable insights into the TNN's functionality. Overall, the TNN outperformed models using only static predictors or cumulative rainfall, offering a robust framework for understanding and predicting landslide susceptibility in diverse scenarios.
How to cite: Rani, R., Dahal, A., Lombardo, L., Tanyas, H., and Berti, M.: Unveiling the Complex Relationship Between Rainfall and Landslide Types Using a Transformer Neural Network, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-6943, https://doi.org/10.5194/egusphere-egu25-6943, 2025.
Landslide susceptibility is traditionally determined by analyzing various topographic, geological, and hydrological factors, which influence the probability of landslide occurrence. Recent research in Italy and Nepal has shown that landslide susceptibility is also controlled by landslide path dependency (LPD), where previous landslides locally and temporally influence the future landslide susceptibility. Our study focusses on Collazzone (Italy), a region predominantly affected by shallow landslides, supported by multi-temporal landslide inventory from 1939 to 2014. Here, we are comparing the impact of earlier landslides on landslide susceptibility in the downslope and lateral directions. We hypothesize that the LPD has more impact in downslope direction than in lateral direction due to the crucial role played by formation of positive feedback loop of the soil-landslide system. In the downslope direction, landslides can create weakly permeable soil layers that increase the water saturation, thus increasing the probability of subsequent landslides. In contrast, the lateral direction lacks this feedback mechanism, making subsequent landslides less likely.
For testing our hypothesis, we used simulated annealing to make artificial landslide inventories which approximate the real landslide inventory in terms of topographic positioning, but that lack any LPD. After that we calculate Ripley’s K by using a space-time cuboid for these control inventories and for the real inventory. Generalized additive models (GAM) were used to analyze the ratio between real and control Ripley’s K values. GAM results indicate that there is a nonlinear relationship between ratio of real to control Ripley’s K and time difference (dT), lateral (dL) and downslope distance (dD) between consecutive landslides. This ratio has a negative relationship with dL, and dD, while the relationship with dT is weak. Moreover, we found that earlier landslides have a stronger impact on future occurrence of landslides in downslope direction than in lateral direction. Our results provide clear evidence that downslope direction plays a significant role in landslide path dependency.
How to cite: Sodhi, H. S., Temme, A., Samia, J., Rossi, M., and Ardizzone, F.: Comparing the strength of Landslide Path Dependency in the downslope and lateral directions using simulated annealing and space-time clustering , EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-4137, https://doi.org/10.5194/egusphere-egu25-4137, 2025.
Landslide susceptibility assessment at scales wider than the single slope has been so far carried out mainly through heuristical/geomorphological and/or statistical methods, except for applications limited to shallow landslide predictions by means of infinite-slope limit equilibrium models (Godt et al. 2008). Owing to the complexity of developing quantitative deterministic susceptibility models at wide scales, taking into account also deep and complex landslide mechanisms, the application of limit equilibrium methods as well as numerical stress-strain methods have been historically limited to the scale of the single slope. However, the increased availability of powerful computational tools as well as the existence of detailed geological and geotechnical databases at scales that are intermediate between the single slope and the regional scales, as for example the scale of a single urban centre, allow for extending the application of three-dimensional limit equilibrium analysis to the assessment of landslide susceptibility at such scale, also taking into account the failure susceptibility of deep and complex landslide mechanisms. This contribution presents a physically-based methodology aimed at assessing landslide susceptibility at the urban area scale, for both shallow and deep instability processes involving urbanized areas that are diffusely affected by landsliding processes (Ugenti et al. 2025). The proposed methodology has been applied to the municipality of Carlantino (Daunia Apennines, Southern Italy) as a test case study, using the available geological and geomorphological datasets as well as the soil geotechnical property data. Based on a three-dimensional geotechnical model, 2.5 km2 wide, a three-dimensional limit equilibrium analysis has been develop to obtain a mechanically-based map of the safety factors at the urban area scale, assuming different scenarios related to the groundwater table depth, which has been validated against geomorphological evidence and remote sensing data. The proposed approach, which is supposed to be exportable to other geological environments, provides a valuable tool for quantitative assessment of the slope stability conditions of an overall urban area to be used for a more rational approach of urban planning policies and risk management activities.
References:
Godt J.W., Baum R.L., Savage W.Z., Salciarini D., Schulz W.Z., Harp E.L. (2008). Transient deterministic shallow landslide modeling: requirements for susceptibility and hazard assessments in a GIS framework. Engineering Geology, 102 (3-4), 214-226.
Ugenti A., Mevoli F.A., de Lucia D., Lollino P., Fazio N.L. (2025). Moving beyond single slope quantitative analysis: a 3D slope stability assessment at urban scale. Engineering Geology, 344, 107841, doi: 10.1016/j.enggeo.2024.107841.
How to cite: Lollino, P., Ugenti, A., Mevoli, F. A., de Lucia, D., and Fazio, N. L.: 3D limit equilibrium analysis: an opportunity for quantitative landslide susceptibility assessment at the scale of the urban area, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-16171, https://doi.org/10.5194/egusphere-egu25-16171, 2025.
The increase in localized heavy rainfall and intense storms due to climate change has led to a continuous rise in landslide damages in South Korea, including slope failures and debris flows. While post-landslide recovery and damage site assessments are crucial, it is equally important to develop proactive and systematic landslide adaptation strategies to predict and prepare for landslides in advance. This study aims to develop an interpretable machine learning-based landslide susceptibility model and analyze landslide-prone areas under future climate change scenarios. Through this approach, it seeks to clearly identify the impact of forest management factors on landslides and establish effective adaptation strategies tailored to climate change scenarios. A dataset comprising 6,517 recorded landslide events from 2011 to 2024 was utilized. Various external and internal conditioning factors were obtained and constructed with a resolution of 100 meters. Climate scenario analysis employed SSP 1-2.6, SSP 2-4.5, and SSP 5-8.5, with extreme climate factors including rainfall intensity, the number of heavy rain days, daily rainfall, and 5day cumulative rainfall. Notably, changes in stand age class, DBH class, and growing stock under future forest management scenarios were calculated and integrated into the landslide model, enabling an evaluation of how management strategies affect landslide susceptibility. Results were validated by comparing past actual occurrence data. The SSP 5-8.5 scenario indicates a significant increase in landslide occurrences. These findings provide valuable insights into the effects of climate change on landslide susceptibility in South Korea and examine the potential of future landslide management strategies to mitigate associated susceptibility.
How to cite: Kim, U., Lee, S., Roh, M., Kim, S., and Lee, W.: Spatial Prediction of the Future Landslide Susceptibility under the SSP Scenario Using Machine learning Algorithms, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-14735, https://doi.org/10.5194/egusphere-egu25-14735, 2025.
An active-layer detachment slide (ALDS) occurred on September 21, 2018, in the Fenghuoshan mountains of the Qinghai-Tibet Plateau (QTP) (34◦39.1′N, 92◦53.5′E). With the Sentinel-1A image from Copernicus Open Access Hub, we use small baseline subset to achieve the time series deformation map to analyze the thermo-spatial creep feature, motion pattern, trigger mechanism, and correlation of environmental changes in the ALDS. The SBAS (the Small Baselines Subset) results show that the trailing part of ALDS has the largest downward deformation rate; however, the leading area was small, and the creep feature shows a clear seasonal change corresponding to the freeze-thaw cycle. We also divide the motion pattern into three stages: moderate creep, steady creep, and rapid collapse, based on the deformation rate. Meteorological observation and reanalysis data, as well as borehole data, show that heavy precipitation in the summer of 2017 and 2018 promote the formation of underground ice, while high air temperatures allow the thaw plane to reach the ice-rich zone, and confined water generated by the two-way freezing process result in ALDS. Moreover, there exists a temporal delay of approximately one month in the association between deformation rate and both precipitation and temperature. Furthermore, there is a clear correlation between variations in thawing depth and deformation, which serves as the primary catalyst for ALDS in permafrost regions. Finally, we also identify that ALDS is a mixed-type landslide and that cumulative deformation and creep damage play the main roles in triggering ALDS.
How to cite: wen, Z. and wang, F.: Creep features and mechanism of active-layer detachment slide on the Qinghai-Tibet Plateau by InSAR, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-2110, https://doi.org/10.5194/egusphere-egu25-2110, 2025.
Posters on site: Mon, 28 Apr, 10:45–12:30 | Hall X3
From September 26 to 28, 2024, Nepal experienced exceptionally heavy rainfall, severely impacting large areas, particularly the Kathmandu Valley and its surrounding districts, triggering flash floods and landslides. This study presents preliminary findings from an assessment conducted approximately two months after the event, focusing on the upstream region of the Nallu Khola watershed in Lalitpur District, one of the areas most severely impacted. The event recorded a cumulative rainfall total of 518 mm at the Lele AWS Station (Department of Hydrology and Meteorology, Nepal), located approximately 2 km NW of the study area. This rainfall was about 4.3 times the total monthly rainfall for September 2023. The maximum hourly rainfall, observed at 5:00 AM on September 28, reached 39.8 mm, while the highest 24-hour rainfall was an extraordinary 441.2 mm. The rainfall triggered a series of compound sediment disasters, including raising the river level by approximately 3 m above the riverbed, along with numerous landslides and debris flows. The landslides predominantly consist of shallow failures, primarily occurring along roads and in areas associated with cultivated land, while areas covered with forest exhibit relatively few failures. Debris flows are predominantly concentrated in creeks, with a comparable event having occurred on September 30, 1981. Following that event, debris flow mitigation engineering measures (e.g. gabion check dams and channel works) were implemented in some creeks and the impacts of the 2024 event appear to have been largely confined to these mitigated creeks. This underscores the importance of implementing and maintaining effective mitigation measures to manage debris flow hazards in vulnerable areas.
How to cite: Tsou, C.-Y., Bhusal, Z., Kakinuma, H., Kawakami, R., Higaki, D., Bhusal, J. K., Dhakal, S., and Amatya, S. C.: Sediment disasters induced by the 26-28 September 2024 extreme rainfall event in Nallu Khola watershed of Central Nepal, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-2400, https://doi.org/10.5194/egusphere-egu25-2400, 2025.
Landsliding problems in slopes surrounding reservoir lakes are sometimes induced or reactivated by reservoir operation activities (Xia et al., 2015). Regular landslide susceptibility assessments are essential for safeguarding lives, infrastructure, and the environment, this being highly amplified for reservoir basins. Landslide assessments are commonly done through heuristic, statistical and physically-based quantitative methods such as limit equilibrium (LE) analysis. However, quantitative LE analyses have been historically carried out in 2D and at single-slope scale due to the need of reducing computational requirements, although realistically slope failures are 3D in nature; hence, using 3D methods can likely yield more accurate results and is more suitable for the understanding of the landsliding processes. Nowadays, with increased computational capability, it is possible to move to more representative 3D approaches and even attempt to extend the scale of application not only for shallow landslides, but also deep and complex landsliding processes. Since most 3D LE analyses are performed at slope scale, this study aimed at moving from slope to reservoir basin scale to assess the overall susceptibility to slope failure at the San Pietro Dam. The adopted methodology used Slide3 Software and involved generation of study area's 3D geometry from a 10-m resolution DEM. Then, stratigraphic borehole data, along with stratigraphic sections obtained from geological reports for the area were used to reconstruct 3D geological schematization. Geotechnical strength parameters between residual and peak strength derived from literature were used as inputs for stability analysis. Specifically, 3D extension of the Morgenstern and Price method, which divides the potential failure surface into columns based on Cheng & Yip (2007) formulation for asymmetrical slopes was used. Results indicate that the approach is able to provide distribution of potential areas susceptible to slope instability as safety factor (SF) values which were in good agreement with field observations and the landslide inventory map. In particular, many landslides fall in marginally stable pixels of the SF map and can reactivate depending on the increase of water table levels along the slopes. Effect of potential rapid drawdown of the reservoir level on the stability of surrounding slopes was also investigated. The results shed light on possible extension of 3D LEM to scales larger than a slope so that it can become a useful tool for landslide risk management in reservoir environments.
References
Cheng, Y. M., & Yip, C. J. (2007). Three-Dimensional Asymmetrical Slope Stability Analysis Extension of Bishop’s, Janbu’s, and Morgenstern–Price’s Techniques. Journal of Geotechnical and Geoenvironmental Engineering, 133(12), 1544–1555. https://doi.org/10.1061/(asce)1090-0241(2007)133:12(1544)
Xia, M., Ren, G. M., Zhu, S. S., & Ma, X. L. (2015). Relationship between landslide stability and reservoir water level variation. Bulletin of Engineering Geology and the Environment, 74(3), 909–917. https://doi.org/10.1007/s10064-014-0654-0
How to cite: Chikalamo, E., Lollino, P., and Mavroulli, O.: Reservoir-Scale Landslide Susceptibility Analysis of Slopes Surrounding Artificial Impoundments by Three-Dimensional (3D) Limit Equilibrium Models: A Case Study of San Pietro Dam, Avellino Province, Italy., EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-2781, https://doi.org/10.5194/egusphere-egu25-2781, 2025.
The Ixhuatlan de Madero area is located (in the geology of Mexico) between the Gulf of Mexico and the micro-continent Oaxaquia. The regional stratigraphy comprises the Paleocene Chicontepec Formation (chiefly sandstones and shales), overlying the Cretaceous Mendez and Tamaulipas formations, respectively composed of shales and limestones. Analysis of the structural data collected in the field indicates five stages of deformation. The first stage is characterized by upright folds plunging to the NE and SW. The second stage corresponds to the Laramide orogeny (i.e. ~ 40 Ma) and involves NE-vergent folds. The folding produced south-westwards shallow-dipping layers (i.e. less than 30°) and overturning of the first stage folds to the NE. The third stage is marked by reverse faults compatible with NE-SW compression as observed in the village of Cantollano. In contrast, NW-SE normal faults observed to the NW of Ixhuatlan reveal a fourth stage characterised by an extensional regime. The fifth stage involves NW-SE and NE-SW fractures present mainly west of Ixhuatlan de Madero. The latter fractures represent pronounced weakness zones within the rock mass and are further opened by plant roots and excavated by the tropical rains of the region. The control local disintegration of the rock and lead eventually to landslides. The landslides promote mass transport towards the NE and SW dominantly and, furthermore, the building of houses and human infrastructures amplify them.
How to cite: Pascal, C., Zarate Velazquez, O. D., and Lopez Doncel, R. A.: Geological Controls on natural hazards in Ixhuatlan de Madero, Veracruz, Mexico, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-6333, https://doi.org/10.5194/egusphere-egu25-6333, 2025.
The mapping unit is a classification of land surfaces based on specific criteria, serving as the fundamental unit for spatial data extraction in landslide susceptibility analysis. In the landslide susceptibility analyses, a grid unit is frequently employed due to its ease of generation as uniform grid cells of a designated size. However, the utilization of grid units does have certain limitations. Specifically, these units often fail to accurately represent the actual site topography. Consequently, they result in an incomplete consideration of valley and drainage lines, which are critical factors in landslide occurrence. In contrast, slope units, delineated based on hydrological criteria (e.g., ridges, valleys), offer a more topographically accurate representation. This is due to the fact that they integrate spatial data and topographic factors more effectively into the analysis than grid units.
This study aims to compare the impact and performance of grid units and slope units in landslide susceptibility analysis. To this end, the study utilizes various analytical techniques to evaluate the influence of conditioning factors across these mapping units. The study area, designated as Jecheon-si, Chungcheongbuk-do, Republic of Korea, was selected to assess the impact of mapping units due to its experience with several landslides in August 2020. The analysis incorporated a range of conditioning factors, including elevation, slope aspect, slope angle, standard curvature, planar curvature, profile curvature, Specific Catchment Area (SCA), Topographic Wetness Index (TWI), Stream Power Index (SPI), forest type, forest density, forest stand height, timber diameter, timber age, soil texture, soil depth, slope shape, topography, land use, and lithology. In order to assess the significance and contribution of these factors, visualization techniques were employed, including SHAP (Shapley Additive Explanations) plots, summary plots, and dependence plots. These methods facilitated a comparative analysis of factor importance and influence on landslide susceptibility using the two mapping units. Additionally, correlation analysis among the conditioning factors and trend identification within each unit were conducted to enhance the accuracy and interpretability of the results. The landslide susceptibility analysis was implemented using a Multi-Layer Perceptron model, and the performance of the model was evaluated using the Area Under the Curve (AUC). Finally, the results of the study were analyzed to compare and evaluate the relative advantages and limitations of the slope unit and the grid unit in landslide susceptibility assessment.
Acknowledgments
This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (RS-2023-00222563) and by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (RS-2024-00463587).
How to cite: Jo, E.-B., Lee, J.-H., and Park, H.-J.: Analysis of the Impact of Mapping Units on Landslide Susceptibility: A Comparative Study of Grid Units and Slope Units, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-7844, https://doi.org/10.5194/egusphere-egu25-7844, 2025.
Landslides, a major natural disaster in South Korea, are primarily triggered by heavy rainfall associated with global climate anomalies. In particular, the years 2020 and 2022 witnessed unprecedented torrential rains during the summer, resulting in the most severe landslide damages recorded in recent history, with significant human and economic losses.
Landslide susceptibility assessment involves the spatial analysis of direct triggering factors, such as rainfall, and conditioning factors both internal and external to slopes, to predict the likelihood and impact of landslide occurrences. Based on the mechanisms considered, assessment methodologies are typically classified into physically-based models and data-driven models. Physically-based models, which have been extensively studied globally, are well-suited for landslide susceptibility analysis in South Korea as they allow for the integration of engineering principles to address rainfall and internal slope conditions. However, their limitations in addressing the multifaceted interactions among diverse influencing factors necessitate the incorporation of data-driven approaches.
This study seeks to integrate physically-based models with data-driven models to capture both the engineering mechanisms of rainfall-induced landslides and the complex interrelationships among diverse influencing factors. Since these models operate as independent frameworks, a fusion approach is adopted to combine their outputs effectively. Fusion methodologies vary depending on the stage at which data or information is integrated. In this research, decision-level fusion is employed, which aggregates the independent decisions or outputs of multiple models to produce the final result. Specifically, serial decision fusion and parallel decision fusion, two representative decision-level fusion techniques, are implemented. The study evaluates the performance and applicability of the fusion models by comparing the outcomes of different fusion strategies.
Acknowledgements
This work was supported by the National Research Foundation of Korea(NRF) grant funded by the Korea government(MSIT) (No. RS-2024-00358026).
How to cite: Lee, J.-H., Park, H.-J., and Lee, Y.-J.: A Study on Landslide Susceptibility Fusion Models Using Decision-Level Fusion, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-7971, https://doi.org/10.5194/egusphere-egu25-7971, 2025.
Typhoon-induced shallow landslides have caused significant economic losses and casualties in China's coastal regions. Accurate prediction and hazard assessment of typhoon-induced landslides are crucial for effective geohazard prevention and management. However, providing accurate hazard evaluation remains challenging due to limited data on rainfall triggers and relevant geological parameters. Therefore, our study integrates the effective rainfall model and the probabilistic physical model TRIGRS to analyze the early warning of regional shallow landslides. In this study, we selected Daoshi Town, in Zhejiang Province of China, which was heavily impacted by Super Typhoon Lekima on August 10, 2019. To find out the distribution and regularity of landslides after typhoon rainfall, we identified a total of 190 shallow landslides through field surveys and remote sensing interpretation. The soil thickness of the study area was simulated using the random forest algorithm based on the soil thickness dataset from the field survey. Rainfall characteristics and thresholds were established using an effective rainfall model that accounts for the 6-hour rainfall on the day of analysis and the cumulative rainfall over the preceding three days. To assess slope stability under different rainfall scenario, TRIGRS was employed, considering key parameters of different soil types such as cohesion and internal friction angle. The results indicate that 90% of the landslides occurred in areas classified as highly unstable. Validation using landslide data from 2020 and 2021 showed that 81% of new landslides occurred in unstable areas, demonstrating the reliability of the proposed early warning approach. It shows that our results are reliable and can provide reference for the hazard assessment and management of typhoon-induced shallow landslides in coastal regio.
How to cite: Liao, Y., Chen, L., Li, Y., and Yin, K.: Hazard assessment for typhoon-induced shallow landslides based on rainfall thresholds and physical modeling, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-8942, https://doi.org/10.5194/egusphere-egu25-8942, 2025.
Shallow landslides, typically occurring on steep slopes, are often triggered by intense, short-duration rainfall or extended periods of lighter rainfall. These events present severe hazards in mountainous regions, causing substantial soil loss, fatalities, and economic damage (Tonini and Cama, 2019). Accurate prediction and early warning systems are essential for mitigating such impacts. To address these challenges, previous studies in Switzerland have examined rainfall thresholds related to landslide triggering by regionalizing landslide occurrences according to geomorphological factors (Leonarduzzi et al., 2017). To enhance the overall accuracy of such predictions, it is essential to utilize datasets with higher temporal and spatial resolution.
This work adapts a robust deep learning approach initially developed by Mondini et al. (2023) for Italy to the case of Switzerland. Unlike previous studies that relied solely on rain gauge data, which is often highly variable, we use the CombiPrecip product from the Swiss Federal Office of Meteorology. This product integrates radar measurements with rain gauge data to provide a kilometer-scale, hourly precipitation dataset covering the past 20 years. The landslide input dataset comes from the Swiss Federal Institute for Forest, Snow, and Landscape Research (WSL), which has systematically collected data on damage caused by naturally triggered floods, debris flows, and landslides since 1972 (Hilker et al., 2009).
To compensate for the relative sparsity of landslide events in our training set, we carry out an ensemble approach where we train 24 classifiers, thus resulting in increased robustness and a probabilistic outcome. The ultimate goal of this research is to compare various classification algorithms and evaluate their integration into an early warning system that leverages susceptibility maps and geological factors.
REFERENCES
- Tonini M, Cama M (2019). Spatio-temporal pattern distribution of landslides causing damage in Switzerland. Landslides 16, 2103–2113.
- Leonarduzzi E, Molnar P, McArdell BW (2017). Predictive performance of rainfall thresholds for shallow landslides in Switzerland from gridded daily data. Water Resources Research, 53(8), 6612‑6625.
- Mondini AC, Guzzetti F, Melillo M (2023). Deep learning forecast of rainfall-induced shallow landslides. Nature Communications, 14(1), 2466.
- Hilker N, Badoux A, Hegg C (2009). The Swiss flood and landslide damage database 1972-2007. Nat Hazards Earth Syst Sci 9:913–925.
How to cite: Soutter, J., Dunand, M., and Tonini, M.: Enhancing rainfall-triggered landslide forecasting in Switzerland using ensemble learning, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-9871, https://doi.org/10.5194/egusphere-egu25-9871, 2025.
Shallow landslides are characterized by a superficial sliding surface whose depth is at least five meters below the ground. Their occurrence has increased in recent decades due to climate change, especially in Northern Italy where extreme meteorological events (the main triggering factors) have been reported to increase in intensity. Since shallow landslides are a very common geohazard in mountain and hilly areas, whose consequences can be catastrophic both for people and the natural environment, new methodologies that aim to better estimate landslide susceptibility have been explored in the literature. Here a new tool called BEMSL (“Basin Ensemble Models for Shallow Landslides”) has been developed to forecast effectively shallow landslides at the basin scale.
The BEMSL is a model that considers an ensemble approach for susceptibility mapping, and it is conceptually divided into three parts. Primarily, it includes different limit equilibrium and infinite slope formulations that describe the stability of a slope in terms of safety factor (FS), which is defined as the ratio between stabilizing and destabilizing actions. Even if, theoretically, the actions acting on a slope should be always the same, many authors in this field have proposed different FS equations, trying to choose the most relevant acting actions depending on the local geology, soil composition and other predisposing factors. Consequently, it is difficult to choose the most suitable FS formulation that fits best to the considered situation. To provide a unique answer, the second part of BEMSL includes the Random Forest (RF) approach that creates a model ensemble able to merge the outputs from the implemented FS formulations. Since RF is a machine-learning algorithm that works autonomously on FS data provided, countermeasures to avoid overfitting have been considered. In the last part, the output validation was assessed using the ROC (“Receiver Operation Characteristics”) approach, which essentially consists of the quantification of how many true positives (TP), false positives (FP), true negatives (TN) and false negatives (FN) compared to the available landslide census.
BEMSL was applied to retrieve dynamic landslide susceptibility maps, giving site-specific insight into the probability of shallow terrain failures. The reliability of this BEMSL tool was tested considering the event that happened in July 1987 in Tartano Valley (Sondrio province, located in Northern Italy). In the late afternoon of 18 July 1987, an extreme storm triggered several shallow landslides across Tartano Valley, which evolved into a catastrophic debris flow, resulting in 21 casualties and extensive infrastructure damages. In this case study, the risk of failure of punctual and linear electrical powerlines was investigated using the BEMSL. A dependence on the risk of failure due to the rainfall intensity temporal evolution has shown the vulnerabilities of the Tartano Valley electrical infrastructures developed during the extreme geo-hydrological event.
How to cite: Abbate, A., Scaioli, A., Corti, M., Papini, M., and Longoni, L.: A new tool for studying shallow landslides at the basin scale: BEMSL, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-10083, https://doi.org/10.5194/egusphere-egu25-10083, 2025.
How to cite: Fulmer, E., Leshchinsky, B., Senogles, A., Olsen, M., Mohney, C., and Glover-Cutter, K.: Improving Prediction, Response, and Safety through Real-Time Surface Monitoring with RTK-GNSS Arrays, Oregon, USA, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-13930, https://doi.org/10.5194/egusphere-egu25-13930, 2025.
Rain-induced landslides pose a global threat, resulting in significant casualties and infrastructure damage annually. Such impacts can be reduced utilizing efficient early warning systems to plan mitigation measures and protect vulnerable elements. This study presents an innovative geophysical monitoring approach that combines electrical resistivity tomography (ERT) and quasi-distributed opto-electronic sensing (OES), deployed on a clay rich slope typical of thousands in the Greater Bay Area, China. ERT is used to generate detailed dynamic resistivity maps, combined with OES-indicated moisture content, highlighting the spatial-temporal distribution of slope-scale moisture content. The relationship between the analytical solution of Factor of safety informed by ERT-derived dynamic moisture maps and contemporaneous landslide displacement is confirmed by quasi-distributed OES strain measurements. By revealing the connection between landslide movement and ERT-OES-informed slope stability, this combined ERT and OES monitoring approach offers new insights into landslide mechanisms. Our study demonstrates the importance of relying on multi-source observations to develop effective landslide risk management strategies and accents the advantages of incorporating subsurface geophysical monitoring techniques to enhance landslide early warning approaches.
How to cite: Shen, P.: In-situ Time-lapse Geophysical Monitoring for Rain-induced Landslide Forecasting, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-16837, https://doi.org/10.5194/egusphere-egu25-16837, 2025.
Due to the impact of climate change, the increasing frequency of extreme rainfall events, with concentrated rainfalls, commonly cause landslide hazard in the mountain areas of Taiwan. However, there are uncertainties for the predicted rainfall as well as the landslide susceptibility analysis. This study employs machine learning approached, including the logistic regression method LR to analyze the landslide susceptibilities. Together with the predicted temporal rainfall, the predictive analysis of landslide susceptibility was performed in the adopted study area in Central Taiwan. The uncertainties within the rainfall prediction was firstly investigated before applied to the landslide susceptibility analysis. To assess the susceptibility of the landslides, logistic regression method LR was applied. The results of predictive analysis, with the discussions on the accuracy and uncertainties, can be applied for the landslide hazard management.
How to cite: Shou, K.-J.: Spatial and Temporal Analysis of Landslide Susceptibility– for the Case in Taiwan, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-2117, https://doi.org/10.5194/egusphere-egu25-2117, 2025.
Landslide susceptibility assessment is crucial for preventing landslide risks. However, existing methods only consider local environmental features related to landslides, neglecting remote yet interconnected geographical features, leading to unreliable landslide susceptibility maps. This study fully considers the complex terrain and landform features of mountainous areas where landslides occur. From the perspectives of mapping units and susceptibility assessment models, it introduces geographical environmental correlations to achieve a comprehensive association between landslides and affected environments, thereby improving the accuracy of landslide susceptibility assessments. At the same time, since the world's first scientific satellite dedicated to serving the United Nations 2030 Agenda for Sustainable Development, the Sustainable Development Goals Scientific Satellite 1 (SDGSAT-1), was launched in 2021, its potential in monitoring and assessing landslide disasters remains to be developed. Therefore, this study innovatively applies SDGSAT-1 data in the field of landslide research and conducts landslide susceptibility assessment in Jiulong County, Ganzi, based on the optimal scale slope units and Graph Neural Networks (GNN).
We propose the following method: First, establish appropriately sized slope units using R.Slopeunits to simulate complex mountainous terrain. Second, extract various landslide influencing factors using SDGSAT-1 satellite imagery data. Then, select the most representative graph nodes by constraining environmental similarity and influencing factor feature similarity, constructing a graph structure. Finally, perform landslide susceptibility assessment in the study area using the GraphSage model, which includes environmental information aggregation.
This study's distinctive feature lies in fully considering the complex terrain and landform characteristics of mountainous areas where landslides occur. From the perspectives of mapping units and evaluation models, it introduces geographical environmental correlations to achieve a comprehensive association between landslides and affected environments. Furthermore, to validate the effectiveness of the proposed method, we selected raster units and the classic Artificial Neural Network (ANN) model as control experiments. Simultaneously, we conducted comparative experiments using Landsat and SDGSAT-1 satellite imagery, analyzing differences from two aspects: landslide influencing factors and landslide susceptibility evaluation results.
The results indicate that: (1) Compared to the commonly used Landsat series satellite data in previous studies, SDGSAT-1 satellite imagery offers higher spatial resolution, capturing more spectral information with richer hue and detail. Additionally, it can generate more angles of landslide influencing factors compared to Landsat satellite data. (2) Employing global heterogeneity evaluation metrics allows for reasonable determination of slope unit scales, thereby maximizing internal consistency and external heterogeneity control within slope units. (3) By utilizing the Graph Neural Network (GNN) model that incorporates environmental information aggregation for landslide susceptibility assessment in the study area, it can, to some extent, overcome spatial limitations and integrate complex mountainous environmental information, facilitating the induction of reliable landslide characteristics.
How to cite: Lei, X., Liu, H., Chen, Z., Li, S., Chen, H., Zeng, S., Wang, X., Bai, W., Li, W., and Picco, L.: Investigating the Landslide Susceptibility Assessment Methods for Multi-Scale Slope Units Based on SDGSAT-1 and Graph Neural Networks., EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-4345, https://doi.org/10.5194/egusphere-egu25-4345, 2025.
