Refine
Year of publication
Document Type
- Article (48)
- Other (7)
- Postprint (2)
- Conference Proceeding (1)
- Doctoral Thesis (1)
- Review (1)
Keywords
- Erdbeben (5)
- Seismologie (4)
- Seismology (4)
- Site effects (4)
- Arrayseismologie (3)
- Earthquake (3)
- Surface waves and free oscillations (3)
- Tsunami (3)
- Wave propagation (3)
- array seismology (3)
Institute
Volcanic eruptions are often preceded by seismic activity that can be used to quantify the volcanic activity. In order to allow consistent inference of the volcanic activity state from the observed seismicity patterns, objective and time-invariant classification results achievable by automatic systems should be preferred. Most automatic classification approaches need a large preclassified data set for training the system. However, in case of a volcanic crisis, we are often confronted with a lack of training data due to insufficient prior observations. In the worst case (e. g., volcanic crisis related reconfiguration of stations), there are even no prior observations available. Finally, due to the imminent crisis there might be no time for the time-consuming process of preparing a training data set. For this reason, we have developed a novel seismic-event spotting technique in order to be less dependent on previously acquired data bases and classification schemes. We are using a learning-while-recording approach based on a minimum number of reference waveforms, thus allowing for the build-up of a classification scheme as early as interesting events have been identified. First, short-term wave-field parameters (here, polarization and spectral attributes) are extracted from a continuous seismic data stream. The sequence of multidimensional feature vectors is then used to identify a fixed number of clusters in the feature space. Based on this general description of the overall wave field by a mixture of multivariate Gaussians, we are able to learn particular event classifiers (here, hidden Markov models) from a single waveform example. To show the capabilities of this new approach we apply the algorithm to a data set recorded at Soufriere Hills volcano, Montserrat. Supported by very high classification rates, we conclude that the suggested approach provides a valuable tool for volcano monitoring systems.
Constructing a hidden Markov Model based earthquake detector: application to induced seismicity
(2012)
The triggering or detection of seismic events out of a continuous seismic data stream is one of the key issues of an automatic or semi-automatic seismic monitoring system. In the case of dense networks, either local or global, most of the implemented trigger algorithms are based on a large number of active stations. However, in the case of only few available stations or small events, for example, like in monitoring volcanoes or hydrothermal power plants, common triggers often show high false alarms. In such cases detection algorithms are of interest, which show reasonable performance when operating even on a single station. In this context, we apply Hidden Markov Models (HMM) which are algorithms borrowed from speech recognition. However, many pitfalls need to be avoided to apply speech recognition technology directly to earthquake detection. We show the fit of the model parameters in an innovative way. State clustering is introduced to refine the intrinsically assumed time dependency of the HMMs and we explain the effect coda has on the recognition results. The methodology is then used for the detection of anthropogenicly induced earthquakes for which we demonstrate for a period of 3.9 months of continuous data that the single station HMM earthquake detector can achieve similar detection rates as a common trigger in combination with coincidence sums over two stations. To show the general applicability of state clustering we apply the proposed method also to earthquake classification at Mt. Merapi volcano, Indonesia.
Forecasting seismo-volcanic activity by using the dynamical behavior of volcanic earthquake rates
(2012)
We present a novel approach for short-term forecasting of volcano seismic activity. Volcanic earthquakes can be seen as a response mechanism of the earth crust to stresses induced by magma injection. From this point of view the temporal evolution of seismicity can be represented as a diffusion process which compensates pressure differences. By means of this dynamical approach we are able to estimate the system behavior in the near future which in turn allows us to forecast the evolution of the earthquake rate for the next time span from actual and past observations. For this purpose we model the earthquake rate as a random walk process embedded in a moving and deforming potential function. The center of the potential function is given by a moving average of the random walk's trace. We successfully apply this procedure to estimate the next day seismicity at Soufriere Hills volcano, Montserrat, over a time period of six years. When comparing the dynamical approach to the well known method of material failure forecast we find much better predictions of the critical stages of volcanic activity using the new approach.
Other than commonly assumed in seismology, the phase velocity of Rayleigh waves is not necessarily a single-valued function of frequency. In fact, a single Rayleigh mode can exist with three different values of phase velocity at one frequency. We demonstrate this for the first higher mode on a realistic shallow seismic structure of a homogeneous layer of unconsolidated sediments on top of a half-space of solid rock (LOH). In the case of LOH a significant contrast to the half-space is required to produce the phenomenon. In a simpler structure of a homogeneous layer with fixed (rigid) bottom (LFB) the phenomenon exists for values of Poisson's ratio between 0.19 and 0.5 and is most pronounced for P-wave velocity being three times S-wave velocity (Poisson's ratio of 0.4375). A pavement-like structure (PAV) of two layers on top of a half-space produces the multivaluedness for the fundamental mode. Programs for the computation of synthetic dispersion curves are prone to trouble in such cases. Many of them use mode-follower algorithms which loose track of the dispersion curve and miss the multivalued section. We show results for well established programs. Their inability to properly handle these cases might be one reason why the phenomenon of multivaluedness went unnoticed in seismological Rayleigh wave research for so long. For the very same reason methods of dispersion analysis must fail if they imply wave number k(l)(omega) for the lth Rayleigh mode to be a single-valued function of frequency.. This applies in particular to deconvolution methods like phase-matched filters. We demonstrate that a slant-stack analysis fails in the multivalued section, while a Fourier-Bessel transformation captures the complete Rayleigh-wave signal. Waves of finite bandwidth in the multivalued section propagate with positive group-velocity and negative phase-velocity. Their eigenfunctions appear conventional and contain no conspicuous feature.
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/LTAtrigger 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 system.
Surface wave methods gained in the past decades a primary role in many seismic projects. Specifically, they are often used to retrieve a 1D shear wave velocity model or to estimate the V-s,V-30 at a site. The complexity of the interpretation process and the variety of possible approaches to surface wave analysis make it very hard to set a fixed standard to assure quality and reliability of the results. The present guidelines provide practical information on the acquisition and analysis of surface wave data by giving some basic principles and specific suggestions related to the most common situations. They are primarily targeted to non-expert users approaching surface wave testing, but can be useful to specialists in the field as a general reference. The guidelines are based on the experience gained within the InterPACIFIC project and on the expertise of the participants in acquisition and analysis of surface wave data.
In the present study, we investigated the dispersion characteristics of medium-to-long period Rayleigh waves (2 s < T < 20 s) using both single-station techniques (multiple-filter analysis, and phase-match filter) and multichannel techniques (horizontal slowness [p] and angular frequency [omega] stack, and cross-correlation) to determine the velocity structure for the Mt. Etna volcano. We applied these techniques to a dataset of teleseisms, as regional and local earthquakes recorded by two broad-band seismic arrays installed at Mt. Etna in 2002 and 2005, during two seismic surveys organized by the Istituto Nazionale di Geofisica e Vulcanologia (INGV), sezione di Napoli. The dispersion curves obtained showed phase velocities ranging from 1.5 km/s to 4.0 km/s in the frequency band 0.05 Hz to 0.45 Hz. We inverted the average phase velocity dispersion curves using a non-linear approach, to obtain a set of shear-wave velocity models with maximum resolution depths of 25 km to 30 km. Moreover, the presence of lateral velocity contrasts was checked by dividing the whole array into seven triangular sub-arrays and inverting the dispersion curves relative to each triangle.
Ambient vibration measurements with small, temporary arrays that produce estimates of surface wave dispersion have become increasingly popular as a low-cost, non-invasive tool for site characterisation. An important requirement for these measurements to be meaningful, however, is the temporal consistency and repeatability of the resulting dispersion and spatial autocorrelation curve estimates. Data acquired within several European research projects (NERIES task JRA4, SESAME, and other multinational experiments) offer the chance to investigate the variability of the derived data products. The dataset analysed here consists of repeated array measurements, with several years of time elapsed between them. The measurements were conducted by different groups in different seasons, using different instrumentations and array layouts, at six sites in Greece and Italy. Ambient vibration amplitude spectra and locations of dominant sources vary between the two measurements at each location. Still, analysis indicates that this does not influence the derived dispersion information, which is stable in time and neither influenced by the instrumentation nor the analyst. The frequency range over which the dispersion curves and spatial autocorrelation curves can be reliably estimated depends on the array dimensions (minimum and maximum aperture) used in the specific deployment, though, and may accordingly vary between the repeated experiments. The relative contribution of Rayleigh and Love waves to the wavefield can likewise change between repeated measurements. The observed relative contribution of Rayleigh waves is generally at or below 50%, with especially low values for the rural sites. Besides, the visibility of higher modes depends on the noise wavefield conditions. The similarity of the dispersion and autocorrelation curves measured at each site indicates that the curves are stable, mainly determined by the sub-surface structure, and can thus be used to derive velocity information with depth. Differences between velocity models for the same site derived from independently determined dispersion and autocorrelation curves-as observed in other studies-are consequently not adequately explained by uncertainties in the measurement part.
This study presents an unsupervised feature selection and learning approach for the discovery and intuitive imaging of significant temporal patterns in seismic single-station or network recordings. For this purpose, the data are parametrized by real-valued feature vectors for short time windows using standard analysis tools for seismic data, such as frequency-wavenumber, polarization, and spectral analysis. We use Self-Organizing Maps (SOMs) for a data-driven feature selection, visualization and clustering procedure, which is in particular suitable for high-dimensional data sets. Our feature selection method is based on significance testing using the Wald-Wolfowitz runs test for-individual features and on correlation hunting with SOMs in feature subsets. Using synthetics composed of Rayleigh and Love waves and real-world data, we show the robustness and the improved discriminative power of that approach compared to feature subsets manually selected from individual wavefield parametrization methods. Furthermore, the capability of the clustering and visualization techniques to investigate the discrimination of wave phases is shown by means of synthetic waveforms and regional earthquake recordings.
The inversion of surface-wave dispersion curve to derive shear-wave velocity profile is a very delicate process dealing with a nonunique problem, which is strongly dependent on the model space parameterization. When independent and reliable information is not available, the selection of most representative models within the ensemble produced. by the inversion is often difficult. We implemented a strategy in the inversion of dispersion curves able to investigate the influence of the parameterization of the model space and to select a "best" class of models. We analyzed surface-wave dispersion curves measured at 14 European strong..-motion sites within the NERIES EC-Project. We focused on the inversion task exploring the model space by means of four distinct pararneterization classes composed of layers progressively added over a half-space. The classes differ in the definition of the shear-wave velocity profile; we considered models with uniform velocity as well as models with increasing velocity with depth. At each site and for each model parameterization, we performed an extensive surface-wave inversion (200,100 models for five seeds) using the conditional neighborhood algorithm. We addressed the model evaluation following the corrected Akaike's information criterion (AlCc) that combines the concept of misfit to the number of degrees of freedom of the system. The misfit was computed as least-squares estimation between theoretical and observed dispersion curve. The model complexity was accounted in a penalty term by AlCc. By applying such inversion strategy on 14 strong-motion sites, we found that the best parameterization of the model space is mostly three to four layers over a half-space: where the shear-wave velocity of the uppermost layers can follow uniform or power-law dependence with depth. The shear-wave velocity profiles derived by inversion agree with shear-wave velocity profiles provided by borehole surveys at approximately 80% of the sites.
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.
We apply and evaluate a recent machine learning method for the automatic classification of seismic waveforms. The method relies on Dynamic Bayesian Networks (DBN) and supervised learning to improve the detection capabilities at 3C seismic stations. A time-frequency decomposition provides the basis for the required signal characteristics we need in order to derive the features defining typical "signal" and "noise" patterns. Each pattern class is modeled by a DBN, specifying the interrelationships of the derived features in the time-frequency plane. Subsequently, the models are trained using previously labeled segments of seismic data. The DBN models can now be compared against in order to determine the likelihood of new incoming seismic waveform segments to be either signal or noise. As the noise characteristics of seismic stations varies smoothly in time (seasonal variation as well as anthropogenic influence), we accommodate in our approach for a continuous adaptation of the DBN model that is associated with the noise class. Given the difficulty for obtaining a golden standard for real data (ground truth) the proof of concept and evaluation is shown by conducting experiments based on 3C seismic data from the International Monitoring Stations, BOSA and LPAZ.
The main scope of the InterPACIFIC (Intercomparison of methods for site parameter and velocity profile characterization) project is to assess the reliability of in-hole and surface-wave methods, used for estimating shear wave velocity. Three test-sites with different subsurface conditions were chosen: a soft soil, a stiff soil and a rock outcrop. This paper reports the surface-wave methods results. Specifically 14 teams of expert users analysed the same experimental surface-wave datasets, consisting of both passive and active data. Each team adopted their own strategy to retrieve the dispersion curve and the shear-wave velocity profile at each site. Despite different approaches, the dispersion curves are quite in agreement with each other. Conversely, the shear-wave velocity profiles show a certain variability that increases in correspondence of major stratigraphic interfaces. This larger variability is mainly due to non-uniqueness of the solution and lateral variability. As expected, the observed variability in V-s,V-30 estimatesis small, as solution non-uniqueness plays a limited role. (C) 2015 Elsevier Ltd. All rights reserved.
The resonance frequency of the transmission response in layered half-space model is important in the study of site effect because it is the frequency where the shake-ability of the ground is enhanced significantly. In practice, it is often determined by the H/V ratio technique in which the peak frequency of recorded H/V spectral ratio is interpreted as the resonance frequency. Despite of its importance, there has not been any formula of the resonance frequency of the layered half-space structure. In this paper, a simple approximate formula of the fundamental resonance frequency is presented after an exact formula in explicit form of the response function of vertically SH incident wave is obtained. The formula is in similar form with the one used in H/V ratio technique but it reflects several major effects of the model to the resonance frequency such as the arrangement of layers, the impedance contrast between layers and the half-space. Therefore, it could be considered as an improved formula used in H/V ratio technique. The formula also reflects the consistency between two approaches of the H/V ratio technique based on SH body waves or Rayleigh surface waves on the peak frequency under high impedance contrast condition. This formula is in explicit form and, therefore, may be used in the direct and inverse problem efficiently. A numerical illustration of the improved formula for an actual layered half-space model already investigated by H/V ratio technique is presented to demonstrate its new features and its improvement to the currently used formula.
In this paper we present densely sampled fumarole temperature data, recorded continuously at a high-temperature fumarole of Mt. Merapi volcano (Indonesia). These temperature time series are correlated with continuous records of rainfall and seismic waveform data collected at the Indonesian - German multi-parameter monitoring network. The correlation analysis of fumarole temperature and precipitation data shows a clear influence of tropical rain events on fumarole temperature. In addition, there is some evidence that rainfall may influence seismicity rates, indicating interaction of meteoric water with the volcanic system. Knowledge about such interactions is important, as lava dome instabilities caused by heavy-precipitation events may result in pyroclastic flows. Apart from the strong external influences on fumarole temperature and seismicity rate, which may conceal smaller signals caused by volcanic degassing processes, the analysis of fumarole temperature and seismic data indicates a statistically significant correlation between a certain type of seismic activity and an increase in fumarole temperature. This certain type of seismic activity consists of a seismic cluster of several high-frequency transients and an ultra-long-period signal (< 0.002 Hz), which are best observed using a broadband seismometer deployed at a distance of 600 m from the active lava dome. The corresponding change in fumarole temperature starts a few minutes after the ultra-long-period signal and simultaneously with the high-frequency seismic cluster. The change in fumarole temperature, an increase of 5 degreesC on average, resembles a smoothed step. Fifty-four occurrences of simultaneous high-frequency seismic cluster, ultra-long period signal and increase of fumarole temperature have been identified in the data set from August 2000 to January 2001. The observed signals appear to correspond to degassing processes in the summit region of Mt. Merapi. (C) 2004 Elsevier B.V. All rights reserved
Earthquake rupture length and width estimates are in demand in many seismological applications. Earthquake magnitude estimates are often available, whereas the geometrical extensions of the rupture fault mostly are lacking. Therefore, scaling relations are needed to derive length and width from magnitude. Most frequently used are the relationships of Wells and Coppersmith (1994) derived on the basis of a large dataset including all slip types with the exception of thrust faulting events in subduction environments. However, there are many applications dealing with earthquakes in subduction zones because of their high seismic and tsunamigenic potential. There are no well-established scaling relations for moment magnitude and length/width for subduction events. Within this study, we compiled a large database of source parameter estimates of 283 earthquakes. All focal mechanisms are represented, but special focus is set on (large) subduction zone events, in particular. Scaling relations were fitted with linear least-square as well as orthogonal regression and analyzed regarding the difference between continental and subduction zone/oceanic relationships. Additionally, the effect of technical progress in earthquake parameter estimation on scaling relations was tested as well as the influence of different fault mechanisms. For a given moment magnitude we found shorter but wider rupture areas of thrust events compared to Wells and Coppersmith (1994). The thrust event relationships for pure continental and pure subduction zone rupture areas were found to be almost identical. The scaling relations differ significantly for slip types. The exclusion of events prior to 1964 when the worldwide standard seismic network was established resulted in a remarkable effect on strike-slip scaling relations: the data do not show any saturation of rupture width of strike- slip earthquakes. Generally, rupture area seems to scale with mean slip independent of magnitude. The aspect ratio L/W, however, depends on moment and differs for each slip type.
Tsunami early warning (TEW) is a challenging task as a decision has to be made within few minutes on the basis of incomplete and error-prone data. Deterministic warning systems have difficulties in integrating and quantifying the intrinsic uncertainties. In contrast, probabilistic approaches provide a framework that handles uncertainties in a natural way. Recently, we have proposed a method using Bayesian networks (BNs) that takes into account the uncertainties of seismic source parameter estimates in TEW. In this follow-up study, the method is applied to 10 recent large earthquakes offshore Sumatra and tested for its performance. We have evaluated both the general model performance given the best knowledge we have today about the source parameters of the 10 events and the corresponding response on seismic source information evaluated in real-time. We find that the resulting site-specific warning level probabilities represent well the available tsunami wave measurements and observations. Difficulties occur in the real-time tsunami assessment if the moment magnitude estimate is severely over- or underestimated. In general, the probabilistic analysis reveals a considerably large range of uncertainties in the near-field TEW. By quantifying the uncertainties the BN analysis provides important additional information to a decision maker in a warning centre to deal with the complexity in TEW and to reason under uncertainty.