Seismic characterization and reconstruction of reference ground motion at accelerometric sites of the Italian national accelerometric network (RAN)

We present a field procedure that has been extensively used in Italy to characterize local seismic response at accelerometric sites and to retrieve ground motion at reference soil conditions by deconvolution analysis. To allow a generalized application to large areas where borehole data are generally lacking or inadequate for the seismic characterization for soils down to the reference seismic bedrock, cost-effectiveness of the considered procedures is a main issue. Thus, major efforts have been devoted to optimize available information and exploit fast and cheap surface geophysical prospecting. In particular, geological/geomorphological survey and passive seismic prospecting (both in single- and multi-station configurations) were jointly considered to reconstruct seismo-stratigraphical site conditions. This information was then used to feed numerical modeling aiming at computing the local seismic response and performing a deconvolution analysis to reconstruct ground motion at reference soil conditions. Major attention was devoted to evaluate and manage uncertainty involved in the procedure and to quantify its effect on final outcomes. An application of this procedure to a set of sites included in the Italian Accelerometric Network is presented.

     

Geo-Characterization at selected accelerometric stations in Crete (Greece) and comparison of earthquake data recordings with EC8 elastic spectra

To estimate the seismic response according to Eurocode (EC8) and almost all other national codes, site conditions have to be properly characterized so that soil amplification and the corresponding peak ground motion can be calculated. In this work, different geophysical and geotechnical methods are combined in order to define the detailed ground conditions in selected sites of the Hellenic Accelerometric Network (HAN) in Crete. For this purpose, the geological information of the sites and shear wave velocity, calculated from surface wave measurements, is used. Additionally, ground acceleration data recorded through HAN have been utilized from intermediate depth earthquakes in the broader area of South Aegean Sea. Using the recorded ground motion data and the procedure defined in EC8, the corresponding elastic response spectrum is calculated for the selected sites. The resulting information is compared to the values defined in the corresponding EC8 spectrum for the seismic zone that includes the island of Crete. The comparison shows that accurate definition of ground type through geological, geotechnical and geophysical investigations is important. However, our current comparison focuses on the distribution of values rather than the absolute values of EC8-prescribed spectra, and the results should be considered in this context.     

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.     

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.     

Analysis of soil effects and distribution of damage from the Pyrgos 1993 (Greece) earthquake

Summary On 26 March, 1993, a moderate magnitude earthquake (M _s=5.5) occurred at 3 km epicentral distance from the town of Pyrgos, in Southern Greece, causing extensive damage to masonry houses. To explain the variability of seismic intensity over the town and to propose measures against future seismic activity, a microzonation study was undertaken which combined geological, geophysical and geotechnical investigations, site specific analyses of seismic ground response and detailed recording of structural damage. The analytical predictions of ground response are correlated to soil conditions and then used to identify (micro-)zones of sites with similar seismic response. Furthermore, they are compared to quantitative estimates of damage distribution over the town. It is concluded that the peak ground acceleration, normalized against the input peak seismic acceleration, is a function of the local soil conditions as well as the seismic excitation characteristics. Hence, it cannot be defined uniquely at a site, without reference to the seismic excitation. However, the normalized peak ground velocity and the acceleration response spectra are mainly functions of the soil conditions and can be used as criteria for the practical definition of (micro-)zones. The distribution of damage in various parts of the town is at least partially attributable to local soil effects. The small epicentral distance of the earthquake, connected with the direction of the fault rupture, as well as the quality and techniques of construction, are additional factors that may have influenced the extent and distribution of damage.     

Site-Specific Analysis of Strong Motion Data from the September 7, 1999 Athens,Greece Earthquake

The strong ground motion from Athens, Greece 07/09/1999 earthquake has been recorded by eighteen (18) stations, fourteen (14) within the central Athens area and four (4) at the centers of nearby towns. The ground conditions for most of the recording sites were identified, based on previous geotechnical investigations carried out in the wider area of the sites, and consequently correlated to the seismic motion characteristics. Hence, it has been possible to evaluate the accuracy of different seismological methods for site characterization and also estimate soil effects on peak ground acceleration and elastic response spectra. In addition, preliminary estimates are drawn for the seismic motion characteristics at the epicentral area, where no strong motion recordings are available. The detailed soil profiles at the recordingsites are placed in the Appendix.     

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

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

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.     

Effects of seismic wave propagation on a long and narrow building: Body wave and surface wave propagation

In the case of the dynamic analysis of the structures using the recorded earthquake ground motions, it is usually assumed that the ground motion consists of body waves propagating vertically. However, the response of a long and narrow structure may be influenced by the oblique propagation of body waves and the dispersion of surface waves. In this paper, the effects of the seismic wave propagation on the response of this kind of structure are investigated. The characteristics of the wave propagation were verified using the recorded motions and soil information at the building site. The ground motion at every input point of the building was evaluated using the difference of arrival time of seismic waves calculated by assuming the velocity and the direction of the wave propagation. Using these ground motions, response analyses of the lumped mass model for the structure were performed. By considering the characteristics of the seismic wave propagation, the average response decreased but the local response increased around the end of the roof. Further studies of the structure were also performed in order to restrain the response around the end of the roof.     

Estimation of seismic spectral acceleration in Peninsular India

Peninsular India (PI), which lies south of 24°N latitude, has experienced several devastating earthquakes in the past. However, very few strong motion records are available for developing attenuation relations for ground acceleration, required by engineers to arrive at rational design response spectra for construction sites and cities in PI. Based on a well-known seismological model, the present paper statistically simulates ground motion in PI to arrive at an empirical relation for estimating 5% damped response spectra, as a function of magnitude and source to site distance, covering bedrock and soil conditions. The standard error in the proposed relationship is reported as a function of the frequency, for further use of the results in probabilistic seismic hazard analysis.     

Seismic hazard and site-specific ground motion for typical ports of Gujarat

Economic importance of major ports is well known, and if ports are located in seismically active regions, then site-specific seismic hazard studies are essential to mitigate the seismic risk of the ports. Seismic design of port sites and related structures can be accomplished in three steps that include assessment of regional seismicity, geotechnical hazards, and soil structure interaction analysis. In the present study, site-specific probabilistic seismic hazard analysis is performed to identify the seismic hazard associated with four typical port sites of Gujarat state (bounded by 20°–25.5°N and 68°–75°E) of India viz. Kandla, Mundra, Hazira, and Dahej ports. The primary aim of the study is to develop consistent seismic ground motion for the structures within the four port sites for different three levels of ground shaking, i.e., operating level earthquake (72 years return period), contingency level earthquake (CLE) (475 year return period), and maximum considered earthquake (2,475 year return period). The geotechnical characterization for each port site is carried out using available geotechnical data. Shear wave velocities of the soil profile are estimated from SPT blow counts using various empirical formulae. Seismicity of the Gujarat region is modeled through delineating the 40 fault sources based on the seismotectonic setting. The Gujarat state is divided into three regions, i.e., Kachchh, Saurashtra, and Mainland Gujarat, and regional recurrence relations are assigned in the form of Gutenberg-Richter parameters in order to calculate seismic hazard associated with each port site. The horizontal component of ground acceleration for three levels of ground shaking is estimated by using different ground motion attenuation relations (GMAR) including one country-specific GMAR for Peninsular India. Uncertainty in seismic hazard computations is handled by using logic tree approach to develop uniform hazard spectra for 5% damping which are consistent with the specified three levels of ground shaking. Using recorded acceleration time history of Bhuj 2001 earthquake as the input time motion, synthetic time histories are generated to match the developed designed response spectra to study site-specific responses of port sites during different levels of ground shaking. It is observed that the Mundra and Kandla port sites are most vulnerable sites for seismic hazard as estimated CLE ground motion is in order of 0.79 and 0.48 g for Mundra and Kandla port sites, respectively. Hazira and Dahej port sites have comparatively less hazard with estimated CLE ground motion of 0.17 and 0.11 g, respectively. The ground amplification factor is observed at all sites which ranges from 1.3 to 2.0 for the frequency range of 1.0–2.7 Hz. The obtained spectral accelerations for the three levels of ground motions and obtained transfer functions for each port sites are compared with provisions made in Indian seismic code IS:1893-Part 1 (2002). The outcome of present study is recommended for further performance-based design to evaluate the seismic response of the port structures with respect to various performance levels.     

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.     

大型地下洞室群地震动输入机制探讨

当前地下工程地震稳定性分析中对地震响应规律研究较多,而对分析中地震动输入机制研究则明显不足。本文针对这一现状,对大型地下洞室群地震响应分析中的地震动输入机制进行了一定的讨论。首先提出将分析中所涉及的地震动分为指定的实测强震记录、优选的实测强震记录、人工地震波3类,分别讨论了每种地震动的选取方法,形成了一套实用性较强的选取流程。随后提出了对设计地震动基准面(点)的建议,并研究了基于基准面(点)的人工边界输入地震动量值的确定方法,形成了适用于高山峡谷地形下大型地下洞室群地震动输入方法。将上述地震波的合理选取方法和高山峡谷地形下大型地下洞室群地震动的输入方法相组合,最终形成了一套可操作性较强的高山峡谷地形下大型地下洞室群地震动的输入机制。最后以白鹤滩水电工程地下洞室群为工程实例,对本文提出的地震动合理选取方法和高山峡谷地形下大型地下洞室群地震动的输入方法进行了工程应用,表明本文的方法为工程抗震设防分析提供了优质的地震输入参数。     

A study of the seismic response of the city of Benevento (Southern Italy) through a combined analysis of seismological and geological data

A previous analysis [Improta, L., G. Di Giulio, and A. Rovelli (2005). Variations of local seismic response in Benevento (Southern Italy) using earthquakes and ambient noise recordings, J. Seism. 9, 191–210.] of small magnitude earthquakes recorded at 12 sites within the city of Benevento has stressed the significant role played by near-surface geology in causing variability of the ground motion. In this paper, we extend the study of the seismic response from 12 sites to the entire urban area. Based on inferences from the comparison at the 12 sites between earthquake and ambient vibration results, we have collected ambient noise at about 100 sites within the city, intensifying measurements across the main shallow geological variations. We use borehole data to interpret ambient noise H/V spectral ratios in terms of near-surface geology comparing H/V curves to theoretical transfer functions of 1D models along five well-constrained profiles.

On the basis of geological, geotechnical, and seismic data, we identify three main typologies of seismic response in the city. Each type of response is associated to zones sharing common soil conditions and similar soil classes according to building codes for seismic design. Moreover, we find that the spatial variation of the seismic response in the ancient town area is consistent with the damage pattern produced by a very destructive, well-documented historical earthquake that struck the city in 1688, causing MCS intensity of IX–X in Benevento.

Finally, we use ground motions recorded during the experiment by Improta et al. [Improta, L., G. Di Giulio, and A. Rovelli (2005). Variations of local seismic response in Benevento (Southern Italy) using earthquakes and ambient noise recordings, J. Seism. 9, 191–210.] to generate synthetic seismograms of moderate to strong (Mw 5.7, Molise 2002 and Ms 6.9, 1980 Irpinia) earthquakes. We calibrate the random summation technique by Ordaz et al. [Ordaz, M., J. Arboleda, and S.K. Singh (1995). A scheme of random summation of an Empirical Green's Function to estimate ground motions for future large earthquakes, Bull. Seism. Soc. Am. 85, 1635–1647.] using recordings of these earthquakes available in Benevento. After a satisfactory fit between observed and synthetic seismograms, we compute response spectra at different sites and speculate on effects of the geology class at large level of shaking, including soil nonlinearity. We find that large discrepancies from design spectra prescribed by seismic codes can occur for a wide sector of Benevento, especially for periods < 0.5 s.     

大连开发区场地地震设计参数

根据工作区潜在震源区划分、潜在震源区地震活动性参数和地震动衰减关系, 进行场地地震危险性分析, 得到不同概率水平下场区相应地震烈度和基岩水平加速度峰值及其反应谱。根据场地工程地震条件划分不同地质单元及相应的场地类别, 进行不同概率水平的地震反应分析计算, 确定地震动设计参数。     

Nonlinear Analysis of Local Site Effects on Seismic Ground Response in the Bam Earthquake

The influence of local geologic and soil conditions on the intensity of ground shaking is addressed in this study. The amplification of the ground motion due to local site effects resulted in severe damage to dwellings in the Bam area during the 2003 Bam Earthquake. A unique set of strong motion acceleration recordings was obtained at the Bam accelerograph station. Although the highest peak ground acceleration recorded was the vertical component (nearly 1 g), the longitudinal component (fault-parallel motion) clearly had the largest maximum velocity as well as maximum ground displacement. Subsurface geotechnical and geophysical (down-hole) data in two different sites have been obtained and used to estimate the local site condition on earthquake ground motion in the area. The ground response analyses have been conducted considering the nonlinear behavior of the soil deposits using both equivalent linear and nonlinear approaches. The fully nonlinear method embodied in FLAC was used to evaluate the nonlinear soil properties on earthquake wave propagation through the soil layer, and compare with the response from the equivalent linear approach. It is shown that thick alluvium deposits amplified the ground motion and resulted in significant damage in residential buildings in the earthquake stricken region. The comparison of results indicated similar response spectra of the motions for both equivalent and nonlinear analyses, showing peaks in the period range of 0.3–1.5 s. However, the amplification levels of nonlinear analysis were less than the equivalent linear method especially in long periods. The observed response spectra are shown to be above the NEHRP building code design requirements, especially at high frequencies.     

An effective procedure for seismic hazard analysis including nonlinear soil response

This paper presents probabilistic seismic hazard analysis (PSHA) of Tehran, Iran, accounting the effect of nonlinear soil response. It is well-known that soil nonlinearity and its accurate prediction could play important role in seismic hazard study. For this purpose, two different approaches have been carried out for predicting the hazard curves by (1) applying site modifications to the ground motion prediction equation based on generic site classes and use of constant coefficients (2) using a close-form solution that modifies the hazard results at the rock level. Also, efficiency of the Monte Carlo method in modeling of amplification function for the six selected sites in the study area was examined. Results showed important effect of nonlinear soil response mainly for frequencies lower than 8?Hz, which should be considered properly in hazard estimation. As an interesting subject, influence of soil plasticity index (PI) on hazard estimation of clayey sites including the nonlinear soil response was evaluated.     

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.     

西安地裂缝场地地震效应分析

地裂缝是由内、外营力以及人类活动等因素的作用引致发生的一种地面破坏现象。在这种不良的地质条件下,必然会对场地的地震响应带来一定程度的影响。针对地裂缝场地地震效应研究和工程评价的不足,采用有限元法研究了地震荷载作用下,地裂缝对场地地震动效应的影响问题。分析得出了地裂缝对场地地表水平峰值加速度及反应谱的影响特征和影响范围,为地裂缝场地工程结构的抗震设防和地裂缝灾害避让提供了新的参考。     

Influence of seismic motions on behavior of piles in liquefied soils

In the present study an analytical procedure based on finite element technique is proposed to investigate the influence of vertical load on deflection and bending moment of a laterally loaded pile embedded in liquefiable soil, subjected to permanent ground displacement. The degradation of subgrade modulus due to soil liquefaction and effect of nonlinearity are also considered. A free headed vertical concrete elastic nonyielding pile with a floating tip subjected to vertical compressive loading, lateral load, and permanent ground displacement due to earthquake motions, in liquefiable soil underlain by nonliquefiable stratum, is considered. The input seismic motions, having varying range of ground motion parameters, considered here include 1989 Loma Gilroy, 1995 Kobe, 2001 Bhuj, and 2011 Sikkim motions. It is calculated that maximum bending moment occurred at the interface of liquefiable and nonliquefiable soil layers and when thickness of liquefiable soil layer is around 60% of total pile length. Maximum bending moment of 1210 kNm and pile head deflection of 110 cm is observed because of 1995 Kobe motion, while 2001 Bhuj and 2011 Sikkim motions amplify the pile head deflection by 14.2 and 14.4 times and bending moment approximately by 4 times, when compared to nonliquefiable soil. Further, the presence of inertial load at the pile head increases bending moment and deflection by approximately 52% when subjected to 1995 Kobe motion. Thus, it is necessary to have a proper assessment of both kinematic and inertial interactions due to free field seismic motions and vertical loads for evaluating pile response in liquefiable soil.