Analysis of Small Earthquake Activities before Medium-strong Earthquakes in Northeast China

13 earthquakes with M≥5.0 have been recorded in the northeast region of China since1970,among which eight medium-strong earthquakes are independent in space and time.Studies of seismicity before these eight medium-strong earthquakes show that small earthquake activity was enhanced before the occurrences.Though seismicity increase is a common phenomenon in the northeast China region,we have difficulty in predicting the medium-strong earthquakes by this phenomenon alone.In order to predict medium-strong earthquakes through seismicity increase,this paper tries to propose a new method that calculates small earthquake frequency through the changing pattern of small earthquake activities based on the characteristics of small earthquake activity in the northeast China region.The results show that we can get the obvious anomaly frequency of small earthquakes before medium-strong earthquakes through the new method,and can obtain a medium to short term anomaly index for the northeast China region.

The MW5.5 earthquake on August 6,2023,in Pingyuan,Shandong,China:A rupture on a buried fault

On August 6,2023,a magnitude MW5.5 earthquake struck Pingyuan County,Dezhou City,Shandong Province,China.This event was significant as no large earthquakes had been recorded in the region for over a century,and no active fault had been previously identified.This study collects 1309 P-wave arrival times and 866 S-wave arrival times from 74 seismic stations less than 200 km to the epicenter to constrain the spatial distribution of the mainshock and its 125 early aftershocks by the double difference earthquake relocation method,and selects 864 P-waveforms from 288 stations located within 800 km of the epicenter to constrain the focal mechanism solution of the mainshock through centroid moment tensor inversion.The relocation and the inversion indicate,the Pingyuan MW5.5 earthquake was caused by a rupture on a buried fault,likely an extensive segment of the Gaotang fault.This buried fault exhibited a dip of approximately 75°to the northwest,with a strike of 222°,similar to the Gaotang fault.The rupture initiated at the depth of 18.6 km and propagated upward and northeastward.However,the ground surface was not broken.The total duration of the rupture was~6.0 s,releasing the scalar moment of 2.5895×1017 N·m,equivalent to MW5.54.The moment rate reached the maximum only 1.4 seconds after the rupture initiation,and the 90%scalar moment was released in the first 4.6 s.In the first 1.4 seconds of the rupture process,the rupture velocity was estimated to be 2.6 km/s,slower than the local S-wave velocity.As the rupture neared its end,the rupture velocity decreased significantly.This study provides valuable insights into the seismic characteristics of the Pingyuan MW5.5 earthquake,shedding light on the previously unidentified buried fault responsible for the seismic activity in the region.Understanding the behavior of such faults is crucial for assessing seismic hazards and enhancing earthquake preparedness in the future.

Characterization and spatial analysis of coseismic landslides triggered by the Luding Ms 6.8 earthquake in the Xianshuihe fault zone, Southwest China

On September 5, 2022, a magnitude Ms 6.8 earthquake occurred along the Moxi fault in the southern part of the Xianshuihe fault zone located in the southeastern margin of the Tibetan Plateau,resulting in severe damage and substantial economic loss. In this study, we established a coseismic landslide database triggered by Luding Ms 6.8 earthquake, which includes 4794 landslides with a total area of 46.79 km^(2). The coseismic landslides primarily consisted of medium and small-sized landslides, characterized by shallow surface sliding. Some exhibited characteristics of high-position initiation resulted in the obstruction or partial obstruction of rivers, leading to the formation of dammed lakes. Our research found that the coseismic landslides were predominantly observed on slopes ranging from 30° to 50°, occurring at between 1000 m and 2500 m, with slope aspects varying from 90° to 180°. Landslides were also highly developed in granitic bodies that had experienced structural fracturing and strong-tomoderate weathering. Coseismic landslides concentrated within a 6 km range on both sides of the Xianshuihe and Daduhe fault zones. The area and number of coseismic landslides exhibited a negative correlation with the distance to fault lines, road networks, and river systems, as they were influenced by fault activity, road excavation, and river erosion. The coseismic landslides were mainly distributed in the southeastern region of the epicenter, exhibiting relatively concentrated patterns within the IX-degree zones such as Moxi Town, Wandong River basin, Detuo Town to Wanggangping Township. Our research findings provide important data on the coseismic landslides triggered by the Luding Ms 6.8 earthquake and reveal the spatial distribution patterns of these landslides. These findings can serve as important references for risk mitigation, reconstruction planning, and regional earthquake disaster research in the earthquake-affected area.

The 2023 Turkey earthquake doublet: Earthquake relocation, seismic tomography, and stress field inversion

On February 6,2023,two earthquakes with magnitudes of M_(W) 7.8 and M_(W) 7.5 struck southeastern Turkey,causing significant casualties and economic losses.These seismic events occurred along the East Anatolian Fault Zone,a convergent boundary between the Arabian Plate and the Anatolian Subplate.In this study,we analyze the M_(W) 7.8 and M_(W) 7.5 earthquakes by comparing their aftershock relocations,tomographic images,and stress field inversions.The earthquakes were localized in the upper crust and exhibited steep dip angles.Furthermore,the aftershocks occurred either close to the boundaries of low and high P-wave velocity anomaly zones or within the low P-wave velocity anomaly zones.The East Anatolia Fault,associated with the M_(W) 7.8 earthquake,and the SürgüFault,related to the M_(W) 7.5 earthquake,predominantly experienced shear stress.However,their western sections experienced a combination of strike-slip and tensile stresses in addition to shear stress.The ruptures of the M_(W) 7.8 and M_(W) 7.5 earthquakes appear to have bridged a seismic gap that had seen sparse seismicity over the past 200 years prior to the 2023 Turkey earthquake sequence.

Rapid report of source parameters of 2023 M6.2 Jishishan,Gansu earthquake sequence

The M6.2 earthquake in Jishishan,Gansu Province,on December 18,2023,caused extraordinary earthquake disasters.It was located in the northern part of the north−south seismic zone,which is a key area for earthquake monitoring in China.The newly built dense strong motion stations in this area provide unprecedented conditions for high-precision earthquake relocation,especially the earthquake focal depth.This paper uses the newly built strong motion and traditional broadband seismic networks to relocate the source locations of the M3.0 and above aftershocks and to invert their focal mechanisms.The horizontal error of earthquake location is estimated to be 0.5−1 km,and the vertical error is 1−2 km.The focal depth range of aftershocks is 9.6−14.6 km,distributed in a 12-km-long strip with SSE direction.Aftershocks in the south are more concentrated horizontally and vertically,while aftershocks in the north are more scattered.The focal mechanisms of the main shock and aftershocks are relatively consistent,and the P-axis orientation is consistent with the regional strain direction.There is a seismic blank area of M3.0 and above,about 3−5 km between the main shock and aftershocks.It is suggested that the energy released by the main shock rupture is concentrated in this area.Based on the earthquake location and focal mechanism of the main shock,it is inferred that the Northern Lajishan fault zone is the seismogenic structure of the main shock,and the main shock did not occur on the main fault,but on a secondary fault.The initial rupture depth and centroid depth of the main shock were 12.8 and 14.0 km,respectively.The source rupture depth may not be the main reason for the severe earthquake disaster.

Key technologies of earthquake early warning system for China’s high-speed railway

Purpose–The purpose of this study is to introduce the top-level design ideas and the overall architecture of earthquake early-warning system for high speed railways in China,which is based on P-wave earthquake early-warning and multiple ways of rapid treatment.Design/methodology/approach–The paper describes the key technologies that are involved in the development of the system,such as P-wave identification and earthquake early-warning,multi-source seismic information fusion and earthquake emergency treatment technologies.The paper also presents the test results of the system,which show that it has complete functions and its major performance indicators meet the design requirements.Findings–The study demonstrates that the high speed railways earthquake early-warning system serves as an important technical tool for high speed railways to cope with the threat of earthquake to the operation safety.The key technical indicators of the system have excellent performance:The first report time of the P-wave is less than three seconds.From the first arrival of P-wave to the beginning of train braking,the total delay of onboard emergency treatment is 3.63 seconds under 95%probability.The average total delay for power failures triggered by substations is 3.3 seconds.Originality/value–The paper provides a valuable reference for the research and development of earthquake early-warning system for high speed railways in other countries and regions.It also contributes to the earthquake prevention and disaster reduction efforts.

Identification of earthquake induced structural damage based on synchroextracting transform

Several popular time-frequency techniques,including the Wigner-Ville distribution,smoothed pseudo-Wigner-Ville distribution,wavelet transform,synchrosqueezing transform,Hilbert-Huang transform,and Gabor-Wigner transform,are investigated to determine how well they can identify damage to structures.In this work,a synchroextracting transform(SET)based on the short-time Fourier transform is proposed for estimating post-earthquake structural damage.The performance of SET for artificially generated signals and actual earthquake signals is examined with existing methods.Amongst other tested techniques,SET improves frequency resolution to a great extent by lowering the influence of smearing along the time-frequency plane.Hence,interpretation and readability with the proposed method are improved,and small changes in the time-varying frequency characteristics of the damaged buildings are easily detected through the SET method.

Identifying drivers of urban landuse changes in the Wenchuan earthquake- affected area by using night-time light data

To learn the process of urban land evolution before and after an earthquake is vital to formulate the urban reconstruction control policies and recovery measures in the earthquake-stricken areas.However,spatiotemporal evolution and its driving factors of urban land in earthquake-prone areas remains limited due to the scarcity of ground observation data.This research,leveraging night-time light remote sensing imagery and land cover data,conducted a comprehensive analysis of the long-term evolution characteristics of urban land in earthquake-prone areas.It introduced methodologies for assessing the socio-economic impact and the primary natural environmental factors driving urban land evolution in these regions.To validate the proposed methods,the 2008 Wenchuan earthquake-affected area in China was selected as a representative study area.The results indicated that the average Digital Number(DN)values in socio-economically impacted areas showed a trend of rising,falling,and then rising again after the earthquake.DN values in three types of damaged areas including Type Ⅱ,Type Ⅲ,and Type Ⅳ exceeded pre-earthquake levels.The analysis of determinative factors influencing urban land evolution revealed that slope and elevation were key elements in controlling urban land expansion before the earthquake,whereas factors such as slope,elevation,lithology,and faults had a stronger influence on urban land expansion after the earthquake.It can be seen that,in view of the differences in the natural conditions of regions for post-disaster reconstruction,the local government need to actively adjust and adapt to urban spatial planning,so as to leverage the scale effect of large-scale inputs of funds,facilities,human resources and other factors after the disaster,thus enhancing resilience and recovery efficiency in response to disaster impacts.

Kinematic deformation and intensity assessment of the 2021 Maduo M_(S)7.4 earthquake in Qinghai revealed by high-frequency GNSS

Rapid acquisition of the kinematic deformation field and seismic intensity distribution of large earthquakes is crucial for postseismic emergency rescue,disaster assessment,and future seismic risk research.The advancement of GNSS observation and data processing makes it play an important role in this field,especially the high-frequency GNSS.We used the differential positioning method to calculate the 1 HZ GNSS data from 98 sites within 1000 km of the M_(S)7.4 Maduo earthquake epicenter.The kinematic deformation field and the distribution of the seismic intensity by using the peak ground velocity derived from displacement waveforms were obtained.The results show that:1)Horizontal coseismic response deformation levels ranging from 25 mm to 301 mm can be observed within a 1000 km radius from the epicenter.Coseismic response deformation on the east and west sides shows bilateral asymmetry,which markedly differs from the symmetry presented by surface rupture.2)The seismic intensity obtained through high-frequency GNSS and field investigations exhibits good consistency of the scope and orientation in the high seismic intensity area,although the former is generally slightly smaller than the latter.3)There may exist obstacles on the eastern side of the seismogenic fault.The Maduo earthquake induced a certain tectonic stress loading effect on the western Kunlun Pass-Jiangcuo fault(KPJF)and Maqin-Maqu segment,resulting in higher seismic risk in the future.

Analytical Solution and Simplified Formula for Added Mass of Horizontal and Vertical Motions of Truncated Cylinders Under Earthquake Action

This paper investigates the hydrodynamic characteristics of floating truncated cylinders undergoing horizontal and vertical motions due to earthquake excitations in the finite water depth.The governing equation of the hydrodynamic pressure acting on the cylinder is derived based on the radiation theory with the inviscid and incompressible assumptions.The governing equation is solved by using the method of separating variables and analytical solutions are obtained by assigning reasonable boundary conditions.The analytical result is validated by a numerical model using the exact artificial boundary simulation of the infinite water.The main variation and distribution characteristics of the hydrodynamic pressure acting on the side and bottom of the cylinder are analyzed for different combinations of wide-height and immersion ratios.The added mass coefficient of the cylinder is calculated by integrating the hydrodynamic pressure and simplified formulas are proposed for engineering applications.The calculation results show that the simplified formulas are in good agreement with the analytical solutions.

Reactive Transport Process of Earthquake-induced Hydrochemical Changes in Guanding Thermal Spring,Western Sichuan,China

Earthquake-related hydrochemical changes in thermal springs have been widely observed;however,quantitative modeling of the reactive transport process is absent.In the present study,we apply reactive transport simulation to capture the hydrochemical responses in a thermal spring following the Wenchuan Ms 8.0 and Lushan Ms 7.0 earthquakes.We first constrain deep reservoir geothermal fluid compositions and temperature by multicomponent geothermometry,and then a reactive geochemical transport model is constructed to reproduce the hydrochemical evolution process.The results show that the recharge from the shallow aquifer increases gradually until it reaches a peak because of the permeability enhancement caused by the Lushan earthquake,which may be the mechanism to explain the earthquake-related hydrochemical responses.In contrast to the postseismic effect of the Wenchuan earthquake,the chemical evolution can be considered as hydrochemical anomalies related to the Lushan earthquake.This study proves that the efficient simulation of reactive transport processes is useful for investigating earthquake-related signals in hydrochemical time series.

The Formation of Oscillation Patterns Based on the Planetary Gravitational Field and Their Suitability for Earthquake Prediction

The fluctuating planetary gravitational field influences not only activities on the Sun but also on the Earth. A special correlation function describes the harmonics of these fluctuations. Groups of earthquakes form oscillation patterns that differ significantly from randomly chosen control groups. These patterns are suitable as an element of an AI for the probability of earthquakes.

Lithospheric strength of the Anatolian plateau and implications for strong earthquake activity in Turkey

On February 6, 2023, the doublet earthquake including two main shocks with magnitudes MW7.8 and MW7.5, occurred near the western side of the East Anatolian Fault at the southeast boundary of the Anatolian Plateau in Turkey. Based on the WGM2012 Bouguer gravity anomaly data and the Etopo1topography data, this study first introduced a joint inversion of admittance and coherence functions and used the Bayesian optimal parameter estimation method to obtain the effective elastic thickness Teand loading ratio F of the lithosphere for various tectonic units in the Anatolian Plateau. Secondly, we discussed the characteristics and influencing factors of the lithospheric mechanical strength and analyzed its relationship with seismic activity. The lithospheric mechanical strength of the Anatolian Plateau showed clear lateral heterogeneity and a "weak-strong-weak" spatial pattern from east to west,reflecting various tectonic processes. At last, the strong seismic activity was found where the lithospheric strength was low in the Anatolian Plate. We also incorporated GPS strain rate and other results to investigate the tectonic background and primary causes of the MW7.8 and MW7.5 doublet earthquakes in Turkey. The results have a good insight into urban safety design in the Turkish region, including postdisaster rehabilitation, earthquake hazard assessment, and loss reduction.

Stress triggering effect on the 2022 Honghe M_(S)5.0 earthquake with historical strong earthquakes

The 2022 Honghe M_(S)5.0 seismic event is intriguing due to its occurrence in the south of the Red River Fault,an area historically lacking seismic activities greater than M_(S)5.0.To elucidate the seismogenic mechanism and scrutinize stress-triggered interactions,we calculated co-seismic and post-seismic Coulomb stress alterations induced by nine historical seismic events(M≥6.0).The analysis reveals that these substantial seismic events provoked co-seismic stress augmentations of 1.409 bar and postseismic stress increments of 0.159 bar.Noteworthy seismic events,such as the 1833 Songming,1877Shiping,1913 Eshan,and 1970 Tonghai earthquakes,catalyzed the occurrence of the Honghe earthquake.Areas of heightened future seismic risk include the southern region of the Red River Fault and the eastern segments of the Shiping-Jianshui and Qujiang faults.Additionally,we assessed the correlation between the spatial distribution of aftershocks and the Coulomb stress shift triggered by the mainshock,taking into account the influence of calculation parameter settings.

Investigating the reactivation of historical landslides during the 2022 Luding M_(S)6.8 earthquake

On September 5,2022,a strong earthquake with a magnitude of MS6.8 struck Luding County in Sichuan Province,China,triggering thousands of landslides along the Dadu River in the northwest-southeast(NW-SE)direction.We investigated the reactivation characteristics of historical landslides within the epicentral area of the Luding earthquake to identify the initiation mechanism of earthquake-induced landslides.Records of the two newly triggered and historical landslides were analyzed using manual and threshold methods;the spatial distribution of landslides was assessed in relation to topographical and geological factors using remote sensing images.This study sheds light on the spatial distribution patterns of landslides,especially those that occur above historical landslide areas.Our results revealed a similarity in the spatial distribution trends between historical landslides and new ones induced by earthquakes.These landslides tend to be concentrated within a range of 0.2 km from the river and 2 km from the fault.Notably,both rivers and faults predominantly influenced the reactivation of historical landslides.Remarkably,the reactivated landslides are characterized by their small to medium size and are predominantly situated in historical landslide zones.The number of reactivated landslides surpassed that of previously documented historical landslides within the study area.We provide insights into the critical factors responsible for historical landslides during the 2022 Luding earthquake,thereby enhancing our understanding of the potential implications for future co-seismic hazard assessments and mitigation strategies.

Seismic response of a mid-story isolated structure considering SSI in mountainous areas under long-period earthquakes

At present,there is not much research on mid-story isolated structures in mountainous areas.In this study,a model of a mid-story isolated structure considering soil-structure interaction(SSI)in mountainous areas is established along with a model that does not consider SSI.Eight long-period earthquake waves and two ordinary earthquake waves are selected as inputs for the dynamic time history analysis of the structure.The results show that the seismic response of a mid-story isolated structure considering SSI in mountainous areas can be amplified when compared with a structure that does not consider SSI.The structure response under long-period earthquakes is larger than that of ordinary earthquakes.The structure response under far-field harmonic-like earthquakes is larger than that of near-fault pulse-type earthquakes.The structure response under near-fault pulse-type earthquakes is larger than that of far-field non-harmonic earthquakes.When subjected to long-period earthquakes,the displacement of the isolated bearings exceeded the limit value,which led to instability and overturning of the structure.The structure with dampers in the isolated story could adequately control the nonlinear response of the structure,effectively reduce the displacement of the isolated bearings,and provide a convenient,efficient and economic method not only for new construction but also to retrofit existing structures.

Co-seismic surface displacement of the June 21, 2022 M_W6 Khōst,Afghanistan earthquake from InSAR observations

A robust estimation of the earthquake location, seismic moment, and fault geometry is essential for objective seismic hazard assessment. Seismic events in a remote location, specifically in the absence of seismic and GNSS networks, can be investigated effectively using the In SAR-based technique. This study adopts the Differential Interferometric SAR(DIn SAR) technique to quantify the co-seismic surface displacement caused by the June 21, 2022, Khōst M_(W)6 earthquake that occurred along the western plate boundary between the Indian and Eurasian plate. The interferograms show that the maximum surface deformation occurred on the northwest and southwest of the fault line. From coherence, the Line of Sight(LOS) displacement, and the co-seismic surface displacement analysis, it has been observed that surface deformation was most pronounced in the southwest region of the fault line, and the surface has moved to the opposite direction along the fault line, which indicates a sinistral slightly oblique strike-slip movement. This In SAR-based observation appears consistent with the seismic waveforms derived from co-seismic surface displacements. Further, it has been argued that the slip deficit accumulated during the period of the last about 48 years along the frontal region of the northward extension of the Suleiman range and associated fault zone is qualitatively estimated at about 1.5 m, which is consistent with the seismic waveforms derived finite slip model.

Strong ground motion characteristics observed in the February 6,2023 M_(W)7.7 Türkiye earthquake

Türkiye is located in a seismically active region,where the Anatolian,African,and Arabian tectonic plates converge.High seismic hazards cause the region to be struck repeatedly by major earthquakes.On February 06,2023,a devastating M_(W)7.7 earthquake struck Türkiye at 01:17 am local time(01:17 UTC).In this regard,near and far-field ground motion data within the distance of 120 km are compiled and later characterized to identify the key ground motion intensity measures.Additionally,the vertical components of ground motions were examined to capture the complete three-dimensional nature of the seismic event.Moreover,the effect of Pulse-Like(PL)and Non-Pulse-Like(NPL)ground motion on a representative RC frame structure built as per the Türkiye code was investigated.The results indicate that PL behavior was observed in both horizontal and vertical components of ground motions and PL behavior were noted both near the epicenter and at higher distances from the epicenter.Moreover,the ratio of the peak vertical acceleration to peak horizontal acceleration at certain stations was found to be close to 1.Finally,the non-linear time history analysis of the representative reinforced concrete frame structure for ground motions recorded at stations located equidistant from the epicenter,indicated that PL ground motions led to more significant damage compared to NPL ground motions.

Coseismic deformation and fault slip distribution of the 2023 M_(W)7.8 and M_(W)7.6 earthquakes in Türkiye

On February 6,2023,a devastating earthquake with a moment magnitude of M_(W)7.8 struck the town of Pazarcik in south-central Türkiye,followed by another powerful earthquake with a moment magnitude of M_(W)7.6 that struck the nearby city of Elbistan 9 h later.To study the characteristics of surface deformation caused by this event and the influence of fault rupture,this study calculated the static coseismic deformation of 56 stations and dynamic displacement waveforms of 15 stations using data from the Turkish national fixed global navigation satellite system(GNSS)network.A maximum static coseismic displacement of 0.38 m for the M_(W)7.8 Kahramanmaras earthquake was observed at station ANTE,36 km from the epicenter,and a maximum dynamic coseismic displacement of 4.4 m for the M_(W)7.6 Elbistan earthquake was observed at station EKZ1,5 km from the epicenter.The rupture-slip distributions of the two earthquakes were inverted using GNSS coseismic deformation as a constraint.The results showed that the Kahramanmaras earthquake rupture segment was distinct and exposed on the ground,resulting in significant rupture slip along the Amanos and Pazarcik fault segments of the East Anatolian Fault.The maximum slip in the Pazarcik fault segment was 10.7 m,and rupture occurred at depths of 0–15 km.In the Cardak fault region,the Elbistan earthquake caused significant ruptures at depths of 0–12 km,with the largest amount of slip reaching 11.6 m.The Coulomb stress change caused by the Kahramanmaras earthquake rupture along the Cardak fault segment was approximately 2 bars,and the area of increased Coulomb stress corresponded to the subsequent rupture region of the M_(W)7.6 earthquake.Thus,it is likely that the M_(W)7.8 earthquake triggered or promoted the M_(W)7.6 earthquake.Based on the cumulative stress impact of the M_(W)7.8 and M_(W)7.6 events,the southwestern segment of the East Anatolian Fault,specifically the Amanos fault segment,experienced a Coulomb rupture stress change exceeding 2 bars,warranting further attention to assess its future seismic hazard risk.

Variations of shear-wave splitting parameters in the source region of the 2023 Türkiye doublet earthquakes

In this study,the shear-wave splitting parameters of local seismic events from the source regions of the 2023 Türkiye MW7.7 and MW7.6 doublet earthquakes(event 1 and event 2,respectively)were measured from June 1,2022,to April 25,2023,and their spatiotemporal characteristics were analyzed.The results revealed clear spatial and temporal differences.Spatially,the dominant fast-wave polarization direction at each station shows a strong correlation with the direction of the maximum horizontal principal compressive stress,as characterized by focal mechanism solutions of seismic events(MW≥3.5)near the station.The dominant fast-wave polarization direction and the regional stress field also showed a strong correlation with the intermovement of the Arabian Plate,African Plate,and Anatolian Block.Along the Nurdagi-Pazarcik fault zone,the seismic fault of event 1,stations closer to the middle of the fault where the mainshock occurred exhibited notably greater delay times than stations located towards the ends of the fault and far from the mainshock.In addition,the stations located to the east of the Nurdagi-Pazarcik fault and to the north of the Sürgüfault also exhibited large delay times.The spatial distribution of shear-wave splitting parameters obtained from each station indicates that the upper-crust anisotropy in the source area is mainly controlled by the regional stress field,which is closely related to the state of the block motion.During the seismogenic process of the MW7.7 earthquake,more stress accumulated in the middle of the Nurdagi-Pazarcik fault than at either end of the fault.Under the influence of the MW7.7 and MW7.6 events,the stress that accumulated during the seismogenic process of the earthquake doublet may have migrated towards some areas outside the aftershock intensive area after the earthquakes,and the crustal stress and its adjustment range near the outer stations increased significantly.With the exception of two stations with few effective events,all stations showed a consistent change in shear-wave splitting parameters over time.In particular,each station showed a decreasing trend in delay times after the doublet earthquakes,reflecting the obvious intensification of crustal stress adjustment in the seismogenic zone after the doublet earthquakes.With the occurrence of the earthquake doublet and a large number of aftershocks,the stress accumulated during the seismogenic process of the doublet earthquakes is gradually released,and then the adjustment range of crustal stress is also gradually reduced.