Tarawera volcano (New Zealand) is volumetrically dominated by rhyolitic lavas and pyroclastic deposits, but the most recent event in AD 1886 was a basaltic Plinian fissure eruption. In March 2019 a swarm of at least 64 earthquakes occurred to the NE of Tarawera volcano, as recorded by the New Zealand Geohazard Monitoring Network (GeoNet). We use seismological analysis to show that this swarm was most likely caused by a dyke that intruded into the brittle crust between depths of 8–10 km and propagated toward Tarawera volcano for 2 km at a rate of 0.3–0.6 m s−1. We infer that this was a dyke of basaltic composition that was stress-guided toward Tarawera volcano by the topographic load of the volcanic edifice. Dyke intrusions of this nature are most likely a common occurrence but a similar process may have occurred during the 1886 eruption with a dyke sourced from some lateral distance away from the volcano. The 2019 intrusion was not detected by InSAR geodesy and we use synthetic models to show that geodetic monitoring could only detect a ≥6 m wide dyke at these depths. Improvements to geodetic monitoring, combined with detailed seismological analysis, could better detect future magmatic intrusions in the region and serve to help assess ongoing changes in the magmatic system and the associated possibilities of a volcanic event.

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.

No abstract

No abstract

No abstract