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.     

Cenozoic volcanism and tectonic evolution of the Tibetan plateau

Cenozoic volcanism on the Tibetan plateau, which shows systematic variations in space and time, is the volcanic response to the India–Asia continental collision. The volcanism gradually changed from Na-rich + K-rich to potassic–ultrapotassic + adakitic compositions along with the India–Asia collision shifting from contact-collision (i.e. “soft collision” or “syn-collision”) to all-sided collision (i.e. “hard collision”). The sodium-rich and potasium-rich lavas with ages of 65–40 Ma distribute mainly in the Lhasa terrane of southern Tibet and subordinately in the Qiangtang terrane of central Tibet. The widespread potassic–ultrapotassic lavas and subordinate adakites were generated from ~ 45 to 26 Ma in the Qiangtang terrane of central Tibet. Subsequent post-collisional volcanism migrated southwards, producing ultrapotassic and adakitic lavas coevally between ~ 26 and 8 Ma in the Lhasa terrane. Then potassic and minor adakitic volcanism was renewed to the north and has become extensive and semicontinuous since ~ 20 Ma in the western Qiangtang and Songpan–Ganze terranes. Such spatial–temporal variations provide important constraints on the geodynamic processes that evolved at depth to form the Tibetan plateau. These processes involve roll-back and break-off of the subducted Neo-Tethyan slab followed by removal of the thickened Lhasa lithospheric root, and consequently northward underthrusting of the Indian lithosphere. The Tibetan plateau is suggested to have risen diachronously from south to north. Whereas the southern part of the plateau may have been created and maintained since the late-Oligocene, the northern plateau would have not attained its present-day elevation and size until the mid-Miocene when the lower part of the western Qiangtang and Songpan–Ganze lithospheres began to founder and detach owing to the persistently northward push of the underthrust Indian lithosphere.     

Detrital zircon U–Pb geochronology,trace-element and Hf isotope geochemistry of the metasedimentary rocks in the Eastern Himalayan syntaxis: Tectonic and paleogeographic implications

The origin of the Greater Himalayan Sequence in the Himalaya and the paleogeographic position of the Lhasa terrane within Gondwanaland remain controversial. In the Eastern Himalayan syntaxis, the basement complexes of the northeastern Indian plate (Namche Barwa Complex) and the South Lhasa terrane (Nyingchi Complex) can be studied to explore these issues. Detrital zircons from the metasedimentary rocks in the Namche Barwa Complex and Nyingchi Complex yield similar U–Pb age spectra, with major age populations of 1.00–1.20 Ga, 1.30–1.45 Ga, 1.50–1.65 Ga and 1.70–1.80 Ga. The maximum depositional ages for their sedimentary protoliths are ~ 1.0 Ga based on the mean ages of the youngest three detrital zircons. Their minimum depositional ages are ~ 477 Ma for the Namche Barwa Complex and ~ 499 Ma for the Nyingchi Complex. Detrital zircons from the Namche Barwa Complex and Nyingchi Complex also display similar trace-element signatures and Hf isotopic composition, indicating that they were derived from common provenance. The trace-element signatures of 1.30–1.45 Ga detrital zircons indicate that the 1.3–1.5 Ga alkalic and mafic rocks belt in the southeastern India is a potential provenance. Most 1.50–1.65 Ga zircons have positive ε_Hf(t) values (+ 1.2 to + 9.0), and most 1.70–1.80 Ga zircons have negative ε_Hf(t) values (− 7.1 to − 1.9), which are compatible with those of the Paleo- to Mesoproterozoic orthogneisses in the Namche Barwa Complex. Provenance analysis indicates that the southern Indian Shield, South Lhasa terrane and probably Eastern Antarctica were the potential detrital sources. Combined with previous studies, our results suggest that: (1) the Namche Barwa Complex is the northeastern extension of the Greater Himalaya Sequence; (2) the metasedimentary rocks in the Namche Barwa Complex represent distal deposits of the northern Indian margin relative to the Lesser Himalaya; (3) the South Lhasa terrane was tectonically linked to northern India before the Cambrian.     

Metamorphism and tectonic evolution of the Lhasa terrane,Central Tibet

The Lhasa terrane in southern Tibet is composed of Precambrian crystalline basement, Paleozoic to Mesozoic sedimentary strata and Paleozoic to Cenozoic magmatic rocks. This terrane has long been accepted as the last crustal block to be accreted with Eurasia prior to its collision with the northward drifting Indian continent in the Cenozoic. Thus, the Lhasa terrane is the key for revealing the origin and evolutionary history of the Himalayan–Tibetan orogen. Although previous models on the tectonic development of the orogen have much evidence from the Lhasa terrane, the metamorphic history of this terrane was rarely considered. This paper provides an overview of the temporal and spatial characteristics of metamorphism in the Lhasa terrane based mostly on the recent results from our group, and evaluates the geodynamic settings and tectonic significance. The Lhasa terrane experienced multistage metamorphism, including the Neoproterozoic and Late Paleozoic HP metamorphism in the oceanic subduction realm, the Early Paleozoic and Early Mesozoic MP metamorphism in the continent–continent collisional zone, the Late Cretaceous HT/MP metamorphism in the mid-oceanic ridge subduction zone, and two stages of Cenozoic MP metamorphism in the thickened crust above the continental subduction zone. These metamorphic and associated magmatic events reveal that the Lhasa terrane experienced a complex tectonic evolution from the Neoproterozoic to Cenozoic. The main conclusions arising from our synthesis are as follows: (1) The Lhasa block consists of the North and South Lhasa terranes, separated by the Paleo-Tethys Ocean and the subsequent Late Paleozoic suture zone. (2) The crystalline basement of the North Lhasa terrane includes Neoproterozoic oceanic crustal rocks, representing probably the remnants of the Mozambique Ocean derived from the break-up of the Rodinia supercontinent. (3) The oceanic crustal basement of North Lhasa witnessed a Late Cryogenian (~ 650 Ma) HP metamorphism and an Early Paleozoic (~ 485 Ma) MP metamorphism in the subduction realm associated with the closure of the Mozambique Ocean and the final amalgamation of Eastern and Western Gondwana, suggesting that the North Lhasa terrane might have been partly derived from the northern segment of the East African Orogen. (4) The northern margin of Indian continent, including the North and South Lhasa, and Qiangtang terranes, experienced Early Paleozoic magmatism, indicating an Andean-type orogeny that resulted from the subduction of the Proto-Tethys Ocean after the final amalgamation of Gondwana. (5) The Lhasa and Qiangtang terranes witnessed Middle Paleozoic (~ 360 Ma) magmatism, suggesting an Andean-type orogeny derived from the subduction of the Paleo-Tethys Ocean. (6) The closure of Paleo-Tethys Ocean between the North and South Lhasa terranes and subsequent terrane collision resulted in the formation of Late Permian (~ 260 Ma) HP metamorphic belt and Triassic (220 Ma) MP metamorphic belt. (7) The South Lhasa terrane experienced Late Cretaceous (~ 90 Ma) Andean-type orogeny, characterized by the regional HT/MP metamorphism and coeval intrusion of the voluminous Gangdese batholith during the northward subduction of the Neo-Tethyan Ocean. (8) During the Early Cenozoic (55–45 Ma), the continent–continent collisional orogeny has led to the thickened crust of the South Lhasa terrane experiencing MP amphibolite-facies metamorphism and syn-collisional magmatism. (9) Following the continuous continent convergence, the South Lhasa terrane also experienced MP metamorphism during Late Eocene (40–30 Ma). (10) During Mesozoic and Cenozoic, two different stages of paired metamorphic belts were formed in the oceanic or continental subduction zones and the middle and lower crust of the hanging wall of the subduction zone. The tectonic imprints from the Lhasa terrane provide excellent examples for understanding metamorphic processes and geodynamics at convergent plate boundaries.     

Crustal shear-wave velocity structure beneath northeast India from teleseismic receiver function analysis

We investigated the seismic shear-wave velocity structure of the crust beneath nine broadband seismological stations of the Shillong–Mikir plateau and its adjoining region using teleseismic P-wave receiver function analysis. The inverted shear wave velocity models show ∼34–38 km thick crust beneath the Shillong Plateau which increases to ∼37–38 km beneath the Brahmaputra valley and ∼46–48 km beneath the Himalayan foredeep region. The gradual increase of crustal thickness from the Shillong Plateau to Himalayan foredeep region is consistent with the underthrusting of Indian Plate beyond the surface collision boundary. A strong azimuthal variation is observed beneath SHL station. The modeling of receiver functions of teleseismic earthquakes arriving the SHL station from NE backazimuth (BAZ) shows a high velocity zone within depth range 2–8 km along with a low velocity zone within ∼8–13 km. In contrast, inversion of receiver functions from SE BAZ shows high velocity zone in the upper crust within depth range ∼10–18 km and low velocity zone within ∼18–36 km. The critical examination of ray piercing points at the depth of Moho shows that the rays from SE BAZ pierce mostly the southeast part of the plateau near Dauki fault zone. This observation suggests the effect of underthrusting Bengal sediments and the underlying oceanic crust in the south of the plateau facilitated by the EW-NE striking Dauki fault dipping 30⁰ toward northwest.     

The longest voyage: Tectonic,magmatic, and paleoclimatic evolution of the Indian plate during its northward flight from Gondwana to Asia

The tectonic evolution of the Indian plate, which started in Late Jurassic about 167 million years ago (~ 167 Ma) with the breakup of Gondwana, presents an exceptional and intricate case history against which a variety of plate tectonic events such as: continental breakup, sea-floor spreading, birth of new oceans, flood basalt volcanism, hotspot tracks, transform faults, subduction, obduction, continental collision, accretion, and mountain building can be investigated. Plate tectonic maps are presented here illustrating the repeated rifting of the Indian plate from surrounding Gondwana continents, its northward migration, and its collision first with the Kohistan–Ladakh Arc at the Indus Suture Zone, and then with Tibet at the Shyok–Tsangpo Suture. The associations between flood basalts and the recurrent separation of the Indian plate from Gondwana are assessed. The breakup of India from Gondwana and the opening of the Indian Ocean is thought to have been caused by plate tectonic forces (i.e., slab pull emanating from the subduction of the Tethyan ocean floor beneath Eurasia) which were localized along zones of weakness caused by mantle plumes (Bouvet, Marion, Kerguelen, and Reunion plumes). The sequential spreading of the Southwest Indian Ridge/Davie Ridge, Southeast Indian Ridge, Central Indian Ridge, Palitana Ridge, and Carlsberg Ridge in the Indian Ocean were responsible for the fragmentation of the Indian plate during the Late Jurassic and Cretaceous times. The Réunion and the Kerguelen plumes left two spectacular hotspot tracks on either side of the Indian plate. With the breakup of Gondwana, India remained isolated as an island continent, but reestablished its biotic links with Africa during the Late Cretaceous during its collision with the Kohistan–Ladakh Arc (~ 85 Ma) along the Indus Suture. Soon after the Deccan eruption, India drifted northward as an island continent by rapid motion carrying Gondwana biota, about 20 cm/year, between 67 Ma to 50 Ma; it slowed down dramatically to 5 cm/year during its collision with Asia in Early Eocene (~ 50 Ma). A northern corridor was established between India and Asia soon after the collision allowing faunal interchange. This is reflected by mixed Gondwana and Eurasian elements in the fossil record preserved in several continental Eocene formations of India. A revised India–Asia collision model suggests that the Indus Suture represents the obduction zone between India and the Kohistan–Ladakh Arc, whereas the Shyok-Suture represents the collision between the Kohistan–Ladakh arc and Tibet. Eventually, the Indus–Tsangpo Zone became the locus of the final India–Asia collision, which probably began in Early Eocene (~ 50 Ma) with the closure of Neotethys Ocean. The post-collisional tectonics for the last 50 million years is best expressed in the evolution of the Himalaya–Tibetan orogen. The great thickness of crust beneath Tibet and Himalaya and a series of north vergent thrust zones in the Himalaya and the south-vergent subduction zones in Tibetan Plateau suggest the progressive convergence between India and Asia of about 2500 km since the time of collision. In the early Eohimalayan phase (~ 50 to 25 Ma) of Himalayan orogeny (Middle Eocene–Late Oligocene), thick sediments on the leading edge of the Indian plate were squeezed, folded, and faulted to form the Tethyan Himalaya. With continuing convergence of India, the architecture of the Himalayan–Tibetan orogen is dominated by deformational structures developed in the Neogene Period during the Neohimalayan phase (~ 21 Ma to present), creating a series of north-vergent thrust belt systems such as the Main Central Thrust, the Main Boundary Thrust, and the Main Frontal Thrust to accommodate crustal shortening. Neogene molassic sediment shed from the rise of the Himalaya was deposited in a nearly continuous foreland trough in the Siwalik Group containing rich vertebrate assemblages. Tomographic imaging of the India–Asia orogen reveals that Indian lithospheric slab has been subducted subhorizontally beneath the entire Tibetan Plateau that has played a key role in the uplift of the Tibetan Plateau. The low-viscosity channel flow in response to topographic loading of Tibet provides a mechanism to explain the Himalayan–Tibetan orogen. From the start of its voyage in Southern Hemisphere, to its final impact with the Asia, the Indian plate has experienced changes in climatic conditions both short-term and long-term. We present a series of paleoclimatic maps illustrating the temperature and precipitation conditions based on estimates of Fast Ocean Atmospheric Model (FOAM), a coupled global climate model. The uplift of the Himalaya–Tibetan Plateau above the snow line created two most important global climate phenomena—the birth of the Asian monsoon and the onset of Pleistocene glaciation. As the mountains rose, and the monsoon rains intensified, increasing erosional sediments from the Himalaya were carried down by the Ganga River in the east and the Indus River in the west, and were deposited in two great deep-sea fans, the Bengal and the Indus. Vertebrate fossils provide additional resolution for the timing of three crucial tectonic events: India–KL Arc collision during the Late Cretaceous, India–Asia collision during the Early Eocene, and the rise of the Himalaya during the Early Miocene.     

Late Triassic crustal growth in southern Tibet: Evidence from the Gangdese magmatic belt

The Gangdese magmatic belt, located in the southern margin of the Lhasa terrane and carrying significant copper and polymetallic mineralization, preserves important information relating to the tectonics associated with Indian–Eurasian collision and the crustal growth of southern Tibet. Here we investigate the Quxu batholith in the central domain of the Gangdese magmatic belt and report the occurrence of hornblende gabbros for the first time. We present petrologic, zircon U–Pb–Hf isotopic and bulk-rock chemistry data on these rocks. The hornblende gabbros display sub-alkaline features, and correspond to tholeiite composition. They also show medium K calc-alkaline to low K affinity. The rocks show enrichment in LILEs and LREEs, but are depleted in HFSEs, indicating a subduction-related active continental margin setting for the magma genesis. Our computations show that the gabbroic pluton was emplaced in the middle-lower crustal depth of ca. 18 km. Zircons from the hornblende gabbros yield crystallization age of ca. 210 Ma, revealing a late Triassic magmatic event. Combined with available data from the Gangdese magmatic belt, our study suggests that the northward subduction of the Neo-Tethys oceanic crust beneath the southern margin of the Lhasa terrane might have been initiated not later than the Norian period of Triassic. Zircons from the hornblende gabbro show positive ε_Hf(t) values of 9.56 to 14.75 (mean value 12.44), corresponding to single stage model ages (T_DM1) in the range of 256 Ma to 459 Ma, attesting to crustal growth in the southern Lhasa terrane associated with the subduction of the Neo-Tethys oceanic crust.     

Post-collisional potassic magmatism in the eastern Lhasa terrane,South Tibet: Products of partial melting of mélanges in a continental subduction channel

Post-collisional, potassic magmatic rocks widely distributed in the eastern Lhasa terrane provide significant information for comprehensive understanding of geodynamic processes of northward subduction of the Indian lithosphere and uplift of the Tibetan Plateau. A combined dataset of whole-rock major and trace elements, Sr–Nd–Pb isotopes, and in situ zircon U–Pb dating and Hf–O isotopic analyses are presented for the Yangying potassic volcanic rocks (YPVR) in the eastern part of the Lhasa terrane, South Tibet. These volcanic rocks consist of trachytes, which are characterized by high K₂O (5.46–9.30 wt.%), SiO₂ (61.34–68.62 wt.%) and Al₂O₃ (15.06–17.36 wt.%), and relatively low MgO (0.47–2.80 wt.%) and FeO_t (1.70–4.90 wt.%). Chondrite-normalized rare earth elements (REE) patterns display clearly negative Eu anomalies. Primitive mantle-normalized incompatible trace elements diagrams exhibit strong enrichment in large ion lithophile elements (LILE) relative to high field strength elements (HFSE) and display significantly negative Nb–Ta–Ti anomalies. Initial isotopic compositions indicate relatively radiogenic Sr [(⁸⁷Sr/⁸⁶Sr)_i = 0.711978–0.712090)] and unradiogenic Nd [(¹⁴³Nd/¹⁴⁴Nd)_i = 0.512121–0.512148]. Combined with their Pb isotopic compositions [(²⁰⁶Pb/²⁰⁴Pb)_i = 18.615–18.774, (²⁰⁷Pb/²⁰⁴Pb)_i = 15.708–15.793, (²⁰⁸Pb/²⁰⁴Pb)_i = 39.274–39.355)], these data are consistent with the involvement of component from subducted continental crustal sediment in their source region. The whole-rock Sr–Nd–Pb isotopic compositions exhibit linear trends between enriched mantle-derived mafic ultrapotassic magmas and relatively depleted crustal contaminants from the Lhasa terrane. The enrichment of the upper mantle below South Tibet is considered to result from the addition of components derived from subducted Indian continental crust to depleted MORB-source mantle during northward underthrusting of the Indian continental lithosphere beneath the Lhasa terrane since India–Asia collision at ~ 55 Ma. Secondary Ion Mass Spectrometry (SIMS) U–Pb zircon analyses yield the eruptive ages of 10.61 ± 0.10 Ma and 10.70 ± 0.18 Ma (weighted mean ages). Zircon Hf isotope compositions [Ɛ_Hf(t) = −4.79 to −0.17], combined with zircon O isotope ratios (5.51–7.22‰), imply an addition of crustal material in their petrogenesis. Clinopyroxene-liquid thermobarometer reveals pressure (2.5–4.1 kbar) and temperature (1029.4–1082.9 °C) of clinopyroxene crystallization, suggesting that depth of the magma chamber was 11.6–16.4 km. Energy-constrained assimilation and fractional crystallization (EC–AFC) model calculation indicates depth of assimilation and fractional crystallization in the region of 14.40–18.75 km underneath the Lhasa terrane, which is in consistent with depth of the magma chamber as suggested by clinopyroxene-liquid thermobarometer. Based on the whole-rock major and trace elements and Sr–Nd–Pb isotope compositions, combined with EC–AFC modeling simulations and zircon Hf–O isotope data, we propose that the YPVR resulted from assimilation and fractional crystallization (AFC) process of the K-rich mafic primitive magmas, which were caused by partial melting of the Indian continental subduction-induced mélange rocks.     

Rapid Eocene erosion,sedimentation and burial in the eastern Himalayan syntaxis and its geodynamic significance

The lower Bomi Group of the eastern Himalayan syntaxis comprises a lithological package of sedimentary and igneous rocks that have been metamorphosed to upper amphibolite-facies conditions. The lower Bomi Group is bounded to the south by the Indus–Yarlung Suture and to the north by unmetamorphosed Paleozoic sediments of the Lhasa terrane. We report U–Pb zircon dating, geochemistry and petrography of gneiss, migmatite, mica schist and marble from the lower Bomi Group and explore their geological implications for the tectonic evolution of the eastern Himalaya. Zircons from the lower Bomi Group are composite. The inherited magmatic zircon cores display ²⁰⁶Pb/²³⁸U ages from ~ 74 Ma to ~ 41.5 Ma, indicating a probable source from the Gangdese magmatic arc. The metamorphic overgrowth zircons yielded ²⁰⁶Pb/²³⁸U ages ranging from ~ 38 Ma to ~ 23 Ma, that overlap the anatexis time (~ 37 Ma) recorded in the leucosome of the migmatites. Our data indicate that the lower Bomi Group do not represent Precambrian basement of the Lhasa terrane. Instead, the lower Bomi Group may represent sedimentary and igneous rocks of the residual forearc basin, similar to the Tsojiangding Group in the Xigaze area, derived from denudation of the hanging wall rocks during the India–Asia continental collision. We propose that following the Indian–Asian collision, the forearc basin was subducted, together with Himalayan lithologies from the Indian continental slab. The minimum age of detrital magmatic zircons from the supracrustal rocks is ~ 41.5 Ma and their metamorphism had happened at ~ 37 Ma. The short time interval (< 5 Ma) suggests that the tectonic processes associated with the eastern Himalayan syntaxis, encompassing uplift and erosion of the Gangdese terrane, followed by deposition, imbrication and subduction of the forearc basin, were extremely rapid during the Late Eocene.     

Exhumation history of the NW Indian Himalaya revealed by fission track and 40Ar/39Ar ages

New fission track and Ar/Ar geochronological data provide time constraints on the exhumation history of the Himalayan nappes in the Mandi (Beas valley) – Tso Morari transect of the NW Indian Himalaya. Results from this and previous studies suggest that the SW-directed North Himalayan nappes were emplaced by detachment from the underthrusted upper Indian crust by 55 Ma and metamorphosed by ca. 48–40 Ma. The nappe stack was subsequently exhumed to shallow upper crustal depths (<10 km) by 40–30 Ma in the Tso Morari dome (northern section of the transect) and by 30–20 Ma close to frontal thrusts in the Baralacha La region. From the Oligocene to the present, exhumation continued slowly.Metamorphism started in the High Himalayan nappe prior to the Late Oligocene.High temperatures and anatexis of the subducting upper Indian crust engendered the buoyancy-driven ductile detachment and extrusion of the High Himalayan nappe in the zone of continental collision. Late extrusion of the High Himalayan nappe started about 26 Ma ago, accompanied by ductile extensional shearing in the Zanskar shear zone in its roof between 22 and 19 Ma concomitant with thrusting along the basal Main Central Thrust to the south. The northern part of the nappe was then rapidly exhumed to shallow depth (<10 km) between 20 and 6 Ma, while its southern front reached this depth at 10–5 Ma.     

喜马拉雅西部雅鲁藏布江缝合带地壳尺度的构造叠置

陆陆碰撞过程是板块构造缺失的链条。印度板块与亚洲板块的碰撞造就了喜马拉雅造山带和青藏高原的主体。然而,人们对印度板块在大陆碰撞过程中的行为尚不了解。如大陆碰撞及其碰撞后的大陆俯冲是如何进行的、印度板块是俯冲在青藏高原之下还是回转至板块上部(喜马拉雅造山带内)以及两者比例如何,这些仍是亟待解决的问题。印度板块低角度沿喜马拉雅主逆冲断裂(MHT)俯冲在低喜马拉雅和高喜马拉雅之下已经被反射地震图像很好地揭示。然而,关于MHT如何向北延伸,前人的研究仅获得了分辨率较低的接收函数图像。因而,MHT和雅鲁藏布江缝合带之间印度板块的俯冲行为仍是一个谜。喜马拉雅造山楔增生机制,也就是印度地壳前缘的变形机制,反映出物质被临界锥形逆冲断层作用转移到板块上部,或是以韧性管道流的样式向南溢出。在本次研究中,我们给出在喜马拉雅造山带西部地区横过雅鲁藏布江缝合带的沿东经81.5°展布的高分辨率深地震反射剖面,精细揭示了地壳尺度结构构造。剖面显示,MHT以大约20°的倾斜角度延伸至大约60 km深度,接近埋深为70~75 km的Moho面。越过雅鲁藏布江缝合带运移到北面的印度地壳厚度已经不足15 km。深地震反射剖面还显示中地壳逆冲构造反射发育。我们认为,伴随着印度板块俯冲,地壳尺度的多重构造叠置作用使物质自MHT下部的板块向其上部板块转移,这一过程使印度地壳厚度减薄了,同时加厚了喜马拉雅地壳。     

Isotopic and structural constraints on the late Miocene to Pliocene evolution of the Namche Barwa area,eastern Himalayan syntaxis,SE Tibet

Within the Namche Barwa area, SE Tibet, the Indus–Yarlung suture zone separates the Lhasa terrain in the north from the Himalayan unit including the Tethyan (sedimentary and volcanic rocks), Dongjiu (greenschist to lower amphibolite facies), Namche Barwa (granulite facies), Pei (amphibolite facies) and Laiguo (greenschist facies) sequences in the south. Two fault systems were distinguished in the Namche Barwa area. The former includes a top-down-to-the-north normal fault in the north and two top-to-the-south thrust zones in the south named as Upper and Lower Thrusts, respectively. The Namche Barwa and Pei sequences were exhumed southwards from beneath the Dongjiu sequence by these faults. Thus, the fault system is regarded as a southward extrusion structure. Subsequently, the exposed Dongjiu, Namche Barwa, Pei and Laiguo sequences were displaced northwards onto the Lhasa terrain by the top-to-the-north fault system, thus, marking it as northward indentation structure. Monazite TIMS U–Pb dating demonstrates that the normal fault and the Lower Thrust from the southward extrusion system were probably active at ~ 6 Ma and ~ 10 Ma, respectively. Zircon U–Pb SHRIMP and phlogopite K–Ar ages further suggest that the Upper Thrust was active between 6.2 ± 0.2 Ma and 5.5 ± 0.2 Ma. The northward indentation structures within the core portion of the eastern Himalayan syntaxis were perhaps active between 3.0 Ma and 1.5 Ma, as inferred by published zircon U–Pb SHRIMP and hornblende Ar–Ar ages. The monazite from upper portions of the Pei sequence dated by U–Pb TIMS indicates that the precursor sediments of this sequence were derived from Proterozoic source regions. Nd isotopic data further suggest that all the metamorphic rocks within eastern Himalaya (εNd = ? 13 to ? 19) correlate closely with those from the Greater Himalayan Sequences, whereas the western Himalayan syntaxis is mainly comprised of Lesser Himalayan Sequences. The two indented corners of the Himalaya are, thus, different.     

Crustal density structure from 3D gravity modeling beneath Himalaya and Lhasa blocks,Tibet

The Himalaya and Lhasa blocks act as the main belt of convergence and collision between the Indian and Eurasian plates. Their crustal structures can be used to understand the dynamic process of continent–continent collision. Herein, we present a 3D crustal density model beneath these two tectonic blocks constrained by a review of all available active seismic and passive seismological results on the velocity structure of crust and lower lithosphere. From our final crustal density model, we infer that the present subduction-angle of the Indian plate is small, but presents some variations along the west–east extension of the orogenic belt: The dip angle of the Moho interface is about 8–9° in the eastern and western part of the orogenic belt, and about 16° in the central part. Integrating crustal P-wave velocity distribution from wide-angle seismic profiling, geothermal data and our crustal density model, we infer a crustal composition model, which is composed of an upper crust with granite–granodiorite and granite gneiss beneath the Lhasa block; biotite gneiss and phyllite beneath the Himalaya, a middle crust with granulite facies and possible pelitic gneisses, and a lower crust with gabbro–norite–troctolite and mafic granulite beneath the Lhasa block. Our density structure (<3.2 g/cm³) and composition (no fitting to eclogite) in the lower crust do not be favor to the speculation of ecologitized lower crust beneath Himalaya and the southern of Lhasa block.     

Petrogenesis of the Guangtoushan granitoid suite,central China: Implications for Early Mesozoic geodynamic evolution of the Qinling Orogenic Belt

The Qinling Orogenic Belt marks the link between the South China and North China Blocks and is an important region to understand the geological evolution of the Chinese mainland as well as the Asian tectonic collage. However, the tectonic affinity and geodynamic evolution of the South Qinling Tectonic Belt (SQTB), a main unit of the Qinling Orogenic Belt, remains debated. Here we present detailed geological, geochemical and zircon U–Pb–Hf isotopic studies on the Zhangjiaba, Xinyuan, Jiangjiaping, Guangtoushan and Huoshaodian plutons from the Guangtoushan granitoid suite (GGS) in the western segment of the SQTB. Combining geology, geochronology and whole-rock geochemistry, we identify four distinct episodes of magmatism as: (1) ~ 230–228 Ma quartz diorites and granodiorites, (2) ~ 224 Ma fine-grained granodiorites and monzogranites, (3) ~ 218 Ma porphyritic monzogranites and (4) ~ 215 Ma high-Mg# quartz diorites and granodiorites as well as coeval muscovite monzogranites. The ~ 230–228 Ma quartz diorites and granodiorites were generated by magma mixing between a mafic melt from mantle source and a granodioritic melt derived from partial melting of Neoproterozoic rocks in the lower continental crust related to a continental arc regime. The ~ 224 Ma fine-grained granodiorites and monzogranites were formed through partial melting of a transitional source with interlayers of basaltic rocks and greywackes in the deep zones of the continental arc. The ~ 218 Ma porphyritic monzogranites originated from partial melting of metamorphosed greywackes in lower crustal levels, suggesting underthrusting of middle or upper crustal materials into lower crustal depths. The ~ 215 Ma high-Mg# quartz diorites and granodiorites (with Mg# values higher than 60) were derived from an enriched mantle altered by sediment-derived melts. Injection of hot mantle-derived magmas led to the emergence of the ~ 215 Ma S-type granites at the final stage.Integrating our studies with previous data, we propose that the Mianlue oceanic crust was still subducting beneath the SQTB during ~ 248–224 Ma, and final closure of the Mianlue oceanic basin occurred between ~ 223 Ma and ~ 218 Ma. After continental collision between the South China Block and the SQTB, slab break-off occurred, following which the SQTB transformed into post-collisional extension setting.     

The exhumation history of collision-related mineralizing systems in Tibet: Insights from thermal studies of the Sharang and Yaguila deposits,central Lhasa

The large, newly discovered Sharang porphyry Mo deposit and nearby Yaguila skarn Pb–Zn–Ag (–Mo) deposit reside in the central Lhasa terrane, northern Gangdese metallogenic belt, Tibet. Multiple mineral chronometers (zircon U–Pb, sericite ⁴⁰Ar–³⁹Ar, and zircon and apatite (U–Th)/He) reveal that ore-forming porphyritic intrusions experienced rapid cooling (> 100 °C/Ma) during a monotonic magmatic–hydrothermal evolution. The magmatic–hydrothermal ore-forming event at Sharang lasted ~ 6.0 Myr (~ 1.8 Myr for cooling from > 900 to 350 °C and ~ 4.0 Myr for cooling from 350 to 200 °C) whereas cooling was more prolonged during ore formation at Yaguila (~ 1.8 Myr from > 900 to 500 °C and a maximum of ~ 16 Myr from > 900 to 350 °C). All porphyritic intrusions in the ore district experienced exhumation at a rate of 0.07–0.09 mm/yr (apatite He ages between ~ 37 and 30 Ma). Combined with previous studies, this work implies that uplift of the eastern section of the Lhasa terrane expanded from central Lhasa (37–30 Ma) to southern Lhasa (15–12 Ma) at an increasing exhumation rate. All available geochronologic data reveal that magmatic–hydrothermal–exhumation activities in the Sharang–Yaguila ore district occurred within four periods of magmatism with related mineralization. Significant porphyry-type Mo mineralization was associated with Late Cretaceous–Eocene felsic porphyritic intrusions in the central Lhasa terrane, resulting from Neotethyan oceanic subduction and India–Asia continental collision.     

Present-day crustal motion around the Pamir Plateau from GPS measurements

We have used geodetic techniques to improve constraints on the crustal motion of the Pamir Plateau. Three campaigns of Global Position System data acquisition between 2011 and 2015 demonstrate that, in association with the India–Asia collision, a complex pattern of crustal motion exists in the Pamir Plateau. In a north–south direction from the Indian Plate to the Hazak Block, the crust has absorbed ~ 35 mm/yr of shortening, of which ~ 35% is distributed around the Hindu Kush region (~ 12 mm/yr), and another ~ 35% is taken up around the Alai Valley (also ~ 12 mm/yr). Global Position System measurements also show ~ 5 mm/yr of shortening between the Pamir Plateau and the Tajik Basin, whereas between the Pamir and the Tarim Basin, an ~ 10 mm/yr extension rate is observed. With respect to the stable Eurasian Plate, the Pamir rotates counterclockwise at a rate of ~ 1.822°Myr^− 1, with an Euler pole positioned about the west end of the Tajik Basin (37.03 ± 0.74°N, 65.89 ± 0.12°E). The strain rate field calculated from Global Position System velocities reveals that the crustal motion is consistent with localized deformation around the Hindu Kush and the Alai Valley, the latter representing a zone with strong shallow seismic activity.     

The youngest eclogite in central Himalaya: P–T path,U–Pb zircon age and its tectonic implication

Relict omphacite inclusions have been discovered in mafic granulite at Dinggye of China, confirming the existence of eclogite in central Himalayan orogenic belt. Detailed petrological studies show that relict omphacite occur as inclusions in both garnets and zircons, and the peak mineral assemblage of eclogite-facies should be garnet, omphacite, rutile, muscovite and quartz which was strongly overprinted by granulite-facies minerals during the exhumation. Phase equilibria modeling and associated geothermometer predict that the minimum P–T conditions for peak eclogite-facies stage are 720–760 °C and 20–21 kbar, and those of overprinted granulite-facies are 750 °C and 7–9 kbar in water-undersaturated condition. Thus, a near isothermal decompression P–T path for central Himalayan eclogite has been obtained. Zircon SHRIMP U–Pb dating of two studied eclogite samples at Dinggye yields the peak metamorphic ages of 13.9 ± 1.2 Ma and 14.9 ± 0.7 Ma, respectively, which indicates that the Dinggye eclogite should be the youngest eclogite in Himalayan orogenic belt. Geochemical characteristics and zircon analyses show that the protoliths of eclogite in Dinggye are predicted to be continental rift-related basaltic rocks. The eclogite at Dinggye in central Himalaya should be formed by the crustal thickening during the long-lasting continental overthrusting by Indian plate beneath Euro-Asian continent, and its exhumation process may be related with channel flow and orogen-parallel extension. In the middle Miocene (~ 14 Ma), Indian continental crust had reached at least ~ 65 km depth in southern Tibet.     

A new kinematic evolutionary model for the growth of a duplex — an example from the Rangit duplex,Sikkim Himalaya,India

The Lesser Himalayan duplex (LHD) is a prominent structure through much of the Lesser Himalayan fold–thrust belt. In the Darjeeling - Sikkim Himalaya a component of the LHD is exposed in the Rangit window as the Rangit duplex (RD). The RD consists of ten horses of the upper Lesser Himalayan Sequence (Gondwana, Buxa, Upper Daling). The duplex varies from hinterland-dipping in the north, through an antiformal stack in the middle to foreland-dipping in the south. The Ramgarh thrust (RT) is the roof thrust and, based on a balanced cross-section, the Main Himalayan Sole thrust is the floor thrust at a depth of ~ 10 km and with a dip of ~ 3.5° N.Retrodeformation suggests that the RD initiated as a foreland-dipping duplex with the Early Ramgarh thrust as the roof thrust and the RT as the floor thrust. The RT became the roof thrust during continued duplexing by a combination of footwall imbrication and concurrent RT reactivation. This kinematic history best explains the large translation of the overlying MCT sheets. The restoration suggests that RD shortening is ~ 125 km, and the original Gondwana basin extended ~ 142 km northward of its present northernmost exposures within the window.     

Detrital zircon U–Pb age and Hf isotope constrains on the generation and reworking of Precambrian continental crust in the Cathaysia Block,South China: A synthesis

The South China Block, consisting of the Yangtze and the Cathaysia blocks, is one of the largest Precambrian blocks in eastern Asia. However, the early history of the Cathaysia Block is poorly understood due largely to intensive and extensive reworking by Phanerozoic polyphase orogenesis and magmatism which strongly overprinted and obscured much of the Precambrian geological record. In this paper, we use the detrital zircon U–Pb age and Hf isotope datasets as an alternative approach to delineate the early history of the Cathaysia Block. Compilation of published 4041 Precambrian detrital zircon ages from a number of (meta)sedimentary samples and river sands exhibits a broad age spectrum, with three major peaks at ~ 2485 Ma, ~ 1853 Ma and ~ 970 Ma (counting for ~ 10%, ~ 16% and ~ 24% of all analyses, respectively), and four subordinate peaks at ~ 1426 Ma, ~ 1074 Ma, ~ 780 Ma and ~ 588 Ma. Five of seven detrital zircon age peaks are broadly coincident with the crystallisation ages of ~ 1.89–1.83 Ga, ~ 1.43 Ga, ~ 1.0–0.98 Ga and ~ 0.82–0.72 Ga for known igneous rocks exposed in Cathaysia, whereas, igneous rocks with ages of ~ 2.49 Ga and ~ 0.59 Ga have not yet been found. The Hf isotopic data from 1085 detrital zircons yield Hf model ages (T_DM^C) between ~ 4.19 Ga and ~ 0.81 Ga, and the calculated εHf(t) values between − 40.2 and 14.4. The Archean detrital zircons are exclusively oval in shape with complicated internal textures, indicating that they were sourced by long distance transportations and strong abrasion from an exotic Archean continent. In contrast, the majority of detrital zircons in age between ~ 1.9 and ~ 0.8 Ga are euhedral to subhedral crystals, indicative of local derivation by short distance transportations from their sources. The oldest crustal basement rocks in Cathaysia were most likely formed by generation of juvenile crust and reworking of recycled Archean components in Late Paleoproterozoic at ~ 1.9–1.8 Ga, rather than in the Archean as previously speculated. Reworking and recycling of the continental crust are likely the dominant processes for the crustal evolution of Cathaysia during the Mesoproterozoic to Neoproterozoic time, with an intervenient period of significant generation of juvenile crust at ~ 1.0 Ga.Precambrian crustal evolutions of the Cathaysia Block are genetically related to the supercontinent cycles. By comparing detrital zircon data from Cathaysia with those for other continents, and integrating multiple lines of geological evidence, we interpret the Cathaysia Block as an orogenic belt located between East Antarctica, Laurentia and Australia during the assembly of supercontinent Columbia/Nuna at ~ 1.9–1.8 Ga. The Cathaysia Block amalgamated with the Yangtze Block to form the united South China Block during the Sibao Orogeny at ~ 1.0–0.89 Ga. The Laurentia–Cathaysia–Yangtze–Australia–East Antarctica connection gives the best solution to the paleo-position of Cathaysia in supercontinent Rodinia. The significant amount of ~ 0.6–0.55 Ga detrital zircons in Cathaysia and West Yangtze have exclusively high crustal incubation time of > 300 Ma, indicating crystallisation from magmas generated dominantly by crustal reworking. This detrital zircon population compares well with the similar-aged zircon populations from a number of Gondwana-derived terranes including Tethyan Himalaya, High Himalaya, Qiangtang and Indochina. The united South China–Indochina continent was likely once an integral part of Gondwanaland, connected to northern India by a “Pan-African” collisional orogen.     

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.