@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{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{DiGiacomoBindiParolaietal.2011,
  author    = {Di Giacomo, Domenico and Bindi, Dino and Parolai, Stefano and Oth, Adrien},
  title     = {Residual analysis of teleseismic P-wave energy magnitude estimates: inter- and intrastation variability},
  series = {Geophysical journal international},
  volume    = {185},
  journal   = {Geophysical journal international},
  number    = {3},
  publisher = {Wiley-Blackwell},
  address   = {Malden},
  issn      = {0956-540X},
  doi       = {10.1111/j.1365-246X.2011.05019.x},
  pages     = {1444 -- 1454},
  year      = {2011},
  abstract  = {P>Computing the magnitude of an earthquake requires correcting for the propagation effects from the source to the receivers. This is often accomplished by performing numerical simulations using a suitable Earth model. In this work, the energy magnitude M(e) is considered and its determination is performed using theoretical spectral amplitude decay functions over teleseismic distances based on the global Earth model AK135Q. Since the high frequency part (above the corner frequency) of the source spectrum has to be considered in computing M(e), the influence of propagation and site effects may not be negligible and they could bias the single station M(e) estimations. Therefore, in this study we assess the inter- and intrastation distributions of errors by considering the M(e) residuals computed for a large data set of earthquakes recorded at teleseismic distances by seismic stations deployed worldwide. To separate the inter- and intrastation contribution of errors, we apply a maximum likelihood approach to the M(e) residuals. We show that the interstation errors (describing a sort of site effect for a station) are within +/- 0.2 magnitude units for most stations and their spatial distribution reflects the expected lateral variation affecting the velocity and attenuation of the Earth's structure in the uppermost layers, not accounted for by the 1-D AK135Q model. The variance of the intrastation error distribution (describing the record-to-record component of variability) is larger than the interstation one (0.240 against 0.159), and the spatial distribution of the errors is not random but shows specific patterns depending on the source-to-station paths. The set of coefficients empirically determined may be used in the future to account for the heterogeneities of the real Earth not considered in the theoretical calculations of the spectral amplitude decay functions used to correct the recorded data for propagation effects.},
  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{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{BindiSpallarossaPicozzietal.2018,
  author    = {Bindi, Dino and Spallarossa, D. and Picozzi, M. and Scafidi, D. and Cotton, Fabrice},
  title     = {Impact of magnitude selection on aleatory variability associated with ground-motion prediction equations},
  series = {Bulletin of the Seismological Society of America},
  volume    = {108},
  journal   = {Bulletin of the Seismological Society of America},
  number    = {3A},
  publisher = {Seismological Society of America},
  address   = {Albany},
  issn      = {0037-1106},
  doi       = {10.1785/0120170356},
  pages     = {1427 -- 1442},
  year      = {2018},
  abstract  = {In this study, we analyzed 10 yrs of seismicity in central Italy from 2008 to 2017, a period witnessing more than 1400 earthquakes in the magnitude range 2.5≤Mw≤6.5⁠. The data set includes the main sequences that have occurred in the area, including those associated with the 2009 Mw 6.3 L'Aquila earthquake and the 2016-2017 sequence (⁠Mw 6.2 Amatrice, Mw 6.1 Visso, and Mw 6.5 Norcia earthquakes). We calibrated a local magnitude scale, investigating the impact of changing the reference distance at which the nonparametric attenuation is tied to the zero-magnitude attenuation function for southern California. We also developed an attenuation model to compute the radiated seismic energy (⁠Es⁠) from the time integral of the squared ground-motion velocity. Seismic moment (⁠M0⁠) and stress drop (⁠Δσ⁠) were estimated for each earthquake by fitting a ω-square model to the source spectra obtained by applying a nonparametric spectral inversion. The Δσ-values vary over three orders of magnitude from about 0.1 to 10 MPa, the larger values associated with the mainshocks. The Δσ-values describe a lognormal distribution with mean and standard deviation equal to log(Δσ)=(-0.25±0.45) (i.e., the mean Δσ is 0.57 MPa, with a 95\% confidence interval from 0.08 to 4.79 MPa). The Δσ variability introduces a spread in the distribution of seismic energy versus moment, with differences in energy up two orders of magnitudes for earthquakes with the same moment. The variability in the high-frequency spectral levels is captured by the local magnitude (⁠ML⁠), which scales with radiated energy as ML=(-1.59+0.52logEs) for logEs≤10.26 and ML=(-1.38+0.50logEs) otherwise. As the peak ground velocity increases with increasing Δσ⁠, local and energy magnitudes perform better than moment magnitude as predictors for the shaking potential. The availability of different magnitude scales and source parameters for a large earthquake population will help characterize the between-event ground-motion variability in central Italy.},
  language  = {en}
}
@article{BindiKothaWeatherilletal.2018,
  author    = {Bindi, Dino and Kotha, Sreeram Reddy and Weatherill, Graeme and Lanzano, Giovanni and Luzi, Lucia and Cotton, Fabrice},
  title     = {The pan-European engineering strong motion (ESM) flatfile},
  series = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering},
  volume    = {17},
  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-018-0466-x},
  pages     = {583 -- 602},
  year      = {2018},
  abstract  = {We present the results of a consistency check performed over the flatfile extracted from the engineering strong motion (ESM) database. The flatfile includes 23,014 recordings from 2179 earthquakes in the magnitude range from 3.5 to 7.8 that occurred since the 1970s in Europe and Middle East, as presented in the companion article by Lanzano et al. (Bull Earthq Eng, 2018a). The consistency check is developed by analyzing different residual distributions obtained from ad-hoc ground motion prediction equations for the absolute spectral acceleration (SA), displacement and Fourier amplitude spectra (FAS). Only recordings from earthquakes shallower than 40 km are considered in the analysis. The between-event, between-station and event-and-station corrected residuals are computed by applying a mixed-effect regression. We identified those earthquakes, stations, and recordings showing the largest deviations from the GMPE median predictions, and also evaluated the statistical uncertainty on the median model to get insights on the applicable magnitude-distance ranges and the usable period (or frequency) range. We observed that robust median predictions are obtained up to 8 s for SA and up to 20 Hz for FAS, although median predictions for Mw ≥ 7 show significantly larger uncertainties with 'bumps' starting above 5 s for SA and below 0.3 Hz for FAS. The between-station variance dominates over the other residual variances, and the dependence of the between-station residuals on logarithm of Vs30 is well-described by a piece-wise linear function with period-dependent slopes and hinge velocity around 580 m/s. Finally, we compared the between-event residuals obtained by considering two different sources of moment magnitude. The results show that, at long periods, the between-event terms from the two regressions have a weak correlation and the overall between-event variability is dissimilar, highlighting the importance of magnitude source in the regression results.},
  language  = {en}
}
@article{ZoellerUllahBindietal.2017,
  author    = {Z{\"o}ller, Gert and Ullah, Shahid and Bindi, Dino and Parolai, Stefano and Mikhailova, Natalya},
  title     = {The largest expected earthquake magnitudes in Central Asia},
  series = {Seismicity, fault rupture and earthquake hazards in slowly deforming regions},
  volume    = {432},
  journal   = {Seismicity, fault rupture and earthquake hazards in slowly deforming regions},
  publisher = {The Geological Society},
  address   = {London},
  isbn      = {978-1-86239-745-3},
  issn      = {0305-8719},
  doi       = {10.1144/SP432.3},
  pages     = {29 -- 40},
  year      = {2017},
  abstract  = {The knowledge of the largest expected earthquake magnitude in a region is one of the key issues in probabilistic seismic hazard calculations and the estimation of worst-case scenarios. Earthquake catalogues are the most informative source of information for the inference of earthquake magnitudes. We analysed the earthquake catalogue for Central Asia with respect to the largest expected magnitudes m(T) in a pre-defined time horizon T-f using a recently developed statistical methodology, extended by the explicit probabilistic consideration of magnitude errors. For this aim, we assumed broad error distributions for historical events, whereas the magnitudes of recently recorded instrumental earthquakes had smaller errors. The results indicate high probabilities for the occurrence of large events (M >= 8), even in short time intervals of a few decades. The expected magnitudes relative to the assumed maximum possible magnitude are generally higher for intermediate-depth earthquakes (51-300 km) than for shallow events (0-50 km). For long future time horizons, for example, a few hundred years, earthquakes with M >= 8.5 have to be taken into account, although, apart from the 1889 Chilik earthquake, it is probable that no such event occurred during the observation period of the catalogue.},
  language  = {en}
}
@article{GruenthalStromeyerBosseetal.2018,
  author    = {Gr{\"u}nthal, Gottfried and Stromeyer, Dietrich and Bosse, Christian and Cotton, Fabrice 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 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}
}
@article{BindiMarzoratiParolaietal.2009,
  author    = {Bindi, Dino and Marzorati, Simone and Parolai, Stefano and Strollo, Angelo and Jaeckel, Karl-Heinz},
  title     = {Empirical spectral ratios estimated in two deep sedimentary basins using microseisms recorded by short-period seismometers},
  issn      = {0956-540X},
  doi       = {10.1111/j.1365-246X.2008.03958.x},
  year      = {2009},
  abstract  = {In this work, we analyse continuous measurements of microseisms to assess the reliability of the fundamental resonance frequency estimated by means of the horizontal-to-vertical (H/V) spectral ratio within the 0.1-1 Hz frequency range, using short-period sensors (natural period of 1 s). We apply the H/V technique to recordings of stations installed in two alluvial basins with different sedimentary cover thicknesses-the Lower Rhine Embayment (Germany) and the Gubbio Plain (Central Italy). The spectral ratios are estimated over the time-frequency domain, and we discuss the reliability of the results considering both the variability of the microseism activity and the amplitude of the instrumental noise. We show that microseisms measured by short period sensors allow the retrieval of fundamental resonance frequencies greater than about 0.1-0.2 Hz, with this lower frequency bound depending on the relative amplitude of the microseism signal and the self-noise of the instruments. In particular, we show an example where the considered short-period sensor is connected to instruments characterized by an instrumental noise level which allows detecting only fundamental frequencies greater than about 0.4 Hz. Since the frequency at which the peak of the H/V spectral ratio is biased depends upon the seismic signal-to-instrument noise ratio, the power spectral amplitude of instrumental self- noise should be always considered when interpreting the frequency of the peak as the fundamental resonance frequency of the investigated site.},
  language  = {en}
}
@article{PicozziParolaiBindietal.2009,
  author    = {Picozzi, Matteo and Parolai, Stefano and Bindi, Dino and Strollo, Angelo},
  title     = {Characterization of shallow geology by high-frequency seismic noise tomography},
  issn      = {0956-540X},
  doi       = {10.1111/j.1365-246X.2008.03966.x},
  year      = {2009},
  abstract  = {To study the applicability of the passive seismic interferometry technique to near-surface geological studies, seismic noise recordings from a small scale 2-D array of seismic stations were performed in the test site of Nauen (Germany). Rayleigh wave Green's functions were estimated for different frequencies. A tomographic inversion of the traveltimes estimated for each frequency from the Green's functions is then performed, allowing the laterally varying 3-D surfacewave velocity structure below the array to be retrieved at engineering-geotechnical scales. Furthermore, a 2-D S-wave velocity cross-section is obtained by combining 1-D velocity structures derived from the inversion of the dispersion curves extracted at several points along a profile where other geophysical analyses were performed. It is shown that the cross-section from passive seismic interferometry provides a clear image of the local structural heterogeneities that are in excellent agreement with georadar and geoelectrical results. Such findings indicate that the interferometry analysis of seismic noise is potentially of great interest for deriving the shallow 3-D velocity structure in urban areas.},
  language  = {en}
}
@article{StrolloParolaiBindietal.2012,
  author    = {Strollo, Angelo and Parolai, Stefano and Bindi, Dino and Chiauzzi, Leonardo and Pagliuca, Rossella and Mucciarelli, Marco and Zschau, Jochen},
  title     = {Microzonation of Potenza (Southern Italy) in terms of spectral intensity ratio using joint analysis of earthquakes and ambient noise},
  series = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering},
  volume    = {10},
  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-011-9256-4},
  pages     = {493 -- 516},
  year      = {2012},
  abstract  = {A temporary seismic network composed of 11 stations was installed in the city of Potenza (Southern Italy) to record local and regional seismicity within the context of a national project funded by the Italian Department of Civil Protection (DPC). Some stations were moved after a certain time in order to increase the number of measurement points, leading to a total of 14 sites within the city by the end of the experiment. Recordings from 26 local earthquakes (M-l 2.2-3.8 ) were analyzed to compute the site responses at the 14 sites by applying both reference and non-reference site techniques. Furthermore, the Spectral Intensity (SI) for each local earthquake, as well as their ratios with respect to the values obtained at a reference site, were also calculated. In addition, a field survey of 233 single station noise measurements within the city was carried out to increase the information available at localities different from the 14 monitoring sites. By using the results of the correlation analysis between the horizontal-to-vertical spectral ratios computed from noise recordings (NHV) at the 14 selected sites and those derived by the single station noise measurements within the town as a proxy, the spectral intensity correction factors for site amplification obtained from earthquake analysis were extended to the entire city area. This procedure allowed us to provide a microzonation map of the urban area that can be directly used when calculating risk scenarios for civil defence purposes. The amplification factors estimated following this approach show values increasing along the main valley toward east where the detrital and alluvial complexes reach their maximum thickness.},
  language  = {en}
}
@article{PilzCottonZaccarellietal.2019,
  author    = {Pilz, Marco and Cotton, Fabrice and Zaccarelli, Riccardo and Bindi, Dino},
  title     = {Capturing Regional Variations of Hard-Rock Attenuation in Europe},
  series = {Bulletin of the Seismological Society of America},
  volume    = {109},
  journal   = {Bulletin of the Seismological Society of America},
  number    = {4},
  publisher = {Seismological Society of America},
  address   = {Albany},
  issn      = {0037-1106},
  doi       = {10.1785/0120190023},
  pages     = {1401 -- 1418},
  year      = {2019},
  abstract  = {A proper assessment of seismic reference site conditions has important applications as they represent the basis on which ground motions and amplifications are generally computed. Besides accounting for the average S-wave velocity over the uppermost 30 m (V-S30), the parameterization of high-frequency ground motions beyond source-corner frequency received significant attention. kappa, an empirical parameter introduced by Anderson and Hough (1984), is often used to represent the spectral decay of the acceleration spectrum at high frequencies. The lack of hard-rock records and the poor understanding of the physics of kappa introduced significant epistemic uncertainty in the final seismic hazard of recent projects. Thus, determining precise and accurate regional hard-rock kappa(0) values is critical. We propose an alternative procedure for capturing the reference kappa(0) on regional scales by linking thewell-known high-frequency attenuation parameter kappa and the properties of multiple-scattered coda waves. Using geological and geophysical data around more than 1300 stations for separating reference and soft soil sites and based on more than 10,000 crustal earthquake recordings, we observe that kappa(0) from multiple-scattered coda waves seems to be independent of the soil type but correlated with the hard-rock kappa(0), showing significant regional variations across Europe. The values range between 0.004 s for northern Europe and 0.020 s for the southern and southeastern parts. On the other hand, measuring kappa (and correspondingly kappa(0)) on the S-wave window (as classically proposed), the results are strongly affected by transmitted (reflected, refracted, and scattered) waves included in the analyzed window biasing the proper assessment of kappa(0). This effect is more pronounced for soft soil sites. In this way, kappa(coda)(0) can serve as a proxy for the regional hard-rock kappa(0) at the reference sites.},
  language  = {en}
}
@article{DahmHeimannFunkeetal.2018,
  author    = {Dahm, Torsten and Heimann, Sebastian and Funke, Sigward and Wendt, Siegfried and Rappsilber, Ivo and Bindi, Dino and Plenefisch, Thomas and Cotton, Fabrice},
  title     = {Seismicity in the block mountains between Halle and Leipzig, Central Germany},
  series = {Journal of seismology},
  volume    = {22},
  journal   = {Journal of seismology},
  number    = {4},
  publisher = {Springer},
  address   = {Dordrecht},
  issn      = {1383-4649},
  doi       = {10.1007/s10950-018-9746-9},
  pages     = {985 -- 1003},
  year      = {2018},
  abstract  = {On April 29, 2017 at 0:56 UTC (2:56 local time), an M (W) = 2.8 earthquake struck the metropolitan area between Leipzig and Halle, Germany, near the small town of Markranstadt. The earthquake was felt within 50 km from the epicenter and reached a local intensity of I (0) = IV. Already in 2015 and only 15 km northwest of the epicenter, a M (W) = 3.2 earthquake struck the area with a similar large felt radius and I (0) = IV. More than 1.1 million people live in the region, and the unusual occurrence of the two earthquakes led to public attention, because the tectonic activity is unclear and induced earthquakes have occurred in neighboring regions. Historical earthquakes south of Leipzig had estimated magnitudes up to M (W) ae 5 and coincide with NW-SE striking crustal basement faults. We use different seismological methods to analyze the two recent earthquakes and discuss them in the context of the known tectonic structures and historical seismicity. Novel stochastic full waveform simulation and inversion approaches are adapted for the application to weak, local earthquakes, to analyze mechanisms and ground motions and their relation to observed intensities. We find NW-SE striking normal faulting mechanisms for both earthquakes and centroid depths of 26 and 29 km. The earthquakes are located where faults with large vertical offsets of several hundred meters and Hercynian strike have developed since the Mesozoic. We use a stochastic full waveform simulation to explain the local peak ground velocities and calibrate the method to simulate intensities. Since the area is densely populated and has sensitive infrastructure, we simulate scenarios assuming that a 12-km long fault segment between the two recent earthquakes is ruptured and study the impact of rupture parameters on ground motions and expected damage.},
  language  = {en}
}
@article{BindiPicozziSpallarossaetal.2019,
  author    = {Bindi, Dino and Picozzi, Matteo and Spallarossa, Daniele and Cotton, Fabrice and Kotha, Sreeram Reddy},
  title     = {Impact of Magnitude Selection on Aleatory Variability Associated with Ground-Motion Prediction Equations},
  series = {Bulletin of the Seismological Society of America},
  volume    = {109},
  journal   = {Bulletin of the Seismological Society of America},
  number    = {1},
  publisher = {Seismological Society of America},
  address   = {Albany},
  issn      = {0037-1106},
  doi       = {10.1785/0120180239},
  pages     = {251 -- 262},
  year      = {2019},
  abstract  = {We derive a set of regional ground-motion prediction equations (GMPEs) in the Fourier amplitude spectra (FAS-GMPE) and in the spectral acceleration (SA-GMPE) domains for the purpose of interpreting the between-event residuals in terms of source parameter variability. We analyze a dataset of about 65,000 recordings generated by 1400 earthquakes (moment magnitude 2: 5 <= M-w <= 6: 5, hypocentral distance R-hypo <= 150 km) that occurred in central Italy between January 2008 and October 2017. In a companion article (Bindi, Spallarossa, et al., 2018), the nonparametric acceleration source spectra were interpreted in terms of omega-square models modified to account for deviations from a high-frequency flat plateau through a parameter named k(source). Here, the GMPEs are derived considering the moment (M-w), the local (M-L), and the energy (M-e) magnitude scales, and the between-event residuals are computed as random effects. We show that the between-event residuals for the FAS-GMPE implementing M-w are correlated with stress drop, with correlation coefficients increasing with increasing frequency up to about 10 Hz. Contrariwise, the correlation is weak for the FAS-GMPEs implementing M-L and M-e, in particular between 2 and 5 Hz, where most of the corner frequencies lie. At higher frequencies, all models show a strong correlation with k(source). The correlation with the source parameters reflects in a different behavior of the standard deviation tau of the between-event residuals with frequency. Although tau is smaller for the FAS-GMPE using M-w below 1.5 Hz, at higher frequencies, the model implementing either M-L or M-e shows smaller values, with a reduction of about 30\% at 3 Hz (i.e., from 0.3 for M-w to 0.1 for M-L). We conclude that considering magnitude scales informative for the stress-drop variability allows to reduce the between-event variability with a significant impact on the hazard assessment, in particular for studies in which the ergodic assumption on site is removed.},
  language  = {en}
}
@article{KothaCottonBindi2019,
  author    = {Kotha, Sreeram Reddy and Cotton, Fabrice and Bindi, Dino},
  title     = {Empirical models of shear-wave radiation pattern derived from large datasets of ground-shaking observations},
  series = {Scientific reports},
  volume    = {9},
  journal   = {Scientific reports},
  publisher = {Nature Publ. Group},
  address   = {London},
  issn      = {2045-2322},
  doi       = {10.1038/s41598-018-37524-4},
  pages     = {11},
  year      = {2019},
  abstract  = {Shear-waves are the most energetic body-waves radiated from an earthquake, and are responsible for the destruction of engineered structures. In both short-term emergency response and long-term risk forecasting of disaster-resilient built environment, it is critical to predict spatially accurate distribution of shear-wave amplitudes. Although decades' old theory proposes a deterministic, highly anisotropic, four-lobed shear-wave radiation pattern, from lack of convincing evidence, most empirical ground-shaking prediction models settled for an oversimplified stochastic radiation pattern that is isotropic on average. Today, using the large datasets of uniformly processed seismograms from several strike, normal, reverse, and oblique-slip earthquakes across the globe, compiled specifically for engineering applications, we could reveal, quantify, and calibrate the frequency-, distance-, and style-of-faulting dependent transition of shear-wave radiation between a stochastic-isotropic and a deterministic-anisotropic phenomenon. Consequent recalibration of empirical ground-shaking models dramatically improved their predictions: with isodistant anisotropic variations of ±40\%, and 8\% reduction in uncertainty. The outcomes presented here can potentially trigger a reappraisal of several practical issues in engineering seismology, particularly in seismic ground-shaking studies and seismic hazard and risk assessment.},
  language  = {en}
}
@article{BindiCottonSpallarossaetal.2018,
  author    = {Bindi, Dino and Cotton, Fabrice and Spallarossa, Daniele and Picozzi, Matteo and Rivalta, Eleonora},
  title     = {Temporal variability of ground shaking and stress drop in Central Italy},
  series = {Bulletin of the Seismological Society of America},
  volume    = {108},
  journal   = {Bulletin of the Seismological Society of America},
  number    = {4},
  publisher = {Seismological Society of America},
  address   = {Albany},
  issn      = {0037-1106},
  doi       = {10.1785/0120180078},
  pages     = {1853 -- 1863},
  year      = {2018},
  abstract  = {Ground-motion prediction equations (GMPEs) are calibrated to predict the intensity of ground shaking at any given location, based on earthquake magnitude, source-to-site distance, local soil amplifications, and other parameters. GMPEs are generally assumed to be independent of time; however, evidence is increasing that large earthquakes modify the shallow soil conditions and those of the fault zone for months or years. These changes may affect the intensity of shaking and result in time-dependent effects that can potentially be resolved by analyzing between-event residuals (residuals between observed and predicted ground motion for individual earthquakes averaged over all stations). Here, we analyze a data set of about 65,000 recordings for about 1400 earthquakes in the moment magnitude range 2.5-6.5 that occurred in central Italy from 2008 to 2017 to capture the temporal variability of the ground shaking at high frequency. We first compute between-event residuals for each earthquake in the Fourier domain with respect to a GMPE developed ad hoc for the analyzed data set. The between-events show large changes after the occurrence of mainshocks such as the 2009 Mw 6.3 L'Aquila, the 2016 Mw 6.2 Amatrice, and Mw 6.5 Norcia earthquakes. Within the time span of a few months after the mainshocks, the between-event contribution to the ground shaking varies by a factor 7. In particular, we find a large drop in the between-events in the aftermath of the L'Aquila earthquake, followed by a slow positive trend that leads to a recovery interrupted by a new drop at the beginning of 2014. We also quantify the frequency-dependent correlation between the Brune stress drop Δσ and the between-events. We find that the temporal changes of Δσ resemble those of the between-event residuals; in particular, during the period when the between-events show the positive trend, the average logarithm of Δσ increases with an annual rate of 0.19 (i.e., the amplification factor for Δσ is 1.56 per year). Breakpoint analysis located a change in the linear trend coefficients of Δσ versus time in February 2014, although no large earthquakes occurred at that time. Finally, the temporal variability of Δσ mirrors the relative seismic-velocity variations observed in previous studies for the same area and period, suggesting that both crack healing along the main fault system and healing of microcracks distributed at shallow depths throughout the surrounding region might be necessary to explain the wider observations of postearthquake recovery.},
  language  = {en}
}
@misc{DahmHeimannFunkeetal.2018,
  author    = {Dahm, Torsten and Heimann, Sebastian and Funke, Sigward and Wendt, Siegfried and Rappsilber, Ivo and Bindi, Dino and Plenefisch, Thomas and Cotton, Fabrice},
  title     = {Correction to: Seismicity in the block mountains between Halle and Leipzig, Central Germany: centroid moment tensors, ground motion simulation, and felt intensities of two M approximate to 3 earthquakes in 2015 and 2017 (vol 22, pg 985, 2018)},
  series = {Journal of seismology},
  volume    = {22},
  journal   = {Journal of seismology},
  number    = {6},
  publisher = {Springer},
  address   = {Dordrecht},
  issn      = {1383-4649},
  doi       = {10.1007/s10950-018-9773-6},
  pages     = {1669 -- 1671},
  year      = {2018},
  language  = {en}
}
@article{KothaBindiCotton2017,
  author    = {Kotha, Sreeram Reddy and Bindi, Dino and Cotton, Fabrice},
  title     = {From Ergodic to Region- and Site-Specific Probabilistic Seismic Hazard Assessment: Method Development and Application at European and Middle Eastern Sites},
  series = {Earthquake spectra : the professional journal of the Earthquake Engineering Research Institute},
  volume    = {33},
  journal   = {Earthquake spectra : the professional journal of the Earthquake Engineering Research Institute},
  publisher = {Earthquake Engineering Research Institute},
  address   = {Oakland},
  issn      = {8755-2930},
  doi       = {10.1193/081016EQS130M},
  pages     = {1433 -- 1453},
  year      = {2017},
  abstract  = {The increasing numbers of recordings at individual sites allows quantification of empirical linear site-response adjustment factors (delta S2S(s)) from the ground motion prediction equation (GMPE) residuals. The delta S2S(s) are then used to linearly scale the ergodic GMPE predictions to obtain site-specific ground motion predictions in a partially non-ergodic Probabilistic Seismic Hazard Assessment (PSHA). To address key statistical and conceptual issues in the current practice, we introduce a novel empirical region-and site-specific PSHA methodology wherein, (1) site-to-site variability (phi(S2S)) is first estimated as a random-variance in a mixed-effects GMPE regression, (2) delta S2S(s) at new sites with strong motion are estimated using the a priori phi(S2S), and (3) the GMPE site-specific single-site aleatory variability sigma(ss,s) is replaced with a generic site-corrected aleatory variability sigma(0). Comparison of region- and site-specific hazard curves from our method against the traditional ergodic estimates at 225 sites in Europe and Middle East shows an approximate 50\% difference in predicted ground motions over a range of hazard levels-a strong motivation to increase seismological monitoring of critical facilities and enrich regional ground motion data sets.},
  language  = {en}
}