@article{KtenidouRoumeliotiAbrahamsonetal.2018, author = {Ktenidou, Olga-Joan and Roumelioti, Zafeiria and Abrahamson, Norman and Cotton, Fabrice and Pitilakis, Kyriazis and Hollender, Fabrice}, title = {Understanding single-station ground motion variability and uncertainty (sigma)}, series = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering}, volume = {16}, journal = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering}, number = {6}, publisher = {Springer}, address = {Dordrecht}, issn = {1570-761X}, doi = {10.1007/s10518-017-0098-6}, pages = {2311 -- 2336}, year = {2018}, abstract = {Accelerometric data from the well-studied valley EUROSEISTEST are used to investigate ground motion uncertainty and variability. We define a simple local ground motion prediction equation (GMPE) and investigate changes in standard deviation (σ) and its components, the between-event variability (τ) and within-event variability (φ). Improving seismological metadata significantly reduces τ (30-50\%), which in turn reduces the total σ. Improving site information reduces the systematic site-to-site variability, φ S2S (20-30\%), in turn reducing φ, and ultimately, σ. Our values of standard deviations are lower than global values from literature, and closer to path-specific than site-specific values. However, our data have insufficient azimuthal coverage for single-path analysis. Certain stations have higher ground-motion variability, possibly due to topography, basin edge or downgoing wave effects. Sensitivity checks show that 3 recordings per event is a sufficient data selection criterion, however, one of the dataset's advantages is the large number of recordings per station (9-90) that yields good site term estimates. We examine uncertainty components binning our data with magnitude from 0.01 to 2 s; at smaller magnitudes, τ decreases and φ SS increases, possibly due to κ and source-site trade-offs Finally, we investigate the alternative approach of computing φ SS using existing GMPEs instead of creating an ad hoc local GMPE. This is important where data are insufficient to create one, or when site-specific PSHA is performed. We show that global GMPEs may still capture φ SS , provided that: (1) the magnitude scaling errors are accommodated by the event terms; (2) there are no distance scaling errors (use of a regionally applicable model). Site terms (φ S2S ) computed by different global GMPEs (using different site-proxies) vary significantly, especially for hard-rock sites. This indicates that GMPEs may be poorly constrained where they are sometimes most needed, i.e., for hard rock.}, 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{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{ZhuCottonKawaseetal.2022, author = {Zhu, Chuanbin and Cotton, Fabrice and Kawase, Hiroshi and H{\"a}ndel, Annabel and Pilz, Marco and Nakano, Kenichi}, title = {How well can we predict earthquake site response so far?}, series = {Earthquake spectra : the professional journal of the Earthquake Engineering Research Institute}, volume = {38}, journal = {Earthquake spectra : the professional journal of the Earthquake Engineering Research Institute}, number = {2}, publisher = {Sage Publ.}, address = {Thousand Oaks}, issn = {8755-2930}, doi = {10.1177/87552930211060859}, pages = {1047 -- 1075}, year = {2022}, abstract = {Earthquake site responses or site effects are the modifications of surface geology to seismic waves. How well can we predict the site effects (average over many earthquakes) at individual sites so far? To address this question, we tested and compared the effectiveness of different estimation techniques in predicting the outcrop Fourier site responses separated using the general inversion technique (GIT) from recordings. Techniques being evaluated are (a) the empirical correction to the horizontal-to-vertical spectral ratio of earthquakes (c-HVSR), (b) one-dimensional ground response analysis (GRA), and (c) the square-root-impedance (SRI) method (also called the quarter-wavelength approach). Our results show that c-HVSR can capture significantly more site-specific features in site responses than both GRA and SRI in the aggregate, especially at relatively high frequencies. c-HVSR achieves a "good match" in spectral shape at similar to 80\%-90\% of 145 testing sites, whereas GRA and SRI fail at most sites. GRA and SRI results have a high level of parametric and/or modeling errors which can be constrained, to some extent, by collecting on-site recordings.}, language = {en} } @article{KothaCottonBindi2018, author = {Kotha, Sreeram Reddy and Cotton, Fabrice and Bindi, Dino}, title = {A new approach to site classification}, series = {Soil Dynamics and Earthquake Engineering}, volume = {110}, journal = {Soil Dynamics and Earthquake Engineering}, publisher = {Elsevier}, address = {Oxford}, issn = {0267-7261}, doi = {10.1016/j.soildyn.2018.01.051}, pages = {318 -- 329}, year = {2018}, abstract = {With increasing amount of strong motion data, Ground Motion Prediction Equation (GMPE) developers are able to quantify empirical site amplification functions (delta S2S(s)) from GMPE residuals, for use in site-specific Probabilistic Seismic Hazard Assessment. In this study, we first derive a GMPE for 5\% damped Pseudo Spectral Acceleration (g) of Active Shallow Crustal earthquakes in Japan with 3.4 <= M-w <= 7.3 and 0 <= R-JB <= 600km. Using k-mean spectral clustering technique, we then classify our estimated delta S2S(s)(T = 0.01 - 2s) of 588 wellcharacterized sites, into 8 site clusters with distinct mean site amplification functions, and within-cluster site-tosite variability similar to 50\% smaller than the overall dataset variability (phi(S2S)). Following an evaluation of existing schemes, we propose a revised data-driven site classification characterized by kernel density distributions of V-s30, V-s10, H-800, and predominant period (T-G) of the site clusters.}, language = {en} } @article{vonSpechtHeidbachCottonetal.2018, author = {von Specht, Sebastian and Heidbach, Oliver and Cotton, Fabrice and Zang, Arno}, title = {Uncertainty reduction of stress tensor inversion with data-driven catalogue selection}, series = {Geophysical journal international}, volume = {214}, journal = {Geophysical journal international}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggy240}, pages = {2250 -- 2263}, year = {2018}, abstract = {The selection of earthquake focal mechanisms (FMs) for stress tensor inversion (STI) is commonly done on a spatial basis, that is, hypocentres. However, this selection approach may include data that are undesired, for example, by mixing events that are caused by different stress tensors when for the STI a single stress tensor is assumed. Due to the significant increase of FM data in the past decades, objective data-driven data selection is feasible, allowing more refined FM catalogues that avoid these issues and provide data weights for the STI routines. We present the application of angular classification with expectation-maximization (ACE) as a tool for data selection. ACE identifies clusters of FM without a priori information. The identified clusters can be used for the classification of the style-of-faulting and as weights of the FM data. We demonstrate that ACE effectively selects data that can be associated with a single stress tensor. Two application examples are given for weighted STI from South America. We use the resulting clusters and weights as a priori information for an STI for these regions and show that uncertainties of the stress tensor estimates are reduced significantly.}, language = {en} } @article{LetortRetailleauBoueetal.2018, author = {Letort, Jean and Retailleau, Lise and Boue, Pierre and Radiguet, Mathilde and Gardonio, Blandine and Cotton, Fabrice and Campillo, Michel}, title = {Lateral variations of the Guerrero-Oaxaca subduction zone (Mexico) derived from weak seismicity (M(b)3.5+) detected on a single array at teleseismic distance}, series = {Geophysical journal international}, volume = {213}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggy035}, pages = {1002 -- 1012}, year = {2018}, abstract = {Detections of pP and sP phase arrivals (the so-called depth phases) at teleseismic distance provide one of the best ways to estimate earthquake focal depth, as the P-pP and the P-sP delays are strongly dependent on the depth. Based on a new processing workflow and using a single seismic array at teleseismic distance, we can estimate the depth of clusters of small events down to magnitude M-b 3.5. Our method provides a direct view of the relative variations of the seismicity depth from an active area. This study focuses on the application of this new methodology to study the lateral variations of the Guerrero subduction zone (Mexico) using the Eielson seismic array in Alaska (USA). After denoising the signals, 1232 M-b 3.5 + events were detected, with clear P, pP, sP and PcP arrivals. A high-resolution view of the lateral variations of the depth of the seismicity of the Guerero-Oaxaca area is thus obtained. The seismicity is shown to be mainly clustered along the interface, coherently following the geometry of the plate as constrained by the receiver-function analysis along the Meso America Subduction Experiment profile. From this study, the hypothesis of tears on the western part of Guerrero and the eastern part of Oaxaca are strongly confirmed by dramatic lateral changes in the depth of the earthquake clusters. The presence of these two tears might explain the observed lateral variations in seismicity, which is correlated with the boundaries of the slow slip events.}, language = {en} } @article{SocquetValdesJaraetal.2017, author = {Socquet, Anne and Valdes, Jesus Pina and Jara, Jorge and Cotton, Fabrice and Walpersdorf, Andrea and Cotte, Nathalie and von Specht, Sebastian and Ortega-Culaciati, Francisco and Carrizo, Daniel and Norabuena, Edmundo}, title = {An 8month slow slip event triggers progressive nucleation of the 2014 Chile megathrust}, series = {Geophysical research letters}, volume = {44}, journal = {Geophysical research letters}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1002/2017GL073023}, pages = {4046 -- 4053}, year = {2017}, abstract = {The mechanisms leading to large earthquakes are poorly understood and documented. Here we characterize the long-term precursory phase of the 1 April 2014 M(w)8.1 North Chile megathrust. We show that a group of coastal GPS stations accelerated westward 8months before the main shock, corresponding to a M(w)6.5 slow slip event on the subduction interface, 80\% of which was aseismic. Concurrent interface foreshocks underwent a diminution of their radiation at high frequency, as shown by the temporal evolution of Fourier spectra and residuals with respect to ground motions predicted by recent subduction models. Such ground motions change suggests that in response to the slow sliding of the subduction interface, seismic ruptures are progressively becoming smoother and/or slower. The gradual propagation of seismic ruptures beyond seismic asperities into surrounding metastable areas could explain these observations and might be the precursory mechanism eventually leading to the main shock.}, language = {en} } @article{KothaWeatherillBindietal.2022, author = {Kotha, Sreeram Reddy and Weatherill, Graeme and Bindi, Dino and Cotton, Fabrice}, title = {Near-source magnitude scaling of spectral accelerations}, series = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering}, volume = {20}, journal = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering}, number = {3}, publisher = {Springer}, address = {Dordrecht}, issn = {1570-761X}, doi = {10.1007/s10518-021-01308-5}, pages = {1343 -- 1370}, year = {2022}, abstract = {Ground-motion models (GMMs) are often used to predict the random distribution of Spectral accelerations (SAs) at a site due to a nearby earthquake. In probabilistic seismic hazard and risk assessment, large earthquakes occurring close to a site are considered as critical scenarios. GMMs are expected to predict realistic SAs with low within-model uncertainty (sigma(mu)) for such rare scenarios. However, the datasets used to regress GMMs are usually deficient of data from critical scenarios. The (Kotha et al., A Regionally Adaptable Ground-Motion Model for Shallow Crustal Earthquakes in Europe Bulletin of Earthquake Engineering 18:4091-4125, 2020) GMM developed from the Engineering strong motion (ESM) dataset was found to predict decreasing short-period SAs with increasing M-W >= M-h = 6.2, and with large sigma(mu) at near-source distances <= 30km. In this study, we updated the parametrisation of the GMM based on analyses of ESM and the Near source strong motion (NESS) datasets. With M-h = 5.7, we could rectify the M-W scaling issue, while also reducing sigma(mu). at M-W >= M-h. We then evaluated the GMM against NESS data, and found that the SAs from a few large, thrust-faulting events in California, New Zealand, Japan, and Mexico are significantly higher than GMM median predictions. However, recordings from these events were mostly made on soft-soil geology, and contain anisotropic pulse-like effects. A more thorough non-ergodic treatment of NESS was not possible because most sites sampled unique events in very diverse tectonic environments. We provide an updated set of GMM coefficients,sigma(mu), and heteroscedastic variance models; while also cautioning against its application for M-W <= 4 in low-moderate seismicity regions without evaluating the homogeneity of M-W estimates between pan-European ESM and regional datasets.}, 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} } @article{vonSpechtOeztuerkVehetal.2019, author = {von Specht, Sebastian and {\"O}zt{\"u}rk, Ugur and Veh, Georg and Cotton, Fabrice and Korup, Oliver}, title = {Effects of finite source rupture on landslide triggering}, series = {Solid earth}, volume = {10}, journal = {Solid earth}, number = {2}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1869-9510}, doi = {10.5194/se-10-463-2019}, pages = {463 -- 486}, year = {2019}, abstract = {The propagation of a seismic rupture on a fault introduces spatial variations in the seismic wave field surrounding the fault. This directivity effect results in larger shaking amplitudes in the rupture propagation direction. Its seismic radiation pattern also causes amplitude variations between the strike-normal and strike-parallel components of horizontal ground motion. We investigated the landslide response to these effects during the 2016 Kumamoto earthquake (M-w 7.1) in central Kyushu (Japan). Although the distribution of some 1500 earthquake-triggered landslides as a function of rupture distance is consistent with the observed Arias intensity, the landslides were more concentrated to the northeast of the southwest-northeast striking rupture. We examined several landslide susceptibility factors: hillslope inclination, the median amplification factor (MAF) of ground shaking, lithology, land cover, and topographic wetness. None of these factors sufficiently explains the landslide distribution or orientation (aspect), although the landslide head scarps have an elevated hillslope inclination and MAF. We propose a new physics-based ground-motion model (GMM) that accounts for the seismic rupture effects, and we demonstrate that the low-frequency seismic radiation pattern is consistent with the overall landslide distribution. Its spatial pattern is influenced by the rupture directivity effect, whereas landslide aspect is influenced by amplitude variations between the fault-normal and fault-parallel motion at frequencies < 2 Hz. This azimuth dependence implies that comparable landslide concentrations can occur at different distances from the rupture. This quantitative link between the prevalent landslide aspect and the low-frequency seismic radiation pattern can improve coseismic landslide hazard assessment.}, language = {en} } @article{PenaMetzgerHeidbachetal.2022, author = {Pe{\~n}a, Carlos and Metzger, Sabrina and Heidbach, Oliver and Bedford, Jonathan and Bookhagen, Bodo and Moreno, Marcos and Oncken, Onno and Cotton, Fabrice}, title = {Role of poroelasticity during the early postseismic deformation of the 2010 Maule megathrust earthquake}, series = {Geophysical research letters}, volume = {49}, journal = {Geophysical research letters}, number = {9}, publisher = {Wiley}, address = {Hoboken, NJ}, issn = {0094-8276}, doi = {10.1029/2022GL098144}, pages = {11}, year = {2022}, abstract = {Megathrust earthquakes impose changes of differential stress and pore pressure in the lithosphere-asthenosphere system that are transiently relaxed during the postseismic period primarily due to afterslip, viscoelastic and poroelastic processes. Especially during the early postseismic phase, however, the relative contribution of these processes to the observed surface deformation is unclear. To investigate this, we use geodetic data collected in the first 48 days following the 2010 Maule earthquake and a poro-viscoelastic forward model combined with an afterslip inversion. This model approach fits the geodetic data 14\% better than a pure elastic model. Particularly near the region of maximum coseismic slip, the predicted surface poroelastic uplift pattern explains well the observations. If poroelasticity is neglected, the spatial afterslip distribution is locally altered by up to +/- 40\%. Moreover, we find that shallow crustal aftershocks mostly occur in regions of increased postseismic pore-pressure changes, indicating that both processes might be mechanically coupled.}, language = {en} } @article{ZhuCottonKwaketal.2021, author = {Zhu, Chuanbin and Cotton, Fabrice and Kwak, Dong-Youp and Ji, Kun and Kawase, Hiroshi and Pilz, Marco}, title = {Within-site variability in earthquake site response}, series = {Geophysical journal international}, volume = {229}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggab481}, pages = {1268 -- 1281}, year = {2021}, abstract = {The within-site variability in site response is the randomness in site response at a given site from different earthquakes and is treated as aleatory variability in current seismic hazard/risk analyses. In this study, we investigate the single-station variability in linear site response at K-NET and KiK-net stations in Japan using a large number of earthquake recordings. We found that the standard deviation of the horizontal-to-vertical Fourier spectral ratio at individual sites, that is single-station horizontal-to-vertical spectral ratio (HVSR) sigma sigma(HV,s), approximates the within-site variability in site response quantified using surface-to-borehole spectral ratios (for oscillator frequencies higher than the site fundamental frequency) or empirical ground-motion models. Based on this finding, we then utilize the single-station HVSR sigma as a convenient tool to study the site-response variability at 697 KiK-net and 1169 K-NET sites. Our results show that at certain frequencies, stiff, rough and shallow sites, as well as small and local events tend to have a higher sigma(HV,s). However, when being averaged over different sites, the single-station HVSR sigma, that is sigma(HV), increases gradually with decreasing frequency. In the frequency range of 0.25-25 Hz, sigma(HV) is centred at 0.23-0.43 in ln scales (a linear scale factor of 1.26-1.54) with one standard deviation of less than 0.1. sigma(HV) is quite stable across different tectonic regions, and we present a constant, as well as earthquake magnitude- and distance-dependent sigma(HV) models.}, language = {en} } @article{GomezZapataPittoreCottonetal.2022, author = {Gomez-Zapata, Juan Camilo and Pittore, Massimiliano and Cotton, Fabrice and Lilienkamp, Henning and Shinde, Simantini and Aguirre, Paula and Santa Maria, Hernan}, title = {Epistemic uncertainty of probabilistic building exposure compositions in scenario-based earthquake loss models}, series = {Bulletin of Earthquake Engineering}, volume = {20}, journal = {Bulletin of Earthquake Engineering}, number = {5}, publisher = {Springer}, address = {Dordrecht}, issn = {1570-761X}, doi = {10.1007/s10518-021-01312-9}, pages = {2401 -- 2438}, year = {2022}, abstract = {In seismic risk assessment, the sources of uncertainty associated with building exposure modelling have not received as much attention as other components related to hazard and vulnerability. Conventional practices such as assuming absolute portfolio compositions (i.e., proportions per building class) from expert-based assumptions over aggregated data crudely disregard the contribution of uncertainty of the exposure upon earthquake loss models. In this work, we introduce the concept that the degree of knowledge of a building stock can be described within a Bayesian probabilistic approach that integrates both expert-based prior distributions and data collection on individual buildings. We investigate the impact of the epistemic uncertainty in the portfolio composition on scenario-based earthquake loss models through an exposure-oriented logic tree arrangement based on synthetic building portfolios. For illustrative purposes, we consider the residential building stock of Valparaiso (Chile) subjected to seismic ground-shaking from one subduction earthquake. We have found that building class reconnaissance, either from prior assumptions by desktop studies with aggregated data (top-down approach), or from building-by-building data collection (bottom-up approach), plays a fundamental role in the statistical modelling of exposure. To model the vulnerability of such a heterogeneous building stock, we require that their associated set of structural fragility functions handle multiple spectral periods. Thereby, we also discuss the relevance and specific uncertainty upon generating either uncorrelated or spatially cross-correlated ground motion fields within this framework. We successively show how various epistemic uncertainties embedded within these probabilistic exposure models are differently propagated throughout the computed direct financial losses. This work calls for further efforts to redesign desktop exposure studies, while also highlighting the importance of exposure data collection with standardized and iterative approaches.}, language = {en} } @article{DouglasAkkarAmerietal.2014, author = {Douglas, John and Akkar, Sinan and Ameri, Gabriele and Bard, Pierre-Yves and Bindi, Dino and Bommer, Julian J. and Bora, Sanjay Singh and Cotton, Fabrice and Derras, Boumediene and Hermkes, Marcel and Kuehn, Nicolas Martin and Luzi, Lucia and Massa, Marco and Pacor, Francesca and Riggelsen, Carsten and Sandikkaya, M. Abdullah and Scherbaum, Frank and Stafford, Peter J. and Traversa, Paola}, title = {Comparisons among the five ground-motion models developed using RESORCE for the prediction of response spectral accelerations due to earthquakes in Europe and the Middle East}, series = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering}, volume = {12}, journal = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering}, number = {1}, publisher = {Springer}, address = {Dordrecht}, issn = {1570-761X}, doi = {10.1007/s10518-013-9522-8}, pages = {341 -- 358}, year = {2014}, abstract = {This article presents comparisons among the five ground-motion models described in other articles within this special issue, in terms of data selection criteria, characteristics of the models and predicted peak ground and response spectral accelerations. Comparisons are also made with predictions from the Next Generation Attenuation (NGA) models to which the models presented here have similarities (e.g. a common master database has been used) but also differences (e.g. some models in this issue are nonparametric). As a result of the differing data selection criteria and derivation techniques the predicted median ground motions show considerable differences (up to a factor of two for certain scenarios), particularly for magnitudes and distances close to or beyond the range of the available observations. The predicted influence of style-of-faulting shows much variation among models whereas site amplification factors are more similar, with peak amplification at around 1s. These differences are greater than those among predictions from the NGA models. The models for aleatory variability (sigma), however, are similar and suggest that ground-motion variability from this region is slightly higher than that predicted by the NGA models, based primarily on data from California and Taiwan.}, language = {en} } @article{MolkenthinScherbaumGriewanketal.2017, author = {Molkenthin, Christian and Scherbaum, Frank and Griewank, Andreas and Leovey, Hernan and Kucherenko, Sergei and Cotton, Fabrice}, title = {Derivative-Based Global Sensitivity Analysis: Upper Bounding of Sensitivities in Seismic-Hazard Assessment Using Automatic Differentiation}, series = {Bulletin of the Seismological Society of America}, volume = {107}, journal = {Bulletin of the Seismological Society of America}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0037-1106}, doi = {10.1785/0120160185}, pages = {984 -- 1004}, year = {2017}, abstract = {Seismic-hazard assessment is of great importance within the field of engineering seismology. Nowadays, it is common practice to define future seismic demands using probabilistic seismic-hazard analysis (PSHA). Often it is neither obvious nor transparent how PSHA responds to changes in its inputs. In addition, PSHA relies on many uncertain inputs. Sensitivity analysis (SA) is concerned with the assessment and quantification of how changes in the model inputs affect the model response and how input uncertainties influence the distribution of the model response. Sensitivity studies are challenging primarily for computational reasons; hence, the development of efficient methods is of major importance. Powerful local (deterministic) methods widely used in other fields can make SA feasible, even for complex models with a large number of inputs; for example, automatic/algorithmic differentiation (AD)-based adjoint methods. Recently developed derivative-based global sensitivity measures can combine the advantages of such local SA methods with efficient sampling strategies facilitating quantitative global sensitivity analysis (GSA) for complex models. In our study, we propose and implement exactly this combination. It allows an upper bounding of the sensitivities involved in PSHA globally and, therefore, an identification of the noninfluential and the most important uncertain inputs. To the best of our knowledge, it is the first time that derivative-based GSA measures are combined with AD in practice. In addition, we show that first-order uncertainty propagation using the delta method can give satisfactory approximations of global sensitivity measures and allow a rough characterization of the model output distribution in the case of PSHA. An illustrative example is shown for the suggested derivative-based GSA of a PSHA that uses stochastic ground-motion simulations.}, language = {en} } @article{ZhuCottonPilz2019, author = {Zhu, Chuanbin and Cotton, Fabrice and Pilz, Marco}, title = {Testing the Depths to 1.0 and 2.5 km/s Velocity Isosurfaces in a Velocity Model for Japan and Implications for Ground-Motion Modeling}, series = {Bulletin of the Seismological Society of America}, volume = {109}, journal = {Bulletin of the Seismological Society of America}, number = {6}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0037-1106}, doi = {10.1785/0120190016}, pages = {2710 -- 2721}, year = {2019}, abstract = {In the Next Generation Attenuation West2 (NGA-West2) project, a 3D subsurface structure model (Japan Seismic Hazard Information Station [J-SHIS]) was queried to establish depths to 1.0 and 2.5 km/s velocity isosurfaces for sites without depth measurement in Japan. In this article, we evaluate the depth parameters in the J-SHIS velocity model by comparing them with their corresponding site-specific depth measurements derived from selected KiK-net velocity profiles. The comparison indicates that the J-SHIS model underestimates site depths at shallow sites and overestimates depths at deep sites. Similar issues were also identified in the southern California basin model. Our results also show that these underestimations and over-estimations have a potentially significant impact on ground-motion prediction using NGA-West2 ground-motion models (GMMs). Site resonant period may be considered as an alternative to depth parameter in the site term of a GMM.}, language = {en} } @article{ZaccarelliBindiStrolloetal.2019, author = {Zaccarelli, Riccardo and Bindi, Dino and Strollo, Angelo and Quinteros, Javier and Cotton, Fabrice}, title = {Stream2segment: An Open-Source Tool for Downloading, Processing, and Visualizing Massive Event-Based Seismic Waveform Datasets}, series = {Seismological research letters}, volume = {90}, journal = {Seismological research letters}, number = {5}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0895-0695}, doi = {10.1785/0220180314}, pages = {2028 -- 2038}, year = {2019}, abstract = {The task of downloading comprehensive datasets of event-based seismic waveforms has been made easier through the development of standardized webservices but is still highly nontrivial because the likelihood of temporary network failures or subtle data errors naturally increases when the amount of requested data is in the order of millions of relatively short segments. This is even more challenging because the typical workflow is not restricted to a single massive download but consists of fetching all possible available input data (e.g., with several repeated download executions) for a processing stage producing any desired user-defined output. Here, we present stream2segment, a highly customizable Python 2+3 package helping the user in the entire workflow of downloading, inspecting, and processing event-based seismic data by means of a relational database management system as archiving storage, which has clear performance and usability advantages, and an integrated processing subroutine requiring a configuration file and a single Python function to produce user-defined output. Stream2segment can also produce diagnostic maps or user-defined plots, which, unlike existing tools, do not require external software dependencies and are not static images but instead are interactive browser-based applications ideally suited for data inspection or annotation tasks and subsequent training of classifiers in foreseen supervised machine-learning applications. Stream2segment has already been used as a data quality tool for datasets within the European Integrated Data Archive and to create a weak-motion database (in the form of a so-called flat file) for the stable continental region of Europe in the context of the European Ground Shaking Intensity Model service, in turn an important building block for seismic hazard studies.}, language = {en} } @article{LilienkampvonSpechtWeatherilletal.2022, author = {Lilienkamp, Henning and von Specht, Sebastian and Weatherill, Graeme and Caire, Giuseppe and Cotton, Fabrice}, title = {Ground-Motion modeling as an image processing task}, series = {Bulletin of the Seismological Society of America}, volume = {112}, journal = {Bulletin of the Seismological Society of America}, number = {3}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0037-1106}, doi = {10.1785/0120220008}, pages = {1565 -- 1582}, year = {2022}, abstract = {We construct and examine the prototype of a deep learning-based ground-motion model (GMM) that is both fully data driven and nonergodic. We formulate ground-motion modeling as an image processing task, in which a specific type of neural network, the U-Net, relates continuous, horizontal maps of earthquake predictive parameters to sparse observations of a ground-motion intensity measure (IM). The processing of map-shaped data allows the natural incorporation of absolute earthquake source and observation site coordinates, and is, therefore, well suited to include site-, source-, and path-specific amplification effects in a nonergodic GMM. Data-driven interpolation of the IM between observation points is an inherent feature of the U-Net and requires no a priori assumptions. We evaluate our model using both a synthetic dataset and a subset of observations from the KiK-net strong motion network in the Kanto basin in Japan. We find that the U-Net model is capable of learning the magnitude???distance scaling, as well as site-, source-, and path-specific amplification effects from a strong motion dataset. The interpolation scheme is evaluated using a fivefold cross validation and is found to provide on average unbiased predictions. The magnitude???distance scaling as well as the site amplification of response spectral acceleration at a period of 1 s obtained for the Kanto basin are comparable to previous regional studies.}, language = {en} } @article{PilzCotton2019, author = {Pilz, Marco and Cotton, Fabrice}, title = {Does the One-Dimensional Assumption Hold for Site Response Analysis?}, 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 = {2}, publisher = {Earthquake Engineering Research Institute}, address = {Oakland}, issn = {8755-2930}, doi = {10.1193/050718EQS113M}, pages = {883 -- 905}, year = {2019}, abstract = {The one-dimensional (1-D) approach is still the dominant method to incorporate site effects in engineering applications. To bridge the 1-D to multidimensional site response analysis, we develop quantitative criteria and a reproducible method to identify KiK-net sites with significant deviations from 1-D behavior. We found that 158 out of 354 show two-dimensional (2-D) and three-dimensional (3-D) effects, extending the resonance toward shorter periods at which 2-D or 3-D site effects exceed those of the classic 1-D configurations and imposing an additional amplification to that caused by the impedance contrast alone. Such 2-D and 3-D effects go along with a large within-station ground motion variability. Remarkably, these effects are found to be more pronounced for small impedance contrasts. While it is hardly possible to identify common features in ground motion behavior for stations with similar topography typologies, it is not over-conservative to apply a safety factor to account for 2-D and 3-D site effects in ground motion modeling.}, language = {en} }