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.     

Response spectral attenuation relationships for Sumatran-subduction earthquakes and the seismic hazard implications to Singapore and Kuala Lumpur

Singapore and Kuala Lumpur, the capital of Malaysia, may well represent the classic examples of area with low seismic hazard but with high consequence. Both cities are located in a low-seismicity region of Southeast Asia, where active seismic sources are located more than 300 km away. Seismic designs have not been implemented in this seemingly low-hazard region though distant earthquakes in Sumatra had frequently shaken high-rise structures in the two cities. Several studies have been conducted to systematically assess the seismic hazards of Singapore and the Malay Peninsula. The present research particularly addresses issues in deriving a new set of attenuation relationships of peak ground acceleration (PGA), peak ground velocity (PGV) and response spectral acceleration (RSA) for the Sumatran-subduction earthquakes. To be relevant for the seismic hazard assessment of the remote metropolises, the derived attenuation relationships cover a long distance range from 150 to 1500 km. The attenuation relationships are derived using synthetic seismograms that account for source and path effects. The uncertainties in rupture parameters, such as stress drop, strike, dip and rake angles, have been defined according to the regional geological and tectonic settings as well as the ruptures of previous earthquakes. The seismic potential of the Sumatran subduction zone are high in the region from 2°N to 5°S as there has been no recurrence of great thrust events since 1861. A large event with M_w greater than 7.8 in this particular subduction zone may be capable of generating destructive ground motions in Singapore and Kuala Lumpur, even at a distance of 700 km.     

Probabilistic seismic hazard assessment for Romania and sensitivity analysis: A case of joint consideration of intermediate-depth (Vrancea) and shallow (crustal) seismicity

The earthquake risk on Romania is one of the highest in Europe, and seismic hazard for almost half of the territory of Romania is determined by the Vrancea seismic region, which is situated beneath the southern Carpathian Arc. The region is characterized by a high rate of occurrence of large earthquakes in a narrow focal volume at depth from 70 to 160 km. Besides the Vrancea area, several zones of shallow seismicity located within and outside the Romanian territory are considered as seismically dangerous. We present the results of probabilistic seismic hazard analysis, which implemented the “logic tree” approach, and which considered both the intermediate-depth and the shallow seismicity. Various available models of seismicity and ground-motion attenuation were used as the alternative variants. Seismic hazard in terms of macroseismic intensities, peak ground acceleration, and response spectra was evaluated for various return periods. Sensitivity study was performed to analyze the impact of variation of input parameters on the hazard results. The uncertainty on hazard estimates may be reduced by better understanding of parameters of the Vrancea source zone and the zones of crustal seismicity. Reduction of uncertainty associated with the ground-motion models is also very important issue for Romania.     

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.     

南北地震带南段水平向地震动衰减特征

地震动衰减关系作为抗震救灾的重要依据,一直都是现代地震学研究的重点之一.本研究使用南北地震带南段区域2009—2016年共217个地震事件获得随震中距变化的水平向地震动峰值速度(PGV)和地震动峰值加速度(PGA)经验衰减关系,并计算场地响应.研究结果显示PGV衰减关系的拟合效果较PGA更好,两者的距离衰减系数会随事件矩震级增大呈线性减小;相较于大矩震级事件,中矩震级事件在近场可能产生较衰减关系理论值更大的PGV和PGA,同时衰减关系的拟合标准差会随事件矩震级的增大而减小.进行场地响应校正后的PGV和PGA更加符合经验衰减关系,PGA的场地响应影响较PGV更强但两者的趋势一致,并与该区域前人计算得到的地壳Qs值分布对应,表明地壳介质放大或压制地震波振幅和其传递地震波能量的能力是相关联的.本文结果一定程度上揭示了南北地震带南段的地震动强度衰减特征,为未来中国西南部的抗震减灾工作提供了重要的参考.     

湖南地区分布式地震活动性模型在PSHA中的应用研究

本文采用了空间光滑地震活动性模型,该模型无需潜在震源区划分,同时发展了概率地震危险性分析新方法。根据三种地震目录资料建立了三种地震活动性模型,利用高斯光滑函数获得了湖南区域内的比值分布特征,使用了两种典型的衰减模型,计算了50年内超越概率10％的地震动峰值加速度（PGA）分布。其分析结果显示PGA分布特征与中国地震动参数区划图大体一致,部分区域PGA提高,PGA达0．05g的区域显著扩大,其中包括邵阳、湘潭、吉首、怀化等重要城市,而这种PGA分布特征与该地区地震活动性特征是一致的。概率危险性曲线的结果表明常德等地区的潜在地震危险性比湖南区域内其他城市高。表明此模型用于地震危险性计算中是简便易行的,且具有较高的精度。尤其对于地质和地震构造信息缺乏的弱震区和中强震区,该方法作为替代方法并有着广泛的应用价值。     

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

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

Probabilistic seismic hazard analysis in Nepal

The seismic ground motion hazard for Nepal has been estimated using a probabilistic approach. A catalogue of earthquakes has been compiled for Nepal and the surrounding region （latitude 26% N and 31.7% N and longitude 79° E and 90° E） from 1255 to 2011. The distribution of catalogued earthquakes, together with available geological and tectonic information were used to delineate twenty-three seismic source seismic source information and probabilistic earthquake hazard prediction relationship, peak ground accelerations （PGAs） have zones in Nepal and the surrounding region. By using the parameters in conjunction with a selected ground motion been calculated at bedrock level with 63%, 10%, and 2% probability of exceedance in 50 years. The estimated PGA values are in the range of 0.07-0.16 g, 0.21 0.62 g, and 0.38-1.1 g for 63%, 10%, and 2% probability of exceedance in 50 years, respectively. The resulting ground motion maps show different characteristics of PGA distribution, i.e., high hazard in the far-western and eastern sections, and low hazard in southern Nepal. The quantified PGA values at bedrock level provide information for microzonation studies in different parts of the country.     

A new seismic hazard assessment in the region of Andalusia (Southern Spain)

A probabilistic seismic hazard assessment of Andalusia (Southern Spain) 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 was here calculated in terms of magnitude, using published spectral ground-motion models. Moreover, we considered different ground-motion models for the Atlantic sources, since the attenuation of those motions seems to be slower, as evidenced in the case of the extensive macroseismic areas of earthquakes like those occurred in the years 1755, 1969 and 2007. A comprehensive review of the seismic catalogue and of the seismogenic models proposed for the region was carried out, including those for Northern Africa, which is part of the influence area. Hazard calculations were performed following the Probabilistic Seismic Hazard Assessment (PSHA) methodology using a logic tree, which accounts for six different seismic source zonings and five different ground-motion attenuation relationships. Hazard maps in terms of PGA and SA (0.2 s) and SA (1 s) and coefficient of variation (COV) maps, for the 475-year return period were first obtained in rock sites. A geotechnical classification and amplification factors were proposed and new hazard maps including local effects were represented, showing PGA values ranging from 24 to 370 cm/s² for the whole Andalusian territory, with the highest expected values (PGA > 300 cm/s²) in some parts of the Granada Province and in the town of Vélez Málaga. Lowest values (PGA < 50 cm/s²) correspond to some towns of the Huelva and Córdoba provinces. The inclusion of soil effects provides a more detailed picture of the actual hazard the region is subjected to.     

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

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

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.     

Discussion on uncertainties，attenuation of ground motion and aseismic design criterion

Ｄｉｓｃｕｓｓｉｏｎｏｎｕｎｃｅｒｔａｉｎｔｉｅｓ，ａｔｔｅｎｕａｔｉｏｎｏｆｇｒｏｕｎｄ　ｍｏｔｉｏｎａｎｄａｓｅｉｓｍｉｃｄｅｓｉｇｎｃｒｉｔｅｒｉｏｎＴｉａｎ－ＺｈｏｎｇＺＨＡＮＧ（张天中）；Ｙｕｎ－ＳｈｅｎｇＭＡ（马云生）ａｎｄＸｉＳＨＵ（舒...     

Copula joint function and its application in probability seismic hazard analysis

A mature mathematical technique called copula joint function is introduced in this paper, which is commonly used in the financial risk analysis to estimate uncertainty. The joint function is generalized to the n-dimensional Frank’s copula. In addition, we adopt two attenuation models proposed by YU and Boore et al, respectively, and construct a two-dimensional copula joint probabilistic function as an example to illustrate the uncertainty treatment at low probability. The results show that copula joint func...     

地震岩相识别概率表征方法

储层岩相分布信息是油藏表征的重要参数,基于地震资料开展储层岩相识别通常具有较强的不确定性.传统方法仅获取唯一确定的岩相分布信息,无法解析反演结果的不确定性,增加了油藏评价的风险.本文引入基于概率统计的多步骤反演方法开展地震岩相识别,通过在其各个环节建立输入与输出参量的统计关系,然后融合各环节概率统计信息构建地震数据与储层岩相的条件概率关系以反演岩相分布概率信息.与传统方法相比,文中方法通过概率统计关系表征了地震岩相识别各个环节中地球物理响应关系的不确定性,并通过融合各环节概率信息实现了不确定性传递的数值模拟,最终反演的岩相概率信息能够客观准确地反映地震岩相识别结果的不确定性,为油藏评价及储层建模提供了重要参考信息.模型数据和实际资料应用验证了方法的有效性.     

Simulations of Seismic Hazard for the Pacific Northwest of the United States from Earthquakes Associated with the Cascadia Subduction Zone

-- We investigate the impact of different rupture and attenuation models for the Cascadia subduction zone by simulating seismic hazard models for the Pacific Northwest of the U.S. at 2% probability of exceedance in 50 years. We calculate the sensitivity of hazard (probabilistic ground motions) to the source parameters and the attenuation relations for both intraslab and interface earthquakes and present these in the framework of the standard USGS hazard model that includes crustal earthquakes. Our results indicate that allowing the deep intraslab earthquakes to occur anywhere along the subduction zone increases the peak ground acceleration hazard near Portland, Oregon by about 20%. Alternative attenuation relations for deep earthquakes can result in ground motions that differ by a factor of two. The hazard uncertainty for the plate interface and intraslab earthquakes is analyzed through a Monte-Carlo logic tree approach and indicates a seismic hazard exceeding 1 g (0.2 s spectral acceleration) consistent with the U.S. National Seismic Hazard Maps in western Washington, Oregon, and California and an overall coefficient of variation that ranges from 0.1 to 0.4. Sensitivity studies indicate that the paleoseismic chronology and the magnitude of great plate interface earthquakes contribute significantly to the hazard uncertainty estimates for this region. Paleoseismic data indicate that the mean earthquake recurrence interval for great earthquakes is about 500 years and that it has been 300 years since the last great earthquake. We calculate the probability of such a great earthquake along the Cascadia plate interface to be about 14% when considering a time-dependent model and about 10% when considering a time-independent Poisson model during the next 50-year interval.     

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 and uncertainty analysis of the Taiwan area

The objectives of this paper are (1) to obtain estimates on the effect of uncertainties of the hazard model, and (2) to evaluate the seismic hazard in Taiwan for structural analysis and design purposes. The seismic hazard in the Taiwan area is presented in terms of an iso-acceleration map. Such a map is developed for return periods of peak ground acceleration of 225 years and 475 years. The contour map of b-values and mean occurence rates for this region is also presented. Uncertainty analyses of model parameters in hazard analysis are concentrated on the analysis of dispersion of PGA values and the probabilistic modeling of stationary and nonstationary Poisson models of occurrences. Th e overall results are considered to be conservative since for most uncertainty analyses the more conservative values are used.     

Evaluation of the seismic hazard of Lebanon

This paper presents the results of a study undertaken todetermine the seismic hazard of Lebanon. The seismic hazard evaluation wasconducted using probabilistic methods of hazard analysis. Potential sourcesof seismic activities that affect Lebanon were identified and the earthquakerecurrence relationships of these sources were developed from instrumentalseismology data, historical records, and earlier studies undertaken toevaluate the seismic hazard of neighboring countries. The sensitivityof the results to different assumptions regarding the seismic sources in theLebanese segment and choice of the attenuation relationship wasevaluated. Maps of peak ground acceleration contours, based on 10percent of probability of exceedance in 50 years and 100 years time spans,were developed.     

基于性态的液化场地多跨桥梁桩基地震响应分析

针对液化场地多跨简支桩基桥梁体系,考虑地震随机性的不确定性和认知的不确定性,结合地震危险性曲线自身的不确定性,推导性态指标危险性曲线的解析表达式.利用地震动强度指标PGV/PGA,输入不同幅值的地震动,进行液化场地多跨桩基桥梁体系地震反应有限元分析.基于有限元数值分析结果,选取地震过程中关键位置位移和弯矩的最大值作为性...     

Earthquake data review and seismic hazard analysis for the Hong Kong region

This paper presents the results of a review of earthquake data and a study of seismic hazard for Hong Kong. A region of about 660 km × 660 km around Hong Kong was selected for the study. In this study, the earthquake information available for the region was reviewed, and where possible, earthquake magnitudes were calculated. Since there is no strong motion record for any earthquake within the region, Joyner and Boore's attenuation law¹ has been used for the analysis. The results show that the seismic hazard in Hong Kong is relatively small, but because of the uncertainty inherent in the assumed attenuation relationship, accurate prediction of peak ground acceleration is not possible.