10 Jul 2023
 | 10 Jul 2023

Lahar events in the last 2,000 years from Vesuvius eruptions. Part 2: Formulation and validation of a computational model based on a shallow layer approach

Mattia de' Michieli Vitturi, Antonio Costa, Mauro A. Di Vito, Laura Sandri, and Domenico M. Doronzo

Abstract. In this paper we present a new model for the simulation of lahars, based on the depth-averaged code IMEX-SfloW2D with new governing and constitutive equations introduced to better describe the dynamics of lahars. A thorough sensitivity analysis is carried out to identify the critical processes (such as erosion and deposition) and parameters (both numerical and physical) controlling lahar runout, using both synthetic and real cases topographies. In particular, an application of the model to a syn-eruptive lahar from a reference-size eruption from Somma-Vesuvius, affecting the Campanian Plain (Southern Italy), described in Di Vito et al. (this issue), is used in this work for the sensitivity analysis. Effects of erosion and deposition are investigated by comparing simulations with and without these processes. By comparing flow thickness and area covered by the flow and their evolution with time, we show that the modelling of both the processes is important to properly simulate the effects of the bulking and debulking and the associated changes in rheology. From a computational point of view, the comparison of simulations obtained for different numerical grids (from 25 m to 100 m), scheme order, and grain size discretization were useful to find a good compromise between resolution and computational speed. The companion paper by Sandri et al. (this issue) shows an application of the presented model for probabilistic volcanic hazard assessment for lahars from Vesuvius deposits in the Neapolitan area.

Mattia de' Michieli Vitturi et al.

Status: final response (author comments only)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
  • RC1: 'Comment on egusphere-2023-1301', Anonymous Referee #1, 05 Sep 2023
    • AC1: 'Reply on RC1', Mattia de’ Michieli Vitturi, 21 Nov 2023
  • RC2: 'Comment on egusphere-2023-1301', Gordon Woo, 03 Oct 2023
    • AC2: 'Reply on RC2', Mattia de’ Michieli Vitturi, 21 Nov 2023

Mattia de' Michieli Vitturi et al.

Model code and software

IMEX-SfloW2D-v2 Mattia de' Michieli Vitturi

Mattia de' Michieli Vitturi et al.


Total article views: 285 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
195 71 19 285 17 11
  • HTML: 195
  • PDF: 71
  • XML: 19
  • Total: 285
  • BibTeX: 17
  • EndNote: 11
Views and downloads (calculated since 10 Jul 2023)
Cumulative views and downloads (calculated since 10 Jul 2023)

Viewed (geographical distribution)

Total article views: 271 (including HTML, PDF, and XML) Thereof 271 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
Latest update: 11 Dec 2023
Short summary
We present a numerical model for lahar generated by the mobilization of tephra deposits from a Reference-Size eruption at Somma-Vesuvius. The paper present the model (pyhsics and numerics) and a sensitivity analysis to the processes modeled, the numerical schemes and the grid resolution. This work provide the basis for the application to hazard quantification for lahars in the Vesuvian area. To this end, we rely on results of the two companion papers (Part1 on field data, Part3 on hazard maps).