Kimberlites reveal 2.5-billion-year evolution of a deep, isolated mantle reservoir

Primordial neon and the deep mantle origin of kimberlites

The genesis of kimberlites is unclear despite the economic and scientific interest surrounding these diamond-bearing magmas. One critical question is whether they tap ancient, deep mantle domains or the shallow convecting mantle with partial melting triggered by plumes or plate tectonics. To address this question, we report the He-Ne-Ar isotopic compositions of magmatic fluids trapped in olivine from kimberlites worldwide. The kimberlites which have been least affected by addition of deeply subducted or metasomatic components have Ne isotopes less nucleogenic than the upper mantle, hence requiring a deep-mantle origin. This is corroborated by previous evidence of small negative W isotope anomalies and kimberlite location along age-progressive hot-spot tracks. The lack of strong primordial He isotope signatures indicates overprinting by lithospheric and crustal components, which suggests that Ne isotopes are more robust tracers of deep-mantle contributions in intraplate continental magmas. The most geochemically depleted kimberlites may preserve deep remnants of early-Earth heterogeneities.

Global mantle perturbations following the onset of modern plate tectonics

Plate tectonics drives the compositional diversity of Earth's convecting mantle through subduction of lithosphere. In this context, the role of evolving global geodynamics and plate (re)organization on the spatial and temporal distribution of compositional heterogeneities in the convecting mantle is poorly understood. Here, using the geochemical compositions of intracontinental basalts formed over the past billion years, we show that intracontinental basalts with subchondritic initial neodymium-144/neodymium-143 values become common only after 300 million years, broadly coeval with the global appearance of kimberlites with geochemically enriched isotopic signatures. These step changes in the sources of intraplate magmatism stem from a rapid increase in the supply of deeply subducted lithosphere during the protracted formation of Pangea following the widespread onset of "modern" (cold and deep) subduction in the late Neoproterozoic. We argue that the delay (~300 million years) in the appearance of enriched intraplate magmas reflects the time required for the sinking and (re)incorporation of slabs into the sources of mantle-derived magmas.

A common precursor for global hotspot lavas

Hotspot lavas exhibit chemical heterogeneity, much of which is ascribed to heterogeneous deep mantle sources that contain various components with distinct composition, origin and age. However, characterizing primary melt compositions and mantle heterogeneity directly is challenging. Here we investigate a global dataset of hotspot lavas to constrain the incompatible-element composition of their parental melts and sources. Trace-element ratios indicate that the compositional heterogeneity of global hotspot lavas is not primary, but reflects processes that hotspot melts undergo as they ascend to the surface. We find the parental melts of these lavas, as well as of kimberlites and basalts from large igneous provinces, to be uniform in their elemental, and radiogenic and noble-gas isotope, composition. We suggest that the parental melts to all of these lavas derive from a depleted and outgassed mantle reservoir that was replenished with incompatible element-enriched material during the Archaean. This interpretation explains the elemental, radiogenic and noble-gas isotope compositions of hotspot lavas without requiring a heterogeneous lower mantle or the long-term survival of undegassed relics from a primordial Earth.

Kimberlite eruptions driven by slab flux and subduction angle

Kimberlites are sourced from thermochemical upwellings which can transport diamonds to the surface of the crust. The majority of kimberlites preserved at the Earth's surface erupted between 250 and 50 million years ago, and have been attributed to changes in plate velocity or mantle plumes. However, these mechanisms fail to explain the presence of strong subduction signatures observed in some Cretaceous kimberlites. This raises the question whether there is a subduction process that unifies our understanding of the timing of kimberlite eruptions. We develop a novel formulation for calculating subduction angle based on trench migration, convergence rate, slab thickness and density to connect the influx of slab material into the mantle with the timing of kimberlite eruptions. We find that subduction angles combined with peaks in slab flux predict pulses of kimberlite eruptions. High rates of subducting slab material trigger mantle return flow that stimulates fertile reservoirs in the mantle. These convective instabilities transport slab-influenced melt to the surface at a distance inbound from the trench corresponding to the subduction angle. Our deep-time slab dip formulation has numerous potential applications including modelling the deep carbon and water cycles, and an improved understanding of subduction-related mineral deposits.

Perturbation of the deep-Earth carbon cycle in response to the Cambrian Explosion

Earth's carbon cycle is strongly influenced by subduction of sedimentary material into the mantle. The composition of the sedimentary subduction flux has changed considerably over Earth's history, but the impact of these changes on the mantle carbon cycle is unclear. Here, we show that the carbon isotopes of kimberlite magmas record a fundamental change in their deep-mantle source compositions during the Phanerozoic Eon. The ¹³C/¹²C of kimberlites before ~250 Ma preserves typical mantle values, whereas younger kimberlites exhibit lower and more variable ratios-a switch coincident with a recognized surge in kimberlite magmatism. We attribute these changes to increased deep subduction of organic carbon with low ¹³C/¹²C following the Cambrian Explosion when organic carbon deposition in marine sediments increased significantly. These observations demonstrate that biogeochemical processes at Earth's surface have a profound influence on the deep mantle, revealing an integral link between the deep and shallow carbon cycles.

Light oxygen isotopes in mantle-derived magmas reflect assimilation of sub-continental lithospheric mantle material

Oxygen isotope ratios in mantle-derived magmas that differ from typical mantle values are generally attributed to crustal contamination, deeply subducted crustal material in the mantle source or primordial heterogeneities. Here we provide an alternative view for the origin of light oxygen-isotope signatures in mantle-derived magmas using kimberlites, carbonate-rich magmas that assimilate mantle debris during ascent. Olivine grains in kimberlites are commonly zoned between a mantle-derived core and a magmatic rim, thus constraining the compositions of both mantle wall-rocks and melt phase. Secondary ion mass spectrometry (SIMS) analyses of olivine in worldwide kimberlites show a remarkable correlation between mean oxygen-isotope compositions of cores and rims from mantle-like ¹⁸O/¹⁶O to lower 'crustal' values. This observation indicates that kimberlites entraining low-¹⁸O/¹⁶O olivine xenocrysts are modified by assimilation of low-¹⁸O/¹⁶O sub-continental lithospheric mantle material. Interaction with geochemically-enriched domains of the sub-continental lithospheric mantle can therefore be an important source of apparently 'crustal' signatures in mantle-derived magmas.

Tungsten-182 evidence for an ancient kimberlite source

Globally distributed kimberlites with broadly chondritic initial ¹⁴³Nd-¹⁷⁶Hf isotopic systematics may be derived from a chemically homogenous, relatively primitive mantle source that remained isolated from the convecting mantle for much of the Earth's history. To assess whether this putative reservoir may have preserved remnants of an early Earth process, we report ¹⁸²W/¹⁸⁴W and ¹⁴²Nd/¹⁴⁴Nd data for "primitive" kimberlites from 10 localities worldwide, ranging in age from 1,153 to 89 Ma. Most are characterized by homogeneous μ¹⁸²W and μ¹⁴²Nd values averaging -5.9 ± 3.6 ppm (2SD, n = 13) and +2.7 ± 2.9 ppm (2SD, n = 6), respectively. The remarkably uniform yet modestly negative μ¹⁸²W values, coupled with chondritic to slightly suprachondritic initial ¹⁴³Nd/¹⁴⁴Nd and ¹⁷⁶Hf/¹⁷⁷Hf ratios over a span of nearly 1,000 Mya, provides permissive evidence that these kimberlites were derived from one or more long-lived, early formed mantle reservoirs. Possible causes for negative μ¹⁸²W values among these kimberlites include the transfer of W with low μ¹⁸²W from the core to the mantle source reservoir(s), creation of the source reservoir(s) as a result of early silicate fractionation, or an overabundance of late-accreted materials in the source reservoir(s). By contrast, two younger kimberlites emplaced at 72 and 52 Ma and characterized by distinctly subchondritic initial ¹⁷⁶Hf/¹⁷⁷Hf and ¹⁴³Nd/¹⁴⁴Nd have μ¹⁸²W values consistent with the modern upper mantle. These isotopic compositions may reflect contamination of the ancient kimberlite source by recycled crustal components with μ¹⁸²W ≥ 0.

Remnants of early Earth differentiation in the deepest mantle-derived lavas

The noble gas isotope systematics of ocean island basalts suggest the existence of primordial mantle signatures in the deep mantle. Yet, the isotopic compositions of lithophile elements (Sr, Nd, Hf) in these lavas require derivation from a mantle source that is geochemically depleted by melt extraction rather than primitive. Here, this apparent contradiction is resolved by employing a compilation of the Sr, Nd, and Hf isotope composition of kimberlites-volcanic rocks that originate at great depth beneath continents. This compilation includes kimberlites as old as 2.06 billion years and shows that kimberlites do not derive from a primitive mantle source but sample the same geochemically depleted component (where geochemical depletion refers to ancient melt extraction) common to most oceanic island basalts, previously called PREMA (prevalent mantle) or FOZO (focal zone). Extrapolation of the Nd and Hf isotopic compositions of the kimberlite source to the age of Earth formation yields a ¹⁴³Nd/¹⁴⁴Nd-¹⁷⁶Hf/¹⁷⁷Hf composition within error of chondrite meteorites, which include the likely parent bodies of Earth. This supports a hypothesis where the source of kimberlites and ocean island basalts contains a long-lived component that formed by melt extraction from a domain with chondritic ¹⁴³Nd/¹⁴⁴Nd and ¹⁷⁶Hf/¹⁷⁷Hf shortly after Earth accretion. The geographic distribution of kimberlites containing the PREMA component suggests that these remnants of early Earth differentiation are located in large seismically anomalous regions corresponding to thermochemical piles above the core-mantle boundary. PREMA could have been stored in these structures for most of Earth's history, partially shielded from convective homogenization.

Kimberlite genesis from a common carbonate-rich primary melt modified by lithospheric mantle assimilation

Quantifying the compositional evolution of mantle-derived melts from source to surface is fundamental for constraining the nature of primary melts and deep Earth composition. Despite abundant evidence for interaction between carbonate-rich melts, including diamondiferous kimberlites, and mantle wall rocks en route to surface, the effects of this interaction on melt compositions are poorly constrained. Here, we demonstrate a robust linear correlation between the Mg/Si ratios of kimberlites and their entrained mantle components and between Mg/Fe ratios of mantle-derived olivine cores and magmatic olivine rims in kimberlites worldwide. Combined with numerical modeling, these findings indicate that kimberlite melts with highly variable composition were broadly similar before lithosphere assimilation. This implies that kimberlites worldwide originated by partial melting of compositionally similar convective mantle sources under comparable physical conditions. We conclude that mantle assimilation markedly alters the major element composition of carbonate-rich melts and is a major process in the evolution of mantle-derived magmas.