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.     

2019年10月28日甘肃夏河5.7级地震强地面运动特征分析

2019年10月28日甘肃省夏河县发生5.7级地震,中国数字强震动台网的18个专业台站在此次地震中触发。本文处理捕获的54条三分向加速度记录,给出近场台站的地震动参数,绘制了震中附近区域峰值加速度等值线图,其长轴呈WN-ES方向展布。将实际观测数据与几种常用地震动衰减关系对比,发现俞言祥^[1]短轴衰减预测模型能更好地反映此次地震的影响场。将振幅较大的62LBL、62BLX台的反应谱与我国抗震设计反应谱比较,采用最小二乘法拟合出不同震中距5个台站各周期谱加速度衰减特性,总结了此次地震的反应谱基本特征。运用H/V谱比法对4个典型台站进行场地地震反应分析,研究了局部场地条件对峰值加速度和峰值速度的影响过程。     

对二滩水电站坝区场地地面运动的估计

本文以二滩地区为例,介绍如何结合地震危险性分析的结果,采用理论模拟方法估计工作地区近场强地面运动情况,为工程抗震设计提供必需的地面运动动力学参数。 文中采用理论方法计算剪切位错点源的格林函数;根据断裂动力学模型、近场地面峰值加速度衰减曲线和场地附近的加速度谱来标定震源模型;计算了若在坝区附近发生M_s=6.2级地震时场地的综合地震图、加速度傅氏谱和反应谱等。由理论模拟计算的峰值加速度衰减曲线和观测结果符合较好。在震中距为20km处的加速度谱的截止频率f_max与观测值相吻合（约为8 Hz左右）,相应的峰值加速度为211cm/s²,振动持续时间为3.4秒。     

熊本MW7.0地震近场地表与井下地震动对比研究

选取日本熊本M_W7.0地震断层距小于200 km的82个近场KiK-net台站记录到的三分量记录数据进行基线校正后,获得近场地面运动水平向的峰值加速度PGA、峰值速度PGV及周期为0.2,1,2,3,5和10 s的加速度反应谱数据,并与美国NGA-West2的地震动预测模型相比较,研究熊本地震地表和井下地震动峰值及反应谱的衰减特征,通过比较KiK-net台站地表与井下记录结果,探讨浅层场地放大效应的影响。研究结果表明:① 对于井下观测结果,NGA-West2的地震动模型对PGA和短周期0.2 s的反应谱的预测值与井下观测值相比整体偏高,而PGV和较长周期地震动（如1,2和3 s的反应谱）的预测值与井下观测值较为吻合;② 地表观测记录的PGA,PGV和周期为0.2—3 s的反应谱残差整体上随v_S30对数值的增大呈线性减小的趋势,而周期为5 s和10 s的长周期部分,其场地效应的影响很小;③ 相对于井下记录,地表记录的地震动PGA,PGV和周期为0.2,1和2 s的反应谱有明显的放大,这种放大作用随浅层场地剪切波速的增大而减小;周期为3,5和10 s时长周期地震动的放大效应很小。      

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

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

A note on tsunami caused by submarine slides and slumps spreading in one dimension with nonuniform displacement amplitudes

Tsunami created by spreading submarine slides and slumps with spatially variable final uplift are investigated in the near-field using a kinematic model. It is shown that for velocities of spreading comparable to and smaller than the long period tsunami velocity (g is the acceleration due to gravity and h is the ocean depth), the models with spatially uniform final uplift of the accumulation and depletion zones provide good approximation for the tsunami amplitudes in the near-field. For spreading velocities 2–5 times greater than c_T, and for applications that use wavelengths of the order of the source dimensions, the spatial variability of the final uplift has to be considered in estimation of the high-frequency tsunami amplitudes in the near-field.     

A note on tsunami amplitudes above submarine slides and slumps

Tsunami generated by submarine slumps and slides are investigated in the near-field, using simple source models, which consider the effects of source finiteness and directivity. Five simple two-dimensional kinematic models of submarine slumps and slides are described mathematically as combinations of spreading constant or slopping uplift functions. Tsunami waveforms for these models are computed using linearized shallow water theory for constant water depth and transform method of solution (Laplace in time and Fourier in space). Results for tsunami waveforms and tsunami peak amplitudes are presented for selected model parameters, for a time window of the order of the source duration.The results show that, at the time when the source process is completed, for slides that spread rapidly (c_R/c_T≥20, where c_R is the velocity of predominant spreading), the displacement of the free water surface above the source resembles the displacement of the ocean floor. As the velocity of spreading approaches the long wavelength tsunami velocity the tsunami waveform has progressively larger amplitude, and higher frequency content, in the direction of slide spreading. These large amplitudes are caused by wave focusing. For velocities of spreading smaller than the tsunami long wavelength velocity, the tsunami amplitudes in the direction of source propagation become small, but the high frequency (short) waves continue to be present. The large amplification for c_R/c_T1 is a near-field phenomenon, and at distances greater than several times the source dimension, the large amplitude and short wavelength pulse becomes dispersed.A comparison of peak tsunami amplitudes for five models plotted versus L/h (where L is characteristic length of the slide and h is the water depth) shows that for similar slide dimensions the peak tsunami amplitude is essentially model independent.     

Attenuation relations for strong seismic ground motion in the Himalayan region

Strong motion data from various regions of India have been used to study attenuation characteristics of horizontal peak acceleration and velocity. The strong ground motion data base considered in the present work consists of various earthquakes recorded in the northern part of India since 1986 with magnitudes 5.7 to 7.2. Using these data, relations for horizontal peak acceleration and velocity, which are $$\begin{gathered} log_{10} a = 1.14 + 0.31M + 0.65log_{10} R \hfill \\ log_{10} v = 0.571 + 0.41M + 0.768log_{10} R \hfill \\ \end{gathered} $$ have been proposed wherea is the peak horizontal acceleration in cm/sec²,v is the peak horizontal velocity in mm/sec,M is body wave magnitude, andR is the hypocentral distance in km. The proposed relations are in reasonable agreement with the small amount of strong ground motion data available for the northern part of India. The present results will be useful in estimating strong ground motion parameters and in the earthquake resistant design in the Himalayan region.     

A SIMPLIFIED METHOD FOR PREDICTION OF KINEMATIC SOIL–FOUNDATION INTERACTION EFFECTS ON PEAK HORIZONTAL ACCELERATION OF A RIGID FOUNDATION

The kinematic soil–foundation interaction changes the free field ground motion to a different motion at the foundation of a structure. This interaction effect may be expressed by the ratio of the peak horizontal acceleration of a rigid and relatively lightweight foundation to the peak horizontal acceleration at the ground surface in the free field. It is found that the interaction effect can be defined by a simple function of the ratio of the peak horizontal ground velocity and ground acceleration in the free field, the length of the foundation and the shear wave velocity in the soil. Predictive equations for the kinematic soil foundation effect are derived using 350 strong motion records generated by 114 earthquakes world-wide. At the same time, an attenuation relationship is derived for the ratio of the peak horizontal ground velocity and acceleration from the same set of data. Ten case histories are studied; the interaction effects are calculated by using the predictive equations and then compared with measured field values. The results of the comparison illustrate the degree of predictive capability of the method when the foundation mass and the inertial soil–foundation interaction are not considered.     

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.     

Generation of probabilistic displacement response spectra for displacement-based design

Displacement response spectrum (DRS), as the input, is of great significance to the displacement-based design just like the acceleration response spectrum to the traditional force-based design. Although the procedure of performance-based, in particular the displacement-based design has achieved considerable development, there is not a general DRS covering an enough long period range for common seismic design yet. This paper develops a systematic ground motion data processing procedure for the purpose of correcting the noise in the earthquake records and generating consistent DRS for seismic design. An adaptive algorithm is proposed to determine the cutoff frequency of the high-pass digital filter. The DRS of more than 500 recorded earthquake ground motions are generated and they are classified into three groups according to the ratio of the peak ground acceleration to the peak ground velocity (A/V) and/or the ratio of the peak ground velocity to the peak ground displacement (V/D). In each group, all the ground motions are normalized with respect to a selected scaling factor. Their corresponding DRS are obtained and then averaged to get the mean and standard deviation DRS, which can be used for both deterministic and probabilistic displacement-based design.     

A note on scaling peak acceleration, velocity and displacement of strong earthquake shaking by Modified Mercalli Intensity (MMI) and site soil and geologic conditions

It is shown that the overall trends of peak amplitudes of strong motion acceleration, velocity and displacement versus site intensity in the Western United States are consistent with our previous inferences, ¹⁸ but the larger present data base allows us to begin with identification of the additional changes of average peak amplitudes also in terms of the local soil and local geologic site conditions. These trends display a reversal of peak amplitude trends in going from small to large intensities, making it difficult to reconcile the simple scaling of design spectra in terms of peak ground acceleration.     

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

1 Introduction Earthquake disaster investigations show that numerous strong earthquakes were caused by remobilization of active faults. Major casualties and severe damage to buildings as well as signi?cant economic losses resulted from the ground motions of strongearthquakescausedbyactivefaultslocatedbeneath urban areas. Recently, the potential hazard prediction of and its mitigation against active faults located beneath urban areas have become an important research topic for seismologists and…     

Attenuation relations for horizontal peak ground acceleration and response spectrum in northeastern Tibetan Plateau region

Introduction Developing local attenuation relations of ground motion is one of the key steps in seismic hazard assessment. Because of inadequate strong ground motion records in China, the attenuation relations used in China are usually developed by using the transforming method (Hu, Zhang, 1983; HU, ZHANG, 1984). To use this method, we need to have both the attenuation relation of seismic intensity for the studied region and the attenuation relations of seismic intensity and ground mo-tion…     

Duration of earthquake fault motion in California

Duration of high frequency (5–25 Hz) radiation of energy from earthquake sources in California is consistent with the estimates of fault length and with dislocation velocity estimates of 2–3 km/sec. This duration can be described by an exponential function of magnitude for 2·5 < M < 7·5. It is related to the times it takes the dislocation to spread over the fault width (1/f₂), and the fault length (～ 1/f₁), and to reach its ultimate amplitude (T₀). The results in this paper can be used to estimate the range of amplitudes and the duration of long period pulses of strong ground motion near faults, as these long period pulses can be related to the properties of high-frequency radiation from the source. Such pulses must be considered in the analyses of yielding structures, when the average peak acceleration of the pulse exceeds the yield resistance seismic coefficient of the structure.     

Seismic motion and response prediction alternatives

Statistical methods are available which predict the maximum response of simple oscillators given the peak acceleration (A_p), peak velocity (V_p) or peak displacement (D_p) of seismic ground motions. An alternative parameter, namely an ordinate (or ordinates) of the Fourier amplitude spectrum of ground motion acceleration, FS(f), may in fact be a preferred predictor of peak response, especially in a frequency range close to f. Other statistical methods (attenuation laws) use distance R and other parameters such as magnitude (M), Modified Mercalli epicentral Intensity (I_o) and Modified Mercalli site Intensity (MMI or I_s) to predict spectral velocity (S_v(f)), etc. In using such approaches, it is most desirable to know the total uncertainty in the predicted peak response of the system given the starting parameter values. An extensive strong motion data set is used to study these questions, The most direct prediction models are found to be preferable (have lower prediction dispersion) but data may not be available in all regions to permit their use.     

An evaluation of strong motion records and a new parameter A95

From the past studies of acceleration records, it is generally recognised that peak acceleration is not a suitable measure of the ground motion for design purposes. Usually, peak accelerations are difficult to predict. They occur in high frequency pulses containing very little energy and therefore are of little consequence. On the other hand acceleration has direct and useful application in engineering design. A new parameter A₉₅ is introduced in this paper which is defined as that level of acceleration which contains up to 95 per cent of the Arias Intensity. From the study of 135 records, it is noticed that the variation of the fraction of Arias Intensity with acceleration level can be described by a very well defined relationship. It is also observed that A₉₅ and Arias Intensity bear a stable relationship. Furthermore, the Arias Intensity at a site correlates more closely with the earthquake magnitude and the distance of the site from the source in attenuation relationships. From the available data, relationships are derived which show that there is a difference between the near-field and the far-field attenuation laws. The following relationships are obtained: where E_s = Arias Intensity, E_x = fraction of Arias Intensity above an acceleration level A, M = surface wave magnitude, R = nearest distance to the energy source or focal distance in km.     

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.     

Artificial ground motion compatible with specified peak velocity and target spectrum

In this paper, a method, which synthesizes the artificial ground motion compatible with the specified peak velocity as well as the target acceleration response spectrum, was proposed. In this method, firstly, an initial acceleration time history α8^（0） （t）, which satisfies the prescribed peak ground acceleration, the target spectral acceleration ST（ω, ζ）,and the specified intensity envelope, is generated by the traditional method that generates the requency domain; secondly,α8^（0） （t）is further modulated by superimposing narrow-band time histories upon it in the time domain to make its peak velocity, approach the target peak ground velocity, and at the same time to improve its fitting precision to the target spectrum. Numerical examples show that this algorithm boasts high calculation precisions.     

Scaling of peak ground motions from digital recordings of small earthquakes at Campi Flegrei,southern Italy

The dependence of peak ground acceleration and velocity on seismic moment is studied for a set of small earthquakes (0.7<M _L<3.2) recorded digitally at distances of a few km in the Campi Flegrei volcanic area near Naples, Italy, during the ground uplift episode of 1982–1984. Numerical simulations, using the -square spectral model with constant stress drop and ane ^–kf high frequency decay, fit well both the velocity and acceleration data for an averagek=0.015. The observed ground motions in the 1–24 Hz frequency band appear to consist of radiation from simple sources modified only slightly by attenuation effects. Moreover, the scaling of peak values agrees closely with those determined in nonvolcanic areas, once the difference in stress drop is taken into account.