Rheology Of Lithosphere Research Articles

The Effective Elastic Thickness of the Lithosphere Reveals Tectonic Process and Lithospheric Rheology of the Eurasian Basin

ABSTRACTWe obtained the effective elastic thickness of the lithosphere (Te) in the Eurasian Basin using the fan‐shaped wavelet method. For ultraslow spreading ridges, the magmatic heat input is insufficient to form continuous quasisteady‐state magma lenses along the axis, resulting in significant variability in the rheological strength of the lithosphere. We suggested that the tectonic inheritance resulted in the different lithospheric rheology of WVZ and SMZ, consequently leading to distinct tectonic‐magmatic activities. The SMZ with less magmatic activity has a weaker lithosphere than that of the EVZ with extensive magmatism. We propose that the lithospheric rheology and thermal structure of ultraslow spreading ridge may be influenced not by the amount of magma supply but by the mode of magma emplacement.

A numerical study of contemporary crustal deformation partitioning across the Southwestern Tian Shan orogen

The Tian Shan region is a typical example of crustal deformation within an intracontinental environment, where the lithospheric rheological properties are marked by significant spatial heterogeneity. However, the role of lithospheric rheology in crustal strain partitioning within the interior of the Tian Shan orogenic belt remains nebulous. Here, we utilized a two-dimensional viscoplastic model with contact elements constrained by GPS velocities and fault slip rates to investigate the influence of the block's strength on crustal deformation in the southwestern Tian Shan region. Our founding suggests that a weaker lower crust beneath the northern Tian Shan region offers a potential mechanism for the contemporary deformation patterns. Contemporary crustal deformation at Tian Shan is mainly concentrated along the piedmont thrust-and-fold belts and diffuses in the orogen's interior. This pattern of crustal strain partitioning is closely linked to the interplay between fault slips and lithospheric deformation. This finding emphasizes the significant role that the lithospheric rheology and pre-existing faults played in the deformation partitioning in the interior of the Tian Shan orogeny.

Present-day stress field, strain rate field and seismicity of the Pannonian region: overview and integrated analysis

We present an overview of recent findings on the seismicity, stress field and strain rate field of the Pannonian region. We furthermore show new inferences for deformation mechanisms and lithospheric rheology via an integrated analysis of the datasets. The N-NW motion of the Dinarides and the opposite, SW motion of the E-and S-Carpathians induce shortening in the western and central Pannonian Basin while leading to regional shearing around the basin's SE boundary. This induces moderate seismic activity related to strike-slip and reverse faulting. Maximum horizontal stress and geodetic shortening directions generally agree, implying that upper crustal stresses and surface deformation correspond to the same forces. 2D rheological models confirm that the shallow upper crust is the only brittle layer in the Pannonian lithosphere, while the brittle-ductile transition zone could be as shallow as 6-9.5 km in the Danube Basin. Comparison of seismic moment rates and moment rate predictions from geodetic strain rates show that deep Pannonian sub-basins accumulate major seismic deficits. This could be explained by dominantly aseismic deformation of the weak upper crust, as supported by a comparison with case studies from different geodynamic environments. Stronger mountainous regions in the study area show no to moderate seismic deficits.

Spatial variations in the effective elastic thickness of the Indian Ocean lithosphere

The Indian Ocean lithosphere is a complex assemblage of large igneous provinces, seamounts, plateaus and ridges of different loading ages and tectono-thermal evolution. As a proxy for the strength of tectonic plates, effective elastic thickness (Te) illustrates the relationship between surface deformation and lithospheric rheology of the diverse provinces in response to long-term tectonic processes. Mapping the spatial variations in lithospheric rheology can aid in understanding the detailed tectono-thermal history of the Indian Ocean. In this paper, we perform an assessment of the spatial variation of Te for the Indian Ocean from the inversion of the real free-air admittance between free-air gravity anomalies and bathymetry corrected for the effect of density variations within sediments using a continuous wavelet spectral analysis. Incorporating the effect of sediments substantially reduces Te estimates and better corresponds with the tectonic units in the study region. The results show low overall Te over the Indian Ocean attributed to magmatism and temperature during a multistage opening process. We further demonstrate that temperature controls the strength of warm and young oceanic lithosphere, evidenced by the positive correlation between Te and geothermal proxies. Finally, moderately low Te values at the Southwest Indian Ridge suggest a relatively cold ultraslow lithosphere with sparse magmatism compared to typical mid-ocean ridges.

Lithospheric rheology in West Tibet and Pamir constrained from the postseismic deformation of the 2005 Mw7.6 Kashmir (Pakistan) earthquake

We study the postseismic deformation of the 2005 Mw7.6 Kashmir earthquake in Pakistan, with a combination of the near- and far-field geodetic observations from the Kashmir Himalaya and the Tibetan Plateau, respectively. We construct a 3D finite-element model that integrates topographic relief, curved fault plane, and layered lithospheric structures to analyze stress-driven afterslip and estimate the rheological parameters of the Tibetan lithosphere. Our findings reveal that the ten-year afterslips of up to 30 cm were primarily localized at the downdip edges of co-seismic slip patches, causing near-field surface deformation of up to ∼5 cm within the first year. In contrast, viscoelastic relaxation in the lower crust of western Tibet and Pamir was required to reproduce far-field displacements of up to 1–2 cm recorded by GPS during the first 7 years. The model that best fits to the GPS displacements suggests a Maxwell viscosity exceeding 1020 Pa s for the Indian lower crust, while a much lower viscosity of ∼4–10×1018 Pa s is inferred for the lower crust beneath Pamir and western Tibet, indicating a softer Tibet and Pamir in response to the continental collision between India and Eurasia.

Spatiotemporal dominance of afterslip and viscoelastic relaxation revealed by four decades of post-1973 Luhuo earthquake observations

Measuring surface displacements driven by afterslip and viscoelastic relaxation following large earthquakes enables inferring frictional and rheological properties of Earth's outermost layers where hazardous earthquakes occur. However, the two concurrent mechanisms generate similar and intertwined surface deformations, complicating the extraction of physical information from the measurements. Here, we quantify the spatiotemporal dominance of co-evolving afterslip and viscoelastic relaxation following the 1973 Ms 7.6 Luhuo earthquake in eastern Tibet, using 42 years (1976‒2018) of fault-crossing short-baseline observations. The finite-fault slip of this M7+ earthquake was constructed from triangulation data measured in 1961‒1975. Based on the model incorporating two postseismic relaxation mechanisms and calibrated by the measurements over four decades, we show that, in the temporal domain, afterslip may produce geodetically measurable deformations (>1.5 mm/yr) in local near-field regions even 20 years after the mainshock, with spatially broader impact within ∼10 years. Viscoelastic relaxation may last for nearly six decades, imposing a broad influence for three decades. In the spatial domain, afterslip may produce deformations within ∼2 times the seismogenic depth (SD; 20 km) from the fault, but dominantly in the near-field of ∼1 times the SD. In contrast, viscoelastic relaxation may affect a much wider region reaching 200 km. While afterslip and viscoelastic relaxation collectively affected the region ∼2 times the SD from the fault in the early postseismic period, their resulting deformations gradually separated in space with time, with afterslip-induced deformation shrinking toward the fault, and surface deformation caused by viscoelastic relaxation concentrating in the mid-field, ∼2–3 times the SD from the fault. Their spatiotemporal partitioning helps better learn fault frictional properties and lithospheric rheology based on geodetic data acquired from different domains.

4D fault evolution revealed by footwall exhumation modelling: A natural experiment in the Malawi rift

The evolution of normal fault arrays during rift extension reflects paleo-plate boundary conditions and lithospheric rheology, while controlling seismic hazard and the distribution of basin-hosted resources. Yet, constraining their spatiotemporal development is challenging, particularly where geophysical and subsurface data are absent. Here, we test footwall exhumation modelling using thermochronology as a means of elucidating 4D normal fault array evolution, using the Miocene Central Basin of the Malawi Rift as a natural laboratory. Along-strike trends in exhumational cooling recorded by vertical transects of apatite fission-track and (U–Th)/He data from the basin-bounding Usisya fault scarp reveal a diachronous footwall uplift history that closely reflects 4D trends in hangingwall subsidence recorded by previously published seismic and well data. Initially, pronounced footwall exhumation is restricted to the centres of a series of four isolated normal faults, mirroring the distribution of early syn-rift depocentres. The later onset of footwall exhumation in the intervening areas marks subsequent fault segment propagation and linkage as they formed the through-going Usisya fault system. Elsewhere, cumulative exhumation recorded in the Usisya footwall remains low, coinciding with more significant intra-basinal faulting. This study shows that footwall exhumation modelling constrained by thermochronologic data can reveal the spatiotemporal evolution and strain partitioning within normal fault arrays.

Physical mechanism-driven extraction model for post-seismic deformation recorded by continuous GNSS observations

Post-seismic deformation caused by large earthquakes is one of the most effective parameters for investigating the frictional behavior of faults and the rheological characteristics of the lithosphere. This study introduces a modified functional model that accounts for both afterslip and viscoelastic relaxation effects, incorporates a physical mechanism, and effectively disentangles distinct mechanisms of post-seismic deformation captured by continuous GNSS observations.A parameter optimization solution method based on the Bayesian framework was adopted to obtain optimal estimates and corresponding confidence intervals. The effectiveness of proposed model was validated through simulation experiments. The proposed method was applied to study post-seismic deformations from the 2015 Nepal MW 7.8 earthquake, and the different characteristics of post-seismic deformations caused by different mechanisms were discovered at different stations. The afterslip effect has the duration of 2.67 years, whereas the viscoelastic relaxation effect will continue for 150.8 years. Moreover, the calculated equivalent viscosity coefficient in the southern Tibet region was 1.1 ± 0.3 × 1018 Pa·s. Our proposed model can more accurately and objectively separate different mechanisms from the post-seismic displacements observed by GNSS than previous one, and thus improve the accuracy of fault kinematics and regional lithosphere rheology investigations.

The importance of being molten: 100 Myr of synmagmatic shearing in the Yilgarn Craton (Western Australia). Implications for mineral systems

Between the Neoarchean and the early Proterozoic, plate tectonics was gradually established on our planet. How this transition took place, and what kind of geodynamic paradigm dominated before, is a matter of debate, and systematic multidisciplinary studies help improving our understanding of such a transitional period of Earth’s evolution.The Yilgarn Craton of Western Australia exposes upper-crustal granite–greenstone terranes that are juxtaposed along large-scale, Neoarchean shear zones, which mainly developed during the 2730 – 2660 Ma Yilgarn Orogeny. By combining structural, microstructural and geochemical data, we present here a systematic study of the main shear zones occurring in the best-exposed, northwestern half of the craton. Our structural data demonstrate that large granitic magma sheets were emplaced during shearing. Although these shear zones exhibit a variety of fabrics that developed over a wide range of temperatures, all preserve memory of at least one stage of shearing at near-solidus temperatures. Our results allow tracing a c. 100 Myr-long tectonomagmatic evolution of this portion of the Archean crust, with implications for lithosphere rheology, crustal evolution, and mineral systems. The major synorogenic structures can be interpreted as inclined transpressional shear zones, along which the bulk of the 3D deformation was efficiently partitioned between greenstones and crystallizing granitic sheets within the shear zones, with the latter accommodating large amounts of orogen‐parallel strain through suprasolidus viscous flow. The positive feedback between protracted magmatism and transpression promoted the extraction and upward transfer of syntectonic magma originated between the uppermost mantle and the mid-crust, as demonstrated by our geochemical dataset. This synkinematic mode of granitic magma transfer contrasts fundamentally to the diapiric mode that was dominant during the lithospheric extension phase that pre-dated the Yilgarn orogeny. Our study suggests that the shear zone system studied here likely played a major role in controlling magma/fluid pathways throughout the late Archean lithosphere, therefore playing a critical role in controlling the development of near-surface mineral systems.

Integrated Investigation on Heterogeneous Lower Crust Rheology in Kyushu and Afterslip Behavior Following the 2016 Mw7.1 Kumamoto Earthquake

AbstractThe viscoelastic lower crust beneath Kyushu Island, influenced by the volcanic arc, interplays with active crustal faults in this region and helps to shape local tectonics. In this study, we employed a three‐dimensional viscoelastic finite element model to gain insights into the lithospheric rheology and crustal faulting kinematics, through modeling the postseismic deformation processes of the 2016 Mw 7.1 Kumamoto earthquake. Our model reveals a viscosity of 2 × 1020 Pa s for the lower crust and 2 × 1019 Pa s for the upper mantle. A reduced lower crust viscosity of 2 × 1019 Pa s in the volcanic arc area is required for better reproducing the Global Positioning System data. The stress‐driven afterslip decays rapidly over time and is up to 0.3 m within 5 years after the earthquake. We propose additional normal‐component afterslip to better explain the complex postseismic deformation in the near field, which may be due to the interaction between the fault and volcano Aso.

Illite K‐Ar Dating of the Leibo Fault Zone, Southeastern Margin of the Tibetan Plateau: Implications for the Quasi‐Synchronous Far‐Field Tectonic Response to the India‐Asia Collision

AbstractWhether tectonic strain from the early stage India‐Asia collision has synchronously affected the far‐field margin of the Tibetan Plateau is crucial for understanding plateau deformation and growth processes. However, direct evidence for early far‐field deformation remains scarce. Utilizing illite K‐Ar dating of three fault gouge samples, we established the faulting history of the Leibo fault zone (LFZ) at the southeastern margin of the Tibetan Plateau (SEMTP). Consistent authigenic illite ages of 52 ± 2, 54 ± 12 and 55 ± 6 Ma suggest the reactivated thrust faulting of the LFZ in the Early Cenozoic. Positioned ∼700 km east of the collisional boundary and at the intersection of three blocks with distinct lithospheric rheology in strength/viscosity, this event suggests a quasi‐synchronous far‐field tectonic response in the SEMTP to the India‐Asia collision.

Lithospheric Rheology and Crustal Deformation Across the Northeastern Tibet and Their Implications for Plateau Growth

AbstractUnderstanding lithospheric rheology is crucial in investigating tectonic evolution of intra‐continental tectonic boundary. Here, we use geodetic observations to infer lithospheric rheology across the northeastern Tibet based on a 2D viscoelastic model. Our findings reveal a lower‐crust viscosity of <1022 Pa·s underneath its margins, lower than those estimated underneath its vicinities. By comparing deformation patterns and lithospheric rheology here with those observed in the eastern Tibet, we propose that lateral variations in lower‐crust viscosity control deformation patterns and topographic gradients along the Tibet margins. The presence of low viscosity lower‐crust can lead to the development of contrasting topographic gradients and shape the plateau's geomorphology and deformation characteristics during outward growth of the Tibet. We here emphasize the subtle variations in the lower‐crust rheology between deforming blocks and the corresponding mountain ranges, which play an important role in orogeny along the intracontinental convergence boundary.

Lithosphere Removal in the Sierra Nevada de Santa Marta, Colombia

AbstractThe Sierra Nevada de Santa Marta (SNSM) in northwestern Colombia is one of the world's highest coastal mountains, with an elevation above 5 km. Gravity measurements show that the SNSM has a positive Bouguer anomaly (>+130 mGal), indicating that the mountain lacks a crustal root. In this work, we test the hypothesis that these observations can be explained by gravitational removal of the dense lower lithosphere. We use 2D numerical models to examine the dynamics of lithosphere removal and its effect on surface elevation and gravity anomaly. The models consist of continental lithosphere that includes a pre‐thickened crustal region, representing the SNSM. In our preferred model, the dense mantle lithosphere and crustal root are gravitationally unstable and undergo removal as local drips within ∼10 Ma from the onset of foundering. This creates an area of thinned crust (∼38 km) underlain by a buoyant sublithospheric mantle where melting and low seismic velocities are predicted. Subsequent non‐isostatic forces maintain a topography of 3.3 km with a 2D Bouguer gravity anomaly of +103 mGal. Parameter tests show that a strong lower‐crustal rheology provides greater support for the high topography and that a weak mantle lithosphere rheology produces faster removal. The models demonstrate that local lithosphere dynamics can explain the first‐order observations in the SNSM. We propose that lithosphere removal could have occurred recently (∼2 Ma), explaining the localized low seismic velocity zone below the SNSM, or at 56–40 Ma, inducing anomalous short‐lived Eocene magmatism.

Lower Crust Viscosity in Central Tibet Inferred From InSAR Derived Deformation Around Siling Co Lake After Its Rapid Expansion in the 2000s

AbstractConstraining the lithospheric rheology of Tibetan Plateau is important for the physical understanding of its tectonics. Siling Co, the largest high‐altitude lake in the world, has experienced a rapid water level increase in the 2000s. The resulting loading changes stimulate the viscoelastic response of the lower crust, giving access to study the lithospheric rheology in central Tibet. Here, we derive a clear subsidence signal around the lake with peak velocity of 4 mm/yr from two tracks of Sentinel‐1 InSAR images acquired from 2017 to 2022. Our viscoelastic modeling suggests a ductile layer 15 km beneath the surface with a decadal‐scale, steady‐state viscosity of 1–4 × 1019 Pa s. This value is consistent with the viscosity inferred from millennium‐scale shoreline changes, but is about 10 times higher than the viscosities derived from the deformation surrounding Siling Co before 2011, highlighting the role of transient viscosity in controlling the surface deformation.

Rheology of continental lithosphere and seismic anisotropy

Rheology of continental lithosphere and seismic anisotropy

Middle Miocene Onset of the Litang Fault System Records Kinematic Change in Eastern Tibet

AbstractThe ∼400‐km‐long Litang fault system (LTFS) is a major intracontinental strike‐slip fault inside the Chuandian block, eastern Tibet, but its evolution and role in accommodating the India‐Asia convergence remain poorly known. Structural analysis shows that the LTFS splits into 5 strands as a left‐lateral, right‐stepping en‐echelon pattern formed under NW‐directed compression, subsequently reactivated by transtensive faults under NNE‐directed extension. Displaced geological and morphological markers yield a cumulative left‐lateral offset of 28.9–42.8 km. Inverse thermal‐history modeling of thermochronological data of the faulted rocks reveal accelerated cooling at 38–35 Ma, 16–13 Ma, and 7–5 Ma. The late Eocene rapid cooling is ascribed to the reactivation of the Garze‐Litang suture. Rapid cooling events at 16–13 Ma and 7–5 Ma record the onset of transpression and transtension of the LTFS, respectively, yielding a geologic slip rate of 2.6 ± 0.7 mm/yr. Both bifurcated geometry and slow slip rate of the LTFS since 16–13 Ma indicate diffuse deformation inside the Chuandian block, contrasting with strain localized on fast‐slip strike‐slip faults on the block margins. This implies a significant kinematic transition in the middle Miocene, such that the extrusion of the segmented mega‐blocks has been accommodated by both localized and distributed deformation in eastern Tibet. This tectonic transition could be explained by a change in lithospheric rheology from an earlier rigid state to a viscous state underneath the Chuandian block due to thermal weakening of the lower crust. We thus reconcile the end‐member geodynamic models of block extrusion and lower crustal flow in late Cenozoic times.

SMatStack to Enhance Noisy Teleseismic Seismic Phases: Validation and Application to Resolving Depths of Oceanic Transform Earthquakes

AbstractThe depths of most moderate‐sized oceanic earthquakes are poorly constrained because of the lack of local recording stations and noisy teleseismic recordings. This hampers our understanding of slip behaviors along oceanic faults and the mechanical properties of the oceanic lithosphere. In this study, we develop a new method to enhance the weak body‐wave signals, particularly the depth phases, associated with earthquakes on oceanic transform faults using large‐aperture arrays in the teleseismic range. We simulate the enhanced teleseismic signals to refine the centroid depths of moderate‐sized earthquakes. We validate the new approach using synthetic waveforms and show it outperforms conventional methods when dealing with noisy signals. We obtain the depth estimates for three moderate‐sized earthquakes on the Chain transform fault in the equatorial Atlantic Ocean and find two of them are consistent with a local catalog derived using oceanic bottom seismometers. Application of the new method to the past decades of teleseismic recordings of moderate‐sized earthquakes on the large and slow‐slipping transform faults will provide significantly improved constraints on the width of the seismogenic zone, thus advancing our understanding of the rheology of oceanic lithosphere and earthquake processes in oceanic settings and, by comparison, their more dangerous continental counterparts. The new method is not limited to oceanic transform earthquakes, but can be easily adapted to other seismological studies in which noisy but coherent signals could be revisited for better usage.

Weakness of the Indian Lower Crust Beneath the Himalaya Inferred From Postseismic Deformation of the 2015 Mw 7.8 Gorkha Earthquake

AbstractThe 2015 Mw 7.8 Gorkha (Nepal) earthquake induced prolonged postseismic deformation extending northward beyond the Yarlung Zangbo Suture, which provides unique opportunities to better understand the lithospheric rheology in Himalaya and southern Tibet. Here, we used the first 5‐year Global Positioning System observations to study the main postseismic processes following this event, including viscoelastic relaxation and afterslip, based on a three‐dimensional finite element model. We considered a realistic geometry of the underthrusting Indian plate according to various geophysical images. We found that the models with a uniform elastic Indian lower crust fail to fit the vertical displacements. A heterogeneous Indian lower crust with the transition from elastic (high‐viscosity) to low‐viscosity approximately under the Main Central Thrust is required to reproduce the observed postseismic uplift between China‐Nepal border and Peiku Lake, indicating the weakness of the Indian lower crust from the Lesser to High Himalaya. The afterslip simulated using a weak shear zone takes place in the adjacent area downdip of the rupture zone. The preferred model suggested that viscosities of the Tibetan lower crust, weakened Indian lower crust, and shear zone are 3 × 1018, 1019, and 4 × 1018 Pa s, respectively. The viscosity of the underthrusting Indian upper mantle was roughly estimated to be greater than 1021 Pa s. The model results imply that the near‐field deformation is dominated by both afterslip and viscoelastic relaxation of the weakened Indian lower crust, not only afterslip as suggested by previous studies.

Two-phase kinematic evolution of the Qilian Shan, northern Tibetan Plateau: Initial Eocene−Oligocene deformation that accelerated in the mid-Miocene

The Cenozoic growth of the Tibetan Plateau and the distribution of deformation across it are a consequence of India-Asia collision and continued convergence, which have implications for studies of continental tectonics. The spatio-temporal development of Cenozoic deformation along the northern margin of the plateau is an important issue that can be better understood by testing various models of plateau growth. The northern Tibetan Plateau is bounded by the Cenozoic Qilian Shan thrust belt and the Haiyuan left-slip fault. We conducted geologic mapping, field observations, electron spin resonance (ESR) dating, and apatite (U-Th)/He (AHe) and apatite fission-track (AFT) analysis in the Qilian Shan thrust belt to improve our understanding of the timing of brittle faulting and range exhumation in the northern Tibetan Plateau. We document the first direct age constraints for Oligocene deformation within the central Qilian Shan via ESR dating, which correlates with AHe-AFT cooling ages in adjacent ranges. We demonstrate that the Qilian Shan thrust belt experienced a two-phase growth history, including Eocene−Oligocene fault-related uplift shortly after the India-Asia convergence, and mid-Miocene regional overprinting deformation that reactivated the proximal thrust faults. This deformational pattern suggests that the Qilian Shan thrust belt has experienced out-of-sequence development since the Eocene−Oligocene and has persisted as the stationary northeastern boundary of the Himalayan-Tibetan Orogen throughout the Cenozoic. The Paleozoic Qilian suture systems acted as a pre-existing weakness and played a decisive role in controlling the lithospheric rheology, which therefore impacted the timing, pattern, and strain distribution of Cenozoic deformation across the northern Tibetan Plateau.

Post-seismic deformation of the 2008 Wenchuan earthquake reveals a misaligned rheological boundary in the lower crust along the eastern Tibetan Plateau margin

SUMMARY Along the margins of orogenic plateaus, the viscous Earth structure and fault geometries play a primary role in controlling the tectonic evolution and earthquake generation. After the 2008 Mw 7.9 Wenchuan earthquake, the long-standing debate regarding the tectonics producing and maintaining prominent topography across the Longmen Shan reignited. Post-seismic deformation, representing the surface strain history in response to lithospheric stress perturbations, provides important insights into the lithospheric rheology and active structures. Here, we construct a new 3-D post-seismic deformation model for the Wenchuan earthquake, invoking viscoelastic relaxation and afterslip. Our best-fitting model indicates that the steady-state viscosities of the lower crust in the region to the immediately west of the Songpan-Ganzi terrane and beneath the Songpan-Ganzi terrane are estimated to be 4.0 × 1018 and 1.0 × 1018 Pa s, respectively. Our results, combining geophysical and geodetic observations and model analyses, highlight the prevalent parallelism between the rheological and structural boundaries of the lower crust, which diverge northward away from the trend of the Longmen Shan fault at ∼20°. This diverging rheological structure and the partially coupled upper and lower crust have broad implications for the stress build-up, strain partitioning and deformation styles along the eastern Tibetan Plateau margin.