Density structure of the cratonic mantle in Southern Africa: 2. Correlations with kimberlite distribution,seismic velocities,and Moho sharpness

We present a new regional model for the depth-averaged density structure of the cratonic lithospheric mantle in southern Africa constrained on a 30′ × 30′ grid and discuss it in relation to regional seismic models for the crust and upper mantle, geochemical data on kimberlite-hosted mantle xenoliths, and data on kimberlite ages and distribution. Our calculations of mantle density are based on free-board constraints, account for mantle contribution to surface topography of ca. 0.5–1.0 km, and have uncertainty ranging from ca. 0.01 g/cm³ for the Archean terrains to ca. 0.03 g/cm³ for the adjacent fold belts. We demonstrate that in southern Africa, the lithospheric mantle has a general trend in mantle density increase from Archean to younger lithospheric terranes. Density of the Kaapvaal mantle is typically cratonic, with a subtle difference between the eastern, more depleted, (3.31–3.33 g/cm³) and the western (3.32–3.34 g/cm³) blocks. The Witwatersrand basin and the Bushveld Intrusion Complex appear as distinct blocks with an increased mantle density (3.34–3.35 g/cm³) with values typical of Proterozoic rather than Archean mantle. We attribute a significantly increased mantle density in these tectonic units and beneath the Archean Limpopo belt (3.34–3.37 g/cm³) to melt-metasomatism with an addition of a basaltic component. The Proterozoic Kheis, Okwa, and Namaqua–Natal belts and the Western Cape Fold Belt with the late Proterozoic basement have an overall fertile mantle (ca. 3.37 g/cm³) with local (100–300 km across) low-density (down to 3.34 g/cm³) and high-density (up to 3.41 g/cm³) anomalies. High (3.40–3.42 g/cm³) mantle densities beneath the Eastern Cape Fold belt require the presence of a significant amount of eclogite in the mantle, such as associated with subducted oceanic slabs.We find a strong correlation between the calculated density of the lithospheric mantle, the crustal structure, the spatial pattern of kimberlites, and their emplacement ages. (1) Blocks with the lowest values of mantle density (ca. 3.30 g/cm³) are not sampled by kimberlites and may represent the “pristine” Archean mantle. (2) Young (< 90 Ma) Group I kimberlites sample mantle with higher density (3.35 ± 0.03 g/cm³) than the older Group II kimberlites (3.33 ± 0.01 g/cm³), but the results may be biased by incomplete information on kimberlite ages. (3) Diamondiferous kimberlites are characteristic of regions with a low-density cratonic mantle (3.32–3.35 g/cm³), while non-diamondiferous kimberlites sample mantle with a broad range of density values. (4) Kimberlite-rich regions have a strong seismic velocity contrast at the Moho, thin crust (35–40 km) and low-density (3.32–3.33 g/cm³) mantle, while kimberlite-poor regions have a transitional Moho, thick crust (40–50 km), and denser mantle (3.34–3.36 g/cm³). We explain this pattern by a lithosphere-scale (presumably, pre-kimberlite) magmatic event in kimberlite-poor regions, which affected the Moho sharpness and the crustal thickness through magmatic underplating and modified the composition and rheology of the lithospheric mantle to make it unfavorable for consequent kimberlite eruptions. (5) Density anomalies in the lithospheric mantle show inverse correlation with seismic Vp, Vs velocities at 100–150 km depth. However, this correlation is weaker than reported in experimental studies and indicates that density-velocity relationship in the cratonic mantle is strongly non-unique.     

The structure of mantle flows and stress fields in a two-dimensional convection model with non-Newtonian viscosity

The structure of mantle convection and spatial fields of superlitho static pressure and vertical and horizontal stresses in the Earth’s mantle are studied in a 2D numerical model with non-Newtonian viscosity and heat sources. The model demonstrates a jump-like motion of subduction zones and reveals abrupt changes in the stress fields depending on the stage of slab detachment. The stresses decrease dramatically in the areas without slabs. The horizontal stresses o_xx, superlitho static pressure, and vertical stresses o_zz in the part of the mantle lacking intense near-vertical flows are approximately equal, varying within ± 6, ± 8, and ± 10 MPa, respectively. However, these fields are stronger in the areas of descending slabs, where the values of the above parameters are about an order of magnitude higher (± 50 MPa).This result agrees with the current views of the oceanic slabs as the most important gent of mantle convection. We have found significant differences among the o_xx, o_zz, and pressure fields. The pressure field reveals both the vertical and horizontal features of slabs and plumes, clearly showing their long thermal conduits with broader heads. The distributions of o_xx are sensitive to the near-horizontal parts of the flows, whereas the o_zz fields reveal mainly their vertical substructures. The model shows the presence of relatively cold remnants of slabs in the lower mantle above the thermal boundary layer. Numerous hot plumes penetrating through these high-viscosity remnants, as well as the new descending slabs, induce intense stress fields in the lower mantle, which are strongly inhomogeneous in space and time.     

Lithospheric mantle xenoliths sampled by melts from upwelling asthenosphere: The Quaternary Tasse alkaline basalts of southeastern British Columbia,Canada

The Quaternary Tasse basalts are exposed near the north shore of Quesnel Lake in southeastern British Columbia. They host a variety of mantle xenoliths consisting predominantly of spinel lherzolite with minor dunite and pyroxenite. Mineralogically, the xenoliths are composed of olivine, orthopyroxene, clinopyroxene and spinel characterized by forsterite (Fo_87–93), enstatite (En_90–92), diopside (En_45–50–Wo_40–45–Fs₅), and Cr-spinel (6 − 11 wt.% Cr), respectively. All of the mantle xenoliths are coarse-grained and show granoblastic textures. Clinopyroxene and spinel display textural evidence for chemical reactions with percolating melts.The mantle xenoliths are characterized by restricted Mg-numbers (89 − 92) and low abundances of incompatible elements (Ba = 2 − 11 ppm; Sr = 3 − 31 ppm) and Yttrium (1 − 3 ppm). On the basis of REE patterns, the xenoliths are divided into three groups reflecting the various degrees of mantle metasomatism: (1) Group 1 consists of concave-up LREE patterns (La/Sm_cn = 0.48 − 1.16; Gd/Yb_cn = 0.71 − 0.92); (2) Group 2 possesses flat to moderately LREE-enriched patterns (La/Sm_cn = 1.14 − 1.92; Gd/Yb_cn = 0.87 − 1.09); and (3) Group 3 is characterized by strongly LREE-enriched patterns (La/Sm_cn = 1.53 − 2.45; Gd/Yb_cn = 1.00 − 1.32). On MORB-normalized trace element diagrams, the majority of the xenolith samples share the enrichment of LILE (Rb, Ba, K), U, Th, Pb, Sr and the depletion of HFSE (Nb, Ta, Ti, Y) relative to REE. These geochemical characteristics are consistent with a compositionally heterogeneous subcontinental lithospheric mantle source that originated as subarc mantle wedge peridotite at a convergent plate margin.The Tasse basalts have alkaline compositions characterized by low SiO₂ (44 − 46 wt.%) and high alkali (Na₂O + K₂O = 5.1 − 6.6 wt.%) contents. They are strongly enriched in incompatible elements (TiO₂ = 2.4 − 3.1 wt.%; Ba = 580 − 797 ppm; Sr = 872 − 993 ppm) and, display OIB-like trace element patterns (La/Sm_n = 3.15 − 3.85; Gd/Yb_n = 3.42 − 4.61). They have positive εNd (+ 3.8 to + 5.5) values, with 338 − 426 Ma depleted mantle model ages, and display uniform OIB-like Sr (⁸⁷Sr/⁸⁶Sr = 0.703346 − 0.703591) and Pb (²⁰⁶Pb/²⁰⁴Pb = 19.40 − 19.58; ²⁰⁷Pb/²⁰⁴Pb = 15.57 − 15.60; ²⁰⁸Pb/²⁰⁴Pb = 38.99 − 39.14) isotopic compositions. The basalts erupted discontinuously along a > 1000 km long SE-NW-trending linear belt with minimal compositional variation indicative of a homogenous mantle source. The Sr − Nd − Pb isotope and trace element systematics of the alkaline basalts suggests that they originated from partial melting of an upwelling asthenospheric mantle source. Melting of the asthenospheric mantle might have stemmed from extension of the overlying lithosphere in response to the early stages of back-arc basin opening in the Omineca and Intermontane belts. Ridge subduction beneath the Canadian Cordillera might have played an important role in the weakening of the lithospheric mantle prior to its extension. Alternatively, melting of the upwelling asthenosphere in response to the delamination of the lithospheric mantle beneath the Rocky Mountain Trench might have generated the alkaline lavas.     

Lithospheric petrology of the eastern Arabian Plate: Constraints from Al-Ashkhara (Oman) xenoliths

Mafic granulite and spinel lherzolite xenoliths from Cenozoic alkaline basalts near Al-Ashkhara, eastern Oman, have been selected for a systematic mineralogical, geochemical and Sr–Nd–Pb isotopic study. This is the only place in E Arabia where samples of both lower crust and upper mantle can be examined. Lower crustal xenoliths consist of two mineralogically and chemically distinct groups: gabbronorite (subequal abundances of ortho- and clino-pyroxene and plagioclase) and plagioclase pyroxenite (dominant pyroxene and subordinate plagioclase). Temperature estimates for lower crustal xenoliths using the two pyroxene geothermometer (T-Wells) yield 810–865 °C. The mineral assemblage (spinel–pyroxene–plagioclase) and Al content in pyroxene indicate that plagioclase-bearing xenoliths equilibrated at 5–8 kbar (13 and 30 km depth) in the lower crust. εNd and ⁸⁷Sr/⁸⁶Sr calculated at 700 Ma for Al-Ashkhara lower crustal xenoliths (+ 6.4 to + 6.6; ⁸⁷Sr/⁸⁶Sr = 0.7028 to 0.7039) are consistent with the interpretation that juvenile, mafic melts were added to the lower crust during Neoproterozoic time and that there was no discernible contribution from pre-Neoproterozoic crust. Upper mantle xenoliths consist of both dry and hydrous (phlogopite-bearing) lherzolites. These peridotites are more Fe-rich than expected for primitive mantle or melt residues and probably formed by pervasive circulation of melts that have refertilized pre-existing mantle peridotites. Mineral equilibration temperatures range from 990 to 1070 °C. Isotopic compositions calculated at 700 Ma are εNd = + 6.8 to + 7.8 and ⁸⁷Sr/⁸⁶Sr = 0.7016 to 0.7025, indicating depleted upper mantle. Pb isotopic compositions indicate that the metasomatism was relatively recent, perhaps related to Paleogene tectonics and basanite igneous activity. Nd model ages for the spinel peridotite xenoliths range between 0.59 and 0.65 Ga. The xenolith data suggest that eastern Arabian lower crust is of hotspot origin, in contrast to western Arabian lower crust, which mostly formed at a convergent plate margin. Geochemical and isotopic differences between lower crust and upper mantle indicate that these are unrelated, possibly because delamination replaced the E Arabian mantle root in Neoproterozoic time.     

Ultramafic xenoliths from Lake Nyos area,Cameroon volcanic line,West-central Africa: Petrography,mineral chemistry,equilibration conditions and metasomatic features

Spinel-bearing mantle xenoliths have been recovered in the pyroclastic breccia surrounding the Lake Nyos maar. These include spinel lherzolites, spinel harzburgites and olivine websterites. They exhibit coarse granular or protogranular to weakly porphyroclastic textures, and show variations in mineral chemistry, modal compositions and equilibrium temperature. The xenoliths consist of four mineral phases typical of upper mantle origin: olivine (Fo₈₉–Fo_91.5, NiO = 0.29–0.38 wt%, CaO = 0.02–0.17 wt%), enstatite (Mg# = 90–92, Cr₂O₃ = 0.35 ± 0.04 wt%), Cr-diopside (Mg# = 92–98, Cr₂O₃ = 0.7–1.65 wt%, TiO₂ = 0.26–0.6 wt%) and spinel (high Mg# of 70–80, low TiO₂ ≤ 0.4 wt%). Spinels are aluminous (Cr# = 9.7–11) in most lherzolites, and become increasingly chromiferous from websterites (Cr#_Sp = 15.3–19.8) to harzburgites (Cr#_Sp = 19–33.6). The lherzolites are composed of olivine (48–58%), orthopyroxene (22–30%), and clinopyroxene (8–15%). The harzburgites modes are olivine (60–81%), orthopyroxene (11–29%), and clinopyroxene (<5%). The websterites are mainly composed of pyroxene (～62%) with variable amounts of olivine (23–31%). Temperatures of mineral equilibration in the xenoliths have been estimated from the two-pyroxene thermometer of Wells (1977) and range between 850 and 1050 °C, corresponding to about 10–30 kbar at a depth mantle of 30 km at least. These P–T conditions show significant variations between different petrographical types, the maximum conditions being recorded in two spinel lherzolites (NY-05 and NY-23) that have atypical chemical compositions and textures suggesting that they were initially formed in an environment close to the garnet stability field, then re-equilibrated within the spinel stability field prior to their incorporation in the host magma. With the exception of minerals from these two lherzolite nodules, all the minerals exhibit depletion of light REE, a typical feature of abyssal peridotites implying that some xenoliths from the Cameroon volcanic line were probably sampled in a part of the sub-continental mantle that is chemically similar to sub-oceanic mantle. The variations observed in the mineral chemistry and modal compositions of xenoliths suggest that the spinel harzbugite nodules which represent residues of a significant degree of partial melting of lherzolitic mantle were affected by infiltration of alkali-enriched metasomatizing melts (or fluids) within the uppermost mantle to produce pargasitic amphiboles prior to their sampling by the host lava. The features of this metasomatism event occur in the rocks of all three petrographical facies xenoliths from Lake Nyos.     

Isotope (U–Pb,Sm–Nd,Rb–Sr) geochronology of alkaline basic plutons of the Kuznetsk Alatau

On the northeastern slope of the Kuznetsk Alatau, small differentiated alkaline basic intrusive massifs form an isometric area ~ 100 km across. They are composed of subalkalic and alkali gabbroids, basic and ultrabasic foidolites, nepheline and alkali syenites, and carbonatites. Results of complex (U–Pb, Sm–Nd, and Rb–Sr) isotope dating suggest that alkaline basic magmatism developed at two stages, in the Middle Cambrian–Early Ordovician (~ 510–480 Ma) and in the Early–Middle Devonian (~ 410–385 Ma). Finding of accessory zircons (age 1.3–2.0 Ga) in alkaline rocks suggests that the ascent of mantle plume was accompanied by the melting of fragments of Proterozoic mature continental crust composing the basement of the Caledonian orogen of the Kuznetsk Alatau. Probably, parental Cambrian–Ordovician alkaline mafic melts initiated metasomatism and lithosphere erosion. During the next melting of lithosphere substrate in ~ 100 Myr, this caused the generation of magmas of similar composition with inherited isotope parameters (ε_Nd(T) ≈ + 4.8 to + 5.7, T_Nd(DM) ≈ 0.8–0.9 Ga) pointing to the similar nature of their matter sources in the moderately depleted mantle.     

Petrogenesis of high-Ti mafic dykes from Southern Qiangtang,Tibet: Implications for a ca. 290 Ma large igneous province related to the early Permian rifting of Gondwana

The presence and/or generation mechanism of a mantle plume associated with early Permian rifting on the northern margin of Gondwana are topics of debate. Here we report LA–ICP–MS U–Pb zircon ages, whole-rock geochemistry, and Sr–Nd isotope data for high-Ti mafic dykes from southern Qiangtang, Tibet, with the aim of assessing if a mantle plume formed in this region during the early Permian. Zircon U–Pb dating of diabase dykes yielded ages of 290.6 ± 3.5 Ma and 290.1 ± 1.5 Ma, indicating they were emplaced during the early Permian. Whole-rock geochemistry shows that these mafic dykes are alkaline (Nb/Y = 0.73–0.99), have high TiO2 (3.6%–4.8%), and have ocean-island basalt (OIB)-like trace element patterns with enrichments in Nb, Ta, and Ti. Whole-rock Sr–Nd isotope data show a relatively narrow range of εNd(t) (+ 2.29 to + 3.53), similar to basalts produced by a mantle plume (e.g., Emeishan continental flood basalts (ECFB)). Elemental and isotope data suggest that the dykes have undergone fractionation crystallization of mafic minerals and have experienced negligible crustal contamination. These mafic rocks show an affinity to OIB and may have been generated by partial melting of an OIB-type, garnet-bearing asthenospheric mantle source. On the basis of a similar emplacement age to the Panjal Traps basalts in the Himalayas, combined with a tectonic reconstruction of Gondwana in the early to middle Permian, our work suggests that the high-Ti mafic dykes in the Southern Qiangtang terrane and the coeval Panjal Traps basalts in the Himalayas together comprise a ca. 290 Ma large igneous province linked to a mantle plume, which probably played an active role in early Permian rifting on the northern margin of Gondwana and was related to circum-Pangea subduction.     

Lithospheric mantle structure of the Siberian craton inferred from the superlong Meteorite and Rift seismic profiles

Modeling of the seismic, thermal, and density structure of the Siberian craton lithospheric mantle at depths of 100-300 km has been performed along the superlong Meteorite and Rift seismic profiles. The 2D velocity sections reflect the specific features of the internal structure of the craton: lateral inhomogeneities, seismic-boundary relief at depths of ~ 100, 150, 240, and 300 km, velocities of 8.3-8.7 km/s, and the lack of low-velocity zone in the lower lithosphere. Mapping of the thermal state along the Meteorite and Rift profiles shows a significant temperature decrease in the cratonic mantle as compared with the average temperatures of the surrounding Phanerozoic mantle (> 300 °C) estimated from the global reference model AK135. Lateral temperature variations, reflecting the thermal anomalies in the cratonic keel, are observed at depths of < 200 km (with some decrease in temperature in the central part of the craton), whereas at depths of > 200 km, temperature variations are negligible. This suggests the preservation of residual thermal perturbations at the base of the lithosphere, which must lead to the temperature equalization in the transition zone between the lithosphere and the asthenosphere. Variations in chemical composition have a negligible effect on the thermal state but affect strongly the density structure of the mantle. The results of modeling admit a significant fertilization of matter at depths more than 180-200 km and stratification of the cratonic mantle by chemical composition. The thicknesses of chemical (petrologic) and thermal boundary layers beneath the Siberian craton are estimated. The petrologic lithosphere is localized at depths of ~ 200 km. The bottom of the thermal boundary layer is close to the 1450 °C isotherm and is localized at a depth of 300 km, which agrees with heat flow and seismic-tomography data.     

Rhodium,gold and other highly siderophile elements in orogenic peridotites and peridotite xenoliths

The concentrations of Rh, Au and other highly siderophile elements (HSE: Re, Os, Ir, Ru, Pt, Rh, Pd and Au), and ¹⁸⁷Os/¹⁸⁸Os isotope ratios have been determined for samples from peridotite massifs and xenoliths in order to further constrain HSE abundances in the Earth's mantle and to place constraints on the distributions processes accounting for observed HSE variations between fertile and depleted mantle lithologies. Concentrations of Re, Os, Ir, Ru, Pt and Pd were determined by isotope dilution ICP-MS and N-TIMS. The monoisotopic elements Rh and Au were quantified by standardization relative to the concentrations of Ru and Ir, respectively, and were determined from the same digestion aliquot as other HSE. The measurement precision of the concentration data under intermediate precision conditions, as inferred from repeated analyses of 2 g test portions of powdered samples, is estimated to be better than 10% for Rh and better than 15% for Au (1 s).Fertile lherzolites display non-systematic variation of Rh concentrations and constant Rh/Ir of 0.34 ± 0.03 (1 s, n = 57), indicating a Rh abundance for the primitive mantle of 1.2 ± 0.2 ng/g. The data also suggest that Rh behaves as a compatible element during low to moderate degrees of partial melting in the mantle or melt–mantle interaction, but may be depleted at higher degrees of melting. In contrast, Au concentrations and Au/Ir correlate with peridotite fertility, indicating incompatible behaviour of Au during magmatic processes in the mantle. Fertile lherzolites display Au/Ir ranging from 0.20 to 0.65, whereas residual harzburgites have Au/Ir < 0.20. Concentrations of Au and Re are correlated with each other and suggest similar compatibility of both elements. The primitive mantle abundance of Au calculated from correlations displayed by Au/Ir with Al₂O₃ and Au with Re is 1.7 ± 0.5 ng/g (1 s).The depletion of Pt, Pd, Re and Au relative to Os, Ir, Ru and Rh displayed by residual harzburgites, suggests HSE fractionation during partial melting. However, the HSE abundance variations of fertile and depleted peridotites cannot be explained by a simple fractionation process. Correlations displayed by Pd/Ir, Re/Ir and Au/Ir with Al₂O₃ may reflect refertilization of previously melt depleted mantle rocks due to reactive infiltration of silicate melts.Relative concentrations of Rh and Au inferred for the primitive mantle model composition are similar to values of ordinary and enstatite chondrites, but distinct from carbonaceous chondrites. The HSE pattern of the primitive mantle is inconsistent with compositions of known chondrite groups. The primitive mantle composition may be explained by late accretion of a mixture of chondritic with slightly suprachondritic materials, or alternatively, by meteoritic materials mixed into mantle with a HSE signature inherited from core formation.     

Chemical variations of abyssal peridotites in the central Oman ophiolite: Evidence of oceanic mantle heterogeneity

Basal peridotites above the metamorphic sole outcropped around Wadi Sarami in the central Oman ophiolite give us an excellent opportunity to understand the spatial extent of the mantle heterogeneity and to examine peridotites−slab interactions. We recognized two types of basal lherzolites (Types I and II) that change upward to harzburgites. Their pyroxene and spinel compositions display severely variations at small scales over < 0.5 km, and encompass the entire abyssal peridotite trend; clinopyroxenes (Cpxs) show wide ranges of Al₂O₃, Na₂O, Cr₂O₃ and TiO₂ contents. Primary spinels show a large variation of Cr# [= Cr/(Cr + Al)] from 0.04 to 0.53, indicating various degrees of partial melting. Trace-element compositions of peridotites and their pyroxenes also show a large chemical heterogeneity in the base of the Oman mantle section. This heterogeneity mainly resulted from variations of partial-melting degrees due to the change of a mantle thermal regime and a distance from the spreading ridge or the mantle diapir. It was overlapped with subsolidus modification during cooling and fluid metasomatism prior and/or during emplacement. The studied peridotites are enriched in Rb, Cs, Ba, Sr and LREE due to fluid influx during detachment and emplacement stages. Chondrite (CI)-normalized REE patterns for pyroxenes are convex upward with strong LREE depletion due to their residual origin, similar to abyssal peridotites from a normal ridge segment. The Cpxs are enriched in fluid mobile elements (e.g., B, Li, Cs, Pb, Rb) and depleted in HFSE (Ta, Nb, Th, Zr) + LREE, suggesting no effect of melt refertilization. Their HREE contents, combined with spinel compositions, suggest two melting series with 1–5% melting for type II lherzolites, 3– < 10% melting for type I lherzolites and ~ 15% for harzburgites. Hornblendes are enriched in fluid-mobile elements relative to HFSE + U inherited from their precursor Cpx. The clinopyroxenite lens crosscuts the basal lherzolites, forming small-scale (< 5 cm) mineralogical and chemical heterogeneities. It was possibly formed from fractional crystallization of interstitial incremental melt that formed during decompression melting of a normal MORB mantle source. The studied peridotites possibly represent a chemical heterogeneity common to the mantle at an oceanic spreading center.     

Compositional signatures of SSZ-type peridotites from the northern Vourinos ultra-depleted upper mantle suite,NW Greece

The northern Vourinos massif, located in the Dinarides-Hellenides mountain belt in the Balkan Peninsula, forms a section of the so-called Neotethyan ophiolitic belt in the Alpine-Himalayan orogenic system. It is comprised mainly of a well-preserved mantle sequence, dominated by voluminous massive harzburgite with variable clinopyroxene and olivine modal abundances, accompanied by subordinate coarse- and fine-grained dunite. The harzburgite rock varieties are characterized by high Cr# [Cr/(Cr + Al)] values in Cr-spinel (0.47–0.74), elevated Mg# [Mg/(Mg + Fe²⁺)] in olivine (0.90–0.93), low Al₂O₃ content in clinopyroxene (≤1.82 wt.%) and low average bulk-rock concentrations of CaO (0.52 wt.%) and Al₂O₃ (0.40 wt.%), which are indicative of their refractory nature. In addition, dunite-type rocks display even more depleted compositions, containing Cr-spinel and olivine with higher Cr# (0.76–0.84) and Mg# (0.91–0.94), respectively. They also display extremely low average abundances of CaO (0.13 wt.%) and Al₂O₃ (0.15 wt.%). The vast majority of the studied peridotites are also strongly depleted in REE. Simple batch and fractional melting models are not sufficient to explain their ultra-depleted composition. Whole-rock trace element abundances of the northern Vourinos mantle rocks can be modeled by up to 22–31% closed-system non-modal dynamic melting of an assumed primitive mantle (PM) source having spinel lherzolite composition. The highly depleted compositional signatures of the investigated peridotites indicate that they have experienced hydrous melting in the fore-arc mantle region above a SSZ. This intense melting event was responsible for the release of arc-related melts from the mantle. These melts reacted with the studied peridotites causing incongruent melting of pyroxenes followed by considerable olivine and Cr-spinel addition in terms of cryptic metasomatism. This later metasomatic episode has obscured any geochemical fingerprints indicative of an early mantle melting event in a MOR setting. The lack of any MOR-type peridotites in the northern Vourinos depleted mantle suite is quite uncommon for SSZ-type Neotethyan ophiolites.     

Alkaline series related to Early-Middle Miocene intra-continental rifting in a collision zone: An example from Polatlı, Central Anatolia,Turkey

A large volcanic area (∼7600 km²), the Galatean Volcanic Province (GVP), developed in northwest Central Anatolia during the Miocene along the Neo-Tethys Ocean suture zone possibly by post-collisional processes. The GVP mainly comprises 20–14 My old acid to intermediate volcanites with a geochemical signature indicating a mantle source modified by earlier (Late Cretaceous) subduction-related events. 100 km south of the GVP, near Polatlı, Ankara, basaltic rocks that cover large areas are intercalated with the Miocene deposits of the Beypazarı basin, an intra-continental subsidence zone at the southwest of the GVP. Field observations, geochemistry and K–Ar age dating of the Polatlı volcanites show that they are Early (19.9 Ma) to mid (14.1 Ma) Miocene in age, covering an area as large as 215 km². Variations in lava thickness and the thickness of the underlying silicified/baked zones suggest that the basaltic lavas erupted from a southern source, possibly from the Eskişehir fault zone, and flowed northwards. Most Polatlı samples have chemical compositions that indicate derivation from a mantle source with crustal contamination during ascent. They do not display any characteristic to suggest a subductional component. Although the GVP and Polatlı lavas formed close in time and space, they were derived from different mantle sources. Considering the positions of these two magmatic regions with regard to the Tethyan suture zone, we propose that the mantle beneath the GVP and near the suture zone memorised the earlier subduction while the mantle beneath Polatlı that is located about 100 km further from the suture zone remained apparently unchanged. After a significant volume of magma was consumed in the GVP, a later (∼10 My) and last activity (Güvem activity) has produced quantitatively much less basaltic rocks where this subductional signature seems to completely disappear. Considering that the western Anatolian crust is proposed to undergo extension since the Late Oligocene–Early Miocene times, the Early Miocene intra-plate Polatlı activity may have developed within this extensional tectonic regime. Combined with regional data, Polatlı data also provide broad estimations on how long a subductional event continues to modify the mantle after the subduction ceased (at least ∼20 My), how long the subductional signature is preserved during significant magmatism (between 6 and 10 My) and how far the subductional effect disappears laterally on the mantle with respect to the collision zone (<100 km).     

Water contents of the Cenozoic lithospheric mantle beneath the western part of the North China Craton: Peridotite xenolith constraints

Nominally anhydrous phases (clinopyroxene (cpx), orthopyroxene (opx), and olivine (ol)) of peridotite xenoliths hosted by the Cenozoic basalts from Beishan (Hebei province), and Fansi (Shanxi province), Western part of the North China Craton (WNCC) have been investigated by Fourier transform infrared spectrometry (FTIR). The H₂O contents (wt.) of cpx, opx and ol are 30–255 ppm, 14–95 ppm and ~ 0 ppm, respectively. Although potential H-loss during xenolith ascent cannot be excluded for olivine, pyroxenes (cpx and opx) largely preserve the H₂O content of their mantle source inferred from (1) the homogenous H₂O content within single pyroxene grains, and (2) equilibrium H₂O partitioning between cpx and opx. Based on mineral modes and assuming a partition coefficient of 10 for H₂O between cpx and ol, the recalculated whole-rock H₂O contents range from 6 to 42 ppm. In combination with previously reported data for other two localities (Hannuoba and Yangyuan from Hebei province), the H₂O contents of cpx, opx and whole-rock of peridotite xenoliths (43 samples) hosted by the WNCC Cenozoic basalts range from 30 to 654 ppm, 14 to 225 ppm, and 6 to 262 ppm respectively. The H₂O contents of the Cenozoic lithospheric mantle represented by peridotite xenoliths fall in a similar range for both WNCC and the eastern part of the NCC (Xia et al., 2010, Journal of Geophysical Research). Clearly, the Cenozoic lithospheric mantle of the NCC is dominated by much lower water content compared to the MORB source (50–250 ppm). The low H₂O content is not caused by oxidation of the mantle domain, and likely results from mantle reheating, possibly due to an upwelling asthenospheric flow during the late Mesozoic–early Cenozoic lithospheric thinning of the NCC. If so, the present NCC lithospheric mantle mostly represents relict ancient lithospheric mantle. Some newly accreted and cooled asthenospheric mantle may exist in localities close to deep fault.     

Heterogeneous distribution of water in the mantle beneath the central Siberian Craton: Implications from the Udachnaya Kimberlite Pipe

The paper presents new petrographic, major element and Fourier transform infrared (FTIR) spectroscopy data and PT-estimates of whole-rock samples and minerals of a collection of 19 relatively fresh peridotite xenoliths from the Udachnaya kimberlite pipe, which were recovered from its deeper levels. The xenoliths are non-deformed (granular), medium-deformed and highly deformed (porphyroclastic, mosaic-porphyroclastic, mylonitic) lherzolites, harzburgite and dunite. The lherzolites yielded equilibration temperatures (T) and pressures (P) ranging from 913 to 1324 °C and from 4.6 to 6.3 GPa, respectively. The non-deformed and medium-deformed peridotites match the 35 mW/m² conductive continental geotherm, whereas the highly deformed varieties match the 45 mW/m² geotherm. The content of water spans 2 ± 1–95 ± 52 ppm in olivine, 1 ± 0.5–61 ± 9 ppm in orthopyroxene, and 7 ± 2–71 ± 30 ppm in clinopyroxene. The amount of water in garnets is negligible. Based on the modal proportions of mineral phases in the xenoliths, the water contents in peridotites were estimated to vary over a wide range from < 1 to 64 ppm. The amount of water in the mantle xenoliths is well correlated with the deformation degree: highly deformed peridotites show highest water contents (64 ppm) and those medium-deformed and non-deformed contain ca. 1 ppm of H₂O. The high water contents in the deformed peridotites could be linked to metasomatism of relatively dry diamondiferous cratonic roots by hydrous and carbonatitic agents (fluids/melts), which may cause hydration and carbonation of peridotite and oxidation and dissolution of diamonds. The heterogeneous distribution of water in the cratonic mantle beneath the Udachnaya pipe is consistent with the models of mantle plume or veined mantle structures proposed based on a trace element study of similar xenolithic suits. Mantle metasomatism beneath the Siberian Craton and its triggered kimberlite magmatism could be induced by mantle enrichment in volatiles (H₂O, CO₂) supplied by numerous subduction zones which surrounded the Siberian continent in Neoproterozoic-Cambrian time.     

Post-collisional ultrapotassic rocks and mantle xenoliths in the Sailipu volcanic field of Lhasa terrane,south Tibet: Petrological and geochemical constraints on mantle source and geodynamic setting

Post-collisional ultrapotassic magmatic rocks (15.2–18.8 Ma), containing mantle xenoliths, are extensively distributed in the Sailipu volcanic field of the Lhasa terrane in south Tibet. They could be subdivided into high-MgO and low-MgO subgroups based on their petrological and geochemical characteristics. The high-MgO subgroup has olivine-I (Fo_87–92), phlogopite and clinopyroxene as phenocryst phases, while the low-MgO subgroup consists mainly of phlogopite, clinopyroxene and olivine-II (Fo_77–89). These ultrapotassic magmatic rocks have high MgO (4.6–14.5 wt%), Ni (145–346 ppm), Cr (289–610 ppm) contents, and display enrichment in light rare earth element (REE) over heavy REE and enriched large ion lithophile elements (LILE) relative to high field strength elements (HFSE) with strongly negative Nb-Ta-Ti anomalies in primitive mantle-normalized trace element diagrams. They have extremely radiogenic (⁸⁷Sr/⁸⁶Sr)_i (0.7167–0.7274) and unradiogenic (¹⁴³Nd/¹⁴⁴Nd)_i (0.5118–0.5120), high (²⁰⁷Pb/²⁰⁴Pb)_i (15.740–15.816) and (²⁰⁸Pb/²⁰⁴Pb)_i (39.661–39.827) at a given (²⁰⁶Pb/²⁰⁴Pb)_i (18.363–18.790) with high δ¹⁸O values (7.3–9.7‰). Strongly linear correlations between depleted mid-ocean ridge basalt-source mantle (DMM) and the Indian continental crust (HHCS) in Sr-Nd-Pb-O isotopic diagrams indicate that the geochemical features could result from reaction between mantle peridotite and enriched components (fluids and melts) released by the eclogitized Indian continental crust (HHCS) in the mantle wedge. The high-MgO (13.7–14.5 wt%) subgroup displays higher (¹⁴³Nd/¹⁴⁴Nd)_i, lower (⁸⁷Sr/⁸⁶Sr)_i and (²⁰⁶Pb/²⁰⁴Pb)_i ratios and lower δ¹⁸O values compared with the low-MgO (4.6–8.8 wt%) subgroup. High Ni (850–4862 ppm) contents of olivine phenocrysts and high whole-rock SiO₂, NiO, low CaO contents indicate that the low-MgO ultrapotassic magmatic rocks are derived from partial melting of olivine-poor mantle pyroxenite. However, lower Ni concentrations of olivine phenocryst and lower whole-rock SiO₂, NiO, higher CaO contents of the high-MgO ultrapotassic rocks may indicate their peridotite mantle source. This could be attributed to different amounts of silicate-rich components added into the mantle sources of the parental magmas in the mantle wedge caused by the northward subduction of the Indian continental lithosphere. The reaction-formed websterite xenoliths, reported for the first time in this study, are made up of anhedral and interlocking clinopyroxene (45–65 vol%) and orthopyroxene (30–50 vol%) with minor phlogopite (< 3 vol%) and quartz (< 2 vol%) and are suggested to be formed by silicate metasomatism of the mantle peridotite. The harzburgites, another major type of mantle xenolith in south Tibet, have a mineral assemblage of olivine (60–75 vol%), orthopyroxene (20–35 vol%), clinopyroxene (< 3 vol%), phlogopite (< 2 vol%) and spinel (< 2 vol%) and may have experienced subduction-related metasomatism. Combined with two types of ultrapotassic magmas, we propose that compositions of mantle wedge beneath south Tibet may gradually evolve from harzburgite through lherzolite to websterite with strong metasomatism of silicate-rich components in their mantle source region. Partial melting of the enriched mantle sources could be triggered by rollback of Indian continental slab during 25–8 Ma in south Tibet.     

Magmatic evolution of a dying spreading axis: Evidence for the interaction of tectonics and mantle heterogeneity from the fossil Phoenix Ridge,Drake Passage

New ⁴⁰Ar–³⁹Ar ages of 5.6 to 1.3 Ma for lavas from the fossil Phoenix Ridge in the Drake Passage show that magmatism continued for at least 2 Ma after the cessation of spreading at 3.3 ± 0.2 Ma. The Phoenix Ridge lavas are incompatible element-enriched relative to average MORB and show an increasing enrichment with decreasing age, corresponding to progressively decreasing degrees of partial melting of spinel peridotite after spreading stopped. The low-degree partial melts increasingly tap a mantle source with radiogenic Sr and Pb but unradiogenic Nd isotope ratios implying an ancient enrichment. The post-spreading magmas apparently form by buoyant ascent of enriched and easily fusible portions of the upper mantle. Only segments of fossil spreading ridges underlain by such enriched and fertile mantle show post-spreading volcanism frequently forming bathymetric highs. The Phoenix Ridge lavas belong to the Pacific, rather than the Atlantic, mantle domain in regional Sr–Nd–Pb space. Our new data show that the southern Pacific Ocean mantle is heterogeneous containing significant enriched portions that are preferentially tapped at low melt fractions. Isotopic mapping reveals that Pacific-type upper mantle flows eastward through Drake Passage and surrounds the subducting Phoenix Plate beneath the Bransfield Basin.     

Stability and crystal chemistry of iron-bearing dense hydrous magnesium silicates

Dense hydrous magnesium silicates (DHMS) are supposed to be key phases in planetary water cycles because of their ability to carry water to deep mantle regions in subduction slab environments. In order to understand water cycles in iron-enriched planetary systems such as Mars knowledge of the water content and stability of iron-bearing DHMS is required. Iron-bearing DHMS were synthesized based on two starting compositions, MgFeSiO₄ + 9.5 wt% H₂O system and a simple hydrous Martian mantle composition containing Fe, Mg, Al and Si + 12.35 wt% H₂O (hydrous FMAS system). Compared to literature data on phase D, iron-bearing phase D shows analogous variations in water contents as Mg-phase D but appears to be stable at higher temperatures than Mg-phase D for both starting compositions used in this study. Iron-bearing superhydrous phase B contains up to 7 wt% H₂O and shows an extended thermal stability in the hydrous FMAS system. The high-temperature stability of iron-bearing DHMS with a Mars-like bulk composition indicates that these hydrous phases could host significant amounts of water at core-mantle boundary conditions (1500 °C and 23 GPa) in a hydrous Martian mantle.     

Regimes of subduction and lithospheric dynamics in the Precambrian: 3D thermomechanical modelling

Comparing the early Earth to the present day, geological–geochemical evidence points towards higher mantle potential temperature and a different type of tectonics. In order to investigate possible changes in Precambrian tectonic styles, we conduct 3D high-resolution petrological–thermomechanical numerical modelling experiments for oceanic plate subduction under an active continental margin at a wide range of mantle potential temperature T_P (∆ T_P = 0 − 250 K, compared to present day conditions). At present day mantle temperatures (∆ T_P = 0 K), results of numerical experiments correspond to modern-style subduction, whereas at higher temperature conditions important systematic changes in the styles of both lithospheric deformation and mantle convection occur. For ∆ T_P = 50 − 100 K a regime of dripping subduction emerges which is still very similar to present day subduction but is characterised by frequent dripping from the slab tip and a loss of coherence of the slab, which suggests a close relationship between dripping subduction and episodic subduction. At further increasing ∆ T_P = 150 − 200 K dripping subduction is observed together with unstable dripping lithosphere, which corresponds to a transitional regime. For ∆ T_P = 250 K, presumably equivalent to early Archean, the dominating tectonic style is characterised by small-scale mantle convection, unstable dripping lithosphere, thick basaltic crust and small plates. Even though the initial setup is still defined by present day subduction, this final regime shows many characteristics of plume-lid tectonics. Transition between the two end-members, plume-lid tectonics and plate tectonics, happens gradually and at intermediate temperatures elements of both tectonic regimes are present. We conclude, therefore, that most likely no abrupt geodynamic regime transition point can be specified in the Earth's history and its global geodynamic regime gradually evolved over time from plume-lid tectonics into modern style plate tectonics.     

Postdated melting of subcontinental lithospheric mantle by the Emeishan mantle plume: Evidence from the Anyi intrusion,Yunnan, SW China

The Anyi intrusion is located in the central zone of Emeishan large igneous province (ELIP), SW China. It outcrops in an area of about 0.65 km² and ~ 1 km thick and dips to the southwest. The Anyi intrusion consists of a lower clinopyroxenite zone, middle gabbro zone, and an upper monzonite–syenite zone. Up to 400 m thick stratiform disseminated Fe–Ti oxide layer with grades of 16–18 wt.% total Fe is hosted in the lower clinopyroxenite zone. Zircon SHRIMP U–Pb age (247 ± 3 Ma) indicates that the Anyi intrusion represents postdated mafic magmatism resulting from the ~ 260 Ma Emeishan mantle plume. Compared with the typical oxide-bearing intrusions (such as Panzhihua and Baima) formed at ~ 260 Ma in the ELIP, the Anyi intrusion is characterized by high alkaline contents and LREE/HREE ratios, extremely low εNd values (− 6.2 to − 7.6) and moderate high (⁸⁷Sr/⁸⁶Sr)_i values (0.7072 to 0.7086). These characteristics of the Anyi intrusion cannot be explained by fractional crystallization or crustal contamination, but may reflect a unique enriched continental lithospheric mantle source (a mantle source mixed between garnet pyroxenite and spinel peridotite). We propose that the postdated mafic magmatism associated with the formation of the Anyi intrusion and its Fe–Ti oxide ore may be the product of melting of a mantle source mixed between garnet pyroxenite and spinel peridotite in the shallow lithosphere caused by conductive heating combined with lithosphere thinning due to plume–lithosphere interaction.     

Reduction of buried oxidized oceanic crust during subduction

The relationship among subducted oxidized oceanic crust and oxidation state of the subarc mantle, and arc magmas is one of the important aspects to evaluate convergent margin tectonics. However details of the oxidized mass transferred from buried oceanic crust to the overlying subarc mantle wedge remain obscure. Here we investigate the Songduo eclogites from south Tibet formed by the subduction of the paleo-Tethyan oceanic crust, and identify an abrupt decrease in pyrope and increase in almandine contents from the mantle to rim of garnet grains. This is coupled with a decrease in the Fe^3 + content of epidote and Fe^3 +/(Fe^2 ++ Fe^3 +) ratios from garnet core to rim domains, as well as speciation of calcite, a new mineral phase, in the rock matrix. Minor sulfates occur only as inclusions in garnet core domains, whereas sulfides are confined to the matrix as an accessory mineral phase. Aegirine augite occurs as relics or inclusions in garnet and omphacite. These features clearly suggest that oxidized components, Fe^3 + and S^6 +, were reduced as Fe^2 + and S^2 −, respectively, at the subduction zone. Thermodynamic modeling in the P–T-log₁₀f_O2 space using updated Perplex_X programs further revealed that the Songduo eclogites experienced oxygen fugacity variation of up to 8 log₁₀ units, with decreasing pressure. Petrological observations further suggest that the strong redox processes took place, after breaking of garnet, during the initial exhumation of the eclogites. CO₂ and minor sulfur are subsequently transferred from the cold oceanic subduction zone to the overlying mantle wedge, partially released by arc volcanoes to atmosphere. Our study presents a case of C and S recycling between the Earth's exterior and interior.