Longmenshan Fault Research Articles

Assessment of Emergency Water Sources Using Electrical Resistivity Tomography: A Case Study in the Longmen Shan Fault Zone

The Longmen Shan fold and thrust belt, situated on the eastern margin of the Qinghai–Tibet Plateau, is prone to disasters like earthquakes and debris flows. Thus, applying rapid assessment methods for emergency water sources in disaster-affected areas is crucial for local populations and ensuring an effective response to disasters. In this study, we employed electrical resistivity tomography (ERT) to investigate groundwater resources in post-disaster regions. By integrating the results from ERT profiles with geological and borehole information, we determined the lithology and depth of aquifers. Additionally, we analyzed groundwater recharge and discharge patterns relative to surface water during various precipitation periods and generated a hydrogeological profile for the region. Borehole information confirmed our inferred lithology and aquifer depth, thereby ensuring a reliable water supply during emergencies. This study demonstrates the feasibility of ERT for rapidly identifying water resources in geologically complex environments, providing a scientific foundation for water resource management in disaster-prone regions.

Study on the attenuation relationship of the acceleration envelope parameters for the Wenchuan earthquake aftershocks

The seismic attenuation relationship between ground motion parameters (such as peak acceleration and response spectra value) and seismic parameters (such as magnitude and epicentral distance) is an important foundation for seismic hazard analysis and the core of determining seismic input parameters for seismic resistance engineering. The acceleration envelope parameters, which describe the relationship between ground motion intensity and time variation, are primarily used for artificially synthesizing seismic motion. At present, little research has been performed on the attenuation relationship of acceleration envelope parameters in the Longmenshan fault zone on the eastern side of the Tibetan Plateau in China. Therefore, this study selected the Mw4-6 aftershock records of the 2008 Wenchuan earthquake that occurred on the Longmenshan fault zone and established a commonly used three segment envelope model parameter attenuation relationship. We classified aftershock records based on their source mechanisms and obtained attenuation relationship models for thrust slip aftershocks, thrust and strike slip aftershocks, and strike slip aftershocks. The results are as follows: (1) The thrust slip aftershock records had the longest rising stage which is the time difference from the arrival of P-waves to the beginning of the stable sustained stage of seismic motion. Thrust and strike slip aftershocks records had the longest stable sustained stage period and the slowest attenuation at the tail of the record. (2) The attenuation relationship of the acceleration envelope parameters commonly used in Chinese engineering for artificially synthesizing seismic motion is based on the strong earthquake records in the western United States. But, compared to the attenuation characteristics of the acceleration envelope function in the western United States, the Mw4-6 earthquake records on the Longmenshan fault zone had a slower attenuation rate at the tail of the record. So, accurately artificially synthesizing seismic motion through envelope parameter attenuation model requires the use of attenuation model established by earthquake records in this region.

Dynamic changes of gravity field before the Luding MS6.8 earthquake and its crustal material migration characteristics

On September 5, 2022, an earthquake of magnitude MS6.8 occurred in Luding County, Sichuan Province. This earthquake occurred at the key part of the southeast-clockwise extrusion of material on the eastern margin of the Qinghai–Xizang Plateau, the Y-shaped confluence of the Xianshuihe, Longmenshan and Anninghe fault zones. In this study, the three-dimensional dynamic crustal density changes in the earthquake area are obtained by the typical gravity change data from 2019 to 2022 before the earthquake and gravity inversion by growing bodies. The results indicate that gravity changes presented an obvious four-quadrant and gradient belt distribution in the Luding area before the earthquake. The three-dimensional density horizontal slices show that small density changes occurred at the epicenter in the mid-to-upper crust between 2019.9 - 2020.9 and 2019.9–2021.9. At the same time, the surrounding areas exhibited a positive and negative quadrant distribution. These observations indicate that the source region was likely in a stable locked state, with locking-in shear forces oriented in the NW and NE directions. From 2021.9 to 2022.8, the epicentral region showed negative density changes, indicating that the source region was in the expansion stage, approaching a near-seismic state. The three-dimensional density vertical slices reveal a southeastward migration of positive and negative densities near the epicenter and on the western of the Xianshuihe Fault Zone, indicating that the material is flowing out to the southeast. The observed local negative density changes at the epicenter along the Longmenshan Fault Zone are likely associated with the NE-oriented extensional stress shown by the seismic source mechanism. The above results can provide a basis for interpreting pre-earthquake gravity and density changes, thereby contributing to the advancement of earthquake precursor theory.

Shallow crustal structure of the northern Longmen Shan fault zone revealed by a dense seismic array with ambient noise analysis

The Longmen Shan (LMS) fault zone located in the southwest of China is not only the boundary tectonic zone of the Tibetan Plateau and South China block, but also an important part of the north–south seismic belt of China. To investigate the detailed crustal structure of the LMS fault zone, we deployed a dense seismic array of 151 short-period temporary stations in its northern section from March 16 to April 6, 2023. Continuous vertical component ambient noise waveforms recorded by these stations were cross-correlated for all station pairs, enabling the extraction of Rayleigh wave phase velocity dispersion curves across periods ranging from 0.5 s to 6 s. We then inverted these data to derive the crustal shear wave velocity using direct surface wave tomography. Meanwhile, a linear sub-array consisting of 54 stations and crossing the surface rupture of the 2008 Wenchuan earthquake was used to image the fault structure by applying the Transmitted Surface Wave Reverse Time Migration method (TSW-RTM). The ambient noise tomography results show that strong lateral heterogeneity exists in the shallow crust in the study area, with a high shear velocity in the northwest and a low shear velocity in the southeast. The low shear velocity is generally distributed between the Yingxiu-Beichuan fault (YBF) and the Guanxian-Jiangyou fault (GJF). The results of TSW-RTM indicate that the Wenchuan-Maoxian fault (WMF) and the YBF are both characterized with a steep dip to the northwest. Our results could serve as the important basis for the study of segmentation characteristics of the LMS fault and seismic hazard preparation and mitigation.

Stress heterogeneity in the eastern Tibetan Plateau and implications for the present-day plateau expansion

The eastward expansion of the Tibetan Plateau has resulted in different earthquake types in the eastern Tibetan Plateau, but the mechanism remains unclear. Here, we construct a three-dimensional visco-elastoplastic finite element model considering the topography to investigate the influence of fault geometry and rheological heterogeneity on stress fields. In our best-fitting model, the minimum principal stress is nearly vertical around the southern Huya fault zone, which is adjacent to the Longmen Shan fault zone, due to the significant mid-lower crust lateral rheological heterogeneity, and the thrust stress regime accounts for the reverse fault and thrust-dominated earthquakes. In this scenario, the eastward horizontal motion of the mid-lower crust is obstructed and facilitates thrust faulting, suggesting the limited eastward expansion of the Tibetan Plateau. In contrast, the northern Huya fault zone, one of the terminal branches of the East Kunlun fault, accommodates the continuous eastward extrusion of the East Kunlun fault, where the stress regime under a more homogenized crust favors the strike-slip faulting process, along with the dominant strike-slip earthquakes. Moreover, the best-fitting of stress regime explains the thrust-dominated 2008 Ms. 8.0 Wenchuan and 2013 Ms. 7.0 Lushan earthquakes on the Longmen Shan fault zone. Combining geophysical and geodetic observations and model analyses, we propose that the hybrid deformation mode in the eastern Tibetan Plateau is accommodated by upper crustal shear and thrusting deformation and mid-lower crustal thickening driven by the gravitational potential energy gradient. Our results elucidate the mechanism for differences in strong historical earthquakes and, more importantly, isolate the effect of fault geometry from those of heterogeneous viscosity on crustal deformation and stress heterogeneity in the eastern Tibetan Plateau.

Stress modeling for the upper and lower crust along the Anninghe, Xianshuihe, and Longmenshan Faults in southeastern Tibetan plateau

Earthquake occurrence depth in the crust is related to stress, temperature, and brittle–ductile transition, which is also near the transition depth of the upper to lower crust. The composition variation between the upper and lower crust causes remarkable changes of rheological properties and variation in stress distribution. Clarifying the detailed stress distribution in the upper and lower crust is crucial for understanding the brittle–ductile transition and the stress environment of the seismogenic zone. The Southeastern Tibetan Plateau (SETP), with wide spread of active strike−slip faults and clustered earthquakes, provides a natural field for investigating the relationships between crustal stresses, deformation behaviors, and earthquake mechanics. By considering the rheological properties of granite and anorthite, this paper established stress models with different boundary depths (15, 20 and 25 km) between the upper and lower crust along the Anninghe, Xianshuihe, and Longmenshan Faults in the SETP with a horizontal strain of 6 × 10−4 extracted from in situ stress data. The stress model with different geothermal gradients and a boundary depth of 20 km between the upper and lower crust suggests two distinct types of the brittle–ductile transition below these three faults. Simultaneously, the stress model can account for the continuity of earthquake depth distribution below the Longmenshan Fault and the seismic gap below the Anninghe and Xianshuihe Faults. The continuity of earthquake depth distribution or seismic gap below these three faults can be explained by their different geothermal gradients. These findings provide new insights for understanding the stress environment of the seismogenic zone in the SETP. Our model reveals the relationships between differential stress, seismicity, brittle–ductile transition, and boundary depth of the upper and lower crust in the continental crust, and connects the multiple observations from geophysics and geology. Furthermore, our model provides insights for studying multiple processes in the continental crust, such as crustal deformation, fault slip, and earthquake occurring.

Integrated seismotectonic model of the 2013 Mw6.6 and 2022 Mw5.8 earthquake sequences in the southwestern section of the longmenshan fault zone, China, and its implication

Integrated seismotectonic model of the 2013 Mw6.6 and 2022 Mw5.8 earthquake sequences in the southwestern section of the longmenshan fault zone, China, and its implication

Regional stress field in the SE margin of the Tibetan Plateau revealed by the focal mechanisms of small and moderate earthquakes

In this study, we investigate the stress field in the southeastern margin of the Tibetan Plateau. We first determine the focal mechanism solutions of 1537 small and moderate (3.2 ≤ MW ≤ 6.7) regional earthquakes from January 2009 to June 2021, and then use the focal mechanisms to invert for the spatial variation of crustal stress field by a damped linear inversion method. Our result suggests that in the southeastern margin of the Tibetan Plateau the seismogenic zone is in the upper crust above 15-km depth, and the stress field is predominantly strike-slip in the Sichuan-Yunnan Rhombic Block (SYRB). The maximum compressional stress axis is oriented in a fan-shaped pattern, rotating clockwise from nearly east-west in the Songpan-Ganzi Terrain in the north to northwest-southeast in the SYRB to nearly north-south across the Red River Fault in the Indo-China Block (ICB), consistent with the GPS-derived surface strain rate. The stress field around the border of the Tibetan Plateau with high elevation relief appears to be largely caused by gravitational effect with the maximum extensional axis perpendicular to the topography gradient. The stress field in the vicinity of the Longmenshan Fault Zone and in the Yangtze Craton is mainly thrust as a result of the eastward expansion of the Tibetan Plateau and the resistance of the Sichuan Basin. Near the epicenter of the 2008 Wenchuan earthquake and in the northeastern end of the Longmenshan Fault Zone, the thrust stress field shows spatial variations as a result of the perturbation by complex geometry and the post-seismic healing process. Our result provides multi-resolution images of the stress field for better understanding about the mechanisms of seismic activity and crustal deformation in the southeastern margin of Tibetan Plateau.

Three-dimensional fault model and features of chained hazards of the Luding MS 6.8 earthquake, Sichuan Province, China

The MS 6.8 Luding earthquake in 2022 is located on the NNW-trending Moxi segment of the Xianshuihe fault with left-lateral strike-slip behavior. This area is where the Xianshuihe, Anninghe, Daliangshan and Longmenshan faults intersect. China Earthquake Administration has identified that intersection area, among the Moxi segment of the Xianshuihe fault, the Anninghe fault, the Daliangshan fault and the southern part of the Longmenshan fault, as a high-magnitude earthquake hazard area. According to existing data on the Luding earthquake, including the focal parameters, the spatial distribution of re-located aftershocks, dominated azimuth of the earthquake intensities and earthquake-induced ground fissures, we built a 3D earthquake fault model. We found that two discontinuous NNW-trending vertical strike-slip faults with left stepping were the seismogenic faults of the Luding earthquake. Its coseismic left-lateral dislocation triggered transtensional slips and aftershocks on the NW-trending secondary faults at its northernmost tensile area. Meanwhile, local crustal coseismic shortening on the side of Mt. Gongga triggered the aftershocks on the NE- and NW-trending secondary conjugated strike-slip faults, which were confirmed by GNSS observations and InSAR deformation field around the epicenter. This earthquake rupturing pattern also controlled the spatial distribution of the earthquake intensity IX area and earthquake chain hazards. The Coulomb stress calculation shows that the Luding earthquake increases the risk of high-magnitude earthquake occurrence on the southernmost part of the Xianshuihe fault and the Anninghe fault. Finally, we suggested doing good monitoring of the Anninghe fault and the southernmost part of the Xianshuihe fault and avoiding active faults with seismogenic capacity and areas prone to earthquake-chained hazards during the site selection and planning of reconstruction.

Anisotropic tomography of eastern Tibet and its uncertainty from hypocentral errors

SUMMARY The mechanism responsible for the lateral expansion and uplift of the eastern Tibetan Plateau remains a topic of ongoing debate, partly due to discrepancies in the results of seismic velocity and anisotropy. In local earthquake tomography, hypocentral uncertainties can cause significant errors in the tomographic model. However, this issue has received limited attention in previous studies. In this work, we employ the weighted least-squares (WLS) method to solve the tomographic inversion problem. A power exponent coefficient, which is called weighting level, is introduced into the weighting matrix to control the relative contribution of the data with different hypocentral errors to the final tomographic result. Our data set contains high-quality Pg, Pn and Sg arrival times of local earthquakes recorded by the dense Chinese seismic network in eastern Tibet during 2008–2022. We comprehensively analyse the inversion results derived from the WLS inversions with different weighting levels to evaluate the robustness of isotropic velocity anomalies and azimuthal anisotropy. The most robust feature of our results is a striking low-velocity (low-Vp) zone surrounded by high-velocity (high-Vp) anomalies and fault parallel fast-velocity directions (FVDs) of azimuthal anisotropy in the lower crust beneath the western side of the Longmenshan fault zone. Taking into account many previous results of the region, we deem that the low-Vp zone reflects hot and wet upwelling flow from the deep asthenosphere, which ascends to the lower crust along the fault zone. At the NE margin of the Tibetan Plateau, significant low-Vp anomalies exist in the lower crust and the FVDs are consistent with the motion direction of the Tibetan block revealed by GPS (Global Positioning System) observations. We think that lower crustal flow exists beneath NE Tibet, which controls the plateau expansion toward the northeast. A low-Vp anomaly appears at 30 km depth beneath the Sichuan Basin. However, as the weighting level increases, the amplitude of this low-Vp anomaly decreases by more than 6 per cent, suggesting that this low-Vp anomaly has a lower accuracy than the other features.

Spatio-temporal characteristics of seismic strain anomalies reveal seismic risk zones along the Longmenshan fault zone and adjacent areas

Introduction: Identifying and quantifying earthquake precursors, and analyzing their physical mechanisms, continues to be a challenge for earthquake forecasting. In this study, orthogonal functions were developed to effectively identify precursor anomalies, thereby improving the forecasting of strong earthquakes.Methods: To study the spatio-temporal contour anomalies in seismic strain fields, we assessed them for seismic activity variables and natural orthogonal function expansion, in six strong earthquakes near the Longmenshan Fault Zone, China, that have occurred since 2008.Results: We observed that, prior to these earthquakes, the temporal factor (the time variation characteristics of the strain field) displayed anomalies with high/low values exceeding the mean square error within a stable context. The anomalies exhibited multi-component characteristics and were primarily concentrated in the first four-strain fields. Short-term and impending-earthquake anomalies were observed in the temporal factor before the 2008 Wenchuan (M8.0) and 2013 Lushan (M7.0) earthquakes, while medium-term and long-term anomalies appeared before the other four strong earthquakes, without notable short-term anomalies. The temporal evolution of strain field contour anomalies, and the strain contours positive and negative intersection, showed that central areas surrounded by multiple strain field contour anomalies were potential locations for strong earthquakes. This suggests a potential approach for earthquake location forecasting. Since 2009, there have been five strong earthquakes, each affected to varying degrees by anomalous strain fields from the 2008 Wenchuan (M8.0) earthquake.Conclusion: The results of this study corroborate the findings of the focal mechanism’s node shear stress, indicating significant physical implications of the anomalies and the reliability of these conclusion.

Cenozoic tectonic transition within the western segment of the Longmenshan fault, southeast margin of the Tibetan Plateau: Insights from geological and geophysical data

Thrust fault zones around the Tibetan Plateau (TP) record the tectonic evolution between the Plateau and its surrounding terranes, which is helpful for understanding the uplift mechanism and deformational processes of the TP. The Longmenshan fault (LMSF) is the tectonic boundary (TB) between the Yangtze terrane (YT) and Songpan-Garze terrane (SGT), while the TB of its western segment, either the Lijiang-Xiaojinhe fault (LXF) or Jinhe-Qinghe fault (JQF), is controversial. Therefore, we conducted magnetotelluric (MT) imaging, surface structure surveys, and petrologic analysis to further determine the deep–shallow structural relationship of the western LMSF segment. Resistors R1 and R2 revealed by MT imaging may have originated from different magmatism. Among them, R2 may have originated from plume underplating, which is consistent with previous studies, while R1 may have originated from the residue of episodic mafic magma intrusion along the JQF over a broader period. Based on regional geophysics, surface structural patterns and petrologic mineralogy, it is suggested that the JQF may have deformed deep into the lower crust or upper mantle, accommodating the southeast expansion of the TP by thrusting and acting as the TB between the YT and SGT before ∼15 Ma. After ∼15 Ma, due to the activation of the large-scale strike-slip faults, the LXF gradually replaced the JQF to dominate the structural deformation of the western LMSF segment. Our results indicate that the above tectonic transition might be associated with the geodynamic process from centralized deformation to diffuse deformation within the southeast TP during the late Cenozoic.

Topography effect on seismic waveform tomography: a quantitative study

SUMMARY In seismic tomography practices the Earth's surface is sometimes assumed to be either spherical or flat for convenience in forward modelling calculations. The effect of irregular surface topography on seismic wave propagation is thus ignored, resulting in biases in the phases and amplitudes of synthetic seismograms, which contribute to the residuals that are mapped into velocity structures in tomography inversions. In this study, we conduct a series of inversion experiments based on the adjoint waveform tomography method to quantitatively assess the topography effect on waveform-based inversion results. We first employ models with simplified topography to better highlight and quantify the topography effect. Results show that when topography effect is ignored in the forward modelling, it is mapped into velocity perturbations, leading to spurious velocity anomalies in tomography models. The strength of the spurious velocity anomaly is quasi-linearly related to locally averaged topography gradient. Our inversion experiments demonstrate that in places of strong topography variation, such as the Longmenshan Fault Zone region where the 2008 Mw 7.9 Wenchuan earthquake occurred, topography effect can lead to spurious relative velocity anomalies of up to 10 per cent, which cannot be ignored in waveform-based tomography inversions.

Asymmetrical Microfracture Density Across an Active Thrust Fault: Evidence from the Longmen Shan Fault, Eastern Tibet

Abstract Microfracture density in fault damage zones can reflect spatial variability that decays in intensity as a function of distance from the fault, which is crucial in understanding the mechanical, seismological, and fluid-flow properties of the fault system. However, few studies explored the characteristics of fracture density between the two sides of active dip-slip faults due to rare field observations. Here, we measured and modeled microfractures across an active thrust fault associated with the 2008 Mw 7.9 Wenchuan earthquake in the Longmen Shan, eastern Tibetan Plateau. The results showed that the microfracture density at the Qingping site developed more intensely in the hanging wall than in the footwall for an exposed thrust fault, indicating an asymmetrical pattern. The hidden thrust fault at the Jushui site showed that microfractures developed more intensely in vertical planes in the hanging wall than in the footwall, whereas microfractures developed similarly in horizontal planes within the two sides, indicating a quasiasymmetrical pattern. Comparing the data at the two sites with computational modeling, we suggest that fault geometry might exert a first-order control of the asymmetrical microfracture density pattern, which is helpful for revealing different deformational behaviors of rock masses in the fault damage zones and better understanding the hanging-wall effect for evaluating seismic hazards on active thrust faults.

Electrical Structure of Southwestern Longmenshan Fault Zone: Insights into Seismogenic Structure of 2013 and 2022 Lushan Earthquakes

On 1 June 2022, a magnitude 6.1 earthquake struck the southern segment of the Longmenshan fault zone on the eastern edge of the Tibetan Plateau, once again causing casualties and economic losses. Understanding the deep-seated dynamic mechanisms that lead to seismic events in the Lushan earthquake area and assessing the potential hazards in seismic gap areas are of significant importance. In this study, we utilized 118 magnetotelluric datasets collected from the Lushan earthquake area and employed three-dimensional electromagnetic inversion with topographic considerations to characterize the deep-seated three-dimensional resistivity structure of the Lushan earthquake area. The results reveal that the Shuangshi–Dachuan fault in the Lushan earthquake area can be divided into two relatively low-resistivity zones: a western zone dipping southeastward and an eastern zone with a steeper slightly northwestern dip. These two zones intersect at a depth of approximately 20 km, forming an extensional pattern resembling a “Y” shape. The epicenters of both the 2013 and 2022 Lushan earthquakes are primarily located in the upper constricted portion of the pocket-like low-resistivity body at depth. The distribution of seismic aftershocks is confined within the region enclosed by the high-resistivity body, following the pattern of the Y-shaped low-resistivity zone.

Possible layered lithospheric anisotropy around Longmenshan Faults by teleseismic S wave splitting and receiver functions

This study conducts an in-depth analysis of seismic anisotropy around the Longmenshan Faults. Utilizing a dataset of about 7710 earthquake catalogs from July 2007 to March 2023, we applied S wave splitting and receiver function methods to examine Pms and XKS waveforms collected from 12 fixed broadband stations across Gansu and Sichuan provinces. Our analysis revealed significant variations in seismic anisotropy between the crust and lithosphere, marked by distinct fast wave directions and delay times. These characteristics point to the possibility for layered anisotropy within the region. A two-layer anisotropy inversion analysis at key stations further delineated the differential anisotropic behaviors between the crust and the upper mantle, underscoring the impact of local geological structures and mantle dynamics. Crucially, our study posits the existence of layered anisotropy around the Longmenshan Fault Zone, a finding that significantly advances our comprehension of the region’s seismic anisotropy and adds a vital dimension to our understanding of its subsurface structural intricacies and tectonic evolution.

Earthquake monitoring and high-resolution velocity tomography for the central Longmenshan fault zone by a temporary dense seismic array

Earthquake monitoring and high-resolution velocity tomography for the central Longmenshan fault zone by a temporary dense seismic array

Crustal heterogeneity effects on coseismic deformation: numerical simulation of the 2008 MW 7.9 Wenchuan earthquake

Coseismic deformation of large earthquakes causes significant property damages and fatalities, which requires quantitative research of multiple disciplines such as geodesy, geological investigation, seismic tomography, and seismic dislocation theory. The finite element method accounts for material heterogeneity and geometric complexity, making it suitable for studying the coseismic deformation of large earthquakes. This paper develops a parallel elastic finite element program that utilizes split nodes and high-performance parallel computing technology on the FELAC software platform to study the coseismic deformation of large earthquakes. We verify the accuracy of the parallel elastic finite element program by comparing its results with the analytical solutions from seismic dislocation theory for four ideal earthquake cases. Finally, we established parallel elastic finite element models to study the coseismic deformation of the 2008 Wenchuan earthquake. The simulation results are consistent with the GPS and InSAR data. Coseismic surface deformation results are significantly influenced by medium regional heterogeneity with different layered structures besides the Longmenshan fault. The finite element program lay the foundation for the inversion of the coseismic fault rupture process based on the heterogeneous medium model and complex geometric model.

The crustal structure of the Longmenshan fault zone and its implications for seismogenesis: new insight from aeromagnetic and gravity data

Abstract. Although many geophysical models have been proposed in the Longmenshan fault zone (LFZ) and its surrounding areas, the deep structure of the seismic gap and its constraint of the Wenchuan and Lushan earthquakes remain uncertain. Based on the compiled aeromagnetic data and Bouguer gravity data, we have tried to create a more detailed and reasonable magnetic and density model using 2D forward modeling and 3D inversion and made the deep structure of the LFZ visible. The research shows that structure is heterogenous across the LFZ. The earthquake epicenters are located in regions with high-magnetic anomalies and gravity gradients that are associated with rigid blocks that were likely to accumulate stress. However, the seismic gap shows low-magnetic anomalies and transition of gravity anomalies related to a weak zone. The Sichuan Basin has two NE-trending banded high-magnetic blocks extending beneath the LFZ that firmly support the idea that the crust of the Sichuan Basin has subducted downward the LFZ. More importantly, the basement subducts to approximately 33 km west of the Wenchuan–Maoxian fault, with a low dip angle beneath the middle segment of the LFZ, whereas the distance decreases to approximately 17 and 19 km under the southern segment. Thus, the crust of the Sichuan Basin beneath the middle segment extends farther than that beneath the southern segment, with the seismic gap as the transition zone. Therefore, we propose that the structural heterogeneity of the basement on the western margin of the Sichuan Basin may be the main reason for the different focal mechanisms and geodynamics of the Wenchuan and Lushan earthquakes.

Diverse Late Triassic arc magmatism and porphyry Cu deposits in the southeastern Tibetan Plateau: Related to slab tear along a subducting passive margin?

The subduction of the ocean–continent transition (OCT) zone on a passive margin has received limited attention, despite its significant impact on tectono-magmatic processes. This study focuses on the Late Triassic basalts and andesites in the southeastern Tibetan Plateau, where the Yushu-Yidun Arc belt intersects the Longmen Shan Fault system. The goal is to understand the interplay between these linear systems. The basalts exhibit trace element and isotopic signatures resembling those of within-plate enriched mid-ocean ridge basalts or ocean island basalts (EOBs), such as positive NbTa anomalies, high Nb/Th (7.72 to 10.3) and εNd(t) (+2.23 to +6.34) values, and low La/Sm values (2.69 to 4.05), suggesting an asthenospheric mantle source. The contemporaneous andesites display typical arc-like geochemical characteristics (e.g., enrichment in Ba, Th, and U but depletion in Nb, Ta, and Ti), with negative εNd(t) values (−2.41 to −1.08), indicating a mantle wedge source affected by subducted components. Geological evidence indicates that the Longmen Shan Fault zone, previously defined by its magmatic rocks, was an active OCT belt during the Late Triassic, separating the South China Continent from the Garze–Litang Ocean. Consequently, we propose a model suggesting the tearing of a subducted oceanic slab from the continental lithosphere within the OCT belt under the southernmost segment of the Yushu–Yidun Arc belt, driven by strong slab pull in the Garze–Litang Ocean subduction zone. This model effectively explains the intriguing geology of the Songpan–Garze–Hoh Xil Terrane, in which Late Triassic EOBs and normal arc-related igneous rocks coexist and are accompanied by abundant slab-derived mineralized porphyries. Our findings have global implications for investigating geochemically diverse magmatism and associated mineralization on passive continental margins.