@article{KtenidouRoumeliotiAbrahamsonetal.2018,
  author    = {Ktenidou, Olga-Joan and Roumelioti, Zafeiria and Abrahamson, Norman and Cotton, Fabrice Pierre 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{GruenthalStromeyerBosseetal.2018,
  author    = {Gr{\"u}nthal, Gottfried and Stromeyer, Dietrich and Bosse, Christian and Cotton, Fabrice Pierre and Bindi, Dino},
  title     = {The probabilistic seismic hazard assessment of Germany-version 2016, considering the range of epistemic uncertainties and aleatory variability},
  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    = {10},
  publisher = {Springer},
  address   = {Dordrecht},
  issn      = {1570-761X},
  doi       = {10.1007/s10518-018-0315-y},
  pages     = {4339 -- 4395},
  year      = {2018},
  abstract  = {The basic seismic load parameters for the upcoming national design regulation for DIN EN 1998-1/NA result from the reassessment of the seismic hazard supported by the German Institution for Civil Engineering (DIBt). This 2016 version of the national seismic hazard assessment for Germany is based on a comprehensive involvement of all accessible uncertainties in models and parameters and includes the provision of a rational framework for integrating ranges of epistemic uncertainties and aleatory variabilities in a comprehensive and transparent way. The developed seismic hazard model incorporates significant improvements over previous versions. It is based on updated and extended databases, it includes robust methods to evolve sets of models representing epistemic uncertainties, and a selection of the latest generation of ground motion prediction equations. The new earthquake model is presented here, which consists of a logic tree with 4040 end branches and essential innovations employed for a realistic approach. The output specifications were designed according to the user oriented needs as suggested by two review teams supervising the entire project. Seismic load parameters, for rock conditions of nu(S30) = 800 m/s, are calculated for three hazard levels (10, 5 and 2\% probability of occurrence or exceedance within 50 years) and delivered in the form of uniform hazard spectra, within the spectral period range 0.02-3 s, and seismic hazard maps for peak ground acceleration, spectral response accelerations and for macroseismic intensities. Results are supplied as the mean, the median and the 84th percentile. A broad analysis of resulting uncertainties of calculated seismic load parameters is included. The stability of the hazard maps with respect to previous versions and the cross-border comparison is emphasized.},
  language  = {en}
}
@misc{GruenthalStromeyerBosseetal.2018,
  author    = {Gr{\"u}nthal, Gottfried and Stromeyer, Dietrich and Bosse, Christian and Cotton, Fabrice Pierre and Bindi, Dino},
  title     = {Correction to: The probabilistic seismic hazard assessment of Germanyversion 2016, considering the range of epistemic uncertainties and aleatory variability (vol 16, pg 4339, 2018)},
  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    = {10},
  publisher = {Springer},
  address   = {Dordrecht},
  issn      = {1570-761X},
  doi       = {10.1007/s10518-018-0398-5},
  pages     = {4397 -- 4398},
  year      = {2018},
  abstract  = {One paragraph of the manuscript of the paper has been inadvertently omitted in the very final stage of its compilation due to a technical mistake. Since this paragraph discusses the declustering of the used earthquake catalogue and is therefore necessary for the understanding of the seismicity data preprocessing, the authors decided to provide this paragraph in form of a correction. The respective paragraph belongs to chapter 2 of the paper, where it was placed originally, and should be inserted into the published paper before the second to the last paragraph. The omitted text reads as follows:},
  language  = {en}
}
@misc{JaraSanchezReyesSocquetetal.2018,
  author    = {Jara, Jorge and S{\´a}nchez-Reyes, Hugo and Socquet, Anne and Cotton, Fabrice Pierre and Virieux, Jean and Maksymowicz, Andrei and D{\´i}az-Mojica, John and Walpersdorf, Andrea and Ruiz, Javier and Cotte, Nathalie and Norabuena, Edmundo},
  title     = {Corrigendum to: Kinematic study of Iquique 2014 Mw 8.1 earthquake: Understanding the segmentation of the seismogenic zone. - (Earth and planetary science letters. - 503 (2018) S. 131 - 143)},
  series = {Earth and planetary science letters},
  volume    = {506},
  journal   = {Earth and planetary science letters},
  publisher = {Elsevier},
  address   = {Amsterdam},
  issn      = {0012-821X},
  doi       = {10.1016/j.epsl.2018.11.026},
  pages     = {347 -- 347},
  year      = {2018},
  abstract  = {We study the rupture processes of Iquique earthquake 8.1 (2014/04/01) and its largest aftershock 7.7 (2014/04/03) that ruptured the North Chile subduction zone. High-rate Global Positioning System (GPS) recordings and strong motion data are used to reconstruct the evolution of the slip amplitude, rise time and rupture time of both earthquakes. A two-step inversion scheme is assumed, by first building prior models for both earthquakes from the inversion of the estimated static displacements and then, kinematic inversions in the frequency domain are carried out taken into account this prior information. The preferred model for the mainshock exhibits a seismic moment of 1.73 × 1021 Nm ( 8.1) and maximum slip of ∼9 m, while the aftershock model has a seismic moment of 3.88 × 1020 ( 7.7) and a maximum slip of ∼3 m. For both earthquakes, the final slip distributions show two asperities (a shallow one and a deep one) separated by an area with significant slip deficit. This suggests a segmentation along-dip which might be related to a change of the dipping angle of the subducting slab inferred from gravimetric data. Along-strike, the areas where the seismic ruptures stopped seem to be well correlated with geological features observed from geophysical information (high-resolution bathymetry, gravimetry and coupling maps) that are representative of the long-term segmentation of the subduction margin. Considering the spatially limited portions that were broken by these two earthquakes, our results support the idea that the seismic gap is not filled yet.},
  language  = {en}
}
@article{PinaValdesSocquetCotton2018,
  author    = {Pina-Valdes, Jesus and Socquet, Anne and Cotton, Fabrice Pierre},
  title     = {Insights on the Japanese Subduction Megathrust Properties From Depth and Lateral Variability of Observed Ground Motions},
  series = {Journal of geophysical research : Solid earth},
  volume    = {123},
  journal   = {Journal of geophysical research : Solid earth},
  number    = {10},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {2169-9313},
  doi       = {10.1029/2018JB015743},
  pages     = {8937 -- 8956},
  year      = {2018},
  abstract  = {Two ground motion prediction equation models for subduction zones have been tested using a public ground motion database of the KiK-net records obtained by automated processing protocols (Dawood et al., 2016, https://doi.org/10.1193/071214EQS106). The database contains records of more than 700 interface earthquakes that occurred on the Japan subduction between 1998 and 2012. The Zhao et al. (2006, https://doi.org/10.1785/0120050122) ground motion prediction equation was shown to be the best suited model for the region. It was then used as backbone to analyze the variability of ground motion records. The residuals between observed and predicted ground motions have been analyzed to study the spatial variation of the earthquakes' ground motion frequency content on the Japan megathrust. This analysis revealed a depth dependency of generated ground motions consistent with the downdip segmentation proposed for subduction interfaces (Lay et al., 2012, https://doi.org/10.1029/2011JB009133), a regional ground motion dependency that may be related with lateral variations of the mechanical properties of the subduction interface and a high-frequency radiations drop in the earthquake sequence that preceded the Tohoku-Oki earthquake Mw 9.0. The regional ground motion dependency suggests the existence of different domains along trench of the Japan subduction megathrust that control the ground motions and the wave radiation patterns of interface earthquakes. The location of their boundaries is consistent with the extension of the rupture of the 2011 Tohoku-Oki earthquake, with pre-Tohoku interseismic coupling, and with the free air gravity anomalies.},
  language  = {en}
}
@article{JaraSanchezReyesSocquetetal.2018,
  author    = {Jara, Jorge and Sanchez-Reyes, Hugo and Socquet, Anne and Cotton, Fabrice Pierre and Virieux, Jean and Maksymowicz, Andrei and Diaz-Mojica, John and Walpersdorf, Andrea and Ruiz, Javier and Cotte, Nathalie and Norabuena, Edmundo},
  title     = {Kinematic study of Iquique 2014 M-w 8.1 earthquake},
  series = {Earth \& planetary science letters},
  volume    = {503},
  journal   = {Earth \& planetary science letters},
  publisher = {Elsevier},
  address   = {Amsterdam},
  issn      = {0012-821X},
  doi       = {10.1016/j.epsl.2018.09.025},
  pages     = {131 -- 143},
  year      = {2018},
  abstract  = {We study the rupture processes of Iquique earthquake M-w 8.1 (2014/04/01) and its largest aftershock M-w 7.7 (2014/04/03) that ruptured the North Chile subduction zone. High-rate Global Positioning System (GPS) recordings and strong motion data are used to reconstruct the evolution of the slip amplitude, rise time and rupture time of both earthquakes. A two-step inversion scheme is assumed, by first building prior models for both earthquakes from the inversion of the estimated static displacements and then, kinematic inversions in the frequency domain are carried out taken into account this prior information. The preferred model for the mainshock exhibits a seismic moment of 1.73 x 10(21) Nm (M-w 8.1) and maximum slip of similar to 9 m, while the aftershock model has a seismic moment of 3.88 x 10(20) (M-w 7.7) and a maximum slip of similar to 3 m. For both earthquakes, the final slip distributions show two asperities (a shallow one and a deep one) separated by an area with significant slip deficit. This suggests a segmentation along-dip which might be related to a change of the dipping angle of the subducting slab inferred from gravimetric data. Along-strike, the areas where the seismic ruptures stopped seem to be well correlated with geological features observed from geophysical information (high-resolution bathymetry, gravimetry and coupling maps) that are representative of the long-term segmentation of the subduction margin. Considering the spatially limited portions that were broken by these two earthquakes, our results support the idea that the seismic gap is not filled yet. (C) 2018 Elsevier B.V. All rights reserved.},
  language  = {en}
}
@article{DerrasBardCotton2017,
  author    = {Derras, Boumediene and Bard, Pierre-Yves and Cotton, Fabrice Pierre},
  title     = {V-S30, slope, H-800 and f(0): performance of various site-condition proxies in reducing ground-motion aleatory variability and predicting nonlinear site response},
  series = {Earth, planets and space},
  volume    = {69},
  journal   = {Earth, planets and space},
  publisher = {Springer},
  address   = {Heidelberg},
  issn      = {1880-5981},
  doi       = {10.1186/s40623-017-0718-z},
  pages     = {1623 -- 1629},
  year      = {2017},
  abstract  = {The aim of this paper is to investigate the ability of various site-condition proxies (SCPs) to reduce ground-motion aleatory variability and evaluate how SCPs capture nonlinearity site effects. The SCPs used here are time-averaged shear-wave velocity in the top 30 m (V-S30), the topographical slope (slope), the fundamental resonance frequency (f(0)) and the depth beyond which V-s exceeds 800 m/s (H800). We considered first the performance of each SCP taken alone and then the combined performance of the 6 SCP pairs [V-S30-f(0)], [V-S30-H-800], [f(0)-slope], [H-800-slope], [V-S30-slope] and [f(0)-H-800]. This analysis is performed using a neural network approach including a random effect applied on a KiK-net subset for derivation of ground-motion prediction equations setting the relationship between various ground-motion parameters such as peak ground acceleration, peak ground velocity and pseudo-spectral acceleration PSA (T), and Mw, RJB, focal depth and SCPs. While the choice of SCP is found to have almost no impact on the median groundmotion prediction, it does impact the level of aleatory uncertainty. VS30 is found to perform the best of single proxies at short periods (T < 0.6 s), while f(0) and H-800 perform better at longer periods; considering SCP pairs leads to significant improvements, with particular emphasis on [V-S30-H-800] and [f(0)-slope] pairs. The results also indicate significant nonlinearity on the site terms for soft sites and that the most relevant loading parameter for characterising nonlinear site response is the "stiff" spectral ordinate at the considered period.},
  language  = {en}
}
@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{ZaccarelliBindiStrolloetal.2019,
  author    = {Zaccarelli, Riccardo and Bindi, Dino and Strollo, Angelo and Quinteros, Javier and Cotton, Fabrice Pierre},
  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{LetortRetailleauBoueetal.2018,
  author    = {Letort, Jean and Retailleau, Lise and Boue, Pierre and Radiguet, Mathilde and Gardonio, Blandine and Cotton, Fabrice Pierre 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}
}