USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geological Society of America Bulletin”

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 919 records · Page 51Linked to original sources
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↗
USGS2008

Ground-motion modeling of the 1906 San Francisco Earthquake, part II: Ground-motion estimates for the 1906 earthquake and scenario events

We estimate the ground motions produce by the 1906 San Francisco earthquake making use of the recently developed Song et al. (2008) source model that combines the available geodetic and seismic observations and recently constructed 3D geologic and seismic velocity models. Our estimates of the ground motions for the 1906 earthquake are consistent across five ground-motion modeling groups employing different wave propagation codes and simulation domains. The simulations successfully reproduce the main features of the Boatwright and Bundock (2005) ShakeMap, but tend to over predict the intensity of shaking by 0.1-0.5 modified Mercalli intensity (MMI) units. Velocity waveforms at sites throughout the San Francisco Bay Area exhibit characteristics consistent with rupture directivity, local geologic conditions (e.g., sedimentary basins), and the large size of the event (e.g., durations of strong shaking lasting tens of seconds). We also compute ground motions for seven hypothetical scenarios rupturing the same extent of the northern San Andreas fault, considering three additional hypocenters and an additional, random distribution of slip. Rupture directivity exerts the strongest influence on the variations in shaking, although sedimentary basins do consistently contribute to the response in some locations, such as Santa Rosa, Livermore, and San Jose. These scenarios suggest that future large earthquakes on the northern San Andreas fault may subject the current San Francisco Bay urban area to stronger shaking than a repeat of the 1906 earthquake. Ruptures propagating southward towards San Francisco appear to expose more of the urban area to a given intensity level than do ruptures propagating northward.

Bulletin of the Seismological Society of America↗
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↗
USGS2004

Comparing the November 2002 Denali and November 2001 Kunlun earthquakes

Major strike-slip earthquakes recently occurred in Alaska on the central Denali fault (M 7.9) on 3 November 2002, and in Tibet on the central Kunlun fault (M 7.8) on 14 November 2001. Both earthquakes generated large surface waves with Ms [U.S. Geological Survey (USGS)] of 8.5 (Denali) and 8.0 (Kunlun). Each event occurred on an east-west-trending strike-slip fault situated near the northern boundary of an intense deformation zone that is characterized by lateral extrusion and rotation of crustal blocks. Each earthquake produced east-directed nearly unilateral ruptures that propagated 300 to 400 km. Maximum lateral surface offsets and maximum moment release occurred well beyond 100 km from the rupture initiation, with the events exhibiting by far the largest separations of USGS hypocenter and Harvard Moment Tensor Centroid (CMT) for strike-slip earthquakes in the 27-year CMT catalog. In each sequence, the largest aftershock was more than two orders of magnitude smaller than the mainshock. Regional moment release had been accelerating prior to the main shocks. The close proximity in space and time of the 1964 Prince William Sound and 2002 Denali earthquakes, relative to their rupture lengths and estimated return times, suggests that these events may be part of a recurrent cluster in the vicinity of a complex plate boundary.

Bulletin of the Seismological Society of America↗
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↗
USGS1996

Earthquake site response in Santa Cruz, California

Aftershocks of the 1989 Loma Prieta, California, earthquake are used to estimate site response in a 12-km 2 area centered on downtown Santa Cruz. A total of 258 S -wave records from 36 aftershocks recorded at 33 sites are used in a linear inversion for site-response spectra. The inversion scheme takes advantage of the redundancy of the large data set for which several aftershocks are recorded at each site. The scheme decomposes the observed spectra into source, path, and site terms. The path term is specified before the inversion. The undetermined degree of freedom in the decomposition into source and site spectra is removed by specifying the site-response factor to be approximately 1.0 at two sites on crystalline bedrock. The S -wave site responses correlate well with the surficial geology and observed damage pattern of the mainshock. The site-response spectra of the floodplain sites, which include the heavily damaged downtown area, exhibit significant peaks. The largest peaks are between 1 and 4 Hz. Five floodplain sites have amplification factors of 10 or greater. Most of the floodplain site-response spectra also have a smaller secondary peak between 6 and 8 Hz. Residential areas built on marine terraces above the floodplain experienced much less severe damage. Site-response spectra for these areas also have their largest peaks between 1 and 4 Hz, but the amplification is generally below 6. Several of these sites also have a secondary peak between 6 and 8 Hz. The response peaks seen at nearly all sites between 1 and 4 Hz are probably caused by the natural resonance of the sedimentary rock column. The higher amplifications at floodplain sites may be caused by surface waves generated at the basin margins. The secondary peak between 6 and 8 Hz at many sites may be a harmonic of the 1- to 4-Hz peaks.

California↗
USGS2022

High-resolution marine seismic imaging of the Seattle fault zone: Near surface insights into fault zone geometry, Quaternary deformation, and long-term evolution

The Seattle fault zone (SFZ) is a north‐directed thrust fault system that underlies the greater Seattle metropolitan area. Evidence of past land level changes, landslides, liquefaction, and a local tsunami indicate that this 70‐km‐long fault system can host up to M 7–7.5 earthquakes. Both the geometry and earthquake recurrence of the SFZ are debated and surveys of the shallow subsurface have not yet been incorporated into deeper crustal‐scale structural interpretations, especially where the SFZ cuts across marine portions of the Puget Lowland. Here we use a new high‐resolution marine seismic reflection dataset to image fault‐related deformation in Quaternary sediments and Tertiary bedrock throughout Puget Sound and Lake Washington. We use this perspective of shallow geology as a link between existing crustal‐scale geophysical insights into fault geometry at depth and paleoseismological observations of faulting at the surface and propose a refined structural model for the SFZ. We interpret that our new seismic reflection data in the Rich Passage area of Puget Sound images evidence of an inactive, south‐dipping strand of the SFZ, which is overprinted by Quaternary folding and slip along north‐dipping backthrusts within the hanging wall of a blind, south‐dipping fault located 6 km farther north. To explain these results, we propose that the SFZ is a normal sequence fault propagation fold that has stepped northward through time, and we show the plausibility of this model through trishear forward modeling. Growth strata and faulting imaged in Quaternary sediments in Lake Washington and Rich Passage are consistent with the spatial distribution of folding and backthrusting that occurred during an M 7–7.5 earthquake in A.D. 900–930, corroborating existing evidence that the SFZ has been active throughout the Quaternary.

Bulletin of the Seismological Society of America↗
USGS2024

Forecasting the long-term spatial distribution of earthquakes for the 2023 US National Seismic Hazard Model using gridded seismicity

Probabilistic seismic hazard analyses such as the U.S. National Seismic Hazard Model (NSHM) typically rely on declustering and spatially smoothing an earthquake catalog to estimate a long‐term time‐independent (background) seismicity rate to forecast future seismicity. In support of the U.S. Geological Survey’s (USGS) 2023 update to the NSHM, we update the methods used to develop this background or gridded seismicity model component of the NSHM. As in 2018, we use a combination of fixed and adaptive kernel Gaussian smoothing. However, we implement two additional declustering methods to account for the fact that declustering is a nonunique process. These new declustering methods result in different forecasts for the locations of future seismicity, as represented by spatial probability density functions that are later combined with a rate model to produce a full gridded earthquake rate forecast. The method updates, particularly in the separation of the spatial and rate models as well as revised regional boundaries, in some places cause substantive changes to the seismic hazard forecast compared to the previous 2018 NSHM. Additional updates to catalog processing and induced seismicity zones also contribute to changes in the gridded seismicity hazard since 2018. However, these changes are well understood and reflect improvements in our modeling of gridded seismicity hazard.

Continental United States↗
USGS1991

The Marina District, San Francisco, California: Geology, history, and earthquake effects

A northwest-trending valley in the bedrock surface is buried by firm Pleistocene bay clay, a dense Pleistocene sand layer, soft Holocene bay sediments, loose to dense Holocene beach and dune sands, and artificial fill that have an aggregate maximum thickness of about 90 m (300 ft). Artificial filling of a cove at the site of The Marina District proceeded gradually from the late 1860s to 1912, when major hydraulic filling was done for the Panama-Pacific International Exposition. The remains of thousands of piles driven for the Exposition very probably still exist and have had unknown effects on long-term ground settlement and earthquake-related ground displacements. Intensity maps of the 1906 earthquake, and seismic recordings and severe building damage in 1989, reported by others, indicate that ground motion was amplified on both natural and artificial ground. This suggests that the configuration of the bedrock surface and the location and thickness of various clay and sand deposits underlying the fill had an important effect on the shaking. However, most of the settlement and liquefaction and the damage to pipelines, building foundations, streets, side-walks, and curbs occurred in areas of artificial fill consisting mainly of loose sand.

California↗