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.     

Spectral Attenuation Relations at Soft Sites Based on Existing Attenuation Relations for Rock Sites

A simple and general method based on well-known random vibration theory is used to compute spectral attenuation relations at soft sites based on existing spectral attenuation relations at rock sites. The method consists of: (1) computation, for given magnitude and distance, of the expected Fourier amplitude spectrum associated with the median rock response spectrum computed with the attenuation relation; (2) inclusion of site effects characterized by a frequency-dependent, linear, or nonlinear transfer function; and (3) computation of the response spectrum at the soft site.     

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.     

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.     

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.     

A PROCEDURE FOR ESTIMATING INPUT ENERGY SPECTRA FOR SEISMIC DESIGN

In this paper, the damage potential of an earthquake ground motion is evaluated in terms of the total power of the acceleration of the ground motion. By assuming an appropriate spectral shape for the input energy spectrum, and using the well-known Parseval theorem for evaluating the total power of a random signal, the peak amplification factor for the equivalent input energy velocity spectrum can be determined. It is shown that the peak amplification factor for the input energy spectrum depends on the peak-ground-acceleration to peak-ground-velocity ratio and duration of the strong motion phase of the ground motion. Values for the equivalent input energy velocity amplification factor vary from about 2 to 10 for most of the recorded ground motions used in this study. Although a considerable scatter of data is observed in this study, the peak amplification factor predicted by the Fourier amplitude spectrum of the ground acceleration provides a fairly good estimate of the mean value of the peak input energy compared to that determined from inelastic dynamic time history analyses, particularly for systems with high damping and low lateral strength. The peak amplification factor derived in this paper provides a more consistent approach for estimation of seismic demand when compared to an earlier empirical expression used for the formulation of duration-dependent inelastic seismic design spectra, even though only a slight difference in the required lateral strength results from the use of the new formula.     

MODIFICATION OF SITE RESPONSE IN PRESENCE OF LOCALISED SOFT LAYER

Experimental and numerical simulations are performed to evaluate the modification of ground response resulting from either the presence of soft layers or occurrence of partial liquefaction. Results from two densely instrumented dynamic centrifuge tests are presented to show the ambiguous role played by the presence of a soft layer. It was found that the lateral extent of the soft layer has significant influence on the overall response of the layered strata and any structure founded on it. The experimental observations are supported by simplified numerical analysis. The amplification or deamplification of the input motion is found to be a function of the ratio of the width of soft layer to the wave length. Based on the numerical analysis, a general function describing the site amplification is presented which may be used as a guide in seismic design of foundations in such layered strata.     

TORSIONAL EFFECTS IN THE PUSHOVER-BASED SEISMIC ANALYSIS OF BUILDINGS

The general trends of the inelastic behaviour of plan-asymmetric structures have been studied. Systems with structural elements in both orthogonal directions and bi-axial eccentricity were subjected to bi-directional excitation. Test examples include idealised single-storey and multi-storey models, and a three-storey building, for which test results are available. The response in terms of displacements was determined by nonlinear dynamic analyses. The main findings, limited to fairly regular and simple investigated buildings, are: (a) The amplification of displacements determined by elastic dynamic analysis can be used as a rough, and in the majority of cases conservative estimate in the inelastic range, (b) Any favourable torsional effect on the stiff side, which may arise from elastic analysis, may disappear in the inelastic range. These findings can be utilised in the approximate pushover-based seismic analysis of asymmetric buildings, e.g. in the N2 method. It is proposed that the results obtained by pushover analysis of a 3D structural model be combined with the results of a linear dynamic (spectral) analysis. The former results control the target displacements and the distribution of deformations along the height of the building, whereas the latter results define the torsional amplifications. The proposed approach is partly illustrated and evaluated by test examples.     

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.     

THE USE OF SIMPLE PULSES TO ESTIMATE INELASTIC RESPONSE SPECTRA

Inelastic response spectra are estimated for elasto-plastic SDOF systems subjected to strong earthquake ground motions by applying the strength reduction factors determined for a simple pulse to the elastic response spectrum of the ground motion. This approach relies upon similarities in the strength reduction factors computed for earthquake ground motions and for short duration pulses. The accuracy of the estimated inelastic spectra obtained using 24 simple pulse waveforms is assessed in order to identify subsets of just several pulse waveforms that are suited for this purpose. Based upon the ground motions and pulses investigated, this approach appears to be equally applicable to short and long duration ground motions and those having near-fault forward directivity features.     

EVALUATION OF ACCIDENTAL ECCENTRICITY IN BUILDINGS DUE TO ROTATIONAL COMPONENT OF EARTHQUAKE

This study adopts a random procedure in the evaluation of the effect of the rotational component of earthquake on the accidental eccentricity of symmetric and asymmetric buildings. The spectral density function of the rotational component of earthquake acceleration (about the vertical axis) is obtained on the basis of the spectral density function of the horizontal component of earthquake acceleration. The rotational component of an earthquake can increase the response of the structure. The degree of the increase is highly dependent upon the dynamic characteristics of the system and the rotational component of the earthquake. To bring this increase under consideration, seismic codes represent a parameter referred to as accidental eccentricity, as a part of the design eccentricity. The purpose of the present study is to estimate the value of this increase and to make appropriate suggestions based on frequency domain analysis.     

ON THE EARTHQUAKE RESPONSE OF ROTORS

In this paper, the influence of the rotation speed on the earthquake response of rotors is investigated by means of a simple model with only two degrees of freedom. Some worth noticing properties of the eigenmodes and the transfer functions in the presence of solid gyroscopic effect are discussed. In particular, the importance of the relative contribution of the forward and backward modes is evaluated using basic properties of the transfer functions. In the case of a white noise excitation the backward mode contribution dominates the response, for high values of spin velocity. Depending on the frequency content of the excitation and the rotation speed, an important amplification of the response may occur due to gyroscopic effects, as a result of the amplification of the backward mode response. Nonlinear boundary conditions reduce, to some degree, the differences between the responses of rotating and non-rotating systems, but, even in this case, the importance of the contribution of the backward mode in the response appears clearly. Rotors immersed in a dense fluid with moderate confinement exhibit different behaviour and amplification is possible in the case of resonant response of the second forward mode.     

Floor Spectra of Inelastic RC Frame Buildings Considering Ground Motion Characteristics

A range of reinforced concrete frame buildings with different levels of inelasticity as well as periods of vibration is analyzed to study the floor response. The derived floor acceleration response spectra are normalized by peak ground acceleration, peak floor acceleration, and ground response spectrum. The normalization with respect to ground response spectrum leads to the lowest coefficients of variation. Based on this observation as well as previous studies, an amplification function is proposed that can be used to develop design floor spectra from the ground motion spectrum, considering the building’s dynamic characteristics and level of inelasticity.     

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.     

WEST OF THE BEST: ROCK ART AND ARCHAEOLOGICAL DISCOVERIES IN THE DORO !NAWAS REGION OF NORTHWEST NAMIBIA

The Northwest of Namibia bears a concentration of rock art that is among the highest in the world. One of the centres is Twyfelfontein, a cluster of sites surrounding a waterhole at the eastern fringe of the Namib Desert, classified as the prime rock engraving site in southern Africa and designated a UNESCO World Heritage Site. Recent prospections in the desert west of Twyfelfontein have resulted in the discovery of engravings and paintings that exceed those at the UNESCO site, both in quantity and quality. Our contribution to this volume in memory of Klavs Randsborg introduces this newly discovered archaeological landscape and discusses the connection between water ‐ a rare commodity in the desert environment ‐ and the rich array of rock art and other archaeological features found here.     

EXPLORING IDENTITY AND PLACE: AN ANALYSIS OF THE PROVENANCE OF PASSAGE GRAVE STONES ON GUERNSEY AND JERSEY IN THE MIDDLE NEOLITHIC

Summary. This study examines the provenance of rock types used in the construction of Middle Neolithic passage grave stones on the islands of Guernsey and Jersey, focusing on the social dimensions of stone selection. The use of stones in passage grave construction includes both local and non-local rock types, which at some sites are organized in distinctive patterns. It is argued that the choice of stones was bound by concepts of identity, and that the communities which gathered to build these monuments may have used specific rock types to represent their community and their local mythologies. The relationship between identity and stone selection is supported both by analogy and by research into the role of landmarks in the development of local landscapes and ideology. The success of megalith provenance studies on Guernsey and Jersey suggests considerable potential for future research in other geologically diverse regions.     

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.     

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.     

QUICK SURVEY OF THE POSSIBLE CAUSES OF DAMAGE ENHANCEMENT OBSERVED IN SAN GIULIANO AFTER THE 2002 MOLISE,ITALY SEISMIC SEQUENCE

On October 31 and November 1, 2002, two earthquakes of magnitude 5.4 and 5.3 hit the area at the border between the Molise and Puglia regions in Southern Italy. The damage pattern in the epicentral area qualified the quake as an intensity VII MCS event, although providing a notable exception relevant to the small village of San Giuliano di Puglia. Since the first macroseismic survey, it appeared clear that in S. Giuliano the intensity was two degrees higher with respect to three neighbouring villages located within a radius of 3 km. Soon after the quake, our team started a campaign of microtremor HVSR measurements (Horizontal to Vertical Spectral Ratio), then we installed accelerometers and carried out damage and geological surveys. Finally, we performed a geoelectrical tomography and two profiles of Vs velocity with depth using the NASW technique (Noise Analysis of Surface Waves). The preliminary observations indicate that ground motion amplification is present in S. Giuliano within the frequency band that may affect building. A strong velocity contrast 20 m deep causes the predominant peak. More amplification could be due to more complicated, 2D effects. As regards the damage pattern, it divides S. Giuliano in three zones showing different characteristics and seismic behaviour. A building-by-building survey is still under way to better evaluate vulnerability variations in different zones of the village. However, the acquired data so far is sufficient to propose site amplification as a possible cause of the damage enhancement observed in S. Giuliano.     

SITE AMPLIFICATION FACTOR IN ANCHORAGE,ALASKA BASED ON LONG PERIOD MICROTREMORS

Long period microtremors with periods ranging from 0.5 to 10 seconds were measured in the Anchorage metropolitan area. Two horizontal components of motion were recorded at 81 sites uniformly distributed throughout the basin with spatial resolution of about 2 km. Recording at each site was done for 300 seconds with a sampling rate of 20 Hz. Repeated measurements were performed at a bedrock reference site simultaneously with the measurements in the field. The measurements were completed in six days. In addition, multiple recordings were obtained concurrently at the reference bedrock site and a sediment site. Based on these measurements the Fourier spectra were calculated for each of the site. Ground motion amplification is determined in terms of spectral ratio of horizontal spectral amplitudes at a sediment site and the reference bedrock site. Mean spectral ratio contours were evaluated for different period bands. The results show that for period band 3 to 5 seconds the spectral ratio contours agree well with the ground failure susceptiblity map of Anchorage.