PSHA: is it science?

Probabilistic seismic hazard analysis (PSHA) is beginning to be seen as unreliable. The problem with PSHA is that its data are inadequate and its logic is defective. Much more reliable, and more scientific, are deterministic procedures, especially when coupled with engineering judgment.     

Seismic Hazard Analysis — Quo vadis?

The paper is dedicated to the review of methods of seismic hazard analysis currently in use, analyzing the strengths and weaknesses of different approaches. The review is performed from the perspective of a user of the results of seismic hazard analysis for different applications such as the design of critical and general (non-critical) civil infrastructures, technical and financial risk analysis. A set of criteria is developed for and applied to an objective assessment of the capabilities of different analysis methods. It is demonstrated that traditional probabilistic seismic hazard analysis (PSHA) methods have significant deficiencies, thus limiting their practical applications. These deficiencies have their roots in the use of inadequate probabilistic models and insufficient understanding of modern concepts of risk analysis, as have been revealed in some recent large scale studies. These deficiencies result in the lack of ability of a correct treatment of dependencies between physical parameters and finally, in an incorrect treatment of uncertainties. As a consequence, results of PSHA studies have been found to be unrealistic in comparison with empirical information from the real world. The attempt to compensate these problems by a systematic use of expert elicitation has, so far, not resulted in any improvement of the situation. It is also shown that scenario-earthquakes developed by disaggregation from the results of a traditional PSHA may not be conservative with respect to energy conservation and should not be used for the design of critical infrastructures without validation. Because the assessment of technical as well as of financial risks associated with potential damages of earthquakes need a risk analysis, current method is based on a probabilistic approach with its unsolved deficiencies.

Traditional deterministic or scenario-based seismic hazard analysis methods provide a reliable and in general robust design basis for applications such as the design of critical infrastructures, especially with systematic sensitivity analyses based on validated phenomenological models. Deterministic seismic hazard analysis incorporates uncertainties in the safety factors. These factors are derived from experience as well as from expert judgment. Deterministic methods associated with high safety factors may lead to too conservative results, especially if applied for generally short-lived civil structures. Scenarios used in deterministic seismic hazard analysis have a clear physical basis. They are related to seismic sources discovered by geological, geomorphologic, geodetic and seismological investigations or derived from historical references. Scenario-based methods can be expanded for risk analysis applications with an extended data analysis providing the frequency of seismic events. Such an extension provides a better informed risk model that is suitable for risk-informed decision making.     

On the long-term seismic hazard analysis in the Zhangjiakou–Penglai seismotectonic zone, China

The Zhangjiakou–Penglai seismotectonic zone (ZPSZ) lies in the northern part of North China and extends along the Zhangjiakou–Beijing–Tianjin–Bohai Bay–Penglai–Yellow Sea. It is about 900 km long and some 250 km wide in a northwest direction. The great Sanhe-Pinggu (M_S=8.0) earthquake occurred on September 1679 and the Tangshan (M_S=7.8) earthquake on July 1976 caused serious economic and life losses. According to some differences in crust structure and regional tectonic stress field, the ZPSZ is divided into western and eastern segment by the 117°E line for study on long-term seismic hazard analysis. An analysis of Gutenberg–Richter's empirical relation of earthquake-frequency and time process of historic and recent earthquakes along the eastern and western segments shows that the earthquake activity obeys a Poisson process, and these calculations indicate that the earthquake occurrence probability of M_S=6.0–6.9 is 0.77–0.83 in the eastern segment and the earthquake occurrence probability of M_S=7.0–7.9 is 0.78–0.80 in the western segment of the ZPSZ during a period from 2005 to 2015.     

Seismic hazard maps for the U.K.

Past studies of seismic hazard in the U.K. that have used modern probabilistic methods of hazard assessment have been site-specific studies, mostly in connection with nuclear installations. There has been a need for general-purpose maps of seismic hazard to show relative variation of exposure within the U.K. and to give some guidance on absolute values. Such maps have now been produced, incorporating, for the first time, the wealth of new information on historical earthquakes in Britain that has been gathered over the last 15 years. The hazard calculations were undertaken using a new computer code based on the USGS program SEISRISK III, but incorporating a logic tree approach to model variation in the input parameters (e.g. focal depth) or uncertainty in the formulation of the model (e.g. attenuation parameters). An innovative approach was taken to the formulation of seismic source zones, in which two overlapping models were employed. The first of these uses relatively broad source zones based loosely on an interpretation of seismicity and tectonics, while the second uses numerous small zones that reflect the locations of past significant earthquakes. This double approach (using the logic tree methodology) has the merit of both considering the general trend of earthquake activity as well as focusing in on known danger spots. The results show that the areas of highest hazard are western Scotland, north-western England and Wales, where the intensity with 90% probability of non-exceedance in 50 years is 6 EMS.     

地震迁移的类型、特征及机制讨论

地震活动的迁移是指地震沿着某一方向循序地发生,是地震活动总体无序中局部出现的有序结构。通过系统分析区域上典型的地震迁移现象可以发现,迁移可划分为沿断裂走向的纵向迁移、垂直断裂走向的横向迁移与岩石圈尺度的深源迁移三大类。结合具体的实例分析,可初步归纳出纵向迁移（包括单向、双向、反复和跳跃式迁移等常见形式）、横向迁移和深源地震迁移的主要特征,并初步估算出了不同类型迁移的速度值,其中沿全球板块边界纵向迁移平均速度约为V=569 km/a,沿亚板块边界的平均速度约为V=120 km/a,沿大陆内部断裂带平均速度约为V=50 km/a。横向迁移相对比较复杂,初步发现在东亚存在两种速度分别为约20 km/a、50 km/a的迁移现象。而深源地震迁移速度的全球平均值大约为360 km/a。地震的纵向与横向迁移都存在不同层次和级别,也存在多种不同频率、能量与速度的形变波与迁移现象,这很可能是区域上大地震丛集活动过程中断层相互作用、地震应力触发和岩石圈尺度的形变波传播等因素共同作用的结果,而这种大空间尺度上的地震迁移现象的存在及其所具有的规律性特征,显然可为开展区域地震危险性分析提供新的思路和方法参考。      

Seismic Hazard Studies in the U.K.: Source Specification Problems of Intraplate Seismicity

Although the U.K. is in an area of only low to moderate seismicity, the seismic hazard is sufficient to pose a threat to sensitive structures such as chemical plants and nuclear facilities. In quantifying the level of hazard by conventional probabilistic methodology, however, some problems arise in attempting to interpret earthquake data in terms of geological structure and faults. In the U.K., not only is it impossible to identify any demonstrably active faults, but also it is extremely difficult to discern any relationship between the pattern of seismicity and local or regional geological structure.This study discusses the use of two zonation approaches which complement each other in such a way that the general character and trend of seismicity is preserved. In one approach, the zonation is informed by the structural geology, where possible; geological zonation is avoided if it produces sources with heterogeneous seismicity. In the other approach, the record of past earthquakes is divided up into very small zones around individual epicentres or groups of epicentres, the size of each zone usually being proportional to the uncertainty in the epicentral determination of the appropriate event. This zonation preserves an observed tendency of some British earthquakes to repeat themselves. It is suggested that, in intraplate areas such as the U.K., it is often inappropriate to attempt to model individual fault sources. No faults in the U.K. are provably active. Because an earthquake of moderate size can occur on a very short fault segment, it is impractical to restrict fault modelling to major features. Even the two largest U.K. faults, suspected to be active, pose problems in attributing historical seismicity to them as distinct features.     

Development of spectral hazard maps for a proposed revision of the Indonesian Seismic Building Code

The need to revise the current Indonesian Seismic Hazard Map contained in Indonesian Earthquake Resistant Building Code SNI 03-1726-2002 which partially adopts the concept of UBC 1997, was driven among others by the desire to better reflect the potential larger earthquake disasters faced by the nation in the future. The much larger than maximum predicted Aceh Earthquake (M _w 9.0–9.3) of 2004, followed by the destruction observed during the 2005 Nias Earthquake (M _w 8.7) urgently underline to need to consider the new conceptual approach and technological shift shown in the transition of UBC 1997 to IBC 2006. This paper presents research works for developing spectral hazard maps for Indonesia. Some improvements in seismic hazard analysis were implemented using recent seismic records. Seismic sources were modeled by background, fault, and subduction zones by considering a truncated exponential model, pure characteristic model or both models. A logic tree method was performed to account for the epistemic uncertainty and several attenuation functions were selected. Maps of PGA and spectral accelerations for a short period (0.2 s) and for a 1-s period were then developed using a probabilistic approach. The maps will be proposed as a revision for the current seismic hazard map in the Indonesian Seismic Building Code.     

Seismic hazard evaluation for dams in northern Colorado, U.S.A.

A seismic hazard evaluation for three dams in the Rocky Mountains of northern Colorado is based on a study of the historical seismicity. To model earthquake occurrence as a random process utilizing a maximum likelihood method, the catalog must exhibit random space-time characteristics. This was achieved using a declustering procedure and correction for completeness of recording. On the basis of the resulting a- and b-values, probabilistic epicentral distances for a 2 × 10^–5 annual probability were calculated. For a random earthquake of magnitude M _L 6.0–6.5, this distance is 15 km. Suggested ground motion parameters were estimated using a probabilistic seismic hazard analysis. Critical peak horizontal accelerations at the dams are 0.22g if median values are assumed and 0.39g if variable attenuation and seismicity rates are taken into account. For structural analysis of the dams, synthetic acceleration time series were calculated to match the empirical response spectra. In addition, existing horizontal strong motion records from two Mammoth Lakes, California earthquakes were selected and scaled to fit the target horizontal acceleration response spectra.     

地震危险性、地震危害性和地震易损性

做好城市防震减灾工作 ,关键是对未来可能遭遇地震灾害定量化预测。本文论述地震灾害定量化的三要素 :地震危险性 ,地震危害性和地震易损性 ,它们概念完全不同 ,且很容易混淆 ,但又存在因与果的关系。本文还简要介绍厦门市地震科技工作者开展的“闽南地区综合防震减灾示范工程”。     

Seismic hazard assessment in Jordan and its vicinity

The earthquake hazard in Jordan and its vicinity is assessed on the basis of probabilistic methods. For this purpose, an updated earthquake catalog is compiled which covers the period between AD 1–1989. The earthquakes lie between latitudes 27.0°-35.5° N and longitudes 32.0°-39.0° E. Thirteen seismic zones are defined on a regional seismic and tectonic map presented for the area. Point-source and line-source models are used. The seismic hazard parameters, namely, theb-parameter (of the Gutenberg-Richter relation),m ₁ (the upper bound magnitude), and ₄ (the annual rate of occurrence of earthquakes with local magnitudeM _L 4.0) are calculated for each zone. The results of the seismic hazard assessment are displayed as iso-acceleration contours expected to be exceeded during typical economic life times of structures, i.e. 50 and 100 years. For each model, two seismic hazard maps are derived. In order to determine the importance of the South-eastern Mediterranean zone and the north part of the Red Sea zone from a seismic hazard point of view for Jordan, one seismic hazard map which corresponds to 50 years' economic life for every model, excluding the seismicity of these zones, is derived.     

The importance of paleoseismology in seismic hazard studies for critical facilities

This paper is intended to provide a perspective on the use of paleoseismological studies in the seismic hazard assessment of critical facilities, such as dams, chemical/petrochemical facilities and nuclear power plants. In particular, the use of data obtained from paleoseismological studies for probabilistic seismic hazard analyses, when the required probabilities of exceedance are very low (e.g. 10^{− 6}–10^{− 7}) is considered. Recent revisions to the IAEA Safety Standards that provide guidance to Member States in their work related to the seismic safety of nuclear power plants are presented to illustrate the importance of this emerging discipline.     

Overview of Taiwan Earthquake Loss Estimation System

The National Science Council (NSC) of Taiwan started the HAZ-Taiwan project in 1998 to promote researches on seismic hazard analysis, structural damage assessment, and socio-economic loss estimation. The associated application software, “Taiwan Earthquake Loss Estimation System (TELES)”, integrates various inventory data and analysis modules to fulfill three objectives. First, it helps to obtain reliable estimates of seismic hazards and losses soon after occurrence of large earthquakes. Second, it helps to simulate earthquake scenarios and to provide useful estimates for local governments or public services to propose their seismic disaster mitigation plans. Third, it helps to provide catastrophic risk management tools, such as proposing the seismic insurance policy for residential buildings. This paper focuses on the development and application of analysis modules used in early loss estimation system. These modules include assessments of ground motion intensity, soil liquefaction potential, building damage and casualty.     

Seismic Hazard Assessment in the Area of Mystras-Sparta,South Peloponnesus,Greece, Based on Local Seismotectonic,Seismic, Geologic Information and on Different Models of Rupture Propagation

Strong seismic events once again confirm the view that great destructive earthquakes are produced by the reactivation of pre-existing faults although they have usually remained inactive for many, perhaps thousands of years. It is evident that such active seismogenic zones, with little or no seismicity, have presumably been ignored in the determination of the region's seismic hazard. At south Peloponnesus, Greece, is situated at Taygetos mountain. At its eastern front lies a large normal fault system, the southern segment being the Sparta fault. This area has been characterized by low seismicity for the last 25 centuries. However, during the 6th and 5th centuries B.C. several destructive earthquakes have been reported. That of 464 B.C., was the most destructive and devastated the city of Sparta. Detailed morphotectonic observations of this area, suggest that the earthquake of 464 B.C. could be related to the most recent reactivation of this fault. The ground accelerations that would be produced by a future activation of the Sparta fault, were calculated, by applying a method which takes into account information mainly from the seismotectonic parameters of the Sparta fault, the rupture pattern, the properties of the propagation medium and the local ground conditions. Moreover, these results were compared with those of other independent studies based mainly on the seismic data of the area. This method estimated greater expected values of ground acceleration than those computed by the conventional seismic hazard methods. The highest values correspond to the activation of the Sparta fault either in a unilateral rupture, which would start from the southernmost point of the fault, or in a circular one. Furthermore, an increase is observed of the order of 50% in the ground acceleration values in unconsolidated soft ground in relation to the corresponding values of hard ground. This revised version was published online in June 2006 with corrections to the Cover Date.     

基于MAPGIS的概率性地震危险性分析

对工程场地的地震危险性分析是地震安全性评价的主要方法, 并且为工程师提供抗震参数.在分析国内概率性地震危险性分析(PSHA) 方法基础上, 提出了基于GIS的概率性地震危险性分析的可行性方案.通过MAPGIS二次开发编写了地震危险性分析程序, 并以三峡坝区某工程场地为例对程序进行了测试.基于MAPGIS的概率性地震危险性分析程序提供友好的人机交互界面, 提高了地震危险性分析的可操作性, 更重要的是帮助用户从空间数据中挖掘更多的信息.      

First Order Seismic Microzonation of Delhi,India Using Geographic Information System (GIS)

A first order seismic microzonation map of Delhi is prepared using five thematic layers viz., Peak Ground Acceleration (PGA) contour, different soil types at 6 m depth, geology, groundwater fluctuation and bedrock depth, integrated on GIS platform. The integration is performed following a pair-wise comparison of Analytical Hierarchy Process (AHP), wherein each thematic map is assigned weight in the 5-1 scale: depending on its contribution towards the seismic hazard. Following the AHP, the weightage assigned to each theme are: PGA (0.333), soil (0.266), geology (0.20), groundwater (0.133) and bedrock depth (0.066). The thematic vector layers are overlaid and integrated using GIS. On the microzonation theme, the Delhi region has been classified into four broad zones of vulnerability to the seismic hazard. They are very high (> 52%), high (38–52%), moderate (23–38%) and less ( < 23%) zones of seismic hazard. The “very high” seismic hazard zone is observed where the maximum PGA varies from 140 to 210 gal for a finite source model of M_w 8.5 in the central seismic gap. A site amplification study from local and regional earthquakes for Delhi region using Delhi Telemetry Network data shows a steeper site response gradient in the eastern side of the Yamuna fluvial deposits at 1.5 Hz. The ‘high’ seismic hazard zone occupies most of the study area where the PGA value ranges from 90 to 140 gal. The ‘moderate’ seismic hazard zone occurs on either side of the Delhi ridge with PGA value varying from 60 to 90 gal. The ‘less’ seismic hazard zone occurs in small patches distributed along the study area with the PGA value less than 60 gal. Site response studies, PGA distribution and destruction pattern of the Chamoli earthquake greatly corroborate the seismic hazard zones estimated through microzonation on GIS platform and also establishes the methodology incorporated in this study.     

新构造、活动构造与地震地质

新构造、活动构造和地震地质研究都是开展地震危险性评价的重要基础性地质工作。在综述新构造、活动构造和地震地质的基本含义、相互联系与区别、主要工作内容及方法的基础上,简要回顾了国内外在相关研究领域的主要进展,提出了中国活动构造与地震地质工作中应注意的主要问题和对未来工作的几点建议。最后,重点介绍了青藏高原东南缘开展活动构造体系和玉树地区活动断裂与地震地质调查研究工作所取得的主要进展与成果。     

四川康定城市地震地质及地质灾害风险评价

四川康定处于鲜水河、龙门山及安宁河三条地震带的＂Y＂字形交汇部位,境内地震活动频繁,地震基本烈度达IX度。康定城区位于鲜水河断裂南东段的雅拉河-康定断裂,城区几乎处于康定断裂的次级断层围限之中,断裂发育,岩体破碎,在降雨、地震和人类经济活动的影响下,地质灾害十分发育。康定城发育地质灾害25处,其中滑坡14处,崩塌3处,泥石流5处,不稳定斜坡3处。地质灾害多沿城区河流两岸斜坡和支沟分布,对城区居民的生命财产安全构成巨大的威胁,一旦发生强烈地震及其产生的次生地质灾害,将对康定城市带来灾难性的毁坏。近年来随着城镇化建设速度的加快,城市规模不断扩大,城市容量空前扩张,加之活动断裂活动性增加,城市安全风险增大,正确合理规划城市规模,科学地开展城市环境地质调查,制定有效的减灾防灾预案是确保城市安全的重要方法。     

Probabilistic Analysis of Tsunami Hazards*

Determining the likelihood of a disaster is a key component of any comprehensive hazard assessment. This is particularly true for tsunamis, even though most tsunami hazard assessments have in the past relied on scenario or deterministic type models. We discuss probabilistic tsunami hazard analysis (PTHA) from the standpoint of integrating computational methods with empirical analysis of past tsunami runup. PTHA is derived from probabilistic seismic hazard analysis (PSHA), with the main difference being that PTHA must account for far-field sources. The computational methods rely on numerical tsunami propagation models rather than empirical attenuation relationships as in PSHA in determining ground motions. Because a number of source parameters affect local tsunami runup height, PTHA can become complex and computationally intensive. Empirical analysis can function in one of two ways, depending on the length and completeness of the tsunami catalog. For site-specific studies where there is sufficient tsunami runup data available, hazard curves can primarily be derived from empirical analysis, with computational methods used to highlight deficiencies in the tsunami catalog. For region-wide analyses and sites where there are little to no tsunami data, a computationally based method such as Monte Carlo simulation is the primary method to establish tsunami hazards. Two case studies that describe how computational and empirical methods can be integrated are presented for Acapulco, Mexico (site-specific) and the U.S. Pacific Northwest coastline (region-wide analysis). * The U.S. Government’s right to retain a non-exclusive, royalty-free license in and to any copyright is acknowledged.