SEISMIC DESIGN REGULATION CODES: CONTRIBUTION OF K-NET DATA TO SITE EFFECT EVALUATION

The purpose of this study is to compare the site effect section of building codes (EC8 and UBC97) with the set of data provided by the Kyoshin network. In order to obtain a set of site coefficients and spectral shapes, we have first deduced an attenuation law for both horizontal and vertical motion. Site conditions are represented by the shear velocity averaged over the upper 30 m (V _s ³⁰). Our site classification (4 categories similar to those proposed in the new ECS and the UBC97) is based on borehole investigations at every station. This classification has permitted to distinguish clearly four response spectra which demonstrates the efficiency of V _s ³⁰ as characterising site conditions. Our law is then used to test site coefficients and spectral shapes of building codes ECS and UBC97. Concerning spectral shapes and site coefficients, our results are found to be in good agreement with EC8 and UBC97 only if category B (400<V _s ³⁰<800 m/s) is taken as reference. We also conclude that a site which is characterised as “rock” on geological criteria can not generally be classified in category A (V _s ³⁰>800 m/s). This suggests that classification in category A should be based only on field measurements. Concerning vertical motion, our analysis of the K-NET data shows that the ratio av/ah (vertical peak ground acceleration over horizontal peak ground acceleration) is between 0.50 and 0.68.     

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.     

ON THE RESOLUTION OF H/V MEASUREMENTS TO DETERMINE SEDIMENT THICKNESS,A CASE STUDY ACROSS A NORMAL FAULT IN THE LOWER RHINE EMBAYMENT,GERMANY

In recent years, H/V measurements have been increasingly used to map the thickness of sediment fill in sedimentary basins in the context of seismic hazard assessment. This parameter is believed to be an important proxy for the site effects in sedimentary basins (e.g. in the Los Angeles basin). Here we present the results of a test using this approach across an active normal fault in a structurally well known situation. Measurements on a 50 km long profile with 1 km station spacing clearly show a change in the frequency of the fundamental peak of H/V ratios with increasing thickness of the sediment layer in the eastern part of the Lower Rhine Embayment. Subsequently, a section of 10 km length across the Erft-Sprung system, a normal fault with ca. 750 m vertical offset, was measured with a station distance of 100 m. Frequencies of the first and second peaks and the first trough in the H/V spectra are used in a simple resonance model to estimate depths of the bedrock. While the frequency of the first peak shows a large scatter for sediment depths larger than ca. 500 m, the frequency of the first trough follows the changing thickness of the sediments across the fault. The lateral resolution is in the range of the station distance of 100 m. A power law for the depth dependence of the S-wave velocity derived from down hole measurements in an earlier study [Budny, 1984] and power laws inverted from dispersion analysis of micro array measurements [Scherbaum et aL, 2002] agree with the results from the H/V ratios of this study.     

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.     

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.     

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.     

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.     

Seismic Site Response Modeling for Three Missouri River Highway Bridges

Three highway bridges spanning the Missouri River flood plain were selected for evaluation of seismic site response for moderate size earthquakes emanating from the New Madrid Seismic Zone (NMSZ) in the Midwestern United States. The NMSZ is known to be capable spawning earthquakes larger than magnitude (M) 7.0, four of which occurred in a three-month period between 1811 and 1812, and the M_w 6.0 earthquake of October 1895 centered near Charleston, Missouri. This study evaluated the likely impacts of long period motion of these historic earthquakes on three long-span highway bridges using geotechnical data obtained from recent investigations. Our results suggest site amplification between 6× and 9×, depending on the magnitude and epicentral distance. We believe that threshold magnitude for serious foundation failure and damage to these bridges is between M_w 6.5 and 6.6. Above these magnitudes widespread liquefaction is predicted, which would effect the peak horizontal acceleration and spectral accelerations, causing the ground motions to be different than predicted. Increase in amplification of the response spectra also should be expected where the periods are higher than 1.0 sec. Therefore, M_w 6.5+ earthquakes at ranges 210–260 km could be expected to engender resonant frequency problems for multiple span bridges and tall buildings (10 to 25 stories) in channel corridors containing 20 to 46 m of unconsolidated sediment.     

EUROCODE 8 DESIGN RESPONSE SPECTRA EVALUATION USING THE K-NET JAPANESE DATABASE

The Eurocode 8 (EC8) currently proposes two standard shapes for the design response spectra. Type 1 spectra are enriched in long period and are suggested for high seismicity regions. Conversely, Type 2 spectra are proposed for low to moderate seismicity areas (like France), and exhibit both a larger amplification at short period, and a much smaller long period contents, with respect to Type 1 spectra. These propositions, however, were constrained using few events mostly recorded on analogical instruments. In the present study, we use the Japanese high quality digital K-net array in order to evaluate the proposed ECS response spectra. Furthermore, all K-net stations have geotechnical characterisation. We first constructed a database of shallow events, depth less than 25 km, to avoid subduction related records. The database spans six years of seismicity from 1996 until 2003. Thus, 591 events were selected with moment magnitude between 4 and 7.3, recorded at 691 stations, giving a total of 6812 two horizontal components accelerograms. Using these records, we computed spectral ground-motion prediction equations and we used them to review the shape of the proposed EC8 spectra. In particular, we studied the plateau-PGA ratio level, the period interval where this plateau is constant, and site amplification effects. The results show surprisingly that the Type 2 rock better envelope the Japanese data. Another interesting observation is that the K-net data corresponding to all soil classes are rich in short periods around 0.1 s. This characteristic has not been observed in other worldwide databases. Normalised empirical predictions show a widening of the plateau as the soil conditions degrade. This suggests that the Type 2 EC8 spectra do not cover enough the long periods for EC3-soil classes C, D and E. Finally, the computed ground-motion prediction equations show that the peak ground acceleration (PGA) is nearly invariant to the soil conditions. Soil effects are mainly seen in the shape and plateau level.     

FINITE FAULT MODELLING OF NEAR-FIELD ROCK MOTIONS IN THE NEW MADRID SEISMIC ZONE

Due to lack of strong motion records, point-source and finite-fault models have been used to simulate far-field motions at Memphis and St. Louis Cities from earthquake events in the New Madrid Seismic Zone. However, near-field rock motions and their associated uncertainties have never been studied within this zone. The objectives of this study are to develop a simple procedure to account for the uncertainty effect of earthquake source parameters, to analyze the sensitivity of near-field rock motions to input source parameters, and finally, to generate rock motions at two sites located within 11 km from the southwestern segment (strike-fault) and a third site bove the Reelfoot Rift (reverse fault) using a well-validated finite-fault simulation program; FINSIM. An equal-weight logic tree was developed to ensure that the assumed uncertainties are within physical, geological, and seismological constraints. For each site, 100 acceleration time histories with various combinations of parameter uncertainties were respectively simulated for an earthquake of M _w 7.0, 7.5, and 8.0 from each of the two faults. Their average spectral accelerations were in good agreement with those derived from the attenuation relation-ships representative to the Central and Eastern United States. Numerical simulations indicated that spectral accelerations are sensitive to the slip velocity, depth to top of fault, fault strike, slip distribution, and hypocentre location along the strike.     

NEW EMPIRICAL RESPONSE SPECTRAL ATTENUATION LAWS FOR MODERATE EUROPEAN EARTHQUAKES

This paper presents response spectral attenuation laws used in the new French Safety Rule, which is the reference for nuclear safety studies in France. Attenuation laws were derived from 965 horizontal and 485 vertical components from a two-step inversion method and accounts for geometrical spreading, anelastic attenuation and geological site condition. The datasets are mainly constituted of European strong motion records (83%) recently collected and homogeneously processed. In order to complete the distribution data beyond magnitude 6, a few American records were added, representing 17% of the datasets. The magnitude type and source-to-site distance definitions chosen to derive the laws are tested with respect to other definitions. These parametric tests induce a conservative law, for some magnitude and distance ranges of interest. The residual values between observed and predicted spectral accelerations are studied and do not exhibit any bias. The inferred laws are in good agreement with classical strong motion attenuation laws.     

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     

EXPERIMENTAL STUDY OF TOPOGRAPHIC AMPLIFICATION USING THE ISRAEL SEISMIC NETWORK

Earthquakes and microtremor records are used for estimating the site response of hard rock sites comprising four three-component seismic stations which operate as part of the Israel Seismic Network. The response functions are determined by implementing the horizontal-to-vertical component spectral ratio of earthquake shear-waves (receiver function estimates) and microtremors (Nakamura's estimate) observed simultaneously at the site. The sites of seismic stations ATZ (Mt. Atzmon), MBH (Mt. Berech) and MRNI (Mt. Meron) exhibit amplification attributed to topography effects. At ATZ, within the 1.3–2.0 Hz range, the amplification is in the order of factor 4. At MBH amplification levels of 3.0–3.5 are observed in the frequency range 1.5–4.0 Hz. Station MRNI exhibits a relatively strong amplification effect (up to 4) in the frequency range of about 2.5 to 3.5 Hz. Slight amplification around 5 Hz is observed at ATR (the proposed site for a nuclear power plant). These effects were correlated with the thickness of the weathered layer above unweathered chalk. A comparison between the amplification factor observed during earthquakes and those inferred from microtremors shows that these are, in general, in agreement. However, details of the spectral ratios from different microtremor recordings are not exactly the same. Differences appear mainly in the frequency at which the maximum amplification occurs. These observations demonstrate the usefulness of non-reference technique in estimating the topographical effects of ground shaking. These methods may be used in the process of seismic hazard assessment for ridges and mountain tops, common sites for settlements, communication relay stations, bridges, rope-drive and power transmission towers.     

Development and Verification of a New Site Classification System and Site Coefficients for Regions of Shallow Bedrock in Korea

In this article, site response analyses for 124 sites collected in Korea were performed to evaluate earthquake ground motions in regions of shallow bedrock in Korea. Based on the results, a new two-parameter site classification system was developed for use in regions of shallow bedrock. The system incorporates depth to bedrock and mean V_S of soil above bedrock as parameters for site classification instead of V_S30. Soil sites were divided into seven site classes and the corresponding site coefficients were proposed for each site category. Some verification work demonstrates the superiority of the proposed system over the current seismic codes.     

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.     

SEISMIC HAZARD ASSESSMENT OF UNITED ARAB EMIRATES AND ITS SURROUNDINGS

This paper presents the seismic hazard assessment and seismic zoning of the United Arab Emirates (UAE) and its surroundings based on the probabilistic approach. The area that has been studied lies between 50°E-60°E and 20°N-30°N and spans several Gulf countries. First, the tectonics of the area and its surroundings is reviewed. An updated catalogue, containing both historical and instrumental events is used. Seismic source regions are modelled and relationships between earthquake magnitude and earthquake frequency is established. A modified attenuation relation for Zagros region is adopted. Seismic hazard assessment is then carried out for 20 km interval grid points. Seismic hazard maps of the studied area based on probable Peak Ground Acceleration (PGA) for 10% probability of exceedance for time-spans of 50, 100 and 200 years are shown. A seismic zone map is also shown for a 475-year return period. Although the results of the seismic hazard assessment indicated that UAE has moderate to low seismic hazard levels, nevertheless high seismic activities in the northern part of UAE warrant attention. The northern Emirates region is the most seismically active part of UAE. The PGA on bedrock in this region ranges between 0.22 g for a return period of 475 years to 0.38 g for a return period of 1900 years. This magnitude of PGA, together with amplification from local site effect, can cause structural damage to key structures and lifeline systems.     

Mapping of an Archaeological Site at Olduvai Gorge,Tanzania

Abstract

An archaeological site at Olduvai Gorge, Tanzania, consists of a roughly quadrilateral surface of consolidated silt, about 25 feet by 47 feet. It shows a complex microtopography of pits and grooves. In order to record all these features a plan was made with a horizontal scale of 1:4 and contour interval of 0.1 foot (1.2 inches). The contours were interpolated between a large number of spot-heights. The positions of the spot-heights in the horizontal plane were fixed by distance measurement along the lines in a rectangular grid laid out from points on the margins of the surveyed area. The elevations of the spot-heights were determined by levelling with an automatic level. When the spot-heights had been plotted the map was mounted and taken into the field. Constant comparison was necessary to ensure that the interpolated contours provided an accurate representation of the complex topography.     

Saturated gravity wave spectra derived from foil chaff experiments performed during campaign ‘DYANA’ at 44°N (Biscarrosse,Southwest France)

Two series of measurements of horizontal wind velocity, one performed in winter 1990 at 44°N latitude, and the other in summer 1988 at 69°N latitude, were subjected to a spectral analysis. The two series covered a height interval of 11–15 km centred at about the same height (92.5km in winter, 91.5km in summer). In both cases the slopes of the spectra were close to −3. A dominant wavelength (if it existed) must have been larger than the height interval covered.     

SEGMENTAL CROSS-SPECTRUM AS A NEW TECHNIQUE IN SITE RESPONSE ESTIMATION USING SPECTRAL RATIO ANALYSIS

Different aspects of spectral analysis for site response evaluation are investigated in this study. The segmental cross-spectrum is proposed in spectral analysis of earthquake ground motions. The performance of segmental cross-spectrum in contrast with the conventional methods is investigated through the mathematical modelling, numerical analysis and application to earthquake data recorded at Chiba and Shinfuji downhole arrays in Japan. In analysis of earthquake data, the soil amplification function is identified using both uphole/downhole (U/D) and H/V spectral ratios. The advantage of seg-mental cross-spectrum is assessed by comparing identified amplification functions using different spectral methods and theoretical soil response. The reliability of site response estimations obtained by H/V spectral ratio using segmental cross- and Fourier spectra is also examined by means of cross-validation with the U/D spectral ratio of earthquake motion and theoretical soil response. Furthermore, the application of segmental cross-spectrum in nonlinear soil response is examined by comparing the amplification function of weak and strong motions for both methods. The results validate the advantage of segmental cross-spectrum in both linear and nonlinear soil response, particularly, when it used with H/V technique.     

ATTENUATION RELATION FOR WEST EURASIA DETERMINED WITH RECENT NEAR-FAULT RECORDS FROM CALIFORNIA,JAPAN AND TURKEY

Strong ground motion close to a fault can be expected to be very large, so its estimation is essential for human safetv. Although a few strong-motion data exist for the west Eurasian region, we proposed in a previous work [Berge-Thierry et al., 2003] an attenuation relation for spectral acceleration using strong-motion data recorded in west Eurasia (mainly in Europe) and some in the western United States: this relationship was derived for the French Safety Rule, which is applied for seismic hazard assessment at nuclear power plants. In this study, we propose a constraining of the amplitude saturation term related to the proximity of the fault, and an adding of an amplitude saturation term in the regression model. We add, to the data-set previously used to derive the west Eurasian attenuation relationship strong-motions recorded during recent large earth-quakes: the 1995 Hyogo-ken Nanbu (Kobe) event in Japan and the 1999 Kocaeli (Izmit) event in Turkey. The regression analysis, adopted from Fukushima and Tanaka [1990], is non-linear, so an iterative procedure is applied. The determined regression coefficients lead to a prediction of a peak ground acceleration of about 0.7 g for soil site conditions at a fault distance of 0.5 km. The Q coefficient deduced, from the distance coefficient is in agreement with scattering Q models. The introduction of the saturation term leads to significantly lower predictions of average spectral accelerations at short distances as compared with using the Berge-Thierry et al. [2003] empirical model.