Seismic risk assessment for developing countries: Pakistan as a case study

Modern Earthquake Risk Assessment (ERA) methods usually require seismo-tectonic information for Probabilistic Seismic Hazard Assessment (PSHA) that may not be readily available in developing countries. To bypass this drawback, this paper presents a practical event-based PSHA method that uses instrumental seismicity, available historical seismicity, as well as limited information on geology and tectonic setting. Historical seismicity is integrated with instrumental seismicity to determine the long-term hazard. The tectonic setting is included by assigning seismic source zones associated with known major faults. Monte Carlo simulations are used to generate earthquake catalogues with randomized key hazard parameters. A case study region in Pakistan is selected to demonstrate the effectiveness of the method. The results indicate that the proposed method produces seismic hazard maps consistent with previous studies, thus being suitable for generating such maps in regions where limited data are available. The PSHA procedure is developed as an integral part of an ERA framework named EQRAM. The framework is also used to determine seismic risk in terms of annual losses for the study region.     

Probabilistic seismic hazard assessment for South-Eastern France

The accurate evaluation and appropriate treatment of uncertainties is of primary importance in modern probabilistic seismic hazard assessment (PSHA). One of the objectives of the SIGMA project was to establish a framework to improve knowledge and data on two target regions characterized by low-to-moderate seismic activity. In this paper, for South-Eastern France, we present the final PSHA performed within the SIGMA project. A new earthquake catalogue for France covering instrumental and historical periods was used for the calculation of the magnitude-frequency distributions. The hazard model incorporates area sources, smoothed seismicity and a 3D faults model. A set of recently developed ground motion prediction equations (GMPEs) from global and regional data, evaluated as adequately representing the ground motion characteristics in the region, was used to calculate the hazard. The magnitude-frequency distributions, maximum magnitude, faults slip rate and style-of-faulting are considered as additional source of epistemic uncertainties. The hazard results for generic rock condition (Vs30 = 800 m/s) are displayed for 20 sites in terms of uniform hazard spectra at two return periods (475 years and 10,000 years). The contributions of the epistemic uncertainties in the ground motion characterizations and in the seismic source characterization to the total hazard uncertainties are analyzed. Finally, we compare the results with existing models developed at national scale in the framework of the first generation of models supporting the Eurocode 8 enforcement, (MEDD 2002 and AFPS06) and at the European scale (within the SHARE project), highlighting significant discrepancies at short return periods.     

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.     

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

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

Seismic Inventory Simulation and Conformance Test Based on the Monte Carlo Method——A Case Study of the Fenhe-Weihe Belt

Artificial earthquake catalogue simulation is one of the ways to effectively improve the incompleteness of the existing earthquake catalogue,the scarcity of large earthquake records and the improvement of seismological research.Based on the Poisson distribution model of seismic activity and the Gutenberg-Richter magnitude-frequency relationship, the Monte Carlo method which can describe the characteristics of the stochastic nature and the physical experiment process is used.This paper simulates the future seismic catalogues of the Fenhe-Weihe seismic belt of different durations and conducts statistical tests on them. The analysis shows that the simulation catalogue meets the set seismic activity parameters and meets the Poisson distribution hypothesis,which can obtain a better simulated earthquake catalogues that meets the seismic activity characteristics.According to the simulated earthquake catalogues,future earthquake trends in this region are analyzed to provide reference for seismic hazard analysis.     

A probabilistic seismic hazard assessment for the Turkish territory: part II—fault source and background seismicity model

Over the years, several local and regional seismic hazard studies have been conducted for the estimation of the seismic hazard in Turkey using different statistical processing tools for instrumental and historical earthquake data and modeling the geologic and tectonic characteristics of the region. Recently developed techniques, increased knowledge and improved databases brought the necessity to review the national active fault database and the compiled earthquake catalogue for the development of a national earthquake hazard map. A national earthquake strategy and action plan were conceived and accordingly with the collaboration of the several institutions and expert researchers, the Revision of Turkish Seismic Hazard Map Project (UDAP-Ç-13-06) was initiated, and finalized at the end of 2014. The scope of the project was confined to the revision of current national seismic hazard map, using the state of the art technologies and knowledge of the active fault, earthquake database, and ground motion prediction equations. The following two seismic source zonation models are developed for the probabilistic earthquake hazard analysis: (1) Area source model, (2) Fault and spatial smoothing seismic source model (FSBCK). In this study, we focus on the development and the characterization of the Fault Source model, the background spatially smoothed seismicity model and intrinsic uncertainty on the earthquake occurrence-rates-estimation. Finally, PSHA results obtained from the fault and spatial smoothed seismic source model are presented for 43, 72, 475 and 2475 years return periods (corresponding to 69, 50, 10, and 2% probability of exceedance in 50 years) for PGA and 5% damped spectral accelerations at 0.2 and 1.0 s.     

Seismogenic potential and earthquake hazard assessment in the Tell Atlas of Algeria

Seismotectonic zonation studies in the Tell Atlas of Algeria, a branch of the Africa-Eurasia plate boundary, provide a valuable input for deterministic seismic hazard calculations. We delineate a number of seismogenic zones from causal relationships established between geological structures and earthquakes and compile a working seismic catalogue mainly from readily available sources. To this catalogue, for a most rational and best-justified hazard analysis, we add estimates of earthquake size translated from active faulting characteristics. We assess the regional seismic hazard using a deterministic procedure based on the computation of complete synthetic seismograms (up to 1 Hz) by the modal summation technique. As a result, we generate seismic hazard maps of maximum velocity, maximum displacement, and design ground acceleration that blend information from geology, historical seismicity and observational seismology, leading to better estimates of the earthquake hazard throughout northern Algeria. Our analysis and the resulting maps illustrate how different the estimate of seismic hazard is based primarily on combined geologic and seismological data with respect to the one for which only information from earthquake catalogues has been used.     

Effects of seismic declustering on seismic hazard assessment: a sensitivity study using the Turkish earthquake catalogue

Earthquakes trigger other earthquakes and build up clusters in space and time that in turn create a bias in seismic catalogues. Therefore, declustering is considered as a prerequisite in seismic studies, particularly for probabilistic seismic hazard analysis, not only to eliminate the bias but also to decouple mainshocks and triggered events. However, a declustering process is not a straightforward task due to the complex nature of earthquake phenomena. There exist several declustering methods that mostly employ subjective rules to distinguish between background seismicity and offsprings. Eventually, the final declustered catalogues usually deviate significantly according to the employed method. This issue is raising some concerns, such as how to select the most suitable declustering algorithm, or to assess how this selection affects seismic hazard assessment. In consequence, the main goal of this paper is to quantify the sensitivity of seismic hazard assessments to different declustering techniques. Accordingly, the recently compiled Turkish earthquake catalogue was declustered by making use of three declustering algorithms. A total of six declustered catalogues, two catalogues per method, one by implementing the default input parameters, and one by altering the free input parameters of the employed methods, were produced. The clusters of selected earthquakes were studied in terms of the spatial–temporal distribution of earthquake sequences. A sensitivity analysis was conducted through the major steps of seismic hazard assessment for Istanbul metropolitan city. The seismicity of Istanbul and surroundings was modeled on the basis of four areal source zones. Comparative studies showed that, while the selected declustering algorithm did not significantly affect the completeness periods of moderate to large size earthquakes, it considerably altered those of small magnitude events (e.g. M_w 4.3–5.2) and consequently the recurrence parameters of the source zones. Depending on the declustering algorithm and input parameters, the activity rate was observed to vary up to a factor of two. The differences in the declustered catalogues obtained from different declustering approaches resulted in considerable variations in seismic hazard estimations. The hazard maps at return periods of 475 and 2475 years indicated that peak ground acceleration values may vary up to 20% at some locations. Moreover, the differences in 5% damped elastic spectral accelerations at T = 0.2 for the return periods of 475 and 2475 years are about 18 and 12%, respectively, on the southern shores of Istanbul where the highest hazard levels are observed.     

Seismic source identification and characterisation for probabilistic seismic hazard analyses conducted in the Buller–NW Nelson Region, South Island, New Zealand

A seismic source model is presented for use in probabilistic seismic hazard analyses to be conducted for sites within the Buller–NW Nelson region of New Zealand. The application of common probabilistic seismic hazard analysis (PSHA) methodology for sites in this region has been complicated by the long-held suspicion that the observed rates of seismic activity are high and not representative of long-term earthquake activity. However, recent analyses of geological, seismicity and geodetic data indicate that the extent of this anomaly may have been overestimated and that current rates of seismic activity within this region are likely to continue into the foreseeable future. Probable bounds for the most appropriate long-term rates of seismic activity are estimated after considering all available sources of constraint. These include geodetic analyses, plate-motion modelling, finite element modelling, structural geological considerations, paleoseismic information, tree-ring analyses, precarious rock information, observed seismicity and fundamental mechanics. A suite of fault sources is identified, and the observed seismicity is partitioned between these sources and a background source using Bayesian inference, and then analysed to obtain a magnitude–frequency distribution for each seismic source. The annual moment release rate for the region, resulting from the identified and characterised sources, is shown to be consistent with available constraints. Consequently, it is demonstrated that the observed seismicity in the Buller–NW Nelson region can be used to model future earthquake occurrence within the region and that standard PSHA may therefore be implemented within the region.     

Uncertainty in recurrence rates of large magnitude events due to short historic catalogs

Seismic hazard analysis requires knowledge of the recurrence rates of large magnitude earthquakes that drive the hazard at low probabilities of interest for seismic design. Earthquake recurrence is usually determined through studies of the historic earthquake catalogue for a given region. Reliable historic catalogues generally span time periods of 100–200 years in North America, while large magnitude events (M?≥?7) have recurrence rates on the order of hundreds or thousands of years in many areas, resulting in large uncertainty in recurrence rates for large events. Using Monte Carlo techniques and assuming typical recurrence parameters, we simulate earthquake catalogues that span long periods of time. We then split these catalogues into smaller catalogues spanning 100–200 years that mimic the length of historic catalogues. For each of these simulated “historic” catalogues, a recurrence rate for large magnitude events is determined. By comparing recurrence rates from one historic-length catalogue to another, we quantify the uncertainty associated with determining recurrence rates from short historic catalogues. The use of simulations to explore the uncertainty (rather than analytical solutions) allows us flexibility to consider issues such as the relative contributions of aleatory versus epistemic uncertainty, and the influence of fitting method, as well as lending insight into extreme-event statistics. The uncertainty in recurrence rates of large (M?>?7) events is about a factor of two in regions of high seismicity, due to the shortness of historic catalogues. This uncertainty increases greatly with decreasing seismic activity. Uncertainty is dependent on the length of the catalogue as well as the fitting method used (least squares vs. maximum likelihood). Examination of 90th percentile recurrence rates reveals that epistemic uncertainty in the true parameters may cause recurrence rates determined from historic catalogues to be uncertain by a factor greater than 50.     

Seismicity assessment using earthquake catalogues with uncertain and incomplete data: probabilistic formulation

A novel generalized probabilistic formulation is proposed to assess seismicity using earthquake catalogues with uncertain and incomplete data. The seismicity, described by the complete exceedance rate of magnitudes, is estimated starting from a consistent incomplete exceedance rate which is rationally linked to the catalogue data. Complete and incomplete exceedance rates are represented by similar functional forms and they are related by a completeness function, which expresses the probability that an event is included in a data set. Completeness is considered uncertain and it is defined by a suitable, continuous, analytical, magnitude dependent function. The importance of this work lies on its applicability because it can be useful in seismic zones where information about seismic activity is scarce or simply when the catalogue is incomplete in a range of magnitudes that can have a significant influence on the seismic hazard analysis and on the resulting seismic risk assessment. Moreover, it can also be applied in the common case when the catalogue is considered complete above a given magnitude threshold. Numerical examples are presented to illustrate the influence of catalogue incompleteness on the complete exceedance rate estimations. In companion papers, attention is focused on the estimation of completeness probabilities of available catalogues and on parameter estimation of the exceedance rate functions.     

Rating of seismicity and reconstruction of the fault geometries in northern Tien Shan within the project “Seismic Hazard Assessment for Almaty”

The project “Seismic Hazard Assessment for Almaty” has a main objective to improve existing seismic hazard maps for the region of northern Tien Shan and especially for the surroundings of Almaty and to generate a new geodynamic model of the region.In the first step a composite seismic catalogue for the northern Tien Shan region was created, which contains about 20,000 events and is representative for strong earthquakes for the period back to the year 500. For the period of instrumental observation 1911–2006 the catalogue contains data for earthquakes with a body wave magnitude larger than 4. For smaller events with magnitudes up to 2.2 the data are only available since 1980. The composite catalogue was created on the basis of several catalogues from the United States Geologic Survey (USGS), local catalogues from the Kazakh National Data Centre (KNDC) and the USSR earthquake catalogue. Due to the different magnitudes used in several catalogues a magnitude conversion was necessary.Event density maps were created to rate the seismicity in the region and to identify seismic sources. Subsurface fault geometries were constructed using tectonic model which uses fault parallel material flow and is constrained by GPS data. The fault geometry should improve the estimation of the expected seismic sources from seismic density maps.First analysis of the earthquake catalogue and the density maps has shown that nearly all large events are related to fault systems. Annual seismicity distribution maps suggest different processes as the cause for the seismic events. Apart from tectonics, also fluids play a major part in triggering of the earthquakes.Beneath the Issyk-Kul basin the absence of strong seismic activity suggests aseismic sliding at the flat ramp in a ductile crust part and low deformation within the stable Issyk-Kul micro-continent which underthrust the northern ranges of Tien Shan. First results suggest a new partition of the region in tectonic units, whose bounding faults are responsible for most of the seismic activity.     

K-means cluster analysis and seismicity partitioning for Pakistan

Pakistan and the western Himalaya is a region of high seismic activity located at the triple junction between the Arabian, Eurasian and Indian plates. Four devastating earthquakes have resulted in significant numbers of fatalities in Pakistan and the surrounding region in the past century (Quetta, 1935; Makran, 1945; Pattan, 1974 and the recent 2005 Kashmir earthquake). It is therefore necessary to develop an understanding of the spatial distribution of seismicity and the potential seismogenic sources across the region. This forms an important basis for the calculation of seismic hazard; a crucial input in seismic design codes needed to begin to effectively mitigate the high earthquake risk in Pakistan. The development of seismogenic source zones for seismic hazard analysis is driven by both geological and seismotectonic inputs. Despite the many developments in seismic hazard in recent decades, the manner in which seismotectonic information feeds the definition of the seismic source can, in many parts of the world including Pakistan and the surrounding regions, remain a subjective process driven primarily by expert judgment. Whilst much research is ongoing to map and characterise active faults in Pakistan, knowledge of the seismogenic properties of the active faults is still incomplete in much of the region. Consequently, seismicity, both historical and instrumental, remains a primary guide to the seismogenic sources of Pakistan. This study utilises a cluster analysis approach for the purposes of identifying spatial differences in seismicity, which can be utilised to form a basis for delineating seismogenic source regions. An effort is made to examine seismicity partitioning for Pakistan with respect to earthquake database, seismic cluster analysis and seismic partitions in a seismic hazard context. A magnitude homogenous earthquake catalogue has been compiled using various available earthquake data. The earthquake catalogue covers a time span from 1930 to 2007 and an area from 23.00° to 39.00°N and 59.00° to 80.00°E. A threshold magnitude of 5.2 is considered for K-means cluster analysis. The current study uses the traditional metrics of cluster quality, in addition to a seismic hazard contextual metric to attempt to constrain the preferred number of clusters found in the data. The spatial distribution of earthquakes from the catalogue was used to define the seismic clusters for Pakistan, which can be used further in the process of defining seismogenic sources and corresponding earthquake recurrence models for estimates of seismic hazard and risk in Pakistan. Consideration of the different approaches to cluster validation in a seismic hazard context suggests that Pakistan may be divided into K?=?19 seismic clusters, including some portions of the neighbouring countries of Afghanistan, Tajikistan and India.     

新版地震区划图地震活动性模型与参数确定

地震活动性模型和地震动预测模型是概率地震危险性分析的两个核心。在新版地震区划图中,依据板内地震活动空间不均匀性分布的特点,在概率地震危险性分析方法(CPSHA)中采用了由地震统计区、背景潜在震源区和构造潜在震源区构成的三级层次性潜在震源区模型,并构建了相应的地震活动性模型。本文在论述CPSHA方法及其地震活动性模型基本概念的基础上,重点介绍了新版地震区划图地震活动性模型的三级潜在震源区模型的构成、地震活动性假定和基本特点,同时,也对新版地震区划图地震活动性模型的重要参数确定思路、方法与结果进行了介绍。本文将为更好地认识与理解我国新版地震动参数区划图提供有益的参考。     

Evaluation of the seismic hazard for 20 cities in Romania using Monte Carlo based simulations

This work focuses on the evaluation of the seismic hazard for Romania using earthquake catalogues generated by a Monte Carlo approach. The seismicity of Romania can be attributed to the Vrancea intermediate-depth seismic source and to 13 other crustal seismic sources. The recurrence times of large magnitude seismic events (both crustal and subcrustal), as well as the moment release rates are computed using simulated earthquake catalogues. The results show that the largest contribution to the overall moment release for the crustal seismic sources is from the seismic regions in Bulgaria, while the seismic regions in Romania contribute less than 5% of the overall moment release. In addition, the computations show that the moment release rate for the Vrancea subcrustal seismic source is about ten times larger than that of all the crustal seismic sources. Finally, the Monte Carlo approach is used to evaluate the seismic hazard for 20 cities in Romania with populations larger than 100,000 inhabitants. The results show some differences between the seismic hazard values obtained through Monte-Carlo simulation and those in the Romanian seismic design code P100-1/2013, notably for cities situated in the western part of Romania that are influenced by local crustal seismic sources.     

Geodetic Strain Observations and Return Period of the Strongest Earthquakes of a Given Seismic Source Zone

We present the basis for a method for estimating the return period of large and medium earthquakes that is independent of current deterministic and probabilistic approaches. The two standard techniques of seismic hazard assessment??probabilistic seismic hazard assessment (PSHA) and deterministic seismic hazard assessment (DSHA)??suffer from limited knowledge of seismic prehistory. A further weakness of PSHA is its requirement of homogeneous seismic activity within a seismic zone. Moreover, PSHA and DSHA were developed for seismically active areas and, thus, cannot reliably be used in areas of medium and low activity. In this paper we propose the combined use of geodetic strain rate data and the seismic moment data set determined for past seismic events. This combination represents a new and independent approach to estimation of future seismic activity. Using a modified version of Kostrov??s (Phys Solid Earth 1:23?C40, 1974) equation and the catalogue of seismic moments, the minimum return period of the strongest earthquakes of a source area is estimated.     

Seismic hazard assessment of Tehran,Iran with emphasis on near-fault rupture directivity effects

Many destructive earthquakes happened in Tehran, Iran in the last centuries. The existence of active faults like the North Tehran is the main cause of seismicity in this city. According to past investigations, it is estimated that in the scenario of activation of the North Tehran fault, many structures in Tehran will collapse. Therefore, it is necessary to incorporate the near field rupture directivity effects of this fault into the seismic hazard assessment of important sites in Tehran. In this study, using calculations coded in MATLAB,Probabilistic Seismic Hazard Analysis(PSHA) is conducted for an important site in Tehran. Following that, deaggregation technique is performed on PSHA and the contribution of seismic scenarios to hazard is obtained in the range of distance and magnitude. After identifying the North Tehran fault as the most hazardous source affecting the site in 10000-year return period, rupture directivity effects of this fault is incorporated into the seismic hazard assessment using Somerville et al.(1997) model with broadband approach and Shahi and Baker(2011) model with narrowband approach. The results show that the narrowband approach caused a 27% increase in the peak of response spectrum in 10000-year return period compared with the conventional PSHA. Therefore, it is necessary to incorporate the near fault rupture directivity effects into the higher levels of seismic hazard assessment attributed to important sites.     

基于特大地震发生率的川滇地区地震危险性预测新模型

川滇地区是我国地震危险性较高的地区之一.本文基于对特大强震的风险性考虑,使用全球地震模型OpenQuake软件,建立了川滇地区地震危险性预测新模型.首先根据构造特征划分多个震源分区,并整理出这些震源分区内断层活动特征与滑动速率;基于震源分区和断层模型,使用GPS应变率转换成的锥形古登堡-里克特关系作为整个区域的地震积累率,并允许超过历史最大震级的特大地震的出现,结合活动断层滑动速率所积累的地震发生率,给出震源分区内断层地震源和背景地震源的地震发生率的比率分配关系;在活动断层分段上,保留了大型断裂或其主要部分,没有根据小的阶区来对断层进行详细分段,以便分配特大地震发生率;并使用地震率平滑方法分配背景地震发生率.最后在OpenQuake中加入地震动预测方程,计算出了川滇地区的PGA分布图,为区域地震危险性提供科学依据.     

Updated earthquake catalogue for seismic hazard analysis in Pakistan

A reliable and homogenized earthquake catalogue is essential for seismic hazard assessment in any area. This article describes the compilation and processing of an updated earthquake catalogue for Pakistan. The earthquake catalogue compiled in this study for the region (quadrangle bounded by the geographical limits 40–83° N and 20–40° E) includes 36,563 earthquake events, which are reported as 4.0–8.3 moment magnitude (M_W) and span from 25 AD to 2016. Relationships are developed between the moment magnitude and body, and surface wave magnitude scales to unify the catalogue in terms of magnitude M_W. The catalogue includes earthquakes from Pakistan and neighbouring countries to minimize the effects of geopolitical boundaries in seismic hazard assessment studies. Earthquakes reported by local and international agencies as well as individual catalogues are included. The proposed catalogue is further used to obtain magnitude of completeness after removal of dependent events by using four different algorithms. Finally, seismicity parameters of the seismic sources are reported, and recommendations are made for seismic hazard assessment studies in Pakistan.     

Seismic Hazard Assessment: Issues and Alternatives

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