skip to main content

Title: ReversePetrogen : A Multiphase Dry Reverse Fractional Crystallization‐Mantle Melting Thermobarometer Applied to 13,589 Mid‐Ocean Ridge Basalt Glasses

We present a new algorithm,ReversePetrogen (RevPet), to infer mantle melting conditions (pressure, temperature, source composition) using evolved basalts that have experienced multiphase fractional crystallization.RevPetmeasures and minimizes the compositional distance between experimentally predicted phase saturation boundaries and an erupted basalt and the more primitive liquids that return it to a primary melt. We useRevPetto investigate mantle melting conditions at mid‐ocean ridges (MORs) using a global data set of 13,589 basaltic glasses. We find that their average apparent mantle potential temperature (TP*) is 1322°C ± 56°C with melting pressures of 13.0 ± 5 kbars. Inferring the true initial (pre‐melted)TPfromTP* requires knowing the style and degree of melting of the input basalts. If MORB glasses were entirely produced by near‐fractional melting of a homogenous source, they would record the cooling of the mantle during melting from an initialTP = ∼1380°C (ΔTP = 0°C) down toTP = ∼1270°C. If, instead, they were all fully pooled near‐fractional melts of the same source, they would record variations in ambient MORTPfrom ∼1300°C to 1450°C (ΔTP = 150°C). However, because MOR basalts are thought to be both near‐fractional and variably pooled melts of variable sources, MORTPmust be intermediate between these two extremes. Our best estimate, consistent with MOR crustal thickness, is that ambient MORTPis homogenous (∼1350°C–1400°C) except near hotspots whereTPreaches ∼1600°C. Some primitive glasses found near slow‐spreading ridges and back‐arcs record very low temperatures (TP* < 1250°C) and pressures of melting (<10 kbar) and reflect mantle cooling during melting and melt equilibration in the mantle lithosphere.

more » « less
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The Mado Megamullion is an oceanic core complex (OCC) in the Shikoku back‐arc basin within the Philippine Sea Plate. Mantle peridotites (serpentinized) recovered by six dredge and submersible cruises exhibit signatures of extensive deformation. Amorphous pseudomorphs after plagioclase in many of the samples, as well as plagioclase‐spinel intergrowths, are clear evidence of melt stagnation and mantle reaction. Spinels show a wide range of compositions in terms of their Cr#, Mg#, and TiO2content. The presence of apparently magmatic high‐temperature pargasitic amphibole in veins and as replacement of clinopyroxene suggests that it may be a primary or near‐primary mineral crystallized from a hydrous melt which is unusual for abyssal peridotites. Two trace‐element populations of clinopyroxenes are in equilibrium with depleted and enriched basaltic melts, respectively. Rare‐earth element (REE) in the most depleted clinopyroxenes are modeled by 10% fractional melting except for a ubiquitous La‐Ce “kick.” Multiple models of open system melting combined with subsequent mixing of an enriched melt can explain the REE data. Broadly it appears that the peridotites underwent variable degrees of partial melting with moderate influx of enriched melts, which agrees with the other textural and chemical evidence of melt‐rock reaction and re‐fertilization. The compositions of the accumulated melts simulated by the open system models reproduce the enrichments in fluid mobile elements (Ba, U, and Pb) observed in basalts dredged from the Shikoku basin. Back‐arc basin peridotites at Mado Megamullion appear to have a unique petrographic and geochemical character that is distinct from those of peridotites exposed at the seafloor after formation from mid‐ocean ridges.

    more » « less
  2. null (Ed.)
    Oceanic island basalts are targeted for geochemical study because they provide a direct window into mantle composition and a wealth of information on the dynamics and timescales associated with Earth mixing. Previous studies mainly focused on the shield volcanic stage of oceanic islands and the more fusible, enriched mantle components that are easily distinguished in those basalts. Mantle depleted compositions are typically more difficult to resolve unless large amounts of this material participated in mantle melting (e.g., mid-ocean ridges), or unique processes allow for their compositions to be erupted undiluted, such as very small degrees of melting of a source with minimal fusible enriched components (e.g., rejuvenated basalts) or as xenoliths (e.g., abyssal peridotites). Mantle depleted components, defined here as material with low time-integrated Rb/Sr (low 87Sr/86Sr) and high time-integrated Sm/Nd and Lu/Hf ratios (high 143Nd/144Nd and 176Hf/177Hf) relative to primitive mantle, derive from a potentially very large volume reservoir (up to 80% of the mantle), and therefore need adequate characterization in order estimate the composition of the Earth and mantle-derived melts. This review focuses on mantle depleted compositions in oceanic island basalts using the Hawaiian-Emperor chain as a case study. The Hawaiian-Emperor chain is the ∼6000 km long geological record of the deeply sourced Hawaiian mantle plume, active for>81 Myr. Hawaiian volcanism evolves through four volcanic stages as a volcano traverses the Hawaiian plume: alkalic preshield, tholeiitic shield (80–90% volcano volume), alkalic postshield (∼1%), and silica undersaturated rejuvenated (< 0.1%). We report Pb-Sr-Nd-Hf isotope compositions and trace element concentrations of three rejuvenated Northwest Hawaiian Ridge basalts and compare them to an exhaustive compiled dataset of basalts from the Hawaiian Islands to the Emperor Seamounts. The Northwest Hawaiian Ridge (NWHR) includes 51 volcanoes spanning ∼42 m.y. between the bend in the Hawaiian-Emperor chain and the Hawaiian Islands where there is no high-precision isotopic data published on the rejuvenated-stage over ∼47% of the chain. NWHR and Hawaiian Island rejuvenated basalts are geochemically similar, indicating a consistent source for rejuvenated volcanism over ∼12.5 million years. In contrast, shield-stage basalts from the oldest Emperor Seamounts are more depleted in isotopic composition (i.e., higher 176Hf/177Hf, and 143Nd/144Nd with lower 87Sr/86Sr and 208Pb*/206Pb*) and trace element concentrations (i.e., much lower concentrations of highly incompatible elements) than all other depleted Hawaiian basalts younger than the bend, including NWHR rejuvenated basalts. The strongly depleted source for the oldest Emperor Seamounts (> 70 Ma) was likely related to interaction with the Kula-Pacific-Izanagi mid-ocean ridge spreading system active near the Hawaiian plume in the Late Cretaceous. In contrast, the incompatible trace element ratios of NWHR rejuvenated basalts require a distinct source in the Hawaiian mantle plume that was imprinted by ancient (> 1 Ga) partial melting, likely ancient recycled oceanic lithosphere. This review of the geochemistry of Hawaiian depleted components documents the need for the sampling of multiple distinctive depleted compositions, each preferentially melted during specific periods of Hawaiian plume activity. This suggests that the composition of depleted components can evolve during the lifetime of the mantle plume, as observed for enriched components in the Hawaiian mantle plume. Changes in the composition of depleted components are dominantly controlled by the upper mantle tectonic configurations at the time of eruption (i.e., proximity to a mid-ocean ridge), as this effect overwhelms the signal imparted by potentially sampling different lower mantle components through time. 
    more » « less
  3. Abstract

    Young mafic lavas from the East African Western Rift record melting of subcontinental lithospheric mantle that was metasomatically modified by multiple tectonic events. We report new isotope data from monogenetic cinder cones near Bufumbira, Uganda, in the Virunga Volcanic Field:87Sr/86Sr = 0.7059–0.7079,εNd = −6.5 to −1.3,εHf = −6.3 to +0.9,208Pb/204Pb = 40.1–40.7,207Pb/204Pb = 15.68–15.75, and206Pb/204Pb = 19.27–19.45. Olivine phenocrysts from the Bufumbira lavas have3He/4He = 6.0–7.4RA. The isotopic data, in conjunction with major and trace element systematics, indicate that primitive Bufumbira magmas are derived from two different metasomatized lithospheric source domains. Melts generated by lower degrees of melting record greater contributions from ∼1 to 2 Ga isotopically enriched garnet‐amphibole‐phlogopite pyroxenite veins within the lithosphere. As melting progresses, these vein melts become increasingly diluted by melts that originate near the lithosphere/asthenosphere boundary, shifting the isotopic compositions toward the common lithospheric mantle (CLM) proposed by Furman and Graham (1999, This ∼450–500 Ma source domain appears to underlie all Western Rift volcanic provinces and is characterized by87Sr/86Sr ∼ 0.705,εNd∼ 0,εHf∼ +1 to +3,206Pb/204Pb ∼ 19.0–19.2,208Pb/204Pb ∼ 39.7, and3He/4He ∼ 7RA. Basal portions of the dense subcontinental lithospheric mantle may become gravitationally unstable and founder into underlying warmer asthenosphere, exposing surfaces where melting of locally heterogeneous veins produces small‐volume, alkaline mafic melts. Mafic lavas from all Western Rift volcanic provinces record mixing between the CLM and locally variable metasomatized source domains, suggesting this style of melt generation is fundamental to the development of magma‐poor rifts.

    more » « less
  4. 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
  5. Abstract

    The Izu‐Bonin‐Mariana (IBM) fore arc preserves igneous rock assemblages that formed during subduction initiation circa 52 Ma. International Ocean Discovery Program (IODP) Expedition 352 cored four sites in the fore arc near the Ogasawara Plateau in order to document the magmatic response to subduction initiation and the physical, petrologic, and chemical stratigraphy of a nascent subduction zone. Two of these sites (U1440 and U1441) are underlain by fore‐arc basalt (FAB). FABs have mid‐ocean ridge basalt (MORB)‐like compositions, however, FAB are consistently lower in the high‐field strength elements (TiO2, P2O5, Zr) and Ni compared to MORB, with Na2O at the low end of the MORB field and FeO* at the high end. Almost all FABs are light rare earth element depleted, with low total REE, and have low ratios of highly incompatible to less incompatible elements (Ti/V, Zr/Y, Ce/Yb, and Zr/Sm) relative to MORB. Chemostratigraphic trends in Hole U1440B are consistent with the uppermost lavas forming off axis, whereas the lower lavas formed beneath a spreading center axis. Axial magma of U1440B becomes more fractionated upsection; overlying off‐axis magmas return to more primitive compositions. Melt models require a two‐stage process, with early garnet field melts extracted prior to later spinel field melts, with up to 23% melting to form the most depleted compositions. Mantle equilibration temperatures are higher than normal MORB (1,400 °C–1,480 °C) at relatively low pressures (1–2 GPa), which may reflect an influence of the Manus plume during subduction initiation. Our data support previous models of FAB origin by decompression melting but imply a source more depleted than normal MORB source mantle.

    more » « less