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 1,405 records · Page 78Linked to original sources

Illumination of rheological mantle heterogeneity by the M7.2 2010 El Mayor-Cucapah earthquake

Major intracontinental strike-slip faults tend to mark boundaries between lithospheric blocks of contrasting mechanical properties along much of their length. Both crustal and mantle heterogeneities can form such boundaries, but the role of crustal versus mantle strength contrasts for localizing strain sufficiently to generate major faults remains unclear. Using the crustal velocity field observed through the Global Positioning System (GPS) in the epicentral area of the M7.2 2010 El Mayor-Cucapah earthquake, Baja California, we find that transient deformation observed after the event is anomalously small in areas of relatively high seismic velocity in the shallow upper mantle (∼50 km depth). This pattern is best explained with a laterally heterogeneous viscoelastic structure that mimics the seismic structure. The mantle of the Southern Colorado River Desert (SCRD) and Peninsular Ranges (PR), which bound the fault system to its east and west, respectively, have anomalously high viscosity and seismic velocity. We hypothesize that compared with the rest of the San Andreas fault (SAF) system to its north, the strike-slip fault system in northern Baja California is narrow because of the presence of the PR and SCRD high-viscosity regions which bound it.

Arizona, Baja California, California, Sonora↗

Injection‐induced earthquakes near Milan, Kansas, controlled by Karstic Networks

Induced earthquakes from waste disposal operations in otherwise tectonically stable regions significantly increases seismic hazard. It remains unclear why injections induce large earthquakes on non‐optimally oriented faults kilometers below the injection horizon, particularly since fluids are not injected under pressure, but rather poured, into the well as observed in the Milan, Kansas area. Here we propose a mechanism for induced earthquakes whereby the karstic lower Arbuckle provides the short‐circuit that establishes a tens of MPa stepwise fluid pressure increase within the basement upon arrival of the hydraulic connection to the free surface and ultimately induce slip on the deeper fault. We investigate this scenario through modeling and mechanical analysis and show that earthquakes near Milan are likely induced by large (and sudden) fluid pressure changes when the karst network links two previously isolated hydrological systems.

Kansas↗

Brittle faulting at elevated temperature and vanishing effective stress

If brittle fault strength depends only on friction, slip instability is discouraged at low effective normal stress, σ . Stress drop and the critical stiffness necessary for unstable sliding both vanish with σ ; small earthquakes cannot occur. Very low σ is inferred in the source region of low-frequency earthquakes (LFEs) on the San Andreas fault (SAF). Moreover, if pore pressure, p , is undrained at low σ , then instabilities are prevented at all scales. This is due to dilatant strengthening which arises due to a dependence of porosity on strain rate. Dilatant strengthening is σ -independent and dominates at low σ . Undrained p is inferred over time scales of less than a few days for the SAF LFEs. Based on experiments that measure rapid contact overgrowth between 350 and 530°C at very low σ , fault failure controlled by time-dependent cementation is invoked as an explanation for the SAF LFEs. Because this “cohesion” is σ -independent, stress drops can occur at σ = 0. If in addition cohesion exceeds any dilatant strengthening during slip, cohesion dominates strength at low σ . Dilatancy measured in prior faulting and shear experiments indicate that at all stress levels steady-state porosity depends on σ in addition to strain rate. Moreover, porosity at low σ depends elastically on the confining and differential stresses. A model with these additional pore pressure effects, friction, and time-dependent cohesion, applied to the SAF LFEs produces stress drops, slip speeds, and durations that are consistent with the observations, when the shear-induced dilatancy is not extreme.

JGR Solid Earth↗

Probability of one or more M ≥7 earthquakes in southern California in 30 years

Eight earthquakes of magnitude greater than or equal to seven have occurred in southern California in the past 200 years. If one assumes that such events are the product of a Poisson process, the probability of one or more earthquakes of magnitude seven or larger in southern California within any 30 year interval is 67% ± 23% (95% confidence interval). Because five of the eight M ≥ 7 earthquakes in southern California in the last 200 years occurred off of the San Andreas fault system, the probability of one or more M ≥ 7 earthquakes in southern California but not on the San Andreas fault system occurring within 30 years is 52% ± 27% (95% confidence interval).

Southern California↗

Images of crust beneath southern California will aid study of earthquakes and their effects

The Whittier Narrows earthquake of 1987 and the Northridge earthquake of 1991 highlighted the earthquake hazards associated with buried faults in the Los Angeles region. A more thorough knowledge of the subsurface structure of southern California is needed to reveal these and other buried faults and to aid us in understanding how the earthquake-producing machinery works in this region.

California↗

Geometry of the 1954 Fairview Peak-Dixie Valley earthquake sequence from a joint inversion of leveling and triangulation data

In 1954, four earthquakes greater than M s = 6.0 occurred within a 30-km radius and in a period of 6 months. Elevation and angle changes calculated from repeated leveling and triangulation surveys which span the coseismic period provide constraints on the fault geometries and coseismic slip of the faults which were activated. The quality of the coseismic geodetic data is assessed. Corrections are applied to the leveling data for subsidence due to groundwater withdrawal in the Fallon area, and a rod miscalibration error of 150 ± 30 ppm is isolated in leveling surveys made in 1967. The leveling and triangulation observations are then simultaneously inverted using the single value decomposition (SVD) inversion method to determine fault geometries and coseismic slip. Using SVD, it is possible to determine on which faults slip is resolvable given the data distribution. The faults are found to dip between 50° and 80° and extend to depths of 5 to 14 km. The geodetically derived slip values are generally equal to, or greater than, the maximum observed displacement along the surface scarps. Where slip is resolvable the geodetic data indicates the 1954 sequence contained a significant component of right-lateral slip. This is consistent with the N15°W trending shear zone which geodetic surveys have detected in western Nevada.

Journal of Geophysical Research B: Solid Earth↗

Earthquake triggering by transient and static deformations

Observational evidence for both static and transient near-field and far-field triggered seismicity are explained in terms of a frictional instability model, based on a single degree of freedom spring-slider system and rate- and state-dependent frictional constitutive equations. In this study a triggered earthquake is one whose failure time has been advanced by Δ t (clock advance) due to a stress perturbation. Triggering stress perturbations considered include square-wave transients and step functions, analogous to seismic waves and coseismic static stress changes, respectively. Perturbations are superimposed on a constant background stressing rate which represents the tectonic stressing rate. The normal stress is assumed to be constant. Approximate, closed-form solutions of the rate-and-state equations are derived for these triggering and background loads, building on the work of Dieterich [1992, 1994]. These solutions can be used to simulate the effects of static and transient stresses as a function of amplitude, onset time t 0 , and in the case of square waves, duration. The accuracies of the approximate closed-form solutions are also evaluated with respect to the full numerical solution and t 0 . The approximate solutions underpredict the full solutions, although the difference decreases as t 0 approaches the end of the earthquake cycle. The relationship between Δ t and t 0 differs for transient and static loads: a static stress step imposed late in the cycle causes less clock advance than an equal step imposed earlier, whereas a later applied transient causes greater clock advance than an equal one imposed earlier. For equal Δ t , transient amplitudes must be greater than static loads by factors of several tens to hundreds depending on t 0 . We show that the rate-and-state model requires that the total slip at failure is a constant, regardless of the loading history. Thus a static load applied early in the cycle, or a transient applied at any time, reduces the stress at the initiation of failure, whereas static loads that are applied sufficiently late raise it. Rate-and-state friction predictions differ markedly from those based on Coulomb failure stress changes (ΔCFS) in which Δ t equals the amplitude of the static stress change divided by the background stressing rate. The ΔCFS model assumes a stress failure threshold, while the rate-and-state equations require a slip failure threshold. The complete rate-and-state equations predict larger Δ t than the ΔCFS model does for static stress steps at small t 0 , and smaller Δ t than the ΔCFS model for stress steps at large t 0 . The ΔCFS model predicts nonzero Δ t only for transient loads that raise the stress to failure stress levels during the transient. In contrast, the rate-and-state model predicts nonzero Δ t for smaller loads, and triggered failure may occur well after the transient is finished. We consider heuristically the effects of triggering on a population of faults, as these effects might be evident in seismicity data. Triggering is manifest as an initial increase in seismicity rate that may be followed by a quiescence or by a return to the background rate. Available seismicity data are insufficient to discriminate whether triggered earthquakes are “new” or clock advanced. However, if triggering indeed results from advancing the failure time of inevitable earthquakes, then our modeling suggests that a quiescence always follows transient triggering and that the duration of increased seismicity also cannot exceed the duration of a triggering transient load. Quiescence follows static triggering only if the population of available faults is finite.

Journal of Geophysical Research B: Solid Earth↗

Bimodal distribution of creep event amplitudes on the San Andreas fault, California

EPISODIC fault creep, at several instrument sites along the San Andreas and associated faults in central California consists of a few small and large slip events per year generally superimposed on a background of gradual yielding at low rates1-3. Most of the events are aseismic, but a few minor displacement steps have occurred in association with local earthquakes 12. After removal of earthquake steps, event lists for several sites include significant numbers of small events about an order or magnitude below the typical 1-4-mm amplitude range for large events1, 3. Recent experimental rock-deformation results demonstrate that under biaxial loading some rocks show episodic slip on pre-cut surfaces9,10. It is not yet clear how the laboratory and field observations are related, but the data presented here indicate that episodic fault creep in nature may be more complex than previously realised. In light of the laboratory results, it is more important than ever to consider all the details of the field data concerning fault creep.

Nature↗

Gravitational body forces focus North American intraplate earthquakes

Earthquakes far from tectonic plate boundaries generally exploit ancient faults, but not all intraplate faults are equally active. The North American Great Plains exemplify such intraplate earthquake localization, with both natural and induced seismicity generally clustered in discrete zones. Here we use seismic velocity, gravity and topography to generate a 3D lithospheric density model of the region; subsequent finite-element modelling shows that seismicity focuses in regions of high-gravity-derived deviatoric stress. Furthermore, predicted principal stress directions generally align with those observed independently in earthquake moment tensors and borehole breakouts. Body forces therefore appear to control the state of stress and thus the location and style of intraplate earthquakes in the central United States with no influence from mantle convection or crustal weakness necessary. These results show that mapping where gravitational body forces encourage seismicity is crucial to understanding and appraising intraplate seismic hazard.

Nature Communications↗

The 1933 Long Beach Earthquake (California, USA): Ground motions and rupture scenario

We present a synoptic analysis of the ground motions from the 11 March 1933 Mw 6.4 Long Beach, California, earthquake, the largest known earthquake within the central Los Angeles Basin region. Our inferred shaking intensity pattern supports the association of the earthquake with the Newport-Inglewood fault; it further illuminates the concentration of severe damage in the town of Compton, where accounts suggest vertical ground motions exceeding 1 g . We use a broadband simulation approach to develop a rupture scenario for this earthquake, informed by the damage distribution. The predicted shaking for a 25-km-long fault matches the intensity distribution, with an indication that non-linear site response on soft sediments in some near-field regions was stronger than predicted using a simple model to account for non-linearity. Our results suggest that the concentration of damage near Compton can be explained by a combination of local site amplification, source-controlled directivity, and three-dimensional basin effects whereby energy was channeled towards the deepest part of the Los Angeles Basin.

California↗

The temporal distribution of seismic radiation during deep earthquake rupture

The time history of energy release during earthquakes illuminates the process of failure, which remains enigmatic for events deeper than about 100 kilometers. Stacks of teleseismic records from regional arrays for 122 intermediate (depths of 100 to 350 kilometers) and deep (depths of 350 to 700 kilometers) earthquakes show that the temporal pattern of short-period seismic radiation has a systematic variation with depth. On average, for intermediate depth events more radiation is released toward the beginning of the rupture than near the end, whereas for deep events radiation is released symmetrically over the duration of the event, with an abrupt beginning and end of rupture. These findings suggest a variation in the style of rupture related to decreasing fault heterogeneity with depth.The time history of energy release during earthquakes illuminates the process of failure, which remains enigmatic for events deeper than about 100 kilometers. Stacks of teleseismic records from regional arrays for 122 intermediate (depths of 100 to 350 kilometers) and deep (depths of 350 to 700 kilometers) earthquakes show that the temporal pattern of short-period seismic radiation has a systematic variation with depth. On average, for intermediate depth events more radiation is released toward the beginning of the rupture than near the end, whereas for deep events radiation is released symmetrically over the duration of the event, with an abrupt beginning and end of rupture. These findings suggest a variation in the style of rupture related to decreasing fault heterogeneity with depth.

Science↗

Geology and neotectonism in the epicentral area of the 2011 M5.8 Mineral, Virginia, earthquake

This fi eld guide covers a two-day west-to-east transect across the epicentral region of the 2011 M5.8 Mineral, Virginia, earthquake, the largest ever recorded in the Central Virginia seismic zone. The fi eld trip highlights results of recent bedrock and surficial geologic mapping in two adjoining 7.5-min quadrangles, the Ferncliff and the Pendleton, which together encompass the epicenter and most of the 2011–2012 aftershocks. Tectonic history of the region includes early Paleozoic accretion of an island arc (Ordovician Chopawamsic Formation) to Laurentia, intrusion of a granodiorite pluton (Ordovician Ellisville pluton), and formation of a post-Chopawamsic successor basin (Ordovician Quantico Formation), all accompanied by early Paleozoic regional deformation and metamorphism. Local transpressional faulting and retrograde metamorphism occurred in the late Paleozoic, followed by diabase dike intrusion and possible local normal faulting in the early Mesozoic. The overall goal of the bedrock mapping is to determine what existing geologic structures might have been reactivated during the 2011 seismic event, and surfi cial deposits along the South Anna River are being mapped in order to determine possible neotectonic uplift. In addition to bedrock and surfi cial studies, we have excavated trenches in an area that contains two late Paleozoic faults and represents the updip projection of the causative fault for the 2011 quake. The trenches reveal faulting that has offset surfi cial deposits dated as Quaternary in age, as well as numerous other brittle structures that suggest a geologically recent history of neotectonic activity.

Virginia↗

Rapid mapping of ultrafine fault zone topography with structure from motion

Structure from Motion (SfM) generates high-resolution topography and coregistered texture (color) from an unstructured set of overlapping photographs taken from varying viewpoints, overcoming many of the cost, time, and logistical limitations of Light Detection and Ranging (LiDAR) and other topographic surveying methods. This paper provides the first investigation of SfM as a tool for mapping fault zone topography in areas of sparse or low-lying vegetation. First, we present a simple, affordable SfM workflow, based on an unmanned helium balloon or motorized glider, an inexpensive camera, and semiautomated software. Second, we illustrate the system at two sites on southern California faults covered by existing airborne or terrestrial LiDAR, enabling a comparative assessment of SfM topography resolution and precision. At the first site, an ∼0.1 km 2 alluvial fan on the San Andreas fault, a colored point cloud of density mostly >700 points/m 2 and a 3 cm digital elevation model (DEM) and orthophoto were produced from 233 photos collected ∼50 m above ground level. When a few global positioning system ground control points are incorporated, closest point vertical distances to the much sparser (∼4 points/m 2 ) airborne LiDAR point cloud are mostly <3 cm. The second site spans an ∼1 km section of the 1992 Landers earthquake scarp. A colored point cloud of density mostly >530 points/m 2 and a 2 cm DEM and orthophoto were produced from 450 photos taken from ∼60 m above ground level. Closest point vertical distances to existing terrestrial LiDAR data of comparable density are mostly <6 cm. Each SfM survey took ∼2 h to complete and several hours to generate the scene topography and texture. SfM greatly facilitates the imaging of subtle geomorphic offsets related to past earthquakes as well as rapid response mapping or long-term monitoring of faulted landscapes.

California↗

1857 slip on the San Andreas fault Southeast of Cholame, California

Sieh and Jahns (1984) forecasted that the next moderate Parkfield earthquake might trigger a major earthquake along a fault segment greater than 30 km long southeast of Cholame. Their forecast assumed (1) the slip was 3-4 m in 1857 and characteristic of the segment; (2) a slip rate of 3.4 cm/yr; and (3) full strain release in earthquakes. This study represents an independent measurement of channel offsets, on 1:2400-scale low-sun aerial photographs and by field investigation, to estimate the amount of 1857 slip. Although rainfall is only moderate (30 cm/yr), few reliable offsets of less than 20 m persist here because cattle grazing and agricultural disking of soft sediments on the steep terrain greatly aggravate erosion. Reconstruction of offset geometry and size depends heavily on assumptions made about the post-1857 erosion. Most of the apparent 3- to 4-m offsets of Sieh and Jahns (1984) can also be measured as 2 to 3 m larger with equal or lower uncertainty. The four offsets judged as most reliable range between 5.4 and 6.7 m, and the 11 offsets of medium-high reliability average 5.8 ?? 0.3 m. Data are too sparse and ambiguous to resolve details of the 1857 slip for this segment but it is distinctly less than the 9 m of the Carrizo Plain and more than the 3-4 m previously estimated. Further trenching may refine some measurements, but probability calculations for a Cholame segment earthquake must allow for large observer-dependent uncertainty in the 1857 slip. Although the probability of an M ???7 Cholame event seems less than that suggested by a 3.5-m characteristic earthquake model, it remains among the highest in the state.

Bulletin of the Seismological Society of America↗

Transient stress-coupling between the 1992 Landers and 1999 Hector Mine, California, earthquakes

A three-dimensional finite-element model (FEM) of the Mojave block region in southern California is constructed to investigate transient stress-coupling between the 1992 Landers and 1999 Hector Mine earthquakes. The FEM simulates a poroelastic upper-crust layer coupled to a viscoelastic lower-crust layer, which is decoupled from the upper mantle. FEM predictions of the transient mechanical behavior of the crust are constrained by global positioning system (GPS) data, interferometric synthetic aperture radar (InSAR) images, fluid-pressure data from water wells, and the dislocation source of the 1999 Hector Mine earthquake. Two time-dependent parameters, hydraulic diffusivity of the upper crust and viscosity of the lower crust, are calibrated to 10 –2 m 2 ·sec –1 and 5 × 10 18 Pa·sec respectively. The hydraulic diffusivity is relatively insensitive to heterogeneous fault-zone permeability specifications and fluid-flow boundary conditions along the elastic free-surface at the top of the problem domain. The calibrated FEM is used to predict the evolution of Coulomb stress during the interval separating the 1992 Landers and 1999 Hector Mine earthquakes. The predicted change in Coulomb stress near the hypocenter of the Hector Mine earthquake increases from 0.02 to 0.05 MPa during the 7-yr interval separating the two events. This increase is primarily attributed to the recovery of decreased excess fluid pressure from the 1992 Landers coseismic (undrained) strain field. Coulomb stress predictions are insensitive to small variations of fault-plane dip and hypocentral depth estimations of the Hector Mine rupture.

Bulletin of the Seismological Society of America↗

The 2008 Wells, Nevada earthquake sequence: Source constraints using calibrated multiple event relocation and InSAR

The 2008 Wells, NV earthquake represents the largest domestic event in the conterminous U.S. outside of California since the October 1983 Borah Peak earthquake in southern Idaho. We present an improved catalog, magnitude complete to 1.6, of the foreshock-aftershock sequence, supplementing the current U.S. Geological Survey (USGS) Preliminary Determination of Epicenters (PDE) catalog with 1,928 well-located events. In order to create this catalog, both subspace and kurtosis detectors are used to obtain an initial set of earthquakes and associated locations. The latter are then calibrated through the implementation of the hypocentroidal decomposition method and relocated using the BayesLoc relocation technique. We additionally perform a finite fault slip analysis of the mainshock using InSAR observations. By combining the relocated sequence with the finite fault analysis, we show that the aftershocks occur primarily updip and along the southwestern edge of the zone of maximum slip. The aftershock locations illuminate areas of post-mainshock strain increase; aftershock depths, ranging from 5 to 16 km, are consistent with InSAR imaging, which shows that the Wells earthquake was a buried source with no observable near-surface offset.

Nevada↗

Kinematic rupture and 3D wave propagation simulations of the 2019 Mw 7.1 Ridgecrest, California, earthquake

We model the kinematic rupture process of the 2019 M w "> M w Mw 7.1 Ridgecrest, California, earthquake using numerical simulations to reproduce the elastodynamic wave field observed by inertial seismometers, high‐rate Global Navigation Satellite System stations, and borehole strainmeters. This was the largest earthquake in Southern California in 20 yr and was widely felt throughout the region. The M w "> M w Mw 7.1 mainshock was part of a large sequence of &#x223C; 30 , 000 "> ∼ 30 , 000 ∼30,000 aftershocks and was notably preceded by an M w "> M w Mw 6.4 foreshock by 34 hr on fault structures that were once poorly understood. A large number of seismic and geodetic instruments measured the rupture process for both events, with many stations located in the near field. Hence, this is a rare opportunity to better understand complex earthquake processes that arise in an immature fault zone using advanced computing. Of the kinematic rupture models that we tested, our preferred is the simplest one that reproduces signals recorded by the three different geophysical datasets; it is composed of four distinct ruptures that progressively migrate to the southeast with delayed initiation times, and typical rupture speeds. This type of model does a better job at matching the recorded ground motions and deformations than does one composed of a continuous rupture with very low‐rupture velocity, as proposed in other studies of this earthquake.

California↗

Enumerating plausible multifault ruptures in complex fault systems with physical constraints

We propose a new model for determining the set of plausible multifault ruptures in an interconnected fault system. We improve upon the rules used in the Third Uniform California Earthquake Rupture Forecast (UCERF3) to increase connectivity and the physical consistency of ruptures. We replace UCERF3’s simple azimuth change rules with new Coulomb favorability metrics and increase the maximum jump distance to 15 km. Although the UCERF3 rules were appropriate for faults with similar rakes, the Coulomb calculations used here inherently encode preferred orientations between faults with different rakes. Our new rules are designed to be insensitive to discretization details and are generally more permissive than their UCERF3 counterparts; they allow more than twice the connectivity compared to UCERF3, yet heavily penalize long ruptures that take multiple improbable jumps. The set of all possible multifault ruptures in the California fault system is near-infinite, but our model produces a tractable set of 326,707 ruptures (a modest 29% increase over UCERF3, despite the greatly increased connectivity). Inclusion in the rupture set does not dictate that a rupture receives a significant rate in the final model; rupture rates are subsequently determined by data constraints used in an inversion. We describe the rupture building algorithm and its components in detail and provide comparisons with ruptures generated by a physics-based multicycle earthquake simulator. We find that greater than twice as many ruptures generated by the simulator violate the UCERF3 rules than violate our proposed model.

California↗