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

At least 1,369 records · Page 76Linked to original sources

USGS2000

Rupture history of the 1997 Cariaco, Venezuela, earthquake from teleseismic P waves

A two-step finite-fault waveform inversion scheme is applied to the broadband teleseismic P waves recorded for the strike-slip, Cariaco, Venezuela, earthquake of 9 July 1997 to recover the distribution of mainshock slip. The earthquake is first analyzed using a long narrow fault with a maximum rise time of 20 sec. This line-source analysis indicates that slip propagated to the west with a constant rupture velocity and a relatively short rise time. The results are then used to constrain a second inversion of the P waveforms using a 60-km by 20-km two-dimensional fault. The rupture shows a zone of large slip (1.3-m peak) near the hypocenter and a second, broader source extending updip and to the west at depths shallower than 5 km. The second source has a peak slip of 2.1 meters and accounts for most of the moment of 1.1 × 10 26 dyne-cm (6.6 M w w) estimated from the P waves. The inferred rupture pattern is consistent with macroseismic effects observed in the epicentral area.

Geophysical Research Letters↗

USGS2009

Source analysis using regional empirical Green's functions: The 2008 Wells, Nevada, earthquake

We invert three-component, regional broadband waveforms recorded for the 21 February 2008 Wells, Nevada, earthquake using a finite-fault methodology that prescribes subfault responses using eight M W &sim;4 aftershocks as empirical Green's functions (EGFs) distributed within a 20-km by 21.6-km fault area. The inversion identifies a seismic moment of 6.2 x 10 24 dyne-cm (5.8 M W ) with slip concentrated in a compact 6.5-km by 4-km region updip from the hypocenter. The peak slip within this localized area is 88 cm and the stress drop is 72 bars, which is higher than expected for Basin and Range normal faults in the western United States. The EGF approach yields excellent fits to the complex regional waveforms, accounting for strong variations in wave propagation and site effects. This suggests that the procedure is useful for studying moderate-size earthquakes with limited teleseismic or strong-motion data and for examining uncertainties in slip models obtained using theoretical Green's functions.

Nevada↗

USGS2011

Tremor reveals stress shadowing, deep postseismic creep, and depth-dependent slip recurrence on the lower-crustal San Andreas fault near Parkfield

The 2003 magnitude 6.5 San Simeon and the 2004 magnitude 6.0 Parkfield earthquakes induced small, but significant, static stress changes in the lower crust on the central San Andreas fault, where recently detected tectonic tremor sources provide new constraints on deep fault creep processes. We find that these earthquakes affect tremor rates very differently, consistent with their differing transferred static shear stresses. The San Simeon event appears to have cast a "stress shadow" north of Parkfield, where tremor activity was stifled for 3-6 weeks. In contrast, the 2004 Parkfield earthquake dramatically increased tremor activity rates both north and south of Parkfield, allowing us to track deep postseismic slip. Following this event, rates initially increased by up to two orders of magnitude for the relatively shallow tremor sources closest to the rupture, with activity in some sources persisting above background rates for more than a year. We also observe strong depth dependence in tremor recurrence patterns, with shallower sources generally exhibiting larger, less-frequent bursts, possibly signaling a transition toward steady creep with increasing temperature and depth. Copyright 2011 by the American Geophysical Union.

California↗

USGS2024

Afterslip and creep in the rate-dependent framework: Joint inversion of borehole strain and GNSS displacements for the Mw 7.1 Ridgecrest earthquake

The elusive transition toward afterslip following an earthquake is challenging to capture with typical data resolution limits. A dense geodetic network recorded the Mw 7.1 Ridgecrest earthquake, including 16 Global Navigation Satellite System (GNSS) stations and 3 borehole strainmeters (BSM). The sub-nanostrain precision and sub-second sampling rate of BSMs bridges a gap between conventional seismologic and geodetic methods, exemplified by atypical postseismic shear strain reversals observed at nearfield (<2 km) station B921 that remain unexplained. We jointly invert GNSS displacements and BSM strains for coseismic and postseismic slip spanning hours to months over 7 independent periods. Cosiesmically, our model resolves the largest slip magnitudes of up to 6.6 m on the mainshock rupture plane, with similar patterns to other inferred slip distributions. The foreshock fault appears to slip coincidently with mainshock, revealing potential asperities activated during the preceding Mw 6.4 event. Postseismically, the best-fitting models adhere to mechanical rate-and-state expectations of logarithmically decaying slip adjacent to the coseismic rupture terminus, and where deep rheologic conditions favor creep. Most spatial variation occurs in the early postseismic timeframe (<1–2 weeks), with evidence for regional rheologic control and static stress dependence. Triggered creep on the neighboring Garlock Fault unexpectedly persists for >178 days—further highlighting the importance of fault networks in postseismic stress redistribution, critical to assessing future hazard.

Journal of Geophysics Research↗

USGS2025

Evidence for fluid pressurization of fault zones and persistent sensitivity to injection rate beneath the Raton Basin

Subsurface wastewater injection has increased the seismicity rate within the Raton Basin over more than two decades, with the basin-wide injection rate peaked between 2009-2015. To understand the evolution of injection-induced earthquakes, we systematically analyzed 2016-2024 broadband recordings with a machine-learning-based phase picker and constructed a catalog with 95,993 earthquakes (-1≤ M L ≤4.3). We then inverted for full centroid moment tensors (CMT) for 90 M L ≥ 2 events, with a special interest in constraining the non-double-couple components via probabilistic metrics. Both relocations and CMT solutions support basement-rooted normal faults, including graben and half-graben structures. Furthermore, we observe the non-double-couple components that imply elevated pore pressure in the fault zones. An earthquake cluster emerged in the north-central basin in 2023, preceded by ~1-yr of increased injection volume from wells within 15km. Despite a basin-wide decrease in the injection volume, we highlights the persistence of seismicity that remains to sensitive to injection rates within the Raton Basin.

Colorado, New Mexico↗

USGS2012

Tectonic relations of carbon dioxide discharges and earthquakes

CO 2 ‐rich springs occur worldwide along major zones of seismicity. They are mostly in young orogenic belts, but some are in areas of rifting continental platforms. Analyses of 13 C content indicate that much of the CO 2 is derived from the mantle and that other important sources are the metamorphism of marine carbonate‐bearing sedimentary rocks and the degradation of organic material. The presence of calc‐silicate minerals, such as pumpellyite in metagraywacke, is evidence of former conversion of carbonate‐bearing rocks into calc‐silicate minerals and release of CO 2 . The CO 2 pressure in fractured rocks of a fault zone reduces the effective normal stress and, if it is sufficiently great, allows the fault to slip. If the pressure were maintained at a sufficiently high level, the fault behavior might be characterized by frequent small earthquakes and aseismic creep such as occur along active segments of the San Andreas system. The presence of CO 2 ‐rich springs may indicate a potentially hazardous seismic region. Monitoring of CO 2 discharges could be useful in earthquake prediction.

Journal of Geophysical Research B: Solid Earth↗

USGS2010

The ShakeOut earthquake scenario: Verification of three simulation sets

This paper presents a verification of three simulations of the ShakeOut scenario, an M w 7.8 earthquake on a portion of the San Andreas fault in southern California, conducted by three different groups at the Southern California Earthquake Center using the SCEC Community Velocity Model for this region. We conducted two simulations using the finite difference method, and one by the finite element method, and performed qualitative and quantitative comparisons between the corresponding results. The results are in good agreement with each other; only small differences occur both in amplitude and phase between the various synthetics at ten observation points located near and away from the fault-as far as 150 km away from the fault. Using an available goodness-of-fit criterion all the comparisons scored above 8, with most above 9.2. This score would be regarded as excellent if the measurements were between recorded and synthetic seismograms. We also report results of comparisons based on time-frequency misfit criteria. Results from these two criteria can be used for calibrating the two methods for comparing seismograms. In those cases in which noticeable discrepancies occurred between the seismograms generated by the three groups, we found that they were the product of inherent characteristics of the various numerical methods used and their implementations. In particular, we found that the major source of discrepancy lies in the difference between mesh and grid representations of the same material model. Overall, however, even the largest differences in the synthetic seismograms are small. Thus, given the complexity of the simulations used in this verification, it appears that the three schemes are consistent, reliable and sufficiently accurate and robust for use in future large-scale simulations. ?? 2009 The Authors Journal compilation ?? 2009 RAS.

Geophysical Journal International↗

USGS1986

Slip deficit on the San Andreas fault at Parkfield, California, as revealed by inversion of geodetic data

A network of geodetic lines spanning the San Andreas fault near the rupture zone of the 1966 Parkfield, California, earthquake (magnitude M = 6) has been repeatedly surveyed since 1959. In the study reported here the average rates of line-length change since 1966 were inverted to determine the distribution of interseismic slip rate on the fault. These results indicate that the Parkfield rupture surface has not slipped significantly since 1966. Comparison of the geodetically determined seismic moment of the 1966 earthquake with the interseismic slip-deficit rate suggests that the strain released by the latest shock will most likely be restored between 1984 and 1989, although this may not occur until 1995. These results lend independent support to the earlier forecast of an M = 6 earthquake near Parkfield within 5 years of 1988.

California↗

USGS2014

Deformation from the 1989 Loma Prieta earthquake near the southwest margin of the Santa Clara Valley, California

Damage to pavement and near-surface utility pipes, caused by the 17 October 1989, Loma Prieta earthquake, provides evidence for ground deformation in a 663 km 2 area near the southwest margin of the Santa Clara Valley, California (USA). A total of 1427 damage sites, collected from more than 30 sources, are concentrated in four zones, three of which lie near previously mapped faults. In one of these zones, the channel lining of Los Gatos Creek, a 2-km-long concrete strip trending perpendicular to regional geologic structure, was broken by thrusts that were concentrated in two belts, each several tens of meters wide, separated by more than 300 m of relatively undeformed concrete. To gain additional measurement of any permanent ground deformation that accompanied this damage, we compiled and conducted post-earthquake surveys along two 5-km lines of horizontal control and a 15-km level line. Measurements of horizontal distortion indicate approximately 0.1 m shortening in a NE-SW direction across the valley margin, similar to the amount measured in the channel lining. Evaluation of precise leveling by the National Geodetic Survey showed a downwarp, with an amplitude of >0.1 m over a span of >12 km, that resembled regional geodetic models of coseismic deformation. Although the leveling indicates broad, regional warping, abrupt discontinuities characteristic of faulting characterize both the broad-scale distribution of damage and the local deformation of the channel lining. Reverse movement largely along preexisting faults and probably enhanced significantly by warping combined with enhanced ground shaking, produced the documented coseismic ground deformation.

California↗

USGS2024

The importance of nodal plane orientation diversity for earthquake focal mechanism stress inversions

Inversions of earthquake focal mechanisms are among the most accessible and reliable methods for determining crustal stress. However, the use of this method varies widely, and assumptions that underpin it are often violated, potentially compromising stress estimates. We investigate the consequences of violating the little-studied assumption that the focal mechanisms have diverse orientations. Our approach is to employ data-informed synthetic mechanisms, with nodal plane orientations defined by recent earthquake lineaments in the Midland Basin, western Texas, and rakes consistent with slip in the mapped stress field. Using both the traditional stress inversion method that assumes constant shear stress magnitudes on the causative faults as well as a recently published variable shear stress method, we show that low fault plane diversity can cause maximum horizontal stress ( S Hmax ) orientation and relative principal stress magnitude (faulting regime) estimates to differ markedly from the true values. This problem is compounded for catalogs with even modest amounts of noise (≤15°) or few (e.g., 20) mechanisms. Significantly, traditional approaches for quantifying uncertainty such as the bootstrap can severely underestimate the true uncertainty under these circumstances. To remedy this, we provide simple tools to quantify nodal plane orientation diversity and stress inversion reliability.

Geological Society of London Special Publications↗

USGS2015

Introduction to the special issue on the 2012 Haida Gwaii and 2013 Craig earthquakes at the Pacific–North America plate boundary (British Columbia and Alaska)

The 27 October 2012 M w 7.8 Haida Gwaii thrust earthquake and the 5 January 2013 M w 7.5 Craig strike‐slip earthquake are the focus of this special issue. They occurred along the transform boundary between the Pacific and North American plates (Fig. 1). The most identifiable feature of the plate boundary, the strike‐slip Queen Charlotte fault, might be viewed as typical of continent–ocean transform faults because it separates the continental crust of the North American plate from oceanic crust of the Pacific plate for most of its length. However, the current relative plate motion of about 5 cm/yr is highly oblique to the Queen Charlotte fault, causing a transpressive plate boundary in the region.

Alaska, British Columbia↗

USGS2020

Surface displacement distributions for the July 2019 Ridgecrest, California earthquake ruptures

Surface rupture in the 2019 Ridgecrest, California, earthquake sequence occurred along two orthogonal cross faults and includes dominantly left‐lateral and northeast‐striking rupture in the M w "> M w 6.4 foreshock and dominantly right‐lateral and northwest‐striking rupture in the M w "> M w 7.1 mainshock. We present > 650 "> > 650 field‐based, surface‐displacement observations for these ruptures and synthesize our results into cumulative along‐strike displacement distributions. Using these data, we calculate displacement gradients and compare our results with historical strike‐slip ruptures in the eastern California shear zone. For the M w "> M w 6.4 rupture, we report 96 displacements measured along 18 km of northeast‐striking rupture. Cumulative displacement curves for the rupture yield a mean left‐lateral displacement of 0.3–0.5 m and maximum of 0.7–1.6 m. Net mean vertical displacement based on the difference of down‐to‐the‐west (DTW) and down‐to‐the‐east (DTE) displacement curves is close to zero (0.02 m DTW). The M w "> M w 6.4 displacement distribution shows that the majority of displacement occurred southwest of the intersection with the M w "> M w 7.1 rupture. The M w "> M w 7.1 rupture is northwest‐striking and 50 km long based on 576 field measurements. Displacement curves indicate a mean right‐lateral displacement of 1.2–1.7 m and a maximum of 4.3–7.0 m. Net vertical displacement in the rupture averages 0.3 m DTW. The M w "> M w 7.1 displacement distributions demonstrate that maximum displacement occurred along a 12‐km‐long portion of the fault near the M w "> M w 7.1 epicenter, releasing 66% of the geologically based seismic moment along 24% of the total rupture length. Using our displacement distributions, we calculate kilometer‐scale displacement gradients for the M w "> M w 7.1 rupture. The steepest gradients ( ⁠ ∼ 1 – 3    m / km "> ∼ 1 – 3 m / km ⁠ ) flank the 12‐km‐long region of maximum displacement. In contrast, gradients for the 1992 M w "> M w 7.3 Landers and 1999 M w "> M w 7.1 Hector Mine earthquakes are < 0.6    m / km "> < 0.6 m / km . Our displacement distributions are important for understanding the influence of cross‐fault rupture on M w "> M w 6.4 and 7.1 rupture length and displacement and will facilitate comparisons with distributions generated remotely and at broader scales.

California↗

USGS2021

Post-seismic relaxation following the 2019 Ridgecrest, California, earthquake sequence

The 2019 Ridgecrest, California, earthquake sequence involved predominantly right‐lateral strike slip on a northwest–southeast‐trending subvertical fault in the 6 July M 7.1 mainshock, preceded by left‐lateral strike slip on a northeast–southwest‐trending subvertical fault in the 4 July M 6.4 foreshock. To characterize the postseismic deformation, we assemble displacements measured by Global Positioning System (GPS) and Interferometric Synthetic Aperture Radar. The geodetic measurements illuminate vigorous postseismic deformation for at least 21 months following the earthquake sequence. The postseismic transient deformation is particularly well constrained from survey‐mode GPS (sGPS) in the epicentral region carried out during the weeks after the mainshock. We interpret these observations with mechanical models including afterslip and viscoelastic relaxation of the lower crust and mantle asthenosphere. During the first 21 months, up to several centimeters of horizontal motions are measured at continuous GPS and sGPS sites, with amplitude that diminishes slowly with distance from the mainshock rupture, suggestive of deeper afterslip or viscoelastic relaxation. We find that although afterslip involving right‐lateral strike slip along the mainshock fault traces and their deeper extensions reach a few decimeters, most postseismic deformation is attributable to viscoelastic relaxation of the lower crust and mantle. Within the Basin and Range crust and mantle, we infer a transient lower crust viscosity several times that of the mantle asthenosphere. The transient mantle asthenosphere viscosity is ∼ 1.3 × 10 17    Pa   s ">

California↗

USGS2015

Ground motion observations of the 2014 South Napa earthquake

Ground motions of the South Napa earthquake (24 August 2014; M 6.0) were recorded at 19 stations within 20 km and 292 stations within 100 km of the rupture surface trace, generating peak ground motions in excess of 50% g and 50  cm/s in and near Napa Valley. This large dataset allows us to compare the ground motion from the earthquake to existing ground‐motion prediction equations (GMPEs) in considerable detail. Using the ground‐motion data compiled and reported by ShakeMap ( Wald et al. , 2000 ), we examine the peak ground acceleration (PGA) and peak ground velocity (PGV), as well as the pseudospectral acceleration (PSA) at periods of 0.3, 1.0, and 3.0 s. At the higher frequencies, especially PGA, data recorded at close distances (within &sim;20  km) are very consistent with the GMPEs, implying a stress drop for this event similar to the median for California, that is, 5 MPa ( Baltay and Hanks, 2014 ). At all frequencies, the attenuation with distance is stronger than the GMPEs would predict, which suggests the attenuation in the Napa and San Francisco Bay delta region is stronger than the average attenuation in California. The spatial plot of the ground‐motion residuals is positive to the north, in both Napa and Sonoma Valleys, consistent with increases in amplitude expected from both the directivity and basin effects. More interestingly, perhaps, there is strong ground motion to the south in the along‐strike direction, particularly for PSA at 1.0 s. These strongly positive residuals align with an older, Quaternary fault structure associated with the Franklin or Southampton fault, potentially indicating a fault‐zone‐guided wave.

California↗

USGS2020

Rapid geodetic observations of spatiotemporally varying postseismic deformation following the Ridgecrest earthquake sequence: The U.S. Geological Survey response

The U.S. Geological Survey’s geodetic response to the 4–5 July 2019 (Pacific time) Ridgecrest earthquake sequence comprised primarily the installation and/or reoccupation of Global Navigation Satellite System (GNSS) monumentation. Our response focused primarily on the United States’ Navy’s China Lake Naval Air Weapons Station base (NAWSCL). This focus was because much of the surface rupture occurred on the NAWSCL and because of NAWSCL access restrictions only permitting Federal and State of California personnel. In total, we measured or are still measuring at 24 sites, 14 of which were on the NAWSCL and, as of this writing, operational. The majority of sites were set up as continuous stations logging at either 1 sample per second or 1 sample per 15 s. Two stations were recording a 200 m cross‐rupture aperture starting ∼10 hr after the M 6.4 event, and they recorded the coseismic displacements of the M 7.1. Approximately, 1 hr after the M 7.1 event, two new stations were recording a ∼200 m cross‐rupture aperture of the surface rupture. In the days following, we established the rest of the stations ranging to a distance of ∼15 km from the M 7.1 principal rupture trace. The lack of differential displacement across the M 6.4 rupture during the M 7.1 event suggests that it did not reactivate the M 6.4 plane. The lack of differential cross‐fault displacement for both events suggests that rapid shallow afterslip did not occur at those two locations. The postseismic time series from these stations shows centimeters of horizontal displacement over periods of a few months. They record a mixture of fault‐parallel and fault‐normal displacements that, in conjunction with analysis of more spatially complete Interferometric Synthetic Aperture Radar displacement fields, suggest that both poroelastic and afterslip phenomena occur along the M 6.4 and 7.1 rupture planes. Using preliminary data from these and other regional stations, we also explore the Ridgecrest sequence’s effect on regional GNSS time series and the differentiation of long‐term postseismic motions and secular deformation rates. We find that redefining a common‐mode noise filter using different GNSS stations that are assumed to be unaffected by the earthquakes results in small but systematic differences in the regional velocity field estimate.

Seismological Research Letters↗

USGS1976

Effects of local geological conditions in the San Francisco Bay region on ground motions and the intensities of the 1906 earthquake

Measurements of ground motion generated by nuclear explosions in Nevada have been completed for 99 locations in the San Francisco Bay region, California. The recordings show marked amplitude variations in the frequency band 0.25 to 3.0 Hz that are consistently related to the local geological conditions of the recording site. The average spectral amplifications observed for vertical and horizontal ground motions are, respectively: (1, 1) for granite, (1.5, 1.6) for the Franciscan Formation, (3.0, 2.7) for the Santa Clara Formation, (3.3, 4.4) for alluvium, and (3.7, 11.3) for bay mud. Spectral amplification curves define predominant ground frequencies in the band 0.25 to 3.0 E for bay mud sites and for some alluvial sites. Amplitude spectra computed from recordings of seismic background noise at 50 sites do not generally define predominant ground frequencies. The intensities ascribed to various sites in the San Francisco Bay region for the California earthquake of April 18, 1906, are strongly dependent on distance from the zone of surface faulting and the geological character of the ground. Considering only those sites (approximately one square city block in size) for which there is good evidence for the degree of ascribed intensity, the intensities for 917 sites on Franciscan rocks generally decrease with the logarithm of distance as Intensity = 2.69 - 1.90 log (Distance in kilometers). (1) For sites on other geological units, intensity increments, derived from this empirical relation, correlate strongly with the Average Horizontal Spectral Amplifications (AHSA) according to the empirical relation Intensity Increment = 0.27 + 2.70 log (AHSA). (2) Average intensity increments predicted for the various geological units are −0.3 for granite, 0.2 for the Franciscan Formation, 0.6 for the Great Valley sequence, 0.8 for the Santa Clara Formation, 1.3 for alluvium, and 2.4 for bay mud. The maximum intensity map predicted on the basis of these data delineates areas in the San Francisco Bay region of potentially high intensity for large earthquakes on either the San Andreas fault or the Hayward fault. The map provides a crude form of seismic zonation for the region and may be useful for certain general types of land-use zonation.

California↗

USGS2002

Felt reports and intensity assignments for aftershocks and triggered events of the great 1906 California earthquake

The San Andreas fault is the longest fault in California and one of the longest strikeslip faults in the world, yet little is known about the aftershocks following the most recent great event on the San Andreas, the M 7.8 San Francisco earthquake, on 18 April 1906. This open-file report is a compilation of first-hand accounts (felt reports) describing aftershocks and triggered events of the 1906 earthquake, for the first twenty months of the aftershock sequence (through December 1907). The report includes a chronological catalog. For the larger events, Modified Mercalli intensities (MMIs) have been assigned based on the descriptions judged to be the most reliable.

California↗

USGS2004

Effect of structural heterogeneity and slip distribution on coseismic vertical displacement from rupture on the Seattle Fault

Workshops in 2001 and 2002 were convened to determine critical issues in the development of tsunami inundation maps for the Puget Sound region. The Tsunami Inundation Mapping Effort (TIME) is conducted under the multi-agency National Tsunami Hazard Mitigation Program (NTHMP). The Puget Sound Tsunami/Landslide Workshop in 2001 focused on integrated tsunami research involving a wide range of research studies and tsunami hazard mitigation issues. The 2002 Puget Sound Tsunami Sources workshop (González et al., 2003) made specific recommendations for tsunami source modeling and improving our state of knowledge for sources in the Puget Sound region. One of the recommendations stated in González et al. (2003) is "Develop methods to assess the sensitivity of coastal areas to tsunami inundation, based on multiple simulations that reflect the possible range of variations in the source parameters." Tsunami inundation models rely heavily on the imposed initial conditions which, for an earthquake source, is the coseismic vertical displacement field. For example, Koshimura et al. (2002) use the geologic uplift observations (Buknam et al., 1992) to constrain the slip distribution for the event that occurred 1100 years ago, resulting in an average slip of 3.7 m and a magnitude of 7.6. Walsh et al. (2003) develop a tsunami inundation map for Elliot Bay based on a M 7.3 earthquake and the geologic uplift observations from the 1100 y.b.p. event as in Koshimura et al. (2002), though they use a constant fault dip of 60° rather than different dips for deep and shallow segments. The objective of this report is to examine how coseismic vertical displacement from a smaller M 6.5 Seattle Fault earthquake (as in Hartzell et al., 2002) is affected by structural heterogeneity and different slip distribution patterns.

Washington↗