A simplified method for simulation of strong ground motion using finite rupture model of the earthquake source

We present a simplified method to simulate strong ground motion for a realistic representation of a finite earthquake source burried in a layered earth. This method is based on the stochastic simulation method of Boore (Boore, D. M., 1983, Bull. Seism. Soc. Am. 73, 1865–1894) and the Empirical Greens Function (EFG) method of Irikura (Irikura, K., 1986, Proceedings of the 7th Japan Earthquake symposium, pp. 151–156). The rupture responsible for an earthquake is represented by several subfaults. The geometry of subfaults and their number is decided by the similarity relationships. For simulation of ground motion using the stochastic simulation technique we used the shapping window based on the kinetic source model of the rupture plane. The shaping window deepens on the geometry of the earthquake source and the propagation characteristics of the energy released by various subfaults. The division of large fault into small subfaults and the method for accounting their contribution at the surface is identical to the EGF. The shapping window has been modified to take into account the effect of the transmission of energy released form the finite fault at various boundaries of the layered earth model above the source. In the present method we have applied the correction factor to adjust slip time function of small and large earthquakes. The correction factor is used to simulate strong motion records having basic spectral shape of ² source model in broad frequency range. To test this method we have used the strong motion data of the Geiyo earthquake of 24th March 2001, Japan recorded by KiK network. The source of this earthquake is modelled by a simple rectangular rupture of size 24 × 15 km, burried at a depth of 31 km in a multilayered earth model. This rupture plane is divided into 16 rectangular subfaults of size 6.0 × 3.75 km each. Strong motion records at eight selected near-field stations were simulated and compared with the observed records in terms of the acceleration and velocity records and their response spectrum. The comparison confirms the suitability of proposed rupture model responsible for this earthquake and the efficacy of the approach in predicting the strong motion scenario of earthquakes in the subduction zone. Using the same rupture model of the Geiyo earthquake, we compared the simulated records from our and the EGF techniques at one near-field station. The comparison shows that this technique gives records which matches in a wide frequency range and that too from simple and easily accessible parameters of burried rupture.     

Simulating strong ground motion with the `k -2 ' kinematic source model: An application to the seismic hazard in the Erzincan basin, Turkey

We present a seismic hazard application of a kinematic broad-bandrupture model. This model is based on the k-square dislocation distribution concept (Herrero and Bernard, 1994).Synthetic seismograms are calculated in the far-field approximation with alayered velocity medium for the 13 March 1992 Erzincan earthquake.With a parametrization of the source constrained by other studies,the far-field contribution correctly fits the recorded strong ground motion, which presents a 0.5 to 2 Hz dominant frequency range.As the k-square model is particularly well adapted to synthetize realistic strong-motion at short distances from the fault,it is a reliable tool for calculating seismic hazard maps around active faults. We thus present a synthetic peak ground acceleration map associatedwith the 13 March 1992 activated fault, for a 60 km × 60 km regionaround the epicenter taking into account a smoothed velocity structure ofthe basin in agreement with the absence of significant site effects related to1D resonance deduced from the aftershock records study.This map is compared with several post seismic reports: macroseismicintensities, detailed distribution of damage, and soil cracking andliquefaction. Our model shows that the values of the peak acceleration andvelocity can explain the dominant spatial distributionof these effects, which concentrates in a narrow band along theactivated segment fault, and in particular at its southernextremity. These results enable us to present such maps forhypothetical future earthquake ruptures, located on the major visible orinferred active fault segments in and around the basin. The effects ofthese potential sources are analyzed in relation to the 1992 eventeffects in order to eliminate unknown site responses. We show that thesouthern part of the basin is particularly exposed because of thepresence of strike-slip faults,and that the western part of the basin would suffer a significantlyhigher strong motion levelthan during the 1992 event with the activation of moderate sizednormal faults evidenced on thesouthwestern edge of the basin.     

2008年汶川地震近断层竖向与水平向地震动特征

选取分布在北川-映秀中央断裂两侧断层距120 km以内的40个强震动台站的记录,对汶川地震近断层地震动竖向和水平向加速度峰值、速度峰值、竖向和水平向加速度反应谱及谱比值进行了统计分析.研究表明:(1)地震动加速度峰值有显著的上盘效应,经验衰减模型的结果表明,在距地表破裂3~60 km的范围内,龙门山发震断层上盘一侧竖向与水平向的加速度峰值要比衰减模型得到的平均值大30%~40%.上盘的加速度峰值残差大部分是正值,而断层下盘残差大部分为负;水平地震动的东西分量幅值总体要大于南北分量,东西分量衰减相对较慢.(2)地震动长周期成分较弱,加速度反应谱值随周期增大而迅速减小,在周期1.0 s 时,即使在靠近中央断裂的最大加速度反应谱值也只有0.5 g;地震动加速度反应谱谱比值(竖向/水平向)沿龙门山断层周围的分布,在较长周期(T=0.2 s, 0.5 s, 1.0 s)与短周期(T=0.05 s, 0.1 s)有明显的不同.(3)近断层竖向地震动显著,地震动加速度峰值比在(竖向/水平向)可达1.4.在龙门山发震断层的上盘,地震动加速度峰值比整体上比下盘要大,竖向地震动尤为剧烈.部分近断层记录的地震动谱比值(竖向/水平向)在短周期(< 0.1 s)甚至超过1.5,统计分析还表明谱比值在短周期段(< 0.1 s)随断层距的增大而减小.     

Simulation of accelerograms from simplified deterministic approach for the 23rd October 2004 Niigata-ken Chuetsu,Japan earthquake

A simple hybrid approach for the simulation of strong ground motion is presented in this paper. This approach is based on the deterministic modelling of rupture plane initially started by Midorikawa, Tectonophysics 218:287–295, (1993) and further modified by Joshi, Pure Appl Geophys (PAGEOPH) 8:161, (2004). In this technique, the finite rupture plane of the target event is divided into several subfaults, which satisfy scaling relationship. In this paper, simulation of strong ground motion due to a rupture buried in a earth medium consisting of several layers of different velocities and thicknesses is made by considering (1) transmission of energy at each layer; (2) frequency filtering properties of medium and earthquake source; (3) correction factor for slip of large and small magnitude earthquakes and (4) site amplification ratio at various stations. To test the efficacy of the developed technique, strong motion records were simulated at different stations that have recorded the 2004 Niigata-ken Chuetsu, Japan earthquake (M _s 7.0). Comparison is made between the simulated and observed velocity and acceleration records and their response spectra. Distribution of peak ground acceleration, velocity and displacement surrounding the rupture plane is prepared from simulated and observed records and are compared with each other. The comparison of synthetic with the observed records over wide range of frequencies shows that the present technique is effective to predict various strong motion parameters from simple deterministic model which is based on simple regression relations and modelling parameters.     

the faulting characteristics of 2008 wenchuan ms8.0 earthquake and its relation with strong ground motion

The 2008 M_S8 Wenchuan earthquake occurs on a high angle listric thrust fault. It is the first time that the near and far field strong ground motion was observed for such special type thrust earthquake. This paper jointly interprets the distribution of peak acceleration of ground motion data with seismogenic structure and slip propagating process to investigate how high angle listric thrust fault controls the pattern of strong ground motion. We found that the distribution of peak acceleration of strong ground motion during the Wenchuan earthquake has four distinctive features: 1)The peak acceleration of ground motion inside the Longmenshan fault zone is large, that is, nearly twice as strong as that outside the fault zone; 2)This earthquake produces significant vertical ground motion, prevailing against horizontal components in the near field; 3)The far field records show that the peak acceleration is generally higher and attenuates slower versus station-fault distance in the hanging wall. It is doubtful that the attenuation of horizontal components also has the hanging wall effect since no evidence yet proving that the unexpected high value at long distance need be omitted; 4)As to the attenuation in directions parallel to the source fault(Yingxiu-Beichuan Fault), the far field records also exhibit azimuthal heterogeneity that the peak acceleration of horizontal components decreases slower in the north-northeastern direction in which the co-seismic slip propagates than that in the backward way. However, the attenuation of vertical component displays very weak heterogeneity of this kind. Synthetically considered with shallow dislocation, high dip angle, and prevailing vertical deformation during co-seismic process of the Wenchuan earthquake, our near and far field ground motion records reflect the truth that the magnitude of ground motion is principally determined by slip type of earthquake and actual distance between the slipping source patches and stations. As a further interpretation, the uniqueness of high angle listric thrust results in that the ground motion effects of the Wenchuan earthquake are similar to that due to a common thrust earthquake in some components while differ in the others.     

Simulating time-histories and pseudo-spectral accelerations from the 1992 Cairo earthquake at the proposed El-Fayoum New City Site,Egypt

El-Fayoum New City represents one of the new urban settlements that are recently erected all over Egypt. Because seismic recordings are not available, I used the stochastic method to simulate the largest damaging earthquake from the closest seismic source to the proposed area of the city. To verify the method and its computed results in Egypt, a study termed “method verification” was performed. I found that the October 12, 1992, earthquake (M _b = 5.8) that occurred southwest of Cairo in the vicinity of the Dahshour region, at the coordinates 29.77°N, 31.07°E, is a significant earthquake to the city. The parameters of the path from the hypocenter of the event to the city were taken into consideration. To determine the site parameters, a shallow seismic refraction survey was carried out in the studied area. Accordingly, I simulated time-histories and pseudo-spectral accelerations from the October 12, 1992, earthquake at the location of seismic profiles. Finally, it is demonstrated that the site is characterized by high ground motion amplification factors, producing a high ground motion acceleration value.     

2015年尼泊尔Gorkha地震强地面运动记录分析

2015年4月25日在尼泊尔Gorkha地区发生M_W7.8地震,距离发震断层约11 km的KATNP台站完整记录了主震的加速度时程.本文根据KATNP台站记录的加速度数据分析了Gorkha地震的地震动特征.结果表明Gorkha地震在KATNP台站处产生的水平向峰值加速度为0.17 g,竖直向峰值加速度为0.19 g,该数值小于科学家们对如此大规模地震产生的地震动的预期,初步推测这可能是由加德满都山谷产生的非线性响应造成的（Dixit et al.,2015）;地震在KATNP台站处产生了地表永久位移,其中竖向永久位移为131.9 cm,水平向永久位移的绝对值为159.2 cm,方向为南偏西19°（199°）,据此可简单推算出断层走向约为289°（109°）.地震产生了脉冲型地震动,影响因素有盆地效应、地震破裂的向前的方向性效应以及滑冲效应,其中盆地效应的周期约为5 s左右,方向性效应产生的速度脉冲的周期约为8 s左右.加速度反应谱显示在0.5 s和5.0 s左右各有一个峰值,前者是由地震破裂的脉冲式滑移产生的大量高频地震动造成的,后者是由于盆地效应和地震破裂的方向性效应造成的.基于阿里亚斯烈度计算的地震动持时约在36~46 s之间,小于与其规模相当的地震产生的地震动持时,并且不同方向上的地震动持时可能与地震破裂方向有关.阿里亚斯烈度随时间的变化比较简单,也反映了Gorkha地震是一次连续的、能量释放相对简单的地震事件.     

台湾集集地震近场地震动的上盘效应

1999年9月21日（当地时间）台湾集集7.6级地震是一个逆断层型地震．用回归分析法对台湾集集地震的加速度峰值数据进行分析，得出了这次地震的水平与垂直向的加速度峰值衰减关系．从残差分布上看，位于断层上盘和下盘上的加速度峰值与从衰减关系所得到的结果相比存在不同的系统偏差，断层上盘地表的加速度峰值较高，而下盘地表的加速度峰值较低．从这次地震的加速度峰值分布等值线图上也可以看出，加速度峰值的分布相对于断层呈现明显的不对称性，上盘衰减较慢而下盘衰减较快．在近断层强地面运动研究、地震危险性分析、设定地震研究与震害预测等工作中，应考虑可能地震的震源机制特点，以便使所用的衰减模型更能反映不同地震环境地区的地震动分布特征．      

Numerical Green’s function method:Application to quantifying ground motion variations of M7 earthquakes

The concept of "numerical Green’s functions" (NGF or Green’s function database) is developed. The basic idea is: a large seismic fault is divided into subfaults of appropriate size, for which synthetic Green’s functions at the surface (NGF) are calculated and stored. Consequently, ground motions from arbitrary kinematic sources can be simulated, rapidly, for the whole fault or parts of it by superposition. The target fault is a simplified, vertical model of the Newport-Inglewood fault in the Los Angeles basin. This approach and its functionality are illustrated by investigating the variations of ground motions (e.g. peak ground velocity and synthetic seismograms) due to the source complexity. The source complexities are considered with two respects: hypocenter location and slip history. The results show a complex behavior, with dependence of absolute peak ground velocity and their variation on source process directionality, hypocenter location, local structure, and static slip asperity location. We concluded that combining effect due to 3-D structure and finite-source is necessary to quan- tify ground motion characteristics and their variations. Our results will facilitate the earthquake hazard assessment projects.     

Conversion between seismic intensity and peak ground shaking

In this paper, the relationships between seismic intensity and peak ground shaking are studied under no specific condition, separately or simultaneously considering the number of building storey and site category, based on data of mean peak values of horizontal ground motion recorded during strong earthquakes. Then, according to the statistical results, the variation of mean peak value with intensity rating is discussed, and schemes of peak ground velocity, peak ground acceleration or response spectrum of an designed earthquake converted from intensity rating are recommended. Finally, a methodology of converting seismic intensity from response spectrum of design earthquake is also discussed, and the conversion scheme is recommended. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 32–40, 1991. This paper is sponsored by the Chinese Joint Seismological Science Foundation.     

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 CHARACTERISTICS OF GAOYOU-BAOYING M4.9 EARTHQUAKE

An M4.9 earthquake occurred at the junction of Gaoyou and Baoying on July 20, 2012. In this paper, 43 sets of strong motion records of the main shock are analyzed. With these data, we analyzed the characteristics of the peak ground motion value, attenuation relation, duration and acceleration response spectrum. We draw the peak acceleration contour map of the region near the epicenter. The contour line is smooth and the trend of long axis is northwest-southeast. Distribution of peak acceleration of the observed records is basically consistent with the real intensity distribution. Compared with the predicted result based on the seismic attenuation relation proposed by Yu Yanxiang and Wang Suyunon for eastern China and the Fifth-generation ground motion zonation map, the horizontal PGA and PGV of Gaoyou-Baoying earthquake are higher than the predicted results that are based on the model of Fifth-generation ground motion zonation map, while the PGV is similar with the predicted results which are based on Yu Yanxiang and Wang Suyun's model. We regressively analyzed the spatial-temporal change curves of the two types of relative ground motion durations. Compared with the predicted results proposed by Bommer et al. (2009) based on the NGA strong motion records, the durations of all the three components of this earthquake are higher. 10 typical recordings' acceleration response spectra with 5% damping are calculated, their peak periods are around 0.1~0.3s. The acceleration response spectrum of the station 32BYT, which has the largest amplitude, is considerably larger than the Chinese code design spectra, while it becomes notably smaller when the period is larger than 0.4s. Compared with the horizontal bedrock acceleration response spectrum predicted by the attenuation relationship for the eastern part of China, the observed response spectrum shape is similar with the predicted ones, while almost all the observed response spectrum values (except station 32YCT)are smaller than the predicted bedrock acceleration response spectrum. These phenomena suggest that this earthquake has a weak impact on the seismic fortification standards in this area. Using H/V single-station spectral ratio method, amplitude and site amplification effect of the two typical stations are calculated, and the results show the H/V values are obviously larger than that of ground microtremor. This suggests that the site of the station has obvious amplification effect on ground motion.     

Estimation of Strong Ground Motions for 2001 Bhuj (M w 7.6), India Earthquake

The strong ground motions for the 2001 Bhuj (M _w 7.6) India earthquake have been estimated on hard rock and B/C boundary (NEHRP) levels using a recently modified version of stochastic finite fault modeling based on dynamic corner frequency (Motazedian and Atkinson in Bull Seismol Soc Am 95, 995–1010 2005). Incorporation of dynamic corner frequency removes the limitations of earlier stochastic methods. Simulations were carried out at 13 sites in Gujarat where structural response recorder (SRR) recordings are available. In addition, accelerograms were simulated at the B/C boundary at a large number of points distributed on a grid. The corresponding response spectra have also been estimated. The values of peak ground accelerations and spectral accelerations at three periods (0.4, 0.75 and 1.25 s) are presented in the form of contour maps. The maximum value of peak ground acceleration (PGA) in the center of meizoseismal zone is 550 cm/s². The response spectral acceleration in same zone is 900 cm/s² (T = 0.4 s), 600 cm/s² (T = 0.75 s) and 300 cm/s² (T = 1.25 s). The innermost PGA contour is on the fault plane. A comparison of the PGA values obtained at 13 sites in this study with those obtained in earlier studies on the same sites, but employing different methods, show that the present PGA values are comparable at most of the sites. The rate of decay of PGA values is fast at short distances as compared to that at longer distances. The PGA values obtained here put some constraints on the expected values from a similar earthquake in the region. A synthetic intensity map has been prepared from the estimated values of PGA using an empirical relation. A comparison with the reported intensity map of the earthquake shows the synthetic MMI values, as expected, are lower by 1 unit compared to reported intensity map. The contour map of PGA along with the contour maps of spectral acceleration at various periods permit the assessment of damage potential to various categories of houses and other structures. Such information will be quite important in planning of mitigation and disaster management programs in the region.     

2022年1月8日青海门源6.9级强地震动模拟研究

为满足数字化时代应急处置的新要求,基于随机有限断层法的发展,利用青海门源地区地下三维速度结构的模型和v_S30数据,通过逐步迭代射线追踪法、格林函数位移解析相位谱和有限断层法的联合计算,以青海门源6.9级地震为例,获得具有地表土层放大效应的强地震动模拟数据,进而绘制出研究区内PGA和烈度分布模拟图;并与实际台站记录的PGA和实地调查烈度结果相比,其烈度区划范围基本一致,同时也验证此联合计算方法可用于未来地震灾害的快速评估,并为灾后应急救援提供参考。     

2018年陕西宁强5.3级地震强地面运动特征及局部场地效应分析

2018年9月12日陕西省宁强县发生5.3级地震,中国数字强震动台网的39个专业台站在此次地震中触发。文章中通过处理捕获的117条三分向加速度记录,给出近场台站的地震动参数,绘制震中附近区域峰值加速度等值线图,其长轴呈西南-东北方向展布。采用实际观测数据与几种常用地震动衰减关系对比,发现霍俊荣衰减预测模型能更好地反映此次地震的影响场。将振幅最大的51GYD台的反应谱与我国抗震设计反应谱比较,采用最小二乘法拟合出不同震中距5个台站各周期谱加速度衰减特性,总结出此次地震的反应谱基本特征。运用H/V谱比法对51GYD土层台和62ZM台阵进行局部场地地震反应分析,研究覆盖土层对地震动的放大作用,及局部地形对峰值加速度和峰值速度的影响过程。     

基于《建筑抗震设计规范》强度包线函数参数取值研究

利用《中国地震动参数区划图》采用的地震动参数衰减关系,以及《中国地震动参数区划图》中地震动峰值加速度和地震动加速度反应谱特征周期反推不同设防烈度和设计地震分组对应的震级和震中距,再根据《建筑抗震设计规范》中各设防水准的峰值加速度确定对应的震级和震中距,进而根据地震动强度包线参数与震级和震中距关系计算地震动强度包线参数的取值,为基于强度包线函数生成人工地震动提供参考,并讨论强度包线参数的取值规律：(1)随着设防烈度的提高,加速度时程曲线上升段持续时间t₁和平稳段持续时间t_s减小,下降段衰减指数c增大;(2)随着地震水准和设计地震分组的提高,加速度时程曲线上升段持续时间t₁和平稳段持续时间t_s增加,下降段衰减指数c减小;(3)在生成人工地震动时,除考虑峰值加速度和设计地震分组影响外,还需要考虑设防烈度影响。     

Modelling of rupture planes for peak ground accelerations and its application to the isoseismal map of MMI scale in Indian region

The rupture plane for an earthquake has been modelledby using the semi empirical technique of Midorikawa(1993). This technique estimates ground accelerationby modelling the rupture process during an earthquake.Modifications in this technique have been made for itsapplication to the Indian region. This has been tested forthe Uttarkashi earthquake of 20th Oct, 1991, India, whichwas well recorded at thirteen stations of installedstrong motion array in this region. After testingseveral possible rupture models, a final model has beenselected and peak ground acceleration due to thismodel is simulated at thirteen different stations.Dependency of methodology on model parameters, e.g.dip and mode of rupture propagation have also beenstudied in detail.Using this technique synthetic isoseismal maps wereprepared by converting peak ground acceleration intoMMI scale. Dependency of rupture models on syntheticisoseismals has also been studied in detail. Usingthis method, peak ground acceleration for the Laturearthquake of Sept 30, 1993 has been obtained atvarious places within meisoseismal area. Synthetic andfield intensity was compared at various well-knownsites. Since the region was not covered by anyinstrumental array during Latur earthquake, thesimulated peak ground accelerations are expected toserve basis of design criteria in this region.     

夏垫断裂MW≥7.5地震动的预测

针对夏垫断裂开展了M_W≥7.5地震动预测研究。首先基于全破裂模式设定震源（使其尽可能涵盖夏垫断裂的未知信息）模拟得到夏垫断裂发生M_W≥7.5地震时研究区域内的地面地震动场,进而依据分位数筛选出各场点的地震动空间分布,讨论了包含不确定震源下的加速度峰值和速度峰值的分布特征,结果显示当夏垫断裂发生M_W7.9地震时,通州城区、北京中心城区均会发生强烈的运动。之后对比讨论了仿真震源下M_W7.5地震所引起的地面运动场的空间变化,结果显示对于同等震级而言,两种震源的模拟结果可以相互印证。      

Housner谱烈度的一种快速近似算法

Housner谱烈度及修正谱烈度作为基于加速度记录时程直接得到的地震动强度参数,与建筑结构破坏及地震宏观烈度存在较高的相关性,是可靠的地震仪器烈度物理参数指标。然而,相对于地面加速度峰值、地面速度峰值等地震动峰值参数,三分量加速度记录对应的谱烈度计算过程较为复杂,耗时相对较长,影响了利用谱烈度确定地震仪器烈度的时效性。基于对强震动加速度记录的统计分析,本文提出了谱烈度的快速近似算法,仅计算4个方向上的谱烈度值,采用其中3点作圆即可获得水平面内谱烈度迹线的近似最大值,使计算速度提高了45倍,且保持了谱烈度作为地震仪器烈度物理指标的精度。利用在汶川M_S 8.0地震等386次M_S > 3.0地震中获取的2701组强震动加速度记录,经可靠性检验,结果表明所提出的Housner谱烈度快速近似算法的计算误差在±4.5%以内,可以同时满足地震仪器烈度速报的可靠性和时效性需求。     

Effects of the hanging wall and footwall on peak acceleration during the Jiji (Chi-Chi), Taiwan Province, earthquake

Introduction In the study of ground motion attenuation model the considered parameters are generally sim-plified as ground motion parameters (acceleration, velocity, displacement, response spectrum, dura-tion, etc), earthquake magnitude, distance and site condition. As the accumulation of ground motion records, it was found that the characteristic of ground motion attenuation in tectonicly compressional region was different from that in tectonicly tensional region and the peak ground accelerat…