Search

New geochemical data, including Sr–Nd–Pb isotopes for whole-rock and groundmass samples, are reported for edifice-forming eruptives at Tongariro volcano, New Zealand, which span its ~ 350 ka to late Holocene history. Tongariro eruptives are medium-K basaltic-andesites to dacites (53.0–66.2 wt% SiO2) that evolved via assimilation-fractional crystallisation (AFC) processes partly or mostly in the uppermost 15 km of the crust. When ordered chronologically using a high-resolution 40Ar/39Ar-dated eruptive stratigraphy, the compositional data show systematic 10–130 kyr cycles. Mafic replenishment events inferred from MgO values occurred at ~ 230, ~ 151, ~ 88 and ~ 56 ka and in the late Holocene, with high-MgO flank vents erupting at ~ 160, ~ 117, ~ 35 and ~ 17.5 ka. Cycles in Sm/Nd, 87Sr/86Sr, 143Nd/144Nd and Pb isotopic ratios, which are decoupled from MgO, K2O and Rb/Sr cycles, indicate periods of prolonged crustal residence of magmas from ~ 230 to ~ 100 ka and ~ 95 to ~ 30 ka. AFC modelling shows that intermediate and silicic melt compositions, with r-values between 0.1 and 1, are needed to reproduce Tongariro compositional arrays. AFC models also indicate that ~ 20% of the average Tongariro magma comprises assimilated (meta)sedimentary basement material. Locally, Tongariro and adjacent Ruapehu volcanoes attain their most crust-like 87Sr/86Sr and 143Nd/144Nd compositions at ~ 100 and ~ 30 ka, paralleling with zircon model-age crystallisation modes at the rhyolitic Taupō volcano ~ 50 km to the NNE. These coincidences suggest that the timing and tempo of magma assembly processes at all three volcanoes were contemporaneous and may have been tectonically influenced since at least 200 ka.

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.

This work has been supported by the Marsden Fund grant VUW0813 (Royal Society of New Zealand to C.J.N.W.), a James Cook Fellowship (Royal Society of New Zealand) to C.J.N.W., and the ECLIPSE Programme, funded by the N.Z. Ministry of Business, Innovation and Employment. G.F.C. is supported by a NERC Standard Grant (NE/T000317/1), M.L.M. is supported by an NSF CAREER grant (EAR 2042662) and S.J.B. acknowledges Marsden Fund grant VUW1627.

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 ( 130 Tg of S (390 Tg sulfate) and up to ~ 1800 Tg of Cl, with potentially global impacts on climate and stratospheric ozone.

We consider here the validity, accuracy, and precision of rhyolite-MELTS modelling in inferring the pre-eruptive magma storage conditions for the caldera-forming 25.5 ka Oruanui and 232 CE (1.72 ka) Taupō eruptions at Taupō volcano, New Zealand as proposed by Pamukcu et al. (2020: Contrib Mineral Petrol 175: 48). There are four major issues with the modelling as used. First, the model is inappropriately applied for the Taupō event as this magma does not contain phenocrystic quartz. Second, the products of both eruptions, as is typical for virtually all Taupō Volcanic Zone rhyolites, contain plagioclase but lack sanidine, leading to increased model errors beyond those stated. Third, temperatures utilised for modelling crystallisation histories for both magmas are in conflict with independent measures of magmatic temperatures for these rocks. Fourth, glass compositions for each of these eruptions individually fall within analytical uncertainty yet yield model pressures that vary over hundreds of MPa, in part due to contrasting analytical protocols. We conclude that rhyolite-MELTS is too imprecise (errors of ± ~ 100 MPa or more for the methods applied) for quartz + 1 feldspar systems like at Taupō volcano specifically and the Taupō Volcanic Zone in general to be applied with confidence to deriving relative or absolute depths of magma storage.

The following statement was omitted from the acknowledgments section of the manuscript but is included now in accordance with U.S. Geological Survey policy. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government.

The youngest major caldera-forming event at Yellowstone was the ~ 630-ka eruption of the Lava Creek Tuff. The tuff as mapped consists of two major ignimbrite packages (members A and B), linked to widespread coeval fall deposits and formation of the Yellowstone Caldera. Subsequent activity included emplacement of numerous rhyolite flows and domes, and development of two structurally resurgent domes (Mallard Lake and Sour Creek) that accommodate strain due to continual uplift/subsidence cycles. Uplifted lithologies previously mapped on and adjacent to Sour Creek dome were thought to include the ~ 2.08-Ma Huckleberry Ridge Tuff, cropping out beneath Lava Creek Tuff members A and B. Mapped outcrops of this Huckleberry Ridge Tuff material were sampled as welded ignimbrite (sample YR345) on Sour Creek dome, and at nearby Bog Creek as welded ignimbrite (YR311) underlain by an indurated lithic lag breccia containing blocks of another welded ignimbrite (YR324). Zircon near-rim U–Pb analyses from these samples yield weighted mean ages of 661 ± 13 ka (YR345: 95% confidence), 655 ± 11 ka (YR311), and 664 ± 15 ka (YR324) (combined weighted mean of 658.8 ± 6.6 ka). We also studied two samples of ignimbrite previously mapped as Huckleberry Ridge Tuff on the northeastern perimeter of the Yellowstone Caldera, ~ 12 km ENE of Sour Creek dome. Sanidines from these samples yield 40Ar/39Ar age estimates of 634.5 ± 6.8 ka (8YC-358) and 630.9 ± 4.1 ka (8YC-359). These age data show that all these units represent previously unrecognized parts of the Lava Creek Tuff and do not have any relationship to the Huckleberry Ridge Tuff. Our observations and data imply that the Lava Creek eruption was more complex than is currently assumed, incorporating two tuff units additional to those currently mapped, and which themselves are separated by a time break sufficient for cooling and some reworking. The presence of a lag breccia suggests that a source vent lay nearby (

Complexities in the nature of large-scale silicic eruptions and their magmatic systems can be discerned through micro-analytical geochemical studies. We present high-resolution, stratigraphically constrained compositional data on glassy matrix material and feldspar crystals from the initial fall deposits and earliest ignimbrite (base of member A) of the 2.08 Ma, ~ 2500 km3 Huckleberry Ridge Tuff (HRT), Yellowstone. We use these data to document the nature of the magmatic system and compositional changes related to the transition from fall to widespread ignimbrite deposition, inferred to mark the onset of caldera collapse. Although major element glass compositions are relatively uniform, trace elements span a large range (e.g. Ba 10–900 ppm, Sr/Rb = 0.005–0.09), with highly evolved glasses dominating in the fall deposits. Several trace elements (e.g. Ba and light rare earth elements) in the glass samples serve to define statistically significant compositional clustering in the fall deposits and basal ignimbrite. These clusters are inferred to reflect melt compositions controlled by fractional crystallisation processes and are interpreted to represent multiple, discrete melt-dominant domains that were tapped by multiple vents. The onset of widespread ignimbrite deposition is marked by an increase in the number of erupted melt compositional clusters from four in the fall deposits to eight, representing nine melt-dominant domains. There is an absence of geographical variation of glass compositions within the basal ignimbrite, with samples from proximal to distal localities north, west and south of the HRT caldera exhibiting similar variability. Pairing of glass analyses with sanidine major and minor element compositional data suggests that the nine melt compositional domains converged at depth into two compositionally distinct upper-crustal magmatic lineages that were both active during these early stages of the eruption. Our data collectively indicate the evacuation of an exceptionally complex and heterogeneous magma system. The onset of widespread ignimbrite deposition, inferred to relate to caldera collapse, occurred after ~ 50 km3 of magma had been discharged. Although external controls were important as an eruption trigger, depressurisation of the system led to caldera collapse with the eruption of numerous discrete melt-dominant domains.

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.

Multiple, independent time markers are essential to correlate sediment and ice cores from the terrestrial, marine and glacial realms. These records constrain global paleoclimate reconstructions and inform future climate change scenarios. In the Northern Hemisphere, sub-visible layers of volcanic ash (cryptotephra) are valuable time markers due to their widespread dispersal and unique geochemical fingerprints. However, cryptotephra are not as widely identified in the Southern Hemisphere, leaving a gap in the climate record, particularly during the Last Glacial Maximum (LGM). Here we report the first identification of New Zealand volcanic ash in Antarctic ice. The Oruanui supereruption from Taupo volcano (25,580 ± 258 cal. a BP) provides a key time marker for the LGM in the New Zealand sector of the SW Pacific. This finding provides a high-precision chronological link to mid-latitude terrestrial and marine sites, and sheds light on the long-distance transport of tephra in the Southern Hemisphere. As occurred after identification of the Alaskan White River Ash in northern Europe, recognition of ash from the Oruanui eruption in Antarctica dramatically increases the reach and value of tephrochronology, providing links among climate records in widely different geographic areas and depositional environments.

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 km3 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 km3) 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.

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 km3 (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-SiO2 (>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.

Many submarine caldera volcanoes are blanketed with deposits of highly vesicular pumice, typically attributed to vigorous explosive activity1, 2, 3, 4. However, it is challenging to relate volcanic products to specific eruptive styles in submarine volcanism5, 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.

Trace concentrations of Ti in quartz are used to indicate the pressure and temperature conditions of crystallization in the ‘TitaniQ’ geothermobarometer of Thomas et al. (Contrib Miner Petrol 160:743–759, 2010). It utilises the partitioning of Ti into quartz as an indicator of the pressures and/or temperatures of crystal growth. For a given value of TiO2 activity in the system, if temperatures are inferred to ±20 °C, pressure is constrained to ±1 kbar and vice versa. There are significant contrasts, however, between the conclusions from TitaniQ and those for natural quartz (as well as other mineral phases) in volcanic rocks. Application of the TitaniQ model to quartz from the 27 ka Oruanui and 760 ka Bishop high-silica rhyolites, where the values of T, P and TiO2 activity are constrained by other means (Fe–Ti oxide equilibria, melt inclusion entrapment pressures in gas-saturated melts, melt and amphibole compositions), yields inconsistent results. If realistic values are given to any two of these three parameters, then the value of the third is wholly unrealistic. The model yields growth temperatures at or below the granite solidus, pressures in the lower crust or upper mantle, or TiO2 activities inconsistent with the mineralogical and chemical compositions of the magmas. CL imagery and measurements of Ti (and other elements) in quartz are of great value in showing the growth histories and changes in conditions experienced by crystals, but direct linkages to P, T conditions during crystal growth cannot be achieved.

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.

Cathodoluminescence (CL) zoning in quartz crystals from rhyolitic pumices in two ignimbrite members of the ~340-ka Whakamaru super-eruption deposits, Taupo Volcanic Zone, New Zealand, is investigated in conjunction with the analysis of Ti concentration in quartz to reconstruct the history of changing magma chamber conditions and to elucidate the eruption-triggering processes. CL intensity images are taken as a proxy for Ti concentration and thus temperature and/or pressure and/or compositional variations during crystal growth history. Estimates of the maximum temperature changes (i.e., assuming other factors influencing Ti uptake remain constant) are made using the TitaniQ geothermometer based on the Ti concentration in quartz. These results are reviewed in comparison with Fe–Ti oxide, feldspar-melt and amphibole geothermometry. Core-to-rim quartz Ti profiles record a marked change in conditions (temperature increase and/or pressure decrease and/or change in melt composition) causing and then following a significant resorption horizon in the outer parts of the crystals. Two alternative models that could explain the quartz Ti zonation invoke a temperature increase caused by mafic recharge and/or a pressure decrease due to magma ponding and re-equilibration at shallow crustal levels. Concomitant changes in melt composition and Ti activity may, however, also have strongly influenced Ti uptake into the quartz. Some crystals also show other marked increases in CL brightness internally, but any accompanying magmatic changes did not result in eruption. Diffusion modelling indicates that this significant change in conditions occurred over ~10–85 years prior to caldera-forming eruption. This rapid thermal pulse or pressure change is interpreted as evidence for open-system processes, and appears to record a magma chamber recharge event that rejuvenated the Whakamaru magma system (melt-dominant magma plus crystal mush), and potentially acted as a trigger for processes that led to eruption.

Rhyolitic pumices in the 26.5 ka Oruanui eruption (Taupo volcano, New Zealand) contain an average of 10 wt% crystals. About 2 wt% of the crystal population is feldspar crystals that display bluish–grey cloudy cores, the colour being imparted by exsolved needles of rutile. The volume of cloudy-cored feldspars thus amounts to ~1.0 km3 in a total magma volume of ~530 km3. The cored feldspars show great variability in detail, but in general have a rounded cloudy core bounded by a zone rich in glass and mineral inclusions, that was then overgrown by a euhedral clear rim. Sr-isotopic variations in eight representative crystals were measured on micromilled samples of selected growth zones in the cores and rims, and linked to feldspar compositions through microprobe traverses. The cloudy cores range from 87Sr/86Sr = 0.70547 to 0.70657, with compositions of An43 to An78. The overgrowth rims display wider variations: inner parts show extreme ranges in composition (maxima 87Sr/86Sr = 0.70764 and An78), while outer parts in seven of eight crystals are zoned, with outward-decreasing Sr-isotopic and An values to figures that are in accord with the bulk pumice and other, clear-feldspar values, respectively. The three parts of the crystals represent distinct regimes. The cloudy cores are inherited from an intermediate plutonic protolith that has been subjected to melting. The inner overgrowth rims were crystallised from a high temperature, relatively radiogenic melt derived from Mesozoic-Palaeozoic metasedimentary rocks (“greywacke”). The outer euhedral rims reflect mixing into and continued growth within the growing Oruanui magma body. The cloudy-cored feldspars also contain rare zircon inclusions. Twenty one zircons were recovered by HF digestion of a bulk sample of cloudy feldspars and analysed by SHRIMP for U–Th isotopes with which to calculate model ages. Eighteen of 21 crystals returned finite ages, the model-age spectrum of which is similar to the age spectra from free zircons in Oruanui pumices. Assembly of the Oruanui magma body was not only rapid (over ~40 kyr, as shown by other data) but involved a wide open system, with significant contributions from partly-melted intermediate-composition igneous intrusions (cloudy cores) and greywacke melts (inner overgrowths) being introduced into the magma body up to the point of eruption. Such open system behaviour contrasts with that proposed in models for comparably voluminous silicic magmas derived dominantly by fractionation (such as the Bishop Tuff) where the magma and its crystal cargo were better insulated thermally and chemically from country-rock interaction.

Crater Hill is a small volume alkali olivine basalt volcano in the Auckland volcanic field. Crater Hill consists of a sequence of pyroclastic and effusive eruptive units of which the earliest have low silica, ferromagnesian elements and Mg/Fe ratios, high incompatible elements and are more silica undersaturated while the last material to be erupted has higher silica, ferromagnesian elements and Mg/Fe ratios but relatively low incompatible elements. Through the sequence, Mg-number changes from 59 to 67 and LaN/LuN decreases by a factor of 3. This systematic compositional variation is interpreted to be the result of clinopyroxene ± spinel fractionation at pressures of at least 1.4–1.9 GPa, from a primary magma generated by small-degree partial melting in the garnet peridotite stability field (>2.5 GPa). Fractionation occurred where early crystals grew and accumulated along the conduit walls. The rising magma evolved along a polybaric liquid line of descent until it encountered lithosphere cold enough to chill the dike margin. Above this depth, further cooling resulted only in the growth of suspended phenocrysts in a magma separated from the country rock by a chilled margin. This process is observed in the Auckland volcanic field because the rate of magma production is very small allowing compositional features to be preserved that would be overwhelmed in a larger scale magmatic system.

The 1886 Plinian eruption of Tarawera, New Zealand, is a unique basaltic fissure-fed eruption with exceptionally well preserved fall deposits to within 200 meters of the source vents. These proximal deposits form a series of spatter/cinder half-cones along the northeastern 8-km-long segment of the 1886 fissure. Here we examine these deposits using grain size and clast componentry techniques. We contrast the products of the phreatomagmatic (phases I and III) and Plinian (Phase II) stages of the eruption and examine deposit variability as a function of contrasting eruptive intensity within the climactic phase (II) of the eruption. The opening phreatomagmatic phase I of the eruption involved gas-rich magma interacting with water and fragmenting at least 300 meters below the surface. The deposits of the climactic phase that followed have relatively uniform grain size but marked contrasts in the relative abundance of juvenile and wall rock (lithic) clasts. Deposits linked to vents associated with the high Plinian plume are more uniform than those characterized by a weaker cone-forming eruption style. During the third, and closing, phase of the eruption, magma withdrawal accompanied the onset of decoupling of the exsolved gas phase, leading to fragmentation at increasingly greater depths and significant wall rock collapse into the erupting vents. Variability in eruptive style during phase II along the fissure appears to be a function of shallow seated controls, in particular the variable extent of incorporation of lithic wall rock into the erupting jet, as a consequence of vent wall collapse. Widely dispersed beds centralized around Plinian sources along the fissure have very low lithic content; cone-forming beds at other craters that contain very high lithic contents. This incorporation led to a significant reduction of the velocity and stability of the jet at these latter steep-walled craters, and induced episodicity in the form of vent-clearing explosions. The result is a large reduction of the physical and thermal ability of these vents to contribute to a stable high eruptive plume. Instead large volumes of ejecta were sedimented prematurely from shallow heights at rates an order of magnitude greater than for historical Strombolian, Hawaiian and subPlinian eruptions. This study illustrates that sustained powerful Plinian eruptions can be accompanied by heterogeneities and instabilities of the eruptive jet. At Tarawera, the record of complex proximal transport and deposition processes in the eruptive jet cannot be inferred from the eruption products at distances greater than 400 m from the eruptive fissure. We suggest that study of ultraproximal deposits, as seen at Mt Tarawera, provides the only opportunity to document the complex, dynamic behavior of the jet region of Plinian eruptions.

Supereruptions (>1015 kg ≈ 450 km3 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 km3) 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 km3 and possibly up to 1000 km3 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.

Ice exerts a first-order control over the distribution and preservation of eruptive products on glaciated volcanoes. Defining the temporal and spatial distributions of ice-marginal lava flows provides valuable constraints on past glacial extents and is crucial for understanding the eruptive histories of such settings. Ice-marginal lava flows are well displayed on Ruapehu, a glaciated andesite-dacite composite cone in the southern Taupo Volcanic Zone, New Zealand. Flow morphology, fracture characteristics and 40Ar/39Ar geochronological data indicate that lavas erupted between ~51 and 15 ka interacted with large valley glaciers on Ruapehu. Ice-marginal lava flows exhibit grossly overthickened margins adjacent to glaciated valleys, are intercalated with glacial deposits, display fine-scale fracture networks indicative of chilling against ice, and are commonly ridge-capping due to their exclusion from valleys by glaciers. New and existing 40Ar/39Ar eruption ages for ice-marginal lava flows indicate that glaciers descending to 1300 m above sea level were present on Ruapehu between ~51–41 and ~27–15 ka. Younger lava flows located within valleys are characterised by blocky flow morphologies and fracture networks indicative of only localised and minor interaction with ice and/or snow, mainly in their upper reaches at elevations of ~2600–2400 m. An 40Ar/39Ar eruption age of 9 ± 3 ka (2σ error) determined for a valley-filling flow on the northern flank of Ruapehu indicates that glaciers had retreated to near-historical extents by the time of emplacement for this lava flow. The applicability of 40Ar/39Ar dating to ice-marginal flows on glaciated andesite-dacite composite volcanoes makes this technique an additional proxy for paleoclimate reconstructions.

The 1886 eruption of Tarawera, New Zealand, was unusual for a Plinian eruption because it involved entirely basaltic magma, originated in a 17-km-long fissure, and produced extremely overthickened proximal deposits with a complex geometry. This study focuses on an 8-km-long segment cutting across Mount Tarawera where over 50 point-source vents were active during the 5.5-h eruption. A detailed characterization of the proximal deposits is developed and used to interpret the range of styles and intensities of the vents, including changes with time. We identify the four vents that contributed most heavily to the Plinian fall and evaluate the extent to which current volcanic plume models are compatible with the depositional patterns at Tarawera. Three proximal units are mapped that have phreatomagmatic, magmatic, and phreatomagmatic characteristics, respectively. Within the magmatic proximal unit, beds of like character are grouped into packages and delineated on scaled cross sections. Package dispersal is quantified by measuring the linear thickness half-distance (t1/2) in the planes of the fissure walls. Most packages have localized dispersals (low t1/2), indicating that Strombolian-style activity dominated most vents. The more widely dispersed packages (high t1/2) reflect contributions from additional transport regimes that were more vigorous but still contributed considerable material to the proximal region. We conclude that the geometry of the observed proximal deposits requires three modes of fall transport: (1) fallout from the upper portions of the Plinian plumes produced principally by vents in four craters; (2) sedimentation from the margins of the lower portions of the Plinian plumes including the jets and possibly the lower convective regions; and (3) ejection by weak Strombolian-style explosions from vents that did not contribute significant volumes of particles to the high plume. We suggest that the curvature of the velocity profile across the jet region of each plume (1–4 km height) was important, and that the lower velocity at the margins allowed proximal deposition of a large volume of material with a wide grain-size range.

Large pyroclastic rhyolites are snapshots of evolving magma bodies, and preserved in their eruptive pyroclasts is a record of evolution up to the time of eruption. Here we focus on the conditions and processes in the Oruanui magma that erupted at 26.5 ka from Taupo Volcano, New Zealand. The 530 km3 (void-free) of material erupted in the Oruanui event is comparable in size to the Bishop Tuff in California, but differs in that rhyolitic pumice and glass compositions, although variable, did not change systematically with eruption order. We measured the concentrations of H2O, CO2 and major and trace elements in zoned phenocrysts and melt inclusions from individual pumice clasts covering the range from early to late erupted units. We also used cathodoluminescence imaging to infer growth histories of quartz phenocrysts. For quartz-hosted inclusions, we studied both fully enclosed melt inclusions and reentrants (connecting to host melt through a small opening). The textures and compositions of inclusions and phenocrysts reflect complex pre-eruptive processes of incomplete assimilation/partial melting, crystallization differentiation, magma mixing and gas saturation. ‘Restitic’ quartz occurs in seven of eight pumice clasts studied. Variations in dissolved H2O and CO2 in quartz-hosted melt inclusions reflect gas saturation in the Oruanui magma and crystallization depths of ∼3.5–7 km. Based on variations of dissolved H2O and CO2 in reentrants, the amount of exsolved gas at the beginning of eruption increased with depth, corresponding to decreasing density with depth. Pre-eruptive mixing of magma with varying gas content implies variations in magma bulk density that would have driven convective mixing.

Palaeomagnetic data from lithic clasts collected at 46 sites within layers 1 and 2 of the 1.8-ka Taupo ignimbrite, New Zealand, have been used to determine the palaeotemperatures and thermal structure of the deposit on its emplacement. Equilibrium temperatures from sites less than 30–40 km from vent are 150–300 °C, whereas at greater distances site equilibrium temperatures increase up to 400–500 °C. This variation is seen in both layer 1 and 2 deposits, with values for layer 1 being somewhat cooler, and with its increase in temperature occurring at a greater distance from vent. A temperature maximum at ~50 km from vent coincides with a zone of pink thermal-oxidation colouration of pumices previously inferred to reflect higher emplacement temperatures. Additional palaeomagnetic data collected by us and others from pumice clasts show comparable temperature variations, but these temperature estimates are shown here to be due to a chemical remanence and unreliable for accurate temperature estimates. Cooler temperatures in proximal parts of the ignimbrite are consistent with admixture of >20% cold lithic clasts at source and interaction with the pre-eruption Lake Taupo. The similar, but offset, increases in equilibrium temperatures for medial and distal layers 1 and 2 are consistent with both layers being deposited from the same flow. However, any proximal deposits left by the later, hotter material must have been subsequently eroded, or be so thin that our collection failed to sample them. Radial asymmetries in equilibrium temperatures as well as other physical parameters suggest that the deposit emplacement temperature is primarily determined at source, rather than by interaction with air during transport. These data support previous interpretations that a concentrated basal flow played a dominant role in emplacement and deposition of the Taupo ignimbrite.

Proximal (

Critical to understanding explosive eruptions is establishing how accurately representative pyroclasts are of processes during magma vesiculation and fragmentation. Here, we present data on densities, and vesicle size and number characteristics, for representative pyroclasts from six silicic eruptions of contrasting size and style from Raoul volcano (Kermadec arc). We use these data to evaluate histories of bubble nucleation, coalescence, and growth in explosive eruptions and to provide comparisons with pumiceous dome carapace material. Density/vesicularity distributions show a scarcity of pyroclasts with ∼65–75 % vesicularity; however, pyroclasts closest to this vesicularity range have the highest bubble number density (BND) values regardless of eruptive intensity or style. Clasts with vesicularities greater than this 65–75 % “pivotal” vesicularity range have decreasing BNDs with increasing vesicularities, interpreted to reflect continuing bubble growth and coalescence. Clasts with vesicularities less than the pivotal range have BNDs that decrease with decreasing vesicularity and preserve textures indicative of processes such as stalling and open system degassing prior to vesiculation in a microlite-rich magma, or vesiculation during slow ascent of degassing magma. Bubble size distributions (BSDs) and BNDs show variations consistent with 65–75 % representing the vesicularity at which vesiculating magma is most likely to undergo fragmentation, consistent with the closest packing of spheres. We consider that the observed vesicularity range may reflect the development of permeability in the magma through shearing as it flows through the conduit. These processes can act in concert with multiple nucleation events, generating a situation of heterogeneous bubble populations that permit some regions of the magma to expand and bubbles to coalesce with other regions in which permeable networks are formed. Fragmentation preserves the range in vesicularity seen as well as any post-fragmentation/pre-quenching expansion which may have occurred. We demonstrate that differing density pyroclasts from a single eruption interval can have widely varying BND values corresponding to the degree of bubble maturation that has occurred. The modal density clasts (the usual targets for vesicularity studies) have likely undergone some degree of bubble maturation and are therefore may not be representative of the magma at the onset of fragmentation.

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.

A series of alternating phreatomagmatic ("wet") and magmatic ("dry") basaltic pyroclastic deposits forming the Crater Hill tuff ring in New Zealand contains one unit (M1) which can only be interpreted as the products of mixing of ejecta from simultaneous wet and dry explosions at different portions of a multiple vent system. The principal characteristics of M1 are (a) rapid lateral changes in the thicknesses of, and proportions in juvenile components in individual beds, and (b) wide ranges of juvenile clast densities in every sample. M1 appears to have been associated with an elongate source of highly variable and fluctuating magma : water ratios and magma discharge rates. This contrasts with the only other documented mixed (wet and dry) basaltic pyroclastic deposits where mixing from two point sources of quite different but stable character has been inferred.

Two groups of poorly sorted ash-rich beds, previously interpreted as rain-flushed ashes, occur in the ca. AD 180 Hatepe Plinian pumice fall deposit at Taupo volcano, New Zealand. Two ash beds with similar dispersal patterns and an aggregate thickness of up to 13 cm make up the lowermost group (A). Group A beds extend 45 km north-east of the vent and cover 290 km2. In the southern part of the group A distribution area, a coarse ash to lapilli-size Plinian pumice bed (deposit B) separates the two group A beds. The scarcity of lapilli (material seen elsewhere from the still-depositing pumice fall) in group A beds indicates that they were rapidly transported and deposited. However, this rapid transportation and deposition did not produce cross-bedding, nor did it erode the underlying deposits. It is proposed that thick (>600 m) but dilute gravity currents generated from the collapsing outer margin of the otherwise buoyant Hatepe Plinian eruption column deposited the group A beds. The upper ash beds (group C) consist of one to seven layers, attain an aggregate thickness of ≤35 cm, and vary considerably in thickness and number of beds with respect to distance from vent. Group C beds contain variable amounts of ash mixed with angular Plinian pumices and are genuine rain-flushed ashes. Several recent eruptions at other volcanoes (Ukinrek Maars, Vulcan, Rabaul, La Soufrere de Guadeloupe and Soufriere, St Vincent) have produced gravity currents similar in style, but much smaller than those envisaged for group A deposits. The overloaded margins of otherwise buoyant eruption plumes generated these gravity currents. Laboratory studies have produced experimental gravity current analogues. Hazards from dilute gravity currents are considerable but often overlooked, thus the recognition of gravity current deposits will contribute to more thorough volcanic hazard assessment of prehistoric eruption sequences.

The 14 ka Puketarata eruption of Maroa caldera in Taupo Volcanic Zone was a dome-related event in which the bulk of the 0.25 km3 of eruption products were emplaced as phreatomagmatic fall and surge deposits. A rhyolitic dike encountered shallow groundwater during emplacement along a NE-trending normal fault, leading to shallow-seated explosions characterised by low to moderate water/magma ratios. The eruption products consist of two lava domes, a proximal tuff ring, three phreatic collapse craters, and a widespread fall deposit. The pyroclastic deposits contain dominantly dense juvenile clasts and few foreign lithics, and relate to very shallow-level disruption of the growing dome and its feeder dike with relatively little involvement of country rock. The distal fall deposit, representing 88% of the eruption products is, despite its uniform appearance and apparently subplinian dispersal, a composite feature equivalent to numerous discrete proximal phreatomagmatic lapilli fall layers, each deposited from a short-lived eruption column. The Puketarata products are subdivided into four units related to successive phases of:(A) shallow lava intrusion and initial dome growth; (B) rapid growth and destruction of dome lobes; (C) slower, sustained dome growth and restriction of explosive disruption to the dome margins; and (D) post-dome withdrawal of magma and crater-collapse. Phase D was phreatic, phases A and C had moderate water: magma ratios, and phase B a low water: magma ratio. Dome extrusion was most rapid during phase B, but so was destruction, and hence dome growth was largely accomplished during phase C. The Puketarata eruption illustrates how vent geometry and the presence of groundwater may control the style of silicic volcanism. Early activity was dominated by these external influences and sustained dome growth only followed after effective exclusion of external water from newly emplaced magma.

PUMiCE-rich pyroclastic flows leave deposits called ignimbrites1, which are common in the eruption records of arc and intracontinental volcanoes2,3. Pyroclastic flows represent the most dangerous manifestation of volcanism4 because of their volumes, rapid emplacement at high temperatures, and capacity to extend >100 km from the source. However, maximum distances travelled and extents of areas destroyed by pyroclastic flows remain poorly known. Historic examples have been relatively small, the largest reaching to only ∼35 km from the vent5, but prehistoric ignimbrites over 150 km across are known3,6–10. The largest ignimbrites are invariably partly buried or eroded, and only minimum estimates of their size (and hence of their destructive capacity) can be made. Here we describe the ∼l-Myr-old Kidnappers ignimbrite, erupted from the Taupo Volcanic Zone in New Zealand, which we have correlated for ≥385 km and which thus represents the most wide-spread ignimbrite yet known. We suggest that this wide distribution primarily reflects a coupling of high initial flow velocities with large mass fluxes, and that favourable flow paths and emplacement of flows over a wet substrate may have aided their mobility.

The vesicularity of juvenile clasts in pyroclastic deposits gives information on the relative timing of vesiculation and fragmentation, and on the role of magmatic volatiles versus external water in driving explosive eruptions. The vesicularity index and range are defined as the arithmetic mean and total spread of vesicularity values, respectively. Clast densities are measured for the 16–32 mm size fraction by water immersion techniques and converted to vesicularities using measured dense-rock equivalent densities. The techniques used are applied to four case studies involving magmas of widely varying viscosities and discharge rates: Kilauea Iki 1959 (basalt), Eifel tuff rings (basanite), Mayor Island cone-forming deposits (peralkaline rhyolite) and Taupo 1800 B.P. (calc-alkaline rhyolite). Previous theoretical studies suggested that a spectrum of clast vesicularities should be seen, depending on the magma viscosity, eruption rate, and the presence and timing of magma: water interaction. The new data are consistent with these predictions. In magmatic “dry” eruptions the vesicularity index lies uniformly in the range 70%–80% regardless of magma viscosity. For high viscosities and eruption rates the vesicularity ranges are narrow (< 25%), but broaden to between 30% and 50% as the viscosity and eruption rates are lowered and the volatiles and magma can de-couple. In phreatomagmatic “wet” eruptions, widely varying clast vesicularities reflect complex variations in the relative timing of vesiculation and water-induced fragmentation. Magma:water interaction at an early stage greatly reduces the vesicularity indices (< 40%) and broadens the ranges (as high as 80%), whereas late-stage interaction has only a minor effect on the index and broadens the range to a limited extent. Clast vesicularity represents a useful third parameter in addition to dispersal and fragmentation to characterise pyroclastic deposits.

The aureole materials that form an annulus of corrugated terrain surrounding Olympus Mons are considered to be the product of mass movement. The scarp at the mountain's foot formed as a result of this massive removal of material from the volcano's outer flanks. This interpretation is supported by a comparison of the amount of material originally available before scarp formation, and the present volume of aureole materials. On the basis of distribution, surface textures and theoretical considerations it is considered that the aureole was produced by a series of megaslides, rather than by a flow mechanism. Production of the megaslides may have been assisted by a period of widespread melting of permafrost.

The Taupo eruption was a complex volcanic event that represents the most recent activity at the Taupo Volcanic Centre1 in the North Island of New Zealand. The average2 of many radiocarbon ages dates the eruption at ∼AD 130. The records of ancient China and Rome refer to events in ∼AD 186 of the type which follow a major volcanic eruption, and we present here reasons for considering that these events were due to the Taupo eruption, and that the true date of the eruption is ∼AD 186.

Conventionally1, three important characteristics of ignimbrites are that they show a very pronounced tendency to pond in topographic depressions, possess a flat, horizontal or gently sloping upper depositional surface, and have a thickness generally between 10 and 1,000 m. The youngest major ignimbrite, that of 1912 in the Valley of Ten Thousand Smokes (Alaska)2–4 shows these characteristics well, being ponded in a valley 25 km long with a flat upper surface sloping down-valley at an average of 1.3°, and having an estimated thickness exceeding 100 m in places. AU three characteristics are commonly used as field criteria to distinguish ignimbrites from other pyroclastic rock bodies. Here we look at certain ignimbrites which depart significantly from this conventional form in that a major part of them occurs as a thin layer mantling the landscape, resting on slopes as steep as 30°, and with an upper surface sensibly parallel with the underlying surface. We have studied two examples, the 1,800-yr-old5 Taupo ignimbrite (New Zealand), and the 1,400-yr-old Rabual ignimbrite (New Britain)6,7. We also cite several others. These ignimbrites have a remarkably low aspect ratio. This ratio, previously applied to lava extrusions8,9, provides a useful means of quantifying the overall geometry of rock bodies. The ratio is that of average thickness to lateral spread; conveniently the latter is taken as the diameter of the circle which covers the same area as the rock body.