Basaltic caldera-rift dike intrusions and island-scale subsidence

Shreve, T., Grandin, R., Boichu, M. et al. (2019). From prodigious volcanic degassing to caldera subsidence and quiescence at Ambrym (Vanuatu): the influence of regional tectonics. Sci Rep 9, 18868 . https://doi.org/10.1038/s41598-019-55141-7.

A conceptual model of the December 2018 submarine eruption at Ambrym volcano (Vanuatu). The inset is a zoom on Ambrym’s tectonic setting, emphasizing that Ambrym’s rift zone orientation is sub-parallel to the regional maximum compressive stress, allowing us to interpret Ambrym as a large tension fracture.

Caldera ring-fault activation measured with InSAR

Shreve, T.,  Grandin, R.,  Smittarello, D.,  Cayol, V.,  Pinel, V.,  Boichu, M., &  Morishita, Y. (2021).  What triggers caldera ring-fault subsidence at Ambrym volcano? Insights from the 2015 dike intrusion and eruption. Journal of Geophysical Research: Solid Earth,  126, e2020JB020277. https://doi.org/10.1029/2020JB020277.

Shreve, T., &  Delgado, F. (2023).  Trapdoor fault activation: A step toward caldera collapse at Sierra Negra, Galápagos, Ecuador. Journal of Geophysical Research: Solid Earth,  128, e2023JB026437. https://doi.org/10.1029/2023JB026437.


At Ambrym volcano (Vanuatu) in 2015, a reservoir (blue spheroid) at depths of 4–5 km b.s.l. fed a dike intrusion and intracaldera fissure eruption. An extraction of less than 7% of magma from the reservoir was sufficient to unclamp the western portion of Ambrym's caldera ring-fault (light green rectangle). We hypothesize that the fault was prestressed by a previous event associated with a depressurizing reservoir (due to an eruption, degassing, etc.), and the reservoir deflation in 2015 brought the fault to failure. 

Degassing-driven subsidence and "stealthy" magma influx

Shreve, T., Grandin, R., Boichu, M. (2022). Reservoir depressurization driven by passive gas emissions at Ambrym volcano. Earth and Planetary Science Letters, Volume 584, 117512, ISSN 0012-821X. https://doi.org/10.1016/j.epsl.2022.117512.


A schematic showing a theoretical model used to explain the end of subsidence at Ambrym volcano in 2017, after 2 years of deformation (modified after Girona et al. (2014)). This model assumes a reservoir is connected to the surface by an open conduit, and that the pressure in the reservoir equilibrates the weight of the magma in the conduit. Consequently, mass changes within the system (either in the conduit or reservoir) due to degassing result in reservoir depressurization, causing ground subsidence.

Evidence for multiple magma storage regions from InSAR time series

Shreve, T.,  Zhan, Y.,  Le Mével, H.,  Roman, D., &  Moussallam, Y. (2023).  Two distinct magma storage regions at Ambrym volcano detected by satellite geodesy. Geophysical Research Letters,  50, e2023GL102925. https://doi.org/10.1029/2023GL102925.


An interpretation of Ambrym's magmatic plumbing system after the December 2018 eruption. The annotated sketch includes variables from a theoretical model discussed in this study. Using geodetic modelling and a data assimilation approach, we estimate values of pressure change over time for two separate reservoirs. Furthermore, we apply a theoretical model to demonstrate that the reservoirs may not currently be hydraulically connected, despite evidence of physical mixing of magma derived from each reservoir during the December 2018 eruption. 

Heat and gas transfer through a quiescent lava dome

Shreve, T., Girona, T., & Anderson, K (2023).  Time-dependent Numerical Modelling of Gas and Heat Transfer Through a Solidified Lava Dome. American Geophysical Union Fall Meeting, December 2023. San Francisco, USA. Poster Presentation. https://agu.confex.com/agu/fm23/meetingapp.cgi/Paper/1347661.

This current work applies the finite element method to investigate how changes in surface temperature relate to pore pressure conditions beneath and within a quiescent lava dome. Preliminary results find that high surface temperatures alone do not necessarily indicate high overpressures and pre-eruptive conditions, and that distributed regions of moderately to highly elevated temperatures may indicate increased dome overpressure. Time-dependent models also indicate that a dome permeability decrease may cause a local and temporary decrease in surface temperature.