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,531 records · Page 85Linked to original sources

USGS1988

The 1987 Whittier Narrows earthquake in the Los Angeles metropolitan area, California

The Whittier Narrows earthquake sequence (local magnitude, M L = 5.9), which caused over 358-million dollars damage, indicates that assessments of earthquake hazards in the Los Angeles metropolitan area may be underestimated. The sequence ruptured a previously unidentified thrust fault that may be part of a large system of thrust faults that extends across the entire east-west length of the northern margin of the Los Angeles basin. Peak horizontal accelerations from the main shock, which were measured at ground level and in structures, were as high as 0.6 g (where g is the acceleration of gravity at sea level) within 50 kilometers of the epicenter. The distribution of the modified Mercalli intensity VII reflects a broad north-south elongated zone of damage that is approximately centered on the main shock epicenter.

California↗

USGS1993

The Cape Mendocino, California, earthquakes of April 1992: Subduction at the triple junction

The 25 April 1992 magnitude 7.1 Cape Mendocino thrust earthquake demonstrated that the North America—Gorda plate boundary is seismogenic and illustrated hazards that could result from much larger earthquakes forecast for the Cascadia region. The shock occurred just north of the Mendocino Triple Junction and caused strong ground motion and moderate damage in the immediate area. Rupture initiated onshore at a depth of 10.5 kilometers and propagated up-dip and seaward. Slip on steep faults in the Gorda plate generated two magnitude 6.6 aftershocks on 26 April. The main shock did not produce surface rupture on land but caused coastal uplift and a tsunami. The emerging picture of seismicity and faulting at the triple junction suggests that the region is likely to continue experiencing significant seismicity.

California↗

USGS2000

Heightened odds of large earthquakes near Istanbul: an interaction-based probability calculation

We calculate the probability of strong shaking in Istanbul, an urban center of 10 million people, from the description of earthquakes on the North Anatolian fault system in the Marmara Sea during the past 500 years and test the resulting catalog against the frequency of damage in Istanbul during the preceding millennium, departing from current practice, we include the time-dependent effect of stress transferred by the 1999 moment magnitude M = 7.4 Izmit earthquake to faults nearer to Istanbul. We find a 62 ± 15% probability (one standard deviation) of strong shaking during the next 30 years and 32 ± 12% during the next decade.

Science↗

USGS1989

Earthquakes triggered by fluid extraction

Seismicity is correlated in space and time with production from some oil and gas fields where pore pressures have declined by several tens of megapascals. Reverse faulting has occurred both above and below petroleum reservoirs, and normal faulting has occurred on the flanks of at least one reservoir. The theory of poroelasticity requires that fluid extraction locally alter the state of stress. Calculations with simple geometries predict stress perturbations that are consistent with observed earthquake locations and focal mechanisms. Measurements of surface displacement and strain, pore pressure, stress, and poroelastic rock properties in such areas could be used to test theoretical predictions and improve our understanding of earthquake mechanics.

Geology↗

USGS2017

Influence of the megathrust earthquake cycle on upper-plate deformation in the Cascadia forearc of Washington State, USA

The influence of subduction zone earthquake cycle processes on permanent forearc deformation is poorly understood. In the Cascadia subduction zone forearc of Washington State, USA, deformed and incised fluvial terraces serve as archives of longer-term (10 3 –10 4 yr) strain manifest as both fluvial incision and slip on upper-plate faults. We focus on comparing these geomorphic records in the Wynoochee River valley in the southern Olympic Mountains with short-term (10 1 yr) deformation driven by interseismic subduction zone coupling. We use optically stimulated luminescence dating and high-resolution elevation data to characterize strath terrace incision and differential uplift across the Canyon River fault, which cuts Wynoochee River terraces. This analysis demonstrates reverse slip rates of ∼0.1–0.3 mm/yr over the past ∼12–37 k.y., which agree with rates predicted by a GPS-constrained boundary element model of interseismic stress from Cascadia subduction zone coupling. Similarly, model-predicted patterns of interseismic uplift mimic the overall pattern of incision in the lower Wynoochee River valley, as revealed by strath elevations dated at 14.1 ± 1.2 ka. Agreement between modeled short-term and observed long-term records of forearc strain suggests that interseismic stress drives slip on upper-plate faults and fluvial incision in Cascadia. Consistency over multiple time scales may indicate relative stability in spatial patterns of subduction zone coupling over at least ∼10 4 yr intervals.

Washington↗

USGS2023

Rapid active thrust faulting at the northern Alaska Range front

Plate convergence rates strongly influence seismicity and mountain building inboard of convergent margins, but the distribution and kinematics of structures accommodating farfield convergence can be elusive. In interior Alaska, Yakutat microplate convergence drives late Pleistocene–recent right slip on the Denali fault, but westward-decreasing slip rates leave substantial residual Yakutat motion unaccounted for. Here, we show that Northern Foothills thrust slip beneath the northern Alaska Range absorbs a modern 4.4 mm/yr geodetic velocity gradient equivalent to ~78% of the 5.6 mm/yr residual Yakutat convergence along the central Denali fault. Infrared-stimulated luminescence ages of strath terrace deposits (67–4 ka; six sites) quantify Totatlanika River bedrock incision across the 1947 M w 7.1 thrust earthquake epicentral region. Incision rates increase abruptly from <1 mm/yr to 4.8–5.6 mm/yr above the blind thrust tip near the range front. Rapid slip at 6.7 mm/yr on a steep thrust ramp beneath the northern Alaska Range front accommodates the geodetic gradient, drives rock uplift at rates matching measured incision rates, and implies that large earthquakes like the 1947 event may recur with 500–1400 yr frequency. Results illuminate focused seismogenic strain inboard of a complex convergent margin and prompt reevaluation of Alaska’s neotectonic framework.

Alaska↗

USGS2011

No correlation between Anderson Reservoir stage level and underlying Calaveras fault seismicity despite calculated differential stress increases

Concerns have been raised that stresses from reservoir impoundment may trigger damaging earthquakes because rate changes have been associated with reservoir impoundment or stage-level changes globally. Here, the idea is tested blindly using Anderson Reservoir, which lies atop the seismically active Calaveras fault. The only knowledge held by the author going into the study was the expectation that reservoir levels change cyclically because of seasonal rainfall. Examination of seismicity rates near the reservoir reveals variability, but no correlation with stage-level changes. Three-dimensional finite-element modeling shows stress changes sufficient for earthquake triggering along the Calaveras fault zone. Since many of the reported cases of induced triggering come from low-strain settings, it is speculated that gradual stressing from stage-level changes in high-strain settings may not be significant. From this study, it can be concluded that reservoirs are not necessarily risky in active tectonic settings.

California↗

USGS2000

Probability of rupture of multiple fault segments

Fault segments identified from geologic and historic evidence have sometimes been adopted as features limiting the likely extends of earthquake ruptures. There is no doubt that individual segments can sometimes join together to produce larger earthquakes. This work is a trial of an objective method to determine the probability of multisegment ruptures. The frequency of occurrence of events on all conjectured combinations of adjacent segments in northern California is found by fitting to both geologic slip rates and to an assumed distribution of event sizes for the region as a whole. Uncertainty in the shape of the distribution near the maximum magnitude has a large effect on the solution. Frequencies of individual events cannot be determined, but it is possible to find a set of frequencies to fit a model closely. A robust conclusion for the San Francisco Bay region is that large multisegment events occur on the San Andreas and San Gregorio faults, but single-segment events predominate on the extended Hayward and Calaveras strands of segments.

Bulletin of the Seismological Society of America↗

USGS2001

Data files from “CWB Free-Field Strong-Motion Data from the 21 September Chi-Chi, Taiwan, Earthquake”

The Central Weather Bureau (CWB) of Taiwan completed a deployment of 1200 modern digital strong-motion instruments in 1996 at free-field sites and in buildings and bridges. Consequently, a very extensive set of strong-motion records were obtained for the M W = 7.6 Chi-Chi earthquake in 1999, including over 60 near-field records within 20 km of the fault ruptures. For documentation purposes, we included all relevant data files on the attached CD-ROM from our article, “CWB Free-Field Strong-Motion Data from the 21 September Chi-Chi, Taiwan, Earthquake” in this issue. We presented the data in four different ways to make the data more user friendly: (1) the original recorded data by accelerograph type; (2) the processed data in ASCII text format; (3) the processed data in SUDS format; and (4) the processed data in SAC format. (1)-(3) are exactly the same as in Lee et al. ( 2001 ), except the files are zipped (or compressed). (4) is provided for users using the SAC software on a Unix computer, and it also contains waveform plots of the data in PDF files for ease of viewing. The attached CD contains the full report of Lee et al. ( 2001 ) in 25 PDF files for viewing or printing. This full report described the data set and the data processing performed in details. It also contains plots of all processed data with an index map showing the fault, mainshock location, and the station location. We also provided a readme.txt on the CDROM to augment this Short Note with further explanations. We classify the recorded data files into four quality groups. Since the recorded strong-motion data are based on a triggering algorithm, a main concern is the amount of pre-event data and whether or not the record is long enough to cover the entire duration of the ground shaking at that station. In addition, we are concerned whether or not the record has defects (e.g., spikes, or a component was not recorded), and whether or not the record has absolute timing. In general, QA-class records are the best and can be used for any studies. QB-class records are the next best, except most of them do not have absolute timing. QC-class records cover the principal strong motions but may not have adequate pre-event data and/or post strong-motion data. QD-class records have some defects; they are included here for completeness and should not be used for most studies. Since most A800 accelerographs have a colocated A900 or A900A accelerograph, the user should use the A900 or A900A data files whenever possible. There are a readme file and six subdirectories containing the data files described in our article on the attached CDROM, under the directory of \LeeWHK: (1) readme.txt: an ASCII text file of this Short Note and additional explanations. (2) DOCfiles (subdirectory): documentary files in PDF format for viewing, and some files in ASCII text format for use in data processing. User should view or print the Reportxx.pdf files which describe the Chi-Chi strong-motion data set in details (xx is a number from 01 through 25). (3) Software (subdirectory): containing computer programs from the accelerograph manufacturers for viewing and converting the original recorded data; they are grouped by manufacturer. (4) ORGfiles (subdirectory): the original recorded data files grouped by accelerograph types; these files are for archival purposes only. (5) SUDfiles (subdirectory): the processed strong-motion data files in PC-SUDS format, and are grouped by quality class. (6) ASCfiles (subdirectory): the processed strong-motion data files are in ASCII text format and are grouped by quality class. (7) SACfiles (subdirectory): the processed strong-motion data files are in SAC format and are grouped by quality class. Within each quality class, the data are given by station. Each station has two files—xxx.pdf for viewing the three-component strong-motion data, and xxx_SAC.zip, a WINZIP file containing three SAC data files corresponding to the three components of the recorded acceleration. We use “xxx” here to denote the station name.

Bulletin of the Seismological Society of America↗

USGS2002

The Emerson Lake Body: A link between the Landers and Hector Mine earthquakes, southern California, as inferred from gravity and magnetic anomalies

Gravity and magnetic data indicate a mafic crustal heterogeneity that lies between the Hector Mine 16 October 1999 ( M w 7.1) and Landers 28 June 1992 ( M w 7.3) epicenters. The aftershocks and ruptures of these two events avoided the interior of the body. Two- and three-dimensional modeling of the potential-field anomalies shows that the source, here named the Emerson Lake body (ELB), extends to a depth of approximately 15 km. The source of the gravity and magnetic anomaly is most likely Jurassic diorite because exposures of these rocks coincide with both gravity and magnetic highs west of Emerson Lake. Seismic tomography also shows higher velocities within the region of the ELB. We propose that the ELB was an important influence on the rupture geometry of the Landers and Hector Mine ruptures and that the ELB may have played a role in transferring of stress from the Landers earthquake to the Hector Mine hypocenter. Seismicity before the Landers earthquake also tended to avoid the ELB, suggesting that the ELB affects how strain is distributed in this part of the Mojave Desert. Thus, faults within the body should have limited rupture sizes and lower seismic hazard than faults bounding or outside this mafic crustal heterogeneity.

California↗

USGS2002

The 26 January 2001 M 7.6 Bhuj, India, earthquake: Observed and predicted ground motions

Although local and regional instrumental recordings of the devastating 26, January 2001, Bhuj earthquake are sparse, the distribution of macroseismic effects can provide important constraints on the mainshock ground motions. We compiled available news accounts describing damage and other effects and interpreted them to obtain modified Mercalli intensities (MMIs) at >200 locations throughout the Indian subcontinent. These values are then used to map the intensity distribution throughout the subcontinent using a simple mathematical interpolation method. Although preliminary, the maps reveal several interesting features. Within the Kachchh region, the most heavily damaged villages are concentrated toward the western edge of the inferred fault, consistent with western directivity. Significant sediment-induced amplification is also suggested at a number of locations around the Gulf of Kachchh to the south of the epicenter. Away from the Kachchh region, intensities were clearly amplified significantly in areas that are along rivers, within deltas, or on coastal alluvium, such as mudflats and salt pans. In addition, we use fault-rupture parameters inferred from teleseismic data to predict shaking intensity at distances of 0-1000 km. We then convert the predicted hard-rock ground-motion parameters to MMI by using a relationship (derived from Internet-based intensity surveys) that assigns MMI based on the average effects in a region. The predicted MMIs are typically lower by 1-3 units than those estimated from news accounts, although they do predict near-field ground motions of approximately 80%g and potentially damaging ground motions on hard-rock sites to distances of approximately 300 km. For the most part, this discrepancy is consistent with the expected effect of sediment response, but it could also reflect other factors, such as unusually high building vulnerability in the Bhuj region and a tendency for media accounts to focus on the most dramatic damage, rather than the average effects. The discrepancy may also be partly attributable to the inadequacy of the empirical relationship between MMI and peak ground acceleration (PGA), when applied to India. The MMI-PGA relationship was developed using data from California earthquakes, which might have a systematically different stress drop and therefore, a different frequency content than intraplate events. When a relationship between response spectra and MMI is used, we obtain larger predicted MMI values, in better agreement with the observations.

Bulletin of the Seismological Society of America↗

USGS2005

Orientation of three-component geophones in the San Andreas Fault observatory at depth Pilot Hole, Parkfield, California

To identify and constrain the target zone for the planned SAFOD Main Hole through the San Andreas Fault (SAF) near Parkfield, California, a 32-level three-component (3C) geophone string was installed in the Pilot Hole (PH) to monitor and improve the locations of nearby earthquakes. The orientation of the 3C geophones is essential for this purpose, because ray directions from sources may be determined directly from the 3D particle motion for both P and S waves. Due to the complex local velocity structure, rays traced from explosions and earthquakes to the PH show strong ray bending. Observed azimuths are obtained from P-wave polarization analysis, and ray tracing provides theoretical estimates of the incoming wave field. The differences between the theoretical and the observed angles define the calibration azimuths. To investigate the process of orientation with respect to the assumed velocity model, we compare calibration azimuths derived from both a homogeneous and 3D velocity model. Uncertainties in the relative orientation between the geophone levels were also estimated for a cluster of 36 earthquakes that was not used in the orientation process. The comparison between the homogeneous and the 3D velocity model shows that there are only minor changes in these relative orientations. In contrast, the absolute orientations, with respect to global North, were significantly improved by application of the 3D model. The average data residual decreased from 13?? to 7??, supporting the importance of an accurate velocity model. We explain the remaining residuals by methodological uncertainties and noise and with errors in the velocity model.

Bulletin of the Seismological Society of America↗

USGS2008

The distribution of modified mercalli intensity in the 18 April 1906 San Francisco earthquake

We analyze Boatwright and Bundock's (2005) modified Mercalli intensity (MMI) map for the 18 April 1906 San Francisco earthquake, reviewing their interpretation of the MMI scale and testing their correlation of 1906 cemetery damage with MMI intensity. We consider in detail four areas of the intensity map where Boatwright and Bundock (2005) added significantly to the intensity descriptions compiled by Lawson (1908). We show that the distribution of off-fault damage in Sonoma County suggests that the rupture velocity approached the P-wave velocity along Tomales Bay. In contrast, the falloff of intensity with distance from the fault appears approximately constant throughout Mendocino County. The intensity in Humboldt County appears somewhat higher than the intensity in Mendocino County, suggesting that the rupture process at the northern end of the rupture was relatively energetic and that there was directivity consistent with a subsonic rupture velocity on the section of the fault south of Shelter Cove. Finally, we show that the intensity sites added in Santa Cruz County change the intensity distribution so that it decreases gradually along the southeastern section of rupture from Corralitos to San Juan Bautista and implies that the stress release on this section of rupture was relatively low.

Conference Paper↗

USGS2013

The Greenville Fault: preliminary estimates of its long-term creep rate and seismic potential

Once assumed locked, we show that the northern third of the Greenville fault (GF) creeps at 2 mm/yr, based on 47 yr of trilateration net data. This northern GF creep rate equals its 11-ka slip rate, suggesting a low strain accumulation rate. In 1980, the GF, easternmost strand of the San Andreas fault system east of San Francisco Bay, produced a M w 5.8 earthquake with a 6-km surface rupture and dextral slip growing to ≥2 cm on cracks over a few weeks. Trilateration shows a 10-cm post-1980 transient slip ending in 1984. Analysis of 2000-2012 crustal velocities on continuous global positioning system stations, allows creep rates of ~2 mm/yr on the northern GF, 0-1 mm/yr on the central GF, and ~0 mm/yr on its southern third. Modeled depth ranges of creep along the GF allow 5-25% aseismic release. Greater locking in the southern two thirds of the GF is consistent with paleoseismic evidence there for large late Holocene ruptures. Because the GF lacks large (>1 km) discontinuities likely to arrest higher (~1 m) slip ruptures, we expect full-length (54-km) ruptures to occur that include the northern creeping zone. We estimate sufficient strain accumulation on the entire GF to produce M w 6.9 earthquakes with a mean recurrence of ~575 yr. While the creeping 16-km northern part has the potential to produce a M w 6.2 event in 240 yr, it may rupture in both moderate (1980) and large events. These two-dimensional-model estimates of creep rate along the southern GF need verification with small aperture surveys.

California↗

USGS2017

Uncertainty, variability, and earthquake physics in ground‐motion prediction equations

Residuals between ground‐motion data and ground‐motion prediction equations (GMPEs) can be decomposed into terms representing earthquake source, path, and site effects. These terms can be cast in terms of repeatable (epistemic) residuals and the random (aleatory) components. Identifying the repeatable residuals leads to a GMPE with reduced uncertainty for a specific source, site, or path location, which in turn can yield a lower hazard level at small probabilities of exceedance. We illustrate a schematic framework for this residual partitioning with a dataset from the ANZA network, which straddles the central San Jacinto fault in southern California. The dataset consists of more than 3200 1.15≤ M ≤3 earthquakes and their peak ground accelerations (PGAs), recorded at close distances ( R ≤20 km). We construct a small‐magnitude GMPE for these PGA data, incorporating V S 30 site conditions and geometrical spreading. Identification and removal of the repeatable source, path, and site terms yield an overall reduction in the standard deviation from 0.97 (in ln units) to 0.44, for a nonergodic assumption, that is, for a single‐source location, single site, and single path. We give examples of relationships between independent seismological observables and the repeatable terms. We find a correlation between location‐based source terms and stress drops in the San Jacinto fault zone region; an explanation of the site term as a function of kappa, the near‐site attenuation parameter; and a suggestion that the path component can be related directly to elastic structure. These correlations allow the repeatable source location, site, and path terms to be determined a priori using independent geophysical relationships. Those terms could be incorporated into location‐specific GMPEs for more accurate and precise ground‐motion prediction.

Bulletin of the Seismological Society of America↗

USGS2021

Refinements to the Graves–Pitarka kinematic rupture generator, including a dynamically consistent slip‐rate function, applied to the 2019 Mw 7.1 Ridgecrest earthquake

The main objective of this study is to develop physics‐based constraints on the spatiotemporal variation of the slip‐rate function using a simplified dynamic rupture model. First, we performed dynamic rupture modeling of the 2019 M w 7.1 Ridgecrest, California, earthquake, to analyze the effects of depth‐dependent stress and material friction on slip rate. Then, we used our modeling results to guide refinements to the slip‐rate function that were implemented in the Graves–Pitarka kinematic rupture generation technique. The dynamic ruptures were computed on a surface‐rupturing, planar strike‐slip fault that includes a weak (negative to low‐stress‐drop) zone in the upper 4 km of the crust. Below the weak zone, we placed high‐stress‐drop patches designed to mirror the large‐slip areas seen in various rupture model inversions of the event. The locations of the high‐stress‐drop patches and the hypocenter were varied in multiple realizations to investigate how changing the dynamic conditions affected the resulting rupture kinematics, in particular, the slip rate. From these simulations, we observed a systematic change in the shape of the slip‐rate function from Kostrov type below the weak zone to a predominantly symmetric shape within the weak zone, along with a depth‐dependent reduction of peak slip rate. We generalized these shallow rupture features into a depth‐dependent parametric variation of the slip‐rate function and implemented it in the Graves–Pitarka kinematic rupture model generator. The performance of the updated kinematic approach was then verified in 0–4 Hz simulations of the M w 7.1 Ridgecrest earthquake, which showed that incorporating the depth‐dependent variation in the shape of the slip‐rate function improves the fit to the observed near‐fault ground motions in the 0.5–3 s period range.

California↗

USGS2025

Magnitude, depth and methodological variations of spectral stress drop within the SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest Earthquake Sequence

We present the first ensemble analysis of the 56 different sets of results submitted to the ongoing Community Stress Drop Validation Study using the 2019 Ridgecrest, California, earthquake sequence. Different assumptions and methods result in different estimation of the source contribution to recorded seismograms, and hence to the source parameters (principally corner frequency, f c ⁠ , spectral stress drop, Δσ, and seismic moment, M 0 ⁠ ) obtained from modeling calculated source spectra. For earthquakes smaller than magnitude (M) 2.5 there is negligible correlation between the f c values obtained by different studies, implying that no present method is reliable using available data. For larger magnitude events, correlation between f c measurements of different studies, within even a small M range is always higher than spectral ⁠Δσ , because the f c measurements simply reflect the underlying physical decrease in f c with increasing M. We model the observed trends of submitted f c with both magnitude and depth. Most methods report an increase in spectral Δσ with M, although a magnitude‐invariant spectral Δσ is within the confidence limits. The depth dependence is smaller and depends on whether a study allows attenuation to vary with source depth; a combination of depth‐dependent attenuation correction, and depth‐dependent shear‐wave velocity can compensate for reported depth trends. We model the submitted values to remove differing M and depth variation to investigate the relative interevent variability. We find consistent relative variation between individual events, and also lower relative spectral Δσ in the northwest of the aftershock sequence, and higher on the cross fault and in the region of main fault intersection. This large‐scale comparison implies that absolute spectral Δσ estimates are dependent on the methods used; studies of different regions or using different methods should not be directly compared and improved constraints on path and site corrections are needed to resolve these absolute spectral Δσ differences.

California↗

USGS2020

The U.S. Geological Survey’s Rapid Seismic Array Deployment for the 2019 Ridgecrest Earthquake Sequence

Rapid seismic deployments following large earthquakes capture ephemeral near‐field recordings of aftershocks and ambient noise that can provide valuable data for seismological studies. The U.S. Geological Survey installed 19 temporary seismic stations following the 4 July 2019 M w 6.4 and 6 July 2019 (UTC) M w 7.1 earthquakes near the city of Ridgecrest, California. The stations record the aftershock sequence beginning two days after the mainshock and are expected to remain in the field through approximately January 2020. The deployment augments the permanent seismic network in the area to improve azimuthal coverage and provide additional near‐field observations. This article summarizes the motivation and goals of the deployment; details of station installation, instrumentation, and configurations; and initial data quality and observations from the network. We expect these data to be useful for a range of studies including detailing near‐field variability in strong ground motions, determining stress drops and rupture directivity of small events, imaging the fault zone, documenting the evolution of crustal properties within and outside of the fault zone, and others.

California↗