Imaging the Indian subcontinent beneath the Himalaya

Seismic structure of the 2015 Mw7.8 Gorkha earthquake revealed by ambient seismic noise and teleseismic surface wave tomography

The destructive 2015 Mw7.8 Gorkha earthquake occurred in the Main Himalayan Thrust due to the collision of the Indian and Asian plates, which provides a unique opportunity to understand the deep dynamic processes and seismogenic mechanisms of strong earthquakes. We construct a regional-scale shear-wave velocity model of the crust and uppermost mantle using ambient seismic noise and teleseismic surface wave at periods of 5-100 s around the Gorkha earthquake region. The new shear-wave velocity model exhibits prominently lateral heterogeneities in the Gorkha earthquake areas. We observe a high-velocity (high-V) zone around the Gorkha main shock in the Main Himalayan Thrust, indicating the existence of a high-strength asperity that sustains the stress accumulating. The aftershocks are primarily located in the low-velocity (low-V) anomalies and enclosed by two high-V anomalies, which appear to act as structural barriers that influence the spread of the aftershocks. Prominent low-Vanomalies from the lower crust to the mantle lithosphere are observed along the north-south trending rifts, suggesting the hot materials upwelling due to the tearing of the northward subducting Indian lithosphere. These observations may indicate that seismic velocity heterogeneity could play an essential role in earthquake initiation and the rupture process.

Geodynamic models of Indian continental flat slab subduction with implications for the topography of the Himalaya-Tibet region

The slab structure and high elevation of the Himalaya-Tibet region and their underlying mechanisms have been widely discussed. Many studies interpret a flat slab segment of Indian continental lithosphere located below the overriding plate, but interpretations of the northward extent of the flat slab differ substantially, with minimum estimates placing the boundary at the northern margin of the Himalaya (Indus-Yarlung Tsangpo suture), and maximum estimates placing it at the northern boundary of Tibet. In this study, we investigate for the first time if a flat slab segment of subducted buoyant Indian continental lithosphere below the Himalaya-Tibet region is geodynamically feasible and we quantify its northward extent, as well as its contribution to the high topography of the region. We conduct three large-scale fully-dynamic (buoyancy-driven) analogue experiments to simulate the subduction of the Indian continent. Our preferred, and geodynamically most feasible, model shows a continental flat slab extending northward up to ~ 320 km from the Himalayan thrust front, in agreement with recent estimates. Furthermore, it suggests that the positively buoyant flat slab segment of the Indian continent contributes some ~ 1.5-2 km to the high topography of the Himalaya-Southern Tibet region by providing an upward force to elevate the overriding Eurasian plate.

Southern Tibetan rifting since late Miocene enabled by basal shear of the underthrusting Indian lithosphere

Syncontractional extension is prominent in present-day Tibet, but its origin remains vigorously debated. Several deep-seated geodynamic processes (e.g., Indian underthrusting, horizontal flow, and mantle upwelling) have been linked to Tibetan rifting. Indian underthrusting is a good candidate because it can well explain why surface rifts are more prominent south of the Bangong-Nujiang suture; however, how Indian underthrusting causes extension is not well understood and lacks observational constraints. Seismic anisotropy, measured by exploiting the birefringence effect of shear waves, can be indicative of the deformation styles within the crust. Here, we unveil the dominant convergence-parallel alignment of anisotropic fabrics in the deep crust of the southern Tibetan rifts using seismic recordings collected from our recently deployed and existing seismic stations. This finding suggests that the strong north-directed shearing exerted by the underthrusting Indian plate is key to enabling present-day extension in southern Tibet.

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.

Modelling of crustal composition and Moho depths and their Implications toward seismogenesis in the Kumaon-Garhwal Himalaya

We image the lateral variations in the Moho depths and average crustal composition across the Kumaon–Garhwal (KG) Himalaya, through the H–K stacking of 1400 radial PRFs from 42 three-component broadband stations. The modelled Moho depth, average crustal Vp/Vs, and Poisson’s ratio estimates vary from 28.3 to 52.9 km, 1.59 to 2.13 and 0.17 to 0.36, respectively, in the KG Himalaya. We map three NS to NNE trending transverse zones of significant thinning of mafic crust, which are interspaced by zones of thickening of felsic crust. These mapped transverse zones bend toward the north to form a NE dipping zone of maximum changes in Moho depths, below the region between Munsiari and Vaikrita thrusts. The 1991 M_w6.6 Uttarakashi and 1999 M_w6.4 Chamoli earthquakes have occurred on the main Himalayan thrust (MHT), lying just above the mapped zone of maximum changes in Moho depths. Modelled large values of average crustal Vp/Vs (> 1.85) could be attributed to the high fluid (metamorphic fluids) pressure associated with the mid-crustal MHT. Additionally, the serpentinization of the lowermost crust resulted from the continent–continent Himalayan collision process could also contribute to the increase of the average crustal Vp/Vs ratio in the region.

Crustal deformation rates in Kashmir valley and adjoining regions from continuous GPS measurements from 2008 to 2019

We present GPS velocities in Kashmir valley and adjoining regions from continuous Global Positioning System (cGPS) network during 2008 to 2019. Results indicate total arc normal shortening rates of ~ 14 mm/year across this transect of Himalaya that is comparable to the rates of ~ 10 to 20 mm/year reported else-where in the 2500 km Himalaya Arc. For the first time in Himalayas, arc-parallel extension rate of ~ 7 mm/year was recorded in the Kashmir valley, pointing to oblique deformation. Inverse modeling of the contemporary deformation rates in Kashmir valley indicate oblique slip of ~ 16 mm/year along the decollement with locking depth of ~ 15 km and width of ~ 145 km. This result is consistent with the recorded micro-seismicity and low velocity layer at a depth of 12 to 16 km beneath the Kashmir valley obtained from collocated broadband seismic network. Geodetic strain rates are consistent with the dislocation model and micro-seismic activity, with high strain accumulation (~ 7e−08 maximum compression) to the north of Kashmir valley and south of Zanskar ranges. Assuming the stored energy was fully released during 1555 earthquake, high geodetic strain rate since then and observed micro-seismicity point to probable future large earthquakes of Mw ~ 7.7 in Kashmir seismic gap.

Localized foundering of Indian lower crust in the India-Tibet collision zone

The deep structure of the continental collision between India and Asia and whether India's lower crust is underplated beneath Tibet or subducted into the mantle remain controversial. It is also unknown whether the active normal faults that facilitate orogen-parallel extension of Tibetan upper crust continue into the lower crust and upper mantle. Our receiver-function images collected parallel to the India-Tibet collision zone show the 20-km-thick Indian lower crust that underplates Tibet at 88.5-92°E beneath the Yarlung-Zangbo suture is essentially absent in the vicinity of the Cona-Sangri and Pumqu-Xainza grabens, demonstrating a clear link between upper-crustal and lower-crustal thinning. Satellite gravity data that covary with the thickness of Indian lower crust are consistent with the lower crust being only ∼30% eclogitized so gravitationally stable. Deep earthquakes coincide with Moho offsets and with lateral thinning of the Indian lower crust near the bottom of the partially eclogitized Indian lower crust, suggesting the Indian lower crust is locally foundering or stoping into the mantle. Loss of Indian lower crust by these means implies gravitational instability that can result from localized rapid eclogitization enabled by dehydration reactions in weakly hydrous mafic granulites or by volatile-rich asthenospheric upwelling directly beneath the two grabens. We propose that two competing processes, plateau formation by underplating and continental loss by foundering or stoping, are simultaneously operating beneath the collision zone.

Seismological evidence for the earliest global subduction network at 2 Ga ago

The earliest evidence for subduction, which could have been localized, does not signify when plate tectonics became a global phenomenon. To test the antiquity of global subduction, we investigated Paleoproterozoic time, for which seismic evidence is available from multiple continents. We used a new high-density seismic array in North China to image the crustal structure that exhibits a dipping Moho bearing close resemblance to that of the modern Himalaya. The relict collisional zone is Paleoproterozoic in age and implies subduction operating at least as early as ~2 billion years (Ga) ago. Seismic evidence of subduction from six continents at this age is interpreted as the oldest evidence of global plate tectonics. The sutures identified can be linked in a plate network that resulted in the assembly of Nuna, likely Earth's first supercontinent. Global subduction by ~2 Ga ago can explain why secular planetary cooling was not appreciable until Proterozoic time.

Strain-Induced Fabric Transition of Chlorite and Implications for Seismic Anisotropy in Subduction Zones

Seismic anisotropy of S-wave, trench-parallel or trench-normal polarization direction of fast S-wave, has been observed in the fore-arc and back-arc regions of subduction zones. Lattice preferred orientation (LPO) of elastically anisotropic chlorite has been suggested as one of the major causes of seismic anisotropy in subduction zones. However, there are two different LPOs of chlorite reported based on the previous studies of natural chlorite peridotites, which can produce different expression of seismic anisotropy. The mechanism for causing the two different LPOs of chlorite is not known. Therefore, we conducted deformation experiments of chlorite peridotite under high pressure–temperature conditions (P = 0.5–2.5 GPa, T = 540–720 °C). We found that two different chlorite LPOs were developed depending on the magnitude of shear strain. The type-1 chlorite LPO is characterized by the [001] axes aligned subnormal to the shear plane, and the type-2 chlorite LPO is characterized by a girdle distribution of the [001] axes subnormal to the shear direction. The type-1 chlorite LPO developed under low shear strain (γ ≤ 3.1 ± 0.3), producing trench-parallel seismic anisotropy. The type-2 chlorite LPO developed under high shear strain (γ ≥ 5.1 ± 1.5), producing trench-normal seismic anisotropy. The anisotropy of S-wave velocity (AVs) of chlorite was very strong up to AVs = 48.7% so that anomalous seismic anisotropy in subduction zones can be influenced by the chlorite LPOs.

Probing crustal anisotropy by receiver functions at the deep continental drilling site KTB in Southern Germany

Seismic anisotropy is a unique observational tool for remotely studying deformation and stress within the Earth. Effects of anisotropy can be seen in seismic data; they are due to mineral alignment, fractures or layering. Seismic anisotropy is linked to local stress and strain, allowing modern geophysics to derive geomechanical properties from seismic data for supporting well planning and fracking. For unravelling anisotropic properties of the crust, the teleseismic receiver functions methodology has started to be widely applied recently due to its ability in retrieving the three-dimensional characteristics of the media sampled by the waves. The applicability of this technique is tested here by a field test carried out around the Kontinental Tiefbohrung site in southeastern Germany. We compare our results to previous investigations of the metamorphic rock pile of the Zone Erbendorf-Vohenstrauss, drilled down to 9 km depth, which sampled an alternating sequence of paragneiss and amphibolite, in which a strong foliation has been produced by ductile deformation. The application of the receiver functions reveals the presence of two distinct anisotropic layers within the metamorphic rock pile at 0-4 km and below 6 km depth, with up to 8% anisotropy; the depth of these two layers corresponds to the location of mica-rich paragneiss which show intense foliation, and finally proves the relation between the signal in the receiver functions, rock texture and presence of cracks. We have now the capability of providing insights from passive seismic data on geomechanical properties of the rocks, useful for geological exploration and engineering purposes, which will help influencing expensive drilling decisions thanks to future application of this seismic technique.

Lateral variation of the Main Himalayan Thrust controls the rupture length of the 2015 Gorkha earthquake in Nepal

The Himalaya orogenic belt produces frequent large earthquakes that affect population centers along a length of over 2500 km. The 2015 Gorkha, Nepal earthquake (M _w 7.8) ruptured the Main Himalayan Thrust (MHT) and allows direct measurements of the behavior of the continental collision zone. We study the MHT using seismic waveforms recorded by local stations that completely cover the aftershock zone. The MHT exhibits clear lateral variation along geologic strike, with the Lesser Himalayan ramp having moderate dip on the MHT beneath the mainshock area and a flatter and deeper MHT beneath the eastern end of the aftershock zone. East of the aftershock zone, seismic wave speed increases at MHT depths, perhaps due to subduction of an Indian basement ridge. A similar magnitude wave speed change occurs at the western end of the aftershock zone. These gross morphological structures of the MHT controlled the rupture length of the Gorkha earthquake.

Lower-crustal earthquakes in southern Tibet are linked to eclogitization of dry metastable granulite

Southern Tibet is the most active orogenic region on Earth where the Indian Plate thrusts under Eurasia, pushing the seismic discontinuity between the crust and the mantle to an unusual depth of ~80 km. Numerous earthquakes occur in the lower portion of this thickened continental crust, but the triggering mechanisms remain enigmatic. Here we show that dry granulite rocks, the dominant constituent of the subducted Indian crust, become brittle when deformed under conditions corresponding to the eclogite stability field. Microfractures propagate dynamically, producing acoustic emission, a laboratory analog of earthquakes, leading to macroscopic faults. Failed specimens are characterized by weak reaction bands consisting of nanometric products of the metamorphic reaction. Assisted by brittle intra-granular ruptures, the reaction bands develop into shear bands which self-organize to form macroscopic Riedel-like fault zones. These results provide a viable mechanism for deep seismicity with additional constraints on orogenic processes in Tibet.

The triggering mechanism of deep seismicity in Tibet remains unclear. Here the authors use experiments to show that granulite when deformed becomes brittle as it passes into the ecologite stability field developing macroscopic riedel fault zones thus providing an explanation for deep seismicity in Southern Tibet.

Nonuniform subduction of the Indian crust beneath the Himalayas

Himalayan tectonic activity is triggered by downward penetration of the Indian plate beneath the Asian plate. The subsurface geometry of this interaction has not been fully investigated. This study presents novel constraints on this geometry provided by two newly obtained, deep seismic reflection profiles. The profiles cover 100- and 60-km transects across the Yarlung-Zangbo suture of the Himalaya-Tibet orogen at c. 88°E. Both profiles show a crustal-scale outline of the subducting Indian crust. This outline clearly shows Indian understhrusting southern Tibet, but only to a limited degree. When combined with a third seismic reflection profile of the western Himalayas, the new profiles reveal progressive, eastward steepening and shortening in the horizontal advance of the subducting Indian crust.

Growth of the northeastern margin of the Tibetan Plateau by squeezing up of the crust at the boundaries

In classic orogenic models, the mountain range is underlain by a deep crustal root. Here we present the crustal and upper mantle structures along two receiver function profiles across Qilian, an orogen experiencing recent growth at the northern margin of the Tibetan plateau. Opposite to an expected crustal root beneath the orogen, the Moho beneath Qilian is arch-like, shallower beneath the center and deepens by up to 10 km beneath its southern and northern boundaries. Additional velocity interfaces sub-parallel to the Moho are observed in the lower crust of the basins south of Qilian, which we interpret as the top of a mechanically strong lower crust thrusting several tens of kilometers underneath Qilian. In the north, the small lateral offset between the surface and mantle traces of the thrust system reveals a steep boundary, indicating that the North China cratonic crust acts as a strong resistance to the northward growth of the plateau, forcing the development of the left-lateral strike-slip Haiyuan fault south of the northern Qilian suture. The young Qilian orogen thus has been rising and growing progressively from the boundaries to the center, squeezed up by more rigid tectonic blocks in the north and south.

Widespread ground motion distribution caused by rupture directivity during the 2015 Gorkha, Nepal earthquake

The ground motion and damage caused by the 2015 Gorkha, Nepal earthquake can be characterized by their widespread distributions to the east. Evidence from strong ground motions, regional acceleration duration, and teleseismic waveforms indicate that rupture directivity contributed significantly to these distributions. This phenomenon has been thought to occur only if a strike-slip or dip-slip rupture propagates to a site in the along-strike or updip direction, respectively. However, even though the earthquake was a dip-slip faulting event and its source fault strike was nearly eastward, evidence for rupture directivity is found in the eastward direction. Here, we explore the reasons for this apparent inconsistency by performing a joint source inversion of seismic and geodetic datasets, and conducting ground motion simulations. The results indicate that the earthquake occurred on the underthrusting Indian lithosphere, with a low dip angle, and that the fault rupture propagated in the along-strike direction at a velocity just slightly below the S-wave velocity. This low dip angle and fast rupture velocity produced rupture directivity in the along-strike direction, which caused widespread ground motion distribution and significant damage extending far eastwards, from central Nepal to Mount Everest.

GEOPHYSICS. Layered deformation in the Taiwan orogen

The underthrusting of continental crust during mountain building is an issue of debate for orogens at convergent continental margins. We report three-dimensional seismic anisotropic tomography of Taiwan that shows a nearly 90° rotation of anisotropic fabrics across a 10- to 20-kilometer depth, consistent with the presence of two layers of deformation. The upper crust is dominated by collision-related compressional deformation, whereas the lower crust of Taiwan, mostly the crust of the subducted Eurasian plate, is dominated by convergence-parallel shear deformation. We interpret this lower crustal shearing as driven by the continuous sinking of the Eurasian mantle lithosphere when the surface of the subducted plate is coupled with the orogen. The two-layer deformation clearly defines the role of subduction in the formation of the Taiwan mountain belt.

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.

Great Himalayan earthquakes and the Tibetan plateau

It has been assumed that Himalayan earthquakes are driven by the release of compressional strain accumulating close to the Greater Himalaya. However, elastic models of the Indo-Asian collision using recently imaged subsurface interface geometries suggest that a substantial fraction of the southernmost 500 kilometres of the Tibetan plateau participates in driving great ruptures. We show here that this Tibetan reservoir of elastic strain energy is drained in proportion to Himalayan rupture length, and that the consequent growth of slip and magnitude with rupture area, when compared to data from recent earthquakes, can be used to infer a approximately 500-year renewal time for these events. The elastic models also illuminate two puzzling features of plate boundary seismicity: how great earthquakes can re-rupture regions that have already ruptured in recent smaller earthquakes and how mega-earthquakes with greater than 20 metres slip may occur at millennia-long intervals, driven by residual strain following many centuries of smaller earthquakes.