Probabilistic seismic hazard assessment for South-Eastern France

The accurate evaluation and appropriate treatment of uncertainties is of primary importance in modern probabilistic seismic hazard assessment (PSHA). One of the objectives of the SIGMA project was to establish a framework to improve knowledge and data on two target regions characterized by low-to-moderate seismic activity. In this paper, for South-Eastern France, we present the final PSHA performed within the SIGMA project. A new earthquake catalogue for France covering instrumental and historical periods was used for the calculation of the magnitude-frequency distributions. The hazard model incorporates area sources, smoothed seismicity and a 3D faults model. A set of recently developed ground motion prediction equations (GMPEs) from global and regional data, evaluated as adequately representing the ground motion characteristics in the region, was used to calculate the hazard. The magnitude-frequency distributions, maximum magnitude, faults slip rate and style-of-faulting are considered as additional source of epistemic uncertainties. The hazard results for generic rock condition (Vs30 = 800 m/s) are displayed for 20 sites in terms of uniform hazard spectra at two return periods (475 years and 10,000 years). The contributions of the epistemic uncertainties in the ground motion characterizations and in the seismic source characterization to the total hazard uncertainties are analyzed. Finally, we compare the results with existing models developed at national scale in the framework of the first generation of models supporting the Eurocode 8 enforcement, (MEDD 2002 and AFPS06) and at the European scale (within the SHARE project), highlighting significant discrepancies at short return periods.     

Probabilistic seismic hazard assessment of the Pyrenean region

A unified probabilistic seismic hazard assessment (PSHA) for the Pyrenean region has been performed by an international team composed of experts from Spain and France during the Interreg IIIA ISARD project. It is motivated by incoherencies between the seismic hazard zonations of the design codes of France and Spain and by the need for input data to be used to define earthquake scenarios. A great effort was invested in the homogenisation of the input data. All existing seismic data are collected in a database and lead to a unified catalogue using a local magnitude scale. PSHA has been performed using logic trees combined with Monte Carlo simulations to account for both epistemic and aleatory uncertainties. As an alternative to hazard calculation based on seismic sources zone models, a zoneless method is also used to produce a hazard map less dependant on zone boundaries. Two seismogenic source models were defined to take into account the different interpretations existing among specialists. A new regional ground-motion prediction equation based on regional data has been proposed. It was used in combination with published ground-motion prediction equations derived using European and Mediterranean data. The application of this methodology leads to the definition of seismic hazard maps for 475- and 1,975-year return periods for spectral accelerations at periods of 0 (corresponding to peak ground acceleration), 0.1, 0.3, 0.6, 1 and 2 s. Median and percentiles 15% and 85% acceleration contour lines are represented. Finally, the seismic catalogue is used to produce a map of the maximum acceleration expected for comparison with the probabilistic hazard maps. The hazard maps are produced using a grid of 0.1°. The results obtained may be useful for civil protection and risk prevention purposes in France, Spain and Andorra.     

Assessment of ground motion variability and its effects on seismic hazard analysis: a case study for iceland

Probabilistic seismic hazard analysis (PSHA) generally relies on the basic assumption that ground motion prediction equations (GMPEs) developed for other similar tectonic regions can be adopted in the considered area. This implies that observed ground motion and its variability at considered sites could be modelled by the selected GMPEs. Until now ground-motion variability has been taken into account in PSHA by integrating over the standard deviation reported in GMPEs, which significantly affects estimated ground motions, especially at very low probabilities of exceedance. To provide insight on this issue, ground-motion variability in the South Iceland Seismic Zone (SISZ), where many ground-motion records are available, is assessed. Three statistical methods are applied to separate the aleatory variability into source (inter-event), site (inter-site) and residual (intra-event and intra-site) components. Furthermore, the current PSHA procedure that makes the ergodic assumption of equality between spatially and temporal variability is examined. In contrast to the ergodic assumption, several recent studies show that the observed ground-motion variability at an individual location is lower than that implied by the standard deviation of a GMPE. This could imply a mishandling of aleatory uncertainty in PSHA by ignoring spatial variability and by mixing aleatory and epistemic uncertainties in the computation of sigma. Station correction coefficients are introduced in order to capture site effects at different stations. The introduction of the non-ergodic assumption in PSHA leads to larger epistemic uncertainty, although this is not the same as traditional epistemic uncertainty modelled using different GMPEs. The epistemic uncertainty due to the site correction coefficients (i.e. mean residuals) could be better constrained for future events if more information regarding the characteristics of these seismic sources and path dependence could be obtained.     

A probabilistic seismic hazard assessment for the Turkish territory: part II—fault source and background seismicity model

Over the years, several local and regional seismic hazard studies have been conducted for the estimation of the seismic hazard in Turkey using different statistical processing tools for instrumental and historical earthquake data and modeling the geologic and tectonic characteristics of the region. Recently developed techniques, increased knowledge and improved databases brought the necessity to review the national active fault database and the compiled earthquake catalogue for the development of a national earthquake hazard map. A national earthquake strategy and action plan were conceived and accordingly with the collaboration of the several institutions and expert researchers, the Revision of Turkish Seismic Hazard Map Project (UDAP-Ç-13-06) was initiated, and finalized at the end of 2014. The scope of the project was confined to the revision of current national seismic hazard map, using the state of the art technologies and knowledge of the active fault, earthquake database, and ground motion prediction equations. The following two seismic source zonation models are developed for the probabilistic earthquake hazard analysis: (1) Area source model, (2) Fault and spatial smoothing seismic source model (FSBCK). In this study, we focus on the development and the characterization of the Fault Source model, the background spatially smoothed seismicity model and intrinsic uncertainty on the earthquake occurrence-rates-estimation. Finally, PSHA results obtained from the fault and spatial smoothed seismic source model are presented for 43, 72, 475 and 2475 years return periods (corresponding to 69, 50, 10, and 2% probability of exceedance in 50 years) for PGA and 5% damped spectral accelerations at 0.2 and 1.0 s.     

Application-driven ground motion prediction equation for seismic hazard assessments in non-cratonic moderate-seismicity areas

We present a ground motion prediction equation (GMPE) for probabilistic seismic hazard assessments (PSHA) in low-to-moderate seismicity areas, such as Germany. Starting from the NGA-West2 flat-file (Ancheta et al. in Earthquake Spectra 30:989–1005, 2014), we develop a model tailored to the hazard application in terms of data selection and implemented functional form. In light of such hazard application, the GMPE is derived for hypocentral distance (along with the Joyner-Boore one), selecting recordings at sites with vs30 ≥ 360 m/s, distances within 300 km, and magnitudes in the range 3 to 8 (being 7.4 the maximum magnitude for the PSHA in the target area). Moreover, the complexity of the considered functional form is reflecting the availability of information in the target area. The median predictions are compared with those from the NGA-West2 models and with one recent European model, using the Sammon’s map constructed for different scenarios. Despite the simplification in the functional form, the assessed epistemic uncertainty in the GMPE median is of the order of those affecting the NGA-West2 models for the magnitude range of interest of the hazard application. On the other hand, the simplification of the functional form led to an increment of the apparent aleatory variability. In conclusion, the GMPE developed in this study is tailored to the needs for applications in low-to-moderate seismic areas and for short return periods (e.g., 475 years); its application in studies where the hazard is involving magnitudes above 7.4 and for long return periods is not advised.     

Seismic hazard map of the Middle East

The collaborative project Earthquake Model of the Middle East (EMME, 2010–2015) brought together scientists and engineers from the leading research institutions in the region and delivered state-of-the-art seismic hazard assessment covering Afghanistan, Armenia, Azerbaijan, Cyprus, Georgia, Iran, Iraq, Jordan, Lebanon, Palestine, Pakistan, Syria and Turkey. Their efforts have been materialized in the first homogenized seismic hazard model comprising earthquake catalogues, mapped active faults, strong motions databank, ground motion models and the estimated ground motion values for various intensity measure types and relevant return periods (e.g. 475–5000 years). The reference seismic hazard map of the Middle East, depicts the mean values of peak ground acceleration with a 10% chance of exceedance in 50 years, corresponding to a mean return period of 475 years. A full resolution poster is provided with this contribution.     

Seismic hazard assessment for Iceland in terms of macroseismic intensity using a site approach

We present the results of probabilistic seismic hazard assessment for Iceland in the framework of the EU project UPStrat-MAFA using the so-called site approach implemented in the SASHA computational code. This approach estimates seismic hazard in terms of macroseismic intensity by basically relying on local information about documented effects of past seismic events in the framework of a formally coherent and complete treatment of intensity data. In the case of Iceland, due to the lack of observed intensities for past earthquakes, local seismic histories were built using indirect macroseismic estimates deduced from epicentral information through an empirical attenuation relationship in probabilistic form. Seismic hazard was computed for four exceedance probabilities for an exposure time of 50 years, equivalent to average return periods of 50, 200, 475 and 975 years. For some localities, further return periods were examined and deaggregation analysis was performed. Results appear significantly different from previous seismic hazard maps, though just a semi-qualitative comparison is possible because of the different shaking measure considered (peak ground acceleration versus intensity), and the different computational methodology and input data used in these studies.     

Effects of seismic declustering on seismic hazard assessment: a sensitivity study using the Turkish earthquake catalogue

Earthquakes trigger other earthquakes and build up clusters in space and time that in turn create a bias in seismic catalogues. Therefore, declustering is considered as a prerequisite in seismic studies, particularly for probabilistic seismic hazard analysis, not only to eliminate the bias but also to decouple mainshocks and triggered events. However, a declustering process is not a straightforward task due to the complex nature of earthquake phenomena. There exist several declustering methods that mostly employ subjective rules to distinguish between background seismicity and offsprings. Eventually, the final declustered catalogues usually deviate significantly according to the employed method. This issue is raising some concerns, such as how to select the most suitable declustering algorithm, or to assess how this selection affects seismic hazard assessment. In consequence, the main goal of this paper is to quantify the sensitivity of seismic hazard assessments to different declustering techniques. Accordingly, the recently compiled Turkish earthquake catalogue was declustered by making use of three declustering algorithms. A total of six declustered catalogues, two catalogues per method, one by implementing the default input parameters, and one by altering the free input parameters of the employed methods, were produced. The clusters of selected earthquakes were studied in terms of the spatial–temporal distribution of earthquake sequences. A sensitivity analysis was conducted through the major steps of seismic hazard assessment for Istanbul metropolitan city. The seismicity of Istanbul and surroundings was modeled on the basis of four areal source zones. Comparative studies showed that, while the selected declustering algorithm did not significantly affect the completeness periods of moderate to large size earthquakes, it considerably altered those of small magnitude events (e.g. M_w 4.3–5.2) and consequently the recurrence parameters of the source zones. Depending on the declustering algorithm and input parameters, the activity rate was observed to vary up to a factor of two. The differences in the declustered catalogues obtained from different declustering approaches resulted in considerable variations in seismic hazard estimations. The hazard maps at return periods of 475 and 2475 years indicated that peak ground acceleration values may vary up to 20% at some locations. Moreover, the differences in 5% damped elastic spectral accelerations at T = 0.2 for the return periods of 475 and 2475 years are about 18 and 12%, respectively, on the southern shores of Istanbul where the highest hazard levels are observed.     

Seismic hazard assessment for the Itoiz dam site (Western Pyrenees,Spain)

This work summarises the seismic hazard analysis performed for the complete characterisation of strong ground-motion at the site of the Itoiz dam (Western Pyrenees, Spain). The hazard analysis includes the compilation of a composite catalogue from French and Spanish agencies, the definition of an original hybrid seismogenic source model (including zones and major faults) and the selection of ground motion prediction equations (GMPEs). Hazard results are provided as hazard curves and acceleration response spectra on rock for the 1000- and 5000-year return periods, which correspond respectively to the operating basis earthquake (OBE) and safety evaluation earthquake (SEE). The impact of truncating GMPEs at a number of standard deviations (epsilon) has been found not critical here for the return periods targeted. Subsequently, an analysis of the contribution of each source to total hazard and a hazard disaggregation analysis are performed in order to establish the earthquake-source parameters for both the OBE and SEE scenarios consistently with the seismotectonics of the region. The European Strong Motion database is then searched and a selection of records is proposed for each of the scenarios. Our results suggest that seismic hazard in the region is underestimated by the official Spanish seismic hazard map included in the current version of the code (NCSE-02), which is the reference document for the definition of seismic actions for dam projects in the whole Pyrenees.     

Application of Gumbel I and Monte Carlo methods to assess seismic hazard in and around Pakistan

A proper assessment of seismic hazard is of considerable importance in order to achieve suitable building construction criteria. This paper presents probabilistic seismic hazard assessment in and around Pakistan (23° N–39° N; 59° E–80° E) in terms of peak ground acceleration (PGA). Ground motion is calculated in terms of PGA for a return period of 475 years using a seismogenic-free zone method of Gumbel’s first asymptotic distribution of extreme values and Monte Carlo simulation. Appropriate attenuation relations of universal and local types have been used in this study. The results show that for many parts of Pakistan, the expected seismic hazard is relatively comparable with the level specified in the existing PGA maps.     

Ground-motion characterization for the probabilistic seismic hazard assessment in Turkey

This study describes the methodology implemented to establish the ground-motion logic-tree for national probabilistic seismic hazard map of Turkey for shallow active crustal regions. The presented procedure provides quantitative information to guide the hazard experts while establishing the logic tree to capture the epistemic uncertainty in ground-motion characterization. It uses non-data-driven and data-driven testing methods to identify and rank candidate ground-motion prediction equations (GMPEs) under a specific ground-motion database. The candidate GMPEs are subjected to visual inspection and are classified into center, body and range (CBR) spectral estimates for a proper consideration of epistemic uncertainty. The GMPEs classified into CBR are then used in a suite of seismic hazard sensitivity analysis to establish the most suitable GMPE logic-tree whose spectral estimates are not biased by any one of the GMPEs in the logic-tree structure. The sensitivity analysis considers normalized spectral ordinates and is not manipulated by the spectral amplitudes. The proposed procedure is inherited from the relevant studies of the Earthquake Model of the Middle East (EMME; www.efehr.org:8080/jetspeed/portal/emme.psml) regional seismic hazard project. This paper also highlights the similarities and differences in ground-motion characterization between EMME and our approach.     

余震时空丛集对概率地震危险性分析的影响

本文以东北、华北及川滇地区为例,系统研究了余震时空丛集对概率地震危险性分析的影响.采用基于传染型余震序列模型(ETAS)的蒙特卡罗模拟方法,模拟了包含余震和不包含余震的两套地震序列,然后以模拟地震目录为基础输入,采用基于空间光滑地震活动性模型的地震危险性分析方法计算了两套地震危险性结果——PGA(Peak Ground Acceleration,峰值加速度),通过分析比较这两套PGA的绝对差值和相对差值来研究余震时空丛集对概率地震危险性分析的影响.研究结果表明余震对50年超越概率10%地震危险性计算结果的影响均值为6%左右,最大可达10%,并且随着超越概率水平的提高,余震影响也越大.弱地震活动区余震对概率地震危险性分析的影响要高于强地震活动区.研究结果还进一步揭示两套PGA结果绝对差值的最大值约为15 cm·s~(-2),且出现在高PGA区,这意味着余震对概率地震危险性计算结果不会产生显著影响.因此在地震区划或一般性地震危险性分析中可考虑不用删除余震.     

Spatial sensitivity of seismic hazard results to different background seismic activity and temporal earthquake occurrence models

Spatial sensitivity of seismic hazard results to different models with respect to background seismic activity and earthquake occurrence in time is investigated. For the contribution of background seismic activity to seismic hazard, background area source with uniform seismicity and spatially smoothed seismicity models are taken into consideration. For the contribution of faults, through characteristic earthquakes, both the memoryless Poisson and the time dependent renewal models are utilized. A case study, involving the assessment of seismic hazard for the Bursa province in Turkey, is conducted in order to examine quantitatively the influence of these models on seismic hazard results. The spatial variation of the difference in Peak Ground Acceleration (PGA) values obtained from these different models is presented in the form of difference maps for return periods of 475 and 2475 years. Best estimate seismic hazard maps for PGA and Spectral Accelerations (SA) at 0.2 and 1.0 s are obtained by using the logic tree method.     

地震危险性概率分析方法:存在的问题和纠正

地震危险性概率分析(PSHA)是目前最广泛应用于地震灾害与风险性评估的方法。然而它在计算中却存在着一个错误:把强地面运动衰减关系(一个函数)的条件超越概率等同于强地面运动误差(一个变量)的超越概率。这个错误导致了运用强地面运动误差(空间分布特征)去外推强地面运动的发生(时间分布特征)或称之为遍历性假设,同时也造成了对PSHA理解和应用上的困难。本文推导出新的灾害计算方法(称之为KY-PSHA)来纠正这种错误。     

Site specific probabilistic seismic hazard analysis at Dubai Creek on the west coast of UAE

A probabilistic seismic hazard analysis (PSHA) was conducted to establish the hazard spectra for a site located at Dubai Creek on the west coast of the United Arab Emirates (UAE). The PSHA considered all the seismogenic sources that affect the site, including plate boundaries such as the Makran subduction zone, the Zagros fold-thrust region and the transition fault system between them; and local crustal faults in UAE. PSHA indicated that local faults dominate the hazard. The peak ground acceleration (PGA) for the 475-year return period spectrum is 0.17 g and 0.33 g for the 2,475-year return period spectrum. The hazard spectra are then employed to establish rock ground motions using the spectral matching technique.     

Toward a ground-motion logic tree for probabilistic seismic hazard assessment in Europe

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?C3247, 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.     

A probabilistic seismic hazard assessment for the Turkish territory—part I: the area source model

The seismic zoning map of Turkey that is used in connection with the national seismic design code (versions issued both in 1997 and 2007) is based on a probabilistic seismic hazard assessment study conducted more than 20 years ago (Gülkan et al. in En son verilere göre haz?rlanan Türkiye deprem bölgeleri haritas?, Report No: METU/EERC 93-1, 1993). In line with the efforts for the update of the seismic design code, the need aroused for an updated seismic hazard map, incorporating recent data and state-of-the-art methodologies and providing ground motion parameters required for the construction of the design spectra stipulated by the new Turkish Earthquake Design Code. Supported by AFAD (Disaster and Emergency Management Authority of Turkey), a project has been conducted for the country scale assessment of the seismic hazard by probabilistic methods. The present paper describes the probabilistic seismic hazard assessment study conducted in connection with this project, incorporating in an area source model, all recently compiled data on seismicity and active faulting, and using a set of recently developed ground motion prediction equations, for both active shallow crustal and subduction regimes, evaluated as adequately representing the ground motion characteristics in the region. The area sources delineated in the model are fully parameterized in terms of maximum magnitude, depth distribution, predominant strike and dip angles and mechanism of possible ruptures. Resulting ground motion distributions are quantified and presented for PGA and 5 % damped spectral accelerations at T = 0.2 and 1.0 s, associated with return periods of 475 and 2475 years. The full set of seismic hazard curves was also made available for the hazard computation sites. The second part of the study, which is based on a fault source and smoothed seismicity model is covered in Demircioglu et al. in Bull Earthq Eng, (2016).     

Comparing predicted and observed ground motions from subduction earthquakes in the Lesser Antilles

This brief article presents a quantitative analysis of the ability of eight published empirical ground-motion prediction equations (GMPEs) for subduction earthquakes (interface and intraslab) to estimate observed earthquake ground motions on the islands of the Lesser Antilles (specifically Guadeloupe, Martinique, Trinidad, and Dominica). In total, over 300 records from 22 earthquakes from various seismic networks are used within the analysis. It is found that most of the GMPEs tested perform poorly, which is mainly due to a larger variability in the observed ground motions than predicted by the GMPEs, although two recent GMPEs derived using Japanese strong-motion data provide reasonably good predictions. Analyzing separately the interface and intraslab events does not significant modify the results. Therefore, it is concluded that seismic hazard assessments for this region should use a variety of GMPEs in order to capture this large epistemic uncertainty in earthquake ground-motion prediction for the Lesser Antilles.     

Seismic hazard assessment of the Province of Murcia (SE Spain): analysis of source contribution to hazard

A probabilistic seismic hazard assessment of the Province of Murcia in terms of peak ground acceleration (PGA) and spectral accelerations [SA(T)] is presented in this paper. In contrast to most of the previous studies in the region, which were performed for PGA making use of intensity-to-PGA relationships, hazard is here calculated in terms of magnitude and using European spectral ground-motion models. Moreover, we have considered the most important faults in the region as specific seismic sources, and also comprehensively reviewed the earthquake catalogue. Hazard calculations are performed following the Probabilistic Seismic Hazard Assessment (PSHA) methodology using a logic tree, which accounts for three different seismic source zonings and three different ground-motion models. Hazard maps in terms of PGA and SA(0.1, 0.2, 0.5, 1.0 and 2.0 s) and coefficient of variation (COV) for the 475-year return period are shown. Subsequent analysis is focused on three sites of the province, namely, the cities of Murcia, Lorca and Cartagena, which are important industrial and tourism centres. Results at these sites have been analysed to evaluate the influence of the different input options. The most important factor affecting the results is the choice of the attenuation relationship, whereas the influence of the selected seismic source zonings appears strongly site dependant. Finally, we have performed an analysis of source contribution to hazard at each of these cities to provide preliminary guidance in devising specific risk scenarios. We have found that local source zones control the hazard for PGA and SA(T ≤ 1.0 s), although contribution from specific fault sources and long-distance north Algerian sources becomes significant from SA(0.5 s) onwards.