Search USGS⌕ Search

Geology topics

Edward H. Field

Publications and source records attributed to Edward H. Field.

At least 19 recordsLinked to original sources

Testing characteristic magnitude distributions in modern PSHA models

The characteristic magnitude distribution hypothesis predicts a higher rate of large earthquakes than a Gutenberg–Richter extrapolation of the small‐earthquake rate would imply. Characteristic magnitude distributions have been commonly applied to faults in probabilistic seismic hazard analysis (PSHA), and in modern models they can emerge from the way short‐term seismicity constraints are combined with long‐term geologic and geodetic constraints. We test the characteristic magnitude distribution hypothesis by comparing the fault‐based magnitude distributions from the 2023 update to the National Seismic Hazard Model (NSHM23) in the Western United States with observed seismicity over the past 93 yr. We find that observed magnitude distributions fall outside the model‐predicted confidence bounds in regions where NSHM23 produces characteristic magnitude distributions: in these regions, the model predicts higher rates of large earthquakes than are observed. An analysis of the earlier California model (Uniform California Earthquake Rupture Forecast, version 3) also reveals discrepancies between the modeled and observed magnitude distributions. In addition, we find that observed magnitude distributions near modeled faults are not significantly different from those in background regions. These results challenge the prevalence of characteristic magnitude distributions in fault‐based seismic hazard models and call for a reassessment of how disparate data sets are integrated in PSHA.

western United States↗

2025 USGS National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands: Overview of model and hazard results

The U.S. Geological Survey recently updated the National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). The first version of the PRVI NSHM was released in 2003, and therefore this 2025 update includes over 20 years of new geologic, geophysical, and engineering data, methods, and models. Updates follow similar efforts performed in the recent 2023 50-state NSHM. However, this is the first NSHM in which we: (1) apply an inversion methodology to subduction interface fault sources in the earthquake rupture forecast (ERF) model; (2) develop scaled backbone median ground-motion models and independent aleatory variability models that are applied in the ground-motion characterization (GMC) model; and (3) calculate epistemic uncertainty related to alternative scenarios in the ERF and GMC models for all grid points in the study region. Long-term time-independent mean hazard calculations were performed for peak ground acceleration and 5%-damped pseudospectral acceleration at 21 spectral periods from 0.01- to 10.0-s, for eight National Earthquake Hazards Reduction Program site conditions ranging from V S30 = 150 to 1500 m/s, and for 2%, 5%, and 10% in 50-year probabilities of exceedance (return periods of 2475, 975, and 475 years, respectively). Epistemic uncertainty, in the form of selected percentiles, is also provided for a suite of test sites and all grid points in the study region for limited periods, site conditions, and probabilities of exceedance. Selected results, including comparisons with the 2003 PRVI NSHM, are shown and discussed for selected periods, site conditions, and probabilities of exceedance. When comparing the 2025 PRVI NSHM with the 2003 PRVI NSHM, hazard is generally higher at shorter periods and lower at longer periods, as a result of updates in both ERF and GMC models. The 2025 PRVI NSHM is applicable for return periods greater than ∼475 or less than ∼10,000 years.

Puerto Rico, U.S. Virgin Islands↗

The U.S. Geological Survey 2025 Puerto Rico and U.S. Virgin Islands time-independent earthquake rupture forecast

We present the 2025 U.S. Geological Survey Puerto Rico and U.S. Virgin Islands (PRVI) time‐independent earthquake rupture forecast (ERF), developed for the 2025 update to the National Seismic Hazard Model (NSHM) for PRVI. The updated ERF improves upon a prior model from 2003, including an expanded fault inventory with slip‐rate estimates, updated seismicity catalogs, and refined subduction zone geometries and deformation models. It applies the fault‐system inversion methodology to solve for rates of ruptures on modeled faults, adapted from the 2023 NSHM (NSHM23) for the western United States, including the first application of the inversion to model rates on a U.S. subduction interface. Off‐fault and intraslab seismicity are constrained by observed seismicity and use updated methods developed for NSHM23. Uncertainties in model components are substantial, and the ERF represents epistemic uncertainties through a comprehensive logic tree consisting of 1.7 billion logic‐tree branches combined across all sources.

Puerto Rico, U.S. Virgin Islands↗

An entropic explanation for Gutenberg-Richter scaling

We develop a simple explanation for Gutenberg-Richter (G-R) size scaling of earthquakes on a single fault. We discretize the fault and consider all possible contiguous ruptures at that level of discretization. In this static model, we assume that slip scales with rupture length, and that the rupture rates at each point along the fault are consistent with an a priori long-term slip rate. These simple assumptions define an (under-determined) non-negative least-squares inverse problem. Each solution to this inverse problem is a set of earthquake rates that matches the slip-rate constraint. We use a Markov Chain Monte Carlo (MCMC) algorithm to uniformly sample the solution space assuming constant slip rates along the fault. At finer discretizations, deviations from G-R behavior decrease, which is consistent with an entropic pressure towards G-R solutions. When the fault is discretized into 10 or more segments, random solutions found by the MCMC algorithm have G-R size scaling, even though there are trivial solutions that, for example, have earthquakes of only one size. This is because there are simply far more solutions that have G-R scaling; as the problem size increases, the strong degeneracy of GR solutions results in other solutions becoming improbably rare. Also, the entropically favored G-R distribution has a b -value of approximately 1, which agrees with measured b -values in real earthquake catalogs.

JGR Solid Earth↗

A scientific vision and roadmap for earthquake rupture forecast developments, a USGS perspective

We articulate a scientific vision and roadmap for the development of improved Earthquake Rupture Forecast models, which are one of the two main modeling components used in modern seismic hazard and risk analysis. One primary future objective is to provide fully time-dependent models that include both elastic rebound and spatiotemporal clustering nationwide, which is particularly important for shorter-term hazard and risk considerations (e.g., earthquake insurance products). We also discuss the importance and perennial challenges associated with quantifying epistemic uncertainties, including those associated with deformation-model slip rates, un-quantified sampling errors with respect to off-fault seismicity, and any spatial covariances. The need for more physics-based approaches is also emphasized, as is the benefit of adding model valuation (quantifying usefulness) to our verification and validation protocols. Given the multidisciplinary and system-level nature of this activity, modular design is critical. Future updates will also draw from best-available science by both the United States Geological Survey and the external community. The primary goal of this paper is to highlight plans that guide research and facilitate community engagement with model development, especially with respect to lowering the entry barrier for early career scientists and engineers. The paper is written so readers can focus on the sections that interest them most (see table of contents), with the Introduction and Discussion providing a stand-alone overview and summary.

Bulletin of the Seismological Society of America↗

Risk implications of Poisson assumptions and declustering inferred from a fully time-dependent earthquake forecast

We use the Third Uniform California Earthquake Rupture Forecast Epidemic Type Aftershock Sequence model, which is fully time-dependent in terms of including spatiotemporal clustering, to evaluate the effects of the Poisson assumption and declustering algorithms on statewide loss exceedance curves. The model is simulation based, meaning it produces synthetic catalogs that exhibit realistic behavior with respect to aftershocks and multi-fault earthquakes. A Poisson version of the model was constructed by randomizing event times, and the influence of two declustering algorithms was examined as well. We demonstrate that the probability of one-or-more loss exceedances (occurrence exceedance probability) is greater for the Poisson model because it has fewer seismically quiet time windows. The discrepancy between dollar loss estimates with a given exceedance probability is up to a factor of 32% but varies depending on the loss threshold (the x-axis value) and the forecast duration (we examined a range between 24 h and 50 years, with the discrepancy for the latter being negligible). We discuss how the one-or-more loss exceedance metric is questionable because it ignores all but the maximum loss experienced in each timeframe. An alternative metric based on total aggregate loss in each time window (aggregate exceedance probability) was therefore also examined, for which the Poisson model again implies higher risk at intermediate losses but lower risk at higher losses (because large, triggered events now contribute to total aggregate losses for the fully time-dependent model). We also argue that declustering is not a scientifically justifiable way to deal with full time dependence, in agreement with a chorus from other recent studies. It is difficult to draw generally applicable conclusions from our study, in part because application specific details will likely be important, but our results highlight how full time dependence can be reckoned with once authoritative forecast models are made available.

California↗

U.S. Geological Survey Earthquake Hazards Program decadal science strategy, 2024–33

Executive Summary Earthquakes represent one of our Nation’s most significant and costly natural hazards, with estimated annual loses from earthquakes close to $15 billion in 2023. Over the past two centuries, 37 U.S. States have experienced an earthquake exceeding a magnitude of 5, and 50 percent of States have a significant potential for future damaging shaking; these statistics speak to the need for nationwide interest and investment in earthquake hazard characterization and risk reduction. Authorized under the Earthquake Hazards Reduction Authorization Act, the U.S. Geological Survey (USGS) Earthquake Hazards Program (EHP) provides the scientific information, situational awareness, and knowledge necessary to reduce deaths, injuries, and economic losses from earthquakes and earthquake-induced tsunamis, landslides, and soil liquefaction. The EHP supports activities in three focused topical areas: (1) earthquake monitoring, (2) hazard assessment, and (3) applied research, using the results of each—and the coordination among them—to further support risk translation and communication in regions at risk nationwide. For earthquake monitoring, the Advanced National Seismic System (ANSS), a cooperative effort of USGS networks, university partner regional seismic networks, and real-time geodetic networks, collects and analyzes data on earthquakes; issues timely, reliable notifications of their occurrence and impacts; and provides data for earthquake research, hazard, and risk assessment as a foundation for building an earthquake-resilient Nation. The USGS-operated ShakeAlert Earthquake Early Warning system is a recent addition to EHP’s ANSS infrastructure. In the realm of earthquake hazard assessment, the EHP contributes to earthquake risk mitigation strategies by developing the National Seismic Hazard Model and maps, and other related products, that describe the likelihood and potential effects of earthquakes nationwide, especially in the urban areas of highest risk. The EHP also conducts research on the causes, characteristics, and effects of earthquakes and prioritizes work that directly increases the accuracy and precision of earthquake hazards assessments, earthquake forecasts, and earthquake monitoring and situational-awareness products and that supports the Nation’s earthquake mitigation practices. Bridging the EHP’s efforts across research, hazard assessments, and earthquake monitoring is a broad and comprehensive collection of earthquake information products, including the National Seismic Hazard Model, ShakeAlert, and other products describing impact, such as ShakeMap and PAGER (Prompt Assessment of Global Earthquakes for Response), which have been developed and integrated into EHP’s real-time monitoring systems. EHP funds external partners to carry out many important collaborative activities through an active external grants program—one of the largest in the USGS—and through cooperative agreements with other partners such as the university-operated regional seismic networks, funded as part of the ANSS. To continue its support of earthquake hazard characterization and risk reduction, the EHP aims to strengthen its foundational products and practices while positioning itself to respond to the evolving needs of the Nation and follow best practices of the scientific community. This document describes a strategy for the program to ensure it can meet these demands. The foundational priorities outlined in this strategy represent those activities that remain critical to the core functionality of the program and those that can be supported under current fiscal year 2024-level appropriations. Priorities described as aspirational are important for future growth, and to maintain the program’s position as a leading global resource in earthquake science, but would require increases in appropriated funding to be fully realized. Across the program’s portfolio of activities, several major themes have been identified as the most critical activities to advance EHP science over the coming decade. Together, these activities provide the framework necessary to integrate critical hazard characterization and risk reduction activities across the program. They provide the structure for research to advance the understanding of where, when, and why earthquakes occur and how we can use improved knowledge to drive short-term and actionable forecasts of seismic activity. They expand the usefulness of critical earthquake products and advance the sophistication of those products to keep pace with the rapidly evolving needs of an ever-expanding user base while maintaining the position of the USGS as a global leader in earthquake science. Focus on system-level science. Establish an automated earthquake-processing pipeline. Enhance the accuracy and reliability of the ShakeAlert Earthquake Early Warning system and plan for extension to other regions. Implement time-dependent earthquake forecasting. Develop physically realistic models. Expand computational capacity. This science strategy is organized into three primary sections. The first section provides an overview of the EHP and its budget, governance, and program council. Readers familiar with the program may wish to focus on the second section, which describes the core of the science strategy, including priorities across each of the EHP’s major program activities in monitoring, hazard assessment, and targeted research. The third section outlines science priorities that cut across program activities, including those involving collaborations external to the EHP.

American Samoa, Guam, Northern Marianas Islands, U↗

Trimming the UCERF3-TD logic tree: Model order reduction for an earthquake rupture forecast considering loss exceedance

The Uniform California Earthquake Rupture Forecast version 3-Time Dependent depicts California’s seismic faults and their activity. Its logic tree has 5760 leaves. Considering 30 more model combinations related to ground motion produces 172,800 distinct models representing so-called epistemic uncertainties. To calculate risk to a portfolio of buildings, one also considers millions of earthquakes and spatially correlated ground-motion variability. We offer a tree-trimming technique that retains the probability distribution of portfolio loss and identifies the leading sources of uncertainty for further study. We applied it to a California statewide building portfolio and various levels of nonexceedance probability between one in 100 and one in 2500. We trimmed the logic tree from 172,800 leaves to as few as 15. The result: a supercomputer that would otherwise run 24 h to estimate the distribution of one-in-250-year loss can calculate it in moments with the reduced-order model. Others can use the reduced-order model to calculate risk to different California portfolios, and scientists can prioritize study to reduce the remaining epistemic uncertainty.

Earthquake Spectra↗

The 2023 US 50-State National Seismic Hazard Model: Overview and implications

The US National Seismic Hazard Model (NSHM) was updated in 2023 for all 50 states using new science on seismicity, fault ruptures, ground motions, and probabilistic techniques to produce a standard of practice for public policy and other engineering applications (defined for return periods greater than ∼475 or less than ∼10,000 years). Changes in 2023 time-independent seismic hazard (both increases and decreases compared to previous NSHMs) are substantial because the new model considers more data and updated earthquake rupture forecasts and ground-motion components. In developing the 2023 model, we tried to apply best available or applicable science based on advice of co-authors, more than 50 reviewers, and hundreds of hazard scientists and end-users, who attended public workshops and provided technical inputs. The hazard assessment incorporates new catalogs, declustering algorithms, gridded seismicity models, magnitude-scaling equations, fault-based structural and deformation models, multi-fault earthquake rupture forecast models, semi-empirical and simulation-based ground-motion models, and site amplification models conditioned on shear-wave velocities of the upper 30 m of soil and deeper sedimentary basin structures. Seismic hazard calculations yield hazard curves at hundreds of thousands of sites, ground-motion maps, uniform-hazard response spectra, and disaggregations developed for pseudo-spectral accelerations at 21 oscillator periods and two peak parameters, Modified Mercalli Intensity, and 8 site classes required by building codes and other public policy applications. Tests show the new model is consistent with past ShakeMap intensity observations. Sensitivity and uncertainty assessments ensure resulting ground motions are compatible with known hazard information and highlight the range and causes of variability in ground motions. We produce several impact products including building seismic design criteria, intensity maps, planning scenarios, and engineering risk assessments showing the potential physical and social impacts. These applications provide a basis for assessing, planning, and mitigating the effects of future earthquakes.

Earthquake Spectra↗

The USGS 2023 Conterminous U.S. time‐independent earthquake rupture forecast

We present the 2023 U.S. Geological Survey time‐independent earthquake rupture forecast for the conterminous United States, which gives authoritative estimates of the magnitude, location, and time‐averaged frequency of potentially damaging earthquakes throughout the region. In addition to updating virtually all model components, a major focus has been to provide a better representation of epistemic uncertainties. For example, we have improved the representation of multifault ruptures, both in terms of allowing more and less fault connectivity than in the previous models, and in sweeping over a broader range of viable models. An unprecedented level of diagnostic information has been provided for assessing the model, and the development was overseen by a 19‐member participatory review panel. Although we believe the new model embodies significant improvements and represents the best available science, we also discuss potential model limitations, including the applicability of logic tree branch weights with respect different types of hazard and risk metrics. Future improvements are also discussed, with deformation model enhancements being particularly worthy of pursuit, as well as better representation of sampling errors in the gridded seismicity components. We also plan to add time‐dependent components, and assess implications with a wider range of hazard and risk metrics.

Bulletin of the Seismological Society of America↗

A comprehensive fault system inversion approach: Methods and application to NSHM23

We present updated inversion‐based fault‐system solutions for the 2023 update to the National Seismic Hazard Model (NSHM23), standardizing earthquake rate model calculations on crustal faults across the western United States. We build upon the inversion methodology used in the Third Uniform California Earthquake Rupture Forecast (UCERF3) to solve for time‐independent rates of earthquakes in an interconnected fault system. The updated model explicitly maps out a wide range of fault recurrence and segmentation behavior (epistemic uncertainty), more completely exploring the solution space of viable models beyond those of UCERF3. We also improve the simulated annealing implementation, greatly increasing computational efficiency (and thus inversion convergence), and introduce an adaptive constraint weight calculation algorithm that helps to mediate between competing constraints. Hazard calculations show that ingredient changes (especially fault and deformation models) are the primary driver of hazard changes between NSHM23 and UCERF3. Updates to the inversion methodology are also consequential near faults in which the slip rate in UCERF3 was poorly fit or was satisfied primarily using large multifault ruptures that are now restricted by explicit b ‐value and segmentation constraints.

Bulletin of the Seismological Society of America↗

Western U.S. deformation models for the 2023 update to the U.S. National Seismic Hazard Model

This report describes geodetic and geologic information used to constrain deformation models of the 2023 update to the National Seismic Hazard Model (NSHM), a set of deformation models to interpret these data, and their implications for earthquake rates in the western United States. Recent updates provide a much larger data set of Global Positioning System crustal velocities than used in the 2014 NSHM, as well as hundreds of new faults considered as active sources for the 2023 NSHM. These data are interpreted by four geodetic models of deformation that estimate fault slip rates and their uncertainties together with off‐fault moment release rates. Key innovations in the 2023 NSHM relative to past practice include (1) the addition of two new (in addition to two existing) deformation models, (2) the revision and expansion of the geologic slip rate database, (3) accounting for fault creep through development of a creep‐rate model that is employed by the four deformation models, and (4) accounting for time‐dependent earthquake‐cycle effects through development of viscoelastic models of the earthquake cycle along the San Andreas fault and the Cascadia subduction zone. The effort includes development of a geologic deformation model that complements the four geodetic models. The current deformation models provide a new assessment of outstanding discrepancies between geologic and geodetic slip rates, at the same time highlighting the need for both geologic and geodetic slip rates to robustly inform the earthquake rate model.

Seismological Research Letters↗

Simplifying complex fault data for systems-level analysis: Earthquake geology inputs for U.S. NSHM 2023

As part of the U.S. National Seismic Hazard Model (NSHM) update planned for 2023, two databases were prepared to more completely represent Quaternary-active faulting across the western United States: the NSHM23 fault sections database (FSD) and earthquake geology database (EQGeoDB). In prior iterations of NSHM, fault sections were included only if a field-measurement-derived slip rate was estimated along a given fault. By expanding this inclusion criteria, we were able to assess a larger set of faults for use in NSHM23. The USGS Quaternary Fault and Fold Database served as a guide for assessing possible additions to the NSHM23 FSD. Reevaluating available data from published sources yielded an increase of fault sections from ~650 faults in NSHM18 to ~1,000 faults proposed for use in NSHM23. EQGeoDB, a companion dataset linked to NSHM23 FSD, contains geologic slip rate estimates for fault sections included in FSD. Together, these databases serve as common input data used in deformation modeling, earthquake rupture forecasting, and additional downstream uses in NSHM development.

Scientific Data↗

Enumerating plausible multifault ruptures in complex fault systems with physical constraints

We propose a new model for determining the set of plausible multifault ruptures in an interconnected fault system. We improve upon the rules used in the Third Uniform California Earthquake Rupture Forecast (UCERF3) to increase connectivity and the physical consistency of ruptures. We replace UCERF3’s simple azimuth change rules with new Coulomb favorability metrics and increase the maximum jump distance to 15 km. Although the UCERF3 rules were appropriate for faults with similar rakes, the Coulomb calculations used here inherently encode preferred orientations between faults with different rakes. Our new rules are designed to be insensitive to discretization details and are generally more permissive than their UCERF3 counterparts; they allow more than twice the connectivity compared to UCERF3, yet heavily penalize long ruptures that take multiple improbable jumps. The set of all possible multifault ruptures in the California fault system is near-infinite, but our model produces a tractable set of 326,707 ruptures (a modest 29% increase over UCERF3, despite the greatly increased connectivity). Inclusion in the rupture set does not dictate that a rupture receives a significant rate in the final model; rupture rates are subsequently determined by data constraints used in an inversion. We describe the rupture building algorithm and its components in detail and provide comparisons with ruptures generated by a physics-based multicycle earthquake simulator. We find that greater than twice as many ruptures generated by the simulator violate the UCERF3 rules than violate our proposed model.

California↗

STEPS: Slip time earthquake path simulations applied to the San Andreas and Toe Jam Hill Faults to redefine geologic slip rate uncertainty

Geologic slip rates are a time-averaged measurement of fault displacement calculated over hundreds to million-year time scales and are a primary input for probabilistic seismic hazard analyses, which forecast expected ground shaking in future earthquakes. Despite their utility for seismic hazard calculations, longer-term geologic slip rates represent a time-averaged measure of the tempo of strain release and do not measure variability across earthquake cycles. We have developed a numerical approach called STEPS (Slip Time Earthquake Path Simulations), which is built upon field-based observations and explicitly incorporates realistic variations in displacement per event and variability in the recurrence interval between earthquakes. The STEPS approach, which simulates strain release through time, relies on representing earthquake cycles as stairsteps, rather than straight-line paths, connecting per earthquake time-displacement coordinates. We simulate earthquake histories based on these input constraints using two examples: the Carrizo section of the San Andreas fault and the Toe Jam Hill fault of the Seattle fault zone. We find that modeled slip rate distributions agree with slip rates reported for the sites of interest by the original investigators, while providing a slip rate distribution that reflects the variability of earthquake frequency and displacement. The STEPS approach provides an estimate of fault slip rate uncertainty based on a synthetic suite of plausible time-displacement paths resulting from individual earthquakes, rather than measurement uncertainties associated with offset features. When considering this simulated earthquake behavior between measurements, the uncertainty associated with earthquake paths is greater than that calculated by the long-term rate.

California, Washington↗

The seismic hazard implications of declustering and poisson assumptions inferred from a fully time‐dependent model

We use the Third Uniform California Earthquake Rupture Forecast (UCERF3) epidemic‐type aftershock sequence (ETAS) model (UCERF3‐ETAS) to evaluate the effects of declustering and Poisson assumptions on seismic hazard estimates. Although declustering is necessary to infer the long‐term spatial distribution of earthquake rates, the question is whether it is also necessary to honor the Poisson assumption in classic probabilistic seismic hazard assessment. We use 500,000 yr, M ≥ 2.5 synthetic catalogs to address this question, for which UCERF3‐ETAS exhibits realistic spatiotemporal clustering effects (e.g., aftershocks). We find that Gardner and Knopoff (1974) declustering, used in the U.S. Geological Survey seismic hazard models, lowers 2% in 50 yr and risk‐targeted ground‐motion hazard metrics by about 4% on average (compared with the full time‐dependent [TD] model), with the reduction being 5% at 40% in 50 yr ground motions. Keeping all earthquakes and treating them as a Poisson process increases these same hazard metrics by about 3%–12%, on average, due to the removal of relatively quiet time periods in the full TD model. In the interest of model simplification, bias minimization, and consideration of the probabilities of multiple exceedances, we agree with others ( Marzocchi and Taroni, 2014 ) that we are better off keeping aftershocks and treating them as a Poisson process rather than removing them from hazard consideration via declustering. Honoring the true time dependence, however, will likely be important for other hazard and risk metrics, and this study further exemplifies how this can now be evaluated more extensively.

California↗

Improvements to the Third Uniform California Earthquake Rupture Forecast ETAS Model (UCERF3‐ETAS)

We describe recent improvements to the Third Uniform California Earthquake Rupture Forecast ETAS Model (UCERF3‐ETAS), which continues to represent our most advanced and complete earthquake forecast in terms of relaxing segmentation assumptions and representing multifault ruptures, elastic‐rebound effects, and spatiotemporal clustering (the latter to represent aftershocks and otherwise triggered events). The two main improvements include adding aleatory variability in aftershock productivity and the option to represent off‐fault events with finite‐rupture surfaces. We also summarize the studies that led to these modifications, and reflect on how past and future uses of the model can improve our understanding of earthquake processes and the hazards and risks they pose.

California↗