Intensity attenuation in the Campania region,Southern Italy

Ground motion prediction equations (GMPE) in terms of macroseismic intensity are a prerequisite for intensity-based shake maps and seismic hazard assessment and have the advantage of direct relation to earthquake damage and good data availability also for historical events. In this study, we derive GMPE for macroseismic intensity for the Campania region in southern Italy. This region is highly exposed to the seismic hazard related to the high seismicity with moderate- to large-magnitude earthquakes in the Appenninic belt. The relations are based on physical considerations and are easy to implement for the user. The uncertainties in earthquake source parameters are accounted for through a Monte Carlo approach and results are compared to those obtained through a standard regression scheme. One relation takes into account the finite dimensions of the fault plane and describes the site intensity as a function of Joyner–Boore distance. Additionally, a relation describing the intensity as a function of epicentral distance is derived for implementation in cases where the dimensions of the fault plane are unknown. The relations are based on an extensive dataset of macroseismic intensities for large earthquakes in the Campania region and are valid in the magnitude range M _w = 6.3–7.0 for shallow crustal earthquakes. Results indicate that the uncertainties in earthquake source parameters are negligible in comparison to the spread in the intensity data. The GMPE provide a good overall fit to historical earthquakes in the region and can provide the intensities for a future earthquake within 1 intensity unit.     

Convex set theory-based seismic hazard analysis of low seismicity area

Historical seismic data and seismogenic information are quite scarce for the low seismicity region, and modeling the parameters uncertainties based on probabilistic model is suspicious. The convex set theory-based seismic hazard analysis approach is proposed. The uncertainties of b value, the annual occurrence rate v and the upper bound magnitude M_u are described by the envelop bound convex model and the ellipsoidal bound convex model. Convex analysis method and China probabilistic seismic hazard analysis methodology are combined to perform a bound seismic hazard analysis for Ningbo city, China. The seismic intensity interval obtained using the bound seismic hazard analysis is compared with that calculated using China probabilistic seismic hazard analysis methodology. The sensitivity analysis indicates that the interval of seismic intensity is most sensitive to the annual occurrence rate v. Furthermore, the different convex models have little effect on the interval of seismic intensity.     

Generation of hazard-consistent ground motion

This paper presents an integrated approach for evaluating seismic hazard and establishing ground motion at a site. In this approach, we combine the advantage of probabilistic and deterministic seismic hazard analyses and generate synthetic ground motion by considering the characteristics of seismic source, path attenuation, and local soil condition. Furthermore, uncertainties in seismic and soil parameters are taken into account. The proposed approach can be used to establish site-specific ground motion for engineering applications.     

Building damage scenarios based on exploitation of Housner intensity derived from finite faults ground motion simulations

In this paper earthquake damage scenarios for residential buildings (about 4200 units) in Potenza (Southern Italy) have been estimated adopting a novel probabilistic approach that involves complex source models, site effects, building vulnerability assessment and damage estimation through Damage Probability Matrices. Several causative faults of single seismic events, with magnitude up to 7, are known to be close to the town. A seismic hazard approach based on finite faults ground motion simulation techniques has been used to identify the sources producing the maximum expected ground motion at Potenza and to generate a set of ground motion time histories to be adopted for building damage scenarios. Additionally, site effects, evaluated in a previous work through amplification factors of Housner intensity, have been combined with the bedrock values provided by hazard assessment. Furthermore, a new relationship between Housner and EMS-98 macroseismic intensity has been developed. This relationship has been used to convert the probability mass functions of Housner intensity obtained from synthetic seismograms amplified by the site effects coefficients into probability mass function of EMS-98 intensity. Finally, the Damage Probability Matrices have been applied to estimate the damage levels of the residential buildings located in the urban area of Potenza. The proposed methodology returns the full probabilistic distribution of expected damage, thus avoiding average damage index or uncertainties expressed in term of dispersion indexes.     

潜在地震滑坡危险区区划方法

不同地区地震活动的强度和频率是不同的.基于地震危险性分析的地震滑坡危险研究在综合了地震烈度、位置、复发时间等因素的基础上,考虑了地震动峰值加速度时空分布的特点,可以有效地应用于潜在地震滑坡危险区区划.以汶川地震灾区为研究对象,根据研究区的地质构造、地震活动特点等划分出灾区的潜在震源区,对该区进行地震危险性分析,并在此基础上采用综合指标法做出基于地震危险性分析的地震滑坡危险性区划.所得地震滑坡危险性区划按照滑坡危险程度分为高危险、较高危险、较低危险和低危险四级,表示未来一段时间内研究区在遭受一定超越概率水平的地震动作用下,不同地区地震滑坡发生的可能程度. 本文给出的地震滑坡危险性区划结果中,汶川地震滑坡崩塌较发育的汶川、北川、茂县等部分区域均处于高危险或较高危险区域;在对具有较高DEM精度的北川擂鼓镇地区所作的地震滑坡危险性区划中,汶川地震中实际发生的地震滑坡灾害与地震滑坡危险区划结果表现出较好的一致性.对区域范围而言,基于地震危险性分析的地震滑坡区划,可为初期阶段的土地规划使用及重大工程选址提供参考.     

Incremental dynamic analysis of wood-frame houses in Canada: Effects of dominant earthquake scenarios on seismic fragility

Seismic fragility can be assessed by conducting incremental dynamic analysis (IDA). This study extends the current conditional mean spectrum (CMS)-based record selection approach for IDA by taking into account detailed seismic hazard information. The proposed method is applied to conventional wood-frame houses in Canada, across which dominant earthquake scenarios and associated hazard levels vary significantly. Effects due to different seismic environments, site conditions, CMS-based record selection methods, and house models are investigated by comparing various seismic fragility models. Moreover, relative impact of the key characteristics is evaluated in terms of seismic loss curve for a group of wood-frame houses. Importantly, a close examination of regional seismic hazard characteristics using seismic hazard curve and seismic deaggregation facilitates the deeper understanding of the impact of ground motion characteristics on seismic fragility. A comprehensive and systematic assessment of key uncertainties associated with seismic fragility is provided.     

From ambient noise recordings to site effect assessment: The case study of Barcelona microzonation

As local site effects have a drastic influence on seismic hazard, it is a major issue to characterize them in vulnerable areas such as highly urbanized zones, like Barcelona city. The aim of this work is to improve the knowledge of geophysical characteristics of Barcelona in the perspective of a seismic microzonation that takes into account site effect. The first step was to gather the existing data from geological, geotechnical, geophysical, and seismological investigations, bringing us to keep the four zones proposed by previous work as a base of zonation. The second step was to characterize each zone by time-averaged shear-wave velocity and fundamental resonance frequency, with ambient noise techniques over 17 sites, providing new knowledge about the soil of Barcelona. The third step was to propose an amplification function between an average soil for each zone and a standard reference rock site, using empirically based propositions and to compare them to previous numerical approaches.     

Accounting for source location errors in the bayesian analysis of seismicity and seismic hazard

A method is presented for incorporating the uncertainties associated with hypocentral locations in the formulation of probabilistic models of the time and space distributions of the activity of potential seismic sources, as well as of the resulting seismic hazard functions at sites in their vicinity. For this purpose, a bayesian framework of analysis is adopted, where the probabilistic models considered are assumed to have known forms and uncertain parameters, the distribution of the latter being the result of an a priori assessment and its updating through the incorporation of the direct statistical information, including the uncertainty associated with the relations between the actual hypocentral locations and the reported data. This uncertainty is incorporated in the evaluation of the likelihood function of the parameters to be estimated for a given sample of recorded locations. For the purpose of illustration, the method proposed is applied to the modelling of the seismic sources near a site close to the southern coast of Mexico. The results of two alternate algorithms for the incorporation of location uncertainties are compared with those arising from neglecting those uncertainties. One of them makes use of Monte Carlo simulation, while the other is based on a closed-form analytical integration following the introduction of some simplifying assumptions. For the particular case studied, accounting for location uncertainties gives place to significant changes in the probabilistic models of the seismic sources. Deviations of the same order of magnitude can be ascribed to differences in the mathematical and/or numerical tools used in the uncertainty analysis. The resulting variability of the seismic hazard at the site of interest is less pronounced than that affecting the estimates of activity of individual seismic sources.     

"Parametric-historic" Procedure for Probabilistic Seismic Hazard Analysis Part II: Assessment of Seismic Hazard at Specified Site

—A new methodology for probabilistic seismic hazard analysis is described. The approach combines the best features of the "deductive" (Cornell, 1968) and "historic" (Veneziano et al., 1984) procedures. It can be called a "parametric-historic" procedure. The maximum regional magnitude mmax is of paramount importance in this approach and Part I of the authors’ work (Kijko and Graham, 1998) was dedicated to developing efficient statistical procedures that can be used for the evaluation of this parameter. In Part II the approach of a probabilistic seismic hazard assessment at a given site is described. The approach permits the utilization of incomplete earthquake catalogues. It is assumed that a typical catalogue contains two types of information historical macroseismic events that occurred over a period of a few hundred years and recent, instrumental data. The historical part of the catalogue contains only the strongest events, whereas the complete part can be divided into several subcatalogues, each assumed complete above a specified threshold of magnitude. The author’s approach also takes into account uncertainty in the determination of the earthquake magnitude. The technique has been developed specifically for the estimation of seismic hazard at individual sites, without the subjective judgment involved in the definition of seismic source zones, in which specific active faults have not been mapped and identified, and where the causes of seismicity are not well understood. As an example of the application of the new technique, the results of a typical hazard analysis for a hypothetical engineering structure located in the territory of South Africa are presented. It was assumed that the only reliable information in the assessment of the seismic hazard parameters in the vicinity of the selected site comes from a knowledge of past seismicity. The procedure was applied to seismic data that were divided into an incomplete part, containing only the largest events, and two complete parts, containing information obtained from instruments. The simulation experiments described in Part I of our study have shown that the Bayesian estimator K-S-B tends to perform very well, especially in the presence of inevitable deviations from the simple Gutenberg–Richter model. In the light of this fact value &gif1; = 6.66 - 0.44, which was obtained from the K-S-B technique, was regarded as the best choice. At an exceedance probability of 10^х per annum, the median value of peak ground acceleration on rock at the site is 0.31g, and at an exceedance probability of 10^ц per annum, the median peak ground acceleration at the site is 0.39g. The median value of the maximum possible acceleration at the site is 0.40g, which was calculated from attenuation formulae by assuming the occurrence of the strongest possible earthquake, e.g., with magnitude &gif1; = 6.66 at distance 10 km.     

A Bayesian Approach to Seismic Hazard Estmation: Maximum Values of Magnitudes and Peak Ground Accelerations

Earthquake Researeh in Ch一na461 .METHODLet R be some value measured or estimated as a sequenee ina’‘Past”time interval(一丁,O)(I)万‘月,=(RI,…,R。),R,之R。,R=nlaX l二f匕11(RI,…,R,,) Values(l)eould have an arbitrary Physieal nature.BelowweshalleonsiderEq.(l)asearthquakemagnitudes in a given seismic aetive region or logarithms of seismie Peak ground aeeelerations at习given site.Ro isa而nimum eutoff value;it 15 defined by Possibilities of registration systems or wasehosen as the …     

青藏高原东北缘地震活动性广义帕累托模型的全域敏感性分析

基于广义帕累托分布构建地震活动性模型,因其输入参数取值难以避免不确定性,导致依据该模型所得的地震危险性估计结果具有不确定性。鉴于此,本文选取青藏高原东北缘为研究区,提出了基于全域敏感性分析的地震危险性估计的不确定性分析流程和方法。首先,利用地震活动性广义帕累托模型,进行研究区地震危险性估计;然后,选取地震记录的起始时间和震级阈值作为地震活动性模型的输入参数,采用具有全域敏感性分析功能的E-FAST方法,对上述两个参数的不确定性以及两参数之间的相互作用对地震危险性估计不确定性的影响进行定量分析。结果表明：地震危险性估计结果（不同重现期的震级重现水平、震级上限及相应的置信区间）对两个输入参数中的震级阈值更为敏感;不同重现期的地震危险性估计结果对震级阈值的敏感程度不同;对不同的重现期而言,在影响地震危险性估计结果的不确定性上,两个输入参数之间存在非线性效应,且非线性效应程度不同。本文提出的不确定性分析流程和方法,可以推广应用于基于其它类型地震活动性模型的地震危险性估计不确定性分析。     

Comment on Sousa, M. L. and Costa, A. C., ��Ground motion scenarios consistent with probabilistic seismic hazard disaggregation analysis. Application to Mainland Portugal��

Disaggregation of the seismic hazard has become a popular technique to convey information on the main sources contributing to the hazard at a particular site. Recently published work adopts geographic disaggregation analysis of seismic hazard as a tool to identify dominant rupture scenarios for municipalities in Mainland Portugal. The authors conclude that the seismic hazard in South and Central Portugal is dominated by the seismicity that takes place offshore, around 70km WSW of Cape S. Vicente, both for the return periods of 475 years and 975 years. Whilst recognizing the merits of the approach taken and the utility of the tools developed, we take issue with this last conclusion. We consider that the proposed disaggregation returns a picture of the biases in the hazard analysis, more than any real feature of the distribution of relevant seismogenic sources.     

A generalized conditional intensity measure approach and holistic ground‐motion selection

A generalized conditional intensity measure (GCIM) approach is proposed for use in the holistic selection of ground motions for any form of seismic response analysis. The essence of the method is the construction of the multivariate distribution of any set of ground‐motion intensity measures conditioned on the occurrence of a specific ground‐motion intensity measure (commonly obtained from probabilistic seismic hazard analysis). The approach therefore allows any number of ground‐motion intensity measures identified as important in a particular seismic response problem to be considered. A holistic method of ground‐motion selection is also proposed based on the statistical comparison, for each intensity measure, of the empirical distribution of the ground‐motion suite with the ‘target’ GCIM distribution. A simple procedure to estimate the magnitude of potential bias in the results of seismic response analyses when the ground‐motion suite does not conform to the GCIM distribution is also demonstrated. The combination of these three features of the approach make it entirely holistic in that: any level of complexity in ground‐motion selection for any seismic response analysis can be exercised; users explicitly understand the simplifications made in the selected suite of ground motions; and an approximate estimate of any bias associated with such simplifications is obtained. Copyright © 2010 John Wiley & Sons, Ltd.     

Validation through HVSR measurements of a method for the quick detection of site amplification effects from intensity data: an application to a seismic area in Northern Italy

An instrumental validation is attempted of an innovative approach devoted to the quick individuation, from macroseismic data, of site amplification phenomena able to significantly modify seismic hazard levels expected on the basis of average propagation effects only. According to this methodology, two evaluations of hazard are performed at each investigated locality: the former, obtained by epicentral intensity data ‘reduced’ at the site through a probabilistic attenuation function and, the latter, computed by integrating such data with seismic effects actually observed at the site during past earthquakes. The comparison, for each locality, between these two hazard estimates allow to orientate the identification of those sites where local amplifications of earthquake ground motion could be significant. In order to check such methodology, indications obtained in this way from macroseismic data are compared with the estimates of transfer functions performed through the HVSR technique applied to microtremors. Results concerning municipalities located in a seismic area of Northern Italy indicate a good agreement between macroseismic and instrumental estimates.     

A method to estimate intensity occurrence probabilities in low seismic activity regions

The estimation of site intensity occurrence probabilities in low seismic activity regions has been studied from different points of view. However, no method has been definitively established because several problems arise when macroseismic historical data are incomplete and the active zones are not well determined. The purpose of this paper is to present a method that estimates site occurrence probabilities and at the same time measures the uncertainties inherent in these probabilities in low activity regions. The region to be studied is divided into very broad seismic zones. An exponential intensity probability law is adjusted for each zone and the degree of uncertainty in the assumed incompleteness of the catalogue is evaluated for each intensity. These probabilities are used to establish what may be termed ‘prior site occurrence models’. A Bayesian method is used to improve ‘prior models’ and to obtain the ‘posterior site occurrence models’. Epicentre locations are used to recover spatial information lost in the prior broad zoning. This Bayesian correction permits the use of specific attenuation for different events and may take into account, by means of conservative criteria, epicentre location errors. Following Bayesian methods, probabilities are assumed to be random variables and their distribution may be used to estimate the degree of uncertainty arising from (a) the statistical variance of estimators, (b) catalogue incompleteness and (c) mismatch of data to prior assumptions such as Poisson distribution for events and exponential distribution for intensities. The results are maps of probability and uncertainty for each intensity. These maps exhibit better spatial definition than those obtained by means of simple, broad zones. Some results for Catalonia (NE of Iberian Peninsula) are shown.     

A GIS-based time-dependent seismic source modeling of Northern Iran

The first step in any seismic hazard study is the definition of seismogenic sources and the estimation of magnitude-frequency relationships for each source. There is as yet no standard methodology for source modeling and many researchers have worked on this topic. This study is an effort to define linear and area seismic sources for Northern Iran. The linear or fault sources are developed based on tectonic features and characteristic earthquakes while the area sources are developed based on spatial distribution of small to moderate earthquakes. Time-dependent recurrence relationships are developed for fault sources using renewal approach while time-independent frequency-magnitude relationships are proposed for area sources based on Poisson process. GIS functionalities are used in this study to introduce and incorporate spatial-temporal and geostatistical indices in delineating area seismic sources. The proposed methodology is used to model seismic sources for an area of about 500 by 400 square kilometers around Tehran. Previous researches and reports are studied to compile an earthquake/fault catalog that is as complete as possible. All events are transformed to uniform magnitude scale; duplicate events and dependent shocks are removed. Completeness and time distribution of the compiled catalog is taken into account. The proposed area and linear seismic sources in conjunction with defined recurrence relationships can be used to develop time-dependent probabilistic seismic hazard analysis of Northern Iran.     

Seismic hazard analysis for engineering sites based on the stochastic finite-fault method

Seismic hazard analyses are mainly performed using either deterministic or probabilistic methods. However, there are still some defects in these statistical model-based approaches for regional seismic risk assessment affected by the near-field of large earthquakes. Therefore, we established a deterministic seismic hazard analysis method that can characterize the entire process of ground motion propagation based on stochastic finite-fault simulation, and we chose the site of the Xiluodu dam to demonstrate the method. This method can characterize earthquake source properties more realistically than other methods and consider factors such as the path and site attenuation of seismic waves. It also has high computational efficiency and is convenient for engineering applications. We first analyzed the complexity of seismogenic structures in the Xiluodu dam site area, and then an evaluation system for ground motion parameters that considers various uncertainties is constructed based on a stochastic finite-fault simulation. Finally, we assessed the seismic hazard of the dam site area comprehensively. The proposed method was able to take into account the complexity of the seismogenic structures affecting the dam site and provide multi-level parameter evaluation results corresponding to different risk levels. These results can be used to construct a dam safety assessment system of an earthquake in advance that provides technical support for rapidly and accurately assessing the post-earthquake damage state of a dam, thus determining the influence of an earthquake on dam safety and mitigating the risk of potential secondary disasters.     

Macroseismic intensity attenuation in Iran

Macroseismic intensity data plays an important role in the process of seismic hazard analysis as well in developing of reliable earthquake loss models. This paper presents a physical-based model to predict macroseismic intensity attenuation based on 560 intensity data obtained in Iran in the time period 1975–2013. The geometric spreading and energy absorption of seismic waves have been considered in the proposed model. The proposed easy to implement relation describes the intensity simply as a function of moment magnitude, source to site distance and focal depth. The prediction capability of the proposed model is assessed by means of residuals analysis. Prediction results have been compared with those of other intensity prediction models for Italy, Turkey, Iran and central Asia. The results indicate the higher attenuation rate for the study area in distances less than 70 km.     

Global sensitivity analysis for stochastic ground motion modeling in seismic-risk assessment

Seismic risk assessment requires adoption of appropriate models for the earthquake hazard, the structural system and for its performance, and quantification of the uncertainties involved in these models through appropriate probability distributions. Characterization of the seismic hazard comprises undoubtedly the most critical component of this process, the one associated with the largest amount of uncertainty. For applications involving dynamic analysis this hazard is frequently characterized through stochastic ground motion models. This paper discusses a novel, global sensitivity analysis for the seismic risk with emphasis on such a stochastic ground motion modeling. This analysis aims to identify the overall (i.e. global) importance of each of the uncertain model parameters, or of groups of them, towards the total risk. The methodology is based on definition of an auxiliary density (distribution) function, proportional to the integrand of the integral quantifying seismic risk, and on comparison of this density to the initial probability distribution for the model parameters of interest. Uncertainty in the rest of the model parameters is explicitly addressed through integration of their joint auxiliary distribution to calculate the corresponding marginal distributions. The relative information entropy is used to quantify the difference between the compared density functions and an efficient approach based on stochastic sampling is introduced for estimating this entropy for all quantities of interest. The framework is illustrated in an example that adopts a source-based stochastic ground motion model, and valuable insight is provided for its implementation within structural engineering applications.     

An Anisotropic Attenuation Law of Macroseismic Intensity Performed on Virtual Intensity Distribution of Seismogenic Zones

An anisotropic attenuation law, based on an anisotropic characterization of intensity distribution for seismogenic zones, is proposed. This approach, that distinguishes itself for its consistency to the observed data, initially reconfigured by filtering procedures, is particularly suitable for seismic hazard evaluation.