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 217 records · Page 12Linked to original sources

USGS2003

The 2002 Denali fault earthquake, Alaska: A large magnitude, slip-partitioned event

The MW (moment magnitude) 7.9 Denali fault earthquake on 3 November 2002 was associated with 340 kilometers of surface rupture and was the largest strike-slip earthquake in North America in almost 150 years. It illuminates earthquake mechanics and hazards of large strike-slip faults. It began with thrusting on the previously unrecognized Susitna Glacier fault, continued with right-slip on the Denali fault, then took a right step and continued with right-slip on the Totschunda fault. There is good correlation between geologically observed and geophysically inferred moment release. The earthquake produced unusually strong distal effects in the rupture propagation direction, including triggered seismicity.

Science↗

USGS2010

Ground-motion modeling of Hayward fault scenario earthquakes, part I: Construction of the suite of scenarios

We construct kinematic earthquake rupture models for a suite of 39 M w 6.6-7.2 scenario earthquakes involving the Hayward, Calaveras, and Rodgers Creek faults. We use these rupture models in 3D ground-motion simulations as discussed in Part II (Aagaard et al. , 2010) to provide detailed estimates of the shaking for each scenario. We employ both geophysical constraints and empirical relations to provide realistic variation in the rupture dimensions, slip heterogeneity, hypocenters, rupture speeds, and rise times. The five rupture lengths include portions of the Hayward fault as well as combined rupture of the Hayward and Rodgers Creek faults and the Hayward and Calaveras faults. We vary rupture directivity using multiple hypocenters, typically three per rupture length, yielding north-to-south rupture, bilateral rupture, and south-to-north rupture. For each rupture length and hypocenter, we consider multiple random distributions of slip. We use two approaches to account for how aseismic creep might reduce coseismic slip. For one subset of scenarios, we follow the slip-predictable approach and reduce the nominal slip in creeping regions according to the creep rate and time since the most recent earthquake, whereas for another subset of scenarios we apply a vertical gradient to the nominal slip in creeping regions. The rupture models include local variations in rupture speed and use a ray-tracing algorithm to propagate the rupture front. Although we are not attempting to simulate the 1868 Hayward fault earthquake in detail, a few of the scenarios are designed to have source parameters that might be similar to this historical event.

Bulletin of the Seismological Society of America↗

USGS2012

Probabilistic estimates of surface coseismic slip and afterslip for Hayward fault earthquakes

We examine the partition of long‐term geologic slip on the Hayward fault into interseismic creep, coseismic slip, and afterslip. Using Monte Carlo simulations, we compute expected coseismic slip and afterslip at three alinement array sites for Hayward fault earthquakes with nominal moment magnitudes ranging from about 6.5 to 7.1. We consider how interseismic creep might affect the coseismic slip distribution as well as the variability in locations of large and small slip patches and the magnitude of an earthquake for a given rupture area. We calibrate the estimates to be consistent with the ratio of interseismic creep rate at the alinement array sites to the geologic slip rate for the Hayward fault. We find that the coseismic slip at the surface is expected to comprise only a small fraction of the long‐term geologic slip. The median values of coseismic slip are less than 0.2 m in nearly all cases as a result of the influence of interseismic creep and afterslip. However, afterslip makes a substantial contribution to the long‐term geologic slip and may be responsible for up to 0.5–1.5 m (median plus one standard deviation [S.D.]) of additional slip following an earthquake rupture. Thus, utility and transportation infrastructure could be severely impacted by afterslip in the hours and days following a large earthquake on the Hayward fault that generated little coseismic slip. Inherent spatial variability in earthquake slip combined with the uncertainty in how interseismic creep affects coseismic slip results in large uncertainties in these slip estimates.

California↗

USGS2020

The normal faulting 2020 Mw5.8 Lone Pine, Eastern California earthquake sequence

The 2020 M w "> M w 5.8 Lone Pine earthquake, the largest earthquake on the Owens Valley fault zone, eastern California, since the nineteenth century, ruptured an extensional stepover in that fault. Owens Valley separates two normal‐faulting regimes, the western margin of the Great basin and the eastern margin of the Sierra Nevada, forming a complex seismotectonic zone, and a possible nascent plate boundary. Foreshocks began on 22 June 2020; the largest M w "> M w 4.7 foreshock occurred at ∼ 6    km "> ∼6 km depth, with primarily normal faulting, followed ∼ 40    hr "> ∼ 40 hr later on 24 June 2020 by an M w "> M w 5.8 mainshock at ∼ 7    km "> ∼ 7 km depth. The sequence caused overlapping ruptures across a ∼ 0.25    km 2 "> ∼0.25 km 2 area, extended to ∼ 4    km 2 "> ∼ 4 km 2 , and culminated in an ∼ 25    km 2 "> ∼25 km 2 aftershock area. The mainshock was predominantly normal faulting, with a strike of 330° (north‐northwest), dipping 60°–65° to the east‐northeast. Comparison of background seismicity and 2020 Ridgecrest aftershock rates showed that this earthquake was not an aftershock of the Ridgecrest mainshock. The M w "> M w M w – m B "> – m B relationship and distribution of ground motions suggest typical rupture speeds. The aftershocks form a north‐northwest‐trending, north‐northeast‐dipping, 5 km long distribution, consistent with the rupture length estimated from analysis of regional waveform data. No surface rupture was reported along the 1872 scarps from the 2020 M w "> M w 5.8 mainshock, although, the dipping rupture zone of the M w "> M w 5.8 mainshock projects to the surface in the general area. The mainshock seismic energy triggered rockfalls at high elevations ( ⁠ > 3.0    km "> > 3.0 km ⁠ ) in the Sierra Nevada, at distances of 8–20 km, and liquefaction along the western edge of Owens Lake. Because there were ∼ 30 % "> ∼ 30 % fewer aftershocks than for an average southern California sequence, the aftershock forecast probabilities were lower than expected. ShakeAlert, the earthquake early warning system, provided first warning within 9.9 s, as well as subsequent updates.

California↗

USGS1996

The 25 March 1993 Scotts Mills, Oregon, earthquake and aftershock sequence: Spatial distribution, focal mechanisms, and the mount angel fault

The 25 March 1993 M L = 5.7 crustal earthquake near Scotts Mills, Oregon, was the largest earthquake to occur in the Pacific Northwest in over a decade. The mainshock was located at 45.033°N, 122.586°W and at a depth of about 15.1 km, based on arrival time data from the short-period Pacific Northwest Seismograph Network. Beginning about 12 h after the mainshock, investigators from the U.S. Geological Survey deployed 22 digital seismographs to record aftershocks. Using data from the temporary and permanent stations, we analyzed a subset of 50 aftershocks with quality locations. Hypocenters of these aftershocks lie on a northwest-trending steeply dipping plane (strike 290 ± 10°, dipping 60 ± 5° to the north-northeast), in agreement with the preferred slip plane of the mainshock focal mechanism solution (strike 294°, dipping 58° to the north-northeast). The planar structure defined by the aftershock locations may be a southeast continuation of the Mount Angel Fault, a reverse fault identified from both surface and subsurface evidence. The mapped southeast extent of the Mount Angel Fault is located less than 10 km west of the Scotts Mills epicentral region. In addition, the mainshock focal mechanism solution, with a combination of reverse motion and right-lateral strike slip, has a geometry and sense of motion consistent with the Mount Angel Fault. While aftershock focal mechanisms are varied, P axes are consistently oriented in a subhorizontal north-south direction. This earthquake sequence, together with the geological and geophysical evidence for the Mount Angel Fault, suggests a significant crustal earthquake hazard for this region of northwest Oregon.

Oregon↗

USGS2020

Seismic analysis of the 2020 Magna, Utah, earthquake sequence: Evidence for a listric Wasatch fault

The 18 March 2020 M w 5.7 Magna earthquake near Salt Lake City, Utah, offers a rare glimpse into the subsurface geometry of the Wasatch fault system—one of the world's longest active normal faults and a major source of seismic hazard in northern Utah. We analyze the Magna earthquake sequence and resolve oblique-normal slip on a shallow (30–35°) west-dipping fault at ~9- to 12-km depth. Combined with near-surface geological observations of steep dip (~70°), our results support a curved, or listric, fault shape. High-precision aftershock locations show the activation of multiple, low-angle (<30–35°) structures, indicating the existence of a complicated fault system. Our observations constrain the deep structure of the Wasatch fault system and suggest that ground shaking in the Salt Lake City region in future Wasatch fault earthquakes may be higher than previously estimated.

Utah↗

USGS1991

Mantle phase changes and deep-earthquake faulting in subducting lithosphere

Inclined zones of earthquakes are the primary expression of lithosphere subduction. A distinct deep population of subduction-zone earthquakes occurs at depths of 350 to 690 kilometers. At those depths ordinary brittle fracture and frictional sliding, the faulting processes of shallow earthquakes, are not expected. A fresh understanding of these deep earthquakes comes from developments in several areas of experimental and theoretical geophysics, including the discovery and characterization of transformational faulting, a shear instability connected with localized phase transformations under nonhydrostatic stress. These developments support the hypothesis that deep earthquakes represent transformational faulting in a wedge of olivine-rich peridotite that is likely to persist metastably in coldest plate interiors to depths as great as 690 km. Predictions based on this deep structure of mantle phase changes are consistent with the global depth distribution of deep earthquakes, the maximum depths of earthquakes in individual subductions zones, and key source characteristics of deep events.

Science↗

USGS2005

Viscoelasticity, postseismic slip, fault interactions, and the recurrence of large earthquakes

The Brownian Passage Time (BPT) model for earthquake recurrence is modified to include transient deformation due to either viscoelasticity or deep post seismic slip. Both of these processes act to increase the rate of loading on the seismogenic fault for some time after a large event. To approximate these effects, a decaying exponential term is added to the BPT model's uniform loading term. The resulting interevent time distributions remain approximately lognormal, but the balance between the level of noise (e.g., unknown fault interactions) and the coefficient of variability of the interevent time distribution changes depending on the shape of the loading function. For a given level of noise in the loading process, transient deformation has the effect of increasing the coefficient of variability of earthquake interevent times. Conversely, the level of noise needed to achieve a given level of variability is reduced when transient deformation is included. Using less noise would then increase the effect of known fault interactions modeled as stress or strain steps because they would be larger with respect to the noise. If we only seek to estimate the shape of the interevent time distribution from observed earthquake occurrences, then the use of a transient deformation model will not dramatically change the results of a probability study because a similar shaped distribution can be achieved with either uniform or transient loading functions. However, if the goal is to estimate earthquake probabilities based on our increasing understanding of the seismogenic process, including earthquake interactions, then including transient deformation is important to obtain accurate results. For example, a loading curve based on the 1906 earthquake, paleoseismic observations of prior events, and observations of recent deformation in the San Francisco Bay region produces a 40% greater variability in earthquake recurrence than a uniform loading model with the same noise level.

Bulletin of the Seismological Society of America↗

USGS1991

The 1954 and 1980 Algerian earthquakes: Implications for the characteristic-displacement model of fault behavior

Joint epicenter determination of earthquakes that occurred in northern Algeria near Ech Cheliff (named Orleansville in 1954 and El Asnam in 1980) shows that the earthquake of 9 September 1954 (M=6.5) occurred at nearly the same location as the earthquake of 10 October 1980 (M=7.3). The 1954 main shock and earliest aftershocks were concentrated close to the boundaries of segment B (nomenclature of Deschamps et al., 1982; King and Yielding, 1984) of the 1980 fault system, which was to experience approximately 8 m of slip in the 1980 earthquake. Later aftershocks of the 1954 earthquake were spread over a broad area, notably in a region north of the 1980 fault system that also experienced many aftershocks to the 1980 earthquake. The closeness of the 1954 main shock and earliest aftershocks to the 1980 segment B implies that the 1954 earthquake involved either 1) rupture of segment B proper, or 2) rupture of a distinct fault in the hanging wall of footwall block of segment B.

Bulletin of the Seismological Society of America↗

USGS2022

Simplifying complex fault data for systems-level analysis: Earthquake geology inputs for U.S. NSHM 2023

As part of the U.S. National Seismic Hazard Model (NSHM) update planned for 2023, two databases were prepared to more completely represent Quaternary-active faulting across the western United States: the NSHM23 fault sections database (FSD) and earthquake geology database (EQGeoDB). In prior iterations of NSHM, fault sections were included only if a field-measurement-derived slip rate was estimated along a given fault. By expanding this inclusion criteria, we were able to assess a larger set of faults for use in NSHM23. The USGS Quaternary Fault and Fold Database served as a guide for assessing possible additions to the NSHM23 FSD. Reevaluating available data from published sources yielded an increase of fault sections from ~650 faults in NSHM18 to ~1,000 faults proposed for use in NSHM23. EQGeoDB, a companion dataset linked to NSHM23 FSD, contains geologic slip rate estimates for fault sections included in FSD. Together, these databases serve as common input data used in deformation modeling, earthquake rupture forecasting, and additional downstream uses in NSHM development.

Scientific Data↗

USGS2015

On- and off-fault deformation associated with the September 2013 Mw7.7 Balochistan earthquake: Implications for geologic slip rate measurements

The 24 September 2013 M w 7.7 Balochistan, Pakistan earthquake ruptured a ~ 200 km-long stretch of the Hoshab fault in southern Pakistan and produced the second-largest lateral surface displacement observed for a continental strike-slip earthquake. We remotely measured surface deformation associated with this event using high-resolution (0.5 m) pre- and post-event satellite optical imagery. We document left lateral, near-field, on-fault offsets (10 m from fault) using 309 laterally offset piercing points, such as streams, terrace risers, and roads. Peak near-field displacement is 13.6 + 2.5 / − 3.4 m. We characterize off-fault deformation by measuring medium- (< 350 m from fault) and far-field (> 350 m from fault) displacement using manual (259 measurements) and automated image cross-correlation methods, respectively. Off-fault peak lateral displacement values are ~ 15 m and exceed on-fault displacement magnitudes for ~ 85% of the rupture length. Our observations suggest that for this rupture, coseismic surface displacement typically increases with distance away from the surface trace of the fault; however, nearly 100% of total surface displacement occurs within a few hundred meters of the primary fault trace. Furthermore, off-fault displacement accounts for, on average, 28% of the total displacement but exhibits a highly heterogeneous along-strike pattern. The best agreement between near-field and far-field displacements generally corresponds to the narrowest fault zone widths. Our analysis demonstrates significant and heterogeneous mismatches between on- and off-fault coseismic deformation, and we conclude that this phenomenon should be considered in hazard models based on geologically determined on-fault slip rates.

Tectonophysics↗

USGS2012

Seismic guided waves trapped in the fault zone of the Landers, California, earthquake of 1992

A mobile seismic array of seven stations was deployed at 11 sites along the fault trace of the M 7.4 Landers earthquake of June 28, 1992, with a maximum offset of 1 km from the trace. We found a distinct wave train with a relatively long period following the S waves that shows up only when both the stations and the events are close to the fault trace. This wave train is interpreted as a seismic guided wave trapped in a low‐velocity fault zone. To study the distribution of amplitude of the guided waves with distance from the fault trace and also their attenuation with travel distance along the fault zone, we eliminated source and recording site effects by the coda normalization method. The normalized amplitudes of guided waves show a spectral peak at 3–4 Hz, which decays sharply with distance from the fault trace. Spectral amplitudes at high frequencies (8–15 Hz) show an opposite trend, increasing with distance from the fault trace. The normalized amplitudes of guided waves at 3–4 Hz also show a systematic decrease with hypocentral distance along the fault zone, from which we infer an apparent Q of 50. In order to confirm the existence of the guided waves, a dense array of 31 stations was deployed at one of the 11 sites. The resultant records revealed unequivocal evidence for the existence of guided waves associated with the fault zone. By modeling the waveforms as S waves trapped in a low‐velocity waveguide sandwiched between two homogeneous half‐spaces with velocity V s = 3.0 km/s, we infer a waveguide width of about 180 m, a shear velocity of 2.0–2.2 km/s, and a Q of ∼50. Hypocenters of aftershocks with clear guided waves show a systematic distribution both laterally and with depth delineating the extent of the low‐velocity fault zone in three dimensions. We find that the zone extends to a depth of at least 10 km. This zone apparently continues to the south across the Pinto Mountain fault because guided waves are observed at stations north of the Pinto Mountain fault for earthquakes with epicenters south of it. On the other hand, the zone appears to be discontinuous at the fault bend located about 20 km north of the mainshock epicenter; guided waves were observed for stations and epicenters which are located on the same sides of the fault bend but not for those on the opposite sides.

California↗

USGS1998

Faulting apparently related to the 1994 Northridge, California, earthquake and possible co-seismic origin of surface cracks in Potrero Canyon, Los Angeles County, California

Apparent southward-dipping, reverse-fault zones are imaged to depths of about 1.5 km beneath Potrero Canyon, Los Angeles County, California. Based on their orientation and projection to the surface, we suggest that the imaged fault zones are extensions of the Oak Ridge fault. Geologic mapping by others and correlations with seismicity studies suggest that the Oak Ridge fault is the causative fault of the 17 January 1994 Northridge earthquake (Northridge fault). Our seismically imaged faults may be among several faults that collectively comprise the Northridge thrust fault system. Unusually strong shaking in Potrero Canyon during the Northridge earthquake may have resulted from focusing of seismic energy or co-seismic movement along existing, related shallow-depth faults. The strong shaking produced ground-surface cracks and sand blows distributed along the length of the canyon. Seismic reflection and refraction images show that shallow-depth faults may underlie some of the observed surface cracks. The relationship between observed surface cracks and imaged faults indicates that some of the surface cracks may have developed from nontectonic alluvial movement, but others may be fault related. Immediately beneath the surface cracks, P-wave velocities are unusually low (<400 m/sec), and there are velocity anomalies consistent with a seismic reflection image of shallow faulting to depths of at least 100 m. On the basis of velocity data, we suggest that unconsolidated soils (<800 m/sec) extend to depths of about 15 to 20 m beneath our datum (<25 m below ground surface). The underlying rocks range in velocity from about 1000 to 5000 m/sec in the upper 100 m. This study illustrates the utility of high-resolution seismic imaging in assessing local and regional seismic hazards.

California↗

USGS2004

A teleseismic study of the 2002 Denali fault, Alaska, earthquake and implications for rapid strong-motion estimation

Slip histories for the 2002 M7.9 Denali fault, Alaska, earthquake are derived rapidly from global teleseismic waveform data. In phases, three models improve matching waveform data and recovery of rupture details. In the first model (Phase I), analogous to an automated solution, a simple fault plane is fixed based on the preliminary Harvard Centroid Moment Tensor mechanism and the epicenter provided by the Preliminary Determination of Epicenters. This model is then updated (Phase II) by implementing a more realistic fault geometry inferred from Digital Elevation Model topography and further (Phase III) by using the calibrated P-wave and SH-wave arrival times derived from modeling of the nearby 2002 M6.7 Nenana Mountain earthquake. These models are used to predict the peak ground velocity and the shaking intensity field in the fault vicinity. The procedure to estimate local strong motion could be automated and used for global real-time earthquake shaking and damage assessment. ?? 2004, Earthquake Engineering Research Institute.

Earthquake Spectra↗

USGS2011

Estimating surface faulting impacts from the shakeout scenario earthquake

An earthquake scenario, based on a kinematic rupture model, has been prepared for a Mw 7.8 earthquake on the southern San Andreas Fault. The rupture distribution, in the context of other historic large earthquakes, is judged reasonable for the purposes of this scenario. This model is used as the basis for generating a surface rupture map and for assessing potential direct impacts on lifelines and other infrastructure. Modeling the surface rupture involves identifying fault traces on which to place the rupture, assigning slip values to the fault traces, and characterizing the specific displacements that would occur to each lifeline impacted by the rupture. Different approaches were required to address variable slip distribution in response to a variety of fault patterns. Our results, involving judgment and experience, represent one plausible outcome and are not predictive because of the variable nature of surface rupture. ?? 2011, Earthquake Engineering Research Institute.

Earthquake Spectra↗

USGS2004

Inverse kinematic and forward dynamic models of the 2002 Denali fault earthquake, Alaska

We perform inverse kinematic and forward dynamic models of the M 7.9 2002 Denali fault, Alaska, earthquake to shed light on the rupture process and dynamics of this event, which took place on a geometrically complex fault system in central Alaska. We use a combination of local seismic and Global Positioning System (GPS) data for our kinematic inversion and find that the slip distribution of this event is characterized by three major asperities on the Denali fault. The rupture nucleated on the Susitna Glacier thrust fault, and after a pause, propagated onto the strike-slip Denali fault. Approximately 216 km to the east, the rupture abandoned the Denali fault in favor of the more southwesterly directed Totschunda fault. Three-dimensional dynamic models of this event indicate that the abandonment of the Denali fault for the Totschunda fault can be explained by the Totschunda fault's more favorable orientation with respect to the local stress field. However, a uniform tectonic stress field cannot explain the complex slip pattern in this event. We also find that our dynamic models predict discontinuous rupture from the Denali to Totschunda fault segments. Such discontinuous rupture helps to qualitatively improve our kinematic inverse models. Two principal implications of our study are (1) a combination of inverse and forward modeling can bring insight into earthquake processes that are not possible with either technique alone, and (2) the stress field on geometrically complex fault systems is most likely not due to a uniform tectonic stress field that is resolved onto fault segments of different orientations; rather, other forms of stress heterogeneity must be invoked to explain the observed slip patterns.

Alaska↗

USGS2008

12 May 2008 M = 7.9 Wenchuan, China, earthquake calculated to increase failure stress and seismicity rate on three major fault systems

The Wenchuan earthquake on the Longmen Shan fault zone devastated cities of Sichuan, claiming at least 69,000 lives. We calculate that the earthquake also brought the Xianshuihe, Kunlun and Min Jiang faults 150-400 km from the mainshock rupture in the eastern Tibetan Plateau 0.2-0.5 bars closer to Coulomb failure. Because some portions of these stressed faults have not ruptured in more than a century, the earthquake could trigger or hasten additional M > 7 earthquakes, potentially subjecting regions from Kangding to Daofu and Maqin to Rangtag to strong shaking. We use the calculated stress changes and the observed background seismicity to forecast the rate and distribution of damaging shocks. The earthquake probability in the region is estimated to be 57-71% for M ??? 6 shocks during the next decade, and 8-12% for M ??? 7 shocks. These are up to twice the probabilities for the decade before the Wenchuan earthquake struck. Copyright 2008 by the American Geophysical Union.

Geophysical Research Letters↗

USGS1996

The 1954 Rainbow Mountain-Fairview Peak-Dixie Valley earthquakes: A triggered normal faulting sequence

In 1954, four earthquakes of M > 6.0 occurred within a 30 km radius in a period of six months. The Rainbow Mountain-Fairview Peak-Dixie Valley earthquakes are among the largest to have been recorded geodetically in the Basin and Range province. The Fairview Peak earthquake (M = 7.2, December 12, 1954) followed two events in the Rainbow Mountains ( M = 6.2, July 6, and M = 6.5, August 24, 1954) by 6 months. Four minutes later the Dixie Valley fault ruptured ( M = 6.7, December 12, 1954). The changes in static stresses caused by the events are calculated using the Coulomb-Navier failure criterion and assuming uniform slip on rectangular dislocations embedded in an elastic half-space. Coulomb stress changes are resolved on optimally oriented faults and on each of the faults that ruptured in the chain of events. These calculations show that each earthquake in the Rainbow Mountain-Fairview Peak-Dixie Valley sequence was preceded by a static stress change that encouraged failure. The magnitude of the stress increases transferred from one earthquake to another ranged from 0.01 MPa (0.1 bar) to over 0.1 MPa (1 bar). Stresses were reduced by up to 0.1 MPa over most of the Rainbow Mountain-Fairview Peak area as a result of the earthquake sequence.

Journal of Geophysical Research B: Solid Earth↗