Stochastic ground-motion simulation of two Himalayan earthquakes: seismic hazard assessment perspective

The earthquakes in Uttarkashi (October 20, 1991, M _w 6.8) and Chamoli (March 8, 1999, M _w 6.4) are among the recent well-documented earthquakes that occurred in the Garhwal region of India and that caused extensive damage as well as loss of life. Using strong-motion data of these two earthquakes, we estimate their source, path, and site parameters. The quality factor (Q _β ) as a function of frequency is derived as Q _β (f) = 140f ^1.018. The site amplification functions are evaluated using the horizontal-to-vertical spectral ratio technique. The ground motions of the Uttarkashi and Chamoli earthquakes are simulated using the stochastic method of Boore (Bull Seismol Soc Am 73:1865–1894, 1983). The estimated source, path, and site parameters are used as input for the simulation. The simulated time histories are generated for a few stations and compared with the observed data. The simulated response spectra at 5% damping are in fair agreement with the observed response spectra for most of the stations over a wide range of frequencies. Residual trends closely match the observed and simulated response spectra. The synthetic data are in rough agreement with the ground-motion attenuation equation available for the Himalayas (Sharma, Bull Seismol Soc Am 98:1063–1069, 1998).     

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 uniform catalog of earthquakes for seismic hazard asesment in Iran

ＡｕｎｉｆｏｒｍｃａｔａｌｏｇｏｆｅａｒｔｈｑｕａｋｅｓｆｏｒｓｅｉｓｍｉｃｈａｚａｒｄａｓｅｓｍｅｎｔｉｎＩｒａｎＮｏｏｒｂａｋｈｓｈＭｉｒｚａｅｉ１，２）ＭＥＮＧ－ＴＡＮＧＡＯ１）（高孟谭）ＹＵＮ－ＴＡＩＣＨＥＮ１）（陈运泰）ＪＩＡＮＷＡＮＧ１...     

Treatment of ground-motion predictive relationships for the reference seismic hazard map of Italy

In the framework of the 2004 reference seismic hazard map of Italy the amplitude of the strong-motion (expressed in terms of Peak Horizontal Acceleration with 10% probability of non-exceedence in 50 years, referred to average hard ground conditions) was computed using different predictive relationships. Equations derived in Italy and in Europe from strong-motion data, as well as a set of weak and strong-motion based empirical predictive relationships were employed in a logic tree procedure, in order to capture the epistemic uncertainty affecting ground-motion attenuation. This article describes the adjustments and conversions required to eliminate the incompatibilities amongst the relations. Particularly significant are distance conversions and style-of-faulting adjustments, as well as the problems related to the use of regional relations, such as the selection of a reference depth, the quantification of random variability and the strong-motion prediction. Moreover, a regional attenuation relationship specific for volcanic areas was also employed, allowing a more realistic evaluation of seismic hazard, as confirmed by the attenuation of macroseismic intensities.     

Influence of ground-motion correlation on probabilistic assessments of seismic hazard and loss: sensitivity analysis

Recent studies have shown that the proper treatment of ground-motion variability and, particularly, the correlation of ground motion are essential for the estimation of the seismic hazard, damage and loss for distributed portfolios. In this work we compared the effects of variations in the between-earthquake correlation and in the site-to-site correlation on probabilistic estimations of seismic damage and loss for the extended objects (hypothetical portfolio) and critical elements (e.g. bridges) of a network. Taiwan Island has been chosen as a test case for this study because of relatively high seismicity and previous experience in earthquake hazard modelling. The hazard and loss estimations were performed using Monte Carlo approach on the basis of stochastic catalogues and random ground-motion fields. We showed that the influence of correlation on parameters of seismic hazard, characteristics of loss distribution and the probability of damage depend, on one hand, on level of hazard and probability level of interest (return period) and, on the other hand, the relative influence of each type of correlation is not equal.     

Updated seismic hazard assessment of Tunisia

Modern earthquake loss models make use of earthquake catalogs relevant to the seismic hazard assessment upon seismicity and seismotectonic analysis. The main objective of this paper is to investigate a recently compiled catalog (National Institute of Meteorology or INM catalog: 412-2011) and to generate seismic hazard maps through classical probabilistic seismic hazard assessment (PSHA) and smoothed-gridded seismicity models for Tunisia. It is now established with the local earthquake bulletin that the recent seismicity of Tunisia is sparse and moderate. Therefore, efforts must be undertaken to elaborate a robust hazard analysis for risk assessment and seismic design purposes. These recommendations follow the recently published reports by the World Bank that describe the seismic risk in Tunis City as being beyond a tolerable level with an MSK intensity level of VII. Some attempts were made during the past two decades to assess the seismic hazard for Tunisia and they have mostly failed to properly investigate the historical and instrumental seismicity catalog. This limitation also exists for the key aspect of epistemic and random uncertainties impact on the final seismic hazard assessment. This study also investigates new ground motion prediction equations suitable for use in Tunisia. The methodology applied herein uses, for the first time in PSHA of Tunisia, seismicity parameters integrated in logic tree framework to capture epistemic uncertainties through three different seismic source models. It also makes use of the recently released version of OpenQuake engine; an open-source tool for seismic hazard and risk assessment developed in the framework of the Global Earthquake Model.     

Update of likelihood-based ground-motion model selection for seismic hazard analysis in western central Europe

Scherbaum et al. [(2004) Bull Seismolo Soc Am 94(6): 2164–2185] proposed a likelihood-based approach to select and rank ground-motion models for seismic hazard analysis in regions of low-seismicity. The results of their analysis were first used within the PEGASOS project [Abrahamson et al. (2002), In Proceedings of the 12 ECEE, London, 2002, Paper no. 633] so far the only application of a probabilistic seismic hazard analysis (PSHA) in Europe which was based on a SSHAC Level 4 procedure [(Budnitz et al. 1997, Recommendations for PSHA: guidance on uncertainty and use of experts. No. NUREG/CR-6372-V1). The outcome of this project have generated considerable discussion (Klügel 2005, Eng Geol 78:285–307, 2005b) Eng Geol 78: 285–307, (2005c) Eng Geol 82: 79–85 Musson et al. (2005) Eng Geol 82(1): 43–55]; Budnitz et al. (2005), Eng Geol 78(3–4): 285–307], a central part of which is related to the issue of ground-motion model selection and ranking. Since at the time of the study by Scherbaum et al. [(2004.) Bull Seismolo Soc Am 94(6): 2164–2185], only records from one earthquake were available for the study area, here we test the stability of their results using more recent data. Increasing the data set from 12 records of one earthquake in Scherbaum et al. [(2004) Bull Seismolo Soc Am 94(6): 2164–2185] to 61 records of 5 earthquakes, which have mainly occurred since the publication of the original study, does not change the set of the three top-ranked ground-motion models [Abrahamson and Silva (1997) Seismolo Res Latt 68(1): 94–127; Lussou et al. (2001) J Earthquake Eng 5(1):13–33; Berge-Thierry et al. (2003) Bull Seismolog Soc Am 95(2): 377–389. Only for the lower-ranked models do we obtain modifications in the ranking order. Furthermore, the records from the Waldkirch earthquake (Dec, 5th, 2004, M _w = 4.9) enabled us to develop a new stochastic model parameter set for the application of Campbell’s [(2003) Bull Seismolo Soc Am 93(3): 1012–1033] hybrid empirical model to SW Germany and neighbouring regions.     

From regional- to local-scale seismic hazard assessment: examples from Southern Spain

An approach that relates results from a regional seismic hazard assessment study with local-scale site-effect characterizations in an area of low-to-moderate seismic activity such as Andalusia (southern Spain), is presented. Results of a previous probabilistic seismic hazard analysis of Andalusia on rock conditions are disaggregated to infer hazard controlling earthquakes for different target motions. A collection of controlling magnitude-distance pairs and the corresponding site-specific response spectra at main capital cities of the region are obtained. These spectra are first-order approximations to expected seismic actions required in local earthquake risk assessments. In addition, results of independent, local-scale studies developed in Almeria City (SE Andalusia) are used to derive an updated seismic zonation of the city. These include predominant soil period estimates and shear-wave velocity profiles at different locations. If a local seismic risk assessment study or an earthquake-resistant structural design is to be developed, it may be recommended the use of different seismic actions on sites characterized by distinct response to seismic shaking (as derived from the seismic zonation). The seismic action related to worst-case scenarios may be modeled through a hazard-consistent response spectrum, obtained by hazard disaggregation at the spectral acceleration with period matching the prevailing resonant period of the target site or structure.     

Application of SASHA to seismic hazard assessment for Portugal mainland

In the frame of the UPStrat-MAFA “Urban Disaster Prevention Strategies Using MAcroseismic Fields and FAult Sources” project, seismic hazard has been assessed in Portugal in terms of macroseismic intensity. Assessment has been performed by using a probabilistic approach based on the statistical analysis of local seismic histories (i.e., the record of seismic effects at each locality) performed by a new version of the SASHA code (D’Amico and Albarello in Res Lett 79(5):663–671, 2008) on purpose modified to account for this specific area of study. Local seismic histories are reconstructed by considering documented effects or indirect estimates deduced from epicentral information or numerical simulations. All these pieces of evidence are combined taking into account relevant uncertainty and statistical completeness of information locally available. Distribution of expected maximum intensity (i.e., the maximum intensity characterized by a fixed exceedance probability for a exposure time of 50 years) has been obtained and compared with the one deduced from alternative approaches.     

A Bayesian methodology to update the probabilistic seismic hazard assessment

This paper presents a Bayesian methodology for updating the seismic hazard curves. The methodology is based on the comparison of predictive exceedance rates of a fixed acceleration level (given by the seismic hazard curves) and the observed exceedance rates in some selected sites. The application of the methodology needs, firstly, the definition of a prior probabilistic seismic hazard assessment based in a logic tree. Each main branch corresponds to a probabilistic model of calculus of seismic hazard. The method considers that, initially (or a priori), the weights of all branches of the logic tree are equivalent. Secondly, the method needs to compile the observations in the region. They are introduced in a database containing the recorded acceleration data (during the instrumental period). Nevertheless, the instrumental period in stable zones (as France) shows only very low acceleration levels recorded during a short observation period. Then, a method to enlarge the REX (number of observations) is presented taking into account the historical data and defining “synthetic” accelerations in the sites of observation. The synthetic REX allows to expand the period of observation and to increase the acceleration thresholds used in the Bayesian updating process. The application of the Bayesian approach leads to a new and more objective definition of the weights of each branch of the logic tree and, therefore, to new seismic hazard curves (mean and centiles). The Bayesian approach doesn’t change the probabilistic models (seismic hazard curves). It only modifies the weights of each branch of the logic tree.     

Interpreting intraplate tectonics for seismic hazard: a UK historical perspective

It is notoriously difficult to construct seismic source models for probabilistic seismic hazard assessment in intraplate areas on the basis of geological information, and many practitioners have given up the task in favour of purely seismicity-based models. This risks losing potentially valuable information in regions where the earthquake catalogue is short compared to the seismic cycle. It is interesting to survey how attitudes to this issue have evolved over the past 30 years. This paper takes the UK as an example, and traces the evolution of seismic source models through generations of hazard studies. It is found that in the UK, while the earliest studies did not consider regional tectonics in any way, there has been a gradual evolution towards more tectonically based models. Experience in other countries, of course, may differ.     

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.     

Probabilistic seismic hazard assessment of Eastern Marmara Region

The objective of this study is to evaluate the seismic hazard in Eastern Marmara Region using an improved probabilistic seismic hazard assessment methodology. Two significant improvements over the previous seismic hazard assessment practices are accomplished in this study: advanced seismic source characterization models in terms of source geometry and recurrence relationships are developed, and improved global ground motion models (NGA-W1 models) are employed to represent the ground motion variability. Planar fault segments are defined and a composite magnitude distribution model is used for all seismic sources in the region to properly represent the characteristic behavior of the North Anatolian Fault without the need for an additional background zone. Multi-segment ruptures are considered using the rupture model proposed by the Working Group on California Earthquake Probabilities (2003). Events in the earthquake catalogue are attributed to the fault zones and scenario weights are determined by releasing the accumulated seismic energy. The uniform hazard spectra at 10 % probability of exceedance in 50 years hazard level for different soil conditions (soil and rock) are revealed for specific locations in the region (Adapazar?, Düzce, Gemlik, Izmit, Iznik and Sapanca). Hazard maps of the region for rock site conditions at the selected hazard levels are provided to allow the readers perform site-specific hazard assessment and develop site-specific design spectrum for local site conditions.     

Probabilistic seismic hazard assessment of Switzerland: best estimates and uncertainties

We present the results of a new genera tion of probabilistic seismic hazard assessment for Switzerland. This study replaces the previous intensity-based generation of national hazard maps of 1978. Based on a revised moment-magnitude earthquake catalog for Switzerland and the surrounding regions, covering the period 1300–2003, sets of recurrence parameters (a and b values, M _max ) are estimated. Information on active faulting in Switzerland is too sparse to be used as source model. We develop instead two models of areal sources. The first oriented towards capturing historical and instrumental seismicity, the second guided largely by tectonic principles and express ing the alterative view that seismicity is less stationary and thus future activity may occur in previously quiet regions. To estimate three alterna tive a and b value sets and their relative weighting, we introduce a novel approach based on the modified Akaike information criterion, which allows us to decide when the data in a zone deserves to be fitted with a zone-specific b value. From these input parameters, we simulate synthetic earthquake catalogs of one-million-year duration down to magnitude 4.0, which also reflect the difference in depth distribution between the Alpine Foreland and the Alps. Using a specific predictive spectral ground motion model for Switzerland, we estimate expected ground motions in units of the 5% damped acceleration response spectrum at frequencies of 0.5–10 Hz for all of Switzerland, referenced to rock sites with an estimated shear wave velocity of 1,500 m/s² in the upper 30 m. The highest hazard is found in the Wallis, in the Basel region, in Graubünden and along the Alpine front, with maximum spectral accelerations at 5 Hz frequency reaching 150 cm/s² for a return period of 475 years and 720 cm/s² for 10,000 years. Electronic supplementary material The online version of this article (doi: ) contains supplementary material, which is available to authorized users.     

近断层地震下可液化场地土反应的数值模拟

液化条件下场地土大变形是造成工程结构失效的主要原因之一。文中以某高速公路特大桥的可液化岸坡近场场地为研究对象,基于PL-Finn液化本构模型,利用FLAC 3D程序对其在地震作用下的动响应全过程进行了数值模拟分析。结果表明,PL-Finn模型可较好地反映地震过程中孔隙水压、液化区的变化规律,并能较好的预测液化后场地土的变形规律。液化引起的地基土大变形对桥梁桩基的影响需引起重视。     

Seismic waves in 3-D: from mantle asymmetries to reliable seismic hazard assessment

A global cross-section of the Earth parallel to the tectonic equator (TE) path, the great circle representing the equator of net lithosphere rotation, shows a difference in shear wave velocities between the western and eastern flanks of the three major oceanic rift basins. The low-velocity layer in the upper asthenosphere, at a depth range of 120 to 200 km, is assumed to represent the decoupling between the lithosphere and the underlying mantle. Along the TE-perturbed (TE-pert) path, a ubiquitous LVZ, about 1,000-km-wide and 100-km-thick, occurs in the asthenosphere. The existence of the TE-pert is a necessary prerequisite for the existence of a continuous global flow within the Earth. Ground-shaking scenarios were constructed using a scenario-based method for seismic hazard analysis (NDSHA), using realistic and duly validated synthetic time series, and generating a data bank of several thousands of seismograms that account for source, propagation, and site effects. Accordingly, with basic self-organized criticality concepts, NDSHA permits the integration of available information provided by the most updated seismological, geological, geophysical, and geotechnical databases for the site of interest, as well as advanced physical modeling techniques, to provide a reliable and robust background for the development of a design basis for cultural heritage and civil infrastructures. Estimates of seismic hazard obtained using the NDSHA and standard probabilistic approaches are compared for the Italian territory, and a case-study is discussed. In order to enable a reliable estimation of the ground motion response to an earthquake, three-dimensional velocity models have to be considered, resulting in a new, very efficient, analytical procedure for computing the broadband seismic wave-field in a 3-D anelastic Earth model.     

On unified analysis of uncertainty in seismic hazard assessment

（胡聿贤，陈汉尧）Ｏｎｕｎｉｆｉｅｄａｎａｌｙｓｉｓｏｆｕｎｃｅｒｔａｉｎｔｙｉｎｓｅｉｓｍｉｃｈａｚａｒｄａｓｓｅｓｓｍｅｎｔ￥Ｙｕ－ＸｉａｎＨＵａｎｄＨａｎ－ＹａｏＣＨＥＮ（ＩｎｓｔｉｔｕｔｅｏｆＧｅｏｐｈｙｓｉｃｓ，ＳｔａｔｅＳｅｉｓｍｏｌｏｇ...     

Probabilistic seismic hazard assessment of southern part of Ghana

This paper presents a seismic hazard map for the southern part of Ghana prepared using the probabilistic approach, and seismic hazard assessment results for six cities. The seismic hazard map was prepared for 10% probability of exceedance for peak ground acceleration in 50 years. The input parameters used for the computations of hazard were obtained using data from a catalogue that was compiled and homogenised to moment magnitude (Mw). The catalogue covered a period of over a century (1615–2009). The hazard assessment is based on the Poisson model for earthquake occurrence, and hence, dependent events were identified and removed from the catalogue. The following attenuation relations were adopted and used in this study—Allen (for south and eastern Australia), Silva et al. (for Central and eastern North America), Campbell and Bozorgnia (for worldwide active-shallow-crust regions) and Chiou and Youngs (for worldwide active-shallow-crust regions). Logic-tree formalism was used to account for possible uncertainties associated with the attenuation relationships. OpenQuake software package was used for the hazard calculation. The highest level of seismic hazard is found in the Accra and Tema seismic zones, with estimated peak ground acceleration close to 0.2 g. The level of the seismic hazard in the southern part of Ghana diminishes with distance away from the Accra/Tema region to a value of 0.05 g at a distance of about 140 km.     

Ancient destructive earthquakes in Chersonesus Taurica and their importance for a long-term seismic hazard assessment of the Southwestern Crimea region

Continuing the author’s earlier works, the signs of ancient destructive earthquakes in the southwest Crimea have been collected and analyzed. These signs are considered using archeoseismological method based on the data of multiannual archeological excavations at Chersonesus Taurica and agricultural hinterland around it (chora) separately on each period of total destruction, reconstruction and new construction, and on structures of different purposes and degrees of earthquake-resistance (defensive structures, large civil buildings, dwellings, fortified estates, necropolises, and burials). Indirect signs are also involved. Both direct and indirect signs unambiguously contribute to form the general scenarios of considerable destruction in the polis and its chora, with the disturbance of entire socioeconomic life, which occurred in 340(±20) and 105(±5) BC as a result of strong earthquakes with I = VIII–IX (on MSK-64 scale). Based on the results of earlier works on destructive earthquakes of a later historical period in mind, the seismic potential of the Sevastopol source zone has been estimated.