Cratonic Lithospheric Mantle Research Articles

Multi-stage magmatic history of olivine–leucite lamproite dykes from Banganapalle, Dharwar craton, India: evidence from compositional zoning of spinel

Mesoproterozoic lamproite dykes occurring in the Banganapalle Lamproite Field of southern India show extensive hydrothermal alteration, but preserve fresh spinel, apatite and rutile in the groundmass. Spinels belong to three genetic populations. Spinels of the first population, which form crystal cores with overgrowth rims of later spinels, are Al-rich chromites derived from disaggregated mantle peridotite. Spinels of the second population include spongy-textured grains and alteration rims of titanian magnesian aluminous chromites that formed by metasomatic interactions between mantle wall-rocks and precursor lamproite melts before their entrainment into the erupting lamproite magma. Spinels that crystallised directly from the lamproite magma constitute the third population and show five distinct compositional subtypes (spinel-IIIa to IIIe), which represent discrete stages of crystal growth. First stage magmatic spinel (spinel-IIIa) includes continuously zoned macrocrysts of magnesian aluminous chromite, which formed together with Al–Cr-rich phlogopite macrocrysts from an earlier pulse of lamproite magma at mantle depth. Crystallisation of spinel during the other four identified stages occurred during magma emplacement at crustal levels. Titanian magnesian chromites (spinel-IIIb) form either discrete crystals or overgrowth rims on spinel-IIIa cores. Further generations of overgrowth rims comprise titanian magnesian aluminous chromite (spinel-IIIc), magnetite with ulvospinel component (spinel-IIId) and lastly pure magnetite (spinel-IIIe). Abrupt changes of the compositions between successive zones of magmatic spinel indicate either a hiatus in the crystallisation history or co-crystallisation of other groundmass phases, or possibly magma mixing. This study highlights how different textural and compositional populations of spinel provide important insights into the complex evolution of lamproite magmas including clues to elusive precursor metasomatic events that affect cratonic mantle lithosphere.

Clinopyroxene and Garnet Mantle Cargo in Kimberlites as Probes of Dharwar Craton Architecture and Geotherms, with Implications for Post-1·1 Ga Lithosphere Thinning Events Beneath Southern India

AbstractThe Wajrakarur Kimberlite Field (WKF) on the Eastern Dharwar Craton in southern India hosts several occurrences of Mesoproterozoic kimberlites, lamproites and ultramafic lamprophyres, for which mantle-derived xenoliths are rare and only poorly preserved. The general paucity of mantle cargo has hampered the investigation of the nature and evolution of the continental lithospheric mantle (CLM) beneath cratonic southern India. We present a comprehensive study of the major and trace element compositions of clinopyroxene and garnet xenocrysts recovered from heavy mineral concentrates for three c.1·1 Ga old WKF kimberlite pipes (P7, P9, P10), with the goal to improve our understanding of the cratonic mantle architecture and its evolution beneath southern India. The pressure-temperature conditions recorded by peridotitic clinopyroxene xenocrysts, estimated using single-pyroxene thermobarometry, suggest a relatively moderate cratonic mantle geotherm of 40 mW/m2 at 1·1 Ga. Reconstruction of the vertical distribution of clinopyroxene and garnet xenocrysts, combined with some rare mantle xenoliths data, reveals a compositionally layered CLM structure. Two main lithological horizons are identified and denoted as layer A (∼80–145 km depth) and layer B (∼160–190 km depth). Layer A is dominated by depleted lherzolite with subordinate amounts of pyroxenite, whereas layer B comprises mainly refertilised and Ti-metasomatized peridotite. Harzburgite occurs as a minor lithology in both layers. Eclogite stringers occur within the lower portion of layer A and at the bottom of layer B near the lithosphere–asthenosphere boundary at 1·1 Ga. Refertilisation of layer B is marked by garnet compositions with enrichment in Ca, Ti, Fe, Zr and LREE, although Y is depleted compared to garnet in layer A. Garnet trace element systematics such as Zr/Hf and Ti/Eu indicate that both kimberlitic and carbonatitic melts have interacted with and compositionally overprinted layer B. Progressive changes in the REE systematics of garnet grains with depth record an upward percolation of a continuously evolving metasomatic agent. The intervening zone between layers A and B at ∼145–160 km depth is characterized by a general paucity of garnet. This ‘garnet-paucity’ zone and an overlying type II clinopyroxene-bearing zone (∼115–145 km) appear to be rich in hydrous mineral assemblages of the MARID- or PIC kind. The composite horizon between ∼115–160 km depth may represent the product of intensive melt/rock interaction by which former garnet was largely reacted out and new metasomatic phases such as type II clinopyroxene and phlogopite plus amphibole were introduced. By analogy with better-studied cratons, this ‘metasomatic horizon’ may be a petrological manifestation of a former mid-lithospheric discontinuity at 1·1 Ga. Importantly, the depth interval of the present-day lithosphere–asthenosphere boundary beneath Peninsular India as detected in seismic surveys coincides with this heavily overprinted metasomatic horizon, which suggests that post-1·1 Ga delamination of cratonic mantle lithosphere progressed all the way to mid-lithospheric depth. This finding implies that strongly overprinted metasomatic layers, such as the ‘garnet-paucity’ zone beneath the Dharwar craton, present structural zones of weakness that aid lithosphere detachment and foundering in response to plate tectonic stresses.

Plume activity and carbonated silicate melt metasomatism in Dharwar cratonic lithosphere: Evidence from peridotite xenoliths in Wajrakarur kimberlites

We report petrography, mineralogy, major- and trace-element compositions of a rare selection of spinel- and garnet-bearing peridotite xenoliths and single crystals separated from peridotites hosted in the Mesoproterozoic Wajrakarur kimberlites from the Eastern Dharwar craton (EDC), India. These ultramafic xenoliths consist of olivine (modal 74–82 vol. %) with Fo92–93, clinopyroxene, orthopyroxene, spinel, garnet, and/or ilmenite. Calculated equilibrium pressure and temperature conditions are 2.5–5.0 GPa and 710–1179 °C for these peridotites, which suggests residence depths >160 km near the base of the Dharwar cratonic lithospheric mantle. Garnet in these ultramafic rocks [with Mg# = molar (Mg/(Mg + Fetotal) × 100 of 80–88] displays either “sinuous” LREE-enriched patterns with depletion in Gd and Er for harzburgites or “normal” LREE-depleted, HREE-enriched patterns for lherzolites. Two groups of clinopyroxenes (group-I and group-II) were also observed with high LREE (LaN > 10) and low LREE (LaN < 10), respectively. The Yb vs. Zr contents in peridotitic garnet are characterised by two distinct trends: one with low-T metasomatism for harzburgitic garnets and the other with high-T metasomatism for lherzolitic garnet, which suggests metasomatism from fluids circulating with in the continental lithospheric mantle, resulting in refertilization from harzburgite to lherzolite. REE concentrations of hypothetical melts calculated to be in equilibrium with peridotitic garnet resemble small volume carbonated silicate melts closely similar to natural kimberlite and lamproites from the Dharwar craton. Modeling suggest that olivine-rich peridotites of the Dharwar craton are the residues of after ~35–50% of melt extraction at ~1500–1540 °C and ≤ 2.0 GPa (≤60 km depth). Residues of these high-Mg# (~92–93) peridotites were generated by shallow and a hot plume melting environment. These processes formed a thick depleted mantle residue of primordial cratonic nucleus of Eastern and Western Dharwar. Later, these two blocks collided and were amalgamated, which further thickened the lithospheric mantle beneath the Dharwar craton. The peridotites later were subjected to metasomatic modification, possibly related to proto-kimberlite activity prior to being transported to the surface by the Mesoproterozoic Wajrakarur kimberlites.

Metasomatic control of hydrogen contents in the layered cratonic mantle lithosphere sampled by Lac de Gras xenoliths in the central Slave craton, Canada

Whether hydrogen incorporated in nominally anhydrous mantle minerals plays a role in the strength and longevity of the thick cratonic lithosphere is a matter of debate. In particular, the percolation of hydrogen-bearing melts and fluids could potentially add hydrogen to the mantle lithosphere, weaken its olivines (the dominant mineral in mantle peridotite), and cause delamination of the lithosphere's base. The influence of metasomatism on hydrogen contents of cratonic mantle minerals can be tested in mantle xenoliths from the Slave Craton (Canada) because they show extensive evidence for metasomatism of a layered cratonic mantle. Minerals from mantle xenoliths from the Diavik mine in the Lac de Gras kimberlite area located at the center of the Archean Slave craton were analyzed by FTIR for hydrogen contents. The 18 peridotites, two pyroxenites, one websterite and one wehrlite span an equilibration pressure range from 3.1 to 6.6 GPa and include samples from the shallow (≤145 km), oxidized ultra-depleted layer; the deeper (∼145–180 km), reduced less depleted layer; and an ultra-deep (≥180 km) layer near the base of the lithosphere. Olivine, orthopyroxene, clinopyroxene and garnet from peridotites contain 30–145, 110–225, 105–285, 2–105 ppm H2O, respectively. Within each deep and ultra-deep layer, correlations of hydrogen contents in minerals and tracers of metasomatism (for example light over heavy rare-earth-element ratio (LREE/HREE), high-field-strength-element (HFSE) content with equilibration pressure) can be explained by a chromatographic process occurring during the percolation of kimberlite-like melts through garnet peridotite. The hydrogen content of peridotite minerals is controlled by the compositions of the evolving melt and of the minerals and by mineral/melt partition coefficients. At the beginning of the process, clinopyroxene scavenges most of the hydrogen and garnet most of the HFSE. As the melt evolves and becomes enriched in hydrogen and LREE, olivine and garnet start to incorporate hydrogen and pyroxenes become enriched in LREE. The hydrogen content of peridotite increases with decreasing depth, overall (e.g., from 75 to 138 ppm H2O in the deep peridotites). Effective viscosity calculated using olivine hydrogen content for the deepest xenoliths near the lithosphere-asthenosphere boundary overlaps with estimates of asthenospheric viscosities. These xenoliths cannot be representative of the overall cratonic root because the lack of viscosity contrast would have caused basal erosion of lithosphere. Instead, metasomatism must be confined in narrow zones channeling kimberlite melts through the lithosphere and from where xenoliths are preferentially sampled. Such localized metasomatism by hydrogen-bearing melts therefore does not necessarily result in delamination of the cratonic root.

Eclogites of the North Atlantic Craton: insights from the Chidliak eclogite xenoliths (S. Baffin Island, Canada)

The 156–138 Ma Chidliak kimberlites on the Eastern Hall peninsula (EHP) of Baffin Island entrained mantle xenoliths interpreted to have been a part of the Archean North Atlantic Craton (NAC) lithospheric mantle. We studied 19 Chidliak eclogite xenoliths that comprise 10 bimineralic, 5 rutile-bearing, 3 orthopyroxene-bearing and 1 kyanite-bearing eclogites. We report major and trace element compositions of the minerals, calculated bulk compositions, pressures and temperatures of the rock formation and model melt extraction from viable protoliths. The eclogite samples are classified into three groups of HREE-enriched, LREE-depleted and metasomatized based on their reconstructed whole-rock REE patterns. PT parameters of the eclogites were calculated by projecting garnet–clinopyroxene temperatures onto the local P–T arrays for 65 Chidliak peridotite xenoliths. All Chidliak eclogites are equilibrated in the diamond P–T field and cluster in two groups, low-temperature (n = 5, 840–990 °C at 4.1–5.0 GPa) and high-temperature (n = 11, T > 1320 °C at P > 7.0 GPa). The reconstructed Mg-rich major element bulk compositions and trace elements patterns are similar to Archean basalts from the North Atlantic and Superior cratons and the oceanic gabbros. The LREE-depleted Chidliak eclogites could be residues after 15–55% partial melting of Archean basalt at the eclogite facies of metamorphism that led to extraction of a tonalite–trondhjemite–granodiorite melt from the EHP. The HREE-depleted eclogites may have experienced a lower degree ( 95%. The former are more magnesian, less ferrous and calcic, contain more magnesian and less calcic garnets, and lower proportions of group C eclogites. The contrast may relate to the stronger NAC metasomatism by silicate–carbonate melt observed in Chidliak peridotitic mantle, or to the different formation ages of the eclogites beneath the two cratons.

Archean lithospheric differentiation: Insights from Fe and Zn isotopes

Abstract The Archean continental lithosphere consists of a dominantly felsic continental crust, made of tonalite-trondhjemite-granodiorite (TTG) and subordinate granitoids, and a cratonic lithospheric mantle, made of highly refractory peridotites. Whether they stemmed from the same process of differentiation from the primitive mantle, or were two distinct components that were physically juxtaposed, remains debated. Metal stable isotope ratios are sensitive to magmatic and metamorphic processes and do not evolve with time. Therefore, stable isotope ratios are complementary to radiogenic isotope ratios, and they allow direct comparisons to be made between different terrestrial components without age corrections. Isotopes of iron and zinc, metals ubiquitous in Earth’s lithosphere, can be tracers of lithospheric formation and evolution because they are affected by partial melting (Fe, Zn), redox state (Fe), and the presence of sulfides (Fe, Zn). Here, using stable Fe and Zn isotopic data from Archean samples of the lithospheric mantle and the continental crust, we show that Fe and Zn isotopes define a linear array, best explained by their coupled fractionation behavior during magmatic processes. Our data show that high degrees of partial melting (&amp;gt;30%) during the formation of the cratonic mantle and mafic protocrust, and reworking of the early crust significantly fractionate Fe and Zn isotopes. Conversely, Fe and Zn isotope ratios in the TTG are similar to those in Archean mafic rocks, suggesting an origin by fractional crystallization of basalt, and implying limited Fe and Zn isotopic fractionation, instead of partial melting of mafic crust. Moreover, the absence of Fe and Zn isotope decoupling due to redox effects, melt (fluid)–rock or sediment-rock interaction, and decarbonation indicates that subduction, at least as we understand it now, is not required to explain the Fe and Zn isotope composition of the Archean lithosphere.

The age and origin of cratonic lithospheric mantle: Archean dunites vs. Paleoproterozoic harzburgites from the Udachnaya kimberlite, Siberian craton

Cratonic lithospheric mantle is believed to have been formed in the Archean, but kimberlite-hosted coarse peridotites from Udachnaya in the central Siberian craton typically yield Paleoproterozoic Re-depletion Os isotope ages (TRD). By comparison, olivine megacrysts from Udachnaya, sometimes called “megacrystalline peridotites”, often yield Archean TRD ages, but the nature of these rare materials remains enigmatic. We provide whole-rock (WR) Re-Os isotope and PGE analyses for 24 olivine-rich xenoliths from Udachnaya as well as modal and petrographic data, WR and mineral major and trace element compositions. The samples were selected based on (a) high olivine abundances in hand specimens and (b) sufficient freshness and size to yield representative WR powders. They comprise medium- to coarse-grained (olivine < 1 cm) dunites, a megacrystalline (olivine > 1 cm) dunite, olivine megacrysts and low-orthopyroxene (11–21% opx) harzburgites equilibrated at 783–1154 °C and 3.9–6.5 GPa; coarse dunites have not been previously reported from Udachnaya; two xenoliths contain ilmenite. The harzburgites and dunites have similar WR variation ranges of Ca, Al, Fe, Cr and Mg# (0.917–0.934) typical of refractory cratonic peridotites, but the dunites tend to have higher MgO, NiO and Mg/Si. Mineral abundances and those of Ca and Al are not correlated with Mg#WR; they are not due to differences in melting degrees but are linked to metasomatism. Several samples with high 187Re/188Os show a positive linear correlation with 187Os/188Os with an apparent age of 0.37 Ga, same as eruption age of host kimberlite. Robust TRD ages were obtained for 16 xenoliths with low 187Re/188Os (0.02–0.13). TRD ages for low-opx harzburgites (1.9–2.1 Ga; average 2.0 ± 0.1 Ga, 1 σ) are manifestly lower than for dunites and megacrysts (2.4–3.1 Ga); the latter define two subsets with average TRD of 2.6 ± 0.1 Ga and 3.0 ± 0.1 Ga, and TMA of 3.0 ± 0.2 Ga and 3.3 ± 0.1 Ga, respectively. Differences in olivine grain size (coarse vs. megacrystalline) are not related to age. The age relations suggest that the dunites and megacrysts could not be produced by re-melting of harzburgites, e.g. in arc settings, nor be melt channel materials in harzburgites. Instead, they are relict fragments of lithospheric mantle formed in the Archean (likely in two events at or after 2.6 Ga and 3.0 Ga) that were incorporated into cratonic lithosphere during the final assembly of the Siberian craton in the Paleoproterozoic. A multi-stage formation of the Siberian lithospheric mantle is consistent with crustal basement ages from U-Pb dating of zircons from crustal xenoliths at Udachnaya and detrital zircons from the northern Siberian craton (1.8–2.0, 2.4–2.8 and 3.0–3.4 Ga). The new data from the Siberian and other cratons suggest that the formation of strongly melt-depleted cratonic lithosphere (e.g. Mg# ≥0.92) did not stop at the Archean-Proterozoic boundary as is commonly thought, but continued in the Paleoproterozoic. The same may be valid for the transition from the ‘Archean’ (4–2.5 Ga) to modern tectonic regimes.

Reduced methane-bearing fluids as a source for diamond

Diamond formation in the Earth has been extensively discussed in recent years on the basis of geochemical analysis of natural materials, high-pressure experimental studies, or theoretical aspects. Here, we demonstrate experimentally for the first time, the spontaneous crystallization of diamond from CH4-rich fluids at pressure, temperature and redox conditions approximating those of the deeper parts of the cratonic lithospheric mantle (5–7 GPa) without using diamond seed crystals or carbides. In these experiments the fluid phase is nearly pure methane, even though the oxygen fugacity was significantly above metal saturation. We propose several previously unidentified mechanisms that may promote diamond formation under such conditions and which may also have implications for the origin of sublithospheric diamonds. These include the hydroxylation of silicate minerals like olivine and pyroxene, H2 incorporation into these phases and the “etching” of graphite by H2 and CH4 and reprecipitation as diamond. This study also serves as a demonstration of our new high-pressure experimental technique for obtaining reduced fluids, which is not only relevant for diamond synthesis, but also for investigating the metasomatic origins of diamond in the upper mantle, which has further implications for the deep carbon cycle.

Mesoproterozoic (~1.32 Ga) modification of lithospheric mantle beneath the North China craton caused by break-up of the Columbia supercontinent

The Mesoproterozoic Yanliao diabase sills are widespread at the northern margin of the eastern North China Craton (NCC) and provide a means to examine the evolution of the NCC lithospheric mantle from amalgamation to breakup and dispersal of the Columbia supercontinent. Here, we report new baddeleyite U–Pb ages, whole-rock geochemical and Sr–Nd–Hf isotopic data, and baddeleyite Hf isotopes for diabase rocks from different sills, which provide insights into the source and petrogenesis of the magmas and the Mesoproterozoic evolution of the NCC. In situ baddeleyite U–Pb dating yields consistent ages of ca. 1320 Ma for these sills, suggesting coeval emplacement. The diabase samples are all tholeiitic and have high TiO2 and total Fe2O3, and low alkali contents. They have variable Mg# values and major- and trace-element contents but uniform primitive-mantle-normalized trace-element patterns, with weak enrichment in light rare-earth elements (LREEs) and large-ion lithophile elements (LILEs), and depletion in high-field-strength elements (HFSEs). The whole-rock and in situ isotopic contents of the diabase samples are homogeneous, with initial 87Sr/86Sr ratios of 0.7032–0.7049, εNd(t) values of −0.7 to +3.1, εHf(t) values of +0.8 to +3.7, and baddeleyite εHf(t) values of +1.9 to +7.2. All of these geochemical features suggest that these rocks were derived from high-degree partial melting of a fertile lithospheric mantle source at the depth of spinel stability, with subsequent olivine, clinopyroxene, and plagioclase fractionation. This lithospheric mantle was distinct from the previous Archean to Paleoproterozoic cratonic lithospheric mantle beneath the NCC, indicating a refertilization process by upwelling asthenospheric materials during the Mesoproterozoic. This upwelling, together with other Mesoproterozoic magmatism at the craton margin, would have significantly modified the subcontinental lithospheric mantle of the eastern NCC during the Mesoproterozoic, thereby increasing the susceptibility of the craton to subsequent Mesozoic decratonization.

Sulphur-rich mantle metasomatism of Kaapvaal craton eclogites and its role in redox-controlled platinum group element mobility

Eclogite mantle xenoliths from various kimberlite occurrences on the Kaapvaal craton show evidence for depth- and redox-dependent metasomatic events that led to variable base metal sulphide and incompatible element enrichments. Eclogite xenoliths from the Roberts Victor, Jagersfontein, Kimberley (Kamfersdam) and Premier kimberlites were investigated for their silicate and base metal sulphide geochemistry, stable oxygen isotope compositions and oxybarometry. The variably metasomatised eclogites had basaltic, picritic and gabbroic protolith compositions and have garnet δ18O values that range from +3.3 to +7.9‰, which, when coupled with the trace element characteristics, indicate oceanic lithosphere protoliths that had undergone variable degrees of seawater alteration. The deepest equilibrated eclogites (175–220 km depth) from near the base of the Kaapvaal craton lithosphere are the most refractory and feature significant light rare earth element (LREE) depletions. They show the most oxidised redox compositions with ∆logƒO2 values of FMQ-3.9 to FMQ-1.5. Subtle metasomatic overprinting of these eclogites resulted in base metal sulphide formation with relatively depleted and highly fractionated HSE compositions. These deepest eclogites and their included base metal sulphides suggest interaction with relatively oxidised melts or fluids, which, based on their HSE characteristics, could be related to precursor kimberlite metasomatism that was widespread within the Kaapvaal craton mantle lithosphere. In contrast, eclogites that reside at shallower, “mid-lithospheric” depths (140–180 km) have been enriched in LREE and secondary diopside/phlogopite. Importantly, they host abundant metasomatic base metal sulphides, which have higher HSE contents than those in the deeper eclogites at the lithosphere base. The mid-lithospheric eclogites have more reducing redox compositions (∆logfO2 = FMQ-5.3 − FMQ-3.3) than the eclogites from the lowermost Kaapvaal lithosphere. The compositional overprint of the shallower mantle eclogites resembles basaltic rather than kimberlitic/carbonatitic metasomatism, which is also supported by their relatively reducing redox state. Base metal sulphides from the mid-lithospheric eclogites have HSE abundances and distributions that are similar to Karoo flood basalts from southern Africa, suggesting a link between the identified shallow mantle metasomatism of the Kaapvaal cratonic lithosphere and the Karoo large igneous event during the Mesozoic. The sulphide-hosted platinum group element abundances of the mid-lithospheric eclogites are higher compared with their analogues from the deeper lithospheric eclogites, which in combination with their contrasting oxidation states, may imply redox-controlled HSE mobility during sulphur-rich metasomatism of continental mantle lithosphere.

Trace element mapping of high-pressure, high-temperature experimental samples with laser ablation ICP time-of-flight mass spectrometry – Illuminating melt-rock reactions in the lithospheric mantle

Laser ablation inductively coupled plasma time-of-flight mass spectrometry (LA-ICP-TOFMS) is a fast, relatively high-spatial resolution method which allows to measure and visualize major, minor and trace elements in fine-grained samples. Thus, it is a promising tool for the investigation of high-pressure, high-temperature (HP-HT) reaction experiments that typically produce small grain sizes. Here, we apply LA-ICP-TOFMS mapping to HP-HT reaction experiments by Gervasoni et al. 2017 between a) hydrous eclogite and fertile peridotite, which simulates metasomatism in subduction zones, and b) ultramafic silicate-carbonate melt and peridotite, which simulates metasomatism at the base of cratonic lithospheric mantle. A spatial resolution (i.e., laser spot size) of 5 μm and an acquisition rate of 20 multi-element pixels per second (i.e., laser repetition rate of 20 Hz) were used for mapping of the experimental samples. We demonstrate the applicability of LA-ICP-TOFMS by comparison of our results to maps created with electron probe micro-analysis (EPMA) and elemental concentrations measured by conventional (i.e., crater-drilling) LA-ICP-MS. Our new data suggest that metasomatism in subduction zones is expected to result in the formation of Al–Ti-rich amphiboles that strongly fractionate the rare earth elements (REE) and high-field strength elements (HFSE) from other trace elements. In contrast, metasomatism of the cratonic lithosphere involves fractionation in the light and middle REE, large ion lithophile elements (LILE) and HFSE, leaving Ni and heavy REE behind in olivine and garnet, respectively. The trace element enrichment patterns of the resulting melt bear a general resemblance to those of natural Group I kimberlites.

An experimental study of peridotite dissolution in eclogite-derived melts: Implications for styles of melt-rock interaction in lithospheric mantle beneath the North China Craton

Interaction between eclogite-derived melt and peridotite has played an important role in modifying the cratonic lithospheric mantle and generation of high-Mg igneous rocks. To better understand the effect of peridotite physical state on melt-peridotite interaction, we conducted reaction experiments of lherzolite with two basaltic andesites and a ferro-basalt at temperatures of 1300 °C, 1375 °C, and 1425 °C and a pressure of 2 GPa using the reaction couple method. At 1300 °C, when lherzolite is subsolidus, its dissolution rate is slow, and it is mineralogically and texturally unchanged. Garnet and clinopyroxene precipitate at the melt-rock interface. The low-temperature reaction enriches the melt with SiO2 and Na2O, and depletes the melt with Al2O3, FeO, and CaO. At 1375 °C and 1425 °C, when the lherzolite is partially molten, the dissolution and re-equilibration rates are fast due to grain-scale process that involves dissolution, precipitation, and reprecipitation. Dissolution of olivine and precipitation of orthopyroxene produce a melt-bearing orthopyroxene-rich lithology at the melt-rock interface followed by a melt-bearing harzburgite and a melt-bearing lherzolite. The lithology near the interface can be an orthopyroxene-rich harzburgite or an orthopyroxenite depending on the reacting melt composition. The high-temperature reaction produces melt with increased MgO, FeO, CaO, and Mg#, and decreased SiO2, Al2O3, and Na2O. The dissolution rate obtained from the experiments is used to assess the survival of mantle xenoliths in the early Cretaceous high-Mg diorites. Combined with field observations from mantle samples, the experimental results provide insight into the style of eclogite-derived melt and peridotite interaction in the lithospheric mantle beneath the North China Craton. At the lithosphere-asthenosphere boundary, the reaction is dominated by the high-temperature regime, which produces orthopyroxene-rich lithologies such as orthopyroxene-rich harzburgite and orthopyroxenite-veined peridotite. Above the lithosphere-asthenosphere boundary where peridotite is subsolidus, the reaction is characterized by the low-temperature regime, which produces garnet-bearing or garnet-rich lithologies such as garnet pyroxenite, garnetite, and high-Mg granulite. Interactions between eclogite-derived melts and peridotite in the two regimes are responsible for the geochemical features of the early Cretaceous high-Mg igneous rocks from the North China Craton.

Heterogeneous kimberlite metasomatism revealed from a combined He-Os isotope study of Siberian megacrystalline dunite xenoliths

The Siberian sub-continental lithospheric mantle (SCLM) is regionally heterogeneous due to the complex modification of ancient cratonic material by various metasomatic fluids and interaction with the Siberian plume at ∼250 Ma. Here, we assess the extent and origin of this heterogeneity by analysis of helium isotopes, rhenium-osmium isotopes, trace-element, and highly siderophile element [HSE: Os, Ir, Ru, Pt, Pd, Re] abundances in a suite of mantle-derived megacrystalline dunite xenoliths from the 360 Ma Udachnaya East kimberlite pipe, Siberia. This enables assessment of the style and extent of metasomatism acting to modify the volatile budget of the Siberian SCLM. The olivine megacrysts have 3He/4He values that range from 2.3 to 7.4 RA (where RA = the 3He/4He of air), outside the canonical range for the lithospheric mantle (6.1 ± 2.1 RA; Day et al., 2015). High [He] and low U + Th concentrations in these olivine megacrysts indicate that the Udachnaya megacrystalline dunite xenolith suite has undergone minimal post-eruptive modification of He isotopes by 4He recoil. We instead interpret He isotope variations to reflect pre- or syn-eruptive metasomatic signatures during kimberlite emplacement. The dunites can be divided into two groups, ultimately reflecting modification by two distinct metasomatic styles. Group 1 dunites are characterized by highly unradiogenic Os isotope compositions (γOs360Ma −16 to −13) and ancient melt depletion ages (∼3 Ga) typical of ancient cratonic lithospheric mantle; yet display a range of He isotope ratios (2.3–5.9 RA). This indicates that group 1 dunites have been modified by gas-rich metasomatic fluids acting to change He abundances and isotope ratios, while preserving unradiogenic Os isotope and HSE abundance systematics. Group 2 dunites extend to significantly more radiogenic Os isotope signatures and higher He isotope values (γOs360Ma > −5 to +53; 3He/4He from 3.7 to 7.4 RA). These two groups define a negative hyperbolic array between Os and He isotopes which we interpret to result from metasomatism by variable mixtures between: (1) a gas-rich and HSE-poor fluid (similar to the fluids acting to modify Group 1 dunites) with helium isotope compositions above the MORB range (>8 RA) derived from the asthenosphere, and: (2) strongly radiogenic fluids (∼2 RA; 187Os/188Os > 0.25) likely sourced from within the SCLM. The Udachnaya megacrystalline dunite xenolith suite provides insight into how mantle-plume derived fluids interact with metasomatic fluids within the SCLM.

Fragments of asthenosphere incorporated in the lithospheric mantle underneath the Subei Basin, eastern China: Constraints from geothermobarometric results and water contents of peridotite xenoliths in Cenozoic basalts

Anhydrous, medium/coarse-grained spinel bearing mantle xenoliths from the Subei Basin, Eastern China have mineral arrangements that reflect low energy geometry. Because of clinopyroxene modal contents, they are grouped into cpx-rich lherzolites (cpx ≥ 14percentage), lherzolites (8 < cpx < 14%), cpx-poor lherzolites (6 < cpx ≤ 8%) and harzburgites (cpx up to 4%), without relevant textural differences from the most fertile to the most depleted lithotypes. The cpx-rich lherzolites have mineral chemistry close in composition to Primitive Mantle (PM), whereas cpx-poor lherzolite (and lherzolite) and harzburgite groups cannot be considered as a result of direct melting of the PM source. In addition, the high LREE, Th and U contents, coupled with a Sr enrichment of clinopyroxenes in these lithotypes, indicate the circulation of a silicate melt (with crustal components) in a variably depleted mantle sector well before the entrapment of the xenoliths by the host basalt.Despite the large differences in refractory lithophile element contents (i.e. Ca, Al and REE), the equilibrium temperatures never exceed 1021 °C with a constant difference (<200 °C) from harzburgites to cpx-rich lherzolites. Measured mineral water contents indicate that the whole rock contains, on average, 19 ± 7 ppm of H2O without any systematic variation among rock types nor correlation with Al2O3, light-REE and Yb (or Y) contents of cpx. The cpx H2O contents of cpx-rich lherzolites (41–96 ppm) are, on average, one order of magnitude lower than those theoretically expected (214–530 ppm) for a residuum after a maximum of 3% of PM fractional (≈ bulk) melting in the spinel stability field.The proposed dehydration model suggests that the cold highly refractory harzburgites and cpx-poor lherzolites (and lherzolites?) may represent old cratonic lithospheric mantle modified at depth by the interaction with silicate melts, which may also have involved crustal components. In turn, cpx-rich lherzolites constitute fragments of upwelling fertile asthenosphere, which caused the removal/erosion of the lowermost part of the lithospheric mantle. This asthenosphere portion may have been incorporated in the lithospheric region since the Jurassic and it may have progressively cooled down after one (or more) partial melting episodes. The water depletion can be accounted for a continuous loss by diffusion during the subsolidus chemical-physical readjustment, well after (>5My, based on modelled H2O solid-solid diffusion rate) the occurrence of the last melting episode.

Olivine trace element compositions in diamondiferous lamproites from India: Proxies for magma origins and the nature of the lithospheric mantle beneath the Bastar and Dharwar cratons

The ~1100 Ma CC2 and P13 lamproite dykes in the Wajrakarur Kimberlite Field (WKF), Eastern Dharwar Craton, and ~65 Ma Kodomali and Behradih lamproite diatremes in the Mainpur Kimberlite Field (MKF), Bastar Craton share a similar mineralogy, although the proportions of individual mineral phases vary significantly. The lamproites contain phenocrysts, macrocrysts and microcrysts of olivine set in a groundmass dominated by diopside and phlogopite with a subordinate amount of spinel, perovskite, apatite and serpentine along with rare barite. K-richterite occurs as inclusion in olivine phenocrysts in Kodomali, while it is a late groundmass phase in Behradih and CC2. Mineralogically, the studied intrusions are classified as olivine lamproites. Based on microtextures and compositions, three distinct populations of olivine are recognised. The first population comprises Mg-rich olivine macrocrysts (Fo89–93), which are interpreted to be xenocrysts derived from disaggregated mantle peridotites. The second population includes Fe-rich olivine macrocrysts (Fo82–89), which are suggested to be the product of metasomatism of mantle wall-rock by precursor lamproite melts. The third population comprises phenocrysts and overgrowth rims (Fo83–92), which are clearly of magmatic origin. The Mn and Al systematics of Mg-rich olivine xenocrysts indicate an origin from diverse mantle lithologies including garnet peridotite, garnet–spinel peridotite and spinel peridotite beneath the WKF, and mostly from garnet peridotite beneath the MKF. Modelling of temperatures calculated using the Al-in-olivine thermometer for olivine xenocrysts indicates a hotter palaeogeotherm of the SCLM beneath the WKF (between 41 and 43 mW/m2) at ~1100 Ma than beneath the MKF (between 38 and 41 mW/m2) at ~65 Ma. Further, a higher degree of metasomatism of the SCLM by precursor lamproite melts has occurred beneath the WKF compared to the MKF based on the extent of CaTi enrichment in Fe-rich olivine macrocrysts. For different lamproite intrusions within a given volcanic field, lower Fo olivine overgrowth rims are correlated with higher phlogopite plus oxide mineral abundances. A comparison of olivine overgrowth rims from the two fields shows that WKF olivines with lower Fo content than MKF olivines are associated with increased XMg in spinel and phlogopite and vice versa. Melt modelling indicates relatively Fe-rich parental melt for WKF intrusions compared to MKF intrusions. The Ni/Mg and Mn/Fe systematics of magmatic olivines indicate derivation of the lamproite melts from mantle source rocks with a higher proportion of phlogopite and/or lower proportion of orthopyroxene for the WKF on the Eastern Dharwar Craton compared to those for the MKF on the Bastar Craton. This study highlights how olivine cores provide important insights into the composition and thermal state of cratonic mantle lithosphere as sampled by lamproites, including clues to elusive precursor metasomatic events. Variable compositions of olivine rims testify to the complex interplay of parental magma composition and localised crystallisation conditions including oxygen fugacity variations, co-crystallisation of groundmass minerals, and assimilation of entrained material.

Craton Destruction 2: Evolution of Cratonic Lithosphere After a Rapid Keel Delamination Event

AbstractCratonic lithosphere beneath the eastern North China Craton has undergone extensive destruction since early Jurassic times (approximately 190 Ma). This is recorded in its episodic tectonic and magmatic history. In this time, its lithosphere changed thickness from approximately 200 km to &lt;60 km. This change was associated with a peak time (approximately 120 Ma) of lithospheric thinning and magmatism that was linked with high surface heat flow recorded in rift basins. We believe that these records are best explained by a two‐stage evolutionary process. First, approximately 100 km of cratonic “keel” underlying a weak midlithospheric discontinuity layer (approximately 80–100 km) was rapidly removed in &lt;10–20 Ma. This keel delamination stage was followed by a protracted (approximately 50–100 Ma) period of convective erosion and/or lithospheric extension that thinned the remaining lithosphere and continuously reworked the former cratonic lithospheric mantle. This study focuses on numerical exploration of the well‐recorded second stage of the eastern North China Craton's lithospheric evolution. We find that (1) lithospheric mantle capped by thick crust can be locally replaced by deeper mantle material in 100 Ma due to small‐scale convective erosion; (2) asthenospheric upwelling and related extension can replace lithospheric mantle over horizontal length scales of ~50–150 km, and account for observed “mushroom‐shaped” low‐velocity structures; (3) modeling shows conditions that could lead to the multiple eastern North China Craton magmatic pulses between 190 and 115 Ma that are associated with temporal and spatial changes in magma source petrology and a magmatic hiatus; and (4) a “wet” midlithospheric discontinuity layer provides a potential source material for on‐craton magmatism.

Making Cratonic Lithospheric Mantle

AbstractThe origin of cratonic lithospheric mantle has been attributed to either high‐pressure (5–7 GPa) melting in hot mantle plumes or low‐pressure (&lt;5 GPa) melting in mid‐ocean ridges or suprasubduction zones. To resolve this long‐standing debate, it is necessary to confirm under what depths the incipient cratonic mantle melted. Compared with most cratonic mantle xenoliths and xenocrysts, which commonly experienced metasomatic modification after melt extraction, diamond inclusions with predominantly harzburgitic compositions can better track the compositional signature of pristine cratonic mantle. This paper presents thermodynamic calculations performed with THERMOCALC software to establish a quantitative relationship between garnet and cratonic peridotite compositions. Along the normal cratonic geotherm (40 mW/m2), the XCa [atomic Ca/(Ca + Mg + Fe2+)] and Cr# [atomic Cr/(Cr + Al) × 100] values in garnet depend mainly on the bulk CaO/Al2O3 and Cr# values, respectively. Furthermore, mantle melting modeling shows that the Cr# of the residue displays a negative correlation with melting pressure at melt fractions of greater than ~15%. Therefore, the high Cr# values (mostly 12–36) of global garnet inclusions in diamond support a low‐pressure (&lt;4 GPa) origin for cratonic lithospheric mantle. More importantly, the incipient cratonic mantle calculated from the garnet inclusions exhibits similar bulk CaO/Al2O3 and Cr# compositions to oceanic lithosphere but different values from arc lithosphere with higher CaO/Al2O3 and Cr# values. These results imply that cratonic mantle was initially formed by the extensive melting of hot ambient mantle in the Archean ocean ridge environments. The shallow oceanic lithosphere was subsequently stacked to generate the thick and stable cratonic lithospheric mantle.

Diamondiferous Paleoproterozoic mantle roots beneath Arctic Canada: A study of mantle xenoliths from Parry Peninsula and Central Victoria Island

While the mantle roots directly beneath Archean cratons have been relatively well studied because of their economic importance, much less is known about the genesis, age, composition and thickness of the mantle lithosphere beneath the regions that surround the cratons. Despite this knowledge gap, it is fundamentally important to establish the nature of relationships between this circum-cratonic mantle and that beneath the cratons, including the diamond potential of circum-cratonic regions. Here we present mineral and bulk elemental and isotopic compositions for kimberlite-borne mantle xenoliths from the Parry Peninsula and Central Victoria Island, Arctic Canada. These xenoliths provide key windows into the lithospheric mantle underpinning regions to the North and Northwest of the Archean Slave craton, where the presence of cratonic material has been proposed. The mantle xenolith data are supplemented by mineral concentrate data obtained during diamond exploration. The mineral and whole rock chemistry of peridotites from both localities is indistinguishable from that of typical cratonic mantle lithosphere. The cool mantle paleogeotherms defined by mineral thermobarometry reveal that the lithospheric mantle beneath the Parry Peninsula and Central Victoria Island terranes extended well into the diamond stability field at the time of kimberlite eruption, and this is consistent with the recovery of diamonds from both kimberlite fields. Bulk xenolith Se and Te contents, and highly siderophile element (including Os, Ir, Pt, Pd and Re) abundance systematics, plus corresponding depletion ages derived from Re-Os isotope data suggest that the mantle beneath these parts of Arctic Canada formed in the Paleoproterozoic Era, at ∼2 Ga, rather than in the Archean. The presence of a diamondiferous Paleoproterozoic mantle root is part of the growing body of global evidence for diamond generation in mantle roots that stabilized well after the Archean. In the context of regional tectonics, we interpret the highly depleted mantle compositions beneath both studied regions as formed by mantle melting associated with hydrous metasomatism in the major Paleoproterozoic Wopmay-Great Bear-Hottah arc systems. These ∼2 Ga arc systems were subsequently accreted along the margin of the Slave craton to form a craton-like thick lithosphere with diamond potential thereby demonstrating the importance of subduction accretion in building up Earth’s long-lived continental terranes.

Multidisciplinary Constraints on the Abundance of Diamond and Eclogite in the Cratonic Lithosphere

AbstractSome seismic models derived from tomographic studies indicate elevated shear‐wave velocities (≥4.7 km/s) around 120–150 km depth in cratonic lithospheric mantle. These velocities are higher than those of cratonic peridotites, even assuming a cold cratonic geotherm (i.e., 35 mW/m2 surface heat flux) and accounting for compositional heterogeneity in cratonic peridotite xenoliths and the effects of anelasticity. We reviewed various geophysical and petrologic constraints on the nature of cratonic roots (seismic velocities, lithology/mineralogy, electrical conductivity, and gravity) and explored a range of permissible rock and mineral assemblages that can explain the high seismic velocities. These constraints suggest that diamond and eclogite are the most likely high‐Vs candidates to explain the observed velocities, but matching the high shear‐wave velocities requires either a large proportion of eclogite (&gt;50 vol.%) or the presence of up to 3 vol.% diamond, with the exact values depending on peridotite and eclogite compositions and the geotherm. Both of these estimates are higher than predicted by observations made on natural samples from kimberlites. However, a combination of ≤20 vol.% eclogite and ~2 vol.% diamond may account for high shear‐wave velocities, in proportions consistent with multiple geophysical observables, data from natural samples, and within mass balance constraints for global carbon. Our results further show that cratonic thermal structure need not be significantly cooler than determined from xenolith thermobarometry.

Seismic evidence for depth-dependent metasomatism in cratons

The long-term stability of cratons has been attributed to low temperatures and depletion in iron and water, which decrease density and increase viscosity. However, steady-state thermal models based on heat flow and xenolith constraints systematically overpredict the seismic velocity-depth gradients in cratonic lithospheric mantle. Here we invert for the 1-D thermal structure and a depth distribution of metasomatic minerals that fit average Rayleigh-wave dispersion curves for the Archean Kaapvaal, Yilgarn and Slave cratons and the Proterozoic Baltic Shield below Finland. To match the seismic profiles, we need a significant amount of hydrous and/or carbonate minerals in the shallow lithospheric mantle, starting between the Moho and 70 km depth and extending down to at least 100–150 km. The metasomatic component can consist of 0.5–1 wt% water bound in amphibole, antigorite and chlorite, ∼0.2 wt% water plus potassium to form phlogopite, or ∼5 wt% CO2 plus Ca for carbonate, or a combination of these. Lithospheric temperatures that fit the seismic data are consistent with heat flow constraints, but most are lower than those inferred from xenolith geothermobarometry. The dispersion data require differences in Moho heat flux between individual cratons, and sublithospheric mantle temperatures that are 100–200 °C less beneath Yilgarn, Slave and Finland than beneath Kaapvaal. Significant upward-increasing metasomatism by water and CO2-rich fluids is not only a plausible mechanism to explain the average seismic structure of cratonic lithosphere but such metasomatism may also lead to the formation of mid-lithospheric discontinuities and would contribute to the positive chemical buoyancy of cratonic roots.