Oxygen isotope and trace element evidence for three-stage petrogenesis of the youngest episode (260–79 ka) of Yellowstone rhyolitic volcanism

We report the first high-precision δ¹⁸O analyses of glass, δ¹⁸O of minerals, and trace element concentrations in glass and minerals for the 260–79 ka Central Plateau Member (CPM) rhyolites of Yellowstone, a >350 km³ cumulative volume of lavas erupted inside of 630 ka Lava Creek Tuff (LCT) caldera. The glass analyses of these crystal-poor rhyolites provide direct characterization of the melt and its evolution through time. The δ¹⁸O_glass values are low and mostly homogeneous (4.5 ± 0.14 ‰) within and in between lavas that erupted in four different temporal episodes during 200 ka of CPM volcanism with a slight shift to lower δ¹⁸O in the youngest episode (Pitchstone Plateau). These values are lower than Yellowstone basalts (5.7–6 ‰), LCT (5.5 ‰), pre-, and extracaldera rhyolites (~7–8 ‰), but higher than the earliest 550–450 ka post-LCT rhyolites (1–2 ‰). The glass δ¹⁸O value is coupled with new clinopyroxene analyses and previously reported zircon analyses to calculate oxygen isotope equilibration temperatures. Clinopyroxene records >900 °C near-liquidus temperatures, while zircon records temperatures <850 °C similar to zircon saturation temperature estimates. Trace element concentrations in the same glass analyzed for oxygen isotopes show evidence for temporal decreases in Ti, Sr, Ba, and Eu—related to Fe–Ti oxide and sanidine (±quartz) crystallization control, while other trace elements remain similar or are enriched through time. The slight temporal increase in glass Zr concentrations may reflect similar or higher temperature magmas (via zircon saturation) through time, while previosuly reported temperature decreases (e.g., Ti-in-quartz) might reflect changing Ti concentrations with progressive melt evolution. Multiple analyses of glass across single samples and in profiles across lava flow surfaces document trace element heterogeneity with compatible behavior of all analyzed elements except Rb, Nb, and U. These new data provide evidence for a three-stage geochemical evolution of these most recent Yellowstone rhyolites: (1) repeated batch melting events at the base of a homogenized low-δ¹⁸O intracaldera fill resulting in liquidus rhyolite melt and a refractory residue that sequesters feldspar-compatible elements over time. This melting may be triggered by conductive "hot plate" heating by basaltic magma intruding beneath the Yellowstone caldera resulting in contact rhyolitic melt that crystallizes early clinopyroxene and/or sanidine at high temperature. (2) Heterogeneity within individual samples and across flows reflects crystallization of these melts during preeruptive storage of magma at at lower, zircon-saturated temperatures. Compatible behavior and variations of most trace elements within individual lava flows are the result of sanidine, quartz, Fe–Ti oxide, zircon, and chevkinite crystallization at this stage. (3) Internal mixing immediately prior to and/or during eruption disrupts, these compositional gradients in each parental magma body that are preserved as melt domains distributed throughout the lava flows. These results based on the most recent and best-preserved volcanic products from the Yellowstone volcanic system provide new insight into the multiple stages required to generate highly fractionated hot spot and rift-related rhyolites. Our proposed model differs from previous interpretations that extreme Sr and Ba depletion result from long-term crystallization of a single magma body—instead we suggest that punctuated batch melting events generated a sanidine-rich refractory residue and a melt source region progressively depleted in Sr and Ba.     

Magma mixing and the generation of isotopically juvenile silicic magma at Yellowstone caldera inferred from coupling 238U–230Th ages with trace elements and Hf and O isotopes in zircon and Pb isotopes in sanidine

The nature of compositional heterogeneity within large silicic magma bodies has important implications for how silicic reservoirs are assembled and evolve through time. We examine compositional heterogeneity in the youngest (~170 to 70 ka) post-caldera volcanism at Yellowstone caldera, the Central Plateau Member (CPM) rhyolites, as a case study. We compare ²³⁸U–²³⁰Th age, trace-element, and Hf isotopic data from zircons, and major-element, Ba, and Pb isotopic data from sanidines hosted in two CPM rhyolites (Hayden Valley and Solfatara Plateau flows) and one extracaldera rhyolite (Gibbon River flow), all of which erupted near the caldera margin ca. 100 ka. The Hayden Valley flow hosts two zircon populations and one sanidine population that are consistent with residence in the CPM reservoir. The Gibbon River flow hosts one zircon population that is compositionally distinct from Hayden Valley flow zircons. The Solfatara Plateau flow contains multiple sanidine populations and all three zircon populations found in the Hayden Valley and Gibbon River flows, demonstrating that the Solfatara Plateau flow formed by mixing extracaldera magma with the margin of the CPM reservoir. This process highlights the dynamic nature of magmatic interactions at the margins of large silicic reservoirs. More generally, Hf isotopic data from the CPM zircons provide the first direct evidence for isotopically juvenile magmas contributing mass to the youngest post-caldera magmatic system and demonstrate that the sources contributing magma to the CPM reservoir were heterogeneous in ¹⁷⁶Hf/¹⁷⁷Hf at ca. 100 ka. Thus, the limited compositional variability of CPM glasses reflects homogenization occurring within the CPM reservoir, not a homogeneous source.     

Genesis of the post-caldera eastern Upper Basin Member rhyolites, Yellowstone, WY: from volcanic stratigraphy, geochemistry, and radiogenic isotope modeling

An array of samples from the eastern Upper Basin Member of the Plateau Rhyolite (EUBM) in the Yellowstone Plateau, Wyoming, were collected and analyzed to evaluate styles of deposition, geochemical variation, and plausible sources for low δ¹⁸O rhyolites. Similar depositional styles and geochemistry suggest that the Tuff of Sulphur Creek and Tuff of Uncle Tom’s Trail were both deposited from pyroclastic density currents and are most likely part of the same unit. The middle unit of the EUBM, the Canyon flow, may be composed of multiple flows based on a wide range of Pb isotopic ratios (e.g., ²⁰⁶Pb/²⁰⁴Pb ranges from 17.54 to 17.86). The youngest EUBM, the Dunraven Road flow, appears to be a ring fracture dome and contains isotopic ratios and sparse phenocrysts that are similar to extra-caldera rhyolites of the younger Roaring Mountain Member. Petrologic textures, more radiogenic ⁸⁷Sr/⁸⁶Sr in plagioclase phenocrysts (0.7134–0.7185) than groundmass and whole-rock ratios (0.7099–0.7161), and δ¹⁸O depletions on the order of 5‰ found in the Tuff of Sulphur Creek and Canyon flow indicate at least a two-stage petrogenesis involving an initial source rock formed by assimilation and fractional crystallization processes, which cooled and was hydrothermally altered. The source rock was then lowered to melting depth by caldera collapse and remelted and erupted. The presence of a low δ¹⁸O extra-caldera rhyolite indicates that country rock may have been hydrothermally altered at depth and then assimilated to form the Dunraven Road flow.     

Rapid extraction of discrete magma batches from a large differentiating magma chamber: the Central Plateau Member rhyolites, Yellowstone Caldera, Wyoming

The Central Plateau Member rhyolites have been erupted between 173 and 70 ka and are the youngest Yellowstone intracaldera rhyolites. They mostly comprise very voluminous lava flows totaling ~600 km³ in volume. Their eruptive vents define two NNW-trending lineaments which are aligned with regional faults. We present new whole rock, glass, and mineral analyses and propose a petrogenetic and volcano-tectonic model for these rhyolites. At a caldera-wide scale, there is a temporal enrichment in elements such as Nb, Y and HREE, and a depletion in Sr, Ba, and Ce/Yb. Simultaneously, clinopyroxene becomes less magnesian while Ti contents in quartz decrease. By contrast, quartz in all rhyolites is rounded and bears long glass re-entrants, suggesting heating. Based on these data and observations, we propose that the Central Plateau Member rhyolites have been generated as follows. A hydrothermally altered low-δ¹⁸O rhyolitic protolith beneath the Mallard Lake Resurgent Dome in the southwestern part of the caldera started to melt at ~250 ka. Repeated heating pulses caused the melting front to expand radially, and a large crystal mush formed beneath much of the caldera. The mush was able to differentiate but not erupt due to its high crystallinity and viscosity. Further inputs of heat and silicic magma in this mush increased the degree of melting, forming crystal-poor magma batches which erupted a few hundred to a few thousand years later through regional faults to form the Central Plateau Member rhyolites.     

Initiation of large-volume silicic centers in the Yellowstone hotspot track: insights from H<Subscript>2</Subscript>O- and F-rich quartz-hosted rhyolitic melt inclusions in the Arbon Valley Tuff of the Snake River Plain

During the onset of caldera cluster volcanism at a new location in the Snake River Plain (SRP), there is an increase in basalt fluxing into the crust and diverse silicic volcanic products are generated. The SRP contains abundant and compositionally diverse hot, dry, and often low-δ¹⁸O silicic volcanic rocks produced through time during the formation of individual caldera clusters, but more H₂O-rich eruptive products are rare. We report analyses of quartz-hosted melt inclusions from pumice clasts from the upper and lower Arbon Valley Tuff (AVT) to gain insight into the initiation of caldera cluster volcanism. The AVT, a voluminous, caldera-forming rhyolite, represents the commencement of volcanism (10.44 Ma) at the Picabo volcanic field of the Yellowstone hotspot track. This is a normal δ¹⁸O rhyolite consisting of early and late erupted members (lower and upper AVT, respectively) with extremely radiogenic Sr isotopes and unradiogenic Nd isotopes, requiring that ~50 % of the mass of these elements is derived from melts of Archean upper crust. Our data reveal distinctive features of the early erupted lower AVT melt including: variable F concentrations up to 1.4 wt%, homogenous and low Cl concentrations (~0.08 wt%), H₂O contents ranging from 2.3 to 6.4 wt%, CO₂ contents ranging from 79 to 410 ppm, and enrichment of incompatible elements compared to the late erupted AVT, subsequent Picabo rhyolites, SRP rhyolites, and melt inclusions from other metaluminous rhyolites (e.g., Bishop Tuff, Mesa Falls Tuff). We couple melt inclusion data with Ti measurements and cathodoluminescence (CL) imaging of the host quartz phenocrysts to elucidate the petrogenetic evolution of the AVT rhyolitic magma. We observe complex and multistage CL zoning patterns, the most critical being multiple truncations indicative of several dissolution–reprecipitation episodes with bright CL cores (higher Ti) and occasional bright CL rims (higher Ti). We interpret the high H₂O, F, F/Cl, and incompatible trace element concentrations in the context of a model involving melting of Archean crust and mixing of the crustal melt with basaltic differentiates, followed by multiple stages of fractional crystallization, remelting, and melt extraction. This multistage process, which we refer to as distillation, is further supported by the complex CL zoning patterns in quartz. We interpret new Δ¹⁸O(Qz-Mt) isotope measurements, demonstrating a 0.4 ‰ or ~180 °C temperature difference, and strong Sr isotopic and chemical differences between the upper and lower AVT to represent two separate eruptions. Similarities between the AVT and the first caldera-forming eruptions of other caldera clusters in the SRP (Yellowstone, Heise and Bruneau Jarbidge) suggest that the more evolved, lower-temperature, more H₂O-rich rhyolites of the SRP are important in the initiation of a caldera cluster during the onset of plume impingement.     

Geochronological and isotopic records of crustal storage and assimilation in the Wolverine Creek–Conant Creek system,Heise eruptive centre,Snake River Plain

Understanding the processes of differentiation of the Yellowstone–Snake River Plain (YSRP) rhyolites is typically impeded by the apparent lack of erupted intermediate compositions as well as the complex nature of their shallow interaction with the surrounding crust responsible for their typically low O isotopic ratios. A pair of normal-δ¹⁸O rhyolitic eruptions from the Heise eruptive centre in eastern Idaho, the Wolverine Creek Tuff and the Conant Creek Tuff, represent unique magmatic products of the Yellowstone hotspot preserving abundant vestiges of the intermediate differentiation steps leading to rhyolite generation. We address both shallow and deep processes of magma generation and storage in the two units by combining high-precision ID–TIMS U–Pb zircon geochronology, trace element, O and Hf isotopic studies of zircon, and Sr isotopic analyses of individual high-Mg# pyroxenes inherited from lower- to mid-crustal differentiation stages. The zircon geochronology confirms the derivation of both tuffs from the same rhyolitic magma reservoir erupted at 5.5941 ± 0.0097 Ma, preceded by at least 92 ± 14 ky of continuous or intermittent zircon saturation approximating the length of pre-eruptive magma accumulation in the upper crust. Some low-Mg# pyroxenes enclosing zircons predate the eruption by at least 45 ± 27 ky, illustrating the co-crystallisation of major and accessory phases in the near-liquidus rhyolitic melts of the YSRP over a significant period of time. Coeval zircon crystals are isotopically heterogeneous (two populations at εHf ~?5 and ?13), requiring the assembly of isotopically distinct melt pockets directly prior to, or during, the eruption. The primitive Mg# 60–90 pyroxenes are out of isotopic equilibrium with the host rhyolitic melt (⁸⁷Sr/⁸⁶Sr_i = 0.70889), covering a range of ⁸⁷Sr/⁸⁶Sr_i = 0.70705–0.70883 corresponding to ratios typical of the most radiogenic YSRP basalts to the least radiogenic YSRP rhyolites. Together with the low εHf in zircon, the Sr isotopic ratios illustrate limited assimilation dominated by radiogenic Archean crustal source materials incorporated into variably evolved YSRP melts as they progress towards rhyolitic compositions by assimilation–fractional crystallisation.     

Evolution of silicic magmas in the Kos-Nisyros volcanic center,Greece: a petrological cycle associated with caldera collapse

Multiple eruptions of silicic magma (dacite and rhyolites) occurred over the last ~3 My in the Kos-Nisyros volcanic center (eastern Aegean sea). During this period, magmas have changed from hornblende-biotite-rich units with low eruption temperatures (≤750–800°C; Kefalos and Kos dacites and rhyolites) to hotter, pyroxene-bearing units (>800–850°C; Nisyros rhyodacites) and are transitioning back to cooler magmas (Yali rhyolites). New whole-rock compositions, mineral chemistry, and zircon Hf isotopes show that these three types of silicic magmas followed the same differentiation trend: they all evolved by crystal fractionation and minor crustal assimilation (AFC) from parents with intermediate compositions characterized by high Sr/Y and low Nb content, following a wet, high oxygen fugacity liquid line of descent typical of subduction zones. As the transition between the Kos-Kefalos and Nisyros-type magmas occurred immediately and abruptly after the major caldera collapse in the area (the 161 ka Kos Plateau Tuff; KPT), we suggest that the efficient emptying of the magma chamber during the KPT drew out most of the eruptible, volatile-charged magma and partly solidified the unerupted mush zone in the upper crust due to rapid unloading, decompression, and coincident crystallization. Subsequently, the system reestablished a shallow silicic production zone from more mafic parents, recharged from the mid to lower crust. The first silicic eruptions evolving from these parents after the caldera collapse (Nisyros units) were hotter (up to >100°C) than the caldera-forming event and erupted from reservoirs characterized by different mineral proportions (more plagioclase and less amphibole). We interpret such a change as a reflection of slightly drier conditions in the magmatic column after the caldera collapse due to the decompression event. With time, the upper crustal intermediate mush progressively transitioned into the cold-wet state that prevailed during the Kefalos-Kos stage. The recent eruptions of the high-SiO₂ rhyolite on Yali Island, which are low temperature and hydrous phases (sanidine, quartz, biotite), suggest that another large, potentially explosive magma chamber is presently building under the Kos-Nisyros volcanic center.     

U–Pb zircon geochronology of silicic tuffs from the Timber Mountain/Oasis Valley caldera complex, Nevada: rapid generation of large volume magmas by shallow-level remelting

Large volumes of silicic magma were produced on a very short timescale in the nested caldera complex of the SW Nevada volcanic field (SWNVF). Voluminous ash flows erupted in two paired events: Topopah Spring (TS, >1,200 km³, 12.8 Ma)–Tiva Canyon (TC, 1,000 km³, 12.7 Ma) and Rainier Mesa (RM, 1,200 km³, 11.6 Ma)–Ammonia Tanks (AT, 900 km³, 11.45 Ma; all cited ages are previously published ⁴⁰Ar/³⁹Ar sanidine ages). Within each pair, eruptions are separated by only 0.1–0.15 My and produced tuffs with contrasting isotopic values. These events represent nearly complete evacuation of sheet-like magma chambers formed in the extensional Basin and Range environment. We present ion microprobe ages from zircons in the zoned ash-flow sheets of TS, TC, RM, and AT in conjunction with δ¹⁸O values of zircons and other phenocrysts, which differ dramatically among subsequently erupted units. Bulk zircons in the low-δ¹⁸O AT cycle were earlier determined to exhibit ∼1.5‰ core-to-rim oxygen isotope zoning; and high-spatial resolution zircon analyses by ion microprobe reveal the presence of older grains that are zoned by 0.5–2.5‰. The following U–Pb isochron ages were calculated after correcting for the initial U–Pb disequilibria: AT (zircon rims: 11.7 ± 0.2 Ma; cores: 12.0 ± 0.1 Ma); pre-AT rhyolite lava: (12.0 ± 0.3 Ma); RM: 12.4 ± 0.3); TC: (13.2 ± 0.15 Ma); TS: (13.5 ± 0.2). Average zircon crystallization ages calculated from weighted regression or cumulative averaging are older than the Ar–Ar stratigraphy, but preserve the comparably short time gaps within each of two major eruption cycles (TS/TC, RM/AT). Notably, every sample yields average zircon ages that are 0.70–0.35 Ma older than the respective Ar–Ar eruption ages. The Th/U ratio of SWNVF zircons are 0.4–4.7, higher than typically found in igneous zircons, which correlates with elevated Th/U of the whole rocks (5–16). High Th/U could be explained if uranium was preferentially removed by hydrothermal solutions or is retained in the protolith during partial melting. For low-δ¹⁸O AT-cycle magmas, rim ages from unpolished zircons overlap within analytical uncertainties with the ⁴⁰Ar/³⁹Ar eruption age compared to core ages that are on average ∼0.2–0.3 My older than even the age of the preceding caldera forming eruption of RM tuff. This age difference, the core-to-rim oxygen isotope zoning in AT zircons, and disequilibrium quartz–zircon and melt-zircon isotopic fractionations suggest that AT magma recycled older zircons derived from the RM and older eruptive cycles. These results suggest that the low-δ¹⁸O AT magmas were generated by melting a hydrothermally-altered protolith from the same nested complex that erupted high-δ¹⁸O magmas of the RM cycle only 0.15 My prior to the eruption of the AT, the largest volume low-δ¹⁸O magma presently known.Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Accessory mineral U–Th–Pb ages and <Superscript>40</Superscript>Ar/<Superscript>39</Superscript>Ar eruption chronology,and their bearing on rhyolitic magma evolution in the Pleistocene Coso volcanic field,California

We determined Ar/Ar eruption ages of eight extrusions from the Pleistocene Coso volcanic field, a long-lived series of small volume rhyolitic domes in eastern California. Combined with ion-microprobe dating of crystal ages of zircon and allanite from these lavas and from granophyre geothermal well cuttings, we were able to track the range of magma-production rates over the past 650 ka at Coso. In ≤230 ka rhyolites we find no evidence of protracted magma residence or recycled zircon (or allanite) from Pleistocene predecessors. A significant subset of zircon in the ~85 ka rhyolites yielded ages between ~100 and 200 Ma, requiring that generation of at least some rhyolites involves material from Mesozoic basement. Similar zircon xenocrysts are found in an ~200 ka granophyre. The new age constraints imply that magma evolution at Coso can occur rapidly as demonstrated by significant changes in rhyolite composition over short time intervals (≤10’s to 100’s ka). In conjunction with radioisotopic age constraints from other young silicic volcanic fields, dating of Coso rhyolites highlights the fact that at least some (and often the more voluminous) rhyolites are produced relatively rapidly, but that many small-volume rhyolites likely represent separation from long-lived mushy magma bodies.     

Origin and Evolution of Silicic Magmatism at Yellowstone Based on Ion Microprobe Analysis of Isotopically Zoned Zircons

The origin of large-volume Yellowstone ignimbrites and smaller-volumeintra-caldera lavas requires shallow remelting of enormous volumesof variably ¹⁸O-depleted volcanic and sub-volcanic rocks alteredby hydrothermal activity. Zircons provide probes of these processesas they preserve older ages and inherited ¹⁸O values. This studypresents a high-resolution, oxygen isotope examination of volcanismat Yellowstone using ion microprobe analysis with an averageprecision of ± 0·2 and a 10 µm spot size.We report 357 analyses of cores and rims of zircons, and isotopeprofiles of 142 single zircons in 11 units that represent majorYellowstone ignimbrites, and post-caldera lavas. Many zirconsfrom these samples were previously dated in the same spots bysensitive high-resolution ion microprobe (SHRIMP), and all zirconswere analyzed for oxygen isotope ratios in bulk as a functionof grain size by laser fluorination. We additionally reportoxygen isotope analyses of quartz crystals in three units. Theresults of this work provide the following new observations.(1) Most zircons from post-caldera low-¹⁸O lavas are zoned,with higher ¹⁸O values and highly variable U–Pb ages inthe cores that suggest inheritance from pre-caldera rocks exposedon the surface. (2) Many of the higher-¹⁸O zircon cores in theselavas have U–Pb zircon crystallization ages that postdatecaldera formation, but pre-date the eruption age by 10–20kyr, and represent inheritance of unexposed post-caldera sub-volcanicunits that have ¹⁸O similar to the Lava Creek Tuff. (3) Youngand voluminous 0·25–0·1 Ma intra-calderalavas, which represent the latest volcanic activity at Yellowstone,contain zircons with both high-¹⁸O and low-¹⁸O cores surroundedby an intermediate-¹⁸O rim. This implies inheritance of a varietyof rocks from high-¹⁸O pre-caldera and low-¹⁸O post-calderaunits, followed by residence in a common intermediate-¹⁸O meltprior to eruption. (4) Major ignimbrites of Huckleberry Ridge,and to a lesser extent the Lava Creek and Mesa Falls Tuffs,contain zoned zircons with lower-¹⁸O zircon cores, suggestingthat melting and zircon inheritance from the low-¹⁸O hydrothermallyaltered carapace was an important process during formation ofthese large magma bodies prior to caldera collapse. (5) The¹⁸O zoning in the majority of zircon core–rim interfacesis step-like rather than smoothly inflected, suggesting thatprocesses of solution–reprecipitation were more importantthan intra-crystalline oxygen diffusion. Concave-downward zirconcrystal size distributions support dissolution of the smallercrystals and growth of rims on larger crystals. This study suggeststhat silicic magmatism at Yellowstone proceeded via rapid, shallow-levelremelting of earlier erupted and hydrothermally altered Yellowstonesource rocks and that pulses of basaltic magma provided theheat for melting. Each post-caldera Yellowstone lava representsan independent homogenized magma batch that was generated rapidlyby remelting of source rocks of various ages and ¹⁸O values.The commonly held model of a single, large-volume, super-solidus,mushy-state magma chamber that is periodically reactivated andproduces rhyolitic offspring is not supported by our data. Rather,the source rocks for the Yellowstone volcanism were cooled belowthe solidus, hydrothermally altered by heated meteoric watersthat caused low ¹⁸O values, and then remelted in distinct pocketsby intrusion of basic magmas. Each packet of new melt inheritedzircons that retained older age and ¹⁸O values. This interpretationmay have significance for interpreting seismic data for crustallow-velocity zones in which magma mush and solidified areasexperiencing hydrothermal circulation occur side by side. Newbasalt intrusions into this solidifying batholith are requiredto form the youngest volcanic rocks that erupted as independentrhyolitic magmas. We also suggest that the Lava Creek Tuff magmawas already an uneruptable mush by the time of the first post-calderaeruption after 0·1 Myr of the climactic caldera-formingeruption. KEY WORDS: Yellowstone; oxygen isotopes; geochronology; isotope zoning; zircon; U–Pb dating; caldera; rhyolite; ion microprobe     

Origins and evolution of rhyolitic magmas in the central Snake River Plain: insights from coupled high-precision geochronology,oxygen isotope,and hafnium isotope analyses of zircon

We present new high-precision CA-ID-TIMS and in situ U–Pb ages together with Hf and O isotopic analyses (analyses performed all on the same grains) from four tuffs from the 15?10 Ma Bruneau–Jarbidge center of the Snake River Plain and from three rhyolitic units from the Kimberly borehole in the neighboring 10?6 Ma Twin Falls volcanic center. We find significant intrasample diversity in zircon ages (ranges of up to 3 Myr) and in δ¹⁸O (ranges of up to 6‰) and ε_Hf (ranges of up to 24 ε units) values. Zircon rims are also more homogeneous than the associated cores, and we show that zircon rim growth occurs faster than the resolution of in situ dating techniques. CA-ID-TIMS dating of a subset of zircon grains from the Twin Falls samples reveals complex crystallization histories spanning 10⁴–10⁶ years prior to some eruptions, suggesting that magma genesis was characterized by the cyclic remelting of buried volcanic rocks and intrusions associated with previous magmatic episodes. Age-dependent trends in zircon isotopic compositions show that rhyolite production in the Yellowstone hotspot track is driven by the mixing of mantle-derived melts (normal δ¹⁸O and ε_Hf) and a combination of Precambrian basement rock (normal δ¹⁸O and ε_Hf down to ??60) and shallow Mesozoic and Cenozoic age rocks, some of which are hydrothermally altered (to low δ¹⁸O values) by earlier stages of Snake River Plain magmatism. These crustal melts hybridize with juvenile basalts and rhyolites to produce the erupted rhyolites. We also observe that the Precambrian basement rock is only an important component in the erupted magmas in the first eruption at each caldera center, suggesting that the accumulation of new intrusions quickly builds an upper crustal intrusive body which is isolated from the Precambrian basement and evolves towards more isotopically juvenile and lower-δ¹⁸O compositions over time.     

The dark side of zircon: textural,age, oxygen isotopic and trace element evidence of fluid saturation in the subvolcanic reservoir of the Island Park-Mount Jackson Rhyolite,Yellowstone (USA)

The Island Park-Mount Jackson series in the Yellowstone volcanic field, Wyoming (USA), is a suite of rhyolitic domes and lavas that erupted between the caldera-forming eruptions of the Mesa Fall Tuff (1.3 Ma) and the Lava Creek Tuff (0.6 Ma). Combined zircon U/Pb geochronology, Raman spectroscopy, oxygen isotopic and trace elemental compositions document storage conditions of these magmas between consecutive supereruptions. Based on comparison with co-erupted melt compositions and textural criteria, four zircon compositional groups are identified that record different stages along a continuous magmatic evolution from trace element-poor rhyolite at high temperatures to extremely fractionated rhyolite where zircon trace elements are highly enriched (e.g., >?1000 ppm U). These latter zircon domains are dark in cathodoluminescence images and show broadened Raman peaks relative to near-endmember zircon, indicating that substitution of non-stoichiometric trace elements into zircon leads to distortion of the crystal lattice. Some of these zircon domains contain inclusions of U-Th-REE-phases, likely originating from coupled dissolution–reprecipitation of metastable trace element-rich zircon in the presence of a fluid phase. Rhyolite-MELTS simulations indicate that at the conditions required to produce the observed enrichment in trace elements, a fluid phase is likely present. These findings illustrate that zircons can be assembled from a variety of co-existing magmatic environments in the same magma reservoir, including near-solidus volatile-rich melts close to the magmatic–hydrothermal transition.     

Do S-type granites commonly sample infracrustal sources? New results from an integrated O, U–Pb and Hf isotope study of zircon

In contrast to I-type granites, which commonly comprise infracrustal and supracrustal sources, S-type granites typically incorporate predominantly supracrustal sources. The initial aim of this study was to identify the sources of three Scottish Caledonian (~460 Ma) S-type granites (Kemnay, Cove and Nigg Bay) by conducting oxygen, U–Pb and Hf isotope analyses in zircon in order to characterise one potential end-member magma involved in the genesis of the voluminous late Caledonian (~430–400 Ma) I-type granites. Field, whole-rock geochemical and isotopic data are consistent with the generation of the S-type granites by melting their Dalradian Supergroup country rocks. While Hf isotope compositions of magmatic zircon, U–Pb data of inherited zircons, and high mean zircon δ¹⁸O values of 9.0 ± 2.7‰ (2SD) and 9.8 ± 2.0‰ for the Kemnay and Cove granites support this model, the Nigg Bay Granite contains zircons with much lower δ¹⁸O values (6.8 ± 2.1‰), similar to those found in Scottish I-type granites. This suggests that the Nigg Bay Granite contains low-δ¹⁸O material representing either altered supracrustal material, or more likely, an infracrustal source component with mantle-like δ¹⁸O. Mixing trends in plots of δ¹⁸O vs. εHf for S-type granite zircons indicate involvement of at least two sources in all three granites. This pilot study of Scottish Caledonian S-type granites demonstrates that, while field and whole-rock geochemical data are consistent with local melting of only supracrustal sources, the oxygen isotopic record stored in zircon reveals a much more complex petrogenetic evolution involving two or more magma sources.     

Oxygen isotope study of the Long Valley magma system, California: isotope thermometry and convection in large silicic magma bodies

Products of voluminous pyroclastic eruptions with eruptive draw-down of several kilometers provide a snap-shot view of batholith-scale magma chambers, and quench pre-eruptive isotopic fractionations (i.e., temperatures) between minerals. We report analyses of oxygen isotope ratio in individual quartz phenocrysts and concentrates of magnetite, pyroxene, and zircon from individual pumice clasts of ignimbrite and fall units of caldera-forming 0.76 Ma Bishop Tuff (BT), pre-caldera Glass Mountain (2.1-0.78 Ma), and post-caldera rhyolites (0.65-0.04 Ma) to characterize the long-lived, batholith-scale magma chamber beneath Long Valley Caldera in California. Values of '¹⁸O show a subtle 1‰ decrease from the oldest Glass Mountain lavas to the youngest post-caldera rhyolites. Older Glass Mountain lavas exhibit larger (~1‰) variability of '¹⁸O(quartz). The youngest domes of Glass Mountain are similar to BT in '¹⁸O(quartz) values and reflect convective homogenization during formation of BT magma chamber surrounded by extremely heterogeneous country rocks (ranging from 2 to +29‰). Oxygen isotope thermometry of BT confirms a temperature gradient between "Late" (815 °C) and "Early" (715 °C) BT. The '¹⁸O(quartz) values of "Early" and "Late" BT are +8.33 and 8.21‰, consistent with a constant '¹⁸O(melt)=7.8ǂ.1‰ and 100 °C temperature difference. Zircon-melt saturation equilibria gives a similar temperature range. Values of '¹⁸O(quartz) for different stratigraphic units of BT, and in pumice clasts ranging in pre-eruptive depths from 6 to 11 km (based on melt inclusions), and document vertical and lateral homogeneity of '¹⁸O(melt). Worldwide, five other large-volume rhyolites, Lava Creek, Lower Bandelier, Fish Canyon, Cerro Galan, and Toba, exhibit equal '¹⁸O(melt) values of earlier and later erupted portions in each of the these climactic caldera-forming eruptions. We interpret the large-scale '¹⁸O homogeneity of BT and other large magma chambers as evidence of their longevity (>10⁵ years) and convection. However, remaining isotopic zoning in some quartz phenocrysts, trace element gradients in feldspars, and quartz and zircon crystal size distributions are more consistent with far shorter timescales (10²-10⁴ years). We propose a sidewall-crystallization model that promotes convective homogenization, roofward accumulation of more evolved and stagnant, volatile-rich liquid, and develops compositional and temperature gradients in pre-climactic magma chamber. Crystal + melt + gas bubbles mush near chamber walls of variable '¹⁸O gets periodically remobilized in response to chamber refill by new hotter magmas. One such episode of chamber refill by high-Ti, Sr, Ba, Zr, and volatile-richer magma happened 10³-10⁴ years prior to the 0.76-Ma caldera collapse that caused magma mixing at the base, mush thawing near the roof and walls, and downward settling of phenocrysts into this hybrid melt.     

Decoupled crystallization and eruption histories of the rhyolite magmatic system at Tarawera volcano revealed by zircon ages and growth rates

We obtained U–Th disequilibrium age data on zircons from each of the four rhyolite eruptions that built Tarawera volcano in the last 22 ka within the Okataina Volcanic Center (OVC), caldera, New Zealand. Secondary ion mass spectrometry analyses on unpolished euhedral crystal faces that lack resorption features show that crystal growth variously terminated from near-eruption age to ~100 ka prior to eruption. Age-depth profiling of crystals reveals long periods of continuous (~34 ka) and discontinuous growth (~90 ka). Growth hiatuses of up to ~40 ka duration occur, but do not all relate to obvious resorption surfaces. Age differences up to similar magnitude are found on opposing faces of some crystals suggesting episodes of partial exposure to melts. These features are best explained by periodic, complete, or partial, sub-solidus storage and/or inclusion in larger crystal phases, followed by rapid liberation prior to eruption. This is supported by high abundances of U and Th (~500 − >2,000 ppm) in some zircons consistent with periods of high crystallinity (>70%) in the magmatic system, based on crystal/melt partitioning. Contemporaneous but contrasting rim-ward trends of these elements within crystals, even in the same lava hand sample, require synchronous growth in separate melt bodies and little connectivity within the system, but also significant crystal transport and mixing prior to eruption. Many crystals record continuity of growth through the preceding ~60 ka OVC caldera-collapse and subsequent eruptions from Tarawera. This demonstrates a decoupling between eruption triggers, such as shallow crustal extension and mafic intrusion, and the crystallization state of the OVC silicic magmatic system. The data highlights the need to distinguish between the time for accumulation of eruptible magma and the long-term magma residence time based on the age of crystals with high closure temperatures, when assessing the potential for catastrophic eruptions.     

The behavior of chlorine and sulfur during differentiation of the Mt. Somma-Vesuvius magmatic system

Summary Reheated silicate melt inclusions in volcanic rock samples from Mt. Somma-Vesuvius, Italy, have been analyzed for 29 constituents including H₂O, S, Cl, F, B, and P₂O₅. This composite volcano consists of the older Mt. Somma caldera, formed between 14 and 3.55 ka before present, and the younger Vesuvius cone. The melt inclusion compositions provide important constraints on pre-eruptive magma geochemistry, identify relationships that relate to eruption behavior and magma evolution, and provide extensive evidence for magmatic fluid exsolution well before eruption. The melt inclusion data have been categorized by groups that reflect magma compositions, age, and style of eruptions. The data show distinct differences in composition for eruptive products older than 14.0 ka (pre-caldera rocks) versus eruptive products younger than 3.55 ka. Moreover, pre-caldera eruptions were associated with magmas relatively enriched in SiO₂, whereas eruptions younger than 3.55 ka (i.e., the syn- and post-caldera magmas which generated the Somma caldera and the Vesuvius cone) were derived from magmas comparatively enriched in S, Cl, CaO, MgO, P₂O₅, F, and many lithophile trace elements. Melt inclusion data indicate that eruptive behavior at Vesuvius correlates with pre-eruptive volatile enrichments. Most magmas associated with explosive plinian and subplinian events younger than 3.55 ka contained more H₂O, contained significantly more S, and exhibited higher (S/Cl) ratios than syn- and post-caldera magmas which erupted during relatively passive interplinian volcanic phenomena. Received January 10, 2000 Revised version accepted July 17, 2000     

Extraction of crystal-poor rhyolite from a hornblende-bearing intermediate mush: a case study of the caldera-forming Matahina eruption,Okataina volcanic complex

The Matahina Ignimbrite (~160 km³ rhyolite magma, 330 ka) was deposited during a caldera-forming eruption from the Okataina Volcanic Centre, Taupo Volcanic Zone (TVZ), New Zealand. Juvenile clasts are divided into three groups: Group (1) the dominant crystal-poor rhyolite type, Group (2) a minor coarse-grained, mingled/mixed intermediate type, and Group (3) a rare fine-grained basalt. The ignimbrite consists of the Group 1 type and is divided into three members: a lower and middle member, which is high-silica, crystal-poor (<10 vol.%) rhyolite, and the upper member, which is low-silica and slightly more crystal-rich (up to 21 vol.%). Cognate, crystal-rich (up to 50 vol.%) basalt to intermediate pumice occurs on top of lag breccias and within lithic-rich pyroclastic density current deposits along the caldera margin (Groups 2 and 3). Several lines of evidence indicate that the intermediate clasts represent the cumulate complement to the melt-rich rhyolite: (1) continuity in the compositions of plagioclase, orthopyroxene, hornblende, and oxides and normal zoning of individual phenocrysts; (2) the silicic glass from the intermediate magma (interstitial melt) overlaps compositionally with the bulk rock rhyolite and glass; (3) high Zr and a slight positive Eu anomaly in the intermediate magma relative to quenched enclaves from other intermediate TVZ eruptions indicates zircon and plagioclase accumulation, respectively; (4) an increase in the Cl contents in glass from the least evolved to most evolved is consistent with the concentration of volatiles during magma evolution. Most of the compositional variations in the low- to high-silica rhyolites can be accounted for by continued Rayleigh fractionation (up to 15%), following melt extraction from the underlying mush, under varying fO₂–fH₂O conditions to form a slightly compositionally zoned rhyolitic cap. This link to the varying fO₂–fH₂O conditions is evidenced by the strong correlation between key geochemical parameters (e.g. Dy, Y), that qualitatively reflect fH₂O conditions (presence or absence of hornblende/biotite), and fO₂ estimated from Fe–Ti oxide equilibrium. Magma mingling/mixing between the basalt–andesite and the main slightly compositionally zoned rhyolitic magma occurred during caldera-collapse, modifying the least-evolved rhyolite at the lower portion of the reservoir and effectively destroying any pre-eruptive gradients.     

The origin of high δ<Superscript>18</Superscript>O zircons: marbles,megacrysts, and metamorphism

The oxygen isotope ratios (δ¹⁸O) of most igneous zircons range from 5 to 8‰, with 99% of published values from 1345 rocks below 10‰. Metamorphic zircons from quartzite, metapelite, metabasite, and eclogite record δ¹⁸O values from 5 to 17‰, with 99% below 15‰. However, zircons with anomalously high δ¹⁸O, up to 23‰, have been reported in detrital suites; source rocks for these unusual zircons have not been identified. We report data for zircons from Sri Lanka and Myanmar that constrain a metamorphic petrogenesis for anomalously high δ¹⁸O in zircon. A suite of 28 large detrital zircon megacrysts from Mogok (Myanmar) analyzed by laser fluorination yields δ¹⁸O from 9.4 to 25.5‰. The U–Pb standard, CZ3, a large detrital zircon megacryst from Sri Lanka, yields δ¹⁸O = 15.4 ± 0.1‰ (2 SE) by ion microprobe. A euhedral unzoned zircon in a thin section of Sri Lanka granulite facies calcite marble yields δ¹⁸O = 19.4‰ by ion microprobe and confirms a metamorphic petrogenesis of zircon in marble. Small oxygen isotope fractionations between zircon and most minerals require a high δ¹⁸O source for the high δ¹⁸O zircons. Predicted equilibrium values of Δ¹⁸O(calcite-zircon) = 2–3‰ from 800 to 600°C show that metamorphic zircon crystallizing in a high δ¹⁸O marble will have high δ¹⁸O. The high δ¹⁸O zircons (>15‰) from both Sri Lanka and Mogok overlap the values of primary marine carbonates, and marbles are known detrital gemstone sources in both localities. The high δ¹⁸O zircons are thus metamorphic; the 15–25‰ zircon values are consistent with a marble origin in a rock-dominated system (i.e., low fluid_(external)/rock); the lower δ¹⁸O zircon values (9–15‰) are consistent with an origin in an external fluid-dominated system, such as skarn derived from marble, although many non-metasomatized marbles also fall in this range of δ¹⁸O. High δ¹⁸O (>15‰) and the absence of zoning can thus be used as a tracer to identify a marble source for high δ¹⁸O detrital zircons; this recognition can aid provenance studies in complex metamorphic terranes where age determinations alone may not allow discrimination of coeval source rocks. Metamorphic zircon megacrysts have not been reported previously and appear to be associated with high-grade marble. Identification of high δ¹⁸O zircons can also aid geochronology studies that seek to date high-grade metamorphic events due to the ability to distinguish metamorphic from detrital zircons in marble.     

Temporal and chemical connections between plutons and ignimbrites from the Mount Princeton magmatic center

The Mount Princeton magmatic center, located in central Colorado, consists of the epizonal Mount Princeton batholith, the nested Mount Aetna caldera, and volumetrically minor leucogranites. New CA-TIMS U/Pb zircon ages indicate the majority of the Mount Princeton batholith was emplaced during a period of regional ignimbrite quiescence. The structurally highest unit of quartz monzonite yields a ²⁰⁶Pb/²³⁸U age of 35.80 ± 0.10 Ma, and the youngest dated unit of the quartz monzonite is a porphyritic unit that yields a ²⁰⁶Pb/²³⁸U age of 35.37 ± 0.10 Ma. Using the exposed, dated volume of the quartz monzonite and new geochronology yields an estimated pluton filling rate of ~0.002 km³/a. This rate is comparable to the accumulation rates published for other plutons, and at least an order of magnitude slower than fluxes necessary to support accumulation of large eruptible magma volumes. Geochronology for the two large ignimbrites spatially associated with the batholith indicates a temporal disconnect between the vast majority of pluton building and explosive eruption of magma. The Wall Mountain Tuff erupted from a source in the same geographic area as the Mount Princeton batholith at 37.3 Ma (Ar/Ar sanidine), but no structural evidence of a caldera or temporally associated plutonic rocks is known. The Badger Creek Tuff erupted at 34.3 Ma (Ar/Ar sanidine) during the formation of the Mount Aetna caldera in the southern portion of the batholith. Our ²⁰⁶Pb/²³⁸U age for the Badger Creek Tuff is 34.47 ± 0.05. The only analyzed plutonic rocks of similar age to the Badger Creek Tuff are an extra-caldera dike with a ²⁰⁶Pb/²³⁸U age of 34.57 ± 0.08 Ma, a ring dike with a ²⁰⁶Pb/²³⁸U age of 34.48 ± 0.09 Ma, and a portion of the Mount Aetna pluton with a ²⁰⁶Pb/²³⁸U age of 34.60 ± 0.13 Ma. The small volume intrusions related to the eruption of the Badger Creek Tuff are chemically similar to the ignimbrite and show no signature of crystal–liquid separation in the shallow crust.     

Development of the Long Valley,California, magma chamber recorded in precaldera rhyolite lavas of Glass Mountain

Glass Mountain, California, consists of >50 km³ of high-silica rhyolite lavas and associated pyroclastic deposits that erupted over a period of >1 my preceding explosive eruption of the Bishop Tuff and formation of the Long Valley caldera at 0.73 Ma. These “minimum-melt” rhyolites yield Fe-Ti-oxide temperatures of 695–718°C and contain sparse phenocrysts of plagioclase+quartz+magnetite+apatite±sanidine, biotite, ilmenite, allanite, and zircon. Incompatible trace elements show similar or larger ranges within the Glass Mountain suite than within the Bishop Tuff, despite a much smaller range of major-element concentrations, largely due to variability among the older lavas (erupted between 2.1 and 1.2 Ma). Ratios of the most incompatible elements have larger ranges in the older lavas than in the younger lavas (1.2–0.79 Ma), and concentrations of incompatible elements span wide ranges at nearly constant Ce/Yb, suggesting that the highest concentrations of these elements are not the result of extensive fractional crystallization alone; rather, they are inherited from parental magmas with a larger proportion of crustal partial melt. Evidence for the nature of this crustal component comes from the presence of scarce, tiny xenocrysts derived from granitic and greenschist-grade metamorphic rocks. The wider range of chemical and isotopic compositions in the older lavas, the larger range in phenocryst modes, the eruption of magmas with different compositions at nearly the same time in different parts of the field, and the smaller volume of individual lavas suggest either that more than one magma body was tapped during eruption of the older lavas or that a single chamber tapped by all lavas was small enough that the composition of its upper reaches easily affected by new additions of crustal melts. We interpret the relative chemical, mineralogical, and isotopic homogeneity of the younger Glass Mountain lavas as reflecting eruptions from a large, integrated magma chamber. The small number of cruptions between 1.4 and 1.2 ma may have allowed time for a large magma body to coalesce, and, as the chamber grew, its upper reaches became less affected by new inputs of crustal melts, so that trace-element trends in magmas erupted after 1.2 Ma are largely controlled by fractional crystallization. The extremely low Sr concentrations of Glass Mountain lavas imply extensive crystallization in chambers at least hundreds of cubic kilometers in volume. The close similarity in Sr, Nd, and Pb isotopic ratios between the younger Glass Mountain lavas and unaltered Bishop Tuff indicates that they tapped the same body of magma, which had become isotopically homogenous by 1.2 Ma but continued to differentiate after that time. From 1.2 to 0.79 Ma, volumetric eruptive rates may have exceeded rates of differentiation, as younger Glass Mountain lavas become slightly less evolved with time. Early-erupted Bishop Tuff is more evolved than the youngest of the Glass Mountain lavas and is characterized by slightly different trace element ratios. This suggests that although magma had been present for 0.5 my, the composiional gradient exhibited by the Bishop Tuff had not been a long-term, steady-state condition in the Long Valley magma chamber, but developed at least in part during the 0.06-my hiatus between extrusion of the last Glass Mountain lava and the climactic eruption.