Mapping crustal structure beneath southern Tibet: Seismic evidence for continental crustal underthrusting

Receiver function imaging along a temporary seismic array (ANTILOPE-2) reveals detailed information of the underthrusting of the Indian crust in southern Tibet. The Moho dips northward from ~ 50 km to 80 km beneath the Himalaya terrane, and locally reaches ~ 85 km beneath the Indus–Yalung suture. It remains at ~ 80 km depth across the Lhasa terrane, and shallows to ~ 70 km depth under the Qiangtang terrane. An intra-crustal interface at ~ 60 km beneath the Lhasa terrane can be clearly followed southward through the Main Himalaya Thrust and connects the Main Boundary Thrust at the surface, which represents the border of the Indian crust that is underthrusting until south of the Bangong–Nujiang Suture. A mid-crustal low velocity zone is observed at depths of 14–30 km beneath the Lhasa and Himalaya terranes probably formed by partial melt and/or aqueous fluids.     

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.     

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.     

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.     

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.     

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.     

Subduction of the Indian lower crust beneath southern Tibet revealed by the post-collisional potassic and ultrapotassic rocks in SW Tibet

New major and trace elemental, Sr–Nd–Pb isotope, and zircon U–Pb geochronological and Hf–O isotope data of post-collisional potassic and ultrapotassic volcanic rocks (PVRs and UPVs, respectively) along with geochemical data of PVRs, UPVs, and Mg-rich potassic rocks (MPRs) in the literature are used to constrain their mantle source and genesis. The PVRs, UPVs, and MPRs share similar geochemical features but with some discrepancies, suggesting that they were derived from subcontinental lithospheric mantle (SCLM) with isotopic heterogeneity resulting from the varying contributions of subducted Indian lower crust into the mantle source (ca. 6–20%, ca. 8–30%, and ca. 9–30%, respectively). The zircon Hf–O isotopic compositions of these rocks can be classified into two groups, including Group I rocks with high δ¹⁸O (6.7–11.3‰), low ε_Hf(t) (− 17.0 to − 12.0), and old Hf crustal model ages (1.87–2.19 Ga) that indicate an ancient SCLM source, and Group II rocks with δ¹⁸O values of 6.8–10.7‰, ε_Hf(t) values of − 11.8 to − 6.3, and younger Hf crustal model ages (1.50–1.86 Ga). The negative correlation defined by δ¹⁸O and ε_Hf(t) of Group II samples suggests a two-component mixing between mantle- and crust-derived melts, in which the latter would be the subducted Indian lower crust as indicated by the similar negative ε_Hf(t) values between Group II samples (− 11.8 to − 6.3) and the High Himalayan gneiss (− 14.2 to + 0.3). Thus we propose two enrichment events to account for the Hf–O isotopic compositions of the PVRs and UPVs/MPRs: the first involves the enrichment of the overlying SCLM that was metasomatized by fluids derived from dehydration of the subducted Indian lower crust, and the second invokes the enrichment of the overlying SCLM metasomatized by melts of the already dehydrated different proportions of the Indian lower crust. We argue that break-off of the northwards subducted Indian Plate in the early Miocene caused the asthenospheric upwelling under the Indian plate through slab window, resulting in varying degrees of partial melting of the overlying metasomatized heterogeneous SCLM to produce the primitive magmas of the PVRs, UPVs, and MPRs in an extensional setting. These observations and interpretations imply that the Indian lower crust was subducted beneath the Lhasa terrane in the Early–Middle Miocene.     

3D thermal structure of the continental lithosphere beneath China and adjacent regions

Based on the Crust2.0 model and the topography data of Chinese continent and its adjacent regions, a three-dimensional finite element model is constructed in terms of the spherical coordinate system. In our numerical model, the average annual ground temperature from 195 meteorological stations and temperature of upper mantle derived from the seismic velocities are adopted as the top and bottom boundary conditions, respectively. The observed thermal conductivity and heat production from P wave velocity based on empirical formula are employed in our numerical model as well. The comparison between the calculated and observed surface heat flow proved that our results are reliable. The temperature beneath the Precambrian cratons is lower than that of other areas for 100–300 °C also. The typical temperature rang at the Moho is estimated to be 800–1000 °C beneath the Tibetan plateau and 500–700 °C beneath the Precambrian cratons (such as Indian plate, Sichuan basin, South China, North China and Tarim), respectively. The thermal state in the eastern part of Sino-Korean craton at the depth deeper than 60 km indicates that it was destructed. The thermal structure in center of Tibetan plateau (especially beneath Qiangtang area) supports the proposed flow of lower crustal or upper mantle material to the east. Generally, the distribution of volcanoes in Chinese continent is consistent with the high temperature areas in the crust or upper mantle. There are many obvious thermal transition zones across the orogenic belts. The thermal transition zone between eastern and western parts in the crust of Chinese continent is consistent with the north–south seismic zone.     

西藏及其周围地区地壳、地幔地震层析成像——印度板块大规模俯冲于西藏高原之下的证据

文中介绍有关西藏—喜马拉雅碰撞带的一项地震层析成像研究。根据一个用天然地震数据产生的全球波速模型 ,印度板块有可能以近水平状俯冲于整个西藏高原之下至 16 5～ 2 6 0km深度。西藏岩石圈具有低波速地壳和高波速下岩石圈 (75～ 12 0km深 )。在 12 0～ 16 5km深度范围 ,西藏岩石圈与俯冲的印度板块之间有一层低速软流圈物质。高原中部从地表到 310km深处有一低速体 ,说明地幔物质有可能穿过俯冲板块的脆弱部位上隆。这些结果以及野外实测的地壳缩短值说明高原的抬升得助于印度板块的近水平俯冲。我们推论俯冲印度板块的升温上浮以及上覆软流层的存在是造成西藏高原高海拔抬升以及内部地表仍相对平坦的主要原因。2 0 0 1年 1月 2 6日在印度西部发生的毁灭性大地震有可能是俯冲应力在印度板块后缘薄弱处引发的岩石圈大断裂。     

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.     

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.     

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.     

Upper mantle P-wave tomography across the Longmenshan fault belt from passive-source seismic observations along Aba-Longquanshan profile

A passive seismic experiment across the Longmenshan (LMS) fault belt had been conducted between August 2006 and July 2007 for the understanding of geodynamic process between the Eastern Tibet and Sichuan basin. We herein collected 3677 first P arrival times with high precision from seismograms of 288 teleseismic events so as to reconstruct the upper mantle velocity structure. Our results show that the depth of the Lithosphere–asthenosphere boundary (LAB) changes from 70 km beneath Eastern Tibet to about 110 km beneath Longquanshan, Sichuan Basin, which is consistent with the receiver function imaging results. The very thin mantle part of the lithosphere beneath Eastern Tibet may suggest the lithosphere delamination due to strong interaction between the Tibetan eastward escaping flow and the rigid resisting Sichuan basin, which can be further supported by the existences of two high-velocity anomalies beneath LAB in our imaging result. We also find there are two related low-velocity anomalies beneath the LMS fault belt, which may indicate magmatic upwelling from lithosphere delamination and account for the origin of tremendous energy needed by the devastating Wenchuan earthquake.     

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.     

青藏高原的地幔结构:地幔羽、地幔剪切带及岩石圈俯冲板片的拆沉

通过横穿青藏高原近 80 0 0km长的 4条天然地震层析剖面 ,获得 4 0 0km深度以上的地壳和地幔速度图像及地震波各向异性 ,揭示了青藏高原 4 0 0km深度范围内的地壳和地幔结构特征。地幔速度图像显示 ,青藏高原腹地的深地幔中存在以大型低速异常体为特征的地幔羽 ,其可能通过热通道与大面积分布的可可西里新生代高钾碱性火山作用有成因联系 ;阿尔金、康西瓦、金沙江、嘉黎及雅鲁藏布江等走滑断裂可下延至 30 0～ 4 0 0km深度 ,显示了低速高热物质组成的垂向低速异常带特征及大型超岩石圈或地幔剪切带的产出 ;发现康西瓦、东昆仑—金沙江、班公湖—怒江和雅鲁藏布缝合带下部存在不连续的高速异常带 ,可以解释为青藏高原地体拼合及碰撞过程中可能保留的加里东、古特提斯和中特提斯大洋岩石圈“化石”残片 ,是“拆沉”的地球物理证据。印度大陆岩石圈的巨厚俯冲板片以 15～ 2 0°倾角向北插入唐古拉山下 30 0km深处 ,并被高热物质组成的地幔剪切带分开。结合新的横穿喜马拉雅及青藏高原的地幔层析资料 ,提出青藏高原碰撞动力学新模式 :青藏高原南部印度岩石圈板片的翻卷式陆内超深俯冲 ,北缘克拉通向南的陆内俯冲 ,腹地深部的地幔羽上涌 ,以及地幔范围内的高原“右旋隆升”及物质向东及北东方向运动及挤出。     

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.     

Combined U-Pb,Lu-Hf,Sm-Nd and Ar-Ar multichronometric dating on the Bailang eclogite constrains the closure timing of the Paleo-Tethys Ocean in the Lhasa terrane,Tibet

The Lhasa terrane, the main tectonic component of the Himalayan–Tibetan orogen, has received much attention as it records the entire history of the orogeny. The occurrence of Permian to Triassic high-pressure eclogites has a significant bearing on the understanding of the Paleo-Tethys subduction and plate suturing processes in this area. An eclogite from the Bailang, eastern Lhasa terrane, was investigated with a combined metamorphic P–T and U–Pb, Lu–Hf, Sm–Nd and Ar–Ar multichronometric approach. Pseudosection modeling combined with thermobarometric calculations indicate that the Bailang eclogite equilibrated at peak P–T conditions of ~ 2.6 GPa and 465–503 °C, which is much lower than those of Sumdo and Jilang eclogites in this area. Garnet–whole rock–omphacite Lu–Hf and Sm–Nd ages of 238.1 ± 3.6 Ma and 230.0 ± 4.7 Ma were obtained on the same sample, which are largely consistent with the corresponding U–Pb age of 227.4 ± 6.4 Ma for the metamorphic zircons within uncertainty. The peak metamorphic temperature of the sample is lower than the Lu–Hf and Sm–Nd closure temperatures in garnet. This, combined with the core-to-rim decrease in Mn and HREE concentrations, the slightly U-shaped Sm zonation across garnet and the exclusive occurrence of omphacite inclusion in garnet rim, are consistent with the Lu–Hf system skewing to the age of the garnet core and the Sm–Nd system favoring the rim age. The Sm–Nd age was thus interpreted as the age of eclogite-facies metamorphism and the Lu–Hf age likely pre-dated the eclogite-facies metamorphism. ⁴⁰Ar/³⁹Ar dating of hornblende from the eclogite yielded ages about 200 Ma, which is interpreted as a cooling age and is probably indicative of the time of exhumation to the middle crust. The difference of peak eclogite-facies metamorphic conditions and the distinct metamorphic ages for the Bailang eclogite (~ 2.6 GPa and ~ 480 °C; ca. 230 Ma), the Sumdo eclogite (~ 3.4 GPa and ~ 650 °C; ca. 262 Ma) and Jiang eclogite (~ 3.6 GPa and ~ 750 °C; ca. 261 Ma) in the same (ultra)-high-pressure belt indicate that this region likely comprises different slices that had distinct P–T histories and underwent (U)HP metamorphism at different times. The initiation of the opening the Paleo-Tethys Ocean in the Lhasa terrane could trace back to the early Permian. The ultimate closure of the Paleo-Tethys Ocean in the Lhasa terrane was no earlier than ca. 230 Ma.     

Zedong terrane revisited: An intra-oceanic arc within Neo-Tethys or a part of the Asian active continental margin?

Precise timing of the India-Asia collision is important to constrain the evolution history of both the Himalayan orogen and the Tibetan Plateau. It has been proposed that the Indian plate first collided with an intra-oceanic arc at ∼55 Ma, and then the composite terrane collided with the Asian continent at ∼35 Ma. The Zedong terrane has been suggested to represent the vestige of such an intra-oceanic arc developed within the Neo-Tethys Ocean, as some volcanic rocks with high K₂O have been classified as shoshonites. In this study, we present detailed geochemical and geochronological data of various types of magmatic rocks (including volcanic, cumulate and granitic rocks) widely exposed in the Zedong terrane to constrain the formation age and tectonic setting of the Zedong terrane. We found that the Zedong volcanic rocks belong to calc-alkaline series rather than shoshonites and high K₂O contents in some volcanic rocks resulted from alteration. The basalts are highly enriched in LREE and LILE, but strongly depleted in HFSE, indicating they were derived from a metasomatized mantle. Presence of hornblende phenocryst in both gabbros and hornblendites indicates that the cumulates were produced from hydrous basalts through crystallization. The granitic rocks have adakite-like compositional characteristics, i.e., high Sr/Y ratios but low Y contents, which were formed by melting of a thickened lower crust. Zircons from six samples, including a volcanic rock (an andesite), three cumulates (a hornblendites, a hornblende-bearing gabbro and a gabbro) and two granitic (a tonalite and a granodiorite) rocks, have been dated to yield identical ages of ∼155–160 Ma. This suggests that the volcanic eruption and plutonic emplacement were coevally developed in the Zedong terrane. Zircons from both the andesite and the cumulates have similar positive ε_Hf(t) values (∼+11.6 to +16.7), indicating they were stemmed from similarly depleted mantle sources. Meanwhile, zircons from the granitic rocks also have positive ε_Hf(t) values of ∼+12.6 to +15.2, implying their derivation from a juvenile lower crust. Therefore, we proposed that the basalts in the Zedong terrane were formed through partial melting of the mantle wedge metasomatized by slab-released fluids/melts. A part of hydrous basalts were underplated in the thickened lower crust beneath the Zedong terrane, which gave rise to the cumulate and granitic rocks. By comparison, magmatic rocks in the Zedong terrane show compositional similarities with the Jurassic rocks exposed in the Gangdese arc. This suggests that the Zedong terrane represents a slice of the active continental margin developed on the southern margin of the Lhasa terrane as a result of the northward subduction of the Neo-Tethys Ocean during the Late Jurassic, rather than the vestige of an intra-oceanic arc.     

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.     

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.