Robust attenuation relations for peak time-domain parameters of strong ground motions

This study presents new attenuation models for the estimation of peak ground acceleration (PGA), peak ground velocity (PGV), and peak ground displacement (PGD) using a hybrid method coupling genetic programming and simulated annealing, called GP/SA. The PGA, PGV, and PGD were formulated in terms of earthquake magnitude, earthquake source to site distance, average shear-wave velocity, and faulting mechanisms. A worldwide database of strong ground motions released by Pacific Earthquake Engineering Research Center (PEER) was employed to establish the models. A traditional genetic programming analysis was performed to benchmark the proposed models. For more validity verification, the GP/SA models were employed to predict the ground-motion parameters of the Iranian plateau earthquakes. Sensitivity and parametric analyses were carried out and discussed. The results show that the GP/SA attenuation models can offer precise and efficient solutions for the prediction of estimates of the peak time-domain characteristics of strong ground motions. The performance of the proposed models is better than or comparable with the attenuation relationships found in the literature.     

Hybrid attenuation model for estimation of peak ground accelerations in the Kutch region, India

The Kutch region of Gujarat in India is the locale of one of the most devastating earthquake of magnitude (M _w) 7.7, which occurred on January 26, 2001. Though, the region is considered as seismically active region, very few strong motion records are available in this region. First part of this paper uses available data of strong motion earthquakes recorded in this region between 2006 and 2008 years to prepare attenuation relation. The developed attenuation relation is further used to prepare synthetic strong motion records of large magnitude earthquakes using semiempirical simulation technique. Semiempirical simulation technique uses attenuation relation to simulate strong ground motion records of any target earthquake. The database of peak ground acceleration obtained from simulated records is used together with database of peak ground acceleration obtained from observed record to develop following hybrid attenuation model of wide applicability in the Kutch region: $$ \begin{aligned} \ln \left( {\text{PGA}} \right) & = - 2.56 + 1.17 \, M_{\text{w}} - \, 0.015R - 0.0001\ln \left( {E + 15} \right) \\ &\quad 3.0 \le M_{\text{w}} \le 8.2;\quad 12 \le R \le 120;\quad {\text{std}} . {\text{ dev}}.(\sigma ): \pm 0.5 \\ \end{aligned} $$ ln ( PGA ) = ? 2.56 + 1.17 M w ? 0.015 R ? 0.0001 ln ( E + 15 ) 3.0 ≤ M w ≤ 8.2 ; 12 ≤ R ≤ 120 ; std . dev . ( σ ) : ± 0.5 In the above equation, PGA is maximum horizontal ground acceleration in gal, M _w is moment magnitude of earthquake, R is hypocentral distance, and E is epicentral distance in km. The standard deviation of residual of error in this relation is 0.5. This relation is compared with other available relations in this region, and it is seen that developed relation gives minimum root mean square error in comparison with observed and calculated peak ground acceleration from same data set. The applicability of developed relation is further checked by testing it with the observed peak ground acceleration from earthquakes of magnitude (M _w), 3.6, 4.0, 4.4, and 7.7, respectively, which are not included in the database used for regression analysis. The comparison demonstrates the efficacy of developed hybrid attenuation model for calculating peak ground acceleration values in the Kutch region.     

Generation of isoseismal maps for Mexico using velocity and acceleration data

Acceleration and velocity data of large earthquakes recorded by the Mexican National Broadband Seismological Network were used to generate isoseismal maps for Mexico. The seismic data consist of 99 events recorded at 27 seismic stations located in the southern and northern regions of Mexico from 2004 to 2009. The magnitude (Mw) of these events ranged from 4.1 to 7.1. Peak ground velocity values (PGV) and peak ground accelerations (PGA) were estimated, and velocity-derived peak ground accelerations (PGAv) were calculated. No important differences between PGAv and PGA values were found; thus, both parameters were used in the isoseismal determination. The generated synthetic isoseismal maps were compared with those from an existing catalog for large earthquakes in Mexico. Using empiric relations between PGA and MMI (modified Mercalli intensity) and PGV and MMI, the obtained simulated intensity maps showed similar behavior to those reported in the catalog. The results indicate that the PGAv values can be used to determine intensities when acceleration records (PGA) are unavailable.     

Surface ruptures induced by the devastating 1068 AD earthquake in the southern Arava valley, Dead Sea Rift, Israel

The Elat fault (a segment of the Dead Sea Transform) runs along the southern Arava valley (part of the Dead Sea Rift, Israel) forming a complex fault zone that displays a time-dependent seismic behaviour. Paleoseismic evidence shows that this fault zone has generated at least 15 earthquakes of magnitude larger than M 6 during the late Pleistocene and the Holocene. However, at present the Elat fault is one of the quietest segments of the Dead Sea Transform, lacking even microsesimicity. The last event detected in the southern Arava valley occurred in the Avrona playa and was strong enough to have deformed the playa and to change it from a closed basin with internal drainage into an open basin draining to the south.Paleoseismological, geophysical and archaeological evidences indicate that this event was the historical devastating earthquake, which occurred in 1068 AD in the eastern Mediterranean region. According to the present study this event was strong enough to rupture the surface, reactivate at least two fault branches of the Elat fault and vertically displace the surface and an early Islamic irrigation system by at least 1 m. In addition, the playa area was uplifted between 2.5 and 3 m along the eastern part of the Elat fault shear zone. Such values are compatible with an earthquake magnitude ranging between M 6.6 and 7. Since the average recurrence interval of strong earthquakes during the Holocene along the Elat fault is about 1.2 ± 0.3 ky and the last earthquake occurred more about 1000 years ago, the possibility of a very strong earthquake in this area in the future should be seriously considered in assessing seismic hazards.     

Recent Advances in Prediction and Processing of Strong Ground Motions

We present an overview of our recent results on utilizing small earthquakes in the earthquake engineering practice. Site-specific ground motion time-histories of large earthquakes can be successfully simulated using recordings of small earthquakes which are often referred to as 'empirical Green's functions' in seismology. Another important practical problem is whether and how these observations can be used in seismic risk studies which are based on empirical attenuation relations for ground motion parameters. We study a possibility of extrapolating attenuation relations for small earthquakes, to larger magnitudes using the data from the Garner Valley downhole array in Southern California. Finally we introduce efficient ground motion processing techniques in frequency- and time-domains and apply them to site response estimation.     

Sensitivity analysis of seismic hazard for the northwestern portion of the state of Gujarat, India

We test the sensitivity of seismic hazard to three fault source models for the northwestern portion of Gujarat, India. The models incorporate different characteristic earthquake magnitudes on three faults with individual recurrence intervals of either 800 or 1600 years. These recurrence intervals imply that large earthquakes occur on one of these faults every 266–533 years, similar to the rate of historic large earthquakes in this region during the past two centuries and for earthquakes in intraplate environments like the New Madrid region in the central United States. If one assumes a recurrence interval of 800 years for large earthquakes on each of three local faults, the peak ground accelerations (PGA; horizontal) and 1-Hz spectral acceleration ground motions (5% damping) are greater than 1 g over a broad region for a 2% probability of exceedance in 50 years' hazard level. These probabilistic PGAs at this hazard level are similar to median deterministic ground motions. The PGAs for 10% in 50 years' hazard level are considerably lower, generally ranging between 0.2 g and 0.7 g across northwestern Gujarat. Ground motions calculated from our models that consider fault interevent times of 800 years are considerably higher than other published models even though they imply similar recurrence intervals. These higher ground motions are mainly caused by the application of intraplate attenuation relations, which account for less severe attenuation of seismic waves when compared to the crustal interplate relations used in these previous studies. For sites in Bhuj and Ahmedabad, magnitude (M) 7 3/4 earthquakes contribute most to the PGA and the 0.2- and 1-s spectral acceleration ground motion maps at the two considered hazard levels.     

The first peak ground motion attenuation relationships for North of Vietnam

The first attenuation relationships of peak ground acceleration (PGA) and peak ground velocity (PGV) for northern Vietnam are obtained in this study. Ground motion data are collected by a portable broadband seismic network in northern Vietnam as a part of cooperation between the Institute of Geophysics, Vietnamese Academy of Science and Technology, Vietnam and Institute of Earth Sciences, Academia Sinica, Taiwan. The database comprises a total of 330 amplitude records by 14 broadband stations from 53 shallow earthquakes, which were occurred in and around northern Vietnam in the period between 01/2006 and 12/2009. These earthquakes are of local magnitudes between 1.6 and 4.6, focal depths less than 30 km, and epicentral distances less than 500 km. The new attenuation relationships for PGA and PGV are:

log₁₀(PGA)=-0.987+0.7521M_L-log₁₀(R)-0.00475R,     

基于复合方法的天水盆地宽频带地震动模拟

本文采用有限差分和随机振动合成结合的复合方法,模拟了当礼县—罗家堡断裂发生矩震级Mw7.7级大地震时,在天水盆地产生的宽频带地震动场,分析了在设定地震条件下盆地内的地震动分布特征,为该区黄土地震滑坡分析提供了地震动参数结果。结果显示:(1)有限差分法和随机振动合成法可以很好地互补,得到盆地内地表宽频带地震动;(2)地震在盆地区域产生了强烈地震动,PGA(峰值加速度)介于150~900 gal,离断层较近的区域东南角的PGA最大,随着断层距的增加,PGA逐渐减小。河谷南侧的PGA值相比北侧较大,具备诱发滑坡的强大动力条件;(3)盆地区域PGV(峰值速度)最大为120 cm/s。受第四系覆盖层放大效应和地形放大效应共同影响,水平向地震动在盆地区域东侧和中部具有较大PGV,而西侧PGV相对较小。竖向地震动在盆地区域东侧较弱,而在中部和西侧较强,特别是最西侧陡峭的山坡上,PGV达到了最大值。此外,竖向地震动明显受到覆盖层厚度的影响,譬如在盆地区域南侧的中间部位,也具有较大的PGV。     

A new probabilistic seismic hazard study of Spain

In this paper, we present a probabilistic seismic hazard analysis (PSHA) for mainland Spain that takes into account recent new results in seismicity, seismic zoning, and strong ground attenuation not considered in the latest PSHA of the Spanish Building Code. Those new input data have been obtained as a three-step project carried out in order to improve the existing hazard map for mainland Spain. We have produced a new earthquake catalogue for the area, in which the earthquakes are given in moment magnitude through specific deduced relationships for our territory based on intensity data (Mezcua et al. in Seismol Res Lett 75:75–81, 2004). In addition, we included a new seismogenetic zoning based on the recent partial zoning studies performed by different authors. Finally, as we have developed a new strong ground motion model for the area García Blanco (2009), it was considered in the hazard calculation together with other attenuations gathered from different authors using data compatible with our region. With this new data, a logic tree process is defined to quantify the epistemic uncertainty related to those parts of the process. A sensitivity test has been included in order to analyze the different models of ground motion and seismotectonic zonation used in this work. Finally, after applying a weighting scheme, a mean hazard map for PGA, based on rock type condition for 10% exceedance probability in 50 years, is presented, including 15th and 85th percentile hazard maps. The main differences with the present official building code hazard map are analyzed.     

Probabilistic assessment of seismic hazard of dam sites in Jordan

A Probabilistic method is used to evaluate the seismic hazard of nineteen embankment dam sites in Jordan. A line source model developed by McGuire (1978) is used in this study. An updated earthquake catalogue covering the period from 1 A.D. to 1991 A.D. is used for this purpose. This catalogue includes all earthquakes that occurred in Jordan and adjacent areas, more specifically between latitudes 27.0°–35.5° N and longitudes 32.0°–39.0° E.Nine distinct seismic sources of potential seismic activities are identified. The seismic hazard parameters are determined using the method suggested by Kijko and Sellevoll (1989).The Peak Ground Acceleration (PGA) is selected as a measure of ground motion severity. Esteva (1974) attenuation relationship is used in evaluating PGA values at each dam site. Analysis is carried out for 50%, 90%, and 95% probability that is not being exceeded in a life time of 50, 100, and 200 years.Results of analysis indicate that PGA values are higher for dam sites closer to the Dead Sea Fault. This fault is believed to be responsible for most earthquake activities in Jordan and vicinity. The highest PGA value is found to be for Al-Karama dam site.     

http://www.sciencedirect.com/science/article/pii/S1674987114001376

A new model is derived to predict the peak ground acceleration(PGA) utilizing a hybrid method coupling artificial neural network(ANN) and simulated annealing(SA), called SA-ANN. The proposed model relates PGA to earthquake source to site distance, earthquake magnitude, average shear-wave velocity,faulting mechanisms, and focal depth. A database of strong ground-motion recordings of 36 earthquakes,which happened in Iran's tectonic regions, is used to establish the model. For more validity verification,the SA-ANN model is employed to predict the PGA of a part of the database beyond the training data domain. The proposed SA-ANN model is compared with the simple ANN in addition to 10 well-known models proposed in the literature. The proposed model performance is superior to the single ANN and other existing attenuation models. The SA-ANN model is highly correlated to the actual records(R=0.835 and r =0.0908) and it is subsequently converted into a tractable design equation.     

Tectonic stress,seismicity, and seismic hazard in the southeastern Carpathians

Intermediate-depth earthquakes in the Vrancea region occur in response to stress generation due to descending lithosphere beneath the southeastern Carpathians. In this article, tectonic stress and seismicity are analyzed in the region on the basis of a vast body of observations. We show a correlation between the location of intermediate-depth earthquakes and the predicted localization of maximum shear stress in the lithosphere. A probabilistic seismic hazard assessment (PSHA) for the region is presented in terms of various ground motion parameters on the utilization of Fourier amplitude spectra used in engineering practice and risk assessment (peak ground acceleration, response spectra amplitude, and seismic intensity). We review the PSHA carried out in the region, and present new PSHA results for the eastern and southern parts of Romania. Our seismic hazard assessment is based on the information about the features of earthquake ground motion excitation, seismic wave propagation (attenuation), and site effect in the region. Spectral models and characteristics of site-response on earthquake ground motions are obtained from the regional ground motion data including several hundred records of small and large earthquakes. Results of the probabilistic seismic hazard assessment are consistent with the features of observed earthquake effects in the southeastern Carpathians and show that geological factors play an important part in the distribution of the earthquake ground motion parameters.     

Probabilistic seismic hazard analysis for Sumatra, Indonesia and across the Southern Malaysian Peninsula

The ground motion hazard for Sumatra and the Malaysian peninsula is calculated in a probabilistic framework, using procedures developed for the US National Seismic Hazard Maps. We constructed regional earthquake source models and used standard published and modified attenuation equations to calculate peak ground acceleration at 2% and 10% probability of exceedance in 50 years for rock site conditions. We developed or modified earthquake catalogs and declustered these catalogs to include only independent earthquakes. The resulting catalogs were used to define four source zones that characterize earthquakes in four tectonic environments: subduction zone interface earthquakes, subduction zone deep intraslab earthquakes, strike-slip transform earthquakes, and intraplate earthquakes. The recurrence rates and sizes of historical earthquakes on known faults and across zones were also determined from this modified catalog. In addition to the source zones, our seismic source model considers two major faults that are known historically to generate large earthquakes: the Sumatran subduction zone and the Sumatran transform fault. Several published studies were used to describe earthquakes along these faults during historical and pre-historical time, as well as to identify segmentation models of faults. Peak horizontal ground accelerations were calculated using ground motion prediction relations that were developed from seismic data obtained from the crustal interplate environment, crustal intraplate environment, along the subduction zone interface, and from deep intraslab earthquakes. Most of these relations, however, have not been developed for large distances that are needed for calculating the hazard across the Malaysian peninsula, and none were developed for earthquake ground motions generated in an interplate tectonic environment that are propagated into an intraplate tectonic environment. For the interplate and intraplate crustal earthquakes, we have applied ground-motion prediction relations that are consistent with California (interplate) and India (intraplate) strong motion data that we collected for distances beyond 200 km. For the subduction zone equations, we recognized that the published relationships at large distances were not consistent with global earthquake data that we collected and modified the relations to be compatible with the global subduction zone ground motions. In this analysis, we have used alternative source and attenuation models and weighted them to account for our uncertainty in which model is most appropriate for Sumatra or for the Malaysian peninsula. The resulting peak horizontal ground accelerations for 2% probability of exceedance in 50 years range from over 100% g to about 10% g across Sumatra and generally less than 20% g across most of the Malaysian peninsula. The ground motions at 10% probability of exceedance in 50 years are typically about 60% of the ground motions derived for a hazard level at 2% probability of exceedance in 50 years. The largest contributors to hazard are from the Sumatran faults.     

Modeling of strong ground motions for 1991 Uttarkashi, 1999 Chamoli earthquakes, and a hypothetical great earthquake in Garhwal–Kumaun Himalaya

In this study, stochastic finite fault modeling is used to simulate Uttarkashi (1991) and Chamoli (1999) earthquakes using all available source, path, and site parameters available for the region. These two moderate earthquakes are recorded at number of stations of a strong motion network. The predicted peak ground accelerations at these stations are compared with the observed data and the ground motion parameters are constrained. The stress drop of Uttarkashi and Chamoli earthquakes is constrained at 77 and 65?bars, respectively, whereas the quality factor Q _C is 112 $ f^{0.97} $ and 149 $ f^{0.95} $ for these two regions. The high-frequency attenuation parameter Kappa is in the range 0.04?C0.05. The constrained ground motion parameters are then used to simulate Mw 8.5 earthquake in central seismic gap region of Himalaya. Two scenarios are considered with epicenter of future great earthquake at locations of Uttarkashi and Chamoli earthquakes using above constrained parameters. The most vulnerable towns are the towns of Dehradun and Almora where expected PGA is in excess of 600?cm/s² at V_S30 520?m/s when the epicenter of the great earthquake is at the location of Uttarkashi (1991) earthquake. The towns of Shimla and Chandigarh can expect PGA close to 200?cm/s². Whereas when the epicenter of the great earthquake is at the location of Chamoli (1999) earthquake, the towns of Dehradun and Almora can expect PGA of around 500 and 400?cm/s², respectively, at V_S30 620?m/s. The National Capital Region, Delhi can expect accelerations of around 80?cm/s² in both the cases. The PGA contour maps obtained in this study can be used to assess the seismic hazard of the region and identify vulnerable areas in and around central Himalaya from a future great earthquake.     

An attenuation relationship based on Turkish strong motion data and iso-acceleration map of Turkey

This paper presents an attenuation relationship of peak ground acceleration (PGA) derived from Turkish strong motion data for rock, soil and soft soil sites and an iso-acceleration map of Turkey based on this relationship. For the purpose, among all the three-component accessible records, 221 records from 122 earthquakes that occurred in Turkey between 1976 and November 2003 were selected. The database was compiled for earthquakes with moment magnitudes (M_w) and PGA values ranging between 4.1 and 7.5, and 20 and 806 gal, and distances to epicenter considered in the database were between 5 and 100 km. From the regression analysis of the data, an attenuation equation of PGA considering rock, soil and soft soil conditions was developed. The PGA values predicted from the equation suggested in this study and those both from a few domestic equations and some imported equations were compared. In addition, an iso-acceleration map of Turkey was constructed using the suggested attenuation equation and considering both known active faults and epicenter locations of the earthquakes that have occurred in Turkey.     

Strong ground motion attenuation in Aswan area, Egypt

Recent and paleo seismicity indicate that moderate seismic activity is relatively large for Aswan area. This is a warning on the possibility of occurrence of earthquakes in the future too. No strong motion records are available in Aswan area for engineers to rely upon. Consequently, the seismological modeling is an alternative approach till sufficient instrumental records around Aswan become available. In the present study, we have developed new ground motion attenuation relationship for events spanning 4.0?≤? M _w?≤?7.0 and distance to the surface projection of the fault up to 100 km for Aswan based on a statistically simulated seismological model. We generated suites of ground motion time histories using stochastic technique. The ground motion attenuation relation describes the dependence of the strength of the ground motions on the earthquake magnitude and distance from the earthquake. The proposed equation for peak ground acceleration (PGA) for the bed rock is in the form of: $ {\mathbf{log}}{\text{ }}\left( {{\mathbf{PGA}}/{\mathbf{gal}}} \right){\text{ }} = {\mathbf{1}}.{\mathbf{24}} + {\mathbf{0}}.{\mathbf{358}}{M_{\mathbf{w}}} - {\text{ }}{\mathbf{log}}\left( {\mathbf{R}} \right){\text{ }}-{\text{ }}{\mathbf{0}}.{\mathbf{008}}{\text{ }}{\mathbf{R}}{\text{ }} + {\text{ }}{\mathbf{0}}.{\mathbf{22}}{\text{ }}{\mathbf{P}} $ . Where PGA is the peak ground acceleration in gal (cm/s²); M_w, its moment magnitude; R is the closest distance between the rupture projection and the site of interest; and the factor P is a dummy variable. It is observed that attenuation of strong motion in Aswan is correlated with those used before in Egypt.     

四川汶川地震诱发滑坡与峰值速度的关系

利用汶川地震诱发的滑坡资料和得到的地震动数据,对地震诱发滑坡与峰值速度的关系进行了分析,得到一些基本认识：①地震诱发滑坡与地震动峰值速度(PGV)存在明显的正相关性,可以利用PGV作为判别地震滑坡的判据;②汶川地震龙门山震区可以触发地震滑坡的PGV下限约为0.5m/s;③龙门山震区可以触发地震滑坡的PGV上限约为1.5m/s;④建立了汶川地震龙门山震区PGV衰减关系,可用以估计地震滑坡致灾的范围。利用地震动参数研究地震诱发滑坡具有很好的一致性,可以克服应用烈度时存在的不足。上述研究结果可以应用于震后灾害快速评估工作中,为应急救援方案的制定提供参考信息。     

Seismicity, seismotectonics and seismic hazard of Italy

A first generation of probabilistic seismic hazard maps of the Italian country are presented. They are based on seismogenic zoning deriving from a kinematic model of the structural tectonic units and on an earthquake catalogue with the foreshock and aftershock events filtered out. The following ground motion parameters have been investigated and mapped using attenuation equations based on strong-motion recordings of Italian earthquakes: peak ground acceleration and velocity; Arias intensity; strong motion duration; and the pseudovelocity and pseudoacceleration spectral values at 14 fixed frequencies both for the vertical and the largest horizontal component. A Poissonian model of earthquake occurrence is assumed as a default and the hazard maps are presented in terms of ground motion values expected to be exceeded at a 10% probability level in 50 years (return period 475 years) according to the requirement of Eurocode 8 for the seismic classification of national territories, as well as in terms of exceedance probabilities of selected ground motion values. Finally, as a tentative study, the use of hybrid methods (implementing both seismogenic zones and structures), renewal processes (including earthquake forecasting) and the influence of site effects (as the basis for the planning of earthquake scenarios) were explored.     

四川汶川地震诱发滑坡与峰值速度的关系

利用汶川地震诱发的滑坡资料和得到的地震动数据,对地震诱发滑坡与峰值速度的关系进行了分析,得到一些基本认识：①地震诱发滑坡与地震动峰值速度(PGV)存在明显的正相关性,可以利用PGV作为判别地震滑坡的判据;②汶川地震龙门山震区可以触发地震滑坡的PGV下限约为0.5m/s;③龙门山震区可以触发地震滑坡的PGV上限约为1.5m/s;④建立了汶川地震龙门山震区PGV衰减关系,可用以估计地震滑坡致灾的范围。利用地震动参数研究地震诱发滑坡具有很好的一致性,可以克服应用烈度时存在的不足。上述研究结果可以应用于震后灾害快速评估工作中,为应急救援方案的制定提供参考信息。