Field survey around strong motion stations and its implications on the seismic intensity in the Lushan earthquake on April 20, 2013

The M _s7.0 Lushan earthquake on April 20, 2013 is another destructive event in China since the M _s8.0 Wenchuan earthquake in 2008 and M _s7.1 Yushu earthquake in 2010. A large number of strong motion recordings were accumulated by the National Strong Motion Observation Network System of China. The maximum peak ground acceleration (PGA) at Station 51BXD in Baoxing Country is recorded as ?1,005.3 cm/s², which is even larger than the maximum one in the Wenchuan earthquake. A field survey around three typical strong motion stations confirms that the earthquake damage is consistent with the issued map of macroseismic intensity. For the oscillation period 0.3–1.0 s which is the common natural period range of the Chinese civil building, a comparison shows that the observed response spectrums are considerably smaller than the designed values in the Chinese code and this could be one of the reasons that the macroseismic intensity is lower than what we expected despite the high amplitude of PGAs. The Housner spectral intensities from 16 stations are also basically correlated with their macroseismic intensities, and the empirical distribution of spectral intensities from Lushan and Wenchuan Earthquakes under the Chinese scale is almost identical with those under the European scale.     

Earthquake prediction on the basis of Vp/Vs variations — a case history

On November 7, 1976, an earthquake of the strike-slip fault type, and of magnitude 6.7, occurred at two mountainous localities, Yanyuan and Ninglang, in the border region of the Szechuan and Yunnan Provinces, China. One year before the earthquake, a prediction was made as to the magnitude and location of the impending earthquake by the present authors, on the basis of the results of a general survey of an area of about 20,000 km² for V_p/V_s ratio variations. The prediction of time of occurrence was made afterwards by the combined analysis of results of observations of some other precursory phenomena. The actual occurrence of the event was generally considered as being in agreement with the prediction.The present study may be taken as a new example for the detection of V_p/V_s ratio variations prior to an earthquake of the strike-slip fault type. By considering the difference between the shear rupture and stick-slip motion and the anisotropy induced by dilatancy, a preliminary discussion is made concerning the related results.     

Reconnaissance report and preliminary ground motion simulation of the 12 May 2008 Wenchuan earthquake

The M _w = 8.0 Wenchuan earthquake of May 12, 2008, caused destruction over a wide area. The earthquake cost more than 69,000 lives and the damage is reported to have left more than 5 million people homeless. It is estimated that 5.36 million buildings were destroyed and 21 million buildings were damaged in Sichuan and the nearby provinces. Economic losses due to the event are estimated to be 124 billion USD. From a field reconnaissance trip conducted in October 2008, it is evident that the combination of several factors, including mountainous landscape, strong ground shaking, extensive landslides and rock-falls, has exacerbated the human and economic consequences of this earthquake. Extensive damage occurred over a wide area due to the shear size of the earthquake rupture combined with poor quality building construction. In order to investigate the ground shaking during the earthquake, we have conducted a strong ground motion simulation study, applying a hybrid broadband frequency technique. The preliminary results show large spatial variation in the ground shaking, with the strongest ground motions along the fault plane. The simulation results have been calibrated against the recorded ground motion from several near-field stations in the area, and acceleration values of the order of 1 g are obtained, similar to what was recorded during the event. Comparison with the damage distribution observed in the field confirms that the effect of fault rupture complexity on the resulting ground motion distribution also controls to a large extent the damage distribution. The applied simulation technique provides a promising platform for predictive studies.     

Critical evaluation of strong motion in Kocaeli and Düzce (Turkey) earthquakes

Turkey was struck by two major events on August 17th and November 12th, 1999. Named Kocaeli (M_w=7.4) and Düzce (M_w=7.2) earthquakes, respectively, the two earthquakes provided the most extensive strong ground motion data set ever recorded in Turkey. The strong motion stations operated by the General Directorate of Disaster Affairs, the Kandilli Observatory and Earthquake Research Institute of Bogazici University and Istanbul Technical University have produced at least 27 strong motion records for the Kocaeli earthquake within 200 km of the fault. Kocaeli earthquake has generated six motions within 20 km of the fault adding significantly to the near-field database of ground motions for M_w>=7.0 strike–slip earthquakes. The paper discusses available strong motion data, studies their attenuation characteristics, analyses time domain, as well as spectral properties such as spectral accelerations with special emphasis on fault normal and fault parallel components and the elastic attenuation parameter, kappa. A simulation of the Kocaeli earthquake using code FINSIM is also presented.     

Long period microtremors, microseisms and earthquake damage: Northridge, CA, earthquake of 17 January 1994

Three studies of site amplification factors, based on the recorded aftershocks, and one study based on strong motion data, are compared one with another and with the observed distribution of damage from the Northridge, CA, earthquake of 17 January 1994 (M_L=6.4). In the epicentral area, when the peak ground velocities are larger than v_m≈15 cm/s, nonlinear response of soil begins to distort the amplification factors determined from small amplitude (linear) wave motion. Moving into the area of near-field and strong ground motion (v_m>30 cm/s), the site response becomes progressively more affected by the nonlinear soil response. Based on the published results, it is concluded that site amplification factors determined from small amplitude waves (aftershocks, small earthquakes, coda waves) and their transfer-function representation may be useful for small and distant earthquake motions, where soils and structures respond to earthquake waves in a linear manner. However in San Fernando Valley, during the Northridge earthquake, the observed distribution of damage did not correlate with site amplification determined from spectra of recorded weak motions. Mapping geographical distribution of site amplification using other than very strong motion data, therefore appears to be of little use for seismic hazard analyses.     

Strong ground motion from the November 12, 2017, <Emphasis Type="Bold">M</Emphasis> 7.3 Kermanshah earthquake in western Iran

In this paper, we analyzed the strong ground motion from the November 12, 2017, Kermanshah earthquake in western Iran with moment magnitude (M) of 7.3. Nonlinear and linear amplification of ground motion amplitudes were observed at stations with soft soil condition at hypocentral distances below and above 100 km, respectively. Observation of large ground motion amplitudes dominated with long-period pulses on the strike-normal component of the velocity time series suggests a right-lateral component of movement and propagation of rupture towards southeast. Comparison of the horizontal peak ground acceleration (PGA) from the M 7.3 earthquake with global PGA values showed a similar decay in ground motion amplitudes, although it seems that PGA from the M 7.3 Kermanshah earthquake is higher than global values for NEHRP site class B. We also found that the bracketed duration (D_b) was higher in the velocity domain than in the acceleration domain for the same modified Mercalli intensity (MMI) threshold. For example, D_b reached ~?30 s at the maximum PGA while it was ~?50 s at the maximum peak ground velocity above the threshold of MMI?=?5. Although the standard design spectrum from Iranian Code of Practice for Seismic Resistant Design of Buildings (standard No. 2800) seems to include appropriate values for the design of structures with fundamental period of 1 s and higher, it is underestimated for near-field ground motions at lower periods.     

Imaging the rupture process of the 10 January 2018 MW7.5 Swan island,Honduras earthquake*

The 10 January 2018 M_W7.5 Swan island, Honduras earthquake occurred on the Swan island fault, which is a transform plate boundary between the North American and Caribbean plates. Here we back-project the rupture process of the earthquake using dense seismic stations in Alaska, and find that the earthquake ruptured at least three faults (three stages) for a duration of ~40 s. The rupture speed for the longest fault (stage 3) is as fast as 5 km/s, which is much faster than the local shear wave velocity of ~4 km/s. Supershear rupture was incidentally observed on long and straight strike-slip faults. This study shows a supershear rupture that occured on a strike-slip fault with moderate length, implying that supershear rupture might commonly occur on large strike-slip earthquakes. The common occurrence of supershear rupture on strike-slip earthquakes will challenge present understanding of crack physics, as well as strong ground motion evaluation in earthquake engineering.     

Strong ground motion simulation of the 2003 Bam, Iran, earthquake using the empirical Green’s function method

The 2003 Bam, Iran, earthquake caused catastrophic damage to the city of Bam and neighboring villages. Given its magnitude (M _w ) of 6.5, the damage was remarkably large. Large-amplitude ground motions were recorded at the Bam accelerograph station in the center of Bam city by the Building and Housing Research Center (BHRC) of Iran. We simulated the Bam earthquake acceleration records at three BHRC strong-motion stations—Bam, Abaraq, and Mohammad-Abad—by the empirical Green’s function method. Three aftershocks were used as empirical Green’s functions. The frequency range of the empirical Green’s function simulations was 0.5–10 Hz. The size of the strong motion generation area of the mainshock was estimated to be 11 km in length by 7 km in width. To estimate the parameters of the strong motion generation area, we used 1D and 2D velocity structures across the fault and a combined source model. The empirical Green’s function method using a combination of aftershocks produced a source model that reproduced ground motions with the best fit to the observed waveforms. This may be attributed to the existence of two distinct rupture mechanisms in the strong motion generation area. We found that the rupture starting point for which the simulated waveforms best fit the observed ones was near the center of the strong motion generation area, which reproduced near-source ground motions in a broadband frequency range. The estimated strong motion generation area could explain the observed damaging ground motion at the Bam station. This suggests that estimating the source characteristics of the Bam earthquake is very important in understanding the causes of the earthquake damage.     

Source mechanism of the 22/11/1995 Gulf of Aqaba earthquake and its aftershock sequence

The Gulf of Aqaba earthquake occurred on 22/11/1995 with M _W = 7.2 and was the largest event to occur along the Dead Sea Transform in at least a century. We determined the source mechanism of the event based on the inversion of the P and SH waveforms observed by teleseismic stations. Our solution consists of 2 similar subevents, where the first exhibits a left-lateral strike-slip motion with a normal component and the second subevent shows an almost pure left-lateral strike-slip motion along the gulf major fault system. The total seismic moment is 7.7 × 10²⁶ dyne-cm, with source duration of 15 seconds. The aftershock sequence was recorded by the regional broad band stations BGIO and JER, of the Israel Seismic Network, and KEG, of the MEDNET regional network for earthquakeswith M _W > 4.0. Aftershock epicenter seems to cluster mainly in the Eilat basin and the Aragonese basin. The dominant mechanism indicates left-lateral strike slip motion along the gulf trend, similar to the mainshock.     

Estimation of strong ground motions from hypothetical earthquakes on the Cascadia subduction zone,Pacific Northwest

Strong ground motions are estimated for the Pacific Northwest assuming that large shallow earthquakes, similar to those experienced in southern Chile, southwestern Japan, and Colombia, may also occur on the Cascadia subduction zone. Fifty-six strong motion recordings for twenty-five subduction earthquakes ofM _s7.0 are used to estimate the response spectra that may result from earthquakesM _w<81/4. Large variations in observed ground motion levels are noted for a given site distance and earthquake magnitude. When compared with motions that have been observed in the western United States, large subduction zone earthquakes produce relatively large ground motions at surprisingly large distances. An earthquake similar to the 22 May 1960 Chilean earthquake (M _w 9.5) is the largest event that is considered to be plausible for the Cascadia subduction zone. This event has a moment which is two orders of magnitude larger than the largest earthquake for which we have strong motion records. The empirical Green's function technique is used to synthesize strong ground motions for such giant earthquakes. Observed teleseismicP-waveforms from giant earthquakes are also modeled using the empirical Green's function technique in order to constrain model parameters. The teleseismic modeling in the period range of 1.0 to 50 sec strongly suggests that fewer Green's functions should be randomly summed than is required to match the long-period moments of giant earthquakes. It appears that a large portion of the moment associated with giant earthquakes occurs at very long periods that are outside the frequency band of interest for strong ground motions. Nevertheless, the occurrence of a giant earthquake in the Pacific Northwest may produce quite strong shaking over a very large region.     

Simulation of strong earthquake characteristics of a scenario earthquake (MS7.5) based on the enlightenment of 2022 MS6.9 earthquake in Menyuan

The Menyuan area is an important transportation hub in the Hexi Corridor. The Menyuan M_S6.9 earthquake that occurred on January 8, 2022 had a major impact on the local infrastructure and transportation of this region. Due to the high possibility of similar strong earthquakes occurring in this area in the future, preliminary assessment of the seismic intensity characteristics of destructive earthquakes in this region is essential for effective disaster control. This paper uses the empirical Green′s function (EGF) method as a numerical simulation tool to predict the ground motion intensity of Datong Autonomous County under the action of the scenario earthquake (M_S7.5). Seismic records of aftershocks of the 2016 Menyuan M_S6.4 earthquake were used as Green’s functions for this simulation. The uncertainties associated with various source parameters were considered, and 36 possible earthquake scenarios were simulated to obtain 72 sets of horizontal ground motions in Datong County. The obtained peak ground acceleration (PGA) vs. time histories of the horizontal ground motion were screened using the attenuation relationships provided by the fifth-edition of China's Seismic Ground Motion Parameter Zoning Map and the NGA-West2 dataset. Ultimately, 32 possible acceleration-time histories were selected for further analysis. The screened PGA values ranged from 78.8 to 153 cm/s². The uncertainty associated with the initial rupture point was found to greatly affect the results of the earthquake simulation. The average acceleration spectrum of the selected acceleration-time history exceeded the expected spectrum of a intermediate earthquake, which means that buildings in Datong County might sustain some damage should the scenario earthquake occur. This research can provide reliable ground motion input for urban earthquake damage simulation and seismic design in Datong County. Growing the dataset of small earthquakes recorded in this region will facilitate the large-scale simulation of ground motions under different earthquake scenarios.     

Rapid report of June 1, 2022 MW 5.9 Lushan earthquake,China with geodetic and teleseismic data

Timely response to earthquake characterization can facilitate earthquake emergency rescue and further scientific investigations. On June 1, 2022, M_W 5.9 earthquake occurred in the southern area of the Longmenshan fault zone. This event also happened at the south end of the Dayi seismic gap and is the largest earthquake that has occurred in this seismic gap since the 1970 M 6.2 event. The slip-distribution model constrained by the seismic waveforms suggests a thrust-dominated faulting mechanism. The main slip occurs at a depth of ～14 ?km, and the cumulative energy is released in the first 6 ?s. The variations of Coulomb stress caused by the mainshock show a positive change in the southwest area of the Dayi seismic gap, indicating possible activation of future earthquakes. In addition, we emphasize the importance of rapid estimation of deformation for near-field hazard delineation, especially when interferometric radar fails to image coseismic deformation in a high relief terrain.     

Strong motion characteristics of the <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> 6.6 Lushan earthquake,Sichuan, China — an insight into the spatial difference of a typical thrust fault earthquake

Near-field strong ground motions are useful for engineering seismology studies and seismic design, but dense observation networks of damaging earthquakes are still rare. In this study, based on the strong-motion data from the M w 6.6 Lushan earthquake, the ground motion parameters in different spatial regions are systematically analyzed, and the contributions from different effects, like the hanging-wall effect, directivity effect, and attenuation effect are separated to the extent possible. Different engineering parameters from the observed ground motions are compared with the local design response spectra and a new attenuation relation of Western China. General results indicate that the high frequency ground motion, like the peak ground acceleration, on two sides of the fault plane is sensitive to the hanging-wall effect, whereas the low frequency ground motion, like the long period spectral acceleration, in the rupture propagation direction is affected by the directivity effect. Moreover, although the M w 6.6 Lushan earthquake is not a large magnitude event, the spatial difference of ground motion is still obvious; thus, for a thrust faulting earthquake, in addition to the hanging effect, the directivity effect should also be considered.     

基于随机振动模型的强地面运动数值模拟及其在广西灵山M6(3/4)地震中的应用

1936年4月1日广西灵山县平山镇东南罗阳山附近发生M6(3/4)地震,该地震是华南大陆自有地震记载以来发生的最大地震。本研究收集整理了灵山M6(3/4)地震的地质资料、活动断层探测最新成果等,选取适当的研究区域,利用随机振动有限断层模型计算区域内网格点的峰值加速度(PGA)及峰值速度(PGV)等参数,并且加入浅层横波速度结构V_s³⁰对地震动参数的影响,最终得到此次地震的地震动分布并分析了地震动特征。本研究将模拟结果与野外调查烈度数据和地震动衰减关系进行对比,结果显示模拟结果与调查烈度值和衰减关系在整体衰减特征、极震区的分布等方面均符合较好,模拟结果可为该地区未来地震危害性评估提供依据。本研究所使用的方法流程亦可应用于本地区地震烈度速报,为震后应急救援及决策指挥提供帮助。     

Modeling the strong ground motion and rupture characteristics of the March 31, 2006, Darb-e-Astane earthquake,Iran, using a hybrid of near-field SH-wave and empirical Green’s function method

We analyze the strong motion accelerograms of the moderate (M _w = 6.1), March 31, 2006, Darb-e-Astane earthquake of western Iran and also those of one of its prominently recorded, large (M _w = 5.1) foreshock and (M _w = 4.9) aftershock. (1) Using derived SH-wave spectral data, we first objectively estimate the parameters W _o\mathit{\Omega} _{\rm o} (long period spectral level), f _c (corner frequency) and Q(f) (frequency dependent, average shear wave quality factor), appropriate for the best-fit Brune ω ^− 2 spectrum of each of these three events. We then perform a non-linear least square analysis of the SH-wave spectral data to provide approximate near-field estimates of the strike, dip, and rake of the causative faults and also the seismic moment, moment magnitude, source size, and average stress drop of these three events. (2) In the next step, we use these approximate values and an empirical Green’s function approach, in an iterative manner, to optimally model the strong ground motion and rupture characteristics of the main event in terms of peak ground acceleration/velocity/displacement and duration of ground shaking and thereby provide improved, more reliable estimates of the causative fault parameters of the main event and its asperities. Our near-field estimates for both the main moderate event and the two smaller events are in good conformity with the corresponding far-field estimates reported by other studies.     

Earthquake strong motions recorded by a large near-source array of digital seismographs

Specially designed arrays of strong-motion seismographs near to the earthquake source are required for seismological and engineering studies of the generation and near-field properties of seismic waves. The first such large digital array, called SMART 1 (with radius 2 km and 37 accelerometers), to record substantial ground motion (up to 0·24g horizontal acceleration) became operational in late 1980 in a highly seismic region of Taiwan. During the first 6 months of operation, SMART 1 recorded nine earthquakes with magnitudes ranging from M_L 3·8 to 6·9. Three were located directly below the array at focal depths of 59 to 76 km. The remaining six had shallow depths and epicentral distances from 7 to 193 km. Digital records from 27 three-component accelerographs were obtained from a magnitude 6·9 (M_L) local earthquake on 29 January 1981. Representative measurements are described of seismic wave coherency and power spectrum as a function of wave number, frequency, azimuth of propagation and wave type. Acceleration waveforms varied significantly across the array for each event. On average, peak acceleration of horizontal components was about three times that of the vertical component. Relative spectral changes from earthquake to earthquake were large.     

Earthquake ground motion predictive equations for Garhwal Himalaya,India

Predictive equations based on the stochastic approach are developed for earthquake ground motions from Garhwal Himalayan earthquakes of 3.5≤M_w≤6.8 at a distance of 10≤R≤250 km. The predicted ground motion parameters are response spectral values at frequencies from 0.25 to 20 Hz, and peak ground acceleration (PGA). The ground motion prediction equations (GMPEs) are derived from an empirically based stochastic ground motion model. The GMPEs show a fair agreement with the empirically developed ground motion equations from Himalaya as well as the NGA equation. The proposed relations also reasonably predict the observed ground motion of two major Himalayan earthquakes from Garhwal Himalayan region. For high magnitudes, there is insufficient data to satisfactorily judge the relationship; however it reasonably predicts the 1991 Uttarkashi earthquake (M_w=6.8) and 1999 Chamoli earthquake (M_w=6.4) from Garhwal Himalaya region.     

Topographic site effects of Xishan Park ridge in Zigong city,Sichuan considering epicentral distance

Ground motion amplification induced by topography plays a vital role in engineering seismology. A topographic array of 8 accelerographs has been operating along the ridge in Xishan Park since 2007. The topographic site effects in Zigong city are studied based on the strong motion data of 2008 M_s 8.0 Wenchuan earthquake (the epicentral distance?=?225 km) and 2019 M_s 5.2 Zizhong earthquake (the epicentral distance?=?29 km). We compare the peak ground acceleration (PGA) of the two earthquakes and find that the PGA of Station 7#, which locates on a relatively steep slope, is amplified by 4.41 times comparing with the reference station in Zizhong earthquake, while this value is only 1.62 in Wenchuan earthquake. Fourier amplitude spectrum shows that the high frequency content of Zizhong earthquake is more abundant because of its smaller epicentral distance. By using the standard spectral ratio (SSR) method, we conclude that the amplification occurs because high-frequency ground motion is likely to resonate at small-scale features. Finally, the 3D numerical simulations are used to verify these conclusions. Our work indicates that more sophisticated numerical models need to be established for more accurate topographic site effects quantification. In addition, the influence of nearby topographic features should be considered when selecting reference stations.