Explore the vast range of titles available.

Supereruptions (>10 15 kg ≈ 450 km 3 of ejected magma) have received much attention because of the challenges in explaining how and over what time intervals such large volumes of magma are accumulated, stored and erupted. However, the processes that follow supereruptions, particularly those focused around magmatic recovery, are less fully documented. We present major and trace-element data from whole-rock, glass and mineral samples from eruptive products from Taupo volcano, New Zealand, to investigate how the host magmatic system reestablished and evolved following the Oruanui supereruption at 25.4 ka. Taupo’s young eruptive units are precisely constrained chronostratigraphically, providing uniquely fine-scale temporal snapshots of a post-supereruption magmatic system. After only ~5 kyr of quiescence following the Oruanui eruption, Taupo erupted three small volume (~0.1 km 3 ) dacitic pyroclastic units from 20.5 to 17 ka, followed by another ~5-kyr-year time break, and then eruption of 25 rhyolitic units starting at ~12 ka. The dacites show strongly zoned minerals and wide variations in melt-inclusion compositions, consistent with early magma mixing followed by periods of cooling and crystallisation at depths of >8 km, overlapping spatially with the inferred basal parts of the older Oruanui silicic mush system. The dacites reflect the first products of a new silicic system, as most of the Oruanui magmatic root zone was significantly modified in composition or effectively destroyed by influxes of hot mafic magmas following caldera collapse. The first rhyolites erupted between 12 and 10 ka formed through shallow (4–5 km depth) cooling and fractionation of melts from a source similar in composition to that generating the earlier dacites, with overlapping compositions for melt inclusions and crystal cores between the two magma types. For the successively younger rhyolite units, temporal changes in melt chemistry and mineral phase stability are observed, which reflect the development, stabilisation and maturation of a new, probably unitary, silicic mush system. This new mush system was closely linked to, and sometimes physically interacted with, underlying mafic melts of similar composition to those involved in the Oruanui supereruption. From the inferred depths of magma storage and geographical extent of vent sites, we consider that a large silicic mush system (>200 km 3 and possibly up to 1000 km 3 in volume) is now established at Taupo and is capable of feeding a new episode or cycle of volcanism at any stage in the future.

Chemical anomalies in polar ice core records are frequently linked to volcanism; however, without the presence of (crypto)tephra particles, links to specific eruptions remain speculative. Correlating tephras yields estimates of eruption timing and potential source volcano, offers refinement of ice core chronologies, and provides insights into volcanic impacts. Here, we report on sparse rhyolitic glass shards detected in the Roosevelt Island Climate Evolution (RICE) ice core (West Antarctica), attributed to the 1.8 ka Taupō eruption (New Zealand)—one of the largest and most energetic Holocene eruptions globally. Six shards of a distinctive geochemical composition, identical within analytical uncertainties to proximal Taupō glass, are accompanied by a single shard indistinguishable from glass of the ~25.5 ka Ōruanui supereruption, also from Taupō volcano. This double fingerprint uniquely identifies the source volcano and helps link the shards to the climactic phase of the Taupō eruption. The englacial Taupō-derived glass shards coincide with a particle spike and conductivity anomaly at 278.84 m core depth, along with trachytic glass from a local Antarctic eruption of Mt. Melbourne. The assessed age of the sampled ice is 230 ± 19 CE (95% confidence), confirming that the published radiocarbon wiggle-match date of 232 ± 10 CE (2 SD) for the Taupō eruption is robust.

Key processes influencing the aggregation of volcanic ash and hydrometeors are examined with an experimental method employing vibratory pan aggregation. Mechanisms of aggregation in the presence of hail and ice pellets, liquid water (≤30 wt%), and mixed water phases are investigated at temperatures of 18 and −20 °C. The experimentally generated aggregates, examined in hand sample, impregnated thin sections, SEM imagery, and X-ray microtomography, closely match natural examples from phreatomagmatic phases of the 27 ka Oruanui and 2010 Eyjafjallajökull eruptions. Laser diffraction particle size analysis of parent ash and aggregates is also used to calculate the first experimentally derived aggregation coefficients that account for changing liquid water contents and subzero temperatures. These indicate that dry conditions (<5–10 wt% liquid) promote strongly size selective collection of sub-63 μm particles into aggregates (given by aggregation coefficients >1). In contrast, liquid-saturated conditions (>15–20 wt% liquid) promote less size selective processes. Crystalline ice was also capable of preferentially selecting volcanic ash <31 μm under liquid-free conditions in a two-stage process of electrostatic attraction followed by ice sintering. However, this did not accumulate more than a monolayer of ash at the ice surface. These quantitative relationships may be used to predict the timescales and characteristics of aggregation, such as aggregate size spectra, densities, and constituent particle size characteristics, when the initial size distribution and water content of a volcanic cloud are known. The presence of an irregularly shaped, millimeter-scale vacuole at the center of natural aggregates was also replicated during interaction of ash and melting ice pellets, followed by sublimation. Fine-grained rims were formed by adding moist aggregates to a dry mixture of sub-31 μm ash, which adhered by electrostatic forces and sparse liquid bridges. From this, we infer that the fine-grained outer layers of natural aggregates reflect recycled exposure of moist aggregates to regions of volcanic clouds that are relatively dry and dominated by <31 μm ash.

The magmatic systems that give rise to voluminous crystal-poor rhyolite magma bodies can be considered to operate on two contrasting timescales: Those governed by longer-term processes by which a magma acquires its chemical and isotopic characteristics (e.g., fractional crystallisation and assimilation), and those operating at shorter timescales during the physical accumulation of the melt-dominant magma body that finally erupts. We explore the compositional and textural relationships between amphibole and orthopyroxene crystals from the 25.4 ka, 530 km 3 (magma) Oruanui eruption products (Taupo volcano, New Zealand) to investigate how processes related to the physical assembly of the pre-eruptive magma body are represented in the crystal record. Over 90 % of orthopyroxenes from the volumetrically dominant high-SiO 2 (>74 wt%) rhyolite pumices record textural evidence for a significant disequilibrium event (partial dissolution ± resorption of cores and interiors) prior to the growth of 40–500 μm thick rim zones. This dissolution/regrowth history of orthopyroxene is recorded in the chemistry of co-crystallising amphiboles as a prominent inflection in the concentrations of Mn and Zn, two elements notably enriched in orthopyroxene relative to amphibole. Textural and chemical features, linked with in situ thermobarometric estimates, indicate that a major decompression event preceded the formation of the melt-dominant body. The decompression event is inferred to represent the extraction of large volumes of melt plus crystals from the Oruanui crystal mush/source zone at pressures of 140–300 MPa (~6–12 km depth). Orthopyroxene underwent partial dissolution during ascent before reestablishing in the melt-dominant magma body at pressures of 90–140 MPa (~3.5–6 km). We model Fe–Mg diffusion across the core-rim boundaries along the crystallographic a or b -axes to constrain the timing of this decompression event, which marked establishment of the melt-dominant magma body. Maximum modelled ages indicate that this event did not begin until ~1,600 years before eruption, consistent with constraints from zircon model-age spectra. Once extraction began, it underwent runaway acceleration with a peak extraction age of ~230 years, followed by an apparent period of stasis of ~60 years prior to eruption. The rapidity of the extraction and accumulation processes implies the involvement of a dynamic driving force which, in the rifted continental arc setting of the Taupo Volcanic Zone, seems likely to be represented by magma-assisted extensional tectonic processes.

We use comprehensive geochemical and petrological records from whole-rock samples, crystals, matrix glasses and melt inclusions to derive an integrated picture of the generation, accumulation and evacuation of 530 km 3 of crystal-poor rhyolite in the 25.4 ka Oruanui supereruption (New Zealand). New data from plagioclase, orthopyroxene, amphibole, quartz, Fe–Ti oxides, matrix glasses, and plagioclase- and quartz-hosted melt inclusions, in samples spanning different phases of the eruption, are integrated with existing data to build a history of the magma system prior to and during eruption. A thermally and compositionally zoned, parental crystal-rich (mush) body was developed during two periods of intensive crystallisation, 70 and 10–15 kyr before the eruption. The mush top was quartz-bearing and as shallow as ~3.5 km deep, and the roots quartz-free and extending to >10 km depth. Less than 600 year prior to the eruption, extraction of large volumes of ~840 °C low-silica rhyolite melt with some crystal cargo (between 1 and 10%), began from this mush to form a melt-dominant (eruptible) body that eventually extended from 3.5 to 6 km depth. Crystals from all levels of the mush were entrained into the eruptible magma, as seen in mineral zonation and amphibole model pressures. Rapid translation of crystals from the mush to the eruptible magma is reflected in textural and compositional diversity in crystal cores and melt inclusion compositions, versus uniformity in the outermost rims. Prior to eruption the assembled eruptible magma body was not thermally or compositionally zoned and at temperatures of ~790 °C, reflecting rapid cooling from the ~840 °C low-silica rhyolite feedstock magma. A subordinate but significant volume (3–5 km 3 ) of contrasting tholeiitic and calc-alkaline mafic material was co-erupted with the dominant rhyolite. These mafic clasts host crystals with compositions which demonstrate that there was some limited pre-eruptive physical interaction of mafic magmas with the mush and melt-dominant body. However, the mafic magmas do not appear to have triggered the eruption or controlled magmatic temperatures in the erupted rhyolite. Integration of textural and compositional data from all available crystal types, across all dominant and subordinate magmatic components, allow the history of the Oruanui magma body to be reconstructed over a wide range of temporal scales using multiple techniques. This history spans the tens of millennia required to grow the parental magma system (U–Th disequilibrium dating in zircon), through the centuries and decades required to assemble the eruptible magma body (textural and diffusion modelling in orthopyroxene), to the months, days, hours and minutes over which individual phases of the eruption occurred, identified through field observations tied to diffusion modelling in magnetite, olivine, quartz and feldspar. Tectonic processes, rather than any inherent characteristics of the magmatic system, were a principal factor acting to drive the rapid accumulation of magma and control its release episodically during the eruption. This work highlights the richness of information that can be gained by integrating multiple lines of petrologic evidence into a holistic timeline of field-verifiable processes.

Pyroclastic deposits from four caldera volcanoes in the Kermadec arc have been sampled from subaerial sections (Raoul and Macauley) and by dredging from the submerged volcano flanks (Macauley, Healy, and the newly discovered Raoul SW). Suites of 16–32 mm sized clasts have been analyzed for density and shape, and larger clasts have been analyzed for major element compositions. Density spectra for subaerial dry-type eruptions on Raoul Island have narrow unimodal distributions peaking at vesicularities of 80–85%, whereas ingress of external water (wet-type eruption) or extended timescales for degassing generate broader distributions, including denser clasts. Submarine-erupted pyroclasts show two different patterns. Healy and Raoul SW dredge samples and Macauley Island subaerial-emplaced samples are dominated by modes at ~80–85%, implying that submarine explosive volcanism at high eruption rates can generate clasts with similar vesicularities to their subaerial counterparts. A minor proportion of Healy and Raoul SW clasts also show a pink oxidation color, suggesting that hot clasts met air despite 0.5 to >1 km of intervening water. In contrast, Macauley dredged samples have a bimodal density spectrum dominated by clasts formed in a submarine-eruptive style that is not highly explosive. Macauley dredged pyroclasts are also the mixed products of multiple eruptions, as shown by pumice major-element chemistry, and the sea-floor deposits reflect complex volcanic and sedimentation histories. The Kermadec calderas are composite features, and wide dispersal of pumice does not require large single eruptions. When coupled with chemical constraints and textural observations, density spectra are useful for interpreting both eruptive style and the diversity of samples collected from the submarine environment.

Deposits of highly vesicular pumice that blanket submarine volcanoes are often attributed to explosive eruptions. Density and textural analysis of clasts dredged from the submarine Macauley Volcano, southwest Pacific Ocean, however, reveal an eruptive style that is neither explosive nor effusive, with clasts instead forming from buoyant detachment of a magma foam. Many submarine caldera volcanoes are blanketed with deposits of highly vesicular pumice, typically attributed to vigorous explosive activity 1 , 2 , 3 , 4 . However, it is challenging to relate volcanic products to specific eruptive styles in submarine volcanism 5 , 6 . Here we document vesicularity and textural characteristics of pumice clasts dredged from the submarine Macauley volcano in the Kermadec arc, southwest Pacific Ocean. We find that clasts show a bimodal distribution, with corresponding differences in vesicle abundances and shapes. Specifically, we find a sharp mode at 91% vesicularity and a broad mode at 65–80%. Subordinate clasts show gradients in vesicularity. We attribute the bimodality to a previously undocumented eruptive style that is neither effusive nor explosive. The eruption rate is insufficient to cause magma to fragment explosively, yet too high to passively feed a lava dome. Instead, the magma foam buoyantly detaches at the vent and rises as discrete magma parcels, or blebs, while continuing to vesiculate internally. The blebs are widely distributed by ocean currents before they disintegrate or become waterlogged. This disintegration creates individual clasts from interior and rim fragments, yielding the bimodal vesicularity characteristics. We conclude that the generation and widespread dispersal of highly vesicular pumice in the marine environment does not require highly explosive activity.

Important clues to the initiation and early behavior of large (super-) eruptions lie in the records of degassing during magma ascent. Here we investigate the timescales of magma ascent for three rhyolitic supereruptions that show field evidence for contrasting behavior at eruption onset: (1) 650 km3, 0.767 Ma Bishop Tuff, Long Valley; (2) 530 km3, 25.4 ka Oruanui eruption, Taupo; and (3) 2500 km3, 2.08 Ma Huckleberry Ridge Tuff, Yellowstone. During magma ascent, decompression causes volatile exsolution from the host melt into bubbles, leading to H2O and CO2 gradients in quartz-hosted re-entrants (REs; unsealed inclusions). These gradients are modeled to estimate ascent rates. We present best-fit modeled ascent rates for H2O and CO2 profiles for REs in early-erupted fall deposits from Bishop (n = 13), Oruanui (n = 9), and Huckleberry Ridge (n = 9). Using a Matlab script that includes an error minimization function, Bishop REs yield ascent rates of 0.6-13 m/s, overlapping with and extending beyond those of the Huckleberry Ridge (0.3-4.0 m/s). Re-entrants in Oruanui quartz crystals from the first two eruptive phases (of 10) yield the slowest ascent rates determined in this study (0.06-0.48 m/s), whereas those from phase three, which has clear field evidence for a marked increase in eruption intensity, are uniformly higher (1.4-2.6 m/s). For all three eruptions, the interiors of most REs appear to have re-equilibrated to lower H2O and CO2 concentrations when compared to co-erupted, enclosed melt inclusions in quartz. Such reequilibration implies the presence of an initial period of slower ascent, likely resulting from movement of magma from storage into a developing conduit system, prior to the faster (<1-2.5 h) final ascent of magma to the surface. This slower initial movement represents hours to several days of reequilibration, invalidating any assumption of constant decompression conditions from the storage region. However, the number of REs with deeper starting depths increases with stratigraphic height in all three deposits (particularly the Bishop Tuff), suggesting progressive elimination of the deep, sluggish ascent stage over time, which we interpret to be the result of maturing of the conduit system(s). Our results agree well with ascent rates estimated using theoretical approximations and numerical modeling for plinian rhyolitic eruptions (0.7-30 m/s), but overlap more with the slower estimates.

The transport and degassing pathways of volatiles through large silicic magmatic systems are central to understanding geothermal fluid compositions, ore deposit genesis, and volcanic eruption dynamics and impacts. Here, we document sulfur (S), chlorine (Cl), and fluorine (F) concentrations in a range of host materials in eruptive deposits from Taupō volcano (New Zealand). Materials analysed are groundmass glass, silicic melt inclusions, and microphenocrystic apatite that equilibrated in shallow melt-dominant magma bodies; silicic melt and apatite inclusions within crystal cores inferred to be sourced from deeper crystal mush; and olivine-hosted basaltic melt inclusions from mafic enclaves that represent the most primitive feedstock magmas. Sulfur and halogen concentrations each follow distinct concentration pathways during magma differentiation in response to changing pressures, temperatures, oxygen fugacities, crystallising mineral phases, the effects of volatile saturation, and the presence of an aqueous fluid phase. Sulfur contents in the basaltic melt inclusions (~ 2000 ppm) are typical for arc-type magmas, but drop to near detection limits by dacitic compositions, reflecting pyrrhotite crystallisation at ~ 60 wt. % SiO2 during the onset of magnetite crystallisation. In contrast, Cl increases from ~ 500 ppm in basalts to ~ 2500 ppm in dacitic compositions, due to incompatibility in the crystallising phases. Fluorine contents are similar between mafic and silicic compositions (< 1200 ppm) and are primarily controlled by the onset of apatite and/or amphibole crystallisation and then destabilisation. Sulfur and Cl partition strongly into an aqueous fluid and/or vapour phase in the shallow silicic system. Sulfur contents in the rhyolite melts are low, yet the Oruanui supereruption is associated with a major sulfate peak in ice core records in Antarctica and Greenland, implying that excess S was derived from a pre-eruptive gas phase, mafic magma recharge, and/or disintegration of a hydrothermal system. We estimate that the 25.5 ka Oruanui eruption ejected > 130 Tg of S (390 Tg sulfate) and up to ~ 1800 Tg of Cl, with potentially global impacts on climate and stratospheric ozone.

Late Pleistocene diatomaceous sediment was widely dispersed along with volcanic ash (tephra) across and beyond New Zealand by the 25.4 ka Oruanui supereruption from Taupo volcano. We present a detailed analysis of the diatom populations in the Oruanui tephra and the newly discovered floras in two other eruptions from the same volcano: the 28.6 ka Okaia and 1.8 ka Taupo eruptions. For comparison, the diatoms were also examined in Late Pleistocene and Holocene lake sediments from the Taupo Volcanic Zone (TVZ). Our study demonstrates how these microfossils provide insights into the lake history of the TVZ since the Last Glacial Maximum. Morphometric analysis of Aulacoseira valve dimensions provides a useful quantitative tool to distinguish environmental and eruptive processes within and between individual tephras. The Oruanui and Okaia diatom species and valve dimensions are highly consistent with a shared volcanic source, paleolake and eruption style (involving large-scale magma–water interaction). They are distinct from lacustrine sediments sourced elsewhere in the TVZ. Correspondence analysis shows that small, intact samples of erupted lake sediment (i.e., lithic clasts in ignimbrite) contain heterogeneous diatom populations, reflecting local variability in species composition of the paleolake and its shallowly buried sediments. Our analysis also shows a dramatic post-Oruanui supereruption decline in Cyclostephanos novaezelandiae, which likely reflects a combination of (1) reorganisation of the watershed in the aftermath of the eruption, and (2) overall climate warming following the Last Glacial Maximum. This decline is reflected in substantially lower proportions of C. novaezelandiae in the 1.8 ka Taupo eruption deposits, and even fewer in post-1.8 ka sediments from modern (Holocene) Lake Taupo. Our analysis highlights how the excellent preservation of siliceous microfossils in volcanic tephra may fingerprint the volcanic source region and retain a valuable record of volcanically influenced environmental change.