Subduction is a component of plate tectonics, which is widely accepted as having operated in a manner similar to the present-day back through the Phanerozoic Eon. However, whether Earth always had ...plate tectonics or, if not, when and how a globally linked network of narrow plate boundaries emerged are matters of ongoing debate. Earth's mantle may have been as much as 200-300 °C warmer in the Mesoarchean compared to the present day, which potentially required an alternative tectonic regime during part or all of the Archean Eon. Here we use a data set of the pressure (P), temperature (T), and age of metamorphic rocks from 564 localities that vary in age from the Paleoarchean to the Cenozoic to evaluate the petrogenesis and secular change of metamorphic rocks associated with subduction and collisional orogenesis at convergent plate boundaries. Based on the thermobaric ratio (T/P), metamorphic rocks are classified into three natural groups: high T/P type (T/P > 775 °C/GPa, mean T/P ∼1105 °C/GPa), intermediate T/P type (T/P between 775 and 375 °C/GPa, mean T/P ∼575 °C/GPa), and low T/P type (T/P < 375 °C/GPa, mean T/P ∼255 °C/GPa). With reference to published thermal models of active subduction, we show that low T/P oceanic metamorphic rocks preserving peak pressures >2.5 GPa equilibrated at P-T conditions similar to those modeled for the uppermost oceanic crust in a wide range of active subduction environments. By contrast, those that have peak pressures <2.2 GPa may require exhumation under relatively warm conditions, which may indicate subduction of young oceanic lithosphere or exhumation during the initial stages of subduction. However, low T/P oceanic metamorphic rocks with peak pressures of 2.5-2.2 GPa were exhumed from depths where, in models of active subduction, the slab and overriding plate change from being decoupled (at lower P) to coupled (at higher P), possibly suggesting a causal relationship. In relation to secular change, the widespread appearance of low T/P metamorphism in the Neoproterozoic represents a 'modern' style of cold collision and deep slab breakoff, whereas rare occurrences of low T/P metamorphism in the Paleoproterozoic may reveal atypical localized regions of cold collision. Low T/P metamorphism is not known from the Archean geological record, but the absence of blueschists in particular is unlikely to reflect secular change in the composition of the oceanic crust. In addition, the premise that the formation of lawsonite requires abnormally low thermal gradients and the postulate that oceanic subduction-related rocks register significantly lower maximum pressures than do continental subduction-related rocks, and imply different mechanisms of exhumation, are not supported. The widespread appearance of intermediate T/P and high T/P metamorphism at the beginning of the Neoarchean, and the subsequent development of a clear bimodality in tectono-thermal environments are interpreted to be evidence of the stabilization of subduction during a transition to a globally linked network of narrow plate boundaries and the emergence of plate tectonics.
The zinc (Zn) stable isotope system has great potential for tracing planetary formation and differentiation processes due to its chalcophile, lithophile and moderately volatile character. As an ...initial approach, the terrestrial mantle, and by inference, the bulk silicate Earth (BSE), have previously been suggested to have an average δ66Zn value of ∼+0.28‰ (relative to JMC 3-0749L) primarily based on oceanic basalts. Nevertheless, data for mantle peridotites are relatively scarce and it remains unclear whether Zn isotopes are fractionated during mantle melting. To address this issue, we report high-precision (±0.04‰; 2SD) Zn isotope data for well-characterized peridotites (n=47) from cratonic and orogenic settings, as well as their mineral separates. Basalts including mid-ocean ridge basalts (MORB) and ocean island basalts (OIB) were also measured to avoid inter-laboratory bias. The MORB analyzed have homogeneous δ66Zn values of +0.28±0.03‰ (here and throughout the text, errors are given as 2SD), similar to those of OIB obtained in this study and in the literature (+0.31±0.09‰). Excluding the metasomatized peridotites that exhibit a wide δ66Zn range of −0.44‰ to +0.42‰, the non-metasomatized peridotites have relatively uniform δ66Zn value of +0.18±0.06‰, which is lighter than both MORB and OIB. This difference suggests a small but detectable Zn isotope fractionation (∼0.1‰) during mantle partial melting. The magnitude of inter-mineral fractionation between olivine and pyroxene is, on average, close to zero, but spinels are always isotopically heavier than coexisting olivines (Δ66ZnSpl-Ol=+0.12±0.07‰) due to the stiffer Zn-O bonds in spinel than silicate minerals (Ol, Opx and Cpx). Zinc concentrations in spinels are 11–88 times higher than those in silicate minerals, and our modelling suggests that spinel consumption during mantle melting plays a key role in generating high Zn concentrations and heavy Zn isotopic compositions of MORB. Therefore, preferential melting of spinel in the peridotites may account for the Zn isotopic difference between spinel peridotites and basalts. By contrast, the absence of Zn isotope fractionation between silicate minerals suggests that Zn isotopes are not significantly fractionated during partial melting of spinel-free garnet-facies mantle. If the studied non-metasomatized peridotites represent the refractory upper mantle, mass balance calculation shows that the depleted MORB mantle (DMM) has a δ66Zn value of +0.20±0.05‰ (2SD), which is lighter than the primitive upper mantle (PUM) estimated in previous studies (+0.28±0.05‰, 2SD, Chen et al., 2013b; +0.30±0.07‰, 2SD, Doucet et al., 2016). This indicates that the Earth’s upper mantle has a heterogeneous Zn isotopic composition vertically, which is probably due to shallow mantle melting processes.
Knowledge about the properties of silicate melts is needed by volcanologists and petrologists to evaluate the dynamics of volcanic eruptions and magmatic processes. These properties include the ...solubility and diffusivity of volatile components in silicate melts, silicate melt viscosity, and the fragmentation condition. Data and models of each property are reviewed and assessed. For rhyolitic melts many properties are sufficiently well known to allow realistic modeling of volcanic and magmatic processes. One interesting example is the role of speciation in the solubility and diffusivity of H2O and CO2. Even though both H2O and CO2 are present in silicate melts as at least two species, the complexity in the solubility and diffusion behavior of H2O and the simplicity of CO2 are due to differences in the speciation reaction: For the H2O component the stoichiometric coefficient is one for one hydrous species (molecular H2O) but is two for the other hydrous species (OH) in the species interconversion reaction, whereas for CO2 the stoichiometric coefficients for all carbon species are one. The investigation of the species reaction not only helps in understanding the solubility and diffusion behavior, but the reaction among the hydrous species also serves as a geospeedometer (cooling rate indicator) for hydrous rhyolitic pyroclasts and glass and provides a method to infer viscosity. For melts other than rhyolite, a preliminary description of their properties is also available, but much more experimental and modeling work is necessary to quantify these properties more accurately.
The Neoarchean Gadwal greenstone belt in the eastern Dharwar craton, India, hosts a well preserved metavolcanic sequence that is dominated by tholeiitic and calc-alkaline ...basalt-andesite-dacite-rhyolite series, which includes boninitic geochemical varieties. Bulk-rock Lu–Hf and Sm–Nd isotope systematics of these apparently arc-related volcanic rocks yield indistinguishable ages of 2.701±0.024Ga and 2.702±0.026Ga, respectively. On the basis of the close spatial association and identical ages of the different rock types we suggest 2.70±0.03Ga as the age of crystallization of the different rock types within the Gadwal metavolcanic sequence. In contrast, bulk-rock Pb–Pb isotope systematics of the same samples yield a significantly younger and less precise age of 2.466Ga (+0.068/−0.110Ga). We tentatively interpret this younger age to represent a metallogenic and crustal reworking event in the Dharwar craton, which disturbed the U–Pb system but not the Lu–Hf or Sm–Nd systems. The Gadwal metavolcanic rocks have positive initial εHf(2.70Ga)=+1.6 to+8.7 and slightly negative to positive εNd(2.70Ga)=−0.1 to+3.0 values, consistent with an origin from a long term depleted source relative to a chondritic reservoir at ∼2.7Ga. Lack of correlation between initial isotopic compositions and major or trace element indices of fractionation and alteration suggest that the observed isotope variability probably reflects compositional variation in the Gadwal source, similar to that observed in modern day island arcs. Two boninitic samples of the Gadwal sequence have εHf∼8.3 and 8.7, and are more radiogenic than average depleted mantle for the time period 3.2 to 2.5Ga (εHf=4 to 6). Early (perhaps Hadean) differentiation events that led to a depleted and heterogeneous mantle are apparent in the Nd and Hf isotope systematics of 3.7–3.8Ga Isua supracrustal rocks. The radiogenic Hf isotopes of the Gadwal boninites and the Hf, Nd isotope systematics of rocks from other locations in the 3.4 to 2.5Ga time period are consistent with the survival of fragments of an early depleted mantle later in the Archean. From ∼2.0Ga to present, the time-integrated 176Lu/177Hf and 147Sm/144Nd of the depleted mantle appears nearly constant and similar to the present day average MORB source. These data indicate that progressive elimination of early (>4.5Ga) formed heterogeneities in the depleted mantle dominated the history of the Archean mantle, and that portions of early depleted reservoirs survived through the Mesoarchean. These results have implications for the mixing scales for the early terrestrial mantle and the timing of the initiation of present day plate tectonics.
Alkaline magmatism is an important chemical end-member of magmatic activity that typically occurs in response to small volume melting of asthenospheric- and/or lithospheric mantle material in ...intra-continental settings. Understanding trace element partitioning and phase equilibria during alkaline magmatism can therefore provide constraints on intra-continental geodynamic settings. However, the partitioning of trace elements between alkaline melts and their dominant equilibrium mineral phases remains poorly constrained. Feldspathoids in particular have received limited attention with regards to their trace element contents, hampering our ability to interpret geochemical trends in alkaline magmatic systems. In this study, we performed a series of 1 atmosphere experiments in a gas-mixing furnace using a variety of highly alkaline (Na2O + K2O = 4.15–14.97 wt%) and silica-undersaturated (SiO2 = 36.73–45.96 wt%) lava compositions from Nyiragongo, Democratic Republic of Congo, in order to investigate the partitioning behaviour of trace elements in minerals from alkaline magmas. Experimental runs were performed with oxygen fugacity buffered at both QFM (quartz-fayalite-magnetite equilibrium) and QFM + 1 and cover a range of geologically-relevant temperatures (1025–1200 °C). The quenched products of these experiments contained leucite, nepheline, melilite, clinopyroxene, olivine, and rhönite crystals, of which glass-crystal pairs were analysed for rare earth elements, large-ion lithophile elements, and high-field-strength elements. Leucite and nepheline host considerable quantities of large-ion lithophile elements but take up negligible amounts of more highly charged cations. Åkermanitic melilite readily incorporates mono- to trivalent cations with a preference for light over heavy rare earth elements, but incorporates only select divalent cations. Rhönite and clinopyroxene have analogous partitioning behaviours, with a strong preference for heavy over light rare earth elements. Fractionation modelling using the reported partitioning behaviours reproduces the 2021 eruption products of Nyiragongo, with 48% fractionation from an olivine-melilitic parental melt composition. Crystallization of trace-element poor feldspathoid amplifies pre-existing high LREE/MREE ratios of the parental magma and progressively increase trace element abundances for all but monovalent cations.
This paper introduces Petrolog3, software for modeling (1) fractional and equilibrium crystallization, (2) reverse fractional crystallization at variable pressure, melt oxidation state and melt H2O ...contents, and (3) postentrapment reequilibration of melt inclusions in olivine. Petrolog3 offers an algorithm that allows calculations with a potentially unlimited number of (1) mineral‐melt equilibrium models for major and trace elements and (2) models describing melt physical parameters such as density and viscosity, melt oxidation state, and solubility of fluid components in silicate melts. The current version of the software incorporates 46 mineral‐melt equilibrium models for 8 minerals; 3 models describing distribution of trace elements between minerals and melt; 4 models of melt oxidation state; 1 model for H2O solubility in silicate melts; and 4 models describing melt density and viscosity. The idea behind the program is to provide the community of igneous petrologists and geochemists with a user‐friendly interface for using any combinations of available mineral‐melt equilibrium models for computer simulation of the crystallization process.
Key Points
New algorithm for crystallization modeling
User friendly software interface
Ability to compare available models
Abstract
Serpentinization is a metamorphic process that can stabilize highly reduced hydrogen-rich fluids. Previous measurements of elevated CH4 and H2 concentrations in ultramafic-hosted submarine ...springs indicate that active serpentinization occurs along mid-ocean ridge systems at seafloor pressures (∼<500 bar) and temperatures (∼<350 °C). Serpentinites also exist at higher pressures in subduction zones; for example, during retrograde hydration of the forearc mantle wedge and during prograde deserpentinization within the subducted slab. However, many studies demonstrating the thermodynamic stability of reduced serpentinite fluids have been limited to terrestrial seafloor conditions. To investigate the redox state of serpentinite fluids at elevated pressures, phase equilibria and fluid compositions were computed for 100–700 °C and 1–20 kbar using aqueous silica activity (aSiO2(aq)) as a governing parameter. Silica-sensitive, redox-buffering assemblages were selected to be consistent with previously proposed reactions: SiO2(aq)–fayalite–magnetite (QFM), SiO2(aq)–Fe-brucite–cronstedtite, SiO2(aq)–Fe-brucite–Fe3+-serpentine, plus the silica-free buffer Fe-brucite–magnetite. Fluid species are limited to simple, zerovalent compounds. For silica-bearing redox reactions, aSiO2(aq) is buffered by coexisting ultramafic mineral assemblages in the system MgO–SiO2–H2O. Silica activity and fO2 are directly correlated, with the most reduced fluids stabilized by the least siliceous assemblages. Silica activity and fO2 increase with pressure, but are more strongly dependent on temperature, leading to greater silica enrichment and more oxidized conditions along shallow, warm subduction paths than along steeper, colder paths. Reduced fluids with mCH4/mCO2 > 1 and fO2 below QFM are present only when serpentine is stable, and are favored along all subduction trajectories except shallow P–T paths at eclogite-grade. Values of mH2 and mCO/mCO2 depend strongly on P and T, but also on the choice of redox buffer, especially whether the Fe-serpentine component is cronstedtite or Fe3+-serpentine. Methane and H2S production are thermodynamically favored throughout the P–T range of the serpentinized forearc mantle and in other settings with similar conditions; for example, deep planetary seafloors. The model offers a generalized technique for estimating the redox state of a fluid-saturated serpentinite at elevated P and T, and yields results consistent with previous petrographic and thermodynamic analyses. High-pressure serpentinization may be an important source of reduced species that could influence prebiotic chemistry, support microbial life in the deep biosphere or in deep planetary oceans, or promote greenhouse warming on early Earth.