Impact of Higher Frequency Content of Input Motion Upon Equivalent Linear Site Response Analysis for the Study Area of Delhi

While travelling through the subsoil layers, earthquake generated bedrock motions get modified significantly due to local soil and should be quantified using ground response analysis. Present study concentrates on equivalent linear method of site response analysis in SHAKE2000 software. It is a frequency based analysis tool having default frequency set to 15 Hz. While due consideration is given to amplitude, no to very limited information about the frequency content of the input motion to be considered in ground response analysis is available. In the present work, the effect of the maximum frequency of ground motion in site response analysis using SHAKE2000 is examined. Two sets of analyses are carried out in this work based on 30 globally recorded input motions. In the first analyses, input motion up to 15 Hz maximum frequency, which is a default value in SHAKE2000 is considered while second analyses are based on considering each of the 30 input motions up to the Nyquist frequency. Comparing the results from the two sets of analyses highlight that selection of maximum frequency in SHAKE2000 has considerable effect in ground motion amplification at different depths. As a result, even the peak ground acceleration which controls the building behavior and damage scenario, is going to change considerably even in case same input motion is used in the analysis.     

Evaluating the effects of ground motion parameters on response spectra in Uttarakhand Himalayas,India

Uttarakhand, a state of India, is located in seismically active Himalayan region and in the proximity of plate boundaries. The effects of important ground motion parameters like magnitude, distance, and local geology and site conditions on acceleration response spectra are examined in Uttarakhand Himalayas in this work. A total of 447 strong ground motion histories (horizontal and vertical) from 42 earthquakes were selected. The results show that the shape of the acceleration response spectra is influenced by the local site conditions and regional geology. The studies are carried out for two categories of sites, i.e., rock sites and soft soil sites. The maximum average horizontal spectral amplification for rock sites is 2.7 at 0.1 s, while for soft soil sites, it is found to be 3.2 at 0.2 s. In the same way, the maximum average vertical spectral amplification for rock is found to be 2.7 at 0.1 s, while for soft soil, it is found to be 2.95 at 0.1 s. The average spectral amplification in vertical component also shifts from low period (rock) to high period (soft soil). The level of spectra increases with decrease in distance for rock sites as well as soft soil sites. When comparing different magnitude earthquakes in different geological conditions, the response spectra are found to follow each other up to 0.04 s, while for period greater than 0.04 s, the spectra of higher magnitude earthquake is observed on the higher side. For soft soil sites, spectra from different magnitude earthquakes are observed to follow each other up to 0.1 s, beyond which they get separated.     

Comparative study on the equivalent linear and the fully nonlinear site response analysis approaches

Seismic site effect has been a major issue in the field of earthquake engineering due to the large local amplification of the seismic motion. This paper presents the importance of an appropriate soil behavior model to simulate earthquake site response and gives an overview of the field of site response analysis. Some of the well-known site response analysis methods are discussed. The objective of this paper is to investigate the influences of nonlinearity on the site response analysis by means of a more precise numerical model. In this respect, site responses of four different types of one-layered soil deposit, based on various shear wave velocities with the assumption of linear and rigid base bedrock, were analyzed by using the equivalent linear and fully nonlinear approaches. Nonlinear analyses?? results were compared with those of the linear method, and both of the similarities and differences are discussed. It is concluded that in the case of nonlinearity of soil under strong ground motions, 1-D equivalent linear modeling overestimates the amplification patterns in terms of absolute amplification level, and cannot correctly account for resonant frequencies and hysteric soil behavior. Therefore, more practical and appropriate numerical techniques for ground response analysis should be surveyed.     

Site-specific earthquake response study for hazard assessment in Kolkata city,India

The assessment of local site effects on seismic ground motions is of great importance in earthquake engineering practice. Several destructive earthquakes in the past have demonstrated that the amplification of ground motion and associated damage to structures due to local site conditions is a significant consideration in earthquake hazard analysis. A recent paper published in this journal highlights the hazard posed by earthquakes in the megacity of Kolkata in India due to its seismic and geological settings. The seismic hazard assessment study speculates that the deep alluvial deposit in the city may increase the seismic hazard probably due to the amplification of the seismic energies. This paper focuses on the seismic response studies of the various soil strata (i.e. for local subsurface conditions) obtained from various construction sites in the city for predicted earthquake. It is very well recognized that site response studies (a part of seismic microhazard zonation for urban areas) are the first step towards performance-based foundation design or seismic risk analysis and mitigation strategy. One of the problems for carrying out site-specific study in Kolkata is the lack of recorded strong motion data in the city. Hence, this paper outlines a methodology to carry out site-specific study, where no strong motion data or seismic data are available. The methodology uses wavelet-based spectrum compatibility approach to generate synthetic earthquake motions and equivalent linear method for seismic site response analysis. The Mega City of Kolkata has been considered to explain the methodology. Seismic hazard zonation map by the Bureau of Indian Standards classifies the City of Kolkata as moderate seismic zone (Zone III) with a zone factor 0.16. On the other hand, GSHAP(Global Seismic Hazard Assessment Program) map which is based on 10% probability of exceedance in 50 years specifies a maximum peak ground acceleration (PGA) of 1.6 m/s² (0.163 g) for this region. In the present study, the seismic response has been carried out based on GSHAP. The results of the analysis indicate the amplification of ground motion in the range of 4.46–4.82 with the fundamental period ranging from 0.81 to 1.17 s. Furthermore, the maximum spectral accelerations vary in the range of 0.78–0.95 g.     

Equivalent-Linear Seismic Ground Response Analysis of Some Typical Sites in Mumbai

Any earthquake event is associated with a rupture mechanism at the source, propagation of seismic waves through underlying rock and finally these waves travel through the soil layers to the particular site of interest. The bedrock motion is significantly modified at the ground surface due to the presence of local soil layers above the bedrock beneath the site of interest. The estimation of the amplifications in ground response due to the local soil sites is a complex problem to the designers and the problem is more important for mega cities like Mumbai in India, where huge population may get affected due to devastations of earthquake. In the present study, the effect of local soil sites in modifying ground response is studied by performing one dimensional equivalent-linear ground response analysis for some of the typical Mumbai soil sites. Field borelog data of some typical sites in Mumbai city viz. Mangalwadi site, Walkeswar site, BJ Marg near Pandhari Chawl site are considered in this study. The ground responses are observed for range of input motions and the results are presented in terms of surface acceleration time history, ratio of shear stress to vertical effective stress versus time, acceleration response spectrum, Fourier amplitude ratio versus frequency etc. The typical amplifications of ground accelerations considering four strong ground motions with wide variation of low to high MHA, frequency contents and durations are obtained. Results show that MHA, bracketed duration, frequency content have significant effects on the amplification of seismic accelerations for typical 2001 Bhuj motion. The peak ground acceleration amplification factors are found to be about 2.50 for Mangalwadi site, 2.60 for Walkeswar site and 3.45 for BJ Marg site using 2001 Bhuj input motion. The response spectrum along various soil layers are obtained which will be useful for designers for earthquake resistant design of geotechnical structures in Mumbai for similar sites in the absence of site specific data.     

Evaluation of ground motion characteristics,effects of local geology and liquefaction susceptibility: the case of Itea,Corinth Gulf (Greece)

In the present paper we analyze the effect of local geology on ground motion by means of numerical calculations (numerical models) using total (TS) and effective stress (ES) methods. These numerical calculations have been applied to the site of Itea, Corinth Gulf, which was chosen based on liquefaction susceptibility criteria and field inspection. Data from seismic refraction experiments and cone penetration test N-values as well as selected records of ground motion in nearby areas were used to construct the input file for the numerical model. By means of␣dynamic analysis such characteristics of ground motion as acceleration time histories, response spectra, and amplification function were evaluated. A one-dimensional soil amplification effect was clearly shown. Liquefaction probability at the Itea site was predicted based on the safety factor and the calculation of the induced settlement at the test site. Results of the TS and ES modeling lead us to conclude: (1) the presence of soft soil at Itea caused significant amplification (almost 2.5-fold higher magnitude) of the underlying bedrock motion and, therefore, can contribute to damage; (2) the area of Itea is highly susceptible to liquefaction due to presence of silty sand deposits at depths between 2.48 m (the position of the water table) and 12 m that demonstrate the rapid growth of the excess pore water pressure (EPWP) ratio with an increase in peak ground acceleration values.     

远场大地震作用下大尺度深软场地的非线性地震效应分析

基于ABAQUS软件自行研制的并行计算显式算法集群平台,针对苏州城区典型地层剖面,建立了大尺度深软场地的二维精细化非线性有限元分析模型,对人工地震波和大地震远场地震动作用下深软场地的非线性地震效应进行了比较研究。研究结果表明：（1）与人工地震波作用时深软场地的地表峰值加速度放大效应相比,大地震远场地震波作用时的放大效应尤为显著,由于土介质的横向不均匀性及其非线性,使不同地表的峰值加速度放大效应存在明显的变异性。（2）深软场地对周期小于0.3 s的高频地震波均具有显著的滤波效应;大地震远场地震波作用时,深软场地对周期0.85～1.65 s的长周期地震波的放大效应非常显著,但对2.5～7.0 s的长周期地震波呈现出明显的滤波效应。（3）地震动峰值加速度PGA值沿土层深度和横向的分布形态呈现出明显的高低起伏现象,在不同成因的土层更迭面附近及土介质横向不均匀性显著的区域,地震波的局部聚焦放大和过滤减小现象尤为明显,且大地震远场地震动作用时,20 m以浅土层的PGA值呈现出非常显著的放大效应。（4）地震波的频谱特性、土层的横向不均匀性对深软场地地表加速度反应谱? 谱的谱形有显著影响;给出了描述加速度反应谱沿土层深度变化特征的三维谱形曲线,可以直观地展示出深软场地中细长地下结构地震反应可能存在类共振现象的土层深度。     

Dynamic analysis of ground with rigorous use of strain dependency and its application to seismic microzonation of alluvial plane

Seismic microzonation is one of the most important measures to mitigate earthquake hazards in urban areas. Because the ground motion varies significantly with the subsurface geology, it is needed for microzonation to account as much as possible for the local soil conditions. Noteworthy is that nonlinear deformation properties of soil play essential roles in amplification of strong ground motion. It is desired furthermore to focus on the expected damage extent in addition to the calculated maximum acceleration and/or velocity. The present study first developed a computer code for one-dimensional response analysis of ground that reasonably takes into account nonlinear dynamic soil properties. Second, correlations between the calculated ground motion and damage extent were obtained by examining seismic damages during the past earthquakes. By combining these two issues, seismic microzonation was carried out, and detailed damage distribution was assessed. The product of this study covers not only the damage caused by ground shaking but also liquefaction problem and lifeline damage.     

Quantifying sensitivity of local site response models to statistical variations in soil properties

We perform a combined stochastic-deterministic analysis of local site response using two computer codes, an equivalent linear analysis program SHAKE and a fully nonlinear finite element code SPECTRA. Our goal is to compare the relative sensitivity of the two codes to statistical variations in soil properties. For the case studies, we re-analyze two ground motion records in Lotung, Taiwan, and one ground motion record in Gilroy, California, utilizing the recorded ground motions at the site deterministically as input into the two codes while treating the uncertain soil parameters as random variables. We then obtain empirical cumulative distribution functions of the Arias intensity and acceleration spectrum intensity, two measures of cumulative damage, to compare the relative sensitivity of the two codes to variations in model parameters. We show that the two codes exhibit comparable sensitivities to statistical parameter variations, indicating that even in the presence of fluctuations in the soil parameter values it is possible to pursue a fully nonlinear site response analysis with SPECTRA and benefit from its superior accuracy.     

汶川地震远场黄土场地地震动场地放大效应机制研究

通过汶川地震震后科学考察发现,汶川地震局部场地的震害和地震动放大效应明显,地震对甘肃省境内远离震中的黄土地区造成了较为严重的破坏。在此基础上,结合钻孔波速测试,运用二维等价线性时程响应动分析法对甘肃省平凉市典型黄土塬进行了地震动力响应计算,分析了该地区局部场地条件对地震动放大效应的影响。研究发现,黄土塬具有地震动放大效应,随着塬高的增加,加速度、速度、位移均出现放大效应;随着地震动在黄土覆盖层中的传递,卓越频率向低频移动,高频成分被吸收,加速度反应谱中长周期分量逐渐增大。数值计算结果与实际震害基本吻合,其结果对黄土地区进行合理的抗震设防具有一定的科学意义。     

Seismic hazard assessment of Chennai city considering local site effects

Chennai city suffered moderate tremors during the 2001 Bhuj and Pondicherry earthquakes and the 2004 Sumatra earthquake. After the Bhuj earthquake, Indian Standard IS: 1893 was revised and Chennai city was upgraded from zone II to zone III which leads to a substantial increase of the design ground motion parameters. Therefore, a comprehensive study is carried out to assess the seismic hazard of Chennai city based on a deterministic approach. The seismicity and seismotectonic details within a 100 km radius of the study area have been considered. The one-dimensional ground response analysis was carried out for 38 representative sites by the equivalent linear method using the SHAKE91 program to estimate the ground motion parameters considering the local site effects. The shear wave velocity profile was inferred from the corrected blow counts and it was verified with the Multichannel Analysis of Surface Wave (MASW) test performed for a representative site. The seismic hazard is represented in terms of characteristic site period and Spectral Acceleration Ratio (SAR) contours for the entire city. It is found that structures with low natural period undergo significant amplification mostly in the central and southern parts of Chennai city due to the presence of deep soil sites with clayey or sandy deposits and the remaining parts undergo marginal amplification.     

刚度和阻尼频率相关的等效线性化方法

在实际工程中发现,在强地震动作用下等效线性化方法低估了地表地震动的峰值加速度反应。通过与场地效应竖向台阵观测记录的对比分析证实了这一现象,其主要原因是等效剪应变水平的取值高于高频段实际剪应变水平,从而抑制了土层地震反应中的高频成份。若干分析表明,动力反应过程中土体的剪应变是随频率而变化的,且试验表明土体的刚度和阻尼是与剪应变相关的,据此假设土体刚度和阻尼是与频率相关的。基于这一假设,对传统等效剪应变的选取如何使得土层对基底地震动输入中的高频成分产生滤波效应的原理进行探讨,且通过频率相关的等效线性化方法对该原理进行证明。将该方法应用于Port Island场地效应竖向台阵的分析,并与实际地震动记录进行了对比,结果表明,较之传统等效线性化方法,能更加合理地考虑地震动输入中的高频成分。     

A comparative study on linear and nonlinear site response analysis

Site response analysis is usually the first step of any seismic soil-structure study. Geotechnical earthquake engineers and engineering geologist have been trying to find both practical and most appropriate solution techniques for ground response analysis under earthquake loadings. The paper attempts to give a critical overview of the field of site response analysis. In this paper, the influences of nonlinearity on the site response analysis summarized and were evaluated with a numerical example. Site response of a two layered soil deposit with the assumption of linear and rigid base bedrock (or viscoelastic half-space) was analyzed by using linear and nonlinear approaches. The amplification spectrum of the soil column is computed between the top and the bottom of this soil deposit. Nonlinear analysis was compared with the linear method of analysis. Steps involved in ground response analyses to develop site-specific response spectra at a soil site are briefly summarized. Some of the well-known site response analysis methods are summarized and similarities and differences between linear and nonlinear methods are compared by a numerical example.     

场地土对基岩峰值加速度放大效应分析

通过实际土层地震反应结果的统计分析和强震加速度观测结果的对比, 讨论了不同场地条件对基岩峰值加速度的放大效应及其特点。该分析可为地震动参数区划图编制和地震安全性评价中场地效应的估计、由基岩地震动估算场地地面地震动提供参考。     

Nonlinear soil response to ground motions during different earthquakes in Nepal,to arrive at surface response spectra

Catastrophic damages reported during an earthquake include building damages, excessive ground shaking, uneven settlements and liquefaction. While most of the seismic hazard studies map the probable level of ground shaking at the bedrock level, their use in assessing the above damages is very limited until the response of the local soil is also taken into account. Determination of the local soil response needs regionally recorded ground motions, dynamic soil properties, in situ geotechnical details, etc., which most of the time are not readily available for the region under study. In the present work, the response of local soil for Nepal has been studied indirectly taking into account the surface level of ground shaking during various past as well as recent EQs observed at various locations. Based on the present analysis, a low value of amplification factor for high peak horizontal acceleration and vice versa is observed in central, western as well as southern parts of Nepal. These observations suggest nonlinear soil behavior and are in accordance with the available literature. Further, the ground motion records during 2015 Nepal EQ show maximum soil response at 0.3 s which is exactly matching with the site class C obtained from in situ data for the above locations. Based on the above observations, various correlations between the high peak horizontal acceleration and the surface spectral acceleration are proposed to obtained site specific surface response spectrum for Nepal.     

Estimation of Response Spectra in the Severely Damaged Area During the 1995 Hyogo-Ken Nambu Earthquake

We calculated theoretical acceleration response spectra at the ground surface in the region near the 1995 Hyogo-ken Nambu earthquake source based on the spectral moment method. To estimate earthquake motion on the ground surface, a formula of earthquake motion at base rock level was derived. The amplification effect of the ground was introduced by using multiple reflection theory. Theoretically estimated response spectrum were modified by the response spectra calculated using observed earthquake motions.     

Ground motion estimation at Guwahati city for an Mw 8.1 earthquake in the Shillong plateau

In this paper, the ground motion at Guwahati city for an 8.1 magnitude earthquake on Oldham fault in the Shillong plateau has been estimated by stochastic finite-fault simulation method. The corresponding acceleration time histories on rock level at several sites in the epicentral region have been computed. These results are validated by comparing them with the estimates obtained from Medvedev–Sponheuer–Karnik (MSK) intensity observations of 1897 Shillong earthquake. Using the local soil parameters, the simulated rock level acceleration time history at Guwahati city is further amplified up to the ground surface by nonlinear site response analysis. The results obtained are presented in the form of peak ground acceleration (PGA) contour map. The maximum amplification for PGA over Guwahati city is as high as 2.5. Based on the simulated PGA, the liquefaction susceptibility at several locations in the city has been estimated. The results are presented in the form of contours of factor of safety against liquefaction at different depths below the ground surface. It is observed that over a large part of the Guwahati city, the factor of safety against liquefaction is less than one, indicating that the city is highly vulnerable to liquefaction in the event of this earthquake. The contour maps obtained can be used in identifying vulnerable areas and disaster mitigation.     

脉冲地震动中竖向加速度对场地液化的影响

为研究近断层脉冲地震动中竖向加速度对砂土场地液化的影响,基于有限元平台OpenSees开发的边界面塑性本构模型,建立了动单剪单元试验模型和饱和砂土三维有限元模型。选取台湾Chi-Chi地震中10条具有速度脉冲特性的地震波,对比分析了水平双向脉冲波与三向脉冲波作用下土柱竖向位移、循环应力比、孔压比及等效循环周数的差异性,继而明确了脉冲地震动中竖向加速度对砂土液化的影响规律。研究表明,三向脉冲地震波中竖向加速度分量对场地永久位移值影响较小,但使永久位移的发展持时明显增大;土柱循环应力比受竖向地震动影响较小,因此分析脉冲地震动对场地剪切特性的影响时,可将三向脉冲地震动简化为水平双向地震动;考虑竖向地震动的三向脉冲地震波引起的孔压比变化幅度较大,孔压消散时间较长;三向脉冲地震波对应的等效循环周数较大,地震动发展持时长,可认为竖向加速度对场地液化有促进作用。     

Equivalent-linear site response analysis on the site of the historical Trakošćan Castle,Croatia, using HVSR method

This paper presents a major extension of seismic vulnerability research project on the site of Trako??an Castle based on the initial horizontal-to-vertical-spectral-ratio (HVSR) results from Stanko et al. (2016). The estimated HVSR site frequencies and HV amplification at Trako??an Castle can only be used as an indication of the initial soil site frequency and amplification, so-called natural soil model, corresponding to the subsoil profile without the influence of an earthquake. The equivalent-linear (EQL) site response analysis has been carried out for different earthquake scenarios for a maximum input rock peak ground acceleration (PGA_ROCK) that corresponds to return periods of 95 (0.08 g), 475 (0.18 g) and 1000 years (0.31 g). The aim of the research is to evaluate structural seismic design responses and to determine type and degree of damage caused by local site effect, which is the result of an alluvial basin and topographic influences. The main objective of this research is the formation of local microseismic zones based on an EQL analysis: surface spectral acceleration and amplification maps at the predominant frequency. Based on the HVSR frequency response of the core structure of Trako??an Castle and the Tower itself (fundamental and higher frequency modes), maps of surface spectral acceleration and soil amplification at different frequencies (3, 5 and 10 Hz) are developed for different input PGA_ROCK levels (0.08, 0.18 and 0.31 g) to evaluate seismic response of the Castle. Observed amplifications are correlated with ground motion polarization and directionality of the ground motion from the alluvial basin to the hilltop. Shortening of predominant frequencies (lengthening of the period), particularly in the alluvial basin, has been observed with higher input PGA_ROCK in the EQL analysis. This effect is not manifested in the Trako??an hill, and predominant frequencies match HVSR frequencies. The use of certain geophysical survey methods at historical sites is a big problem, because terrain features (e.g. steep hills, mountains, ridges, slopes, cliffs) create lack of space and make it impossible to carry out geophysical investigation. Microtremor measurements at historical sites can overcome this limitation and provide local seismic response and vulnerability behaviour of historical monuments without destroying their authenticity. Also, computational modelling can greatly improve the results. The EQL site response analysis on the site of Trako??an Castle has confirmed and improved the results of seismic response and vulnerability based on HVSR method.     

黄土塬地震动响应特征分析

2008年汶川Ms8.0级强烈地震对远离震中的黄土塬地区造成了较为严重的破坏,局部场地的震害和地震动放大效应显著。以典型黄土塬场地为对象,应用地形台阵流动观测和有限元计算方法,系统开展强震动作用下黄土塬边坡场地动力响应特征研究;应用大型振动台模拟试验和数值计算方法,重点探讨黄土塬平台场地在汶川地震作用下地表加速度响应随覆盖层厚度、地震动强度的变化规律以及对建筑结构的潜在影响。结果表明：黄土塬边坡顶部存在低频但较高放大系数的现象,可能与斜坡高差与入射波波长之比密切相关;大角度（60°～70°）黄土塬边坡对地震动的放大效应十分显著,坡顶加速度峰值（PGA）放大可达2倍,反应谱卓越周期放大可达5倍;较厚黄土塬平台场地加速度放大可达2倍以上,地震烈度增加1度,对于场地上固有周期0.7～2.0 s和周期大于3.0 s的建筑物地震反应将显著增加。