USGS Science⌕ Search

Geology topics

Faults & earthquakes

Explore USGS earthquake and fault publications. Search fault names, locations and seismic hazards; follow the original report for maps, methods and publication context.

Follow this topic with RSS
336 records · Page 7Linked to original sources
USGS2025

Three-dimensional seismic velocity models for the San Francisco Bay region, California from joint body-wave and surface-wave tomography validated by waveform simulation

A high-resolution seismic velocity model for both P and S waves is essential for the San Francisco Bay (SFB) region to accurately simulate earthquake ground motions and assess seismic hazards in this high-risk region. We present new three-dimensional P-wave and S-wave velocity (Vp and Vs) models of the SFB region developed using joint tomographic inversion with a rich set of earthquake and active-source body-wave data and ambient-noise surface-wave data. We designed a dynamic balance strategy to effectively balance the contributions of body-wave and surface-wave data during each iteration of the inversion. Owing to the complementary resolution of body-wave and surface-wave data, our new models are significantly improved compared to previous tomographic models. Major geologic features in the SFB region, including low-velocity basins, high-velocity Franciscan Complex, granite, and gabbro rocks, and across-fault velocity contrasts, are imaged in our Vp and Vs models. Compared to the geology-based U.S. Geological Survey velocity model, our new model has better defined structures for some important areas that are of great concern for seismic hazard assessment, such as the Sacramento-San Joaquin Delta, and the 2014 magnitude 6.0 South Napa earthquake area. We validated our model by simulating waveforms for a set of moderate-magnitude earthquakes. The good performance of our model in waveform simulation indicates its potential for improving ground motion simulation and seismic hazard assessment in the SFB region.

California↗
USGS2026

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

Magnitude conversion relations create substantial differences in seismic hazard models

Earthquake catalogs are essential data inputs for seismic hazard modeling. Because earthquake magnitudes are reported in a variety of types (e.g., local magnitudes and moment magnitudes), magnitude conversion relationships must be used to convert the different magnitude types present in a catalog to a uniform magnitude type to avoid biases in the hazard computation. However, these conversion relationships are often uncertain and have been shown to sometimes perform poorly. Here, we investigate the sensitivity of the gridded seismicity component of the National Seismic Hazard Model (NSHM) to the catalog conversion equations in the Eastern United States. In the 2023 NSHM, magnitudes of various types were converted to moment magnitudes using equations developed by the Central and Eastern United States Seismic Source Characterization for Nuclear Facilities (CEUS‐SSCn), based on least‐squares (LS) regressions made using data from a catalog containing events up through 2008. We recompute these equations using events in the Advanced National Seismic System Comprehensive Earthquake Catalog with multiple magnitudes from 2000 to 2023. Although we prefer the use of orthogonal regressions for our datasets, LS regressions produce broadly similar results, with both approaches exhibiting large deviations from the CEUS‐SSCn conversions, especially at smaller magnitudes. We compare the spatial distribution of annual rates using three different models: (1) the 2023 NSHM conversions, (2) our updated conversions, and (3) no conversions. We find that the choice of conversions leads to substantial differences in the rate forecasts, which can greatly impact the seismic hazard model, particularly in regions with low‐seismicity rates such as the Eastern United States, where the hazard is dominated by gridded seismicity rather than a fault model.

Seismological Research Letters↗
USGS2025

Constraining large magnitude event source and path effects using ground motion simulations

The purpose of this study is to use ground motion simulations to investigate ways in which source and path effects for large magnitude events can be represented in non-ergodic GMMs. While we initially developed computation techniques using CyberShake simulations, the range of magnitudes and source-site combinations is not adequate to replicate what is observed empirically. We therefore designed a new ground motion simulation study, which includes earthquakes with a large range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of simulations (M4-M7), we then develop a non-ergodic GMM with the simulation data. We find that the within-site residuals are dominated by the radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within-site residuals among multiple source realizations. Finally, after removing the source effects from the within-site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site with the large events, while the second approach considers all small events on the fault plane. The results indicate that the path effects of large events cannot be satisfactorily approximated with that of small events using either approach.

Conference Paper↗
USGS2025

Long-term trends in microseismicity during operational shut-ins at the Coso Geothermal Field, California

Pausing injection and production can lead to induced seismicity in a variety of settings, with some of the largest events occurring during these so-called shut-ins. In geothermal fields, shut-ins are periodically conducted for maintenance on wells and surface infrastructure, thereby offering recurring means of estimating stress changes in the subsurface that lead to increased seismicity rates. Here, we examine daily production and microseismicity data from the Coso Geothermal Field (CGF) in California between 1996 and 2010 to investigate the repetitive trends of operational shut-in microseismicity. Using 24 local seismic stations, we first analyze spatial and temporal trends of over 60,000 earthquakes with magnitudes between -0.4 to 3.8. We find that the northern region exhibits no significant seismicity changes during shut-ins, whereas the rest of the field experiences induced seismicity during almost every shut-in with an increasing intensity towards the southern and eastern portions of the field, highlighting local differences in stress within the CGF. Additionally, we cluster the seismicity using waveform cross-correlation, revealing several earthquake clusters primarily occurring during shut-in periods. These observations suggest that certain fracture and fault sections respond quicker to changes in pore pressure and poroelastic stresses within the geothermal system, possibly highlighting main fluid pathways.

California↗
USGS2025

Geologic map of the Buckner 7.5-Minute quadrangle, Louisa County, Virginia

The Buckner 7.5-minute quadrangle straddles three terrane boundaries in the Piedmont Physiographic Province in central Virginia: the Chopawamsic terrane, the Elk Hill Complex, and the Goochland terrane. In much of the map area, the Elk Hill Complex separates the Chopawamsic and Goochland terranes. Rocks of the Chopawamsic terrane include Ordovician metavolcanic and metasedimentary rocks of the Chopawamsic Formation, Ordovician to Silurian granodiorite sheet intrusions, and Paleozoic mafic intrusions. Silurian to Devonian rocks of the Quantico Formation, mostly garnet-mica schist, crop out in the northwesternmost part of the map area, and are in unconformable contact with rocks of the Chopawamsic Formation on the southeastern limb of the Quantico synclinorium. The main map unit in the Elk Hill Complex is Neoproterozoic mica gneiss, which is in pre-metamorphic fault contact with rocks of the Chopawamsic Formation to the west. The main map unit of the Goochland terrane is the Maidens Gneiss. Except for Jurassic diabase dikes, all rocks on the Buckner 7.5-minute quadrangle were metamorphosed to amphibolite facies during the Alleghanian orogeny and preserve multiple compositional and phyllosilicate penetrative foliations. Evidence of amphibolite-facies metamorphism during the Taconic orogeny is preserved in rocks of the Elk Hill Complex. The entire width of the Maidens Gneiss on the Buckner 7.5-minute quadrangle is within the Spotsylvania high-strain zone and amphibolite-facies mylonitic textures are pervasive. Quartz veins and Jurassic diabase dikes crosscut all older rocks of the quadrangle. Multiple levels of terrace deposits are present along and near the major streams of the quadrangle. The lower terrace deposits are likely remnants of former positions of the Little River on the landscape, whereas higher deposits may be remnants of former deposits of the Atlantic Coastal Plain that covered this portion of the Piedmont Province. A linear cluster of aftershocks from the magnitude 5.8 earthquake that occurred near Mineral, Virginia, in 2011 defines the Fredericks Hall fault, which is at depth on the Buckner quadrangle. Most of the aftershocks occurred in the core of the Elk Creek antiform and have no relation to faults mapped at the surface. Several abandoned crushed stone and building stone quarries, as well as a mica prospect, exist in the quadrangle.

Virginia↗
USGS2025

Surficial geology and Quaternary tectonics of the Madison Valley and fault zone, Madison, Gallatin, and Beaverhead Counties, southwest Montana

The north-northwest-striking Madison fault is approximately 95 kilometers in length, lying at the confluence of the northeastern Basin and Range province and the Yellowstone tectonic parabola. The fault zone consists primarily of west-dipping normal faults that have east-dipping antithetic faults, which create the Madison Valley graben and several northeast-trending intrabasin faults. The Madison fault and associated sections discussed herein refer to the main west-dipping, range-bounding fault along the eastern side of the valley. Detailed geologic mapping (1:12,000 scale) of the entire fault zone and fault scarp profiling (total of 102 profiles) of the Madison fault reveal greater late Quaternary paleoseismic activity towards the south, including at least three paleoevents along the southern part of the fault that postdate Pinedale glaciation. Early to middle Holocene alluvial fans have vertical surface offsets that average between 2.0 and 3.0 meters and define the characteristic single-event surface offset. Pinedale lateral moraines have vertical surface offsets as great as 12.0 meters. Late Pleistocene to Holocene multiple-event fault scarps show little evidence of beveling, suggesting short seismic recurrence intervals and potential late Pleistocene and Holocene temporal clustering. Long-term average tectonic activity rates indicate slip rates ranging from 0.18–0.6 millimeters per year. Based on a comparison of fault-scarp height versus maximum slope angle of known regression lines developed from other paleoseismic investigations, the most recent event ranges from 5–1 ka. The northern section of the fault zone is defined by multiple normal faults, which detached the hanging walls of Laramide thrust faults within the Paleozoic and Mesozoic strata. This resulted in the partitioning of extension along multiple preexisting structures and less displacement along individual normal fault strands. Structural controls on lateral propagation of individual paleoevents involve the position of lateral ramps along preexisting Laramide contractional faults. This resulted in greater displacement within the larger basement-cored structures along the southern section of the fault zone, where extension is accommodated by one inferred principal basement-involved normal fault. Inferred east-northeast trending, intrabasin, normal faults within the southern half of the fault zone have no late Pleistocene displacement.

Idaho, Montana, Wyoming↗
USGS2026

Preliminary geology of the North Meadow Creek Area, Tobacco Root Mountains, southwest Montana: The North Meadow Creek Fault

This report documents a previously unmapped fault, informally referred to herein as the North Meadow Creek fault, on the east flank of the Tobacco Mountains of southwestern Montana. This fault has an apparently long and complex history, including Quaternary offset of an older alluvial terrace, offset of a debris flow of presumed Pleistocene age, offset of rhyolite lava flows of presumed Tertiary age from their presumed source several kilometers east, on the opposite side of the fault, and offset and possible rotation of blocks of Archean gneiss on northeast side of the fault. The western part of the fault is covered by glacial moraines of two apparent ages, both Pleistocene.

Montana↗
USGS2025

Map of topographic lineaments interpreted as recent surface ruptures along the Bennett Valley and Southern Maacama Fault Zones, Sonoma County, California

This study documents lidar-illuminated topographic lineaments interpreted as evidence of recent surface fault ruptures and surface ruptures related to distributed deformation along the Bennett Valley Fault Zone and the southernmost Maacama Fault Zone in the northern San Francisco Bay area (fig. 1, on map sheet). Together, these fault zones form a structural connection across a right stepover between the main Maacama Fault Zone and the Rodgers Creek Fault, overlapping principal strands of the San Andreas plate boundary system north of San Francisco Bay (figs. 1, 2, on map sheet) that accommodate about a quarter of the ~40 millimeters per year of regional dextral tectonic slip (Parsons and others, 2013). Although much of the Bennett Valley Fault Zone is seismically active (McLaughlin and others, 2012; Sowers and others, 2010), only the north end of the fault zone (the “Spring Valley strand,” indicated on southern map panel of map sheet and labeled “D” in figure 2, on map sheet) was previously known to have youthful geomorphic expression and to displace Holocene deposits (McLaughlin and others, 2008; Sowers and others, 2010; Sowers and others, 2016). Holocene activity along the length of the Bennett Valley Fault Zone and southernmost Maacama Fault Zone, with possible implications for rupture propagation, continuity, and slip transfer, had not been identified. However, fault splays that project eastward from the Rodgers Creek Fault have been mapped and hypothesized to accommodate slip transfer to the Bennett Valley Fault Zone (Hecker and Randolph Loar, 2018). The generally subtle and distributed nature of surface-rupture evidence along the Bennett Valley and southernmost Maacama Fault Zones, and extensive vegetation cover, had left recent faulting previously unmapped along most of the zone. The map presented here represents a new compilation of inferred surface-rupture features detected using high-resolution topographic lidar data from an airborne lidar survey of Sonoma County, California (OpenTopography, 2014). These data, which enable subtle topographic features to be discerned, indicate that recent (likely Holocene) surface ruptures extend throughout the Bennett Valley and southernmost Maacama Fault Zones. The rupture-lineament map was created using a GIS (geographic information system) framework and is included herein as an image map at a scale of 1:36,000 and as digital datasets (included as supplemental information to this report). The mapping is intended to lay the groundwork for future studies designed to better characterize how plate-boundary slip is accommodated on this important and complex system of faults.

California↗
USGS2025

A newly identified creeping strand of the Concord fault, San Francisco Bay Area

The Concord fault constitutes a major branch of the Pacific–North America transform plate boundary in Northern California, bridging the strike‐slip Bartlett Springs ‐ Green Valley Fault system to the north with the Greenville and Calaveras Faults to the south. Like many faults in the San Francisco Bay Area its long‐term slip is partially accommodated by aseismic slip (creep). Although creep has been recognized and monitored on the northern half of the fault for decades, the precise location of the southern half of the fault and its slip rate—whether accommodated seismically or aseismically—has remained enigmatic. How slip transfers between the Concord and Greenville or Calaveras faults to the south remains an outstanding question. New field observations presented here indicate that the active trace of the fault south of downtown Concord is not where previously interpreted and is indeed actively creeping. We report observations of shallow creep continuing >7 km farther south along the Concord fault than previously reported, along a fault strand not previously recognized for most of its length. This is evident as right‐laterally deflected concrete curbs and sidewalk slabs on both sides of every street that crosses the fault at a high angle in southeast Concord and northeast Walnut Creek. We document the magnitude and location of these deflections to estimate accumulated right‐lateral aseismic slip expressed in engineered structures. Offsets of these piercing lines range from 8 to 18 cm, over widths varying from narrow breaks along centimeter‐scale concrete joints to 10‐m‐wide zones of deflection. Significantly, this active trace is ∼400 m west of where the Quaternary active trace has previously been inferred, placing it within—rather than bounding—the built area of suburban Concord. Slip along the fault has already caused infrastructure damage. These results revise our understanding of the southern Concord fault and help constrain its seismic potential.

California↗
USGS2025

Late Pleistocene kinematics of the Great Southern Puerto Rico Fault Zone, Puerto Rico

Several onshore faults in southern Puerto Rico have recently been recognized as Quaternary active. However, the kinematics of these faults, particularly any lateral component, remain largely unconstrained. It is difficult to characterize low strain‐rate faults, partially due to extensive erosional and anthropogenic landscape modification, steep relief, and frequent landsliding, limiting the preservation of geomorphic features that could serve as recorders of fault motion. Here, we constrain the kinematics along sections of the Great Southern Puerto Rico Fault Zone (GSPRFZ) on the southern coastal plain of Puerto Rico. We integrate ∼1‐m‐resolution light detection and ranging (lidar)‐derived topography, historical air photos, and field mapping to identify a series of ∼50–1200‐m‐long fault scarps and lineaments that trend northwest–southeast and extend for ≥25 km across the southern coastal plain. Fault scarps are primarily south facing, cut across topography, and displace Quaternary deposits and landforms. We document multiple offset geomorphic markers, including channel thalwegs and interfluves formed in deposits previously mapped as Quaternary piedmont alluvial plain. We observe both vertical (south‐side‐down) and right‐lateral meter‐scale displacements, which indicate that the GSPRFZ accommodates right‐lateral oblique motion in the late Pleistocene, consistent with northeast motion of the Puerto Rico and the Virgin Islands microplate away from the Hispaniola block.

Puerto Rico↗
USGS2025

Anomalous shear stress variation in wet granular medium: Implications for landslide lateral faults

Landslide assessments typically focus on the mechanical properties of the basal shear zone, but lateral faults are frequently overlooked, possibly due to their lower normal stresses and variably saturated conditions. Using double-cylinder shear experiments on wet granular systems as analogs for landslide lateral faults, we observe anomalous shear stress variations with fluid volume fractions, defying an expected unimodal relationship associated with capillary cohesion. At low fluid volume fractions, shear strength weakens as the wet grain assembly experiences reduced lateral pressure and increased boundary slip. This boundary slip subsequently vanishes, with an abrupt strengthening due to the dilation of the grain assembly against fluid surface tension as saturation approaches. Strike-slip motion and confinement in this system explain the strength anomaly, highlighting a critical role of lateral faults in landslide stability, particularly in cases where dynamics cannot be adequately explained by monitored pore-water pressure or basal friction.

Geophysical Research Letters↗
USGS2026

Structural evolution and slip rate variations through time of the Puente-Hills blind-thrust fault beneath Los Angeles: Implications for seismic hazard and folding kinematics

Using seismic reflection profiles, historical well logging data, and luminescence and radiocarbon ages, we determine a Pleistocene-Holocene slip history for the central, Santa Fe Springs segment of the Puente Hills blind-thrust fault (PHT), a major seismogenic fault situated beneath the urbanized Los Angeles metropolitan region. We analyze the geometry of correlative stratigraphic units in the forelimb and backlimb of the overlying growth-fold, the Santa Fe Springs anticline, and determine the uplift of seven age-correlative markers. Uplift measurements are converted to thrust displacements on the underlying PHT using a structural method laid out by Don et al. (2022, https://doi.org/10.1785/0120220048 ) that accounts for the geometry of the fault. These data indicate that deep thrust displacement on the PHT is partially consumed updip in the creation of a hanging-wall fault-bend fold, with forelimb growth strata recording <80% of the slip documented within the backlimb. Chronological data from growth strata yield age constraints for folding and faulting on the underlying PHT, providing a detailed incremental slip history derived from both the forelimb and backlimb folding for seven discrete growth horizons, spanning the past 1.4 million years. The resulting six incremental slip rates demonstrate that fault slip has varied through time from the middle Pleistocene to the Holocene. Moreover, these results reveal synchronous acceleration of both the central, Santa Fe Springs and western, Los Angeles segments of the fault system since late Pleistocene time (after 200 ka) and slip rates of greater than 2 mm/yr on the downdip, backlimb, fault-ramp below the anticline.

California↗
USGS2025

Comparative properties of saponitic fault gouge and serpentinite muds cored from mud volcanoes of the Mariana subduction zone

We obtained 12 core samples for physical and chemical characterization from three serpentinite mud volcanoes (Yinazao, Asùt Tesoru, and Fantangisña) located on the forearc of the Mariana subduction system, that were drilled during International Ocean Discovery Program Expedition 366. Two samples from the Fantangisña mud volcano are interpreted to be clay-rich fault gouges derived from the subduction channel. Their bulk compositions are intermediate between the serpentinites and oceanic basalts. The oceanic crustal materials in the gouges have been thoroughly metasomatized and the serpentinites extensively altered to the trioctahedral, Mg-rich smectite clays saponite and corrensite. The only relict phases in clasts of crustal rock are accessory Ti- and P-bearing minerals. The two fault gouge samples have lower frictional strengths ( μ < 0.2) than the serpentinites ( μ = 0.2–0.4), and their measured permeabilities are also somewhat lower. Their physical and compositional properties correspond to saponitic gouges from other faults that juxtapose serpentinite against crustal rocks, in particular gouges from the two creeping traces of the San Andreas Fault recovered in the core from the San Andreas Fault Observatory at Depth. The décollement beneath Fantangisña mud volcano is thus expected to be very weak and likely characterized by stable slip.

Geophysics, Geochemistry, Geosystems↗
USGS2025

Memory and jamming in fault zone sediments

Many subsurface processes involve transitions in granular material states, from arrested to creeping to flowing. Experiments and frameworks for idealized systems reveal that granular fabrics develop during shearing, co-evolve with applied stress, and govern such transitions. We use microtomography to test whether fabrics at two San Andreas fault sites reflect slip history and whether idealized frameworks extend to nature. Near-surface sediments within the fault zone transition between deformation patterns over the seismic cycle, including bulk/localized grain re-arrangements, individual grain fracturing, and localized zones of fracturing. Aseismic and co-seismic shearing produce distinct preferred grain orientations. Co-seismic fabrics can be preserved after centuries of aseismic strain, aseismic fabrics may be overprinted, and grain size and coordination number influence the fabrics. Idealized frameworks, namely anisotropic critical state theory, frictional jamming, and material memory, can explain our observations, and fault zone sediments likely undergo cycles of memory creation and erasure that influence rigidity spatiotemporally.

California↗
USGS2026

Discrete element investigation of the influence of shallow soil density on the manifestations of strike-slip surface fault rupture

This study investigates the influence of soil relative density on strike-slip surface fault rupture manifestations using three-dimensional numerical simulations performed with the discrete element method (DEM). The simulations capture the formation of distinctive fault strands within complex flower structures using tens of millions of grains. The tendency for dense soils to localize shear manifests as multiple localized shear bands within wide zones of deformation in strike-slip fault rupture, whereas diffuse shear deformation develops within narrow zones in loose soils. The spatial extents of soil deformation are consistent between simulations having similar relative densities but with different quantities of grains. However, individual shears are more distinguishable in assemblages of finer grains than in assemblages of coarser grains. The simulations show the progressive development of new shears within the bounds of previously developed shears. Shear activity transitions inward as fault activity diminishes along the outermost shears and continues along newly developed shears until a vertically dipping throughgoing shear structure develops that accommodates most of the fault displacement thereon. The throughgoing fault develops at smaller fault displacements in looser soils because the first shear rupture propagates closer to the vertical direction and does not undergo as much inward translation of shear activity, as is observed in denser soils. In all simulations, ground surface uplift develops between nonintersecting active shears, and ground surface subsidence tends to develop where new shears intersect previous shears. The surface traces in these simulations are shown to be consistent with analog models and case histories of surface fault rupture occurring in different shallow subsurface materials. Although computationally costly, these modeling results are valuable for providing a strong numerical supplement to traditional analog models used to represent the mechanics of strike-slip zones in soil, and they provide quantifiable stresses and large-strain deformations throughout the model domain.

Journal of Geotechnical and Geoenvironmental Engin↗