ORIGINAL RESEARCH article

Front. Earth Sci., 30 August 2023

Sec. Structural Geology and Tectonics

Volume 11 - 2023 | https://doi.org/10.3389/feart.2023.1252670

Lithospheric structures of the eastern Himalayas as revealed by receiver function analysis

  • 1. School of Earth Sciences and Engineering, Sun Yat-sen University, Zhuhai, China

  • 2. Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai, China

Abstract

Although the Himalayan-Tibetan orogen is a result of compressional tectonics, the orogen also hosts active rifts accommodating east-west extension orthogonal to the north-south India-Asia convergence. In this study we address the question of how the north-trending rifts were formed by conducting high-resolution seismic imaging survey across southeastern Tibet where the Cona rift is exposed. Our work shows that the crustal structures of this youngest rift in southern Tibet was constructed by multiple-scale structures that are decoupled with depth and long rift trend. We suggest this deformation style to have resulted from eastward extrusion of the middle and upper crust with increasing speeds to the north towards the Yarlunbg-Zangpo suture. The differential eastward extrusion in turn may have contributed to the formation and evolution of the eastern Himalayan syntaxis.

1 Introduction

Active north-trending rifts in southern Tibet and the northern Himalaya across the Yarlung-Zangpo suture (YZS) zone between India and Asia have been known since the pioneer work of (; Yin and Harrison, 2000) (Figure 1). However, whether the initiation of the rifts was a result of 1) gravitational collapse () when the Tibetan plateau reached its maximum gravity-supported elevation (), 2) large-scale strike-slip faulting due to lateral extrusion of northern Tibet (), 3) convective removal of the Tibetan mantle lithosphere (), 4) plate-boundary-enforced continental-scale east-west extension (Yin, 2000), or 5) gravitational spreading of the Tibetan lithosphere assisted by lateral asthenospheric flow (Yin and Taylor, 2011) has been debated. Studies aiming at differentiating the above competing hypotheses helped clarified the possible dynamic causes, timing, and potential processes during the rift formation and evolution (; ; ). Despite these efforts, little progress has been made for diagnostic tests of the above models due to the lack of high-resolution imaging of the lithospheric structures across the Tibetan rifts. In this study, we address this issue by presenting newly acquired short-period seismological data to better document the rift structures relevant to differentiating the competing models.

FIGURE 1

). Dark green circles mark the earthquake hypocenters determined by the Global Centroid-Moment-Tensor (CMT) project (). (C) The receiver function image of that shows the main structures within the crust along latitude 29.2 °N. See (B) for the profile location. Red and gray colors represent positive and negative increases in seismic velocities with depth, respectively. The locations of focal mechanisms are shown in (B) (). Yellow rectangle is the distribution of ultrapotassic rocks.

2 Geological and geophysical background

East-west extension in southern Tibet may have been initiated between the middle and late Miocene (; ; ). According to the distribution and focal mechanisms of local earthquakes (Figures 1B, C), the middle crust was an aseismic zone and the normal faults were generally associated with north-trending rifts (). The earthquakes in lower crust were triggered by the strong interactions between the crust and mantle (; ; ).

Cenozoic ultrapotassic volcanic rocks and leucogranites are widely developed in the Lhasa Terrane and Tethyan Himalayan Sequence (THS) (; ); their ages range from the late Oligocene to the Miocene (; ; ). The origin of the ultrapotassic volcanic rocks may constrain the evolution of continent-continent collision zones (; ). Geochemical analyses suggest that the ultrapotassic melt was derived from the interaction between the Indian lithospheric mantle and upwelling asthenosphere (Xu et al., 2017; ). The asthenospheric upwelling was generally observed as a factor strongly related to mantle-lithosphere delamination (Zhao et al., 2009) and tearing ().

The Himalayan leucogranites are mainly exposed within the Tethyan Himalayan sequence in north and the Greater Himalayan Crystalline Complex (GHC) in south (; ; ; ). The emplacement age of leucogranites was from 44 to 7 Ma, and most leucogranites were emplaced during the Neo-Himalayan period (25-14 Ma) in both belts (). The Eo-Himalayan stage (46–25 Ma) had a few leucogranite emplacements in the northern belt (; ). proposed a Post-Himalayan stage (<14 Ma), leucogranites of this stage were sparse and were mainly discovered along N–S-trending rifts (). The in situ partial melting of the subducted GHC at depth or the fractional crystallization that occurred with the exhumation of this system are used to explain the series of characteristics of the leucogranites (Wu et al., 2015; Zhang et al., 2017). Overall, these studies showed that the evolution of the Himalayas was significantly associated with the exhumation of GHC (; ).

The South Tibetan Detachment System (STDS) is a low-angle normal fault system parallel to the entire Himalayan orogenic belt (). The STDS separated the Tethyan Himalayan Sequence from the underlying the GHC in the north-central Himalaya (), but merges to the south with the Main Central Thrust (MCT) in the western Himalaya and the southernmost Nepal Himalaya (; ; ). The initiation age of the STDS as a normal fault is estimated between 35 Ma and 23 Ma (; ). The cessation time of the STDS is mainly determined by the north-trending rifts (; Xu et al., 2013). Motion along the STDS has been related to the exhumation of the GHC (Yin, 2006).

Evidence from several geophysical studies suggests that the subducted Indian lithospheric slab was torn into pieces with different angles (; ; ). The lateral variation of the subducted slab may cause the asthenospheric upwelling, it appears to be positively related to the partial melting of the heated Indian lithospheric mantle and lower crust (). What we know about the partial eclogitization of Indian lower crust was largely based upon gravity anomalies and seismic data (; ; ; ; ; ). The anatexis of granulite and eclogite took place at ∼50 km depth (800°C–850°C) (; ) during which magmatic melts were transported upwards through the crust. The low velocity and resistivity layer, interpreted as a partially molten zone in the middle and lower crust of the Himalaya, was identified by tomographic and magnetotelluric images (; ; ; ; Xu et al., 2017; ; ). Numerical simulations support the suggestion that the underlying low-viscosity layer in the Himalaya may have assisted GHC exhumation via erosion-induced lower-crustal extrusion (). The receiver function (RF) illustrated the anisotropic shear zone at shallow depths by exhumation of GHC (). found that the Indian lower crust was gradually thinning towards the east (Figure 1C), and the rifts and Moho offset occur in pairs.

The Cona rift is located in southeast Tibet that cuts across the Tethyan Himalayan Sequence (Figure 1B). The formation mechanism of the Cona rift is controversial (; ). This study area (Figure 2A) is located in the southeast of the THS. Two gneiss domes are exposed on the east side of the Cona rift, Yardoi and Cuonadong domes. The study area is dominated by north-dipping thrust faults, from north to south, the Lhunze thrust fault, Rongbu-Gudui Fault, Lhozhag thrust and Main central fault (MCT) (). The deformation age of these tectonic units is mainly in the Oligocene and Miocene (; ).

FIGURE 2

). The stations are divided into two segments (Profiles 1 and 2). The receiver function (RF) stations are located northwest of point “P” along profile 1 and east of point “P” along profile 2. Distribution of piercing points at a depth of 70 km is shown: the purple piercing points are related to the RFs used in profile 1 and the orange dots are in profile 2. Profile 3 shows the RFs of . The labeled boxes show the range of piercing points. (B) Distribution of teleseismic events used in this study. (C) The 1-D S-wave velocity model used for CCP stacking.

To date, geophysical and petrological studies have accumulated a wealth of information about the Cenozoic evolution of the Himalaya. However, the current resolution of the available seismic data is insufficient to conform or reject the competing models for the Tibetan rift formation. In this study, we provide the highest-resolution images of the crustal structures across the Cona rift.

3 Data and methods

A receiver function (RF) is a time series containing the P-to-S converted waves and multiples formed by velocity discontinuities (; Zhu and Kanamori, 2000b). It can reveal the fine structures of the crust (). Since the common conversion point (CCP) stacking (Zhu, 2000a) was proposed, as one of the receiver function migration imaging methods, it has been widely used in the study of velocity discontinuities in the crust and upper mantle (Zhang et al., 2014).

The short-period, dense-nodal array has a higher lateral resolution of the RF image at the crustal scale compared to the broadband-station array (). Due to the high integration of node instruments, high-density deployment becomes possible (). These new instruments can be used to study the geological bodies at different scales.

The teleseismic waveforms used in this study were recorded by 148 5-Hz nodal geophones (Fairfield Zland 3C) with a spacing of 1 km and an observation period of ∼40 days (end of the year 2019). The survey line is roughly in the northwest-southeast direction (Figure 2A), crossing the main part of the Cona rift. A total of 24 teleseismic events with magnitudes ≥5.4 were observed for subsequent P-wave RF analysis. The distribution range of epicenter distance is 30°–90°, and most of the earthquakes occurred in the southeast of our study area (Figure 2B). Finally, the 1,604 high-quality RFs were calculated for imaging the crustal structures (Please see the Supplementary Material for the details of the data processing). In general, care was taken to avoid higher Gaussian filter values, because they make it more difficult to separate noises from signals (). Dense nodal array is closely spaced and can track continuous arrival times through multiple geophones (). To enhance the imaged crustal structures, we stacked the RFs by using spatial () and CCP stacking (Xu et al., 2018) methods. The CCP stacking process involves two main steps. Initially, every amplitude on the receiver function subsequent to the direct P-wave is presumed to originate from a solitary P-to-S conversion occurring along the ray path. This amplitude is then projected to the conversion point utilizing a reference velocity model. Subsequently, the volume of the crust is partitioned into bins of specified dimensions, and all amplitudes within a given bin are averaged to generate a structural image. The 1-D S-wave velocity model used for CCP stacking from Xin et al. (2019) (Figure 2C). Our CCP image of crustal structure used bins 3-km long and 1-km thick.

4 Results

Our survey line intersects obliquely the regional tectonic trends (Figure 2A). In order to separate east-trending compressional structures from north-trending extensional structures along the rifting zone, our imaging results were projected onto two separate profiles (Figure 2A), respectively. This helped reveal the relationship between structures responsible for GHC exhumation and structures associated with north-trending rifting.

The continuity of the crustal structures (grey thin dash lines in Figure 3) increases with depth in profile 1. The deformation of the middle and upper crust was strong, and the structure of the lower crust is more continuous (Figures 3A, C). The Moho discontinuity is traceable north of the Cona rift, the projected Indian crust dips to the north, and the crustal thickness increases from ∼70 km to ∼80 km (Figures 3A, C). The lower crust is overlain by a set of wide and gentle fold-thrust structures, and the upper crust displays a series of north-dipping thrusts associated closely with fault-bend folds (Figure 3C). There are two scales of folds above and below a 20-km deep interface (Figure 3C), which can be traced to the STDS (thick red dash line in Figures 4A, B). The undeformed lower crust was decoupled from the overlying fold-thrust structures (dark blue zone in Figures 4A, C), there was a decollement (thin pink dash line in Figures 4A, C) between them. Based on previous studies (; ), this interface is interpreted as the Main Himalayan Thrust (MHT). Extensively developed thrusts indicate that the seismically imaged area has been strongly affected by north-south shortening. A possible normal fault was interpreted in the Cona rift zone, but few obvious normal faults can be established in this study. The short-wavelength folds beneath the STDS show the characteristics of low viscosity. The low-velocity zones were highlighted in Figures 3A, C, and the leucogranites were widely exposed in this area. Hence, the low-velocity zones were possible channels for exhumation of high-grade metamorphic rocks. A continuous interface appears in the northern part of Profile 1 at a depth of 70 km (Figure 3C), which might be the partial eclogitization of Indian lower crust by comparing previous study () in the same location.

FIGURE 3

FIGURE 4

) at the approximate locations are projected on (A, B), (D, E). The isotherms in (C) refer to .

The Cona rift is north-trending normal faults, Profile 2 had been imaged to get a better view of its deep structures (Figures 3B, D). The result of low Gaussian coefficient (2.5) is relatively simple, and the structures located in the middle crust are overall east-dipping (Figure 3B). The structures shown in high Gaussian coefficient (5.5) were a series of short flat layers (Figure 3D). Almost all strata remain horizontal, without any dip to the east or west (Figure 3D). However, there were a large number of offsets in these strata. The strata show a gradual stepping towards the east. This is also consistent with the overall eastward dipping structure shown by the low Gaussian coefficient result.

Since the Moho was rather heterogeneous in the southern section of the survey line, we conducted partition analysis for different locations of the survey line (the boxes marked with numbers, and please see Figure 2 for locations). There was a data intersection between the south section of Profile 1 and the west side of Profile 2 (boxes marked 4, 5 and 6 in Figure 2A). There is a very strong positive amplitude at a depth of 60–80 km at the southern end of Profile 1 (Figure 3C), the west side of Profile 2 did not have the similar amplitude. There were also several strong positive amplitudes on the west side of Profile 2 (dash lines in boxes marked 4 and 5 in Figure 3D). By comparing the interface thickness (green and yellow bars in Figures 3C, D) and the continuity of Moho in profile 1, we can confirm that the western side of profile 2 has the Moho feature at about 70 km depth (Figure 3D), it helped us identify the Moho at the southern end of Profile 1. There is a robust Moho feature on the east side of Profile 2 (boxes marked 1, 2 and 3 in Figure 3D).

The comparison between Profile 2 (Figure 3D) and Profile 3 (Figure 3E) shows that the crust has the following characteristics in the east-west direction. The depth of Moho is deeper in the west than in the east, and there are offsets in the middle. Taking the offset as the reference point, the surface position of the Cona Rift in the northern section is migrated eastward. The lower crust is unevenly thick on both sides of the offset, and thicker on the western side. The images of the low Gaussian coefficient (Figures 3B, E) indicate that the upper and middle crust is tilted to the east.

5 Discussion

By comparing our results with the finding of , the crustal thickness on the west side in our result is thicker than that on the east side (Figure 4E). A large number of regularly arranged folds located above the MHT (Figures 4A, C) are interpreted as shear zones (; ). The fold trend shows the exhumation and shortening of the upper and middle crust in the north-south direction.

There is a distinct wedge in the upper and middle crust to the west (Figure 4B), the wedge is actually due to an increase in the number of strata to the west (light blue area and light pink area Supplementary Figure S10B). The horizontal stratum had also been imaged in profile 1 (light pink area in Supplementary Figure S10A), the ray coverage area (blue box marked 6 in Figure 2A) was similar to profile 2 (black boxes marked 4 and 5 in Figure 2A). The extra stratum on the west side was likely to be the exhumation of the GHC. The uplift of upper crustal (above the STDS) folds was produced by crustal shortening, the normal faults indicate the relative movements of extrusion (Figure 4D). But the upper crustal folds are gentler in the east-west direction, and the middle (GHC) and lower curst is flat with offsets in the east direction (Figures 4B, D). Comparing the depth and thickness of the middle crustal structures, the fold-thrust structures thicken the middle crust and the offsets of the Moho are in response to the lateral inhomogeneity of the middle crustal thickening.

The GHC generally has experienced multiple phases of partial melting expressed by the emplacement of leucogranites (; ). The partial melting weakens the mid-to-lower crust (), favoring flow and exhumation under the combined action of plateau-induced pressure gradients and the density difference between molten and un-molten rocks (). Our results (Figures 4C, D) shows that the subducted Indian crust can reach the depth over 60km, the P-T condition at this depth made the eclogitization and partial melting of lower and middle crust happen (; ; ; ; ). The partial melting of granulite- and eclogite-facies metamorphism initiated under the P-T conditions of 1.2 GPa and 800°C–850°C, respectively (Figure 4C) (Zhang et al., 2017; ). The partial melting of subducted GHC due to heating from below triggered the exhumation of high-grade metamorphosed rocks (). The MHT was the interface between the exhumation and the underthrusting Indian lower crust.

The exhumation process may have been partitioned (Figure 4C): the lower GHC section was by tectonic wedging () while the upper GHC section was by channel-flow tunnel (; ). The leucogranites were the products of partial melting of subducted the GHC and subsequent highly fractionated magma (; ; Wu et al., 2015). Although we are unable to identify rock lithology by our seismic imaging techniques, the inferred partial melting in the models of could be represented by the presence of low velocity zones (; Yang et al., 2012; ). The likely low velocity zone may be located by the polarity converter (negative-over-positive) in the upper crust and the upper part of the middle crust (Figures 3A–D). The migration of the melting magmas went through a long distance (Figures 4A, C), it helps researchers to explain the different stages and locations of exposed leucogranites, and also the formation of the leucogranites through a high degree of fractionation along the STDS ().

A significant finding of this study is that the deformation of the Indian crust is lateral variation (Figures 4B, D). The crustal thickening was caused by the crustal shortening and the exhumation of the GHC, but the time of exhumation can be different. The time was controlled by the dipping angles of the Indian subducting continental lithosphere; higher angle made the GHC reach the high P-T condition earlier. The Indian continental lithosphere began to subduct under the pulling of the subducting Neo-Tethys Oceanic lithosphere (Zhang et al., 2012). The age of the subducted Neo-Tethys oceanic slab breakoff is relatively consistent, which is 45–30 Ma (). The subducting Indian continental lithosphere without the slab pull continues to subduct beneath the Eurasian plate due to inertia (). Evidence from the Southern Tibet Rift System distribution, magmatic rock composition, or geophysical observations all show that the subducting Indian continental lithospheric slab is torn in a north-south direction (; Yin and Taylor, 2011; ). The change in dip of the subducting slab may have controlled the stress state in the overlying plate (). Lateral variation in the angle of the subducted slab may have caused slab tearing (Figures 5B, C). Different from the previous slab tearing model (Yin, 2000), the slab window (Figure 5C) formed by the slab tearing allows the hot asthenosphere material to erode and destroy the upper lithosphere, which promotes the development of the north-trending rifts. Slab tear started at 25 Ma (), and it was torn into at least four parts (). The inhomogeneity of subducting slab and mantle, lateral slab breakoff, and subduction duration are the reasons why the slab was torn (; ). Furthermore, the subduction angle of the subducting Indian continental lithosphere increases from west to east (). The regional seismic tomography shows the geometry of the subducting Indian continental lithosphere (). Based on the distribution and age of the surface outcrops of ultramafic magmatism () and the locations of the slab tears, we infer that the slab tears occur sequentially from west to east. Variation in the regional vertical stress was caused by changes in the angle of subduction of the Indian lithospheric mantle (Yin, 2006). The contact between the overlying Indian crust and the upwelling high temperature mantle promoted the partial melting of the GHC().

FIGURE 5

). (B) The crustal shortening and exhumation of GHC led to the lateral variation of crustal thickening. The middle and upper crust extruded eastward, and the Cona rift migrated eastward with it. (C) The slab tearing allowed the hot asthenosphere material to upwell and melt. The eastward escape of the exhumation of the GHC reshaped the crustal structure in the rift region.

Helium isotope data from hot springs in southern Tibet have been used to infer the location of the mantle suture between India and Asia () (Figure 5A). Surface-collected helium isotopes sourced from the mantle must have been transported through the crust via complex conduit networks that may have been strongly controlled by the north-trending rifts due to the spatial correlation of the springs and rifts. In the southern part of the Cona rift, the crustal structure of the rift is relatively simple, with the structures traceable continuously from the Moho to the surface. It shows an east-dipping structure that may provide a conduit for helium isotope to rise to the surface (Figure 5B). In the northern segment, the structure of the rift at mid-crustal depth is replaced by a series of east-dipping structures, and the location of the rift at the surface is farther east than the Moho offset (Figure 5C). This relationship may imply that the rift structures did not form a conduit network for helium isotopes to rise from the mantle to the surface. From south to north, the Cona rift changes its trend from north to northeast direction (Figure 5A). This may indicate that the upper and middle crust in the eastern Tethyan Himalaya may have been extruded eastward with an increasing eastward velocity towards the Yarlung-Zangpo suture zone.

6 Conclusion

A high-resolution receiver-function study in southeastern Tibet reveals geological structures at different scales across a north-trending rift zone (i.e., the Cona rift). The intensity of the crustal deformation decreases with increasing depth, which may have been resulted from the formation of upper and lower decollements as presented by the STDS and MHT, respectively. The unusually thick crustal thickness was contributed to crustal shortening and exhumation of the GHC. The GHC exhumation is the tectonic wedging accompanied by channel tunneling model. The distribution of exposed leucogranites were traceable to the migration of the partial melts. Stress changed due to slab tearing; the exhumation was not synchronized from west to east. At the present, it is possible that extensional structures along the Cona rift are decoupled with depth and rift strike. The eastward extrusion of the middle and upper crust may have different velocities, and this difference may contribute to the evolution of eastern Himalayan syntaxis.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

BX, XX, XT, and XG planned the field observation scheme. BX, XX, XT, JW, CL, XL, and JY carried out the field work and collected the waveforms. BX, XX, XT, JW and XG analyzed the data and wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

Lithospheric Structures of the Eastern Himalayas as Revealed by Receiver Function Analysis has received funding from the National Natural Science Foundation of China (Grant No 41874102, 42325402, 42274120, 41974097).

Acknowledgments

The authors appreciate the editors and reviewers for their valuable comments and suggestions which greatly helped us to improve the present manuscript. We would like to thank professor An Yin for his help with the English expression.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2023.1252670/full#supplementary-material

References

  • 1

    ArmijoR.TapponnierP.HanT. (1989). Late Cenozoic right-lateral strike-slip faulting in southern Tibet. J. Geophys. Res. Solid Earth94, 27872838. 10.1029/JB094iB03p02787

  • 2

    BeaumontC.JamiesonR. A.NguyenM. H.LeeB. (2001). Himalayan tectonics explained by extrusion of a low-viscosity crustal channel coupled to focused surface denudation. Nature414, 738742. 10.1038/414738a

  • 3

    BeaumontC.JamiesonR. A.NguyenM. H.MedvedevS. (2004). Crustal channel flows: 1. Numerical models with applications to the tectonics of the himalayan-Tibetan orogen. J. Geophys. Res. Solid Earth109. 10.1029/2003jb002809

  • 4

    BianS.GongJ.ZuzaA. V.YangR.TianY.JiJ.et al (2020). Late Pliocene onset of the Cona rift, eastern Himalaya, confirms eastward propagation of extension in Himalayan-Tibetan orogen. Earth Planet. Sci. Lett.544, 116383. 10.1016/j.epsl.2020.116383

  • 5

    BurchfielB. C.RoydenL. H. (1985). North-south extension within the convergent Himalayan region. Geology13, 679682. 10.1130/0091-7613(1985)13<679:NEWTCH>2.0.CO;2

  • 6

    CarosiR.LombardoB.MolliG.MusumeciG.PertusatiP. (1998). The South Tibetan detachment system in the Rongbuk valley, Everest region. Deformation features and geological implications. J. Asian Earth Sci.16, 299311. 10.1016/S0743-9547(98)00014-2

  • 7

    CarosiR.MontomoliC.IaccarinoS.MassonneH. J.RubattoD.LangoneA.et al (2016). Middle to late eocene exhumation of the greater himalayan sequence in the central Himalayas: progressive accretion from the Indian plate. Geol. Soc. Am. Bull.128, 15711592. 10.1130/B31471.1

  • 8

    ChenJ.GaillardF.VillarosA.YangX.LaumonierM.JolivetL.et al (2018). Melting conditions in the modern Tibetan crust since the Miocene. Nat. Commun.9, 3515. 10.1038/s41467-018-05934-7

  • 9

    ChenW.-P.KaoH. J. W.GeologyR. (1996). Seismotectonics of Asia. Some recent Prog.1, 3762.

  • 10

    ChenY.LiW.YuanX.BadalJ. H.TengJ. W. (2015). Tearing of the Indian lithospheric slab beneath southern Tibet revealed by SKS-wave splitting measurements. Earth Planet. Sci. Lett.413, 1324. 10.1016/j.epsl.2014.12.041

  • 11

    ChungS.-L.ChuM.-F.ZhangY.XieY.LoC.-H.LeeT.-Y.et al (2005). Tibetan tectonic evolution inferred from spatial and temporal variations in post-collisional magmatism. Earth-Science Rev.68, 173196. 10.1016/j.earscirev.2004.05.001

  • 12

    ChungS.-L.LiuD.JiJ.ChuM.-F.LeeH.-Y.WenD.-J.et al (2003). Adakites from continental collision zones: melting of thickened lower crust beneath southern Tibet. Geology31, 1021. 10.1130/g19796.1

  • 13

    DoglioniC.CarminatiE.CuffaroM.ScroccaD. (2007). Subduction kinematics and dynamic constraints. Earth-Science Rev.83, 125175. 10.1016/j.earscirev.2007.04.001

  • 14

    DongX.LiW.LuZ.HuangX.GaoR. (2020). Seismic reflection imaging of crustal deformation within the eastern Yarlung-Zangbo suture zone. Tectonophysics780, 228395. 10.1016/j.tecto.2020.228395

  • 15

    EdwardsM.HarrisonT. (1997). When did the roof collapse? Late Miocene north-south extension in the high Himalaya revealed by Th-Pb monazite dating of the khula kangri granite. Geology25, 543546. 10.1130/0091-7613(1997)025<0543:WDTRCL>2.3.CO;2

  • 16

    EnglandP.HousemanG. (1989). Extension during continental convergence, with application to the Tibetan Plateau. J. Geophys. Res.94, 17561. 10.1029/JB094iB12p17561

  • 17

    FuJ.LiG.WangG.ZhangL.LiangW.ZhangX.et al (2020). Structural analysis of sheath folds and geochronology in the Cuonadong Dome, southern Tibet, China: new constraints on the timing of the South Tibetan detachment system and its relationship to North Himalayan Gneiss Domes. Terra nova.32, 300323. 10.1111/ter.12462

  • 18

    GrujicD.CaseyM.DavidsonC.HollisterL. S.KündigR.PavlisT.et al (1996). Ductile extrusion of the higher himalayan crystalline in Bhutan: evidence from quartz microfabrics. Tectonophysics260, 2143. 10.1016/0040-1951(96)00074-1

  • 19

    GuoZ.GaoX.YaoH.LiJ.WangW. (2009). Midcrustal low-velocity layer beneath the central Himalaya and southern Tibet revealed by ambient noise array tomography. Geochem. Geophys. Geosystems10. 10.1029/2009gc002458

  • 20

    GuoZ.WilsonM. (2012). The Himalayan leucogranites: constraints on the nature of their crustal source region and geodynamic setting. Gondwana Res.22, 360376. 10.1016/j.gr.2011.07.027

  • 21

    HackerB. R.KelemenP. B.BehnM. D. (2015). Continental lower crust. Annu. Rev. Earth Planet. Sci.43, 167205. 10.1146/annurev-earth-050212-124117

  • 22

    HarrisonT. M.CopelandP.KiddW.LoveraO. M. (1995). Activation of the Nyainqentanghla shear zone: implications for uplift of the southern Tibetan Plateau. Tectonics14, 658676. 10.1029/95TC00608

  • 23

    HazarikaD.SenK.KumarN. (2014). Characterizing the intracrustal low velocity zone beneath northwest India–Asia collision zone. Geophys. J. Int.199, 13381353. 10.1093/gji/ggu328

  • 24

    HetényiG.VergneJ.BollingerL.CattinR. (2011). Discontinuous low-velocity zones in southern Tibet question the viability of the channel flow model. Geol. Soc. Lond. Spec. Publ.353, 99108. 10.1144/SP353.6

  • 25

    HouZ. Q.GaoY. F.QuX. M.RuiZ. Y.MoX. X. (2004). Origin of adakitic intrusives generated during mid-Miocene east–west extension in southern Tibet. Earth Planet. Sci. Lett.220, 139155. 10.1016/s0012-821x(04)00007-x

  • 26

    HouZ. Q.ZhaoZ. D. D.GaoY. F.YangZ. M.JiangW. (2006). Tearing and dischronal subduction of the Indian continental slab: evidence from cenozoic Gangdese volcano-magmatic rocks in south Tibet. Acta Petrol Sin.22, 761774. 10.1016/j.sedgeo.2005.11.021

  • 27

    JamiesonR. A.BeaumontC.NguyenM.GrujicD. (2006). Provenance of the greater himalayan sequence and associated rocks: predictions of channel flow models. Geol. Soc. Lond. Spec. Publ.268, 165182. 10.1144/GSL.SP.2006.268.01.07

  • 28

    KindR.YuanX.SaulJ.NelsonD.SobolevS. V.MechieJ.et al (2002). Seismic images of crust and upper mantle beneath Tibet: evidence for Eurasian plate subduction. Science298, 12191221. 10.1126/science.1078115

  • 29

    KlempererS. L.ZhaoP.WhyteC. J.DarrahT. H.CrosseyL. J.KarlstromK. E.et al (2022). Limited underthrusting of India below Tibet: (3)He/(4)He analysis of thermal springs locates the mantle suture in continental collision. Proc. Natl. Acad. Sci. U. S. A.119, e2113877119. 10.1073/pnas.2113877119

  • 30

    LangstonC. A. (1979). Structure under Mount Rainier, Washington, inferred from teleseismic body waves. J. Geophys. Res.84, 4749. 10.1029/JB084iB09p04749

  • 31

    LeahyG. M.SaltzerR. L.SchmedesJ. (2012). Imaging the shallow crust with teleseismic receiver functions. Geophys. J. Int.191, 627636. 10.1111/j.1365-246X.2012.05615.x

  • 32

    LeeJ.WhitehouseM. J. (2007). Onset of mid-crustal extensional flow in southern Tibet: evidence from U/Pb zircon ages. Geology35, 4548. 10.1130/G22842A.1

  • 33

    LiC.van der HilstR. D.MeltzerA. S.EngdahlE. R. (2008). Subduction of the Indian lithosphere beneath the Tibetan plateau and Burma. Earth Planet. Sci. Lett.274, 157168. 10.1016/j.epsl.2008.07.016

  • 34

    LiJ.SongX. (2018). Tearing of Indian mantle lithosphere from high-resolution seismic images and its implications for lithosphere coupling in southern Tibet. Proc. Natl. Acad. Sci. U. S. A.115, 82968300. 10.1073/pnas.1717258115

  • 35

    LinF.-C.LiD.ClaytonR. W.HollisD. (2013). High-resolution 3D shallow crustal structure in Long Beach, California: application of ambient noise tomography on a dense seismic array. Geophysics78, Q45Q56. 10.1190/geo2012-0453.1

  • 36

    LiuZ.-C.WuF.-Y.DingL.LiuX.-C.WangJ.-G.JiW.-Q. (2016). Highly fractionated late eocene (∼35 Ma) leucogranite in the xiaru dome, tethyan Himalaya, south Tibet. Lithos240-243, 337354. 10.1016/j.lithos.2015.11.026

  • 37

    LiuZ.-C.WuF.-Y.JiW.-Q.WangJ.-G.LiuC.-Z. (2014). Petrogenesis of the Ramba leucogranite in the Tethyan Himalaya and constraints on the channel flow model. Lithos208-209, 118136. 10.1016/j.lithos.2014.08.022

  • 38

    LiuZ.-C.WuF.-Y.QiuZ.-L.WangJ.-G.LiuX.-C.JiW.-Q.et al (2017a). Leucogranite geochronological constraints on the termination of the South Tibetan Detachment in eastern Himalaya. Tectonophysics721, 106122. 10.1016/j.tecto.2017.08.019

  • 39

    LiuZ.TianX.GaoR.WangG.WuZ.ZhouB.et al (2017b). New images of the crustal structure beneath eastern Tibet from a high-density seismic array. Earth Planet. Sci. Lett.480, 3341. 10.1016/j.epsl.2017.09.048

  • 40

    LiuZ.TianX.YuanX.LiangX.ChenY.ZhuG.et al (2020a). Complex structure of upper mantle beneath the Yadong-Gulu rift in Tibet revealed by S-to-P converted waves. Earth Planet. Sci. Lett.531, 115954. 10.1016/j.epsl.2019.115954

  • 41

    LiuZ.WangJ. G.LiuX. C.LiuY.LaiQ. Z. (2020b). Middle Miocene ultrapotassic magmatism in the Himalaya: A response to mantle unrooting process beneath the orogen. Terra nova.33, 240251. 10.1111/ter.12507

  • 42

    MahéoG.GuillotS.Blichert-ToftJ.RollandY.PêcherA. (2002). A slab breakoff model for the Neogene thermal evolution of South Karakorum and South Tibet. Earth Planet. Sci. Lett.195, 4558. 10.1016/S0012-821X(01)00578-7

  • 43

    Mark HarrisonT.LoveraO. M.GroveM. (1997). New insights into the origin of two contrasting Himalayan granite belts. Geology25, 899902. 10.1130/0091-7613(1997)025<0899:NIITOO>2.3.CO;2

  • 44

    MitsuishiM.WallisS. R.AoyaM.LeeJ.WangY. (2012). E–W extension at 19 Ma in the Kung Co area, S. Tibet: evidence for contemporaneous E–W and N–S extension in the Himalayan orogen. Earth Planet. Sci. Lett.325-326, 1020. 10.1016/j.epsl.2011.11.013

  • 45

    MolnarP.ChenW.-P. (1983). Focal depths and fault plane solutions of earthquakes under the Tibetan Plateau. J. Geophys. Res.88, 1180. 10.1029/JB088iB02p01180

  • 46

    MolnarP.TapponnierP. (1978). Active tectonics of Tibet. J. Geophys. Res.83, 5361. 10.1029/JB083iB11p05361

  • 47

    MurphyM.YinA.HarrisonT.DurrS.RyersonF.KiddW.et al (1997). Did the Indo-Asian collision alone create the Tibetan plateau?Geology25, 719722. 10.1130/0091-7613(1997)025<0719:DTIACA>2.3.CO;2

  • 48

    NabelekJ.HetenyiG.VergneJ.SapkotaS.KafleB.JiangM.et al (2009). Underplating in the Himalaya-Tibet collision zone revealed by the Hi-CLIMB experiment. Science325, 13711374. 10.1126/science.1167719

  • 49

    NelsonK. D.ZhaoW.BrownL.KuoJ.CheJ.LiuX.et al (1996). Partially molten middle crust beneath southern Tibet: synthesis of project INDEPTH results. science274, 16841688. 10.1126/science.274.5293.1684

  • 50

    O'BrienP. J. (2018). Eclogites and other high-pressure rocks in the Himalaya: A review. Geol. Soc. Lond. Spec. Publ.483, 183213. 10.1144/sp483.13

  • 51

    PangY. J.ZhangH.GeryaT. V.LiaoJ.ChengH. H.ShiY. L. (2018). The mechanism and dynamics of N-S rifting in southern Tibet: insight from 3-D thermomechanical modeling. J. Geophys Res-Sol Ea.123, 859877. 10.1002/2017jb014011

  • 52

    RubattoD.ChakrabortyS.DasguptaS. (2012). Timescales of crustal melting in the Higher Himalayan Crystallines (Sikkim, Eastern Himalaya) inferred from trace element-constrained monazite and zircon chronology. Contributions Mineralogy Petrology165, 349372. 10.1007/s00410-012-0812-y

  • 53

    Schulte-PelkumV.MonsalveG.SheehanA.PandeyM. R.SapkotaS.BilhamR.et al (2005). Imaging the Indian subcontinent beneath the Himalaya. Nature435, 12221225. 10.1038/nature03678

  • 54

    ShiD.KlempererS. L.ShiJ.WuZ.ZhaoW. (2020). Localized foundering of Indian lower crust in the India-Tibet collision zone. Proc. Natl. Acad. Sci. U. S. A.117, 2474224747. 10.1073/pnas.2000015117

  • 55

    ShiD.WuZ.KlempererS. L.ZhaoW.XueG.SuH. (2015). Receiver function imaging of crustal suture, steep subduction, and mantle wedge in the eastern India–Tibet continental collision zone. Earth Planet. Sci. Lett.414, 615. 10.1016/j.epsl.2014.12.055

  • 56

    SorcarN.HoppeU.DasguptaS.ChakrabortyS. (2014). High-temperature cooling histories of migmatites from the high himalayan crystallines in Sikkim, India: rapid cooling unrelated to exhumation?Contributions Mineralogy Petrology167, 957. 10.1007/s00410-013-0957-3

  • 57

    StyronR. H.TaylorM. H.SundellK. E.StockliD. F.OalmannJ. A. G.MöllerA.et al (2013). Miocene initiation and acceleration of extension in the South Lunggar rift, western Tibet: evolution of an active detachment system from structural mapping and (U-Th)/He thermochronology. Tectonics. 10.1002/tect.20053

  • 58

    SubediS.HetényiG.VergneJ.BollingerL.Lyon‐CaenH.FarraV.et al (2018). Imaging the Moho and the main himalayan thrust in western Nepal with receiver functions. Geophys. Res. Lett.45. 10.1029/2018gl080911

  • 59

    TapponnierP.ZhiqinX.RogerF.MeyerB.ArnaudN.WittlingerG.et al (2001). Oblique stepwise rise and growth of the Tibet plateau. Science294, 16711677. 10.1126/science.105978

  • 60

    TaylorM.YinA.RyersonF. J.KappP.DingL. (2003). Conjugate strike-slip faulting along the Bangong-Nujiang suture zone accommodates coeval east-west extension and north-south shortening in the interior of the Tibetan Plateau. Tectonics22. 10.1029/2002tc001361

  • 61

    TianX. B.ChenY.TsengT. L.KlempererS. L.ThyboH.LiuZ.et al (2015). Weakly coupled lithospheric extension in southern Tibet. Earth Planet. Sci. Lett.430, 171177. 10.1016/j.epsl.2015.08.025

  • 62

    TurnerS.ArnaudN.LiuJ.RogersN.HawkesworthC.HarrisN.et al (1996). Post-collision, shoshonitic volcanism on the Tibetan plateau: implications for convective thinning of the lithosphere and the source of ocean island basalts. J. Petrology37, 4571. 10.1093/petrology/37.1.45

  • 63

    UnsworthM. J.JonesA. G.WeiW.MarquisG.GokarnS. G.SprattJ. E.et al (2005). Crustal rheology of the Himalaya and Southern Tibet inferred from magnetotelluric data. Nature438, 7881. 10.1038/nature04154

  • 64

    VanceD.HarrisN. (1999). Timing of prograde metamorphism in the zanskar Himalaya. Geology27, 395398. 10.1130/0091-7613(1999)027<0395:Topmit>2.3.Co;2

  • 65

    WangJ.-M.WuF.-Y.RubattoD.LiuK.ZhangJ.-J.LiuX.-C. (2018). Early Miocene rapid exhumation in southern Tibet: insights from P–T–t–D–magmatism path of Yardoi dome. Lithos304-307, 3856. 10.1016/j.lithos.2018.02.003

  • 66

    WardK. M.LinF.SchmandtB. (2018). High‐resolution receiver function imaging across the cascadia subduction zone using a dense nodal array. Geophys. Res. Lett.45 (12), 218212. 10.1029/2018gl079903

  • 67

    WebbA. A. (2011b). Cenozoic tectonic history of the Himachal Himalaya (northwestern India) and its constraints on the formation mechanism of the Himalayan orogen. Geosphere7, 1013. 10.1130/ges00627.1

  • 68

    WebbA. A. G.GuoH.CliftP. D.HussonL.MüllerT.CostantinoD.et al (2017). The Himalaya in 3D: slab dynamics controlled mountain building and monsoon intensification. Lithosphere, L636.1. 10.1130/l636.1

  • 69

    WebbA. A. G.SchmittA. K.HeD.WeigandE. L. (2011a). Structural and geochronological evidence for the leading edge of the Greater Himalayan Crystalline complex in the central Nepal Himalaya. Earth Planet. Sci. Lett.304, 483495. 10.1016/j.epsl.2011.02.024

  • 70

    WebbA. A. G.YinA.HarrisonT. M.CélérierJ.BurgessW. P. (2007). The leading edge of the Greater Himalayan Crystalline complex revealed in the NW Indian Himalaya: implications for the evolution of the Himalayan orogen. Geology35, 955. 10.1130/g23931a.1

  • 71

    WilliamsH. M. (2004). Nature of the source regions for post-collisional, potassic magmatism in southern and northern Tibet from geochemical variations and inverse trace element modelling. J. Petrology45, 555607. 10.1093/petrology/egg094

  • 72

    WilliamsH.TurnerS.KelleyS.HarrisN.,(2001). Age and composition of dikes in Southern Tibet: new constraints on the timing of east-west extension and its relationship to postcollisional volcanism. Geology. 29, 339342. 10.1130/0091-7613(2001)029<0339:Aacodi>2.0.Co;2

  • 73

    WittlingerG.FarraV.HetényiG.VergneJ.NábělekJ. (2009). Seismic velocities in southern Tibet lower crust: A receiver function approach for eclogite detection. Geophys. J. Int.177, 10371049. 10.1111/j.1365-246X.2008.04084.x

  • 74

    WuC.TianX.XuT.LiangX.ChenY.TaylorM.et al (2019). Deformation of crust and upper mantle in central Tibet caused by the northward subduction and slab tearing of the Indian lithosphere: new evidence based on shear wave splitting measurements. Earth Planet. Sci. Lett.514, 7583. 10.1016/j.epsl.2019.02.037

  • 75

    WuF.-Y.LiuX.-C.LiuZ.-C.WangR.-C.XieL.WangJ.-M.et al (2020). Highly fractionated Himalayan leucogranites and associated rare-metal mineralization. Lithos352-353, 105319105353. 10.1016/j.lithos.2019.105319

  • 76

    WuF.LiuZ.LiuX.JiW. (2015). Himalayan leucogranite: petrogenesis and implications to orogenesis and plateau uplift. Acta Petrol Sin.31, 136. 10.1016/j.jseaes.2014.10.014

  • 77

    XuB.GriffinW. L.XiongQ.HouZ. Q.O'ReillyS. Y.GuoZ.et al (2017a). Ultrapotassic rocks and xenoliths from South Tibet: contrasting styles of interaction between lithospheric mantle and asthenosphere during continental collision. Geology45, 5154. 10.1130/G38466.1

  • 78

    XuM.HuangH.HuangZ.WangP.WangL.XuM.et al (2018). Insight into the subducted Indian slab and origin of the Tengchong volcano in SE Tibet from receiver function analysis. Earth Planet. Sci. Lett.482, 567579. 10.1016/j.epsl.2017.11.048

  • 79

    XuQ.ZhaoJ.YuanX.LiuH.PeiS. (2017b). Detailed configuration of the underthrusting Indian lithosphere beneath western Tibet revealed by receiver function images. J. Geophys. Res. Solid Earth122, 82578269. 10.1002/2017jb014490

  • 80

    XuZ.WangQ.PêcherA.LiangF.QiX.CaiZ.et al (2013). Orogen‐parallel ductile extension and extrusion of the Greater Himalaya in the late Oligocene and Miocene. Tectonics32, 191215. 10.1002/tect.20021

  • 81

    YangY.RitzwollerM. H.ZhengY.ShenW.LevshinA. L.XieZ. (2012). A synoptic view of the distribution and connectivity of the mid-crustal low velocity zone beneath Tibet. J. Geophys. Res. Solid Earth117. 10.1029/2011jb008810

  • 82

    YinA. (2006). Cenozoic tectonic evolution of the Himalayan orogen as constrained by along-strike variation of structural geometry, exhumation history, and foreland sedimentation. Earth-Science Rev.76, 1131. 10.1016/j.earscirev.2005.05.004

  • 83

    YinA.HarrisonT. M. (2000). Geologic evolution of the Himalayan-Tibetan orogen. Annu. Rev. Earth Planet. Sci.28, 211280. 10.1146/annurev.earth.28.1.211

  • 84

    YinA. (2000). Mode of Cenozoic east-west extension in Tibet suggesting a common origin of rifts in Asia during the Indo-Asian collision. J. Geophys. Res. Solid Earth105, 2174521759. 10.1029/2000jb900168

  • 85

    YinA.TaylorM. H. (2011). Mechanics of V-shaped conjugate strike-slip faults and the corresponding continuum mode of continental deformation. Geol. Soc. Am. Bull.123, 17981821. 10.1130/B30159.1

  • 86

    ZhangK. J.ZhangY. X.TangX. C.XiaB. (2012). Late Mesozoic tectonic evolution and growth of the Tibetan plateau prior to the Indo-Asian collision. Earth-Science Rev.114, 236249. 10.1016/j.earscirev.2012.06.001

  • 87

    ZhangR.WuQ.SunL.HeJ.GaoZ. (2014). Crustal and lithospheric structure of Northeast China from S-wave receiver functions. Earth Planet. Sci. Lett.401, 196205. 10.1016/j.epsl.2014.06.017

  • 88

    ZhangZ.XiangH.DongX.LiW.DingH.GouZ.et al (2017). Oligocene HP metamorphism and anatexis of the higher himalayan crystalline sequence in yadong region, east-central Himalaya. Gondwana Res.41, 173187. 10.1016/j.gr.2015.03.002

  • 89

    ZhaoZ. D.MoX. X.DilekY.NiuY. L.DePaoloD. J.RobinsonP.et al (2009). Geochemical and Sr-Nd-Pb-O isotopic compositions of the post-collisional ultrapotassic magmatism in SW Tibet: petrogenesis and implications for India intra-continental subduction beneath southern Tibet. Lithos113, 190212. 10.1016/j.lithos.2009.02.004

  • 90

    ZhuL. (2000a). Crustal structure across the San Andreas Fault, southern California from teleseismic converted waves. Earth Planet. Sci. Lett.179, 183190. 10.1016/S0012-821x(00)00101-1

  • 91

    ZhuL.KanamoriH. (2000b). Moho depth variation in southern California from teleseismic receiver functions. J. Geophys. Res. Solid Earth105, 29692980. 10.1029/1999jb900322

Summary

Keywords

north-trending rifts, receiver function, Cona rift, dense seismic array, Moho offset

Citation

Xiang B, Xu X, Tong X, Wu J, Li C, Yu J, Luo X and Guo X (2023) Lithospheric structures of the eastern Himalayas as revealed by receiver function analysis. Front. Earth Sci. 11:1252670. doi: 10.3389/feart.2023.1252670

Received

04 July 2023

Accepted

21 August 2023

Published

30 August 2023

Volume

11 - 2023

Edited by

Tao Xu, Chinese Academy of Sciences (CAS), China

Reviewed by

Wenhui Li, Chinese Academy of Geological Sciences (CAGS), China

Zhen Guo, Southern University of Science and Technology, China

Updates

Copyright

*Correspondence: Xiao Xu,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics