USGS Science⌕ Search

Geology topics

Faults & earthquakes

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

At least 271 records · Page 15Linked to original sources

Buried shallow fault slip from the South Napa earthquake revealed by near-field geodesy

Earthquake-related fault slip in the upper hundreds of meters of Earth’s surface has remained largely unstudied because of challenges measuring deformation in the near field of a fault rupture. We analyze centimeter-scale accuracy mobile laser scanning (MLS) data of deformed vine rows within ±300 m of the principal surface expression of the M (magnitude) 6.0 2014 South Napa earthquake. Rather than assuming surface displacement equivalence to fault slip, we invert the near-field data with a model that allows for, but does not require, the fault to be buried below the surface. The inversion maps the position on a preexisting fault plane of a slip front that terminates ~3 to 25 m below the surface coseismically and within a few hours postseismically. The lack of surface-breaching fault slip is verified by two trenches. We estimate near-surface slip ranging from ~0.5 to 1.25 m. Surface displacement can underestimate fault slip by as much as 30%. This implies that similar biases could be present in short-term geologic slip rates used in seismic hazard analyses. Along strike and downdip, we find deficits in slip: The along-strike deficit is erased after ~1 month by afterslip. We find no evidence of off-fault deformation and conclude that the downdip shallow slip deficit for this event is likely an artifact. As near-field geodetic data rapidly proliferate and will become commonplace, we suggest that analyses of near-surface fault rupture should also use more sophisticated mechanical models and subsurface geomechanical tests.

California↗

Evidence for large Holocene earthquakes along the Denali fault in southwest Yukon, Canada

The Yukon–Alaska Highway corridor in southern Yukon is subject to geohazards ranging from landslides to floods and earthquakes on faults in the St. Elias Mountains and Shakwak Valley. Here we discuss the late Holocene seismic history of the Denali fault, located at the eastern front of the St. Elias Mountains and one of only a few known seismically active terrestrial faults in Canada. Holocene faulting is indicated by scarps and mounds on late Pleistocene drift and by tectonically deformed Pleistocene and Holocene sediments. Previous work on trenches excavated against the fault scarp near the Duke River reveals paleoseismic sediment disturbance dated to ∼300–1,200, 1,200–1,900, and 3,000 years ago. Re-excavation of the trenches indicates a fourth event dated to 6,000 years ago. The trenches are interpreted to show a negative flower structure produced by extension of sediments by dextral strike-slip fault movement. Nearby Crescent Lake is ponded against the fault scarp. Sediment cores reveal four abrupt sediment and diatom changes reflecting seismic shaking at ∼1,200–1,900, 1,900–5,900, 5,900–6,200, and 6,500–6,800 years ago. At the Duke River, the fault offsets sediments, including two White River tephra layers (∼1,900 and 1,200 years old). Late Pleistocene outwash gravel and overlying Holocene aeolian sediments show in cross section a positive flower structure indicative of post-glacial contraction of the sediments by dextral strike-slip movement. Based on the number of events reflecting ∼M6, we estimate the average recurrence of large earthquakes on the Yukon part of the Denali fault to be about 1,300 years in the past 6,500–6,800 years.

Yukon↗

Field and Laboratory Data From an Earthquake History Study of Scarps in the Hanging Wall of the Tacoma Fault, Mason and Pierce Counties, Washington

As part of the effort to assess seismic hazard in the Puget Sound region, we map fault scarps on Airborne Laser Swath Mapping (ALSM, an application of LiDAR) imagery (with 2.5-m elevation contours on 1:4,000-scale maps) and show field and laboratory data from backhoe trenches across the scarps that are being used to develop a latest Pleistocene and Holocene history of large earthquakes on the Tacoma fault. We supplement previous Tacoma fault paleoseismic studies with data from five trenches on the hanging wall of the fault. In a new trench across the Catfish Lake scarp, broad folding of more tightly folded glacial sediment does not predate 4.3 ka because detrital charcoal of this age was found in stream-channel sand in the trench beneath the crest of the scarp. A post-4.3-ka age for scarp folding is consistent with previously identified uplift across the fault during AD 770-1160. In the trench across the younger of the two Stansberry Lake scarps, six maximum 14C ages on detrital charcoal in pre-faulting B and C soil horizons and three minimum ages on a tree root in post-faulting colluvium, limit a single oblique-slip (right-lateral) surface faulting event to AD 410-990. Stratigraphy and sedimentary structures in the trench across the older scarp at the same site show eroded glacial sediments, probably cut by a meltwater channel, with no evidence of post-glacial deformation. At the northeast end of the Sunset Beach scarps, charcoal ages in two trenches across graben-forming scarps give a close maximum age of 1.3 ka for graben formation. The ages that best limit the time of faulting and folding in each of the trenches are consistent with the time of the large regional earthquake in southern Puget Sound about AD 900-930.

Scientific Investigations Map↗

Near-field ground motion of the 2002 Denali fault, Alaska, earthquake recorded at pump station 10

A free-field recording of the Denali fault earthquake was obtained by the Alyeska Pipeline Service Company 3 km from the surface rupture of the Denali fault. The instrument, part of the monitoring and control system for the trans-Alaska pipeline, was located at Pump Station 10, approximately 85 km east of the epicenter. After correction for the measured instrument response, we recover a seismogram that includes a permanent displacement of 3.0 m. The recorded ground motion has relatively low peak acceleration (0.36 g) and very high peak velocity (180 cm/s). Nonlinear soil response may have reduced the peak acceleration to this 0.36 g value. Accelerations in excess of 0.1 g lasted for 10 s, with the most intense motion occurring during a 1.5-s interval when the rupture passed the site. The low acceleration and high velocity observed near the fault in this earthquake agree with observations from other recent large-magnitude earthquakes. ?? 2004, Earthquake Engineering Research Institute.

Earthquake Spectra↗

Tidal triggering of low frequency earthquakes near Parkfield, California: Implications for fault mechanics within the brittle-ductile transition

Studies of nonvolcanic tremor (NVT) have established the significant impact of small stress perturbations on NVT generation. Here we analyze the influence of the solid earth and ocean tides on a catalog of ∼550,000 low frequency earthquakes (LFEs) distributed along a 150 km section of the San Andreas Fault centered at Parkfield. LFE families are identified in the NVT data on the basis of waveform similarity and are thought to represent small, effectively co-located earthquakes occurring on brittle asperities on an otherwise aseismic fault at depths of 16 to 30 km. We calculate the sensitivity of each of these 88 LFE families to the tidally induced right-lateral shear stress (RLSS), fault-normal stress (FNS), and their time derivatives and use the hypocentral locations of each family to map the spatial variability of this sensitivity. LFE occurrence is most strongly modulated by fluctuations in shear stress, with the majority of families demonstrating a correlation with RLSS at the 99% confidence level or above. Producing the observed LFE rate modulation in response to shear stress perturbations requires low effective stress in the LFE source region. There are substantial lateral and vertical variations in tidal shear stress sensitivity, which we interpret to reflect spatial variation in source region properties, such as friction and pore fluid pressure. Additionally, we find that highly episodic, shallow LFE families are generally less correlated with tidal stresses than their deeper, continuously active counterparts. The majority of families have weaker or insignificant correlation with positive (tensile) FNS. Two groups of families demonstrate a stronger correlation with fault-normal tension to the north and with compression to the south of Parkfield. The families that correlate with fault-normal clamping coincide with a releasing right bend in the surface fault trace and the LFE locations, suggesting that the San Andreas remains localized and contiguous down to near the base of the crust. The deep families that have high sensitivity to both shear and tensile normal stress perturbations may be indicative of an increase in effective fault contact area with depth. Synthesizing our observations with those of other LFE-hosting localities will help to develop a comprehensive understanding of transient fault slip below the “seismogenic zone” by providing constraints on parameters in physical models of slow slip and LFEs.

California↗

Coseismic deformation of the Wrights tunnel during the 1906 San Francisco earthquake: A key to understanding 1906 fault slip and 1989 surface ruptures in the southern Santa Cruz Mountains, California

The Wrights tunnel is an abandoned railroad tunnel that crosses the San Andreas fault in the southern Santa Cruz Mountains in the vicinity of the 1989 Loma Prieta earthquake. The tunnel was damaged and deformed during the 1906 San Francisco earthquake and a plot showing postearthquake measurements made in the tunnel is given by Lawson [1908]. The amount of offset shown on this plot (1.5 m) has been used in several studies as being representative of the amount of fault offset along this segment of the San Andreas fault in 1906. Our historical research shows that different observers reported different amounts of fault offset in the tunnel and that the 1.5 m given on the plot is not a surveyed measurement. In addition, the plot of the tunnel has been interpreted in several previous studies as evidence of a broad (1.5 km) zone of faulting beneath Summit Ridge. Our analysis shows that this plot need not indicate a broad zone of deformation. Our historical research and modeling of the tunnel measurements indicate that faulting was confined to a zone less than 400 m wide and that 60–85% of the coseismic slip occurred across a single fault plane. There is no evidence for offset across a second shear zone beneath Summit Ridge in 1906. This implies that surface fractures reported on Summit Ridge in 1906 were not associated with significant deformation of the tunnel, implying that they were shallow, surficial features. By analogy, the very similar fractures that occurred on Summit Ridge in 1989 were also probably the result of shallow gravitational, rather than deep-seated tectonic, processes. Our modeling also indicates that total coseismic, near-surface slip across the San Andreas fault zone in the Wrights tunnel in 1906 was at least 1.7–1.8 m.

Journal of Geophysical Research B: Solid Earth↗

Surface rupture map of the 2002 M7.9 Denali fault earthquake, Alaska: Digital data

The November 3, 2002, Mw7.9 Denali Fault earthquake produced about 340 km of surface rupture along the Susitna Glacier Thrust Fault and the right-lateral, strike-slip Denali and Totschunda Faults. Digital photogrammetric methods were primarily used to create a 1:500-scale, three-dimensional surface rupture map, and 1:6,000-scale aerial photographs were used for three-dimensional digitization in ESRI's ArcMap GIS software, using Leica's StereoAnalyst plug in. Points were digitized 4.3 m apart, on average, for the entire surface rupture. Earthquake-induced landslides, sackungen, and unruptured Holocene fault scarps on the eastern Denali Fault were also digitized where they lay within the limits of air photo coverage. This digital three-dimensional fault-trace map is superior to traditional maps in terms of relative and absolute accuracy, completeness, and detail and is used as a basis for three-dimensional visualization. Field work complements the air photo observations in locations of dense vegetation, on bedrock, or in areas where the surface trace is weakly developed. Seventeen km of the fault trace, which broke through glacier ice, were not digitized in detail due to time constraints, and air photos missed another 10 km of fault rupture through the upper Black Rapids Glacier, so that was not mapped in detail either.

Alaska↗

Far-field pressurization likely caused one of the largest injection induced earthquakes by reactivating a large pre-existing basement fault structure

The M w 5.1 Fairview, Oklahoma, earthquake on 13 February 2016 and its associated seismicity produced the largest moment release in the central and eastern United States since the 2011 M w 5.7 Prague, Oklahoma, earthquake sequence and is one of the largest earthquakes potentially linked to wastewater injection. This energetic sequence has produced five earthquakes with M w 4.4 or larger. Almost all of these earthquakes occur in Precambrian basement on a partially unmapped 14 km long fault. Regional injection into the Arbuckle Group increased approximately sevenfold in the 36 months prior to the start of the sequence (January 2015). We suggest far-field pressurization from clustered, high-rate wells greater than 12 km from this sequence induced these earthquakes. As compared to the Fairview sequence, seismicity is diffuse near high-rate wells, where pressure changes are expected to be largest. This points to the critical role that preexisting faults play in the occurrence of large induced earthquakes.

Geophysical Research Letters↗

Fault structure and kinematics of the Long Valley Caldera region, California, revealed by high-accuracy earthquake hypocenters and focal mechanism stress inversions

We have determined high-resolution hypocenters for 45,000+ earthquakes that occurred between 1980 and 2000 in the Long Valley caldera area using a double-difference earthquake location algorithm and routinely determined arrival times. The locations reveal numerous discrete fault planes in the southern caldera and adjacent Sierra Nevada block (SNB). Intracaldera faults include a series of east/west-striking right-lateral strike-slip faults beneath the caldera's south moat and a series of more northerly striking strike-slip/normal faults beneath the caldera's resurgent dome. Seismicity in the SNB south of the caldera is confined to a crustal block bounded on the west by an east-dipping oblique normal fault and on the east by the Hilton Creek fault. Two NE-striking left-lateral strike-slip faults are responsible for most seismicity within this block. To understand better the stresses driving seismicity, we performed stress inversions using focal mechanisms with 50 or more first motions. This analysis reveals that the least principal stress direction systematically rotates across the studied region, from NE to SW in the caldera's south moat to WNW-ESE in Round Valley, 25 km to the SE. Because WNW-ESE extension is characteristic of the western boundary of the Basin and Range province, caldera area stresses appear to be locally perturbed. This stress perturbation does not seem to result from magma chamber inflation but may be related to the significant (???20 km) left step in the locus of extension along the Sierra Nevada/Basin and Range province boundary. This implies that regional-scale tectonic processes are driving seismic deformation in the Long Valley caldera.

California↗

Fault interactions and large complex earthquakes in the Los Angeles area

Faults in complex tectonic environments interact in various ways, including triggered rupture of one fault by another, that may increase seismic hazard in the surrounding region. We model static and dynamic fault interactions between the strike-slip and thrust fault systems in southern California. We find that rupture of the Sierra Madre-Cucamonga thrust fault system is unlikely to trigger rupture of the San Andreas or San Jacinto strike-slip faults. However, a large northern San Jacinto fault earthquake could trigger a cascading rupture of the Sierra Madre-Cucamonga system, potentially causing a moment magnitude 7.5 to 7.8 earthquake on the edge of the Los Angeles metropolitan region.

Science↗

Earthquake hypocenters and focal mechanisms in central Oklahoma reveal a complex system of reactivated subsurface strike-slip faulting

The sharp increase in seismicity over a broad region of central Oklahoma has raised concern regarding the source of the activity and its potential hazard to local communities and energy industry infrastructure. Since early 2010, numerous organizations have deployed temporary portable seismic stations in central Oklahoma in order to record the evolving seismicity. In this study, we apply a multiple-event relocation method to produce a catalog of 3,639 central Oklahoma earthquakes from late 2009 through 2014. RMT source parameters were determined for 195 of the largest and best-recorded earthquakes. Combining RMT results with relocated seismicity enabled us to determine the length, depth and style-of-faulting occurring on reactivated subsurface fault systems. Results show that the majority of earthquakes occur on near vertical, optimally oriented (NE-SW and NW-SE), strike-slip faults in the shallow crystalline basement. These are necessary first order observations required to assess the potential hazards of individual faults in Oklahoma.

Oklahoma↗

Landslides triggered by the 2002 M 7.9 Denali Fault earthquake, Alaska, USA

The 2002 M 7.9 Denali earthquake in Alaska, USA, was the largest inland earthquake in North America in nearly 150 years. The earthquake involved oblique thrusting but mostly strike-slip motion, and faults ruptured the ground surface over 330 km. Fault rupture occurred in a rugged, mountainous, subarctic environment with extensive permafrost and variable glaciation, geology, and groundwater presence, and many triggered landslides mobilized into avalanches that traversed varied physiographic settings, some moving farther than 11 km. However, only several thousand landslides were triggered, and these occurred in a narrow zone along the fault rupture. These characteristics of the event provide ample opportunities to improve understanding of controls on coseismic landsliding and avalanching mechanisms. The paucity and limited extent of landslides likely resulted from high directivity of seismic energy and relatively low levels of high-amplitude, high-frequency ground motion; glacial damping of seismic energy was likely not a factor. Landslides preferentially occurred on hillslopes steeper and higher than average. Meteorological conditions were near historical averages at the time of the event, although the region has experienced gradual warming historically that appears to have resulted in increased landslide occurrence in other parts of Alaska during recent years. Inherent susceptibility of geological formations to landsliding was not apparent with the available data, although discontinuities created dip-slope conditions for the five largest slides. Historical thinning of glaciers and consequent slope debuttressing may have been a factor in aiding occurrence of some of the earthquake-induced landslides, particularly some of the largest. Mobility of the largest avalanches was above average compared to global data, and mobility of all sizes of avalanches was apparently aided by movement over glacial ice. Avalanche deposits displayed characteristics indicative of turbulent flow on steeper slopes and laminar plug flow in flatter areas, where sliding also likely occurred.

Alaska↗

Numerical simulations of large earthquakes: Dynamic rupture propagation on heterogeneous faults

Our current conceptions of earthquake rupture dynamics, especially for large earthquakes, require knowledge of the geometry of the faults involved in the rupture, the material properties of the rocks surrounding the faults, the initial state of stress on the faults, and a constitutive formulation that determines when the faults can slip. In numerical simulations each of these factors appears to play a significant role in rupture propagation, at the kilometer length scale. Observational evidence of the earth indicates that at least the first three of the elements, geometry, material, and stress, can vary over many scale dimensions. Future research on earthquake rupture dynamics needs to consider at which length scales these features are significant in affecting rupture propagation. ?? Birkha??user Verlag, Basel, 2004.

Pure and Applied Geophysics↗

Three Mw ≥ 4.7 earthquakes within the Changning (China) shale gas field ruptured shallow faults intersecting with hydraulic fracturing wells

From 2017 to 2019, three destructive earthquakes (27 January 2017 Mw 4.7, 16 December 2018 Mw 5.2, and 3 January 2019 Mw 4.8) occurred in the Changning shale gas field in the southwest Sichuan Basin, China. Previous seismological studies attributed these events to hydraulic fracturing (HF), but were unable to identify the causative seismogenic faults and their slip behaviors. Here, we use Sentinel-1 synthetic aperture radar data to measure surface deformation triggered by the three events and conduct geodetic inversions to characterize their rupture models. The resulting coseismic interferograms show prominent surface deformation with the maximum line-of-sight displacements of up to 4 cm. The inversion results show that all three earthquakes mainly ruptured sedimentary formations above the shale gas bed, in the upper 3 km of the crust, with slip magnitudes ranging from 8.5 to 15 cm, and stress drops ranging from ∼1.8 to ∼3.3 MPa. Their source faults intersect with horizontal HF wells, but do not root in the crystalline basement. Combined with the reported difficulty of increasing HF operation pressures prior to the three events, we argue that they were most likely induced by direct injection of pressurized fluids into the fault zones. Crustal deformation patterns inferred from regional topography and GPS velocities highlight that the Changning field is located within a triple junction region near the southeastern margin of the Tibetan Plateau with large deformation gradients; such conditions are not only favorable to the development of critically stressed faults, but also facilitate the occurrence of at least moderate magnitude earthquakes.

Changning shale gas field, Sichuan Basin↗

Finite-fault analysis of the 1979 March 14 Petatlan, Mexico, earthquake using teleseismic P waveforms

Vertical, teleseismic P waves recorded for the 1979 March 14 Petatlan, Mexico, earthquake were used to derive the distribution of coseismic slip using a linear finite-fault inversion scheme that solves for the amount of slip in each of a series of consecutive time windows. Data recorded by six stations of the Global Digital Seismograph Network were inverted in addition to digitized analogue long-period recordings available from nine Worldwide Standardized Seismograph Network stations. The digital data include four broad-band and short-period velocity waveforms reconstructed from the short- and long-period components. The time-window approach allows for a variable rise time on the fault and accounts for the source multiplicity evident in the recorded P waveforms. Synthetic tests conducted using the inversion method on the limited data set, however, reveal that the data are insufficient to identify the exact dislocation duration on the fault. The method is thus implemented by prescribing the fault rise time using five consecutive 1 s time windows. The coseismic slip inferred from the P waves shows a small 70 cm peak near the earthquake hypocentre and a large zone of dislocation (1.2 m maximum) further south-east. The slip pattern covers depths from 3 to 25 km and is located south-east of other recent large interplate ruptures on the Michoacan segment of the Mexican subduction zone. This result indicates that the 1979 Petatlan earthquake broke an independent, adjacent portion of the Cocos-North America plate boundary. The seismic moment of 1.5 × 10 27 dyn cm inferred from the P waves is approximately one-half the long-period moment estimated by other investigators from the observed surface waves. Although the discrepancy is within the uncertainty of the seismic-moment estimates, it may suggest the presence of a component of slow interplate motion that did not radiate significant P -wave energy.

Geophysical Journal International↗

Earthquake swarm along the San Andreas fault near Palmdale, Southern California, 1976 to 1977

Between November 1976 and November 1977 a swarm of small earthquakes (local magnitude ≤ 3) occurred on or near the San Andreas fault near Palmdale, California. This swarm was the first observed along this section of the San Andreas since cataloging of instrumental data began in 1932. The activity followed partial subsidence of the 35-centimeter vertical crustal uplift known as the Palmdale bulge along this "locked" section of the San Andreas, which last broke in the great (surface-wave magnitude = 8¼;+) 1857 Fort Tejon earthquake. The swarm events exhibit characteristics previously observed for some foreshock sequences, such as tight clustering of hypocenters and time-dependent rotations of stress axes inferred from focal mechanisms. However, because of our present lack of understanding of the processes that precede earthquake faulting, the implications of the swarm for future large earthquakes on the San Andreas fault are unknown.

California↗

Rupture models for the A.D. 900-930 Seattle fault earthquake from uplifted shorelines

A major earthquake on the Seattle fault, Washington, ca. A.D. 900-930 was first inferred from uplifted shorelines and tsunami deposits. Despite follow-up geophysical and geological investigations, the rupture parameters of the earthquake and the geometry of the fault are uncertain. Here we estimate the fault geometry, slip direction, and magnitude of the earthquake by modeling shoreline elevation change. The best fitting model geometry is a reverse fault with a shallow roof ramp consisting of at least two back thrusts. The best fitting rupture is a SW-NE ohlique reverse slip with horizontal shortening of 15 m, rupture depth of 12.5 km, and magnitude Mw = 7.5. ?? 2006 Geological Society of America.

Washington↗

Deformed submarine terraces in Puget Sound, Pacific Northwest, indicate only one M >~7.5 earthquake on the Seattle fault zone in the past 11,000 yr

Submerged marine terraces in Puget Sound, deformed across the Seattle fault zone (SFZ), indicate that only one earthquake as large as M~7.5 has occurred in at least the past 11 kyr. Previous paleoseismic studies document a M~7.5 earthquake between 923–4 CE, which uplifted coastal marine terraces by as much as 8 m. We demonstrate that this earthquake was the only such event since ~11 ka by mapping and quantifying deformation of older marine terraces, now submerged in Puget Sound. The submerged terraces, attributed to a late-glacial sea-level lowstand, record both glacial isostatic rebound and tectonic deformation. Vertical offset of the ~11 ka terraces within the SFZ is comparable to that of the marine terraces uplifted in 923 CE, implying no additional large (M>~7.5) earthquake on the SFZ since ~11 ka. This result implies a longer recurrence interval than current hazard estimates, which assumes recurrence of M>7.1 events every 5 kyr. Our mapping of SFZ deformation since ~11 ka also supports fault segmentation and contiguous block uplift between the Seattle and Tacoma fault zones.

Washington↗