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.     

Long wavelength gravity anomalies over India: Crustal and lithospheric structures and its flexure

Long wavelength gravity anomalies over India were obtained from terrestrial gravity data through two independent methods: (i) wavelength filtering and (ii) removing crustal effects. The gravity fields due to the lithospheric mantle obtained from two methods were quite comparable. The long wavelength gravity anomalies were interpreted in terms of variations in the depth of the lithosphere–asthenosphere boundary (LAB) and the Moho with appropriate densities, that are constrained from seismic results at certain points. Modeling of the long wavelength gravity anomaly along a N–S profile (77°E) suggest that the thickness of the lithosphere for a density contrast of 0.05 g/cm³ with the asthenosphere is maximum of ∼190 km along the Himalayan front that reduces to ∼155 km under the southern part of the Ganga and the Vindhyan basins increasing to ∼175 km south of the Satpura Mobile belt, reducing to ∼155–140 km under the Eastern Dharwar craton (EDC) and from there consistently decreasing south wards to ∼120 km under the southernmost part of India, known as Southern Granulite Terrain (SGT).The crustal model clearly shows three distinct terrains of different bulk densities, and thicknesses, north of the SMB under the Ganga and the Vindhyan basins, and south of it the Eastern Dharwar Craton (EDC) and the Southern Granulite Terrain (SGT) of bulk densities 2.87, 2.90 and 2.96 g/cm³, respectively. It is confirmed from the exposed rock types as the SGT is composed of high bulk density lower crustal rocks and mafic/ultramafic intrusives while the EDC represent typical granite/gneisses rocks and the basement under the Vindhyan and Ganga basins towards the north are composed of Bundelkhand granite massif of the lower density. The crustal thickness along this profile varies from ∼37–38 km under the EDC, increasing to ∼40–45 km under the SGT and ∼40–42 km under the northern part of the Ganga basin with a bulge up to ∼36 km under its southern part. Reduced lithospheric and crustal thicknesses under the Vindhyan and the Ganga basins are attributed to the lithospheric flexure of the Indian plate due to Himalaya. Crustal bulge due to lithospheric flexure is well reflected in isostatic Moho based on flexural model of average effective elastic thickness of ∼40 km. Lithospheric flexure causes high heat flow that is aided by large crustal scale fault system of mobile belts and their extensions northwards in this section, which may be responsible for lower crustal bulk density in the northern part. A low density and high thermal regime in north India north of the SMB compared to south India, however does not conform to the high S-wave velocity in the northern part and thus it is attributed to changes in composition between the northern and the southern parts indicating a reworked lithosphere. Some of the long wavelength gravity anomalies along the east and the west coasts of India are attributed to the intrusives that caused the breakup of India from Antarctica, and Africa, Madagascar and Seychelles along the east and the west coasts of India, respectively.     

Lithospheric structure of the North China Craton: Integrated gravity,geoid and topography data

The lithospheric structure of ancient cratons provides important constraints on models relating to tectonic evolution and mantle dynamics. Here we present the 3D lithospheric structure of the North China Craton (NCC) from a joint inversion of gravity, geoid and topography data. The NCC records a prolonged history of Archean and Paleoproterozoic accretion of crustal blocks through subduction and collision building the cratonic architecture, which was subsequently differentially destroyed during Mesozoic through extensive magmatism. The thermal structure obtained in our study is considered to define the lithosphere-asthenosphere boundary (LAB) of the NCC, and reflects the density variations within the mantle lithosphere. Employing the Moho depths from deep seismic sounding profiles for the inversion, and based on repeated computations using different parameters, we estimate the Moho depth, LAB depth and average crustal density of the craton. The Moho depth varies from 28 to 50 km and the LAB depth varies from 105 to 205 km. The LAB and Moho show concordant thinning from West to East of the NCC. The average crustal density is 2870 kg m^− 3 in the western part of the NCC, higher than that in the eastern part (2750 kg m^− 3). The results of joint inversion in our study yielded LAB depth and lithospheric thinning features similar to those estimated from thermal and seismic studies, although our results show different depth and variations in the thickness. The lithosphere gently thins from 145 to 105 km in the eastern NCC, where as the thinning is much less pronounced in the western NCC with average depth of about 175 km. The joint inversion results in this study provide another perspective on the lithospheric structure from the density properties and corresponding geophysical responses in an ancient craton.     

Cumberland batholith,Trans-Hudson Orogen,Canada: Petrogenesis and implications for Paleoproterozoic crustal and orogenic processes

Large volume, plutonic belts, such as the ～ 221,000 km², ca. 1.865–1.845 Ga Cumberland batholith (CB) of the Trans-Hudson Orogen in Canada, are major components of Paleoproterozoic orogenic belts. In many cases, they have been interpreted as continental arc batholiths. The petrogenesis and tectonic context of the CB and implications for crustal growth and recycling are interpreted herein based on a 900 km geochemical-isotopic (Nd–O) transect across it and into granitoid plutons within bounding Archean cratons in central and southern Baffin Island.The mainly granulite grade CB, emplaced over an age span of between 14 and 24 Ma, consists mainly of high-K to shoshonitic monzogranite and granodiorite, but also includes low- and medium-K granitoid rocks. Metaluminous to slightly peraluminous compositions and δ¹⁸O (VSMOW) values (+ 6 to + 10‰) indicate derivation from infracrustal (I-type) sources. ε_Nd 1.85 Ga signatures (? 12 to ? 2) of both mafic and felsic units suggest a dominance of evolved sources. Isotopic signatures in the interior of the CB (? 2 to ? 7) are more radiogenic than those within Archean domains in central (? 8 to ? 15) and southern (? 5 to ? 19) Baffin Island. The isotopic transect is interpreted as ‘imaging’ an accreted microcontinental block (Meta Incognita) and bounding Archean cratons. The CB includes granites of arc, within-plate (A-type) and post-collisional affinity and volumetrically minor mafic rocks with both arc and non-arc features. (La/Yb)_CN and Sr/Y values range from < 1 to 225 and < 1 to 611, respectively. In these respects, some CB granitoid rocks resemble Paleozoic adakitic granites, interpreted as partial melts of greatly thickened crust within post-collisional settings, such as Tibet. Thus, the CB likely encompasses various non-consanguineous magmatic suites generated at deep- to mid-crustal depths. Although CB granitoid rocks undoubtedly had important crustal sources, it is hard to assess the relative contribution of mantle-derived magmas.The CB is best interpreted as a post-accretion batholith resulting from large-scale lithospheric mantle delamination followed by the upwelling of hot asthenospheric mantle leading to voluminous crustal partial melting. Contributors to crustal instability which may have facilitated such delamination included: (a) a collage of recently assembled small cratons underlain by hot, weak lithosphere with mantle-depth structural breaks within this segment of the Trans-Hudson Orogen; (b) the gabbro-eclogite phase transformation, and (c) a greatly thickened crustal section (> 60 km), as evidenced by adakitic granites.     

Highly heterogeneous lithospheric mantle beneath the Central Zone of the North China Craton evolved from Archean mantle through diverse melt refertilization

High-Mg# peridotite xenoliths in the Cenozoic Hebi basalts from the North China Craton have refractory mineral compositions (Fo > 91.5) and highly heterogeneous Sr–Nd isotopic compositions (⁸⁷Sr/⁸⁶Sr = 0.7031–0.7048, ¹⁴³Nd/¹⁴⁴Nd = 0.5130–0.5118) ranging from MORB-like to EM1-type mantle, which are similar to those of peridotites from Archean cratons. Thus, the high-Mg# peridotites may represent relics of the ancient lithospheric mantle. Published Re–Os isotopic data for Cenozoic basalt-borne xenoliths show T_RD ages of 3.0–1.5 Ga for the peridotites from Hebi (the center of the craton), 2.2–0 Ga for those from Hannuoba and Jining (north margin of the craton), and 2.6–0 Ga for those from Fanshi and Yangyuan (midway between the center and north margin of the craton). In situ Re–Os data of sulfides in Hannuoba peridotites suggest that whole-rock Re–Os model ages represent mixtures of multiple generations of sulfides with varying Os isotopic compositions. These observations indicate that initial lithospheric mantle beneath the Central Zone of the North China Craton formed during the Archean and was refertilized by multiple melt additions after its formation. The refertilization became more intensive from the interior to the margin of the craton, leading to the high heterogeneity of the lithospheric mantle: more ancient and refractory peridotites with highly variable Sr–Nd isotopic compositions in the interior, and more young and fertile peridotites with depleted Sr–Nd isotopic composition in the margin. Our data, coupled with published petrological and geochemical data of peridotites from the Central Zone of the North China Craton, suggest that the lithospheric mantle beneath this region is highly heterogeneous, likely produced by refertilization of Archean mantle via multiple additions of melts/fluids, which were closely related to the Paleoproterozoic collision between the Eastern and the Western Blocks and subsequent circum-craton subduction events.     

Inhomogeneous lithospheric thinning in the central North China Craton: Zircon U–Pb and S–He–Ar isotopic record from magmatism and metallogeny in the Taihang Mountains

The large scale Mesozoic magmatism and related metallogeny in the Taihang Mountains (TM) provide important clues for the lithospheric thinning of the North China Craton (NCC). Among the ore deposits, the vein gold mineralization of Shihu in the Fuping region and the skarn ore deposit of Xishimen in the Wu'an region represent typical Mesozoic metallogeny in the TM. In the Shihu gold mine, the Mapeng batholith is dominantly composed of monzogranite and granodiorite, whereas, the Wu'an pluton in the Xishimen iron mine mainly comprises monzonite and diorite. Here we present zircon LA–ICP-MS U–Pb data from 8 samples which reveal the timing of magmatism in the TM as ca. 130 Ma, which is contemporaneous with the large-scale metallogeny in the margins of the NCC. The δ³⁴S values recorded in the sulfide minerals from the Shihu gold deposit and the Xishimen skarn iron deposit show a range of 2.2‰–5.0‰, and 11.6‰–18.7‰, respectively. Helium isotopic compositions of fluid inclusions in pyrite from the Shihu gold deposit vary from 0.12 to 1.98 Ra (where Ra is the ³He/⁴He ratio of air = 1.39 × 10^? 6), with calculated mantle helium values of 1.4%–25%, whereas, those of the Xishimen skarn iron deposit range from 0.06 to 0.19 Ra, with calculated mantle helium of 0.7%–2.2%. The S–He–Ar isotopic data suggest a lower crustal origin for the ore-forming components, with variable inputs of mantle source. The large population of inherited zircons in our samples, with ²⁰⁷Pb/²⁰⁶Pb ages ranging between 2500 Ma and 1800 Ma, also supports crustal participation. Our data reveal that the Shihu gold deposit witnessed greater mantle input than the Xishimen skarn iron deposit, suggesting that the continental lithosphere is markedly thinner under the Fuping region than that under the Wu'an region. Our interpretation is also supported by published data from two ultra-broadband high-precision magnetotelluric sounding profiles across the TM region showing a variation in the lithosphere thickness from 155 km to 70 km while moving from the south (Wu'an region) to the north (Fuping region). Our study suggests that inhomogeneous lithospheric thinning in the central NCC occurred at least as early as ca. 130 Ma ago.     

Subduction of Indian continent beneath southern Tibet in the latest Eocene (~ 35 Ma): Insights from the Quguosha gabbros in southern Lhasa block

Geophysical data illustrate that the Indian continental lithosphere has northward subducted beneath the Tibet Plateau, reaching the Bangong–Nujiang suture in central Tibet. However, when the Indian continental lithosphere started to subduct, and whether the Indian continental crust has injected into the mantle beneath southern Lhasa block, are not clear. Here we report new results from the Quguosha gabbros of southern Lhasa block, southern Tibet. LA-ICP-MS zircon U–Pb dating of two samples gives a ca. 35 Ma formation age (i.e., the latest Eocene) for the Quguosha gabbros. The Quguosha gabbro samples are geochemically characterized by variable SiO₂ and MgO contents, strongly negative Nb–Ta–Ti and slightly negative Eu anomalies, and uniform initial ⁸⁷Sr/⁸⁶Sr (0.7056–0.7058) and ε_Nd(t) (− 2.2 to − 3.6). They exhibit Sr–Nd isotopic compositions different from those of the Jurassic–Eocene magmatic rocks with depleted Sr–Nd isotopic characteristics, but somewhat similar to those of Oligocene–Miocene K-rich magmatic rocks with enriched Sr–Nd isotopic characteristics. We therefore propose that an enriched Indian crustal component was added into the lithospheric mantle beneath southern Lhasa by continental subduction at least prior to the latest Eocene (ca. 35 Ma). We interpret the Quguosha mafic magmas to have been generated by partial melting of lithospheric mantle metasomatized by subducted continental sediments, which entered continental subduction channel(s) and then probably accreted or underplated into the overlying mantle during the northward subduction of the Indian continent. Continental subduction likely played a key role in the formation of the Tibetan plateau at an earlier date than previously thought.     

Upper mantle structure of the Saharan Metacraton

The ∼500,000 km² Saharan Metacraton in northern Africa (metacraton refers to a craton that has been mobilized during an orogenic event but that is still recognisable through its rheological, geochronological and isotopic characteristics) is an Archean–Paleoproterozoic cratonic lithosphere that has been destabilized during the Neoproterozoic. It extends from the Arabian–Nubian Shield in the east to the Trans-Saharan Belt in the west, and from the Oubanguides Orogenic Belt in the south to the Phanerozoic cover of North Africa. Here, we show that there are high S-wave velocity anomalies in the upper 100 km of the mantle beneath the metacraton typical of cratonic lithosphere, but that the S-wave velocity anomalies in the 175–250 km depth are much lower than those typical of other cratons. Cratons have possitive S-wave velocity anomalies throughout the uppermost 250 km reflecting the presence of well-developed cratonic root. The anomalous upper mantle structure of the Saharan Metacraton might be due to partial loss of its cratonic root. Possible causes of such modification include mantle delamination or convective removal of the cratonic root during the Neoproterozoic due to collision-related deformation. Partial loss of the cratonic root resulted in regional destabilization, most notably in the form of emplacement of high-K calc-alkaline granitoids. We hope that this work will stimulate future multi-national research to better understand this part of the African Precambrian. Specifically, we call for efforts to conduct systematic geochronological, geochemical, and isotopic sampling, deploy a reasonably-dense seismic broadband seismic network, and conduct systematic mantle xenoliths studies.     

Geochronology and geochemistry of Mesozoic mafic–ultramafic complexes in the southern Liaoning and southern Jilin provinces,NE China: Constraints on the spatial extent of destruction of the North China Craton

LA–ICP–MS zircon U–Pb ages, geochemical and Sr–Nd–Pb isotope data are presented for mafic–ultramafic complexes from the southern Liaoning–southern Jilin area with the aim of determining the nature of the Mesozoic lithospheric mantle and to further constrain the spatial extent of destruction of the North China Craton (NCC). The complexes consist of olivine-websterite, gabbro, dolerite, and gabbro-diorite. Zircons from the complexes show typical zoning absorption, are euhedral–subhedral in shape, and yield high Th/U ratios (1.23–2.87), indicating a magmatic origin. Zircon U–Pb age data indicate that they formed in the Early Cretaceous (129–137 Ma). Geochemically, they have SiO₂ = 44.3–49.8%, MgO = 6.8–26.5%, Cr = 102–3578 ppm, and Ni = 31–1308 ppm, and are characterized by enrichment in large ion lithophile elements (LILEs) and light rare earth elements (LREEs), and depletion in high field strength elements (HFSEs) and heavy rare earth elements (HREEs), as well as a wide range of Sr–Nd–Pb isotopic compositions [(⁸⁷Sr/⁸⁶Sr)i = 0.70557–0.71119; ε_Nd (t) = ?5.4 to ?20.1; (²⁰⁶Pb/²⁰⁴Pb)i = 15.13–17.85; Δ7/4 = ?11.49 to 16.00; Δ8/4 = 102.64–203.48]. Compared with the southern Liaoning mafic–ultramafic rocks, the southern Jilin mafic–ultramafic rocks have high TiO₂ and Al₂O₃ contents, high ε_Nd (t) values, low (La/Yb)_N values, low initial ⁸⁷Sr/⁸⁶Sr ratios, and low radiogenic Pb isotopic compositions. These findings indicate that the primary magmas of the southern Jilin complexes were derived from lithospheric mantle that was previously metasomatized by a melt derived from the delaminated ancient lower crust, whereas the primary magmas of the southern Liaoning complexes originated from partial melting of a lithospheric mantle source that was previously modified by melt derived from the broken-off Yangtze slab. Therefore, the lateral extent of the NCC destruction should include the southern Liaoning–southern Jilin area.     

Origin of high-velocity anomalies beneath the Siberian craton: A fingerprint of multistage magma underplating since the Neoarchean

Despite the violent eruption of the Siberian Traps at ~ 250 Ma, the Siberian craton has an extremely low heat flow (18–25 mW/m²) and a very thick lithosphere (300–350 km), which makes it an ideal place to study the influence of mantle plumes on the long-term stability of cratons. Compared with seismic velocities of rocks, the lower crust of the Siberian craton is composed mainly of mafic granulites and could be rather heterogeneous in composition. The very high V_p (> 7.2 km/s) in the lowermost crust can be fit by a mixture of garnet granulites, two-pyroxene granulites, and garnet gabbro due to magma underplating. The high-velocity anomaly in the upper mantle (V_p = 8.3-8.6 km/s) can be interpreted by a mixture of eclogites and garnet peridotites. Combined with the study of lower crustal and mantle xenoliths, we recognized multistage magma underplating at the crust-mantle boundary beneath the Siberian craton, including the Neoarchean growth and Paleoproterozoic assembly of the Siberian craton beneath the Markha terrane, the Proterozoic collision along the Sayan-Taimyr suture zone, and the Triassic Siberian Trap event beneath the central Tunguska basin. The Moho becomes a metamorphism boundary of mafic rocks between granulite facies and eclogite facies rather than a chemical boundary that separates the mafic lower crust from the ultramafic upper mantle. Therefore, multistage magma underplating since the Neoarchean will result in a seismic Moho shallower than the petrologic Moho. Such magmatism-induced compositional change and dehydration will increase viscosity of the lithospheric mantle, and finally trigger lithospheric thickening after mantle plume activity. Hence, mantle plumes are not the key factor for craton destruction.     

Mineral chemistry and zircon geochronology of xenocrysts and altered mantle and crustal xenoliths from the Aries micaceous kimberlite: Constraints on the composition and age of the central Kimberley Craton,Western Australia

The Neoproterozoic (∼ 820 Ma) Aries micaceous kimberlite intrudes the central Kimberley Basin, northern Western Australia, and has yielded a suite of 27 serpentinised ultramafic xenoliths, including spinel-bearing and rare, metasomatised, phlogopite–biotite and rutile-bearing types, along with minor granite xenoliths. Proton-microprobe trace-element analysis of pyrope and chromian spinel grains derived from heavy mineral concentrates from the kimberlite has been used to define a ∼ 35–40 mW/m² Proterozoic geotherm for the central Kimberley Craton. Lherzolitic chromian pyrope highly depleted in Zr and Y, and Cr-rich magnesiochromite xenocrysts (class 1), probably were derived from depleted garnet peridotite mantle at ∼ 150 km depth. Sampling of shallower levels of the lithospheric mantle by kimberlite magmas in the north and north-extension lobes entrained high-Fe chromite xenocrysts (class 2), and aluminous spinel-bearing xenoliths, where both spinel compositions are anomalously Fe-rich for spinels from mantle xenoliths. This Fe-enrichment may have resulted from Fe–Mg exchange with olivine during slow cooling of the peridotite host rocks. Fine exsolution rods of aluminous spinel in diopside and zircon in rutile grains in spinel- and rutile-bearing serpentinised ultramafic xenoliths, respectively, suggest nearly isobaric cooling of host rocks in the lithospheric mantle, and indicate that at least some aluminous spinel in spinel-facies peridotites formed through exsolution from chromian diopside. Fe–Ti-rich metasomatism in the spinel-facies Kimberley mantle probably produced high-Ti phlogopite–biotite + rutile and Ti, V, Zn, Ni-enriched aluminous spinel ± ilmenite associations in several ultramafic xenoliths. U–Pb SHRIMP ²⁰⁷Pb/²⁰⁶Pb zircon ages for one granite (1851 ± 10 Ma) and two serpentinised ultramafic xenoliths (1845 ± 30 Ma; 1861 ± 31 Ma) indicate that the granitic basement and lower crust beneath the central Kimberley Basin are at least Palaeoproterozoic in age. However, Hf-isotope analyses of the zircons in the ultramafic xenoliths suggest that the underlying lithospheric mantle is at least late Archean in age.     

Density heterogeneity of the cratonic lithosphere: A case study of the Siberian Craton

Using free-board modeling, we examine a vertically-averaged mantle density beneath the Archean–Proterozoic Siberian Craton in the layer from the Moho down to base of the chemical boundary layer (CBL). Two models are tested: in Model 1 the base of the CBL coincides with the LAB, whereas in Model 2 the base of the CBL is at a 180 km depth. The uncertainty of density model is < 0.02 t/m³ or < 0.6% with respect to primitive mantle. The results, calculated at in situ and at room temperature (SPT) conditions, indicate a heterogeneous density structure of the Siberian lithospheric mantle with a strong correlation between mantle density variations and the tectonic setting. Three types of cratonic mantle are recognized from mantle density anomalies. ‘Pristine’ cratonic regions not sampled by kimberlites have the strongest depletion with density deficit of 1.8–3.0% (and SPT density of 3.29–3.33 t/m³ as compared to 3.39 t/m³ of primitive mantle). Cratonic mantle affected by magmatism (including the kimberlite provinces) has a typical density deficit of 1.0–1.5%, indicative of a metasomatic melt-enrichment. Intracratonic sedimentary basins have a high density mantle (3.38–3.40 t/m³ at SPT) which suggests, at least partial, eclogitization. Moderate density anomalies beneath the Tunguska Basin imply that the source of the Siberian LIP lies outside of the Craton. In situ mantle density is used to test the isopycnic condition of the Siberian Craton. Both CBL thickness models indicate significant lateral variations in the isopycnic state, correlated with mantle depletion and best achieved for the Anabar Shield region and other intracratonic domains with a strongly depleted mantle. A comparison of synthetic Mg# for the bulk lithospheric mantle calculated from density with Mg# from petrological studies of peridotite xenoliths from the Siberian kimberlites suggests that melt migration may produce local patches of metasomatic material in the overall depleted mantle.     

Thermochemical structure of the North China Craton from multi-observable probabilistic inversion: Extent and causes of cratonic lithosphere modification

We present a 3D thermochemical model of the North China Craton (NCC) from the surface down to 350 km by jointly inverting surface wave phase velocity data, geoid height, surface heat flow and absolute elevation with a multi-observable probabilistic inversion method. Our model reveals a thin (~ 65–100 km) and chemically fertile lithosphere (87 < Mg# < 90) beneath the Eastern NCC, consistent with independent results from mantle xenoliths, and supports the idea that the Eastern NCC experienced significant lithospheric destruction and refertilization during the Phanerozoic. In contrast, beneath the Trans-North China Orogen, Inner Mongolia Suture Zone and Yinshan belt, we observe a more heterogeneous (chemically and thermally) lithosphere, indicating that these areas have been partly involved in lithospheric modification and mechanical erosion at multiple scales. A cold and chemically refractory (Mg# > 90) lithospheric mantle is imaged beneath the central TNCO and Ordos Block, reaching depths > 260 km. This lithospheric “keel” is surrounded to the east by a high-temperature sublithospheric anomaly that originates at depths > 280 km. The spatial distribution of this anomaly and its correlation with the location of recent volcanism in the region suggest that the anomaly represents a deep mantle upwelling being diverted by the cratonic keel and spreading onto regions of shallow lithosphere. Our results indicate that the present-day thermochemical structure beneath the NCC is the result of a complex interaction between a large-scale return flow associated with the subduction of the Pacific slab and the shallow lithospheric structure.     

Anomalous crustal and lithospheric mantle structure of southern part of the Vindhyan Basin and its geodynamic implications

Tectonically active Vindhyan intracratonic basin situated in central India, forms one of the largest Proterozoic sedimentary basins of the world. Possibility of hydrocarbon occurrences in thick sediments of the southern part of this basin, has led to surge in geological and geophysical investigations by various agencies. An attempt to synthesize such multiparametric data in an integrated manner, has provided a new understanding to the prevailing crustal configuration, thermal regime and nature of its geodynamic evolution. Apparently, this region has been subjected to sustained uplift, erosion and magmatism followed by crustal extension, rifting and subsidence due to episodic thermal interaction of the crust with the hot underlying mantle. Almost 5–6 km thick sedimentation took place in the deep faulted Jabera Basin, either directly over the Bijawar/Mahakoshal group of mafic rocks or high velocity-high density exhumed middle part of the crust. Detailed gravity observations indicate further extension of the basin probably beyond NSL rift in the south. A high heat flow of about 78 mW/m² has also been estimated for this basin, which is characterized by extremely high Moho temperatures (exceeding 1000 °C) and mantle heat flow (56 mW/m²) besides a very thin lithospheric lid of only about 50 km. Many areas of this terrain are thickly underplated by infused magmas and from some segments, granitic–gneissic upper crust has either been completely eroded or now only a thin veneer of such rocks exists due to sustained exhumation of deep seated rocks. A 5–8 km thick retrogressed metasomatized zone, with significantly reduced velocities, has also been identified around mid to lower crustal transition.     

Refertilization of lithospheric mantle beneath the Yangtze craton in south-east China: Evidence from noble gases geochemistry

The Yangtze craton (YC), in eastern China, is one of the oldest cratons in the world and is characterized by a complex tectonic and geodynamic evolution. This evolution regards most of the eastern China craton, which since Mesozoic time has undergone significant thinning (> 200 km) of Archean lithosphere. This thinning favored the refertilization of the old refractory subcontinental lithospheric mantle (SCLM) by the upwelling of younger fertile asthenosphere. Whether this feature is localized only beneath certain areas of eastern China or is a more widespread characteristic of the mantle, including the YC, is a matter of debate.In order to constrain the history of the YC SCLM, we have measured the He- and Ar-isotopic compositions of fluid inclusions hosted in mantle xenoliths in the Lianshan area, which is part of the poorly investigated YC in south-east China. We also report new mineral chemistry and trace element compositions of clinopyroxenes from the same suite of samples, for comparison with noble gases. Two distinct types of xenoliths can be identified: Type 1, characterized by mantle-like He-isotopic (³He/⁴He) ratios (up to 9.1 Ra), represents fragments of a fertile lithospheric mantle; Type 2, showing ³He/⁴He values in the SCLM range (³He/⁴He < 7 Ra), represents shallow relicts of a refractory mantle. The patterns of rare-earth elements as well as the Y and Yb concentrations in the clinopyroxenes normalized to primitive mantle (Y_N and Yb_N, respectively) indicate that fractional partial melting might have affected the local mantle by < 3% in Type 1 and up to 20% in Type 2 xenoliths from Lianshan, respectively. The range of ⁴He/⁴⁰Ar* (⁴⁰Ar* is corrected for atmospheric contamination) ranges from 4.9 × 10^− 4 to 3.6 × 10^− 1, which is below the typical production ratio of the mantle (⁴He/⁴⁰Ar* = 1–5); this range is however compatible with this fractional partial melting. The variable ³He/⁴He and ⁴He/⁴⁰Ar* values in Lianshan xenoliths suggest that the local mantle source was also influenced by kinetic fractionation, possibly triggered by metasomatic melts. Metasomatism associated with carbonatitic melts, together with fluxing by CO₂-rich fluids, have permeated the mantle beneath Lianshan, generating the observed decoupling between noble gases and trace elements. The interpretative framework is also applicable for other mantle xenoliths from eastern China, indicating that the refertilization of the SCLM by ascending mantle-like melts is common also to YC, which can be identified using noble gases.     

Coexistence of the moderately refractory and fertile mantle beneath the eastern Central Asian Orogenic Belt

Relative to the North China Craton, the subcontinental lithospheric mantle (SCLM) beneath the Central Asian Orogenic Belt is little known. Mantle-derived peridotite xenoliths from the Cenozoic basalts in the Xilinhot region, Inner Mongolia, provide samples of the lithospheric mantle beneath the eastern part of the belt. The xenoliths are predominantly lherzolites with minor harzburgites, and can be subdivided into three groups, based on the REE patterns of clinopyroxenes. Group 1 peridotites (LREE-enriched), with low modal Cpx (3–7%), high Mg^# in olivine (> 90.6) and Cr^# in spinel (> 43.8), low whole-rock CaO + Al₂O₃ contents (1.62–3.22 wt.%) and estimated temperatures of 1043–1126 °C, represent moderately refractory SCLM that has experienced carbonatite-related metasomatism. Group 2 peridotites (LREE-depleted), with high modal Cpx (9–13%), low Mg^# in olivine (< 90.6) and Cr^# in spinel (< 20.0), high whole-rock CaO + Al₂O₃ contents (4.93–6.37 wt.%) and estimated temperatures of 814–970 °C, show affinity with Phanerozoic fertile SCLM that has undergone silicate-related metasomatism. Group 3 peridotites (convex-upward REE patterns), show wide ranges of olivine-Mg^# (88.4–90.6), spinel-Cr^# (11.5–47.6), and modal Cpx (3–14%) that overlap Groups 1 and 2. Their spinels have high TiO₂ contents (> 0.41 wt.%), implying involvement of reactions between melt and peridotites. The estimated temperatures of Group 3 (1033–1156 °C) are similar to those of Group 1. We suggest that the pre-existing moderately refractory lithospheric mantle (i.e., Group 1) beneath the eastern part of the Central Asian Orogenic Belt was strongly penetrated by upwelling asthenospheric material, and the cooling of this material produced fertile lithospheric mantle (i.e., Group 2). The present lithospheric mantle of this area consists of interspersed volumes of younger fertile and older more refractory lithosphere, with the fertile type dominating the shallower levels of the mantle.     

Eocene potassic and ultrapotassic volcanism in south Tibet: New constraints on mantle source characteristics and geodynamic processes

In the Yangbajing area, southern Tibet, several monogenic volcanoes were conformably superimposed on the Linzizong calc-alkaline volcanic successions. According to their petrologic and geochemical characteristics, these monogenic volcanoes are composed of three rock varieties: tephritic phonolitic plugs and shoshonitic and trachytic lavas. Their geochemical systematics reveals that low-pressure evolutionary processes in the large voluminous Linzizong calc-alkaline magmas were not responsible for the generation of these potassic–ultrapotassic rocks, but the significant change in petrologic and geochemical characteristics from the Linzizong calc-alkaline to potassic–ultrapotassic magma is likely accounted for the change of metasomatic agents in the southern Tibetan lithospheric mantle source during the Paleocene to Eocene. The tephritic phonolites containing both leucite and plagioclase show primary ultrapotassic character similar to that of Mediterranean plagioleucititic magmas. Radiogenic Sr increases with SiO₂ in the xenolith-bearing trachytes strongly suggesting significant crustal assimilation in the shoshonitic magmas. The Yangbajing ultrapotassic rocks have high K₂O and Al₂O₃, and show depletion of high field strength elements (HFSEs) with respect to large ion lithophile elements. In primitive mantle-normalized element diagrams, all samples are characterized by positive spikes at Th (U) and Pb with negative anomalies at Ba, Nb–Ta and Ti, reflecting the orogenic nature of the ultrapotassic rocks. They are characterized by highly radiogenic ⁸⁷Sr/⁸⁶Sr_(i) ratios (0.7061–0.7063) and unradiogenic ¹⁴³Nd/¹⁴⁴Nd_(i) (0.5125), and Pb isotopic compositions (²⁰⁶Pb/²⁰⁴Pb = 18.688–18.733, ²⁰⁷Pb/²⁰⁴Pb = 15.613–15.637, and ²⁰⁸Pb/²⁰⁴Pb = 38.861–38.930) similar to the global subducting sediment. Strong enrichment of incompatible trace elements and high Th fractionation from the other HFSEs (such as Nb and U) clearly indicate that the Th-enriched sedimentary component in a network veined mantle source was mainly introduced by sediment-derived melts. In addition, the ultrapotassic rocks have significant Ce (Ce/Ce* = 0.77–0.84) and Eu (Eu/Eu* = 0.72–0.75) anomalies, suggesting a subduction sediment input into the southern Tibetan lithospheric mantle source. In contrast, high U/Th (> 0.20) and Ba/Th (> 32) and low Th/La (< 0.3) in the shoshonites indicate that the Eocene potassic magma originated from partial melting of the surrounding peridotite mantle pervasively affected by slab-related fluid addition from the dehydration of either the subducting oceanic crust or the sediment. Thus, at least two different subduction-related metasomatic agents re-fertilized the upper mantle. According to the radiometric ages and spatial distribution, the Gangdese magmatic association shows a temporal succession from the Linzizong calc-alkaline to ultrapotassic magmas. This indicates a late arrival of recycled sediments within the Tibetan lithospheric mantle wedge. The most diagnostic signatures for the involvement of continent-derived materials are the super-chondritic Zr/Hf (45.5–49.2) and elevated Hf/Sm values (0.81–0.91) in the ultrapotassic rocks. Therefore, the occurrence of orogenic magmatism in the Gangdese belt likely represents the volcanic expression of the onset of the India–Asia collision, preceding the 10 Ma Neo-Tethyan slab break-off process at 42–40 Ma. The absence of residual garnet in the mantle source for the ultrapotassic volcanism seems to imply that the southern Tibetan lithosphere was not been remarkably thickened until the Eocene (～ 50 Ma).     

Source of metals in the Guocheng gold deposit,Jiaodong Peninsula,North China Craton: Link to early Cretaceous mafic magmatism originating from Paleoproterozoic metasomatized lithospheric mantle

Widespread Mesozoic Au and other hydrothermal polymetal (Zn–Pb–Cu–Mo–Ag–W–Fe–REE) deposits or smaller prospects occur in association with ancient mobile belts surrounding and cutting through the North China Carton (NCC). Among these, the gold ores of the Jiaodong Peninsula, Shandong Province, eastern NCC, represent the largest gold district in China. However, the genesis of these important gold mineralizations has remained controversial, notably their relationships to widespread mafic magmatism of alkaline affinity.The ore bodies of the Guocheng gold deposit on the Jiaodong Peninsula are fracture-controlled, sulfide-rich veins and disseminations, formed contemporaneously with abundant dolerite, lamprophyre and monzonite dikes at ca. 120 Ma. Dolerite dikes possess mantle-like major element compositions and alkaline affinity, associated with prominent subduction-type trace element enrichments. The dikes show petrographic and chemical evidence of magma mixing that triggered exsolution of magmatic sulfide and anhydrite crystallization, preserved as primary inclusions in phenocrysts. LA-ICP-MS analysis of magmatic sulfide inclusions demonstrates that metal abundance ratios (Ag, As, Au, Bi, Co, Cu, Mo, Ni, Pb, Sb, Zn) largely correspond to those of both unaltered bulk rock and bulk ore. Together with identical Pb isotope ratios of dolerite and bulk ore, this demonstrates that gold mineralization and dolerite dikes share a common source.Lead isotope signatures of the ore sulfides are much less radiogenic (17.08 < ²⁰⁶Pb/²⁰⁴Pb < 17.25, 15.41 <²⁰⁷Pb/²⁰⁴Pb < 15.45, 37.55 < ²⁰⁸Pb/²⁰⁴Pb < 37.93) relative to the Pb signature of Phanerozoic convecting mantle and plot to the left of the Geochron and above the MORB-source mantle Pb evolution line. Forward Monte Carlo simulations indicate three events for the U–Th–Pb isotope evolution: (1) late Archean formation of juvenile crust is followed by (2) subduction of this aged crust at ca. 1.85 Ga along with the assembly of Jiao–Liao–Ji mobile belt (suture within Columbia supercontinent). This late-Archean subducted crust released fluids with drastically reduced U/Pb that metasomatized the overlying depleted mantle, which formed cratonic lithospheric mantle. This metasomatized lithospheric mantle was (3) tapped in response to early Cretaceous extensional tectonics affecting notably the eastern margin of the NCC to generate mafic magmas and associated gold mineralization at Guocheng. Similarly non-radiogenic uranogenic Pb isotope data characterize the contemporaneous mafic dikes and gold deposits in the entire Jiaodong Peninsula, suggesting that our genetic model applies to the entire Jiaodong gold district.We propose that early Cretaceous melting of subcontinental lithospheric mantle metasomatized by subduction fluids during Paleoproterozoic amalgamation of terranes to the eastern NCC along with Columbia supercontinent assembly generated mafic magmatism and associated gold deposits. Given the conspicuous association of Phanerozoic hydrothermal ore deposits associated with reactivated Paleoproterozoic mobile belts, we envisage that our genetic model, which largely corresponds to that which is proposed for the Bingham porphyry-Cu–Au–Mo deposit, USA, may explain much of the magmatic-hydrothermal activity and associated ore formation all around the NCC.     

The alkaline-carbonatite complex of Jacupiranga (Brazil): Magma genesis and mode of emplacement

A comprehensive study including new field, petrological and geochemical data is reported on the Jacupiranga alkaline-carbonatite complex (133–131 Ma) which, together with other alkaline complexes, occurs in southern Brazil and is coeval with the Paraná CFB province. It consists of a shallow intrusion (ca. 65 km²) in the Precambrian crystalline basement, and can be subdivided in two main diachronous plutonic bodies: an older dunite-gabbro-syenite in the NW and a younger clinopyroxenite-ijolite (s.l.) in the SE, later injected by a carbonatitic core (ca. 1 km²). An integrated petrogenetic model, based on bulk rock major and trace element analyses, mineral chemistry and Sr-Nd-Pb-C isotopic data, suggests that the two silicate intrusions generated from different mantle-derived magmas that evolved at shallow level (2–3 km depth) in two zoned cup-shaped plutonic bodies growing incrementally from independent feeding systems. The first intrusion was generated by OIB-like alkaline to mildly alkaline parental basalts that initially led to the formation of a dunitic adcumulate core, discontinuously surrounded by gabbroic cumulates, in turn injected by subanular syenite intrusive and phonolite dykes. Nephelinitic (± melilite) melts – likely generated deep in the lithosphere at ≥ 3 GPa – were the parental magmas of the second intrusion and gave rise to large coarse-grained clinopyroxenite ad- to meso-cumulates, in turn surrounded, and partially cut, by semi-annular fine-layered melteigite-ijolite-urtite ortho-cumulates. The available isotopic data do not evidence genetic links between carbonatites and the associated silicate intrusions, thus favouring an independent source from the mantle. Moreover, it may be suggested that, unlike gabbro-syenites and carbonatites, mostly generated from lithospheric mantle sources, the parental magmas of the ijolite-clinopyroxenite intrusion also record the influence of sublithospheric (plume-related?) geochemical components.     

Amount of Asian lithospheric mantle subducted during the India/Asia collision

Body wave seismic tomography is a successful technique for mapping lithospheric material sinking into the mantle. Focusing on the India/Asia collision zone, we postulate the existence of several Asian continental slabs, based on seismic global tomography. We observe a lower mantle positive anomaly between 1100 and 900 km depths, that we interpret as the signature of a past subduction process of Asian lithosphere, based on the anomaly position relative to positive anomalies related to Indian continental slab. We propose that this anomaly provides evidence for south dipping subduction of North Tibet lithospheric mantle, occurring along 3000 km parallel to the Southern Asian margin, and beginning soon after the 45 Ma break-off that detached the Tethys oceanic slab from the Indian continent. We estimate the maximum length of the slab related to the anomaly to be 400 km. Adding 200 km of presently Asian subducting slab beneath Central Tibet, the amount of Asian lithospheric mantle absorbed by continental subduction during the collision is at most 600 km. Using global seismic tomography to resolve the geometry of Asian continent at the onset of collision, we estimate that the convergence absorbed by Asia during the indentation process is ~ 1300 km. We conclude that Asian continental subduction could accommodate at most 45% of the Asian convergence. The rest of the convergence could have been accommodated by a combination of extrusion and shallow subduction/underthrusting processes. Continental subduction is therefore a major lithospheric process involved in intraplate tectonics of a supercontinent like Eurasia.