skip to main content


Title: Origin of Compositional Gradients with Temperature in the High-SiO2 Rhyolite Portion of the Bishop Tuff: Constraints on Mineral–Melt–Fluid Reactions in the Parental Mush
Abstract The Bishop Tuff (BT), erupted from the Long Valley caldera in California, displays two types of geochemical gradients with temperature: one is related to magma mixing, whereas the other is found in the high-SiO2 rhyolite portion of the Bishop Tuff and is characterized by twofold or lower concentration variations in minor and trace elements that are strongly correlated with temperature. It is proposed that the latter zonation, which preceded phenocryst growth, developed as a result of mineral–melt partitioning between interstitial melt and surrounding crystals in a parental mush, from which variable melt fractions were segregated. To test this hypothesis, trends of increasing vs decreasing element concentrations with temperature (as a proxy for melt fraction), obtained from published data on single-clast pumice samples from the high-SiO2 rhyolite portion of the Bishop Tuff, were used to infer their relative degrees of incompatibility vs compatibility between crystals and melt in the parental mush. Relative compatibility values (RCVi) for all elements i, defined as the concentration slope with temperature divided by average concentration, are shown to be linearly correlated with their respective bulk partition coefficients (bulk Di). Mineral–melt partition coefficients from the literature were used to constrain the average stoichiometry of the crystallization/melting reaction in the parental mush: 32 % quartz + 34 % plagioclase + 31 % K-feldspar + 1·60 % biotite + 0·42 % titanomagnetite + 0·34 % ilmenite + 0·093 % allanite + 0·024 % zircon + 0·025 % apatite = 100 % liquid. The proportions of tectosilicates in the reaction (i.e. location of eutectic) are consistent with depths of melt segregation of ~400–550 MPa and an activity of H2O of ~0·4–0·6. Temperatures of <770–780 °C are constrained by allanite in the reaction. Evidence that a fluid phase was present in the parental mush is seen in the decreasing versus increasing H2O and CO2 contents with temperature in the segregated interstitial melt that formed the high-SiO2 rhyolite portion of the Bishop Tuff. The presence of an excess fluid phase, which strongly partitions CO2 relative to the melt, is required to explain the compatible behavior of CO2, whereas the fluid abundance must have been low to explain the incompatible behavior of H2O. Calculated degassing paths for interstitial melts, which segregated from the parental mush and ascended to shallower depths to grow phenocrysts, match published volatile analyses in quartz-hosted melt inclusions and constrain fluid abundances in the mush to be ≤1 wt%. The source of volatiles in the parental mush, irrespective of whether it formed by crystallization or partial melting, must have been primarily from associated basalts, as granitoid crust is too volatile poor. Approximately twice as much basalt as rhyolite is needed to provide the requisite volatiles. The determination of bulk Di for several elements gives the bulk composition of the parental leucogranitic mush and shows that it is distinct from Mesozoic Sierran arc granitoids, as expected. Collectively, the results from this study provide new constraints for models of the complex, multi-stage processes throughout the Plio-Quaternary, involving both mantle-derived basalt and pre-existing crust, that led to the origin of the parental body to the Bishop Tuff.  more » « less
Award ID(s):
1855751
NSF-PAR ID:
10354017
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Journal of Petrology
Volume:
62
Issue:
12
ISSN:
0022-3530
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Many lines of evidence from high P–T experiments, thermodynamic models, and natural observations suggest that slab-derived aqueous fluids, which flux mantle wedges contain variable amounts of dissolved carbon. However, constraints on the effects of H2O–CO2 fluids on mantle melting, particularly at mantle wedge P–T conditions, are limited. Here, we present new piston cylinder experiments on fertile and depleted peridotite compositions with 3.5 wt.% H2O and XCO2 [= molar CO2 / (CO2 + H2O)] of 0.04–0.17. Experiments were performed at 2–3 GPa and 1350°C to assess how temperature, peridotite fertility, and XCO2 of slab-derived fluid affects partial melting in mantle wedges. All experiments produce olivine + orthopyroxene +7 to 41 wt.% partial melt. Our new data, along with previous lower temperature data, show that as mantle wedge temperature increases, primary melts become richer in SiO2, FeO*, and MgO and poorer CaO, Al2O3, and alkalis when influenced by H2O–CO2 fluids. At constant P–T and bulk H2O content, the extent of melting in the mantle wedge is largely controlled by peridotite fertility and XCO2 of slab-fluid. High XCO2 depleted compositions generate ~7 wt.% melt, whereas, at identical P–T, low XCO2 fertile compositions generate ~30 to 40 wt.% melt. Additionally, peridotite fertility and XCO2 have significant effects on peridotite partial melt compositions. At a constant P–T–XCO2, fertile peridotites generate melts richer in CaO and Al2O3 and poorer in SiO2, MgO + FeO, and alkalis. Similar to previous experimental studies, at a constant P–T fertility condition, as XCO2 increases, SiO2 and CaO of melts systematically decrease and increase, respectively. Such distinctive effects of oxidized form of dissolved carbon on peridotite partial melt compositions are not observed if the carbon-bearing fluid is reduced, such as CH4-bearing. Considering the large effect of XCO2 on melt SiO2 and CaO concentrations and the relatively oxidized nature of arc magmas, we compare the SiO2/CaO of our experimental melts and melts from previous peridotite + H2O ± CO2 studies to the SiO2/CaO systematics of primitive arc basalts and ultra-calcic, silica-undersaturated arc melt inclusions. From this comparison, we demonstrate that across most P–T–fertility conditions predicted for mantle wedges, partial melts from bulk compositions with XCO2 ≥ 0.11 have lower SiO2/CaO than all primitive arc melts found globally, even when correcting for olivine fractionation, whereas partial melts from bulk compositions with XCO2 = 0.04 overlap the lower end of the SiO2/CaO field defined by natural data. These results suggest that the upper XCO2 limit of slab-fluids influencing primary arc magma formation is 0.04 < XCO2 < 0.11, and this upper limit is likely to apply globally. Lastly, we show that the anomalous SiO2/CaO and CaO/Al2O3 signatures observed in ultra-calcic arc melt inclusions can be reproduced by partial melting of either CO2-bearing hydrous fertile and depleted peridotites with 0 < XCO2 < 0.11 at 2–3 GPa, or from nominally CO2-free hydrous fertile peridotites at P > 3 GPa.

     
    more » « less
  2. The granitic water-saturated solidus (G-WSS) is the lower temperature limit of magmatic mineral crystallization. The accepted water-saturated solidus for granitic compositions was largely determined >60 years ago1. More recent advances in experimental petrology, improved analytical techniques, and recent observations that granitic systems can remain active or spend a significant proportion of their lives at conditions below the traditional G-WSS2–5 necessitate a careful experimental investigation of the near-solidus regions of granitic systems. Natural and synthetic starting materials were melted at 10 kbar and 900°C with 48 wt% H2O to produce hydrous glasses for subsequent experiments at lower PT conditions used to locate the G-WSS. We performed crystallization experiments and melting experiments at temperatures ranging from 575 to 800°C and 1, 6, 8, and 10 kbar on 12 granitoid compositions. First, we ran a series of isothermal crystallization experiments along each isobar at progressively lower temperatures until runs completely crystallized to identify apparent solidus temperatures. Geochemical analyses of quenched glass compositions demonstrate that progressive crystallization drives all starting compositions towards silica-rich, water-saturated rhyolitic/granitic melts (e.g., ~7578 wt% SiO2). After identifying the apparent solidus temperatures at which the various compositions crystallized, we then ran series of reversal-type melting experiments. With the goal of producing rocks with hydrous equilibrium microstructures, we crystallized compositions at temperatures ~10°C below the apparent solidus identified in crystallization experiments, and then heated isobarically to conditions that produced ~20% melt during the crystallization experiments. Importantly, crystallization experiments and heating experiments at the same PT conditions produced similar proportions of melt, crystals, and vapor. A time-series of experiments 230 days at PT conditions previously identified to produce ~10% to 20% melt did not reveal any kinetic effects on melt crystallization. Experiments at 6 to 10 kbar crystallized/melted at temperatures close to the published G-WSS. However, at lower pressures where the published G-WSS is strongly curved in PT space, all compositions investigated contained melt to temperatures ~75 to 100°C below the accepted G-WSS. The similarity of crystallization temperatures for the higher-pressure experiments to previously published results, similar phase proportions in melting and crystallization experiments, and the lack of kinetic effects on crystallization collectively suggest that our lower pressure constraints on the G-WSS are accurate. The new experimental results demonstrating that the lower-pressure G-WSS is significantly lower than unanimously accepted estimates will help us to better understand the storage conditions, evolution, and potential for eruption in mid- to upper-crustal silicic magmatic systems. (1) Tuttle, O.; Bowen, N. Origin of Granite in the Light of Experimental Studies in the System NaAlSi3O8–KAlSi3O8–SiO2–H2O; Geological Society of America Memoirs; Geological Society of America, 1958; Vol. 74. https://doi.org/10.1130/MEM74. (2) Rubin, A. E.; Cooper, K. M.; Till, C. B.; Kent, A. J. R.; Costa, F.; Bose, M.; Gravley, D.; Deering, C.; Cole, J. Rapid Cooling and Cold Storage in a Silicic Magma Reservoir Recorded in Individual Crystals. Science 2017, 356 (6343), 1154–1156. https://doi.org/10.1126/science.aam8720. (3) Andersen, N. L.; Jicha, B. R.; Singer, B. S.; Hildreth, W. Incremental Heating of Bishop Tuff Sanidine Reveals Preeruptive Radiogenic Ar and Rapid Remobilization from Cold Storage. Proceedings of the National Academy of Sciences 2017, 114 (47), 12407–12412. https://doi.org/10.1073/pnas.1709581114. (4) Ackerson, M. R.; Mysen, B. O.; Tailby, N. D.; Watson, E. B. Low-Temperature Crystallization of Granites and the Implications for Crustal Magmatism. Nature 2018, 559 (7712), 94–97. https://doi.org/10.1038/s41586-018-0264-2. (5) Glazner, A. F.; Bartley, J. M.; Coleman, D. S.; Lindgren, K. Aplite Diking and Infiltration: A Differentiation Mechanism Restricted to Plutonic Rocks. Contributions to Mineralogy and Petrology 2020, 175 (4). https://doi.org/10.1007/s00410-020-01677-1. 
    more » « less
  3. Volcano inflation, for durations of months to years immediately following an eruption, has been observed at a number of volcanoes, including the 2011/12 eruption of Cordón Caulle, Chile. Such reinflation is often explained by replenishment of the magma reservoir from a deeper source. Whether and why that is the case remains uncertain in most instances, but the implications for renewed eruptive potential may be profound. Here, we posit redistribution of melt within a zoned magma reservoir consisting of a crystal rich mush overlain by an eruptible layer of crystal poor rhyolite as an alternate mechanism for reinflation. Such a zoned magma body is consistent with conceptual models for how crystal poor rhyolites form and with the presence of mafic enclaves within the Cordón Caulle rhyolite. The enclaves can be interpreted as pieces of mush entrained into the overlying rhyolite during its withdrawal from the reservoir. We test the hypothesis that melt from the inter-crystalline pores of the mush can redistribute by porous flow into the overlying crystal poor rhyolite, causing inflation after an eruption. We simulate the flow of melt within the zoned reservoir during and after eruption with a numerical model. As crystal poor rhyolite is erupted, magma pressure within the rhyolite layer above the mush decreases. Consequently, interstitial melt flows upward within the mush, toward the reduced pressure at the interface of mush and crystal poor rhyolite. The mush is treated as a poroelastic material, with interstitial melt flow governed by Darcy's law. Thus, the change in pressure caused by withdrawal from the overlying rhyolite diffuses downward into the mush as the interstitial melt flows upward. The change in pore pressure results in an elastic deformation of the mush matrix. Because pore pressure diffusivity is small, melt redistribution can persist for years after eruptive activity ends, leading to slow inflation compared to fast eruptive deflation. We predict a partial recovery of volume lost from the eruption. Reinflation occurs because the expansion of decompressing melt flowing from the mush into the crystal poor rhyolite exceeds compression of the poroelastic mush. For cases where the interstitial melt is moderately compressible due to exsolved volatiles, our model reproduces the deformation observed at Cordón Caulle. 
    more » « less
  4. Abstract

    Interpreting unrest at silicic volcanoes requires knowledge of the magma storage conditions and dynamics that precede eruptions. The Laguna del Maule volcanic field, Chile, has erupted ~40 km3of rhyolite over the last 20 ka. Astonishing rates of sustained surface inflation at >25 cm/year for >12 years reveal a large, restless system. Integration of geochronologic, petrologic, geomorphic, and geophysical observations provides an unusually rich context to interpret ongoing and prehistoric processes. We present new volatile (H2O, CO2, S, F, and Cl), trace element, and major element concentrations from 109 melt inclusions hosted in quartz, plagioclase, and olivine from seven eruptions. Silicic melts contain up to 8.0 wt. % H2O and 570 ppm CO2. In rhyolites melt inclusions track decompression‐driven fractional crystallization as magma ascended from ~14 to 4 km. This mirrors teleseismic tomography and magnetotelluric findings that reveal a domain containing partial melt spanning from 14 to 4 km. Ce and Cl contents of rhyolites support the generation of compositionally distinct domains of eruptible rhyolite within the larger reservoir. Heat, volatiles, and melt derived from episodic mafic recharge likely incubate and grow the shallow reservoir. Olivine‐hosted melt inclusions in mafic tephra contain up to 2.5 wt. % H2O and 1,140 ppm CO2and proxy for the volatile load delivered via recharge into the base of the silicic mush at ~14 to 8 km. We propose that mafic recharge flushes deeper reaches of the magma reservoir with CO2that propels H2O exsolution, upward accumulation of fluid, pressurization, and triggering of rhyolitic eruptions.

     
    more » « less
  5. Abstract

    Quartz‐hosted melt inclusions are windows into the assembly and storage of compositionally zoned magmas systems like the giant magma body responsible for the 2,800 km3Youngest Toba Tuff (YTT), Indonesia. Feldspar‐sensitive element concentrations and interelement ratios in YTT quartz‐hosted melt inclusions are found to cluster into three discrete melt inclusion populations. Each population represents a small fraction of the overall range in YTT melt compositions and no more than ∼12% to ∼21% fractional crystallization at mostly eutectoid conditions accounts for the trace elements variations within them. Kinships between host matrix glasses and co‐hosted phases are interpreted as evidence of spatially discrete chemical domains of magma within a more broadly zoned magma system. Quartz growth in the dominant portion of the system apparently occurred after extraction of melt from mush, whereas chemically distinct crystal cores in the more crystal‐rich, low‐silica magmas may have been derived from mush during extraction or remobilized by mush rejuvenation. Collectively, the chemical zoning of the YTT appears to have developed prior to most quartz growth, likely “bottom‐up” due to mush‐derived melt heterogeneity.

     
    more » « less