Experimental study of seismic behavior of two hilly sites in Tehran and comparison with 2D and 3D numerical modeling

Two hilly sites were selected to study seismic site response due to topography effects. The sites were selected in a manner to be as much as possible homogenous and free of the soft soil layers effects. The hills were instrumented by nine velocimetric stations to record microtremors and the obtained data were analyzed using horizontal to vertical spectral ratios. Some standard spectral ratio tests were performed on noise as well. Then the instrumented hills were modeled (both 2D and 3D) assuming a linear elastic constitutive behavior subjected to vertically propagating SV and P Ricker wavelets. All calculations were performed in time domain using direct boundary element method. Different transfer function components, amplification patterns and spectral ratios were calculated in frequency domain. The frequency of vibration, obtained by experimental studies, is between 4 and 5 Hz for both of the hills. The spectral ratios derived by numerical simulations were compared with the observed spectral ratios. They show relatively good similarities between the results of these two methods. The frequencies of vibration derived from different methods seem to be nearly identical. The agreement in term of resonance frequency between microtremors and numerical modeling suggests that noise measurements could represent a simple, even if preliminary, tool in order to identify possible topographic amplification.     

Microtremor Measurements for the Microzonation of Dinar

v--vThe geotechnical site conditions in Dinar town located in western Turkey were investigated after the 1995 Dinar earthquake based on borings, in situ penetration tests, seismic wave velocity measurements, and microtremor records. The variation of damage distribution within the town was evaluated with respect to 23 district damage ratios calculated, based on the detailed damage survey conducted by the General Directorate of Disaster Affairs. Site amplifications were estimated from microtremor spectral ratios and microzonation was performed using a GIS methodology. The results of in situ penetration tests and seismic wave velocity measurements as well as the damage distribution were compared with the amplification zonation obtained from microtremor records. The results indicate the applicability of microtremor spectral ratios for assessing the local site conditions and site amplifications.     

Near surface geophysical surveys with a high frequency mutual impedance measuring system

Relationships for the mutual impedance of horizontal coplanar, vertical coplanar, perpendicular and vertical coaxial loop antenna arrays operating at a high frequency, located over a layered conducting medium, have been derived. The relationships take into account the occurrence of displacement currents in both the free space and the conducting medium. The relations can be used for calculating the earth parameters in geophysical prospecting. The use of high frequency antenna array makes it possible to determine distributions of apparent resistivity as well as apparent electrical permittivity. The potential application of the measurements of mutual impedance of loop antennas operating at a high frequency to the geophysical examinations is presented. The measuring system developed for this purpose is described. A simple method of interpreting the mutual impedance measurements to obtain the traces of apparent resistivity and apparent permittivity of the investigated earth is given.     

Mapping the spatial variation of soil water content at the field scale with different ground penetrating radar techniques

Two ground penetrating radar (GPR) techniques were used to estimate the shallow soil water content at the field scale. The first technique is based on the ground wave velocity measured with a bistatic impulse radar connected to 450 MHz ground-coupled antennas. The second technique is based on inverse modeling of an off-ground monostatic TEM horn antenna in the 0.8–1.6 GHz frequency range. Data were collected on a 8 by 9 m partially irrigated intensive research plot and along four 148.5 m transects. Time domain reflectometry, capacitance sensors, and volumetric soil samples were used as reference measurements. The aim of the study was to test the applicability of the ground wave method and the off-ground inverse modeling approach at the field scale for a soil with a silt loam texture. The results for the ground wave technique were difficult to interpret due to the strong attenuation of the GPR signal, which is related to the silt loam texture at the test site. The root mean square error of the ground wave technique was 0.076 m³ m⁻³ when compared to the TDR measurements and 0.102 m³ m⁻³ when compared with the volumetric soil samples. The off-ground monostatic GPR measured less within-field soil water content variability than the reference measurements, resulting in a root mean square error of 0.053 m³ m⁻³ when compared with the TDR measurements and an error of 0.051 m³ m⁻³ when compared with the volumetric soil samples. The variability between the two GPR measurements was even larger with a RSME of 0.115 m³ m⁻³. In summary, both GPR methods did not provide adequate spatial information on soil water content variation at the field scale. The main reason for the deviating results of the ground wave method was the poor data quality due to high silt and clay content at the test site. Additional reasons were shallow reflections and the dry upper soil layer that cannot be detected by the ground wave method. In the case of off-ground GPR, the high sensitivity to the dry surface layer is the most likely reason for the observed deviations. The off-ground GPR results might be improved by using a different antenna that allows data acquisition in a lower frequency range.     

Application of microtremor measurements to the estimation of earthquake ground motions in Kushiro city during the Kushiro-Oki earthquake of 15 january 1993

Microtremor measurements were conducted in Kushiro City, Hokkaido Island of Japan. The thickness of alluvial deposits in the central area of Kushiro City abruptly changes from 0 to 80 m. Predominant ground motion at a microtremor measuring site is represented by the spectral ratio, the ratio of Fourier amplitude spectrum of microtremor in the horizontal direction to that in the vertical direction. The peak frequency of the spectral ratio corresponds well to the predominant frequency estimated from the thickness of alluvial layer at each site. Based on a hypothesis that the spectral ratio can be regarded identical with a half of the amplification factor from diluvial bed rock to the ground surface, the strong ground motion due to the Kushiro-oki earthquake at each site is estimated. The method of estimation is verified by the comparison of the synthesized accelerogram with the recorded one at West Port of Kushiro. The results satisfactorily explain major damages and vibrations felt in Kushiro City during the earthquake.     

瑞利阻尼矩阵对深厚场地地震反应的影响分析

介绍场地时域分析基本理论和几种常见的阻尼矩阵确定方法,并基于有限差分方法,采用实际地震记录,考虑输入地震动强度和频谱特性,分别选择5种不同的瑞利阻尼矩阵形式,比较分析其对深厚场地时域动力反应的影响。结果表明:不同阻尼矩阵对加速度和反应谱均有影响,影响程度和输入地震动的频谱特性有关;当土层基频接近输入地震动卓越频率时,只根据土层基频确定阻尼矩阵是可行的,但当输入地震动卓越频率高于或远高于土层基频时,该方法会严重低估地表加速度峰值和反应谱谱值,时域分析中阻尼的确定应综合考虑土层特性和地震动频谱特性。     

Evaluation of local site effect from microtremor measurements in Babol City,Iran

Every year, numerous casualties and a large deal of financial losses are incurred due to earthquake events. The losses incurred by an earthquake vary depending on local site effect. Therefore, in order to conquer drastic effects of an earthquake, one should evaluate urban districts in terms of the local site effect. One of the methods for evaluating the local site effect is microtremor measurement and analysis. Aiming at evaluation of local site effect across the city of Babol, the study area was gridded and microtremor measurements were performed with an appropriate distribution. The acquired data was analyzed through the horizontal-to-vertical noise ratio (HVNR) method, and fundamental frequency and associated amplitude of the H/V peak were obtained. The results indicate that fundamental frequency of the study area is generally lower than 1.25 Hz, which is acceptably in agreement with the findings of previous studies. Also, in order to constrain and validate the seismostratigraphic model obtained with this method, the results were compared with geotechnical, geological, and seismic data. Comparing the results of different methods, it was observed that the presented geophysical method can successfully determine the values of fundamental frequency across the study area as well as local site effect. Using the data obtained from the analysis of microtremor, a microzonation map of fundamental frequency across the city of Babol was prepared. This map has numerous applications in designing high-rise building and urban development plans.     

The effects of torsion and motion coupling in site response estimation

Soil amplification characteristics are investigated using data from the Chibaken‐Toho‐Oki earthquake and its aftershocks recorded at Chiba dense array in Japan. The frequency‐dependent amplification function of soil is calculated using uphole‐to‐downhole spectral ratio analysis, considering the horizontal components of shear wave. The identified spectral ratios consistently demonstrate the splitting of peaks in their resonance frequencies and low amplification values in comparison with a 1D model. The torsional behaviour and horizontal ground motion coupling are clarified as the reasons for these phenomena at the site. To prove the hypothesis, the torsional motion is directly evaluated using the data of the horizontal dense array in different depths at the site. The comparison between Fourier spectra of torsional motion and identified transfer functions reveals the peaks at the same frequencies. The wave equation including torsion and horizontal motion coupling is introduced and solved for the layered media by applying wave propagation theory. Using the developed model, the effects of torsional motion with horizontal motion coupling on soil transfer function are numerically examined. Splitting and low amplification at resonance frequencies are confirmed by the results of numerical analysis. Furthermore, the ground motion in two horizontal directions at the site is simulated using site geotechnical specification and optimizing the model parameters. The simulated and recorded motions demonstrate good agreement that is used to validate the hypothesis. In addition, the spectral density of torsional ground motions are compared with the calculated one and found to be well predicted by the model. Finally, the results are used to explain the overestimation of damping in back‐calculation of dynamic soil properties using vertical array data in small strain level. Copyright © 2003 John Wiley & Sons, Ltd.     

Seismic vulnerability estimation of the building considering seismic environment and local site condition

Introduction The estimation of damage probability distribution among different damage states of rein-forced concrete buildings is a key component of earthquake loss estimation for modern city or a group of cities. With the development of city, the reinforced concrete buildings are major compo-nent parts of modern cities. Vulnerability estimates for these kinds of buildings are of importance to those responsible for civil protection, relief, and emergency services to enable adequate contin-genc…     

Bayesian inference of the Cole–Cole parameters from time- and frequency-domain induced polarization

The inversion of induced‐polarization parameters is important in the characterization of the frequency electrical response of porous rocks. A Bayesian approach is developed to invert these parameters assuming the electrical response is described by a Cole–Cole model in the time or frequency domain. We show that the Bayesian approach provides a better analysis of the uncertainty associated with the parameters of the Cole–Cole model compared with more conventional methods based on the minimization of a cost function using the least‐squares criterion. This is due to the strong non‐linearity of the inverse problem and non‐uniqueness of the solution in the time domain. The Bayesian approach consists of propagating the information provided by the measurements through the model and combining this information with a priori knowledge of the data. Our analysis demonstrates that the uncertainty in estimating the Cole–Cole model parameters from induced‐polarization data is much higher for measurements performed in the time domain than in the frequency domain. Our conclusion is that it is very difficult, if not impossible, to retrieve the correct value of the Cole–Cole parameters from time‐domain induced‐polarization data using standard least‐squares methods. In contrast, the Cole–Cole parameters can be more correctly inverted in the frequency domain. These results are also valid for other models describing the induced‐polarization spectral response, such as the Cole–Davidson or power law models.     

New relationship between fundamental site frequency and thickness of soft sediments from seismic ambient noise

In this contribution, new relationship between the fundamental site frequency and the thickness of soft sediments is obtained for many sites in Egypt. The Horizontal-to-Vertical Spectral Ratio (“H/V”) technique (known as Nakamura technique) can be used as a robust tool to determine the thickness of soft sediments layers overlaying bedrock from observations and measurements of seismic ambient noise data. In Egypt, numerous seismic ambient noise measurements have been conducted in several areas to determine the dynamic properties of soft soil for engineering purposes. At each site in each studied area, the fundamental site frequency was accurately estimated from the main peak in the spectral ratio between the horizontal and vertical component. Consequently, an extensive database of microtremor measurements, well logging data, and shallow seismic refraction data have been configured and assembled for the studied areas. New formula between fundamental site frequency (f₀) and thickness of soft sediments (h) is established. The new formula has been validated and compared with other formulas of earlier scientists, and the results indicate that the calculated depth and geometry of the bedrock surface using new formula are in a good agreement with well logs data and previously published seismic refraction surveys in the investigated sites.     

A comparison of surface and underground array measurements of ambient noise recorded in Naples (Italy)

In this study, we describe two experiments of seismic noise measurements carried out in Naples, Italy. The site allowed measurements to be obtained both at the surface and in a tunnel that is 120-m-deep. The main goal was to compare the seismic response evaluated at the surface to the in-tunnel response, through spectral, polarization, and resonance directivity analyses. In the 1 to 20 Hz frequency band, the noise level was up to 15 dB higher at the surface than in the tunnel. The polarization properties and horizontal-to-vertical spectral ratios appear not to be influenced by the tunnel geometry or by the topography. Some preferential alignments were observed in the polarization azimuths computed at the surface, which are likely to be due to local sources, rather than morphological features. The absence of directivity effects and the low noise levels in the tunnel make this site suitable for installing seismic stations. We also studied how the subsoil structure affects the seismic motion at the surface. The dispersive properties of the Rayleigh waves were investigated using the spatial autocorrelation method. A joint inversion of the dispersion data and the horizontal-to-vertical spectral ratios provided the subsurface Vs profile. The derived model has a low velocity contrast at depth, such as to generate moderate and broad H/V spectral ratio peak amplitude. The normalized spectral ratio appears more appropriate to identify the soil-resonance frequencies.     

Spectral bandwidth enhancement of GPR profiling data using multiple-frequency compositing

The amplitude spectrum of ground penetrating radar (GPR) reflection data acquired with a particular antenna set is normally concentrated over a spectral bandwidth of a single octave, limiting the resolving power of the GPR wavelet. Where variously-sized GPR targets are located at numerous depths in the ground, it is often necessary to acquire several profiles of GPR data using antennas of different nominal frequencies. The most complete understanding of the subsurface is obtained when those frequency-limited radargrams are jointly interpreted, since each frequency yields a particular response to subsurface reflectivity. The application of deconvolution to GPR data could improve image quality, but is often hindered by limited spectral bandwidth.We present multiple-frequency compositing as a means of combining data from several frequency-limited datasets and improving the spectral bandwidth of the GPR profile. A multiple-frequency composite is built by summing together a number of spatially-coincident radargrams, each acquired with antennae of different centre frequency. The goal of the compositing process is therefore to produce a composite radargram with balanced contributions from frequency-limited radargrams and obtain a composite wavelet that has properties approximating a delta function (i.e. short in duration and having a broad, uniform spectral bandwidth).A synthetic investigation of the compositing process was performed using Berlage wavelets as proxies for GPR source pulses. This investigation suggests that a balanced, broad bandwidth, effective source pulse is obtained by a compositing process that equalises the spectral maxima of frequency-limited wavelets prior to summation into the composite. The compositing of real GPR data was examined using a set of 225, 450 and 900 MHz GPR common offset profiles acquired at a site on the Waterloo Moraine in Ontario, Canada. The most successful compositing strategy involved derivation of scaling factors from a time-variant least squares analysis of the amplitude spectra of each frequency-limited dataset. Contributions to the composite from each nominal acquisition frequency are clear, and the trace averaged amplitude spectrum of the corresponding composite is broadened uniformly over a bandwidth approaching two-octaves. Improvements to wavelet resolution are clear when a composite radargram is treated with a spiking deconvolution algorithm. Such improvement suggests that multiple-frequency compositing is a useful imaging tool, and a promising foundation for improving deconvolution of GPR data.     

Assessing the state of the atmospheric acoustic channel using the IDEAS data on long-distance microbarom propagation

This paper discusses the technique for assessing the state of atmospheric acoustic channels (AACs) for the long-distance propagation of microbaroms. We calculated two possible microbarom propagation paths to station “Badary” from the sources located in (1) the North Atlantic at an azimuth of 320° and (2) the Northwestern Pacific at an azimuth of 60°. We investigate the spatio-temporal structure of the AAC. The experimental data are compared with modeling results.     

统一考虑地震环境和局部场地影响的建筑物易损性研究

提出了一种综合考虑地震环境和场地影响的钢筋混凝土房屋地震易损性分析方法. 将地震环境、局部场地和工程结构作为一个整体,以概率地震危险性分析的方式考虑地震环境的影响,在此基础上详细考虑了随局部场地而变化的反应谱形状对结构地震反应及其破坏概率分布的影响. 此外,还提出了另一种表述结构易损性的方式,以对应于不同超越概率地震危险水平的方式, 提供结构地震破坏概率分布的信息.      

Evaluation of the influence of experimental conditions on H/V results from ambient noise recordings

The H/V-noise technique is now widely used to estimate site effect parameters (fundamental frequency and sometimes the associated soil amplification), and many surveys using this technique have provided convincing results. However, a general agreement on a methodology for data acquisition, data processing and result interpretation has yet to be found. H/V measurements from ambient noise recordings imply both reliability of the results and rapidity of data collection. It is therefore important to understand which experimental conditions (1) influence data quality and reliability, and (2) can help speeding up the recording process. Within the framework of the SESAME European project, a specific task was defined to investigate the reliability of the H/V spectral ratio technique in assessing the site effects. The aim of WP02, one specific Work Package of the SESAME project, is to study the effects of experimental conditions on both stability and reproducibility of H/V results. This study has been conducted in a purely experimental way, by testing the possible influence of various experimental conditions on H/V results both on the frequency peak value and on its amplitude. WP02 results help setting up the experimental conditions under which ambient noise recordings have to be performed in order to provide reproducible, reliable and meaningful H/V results. In this paper we present the results of the WP02 SESAME project concerning the evaluation of the influence of experimental conditions of ambient noise recording on H/V results.     

Use of microtremor measurement for assessing site effects in Northern Belgium – interpretation of the observed intensity during the MS = 5.0 June 11 1938 earthquake

To provide quantitative information on site effects in the northern part of Belgium, forty-seven H/V microtremor measurements were performed with 5 second seismometers over an area of about 15.000 km². Most of the results show a northward regular increase of the fundamental period in agreement with the augmentation of the Mesozoic and Cainozoic soft sedimentsthickness from a few meters 40 km south of Brussels to 900 m at the Netherlands-Belgium border. The measured resonance frequency values wereconsistent with theoretical computations performed at different sites onthe basis of existing information and shallow seismic experiments. At oneparticular site (Uccle) where borehole data were available, microtremor measurements using an array of four seismological stations with differentapertures allowed to obtain the low frequency part of the Rayleigh wave dispersion curve, extending the range covered by the analysis of surfacewaves artificially generated. The Vs profile derived from the surface waveinversion corroborates the 1 Hz natural frequency of the site. Comparison of these results with the macroseismic information concerning the M_S = 5.0 1938 earthquake which occurred 50 km west of Brussels,confirmed the hypothesis that the geological structure of the Brabant massifis likely to control damage distribution during such an earthquake. Comparisonbetween the intensity map of the 1938 earthquake and the resonance period ofsediments obtained by our microtremor study shows a clear relation betweenthe two parameters. During the 1938 earthquake, site effects played a prominent role due to the dimension of the source whose corner frequency wasabout 1 Hz.     

Determination of S-wave Site Response in Anchorage, Alaska in the 1–9 Hz Frequency Band

—?Site response was estimated at 19 sites in the Anchorage basin in south-central Alaska, using 15 local earthquakes recorded with good signal-to-noise ratio by a temporary weak motion network. The receiver-function-type horizontal-to-vertical spectral ratios (HVSR) were computed at 1–9 Hz frequency band and the resulting HVSR contour maps at 1, 5 and 9 Hz are presented here. The spatial site response distribution shows considerable variation from the foothills of the Chugach Mountains in the east to the western part of Anchorage. The site response increases by a factor of 3 and 3.5 at 1 and 5 Hz, respectively, from the area of older glacial deposits in the eastern part of the city to the area occupied by the Bootlegger Cove formation, particularly in the section adjoining Knik Arm. At 9 Hz, the variation of HVSR from the east to the west is smaller, approximately by a factor of 2. Moreover, the trend of the HVSR variation at 1 and 5 Hz shows good correlation with that of the soil class obtained from surface measurements of S-wave velocity in the 0–30?m depth range and available results on ground failure susceptibility of Anchorage.