Temporal evolution and compositional signatures of two supervolcanic systems recorded in zircons from Mangakino volcanic centre, New Zealand

Mangakino, the oldest rhyolitic caldera centre delineated in the Taupo Volcanic Zone of New Zealand, generated two very large (super-sized) ignimbrite eruptions, the 1.21 ± 0.04 Ma >500 km3 Ongatiti and ~1.0 Ma ~1,200 km3 Kidnappers events, the latter of which was followed after a short period of erosion by the ~200 km3 Rocky Hill eruption. We present U/Pb ages and trace-element analyses on zircons from pumice clasts from these three eruptions by Secondary Ion Mass Spectrometry (SIMS) using SHRIMP-RG instruments to illustrate the evolution of the respective magmatic systems. U–Pb age spectra from the Ongatiti imply growth of the magmatic system over ~250 kyr, with a peak of crystallisation around 1.32 Ma, ~100 kyr prior to eruption. The zircons are inferred to have then remained stable in a mush with little crystallisation and/or dissolution before later rejuvenation of the system at the lead-in to eruption. The paired Kidnappers and Rocky Hill eruptions have U–Pb zircon ages and geochemical signatures that suggest they were products of a common system grown over ~200 kyr. The Kidnappers and Rocky Hill samples show similar weakly bimodal age spectra, with peaks at 1.1 and 1.0 Ma, suggesting that an inherited antecrystic population was augmented by crystals grown at ages within uncertainty of the eruption age. In the Kidnappers, this younger age peak is dominantly seen in needle-shaped low U grains with aspect ratios of up to 18. In all three deposits, zircon cores show larger ranges and higher absolute concentrations of trace elements than zircon rims, consistent with zircon crystallisation from evolving melts undergoing crystal fractionation involving plagioclase and amphibole. Abundances and ratios of many trace elements frequently show variations between different sectors within single grains, even where there is no visible sector zoning in cathodoluminescence (CL) imaging. Substitution mechanisms, as reflected in the molar (Sc + Y + REE3+)/P ratio, differ in the same growth zone between the sides (along a-axis and b-axis: values approaching 1.0) and tips (c-axis: values between 1.5 and 5.0) of single crystals. These observations have implications for the use of zircons for tracking magmatic processes, particularly in techniques where CL zonation within crystals is not assessed and small analytical spot sizes cannot be achieved. These observations also limit applicability of the widely used Ti-in-zircon thermometer. The age spectra for the Ongatiti and Kidnappers/Rocky Hill samples indicate that both magmatic systems were newly built in the time-breaks after respective previous large eruptions from Mangakino. Trace element variations defining three-component mixing suggest that zircons, sourced from multiple melts, contributed to the population in each system.