GC11-solidearth-13, updated on 16 Mar 2023
Galileo Conference: Solid Earth and Geohazards in the Exascale Era
© Author(s) 2023. This work is distributed under
the Creative Commons Attribution 4.0 License.

Ash fallout long term probabilistic volcanic hazard assessment for Neapolitan volcanoes: an example of what Earth Scientists can do with HPC resources

Manuel Stocchi1, Silvia Massaro1,2, Beatriz Martínz Montesinos2, Laura Sandri2, Jacopo Selva3, Roberto Sulpizio1,2,4, Biagio Giaccio4, Massimiliano Moscatelli4, Edoardo Peronace4, Marco Nocentini4,5, Roberto Isaia6, Manuel Titos Luzón7, Pierfrancesco Dellino1, Giuseppe Naso8, and Antonio Costa2
Manuel Stocchi et al.
  • 1Università degli Studi Aldo Moro di Bari, Department of Earth and geo-environmental sciences, Italy (manuel.stocchi@uniba.it)
  • 2Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Bologna, Italy
  • 3Dipartimento di Scienze della Terra, dell’Ambiente e delle Risorse, Università Federico II, Napoli, Italy
  • 4Istituto di Geologia Ambientale e Geoingegneria, Consiglio Nazionale delle Ricerche, Roma, Italy
  • 5Dipartimento Servizio Geologico d’Italia, Istituto Superiore per la Protezione e Ricerca Ambientale, Rome, Italy
  • 6Istituto Nazionale di Geofisica e Vulcanologia, Osservatorio Vesuviano, Napoli, Italy
  • 7Signal Processing and Machine Learning, University of Granada, Spain
  • 8Dipartimento di Protezione Civile, Rome, Italy

The creation of hazard maps relative to volcanic phenomena requires taking into account the intrinsic complexity and variability of eruptions. Here we present an example of how HPC can allow producing a high resolution multi-source probabilistic hazard assessment due to tephra fallout over a domain covering Southern Italy.  

The three active volcanoes in the Neapolitan area, Somma-Vesuvius, Campi Flegrei and Ischia, were considered as volcanic sources. For each one, we explored  three explosive size classes (Small, Medium and Large) for Somma Vesuvius and Campi Flegrei, and one explosive size class (Large) for Ischia. For each size class, we performed 1500 numerical simulations of ash dispersion (total of 10500) using the  Fall3D (V8.0) model over a computational domain covering  Southern Italy with a 0.03° ⨉ 0.03° (~3 km ⨉ 3 km) resolution. Within each size class, the eruptive parameters have been randomly sampled from well-suited probability distributions and with different meteorological conditions, obtained by randomly sampling a day between 1990 and 2020 and retrieving the relative data from the ECMWF ERA5 database. This allowed exploring the intra-class variability and to quantify aleatoric uncertainty. The results of these simulations have been post-processed with a statistical approach by assigning a weight to each eruption (based on its eruption magnitude) and the annual eruption rate of each size class. For the case of Campi Flegrei, the variability in the eruptive vent position has also been explored by constructing a grid of possible vent locations with different spatial probability. By merging the results obtained for each source and size class we produced a portfolio of hazard maps showing the expected mean annual frequency of overcoming selected thresholds in  ground tephra load. A disaggregation analysis has also been performed in order to understand which particular source and/or size class had the greater impact on a particular area. 

The completion of this work, considering both numerical simulations and the statistical elaboration of the results has required a total of more than 5000 core hours and the processing of more than 2TB of data, an effort that wouldn’t have been possible without the access to high level HPC resources.

How to cite: Stocchi, M., Massaro, S., Martínz Montesinos, B., Sandri, L., Selva, J., Sulpizio, R., Giaccio, B., Moscatelli, M., Peronace, E., Nocentini, M., Isaia, R., Titos Luzón, M., Dellino, P., Naso, G., and Costa, A.: Ash fallout long term probabilistic volcanic hazard assessment for Neapolitan volcanoes: an example of what Earth Scientists can do with HPC resources, Galileo Conference: Solid Earth and Geohazards in the Exascale Era, Barcelona, Spain, 23–26 May 2023, GC11-solidearth-13, https://doi.org/10.5194/egusphere-gc11-solidearth-13, 2023.