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.

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.

Spatial correlations in ground motion intensity measuring from the 2023 Türkiye earthquake

When evaluating an area's seismic risk or resilience,it is necessary to use the spatial correlation to analyze the ground motion parameters of multiple sites together in an earthquake.These two large earthquakes in Türkiye provided the possibility for spatial correlation analysis of ground motion intensity measurements in this area.Based on the strong motion records provided by The Disaster and Emergency Management Authority of Türkiye(AFAD),this study uses the local ground motion prediction equation in Türkiye to give spatial correlation analysis of Intensity Measurements.This study gives an exponential model based on a semivariogram and compares it with the correlation model obtained from previous studies.

Preliminary report of coseismic surface rupture(part)of Türkiye's M_(W)7.8 earthquake by remote sensing interpretation

Both M_(W) 7.8 and M_(W) 7.5 earthquakes occurred in southeastern Türkiye on February 6,2023,resulting in numerous buildings collapsing and serious casualties.Understanding the distribution of coseismic surface ruptures and secondary disasters surrounding the epicentral area is important for post-earthquake emergency and disaster assessments.High-resolution Maxar and GF-2 satellite data were used after the events to extract the location of the rupture surrounding the first epicentral area.The results show that the length of the interpreted surface rupture zone(part of)is approximately 75 km,with a coseismic sinistral dislocation of 2-3 m near the epicenter;however,this reduced to zero at the tip of the southwest section of the East Anatolia Fault Zone.Moreover,dense soil liquefaction pits were triggered along the rupture trace.These events are in the western region of the Eurasian Seismic Belt and result from the subduction and collision of the Arabian and African Plates toward the Eurasian Plate.The western region of the Chinese mainland and its adjacent areas are in the eastern section of the Eurasian Seismic Belt,where seismic activity is controlled by the collision of the Indian and Eurasian Plates.Both China and Türkiye have independent tectonic histories.

Rapid source inversions of the 2023 SE Türkiye earthquakes with teleseismic and strong-motion data

We conducted rapid inversions of rupture process for the 2023 earthquake doublet occurred in SE Türkiye,the first with a magnitude of M_(W)7.8 and the second with a magnitude of M_(W)7.6,using teleseismic and strong-motion data.The teleseismic rupture models of the both events were obtained approximately 88 and 55 minutes after their occurrences,respectively.The rupture models indicated that the first event was an asymmetric bilateral event with ruptures mainly propagating to the northeast,while the second one was a unilateral event with ruptures propagating to the west.This information could be useful in locating the meizoseismal areas.Compared with teleseismic models,the strong-motion models showed relatively higher resolution.A noticeable difference was found for the M_(W)7.6 earthquake,for which the strong-motion models shows a bilateral event,rather than a unilateral event,but the dominant rupture direction is still westward.Nevertheless,all strong-motion models are consistent with the teleseismic models in terms of magnitudes,durations,and dominant rupture directions.This suggests that both teleseismic and strong-motion data can be used for fast determination of major source characteristics.In contrast,the strong-motion data would be preferable in future emergency responses since they are recorded earlier and have a better resolution ability on the source ruptures.

Geophysical Reactions to Remote 2022 Tonga Eruption and to Türkiye Earthquakes in Georgia (Caucasus): Hydrogeology, Geomagnetics and Seismicity

The paper is devoted to analysis of hydrogeological, geomagnetic and seismic response to the two great remote geophysical events, 2022 Tonga volcano eruption and 2020-2023 Türkiye earthquakes in Georgia (Caucasus). The geophysical observation system in Georgia, namely, water level stations in the network of deep wells, atmospheric pressure and the geomagnetic sensors of the Dusheti Geophysical Observatory (DGO) as well as seismic data in Garni Observatory (Armenia) respond to the Tonga event by anomalies in the time series. These data show that there are two types of respond: infrasound disturbances in atmospheric pressure and seismic waves in the Earth generated by the eruption. After Tonga eruption January 15 at 04:21 UTC three groups of N-shaped waveforms were registered in the water level corresponding to the global propagation characteristics of the N-shaped waveform of infrasound signals on the barograms generated by eruption at the distance ~15,700 km: they were identified as the Lamb wave, a surface wave package running in the atmosphere with a velocity around ~314 m/s. The paper also presents the WL reactions to three strong EQs that occur in Türkiye 2020-2023, namely Elazığ, Van and Türkiye-Syria EQs. WL in Georgian well network reacts to these events by anomalies of different intensity, which points to the high sensitivity of hydrosphere to remote (several hundred km) strong EQs. The intensity and character of WL reactions depend strongly on the local hydrogeological properties of rocks, surrounding the well.

Tectonic geomorphology of Türkiye and its insights into the neotectonic deformation of the Anatolian Plate

Quantitative geomorphic analyses are usually powerful in identifying active tectonics across global orogenic belts.Our present study will focus on the Anatolian Plate which hosts a lot of recent catastrophic earthquakes in Türkiye.Six geomorphic indices for 100 sub-basins around Türkiye have been computed including local relief,slope,normalized steepness index(k_(Sn)),hypsometric curve and integral(HI),transverse topographic symmetry factor(Tf),and the basin asymmetry factor(Af).The averaged kSnand Af values have shown four high-value anomalous zones,suggesting relatively high uplift rates featured by high river incision and regional tilting.The values of 0.35≤HI<0.6 for basins with S-shaped curves imply intensive tectonic activities along the eastern part of the North Anatolian Fault Zone(NAFZ),the Northeast Anatolian Fault Zone(NEAFZ),the East Anatolian Fault Zone(EAFZ),and the Central Anatolian Fault Zone(CAFZ).All results of the geomorphic indices analysis suggest a relatively high degree of tectonic activity in the following four areas,the Isparta Angle,the Eastern Black Sea Mountains,the South-eastern Anatolia Region,and the Central Anatolian fault zone.We further suggest that the eastern part of the NAFZ,NEAFZ,EAFZ,and CAFZ will be more active in tectonic activities,with a greater potential for strong earthquake occurrence.

High-resolution seismicity imaging and early aftershock migration of the 2023 Kahramanmara?(SE Türkiye)M_W7.9&7.8 earthquake doublet

We build a high-resolution early aftershock catalog for the 2023 SE Türkiye seismic sequence with PALM,a seamless workflow that sequentially performs phase picking,association,location,and matched filter for continuous data.The catalog contains 29,519 well-located events in the two mainshocks rupture region during 2023-02-01–2023-02-28,which significantly improves the detection completeness and relocation precision compared to the public routine catalog.Employing the new PALM catalog,we analyze the structure of the seismogenic fault system.We find that the Eastern Anatolian Fault(EAF)that generated the first M_(W)7.9 mainshock is overall near-vertical,whereas complexities are revealed in a small-scale,such as subparallel subfaults,unmapped branches,and stepovers.The seismicity on EAF is shallow(<15 km)and concentrated in depth distribution,indicating a clear lock-creep transition.In contrast,the SürgüFault(SF)that is responsible for the second M_(W)7.8 mainshock is shovel-shaped for the nucleation segment and has overall low dip angles(~40°–80°).Aftershocks on the SF distribute in a broad range of depth,extending down to~35 km.We also analyze the temporal behavior of seismicity,discovering no immediate foreshocks within~5 days preceding the first mainshock,and no seismic activity on the SF before the second mainshock.

Tor topography in Giresun mountains,Türkiye

Tors are mound-like rock landforms,resistant to erosion,that may sometimes reach the size of several storeys of house by volume.They may develop in different climate regions and different rock types,led by granitic rocks.In this study,tors developed in granitic rocks around Yaylal?village linked to Kürtün county in Gümü?hane were investigated.The study area is located in the eastern part of the Giresun mountain range in northeastern Turkey,2000-2300 m above sea level.The aim of the study is to determine the natural environmental characteristics in tor formation and to reveal the effect of bedrock on tor formation.During field studies,18 tors were selected by noting location,size and structural features,and dimensions were measured.Thin sections were first made from the rock samples taken from the tors and then the mineralogical composition was investigated by modal analysis method.As a result of this analysis,it was determined that the bedrock forming the tors is granite.During the field observations,it was determined that the intersecting vertical joints supported the formation of castle-like tors,while the horizontal joints supported the formation of pita pileshaped tors.There are also pseudokarstic shapes resembling kamenitza and karst pavement.In conclusion,joint sets causing local resistance differences directed the formation of tors,while excess slope affected surface erosion and suitability of climate conditions-controlled weathering rate.

从“Turkey”到“Türkiye”——土耳其英语国名的改变考析

国名是具有多重意义的政治符号,考察土耳其国名英语名称的改变,既可以从一个侧面反映出土耳其人对自身定位的深化历程,同时由于英语国家人士参与了土耳其国名的修改,因此也能从中了解其对土耳其政治文化的认识程度的变化。当一个在联合国已经有百多年历史的原国家名字修改时,该如何处理才合适?该文试图从三个角度加以梳理与分析,在土耳其与英语语言国家的交流中,国家文化背景与传统准确的引导作用发挥到了极致,显示出国名作为政治符号,在对外传播中具有多重意义和重要性。

Changes in snow cover extent in the Central Taurus Mountains from 1981 to 2021 in relation to temperature, precipitation, and atmospheric teleconnections

The snow cover over the Taurus Mountains affects water supply, agriculture, and hydropower generation in the region. In this study, we analyzed the monthly Snow Cover Extent(SCE) from November to April in the Central Taurus Mountains(Bolkar, Aladaglar, Tahtali and Binboga Mountains) from 1981 to 2021. Linear trends of snow cover season(November to April) over the last 41 years showed decreases in SCE primarily at lower elevations. The downward trend in SCE was found to be more pronounced and statistically significant for only November and March. SCE in the Central Taurus Mountains has declined about-6.3% per decade for 2500-3000 m in November and about-6.0% per decade for 1000-1500 m and 3000+ m in March over the last 41 years. The loss of SCE has become evident since the 2000s, and the lowest negative anomalies in SCE have been observed in 2014, 2001, and 2007 in the last 41 years, which are consistent with an increase in air temperature and decreased precipitation. SCE was correlated with both mean temperature and precipitation, with temperature having a greater relative importance at all elevated gradients. Results showed that there is a strong linear relationship between SCE and the mean air temperature(r =-0.80) and precipitation(r = 0.44) for all elevated gradients during the snow season. The Arctic Oscillation(AO), the North Atlantic Oscillation(NAO), and the Mediterranean Oscillation(MO) winter indices were used to explain the year-to-year variability in SCE over the Central Taurus Mountains. The results showed that the inter-annual variability observed in the winter SCE on the Central Taurus Mountains was positively correlated with the phases of the winter AO, NAO and MO, especially below 2000 m elevation.

Geoelectrical survey over perched aquifers in the northern part of Upper Sakarya River Basin,Türkiye

In this study,a groundwater exploration survey was conducted using the DC Resistivity(DCR)method in a hydrogeological setting containing a perched aquifer.DCR data were gathered and an electrical tomography section was recovered using conventional four-electrode instruments with a Schlumberger array and a two-dimensional(2D)inversion scheme.The proposed scheme was tested over a synthetic threedimensional(3D)subsurface model before deploying it in a field situation.The proposed method indicated that gathering data with simple four-electrode instruments at stations along a line and 2D inversion of datasets at multiple stations can recover depth intervals of the studied aquifer in the hydrogeological setting even if it has a 3D structure.In this study,2D inversion of parallel profiles formed a pseudo-3D volume of the subsurface resistivity structures and mapped out multiple resistive(>25 ohm·m)bodies at shallow(between 50-100 m)and deep sections(>150 m).In general,the proposed method is convenient to encounter geological units that have limited vertical and spatial extensions in any direction and presents resistivity contrast from groundwater-bearing geologic materials.

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.

Representative value of cross-fault in the northeastern margin of the Qinghai-Tibet block and case analysis of the 2016 Menyuan Ms6.4 earthquake

The equation for determining cross-fault representative value is calculated based on hanging wall and foot wall reference level surfaces. The cross-fault data reliability are analyzed base on the stability of reference datum and observation points, thereby facilitating plotting of the representative value curves after removing interference. The spatial and temporal characteristics of fault deformation abnormalities before the 2016 Menyuan Ms6.4 earthquake, as well as the fault-movement characteristics reflected by representative value, are summarized. The results show that many site trends had changed 1-3 years before the Menyuan Ms6.4 earthquake in the Qilian Fault, reflecting certain background abnormalities. The short-term abnormalities centrally had appeared in the 6 months to 1 year period before the earthquake near and in the neighborhood of the source region,demonstrating a significantly increased number of short-term abnormalities. Many sites near and in the neighborhood of the source region had strengthened inverse activities or had changed from positive to inverse activities in the most recent 2-3 years, which reflect stress-field enhancements or adjustment features.

The 2023 Earthquake in Türkiye and Implications for China's Response to Catastrophe

On 6 February 2023,two 7.8 magnitude earthquakes consecutively hit south-central Türkiye,causing great concern from all governments,the United Nations,academia,and all sectors of society.Analyses indicate that there is also a high possibility of strong earthquakes with a magnitude of 7.8 or above occurring in the western region of China in the coming years.China is a country that is highly susceptible to catastrophic disasters such as earthquakes,floods,and other natural calamities,which can cause significant damages to both human life and property,as well as widespread impacts on the society.Currently,China's capacity for disaster prevention and control is still limited.In order to effectively reduce the impact of catastrophic disasters,ensure the safety of people's lives and property to the greatest extent possible,maintain social stability in high-risk areas,and ensure high-quality and sustainable regional development,it is urgent to improve the seismic resistance level of houses and critical infrastructure in high earthquake risk zones and increase the earthquake-resistant design level of houses in high-risk fault areas with frequent seismic activities;significantly enhance the ability to defend against extreme weather and ocean disasters in economically developed areas along the southeastern coast,as well as the level of fortification in response to extreme meteorological and hydrological disasters of coastal towns/cities and key infrastructure;vigorously enhance the emergency response capacity and disaster risk prevention level in western and ethnic minority regions;comprehensively improve the defense level of residential areas and major infrastructure in high geological hazard risk zones with flash floods,landslides,and mudslides;systematically promote national disaster prevention and mitigation education;and greatly enhance the societal disaster risk reduction ability,including catastrophic insurance.

苍山5.2级地震前的R(t)图异常特征及其时间变化

研究了１９９５ 年苍山５２ 级地震前地震时空分布参数Ｒ （ｔ） 图的异常变化， 由Ｒ （ｔ）图清楚的反映出苍山５２ 级地震前地震活动的三个阶段性发展特征。即早期地震活动增强，形成大的孕震空区； 中期在上述背景上地震活动再次向外围扩散并收缩， 地震活动逐渐减弱， 孕震空区缩小； 晚期地震活动出现短平静（１８ 个月无ＭＬ≥２５ 级地震） ； 临震前１８ 天孕震空区内出现３ 个特征前震，

印度洋T波与地震震源参数的相关性

从2004-12-26苏门答腊(Mw=9.3)大地震发生后,在该区域诱发的137个地震的波形记录中提取出了清晰的T波,并分析研究了T波波形特征与地震震源参数的相关性.通过分析计算得到如下结论:(Ⅰ)对于同一震级的地震,走滑型地震在海水声纳层中激发的T波能量大于正逆断层型地震;(Ⅱ)随着深度的增加,激发的T波与P波能量的比值(T/P)在逐渐减小;(Ⅲ)T与P,S波峰值到时差与地震震中距呈现较好的线性关系;(Ⅳ)在纬度1.0至2.5的区域,存在一个地震波向声波转换效率高的区域.

A Hybrid Multi-Hazard Susceptibility Assessment Model for a Basin in Elazig Province,Türkiye

Preparation of accurate and up-to-date susceptibility maps at the regional scale is mandatory for disaster mitigation,site selection,and planning in areas prone to multiple natural hazards.In this study,we proposed a novel multi-hazard susceptibility assessment approach that combines expert-based and supervised machine learning methods for landslide,flood,and earthquake hazard assessments for a basin in Elazig Province,Türkiye.To produce the landslide susceptibility map,an ensemble machine learning algorithm,random forest,was chosen because of its known performance in similar studies.The modified analytical hierarchical process method was used to produce the flood susceptibility map by using factor scores that were defined specifically for the area in the study.The seismic hazard was assessed using ground motion parameters based on Arias intensity values.The univariate maps were synthesized with a Mamdani fuzzy inference system using membership functions designated by expert.The results show that the random forest provided an overall accuracy of 92.3%for landslide susceptibility mapping.Of the study area,41.24%were found prone to multi-hazards(probability value>50%),but the southern parts of the study area are more susceptible.The proposed model is suitable for multi-hazard susceptibility assessment at a regional scale although expert intervention may be required for optimizing the algorithms.