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.

Algebraic processing technique for extracting frequency-dependent shear-wave splitting parameters in an anisotropic medium

Based on the dual source cumulative rotation technique in the time-domain proposed by Zeng and MacBeth（1993）,a new algebraic processing technique for extracting shear-wave splitting parameters from multi-component VSP data in frequency-dependent medium has been developed.By using this dual source cumulative rotation technique in the frequency-domain（DCTF）,anisotropic parameters,including polarization direction of the shear-waves and timedelay between the fast and slow shear-waves,can be estimated for each frequency component in the frequency domain.It avoids the possible error which comes from using a narrow-band filter in the current commonly used method.By using synthetic seismograms,the feasibility and validity of the technique was tested and a comparison with the currently used method was also given.The results demonstrate that the shear-wave splitting parameters frequency dependence can be extracted directly from four-component seismic data using the DCTF.In the presence of larger scale fractures,substantial frequency dependence would be found in the seismic frequency range,which implies that dispersion would occur at seismic frequencies.Our study shows that shear-wave anisotropy decreases as frequency increases.

Shear-wave splitting in the crust:Regional compressive stress from polarizations of fast shear-waves

When propagating through anisotropic rocks in the crust, shear-waves split into faster and slower components with almost orthogonal polarizations. For nearly vertical propagation the polarization of fast shear- wave （PFS） is parallel to both the strike of the cracks and the direction of maximum horizontal stress, therefore it is possible to use PFS to study stress in the crust. This study discusses several examples in which PFS is applied to deduce the compressive stress in North China, Longmenshan fault zone of east edge of Tibetan plateau and Yunnan zone of southeast edge of Tibetan plateau, also discusses temporal variations of PFS orientations of 1999 Xiuyan earthquake sequences of northeastern China. The results are consistent to those of other independent traditional stress measurements. There is a bridge between crustal PFS and the crustal principal compressive stress although there are many unclear disturbance sources. This study suggests the PFS results could be used to deduce regional and in situ principal compressive stress in the crust only if there are enough seismic stations and enough data. At least, PFS is a useful choice in the zone where there are a large number of dense seismic stations.

Variational characteristics of shear-wave splitting on the 2001 Shidian earthquakes in Yunnan, China

In 2001 three earthquakes occurred in Shidian in Yunnan Province, which were the MS=5.2 on April 10, the MS=5.9 on April 12 and the MS=5.3 on June 8. Based on the data from the station Baoshan of Yunnan Telemetry Digital Seismograph Network, the variational characteristics of shear-wave splitting on these series of strong earthquakes has been studied by using the systematic analysis method (SAM) of shear-wave splitting. The result shows the time delays of shear-wave splitting basically increase with earthquake activity intensifying. However the time delays abruptly decrease immediately before strong aftershocks. It accords with the stress relaxation before earthquakes, which was found recently in study on shear-wave splitting. The result suggests it is significant for reducing the harm degree of earthquakes to develop the stress-forecasting on earthquake in strong active tectonic zones and economic developed regions or big cities under the danger of strong earthquakes.

Preliminary analysis of the shear-wave splitting observations from the Qiaojia seismic array

A total of 351 shear-wave splitting results at 25 stations were obtained from the seismic data recorded in period of January,2013 to December,2016,by a broadband seismic array deployed in the northern segment of Xiaojiang Fault Zone(n-XJFZ).Meanwhile,the stress field of the n-XJFZ was determined by inverting 140 focal mechanism solutions of the small earthquakes within the study area which were recorded in the same period.This determination confirmed a compressive stress in NW-SE orientation and an extensional stress in the NE-SW orientation,with little difference from those released by previous studies.The shear-wave splitting results show a spatial complexity in polarization orientation,different from one site to another.How-ever,the polarization orientations integrated for the subareas suggest that the fault trends seemingly played important roles.All the subareas bear two dominant orientations,N10°E and N90°E,both of which are different from the azimuths of the principal compressive stress,due to the fault distribution.The time delay averaged over the entire region is 4.56 ms/km,close to that of the upper boundary of the generally accepted interval worldwide but larger than those in most of the investigated regions in the Chinese mainland,which probably implies an alignment of more micro-cracks in the n-XJFZ.Interestingly,the 2014 Ms6.5 Ludian earthquake was found to have caused a variation in the time delays of the slow shear waves within the study area though its epicenter was outside.This earthquake resulted in an evident drop of the time delays remaining for 4 months,however,lifted a bit the time-delay level with respect to that prior to the earthquake.

Systematic Analysis Method of Shear-Wave Splitting:SAM Software System

In order to make a more effective use of the data from regional digital seismograph networks and to promote the study on shear wave splitting and its application to earthquake stress-forecasting, SAM software system, i.e., the software on systematic analysis method of shear wave splitting has been developed. This paper introduces the design aims, system structure, function and characteristics about the SAM software system and shows some graphical interfaces of data input and result output. Lastly, it discusses preliminarily the study of shear wave splitting and its application to earthquake forecasting.

The Character of Shear-wave Splitting in Marble in the Critical State of Rupture

This paper mainly observed and analyzed the character of shear-wave splitting in rock specimens while they were in the critical state of rupture. The rock specimens for study are made of Laizhou marble from Shandong, China. A series of records were obtained from two rock specimens when they were in the critical state of rupture. The result shows that, in the critical state just before rock rupture, there may be the phenomenon of rise and fall in the time delay of shear-wave splitting, even though the load was kept constant. That is to say, the time delay of shear-wave splitting may have a falling process before rock rupture.

Numerical simulation of influence of material parameters on splitting spinning of aluminum alloy

During the splitting spinning process, the material parameters of disk blank have a significant effect on the determination of forming parameters and the quality of deformed blank. The influence laws of material parameters, including yield stress, hardening exponent and elastic modulus, on splitting spinning force, splitting spinning moment, degree of inhomogeneous deformation and quality of flange (average thickness and average deviation angle) were investigated by 3D-FE numerical simulation based on elasto-plastic dynamic explicit FEM under ABAQUS/Explicit environment. The results show that, the splitting spinning force and the splitting spinning moment increase with the increase of yield stress, hardening exponent and elastic modulus. The degree of inhomogeneous deformation increases with the decrease of yield stress and hardening exponent and the increase of elastic modulus. The average thickness of flange increases with the decrease of yield stress and the increase of hardening exponent and elastic modulus. The average deviation angle of upper surface increases with the increase of yield stress and the decrease of hardening exponent and elastic modulus. The average deviation angle of lower surface increases with the decrease of yield stress, hardening exponent and elastic modulus. Meanwhile, the corresponding variation ranges are given. The achievements may serve as an important guide for selecting the reasonable processing parameters of splitting spinning based on different aluminum alloys, and are very significant for optimum design and precision control of the splitting spinning process.

A seismological evidence for the northwestward movement of Africa with respect to Iberia from shear-wave splitting

Seismic anisotropy and its main features along the convergent boundary between Africa and Iberia are detected through the analysis of teleseismic shear-wave splitting. Waveform data generated by 95 teleseismic events recorded at 17 broadband stations deployed in the western Mediterranean region are used in the present study. Although the station coverage is not uniform in the Iberian Peninsula and north- west Africa, significant variations in the fast polarization directions and delay times are observed at stations located at different tectonic domains. Fast polarization directions are oriented predominantly NW-SE at most stations which are close to the plate boundary and in central Iberia; being consistent with the absolute plate motion in the region. In the northern part of the Iberian Peninsula, fast velocity direc- tions are oriented nearly E--W; coincident with previous results. Few stations located slightly north of the plate boundary and to the southeast of lberia show E--W to NE-SW fast velocity directions, which may be related to the Alpine Orogeny and the extension direction in lberia. Delay times vary significantly between 0.2 and 1.9 s for individual measurements, reflecting a highly anisotropic structure beneath the recording stations. The relative motion between Africa and lberia represents the main reason for the observed NW-SE orientations of the fast velocity directions. However, different causes of anisotropy have also to be considered to explain the wide range of the splitting pattern observed in the western Mediterranean region. Many geophysical observations such as the low Pn velocity, lower lithospheric Q values, higher heat flow and the presence of high conductive features support the mantle flow in the western Mediterranean, which may contribute and even modify the splitting pattern beneath the studied region.

Zero-field splitting parameters and local structures for tetragonal Cr^(2+) centers in Cr^(2+)-doped ZnSe semiconductors

The inter-relation between zero-field splitting （ZFS） parameters and local lattice structures of the （CrSe4）6 clusters in ZnSe semiconductors has been established by using the complete diagonalization （of the energy matrix） method. On the basis of this, the local lattice distortions, the ZFS parameters D, a, F and the optical spectrum for Cr2＋ ions doped into ZnSe are theoretically investigated, and the contributions of the spin singlets have been taken into account. The calculated ZFS parameters are in good agreement with the experimental values. From our calculations, the tetragonal distortion parameters AR = 0.091A and Aθ = 4.28° of Cr2＋ in ZnSe are acquired, and the results suggest that there exists a tetragonal expansion distortion for the local lattice structure of （CrSe4）6- clusters in ZnSe crystals. The influence of the spin singlets on ZFS parameters is also discussed, indicating that the contributions to ZFS parameters a and F cannot be ignored.

Shear-wave Splitting in the Crust beneath Shandong and Its Adjacent Area

Based on the seismic data of near-field source from the Shandong Digital Seismic Network,this study obtained the shear-wave splitting parameters from multiple stations. The results show that dominant directions of polarizations of fast shear-waves reflect the spatial distribution characteristics of tectonic stress in this area,CHD and LIS stations show 2 dominant directions,and reveals that the crustal seismic anisotropy in the Shandong area is constrained by the regional stress background,local failure and structure.

The Temporal Evolution of Shear-wave Splitting and Geoelectrical Anisotropy During Yongdeng Earthquake

The temporal evolution of shear-wave splitting and geoelectrical anisotropy of earth media during the preparation and occurrence of the Yongdeng MS5.8 earthquake is studied based on the digital seismic data and continuous magnetotelluric (MT) records. The results show that the direction of polarization of the fast S-wave gradually rotated from north by east to north by west before the Yongdeng earthquake and returned to north by east after the earthquake. Moreover, the time delay between the fast S wave and slow S wave increased to the largest until the earthquake occurred. On the other hand, the directions of geoelectrical principal axes also changed before earthquake, and exhibited identical characteristics in the frequency domain. The direction of the axis (pxy), which was originally perpendicular to the fault near station, varied from N17°to N15°before earthquake, and returned to north by east after the earthquake. Correspondingly, the change of pxy was most obvious in the period range of 160 s to 226 s, but the variation of pyx was less obvious. The focal mechanism solution of the earthquake showed that the direction of P-axis was N15°. Conclusions are: (1) The fast S-wave polarization, one of the geoelectrical principal axes and the P-axis were in agreement in orientation during earthquake preparation, which is the result of the alignment of EDA cracks; (2) The fact that the geoelectrical axes identically varied with frequency demonstrated that EDA is really widespread in the crust; (3) The variation of MT apparent resistivity also showed the existence of anisotropy and has its own features: The static anisotropy of resistivity is controlled by the tectonic conditions of the station, while the dynamic anisotropy is caused and affected by the seismic stress field; The resistivity along the axis parallel to the P-axis showed the most obvious variation, while the difference in variation with frequency shows that the variation of resistivity depends not only on direction but also on buried depth.

Nonlinear Algebraic Equations Solved by an Optimal Splitting-Linearizing Iterative Method

How to accelerate the convergence speed and avoid computing the inversion of a Jacobian matrix is important in the solution of nonlinear algebraic equations(NAEs).This paper develops an approach with a splitting-linearizing technique based on the nonlinear term to reduce the effect of the nonlinear terms.We decompose the nonlinear terms in the NAEs through a splitting parameter and then linearize the NAEs around the values at the previous step to a linear system.Through the maximal orthogonal projection concept,to minimize a merit function within a selected interval of splitting parameters,the optimal parameters can be quickly determined.In each step,a linear system is solved by the Gaussian elimination method,and the whole iteration procedure is convergent very fast.Several numerical tests show the high performance of the optimal split-linearization iterative method(OSLIM).

Multi-Modal Medical Image Fusion Based on Improved Parameter Adaptive PCNN and Latent Low-Rank Representation

Multimodal medical image fusion can help physicians provide more accurate treatment plans for patients, as unimodal images provide limited valid information. To address the insufficient ability of traditional medical image fusion solutions to protect image details and significant information, a new multimodality medical image fusion method(NSST-PAPCNNLatLRR) is proposed in this paper. Firstly, the high and low-frequency sub-band coefficients are obtained by decomposing the source image using NSST. Then, the latent low-rank representation algorithm is used to process the low-frequency sub-band coefficients;An improved PAPCNN algorithm is also proposed for the fusion of high-frequency sub-band coefficients. The improved PAPCNN model was based on the automatic setting of the parameters, and the optimal method was configured for the time decay factor αe. The experimental results show that, in comparison with the five mainstream fusion algorithms, the new algorithm has significantly improved the visual effect over the comparison algorithm,enhanced the ability to characterize important information in images, and further improved the ability to protect the detailed information;the new algorithm has achieved at least four firsts in six objective indexes.

Effect of Lanthanum on Rice Growth and Phy siological Parameters with Split-Root Nutrient Solution Culture

Split-root solution culture was used to study the promoting effect of lanthanum on rice (Oryza sativa) growth and its physiological mechanisms. Results sho w that low concentration (0.05～1.5 mg·L -1) increases rice yield an d grain numbers. High concentration depresses grain formation (9～30 mg·L -1 ) and root elongation (1.5～30 mg·L -1). No significant influence on str aw dry weight was found over the whole concentration range except the 0.05 mg·L -1 treatment. With the increase of La concentration from 0.05 to 0.75 mg· L -1, catalase (CAT) activity in the first fully expandeing leaves and root s decreases. When La concentration is greater than 0.75 mg·L -1 or less than 9 mg·L -1, it significantly decreases superoxide dismutase activity ( SOD) in the leaves and roots. No significant effects were found on chlorophyll, protein and malondialdehyde (MDA) content. Possible mechanisms of La′s promotin g effect on rice growth and reduction effect of ·O- 2 were discussed.

Upper mantle anisotropy inferred from shear wave splitting beneath the Eastern Indian Shield region

We estimate the shear wave splitting parameters vis-a-vis the thicknesses of the continental lithosphere beneath the two permanent seismic broadband stations located at Dhanbad(DHN) and Bokaro(BOKR) in the Eastern Indian Shield region. Broadband seismic data of 146 and 131 teleseismic earthquake events recorded at DHN and BOKR stations during 2007-2014 were analyzed for the present measurements.The study is carried out using rotation-correlation and transverse component minimization methods. We retain our "Good", "Fair" and "Null" measurements, and estimate the splitting parameters using 13 "Good" results for DHN and 10 "Good" results for BOKR stations. The average splitting parameters(φ,δt)for DHN and BOKR stations are found to be 50.76°±5.46° and 0.82±0.2 s and 56.30°±5.07° and0.95 ± 0.17 s, and the estimated average thicknesses of the anisotropic layers beneath these two stations are ～94 and ～109 km,respectively. The measured deviation of azimuth of the fast axis direction(0)from the absolute motion of the Indian plate ranges from ～8° to 14°. The measured deviation of azimuth of the fast axis direction(0) from the absolute motion of the Indian plate ranges from ～8° to 14°. The eastward deviation of the fast axis azimuths from absolute plate motion direction is interpreted to be caused by induced outflow from the asthenosphere. Further, the delay time found in the present analysis is close to the global average for continental shield areas, and also coherent with other studies for Indian shield regions. The five "Null" results and the lower delay time of ～0.5-0.6 s might be indicating multilayer anisotropy existing in the mantle lithosphere beneath the study area.

Experimental study on tensile failure process of double-K fracture parameters in roller compacted concrete layer

Study on failure of soft stratum of roller compacted concrete (RCC) is an important aspect of stability of high RCC dam. Six kinds of specimens with different interfaces were investigated by wedge splitting method. Double-K fracture parameters (initial fracture parameter and unsteady fracture parameter) were calculated by the concrete double-K fracture theory. It is indicated that the approach of construction joint or old joint after RCC final set is the most efficient among the six cases, and its fracture parameter is the largest among them. Its propagation path is sinuous. Its failure surface is scraggly. Bedding plane crack fails at the underside of the concrete surface, bond course or the surface between them for each approach. So disturbance on the underside of the concrete surface should be avoided or decreased at best during RCC construction.

Study on Correlation and Quantiative Error Estimation Method Among the Splitting Shear Wave Identification Methods

Based on the propagation characteristics of shear wave in the anisotropic layers,thecorrelation among several splitting shear-wave identification methods hasbeen studied.Thispaper puts forward the method estimating splitting shear-wave phases and its reliability byusing of the assumption that variance of noise and useful signal data obey normaldistribution.To check the validity of new method,the identification results and errorestimation corresponding to 95% confidence level by analyzing simulation signals have beengiven.

Study on the S-wave Splitting Characteristics in the Longtan Reservoir Area,Guangxi,China

In this paper,the shear wave splitting features of the Longtan reservoir area are studied by adopting the traditional cross-correlation coefficient method and polarization analysis,using the data recorded by the seismic network founded by a project under the National Science and Technology Pillar Program from April 2009 to April 2010.We found that most of polarization directions at seismic stations are consistent with the direction of the overall regional stress field,but local structures and faults may control or influence the fast shear-wave polarization direction.The time-delay normalized to source-station path is between 10 to 25ms/km,and among them,the time-delay is about 10ms/km at the LIL and XIL sites,which are farther away from the dam.The water depth is relatively shallow and seismic activity relatively weak after water storage,indicating the effect of reservoir water penetration,or loading,on the state of cracks in the reservoir area.We also found that the delay time changes consistently with the water level at stations DPD and GAL.It may be related to crack expansion and water penetration caused by the reservoir impoundment.

On Shear Wave Splitting in the Zhangjiakou-Bohai Sea Seismic Belt

The Zhangjiakou-Bohai Sea seismic belt is an important seismic zone in North China. The direction of principal compressive stress is near EW in this region. According to digital seismic data from the Capital Region from September 2005 to September 2010,and using the SAM method,the spatial distribution of the crust anisotropy characteristics are studied and discussed in the middle-eastern part of the Zhangjiakou-Bohai Sea seismic belt. The principal polarization direction,which is near EW direction,is obvious in the middleeastern Zhangjiakou-Bohai Sea seismic belt. The spatial distribution of polarization direction crossing the Zhangjiakou-Bohai Sea seismic belt shows that there is little difference among the Yanshan uplift area,inside of the seismic zone and North China basin,and the principal polarization direction is near EW.