@article{WeatherillKothaCotton2020, author = {Weatherill, Graeme and Kotha, Sreeram Reddy and Cotton, Fabrice Pierre}, title = {Re-thinking site amplification in regional seismic risk assessment}, series = {Earthquake spectra : the professional journal of the Earthquake Engineering Research Institute}, volume = {36}, journal = {Earthquake spectra : the professional journal of the Earthquake Engineering Research Institute}, number = {1_SUPPL}, publisher = {Sage Publishing}, address = {Thousand Oaks, CA}, issn = {8755-2930}, doi = {10.1177/8755293019899956}, pages = {274 -- 297}, year = {2020}, abstract = {Probabilistic assessment of seismic hazard and risk over a geographical region presents the modeler with challenges in the characterization of the site amplification that are not present in site-specific assessment. Using site-to-site residuals from a ground motion model fit to observations from the Japanese KiK-net database, correlations between measured local amplifications and mappable proxies such as topographic slope and geology are explored. These are used subsequently to develop empirical models describing amplification as a direct function of slope, conditional upon geological period. These correlations also demonstrate the limitations of inferring 30-m shearwave velocity from slope and applying them directly into ground motion models. Instead, they illustrate the feasibility of deriving spectral acceleration amplification factors directly from sets of observed records, which are calibrated to parameters that can be mapped uniformly on a regional scale. The result is a geologically calibrated amplification model that can be incorporated into national and regional seismic hazard and risk assessment, ensuring that the corresponding total aleatory variability reflects the predictive capability of the mapped site proxy.}, language = {en} } @article{KothaBazzurroPagani2018, author = {Kotha, Sreeram Reddy and Bazzurro, Paolo and Pagani, Marco}, title = {Effects of epistemic uncertainty in seismic hazard estimates on building portfolio losses}, series = {Earthquake spectra : the professional journal of the Earthquake Engineering Research Institute}, volume = {34}, journal = {Earthquake spectra : the professional journal of the Earthquake Engineering Research Institute}, number = {1}, publisher = {Earthquake Engineering Research Institute}, address = {Oakland}, issn = {8755-2930}, doi = {10.1193/020515EQS020M}, pages = {217 -- 236}, year = {2018}, abstract = {In catastrophe risk modeling, a defensible estimation of impact severity and its likelihood of occurrence to a portfolio of assets can only be made through a rigorous treatment of uncertainty and the consideration of multiple alternative models. This approach, however, requires repeating lengthy calculations multiple times. To limit the demand on computational time and resources, a frequent practice in the industry is to estimate the distribution of earthquake-induced portfolio losses using a simulated catalog of events from a single representative mean ground motion hazard model for the region. This simplified approach is faster but may provide biased estimates of the likelihood of occurrence of the large and infrequent losses that drive many risk mitigation decisions. Investigation through case studies of different portfolios of assets located in the San Francisco Bay Region shows the potential for both a bias in the mean loss estimates and an underestimation of their central 70\% interpercentile. We propose a simplified and computationally practical approach that reduces the bias in the mean portfolio loss estimates. This approach does not improve the estimate of the interpercentile range, however, a quantity of no direct practical use.}, language = {en} } @article{KothaCottonBindi2018, author = {Kotha, Sreeram Reddy and Cotton, Fabrice Pierre 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{KothaCottonBindi2019, author = {Kotha, Sreeram Reddy and Cotton, Fabrice Pierre 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{BindiKothaWeatherilletal.2018, author = {Bindi, Dino and Kotha, Sreeram Reddy and Weatherill, Graeme and Lanzano, Giovanni and Luzi, Lucia and Cotton, Fabrice Pierre}, 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{BindiPicozziSpallarossaetal.2019, author = {Bindi, Dino and Picozzi, Matteo and Spallarossa, Daniele and Cotton, Fabrice Pierre 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{KothaBindiCotton2017, author = {Kotha, Sreeram Reddy and Bindi, Dino and Cotton, Fabrice Pierre}, 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} } @article{KothaBindiCotton2017, author = {Kotha, Sreeram Reddy and Bindi, Dino and Cotton, Fabrice Pierre}, title = {Site-Corrected Magnitude- and Region-Dependent Correlations of Horizontal Peak Spectral Amplitudes}, 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/091416EQS150M}, pages = {1415 -- 1432}, year = {2017}, abstract = {Empirical correlations of horizontal peak spectral amplitudes (PSA) are modeled using the total-residuals obtained in a ground motion prediction equation (GMPE) regression. Recent GMPEs moved toward partially non-ergodic region-and site-specific predictions, while the residual correlation models remained largely ergodic. Using mixed-effects regression, we decompose the total-residuals of a pan-European GMPE into between-event, between-site, and event-and-site corrected residuals to investigate the ergodicity in empirical PSA correlations. We first observed that the between-event correlations are magnitude-dependent, partially due to the differences in source spectra, and influence of stress-drop parameter on small and large events. Next, removing the between-site residuals from within-event residuals yields the event-and-site corrected residuals which are found to be region-dependent, possibly due to the regional differences in distance-decay of short period PSAs. Using our site-corrected magnitude- and region-dependent correlations, and the between-site residuals as empirical site-specific ground motion adjustments, we compute partially non-ergodic conditional mean spectra at four well-recorded sites in Europe and Middle Eastern regions.}, language = {en} }