Modelling Pre-eruptive Progressive Damage in Basaltic Volcanoes: Consequences for the Pre-eruptive Process
Abstract
At Grimsvötn volcano, high-quality earthquake and continuous GPS data were recorded by the Icelandic Meteorological Office during its 2004-2011 inter-eruptive period and exhibited remarkable patterns : acceleration of the cumulated earthquake number, and a 2-year exponential decrease in displacement rate followed by a 4-year constant inflation rate. We proposed a model with one magma reservoir in a non-linear elastic damaging edifice, with incompressible magma and a constant pressure at the base of the magma conduit. We first modelled seismicity rate and damage as a function of time, and show that Kachanov's elastic brittle damage law may be used to express the decrease of the effective shear modulus with time. We then derived simple analytical expressions for the magma reservoir overpressure and the surface displacement as a function of time. We got a very good fit of the seismicity and surface displacement data by adjusting only three phenomenological parameters and computed magma reservoir overpressure, magma flow and strain power as a function of time. Overpressure decrease is controlled by damage and shear modulus decrease. Displacement increases, although overpressure is decreasing, because shear modulus decreases more than overpressure. Normalized strain power reaches a maximum 0.25 value. This maximum is a physical limit, after which the elasticity laws are no longer valid, earthquakes cluster, cumulative number of earthquakes departs from the model. State variable extrema provide four reference times that may be used to assess the mechanical state and dynamics of the volcanic edifice. We also performed the spatial modelling of the progressive damage and strain localization around a pressurized magma reservoir. We used Kachanov's damage law and finite element modelling of an initially elastic volcanic edifice pressurized by a spherical magma reservoir, with a constant pressure in the reservoir and various external boundary conditions. At each node of the model, Young's modulus is decreased if deviatoric stress locally reaches the Mohr-Coulomb plastic threshold. For a compressive horizontal stress, the result shows a complex strain localization pattern, showing reverse and normal faulting very similar to what is obtained from analog modelling and observed at volcanic resurgent domes.
- Publication:
-
AGU Fall Meeting Abstracts
- Pub Date:
- December 2017
- Bibcode:
- 2017AGUFM.V34C..02G
- Keywords:
-
- 1207 Transient deformation;
- GEODESY AND GRAVITY;
- 8159 Rheology: crust and lithosphere;
- TECTONOPHYSICS;
- 8434 Magma migration and fragmentation;
- VOLCANOLOGY;
- 8485 Remote sensing of volcanoes;
- VOLCANOLOGY