Technical Papers and Presentations

Campi Flegrei Caldera is an active volcano located westward of Naples, Italy. Despite its volcanic activity spanning over 39,000 years, the area is densely populated and poses a significant threat to the inhabitants due to ongoing seismicity, ground uplift, and hydrothermal activity resulting from increased pressure due to fluid injections and thermal variations at depth. These factors highlight the potential for a major volcanic eruption. Being able to model the processes that lead to ground deformation could be of vital importance in predicting an imminent eruption and enabling the evacuation of the resident population in the area. Such models are being schematized and analyzed using specific Finite Element Method (FEM) software capable of solving complex mathematical problems involving partial differential equations. In this manuscript, we utilized the COMSOL Multiphysics® software, configured with the Structural Mechanics Module, to examine how deep changes in pressure and temperature influence the observed surface deformation field within the Campi Flegrei caldera. To simulate these deep changes, we employed the Tough software, an open-source numerical simulation program designed for multi-dimensional fluid and heat flows of multiphase and multicomponent fluid mixtures in porous and fractured media. The output of Tough served as an input for the COMSOL® model, representing the source at depth. By using the Structural Mechanics Module, we were able to assess the accuracy of the proposed model in comparison to analytical solutions. Furthermore, we were able to model the geometry of the deep source in more detail and verify that the surface deformation pattern aligned with the measurements obtained from sensors. The surface deformations obtained through simulation corresponded with those observed through GPS/GNSS time series. By leveraging COMSOL Multiphysics®, therefore, we have constructed a mathematical model that accurately captures the intricate interplay of fluid injections, thermal variations, and rock mechanics, enabling us to simulate volcanic crustal deformations with remarkable fidelity. Additionally, we provide an explanation for the mini-uplifts recorded by various geodetic networks, managed by Osservatorio Vesuviano, during the period from 2006 to 2022.