Lithospheric structure and crust–mantle decoupling in the southeast edge of the Tibetan Plateau

Convergence between the Indian plate and the Eurasian plate has resulted in the uplift of the Tibetan Plateau, and understanding the associated dynamical processes requires investigation of the structures of the crust and the lithosphere of the Tibetan Plateau. Yunnan is located in the southwest edge of the plateau and adjacent to Myanmar to the west. Previous observations have confirmed that there is a sharp transition in mantle anisotropy in this area, as well as clockwise rotations of the surface velocity, surface strain, and fault orientation. We use S receiver functions from 54 permanent broad-band stations to investigate the structures of the crust and the lithosphere beneath Yunnan. The depth of the Moho is found to range from 36 to 40 km beneath southern Yunnan and from 55 to 60 km beneath northwestern Yunnan, with a dramatic variation across latitude 25–26°N. The depth of the lithosphere–asthenosphere boundary (LAB) ranges from 180 km to less than 70 km, also varying abruptly across latitude 25–26°N, which is consistent with the sudden change of the fast S-wave direction (from NW–SE to E–W across 26–28°N). In the north of the transition belt, the lithosphere is driven by asthenospheric flow from Tibet, and the crust and the upper mantle are mechanically coupled and moving southward. Because the northeastward movement of the crust in the Burma micro-plate is absorbed by the right-lateral Sagaing Fault, the crust in Yunnan keeps the original southward movement. However, in the south of the transition belt, the northeastward mantle flow from Myanmar and the southward mantle flow from Tibet interact and evolve into an eastward flow (by momentum conservation) as shown by the structure of the LAB. This resulting mantle flow has a direction different from that of the crustal movement. It is concluded that the Sagaing Fault causes the west boundary condition of the crust to be different from that of the lithospheric mantle, thus leading to crust–mantle decoupling in Yunnan.     

P-wave tomography and dynamics of the crust and upper mantle beneath western Tibet

The continental collision between the Indian and Asian plates plays a key role in the geologic and tectonic evolution of the Tibetan plateau. In this article we present high-resolution tomographic images of the crust and upper mantle derived from a large number of high-quality seismic data from the ANTILOPE project in western Tibet. Both local and distant earthquakes were used in this study and 35,115 P-wave arrival times were manually picked from the original seismograms. Geological and geochemical results suggested that the subducting Indian plate has reached northward to the Lhasa terrane, whereas our new tomography shows that the Indian plate is currently sub-horizontal and underthrusting to the Jinsha river suture at depths of ~ 100 to ~ 250 km, suggesting that the subduction process has evolved over time. The Asian plate is also imaged clearly from the surface to a depth of ~ 100 km by our tomography, and it is located under the Tarim Basin north of the Altyn Tagh Fault. There is no obvious evidence to show that the Asian plate has subducted beneath western Tibet. The Indian and Asian plates are separated by a prominent low-velocity zone under northern Tibet. We attribute the low-velocity zone to mantle upwelling, which may account for the warm crust and upper mantle beneath that region, and thus explain the different features of magmatism between southern and northern Tibet. But the upwelling may not penetrate through the whole crust. We propose a revised geodynamic model and suggest that the high-velocity zones under Lhasa terrane may reflect a cold crust which has interrupted the crustal flow under the westernmost Tibetan plateau.     

Seismic structure and rheology of the crust under mainland China

The crust and upper mantle in mainland China were relatively densely probed with wide-angle seismic profiling since 1958, and the data have provided constraints on the amalgamation and lithosphere deformation of the continent. Based on the collection and digitization of crustal P-wave velocity models along related wide-angle seismic profiles, we construct several crustal transects across major tectonic units in mainland China. In our study, we analyzed the seismic activity, and seismic energy releases during 1970 and 2010 along them. We present seismogenic layer distribution and calculate the yield stress envelopes of the lithosphere along the transects, yielding a better understanding of the lithosphere rheology strength beneath mainland China. Our results demonstrate that the crustal thicknesses of different tectonic provinces are distinctively different in mainland China. The average crustal thickness is greater than 65 km beneath the Tibetan Plateau, about 35 km beneath South China, and about 36–38 km beneath North China and Northeastern China. For the basins, the thickness is ~ 55 km beneath Qaidam, ~ 50 km beneath Tarim, ~ 40 km beneath Sichuan and ~ 35 km beneath Songliao. Our study also shows that the average seismic P-wave velocity is usually slower than the global average, equivalent with a more felsic composition of crust beneath the four tectonic blocks of mainland China resulting from the complex process of lithospheric evolution during Triassic and Cenozoic continent–continent and Mesozoic ocean–continent collisions. We identify characteristically different patterns of seismic activity distribution in different tectonic blocks, with bi-, or even tri-peak distribution of seismic concentration in South Tibet, which may suggest that crustal architecture and composition exert important control role in lithosphere deformation. The calculated yield stress envelopes of lithosphere in mainland China can be divided into three groups. The results indicate that the lithosphere rheology structure can be described by jelly sandwich model in eastern China, and crème brulee models with weak and strong lower crust corresponding to lithosphere beneath the western China and Kunlun orogenic belts, respectively. The spatial distribution of lithospheric rheology structure may provide important constraints on understanding of intra- or inter-plate deformation mechanism, and more studies are needed to further understand the tectonic process(es) accompanying different lithosphere rheology structures.     

An overview of the crustal structure of the Tibetan plateau after 35 years of deep seismic soundings

Since the pioneer wide-angle seismic profile along the Yadong–Gulu rift acquired in 1974 by the ex-Institute of Geophysics, Chinese Academy of Sciences (CAS), several research programs aimed to deep geophysics, performed thanks to the participation of Chinese national and international institutions, have been developed during last 35 years, including 23 wide-angle seismic profiles with total length of about 6000 km. These profiles are unevenly distributed, most of them in eastern Tibet and few profiles in western Tibet. In this paper, we make a summarized presentation of all these wide-angle seismic profiles and provide an overall view of the seismic velocity structure of the crust beneath the broad Tibetan plateau, which is the product of the continuous convergence and collision of the Indian and Eurasian plates since about 50 Ma ago. Different patterns of crustal thickness variation related to the tectonic blocks and along suture zones of the region are displayed. The crust thickness is confirmed to be about 70–75 km under southern Tibet, and 60–65 km under northern, northeastern and southeastern Tibet. The leading edge of the subducted lithosphere reaches the northern margin of the plateau and directly contacts with Tarim Basin. Westward of the 90°E boundary, the Indian crust is moving towards the northern edge of the plateau and collides with Tarim Basin at 80°E while reach the Bangong–Nujiang suture belt at 88°E; eastward of the 90°E boundary, the northern edge of the crust should be at 50–100 km south of Bangong–Nujiang suture. The results supply helpful constrains to understand the mechanism of the continent–continent collision and its consequences in the plateau and neighbouring areas.     

S-wave velocity structure of the North China from inversion of Rayleigh wave phase velocity

We constructed the S-wave velocity structure of the crust and uppermost mantle (10–100 km) beneath the North China based on the teleseismic data recorded by 187 portable broadband stations deployed in this region. The traditional two-step inversion scheme was adopted. Firstly, we measured the interstation fundamental Rayleigh wave phase velocity of 10–60 s and imaged the phase velocity distributions using the Tarantola inversion method. Secondly, we inverted the 1-D S-wave velocity structure with a grid spacing of 0.25° × 0.25° and constructed the 3-D S-wave velocity structure of the North China. The 3-D S-wave velocity model provides valuable information about the destruction mechanism and geodynamics of the North China Craton (NCC). The S-wave velocity structures in the northwestern and southwestern sides of the North–South Gravity Lineament (NSGL) are obviously different. The southeastern side is high velocity (high-V) while the northeastern side is low velocity (low-V) at the depth of 60–80 km. The upwelling asthenosphere above the stagnated Pacific plate may cause the destruction of the Eastern Block and form the NSGL. A prominent low-V anomaly exists around Datong from 50 to 100 km, which may due to the upwelling asthenosphere originating from the mantle transition zone beneath the Western Block. The upwelling asthenosphere beneath the Datong may also contribute to the destruction of the Eastern Block. The Zhangjiakou-Penglai fault zone (ZPFZ) may cut through the lithosphere and act as a channel of the upwelling asthenosphere. A noticeable low-V zone also exists in the lower crust and upper mantle lid (30–50 km) beneath the Beijing–Tianjin–Tangshan (BTT) region, which may be caused by the upwelling asthenosphere through the ZPFZ.     

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.     

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.     

The Mesoproterozoic mantle plume beneath the northern part of the Siberian craton

The study of the Mesoproterozoic (1473 ± 24 Ma) dolerites of the Olenek uplift of the Siberian craton basement has shown their petrologic and geochemical similarity to typical OIB produced with participation of a mantle plume. The dolerites are characterized by variations in the geochemical composition explained by different degrees of melting of the same source. A conclusion is drawn that the parental melts of the rocks were slightly modified by crustal contamination, as evidenced from their Nd isotope composition (£_Nd(T) = + 0.6 to − 0.8) and the presence of inherited zircons of four ages (2564, 2111, 2053, and 1865 Ma). Since the Siberian craton in the structure of the Nuna supercontinent (Columbia) was located relatively close to the Baltic continent and the Congo and Sao Francisco cratons, we assume that the Early Mesoproterozoic mafic intrusions (1500–1470 Ma) of all these cratons belong to the same large igneous province (LIP). The province formation was related to the activity of superplume (or mantle hot field), which supplied mantle matter to the lithosphere basement. The superplume core was probably located beneath the northern part of the Siberian craton, where basites are compositionally most similar to the primary mantle source.     

Geophysical constraints on geodynamic processes at convergent margins: A global perspective

Convergent margins, being the boundaries between colliding lithospheric plates, form the most disastrous areas in the world due to intensive, strong seismicity and volcanism. We review global geophysical data in order to illustrate the effects of the plate tectonic processes at convergent margins on the crustal and upper mantle structure, seismicity, and geometry of subducting slab. We present global maps of free-air and Bouguer gravity anomalies, heat flow, seismicity, seismic Vs anomalies in the upper mantle, and plate convergence rate, as well as 20 profiles across different convergent margins. A global analysis of these data for three types of convergent margins, formed by ocean–ocean, ocean–continent, and continent–continent collisions, allows us to recognize the following patterns. (1) Plate convergence rate depends on the type of convergent margins and it is significantly larger when, at least, one of the plates is oceanic. However, the oldest oceanic plate in the Pacific ocean has the smallest convergence rate. (2) The presence of an oceanic plate is, in general, required for generation of high-magnitude (M > 8.0) earthquakes and for generating intermediate and deep seismicity along the convergent margins. When oceanic slabs subduct beneath a continent, a gap in the seismogenic zone exists at depths between ca. 250 km and 500 km. Given that the seismogenic zone terminates at ca. 200 km depth in case of continent–continent collision, we propose oceanic origin of subducting slabs beneath the Zagros, the Pamir, and the Vrancea zone. (3) Dip angle of the subducting slab in continent–ocean collision does not correlate neither with the age of subducting oceanic slab, nor with the convergence rate. For ocean–ocean subduction, clear trends are recognized: steeply dipping slabs are characteristic of young subducting plates and of oceanic plates with high convergence rate, with slab rotation towards a near-vertical dip angle at depths below ca. 500 km at very high convergence rate. (4) Local isostasy is not satisfied at the convergent margins as evidenced by strong free air gravity anomalies of positive and negative signs. However, near-isostatic equilibrium may exist in broad zones of distributed deformation such as Tibet. (5) No systematic patterns are recognized in heat flow data due to strong heterogeneity of measured values which are strongly affected by hydrothermal circulation, magmatic activity, crustal faulting, horizontal heat transfer, and also due to low number of heat flow measurements across many margins. (6) Low upper mantle Vs seismic velocities beneath the convergent margins are restricted to the upper 150 km and may be related to mantle wedge melting which is confined to shallow mantle levels.     

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.     

Formation and evolution of Precambrian continental lithosphere in South China

An overview is presented for the formation and evolution of Precambrian continental lithosphere in South China. This is primarily based on an integrated study of zircon U–Pb ages and Lu–Hf isotopes in crustal rocks, with additional constraints from Re–Os isotopes in mantle-derived rocks. Available Re–Os isotope data on xenolith peridotites suggest that the oldest subcontinental lithospheric mantle beneath South China is primarily of Paleoproterozoic age. The zircon U–Pb ages and Lu–Hf isotope studies reveal growth and reworking of the juvenile crust at different ages. Both the Yangtze and Cathaysia terranes contain crustal materials of Archean U–Pb ages. Nevertheless, zircon U–Pb ages exhibit two peaks at 2.9–3.0 Ga and ~ 2.5 Ga in Yangtze but only one peak at ~ 2.5 Ga in Cathaysia. Both massive rocks and crustal remnants (i.e., zircon) of Archean U–Pb ages occur in Yangtze, but only crustal remnants of Archean U–Pb ages occur in Cathaysia. Zircon U–Pb and Lu–Hf isotopes in the Kongling complex of Yangtze suggest the earliest episode of crustal growth in the Paleoarchean and two episodes of crustal reworking at 3.1–3.3 Ga and 2.8–3.0 Ga. Both negative and positive ε_Hf(t) values are associated with Archean U–Pb ages of zircon in South China, indicating both the growth of juvenile crust and the reworking of ancient crust in the Archean. Paleoproterozoic rocks in Yangtze exhibit four groups of U–Pb ages at 2.1 Ga, 1.9–2.0 Ga, ~ 1.85 Ga and ~ 1.7 Ga, respectively. They are associated not only with reworking of the ancient Archean crust in the interior of Yangtze, but also with the growth of the contemporaneous juvenile crust in the periphery of Yangtze. In contrast, Paleoproterozoic rocks in Cathaysia were primarily derived from reworking of Archean crust at 1.8–1.9 Ga. The exposure of Mesoproterozoic rocks are very limited in South China, but zircon Hf model ages suggest the growth of juvenile crust in this period due to island arc magmatism of the Grenvillian oceanic subduction. Magmatic rocks of middle Neoproterozoic U–Pb ages are widespread in South China, exhibiting two peaks at about 830–800 Ma and 780–740 Ma, respectively. Both negative and positive ε_Hf(t) values are associated with the middle Neoproterozoic U–Pb ages of zircon, suggesting not only growth and reworking of the juvenile Mesoproterozoic crust but also reworking of the ancient Archean and Paleoproterozoic crust in the middle Neoproterozoic. The tectonic setting for this period of magmatism would be transformed from arc–continent collision to continental rifting with reference to the plate tectonic regime in South China.     

Shear velocity and Vp/Vs ratio structure of the crust beneath the southern margin of South China continent

We here present the results of the inverse modeling of crustal S-phases recorded from a 400-km-long seismic profile, with azimuth nearly N30W, from Lianxian, near Hunan Province, to Gangkou Island, near Guangzhou City, Guangdong Province, in the southern margin of South China continent. The finding in this case is that many shot gathers provided by this wide-angle seismic experiment show relatively strong reflected and refracted S-phases, in particular some crustal refractions (Sg waves) and Moho reflections (SmS waves or simply Sm waves). The P-wave velocity structure of the crust and uppermost mantle was already obtained through the interpretation of vertical-component shot gathers. Now, with constraints introduced by the P-wave velocity architecture and after picking up S-wave traveltime data on the seismograms, we have obtained the S-velocity model of the crust by adjusting these traveltimes but keeping the geometry of the crustal reflectors. Our results demonstrate: (1) the average crustal S-velocity is about 3.64 km/s to the northwest of the Wuchuan-Sihui fault, and 3.62 km/s to the southeast of this fault; (2) relatively constant S-velocity of about 3.42 km/s for the upper crust, 3.55 km/s for the middle crust and laterally varying shear velocity around 3.82 km/s for the lower crust; (3) correspondingly, Vp/Vs ratio is 1.73 for the upper crust, 1.71 for the middle crust and 1.74 for the lower crust. Both shear velocities and Vp/Vs ratio correlate well with the major active faults that break the study area, and show significant changes especially in the upper crust. High Poisson’s ratio (1.8) is observed at shallow depth beneath the Minzhong depression to the southeast of the Wuchuan-Sihui fault and the Huiyuan depression in the southern margin of South China continent. In contrast, a very low Vp/Vs ratio (1.68) is observed between 8 and 14 km depth beneath Huiyuan. At deeper depth, a high Vp/Vs ratio (1.76) is observed in the lower crust beneath the Minzhong depression.     

Lithospheric architecture and tectonic evolution of the Hudson Bay region

Hudson Bay conceals several fundamental tectonic elements of the North American continent, including most of the ca. 1.9–1.8 Ga Trans-Hudson orogen (THO) and the Paleozoic Hudson Bay basin. Formed due to a collision between two cratons, the THO is similar in scale and tectonic style to the modern Himalayan–Karakorum orogen. During collision, the lobate shape of the indentor (Superior craton) formed an orogenic template that, along with the smaller Sask craton, exerted a persistent influence on the tectonic evolution of the region resulting in anomalous preservation of juvenile Proterozoic crust. Extensive products of 2.72–2.68 Ga and 1.9–1.8 Ga episodes of subduction are preserved, but the spatial scale of corresponding domains increases by roughly an order-of-magnitude (to 1000 km, comparable to modern subduction environments) from the Archean to the Proterozoic. Based on analysis of gravity and magnetic data and published field evidence, we propose a new tectonic model in which Proterozoic crust in the southeastern third of Hudson Bay formed within an oceanic or marginal-basin setting proximal to the Superior craton, whereas the northwestern third is underlain by Archean crust. An intervening central belt truncates the southeastern domains and is interpreted to be a continental magmatic arc.Thick, cold and refractory lithosphere that underlies the Bay is well imaged by surface-wave studies and comprises a large component of the cratonic mantle keel beneath North America. The existence of an unusually thick mantle root indicates that subduction and plate collision during the Trans-Hudson orogeny were ‘root-preserving’ (if not ‘root-forming’) processes. Although the Hudson Bay basin is the largest by surface area of four major intracratonic basins in North America, it is also the shallowest. Available evidence suggests that basin subsidence may have been triggered by eclogitization of lower-crustal material. Compared to other basins of similar age in North America, the relatively stiff lithospheric root may have inhibited subsidence of the Hudson Bay basin.     

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.     

Density structure of the cratonic mantle in southern Africa: 1. Implications for dynamic topography

The origin of high topography in southern Africa is enigmatic. By comparing topography in different cratons, we demonstrate that in southern Africa both the Archean and Proterozoic blocks have surface elevation 500–700 m higher than in any other craton worldwide, except for the Tanzanian Craton. An unusually high topography may be caused by a low density (high depletion) of the cratonic lithospheric mantle and/or by the dynamic support of the mantle with origin below the depth of isostatic compensation (assumed here to be at the lithosphere base). We use free-board constraints to examine the relative contributions of the both factors to surface topography in the cratons of southern Africa. Our analysis takes advantage of the SASE seismic experiment which provided high resolution regional models of the crustal thickness.We calculate the model of density structure of the lithospheric mantle in southern Africa and show that it has an overall agreement with xenolith-based data for lithospheric terranes of different ages. Density of lithospheric mantle has significant short-wavelength variations in all tectonic blocks of southern Africa and has typical SPT values of ca. 3.37–3.41 g/cm³ in the Cape Fold and Namaqua–Natal fold belts, ca. 3.34–3.35 g/cm³ in the Proterozoic Okwa block and the Bushveld Intrusion Complex, ca. 3.34–3.37 g/cm³ in the Limpopo Belt, and ca. 3.32–3.33 g/cm³ in the Kaapvaal and southern Zimbabwe cratons.The results indicate that 0.5–1.0 km of surface topography, with the most likely value of ca. 0.5 km, cannot be explained by the lithosphere structure within the petrologically permitted range of mantle densities and requires the dynamic (or static) contribution from the sublithospheric mantle. Given a low amplitude of regional free air gravity anomalies (ca. + 20 mGal on average), we propose that mantle residual (dynamic) topography may be associated with the low-density region below the depth of isostatic compensation. A possible candidate is the low velocity layer between the lithospheric base and the mantle transition zone, where a temperature anomaly of 100–200 °C in a ca. 100–150 km thick layer may explain the observed reduction in Vs velocity and may produce ca. 0.5–1.0 km to the regional topographic uplift.     

Detrital zircons from Neoproterozoic sedimentary rocks in the Yili Block: Constraints on the affinity of microcontinents in the southern Central Asian Orogenic Belt

The Yili Block is one of the Precambrian microcontinents dispersed in the Central Asian Orogenic Belt (CAOB). Detrital zircon U–Pb ages and Hf isotopic data of Neoproterozoic meta-sedimentary rocks (the Wenquan Group) are presented to constrain the tectonic affinity and early history of the Yili Block. The dating of detrital zircons indicates that both the lower and upper Wenquan Groups have two major populations with ages at 950–880 Ma and 1600–1370 Ma. Moreover, the upper Wenquan Group has two minor populations at ~ 1100 Ma and 1850–1720 Ma. According to the youngest age peaks of meta-sedimentary rocks and the ages of related granitoids, the lower Wenquan Group is considered to have been deposited during the early Neoproterozoic (900–845 Ma), whereas the upper Wenquan Group was deposited at 880–857 Ma. The zircon ε_Hf (t) values suggest that the 1.85–1.72 Ga source rocks for the upper Wenquan Group were dominated by juvenile crustal material, whereas those for the lower Wenquan Group involved more ancient crustal material. For the 1.60–1.37 Ga source rocks, however, juvenile material was a significant input into both the upper and lower Wenquan Groups. Therefore, two synchronous crustal growth and reworking events were identified in the northern Yili Block at ca. 1.8–1.7 Ga and 1.6–1.3 Ga, respectively. After the last growth and reworking event, continuous crustal reworking took place in the northern Yili Block until the early Neoproterozoic. Comparing the age patterns and Hf isotopic compositions of detrital zircons from the Yili Block and the surrounding tectonic units indicates that the Yili Block has a close tectonic affinity to the Chinese Central Tianshan Block in the Precambrian. The Precambrian crustal evolution of the Yili Block is distinct from that of the Siberian, North China and Tarim Cratons. Such difference therefore suggests that the Yili Block and the Chinese Central Tianshan Block may have been united in an isolated Precambrian microcontinent within the CAOB rather than representing two different blocks rifted from old cratons on both sides of the Paleo-Asian Ocean.     

Whole-rock Nd–Hf isotopic study of I-type and peraluminous granitic rocks from the Chinese Altai: Constraints on the nature of the lower crust and tectonic setting

The nature of the lower crust and tectonic setting of the Chinese Altai in the early to middle Paleozoic are still hotly debated. Decoupling between zircon Hf and whole-rock Nd isotopic systems for granites results in different interpretations for the above issues. In order to solve the problem, whole-rock Nd–Hf isotopic analyses were conducted on representative early to middle Paleozoic I-type granite and strongly peraluminous granites and rhyolites from the Chinese Altai. The I-type granites show metaluminous to weakly peraluminous feature and have ε_Nd(t) values ranging from − 2.2 to + 0.8 and ε_Hf(t) from + 3.9 to + 12.9, respectively. The strongly peraluminous granites and rhyolites have similar ε_Nd(t) and ε_Hf(t) values ranging from − 3.0 to + 1.7 and from + 2.1 to + 10.4, respectively. All samples plot above the Terrestrial Array on Nd–Hf isotopic diagram, indicating significant Nd–Hf isotopic decoupling in the magma sources. These samples show flatten HREE pattern and have Lu/Hf ratios similar to the average crust, suggesting that Nd–Hf isotopic decoupling was not originated from an ancient basement with elevated Lu/Hf ratios. The observed isotopic decoupling is similar to those modern island arcs, such as the Lesser Antilles and Sunda, where Nd selectively enriched over Hf due to metasomatism in the mantle wedge and consequently resulted in decoupling between the Sm–Nd and Lu–Hf isotopic systems. Our results, combined with the available data, show that prolonged subduction and crust–mantle interaction caused the Nd–Hf isotopic decoupling in the lithospheric mantle beneath the Chinese Altai. The crust of the Chinese Altai was extracted from the lithospheric mantle and inherited the Nd–Hf isotopic decoupling feature. Therefore, the Hf, rather than Nd, isotopic data more faithfully reflect the nature of the lower crust that was quite juvenile in the Paleozoic, and the Chinese Altai represents an early Paleozoic magmatic arc possibly built near western Mongolia.     

Nature and processes of the lithospheric mantle beneath the western Qinling: Evidence from deformed peridotitic xenoliths in Cenozoic kamafugite from Haoti,Gansu Province,China

The Cenozoic Haoti kamafugite field (23 Ma) is situated at the western Qinling Orogen, Gansu Province in China, which is a conjunction region of the North China Craton, the Yangtze Craton and the Tibetan Plateau. Fresh peridotitic xenoliths entrained in these volcanic rocks provide an opportunity to study the nature and processes of the lithospheric mantle beneath the western Qinling. These xenoliths can be divided into two groups based on the petrological features and mineral compositions, type 1 and type 2. Type 1 xenoliths with strongly deformed texture have higher Fo (90–92.5) contents in olivines, Mg# (91–94) and Cr# (15–35) of clinopyroxenes, and Cr# (36–67) of spinels than the weakly deformed type 2 xenoliths, which have the corresponding values of 89–90, 89–91.5, 10–15 and 5–15 in minerals, respectively. CaO contents in fine-grained olivines are slightly higher than 0.10 wt% compared with coarse-grained ones (less than 0.10 wt%). Fine-grained clinopyroxenes have low Al₂O₃ + CaO contents (generally <23 wt%) relative to coarse-grained ones (>23 wt%). Fo contents in fine-grained olivines mainly in the melt pocket of the type 1 xenoliths are higher than those in coarse-grained ones, which is somewhat contrary to the type 2 xenoliths without melt pocket. Clinopyroxenes of the type 2 display higher Na₂O contents (1.7–1.9 wt%) than those of the type 1 (<1.4 wt%). P–T estimations reveal that the type 1 xenoliths give temperature in range of 1106–1187 °C and pressure of 21–26 kbar and that relatively low temperature (907 and 1022 °C) and pressure (19.0 and 18.5 kbar) for the type 2 xenoliths. The type 1 xenoliths are characterized by depletion due to high degree of partial melting (>10%), modal metasomatic and deformed characteristics, and may represent the old refractory lithospheric mantle. In contrast, the type 2 peridotites show fertile features with low degree of partial melting (<5%) and may represent the newly-accreted lithospheric mantle. The lithospheric mantle beneath the western Qinling underwent partial melting, recrystallization, deformation and metasomatism due to asthenospheric upwelling and the latest decompression responding to the Cenozoic extensive tectonic environment. These processes perhaps are closely related to the evolution of Tibetan Plateau caused by the India-Asian collision.     

Subduction-related metasomatic mantle source in the eastern Central Asian Orogenic Belt: Evidence from amphibolites in the Xilingol Complex,Inner Mongolia,China

The Central Asian Orogenic Belt (CAOB) formed mainly in the Paleozoic due to the closure of the Paleo-Asian oceanic basins and accompanying prolonged accretion of pelagic sediments, oceanic crust, magmatic arcs, and Precambrian terranes. The timing of subduction–accretion processes and closure of the Paleo-Asian Ocean has long been controversial and is addressed in a geochemical and isotopic investigation of mafic rocks, which can yield important insight into the geodynamics of subduction zone environments. The Xilingol Complex, located on the northern subduction–accretion zone of the CAOB, mainly comprises strongly deformed quartzo-feldspathic gneisses with intercalated lenticular or quasi-lamellar amphibolite bodies. An integrated study of the petrology, geochemistry, and geochronology of a suite of amphibolites from the complex constrains the nature of the mantle source and the tectono-metamorphic events in the belt. The protoliths of these amphibolites are gabbros and gabbroic diorites that intruded at ca. 340–321 Ma with positive εHf_(t) values ranging from + 2.89 to + 12.98. Their T_DM1 model ages range from 455 to 855 Ma and peak at 617 Ma, suggesting that these mafic rocks are derived from a depleted continental lithospheric mantle. The primitive magma was generated by variable degrees of partial melting of spinel-bearing peridotites. Fractionation of olivine, clinopyroxene and hornblende has played a dominant role during magma differentiation with little or no crustal contamination. The mafic rocks are derived from a Late Neoproterozoic depleted mantle source that was subsequently enriched by melts affected by slab-derived fluids and sediments, or melts with a sedimentary source rock. The Carboniferous mafic rocks in the northern accretionary zone of the CAOB record a regional extensional event after the Early Paleozoic subduction of the Paleo-Asian Ocean. Both addition of mantle-derived magmas and recycling of oceanic crust played key roles in significant Late Carboniferous (ca. 340–309 Ma) vertical crustal growth in the CAOB. Amphibolite–facies metamorphism (P = 0.34–0.52 GPa, T = 675–708 °C) affected these mafic rocks in the Xilingol Complex at ca. 306–296 Ma, which may be related to the crustal thickening by northward subduction of a forearc oceanic crust beneath the southern margin of the South Mongolian microcontinent. The final formation of the Solonker zone may have lasted until ca. 228 Ma.     

Constraints and deception in the isotopic record; the crustal evolution of the west Musgrave Province,central Australia

The Hf and Nd isotopic evolution of the Musgrave Province, central Australia, is used to constrain the timing of crust formation and lithospheric organisation of Proterozoic Australia. The dataset from this region challenges two widely held tenets of Hf and Nd isotope systematics, namely; that crust formation events can only be identified as periods when crystallisation ages correspond to model ages, and that linear Hf evolution arrays away from depleted mantle (along crustal Lu/Hf or Sm/Nd slopes) reflect reworking of the same source.Hf isotopes in Musgrave Province zircon crystals indicate two major crust formation events at c. 1900 Ma and at 1600–1550 Ma. Although no juvenile rocks or crystals are known from c. 1900 Ma, radiogenic addition into the crust at this time is required to account for consistent Nd and Hf evolution patterns, which show no indication of an initially heterogeneous source. Oxygen isotopes in zircon grains confirm that much of the c. 1900 Ma Hf isotopic signal is not compromised by mixtures. Furthermore, the correspondence between mantle extraction and the commencement of reworking of Archean material supports new crust generation at c. 1900 Ma and a coupling between lower and upper crustal processes. The c. 1900 Ma timing of juvenile addition is dissimilar to that in the Albany–Fraser and Arunta Orogens and may reflect continental arc development on the margin of a southern continent.The general Hf isotopic evolution trend of the Musgrave Province apparently reflects reworking from a dominant c. 1900 Ma source with some additional unradiogenic and radiogenic input through time. However, in the 1220–1050 Ma interval this apparent isotopic evolution contrasts with geological observations that indicate input of voluminous mantle-derived material. Intracontinental rifts and other regions with sustained very-high temperature crustal recycling processes generate magmatic provinces with extreme HFSE-enrichment. This can have a profound influence on isotopic evolution trends, suppressing typical juvenile addition patterns. Isotopic mixture modelling indicates that a significant volume of mantle derived material can be accommodated within HFSE enriched magmas without diverging isotopic signatures from apparent reworking trends. In the Musgrave Province, the crust had become so HFSE enriched during the prolonged Musgrave Orogeny (1220–1150 Ma) that it was insensitive to mantle input, which is estimated to have been as much as 85% during this event.