More Like this

1. The 3.1 Ma to 100 ka evolution of the Goat Rocks Volcanic Complex: Persistence and evolution of magmatism at a long-lived major andesite locus on the Cascades Arc

   Wall, Kellie T. ( July 2022 , Theses and dissertations Oregon State University)

   PhD Dissertation Abstract: The imposing andesite stratovolcano is the characteristic expression of subduction zone magmatism, posing hazards to coastal populations and bearing insight into deep Earth processes. On a map of a typical volcanic arc, one can easily distinguish the approximately linear alignment and regular spacing of these major edifices that stand out from a diffuse distribution of mafic volcanoes (e.g. the Quaternary Cascades; Hildreth, 2007). The andesitic composite volcanoes have a reputation for being complex, open systems: crystal zoning “stratigraphies,” diverse crystal cargoes including antecrysts or xenocrysts, quenched magmatic inclusions, and variations in isotopic signatures are among the many lines of evidence that these systems involve a variety of igneous processes and melt sources. To investigate the development and evolution of such transcrustal magma factories, I have conducted a detailed temporal, spatial, and geochemical characterization of a long-lived arc volcanic center in the southern Washington Cascades, the Goat Rocks volcanic complex. Results from ⁴⁰Ar/³⁹Ar and U/Pb geochronology constrain the lifespan of the Goat Rocks volcanic complex from ~3.1 Ma to ~100 ka. During this time, four major composite volcanoes were built (as well as several smaller volcanoes). From oldest to youngest, these are Tieton Peak, Bear Creek Mountain, Lake Creek volcano, and Old Snowy Mountain. Four volcanic stages are defined based on the lifespans of these centers and distinct compositional changes that occur from one to the next: Tieton Peak stage (3.1-2.6 Ma), Bear Creek Mountain stage (1.6-1.1 Ma), Lake Creek stage (1.1 Ma to 456 ka), and Old Snowy Mountain stage (440 ka to 115 ka). Two lava flow remnants also have ages in the interim between Tieton Peak stage and Bear Creek Mountain stage (2.3 Ma and 2.1 Ma), and their sources are not yet identified. The ages of the Bear Creek Mountain and Lake Creek stages in fact overlap, and the gap between Lake Creek stage and Old Snowy Mountain stage is only on the order of 10⁴ years. Based on supporting compositional evidence, the Bear Creek Mountain, Lake Creek, and Old Snowy Mountain stage volcanoes are considered to be the migrating surface expressions of a continuous magmatic system that was active over at least ~1.5 million years. It remains uncertain whether the gaps between the Tieton Peak stage, scattered early Pleistocene andesites, and Bear Creek Mountain stage are due to incomplete exposure/sampling or real quiescent periods earlier in the development of the Goat Rocks volcanic complex. Throughout the construction of the andesitic complex, mafic volcanoes were active on its periphery. These include the Miriam Creek volcano (~3.6-3.1 Ma), Devils Washbasin volcano (3.0-2.7 Ma), Hogback Mountain (1.1 Ma – 891 ka), Lakeview Mountain (194 ka), and Walupt Lake volcano (65 ka). Two basalt and basaltic andesite units (Qob₁ and Qob₂, 1.4 and 1.3 Ma; Hammond, 2017) also erupted from the Goat Rocks area, likely an older incarnation of Hogback Mountain. The suite of mafic magmas erupted in this region are all calcalkaline basalt (or basaltic andesite; CAB), but two compositional groups emerge from the trace element and isotopic data. Group 1 is LILE and LREE-enriched, with higher ⁸⁷Sr/⁸⁶Sr isotopes, and includes compositions from Devils Washbasin, Lower Hogback Mountain, and Lakeview Mountain. Group 2 is less enriched in LILE and LREE and lower in ⁸⁷Sr/⁸⁶Sr, and includes the compositions of Miriam Creek, Qob1, Upper Hogback Mountain, Walupt Lake, and Coleman Weedpatch. The two groups are recurrent through time and with no geographic distinction; in fact, both types were tapped by the Hogback Mountain volcano. Together both of these groups, alongside CABs from Mount Adams and various basalts from Mount St. Helens, form a compositional array between the basalts of the High Cascades and the intraplate-type basalts (IPB) of Mount Adams and Simcoe volcanic field. These results lead to three conclusions. 1) Variably subduction-modified mantle is distributed across the region, perhaps either as stratified layers or a web-like network of fluid pathways amongst less metasomatized mantle. 2) Transitional compositions between the IPBs and typical “High Cascades” CAB/HAOT signature suggest a broader influence of the mantle domain that feeds IPBs—if asthenospheric mantle through a slab window, as suggested by Mullen et al. (2017), then perhaps it bleeds in smaller quantities over a broader area. This compositional trend solidifies the interpretation of the southern Washington Cascades as a unique and coherent “segment” of the arc (the Washington segment of Pitcher and Kent, 2019). 3) The recurrence of variable mafic magma types through time, and with no geographic boundaries, indicates that the compositional evolution of the Goat Rocks volcanic complex was not likely driven by a change in mafic input. Indeed, the Sr, Nd, Hf, and Pb isotope ratios of the intermediate to felsic suite are closely aligned with the local basalts and suggest a limited role of crustal assimilation. Importantly, several mineral thermometers (zircon, ilmenite-magnetite pairs, and amphibole) align in recording higher crystallization temperatures in Bear Creek Mountain to early Lake Creek time, a cooling trend through the Lake Creek stage, and a more diverse range of temperatures in the transition to Old Snowy Mountain stage. Thus, it is suggested that the compositional evolution at Goat Rocks represents a thermal cycle of waxing and waning magmatic flux: where the period of Bear Creek Mountain to early Lake Creek volcanism was the climactic phase of a vertically extensive magma homogenization factory, then the system waned and cooled, ultimately losing its ability to filter, homogenize, and enrich magmas. 

   more » « less
2. The Systematics of Olivine CaO + Cr-Spinel in High-Mg# Arc Volcanic Rocks: Evidence for in-Situ Mantle Wedge Depletion at the Arc Volcanic Front

   https://doi.org/10.1093/petrology/egad085  

   Straub, Susanne M. ; Batanova, Valentina ; Sobolev, Alexander ; Gómez-Tuena, Arturo ; Espinasa-Perena, Ramon ; Fleming, W. Lindsey ; Bindeman, Ilya N. ; Stuart, Finlay M. ; Widom, Elisabeth ; Iizuka, Yoshiyuki ( December 2023 , Journal of Petrology)

   We investigated the state of the arc background mantle (i.e. mantle wedge without slab component) by means of olivine CaO and its Cr-spinel inclusions in a series of high-Mg# volcanic rocks from the Quaternary Trans-Mexican Volcanic Belt. Olivine CaO was paired with the Cr# [molar Cr/(Cr + Al) *100] of Cr-spinel inclusions, and 337 olivine+Cr-spinel pairs were obtained from 33 calc-alkaline, high-K and OIB-type arc front volcanic rocks, and three monogenetic rear-arc basalts that lack subduction signatures. Olivine+Cr-spinels display coherent elemental and He–O isotopic systematics that contrast with the compositional diversity of the bulk rocks. All arc front olivines have low CaO (0.135 ± 0.029 wt %) relative to rear-arc olivines which have the higher CaO (0.248 ± 0.028 wt %) of olivines from mid-ocean ridge basalts. Olivine 3He/4He–δ18O isotope systematics confirm that the olivine+Cr-spinels are not, or negligibly, affected by crustal basement contamination, and thus preserve compositional characteristics of primary arc magmas. Variations in melt H2O contents in the arc front series and the decoupling of olivine CaO and Ni are inconsistent with controls on the olivine CaO by melt water and/or secondary mantle pyroxenites. Instead, we propose that low olivine CaO reflects the typical low melt CaO of high-Mg# arc magmas erupting through thick crust. We interpret the inverse correlation of olivine CaO and Cr-spinel Cr# over a broad range of Cr# (~10–70) as co-variations of CaO, Al and Cr of their (near) primary host melts, which derived from a mantle that has been variably depleted by slab-flux driven serial melt extraction. Our results obviate the need for advecting depleted residual mantle from rear- and back-arc region, but do not upset the larger underlying global variations of melt CaO high-Mg# arc magmas worldwide, despite leading to considerable regional variations of melt CaO at the arc front of the Trans-Mexican Volcanic Belt.

    

   more » « less
3. 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

   https://doi.org/10.1093/petrology/egab087  

   Jolles, Jameson S ; Lange, Rebecca A ( December 2021 , Journal of Petrology)

   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
4. Application of an Improved Olivine-Melt Thermometer/Hygrometer to the Colima Cone Basanites and Minettes of Western Mexico: Implications for the Mantle Source of Unusually High-MgO Melts

   https://doi.org/10.1093/petrology/egad064  

   Pu, Xiaofei ; Lange, Rebecca A. ; Moore, Gordon M. ( August 2023 , Journal of Petrology)

   A collection of quaternary, high-MgO (≤13.4 wt%) basanite and minette cinder and lava cones, with an enhanced arc geochemical signature, are located along the northern margin of the N–S Colima rift in western Mexico. The Colima rift overlies the lithospheric suture between the Jalisco block and Guerrero terrane, as well as the tear between the Rivera and Cocos subducting oceanic plates. From the literature, volatile analyses of olivine-hosted melt inclusions in the Colima cone samples document notably high concentrations of dissolved H2O in the melt (≤ 7.0 wt%) as well as degassing-induced phenocryst growth over a range of depths ≤25 km. In this study, it is shown that the high-MgO character of the Colima suite reflects liquid compositions, consistent with evidence for their rapid transit to the surface, without stalling in a crustal magma chamber. The most Mg-rich olivine analyzed in each sample matches the equilibrium composition at the liquidus based on olivine-melt Mn–Mg and Fe2+–Mg exchange coefficients. Application of both a Mg- and Ni-based olivine-melt thermometer, calibrated on the same experimental data set, to the Colima cone suite provides the temperature and dissolved H2O content at the liquidus. Because the Ni thermometer is insensitive to dissolved H2O in the melt, it gives the actual temperature at the onset of olivine phenocryst growth. For the nine Colima samples that range from 13.4–9.2 wt% MgO, resulting temperatures range from 1221°C to 1056°C (± 6–11°C). In contrast, the Mg thermometer is sensitive to dissolved H2O in the melt, and its application (without a correction of H2O) gives the temperature of olivine crystallization under anhydrous conditions. When the Mg- and Ni-based temperatures are paired, the depression of the liquidus (∆T = TMg–TNi) due to dissolved H2O can be obtained. For the high-MgO (>9 wt%) Colima samples, ∆T values range from 188°C to 109°C. Corrections for the effect of pressure (i.e. evidence that phenocryst growth began at ~700 MPa), increase ∆T by ~21°C. An updated model calibration (on experiments from the literature) that relates ∆T with dissolved H2O in the melt shows that the best fit (R2 = 0.95) is linear, wt% H2O = 0.047*∆T, with a standard error of ±0.5 wt%. Although the experimental data set spans a wide range of melt composition (e.g. 47–58 wt% SiO2, 4.4–10.2 wt% MgO, 1.3–4.9 wt% Na2O, 0.1–5.0 wt% K2O, 0.3–5.3 wt% H2O), no dependence on anhydrous melt composition is resolved. Application of this updated model to the Colima suite gives H2O contents of 5.1–8.8 wt% H2O, consistent with those analyzed in olivine-hosted MIs from the literature. When the thermometry and hygrometry results for the Colima cone suite are compared to those for the adjacent calc–alkaline basalts from the Tancítaro Volcanic Field (TVF) in Michoacán, two distinct linear trends in a plot of wt% H2O vs. temperature are found, indicative of different mantle sources. It is proposed that the high-MgO (>11 wt%) Colima cone melts were derived from a phlogopite-bearing harzburgitic mantle at the base of the Jalisco block lithosphere, whereas both TVF and Colima melts with ≤10 wt% MgO were derived from the asthenosphere (i.e. arc mantle wedge). In both mantle sources, slab-derived fluids were an important source of H2O.

    

   more » « less
5. Effects of H2O–CO2 Fluids, Temperature, and Peridotite Fertility on Partial Melting in Mantle Wedges and Generation of Primary Arc Basalts

   https://doi.org/10.1093/petrology/egad047  

   Lara, Michael ; Dasgupta, Rajdeep ( July 2023 , Journal of Petrology)

   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