@article{BoraCottonScherbaumetal.2017,
  author    = {Bora, Sanjay Singh and Cotton, Fabrice Pierre 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{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}
}
@article{ScherbaumCottonStaedtke2006,
  author    = {Scherbaum, Frank and Cotton, Fabrice Pierre 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{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{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{DouglasBungumScherbaum2006,
  author    = {Douglas, John and Bungum, Hilmar and Scherbaum, Frank},
  title     = {Ground-motion prediction equations for Southern Spain and Southern Norway obtained using the composite model perspective},
  issn      = {1363-2469},
  doi       = {10.1080/1363246060935058},
  year      = {2006},
  abstract  = {In this paper, two sets of earthquake ground-motion relations to estimate peak ground and response spectral acceleration are developed for sites in southern Spain and in southern Norway using a recently published composite approach. For this purpose seven empirical ground-motion relations developed from recorded strong-motion data from different parts of the world were employed. The different relations were first adjusted based on a number of transformations to convert the differing choices of independent parameters to a single one. After these transformations, which include the scatter introduced, were performed, the equations were modified to account for differences between the host and the target regions using the stochastic method to compute the host-to-target conversion factors. Finally functions were fitted to the derived ground-motion estimates to obtain sets of seven individual equations for use in probabilistic seismic hazard assessment for southern Spain and southern Norway. The relations are compared with local ones published for the two regions. The composite methodology calls for the setting up of independent logic trees for the median values and for the sigma values, in order to properly separate epistemic and aleatory uncertainties after the corrections and the conversions},
  language  = {en}
}
@article{EwaldIgelHinzenetal.2006,
  author    = {Ewald, Michael and Igel, Heiner and Hinzen, Klaus-G{\"u}nther and Scherbaum, Frank},
  title     = {Basin-related effects on ground motion for earthquake scenarios in the Lower Rhine Embayment},
  issn      = {0956-540X},
  doi       = {10.1111/j.1365-246X.2006.02909.x},
  year      = {2006},
  abstract  = {The deterministic calculation of earthquake scenarios using complete waveform modelling plays an increasingly important role in estimating shaking hazard in seismically active regions. Here we apply 3-D numerical modelling of seismic wave propagation to M 6+ earthquake scenarios in the area of the Lower Rhine Embayment, one of the seismically most active regions in central Europe. Using a 3-D basin model derived from geology, borehole information and seismic experiments, we aim at demonstrating the strong dependence of ground shaking on hypocentre location and basin structure. The simulations are carried out up to frequencies of ca. 1 Hz. As expected, the basin structure leads to strong lateral variations in peak ground motion, amplification and shaking duration. Depending on source-basin-receiver geometry, the effects correlate with basin depth and the slope of the basin flanks; yet, the basin also affects peak ground motion and estimated shaking hazard thereof outside the basin. Comparison with measured seismograms for one of the earthquakes shows that some of the main characteristics of the wave motion are reproduced. Cumulating the derived seismic intensities from the three modelled earthquake scenarios leads to a predominantly basin correlated intensity distribution for our study area},
  language  = {en}
}
@article{MussonToroCoppersmithetal.2005,
  author    = {Musson, R. M. W. and Toro, G. R. and Coppersmith, Kevin J. and Bommer, Julian J. and Deichmann, N. and Bungum, Hilmar and Cotton, Fabrice Pierre and Scherbaum, Frank and Slejko, Dario and Abrahamson, Norman A.},
  title     = {Evaluating hazard results for Switzerland and how not to do it : a discussion of "Problems in the application of the SSHAC probability method for assessing earthquake hazards at Swiss nuclear power plants" by J-U Klugel},
  year      = {2005},
  abstract  = {The PEGASOS project was a major international seismic hazard study, one of the largest ever conducted anywhere in the world, to assess seismic hazard at four nuclear power plant sites in Switzerland. Before the report of this project has become publicly available, a paper attacking both methodology and results has appeared. Since the general scientific readership may have difficulty in assessing this attack in the absence of the report being attacked, we supply a response in the present paper. The bulk of the attack, besides some misconceived arguments about the role of uncertainties in seismic hazard analysis, is carried by some exercises that purport to be validation exercises. In practice, they are no such thing; they are merely independent sets of hazard calculations based on varying assumptions and procedures, often rather questionable, which come up with various different answers which have no particular significance. (C) 2005 Elsevier B.V. All rights reserved},
  language  = {en}
}
@article{ScherbaumBommerBungumetal.2005,
  author    = {Scherbaum, Frank and Bommer, Julian J. and Bungum, Hilmar and Cotton, Fabrice Pierre and Abrahamson, Norman A.},
  title     = {Composite ground-motion models and logic trees: Methodology, sensitivities, and uncertainties},
  issn      = {0037-1106},
  year      = {2005},
  abstract  = {Logic trees have become a popular tool in seismic hazard studies. Commonly, the models corresponding to the end branches of the complete logic tree in a probabalistic seismic hazard analysis (PSHA) are treated separately until the final calculation of the set of hazard curves. This comes at the price that information regarding sensitivities and uncertainties in the ground-motion sections of the logic tree are only obtainable after disaggregation. Furthermore, from this end-branch model perspective even the designers of the logic tree cannot directly tell what ground-motion scenarios most likely would result from their logic trees for a given earthquake at a particular distance, nor how uncertain these scenarios might be or how they would be affected by the choices of the hazard analyst. On the other hand, all this information is already implicitly present in the logic tree. Therefore, with the ground-motion perspective that we propose in the present article, we treat the ground-motion sections of a complete logic tree for seismic hazard as a single composite model representing the complete state-of-knowledge-and-belief of a particular analyst on ground motion in a particular target region. We implement this view by resampling the ground-motion models represented in the ground-motion sections of the logic tree by Monte Carlo simulation (separately for the median values and the sigma values) and then recombining the sets of simulated values in proportion to their logic-tree branch weights. The quantiles of this resampled composite model provide the hazard analyst and the decision maker with a simple, clear, and quantitative representation of the overall physical meaning of the ground-motion section of a logic tree and the accompanying epistemic uncertainty. Quantiles of the composite model also provide an easy way to analyze the sensitivities and uncertainties related to a given logic-tree model. We illustrate this for a composite ground- motion model for central Europe. Further potential fields of applications are seen wherever individual best estimates of ground motion have to be derived from a set of candidate models, for example, for hazard rnaps, sensitivity studies, or for modeling scenario earthquakes},
  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}
}