Defining a 3D physical model for the hydrothermal circulation at Campi Flegrei caldera (Italy)

Our study is aimed to develop a 3D physical model of the Campi Flegrei geothermal system, in order to achieve a more accurate and comprehensive representation of the hydrothermal processes occurring in the caldera. The new model, developed by using the TOUGH2 code simulator, accounts for the caldera rocks' physical properties, bathymetry and water table topography. In particular, the computational domain is constrained by density values obtained by tomographic inversion of gravity data collected during several surveys at CF both onshore and offshore the caldera. Empirical relations between density and porosity and between porosity and permeability, derived by published data on samples cored in deep wells or collected in outcrops, allowed us to characterize the main rocks physical parameters controlling the dynamic of the CF geothermal system. We have performed and compared several simulations investigating the effects of the injection at depth, underneath Solfatara crater, of a hot gaseous mixture rich in CO2. We show that, with respect to the available literature on 2D axisymmetric models, the effects of the water table topography together with the bathymetry and the heterogeneous distribution of the rocks' physical properties, lead to important differences in the hydrothermal circulation of fluids at CF. These constraints allow the activation of convective cells with different behaviors, which produce variable patterns of temperature inside the hydrothermal system. As a consequence, the predominant effect is again represented by a central plume below the Solfatara crater, but with a non-axisymmetric structure and a wider extension. Additionally, high temperature zones are present near the coastline and in the middle part of the submerged area of the caldera with a SE–NW alignment. Moreover, our results indicate that, the submerged part of the CF caldera would deserve a more accurate study and survey, being affected by phenomenon of heating and degassing. This information could be very useful in terms of hazard assessment and volcanic risk mitigation in such an active and densely inhabited volcanic and geothermal area