Numerical modeling of debris avalanche propagation from collapse of volcanic edifices

Debris avalanches produced from the collapse of volcanic edifices are destructive events that involve volumes up to two orders of magnitude larger (cubic kilometer) than most non-volcanic rock and debris avalanches. We replicate the motion and spreading of several volcanic collapses by means of a depth-averaged quasi-3D numerical code. The model assumes a frictional internal rheology and a variable basal rheology (i.e frictional, Voellmy and plastic). We back analyzed seven case-studies against observations reported in the literature to provide a set of calibrated cases. The ASTER and SRTM satellite-derived digital elevation models were used as topographic data. The numerical model captures the main features of the propagation process, including travel distance, lateral spreading and run up. At varying triggering factors and material characteristics the best fitting parameters span in a narrow interval and differ from those typical of non-volcanic rock and debris avalanches. The bulk basal friction angles (the sole parameter required in the frictional rheology) range within 3° and 7.5° whereas typical values for non-volcanic debris avalanches vary from 11° to 31°. The consistency of the back analyzed parameters is encouraging for a possible use of the model in the perspective of hazard mapping. The reconstruction of the pre-event topography is critical, and it is associated to large uncertainty. The quality of the terrain data, more than the resolution of the DEMs used, is relevant for the modeling. Resampling the original square grid to larger cell sizes determines a low increase in the back analyzed rheological parameters, as a result of the lower roughness of the terrain.