Rapid pre-eruptive thermal rejuvenation in a large silicic magma body: the case of the Masonic Park Tuff,Southern Rocky Mountain volcanic field,CO, USA

Determining the mechanisms involved in generating large-volume eruptions (>100 km³) of silicic magma with crystallinities approaching rheological lock-up (~50 vol% crystals) remains a challenge for volcanologists. The Cenozoic Southern Rocky Mountain volcanic field, in Colorado and northernmost New Mexico, USA, produced ten such crystal-rich ignimbrites within 3 m.y. This work focuses on the 28.7 Ma Masonic Park Tuff, a dacitic (~62–65 wt% SiO₂) ignimbrite with an estimated erupted volume of ~500 km³ and an average of ~45 vol% crystals. Near-absence of quartz, titanite, and sanidine, pronounced An-rich spikes near the rims of plagioclase, and reverse zoning in clinopyroxene record the reheating (from ~750 to >800?°C) of an upper crustal mush in response to hotter recharge from below. Zircon U–Pb ages suggest prolonged magmatic residence, while Yb/Dy vs temperature trends indicate co-crystallization with titanite which was later resorbed. High Sr, Ba, and Ti concentrations in plagioclase microlites and phenocryst rims require in-situ feldspar melting and concurrent, but limited, mass addition provided by the recharge, likely in the form of a melt-gas mixture. The larger Fish Canyon Tuff, which erupted from the same location ~0.7 m.y. later, also underwent pre-eruptive reheating and partial melting of quartz, titanite, and feldspars in a long-lived upper crustal mush following the underplating of hotter magma. The Fish Canyon Tuff, however, records cooler pre-eruptive temperatures (~710–760?°C) and a mineral assemblage indicative of higher magmatic water contents (abundant resorbed sanidine and quartz, euhedral amphibole and titanite, and absence of pyroxene). These similar pre-eruptive mush-reactivation histories, despite differing mineral assemblages and pre-eruptive temperatures, indicate that thermal rejuvenation is a key step in the eruption of crystal-rich silicic volcanics over a wide range of conditions.     

Magmatic processes in the Bishop Tuff rhyolitic magma based on trace elements in melt inclusions and pumice matrix glass

To investigate the origin of compositional zonation in the Bishop Tuff magma body, we have analyzed trace elements in the matrix glass of pumice clasts and in quartz-hosted melt inclusions. Our results show contrasting patterns for quartz in different parts of the Bishop Tuff. In all samples from the early part of the eruption, trace element compositions of matrix glasses are similar to but slightly more evolved than quartz-hosted melt inclusions. This indicates a cogenetic relationship between quartz crystals and their surrounding matrix glass, consistent with in situ crystallization. The range of incompatible element concentrations in melt inclusions and matrix glass from single pumice clasts requires 16–20 wt% in situ crystallization. This is greater than the actual crystal content of the pumices (<15 % crystals). In contrast to the pattern for the early pumices, pyroclastic flow samples from the middle part of the eruption show contrasting trends: In some clasts, the matrix is more evolved than the inclusions, whereas in other clasts, the matrix is less evolved. In the late Bishop Tuff, all crystal-rich samples have matrix glasses that are less evolved than the melt inclusions. Trace element abundances indicate that the cores of quartz in the late Bishop Tuff crystallized from more differentiated rhyolitic magma that was similar in many ways, yet distinct from the early-erupted Bishop Tuff. Our results are compatible with a model of secular incremental zoning (Hildreth and Wilson in Compositional zoning of the Bishop Tuff. J Petrol 48(5):951–999, 2007), in which melt batches from underlying crystal mush rise to various levels in a growing magma body according to their buoyancy. Early- and middle-erupted quartz crystallized from highly evolved rhyolitic melt, but then some parts of the middle-erupted magma were invaded by less differentiated rhyolite such that the matrix melt at the time of eruption was less evolved than the melt inclusions. A similar process occurred but to a greater extent in magma that erupted to form the late Bishop Tuff. In addition, there was a final, major magma mixing event in the late magma that formed Ti-rich rims on quartz and Ba-rich rims on sanidine, trapped less evolved rhyolitic melt inclusions, and resulted in dark and swirly crystal-poor pumice that is a rare type throughout much of the Bishop Tuff.     

Experimental constraints on rhyolite-MELTS and the Late Bishop Tuff magma body

Thermodynamic models are vital tools to evaluate magma crystallization and storage conditions. Before their results can be used independently, however, they must be verified with controlled experimental data. Here, we use a set of hydrothermal experiments on the Late-erupted Bishop Tuff (LBT) magma to evaluate the rhyolite-MELTS thermodynamic model, a modified calibration of the original MELTS model optimized for crystallization of silicic magmas. Experimental results that are well captured by rhyolite-MELTS include a relatively narrow temperature range separating the crystallization of the first felsic mineral and the onset of the ternary minimum (quartz plus two feldspars), and extensive crystallization over a narrow temperature range once the ternary minimum is reached. The model overestimates temperatures by ~40 °C, a known limitation of rhyolite-MELTS. At pressures below 110 MPa, model and experiments differ in the first felsic phase, suggesting that caution should be exercised when applying the model to very low pressures. Our results indicate that for quartz, sanidine, plagioclase, magnetite, and ilmenite to crystallize in equilibrium from LBT magma, magma must have been stored at ≤740 °C, even when a substantial amount of CO₂ occurs in the coexisting fluid. Such temperatures are in conflict with the hotter temperatures retrieved from magnetite–ilmenite compositions (~785 °C for the sample used in the experiments). Consistent with other recent studies, we suggest that the Fe–Ti oxide phases in the Late Bishop Tuff magma body are not in equilibrium with the other minerals and thus the retrieved temperature and oxygen fugacity do not reflect pre-eruptive storage conditions.     

Petrological imaging of an active pluton beneath Cerro Uturuncu,Bolivia

Uturuncu is a dormant volcano in the Altiplano of SW Bolivia. A present day ~70 km diameter interferometric synthetic aperture radar (InSAR) anomaly roughly centred on Uturuncu’s edifice is believed to be a result of magma intrusion into an active crustal pluton. Past activity at the volcano, spanning 0.89 to 0.27 Ma, is exclusively effusive and almost all lavas and domes are dacitic with phenocrysts of plagioclase, orthopyroxene, biotite, ilmenite and Ti-magnetite plus or minus quartz, and microlites of plagioclase and orthopyroxene set in rhyolitic groundmass glass. Plagioclase-hosted melt inclusions (MI) are rhyolitic with major element compositions that are similar to groundmass glasses. H₂O concentrations plotted versus incompatible elements for individual samples describe a trend typical of near-isobaric, volatile-saturated crystallisation. At 870 °C, the average magma temperature calculated from Fe–Ti oxides, the average H₂O of 3.2 ± 0.7 wt% and CO₂ typically <160 ppm equate to MI trapping pressures of 50–120 MPa, approximately 2–4.5 km below surface. Such shallow storage precludes the role of dacite magma emplacement into pre-eruptive storage regions as being the cause of the observed InSAR anomaly. Storage pressures, whole-rock (WR) chemistry and phase assemblage are remarkably consistent across the eruptive history of the volcano, although magmatic temperatures calculated from Fe–Ti oxide geothermometry, zircon saturation thermometry using MI and orthopyroxene-melt thermometry range from 760 to 925 °C at NNO ± 1 log. This large temperature range is similar to that of saturation temperatures of observed phases in experimental data on Uturuncu dacites. The variation in calculated temperatures is attributed to piecemeal construction of the active pluton by successive inputs of new magma into a growing volume of plutonic mush. Fluctuating temperatures within the mush can account for sieve-textured cores and complex zoning in plagioclase phenocrysts, resorption of quartz and biotite phenocrysts and apatite microlites. That Fe–Ti oxide temperatures vary by ~50–100 °C in a single thin section indicates that magmas were not homogenised effectively prior to eruption. Phenocryst contents do not correlate with calculated magmatic temperatures, consistent with crystal entrainment from the mush during magma ascent and eruption. Microlites grew during ascent from the magma storage region. Variability in the proportion of microlites is attributed to differing ascent and effusion rates with faster rates in general for lavas >0.5 Ma compared to those <0.5 Ma. High microlite contents of domes indicate that effusion rates were probably slowest in dome-forming eruptions. Linear trends in WR major and trace element chemistries, highly variable, bimodal mineral compositions, and the presence of mafic enclaves in lavas demonstrate that intrusion of more mafic magmas into the evolving, shallow plutonic mush also occurred further amplifying local temperature fluctuations. Crystallisation and resorption of accessory phases, particularly ilmenite and apatite, can be detected in MI and groundmass glass trace element covariation trends, which are oblique to WRs. Marked variability of Ba, Sr and La in MI can be attributed to temperature-controlled, localised crystallisation of plagioclase, orthopyroxene and biotite within the evolving mush.     

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.     

Evolution of Bishop Tuff Rhyolitic Magma Based on Melt and Magnetite Inclusions and Zoned Phenocrysts

The evolution of large bodies of silicic magma is an importantaspect of planetary differentiation. Melt and mineral inclusionsin phenocrysts and zoned phenocrysts can help reveal the processesof differentiation such as magma mixing and crystal settling,because they record a history of changing environmental conditions.Similar major element compositions and unusually low concentrationsof compatible elements (e.g. 0·45–4·6 ppmBa) in early-erupted melt inclusions, matrix glasses and bulkpumice from the Bishop Tuff, California, USA, suggest eutectoidfractional crystallization. On the other hand, late-eruptedsanidine phenocrysts have rims rich in Ba, and late-eruptedquartz phenocrysts have CO₂-rich melt inclusions closest tocrystal rims. Both features are the reverse of in situ crystallizationdifferentiation, and they might be explained by magma mixingor crystal sinking. Log(Ba/Rb) correlates linearly with log(Sr/Rb)in melt inclusions, and this is inconsistent with magma mixing.Melt inclusion gas-saturation pressure increases with CO₂ fromphenocryst core to rim and suggests crystal sinking. Some inclusionsof magnetite in late-erupted quartz are similar to early-eruptedmagnetite phenocrysts, and this too is consistent with crystalsinking. We argue that some large phenocrysts of late-eruptedquartz and sanidine continued to crystallize as they sank severalkilometers through progressively less differentiated melts.Probable diffusive modification of Sr in sanidine phenocrystsand the duration of crystal sinking are consistent with an evolutionaryinterval of some 100 ky or more. Crystal sinking enhanced thedegree of differentiation of the early-erupted magma and pointsto the importance of H₂O (to diminish viscosity and enhancethe rate of crystal sinking) in the evolution of silicic magmas. KEY WORDS: crystal settling; differentiation; melt inclusions; rhyolite; trace elements     

The Blacktail Creek Tuff: an analytical and experimental study of rhyolites from the Heise volcanic field,Yellowstone hotspot system

The magma storage conditions of the 6.62 Ma Blacktail Creek Tuff eruption, belonging to the Heise volcanic field (6.62–4.45 Ma old) of the Yellowstone hotspot system, have been investigated by combining thermobarometric and experimental approaches. The results from different geothermometers (e.g., Fe–Ti oxides, feldspar pairs, apatite and zircon solubility, and Ti in quartz) indicate a pre-eruptive temperature in the range 825–875 °C. The temperature estimated using two-pyroxene pairs varies in a range of 810–950 °C, but the pyroxenes are probably not in equilibrium with each other, and the analytical results of melt inclusion in pyroxenes indicate a complex history for clinopyroxene, which hosts two compositionally different inclusion types. One natural Blacktail Creek Tuff rock sample has been used to determine experimentally the equilibrium phase assemblages in the pressure range 100–500 MPa and a water activity range 0.1–1.0. The experiments have been performed at fluid-present conditions, with a fluid phase composed of H₂O and CO₂, as well as at fluid-absent conditions. The stability of the quartzo-feldspathic phases is similar in both types of experiments, but the presence of mafic minerals such as biotite and clinopyroxene is strongly dependent on the experimental approach. Possible explanations are given for this discrepancy which may have strong impacts on the choice of appropriate experimental approaches for the determination of magma storage conditions. The comparison of the composition of natural phases and of experimentally synthesized phases confirms magma storage temperatures of 845–875 °C. Melt water contents of 1.5–2.5 wt% H₂O are required to reproduce the natural Blacktail Creek Tuff mineral assemblage at these temperatures. Using the Ti-in-quartz barometer and the Qz–Ab–Or proportions of natural matrix glasses, coexisting with quartz, plagioclase and sanidine, the depth of magma storage is estimated to be in a pressure range between 130 and 250 MPa.     

Building zoned ignimbrites by recycling silicic cumulates: insight from the 1,000 km<Superscript>3</Superscript> Carpenter Ridge Tuff,CO

The ~1,000 km³ Carpenter Ridge Tuff (CRT), erupted at 27.55 Ma during the mid-tertiary ignimbrite flare-up in the western USA, is among the largest known strongly zoned ash-flow tuffs. It consists primarily of densely welded crystal-poor rhyolite with a pronounced, highly evolved chemical signature (high Rb/Sr, low Ba, Zr, Eu), but thickly ponded intracaldera CRT is capped by a more crystal-rich, less silicic facies. In the outflow ignimbrite, this upper zone is defined mainly by densely welded crystal-rich juvenile clasts of trachydacite composition, with higher Fe–Ti oxide temperatures, and is characterized by extremely high Ba (to 7,500 ppm), Zr, Sr, and positive Eu anomalies. Rare mafic clasts (51–53 wt% SiO₂) with Ba contents to 4,000–5,000 ppm and positive Eu anomalies are also present. Much of the major and trace-element variations in the CRT juvenile clasts can be reproduced via in situ differentiation by interstitial melt extraction from a crystal-rich, upper-crustal mush zone, with the trachydacite, crystal-rich clasts representing the remobilized crystal cumulate left behind by the melt extraction process. Late recharge events, represented by the rare mafic clasts and high-Al amphiboles in some samples, mixed in with parts of the crystal cumulate and generated additional scatter in the whole-rock data. Recharge was important in thermally remobilizing the silicic crystal cumulate by partially melting the near-solidus phases, as supported by: (1) ubiquitous wormy/sieve textures and reverse zoning patterns in feldspars and biotites, (2) absence of quartz in this very silicic unit stored at depths of >4–5 km, and (3) heterogeneous melt compositions in the trachydacite fiamme and mafic clasts, particularly in Ba, indicating local enrichment of this element due mostly to sanidine and biotite melting. The injection of hot, juvenile magma into the upper-crustal cumulate also imparted the observed thermal gradient to the deposits and the mixing overprint that partly masks the in situ differentiation process. The CRT provides a particularly clear perspective on processes of in situ crystal-liquid separation into a lower crystal-rich zone and an upper eruptible cap, which appears common in incrementally built upper-crustal magma reservoirs of high-flux magmatic provinces.     

The Bishop Tuff giant magma body: an alternative to the Standard Model

The Bishop Tuff, one of the most extensively studied high-silica rhyolite bodies in the world, is usually considered as the archetypical example of a deposit formed from a magma body characterized by thermal and compositional vertical stratification—what we call the Standard Model for the Bishop magma body. We present here new geothermometry and geobarometry results derived using a large database of previously published quartz-hosted glass inclusion compositions. Assuming equilibrium between melt and an assemblage composed of quartz, ±plagioclase, ±sanidine, +zircon, ±fluid, we use Zr contents in glass inclusions to derive quartz crystallization temperatures, and we use (1) silica contents in glass, (2) projection of glass compositions onto the haplogranitic (quartz-albite-orthoclase) ternary, and (3) phase equilibria calculations using rhyolite-MELTS, to constrain crystallization pressures. We find crystallization temperatures of ~740–750 °C for all inclusions from both early- and late-erupted pumice. Crystallization pressures for both early- and late-erupted inclusions are also very similar to each other, with averages of ~175–200 MPa. We find no evidence of late-erupted inclusions having been entrapped at higher temperatures or pressures than early-erupted inclusions, as would be expected by the Standard Model. We argue that the thermal gradient inferred from Fe–Ti oxides—the backbone of the Standard Model—does not reflect equilibrium pre-eruptive conditions; we also note that H₂O–CO₂ systematics of glass inclusions yields overlapping pressure ranges for early- and late-erupted inclusions, similar to the results presented here; and we show that glass inclusion and phenocryst compositions show bimodal distributions, suggestive of compositional separation between early- and late-erupted populations. These findings are inconsistent with the Standard Model. The similarity in crystallization conditions and the compositional separation between early- and late-erupted magmas suggest that two laterally juxtaposed independent magma reservoirs existed in the same region at the same time and co-erupted to form the Long Valley Caldera and the Bishop Tuff. This hypothesis would explain the lack of mixing between early- and late-erupted crystal populations in pumice clasts; it could also explain the inferred eruption pattern—which resulted in early-erupted magmas being deposited only to the south of the caldera—if the early-erupted magma body resided to the south and the late-erupted magma body was located to the north. Our alternative model is consistent with the patchy distribution of thermal anomalies and the inference of co-eruption of distinct magma types in active volcanic areas such as the central Taupo Volcanic Zone.     

A New LA‐ICP‐MS Method for Ti in Quartz: Implications and Application to High Pressure Rutile‐Quartz Veins from the Czech Erzgebirge

Experimental determination of the pressure and temperature controls on Ti solubility in quartz provides a calibration of the Ti‐in‐quartz (TitaniQ) geothermometer applicable to geological conditions up to ~ 20 kbar. We present a new method for determining ⁴⁸Ti mass fractions in quartz by LA‐ICP‐MS at the 1 μg g^?1 level, relevant to quartz in HP‐LT terranes. We suggest that natural quartz such as the low‐CL rims of the Bishop Tuff quartz (determined by EPMA; 41 ± 2 μg g^?1 Ti, 2s) is more suitable than NIST reference glasses as a reference material for low Ti mass fractions because matrix effects are limited, Ca isobaric interferences are avoided, and polyatomic interferences at mass 48 are insignificant, thus allowing for the use of ⁴⁸Ti as a normalising mass. Average titanium mass fraction from thirty‐three analyses of low temperature quartz from the Czech Erzgebirge is 0.9 ± 0.2 μg g^?1 (2s) using ⁴⁸Ti as a normalising mass and Bishop Tuff quartz rims as a reference material. The 2s average analytical uncertainty for individual analyses of ⁴⁸Ti is 8% for 50 μm spots and 7% for 100 μm spots, which offers much greater accuracy than the 21–41% uncertainty (2s) incurred from using ⁴⁹Ti as an analyte.     

Compositional gradients in large reservoirs of silicic magma as evidenced by ignimbrites versus Taylor Creek Rhyolite lava domes

The Taylor Creek Rhyolite of southwest New Mexico consists of 20 lava domes and flows that were emplaced during a period of a few thousand years or less in late Oligocene time. Including genetically associated pyroclastic deposits, which are about as voluminous as the lava domes and flows, the Taylor Creek Rhyolite represents roughly 100 km³ of magma erupted from vents distributed throughout an area of several hundred square kilometers. Major-element composition is metaluminous to weakly peraluminous high-silica rhyolite and is nearly constant throughout the lava field. The magma reservoir for the Taylor Creek Rhyolite was vertically zoned in trace elements, ⁸⁷Sr/⁸⁶Sr, and phenocryst abundance and size. Mean trace-element concentrations, ranges in concentrations, and element-pair correlations are similar to many subalkaline silicic ignimbrites. However, the polarity of the zonation was opposite that in reservoirs for ignimbrites, for most constituents. For example, compared to the Bishop Tuff, only ⁸⁷Sr/⁸⁶Sr and Sc increased upward in both reservoirs. Quite likely, a dominant but nonerupted volume of the magma reservoir for the Taylor Creek Rhyolite was zoned like that for the Bishop Tuff, whereas an erupted, few-hundred-meter-thick cap on the magma body was variably contaminated by roof rocks whose contribution to this part of the magma system moderated relatively extreme trace-element concentrations of uncontaminated Taylor Creek Rhyolite but did not change the sense of correlation for most element pairs. The contaminant probably was a Precambrian rock of broadly granitic composition and with very high ⁸⁷Sr/⁸⁶Sr. Although examples apparently are not yet reported in the literature, evidence for a similar thin contaminated cap on reservoirs for large-volume silicic ignimbrites may exist in the bottom few meters of ignimbrites or perhaps only in the pumice fallout that normally immediately precedes ignimbrite emplacement. ⁸⁷Sr/⁸⁶Sr in sanidine phenocrysts of the Taylor Creek Rhyolite is higher than that of their host whole rocks. Covariation of this isotope ratio with sanidine abundance and size indicates positive correlations for all three features with decreasing distance to the roof of the magma reservoir. The sanidine probably is more radiogenic than host whole rock because growing phenocrysts partly incorporated Sr from the first partial melt of roof rocks, which contained the highly radiogenic Sr of Precambrian biotite ± hornblende, whereas diffusion was too slow for sanidine to incorporate much of the Sr from subsequently produced less radiogenic partial melt of roof rocks, before eruption quenched the magma system. Disequilibrium between feldspar phenocrysts and host groundmass is fairly common for ignimbrites, and a process of contamination similar to that for the Taylor Creek Rhyolite may help explain some of these situations.     

Multiple timescale constraints for high-flux magma chamber assembly prior to the Late Bronze Age eruption of Santorini (Greece)

The rhyodacitic magma discharged during the 30–80 km³ DRE (dense rock equivalent) Late Bronze Age (LBA; also called ‘Minoan’) eruption of Santorini caldera is known from previous studies to have had a complex history of polybaric ascent and storage prior to eruption. We refine the timescales of these processes by modelling Mg–Fe diffusion profiles in orthopyroxene and clinopyroxene crystals. The data are integrated with previously published information on the LBA eruption (phase equilibria studies, melt inclusion volatile barometry, Mg-in-plagioclase diffusion chronometry), as well as new plagioclase crystal size distributions and the established pre-LBA history of the volcano, to reconstruct the events that led up to the assembly and discharge of the LBA magma chamber. Orthopyroxene, clinopyroxene and plagioclase crystals in the rhyodacite have compositionally distinct rims, overgrowing relict, probably source-derived, more magnesian (or calcic) cores, and record one or more crystallization (plag???opx?>?cpx) events during the few centuries to years prior to eruption. The crystallization event(s) can be explained by the rapid transfer of rhyodacitic melt from a dioritic/gabbroic region of the subcaldera pluton (mostly in the 8–12 km depth range), followed by injection, cooling and mixing in a large melt lens at 4–6 km depth (the pre-eruptive magma chamber). Since crystals from all eruptive phases yield similar timescales, the melt transfer event(s), the last of which took place less than 2 years before the eruption, must have involved most of the magma that subsequently erupted. The data are consistent with a model in which prolonged generation, storage and segregation of silicic melts were followed by gravitational instability in the subcaldera pluton, causing the rapid interconnection and amalgamation of melt-rich domains. The melts then drained to the top of the pluton, at fluxes of up to 0.1–1 km³ year^??1, where steep vertical gradients of density and rheology probably caused them to inject laterally, forming a short-lived holding chamber prior to eruption. This interpretation is consistent with growing evidence that some large silicic magma chambers are transient features on geological timescales. A similar process preceded at least one earlier caldera-forming eruption on Santorini, suggesting that it may be a general feature of this rift-hosted magmatic system.     

A cascade of magmatic events during the assembly and eruption of a super-sized magma body

We use comprehensive geochemical and petrological records from whole-rock samples, crystals, matrix glasses and melt inclusions to derive an integrated picture of the generation, accumulation and evacuation of 530 km³ of crystal-poor rhyolite in the 25.4 ka Oruanui supereruption (New Zealand). New data from plagioclase, orthopyroxene, amphibole, quartz, Fe–Ti oxides, matrix glasses, and plagioclase- and quartz-hosted melt inclusions, in samples spanning different phases of the eruption, are integrated with existing data to build a history of the magma system prior to and during eruption. A thermally and compositionally zoned, parental crystal-rich (mush) body was developed during two periods of intensive crystallisation, 70 and 10–15 kyr before the eruption. The mush top was quartz-bearing and as shallow as ~3.5 km deep, and the roots quartz-free and extending to >10 km depth. Less than 600 year prior to the eruption, extraction of large volumes of ~840 °C low-silica rhyolite melt with some crystal cargo (between 1 and 10%), began from this mush to form a melt-dominant (eruptible) body that eventually extended from 3.5 to 6 km depth. Crystals from all levels of the mush were entrained into the eruptible magma, as seen in mineral zonation and amphibole model pressures. Rapid translation of crystals from the mush to the eruptible magma is reflected in textural and compositional diversity in crystal cores and melt inclusion compositions, versus uniformity in the outermost rims. Prior to eruption the assembled eruptible magma body was not thermally or compositionally zoned and at temperatures of ~790 °C, reflecting rapid cooling from the ~840 °C low-silica rhyolite feedstock magma. A subordinate but significant volume (3–5 km³) 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.     

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.     

Ultrahigh-temperature metamorphism from an unusual corundum+orthopyroxene intergrowth bearing Al–Mg granulite from the Southern Marginal Zone,Limpopo Complex,South Africa

The Southern Marginal Zone of the Limpopo Complex is composed of granite-greenstone cratonic rocks reworked by a Neoarchean high-grade tectono-metamorphic event. Petrographic and mineral chemical characterization of an Al–Mg granulite from this zone is presented here. The granulite has a gneissic fabric with distinct Al-rich and Si-rich layers, with the former preserving the unusual lamellar (random and regular subparallel) intergrowths of corundum and symplectic intergrowth of spinel with orthopyroxene. The Al-rich layer preserves mineral assemblages such as rutile with orthopyroxene + sillimanite ± quartz, Al-rich orthopyroxene (~11 wt%), spinel + quartz, and corundum in possible equilibrium with quartz, while the Si-rich layer preserves antiperthites and orthopyroxene + sillimanite ± quartz, all considered diagnostic of ultrahigh-temperature metamorphism. Application of Al-in-opx thermometry, ternary feldspar thermometry and construction of suitable pressure–temperature phase diagrams, compositional and model proportion isopleth results indicate P–T conditions as high as ~1,050–1,100 °C, and ~10–12 kbars for the Al–Mg granulite. Our report of ultrahigh-temperature conditions is significant considering that the very high temperature was reached during decompression of an otherwise high-pressure granulite complex (clockwise P–T path), whereas most other ultrahigh-temperature granulites are linked to magma underplating at the base of the crust (counterclockwise P–T path).     

Microlite transfer by disaggregation of mafic inclusions following magma mixing at Soufrière Hills volcano,Montserrat

The Soufrière Hills volcano on Montserrat has for the past 12 years been erupting andesite with basaltic to basaltic–andesite inclusions. The andesite contains a wide variety of phenocryst textures and strongly zoned microlites. Analysis of minor elements in both phenocrysts and microlites allows us to put detailed constraints on their origins. Compositions of clinopyroxene, from overgrowth rims on quartz and orthopyroxene and coarse-grained breakdown rims on hornblende, are identical to those from the mafic inclusions, indicating that these rims form during interaction with mafic magma. In contrast, resorbed quartz and reversely zoned orthopyroxenes form during heating. Microlites of plagioclase and orthopyroxene are chemically distinct from the phenocrysts, being enriched in Fe and Mg, and Al and Ca respectively. However, microlites of plagioclase, orthopyroxene and clinopyroxene are indistinguishable from the compositions of these phases in the mafic inclusions. We infer that the inclusions disaggregated under conditions of high shear stress during ascent in the conduit, transferring mafic material into the andesite groundmass. The mafic component of the system is therefore greater than previously thought. The presence of mafic-derived microlites in the andesite groundmass also means that care must be taken when using this as a starting material for phase equilibrium experiments.     

Magmatic residence times of zoned phenocrysts: introduction and application of the binary element diffusion modelling (BEDM) technique

This paper describes a general technique, binary element diffusion modelling (BEDM), for determining single-crystal residence times in magmas that relies on modelling the diffusion of two or more elements in the crystal. BEDM has the advantage over other diffusion-based models in that it does not need a precisely defined initial compositional profile for the crystal at “zero time”, and instead requires that the concentrations of the two elements are correlated during crystallisation. Any differences subsequently observed between the two elements are caused by intracrystalline diffusion during residence in hot magma. These differences are removed by artificially ageing the slower-diffusing of the two elements, and the amount of time taken to “undo” the difference between the elements is simply related to the crystal residence time (=decoupling time) at high temperatures. The BEDM principle is demonstrated using artificial data and is then applied to literature data for Sr and Ba in a zoned sanidine crystal from the Bishop Tuff (Anderson et al., in J. Petrol 41(3):449–473, 2000). For this crystal, the method gives a residence time estimate of 114 ka at 800°C, which is then compared with estimates from other methods. In theory, the method can be further expanded for use as a geothermometer as well as geochronometer. However, this is not easily possible with the diffusivity data currently available.     

Compositional Zoning of the Bishop Tuff

Compositional data for >400 pumice clasts, organized accordingto eruptive sequence, crystal content, and texture, providenew perspectives on eruption and pre-eruptive evolution of the>600 km³ of zoned rhyolitic magma ejected as the Bishop Tuffduring formation of Long Valley caldera. Proportions and compositionsof different pumice types are given for each ignimbrite packageand for the intercalated plinian pumice-fall layers that eruptedsynchronously. Although withdrawal of the zoned magma was lesssystematic than previously realized, the overall sequence displaystrends toward greater proportions of less evolved pumice, morecrystals (0·5–24 wt %), and higher FeTi-oxide temperatures(714–818°C). No significant hiatus took place duringthe 6 day eruption of the Bishop Tuff, nearly all of which issuedfrom an integrated, zoned, unitary reservoir. Shortly beforeeruption, however, the zoned melt-dominant portion of the chamberwas invaded by batches of disparate lower-silica rhyolite magma,poorer in crystals than most of the resident magma but slightlyhotter and richer in Ba, Sr, and Ti. Interaction with residentmagma at the deepest levels tapped promoted growth of Ti-richrims on quartz, Ba-rich rims on sanidine, and entrapment ofnear-rim melt inclusions relatively enriched in Ba and CO₂.Varied amounts of mingling, even in higher parts of the chamber,led to the dark gray and swirly crystal-poor pumices sparselypresent in all ash-flow packages. As shown by FeTi-oxide geothermometry,the zoned rhyolitic chamber was hottest where crystal-richest,rendering any model of solidification fronts at the walls orroof unlikely. The main compositional gradient (75–195ppm Rb; 0·8–2·2 ppm Ta; 71–154 ppmZr; 0·40–1·73% FeO*) existed in the melt,prior to crystallization of the phenocryst suite observed, whichincluded zircon as much as 100 kyr older than the eruption.The compositions of crystals, though themselves largely unzoned,generally reflect magma temperature and the bulk compositionalgradient, implying both that few crystals settled or were transportedfar and that the observed crystals contributed little to establishingthat gradient. Upward increases in aqueous gas and dissolvedwater, combined with the adiabatic gradient (for the 5 km depthrange tapped) and the roofward decline in liquidus temperatureof the zoned melt, prevented significant crystallization againstthe roof, consistent with dominance of crystal-poor magma earlyin the eruption and lack of any roof-rind fragments among theBishop ejecta, before or after onset of caldera collapse. Amodel of secular incremental zoning is advanced wherein numerousbatches of crystal-poor melt were released from a mush zone(many kilometers thick) that floored the accumulating rhyoliticmelt-rich body. Each batch rose to its own appropriate levelin the melt-buoyancy gradient, which was self-sustaining againstwholesale convective re-homogenization, while the thick mushzone below buffered it against disruption by the deeper (non-rhyolitic)recharge that augmented the mush zone and thermally sustainedthe whole magma chamber. Crystal–melt fractionation wasthe dominant zoning process, but it took place not principallyin the shallow melt-rich body but mostly in the pluton-scalemush zone before and during batchwise melt extraction. KEY WORDS: Bishop Tuff; ignimbrite; magma zonation; mush model; rhyolite     

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.