Strong Earthquakes Research Articles

Seismic moment of local earthquakes in the central part of the Baikal rift by the coda envelope inversion

The article presents an estimation of scalar seismic moment by the coda envelope inversion of local earthquakes recordings in the central part of the Baikal rift zone. In addition to the earthquake source parameters, the method allows us to simultaneously estimate seismic energy loss and site amplification factors for frequency bands from 0.53 to 34 Hz. Because of the compromise between spectral source energy and site responses, we estimated amplification factors for used stations relative to the "UlanUde" (UUDB) reference site located in crystalline rock. In order to estimate intrinsic and scattering attenuation we used events recorded at three broadband seismic stations and epicentral distances between 40 and 120 km. The result suggests that intrinsic absorption is dominant over scattering attenuation in the central part of the Baikal rift zone for most frequency bands, but the seismic albedo B0, expressing the contribution of scattering to the total attenuation, showed variations from 0.23 to 0.6 for frequencies below 1 Hz, with mean B0 value 0.33. Attenuation and siteamplification factors were used to estimate the seismic moment and moment magnitudes Mw of local earthquakes. The resultant moment magnitudes exhibit a good agreement with routinely reported local magnitude (ML) estimates for the study area. The coda inversion estimates of seismic moment provide stable, unbiased moment magnitudes for events that are too small to be seen at teleseismic distances.

Detection of Small Earthquakes by Waveform Envelope Using Machine Learning

ABSTRACT Analyzing seismic data has helped reveal most of our knowledge about the Earth’s interior. However, the number of moderate-to-large earthquakes is limited; and most small earthquakes, although essential for monitoring the dynamic process between major earthquakes, are difficult to detect from earthquake recordings. In this study, we present a detection method that catches the envelope pattern of seismic data in machine learning. The waveform envelope presents reliable features and a low-frequency pattern for any event rather than being dominated by details of noise and signals. We first apply a fully connected neural network to extract data envelopes from recordings and then a convolutional neural network to detect events using the envelopes as input. Our method is tested against the previously published method on seismic data in Japan. The new approach identifies the largest number of small events. In the testing dataset, the precision and recall of the approach for events are 98.73% and 96.54%, respectively, and those for noise are 97.41% and 99.05%, respectively. We demonstrate that the approach performs well in different signal-to-noise ratios (SNRs) and filter frequency band tests. For events with low SNRs (0–6 dB), the detection accuracy of the approach is approximately 10% higher than a published deep learning method. Besides most of the events in the earthquake catalog, the process also detects more events from continuous data that are not documented in the catalog.

Possible Identification of Precursor ELF Signals on Recent EQs That Occurred Close to the Recording Station

The Lithospheric–Atmospheric–Ionospheric Coupling (LAIC) mechanism stands as the leading model for the prediction of seismic activities. It consists of a cascade of physical processes that are initiated days before a major earthquake. The onset is marked by the discharge of ionized gases, such as radon, through subterranean fissures that develop in the lead-up to the quake. This discharge augments the ionization at the lower atmospheric layers, instigating disturbances that extend from the Earth’s surface to the lower ionosphere. A critical component of the LAIC sequence involves the distinctive perturbations of Extremely Low Electromagnetic Frequencies (ELF) within the Schumann Resonances (SR) spectrum of 2 to 50 Hz, detectable days ahead of the seismic event. Our study examines 10 earthquakes that transpired over a span of 3.5 months—averaging nearly three quakes monthly—which concurrently generated 45 discernible potential precursor seismic signals. Notably, each earthquake originated in Southern Greece, within a radius of 30 to 250 km from the observatory on Mount Parnon. Our research seeks to resolve two important issues. The first concerns the association between specific ELF signals and individual earthquakes—a question of significant importance in seismogenic regions like Greece, where earthquakes occur frequently. The second inquiry concerns the parameters that determine the detectability of an earthquake by a given station, including the requisite proximity and magnitude. Initial findings suggest that SR signals can be reliably linked to a particular earthquake if the observatory is situated within the earthquake’s preparatory zone. Conversely, outside this zone, the correlation becomes indeterminate. Additionally, we observe a differentiation in SR signals based on whether the earthquake took place over land or offshore. The latter category exhibits unique signal behaviors, potentially attributable to the water layers above the epicenter acting as a barrier to the ascending gases, thereby affecting the atmospheric–ionospheric ionization process.

Study on the Effect of Liquefiable Overburden Foundations of Rockfill Dams Based on a Pore Pressure Model

China’s southwestern region boasts abundant hydropower resources. However, the area is prone to frequent strong earthquakes. The areas surrounding dam sites typically have deep overburden, and the liquefaction of saturated sand foundations by earthquakes poses significant safety risks to the construction of high dams in the southwest. The effects of liquefaction and reinforcing measures on the foundations of rockfill dams on liquefiable overburden under seismic action are currently the subject of somewhat unsystematic investigations. The paper utilizes the total stress and effective stress methods, based on the equivalent linear model, to perform numerical simulations on the overburden foundations of rockfill dams. The study explores how factors such as dam height, overburden thickness, liquefiable layer depth, liquefiable layer thickness, ground motion intensity, and seismic wave characteristics affect the liquefaction of the overburden foundations. Additionally, it examines how rockfill dams impact the dynamic response, considering the liquefaction effects in the overburden. The results show that although the total stress method, which ignores the cumulative evolution of pore pressure during liquefaction, can reveal the basic response trend of the dam, its results in predicting the acceleration response are significantly biased compared to those of the effective stress method, which comprehensively considers the cumulative changes in liquefaction pore pressure. Specifically, when the effect of soil liquefaction is considered, the predicted acceleration response is reduced compared to that when liquefaction is not considered, with the reduction ranging from 4% to 30%; with increases in the thickness and burial depth of the liquefiable layer, the effective stress method considering liquefaction significantly reduces the predicted peak acceleration; the effect of liquefiable soil on the attenuation of the speed response is more sensitive to the low-frequency portion of the seismic wave. The study’s findings are a significant source of reference for the planning and building of rockfill dams on liquefiable overburden.

Relative Moment Tensor Inversion for Microseismicity: Application to Clustered Earthquakes in the Cascadia Forearc

The relative abundance of small earthquakes affords significant opportunities for improved understanding of regional seismotectonics; however, determining moment tensors for such events recorded on regional networks is complicated by low signal-to-noise ratios, sparse station sampling and complex wave propagation at short periods. We build upon previous work in designing a multiple-event, simultaneous moment tensor inversion scheme for small earthquakes that employs constraints from P-wave polarities, relative amplitudes of P- and S-waves recorded at common stations, and local magnitude estimates. Our method does not require a priori knowledge of a reference moment tensor. High-fidelity polarity and relative amplitude data are recovered using principal component decomposition of clustered-event waveforms. These data are employed within a multi-stage iterative framework to invert for moment tensors and incorporate local magnitude information. Synthetic examples employing as few as four high-quality and spatially-distributed stations yield accurate moment tensor estimates. We demonstrate our approach on a cluster of seismicity near San Juan Island, Washington, USA, within the Cascadia forearc. Our results are consistent with previous characterization of the local stress regime, and support an interpretation of swarm behaviour resulting from migration of fluids originating from dehydration of the subducting Juan de Fuca plate.

Review of Earthquake Resistance and Mitigation of Underground Structures

The rapid development of underground structural engineering in China has been accompanied by pressing issues in the research of earthquake resistance and mitigation. Typically, underground structures exhibit good seismic performance with relatively few disasters. However, once they are damaged by earthquakes, the consequences can be severe and difficult to repair. This paper first elaborates on the prototype observation, theoretical analysis, model testing, and numerical simulation involved in seismic response analysis of underground structures. It comprehensively analyzes the advantages, disadvantages, and applicability of different research methods. Finally, it summarizes the research achievements on technical measures to mitigate earthquake damage to underground structures in strong earthquake regions, both domestically and internationally.

Increase in short-term and long-term stress-associated mental illness after Jiji earthquake in Taiwan: A twenty-year longitudinal population-based cohort study (1999–2019)

BackgroundEarthquakes have caused profound physical and mental health impacts in human history. The Jiji earthquake, which had a magnitude of 7.6 on the Richter scale, occurred on 21 September 1999 in Taiwan. A close follow-up on the mental health status of affected adults after major natural disasters to construct the short-term and long-term risk and prevalence of stress-associated mental illnesses has not been performed by using the nationwide health databases. MethodsThis population-based cohort study included 468,804 adults affected by Jiji earthquake spanning from 2000 to 2019 who were matched at a 1:4 ratio with unaffected individuals based on age and sex (n = 1,875,216). Employing a subdistribution hazard regression analysis, we assessed the incidence of sleep, anxiety, and depressive disorders after Jiji earthquake. Corrections for multiple comparisons were carried out using the Benjamini–Hochberg procedure. ResultsAffected adults experienced an increased incidence of short-term (approximately twice) stress-associated psychiatric disorders. The risk of the post-traumatic stress disorder (PTSD) is significantly higher in the affected adults (40–64 years: aSHR: 92.0; ≥65 years: aSHR: 96.7, p < 0.0001). Middle-aged (aged 40–64 years) male adults presented with significantly more short-term (< one year) and long-term (up to 20 years) stress-related mental illnesses, i.e., insomnia, anxiety, and depressive disorders, than individuals in the control group. ConclusionsAn earthquake has significant short and long-term effects on sleep quality, anxiety, and depressive disorders in affected adults. Optimal short and long-term close monitoring is needed to deploy medical resources and socioeconomic support to relieve mental stress burdens.

Seismic Upgrading of the Heritage-Protected Reinforced Concrete Warehouse in Rijeka, Croatia

Despite Croatia experiencing two strong earthquakes in 2020, Rijeka was not directly affected, underscoring the importance of proactive seismic assessment and strengthening in all seismic regions. This paper presents a comprehensive case study on the seismic strengthening of a 20th-century concrete building located in Rijeka, Croatia, originally designed according to Austro-Hungarian construction norms and practices. As a heritage-protected structure, the building’s architectural features and construction practices were examined and contextualized within its historical background. The assessment and renovation phases of this project are discussed in detail, demonstrating the practical application of modern seismic strengthening techniques while preserving the building’s historical integrity. This case study aims to highlight the need for such measures to protect heritage structures and to show the implementation of rapid and new (ad hoc) norms for earthquake-damaged buildings in Croatia. This study serves as a reference for engineers, architects, and conservationists involved in the preservation of heritage buildings, demonstrating that it is possible to enhance their structural safety without compromising their architectural authenticity.

Seismic assessment of glulam frames with dual-tube self-centering buckling-restrained braces

The development of resilient lateral load-resisting systems is essential for multi-story and high-rise timber buildings in earthquake-prone areas. In this paper, an attempt is made to integrate dual-tube self-centering buckling-restrained braces (DT-SCBRBs) to glulam frames, aiming to improve their structural resilience to earthquakes. To assess the seismic effectiveness of such a DT-SCBRB glulam frame, nonlinear time-history analyses (NLTHAs) were conducted on a series of prototype DT-SCBRB glulam frames with different building heights and design parameters of DT-SCBRBs. The seismic performance of these prototype structures was evaluated in terms of maximum inter-story drift ratios (MaxISDRs), residual inter-story drift ratios (ResISDRs), and inter-story drift uniformity. The results show that the frame with lower post-tensioning-to-yielding ratios of DT-SCBRBs tended to exhibit lower MaxISDRs but could sometimes lead to higher ResISDRs under major earthquakes. Low deformability of tendons could result in tendon fracture under major earthquakes. The ResISDRs of all the prototype frames were lower than 0.5 %, the drift limit for the “Repairable” performance level. Compared to minor and moderate earthquakes, the DT-SCBRB glulam frames deformed more uniformly during major earthquakes. Incremental dynamic analyses (IDAs) were also conducted on these prototype structures, developing a database to quantify performance levels of DT-SCBRB glulam frames. The MaxISDRs of 0.5 %, 0.7 %, and 2.6 % were recommended for immediate occupancy (IO), life safety (LS), and collapse prevention (CP) limit states of DT-SCBRB glulam frames.

Travel distance estimation of landslide-induced debris flows by machine learning method in Nepal Himalaya after the Gorkha earthquake

Debris flows are more likely to be triggered in the earthquake-strike areas with a widespread presence of unstable slopes, causing severe casualties and changing the surrounding natural topography. In such scenario, estimating the travel distance of debris flows becomes crucial to understand the hazardous areas. Therefore, a hybrid machine learning model (GA-XGBoost) was employed to achieve a reliable estimation of debris-flow travel distance. This model was applied to the Nepal Himalayas, the site of the 2015 Gorkha earthquake. We selected four geomorphological factors for travel distance estimation. They are the volume of failure mass (VL), the height difference between the material source center and end point of movement mass (H), the mean gradient of the travel path (J), and the mean curvature of the travel path (C). Furthermore, to eliminate the noise information and enhance stability of input data, a principal component analysis (PCA) was used to generate three principal components (PC1, PC2, and PC3) from the selected factors to serve as input variables of model development. The performance of this model was evaluated using the assessment indexes, resulting in a mean absolute percentage error (MAPE) of 8.71%, a root mean square error (RMSE) of 144.3 m, and a mean absolute error (MAE) of 86.1 m. Four empirical approaches were also introduced for comparison analysis. Our proposed model has proven to be superior and effective, as the estimated results closely match the actual values. All the results affirm the suitability of our developed model for estimating the travel distance of landslide-induced debris flows following a strong earthquake.

Experiment and analysis on seismic performance of a self-centering Y-eccentrically braced frames structure

Conventional eccentrically braced frames (EBFs) exhibit significant residual deformations following major earthquakes, necessitating costly repairs for the damaged structures. This study proposed a novel self-centering Y-eccentrically braced frames (SC-YEBFs) system to enhance the seismic resilience. A theoretical mechanical model which takes into account its components is formulated to predict the lateral load behavior of SC-YEBFs. The pseudo-static experiments were conducted on four specimens, and the resulting observations and analyses demonstrate that the shear links effectively function as ductile fuses, dissipating a significant portion of the input energy to safeguard the main frame from damage. Furthermore, the self-centering joint exhibits a hinge-like behavior with remarkable self-centering characteristics. The structural reparability of all specimens was found to be exceptional during major seismic events. By increasing the initial prestress and sectional area of steel strands, the self-centering performance can be further enhanced. The damaged shear links could be easily detached by loosening bolts, and the consideration of prestress loss is nonnegligible to ensure an exceptional self-centering performance. The proposed theoretical model predicted the mechanical properties of SC-YEBFs, including lateral stiffness, energy dissipation, and self-centering mechanism, through a comparison with experimental results. Moreover, the theoretical model was utilized to propose an equivalent simulation method for simulating self-centering joint and shear link using Connector function in ABAQUS software, while establishing a simplified frame finite element model for further analysis.

Stone column strategies for mitigating liquefaction-induced uplift of tunnels

Tunnels embedded in liquefiable soil are frequently subjected to uplift and sustain serious damage during major earthquakes. Mitigation methods to prevent the flotation of these tunnels must be developed and implemented in the natural environment. For this purpose, this study proposes a new method that uses stone columns to enhance soil drainage and mitigate soil liquefaction around tunnels. A circular tunnel in liquefied soil is simulated using a 2D finite element model, PLAXIS 2D, and subjected to a sinusoidal input motion. The tunnel's pre-construction and post-construction scenarios are examined. Parametric studies are carried out to investigate the effect of changing several parameters, such as the distance between stone columns and tunnel springing, the diameter of stone columns, the spacing between stone columns, the number of stone column rows, and the stone column arrangement patterns on the effectiveness of liquefaction mitigation. The study reveals that liquefaction mitigation is enhanced by using stone columns closer to tunnel springing, with larger diameters, less spacing, and more rows of stone columns that are arranged in a square pattern. It also emphasizes the importance of timely implementation of stone columns for maximum benefit.

Identifying the Characteristics and Implications of Post-Earthquake Debris Flow Events Based on Multi-Source Remote Sensing Images

Strong earthquakes often bring amounts of loose material, disrupting the balance of material transportation within a watershed and severely impacting the restoration of the ecological environment and human safety downstream. Therefore, it is crucial to identify the frequency and scale of these debris flow events, as well as to explore their long-term development and impact on internal and external channels. Using multi-source remote sensing images from four perspectives, hillslope, channel, accumulation fan, and their relationship with the mainstream, we reconstructed a debris flow event dataset from 2008 to 2020, explored a method for identifying these events, and analyzed their impacts on channels and accumulation fans in Mozi Gully affected by the Wenchuan earthquake. Loose matter was predominantly found in areas proximate to the channel and at lower elevations during debris flow events. Alterations in channel width, accumulation fans downstream, and their potential to obstruct rivers proved to be vital for identifying the large scale of debris flow event. Finally, we encapsulated the evolution patterns and constraints of post-earthquake debris flows. Determination in frequency and scale could offer valuable supplementary data for scenario hypothesis parameters in post-earthquake disaster engineering prevention and control.

Estimation of seismometer clock time offsets using Kalman Filter towards accurate seismic velocity change

SUMMARY Monitoring seismic velocity changes obtained from ambient noise correlations is widely used to understand changes in rock properties in response to earthquakes, volcanic activities and environmental changes. Since continuous seismic data have been accumulated, this method can estimate long-term changes in seismic velocity, such as crustal recovery after a major earthquake and temporal variations in seismic velocity related to long-term environmental change. Changes in seismic velocity can be estimated with a high temporal resolution by measuring the phase differences of ambient noise correlations based on a seismic interferometry method. Still, these phase differences are influenced not only by seismic wave velocity changes but also by errors in clock timing in seismometers. The clock drift occurs due to out-of-synchronization with the GPS clock and the drift of the internal clock. Therefore, to accurately monitor temporal changes in crustal structure by measuring the phase differences of noise correlations, it is crucial to evaluate the contribution of errors in clock timing to the phase differences. Recently, a method using an extended Kalman filter based on a state-space model was developed for reliable detection of temporal changes in the waveforms of ambient noise correlations, with the state-space model offering the advantage of flexible modelling of time-series data. In this study, we incorporated the time-shifts caused by clock time errors of the seismometer into the state-space model of the temporal changes in ambient noise correlations. We estimated seismic velocity changes, amplitude changes of noise correlations and clock time errors from 2010 April to 2021 September at seismic stations around the Shinmoe-dake volcano in Japan, which experienced eruptions in 2011 and 2018, respectively. Several stations exhibited clear clock time offsets, and the occurrence of clock time-shifts coincided with the dates when the data logger was turned off for seismic station maintenance or replacement of the seismometer. The proposed method provides stable estimations with respect to the signal-to-noise ratio of the waveform, and this stable estimation facilitates accurate timing of seismic recordings, enabling precise analysis of seismic phase arrival times.

Analisis Deformasi di Lampung dan Selat Sunda berdasarkan Data GNSS tahun 2018 hingga 2021

The Lampung Province and Sunda Strait have a seismic gap zone with the potential for major earthquakes in the future. This study analyzes the deformation occurring in this region using continuous Global Navigation Satellite System (GNSS) station data from Indonesia Continuously Operating Reference Station (InaCORS) and Sumatran GPS Array (SuGAr) from 2018 to 2021.5. The GNSS data was processed using the Bernese 5.2 scientific software, applying least squares for velocity changes and statistical tests to analyze significance. The data processing was carried out in two schemes: the first scheme covering 2018-2020, and the second covering 2019-2021. The results of the deformation analysis from 2018 to 2021, using two continuous GNSS data processing schemes, showed velocity changes relative to the Sundaland Plate ranging from ~2 mm/year to ~20 mm/year. In the eastern region of the Sumatra fault, the velocity changes were smaller, around ~5 mm/year, due to the minimal influence of tectonic activity. However, in the Sunda Strait region, the deformation was influenced by volcanic activity. The deformation occurring in Lampung Province and the Sunda Strait, based on GNSS velocity changes, significantly contributes to tectonic and volcanic activities.

Evaluation of Earthquake Performance of Reinforced Concrete Buildings with Fuzzy Logic Method

Reinforced concrete building is a type of building whose structural system consists of reinforced concrete columns, beams, wall wall, slabs and foundations. Due to the strong earthquakes in our country, it has been observed that reinforced concrete buildings have been severely damaged or collapsed from past to present. The evaluation of buildings after earthquakes and their performances are quite important for the safety of life and property. In the literature, different methods have been developed for post-earthquake evaluation of buildings, either low-cost and fast or slow and high-cost and tool-demanding. In this study, in order to overcome the gap in the literature, the evaluation of earthquake performance of buildings with fuzzy logic method is discussed. In this context, the performance of the building was evaluated by taking into account the concrete compressive strength, number of floors, ground floor area, area of column and shear walls, ground floor area and architectural parameters. The data of 18, 28 and 146 buildings affected by the earthquakes in Afyon, Bingöl and Van provinces in 2002, 2003 and 2011, respectively, were used in the study. Out of the total 192 building data, 94 buildings were processed as data and fuzzy logic rules were applied. The remaining 98 buildings were tested with this method. The buildings considered are light, moderate and severely damaged or collapsed. Matlab program was used in the study. The method and result used in the study are in the class of fast and second stage methods in the literature, which are reliable but do not reflect the final result. With the method revealed in the study, the damage status of reinforced concrete buildings can be determined quickly and reliably. In this context, the buildings were evaluated as no damage/slight, moderate and severely damaged or collapsed. According to the results, the correct estimation rate of the damage status of 98 buildings was found to be 88%. It is seen that the fuzzy logic method can predict building performances at a quite high rate. Compared to the existing second stage assessment methods in the literature, this result is more conservative.

Imaging microearthquake rupture processes using a dense array in Oklahoma

Both large and small earthquakes rupture in complex ways. However, microearthquakes are often simplified as point sources and their rupture properties are challenging to resolve. We leverage seismic wavefields recorded by a dense array in Oklahoma to image microearthquake rupture processes. We construct machine-learning enabled catalogs and identify four spatially disconnected seismic clusters. These clusters likely delineate near-vertical strike-slip faults. We develop a new approach to use the maximum absolute SH-wave amplitude distributions (S-wave wavefields) to compare microearthquake rupture processes. We focus on one cluster with earthquakes located beneath the dense array and have a local magnitude range of -1.3 to 2.3. The S-wave wavefields of single earthquakes are generally coherent but differ slightly between the low-frequency (&lt;12 Hz) and high-frequency (&gt;12 Hz) bands. The S-wave wavefields are coherent between different earthquakes at low frequencies with average correlation coefficients greater than 0.95. However, the wavefield coherence decreases with increasing frequency for different earthquakes. This reduced coherence is likely due to the rupture differences among individual earthquakes. Our results suggest that earthquake slip of the microearthquakes dominates the radiated S-wave wavefields at higher frequencies. Our method suggests a new direction in resolving small earthquake source attributes using dense seismic arrays without assuming a rupture model.

Coseismic deformation analysis of the 2017 Milin Ms 6.9 earthquake in the Namche Barwa Syntaxis: Implications for regional tectonics

The 2017 Milin earthquake occurred near Namche Barwa on the southeastern Tibetan Plateau. This research employed Sentinel-1 imagery to capture the corresponding coseismic deformation field and adopted a joint inversion approach using geodetic measurements and teleseismic waveform data to assess earthquake mechanics. The coseismic deformation indicated distinct differential movements, indicating a thrust movement of this fault with uplift in northeast side. The calculated seismic moment magnitude was Mw 6.5, with the fault striking of 305° and a dip angle of 72° The earthquake hypocenter was located at a depth of approximately 10 km, and the maximum fault slip was approximately 1.08 m. The event was attributed to the northwestern segment of the Xixingla fault. Numerical simulations indicated a peak ground velocity of approximately 0.3 m/s. The northwestern segment of the Xixingla fault was identified as the seismogenic fault of the event. Furthermore, the Coulomb stress analysis indicates that the earthquake induced stress loading on the Dongju–Milin and Jiali faults. The surface strain analysis reveals a region of high strain adjacent to the fault. The 2017 Milin earthquake occurred because the Xixingla fault underwent predominantly thrust-type sliding due to the continuous tectonic stress induced within the Tibetan Plateau. Notably, the middle segment of the Xixingla fault currently experiences relatively low seismic activity and is in a state of stress and strain accumulation. Thus, strong earthquakes may occur within the middle segment of the Xixingla fault.

The rupture process of the Hualien M7.3 sequence on April 3, 2024

The Hualien M7.3 earthquake on April 3, 2024, was a significant and strong earthquake in Taiwan, China in the past two decades. The rupture process of the main shock and strong aftershocks is of great significance to the subsequent seismic activity and seismogenic tectonic research. Based on local strong-motion data, we used the IDS (Iterative Deconvolution and Stacking) method to obtain the rupture process of the mainshock and two strong aftershocks on the 23rd. The rupture of the mainshock was mainly unilateral, lasting 31 s, with a maximum slip of 2 m, and the depth of the large slip zone is about 41-49 km. There is a clear difference between the rupture depth of the main shock and the two strong aftershocks. The depths of the large slip zones of the latter two are 3-9 km and 8-10 km, respectively. There is also a significant difference in the seismogenic fault between the mainshock and the aftershocks, and we believe that there are two seismogenic fault zones in the study area, the deep and the shallow fault zone. The slip of the deep faults activates the shallow faults.

Resolution analysis of the gravity survey network in the middle and south sections of Tan-Lu fault and recent changes in the gravity field

In this paper, we computed the fractal dimension of three survey areas within the central and southern sections of the Tan-Lu fault zone using fractal analysis. Subsequently, simulations were conducted to analyze the gravity response under a forward model of equivalent density changes. Additionally, we thoroughly investigated the seismic monitoring capabilities of the gravity network in the central and southern regions of the Tan-Lu fault. Expanding on these analyses. Recent gravity field variations were examined in the mid-southern segment of the Tan-Lu fault zone and its surrounding areas from 2013 to 2023. The results indicate that the observation capabilities of the northern network in the study area outperform those of the southern gravity network, with the northern network demonstrating a more evenly distributed coverage. The optimal gravity anomaly recovery effect for the entire study area is achieved at a resolution of 0.5° × 0.5°. With an equivalent observable signal in the range of 30 × 10−8 m/s2 to 40 × 10−8 m/s2, the spatial resolution of the gravity network's field source is estimated to be approximately 55 km. From 2013 to 2023, a significant positive change has been observed in the gravity field within the study area. The Tan-Lu fault zone plays a crucial role in governing the crustal movement in this region, with the dextral strike-slip movement trend of the fault persisting. Small earthquakes occur more frequently in the southern section of the fault zone, while strong earthquakes are less common. The alignment of gravity field changes with the fault strike indicates ongoing activity in the fault zone without any signs of locking. In the central segment of the Tan-Lu fault zone in the Shandong region, there appears to be a weaker correlation between gravity field changes and fault trends. This discrepancy may suggest that the area is locked, resulting in the accumulation of stress and strain. It is imperative to monitor the continuous evolution of the gravity field in this region to gain insights into potential seismic risks.