Rayleigh wave ellipticity as a function of frequency is closely linked to underground structure, i.e., shear wave velocity profile and sediment thickness. The possibility to calculate these ...underground properties by inverting ellipticity curves has recently been shown. We propose a new technique enabling the Rayleigh wave ellipticity to be recovered over a wide frequency range by using ambient noise recordings. Based on the random decrement technique commonly used to characterize dynamic parameters of buildings, this method eliminates all wave types except Rayleigh waves. We apply the method to noise synthetics simulated for different underground structures and show its applicability to real seismic noise data.
The single-station microtremor horizontal-to-vertical spectral ratio (MHVSR) method was initially proposed to retrieve the site amplification function and its resonance frequencies produced by ...unconsolidated sediments overlying high-velocity bedrock. Presently, MHVSR measurements are predominantly conducted to obtain an estimate of the fundamental site frequency at sites where a strong subsurface impedance contrast exists. Of the earthquake site characterization methods presented in this special issue, the MHVSR method is the furthest behind in terms of consensus towards standardized guidelines and commercial use. The greatest challenges to an international standardization of MHVSR acquisition and analysis are (1) the
what
— the underlying composition of the microtremor wavefield is site-dependent, and thus, the appropriate theoretical (forward) model for inversion is still debated; and (2) the
how
— many factors and options are involved in the data acquisition, processing, and interpretation stages. This paper reviews briefly a historical development of the MHVSR technique and the physical basis of an MHVSR (the
what
). We then summarize recommendations for MHVSR acquisition and analysis (the
how
). Specific sections address MHVSR interpretation and uncertainty assessment.
The Avignonet landslide affects a 2 by 2 km area covered by clayey deposits. This paper presents the use of the seismic ambient noise cross‐correlation technique to retrieve a 3‐D model of the shear ...wave velocity of the area. Seismic ambient noise was recorded during 15 days at 13 stations located on the landslide. Cross correlations computed between the vertical components of all station pairs allow the retrieval of the Rayleigh wave Green's functions and the estimation of their group velocity dispersion curves in the 1.7–5 Hz frequency range. At frequencies lower than 1.5 Hz, the anisotropy of the wavefield strongly influences the apparent Rayleigh wave velocities. Moreover, the analysis of the convergence of the cross correlations shows that at frequencies higher than 5 Hz, the recording time length was not sufficient for the cross correlation to be stable. These 1.7–5 Hz passive group dispersion curves are complementary to the ones computed from shot signals in the 3–7 Hz frequency range. A tomographic inversion of the resulting 1.7–7 Hz Rayleigh wave group dispersion curves provides local group dispersion curves at each cell of the tomographic grid. These are inverted with a neighborhood algorithm to retrieve the 3‐D model of the landslide. Despite the complex wave propagation in the eastern part of the landslide and the sparse ray coverage, estimated velocities and first‐order features are in good agreement with previous investigations.
The knowledge of the local soil structure is important for the assessment of seismic hazards. A widespread, but time-consuming technique to retrieve the parameters of the local underground is the ...drilling of boreholes. Another way to obtain the shear wave velocity profile at a given location is the inversion of surface wave dispersion curves. To ensure a good resolution for both superficial and deeper layers, the used dispersion curves need to cover a wide frequency range. This wide frequency range can be obtained using several arrays of seismic sensors or a single array comprising a large number of sensors. Consequently, these measurements are time-consuming. A simpler alternative is provided by the use of the ellipticity of Rayleigh waves. The frequency dependence of the ellipticity is tightly linked to the shear wave velocity profile. Furthermore, it can be measured using a single seismic sensor. As soil structures obtained by scaling of a given model exhibit the same ellipticity curve, any inversion of the ellipticity curve alone will be ambiguous. Therefore, additional measurements which fix the absolute value of the shear wave velocity profile at some points have to be included in the inversion process. Small-scale spatial autocorrelation measurements or MASW measurements can provide the needed data. Using a theoretical soil structure, we show which parts of the ellipticity curve have to be included in the inversion process to get a reliable result and which parts can be omitted. Furthermore, the use of autocorrelation or high-frequency dispersion curves will be highlighted. The resulting guidelines for inversions including ellipticity data are then applied to real data measurements collected at 14 different sites during the European NERIES project. It is found that the results are in good agreement with dispersion curve measurements. Furthermore, the method can help in identifying the mode of Rayleigh waves in dispersion curve measurements.
The collective excitation of city structures by a seismic wavefield and the subsequent multiple Structure-Soil-Structure Interactions (SSSIs) between the buildings are usually disregarded in ...conventional seismology and earthquake engineering practice. The objective here is to qualify and quantify these complex multiple SSSIs through the design of an elementary study case, which serves as a benchmark for theoretical, numerical and experimental crossed-analysis. The experimental specimen consists of an idealized site-city setup with up to 37 anisotropic resonant structures arranged at the top surface of an elastic layer and in co-resonance with it. The experimental data from shaking table measurements is compared with the theoretical and numerical results provided respectively by an equivalent city-impedance model derived analytically from homogenization in the long-wavelength approximation and a model based on boundary elements. The signatures of the site-city interactions are identified in the frequency, time and space domain, and in particular consist of a frequency-dependent free/rigid switch in the surface condition at the city resonance, beatings in the records and the depolarization of the wavefield. A parametric study on the city density shows that multiple SSSIs among the city structures (five are sufficient) can have significant effects on both the seismic response of its implantation site and that of the buildings. Key parameters are provided to assess site-city interactions in the low seismic frequency range: They involve the mass and rigidity of the city compared to those of the soil and the damping of the building.
Site-specific seismic hazard studies involving detailed account of the site response require the prior estimate of the hazard at the local reference bedrock level. As the real characteristics of such ...local bedrock often correspond to “hard-rock” with S-wave velocity exceeding 1.5 km/s, “standard rock” PSHA estimates should be adjusted in order to replace the effects of “standard-rock” characteristics by those corresponding to the local bedrock. The current practice involves the computation of scaling factors determined on the basis of V
S
(S-wave velocity) and “κ
0
” (site specific, high-frequency attenuation parameter) values, and generally predicts larger high-frequency motion on hard rock compared to standard rock. However, it also proves to be affected by large uncertainties (Biro and Renault, Proceedings of the 15th world conference on earthquake engineering, 24–28,
2012
; Al Atik et al., Bull Seism Soc Am 104(1):336–346
2014
), mainly attributed to (i) the measurement of host and target parameters, and (ii) the forward and inverse conversions from the response spectrum domain to the Fourier domain to apply the V
S
and κ
0
adjustments. Moreover, recent studies (Ktenidou and Abrahamson, Seismol Res Lett 87(6):1465–1478,
2016
) question the appropriateness of current V
S
− κ
0
scaling factors, so that the significant amplification of high frequency content for hard-rock with respect to standard-rock seems overestimated. This paper discusses the key aspects of a few, recently proposed, alternatives to the standard approach. The calibration of GMPEs directly in the Fourier domain rather than in the response spectrum domain is one possibility (Bora et al., Bull Seism Soc Am 105(4):2192–2218,
2015
, Bull Earthq Eng 15(11):4531–4561,
2017
). Another possibility is the derivation of GMPEs which be valid also for hard-rock conditions (e.g. Laurendeau et al., Bull Earthq Eng 16(6):2253–2284,
2018
). In this latter case the host site response is first removed using theoretical site response analyses (and site velocity profile), or generalized inversions techniques. A third possibility is to use existing hard rock surface recordings to derive purely empirical scaling models from standard rock to hard rock (Ktenidou et al., PEER Report, Pacific Earthquake Engineering Research Center, Berkeley,
2016
). Finally, when a sufficient amount of records are available at a given site, generic GMPEs can be scaled to the site-specific ground motion using empirical site residual (δ
S2Ss
) (Kotha et al., Earthq Spectra 33(4):1433–1453,
2017
; Ktenidou et al., Bulletin of Earthquake Engineering 16(6):2311–2336,
2018
). Such alternative approaches present the advantage of a significant simplification with respect to the current practice (with thus a reduced number of uncertainty sources); their generalization calls however for high-quality recordings (including high-quality site metadata) for both host regions and target sites, especially for small to moderate magnitude events. Our answer to the question in the title is thus “No, alternative approaches exist and are promising; though, their routine implementation requires additional work regarding systematic site characterization (for the host regions) and high-quality site characterization/instrumentation (for the target site), and so do also the needed improvements of the existing HTTA procedure”.
Elastic fundamental frequency is a key-parameter of simplified seismic design and vulnerability assessment methods. Empirical relationships exist in codes to estimate this frequency but they miss ...experimental data to validate them accounting for national feature of building design and, above all, corresponding uncertainties. Even if resonance frequency extracted from ambient vibrations may be larger than the elastic frequency (at yield) generally used in earthquake engineering, ambient vibration recordings may provide a large set of data for statistical analysis of periods versus building characteristics relationships. We recorded ambient vibrations and estimated the fundamental frequency of about 60 buildings of various types (RC and masonry) in Grenoble City (France). These data complete the set existing yet, made of 26 RC-buildings of Grenoble (Farsi and Bard 2004) and 28 buildings in Nice (France) (Dunand 2005). Statistical analysis of these experimental data was performed for fundamental frequencies of RC shear wall structures and the results are compared with existing relationships. Only building height or number of stories has a statistical relevancy to estimate the resonance frequency but the variability associated to the proposed relationships is large. Moreover, we compared the elastic part of capacity curves of RC and masonry buildings used in the European Risk-UE method for vulnerability assessment with the experimental frequencies. The variability is also large and the curves may not be consistent with French existing buildings.
The InterPACIFIC project was aimed at assessing the reliability, resolution, and variability of geophysical methods in estimating the shear-wave velocity profile for seismic ground response analyses. ...Three different subsoil conditions, which can be broadly defined as soft-soil, stiff-soil, and hard-rock, were investigated. At each site, several participants performed and interpreted invasive measurements of shear wave velocity (Vs) and compression wave velocity (Vp) in the same boreholes. Additionally, participants in the project analysed a common surface-wave dataset using their preferred strategies for processing and inversion to obtain Vs profiles. The most significant difference between the invasive borehole methods and non-invasive surface wave methods is related to resolution of thin layers and abrupt contrasts, which is inherently better for invasive methods. However, similar variability is observed in the estimated invasive and non-invasive Vs profiles, underscoring the need to account for such uncertainty in site response studies. VS,30 estimates are comparable between invasive and non-invasive methods, confirming that the higher resolution provided by invasive methods is quite irrelevant for computing this parameter.
•Reliability of Vs models from geophysical tests is studied with a blind test.•Three subsoil conditions are considered: soft soil, stiff soil, rock outcrop.•In-hole methods and surface wave methods are considered.•In-hole methods provide high resolution and accuracy, but the precision is similar.There is a good match between estimates of VS,30, with comparable variability.
A key component in seismic hazard assessment is the estimation of ground motion for hard rock sites, either for applications to installations built on this site category, or as an input motion for ...site response computation. Empirical ground motion prediction equations (GMPEs) are the traditional basis for estimating ground motion while V
S30
is the basis to account for site conditions. As current GMPEs are poorly constrained for V
S30
larger than 1000 m/s, the presently used approach for estimating hazard on hard rock sites consists of “host-to-target” adjustment techniques based on V
S30
and κ
0
values. The present study investigates alternative methods on the basis of a KiK-net dataset corresponding to stiff and rocky sites with 500 < V
S30
< 1350 m/s. The existence of sensor pairs (one at the surface and one in depth) and the availability of P- and S-wave velocity profiles allow deriving two “virtual” datasets associated to outcropping hard rock sites with V
S
in the range 1000, 3000 m/s with two independent corrections: 1/down-hole recordings modified from within motion to outcropping motion with a depth correction factor, 2/surface recordings deconvolved from their specific site response derived through 1D simulation. GMPEs with simple functional forms are then developed, including a V
S30
site term. They lead to consistent and robust hard-rock motion estimates, which prove to be significantly lower than host-to-target adjustment predictions. The difference can reach a factor up to 3–4 beyond 5 Hz for very hard-rock, but decreases for decreasing frequency until vanishing below 2 Hz.