Structure, soil–structure response and effects of damage based on observations of horizontal-to-vertical spectral ratios of microtremors

The microtremor horizontal-to-vertical-spectral-ratio (HVSR) technique is widely used in the urban environment to assess the fundamental frequency response of the ground. Extensive literature exists about case histories using HVSR for microzonation in several cities, but no systematic studies have been devoted to check the presence of soil–structure interaction effects, and even less attention to study building behaviour after earthquake damage. To evaluate the above-mentioned effects, a series of experiments are reported in this article.We first made a series of microtremor measurements on buildings and civil structures to evaluate the reliability of fundamental frequency determinations. Then, we considered several case studies to evaluate the effect of soil–structure interaction in estimates of site response in the presence of tall buildings. Finally, an experiment on the frequency change due to damage was performed. It was possible to confirm that HVSR is able to detect building fundamental modes and once known the building frequency, it is also possible to detect the presence of soil–structure interaction. Thus, once the presence of the building natural frequency is identified, it is possible to infer the site response from free field measurements. We also found that the HVSR technique is equally useful for detecting structural damage by determining the frequency shift of the buildings.     

Soil-structure interaction effects on an instrumented building

The paper highlights the use of fem and bi-directional lumped-mass-storey-stiffness numerical models for the study of the soil–structure interaction (ssi) effects on an instrumented building. Data on the structural response have been obtained through the project for seismic instrumentation of a 16-storey r/c cast-in-place dwelling building (Chisinau, Republic of Moldova) during a series of earthquakes (Gutenberg–Richter M _GR = 5.0−6.7). The effect of soil–structure interaction is clearly observed comparing the responses recorded on foundation and free-field. ssi becomes more pronounced for higher level of ground shaking amplifying the natural period of the structure and slightly suppressing high frequences on the foundation in comparison with the free-field motion.     

Wave propagation in a seven-story reinforced concrete building: III. Damage detection via changes in wavenumbers

This paper presents low frequency wavenumbers in a seven-storey reinforced concrete building estimated from its recorded response to eleven earthquakes, one of which (1994 Northridge) caused visible structural damage, and two of which are its aftershocks. The wavenumbers, K_i,j(f), are estimated from pairs (i,j) of records at neighboring recording sites in the building, distributed vertically or horizontally. Changes in K_i,j(f) from one event to another are compared in the undamaged (lower) and in the damaged (upper) part of the building, with the aim to find whether trends in K_i,j(f) can indicate damage. The results suggest significant and permanent increase of the wavenumbers in the damaged parts for the 1994 Northridge earthquake and its aftershocks, which is not the case for the other events in the damaged parts, and for all eleven events in the undamaged parts of the building. This increase in wavenumbers in the damaged parts can be explained by reduced wave velocities through the damaged structural members, and by scattering of waves from the discontinuities created by the damage. It is concluded from this qualitative analysis that wavenumbers estimated from strong motion recordings in a building can indicate location of damage, and that it would be useful to refine further this method (extend it to higher frequencies, and add the capability to quantify the damage). However, this would require more dense strong motion instrumentation in buildings than currently available. Deployment of dense arrays in selected buildings would provide data for further work on this subject.     

Ambient vibration tests of a seven-story reinforced concrete building in Van Nuys, California, damaged by the 1994 Northridge earthquake

Results of two detailed ambient vibration surveys of a 7-story reinforced concrete building in Van Nuys, California, are presented. Both surveys were conducted after the building was severely damaged by the 17 January 1994, Northridge earthquake (M_L=5.3, epicenter 1.5 km west from the building site) and its early aftershocks. The first survey was conducted on 4 and 5 February 1994, and the second one on 19 and 20 April 1994, about one month after the 20 March aftershock (M_L=5.3, epicenter 1.2 km north–west from the building site). The apparent frequencies and two- and three-dimensional mode shapes for longitudinal, transverse and vertical vibrations were calculated. The attempts to detect the highly localized damage by simple spectral analyses of the ambient noise data were not successful. It is suggested that very high spatial resolution of recording points is required to identify localized column and beam damage, due to the complex building behavior, with many interacting structural components. The loss of the axial capacity of the damaged columns could be seen in the vertical response of the columns, but similar moderate or weak damage typically would not be noticed in ambient vibration surveys. Previous analysis of the recorded response of this building to 12 earthquakes suggests that, during large response of the foundation and piles, the soil is pushed sideways and gaps form between the foundation and the soil. These gaps appear to be closing during “dynamic compaction” when the building site is shaken by many small aftershocks. The apparent frequencies of the soil–foundation–structure system appear to be influenced significantly by variations in the effective soil–foundation stiffness. These variations can be monitored by a sequence of specialized ambient vibration tests.     

Wave propagation in a seven-story reinforced concrete building: I. Theoretical models

For transient, high frequency, and pulse like excitation of structures in the near field of strong earthquakes, the classical design approach based on relative response spectrum and mode superposition may not be conservative. For such excitations, it is more natural to use wave propagation methods. In this paper (Part I), we review several two-dimensional wave propagation models of buildings and show results for theoretical dispersion curves computed for these models. We also estimate the parameters of these models that would correspond to a seven-story reinforced concrete building in Van Nuys, California. Ambient vibration tests data for this building imply vertical shear wave velocity β_z=112 m/s and anisotropy factor β_x/β_z=0.55 for NS vibrations, and β_z=88 m/s and β_x/β_z=1 for EW vibrations. The velocity of shear waves propagating through the slabs is estimated to be about 2000 m/s. In the companion paper (Part II), we estimate phase velocities of vertically and horizontally propagating waves between seven pairs of recording points in the building using recorded response to four earthquakes.     

Evidence of nonlinear site response in horizontal-to-vertical spectral ratio from near-field earthquakes

The issue is addressed as to whether the horizontal-to-vertical spectral ratio (HVSR) method is sensitive to the amplitude of ground motion from near-field earthquakes. Twenty-one three-component accelerograms from two closely located similar soil sites in the town of Lefkas are used. The recordings represent 17 earthquakes covering a wide range of magnitudes, epicentral distances and azimuths. Peak horizontal accelerations (PGA) and velocities (PGV) lie in the ranges 20–540 cm/s² and 1.4–55.2 cm/s. For each HVS ratio, the site's fundamental-resonance frequency, f_res, is determined visually. Linear correlation analysis shows that f_res is strongly (negatively) correlated to PGA and PGV (r between −0.7 and −0.8); no correlation is found with resonance amplitude or epicentral distance. We show that the observed correlation is attributable to soil nonlinearity and indicate how weak-motion estimates of f_res can be corrected for use in assessing site response during strong shaking.     

An inversion of site response and Lg attenuation using Lg waveform

Based on spectral ratio method, a joint inversion method was used to obtain parameters of Lg wave attenuation and site response. The inversion method allows simple and direct (two-parameter) determination of Lg wave attenua- tion and site response from sparse spectral data, which are not affected by radiation pattern factor and different response of same instrument after geometrical spreading. The method was used successfully for estimating site re- sponse of stations of Zhejiang Seismic Network and measuring Lg wave attenuation. The study is based on 20 earth- quakes occurred in northeast of Taiwan with magnitude MS5.0~6.7 and 960 seismic wave records from 16 stations in Zhejiang area from 2002 to 2005. The parameters of site response and Lg attenuation were calculated with a fre- quency interval of 0.2 Hz in the range of 0.5 Hz to 10 Hz. Lg wave attenuation coefficient corresponding to U-D, E-W and N-S components are γ ( f )=0.00175 f 0.43485, γ ( f )=0.00145f 0.48467 and γ ( f )=0.0021f 0.41241, respectively. It is found that the site response is component-independent. It is also found that the site response of QIY station is significant above the frequency of 1.5 Hz, and that the site response of NIB station is low for most frequency     

Scattering and anelastic attenuation of seismic energy in northern Greece

The relative contribution of scattering (Q _s ^–1 ) and intrinsic (Q _i ^–1 ) attenuation to the totalS-wave attenuation for the frequencies of 1.5, 3.0, 6.0 and 12.0 Hz has been studied by applying the radiative energy transfer theory, Data of local earthquakes which occurred in northern Greece and were recorded by the permanent telementered network of the Geophysical Laboratory of the University of Thessaloniki have been used. The results show that in this area the scattering attenuation is dominant over all frequencies while intrinsic attenuation is significantly lower. The estimatedQ _s ^–1 andQ _i ^–1 values have frequency dependences off ^–0.72 andf ^–0.45, respectively. The frequency dependence ofQ _s ^–1 is the same as that of the codaQ _c ^–1 , obtained by applying the single scattering model, which probably implies that the frequency dependence of the coda wave attenuation is attributed to the frequency dependence of the scattering attenuation.Q _c ^–1 values are very close to scattering attenuation for short lapse times, (10–20 sec), and intermediate between scattering and intrinsic attenuation for the longer lapse times, (50–100 sec). This difference is explained as the result of the depth-dependent attenuation properties and the multiple scattering effects.     

Duration of earthquake fault motion in California

Duration of high frequency (5–25 Hz) radiation of energy from earthquake sources in California is consistent with the estimates of fault length and with dislocation velocity estimates of 2–3 km/sec. This duration can be described by an exponential function of magnitude for 2·5 < M < 7·5. It is related to the times it takes the dislocation to spread over the fault width (1/f₂), and the fault length (～ 1/f₁), and to reach its ultimate amplitude (T₀). The results in this paper can be used to estimate the range of amplitudes and the duration of long period pulses of strong ground motion near faults, as these long period pulses can be related to the properties of high-frequency radiation from the source. Such pulses must be considered in the analyses of yielding structures, when the average peak acceleration of the pulse exceeds the yield resistance seismic coefficient of the structure.     

Internal soil moisture and piezometric responses to rainfall-induced shallow slope failures

Better knowledge regarding internal soil moisture and piezometric responses in the process of rainfall-induced shallow slope failures is the key to an effective prediction of the landslide and/or debris flow initiation. To this end, internal soil moisture and piezometric response of 0.7-m-deep, 1.5-m-wide, 1.7-m-high, and 3.94-m-long semi-infinite sandy slopes rested on a bi-linear impermeable bedrock were explored using a chute test facility with artificial rainfall applications. The internal response time defined by the inflection point of the soil moisture and piezometric response curves obtained along the soil–bedrock interface were closely related to some critical failure states, such as the slope toe failure and extensive slope failures. It was also found that the response times obtained at the point of abrupt bedrock slope decrease can be used as indicators for the initiation of rainfall-induced shallow slope failures. An investigation of spatial distributions of soil water content, ω (or degrees of saturation, S_r), in the slope at critical failure states shows that the 0.2 m – below – surface zone remains unsaturated with S_r 40–60%, regardless of their distances from the toe and the rainfall intensity. Non-uniform distributions of ω (or S_r) along the soil–bedrock interface at critical failure states were always associated with near-saturation states (S_r 80–100%) around the point of bedrock slope change or around the transient ‘toe’ upstream of the slumped mass induced by the retrogressive failure of the slope. These observations suggest the important role of the interflow along the soil–bedrock interface and the high soil water content (or high porewater pressure) around the point of bedrock slope deflection in the rainfall-induced failure of sandy slopes consisting of shallow impermeable bedrocks. The present study proposes an ‘internal response time’ criterion to substantiate the prediction of rainfall-induced shallow slope failures. It is believed that the ‘internal response time’ reflects the overall characteristics of a slope under rainfall infiltration and can be as useful as the conventional meteorology-based threshold times. The ‘internal response time’ theory can be generalized via numerical modeling of slope hydrology, slope geology and slope stability in the future.     

Long period microtremors, microseisms and earthquake damage: Northridge, CA, earthquake of 17 January 1994

Three studies of site amplification factors, based on the recorded aftershocks, and one study based on strong motion data, are compared one with another and with the observed distribution of damage from the Northridge, CA, earthquake of 17 January 1994 (M_L=6.4). In the epicentral area, when the peak ground velocities are larger than v_m≈15 cm/s, nonlinear response of soil begins to distort the amplification factors determined from small amplitude (linear) wave motion. Moving into the area of near-field and strong ground motion (v_m>30 cm/s), the site response becomes progressively more affected by the nonlinear soil response. Based on the published results, it is concluded that site amplification factors determined from small amplitude waves (aftershocks, small earthquakes, coda waves) and their transfer-function representation may be useful for small and distant earthquake motions, where soils and structures respond to earthquake waves in a linear manner. However in San Fernando Valley, during the Northridge earthquake, the observed distribution of damage did not correlate with site amplification determined from spectra of recorded weak motions. Mapping geographical distribution of site amplification using other than very strong motion data, therefore appears to be of little use for seismic hazard analyses.     

Estimation of frequency dependent coda wave attenuation structure at the vicinity of Cairo Metropolitan Area

Estimation of seismic wave attenuation in the shallow crust in terms of coda wave Q structure previously investigated in the vicinity of Cairo Metropolitan Area was improved using seismograms of local earthquakes recorded by the Egyptian National Seismic Network. The seismic wave attenuation was measured from the time decay of coda wave amplitudes on narrow bandpass filtered seismograms based on the single scattering theory. The frequency bands of interest are from 1.5 to 18 Hz. In general, the values obtained for various events recorded at El-Fayoum and Wadi Hagul stations are very similar for all frequency bands. A regional attenuation law Q _c = 85.66 f ^0.79 was obtained.     

Direct estimation of catchment response time parameters in medium to large catchments using observed streamflow data

In single‐event deterministic design flood estimation methods, estimates of the peak discharge are based on a single and representative catchment response time parameter. In small catchments, a simplified convolution process between a single‐observed hyetograph and hydrograph is generally used to estimate time parameters such as the time to peak (T_P), time of concentration (T_C), and lag time (T_L) to reflect the “observed” catchment response time. However, such simplification is neither practical nor applicable in medium to large heterogeneous catchments, where antecedent moisture from previous rainfall events and spatially non‐uniform rainfall hyetographs can result in multi‐peaked hydrographs. In addition, the paucity of rainfall data at sub‐daily timescales further limits the reliable estimation of catchment responses using observed hyetographs and hydrographs at these catchment scales. This paper presents the development of a new and consistent approach to estimate catchment response times, expressed as the time to peak (T_Px) obtained directly from observed streamflow data. The relationships between catchment response time parameters and conceptualised triangular‐shaped hydrograph approximations and linear catchment response functions are investigated in four climatologically regions of South Africa. Flood event characteristics using primary streamflow data from 74 flow‐gauging stations were extracted and analysed to derive unique relationships between peak discharge, baseflow, direct runoff, and catchment response time in terms of T_Px. The T_Px parameters are estimated from observed streamflow data using three different methods: (a) duration of total net rise of a multipeaked hydrograph, (b) triangular‐shaped direct runoff hydrograph approximations, and (c) linear catchment response functions. The results show that for design hydrology and for the derivation of empirical equations to estimate catchment response times in ungauged catchments, the catchment T_Px should be estimated from both the use of an average catchment T_Px value computed using either Methods (a) or (b) and a linear catchment response function as used in Method (c). The use of the different methods in combination is not only practical but is also objective and has consistent results.     

Comparison of amplitude decay rates in reflection,refraction, and local earthquake records

Three types of seismic data recorded near Coalinga, California were analyzed to study the behavior of scattered waves: 1) aftershocks of the May 2, 1983 earthquake, recorded on verticalcomponent seismometers deployed by the USGS; 2) regional refraction profiles using large explosive sources recorded on essentially the same arrays above; 3) three common-midpoint (CMP) reflection surveys recorded with vibrator sources over the same area. Records from each data set were bandpassed filtered into 5 Hz wide passbands (over the range of 1–25 Hz), corrected for geometric spreading, and fit with an exponential model of amplitude decay. Decay rates were expressed in terms of inverse codaQ (Q _c ^–1 ).Q _c ^–1 values for earthquake and refraction data are generally comparable and show a slight decrease with increasing frequency. Decay rates for different source types recorded on proximate receivers show similar results, with one notable exception. One set of aftershocks shows an increase ofQ _c ^–1 with frequency.Where the amplitude decay rates of surface and buried sources are similar, the coda decay results are consistent with other studies suggesting the importance of upper crustal scattering in the formation of coda. Differences in the variation ofQ _c ^–1 with frequency can be correlated with differences in geologic structure near the source region, as revealed by CMP-stacked reflection data. A more detailed assessment of effects such as the depth dependence of scattered contributions to the coda and the role of intrinsic attenuation requires precise control of source-receiver field geometry and the study of synthetic seismic data calculated for velocity models developed from CMP reflection data.     

CalculationoflongperiodresponsespectratoearthquakegroundmotionfromseismogramsofType513seismographs

ＣａｌｃｕｌａｔｉｏｎｏｆｌｏｎｇｐｅｒｉｏｄｒｅｓｐｏｎｓｅｓｐｅｃｔｒａｔｏｅａｒｔｈｑｕａｋｅｇｒｏｕｎｄｍｏｔｉｏｎｆｒｏｍｓｅｉｓｍｏｇｒａｍｓｏｆＴｙｐｅ５１３ｓｅｉｓｍｏｇｒａｐｈｓＹＡＮＸＩＡＮＧＹＵ（俞言祥）ＳＵＹＵＮＷＡＮ...     

Simulation of near source ground motion and its characteristics

It is fact that the severe ground motions of shear waves have a strong effect on the dynamic behavior of buildings and civil structures. We simulate near source strong motions of a pure shear wave and synthesize small motions, using the parameters based on the recorded accelerograms at the site that is regarded as a base rock in the Osaka basin, Japan. By making use of a stochastic technique, we can easily introduce higher frequency contents in the motions and apply the technique to the synthesis of small waves regarding as a green function. We also introduce to the analysis the useful relationships among the time duration T_d, the seismic moment M₀, the corner frequency f_c and the high cutoff frequency f_max which were regressed by a simple representation scheme. Considering two active faults that may affect severe damage on buildings and civil structures, we try to predict strong ground motions in Osaka basin and show the characteristics of them.     

International standard (ISO) on seismic actions for designing geotechnical works–An overview

The seismic performance of geotechnical works is significantly affected by ground displacement. In particular, soil–structure interaction and effects of liquefaction play major roles and pose difficult problems for engineers. An International Standard, ISO23469, is being developed for addressing these issues in a systematic manner within a consistent framework. The objective of this paper is to give an overview of this International Standard.In this International Standard, the seismic actions are determined through two stages. The first stage determines basic seismic action variables, including the earthquake ground motion at the site, the potential for earthquake-associated phenomena such as liquefaction and induced lateral ground displacement. These basic variables are used, in the second stage, for specifying the seismic actions for designing geotechnical works. In the second stage, the soil–structure interaction plays a major role. Types of analyses are classified based on a combination of static/dynamic analyses and the procedure for soil–structure interaction classified as follows:

– simplified: soil–structure interaction of a global system is modeled as an action on a substructure;

– detailed: soil–structure interaction of a global system is modeled as a coupled system.

Earthquake damage detection in the Imperial County Services Building I: The data and time–frequency analysis

The former Imperial County Services Building was a six-story reinforced concrete structure in the El Centro, California, severely damaged by the 1979 Imperial Valley earthquake. It represents a rare case of an instrumented building that has been damaged, and thus can serve as a full-scale benchmark to evaluate and further develop structural health monitoring methods. This paper presents an analysis of inter-story drifts, and of changes in the first NS and EW system frequencies (estimated from the ridge of the Gabor transform of the relative roof displacement response) as indicators of the occurrence of damage. The detected initial decreases of system frequency, of about 28% for NS and 24% for EW motions, are not believed to be due to severe damage. The subsequent decreases, of about 44% for NS and 21% for EW motions, are attributed to damage. Near the end of shaking, increases of about 35% for the NS and 36% for EW motions were detected. (These percentage changes were computed with respect to the previous value for particular interval, rather than a fixed reference). During the most severe shaking, the inter-story drifts exceeded 0.5% for NS and 1.5% for EW motions.     

An attenuation study using earthquakes from the 1997 Umbria-Marche sequence

We studied spatial and temporal characteristics of seismic attenuation inCentral Italy using S- and coda- waves recorded by the MarchesanSeismograph Network from earthquakes located in the epicentral area ofthe 1997 Umbria-Marche sequence. The amplitude decay of the S waveswith distance was defined calculating empirical attenuation functions at 15frequencies between 1 and 25 Hz. We analyzed separately foreshocks andaftershocks and we found the same attenuation functions, suggesting thatthe possible temporal variations could be confined in a small area. Thefrequency dependence of Q _S was approximated by the equation Q _S=18 · f ^2.0between 1 and 10 Hz. At higher frequencies (10–25 Hz), the frequencydependence of Q _s weakens, having an average value of Q _S=990. We also estimated Q from coda waves (Q _C) using the single-scattering models of Aki andChouet (1975) and Sato (1977). We found that Q _C=77 · f ^0.6, (between 2 and 20Hz) at the western side of the mountain chain, using either foreshocks oraftershocks. This relation is consistent with previous estimates of Q _Creported for the Central Apennines. For a volume sampling the Colfioritobasin, the Apennines and the Marche region we found that Q _C=55 · f ^0.8,indicating highattenuation below the mountain belt. To detect small temporal changes ofQ, we calculated spectral ratios of 5 temporal doublets located in theepicentral area and recorded at the closest station. We found temporalchanges of Q that vary from 27% to 56%, depending on the locationof the doublets. This variability suggests that the temporal change ofattenuation may depend on the spatial variation of Q and perhaps on thespatial distribution of tectonic stress in the epicentral area.     

Source Parameters for 1999 North Anatolian Fault Zone Aftershocks

We develop a data set of aftershock recordings of the 1999, M = 7.4 Izmit and M = 7.2 Duzce (Turkey) earthquakes to study their source parameters. We combined seismograms from 44 stations maintained by several sources (organizations) to obtain a unified data set of events (2.1 ≤ M_w ≤ 5.5). We calculate source parameters of these small earthquakes by two methods that use different techniques to address the difficulty in obtaining source spectra for small earthquakes subject to interference from site response. One method (program NetMoment (NM), Hutchings, 2004) uses spectra of direct S waves in a simultaneous inversion of local high-frequency network data to estimate seismic moment, source corner frequency (f_c), site attenuation (k) and whole-path Q. This approach takes advantage of the source commonality in all recordings for a particular earthquake by fitting a common Brune source spectrum to the data with a and individual k. The second approach (Mayeda et al., 2003) uses the coda method (CM) to obtain “nonmodel-based” source spectra and moment estimates from selected broadband recording sites. We found that both methods do well for events that allow the comparison with seismic moment estimates derived from waveform modeling. Also, source spectra obtained from the two methods are very closely matched for most of the events they have in common. We use an F test to examine the trade-off between k and f_c picks identified by the direct S-wave method. About half of the events could be constrained to have less than a 50% average uncertainty in f_c and k. We used these source spectra solutions to calculate energy and apparent stress and compare these to estimates from the selected “good quality” source spectra from CM. Both studies have values mutually consistent and show a similar increase in apparent stress with increasing moment. This result has added merit due to the independent approaches to calculate apparent stress. We conclude that both methods are at least partially validated by our study, and they both have usefulness for different circumstances of recording local small earthquakes. CM would work well in studies for which there is a broad magnitude range of events and NM works well for local events recorded by band-limited recorders.