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.

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.

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.

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.

Strong Earthquake Model of North China Craton: A Case of 2023 Mw 5.5 Earthquake in Pingyuan County, Shandong Province, China

North China is one of the high-risk areas for destructive strong earthquakes in China's Mainland, with a history of numerous significant seismic events. On August 6, 2023, an Mw5.5 earthquake struck Pingyuan County, Dezhou City, in Shandong Province, China. This earthquake was the largest in the eastern North China Craton(NCC) since the Tangshan earthquake of 1976. Due to the absence of surface ruptures, the fault responsible for the Pingyuan Mw5.5 earthquake remains unclear. To reveal the subsurface geological structure near the earthquake epicenter, this study utilized highresolution two-dimensional seismic reflection profiles to interpret pre-existing faults.

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.

Research on the Extraction of OLR Anomaly Prior to Ms 7.5 Sand Point,Alaska Earthquake based on IPI Method

Based on the understanding that the seismic fault system is a nonlinear complex system,Rundle(1995)introduced the nonlinear threshold system used in meteorology to analyze the ocean-atmosphere interface and the El Ni?o Southern Oscillation into the study of seismic activity changes,and then proposed the PI method(Rundle et al.,2000a,b).Wu et al.(2011)modified the Pattern Informatics Method named MPI to extract the ionospheric anomaly by using data from DEMETER satellites which is suitable for 1–3 months short-term prediction.

Myths about Earthquakes:Quo vadis?

We are living in a world of numbers and calculations with enormous amount of pretty fast user-friendly software ready for an automatic output that may lead to a discovery or,alternatively,mislead to a deceptive conclusion,erroneous claims and predictions.As a matter of fact,nowadays,Science can disclose Natural Hazards,assess Risks,and deliver the state-of-the-art Knowledge of looming disaster in advance catastrophes along with useful Recommendations on the level of risks for decision making regarding engineering design,insurance,and emergency management.

Seismicity patterns before the 2021 Fin (Iran) doublet earthquakes using the region-time-length and time-to-failure methods

Knowledge regarding earthquake hazards and seismicity is crucial for crisis management, and the occurrence of foreshocks, seismic activity patterns, and spatiotemporal variations in seismic activity have been studied. Furthermore, the estimation of the region-time-length (RTL) parameter has been proposed to detect seismic quiescence before the occurrence of a large earthquake. In addition, the time-to-failure method has been used to estimate the time occurrence of large earthquakes. Hence, in this study, to gain deeper insight into seismic activity in the southern Zagros region, we utilized the RTL algorithm to identify the quiescence and activation phases leading to the Fin doublet earthquakes. Temporal variations in the RTL parameter showed two significant anomalies. One corresponded to the occurrence time of the first earthquake (2017-12-12);the other anomaly was associated with the occurrence time of the second event (2021-11-14). Based on a negative value of the RTL parameter observed in the vicinity of the Fin epicenters (2021), seismic quiescence (a decrease in seismicity compared to the preceding background rate) was identified. The spatial distribution of the RTL prognostic parameters confirms the appearance of seismic quiescence surrounding the epicenter of the Fin doublet earthquakes (2021). The time-to-failure method was designed using precursory events that describe the acceleration of the seismic energy release before the mainshock. Using the time-to-failure method for the earthquake catalog, it was possible to estimate both the magnitude and time of failure of the Fin doublet. Hence, the time-tofailure technique can be a useful supplementary method to the RTL algorithm for determining the characteristics of impending earthquakes.

Microearthquake reveals the lithospheric structure at midocean ridges and oceanic transform faults

Mid-ocean ridge and oceanic transforms are among the most prominent features on the seafloor surface and are crucial for understanding seafloor spreading and plate tectonic dynamics,but the deep structure of the oceanic lithosphere remains poorly understood.The large number of microearthquakes occurring along ridges and transforms provide valuable information for gaining an indepth view of the underlying detailed seismic structures,contributing to understanding geodynamic processes within the oceanic lithosphere.Previous studies have indicated that the maximum depth of microseismicity is controlled by the 600-℃isotherm.However,this perspective is being challenged due to increasing observations of deep earthquakes that far exceed this suggested isotherm along mid-ocean ridges and oceanic transform faults.Several mechanisms have been proposed to explain these deep events,and we suggest that local geodynamic processes(e.g.,magma supply,mylonite shear zone,longlived faults,hydrothermal vents,etc.)likely play a more important role than previously thought.

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.

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.

The 2024 M_(j) 7.6 Noto Peninsula, Japan earthquake caused by the fluid flow in the crust

On January 1, 2024 at 16:10:09 JST, an M_(j) 7.6 earthquake struck the Noto Peninsula in the southern part of the Sea of Japan. This location has been experiencing an earthquake swarm for more than three years. Here, we provide an overview of this earthquake, focusing on the slip distribution of the mainshock and its relationship with the preceding swarm. We also reexamined the source areas of other large earthquakes that occurred around the Sea of Japan in the past and compared them with the Matsushiro earthquake swarm in central Japan from1964 to 1968. The difference between the Matsushiro earthquake swarm and the Noto earthquake swarm is the surrounding stress field. The Matsushiro earthquake swarm was a strike-slip stress field, so the cracks in the crust were oriented vertically. This allowed fluids seeped from the depths to rise and flow out to the surface. On the other hand, the Noto area was a reverse fault stress field. Therefore, the cracks in the earth's crust were oriented horizontally. Fluids flowing underground in deep areas could not rise and spread over a wide area in the horizontal plane. This may have caused a large amount of fluid to accumulate underground, triggering a large earthquake. Although our proposed mechanism does not take into account other complex geological conditions into consideration, it may provide a simple way to explain why the Noto swarm is followed by a large earthquake while other swarms are not.

Rapid report of the December 18,2023 M_(S)6.2 Jishishan earthquake,Gansu,China

On December 18,2023,the Jishishan area in Gansu Province was jolted by a M_(S) 6.2 earthquake,which is the most powerful seismic event that occurred throughout the year in China.The earthquake occurred along the NWtrending Lajishan fault(LJSF),a large tectonic transformation zone.After this event,China Earthquake Networks Center(CENC)has timely published several reports about seismic sources for emergency responses.The earthquake early warning system issued the first alert 4.9 s after the earthquake occurrence,providing prompt notification that effectively mitigated panics,injuries,and deaths of residents.The near real-time focal mechanism solution indicates that this earthquake is associated with a thrust fault.The distribution of aftershocks,the rupture process,and the recorded amplitudes from seismic monitoring and GNSS stations,all suggest that the mainshock rupture predominately propagates to the northwest direction.The duration of the rupture process is~12 s,and the largest slip is located at approximately 6.3 km to the NNW from the epicenter,with a peak slip of 0.12 m at~8 km depth.Seismic station N0028 recorded the highest instrumental intensity,which is 9.4 on the Mercalli scale.The estimated intensity map shows a seismic intensity reaching up to IX near the rupture area,consistent with field survey results.The aftershocks(up to December 22,2023)are mostly distributed in the northwest direction within~20 km of the epicenter.This earthquake caused serious casualties and house collapses,which requires further investigations into the impact of this earthquake.

Rapid rupture characterization for the 2023 M_(S)6.2 Jishishan earthquake

On December 18, 2023, the M_(S)6.2 Jishishan earthquake occurred in the northeastern region of the QinghaiXizang Plateau, causing heavy casualties and property damage in Gansu and Qinghai Provinces. In this study,we integrate space imaging geodesy, finite fault inversion, and back-projection methods to decipher its rupture property, including fault geometry, coseismic slip distribution, rupture direction, and propagation speed. The results reveal that the seismogenic fault dips to the southwest at an angle of 29°. The major slip asperity is dominated by reverse slip and is concentrated within a depth range of 7–16 km, which explains the significant uplift near the epicenter observed by both the Sentinel-1 ascending and descending In SAR data. Moreover, the teleseismic array waveforms indicate a northwest propagating rupture with an overall slow rupture velocity of~1.91 km/s(AK array) or 1.01 km/s(AU array).

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.

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.

Analysis of faulting destruction and water supply pipeline damage from the first mainshock of the February 6,2023 Türkiye earthquake doublet

In 2023,two consecutive earthquakes exceeding a magnitude of 7 occurred in Türkiye,causing severe casualties and economic losses.The damage to critical urban infrastructure and building structures,including highways,railroads,and water supply pipelines,was particularly severe in areas where these structures intersected the seismogenic fault.Critical infrastructure projects that traverse active faults are susceptible to the influence of fault movement,pulse velocity,and ground motions.In this study,we used a unique approach to analyze the acceleration records obtained from the seismic station array(9 strong ground motion stations)located along the East Anatolian Fault(the seismogenic fault of the MW7.8 mainshock of the 2023 Türkiye earthquake doublet).The acceleration records were filtered and integrated to obtain the velocity and displacement time histories.We used the results of an on-site investigation,jointly conducted by China Earthquake Administration and Türkiye’s AFAD,to analyze the distribution of PGA,PGV,and PGD recorded by the strong motion array of the East Anatolian Fault.We found that the maximum horizontal PGA in this earthquake was 3.0 g,and the maximum co-seismic surface displacement caused by the East Anatolian Fault rupture was 6.50 m.As the fault rupture propagated southwest,the velocity pulse caused by the directional effect of the rupture increased gradually,with the maximum PGA reaching 162.3 cm/s.We also discussed the seismic safety of critical infrastructure projects traversing active faults,using two case studies of water supply pipelines in Türkiye that were damaged by earthquakes.We used a three-dimensional finite element model of the PE(polyethylene)water pipeline at the Islahiye State Hospital and fault displacement observations obtained through on-site investigation to analyze pipeline failure mechanisms.We further investigated the effect of the fault-crossing angle on seismic safety of a pipeline,based on our analysis and the failure performance of the large-diameter Thames Water pipeline during the 1999 Kocaeli earthquake.The seismic method of buried pipelines crossing the fault was summarized.

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.