Spectral parameters of ground motion in different regions: comparison of empirical models

Fourier-amplitude spectrum is one of the most important parameters describing earthquake ground motion, and it is widely used for strong ground motion prediction and seismic hazard estimation. The relationships between Fourier-acceleration spectra, earthquake magnitude and distance were analysed for different seismic regions (the Caucasus and Taiwan island) on the basis of ground motion recordings of small to moderate (3.5≤M_L≤6.5) earthquakes. It has been found that the acceleration spectra of the most significant part of the records, starting from S-wave arrival, can be modelled accurately by the Brune's “ω-squared” point-source model. Parameters of the model are found to be region-dependent. Peak ground accelerations and response spectra for condition of rock sites were calculated using stochastic simulation technique and obtained models of source spectra. The modelled ground-motion parameters are compared with those predicted by recent empirical attenuation relationship for California.     

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.     

Canyon topography effects on ground motion at Feitsui damsite

The free-field accelerograms along Feitsui Canyon are analyzed and modeled by a numerical scheme to study the effect of canyon topography. Since six strong-motion accelerometers (SC1–SC6) were deployed along the Feitsui Canyon in 1991; there are 14 earthquakes (4.9≤M_L≤6.6) recorded by these stations until June 1996, but only five triggered all six stations. The maximum PGA value is 68.6 cm s⁻² recorded at station SC1. According to the present data, the effect of the dam on the ground motions at canyon stations can be negligible. The amplitude of ground motion on the slopes of the canyon is bigger than that at its trough. The integral equation method is applied to a two dimensional model of Feitsui Canyon to study the effects of the canyon topography. We choose the ground motion of SC3 or SC4 station at the trough of the canyon as the input motion for the model, which is then used to predict the ground motion at the other five stations. Apart from the earthquake close to the damsite, the simple model can reproduce the observed accelerations at all frequencies below 4 Hz. Overall, the numerical method can well predict the ground motion along the canyon, although the high-frequency simulation is underestimated.     

Stochastic Finite Fault Modeling of Subduction Zone Earthquakes in Northeastern India

In this article, a stochastic finite fault source model is calibrated to estimate ground motion in northeastern India for intermediate depth events originating in the Indo-Burmese tectonic domain. A total of 47 three-component accelerograms from eight events with magnitudes ranging from M _w 4.8–6.4 are used to estimate the input source and site parameters of the finite fault source model. Key seismic parameters such as stress drop (Δσ) and site amplification function are determined from the recorded strong motion data. The obtained stress drop of the eight recorded events lies in between 105 and 165 bars.     

Ground-Motion Scaling in the Western Alps

In order to empirically obtain the scaling relationships for the high-frequency ground motion in the Western Alps (NW Italy), regressions are carried out on more than 7500 seismograms from 957 regional earthquakes. The waveforms were selected from the database of 6 three-component stations of the RSNI (Regional Seismic network of Northwestern Italy). The events, M _W ranging between 1.2 and 4.8, were recorded within a hypocentral distance of 200 km during the time period: 1996–2001. The peak ground velocities are measured in selected narrow-frequency bands, between 0.5 and 14 Hz. Results are presented in terms of a regional attenuation function for the vertical ground motion, a set of vertical excitation terms at the reference station STV2 (hard-rock), and a set of site terms (vertical and horizontal), all relative to the vertical component of station STV2.The regional propagation of the ground motion is modeled after quantifying the expected duration of the seismic motion as a function of frequency and hypocentral distance. A simple functional form is used to take into account both the geometrical and the anelastic attenuation: a multi-variable grid search yielded a quality factor Q(f) = 310f ^0.20, together with a quadri-linear geometrical spreading at low frequency. A simpler, bi-linear geometrical spreading seems to be more appropriate at higher frequencies (f > 1.0 Hz). Excitation terms are matched by using a Brune spectral model with variable, magnitude-dependent stress drop: at M _w 4.8, we used Δσ = 50 MPa. A regional distance-independent attenuation parameter is obtained (κ₀ = 0.012 s) by modelling the average spectral decay at high frequency of small earthquakes.In order to predict the absolute levels of ground shaking in the region, the excitation/attenuation model is used through the Random Vibration Theory (RVT) with a stochastic point-source model. The expected peak-ground accelerations (PGA) are compared with the ones derived by Ambraseys et al. (1996) for the Mediterranean region and by Sabetta and Pugliese (1996) for the Italian territory.     

Ground Motion Relations for Active Regions in India

In this article, a study on development of ground motion prediction equations (GMPEs) is undertaken for seismically active regions in India. To derive the equations, the seismically active regions are divided into four units based on seismotectonic setting and geology. Due to lack of strong motion data, a stochastic finite-fault simulation method is used for generating a complete synthetic database with respect to magnitude and distance. The input parameters in the stochastic seismological model, such as site amplification and stress drop, are first derived from the past strong-motion data. A total of 236 three-component records from 62 earthquakes with magnitudes ranging from M _w 3.4 to 7.8 are used to calibrate the seismological model. The obtained stress drops of these 62 events lie in between 60 and 165 bars. With the help of a large synthetic database generated from the calibrated seismological model, ground motion relations for 5 % damped spectral acceleration are obtained by regression analysis. The developed ground motion relations are compared with the existing GMPEs of the other active regions in the world. Although the proposed equations have trends similar to those of the existing relations, there are some differences attributed to stress drop and the quality factor of active regions in India. These relations will be useful to prepare spectral acceleration hazard maps of India for a given annual probability of exceedance.     

Evaluation of models for Fourier amplitude spectra for the Taiwan region

In the Taiwan region, the empirical spectral models for estimating ground-motion parameters were obtained recently on the basis of recordings of small to moderate (5.0≤M_L≤6.5) earthquakes. A large collection of acceleration records from the M_L=7.3 Chi-Chi earthquake (21 September, 1999) makes it possible to test the applicability of the established relationships in the case of larger events. The comparison of ground-motion parameters (Fourier amplitude spectra, peak accelerations and response spectra), which were calculated using the models, and the observed data demonstrates that the models could provide an accurate prediction for the case of the Chi-Chi earthquake and the largest aftershocks. However, there are some peculiarities in the ground-motion frequency content and attenuation that, most probably, are caused by the features of the rupture process of the large shallow earthquake source.     

Empirical model for estimating Fourier amplitude spectra of ground acceleration in Taiwan region

A collection of ground‐motion recordings (1070 acceleration records) of moderate (5.1⩽M_L⩽6.5) earthquakes obtained during the execution of the Taiwan Strong Motion Instrumentation Program (TSMIP) since 1991 was used to study source scaling model and attenuation relations for a wide range of earthquake magnitudes and distances and to verify the models developed recently for the Taiwan region. The results of the analysis reveal that the acceleration spectra of the most significant part of the records, starting from S‐wave arrival, can be modelled accurately using the Brune's ω‐squared source model with magnitude‐dependent stress parameter Δσ, that should be determined using the recently proposed regional relationships between magnitude (M_L) and seismic moment (M₀) and between M₀ and Δσ. The anelastic attenuation Q of spectral amplitudes with distance may be described as Q=225 ƒ^1.1 both for deep (depth more than 35 km) and shallow earthquakes. The source scaling and attenuation models allow a satisfactory prediction of the peak ground acceleration for magnitudes 5.1⩽M⩽6.5 and distances up to about 200 km in the Taiwan region, and may be useful for seismic hazard assessment. Copyright © 2000 John Wiley & Sons, Ltd.     

Estimating ground motions using recorded accelerograms

A procedure for estimating ground motions using recorded accelerograms is described. The premise of the study is the assumption that future ground motions will be similar to those observed for similar site and tectonic situations in the past. Direct techniques for scaling existing accelerograms have been developed, based on relative estimates of local magnitude,M _L . Design events are described deterministically in terms of fault dimension, tectonic setting (stress drop), fault distance, and site conditions. A combination of empirical and theoretical arguments is used to develop relationships betweenM _L and other earthquake magnitude scales. In order to minimize scaling errors due to lack of understanding of the physics of strong ground motion, the procedure employs as few intermediate scaling laws as possible. The procedure conserves a meaningful measure of the uncertainty inherent when predicting ground motions from simple parameterizations of earthquake sources and site conditions.     

Comparison between different earthquake magnitudes determined by China Seismograph Network

By linear regression and orthogonal regression methods, comparisons are made between different magnitudes (lo-cal magnitude ML, surface wave magnitudes MS and MS7, long-period body wave magnitude mB and short-period body wave magnitude mb) determined by Institute of Geophysics, China Earthquake Administration, on the basis of observation data collected by China Seismograph Network between 1983 and 2004. Empirical relations between different magnitudes have been obtained. The result shows that: 1 As different magnitude scales reflect radiated energy by seismic waves within different periods, earthquake magnitudes can be described more objectively by using different scales for earthquakes of different magnitudes. When the epicentral distance is less than 1 000 km, local magnitude ML can be a preferable scale; In case M<4.5, there is little difference between the magnitude scales; In case 4.5MS, i.e., MS underestimates magnitudes of such events, therefore, mB can be a better choice; In case M>6.0, MS>mB>mb, both mB and mb underestimate the magnitudes, so MS is a preferable scale for deter-mining magnitudes of such events (6.08.5, a saturation phenomenon appears in MS, which cannot give an accurate reflection of the magnitudes of such large events; 2 In China, when the epicentral distance is less than 1 000 km, there is almost no difference between ML and MS, and thus there is no need to convert be-tween the two magnitudes in practice; 3 Although MS and MS7 are both surface wave magnitudes, MS is in general greater than MS7 by 0.2~0.3 magnitude, because different instruments and calculation formulae are used; 4 mB is almost equal to mb for earthquakes around mB4.0, but mB is larger than mb for those of mB≥4.5, because the periods of seismic waves used for measuring mB and mb are different though the calculation formulae are the same.     

Interrelation between seismicity parameters and delimiting potential seismic sources in a seismic statistical region and its influence on seismic risk estimation

In the paper, we have discovered the abnormal area distribution features of maximum variation values of ground motion parameter uncertainty with different probabilities of exceedance in 50 years within the range of 100°～120°E,29°～42°N for the purpose to solve the problem that abnormal areas of maximum variation values of ground motion parameter uncertainties emerge in a certain cities and towns caused by seismicity parameter uncertainty in a seismic statistical region in an inhomogeneous distribution model that considers tempo-spatial nonuniformity of seismic activity. And we have also approached the interrelation between the risk estimation uncertainty of a site caused by seismicity parameter uncertainty in a seismic statistical region and the delimitation of potential sources, as well as the reasons for forming abnormal areas. The results from the research indicate that the seismicity parameter uncertainty has unequal influence on the uncertainty of risk estimation at each site in a statistical region in the inhomogeneous distribution model, which relates to the scheme for delimiting potential sources. Abnormal areas of maximum variation values of ground motion parameter uncertainty often emerge in the potential sources of Mu≥8 (Mu is upper limit of a potential source) and their vicinity. However, this kind of influence is equal in the homogeneous distribution model. The uncertainty of risk estimation of each site depends on its seat. Generally speaking, the sites located in the middle part of a statistical region are only related to the seismicity parameter uncertainty of the region, while the sites situated in or near the juncture of two or three statistical regions might be subject to the synthetic influences of seismicity parameter uncertainties of several statistical regions.     

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.     

On the use of JMA intensity in earthquake early warning systems

The estimation of strength of shaking at a site from the initial P-wave portion of ground motion is the key problems for shortening the alert time of the earthquake Early Warning (EEW). The most of the techniques proposed for the purpose utilize (a) ground motion models based on the estimated magnitude and hypocentral distance, or (b) the interim proxies, such as initial vertical displacement P _d . We suggest the instrumental Japan Meteorological Agency (JMA) intensity (JMA_I) as a characteristic for fast estimation of damage potential in the EEW systems. We investigated the scaling relations between JMA_I measured using the whole earthquake recordings (overall intensity) and using particular time intervals of various duration (2.0–8.0 s) starting from the P-wave arrival (preliminary intensity). The dataset included 3,660 records (K-NET and the KiK-net networks) from 55 events (M _W 4.1–7.4) occurred in 1999–2008 in Japan. We showed that the time interval of 4–5 s from the P-wave arrival can be used for reliable estimations of the overall intensity with the average standard error of about 0.5 JMA units. The uncertainty in the prediction may be reduced by consideration of local site conditions or by development of the station-specific models.     

Characteristics of strong ground motion during the 1999 Chi-Chi earthquake (Taiwan) and large aftershocks: comparison with the previously established models

In the Taiwan region, the empirical spectral models for estimating ground motion parameters were obtained recently on the basis of recordings of small to moderate (5.0≤M_L≤6.5) earthquakes. A large collection of acceleration records from the recent M_L 7.3 (M_W 7.6) Chi-Chi earthquake (20 September, 1999) makes it possible to test the applicability of the established relationships in the case of larger events. The comparison of ground motion parameters (peak accelerations and response spectra), which were calculated using the stochastic approach based on the modeled Fourier amplitude spectra, and the observed data demonstrates that the models may be successfully used for ground motion prediction for earthquakes of magnitudes up to M_L=6.8–7.0 and hypocentral depth more than 10 km. To satisfy to the peculiarities of ground motion during shallow (depth less than 10 km) and larger (M_L>7.0) events, the models were revised.     

ML of the 1989 Loma Prieta earthquake

The use of regional attenuation in computing the local magnitude, M_L, from strong motion data gathered at distances less than 100 km may lead to systematic underestimates approaching 0·5 magnitude units (Trifunac & Herak, Soil Dynamics and Earthquake Engineering, 1992, 18, 229-41). The use of the attenuation law Att(Δ), for example, with synthetic estimates of Wood-Anderson seismometer response, during the Loma Preita earthquake, leads to estimates of M_L which agree with the surface wave and moment magnitudes, and which are essentially distance-independent.     

Attenuation of modified mercalli intensity in New Zealand earthquakes

This paper describes a comprehensive study of the attenuation of Modified Mercalli intensity in New Zealand earthquakes, which has resulted in significantly different and new findings when compared with those of earlier work. The current study used recently revised magnitudes for the 30 events in the carefully selected data set. Magnitudes ranged from M_L = 5.0 to M_s = 7.8. Special effort was also put into establishing the depths of the events, which ranged from very shallow to 65 km. An expression of a form also used in peak ground acceleration attenuation studies was adopted, I = a + bM + cr + dlog₁₀r with r being taken as the mean distance from the centre of the fault rupture surface to each isoseismal. A two-step stratified regression analysis was used because it results in more realistic estimates of standard errors, although its mean curves were virtually the same as those derived from the standard one-step method. It was found possible to model accurately the attenuation of the deepest events from regressions of the data of the shallowest events, but the inverse was not true. In a study of the influence of source mechanisms it was found that the attenuation was the same for events with normal and strike-slip faulting. In contrast, events with reverse fault mechanisms were found to have higher intensities (for the same M and r) than normal and strike-slip events, by a factor corresponding to that found by Campbell in a study of peak ground accelerations, although the statistical test was narrowly short of the conventional significance level.     

Selection and ranking of ground motion models for seismic hazard analysis in the Pyrenees

The issue addressed in this paper is the objective selection of appropriate ground motion models for seismic hazard assessment in the Pyrenees. The method of Scherbaum et al. (2004a) is applied in order to rank eight published ground motion models relevant to intraplate or to low deformation rate contexts. This method is based on a transparent and data-driven process which quantifies the model fit and also measures how well the underlying model assumptions are met. The method is applied to 15 accelerometric records obtained in the Pyrenees for events of local magnitude between 4.8 and 5.1, corresponding to moment magnitudes ranging from 3.7 to 3.9. Only stations at rock sites are considered. A total of 720 spectral amplitudes are used to rank the selected ground motion models. Some control parameters of these models, such as magnitude and distance definitions, may vary from one model to the other. It is thus important to correct the selected models for their difference with respect to the magnitude and distance definitions used for the Pyrenean data. Our analysis shows that, with these corrections, some of the ground motion models successfully fit the data. These are the Lussou et al. (2001) and the Berge-Thierry et al. (2003) models. According to the selected ground motion models, a possible scenario of a magnitude 6 event is proposed; it predicts response spectra accelerations of 0.08–0.1 g at 1 Hz at a hypocentral distance of 10 km.     

Characteristics of strong ground motion for typical Australian intra-plate earthquakes and their relationship with the recommended response spectra

Global epicentre maps show that the majority of earthquakes are inter-plate, although moderate to large earthquakes do occur intra-plate, i.e. within the plates. The seismicity of the Australian continent is typical of intra-plate environments and a magnitude M_L 6 earthquake has an average return period of about 5 years. Recordings of Australian intra-plate earthquakes are investigated here to characterise their frequency content, peak acceleration and duration.Due to lack of quality strong motion records of large intra-plate earthquakes at short distances, synthetic seismograms are commonly used for testing structural behaviour. An empirical Green's Function method (Geophys. Res. Lett., 5 (1978), 1–4; Proceedings of the Third International Microzonation Conference, Seattle, USA, vol. 1, (1982), pp. 447–458.) is chosen to simulate a large earthquake by summation in time of a number of smaller earthquakes or sub-events, each given a slightly different origin time to represent more realistically the propagation of a rupture along an assumed fault plane. In the first instance, recordings on rock of the magnitude M_L 2.3 aftershock of the 29 December 1989 Newcastle earthquake were used as sub-events to simulate the main shock of magnitude M_L 5.6. Validation studies for events recorded elsewhere in Australia are also considered.The response spectra of such synthetic events will be compared with the recommended spectra developed empirically from a statistical analysis of strong motion data for magnitude 5.4–6.5 intra-plate earthquakes recorded in other parts of the world and normalised to a peak ground velocity of 50 mm/s which is typical for a return period of 500 years in Australia (Australasian Structural Engineering Conference, Auckland, New Zealand, (1998), pp. 439–444.). Preliminary results from this comparison with the response spectra recommended for the Building Code of Australia show that the synthetic waveforms produced by this method are realistic and can be used to represent ground motion during typical Australian intra-plate earthquakes.     

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.

     

A spatial correlation model of peak ground acceleration and response spectra based on data of the Istanbul Earthquake Rapid Response and Early Warning System

Ground motion intensity measures such as the peak ground acceleration (PGA) and the pseudo-spectral acceleration (PSA) at two sites due to the same seismic event are correlated. The spatial correlation needs to be considered when modeling ground-motion fields for seismic loss assessments, since it can have a significant influence on the statistical moments and probability distribution of aggregated seismic loss of a building portfolio.Empirical models of spatial correlation of ground motion intensity measures exist only for a few seismic regions in the world such as Japan, Taiwan and California, since for this purpose a dense observation network of earthquake ground motion is required. The Istanbul Earthquake Rapid Response and Early Warning System (IERREWS) provides one such dense array with station spacing of typically 2 km in the urban area of Istanbul. Based on the records of eight small to moderate (M_w3.5–M_w5.1) events, which occurred since 2003 in the Marmara region, we establish a model of intra-event spatial correlation for PGA and PSA up to the natural period of 1.0 s.The results indicate that the correlation coefficients of PGA and short-period PSA decay rapidly with increasing interstation distance, resulting in correlation lengths of approximately 3–4 km, while correlation lengths at longer natural periods (above 0.5 s) exceed 6 km. Finally, we implement the correlation model in a Monte Carlo simulation to evaluate economic loss in Istanbul's district Zeytinburnu due to a M_w7.2 scenario earthquake.