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
Experimental Constraints on Dacite Magma Storage beneath Volcán Quizapu, Chile
Abstract Volcán Quizapu, Chile, is an under-monitored volcano that was the site of two historical eruptions: an effusive eruption in 1846–1847 and a Plinian eruption in 1932, both of which discharged ∼5 km3 (dense rock equivalent) of lava and/or tephra. The majority of material erupted in both cases is trachydacite, nearly identical for each event. We present H2O-saturated, phase equilibrium experiments on this end-member dacite magma, using a pumice sample from the 1932 eruption as the main starting material. At an oxygen fugacity (fO2) of ∼NNO + 0·2 (where NNO is the nickel–nickel oxide buffer), the phase assemblage of An25–30 plagioclase + amphibole + orthopyroxene, without biotite, is stable at 865 ± 10 °C and 110 ± 20 MPa H2O pressure (PH2O), corresponding to ∼4 km depth. At these conditions, experiments also reproduce the quenched glass composition of the starting pumice. At slightly higher PH2O and below 860 °C, biotite joins the equilibrium assemblage. Because biotite is not part of the observed Quizapu phase assemblage, its presence places an upper limit on PH2O. At the determined storage PH2O of ∼110 MPa, H2O undersaturation of the magma with XH2Ofluid = 0·87 would align Ptotal to mineral-based geobarometry estimates of ∼130 MPa. However, XH2Ofluid < 1 is not required to reproduce the Quizapu dacite phase assemblage and compositions. A second suite of experiments at lower fO2 shows that the stability fields of the hydrous silicates (amphibole and biotite) are significantly restricted at NNO – 2 relative to NNO + 0·2. Additional observations of Quizapu lava and pumice samples support the existing hypothesis that rapid pre-eruptive heating drove the effusive 1846–1847 eruption, with important refinements. We demonstrate that microlites in the end-member dacite lavas are consistent with in situ crystallization (during ascent), rather than transfer from an andesite. In one end-member dacite lava, newly identified reverse zoning in orthopyroxene and incipient destabilization of amphibole are consistent with small degrees of heating. Our work articulates a clear direction for future Quizapu studies, which are warranted given the active nature of the Cerro Azul–Descabezado Grande volcanic axis.
more »
« less
- Award ID(s):
- 1347880
- NSF-PAR ID:
- 10273230
- Date Published:
- Journal Name:
- Journal of Petrology
- Volume:
- 62
- Issue:
- 5
- ISSN:
- 0022-3530
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Volcanic eruptions of rhyolitic magma often show shifts from powerful (Vulcanian to Plinian) explosive episodes to a more gentle effusion of viscous lava forming obsidian flows. Another prevailing characteris-tic of these eruptions is the presence of pyroclastic obsidians intermingled with the explosive tephra. This dense, juvenile product is similar to the tephra and obsidian flow in composition, but is generally less degassed than its flow counterpart. The formation mechanism(s) of pyroclastic obsidians and the information they can provide concerning the extent to which magma degassing modulates the eruptive style of rhyolitic eruptions are currently subject to active research. Porous tephra and pyroclastic and flow obsidians from the 1060CE Glass Mountain rhyolitic eruption at Medicine Lake Volcano (California) were analyzed for their porosity, φ, water content, H2O, and hydrogen isotopic composition, δD. H2O in porous pyroclasts is correlated negatively with δD and positively with φ, indicating that the samples were affected by post-eruptive rehydration. Numerical modeling suggests that this rehydration occurred at an average rate of 10−23.5±0.5m2s−1during the ∼960 years since the eruption, causing some pyroclasts to gain up to 1 wt%of meteoric water. Pyroclastic and flow obsidians were not affected by rehydration due to their very low porosity. Comparison between modeled δD-H2O relationships in degassing magma and values measured in the Glass Mountain samples supports the idea that rhyolitic magma degasses in closed-system until its porosity reaches a value of about 65±5%, beyond which degassing occurs in open-system until quench. During the explosive phase, rapidly ascending magma fragments soon after it becomes permeable, creating porous lapilli and ash that continue to degas in open-system within an expanding gas phase. As suggested by recent studies, some ash may aggregate and sinter on the conduit sides at different depths above the fragmentation level, partly equilibrating with the continuously fluxing heavier magmatic vapor, explaining the wide range of H2O contents and high variability in δD measured in the pyroclastic obsidians. Using only H2O and δD, it is impossible to rule out the possibility that pyroclastic obsidians may also form by permeable foam collapse, either syn-explosively near the conduit sides below the fragmentation level or during more effusive periods interspersed in the explosive phase. During the final effusive phase of the eruption, slowly ascending magma degasses in open-system until it reaches the surface, creating flows with low H2O and δD. This study shows that H2O measured in highly porous pyroclasts of a few hundred years or more cannot be used to infer syn-eruptive magma degassing pathways, unless careful assessment of post-eruptive rehydration is first carried out. If their mechanism of formation can be better understood, detailed analysis of the variations in texture and volatile content of pyroclastic obsidians throughout the explosive phase may help decipher the reasons why rhyolitic eruptions commonly shift from explosive to effusive phases.more » « less
-
Abstract Ignimbrite flare-ups are rare periods of intense silicic volcanism during which the pyroclastic volume and eruptive frequency is more than an order of magnitude higher than background activity. Investigating the compositional differences between flare-up and steady-state magmas provides critical constraints on the petrogenetic causes for the event and can offer unique opportunities to investigate the role of large-scale tectonic or geodynamic processes in arc magmatism. In this study, we focus on the bimodal Deschutes Formation ignimbrite flare-up of Central Oregon, which erupted unusually high volumes of pyroclastic material 6.25–5.45 Ma from a new axis of volcanism in the Cascades arc. This episode is marked by increased eruption rates and eruption of more silicic compositions relative to the Quaternary Cascade arc, which rarely erupts rhyolites. Ignimbrites are crystal-poor (<10%) dacite to rhyolites (mostly 65–77 wt.% SiO2) with anhydrous mineral assemblages and higher FeO/MgO, Y, Eu/Eu*, MREE and Zr/Sr, indicating drier magmatic evolution compared to the Quaternary arc, and are more similar to those from the rear-arc High Lava Plains (HLP) province that lies to the east. Magnetite-ilmenite oxybarometry indicates that Deschutes Formation felsic magmas tend to be hotter and more reduced (NNO-1 to NNO) than the Quaternary arc (NNO to NNO + 1.5). Rhyolite-MELTS geobarometry suggests complex storage of diverse Deschutes Formation magmas within the shallow crust (50–250 MPa), and the common co-eruption of multiple plagioclase populations, pumice compositions, and compositionally banded pumice suggest variable degrees of mixing and mingling of distinct magmas. Deschutes magmas also have low δ18Oplagioclase values that indicate partial melting and assimilation of hydrothermally altered shallow crust. Trace element systematics and rhyolite-MELTS modeling suggests that felsic pumice cannot be produced by simple fractionation of co-erupted mafic pumice or basaltic lavas, and requires a crustal melting origin, and trace elements and Pb isotopes suggest that young mafic crust may have been the primary protolith. We suggest that partial melting produced low-Si rhyolite melt (~72 wt.%) that acted as both a parent for the most evolved rhyolites, and as a mixing endmember to create the dacite to rhyodacite magmas with heterogenous plagioclase populations. Unlike the predominantly calc-alkaline basalts erupted in the Quaternary Cascade arc, Deschutes Formation primary basalts are mostly low-K tholeiites, indicative of decompression melting. These are similar to the compositions erupted during a contemporaneous pulse of low-K tholeiite volcanism across the whole HLP that reached into the Cascades rear-arc. We suggest that intra-arc extension focused decompression melts from the back-arc into the arc and that tensional stresses allowed this high flux of hot-dry-reduced basalt throughout the crustal column, causing partial melting of mafic protoliths and the production of hot-dry-reduced rhyolite melts. Depletion of incompatible elements in successive rhyolites implies progressive depletion in fertility of the protolith. Extension also allowed for the establishment of a robust hydrothermal system, and assimilation of hydrothermally-altered rocks by magmas residing in a shallow, complex storage network lead to low δ18O melts. Our findings suggest the integral role that extensional tectonics played in producing an unusual ignimbrite flare-up of hot-dry-reduced rhyolite magmas that are atypical of the Cascades arc and may be an important contributor to flare-ups at arcs worldwide.more » « less
-
Hydrogen is a rapidly diffusing monovalent cation in nominally anhydrous minerals (NAMs, such as olivine, orthopyroxene, and clinopyroxene), which is potentially re-equilibrated during silicate melt-rock and aqueous fluid-rock interactions in massif and abyssal peridotites. We apply a 3D numerical diffusion modeling technique to provide first-order timescales of complete hydrogen re-equilibration in olivine, clinopyroxene, and orthopyroxene over the temperature range 600-1200°C. Model crystals are 1-3 mm along the c-axis and utilize H+ diffusion coefficients appropriate for Fe-bearing systems. Two sets of models were run with different boundary compositions: 1) “low-H models” are constrained by mineral-melt equilibrium partitioning with a basaltic melt that has 0.75 wt% H2O and 2) “high-H models,” which utilize the upper end of the estimated range of mantle water solubility for each phase. Both sets of models yield re-equilibration timescales that are identical and are fast for all phases at a given temperature. These timescales have strong log-linear trends as a function of temperature (R2 from 0.97 to 0.99) that can be used to calculate expected re-equilibration time at a given temperature and grain size. At the high end of the model temperatures (1000-1200°C), H+ completely re-equilibrates in olivine, orthopyroxene, and clinopyroxene within minutes to hours, consistent with previous studies. These short timescales indicate that xenolith NAM mantle water contents are likely to be overprinted prior to eruption. The models also resolve the decoupled water-trace element relationship in Southwest Indian Ridge peridotites, in which peridotite REE abundances are reproduced by partial melting models whereas the relatively high NAM H2O contents require later re-equilibration with melt. At temperatures of 600-800°C, which correspond to conditions of hydrothermal alteration of pyroxene to amphibole and talc, H+ re-equilibration typically occurs over a range of timescales spanning days to years. These durations are well within existing estimates for the duration of fluid flow in oceanic hydrothermal systems, suggesting that peridotite NAM water contents are susceptible to diffusive overprinting during higher temperature hydrothermal alteration. Thus, diffusion during aqueous fluid-rock interactions may also explain NAM H2O contents that are too high to reflect residues of melting. These relatively short timescales at low temperatures suggest that the origin of water contents measured in peridotite NAMs requires additional constraints on sample petrogenesis, including petrographic and trace element analyses. Our 3D model results also hint that H+ may diffuse appreciably during peridotite serpentinization, but diffusion coefficients at low temperature are unconstrained and additional experimental investigations are needed.more » « less
-
Abstract We present microbeam major- and trace-element data from 14 monzodiorites collected from the Malaspina Pluton (Fiordland, New Zealand) with the goal of evaluating processes involved in the production of andesites in lower arc crust. We focus on relict igneous assemblages consisting of plagioclase and amphibole with lesser amounts of clinopyroxene, orthopyroxene, biotite and quartz. These relict igneous assemblages are heterogeneously preserved in the lower crust within sheeted intrusions that display hypersolidus fabrics defined by alignment of unstrained plagioclase and amphibole. Trace-element data from relict igneous amphiboles in these rocks reveal two distinct groups: one relatively enriched in high field strength element concentrations and one relatively depleted. The enriched amphibole group has Zr values in the range of ∼25–110 ppm, Nb values of ∼5–32 ppm, and Th values up to 2·4 ppm. The depleted group, in contrast, shows Zr values <35 ppm and Nb values <0·25 ppm, and Th is generally below the level of detection. Amphibole crystallization temperatures calculated from major elements range from ∼960 to 830 °C for all samples in the pluton; however, we do not observe significant differences in the range of crystallization temperatures between enriched (∼960–840 °C) and depleted groups (∼940–830 °C). Bulk-rock Sr and Nd isotopes are also remarkably homogeneous and show no apparent difference between enriched (εNdi = 0·1 to –0·1; 87Sr/86Sri = 0·70420–0·70413) and depleted groups (εNdi = 0·3 to –0·4; 87Sr/86Sri = 0·70424–0·70411). Calculated amphibole-equilibrium melt compositions using chemometric equations indicate that melts were highly fractionated (molar Mg# <50), andesitic to dacitic in composition, and were much more evolved than bulk lower continental crust or primitive basalts and andesites predicted to have formed from hydrous melting of mantle-wedge peridotite beneath an arc. We suggest that melts originated from a common, isotopically homogeneous source beneath the Malaspina Pluton, and differences between enriched and depleted trace-element groups reflect varying contributions from subducted sediment-derived melt and sediment-derived fluid, respectively. Our data demonstrate that andesites and dacites were the dominant melts that intruded the lower crust, and their compositions mirror middle and upper bulk-continental crust estimates. Continental crust-like geochemical signatures were acquired in the source region from interaction between hydrous mantle-wedge melts and recycled subducted sediment rather than assimilation and/or remelting of pre-existing lower continental crust.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/petrology/egab027
