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

USGS1999

Mapping the rupture process of moderate earthquakes by inverting accelerograms

We present a waveform inversion method that uses recordings of small events as Green's functions to map the rupture growth of moderate earthquakes. The method fits P and S waveforms from many stations simultaneously in an iterative procedure to estimate the subevent rupture time and amplitude relative to the Green's function event. We invert the accelerograms written by two moderate Parkfield earthquakes using smaller events as Green's functions. The first earthquake ( M = 4.6) occurred on November 14, 1993, at a depth of 11 km under Middle Mountain, in the assumed preparation zone for the next Parkfield main shock. The second earthquake ( M = 4.7) occurred on December 20, 1994, some 6 km to the southeast, at a depth of 9 km on a section of the San Andreas fault with no previous microseismicity and little inferred coseismic slip in the 1966 Parkfield earthquake. The inversion results are strikingly different for the two events. The average stress release in the 1993 event was 50 bars, distributed over a geometrically complex area of 0.9 km 2 . The average stress release in the 1994 event was only 6 bars, distributed over a roughly elliptical area of 20 km 2 . The ruptures of both events appear to grow spasmodically into relatively complex shapes: the inversion only constrains the ruptures to grow more slowly than the S wave velocity but does not use smoothness constraints.

Journal of Geophysical Research B: Solid Earth↗

USGS1996

Earthquake prediction: The interaction of public policy and science

Earthquake prediction research has searched for both informational phenomena, those that provide information about earthquake hazards useful to the public, and causal phenomena, causally related to the physical processes governing failure on a fault, to improve our understanding of those processes. Neither informational nor causal phenomena are a subset of the other. I propose a classification of potential earthquake predictors of informational, causal, and predictive phenomena, where predictors are causal phenomena that provide more accurate assessments of the earthquake hazard than can be gotten from assuming a random distribution. Achieving higher, more accurate probabilities than a random distribution requires much more information about the precursor than just that it is causally related to the earthquake.

Conference Paper↗

USGS2017

Accounting for orphaned aftershocks in the earthquake background rate

Aftershocks often occur within cascades of triggered seismicity in which each generation of aftershocks triggers an additional generation, and so on. The rate of earthquakes in any particular generation follows Omori's law, going approximately as 1/ t . This function decays rapidly, but is heavy-tailed, and aftershock sequences may persist for long times at a rate that is difficult to discriminate from background. It is likely that some apparently spontaneous earthquakes in the observational catalogue are orphaned aftershocks of long-past main shocks. To assess the relative proportion of orphaned aftershocks in the apparent background rate, I develop an extension of the ETAS model that explicitly includes the expected contribution of orphaned aftershocks to the apparent background rate. Applying this model to California, I find that the apparent background rate can be almost entirely attributed to orphaned aftershocks, depending on the assumed duration of an aftershock sequence. This implies an earthquake cascade with a branching ratio (the average number of directly triggered aftershocks per main shock) of nearly unity. In physical terms, this implies that very few earthquakes are completely isolated from the perturbing effects of other earthquakes within the fault system. Accounting for orphaned aftershocks in the ETAS model gives more accurate estimates of the true background rate, and more realistic expectations for long-term seismicity patterns.

Geophysical Journal International↗

USGS2015

Crowdsourced earthquake early warning

Earthquake early warning (EEW) can reduce harm to people and infrastructure from earthquakes and tsunamis, but it has not been implemented in most high earthquake-risk regions because of prohibitive cost. Common consumer devices such as smartphones contain low-cost versions of the sensors used in EEW. Although less accurate than scientific-grade instruments, these sensors are globally ubiquitous. Through controlled tests of consumer devices, simulation of an M w (moment magnitude) 7 earthquake on California’s Hayward fault, and real data from the M w 9 Tohoku-oki earthquake, we demonstrate that EEW could be achieved via crowdsourcing.

Science Advances↗

USGS2017

Nucleation speed limit on remote fluid induced earthquakes

Earthquakes triggered by other remote seismic events are explained as a response to long-traveling seismic waves that temporarily stress the crust. However, delays of hours or days after seismic waves pass through are reported by several studies, which are difficult to reconcile with the transient stresses imparted by seismic waves. We show that these delays are proportional to magnitude and that nucleation times are best fit to a fluid diffusion process if the governing rupture process involves unlocking a magnitude-dependent critical nucleation zone. It is well established that distant earthquakes can strongly affect the pressure and distribution of crustal pore fluids. Earth’s crust contains hydraulically isolated, pressurized compartments in which fluids are contained within low-permeability walls. We know that strong shaking induced by seismic waves from large earthquakes can change the permeability of rocks. Thus, the boundary of a pressurized compartment may see its permeability rise. Previously confined, overpressurized pore fluids may then diffuse away, infiltrate faults, decrease their strength, and induce earthquakes. Magnitude-dependent delays and critical nucleation zone conclusions can also be applied to human-induced earthquakes.

Science Advances↗

USGS2025

Rapid fault healing from cementation controls the dynamics of deep slow slip and tremor

Despite its status as one of the most important discoveries in geophysics, the physical mechanism(s) responsible for slow slip events (SSEs) are not well understood. Here, we synthesize observations of deep SSEs in the Cascadia Subduction Zone and argue that rapid, cohesive fault strengthening may control the dynamics of deep SSEs. Cohesive strength is frequently ignored in constitutive laws used to describe fault rheology in numerical simulations of earthquakes and SSEs alike. To demonstrate its importance, we perform and analyze a suite of petrological experiments that simulate fault healing under representative pressure and temperature conditions. We show that significant cohesive strength recovery caused by dissolution-precipitation processes occurs on timescales of just a few hours. Together, our experimental and observational results support the idea that cohesion is a key component of fault strength under SSE conditions and highlight the need for its inclusion in both future experiments and numerical models of fault slip.

British Columbia, Washington↗

USGS2001

Steep-dip seismic imaging of the shallow San Andreas Fault near Parkfield

Seismic reflection and refraction images illuminate the San Andreas Fault to a depth of 1 kilometer. The prestack depth-migrated reflection image contains near-vertical reflections aligned with the active fault trace. The fault is vertical in the upper 0.5 kilometer, then dips about 70° to the southwest to at least 1 kilometer subsurface. This dip reconciles the difference between the computed locations of earthquakes and the surface fault trace. The seismic velocity cross section shows strong lateral variations. Relatively low velocity (10 to 30%), high electrical conductivity, and low density indicate a 1-kilometer-wide vertical wedge of porous sediment or fractured rock immediately southwest of the active fault trace.

California↗

USGS2015

Stress orientations in subduction zones and the strength of subduction megathrust faults

Subduction zone megathrust faults produce most of the world’s largest earthquakes. Although the physical properties of these faults are difficult to observe directly, their frictional strength can be estimated indirectly by constraining the orientations of the stresses that act on them. A global investigation of stress orientations in subduction zones finds that the maximum compressive stress axis plunges systematically trenchward, consistently making a 45°-60° angle to the subduction megathrust fault. These angles indicate that the megathrust fault is not substantially weaker than its surroundings. Together with several other lines of evidence, this implies that subduction zone megathrusts are weak faults in a low-stress environment. The deforming outer accretionary wedge may decouple the stress state along the megathrust from the constraints of the free surface.

Science↗

USGS1978

Mapped offset on the right-lateral Kern Canyon fault, southern Sierra Nevada, California

The north-trending Kern Canyon fault, the longest fault within the southern Sierra Nevada, has been mapped from lat 36°00′N to its northern end near lat 36°40′N. The fault cuts and offsets granitic plutons as young as 80 m.y., but despite the fact that many recent earthquake foci plot close to the fault, it does not appear to offset an overlying 3.5-m.y.-old basalt flow. Seven granitic plutons are clearly offset by the fault in a right-lateral sense. In the area mapped, offset of plutonic contacts is 6.5 to 13 km and increases southward by 0.2 km/km.

California↗

USGS1996

Sand boils induced by the 1993 Mississippi River flood: Could they one day be misinterpreted as earthquake-induced liquefaction?

In areas that are seismically active but lacking clear surficial faulting, many paleoearthquake studies depend on the interpretation of ancient liquefaction features (sand blows) as indicators of prehistoric seismicity. Sand blows, however, can be mimicked by nonseismic sand boils formed by water seeping beneath levees during floods. We examined sand boils induced by the Mississippi River flood of 1993 in order to compare their characteristics with sand blows of the New Madrid earthquakes of 1811–1812. We found a number of criteria that allow a distinction between the two types of deposits. (1) Earthquake-induced liquefaction deposits are broadly distributed about an epicentral area, whereas flood-induced sand boils are limited to a narrow band along a river's levee. (2) The conduits of most earthquake-induced sand blows are planar dikes, whereas the conduits of flood-induced sand boils are most commonly tubular. (3) Depression of the preearthquake ground surface is usual for sand blows, not for sand boils. (4) Flood-induced sand boils tend to be better sorted and much finer than sand-blow deposits. (5) Source beds for earthquake-induced deposits occur at a wide range of depths, whereas the source bed for sand boils is always near surface. (6) Materials removed from the walls surrounding the vent of a sand blow are seen inside sand blows, but are rarely seen inside sand boils. In general, flood-induced sand boils examined are interpreted to represent a less-energetic genesis than earthquake-induced liquefaction.

Geology↗

USGSDate not supplied

Remnant damage from the 1906 San Francisco earthquake

This field trip consists of two stops at locations where it is possible to see damage from the 1906 earthquake and to gauge the intensity of the ground shaking that caused the damage. The first stop is at a cemetery in Colma, where the damage to monuments and headstones was photographed and roughly quantified in the Report of the State Earthquake Investigation Commission , Lawson (1908), commonly referred to as the “Lawson Report.” The Lawson Report represents the formal study of the earthquake and consists of a compilation of the reports of many investigators who gathered information about faulting, ground failure, and damage due to the 1906 earthquake. The second stop is at a brick office building at the southern limit of San Francisco that was damaged by the earthquake but repaired in such a fashion that the damage is still clearly evident.

California↗

USGS2013

Mw 8.6 Sumatran earthquake of 11 April 2012: rare seaward expression of oblique subduction

The magnitude 8.6 and 8.2 earthquakes off northwestern Sumatra on 11 April 2012 generated small tsunami waves that were recorded by stations around the Indian Ocean. Combining differential travel-time modeling of tsunami waves with results from back projection of seismic data reveals a complex source with a significant trench-parallel component. The oblique plate convergence indicates that ~20-50 m of trench-parallel displacement could have accumulated since the last megathrust earthquake, only part of which has been taken up by the Great Sumatran fault. This suggests that the remaining trench-parallel motion was released during the magnitude 8.6 earthquake on 11 April 2012 within the subducting plate. The magnitude 8.6 earthquake is interpreted to be a result of oblique subduction as well as a reduction in normal stress due to the occurrence of the Sumatra-Andaman earthquake in 2004.

Geology↗

USGS2013

Intermediate-depth earthquakes facilitated by eclogitization-related stresses

Eclogitization of the basaltic and gabbroic layer in the oceanic crust involves a volume reduction of 10%–15%. One consequence of the negative volume change is the formation of a paired stress field as a result of strain compatibility across the reaction front. Here we use waveform analysis of a tiny seismic cluster in the lower crust of the downgoing Pacific plate and reveal new evidence in favor of this mechanism: tensional earthquakes lying 1 km above compressional earthquakes, and earthquakes with highly similar waveforms lying on well-defined planes with complementary rupture areas. The tensional stress is interpreted to be caused by the dimensional mismatch between crust transformed to eclogite and underlying untransformed crust, and the earthquakes are probably facilitated by reactivation of fossil faults extant in the subducting plate. These observations provide seismic evidence for the role of volume change–related stresses and, possibly, fluid-related embrittlement as viable processes for nucleating earthquakes in downgoing oceanic lithosphere.

Geology↗

USGS2025

Evidence for nonlocal sediment transport on hillslopes from fault scarp morphology

Hillslope sediment transport processes such as bioturbation, rainsplash, and granular mechanics occur across the entire planet. Yet, it remains uncertain how these small-scale processes act together to shape landscapes. Longstanding hillslope diffusion theory posits that hillslope processes are spatially limited, whereas new concepts of nonlocal sediment transport argue otherwise. However, each theory produces subtly different, but distinct, predictions for the evolution of fault scarps. We use the topographic change of fault scarps to demonstrate that hillslope processes produce nonlocal sediment transport. Analysis of a global compilation of 340 dated single-earthquake scarp profiles reveals a statistically significant ( p < 0.05 ) relationship between scarp age and scarp asymmetry, here defined as the ratio of imaginary to real components of the Fourier transform of absolute slope. Numerical simulations show that nonlocal models predict this relationship, whereas hillslope diffusion models do not. To further investigate this result, we examined the depositional geometry of a well-exposed colluvial wedge along the Wasatch fault in central Utah, United States. Our quantitative comparison between the exposure and numerical simulations reveals better agreement with the nonlocal model. Nonlocal sediment transport theory appears to best capture the physics of how hillslope processes shape fault scarps, yet hillslope diffusion provides a useful approximation in many cases. As the processes that act on fault scarps are nearly identical to those acting on hillslopes, our results provide evidence supporting nonlocality as a generalized model of hillslope sediment transport.

Utah↗

USGS2003

Signatures of the seismic source in EMD-based characterization of the 1994 Northridge, California, earthquake recordings

In this article we use empirical mode decomposition (EMD) to characterize the 1994 Northridge, California, earthquake records and investigate the signatures carried over from the source rupture process. Comparison of the current study results with existing source inverse solutions that use traditional data processing suggests that the EMD-based characterization contains information that sheds light on aspects of the earthquake rupture process. We first summarize the fundamentals of the EMD and illustrate its features through the analysis of a hypothetical and a real record. Typically, the Northridge strong-motion records are decomposed into eight or nine intrinsic mode functions (IMF's), each of which emphasizes a different oscillation mode with different amplitude and frequency content. The first IMF has the highest-frequency content; frequency content decreases with an increase in IMF component. With the aid of a finite-fault inversion method, we then examine aspects of the source of the 1994 Northridge earthquake that are reflected in the second to fifth IMF components. This study shows that the second IMF is predominantly wave motion generated near the hypocenter, with high-frequency content that might be related to a large stress drop associated with the initiation of the earthquake. As one progresses from the second to the fifth IMF component, there is a general migration of the source region away from the hypocenter with associated longer-period signals as the rupture propagates. This study suggests that the different IMF components carry information on the earthquake rupture process that is expressed in their different frequency bands.

California↗

USGS2010

Intensity, magnitude, location and attenuation in India for felt earthquakes since 1762

A comprehensive, consistently interpreted new catalog of felt intensities for India (Martin and Szeliga, 2010, this issue) includes intensities for 570 earthquakes; instrumental magnitudes and locations are available for 100 of these events. We use the intensity values for 29 of the instrumentally recorded events to develop new intensity versus attenuation relations for the Indian subcontinent and the Himalayan region. We then use these relations to determine the locations and magnitudes of 234 historical events, using the method of Bakun and Wentworth (1997). For the remaining 336 events, intensity distributions are too sparse to determine magnitude or location. We evaluate magnitude and location accuracy of newly located events by comparing the instrumental- with the intensity-derived location for 29 calibration events, for which more than 15 intensity observations are available. With few exceptions, most intensity-derived locations lie within a fault length of the instrumentally determined location. For events in which the azimuthal distribution of intensities is limited, we conclude that the formal error bounds from the regression of Bakun and Wentworth (1997) do not reflect the true uncertainties. We also find that the regression underestimates the uncertainties of the location and magnitude of the 1819 Allah Bund earthquake, for which a location has been inferred from mapped surface deformation. Comparing our inferred attenuation relations to those developed for other regions, we find that attenuation for Himalayan events is comparable to intensity attenuation in California (Bakun and Wentworth, 1997), while intensity attenuation for cratonic events is higher than intensity attenuation reported for central/eastern North America (Bakun et al. , 2003). Further, we present evidence that intensities of intraplate earthquakes have a nonlinear dependence on magnitude such that attenuation relations based largely on small-to-moderate earthquakes may significantly overestimate the magnitudes of historical earthquakes.

Bulletin of the Seismological Society of America↗

USGS2010

Rupture directivity of moderate earthquakes in northern California

We invert peak ground velocity and acceleration (PGV and PGA) to estimate rupture direction and rupture velocity for 47 moderate earthquakes (3.5≥M≥5.4) in northern California. We correct sets of PGAs and PGVs recorded at stations less than 55–125 km, depending on source depth, for site amplification and source–receiver distance, then fit the residual peak motions to the unilateral directivity function of Ben-Menahem (1961). We independently invert PGA and PGV. The rupture direction can be determined using as few as seven peak motions if the station distribution is sufficient. The rupture velocity is unstable, however, if there are no takeoff angles within 30° of the rupture direction. Rupture velocities are generally subsonic (0.5β–0.9β); for stability, we limit the rupture velocity at v=0.92β, the Rayleigh wave speed. For 73 of 94 inversions, the rupture direction clearly identifies one of the nodal planes as the fault plane. The 35 strike-slip earthquakes have rupture directions that range from nearly horizontal (6 events) to directly updip (5 events); the other 24 rupture partly along strike and partly updip. Two strike-slip earthquakes rupture updip in one inversion and downdip in the other. All but 1 of the 11 thrust earthquakes rupture predominantly updip. We compare the rupture directions for 10 M≥4.0 earthquakes to the relative location of the mainshock and the first two weeks of aftershocks. Spatial distributions of 8 of 10 aftershock sequences agree well with the rupture directivity calculated for the mainshock.

California↗

USGS2012

A global earthquake discrimination scheme to optimize ground-motion prediction equation selection

We present a new automatic earthquake discrimination procedure to determine in near-real time the tectonic regime and seismotectonic domain of an earthquake, its most likely source type, and the corresponding ground-motion prediction equation (GMPE) class to be used in the U.S. Geological Survey (USGS) Global ShakeMap system. This method makes use of the Flinn–Engdahl regionalization scheme, seismotectonic information (plate boundaries, global geology, seismicity catalogs, and regional and local studies), and the source parameters available from the USGS National Earthquake Information Center in the minutes following an earthquake to give the best estimation of the setting and mechanism of the event. Depending on the tectonic setting, additional criteria based on hypocentral depth, style of faulting, and regional seismicity may be applied. For subduction zones, these criteria include the use of focal mechanism information and detailed interface models to discriminate among outer-rise, upper-plate, interface, and intraslab seismicity. The scheme is validated against a large database of recent historical earthquakes. Though developed to assess GMPE selection in Global ShakeMap operations, we anticipate a variety of uses for this strategy, from real-time processing systems to any analysis involving tectonic classification of sources from seismic catalogs.

Bulletin of the Seismological Society of America↗

Refine your search to explore more results.