10th International Conference on Geomorphology
© Author(s) 2022. This work is distributed under
the Creative Commons Attribution 4.0 License.

Applying a hierarchical Generalized Additive Model to integrate predisposing, preparatory and triggering factors for landslide prediction

Stefan Steger1, Mateo Moreno1,2, Alice Crespi1, Peter James Zellner1, Robin Kohrs1,3, Jason Goetz3, Stefano Luigi Gariano4, Maria Teresa Brunetti4, Massimo Melillo4, Silvia Peruccacci4, Lotte de Vugt5, Thomas Zieher6, Martin Rutzinger5, Volkmar Mair7, and Massimiliano Pittore1
Stefan Steger et al.
  • 1Eurac Research, Institute for Earth Observation, Bolzano-Bozen, Italy (stefan.steger@eurac.edu)
  • 2University of Twente, Faculty of Geo-Information Science and Earth Observation (ITC), Enschede, The Netherlands
  • 3Department of Geography, Friedrich Schiller University Jena, Germany
  • 4CNR IRPI, Perugia, Italy
  • 5Department of Geography, University of Innsbruck, Innrain 52f, 6020 Innsbruck, Austria
  • 6Institute for Interdisciplinary Mountain Research, Austrian Academy of Sciences, Innrain 25, 6020 Innsbruck, Austria
  • 7Office for Geology and Building Materials Testing, Autonomous Province of Bolzano-South Tyrol, Cardano, Italy

Shallow landslides are frequently caused by an interplay of static predisposing factors (e.g., topography), mid-term preparatory factors (e.g., prolonged rainfall, seasonal changes of vegetation, snow melt), and short-term triggers (e.g., heavy rainfall). For large-area assessments, statistical analyses and data-driven approaches are often used to model landslide susceptibility based on spatial environmental variables or to derive critical landslide-triggering rainfall conditions through the definition of empirical rainfall thresholds. Attempts to integrate the spatial and temporal domains in the context of quantitative regional-scale landslide prediction are still rare.

This contribution focuses on the landslide-prone area of South Tyrol, northern Italy (7,400 km²) and presents a novel data-driven modelling procedure that integrates spatial predisposing factors and dynamic preparatory and triggering factors to predict the probability of landslide occurrence in space and time. The approach is based on time-stamped landslide inventory data from 2000 to 2021, high-resolution gridded daily precipitation observations for the same period, and a set of relevant static environmental variables (e.g., including topographic indices, lithology). Data preparation included an initial filtering of rainfall-induced landslide presence observations and a rule-based stratified sampling of landslide absence observations at landslide locations and at non-landslide locations. Cross-validation was implemented in the model developing stage to select optimal time windows to represent pre-landslide preparatory and triggering cumulative rainfall conditions. Modelling was based on a binomial Hierarchical Generalized Additive Model (HGAM) that considered the non-linear influence and interactions (i.e., via tensor products) of static and dynamic environmental variables on landslide occurrence (presence vs. absence) while simultaneously accounting for the nested data structure (i.e., multiple considerations of each location) and seasonal effects. The study also considered different biases inherent in the input data, such as the known underrepresentation of landslide data in locations far from infrastructure or potential reporting biases across years, by averaging-out associated random effect variables. The results were validated quantitatively (spatial and temporal cross-validation) and qualitatively (geomorphic plausibility) and visualized as maps and probability surface plots.

The assessment and validation confirmed the high generalizability and predictive performance of the model. A closer look at the derived relationships allowed to uncover the effects of predisposing, preparatory and triggering factors on landslide occurrence as well as the associated season-dependent variations. From our perspective, this new approach represents a good compromise between the high model flexibility for landslide prediction purposes and the high interpretability for understanding the underlying relationships.

The research leading to these results is related to the PROSLIDE project that received funding from the research program Research Südtirol/Alto Adige 2019 of the Autonomous Province of Bozen/Bolzano – Südtirol/Alto Adige.

How to cite: Steger, S., Moreno, M., Crespi, A., Zellner, P. J., Kohrs, R., Goetz, J., Gariano, S. L., Brunetti, M. T., Melillo, M., Peruccacci, S., de Vugt, L., Zieher, T., Rutzinger, M., Mair, V., and Pittore, M.: Applying a hierarchical Generalized Additive Model to integrate predisposing, preparatory and triggering factors for landslide prediction, 10th International Conference on Geomorphology, Coimbra, Portugal, 12–16 Sep 2022, ICG2022-388, https://doi.org/10.5194/icg2022-388, 2022.