James Hickey

Estimating volcanic deformation source parameters with a finite element inversion: The 2001–2002 unrest at Cotopaxi volcano, Ecuador

Deformation at Cotopaxi was observed between 2001 and 2002 along with recorded seismicity beneath the northeast (NE) flank, despite the fact that the last eruption occurred in 1942. We use electronic distance meter deformation data along with the patterns of recorded seismicity to constrain the cause of this unrest episode. To solve for the optimum deformation source parameters we employ inverse finite element (FE) models that account for material heterogeneities and surface topography. For a range of source shapes the models converge on a shallow reservoir beneath the southwest (SW) flank. The individual best fit model is a small oblate-shaped source, approximately 4-5 km beneath the summit, with a volume increase of roughly 20 × 10 6 m 3 . This SW source location contrasts with the NE seismicity locations. Subsequently, further FE models that additionally account for temperature-dependent viscoelasticity are used to reconcile the deformation and seismicity simultaneously. Comparisons of elastic and viscous timescales allude to aseismic pressurization of a small magma reservoir in the SW. Seismicity in the NE is then explained through a mechanism of fluid migration from the SW to the NE along fault systems. We extend our analyses to further show that if future unrest crises are accompanied by measurable seismicity around the deformation source, this could indicate a higher magma supply rate and increased likelihood of a forthcoming eruption.

Thermomechanical controls on magma supply and volcanic deformation: application to Aira caldera, Japan.

Ground deformation often precedes volcanic eruptions, and results from complex interactions between source processes and the thermomechanical behaviour of surrounding rocks. Previous models aiming to constrain source processes were unable to include realistic mechanical and thermal rock properties, and the role of thermomechanical heterogeneity in magma accumulation was unclear. Here we show how spatio-temporal deformation and magma reservoir evolution are fundamentally controlled by three-dimensional thermomechanical heterogeneity. Using the example of continued inflation at Aira caldera, Japan, we demonstrate that magma is accumulating faster than it can be erupted, and the current uplift is approaching the level inferred prior to the violent 1914 Plinian eruption. Magma storage conditions coincide with estimates for the caldera-forming reservoir ~29,000 years ago, and the inferred magma supply rate indicates a ~130-year timeframe to amass enough magma to feed a future 1914-sized eruption. These new inferences are important for eruption forecasting and risk mitigation, and have significant implications for the interpretations of volcanic deformation worldwide.

The Ups and Downs of Volcanic Unrest: Insights from Integrated Geodesy and Numerical Modelling

Volcanic eruptions are often preceded by small changes in the shape of the volcano. Such volcanic deformation may be measured using precise surveying techniques and analysed to better understand volcanic processes. Complicating the matter is the fact that deformation events (e.g., inflation or deflation) may result from magmatic, non-magmatic or mixed/hybrid sources. Using spatial and temporal patterns in volcanic deformation data and mathematical models it is possible to infer the location and strength of the subsurface driving mechanism. This can provide essential information to inform hazard assessment, risk mitigation and eruption forecasting. However, most generic models over-simplify their representation of the crustal conditions in which the deformation source resides. We present work from a selection of studies that employ advanced numerical models to interpret deformation and gravity data. These incorporate crustal heterogeneity, topography, viscoelastic rheology and the influence of temperature, to constrain unrest source parameters at Uturuncu (Bolivia), Cotopaxi (Ecuador), Soufriere Hills (Montserrat), and Teide (Tenerife) volcanoes. Such model complexities are justified by geophysical, geological, and petrological constraints. Results highlight how more realistic crustal mechanical conditions alter the way stress and strain are partitioned in the subsurface. This impacts inferred source locations and magmatic pressures, and demonstrates how generic models may produce misleading interpretations due to their simplified assumptions. Further model results are used to infer quantitative and qualitative estimates of magma supply rate and mechanism, respectively. The simultaneous inclusion of gravity data alongside deformation measurements may additionally allow the magmatic or non-magmatic nature of the source to be characterised. Together, these results highlight how models with more realistic, and geophysically consistent, components can improve our understanding of the mechanical processes affecting volcanic unrest and geodetic eruption precursors, to aid eruption forecasting, hazard assessment and risk mitigation.