The role of magma chamber-fault interaction in caldera forming eruptions

• Published in: Bulletin of volcanology 2010, Vol.72 (1), p.85-101
• Main Authors: Simakin, A. G., Ghassemi, A.
• Format: Article
• Language: English
• Subjects: Calderas; Deformation; Earth and Environmental Science; Earth Sciences; Geology; Geophysics; Geophysics/Geodesy; Magma; Mechanical failure; Mineralogy; Plate tectonics; Research Article; Rheology; Rocks; Sedimentology; Strain; Volcanoes; Volcanology
• ISSN: 0258-8900; 1432-0819
• DOI: 10.1007/s00445-009-0306-6

This paper examines the role of the position and orientation of a regional fault in the roof of a magma chamber on stress distribution, mechanical failure, and dyking using 2D finite element numerical simulations. The study pertains to the magma chamber behavior in the relatively short time intervals of several hundreds to thousand of years. The magma chamber is represented as an elliptical inclusion (eccentricity, a/b  = 0.12) at a relative depth, H/a , of 0.9. The fault has a 45° dip and is represented by a frictionless fracture. The temperature field in the host rock is calculated assuming a quasisteady-state thermal regime that develops through periodic episodes of magma supply. The rheology of the surrounding rocks is treated using viscoelasticity with temperature activated strain-rate dependent viscosity. Strain weakening of the rocks in the ductile zone is described within the frame of the Dynamic Power Law model . The magma pressure is coupled with the deformation of the rock mass hosting the chamber, including the fault. The variation of magma pressure in response to magma supply and chamber deformation is calculated in the elastic and viscoelastic regimes. The latter corresponds to slow filling, while the former represents a filling time much less than the viscous relaxation time scale. The resulting “equation of state” for the magma chamber couples the magma pressure with the chamber volume in the elastic regime, and with the filling rate for the viscoelastic regime. Analysis of stresses is used to predict dyke propagation conditions, and the mechanical failure of the chamber roof for different fault positions and magma overpressures. Results show that an outward dipping fault located on the periphery of the chamber roof hinders the propagation of dykes to the surface, causing magma to accumulate under the footwall of the fault. At high to moderate overpressures (30–40 MPa), the fault causes localized shear failure and chamber roof collapse that might lead to the first stage of a caldera-forming eruption.