Most Read articles

    Published in last 1 year |  In last 2 years |  In last 3 years |  All

    All
    Please wait a minute...
    For Selected: Toggle Thumbnails
    STUDY ON SEISMOGENIC TECTONICS OF THE 2025 MYANMAR MS7.9 EARTHQUAKE
    XU Bin-bin, ZHANG Yi-peng, LU Le-jun, TIAN Qing-ying, YANG Xue, WANG Yang, ZHANG Pei-zhen
    SEISMOLOGY AND GEOLOGY    2025, 47 (2): 649-670.   DOI: 10.3969/j.issn.0253-4967.2025.02.20250089
    Abstract877)   HTML48)    PDF(pc) (13633KB)(358)       Save

    According to the China Earthquake Networks Center, an MS7.9 earthquake(hereafter referred to as the Myanmar earthquake)struck the Mandalay region of Myanmar(21.85°N, 95.95°E)on March 28, 2025, at a focal depth of 30km. The earthquake occurred along the central segment of the Sagaing Fault and was characterized by a right-lateral strike-slip rupture, generating a ~350km-long surface rupture zone with a maximum coseismic horizontal displacement of 6 meters. The event caused extensive damage to buildings and varying degrees of destruction to infrastructure, including roads and bridges.
    Situated in a critical tectonic region where the Indian Plate obliquely converges with the Eurasian Plate, the Myanmar earthquake offers valuable insights into plate boundary deformation processes. Detailed analysis of this event enhances our understanding of the deformation mechanisms along the Myanmar plate boundary and provides essential constraints for seismic hazard assessment along the southeastern margin of the Eurasian Plate. This research holds scientific significance for elucidating continental lithospheric deformation in response to oblique plate convergence. The findings contribute to regional early warning strategies and disaster mitigation efforts and offer a valuable reference for seismic risk studies in comparable tectonic settings worldwide.
    This study integrates Global Navigation Satellite System(GNSS)data from across the Sagaing Fault region, establishing a comprehensive GNSS velocity field for Myanmar and addressing previous gaps in coverage along the fault’s southern segment. Using multiscale spherical wavelet analysis and GNSS velocity profiles, we examine the deformation characteristics of the region. We calculate the slip rate deficit distribution along the Sagaing Fault and assess postseismic Coulomb stress changes. Combined with historical seismicity data, we investigate the seismogenic structure and stress perturbations in surrounding areas. The key findings are as follows:
    (1)The Myanmar MS7.9 earthquake was a right-lateral strike-slip event along the central Sagaing fault. The region is affected by the northeastward oblique convergence of the Indian Plate and southeastward extrusion of crustal material from the Tibetan plateau, resulting in strong north-south shear and east-west shortening. The Sagaing fault accommodates most of this deformation, with a rapid right-lateral slip rate of approximately 21~22mm/a.
    (2)High-resolution GNSS velocity profiles indicate significant fault locking at depths of 15~25km along the Sagaing Fault. The slip rate deficit analysis reveals a high locking ratio across the fault, indicating elevated seismic potential. Notably, the central segment shows lower seismic moment accumulation compared to the northern and southern segments, forming a ~300km-long seismic gap since 1900, capable of generating earthquakes exceeding magnitude 7.5.
    (3)Coulomb stress modeling suggests that the earthquake significantly altered the regional stress field. Stress accumulation zones were identified at both ends of the Sagaing fault and in the central Shan Plateau to the east. These regions of increased stress transfer and loading exhibit heightened potential for future large earthquakes, underscoring the need for enhanced seismic monitoring.

    Table and Figures | Reference | Related Articles | Metrics
    CRUSTAL DEFORMATION CHARACTERISTICS AND PROBABILITY PREDICTION OF STRONG EARTHQUAKE RISK IN XINJIANG AND ITS ADJACENT REGION
    CHEN Chang-yun, YIN Hai-quan
    SEISMOLOGY AND GEOLOGY    2025, 47 (2): 384-404.   DOI: 10.3969/j.issn.0253-4967.2025.02.20240151
    Abstract653)   HTML38)    PDF(pc) (9921KB)(205)       Save

    Based on the GNSS velocity field, we analyzed the present-day crustal deformation characteristics of Xinjiang and its adjacent regions before the January 23, 2024, Wushi earthquake using the spherical least squares configuration method. Based on the fundamental concept of active blocks, Xinjiang and its neighboring areas were divided into 17 active blocks by integrating regional seismic geological data. The slip rates of the boundary fault zones of these active blocks were then calculated using a three-dimensional elastic block model. Based on this block delineation, the study area was further divided into 91 potential seismic hazard zones. We incorporated geodetic observations, including fault strike-slip rates and dilatational strain rates derived from the GNSS velocity field, into the conventional probabilistic forecasting of significant seismic hazards. The relationship between the Wushi 7.1 earthquake and the regional crustal deformation characteristics, as well as the prediction results of the probability of strong earthquakes, was comprehensively analyzed.
    The direction and magnitude of GNSS velocity field motions in the study area under different dynamical backgrounds are distinctly characterized by zoning. GNSS stations east of longitude 76° move toward the north or northeast, while those on the west side move north or northwest. The velocity field difference is primarily evident on the north and south sides of the Tianshan Mountain. The velocity change of GNSS sites from the Tibetan plateau to the Tarim Basin is minor, but it significantly decreases after crossing the Tianshan Mountain, indicating that the Tianshan Mountain tectonics absorb most of the remote effects of the Indo-Eurasian plate collision. The principal strains obtained from the least-squares configuration results reveal that the extrusion characteristics are most prominent near the western section of the Southern Tianshan Mountains and the Altyn Tagh fault zone. In the western section of the Southern Tianshan Mountains, the direction of the principal compressive strains is perpendicular to the tectonic direction in the region. This suggests that the region is primarily influenced by a force perpendicular to the tectonic direction, resulting in the predominant retrograde movement of major fault zones in the area. Moreover, in addition to the western section of the Southern Tianshan Mountains, the region with a greater main compressive strain is the Altyn Tagh fault zone. The direction of the main compressive strain intersects obliquely with the Altyn Tagh fault zone, suggesting that the force background is linked to the left-lateral reverse slip movement of the Altyn Tagh fault zone.
    The results derived from the inversion of the three-dimensional elastic dislocation model reveal that the motion features of the main active faults within the study area indicate a predominance of dextral slip along the northwest-trending fault zones in the Tianshan region, while sinistral slip motion is primarily observed along the northeast-trending or northeastern fault zones. Apart from the Altyn Tagh fault zone, the strike slip rate of the dextral-slip fault zones is notably higher than that of the sinistral-slip fault zones. The fault zones at the northern edge of the Junggar Basin and the major faults in the Tianshan region are primarily characterized by extrusion movements. The extrusion rate of the fracture zones in the South Tianshan Mountain is higher than that in the North Tianshan Mountain. Specifically, the west section of the Keping fault zone and the west section of the Nalati fault zone exhibit the most prominent extrusion movements.
    The study area has been divided into 91 potential seismic hazard zones based on block delineation and the findings from previous research. The slip rates and regional surface strains of the major faults obtained from the inversion are utilized in classical probabilistic predictions to derive quantitative results regarding the strong earthquake hazard in the study area over the next 50 years. The results indicate that the strong earthquake hazard is primarily concentrated in the western section of the Southern Tianshan region. This includes areas such as the north-east or nearly east-west-trending Maidan fault, the Nalati fault, and the Usun Ridge Fault, as well as the north-west-trending Talas-Fergana fault zone and the north-west section of the Kyzyltau fault zone. Furthermore, the probability of strong earthquakes is elevated in the northern Luntai fault compared to the surrounding faults. In the northern Tianshan region, areas with relatively high probabilities of strong earthquakes include the Fukang fault and the western section of the Bogda fault.
    The seismic mechanism solution reveals that the M7.1 Wushi earthquake was a thrust earthquake, aligning with the characteristics of the Maidan fault zone. In this zone, the seismogenic fault is primarily influenced by extrusion motion. The Wushi earthquake occurred in the Maidan fault zone, situated at the border of the high shear strain rate zone and the high probability hazard zone. This occurrence validates the effectiveness and accuracy of the probabilistic prediction method.

    Table and Figures | Reference | Related Articles | Metrics
    EARTHQUAKE SEQUENCE RELOCATION AND SEISMOGENIC STRUCTURE OF THE 2024 MS7.1 WUSHI EARTHQUAKE ON JANUARY 23, 2024, XINJIANG
    WANG Xue-zhu, WU Chuan-yong, LIU Jian-ming, ZANG Ke-zhi, YUAN Hai-yang, GAO Zhan, ZHANG Jin-shuo, MA Yun-xiao
    SEISMOLOGY AND GEOLOGY    2025, 47 (2): 488-506.   DOI: 10.3969/j.issn.0253-4967.2025.02.20240153
    Abstract568)   HTML33)    PDF(pc) (11676KB)(154)       Save

    At 2:09 on January 23, 2024(Beijing time), an earthquake of MS7.1 occurred in Wushi County, Aksu region, Xinjiang Uyghur Autonomous Region. This earthquake is the first earthquake with a magnitude greater than 7 in the Tian Shan seismic belt after the 1992 Susamer M7.3 earthquake. This earthquake occurred near the range-front Maidan fault which is the boundary structure between the southern Tian Shan and Tarim Basin. In fault geometry, the Maidan fault shows a complex geometry. On the plane, the Maidan fault comprises multiple roughly parallel secondary faults. On the profile, the gently dipping mountain front fault controls the folding deformation of the late Cenozoic. The front and back mountain faults converge in the deep part. This causes the Tarim Plate to deeply subduct along a large subduction zone under the Tianshan. The thrust-nappe structural system is often composed of multiple fault zones that can generate strong earthquakes because these fault zones converge to the same detachment surface in the deep part. Because a fault zone triggers or suppresses strong seismic activity in adjacent faults after a strong earthquake.
    Because the main earthquake did not generate obvious coseismic surface ruptures, there is still a great controversy and uncertainty on its seismogenic structure. Determination of the seismogenic structure is crucial for analyzing the potential location of the next major earthquake. It is of great significance for evaluating the future risk of strong earthquakes in the region and the stress loading and triggering relationships between different faults.
    According to the analysis of Google’s high-precision satellite image interpretation results and field inspections, the geometric structure of Maidan is more complicated. The fault is roughly near the ancient stream of Yu Shangu Xi River, and it can be divided into two branches, east and west of the left steps. The faults all show clear signs of activity at the surface. In this study, we utilize the earthquake relocation results to determine the seismogenic structure of the 2024 Wushi event. Our results show that the main shock began to break in the deep and the aftershocks are extended from deep to the part. The deeper focal depth may be an important factor in preventing the coseismic sliding of the Wushi MS7.1 earthquake from being transmitted to the surface. The seismic sequence is exhibited by the northeast-southwest. The long axis direction is about 55°, and the total length is about 85km. The aftershock sequence is divided into the northeast, middle section, and southwest section. The rupture range of the Wushi MS7.1 earthquake is about 35km, which is the middle section of the aftershock sequence. The aftershocks in the northeast section are mainly distributed along F1-1, the aftershocks in the southeast section are mainly distributed along F2-2 and the aftershocks in the middle section are mainly distributed along F1-2. The aftershocks of different sections are distributed on different branches, which means that the strong earthquake triggered the adjacent seismic activity, which belongs to a more complicated grade joint rupture earthquake. At the same time, the CAP waveform counter and method was adopted to obtain the focal mechanism. The nodal planes parameters of the best double-couple focal mechanisms are: strike 115°, dip 52° and rake 132° for nodal plane I, and strike 240°, dip 54° and rake 49° for nodal plane Ⅱ, the depth of the centroid is about 17km. Based on the inversion results of the focal mechanism and the spatial distribution characteristics of the earthquake sequence, we believe that nodal plane II is the seismogenic fault plane. Based on the spatial distribution characteristics of the earthquake sequence, the focal mechanism solution, and the geological structure data of the earthquake area, we suggest that the Biedieli-goukou fault is the piedmont branch of the Maidan fault. The upper part of the fault(4~5km)has no rupture, and it still has a strong risk of strong earthquakes in the future. The Wushi MS7.1 earthquake triggered the Biedieli fault and the Aheqi fault. During the Late Quaternary period, these faults have repeatedly ruptured to the surface of strong earthquake incidents.

    Table and Figures | Reference | Related Articles | Metrics
    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
    SEISMOLOGY AND GEOLOGY    2025, 47 (3): 671-688.   DOI: 10.3969/j.issn.0253-4967.2025.03.20250030
    Abstract552)   HTML48)    PDF(pc) (7098KB)(344)       Save

    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.

    Table and Figures | Reference | Related Articles | Metrics
    SEDIMENTARY RECORDS AND SOURCE ANALYSIS OF A TSUNAMI EVENT ABOUT ~1 000 YEARS AGO IN THE PEARL RIVER ESTUARY ALONG THE COAST OF SOUTH CHINA
    WANG Wei-tao, YANG Xiao-qiang, SHU Peng, ZHANG Yu-hao, LIANG Hao, LI Lin-lin, LI Zhi-gang, WANG Da-wei, ZHANG Pei-zhen
    SEISMOLOGY AND GEOLOGY    2025, 47 (4): 999-1019.   DOI: 10.3969/j.issn.0253-4967.2025.04.20240170
    Abstract549)   HTML25)    PDF(pc) (7174KB)(150)       Save

    Since the Cenozoic time, the South China Sea(SCS)has formed one of the largest semi-enclosed marginal basins along the East Asian continental margin through the geological processes such as South China Sea plate rifting, seafloor spreading, and plate subduction. In the South China Sea Basin and its surrounding regions, a series of active geological structures have developed, for example, the Manila subduction zones, the Littoral(Binhai) fault zone, and the Continental Slope(Lupo) fault zone. The activity of these tectonic zones is highly prone to triggering extreme natural disasters such as earthquakes and tsunamis. Along the coastal zone of the northern part of the South China Sea(the South China continental margin), there are densely populated large cities with critical infrastructure, which are also regions severely affected by extreme natural disasters like earthquakes and tsunamis. Therefore, identifying the sedimentary records of large-scale paleotsunami events in the northern part of the South China Sea and analyzing their potential triggering mechanisms are of great significance for seismic-tsunami hazard assessment.

    This study focuses on the sedimentary strata of boreholes E15 and E12 from the Pearl River Estuary along the South China coast. Based on AMS 14C dating and a comparison of magnetic susceptibility data between boreholes E15 and E12, a high-precision chronological sequence was established for the core of borehole E15, spanning approximately the last 1000 years to the present. The core E15 is 5.9m long, with a progressively younger AMS 14C age sequence in the upper part of the core section from 4.24~0m. However, AMS 14C ages of the sediments in the lower part of core E15, from 5.9 to 4.24m, are sometimes reversed. The reversal ages may be attributed to the reworking or recycling of the sediments in the lower part of core E15.

    To reveal the depositional processes of borehole E15, we conducted detailed analyses of sedimentary grain size, sedimentary color, and geochemical elemental composition. The lower part of the core section(5.9~4.24m) for E15, consists of dark gray to grayish-black medium-to-coarse sand layers with poor sorting and contains abundant marine biodetritus. In contrast, the upper part of the core section(4.24~0m) for E15 is composed of dark gray to grayish-brown silty mud, fine sandy silt, and delicate sand layers. The upper core section exhibits finer grain size, lighter color, faint horizontal bedding, and higher terrestrial-derived elemental content, representing typical delta-shallow marine depositional environment. The lower part core section is characterized by a coarser grain size(medium to coarse) of sands, which lack clear sedimentary structures and exhibit higher offshore marine-derived elemental content, but relatively lower terrestrial-derived elemental content.

    Based on the sedimentary features and geochemical composition, the sands from 5.9~4.24m within the borehole E15 are completely different from the overlying normal, typical shallow sea-delta sediments. Considering the reversal AMS 14C ages, coarser grain size, poor sorting, darker color, higher offshore marine-derived elemental content, and lower terrestrial-derived elemental content in the lower part core of the E15, we propose that the sand layers with abundant marine biodetritus in the lower part of boreholes E15 (5.9~4.24m) may be deposits from an extreme hydrological event occurring approximately 1000 years ago. In fact, tsunami deposits dating back about 1000 years have been widely documented along the northern and western coasts of the South China Sea, the inner islands of the South China Sea, and the northwestern Philippines. Therefore, we suggest that the event deposits in the Pearl River Estuary region, along the northern part of the South China Sea, at ~1000 years ago, may also be the result of a tsunami event.

    Combining sedimentary evidence and numerical simulations, we hypothesize that a strong submarine earthquake may have occurred along the Manila subduction zone in the eastern South China Sea approximately 1000 years ago, triggering a large-scale tsunami. The medium and coarse-grained sand layers in the lower part (5.9~4.24m) of the E15 borehole within the Pearl River Estuary may be the consequence of this tsunami event.

    Table and Figures | Reference | Related Articles | Metrics
    PRELIMINARY STUDY ON THE SEISMOGENIC STRUCTURE OF THE 2024 M5 EARTHQUAKE SWARM IN YULI, TARIM BASIN
    LI Jin, DENG Ming-wen, ZHANG Zhi-guang, SUN Ye-jun, YAO Yuan, XU Kai-chi
    SEISMOLOGY AND GEOLOGY    2025, 47 (2): 463-487.   DOI: 10.3969/j.issn.0253-4967.2025.02.20250064
    Abstract519)   HTML22)    PDF(pc) (9183KB)(144)       Save

    The Tarim Basin, a typical small cratonic basin and the largest inland basin in China, is bordered by three major orogenic belts: the Kunlun, Altun, and Tianshan Mountains. Prior to 2000, seismic activity within the basin was sparse. However, with the ongoing development of oil and gas exploration and extraction after 2000, seismicity began to increase, forming a distinct pattern of earthquake swarms—referred to as the Tarim earthquake swarm. Notably, since 2017, both the frequency and magnitude of earthquakes have risen significantly. Major events include the Kuqa M5.6 earthquake on January 16, 2020, the Shaya M6.1 earthquake on January 30, 2023, and three M5.0 earthquakes in the Yuli-Kuqa area in 2024. The rising seismicity has drawn increasing attention, particularly due to heightened activity in the Yuli region since June 2024, now considered the most seismically active area in both Xinjiang and China.
    This study applies the Cut and Paste(CAP)method to invert focal mechanisms for three M5.0 earthquakes and other events above ML4.0 within the Yuli earthquake swarm. The double-difference relocation method was employed to refine the hypocentral locations of earthquakes above ML2.0. Furthermore, by analyzing the spatial distribution of smaller events and modeling the regional stress field, the study infers plausible parameters for seismogenic fault planes and investigates the faults responsible for the swarm. Among the 44 earthquakes analyzed, 42 exhibit strike-slip focal mechanisms. Cluster analysis reveals the presence of nearly north-south(NS)oriented nodal planes, consistent with regional strike-slip faulting trends in the NS to NNE direction. The relocated seismic sequence shows a clear division into eastern and western branches, both exhibiting linear NNE-trending patterns. Initially, seismic activity was concentrated in the eastern segment, but both branches demonstrated a northward migration over time. Based on the focal mechanisms, relocation data, and inferred fault plane parameters—combined with known fault structures in the region—it is preliminarily concluded that the Yuli earthquake swarm occurred along two nearly parallel, NNE-trending, right-lateral strike-slip faults: F1(east)and F2(west). Both faults are nearly vertical, with F1 exhibiting a steeper dip. The fault depths are estimated at 20~25km. The previously mapped faults F18 and F16 near the Yuli swarm may correspond to the shallow and deep sections of faults F1 and F2, respectively. The eastern fault(F18/F1) maintains consistent dip characteristics from shallow to deep levels, whereas the western fault(F16/F2)displays a gentler dip at depth compared to the surface.
    In addition to the Yuli swarm, the Tarim Basin has recently experienced the 2012 Luopu M6.0 earthquake(thrust)and the 2023 Shaya M6.1 earthquake(strike-slip), reflecting diverse faulting mechanisms driven by a common regional tectonic regime. The Luopu event resulted from thrusting of the Bachu Uplift toward the Awati Depression, induced by the ongoing convergence of the Indian and Eurasian plates. In contrast, the Shaya earthquake and the Yuli swarm likely represent secondary strike-slip structures accommodating crustal deformation from northward-directed compressive stress associated with NW- and EW-trending thrust belts in the basin interior. The Yuli earthquake swarm is located in the Yuke region of eastern Tarim, an area rich in oil and gas resources. In recent years, extraction activities have intensified, raising questions about a potential link between these activities and the increased seismicity. Further investigation is required to clarify the relationship between industrial activity and the occurrence of the Yuli earthquake swarm.

    Table and Figures | Reference | Related Articles | Metrics
    ARCHAEO-SEIMSIC INVESTIGATION REVEALS A DESTRUC-TIVE EARTHQUAKE OCCURRED IN THE HELUO REGION DURING THE HAN DYNASTY
    HU Xiu, LU Peng, WANG Hong-chi, FU Long-teng, MO Duo-wen, LI You-li, ZHANG Pei-zhen, ZHANG Hui-ping, WANG Zhi-shuo, HUI Ge-ge, CHEN Pan-pan, GUO Ai-lun, LUO Quan-xing, ZHAO Xian-gang
    SEISMOLOGY AND GEOLOGY    2025, 47 (4): 1020-1035.   DOI: 10.3969/j.issn.0253-4967.2025.04.20240126
    Abstract517)   HTML33)    PDF(pc) (8253KB)(68)       Save

    The Heluo region, the cradle of ancient Chinese civilization, has historically served as China’s political, economic, and cultural center. Historical records suggest that earthquake hazards frequently affected this region, particularly during the Eastern Han Dynasty. However, the political motivations of historical record-keeping raise questions about the reliability of these accounts. The buried nature of active faults in the plains complicates the identification of earthquake sources and magnitudes. In areas with abundant liquefaction deposits from historic and prehistoric earthquakes, archaeo-seismic investigations provide crucial information about past seismic activity by identifying liquefaction features related to strong ground motion.

    In this study, detailed archaeo-seismic investigations were conducted at the Gucheng archaeological site, located in the Jialu River floodplain(a sub-channel of the Huai River)in central-west Zhengzhou city. Observed liquefaction features above the paleo-cultural surface include sand dikes and sand blows formed by the upward flow of water and entrained sediments. Using AMS 14C and archaeological dating, this study confirms a catastrophic earthquake occurred during the Han Dynasty. By analyzing coseismic deformation types and distribution, and comparing with previous regional paleo-liquefaction studies, constraints on earthquake sources and magnitudes were established.

    The evidence primarily attributes the 119AD Luoyang earthquake to the Fengmenkou-Wuzhiling fault as the most likely seismogenic source, with a minimum magnitude of MS6.8 based on worldwide liquefaction-magnitude relations. This study demonstrates the potential of archaeo-seismic methods to provide reliable insights into prehistorical and historical earthquake hazards, even with limited evidence. The approach of reconstructing regional paleo-seismic events from liquefaction deposits is broadly applicable to basin zones with buried faults worldwide, particularly in areas with fluvial and lacustrine sediments. Additional archaeo-seismic research may further aid in regional seismic risk assessment and evaluating societal impacts.

    Table and Figures | Reference | Related Articles | Metrics
    STUDY ON FAULT DEFORMATION CHARACTERISTICS AND SEISMIC HAZARD IN THE KEPING THRUST BELT
    ZHU Shuang, GUO Nan-nan, PANG Ya-jin
    SEISMOLOGY AND GEOLOGY    2025, 47 (2): 448-462.   DOI: 10.3969/j.issn.0253-4967.2025.02.20240147
    Abstract497)   HTML18)    PDF(pc) (6489KB)(141)       Save

    As one of the most seismically active regions in China, the southwestern Tianshan exhibits the necessary tectonic conditions and energy accumulation for earthquakes exceeding magnitude 6. Investigating the region's tectonic deformation characteristics is crucial for understanding the background of strong earthquakes and assessing future seismic hazards.
    In this study, we compiled and analyzed multi-period GNSS velocity fields from domestic and international sources, generating a fused velocity field. The results indicate that the crustal deformation of the Tianshan seismic belt is spatially heterogeneous, influenced by the combined effects of the clockwise rotation of the Tarim Basin and the northward thrust of the Pamir Plateau. Additionally, three north-south profiles across the study area were analyzed. The western Tianshan(Profile 1)exhibits a higher north-south compression rate, which gradually decreases eastward, consistent with the velocity field trends.
    Using the fused GNSS velocity field, we calculated the strain rate field through the least-squares collocation method, obtaining strain characteristics of the Keping thrust belt and surrounding areas. The Piqiang fault zone serves as a boundary, with compressive strain on the western side of the Keping thrust belt exceeding that on the eastern side. The average surface strain rate is -3.3×10-8/a in the west and -1.2×10-8/a in the east. Similarly, the maximum shear strain follows this pattern, with values of 2.7×10-8/a in the west and 1.1×10-8/a in the east, indicating distinct deformation characteristics on either side of the Keping thrust belt.
    Focal mechanism solutions of earthquakes(M≥4.0)since 1976 were collected, and stress inversion analysis was conducted at various depths. Earthquakes in the Tianshan region predominantly occur within the upper 50km of the crust. The results reveal a nearly north-south principal compressive stress orientation, perpendicular to the Tianshan orogenic belt, aligning with the regional surface principal compressive strain rate. However, due to the presence of conjugate strike-slip faults and strike-slip earthquakes in the northern and southern Tianshan hinterland, local horizontal principal stress directions deviate from the regional trend.
    Integrating InSAR-derived deformation data with GNSS velocity field results, we calculated fault slip deficit rates and slip rates for major faults in the Keping thrust belt using a three-dimensional block model. The western Keping fault zone exhibits a higher compression rate of(3.1±0.3)mm/a, whereas the eastern section has a lower rate of(0.3±0.2)mm/a. The Maidan fault's western section has a compression rate of(2.7±0.5)mm/a, increasing to(3.7±0.4)mm/a in the east. Additionally, the left-lateral strike-slip rate is higher in the western Keping fault zone((1.5±0.3)mm/a)compared to the eastern section((0.5±0.2)mm/a). The Maidan fault follows a similar pattern, with strike-slip rates of(0.9±0.5)mm/a in the west and(2.1±0.4)mm/a in the east. The slip deficit rate distribution indicates high values in the western Keping fault zone, corresponding to the Jiashi earthquake swarm, while the eastern section exhibits lower deficit and slip rates, potentially due to a multilayer nappe structure. In the Maidan fault zone, the western section has a lower slip deficit and slip rate, suggesting weak crustal strength and limited stress accumulation, whereas the eastern section has higher values, correlating with the Wushi earthquake and subsequent energy release.
    Finally, based on the calculated strain field and the global focal mechanism earthquake catalog(1976—2021), we applied the SHIFT_GSRM2f model developed by Bird et al. to predict shallow earthquakes. The highest predicted seismic hazard values are concentrated in the southwestern Tianshan, particularly in the Pamir region, which has a higher risk of strong earthquakes than the Keping thrust belt. However, the Keping thrust belt remains one of the most seismically hazardous areas in the Tianshan region, underscoring the continued seismic risk in southwestern Tianshan.

    Table and Figures | Reference | Related Articles | Metrics
    EARTHQUAKE CASUALTY RISK ANALYSIS UNDER THE RECURRENCE SCENARIO OF THE 1902 ARTUX MS8¼ EARTHQUAKE: A CASE STUDY OF KASHGAR AND ARTUX
    NIE Wen-yu, FAN Xi-wei, LI Hua-yue, QI Yuan-meng, LIU Min
    SEISMOLOGY AND GEOLOGY    2025, 47 (2): 627-648.   DOI: 10.3969/j.issn.0253-4967.2025.02.20240156
    Abstract484)   HTML21)    PDF(pc) (16445KB)(62)       Save

    Despite global efforts to reduce earthquake disaster risk, earthquakes remain one of the most destructive natural disasters in the world. Since precise earthquake prediction remains beyond human capability, conducting accurate earthquake casualty risk assessment has become one of the most feasible and effective strategies to reduce human losses. Kashgar and Artux are located in the southwestern part of the Xinjiang Uygur Autonomous Region, on the western edge of the Tarim Basin. These are two important cities in northwest China, situated in the Pamir tectonic knot region of the India-Eurasia continental collision zone, one of the world’s most seismically active regions due to intracontinental subduction. Historically, this region has experienced multiple major earthquakes, with the most representative example being the 1902 Artux earthquake, which had a magnitude of MS8¼. This earthquake caused severe casualties and highlighted the region’s high susceptibility to seismic hazards. In this study, we used the historical seismic intensity scenario of the 1902 Artux earthquake as a deterministic earthquake scenario to understand the potential impact of a similar event under contemporary conditions. First, using time-series mobile signaling data and machine learning methods, we extracted the function types of buildings in the study area and developed a thematic dataset of building function types. Utilizing the extracted building function types, high-resolution population heatmap data, and a mapping method between grids and individual buildings, we allocated the population within the grids to specific buildings. Subsequently, by considering the relationships among building function types, temporal characteristics, local daily activity patterns, and indoor occupancy rates, we determined the spatial distribution of the indoor population in the study area. Understanding the potential damage levels of buildings under specific seismic scenarios is crucial for predicting casualty risks and formulating effective emergency response strategies. Therefore, a quantitative analysis of building damage levels under the deterministic earthquake scenario was conducted to clarify damage distribution for buildings in specific intensity scenarios. Subsequently, based on building vulnerability analysis methods, we estimated the daytime and nighttime fatality risks under the recurrence of the Artux earthquake at a 30″grid scale. The results indicate that if the 1902 MS8¼ Artux earthquake were to recur, areas with high risk of fatalities would mainly be concentrated in the densely populated urban centers of Kashgar and Artux, where buildings are densely packed, populations are concentrated, and some structures lack adequate seismic resistance. The risk of casualties is higher at night than during the daytime. To evaluate the effectiveness of population heatmap data in earthquake casualty assessment, this study compared the evaluation method based on population heatmap data with the method using the Seventh National Census data. The results demonstrate that traditional census data can only provide a rough estimate of casualties based on administrative divisions. However, its reliance on administrative boundaries means it lacks the spatial detail necessary to accurately depict casualty distribution and presents scale inconsistencies that hinder effective risk comparisons across different areas. In contrast, the assessment method based on population heatmap data enables spatial visualization of fatality risks through optimized scaling and establishes a uniform spatial comparison benchmark, thereby providing strong support for the precise allocation of emergency rescue resources. By introducing high-resolution population heatmap data and a thematic dataset of building function types, this study has enhanced the accuracy and practicality of earthquake casualty risk assessment. The results not only reveal the potential for high human casualties if a historical earthquake were to recur today but also demonstrate the distribution characteristics and temporal differences in high-risk zones, emphasizing the importance of considering temporal factors in disaster mitigation strategies. These insights provide robust scientific support for developing precise disaster mitigation and preparedness strategies in earthquake-prone urban areas.

    Table and Figures | Reference | Related Articles | Metrics
    INSAR COSEISMIC DEFORMATION AND SEISMOGENIC STRUCTURE OF THE 2024 MW7.0 WUSHI EARTHQUAKE
    CHEN Zi-long, LIU Gang, LI Qi, CHEN Wei, ZHAO Xin-yu, LIN Mu, TAO Long-wen, QIAO Xue-jun, NIE Zhao-sheng
    SEISMOLOGY AND GEOLOGY    2025, 47 (2): 429-447.   DOI: 10.3969/j.issn.0253-4967.2025.02.20240142
    Abstract477)   HTML26)    PDF(pc) (11046KB)(174)       Save

    On January 30, 2024, an MW7.0 earthquake struck the Wushi region of the southern Tianshan Mountains, Xinjiang, China. This earthquake, located in a tectonically active zone dominated by intense crustal shortening and thrust faulting, providing a valuable opportunity to investigate fault geometry and rupture mechanisms in the region. We utilized Sentinel-1 InSAR data combined with advanced inversion techniques to analyze the coseismic deformation field, determine fault parameters, and explore the spatial relationship between the mainshock and aftershocks.
    High-resolution coseismic deformation fields were generated using D-InSAR processing of Sentinel-1A ascending and descending orbit data. Nonlinear inversion methods were employed to calculate fault geometry and sliding distribution, with both single-fault and dual-fault models tested to accommodate the complex faulting characteristics. Residual analysis was performed to examine the relationship between the mainshock and the MW5.7 aftershock, and geological surveys were used to validate fault models and rupture characteristics.
    The maximum line-of-sight(LOS)displacement of the coseismic deformation field reached 70cm, displaying an elliptical pattern along the Maidan Fault. The dual-fault model revealed significant geometric complexity: the fault strike rotates clockwise by 20°~25° east of the epicenter, and the dip angle decreases from 60° in the west to 40° in the east. Fault slip was primarily concentrated west of the epicenter, characterized by high-angle thrusting with a left-lateral component, while slip in the eastern segment was lower in magnitude and relatively dispersed. The overall distribution exhibited shallow slip deficit. The correlation between geometric variations and sliding distribution suggests the presence of a geometric barrier east of the epicenter, acting as an obstacle to rupture propagation. The coseismic rupture was confined between the Yushanguxi River and the Wushi depression, with fault steps and structural complexities on the eastern and western boundaries limiting the rupture extent. The MW5.7 aftershock produced a clear LOS deformation field, with the fault strike deviating by ~10° from the mainshock trace and dipping southeast. The surface trace of the aftershock fault closely aligned with mapped surface ruptures, and the shallow slip magnitude matched the observed vertical offsets from geological surveys.
    We also demonstrates that fault geometry plays a significant role in controlling rupture propagation and termination. The geometric barrier east of the epicenter effectively limited eastward rupture propagation, while a wide fault step near the Yushanguxi River constrained the western rupture. The Wushi earthquake was identified as a blind rupture event, with no significant primary surface rupture. The distribution of secondary geological hazards aligned well with fault slip characteristics, and the spatial relationship between aftershock slip and the mainshock highlights fault segmentation within the Ush thrust belt.

    Table and Figures | Reference | Related Articles | Metrics
    GEOLOGICAL HAZARD CHARACTERISTICS AND MACRO-SCOPIC EPICENTER OF WUSHI MS7.1 EARTHQUAKE ON JANUARY 23, 2024, XINJIANG
    CHEN Jian-bo, ZHENG Li-ming, LI Tao, CHEN Jie, YAO Yuan, QIAN Li, ZHANG Bo-xuan, HU Zong-kai, LIU Guan-shen, LI Jin
    SEISMOLOGY AND GEOLOGY    2025, 47 (2): 561-576.   DOI: 10.3969/j.issn.0253-4967.2025.02.20240066
    Abstract468)   HTML22)    PDF(pc) (14230KB)(125)       Save

    The Wushi MS7.1 earthquake is the largest seismic event in Xinjiang in the past decade and the most significant earthquake in the Tianshan seismic belt since the Suusamyr MS7.3 earthquake in Kyrgyzstan in 1992. The microseismic epicenter is located at the transition between the Keping nappe structure belt and the Kuqa nappe belt in the southwestern Tianshan Mountains, near the intersection of the southern Tianshan Mountains and the Wushi Depression. This region is a compressional fault depression valley formed in the context of Tianshan orogeny. The northern boundary of the valley is controlled by the Maidan reverse fault zone, while multiple rows of newly developed linear uplifts and piedmont fault scarps in the northern depression are governed by low-angle overthrust and recoil faults. Focal mechanism analysis indicates that both the main shock and aftershocks exhibit a compressional thrust-type mechanism, consistent with the structural characteristics of the Wushishan front depression. Post-earthquake investigations revealed relatively minor casualties and structural damage, but secondary geological disasters were severe. Using field surveys, large-scale satellite image analysis, and drone aerial photography conducted immediately after the earthquake, this study investigates the spatial distribution, types, and development characteristics of seismic-induced geological hazards in various gullies within the affected region. Although no large-scale surface ruptures were observed at or near the epicenter, extensive seismic geological hazards were identified, providing crucial insights for determining the macroseismic epicenter and assessing the post-disaster intensity and loss in the earthquake area.
    Geological hazards induced by the Wushi earthquake were primarily concentrated in Biediligou, Qialematigou, Cocriu Musugou, and adjacent mountainous areas within the VII-degree seismic intensity zone. The most severe damage occurred in the Biediligou-Qialematigou region, where the macroseismic epicenter was located, with diverse types of secondary disasters, including rockfalls, landslides, rolling stones, bank slope collapses, subgrade settlement cracks, and ground fissures. Some ground fissures were accompanied by displacement. Seismic geological hazards were predominantly observed near the microseismic epicenter and within 10km of the north and south sides of the fault zone, with collapses being the most prominent, followed by pavement cracking and riverbank slope failures caused by ground fissures and uneven subgrade settlement. Post-earthquake investigations also identified numerous concentrated surface tension and shear cracks in the macroseismic epicenter, alongside earthquake-induced collapses, landslides, slope instabilities, subgrade settlements, and structural failures. Coseismic InSAR deformation field and image offset analysis revealed a maximum line-of-sight displacement of approximately 40cm and a minimum displacement of -16cm. The long axis of the deformation zone trends NE, with its strongest deformation occurring between Biediligou and Qialamatigou. This spatial distribution aligns well with the observed ground fissures and seismic deformation zones near the macroseismic epicenter, indicating that structural deformation plays a critical role in controlling the formation of earthquake-induced geological hazards.

    Table and Figures | Reference | Related Articles | Metrics
    CHARACTERISTICS OF SEISMICITY AND TECTONIC STRESS FIELD IN THE NORTHEASTERN REGION OF THE PAMIR PLATEAU
    CUI Hua-wei, YIN Xin-zhong, CHEN Jiu-hui, GUO Biao, LI Tao, YAO Yuan, LI Shi-ying, JIA Zhen
    SEISMOLOGY AND GEOLOGY    2025, 47 (2): 577-596.   DOI: 10.3969/j.issn.0253-4967.2025.02.20250047
    Abstract454)   HTML18)    PDF(pc) (15713KB)(104)       Save

    The northeastern region of the Pamir Plateau marks the convergence of the Pamir Plateau, the Tianshan Mountains, and the Tarim Basin. This area is associated with high levels of seismic activity and has experienced significant casualties due to several major earthquakes. In this study, earthquakes are relocated using the double-difference relocation algorithm, based on seismic observation reports from the China Earthquake Networks Centre(CENC)covering the period from 2008 to 2023. A total of 916 earthquake focal mechanisms were collected from various institutions and studies, including the GCMT, WWSSN, previous research, and the Xinjiang Earthquake Agency. The magnitudes of these events range from M2.8 to M8.2, and the time span extends from 1902 to 2024. The tectonic stress field was inverted using a spatial-temporal stress inversion algorithm with a 1.0°×1.0° grid, based on these historical focal mechanisms.
    The relocated earthquakes are distributed along major tectonic zones situated between the Pamir Plateau, the Tianshan Mountains, and the Tarim Basin, forming several seismic clusters within the study area. Most events have focal depths of less than 30km. Seven cross-sectional profiles were selected to analyze the spatial distribution of seismicity. The regional tectonic stress field is primarily dominated by strike-slip and reverse faulting regimes. The maximum principal stress axes are compressional in a nearly horizontal south-north(S-N)direction, while the minimum principal stress axes are extensional in the east-west(E-W)direction, with either horizontal or vertical plunges. Based on the spatial distribution of seismicity and tectonic features, the study area is divided into several subregions: the central Pamir Plateau and its eastern and western flanks; the collision zone between the Pamir Plateau and the Tianshan Mountains and its surrounding areas; the Tianshan Mountains; and the interior of the Tarim Basin.
    Preliminary results of earthquake relocation show that focal depths in the southeastern and northwestern parts of the Qimugen Arc Tectonic Belt are deeper than those within the belt. This is attributed to the subduction of the eastern and northeastern Pamir beneath the Tarim Basin toward the east and northeast. Accordingly, the Moho depth is also greater in these regions. These observations suggest that the Qimugen Arc Tectonic Belt represents a key tectonic boundary, with both the Moho and earthquake depths shallower within the belt than in its southeastern and northwestern margins. Focal depths along the southern margin of the Tianshan Mountains are deeper than those in the adjacent Tarim Basin, which is consistent with a deeper Moho in the former. Similarly, both the Moho and earthquake depths within the Ferghana fault and Piqiang fault zones are greater than those outside these faults along the southern margin of the Tianshan Mountains, indicating that these faults serve as significant tectonic boundaries. An area of sparse seismicity overlaps with the source region of the 1902 Atushi(Kashgar)MW7.7 earthquake, suggesting that accumulated crustal stress was likely fully released by this event. Seismicity in the Tarim Basin primarily consists of aftershock sequences from the 1996-1999 earthquake clusters and the 2003 Bachu-Jiashi MS6.8 earthquake. The main contributing factors to seismicity in the Tarim Basin are thought to be its low seismic wave velocity, low crustal density, and brittle crustal properties.
    The initial results of the tectonic stress inversion indicate that the central Pamir Plateau is undergoing east-west lateral extension under a strike-slip stress regime, driven by compression from the Indian Plate to the south and obstruction from the Tianshan Mountains to the north. A normal faulting regime in the central grid of the Pamir Plateau is interpreted as a center of expansion or gravitational collapse. The tectonic stress field exhibits significant variation across the flanks of the Pamir Plateau. The R-value, stress regime, and plunge of the maximum principal extension axes show gradual and, in some areas, abrupt changes from the central region to the western and eastern margins, as well as in the collision zones between the Pamir Plateau and the Tianshan Mountains. The Tianshan Mountains continue to uplift due to near S-N compression, displaying a dominant reverse faulting regime. In the Bachu uplift and depression region in the southwestern Tarim Basin, the stress field is characterized by both reverse and strike-slip regimes, likely resulting from differences in crustal velocity, density, and structural composition.

    Table and Figures | Reference | Related Articles | Metrics
    GROWING HISTORY AND GEOMORPHIC RESPONSE OF THE EASTERN TERMINATION OF KASHI ANTICLINE, SOUTHWESTERN TIAN SHAN: AN INTEGRATED ANALYSIS OF GEOLOGY, GEOMORPHOLOGY, SEISMIC REFLECTION PROFILE AND MAGNETOSTRATIGRAPHY
    HE Peng-yu, LI Tao, CHEN Zhu-xin, CHEN Jie
    SEISMOLOGY AND GEOLOGY    2025, 47 (2): 369-383.   DOI: 10.3969/j.issn.0253-4967.2025.02.20240157
    Abstract446)   HTML30)    PDF(pc) (7786KB)(201)       Save

    The Kashi anticline, located along the leading edge of the Kashi foreland thrust system in the southwestern Tian Shan of China, confronts the Atushi Anticline to the north and connects with the Mingyaole Anticline to the west. The Kashi anticline manifests at the surface as an elongated hillock with a nearly EW strike. The morphology of the anticline is roughly box-shaped, with the southern limb being gentler compared to the northern limb. As a significant component of the Tianshan orogenic belt’s frontal zone, the study of the tectonic evolution of the Kashi Anticline is crucial for understanding the Cenozoic tectonic deformation and crustal shortening processes of the Tian Shan. Previous studies on this fold are primarily focused on its surface-expression part, with little or no focus on its lateral termination that is not expressed on the land surface, which has limited the comprehensive understanding of the anticline’s overall evolutionary process. The study employs the depth-relief area method, combined with high-resolution seismic reflection profile data, to conduct a detailed structural analysis of the eastern termination of the Kashi anticline.
    Through meticulous interpretation and quantitative analysis of the seismic profiles, the following key insights have been obtained: Firstly, along the seismic profile, the Cenozoic strata are approximately 6.8km thick, while the Mesozoic-Paleozoic strata are about 2.2km thick. The eastern termination of the anticline detaches along the Paleogene unit, with a depth of ~6.8km. This detachment layer governs the deformation patterns and magnitudes of the overlying strata. Secondly, along the seismic profile, the total shortening of the anticline is estimated to be(882±79)m, of which approximately 94% is attributed to shear shortening, and about 6% is due to curvietric shortening. During the folding process, materials with an excess area of ~3.4km2 enter the cross section. Thirdly, according to the analysis of growth strata and published magnetostratigraphic data, the folding of the eastern termination initiates at the age of ~2.1Ma, which implies that the initiation age of the fold should be much older than 2.1Ma. The shortening rate remains at an approximate constant of ~0.4mm/a since the folding initiation. Fourthly, for the pre-growth strata, the uplift of the anticline gradually increases upward with depth, reaching a maximum of approximately 770m at the top boundary of the pre-growth strata. Analysis of the growth strata indicates that the uplift rate either keeps a constant of ~0.4mm/a, or increases significantly from an earlier rate of ~0.1mm/a to a later rate of 0.4 mm/a at the age point of ~1.6Ma. Notably, because the uplift rate is smaller than the sediment rate, the fold exhibits no expression on the surface of the anticline.
    Our study exemplifies that an analysis of the buried lateral termination of a fold can well determine the detachment level, shortening and uplift histories, initiation age of the fold, as well as whether or not the excess area enters the cross section during the folding process. These constraints provide a completer and more reliable basis for understanding the entire growing history of the fold. Furthermore, the results demonstrate that the analysis of buried structures in the frontal zones of orogenic belts is indispensable for a comprehensive understanding of regional tectonic deformation characteristics and evolutionary history.

    Table and Figures | Reference | Related Articles | Metrics
    COSEISMIC SLIP AND COULOMB STRESS OF THE 2024 MW7.1 HYUGANADA EARTHQUAKE
    YANG Si-yuan, HU Yan
    SEISMOLOGY AND GEOLOGY    2025, 47 (2): 597-609.   DOI: 10.3969/j.issn.0253-4967.2025.02.20240152
    Abstract442)   HTML15)    PDF(pc) (4237KB)(68)       Save

    In Southwest Japan, the Philippine Sea Plate subducts beneath the Eurasian Plate, which makes it one of the most seismically active areas in the world. This area nucleates primarily thrust-type earthquakes over the subduction interface. On August 8, 2024, a MW 7.1 earthquake occurred in the Hyuganada region in southwestern Japan, offering a valuable opportunity to investigate the seismic slip of this event. Understanding the seismic slip plays an essential role in evaluating earthquake hazards, because it helps better understand the fault behaviors and stress interactions among active crustal faults and thus potential seismic impacts.
    In this study, we first derived the coseismic displacements caused by the earthquake through Global Navigation Satellite System(GNSS)data. The maximum horizontal displacement is approximately 15cm, decreasing rapidly as the distance from the rupture zone increases. The vertical displacement shows distinct spatial characteristics. In the southeastern area of the Kyushu Island, near the rupture zone, the subsidence is predominant, with a maximum subsidence of about 6cm. Northwest of this area, the displacement gradually changes to uplift. This result indicates a complex deformation from the rupture area to the farther areas. Based on the elastic half-space model, we inverted the coseismic slip distribution using the derived GNSS displacements. In the model, we have considered a rupture fault with a length of 100km and a width of 60km. This fault is divided into 220 small rectangular elements, with 20 elements in the length direction and 11 in the depth direction. The length of each small rectangular tile is fixed at 5km, while the tile width varies with the depth. The minimum width is 5.2km, and the width increases with depth according to a proportionality coefficient of 1.01. We further use the grid search method to determine the fault dip angle and strike direction. The dip angle and strike direction of the fault are varied within ranges of -50° to 100° and 100° to 300°, respectively. Through these processes, the optimal fault dip angle and strike direction are determined to be 24° and 206°, respectively, consistent with the focal mechanism by the United States Geological Survey(USGS). The inverted coseismic slip distribution shows that the slip is primarily concentrated at 5~20km depths, featuring an elliptical pattern with a maximum slip of 1.47m. Based on the inverted fault slip model, we calculated the surface displacement using Okada’s half-space elastic dislocation theory. The results show that the model reproduces the first-order pattern of GNSS observations in both the horizontal and vertical directions. The total released moment is about 5.49×1019 Nm, corresponding to a moment magnitude of 7.09(assuming a shear modulus of 30GPa), which is close to the USGS of MW7.1.
    We calculated the Coulomb stress over the main fault using the inverted coseismic slip model. During this process, we assumed that the frictional coefficient of the fault, shear modulus, and Poisson’s ratio are to be 0.4, 30GPa, and 0.25, respectively. The fault plane’s strike, dip, and slip angles are set to be 206°, 24°, and 76°, respectively. The results reveal that outside the rupture zone, the Coulomb stress is positive, with a peak value of approximately 0.95MPa. In contrast, the overall stress is negative within the rupture zone, especially at depths between 7km and 17.5km, where it reaches a maximum value of -2.03MPa. Furthermore, the analyses on normal and shear stresses reveal a distinct pattern. In the shallow portion of the coseismic rupture zone, spanning from 0km to 7km, the normal stress is positive, with its maximum value reaching approximately 0.30MPa, whereas the shear stress is negative, peaking at around 0.73MPa. In other regions over the fault, the normal and shear stresses exhibit an increasing or decreasing trend.
    In addition, we have also investigated the aftershocks that occurred within one month after the earthquake, using data from the USGS and JMA(Japan Meteorological Agency). The location of these aftershocks features a depth range of 10~60km. To further investigate the mechanisms of the aftershock occurrence, we analyzed the distribution of maximum shear stress at a depth of 40km. The results show that aftershocks mainly occur in areas where the shear stress is more than 100kPa. The majority of aftershocks concentrate in areas with the shear stress larger than 500kPa. This suggests at least hundreds of kPa are required to induce aftershocks of this event
    To summarize, we have derived the coseismic slip and induced Coulomb stress of the 2024 Hyuganada earthquake. This reveals the seismic slip characteristics as well as impacts on the stress states and aftershock activities, which may contribute to further earthquake hazard assessment and mitigation strategies in the region.

    Table and Figures | Reference | Related Articles | Metrics
    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
    SEISMOLOGY AND GEOLOGY    2025, 47 (3): 789-805.   DOI: 10.3969/j.issn.0253-4967.2025.03.20250043
    Abstract438)   HTML13)    PDF(pc) (18102KB)(105)       Save

    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.

    Table and Figures | Reference | Related Articles | Metrics
    STUDY ON THE INFLUENCE OF LARGE LOW DIP THRUST FAULTS AND THICK BASINS ON STRONG GROUND MOTION: A STUDY IN KASHGAR AREA, XINJIANG, China
    SHEN Zhu-yue, ZHAO Shun, SUN Yao-chong, XU Du-yuan, YU Hou-yun, ZANG Nan, WANG Chun-jing
    SEISMOLOGY AND GEOLOGY    2025, 47 (2): 610-626.   DOI: 10.3969/j.issn.0253-4967.2025.02.20240158
    Abstract418)   HTML51)    PDF(pc) (9504KB)(104)       Save

    Since the Cenozoic era, the Xinjiang region has undergone significant geological evolution, resulting in numerous large-scale low-angle faults and extremely thick sedimentary basins. These unique geological features have profoundly influenced strong ground motion within the area, particularly in regions such as Kashgar, which is situated near the Pamir Plateau. This study focuses on the effects of the Tuomuluo’an Fault and the thick sedimentary layers of the Tarim Basin on seismic activity and the amplification of strong ground motion. Additionally, the Kashgar region’s geographical position south of the Tuomuluo’an front fault further highlights its vulnerability to complex seismic dynamics.
    To comprehensively understand these phenomena, researchers developed a series of low-angle thrust fault models, thick sedimentary basin models, and strong ground motion models based on detailed data from the Tuomuluo’an fault system around Kashgar and sedimentary layer characteristics in the Tarim Basin. Numerical experiments were conducted using dynamic simulations of earthquake sources and seismic wave field propagation to analyze how fault dip angle, sedimentary layer thickness and earthquake source location influence strong ground motion.
    The numerical experiments revealed several critical insights into the seismic behavior of the region. First, it was found that the Tuomuluo’an fault in the Kashgar region exhibits a greater tendency for earthquake ruptures to reach the surface when occurring along the upper ramp segment. This is due to the specific geometric configuration and stress conditions of this portion of the fault system. In contrast, ruptures along the lower slope are less likely to propagate to the surface, indicating that the fault’s dip angle significantly determines the likelihood of surface-breaking earthquakes. Second, the study demonstrated that the thick sedimentary basin of the Tarim Basin significantly amplifies strong ground motion in the Kashgar region. This amplification occurs because the soft, low-velocity sediments in the basin trap seismic waves and cause prolonged shaking, increasing ground motion intensity during an earthquake. To further investigate this phenomenon, a scenario earthquake was simulated with its nucleation location on the eastern segment upper ramp of the Tuomuluo’an fault. The results generated peak ground velocity and intensity maps essential for disaster assessment and planning.
    The simulation revealed that earthquake disasters in the Kashgar area are divided into two distinct regions: one near the seismogenic fault extending along its direction, where strong shaking is expected to be most intense, and another above the Tarim thick sedimentary basin, where ground motion is significantly amplified due to the geological characteristics of basin. Additionally, complex terrain changes influenced earthquake disaster zones, highlighting the importance of topographical factors in determining seismic hazard distribution.
    The study concludes that several key factors significantly affect strong ground motion in the Kashgar region: the dip angle and geometric shape of the Tuomuluo’an fault, the thickness of sedimentary layers, and the regional topography. These findings enhance our understanding of seismic hazard assessment in the area and provide valuable insights for disaster preparedness and mitigation strategies. The geological conditions of the Xinjiang region, characterized by large-scale low-angle faults and thick sedimentary basins since the Cenozoic era, further underscore the importance of these findings for regional earthquake risk evaluation.
    This research provides critical scientific evidence to improve seismic hazard assessment in the Kashgar region. The study offers practical recommendations for enhancing earthquake preparedness and reducing potential disaster impacts in the Xinjiang region by identifying the key factors that influence strong ground motion.

    Table and Figures | Reference | Related Articles | Metrics
    HIGH-RESOLUTION SHALLOW CRUSTAL S-WAVE VELOCITY STRUCTURE IMAGING IN THE KASHGAR, XINJIANG
    HUA Qian, PEI Shun-ping, LI Tao, LIU Han-lin, LIU Wei, LI Lei, LI Jia-wei, YANG Yi-hai
    SEISMOLOGY AND GEOLOGY    2025, 47 (2): 533-546.   DOI: 10.3969/j.issn.0253-4967.2025.02.20240160
    Abstract418)   HTML24)    PDF(pc) (7312KB)(132)       Save

    High-density short-period seismometers are increasingly employed in urban environments and local geological structures to explore the crustal structure, their high-resolution images facilitate the precise identification of subsurface faults, the spatial distribution of mineral resources, and conduct building response analysis. In this study, 101 short-period seismometers were deployed across Kashgar for continuous seismic monitoring. Integrating with ambient noise tomography, high-resolution seismic velocity imaging of the shallow crust within the Kashgar was conducted. This study aims to delineate potential subsurface faults, elucidate their tectonic genesis, and provide critical insights for regional seismic risk assessment.
    The empirical Green’s functions extracted from the cross-correlation of the Z component yielded a total of 1 752 Rayleigh wave phase velocity and group velocity dispersion curves. By applying one-step ambient noise tomography, the three-dimensional S-wave velocity structure of the study area was resolved down to a depth of 5km, achieving a lateral resolution of approximately 0.04°. The horizontal and vertical cross-sections of the S-wave velocity model reveal that the S-wave velocities within the upper 5km of the crust in the study area are generally lower than the global average velocity model. The Kashgar Depression is characterized predominantly by Cenozoic sediments, with continuous Quaternary alluvial deposits reaching thicknesses of up to 10~12km. The relatively weak Cenozoic sedimentary basin likely contributes to the overall low S-wave velocities observed in the region.
    The velocity structures exhibit remarkably consistent patterns in varying depths. Below the depth of 1.2km, three notable low-velocity anomalies(LVAs), labeled L1, L2, and L3, are identified beneath the Kashgar. Among these, L1 and L2 form an approximately 16km long, east-west trending bowtie-shaped LVAs that align with the structural trend of the Kashgar anticline. These anomalies cover much of central Kashgar and extend nearly vertically to depths shallower than 5km, showing variation in shape and size at different depths. L3, located at the central southern edge of Kashgar, appears as a semicircle with a diameter of about 10km. Its extent diminishes gradually with increasing depth, which may indicate lithological variations at different depths.
    Based on the integration of seismic reflection profiles, we infer that the frontal zone of the Southwestern Tianshan fold-and-thrust belt has developed multiple north-dipping thrust structures and south-dipping secondary thrust faults propagating toward the basin. These deformations penetrate the entire sedimentary cover, forming multi-level detachments at varying depths. Notably, the slippage of the mud layer(approximately 4km deep)at the base of the sedimentary cover in the Kashgar Depression represents the shallowest detachment layer identified in previous studies. This suggests that the multi-layered weak slippage zones within the sedimentary sequence of the Kashgar Depression may be responsible for the formation of the bowtie-shaped LVAs. Mechanically weak detachment layers likely play a key role in shaping these anomalies.
    Furthermore, the Kashgar-Atushi fold-and-thrust system has experienced both lateral propagation and along-strike shortening during ongoing tectonic activity, resulting in the progressive advancement of the fold-and-thrust system towards the Kashgar Depression, which lies adjacent to the collapse-reverse fault system, may also have been subjected to intense tectonic action to form similar faults. Consequently, the east-west trending bowtie-shaped LVAs may indicate the presence of a secondary blind fault parallel to the Kashgar anticline. This inferred fault crosses the tectonic boundary between the southwestern Tianshan and Pamir regions, exhibiting a significant east-west structural discontinuity. Geological and geomorphological evidence reveals that the Kizilsu River, the largest river in the region, and its tributaries intersect the LVAs beneath Kashgar. We hypothesize that these LVAs may also reflect high-porosity fluvial sediments and folded scarps associated with paleo-river deformation.
    In summary, high-density short-period seismic array imaging enables the precise detection of shallow subsurface structures in urban environments. This approach provides robust datasets for urban active fault detection, seismic amplification effect evaluation, and subsurface resources and energy exploration and development.

    Table and Figures | Reference | Related Articles | Metrics
    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
    SEISMOLOGY AND GEOLOGY    2025, 47 (3): 869-880.   DOI: 10.3969/j.issn.0253-4967.2025.03.20250044
    Abstract415)   HTML7)    PDF(pc) (7125KB)(60)       Save

    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.

    Table and Figures | Reference | Related Articles | Metrics
    DENSE BENDING MOMENT NORMAL FAULT SCARPS ALONG THE GUMAN ANTICLINE AT THE FOOTHILL OF THE WEST KUNLUN MOUNTAINS
    XU Jian-hong, CHEN Jie, LI Tao, ZHANG Bo-xuan, DI Ning
    SEISMOLOGY AND GEOLOGY    2025, 47 (2): 405-428.   DOI: 10.3969/j.issn.0253-4967.2025.02.20240148
    Abstract410)   HTML16)    PDF(pc) (15306KB)(98)       Save

    Bending-moment fault and flexural-slip fault are two types of fold-related faults in compressional tectonic environments. Historical earthquake records suggest that both fault types may be active simultaneously, with their fault scarps providing crucial insights into strong seismic events. In the northern region of the Guman anticline, located at the foothills of the West Kunlun Mountains, numerous prominent bending-moment normal fault scarps have developed, reaching heights between 0.5m and 16.0m. This study focuses on a fault scarp segment approximately 5.4km long and 4.2km wide. A digital elevation model(DEM)with a 0.2m resolution was generated using drone photogrammetry. A total of 739 cross-fault scarp profiles were extracted, providing key parameters such as scarp height, slope, displacement continuity, and cumulative displacement trends. Data analysis yielded the following findings:
    (1)In the study area, dense bending-moment normal faults align along the active anticline axis, dipping toward the axial plane at angle of 70°~80°, as observed in a trench. Among these faults, more than a dozen dip northward, whereas only 1-2 dip southward, forming asymmetric grabens. This asymmetry may be attributed to the overall northward tilt of the strata and the differing limb structures of the underlying anticline. These faults divide the terrace surfaces into multiple rectangular blocks, 380~650m wide. The blocks exhibit outward tilting relative to the fold axis, with those cut by north-dipping faults tilting southward and those cut by south-dipping faults tilting northward. The degree of tilting and fault displacement is closely related to the thickness of the underlying anticlinal strata and the extent of stratal bending.
    (2)Displacement profiles along the faults reveal a step-like decrease in displacement as terrace surfaces become progressively younger, with maximum slope profiles displaying similar trends. This pattern suggests long-term fault activity. Cumulative displacement data confirm this trend, with displacement values of(54.5±3.3)m for terrace T3c and(19.5±1.1)m for terrace T1c. The total displacement of T3c is approximately 2.8 times that of T1c, and displacement ratios across different terraces range from 1.5 to 5.5. Higher ratios indicate greater displacement accumulation on older terraces, suggesting an earlier onset of fault activity. These displacement rankings imply that an initial framework of faults developed in the region, followed by subsequent fault intrusion. Notably, Fault F8 exhibits a displacement ratio of 5.5, forming a(1.0±0.3)m high fault scarp on the young T1b terrace, indicating that even the earliest-formed faults remain active.
    (3)Seismic reflection profiles reveal that the south flank of the Guman anticline dips 3°~6° northward, while the north flank dips 12°~14° northward. The underlying blind thrust exhibits a lower flat-ramp-upper flat geometry. However, bending-moment normal faults are not visible in the seismic reflection data, suggesting that they are secondary structures associated with anticline deformation. The fault zone aligns with the anticline’s fault-bend axis, indicating ongoing activity in the anticline zone. The bending-moment normal faults are rootless, meaning they are not primary seismogenic faults. Instead, they primarily develop in poorly layered strata and are largely independent of the kinematics of fold growth. Their formation is closely tied to the degree of strata bending and the thickness of overlying beds. Despite their shallow nature, the bending-moment normal faults exhibit long-term activity, providing evidence that the underlying anticline remains active. These findings support the interpretation of the Guman anticline as an active fault-bend fold.

    Table and Figures | Reference | Related Articles | Metrics
    CHARACTERISTICS OF HOLOCENE CATASTROPIC EVENT LAYERS AND THEIR RELATIONSHIP WITH REGIONAL EARTHQUAKES IN THE SOUTHERN QINGSHUIHE BASIN, NINGXIA
    HUANG Ting, WU Fang, XIA Cai-xiang, LI Zhen-hong, DONG Xiao-peng, WU Zhong-hai, KOU Lin-lin
    SEISMOLOGY AND GEOLOGY    2025, 47 (4): 1036-1057.   DOI: 10.3969/j.issn.0253-4967.2025.04.20240110
    Abstract407)   HTML14)    PDF(pc) (15724KB)(50)       Save

    Reconstructing paleoseismic sequences using direct geological evidence from fault zones remains challenging in regions dominated by basement rocks, folded terrains, or loess-covered surfaces. In contrast, far-field sedimentary records—particularly those from fluvial and lacustrine facies—offer key advantages due to their continuity, sensitivity, and high temporal resolution, partially compensating for limitations in instrumental records, historical documents, and trench-based fault studies. As such, these sedimentary archives are essential for advancing the understanding of long-term fault activity and regional seismic hazards.

    This study focuses on the Qingshuihe Basin, located at the northeastern margin of the Tibetan plateau, where the Haiyuan-Liupanshan fault zone intersects the Xiangshan-Tianjingshan fault zone—a region of intense tectonic activity and frequent strong earthquakes. Despite prominent geomorphic expressions, thick loess cover and the absence of continuous marker beds hinder fault activity studies at trench or outcrop scales. To address this, we conducted a detailed survey of Holocene fluvial-lacustrine catastrophic event deposits in the southern Qingshuihe Basin, integrating stratigraphic interpretation, chronological constraints, and regional historical-seismological correlations to explore their links to regional seismic activity.

    Three distinct catastrophic event layers were identified within the fluvial-lacustrine stratigraphy, exhibiting characteristics of sudden deposition, local sediment sourcing, and chaotic accumulation. These are often accompanied by faulting, ruptures, and the development of cultural layers. Small-scale faults at the base of these layers show consistent orientation patterns. Their stratigraphic positions correlate with the formation horizons of purplish-red clay veins in adjacent loess deposits—features widely interpreted as surface expressions of tectonic deformation. Both sets of anomalous deposits align with the strike of basin-margin faults, reinforcing their tectonic origin.

    To constrain the timing of these events, we employed AMS radiocarbon dating and optically stimulated luminescence(OSL)techniques. AMS 14C ages were calibrated using OxCal v4.4.4 and the INTCAL20 calibration curve, yielding a 95.4%confidence interval(2σ). Results indicate three major events since the mid-Holocene: i.e. E1: (6 220±95)to(5 393±49)cal aBP; E2: (3 411±30) to (797±52)cal aBP (approximated near(797±52)cal aBP); E3: (797±52) to (730±26)cal aBP.

    Comparative analysis with regional earthquake records shows a strong correlation as follows:

    E1 likely corresponds to a major earthquake(MW≥7.0)affecting the eastern Haiyuan and Liupanshan faults between(6 600±500)and(5 640±540)cal aBP, producing minor faulting in the study area.

    E2 aligns with the 1219AD Guyuan South Earthquake(M6¾), related to reverse-thrust faulting on the Liupanshan and Guanshan faults, producing widespread collapse deposits and an estimated intensity of Ⅶ in the basin.

    E3 corresponds to the 1306AD Kaicheng Lu Earthquake(MS7.0), which also triggered collapse deposits, though with slightly lower intensity(Ⅵ-Ⅶ)due to differences in magnitude, epicentral location, and distance.

    This study underscores the value of far-field sedimentary archives in reconstructing seismic histories in tectonically complex regions like the northeastern Tibetan plateau. The results provide new paleoseismic constraints and contribute valuable data for regional seismic hazard assessments.

    Table and Figures | Reference | Related Articles | Metrics
    PN WAVE VELOCITY AND ANISOTROPY TOMOGRAPHY IN THE UPPERMOST MANTLE OF PAMIR PLATEAU AND ADJACENT REGIONS
    LIU Jia-xin, PEI Shun-ping, GUO Yi-cun
    SEISMOLOGY AND GEOLOGY    2025, 47 (2): 547-560.   DOI: 10.3969/j.issn.0253-4967.2025.02.20240155
    Abstract373)   HTML7)    PDF(pc) (13201KB)(70)       Save

    Since the Cenozoic era, the ongoing collision between the Indian Plate and the Eurasian Plate has formed the largest and youngest continent-continent orogenic belt on Earth. The Pamir Plateau, located at the western end of the India-Eurasia collision zone, is one of the most tectonically active and structurally complex regions globally. It is characterized by widespread folds and faults, frequent M≥7.0 earthquakes, and numerous intermediate to deep-focus earthquakes, making it an ideal natural laboratory for studying plate tectonics and orogenic processes. Investigating the velocity and anisotropic structure of the uppermost mantle beneath the Pamir Plateau and its surrounding regions is of great significance for understanding the tectonic deformation characteristics and dynamic mechanisms of continent-continent collisions. The Pn-wave travel-time tomography method is an effective approach for studying the physical properties of the uppermost mantle, it offers several advantages: 1)Pn-wave ray paths are concentrated in the uppermost mantle, and there is a relatively abundant record of arrival times; 2)this method can simultaneously obtain mutually constrained velocity and anisotropy structures of the uppermost mantle; and 3)although it only provides lateral velocity and azimuthal anisotropy information for the uppermost mantle, it offers higher resolution and accuracy, and the results are not influenced by deeper mantle structures. In this study, using 363 414 Pn-wave arrival times recorded from 103 190 events at 471 stations, we performed Pn-wave travel-time tomography to obtain higher-resolution seismic velocity and anisotropy images of the uppermost mantle beneath the Pamir Plateau and adjacent regions compared to previous studies. The tomography results exhibit a clear correlation with geological structures, with the following key findings: 1)Lower Pn-wave velocities are observed in tectonically active regions such as the Pamir Plateau, Tienshan, Hindu Kush, West Kunlun, and Altyn Tagh, while higher velocities are found in stable ancient blocks like the Indian Plate, Tarim Basin, Tajik Basin, Junggar Basin, and Fergana Basin. 2)The India-Eurasia collision has generated strong Pn-wave azimuthal anisotropy. On both sides of the Pamir collision zone, in the Indian Plate, Tarim Basin, and Tajik Basin, the anisotropy axes align with the plate motion directions. In contrast, in the central part of the collision zone, the anisotropy axes are nearly perpendicular to the maximum compressive stress and crustal motion directions. A similar anisotropy pattern is observed on both sides of the Tienshan. Pn-wave velocities primarily reflect the properties of the uppermost mantle, which is predominantly composed of peridotite. Temperature has a more significant influence on velocity than pressure, and temperature variations are closely linked to tectonic activity. Typically, stable cratonic regions exhibit higher Pn-wave velocities, while tectonically active or volcanic regions with significant fluid activity show lower velocities. The velocity differences between active and stable blocks are mainly attributed to temperature variations. Comparing the surface boundary between the Eurasian and Indian Plates with the high-velocity anomalies observed in tomography reveals that the high-velocity anomaly of the Indian Plate extends approximately 200km northward beneath the Tibetan plateau. Pn-wave azimuthal anisotropy is generally attributed to the preferred orientation of olivine crystals caused by mantle deformation. The anisotropy direction of the Indian Plate is predominantly north-south, consistent with its GPS motion and maximum compressive stress directions. From the Tajik Basin to the western Tarim Basin, the fast-axis anisotropy direction gradually shifts from NW to NE, aligning well with the maximum compressive stress and GPS directions. This is primarily due to the simple shear between the crust and upper mantle caused by crustal shortening and uplift, while the lithospheric mantle subducts during intense plate collision. Similarly, in the Fergana Basin, northern Tarim Basin, and Junggar Basin, the anisotropy directions are nearly north-south, consistent with GPS and maximum compressive stress directions. In regions of intense deformation, such as the Hindu Kush, Pamir Plateau, West Kunlun, Altyn Tagh, and Tienshan orogenic belts, the anisotropy directions are perpendicular to the maximum compressive stress directions, indicating strong pure shear deformation in the uppermost mantle. This suggests that in continent-continent collision zones, not only does the crust undergo significant shortening and uplift, but the uppermost mantle also experiences substantial compressional deformation. Finally, we propose an improved dynamic model of continent-continent collision to elucidate the collision process between the Indian Plate and the Tarim and Tajik Basins, as well as the mechanisms of anisotropy formation.

    Table and Figures | Reference | Related Articles | Metrics
    STRUCTURAL SETTING AND IMPLICATIONS OF THE 2015 PISHAN MW6.4 EARTHQUAKE IN THE FRONT OF WESTERN KUNLUN
    YANG Wen-xin, LI Tao, CHEN Jie, YAO Yuan
    SEISMOLOGY AND GEOLOGY    2025, 47 (2): 507-532.   DOI: 10.3969/j.issn.0253-4967.2025.02.20250049
    Abstract356)   HTML14)    PDF(pc) (7331KB)(117)       Save

    As the dynamic hub between collisional orogenic belts and foreland basins, the formation and propagation of foreland fold-and-thrust belts are often accompanied by destructive earthquakes. In this process, the detachment, which is regarded as an incompetent layer, plays a significant controlling role in the propagation of the thrust system and the occurrence of strong earthquakes. Therefore, research focused on the structural patterns and physical property architecture of foreland fold-and-thrust belts is conducive to analyzing their activity characteristics and seismogenic patterns, thereby providing a basis for regional seismic hazard assessments.
    The western Kunlun orogenic belt is situated at the Qinghai-Xizang Plateau’s northwestern margin and the Tarim Basin’s southwestern edge. Governed by the convergence of the Indian and Eurasian plates, the compressive collision between the Qinghai-Xizang plateau and the Tarim Basin has led to the formation and propagation of the Western Kunlun foreland fold-and-thrust system, triggering intense and frequent tectonic and seismic activities within the region. A MW6.4 earthquake occurred in Pishan on July 3, 2015 in the front of western Kunlun, which was one of the most destructive seismic events that happened in Xinjiang in the recent decade. The epicentral area of this earthquake experienced an intensity of VIII, with many regions feeling strong tremors. The disaster area was vast, affecting a large population, and was accompanied by hazards such as sand liquefaction, ground fissures, and collapses. Studies on this event will provide an opportunity to better understand active structures, stress state, and future seismic hazard of the front of western Kunlun and the adjacent region.
    Synthesizing geology, geomorphology and seismology studies, Slik Anticline, located on the western Kunlun foreland fold-and-thrust belt, is regarded as the seismogenic fault of this event. The Slik Anticline is located within the Guman fold belt of the Hotan thrust belt, which lies in the eastern segment of the western Kunlun foreland. The Slik ramp, situated beneath the Slik Anticline, is approximately at a depth of 10 to 15km, with a ramp angle of about 20 degrees. The Pishan earthquake ruptured blind ramp under the Slik Anticline. From north and shallow to south and deep level, the fault is revealed with the structural style of lower flat, ramp and upper flat and may merge with other thrust in a deeper part near to the hinterland. The upper and lower flats are corresponded to regional detachment of lowest Paleogene and inner Cambrian incompetent strata, respectively. And the mainshock occurred at the Slik ramp which connecting these two detachment, which implicates that the Pishan earthquake is an event caused by independent structural unit. The postseismic deformation mainly occurred on the lower flat as afterslip, and accumulated gradually over time. The moment released after a slip in 5 months corresponds to a sub-MW6.3 event, and is equivalent to the moment released by the mainshock.
    Worldwide, earthquakes exhibiting similar characteristics include the 2015 Gorkha MW7.8 earthquake in Nepal, the 2005 Kashmir MW7.6 earthquake, the 1999 Chi-Chi MW7.6 earthquake in Taiwan, China, and the 1995 Colima-Jalisco MW8.0 earthquake in Mexico, among others. Conducting detailed analyses of these seismic events provides a valuable entry point for clarifying the structural patterns, stress states, and seismogenic features of foreland fold-and-thrust fault zones. Furthermore, it aids in the in-depth analysis of regional seismogenic models and predicting future seismic trends.
    Focusing on the Pishan earthquake, expect settling slip by the folding strata above the upper flat; this phenomenon may also depend on both matter properties and gravitational differences between the upper and lower flats. 1)During the fault deformation, a detachment with low competence and low friction coefficient is unable to accumulate significant stress, resulting in the absence of notable stress drops during frictional sliding and instead manifests as a state of after-slip activity. 2)When the detachment, especially the shallow one, has an undersized gravitational load, it may fail to effectively play a role in detachment, leading to a geological process where fault activity may be altered or inhibited. Occurring at the lower flat, the afterslip of the Pishan earthquake cooperated these two situations. 1)The deeper detachment, located on the Cambrian gypsum, may have a lower friction coefficient compared with the upper detachment. 2)And apparently, the deeper detachment possesses more gravitational load. 3)Plus, stress cannot accumulate to a high level at a flat part of a seismogenic fault. These points may lead to the afterslip release of accumulated stress, mostly on the lower flat, as the postseismic deformation. Therefore, the risk of a major earthquake occurring at the lower flat of the fault may be relatively low.
    On the one hand, this phenomenon may cause stress loading on the lower fault ramp or the root fault at a deeper level below the fault flat, increasing the seismic risk in the deep structures. On the other hand, if an earthquake was trigged on the lower ramp, the lower flat above it, which had undergone a slip and released a significant amount of energy, would limit the extent of the rupture.
    Furthermore, Cambrian gypsum layers are widely present at depth, which can serve as deep detachment surfaces and lower fault flats at the scale of the study area and its adjacent region. Therefore, whether the seismic activity and energy release pattern of the Slik anticline can represent the seismogenic patterns in the front of western Kunlun, and whether such a pattern can effectively provide a basis for regional seismic hazard assessment, deserves further attention and research.

    Table and Figures | Reference | Related Articles | Metrics
    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
    SEISMOLOGY AND GEOLOGY    2025, 47 (3): 689-706.   DOI: 10.3969/j.issn.0253-4967.2025.03.20250034
    Abstract353)   HTML23)    PDF(pc) (10876KB)(209)       Save

    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.

    Table and Figures | Reference | Related Articles | Metrics
    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
    SEISMOLOGY AND GEOLOGY    2025, 47 (3): 734-746.   DOI: 10.3969/j.issn.0253-4967.2025.03.20250033
    Abstract322)   HTML21)    PDF(pc) (8453KB)(99)       Save

    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.

    Table and Figures | Reference | Related Articles | Metrics
    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
    SEISMOLOGY AND GEOLOGY    2025, 47 (3): 932-948.   DOI: 10.3969/j.issn.0253-4967.2025.03.20250019
    Abstract317)   HTML5)    PDF(pc) (10837KB)(58)       Save

    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.

    Table and Figures | Reference | Related Articles | Metrics
    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
    SEISMOLOGY AND GEOLOGY    2025, 47 (3): 949-968.   DOI: 10.3969/j.issn.0253-4967.2025.03.20250026
    Abstract298)   HTML12)    PDF(pc) (11410KB)(45)       Save

    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.

    Table and Figures | Reference | Related Articles | Metrics
    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
    SEISMOLOGY AND GEOLOGY    2025, 47 (3): 777-788.   DOI: 10.3969/j.issn.0253-4967.2025.03.20250037
    Abstract298)   HTML19)    PDF(pc) (5998KB)(65)       Save

    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.

    Table and Figures | Reference | Related Articles | Metrics
    3-D CRUSTAL S-WAVE VELOCITY STRUCTURE IN AND AROUND THE DATONG VOLCANIC GROUP: CONSTRAINTS FROM DIRECT TOMOGRAPHIC IMAGING OF AMBIENT-NOISE SURFACE WAVES
    LI Ruo-hao, LEI Jian-she, SONG Xiao-yan
    SEISMOLOGY AND GEOLOGY    2025, 47 (4): 1090-1112.   DOI: 10.3969/j.issn.0253-4967.2025.04.20240121
    Abstract296)   HTML28)    PDF(pc) (13539KB)(74)       Save

    The Datong volcanic group is located in the central part of the North China Craton, and it has attracted widespread attention due to its complex geological tectonic environment with a high level of seismicity. In this study, we collect continuous seismic waveforms from January to December 2020 recorded at 56 provincial permanent seismic stations of the China Earthquake Administration. In data processing, we first conduct rigorous preprocessing of the raw waveform data, including mean removal, detrending, and bandpass filtering to ensure data quality. After performing cross-correlation of ambient noise, we manually extract Rayleigh wave dispersions in the periods of 5-30s, which effectively reflect the velocity structural characteristics of the crustal and uppermost mantle. Based on the extracted dispersion data, we apply the direct surface wave tomographic method to construct a three-dimensional S-wave velocity structure model extending to a depth of 40km with a spatial resolution of 0.75°×0.75° in the horizontal directions.

    Our imaging results show that the distribution of S-wave velocities corresponds well with geological structural features in the upper crust. The Shanxi rift zone generally exhibits low-velocity anomalies, reflecting the structural characteristics of the Taiyuan Basin, Xinding Basin, and Datong Basin, which are speculated to be related to the Cenozoic sedimentary layers covering the shallow subsurface in this area. In contrast, the Lüliang Mountains and Taihang Mountains exhibit high-velocity anomalies due to their exposed bedrock. In the middle to lower crust, the low-velocity anomaly beneath the Datong volcanic group extends across the northern part of the Shanxi rift zone to the west of the zone, possibly caused by extensive magma activity in the crust due to the upwelling of hot mantle materials under the region. A discontinuous low-velocity anomaly body exists in the crust beneath the Datong volcanic area, potentially serving as a conduit for magma upwelling, but with a possibly discontinuous magma supply. Combining previous deep mantle imaging studies, we speculate that the crustal low-velocity anomalies reflecting hot materials beneath the Datong volcanoes could be jointly caused by the westward deep subduction of the Pacific slab, the extrusion of asthenospheric mantle materials by the Indo-Eurasian collision, and mantle plume activities. Our findings not only deepen our understanding of the deep structure and dynamics of Datong volcano, but also provide new insights into understanding the tectonic evolution of the North China Craton.

    Table and Figures | Reference | Related Articles | Metrics
    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
    SEISMOLOGY AND GEOLOGY    2025, 47 (3): 718-733.   DOI: 10.3969/j.issn.0253-4967.2025.03.20250012
    Abstract295)   HTML19)    PDF(pc) (16018KB)(79)       Save

    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.

    Table and Figures | Reference | Related Articles | Metrics
    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
    SEISMOLOGY AND GEOLOGY    2025, 47 (3): 835-849.   DOI: 10.3969/j.issn.0253-4967.2025.03.20250028
    Abstract294)   HTML3)    PDF(pc) (5197KB)(18)       Save

    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.

    Table and Figures | Reference | Related Articles | Metrics
    SURFACE RUPTURE CHARACTERISTICS OF THE JISHISHAN MS6.2 EARTHQUAKE ON DECEMBER 18, 2023
    LI Lin-lin, JIANG Wen-liang, LI De-wen, JIAO Qi-song, LUO Yi, LI Yong-sheng, TIAN Yun-feng, LI Ying-ying
    SEISMOLOGY AND GEOLOGY    2025, 47 (4): 1058-1074.   DOI: 10.3969/j.issn.0253-4967.2025.04.20240010
    Abstract291)   HTML22)    PDF(pc) (13335KB)(118)       Save

    At 23:59 on December 18, 2023, an MS6.2 earthquake struck Jishishan County, Linxia Hui Autonomous Prefecture, Gansu Province. The epicenter(35.70°N, 102.79°E) was located in the southeastern segment of the Lajishan fault zone, with a focal depth of approximately 10km. According to the Ministry of Emergency Management, the maximum intensity reached Ⅷ degree, with the NNW-striking long axis of the isoseismal zone consistent with the fault strike. Although moderate in magnitude, the earthquake has caused over 150 fatalities, primarily due to its occurrence at midnight, high population density, poor seismic resistance of housing, and secondary hazards such as debris flows.

    Following the event, a comprehensive field investigation was conducted at the epicentral area. Utilizing UAV imagery, digital surface models(DSMs) derived from GF-7 satellite data, and InSAR analysis using Sentinel-1 SAR data, the characteristics of the surface ruptures and the seismogenic structure were examined.

    The InSAR results from the ascending orbit of Sentinel-1 revealed that coseismic deformation was dominated by uplift, with a maximum line-of-sight(LOS) displacement of approximately 65mm. The primary deformation zone exhibited an elliptical shape, trending NNW, and extended approximately 10km along the fault strike—consistent with the expected rupture length for an event of this magnitude.

    Field surveys and UAV imagery identified a ~1km-long NNW-trending surface rupture east of Yinjiashan Village, which cut across multiple geomorphic units. These ruptures exhibited dominantly thrusting motion with minor right-lateral strike-slip components and were linearly distributed, indicating a tectonic origin rather than landslides or secondary processes. The maximum observed offset reached 8cm vertically and 2cm horizontally.

    Digital surface model interpretation from GF-7 imagery revealed several NNW-trending linear structures along the eastern front of the Jishishan Mountains, forming linear topographic scarps and ridges. The observed surface rupture corresponds with one of these structures(F12). Additionally, two local rivers exhibit sharp deflections toward the NNW, controlled by these structures, supporting the interpretation that they are branch faults of the northern Lajishan fault.

    According to empirical relationships between magnitude, rupture length, and displacement(Wells & Coppersmith, 1994), an MS6.2 event is expected to produce a rupture length of approximately 10km. Similar surface rupture lengths(~15km and ~11km) were observed in the 2021 MS6.1 Biru earthquake. Combined with InSAR-derived deformation extent and the lack of field coverage south of the observed rupture, it is inferred that the surface rupture may extend several kilometers southward along the same structural lineament.

    In conclusion, the seismogenic fault responsible for this event is the northern Lajishan fault, comprising multiple NNW-trending branch faults along the eastern front of the Jishishan Mountains. The surface rupture identified corresponds to one of these secondary faults. Despite of the modest scale of the Lajishan fault zone compared to major structures like the Altyn Tagh, East Kunlun, and Qilian-Haiyuan fault zones, and its lack of historical large earthquakes, this event highlights the potential seismic hazard posed by smaller faults under the influence of ongoing crustal uplift and tectonic extension in the northeastern Tibetan plateau.

    Table and Figures | Reference | Related Articles | Metrics
    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
    SEISMOLOGY AND GEOLOGY    2025, 47 (3): 917-931.   DOI: 10.3969/j.issn.0253-4967.2025.03.20250031
    Abstract288)   HTML8)    PDF(pc) (8399KB)(27)       Save

    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.

    Table and Figures | Reference | Related Articles | Metrics
    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
    SEISMOLOGY AND GEOLOGY    2025, 47 (3): 806-819.   DOI: 10.3969/j.issn.0253-4967.2025.03.20250029
    Abstract285)   HTML20)    PDF(pc) (6193KB)(62)       Save

    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.

    Table and Figures | Reference | Related Articles | Metrics
    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
    SEISMOLOGY AND GEOLOGY    2025, 47 (3): 747-760.   DOI: 10.3969/j.issn.0253-4967.2025.03.20250021
    Abstract276)   HTML14)    PDF(pc) (7688KB)(57)       Save

    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.

    Table and Figures | Reference | Related Articles | Metrics
    RELOCATION AND FORESHOCK SEQUENCE IDENTIFICATION OF DINGRI MS6.8 EARTHQUAKE IN XIZANG
    YIN Xin-xin, ZUO Ke-zhen, ZHAO Cui-ping, CAI Run
    SEISMOLOGY AND GEOLOGY    2025, 47 (3): 850-868.   DOI: 10.3969/j.issn.0253-4967.2025.03.20250025
    Abstract271)   HTML10)    PDF(pc) (7453KB)(73)       Save

    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.

    Table and Figures | Reference | Related Articles | Metrics
    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
    SEISMOLOGY AND GEOLOGY    2025, 47 (3): 984-998.   DOI: 10.3969/j.issn.0253-4967.2025.03.20250010
    Abstract265)   HTML5)    PDF(pc) (9350KB)(50)       Save

    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.

    Table and Figures | Reference | Related Articles | Metrics
    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
    SEISMOLOGY AND GEOLOGY    2025, 47 (3): 881-896.   DOI: 10.3969/j.issn.0253-4967.2025.03.20250042
    Abstract257)   HTML7)    PDF(pc) (5456KB)(55)       Save

    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.

    Table and Figures | Reference | Related Articles | Metrics
    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
    SEISMOLOGY AND GEOLOGY    2025, 47 (3): 820-834.   DOI: 10.3969/j.issn.0253-4967.2025.03.20250040
    Abstract255)   HTML9)    PDF(pc) (6328KB)(44)       Save

    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.

    Table and Figures | Reference | Related Articles | Metrics
    SIMULTANEOUS INVERSION OF CRUSTAL VELOCITY STRUCTURE AND EARTHQUAKE RELOCATION IN THE NORTHWEST OF THE BEIJING AREA
    GONG Meng, ZOU Xian-kun, WANG Xiao-shan, LI Guang, SHENG Shu-zhong, LI Hong-xing, XU Rong-hua, LU Chang-sheng
    SEISMOLOGY AND GEOLOGY    2025, 47 (4): 1132-1151.   DOI: 10.3969/j.issn.0253-4967.2025.04.20240014
    Abstract248)   HTML19)    PDF(pc) (8398KB)(71)       Save

    The northwest Beijing area is located in the northwest of the ancient North China Craton block. Due to the long-term and frequent geological structure evolution and tectonic movement of the North China block, the complex geological structure pattern has been created in this area, with the Yanshan tectonic belt in the north, the North China rift basin in the south, the Shanxi depression belt in the west, and the Bohai Sea in the east. Due to its unique geographical location and frequent seismic activity, this area has long been a key concern for geologists and seismologists.

    We collected P waves’ absolute and relative travel time of 20 442 earthquakes in the northwest of the Beijing area record by 145 seismic stations deployed in Hebei, Shanxi and Neimeng Provinces, during January 2009 to December 2020 and used double-difference seismic tomography joint inverted the Seismic source location parameters and the 3D P-wave velocity structure of the study area. To improve the uniformity of ray coverage and the accuracy of data in inversion, seismic phases were selected based on the following conditions. 1)The P-wave phases of each earthquake are required to be recorded by at least four stations; 2)Seismic phases with error greater than ±0.5s were eliminated by using the epicentral distance-travel time fitting curve of the selected earthquake; 3)The distance between each earthquake pair is required to be less than 10km, and the number of double difference data formed by each earthquake pair is required to be larger than 8.

    In the inversion process, the research area is divided into three dimensional grids according to the station location and earthquake distribution. The horizontal direction is divided into 0.3°×0.3° grids, In vertical depth, nodes are set at 0km, 5km, 10km, 15km, 20km, 25km, 30km, 35km, 42km, 50km, and 60km respectively. The value of the damping coefficient is set to 600, the value of the smoothness factor is set to 40, and the number of iterations is set to 10. Thus, after 10 iterations of inversion, the distribution range of residual travel time of seismic data decreases from ±3s to±1s, and the horizontal and vertical errors of the source location after relocation are 0.1~0.9km and 0.1~1.5km, respectively. In order to ensure the accuracy of velocity structure inversion, the reliability of the results is evaluated by using Differential Weighted Sum of Nodes(DWS) and a detection board. Finally, the P-wave velocity structure at depths less than 50km underground in the study area, along with the relocation source parameters of 17613 earthquakes, are obtained.

    The results show that: 1)The focal depth of earthquakes is mainly distributed in the 5~25km depth range. The relocated earthquakes are more closely clustered near the fault zone and the seismic spatial distribution can better describe the geometric morphology of deep faults. There are several NE-trending and NW-trending faults with deep development and steep dip Angle in the Zhang-Bo earthquake zone. Both the Xiadian fault and the Xinhe fault are nearly vertical deep faults with a high dip Angle. 2)The variation of the P-wave velocity had a good correlation with the topography, geomorphology and tectonic environment. Influenced by the surface sediments, the P-wave velocity in shallow crust of the Shanxi fault depression belt, Hebei plain and inter-mountain basin shows low-velocity anomalies. The P-wave velocity in the crust of the junction area of Shanxi, Hebei, and Mongolia is relatively low. The significant low-velocity anomaly of the Zhang-Bo earthquake belt at depth of 42km underground is related to the Destruction of the North China Craton and the upwelling of deep thermal materials. 3)Based on the P-wave velocity structure and seismic relocation results, the fault is developed in the earthquake-prone areas in the northwest of Beijing and its adjacent areas. Most earthquakes occur in the brittle-ductile transition zone between the brittle upper crust and the ductile middle and lower crust. In summary, the seismicity in the northwest Beijing area is closely related to the development of deep and shallow faults and the velocity structure.

    Table and Figures | Reference | Related Articles | Metrics
    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
    SEISMOLOGY AND GEOLOGY    2025, 47 (3): 897-916.   DOI: 10.3969/j.issn.0253-4967.2025.03.20250038
    Abstract245)   HTML18)    PDF(pc) (16398KB)(72)       Save

    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.

    Table and Figures | Reference | Related Articles | Metrics