Numerical modeling of surface deformation and mechanical stability of Vesuvius volcano, Italy

Abstract

Abstract. This study was undertaken with the aim of contributing to the risk evaluation of Vesuvius. We calculate the surface displacements due to an increase in pressure in a shallow reservoir and assess the mechanical instability of the volcanic edifice and of the feeding system caused by pressure on the reservoir's walls and by regional stresses. We consider axisymmetric models that take into account topography, gravity, homogeneous and heterogeneous country rocks, spheroidal magma chambers with different aspect ratios and variable depths, and subvertical intrusions from the top of the chamber to various depths. We impose both symmetric and asymmetric regional stresses increasing with depth as boundary conditions. The models are static. The ground deformation and the stress distribution are calculated, in the framework of linear elasticity, by a numerical finite element method. The surface displacements refer to an overpressure of 10 MPa. This value is considered an upper limit for the fracture of surrounding rocks. We find the maximum vertical displacements to be of a few centimeters in the most favorable case and the displacement gradients to be at the lower limit of measurability. To evaluate the mechanical stability, we calculate the stress distributions of a prolate ellipsoidal reservoir within a heterogeneous medium. We consider a total hydrostatic magmatic pressure starting from the value of 50 MPa at the reservoir's top and three regional stress regimes from symmetric tensile to asymmetric tensile to strike-slip. In the last two cases we use three-dimensional models. The criteria adopted for instability are (1) the development of tensile tangential stress and (2) the Navier-Coulomb criterion, in compression. In no case have we found instability near or on the wall of the reservoir, whereas the slope of the volcanic edifice exhibits a shear failure instability, which increases with greater regional stress anisotropy.

Effects of the hydrothermal circulation on the strain field of the Campanian Plain (southern Italy)

“…Studies on individual volcanoes (scale of ∼10 1−2 km 2 ) show that magmatic reservoirs and/or hydrothermal systems may induce variations in the near‐field strain (e.g. Russo et al ., 1997). However, little is known about the effects of an extensive hydrothermal circulation on the strain field at a larger (∼10 3−4 km 2 ), regional scale.…”

Precise, absolute earthquake location under Somma-Vesuvius volcano using a new three-dimensional velocity model

“…Russo et al . (1997) studied the mechanical stability of Mt Vesuvius volcano through axisymmetric, numerical, finite element modelling of the stress distribution inside and beneath the volcano edifice due to the combined effects of local and regional stress sources.…”

Evolution of large Venusian volcanoes: Insights from coupled models of lithospheric flexure and magma reservoir pressurization

“…Elaborate numerical models that incorporate contributing tectonic sources, topographic loading, and/or lithospheric layers of different material properties have been used to predict stress orientations near volcanic centers [e.g., McGovern and Solomon , 1998; McGovern et al , 2001; Cailleau et al , 2005; Manconi et al , 2007; Long and Grosfils , 2009]. Authors have also explored the stability of magma chambers due to pressurization [e.g., Sartoris et al , 1990; Parfitt et al , 1993] and volcanic edifice loading [e.g., Russo et al , 1997; Pinel and Jaupart , 2003]. However, recent numerical modeling approaches that incorporate gravitational loading and lithostatic prestress conditions predict significantly different patterns of maximum stresses and strains [e.g., Grosfils , 2007; Hurwitz et al , 2009; Long and Grosfils , 2009] compared to earlier efforts.…”