The impact of terrain model source and resolution on snow avalanche modeling

Natural hazard models need accurate digital elevation models (DEMs) to simulate mass movements on real-world terrain. A variety of platforms (terrestrial, drones, aerial, satellite) and sensor technologies (photogrammetry, lidar, interferometric synthetic aperture radar) are used to generate DEMs at a range of spatial resolutions with varying accuracy. As the availability of high-resolution DEMs continues to increase and the cost to produce DEMs continues to fall, hazard modelers must often choose which DEM to use for their modeling. We use satellite photogrammetry and topographic lidar to generate high-resolution DEMs and test the sensitivity of the Rapid Mass Movement Simulation (RAMMS) software to the DEM source and spatial resolution when simulating a large and complex snow avalanche along Milford Road in Aotearoa/New Zealand. Holding the RAMMS parameters constant while adjusting the source and spatial resolution of the DEM reveals how differences in terrain representation between the satellite photogrammetry and topographic lidar DEMs (2 m spatial resolution) affect the reliability of the simulation estimates (e.g., maximum core velocity, powder pressure, runout length, final debris pattern). At the same time, coarser representations of the terrain (5 and 15 m spatial resolution) simulate avalanches that run too far and produce a powder cloud that is too large, though with lower maximum impact pressures, compared to the actual event. The complex nature of the alpine terrain in the avalanche path (steep, rough, rock faces, treeless) makes it a suitable location to specifically test the model sensitivity to digital surface models (DSMs) where both ground and above-ground features on the topography are included in the elevation model. Considering the nature of the snowpack in the path (warm, deep with a steep elevation gradient) lying on a bedrock surface and plunging over a cliff, RAMMS performed well in the challenging conditions when using the high-resolution 2 m lidar DSM, with 99 % of the simulated debris volume located in the documented debris area.