Peak ground motion predictions in India: an appraisal for rock sites

Proper selection and ranking of Ground Motion Prediction Equations (GMPEs) is critical for successful logic-tree implementation in probabilistic seismic hazard analysis. The present study explores this issue in predicting peak ground accelerations at the rock sites in India. Macroseismic intensity data complemented with limited strong ground-motion recordings are used for the purpose. The findings corroborate the possible conformity between the GMPEs developed for tectonically active shallow crust across the globe. On the other hand, the relevant GMPEs in the intraplate regions cluster into two different groups with the equations of lower ranks catering to higher ground motions. The earthquakes in the subduction zones have significant regional implications. However, affinity in the ground-motion attenuations between the major interface events (M _W > 7.4) in Andaman-Nicobar, Japan and Cascadia, respectively, is noted. This can be also observed for the intraslab events in the Hindukush and Taiwan respectively. Overall, we do not observe any significant advantage with the equations developed using the regional data. These findings are expected to be useful in probabilistic seismic hazard analysis across the study region.     

Consistency of ground-motion predictions from the past four decades

Due to the limited observational datasets available for the derivation of ground-motion prediction equations (GMPEs) there is always epistemic uncertainty in the estimated median ground motion. Because of the increasing quality and quantity of strong-motion datasets it would be expected that the epistemic uncertainty in ground-motion prediction (related to lack of knowledge and data) is decreasing. In this study the predicted median ground motions from over 200 GMPEs for various scenarios are plotted against date of publication to examine whether the scatter in the predictions (a measure of epistemic uncertainty) is decreasing with time. It is found that there are still considerable differences in predicted ground motions from the various GMPEs and that the variation between estimates is not reducing although the ground motion estimated by averaging median predictions is roughly constant. For western North America predictions for moderate earthquakes have show a high level of consistency since the 1980s as do, but to a lesser extent, predictions for moderate earthquakes in Europe, the Mediterranean and the Middle East. A good match is observed between the predictions from GMPEs and the median ground motions based on observations from similar scenarios. Variations in median ground motion predictions for stable continental regions and subduction zones from different GMPEs are large, even for moderate earthquakes. The large scatter in predictions of the median ground motion shows that epistemic uncertainty in ground-motion prediction is still large and that it is vital that this is accounted for in seismic hazard assessments.     

基于数值模拟的设定地震场地谱时程拟合方法

本文将确定性数值模拟方法与地震动预测方程相结合,提出了一种重大水电工程场址设定地震的地震动时程生成方法。该方法基于场址设定地震,首先采用地震动预测方程确定场址的场地相关反应谱;其次建立包含震源和场址的场地模型,通过确定性数值模拟方法生成场址地震动时程;最后对生成的场址地震动时程进行调整,使其反应谱与设计谱相一致,用于工程抗震分析。这一方法生成的地震动时程既考虑了震源机制、传播路径以及局部场地效应等物理背景,又与场地相关的设计地震反应谱保持一致,为重大工程抗震分析与评价提供了一种新的思路。     

Simulations of Seismic Hazard for the Pacific Northwest of the United States from Earthquakes Associated with the Cascadia Subduction Zone

-- We investigate the impact of different rupture and attenuation models for the Cascadia subduction zone by simulating seismic hazard models for the Pacific Northwest of the U.S. at 2% probability of exceedance in 50 years. We calculate the sensitivity of hazard (probabilistic ground motions) to the source parameters and the attenuation relations for both intraslab and interface earthquakes and present these in the framework of the standard USGS hazard model that includes crustal earthquakes. Our results indicate that allowing the deep intraslab earthquakes to occur anywhere along the subduction zone increases the peak ground acceleration hazard near Portland, Oregon by about 20%. Alternative attenuation relations for deep earthquakes can result in ground motions that differ by a factor of two. The hazard uncertainty for the plate interface and intraslab earthquakes is analyzed through a Monte-Carlo logic tree approach and indicates a seismic hazard exceeding 1 g (0.2 s spectral acceleration) consistent with the U.S. National Seismic Hazard Maps in western Washington, Oregon, and California and an overall coefficient of variation that ranges from 0.1 to 0.4. Sensitivity studies indicate that the paleoseismic chronology and the magnitude of great plate interface earthquakes contribute significantly to the hazard uncertainty estimates for this region. Paleoseismic data indicate that the mean earthquake recurrence interval for great earthquakes is about 500 years and that it has been 300 years since the last great earthquake. We calculate the probability of such a great earthquake along the Cascadia plate interface to be about 14% when considering a time-dependent model and about 10% when considering a time-independent Poisson model during the next 50-year interval.     

俯冲带地震动特征及其衰减规律探讨

随着我国南海不断开发建设,海洋工程的抗震问题日益受到重视.我国南海东部区域位于大陆板块与海洋板块共同作用的俯冲带地区,地震活动频繁,震级较大,潜在地震对南海开发建设有重要影响.为了研究俯冲带地震的地震动特征及其衰减规律,本文基于实际俯冲带地震数据,并结合数值模拟方法,分析和探讨了俯冲带板内、板缘地震与浅地壳地震的地震动特征和衰减规律的差异.研究结果表明:俯冲带地震动存在区域性差异,在地震动衰减特征方面,同一区域的俯冲带板缘地震要比浅地壳地震衰减慢,俯冲带板内地震要比浅地壳地震衰减得快;数值模拟分析不同深度海水对海底地震动的影响表明,海底地震动水平分量几乎不受海水介质的影响,但是竖向分量随海水深度的增加有减小的趋势.最终,基于数值模拟和经验关系的混合方法建立了南海俯冲带地震动衰减关系模型,其结果可为海域区划等相关研究和海域工程建设提供参考.     

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

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

Peak ground accelerations from large (<Emphasis Type="Bold">M</Emphasis> ≥ 7.2) shallow crustal earthquakes: a comparison with predictions from eight recent ground-motion models

Strong-motion data from large (M ≥ 7.2) shallow crustal earthquakes invariably make up a small proportion of the records used to develop empirical ground motion prediction equations (GMPEs). Consequently GMPEs are more poorly constrained for large earthquakes than for small events. In this article peak ground accelerations (PGAs) observed in 38 earthquakes worldwide with M ≥ 7.2 are compared with those predicted by eight recent GMPEs. Well over half of the 38 earthquakes were not considered when deriving these GMPEs but the data were identified by a thorough literature review of strong-motion reports from the past 60 years. These data are provided in an electronic supplement for future investigations on ground motions from large earthquakes. The addition of these data provides better constraint of the between-event ground-motion variability in large earthquakes. It is found that the eight models generally provide good predictions for PGAs from these earthquakes, although there is evidence for slight under- or over-prediction of motions by some models (particularly for M > 7.6). The between-event variabilities predicted by most models match the observed variability, if data from two events (2001 Bhuj and 2005 Crescent City) that are likely atypical of earthquakes in active regions are excluded. For some GMPEs there is evidence that they are over-predicting PGAs in the near-source region of large earthquakes as well as over-predicting motions on hard rock. Overall, however, all the considered models, despite having been derived using limited data, provide reliable predictions of PGAs in the largest crustal earthquakes.     

Assessment of ground motion variability and its effects on seismic hazard analysis: a case study for iceland

Probabilistic seismic hazard analysis (PSHA) generally relies on the basic assumption that ground motion prediction equations (GMPEs) developed for other similar tectonic regions can be adopted in the considered area. This implies that observed ground motion and its variability at considered sites could be modelled by the selected GMPEs. Until now ground-motion variability has been taken into account in PSHA by integrating over the standard deviation reported in GMPEs, which significantly affects estimated ground motions, especially at very low probabilities of exceedance. To provide insight on this issue, ground-motion variability in the South Iceland Seismic Zone (SISZ), where many ground-motion records are available, is assessed. Three statistical methods are applied to separate the aleatory variability into source (inter-event), site (inter-site) and residual (intra-event and intra-site) components. Furthermore, the current PSHA procedure that makes the ergodic assumption of equality between spatially and temporal variability is examined. In contrast to the ergodic assumption, several recent studies show that the observed ground-motion variability at an individual location is lower than that implied by the standard deviation of a GMPE. This could imply a mishandling of aleatory uncertainty in PSHA by ignoring spatial variability and by mixing aleatory and epistemic uncertainties in the computation of sigma. Station correction coefficients are introduced in order to capture site effects at different stations. The introduction of the non-ergodic assumption in PSHA leads to larger epistemic uncertainty, although this is not the same as traditional epistemic uncertainty modelled using different GMPEs. The epistemic uncertainty due to the site correction coefficients (i.e. mean residuals) could be better constrained for future events if more information regarding the characteristics of these seismic sources and path dependence could be obtained.     

Characteristics of horizontal ground motion measures along principal directions

<正>Ground motion records are often used to develop ground motion prediction equations(GMPEs) for a randomly oriented horizontal component,and to assess the principal directions of ground motions based on the Arias intensity tensor or the orientation of the major response axis.The former is needed for seismic hazard assessment,whereas the latter can be important for assessing structural responses under multi-directional excitations.However,a comprehensive investigation of the pseudo-spectral acceleration(PSA) and of GMPEs conditioned on different axes is currently lacking.This study investigates the principal directions of strong ground motions and their relation to the orientation of the major response axis, statistics of the PSA along the principal directions on the horizontal plane,and correlation of the PSA along the principal directions on the horizontal plane.For these,three sets of strong ground motion records,including intraplate California earthquakes,inslab Mexican earthquakes,and interface Mexican earthquakes,are used.The results indicate that one of the principal directions could be considered as quasi-vertical.By focusing on seismic excitations on the horizontal plane,the statistics of the angles between the major response axis and the major principal axis are obtained;GMPEs along the principal axes are provided and compared with those obtained for a randomly oriented horizontal component;and statistical analysis of residuals associated with GMPEs along the principal directions is carried out.     

Criteria for Selecting and Adjusting Ground-Motion Models for Specific Target Regions: Application to Central Europe and Rock Sites

A vital component of any seismic hazard analysis is a model for predicting the expected distribution of ground motions at a site due to possible earthquake scenarios. The limited nature of the datasets from which such models are derived gives rise to epistemic uncertainty in both the median estimates and the associated aleatory variability of these predictive equations. In order to capture this epistemic uncertainty in a seismic hazard analysis, more than one ground-motion prediction equation must be used, and the tool that is currently employed to combine multiple models is the logic tree. Candidate ground-motion models for a logic tree should be selected in order to obtain the smallest possible suite of equations that can capture the expected range of possible ground motions in the target region. This is achieved by starting from a comprehensive list of available equations and then applying criteria for rejecting those considered inappropriate in terms of quality, derivation or applicability. Once the final list of candidate models is established, adjustments must be applied to achieve parameter compatibility. Additional adjustments can also be applied to remove the effect of systematic differences between host and target regions. These procedures are applied to select and adjust ground-motion models for the analysis of seismic hazard at rock sites in West Central Europe. This region is chosen for illustrative purposes particularly because it highlights the issue of using ground-motion models derived from small magnitude earthquakes in the analysis of hazard due to much larger events. Some of the pitfalls of extrapolating ground-motion models from small to large magnitude earthquakes in low seismicity regions are discussed for the selected target region.     

不同类型俯冲带地震竖向地震动阻尼修正系数对比研究

俯冲带地区竖向地震动的阻尼修正系数在工程结构抗震设计中起着重要作用。由于俯冲带地区的板块构造复杂,俯冲带地区的地震可划分为浅壳上地幔地震、板内地震和板间地震3种类型。为研究不同类型俯冲带地震的竖向地震动阻尼修正系数间是否具有显著差异而需要分别建立不同的阻尼修正系数模型,采用日本俯冲带地区的地震动数据,通过假设检验和构造差异指标的方式对不同类型地震的竖向加速度和位移反应谱的阻尼修正系数进行两两比较。结果显示:不同地震类型的竖向地震动阻尼修正系数在众多谱周期上存在统计意义和工程实际意义上的显著差异。该研究表明:研究俯冲带地区竖向地震动阻尼修正系数时需要考虑地震类型的影响。     

Inferred ground motions on Guadeloupe during the 2004 Les Saintes earthquake

Accurate estimates of the ground motions that occurred during damaging earthquakes are a vital part of many aspects of earthquake engineering, such as the study of the size and cause of the uncertainties within earthquake risk assessments. This article compares a number of methods to estimate the ground shaking that occurred on Guadeloupe (French Antilles) during the 21st November 2004 (M _w 6.3) Les Saintes earthquake, with the aim of providing more accurate shaking estimates for the investigation of the sources of uncertainties within loss evaluations, based on damage data from this event. The various techniques make differing use of the available ground-motion recordings of this earthquake and by consequence the estimates obtained by the different approaches are associated with differing uncertainties. Ground motions on the French Antilles are affected by strong local site effects, which have been extensively investigated in previous studies. In this article, use is made of these studies in order to improve the shaking estimates. It is shown that the simple methods neglecting the spatial correlation of earthquake shaking lead to uncertainties similar to those predicted by empirical ground-motion models and that these are uniform across the whole of Guadeloupe. In contrast, methods (such as the ShakeMap approach) that take account of the spatial correlation in motions demonstrate that shaking within roughly 10 km of a recording station (covering a significant portion of the investigated area) can be defined with reasonable accuracy but that motions at more distant points are not well constrained.     

Ground-motion prediction equations for the intermediate depth Vrancea (Romania) earthquakes

We present the regional ground-motion prediction equations for peak ground acceleration (PGA), peak ground velocity (PGV), pseudo-spectral acceleration (PSA), and seismic intensity (MSK scale) for the Vrancea intermediate depth earthquakes (SE-Carpathians) and territory of Romania. The prediction equations were constructed using the stochastic technique on the basis of the regional Fourier amplitude spectrum (FAS) source scaling and attenuation models and the generalised site amplification functions. Values of considered ground motion parameters are given as the functions of earthquake magnitude, depth and epicentral distance. The developed ground-motion models were tested and calibrated using the available data from the large Vrancea earthquakes. We suggest to use the presented equations for the rapid estimation of seismic effect after strong earthquakes (Shakemap generation) and seismic hazard assessment, both deterministic and probabilistic approaches.     

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.     

Seismic hazard assessment for the Itoiz dam site (Western Pyrenees,Spain)

This work summarises the seismic hazard analysis performed for the complete characterisation of strong ground-motion at the site of the Itoiz dam (Western Pyrenees, Spain). The hazard analysis includes the compilation of a composite catalogue from French and Spanish agencies, the definition of an original hybrid seismogenic source model (including zones and major faults) and the selection of ground motion prediction equations (GMPEs). Hazard results are provided as hazard curves and acceleration response spectra on rock for the 1000- and 5000-year return periods, which correspond respectively to the operating basis earthquake (OBE) and safety evaluation earthquake (SEE). The impact of truncating GMPEs at a number of standard deviations (epsilon) has been found not critical here for the return periods targeted. Subsequently, an analysis of the contribution of each source to total hazard and a hazard disaggregation analysis are performed in order to establish the earthquake-source parameters for both the OBE and SEE scenarios consistently with the seismotectonics of the region. The European Strong Motion database is then searched and a selection of records is proposed for each of the scenarios. Our results suggest that seismic hazard in the region is underestimated by the official Spanish seismic hazard map included in the current version of the code (NCSE-02), which is the reference document for the definition of seismic actions for dam projects in the whole Pyrenees.     

Post-earthquake estimates of different ground motion intensity measures for the 2012 Emilia earthquake

Estimates of the earthquake ground motion intensity over a geographical area have multiple uses, that is, emergency management, civil protection and seismic fragility assessment. In particular, with reference to fragility assessment, it is of interest to have estimates of the values of different ground-motion intensity measures in order to correlate them with the observed damage. To this purpose, the present paper uses a procedure recently proposed in the literature to estimate the ground-motion intensity for the 2012 Emilia mainshocks, considering different ground motion intensity measures and directionality effects. Ground motion prediction equations based on different site effect models, and spatial correlation models are calibrated for the Emilia earthquakes. The paper discusses the accuracy of the shakemaps obtained using the different soil effect models considered and presents the obtained shakemaps as supplementary material. The procedure presented in the paper is aimed at providing ground motion intensity values for seismic fragility assessment and is not intended as a tool to estimate shakemaps for rapid emergency assessment.     

Sensitivity of ground-motion scenarios to earthquake source parameters in the Tehran metropolitan area,Iran

The northern Tehran fault (NTF) is a principal active fault of the Alborz mountain belt in the northern Iran. The fault is located north of the highly populated Metropolitan Area of Tehran. Historical records and paleoseismological studies have shown that the NTF poses a high seismic risk for the Tehran region and the surrounding cities (e.g. Karaj). A series of ground-motion simulations are carried out using a hybrid kinematic-stochastic model to calculate broadband (0.1–20 Hz) ground-motion time histories for deterministic earthquake scenarios (M7.2) on the NTF. We will describe the source characteristics of the target event to develop a list of scenario earthquakes that are probably similar to a large earthquake on the NTF. The effect of varying different rupture parameters such as rupture velocity and rise time on the resulting broadband strong motions has been investigated to evaluate the range of uncertainty in seismic scenarios. The most significant parameters in terms of ground-shaking level are the rise time and the value of the rupture velocity. For the worst-case scenario, the maximum expected horizontal acceleration, and velocity at rock sites in Tehran range between 128 and 1315 cm/s/s and 11–191 cm/s, respectively. For the lowest scenario, the corresponding values range between 102 and 776 cm/s/s and 12 to 81 cm/s. Nonlinear soil effects may change these results but are not accounted for in this study. The largest variability of ground motion is observed in neighborhood of asperity and also in the direction of rupture propagation. The calculated standard deviation of all ground-motion scenarios is less than 30% of the mean. The capability of the simulation method to synthesize expected ground motions and the appropriateness of the key parameters used in the simulations are confirmed by comparing the synthetic peak ground motions (PGA, PGV and response spectra) with empirical ground-motion prediction equations.     

Investigating strong ground-motion variability using analysis of variance and two-way-fit plots

A statistical method to quantitatively assess the relative importance of unmodelled site and source effects on the observed variability (σ) in ground motions is presented. The method consists of analysis of variance (ANOVA) using the computed residuals with respect to an empirical ground-motion model for strong-motion records of various earthquakes recorded at a common set of stations. ANOVA divides the overall variance (σ ²) into the components due to site and source effects (respectively σ _S ² and σ _E ²) not modelled by the ground-motion model plus the residual variance not explained by these effects (σ _R ²). To test this procedure, four sets of observed strong-motion records: two from Italy (Umbria-Marche and Molise), one from the French Antilles and one from Turkey, are used. It is found that for the data from Italy, the vast majority of the observed variance is attributable to unmodelled site effects. In contrast, the variation in ground motions in the French Antilles and Turkey data is largely attributable, especially at short periods, to source effects not modelled by the ground-motion estimation equations used.     

Ground-motion characterization for the probabilistic seismic hazard assessment in Turkey

This study describes the methodology implemented to establish the ground-motion logic-tree for national probabilistic seismic hazard map of Turkey for shallow active crustal regions. The presented procedure provides quantitative information to guide the hazard experts while establishing the logic tree to capture the epistemic uncertainty in ground-motion characterization. It uses non-data-driven and data-driven testing methods to identify and rank candidate ground-motion prediction equations (GMPEs) under a specific ground-motion database. The candidate GMPEs are subjected to visual inspection and are classified into center, body and range (CBR) spectral estimates for a proper consideration of epistemic uncertainty. The GMPEs classified into CBR are then used in a suite of seismic hazard sensitivity analysis to establish the most suitable GMPE logic-tree whose spectral estimates are not biased by any one of the GMPEs in the logic-tree structure. The sensitivity analysis considers normalized spectral ordinates and is not manipulated by the spectral amplitudes. The proposed procedure is inherited from the relevant studies of the Earthquake Model of the Middle East (EMME; www.efehr.org:8080/jetspeed/portal/emme.psml) regional seismic hazard project. This paper also highlights the similarities and differences in ground-motion characterization between EMME and our approach.