Pre-eruptive rhyolite magma ascent rate is rapid and independent of eruption size: a case study from Ōkataina Volcanic Centre, Aotearoa New Zealand

Abstract

Volatile measurements in mineral-hosted sealed melt inclusions, and open-ended embayments, have previously been used to study magma ascent dynamics in large rhyolitic eruptions. However, despite occurring more frequently, smaller-volume explosive events remain under-studied. We present magmatic volatile data from quartz-hosted melt inclusions and embayments for eight post-25.4 ka rhyolitic eruptions at Ōkataina Volcanic Centre, Aotearoa New Zealand. Seven originated from within the main caldera, and the other erupted from the associated Ōkareka Structural Embayment. Melt inclusions preserve volatile contents of 2.92–5.82 wt% H2O and 13–126 ppm CO2, indicating pre-eruptive storage depths of 4.5–7.4 km, with younger eruptions being more shallow. The lack of correlation between H2O, CO2, inclusion size or distance to the crystal rim suggests magma bodies experienced variable degrees of degassing during magma storage, with some amount of post-entrapment volatile modification prior to and concurrent with final magma ascent. Diffusion modelling of measured H2O gradients in melt embayments indicates ascent rates of 0.10–1.67 m.s−1 over time spans of 20–230 min for the intra-caldera events. In contrast, ascent rates for the eruption from the Ōkareka Structural Embayment may be more rapid, at 1.59–4.4 m.s−1 over a time span of 22–34 min. Our findings imply that the final, pre-eruptive magma movement towards the surface could be less than a few hours. Comparisons with published data for caldera-forming explosive events reveal no clear relationships between final ascent rate, eruption size or initial volatile content, implying that other factors besides eruption volume control rhyolite magma ascent.

Deciphering Clues Regarding Magma Composition Encoded in Quartz‐Hosted Embayments and Melt Inclusions Through Direct Numerical Simulations

“…Previous studies propose that higher volatile contents in melt inclusions relative to embayment interiors in the same crystal are indicative of decompression between the formation of melt inclusions and embayments (Elms et al., 2023; Myers et al., 2016, 2018, 2021; Saalfeld et al., 2021). However, our results demonstrate that for the specific case modeled here, namely embayment formation during quartz growth, different post‐formation connectivity with surrounding melt can contribute to the compositional difference.…”

“…Embayments have been increasingly used as a tool for estimating magma ascent rates (Befus et al, 2023;Elms et al, 2023;Ferguson et al, 2016;Geshi et al, 2021;Hosseini et al, 2023;Humphreys et al, 2008;Liu et al, 2007;Lloyd et al, 2014;Moussallam et al, 2019;Myers et al, 2016Myers et al, , 2018Myers et al, , 2021Newcombe et al, 2020;Saalfeld et al, 2021;Zuccarello et al, 2022). During magma ascent, the volatile concentrations of a melt decrease due to volatile saturation and exsolution during decompression.…”