Particle Sedimentation in Numerical Modelling: A Case Study from the Puyehue-Cordón Caulle 2011 Eruption with the PLUME-MoM/HYSPLIT Models

• Published in: Atmosphere Vol. 13; no. 5; p. 784
• Main Authors: Tadini, Alessandro, Gouhier, Mathieu, Donnadieu, Franck, de’ Michieli Vitturi, Mattia, Pardini, Federica
• Format: Journal Article
• Language: English
• Published: Basel MDPI AG 01-05-2022; MDPI
• Subjects: Aggregation; Atmospheric models; Datasets; Earth Sciences; Eruptions; Fallout; Grain size; Hazard assessment; HYSPLIT; Magma; Mass; Mathematical models; Model testing; Modelling; numerical modelling; Parameters; Parametric analysis; Particle settling; Phase transitions; PLUME-MoM; Procedures; Sciences of the Universe; Sedimentation & deposition; Settling rate; Settling velocity; Tephra; tephra fallout; Volcanic eruptions; Volcanology; numerical modelling; tephra fallout; PLUME-MoM; HYSPLIT
• ISSN: 2073-4433; 2073-4433
• DOI: 10.3390/atmos13050784

Numerical modelling of tephra fallout is a fast-developing research area in volcanology. Several models are currently available both to forecast the dispersion of volcanic particles in the atmosphere and to calculate the particles deposited at different locations on the ground. Data from these simulations can then be used both to manage volcanic crises (e.g., protect air traffic) or perform long-term hazard assessment studies (e.g., through hazard maps). Given the importance of these tasks, it is important that each model is thoroughly tested in order to assess advantages and limitations, and to provide useful information for quantifying the model uncertainty. In this study we tested the coupled PLUME-MoM/HYSPLIT models by applying them to the Puyehue–Cordon Caulle 2011 sub-Plinian eruption. More specifically, we tested new features recently introduced in these well-established models (ash aggregation, external water addition, and settling velocity models), we implemented a new inversion procedure, and we performed a parametric analysis. Our main results reaffirm the pivotal role played by mass eruption rate on the final deposit and show that some choices for the input parameters of the model can lead to the large overestimation in total deposited mass (which can be reduced with our inversion procedure). The parametric analysis suggests a most likely value of the mass eruption rate in the range 2.0–6.3 × 106 kg/s. More studies with a similar approach would be advisable in order to provide final users with useful indications about the parameters that should be carefully evaluated before being used as input for this kind of model.