SEISMIC AMPLIFICATION AT A SOFT SOIL SITE WITH LIQUEFIABLE LAYER

It is well known that local soil conditions play a key role in the amplification of earthquake waves. In particular, a liquefiable shallow soil layer may produce a significant influence on ground motion during strong earthquakes. In this paper, the response of a liquefiable site during the 1995 Kobe earthquake is studied using vertical array records, with particular attention on the effects of nonlinear soil behaviour and liquefaction on the ground motion. Variations of the characteristics of the recorded ground motions are analysed using the spectral ratio technique, and the nonlinearity occurring in the shallow liquefied layer during earthquake is identified. A fully coupled, inelastic finite element analysis of the response of the array site is performed. The calculated stress-strain histories of soils and excess pore water pressures at different depths are presented, and their relations to the characteristics of the ground motions are addressed.     

SITE EFFECTS BY H/V RATIO: COMPARISON OF TWO DIFFERENT PROCEDURES

In this article, H/V ratio are evaluated using different spectral techiniques applied to both earthquake and microtremor data (Nakamura technique). In particular, in order to avoid numerical instability, two different numerical techiniques are taken into account: (a) a smoothing procedure applied to the spectra of the seismogram components and (b) a regularization method applied to the H/V ratio (Landweber scheme). The data set consists of more than 70 earthquake events recorded by three component sensors displaced in the town of Fabriano (Central Italy) during the Umbria-Marche sequence started on September 1997. The local magnitudes range between 2.7 and 4.4, while the epicentral distances range between nearly 30 and 60 km. The stations were set to continuous recording so that a huge amount of microtremors was stored. The results are compared in terms of predominant frequencies and amplification levels in order to point out the influence of the adopted methods. The H/V ratio provides similar results if applied to a smoothed version of both earthquake and microtremor spectra, confirming that Nakamura technique is a cheap and a fast method to collect information on the site amplification effects. Moreover, the results relevant to earthquake data seem not to depend on the method used to stabilize the H/V ratio, whereas those relevant to microtremor data does. The explanation of this fact is suggested by the behaviour of the Landweber filter showing that the predominant frequency detected by means of microtremor data lies in a high instability region of the spectra.     

A COMPARATIVE STUDY ON THE ACTUAL AND ESTIMATED SEISMIC RESPONSE OF KIRALKIZI DAM IN TURKEY

Kiralkizi Dam, a 120 m high earthfill dam located in Diyarbakir city, Turkey, was shaken by a moment magnitude, M _w =4.6 earthquake at an epicentral distance of 8 km, on December 24, 2000, at 13:31 local time. The seismic response of the dam was assessed by using spectral ratios between (i) available crest and foundation records (C/F), (ii) horizontal and vertical components of the recorded motions (H/V), (iii) by performing 2 dimensional finite difference-based seismic response analyses (Flac-2D), and (iv) ID elastic shear beam solutions. First mode of vibration of the dam in the transverse direction by all four methods were estimated in the range of 0.55 to 0.62 second. Similar close agreement was not observed in higher modal periods estimated by H/V technique as compared to the predictions by C/F, Flac-2D, shear beam analysis techniques. Thus, H/V technique was concluded to be useful for the estimation of the fundamental resonance frequency of a soil structure, but not for its higher harmonics as consistent with available limited literature. In the longitudinal direction, natural period of the dam was estimated as 0.28 and 0. 82 second by H/V and C/F techniques, respectively. Such disagreement was explained by (i) differences in the definitions of the estimated periods, (ii) internal impedance contrast of the dam, (iii) contributions of 3D valley effects. Single seismometer record obtained from crest level was found to be inadequate for reliably assessing the response of a dam in the longitudinal direction, and it is recommended to install multiple seismometers both within dambody and the abutments. Last but not least, the results of these analyses were further compared by available accelograms recorded at three earthfill and rocknll dams from Japan. In general, it was concluded that the seismic response of Kiralkizi Dam is comparable and within the prediction ranges of available analyses methods and is consistent with the expected response of a dam this height.     

Site Classification Assessment for Estimating Empirical Attenuation Relationships for Central-Northern Italy Earthquakes

The aim of this article is to investigate the ground motion attenuation of the most industrialized and populated regions of Italy, evaluating the capability of different approaches to estimate site dependent models. The 5.2 local magnitude earthquake on November 24, 2004 shocked the areas of Northern Italy producing damage of about 215 million euros. The data set, including 243 earthquakes of local magnitude up to 5.2, has been collected in the period December 2002–October 2005 by 30 three-component seismic stations managed by Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Milano (INGV-MI). Empirical attenuation relationships have been estimated for horizontal peak ground velocity (PGHV), acceleration (PGHA), displacement (PGHD), and for response spectral acceleration (SA) for periods between 0.1 and 1.5 s. To estimate suitable attenuation models, in particular for sites characterized by thick sedimentary geological formations, a soil discrimination based on EU8 code can lead to wrong evaluations. On the contrary, a classification based on H/V spectral ratios of seismic ambient noise (NHV) allows the models to fit better real and predicted data and to reduce the uncertainties of the process. For each receiver, NHV have been strengthened by additional H/V spectral ratio of earthquake data (EHV), calculated considering different portions of the analysed signals. In order to validate the PGHA attenuation relationship for greater magnitudes, accelerometric records, relative to Central-Northern Italy strong motions occurring in the last 30 years, have been collected and superimposed to our attenuation curves.     

Numerical Evaluation of Seismic Soil Pressures on Rigid Walls Fixed to the Bedrock

Seismic soil pressures developed on a 7 m rigid retaining wall fixed to the bedrock are investigated using a finite element model that engages nonlinear soil intended materials available in OpenSees. This allows incorporation of the inelastic behavior of the soil and wave propagation effects in the soil-wall system seismic response. The nonlinear response of the soil was validated using the well-stablished, frequency-domain, linear-equivalent approach. An incremental dynamic analysis was implemented to comprehensively examine the effect of soil nonlinearity and input motion on the induced seismic pressures and to evaluate current code equations/methodologies at different levels of earthquake intensity. The results show that soil nonlinearity and seismic wave amplification may play an important role in the response of the soil-wall system. Therefore, methodologies that rely only on peak ground acceleration may introduce large bias on the estimated seismic pressures in scenarios where high nonlinearity and site amplification are expected.     

VERTICAL EARTHQUAKE GROUND MOTION RECORDS: AN OVERVIEW

Recent research clearly shows the importance of including the vertical component of earthquake ground motion in seismic analysis and design. In addition, pioneering studies [e.g., Elnashai and Papazoglou (1997)] have explored and documented the characteristics of available near-field vertical ground motion records. As a follow-up, this paper complements earlier studies, and investigates additional far-field records and available downhole array vertical motion records. A total of 111 free-field strong motion records (from California) and available downhole array records are employed. Compared to near-field records, far-field records generally contain more energy at longer periods. Based on the available data, response spectra are presented for near-field and far-field records respectively. The currently scarce downhole-array vertical motion records show that significant amplification may occur within the top 10-20 m soil layers. A simple one-dimensional (ID) vertical wave propagation model did not appear adequate for modelling the observed downhole array response. In using such a simplified model, very high viscous damping in the range of 15-25% was needed to match the recorded downhole vertical response, even for small tremors. Additional data and research are required [Beresnev et al., 2002] towards the development of a rational vertical motion site response analysis procedure.     

Analytical Collapse Study of Lightly Confined Reinforced Concrete Frames Subjected to Northridge Earthquake Ground Motions

Review of older non seismically detailed reinforced concrete building collapses shows that most collapses are triggered by failures in columns, beam-column joints, and slab-column connections. Using data from laboratory studies, failure models have previously been developed to estimate loading conditions that correspond to failure of column components. These failure models have been incorporated in nonlinear dynamic analysis software, enabling complete dynamic simulations of building response including component failure and the progression of collapse. A reinforced concrete frame analytical model incorporating column shear and axial failure elements was subjected to a suite of near-fault ground motions recorded during the 1994 Northridge earthquake. The results of this study show sensitivity of the frame response to ground motions recorded from the same earthquake, at sites of close proximity, and with similar soil conditions. This suggests that the variability of ground motion from site to site (so-called intra-event variability) plays an important role in determining which buildings will collapse in a given earthquake.     

Seismic Microzonation and Damage Assessment of Bam City,Southeastern Iran

This article aims to study the dynamic characteristics of soil in Bam, southeastern Iran. Fundamental frequency and amplification factor of soil sediment were estimated by microtremor measurement. This procedure was performed at 49 sites in the city. Two 120-second data were recorded at each site. Segmental cross spectra were applied to calculate spectrum, Nakamura's method (H/V) was used for analyzing data and fundamental frequency and amplification factor values were derived. Iso-frequency and iso-amplification maps of the city were prepared. Results show that soil type in Bam city is mainly stiff, although amplification factor is relatively large. There is a short period zone in northwest to southeast direction. Sediment depth was estimated using a correlation between fundamental frequency and sediment thickness. Results obtained from microtremors were compared to the geotechnical boreholes and it was shown that microtremors can be used for a rough estimation of sediment thickness. Damage distribution map of Bam due to the Bam earthquake on December 26, 2003 was prepared by both field and aerial investigation. Quality of buildings in Bam was evaluated, as a result a crude zoning map on the building quality was prepared. Finally, it was concluded that several factors have contributed to damage intensity variation; earthquake characteristics, local site effect, and buildings' quality. All mentioned parameters affected the damage variation.     

EMPIRICAL EVALUATION OF SITE EFFECTS BY MEANS OF H/V SPECTRAL RATIOS AT THE LOCATIONS OF STRONG MOTION ACCELEROMETERS IN ISRAEL

Acceleration data from local and regional earthquakes is of prime importance in evaluating the seismic hazard. Consequently, strong motion accelerometers are currently installed at more than 60 locations in Israel. We have explored the possibility of site amplification effects at 10 sites where local earthquakes triggered strong motion accelerometers by integrating empirical and analytical estimations. Implementing H/V spectral ratio techniques using 15 accelerograms from nine earthquakes, 105 seismograms shear-wave records of 35 local and regional earthquakes and seismograms of microtremors were used in the empirical evaluations. The subsurface models were constructed by integrating available geological and geophysical information at the analysed site with empirically evaluated site response functions. Amplification effects of factor 3-6 are observed at various frequencies in the 0.8-6.0 Hz band. Through the analysis process it became evident that the instant availability of many useful time windows of microtremors provides systematic estimations of the fundamental resonance frequency of each site and their associated amplification levels, which are similar to those obtained from H/V spectral ratios of seismograms and accelerograms and to those inferred from the subsurface geology. Analytical transfer functions should be reviewed with respect to empirical site response evaluations. Estimations that are based on only one approach may be totally misleading.     

Design Earthquake Ground Motion Prediction for Perth Metropolitan Area with Microtremor Measurements for Site Characterization

Perth is the largest city in Western Australia and home to three-quarters of the state's residents. In recent decades, there have been a lot of earthquake activities just east of Perth in an area known as the South-West Seismic Zone. Previous numerical results of site response analyses based on limited available geology information for PMA indicated that Perth Basin might amplify the bedrock motion by more than 10 times at some frequencies and at some sites. Hence, more detailed studies on site characterization and amplification are necessary. The microtremor method using spatial autocorrelation (SPAC) processing is a useful tool for gaining thickness and shear wave velocity (SWV) of sediments and has been adopted in many previous studies. In this study, the response spectrum of rock site corresponding to the 475-year return period for PMA is defined according to the probabilistic seismic hazard analysis (PSHA) based on the latest ground motion attenuation model of Southwest Western Australia. Site characterization in PMA is performed using two microtremor measurements, namely SPAC technique and H/V method. The clonal selection algorithm (CSA) is introduced to perform direct inversion of SPAC curves to determine the soil profiles of representative PMA sites investigated in this study. Using the simulated bedrock motion as input, the responses of the soil sites are estimated using numerical method based on the shear-wave velocity vs. depth profiles determined from the SPAC technique. The response spectrum of the earthquake ground motion on surface of each site is derived from the numerical results of the site response analysis, and compared with the respective design spectrum defined in the Australian Earthquake Loading Code. The comparison shows that the code spectra are conservative in the short period range, but may slightly underestimate the response spectrum at some long period range.     

SEISMIC ENERGY DEMAND IN STEEL MOMENT FRAMES

The present paper investigates the seismic energy demand in steel moment-resisting frames. The frames, with 3, 6 and 10 storeys, and 4 and 8 spans, are designed according to current seismic code provisions. The energy response (energy quantities and their distributions) in the frames subjected to an ensemble of six earthquake ground motions recorded on different soil conditions, is investigated by nonlinear time history analysis. The study concludes that (1) the results of energy response can be developed into a rational method of seismic evaluation and design for steel moment-resisting frames; (2) the energy concept based on the single-degree-of-freedom has limitations when extended to the realistic structural system for design purposes; and (3) it is necessary to develop the energy-based approach for seismic evaluation and design based on the seismic response of a realistic multi-degree-of-freedom structural system.     

Effect of Nonlinearity of Frame Buildings on Peak Horizontal Floor Acceleration

Using representative numerical models of eight code-designed steel moment-resisting frame buildings and several ground motions, time-history analyses are performed and a critical evaluation of Peak Horizontal Floor Acceleration (PHFA) demands is conducted. The frames are modeled alternatively as linear and nonlinear systems to isolate the effect of building nonlinearity on PHFA. In most cases, PHFA is reduced when nonlinear behavior of a building is considered; however, in some cases, significant amplification of PHFA is observed. Results from the numerical study provide insight into the trend of modal response modification factors presented taking ground motion spectral shape into account.     

Predicting the Seismic Response of Capacity-Designed Structures by Equivalent Linearization

An equivalent linearization procedure is developed for predicting the inelastic deformations and internal forces of capacity-designed structures under earthquake excitations. The procedure employs response spectrum analysis, and mainly consists of the construction of an equivalent linear system by reducing the stiffness of structural members that are expected to respond in the inelastic range. These members are well defined in structures designed with capacity principles. Maximum modal displacement demands of the equivalent linear system are determined either from the equal displacement rule, or from independent nonlinear response history analysis of SDOF systems representing inelastic modes.

Predictions obtained from the proposed equivalent linearization procedure are evaluated comparatively by using the results of nonlinear response history analysis as benchmark, linear elastic response spectrum analysis and conventional pushover analysis. The deformations and capacity controlled actions of a 12-story symmetrical plan concrete frame and a 6-story unsymmetrical plan concrete frame are obtained by each method under 96 strong ground motions. It is observed that the proposed procedure results in better accuracy in estimating the inelastic seismic displacement response parameters and capacity controlled forces than the other two approximate methods.     

Time-Domain Spectral Matching of Earthquake Ground Motions using Broyden Updating

Time-domain spectral matching of an earthquake ground motion consists of iteratively adding sets of wavelets to an acceleration history until the resulting response spectrum sufficiently matches a target spectrum. The spectral matching procedure is at its core a nonlinear problem because the addition of a wavelet often causes shifting in the time of peak response or creation of a larger second peak at a different time. A modification to existing time-domain spectral matching algorithms is proposed using Broyden updating for solving the set of nonlinear equations. Three wavelet bases are evaluated and the corrected tapered cosine wavelet is selected. The proposed algorithm is then tested and compared with other methods that are commonly used for spectral matching. The results show that the proposed algorithm is able to match the target spectrum while reasonably preserving the spectral nonstationarity, energy development, and the frequency content of the original time histories.     

RESPONSE OF UNREINFORCED MASONRY BUILDINGS

Abstract

A summary of dynamic measurements are presented that illustrate relations between linear seismic demand and true nonlinear response of unreinforced masonry buildings with flexible diaphragms and rocking piers subjected to a series of simulated earthquake motions.     

AN EVALUATION OF THE STRONG GROUND MOTION RECORDED DURING THE MAY 1, 2003 BINGÖL TURKEY,EARTHQUAKE

An important record of ground motion from a M6.4 earthquake occurring on May 1, 2003, at epicentral and fault distances of about 12 and 9 km, respectively, was obtained at a station near the city of Bingöl, Turkey. The maximum peak ground values of 0.55 g and 36 cm/s are among the largest ground-motion amplitudes recorded in Turkey. From simulations and comparisons with ground motions from other earthquakes of comparable magnitude, we conclude that the ground motion over a range of frequencies is unusually high. Site response may be responsible for the elevated ground motion, as suggested from analysis of numerous aftershock recordings from the same station. The mainshock motions have some interesting seismological features, including ramps between the P-and S-wave that are probably due to near- and intermediate-field elastic motions and strong polarisation oriented at about 39 degrees to the fault (and therefore not in the fault-normal direction). Simulations of motions from an extended rupture explain these features. The N10E component shows a high-amplitude spectral acceleration at a period of 0.15 seconds resulting in a site specific design spectrum that significantly overestimates the actual strength and displacement demands of the record. The pulse signal in the N10E component affects the inelastic spectral displacement and increases the inelastic displacement demand with respect to elastic demand for very long periods.     

Predicting Probabilistic-Based Strong Ground Motion Time Series for Citadel of Arg-E-Bam (South-East of Iran)

The main objective of this article is to present a probabilistic-based strong motion compatible with the source-path and site soil condition given the probability of exceedence for citadel of Arg-e-Bam site bed rock (South-East of Iran). A Fourier amplitude spectral attenuation relation for bed rock beneath the site is proposed which permits the estimation of time-histories through a probabilistic seismic hazard analysis procedure. Due to lack of data, the two well-known simulation techniques, point source and finite fault models have been used for generating hundreds of strong motion as input data. Tens of model parameter values such as stress-drop nucleation points were used, in each specified magnitude-distance, to reduce the uncertainty effects inherently existing in seismological/geological parameters. The proposed attenuation relation is validated by comparing the estimated strong motion, in the form of Fourier amplitude spectral, using the proposed attenuation relation with those of recorded ground motion data at three stations far away from the assumed source so that the results would not be influenced by the near source problems such as directivity and fling step. The results of proposed technique is assessed by comparing the estimated response spectra, with 10% probability of exceedence and 5% damping ratio, with those of traditional uniform hazard spectra. The proposed technique is supposed to be used in retrofitting procedure of international historical adobe structures in Arg-e-Bam site, which have been damaged during the destructive Bam earthquake 2003, Iran     

GROUND MOTION ANALYSIS IN THE ANCHORAGE BASIN: 1-D APPROACH

Understanding the dynamic behaviour of soil in the Anchorage basin in southcentral Alaska is essential for seismic hazard assessment of this highly seismically active region. The analysis of site responses for 40 sites from weak-motion and strong-motion data with amplitudes less than 0.1 g showed a strong influence of subsurface geological conditions on the characteristics of ground motion. Particularly, the sites in the central part of the city, including the downtown area, showed prominent resonance peaks around 1 Hz with amplification values up to about 4. The numerical analysis, based on one-dimensional multi-layer soil models shows that site response characteristics, and especially aforementioned peaks, are largely related to the thick, soft layer of Quaternary deposits, particularly cohesive fades of Bootlegger Cove Formation. The computed transfer functions for soil profiles of six representative sites are in accord with the site responses in the frequency range from 1 to 5 Hz. There is no significant change in amplification values below 2 Hz corresponding to large-amplitude (up to 0.38 g) ground motions; however, above 2-3 Hz the amplification values are greatly reduced in this case.     

Site Effects at Long Periods from Digital Strong Motion Records of the KiK-net,Japan

The high-quality digital records of the Japanese KiK-net were examined, with the aim of studying the influence of local site conditions on the displacement spectral ordinates at long periods. The attention was limited to those records for which the velocity profiles up to 100–200 m depth were known, and corresponding surface and borehole accelerograms were available. Based on the available records and with the support of 1D numerical simulations, different aspects that may have an influence on the amplification of long period spectral ordinates were studied, including the bedrock velocity profile, the site classification using V_s,30, and the earthquake magnitude.

Small amplification factors at long periods were found, ranging from 1 to 1.3, with median value from 1.05 to 1.1, for Eurocode 8 site classes B and C, respectively. Only for two records on soft soils (at the boundary between Eurocode 8 classes C and D), from small magnitude earthquakes, large amplification factors were obtained, up to about 4. A good correlation was found of the amplification levels with the response spectral ratio D(T₀)/D(10), where D(T₀) and D(10) are displacement spectral ordinates of the input signal at bedrock, at the fundamental period T₀ of the soil profile and at T = 10 s, respectively.     

Comparison of the Spectral Representation Method to Simulate Spatially Variable Ground Motions

The spectral representation method (SRM) is widely used when simulating spatially variable ground motions. It has mainly two formulas, i.e., the random amplitudes and the random phases formulas. There exist three methods for decomposing the cross spectral density matrix: Cholesky decomposition, eigen decomposition, and root decomposition. Therefore, there are six forms with respect to the different combinations of the simulation formulas and the decomposition methods. To provide researchers and engineers with the guidance on choosing simulation method, the six forms are systematically investigated from five aspects: the power intensity, response spectra, and stochastic error of auto/cross spectral density, Fourier spectra, and difference indexes for Fourier amplitudes and phases. Finally, we give the following advice: the characteristics of the ground motions simulated by the random amplitudes formula are independent of the decomposition method, while the characteristics of the ground motions simulated by random phases formula are dependent of the decomposition method. Furthermore, the root decomposition is strongly recommended when utilizing the random phases formula.