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

Sensor placement for the analysis of seismic surface waves: sources of error, design criterion and array design algorithms

Seismic surface waves can be measured by deploying an array of seismometers on the surface of the earth. The goal of such measurement surveys is, usually, to estimate the velocity of propagation and the direction of arrival of the seismic waves. In this paper, we address the issue of sensor placement for the analysis of seismic surface waves from ambient vibration wavefields. First, we explain in detail how the array geometry affects the mean-squared estimation error of parameters of interest, such as the velocity and direction of propagation, both at low and high signal-to-noise ratios (SNRs). Secondly, we propose a cost function suitable for the design of the array geometry with particular focus on the estimation of the wavenumber of both Love and Rayleigh waves. Thirdly, we present and compare several computational approaches to minimize the proposed cost function. Numerical experiments verify the effectiveness of our cost function and resulting array geometry designs, leading to greatly improved estimation performance in comparison to arbitrary array geometries, both at low and high SNR level

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

A Simplified Classification of the Relative Tsunami Potential in Swiss Perialpine Lakes Caused by Subaqueous and Subaerial Mass-Movements

Historical reports and recent studies have shown that tsunamis can also occur in lakes where they may cause large damages and casualties. Among the historical reports are many tsunamis in Swiss lakes that have been triggered both by subaerial and subaqueous mass movements (SAEMM and SAQMM). In this study, we present a simplified classification of lakes with respect to their relative tsunami potential. The classification uses basic topographic, bathymetric, and seismologic input parameters to assess the relative tsunami potential on the 28 Swiss alpine and perialpine lakes with a surface area >1km2. The investigated lakes are located in the three main regions “Alps,” “Swiss Plateau,” and “Jura Mountains.” The input parameters are normalized by their range and a k-means algorithm is used to classify the lakes according to their main expected tsunami source. Results indicate that lakes located within the Alps show generally a higher potential for SAEMM and SAQMM, due to the often steep surrounding rock-walls, and the fjord-type topography of the lake basins with a high amount of lateral slopes with inclinations favoring instabilities. In contrast, the missing steep walls surrounding lakeshores of the “Swiss Plateau” and “Jura Mountains” lakes result in a lower potential for SAEMM but favor inundation caused by potential tsunamis in these lakes. The results of this study may serve as a starting point for more detailed investigations, considering field data

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

Seismic site characterization of the Kastelli (Kissamos) Basin in northwest Crete (Greece): Assessments using ambient noise recordings

Crete is actively seismic and site response studies are needed for estimating local site conditions subjected to seismic activity. In order to collect basic data, we performed ambient noise recordings to estimate the site response of the surface and near subsurface structure of the small-scale Kastelli Basin in northwest Crete. The spatial horizontal to vertical spectral ratios (HVSR) resonance pattern of the investigated sites in the centre of the Basin consists of either one or two peaks divided into low to high frequency range in different sites as follows: (a) in some sites only one amplified peak at low frequencies (0.6–1.2 Hz), (b) in other sites only one amplified peak at medium frequencies (2.9–8.5 Hz) and (c) in yet other sites two amplified peaks in the low to high frequency range (0.6–15.5 Hz). The investigated sites are amplified in the frequency range 0.6–15.5 Hz, while the amplitude reaches to a factor of 4 in the spectral ratios. The one HVSR amplified peak at low frequencies is related to locally soft or thick Quaternary deposits. Microtremors were measured in the coastal northwest part of the Basin in a well—lithified Cretaceous limestone site characterized by fractures and faults striking predominantly in a sector NNE to NNW. Sites of one amplified peak at medium frequencies are extended from coastal northwest to southwest delineating a structure striking to NNW. The two amplified peaks are attributed to shallow subsurface heterogeneities/irregularities, locally induced by fault zones and to the overlying Quaternary deposits. Spatial HVSR variations in the frequency and HVSR shape delineate four structures striking NNE, NNW and in a sector NW to WNW, crosscutting the dense populated Basin suggesting that microtremors could be a valuable tool for providing a first approximation of fault zone delineation at least for the Kastelli-Kissamos Basin. The Basin is classified into the X soil category of the Greek Seismic Code 2000.This work was implemented through the project entitled “Interdisciplinary Multi-Scale Research of Earth-quake Physics and Seismotectonics at the Front of the Hellenic Arc (IMPACT-ARC)” in the framework of action “ARCHIMEDES III—Support of Research Teams at TEI of Crete” (MIS380353) of the Operational Program “Education and Lifelong Learning” and is co-financed by the European Union (European Social Fund) and Greek national fund

Application of Surface wave methods for seismic site characterization

Surface-wave dispersion analysis is widely used in geophysics to infer a shear wave velocity model of the subsoil for a wide variety of applications. A shear-wave velocity model is obtained from the solution of an inverse problem based on the surface wave dispersive propagation in vertically heterogeneous media. The analysis can be based either on active source measurements or on seismic noise recordings. This paper discusses the most typical choices for collection and interpretation of experimental data, providing a state of the art on the different steps involved in surface wave surveys. In particular, the different strategies for processing experimental data and to solve the inverse problem are presented, along with their advantages and disadvantages. Also, some issues related to the characteristics of passive surface wave data and their use in H/V spectral ratio technique are discussed as additional information to be used independently or in conjunction with dispersion analysis. Finally, some recommendations for the use of surface wave methods are presented, while also outlining future trends in the research of this topic

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