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.     

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 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.     

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.     

Seismic response of high concrete face rockfill dams subjected to non-uniform input motion

Analysis of the seismic response of high CFRDs under non-uniform ground motion input is conducted using a novel non-uniform input motion calculation method combined with nonlinear FEM. The non-uniform input motion calculation method and its basic assumption are validated. The response of CFRDs under uniform and non-uniform input is compared to discuss the necessity of conducting seismic analysis of high CFRDs under realistic non-uniform ground motion input. When the acceleration at the surface of the free field for dynamic simulations with uniform and non-uniform input is kept consistent, the seismic response of CFRDs under non-uniform input is in general significantly smaller, while the dynamic tensile stress around the edges of the concrete face slab is greater. The simulation results suggest that non-uniformity of the ground motion input has important effects on the seismic response of high CFRDs and should be considered in the seismic design of CFRDs. The influence of the incident angle of seismic waves is also investigated, with results indicating that the influence is waveform dependent, while being frequency independent.

     

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.     

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.     

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.     

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.     

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.     

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

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

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.     

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.     

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.     

Geologic site conditions and site coefficients for estimating earthquake ground motions in the inland areas of Korea

Site characterization and site-specific ground response analyses were conducted at two representative inland areas in Korea. In situ tests included 25 boring investigations, 7 crosshole tests, 18 downhole tests and 41 SASW tests, and in the laboratory, resonant column tests were performed. The soil deposits in Korea, which were shallower and stiffer than those in the western US, were examined. The fundamental site periods were distributed in the narrow band ranging from 0.1 to 0.4 s. Most sites were designated as site classes C and D based on the mean shear wave velocity of the upper 30 m from the current Korean seismic design guide. Based on the ratio of the acceleration response spectra of ground surface to rock-outcrop, short-period (0.1–0.5 s) site coefficient, F_a ranged from 1.0 to 2.7, and mid-period (0.4–2.0 s) site coefficient, F_v ranged from 1.0 to 1.6, regardless of the input rock outcrop acceleration levels of 0.05 and 0.14 g. The site coefficients specified in the Korean seismic design guide, which is similar to NEHRP provisions and UBC, underestimate the ground motion in the short-period band and overestimate the ground motion in the mid-period band. These differences can be explained by the differences in the depth to bedrock and the soil stiffness profile between Korea and western US. Also, the site coefficients should be re-evaluated accounting for the local geologic conditions on the Korean peninsula.     

Site response zones and short-period earthquake ground motion projections for the Las Vegas Basin

A deterministic seismic hazard analysis was conducted to address the effect of local soil conditions on earthquake-induced strong ground motion in the Las Vegas Basin, Nevada (US). Using a large geological and geotechnical database, two response units were defined: a fine-grained unit, predominantly clay; and a coarse-grained unit, predominantly gravel. A moderate number of high-quality shallow shear wave velocity measurements were collected from which characteristic shear wave velocity profiles were developed for each response unit. An equivalent-linear one-dimensional site response model was used. The model was calibrated using a basin-wide, small-strain ground motion database. Calibration tests showed that ground motion projections become increasingly conservative with increasing ground-motion amplitude. Projections were overconservative for the coarse-grained response unit, likely due to the sparseness of the velocity database. For the earthquake response analyses, historical ground motions were used to model characteristic ‘bedrock’ motion for earthquakes on 10 faults judged to be critical. Response spectral envelopes were generated for each unit through Monte-Carlo simulations. For the fine-grained response unit, 95th percentile peak ground acceleration, peak spectral acceleration and predominant period were 310 cm/s², 1100 cm/s², and 0.29 s, respectively. With respect to codified design spectra, projections are lower at short periods and higher at long periods. Projections of peak spectral accelerations for the coarse-grained response unit, were more than double that of codified spectra; however, they are believed to be overconservative. Near-fault effects and basin-edge effects, though potentially important, were not considered in these analyses.     

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.