Complex decompression and fragmentation of mingled andesite magmas driving multi-phase Plinian eruptions at Mt. Taranaki, New Zealand

Abstract

Estimating the kinetics of andesite magma vesiculation and crystallization inside volcanic plumbing systems is key for unraveling andesite Plinian eruption dynamics. The conduit kinetics provide the necessary input data for estimating the magma flow rates driving magma ascent and the fragmentation mechanisms controlling shifts in eruption explosivity and style. This information is crucial for increasing knowledge on expected hazards and for developing realistic eruption scenarios. In this work, we estimate conduit magma vesiculation and crystallization kinetics during the 3300 cal BP Upper Inglewood Plinian eruptive episode of Mount Taranaki, New Zealand. This episode comprised (i) low-intensity, conduit-opening phases of dome-collapse PDCs; (ii) pre-climactic, highly explosive phases of diverse PDCs, of up to violent 18-km-runout lateral blasts; (iii) climactic phases of steady 22-km-high Plinian eruption columns; and (iv) waning phases of column-collapse PDCs. By employing synchrotron microtomography, combined with mineral/glass chemistry and electron-microscopy, we quantified 3D vesicle and crystal size and shape distributions in juvenile pyroclasts over time, and corresponding number densities ranging from 1.1 × 105 to 2.5 × 106 mm−3 for vesicles, and from 8.0 × 104 to 5.1 × 106 mm−3 for crystals. Our results indicate that tapping of chemically alike yet rheologically contrasting magmas over a multi-phase andesite eruptive episode is linked to: (a) mafic magma recharge and differentiation in multiple storage reservoirs at distinct crustal levels, (b) stepwise to rapid magma decompression while mingling, producing variable pre- and syn-eruptive degassing and crystallization, and (c) syn-eruptive changes in melt viscosity, strain rate, localized shear deformation, and conduit geometry. The earliest and least explosive eruptive phases (≈ 2 × 106 kg s−1) were produced at the slowest rates of magma decompression (0.3–0.6 MPa s−1), ascent (0.01–0.02 m s−1) and strain (< 0.002 s−1), driven by volatile diffusion and exsolution. All subsequent pre-climactic and Plinian phases (4 × 107–1 × 108 kg s−1) were produced at either rapid or intermittent rates of magma decompression (2.0–6.0 MPa s−1), ascent (0.06–0.2 m s−1) and strain (> 0.003–0.010 s−1), powered by combined magma volatile supersaturation and delayed disequilibrium degassing, decompression-induced microlite crystallization and rapid heterogeneous vesiculation kinetics, shear deformation and magma mingling. These processes enabled complex fragmentation mechanisms of the rheologically most homogeneous magmas.