Seismic imaging of the downwelling Indian lithosphere beneath central Tibet

Correlated crustal and mantle melting documents proto-Tibetan Plateau growth

The mechanism that causes the rapid uplift and active magmatism of the Hoh-Xil Basin in the northern Tibetan Plateau and hence the outward growth of the proto-plateau is highly debated, more specifically, over the relationship between deep dynamics and surface uplift. Until recently the Hoh-Xil Basin remained uncovered by seismic networks due to inaccessibility. Here, based on linear seismic arrays across the Hoh-Xil Basin, we present a three-dimensional S-wave velocity (VS) model of the crust and uppermost mantle structure beneath the Tibetan Plateau from ambient noise tomography. This model exhibits a widespread partially molten crust in the northern Tibetan Plateau but only isolated pockets in the south manifested as low-VS anomalies in the middle crust. The spatial correlation of the widespread low-VS anomalies with strong uppermost mantle low-VS anomalies and young exposed magmatic rocks in the Hoh-Xil Basin suggests that the plateau grew through lithospheric mantle removal and its driven magmatism.

Evidence of low velocity layers at the top and bottom of the Mantle Transition Zone (MTZ) below the Uttarakhand Himalaya, India

The Mantle Transition Zone (MTZ) beneath the Uttarakhand Himalaya has been modelled using Common Conversion Point (CCP) stacking and depth-migration of radial P-receiver functions. In the Uttarakhand Himalaya region, the depths of the 410-km discontinuity (d410) and the 660-km discontinuity (d660) are estimated to be approximately 406 ± 8 km and 659 ± 10 km, respectively. Additionally, the thickness of the mantle transition zone (MTZ) is modelled to be 255 ± 7 km. The average arrival times for d410 and d660 conversions are (44.47 ± 1.33) s and (71.08 ± 1.29) s, respectively, indicating an undisturbed slightly deeper d410 and a deformed noticeably deeper d660 in the area. The model identifies the characteristics of the d410 and d660 mantle discontinuities beneath the Lesser Himalayan region, revealing a thickening of the MTZ towards northeast, which could be due to gradual cooling or thickening of the Indian lithosphere towards its northern limit. We simulate a low-velocity layer (perhaps partially molten) above the d410 discontinuity at depths of 350 to 385 km, indicating the presence of a hydrated MTZ beneath the area. We also interpret a negative phase at d660 as a low-velocity layer between 590 and 640 km depths, which could be attributed to the accumulation of old subducted oceanic materials or increased water content at the bottom of the MTZ. Our results suggest the presence of residues from paleo-subducted lithospheric slabs in and below the mantle transition zone underlying the Uttarakhand Himalayas.

Fragmentation of continental subduction is ending the Himalayan orogeny

After two continents collide, plate convergence and orogenesis are sustained because subducted continental lithosphere continues pulling the surface plate. It remains controversial how, why, and when continental plate convergence and collision slow down and eventually cease. We use an unprecedented data coverage and present a regional-scale seismic tomographic image of the mantle structure beneath the Tibetan Plateau. In the mantle transition zone, we identify multiple high-velocity anomalies and interpret them as detached pieces of the Indian continental slab. Facilitated by internal heterogeneity of the continental lithosphere, piecewise slab detachments could reduce the slab pull force, resulting in the Miocene slowdown of the India-Eurasia convergence and coeval diachronous potassic volcanism in southern Tibet. We propose that slab detachment is a mechanism that eventually will lead to the end of the Indo-Eurasian continental collision and the Himalayan orogeny.

Evidence of structural segmentation of the Uttarakhand Himalaya and its implications for earthquake hazard

The earthquake hazard associated with the Main Himalayan Thrust (MHT) is a critical issue for India and its neighbouring countries in the north. We used data from a dense seismic network in Uttarakhand, India, to model the lateral variations in the depths of MHT (2-6% drop in V_s at 12-21 km depths), Moho (a sharp increase in V_s (by ~ 0.5-0.7 km/s) at 39-50 km depths) and lithosphere (a marked decrease in V_s(~ 1-3%) at 136-178 km depths), across the Himalayan collisional front. Our joint inversion of radial PRFs and group velocity dispersion data of Rayleigh waves detects three NNE trending transverse lithospheric blocks segmenting the lithosphere in Uttarakhand Himalaya, which spatially correlate well with the northward extension of the Delhi -Haridwar Indian basement ridge, an inferred tectonic boundary and great boundary fault, respectively. Our radial receiver function imaging detects highly deformed and segmented crustal and lithospheric structures associated with three mapped transverse lithospheric blocks, suggesting a reduction in rupture lengths of future earthquakes, thereby, reducing earthquake hazards in Uttarakhand.

No mafic layer in 80 km thick Tibetan crust

All models of the magmatic and plate tectonic processes that create continental crust predict the presence of a mafic lower crust. Earlier proposed crustal doubling in Tibet and the Himalayas by underthrusting of the Indian plate requires the presence of a mafic layer with high seismic P-wave velocity (V_p > 7.0 km/s) above the Moho. Our new seismic data demonstrates that some of the thickest crust on Earth in the middle Lhasa Terrane has exceptionally low velocity (V_p < 6.7 km/s) throughout the whole 80 km thick crust. Observed deep crustal earthquakes throughout the crustal column and thick lithosphere from seismic tomography imply low temperature crust. Therefore, the whole crust must consist of felsic rocks as any mafic layer would have high velocity unless the temperature of the crust were high. Our results form basis for alternative models for the formation of extremely thick juvenile crust with predominantly felsic composition in continental collision zones.

Combined Gravimetric-Seismic Moho Model of Tibet

Substantial progress has been achieved over the last four decades to better understand a deep structure in the Himalayas and Tibet. Nevertheless, the remoteness of this part of the world still considerably limits the use of seismic data. A possible way to overcome this practical restriction partially is to use products from the Earth’s satellite observation systems. Global topographic data are provided by the Shuttle Radar Topography Mission (SRTM). Global gravitational models have been derived from observables delivered by the gravity-dedicated satellite missions, such as the Gravity Recovery and Climate Experiment (GRACE) and the Gravity field and steady-state Ocean Circulation Explorer (GOCE). Optimally, the topographic and gravity data should be combined with available results from tomographic surveys to interpret the lithospheric structure, including also a Moho relief. In this study, we use seismic, gravity, and topographic data to estimate the Moho depth under orogenic structures of the Himalayas and Tibet. The combined Moho model is computed based on solving the Vening Meinesz–Moritz (VMM) inverse problem of isostasy, while incorporating seismic data to constrain the gravimetric solution. The result of the combined gravimetric-seismic data analysis exhibits an anticipated more detailed structure of the Moho geometry when compared to the solution obtained merely from seismic data. This is especially evident over regions with sparse seismic data coverage. The newly-determined combined Moho model of Tibet shows a typical contrast between a thick crustal structure of orogenic formations compared to a thinner crust of continental basins. The Moho depth under most of the Himalayas and the Tibetan Plateau is typically within 60–70 km. The maximum Moho deepening of ~76 km occurs to the south of the Bangong-Nujiang suture under the Lhasa terrane. Local maxima of the Moho depth to ~74 km are also found beneath Taksha at the Karakoram fault. This Moho pattern generally agrees with the findings from existing gravimetric and seismic studies, but some inconsistencies are also identified and discussed in this study.

Multi-terrane structure controls the contrasting lithospheric evolution beneath the western and central-eastern Tibetan plateau

The Tibetan plateau is manifested by contrasting along-strike lithospheric structures, but its formation mechanism and the relationship with the heterogeneous multi-terrane configuration is a challenging problem. Here we conduct systematic numerical modeling to explore the roles of width, density, and rheological properties of the multiple terranes in the lithospheric evolution of the Tibetan plateau, which reveals two distinct collision modes. In Mode-I, the lithospheric mantles of both the strong and weak terranes in the Tibetan plate are completely detached, followed by the underthrusting of Indian lithosphere beneath the whole plateau. Alternatively, Mode-II is characterized by full detachment of the weak terranes, but (partial) residue of the strong terranes during collision. These two contrasting modes, broadly consistent with the lithospheric structures of western and central–eastern Tibetan plateau, respectively, are strongly dependent on the along-strike variation of the width of the strong Lhasa–Qiangtang terranes.

The Tibetan plateau is manifested by contrasting along-strike lithospheric structures, but the correlation with multi-terrane configuration remains challenging. Here, the authors show the crucial roles of the original geometric shape of accreted terranes in regulating the lithospheric evolution of Tibetan plateau.

Tearing of Indian mantle lithosphere from high-resolution seismic images and its implications for lithosphere coupling in southern Tibet

What happened to the Indian mantle lithosphere (IML) during the Indian-Eurasian collision and what role it has played on the plateau growth are fundamental questions that remain unanswered. Here, we show clear images of the IML from high-resolution P and S tomography, which suggest that the subducted IML is torn into at least four pieces with different angles and northern limits, shallower and extending further in the west and east sides while steeper in the middle. Intermediate-depth earthquakes in the lower crust and mantle are located almost exclusively in the high-velocity (and presumably strong) part of the Indian lithosphere. The tearing of the IML provides a unified mechanism for Late Miocene and Quaternary rifting, current crustal deformation, and intermediate-depth earthquakes in the southern and central Tibetan Plateau and suggests that the deformations of the crust and the mantle lithosphere are strongly coupled.

Lithospheric foundering and underthrusting imaged beneath Tibet

Long-standing debates exist over the timing and mechanism of uplift of the Tibetan Plateau and, more specifically, over the connection between lithospheric evolution and surface expressions of plateau uplift and volcanism. Here we show a T-shaped high wave speed structure in our new tomographic model beneath South-Central Tibet, interpreted as an upper-mantle remnant from earlier lithospheric foundering. Its spatial correlation with ultrapotassic and adakitic magmatism supports the hypothesis of convective removal of thickened Tibetan lithosphere causing major uplift of Southern Tibet during the Oligocene. Lithospheric foundering induces an asthenospheric drag force, which drives continued underthrusting of the Indian continental lithosphere and shortening and thickening of the Northern Tibetan lithosphere. Surface uplift of Northern Tibet is subject to more recent asthenospheric upwelling and thermal erosion of thickened lithosphere, which is spatially consistent with recent potassic volcanism and an imaged narrow low wave speed zone in the uppermost mantle.

The timing and mechanism of uplift of the Tibetan plateau continues to be a source of debate. Here, the authors present a new tomographic model revealing a T-shaped high wave speed structure beneath South-Central Tibet and interpret this an upper-mantle remnant from lithospheric foundering.

Pliocene-Quaternary crustal melting in central and northern Tibet and insights into crustal flow

There is considerable controversy over the nature of geophysically recognized low-velocity–high-conductivity zones (LV–HCZs) within the Tibetan crust, and their role in models for the development of the Tibetan Plateau. Here we report petrological and geochemical data on magmas erupted 4.7–0.3 Myr ago in central and northern Tibet, demonstrating that they were generated by partial melting of crustal rocks at temperatures of 700–1,050 °C and pressures of 0.5–1.5 GPa. Thus Pliocene-Quaternary melting of crustal rocks occurred at depths of 15–50 km in areas where the LV–HCZs have been recognized. This provides new petrological evidence that the LV–HCZs are sources of partial melt. It is inferred that crustal melting played a key role in triggering crustal weakening and outward crustal flow in the expansion of the Tibetan Plateau.

The role of the low velocity-high conductivity zones (LV–HCZs) in developing the Tibetan Plateau has remained controversial. Here, Wang et al. present new geochemical and petrological data that show the LV–HCZs are sources of partial melt thus giving insight into the development of the Tibetan Plateau.

The boundary between the Indian and Asian tectonic plates below Tibet

The fate of the colliding Indian and Asian tectonic plates below the Tibetan high plateau may be visualized by, in addition to seismic tomography, mapping the deep seismic discontinuities, like the crust-mantle boundary (Moho), the lithosphere-asthenosphere boundary (LAB), or the discontinuities at 410 and 660 km depth. We herein present observations of seismic discontinuities with the P and S receiver function techniques beneath central and western Tibet along two new profiles and discuss the results in connection with results from earlier profiles, which did observe the LAB. The LAB of the Indian and Asian plates is well-imaged by several profiles and suggests a changing mode of India-Asia collision in the east-west direction. From eastern Himalayan syntaxis to the western edge of the Tarim Basin, the Indian lithosphere is underthrusting Tibet at an increasingly shallower angle and reaching progressively further to the north. A particular lithospheric region was formed in northern and eastern Tibet as a crush zone between the two colliding plates, the existence of which is marked by high temperature, low mantle seismic wavespeed (correlating with late arriving signals from the 410 discontinuity), poor Sn propagation, east and southeast oriented global positioning system displacements, and strikingly larger seismic (SKS) anisotropy.

Underplating in the Himalaya-Tibet collision zone revealed by the Hi-CLIMB experiment

We studied the formation of the Himalayan mountain range and the Tibetan Plateau by investigating their lithospheric structure. Using an 800-kilometer-long, densely spaced seismic array, we have constructed an image of the crust and upper mantle beneath the Himalayas and the southern Tibetan Plateau. The image reveals in a continuous fashion the Main Himalayan thrust fault as it extends from a shallow depth under Nepal to the mid-crust under southern Tibet. Indian crust can be traced to 31 degrees N. The crust/mantle interface beneath Tibet is anisotropic, indicating shearing during its formation. The dipping mantle fabric suggests that the Indian mantle is subducting in a diffuse fashion along several evolving subparallel structures.