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
    GEOLOGICAL DISASTERS AND SURFACE RUPTURES OF JANUARY 23, 2024 MS7.1 WUSHI EARTHQUAKE, XINJIANG, CHINA
    ZHANG Bo-xuan, QIAN Li, LI Tao, CHEN Jie, XU Jian-hong, YAO Yuan, FANG Li-hua, XIE Chao, CHEN Jian-bo, LIU Guan-shen, HU Zong-kai, YANG Wen-xin, ZHANG Jun-long, PANG Wei
    SEISMOLOGY AND GEOLOGY    2024, 46 (1): 220-234.   DOI: 10.3969/j.issn.0253-4967.2024.01.013
    Abstract1002)   HTML57)    PDF(pc) (14676KB)(743)       Save

    The MS7.1 earthquake in Wushi, Xinjiang on January 23, 2024, represents the largest earthquake in the Tianshan seismic belt since the 1992 Suusamyr MS7.3 earthquake in Kyrgyzstan. Preliminary precise aftershock localization and initial field investigations indicate an NE-trending aftershock zone with a length of 62km that is concentrated at the mountain-basin transition area. This event produced geological hazards, including slope instability, rockfalls, rolling stones, and ground fissures, primarily within a 30-kilometer radius around the epicenter. The epicenter, located approximately 7 kilometers north of the precise positioning in this study, witnessed a rapid decrease in geological hazards such as collapses, with no discernible fresh activity observed on the steep fault scarp along the mountainfront. Consequently, it is inferred that the causative fault for this main shock may be an NW-dipping reverse fault, with potential rupture not reaching the surface.

    Moreover, a surface rupture zone with a general trend of N60°E, extending approximately 2 kilometers, and displaying a maximum vertical offset of 1m, was identified on the western side of the micro-epicenter at the Qialemati River. This rupture zone predominantly follows the pre-existing fault scarp on higher geomorphic surfaces, indicating that it is not new. Its characteristics are mainly controlled by a southeast-dipping reverse fault, opposite in dip to the causative fault of the main shock. The scale of this 2-kilometer-long surface rupture zone is notably smaller than the aftershock zone of the Wushi MS7.1 earthquake. Further investigation is warranted to elucidate whether or not the MS5.7 aftershock and the relationship between the SE-dipping reverse fault responsible for the surface rupture and the NW-dipping causative fault of the main shock produced it.

    Table and Figures | Reference | Related Articles | Metrics
    PRECISE RELOCATION OF SMALL-TO-MODERATE-SIZED EARTHQUAKES IN THE DATONG VOLCANIC GROUP AND SURROUNDING AREAS
    XU Yong-qiang, LEI Jian-she
    SEISMOLOGY AND GEOLOGY    2024, 46 (2): 336-356.   DOI: 10.3969/j.issn.0253-4967.2024.02.006
    Abstract361)   HTML25)    PDF(pc) (14298KB)(201)       Save

    In the present study we collect a large amount of arrival times from 3 218 earthquakes in the Datong volcanic group and surrounding areas from January 2008 to January 2023 through the China Seismic Network Center and relocated these earthquakes using double-difference location algorithm, finally obtain 2 447 relocate earthquakes. Our result shows that most earthquakes occurred above a depth of 16km, and earthquakes in the basin occurred at depths of 5-16km. There are fewer earthquakes occur near the surface at depths of 0-2km, while 6km and 11km are the dominant depths for earthquakes. The overall strike trending of these earthquake sequences is NE-SW, which is consistent with the regional active faults and controlled grabens and semi-graben-type faulting basins. In addition, these earthquakes are more concentrated near the Kouquan fault zone and in the Datong-Yanggao earthquake zone in the eastern part of the volcanic group. The average location errors of these earthquakes in the east-west, north-south, and vertical directions are about 0.21km, 0.22km, and 0.30km, respectively, with an average travel time residual of 0.14s.

    The earthquakes near the Kouquan fault zone changed from deeper and more concentrated in the south to shallower and more scattered in the north. The earthquake sequences in the northern part of the southern section and the southern part of the middle section of the Kouquan fault zone are deeper along the NE-SW direction, roughly vertically distributed on the Kouquan fault. The earthquake sequences in the northern part of the middle section of the Kouquan fault zone did not occur on the Kouquan fault, and the distribution of earthquakes is relatively scattered, and earthquakes with larger magnitude are mostly concentrated at shallow depth, which may be related to the thick sedimentary coal-bearing strata and mining activities in the area. The strike trending of these earthquakes in the northern section of the Kouquan fault zone is, along the NE-SW direction, roughly distributed on the Kouquan fault. However, there are also earthquakes in the northern part of the Kouquan fault zone, which may suggest that the activity of the Kouquan fault zone has extended there.

    The focal depth in the source areas of the Datong-Yanggao earthquake is mostly concentrated at depths of 3-16km on the hidden fault parallel to the NE-SW trending Dawangcun fault to the east. The hidden fault has a large dip angle and dips towards NW, which intersects with the Tubao fault and the Liulengshan piedmont fault, likely related to the aftershock activity of the Datong-Yanggao earthquake.

    Earthquakes occur frequently in the middle section of the Huairen fault, followed by the southern section, and there are few earthquakes in the northern section. The seismic activity of the Shuiyu fault, the east fault of the Cailiangshan mountains, and the Yanggao-Tianzhen fault is relatively weaker. There are some seismic activities in the central part of the northern margin fault of the Tianzhen-Yanggao Basin. Earthquakes in volcanic areas occurred at the boundaries of volcanic clusters, while the seismicities inside the volcanic group area were not very strong, which suggests that the boundary of volcanic clusters is more prone to stress accumulation and earthquake generation than the interior of volcanic clusters.

    Based on the new seismic results of ambient noise tomography in the area, it is found that earthquakes are not only related to faults, but more importantly, most earthquakes occur near the high-to-low-velocity anomaly boundaries. Furthermore, there are obvious low-velocity anomalies visible beneath most earthquake source areas, which may suggest that the occurrence of these earthquakes is closely related to fluids carried by the upwelling of thermal materials rising to the crust from the mantle and reducing the effective normal stress on the fault planes.

    Table and Figures | Reference | Related Articles | Metrics
    APPLICATION OF DEEP LEARNING IN ACTIVE TECTONICS AND GEOMORPHOLOGY
    LIU Xin, WANG Shi-rou, SHI Xu-hua, SU Cheng, LU Chen-yan, QIAN Xiao-yuan, SUN Qiao-yang, DENG Hong-dan, YANG Rong, CHENG Xiao-gan
    SEISMOLOGY AND GEOLOGY    2024, 46 (2): 277-296.   DOI: 10.3969/j.issn.0253-4967.2024.02.003
    Abstract309)   HTML32)    PDF(pc) (5210KB)(272)       Save

    The research on active tectonics and geomorphology involves extensive sub-topics, including the kinematics of crustal movements, the processes underlying the evolution of landforms, and the associated dynamic mechanisms. These sub-topics are intricately connected with the interactions between the Earth’s endogenic and exogenic processes. In the contemporary realm of the Earth system science, research in active tectonics and geomorphology has become a hot topic for interdisciplinary study. The advancement in big data research coupled with the progressive developments in deep learning technologies has furnished this field of study with a voluminous array of data sources and the requisite analytical tools for technical analysis. In recent years, the application of big data and deep learning technologies in this research field has yielded a series of outstanding results, fostering new research directions, and ushering the discipline into a new phase. In this paper we synthesize existing research to outline the data sources pertinent to the study of active tectonics and geomorphology, including field geological survey, unmanned aerial vehicle (UAV)-based photography, aerial photography, and remote sensing observations. Then, we discuss in-depth examination of the recent innovations progresses in deep learning algorithms, including but not limited to convolutional neural networks(CNNs), deep Gaussian processes, and autoencoders. This article further elaborates on innovative applications of deep learning in the study of active tectonics and geomorphology. These include the identification of changes in glacier extent, monitoring volcanic activity and deformation, recognizing river systems, precise surveillance of landslide events, as well as observations of lithospheric deformation co-seismic surface ruptures.

    Based on the summary of prior studies, this paper showcases a distinct application instance. By employing convolutional neural networks(CNNs)within the realm of deep learning image analysis and utilizing UAV-obtained high-resolution images, we undertake the automated detection of structural fractures in granite rocks in Meizhou island, in the southeast of Fujian province, China. In fault damage zones, structural rock fractures are widely developed, and the study of their orientation, system, and secondary characteristics is of great importance for determining their mechanisms of development and the multi-phase tectonic activity events in the region. Under conventional methodologies, the study of structural fractures in rocks is time-consuming and requires considerable manual effort in conducting exhaustive field surveys and detailed interpretation of cartographic representations. However, the application of deep learning can greatly enhance the efficiency of cartographic work. This application case has improved the classic deep learning framework by developing a CNN model specifically designed for the extraction of complex features and multi-scale rock fractures. This model achieved rapid identification of over 9 000 fractures with varied shapes and complex distributions within 55 minutes, attaining an accuracy of 85% and a recall rate of 89%. These findings demonstrate that deep learning significantly enhances operational efficiency in comparison to manual statistical methods for the automated identification of rock structural fractures, while also maintaining exceptional accuracy in fracture detection. Based on the results identified by deep learning, it can be clearly observed that two sets of fractures, oriented NE and NW, develop on the granite outcrops in the study area. According to previous research and the cross-cutting relationships of the fractures, it is known that NE-oriented fractures formed earlier than NW-oriented fractures, corresponding respectively to the Indosinian Movement and the expansion movement of the South China Sea in the tectonic history of South China. Through the automated extraction of deep learning models, the workload of manual mapping can be greatly reduced, yielding results consistent with actual geomorphological phenomena.

    Table and Figures | Reference | Related Articles | Metrics
    CHARACTERISTICS OF FOCAL MECHANISM AND STRESS FIELD IN THE EASTERN BOUNDARY OF THE SICHUAN-YUNNAN BLOCK
    GUO Xiang-yun, FANG Li-hua, HAN Li-bo, LI Zhen-yue, LI Chun-lai, SU Shan
    SEISMOLOGY AND GEOLOGY    2024, 46 (2): 371-396.   DOI: 10.3969/j.issn.0253-4967.2024.02.008
    Abstract281)   HTML20)    PDF(pc) (11846KB)(244)       Save

    It is important to study the characteristics of the tectonic stress field studies which could provide a deeper understanding of the internal stress environment of the crust. It can provide useful assistance for exploring the relationship between the tectonic stress field and earthquake development. At the same time, it plays an important role in understanding block interactions, fault movement, tectonic deformation, and revealing the dynamic mechanical processes of the continent. The focal mechanism solutions contain abundant information reflecting the stress field.

    In this paper, using the broadband records from 128 permanent and temporary regional stations from the Chinese National Seismic Network(CNSN)deployed in the Sichuan-Yunnan Province and its adjacent, we determined the focal mechanisms of 3 951 earthquakes by the cut-and-paste(CAP)method and the HASH method. The friction coefficient and stress properties of the main active fault and characteristics of the tectonic stress field in this area are analyzed by using two different methods which are the damped inversion method(STASI)and iterative joint inversion method from focal mechanisms.

    The results of the focal mechanisms show that: there are 2 512 strike-slip earthquakes in the study area, accounting for 63.58% of all earthquakes; there are 818 normal fault type and normal strike-slip type earthquakes, accounting for 20.70% of all earthquakes; there are 621 reverse strike slip and reverse thrust earthquakes, accounting for 15.72% of all earthquakes. The most of earthquakes in the study area are distributed in active fault zones, the strike of the fault plane is consistent with the orientation of active fault zones. It revealed predominantly strike-slip faulting characteristics of earthquakes in the Eastern Boundary of the Sichuan-Yunnan Block, while the reverse thrust of earthquakes is mainly concentrated in the Longmenshan fault zone, as well as the NW trending Mabian-Yanjin Fault and the NE trending of Ludian-Zhaotong and Lianfeng faults which lied on the eastern boundary of the Sichuan-Yunnan block. Overall, the characteristics of the source mechanism are consistent with the regional tectonic background.

    Results of the stress field inversion confirmed main active fault in the Eastern Boundary of the Sichuan-Yunnan Block is under a strike-slip stress regime, maximum and minimum compressional stress axes are nearly horizontal. The maximum compressional axes are primarily oriented in NW-SE and NWW-SEE direction, and they experience a clockwise rotation from north to south. Against the strike-slip background, normal faulting stress regimes and reverse faulting stress can be seen in the regional areas. The most prominent is the Daliangshan fault zone, which has obvious differences from the overall characteristics of the stress field with the eastern boundary of the Sichuan Yunnan Block. The maximum horizontal principal stress in the northern section shows a nearly EW direction, with a strike-slip type stress property, and the NW-SE direction in the southern section, with a thrust type stress property. The distribution characteristics of the stress field are consistent with the fault type of sinistral strike-slip and thrust on the eastern boundary of the Sichuan Yunnan block

    The shape ratio R-value varies significantly, the R-value in the Sanchakou area is relatively high, with obvious extrusion characteristics, the R-values of the Xianshuihe fault zone, Anninghe fault zone and Xiaojiang fault zone are all between 0.25-0.5, showing NE-SW compression and NW-SE tension, and the tensile stress may be much less than the compressive stress(strike-slip type). The R values of the northern segment of the Daliangshan fault zone, the southern segment of the Anninghe fault zone, and Zemuhe fault zone are all between 0.5-1, showing NW-SE compression and NE-SW tension, and the compressive stress is greater than the tensile stress. To sum up, the current stress characteristics of the eastern boundary of the Sichuan Yunnan rhombic block are shear strain and local compression or tension.

    There are different friction coefficients of the main faults in the study area: The Anninghe fault zone is 0.60, the Xianshuihe and Zemuhe fault zones are 0.80, the Xiaojiang fault zone is 0.75 and northern and southern sections of the Daliangshan fault zone are 0.65 and 0.85. The friction coefficients of the Xianshuihe Fault, the southern section of the Daliangshan Fault, and the Zemuhe Fault are above 0.75. The high friction coefficients of these fault zones may be because they are strike-slip faults, and the friction coefficients themselves are relatively high. The southern section of the Xiaojiang fault zone may be related to the development of fault gouges in the fault zone.

    Table and Figures | Reference | Related Articles | Metrics
    GEOMORPHIC DATING OF SCARPS AND ITS APPLICATION TO ACTIVE TECTONICS AND GEOMORPHOLOGY
    PANG Zhen-hui, XU Hao-ting, SHI Xu-hua, GE Jin, LI Feng
    SEISMOLOGY AND GEOLOGY    2024, 46 (2): 251-276.   DOI: 10.3969/j.issn.0253-4967.2024.02.002
    Abstract262)   HTML28)    PDF(pc) (3088KB)(292)       Save

    Scarps are typical geomorphic features of tectonics, climatic changes, and erosion processes. On one hand, interpreting geological information encoded in scarps allows for the quantitative constraint of the kinematic and dynamic mechanisms of the active structures. On the other hand, studying the evolution processes of scarps contribute to a better understanding of the couplings among tectonics, erosion, and climate during geomorphic evolution processes. In regions characterized by adverse geological conditions, limited accessibility, and logistical challenges hindering researchers from reaching certain areas, traditional dating methods such as radiocarbon dating, luminescence dating, and cosmogenic nuclide dating often face difficulties in determining the age of scarps. The geomorphic dating method of scarps, however, offers a promising avenue to address the scarcity of chronological samples in research areas where either sample availability is limited or conventional dating techniques are impractical. This paper provides a concise summary of the theoretical evolution of geomorphic dating of scarps. Emphasis is placed on elucidating the slope evolution processes, transport models, and associated computational methodologies integral to this approach. Additionally, the specific applications of these methods in active tectonics and geomorphology are highlighted, accompanied by a case study showcasing their practical implementation.

    The theoretical foundation of geomorphic dating of scarps posits that the evolution of scarps during stable erosion stages can be simulated through models describing the evolution of slope surfaces over time. In practical dating applications, it is essential to determine the theoretical models and computational methods for the evolution of scarps. This necessitates the integration of measured profiles of the scarp to establish boundary and initial conditions, facilitating the determination of the geomorphic age of the studied scarps. On one hand, the related slope evolution model mainly involves processes such as bedrock weathering, sediment transport, and tectonic uplift. Previous studies have proposed dozens of quantitative slope evolution models and geomorphic transport functions(e.g., local linear, local nonlinear, non-local, etc.)based on various slope processes, theoretical assumptions, and numerical simulations. In various transport equations, compared to earlier local linear models, later local nonlinear transport models proposed based on experimental simulations and physical derivations exhibit higher fitting accuracy for real slope evolution. In the past decade, some scientists have proposed nonlocal transport models because of the limitations of traditional transport models, and have applied them in research. This nonlocal model assumes that the distance of sediment movement within a given area follows a probability distribution, thus allowing the simulation of long-distance slope processes over short periods. Additionally, many other transport models have been derived from specific slope processes, such as biotic disturbance and dry ravel. The solution methods for the aforementioned models vary as well. For instance, the analytical solution of a local linear diffusion transport model can be relatively easily obtained, while local nonlinear models and nonlocal models can only be numerically solved through specific approaches. On the other hand, the measured topographic profiles of the studied scarps can be used to determine the practical parameters of slope evolution models, including the present-day morphology of the scarps and their ages since their initial formation. In practical applications, various methods have emerged for the geomorphic dating of scarps, generally classified into two types based on the approach to fitting model calculations with actual topographic profiles: the mid-point slope method and the full slope method. The mid-point slope method uses the mid-point gradient value as the fitting morphological feature, representing an early method for dating scarps, mostly combined with linear diffusion transport functions and requiring numerous profiles for statistical analysis. Due to its low data utilization and limited spatiotemporal precision in statistical methods, the mid-point slope method has gradually been replaced by the full slope method. The full slope method involves fitting the overall shape of actual profile curves using model solutions. With the continuous improvement of observation techniques in the field of Earth sciences and the deepening research on related theories, the application scope of scarps geomorphic dating methods is no longer limited to the study of terraces and simple fault scarp evolution processes but has expanded to more complex geological environments, providing more precise constraints on their formation and evolution history.

    For method application, we systematically present the progress in scarp geomorphic dating research across various geomorphic settings(such as river and coastal terraces, lake shorelines, alluvial fans, marine terraces, and extraterrestrial planets). It employs the geomorphic dating of the northeastern Pamir fault scarp as a case study to further explore and anticipate the developmental trajectory of geomorphic dating of scarps within the field of tectonic geomorphology.

    Table and Figures | Reference | Related Articles | Metrics
    INTEGRATED INTERPRETATION ON THE PRECURSORY PROCESS EVOLUTION IN THE META-INSTABILITY STAGE OF THE EARTHQUAKE: A CASE STUDY ON 2014 LUDIAN MS6.5 EARTHQUAKE
    JIANG Hai-kun, DENG Shi-guang, YAO Qi, SONG Jin, WANG Jin-hong
    SEISMOLOGY AND GEOLOGY    2024, 46 (3): 513-535.   DOI: 10.3969/j.issn.0253-4967.2024.03.001
    Abstract253)   HTML37)    PDF(pc) (5232KB)(167)       Save

    The transition from the metastable state to the meta-instability stage indicates that the seismic fault has entered an irreversible deformation process and will lead to an inevitable instability(Ma Jin et al., 2014). Therefore, identifying the meta-instability stage is helpful for the judgment of short-term earthquake precursor anomalies. Under laboratory conditions, the meta-instability stage can be visually identified through stress-time curves, thus potentially predicting the occurrence of the laboratory earthquake. However, there are significant differences between field conditions and laboratory environments. Firstly, the underground medium and structural conditions in the real earthquake source region are unclear and far more complex than laboratory specimens. Secondly, the distribution of the sensors, sensor density, as well as measurement accuracy are limited by various conditions, making it impossible to construct an ideal observation environment covering the entire region. Thirdly, the loading stress cannot be directly measured, and the current actual stress state of the study area is unknown, which is the most difficult problem to solve. Therefore, under the guidance of meta-instability experiments and theories, it is a beneficial attempt to conduct retrospective studies on typical earthquake cases with relatively good observation conditions in the past, analyze the spatial-temporal evolution of different physical fields at different stages before the earthquakes, compare the observed phenomena with the characteristics and change processes of meta-instability stages obtained from experiments or theoretical research. Its final goal is to find possible characteristics or indirect criteria for meta-instability stages under field observation conditions.

    Therefore, taking the Ludian MS6.5 earthquake as an example and under the guidance of the meta-instability experiments and theories, the paper comprehensively analyzes the relationship between the spatial-temporal evolution of precursory anomalies and the meta-instability process based on the seismic activity and the geophysical observation data prior the earthquake, and combined with numerical simulation results of the earthquake nucleation. The Ludian MS6.5 earthquake occurred on August 3, 2014 in northeastern Yunnan Province, China. The observation conditions in this region were relatively good, with 26 seismometers within a 300-km radius of the epicenter, which were able to basically monitor earthquakes with completeness magnitude ML≥1.5 and locating accuracy of less than 20km. There were 79 fixed geophysical observation stations, including 11 within 100km, 32 between 101~200km, and 36 between 201~300km. The observation terms covered 43 deformation observations(22 tilt observations, 18 borehole strain observations, and 3 gravity observations), 187 underground fluid observations(90 water physical observations such as water level and temperature, 43 material compositions measurements including radon, mercury and so on, 26 gas measurements such as CO2, and 28 ion measurements including bicarbonate, calcium, and magnesium), and 52 electromagnetic observations(36 geomagnetic observations, 16 resistivity and electromagnetic wave observations). There were a large number of credible medium- and short-term precursor anomalies before the Ludian MS6.4 earthquake, a total of 48 precursor anomalies were identified. Among of them, there were 8 seismic anomalies and 40 geophysical anomalies, accounting for approximately 15% of all measurement items. Among these 40 geophysical anomalies, 31 were proposed before the earthquake, and most of them were investigated and verified on-site with reliable changes.(Wu, et al., 2019).

    Based on this abundant precursor abnormally data before the Ludian MS6.4 earthquake and further systematic analysis, a typical earthquake case and relevant observational facts have been provided which can support the viewpoint that during the meta-instability stage, the earthquake nucleation occurred in the epicenter region and the synergy process evolved continuously in surrounding area of the epicenter. The results show that based on large-scale strong earthquake activities and the observation data of the mobile gravity, it can be determined that the concerned area was already in a high-stress state before the Luding earthquake. At that time, the stress level in the large area including the epicenter of the Ludian earthquake was relatively high, and the northeastern Yunnan region and its nearby areas where the Ludian earthquake occurred were already in a critical stress state where strong earthquakes could occur at any time. Under the premise of determining a high-stress state, according to the precursors of seismic activities and geophysical observation precursor anomalies, it can be roughly determined that the meta-instability process of the Ludian earthquake may have begun seven or eight months before the mainshock. The most prominent phenomenon or judgment index is the transition of the fault stress state from accumulation to release, characterized by the active of small earthquakes near the epicenter, as well as the synergistic phenomenon of fault deformation characterized by the significant increase in the number of geophysical observation anomalies, which is related to the expansion process of the core weakening zone in the late period of the earthquake nucleation. After that, until the occurrence of the mainshock, two times should be paying attention to. Firstly, four to five months before the mainshock, the spatial distribution range of the geophysical observation anomalies expands significantly from the epicenter area to the periphery region, indicating accelerated synergistic deformation of the fault. Secondly, after two months before the mainshock, the small earthquake activities near the epicenter began to weaken, and the micro-earthquake activities and the geophysical anomalies showed a migration and contraction towards the epicenter, which is associated with the contraction process of the core weakening zone during the final stage of the earthquake nucleation. The concept of seismic meta-instability proposed from the perspective of stress changes in seismic fault has an explicit physical implication, and the meta-instability stage is associated with the earthquake nucleation process(He et al., 2023). The basic premise for the meta-instability theory to play a role in short-term earthquake prediction lies in how to apply laboratory research results to natural earthquakes, understand whether the regional or fault stress state tends to or enters a meta-instability state through field observations, and further utilize it for practical earthquake prediction.

    Table and Figures | Reference | Related Articles | Metrics
    THE DISASTER MECHANISM OF THE MS6.9 EARTHQUAKE IN MENYUAN, QINGHAI PROVINCE, 2022
    NIU Peng-fei, HAN Zhu-jun, GUO Peng, LI Ke-chang, LÜ Li-xing
    SEISMOLOGY AND GEOLOGY    2024, 46 (4): 761-782.   DOI: 10.3969/j.issn.0253-4967.2024.04.001
    Abstract249)   HTML28)    PDF(pc) (22134KB)(134)       Save

    Earthquake disasters are one of the most significant natural disasters faced by human society. Understanding and mitigating earthquake disasters have always been a key focus of research for seismologists. Conducting investigations on post-earthquake seismic disasters is of great significance for the recovery and reconstruction of disaster-stricken areas, as well as for earthquake prevention and mitigation. Earthquake disasters can be classified into two types based on their mechanisms: one is the destruction caused directly by the seismic vibrations on buildings, lifelines, and other structures; the other is the damage related to geological hazards triggered by earthquakes. The former is mainly related to the density of regional economic layout; the latter seismic geological disasters typically include collapses, landslides, debris flows, ground fissures, ground subsidence, and soil liquefaction. These geological disasters often exacerbate the impact of seismic disasters, posing a more significant threat to human life and property safety. Therefore, it is of great significance to investigate the mechanisms of significant engineering disasters caused by earthquakes, as it can provide important insights for engineering recovery, reconstruction, and site selection. The Qilian-Haiyuan fault zone is an important boundary fault on the northeastern margin of the Qingzang Plateau. It plays a crucial role in absorbing and accommodating the convergence of the Indian Plate towards the Eurasian Plate in a NNE direction. With a total length of approximately 1 000km, it is primarily composed of the Tolaishan fault, the Lenglongling fault, the Jinqianghe fault, the Maomaoshan fault, the Laohushan fault, and the Haiyuan fault, from west to east. On January 8, 2022, a magnitude 6.9 earthquake occurred near the stepover of the Longling and Tuolaishan faults of the Qilian-Haiyuan fault zone. Although the earthquake occurred in uninhabited, sparsely wooded alpine grasslands and did not cause any casualties, it completely destroyed the Liuhuanggou bridge and the south-side Daliang tunnel on the Lanzhou-Xinjiang high-speed railway, a major artery of the Silk Road transportation network in China. This marks the first time that a mainline of the high-speed railway network, which is a showcase of China's economic achievements, has been entirely disrupted by earthquake damage. Based on the high-resolution orthophoto images and digital elevation models(DEMs)obtained through post-earthquake emergency scientific investigations using the unmanned aerial vehicles, this article conducted another field investigation on earthquake disasters in vehicles; this article conducted another field investigation on earthquake disasters in the isoseismal area. First, by investigating geological disasters such as collapses, landslides, and soil liquefaction in the meizoseismal area, as well as the damage to buildings and structures. Then, based on field surveys, a detailed mapping of the reverse-type surface ruptures formed by the Mengyuan earthquake was conducted, identifying the distribution patterns and geometric and kinematic characteristics of the surface ruptures and determining the distribution of coseismic vertical displacements. Additionally, the development of geological disasters caused by this earthquake was analyzed, and the disaster-causing mechanism of the Liuhuang Bridge was discussed. The research indicates that the Liuhuanggou River, located in the isoseismal area, does not exhibit large-scale earthquake landslides and collapses. Instead, only smaller-scale rockfalls and accumulations of rolling stones, as well as localized occurrences of sand liquefaction in certain riverbeds, are observed, which is clearly inconsistent with expectations. In addition to the formation of two strike-slip surface rupture zones, the earthquake also generated a reverse-type surface rupture zone approximately 7.9km long within the Liuhuanggou river on the northern side of the western section of the Lenglongling fault. The rupture zone exhibits an unstable southward trend and is primarily composed of discontinuous arc-shaped compressional ruptures, mole tracks, tensile ruptures, and seismic scarps. Along the surface rupture zone, a total of 35 vertical displacement measurements were obtained, with the minimum displacement of (8±1)cm and the maximum displacement of (49±3)cm. The average vertical displacement is approximately 24cm, and the displacement distribution along the strike is uneven. The surface rupture zone, which cuts nearly vertically across the Lanzhou-Xinjiang high-speed railway Liuhuanggou bridge, has caused extensive surface deformation and displacement. This is the direct cause of the destruction of the Liuhuanggou bridge. This finding suggests that when implementing seismic engineering design measures for major linear projects crossing fault zones, it is important to consider the extensive shear effects of reverse-type surface rupture zones.

    Table and Figures | Reference | Related Articles | Metrics
    OPPORTUNITIES BROUGHT BY 3D GEOSCIENCE MODELING FOR EARTHQUAKE NUMERICAL FORECASTING
    YAO Qi, LU Ren-qi, SU Peng, WANG Hui, ZHU Ya-ling, WANG Li-wei
    SEISMOLOGY AND GEOLOGY    2024, 46 (1): 1-18.   DOI: 10.3969/j.issn.0253-4967.2024.01.001
    Abstract246)   HTML30)    PDF(pc) (1599KB)(226)       Save

    Earthquake prediction and forecasting need to transform from the traditional empirical, qualitative, and semi-quantitative to quantitative. The improvement also calls for multi-disciplinary, highly integrated physical and mechanical simulations rather than only a single discipline. The global development of observation technology and the construction of observation networks have already built a data foundation for earthquake numerical prediction and forecasting to a certain extent. However, the biggest constraint is the difficulty of synthesizing a large amount of observation data and quickly establishing complex numerical models with geological significance for numerical calculation. It is a vital issue restricting experimental research and industry development of earthquake numerical prediction and forecasting. Based on a brief introduction of the concept, development, and research status of earthquake numerical prediction and forecasting, this paper analyzes the difficulties in numerical modeling, which essentially come from the disciplinary differences between active tectonics, structural geology, solid earth, seismology, and numerical simulations. The development of 3D Geoscience Modeling and its application in the earthquake industry can establish a large-scale complex earthquake tectonic model close to the real world with geological significance. It provides a significant opportunity and technical means for developing earthquake numerical prediction and forecasting by solving the problems in numerical modeling. 3D geological modeling has built a bridge for multi-disciplinary geological applications. It can multi-disciplinary data fusion, establish a 3D geological model with geological significance and characteristics in line with geomechanical characteristics, and integrate data, geological model, up to building a numerical model, which advances the efficiency of modeling and simulation. Therefore, the rapid development of 3D geological modeling provides an opportunity to solve the modeling difficulties mentioned above in earthquake numerical prediction. Then, we briefly describe the development of 3D geological modeling technology, its application in the seismic industry, and the construction and application of 3D standard fault models domestically and overseas. Here, we introduced the development and essential contents of the Community Fault Model of Southern California in the United States for the Uniform California earthquake rupture forecast, the New Zealand Community Fault Model from the Institute of Geological and Nuclear Sciences Limited, and the Community Fault Model in Sichuan and Yunnan region in China.

    The prospective future of 3D geological modeling and its potential application in earthquake numerical prediction and forecasting makes it a common concern of researchers in earthquake science. The five future modeling trends are the joint modeling of multi-source and multi-precision heterogeneous data, the integrated modeling of the geological model-attribute model-numerical model, flat fault structure modeling, 3D fault structure modeling, data-model-calculation iteration, and mutual driving. Finally, the paper describes the difficulties of applying the 3D geological modeling technique in earthquake numerical prediction and forecasting, including the industry construction, public approval of the 3D Community Fault Model, and the variations of numerical modeling and applications. 3D geological modeling technology can provide more realistic numerical and geometric models for earthquake numerical prediction, forecasting, and related numerical computing fields, reduce construction periods, create fast iterations, and solve modeling difficulties.

    Table and Figures | Reference | Related Articles | Metrics
    JOINT INVERSION OF THE RUPTURE PROCESS OF 2018 ML5.7 XINGWEN EARTHQUAKE BASED ON SEISMIC AND INSAR OBSERVATIONS
    MIAO Si-yu, ZHANG Hai-jiang, GU Ning, LI Jun-lun, TAN Yu-yang, HUA Si-bo, ZHANG Yong
    SEISMOLOGY AND GEOLOGY    2024, 46 (2): 397-413.   DOI: 10.3969/j.issn.0253-4967.2024.02.009
    Abstract234)   HTML15)    PDF(pc) (5595KB)(168)       Save

    The ML5.7 Xingwen earthquake on December 16, 2018 is very likely induced by shale gas hydraulic fracturing, which caused not only massive landslides and rock collapse, but also some casualties in the surrounding area, with the direct economic loss of about 50 million CNY. It is of great significance to determine the source rupturing process of such an induced earthquake with large magnitude.

    Finite fault inversion is one of the commonly adopted methods to determine coseismic slip displacement distribution. For finite fault inversion, various data have different sensitivities to various aspects of the rupture process. The seismic data can provide the historical information about the earthquake rupture process because it contains the Doppler effect of the space-time rupture behavior on the fault. In comparison, the near-field geodetic data(such as InSAR and GPS)can constrain the fault parameters and the static slip distribution well because they contain the surface motion information. Therefore, the reliability of the inversion for the complex seismic rupture process can greatly be improved by combined use of seismicdata and InSAR data.

    In this study, strong-motion seismic data recorded at 8 near-field stations are chosen and filtered by a band-pass of 0.15-0.60Hz. The same InSAR data used in Wang et al.(2022)is adopted in this joint study. For inversion, a sufficiently large potential fault plane of 15km long and 10km wide is chosen and divided into 15×10 subfaults. Finally, the rupture process is obtained by joint inversion of strong-motion seismic data and InSAR data. The results show that the earthquake is characterzied by a typical unilateral rupture with the rupturing direction nearly towards the north. The duration of the rupture process was 6s, and the energy release was mainly concentrated in the first 5s. The rupture process is segmented and can be divided into two stages. The first stage is distributed from 1-3s and is located in the range of 0~5km from the source; and the 2nd stage is distributed from 3-5s and is located between 6 and 8km from the source. The coseismic slip is mainly concentrated in areas shallower than 5km, with a peak slip of approximately 0.27m. This can be used to explain why the Xingwen earthquake with a magnitude of ML5.7 caused relatively serious damages.

    Combined with the distribution of foreshocks and aftershocks, it can be seen that the foreshocks were mainly concentrated to the eastern edge of the major coseismicslip zone, which are close to some hydraulic fracturing wells. This suggests that these foreshocks occuring at the edge of the main rupture zone has a certain correlation with fluids, and the presence of fluids further leads to the fault weakening of the mainshock due to the increase of pore pressure and the decrease of effective compressive stress, which plays a triggering role in the occurrence of the Xingwen earthquake. The aftershocks are mainly distributed around the main slip zone, which are caused by after slips after the mainshock. The results from seismic inversion, InSAR inversion and joint inversion of the two data types reveal that the Xingwen earthquake is a northward unilateral rupture. The rupture propagation direction and coseismic slip distribution may be related to the physical property changes along the fault plane.

    Compared with the two single inversion results, the joint inversion overcomes the influence of uneven distribution of seismic stations, improves the resolution of slip distribution, and produces results that are more consistent with the real physical process. The slip model obtained by joint inversion in this study can be helpful for further understanding the mechanisms of induced earthquake, the correlation between induced earthquake and geological structure, earthquake disaster assessment and post-earthquake disaster prevention and hazard mitigation.

    Table and Figures | Reference | Related Articles | Metrics
    3D STRUCTURAL MODELLING OF THE ANNINGHE-ZEMUHE-XIAOJIANG FAULT ZONE IN THE EASTERN BOUNDARY OF SICHUAN-YUNNAN BLOCK USING MULTI-DATA AND IMPLICIT MODELING METHODS
    WANG Mao-mao, HU Shun-yang, MA Hao-ran, LIANG Bo-yu, ZHANG Jin-yu, LU Ren-qi
    SEISMOLOGY AND GEOLOGY    2024, 46 (1): 19-34.   DOI: 10.3969/j.issn.0253-4967.2024.01.002
    Abstract211)   HTML21)    PDF(pc) (8144KB)(179)       Save

    The Anninghe-Zemuhe-Xiaojiang fault zone is located at the intersection between the Qinghai-Xizang Plateau and the Yangtze block, representing the eastern boundary of the Sichuan-Yunnan block with frequent seismic activities. Its overall kinematic characteristics involve left-lateral strike-slip motion, and the fault structures along its strike are complex, posing significant challenges in accurately characterizing the 3D structural features of deep faults. The main issues include the structural complexity of the fault surfaces, uncertainties in the intersection relationships of fault systems, spatial constraints of blind faults, and the definition of fault surfaces in regions with weak seismic activity. Traditionally, 3D structural modeling for fault geometry heavily relies on high-resolution seismic reflection profiles, 3D seismic data volumes, and borehole data. It defines the geometric shapes of objects with limited nodes in a triangular mesh, and then simulates the topological structure of objects by connecting these nodes. However, obtaining high-resolution seismic reflection data in active tectonic areas like the eastern boundary of the Sichuan-Yunnan block is challenging, and even when available, it is often sparse in space. Alternatively, a large amount of relocated earthquakes and surface fault traces are generally used to create initial three-dimensional models of active faults. However, this approach overlooks the contributions of focal mechanism solutions in constraining the modeling, with more subjectivity in the selection of relocated seismicity, and does not adopt a differentiated weight strategy for various data sources. In this study, a 3D implicit modeling approach, combining deep and shallow geological and geophysical data that are generally available in active tectonic environments, was used to construct a detailed 3D structural model of the Anninghe-Zemuhe-Xiaojiang fault zone at the eastern boundary of the Sichuan-Yunnan block. The modeling process effectively integrated the fault plane constraints provided by focal mechanism solutions with surface fault traces and relocated seismic data, using a multi-iteration process with differentiated weight to increase the accuracy of the fault models. This approach ultimately represented the 3D complex structural features of the eastern boundary of the Sichuan-Yunnan block using multiple data sources. The modeling results show that the Anninghe-Zemuhe fault zone is characterized by a steep strike-slip fault structure with along-strike geometry variations. The Anninghe Fault shows its steepest dip angle in the central segment and gradually becomes gentler to both ends. Meanwhile, the Zemuhe Fault exhibits several asperities that are perpendicular to the direction of fault slip at a depth of 5~15km. By contrast, the north-to-central segment of the Xiaojiang fault zone is more complex. The western branch of the Xiaojiang Fault, which is an east-dipping, left-lateral strike-slip fault, is characterized by a relatively gentle fault plane with an average dip angle of 76° to 78°. The west-dipping segment of the eastern Xiaojiang Fault has a steeper dip with an average angle of 85°. The detailed 3D structural model of active faults constructed through implicit modeling can be used for analyzing fault roughness and fault system studies, which are crucial for understanding the distribution of asperities on fault planes and conducting seismic rupture simulations. Implementing the implicit modeling approach allows for the development of improved fault surface representations that can contribute to Community Fault Models in active tectonic environments, and support fault system modeling, rupture simulations, and regional hazard assessments.

    Table and Figures | Reference | Related Articles | Metrics
    EFFECTS OF STRIKE-SLIP FAULT GEOMETRIC COMPLEXITY ON EARTHQUAKE RUPTURES PROPAGATION: A REVIEW
    WANG Hui, CAO Jian-ling, YAO Qi, WANG Li-wei, ZHU Ya-ling
    SEISMOLOGY AND GEOLOGY    2024, 46 (1): 201-219.   DOI: 10.3969/j.issn.0253-4967.2024.01.012
    Abstract208)   HTML21)    PDF(pc) (1897KB)(200)       Save

    Active strike-slip fault usually hosts major earthquakes. Therefore, studies on fault segmentation, which is controlled by geometric complexity, are very important for assessing the maximum magnitude of potential earthquakes. Based on previous literature, we summarized the behavior of earthquakes on strike-slip faults related to fault geometry, segmentation, and cascading rupture from the aspects of field investigation and numerical modeling.

    Previous field investigations have shown that geometrically complex sections of a strike-slip fault usually play the role of barriers that can separate earthquake rupture segments and effectively stop the propagation of earthquake rupture. Statistical results based on the limited field investigations of the surface rupture illuminated semi-quantitatively the influence of geometrically complex sections on the rupture behavior of the earthquake. Since not all earthquakes can produce surface rupture zones, the case number of surface ruptures are unlikely to meet statistical requirements in the coming years. In addition, knowledge gained from field investigation is mainly statistical results based on fault geometry and kinematics. They have some consistency but vary greatly, indicating the complexity of seismic rupture modes on strike-slip faults. No simple threshold that can be used as a criterion to refine the capability of earthquake rupture propagation on strike-slip faults with complex geometry. Based on the statistical results of field surveys, geologists have applied the factor, that complex geometricity controls earthquake rupture behavior, in seismic risk assessment. Lacking support from dynamic analysis, it is necessary to develop and integrate physics-based dynamic models to help improve earthquake-rate models and probability models.

    Numerical modeling can not only present the earthquake simulation scenario but also provide insights into the fundamental physics of dynamic rupture propagation. The modeling results improve our understanding of how the geometric complexity of the fault influences the dynamic rupture propagation. Different modeling approaches focus on different aspects of this challenging scientific problem, each with unique advantages and disadvantages. 2D modeling is relatively simple. They allow modelers to consider more complex physical processes, variated parameters, and constraints from the field and laboratory observation, etc. They provide a comparative benchmark on rupture dynamics on a strike-slip fault with complex geometry. 3D modeling can provide closer approximations to realistic faults and more direct comparisons to observations. The simulations of one earthquake rupture process may focus on the influence of single/multiple parameters on the rupture process. While the multicycle earthquake simulation can predict spatiotemporal patterns of earthquake ruptures on a strike-slip fault. Both simple 2D modeling and complex 3D modeling have shown that one of the most important factors affecting rupture behavior on strike-slip faults is the fault geometric complexity. In addition, other dynamic factors, such as the initial stress, the properties of the rock medium, and nucleation location, also affect dynamic fracture propagation on strike-slip faults. It indicates that rupture behavior on a certain strike-slip fault has its unique characteristics that are controlled by dynamic factors such as the regional tectonic environment and the properties of the fault itself. The numerical modeling provides a dynamic perspective on the complexity of rupture behaviors on strike-slip faults given by field investigations.

    In the China Seismic Science Experimental Site, 3D dynamic modeling supported by fault detection, dense geophysical observations, and high-performance computation will provide new insights into the rupture behavior in the complex multi fault system. That is helpful in determining the maximum magnitude of a potential earthquake.

    Table and Figures | Reference | Related Articles | Metrics
    FRICTIONAL PROPERTIES OF SERPENTINE MINERALS UNDER HYDROTHERMAL CONDITIONS
    LIU Shi-min, ZHANG Lei, HE Chang-rong
    SEISMOLOGY AND GEOLOGY    2024, 46 (2): 235-250.   DOI: 10.3969/j.issn.0253-4967.2024.02.001
    Abstract208)   HTML26)    PDF(pc) (2735KB)(154)       Save

    Serpentine minerals are among the minerals commonly found in the Earth’s subduction zones, and their unique physicochemical properties have a significant impact on subducting geodynamics. Friction experimental studies of serpentine minerals are essential to gain a deep understanding of the frictional sliding stability of serpentine-containing faults in subduction zones as well as explaining the complicated misalignment behavior of faults in subduction zone. Previous laboratory research has produced an abundance of results, and this work addresses two main aspects: the stable states of occurrence and interconversion relationships of serpentine minerals, and the parameters affecting the frictional strength and sliding stability of serpentine minerals. First of all, studies on the stable endowment state of serpentine minerals and the interconversion relationship show that different types of serpentines diaplay different stable phases under different conditions. Chrysotile and lizardite are stable at low temperatures, and the stability fields of both chrysotile and lizardite roughly overlap, but chrysotile is in a substable state. Antigorite is stable at high temperature conditions, such as subduction zone mantle wedges containing high pore fluid pressure conditions, and undergoes a transition from lizardite to antigorite with increasing temperature. Secondly, studies on the factors controlling the frictional strength and sliding stability of serpentine minerals have shown that temperature, pore fluid, and the effective normal stress are all critical factors, for example, an increase in temperature can significantly increase the frictional strength of lizardite and chrysotile. In addition, the friction strength of serpentine minerals shows an obvious pressure dependence, and it was found through previous experimental studies that the friction strength of chrysotile exhibits a high-pressure sensitivity, and that the friction strength of antigorite gradually increases with increasing temperature under low fluid pressure conditions, showing an obvious temperature strengthening phenomenon. In contrast, the change in frictional strength of antigorite with temperature under high-pressure fluid pressure conditions is diametrically opposed to the results of low-pore fluid pressure conditions, which shows a clear temperature weakening phenomenon. Previous studies have also found that antigorite-undergoes a dehydration reaction with increasing temperature under lower fluid pressure conditions, and then exhibits unstable velocity weakening phenomenon, while antigorite exhibits velocity weakening phenomenon under low shear deformation rate under high-pressure fluid conditions. By analyzing the variation of friction-slip stability of antigorite with the shear slip rate can help us to better explain the phenomenon of subduction-zone slow-slip. Overall, experimental studies of the friction of serpentine minerals provide a key experimental basis for a deep understanding of subduction zone geologic processes. The results of these studies are scientifically important for predicting earthquakes and explaining the evolution of the Earth’s internal tectonics and subduction zones, providing strong support for research and practice in the field of geosciences.

    Table and Figures | Reference | Related Articles | Metrics
    THE DEEP SEISMIC REFLECTION PROFILE UNVEILS FINE STRUC-TURE AND TECTONIC CHARACTERISTICS OF THE CANGXIAN UPLIFT, HUANGHUA DEPRESSION, AND ADJACENT AREAS
    QIN Jing-jing, LIU Bao-jin, FENG Shao-ying, XU Xi-wei, TIAN Yi-ming, ZHU Guo-jun, ZUO Ying
    SEISMOLOGY AND GEOLOGY    2024, 46 (3): 608-626.   DOI: 10.3969/j.issn.0253-4967.2024.03.006
    Abstract200)   HTML14)    PDF(pc) (13832KB)(99)       Save

    A comprehensive seismic profiling study was conducted across the Jizhong depression, Cangxian uplift, Huanghua depression, and Chengning uplift in the North China Plain to investigate crustal fine structure and the relationship between deep and shallow faults. Two profiles were completed: a deep seismic reflection profile spanning approximately 200km and a middle-shallow seismic reflection profile covering about 66km.

    Our results indicate a crust thickness of approximately 30 to 35km along the section, with a thin distribution in the east and thick in the west. Notably, there is a significant uplift on the Moho surface beneath the Jizhong depression, with an uplift amplitude of about 2 to 3km. The deep seismic reflection profile reveals distinct upper and lower structural characteristics of the crust. The upper crust displays typical sedimentary layer reflection characteristics, marked by alternating depressions and uplifts. Numerous large-scale faults are concealed beneath the North China Plain, playing a pivotal role in uplift and sag formation. The lower crust’s reflection structure comprises events with significant changes in reflection energy, unstable stratification, and complex occurrences, contrasting with the strong reflection energy and good horizontal continuity of the upper crust reflections. The piedmont fault of the Taihang Mountain, identified by the mid-shallow seismic profile and deep seismic reflection section, is a large shovel-shaped normal fault with a low angle, linked to the large detachment structure at the eastern margin of Taihang Mountain. The detachment structure is developed between the basement and the sedimentary cover layer, and is shown on the profile as a reflection zone consisting of 3 to 4 strong reflection phases, lasting 0.3~0.4 seconds. This detachment structure controls the formation of graben and horst structures. The Jizhong depression may have been an extensional tectonic system formed in the upper crust on the basis of the extensional detachment of the Taihang Mountain front fault. The deep seismic reflection section highlights the North China Basin’s structural features, characterized by alternating depressions and uplifts, with boundaries clearly delineated by faults such as the Cangxi, Cangdong, and Chengxi faults. These faults control the formation of graben and horst structures and are considered concealed active faults since the Quaternary period. The Cangxi fault, as the eastern boundary of the Jizhong depression, developed in the weak zone of the front thrust nappe tectonic area of the detachment slip structure. The fault plane resembles a plow shape, steep at the top and gently sloping at the bottom. It utilized or transformed the early thrust section, which is now manifested as a west-dip normal fault, controlling the basement structure and stratigraphic sedimentation on the west side of the Cangxian uplift. The Cangdong fault is the eastern boundary fault of the Cangxian uplift, which appears as an east-dipping shovel shaped normal fault on the profile, cut through the reflection waves of the Carboniferous-Permian strata, the Cambrian-Ordovician strata, the Proterozoic strata, and the crystalline basement. It terminates at the interface of the upper and lower crust at a depth of about 18km. The Chengxi fault is a west-dipping normal fault, which cuts through the Cenozoic sedimentary layer at a depth of about 600~700m in the shallow section. It terminates at the interface between the upper and lower crust in a shovel shaped normal fault downward. The deep seismic reflection section also clearly shows the coexisting structural morphology of uplift and depression. Multiple secondary faults that tilt in the same direction or opposite direction to the main fault are developed inside the depression, causing the depression to be divided into multiple secondary structural units, resulting in the complexity of the entire fault basin structure.

    In conclusion, the development of boundary faults plays a decisive role in the stratigraphic sedimentary and tectonic deformation of the strata within the depression.

    The existing deep and shallow structures and tectonic patterns in the Wuji-Yanshan section of the North China Basin are formed by the “graben-horst” structure developed in the upper crust, the complex fault combination style near the surface, the stratified reflection and the upper and lower superimposed reflection structure developed in the lower crust, and the undulating Moho surface. The findings of this study contribute to the seismological understanding of the dynamic processes occurring in the North China Basin, as well as to the analysis of the structural relationship between deep and shallow structures in the region.

    Table and Figures | Reference | Related Articles | Metrics
    REGIONAL VELOCITY STRUCTURE AND RELOCATION OF THE 2022 DONGTAI EARTHQUAKE SEQUENCE
    YU Yue-ying, LI Zheng-kai, YANG Yun, KANG Qing-qing, QIAN Jia-wei, WANG Jun-fei, QU Min, ZHOU Yu-chen, LI Ying-chun, XU Tian
    SEISMOLOGY AND GEOLOGY    2024, 46 (3): 627-644.   DOI: 10.3969/j.issn.0253-4967.2024.03.007
    Abstract199)   HTML11)    PDF(pc) (9248KB)(56)       Save

    The seismic activity in Dongtai, Jiangsu, suddenly intensified from November to December 2022. The largest earthquake observed during this period was a magnitude MS3.0 event on 25 December, which was felt reported by many nearby residents and caused a certain degree of social impact. The Dongtai area is situated in the central part of Jiangsu Province, within the Dongtai depression in structure, which is a secondary tectonic unit in the southern part of the Subei-South Yellow Sea Basin. Multiple fault zones developed in the region. The prominent known fault zones near the epicenter include the Taizhou fault, the Chenjiabao-Xiaohai Fault, and the Benchahe Fault. Among them, the closest to the epicentral area is the Taizhou fault. Additionally, the Subei Basin has a long history of industrial activity. Its geological conditions are complex, and the resources are extremely scattered and fragmented. The scale of underground resource extraction is predominantly small to medium-sized and has entered the middle to high exploration level. Historically, Dongtai has experienced weak seismic activity with only six earthquakes of MS≥3 within 50 kilometers of the epicenter since 1970. The sudden increase in seismic activity prompts investigation into its cause. Analyzing the structural features of the Dongtai earthquake sequence can enhance understanding of seismic activity and seismogenic mechanisms in the region.

    Previous studies on regional velocity structure have primarily focused on large scales, such as the Tan-Lu fault zone, with no specific research dedicated to the Dongtai earthquake sequence. In this study, we collected earthquake arrival time data recorded by the China Earthquake Networks Center from 2008 to 2022. Employing the double-difference tomography method, we conducted a joint inversion to investigate the velocity structure and earthquake locations in the Subei Basin. The resulting outcomes include the three-dimensional P-wave velocity structure of the area and the relocation results of 22 events within the seismic sequence. Furthermore, utilizing clear P-wave initial motion data from station waveform records, we inverted the focal mechanism solutions of the earthquake sequence using a modified grid search method. By integrating these inversion results with data on fault distribution and local industrial activity, we discussed the earthquake-triggering mechanism and possible seismogenic structures.

    The results indicate that: 1)Following relocation, the seismic sequence exhibits a zonal distribution pattern. The earthquakes are predominantly situated north of the Tai-Zhou fault in a nearly north-south orientation, spanning approximately 15 kilometers in total length, with a predominant depth range of 11 to 16 kilometers. Notably, there is no apparent correlation between the earthquakes and the surrounding known fault structures. 2)The focal mechanism solution parameters for the largest earthquake in the sequence, MS3.0, suggest a strike-slip seismogenic structure with a minor normal component. The direction of the stress axis aligns closely with the current tectonic stress field of the study area. Based on the focal mechanism solution and the distribution of the sequence, it is inferred that a dextral strike-slip hidden structure trending in a NNE-SSW direction may exist beneath the sequence. 3)The velocity structure of the epicenter area exhibits significant heterogeneity. The middle crust displays relatively high velocity, while the lower crust shows relatively low velocity. Notably, a spindle-shaped high-velocity anomaly with a P-wave velocity of 6.25km/s is observed at a depth of approximately 15km. The earthquakes primarily cluster southeast of this anomaly. 4)By examining the relationship between the spatial locations of earthquakes and their occurrence times, it is observed that the epicenters exhibit a seismogenic process extending far from the edge of the anomalous body. This suggests the outward release of accumulated elastic energy within the high-velocity anomaly, indicating a potential relationship between earthquake occurrences and the velocity anomaly. 5)Through on-site investigations of the epicentral area, data regarding local industrial activities have been collected. It was observed that three new wells and multiple industrial operation points have been established in the seismic area. Remarkably, 73% of earthquakes in the seismic sequence occurred within a 4.6km radius of well H1, with the largest earthquake in magnitude located approximately 1km from the well. A notable correspondence is observed between the Wulie-Shiyan-Qindong extraction points, the seismic sequence, and the deep high-velocity anomaly. Additionally, the operational timeframe of newly developed wells in the region closely aligns with the timing of earthquakes. However, the dominant depth of seismicity does not correspond with the drilling depth. A preliminary inference suggests that the occurrence of the earthquake sequence may be linked to the deep heterogeneous velocity structure, while industrial production activities near the epicenters may induce alterations in the regional stress state, leading to stress destabilization and subsequent energy release.

    Table and Figures | Reference | Related Articles | Metrics
    MECHANISM DIFFERENCES BETWEEN SEVERAL TYPICAL PYROCLASTIC ROCKS AND THEIR VOLCANISM SIGNIFICANCE
    WEI Hai-quan, CHEN Zheng-quan, LIU Yong-shun, BAI Zhi-da
    SEISMOLOGY AND GEOLOGY    2024, 46 (2): 297-311.   DOI: 10.3969/j.issn.0253-4967.2024.02.004
    Abstract198)   HTML13)    PDF(pc) (6205KB)(154)       Save

    Pyroclastic rock is the most direct object of physical volcanology and the most important topic of identifying the volcanic explosive fragmentation processes. Some particular species of pyroclastic rocks and equivalents can indicate key characteristics of the volcanism process, which is the basis to estimate the eruptive risks. Volcanic hazard is potential risk related to volcanic eruption, and it is one of the most important types of disasters that human beings face in nature. Volcanic disasters are directly related to the types of volcanic eruptions, among which explosive volcanic eruptions can cause the deadly intensive volcanic risks. The direct product of explosive volcanic eruption is to form various pyroclastic rocks, which represent the different types and intensities of volcanic hazards caused by the eruption process. The primary pyroclasts and secondary fragments reflect the difference of volcanic surface processes during eruptive or intermittent periods, while the distinguish of magmatic, phreatomagmatic and phreatic eruptive deposits marks the systematic development of modern volcanology, which is the leading work in the study of volcanic hazards. 1)Pyroclastic rocks are formed directly by transporting, accumulating and diagenesis of the expelled materials during the eruption. They usually consist of the primary materials such as broken magma, accidental fragments trapped by the volcanic conduit, as well as the epiclasts captured by the volcanic fluid flowing on the surface. Pyroclastic rock, as a direct product of explosive volcanism, has naturally becomes the most important research object in volcanology. The volcanic tephra laminae preserved by fine airfall volcanic ash in basins has been attracted attention because of their good isochron and environmental indication, and the associated rocks may need to be distinguished from different types of volcanic sedimentation such as bedded tuff, sedimentary tuff and tuffaceous mudstone. The autoclastic breccia produced by lava emplacement and the hyaloclatite formed by the quenching of lava under water represent fragmentation that is closely related to the lava flow, rather than those from explosive volcanism. 2)Pyroclast is mainly the product of explosive volcanism, but it can contain a certain amount of normal sedimentation and a small amount of rock fragment near the volcanic channel and the magma chamber roof. Pyroclats are generally defined as the direct products of explosive eruption behavior, while volcaniclastics are formed by volcanic degradation such as slope displacement, avalanche, lahar, and the autoclast generated by lava flowage and quenching. This classification not only emphasizes the difference in the forming process of different volcanic products, but also helps to distinguish the different mechanism in volcanological research and hazard estimation. Different types of pyroclastic rocks are formed with different fragment mechanisms and diagenetic ways, and some specific pyroclastic rocks represent various special types and scales of volcanic hazards. Although they are usually classified as primary clastics, the hazard caused by autoclastic breccia is significantly different. Cryptoexplosive breccia, although we have employed a rock name from pyroclastic rocks, is actually more concerned with its resource economics. 3)When we study the genetic types of pyroclastic rocks, the most important basis for identification is the forming mode of the materials, that is, the type of fragmentation, which include primary volcanism and secondary volcanism. Primary clasts are divided into pyroclast, which is formed by the direct action of volcanic eruption, and autoclast, which is produced by the flow process of lava flows, While secondary(exogenous)volcanism includes various kinds of exogenous clasts(epiclast)formed by volcanic surface processes. According to the proportion of magma and water content at eruptive environment, explosive eruption can be divided into three types: magmatic eruption, phreatomagmatic eruption and phreatic eruption, which represent the most basic process of explosive eruption, and are also the problems of genetic classification and identification often faced in the study of pyroclastic rocks.

    Table and Figures | Reference | Related Articles | Metrics
    CRUSTAL FINE VELOCITY STRUCTURE IN THE ERYUAN AREA, YUNNAN FROM DOUBLE-DIFFERENT TOMOGRAPHY
    CAO Ying, QIAN Jia-wei, HUANG Jiang-pei, ZHOU Qing-yun
    SEISMOLOGY AND GEOLOGY    2024, 46 (1): 162-183.   DOI: 10.3969/j.issn.0253-4967.2024.01.010
    Abstract197)   HTML19)    PDF(pc) (14442KB)(144)       Save

    The Eryuan area is located in the central part of the northwest Yunnan region on the southeastern edge of the Qinghai-Xizang Plateau. The geological structure in the area is complex, including the Weixi-Qiaohou-Weishan Fault and Honghe Fault, as well as the Longpan-Qiaohou Fault and Heqing-Eryuan Fault, which intersect in an “X” shape. The geothermal activity in the area is also active and exhibits strong fault structural characteristics. Earthquake activity is also frequent in the Eryuan area. Since 2013, multiple earthquakes with MS≥5.0 have occurred on the west side of the Weixi-Qiaohou-Weishan Fault. Given the complex geological and geothermal background and seismic activity in the Eryuan area, we use the seismic travel time data recorded by the Yunnan Regional Seismic Network and the Northwest Yunnan Dense Array from January 1, 2008, to July 20, 2023, VP/VS model consistency-constrained DD tomography method to obtain the 3D VP, VS, and VP/VS models and relocation results in the Eryuan and its surrounding areas. The research results indicate that: 1)based on the relocation results, the dense small seismic cluster develops at the intersection of the Weixi-Qiaohou-Weishan Fault, Honghe Fault, and Heqing-eryuan Fault deserves special attention. Four earthquakes with MS≥5.0 that have occurred since 2013 are mainly distributed on the west side of the Weixi-Qiaohou-Weishan Fault, with an NNW-SSE trend distribution. The fault system at the western boundary of the Sichuan-Yunnan block is very complex, and its seismic risk deserves our special attention. 2)Based on the tomography results, the widely distributed low-velocity anomalies in the upper crust of the study area may be related to the crustal material migration pathway caused by the escape of the Sichuan-Yunnan diamond block towards SSE. 3)Combining imaging results with geochemical research results for inference, the high VP/VS below 10km depth beneath the small seismic cluster at the intersection of the Weixi-Qiaohou-Weishan Fault, Honghe Fault, and Heqing-Eryuan Fault may correspond to hot spring geothermal fluids. Moreover, the circulation process of hot spring fluids in complex fault systems may have penetrated to a depth of 7~10km below the seismic cluster, accompanied by some small earthquakes. However, there is no imaging evidence of fluid within a depth of 5~7km where some earthquakes are densely distributed. The occurrence of these earthquakes may be related to the distribution of brittle rocks, but it cannot be ruled out that fluid may have a potential infiltration effect on them in the future. Based on the results of velocity structure, they are combined with pre-existing seismology research results. It is found that the Eryuan earthquake sequences on March 3 and April 17, 2013 were mainly located within low VP, low VS, and low VP/VS anomaly. Generally, if there is fluid present, VS decreases faster than VP, resulting in high VP/VS. So low VP/VS indicates that low VS is not caused by fluids, which may be due to lithology. Therefore, the imaging evidence indicates that there is no fluid in the region where the sequence is located, thus inferring that the occurrence of the sequence is not directly related to the fluid. The MS5.1 Yangbi earthquake sequence on March 27, 2017, is spatially close to the Eryuan sequence and also has the same velocity structure characteristics, so it may not be directly related to fluids. The mainshock and some aftershocks of the MS5.1 Yunlong earthquake sequence on May 18, 2016, were mainly located in high VP, high VS, and relatively high VP/VS regions. High VS indicates that it is not the fluid that causes relatively high VP/VS. It is speculated that there may not be fluid in the region where the mainshock is located, and the fluid did not directly participate in the occurrence process of the sequence.

    Table and Figures | Reference | Related Articles | Metrics
    STUDY ON THE SEISMOGENIC STRUCTURE OF THE 2022 GUJIAO ML4.1 EARTHQUAKE IN SHANXI PROVINCE BASED ON FOCAL MECHANISM AND SEISMIC LOCATION
    DONG Chun-li, ZHANG Guang-wei, LI Xin-wei, WANG Yue-jie, DING Da-ye, GONG Zhuo-hong
    SEISMOLOGY AND GEOLOGY    2024, 46 (2): 414-432.   DOI: 10.3969/j.issn.0253-4967.2024.02.010
    Abstract195)   HTML14)    PDF(pc) (4880KB)(204)       Save

    Understanding the mechanism of earthquake sequence in the mining area is important for the time-dependent hazard assessment. An earthquake of ML4.1 occurred in Gujiao, Taiyuan, Shanxi on February 20th, 2022, which caused strong ground motion in Gujiao and surrounding counties. The epicenter of this earthquake is located in the area of Lvliang uplift, where historical earthquakes are relatively rare. In addition, the coal resources are well developed in the earthquake source area which has attracted much attention from society and local governments.

    To investigate the mechanism and the seismogenic fault of Gujiao ML4.1 earthquake, we first apply the double-difference location method to retrieve highly accurate hypocenter locations. The results show that the earthquakes mainly occur at a depth range of 3~5km, and display a dominant distribution direction nearly EW-trending, which differs significantly from the NE-trending fault distribution pattern in this region. We further collect the broad-band seismic waveforms from the regional network of Shanxi province to perform focal mechanism inversion. The inversion results show that the Gujiao earthquake is a left-slip seismic event with a moment magnitude of MW3.96. The optimal double-couple solution is characterized by a strike of 90°, dip of 80°, and a rake angle of -21° for fault plane Ⅰ, while for the fault plane Ⅱ, the strike is 184°, dip is 69°, and rake angle is -169°. The best centroid depth is estimated to be at 3km. This earthquake shows an extremely shallow focal depth. Moreover, By using cluster analysis method, we obtained the central solution for the seismogenic fault plane of the GuJiao earthquake, with a fault strike of 91°and a dip angle of 70°. The focal solutions show that the earthquake exhibit a strike-slip type, and the orientations of earthquake sequence coincide well with the focal mechanisms.

    In addition, to discuss the effect of Gujiao ML4.1 earthquake on regional stress, we calculate the stress drop of this seismic sequence. The results show that the stress drop is significantly smaller than that of the regional earthquakes, exhibiting at least one order of magnitude lower than that of the background earthquakes in the same region. This phenomenon reflects that the stress level in the focal area of the GuJiao earthquake is not high, suggesting that the background stress enhancement in the focal area is not obvious.

    Based on regional geological structure, we found that the known faults in the region are all high-angle normal faults, and the strike of these faults are inconsistent with the focal mechanism solution of Gujiao earthquake sequence, which suggests that the existing faults are not the seismogenic fault. Taking the regional mining activities into account, we speculated that mining may cause strong disturbance to the stress field, and lead to stress redistribution within the rock mass. Such coal mining activity may generate a high stress disturbance on the hidden fault plane, and then the fault become the carrier of stress transfer. So we conclude that the seismogenic mechanism of the Gujiao-seismic sequence may be related to coal mining activities near the focal area, which leads to local stress changes, thus resulting in the activation of preexisting hidden faults and triggering the occurrence of the Gujiao earthquake.

    Table and Figures | Reference | Related Articles | Metrics
    A WEB-BASED PROTOTYPE SYSTEM FOR THE THREE-DIMENSIONAL FAULT MODELS OF THE CHINA SEISMIC EXPERIMENTAL SITE
    WU Xi-yan, LU Ren-qi, ZHANG Jin-yu, SUN Xiao, XU Fang, CHEN Gui-hua
    SEISMOLOGY AND GEOLOGY    2024, 46 (1): 35-47.   DOI: 10.3969/j.issn.0253-4967.2024.01.003
    Abstract189)   HTML20)    PDF(pc) (6764KB)(153)       Save

    The data of active fault structure and three-dimensional(3D)fault models is essential for seismic risk analysis. With more and more requirement for complex 3D fault models, the demand for data sharing and related research increases dramatically. A web-based display system for three-dimensional fault models would improve data sharing and user experience. Moreover, constructing such a web-based system is also an important issue for data sharing.

    The 3D active fault models are built in a data modeling platform, while the web display system is constructed by the geographic information system(GIS)platform. Because the data structure, type, and content between data modeling and GIS platforms are different, the following questions are critical, for example, how to migrate 3D model data from the modeling platform to the GIS platform?and can the migrated data present the right attributions?In this paper we used the Web AppBuilder of ArcGIS 10.6 Enterprise Edition to build a Web prototype system to display 3D fault models of the China Earthquake Science Experimental Field(Sichuan-Yunnan region). The system implemented the basic functions of a 3D Web application and successfully tested the 3D scene display scheme, user interaction mode, and data migration scheme.

    The prototype system adopted a local scene, which can easily switch between the above-ground and underground viewing angles of the scene. The scene included 2D fault surface traces, 3D fault models, and earthquakes with or without focal depth. After data fusion, the 3D fault models were classified and displayed with active age, having a good visual fusion effect with 2D fault data. Earthquakes with or without focal depth were displayed in different colors. The earthquakes without focal depth were uniformly displayed at 17km depth according to the average focal depth of the earthquakes with focal depth. So the earthquakes without focal depth can be highly consistent with other elements in the 3D scene.

    The user interface interaction mode in the 3D scene of the prototype system was consistent with the common interaction mode of 2D map applications in the following aspects: 1)map browsing; 2)Navigation menu; 3)Geographical inquiry; and 4)Functional interactive tools. The system interface was simple, clear, logical, and unified. Users were easily acquainted with the three-dimensional scene interface according to the two-dimensional map interaction experience. It conformed to the user interface interaction principles of simple, consistent, predictable, and easy feedback.

    The prototype system had the basic functions of 3D scene browsing, zooming in and out, 3D object attribute viewing, geographic query, base map switching, layer control, legend, and distance measurement. However, the prototype system needed further development and more complex functions such as data attribute table browsing, space selection, and space query.

    This paper presented a data migration scheme from the modeling platform to the GIS platform. The data migration of this scheme can be divided into four steps: data format conversion, coordinate system conversion, 2D and 3D attribute information mapping, and 3D data attribute table construction. After transforming the data format and coordination system from the modeling platform to the GIS platform, 2D and 3D data fusion should be carried out to make 3D data and 2D data have the same attribution. The format conversion and coordinate system conversion steps can be automatically completed in batches. Otherwise, mapping the 2D and 3D attribute information and building the 3D data attribute table need manual handling.

    In summary, this paper presents a data migration scheme from the modeling platform to the GIS platform. Practice in reality shows that only after conversing data format and coordination system from the modeling platform, the 2D and 3D data fusion steps are caplable of ensuring a better visual integration of them. The Web-based prototype system of displaying 3D fault models of the China Seismic Experimental Site implements the basic functions of 3D scene application and tests the fused 2D and 3D data visualization. It is friendly and open to users, with a great demonstration significance.

    Table and Figures | Reference | Related Articles | Metrics
    GENESIS OF SONGYUAN EARTHQUAKES BASED ON 3D RESIDUAL DENSITY STRUCTURE
    LIU Wen-yu, CHENG Zheng-pu, NIAN Xiu-qing, CHEN Yan, HU Yu-ling, QIN Zu-jian, SHAO Ming-zheng
    SEISMOLOGY AND GEOLOGY    2024, 46 (2): 462-476.   DOI: 10.3969/j.issn.0253-4967.2024.02.013
    Abstract185)   HTML15)    PDF(pc) (7342KB)(110)       Save

    Since 2006, more than 10 earthquakes with the magnitude of 4 or higher have occurred continuously in Ningjiang district and Qian’an county, Songyuan city of Jilin province, with tens of thousands of aftershocks. The largest earthquake is the magnitude 5.8 earthquake on November 23, 2013 in Qiangorlos Mongolian autonomous county, Songyuan city, Jilin province. Frequent earthquake activities not only cause a large number of damaged buildings, but also trigger geological disasters such as sand liquefaction and slope instability, which have attracted widespread attention from the public and government departments. In order to comprehensively understand the causes of earthquakes in the Songyuan area, much studies on earthquake anomalies, seismogenic structures, and source depths have benn conducted recently. However, there are still some difference opinions in understanding the causes of earthquakes in the Songyuan area, for example, the major different viewpoints include: 1)it is caused by the regional stress release caused by Pacific subduction, mainly based on the consistency between the stress field in the source area and the background stress field of the Pacific plate subducting towards the edge of the east Asian continent; 2)it is related to deep melt or fluid migration, mainly based on the existence of low-speed and low resistance anomalous bodies in the deep part of the seismic area; 3)it is related to long-term oil and gas extraction, the main basis is that the seismic source spectrum exhibits characteristics of early and fast attenuation compared to typical structural earthquakes, and the time-domain and frequency-domain characteristic parameters such as waveform complexity and spectral ratio are also significantly different from typical structural earthquakes. In addition, the seismic source mechanism contains a large number of nondual force source components and shallow seismic source depth. Moreover, there is also significant controversy over the seismogenic structure, with the focus on whether the Qian’an and Ningjiang earthquake regions are controlled by the same fault zones. Some scholars infer that the seismogenic structures in both areas are the NE-trending Fuyu-Zhaodong Fault based on the source mechanism solution. Another group of scholars believe that the Qian’an earthquake area is controlled by NW-trending hidden faults, while the seismic structure in the Ningjiang earthquake area is the second NW-trending Songhuajiang Fault. The reason for the above controversy lies in the lack of necessary constraints on the deep structure of the seismic area.

    To solve above problems, in this article we studied the three-dimensional residual density structure of the Songyuan earthquake area by performing nonlinear conjugate gradient focusing inversion on regional gravity data. Combined with oil drilling and reflection seismic data, some new insights were obtained as follows: 1)The residual density anomaly in the study area shows a high-low alternating strip distribution, with the southern anomaly trending NNW and the northern anomaly transitioning to NNE, This feature reflects the different deep earthquake environments in the Qian’an earthquake area and the Ningjiang earthquake area. The former’s source is located in the high-density anomaly body and its edge of Chaganhua, while the latter’s source is located in the middle of the low-density anomaly zone in Songyuan, indicating different rock types in the two locations. The basement of the Qian’an earthquake area is composed of limestone and metamorphic volcanic rocks, while the basement of the Ningjiang earthquake area is composed of fractured granite magmatic rocks. The stability of the basement structures in both areas is poor; 2)The seismic structures of the two are different. The former is controlled by the Chaganhua Fault and Qian’an Fault, while the latter is controlled by the Songyuan Fault and the Fuyubei Fault; 3)The formation of earthquakes is related to factors such as regional stress, basement structure, deep gas/fluid migration, and long-term oil and gas extraction. Long term water injection and oil recovery have damaged geological structures and stress environments, which may be important triggering factors. 4)In terms of the seismic source mechanism, we concluded that under the sustained action of regional stress in the near EW direction in the Songyuan area, the basement fault structures and fractures in the Qian’an and Ningjiang earthquake areas are extremely developed. Deep gas/fluid migrates upwards along the basement fractures, further reducing the stability of the basement structure. Long term water injection, oil and gas recovery, as well as fracturing, greatly damage the structural integrity of the cap rock and the distribution of deep stress, which results in the compressive and torsional strike slip movements in the relevant structures.

    Table and Figures | Reference | Related Articles | Metrics
    SIMULATION OF THE ROCK SURFACE LUMINESCENCE SIGNALS ON BEDROCK FAULT SCARPS BY STICK-SLIP AND CREEP MOVEMENTS
    LUO Ming, CHEN Jie, QIN Jin-tang, YIN Jin-hui, YANG Hui-li, LIU Jin-feng, GONG Zhi-jun
    SEISMOLOGY AND GEOLOGY    2024, 46 (2): 357-370.   DOI: 10.3969/j.issn.0253-4967.2024.02.007
    Abstract184)   HTML9)    PDF(pc) (3308KB)(130)       Save

    The reconstruct of the stick-slip and creep histories is essential for understanding fault activities and seismic hazard assessment. Distinguishing stick-slip and creep using geodetic technology has become a hot research area in recent years, but distinguishing and estimating seismic slip and creep on geological timescales(e.g., over hundreds of years)is challenging due to the lack of historical, geodetic and remote sensing data extending back more than a few hundred years. This study uses a newly developed dating technique(rock surface optically-stimulated-luminescence(OSL)dating)combined with the OSL decay parameters of granite samples from the Langshan fault in Inner Mongolia to simulate optically stimulated OSL-depth curves and depths of half saturation of luminescence signal under various scenarios such as fault seismic slipping, creeping, and erosion of colluvial wedge. The study compares these OSL-depth profiles, especially the depths of the half saturation, under different slipping modes, and summarizes their features.

    During fault seismic slip, samples at different heights along the fault scarp display a “step-like” distribution pattern at their depths of half saturation. While during creep, however, they exhibit a “slope-like” pattern. Such differences may lie in that the slope during accelerating creeping is steeper than the slope during constant-speed creeping. Correspondingly, the resolution of residual luminescence-depth profile and depth of half saturation is also higher during accelerating creeping. During intra-earthquake creep events between seismic slip occurrences on the bedrock fault scarp, the distribution of half-saturation depth in the samples includes segments resembling both “steps” and “slopes”, which indicate the seismic slip and creep activities of the fault respectively. If the samples at the base of the colluvial wedge have had a sufficiently long last exposure time, the luminescence-depth profile and half-saturation depth distribution due to the erosion of the colluvial wedge would be approximately the same as in the three-phase seismic slip scenario. This indicates that samples previously buried by the colluvial wedge may be considered within the seismic displacement. Conversely, if the last exposure time of the base samples at the base of the colluvial wedge is short, the bleaching depth of the luminescence signal of these base samples will be noticeably shallower than that of the other samples within the seismic displacement, indicating the observed erosion of the colluvial wedge in this case. Furthermore, the seismic displacement ideally should include the buried location of the colluvial wedge. Therefore, when the luminescence curves and half-saturation depth distributions fail to identify the presence of the colluvial wedge, it is acceptable to include the buried location of the colluvial wedge in the seismic displacement calculation. Conversely, the luminescence-depth curves and half-saturation depth distributions document the erosion caused by the colluvial wedge. The simulation results demonstrate that this method can effectively distinguish between fault slipping and creeping, obtain corresponding displacements, and potentially record the erosion of colluvial wedge.

    This study also analyzes the temporal resolution of the method for distinguishing fault activity times and the spatial resolution for quantifying displacements. The specific situation is as follows. When exposure age of the bedrock fault scarp is within a thousand years, the rock surface OSL dating method can easily distinguish types of active slips and seismic displacements for the earthquakes with a recurrence interval of hundreds of years. When exposure age of the bedrock fault scarp is in the range of 100-101ka, the method can easily distinguish types of active slips and seismic displacements for the earthquakes with a recurrence interval exceeding a thousand years. When exposure age of the bedrock fault scarp is over ten-thousand years, the resolution of this method may be significantly reduced. The spatial resolution of seismic displacements using this method depends on interval between sampling and testing samples, typically in 10~30cm.

    Table and Figures | Reference | Related Articles | Metrics
    FULL-WAVEFORM INVERSION BASED ON EXPONENTIAL-PHASE COHERENCY MISFIT FUNCTION
    LIU Jian-huan, CHEN Jian-ye, ZHAO Ji-hai, Deyan Draganov
    SEISMOLOGY AND GEOLOGY    2024, 46 (1): 48-62.   DOI: 10.3969/j.issn.0253-4967.2024.01.004
    Abstract183)   HTML7)    PDF(pc) (3494KB)(75)       Save

    Full waveform inversion(FWI)has emerged as a highly effective approach for obtaining accurate and high-resolution S-wave velocity structure of the shallow subsurface. However, there exist several challenges when applying FWI in the field. The first challenge is the issue of local minima that arises when inaccurate initial models are used, especially when inverting shallow surface seismic data dominated by high-frequency and strongly dispersive surface waves. These local minima are caused by the non-linear misfit function that represents the differences between the measured and simulated data. Moreover, defining an appropriate minimization criterion to reduce the sensitivity of the inversion results to errors in the recorded seismic wave amplitudes is another significant challenge. Amplitude errors can arise from various factors, such as different coupling effects between seismic sources, receivers, and the ground, or variations in the strength of seismic sources excited at different shot locations. If the amplitude information of the recorded seismic wavefield is inaccurate, the reliability of the FWI inversion results will be effected negatively.

    To overcome these problems, we propose a novel misfit function that incorporates exponential-phase coherency, thereby eliminating the reliance on amplitude information from the measured and simulated data. This new misfit function is designed to measure the coherency between the measured and simulated data based on exponential phase. It achieves a balance in extracting valuable information from various amplitude components of the recorded seismic wavefield, such as surface waves, reflections, and scattering waves. The method for computing this coherency is inspired by Phase-Weighted Stacking(PWS)method to detect weak but correlated seismic signals. In PWS, the exponential phase of the data is computed to obtain a phase-dependent coherency that is independent of amplitude. This correlation is then used to enhance the stacking of signals with similar instantaneous phases.

    By utilizing the adjoint-state method, we efficiently calculate the gradient of the misfit function with respect to the model and conduct a thorough analysis of its shape and characteristics. To demonstrate the robustness of our proposed misfit function against random noise, we perform experiments by using simulated data contaminated with varying levels of noise. The results demonstrate that the misfit function based on exponential-phase coherency remains highly robust and reliable, even in the presence of significant random noise. This robustness is particularly crucial in practical applications where noise contamination is a common challenge.

    To evaluate the performance of FWI employing exponential-phase coherency in a real field environment, we employ seismic data collected at an archaeological site as a benchmark dataset. This dataset presents a complex and challenging scenario due to the presence of complex subsurface structures. In addition to the inherent noise in the data, we introduce additional random noise to assess the robustness of our proposed misfit function.

    The inversion results obtained with our novel FWI approach exhibit an impressive resemblance to the known velocity structure of the archaeological site. These results are further validated through independent archaeological excavations, which confirms the accuracy and reliability of our imaging technique. The success of using our method in accurately reconstructing subsurface features under challenging field conditions underscores the significant potential of FWI based on exponential-phase coherency to enhance the accuracy and reliability of shallow subsurface imaging in practical scenarios.

    Table and Figures | Reference | Related Articles | Metrics
    STUDY ON THE RELATIONSHIP BETWEEN LITHOSPHERIC MAGNETIC FIELD AND GEOLOGICAL STRUCTURE AND SEISMIC ACTIVITY: TAKING THE 2021 MS6.4 YANGBI EARTHQUAKE AS AN EXAMPLE
    CHEN Zheng-yu, NI Zhe, ZHOU Si-yuan, JIN Yun-hua, YANG Xin-jun
    SEISMOLOGY AND GEOLOGY    2024, 46 (2): 449-461.   DOI: 10.3969/j.issn.0253-4967.2024.02.012
    Abstract182)   HTML13)    PDF(pc) (5142KB)(143)       Save

    The lithospheric magnetic field is an important part of the earth’s magnetic field, which is affected by many factors, such as rock’s magnetization environment, underground geological structure, rock mineralogical composition, thermal and pressure state, and the deep tectonic evolution process. Most earthquakes occur in the crust and uppermost mantle, known as the lithosphere. The preparation and occurrence of earthquakes are usually accompanied by long-term accumulation and sudden release of energy, which will lead to changes in the thermal and pressure state of rocks, causing magnetic susceptibility variation in the lithosphere. Previous studies show that seismic activities can cause abnormal changes in the geomagnetic field, and there is an obvious correlation between the lithospheric magnetic field and seismic activities. The MS6.4 Yangbi earthquake on May 21, 2021, provided a unique opportunity to study the dynamic evolution of seismo-magnetic anomaly.

    Five-term repeat station vector geomagnetic data observed in Yangbi and surrounding areas from 2018-2021 were used in this paper, the first four terms were observed before the earthquake, and the fifth term was observed after the earthquake. After data processing and model calculation, the lithospheric magnetic fields before and after the earthquake are obtained, lithospheric magnetic field models are established using the Surface Spline method, and annual variations are calculated. Based on the analysis of lithospheric magnetic field combined with the regional geological structure, the Northwest Sichuan Subblock shows negative magnetic anomalies, which is consistent with the geological structural characteristics in the study area, altitude and crustal thickness increase sharply from Central Yunnan Subblock to Western Sichuan Plateau. Small areas of positive and negative magnetic anomalies are alternatively distributed in the Central Yunnan Subblock, which reflects the heterogeneity of deep lithosphere structure. The negative magnetic anomaly in the western boundary of the study area is also consistent with the geological characteristics of the Qingzang Plateau. There is also a correspondence between lithospheric magnetic field anomalies and faults, especially the total intensity. Negative magnetic anomaly strips are distributed along the strike of the Honghe Fault and Lijiang-Xiaojinhe Fault, while Weixi-Qiaohou-Weishan Fault appears at the junction of positive and negative magnetic anomalies. The statistical analysis of the MS6.0 and above seismic events from 1970 to 2021 shows that there is a correlation between lithospheric magnetic anomalies and seismic activities. Most earthquakes occur in the weak magnetic anomaly area, especially near zero contour. The earthquakes tend to be distributed in anomaly gradient belts, and the number of earthquakes in the negative anomaly area is higher than that in the positive anomaly area. Analyzing the characteristics of the pre- and post-seismic changes of declination and total intensity near the epicenter of Yangbi MS6.4, it is found that the epicenter of the Yangbi earthquake is located near the zero-contour-line of declination. During the preparation of the Yangbi earthquake, the total intensity gradually changed from a balanced distribution of positive and negative anomalies to the overall negative changes, and the magnetic anomalies recovered the trend of the balanced distribution of positive and negative changes after the earthquake.

    Table and Figures | Reference | Related Articles | Metrics
    THE NEW FINDINGS OF SURFACE RUPTURE ZONES AND ITS SEISMOLOGICAL SIGNIFICANCE OF THE EASTERN MARGIN OF YUMUSHAN FAULT, NORTHEASTERN MARGIN OF QINGZANG PLATEAU
    CHEN Bai-xu, YU Zhong-yuan, XIAO Peng, DAI Xun-ye, ZHANG Shi-long, ZHENG Rong-ying
    SEISMOLOGY AND GEOLOGY    2024, 46 (3): 589-607.   DOI: 10.3969/j.issn.0253-4967.2024.03.005
    Abstract181)   HTML29)    PDF(pc) (20130KB)(157)       Save

    The Hexi Corridor in northwest China has obvious structural deformation and complex fracture image. With the development of several NW thrust fault zones accompanied by a large number of ancient earthquakes and historical seismic events, the earthquake disaster is relatively serious. The eastern margin of Yumushan fault is one of them. The fault is mainly developed in the east site of Yumushan Mountain, with the exposed fault plane striking NW330° and dipping about 41°~85° to the southwest as a whole. Previous research data show that the Eastern Margin of the Yumushan fault is an important part of the Qilian Mountain active thrust fault system in the northeast margin of the Tibet Plateau. It also constitutes the boundary structure between Hexi Corridor and Yumushan Uplift. Its late Quaternary tectonic deformation and recent activity characteristics reflect the northward extension process of the Qilian Mountains and the remote collision effect of the northward extrusion of the Indian Plate. However, there are still some controversies in the study of the latest activity age and deformation characteristics of the eastern Margin of the Yumushan fault zone, which directly affect the seismic risk assessment along the fault line and the Hexi Corridor, as well as the in-depth understanding of the active structural characteristics of the northeast margin of the Tibetan plateau.

    Combined with remote sensing image interpretation, paleoseismologic excavations, aerial photogrammetry of unmanned aerial vehicles and late Quaternary dating, this study carried out field investigations and newly discovered the surface rupture zone of The Eastern Margin of Yumushan Fault and its activity characteristics. The results show that The Eastern Margin of the Yumushan Fault strikes NW330° combined with obvious thrust movement, which is manifested as a fault scarp landform. That’s revealing than the kinematics property of The eastern margin of Yumushan fault is dominated by thrust. The fault forms the dividing line between the Yumushan uplift and Zhangye Basin, and also the dividing line between pre-quaternary strata and Quaternary strata. The southwest side of the fault is dominated by pre-quaternary bedrock which constitutes a mountain landform. Late Quaternary sediments are exposed on the northeast side, and the Holocene strata are widely distributed around the Heihe River. The results show that there are obvious differences in the activity habits of the faults. With the Heihe River as the boundary, the fault activity difference is obvious on the south and north sides of the Heihe River. The latest surface fracture zone in the late Holocene was found along the Heiheokou segment(F1-1). And the Hongshahesegment(F1-2)showed pre-quaternary fault. It can be seen that the Miocene fine sandstone is in fault contact with the early Pleistocene glutenite and late Ordovician metamorphic andesite, and the fault gouge develops near the fault, which is gray-green and yellow-green with moderate hardness and easy to be wet when encountering water.

    The Heihekou segment(F1-1)starts from Daciyaohe River in the north, passes Xiaociyaokou, and reaches Heihekou in the south. The fracture zone moves towards NW330° and tends to SW, with a length of about 10km and a width of 3~10m. For river terraces, gullies, and platforms with young surface faults, the maximum height of the surface scarp is based on the DEM data generated by UAVs. The height of the T1 terrace fault scarp measured by two profile lines is(1.7±0.1)m to(3.3±0.2)m. In the excavating trenches, obvious evidence of fault activity such as traction bending of strata and directional arrangement of gravel can be seen. The strata consist of late Quaternary alluvial sand, gravel layer, loess layer, and silty layer. The optically stimulated luminescence dating results show that the latest surface rupture event occurred at(0.6±0.07)ka BP.

    According to the empirical formula between maximum vertical displacement(Dmax)and magnitude(M), the magnitude of the latest seismic event is estimated. The magnitude and potential seismic risk of the latest rupture event are evaluated. The results reveal that the maximum vertical displacement of the latest surface rupture event is(3.3±0.2)m. Based on the empirical relationship between magnitude and vertical displacement, it is concluded that a large earthquake rupture occurred in the eastern margin of the Yumushan fault in the late Holocene and the corresponding magnitude is estimated to be M7.5.

    Derived from the analysis of existing data, the fault in the eastern margin of the Yumushan fault may conform to the quasiperiodic earthquake recurrence behavior. And the recurrence interval of strong earthquakes may exceed 1 600a. The time interval between the latest event revealed in this paper and its last seismic event is about 1 800a, which is consistent with the time interval under the fault quasiperiodic earthquake recurrence model.

    The results show that the eastern margin of the Yumushan fault has intensive tectonic deformation in the late Quaternary and a large seismic background of M7 or above. The current kinematic mode of the fault is compressive shortening. Its geodynamic process may be mainly controlled by the northward extension of the Qilian Mountains and the remote collision effect of the northward extrusion of the Indian Plate. The deformation process of the fault may be in line with pre-spreading imbricate thrust deformation and the latest deformation has gradually extended from the basin-mountain boundary to the interior of Zhangye Basin, which provides new data to support the seismic risk assessment of the interior of the basin. At the same time, the latest deformation achievement of the eastern margin of the Yumushan fault has important scientific significance for improving the active tectonic image of the northeastern margin of the Qingzang plateau and discussing the kinematics model of the Qingzang plateau.

    Table and Figures | Reference | Related Articles | Metrics
    GEOCHEMICAL CHARACTERISTICS AND GENESIS OF SOIL GAS IN THE PINGYUAN M5.5 EARTHQUAKE
    SU Shu-juan, CHEN Qi-feng, SUN Hao, LIU Jun, FENG Liang-le, XU Ji-long, YANG Yan-ming, LUO Kun-li
    SEISMOLOGY AND GEOLOGY    2024, 46 (2): 433-448.   DOI: 10.3969/j.issn.0253-4967.2024.02.011
    Abstract180)   HTML19)    PDF(pc) (7314KB)(151)       Save

    At 2:33 am on August 6, 2023, a M5.5 earthquake occurred in Pingyuan county, Dezhou city, Shandong Province. The faults within the epicenter and adjacent areas are deeply buried by the thick Quaternary sediment cover on which human activity is intensive, which makes it difficult to determine the location of the buried active faults from the surface based on geological and geomorphological evidences. It is necessary to detect the location of the buried active faults around earthquake areas and estimate their seismic risk.

    In this study, based on the epicenter distribution direction of major earthquake and aftershocks, seismic and geological data of earthquake areas, and damage degree of local buildings, 4 survey lines with a length of 30km were arranged across the epicenters and adjacent areas, and the concentrations of Rn, CO2 and Hg in soil gas were measured on site, and the results are as follows:

    (1)There are obvious spatial differences in the concentrations of soil gas near the epicenter and its vicinities within the distance of 30km. Gas concentrations are relatively high near the epicenter areas and the east and west ends of 4 arranged survey lines, in contrast to those which are relatively low in other non-structural control regions. The spatial distribution pattern of Rn concentration in soil gas is basically consistent with that of CO2, which may be due to CO2 used as a carrier gas of Rn to migrate to the surface. At the southern end of the Lingxian-Guanxian Fault(F1), the spatial concentration patterns of Rn and CO2 gases exhibit multiple peaks or wide anomalous zones. It is speculated that the deformation zone of the fault rupture at this location is relatively wide, and there may be secondary permeable fracture zones in the west of the F1. The escape form of Rn and CO2 gas indicates that there may indeed be multiple small fault branches near the F1, and the fault structure is relatively complex.

    (2)The spatial concentration distributions of Hg, Rn and CO2 in the epicenter areas are similar to that in its eastern region. However, in the western region of the epicenter areas, the spatial concentration distributions of Hg, Rn and CO2 vary greatly, and the Rn and CO2 concentrations near the Jiucheng Fault(F3) in the west of the epicenter regions are higher than those near epicenters. It is speculated that this phenomenon may be related to the high-concentration gas migration caused by strong seismic tectonic activities and the special nasal geological structure controlled by F3.

    (3)The concentrations of Rn, CO2 and Hg in the soil show high-value anomaly zones near the F1 and F3, and the concentrations of Rn and CO2 in the west of F3 exceed those in the epicenter area. After further earthquake relocation analysis, the spatial distribution of aftershocks exhibit a trend from F1 to F3. Combined with geochemical and geophysical research results, it is inferred that Pingyuan M5.5 earthquake should be related to the deep tectonic activities of F1 and F3.

    Above research results show that the soil gas geochemical method can be applied to define the location and distribution direction of the buried faults with thick overburden, which provides an important criterion for earthquake trend tracking analysis. This study is of greatly scientific significance in determining the dynamic source and genetic mechanism of Pingyuan M5.5 earthquake, identifying potential strong earthquake hazard areas, and assessing the risk of future earthquakes in the study area.

    Table and Figures | Reference | Related Articles | Metrics
    EVIDENCE FOR THE HOLOCENE ACTIVITY OF THE LEIBO FAULT ZONE
    ZHANG Guo-xia, SUN Hao-yue, LI Wei, SUN Wen
    SEISMOLOGY AND GEOLOGY    2024, 46 (1): 141-161.   DOI: 10.3969/j.issn.0253-4967.2024.01.009
    Abstract177)   HTML34)    PDF(pc) (18061KB)(192)       Save

    The Yingjing-Mabian-Yanjin tectonic zone(YMYTZ)is an important boundary structure between the southeastern margin of the Tibet Plateau and the Sichuan Basin. It consists of several small-scale secondary faults with different strikes and is generally characterized by the intersections of north-northwest oriented longitudinal faults and nearly east-west oriented transverse faults. The YMYTZ is seismically very active in the late Quaternary and hosted several moderate-strong earthquakes, including two M≥7 earthquakes since 1216AD, namely the 1216 Mahu earthquake and the 1974 Daguanbei earthquake. After the Daguanbei earthquake, several M≥6 earthquakes and hundreds of M≥5 earthquakes occurred along the YMYTZ to date, implying it is a newly generated seismotectonic belt. Even so, the activity of each fault is still unclear, bringing out great uncertainty in understanding the current crustal deformation pattern and in evaluating the regional seismic potential. Specifically, although several M≥6 earthquakes have occurred along the Leibo fault zone in the southern segment of the YMYTZ, the late Quaternary activity of the fault zone has not been well determined due to insufficient work as well as subsequent lack of solid evidence. The Leibo fault zone strikes NE-SW and spreads on the southeast flank of the Chenqiangyan-Shanzhagang anticline. It starts at the Huanglang township near the Mahu Lake, cuts through the Jingkou Mountain, Lianhuashi, and Leibo, and extends southwestwards to the vicinity of Lianlajue. The latest investigation shows that the Leibo fault zone consists of four subparallel right-lateral strike-slip faults named F1—F4 from the north to the south, respectively. These fault branches together constitute a 43km-long and 10km-wide structural belt. Previous paleoseismic work along the Leibo fault zone found that the faults ruptured the late Pleistocene sedimentary layers with their upward terminations covered by the undeformed Holocene deposits, implying it was active in the late Pleistocene and has not been active since the Holocene. However, the ground surface traces of the Leibo fault zone are the most obvious among the faults in the YMYTZ, and recent seismologic studies show strong seismic activity for the Leibo fault zone, bringing out a controversy about whether it is active in the Holocene or not.

    To address the late Quaternary activity of the Leibo fault zone, we conducted detailed indoor deformed geomorphic feature interpretation on remote sensing imageries like 2m-resolution GF-2 imagery and high-resolution imageries on Google Earth, and further mapped the fault traces in the field using an unmanned aerial vehicle(UAV)derived digital orthographs and digital surface models(DSM). Based on the geological and geomorphological surveys, two trenches were excavated at Pengjiashan and Luohangou along the northern(F2)and southern(F4)branches of the Leibo fault zone respectively. On the trench walls, surface-rupturing paleoearthquakes were identified for each fault according to criteria for faulting events like cut-and-cover structures, scarps, and colluvial wedges. Subsequently, we collected and dated several radiocarbon samples from the sedimentary layers immediately before and after the rupturing events, and finally carried out stratigraphic sequence calibration using the acquired ages with the OxCal 4.4 program to constrain the timings of the revealed paleoearthquakes.

    According to the identification criteria of paleoseismic events, it was revealed 3 paleoearthquakes in the Pengjiashan trench on the northern fault branch(F2)and another 7 rupturing events in the Luohangou trench along the southern fault branch(F4). Radiocarbon sample dating constrain the ages of the paleoearthquakes along F2 to be 21190—20590BC(EP1), 20550—12120BC(EP2), and after 12090BC(EP3), while the latest two paleoseismic events on F4 occurred 9270—5040BC(EL6)and after 5000BC(EL7). Our paleoseismic studies show that the LFZ has experienced several surface-rupturing earthquakes in the Holocene, verifying it is a Holocene active fault zone. Moreover, the ages of the paleoseismic events revealed on two fault branches do not overlap with each other, suggesting they are different paleoearthquakes so that the fault branches in the Leibo fault zone are independent seismogenic structures. By collecting and analyzing the magnitudes of strike-slipping earthquakes that have generated surface ruptures in western China since 1920, it is believed that the minimum magnitudes of the paleoearthquakes determined on the Leibo fault zone are 6.5. Through the empirical relationships between magnitude and surface rupture length, it is estimated that the LFZ has the capability to produce an earthquake with M≥7.

    Table and Figures | Reference | Related Articles | Metrics
    STUDY ON THE EFFECT OF EXCESS TOPOGRAPHY ON LANDSLIDES INDUCED BY LUDING MS6.8 EARTHQUAKE IN 2022
    QIU Heng-zhi, MA Si-yuan, CHEN Xiao-li
    SEISMOLOGY AND GEOLOGY    2024, 46 (4): 783-801.   DOI: 10.3969/j.issn.0253-4967.2024.04.002
    Abstract176)   HTML13)    PDF(pc) (14792KB)(63)       Save

    Strong earthquakes in mountainous regions are prone to triggering severe geological disasters, such as landslides, collapses, and debris flows. These disasters are characterized by their wide distribution, large scale, and high frequency, making them among the most destructive secondary effects of earthquakes. In recent years, the central and eastern parts of the Qinghai-Tibet Plateau have experienced frequent strong earthquakes, leading to varying degrees of earthquake-induced landslide disasters.
    Landscape evolution is significantly influenced by tectonic activities and river incising, which alter the materials of hillslopes and their topographic characteristics. Earthquakes can significantly affect the spatial distribution of co-seismic landslides, particularly in areas with excess topography. In tectonically active zones, rock uplift and river erosion gradually increase slope angles and decrease slope stability. Weathering, freeze-thaw cycles, and seismic vibrations can reduce rock strength, leading to slope erosion and landslides. Once these processes occur, the slope tends to reach a critical state of stability, which is characterized by the presence of excess topography. Landslides can rapidly reduce hillside elevations, limiting terrain relief and impacting landform evolution. Excess topography, defined as rock mass inclined at angles greater than a specified threshold, is closely related to unstable slope masses. The essence of a landslide is the disruption of slope equilibrium, often reflected in the presence of excess topography. However, the influence of excess topography on the distribution of co-seismic landslides is not well understood.
    Earthquake-induced landslides occur when slopes become unstable and slide due to seismic forces. The instability arises when ground motion alters the internal friction angle and cohesion forces along rock mass defects, making them unable to resist the gravitational forces that cause sliding. The weight of the slope material plays a crucial role in this process. As a key component of landscape evolution, landslides significantly shape geomorphic forms, as indicated by the presence of excess topography. The undulating terrain of a region is the result of long-term structural and surface erosion interactions, as well as material migration and distribution. Landslide development is closely related to the local environment, particularly geomorphic conditions. Seismic landslides also play a vital role in shaping and reorganizing active orogenic belts, influencing subsequent landscape evolution.
    On September 5, 2022, a MS6.8 earthquake struck Luding county, Sichuan province, China, with the epicenter in Hailuogou Glacier Forest Park(29.59°N, 102.08°E), at a focal depth of 16km and a maximum intensity of Ⅸ degrees. The earthquake, lasting approximately 20 seconds, was strongly felt across many parts of Sichuan Province and induced numerous landslides, causing significant damage. The affected area, located at the transition between the Qinghai-Tibet Plateau and the Sichuan Basin, features a typical alpine and canyon landscape with steep terrain and river incision, providing favorable conditions for landslides. The long-term and intense tectonic activity in the eastern Qinghai-Tibet Plateau has resulted in complex topography and geomorphology, providing the material basis and external conditions for earthquake and landslide disasters.
    With advancements in science and technology, the Digital Elevation Model(DEM)has become widely used in geoscience research. As DEM accuracy improves, its applications have evolved from qualitative descriptions of geomorphic morphology to semi-quantitative and quantitative analyses of various geomorphic parameters. Geomorphic parameters reveal the structural geomorphic information within topography, essential for understanding regional characteristics and evolution mechanisms. The Luding earthquake serves as a case study for analyzing the influence of topography on the distribution of co-seismic landslides.
    In this study, through post-earthquake remote sensing image analysis we identified 1 485 landslides(covering approximately 14.83km2)and analyzed their spatial distribution. Field surveys revealed that most co-seismic landslides are shallow, with relatively small thicknesses, primarily located along the Dadu River. Excess topography calculations based on the ALOS 12.5m DEM and subsequent quantitative analysis of its correlation with co-seismic landslides indicate a strong relationship: with a 30° threshold, excess topography peaks are found along the Dadu River and its tributaries, coinciding with the majority of landslide occurrences. A total of 91.7% of co-seismic landslides are within areas of varying excess topography heights. However, the average height of excess topography in landslide areas(~80m)significantly exceeds landslide thicknesses, suggesting that the Luding earthquake only mobilized a small fraction of the total excess topography. The remaining excess topography may represent potential unstable slopes for future landslides. Furthermore, the spatial distribution of landslides induced by previous earthquakes, such as the Wenchuan earthquake in 2008, the Lushan earthquake in 2013, and the Ludian earthquake in 2014, shows a high degree of consistency. This underscores the importance of understanding the relationship between excess terrain and landslide distribution to enhance the accuracy of earthquake-induced landslide predictions.

    Table and Figures | Reference | Related Articles | Metrics
    EVIDENCE OF HOLOCENE ACTIVITY OF NALATI FAULT ZONE WITHIN THE TIANSHAN
    WANG Lei, REN Zhi-kun, HE Zhong-tai, JI Hao-min, LIU Jin-rui, GUO Long, LI Xing-ao
    SEISMOLOGY AND GEOLOGY    2024, 46 (4): 821-836.   DOI: 10.3969/j.issn.0253-4967.2024.04.004
    Abstract172)   HTML35)    PDF(pc) (18639KB)(100)       Save

    The Tianshan orogenic belt, extending across the Euro-Asian plates, is one of the most significant intracontinental orogenic belts globally. Spanning over 2 500km, it traverses China, Kazakhstan, Kyrgyzstan, and Uzbekistan from east to west. The belt has been continuously uplifted due to the collision of the India-Eurasia plate during the Cenozoic era. The Tianshan is divided into three segments: North Tianshan, Middle Tianshan, and South Tianshan. The crust in this tectonic region is being shortened in the north-south direction, and a series of NEE-trending or NWW-trending strike-slip faults have developed to accommodate the deformation. The Nalati fault zone serves as the collision suture between the Central Tianshan block and the Tarim block and marks the boundary between Central and South Tianshan. This fault zone trends NEE and extends southwest into Kyrgyzstan, connecting to the Nikolayev line. Its eastern segment is located north of the Dayouludusi basin. The north-south shortening rate is approximately 2.0mm/a, and the horizontal strike-slip rate is about 2.9mm/a. Reports indicate a north-south shortening rate of 0.8-1.1mm/a since the late Quaternary, suggesting it is a significant Holocene active fault zone. However, research on this fault zone's activity is limited, with most studies focused on its eastern segment. Research on other sections remains scarce.
    This study focuses on the middle segment of the Nalati fault zone in Tekes county, Ili Prefecture. The Tekes section trends ENE, starting from Qiongkushitai village in the east, passing through Kalatuori, Ayakeaqia, and Kalawenkeer, and reaching Burili in the west, spanning approximately 55km. Methods employed include remote sensing image interpretation, field geological investigation, UAV aerial surveys, trench excavation, Radiocarbon-14 dating, and semi-automatic horizontal dislocation measurement. The main findings are as follows: 1)The linear geomorphological features of the Tekes segment are prominent, with typical fault geomorphological signs such as fault cliffs, triangles, scarps, bulges, gate ridges, passes, guanmen mountains, and left-lateral dislocation ridges and gullies widely observed.; 2)Small Unmanned Aerial Vehicle Mapping and LaDiCaoz semi-automatic dislocation measurement and analysis indicate a minimum horizontal displacement of approximately 3.4m; 3)Faults are developed in the Proterozoic and Paleozoic strata. A trench 4m long and 1.6m wide excavated at a series of fault reverse scarps revealed a sedimentary event of the sag pond at the hillside, indicating at least four paleo-earthquake events; 4)To date the paleo-earthquake events, we collected 11 sediment samples for Radiocarbon-14 dating at the BETA Analytic laboratory. Results show that the sample at a depth of 2m is about(7.06±0.03)ka BP, and the latest colluvial wedge is about(1.67±0.03)ka BP; 5)Using OxCal age correction, the ages of the four paleo-earthquake events were determined at a 95.4%confidence level: event E1 occurred between 2757BC and 413AD, event E2 between 3581BC and 429BC, event E3 between 4702BC and 3932BC, and event E4 between 5742BC and 5230BC. In summary, we propose that the middle segment of the Nalati fault zone has been active since the Holocene.

    Table and Figures | Reference | Related Articles | Metrics
    CUMULATIVE DISLOCATIONS OF RIVER TERRACES IN THE EASTERN SEGMENT OF THE ALTYN TAGH FAULT
    LI Lu-yao, DING Rui, JIANG Da-wei, ZHANG Shi-min
    SEISMOLOGY AND GEOLOGY    2024, 46 (3): 547-569.   DOI: 10.3969/j.issn.0253-4967.2024.03.003
    Abstract170)   HTML24)    PDF(pc) (18302KB)(133)       Save

    The Altyn Tagh fault zone is a sizeable sinistral strike-slip fault on the northern margin of the Qinghai-Tibet Plateau, and its eastern section is obliquely connected with the northwest Qilian Mountain thrust fault zone. The left-lateral strike-slip action of the Altyn Tagh fault zone and the Qilian Mountain thrust belt constitute a structural transformation relationship. The activity behavior of this fault, especially the amount of sinistral dislocations, segmentation, and slip rate, has always been a hot topic of discussion among scholars. At present, based on geological methods and geodetic research, the slip rates of different sections of the Altyn Tagh Fault have been obtained, with a time scale ranging from tens of thousands of years to decades. The research results generally support the gradual attenuation of the slip rate of the Altyn Tagh fault zone from about 93°E to the east, indicating that its left-lateral attenuation is absorbed by a series of NW-trending thrust faults on its eastern side, and this trend has changed little for decades. The above work provides a framework for us to study the structural transformation relationship between the Altyn Tagh fault zone and the Qilian Mountains thrust belt in time and space. However, due to the limitations of previous observation points, especially the different methods of fault geomorphology measurement and dating used by various authors, there are significant differences in the obtained fault slip rates. Currently, it is not possible to analyze the segmented characteristics of the slip rate in the east section of the Altyn Tagh fault zone further and its spatial relationship with the Qilian Mountain thrust structural zone based on this. In recent years, the application of UAV aerial survey technology has allowed image data to be obtained at the centimeter to millimeter level, making the study of tectonic geomorphology more refined. Researchers can obtain the cumulative displacement or co-seismic displacement of several seismic cycles through micro-faulted landforms and reconstruct the dislocation accumulation process of active faults.

    The transfixion terraces developed across the Altyn Tagh fault are mainly controlled by regional tectonic uplift and climate change, showing regional synchronization in time, which provides convenience for the comparison of regional landforms. Although there are differences in the grading standards of terraces at different sites and the dating methods are not completely consistent, the chronological sequence of the late Quaternary river terraces in the study area generally shows good consistency. From new to old, it is about 3-4ka BP, 6-8ka BP, 10-13ka BP, 20-21ka BP and 40-50ka BP, which provides a research basis for our subsequent comparison of the displacement amount of river terraces. Based on the high-resolution image data obtained by unmanned aerial vehicle photogrammetry(SfM), this paper carried out a detailed interpretation of the 127km section of the eastern section of the Altyn fault zone and measured and counted the dislocations of different levels of risers at 9 typical river terrace dislocation points. Based on the distribution of cumulative displacement of the same terrace, the kinematic segmentation characteristics and tectonic mechanism of the eastern section of the Altyn Tagh fault zone are discussed.

    To define the displacement more accurately, we considered the following factors: 1)Usually, the riser is not a straight line but a curve formed by the free swing of the river bed; 2)Strictly speaking, a riser is a slope with a certain width, consisting of an upper edge, a lower edge, and a middle slope zone. When measuring displacement, we cannot only measure the upper edge or lower edge of the scarp but must consider it comprehensively. Based on the above prerequisites, this paper uses two envelope lines to surround the upper and lower edges of the riser and measures the distance between the corresponding envelope lines on both sides of the fault to obtain the displacement of different levels of scarp in the river terrace. Based on the above measurement methods, four dislocation values, A1-A4 and B1-B4, were obtained from the curve envelope of the upper and lower edge of the terrace scarp on both sides of the fault, respectively. After calculating the corresponding mean value, the standard deviation of the four measured values was estimated to reflect the dispersion degree of different measured values relative to the mean value, which was used as the error range of the measurement results.

    The results show that the late Quaternary left-lateral cumulative displacement tends to decline along the Altyn fault zone eastward. The cumulative dislocations are roughly the same in the same fault segment, which may indicate the consistency of seismic ruptures within the segment. In addition, the dislocation amount of the same-level terrace scarp between adjacent fault segments shows a stepwise decrease, indicating the tectonic transformation relationship of the miter fault and the possible seismic rupture segmentation, which provides a basis for the active segmental research and potential earthquake-generating capacity evaluation of the Altyn fault zone. A horizontal dislocation of 2-3m occurred on the front scarp of the T1 terrace in the Gaoyangou area, indicating that the latest earthquake surface dislocation event may have occurred in the Hongliugou-Shaping section at the easternmost end of the Altyn Tagh fault zone.

    Table and Figures | Reference | Related Articles | Metrics
    REMOTE SENSING INTEPRETATION OF COSEISMIC LAND-SLIDES TRIGGERED BY 1976 LONGLING MS7.3 AND MS7.4 EARTHQUAKES AND THE TECTONIC SIGNIFICANCES
    LI Hao-feng, XU Yue-ren, GUO Ya-li, LIU Han, ZHAO Xin-yu, LU Ling-yu, TANG Jia-cheng
    SEISMOLOGY AND GEOLOGY    2024, 46 (1): 117-140.   DOI: 10.3969/j.issn.0253-4967.2024.01.008
    Abstract165)   HTML12)    PDF(pc) (17130KB)(119)       Save

    The Longling-Lancang seismotectonic belt in southwestern Yunnan is critical for accurately defining the boundaries of active blocks and evaluating seismic risks. Using pre- and post-earthquake high-spatial-resolution satellite imagery to study strong earthquakes retrospectively proves to be a practical method in such study. Strong earthquakes frequently cause secondary effects such as coseismic landslides, collapses, and debris flows, which lead to considerable loss of life and property. These secondary effects, often as the most dramatic manifestations of an earthquake, show geologic signatures providing evidence of historic or prehistoric seismic activities. The use of satellite imagery captured shortly after historic earthquakes to interpret these secondary effects is particularly beneficial in determining the intensity and influence radius of earthquakes, thereby helping study on seismogenic faults of earthquakes.

    On May 29, 1976, two strong earthquakes with MS7.3 and MS7.4 occurred in Longling county, southwest China, followed by intense aftershocks. The seismogenic structure of these earthquakes still remains undetermined to present. These earthquakes triggered numerous coseismic landslides in the regolith of the granitic rock mass. The seismic zone, located in subtropical regions, is characterized by high precipitation and dense vegetation. Apart from the ancient landslides in the northwest and southeast, no records of landslides and debris flows persisted in the epicenter zone for a century, making the occurrence of substantial landslides post the main earthquakes unexpected. Currently, these landslides have undergone reshaping by land surface processes and re-vegetation, which makes them indistinguishable in recent remote sensing images. Using Keyhole satellite images with a resolution of 0.6~1.2m offers a useful means to identify the coseismic landslides of the Longling mainshocks. In this study, we employ these images for a comprehensive visual interpretation of the coseismic landslides. To ensure the accuracy and reliability of the results, we used images captured in 1981(the most recent following the earthquakes)to extract coseismic landslides and substantiated them with images from 1974, field investigation photos from 1976, and relevant records. Finally, we have compiled an exhaustive database of coseismic landslides triggered by the 1976 Longling cases.

    Our results are summarized as follows: 1)A total of 14 448 landslides were interpreted, encompassing an overall area of 17.2km2. The area of individual landslides primarily ranged from several hundreds to one thousand m2, and most were superficial slides in the surface regolith with short sliding distances. The regional stratigraphy is complex, with 70.9% of the landslides occurring in the regolith of granitic rock mass, 15.3%in sandstones or siltstones, and a mere 13.8%in other areas such as limestones. Consequently, these landslides were relatively small compared to those in other regions like the Loess Plateau in north China, where the surface sediment is extremely loose. 2)A strong correlation exists between the intense area of coseismic landslides and the earthquake sequence, which tends to migrate from south to north. Notable aftershocks(e.g., MS6.2 on June 9 and MS6.6 on July 21)particularly exhibited the general NNW distribution direction of the earthquake sequence and triggered scattered landslides outside the epicenter zone. Through synthesizing field surveys, combining other records and the findings of this study, we believed that the two main earthquakes triggered numerous coseismic landslides, and the continuous strong aftershocks led to the destabilized regolith of the granitic rock mass creeping successively, resulting in subsequent landslides. 3)The concentration areas of coseismic landslides do not match the high-earthquake-intensity areas, instead, they are all located on one side of active faults, which suggests that the seismogenic fault is neither the Longling-Ruili Fault nor the Wanding Fault. The spatial distribution of the landslides suggests that the scope of the surface rupture zone is about 30km. The conjugated strong earthquake ruptures in southwestern Yunnan may limit the spatial scale of single strong earthquakes, so it is crucial to pay more attention to the intersection zone of NE and NW trending active faults when assessing regional strong earthquake risk in the future.

    Table and Figures | Reference | Related Articles | Metrics
    THREE-DIMENSIONAL CRUSTAL VELOCITY STRUCTURE AND SEISMOGENIC ENVIRONMENT AROUND THE HUOSHAN EARTHQUAKE SWARM
    JI Guo-qiang, LEI Jian-she, ZHAO Da-peng
    SEISMOLOGY AND GEOLOGY    2024, 46 (3): 665-685.   DOI: 10.3969/j.issn.0253-4967.2024.03.009
    Abstract165)   HTML15)    PDF(pc) (8799KB)(83)       Save

    The Huoshan earthquake swarm is tectonically located at the junction among the North China plate, Yangtze plate, and North Dabie orogenic belt. The geological environment in the region is complex, including the Feixi-Hanbaidu Fault, Meishan-longhekou Fault, Xiaotian-Mozitan Fault and Luoerling-Tudiling Fault, as well as the North Dabie tectonic belt, North Huaiyang tectonic belt and Liu’an basin. In the study region, seismicity is intense, and 9 earthquakes with M≥5.0 occurred along the Luoerling-Tudiling Fault in the history. In recent decades, small-to-moderate earthquakes were frequent, mainly gathering at the intersection of the Xiaotian-Mozitan Fault and Luoerling-Tudiling Fault. Furthermore, the frequency of small earthquakes in the Huoshan region has a significant correspondence to the moderate-to-strong earthquakes in East China and even eastern Tibet, so studying the deep structure can shed new light on the relationship between the Huoshan earthquake swarm and moderate-to-strong earthquakes in mainland China.

    In this study, a total of 17 920 seismic arrival-time data, including 7 706 P, 394 PmP, 9 263 S and 557 SmS arrivals, are hand-picked from 1987 local earthquakes to obtain three-dimensional crustal P-wave velocity(VP), S-wave velocity(VS)and VP/VS ratio models down to 30km depth beneath the Huoshan swarm area. The checkerboard resolution test results show that the imaging spatial resolution in most parts of the regions can reach 0.33°×0.33°, and the North Huaiyang tectonic belt near the Huoshan earthquake swarm has good recovery in the entire crust, and the North Dabie tectonic belt and Lu’an basin also have good recovery at 8-30km depths. Due to the addition of PmP/SmS arrivals, the spatial resolution at 18-30km depths is significantly improved, and the pattern and amplitude of velocity anomalies are better recovered.

    Our tomography results show that a vertical continuous high VP/VS anomaly is observed around the intersection of the Xiaotian-Mozitan Fault and the Luoerling-Tudiling Fault, especially at 18km depth appear broad low-velocity and high VP /VS anomalies. At 30km depth, the areas with high VP/VS are reduced and concentrated on both sides of the Luoerling-Tudiling Fault. There are significant high VP/VS characteristics around the Huoshan earthquake swarm. The high VP/VS anomalies extend to 18~30km depths below the Xiaotian-Mozitan Fault, suggesting that fluids could have migrated upward along the fault to reduce the effective normal stress of the fault planes, triggering the activity of the Huoshan earthquake swarm at the weak intersection between the Xiaotian-Mozitan Fault and the Luoerling-Tudiling Fault.

    Combined with the low-velocity anomalies of the upper mantle revealed by the previous tomographic results, we propose that there may be a channel for upwelling of the wet and hot upper-mantle materials with fluids to the crust along the Xiaotian-Mozitan Fault. The upwelling of the wet and hot materials may be related to the dynamics of the big mantle wedge formed due to the deep subduction of the stagnant Pacific slab down to the mantle transition zone and the eastward extrusion of materials in the upper mantle from eastern Xizang along the Dabie orogenic belt. These factors may jointly affect the seismicity characteristics of the Huoshan earthquake swarm. Our results providea new piece of seismological evidence for the interactions among the tectonic activities in the Huoshan region, Tibetan plateau and East China.

    Table and Figures | Reference | Related Articles | Metrics
    THE SPATIAL AND TEMPORAL CHARACTERISTICS OF PRESENT-DAY SEISMICITY IN NORTHEASTERN LONGMENSHAN FAULT ZONE
    HU Nan, LONG Feng, WANG Ying, XU Liang-xin
    SEISMOLOGY AND GEOLOGY    2024, 46 (4): 856-875.   DOI: 10.3969/j.issn.0253-4967.2024.04.006
    Abstract163)   HTML16)    PDF(pc) (6435KB)(101)       Save

    The Longmenshan fault zone, situated along the eastern margin of the Tibetan plateau, represents a significant thrust tectonic belt characterized by pronounced segmentation. It is delineated into northern and central-southern segments at Beichuan, and along its depth, it features three sub-parallel fault belts: the Houshan fault, the Central fault, and the Qianshan fault, extending from the northwest to the southeast. Geological research indicates that since the Quaternary, the central-southern segments of the Longmenshan fault zone have exhibited considerable seismic activity, whereas the northern segment has shown minimal signs of movement. However, paleo-earthquake studies have identified substantial historical seismic events in the Qingchuan fault, a component of the northern segment, dating back to the Holocene. The devastating 2008 Wenchuan earthquake(MS8.0), which occurred in the middle section of the Longmenshan fault zone, resulted in a 240-km-long surface rupture along the Central fault. A multitude of aftershocks radiated northward from the epicenter, with no discernible surface ruptures observed in the northern segment. This study aims to provide a comprehensive analysis of the kinematic features of the northern segment by re-evaluating the Wenchuan earthquake's aftershocks and employing focal mechanisms derived from previous studies.
    Seismic activity is intrinsically linked to active tectonics, and the precise localization of minor earthquakes can offer critical insights into the underlying seismogenic processes and mechanisms. In this paper, we have compiled early aftershock relocation data and further refined the relocation of small earthquakes using an integrated seismic location technique. Seismic phase data were obtained from the networks in Sichuan, Gansu, and Shaanxi over the past decade, spanning from 2010 to 2020. To mitigate the impact of crustal velocity variations, an optimal one-dimensional velocity model for the study area was initially inverted using the VELEST program. The Hypo2000 program was then utilized to adjust the initial seismic source positions, followed by the application of the double-difference method for the relocation of minor earthquakes. The reliability of the localization outcomes, determined using the LSQR method, was verified by the SVD method. Consequently, 10 653 minor earthquakes were relocated with an average travel time residual of 0.053s, a horizontal location error of 281m, and a vertical location error of 260m.
    In the southern extremity of the study area, the relocated earthquakes are predominantly aligned along the parallel faults flanking the primary rupture zone. In the south-central region, the relocated earthquakes exhibit deviations from the rupture zone, revealing multiple seismic clusters. Towards the northern end, the relocated earthquakes demonstrate a migration from the main rupture towards the Qingchuan fault. The depth profiling of seismic sources reveals that the relocated earthquakes are concentrated between 8-15km deep, all situated above the 500℃ isothermal surface. The depth profile in the southern region continues the characteristics of the main rupture surface of the Wenchuan earthquake, while the dip angle becomes increasingly steep as it progresses northward. The northern end's depth profile suggests an interaction between the rupture surface and the Qingchuan fault. Additionally, the analysis of 32 focal mechanisms exceeding ML4.0 within the study area corroborates the geometrical structures of the fault zone, as revealed by the spatial distribution of the relocated earthquakes, further validating the reliability of relocation.
    A comprehensive analysis suggests that the current seismicity in the northern section of the Longmenshan fault zone is multifaceted, with ongoing activity on the main rupture surface(afterslip)and slip on secondary new rupture surfaces triggered by the mainshock. It is hypothesized that the spatial distribution of the relocated earthquakes retains segmented characteristics. In the southern region of the study area, thrust slip induced by the main rupture continues; in the middle region, new ruptures are concurrently active with the main rupture; and in the northern region, influenced by the high velocity of the upper crust around Ningqiang-Mianxian, the rupture zone vanishes at the surface, with the deep triggering of the Qingchuan fault by stress transfer being evident. In conclusion, the complex spatial characteristics of the current seismic activity in the northern section of the Longmenshan are attributed to the interplay of pre-existing faults, new ruptures, and the main rupture, reflecting the spatially heterogeneous process of stress transfer and adjustment following the Wenchuan earthquake, potentially linked to the complex geological structure of the region.

    Table and Figures | Reference | Related Articles | Metrics
    3D MODELING AND MAXIMUM POTENTIAL SEISMIC ASSESS-MENT OF THE EASTERN MARGIN FAULT OF DAXING UPLIFT
    ZHANG Ya-jing, LI Zheng-fang, ZHOU Ben-gang, XIAO Hai-bo
    SEISMOLOGY AND GEOLOGY    2024, 46 (4): 802-820.   DOI: 10.3969/j.issn.0253-4967.2024.04.003
    Abstract162)   HTML27)    PDF(pc) (6731KB)(92)       Save

    The eastern margin fault of Daxing uplift is an important boundary fault in the southeast of Beijing Plain. The fault is located in the southeast of properthe Xiadian Fault and is distributed in the correct order en echelon with the Xiadian Fault, which controls the development of Langgu secondary depression under the extensional tectonic background. Recent shallow seismic reflection profiles and borehole data have found evidence of Holocene activity in the eastern margin of the Daxing Uplift, which has changed the conclusion in recent decades that it has not been active since the late Quaternary. Because the fault is a right order echelon with the Xiadian Fault, and it is similar to the Xiadian Fault in structure, and the Xiadian Fault had the Sanhe-Pinggu M8 earthquake in 1679, it is inferred that the fault has the risk of a large earthquake. It has essential crucial application value to the seismic hazard survey in Beijing. Also, it poses a new challenge to the upper limit of the maximum potential earthquake magnitude of the fault on the eastern margin of the Daxing Uplift.
    Quaternary sediments cover the fault on the east margin of Daxing Uplift and are in a hidden state, which results in its geometric features and deep and shallow coupling relationships that cannot be visually demonstrated by two-dimensional data two-dimensional data cannot visually demonstrate. It is of great significance to establish a three-dimensional model of hidden active faults for the hazard assessment of seismic active faults. In this paper, by collecting the fine location data of small earthquakes in this area and collating several shallow seismic geophysical profiles and deep seismic reflection profiles, SKUA-GOCAD 3D geological modeling software is used to build 3D models of the eastern margin fault of Daxing Uplift and the Xiadian Fault based on the section modeling method, and the distribution of the two faults in 3D space is simulated. The geometric features and the relationship between the depth and shallow structure of the two faults are revealed, including 1)a three-dimensional fault model and stratigraphic information map; 2)a three-dimensional model diagram of fault distribution according to dip Angle; 3)Three-dimensional model diagram of fault distribution according to depth and a three-dimensional map of small earthquake distribution. The 3D map shows that there are strong structural similarities between the faults on the eastern margin of the Daxing uplift and the Xiadian faults. The contrast map shown by depth shows that both faults are deep and shallow faults, the shallow faults disappear at about 15km underground, and the deep faults extend downward to cut the lower crust and the Moho surface. The contrast diagram displayed by apparentthe dip Angle clearly reflects that the two faults have obviously different dip angles in-depth and shallow. The deep fault is almost steep, and the shallow fault shows obvious differences in different sections. The distribution range of small earthquakes is 0-25km, of which the dominant distribution range is 10-20km. Therefore, it is speculated that the east margin fault of Daxing Uplift may have the seismogenic capacity similar to the Sanhe-Pinggu M8 earthquake in 1679. However, as existing studies have shown that the activity of the Xiadan fault and its southern extension section-eastern margin of the Daxing Uplift in this region gradually weakens from north to south, the maximum potential earthquake magnitude of the east margin fault of the Daxing Uplift is inferred in this paper to be less than Sanhe-Pinggu M8 earthquake in 1679.
    Finally, by using the structural analogy of the Xiadian Fault on the eastern margin of the Daxing Uplift, and based on the structural similarity of the two faults, this paper evaluates the maximum potential earthquake magnitude that may be induced by the Daxing Fault using different experiential relations of magnitude-fault rupture scale fitted by predecessors in North China. The conclusion is as follows: the distribution range of the magnitude of the earthquake is 7.3-7.4.
    Based on the structural analogy with the Xiadian Fault and the empirical relationship between magnitude and rupture scale, the maximum potential earthquake magnitude induced by the eastern margin fault of Daxing uplift is estimated to be magnitude 7.5. This conclusion has important scientific guiding significance for earthquake disaster prevention and control in the capital area, and should be paid attention to and actively take prevention and avoidance measures.

    Table and Figures | Reference | Related Articles | Metrics
    THE CHARACTERISTICS AND MECHANISM OF FLUID ANOMALIES IN THE DAZHAI OBSERVATION WELL OF PU’ER, YUNNAN PROVINCE BEFORE THE M5.9 MOJIANG EARTHQUAKE ON SEPTEMBER 8, 2018
    HU Xiao-jing, FU Hong, ZHANG Xiang, LI Li-bo, HUANG Jiang-pei, LI Qiong, GAO Wen-fei
    SEISMOLOGY AND GEOLOGY    2024, 46 (2): 477-491.   DOI: 10.3969/j.issn.0253-4967.2024.02.014
    Abstract161)   HTML12)    PDF(pc) (6351KB)(166)       Save

    The precursors before earthquakes are very useful to earthquake prediction, and fluid anomalies before earthquakes are very important to precursory observations. This paper reviews the characteristics of hydrochemical ions and well-aquifer permeability anomalies of the Dazhai observation Well in Pu’er, which is in the Yunnan-Southwestern region of China, for all M≥5.5 earthquakes since 2004. We find that both the chemical ions and physical parameters before the Mojiang M5.9 earthquake exhibited the largest magnitude of changes since observation, and the abnormal state was much stronger than that of previous historical earthquakes, but the magnitude of the earthquake was below 6. About 1.5-2a before the M5.9 Mojiang earthquake, the composition of hydrogen and oxygen isotopes of the water samples in the Dazhai observation Well showed a significant deviation, accompanied by a continuously increasing concentration of fluoride ions from sources at deeper depths. This might suggest that the deep material in the earthquake source area began to be active. At the same time, starting one year before the earthquake, the phase lag of the water level in the wellhole changed from negative to positive, indicating that the source and pathway of well water recharge have been changed. In addition, around half a year before the earthquake, the continuously observed water chemical ions at shallow depths in the wellhole began to show a dramatic change. Moreover, macroscopic anomalies of hot spring water volume increased sharply before the earthquake, showing a remarkable evolution process from deep to shallow, from background to short-term, and from micro anomalies to macro anomalies before the earthquake. To investigate the causes and mechanisms of this phenomenon, we attempt to discuss the abnormal evolution process before the M5.9 Mojiang earthquake from the aspects of regional deep material activity and regional stress level. The abnormal concentration of the hydrochemical ions and the change of aquifer permeability observed continuously at the Dazhai observation well before the M5.9 Mojiang earthquake were caused by the continuous increase in shear stress in the region, which caused the aquifer to be compressed, resulting in a vertical fluid recharge and ultimately the alternation and mixing of different aquifer water bodies. In addition to being controlled by the continuous increase in regional vertical shear stress, the abnormal formation process was also accompanied by the intense activity of deep-sourced chemical elements such as helium isotope and fluoride ion. The abnormal evolution process showed a remarkably coupled process of migration from deep to shallow, which may be the reason why the shallow ion anomaly before the M5.9 Mojiang earthquake was the most significant among all the observed cases. Therefore, the evolution process of fluid activity starting from the deep and continuously transmitting to the surface with the accumulation of regional stress is essential to the abnormal evolution of the hydrological phenomenon before the M5.9 Mojiang earthquake. The regional stress and the process of deep material activity are the biggest differences between the M5.9 Mojiang earthquake and other historical earthquake cases in the study area, which will be the two main factors to be considered when similar ion changes occur again in the future. Our study provides insight into a comprehensive understanding of the predictive significance of underground fluid anomalies in the Dazhai well and the coupled evolution process of deep-shallow fluid anomalies before the earthquake.

    Table and Figures | Reference | Related Articles | Metrics
    SEISMOLOGY AND GEOLOGY    2024, 46 (3): 756-759.  
    Abstract159)   HTML40)    PDF(pc) (1214KB)(127)       Save
    Related Articles | Metrics
    STUDY ON SURFACE WAVE TOMOGRAPHY OF THE A'ERSHAN VOLCANOES
    HOU Jie, WU Qing-ju, YU Da-xin, YE Qing-dong
    SEISMOLOGY AND GEOLOGY    2024, 46 (4): 893-915.   DOI: 10.3969/j.issn.0253-4967.2024.04.008
    Abstract157)   HTML8)    PDF(pc) (11942KB)(34)       Save

    Since the Cenozoic era, a series of intraplate volcanic groups have developed along the east and west sides of the Songliao Basin in the eastern part of the Central Asian orogenic belt. The A'ershan volcanic group is one of the Cenozoic intraplate volcanoes in the eastern section of the Central Asian orogenic belt. Further study of this volcanic group is of great significance for exploring and understanding the genesis of intraplate volcanoes in the eastern section of the Central Asian orogenic belt. In the past, the distribution of mobile seismic stations and some fixed stations used for imaging research on the A'ershan volcanic group was relatively sparse and did not fully cover the A'ershan volcanic group. The resolution of the crust-mantle structure obtained in the past was also slightly insufficient for exploring the genesis mechanism of the A'ershan volcanic group. This article utilizes the vertical teleseismic waveforms of 29 broadband mobile seismic stations near the A'ershan volcanic group and 8 fixed stations around them from May 2019 to December 2021. Through frequency-time analysis technology, 11 775 Rayleigh wave phase velocity dispersion between two stations is extracted. After excluding non-monotonic rise and phase velocity dispersion curves that differ significantly from most dispersion distributions, 11 010 high-quality Rayleigh wave fundamental phase velocity dispersion curves were ultimately obtained. Then, based on classical ray theory, the two-dimensional phase velocity distribution with a period of 10-80s and a grid size of 0.5°×0.5° is inverted by using the traditional dual station method. Except for areas not covered by radiation in the edge zone, the lateral spatial resolution of phase velocity in the study area is basically within 50km. The checkerboard test also showed that dividing the grid size of the study area into 0.5°×0.5° is feasible, and anomalies with a central area scale less than 0.5°×0.5° can also be identified. Afterward, the CRUST1.0 model was used as the initial crustal model, and the PREM model was used as the initial mantle model. The crustal thickness results obtained from the receiver function were used to constrain the thickness of each layer in the initial crustal model, further reconstructing the three-dimensional S-wave velocity structure of the crust and upper mantle in the study area. The results show that: within the range of the middle and lower crust, the S-wave velocity in the A'ershan volcanic area exhibits apparent low-velocity anomalies. Based on the characteristics of the high wave velocity ratio in the area, it is speculated that there may be a crustal magma chamber in the A'ershan volcanic group. There are multiple high-velocity anomaly structures within a depth range of 40-150km in the A'ershan volcanic group. The difference in the depth of high-velocity anomalies indicates the heterogeneity of the lithosphere thickness, and it is speculated that the thickness of the lithosphere in the A'ershan volcanic area does not exceed 100km. The deeper distribution of high-velocity anomalies may represent the dismantled lithosphere, while the shallower distribution of high-velocity anomalies may represent the undeveloped lithosphere or residual lithosphere after dismantling, reflecting the possibility of lithospheric detachment and subsidence in the region. There are S-wave low-velocity anomalies in the upper mantle on the north and south sides of the A'ershan volcanic group, connecting the asthenosphere and the exposed positions of the A'ershan volcanic group on the surface. The low-velocity anomalies on the north and south sides merge at a depth of 150km. Based on the high heat flux value, high VP/VS, and crustal thinning characteristics of the surface near the distribution area of the A'ershan volcanic group, as well as the previous conclusion based on remote seismic P-wave and S-wave travel time tomography results that there is a clear connection between the low-velocity anomaly below the A'ershan volcanic group and the southern edge of the Songliao Basin in the deep mantle, it is speculated that this low-velocity anomaly may be caused by the upwelling of asthenosphere material caused by the detachment of the lithosphere in the Songliao Basin.

    Table and Figures | Reference | Related Articles | Metrics
    CONSTRUCTION METHOD OF SEISMOTHERMAL INFRARED BACKGROUND FIELD BASED ON GPR-LSTM
    SONG Dong-mei, ZHANG Man-yu, SHAN Xin-jian, CUI Jian-yong, WANG Bin
    SEISMOLOGY AND GEOLOGY    2024, 46 (2): 492-511.   DOI: 10.3969/j.issn.0253-4967.2024.02.015
    Abstract153)   HTML4)    PDF(pc) (11203KB)(101)       Save

    Seismic monitoring is a very important and challenging task. The continuous development of remote sensing technology has strengthened our ability to monitor the Earth’s surface on a macro scale. Research shows that an abnormal rise in surface temperature usually occurs before an earthquake, so a variety of anomaly extraction algorithms have been applied to the study of seismic thermal anomalies. Among them, the extraction method based on background field is widely used because of its strong mechanism interpretation. However, the previous anomaly methods based on the background field mainly limit the background field to a certain fixed threshold value, and ignore the small range of normal LST fluctuations caused by external factors(non-seismic). Therefore, a method of constructing a seismothermal infrared background field based on GPR-LSTM is proposed in this paper. The main idea of the method is that the LST background field is obtained by adding the established annual variable reference field and the fluctuation range of normal LST. First, the LST exhibits a range of fluctuations due to non-seismic factors such as solar radiation, weather, and human activities. Therefore, it is not reasonable to take a fixed value of the LST background field but it should have a certain fluctuation range. Therefore, the dynamic fluctuation characteristics of the background field should be reflected in the construction process of the background field in this study. Secondly, the reason why this method uses the LSTM model to predict the annual variable reference field of the earthquake period based on the annual variable reference field of the non-earthquake period is that on the one hand, the LSTM model can predict the law of long time series data, so it can better learn the annual variation law of LST in the non-earthquake year. At the same time, the LST trend of increasing or decreasing year by year caused by climate change can be predicted, which is conducive to more truly describing the LST trend of the real background field in the year of the earthquake period.

    The method includes two parts: the establishment of the annual variable reference field, the calculation of the residual fluctuation range of the actual LST and the construction of the background field. Based on the MODIS surface temperature product, the precursory thermal anomaly information of the 2008 Wenchuan earthquake in Sichuan Province and the Yutian earthquake in Xinjiang Province was extracted and analyzed by using the proposed method. Based on the MODIS surface temperature data of the 2008 Wenchuan earthquake in Sichuan Province and the Yutian earthquake in Xinjiang Province, the proposed method is used to detect and analyze the thermal anomalies before the earthquake. The following conclusions are drawn:

    (1)The established LST background field on the Qinghai-Xizang Plateau is consistent with the actual variation law of LST and is relatively stable on the whole. In the time dimension, the background field of LST shows the annual variation of high summer and low winter. In terms of spatial dimension, the established LST background field is generally high in the south, low in the middle and low in the north, and the temperature value of the background field is the lowest at the highest altitude of the Qinghai-Xizang Plateau.

    (2)The thermal anomalies obtained based on the new algorithm are usually distributed along the fault zone of the Qinghai-Xizang Plateau, and the anomaly evolution law is obvious. For the Yutian earthquake in Xinjiang, the trend of thermal anomalies is consistent with that of the fault zone in most cases, and the thermal anomalies are mainly distributed in the northern margin of the fault zone. For the Wenchuan earthquake in Sichuan Province, with the onset time approaching, the thermal anomalies gradually moved southward along the Longmenshan fault zone from the north margin, filled the entire fault zone before the earthquake, and disappeared after the earthquake. This proves the validity of the proposed background field construction method.

    (3)Compared with non-seismic years, the spatial characteristics of thermal anomalies along faults are more obvious in seismic years, and the duration and amplitude of the anomalies are longer. This indicates that the tectonic activity in seismic years is more active than that in non-seismic years, resulting in more significant abnormal warming of surface temperature.

    (4)The evolution law of thermal anomalies of earthquakes with magnitude 7 or above on the Qinghai-Tibet Plateau can be summarized as incubation—disappearance—accumulation—disappearance—earthquake occurrence, and the strength of thermal anomalies appears in the way of repeated changes. This indicates that the thermal anomaly is not continuously enhanced or weakened before the earthquake but shows the characteristics of repeated fluctuations between strong and weak until the earthquake eruption thermal anomaly disappears.

    In conclusion, the proposed method can not only ensure that the established background field is not affected by tectonic factors but also fully consider the small range of normal fluctuation of surface temperature caused by non-tectonic factors, which makes the extraction of thermal infrared anomaly information more accurate. This provides new technical means and ideas for the extraction of earthquake anomaly information and is helpful for the further development of earthquake monitoring.

    Table and Figures | Reference | Related Articles | Metrics
    ESR DATING OF CALCITE VEINS AND IMPLICATIONS FOR THE ACTIVITY OF THE JIANCHUAN SECTION OF THE LONGPAN-QIAOHOU FAULT ZONE
    JI Hao, LIU Chun-ru, WEI Chuan-yi, YANG Hui-li, YIN Gong-ming, CHANG Zu-feng
    SEISMOLOGY AND GEOLOGY    2024, 46 (1): 81-100.   DOI: 10.3969/j.issn.0253-4967.2024.01.006
    Abstract153)   HTML18)    PDF(pc) (9015KB)(101)       Save

    The Sichuan-Yunnan region is located in the southeastern part of the Qinghai-Tibet Plateau. Because of the compression and collision dynamics of the Indian Plate and the Eurasian Plate, the tectonic deformation is strong and seismic activities occur frequently. There have been many earthquakes above magnitude 7.0 in history. A series of active fault zones have developed in the region, among which the Sichuan-Yunnan rhombus block bounded by multiple active faults has attracted great research interests in recent years. The Longpan-Qiaohou fault zone is a boundary fault of the Sichuan-Yunnan rhombus block. The fault zone starts from Longpan in the north, passes through Jiuhe, Jianchuan, and Shaxi in the south, and ends at Qiaohou. It is about 120km long and the fault trend is 15°~20°. This fault zone is large in scale and highly active, with frequent seismic activity, complex mechanical properties, and variable movement patterns. The Mesozoic movement was intense. In the early Cenozoic, compression-thrust movement was dominant, and in the late Cenozoic, tension-strike movement was dominant. Since the Holocene, the fault zone has been characterized by left-lateral strike-slip movement with normal faulting properties, and earthquakes of magnitude 5 or above have occurred many times. Therefore, studying the activity of this fault zone is of great significance for the prediction and evaluation of regional strong earthquake risk. Thick calcite veins are well developed on the Henancun Fault of the Jianchuan section of the Longpan-Qiaohou fault zone, providing very valuable materials for fault dating. Calcite veins are coseismic rapid precipitation formed during seismic activity or syntectonic precipitation that filled along fractures after seismic activity. Therefore, their ages represent the latest time at which seismic activity occurred. Previous studies have shown that tensional fissures formed during coseismic events can close in a short period of time(days to months), suggesting that the filling of calcite veins within fault fissures is a relatively rapid process. This paper uses the ESR method to conduct dating study on the calcite veins in the study area. The results show that the ages of the four calcite veins(HNC-ESR01, HNC-ESR02, HNC-ESR03 and HNC-ESR04) are: (7.1±0.8)ka, (7.1±0.9)ka, (7.3±1.7)ka and (6.9±1.5)ka, respectively. The age results are concentrated, and the average age is(7.1±1.3)ka, indicating that the fault was active no later than(7.1±1.3)ka. The age results are consistent within the error range with the second paleoseismic event time revealed by trenching work in the area(between(6 130±30)a BP and(6 320±40)a BP), indicating that the dating of ESR in the fault zone is an effective dating method for the study of active tectonics and paleo-earthquakes. It is an effective chronological method for research, but it can be seen that compared with 14C and luminescence dating, the error of ESR results is relatively large. For faults with short earthquake recurrence intervals, it is still very challenging to accurately judge their activity. In the follow-up work, it is necessary to further improve the experimental process and reduce experimental errors, including refinement of sample pretreatment, accurate monitoring of irradiation dose, and accurate calculation of dose rate. In addition, by using five fitting functions(LIN, SSE, DSE, EXP+LIN and Dgamma)to calculate the equivalent dose values of calcite vein samples in this study, we found that the SSE function is capable of providing the best fitting effect.

    Table and Figures | Reference | Related Articles | Metrics
    DISCUSSION ON REASON WHY AFTERSHOCKS OF THE 2014 LUDIAN EARTHQUAKE OCCURRED ON THE CONJUGATED FAULTS BASED ON TECTONIC STRESS FIELD
    WAN Yong-ge, SONG Ze-yao, GUAN Zhao-xuan, HUANG Rui-qi, GU Pei-yuan, WANG Run-yan
    SEISMOLOGY AND GEOLOGY    2024, 46 (1): 184-200.   DOI: 10.3969/j.issn.0253-4967.2024.01.011
    Abstract151)   HTML17)    PDF(pc) (8278KB)(120)       Save

    An earthquake of MS6.5 occurred in Ludian county, Zhaotong city, Yunnan province, on August 3, 2014. Aftershock activity appears to be on two conjugated near-vertical Zhaotong-Ludian Fault and Baogunao-Xiaohe Fault, which is unusual phenomenon in aftershock distribution. Usually, most earthquakes occur on a single fault or fracture zone with several end-to-end linked faults. Therefore, the 2014 Zhaotong-Ludian earthquake sequence has attracted great attention from geoscientists, and a lot of work has been done in terms of the precise location of aftershocks, the focal mechanisms of the earthquake sequence, co-seismic slip model and the deep structure of the Earth’s crust. In order to determine the fault on which the mainshock occurred and to understand why the aftershocks occurred on the conjugate faults based on the tectonic stress point of view, we used the new developed fuzzy clustering and fault plane fitting method to calculate the geometry parameters of the 2 fault branches. The result shows that the strike and dip angle of the NNW-SSE branch(Baogunao-Xiaohe fault zone)are 336.67° and 88.41°, respectively, and the strike and dip angle of the near EW branch(Zhaotong-Ludian fault zone)are 266.10° and 86.42°, respectively. Both the Zhaotong-Ludian Fault and Baogunao-Xiaohe Fault are nearly vertical faults, which is consistent with the strike-slip stress regime in this area. In order to investigate the relative magnitude of the shear and normal stress on the two faults generated by the regional stress field, 128 focal mechanisms around the 2014 Ludian earthquake were collected, and types of focal mechanism are classified. It is found that most of the focal mechanisms in this area are strike-slip type, accounting for 61.72%, while the sum of other types of focal mechanisms is only 38.28%. By using the collected focal mechanism data, we estimated the stress tensor in this region. The determined stress tensor shows NW-SE compression and NE-SW extension, also a strike-slip stress regime. It reflects that the Qinghai-Tibet plateau was pushed and uplifted by the Indian plate, and the material flowed eastward, which was blocked by the hard blocks of the Sichuan basin and forced to turn southward, showing the extrusion of NW-SE and the tension of NE-SW. By projecting the stress field onto the two faults, it was observed that the relative shear stress on both the Baoguanya-Xiaohe Fault and the Zhaotong-Ludian Fault are large with 0.99 and 0.77, respectively, indicating that both faults are capable of seismic activity. Moreover, the Baogunao-Xiaohe Fault experiences higher relative shear stress compared to the Zhaotong-Ludian Fault, suggesting that seismic activity on the former is stronger. The observed aftershock seismicity on the Baogunao-Xiaohe Fault is stronger than that on the Zhoatong-Ludian Fault illustrated by frequency, magnitude and overall releasing seismic moment v.s. time, which validates the estimated larger shear stress on the Baogunao-Xiaohe Fault.

    Table and Figures | Reference | Related Articles | Metrics
    STUDY ON THE MAGMATIC PROCESSES OF POST-COLLISIONAL POTASSIC VOLCANIC ROCKS FROM WEST KUNLUN: TAKING THE PULU AND KANGXIWA VOLCANIC ROCKS AS EXAMPLES
    DING Ran, LUAN Peng, YU Hong-mei, WEI Fei-xiang, ZHAO Bo, YANG Wen-jian, XU Jian-dong
    SEISMOLOGY AND GEOLOGY    2024, 46 (2): 312-335.   DOI: 10.3969/j.issn.0253-4967.2024.02.005
    Abstract148)   HTML17)    PDF(pc) (16484KB)(113)       Save

    The West Kunlun region is located in the northwest margin of the Qinghai-Xizang Plateau. Due to the subduction and collision of the Indian plate, this region has many post-collisional potassic volcanic areas of different sizes. Scholars have conducted many volcanic geology and petro-geochemistry work in the West Kunlun volcanic area, mainly focusing on the origin of deep magmas and plate dynamics. However, the magmatic processes of these potassic volcanic rocks are still unclear. To reveal the magmatic processes in the West Kunlun region and understand the mechanism of volcanic eruption, we analyzed the whole-rock major elements, structure and compositions of the phenocrysts, the crystal size distribution(CSD), and magma crystallization temperature and pressure conditions of the Pulu and Kangxiwa volcanic rocks. The results show that the magma sources of the two volcanic regions are close and their trace element characteristics are similar. Still, their rock types and mineral compositions are significantly different. The Pulu volcanic rocks are mainly trachyandesite and basaltic trachyandesite. The phenocrysts are composed of plagioclase, olivine, clinopyroxene and a small amount of orthopyroxene. The Kangxiwa volcanic rocks are mainly phonotephrite, consisting of clinopyroxene, biotite, and a small amount of olivine and plagioclase. The erosion and zonation of plagioclase, olivine and clinopyroxene were observed under a microscope. There are Nb-Ta and Ti negative anomalies in the two regions, with relative enrichment of large ion lithophile elements(LILF)and light rare earth elements(LREE), indicating that the magma source area has the characteristics of island arc magma, which is related to plate subduction. Based on the analysis of previous Sr-Nd datas, we suggest that the magmas from these two volcanic areas originated from enriched sources.

    According to the erosion characteristics, zonal composition data, and the concave CSD pattern, we suggest that the magma in Pulu mixed with acidic magma, whereas the magma in Kangxiwa may only mix with the internal magma, resulting in a large amount of melt erosion of phenocrysts. In Pulu volcanic rocks, the retention time of the smaller size(<5mm)crystal is 190-332a, and that of crystallographic size(>5mm)is 339-860a. The CSD curves of clinopyroxene phenocrysts in Kangxiwa phonotephrites kink at the size of 1.5mm. In Kangxiwa volcanic rocks, the residence time of smaller crystallographic size(<1.5mm)is 5.8-6.4a, and that of larger crystallographic size(>1.5mm)is 9.6-21.2a. The CSD curves of the volcanic rocks from Pulu and Kangxiwa volcanic regions are concave upward, indicating that magma mixing may have occurred both in the two volcanic regions. The An values of the core and the rim of the normal zoning plagioclases and the For value of the normal zoning olivines in the Pulu volcanic rocks vary widely, and the feldspars with low An values at the rim are out of balance with the melt. This indicates that the magma of the Pulu volcanic group mixed with the acidic magma. The Mg# of the normal zoning clinopyroxenes in the Kangxiwa volcanic rocks has a narrow range, and they are all in balance with the melt. The crystallization pressure at the rim was low, and the decompression caused a large number of resorbed phenocrysts to melt. This indicates that the mixing of phenocrysts with different degrees of melting and erosion may result in upward concave CSD curves of clinopyroxenes in the Kangxiwa volcanic rocks, so the Kangxiwa volcanic rocks may only have internal magma mixing.

    The mineral-melt equilibrium thermometers show that the equilibrium temperature and pressure of the Pulu volcanic rocks are 1 035-1 218℃, 5.1-9.9kbar, respectively, and the corresponding depth is 19.4~37.3km. The equilibrium temperature of Kangxiwa volcanic rocks is 1 154-1 282℃, the equilibrium pressure is 1.2-11.6kbar, and the corresponding depth is 4.3~43.7km. The variation range of equilibrium pressure in the Kangxiwa region is large, which may be related to the deep fault zone. In this study, by quantitatively studying the magmatic processes of post-collisional potassic volcanic rocks in the West Kunlun region, we provide CSD calculations of the volcanic rocks, reveal the migration and evolution processes of magmas in the crustal magma reservoir, and provide important information for the volcanic activities in the northwest margin of the Qinghai-Xizang Plateau and its surrounding regions.

    Table and Figures | Reference | Related Articles | Metrics
    LATE QUATERNARY DEPOSITION AND INCISION SEQUENCES OF THE BAQU RIVER AND THEIR EXPERIMENTAL IMPLICATION
    ZHANG Hao, HUANG Wei-liang, XIANG Wen, YANG Qian-hao, LIU Bo
    SEISMOLOGY AND GEOLOGY    2024, 46 (3): 570-588.   DOI: 10.3969/j.issn.0253-4967.2024.03.004
    Abstract148)   HTML18)    PDF(pc) (8597KB)(93)       Save

    River terraces are primarily formed by the erosional action of river incision under the influence of vertical movements of the crust or changes in regional erosion base levels, resulting in layered landforms. As products of the long-term evolution of river systems, the formation, development, and evolution of terraces have always been a focal point in Quaternary research. Climate change and tectonic movements play crucial roles in the evolution of river terraces, providing important evidence for studying a region’s paleoclimate and tectonic history, while also indicating the geomorphic evolution of rivers. The ages and elevations of river terraces serve as a crucial window for understanding climate fluctuations and the intensity of tectonic uplift in a specific area. This role cannot be replaced by any other method. Therefore, accurately defining the incision and deposition ages of river terraces is essential for quantitatively reconstructing the development and evolution of rivers, making it a key data point in current research on surface processes and geomorphic evolution.

    The study area is located at the southeastern margin of the Qinghai-Xizang Plateau, positioned in the main area of the Jinsha River suture zone at the southwestern edge of the Songpan-Ganzi orogenic belt and the eastern part of the Sanjiang orogenic belt. The regional tectonic setting is complex. Since the late Quaternary, the tectonic uplift at the southeastern margin of the Qinghai-Xizang Plateau has intensified, with accelerated plateau uplift in the post-Pleistocene era accompanied by significant tectonic activity. This has led to substantial incision of rivers in the region, forming multiple layers of overlapping terrace landforms on both sides of the river valleys. These landforms are crucial for quantitatively understanding the plateau uplift process and climate change.

    The Jinsha River is one of the main large rivers in the western parts of Sichuan. The river terraces developed in the Jinsha River valley serve as an important evidence for studying the deformation of the plateau crust and climate change. However, there are few Holocene terraces developed in the valley, and their resolution is low. Therefore, current research on the Jinsha River terraces mainly focuses on the orbital time scale(from tens of thousands to millions of years)of climate change and the impact of tectonic uplift, with limited studies on the role of short-term time scales(thousands or hundreds of years)in climate change and tectonic uplift, and a lack of constraints on river incision rates since the late Quaternary. The formation and evolution of river landforms since the Holocene are currently the most important means of studying recent tectonic activities and predicting future climate fluctuations. Therefore, the Baqu River, as a major tributary of the Jinsha River, with the terraces preserved in its valley, has become crucial research material reflecting the climate change and tectonic uplift in the Jinsha River Basin since the Holocene.

    The Batang segment of the Baqu River is situated in the midstream valley of the Jinsha River, characterized by a wide valley floor and gentle riverbed slope. Through drilling and shallow seismic exploration to investigate the valley stratigraphy, it was found that the valley sediments can be divided into four layers from top to bottom. The bottom layer consists of Permian strata mainly composed of weathered crystalline limestone, with a core exposure of 22m without reaching the bottom. The third layer is composed of Middle Pleistocene sediments, 68m thick, mainly consisting of large boulders, small gravel, and calcareous clay. The second layer comprises Late Pleistocene sediments, 30m thick, primarily consisting of large gravel and clay. The first layer is mainly composed of fine-grained clay with a small amount of sand and gravel blocks, 10m thick. This indicates that the valley has experienced at least two significant aggradation stages. Using Electron Spin Resonance dating methods, it was determined that these two aggradation events began at approximately 318ka and 143ka, corresponding to Marine Isotope Stages(MIS)10-9 and MIS 6-5, respectively, during glacial melting phases.

    Four levels of river terraces are developed within the valley, with T1-T3 being aggradational terraces and T4 being a bedrock terrace. T1 has a terrace height of 5~10m, T2 ranges from 15~25m, T3 ranges from 30~40m, and T4 has a terrace height of 120m. The terrace topography is generally parallel to the longitudinal profile of the modern riverbed, with only minor fluctuations, indicating a predominant overall uplift in the area after terrace formation, with consistent tectonic uplift rates and insignificant differential uplift. Combining Optically Stimulated Luminescence dating, Carbon-14 dating, and cosmogenic nuclide dating methods, it was determined that T1-T3 formed between 1~5ka, specifically 1~2ka, (3.1±0.2)ka, and(4.5±0.4)ka, respectively, while T4 formed around 62ka. A comparison of terrace ages with paleoclimate data revealed that the incision times of T1-T3 corresponded to transitions from cold to warm climates. Calculating the incision rates of terraces based on their ages and terrace heights and comparing them with incision rates in different sections of the Jinsha River, it was found that from the Late Pleistocene to the mid-Holocene, the Baqu River incision rate was(1.5±0.3)mm/a, consistent with other sections of the Jinsha River in western Sichuan. From the mid-Holocene to the present, the incision rate increased to(5.5±0.8)mm/a, approximately four times the incision rate during the Late Pleistocene. While there is a lack of quantified results on river incision rates since the Holocene in surrounding rivers, the enhanced incision rate aligns with the current vertical crustal deformation rates, indicating that intensified crustal uplift since the Holocene may be the primary driver of rapid river incision.

    Table and Figures | Reference | Related Articles | Metrics