Simulation of Strong Ground Motion During the 1950 Great Assam Earthquake

In this paper, ground motion during the Independence Day earthquake of August 15, 1950 (M_w 8.6, Ben-Menahem et al., 1974) in the northeastern part of India is estimated by seismological approaches. A hybrid simulation technique which combines the low frequency ground motion simulated from an analytical source mechanism model with the stochastically simulated high-frequency components is used for obtaining the acceleration time histories. A series of ground motion simulations are carried out to estimate the peak ground acceleration (PGA) and spectral accelerations at important cities and towns in the epicentral region. One sample PGA distribution in the epicentral region encompassing the epicenter is also obtained. It is found that PGA in the epicentral region has exceeded 1 g during this earthquake. The estimated PGA’s are validated to the extent possible using the MMI values. The simulated acceleration time histories can be used for the assessment of important engineering structures in northeastern India.     

Ground motion parameters of Shillong plateau: One of the most seismically active zones of northeastern India

Strong ground motion parameters for Shillong plateau of northeastern India are examined. Empirical relations are obtained for main parameters of ground motions as a function of earthquake magnitude, fault type, source depth, velocity characterization of medium and distance. Correlation between ground motion parameters and characteristics of seismogenic zones are established. A new attenuation relation for peak ground acceleration is developed, which predicts higher expected PGA in the region. Parameters of ...     

Strong Ground Motion Estimation During the Kutch, India Earthquake

In the absence of strong motion records, ground motion during the 26^th January, 2001 Kutch, India earthquake, has been estimated by analytical methods. A contour map of peak ground acceleration (PGA) values in the near source region is provided. These results are validated by comparing them with spectral response recorder data and field observations. It is found that very near the epicenter, PGA would have exceeded 0.6 g. A set of three aftershock records have been used as empirical Green's functions to simulate ground acceleration time history and 5% damped response spectrum at Bhuj City. It is found that at Bhuj, PGA would have been 0.31 g–0.37 g. It is demonstrated that source mechanism models can be effectively used to understand spatial variability of large-scale ground movements near urban areas due to the rupture of active faults.     

Probabilistic earthquake hazard assessment for Peninsular India

In this paper, a new probabilistic seismic hazard assessment (PSHA) is presented for Peninsular India. The PSHA has been performed using three different recurrence models: a classical seismic zonation model, a fault model, and a grid model. The development of a grid model based on a non-parameterized recurrence model using an adaptation of the Kernel-based method that has not been applied to this region before. The results obtained from the three models have been combined in a logic tree structure in order to investigate the impact of different weights of the models. Three suitable attenuation relations have been considered in terms of spectral acceleration for the stable continental crust as well as for the active crust within the Gujarat region. While Peninsular India has experienced large earthquakes, e.g., Latur and Jabalpur, it represents in general a stable continental region with little earthquake activity, as also confirmed in our hazard results. On the other hand, our study demonstrates that both the Gujarat and the Koyna regions are exposed to a high seismic hazard. The peak ground acceleration for 10 % exceedance in 50 years observed in Koyna is 0.4 g and in the Kutch region of Gujarat up to 0.3 g. With respect to spectral acceleration at 1 Hz, estimated ground motion amplitudes are higher in Gujarat than in the Koyna region due to the higher frequency of occurrence of larger earthquakes. We discuss the higher PGA levels for Koyna compared Gujarat and do not accept them uncritically.     

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

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

鲁甸 MS6.5级地震强震动记录及震害分析

2014年8月3日发生的鲁甸地震是我国继2013年芦山地震之后的又一次破坏性浅源地震,造成了严重的人员伤亡和工程结构破坏.本文处理了我国数字强震动台网捕获的60余组主震三分量强震动记录,绘制了震中附近区域的峰值加速度等值线图,长轴沿西北-东南方向展布.通过与中国西部常用衰减关系的对比,发现各模型的预测值均不同程度高估了峰值加速度、峰值速度的观测值.最后以水平向峰值加速度较大的三个典型台站为例,详细调查了附近建筑的破坏情况,结合宏观烈度分布结果,分析了地震动特征与震害的相关性.     

A time-domain nonlinear effective-stress non-Masing approach of ground response analysis of Guwahati city,India

The response of subsoil strata subjected to seismic excitations plays an important role in governing the response of the overlying superstructures at any site. Ground response analysis(GRA) helps to assess the influence of soil characteristics on the propagating seismic stress waves from the bedrock level to the ground surface during an earthquake. For the northeastern region of India, located in the highest seismic zone in the country, conducting an extensive GRA study is of prime importance. Conventionally, most of the GRA studies are carried out using the equivalent linear method, which, being a simplistic approach, cannot capture the nonlinear behavior of soil during seismic shaking. This paper presents the outcomes of a one-dimensional effective stress based nonlinear GRA conducted for Guwahati city(located in northeast India) incorporating the non-Masing load/unload/reload characteristics. The various ground response parameters evaluated from this study help in assessing the ground shaking, soil amplification, and site responses expected in this region. 2D contour maps, which are representative of the distribution of some of these parameters throughout Guwahati city, are also developed. The results presented herein can serve as guidelines for the design of foundations and superstructures in this region.     

Estimation of Strong Ground Motion from a Great Earthquake Mw 8.5 in Central Seismic Gap Region, Himalaya (India) Using Empirical Green’s Function Technique

In the present study ground motions for a Mw 8.5 scenario earthquake are estimated at 13 sites in Kumaun-Garhwal region using the empirical Green’s function technique. The recordings of 1991 Uttarkashi earthquake of Mw 6.8 at these sites are used as an element earthquake. A heterogeneous source model consisting of two asperities is considered for simulating the ground motions. The entire central seismic gap (CSG) can expect acceleration in excess of 100 cm/s² with NW portion in excess of 400 cm/s² and SE between 100 and 200 cm/s². The central portion can expect peak ground acceleration (PGA) between 200 and 400 cm/s². It has been observed from simulation of strong ground motion that sites located near the rupture initiation point can expect accelerations in excess of 1g. In the present analysis, Bhatwari and Uttarkashi can expect ground accelerations in excess of 1g. The estimates of the PGA are compared with earlier studies in the same region using different methodologies and it was found that the results are comparable. This has put constrains on the expected PGAs in this region. The obtained PGA values can be used in identifying the vulnerable areas in the central Himalaya, thereby facilitating the planning, design and construction of new structures and strengthening of the existing structures in the region.     

Seismic Zonation of the Delhi Region for Bedrock Ground Motion

An endeavor is made to compute peak ground horizontal accelerations at bedrock level in the Delhi region due to the seismogenic sources present around Delhi. The entire area is divided into six seismogenic sources for which seismic hazard analysis is carried out using the complete and extreme part of the seismicity data. Maximum likelihood estimates of hazard parameters viz., seismic activity rate , b value and maximum probable earthquake M _max are made for each zone. The return periods and the probabilities of occurrence of various magnitudes for return periods of 50, 100 and 1000 years are also computed for each zone. The peak ground acceleration (PGA) values for 20% exceedance in 50 years are then computed for the Delhi region from each zone. The maximum PGA value considering all the zones is 0.34 g, which is due to the Mathura fault zone. The seismogenic zones V and VI, i.e., Mathura fault zone and the Sohna fault zone are observed to be contributing maximum PGA values in the Delhi region governing the isoacceleration contours computed for the region. The seismic zonation map for the PGA values at the bedrock level is obtained for the Delhi region. This can be used directly as input for the microzonation of ground motion at the surface by incorporating the local site conditions.     

Effect of seismic super-shear rupture on the directivity of ground motion acceleration

The effect of seismic super-shear rupture on the directivity of ground motions using simulated accelerations of a vertical strike-slip fault model is the topic of this study. The discrete wave number/finite element method was adopted to calculate the ground motion in the horizontal layered half space. An analysis of peak ground acceleration （PGA） indicates that similar to the sub-shear situation, directivity also exists in the super-shear situation. However, there are some differences as tbllows： （1） The PGA of the fault-normal component decreases with super-shear velocity, and the areas that were significantly affected by directivity in the PGA field changed from a cone-shaped region in the forward direction in a sub-shear situation to a limited near-fault region in a super-shear situation. （2） The PGA of the fault-parallel and vertical component is not as sensitive as the fault-normal component to the increasing super-shear velocity. （3） The PGA of the fault-normal component is not always greater than the fault-parallel component when the rupture velocity exceeds the shear wave velocity.     

Constraining paleo PGA values by numerical analysis of overturned columns

The estimated peak ground acceleration (PGA) for a specific region cannot always be determined reliably from compilation and statistical analyses of instrumental data acquired because of the relatively short time window available for such an approach, typically up to 100 years of instrumentation. We propose a complimentary approach for estimating PGA in a specific region, through back analysis of seismically driven column collapse in historic monuments, using the numerical discrete element discontinuous deformation analysis (DDA) method. Preliminary threshold ‘paleo PGA’ values thus obtained constrain the lower bound of PGA estimates using information from a much broader time window, in the case study presented here of approximately 1200 years. Copyright © 2009 John Wiley & Sons, Ltd.     

川滇地区水平向强地震动衰减关系研究

利用近几年在川滇地区破坏性地震及汶川地震强余震中获得的951条加速度记录, 通过选取合适的衰减关系形式和距离参数, 经统计获得了川滇地区水平向峰值加速度与35个周期点对应的加速度反应谱(阻尼比0.05)衰减关系. 统计结果显示, 震级越大, 地震动随距离的衰减越慢; 加速度反应谱的短周期部分随距离的衰减明显快于长周期部分; 中强地震在近场也会产生较大的峰值加速度和较高的短周期加速度反应谱. 由于统计资料的非完整性, 本文统计结果尚不宜直接应用于相关工程场地的设计地震动参数确定中. 但本文结果与目前由转换方法得到的我国西部地震动参数衰减规律的较大差异提示: 进一步深入研究我国地震动参数衰减规律应是紧迫而重要的任务.      

Evaluation of site effects on ground motions based on equivalent linear site response analysis and liquefaction potential in Chennai,south India

We study local site effects with detailed geotechnical and geophysical site characterization to evaluate the site-specific seismic hazard for the seismic microzonation of the Chennai city in South India. A Maximum Credible Earthquake (MCE) of magnitude 6.0 is considered based on the available seismotectonic and geological information of the study area. We synthesized strong ground motion records for this target event using stochastic finite-fault technique, based on a dynamic corner frequency approach, at different sites in the city, with the model parameters for the source, site, and path (attenuation) most appropriately selected for this region. We tested the influence of several model parameters on the characteristics of ground motion through simulations and found that stress drop largely influences both the amplitude and frequency of ground motion. To minimize its influence, we estimated stress drop after finite bandwidth correction, as expected from an M6 earthquake in Indian peninsula shield for accurately predicting the level of ground motion. Estimates of shear wave velocity averaged over the top 30 m of soil (V_S30) are obtained from multichannel analysis of surface wave (MASW) at 210 sites at depths of 30 to 60 m below the ground surface. Using these V_S30 values, along with the available geotechnical information and synthetic ground motion database obtained, equivalent linear one-dimensional site response analysis that approximates the nonlinear soil behavior within the linear analysis framework was performed using the computer program SHAKE2000. Fundamental natural frequency, Peak Ground Acceleration (PGA) at surface and rock levels, response spectrum at surface level for different damping coefficients, and amplification factors are presented at different sites of the city. Liquefaction study was done based on the V_S30 and PGA values obtained. The major findings suggest show that the northeast part of the city is characterized by (i) low V_S30 values (<?200 m/s) associated with alluvial deposits, (ii) relatively high PGA value, at the surface, of about 0.24 g, and (iii) factor of safety and liquefaction below unity at three sites (no. 12, no. 37, and no. 70). Thus, this part of the city is expected to experience damage for the expected M6 target event.     

基于等损伤的主余震序列型地震动PGA放大系数谱研究

大地震发生后往往会伴随强余震发生,强烈余震会加重结构的破坏程度,在抗震设计及结构损伤评估中考虑余震的影响是一重要课题。采用NGA-West2数据库建立主余震序列型地震动记录数据集,基于损伤等效思想定义主余震作用下PGA放大系数α,针对单自由度体系结构,通过动力时程分析建立不同相对强度、不同场地条件下平均PGA放大系数谱,进一步通过回归分析构建PGA放大系数谱的预测方程,分析统计结果的离散性。结果表明：PGA放大系数谱受场地条件影响较小,受主余震相对强度影响显著;放大系数谱值随周期增大而减小。谱预测方程能够提供目标损伤下结构的主震PGA放大需求,并可以作为设计谱调整系数使用,以实现在结构抗震设计及结构损伤评估中考虑余震的影响。     

中强地震活动区地震动衰减关系的确定

地震动衰减关系是影响地震安全性评价特别是地震区划结果的重要因素.我国现行的地震动衰减关系主要是依据6级以上地震的地面运动资料得到的,并没有考虑中强地震的衰减特性.为此,文中利用现有的烈度资料和其他可供参考的研究成果来建立我国中强地震活动区的地震动衰减关系.收集了我国华中、华南、东北等地区的51次地震的烈度等震线资料,运用单随机变量加权最小二乘回归法得到中强地震活动区烈度衰减关系.然后以美国西部地区为参考地区,运用缺乏地震动参数的地震动估计方法-地震对映射法得到中强地震活动区峰值加速度和有效峰值加速度衰减关系.最后,通过与我国强地震区和中强地震区已有的烈度衰减关系和地震动衰减关系的对比,验证了得出的我国中强地震活动区烈度衰减关系和地震动衰减关系的合理性.     

Ground-Motion Hazard Values for Northern Algeria

This study examines distinctive features of ground motion parameters in northern Algeria. An initial computation of seismic hazard in terms of horizontal peak ground acceleration (PGA) and spectral acceleration (SA) at different periods, damped at 5%, is carried out for three different types of soils (rock, stiff soils and soft soils) for return periods of 100 and 475 years. In addition, uniform hazard spectra (UHS) are computed for these two return periods at several locations in the region. Then, the UHS computed for different soil types are proposed as a starting point to define elastic design spectra for building-code purposes. We have used the well-known Newmark-Hall approach. As proposed in the most recent International Building Codes, the SA (0.2 s) value is used to establish the spectral region for lower periods (region controlled by acceleration), whereas the SA (1.0 s) value is used to establish the spectral region for intermediate periods (region controlled by velocity). We also obtained important relations, dependent on site condition, between SA (0.2 s), SA (1.0 s) or SA_max values, and the PGA, for both return periods of 100 and 475 years. Other relationships between PGA or SA_max values have also been derived for return periods of 100 and 475 years, in this case independent of site condition.     

汶川地震诱发滑坡的地震动加速度评判标准

通过对汶川MS8.0地震强震记录资料的计算分析,给出一个与地震诱发滑坡关系比较密切的地震动加速度参数--最大临界加速度.分析表明,最大临界加速度与地震滑坡的相关性高于峰值加速度,而且可以消除峰值加速度作为判据使用时存在的问题.利用汶川地震强震资料建立了该参数在震区的衰减关系,由衰减关系得到的地震滑坡影响区域与汶川地震滑...     

2020年7月12日唐山古冶MS5.1地震的强震动特征

首都圈强震动观测台相对密集,区域强震动观测能力较强。2020年7月12日6时38分在河北省唐山市古冶区发生M_S5.1地震,强震动台和烈度台共收集到三分量强震动记录363组。通过强震动记录数据处理绘制了PGA和PGV的空间分布图,初步认识到此次地震强烈的地震动主要集中于震中附近区域,远离震中的天津塘沽—宁河—宝坻三角区域也观测到“异常增大趋势”的地震动,这与1976年唐山地震的宏观烈度异常区较为一致。同时分析了强震动记录PGA、PGV随距离的衰减规律,通过与为第五代区划图编制建立的东部强震区衰减关系的计算值对比,可以看出与此次地震中PGA、PGV的衰减具有一致的趋势。      

An explanation of the large PGA value of 2013 M S7.0 Lushan earthquake at 51BXD station through topographic analysis

In M _S7.0 Lushan earthquake, a large amount of strong ground motion recordings were collected. In this paper, we analyze the recordings carefully. The abnormality of ground motion recordings is identified through a log linear regression. In the station of 51BXD, the PGA value has exceeded 1 g, which is the biggest peak ground acceleration (PGA) value obtained from all recordings in this earthquake. The log linear relation shows the PGA value in this station is abnormally large. As this station is located on the footage of a hill, the topographic amplification factor is explored in order to explain this abnormality. Through 3D numerical modeling using spectral element method with transmitting boundary conditions, the amplification factor is quantized. In this station, the topographic amplification is highly polarized in the direction of East–West which agrees with the empirical recordings. This research result suggests us in future directionality of topographic amplification should be considered in the aseismic design.     

面向地震烈度仪的现地地震动预测模型的构建与验证

利用初期P波预警参数构建现地地震动预测模型,使其在达到设定阈值时快速发出报警信息,是现地地震预警系统面临的一个关键问题,直接关系到发布信息的准确性和及时性.针对地震烈度仪基于微机电系统传感器记录到的数据质量较差,通过两次积分获取的位移存在较大偏差,会引起更多的误报和漏报,本文采用不同阶数(1—4阶)的巴特沃斯滤波器,分...