Loading...

Table of Content

    20 June 2025, Volume 47 Issue 3
    THREE-DIMENSIONAL MODEL OF SEISMOGENIC FAULT AND SEISMIC ENVIRONMENT OF XIZANG DINGRI MS6.8 EARTHQUAKE OF JANUARY 7, 2025
    GUO Zhao-wu, LU Ren-qi, ZHANG Jin-yu, FANG Li-hua, LIU Guan-shen, WU Xi-yan, SUN Xiao, QI Shi-miao
    2025, 47(3):  671-688.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250030
    Asbtract ( )   HTML ( )   PDF (7098KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    At 09:05a.m. on January 7, 2025, a magnitude MS6.8 earthquake struck Dingri County, Xizang, China, resulting in 126 fatalities and a maximum seismic intensity of Ⅸ. Occurring within a seismically active and tectonically complex region, this event drew significant attention from both the scientific community and the public. The epicenter was located near the Dengmecuo Fault, which has been preliminarily identified as the seismogenic fault.

    This study utilized publicly available geological survey data, aftershock relocations, and focal mechanism solutions to construct a detailed three-dimensional geometric model of the Dengmecuo Fault. The model was developed using the SKUA-GOCAD 3D modeling platform, enabling a comprehensive analysis of the fault’s geometry. Results reveal pronounced geometric segmentation along the fault plane, with the spatial distribution of these structural features closely correlating with observed seismicity, highlighting the influence of fault geometry on earthquake generation.

    The MS6.8 Dingri earthquake occurred near a prominent structural irregularity on the Dengmecuo Fault, at point P3, where the fault plane bends into an eastward-projecting arc. This three-dimensional structural mutation likely played a role in the nucleation of the event, underscoring the relationship between fault complexity and seismic rupture. The Dengmecuo Fault, situated in the southern Tibetan plateau, is a listric normal fault characterized by a steep upper section and a gentler lower section that terminates within a detachment layer in the upper crust. It does not extend into the deeper lithosphere, indicating that it is part of the region’s shallow normal fault system. The earthquake is interpreted as the release of accumulated stress along this shallow fault structure.

    To evaluate post-earthquake stress transfer and seismic hazard, Coulomb stress modeling(Coulomb 3.4)was conducted. The analysis indicates that several regional faults are now in a state of increased Coulomb stress, including the southern segment of the Dengmecuo Fault, the middle segment of the south Xizang detachment system, the southern segment of the Shenzha-Dingjie Fault, the central Yarlung-Zangbo Fault, and the midsection of the Dajiling-Angren-Renbu Fault. These fault segments are identified as potential sites for future seismic activity and merit heightened monitoring.

    This study presents a detailed characterization of the three-dimensional geometry of the seismogenic fault responsible for the Dingri MS6.8 earthquake and offers a preliminary analysis of regional seismogenic structures. The findings provide valuable insights into the tectonic setting of southern Xizang and contribute to the assessment of regional seismic hazard.

    LATE QUATERNARY THROW RATE OF THE SEISMOGENIC FAULT(DENGMECUO FAULT)OF THE 2025 MS6.8 DINGRI EARTHQUAKE IN SHIGATSE
    GAO Yang, WU Zhong-hai, HAN Shuai, TIAN Ting-ting
    2025, 47(3):  689-706.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250034
    Asbtract ( )   HTML ( )   PDF (10876KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    On January 7, 2025, a significant earthquake occurred in Dingri, Shigatse, China. Both the China Earthquake Networks Center and the USGS provide focal mechanism solutions indicating that the earthquake was a normal fault event. The epicenter was located in the Dengme Cograben, part of the southern segment of the Dinggye-Xainza rift, with the seismogenic fault identified as the bounding normal fault of the Dengme Cograben, the Dengmecuo fault. The late Quaternary throw rate of this fault is crucial not only for regional seismic hazard assessments but also for understanding the east-west extensional deformation within the Tibetan plateau. Previous studies on the kinematics of the Dengmecuo fault report varying throw rates: (0.28±0.04)mm/a since ~56ka, (0.28±0.04)mm/a since ~50ka, and(0.09±0.03)mm/a since ~95ka. However, the observed maximum coseismic vertical displacement of ~3m during the January 2025 earthquake suggests that the maximum cumulative time for this displacement is approximately 37ka, based on the published late Quaternary throw rate. This is inconsistent with the ~5ka elapsed time since the most recent earthquake on this fault, highlighting the uncertainty in the late Quaternary throw rate and limiting our understanding of strain partitioning and seismic risk in the southern Dinggye-Xainza rift.

    To resolve this uncertainty, we combine high-resolution remote sensing image interpretation with field geological and geomorphological surveys to determine the geometric distribution of the Dengmecuo fault. Two study sites were selected along the fault, each featuring clear faulted geomorphology suitable for dating. Low-altitude photogrammetry using an unmanned aerial vehicle(UAV)was combined with optically stimulated luminescence(OSL)and AMS 14C dating to refine the late Quaternary throw rate of the Dengmecuo fault and assess its seismic hazard.

    Our results show that the Dengmecuo fault, which strikes NNW and dips NW or W, is approximately 58km long and can be divided into northern, central, and southern segments. The northern segment is defined by the Laangshuiku area, while the southern boundary is located at the edge of Pum Qu. In the northern study site, vertical offsets of 19.2(+3.5/-2.3)m on the T2 alluvial fan and 8.0(+0.9/-0.7)m on the T1 alluvial fan correspond to formation ages of(28.3±1.4)ka and(12.0±1.5)ka, respectively. By matching the vertical offsets with their respective formation ages, we estimate a throw rate of (0.7±0.1)mm/a. However, the throw rate for the T1 fan is uncertain, as its vertical offset is smaller than the cumulative displacement since its formation. At the southern study site, combining a vertical offset of 5.3(+0.3/-0.5)m with the formation age of the T1 terrace ((9.2±1.0)ka), we calculate a throw rate of(0.6±0.1)mm/a.

    Overall, our results indicate late Quaternary throw rates of(0.7±0.1)mm/a since ~28ka and(0.6±0.1)mm/a since the Holocene. Additionally, using the relationship S=D/Rx, where S is the slip rate, D is displacement, and Rx is the recurrence interval, we estimate a slip rate of(0.5±0.1)mm/a based on the average value of the maximum displacement(2.5~3m)of 2025 Dingri MS6.8(Shigatse)earthquake as 2.75m and the interval of paleoseismic events during the Holocene as(5500±1100)a. This result is consistent with the throw rate of(0.6±0.1)mm/a since the Holocene determined from faulted geomorphic surfaces.

    Finally, combining the throw rate of (0.6±0.1)mm/a since the Holocene with the ~5ka elapsed time since the most recent earthquake, we conclude that the Dengmecuo fault had accumulated 2.5~3.5m of coseismic displacement before the 2025 Dingri MS6.8 earthquake, corresponding to a magnitude of MW6.9-7.0. These parameters align with the 2025 MS6.8 earthquake, indicating that the Dengmecuo fault posed a significant seismic risk prior to the event.

    THE MAXIMUM VERTICAL DISPLACEMENT OF THE MS6.8 EARTHQUAKE IN XIZANG AND ITS SURFACE DEFORMATION STYLE
    ZHANG Da, SHI Feng, LUO Quan-xing, QIAO Jun-xiang, WANG xin, YI Wen-xing, LI Tao, LI an
    2025, 47(3):  707-717.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250041
    Asbtract ( )   HTML ( )   PDF (9462KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    On January 7, 2025, at 9:05 AM, a MS6.8 earthquake occurred in Dingri County, Shigatse City, Xizang Autonomous Region. The epicenter was located at 87.45°E and 28.50°N, with a focal depth of 10 kilometers, as determined by the China Earthquake Networks Center. Regarding the maximum vertical displacement of this earthquake, there are differing opinions due to variations in the reference markers used during field measurements of coseismic displacement or deformation amplitude, as well as divergent understandings of the deformation style of earthquake scarps and pre-existing scarps. Additionally, factors such as the complex structure of the surface rupture zone and the short duration of field investigations contribute to the current discrepancies in understanding the maximum co-seismic displacement. Different scholars have varying perspectives on this matter. To determine the maximum co-seismic displacement and provide data for subsequent research, this study utilized GF1 satellite imagery and drone aerial surveying technology to measure the maximum co-seismic displacement following the earthquake.

    The seismogenic fault of the Dingri earthquake is the Dengmecuo Fault. The surface rupture caused by this earthquake is mainly concentrated near Nixiacuo, at the northern end of the Dengmecuo Fault. The surface rupture trace in the Nixiacuo section is evident, striking northeast and arranged in a linear pattern. It develops along the existing steep scar in front of the mountain, cutting through a series of landforms, including alluvial fans and moraines, and aligns well with the original steep scar. This section of surface rupture is large in scale, developing a series of extensional fractures and fault scarps. The surface rupture phenomenon is most pronounced at a location 800 meters north of Nixiacuo, where the largest co-seismic vertical displacement also occurs. At the location of the maximum co-seismic vertical displacement, the main rupture diverges into several secondary ruptures, which later converge back onto the main rupture. Simultaneously, several extensional fractures develop, and a clear vertical displacement is observed between the hanging wall and the footwall, accompanied by rock collapse. This study employed field observations and UAV aerial survey technology to determine the maximum coseismic surface displacement resulting from the Tingri earthquake. Since the surface rupture generated not only vertical displacement but also ground fissures, leading to some horizontal displacement, the displacement could not be measured directly. To address this, we identified two points on the hanging wall and footwall of the fault, measured the distance l between them, and determined the inclination angle θ of l relative to the ground. The vertical displacement was then calculated using the trigonometric relationship l ×sinθ. Four sets of measurements were taken, yielding a result of(2.47±0.1)m. UAV aerial survey technology was used to capture orthophotos of the location with the maximum coseismic vertical displacement. Three profiles were measured, and the largest recorded coseismic vertical displacement was 0.2m.

    In this study, we collected empirical formulas derived from the maximum moment magnitude of the Dingri earthquake inverted by other scholars. The calculated maximum co-seismic surface displacement ranged from 2.37m to 2.97m, which is consistent with the(2.47±0.1)m observed in the field and the(2.6±0.2)m obtained using drone aerial survey technology. The moment magnitude of the Yutian earthquake is equal to that of the Dingri earthquake, but its maximum co-seismic surface displacement is smaller than that of the Dingri earthquake. Firstly, the Yutian earthquake exhibited both left-lateral strike-slip and normal fault characteristics, whereas the Dingri earthquake was mainly characterized by normal faulting. Secondly, the Yutian earthquake produced relatively continuous surface ruptures, whereas the Dingri earthquake produced highly discontinuous surface ruptures with significant displacement differences between segments.

    LATE QUATERNARY ACTIVITY OF THE XIETONGMEN TO DENGMECUO SEGMENT ALONG THE XAINZA-DINGGYE RIFT IN SOUTHERN QINGHAI-XIZANG PLATEAU
    WANG Duo, CHEN Li-chun, LI Yan-bao, WANG Hu, JIA Yong-shun, GAO Yin-yi, XUE Ke-yi
    2025, 47(3):  718-733.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250012
    Asbtract ( )   HTML ( )   PDF (16018KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    The nearly north-trending rift system and nearly east-west strike-slip faults are the major structures accommodating the east-west extensional deformation within the Qinghai-Xizang Plateau. The rift systems are roughly separated by the nearly east-west striking Indus-Yarlung suture zone into north and south parts. The eastern rift systems are a strongly active major seismic zone in Xizang, which is characterized by faulted landforms of alluvial fans, river terraces, and moraines with large magnitude earthquakes. The Yadong-Gulu rift experienced the 1411 Damxung M8 earthquake and the 1952 Gulu M7.4 earthquake. On January 7, 2025, a MS6.8 earthquake occurred along the southern segment of the Xainza-Dinggye rift in Dingri County, Xigaze, which caused widespread concern about the seismic and Late Quaternary active behaviors along the rift systems. However, few studies on fault activity at the junction of the north and south segments of the rift systems were conducted along both side of the Indus-Yarlung suture zone, which greatly hinders us from understanding the active deformation process and seismic activity of the rift systems in southern Qinghai-Xizang Plateau.

    The Xainza-Dinggye rift can be divided into the north and south segments, named the Xainza and Dinggye rifts, by the Indus-Yarlung suture zone. The Xietongmen and Dengmecuo segments are the northern and southern adjacent segments of Xainza-Dinggye rift, respectively. The Xietongmen segment as the south end of the Xainza rift intersect with the Indus-Yarlung strike-slip fault. There are no strong earthquakes of over M6 recorded near the Xietongmen segment, and the small and medium earthquakes are significantly less than those of other segments. The 2025 Dingri earthquake has caused great social concern about whether a larger earthquake will occur between the Dengmecuo and Xietongmen segment of the Xainza-Dinggye rift, especially in the densely populated area of Xietongmen County. A recent earthquake risk survey has not yet found evidence of Holocene activity along the Xietongmen segment. Meanwhile, previous studies suggest very low activity during the late Quaternary along the Dengmecuo segment, which is in stark contrast to current strong earthquake activities. Therefore, there are still great uncertainties about the fault geometry and activity of the Xainza-Dinggye rift on both sides of the Indus-Yarlung suture zone.

    To determine the Late Quaternary activity of the Xietongmen to Dengmecuo segments of the Xainza-Dinggye rift, we used remote sensing interpretation, field survey, optically stimulated luminescence, and radiocarbon dating methods on the displaced landforms. We found new evidence of the latest activity on the eastern and western branches of the Xietongmen segment at the past millennium. The western branch of the Xietongmen segment has crossed the Yarlung Zangbo River southward, then terminates at the intersection basin with the latest active branches of southern Yarlung Zangbo fault. Our results also suggest the Dengmecuo segment had strong activity with several hundred meters width since the late Quaternary. The latest faulting has extended northward into the mountainous area according to the surface ruptures and positioning aftershocks of the MS6.8 earthquake. The nearly east-west striking Indus-Yarlung suture zone constitutes a north-trending structural gap zone between the two segments. The fault geometries, displaced landforms, and seismic activity of these two segments reveal that they are probably approaching by cutting through this structural gap, and the seismic risk cannot be ignored in the future. Due to the limitations of existing data, further detailed field investigations, geodetic observations, and geophysical deep data are needed to verify and improve our speculation. Our results of the Late Quaternary activity of the Xietongmen to Dengmecuo segments of the Xainza-Dinggye rift provide scientific support for seismic risk assessment of national engineering projects and post-disaster reconstruction of the recent MS6.8 earthquake in the Xigaze region.

    JOINT INVERSION OF THE 2025 DINGRI MS6.8 EARTHQUAKE RUPTURE PROCESS BASED ON TELESEISMIC P WAVES, STRONG-MOTION AND INSAR DATA
    XU Yue-yi, XU Bei-bei, XU Chen-yu, SHAO Zhi-gang, HU Chao-zhong
    2025, 47(3):  734-746.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250033
    Asbtract ( )   HTML ( )   PDF (8453KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    At 09:05 Beijing time on January 7, 2025, a MS6.8 earthquake struck Dingri County, Shigatse City, Xizang Autonomous Region, China, with a focal depth of 10km. The earthquake caused strong ground shaking, reaching a maximum intensity of Ⅸ on the China Seismic Intensity Scale, and resulted in 126 fatalities. This seismic event occurred along the Dengmoco Fault, which is part of the north-south-trending Shenzha-Dingjie normal fault system in the southern Qinghai-Xizang Plateau. The fault extends approximately 60km, dips westward, and has maintained a vertical slip rate of(0.28±0.04)mm/a over the past 1 000 years. Geological investigations indicate that the most recent paleo-earthquake on this fault occurred 4 800-4 968 a BP, with an estimated recurrence interval of(5 500±1 100)a in the Holocene. These characteristics demonstrate that the 2025 Dingri earthquake represents a typical event within the fault’s long-term seismic cycle.

    To determine the source mechanism of this event, we performed a joint inversion of the moment tensor using teleseismic P waves and W-phase. By combining the relatively high-frequency P-waveforms with the low-frequency W-phase records, this approach enables a more robust determination of the focal mechanism while providing improved constraints on the centroid location, particularly the centroid depth. The results indicate a moment magnitude of 7.02, with a centroid time offset of 9.7s and a shallow centroid depth of 6km. The optimal centroid is located at(28.6°N, 87.5°E), about 24° east of north from the epicenter. The best-fitting focal mechanism solution yields two nodal planes: Nodal Plane 1 with a strike of 344°, dip of 48°, and rake of -105°; and Nodal Plane 2 with a strike of 185°, dip of 44°, and rake of -74°. Considering the west-dipping geometry of the Dengmoco Fault, Nodal Plane 2 is interpreted as the likely fault plane responsible for the rupture. Based on Nodal Plane 2, we further conducted a joint rupture process inversion using teleseismic P waves, strong-motion waveforms, and InSAR deformation data. The combination of seismic and geodetic observations provided complementary constraints, enhancing both the spatial and temporal resolution of the dynamic rupture process. The results indicate a predominantly normal faulting mechanism with a minor left-lateral component. The rupture propagated mainly northward, with limited southward extension. The rupture lasted approximately 36 seconds, with the main slip occurring between 8 seconds and 24 seconds and concentrated within depths of 0-10km, generating a significant surface rupture approximately 20km north of Changsuo Township. To further explore the fault geometry, we conducted a grid search using a dual-fault model to identify the optimal strike for the northern segment. The analysis identified 240° as the best-fit strike, which slightly improved the overall data fitting and exhibited better consistency with the surface topography. The resulting slip model retained the main rupture characteristics observed in the single-fault scenario. The rupture process can be divided into three stages: ① 0-7s initial nucleation near Cuoguo Township with a relatively minor slip(<1.3m); ② 8-24s primary rupture occurred approximately 20km north of Changsuo, with a peak slip of 4.3m; ③ 25-36s rapid slip attenuation after propagating beyond the northern fault bend.

    Our study indicates that the complex fault structure plays a critical role in rupture dynamics. The relatively minor slip observed between Cuoguo Township and Changsuo Township, together with the sparse aftershock activity, suggests the existence of an unidentified east-dipping, NW-trending fault segment in this region. Additionally, the abrupt rupture termination near the northern segment of the Dengmoco Fault is likely influenced by an unrecognized NE-trending subsidiary fault. This study underscores the importance of rupture directivity in seismic hazard assessment and reveals the structural complexity of the Dengmoco Fault system. Overall, the results enhance our understanding of the seismogenic mechanism and contribute to more accurate earthquake hazard evaluations in the region.

    PRELIMINARY STUDY ON THE SEISMOGENIC STRUCTURE OF THE 2025 XIZANG DINGRI MS6.8 EARTHQUAKE SEQUENCE
    CHEN Han-lin, WANG Qin-cai, GAO Jin-rui, LI Jun
    2025, 47(3):  747-760.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250021
    Asbtract ( )   HTML ( )   PDF (7688KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    This study investigates the MS6.8 earthquake that occurred in Dingri, Xizang(Tibet), on January 7, 2025, using the azimuth spectrum fast detection method. The moment tensor solution, relative position between the seismogenic point and moment centroid, and other key parameters were determined. The seismogenic nodal plane was further constrained using this approach. The method involves dividing the multi-dimensional parameter space into subspaces and applying a combination of grid search and gradient descent to identify the optimal solution with the least parameter misfit in each subspace. Waveform data were collected from 28 fixed and mobile seismic stations, with relevant station and instrument response data provided by the Institute of Earthquake Forecasting of the China Earthquake Administration. Twelve stations with epicentral distances ranging from 100km to 450km and a filtering range of 0.02-0.05Hz were selected for analysis. A time-domain full-band inversion was employed to incorporate more waveform information and enhance constraints on source parameters such as spatial location and rupture direction. Observed and theoretical three-component waveforms were compared, and records with low signal-to-noise ratios or poor fit were excluded. Final calculations were performed using data from seven stations. The resulting moment tensor solution indicates two nodal planes with strike/dip/rake values of(188°, 46°, -90°)and(9°, 44°, -90°), a focal depth of 8.9km, and a moment magnitude(MW)of 7.0. To test solution stability, multiple initial input combinations(strike, dip, and rake)were examined, yielding consistent inversion results. Comparison with results from the USGS, GFZ, and ZHANG Zhe(China Earthquake Administration)shows discrepancies within 24°, 6°, and 13°, respectively, for strike, dip, and rake. Further analysis of the source geometry yielded a boundary radius of 10km, relative rupture velocity of 0.7, and relative distances of 4.5km along strike and 0km along dip between the seismogenic point and centroid. The rupture propagated both upward and downward from the hypocenter. The seismogenic fault plane was identified with strike/dip/rake parameters of 188°, 46°, -90°.

    To further investigate the earthquake sequence, we analyzed phase reports of 3, 545 events from January 7 to 14, 2025, provided by the China Earthquake Networks Center. The HypoDD method was used for relocation, with events recorded within 300km of the cluster center. Parameters included a maximum inter-event distance of 10km and a minimum of 8 links per pair. A total of 197, 898 P-wave and 252, 651 S-wave differential times were successfully used, representing 80% and 81% of the total available data, respectively. Relocation was performed using the conjugate gradient method and a one-dimensional velocity model by Monsalve et al. for southern Tibet. Quality control parameters ranged from 40 to 80, resulting in successful relocation of 3, 155 events. The relocated hypocenters reveal that the sequence can be divided into three segments—southern, central, and northern. The southern swarm extends NW-SE from Guojia Town to Cuoguo, intersected by the Cuoguo and Dengmecuo faults. The central segment trends NNE-SSE from the mainshock along the Cuoguo fault toward Qiugu Village, with relatively sparse seismicity and a seismic gap near Changsuo Town. The northern segment continues NW-SE toward Xingdang and is intersected by the Nongqu fault. While the central swarm aligns with the Cuoguo fault, the southern and northern segments deviate from mapped fault trends, suggesting the presence of NW-SE-trending subsidiary faults.

    Depth profile analysis indicates that all three swarms occurred on west-dipping fault planes. The southern and central segments show clear layering with focal depths of 6-14km and 20-30km, while the northern segment shows less stratification. The spatial pattern suggests a complex, segmented fault system with a possible Z-shaped branch fault in the Shenzha-Dingjie normal fault zone. The Dingri earthquake sequence is therefore attributed to rupture within a complex fault network.

    HORIZONTAL DEVIATORIC STRESS AND ELASTIC LITHOSPHERE THICKNESS CHARACTERISTICS OF THE EPICENTER AND ITS ADJACENT AREAS OF THE DINGRI MS6.8 EARTHQUAKE, XIZANG, CHINA
    MENG Heng-zhou, YANG Guang-liang, QIN Hai-tao, TAN Hong-bo, LIU Sheng, WANG Jia-pei, HUANG Min-fu, ZHANG Ming-hui
    2025, 47(3):  761-776.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250018
    Asbtract ( )   HTML ( )   PDF (4809KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    According to the China Seismic Network, a magnitude MS6.8 earthquake occurred on January 7, 2025, in Dingri, Tibet. The epicenter was located near the Shenza-Dingjie rift(87.45°E, 28.50°N)at a focal depth of 10km. This earthquake is attributed to compressional forces resulting from the ongoing convergence between the Indian and Eurasian plates, which induce east-west(EW)extensional stress within the Tibetan plateau. In southern Tibet, a series of north-south(NS)trending rift valleys have developed in response to these tectonic processes. In these regions, variations in lithospheric density generate deviatoric stress, which closely correlates with the spatial distribution of seismicity. Areas exhibiting high deviatoric stress tend to experience more frequent tectonic activity. Furthermore, the elastic thickness of the lithosphere(Te), a key indicator of lithospheric strength and stability, is generally lower in seismically active zones, particularly within transition zones between strong and weak lithosphere. Therefore, analyzing deviatoric stress and Te in the study area is essential for understanding the mechanisms behind the Dingri earthquake and related seismic phenomena.

    To investigate the seismo-tectonic background of this event, this study constructs an equilibrium equation for deviatoric stress based on gravity field data. The admittance and coherence method is applied to estimate deviatoric stress at various depths, elastic lithospheric thickness(Te), and load ratio(F), using the WGM2012 gravity field model, ETOPO1 topographic data, and CRUST1.0 crustal structure data. The study further analyzes the coupling between deviatoric stress and regional geological structures, examines the spatial distribution of Te and its tectonic implications, and evaluates the influence of load ratio(F)on deviatoric stress estimation. These analyses form the basis for a comprehensive discussion of the focal characteristics of the Dingri earthquake.

    Our results indicate that deviatoric stress in the study area exhibits a clear south-north gradient, with higher values(>15MPa)concentrated in the southern region, particularly south of the Yarlung Zangbo Fault. In the north, elevated stress values are primarily associated with major fault zones such as the Shenza-Dingjie and Yadong-Gulu faults. Deviatoric stress decreases with depth, showing a marked decline at 50km. The elastic lithosphere thickness is generally greater than 40km across the region, with higher values observed in the central and southern areas, consistent with the subsidence and underthrusting of the Indian plate along the southern margin of the plateau. In contrast, lower Te values in the northeastern part of the study area are likely linked to rifting and lithospheric extension. The load ratio(F)varies between 0 and 1, with surface loads(F<0.4)dominating most of the region. However, higher values are observed in the northern segment of the Yadong-Gulu fault zone, suggesting a significant contribution from lower crustal or upper mantle processes. High load ratios can introduce uncertainties in deviatoric stress estimates, particularly in regions of active deep-seated tectonism.

    The epicenter of the Dingri MS6.8 earthquake is situated within the Shenza-Dingjie rift zone. The stress regime in this area is dominated by strike-slip tectonics, with NNW-SSE compression and NEE-SWW extension. Under this stress configuration, NS-trending faults near the epicenter are susceptible to normal faulting. Deviatoric stress values at crustal depths of 0, 10, 20, 30, and 50km are 11.45MPa, 8.46MPa, 4.36MPa, 2.86MPa, and 1.19MPa, respectively. These results indicate that deviatoric stress is predominantly concentrated within the upper 20km of the crust and is oriented mainly in the NNE direction, consistent with the regional tectonic stress field. Additionally, the epicenter lies within a transitional zone of elastic lithospheric thickness, where stress resistance varies, providing favorable conditions for shallow, NS-oriented normal faulting.

    INVERSION OF THE RUPTURE PROCESS OF THE XIZANG DINGRI MW7.1 EARTHQUAKE IN 2025
    LIU Sheng, TAN Hong-bo, YANG Guang-liang, MENG Heng-zhou, QIN Hai-tao, WANG Jia-pei, HUANG Min-fu
    2025, 47(3):  777-788.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250037
    Asbtract ( )   HTML ( )   PDF (5998KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    According to the United States Geological Survey(USGS), a moment magnitude(MW)7.1 earthquake struck Dingri County, Xigaze City, Xizang(28.65°N, 87.36°E)at 01:05:16 UTC on January 7, 2025(09:05:16 Beijing Time). The earthquake occurred at a focal depth of 10km and resulted in significant casualties: As of the afternoon of the same day, 126 deaths were confirmed, and approximately 61, 500 individuals were affected.

    The Dingri earthquake occurred in the southern Qinghai-Xizang Plateau, a region characterized by intense tectonic activity due to the ongoing subduction of the Indian Plate beneath the Eurasian Plate. This area exhibits the typical seismic pattern of “frequent large earthquakes and persistent smaller events.” The epicenter is situated near the intersection of the Shenzha-Dingjie fault zone, south of the Yarlung Zangbo fault zone, and the South Xizang Detachment fault zone. The Dangra Yongco-Xuru Fault lies to the west, and the Shenzha-Dingjie Fault to the east, the latter exhibiting a north-south extensional structure that divides the South Xizang Detachment Fault into eastern and western segments. GPS observations indicate extension rates of 4~5mm/a for both the Dangra Yongco-Xuru and Yadong-Gulu fault zones, while the Shenzha-Dingjie fault exhibits a slower rate of 1~2mm/a. According to historical USGS records, over 700 earthquakes with magnitudes above M3 have occurred in this region since the 20th century, including 604 events in the M3-M5 range, 101 in the M5-M7 range, and two above M7. Most of these events are concentrated along the Himalayan Orogenic Belt and near the Shenzha-Dingjie fault zone. The occurrence of the Dingri earthquake underscores the region’s seismic complexity and highlights the importance of studying rupture dynamics for understanding earthquake mechanisms and assessing seismic hazards. There is a close correlation between the source rupture process and the earthquake expansion law. In-depth study of the source rupture process is helpful for a comprehensive understanding and analysis of the inducing factors of earthquake rupture, the complexity of the source environment, and its potential impact. Therefore, the inversion of the rupture process of this earthquake can provide a reference for earthquake disaster analysis, earthquake emergency rescue, and post-earthquake seismic trend analysis.

    This study utilizes the rupture process inversion of the Dingri earthquake based on the source mechanism parameters(strike/dip/rake=187°/49°/-78°)provided by USGS and far-field waveform data from 51 stations within 30°~90° epicentral distances, sourced from the IRIS database. The analysis employs the AK135f global 1-D velocity model and the Iterative Deconvolution and Stacking(IDS)method proposed by Zhang (2014). The IDS method integrates advantages of both network-based and back-projection approaches and enables automated rupture process inversion without preset rupture time constraints. It has been successfully applied to events such as the 2015 Nepal and 2017 Jiuzhaigou earthquakes. The inversion results indicate an asymmetric bilateral rupture pattern with shallow rupture propagation. The maximum slip reached approximately 2.3 meters, with the rupture occurring primarily within a 0~9km depth range. The total seismic moment was 5.5×1019 N·m, corresponding to an MW of 7.1. The rupture lasted 29 seconds, peaking in moment release at 16 seconds, with most rupture ceasing by 28 seconds.

    The above results of this study align well with those of other studies, showing a maximum variation in magnitude of 0.1(range: MW7.0-7.2)and a slip difference of less than 1 meter(range: 1.5~3.2m). Despite this agreement, however, debate remains regarding whether the rupture was unilateral or bilateral. Contributing factors include variations in input source parameters, rupture initiation locations, station distribution, and inversion uncertainties. However, the distribution of aftershocks on both sides of the rupture supports the bilateral rupture interpretation. Based on these findings, this earthquake is interpreted as a normal-faulting event with an asymmetric bilateral rupture along the Shenzha-Dingjie fault zone. The concentration of slip near the surface suggests that upper crustal structures are more fragile and play a key role in seismic energy release, potentially explaining the severity of the disaster. These results emphasize the need for closer monitoring of near-surface fault slip potential in this region.

    COSEISMIC SURFACE RUPTURE OF THE MS6.8 DINGRI EARTHQUAKE IN XIZANG, CHINA, BASED ON GF IMAGERY INTERPRETATION
    QIAO Jun-xiang, SHI Feng, LI An, LI Tao, ZHANG Da, WANG Xin, Gesangdanzhen, SUN Hao-yue
    2025, 47(3):  789-805.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250043
    Asbtract ( )   HTML ( )   PDF (18102KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    On January 7, 2025, at 09:05, a MS6.8 earthquake occurred in Dingri County, Shigatse City, Xizang, China, with a focal depth of 10km and an epicenter located at 87.45°E, 28.50°N, as per data from the China Earthquake Networks Center. Rapid identification of coseismic surface ruptures is crucial not only for determining the seismogenic structure but also for post-earthquake damage assessment and emergency response. Based on the objective geological conditions of the Dingri earthquake-affected area, preliminary interpretation of pre- and post-earthquake high-resolution remote sensing data to delineate the spatial distribution and geometric characteristics of coseismic surface rupture zones before field investigations can effectively guide emergency response teams in rapidly and accurately identifying coseismic surface ruptures and conducting subsequent field surveys. This methodology demonstrates high feasibility and necessity for optimizing field workflow efficiency and ensuring targeted structural analysis of seismogenic faults.

    This study utilizes high-resolution imagery from the Gaofen -1 satellite to interpret pre- and post-earthquake images, rapidly obtaining the spatial distribution and geometric characteristics of the coseismic surface ruptures. The seismogenic fault is identified as the Dengmecuo Fault, located in the southwest segment of the Shenzha-Dingjie Rift. The coseismic surface rupture zone is primarily distributed near Gulong Village in the northern and central segments of the Dengmecuo Fault, with discontinuous extensions of approximately 15km, consistent with the location of the pre-existing fault. The coseismic surface rupture zone is tectonically partitioned into three distinct segments based on their geographic distribution: The Nixiacuo and Yangmudingcuo segments situated within the northern fault segment, and the Gurong segment developed along the eastern piedmont front of Gurong Village. These three segments exhibit marked differences in spatial scale and morphological characteristics between remote sensing observations and field investigations. The Nixiacuo segment exhibits linear and continuous coseismic surface ruptures, extending approximately 5km, with prominent linear traces visible in satellite imagery, facilitating clear identification. Field investigations reveal that this segment predominantly develops a series of large-scale tensional fractures and three categories of fault scarps with differential heights, with a maximum coseismic displacement of ~3m recorded. In contrast, both the Yangmudingcuo and Gurong segments exhibit smaller-scale coseismic surface ruptures localized along the eastern graben-bounding basin-range boundary. These secondary ruptures are characterized by minor extensional cracks with limited opening amounts(<1m vertical displacements)and fault scarps, manifesting as dark-gray linear features in remote sensing imagery that coincide with pre-existing rupture traces. Their partial obscuration in spectral signatures renders comprehensive visual interpretation impractical, necessitating field validation for complete delineation. Furthermore, this study identifies an ~10km-long associated surface deformation zone along the eastern Dengmecuo Lake shoreline, exhibiting a structural assemblage of extensional fissures proximal to the mountain front and compressional ridges adjacent to the lakeshore, with concomitant sand boil structures observed within the deformation zone. These extensional features present as dark-toned linear traces paralleling the main surface ruptures, displaying discontinuous arcuate configuration convex toward the mountain front within the mid-fan sector of alluvial fans. The compressional uplift zone west of the extensional belt appears as curvilinear bands with central whitish zones flanked by shadowed margins in imagery, demonstrating enhanced spatial continuity along distal fan margins.

    The conclusions of this study exhibit high consistency with previous research based on submeter-level resolution remote sensing imagery, confirming the reliability of remote sensing interpretation. However, interpretations derived from 2m resolution imagery have inherent limitations, including difficulties in identifying small-scale surface ruptures and distinguishing surface deformations of diverse genetic origins. Based on existing research integrated with remote sensing image interpretation and field investigations, the identification capability of coseismic surface ruptures in 2m-resolution remote sensing imagery is fundamentally governed by the three-dimensional geometric parameters of the rupture zone—specifically, rupture zone width, along-strike continuity length, and fault scarp vertical displacement. From the case study of the Dingri earthquake, we hypothesize that GF-1 satellite imagery with 2m resolution is capable of resolving coseismic surface ruptures characterized by a minimum rupture width of 1m, along-strike continuity exceeding 10m, and vertical offset≥1m Furthermore, the composite deformation pattern observed along the eastern shore of Dengmecuo Lake, characterized by rear tensile fractures and frontal thrusting, is interpreted as shallow detachment sliding induced by seismic shaking, representing a secondary associated surface deformation zone rather than direct fault displacement. Therefore, the identification of coseismic surface ruptures necessitates integrating remote sensing interpretation with field investigations, requiring not only the analysis of their geometric configuration and vertical displacement but also a synthetic evaluation of genetic origin based on the geological environment to achieve comprehensive and accurate determination.

    This study conducted a detailed interpretation of the coseismic surface rupture zone in the earthquake-affected area using GF-1 image data acquired within 24 hours post-earthquake, rapidly delineating the overall geometry and location of the coseismic surface rupture. This approach effectively supported subsequent field investigations and enhanced the efficiency of earthquake emergency response. Furthermore, rigorous field reconnaissance was carried out to validate the remote sensing interpretation results. The findings demonstrate a high consistency between the interpreted segments of the coseismic surface rupture and field observations, confirming the reliability of the remote sensing interpretation. This highlights the potential of domestic high-resolution satellite data for rapid coseismic surface rupture mapping and seismogenic structure identification, providing a feasible, rapid, and efficient methodology for future emergency response to major earthquakes.

    SOURCE CHARACTERISTICS OF THE 2025 DINGRI EARTHQUAKE AND ITS IMPLICATIONS FOR THE ACTIVITY OF THE SHENZHA-DINGJIE RIFT ZONE
    WAN Yong-ge, WANG Run-yan, JIN Zhi-tong, LAN Cong-xin
    2025, 47(3):  806-819.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250029
    Asbtract ( )   HTML ( )   PDF (6193KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    On January 7, 2025, a MS6.8 earthquake occurred in Dingri, Qinghai-Xizang Plateau. The earthquake occurred in the Shenzha-Dingjie Rift zone with a moment magnitude of 7.2, which is relatively rare for such a large earthquake to occur in the rift zone. To understand the source characteristics of the earthquake, the available seismic moment tensor solutions of the earthquake are collected, by averaging the corresponding moment tensor element, we obtained the central seismic moment tensor of the earthquake. By using the central seismic moment tensor to replace the central focal mechanism solution to understand the source characteristics, it is not only considering the results of seismic moment tensors from different sources but also a simpler algorithm than the previous central focal mechanism solution algorithms. By decomposing the central seismic moment tensor into a dislocation source part and a compensated linear vector dipole part, it was found that the dislocation source part occurs in a mechanical state of near vertical compression and near east-west tension. In contrast, the compensated linear vector dipole part exhibits a moment release mode of simultaneous vertical and north-south compression and east-west tension. A comprehensive analysis of previous geological surveys shows that the Shenzha-Dingjie rift zone is a steeply dipping normal fault. Therefore, we speculate that the non-double couple moment tensor of the Dingri earthquake is a comprehensive result of continuous sliding with a steep fault in the shallow crust and gentle low-dip in the deep part of the crust, which formed a shovel-shaped fault, and the sliding angles gradually change from shallow to deep. The focal mechanism of aftershocks of the Dingri earthquake and seismic moment tensor data of surrounding historical earthquakes are also collected. The comprehensive seismic moment tensors for aftershocks and historical earthquakes are obtained by summing elements of the seismic moment of every earthquake. The same analysis of the seismic moment of the mainshock was conducted with the comprehensive seismic moment tensors for aftershocks and historical earthquakes. It was found that the patterns of the dislocation source part and the compensated linear vector dipole part obtained were consistent with that of the main shock, which supported the analysis results of the source characteristics of the main shock. This is the first time that the total seismic moment tensor elements of earthquakes and aftershocks in geological fault zones have been averaged to study the properties or characteristics of fault zones or earthquake sequences. The obtained results are still encouraging. It provides a comprehensive method for analyzing fault zones or aftershock zones and analyzing fault properties or focal rupture characteristics. There are multiple hypotheses regarding the formation mechanism of the north-south rift zone on the Qinghai-Xizang Plateau. The analysis of the source characteristics of the Dingri earthquake supports the related models of magmatic activity/intrusion of the lower crust bottom splitting and simultaneous squeezing of Indian plate material. Based on the central seismic moment tensor solution of the Dingri earthquake and the geometry shape of the rift zone, it is inferred that the non-double-couple in the seismic moment tensor originates from the changes in bending faults and sliding angles, which provides ideas for intuitively explaining the non-double-couple part of the seismic moment tensor.

    RESEARCH ON EARLY AFTERSHOCKS OF THE 2025 DINGRI M6.8 EARTHQUAKE BASED ON THE DEEP-LEARNING-BASED SINGLE-STATION LOCATION METHOD
    ZHI Long-xiang, ZHAO Xu
    2025, 47(3):  820-834.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250040
    Asbtract ( )   HTML ( )   PDF (6328KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    This study adopts the deep learning-based “DiTing” single-station method for near-source earthquake detection and seismic phase picking. At the technical level of single-station localization, a novel multi-feature fusion approach is proposed to accurately estimate the epicentral distance. This method integrates the virtual wave velocity technique, the travel-time table method, and a multi-orbit surface wave approach. To determine the back azimuth, the methodology combines principal component analysis(PCA), the single-value moving average method(SV), the radial amplitude maximum method, and the surface wave polarization technique. Furthermore, it introduces a dynamic weighting mechanism that adjusts the contribution of each method based on their respective uncertainties and the reliability of the phase picking. This approach effectively addresses two longstanding challenges in traditional single-station localization: the 180° azimuthal ambiguity and high uncertainty in back azimuth estimation.

    Using continuous waveform data from the two nearest broadband seismic stations to the 2025 Dingri M6.8 earthquake in Tibet, the study conducted rapid analysis for the period from January 7 to 15, 2025. The analysis identified 2, 255 and 1, 730 aftershocks at each respective station, revealing the spatiotemporal characteristics of the aftershock sequence. The aftershocks display a predominantly north-south spatial distribution that aligns closely with the strike direction of the causative fault. The aftershock zone extends approximately 70 kilometers along the fault, primarily concentrated to the west of the Dengmecuo fault. Moreover, the spatial distribution of aftershocks shows a strong correlation with the pattern of mainshock co-seismic slip, with aftershocks clustering in regions of relatively low slip. This correspondence supports the hypothesis that stress redistribution following the mainshock governs aftershock occurrence. Comparison with aftershock catalogs produced by other researchers further confirms the consistency and reliability of the results obtained in this study.

    Extensive experimental results demonstrate that the integration of multiple localization algorithms significantly enhances the stability and accuracy of single-station solutions. The study establishes a cross-validation framework whereby results from different algorithms are compared, enabling the identification and elimination of erroneous data affected by noise or local anomalies. This approach substantially improves the robustness of single-station localization, particularly in complex seismic environments.

    The single-station localization technique allows for the rapid inversion of key earthquake parameters within seconds of an event, offering a substantial reduction in response time compared to traditional multi-station localization systems. This improvement translates into critical additional seconds for issuing warnings and initiating evacuations, thereby mitigating casualties and property loss in the affected regions. As seismic data processing techniques and single-station localization algorithms continue to evolve, further improvements in localization accuracy are anticipated. With the accumulation of large-scale seismic datasets and the integration of advanced intelligent computing technologies, the potential of single-station localization for near-field strong earthquake early warning is expected to be significantly expanded in the future.

    AFTERSHOCK PROBABILISTIC FORECASTING AND TESTING OF OPERABILITY IN EARTHQUAKE FIELD INVESTIGATION ON-SITE: A CASE OF THE 2025 DINGRI MS6.8 EARTHQUAKE IN XIZANG
    ZHANG Sheng-feng, ZHANG Yong-xian
    2025, 47(3):  835-849.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250028
    Asbtract ( )   HTML ( )   PDF (5197KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    On January 7, 2025, a MS6.8 earthquake struck Dingri, Xizang, causing significant economic losses and casualties. In response, the China Earthquake Administration launched a multidisciplinary scientific investigation, among which the analysis of sequence characterization and the probability forecasting of large aftershocks is an important and meaningful part of the work. This study aims to enhance the understanding of the aftershock sequence and provide timely scientific support for field investigations. To achieve this, we employ a temporal Epidemic-Type Aftershock Sequence(ETAS)model to perform a real-time tracking analysis of the aftershock sequence over the first seven days following the mainshock. The temporal ETAS model was employed to analyze the evolving characteristics of the aftershock sequence at 0.1-day intervals, and short-term aftershock probability forecasts were generated for the subsequent one-day period. Model performance was evaluated using the Brier Score, a metric that quantifies the agreement between probabilistic forecasts and observed aftershock occurrences. The evaluation focused on different magnitude thresholds to assess the consistency and predictive skill of the model.

    Key findings from our study include: 1)The fitting of ETAS model to observed aftershock activity was generally consistent with reality. The fitted model parameters suggest that the overall decay rate of aftershocks aligns closely with typical sequence decay behaviors(p=1.06). Moreover, the proportion of triggered ‘offspring’ events within the sequence is relatively low(α=1.58), indicating that off-spring events did not heavily dominate the primary aftershock activity. The model’s fitting results are consistent with the observed seismic sequence, except for a slight deviation identified around the 220th aftershock, where the observed activity exceeds the expectation based on a homogeneous Poisson process. 2)A time-tracking analysis of the model parameters across varying magnitude thresholds reveals that the parameter estimates begin to stabilize approximately 2.8 days after the mainshock. This suggests that incomplete aftershock recordings during the early phase can impact the reliability of early parameter estimation. Thus, early-stage catalog incompleteness should be carefully accounted for in operational forecasting models. 3)The model also demonstrates high sensitivity to the occurrence of strong aftershocks. When such events occur, they are quickly reflected in the intensity and frequency curves, demonstrating the model’s potential and strong applicability for short-term aftershock forecasting, particularly in a science-based emergency response context. 4)Brier score evaluation further supports the model’s forecasting effectiveness. For aftershocks above magnitude 3.5, 4.0, and 5.0, the forecasting performance consistently exceeds that of a random forecast baseline. Although the model underperforms slightly in forecasting aftershocks above magnitude 4.5 in the early stages, its performance improves over time, especially for magnitude 4.5 and 5.0 events, indicating increasing skill as more data accumulates. These findings highlight the potential of integrating Brier Score evaluation into the temporal ETAS model for assessing probabilistic aftershock forecasts.

    The results demonstrate that the ETAS model provides valuable operational forecasting capabilities for guiding scientific investigations and emergency response following major earthquakes. The study also identifies key challenges for future improvements, including data completeness, parameter stability, and model adaptability to complex sequences of aftershocks. Moving forward, further refinement of hybrid forecasting approaches—integrating multiple models based on statistical and physics-based methods—could enhance the accuracy and reliability of short-term aftershock forecasting. The operational feasibility of the ETAS model, combined with rigorous evaluation metrics, underscores its role in advancing earthquake forecasting methodologies and supporting earthquake disaster risk reduction in China and beyond.

    RELOCATION AND FORESHOCK SEQUENCE IDENTIFICATION OF DINGRI MS6.8 EARTHQUAKE IN XIZANG
    YIN Xin-xin, ZUO Ke-zhen, ZHAO Cui-ping, CAI Run
    2025, 47(3):  850-868.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250025
    Asbtract ( )   HTML ( )   PDF (7453KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    Using seismic data from the Tibet Regional Seismic Network between January 2021 and January 2025, we relocated 7, 951 earthquakes employing the double-difference algorithm. To ensure relocation reliability, we selected phase data with minimal travel-time residuals, constraining earthquake pairs to a maximum separation of 30km and requiring at least eight common phase arrivals per pair. This yielded 4, 370 high-precision relocations, with average relative errors of 0.130km(longitude), 0.131km(latitude), and 0.199km(depth). The relocated mainshock location is(28.501°N, 87.477°E)with a focal depth of 9.3km. The aftershock sequence extends approximately 70km in a nearly north-south direction, with depths mainly concentrated between 3 and 15km.

    We conducted a detailed analysis of the foreshock activity preceding the MS6.8 Dingri earthquake. Due to sparse station coverage near the epicenter, traditional seismic monitoring methods were insufficient for detecting small events. To address this, we applied the deep learning-based PhaseNet model to continuous waveform data from the nearest station(ZHF, ~50km from the epicenter), in combination with a single-station amplitude-magnitude empirical relationship for magnitude estimation. This approach significantly improved catalog completeness. Within the 56-hour window prior to the mainshock, we identified 90 seismic events, of which 80(88.9%)were microearthquakes with magnitudes ML<2.0. In contrast, the regional network recorded only 8 events in the same period. A reliable single-station magnitude calibration was established(log10A=0.77ML+1.36) using 293 aftershocks. For commonly detected events with ML<2.0, the average magnitude difference between the single-station and regional network methods was just 0.09, confirming the accuracy of the single-station approach. Based on the enhanced catalog and using the maximum curvature method accounting for magnitude uncertainty, the completeness magnitude was determined to be ML1.10. The b-value, estimated via the maximum likelihood method, was 0.58±0.07. Relocation results show that foreshocks were spatially clustered within the eventual aftershock zone, approximately 20km from the mainshock epicenter. Eight foreshocks occurred within the final hour before the mainshock, with the largest(ML3.7)occurring approximately one hour prior.

    These findings demonstrate that deep learning-based, single-station detection methods can substantially enhance earthquake monitoring in regions with sparse seismic networks. The spatial and temporal characteristics of the foreshock sequence offer critical insights into earthquake preparation processes. The single-station magnitude estimation method presented here provides a valuable reference for seismic monitoring in similarly data-limited regions.

    COMPARISON OF THE CHARACTERISTICS OF EARTHQUAKE SEQUENCE AND INTENSITY OF THE JANUARY 7, 2025 MS6.8 DINGRI EARTHQUAKE IN XIZANG
    WU Xiao-fei, MENG Ling-yuan
    2025, 47(3):  869-880.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250044
    Asbtract ( )   HTML ( )   PDF (7125KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    At 09:05 on January 7, 2025, a magnitude MS6.8 earthquake struck Dingri County, Shigatse City, in the Xizang Autonomous Region of China. The earthquake had a focal depth of 10km and is the largest event in the region in recent years, resulting in severe damage and a wide area of impact. The disaster caused over 120 fatalities, damaged approximately 27,000 dwellings, and led to the collapse of 3600 structures. This study presents a comprehensive analysis incorporating seismogenic structures, historical seismicity, earthquake sequence evolution, and intensity mapping.

    The epicenter is located within the Lhasa Block of the Tibetan Plateau, in the high-altitude valleys and basins north of the Himalayan mountain range, where elevations exceed 4,000m within a 10km radius. The tectonically active Lhasa Block includes multiple fault systems. The earthquake likely originated on the Dengmecuo fault, a segment of the Shenza-Dinggye Rift fault zone, which trends approximately north-south. Since 1900, the region within 300km of the epicenter has experienced 15 earthquakes of magnitude MS6.0-6.9, and notable larger events include the MS8.1 Nepal earthquake of April 25, 2015, and its aftershocks.

    According to the China Earthquake Networks Center, the event had a moment magnitude of MW7.1 and a centroid depth of 15km. The focal mechanism indicates normal faulting. The radiated seismic energy was approximately 1015 J. The moment and energy magnitudes both exceed the surface-wave magnitude, suggesting an efficient release of seismic moment and energy during rupture. The aftershock sequence extends approximately 75km in a north-south direction, aligned with the inferred fault strike. Aftershocks were concentrated near the mainshock and its northern extent, separated by a relatively quiet central segment.

    Intensity analysis indicates a maximum intensity of Ⅸon both the United States Geological Survey(USGS)simulated intensity map and the measured intensity map from the China Earthquake Administration. The meizoseismal area follows a north-south distribution consistent with the fault trend. The Ⅵ-degree isoseismal zone in the USGS map extends nearly 200km in a north-south orientation, while the measured map shows isoseismal lines trending NNE, with a long axis of about 191km and a short axis of 152km. The orientations of the isoseismal lines and the meizoseismal area in both maps are broadly consistent, though discrepancies exist in areal extent and coverage. Notably, the simulated intensity in the northern part of the aftershock zone is higher than that of observed. This discrepancy likely arises from the fact that the area is sparsely populated, with no significant structural damage or reported casualties.

    ANALYSIS ON THE EVOLUTION CHARACTERISTICS OF LOCAL STRESS FIELD IN THE MAGNITUDE 6.8 EARTHQUAKE SEQUENCE IN DINGRI, XIZANG
    WANG Peng, DAI Zong-hui, KONG Xue, LI Bo, XU Chang-peng, ZHANG Meng-xin
    2025, 47(3):  881-896.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250042
    Asbtract ( )   HTML ( )   PDF (5456KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    Understanding the stress evolution of earthquake sequences is critical for elucidating the physical mechanisms driving earthquake nucleation and rupture. This study investigates the spatiotemporal variations in the stress field following the Dingri MS6.8 earthquake in Tibet, using seismic data from both permanent and temporary mobile stations. Double-difference relocation was performed using HypoDD 2.1, incorporating phase data from temporary seismic networks for improved accuracy. Focal mechanism solutions were determined for 189 events with well-constrained P-wave first-motion polarities and adequate azimuthal coverage(≥8 observations), using the polarity method. The SATSI algorithm was subsequently applied to invert the orientations of the principal stress axes and estimate the stress ratio R.

    The relocation results indicate that the mainshock ruptured the southern segment of the Dengmecuo Fault, with aftershocks propagating northward along the fault’s N-S trending structure. The aftershock distribution reveals a westward-dipping fault geometry. In the central portion of the rupture zone, both eastward- and westward-dipping fault branches are present, while the southern segment exhibits intersecting NW- and NE-striking faults, suggesting multiple rupture planes. The mainshock likely occurred near the eastern boundary between the east- and west-dipping segments, consistent with surface ruptures observed in the field.

    Stress inversion results indicate a normal faulting regime. The maximum principal stress(σ1)has a trend of 142° and a plunge of 67°, while the minimum principal stress(σ3)trends at 110°(W)with a shallow plunge of 7°, and the intermediate stress(σ2)trends at 17° with a plunge of 21°. The optimal stress ratio(R=0.22)suggests a dominantly extensional regime, consistent with the regional tectonic setting of N-S compression and NE-SW extension. Temporally, the orientation of σ1 evolved from 135°(SSE)to 180°(S), and σ3 shifted from NEE to an EW orientation, reflecting a post-seismic adjustment toward a stable regime of NS compression and EW extension. The R-value initially decreased from 0.5 to 0.05, followed by a gradual increase to 0.25, indicating early release of horizontal extensional stress and an increasing influence of vertical σ1—typical of normal faulting sequences. Aftershock activity diminished within seven days and stabilized thereafter, indicating progressive dissipation of residual stress. Spatially, the source region was divided into southern, central, and northern clusters, with respective dominant strike orientations of NNW, NNE, and NNW. The southern cluster, which recorded the most events and the most diverse focal mechanisms, yielded well-constrained stress inversions(narrow confidence intervals). However, the plunge of σ1 in this zone was only 31°, deviating from the near-vertical orientation typical of pure normal faulting. This deviation likely reflects complex fault geometry and secondary fracturing, which may have induced localized strike-slip components. In the central and northern zones, σ3 remained horizontally oriented toward the SWW. In the northern cluster, σ1 rotated to a NE orientation, likely influenced by increased strike-slip activity near the Nongqu Fault. Zone Ⅱ exhibited unstable inversion results, with overlapping σ1-σ2 confidence intervals, indicating a more complex local stress field. A northward increase in R suggests a transition from dominantly extensional to more strike-slip-dominated deformation.

    The region remains in a phase of post-seismic stress adjustment and has not yet returned to its pre-mainshock stress state. Continued seismic monitoring, particularly of the structurally complex southern fault system and the northern strike-slip segments, is essential for assessing future seismic hazard and stress accumulation.

    RESEARCH ON GROUND MOTION SIMULATION OF THE DINGRI MS6.8 EARTHQUAKE IN XIZANG BASED ON DIFFERENT SOURCE MODELS
    YIN Xiao-fei, QIANG Sheng-yin, ZHANG Wei, SHAO Zhi-gang, WANG Wu-xing, YUAN Xiao-xiang, LI Yong-sheng, LIU Hao
    2025, 47(3):  897-916.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250038
    Asbtract ( )   HTML ( )   PDF (16398KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    On January 7, 2025, a magnitude MS6.8 earthquake occurred in Dingri County, Shigatse City, within the Xizang(Tibet)Autonomous Region. This normal-faulting earthquake struck the Lhasa Terrane in the southern Tibetan plateau, a region characterized by a series of nearly north-south trending normal faults and associated tectonic rift valleys—features indicative of the region’s ongoing extensional deformation and potential for future strong seismic events. Given the high seismic hazard in the southern Qinghai-Xizang Plateau and the area’s complex mountainous terrain, which increases the likelihood of secondary disasters such as landslides, assessing strong ground motion is crucial for linking and quantifying seismic hazard and risk. Accordingly, simulating strong ground motion for a hypothetical MS6.8 earthquake in this region holds significant practical value. Such analysis contributes both theoretical insight and practical guidance for regional seismic disaster prevention and mitigation.

    This study simulates the strong ground motion of the Dingri MS6.8 earthquake using two source models derived from joint inversion of InSAR coseismic deformation, teleseismic waveforms, and strong motion recordings. A three-dimensional curvilinear finite-difference method with curved grid meshing is employed to model the seismic wave propagation and ground motion characteristics. By comparing the results of the two source models, the spatial distribution of seismic ground motion and the underlying causative mechanisms are analyzed. The key findings are as follows:

    (1)Simulated ground velocity time histories at four near-field stations, processed with a 0.2Hz low-pass Butterworth non-causal filter, closely match the observed strong motion records, verifying the accuracy and reliability of the simulations.

    (2)Due to a NNE-directed unilateral rupture, peak ground velocities(PGVs)in the forward rupture direction(NNE)are significantly higher than those in the reverse direction(SSW), demonstrating a clear rupture directivity effect.

    (3)A comparison of PGV distributions across the fault shows that values on the upper plate(western side)are significantly higher than those on the lower plate(eastern side), indicating a strong upper-plate effect. Vertical surface displacements on the fault’s upper plate, as simulated by the two models, reach 2.0m and 2.1m, respectively—values that are in close agreement with field measurements from the Dingri earthquake geological survey.

    (4)Both source models simulate a maximum seismic intensity of Ⅸ, with high-intensity zones extending predominantly in the NNE direction. The simulated intensity distributions are generally consistent with field observations, though discrepancies exist in two areas: From northern Dingri County to southwestern Angren County, and in central-southern Gangba County. The intensity distribution produced by source Model 2 shows better agreement with the observed data.

    This study highlights the importance of using source models derived from the joint inversion of InSAR, teleseismic, and strong motion data—or a broader combination of geophysical constraints including GPS—to improve the accuracy of strong ground motion simulations. The results offer an important scientific basis for advancing our understanding of seismic wave propagation and strong ground motion characteristics associated with normal faulting earthquakes in the southern Qinghai-Xizang Plateau.

    PRELIMINARY SIMULATION OF LONG-PERIOD GROUND MOTION OF THE DINGRI MS6.8 EARTHQUAKE ON JANUARY 7, 2025
    JI Zhi-wei, YU Hou-yun, LI Zong-chao, JU Chang-hui, SUN Yao-chong, ZHANG Yong-xian, CHEN Xiao-fei
    2025, 47(3):  917-931.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250031
    Asbtract ( )   HTML ( )   PDF (8399KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    On January 7, 2025, a MS6.8 earthquake struck Dingri County, China. Strong earthquake observation records provide critical insights into ground motion, aiding in macro-intensity assessments, post-earthquake emergency responses, and loss estimations. These records help fill gaps in near-fault strong motion observations and capture the complexity of various near-fault vibrations associated with the earthquake source. Due to the lack of strong motion observation stations near the epicenter, no effective near-field long-period strong ground motions were obtained in this earthquake. The curved grid finite-difference method employs a traction mirror technique derived from the stress mirror method to process free-surface boundary conditions within a curved coordinate system accurately. This method has been extensively utilized for rapid earthquake disaster assessment and simulating strong ground motion. To evaluate the long-period ground motion and velocity pulse distribution of this earthquake, this study applies the curved grid finite-difference method, incorporating the strong earthquake rupture model of the Dingri earthquake and topographic data from the source area. The simulation results illustrate the wavefield propagation process and intensity distribution in the affected region. Furthermore, using a velocity pulse identification method, the study determines the velocity pulse distribution characteristics of the source area. The study accounts for the region’s undulating terrain by first linearly interpolating and downsampling the Shuttle Radar Topography Mission(SRTM)terrain data to align it with the computational grid. The velocity medium model, which significantly influences strong ground motion, is also interpolated and corrected to match the terrain, ensuring compliance with computational requirements. The accuracy and reliability of the simulation results are validated by comparing them with observed waveform and velocity wavefield data. The findings indicate that peak ground velocity(PGV)in the vertical(UD)component is significantly higher than in the east-west(EW)and north-south(NS). This phenomenon is attributed to the normal fault mechanism of the Dingri earthquake. Although some vertical ground motion records exist, near-fault vertical motion data remain scarce. Previous studies suggest that, in near-fault regions, the peak vertical acceleration-to-horizontal acceleration ratio is influenced by factors such as magnitude and epicentral distance, often exceeding the standard 2/3 ratio and sometimes surpassing 1. The maximum simulated intensity in this study is IX, with higher-intensity areas concentrated near the fault’s hanging wall, demonstrating a pronounced hanging wall effect. Due to local topographic influences, the intensity distribution appears irregular. However, the simulated intensity pattern aligns with observed intensity trends, confirming the validity of the long-period earthquake simulation results. Further analysis reveals that near-fault intensity distribution is closely linked to the rupture characteristics of the source area. In the hanging wall region, seismic wave propagation is significantly influenced by fault geometry and surrounding geological conditions. Additionally, the study indicates that intensity distribution varies considerably under different terrain conditions, particularly at the interface between mountains and basins, where seismic wave focusing may locally amplify intensity. The simulated velocity pulses of the EW, NS, and UD components primarily concentrate within the surface projection area of the fault. The EW and NS velocity pulse distribution ranges are narrower than that of the UD component. Since a normal fault caused the Dingri earthquake, velocity pulses induced by rupture directivity effects predominantly appear in the component perpendicular to the fault plane. Compared to other fault types, such as strike-slip faults, normal fault earthquakes are less likely to generate significant velocity pulses. Strike-slip and reverse fault earthquakes, in contrast, tend to produce stronger velocity pulses due to their rupture mechanisms. Normal fault earthquakes are relatively rare, and this event has heightened awareness of potential normal fault seismic hazards in the rift zone. Strengthening research on pulse-type ground motions in normal fault earthquakes is crucial for disaster mitigation. Future studies will collect geometric data of the causative fault and regional stress field information to conduct dynamic rupture simulations. Through numerical analysis, this research aims to further understand pulse-type ground motions in normal fault settings, particularly their spatial distribution and influence on source and site conditions. The findings will enhance our understanding of pulse-type ground motions in normal fault earthquakes and provide a scientific basis for assessing potential seismic impacts and developing disaster prevention strategies.

    COMPARATIVE STUDY ON BUILDING DAMAGE CAUSED BY THE 2025 MS6.8 EARTHQUAKE IN DINGRI, XIZANG, BASED ON REMOTE SENSING AND SEISMIC SIMULATION
    YUAN Xiao-xiang, LIN Xu-chuan, CHEN Zi-feng, ZHANG Jian-long, DOU Ai-xia, XIAO Ben-fu, DU Hao-guo, YU Si-han, DING Xiang, FANG Jie, WANG Shu-min
    2025, 47(3):  932-948.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250019
    Asbtract ( )   HTML ( )   PDF (10837KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    Building damage caused by the earthquake is a significant factor contributing to fatalities in destructive earthquakes. Accurate assessment results of building damage after earthquakes are of great significance for revealing the mechanism of building damage under destructive earthquakes, guiding the research and development of seismic reinforcement technology and post-disaster recovery and reconstruction. On January 7, 2025, the Dingri M6.8 earthquake in Tibet caused many building damages and casualties. To quickly evaluate the seismic damage of buildings in this earthquake and improve the scientificity and timeliness of various methods for assessing seismic damage in the emergency stage, a comparative study on the seismic damage of buildings in this earthquake was conducted based on remote sensing and seismic damage simulation.

    Firstly, the GF-2 remote sensing images and Beijing -3 images obtained on January 8, 2025, were quickly collected as the input of remote sensing building seismic damage information extraction. The GF-2 remote sensing images, collected from October 2024, and actual strong ground motion data from the two nearest stations were used for seismic damage simulation prior to the earthquake. Based on data preprocessing, the physical models of the buildings in the disaster area were quickly extracted using a combination of deep learning and human-computer interaction from pre-earthquake remote sensing images.

    Then, based on constructing the remote sensing structure type characteristics and seismic damage interpretation characteristics of typical buildings, combined with some risk census data and a small amount of field survey information, the physical models of buildings in the study area were adjusted and optimized, and the data obtained on site were used for verification. Based on this information and post-earthquake images, the rapid identification of building seismic damage was carried out using optical remote sensing. At the same time, based on the data from two measured strong earthquake stations, the rapid simulation of earthquake damage was carried out by using the urban seismic simulator(YouSimulator).

    Taking the residential area as the statistical unit, the remote sensing seismic damage index was calculated for the quickly acquired remote sensing building seismic damage. According to the regional approximation principle, the corresponding model was used to calculate the ground equivalent seismic damage index, and the remote sensing intensity was estimated. Referring to the seismic intensity evaluation standard, the intensity of the seismic damage simulation results was calculated. Finally, the results of different methods were compared with the formal intensity.

    The findings indicate that most buildings in the disaster area are dispersed along rivers and valleys in the form of zonal distribution, exhibiting a small but relatively concentrated at some local spatial locations. The structural types are predominantly civil and stone-wood structures, which exhibit pronounced vulnerability under this earthquake. The two methods demonstrated a certain degree of consistency in identifying buildings with high seismic damage ratings in intensity zones above Ⅷ, with an overall assessment error of less than 1 degree. However, in zones Ⅶ and below area, there was a specific error in the evaluation results. This observation indicates that, in the aftermath of a major seismic event, both methods can contribute to emergency response efforts at various stages post-event by providing rapid building damage assessment results, thereby serving as scientific references for earthquake emergency relief and disaster reduction.

    TEXTURE FEATURE DAMAGE DETECTION OF SINGLE BUILD-ING BASED ON DRONE IMAGES AFTER EARTHQUAKE: A CASE STUDY OF 2025 DINGRI MS6.8 EARTHQUAKE IN XIZANG, CHINA
    DU Hao-guo, ZUO Xiao-qing, LIN Xu-chuan, XIAO Ben-fu, LU Yong-kun, HE Shi-fang, ZHANG Fang-hao, YUAN Xiao-xiang, TAO Tian-yan, YE Yang, DENG Shu-rong, ZHAO Zheng-xian, XU Jun-zu, BAI Xian-fu, ZHANG Yuan-shuo, ZHANG Lu-lu
    2025, 47(3):  949-968.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250026
    Asbtract ( )   HTML ( )   PDF (11410KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    Earthquakes, as sudden-onset natural disasters with high destructive potential, not only result in significant casualties but also cause severe damage to infrastructure—particularly buildings—posing major challenges to post-disaster rescue and reconstruction efforts. In emergency response scenarios, the rapid and accurate assessment of building damage is a critical prerequisite for formulating effective rescue strategies and allocating resources efficiently. Traditional manual on-site investigation methods, however, present notable limitations. Disaster-affected areas often experience traffic disruptions and harsh environmental conditions, hindering timely access for investigators. Moreover, manual assessments are time-consuming and generally incapable of meeting the urgent demands of rescue operations within the critical 72-hour post-disaster window. Large-scale manual surveys also involve safety risks, potentially leading to secondary casualties. Therefore, the development of rapid, efficient, and accurate building damage assessment methods holds significant practical and strategic importance.

    In response to this need, the this study we proposes an innovative rapid assessment method for earthquake-induced building damage using unmanned aerial vehicle(UAV)imagery combined with machine learning algorithms. This method leverages the advantages of UAV remote sensing—such as high mobility, flexibility, and high spatial resolution—together with advanced image processing and machine learning techniques to enable intelligent identification and assessment of building damage. The study focuses on the MS6.8 earthquake that struck Dingri County, Xizang, using it as a case study to validate the proposed methodology through a structured technical workflow. The assessment framework comprises three key stages. First, an object-oriented remote sensing classification(OORSC) approach was used to extract individual building features from UAV imagery. By employing rule-based classification strategies, this method effectively eliminates background noise such as trees and roads. After morphological filtering, the completeness of building boundary extraction exceeded 95%, and hole-filling performance was markedly improved, ensuring high-quality input data for subsequent analyses. Second, the study focused on the extraction and optimization of surface texture features. Using algorithms such as the Gray-Level Co-occurrence Matrix(GLCM) and Local Binary Pattern(LBP), critical parameters—including contrast, entropy, and variance—were derived. Experimental data show that, on average, the contrast of collapsed buildings is 26%lower than that of intact buildings, while entropy and variance increase by 32% and 41%, respectively. These features provide robust quantitative indicators for identifying structural damage. Lastly, the study implemented a comparative experimental design incorporating four technical routes, systematically evaluating the performance of classification algorithms such as Support Vector Machines(SVM) and Neural Networks(NN).

    The results demonstrate that the neural network model integrating optimized texture features yields the best performance, achieving an overall accuracy(OA)of 91% and a Kappa coefficient of 0.8. Compared to models excluding texture features, the improvement is significant: the neural network model without texture features achieved an OA of 85% and a Kappa of 0.6, while the SVM-based approach achieved an OA of 82% and a Kappa of 0.6. The recognition accuracy by damage level further reveals that severely damaged buildings are most accurately identified(94%)due to their distinctive visual characteristics, followed by collapsed(87%)and moderately damaged buildings(80%). Misclassification of collapsed structures mainly stems from blurred textures in the imagery. These findings underscore the critical role of texture features in building damage identification and validate the proposed method’s effectiveness in supporting post-disaster emergency response.

    Despite the promising results, the proposed method has several limitations in practical application. At the technical level, complex environmental backgrounds—such as similar roof materials and shadow effects—can interfere with detection accuracy and demand high image quality. At the data level, a lack of sufficient real-time ground truth data may compromise model training accuracy. At the application level, the method’s capacity to detect complex damage types—such as internal structural failures—remains limited. To address these challenges, future research will focus on several directions. From a technical innovation perspective, advanced methods such as deep learning will be explored, particularly the use of three-dimensional convolutional neural networks(3D-CNNs)for capturing volumetric building features. In terms of data integration, the fusion of multi-source data—such as LiDAR point clouds, digital surface models(DSM), and thermal infrared imagery—will be pursued to build a multimodal feature fusion framework. Methodologically, transfer learning and data augmentation will be applied to enhance model generalizability, and adaptive algorithms will be developed to manage complex and dynamic disaster scenarios. On the application front, the establishment of a standardized sample library and evaluation system is proposed to support the broader deployment and engineering application of the method.

    The significance of this study is multidimensional. Theoretically, it introduces a novel approach to building damage identification based on texture features and machine learning, enriching the theoretical framework of remote sensing-based disaster assessment. Technologically, it develops a comprehensive UAV image processing and analysis pipeline, offering a replicable technical route for related research. Practically, the established system can be directly applied to post-earthquake emergency response, enhancing the efficiency and effectiveness of rescue operations. With continued technological advancement, the method holds potential for adaptation to other disaster scenarios, such as typhoons and floods, thereby contributing to integrated disaster risk reduction. Future work will continue to advance research in this area, targeting breakthroughs in key challenges such as multi-source data fusion and intelligent algorithm optimization, with the goal of advancing disaster assessment technologies toward greater intelligence and precision, ultimately contributing to the protection of lives and property and the promotion of sustainable development.

    RAPID ESTIMATION OF PARAMETERS FOR THE M6.8 EARTH-QUAKE ON JANUARY 7, 2025 IN DINGRI(XIZANG, CHINA) BASED ON DATA-DRIVEN METHODS
    ZHAO Qing-xu, RONG Mian-shui, ZHANG Bin, WANG Ji-xin, KONG Xiao-shan, LI Xiao-jun
    2025, 47(3):  969-983.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250014
    Asbtract ( )   HTML ( )   PDF (4319KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    On January 7, 2025, a magnitude 6.8 earthquake struck Dingri County, Shigatse City, Xizang, China, causing significant casualties and property damage. Rapid and accurate estimation of earthquake magnitude, instrumental intensity, and ground motion parameters is essential for seismic hazard assessment and emergency response. Magnitude, as a key indicator of the energy released by an earthquake, lays the foundation for preliminary disaster assessment. Instrumental intensity, calculated from the intensity of ground motion observed by instruments, can be used directly to determine the extent of damage and the severity of disasters. Ground motion parameters, such as PGA and PGV, are widely used in seismic design, disaster assessment, and seismic damage prediction and are important metrics for evaluating the impact of earthquakes on buildings and infrastructure.

    In this study, a data-driven multi-task joint estimation framework is proposed that combines the SeismNet model for rapid magnitude and instrumental intensity estimation with the CRAQuake model for rapid estimation of ground shaking parameters. The framework is applied to the January 7, 2025, Dingri earthquake of magnitude 6.8, where the magnitude, instrumental intensity, and ground motion parameters are estimated and analyzed in parallel. The study starts by filtering and processing the strong motion data obtained, and then estimates the magnitude, instrumental intensity, and ground motion parameters by parallel computation. The results show that: 1)The estimated magnitude provided by SeismNet is 6.17 when the seismic data are input for 3s. With the increase in seismic wave duration, the estimated magnitude gradually approaches the catalog value, and the estimated magnitude is 6.71 at 7 seconds, with a significant reduction in the error. 2)For the instrumental intensity estimation, the results obtained by SeismNet are the same as those of the instrumental intensity flash report when the seismic data are input for 8 to 10 seconds. When data of 6 seconds or longer were used, there were no false alarms or omissions, showing a high degree of accuracy. 3)The estimates of ground-motion parameters provided by the CRAQuake model are in good agreement with observations, providing reliable results within a few seconds, especially for PGA, PGV, and other parameters, with minor and stable errors.

    These results indicate that the data-driven estimation model exhibits strong generalization ability in the Dingri earthquake, particularly in the epicenter region and the early post-earthquake period, providing fast and reliable decision support. With the increase of seismic wave duration, the estimation results of SeismNet and CRAQuake are more stable, the errors are gradually reduced, and the estimation accuracy is significantly improved. Through parallel computing, these two models can estimate multiple seismic parameters at the same time, which not only enhances the estimation efficiency, but also provides efficient and comprehensive technical support for earthquake emergency response. Additionally, data-driven methods offer significant advantages in earthquake emergency response, particularly in large-scale earthquake disasters. These methods can quickly estimate magnitude, instrumental intensity, and ground motion parameters, providing more accurate decision support. The results offer new technical insights and methodological support for future large-scale earthquake emergency response and lay the foundation for the widespread application of data-driven methods in the earthquake field.

    STUDY ON THERMAL INFRARED ANOMALIES OF THE 2025 DINGRI MS6.8 EARTHQUAKE AND SEVERAL EARTHQUAKE CASES IN SOUTHERN XIZANG
    ZHANG Li-feng, ZHONG Mei-jiao, PAN Yu-hang, GUO Ying-xia, SUN Xi-hao, ZHANG Yuan-sheng
    2025, 47(3):  984-998.  DOI: 10.3969/j.issn.0253-4967.2025.03.20250010
    Asbtract ( )   HTML ( )   PDF (9350KB) ( )  
    Figures and Tables | References | Related Articles | Metrics

    On January 7, 2025, a magnitude 6.8 earthquake struck Dingri County in Xizang. To investigate pre-seismic signals, we applied the relative power spectrum variation method and analyzed brightness temperature data from the FY geostationary meteorological satellite to detect thermal infrared(TIR)anomalies preceding the event. Building on this approach, the present study offers a more in-depth analysis of TIR anomalies prior to both the Dingri earthquake and a subsequent MS5.5 earthquake that occurred in Maduo on January 8, 2025. The Dingri MS6.8 earthquake occurred in the southern Tibetan Plateau, a region situated at the forefront of the ongoing collision and compression between the Indian and Eurasian plates. This area features major tectonic structures, including the Himalayan Frontal Thrust, the large-scale Karakoram-Jiali strike-slip fault, and seven nearly north-south-trending rift valleys that developed between them. The Dingri earthquake, a normal-faulting event, occurred within one of these rift valleys. To further understand TIR anomalies in this tectonically active region, we selected three additional earthquakes with similar geological settings and magnitudes for comparative analysis.

    The TIR anomalies associated with the Dingri earthquake were primarily distributed within the region bounded by multiple faults, covering a maximum area of approximately 210 000km2. The anomaly persisted for 80 days without complete dissipation and gradually evolved into a localized, high-intensity anomaly migrating in a northeastward direction. The Maduo MS5.5 earthquake occurred at the edge of this localized anomaly one day after the Dingri event. These two stages—widespread anomaly and localized concentration—are interpreted as part of a continuous anomaly evolution process, with the anomaly migration direction pointing toward the epicenter of the Maduo earthquake. Analysis of the time-series relative power spectrum prior to the Dingri earthquake revealed three significant episodes where the anomaly amplitude exceeded six times the background level. The first two episodes lasted 15 days and 22 days, respectively, while the third, which immediately preceded the Dingri event, persisted for 54 days, indicating a marked difference in duration and intensity. The relative power spectrum peaks were 12.4 for the Dingri event and 12.9 for the Maduo event, occurring 123 and 111 days, respectively, prior to the earthquakes.

    The spatial distribution of TIR anomalies associated with multiple earthquakes in southern Xizang appears closely linked to the extensional rift systems and active tectonic structures of the region. The directional evolution of these anomalies correlates with the eventual earthquake epicenters, which were generally located at the leading edge of the migrating anomaly zones. This finding is consistent with previous studies that have observed similar migration characteristics of TIR anomalies preceding earthquakes. Among the four examined earthquakes in southern Xizang, relative power spectrum peaks ranged from 12 to 18 times of the background level, appearing 28 to 123 days prior to the events. The maximum extent of anomalous areas varied between 170 000 and 210 000km2, with the duration of days exceeding the sixfold threshold ranging from 34 to 54 days. Despite some variation in these parameters, all events displayed common features of high-amplitude, large-area, and persistent anomalies, predominantly occurring during the short-term and imminent pre-seismic periods. Notably, in all four cases, no significant anomalies were observed directly at the epicenters; instead, the epicentral locations were consistently positioned at the margins of pronounced anomalous zones.

    The tectonic regime of the southern and central Tibetan Plateau is characterized by east-west extensional stress, resulting from the regional compressive stress field. This has led to the development of numerous north-south-oriented normal faults, which act as conduits for the upwelling of geothermal fluids. Additionally, the region experiences intense hydrothermal activity and significant CO2 degassing. Drawing on previous research, we propose that the combination of extensional rifting and active hydrothermal systems facilitates the ascent of geothermal fluids and greenhouse gases(including CO2)to the surface. This process likely contributes to the enhanced surface thermal radiation observed in satellite data and may explain the large-scale, fault-aligned TIR anomalies detected prior to these earthquakes in southern Xizang.