USGS Science⌕ Search

SEARCH · USGS Science

Results for “earthquake hazards”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,297 records · Page 72Linked to original sources
USGS2017

A fault‐based model for crustal deformation in the western United States based on a combined inversion of GPS and geologic inputs

We develop a crustal deformation model to determine fault‐slip rates for the western United States (WUS) using the Zeng and Shen (2014) method that is based on a combined inversion of Global Positioning System (GPS) velocities and geological slip‐rate constraints. The model consists of six blocks with boundaries aligned along major faults in California and the Cascadia subduction zone, which are represented as buried dislocations in the Earth. Faults distributed within blocks have their geometrical structure and locking depths specified by the Uniform California Earthquake Rupture Forecast, version 3 (UCERF3) and the 2008 U.S. Geological Survey National Seismic Hazard Map Project model. Faults slip beneath a predefined locking depth, except for a few segments where shallow creep is allowed. The slip rates are estimated using a least‐squares inversion. The model resolution analysis shows that the resulting model is influenced heavily by geologic input, which fits the UCERF3 geologic bounds on California B faults and ±one‐half of the geologic slip rates for most other WUS faults. The modeled slip rates for the WUS faults are consistent with the observed GPS velocity field. Our fit to these velocities is measured in terms of a normalized chi‐square, which is 6.5. This updated model fits the data better than most other geodetic‐based inversion models. Major discrepancies between well‐resolved GPS inversion rates and geologic‐consensus rates occur along some of the northern California A faults, the Mojave to San Bernardino segments of the San Andreas fault, the western Garlock fault, the southern segment of the Wasatch fault, and other faults. Off‐fault strain‐rate distributions are consistent with regional tectonics, with a total off‐fault moment rate of 7.2 × 10 18 "> 7.2 × 10 18 and 8.5 × 10 18    N · m / year "> 8.5 × 10 18 N ⋅ m / year for California and the WUS outside California, respectively.

Bulletin of the Seismological Society of America↗
USGS2018

Characterizing and imaging sedimentary strata using depth-converted spectral ratios: An example from the Atlantic Coastal Plain of the Eastern U.S.

Unconsolidated, near-surface sediments can have a profound influence on the amplitudes and frequencies of ground shaking during earthquakes, and these effects should be accounted for when using amplitude observations for seismic hazard assessments. This study explores methods to use teleseismic arrivals recorded on linear receiver arrays to characterize widespread, shallow sedimentary deposits, including estimation of the velocities, determination of the fundamental resonance peaks, and imaging of the major reflectors. The examples used are the extensive Atlantic Coastal Plain (ACP) and associated Mississippi Embayment (ME) strata of the Central and Eastern United States. The large contrast in material properties at the bedrock surface beneath these sediments produces a strong fundamental resonance peak in the 0.2 to 4 Hz frequency range, which is estimated here by computing spectral ratios at each receiver site relative to bedrock sites at the ends of the receiver arrays. Sediment thicknesses derived from published contour maps made from drill hole data allow for the computation of average velocities to match the observed frequencies of resonance peaks with theoretical values at each receiver site, with the sloping bedrock surface allowing for computation of an average velocity versus depth function if horizontal layers are assumed. The velocity function is then used to convert the spectral ratios from frequency to depth, resulting in an image of the subsurface similar to that of a seismic reflection profile. The results demonstrate the use of teleseismic signals for characterizing and imaging shallow sedimentary strata.

Atlantic Coastal Plain↗
USGS2025

Deterministic physics-based earthquake sequence simulators match empirical ground-motion models and enable extrapolation to data poor regimes: Application to multifault multimechanism ruptures

We use the deterministic earthquake simulator RSQSim to generate complex sequences of ruptures on fault systems used for hazard assessment. We show that the source motions combined with a wave propagation code create surface ground motions that fall within the range of epistemic uncertainties for the Next Generation Attenuation‐West2 set of empirical models. We show the model is well calibrated where there are good data constraints, and has good correspondence in regions with fewer data constraints. We show magnitude, distance, and mechanism dependence all arising naturally from the same underlying friction. The deterministic physics‐based approach provides an opportunity for better understanding the physical origins of ground motions. For example, we find that reduced stress drops in shallow layers relative to constant stress drop with depth lead to peak ground velocities in the near field that better match empirical models. The simulators may also provide better extrapolations into regimes that are poorly empirically constrained by data because physics, rather than surface shaking data parameterizations, is underlying the extrapolations. Having shown the model is credible, we apply it to a problem where observations are lacking. We examine the case of crustal faults above a shallow subduction interface seen to break coseismically in simulations of the New Zealand fault system. These types of events were left out of consideration in the most recent New Zealand national seismic hazard model due to the modeling complexity and lack of observational data to constrain ground‐motion models (GMMs). Here, we show that in the model, by breaking up the coseismic crustal and interface rupturing fault motions into two separate subevents, and then recombining the resulting ground‐motion measures in a square‐root‐of‐sum‐of‐squares incoherent manner, we reproduce well the ground‐motion measures from the full event rupture. This provides a new method for extrapolating GMMs to more complex multifault ruptures.

Seismological Research Letters↗
USGS2014

Pulverization provides a mechanism for the nucleation of earthquakes at low stress on strong faults

An earthquake occurs when rock that has been deformed under stress rebounds elastically along a fault plane ( Gilbert, 1884 ; Reid, 1911 ), radiating seismic waves through the surrounding earth. Rupture along the entire fault surface does not spontaneously occur at the same time, however. Rather the rupture starts in one tiny area, the rupture nucleation zone, and spreads sequentially along the fault. Like a row of dominoes, one bit of rebounding fault triggers the next. This triggering is understood to occur because of the large dynamic stresses at the tip of an active seismic rupture. The importance of these crack tip stresses is a central question in earthquake physics. The crack tip stresses are minimally important, for example, in the time predictable earthquake model ( Shimazaki and Nakata, 1980 ), which holds that prior to rupture stresses are comparable to fault strength in many locations on the future rupture plane, with bits of variation. The stress/strength ratio is highest at some point, which is where the earthquake nucleates. This model does not require any special conditions or processes at the nucleation site; the whole fault is essentially ready for rupture at the same time. The fault tip stresses ensure that the rupture occurs as a single rapid earthquake, but the fact that fault tip stresses are high is not particularly relevant since the stress at most points does not need to be raised by much. Under this model it should technically be possible to forecast earthquakes based on the stress-renewaql concept, or estimates of when the fault as a whole will reach the critical stress level, a practice used in official hazard mapping ( Field, 2008 ). This model also indicates that physical precursors may be present and detectable, since stresses are unusually high over a significant area before a large earthquake.

Frontiers in Earth Science↗
USGS2019

Updates to USGS national seismic hazard model (NSHM) and design ground motion maps for 2020 NEHRP recommended provisions

This presentation summarizes the proposed updates to earthquake design ground motions for the 2020 edition of the NEHRP Recommended Seismic Provisions, expected to be incorporated into the ASCE 7-22 Standard. The implications of these updates on the values of design ground motions for example locations in both conterminous and nonconterminous U.S. cities are shown and discussed.

Conference Paper↗
USGS2013

Toward a time-dependent probabilistic seismic hazard analysis for Alaska

We report on a time-dependent seismic hazard analysis for Alaska and the Aleutians to complement our recently completed time-independent map. Whereas the time-independent map treats all sources as statistically independent, the time-dependent analysis is based on calculations of the conditional probability of occurrence for the next 50 years by using a Brownian Passage Time model for the seismic sources judged to be characteristic. We then consider how those probabilities are modified by coseismic and postseismic stress changes resulting from large regional earthquakes occurring from 1938 to 2002. Recombining the time-dependent probabilities with time-independent truncated Gutenberg–Richter and smoothed seismicity sources leads to our time-dependent probabilistic seismic hazard results. We find that when accounting for time dependence without stress changes, earthquake probabilities can be significantly altered, reducing probabilities to near zero or increasing them to several times the time-independent values. In addition, accounting for coseismic stress changes tends to have a local influence on earthquake probabilities, whereas postseismic effects can be far-reaching in both time and space. In sum, however, since we combine time-dependent and time-independent sources, the modification to seismic hazard is relatively minor, increasing or decreasing hazard adjacent to characteristic faults by about 10%. Most cities, located far from characteristic faults, are little affected.

Alaska↗
USGS1987

Failure of a massive earthquake-induced landslide dam in Papua New Guinea

In many areas of the world, landslides dams are both interesting natural phenomena and significant hazards. A few of the these natural blockages attain heights that rival or exceed those of the largest manmade dams. A landslide dam in its natural state differs from a constructed embankment dam in that it is composed of a heterogeneous mass of poorly consolidated earth debris and has no channeled spillway or other protected outlet for its impoundment. Having no outlet, landslide dams upon overtopping very commonly fail by breaching due to erosion by the overflowing stream. Some of the world's largest and most devastating floods have occurred because of the failure of large landslide dams (Schuster and Costa, 1986). This article discusses the recent occurrence of a large earthquake-induced landslide that dammed the Bairaman River in the interior of the island of New Britian, Papua New Guinea, and the subsequent overtopping and failure of this landslide dam.

New Britain↗
USGS2010

Extension in Mona Passage, Northeast Caribbean

As shown by the recent Mw 7.0 Haiti earthquake, intra-arc deformation, which accompanies the subduction process, can present seismic and tsunami hazards to nearby islands. Spatially-limited diffuse tectonic deformation within the Northeast Caribbean Plate Boundary Zone likely led to the development of the submerged Mona Passage between Puerto Rico and the Dominican Republic. GPS geodetic data and a moderate to high level of seismicity indicate that extension within the region is ongoing. Newly-collected high-resolution multibeam bathymetry and multi-channel seismic reflection profiles and previously-collected samples are used here to determine the tectonic evolution of the Mona Passage intra-arc region. The passage is floored almost completely by Oligocene-Pliocene carbonate platform strata, which have undergone submarine and subaerial erosion. Structurally, the passage is characterized by W- to NNW-trending normal faults that offset the entire thickness of the Oligo-Pliocene carbonate platform rocks. The orientation of these faults is compatible with the NE-oriented extension vector observed in GPS data. Fault geometry best fits an oblique extension model rather than previously proposed single-phase, poly-phase, bending-moment, or rotation extension models. The intersection of these generally NW-trending faults in Mona Passage with the N-S oriented faults of Mona Canyon may reflect differing responses of the brittle upper-crust, along an arc-forearc rheological boundary, to oblique subduction along the Puerto Rico trench. Several faults within the passage, if ruptured completely, are long enough to generate earthquakes with magnitudes on the order of Mw 6.5-7. ?? 2010.

Tectonophysics↗
USGS2020

Probabilistic regional-scale liquefaction triggering modeling using 3D Gaussian processes

Liquefaction is a major cause of coseismic damages, occurring irregularly over hundreds or thousands of square kilometers in large earthquakes. Large variations in the extent and location of liquefaction have been observed in recent earthquakes, motivating the need for prediction methods that consider the spatial heterogeneity of geologic deposits at a regional scale. Contemporary regional-scale liquefaction hazard analyses are typically performed using only surficial data, which does not address the complicated subsurface mechanics and spatial variability associated with artificial fill and natural soil deposits. In this study, we develop a probabilistic, regional-scale, subsurface model using data from hundreds of borings to better understand subsurface conditions that could influence liquefaction. We then use this subsurface sample database to train Gaussian process models, yielding 3D independent random fields of groundwater depth, soil plasticity, and penetration resistance for each geologic unit. We incorporate the Gaussian process models into probabilistic liquefaction triggering procedures, producing 3D estimates of the probability of liquefaction for an example study area in Portland, Oregon. Near sampling locations, the variance of the Gaussian process models approaches the variance of site-specific liquefaction triggering procedures. Conversely, when no sample data are nearby to condition a Gaussian process, the variance approaches the marginal variance of the entire recorded dataset. Thus, the procedure described in this study unifies probabilistic site-specific and regional-scale liquefaction triggering procedures and provides an important step towards quantitative liquefaction hazard assessments for regionally distributed infrastructures, such as levees, pipelines, roadways, and electrical transmission facilities.

Oregon↗
USGS2006

The influence of hazard models on GIS-based regional risk assessments and mitigation policies

Geographic information systems (GIS) are important tools for understanding and communicating the spatial distribution of risks associated with natural hazards in regional economies. We present a GIS-based decision support system (DSS) for assessing community vulnerability to natural hazards and evaluating potential mitigation policy outcomes. The Land Use Portfolio Modeler (LUPM) integrates earth science and socioeconomic information to predict the economic impacts of loss-reduction strategies. However, the potential use of such systems in decision making may be limited when multiple but conflicting interpretations of the hazard are available. To explore this problem, we conduct a policy comparison using the LUPM to test the sensitivity of three available assessments of earthquake-induced lateral-spread ground failure susceptibility in a coastal California community. We find that the uncertainty regarding the interpretation of the science inputs can influence the development and implementation of natural hazard management policies. Copyright ?? 2006 Inderscience Enterprises Ltd.

International Journal of Risk Assessment and Manag↗
USGS2010

Earthquake casualty models within the USGS Prompt Assessment of Global Earthquakes for Response (PAGER) system

Since the launch of the USGS’s Prompt Assessment of Global Earthquakes for Response (PAGER) system in fall of 2007, the time needed for the U.S. Geological Survey (USGS) to determine and comprehend the scope of any major earthquake disaster anywhere in the world has been dramatically reduced to less than 30 min. PAGER alerts consist of estimated shaking hazard from the ShakeMap system, estimates of population exposure at various shaking intensities, and a list of the most severely shaken cities in the epicentral area. These estimates help government, scientific, and relief agencies to guide their responses in the immediate aftermath of a significant earthquake. To account for wide variability and uncertainty associated with inventory, structural vulnerability and casualty data, PAGER employs three different global earthquake fatality/loss computation models. This article describes the development of the models and demonstrates the loss estimation capability for earthquakes that have occurred since 2007. The empirical model relies on country-specific earthquake loss data from past earthquakes and makes use of calibrated casualty rates for future prediction. The semi-empirical and analytical models are engineering-based and rely on complex datasets including building inventories, time-dependent population distributions within different occupancies, the vulnerability of regional building stocks, and casualty rates given structural collapse.

Book chapter↗
USGS2004

Stress interaction between subduction earthquakes and forearc strike-slip faults: Modeling and application to the northern Caribbean plate boundary

Strike-slip faults in the forearc region of a subduction zone often present significant seismic hazard because of their proximity to population centers. We explore the interaction between thrust events on the subduction interface and strike-slip faults within the forearc region using three-dimensional models of static Coulomb stress change. Model results reveal that subduction earthquakes with slip vectors subparallel to the trench axis enhance the Coulomb stress on strike-slip faults adjacent to the trench but reduce the stress on faults farther back in the forearc region. In contrast, subduction events with slip vectors perpendicular to the trench axis enhance the Coulomb stress on strike-slip faults farther back in the forearc, while reducing the stress adjacent to the trench. A significant contribution to Coulomb stress increase on strike-slip faults in the back region of the forearc comes from "unclamping" of the fault, i.e., reduction in normal stress due to thrust motion on the subduction interface. We argue that although Coulomb stress changes from individual subduction earthquakes are ephemeral, their cumulative effects on the pattern of lithosphere deformation in the forearc region are significant. We use the Coulomb stress models to explain the contrasting deformation pattern between two adjacent segments of the Caribbean subduction zone. Subduction earthquakes with slip vectors nearly perpendicular to the Caribbean trench axis is dominant in the Hispaniola segment, where the strike-slip faults are more than 60 km inland from the trench. In contrast, subduction slip motion is nearly parallel to the Caribbean trench axis along the Puerto Rico segment, where the strike-slip fault is less than 15 km from the trench. This observed jump from a strike-slip fault close to the trench axis in the Puerto Rico segment to the inland faults in Hispaniola is explained by different distributions of Coulomb stress in the forearc region of the two segments, as a result of the change from the nearly trench parallel slip on the Puerto Rico subduction interface to the more perpendicular subduction slip beneath Hispaniola. The observations and modeling suggest that subduction-induced strike-slip seismic hazard to Puerto Rico may be smaller than previously assumed but the hazard to Hispaniola remains high. Copyright 2004 by the American Geophysical Union.

Puerto Rico↗
USGS2023

2019 M7.1 Ridgecrest earthquake slip distribution controlled by fault geometry inherited from Independence dike swarm

Faults often form through reactivation of pre-existing structures, developing geometries and mechanical properties specific to the system’s geologic inheritance. Competition between fault geometry and other factors (e.g., lithology) to control slip at Earth’s surface is an open question that is central to our knowledge of fault processes and seismic hazards. Here we use remote sensing data and field observations to investigate the origin of the 2019 M 7.1 Ridgecrest, California, earthquake rupture geometry and test its impact on the slip distribution observed at Earth’s surface. Common geometries suggest the fault system evolved through reactivation of structures within the surrounding Independence dike swarm (IDS). Mechanical models testing a range of fault geometries and stress fields indicate that the inherited rupture geometry strongly controlled the M 7.1 earthquake slip distribution. These results motivate revisiting the development of other large-magnitude earthquake ruptures (1992 M 7.3 Landers, 1999 M 7.1 Hector Mine) and tectonic provinces within the IDS.

California↗
USGS2003

Fault interactions and large complex earthquakes in the Los Angeles area

Faults in complex tectonic environments interact in various ways, including triggered rupture of one fault by another, that may increase seismic hazard in the surrounding region. We model static and dynamic fault interactions between the strike-slip and thrust fault systems in southern California. We find that rupture of the Sierra Madre-Cucamonga thrust fault system is unlikely to trigger rupture of the San Andreas or San Jacinto strike-slip faults. However, a large northern San Jacinto fault earthquake could trigger a cascading rupture of the Sierra Madre-Cucamonga system, potentially causing a moment magnitude 7.5 to 7.8 earthquake on the edge of the Los Angeles metropolitan region.

Science↗
USGS2023

Ambient seismic noise tomography of the Kingdom of Saudi Arabia

Harrat Rahat is a Cenozoic volcanic field in the west-central part of the Kingdom of Saudi Arabia, 150 kilometers east of the Red Sea, and is the site of the most recent eruption in the country (1256 C.E.; 654 in the year of the Hijra). The city of Al Madīnah lies at the north end of Harrat Rahat, and its volcanic and seismic risks are frequently reassessed. In 2009 C.E. an earthquake swarm at Harrat Lunayyir, 200 km west-northwest of Al Madīnah, produced significant shaking and some building damage in nearby communities, which prompted a revision of seismic hazard models for the region. A network of seismic stations on this volcanically active western side of the Arabia Plate was installed, and stations were also added in the tectonically active northern part of the country. Although regional earthquakes may be used to determine the crustal structure of the western Arabia Plate, such crustal models are hindered by insufficient numbers of earthquakes in the stable plate interior. Tomography studies can be used to infer material properties of the subsurface, such as presence of partial melt, and are beneficial for volcanic hazard assessments. We use ambient seismic noise to compute Rayleigh and Love surface-wave dispersion maps between 5 and 12 second periods for a subset of seismic stations within and near northern Harrat Rahat. The surface-wave maps are inverted to produce shear-wave velocities using a neighborhood algorithm and interpolated into a pseudo three-dimensional model. The distributions of surface-wave and shear-wave velocities are heterogenous, varying from ±3–8 percent. However, low velocities are not restricted to Harrat Rahat. We observed a difference between Rayleigh- and Love-wave velocities that extends north of the site of the 1256 C.E. eruption and coincides with a low gravity anomaly. We obtain a shear-wave velocity increase of 10–15 percent between 15- and 25-kilometers depth, which is consistent with the presence of a transition between the felsic upper crust and the mafic lower crust of the Arabian Shield. The average shear-wave velocities of the upper and lower crust are estimated to be 3.64 and 3.95 kilometers per second using Rayleigh waves and 3.53 and 4.16 kilometers per second using Love waves, which are in good agreement with the results of other geophysical surveys in this area. The modest magnitude of the low-velocity anomalies within the crust and their locations extending well beyond the limits of Harrat Rahat indicate that they are not caused by a crustal magma chamber. If magma chambers exist, they are smaller than can be imaged with our seismological method (resolution corresponds to a 15-kilometer wavelength), deeper than 30 kilometers, or shallower than 5 kilometers with a small velocity contrast. The last possibility is discounted by weak-to-absent surface geothermal activity and by a lack of shallow seismicity that characterizes other areas with known shallow magmas, such as Hawaiʻi. We then expanded our work to the entire Arabian Shield with principally the same methodology to look at large-scale regional patterns. We found modest shear-wave velocity deviations on the order of only ±3 percent, which are within expected ranges for lithological variation.

Professional Paper↗
USGS2021

Earthquakes indicated magma viscosity during Kīlauea’s 2018 eruption

Magma viscosity strongly controls the style (for example, explosive versus effusive) of a volcanic eruption and thus its hazard potential, but can only be measured during or after an eruption. The identification of precursors indicative of magma viscosity would enable forecasting of the eruption style and the scale of associated hazards 1 . The unanticipated May 2018 rift intrusion and eruption of Kīlauea Volcano, Hawai‘i 2 displayed exceptional chemical and thermal variability in erupted lavas, leading to unpredictable effusion rates and explosivity. Here, using an integrated analysis of seismicity and magma rheology, we show that the orientation of fault-plane solutions (which indicate a fault’s orientation and sense of movement) for earthquakes preceding and accompanying the 2018 eruption indicate a 90-degree local stress-field rotation from background, a phenomenon previously observed only at high-viscosity eruptions 3 , and never before at Kīlauea 4,5,6,7,8 . Experimentally obtained viscosities for 2018 products and earlier lavas from the Pu‘u ‘Ō‘ō vents tightly constrain the viscosity threshold required for local stress-field reorientation. We argue that rotated fault-plane solutions in earthquake swarms at Kīlauea and other volcanoes worldwide provide an early indication that unrest involves magma of heightened viscosity, and thus real-time monitoring of the orientations of fault-plane solutions could provide critical information about the style of an impending eruption. Furthermore, our results provide insight into the fundamental nature of coupled failure and flow in complex multiphase systems.

Hawaii↗
USGS2013

Reverberations on the watery element: A significant tsunamigenic historical earthquake offshore the Carolina coast

We investigate an early nineteenth-century earthquake that has been previously cataloged but not previously investigated in detail or recognized as a significant event. The earthquake struck at approximately 4:30 a.m. LT on 8 January 1817 and was widely felt throughout the southeastern and mid-Atlantic United States. Around 11:00 a.m. the same day, an eyewitness described a 12-inch tide that rose abruptly and agitated boats on the Delaware River near Philadelphia. We show that the timing of this tide is consistent with the predicted travel time for a tsunami generated by an offshore earthquake 6–7 hours earlier. By combining constraints provided by the shaking intensity distribution and the tsunami observation, we conclude that the 1817 earthquake had a magnitude of low- to mid-M 7 and a location 800–1000 km offshore of South Carolina. Our results suggest that poorly understood offshore source zones might represent a previously unrecognized hazard to the southern and mid-Atlantic coast. Both observational and modeling results indicate that potential tsunami hazard within Delaware Bay merits consideration: the simple geometry of the bay appears to catch and focus tsunami waves. Our preferred location for the 1817 earthquake is along a diffuse northeast-trending zone defined by instrumentally recorded and historical earthquakes. The seismotectonic framework for this region remains enigmatic.

South Carolina↗
USGS2005

Comment on "How can seismic hazard around the New Madrid seismic zone be similar to that in California?" by Arthur Frankel

A recent paper in the September/October 2004 issue of Seismological Research Letters, “How Can Seismic Hazard around the New Madrid Seismic Zone Be Similar to That in California?”, by Arthur Frankel, concludes “I have demonstrated how probabilistic seismic hazard for New Madrid can be greater than that at San Francisco at low probabilities for PGA and similar at low probabilities for 5 Hz S.A. By low probabilities, I mean annual probabilities less than the reciprocal of the return time of the New Madrid characteristic source, that is, 1/500. This is a consequence of the higher ground motions, for PGA and 5 Hz S.A. (and other high-frequency measures), expected for large New Madrid characteristic earthquakes compared to San Andreas earthquakes with similar magnitudes.” His conclusion is confusing and not necessarily true. The fundamental differences in assessing seismic hazard between San Francisco and the New Madrid seismic zone are the uncertainties, especially the uncertainties associated with the ground-motion attenuation relationship. The larger ground-motion uncertainties in both the median and its standard deviation derive a higher hazard in the New Madrid area. The larger standard deviation, in particular, derives a higher hazard at low annual probabilities of exceedance. Moreover, Frankel also equated the hazard (annual probability of exceedance) defined in probabilistic seismic hazard analysis (PSHA) with the average recurrence rate defined in engineering risk analysis. There is a fundamental difference in defining risk between PSHA and the engineering risk analysis. The risks—10%, 5%, and 2% probability of exceedance (PE) in 50 years—defined in PSHA are not interchangeable with those defined in engineering risk analysis. Hence, though Frankel's stated intention was to improve understanding of PSHA among nonpractitioners, the paper may cause more confusion and problems for users of the U.S. Geological Survey national seismic hazard maps.

Seismological Research Letters↗