USGS Science⌕ Search

SEARCH · USGS Science

Results for “Bulletin of the Geological Society of America”

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 775 records · Page 43Linked to original sources
USGS2025

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↗
USGS2016

Uncertainty in V s30 -based site response

Methods that account for site response range in complexity from simple linear categorical adjustment factors to sophisticated nonlinear constitutive models. Seismic‐hazard analysis usually relies on ground‐motion prediction equations (GMPEs); within this framework site response is modeled statistically with simplified site parameters that include the time‐averaged shear‐wave velocity to 30 m ( V S 30 ) and basin depth parameters. Because V S 30 is not known in most locations, it must be interpolated or inferred through secondary information such as geology or topography. In this article, we analyze a subset of stations for which V S 30 has been measured to address effects of V S 30 proxies on the uncertainty in the ground motions as modeled by GMPEs. The stations we analyze also include multiple recordings, which allow us to compute the repeatable site effects (or empirical amplification factors [EAFs]) from the ground motions. Although all methods exhibit similar bias, the proxy methods only reduce the ground‐motion standard deviations at long periods when compared to GMPEs without a site term, whereas measured V S 30 values reduce the standard deviations at all periods. The standard deviation of the ground motions are much lower when the EAFs are used, indicating that future refinements of the site term in GMPEs have the potential to substantially reduce the overall uncertainty in the prediction of ground motions by GMPEs.

Bulletin of the Seismological Society of America↗
USGS2017

The 2008 Wells, Nevada earthquake sequence: Source constraints using calibrated multiple event relocation and InSAR

The 2008 Wells, NV earthquake represents the largest domestic event in the conterminous U.S. outside of California since the October 1983 Borah Peak earthquake in southern Idaho. We present an improved catalog, magnitude complete to 1.6, of the foreshock-aftershock sequence, supplementing the current U.S. Geological Survey (USGS) Preliminary Determination of Epicenters (PDE) catalog with 1,928 well-located events. In order to create this catalog, both subspace and kurtosis detectors are used to obtain an initial set of earthquakes and associated locations. The latter are then calibrated through the implementation of the hypocentroidal decomposition method and relocated using the BayesLoc relocation technique. We additionally perform a finite fault slip analysis of the mainshock using InSAR observations. By combining the relocated sequence with the finite fault analysis, we show that the aftershocks occur primarily updip and along the southwestern edge of the zone of maximum slip. The aftershock locations illuminate areas of post-mainshock strain increase; aftershock depths, ranging from 5 to 16 km, are consistent with InSAR imaging, which shows that the Wells earthquake was a buried source with no observable near-surface offset.

Nevada↗
USGS1988

Effects of three-dimensional velocity structure on the seismicity of the 1984 Morgan Hill, California, aftershock sequence

A three-dimensional velocity model for the area surrounding the 24 April 1984 Morgan Hill earthquake has been developed by simultaneously inverting local earthquake and refraction arrival-time data. This velocity model corresponds well to the surface geology of the region, predominantly showing a low-velocity region associated with the sedimentary sequence to the southwest of the Madrone Springs fault. This low-velocity sequence extends to a depth of approximately 5 km. Previously constructed geologic cross-sections for this region had suggested that these sediments extended to a depth of less than 2 km. Relocation of 2,179 hypocenters with the three-dimensional velocity model did not show great differences when compared to locations done with a one-dimensional velocity model with station corrections. In both cases, the Calaveras fault is steep but has a noticeable dip to the northeast and the seismicity does not intersect the surface along the trace of the Calaveras fault in the southern half of the rupture zone. The greatest differences between the two sets of locations was that the locations using the three-dimensional velocity model were deeper, and fewer earthquakes were located within the low-velocity region. The focal mechanisms were also determined for 946 earthquakes using both the one-dimensional and three-dimensional earth models. Both earth models yield similar focal mechanisms for these earthquakes.

California↗
USGS2026

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↗
USGS2021

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↗
USGS2020

Human behavioral response in the Ridgecrest earthquakes: Assessing immediate actions based on data from “Did You Feel It?”

Human behavioral response to earthquake ground motion has long been a subject of multidisciplinary interest and research. In most versions of seismic intensity scales, human perceptions and behavior are one component of the assignment of intensity. Public health research has shown that actions taken during earthquakes have a significant impact on the incidence of injury or the maintenance of safety. Based on this research, emergency managers and organizations promoting emergency preparedness have advocated strategies such as drop, cover, and hold on (DCHO) and promoted this safety measure through public education and annual drills. The “Did You Feel It?” (DYFI) mapping system (see Data and Resources) based on an online questionnaire developed and maintained by the U.S. Geological Survey has provided opportunities for those who have experienced an earthquake to report this experience worldwide since 2004. The DYFI questionnaire, although designed to assign intensity, also contains questions regarding the behavior in which one has engaged during the earthquake. The questionnaire includes other important information that may elucidate behavioral response to earthquakes, including assigned intensity, emotional reaction, and whether damage occurred at the location where the earthquake was experienced. The very large number of people who completed DYFI questionnaires following the July 2019 Ridgecrest, California, earthquakes provides a robust dataset for analysis and suggests that as intensity and levels of fear increase, behavior becomes more active in terms of physical movement to locations of presumed safety. Among active responses including DCHO, going to a doorway, and running outside, DCHO was the least likely to be implemented. The study provides possible explanations for low participation in DCHO despite active campaigns to promote this strategy.

California↗
USGS2000

A time-dependent probabilistic seismic-hazard model for California

For the purpose of sensitivity testing and illuminating nonconsensus components of time-dependent models, the California Department of Conservation, Division of Mines and Geology (CDMG) has assembled a time-dependent version of its statewide probabilistic seismic hazard (PSH) model for California. The model incorporates available consensus information from within the earth-science community, except for a few faults or fault segments where consensus information is not available. For these latter faults, published information has been incorporated into the model. As in the 1996 CDMG/U.S. Geological Survey (USGS) model, the time-dependent models incorporate three multisegment ruptures: a 1906, an 1857, and a southern San Andreas earthquake. Sensitivity tests are presented to show the effect on hazard and expected damage estimates of (1) intrinsic (aleatory) sigma, (2) multisegment (cascade) vs. independent segment (no cascade) ruptures, and (3) time-dependence vs. time-independence. Results indicate that (1) differences in hazard and expected damage estimates between time-dependent and independent models increase with decreasing intrinsic sigma, (2) differences in hazard and expected damage estimates between full cascading and not cascading are insensitive to intrinsic sigma, (3) differences in hazard increase with increasing return period (decreasing probability of occurrence), and (4) differences in moment-rate budgets increase with decreasing intrinsic sigma and with the degree of cascading, but are within the expected uncertainty in PSH time-dependent modeling and do not always significantly affect hazard and expected damage estimates.

California↗
USGS2021

Diffuse deformation and surface faulting distribution from sub-metric image correlation along the 2019 Ridgecrest ruptures (California, USA)

The 2019 M w "> M w 6.4 and 7.1 Ridgecrest, California, earthquake sequence (July 2019) ruptured consecutively a system of high‐angle strike‐slip cross faults (northeast‐ and northwest‐trending) within 34 hr. The complex rupture mechanism was illuminated by seismological and geodetic data, bringing forward the issue of the interdependency of the two fault systems both at depth and at the surface, and of its effect on the final surface displacement pattern. Here, we use high‐resolution (WorldView and Pleiades) optical satellite image correlation to measure the near‐fault horizontal and vertical surface displacement fields at 0.5 m ground resolution for the two earthquakes. We point out significant differences with previous geodetic‐ and geologic‐based measurements, and document the essential role of distributed faulting and diffuse deformation in producing the observed surface displacement patterns. We derive strain fields from the horizontal displacement maps, and highlight the predominant role of rotation and shear strain in the surface rupture process. We discuss the segmentation of the rupture based on the fault geometry and along‐strike slip variations. We also image several northeast‐trending faults with similar orientation to the deeply embedded shear fabric identified in aftershock studies, and show that these cross faults are present all along the rupture, including at a scale <100 m. Finally, we compare our results to kinematic slip inversions, and show that the surface diffuse deformation is primarily associated with areas of shallow slip deficit; however, this diffuse deformation cannot be explained using elastic modeling. We conclude that inelastic processes play an important role in contributing to the total surface deformation associated with the 2019 Ridgecrest sequence.

California↗
USGS2016

Three ingredients for Improved global aftershock forecasts: Tectonic region, time-dependent catalog incompleteness, and inter-sequence variability

Following a large earthquake, seismic hazard can be orders of magnitude higher than the long‐term average as a result of aftershock triggering. Because of this heightened hazard, emergency managers and the public demand rapid, authoritative, and reliable aftershock forecasts. In the past, U.S. Geological Survey (USGS) aftershock forecasts following large global earthquakes have been released on an ad hoc basis with inconsistent methods, and in some cases aftershock parameters adapted from California. To remedy this, the USGS is currently developing an automated aftershock product based on the Reasenberg and Jones (1989) method that will generate more accurate forecasts. To better capture spatial variations in aftershock productivity and decay, we estimate regional aftershock parameters for sequences within the García et al. (2012) tectonic regions. We find that regional variations for mean aftershock productivity reach almost a factor of 10. We also develop a method to account for the time‐dependent magnitude of completeness following large events in the catalog. In addition to estimating average sequence parameters within regions, we develop an inverse method to estimate the intersequence parameter variability. This allows for a more complete quantification of the forecast uncertainties and Bayesian updating of the forecast as sequence‐specific information becomes available.

California↗
USGS2026

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↗
USGS1971

Geological aspects of the May 31, 1970, Perú earthquake

Geological effects of the destructive May 31, 1970, Perú earthquake (Ms = 7.7) extended over roughly 65,000 km2 of west-central Perú. Earthquake-triggered slope failures of all types that occurred throughout the mountainous parts of the region extensively damaged transportation routes and irrigation canals and temporarily dammed some rivers and lakes. The geologically most important and spectacular of these, a cataclysmic debris avalanche from between 5,500 and 6,400 m altitude on the north peak of Huascarán Mountain, buried the city of Yungay and part of Ranrahirca (elevation about 2,500 m) with a loss of more than about 18,000 lives. The debris avalanche involved 50-100 million m3 of rock, ice, snow, and soil that traveled 14.5 km from the source to Yungay at an average velocity of between 280 and 335 km/hr. Exceptionally rapid movement of the avalanche is indicated by eyewitness accounts, by topographic irregularities as high as 140 m that were overridden, and locally by boulders weighing several tons that were hurled as much as 1,000 m beyond its margins. A pulse of muddy water from the debris avalanche that swept down the Río Santa 160 km to the sea inundated farms and small settlements, buried highway and railroad routes, and destroyed the diversion dam and access bridge to a major hydroelectric plant.

Bulletin of the Seismological Society of America↗
USGS1996

Application of ground-penetrating radar to investigation of near-surface fault properties in the San Francisco Bay region

In many geologic environments, ground-penetrating radar (GPR) provides high-resolution images of near-surface Earth structure. GPR data collection is nondestructive and very economical. The scale of features detected by GPR lies between those imaged by high-resolution seismic reflection surveys and those exposed in trenches and is therefore potentially complementary to traditional techniques for fault location and mapping. Sixty-two GPR profiles were collected at 12 sites in the San Francisco Bay region. Results show that GPR data correlate with large-scale features in existing trench observations, can be used to locate faults where they are buried or where their positions are not well known, and can identify previously unknown fault segments. The best data acquired were on a profile across the San Andreas fault, traversing Pleistocene terrace deposits south of Olema in Marin County; this profile shows a complicated multi-branched fault system from the ground surface down to about 40 m, the maximum depth for which data were recorded.

California↗
USGS2013

A scenario study of seismically induced landsliding in Seattle using broadband synthetic seismograms

We demonstrate the value of utilizing broadband synthetic seismograms to assess regional seismically induced landslide hazard. Focusing on a case study of an M w 7.0 Seattle fault earthquake in Seattle, Washington, we computed broadband synthetic seismograms that account for rupture directivity and 3D basin amplification. We then adjusted the computed motions on a fine grid for 1D amplifications based on the site response of typical geologic profiles in Seattle and used these time‐series ground motions to trigger shallow landsliding using the Newmark method. The inclusion of these effects was critical in determining the extent of landsliding triggered. We found that for inertially triggered slope failures modeled by the Newmark method, the ground motions used to simulate landsliding must have broadband frequency content in order to capture the full slope displacement. We applied commonly used simpler methods based on ground‐motion prediction equations for the same scenario and found that they predicted far fewer landslides if only the mean values were used, but far more at the maximum range of the uncertainties, highlighting the danger of using just the mean values for such methods. Our results indicate that landsliding triggered by a large Seattle fault earthquake will be extensive and potentially devastating, causing direct losses and impeding recovery. The high impact of landsliding predicted by this simulation shows that this secondary effect of earthquakes should be studied with as much vigor as other earthquake effects.

Washington↗
USGS2000

Discrepancy between earthquake rates implied by historic earthquakes and a consensus geologic source model for California

We examine the difference between expected earthquake rates inferred from the historical earthquake catalog and the geologic data that was used to develop the consensus seismic source characterization for the state of California [California Department of Conservation, Division of Mines and Geology (CDMG) and U.S. Geological Survey (USGS) Petersen et al., 1996; Frankel et al., 1996]. On average the historic earthquake catalog and the seismic source model both indicate about one M 6 or greater earthquake per year in the state of California. However, the overall earthquake rates of earthquakes with magnitudes (M) between 6 and 7 in this seismic source model are higher, by at least a factor of 2, than the mean historic earthquake rates for both southern and northern California. The earthquake rate discrepancy results from a seismic source model that includes earthquakes with characteristic (maximum) magnitudes that are primarily between M 6.4 and 7.1. Many of these faults are interpreted to accommodate high strain rates from geologic and geodetic data but have not ruptured in large earthquakes during historic time. Our sensitivity study indicates that the rate differences between magnitudes 6 and 7 can be reduced by adjusting the magnitude-frequency distribution of the source model to reflect more characteristic behavior, by decreasing the moment rate available for seismogenic slip along faults, by increasing the maximum magnitude of the earthquake on a fault, or by decreasing the maximum magnitude of the background seismicity. However, no single parameter can be adjusted, consistent with scientific consensus, to eliminate the earthquake rate discrepancy. Applying a combination of these parametric adjustments yields an alternative earthquake source model that is more compatible with the historic data. The 475-year return period hazard for peak ground and 1-sec spectral acceleration resulting from this alternative source model differs from the hazard resulting from the standard CDMG-USGS model by less than 10% across most of California but is higher (generally about 10% to 30%) within 20 km from some faults.

Bulletin of the Seismological Society of America↗
USGS2026

Use of precariously balanced rocks to constrain postglacial earthquake magnitudes in New England, United States

We use the age and fragility of 13 precariously balanced rocks (PBRs) in Massachusetts, New Hampshire, and southwestern Maine to evaluate the estimated probabilistic seismic hazard (PSH) from the U.S. Geological Survey’s National Seismic Hazard Model (NSHM), and also develop a map of maximum postglacial earthquake magnitudes in the region. The PBRs are erratic boulders left behind by the retreat of the Laurentide ice sheet between 13 and 18 ka, and published ice retreat maps are used to estimate their fragility ages. The fragilities and ages of the rocks appear to be consistent with the NSHM, in that the estimated earthquake magnitudes and resulting ground motions that would topple them are only expected on timescales longer than their fragility ages. This result, on the North American passive margin, contrasts with PBR studies at active plate boundaries, where the PBR constraints imply large reductions to estimated PSH. It may be that the PBRs are not old enough in the eastern United States to reveal the same PSH modeling issues as those revealed in the plate boundary settings. The study provides constraints on maximum magnitudes of postglacial earthquakes across much of Massachusetts, New Hampshire, and southwestern Maine, with the constraints limited to the ages of the PBRs.

Maine. Massachusetts, New Hampshire↗
USGS2010

High-resolution seismic reflection imaging of growth folding and shallow faults beneath the Southern Puget Lowland, Washington State

Marine seismic reflection data from southern Puget Sound, Washington, were collected to investigate the nature of shallow structures associated with the Tacoma fault zone and the Olympia structure. Growth folding and probable Holocene surface deformation were imaged within the Tacoma fault zone beneath Case and Carr Inlets. Shallow faults near potential field anomalies associated with the Olympia structure were imaged beneath Budd and Eld Inlets. Beneath Case Inlet, the Tacoma fault zone includes an &sim;350-m wide section of south-dipping strata forming the upper part of a fold (kink band) coincident with the southern edge of an uplifted shoreline terrace. An &sim;2 m change in the depth of the water bottom, onlapping postglacial sediments, and increasing stratal dips with increasing depth are consistent with late Pleistocene to Holocene postglacial growth folding above a blind fault. Geologic data across a topographic lineament on nearby land indicate recent uplift of late Holocene age. Profiles acquired in Carr Inlet 10 km to the east of Case Inlet showed late Pleistocene or Holocene faulting at one location with &sim;3 to 4 m of vertical displacement, south side up. North of this fault the data show several other disruptions and reflector terminations that could mark faults within the broad Tacoma fault zone. Seismic reflection profiles across part of the Olympia structure beneath southern Puget Sound show two apparent faults about 160 m apart having 1 to 2 m of displacement of subhorizontal bedding. Directly beneath one of these faults, a dipping reflector that may mark the base of a glacial channel shows the opposite sense of throw, suggesting strike-slip motion. Deeper seismic reflection profiles show disrupted strata beneath these faults but little apparent vertical offset, consistent with strike-slip faulting. These faults and folds indicate that the Tacoma fault and Olympia structure include active structures with probable postglacial motion.

Washington↗
USGS2002

The 1999 Izmit, Turkey, earthquake: A 3D dynamic stress transfer model of intraearthquake triggering

Before the August 1999 Izmit (Kocaeli), Turkey, earthquake, theoretical studies of earthquake ruptures and geological observations had provided estimates of how far an earthquake might jump to get to a neighboring fault. Both numerical simulations and geological observations suggested that 5 km might be the upper limit if there were no transfer faults. The Izmit earthquake appears to have followed these expectations. It did not jump across any step-over wider than 5 km and was instead stopped by a narrower step-over at its eastern end and possibly by a stress shadow caused by a historic large earthquake at its western end. Our 3D spontaneous rupture simulations of the 1999 Izmit earthquake provide two new insights: (1) the west- to east-striking fault segments of this part of the North Anatolian fault are oriented so as to be low-stress faults and (2) the easternmost segment involved in the August 1999 rupture may be dipping. An interesting feature of the Izmit earthquake is that a 5-km-long gap in surface rupture and an adjacent 25° restraining bend in the fault zone did not stop the earthquake. The latter observation is a warning that significant fault bends in strike-slip faults may not arrest future earthquakes.

Izmit↗