More Like this

1. Temporal changes in diamond formation by subduction throughout Earth history: thermal modeling, seismology, and petrology evidence

   Shirey, Steven B; Smit, Karen V; Pearson, D Graham; Walter, Michael J; Stachel, Thomas (July 2024, Journal of International Kimberlite Conference Abstracts (JIKCA))

   In the current plate tectonic regime, thermal modeling, petrology, and seismology show that subsurface portions of cold slabs carry some of their volatiles into the deep upper mantle, mantle transition zone, and uppermost lower mantle avoiding the devolatilization occurring with normal arc and wedge subduction. Slab crustal remnants at these depths can melt by intersecting their carbonated solidus whereas slab mantle remnants can devolatilize by warming and metamorphosing to ‘dryer’ mineral assemblages. Since fluid release and earthquake production (“dehydration embrittlement”) operates down to ~300 km depths in all subduction zones, we propose, that deep-focus earthquakes trace the places of fluid release at deeper levels (350 to 750 km). Fluids in faults related to earthquake generation will become diamond-forming as they react with mantle rocks along the fault walls. Diamonds thus formed will record deformation produced by mantle convection and slab buckling during mantle storage. Lithospheric diamonds, stored in static ancient continental keels, lack the connection to this type of geodynamic regime that is evident for sublithospheric diamonds. However, a comparison between the two diamond types suggests a geologic model for lithospheric diamond formation in the ancient past. Lithospheric diamonds and sublithospheric diamonds both contain evidence for the recycling of sediments or surficial rocks that have equilibrated at low temperatures with seawater. The known way to inject these materials into diamond-forming regions is slab subduction. Hence both diamond types may have formed by variants of this same process that differ in depth and style over geologic time. Lithospheric diamonds are different from sublithospheric diamonds in critical ways: higher average N content, ages extending into the Paleoarchean, inclusion assemblages indicating formation at lower pressure, and lack of ubiquitous deformation features. Nitrogen content is the key to relating lithospheric diamonds to the subducting slab. Nitrogen occurs in clays and sediments at the slab surface or uppermost crust. Regardless of whether the slab is hot or cold during subduction, nitrogen will be removed into a mantle wedge if one exists. Additionally, diamonds will not survive in the melts/fluids generated in the wedge under oxidizing conditions. For sublithospheric diamonds, their low to non-existent nitrogen content occurs because they are derived from slab fluids/melts once nitrogen has been largely removed or from rocks deeper in the slab where nitrogen is scarce. The much higher nitrogen in lithospheric diamonds suggests that they formed from fluids/melts derived near the slab surface that contained N. In the Archean, such slabs must have subducted close to the nascent mantle keel with no mantle wedge so the fluids could be directly reduced by the mantle keel. We propose a gradual temporal change from shallow, keel-adjacent, mantle-wedge-poor subduction that produced lithospheric diamonds starting in the Paleoarchean to wedge-avoiding, cold and deep subduction that produced sublithospheric diamonds in the Paleozoic. This temporal change is consistent with many geologic features: an early stagnant lid and a buoyant Archean oceanic lithosphere; the slab-imbrication, advective thickening, and diamond-richness of portions of mantle keels; and anomalously diamond-rich ancient eclogites. 

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

   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
3. DiadFit: An open-source Python3 tool for peak fitting of Raman data from silicate melts and CO2 fluids

   https://doi.org/10.30909/vol.07.01.335359  

   Wieser, Penny E; DeVitre, Charlotte (January 2024, Volcanica)

   We present DiadFit—an open-source Python3 tool for efficient processing of Raman spectroscopy data collected from fluid inclusions, melt inclusions and silicate melts. DiadFit is optimized to fit the characteristic peaks from CO2 fluids (Fermi diads, hot bands, 13C), gas species such as SO2, N2, solid precipitates (e.g. carbonates), and Ne emission lines with easily tweakable background positions and peak shapes. DiadFit's peak fitting functions are used as part of a number of workflows optimized for quantification of CO2 in melt inclusion vapour bubbles and fluid inclusions. DiadFit can also convert between temperature, pressure and density using various CO2 and CO2-H2O equations of state (EOS), allowing calculation of fluid inclusion pressures (and depths in the crust), conversion of homogenization temperatures from microthermometry to CO2 density, and propagation of uncertainties associated with EOS calculations using Monte Carlo methods. There are also functions to quantify the area ratio of the silicate vs. H2O region of spectra collected on silicate glasses to determine H2O contents in glasses and melt inclusions. Documentation and worked examples are available (https://bit.ly/DiadFitRTD, https://bit.ly/DiadFitYouTube). 

   more » « less
4. Molecular H2 in silicate melts

   https://doi.org/10.1016/j.gca.2024.10.020  

   Foustoukos, Dionysis I (January 2025, Geochimica et Cosmochimica Acta)

   A series of hydrothermal diamond anvil cell experiments was conducted to constrain the equilibrium distribution of molecular H2 between H2O-saturated sodium aluminosilicate melts and H2O at elevated temperatures (600 – 800 °C) and pressures (317 – 1265 MPa). The distribution of H2 between the silicate liquid and the aqueous fluid was achieved through real-time monitoring of the H-H stretching vibration under in situ conditions using Raman vibrational spectroscopy. Results show that the solubility of H2 in silicate melts saturated with H2O decreases as the temperature increases, with control exerted by the mole fraction of H2O in the melt. The dissolution of H2 in the hydrous silicate melts appears to follow Henrian behavior, resembling that of an inert, neutral non-polar species. The implications are extended into developing an understanding of the H partitioning between H2-rich atmospheres blanketing magma oceans in the early history of planetary bodies. For example, transferring H from primordial atmospheric envelopes to the interior of rocky exoplanets may be less efficient than previously believed, which should be considered in models of volatile retention. Experimental data also suggest that minimal amounts of solar nebula H2 are likely to dissolve in the molten surface of primitive objects in the protoplanetary disk (~10-5 to 10-9 mole faction H2 in melt), contradicting the highly reducing conditions observed in chondrule mineral compositions. 

   more » « less
5. Slab Transport of Fluids to Deep Focus Earthquake Depths—Thermal Modeling Constraints and Evidence From Diamonds

   https://doi.org/10.1029/2020AV000304  

   Shirey, Steven B.; Wagner, Lara S.; Walter, Michael J.; Pearson, D. Graham; van Keken, Peter E. (May 2021, AGU Advances)

   Abstract The nature and cause of deep earthquakes remain enduring unknowns in the field of seismology. We present new models of thermal structures of subducted slabs traced to mantle transition zone depths that permit a detailed comparison between slab pressure/temperature (P/T) paths and hydrated/carbonated mineral phase relations. We find a remarkable correlation between slabs capable of transporting water to transition zone depths in dense hydrous magnesium silicates with slabs that produce seismicity below ∼300‐km depth, primarily between 500 and 700 km. This depth range also coincides with theP/Tconditions at which oceanic crustal lithologies in cold slabs are predicted to intersect the carbonate‐bearing basalt solidus to produce carbonatitic melts. Both forms of fluid evolution are well represented by sublithospheric diamonds whose inclusions record the existence of melts, fluids, or supercritical liquids derived from hydrated or carbonate‐bearing slabs at depths (∼300–700 km) generally coincident with deep‐focus earthquakes. We propose that the hydrous and carbonated fluids released from subducted slabs at these depths lead to fluid‐triggered seismicity, fluid migration, diamond precipitation, and inclusion crystallization. Deep focus earthquake hypocenters could track the general region of deep fluid release, migration, and diamond formation in the mantle. The thermal modeling of slabs in the mantle and the correlation between sublithospheric diamonds, deep focus earthquakes, and slabs at depth demonstrate a deep subduction pathway to the mantle transition zone for carbon and volatiles that bypasses shallower decarbonation and dehydration processes. 

   more » « less