Stress evolution during caldera collapse

- Distinct Element Method modelling of stress evolution during caldera collapse.
- Explicit fracturing due to gravity-driven stress changes from reservoir depletion.
- Rotation of principal stress orientations during collapse causes changes in ring fault type.
- Four end-member stress paths taken to failure in reservoir roof.
- Reservoir unloading-reloading linked to roof geometric and mechanical properties.

The mechanics of caldera collapse are subject of long-running debate. Particular uncertainties concern how stresses around a magma reservoir relate to fracturing as the reservoir roof collapses, and how roof collapse in turn impacts upon the reservoir. We used two-dimensional Distinct Element Method models to characterise the evolution of stress around a depleting sub-surface magma body during gravity-driven collapse of its roof. These models illustrate how principal stress orientations rotate during progressive deformation so that roof fracturing transitions from initial reverse faulting to later normal faulting. They also reveal four end-member stress paths to fracture, each corresponding to a particular location within the roof. Analysis of these paths indicates that fractures associated with ultimate roof failure initiate in compression (i.e. as shear fractures). We also report on how mechanical and geometric conditions in the roof affect pre-failure unloading and post-failure reloading of the reservoir. In particular, the models show how residual friction within a failed roof could, without friction reduction mechanisms or fluid-derived counter-effects, inhibit a return to a lithostatically equilibrated pressure in the magma reservoir. Many of these findings should be transferable to other gravity-driven collapse processes, such as sinkhole formation, mine collapse and subsidence above hydrocarbon reservoirs.