3D integrated geophysical modeling for the 2008 magma intrusion at Etna: Constraints on rheology and dike overpressure

We present a 3D numerical model based on Finite Element Method (FEM) to jointly evaluate geophysical
changes caused by dislocation and overpressure sources in volcanic areas. A coupled numerical problem
was solved to estimate ground deformation, gravity and magnetic changes produced by stress redistribution
accompanying magma migration within the volcano edifice. We successfully applied the integrated
numerical procedure to image the magmatic intrusion occurring in the northern flank of Etna during
the onset of the 2008 eruption. A multi-layered crustal structure of the volcano constrained by geological
models and geophysical data was considered. Geodetic and gravity data provide information on the
strain field, while piezomagnetic changes give constraints on the stress field. Therefore, the integrated
modeling gives insights on Mt Etna rheology and dike overpressure involved in the magma propagation
and improves understanding of dike emplacement in the northern sector of the volcano. Our FEM-based
approach improves the reliability of model-based inference of geophysical parameters obtained during
monitoring of the onset of Etna lateral intrusions that can prelude to an impending eruption.