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

USGS2018

A proposed rupture scenario for the 1925 Mw 6.5 Santa Barbara, California, earthquake

The 29 June 1925 Santa Barbara earthquake is among the largest 20th century earthquakes in southern California. The earthquake also predated the installation of strong motion and local monitoring instruments in southern California; some instrumental data are, however, available from long-period instruments at regional and teleseismic distances. The current catalog moment magnitude is M W 6.8. Initial intensity magnitudes (M I ) estimated from original Coast and Geodetic Survey intensity assignments were lower (M I 6.3). In this study we assign modified Mercalli intensity values at 239 locations, including 144 specific locations within the city of Santa Barbara for which detailed damage information is available. Comparing the reinterpreted intensities with Did You Feel it? intensities for recent events in California, we estimate M W = 6.5, with a plausible range of 6.3–6.6. We further consider reported instrumental amplitudes to estimate an instrumental moment magnitude of M W = 6.6 ± 0.5. Our preferred final estimate is M W 6.5. Based on available constraints including aftershock locations inferred from data recorded on portable instruments, we propose that the earthquake nucleated east of the city of Santa Barbara, closer to the coast than previously estimated, and ruptured unilaterally ~30 km to the west, possibly along the south-dipping Mesa-Rincon Creek, and the More Ranch fault systems. Contrary to suggestions made in earlier studies (e.g. Willis, 1925a), relatively high intensities ~50 km west of Santa Barbara can then be explained by directivity rather than involvement of the Santa Ynez fault. Finally, we discuss the possibility that the earthquake was triggered by the larger M W = 6.6 Clarkston, Montana earthquake the previous day or induced by oil production in the Summerland oil field.

California↗

USGS2026

Re-examination of the 1958 Huslia earthquake sequence and regional tectonics of the northwestern Koyukuk Basin, Alaska

In 1958 a M w 6.8 (M S = 7.3) mainshock and two M w ∼ 5.9 aftershocks occurred near the village of Huslia in the remote Koyukuk Basin of northern interior Alaska, with no surface faulting recognized afterward. This region lies in a transitional zone between compressional tectonics associated with the collision of the Pacific and Yakutat plates with North America and the extrusion of western Alaska. Since 1969 several earthquake swarms with M w 4–5 events have occurred ∼50 km northwest of Huslia in the western Purcell Mountains, including a swarm in 2019–2022 with over 9500 events. Previous waveform modeling studies conducted 30 years ago suggest the 1958 mainshock occurred on a northeast-trending normal fault. This sense of rupture is inconsistent with focal mechanisms and trends of post-1973 seismicity in the 1958 epicentral region that indicate strike-slip motion along northwest or southeast trending faults. A recently mapped east-west striking fault scarp along the southern front of the Purcell Mountains is also inconsistent with northeast-striking normal faulting. To better characterize the 1958 sequence and its relation to the regional tectonics, we measured scarp profiles across the Purcell Mountains scarp, relocated pre-1974 seismicity, analyzed first motions of pre-1974 events, and modeled potential fields (gravity, magnetics) geophysical data. Our results suggest complex regional deformation occurs along east-west, north-south and northwest-southeast trending strike-slip faults and northwest-southeast trending normal faults. The Purcell Mountains swarm may be related to hydrothermal activity as several hot springs surround the region.

Alaska↗

USGS1998

Local tsunamis and earthquake source parameters

This chapter establishes the relationship among earthquake source parameters and the generation, propagation, and run-up of local tsunamis. In general terms, displacement of the seafloor during the earthquake rupture is modeled using the elastic dislocation theory for which the displacement field is dependent on the slip distribution, fault geometry, and the elastic response and properties of the medium. Specifically, nonlinear long-wave theory governs the propagation and run-up of tsunamis. A parametric study is devised to examine the relative importance of individual earthquake source parameters on local tsunamis, because the physics that describes tsunamis from generation through run-up is complex. Analysis of the source parameters of various tsunamigenic earthquakes have indicated that the details of the earthquake source, namely, nonuniform distribution of slip along the fault plane, have a significant effect on the local tsunami run-up. Numerical methods have been developed to address the realistic bathymetric and shoreline conditions. The accuracy of determining the run-up on shore is directly dependent on the source parameters of the earthquake, which provide the initial conditions used for the hydrodynamic models.

Book chapter↗

USGS2018

Induced earthquake families reveal distinctive evolutionary patterns near disposal wells

The timing of events in seismic sequences can provide insights into the physical processes controlling fault slip. In southern Kansas, the rate of earthquakes rose rapidly starting in 2013 following expansion of energy production into the area, demanding the disposal of large volumes of wastewater into deep wells. Seismicity catalogs that are complete to low magnitudes can provide insights into the physical processes that induce seismicity near wastewater disposal. We develop a catalog of over 130,000 earthquakes recorded in southern Kansas from mid‐March 2014 through December 2017 by applying a matched filter algorithm to an original catalog of 5,831 template earthquakes. Detections have nearly identical waveforms to their associated template event and represent slip on nearly co‐located sections of a fault. We select template events with at least 100 associated detections and examine the characteristics of these prolific families of earthquakes. We find that families located close (<10 km) to areas with significant volumes of injected fluids have near‐Poissonian interevent times and the families remain active over longer durations. Families farther from high‐volume injection wells show strong clustering of interevent times and shorter sequence durations. We conclude that increasing pore fluid pressures from nearby disposal of large volumes of wastewater is the primary driver of these long duration episodes, with earthquake‐earthquake interactions driving sequences at greater distance from the wells.

Kansas↗

USGS2020

Seismic attenuation monitoring of a critically stressed San Andreas fault

We show that seismic attenuation ( ) along the San Andreas fault (SAF) at Parkfield correlates with the occurrence of moderate‐to‐large earthquakes at local and regional distances. Earthquake‐related anomalies are likely caused by changes in permeability from dilatant static stress changes, damage by strong shaking from local sources, and pore unclogging/clogging from mobilization of colloids by dynamic strains. We find that, prior to the 2004 M 6 Parkfield earthquake, prefailure conditions for some local events of moderate magnitude correspond to positive anomalies of on the Pacific side, with local and regional earthquakes producing sharp attenuation reversals. After the 2004 Parkfield earthquake, we see higher anomalies along the SAF, but low sensitivity to local and regional earthquakes, probably because the mainshock significantly altered the permeability state of the rocks adjacent to the SAF, and its sensitivity to earthquake‐induced stress perturbations.

California↗

USGS2012

Deformation from 1973 through 1991 in the epicentral area of the 1992 Landers, California, Earthquake (Ms = 7.5)

Deformation of a 50 × 60 km trilateration network that spans the epicenter of the 1992 Landers earthquake( M s = 7.5) was measured by seven surveys over the 19 years preceding the earthquake. Three moderate earthquakes (1979 Homestead Valley, M s = 5.6; 1986 North Palm Springs, M s = 6.0; and 1992 Joshua Tree, M s = 6.1) occurred within the network during those 19 years. Here we use geodetic and seismic data to construct a dislocation model for each of the three moderate earthquakes. Coseismic changes due to these three moderate earthquakes as predicted by the dislocation models are then removed from the trilateration data. The residual geodetic changes appear to be uniform in time. We take those changes to represent secular strain accumulation. This strain accumulation clearly shows right‐lateral shear across the San Andreas fault but also suggests a northwest‐southeast extension northeast of the eastern end of the “Big Bend” of the San Andreas fault. The Landers earthquake ruptured along a north‐northwest trend across this region of northwest‐southeast extension.

California↗

USGS2012

The accommodation of relative motion at depth on the San Andreas fault system in California

Plate motion below the seismogenic layer along the San Andreas fault system in California is generally assumed to occur by aseismic slip along a deeper extension of the fault. It is also possible that below the seismogenic layer, deformation is distributed laterally over a zone. Several observed features of the San Andreas fault in California have implications about the mode of accommodation of relative motion along the plate boundary beneath the seismogenic zone: the shallow depth of all earthquakes in California, the depth to which coseismic slip occurred during the 1906 San Francisco earthquake, the broad zone of strain accumulation, the broad heat flow anomaly, and the existence of widely separated parallel faults. The observations strongly imply that below the seismogenic zone, relative motion is distributed over a zone and occurs by inelastic flow rather than by aseismic slip on discrete fault planes. The existence of multiple faults further suggests that tractions at the base of the brittle layer are significant over time periods of years to hundreds of years.

Journal of Geophysical Research Solid Earth↗

USGS2012

Active faulting and deformation of the Coalinga anticline as interpreted from three-dimensional velocity structure and seismicity

This work gives a clear picture of the geometry of aftershock seismicity in a large thrust earthquake. Interpretation of hypocenters and fault plane solutions, from the 1983 Coalinga, Coast Range California, earthquake sequence, in combination with the three-dimensional velocity structure shows that the active faulting beneath the fold primarily consists of a set of southwest dipping thrusts uplifting blocks of higher-velocity material. Above the main listric blind thrust there is a conjugate fault, steeply northeast dipping, that provides the western limit of the aftershocks within the Coalinga Anticline and that corresponds in location and spatial extent with the adjacent Pleasant Valley syncline. The character of the seismicity varies with the degree of previous deformation on each section of the anticline. Where the previous uplift was largest, the shallow seismicity shows secondary faulting on either side of the fold with orientations that correspond to the preexisting geologic structure. Diffuse seismicity characterizes the area with the least previous deformation. The mainshock rupture terminated where the fold trend was no longer uniform but had competing north and west trending features. The upward extent of the mainshock rupture ended at the approximate boundary between Franciscan and Great Valley Sequence rocks. Above that depth the main thrust appears to splay into a steeper segment and a near-horizontal segment. Thus the extent of rupture area is limited by the area of uniform structural orientation and by the variation in the type of material. With the three-dimensional velocity model each individual hypocenter moved slightly (0–2 km) in accord with the details of the surrounding velocity structure, so that secondary features in the seismicity pattern are more detailed than with a local one-dimensional model and station corrections. The overall character of the fault plane solutions was not altered by the three-dimensional model, but the more accurate ray paths did result in distinct changes. In particular, the mainshock has a fault plane dipping 30° southwest instead of the 23° obtained with the one-dimensional model.

Journal of Geophysical Research Solid Earth↗

USGS2009

Earth science: lasting earthquake legacy

Earthquakes occur within continental tectonic plates as well as at plate boundaries. Do clusters of such mid-plate events constitute zones of continuing hazard, or are they aftershocks of long-past earthquakes? Early on the morning of 16 December 1811, an earthquake of about magnitude 7 shook the centre of the United States around a small town on the Mississippi called New Madrid. By 7 February 1812, it had triggered three more shocks of similar magnitude. The earthquakes broke a set of faults along the Arkansas, Missouri and Tennessee state boundaries, apparently reactivating an ancient rift in the interior of a continental tectonic plate. On 31 August 1886, a magnitude-7 shock struck Charleston, South Carolina; low-level activity continues there today. One view of seismic hazard is that large earthquakes will return to New Madrid and Charleston at intervals of about 500 years. With expected ground motions that would be stronger than average, that prospect produces estimates of earthquake hazard that rival those at the plate boundaries marked by the San Andreas fault and Cascadia subduction zone. The result is two large 'bull's-eyes' on the US National Seismic Hazard Maps — which, for example, influence regional building codes and perceptions of public safety.

Nature↗

USGS1994

Fine structure of the landers fault zone: Segmentation and the rupture process

Observations and modeling of 3- to 6-hertz seismic shear waves trapped within the fault zone of the 1992 Landers earthquake series allow the fine structure and continuity of the zone to be evaluated. The fault, to a depth of at least 12 kilometers, is marked by a zone 100 to 200 meters wide where shear velocity is reduced by 30 to 50 percent. This zone forms a seismic waveguide that extends along the southern 30 kilometers of the Landers rupture surface and ends at the fault bend about 18 kilometers north of the main shock epicenter. Another fault plane waveguide, disconnected from the first, exists along the northern rupture surface. These observations, in conjunction with surface slip, detailed seismicity patterns, and the progression of rupture along the fault, suggest that several simple rupture planes were involved in the Landers earthquake and that the inferred rupture front hesitated or slowed at the location where the rupture jumped from one to the next plane. Reduction in rupture velocity can tentatively be attributed to fault plane complexity, and variations in moment release can be attributed to variations in available energy.

California↗

USGS1994

Stress triggering of the 1994 M = 6.7 Northridge, California, earthquake by its predecessors

A model of stress transfer implies that earthquakes in 1933 and 1952 increased the Coulomb stress toward failure at the site of the 1971 San Fernando earthquake. The 1971 earthquake in turn raised stress and produced aftershocks at the site of the 1987 Whittier Narrows and 1994 Northridge ruptures. The Northridge main shock raised stress in areas where its aftershocks and surface faulting occurred. Together, the earthquakes with moment magnitude M ≥ 6 near Los Angeles since 1933 have stressed parts of the Oak Ridge, Sierra Madre, Santa Monica Mountains, Elysian Park, and Newport-Inglewood faults by more than 1 bar. Although too small to cause earthquakes, these stress changes can trigger events if the crust is already near failure or advance future earthquake occurrence if it is not.

California↗

USGS2007

Crustal structure of the Alaska Range orogen and Denali fault along the Richardson Highway

A suite of geophysical data obtained along the Richardson Highway crosses the eastern Alaska Range and Denali fault and reveals the crustal structure of the orogen. Strong seismic reflections from within the orogen north of the Denali fault dip as steeply as 25° north and extend downward to depths between 20 and 25 km. These reflections reveal what is probably a shear zone that transects most of the crust and is part of a crustal-scale duplex structure that probably formed during the Late Cretaceous. These structures, however, appear to be relict because over the past 20 years, they have produced little or no seismicity despite the nearby Mw = 7.9 Denali fault earthquake that struck in 2002. The Denali fault is nonreflective, but we interpret modeled magnetotelluric (MT), gravity, and magnetic data to propose that the fault dips steeply to vertically. Modeling of MT data shows that aftershocks of the 2002 Denali fault earthquake occurred above a rock body that has low electrical resistivity (>10 ohm-m), which might signify the presence of fluids in the middle and lower crust.

Special Paper of the Geological Society of America↗

USGS2006

Earthquake and volcano clustering via stress transfer at Yucca Mountain, Nevada

The proposed national high-level nuclear waste repository at Yucca Mountain is close to Quaternary cinder cones and faults with Quaternary slip. Volcano eruption and earthquake frequencies are low, with indications of spatial and temporal clustering, making probabilistic assessments difficult. In an effort to identify the most likely intrusion sites, we based a three-dimensional finite-element model on the expectation that faulting and basalt intrusions are sensitive to the magnitude and orientation of the least principal stress in extensional terranes. We found that in the absence of fault slip, variation in overburden pressure caused a stress state that preferentially favored intrusions at Crater Flat. However, when we allowed central Yucca Mountain faults to slip in the model, we found that magmatic clustering was not favored at Crater Flat or in the central Yucca Mountain block. Instead, we calculated that the stress field was most encouraging to intrusions near fault terminations, consistent with the location of the most recent volcanism at Yucca Mountain, the Lathrop Wells cone. We found this linked fault and magmatic system to be mutually reinforcing in the model in that Lathrop Wells feeder dike inflation favored renewed fault slip. ?? 2006 Geological Society of America.

Geology↗

USGS2016

Pseudotachylyte increases the post-slip strength of faults

Solidified frictional melts, or pseudotachylytes, are observed in exhumed faults from across the seismogenic zone. These unique fault rocks, and many experimental studies, suggest that frictional melting can be an important process during earthquakes. However, it remains unknown how melting affects the post-slip strength of the fault and why many exhumed faults do not contain pseudotachylyte. Analyses of triaxial stick-slip events on Westerly Granite (Rhode Island, USA) sawcuts at confining pressures from 50 to 400 MPa show evidence for frictional heating, including some events energetic enough to generate surface melt. Total and partial stress drops were observed with slip as high as 6.5 mm. We find that in dry samples following melt-producing stick slip, the shear failure strength increased as much as 50 MPa, while wet samples had <10 MPa strengthening. Microstructural analysis indicates that the strengthening is caused by welding of the slip surface during melt quenching, suggesting that natural pseudotachylytes may also strengthen faults after earthquakes. These results predict that natural pseudotachylyte will inhibit slip reactivation and possibly generate stress heterogeneities along faults. Wet samples do not exhibit melt welding, possibly because of thermal pressurization of water reducing frictional heating during slip.

Geology↗

USGS2003

Paleoseismic evidence of characteristic slip on the Western segment of the North Anatolian fault, Turkey

We have conducted a paleoseismic investigation of serial fault rupture at one site along the 110-km rupture of the North Anatolian fault that produced the M w 7.4 earthquake of 17 August 1999. The benefit of using a recent rupture to compare serial ruptures lies in the fact that the location, magnitude, and slip vector of the most recent event are all very well documented. We wished to determine whether or not the previous few ruptures of the fault were similar to the recent one. We chose a site at a step-over between two major strike-slip traces, where the principal fault is a normal fault. Our two excavations across the 1999 rupture reveal fluvial sands and gravels with two colluvial wedges related to previous earthquakes. Each wedge is about 0.8 m thick. Considering the processes of collapse and subsequent diffusion that are responsible for the formation of a colluvial wedge, we suggest that the two paleoscarps were similar in height to the 1999 scarp. This similarity supports the concept of characteristic slip, at least for this location along the fault. Accelerator mass spectrometry (AMS) radiocarbon dates of 16 charcoal samples are consistent with the interpretation that these two paleoscarps formed during large historical events in 1509 and 1719. If this is correct, the most recent three ruptures at the site have occurred at 210- and 280-year intervals.

North Anatolian fault↗

USGS2004

Revisiting the 23 February 1892 Laguna Salada earthquake

According to some compilations, the Laguna Salada, Baja California, earthquake of 23 February 1892 ranks among the largest earthquakes in California and Baja California in historic times. Although surface rupture was not documented at the time of the earthquake, recent geologic investigations have identified and mapped a rupture on the Laguna Salada fault that can be associated with high probability with the 1892 event (Mueller and Rockwell, 1995). The only intensity-based magnitude estimate for the earthquake, M 7.8, was made by Strand (1980) based on an interpretation of macroseismic effects and a comparison of isoseismal areas with those from instrumentally recorded earthquakes. In this study we reinterpret original accounts of the Laguna Salada earthquake. We assign modified Mercalli intensity (MMI) values in keeping with current practice, focusing on objective descriptions of damage rather than subjective human response and not assigning MMI values to effects that are now known to be poor indicators of shaking level, such as liquefaction and rockfalls. The reinterpreted isoseismal contours and the estimated magnitude are both significantly smaller than those obtained earlier. Using the method of Bakun and Wentworth (1997) we obtain a magnitude estimate of M 7.2 and an optimal epicenter less than 15 km from the center of the mapped Laguna Salada rupture. The isoseismal contours are elongated toward the northwest, which is qualitatively consistent with a directivity effect, assuming that the fault ruptured from southeast to northwest. We suggest that the elongation may also thus reflect wave propagation effects, with more efficient propagation of crustal surface (Lg) waves in the direction of the overall regional tectonic fabric.

Bulletin of the Seismological Society of America↗

USGS2026

Earthquake catalog for the Fairbanks region of central Alaska, 2014–2024, based on waveform cross-correlation

The Fairbanks region of central Alaska is part of a broad zone of intraplate crustal deformation, situated north of the Denali fault and north of the ongoing collision and flat‐slab subduction of the Yakutat oceanic plateau. Seismicity in the Fairbanks region occurs both in diffuse areas as well as in well‐defined lineaments, such as the left‐lateral Salcha fault, which hosted the 1937 M 8 7.3 earthquake. Starting with the regional seismicity catalog, we perform waveform cross‐correlation, network‐matched filtering, and relative relocation to obtain an enhanced seismicity catalog over the time period 2014–2024. Based on the relocated catalog, we interpret a set of 15 fault segments, including two conjugate faults and two new faults east of the previously documented fault system. Considering the combined seismicity in the Minto and Fairbanks regions, the median depth of seismicity decreases from east (6 km) to west (20 km). Our interpreted faults provide guidance for future tectonic modeling and assessment of seismic hazards in this region.

Alaska↗

USGS2023

Rapid characterization of the February 2023 Kahramanmaraş, Turkey, earthquake sequence

The 6 February 2023 M w "> M w 7.8 Pazarcık and subsequent M w "> M w 7.5 Elbistan earthquakes generated strong ground shaking that resulted in catastrophic human and economic loss across south‐central Türkiye and northwest Syria. The rapid characterization of the earthquakes, including their location, size, fault geometries, and slip kinematics, is critical to estimate the impact of significant seismic events. The U.S. Geological Survey National Earthquake Information Center (NEIC) provides real‐time monitoring of earthquakes globally, including rapid source characterization and impact estimates. Here, we describe the seismic characterization products generated and made available by the NEIC over the two weeks following the start of the earthquake sequence in southeast Türkiye, their evolution, and how they inform our understanding of regional seismotectonics and hazards. The kinematics of rupture for the two earthquakes was complex, involving multiple fault segments. Therefore, incorporating observations from rupture mapping was critical for characterizing these events. Dense local datasets facilitated robust source characterization and impact assessment once these observations were obtained and converted to NEIC product input formats. We discuss how we may improve the timeliness of NEIC products for rapid assessment of future seismic hazards, particularly in the case of complex ruptures.

The Seismic Record↗