Barium content of Archaean continental crust reveals the onset of subduction was not global

Earth’s earliest continental crust is dominated by tonalite–trondhjemite–granodiorite (TTG) suites, making these rocks key to unlocking the global geodynamic regime operating during the Archaean (4.0–2.5 billion years ago [Ga]). The tectonic setting of TTG magmatism is controversial, with hypotheses arguing both for and against subduction. Here we conduct petrological modeling over a range of pressure–temperature conditions relevant to the Archaean geothermal gradient. Using an average enriched Archaean basaltic source composition, we predict Ba concentrations in TTG suites, which is difficult to increase after magma generated in the source. The results indicate only low geothermal gradients corresponding to hot subduction zones produce Ba-rich TTG, thus Ba represents a proxy for the onset of subduction. We then identify statistically significant increases in the Ba contents of TTG suites worldwide as recording the diachronous onset of subduction from regional at 4 Ga to globally complete sometime after 2.7 Ga.

Only subduction zone can produce Ba-rich TTG, representing a proxy for the onset of subduction. Statistical increases in Ba contents of Archaean TTGs reveal the diachronous onset of subduction from regional at 4 Ga to globally complete after 2.7 Ga

Fossilized Melts in Mantle Wedge Peridotites

The shallow oxidized asthenosphere may contain a small fraction of potassic silicate melts that are enriched in incompatible trace elements and volatiles. Here, to determine the chemical composition of such melt, we analysed fossilized melt inclusions, preserved as multiphase solid inclusions, from an orogenic garnet peridotite in the Bohemian Massif. Garnet-poor (2 vol.%) peridotite preserves inclusions of carbonated potassic silicate melt within Zn-poor chromite (<400 ppm) in the clinopyroxene-free harzburgite assemblage that equilibrated within the hot mantle wedge (Stage 1, > 1180 °C at 3 GPa). The carbonated potassic silicate melt, which has a major element oxide chemical composition of K₂O = 5.2 wt.%, CaO = 17 wt.%, MgO = 18 wt.%, CO₂ = 22 wt.%, and SiO₂ = 20 wt.%, contains extremely high concentrations of large ion lithophile elements, similar to kimberlite melts. Peridotites cooled down to ≅800 °C during Stage 2, resulted in the growth of garnet relatively poor in pyrope content, molar Mg/(Mg + Fe + Ca + Mn), (ca. 67 mol.%). This garnet displays a sinusoidal REE pattern that formed in equilibrium with carbonatitic fluid. Subsequently, subduction of the peridotite resulted in the formation of garnet with a slightly higher pyrope content (70 mol.%) during the Variscan subduction Stage 3 (950 °C, 2.9 GPa). These data suggest the following scenario for the generation of melt in the mantle wedge. Primarily, infiltration of sediment-derived potassic carbonatite melt into the deep mantle wedge resulted in the growth of phlogopite and carbonate/diamond. Formation of volatile-bearing minerals lowered the density and strength of the peridotite. Finally, phlogopite-bearing carbonated peridotite rose as diapirs in the mantle wedge to form carbonated potassic silicate melts at the base of the overriding lithosphere.

Layering of subcontinental lithospheric mantle

Recent seismic studies reveal a sharp velocity drop mostly at ∼70-100km depth within the thick mantle keel beneath cratons, termed the mid-lithosphere discontinuity (MLD). The common presence of the MLD in cratonic regions indicates structural and property layering of the subcontinental lithospheric mantle (SCLM). The nature and origin of the MLD, and many issues associated with the layering of the SCLM are essential to understand the formation and evolution of continents, and have become frontier subjects in the Earth sciences.