A combined inversion of Rayleigh wave dispersion and 2‐D resonance frequencies

Shear wave velocities of the sediment fill of a deep Alpine valley are estimated from ambient noise recorded on linear and circular arrays. We propose a combined inversion of 2‐D resonance frequencies identified from site‐to‐reference spectral ratios and Rayleigh wave dispersion curves obtained from frequency‐wavenumber analysis. The method is tested on synthetic noise data and on noise recorded at three sites in the Rhône valley in Southern Switzerland. Previous studies have shown that 2‐D resonance dominates the ambient vibration wavefield at low frequencies at the investigated sites. Inversion techniques which assume that the noise wavefield consists mainly of horizontally propagating surface waves will, therefore, fail to resolve shear wave velocities at depths below around 500 m. We show that standard techniques lead to an overestimation of shear wave velocities at depth when applied to synthetic and observed ambient noise. The combined inversion is able to resolve the shear wave velocities in the initial velocity model when applied to synthetic noise records. Application of the method to observed ambient noise improves resolution at depth and yields realistic shear wave velocities for the lower part of the sediment fil

Time-Frequency-Wavenumber Analysis of Surface Waves Using the Continuous Wavelet Transform

A modified approach to surface wave dispersion analysis using active sources is proposed. The method is based on continuous recordings, and uses the continuous wavelet transform to analyze the phase velocity dispersion of surface waves. This gives the possibility to accurately localize the phase information in time, and to isolate the most significant contribution of the surface waves. To extract the dispersion information, then, a hybrid technique is applied to the narrowband filtered seismic recordings. The technique combines the flexibility of the slant stack method in identifying waves that propagate in space and time, with the resolution of f-k approaches. This is particularly beneficial for higher mode identification in cases of high noise levels. To process the continuous wavelet transform, a new mother wavelet is presented and compared to the classical and widely used Morlet type. The proposed wavelet is obtained from a raised-cosine envelope function (Hanning type). The proposed approach is particularly suitable when using continuous recordings (e.g., from seismological-like equipment) since it does not require any hardware-based source triggering. This can be subsequently done with the proposed method. Estimation of the surface wave phase delay is performed in the frequency domain by means of a covariance matrix averaging procedure over successive wave field excitations. Thus, no record stacking is necessary in the time domain and a large number of consecutive shots can be used. This leads to a certain simplification of the field procedures. To demonstrate the effectiveness of the method, we tested it on synthetics as well on real field data. For the real case we also combine dispersion curves from ambient vibrations and active measurement

Two-dimensional resonances in Alpine valleys identified from ambient vibration wavefields

Although numerical simulations have for long shown the importance of 2-D resonances in site effect estimations of sediment-filled valleys, this phenomenon is usually not taken into account by current hazard assessment techniques. We present an approach to identify the resonance behaviour of a typical Alpine valley by analysis of ambient noise recorded simultaneously on a dense array. The applicability of the method is evaluated further using synthetic ambient noise acquired with current 3-D numerical simulation techniques. Resonance frequencies of the fundamental mode SV and the fundamental and first higher mode of SH are identified from measured data with the reference station method, verifying results of previous studies. Patterns of spectral amplitude and phase behaviour obtained from observed and synthetic noise correlate well with properties expected for 2-D resonance. Application of a frequency-wavenumber technique shows that the noise wavefield is dominated by standing waves at low frequencies (0.25 to 0.50 Hz). The different 2-D resonance modes are creating prominent peaks in horizontal-to-vertical spectral ratios, which can not be interpreted in terms of 1-D resonance. We conclude that ambient noise records measured simultaneously on a linear array perpendicular to the valley axis may be used for identification of resonance modes in sediment-filled valley

High-frequency ground motion amplification during the 2011 Tohoku earthquake explained by soil dilatancy

Ground motions of the 2011 Tohoku earthquake recorded at Onahama port (Iwaki, Fukushima prefecture) rank among the highest accelerations ever observed, with the peak amplitude of the 3-D acceleration vector approaching 2g. The response of the site was distinctively non-linear, as indicated by the presence of horizontal acceleration spikes which have been linked to cyclic mobility during similar observations. Compared to records of weak ground motions, the response of the site during the Mw 9.1 earthquake was characterized by increased amplification at frequencies above 10 Hz and in peak ground acceleration. This behaviour contrasts with the more common non-linear response encountered at non-liquefiable sites, which results in deamplification at higher frequencies. We simulate propagation of SH waves through the dense sand deposit using a non-linear finite difference code that is capable of modelling the development of excess pore water pressure. Dynamic soil parameters are calibrated using a direct search method that minimizes the difference between observed and simulated acceleration envelopes and response spectra. The finite difference simulations yield surface acceleration time-series that are consistent with the observations in shape and amplitude, pointing towards soil dilatancy as a likely explanation for the high-frequency pulses recorded at Onahama port. The simulations also suggest that the occurrence of high-frequency spikes coincided with a rapid increase in pore water pressure in the upper part of the sand deposit between 145 and 170 s. This sudden increase is possibly linked to a burst of high-frequency energy from a large slip patch below the Iwaki regio

Classifying seismic waveforms from scratch: a case study in the alpine environment

Nowadays, an increasing amount of seismic data is collected by daily observatory routines. The basic step for successfully analyzing those data is the correct detection of various event types. However, the visually scanning process is a time-consuming task. Applying standard techniques for detection like the STA/LTA trigger still requires the manual control for classification. Here, we present a useful alternative. The incoming data stream is scanned automatically for events of interest. A stochastic classifier, called hidden Markov model, is learned for each class of interest enabling the recognition of highly variable waveforms. In contrast to other automatic techniques as neural networks or support vector machines the algorithm allows to start the classification from scratch as soon as interesting events are identified. Neither the tedious process of collecting training samples nor a time-consuming configuration of the classifier is required. An approach originally introduced for the volcanic task force action allows to learn classifier properties from a single waveform example and some hours of background recording. Besides a reduction of required workload this also enables to detect very rare events. Especially the latter feature provides a milestone point for the use of seismic devices in alpine warning systems. Furthermore, the system offers the opportunity to flag new signal classes that have not been defined before. We demonstrate the application of the classification system using a data set from the Swiss Seismological Survey achieving very high recognition rates. In detail we document all refinements of the classifier providing a step-by-step guide for the fast set up of a well-working classification syste

A comparison of observed and simulated site response in the Rhône valley

Site effects in the city of Sion in the Rhône valley are analysed from weak motion signals recorded on a dense temporary array. We simulate the recorded events with a 3-D finite difference method for frequencies up to 4 Hz using a recently developed velocity model of the Sion basin. Site-to-reference Fourier spectral ratios are computed from 16 local and regional events. All sites exhibit amplification factors of up to 12 between 0.5 and 0.6 Hz, which can be reproduced by the numerical simulations. By rotating the weak motion to directions parallel and perpendicular to the valley axis, we show that this low-frequency amplification is caused by the SH00 and SV0 fundamental modes of 2-D resonance. Additional peaks of amplification can be observed at higher frequencies, with amplification factors of up to 20 at some sites. Application of the high-resolution frequency-wavenumber and the multiple signal characterization method to the vertical component of recorded and simulated signals show that edge-generated surface waves arriving from almost all directions dominate the wavefield at 1.25 and 2.50 Hz. Peak ground velocities computed from the simulated ground motion show interference patterns that depend strongly on the incidence direction, and the computed amplification of peak ground velocities are generally in agreement with the observations. We conclude that the complex 3-D geometry of the basin needs to be considered to evaluate site effects up to at least 2.5 H

Estimation of non-linear site response in a deep Alpine valley

We simulate non-linear behaviour of soils during strong ground motion in the Rhône valley in southern Switzerland. Previous studies of the site response using weak ground motion, ambient noise and linear 3-D FD simulations suggest that the 2-D structure of the basin will lead to amplification factors of up to 12 in the frequency band between 0.5 and 10 Hz. To estimate the importance of non-linear soil behaviour during strong ground motion in the Rhône valley we simulate the response of a superficial soft layer with a fully non-linear 1-D finite difference code. The non-linear wave propagator is based on an effective stress constitutive soil model capable of predicting pore pressure evolution due to shear. We determine the required dilatancy parameters from laboratory analysis of soil samples using cyclic triaxial tests. In order to include the effect of the strong 2-D structure in our non-linear analysis synthetic seismograms are convolved with the transfer function of the basin and then propagated through a 1-D non-linear layer. We find that reduced amplification due to soil non-linearity can be expected at rock accelerations above 0.5 ms−2, and that de-amplification occurs at ground motion levels of approximately 2 ms−2. Nevertheless, the spectral accelerations simulated for the valley centre are still exceeding the design spectra at about 0.5 Hz for magnitudes above 6.0, which reflects the strong amplification of ground motion by the deep 2-D resonance of the basin. For frequencies above 1 Hz the design spectra are generally in agreement with the strongest simulated accelerations. We evaluate the occurrence of soil failure using the 5 per cent strain criterion as a function of hypocentral distance and magnitude. Results confirm observations of liquefaction reported after the 1855 Mw 6.4 earthquake of Visp, and they suggest that soil liquefaction may occur at distances beyond those predicted by empirical relations in the valley. Near the basin edge, however, the simulated liquefaction occurrence agrees with the empirical relations. These results suggest that the response of the whole structure needs to be simulated in order to estimate the non-linear seismic response of complex basins like the Rhône valle

Empirical evaluation of microtremor H/V spectral ratio

The objective of this work is to perform a purely empirical assessment of the actual capabilities of the horizontal-to-vertical (H/V) spectral ratio technique to provide reliable and relevant information concerning site conditions and/or site amplification. This objective has been tackled through the homogeneous (re)processing of a large volume of earthquakes and ambient noise data recorded by different research teams in more than 200 sites located mainly in Europe, but also in the Caribbean and in Tehran. The original recordings were first gathered in a specific database with information on both the sites and recorded events. Then, for all sites close to an instrumented reference, average site-to-reference spectral ratios (“spectral ratio method” (SSR)) were derived in a homogeneous way (window selection, smoothing, signal-to-noise ratio threshold, averaging), as well as H/V ratios (“HVSRE–RF”) on earthquake recordings. H/V ratios were also obtained from noise recordings at each site (either specific measurements, or extracted from pre- or post-event noise windows). The spectral curves resulting from these three techniques were estimated reliable for a subset of 104 sites, and were thus compared in terms of fundamental frequency, amplitude and amplification bandwidth, exhibiting agreements and disagreements, for which interpretations are looked for in relation with characteristics of site conditions. The first important result consists in the very good agreement between fundamental frequencies obtained with either technique, observed for 81% of the analyzed sites. A significant part of the disagreements correspond to thick, low frequency, continental sites where natural noise level is often very low and H/V noise ratios do not exhibit any clear peak. The second important result is the absence of correlation between H/V peak amplitude and the actual site amplification measured on site-to-reference spectral ratios. There are, however, two statistically significant results about the amplitude of the H/V curve: the peak amplitude may be considered as a lower bound estimate of the actual amplification indicated by SSR (it is smaller for 79% of the 104 investigated sites), and, from another point of view, the difference in amplitude exhibits a questioning correlation with the geometrical characteristics of the sediment/basement interface: large SSR/HV differences might thus help to detect the existence of significant 2D or 3D effects.Published75-1084.1. Metodologie sismologiche per l'ingegneria sismicaJCR Journalreserve

A web application prototype for the multiscale modelling of seismic input

A web application prototype is described, aimed at the generation of synthetic seismograms for user-defined earthquake models. The web application graphical user interface hides the complexity of the underlying computational engine, which is the outcome of the continuous evolution of sophisticated computer codes, some of which saw the light back in the middle '80s. With the web application, even the non-experts can produce ground shaking scenarios at the local or regional scale in very short times, depending on the complexity of the adopted source and medium models, without the need of a deep knowledge of the physics of the earthquake phenomenon. Actually, it may even allow neophytes to get some basic education in the field of seismology and seismic engineering, due to the simplified intuitive experimental approach to the matter. One of the most powerful features made available to the users is indeed the capability of executing quick parametric tests in near real-time, to explore the relations between each model's parameter and the resulting ground motion scenario. The synthetic seismograms generated through the web application can be used by civil engineers for the design of new seismo-resistant structures, or to analyse the performance of the existing ones under seismic load.Comment: 23 pages, 13 figure

Erratum to: Seismic waves in 3-D: from mantle asymmetries to reliable seismic hazard assessment

Acknowledgments This paper is strongly relying upon the work of many colleagues and collaborators to whom we are greatly thankful, in particular: Carlo Doglioni, Gillian Foulger, Vahid Gholami, Hossein Hamzehloo, Volodya Kossobokov, Cristina La Mura, Anatoly Levshin, Andrea Magrin, Antonella Peresan, Federica Riguzzi, Franco Vaccari, Peter Varga, Tatiana Yanovskaya. Financial support from PRIN 2010-2011 and RITMARE projects, funded by Italian Ministry of University and Research, is gratefully acknowledged. In addition, Fig. 3 should be updated by new one as below