skip to main content


Title: Using Trace Elements in Olivine Phenocrysts to Test Models of Melt Differentiation at Convergent Margins: Constraints from Times Series in Monogenetic Volcanoes in the Central Transmexican Volcanic Belt
Models of subduction zone magmatism ascribe the andesitic composition of arc magmas to crustal processes, such as crustal assimilation and/or fractional crystallization, that basaltic mantle melts experience during their ascent through the upper plate crust. However, results from time series study of olivine-phyric high-Nb basalts and basaltic andesites from two monogenetic arc volcanoes (V. Chichinautzin and Texcal Flow) that are constructed on the ~45 km thick continental basement of the central Transmexican Volcanic Belt (TMVB) are inconsistent with this model. Instead, ratios of radiogenic isotope and incompatible trace elements suggest that these volcanoes were constructed through multiple individual melt batches ascending from a progressively changing mantle source. Moreover, the high Ni contents of the olivine phenocrysts, together with their high mantle-like 3He/4Heoliv =7-8 Ra with high crustal δ18O oliv = +5.5 to +6.5‰ (n=12) point to the presence of secondary ‘reaction pyroxenites’ in the mantle source that create primary silicic arc magmas through melt-rock reaction processes in the mantle [1, 2] . Here we present additional trace element concentration of the high-Ni olivines by electron microprobe (Mn, Ca) and laser-ablation ICPMS (Li, Cr and V) analysis in order to test this model. Olivine Li (2-7 ppm) and Mn (1170- 2810 ppm) increase with decreasing fosterite (Fo89 to Fo75), while Cr (29-364 ppm), V (4-11 ppm) and Ca (825-2390 ppm) decrease. Quantitative modeling shows that these trends in their entirety cannot be controlled by fractional crystallization under variable melt water H2O or oxygen fugacity (fO2), or co-crystallization of Cr-spinel. Instead, the variations support the existence of compositionally distinct melt batches during earliest melt evolution. Moreover, the trace element trends are qualitatively consistent with a model of progressive source depletion by serial melting (shown in olivine Ca, V and Cr) that is triggered by the repetitive addition of silicic slab components (shown by olivine Li). These findings suggest mantle source variations are not eliminated despite the thick crust these magmas pass during ascent. [1] Straub et al. (2013) J Petrol 54 (4): 665-701; [2] Straub et al. (2015) Geochim Cosmochim Acta 166: 29-52.  more » « less
Award ID(s):
1921624
NSF-PAR ID:
10175574
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
AGU Fall Meeting 2019
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    The oxygen fugacity (fO2) of the Earth’s upper mantle and its melting products is an important parameter in the geochemical evolution of arc magmas and their connection with the continental crustal construction and growth. Several works have focused on the fO2 of peridotite xenoliths, primitive melts in relatively young arc settings, and mid-ocean ridge basalts (MORB) but few studies have attempted to examine the early redox history of primitive magmas in mature arcs. Hence, our understanding of the nature and evolution of fO2 during the subduction cycle remains limited. Here, we investigate the basaltic tephra from the Los Hornitos monogenetic cones in central-southern Chile, which are among the most primitive materials reported in the Southern Andes (olivine Mg#  92.5, and Ni  5000 mgg1). These features offer a unique opportunity to explore the fO2 conditions below the Andean arc by studying olivine phenocrysts and their contained crystal and melt inclusions. We integrated EPMA, LA-ICP-MS, and m-XANES analyses to constrain the redox conditions recorded in the basaltic tephra by three different and self-reliant methods. First, we determined the fO2 based on the olivine-spinel equilibrium, yielding average values DFMQ + 1.3 ± 0.4 (1r). Second, we constrained the fO2 conditions of melt inclusions using Fe m-XANES data and the redox dependent olivine-melt vanadium partitioning. After correcting for post-entrapment crystallization and diffusive iron loss, the Fe m-XANES data indicate that the melt inclusions were trapped in average at DFMQ +2.5 ± 0.5 (1r). Results using the olivine-melt vanadium partitioning oxybarometer in melt inclusions are in agreement with Fe m-XANES data, yielding average DFMQ values of +2.6 ± 0.3 (1r). In order to test the potential effects of other postentrapment modifications of the melt inclusions that could have affected the fO2 prior to eruption, we assessed the residence time of these magmas using Mg-Fe interdiffusion modelling in olivine. The short residence times (<200 days) compared to vanadium re-equilibration models strongly suggest that the melt inclusions preserve the prevailing fO2 conditions during their entrapment. Correlations between melt inclusions major element composition and their fO2 determined by Fe m-XANES, as well as V/Sc modelling reveal a case of post-melting oxidation of the LHC magmas. We argue that primitive arc magmas behave as an open system with respect to fO2 during their early geochemical evolution. Our data indicate a complex fO2 early history of primitive melts in the southern Andes and provide a cautionary note on the direct extrapolation of primitive melts fO2 values to that of their mantle source. 
    more » « less
  2. Abstract

    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. Constraining the volatile content of magmas is critical to our understanding of eruptive processes and their deep Earth cycling essential to planetary habitability [R. Dasgupta, M. M. Hirschmann, Earth Planet. Sci. Lett. 298 , 1 (2010)]. Yet, much of the work thus far on magmatic volatiles has been dedicated to understanding their cycling through subduction zones. Further, studies of intraplate mafic volcanism have disproportionately focused on Hawaii [P. E. Wieser et al., Geochem. Geophys. Geosyst. 22 , e2020GC009364 (2021)], making assessments of the overall role of intraplate volcanoes in the global volatile cycles a challenge. Additionally, while mafic volcanoes are the most common landform on Earth and the Solar System [C. A. Wood, J. Volcanol. Geotherm. Res. 7 , 387–413 (1980)], they tend to be overlooked in favor of silicic volcanoes when it comes to their potential for explosivity. Here, we report primitive (olivine-hosted, with host Magnesium number – Mg# 78 to 88%) melt inclusion (MI) data from Fogo volcano, Cabo Verde, that suggest that oceanic intraplate silica-undersaturated explosive eruptions sample volatile-rich sources. Primitive MI (melt Mg# 70 to 71%) data suggest that these melts are oxidized (NiNiO to NiNiO+1) and very high in volatiles (up to 2 wt% CO 2 , 2.8 wt% H 2 O, 6,000 ppm S, 1,900 ppm F, and 1,100 ppm Cl) making Fogo a global endmember. Storage depths calculated from these high volatile contents also imply that magma storage at Fogo occurs at mantle depths (~20 to 30 km) and that these eruptions are fed from the mantle. Our results suggest that oceanic intraplate mafic eruptions are sustained from the mantle by high volatile concentrations inherited from their source and that deep CO 2 exsolution (here up to ~800 MPa) drives their ascent and explosivity. 
    more » « less
  4. 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
  5. Abstract

    Siliceous slab-derived partial melts infiltrate the sub-arc mantle and cause rock-melt reactions, which govern the formation of diverse primary arc magmas and lithological heterogeneities. The effect of bulk water content, composition of reactants, and nature of melt infiltration (porous versus channelized) on the rock-melt reactions at sub-arc conditions have been investigated by previous studies. However, the effect of multiple episodes of rock-melt reactions in such scenarios has not been investigated before. Here, we explore mantle wedge modifications through serial additions of hydrous-silicic slab partial melts and whether such a process may ultimately explain the origin of high-Mg# andesites found in arcs worldwide. A series of piston-cylinder experiments simulate a serial addition of silicic slab melts in up to three stages (I through III) at 3 GPa and 800–1050°C, using rock-melt proportions of 75–25 and 50–50. A synthetic KLB-1 and a natural rhyolite (JR-1) represented the mantle and the slab components, respectively. Right from the first rock-melt interaction, the peridotite mantle transforms into olivine-free mica-rich pyroxenites ± amphibole ± quartz/coesite in equilibrium with rhyolitic-hydrous melts (72–80 wt% SiO2 and 40–90 Mg#). The formation of olivine-free pyroxenite seems to be controlled by complex functions of T, P, rock-melt ratio, wedge composition, and silica activity of the slab-melt. Remarkably, the pyroxenites approach a melt-buffered state with progressive stages of rock-melt reactions, where those rhyolitic melts inherit and preserve the major (alkalis, Fe, Mg, Ca) and trace element slab-signature. Our results demonstrate that lithological heterogeneities such as pyroxenites formed as products of rock-melt reactions in the sub-arc mantle may function as melt ‘enablers,’ implying that they may act as pathways that enable the infiltrating melt to retain their slab signature without undergoing modification. Moreover, the density contrast between the products of rock-melt reaction (melts and residues) and the average mantle wedge (~150 to 400 kg/m3) may help forming instabilities and diapiric rise of the slab components into the mantle wedge. However, the fate of the primitive slab-melts seems to be associated with the length of the pathway of mantle interaction which explains the evident wide magma spectrum as well as their degree of slab garnet-signature dilution. This work and the existence of high-Mg# Mexican-trondhjemites indicates that almost pristine slab-melts can make their way up to crustal levels and contribute to the arc magma diversity.

     
    more » « less