Shallow landslides, including debris slides, debris avalanches and debris flows, involve the failure of soft surficial material at rooting depth. They are overwhelmingly triggered by rainfall, which raises porewater pressure in the soil until the shear strength of the slope is overcome. Under climate change scenarios, heavy precipitation over mid-latitude and tropical mountainous regions is expected to intensify, making reliable predictions of where these failures will occur ever more urgent. Yet the field of landslide susceptibility mapping has a dirty secret: many recent studies report near-perfect prediction accuracies that the Norwegian team argues are almost certainly artefacts of overfitting rather than genuine predictive skill.

The root of the problem lies in space itself. Landslides tend to occur in clusters, and the environmental variables that influence them, from elevation to precipitation, vary smoothly across the landscape. When researchers randomly split their data into training and test sets, as standard machine learning practice dictates, points in the test set sit tantalisingly close to points in the training set, sharing nearly identical conditions. The model appears to predict brilliantly while merely memorising local patterns. The team quantified this effect with variogram analysis, showing that under random cross-validation, the residual spatial correlation of key predictors such as relative precipitation approached unity at the distances separating training and test observations, meaning a held-out point was essentially indistinguishable from its nearest training neighbour.

To confront this challenge head-on, the researchers built an exceptionally clean dataset. Rather than relying on Norway’s crowd-sourced landslide database, which misses many events and clusters along roads, they systematically mapped landslide scars using satellite-derived vegetation change indices, validated with high-resolution imagery, orthophotos and fieldwork. This effort detected 648 landslides, a 246 percent increase over the crowd-sourced record, and reduced spatial bias towards roads by 35 percent. After excluding debris floods, which follow a different release mechanism, the final inventory contained 571 landslide starting points, each characterised by twenty explanatory variables spanning topography, geology, hydrology and forest structure.

The team then ran thirty-two gradient boosted decision tree models, systematically varying three dimensions: whether hyperparameters were tuned with nested or simple cross-validation, whether the model included an explicit spatial random effect through the Gaussian Process Boosting framework, and which of four cross-validation schemes was used, ranging from fully random splits to spatially clustered folds and large and small spatial blocks. The results were sobering. Models with random cross-validation posted the highest performance metrics, with test AUC values reaching 0.90, but left significant spatial autocorrelation in their residuals, a statistical fingerprint of over-optimism. Only spatial models combined with block cross-validation fully eliminated residual autocorrelation, and they paid for this honesty with lower predictive power, achieving a test AUC of 0.77 at best.

Perhaps the most unsettling finding concerned variable importance. The top predictors identified by spatial models differed markedly from those of non-spatial models, and rankings shifted dramatically across cross-validation schemes. The researchers attribute this instability largely to spatial confounding, a phenomenon in which predictors correlate with unmodelled spatial structure and absorb its signal. In non-spatial models, variables with strong spatial footprints, such as total precipitation and bedrock weatherability, dominated, likely acting as proxies for spatial pattern rather than true causal factors. When an explicit spatial term was added, those variables lost much of their apparent importance. No single model configuration emerged as optimal, a candid admission that underscores how much work remains in this field.

Yet across the ensemble of models, a coherent physical story did emerge. Landslides were most likely on south-facing slopes, at elevations below roughly 480 metres, and where water inputs were greatest, whether from intense three-hour rainfall bursts, high flow accumulation, or elevated pre-storm groundwater. The team interprets elevation and aspect as indicators of soil conditions that existing datasets fail to capture directly: lower elevations and south-facing slopes in Norway’s climate experience more weathering, producing thicker and more stratified soils in which impermeable layers can channel water and build the critical pore pressures that trigger failure. Notably, bedrock weatherability, a novel reclassification of Norway’s national bedrock map into five classes, proved influential, with easily weathered shales and phyllites increasing landslide probability while resistant gneisses and granites reduced it, explaining why some drenched regions escaped unscathed.

The forest analysis delivered the study’s most surprising result. In models restricted to forested terrain, tree type consistently ranked among the top predictors, with deciduous stands showing markedly higher landslide susceptibility than spruce or pine forests. This runs counter to simple mechanical expectations, since spruce roots provide substantial reinforcement. The authors suggest the pattern reflects where deciduous forest, often birch-dominated, typically grows in Norway: steep, thin-soiled terrain near the tree line that is unsuitable for production forestry. The findings also revealed that forests suppress landslides on gentler slopes, with landslide probability in forest rising only above slope angles of about 30 degrees, compared with a much wider range in open land. Canopy cover above 75 percent further reduced susceptibility, supporting the role of forests as nature-based solutions for landslide mitigation.

The implications reach well beyond Norway. The team recommends that future machine learning studies of landslide susceptibility test for spatial autocorrelation in both response and predictor data, explore explicitly modelled spatial structure, employ multiple spatial cross-validation techniques, and report methods in full detail. Performance metrics and predictor rankings from studies that ignore spatial structure, they caution, should be interpreted with considerable scepticism, particularly when models are meant to generalise to new locations. As extreme rainfall events grow more frequent in a warming climate, the difference between an honest model and an overfitting one may determine whether hazard maps guide evacuation in time. The Hans storm, which forced the evacuation of around 4,600 people with no loss of life, showed what good early warning can achieve; this study shows how the science of predicting where the ground will fail can be made considerably more trustworthy.

Subject of Research: Spatial machine learning modelling of explanatory factors for rainfall-triggered shallow landslide release in southern Norway

Article Title: Spatial machine learning modelling reveals that soil indicators and tree type best explain shallow landslide release

Article References: Rüther, D. C., Haualand, K. F., Peeters, I. L. J., & Gillespie, M. A. K. (2026). Spatial machine learning modelling reveals that soil indicators and tree type best explain shallow landslide release. Natural Hazards and Earth System Sciences, 26(10), 4675-4721. https://doi.org/10.5194/nhess-26-4675-2026

Image Credits: AI Generated

DOI: 10.5194/nhess-26-4675-2026

Keywords: shallow landslides, machine learning, spatial cross-validation, landslide susceptibility, storm Hans, Norway, soil thickness, bedrock weatherability, forest cover, tree type, gradient boosting, spatial autocorrelation