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 325 records · Page 18Linked to original sources

Faults, lineaments, and earthquake epicenters digital map of the Pahute Mesa 30' x 60' quadrangle, Nevada

This map database, identified as Faults, lineaments, and earthquake epicenters digital map of the Pahute Mesa 30' X 60' quadrangle, Nevada, has been approved for release and publication by the Director of the USGS. Although this database has been subjected to rigorous review and is substantially complete, the USGS reserves the right to revise the data pursuant to further analysis and review. Furthermore, it is released on condition that neither the USGS nor the United States Government may be held liable for any damages resulting from its authorized or unauthorized use. This digital map compilation incorporates fault, air photo lineament, and earthquake epicenter data from within the Pahute Mesa 30' by 60' quadrangle, southern Nye County, Nevada (fig. 1). The compilation contributes to the U.S. Department of Energy's Yucca Mountain Project, established to determine whether or not the Yucca Mountain site is suitable for the disposal of high-level nuclear waste. Studies of local and regional faulting and earthquake activity, including the features depicted in this compilation, are carried out to help characterize seismic hazards and tectonic processes that may be relevant to the future stability of Yucca Mountain. The Yucca Mountain site is located in the central part of the Beatty 30' by 60' quadrangle approximately 15 km south of the south edge of the Pahute Mesa quadrangle (fig. 1). The U.S. Geological Survey participates in studies of the Yucca Mountain site under Interagency Agreement DE-AI08-78ET44802. The map compilation is only available on line as a digital database in ARC/INFO ASCII (Generate) and export formats. The database can be downloaded via 'anonymous ftp' from a USGS system named greenwood.cr.usgs.gov (136.177.48.5). The files are located in a directory named /pub/open-file-reports/ofr-96-0262. This directory contains a text document named 'README.1 ST' that contains database technical and explanatory documentation, including instructions for uncompressing the bundled (tar) file. In displaying the compilation it is important to note that the map data set is considered accurate when depicted at a scale of about 1:100,000; displaying the compilation at scales significantly larger than this may result in distortions and (or) mislocations of the data.

Nevada↗

Evidence of previous faulting along the 2019 Ridgecrest, California earthquake ruptures

The July 2019 Ridgecrest earthquake sequence in southeastern California was characterized as surprising because only ~35% of the rupture occurred on previously mapped faults. Employing more detailed inspection of pre-event high-resolution topography and imagery in combination with field observations, we document evidence of active faulting in the landscape along the entire fault system. Scarps, deflected drainages, and lineaments and contrasts in topography, vegetation, and ground color demonstrate previous slip on a dense network of orthogonal faults, consistent with patterns of surface rupture observed in 2019. Not all of these newly mapped fault strands ruptured in 2019. Outcrop-scale field observations additionally reveal tufa lineaments and sheared Quaternary deposits. Neotectonic features are commonly short (<2 km), discontinuous, and display en echelon patterns along both the M 6.4 and M 7.1 ruptures. These features are generally more prominent and better preserved outside the late Pleistocene lake basins. Fault expression may also be related to deformation style: scarps and topographic lineaments are more prevalent in areas where substantial vertical motion occurred in 2019. Where strike-slip displacement dominated in 2019, the faults are mainly expressed by less prominent tonal and vegetation features. Both the NE- and NW-trending active fault systems are subparallel to regional bedrock fabrics that were established as early as ~150 Ma, and may be reactivating these older structures. Overall, we estimate that 50-70% (i.e., an additional 15-35%) of the 2019 surface ruptures could have been recognized as active faults with detailed inspection of pre-event data. Similar detailed mapping of potential neotectonic features could help improve seismic hazard analyses in other regions of eastern California and elsewhere that have distributed faulting or incompletely mapped faults. In areas where faults cannot be resolved as single thoroughgoing structures, a zone of potential faulting should be used as a hazard model input.

California↗

Reactivated faulting near Cushing, Oklahoma: Increased potential for a triggered earthquake in an area of United States strategic infrastructure

In October 2014 two moderate-sized earthquakes ( Mw 4.0 and 4.3) struck south of Cushing, Oklahoma, below the largest crude oil storage facility in the world. Combined analysis of the spatial distribution of earthquakes and regional moment tensor focal mechanisms indicate reactivation of a subsurface unnamed and unmapped left-lateral strike-slip fault. Coulomb failure stress change calculations using the relocated seismicity and slip distribution determined from regional moment tensors, allow for the possibility that the Wilzetta-Whitetail fault zone south of Cushing, Oklahoma, could produce a large, damaging earthquake comparable to the 2011 Prague event. Resultant very strong shaking levels (MMI VII) in the epicentral region present the possibility of this potential earthquake causing moderate to heavy damage to national strategic infrastructure and local communities.

Oklahoma↗

Spatial variations in fault friction related to lithology from rupture and afterslip of the 2014 South Napa, California, earthquake

Following earthquakes, faults are often observed to continue slipping aseismically. It has been proposed that this afterslip occurs on parts of the fault with rate-strengthening friction that are stressed by the mainshock, but our understanding has been limited by a lack of immediate, high-resolution observations. Here we show that the behavior of afterslip following the 2014 South Napa earthquake varied over distances of only a few kilometers. This variability cannot be explained by coseismic stress changes alone. We present daily positions from continuous and survey GPS sites that we re-measured within 12 hours of the mainshock, and surface displacements from the new Sentinel-1 radar mission. This unique geodetic data set constrains the distribution and evolution of coseismic and postseismic fault slip with exceptional resolution in space and time. We suggest that the observed heterogeneity in behavior is caused by lithological controls on the frictional properties of the fault plane.

California↗

Use of microearthquakes in the study of the mechanics of earthquake generation along the San Andreas fault in central California

A small, dense network of independently recording portable seismograph stations was used to delineate the slip surface associated with the 1966 Parkfield-Cholame earthquake by precise three dimensional mapping of the hypocenters of its aftershocks. The aftershocks were concentrated in a very narrow vertical zone beneath or immediately adjacent to the zone of surf ace fracturing that accompanied the main shock. Focal depths ranged from less than 1 km to a maximum of 15 km. The same type of portable network was used to study microearthquakes associated with an actively creeping section of the San Andreas fault south of Hollister during the summer of 1967. Microearthquake activity during the 6-week operation of this network was dominated by aftershocks of a magnitude-4 earthquake that occurred within the network near Bear Valley on July 23. Most of the aftershocks were concentrated in an equidimensional region about 2 1 2km across that contained the hypocenter of the main shock. The zone of the concentrated aftershocks was centered near the middle of the rift zone at a depth of about 3 1 2km. Hypocenters of other aftershocks outlined a 25 km long zone of activity beneath the actively creeping strand of the fault and extending from the surface to a depth of about 13 km. A continuing study of microearthquakes along the San Andreas, Hayward, and Calaveras faults between Hollister and San Francisco has been under way for about 2 years. The permanent telemetered network constructed for this purpose has grown from about 30 stations in early 1968 to about 45 stations in late 1969. Microearthquakes between Hollister and San Francisco are heavily concentrated in narrow, nearly vertical zones along sections of the Sargent, San Andreas, and Calaveras faults. Focal depths range from less than 1 km to about 14 km.

California↗

Earthquake recurrence on the Calaveras fault east of San Jose, California

Occurrence of small (3 ⩽ M L < 4) earthquakes on two 10-km segments of the Calaveras fault between Calaveras and Anderson reservoirs follows a simple linear pattern of elastic strain accumulation and release. The centers of these independent patches of earthquake activity are 20 km apart. Each region is characterized by a constant rate of seismic slip as computed from earthquake magnitudes, and is assumed to be an isolated locked patch on a creeping fault surface. By calculating seismic slip rates and the amount of seismic slip since the time of the last significant ( M ⩾ 3) earthquake, it is possible to estimate the most likely date of the next ( M ⩾- 3) event on each patch. The larger the last significant event, the longer the time until the next one. The recurrence time also appears to be increased according to the moment of smaller (2 < M L < 3) events in the interim. The anticipated times of future larger events on each patch, on the basis of preliminary location data through May 1977 and estimates of interim activity, are tabulated below with standard errors. The occurrence time for the southern zone is based on eight recurrent events since 1969, the northern zone on only three. The 95% confidence limits can be estimated as twice the standard error of the projected least-squares line. Events of M ⩾ 3 should not occur in the specified zones at times outside these limits. The central region between the two zones was the locus of two events ( M = 3.6, 3.3) on July 3, 1977. These events occurred prior to a window based on the three point, post-1969 slip-time line for the central region.

California↗

Evidence for surface rupture in 1868 on the Hayward fault in north Oakland and major rupturing in prehistoric earthquakes

WGCEP90 estimated the Hayward fault to have a high probability (0.45 in 30 yr) of producing a future M7 Bay Area earthquake. This was based on a generic recurrence time and an unverified segmentation model, because there were few direct observations for the southern fault and none for the northern Hayward fault. To better constrain recurrence and segmentation of the northern Hayward fault, we trenched in north Oakland. Unexpectedly, we observed evidence of surface rupture probably from the M7 1868 earthquake. This extends the limit of that surface rupture 13 km north of the segmentation boundary used in the WGCEP90 model and forces serious re-evaluation of the current two-segment paradigm. Although we found that major prehistoric ruptures have occurred here, we could not radiocarbon date them. However, the last major prehistoric event appears correlative with a recently recognized event 13 km to the north dated AD 1640–1776.

Geophysical Research Letters↗

Earthquake recurrence on the south Hayward fault is most consistent with a time dependent, renewal process

Elastic rebound and stress renewal are important components of earthquake forecasting because if large earthquakes can be shown to be periodic, then rupture probability is time dependent. While renewal models are used in formal forecasts, it has not been possible to exclude the alternate view that repeated large earthquakes can happen in rapid succession without requiring time for stress regeneration. Here a consistency test between time dependent and time independent recurrence distributions is made using a Monte Carlo method to replicate the paleoseismic series on the south Hayward fault. Time dependent distributions with recurrence interval of 210 years and coefficient of variation of 0.6 reproduce the event series on the south Hayward 5 times more often than any exponential distribution: a highly significant difference as determined using a two-tailed Z-test for relative proportions. Therefore large Hayward fault earthquakes are quasi-periodic and are most consistent with a stress renewal process.

Geophysical Research Letters↗

Estimates of stress drop and crustal tectonic stress from the 27 February 2010 Maule, Chile, earthquake: Implications for fault strength

The great 27 February 2010 M w 8.8 earthquake off the coast of southern Chile ruptured a ∼600 km length of subduction zone. In this paper, we make two independent estimates of shear stress in the crust in the region of the Chile earthquake. First, we use a coseismic slip model constrained by geodetic observations from interferometric synthetic aperture radar (InSAR) and GPS to derive a spatially variable estimate of the change in static shear stress along the ruptured fault. Second, we use a static force balance model to constrain the crustal shear stress required to simultaneously support observed fore‐arc topography and the stress orientation indicated by the earthquake focal mechanism. This includes the derivation of a semianalytic solution for the stress field exerted by surface and Moho topography loading the crust. We find that the deviatoric stress exerted by topography is minimized in the limit when the crust is considered an incompressible elastic solid, with a Poisson ratio of 0.5, and is independent of Young's modulus. This places a strict lower bound on the critical stress state maintained by the crust supporting plastically deformed accretionary wedge topography. We estimate the coseismic shear stress change from the Maule event ranged from −6 MPa (stress increase) to 17 MPa (stress drop), with a maximum depth‐averaged crustal shear‐stress drop of 4 MPa. We separately estimate that the plate‐driving forces acting in the region, regardless of their exact mechanism, must contribute at least 27 MPa trench‐perpendicular compression and 15 MPa trench‐parallel compression. This corresponds to a depth‐averaged shear stress of at least 7 MPa. The comparable magnitude of these two independent shear stress estimates is consistent with the interpretation that the section of the megathrust fault ruptured in the Maule earthquake is weak, with the seismic cycle relieving much of the total sustained shear stress in the crust.

Maule↗

A maximum rupture model for the southern San Andreas and San Jacinto Faults California, derived from paleoseismic earthquake ages: Observations and limitations

Paleoseismic rupture histories provide spatiotemporal models of earthquake moment release needed to test numerical models and lengthen the instrumental catalog. We develop a model of the fewest and thus largest magnitude earthquakes permitted by paleoseismic data for the last 1,500 years on the southern San Andreas and San Jacinto Faults, California, USA. The largest geometric complexity appears to regulate the system: Only two ruptures break the San Gorgonio Pass region, followed by episodes of ruptures that could bridge the northern San Jacinto Fault and the San Andreas Fault. When tested against independent data on slip per event, the model produces comparable values indicating the end‐member model does not underpredict rupture rates. Rupture of >85% of the fault length in the historic period between 1800 and 1857 and the subsequent quiescence is similar to epochs of activity in the prehistoric model, suggesting that regional clustering of seismicity could be a trait of the system.

California↗

Mountain rivers reveal the earthquake hazard of geologic faults in Silicon Valley

The 1989, M w = 6.9 Loma Prieta earthquake resulted in tens of lives lost and cost California almost 3% of its gross domestic product. Despite widespread damage, the earthquake did not clearly rupture the surface, challenging the identification and characterization of these hidden hazards. Here, we show that they can be illuminated by inverting fluvial topography for slip-and moment accrual-rates—fundamental components in earthquake hazard assessments—along relief-generating geologic faults. We applied this technique to thrust faults bounding the mountains along the western side of Silicon Valley in the San Francisco Bay Area, and discovered that these structures may be capable of generating a M w = 6.9 earthquake every 250–300 years based on moment accrual rates. This method may be deployed broadly to evaluate seismic hazard in developing regions with limited geological and geophysical information.

California↗

Interpretation of changes in water level accompanying fault creep and implications for earthquake prediction

Quantitative calculations for the effect of a fault creep event on observations of changes in water level in wells provide an approach to the tectonic interpretation of these phenomena. For the pore pressure field associated with an idealized creep event having an exponential displacement versus time curve, an analytic expression has been obtained in terms of exponential-integral functions. The pore pressure versus time curves for observation points near the fault are pulselike; a sharp pressure increase (or decrease, depending on the direction of propagation) is followed by more gradual decay to the normal level after the creep event. The time function of the water level change may be obtained by applying the filter—derived by A. G. Johnson and others to determine the influence of atmospheric pressure on water level—to the analytic pore pressure versus time curves. The resulting water level curves show a fairly rapid increase (or decrease) and then a very gradual return to normal. The results of this analytic model do not reproduce the steplike changes in water level observed by Johnson and others. If the procedure used to obtain the water level from the pore pressure is correct, these results suggest that steplike changes in water level are not produced by smoothly propagating creep events but by creep events that propagate discontinously, by changes in the bulk properties of the region around the well, or by some other mechanism. In addition, simplistic calculations show that significant pressure field variations and water level changes near the surface may be expected to accompany a propagating creep event on a buried fault. Water level changes of as much as several meters may be expected at the surface for a creep event having a dislocation amplitude of l m on a semi-infinite dislocation surface extending downward from a depth of 10 km and propagating horizontally at a rate of 10 km/day. The maximum near-surface effect should be observed at a horizontal distance from the fault about equal to the depth of the top of the dislocation surface. These results are consistent with the observations made in China of large water level changes preceding large earthquakes, if some sort of aseismic creep event at depth precedes these earthquakes.

Journal of Geophysical Research Solid Earth↗

Subduction intraslab-interface fault interactions in the 2022 Mw 6.4 Ferndale, California earthquake sequence

The Mendocino triple junction, the intersection of the Pacific, North American, and Gorda plates, activates a collection of disparate faults that reconcile Cascadia subduction with San Andreas transform motion. The December 20, 2022, Mw 6.4 Ferndale, California earthquake occurred within this complex zone as strike-slip faulting within the subducting Gorda slab. Here, we analyze the seismic and geodetic signatures of the mainshock and aftershock sequence to illuminate its role within complex tectonic surroundings. We find aftershocks on varied fault structures within the uppermost Gorda slab, yet seismicity on the subduction interface itself was notably absent. Nevertheless, we identify small but coherent postseismic deformation that is well modeled by aseismic slip on this interface, likely triggered by stresses generated at the updip limit of coseismic rupture. This sequence demonstrates the potential for interactions between intra-slab earthquakes and slip on the subduction megathrust, highlighting the need to consider this and other subduction zones as coupled systems of interacting faults.

California↗

Fiber-imaged supershear dynamics in the 2024 Mw 7 Mendocino Fault earthquake

Fault structure and rupture physics are deeply intertwined, and observations of this coupling are critical for understanding earthquake behavior. Rupture propagation is observable at fine scales using dense seismic networks. Fiber-optic sensing allows for long-term deployments of ultradense arrays that enable high-resolution measurements of infrequent, large earthquakes. We recorded the 2024 moment magnitude ( M w ) 7 Mendocino Fault earthquake with a nearby fiber-optic array and imaged its behavior with seismic beamforming. The rupture propagated to the east at subshear velocity; stagnated near the Mendocino Triple Junction, a zone of structural complexity; and subsequently transitioned to supershear velocity. The correlation between source physics and structure shows how lithospheric heterogeneity affects first-order characteristics of earthquake ruptures. Our results also demonstrate the potential for fiber-optic sensing to improve real-time estimation of key parameters for early warning.

California↗

Birth of a fault: Connecting the Kern County and Walker Pass, California, earthquakes

A band of seismicity transects the southern Sierra Nevada range between the northeastern end of the site of the 1952 M W (moment magnitude) 7.3 Kern County earthquake and the site of the 1946 M W 6.1 Walker Pass earthquake. Relocated earthquakes in this band, which lacks a surface expression, better delineate the northeast-trending seismic lineament and resolve complex structure near the Walker Pass mainshock. Left-lateral earthquake focal planes are rotated counterclockwise from the strike of the seismic lineament, consistent with slip on shear fractures such as those observed in the early stages of fault development in laboratory experiments. We interpret this seismic lineament as a previously unrecognized, incipient, currently blind, strike-slip fault, a unique example of a newly forming structure.

Geology↗

Fault parameter constraints using relocated earthquakes: A validation of first-motion focal-mechanism data

We estimate the strike and dip of three California fault segments (Calaveras, Sargent, and a portion of the San Andreas near San Jaun Bautistia) based on principle component analysis of accurately located microearthquakes. We compare these fault orientations with two different first-motion focal mechanism catalogs: the Northern California Earthquake Data Center ( NCEDC ) catalog, calculated using the FPFIT algorithm (Reasenberg and Oppenheimer, 1985), and a catalog created using the HASH algorithm that tests mechanism stability relative to seismic velocity model variations and earthquake location (Hardebeck and Shearer, 2002). We assume any disagreement (misfit >30&deg; in strike, dip, or rake) indicates inaccurate focal mechanisms in the catalogs. With this assumption, we can quantify the parameters that identify the most optimally constrained focal mechanisms. For the NCEDC/FPFIT catalogs, we find that the best quantitative discriminator of quality focal mechanisms is the station distribution ratio ( STDR ) parameter, an indicator of how the stations are distributed about the focal sphere. Requiring STDR > 0.65 increases the acceptable mechanisms from 34%&ndash;37% to 63%&ndash;68%. This suggests stations should be uniformly distributed surrounding, rather than aligning, known fault traces. For the HASH catalogs, the fault plane uncertainty ( FPU ) parameter is the best discriminator, increasing the percent of acceptable mechanisms from 63%&ndash;78% to 81%&ndash;83% when FPU &le; 35&deg;. The overall higher percentage of acceptable mechanisms and the usefulness of the formal uncertainty in identifying quality mechanisms validate the HASH approach of testing for mechanism stability.

Bulletin of the Seismological Society of America↗

Calibration of PS09, PS10, and PS11 trans-Alaska pipeline system strong-motion instruments, with acceleration, velocity, and displacement records of the Denali fault earthquake, 03 November 2002

In September, 2003, the Alyeska Pipeline Service Company (APSC) and the U.S. Geological Survey (USGS) embarked on a joint effort to extract, test, and calibrate the accelerometers, amplifiers, and bandpass filters from the earthquake monitoring systems (EMS) at Pump Stations 09, 10, and 11 of the Trans-Alaska Pipeline System (TAPS). These were the three closest strong-motion seismographs to the Denali fault when it ruptured in the MW 7.9 earthquake of 03 November 2002 (22:12:41 UTC). The surface rupture is only 3.0 km from PS10 and 55.5 km from PS09 but PS11 is 124.2 km away from a small rupture splay and 126.9 km from the main trace. Here we briefly describe precision calibration results for all three instruments. Included with this report is a link to the seismograms reprocessed using these new calibrations: http://nsmp.wr.usgs.gov/data_sets/20021103_2212_taps.html Calibration information in this paper applies at the time of the Denali fault earthquake (03 November 2002), but not necessarily at other times because equipment at these stations is changed by APSC personnel at irregular intervals. In particular, the equipment at PS09, PS10, and PS11 was changed by our joint crew in September, 2003, so that we could perform these calibrations. The equipment stayed the same from at least the time of the earthquake until that retrieval, and these calibrations apply for that interval.

Open-File Report↗

Faulting, damage, and intensity in the Canyondam earthquake of May 23, 2013

On Thursday evening, May 23, 2013 (0347 May 24 UTC), a moment magnitude (M w ) = 5.7 earthquake occurred northeast of Canyondam, California. A two-person team of U.S. Geological Survey scientists went to the area to search for surface rupture and to canvass damage in the communities around Lake Almanor. While the causative fault had not been identified at the time of the field survey, surface rupture was expected to have occurred just south of Lake Almanor, approximately 2–4 kilometers south of the epicenter. No surface rupture was discovered. Felt intensity among the communities around Lake Almanor appeared to vary significantly. Lake Almanor West (LAW), Lake Almanor Country Club (LACC), and Hamilton Branch (HB) experienced Modified Mercalli Intensity (MMI) ≥7, whereas other communities around the lake experienced MMI ≤6; the maximum observed intensity was MMI 8, in LAW. Damage in the high intensity areas consisted of broken and collapsed chimneys, ruptured pipes, and some damage to foundations and to structural elements within houses. Although this shaking damage is not usually expected for an M w 5.7 earthquake, the intensities at Lake Almanor Country Club correlate with the peak ground acceleration (38 percent g) and peak ground velocity (30 centimeters per second) recorded by the California Strong Motion Instrumentation Program accelerometer located at the nearby Lake Almanor Fire Station. The intensity distribution for the three hardest hit areas (LAW, LACC, and HB) appears to increase as the azimuth from epicenter to the intensity sites approaches the fault strike. The small communities of Almanor and Prattville on the southwestern shore of Lake Almanor experienced somewhat lower intensities. The town of Canyondam experienced a lower intensity as well, despite its location up-dip of the earthquake rupture. This report contains information on the earthquake itself, the search for surface rupture, and the damage we observed and compiled from other sources.

California↗