Pulse-like ground motion observed during the 6 February 2023 M_(W)7.8 Pazarcık Earthquake(Kahramanmaraş,SE Türkiye)

In this study,we analyzed 100 three-component strong ground motion records observed within 200 km of the causative fault of the 6 February 2023 M_(W)7.8 Pazarcık(Kahramanmaraş)Earthquake in SE Türkiye.The wavelet method was utilized to identify and analyze the characteristics of pulse-like ground motions in the near-fault region,while considering the uncertainty of the pulse orientation during the analysis.Our investigation focused on the effects of the focal mechanism and rupture process on the spatial distribution,pulse orientation,and maximum pulse direction of the observed pulse-like ground motion.We also analyzed the amplitude and period of the observed ground pulses and the effect of long-period amplification on the ground motion response spectra.Our results indicated the following:(1)A total of 21 typical ground velocity pulses were observed during this earthquake,exhibiting complex characteristics due to the influence of the strike-slip mechanism and rupture directivity.Most ground pulses(17 out of 21)were recorded within 20 km of the fault,in a wide range of orientations,including normal and parallel to the fault direction.The waveforms exhibited unidirectional features,indicating the effects of left-lateral fault slip.Distinct pulses observed more than 20 km from the fault were mainly oriented normal to the fault.The waveforms were bidirectional with double-or multi-round trips as a result of rupture directivity.(2)The amplitudes of the observed pulses ranged from 30.5 to 220.0 cm/s,with the largest peak velocity of 220.0 cm/s observed at Station 3138.The pulse periods ranged from 2.3 to 14.5 s,with the longest pulse period of 14.5 s observed at Station 3116.The amplitude and period of the pulses observed during this earthquake were comparable to those of similar-magnitude global earthquakes.The amplitude of the pulses decreased significantly with increasing fault distance,whereas the pulse period was not significantly affected by the fault distance.(3)Compared with non-pulse records,the velocity pulse records had a pronounced amplification effect on the acceleration response spectra near the pulse period,with factors ranging from 2.1 to 5.8.The larger velocity pulses also significantly amplified the velocity response spectra,particularly over the long periods.This significant amplification effect of the pulses on the response spectra leads to empirical models underestimating the long-period earthquake ground motion.

Site dependence of far-source ground motions during the Wenchuan earthquake

This paper aimed to examine the site dependence and evaluate the methods for site analysis of far-source ground motions. This was achieved through the examination of frequency content estimated by different methods based on strong ground motions recorded at twelve far-source stations in Shandong province during the Wenchuan earthquake. The stations were located in sites with soil profiles ranging from code classes Ⅰ to Ⅲ. Approaches used included the Fourier amplitude spectrum （FAS）, the earthquake response spectrum （ERS）, the spectral ratio between the horizontal and the vertical components （H/V）, the spectral ratio between the spectra at the site and at a reference site （SRRS）, and coda wave analysis （CWA）. Results showed that major periods of these ground motions obtained by FAS, ERS and H/V ratio methods were all evidently larger than site dominant periods; the periods were also different from each other and mainly reflected the frequency content of long period components. Prominent periods obtained by the SRRS approach neither illuminated the long period aspect nor efficiently determined site features of the motions. The CWA resulted in a period close to site period for stations with good quality recordings. The results obtained in this study will be useful for the evaluation of far-source effect in constructing seismic design spectra and in selecting methods for ground motion site analysis.

Fragility of disconnect switch subjected to random earthquake ground motions

A fragility calculation scheme is estabtished in this paper for porcelain-type equipments subjected to random earthquake ground motions. All steps of the method are illustrated by the seismic damage analysis of GW4-110 disconnect switch. The model of the equipment is built applying the finite element method with flexible joints, and the seismic response of the equipment is analyzed using elastic time history method. On the base, according to the strength damage index and Monte-Carlo Method, the seismic damage ratios are counted and the seismic fragility curves are presented. Then the seismic damage of GW4-110 disconnect switch can be predicted.

Instantaneous spectrum estimation of earthquake ground motions based on unscented Kalman filter method

Representing earthquake ground： motion as time varying ARMA model, the instantaneous spectrum can only be determined by the time varying coefficients of the corresponding ARMA model. In this paper, unscented Kalman filter is applied to estimate the time varying coefficients. The comparison between the estimation results of unscented Kalman filter and Kalman filter methods shows that unscented Kalman filter can more precisely represent the distribution of the spectral peaks in time-frequency plane than Kalman filter, and its time and frequency resolution is finer which ensures its better ability to track the local properties of earthquake ground motions and to identify the systems with nonlinearity or abruptness. Moreover, the estimation results of ARMA models with different orders indicate that the theoretical frequency resolving power of ARMA model which was usually ignored in former studies has great effect on the estimation precision of instantaneous spectrum and it should be taken as one of the key factors in order selection of ARMA model.

Prediction of near-field strong ground motions for scenario earthquakes on active fault

A method to predict near-field strong ground motions for scenario earthquakes on active faults is proposed. First, macro-source parameters characterizing the entire source area, i.e., global source parameters, including fault length, fault width, rupture area, average slip on the fault plane, etc., are estimated by seismogeology survey, seismicity and seismic scaling laws. Second, slip distributions characterizing heterogeneity or roughness on the fault plane, i.e., local source parameters, are reproduced/evaluated by the hybrid slip model. Finally, the finite fault source model, developed from both the global and local source parameters, is combined with the stochastically synthetic technique of ground motion using the dynamic comer frequency based on seismology. The proposed method is applied to simulate the acceleration time histories on three base-rock stations during the 1994 Northridge earthquake. Comparisons between the predicted and recorded acceleration time histories show that the method is feasible and practicable.

A modified method for simulating non-stationary multi-point earthquake ground motion

A spectral-representation-based algorithm is proposed to simulate non-stationary and stochastic processes with evolutionary power,according to a prescribed non-stationary cross-spectral density matrix. Non-stationary multi-point seismic ground motions at different locations on the ground surface are generated for use in engineering applications. First,a modified iterative procedure is used to generate uniformly modulated non-stationary ground motion time histories which are compatible with the prescribed power spectrum. Then,ground motion time histories are modeled as a non-stationary stochastic process with amplitude and frequency modulation. The characteristic frequency and damping ratio of the Clough-Penzien acceleration spectrum are considered as a function of time in order to study the frequency time variation. Finally,two numerical examples are presented to validate the efficiency of the proposed method,and the results show that this method can be effectively applied to the dynamic seismic analysis of long and large scale structures.

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.

Effects of the probability of pulse-like ground motions on landslide susceptibility assessment in near-fault areas

Earthquake-induced strong near-fault ground motion is typically accompanied by largevelocity pulse-like component,which causes serious damage to slopes and buildings.Although not all near-fault ground motions contain a pulse-like component,it is important to consider this factor in regional earthquake-induced landslide susceptibility assessment.In the present study,we considered the probability of the observed pulse-like ground motion at each site(PP)in the region of an earthquake as one of the conditioning factors for landslide susceptibility assessment.A subset of the area affected by the 1994Mw6.7 Northridge earthquake in California was examined.To explore and verify the effects of PP on landslide susceptibility assessment,seven models were established,consisting of six identical influencing factors(elevation,slope gradient,aspect,distance to drainage,distance to roads,and geology)and one or two factors characterizing the intensity of the earthquake(distance to fault,peak ground acceleration,peak ground velocity,and PP)in logistic regression analysis.The results showed that the model considering PP performed better in susceptibility assessment,with an area under the receiver operating characteristic curve value of 0.956.Based on the results of relative importance analysis,the contribution of the PP value to earthquakeinduced landslide susceptibility was ranked fourth after the slope gradient,elevation,and lithology.The prediction performance of the model considering the pulse-like effect was better than that reported previously.A logistic regression model that considers the pulse-like effect can be applied in disaster prevention,mitigation,and construction planning in near-fault areas.

Influence of earthquake ground motion incoherency on multi-support structures

A linear response history analysis method is used to determine the influence of three factors:geometric incoherency,wave-passage,and local site characteristics on the response of lnulti-support structures subjected to differential ground motions.A one-span frame and a reduced model of a 24-span bridge,located in Las Vegas,Nevada are studied,in which the influence of each of the three factors and their combinations are analyzed.It is revealed that the incoherency of earthquake ground motion can have a dramatic influence on structural response by modifying the dynamics response to uniform excitation and inducing pseudo-static response,which does not exist in structures subjected to uniform excitation.The total response when all three sources of ground motion incoherency are included is generally larger than that of uniform excitation.

Ground motion attenuation of M_s8.0 Wenchuan earthquake

The great Ms8.0 Wenchuan earthquake has been the most destructive earthquake since 1949 in China. The earthquake occurred no more than half a year after the establishment of the National Strong Motion Observation Network System （NSMONS） of China; what is more, the epicenter was located in the area with dense strong motion observation sta- tions so that a large number of strong motion records of the main shock were obtained. In this paper, 501 strong motion re- cords from 167 observation stations are utilized to establish the ground motion attenuation relations in three directions in the range of fault distance less than 600 km. The result shows the difference of seismic motion attenuation in two horizontal di- rections is insignificant. It is the first time that strong-motion records are used to establish the ground motion attenuation relations of the Ms8.0 earthquake in China.

Preliminary research on the method for determination of ground motion parameters of diffuse earthquakes

This paper proposes that on the basis of near-source saturation characteristic of ground motion, the peak ground motion acceleration and the acceleration response spectra of diffuse earthquakes affecting the nuclear facility site can be determined by the statistic method. In view of the problem of insufficient statistic samples due to fewer qualified strong earth- quake records on the bedrock site, it also proposes that according to the features of ground motion response of the soil, the acceleration records on the soil site can be transformed into the equivalent acceleration records on rock site for a simultane- ous statistic analysis together with the records on rock site. Consequently, we have obtained by this method the near-source peak acceleration and the acceleration response spectra of small and moderate earthquakes on the site that could meet the need of nuclear facility and made a comparative analysis to the results from relevant researches. This new method for deter- mining diffuse earthquakes is only preliminarily verified in the paper and further research is needed for practical application.

Stochastic earthquake source model for ground motion simulation

In the analysis and design of important structures with relatively long life spans, there is a need to generate strong motion data for possible large events. The source of an earthquake is characterized by the spatial distribution of slip on the fault plane. For future events, this is unknown. In this paper, a stochastic earthquake source model is developed to address this issue. Here, 1D and 2D stochastic models for slip distribution developed by Lavallée et al.(2006) are used. The random field associated with the slip distribution is heavy-tailed stable distribution which can be used for large events. Using 236 past rupture models, the spectral scaling parameter and the four stable or Levy's parameters against empirical relationship for known quantities like magnitude or fault length are developed. The model is validated with data from 411 stations of 1999 Chi-Chi earthquake. The simulated response spectrum showed good agreement to actual data. Further the proposed model is used to generate ground motion for the 1993 Killari Earthquake where strong motion data is not available. The simulated mean peak ground velocity was in turn related to the intensity(MSK) and compared against values in the literature.

Ground-motion simulation for the MW6.1 Ludian earthquake on 3 August 2014 using the stochastic finite-fault method

The stochastic finite-fault simulation method was applied to synthesize the horizontal ground acceleration seismograms produced by the MW6.1 Ludian earthquake on August 3,2014.For this purpose,we produced first a total of 200 kinematic source models for the Ludian event,which are characterized by the heterogeneous slip on the conjugated ruptured fault and the slip-dependent spreading of the rupture front.The results indicated that the heterogeneous slip and the spatial extent of the ruptured fault play dominant roles in the spatial distribution of ground motions in the near-fault area.The peak ground accelerations(PGAs)and 5%-damped pseudospectral accelerations(PSAs)at periods shorter than 0.5 s estimated on the resulting synthetics generally match well with the observations at stations with Joyner-Boore distances(RJB)greater than 20 km.The synthetic PGVs and PSAs at periods of 0.5 s and 0.75 s are in good agreement with predicted medians by the Yu14 model(Yu et al.,2014).However,the synthetic results are generally much lower than the predicted medians by BSSA14 model(Boore et al.,2014).Moreover,the ground motion variability caused by the randomness in the source rupture process was evaluated by these synthetics.The standard deviations of PSAs on the base-10 logarithmic scale,Sigma[log10(PSA)],are closely dependent on either the spectral period or the RJB.The Sigma[log10(PSA)]remains a constant approximately 0.55 at periods shorter than 0.1 s,and then increase continuously up to^0.13 as the period increases from 0.1 to 2.0 s.The Sigma[log10(PSA)]values at periods of 0.1‒2.0 s show the downward tendency as the RJB values increase.However,the Sigma[log10(PSA)]​values at periods shorter than 0.1 s decrease as the RJB values increase up to^50 km,and then increase with the increasing RJB.Furthermore,we found that the ground-motion variability shows the significant dependence on the azimuth.

Rupture directivity effect on strong ground motion during the 12 May 2008 M_(W)7.9 Wenchuan earthquake

We focus here on the rupture directivity effect on the spatial distribution and attenuation characteristics of near-field ground motions during the 2008 MW7.9 Wenchuan earthquake.We examine the difference between the observed ground motions in and opposite the rupture directions and compare them with Next Generation Attenuation-West2(NGA-West2)ground motion prediction models.The isochrone directivity predictor is used to quantify the band-limited nature of the rupture directivity effect on strong ground motion.Our results show that the observed peak ground velocity(PGV)and spectral accelerations of periods of 1.0 s and longer are significantly amplified in the rupture direction,but de-amplified in the opposite direction affected by rupture directivity effect of this event.In contrast,the effect of rupture directivity on the observed peak ground acceleration(PGA)and periods of shorter than 1.0 s are relatively weak.The rupture directivity of this event shows clear period dependent and band limited characteristics with the strongest effect occurring around the period of 7.5 s.

Estimation of the 2010 Mentawai tsunami earthquake rupture process from joint inversion of teleseismic and strong ground motion data

Joint inversion of teleseismic body-wave data and strong ground motion waveforms was applied to determine the rupture process of the 2010 Mentawai earthquake. To obtain stable solutions, smoothing and non-negative constraints were introduced. A total of 33 teleseismic stations and 5 strong ground motion stations supplied data. The teleseismic and strong ground motion data were separately windowed for 150 s and 250 s and bandpass filtered with frequencies of 0.001e1.0 Hz and 0.005e0.5 Hz, respectively. The finitefault model was established with length and width of 190 km and 70 km, and the initial seismic source parameters were set by referring to centroid moment tensor(CMT) solutions. Joint inversion results indicate that the focal mechanism of this earthquake is thrust fault type, and the strike, dip, and rake angles are generally in accordance with CMT results. The seismic moment was determined as 5.814 1020Nm(Mw7.8) and source duration was about 102 s, which is greater than those of other earthquakes of similar magnitude. The rupture nucleated near the hypocenter and then propagated along the strike direction to the northwest, with a maximum slip of 3.9 m. Large uncertainties regarding the amount of slip retrieved using different inversion methods still exist; however, the conclusion that the majority of slip occurred far from the islands at very shallow depths was found to be robust. The 2010 Mentawai earthquake was categorized as a tsunami earthquake because of the long rupture duration and the generation of a tsunami much larger than was expected for an earthquake of its magnitude.

Orientation effect on ground motion measurement for Mexican subduction earthquakes

The existence of the principal directions of the ground motion based on Arias intensity is well-known. These principal directions do not necessarily coincide with the orientations of recording sensors or with the orientations along which the ground motion parameters such as the peak ground acceleration and the pseudo-spectral acceleration （PSA） are maximum. This is evidenced by the fact that the maximum PSA at different natural vibration periods for horizontal excitations do not correspond to the same orientation. A recent analysis carried out for California earthquake records suggests that an orientation-dependent ground motion measurement for horizontal excitations can be developed. The main objective of this study is to investigate and provide seismic ground motion measurements in the horizontal plane, including bidirectional horizontal ground motions, for Mexican interplate and inslab earthquake records. Extensive statistical analyses of PSA are conducted for the assessment. The analysis results suggest that similar to the case of California records, the average behavior of the ratio of the PSA to the maximum resulting PSA can be approximated by a quarter of an ellipse in one quadrant; and that the ratio can be considered to be independent of the value of the maximum resulting PSA, earthquake magnitude, earthquake distance and the focal depth. Sets of response ratios and attenuation relationships that can be used to represent a bidirectional horizontal ground motion measurement for Mexican interplate and inslab earthquakes were also developed.

Ground motion prediction equations based on shallow crustal earthquakes in Georgia and the surrounding Caucasus

Strong ground motions caused by earthquakes with magnitudes ranging from 3.5 to 6.9 and hypocentral distances of up to 300km were recorded by local broadband stations and three-component accelerograms within Georgia’s enhanced digital seismic network.Such data mixing is particularly effective in areas where strong ground motion data are lacking.The data were used to produce models based on ground-motion prediction equations(GMPEs),one benefit of which is that they take into consideration information from waveforms across a wide range of frequencies.In this study,models were developed to predict ground motions for peak ground acceleration and 5%-damped pseudo-absolute-acceleration spectra for periods between 0.01 and 10 s.Short-period ground motions decayed faster than long-period motions,though decay was still in the order of approximately 1/r.Faulting mechanisms and local soil conditions greatly influence GMPEs.The spectral acceleration(SA)of thrust faults was higher than that for either strike-slip or normal faults but the influence of strike-slip faulting on SA was slightly greater than that for normal faults.Soft soils also caused significantly more amplification than rocky sites.

Discussion on the Relationship between “Anti-great Earthquakes” and Ground Motion with a Low Probability of Exceedance

In this paper we discuss the differences and relationship between"great earthquakes",rare ground motion,and very rare ground motion.Taking the Beichuan-Yingxiu potential seismic source zone in Longmenshan seismic belt as an example,we revealed the relationship between the effects of"great earthquakes"and rare ground and very rare ground motion.After pointing out scientific and technical problems in the current seismic fortification system,we suggest that very rare ground motion should be considered if we want to deduce the potential hazard of great earthquakes in the future.

Variation of earthquake ground motionwith depth

The variation of ground motion with depth is investigated in this paper, which is essential in determining the seismic design ground motion for embedded structures and pipelines. Firstly the earthquake ground motion on surface and in subsoil recorded by Hosokura Mine array of Japan and six California Strong Motion Instrumentation Program geotechnical arrays of the United States from about fifty moderate and strong earthquakes are described. Then the arrays were classified into three different categories according to their site conditions. Finally the variation law of ground and sub-ground motion with depth is established as a result of a nonlinear regression analysis of the above-mentioned data. Through comparing the features in different sites, it is concluded that, in general, the amplitude (acceleration, velocity and displacement) of ground motion are decreasing with depth and the attenuation rate is higher in the shallow strata than that in deeper ones.

Simulation of Ground Motion at Guanting Reservoir Dam Based on the Scenario Earthquake in the Yanqing-Huailai Basin

Research at home and abroad shows that the simulation of ground motion using the 3D finite-difference method might be accurate and feasible. Based on related theories and methods,and using the wave velocity and density model of the crust in the Yanqing-Huailai Basin,this paper makes a simulation of ground motion at Guanting Reservoir Dam based on the scenario earthquake in the Yanqing-Huailai Basin. Comparative analysis shows that the results of 3D finite-difference simulation accord with those of the empirical formula. The parameters such as the velocity-time series of ground motion,PGV and frequency might be referred to for the analysis of seismic protection design of the dam's structure.