Reinterpretation of the post 26ka Taupō rhyolitic magmatic system (New Zealand) as deep and vertically extensive based on Isotope thermometry and measured and modeled zircon destinies

Abstract Taupō volcanic zone, the site of the 26 ka Oruanui supereruption, produced ~70 km3 of new rhyolites since 11 ka, culminating in 50 km3 Taupō eruption 1.8 ka. Major phenocrysts decrease from 4 to 1 vol%, and Oruanui and post-Oruanui ignimbrites all have identical high-δ18Omelt values of 7.39 ± 0.1‰ and lack low-δ18O values despite overlapping calderas. The Δ’17O values are −0.07‰, lower than the mantle and indicate source contamination of high-δ18O, low-Δ’17O metasediments, and limited interaction with high-Δ’17O hydrothermally altered crust. Previously published U-Th-Pb zircon ages demonstrate their diversity spanning 104–105 years for each unit. Zircon crystal size distribution shows a decrease in abundance and the mean size, and some units lack small (~10 um) zircons suggesting that zircons were both growing and dissolving in the coexisting magma generation areas. Isotope thermometry indicates heating of the system from ~812 ± 35°C to 874 ± 36°C past zircon saturation in 1.8 ka eruption. We advocate that a deep vertically continuous and laterally discontinuous silicic magma system at the base of the Taupō rift, rather than a shallow batholith or an evolving mush, drives volcanism at Taupō. To explain the post-Oruanui magma production, rift-base silicic magma origin and moderate (~2 km3/1000 years) rhyodacitic magma flux from a growing and heating liquid magma body creates a sufficient solution for the most recent magmatism.