@article{BeauvalHainzlScherbaum2006, author = {Beauval, C{\´e}line and Hainzl, Sebastian and Scherbaum, Frank}, title = {Probabilistic seismic hazard estimation in low-seismicity regions considering non-Poissonian seismic occurrence}, issn = {0956-540X}, doi = {10.1111/j.1365-246X.2006.02863.x}, year = {2006}, abstract = {In low-seismicity regions, such as France or Germany, the estimation of probabilistic seismic hazard must cope with the difficult identification of active faults and with the low amount of seismic data available. Since the probabilistic hazard method was initiated, most studies assume a Poissonian occurrence of earthquakes. Here we propose a method that enables the inclusion of time and space dependences between earthquakes into the probabilistic estimation of hazard. Combining the seismicity model Epidemic Type Aftershocks-Sequence (ETAS) with a Monte Carlo technique, aftershocks are naturally accounted for in the hazard determination. The method is applied to the Pyrenees region in Southern France. The impact on hazard of declustering and of the usual assumption that earthquakes occur according to a Poisson process is quantified, showing that aftershocks contribute on average less than 5 per cent to the probabilistic hazard, with an upper bound around 18 per cent}, language = {en} } @article{HainzlScherbaumBeauval2006, author = {Hainzl, Sebastian and Scherbaum, Frank and Beauval, C{\´e}line}, title = {Estimating background activity based on interevent-time distribution}, issn = {0037-1106}, doi = {10.1785/0120050053}, year = {2006}, abstract = {The statistics of time delays between successive earthquakes has recently been claimed to be universal and to show the existence of clustering beyond the duration of aftershock bursts. We demonstrate that these claims are unjustified. Stochastic simulations with Poissonian background activity and triggered Omori-type aftershock sequences are shown to reproduce the interevent-time distributions observed on different spatial and magnitude scales in California. Thus the empirical distribution can be explained without any additional long-term clustering. Furthermore, we find that the shape of the interevent-time distribution, which can be approximated by the gamma distribution, is determined by the percentage of main-shocks in the catalog. This percentage can be calculated by the mean and variance of the interevent times and varies between 5\% and 90\% for different regions in California. Our investigation of stochastic simulations indicates that the interevent-time distribution provides a nonparametric reconstruction of the mainshock magnitude-frequency distribution that is superior to standard declustering algorithm}, language = {en} } @article{BeauvalScotti2004, author = {Beauval, C{\´e}line and Scotti, O.}, title = {Quantifying sensitivities of PSHA for France to earthquake catalog uncertainties, truncation of ground-motion variability, and magnitude limits}, issn = {0037-1106}, year = {2004}, abstract = {The results of this study clearly identify four key parameters controlling the estimation of probabilistic seismic hazard assessment (PSHA) in France in the framework of the Cornell-McGuire method. Results in terms of peak ground acceleration demonstrate the equally high impact, at all return periods, of the choice of truncation of the predicted ground-motion distribution (at + 2sigma) and of the choice between two different magnitude-intensity correlations. The choice of minimum magnitude (3.5/4.5) on hazard estimates can have an important impact at small return periods (<1000 years), whereas the maximum magnitude (6.5/7.0), on the other hand, is not a key parameter even at large return periods (10,000 years). This hierarchy of impacts is maintained at lower frequencies down to 5 Hz. Below 5 Hz, the choice of the maximum magnitude has a much greater impact, whereas the impact due to the choice of the minimum magnitude disappears. Moreover, variability due to catalog uncertainties is also quantified; these uncertainties that underly all hazard results can engender as high a variability as the controlling parameters. Parameter impacts, calculated at the centers of each source zone, show a linear trend with the seismicity models of the zone, demonstrating the lack of contributions coming from neighboring zones. Indeed, the region of influence that contributes to the PSHA estimate at a given site decreases with increasing return periods. The resulting overall variability in hazard estimates due to input uncertainties is quantified through a logic tree, obtained coefficients of variation vary between 10\% and 20\%. Until better physical models are obtained, the uncertainty on hazard estimates may be reduced by working on an appropriate magnitude-intensity correlation}, language = {en} } @article{BeauvalHainzlScherbaum2006, author = {Beauval, Celine and Hainzl, Sebastian and Scherbaum, Frank}, title = {The impact of the spatial uniform distribution of seismicity on probabilistic seismic-hazard estimation}, series = {Bulletin of the Seismological Society of America}, volume = {96}, journal = {Bulletin of the Seismological Society of America}, number = {6}, publisher = {GeoScienceWorld}, address = {Alexandria, Va.}, issn = {0037-1106}, doi = {10.1785/0120060073}, pages = {2465 -- 2471}, year = {2006}, abstract = {The first step in the estimation of probabilistic seismic hazard in a region commonly consists of the definition and characterization of the relevant seismic sources. Because in low-seismicity regions seismicity is often rather diffuse and faults are difficult to identify, large areal source zones are mostly used. The corresponding hypothesis is that seismicity is uniformly distributed inside each areal seismic source zone. In this study, the impact of this hypothesis on the probabilistic hazard estimation is quantified through the generation of synthetic spatial seismicity distributions. Fractal seismicity distributions are generated inside a given source zone and probabilistic hazard is computed for a set of sites located inside this zone. In our study, the impact of the spatial seismicity distribution is defined as the deviation from the hazard value obtained for a spatially uniform seismicity distribution. From the generation of a large number of synthetic distributions, the correlation between the fractal dimension D and the impact is derived. The results show that the assumption of spatially uniform seismicity tends to bias the hazard to higher values. The correlation can be used to determine the systematic biases and uncertainties for hazard estimations in real cases, where the fractal dimension has been determined. We apply the technique in Germany (Cologne area) and in France (Alps).}, 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{BeauvalTasanLaurendeauetal.2012, author = {Beauval, Celine and Tasan, Hilal and Laurendeau, Aurore and Delavaud, Elise and Cotton, Fabrice and Gueguen, Philippe and K{\"u}hn, Nicolas}, title = {On the testing of ground-motion prediction equations against small-magnitude data}, series = {Bulletin of the Seismological Society of America}, volume = {102}, journal = {Bulletin of the Seismological Society of America}, number = {5}, publisher = {Seismological Society of America}, address = {El Cerrito}, issn = {0037-1106}, doi = {10.1785/0120110271}, pages = {1994 -- 2007}, year = {2012}, abstract = {Ground-motion prediction equations (GMPE) are essential in probabilistic seismic hazard studies for estimating the ground motions generated by the seismic sources. In low-seismicity regions, only weak motions are available during the lifetime of accelerometric networks, and the equations selected for the probabilistic studies are usually models established from foreign data. Although most GMPEs have been developed for magnitudes 5 and above, the minimum magnitude often used in probabilistic studies in low-seismicity regions is smaller. Disaggregations have shown that, at return periods of engineering interest, magnitudes less than 5 may be contributing to the hazard. This paper presents the testing of several GMPEs selected in current international and national probabilistic projects against weak motions recorded in France (191 recordings with source-site distances up to 300 km, 3:8 <= M-w <= 4:5). The method is based on the log-likelihood value proposed by Scherbaum et al. (2009). The best-fitting models (approximately 2:5 <= LLH <= 3:5) over the whole frequency range are the Cauzzi and Faccioli (2008), Akkar and Bommer (2010), and Abrahamson and Silva (2008) models. No significant regional variation of ground motions is highlighted, and the magnitude scaling could be the predominant factor in the control of ground-motion amplitudes. Furthermore, we take advantage of a rich Japanese dataset to run tests on randomly selected low-magnitude subsets, and confirm that a dataset of similar to 190 observations, the same size as the French dataset, is large enough to obtain stable LLH estimates. Additionally we perform the tests against larger magnitudes (5-7) from the Japanese dataset. The ranking of models is partially modified, indicating a magnitude scaling effect for some of the models, and showing that extrapolating testing results obtained from low-magnitude ranges to higher magnitude ranges is not straightforward.}, language = {en} } @article{YepesAudinAlvaradoetal.2016, author = {Yepes, Hugo and Audin, Laurence and Alvarado, Alexandra and Beauval, Celine and Aguilar, Jorge and Font, Yvonne and Cotton, Fabrice}, title = {A new view for the geodynamics of Ecuador: Implication in seismogenic source definition and seismic hazard assessment}, series = {Tectonics}, volume = {35}, journal = {Tectonics}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0278-7407}, doi = {10.1002/2015TC003941}, pages = {1249 -- 1279}, year = {2016}, abstract = {A new view of Ecuador's complex geodynamics has been developed in the course of modeling seismic source zones for probabilistic seismic hazard analysis. This study focuses on two aspects of the plates' interaction at a continental scale: (a) age-related differences in rheology between Farallon and Nazca plates—marked by the Grijalva rifted margin and its inland projection—as they subduct underneath central Ecuador, and (b) the rapidly changing convergence obliquity resulting from the convex shape of the South American northwestern continental margin. Both conditions satisfactorily explain several characteristics of the observed seismicity and of the interseismic coupling. Intermediate-depth seismicity reveals a severe flexure in the Farallon slab as it dips and contorts at depth, originating the El Puyo seismic cluster. The two slabs position and geometry below continental Ecuador also correlate with surface expressions observable in the local and regional geology and tectonics. The interseismic coupling is weak and shallow south of the Grijalva rifted margin and increases northward, with a heterogeneous pattern locally associated to the Carnegie ridge subduction. High convergence obliquity is responsible for the North Andean Block northeastward movement along localized fault systems. The Cosanga and Pallatanga fault segments of the North Andean Block-South American boundary concentrate most of the seismic moment release in continental Ecuador. Other inner block faults located along the western border of the inter-Andean Depression also show a high rate of moderate-size earthquake production. Finally, a total of 19 seismic source zones were modeled in accordance with the proposed geodynamic and neotectonic scheme.}, language = {en} }