@article{Scherbaum1997, author = {Scherbaum, Frank}, title = {Zero Phase FIR filters in digital seismic acquisition systems : blessing or curse}, year = {1997}, language = {en} } @article{ZaliOhrnbergerScherbaumetal.2021, author = {Zali, Zahra and Ohrnberger, Matthias and Scherbaum, Frank and Cotton, Fabrice and Eibl, Eva P. S.}, title = {Volcanic tremor extraction and earthquake detection using music information retrieval algorithms}, series = {Seismological research letters}, volume = {92}, journal = {Seismological research letters}, number = {6}, publisher = {Seismological Society of America}, address = {Boulder, Colo.}, issn = {0895-0695}, doi = {10.1785/0220210016}, pages = {3668 -- 3681}, year = {2021}, abstract = {Volcanic tremor signals are usually observed before or during volcanic eruptions and must be monitored to evaluate the volcanic activity. A challenge in studying seismic signals of volcanic origin is the coexistence of transient signal swarms and long-lasting volcanic tremor signals. Separating transient events from volcanic tremors can, therefore, contrib-ute to improving upon our understanding of the underlying physical processes. Exploiting the idea of harmonic-percussive separation in musical signal processing, we develop a method to extract the harmonic volcanic tremor signals and to detect tran-sient events from seismic recordings. Based on the similarity properties of spectrogram frames in the time-frequency domain, we decompose the signal into two separate spec-trograms representing repeating (harmonic) and nonrepeating (transient) patterns, which correspond to volcanic tremor signals and earthquake signals, respectively. We reconstruct the harmonic tremor signal in the time domain from the complex spectrogram of the repeating pattern by only considering the phase components for the frequency range in which the tremor amplitude spectrum is significantly contribut-ing to the energy of the signal. The reconstructed signal is, therefore, clean tremor signal without transient events. Furthermore, we derive a characteristic function suitable for the detection of tran-sient events (e.g., earthquakes) by integrating amplitudes of the nonrepeating spectro-gram over frequency at each time frame. Considering transient events like earthquakes, 78\% of the events are detected for signal-to-noise ratio = 0.1 in our semisynthetic tests. In addition, we compared the number of detected earthquakes using our method for one month of continuous data recorded during the Holuhraun 2014-2015 eruption in Iceland with the bulletin presented in Agustsdottir et al. (2019). Our single station event detection algorithm identified 84\% of the bulletin events. Moreover, we detected a total of 12,619 events, which is more than twice the number of the bulletin events.}, language = {en} } @article{KoehlerOhrnbergerScherbaum2009, author = {Koehler, Andreas and Ohrnberger, Matthias and Scherbaum, Frank}, title = {Unsupervised feature selection and general pattern discovery using Self-Organizing Maps for gaining insights into the nature of seismic wavefields}, issn = {0098-3004}, doi = {10.1016/j.cageo.2009.02.004}, year = {2009}, abstract = {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.}, language = {en} } @article{DelavaudCottonAkkaretal.2012, author = {Delavaud, Elise and Cotton, Fabrice and Akkar, Sinan and Scherbaum, Frank and Danciu, Laurentiu and Beauval, Celine and Drouet, Stephane and Douglas, John and Basili, Roberto and Sandikkaya, M. Abdullah and Segou, Margaret and Faccioli, Ezio and Theodoulidis, Nikos}, title = {Toward a ground-motion logic tree for probabilistic seismic hazard assessment in Europe}, series = {Journal of seismology}, volume = {16}, journal = {Journal of seismology}, number = {3}, publisher = {Springer}, address = {Dordrecht}, issn = {1383-4649}, doi = {10.1007/s10950-012-9281-z}, pages = {451 -- 473}, year = {2012}, abstract = {The Seismic Hazard Harmonization in Europe (SHARE) project, which began in June 2009, aims at establishing new standards for probabilistic seismic hazard assessment in the Euro-Mediterranean region. In this context, a logic tree for ground-motion prediction in Europe has been constructed. Ground-motion prediction equations (GMPEs) and weights have been determined so that the logic tree captures epistemic uncertainty in ground-motion prediction for six different tectonic regimes in Europe. Here we present the strategy that we adopted to build such a logic tree. This strategy has the particularity of combining two complementary and independent approaches: expert judgment and data testing. A set of six experts was asked to weight pre-selected GMPEs while the ability of these GMPEs to predict available data was evaluated with the method of Scherbaum et al. (Bull Seismol Soc Am 99:3234-3247, 2009). Results of both approaches were taken into account to commonly select the smallest set of GMPEs to capture the uncertainty in ground-motion prediction in Europe. For stable continental regions, two models, both from eastern North America, have been selected for shields, and three GMPEs from active shallow crustal regions have been added for continental crust. For subduction zones, four models, all non-European, have been chosen. Finally, for active shallow crustal regions, we selected four models, each of them from a different host region but only two of them were kept for long periods. In most cases, a common agreement has been also reached for the weights. In case of divergence, a sensitivity analysis of the weights on the seismic hazard has been conducted, showing that once the GMPEs have been selected, the associated set of weights has a smaller influence on the hazard.}, language = {en} } @book{ScherbaumMzhavanadzeArometal.2020, author = {Scherbaum, Frank and Mzhavanadze, Nana and Arom, Simha and Rosenzweig, Sebastian and M{\"u}ller, Meinard}, title = {Tonal Organization of the Erkomaishvili Dataset: Pitches, Scales, Melodies and Harmonies}, series = {Computational Analysis Of Traditional Georgian Vocal Music}, journal = {Computational Analysis Of Traditional Georgian Vocal Music}, number = {1}, editor = {Scherbaum, Frank}, publisher = {Universit{\"a}tsverlag Potsdam}, address = {Potsdam}, issn = {2702-2641}, doi = {10.25932/publishup-47614}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-476141}, publisher = {Universit{\"a}t Potsdam}, pages = {64}, year = {2020}, abstract = {In this study we examine the tonal organization of a series of recordings of liturgical chants, sung in 1966 by the Georgian master singer Artem Erkomaishvili. This dataset is the oldest corpus of Georgian chants from which the time synchronous F0-trajectories for all three voices have been reliably determined (M{\"u}ller et al. 2017). It is therefore of outstanding importance for the understanding of the tuning principles of traditional Georgian vocal music. The aim of the present study is to use various computational methods to analyze what these recordings can contribute to the ongoing scientific dispute about traditional Georgian tuning systems. Starting point for the present analysis is the re-release of the original audio data together with estimated fundamental frequency (F0) trajectories for each of the three voices, beat annotations, and digital scores (Rosenzweig et al. 2020). We present synoptic models for the pitch and the harmonic interval distributions, which are the first of such models for which the complete Erkomaishvili dataset was used. We show that these distributions can be very compactly be expressed as Gaussian mixture models, anchored on discrete sets of pitch or interval values for the pitch and interval distributions, respectively. As part of our study we demonstrate that these pitch values, which we refer to as scale pitches, and which are determined as the mean values of the Gaussian mixture elements, define the scale degrees of the melodic sound scales which build the skeleton of Artem Erkomaishvili's intonation. The observation of consistent pitch bending of notes in melodic phrases, which appear in identical form in a group of chants, as well as the observation of harmonically driven intonation adjustments, which are clearly documented for all pure harmonic intervals, demonstrate that Artem Erkomaishvili intentionally deviates from the scale pitch skeleton quite freely. As a central result of our study, we proof that this melodic freedom is always constrained by the attracting influence of the scale pitches. Deviations of the F0-values of individual note events from the scale pitches at one instance of time are compensated for in the subsequent melodic steps. This suggests a deviation-compensation mechanism at the core of Artem Erkomaishvili's melody generation, which clearly honors the scales but still allows for a large degree of melodic flexibility. This model, which summarizes all partial aspects of our analysis, is consistent with the melodic scale models derived from the observed pitch distributions, as well as with the melodic and harmonic interval distributions. In addition to the tangible results of our work, we believe that our work has general implications for the determination of tuning models from audio data, in particular for non-tempered music.}, language = {en} } @article{AlAtikAbrahamsonBommeretal.2010, author = {Al Atik, Linda and Abrahamson, Norman A. and Bommer, Julian J. and Scherbaum, Frank and Cotton, Fabrice and Kuehn, Nicolas}, title = {The variability of ground-motion prediction models and its components}, issn = {0895-0695}, doi = {10.1785/gssrl.81.5.794}, year = {2010}, language = {en} } @article{BeauvalHainzlScherbaum2006, author = {Beauval, Celine and Hainzl, Sebastian and Scherbaum, Frank}, title = {The impact of the spatial uniform distribution of seismicity on probabilistic seismic-hazard estimation}, series = {Bulletin of the Seismological Society of America}, volume = {96}, journal = {Bulletin of the Seismological Society of America}, number = {6}, publisher = {GeoScienceWorld}, address = {Alexandria, Va.}, issn = {0037-1106}, doi = {10.1785/0120060073}, pages = {2465 -- 2471}, year = {2006}, abstract = {The first step in the estimation of probabilistic seismic hazard in a region commonly consists of the definition and characterization of the relevant seismic sources. Because in low-seismicity regions seismicity is often rather diffuse and faults are difficult to identify, large areal source zones are mostly used. The corresponding hypothesis is that seismicity is uniformly distributed inside each areal seismic source zone. In this study, the impact of this hypothesis on the probabilistic hazard estimation is quantified through the generation of synthetic spatial seismicity distributions. Fractal seismicity distributions are generated inside a given source zone and probabilistic hazard is computed for a set of sites located inside this zone. In our study, the impact of the spatial seismicity distribution is defined as the deviation from the hazard value obtained for a spatially uniform seismicity distribution. From the generation of a large number of synthetic distributions, the correlation between the fractal dimension D and the impact is derived. The results show that the assumption of spatially uniform seismicity tends to bias the hazard to higher values. The correlation can be used to determine the systematic biases and uncertainties for hazard estimations in real cases, where the fractal dimension has been determined. We apply the technique in Germany (Cologne area) and in France (Alps).}, language = {en} } @article{RietbrockScherbaum1998, author = {Rietbrock, Andreas and Scherbaum, Frank}, title = {The GIANT analysis system (Graphical Interaktive Aftershock Network Toolbox)}, year = {1998}, language = {en} } @article{ScherbaumCottonStaedtke2006, author = {Scherbaum, Frank and Cotton, Fabrice and Staedtke, Helmut}, title = {The estimation of minimum-misfit stochastic models from empirical ground-motion prediction equations}, doi = {10.1785/0120050015}, year = {2006}, abstract = {In areas of moderate to low seismic activity there is commonly a lack of recorded strong ground motion. As a consequence, the prediction of ground motion expected for hypothetical future earthquakes is often performed by employing empirical models from other regions. In this context, Campbell's hybrid empirical approach (Campbell, 2003, 2004) provides a methodological framework to adapt ground-motion prediction equations to arbitrary target regions by using response spectral host-to-target-region-conversion filters. For this purpose, the empirical ground-motion prediction equation has to be quantified in terms of a stochastic model. The problem we address here is how to do this in a systematic way and how to assess the corresponding uncertainties. For the determination of the model parameters we use a genetic algorithm search. The stochastic model spectra were calculated by using a speed-optimized version of SMSIM (Boore, 2000). For most of the empirical ground-motion models, we obtain sets of stochastic models that match the empirical models within the full magnitude and distance ranges of their generating data sets fairly well. The overall quality of fit and the resulting model parameter sets strongly depend on the particular choice of the distance metric used for the stochastic model. We suggest the use of the hypocentral distance metric for the stochastic Simulation of strong ground motion because it provides the lowest-misfit stochastic models for most empirical equations. This is in agreement with the results of two recent studies of hypocenter locations in finite-source models which indicate that hypocenters are often located close to regions of large slip (Mai et al., 2005; Manighetti et al., 2005). Because essentially all empirical ground-motion prediction equations contain data from different geographical regions, the model parameters corresponding to the lowest-misfit stochastic models cannot necessarily be expected to represent single, physically realizable host regions but to model the generating data sets in an average way. In addition, the differences between the lowest-misfit stochastic models and the empirical ground-motion prediction equation are strongly distance, magnitude, and frequency dependent, which, according to the laws of uncertainty propagation, will increase the variance of the corresponding hybrid empirical model predictions (Scherbaum et al., 2005). As a consequence, the selection of empirical ground-motion models for host-to-target-region conversions requires considerable judgment of the ground-motion analyst}, language = {en} } @article{ScherbaumWeberBorm2000, author = {Scherbaum, Frank and Weber, Michael H. and Borm, G.}, title = {The deep seismological lab in the KTB borehole: Status 1999}, year = {2000}, language = {en} } @article{WeberAbuAyyashAbueladasetal.2004, author = {Weber, Michael H. and Abu-Ayyash, Khalil and Abueladas, Abdel-Rahman and Agnon, Amotz and Al-Amoush, H. and Babeyko, Andrey and Bartov, Yosef and Baumann, M. and Ben-Avraham, Zvi and Bock, G{\"u}nter and Bribach, Jens and El-Kelani, R. and Forster, A. and F{\"o}rster, Hans-J{\"u}rgen and Frieslander, U. and Garfunkel, Zvi and Grunewald, Steffen and Gotze, Hans-J{\"u}rgen and Haak, Volker and Haberland, Christian and Hassouneh, Mohammed and Helwig, S. and Hofstetter, Alfons and Jackel, K. H. and Kesten, Dagmar and Kind, Rainer and Maercklin, Nils and Mechie, James and Mohsen, Amjad and Neubauer, F. M. and Oberh{\"a}nsli, Roland and Qabbani, I. and Ritter, O. and Rumpker, G. and Rybakov, M. and Ryberg, Trond and Scherbaum, Frank and Schmidt, J. and Schulze, A. and Sobolev, Stephan Vladimir and Stiller, M. and Th,}, title = {The crustal structure of the Dead Sea Transform}, year = {2004}, abstract = {To address one of the central questions of plate tectonics-How do large transform systems work and what are their typical features?-seismic investigations across the Dead Sea Transform (DST), the boundary between the African and Arabian plates in the Middle East, were conducted for the first time. A major component of these investigations was a combined reflection/ refraction survey across the territories of Palestine, Israel and Jordan. The main results of this study are: (1) The seismic basement is offset by 3-5 km under the DST, (2) The DST cuts through the entire crust, broadening in the lower crust, (3) Strong lower crustal reflectors are imaged only on one side of the DST, (4) The seismic velocity sections show a steady increase in the depth of the crust-mantle transition (Moho) from 26 km at the Mediterranean to 39 km under the Jordan highlands, with only a small but visible, asymmetric topography of the Moho under the DST. These observations can be linked to the left-lateral movement of 105 km of the two plates in the last 17 Myr, accompanied by strong deformation within a narrow zone cutting through the entire crust. Comparing the DST and the San Andreas Fault (SAF) system, a strong asymmetry in subhorizontal lower crustal reflectors and a deep reaching deformation zone both occur around the DST and the SAF. The fact that such lower crustal reflectors and deep deformation zones are observed in such different transform systems suggests that these structures are possibly fundamental features of large transform plate boundaries}, language = {en} } @article{BommerAbrahamsonStrasseretal.2004, author = {Bommer, Julian J. and Abrahamson, Norman A. and Strasser, F. O. and Pecker, Alain and Bard, Pierre-Yves and Bungum, Hilmar and Cotton, Fabrice and F{\"a}h, Donat and Sabetta, F. and Scherbaum, Frank and Studer, Jost}, title = {The challenge of defining upper bounds on earthquake ground motions}, issn = {0895-0695}, year = {2004}, language = {en} } @article{KruegerScherbaum2014, author = {Kr{\"u}ger, Frank and Scherbaum, Frank}, title = {The 29 September 1969, Ceres, South Africa, Earthquake: full waveform moment tensor inversion for point source and kinematic source parameters}, series = {Bulletin of the Seismological Society of America}, volume = {104}, journal = {Bulletin of the Seismological Society of America}, number = {1}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0037-1106}, doi = {10.1785/0120130209}, pages = {576 -- 581}, year = {2014}, abstract = {The Ceres earthquake of 29 September 1969 is the largest known earthquake in southern Africa. Digitized analog recordings from Worldwide Standardized Seismographic Network stations (Powell and Fries, 1964) are used to retrieve the point source moment tensor and the most likely centroid depth of the event using full waveform modeling. A scalar seismic moment of 2.2-2.4 x 10(18) N center dot m corresponding to a moment magnitude of 6.2-6.3 is found. The analysis confirms the pure strike-slip mechanism previously determined from onset polarities by Green and Bloch (1971). Overall good agreement with the fault orientation previously estimated from local aftershock recordings is found. The centroid depth can be constrained to be less than 15 km. In a second analysis step, we use a higher order moment tensor based inversion scheme for simple extended rupture models to constrain the lateral fault dimensions. We find rupture propagated unilaterally for 4.7 s from east-southwest to west-northwest for about 17 km ( average rupture velocity of about 3: 1 km/s).}, language = {en} } @article{DelavaudScherbaumKuehnetal.2012, author = {Delavaud, Elise and Scherbaum, Frank and K{\"u}hn, Nicolas and Allen, Trevor}, title = {Testing the global applicability of ground-motion prediction equations for active shallow crustal regions}, series = {Bulletin of the Seismological Society of America}, volume = {102}, journal = {Bulletin of the Seismological Society of America}, number = {2}, publisher = {Seismological Society of America}, address = {El Cerrito}, issn = {0037-1106}, doi = {10.1785/0120110113}, pages = {707 -- 721}, year = {2012}, abstract = {Large research initiatives such as the Global Earthquake Model (GEM) or the Seismic HAzard haRmonization in Europe (SHARE) projects concentrate a great collaborative effort on defining a global standard for seismic hazard estimations. In this context, there is an increasing need for identifying ground-motion prediction equations (GMPEs) that can be applied at both global and regional scale. With increasing amounts of strong-motion records that are now available worldwide, observational data can provide a valuable resource to tackle this question. Using the global dataset of Allen and Wald (2009), we evaluate the ability of 11 GMPEs to predict ground-motion in different active shallow crustal regions worldwide. Adopting the approach of Scherbaum et al. (2009), we rank these GMPEs according to their likelihood of having generated the data. In particular, we estimate how strongly data support or reject the models with respect to the state of noninformativeness defined by a uniform weighting. Such rankings derived from this particular global dataset enable us to explore the potential of GMPEs to predict ground motions in their host region and also in other regions depending on the magnitude and distance considered. In the ranking process, we particularly focus on the influence of the distribution of the testing dataset compared with the GMPE's native dataset. One of the results of this study is that some nonindigenous models present a high degree of consistency with the data from a target region. Two models in particular demonstrated a strong power of geographically wide applicability in different geographic regions with respect to the testing dataset: the models of Akkar and Bommer (2010) and Chiou et al. (2010).}, language = {en} } @article{BoraCottonScherbaumetal.2017, author = {Bora, Sanjay Singh and Cotton, Fabrice and Scherbaum, Frank and Edwards, Benjamin and Traversa, Paola}, title = {Stochastic source, path and site attenuation parameters and associated variabilities for shallow crustal European earthquakes}, series = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering}, volume = {15}, journal = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering}, publisher = {Springer}, address = {Dordrecht}, issn = {1570-761X}, doi = {10.1007/s10518-017-0167-x}, pages = {4531 -- 4561}, year = {2017}, abstract = {We have analyzed the recently developed pan-European strong motion database, RESORCE-2012: spectral parameters, such as stress drop (stress parameter, Delta sigma), anelastic attenuation (Q), near surface attenuation (kappa(0)) and site amplification have been estimated from observed strong motion recordings. The selected dataset exhibits a bilinear distance-dependent Q model with average kappa(0) value 0.0308 s. Strong regional variations in inelastic attenuation were also observed: frequency-independent Q(0) of 1462 and 601 were estimated for Turkish and Italian data respectively. Due to the strong coupling between Q and kappa(0), the regional variations in Q have strong impact on the estimation of near surface attenuation kappa(0). kappa(0) was estimated as 0.0457 and 0.0261 s for Turkey and Italy respectively. Furthermore, a detailed analysis of the variability in estimated kappa(0) revealed significant within-station variability. The linear site amplification factors were constrained from residual analysis at each station and site-class type. Using the regional Q(0) model and a site-class specific kappa(0), seismic moments (M-0) and source corner frequencies f (c) were estimated from the site corrected empirical Fourier spectra. Delta sigma did not exhibit magnitude dependence. The median Delta sigma value was obtained as 5.75 and 5.65 MPa from inverted and database magnitudes respectively. A comparison of response spectra from the stochastic model (derived herein) with that from (regional) ground motion prediction equations (GMPEs) suggests that the presented seismological parameters can be used to represent the corresponding seismological attributes of the regional GMPEs in a host-to-target adjustment framework. The analysis presented herein can be considered as an update of that undertaken for the previous Euro-Mediterranean strong motion database presented by Edwards and Fah (Geophys J Int 194(2):1190-1202, 2013a).}, language = {en} } @article{FaenzaHainzlScherbaum2009, author = {Faenza, Licia and Hainzl, Sebastian and Scherbaum, Frank}, title = {Statistical analysis of the Central-Europe seismicity}, issn = {0040-1951}, doi = {10.1016/j.tecto.2008.04.030}, year = {2009}, abstract = {The aim of this paper is to characterize the spatio-temporal distribution of Central-Europe seismicity. Specifically, by using a non-parametric statistical approach, the proportional hazard model, leading to an empirical estimation of the hazard function, we provide some constrains on the time behavior of earthquake generation mechanisms. The results indicate that the most conspicuous characteristics of M-w 4.0+ earthquakes is a temporal clustering lasting a couple of years. This suggests that the probability of occurrence increases immediately after a previous event. After a few years, the process becomes almost time independent. Furthermore, we investigate the cluster properties of the seismicity of Central-Europe, by comparing the obtained result with the one of synthetic catalogs generated by the epidemic type aftershock sequences (ETAS) model, which previously have been successfully applied for short term clustering. Our results indicate that the ETAS is not well suited to describe the seismicity as a whole, while it is able to capture the features of the short- term behaviour. Remarkably, similar results have been previously found for Italy using a higher magnitude threshold.}, language = {en} } @inproceedings{HainzlScherbaumZoeller2006, author = {Hainzl, Sebastian and Scherbaum, Frank and Z{\"o}ller, Gert}, title = {Spatiotemporal earthquake patterns}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-7267}, year = {2006}, abstract = {Interdisziplin{\"a}res Zentrum f{\"u}r Musterdynamik und Angewandte Fernerkundung Workshop vom 9. - 10. Februar 2006}, language = {en} } @article{FalsaperlaWassermannScherbaum2002, author = {Falsaperla, Susanna and Wassermann, Joachim and Scherbaum, Frank}, title = {Solid earth - 29. Polarization analyses of broadband seismic data recorded on Stromboli Volcano (Italy) from 1996 to 1999 (DOI 10.1029-2001GLO14300)}, year = {2002}, language = {en} } @misc{RoesslerHiemerBachetal.2009, author = {R{\"o}ßler, Dirk and Hiemer, Stephan and Bach, Christoph and Delavaud, Elise and Kr{\"u}ger, Frank and Ohrnberger, Matthias and Sauer, David and Scherbaum, Frank and Vollmer, Daniel}, title = {Small-aperture seismic array monitors Vogtland earthquake swarm in 2008/09}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-29185}, year = {2009}, abstract = {The most recent intense earthquake swarm in the Vogtland lasted from 6 October 2008 until January 2009. Greatest magnitudes exceeded M3.5 several times in October making it the greatest swarm since 1985/86. In contrast to the swarms in 1985 and 2000, seismic moment release was concentrated near swarm onset. Focal area and temporal evolution are similar to the swarm in 2000. Work hypothysis: uprising upper-mantle fluids trigger swarm earthquakes at low stress level. To monitor the seismicity, the University of Potsdam operated a small aperture seismic array at 10 km epicentral distance between 18 October 2008 and 18 March 2009. Consisting of 12 seismic stations and 3 additional microphones, the array is capable of detecting earthquakes from larger to very low magnitudes (M<-1) as well as associated air waves. We use array techniques to determine properties of the incoming wavefield: noise, direct P and S waves, and converted phases.}, language = {en} } @article{ThomasWeberWicksetal.1999, author = {Thomas, Ch. and Weber, Michael H. and Wicks, Chuck and Scherbaum, Frank}, title = {Small scatterers in the lower mantle observed at German broadband arrays}, year = {1999}, language = {en} } @article{HinzenReamerScherbaum2013, author = {Hinzen, Klaus-G and Reamer, Sharon K. and Scherbaum, Frank}, title = {Slow fourier transform}, series = {Seismological research letters}, volume = {84}, journal = {Seismological research letters}, number = {2}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0895-0695}, doi = {10.1785/0220120139}, pages = {251 -- 257}, year = {2013}, language = {en} } @article{MolkenthinScherbaumGriewanketal.2015, author = {Molkenthin, Christian and Scherbaum, Frank and Griewank, Andreas and K{\"u}hn, Nicolas and Stafford, Peter J. and Leovey, Hernan}, title = {Sensitivity of Probabilistic Seismic Hazard Obtained by Algorithmic Differentiation: A Feasibility Study}, series = {Bulletin of the Seismological Society of America}, volume = {105}, journal = {Bulletin of the Seismological Society of America}, number = {3}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0037-1106}, doi = {10.1785/0120140294}, pages = {1810 -- 1822}, year = {2015}, abstract = {Probabilistic seismic-hazard analysis (PSHA) is the current tool of the trade used to estimate the future seismic demands at a site of interest. A modern PSHA represents a complex framework that combines different models with numerous inputs. It is important to understand and assess the impact of these inputs on the model output in a quantitative way. Sensitivity analysis is a valuable tool for quantifying changes of a model output as inputs are perturbed, identifying critical input parameters, and obtaining insight about the model behavior. Differential sensitivity analysis relies on calculating first-order partial derivatives of the model output with respect to its inputs; however, obtaining the derivatives of complex models can be challenging. In this study, we show how differential sensitivity analysis of a complex framework such as PSHA can be carried out using algorithmic/automatic differentiation (AD). AD has already been successfully applied for sensitivity analyses in various domains such as oceanography and aerodynamics. First, we demonstrate the feasibility of the AD methodology by comparing AD-derived sensitivities with analytically derived sensitivities for a basic case of PSHA using a simple ground-motion prediction equation. Second, we derive sensitivities via AD for a more complex PSHA study using a stochastic simulation approach for the prediction of ground motions. The presented approach is general enough to accommodate more advanced PSHA studies of greater complexity.}, language = {en} } @article{BlaserKruegerOhrnbergeretal.2010, author = {Blaser, Lilian and Kr{\"u}ger, Frank and Ohrnberger, Matthias and Scherbaum, Frank}, title = {Scaling relations of earthquake source parameter estimates with special focus on subduction environment}, issn = {0037-1106}, doi = {10.1785/0120100111}, year = {2010}, abstract = {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.}, language = {en} } @article{BlaserOhrnbergerKruegeretal.2012, author = {Blaser, Lilian and Ohrnberger, Matthias and Kr{\"u}ger, Frank and Scherbaum, Frank}, title = {Probabilistic tsunami threat assessment of 10 recent earthquakes offshore Sumatra}, series = {Geophysical journal international}, volume = {188}, journal = {Geophysical journal international}, number = {3}, publisher = {Wiley-Blackwell}, address = {Malden}, issn = {0956-540X}, doi = {10.1111/j.1365-246X.2011.05324.x}, pages = {1273 -- 1284}, year = {2012}, abstract = {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.}, language = {en} } @article{BeauvalHainzlScherbaum2006, author = {Beauval, C{\´e}line and Hainzl, Sebastian and Scherbaum, Frank}, title = {Probabilistic seismic hazard estimation in low-seismicity regions considering non-Poissonian seismic occurrence}, issn = {0956-540X}, doi = {10.1111/j.1365-246X.2006.02863.x}, year = {2006}, abstract = {In low-seismicity regions, such as France or Germany, the estimation of probabilistic seismic hazard must cope with the difficult identification of active faults and with the low amount of seismic data available. Since the probabilistic hazard method was initiated, most studies assume a Poissonian occurrence of earthquakes. Here we propose a method that enables the inclusion of time and space dependences between earthquakes into the probabilistic estimation of hazard. Combining the seismicity model Epidemic Type Aftershocks-Sequence (ETAS) with a Monte Carlo technique, aftershocks are naturally accounted for in the hazard determination. The method is applied to the Pyrenees region in Southern France. The impact on hazard of declustering and of the usual assumption that earthquakes occur according to a Poisson process is quantified, showing that aftershocks contribute on average less than 5 per cent to the probabilistic hazard, with an upper bound around 18 per cent}, language = {en} } @article{ScherbaumCottonSmit2004, author = {Scherbaum, Frank and Cotton, Fabrice and Smit, P.}, title = {On the use of response spectral-reference data for the selection and ranking of ground-motion models for seismic-hazard analysis in regions of moderate seismicity : the case of rock motion}, issn = {0037-1106}, year = {2004}, abstract = {The use of ground-motion-prediction equations to estimate ground shaking has become a very popular approach for seismic-hazard assessment, especially in the framework of a logic-tree approach. Owing to the large number of existing published ground-motion models, however, the selection and ranking of appropriate models for a particular target area often pose serious practical problems. Here we show how observed around-motion records can help to guide this process in a systematic and comprehensible way. A key element in this context is a new, likelihood based, goodness-of-fit measure that has the property not only to quantify the model fit but also to measure in some degree how well the underlying statistical model assumptions are met. By design, this measure naturally scales between 0 and 1, with a value of 0.5 for a situation in which the model perfectly matches the sample distribution both in terms of mean and standard deviation. We have used it in combination with other goodness-of-fit measures to derive a simple classification scheme to quantify how well a candidate ground-rnotion-prediction equation models a particular set of observed-response spectra. This scheme is demonstrated to perform well in recognizing a number of popular ground-motion models from their rock-site- recording, subsets. This indicates its potential for aiding the assignment of logic-tree weights in a consistent and reproducible way. We have applied our scheme to the border region of France, Germany, and Switzerland where the M-w 4.8 St. Die earthquake of 22 February 2003 in eastern France recently provided a small set of observed-response spectra. These records are best modeled by the ground-motion-prediction equation of Berge-Thierry et al. (2003), which is based on the analysis of predominantly European data. The fact that the Swiss model of Bay et al. (2003) is not able to model the observed records in an acceptable way may indicate general problems arising from the use of weak-motion data for strong-motion prediction}, language = {en} } @article{BommerScherbaumBungumetal.2005, author = {Bommer, Julian J. and Scherbaum, Frank and Bungum, Hilmar and Cotton, Fabrice and Sabetta, F. and Abrahamson, Norman A.}, title = {On the use of logic trees for ground-motion prediction equations in seismic-hazard analysis}, issn = {0037-1106}, year = {2005}, abstract = {Logic trees are widely used in probabilistic seismic hazard analysis as a tool to capture the epistemic uncertainty associated with the seismogenic sources and the ground-motion prediction models used in estimating the hazard. Combining two or more ground-motion relations within a logic tree will generally require several conversions to be made, because there are several definitions available for both the predicted ground-motion parameters and the explanatory parameters within the predictive ground-motion relations. Procedures for making conversions for each of these factors are presented, using a suite of predictive equations in current use for illustration. The sensitivity of the resulting ground-motion models to these conversions is shown to be pronounced for some of the parameters, especially the measure of source-to-site distance, highlighting the need to take into account any incompatibilities among the selected equations. Procedures are also presented for assigning weights to the branches in the ground-motion section of the logic tree in a transparent fashion, considering both intrinsic merits of the individual equations and their degree of applicability to the particular application}, language = {en} } @article{BommerDouglasScherbaumetal.2010, author = {Bommer, Julian J. and Douglas, John and Scherbaum, Frank and Cotton, Fabrice and Bungum, Hilmar and Faeh, Donat}, title = {On the selection of ground-motion prediction equations for seismic hazard analysis}, issn = {0895-0695}, doi = {10.1785/gssrl.81.5.783}, year = {2010}, language = {en} } @article{HinzenWeberScherbaum2004, author = {Hinzen, K. G. and Weber, B. and Scherbaum, Frank}, title = {On the resolution of H/V measurements to determine sediment thickness, a case study across a normal fault in the Lower Rhine Embayment, Germany}, issn = {1363-2469}, year = {2004}, abstract = {In recent years, H/V measurements have been increasingly used to map the thickness of sediment fill in sedimentary basins in the context of seismic hazard assessment. This parameter is believed to be an important proxy for the site effects in sedimentary basins (e.g. in the Los Angeles basin). Here we present the results of a test using this approach across an active normal fault in a structurally well known situation. Measurements on a 50 km long profile with 1 km station spacing clearly show a change in the frequency of the fundamental peak of H/V ratios with increasing thickness of the sediment layer in the eastern part of the Lower Rhine Embayment. Subsequently, a section of 10 km length across the Erft-Sprung system, a normal fault with ca. 750 m vertical offset, was measured with a station distance of 100 m. Frequencies of the first and second peaks and the first trough in the H/V spectra are used in a simple resonance model to estimate depths of the bedrock. While the frequency of the first peak shows a large scatter for sediment depths larger than ca. 500 m, the frequency of the first trough follows the changing thickness of the sediments across the fault. The lateral resolution is in the range of the station distance of 100 m. A power law for the depth dependence of the S-wave velocity derived from down hole measurements in an earlier study [Budny, 1984] and power laws inverted from dispersion analysis of micro array measurements [Scherbaum et al., 2002] agree with the results from the H/V ratios of this study}, language = {en} } @article{TranThanhTuanScherbaumMalischewsky2011, author = {Tran Thanh Tuan, and Scherbaum, Frank and Malischewsky, Peter G.}, title = {On the relationship of peaks and troughs of the ellipticity (H/V) of Rayleigh waves and the transmission response of single layer over half-space models}, series = {Geophysical journal international}, volume = {184}, journal = {Geophysical journal international}, number = {2}, publisher = {Wiley-Blackwell}, address = {Malden}, issn = {0956-540X}, doi = {10.1111/j.1365-246X.2010.04863.x}, pages = {793 -- 800}, year = {2011}, abstract = {One of the key challenges in the context of local site effect studies is the determination of frequencies where the shakeability of the ground is enhanced. In this context, the H/V technique has become increasingly popular and peak frequencies of H/V spectral ratio are sometimes interpreted as resonance frequencies of the transmission response. In the present study, assuming that Rayleigh surface wave is dominant in H/V spectral ratio, we analyse theoretically under which conditions this may be justified and when not. We focus on 'layer over half-space' models which, although seemingly simple, capture many aspects of local site effects in real sedimentary structures. Our starting point is the ellipticity of Rayleigh waves. We use the exact formula of the H/V-ratio presented by Malischewsky \& Scherbaum (2004) to investigate the main characteristics of peak and trough frequencies. We present a simple formula illustrating if and where H/V-ratio curves have sharp peaks in dependence of model parameters. In addition, we have constructed a map, which demonstrates the relation between the H/V-peak frequency and the peak frequency of the transmission response in the domain of the layer's Poisson ratio and the impedance contrast. Finally, we have derived maps showing the relationship between the H/V-peak and trough frequency and key parameters of the model such as impedance contrast. These maps are seen as diagnostic tools, which can help to guide the interpretation of H/V spectral ratio diagrams in the context of site effect studies.}, language = {en} } @article{BoraScherbaumKuehnetal.2016, author = {Bora, Sanjay Singh and Scherbaum, Frank and Kuehn, Nicolas and Stafford, Peter}, title = {On the Relationship between Fourier and Response Spectra: Implications for the Adjustment of Empirical Ground-Motion Prediction Equations (GMPEs)}, series = {Bulletin of the Seismological Society of America}, volume = {106}, journal = {Bulletin of the Seismological Society of America}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0037-1106}, doi = {10.1785/0120150129}, pages = {1235 -- 1253}, year = {2016}, abstract = {The functional form of empirical response spectral ground-motion prediction equations (GMPEs) is often derived using concepts borrowed from Fourier spectral modeling of ground motion. As these GMPEs are subsequently calibrated with empirical observations, this may not appear to pose any major problems in the prediction of ground motion for a particular earthquake scenario. However, the assumption that Fourier spectral concepts persist for response spectra can lead to undesirable consequences when it comes to the adjustment of response spectral GMPEs to represent conditions not covered in the original empirical data set. In this context, a couple of important questions arise, for example, what are the distinctions and/or similarities between Fourier and response spectra of ground motions? And, if they are different, then what is the mechanism responsible for such differences and how do adjustments that are made to Fourier amplitude spectrum (FAS) manifest in response spectra? The present article explores the relationship between the Fourier and response spectrum of ground motion by using random vibration theory (RVT). With a simple Brune (1970, 1971) source model, RVT-generated acceleration spectra for a fixed magnitude and distance scenario are used. The RVT analyses reveal that the scaling of low oscillator-frequency response spectral ordinates can be treated as being equivalent to the scaling of the corresponding Fourier spectral ordinates. However, the high oscillator-frequency response spectral ordinates are controlled by a rather wide band of Fourier spectral ordinates. In fact, the peak ground acceleration, counter to the popular perception that it is a reflection of the high-frequency characteristics of ground motion, is controlled by the entire Fourier spectrum of ground motion. Additionally, this article demonstrates how an adjustment made to FAS is similar or different to the same adjustment made to response spectral ordinates. For this purpose, two cases: adjustments to the stress parameter (Delta sigma) (source term), and adjustments to the attributes reflecting site response (V-S - kappa(0)) are considered.}, language = {en} } @article{ScherbaumSchmedesCotton2004, author = {Scherbaum, Frank and Schmedes, J. and Cotton, Fabrice}, title = {On the conversion of source-to-site distance measures for extended earthquake source models}, issn = {0037-1106}, year = {2004}, abstract = {One of the major challenges in engineering seismology is the reliable prediction of site-specific ground motion for particular earthquakes, observed at specific distances. For larger events, a special problem arises, at short distances, with the source-to-site distance measure, because distance metrics based on a point-source model are no longer appropriate. As a consequence, different attenuation relations differ in the distance metric that they use. In addition to being a source of confusion, this causes problems to quantitatively compare or combine different ground- motion models; for example, in the context of Probabilistic Seismic Hazard Assessment, in cases where ground-motion models with different distance metrics occupy neighboring branches of a logic tree. In such a situation, very crude assumptions about source sizes and orientations often have to be used to be able to derive an estimate of the particular metric required. Even if this solves the problem of providing a number to put into the attenuation relation, a serious problem remains. When converting distance measures, the corresponding uncertainties map onto the estimated ground motions according to the laws of error propagation. To make matters worse, conversion of distance metrics can cause the uncertainties of the adapted ground-motion model to become magnitude and distance dependent, even if they are not in the original relation. To be able to treat this problem quantitatively, the variability increase caused by the distance metric conversion has to be quantified. For this purpose, we have used well established scaling laws to determine explicit distance conversion relations using regression analysis on simulated data. We demonstrate that, for all practical purposes, most popular distance metrics can be related to the Joyner-Boore distance using models based on gamma distributions to express the shape of some "residual function." The functional forms are magnitude and distance dependent and are expressed as polynomials. We compare the performance of these relations with manually derived individual distance estimates for the Landers, the Imperial Valley, and the Chi-Chi earthquakes}, language = {en} } @book{Scherbaum2001, author = {Scherbaum, Frank}, title = {Of poles and zeros : fundamentals of digital seismology}, series = {Modern approaches in geophysics}, volume = {15}, journal = {Modern approaches in geophysics}, edition = {Rev. 2. ed., reprint with corr}, publisher = {Springer}, address = {Dordrecht}, isbn = {0-7923-6834-7}, pages = {265 p.}, year = {2001}, language = {en} } @misc{ZaliReinKruegeretal.2023, author = {Zali, Zahra and Rein, Teresa and Kr{\"u}ger, Frank and Ohrnberger, Matthias and Scherbaum, Frank}, title = {Ocean bottom seismometer (OBS) noise reduction from horizontal and vertical components using harmonic-percussive separation algorithms}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1320}, issn = {1866-8372}, doi = {10.25932/publishup-58882}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-588828}, pages = {15}, year = {2023}, abstract = {Records from ocean bottom seismometers (OBSs) are highly contaminated by noise, which is much stronger compared to data from most land stations, especially on the horizontal components. As a consequence, the high energy of the oceanic noise at frequencies below 1 Hz considerably complicates the analysis of the teleseismic earthquake signals recorded by OBSs. Previous studies suggested different approaches to remove low-frequency noises from OBS recordings but mainly focused on the vertical component. The records of horizontal components, which are crucial for the application of many methods in passive seismological analysis of body and surface waves, could not be much improved in the teleseismic frequency band. Here we introduce a noise reduction method, which is derived from the harmonic-percussive separation algorithms used in Zali et al. (2021), in order to separate long-lasting narrowband signals from broadband transients in the OBS signal. This leads to significant noise reduction of OBS records on both the vertical and horizontal components and increases the earthquake signal-to-noise ratio (SNR) without distortion of the broadband earthquake waveforms. This is demonstrated through tests with synthetic data. Both SNR and cross-correlation coefficients showed significant improvements for different realistic noise realizations. The application of denoised signals in surface wave analysis and receiver functions is discussed through tests with synthetic and real data.}, language = {en} } @article{ZaliReinKruegeretal.2023, author = {Zali, Zahra and Rein, Teresa and Kr{\"u}ger, Frank and Ohrnberger, Matthias and Scherbaum, Frank}, title = {Ocean bottom seismometer (OBS) noise reduction from horizontal and vertical components using harmonic-percussive separation algorithms}, series = {Solid earth}, volume = {14}, journal = {Solid earth}, number = {2}, publisher = {Coepernicus Publ.}, address = {G{\"o}ttingen}, issn = {1869-9529}, doi = {10.5194/se-14-181-2023}, pages = {181 -- 195}, year = {2023}, abstract = {Records from ocean bottom seismometers (OBSs) are highly contaminated by noise, which is much stronger compared to data from most land stations, especially on the horizontal components. As a consequence, the high energy of the oceanic noise at frequencies below 1 Hz considerably complicates the analysis of the teleseismic earthquake signals recorded by OBSs. Previous studies suggested different approaches to remove low-frequency noises from OBS recordings but mainly focused on the vertical component. The records of horizontal components, which are crucial for the application of many methods in passive seismological analysis of body and surface waves, could not be much improved in the teleseismic frequency band. Here we introduce a noise reduction method, which is derived from the harmonic-percussive separation algorithms used in Zali et al. (2021), in order to separate long-lasting narrowband signals from broadband transients in the OBS signal. This leads to significant noise reduction of OBS records on both the vertical and horizontal components and increases the earthquake signal-to-noise ratio (SNR) without distortion of the broadband earthquake waveforms. This is demonstrated through tests with synthetic data. Both SNR and cross-correlation coefficients showed significant improvements for different realistic noise realizations. The application of denoised signals in surface wave analysis and receiver functions is discussed through tests with synthetic and real data.}, language = {en} } @article{BoraCottonScherbaum2019, author = {Bora, Sanjay Singh and Cotton, Fabrice and Scherbaum, Frank}, title = {NGA-West2 Empirical Fourier and Duration Models to Generate Adjustable Response Spectra}, series = {Earthquake spectra : the professional journal of the Earthquake Engineering Research Institute}, volume = {35}, journal = {Earthquake spectra : the professional journal of the Earthquake Engineering Research Institute}, number = {1}, publisher = {Sage Publ.}, address = {Thousand Oaks}, issn = {8755-2930}, doi = {10.1193/110317EQS228M}, pages = {61 -- 93}, year = {2019}, abstract = {Adjustment of median ground motion prediction equations (GMPEs) from one region to another region is one of the major challenges within the current practice of seismic hazard analysis. In our approach of generating response spectra, we derive two separate empirical models for a) Fourier amplitude spectrum (FAS) and b) duration of ground motion. To calculate response spectra, the two models are combined within the random vibration theory (RVT) framework. The models are calibrated on recordings obtained from shallow crustal earthquakes in active tectonic regions. We use a subset of NGA-West2 database with M3.2-7.9 earthquakes at distances 0-300 km. The NGA-West2 database expanded over a wide magnitude and distance range facilitates a better constraint over derived models. A frequency-dependent duration model is derived to obtain adjustable response spectral ordinates. Excellent comparison of our approach with other NGA-West2 models implies that it can also be used as a stand-alone model.}, language = {en} } @article{HiemerRoesslerScherbaum2012, author = {Hiemer, Stefan and R{\"o}ßler, Dirk and Scherbaum, Frank}, title = {Monitoring the West Bohemian earthquake swarm in 2008/2009 by a temporary small-aperture seismic array}, series = {Journal of seismology}, volume = {16}, journal = {Journal of seismology}, number = {2}, publisher = {Springer}, address = {Dordrecht}, issn = {1383-4649}, doi = {10.1007/s10950-011-9256-5}, pages = {169 -- 182}, year = {2012}, abstract = {The most recent intense earthquake swarm in West Bohemia lasted from 6 October 2008 to January 2009. Starting 12 days after the onset, the University of Potsdam monitored the swarm by a temporary small-aperture seismic array at 10 km epicentral distance. The purpose of the installation was a complete monitoring of the swarm including micro-earthquakes (M (L) < 0). We identify earthquakes using a conventional short-term average/long-term average trigger combined with sliding-window frequency-wavenumber and polarisation analyses. The resulting earthquake catalogue consists of 14,530 earthquakes between 19 October 2008 and 18 March 2009 with magnitudes in the range of -aEuro parts per thousand 1.2 a parts per thousand currency signaEuro parts per thousand M (L) a parts per thousand currency signaEuro parts per thousand 2.7. The small-aperture seismic array substantially lowers the detection threshold to about M (c) = -aEuro parts per thousand 0.4, when compared to the regional networks operating in West Bohemia (M (c) > 0.0). In the course of this work, the main temporal features (frequency-magnitude distribution, propagation of back azimuth and horizontal slowness, occurrence rate of aftershock sequences and interevent-time distribution) of the recent 2008/2009 earthquake swarm are presented and discussed. Temporal changes of the coefficient of variation (based on interevent times) suggest that the swarm earthquake activity of the 2008/2009 swarm terminates by 12 January 2009. During the main phase in our studied swarm period after 19 October, the b value of the Gutenberg-Richter relation decreases from 1.2 to 0.8. This trend is also reflected in the power-law behavior of the seismic moment release. The corresponding total seismic moment release of 1.02x10(17) Nm is equivalent to M (L,max) = 5.4.}, language = {en} } @article{SchmedesHainzlReameretal.2005, author = {Schmedes, J. and Hainzl, Sebastian and Reamer, S. K. and Scherbaum, Frank and Hinzen, K. G.}, title = {Moment release in the Lower Rhine Embayment, Germany : seismological perspective of the deformation process}, issn = {0956-540X}, year = {2005}, abstract = {An important task of seismic hazard assessment consists of estimating the rate of seismic moment release which is correlated to the rate of tectonic deformation and the seismic coupling. However, the estimations of deformation depend on the type of information utilized (e.g. geodetic, geological, seismic) and include large uncertainties. We therefore estimate the deformation rate in the Lower Rhine Embayment (LRE), Germany, using an integrated approach where the uncertainties have been systematically incorporated. On the basis of a new homogeneous earthquake catalogue we initially determine the frequency-magnitude distribution by statistical methods. In particular, we focus on an adequate estimation of the upper bound of the Gutenberg-Richter relation and demonstrate the importance of additional palaeoseis- mological information. The integration of seismological and geological information yields a probability distribution of the upper bound magnitude. Using this distribution together with the distribution of Gutenberg-Richter a and b values, we perform Monte Carlo simulations to derive the seismic moment release as a function of the observation time. The seismic moment release estimated from synthetic earthquake catalogues with short catalogue length is found to systematically underestimate the long-term moment rate which can be analytically determined. The moment release recorded in the LRE over the last 250 yr is found to be in good agreement with the probability distribution resulting from the Monte Carlo simulations. Furthermore, the long-term distribution is within its uncertainties consistent with the moment rate derived by geological measurements, indicating an almost complete seismic coupling in this region. By means of Kostrov's formula, we additionally calculate the full deformation rate tensor using the distribution of known focal mechanisms in LRE. Finally, we use the same approach to calculate the seismic moment and the deformation rate for two subsets of the catalogue corresponding to the east- and west-dipping faults, respectively}, language = {en} } @article{KuehnRiggelsenScherbaum2011, author = {K{\"u}hn, Nicolas M. and Riggelsen, Carsten and Scherbaum, Frank}, title = {Modeling the joint probability of earthquake, site, and ground-motion parameters using bayesian networks}, series = {Bulletin of the Seismological Society of America}, volume = {101}, journal = {Bulletin of the Seismological Society of America}, number = {1}, publisher = {Seismological Society of America}, address = {El Cerrito}, issn = {0037-1106}, doi = {10.1785/0120100080}, pages = {235 -- 249}, year = {2011}, abstract = {Bayesian networks are a powerful and increasingly popular tool for reasoning under uncertainty, offering intuitive insight into (probabilistic) data-generating processes. They have been successfully applied to many different fields, including bioinformatics. In this paper, Bayesian networks are used to model the joint-probability distribution of selected earthquake, site, and ground-motion parameters. This provides a probabilistic representation of the independencies and dependencies between these variables. In particular, contrary to classical regression, Bayesian networks do not distinguish between target and predictors, treating each variable as random variable. The capability of Bayesian networks to model the ground-motion domain in probabilistic seismic hazard analysis is shown for a generic situation. A Bayesian network is learned based on a subset of the Next Generation Attenuation (NGA) dataset, using 3342 records from 154 earthquakes. Because no prior assumptions about dependencies between particular parameters are made, the learned network displays the most probable model given the data. The learned network shows that the ground-motion parameter (horizontal peak ground acceleration, PGA) is directly connected only to the moment magnitude, Joyner-Boore distance, fault mechanism, source-to-site azimuth, and depth to a shear-wave horizon of 2: 5 km/s (Z2.5). In particular, the effect of V-S30 is mediated by Z2.5. Comparisons of the PGA distributions based on the Bayesian networks with the NGA model of Boore and Atkinson (2008) show a reasonable agreement in ranges of good data coverage.}, language = {en} } @article{KuleshHolschneiderDialloetal.2005, author = {Kulesh, Michail and Holschneider, Matthias and Diallo, Mamadou Sanou and Xie, Q. and Scherbaum, Frank}, title = {Modeling of wave dispersion using continuous wavelet transforms}, issn = {0033-4553}, year = {2005}, abstract = {In the estimate of dispersion with the help of wavelet analysis considerable emphasis has been put on the extraction of the group velocity using the modulus of the wavelet transform. In this paper we give an asymptotic expression of the full propagator in wavelet space that comprises the phase velocity as well. This operator establishes a relationship between the observed signals at two different stations during wave propagation in a dispersive and attenuating medium. Numerical and experimental examples are presented to show that the method accurately models seismic wave dispersion and attenuation}, language = {en} } @article{ScherbaumDelavaudRiggelsen2009, author = {Scherbaum, Frank and Delavaud, Elise and Riggelsen, Carsten}, title = {Model selection in seismic hazard analysis : an information-theoretic perspective}, issn = {0037-1106}, doi = {10.1785/0120080347}, year = {2009}, abstract = {Although the methodological framework of probabilistic seismic hazard analysis is well established, the selection of models to predict the ground motion at the sites of interest remains a major challenge. Information theory provides a powerful theoretical framework that can guide this selection process in a consistent way. From an information- theoretic perspective, the appropriateness of models can be expressed in terms of their relative information loss (Kullback-Leibler distance) and hence in physically meaningful units (bits). In contrast to hypothesis testing, information-theoretic model selection does not require ad hoc decisions regarding significance levels nor does it require the models to be mutually exclusive and collectively exhaustive. The key ingredient, the Kullback-Leibler distance, can be estimated from the statistical expectation of log-likelihoods of observations for the models under consideration. In the present study, data-driven ground-motion model selection based on Kullback-Leibler-distance differences is illustrated for a set of simulated observations of response spectra and macroseismic intensities. Information theory allows for a unified treatment of both quantities. The application of Kullback-Leibler-distance based model selection to real data using the model generating data set for the Abrahamson and Silva (1997) ground-motion model demonstrates the superior performance of the information-theoretic perspective in comparison to earlier attempts at data- driven model selection (e.g., Scherbaum et al., 2004).}, language = {en} } @article{HaendelvonSpechtKuehnetal.2015, author = {H{\"a}ndel, Annabel and von Specht, Sebastian and Kuehn, Nicolas M. and Scherbaum, Frank}, title = {Mixtures of ground-motion prediction equations as backbone models for a logic tree: an application to the subduction zone in Northern Chile}, series = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering}, volume = {13}, journal = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering}, number = {2}, publisher = {Springer}, address = {Dordrecht}, issn = {1570-761X}, doi = {10.1007/s10518-014-9636-7}, pages = {483 -- 501}, year = {2015}, abstract = {In probabilistic seismic hazard analysis, different ground-motion prediction equations (GMPEs) are commonly combined within a logic tree framework. The selection of appropriate GMPEs, however, is a non-trivial task, especially for regions where strong motion data are sparse and where no indigenous GMPE exists because the set of models needs to capture the whole range of ground-motion uncertainty. In this study we investigate the aggregation of GMPEs into a mixture model with the aim to infer a backbone model that is able to represent the center of the ground-motion distribution in a logic tree analysis. This central model can be scaled up and down to obtain the full range of ground-motion uncertainty. The combination of models into a mixture is inferred from observed ground-motion data. We tested the new approach for Northern Chile, a region for which no indigenous GMPE exists. Mixture models were calculated for interface and intraslab type events individually. For each source type we aggregated eight subduction zone GMPEs using mainly new strong-motion data that were recorded within the Plate Boundary Observatory Chile project and that were processed within this study. We can show that the mixture performs better than any of its component GMPEs, and that it performs comparable to a regression model that was derived for the same dataset. The mixture model seems to represent the median ground motions in that region fairly well. It is thus able to serve as a backbone model for the logic tree.}, language = {en} } @article{GianniotisKuehnScherbaum2014, author = {Gianniotis, Nikolaos and Kuehn, Nicolas and Scherbaum, Frank}, title = {Manifold aligned ground motion prediction equations for regional datasets}, series = {Computers \& geosciences : an international journal devoted to the publication of papers on all aspects of geocomputation and to the distribution of computer programs and test data sets ; an official journal of the International Association for Mathematical Geology}, volume = {69}, journal = {Computers \& geosciences : an international journal devoted to the publication of papers on all aspects of geocomputation and to the distribution of computer programs and test data sets ; an official journal of the International Association for Mathematical Geology}, publisher = {Elsevier}, address = {Oxford}, issn = {0098-3004}, doi = {10.1016/j.cageo.2014.04.014}, pages = {72 -- 77}, year = {2014}, abstract = {Inferring a ground-motion prediction equation (GMPE) for a region in which only a small number of seismic events has been observed is a challenging task. A response to this data scarcity is to utilise data from other regions in the hope that there exist common patterns in the generation of ground motion that can contribute to the development of a GMPE for the region in question. This is not an unreasonable course of action since we expect regional GMPEs to be related to each other. In this work we model this relatedness by assuming that the regional GMPEs occupy a common low-dimensional manifold in the space of all possible GMPEs. As a consequence, the GMPEs are fitted in a joint manner and not independent of each other, borrowing predictive strength from each other's regional datasets. Experimentation on a real dataset shows that the manifold assumption displays better predictive performance over fitting regional GMPEs independent of each other. (C) 2014 Elsevier Ltd. All rights reserved.}, language = {en} } @article{HaneyKummerowLangenbruchetal.2011, author = {Haney, Frank and Kummerow, J. and Langenbruch, C. and Dinske, C. and Shapiro, Serge A. and Scherbaum, Frank}, title = {Magnitude estimation for microseismicity induced during the KTB 2004/2005 injection experiment}, series = {Geophysics}, volume = {76}, journal = {Geophysics}, number = {6}, publisher = {Society of Exploration Geophysicists}, address = {Tulsa}, issn = {0016-8033}, doi = {10.1190/GEO2011-0020.1}, pages = {WC47 -- WC53}, year = {2011}, abstract = {We determined the magnitudes of 2540 microseismic events measured at one single 3C borehole geophone at the German Deep Drilling Site (known by the German acronym, KTB) during the injection phase 2004/2005. For this task we developed a three-step approach. First, we estimated local magnitudes of 104 larger events with a standard method based on amplitude measurements at near-surface stations. Second, we investigated a series of parameters to characterize the size of these events using the seismograms of the borehole sensor, and we compared them statistically with the local magnitudes. Third, we extrapolated the regression curve to obtain the magnitudes of 2436 events that were only measured at the borehole geophone. This method improved the magnitude of completeness for the KTB data set by more than one order down to M = -2.75. The resulting b-value for all events was 0.78, which is similar to the b-value obtained from taking only the greater events with standard local magnitude estimation from near-surface stations, b = 0.86. The more complete magnitude catalog was required to study the magnitude distribution with time and to characterize the seismotectonic state of the KTB injection site. The event distribution with time was consistent with prediction from theory assuming pore pressure diffusion as the underlying mechanism to trigger the events. The value we obtained for the seismogenic index of -4 suggested that the seismic hazard potential at the KTB site is comparatively low.}, language = {en} } @article{MalischewskyScherbaum2004, author = {Malischewsky, Peter G. and Scherbaum, Frank}, title = {Love's formula and H/V-ratio (ellipticity) of Rayleigh waves}, issn = {0165-2125}, year = {2004}, abstract = {The ellipticity of Rayleigh surface waves, which is an important parameter characterizing the propagation medium, is studied for several models with increasing complexity. While the main focus lies on theory, practical implications of the use of the horizontal to vertical component ratio (H/V-ratio) to Study the subsurface structure are considered as well. Love's approximation of the ellipticity for an incompressible layer over an incompressible half-space is critically discussed especially concerning its applicability for different impedance contrasts. The main result is an analytically exact formula of H/V for a 2-layer model of compressible media, which is a generalization of Love's formula. It turns out that for a limited range of models Love's approximation can be used also in the general case. (C) 2003 Elsevier B.V. All rights reserved}, language = {en} } @article{ScherbaumKuehn2011, author = {Scherbaum, Frank and K{\"u}hn, Nicolas M.}, title = {Logic tree branch weights and probabilities summing up to one is not enough}, series = {Earthquake spectra : the professional journal of the Earthquake Engineering Research Institute}, volume = {27}, journal = {Earthquake spectra : the professional journal of the Earthquake Engineering Research Institute}, number = {4}, publisher = {Earthquake Engineering Research Institute}, address = {Oakland}, issn = {8755-2930}, doi = {10.1193/1.3652744}, pages = {1237 -- 1251}, year = {2011}, abstract = {Logic trees have become the most popular tool for the quantification of epistemic uncertainties in probabilistic seismic hazard assessment (PSHA). In a logic-tree framework, epistemic uncertainty is expressed in a set of branch weights, by which an expert or an expert group assigns degree-of-belief values to the applicability of the corresponding branch models. Despite the popularity of logic-trees, however, one finds surprisingly few clear commitments to what logic-tree branch weights are assumed to be (even by hazard analysts designing logic trees). In the present paper we argue that it is important for hazard analysts to accept the probabilistic framework from the beginning for assigning logic-tree branch weights. In other words, to accept that logic-tree branch weights are probabilities in the axiomatic sense, independent of one's preference for the philosophical interpretation of probabilities. We demonstrate that interpreting logic-tree branch weights merely as a numerical measure of "model quality," which are then subsequently normalized to sum up to unity, will with increasing number of models inevitably lead to an apparent insensitivity of hazard curves on the logic-tree branch weights, which may even be mistaken for robustness of the results. Finally, we argue that assigning logic-tree branch weights in a sequential fashion may improve their logical consistency.}, language = {en} } @article{WassermannOhrnbergerScherbaumetal.1998, author = {Wassermann, Joachim and Ohrnberger, Matthias and Scherbaum, Frank and Gossler, J. and Zschau, Jochen}, title = {Kontinuierliche seismologische Netz- und Arraymessungen am Dekadenvulkan Merapi (Java, Indonesien) : ein Zwischenres{\"u}mee = Continuous measurements at Merapi volcano (Java, Indonesia) using anetwork of small-scale seismograph arrays}, issn = {0947-1944}, year = {1998}, language = {de} } @article{DialloKuleshHolschneideretal.2006, author = {Diallo, Mamadou Sanou and Kulesh, Michail and Holschneider, Matthias and Kurennaya, Kristina and Scherbaum, Frank}, title = {Instantaneous polarization attributes based on an adaptive approximate covariance method}, series = {Geophysics}, volume = {71}, journal = {Geophysics}, number = {5}, publisher = {SEG}, address = {Tulsa}, issn = {0016-8033}, doi = {10.1190/1.2227522}, pages = {V99 -- V104}, year = {2006}, abstract = {We introduce a method for computing instantaneous-polarization attributes from multicomponent signals. This is an improvement on the standard covariance method (SCM) because it does not depend on the window size used to compute the standard covariance matrix. We overcome the window-size problem by deriving an approximate analytical formula for the cross-energy matrix in which we automatically and adaptively determine the time window. The proposed method uses polarization analysis as applied to multicomponent seismic by waveform separation and filtering.}, language = {en} } @book{LoosScherbaum1999, author = {Loos, Wolfgang and Scherbaum, Frank}, title = {Inner earth : a seismosonic symphony}, publisher = {Traumton (Indigo Vertrieb)}, address = {[s.l.]}, pages = {Audio-CD}, year = {1999}, language = {en} } @article{DelavaudScherbaumKuehnetal.2009, author = {Delavaud, Elise and Scherbaum, Frank and Kuehn, Nicolas and Riggelsen, Carsten}, title = {Information-theoretic selection of ground-motion prediction equations for seismic hazard analysis : an applicability study using Californian data}, issn = {0037-1106}, doi = {10.1785/0120090055}, year = {2009}, abstract = {Considering the increasing number and complexity of ground-motion prediction equations available for seismic hazard assessment, there is a definite need for an efficient, quantitative, and robust method to select and rank these models for a particular region of interest. In a recent article, Scherbaum et al. (2009) have suggested an information- theoretic approach for this purpose that overcomes several shortcomings of earlier attempts at using data-driven ground- motion prediction equation selection procedures. The results of their theoretical study provides evidence that in addition to observed response spectra, macroseismic intensity data might be useful for model selection and ranking. We present here an applicability study for this approach using response spectra and macroseismic intensities from eight Californian earthquakes. A total of 17 ground-motion prediction equations, from different regions, for response spectra, combined with the equation of Atkinson and Kaka (2007) for macroseismic intensities are tested for their relative performance. The resulting data-driven rankings show that the models that best estimate ground motion in California are, as one would expect, Californian and western U. S. models, while some European models also perform fairly well. Moreover, the model performance appears to be strongly dependent on both distance and frequency. The relative information of intensity versus response spectral data is also explored. The strong correlation we obtain between intensity-based rankings and spectral-based ones demonstrates the great potential of macroseismic intensities data for model selection in the context of seismic hazard assessment.}, language = {en} }