Fault Zone Research Articles

High levels of terrestrial gamma radiation exposure in the Kopili Fault Zone on the eastern wedge of the Shillong Plateau, India

High levels of terrestrial gamma radiation exposure in the Kopili Fault Zone on the eastern wedge of the Shillong Plateau, India

Does the direct effect of friction increase continuously with absolute temperature?

Constitutive models of fault friction form the basis of physics-based simulations of seismic activity. A generally accepted framework for the slip-rate and state dependence of friction involves a thermally activated process, whereby the probability of slip along microasperities adheres to an Arrhenius law. This model, which has become widely adopted among experimentalists and theoreticians, predicts a continuous increase of the direct effect with absolute temperature, but is it observed experimentally? Leveraging comprehensive laboratory data across diverse hydrothermal, barometric, and lithological conditions, we demonstrate that, contrary to the classical view, the direct effect for a given deformation mechanism remains largely temperature-independent. Instead, the incremental shifts in the direct effect often coincide with the brittle to semi-brittle transition, across which distinct deformation mechanisms operate. These considerations challenge the validity of the classical model. Realistic constitutive laws for rock failure within the lithosphere must incorporate the contributions of multiple deformation mechanisms within active fault zones.

The Benkar Fault Zone: An Orogen-Scale Cross Fault in the Eastern Nepal Himalaya

Abstract The Benkar Fault Zone (BFZ) is a recently recognized, NNE-striking, brittle to ductile, cross fault that cuts across the dominant metamorphic fabric of the Greater Himalayan Sequence (GHS) and the Lesser Himalayan Sequence (LHS) in eastern Nepal. 40Ar/39Ar-muscovite cooling ages along a transect across the BFZ in the GHS indicate movement younger than 12 Ma. To understand the mode of genesis, and seismo-tectonic implications of the BFZ, we mapped this fault from the Everest region in the upper Khumbu valley toward the south, across the Main Central Thrust, into the LHS and the Greater Himalayan Nappe. We recognize a series of cross faults segments, which we interpret the BFZ system. The currently mapped section of the BFZ is >100 km long, and its width is up to 4 km in the LHS. The BFZ is semi-ductile in the GHS region but is brittle in the south, where it is expressed as gouge zones, tectonically brecciated zones, sharp fault planes, and segments of nonpenetrative brittle deformation zones. From petrographic and kinematic analysis, we interpret largely a right-lateral, extensional sense of shear. Our work did not continue into the Sub-Himalaya, but the BFZ may continue through this zone into the foreland as documented in other Himalayan cross faults. While several genetic models have been proposed for cross faults in the Himalaya and other convergent orogens, we suggest that the BFZ may be related to extensional structures in Tibet. Understanding cross faults is not only important for the tectonic history of the Himalaya but due to the co-location of cross faults and seismogenic boundaries, there may be a causal relationship. Cross faults also follow many of the north-south river segments of the Himalaya and weakened fault rocks on the valley walls may enhance the landslide hazard in these areas.

Hazard zonation for potential earthquake-induced landslide in the eastern East Kunlun fault zone

Abstract Based on probabilistic seismic hazard analysis, the seismic landslide hazard research considers the spatial and temporal distribution characteristics of seismic peak ground acceleration, which integrates the factors such as seismic intensity, location, and recurrence time. The occurrence of future earthquakes has certain randomness. This article presents the landslide hazard zoning of the eastern Kunlun fault zone and its surrounding faults, which is carried out under the action of horizontal ground motion with certain exceeding probability. According to the geological structure and seismicity characteristics of the study area, the potential source is divided. Based on the seismic hazard analysis and Newmark cumulative displacement evaluation model, the seismic landslide hazard in the study area is analyzed. The landslide probability is taken as the risk index. The seismic landslide hazard can be divided into five grades: extremely low-prone area, low-prone area, medium-prone area, high-prone area, and extremely high-prone area. In the results of seismic landslide risk zoning given in this article, the surrounding areas of Tazang fault and Minjiang fault are high-risk areas, which should be paid attention to.

Quantifying Creep on the Laohushan Fault Using Dense Continuous GNSS

The interseismic behavior of faults (whether they are locked or creeping) and their quantitative kinematic constraints are critical for assessing the seismic hazards of faults and their surrounding areas. Currently, the creep of the eastern segment of the Laohushan Fault in the Haiyuan Fault Zone at the northeastern margin of the Tibetan Plateau, as revealed by InSAR observations, lacks confirmation from other observational methods, particularly high-precision GNSS studies. In this study, we utilized nearly seven years of observation data from a dense GNSS continuous monitoring profile (with a minimum station spacing of 2 km) that crosses the eastern segment of the Laohushan Fault. This dataset was integrated with GNSS data from regional continuous stations, such as those from the Crustal Movement Observation Network of China, and multiple campaign measurements to calculate GNSS baseline change time series across the Laohushan Fault and to obtain a high spatial resolution horizontal crustal velocity field for the region. A comprehensive analysis of this primary dataset indicates that the Laohushan Fault is currently experiencing left-lateral creep, characterized by a partially locked shallow segment and a deeper locked segment. The fault creep is predominantly concentrated in the shallow crustal region, within a depth range of 0–5.7 ± 3.4 km, exhibiting a creep rate of 1.5 ± 0.7 mm/yr. Conversely, at depths of 5.7 ± 3.4 km to 16.8 ± 4.2 km, the fault remains locked, with a loading rate of 3.9 ± 1.1 mm/yr. The shallow creep is primarily confined within 3 km on either side of the fault. Over the nearly seven-year observation period, the creep movement within approximately 5 km of the fault’s near field has shown no significant time-dependent variation, instead demonstrating a steady-state behavior. This steady-state creep appears unaffected by postseismic effects from historical large earthquakes in the adjacent region, although the deeper (far-field) tectonic deformation of the Laohushan Fault may have been influenced by the postseismic effects of the 1920 Haiyuan M8.5 earthquake.

Сейсмотектоническая активность Ляховичского разлома в Припятском прогибе

The main result of the work was the identification of the present-day seismotectonic regime of an active zone of the Lyakhovichi regional fault. The areal and space-time distribution of the earthquake epicenters along the fault line was reconstructed. The coordinates of the seismically active fault zone located within an area of 52.738—53.037° N and 26.830—28.034° E were determined. It was found that the distribution of earthquakes along the fault line occurred non-uniformly in space and time. The period from 2012 till 2021 was the most active one by the quantity of occurred earth quakes (169) and released total seismic energy (4.5929•109 J). The strongest earthquake with a magnitude M = 2.9 was recorded in the eastern part of the fault on September 03, 2021. Our reconstructions determined the value of seismic dislocation along the fault, which corresponded to the kinematic slip fault. To characterize the modern seismogeodynamic activity in the Lyakhovichi fault zone, a comprehensive technique based on long-term monitoring of seismic events was used. Based on the performed studies, the seismotectonic potential of the seismoactive zone of the Lyakhovichi fault was determined.

Use of three-dimensional crossgradient joint inversion of marine controlled-source electromagnetic (CSEM) and seismic refraction data for overburden geohazards mapping.

Offshore energy-transition activities typically require the construction of new subsurface facilities where overburden heterogeneity and deformation pattern may not be fully understood. Anisotropic marine controlled-source electromagnetic (CSEM) and seismic refraction imaging can be adapted to geohazards assessment of the overburden to mitigate safety and costs risks in such areas. This is investigated using an area in deepwater NW Borneo as the type-example. For this proof-of-concept study, a well-established Multiphysics joint inversion code was extended to handle marine seismic refraction data from conventional towed streamers. Realistic field-scale synthetic models with anomalous resistivities and velocities in the shallow section (∼50-200 m below mudline) were built and their CSEM and seismic refraction responses were first inverted separately and then jointly in 3D using the crossgradient method. The joint inversion recovered the anomalous bodies better than the separate inversions. The velocity model from the joint inversion showed improvement from the separate inversion but required a good starting model. As a further test of the methodology, first-break travel times were picked from legacy 3D seismic data and inverted with the available sparse CSEM data. For model verification, resistivity and velocity profiles from separate and joint inversion models were extracted at a well location in the survey area and compared with actual well logs. The joint inversion models matched the well logs better than the separate inversion models. Moreover, known geohazards (complex fault-sets and lateral-inflows) from independent geomatic study were imaged in the joint inversion models. Also, the resistivity anisotropy ratio tracked anomalous zones that correspond to zones of recent complex faulting and potential lateral inflows in the overburden mass transport systems. Thus, joint inversion of CSEM and seismic refraction data has potential to improve the mapping of overburden hazards offshore, with the associated resistivity anisotropy serving as a proxy for active ground deformation.

Co- and postseismic subaquatic evidence for prehistoric fault activity near Coyhaique, Aysén Region, Chile

Abstract. Chilean Patagonia is confronted with several geohazards due to its tectonic setting, i.e., the presence of a subduction zone and numerous fault zones, e.g., the Liquiñe-Ofqui Fault Zone (LOFZ). This region has therefore been the subject of numerous paleoseismological studies. However, this study reveals that the seismic hazard is not limited to these large tectonic structures by identifying past fault activity near Coyhaique in Aysén Region. Mass-wasting deposits in Lago Pollux, a lake located ca. 15 km SW of this region's capital, were identified through analysis of reflection-seismic data and were linked to a simultaneous event recorded in nearby Lago Castor. Furthermore, a coeval ∼50-year-long catchment response was identified in Aysén Fjord based on the multiproxy analysis of a portion of a sediment core. Assuming that this widely recognized event was triggered by an earthquake, ground-motion modeling was applied to derive the most likely magnitude and source fault. The model showed that an earthquake rupture along a local fault, in the vicinity of Lago Pollux and Lago Castor, with a magnitude of 5.6–6.8, is the most likely scenario.

Study on the Susceptibility of Steel Arches with Letting Pressure Nodes Based on Incremental Dynamic Analysis

When soft rock tunnels pass through fractured fault zones, they are particularly susceptible to extrusion and large-scale deformations, especially during seismic events. To address these challenges, this study introduces an innovative yield-support steel arch design featuring a circumferential letting pressure node at its core. This design delivers incremental support resistance within the deformation zone and a susceptibility curve is applied to evaluate the damage probability of the steel arch with a letting pressure node under seismic loading conditions. Measurements of the surrounding rock pressure and structural forces on the steel arch with the letting pressure node were conducted at the Baoshan Jewel Mountain Tunnel in China. The field experiment results revealed a 23% reduction in the surrounding rock pressure and an 11% decrease in the internal forces of the support structure. These findings demonstrate the successful application of the letting pressure node-supported steel arch in mitigating large deformations in soft rock environments. Additionally, using finite element software ANSYS 2022, a seismic time-history analysis was conducted, employing the relative deformation rate of the letting pressure node steel arch as the damage index and the peak ground acceleration (PGA) as the strength parameter to generate the incremental dynamic analysis (IDA) curve. According to the susceptibility curve derived from the incremental dynamic analysis, at the design ground motion level of 8 degrees, the letting pressure node steel arch has a 94% probability of exceeding its normal service life limit and experiencing damage. The findings of this study offer a novel approach to addressing large deformations in soft rock tunnels. The proposed susceptibility curves for steel arches with letting pressure nodes provide a robust foundation for predicting the damage probability of yielding support structures under seismic conditions.

Crustal and uppermost mantle S-velocity structure of the Seoul metropolitan area on the Korean Peninsula from Helmholtz tomography

The Seoul metropolitan area, the most densely populated part of the Korean Peninsula, features complex subsurface structures and seismogenic faults, though their characteristics remain ambiguous due to low seismicity and limitations in fault investigation. High-resolution velocity models can provide constraints for identifying subsurface faults by detecting elongated low-velocity anomalies along fault zones. Recently, a dense seismic network was deployed in this area, facilitating the use of Helmholtz tomography, an array-based method that accounts for finite-frequency effects. Utilizing Helmholtz tomography, we obtained a high-resolution S-wave velocity model down to a depth of 50 km with waveform data recorded at 74 broadband seismic stations. We found that a linear low-velocity anomaly along the Pocheon fault extends to the uppermost mantle, with an increasing width with depth. In contrast, the Dongducheon fault, which traverses Seoul from north to south, is not well imaged, indicating its current weak activity. Another linear low-velocity anomaly extends southwest through Seoul from northern Seoul, potentially representing the extension of the Pocheon fault based on similar strike and dip directions. Additionally, a large lateral low-velocity anomaly is identified in the lower crust beneath the northern part of the Seoul metropolitan area, interpreted as a ductile décollement, connected with the Pocheon, Wangsukcheon, and possibly Gyeonggang faults. This study successfully identified the extensions and orientations of subsurface faults beneath the Seoul metropolitan area down to the uppermost mantle, which is critical for seismic hazard predictions and earthquake simulations in this highly populated area.

3D structural modelling of the Kopili fault zone in North-East India: Seismotectonic analysis utilising focal mechanism solutions of small-to-moderate earthquakes

This work presents new findings from recent seismic data which have never been exploited so far for the weak-to-moderately energetic (ML 2.3–6.4) earthquakes and to depict a 3D scheme for the Kopili Fault (KF) Zone. For the first time, we figure out different fault kinematics in the upper crust such as reverse and strike-slip mechanisms located along the KF zone. Moment Tensor (MT) solutions are successfully obtained through full-waveform inversions for 60 seismic events. The stress inversions are estimated based on the MT results as well as previous findings from 16 small to moderate earthquakes (MW 3.6–6.3). The MT solutions demonstrate that strike-slip kinematics is the most prevalent while a few thrust/normal earthquake events are observed in the area. The seismic cross-sections display that these earthquakes took place at a depth of 40 km and are mostly generated on the eastern side of the fault. According to the 3D spatial and temporal distribution of relocated events, the KF zone is not a completely vertical fault but rather one with a slight eastward dip. Our stress inversion results also indicate that this fault is equivalent to a fault plane with a dip angle of around 80˚, which is compatible with the previously estimated dip angle (75˚). We interpret that the research region is under a transpressional stress regime, which is likely to be the reason for the current activity of the KF zone. We postulate that our findings can be useful to enhance seismic hazard assessment in the densely populated region, which is also covered by vital infrastructure.

Two-dimensional numerical analysis of seepage in an undersea mining area considering hydrochemical corrosion

The complex hydrogeological conditions and unique hydrochemical environment in the seabed rock often cause serious mine water inrush during undersea mining. With the exploitation of the undersea gold mining area in Sanshandao as the research background, a numerical analysis method for the seepage behavior of the surrounding rock in undersea gold mining considering the hydrochemical corrosion was put forward. A numerical simulation analysis of the undersea mining area in Sanshandao was conducted by the established method, and the results showed that mining the submarine ore body significantly changed the seepage field of the surrounding rock around the stope and water levels dropped obviously in the ore-controlling fault zone. Differences in the hydrochemical environment significantly influenced the distribution of the seepage field of the surrounding rock. The calculation results based on three fluids with different pH values show that the head variation range is the largest when pH = 2, followed by pH = 4, and pH = 7 has the smallest head variation range. The water inflow in the stope was also directly affected by the hydrochemical environment. With the change of pH, the water inflow at each boundary of the stope increases exponentially.

On the Billefjorden fault zone in Garmdalen, central Spitsbergen: implications for the mapping of major fault zones during geological fieldwork and for the tectonic history of Svalbard.

The present contribution reexamines the geometry of a segment of a presumably long-lived fault in Svalbard, the Balliolbreen Fault segment of the Billefjorden Fault Zone, along which presumably two basement terranes of Svalbard accreted in the early-mid Paleozoic after thousands of kilometers strike-slip displacement. We performed structural fieldwork to Billefjorden in central Spitsbergen and interpreted satellite images. Field observations demonstrate that the Balliolbreen Fault formed as a top-west thrust fault in the early Cenozoic and that weak sedimentary units such as shales of the Lower Devonian Wood Bay Formation and coals of the uppermost Devonian-Mississippian Billefjorden Group partitioned deformation, resulting in significant contrast in deformation intensity between stratigraphic units. For example, tight early Cenozoic folds are localized in shales of the Wood Bay Formation and contemporaneous top-west brittle-ductile thrusts within coals of the Billefjorden Group, whereas Pennsylvanian deposits of the Hultberget (and/or Ebbadalen?) Formation are simply folded into gentle open folds. Rheological contrasts also resulted in the development of décollements locally, e.g., between tightly folded strata of the Wood Bay Formation and Billefjorden Group and flat-lying, brecciated limestone-dominated strata of the Wordiekammen Formation. Despite the limited quality and continuity of outcrops in the area, the eastward-thickening character (i.e., away from the fault) of Pennsylvanian deposits of the Hultberget, Ebbadalen, and Minkinfjellet formations suggests that the fault did not act as a normal fault in Pennsylvanian times. The study suggests that strain partitioning of early Cenozoic Eurekan contraction alone may explain the deformation patterns in Paleozoic rock units in central Spitsbergen, i.e., that Late Devonian Svalbardian contraction is not required, and that a major segment of the Billefjorden Fault Zone formed in the early Cenozoic. The present work illustrates the crucial need for interdisciplinary approaches and composite educational backgrounds in science.

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

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

Very‐Near‐Field Coseismic Fault Pressure Drop and Delayed Postseismic Cross‐Fault Flow Induced by Fault Damage From the 2018 M6.3 Hualien, Taiwan Earthquake

AbstractEarthquakes can produce rock damage, poroelastic deformation, and ground shaking that modify fault zone hydrogeologic properties. Coseismic and postseismic hydrologic response to the ruptured fault can serve as constraints on hydrogeologic property changes. Here, we document fluid pressure responses to the 2018 M6.3 Hualien, Taiwan earthquake and model the postseismic fault zone hydrology inferred from very‐near‐fault data. Two groundwater wells located ∼180–250 m from the fault damage zone experienced 10–15 m of groundwater level decline followed by a prolonged (>6 months) recovery. Poroelastic models indicate that strain dilation in the fault's hanging wall can produce water level reductions of several meters. However, the model cannot explain groundwater level change in the footwall nor the delayed postseismic recovery, suggesting fault zone damage played a primary role in both the coseismic and postseismic water level reductions. Fluid flow models incorporating the effects of coseismic fault damage on the temporal evolution of hydrologic fault properties show enhanced hydraulic anisotropy after the earthquake. The best‐fit models show that while both vertical hydraulic conductivity and specific storage likely increased within the fault zone, the horizontal hydraulic conductivity is low, suggesting the fault zone behaves as a hydraulic barrier to cross‐fault flow with modeled diffusivities as low as ∼10−3 m2/s. High‐fidelity hydrologic measurements in the very‐near‐field of fault zones (e.g., within a few hundred meters) may be a useful tool to constrain the spatial and temporal evolution of fault zone properties before, during, and after rupture.

Numerical simulation of seismic stress evolution in the Lajishan Fault zone on the northeastern margin of the Qinghai-Xizang Plateau

To analyze the triggering relationship of historical earthquakes to the 2023 Jishishan earthquake and estimate the seismic risk in the adjacent region after the earthquake, we calculate the evolution process of seismic stress generated by 11 earthquakes with MS ≥ 6.0 and 1 earthquake with MS5.7 near the source in the past 100 years. The results show that before the 2023 Jishishan earthquake, the entire Lajishan Northern Edge Fault zone was under the stress shadow generated by the 12 earthquakes. The stress reduction at the 2023 Jishishan earthquake epicenter was approximately 15.1 kPa, and the contribution of this stress shadow mainly comes from the 1920 MS8.5 Haiyuan earthquake. The coseismic stress drop at the epicenter was around −292.1 kPa, and the maximum stress drop for the 2023 Jishishan earthquake was around −346.2 kPa. The Jishishan earthquake increased the stress of the two seismic gaps in the Lajishan Northern Edge Fault by 14.7 kPa and 59.7 kPa respectively, the stress on the Lajishan Southern Edge Fault zone increased by 10.9 kPa, and the stress on the east section of Xunhua Nanshan Fault increased by 11.1 kPa. These stress increments potentially elevate the seismic hazard along these fault zones. In the future earthquake prevention and disaster reduction, attention should be paid to the seismic risk of these fault zones.

Distribution rules of remaining oil by bottom water flooding and potential exploitation strategy in fault- controlled fractured-vuggy reservoirs

Based on the tectonic genesis and seismic data of fault-controlled fractured-vuggy reservoirs, the typical fractured-vuggy structure features were analyzed. A 3D large-scale visual physical model of “tree-like” fractured-vuggy structure was designed and made. The experiments of bottom-water flooding and multi-media synergistic oil displacement after bottom-water flooding were conducted with different production rates and different well-reservoir configuration relationships. The formation mechanisms and distribution rules of residual oil during bottom-water flooding under such fractured-vuggy structure were revealed. The producing characteristics of residual oil under different production methods after bottom-water flooding were discovered. The results show that the remaining oil in “tree-like” fractured-vuggy structure after bottom-water flooding mainly include the remaining oil of non-well controlled fault zones and the attic remaining oil at the top of well controlled fault zones. There exists obvious water channeling of bottom-water along the fault at high production rate, but intermittent drainage can effectively weaken the interference effect between fault zones to inhibit water channeling. Compared with the vertical well, horizontal well can reduce the difference in flow conductivity between fault zones and show better resistance to water channeling. The closer the horizontal well locates to the upper part of the “canopy”, the higher the oil recovery is at the bottom-water flooding stage. However, comprehensive consideration of the bottom-water flooding and subsequent gas injection development, the total recovery is higher when the horizontal well locates in the middle part of the “canopy” and drills through a large number of fault zones. After bottom water flooding, the effect of gas huff and puff is better than that of gas flooding, and the effect of gas huff and puff with large slug is better than that of small slug. Because such development method can effectively develop the remaining oil of non-well controlled fault zones and the attic remaining oil at the top of well controlled fault zones transversely connected with oil wells, thus greatly improving the oil recovery.

Fault zone architecture and hydraulic properties characterization using earth tide analysis

Fault zone architecture and hydraulic properties, which determine subsurface fluid flow and storage, have received long-term research attention from hydrogeologists. Current field experimental techniques and numerical simulation methods are limited by the scale and complexity of the investigations and are inadequate for deep subsurface fault zones. Estimating the internal rock properties and structures of fracture zones requires a combination of methods from many different fields. In recent years, groundwater tidal response methods have been widely used to characterize subsurface hydrogeomechanical properties because of their advantages of being in situ and low cost. In this study, tidal analysis was employed for a well within the Blue Ridge Fault Zone, and we integrated the fracture zone properties revealed by a single well. Well W-03 revealed an aquifer with an overall permeability of ∼ 10–14 m2 and hydraulic diffusivity of 0.24 m2⋅s−1. The fault damage zone had a permeability of 10–14 ∼ 10–13 m2 while the host rock had a permeability of 10–16 ∼ 10–14 m2. In addition, we found that the heterogeneity of specific storage affects the assessment of poroelasticity parameters, we offer a method to estimate these parameters without relying on specific storage, yielding porosity estimates ranging from 0.05 to 0.06, aligning with field test results. Furthermore, the estimated dip and strike of the fault zone were also consistent with the outcrop and borehole observations.

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

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

Tectonic-Paleoseismological Characteristics and Quaternary Activity of Maymundağı Fault (Northern Acıgöl Graben)

The Aegean region and its graben system constitute one of Turkey’s most significant seismic zones. The faults within the Aegean graben generate numerous earthquakes, leading to various human and economic losses. To better understand the seismicity of western Anatolia, it is necessary to obtain concrete findings regarding the seismic history of earthquake-producing graben faults. This can be achieved through paleoseismological studies and other relevant disciplines. This study focuses on paleoseismological investigations along the northern boundary fault of the Acıgöl graben, located east of the Aegean graben system. The Maymundağı fault zone has been examined in two separate segments: east and west. The Dazkırı segment to the east shows evidence of movement dating back at least 10,000 years, with subsequent intensified activity observed later on the western Bozkurt segment. An earthquake occurred approximately 2370 years ago east of the Bozkurt segment, followed by movements migrating westward, resulting in earthquakes approximately 1322 and 598 years ago. Further analysis of the western segment indicates an average recurrence interval of 724 years for earthquakes, with a slip rate of 0.58 mm/year. Based on these findings, a future earthquake can be expected in this region around 2028–2129 AD.