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.     

Probabilistic seismic hazard analysis: Early history

Probabilistic seismic hazard analysis (PSHA) is the evaluation of annual frequencies of exceedence of ground motion levels (typically designated by peak ground acceleration or by spectral accelerations) at a site. The result of a PSHA is a seismic hazard curve (annual frequency of exceedence vs ground motion amplitude) or a uniform hazard spectrum (spectral amplitude vs structural period, for a fixed annual frequency of exceedence). Analyses of this type were first conceived in the 1960s and have become the basis for the seismic design of engineered facilities ranging from common buildings designed according to building codes to critical facilities such as nuclear power plants. This Historical Note traces the early history of PSHA. Copyright © 2007 John Wiley & Sons, Ltd.     

Comparison of conventional and Monte Carlo simulation-based probabilistic seismic hazard analyses for Shiraz city,southern Iran

Estimation of ground-motion amplitudes of different hazard levels is of paramount importance in planning of urban development of any metropolis. Such estimation can be computed through a probabilistic seismic hazard analysis (PSHA). This paper concentrates on the PSHA of an area located in Shiraz city, southern Iran. The area includes whole of Shiraz city (i.e., one of the largest and most populous cities of Iran) and its outskirts. Conventional and Monte Carlo simulation-based approaches are utilized to perform the PSHA of the studied area. Two areal seismic source models are delineated, and thence seismicity parameters of all zones associated with their corresponding uncertainties are computed. Uncertainties in ground-motion prediction are accounted for via three ground-motion prediction equations (GMPEs) within the logic tree framework. These GMPEs are applied to estimate bedrock ground shaking (V_s30?=?760 m/s) for several return periods (i.e., 75, 475, 975, and 2475 years). In general, the results of the two abovementioned PSHA approaches show relatively similar results. However, the Monte Carlo simulation-based approach overpredicts bedrock spectral accelerations at periods of 0.4–2.5 s compared to the conventional PSHA approach for return periods of 475, 975, and 2475 years.     

Defining a model of 3D seismogenic sources for Seismic Hazard Assessment applications: The case of central Apennines (Italy)

Geology-based methods for Probabilistic Seismic Hazard Assessment (PSHA) have been developing in Italy. These methods require information on the geometric, kinematic and energetic parameters of the major seismogenic faults. In this paper, we define a model of 3D seismogenic sources in the central Apennines of Italy. Our approach is mainly structural-seismotectonic: we integrate surface geology data (trace of active faults, i.e. 2D features) with seismicity and subsurface geological–geophysical data (3D approach). A fundamental step is to fix constraints on the thickness of the seismogenic layer and deep geometry of faults: we use constraints from the depth distribution of aftershock zones and background seismicity; we also use information on the structural style of the extensional deformation at crustal scale (mainly from seismic reflection data), as well as on the strength and behaviour (brittle versus plastic) of the crust by rheological profiling. Geological observations allow us to define a segmentation model consisting of major fault structures separated by first-order (kilometric scale) structural-geometric complexities considered as likely barriers to the propagation of major earthquake ruptures. Once defined the 3D fault features and the segmentation model, the step onward is the computation of the maximum magnitude of the expected earthquake (M _max). We compare three different estimates of M _max: (1) from association of past earthquakes to faults; (2) from 3D fault geometry and (3) from geometrical estimate corrected by earthquake scaling laws. By integrating all the data, we define a model of seismogenic sources (seismogenic boxes), which can be directly used for regional-scale PSHA. Preliminary applications of PSHA indicate that the 3D approach may allow to hazard scenarios more realistic than those previously proposed.     

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.     

The effect of causal parameter bounds in PSHA‐based ground motion selection

In this paper the effect of causal parameter bounds (e.g. magnitude, source‐to‐site distance, and site condition) on ground motion selection, based on probabilistic seismic hazard analysis (PSHA) results, is investigated. Despite the prevalent application of causal parameter bounds in ground motion selection, present literature on the topic is cast in the context of a scenario earthquake of interest, and thus specific bounds for use in ground motion selection based on PSHA, and the implications of such bounds, is yet to be examined. Thirty‐six PSHA cases, which cover a wide range of causal rupture deaggregation distributions and site conditions, are considered to empirically investigate the effects of various causal parameter bounds on the characteristics of selected ground motions based on the generalized conditional intensity measure (GCIM) approach. It is demonstrated that the application of relatively ‘wide’ bounds on causal parameters effectively removes ground motions with drastically different characteristics with respect to the target seismic hazard and results in an improved representation of the target causal parameters. In contrast, the use of excessively ‘narrow’ bounds can lead to ground motion ensembles with a poor representation of the target intensity measure distributions, typically as a result of an insufficient number of prospective ground motions. Quantitative criteria for specifying bounds for general PSHA cases are provided, which are expected to be sufficient in the majority of problems encountered in ground motion selection for seismic demand analyses. Copyright © 2016 John Wiley & Sons, Ltd.     

Potential rupture surface model and its application on probabilistic seismic hazard analysis

Potential sources are simplified as point sources or linear sources in current probabilistic seismic hazard analysis (PSHA) methods. Focus size of large earthquakes is considerable, and fault rupture attitudes may have great influ-ence upon the seismic hazard of a site which is near the source. Under this circumstance, it is unreasonable to use the simplified potential source models in the PSHA, so a potential rupture surface model is proposed in this paper. Adopting this model, we analyze the seismic hazard near the Chelungpu fault that generated the Chi-Chi (Jiji) earthquake with magnitude 7.6 and the following conclusions are reached. 1 This model is reasonable on the base of focal mechanism, especially for sites near potential earthquakes with large magnitude; 2 The attitudes of poten-tial rupture surfaces have great influence on the results of probabilistic seismic hazard analysis and seismic zoning.     

Peak-over-threshold: Quantifying ground motion beyond design

In performance-based seismic design, as adopted by several building codes worldwide, the structural performance is verified against ground motions that have predetermined exceedance return periods at the site of interest. Such a return period is evaluated by means of probabilistic seismic hazard analysis (PSHA), and the corresponding ground motion is often represented by the uniform hazard spectrum (UHS). The structural performance for ground motions larger than those considered in this design approach is, typically, not explicitly controlled under the assumption that they are sufficiently rare. On one hand, this does not achieve uniform safety at sites characterized by different design ground motions corresponding to the same return period; on the other hand, exceedances of the design spectra are systematically observed over large areas, for example in Italy. The latter issue is because of the nature of UHS, the exceedance of which is likely-to-almost-certain when the construction site is in the epicentral area of moderate-to-high magnitude earthquakes (ie, the design spectrum may be not conservative at these locations), especially if PSHA is based on seismic source zones. The former is partially because of the systematic difference of ground motions for return periods larger than the design one at the different sites. Quantification of the expected ground motion given the exceedance of the design ground motions (ie, the recently introduced as the expected peak-over-threshold or POT) can be of help in quantitatively assessing these issues. In the study, a procedure to compute the POT distribution is derived first; second, POT spectra are introduced and used to help understanding why and how seismic structural reliability of code-conforming structures decreases as the seismic hazard of the site increases; third, expected and 95^th percentile POT maps are shown for Italy to discuss how much high hazard sites are exposed to much larger peak-over-threshold with respect to mid-hazard and low-hazard sites; finally the POT is discussed with respect to the slope of the hazard curve (in log-log scale) at the threshold, a known proxy for ground motion beyond design. All data presented in the maps are made available for the interested reader as a supplemental archive.     

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.     

Near‐source seismic hazard and design scenarios

Earthquakes damage engineering structures near, relatively to the rupture's size, to the source. In this region, the fault's dynamics affect ground motion propagation differently from site to site, resulting in systematic spatial variability known as directivity. Although a number of researches recommend that records with directivity‐related velocity pulses should be explicitly taken into account when defining design seismic action on structures, probabilistic seismic hazard analysis (PSHA), in its standard version, seems inadequate for the scope. In the study, it is critically reviewed why, from the structural engineering point of view, hazard assessment should account for near‐source effects (i.e., pulse‐like ground motions), and how this can be carried out adjusting PSHA analytically via introduction of specific terms and empirically calibrated models. Disaggregation analysis and design scenarios for near‐source PSHA are also formulated. The analytical procedures are then applied to develop examples of hazard estimates for sites close to strike–slip or dip–slip faults and to address differences with respect to the ordinary case, that is, when pulse‐like effects are not explicitly accounted for. Significant increase of hazard for selected spectral ordinates is found in all investigated cases; increments depend on the fault‐site configuration. Moreover, to address design scenarios for seismic actions on structures, disaggregation results are also discussed, along with limitations of current design spectra to highlight the pulse‐like effects of structural interest. Finally, an attempt to overcome these, by means of disaggregation‐based scenarios specific for the pulse occurrence case, is presented. Copyright © 2012 John Wiley & Sons, Ltd.     

Strong ground acceleration seismic hazard in Greece and neighboring regions

In an early paper [Tectonophysics 117 (1985) 259] seismic hazard in Greece was analyzed using a relatively homogeneous earthquake catalogue spanning 1900–1978 and a strong motion attenuation relationship adapted to use in Greece. Improved seismic hazard analyses are obtained here using Gumbel's asymptotic extreme value distribution applied to peak horizontal ground acceleration occurrence, but now taking into account the increased length and quality of earthquake catalogue data spanning 1900–1999 and the burgeoning information on earthquake strong motion data and attenuation relationships appropriate for Europe and, explicitly, Greece. Seismic acceleration hazard results tabulated for six cities reveal (e.g. using arbitrarily the 50-year p.g.a. with 90% probability of not being exceeded) changes of about 10% in the new calculated values: two cities show an increase and four a decrease. These are relatively small and reassuring adjustments.Inspection of the available attenuation relationships leads to a preference for the models of Theodulidis and Papazachos, particularly with the model modification to produce a ‘stiff soil’ site relationship, as these relationships explicitly exploit the Greek strong motion database. Isoacceleration maps are produced for Greece as a whole from each attenuation relationship inspected. The final set of maps based on the Theodulidis and Papazachos models provide a foundation for comparison with the Seismic Hazard Zones adopted in the New Greek Seismic Code where scope can be found to modify zone shape and the level at which p.g.a.s are set. It should be noted that the generation of the present isoacceleration maps is based on a seismogenic zone-free methodology, independent of any Euclidean zoning assumptions.     

Seismic Hazard Assessment for the Main Seismogenic Zones in North Algeria

—A seismogenic zoning for north Algeria is proposed. The characteristics of each zone are studied in detail and different statistical analyses are performed in order to assess the slope of the Gutenberg-Richter relationship, the mean annual activity rate of earthquakes, the average return period, the probability of non-exceedence during a time period and the maximum expected magnitude. The slope of the Gutenberg-Richter relationship determined by traditional methods, physical strain energy release and corresponding magnitude, together with a new approach based on both historical and instrumental data, are used to estimate the seismic hazard parameters.¶The seismic hazard parameters obtained described well the characteristics of each proposed seismogenic zone. The b-value of the Gutenberg-Richter relationship concurs especially well with the most common seismotectonic interpretation of each zone, for example 0.70 ± 0.06 at El Asnam region and 0.63 ± 0.09 around Algiers. The average return period obtained (approximately 300 years for magnitude 7.5 at El Asnam region, exceeding 1000 years in other regions), corroborate those obtained through Paleoseismic studies and others using a digital elevation model.     

Application of a Monte‐Carlo simulation approach for the probabilistic assessment of seismic hazard for geographically distributed portfolio

The conventional integral approach is very well established in probabilistic seismic hazard assessment (PSHA). However, Monte‐Carlo (MC) simulations can become an efficient and flexible alternative against conventional PSHA when more complicated factors (e.g. spatial correlation of ground shaking) are involved. This study aims at showing the implementation of MC simulation techniques for computing the annual exceedance rates of dynamic ground‐motion intensity measures (GMIMs) (e.g. peak ground acceleration and spectral acceleration). We use multi‐scale random field technique to incorporate spatial correlation and near‐fault directivity while generating MC simulations to assess the probabilistic seismic hazard of dynamic GMIMs. Our approach is capable of producing conditional hazard curves as well. We show various examples to illustrate the potential use of the proposed procedures in the hazard and risk assessment of geographically distributed structural systems. Copyright © 2015 John Wiley & Sons, Ltd.     

Seismic Hazard Assessment: Issues and Alternatives

Seismic hazard and risk are two very important concepts in engineering design and other policy considerations. Although seismic hazard and risk have often been used interchangeably, they are fundamentally different. Furthermore, seismic risk is more important in engineering design and other policy considerations. Seismic hazard assessment is an effort by earth scientists to quantify seismic hazard and its associated uncertainty in time and space and to provide seismic hazard estimates for seismic risk assessment and other applications. Although seismic hazard assessment is more a scientific issue, it deserves special attention because of its significant implication to society. Two approaches, probabilistic seismic hazard analysis (PSHA) and deterministic seismic hazard analysis (DSHA), are commonly used for seismic hazard assessment. Although PSHA has been proclaimed as the best approach for seismic hazard assessment, it is scientifically flawed (i.e., the physics and mathematics that PSHA is based on are not valid). Use of PSHA could lead to either unsafe or overly conservative engineering design or public policy, each of which has dire consequences to society. On the other hand, DSHA is a viable approach for seismic hazard assessment even though it has been labeled as unreliable. The biggest drawback of DSHA is that the temporal characteristics (i.e., earthquake frequency of occurrence and the associated uncertainty) are often neglected. An alternative, seismic hazard analysis (SHA), utilizes earthquake science and statistics directly and provides a seismic hazard estimate that can be readily used for seismic risk assessment and other applications.     

A unified framework for earthquake risk assessment of transportation networks and gross regional product

Natural threats like earthquakes, hurricanes or tsunamis have had serious impacts on communities. In the past, major earthquakes in the United States like Loma Prieta 1989, Northridge 1994, or recent events in Italy like L’Aquila 2009 or Emilia 2012 emphasized the importance of preparedness and awareness to reduce social impacts. In addition to that, earthquake damaged businesses dramatically reduced the gross regional product. Generating scenario earthquakes in a proper way is important to suitably assess the risk in bridge networks and social losses in terms of gross regional product reduction. Seismic hazard is traditionally assessed by means of probabilistic seismic hazard analysis (PSHA). Although PSHA well represents the hazard at a specific location it is not suitable for spatially distributed systems. Scenario earthquakes can overcome this problem; they represent the actual distribution of ground shaking for a spatially distributed system while being hazard consistent. In this work a methodology to generate scenario earthquakes has been proposed using a novel approach with the aim of being the basic step for investigating possible earthquake consequences in seismic areas and contributing to reduce losses.     

Building vulnerability and seismic risk analysis in the urban area of Mt. Etna volcano (Italy)

The tectonic system of the eastern flank of Mt. Etna volcano (Sicily, Italy) is the source of most of the strongest earthquakes occurring in the area over the last 205 years. A total of 12 events with epicentre intensities ≥VIII EMS have occurred at Mt. Etna, 10 of which were located on the eastern flank. This indicates a mean recurrence time of about 20 years. This area is highly urbanised, with many villages around the volcano at altitudes up to 700 m a.s.l. The southern and eastern flanks are particularly highly populated areas, with numerous villages very close to each other. The probabilistic seismic hazard due to local faults for Mt. Etna was calculated by adopting a site approach to seismic hazard assessment. Only the site histories of local volcano-tectonic earthquakes were considered, leaving out the effects due to strong regional earthquakes that occurred in north-eastern and south-eastern Sicily. The inventory used in this application refers to residential buildings. These data were extracted from the 1991 census of the Italian National Institute of Statistics, and are grouped according to the census sections. The seismic vulnerability of the elements at risk belonging to a given building typology is described by a vulnerability index, in accordance with a damage model based on macroseismic intensities. For the estimation of economic losses due to physical damage to buildings, an integrated impact indicator was used, which is equivalent to the lost building volume. The expected annualised economic earthquake losses were evaluated both in absolute and in relative terms, and were compared with the geographical distribution of seismic hazard and with similar evaluations of losses for other regions.     

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.     

The probabilistic seismic hazard assessment of Germany—version 2016, considering the range of epistemic uncertainties and aleatory variability

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 \(v_{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.