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

Tectonic framework and Late Cenozoic tectonic history of the northern part of Cyprus: Implications for earthquake hazards and regional tectonics

Located near the triple junction of the African, Arabia, and Eurasian Plate, Cyprus has had an active and complex neotectonic history, which includes devastating historical earthquakes. Investigations into the tectonic framework of the northern part of Cyprus provide important insights into regional tectonism of the Eastern Mediterranean and Middle East. The northern part of Cyprus is divided into two tectono-stratigraphic terranes, the boundary of which is the Ovgos fault zone. In the Middle Miocene, the Ovgos fault zone was a marine platform margin, which separated open-marine platform carbonates from deep-marine turbidites. Transpressive movement along easterly and northeasterly trending structures dominated the Late Miocene; deposition and preservation of Messinian evaporites occurred in grabens at intersections of these trends. N–S compression began in the Early Pliocene and produced contractional tectonism along east–west trends, including major thrusting of allochthonous rocks in the Kyrenia Range. Quaternary deformation has been dominated by strike-slip faults along northeast and northwest trends; movement during the Pleistocene occurred on several of these faults; Holocene movement is documented on one of these faults. A seismic hazard is implied for the Nicosia area because its proximity to the Quaternary faults. Since the Miocene, Cyprus has been continuously uplifted, but the tectonic setting is controversial, as some researchers invoke a subduction zone setting and others ascribe to a regime of strike-slip tectonics. Our neotectonic framework is consistent with a restraining bend model for Cyprus in a regional strike-slip regime.

Journal of Asian Earth Sciences↗

The 1923 Kanto earthquake reevaluated using a newly augmented geodetic data set

This study revisits the mechanism of the 1923 Ms = 7.9 Kanto earthquake in Japan. We derive a new source model and use it to assess quantitative and qualitative aspects of the accommodation of plate motion in the Kanto region. We use a new geodetic data set that consists of displacements from leveling and angle changes from triangulation measurements obtained in surveys between 1883 and 1927. Two unique aspects of our analysis are the inclusion of a large number of second-order triangulation measurements and the application of a correction to remove interseismic deformation. The geometry of the fault planes is adopted from a recent seismic reflection study of the Kanto region. We evaluate the minimum complexity necessary in the model to fit the data optimally. Our final uniform-slip elastic dislocation model consists of two adjacent ???20?? dipping low-angle planes accommodating reverse dextral slip of 6.0 in on the larger, eastern plane and 9.5 m on the smaller, western plane with azimuths of 163?? and 121??, respectively. The earthquake was located in the Sagami trough, where the Philippine Sea plate subducts under Honshu. Compared to the highly oblique angle of plate convergence, the coseismic slip on the large fault plane has a more orthogonal orientation to the strike of the plate boundary, suggesting that slip partitioning plays a role in accommodation of plate motion. What other structure is involved in the partitioning is unclear. Uplift records of marine coastal terraces in Sagami Bay document 7500 years of earthquake activity and predict average recurrence intervals of 400 years for events with vertical displacement profiles similar to those of the 1923 earthquake. This means that the average slip deficit per recurrence interval is ???50% of the relative plate convergence. These findings of plate motion partitioning and slip deficit lead us to suggest that instead of a simple recurrence model with characteristic earthquakes, additional mechanisms are necessary to describe the accommodation of deformation in the Kanto region. So far, obvious candidates for these alternative mechanisms have not been discovered. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Remotely triggered microearthquakes and tremor in central California following the 2010 M w 8.8 Chile earthquake

We examine remotely triggered microearthquakes and tectonic tremor in central California following the 2010 M w 8.8 Chile earthquake. Several microearthquakes near the Coso Geothermal Field were apparently triggered, with the largest earthquake ( M l 3.5) occurring during the large-amplitude Love surface waves. The Chile mainshock also triggered numerous tremor bursts near the Parkfield-Cholame section of the San Andreas Fault (SAF). The locally triggered tremor bursts are partially masked at lower frequencies by the regionally triggered earthquake signals from Coso, but can be identified by applying high-pass or matched filters. Both triggered tremor along the SAF and the M l 3.5 earthquake in Coso are consistent with frictional failure at different depths on critically-stressed faults under the Coulomb failure criteria. The triggered tremor, however, appears to be more phase-correlated with the surface waves than the triggered earthquakes, likely reflecting differences in constitutive properties between the brittle, seismogenic crust and the underlying lower crust.

California↗

Delayed dynamic triggering of disposal-induced earthquakes observed by a dense array in Northern Oklahoma

Recent increases in earthquake occurrence rates in Oklahoma have been linked to the injection of large volumes of saltwater, a byproduct of oil and gas extraction. Here we present a detailed study of remote earthquake triggering in an area of active injection‐induced seismicity in northern Oklahoma using data from the LArge‐n Seismic Survey in Oklahoma (LASSO) temporary array and nearby permanent broadband seismic stations. We estimate changes in earthquake rates and calculate the Coulomb failure stress changes on potential receiver faults due to passing teleseismic surface waves. A statistically significant increase in seismicity is observed ∼8 hr after the 16 April 2016 M w 7.8 Ecuador earthquake. The Coulomb stress changes associated with the Ecuador earthquake are on the order of ∼1 kPa. Physical mechanisms consistent with the observed dynamic stress threshold include failure driven by activation of aseismic slip or hydrological response of the fault system.

Oklahoma↗

The 2020 Westmorland, California earthquake swarm as aftershocks of a slow slip event sustained by fluid flow

Swarms are bursts of earthquakes without an obvious mainshock. Some have been observed to be associated with transient aseismic fault slip, while others are thought to be related to fluids. However, the association is rarely quantitative due to insufficient data quality. We use high-quality GPS/GNSS, InSAR, and relocated seismicity to study a swarm of >2,000 earthquakes which occurred between 30 September and 6 October 2020, near Westmorland, California. Using 5 min sampled Global Positioning System (GPS) supplemented with InSAR, we document a spontaneous shallow M w 5.2 slow slip event that preceded the swarm by 2–15 hr. The earthquakes in the early phase were predominantly non-interacting and driven primarily by the slow slip event resulting in a nonlinear expansion. A stress-driven model based on the rate-and-state friction successfully explains the overall spatial and temporal evolution of earthquakes, including the time lag between the onset of the slow slip event and the swarm. Later, a distinct back front and a square root of time expansion of clustered seismicity on en-echelon fault structures suggest that fluids helped sustain the swarm. Static stress triggering analysis using Coulomb stress and statistics of interevent times suggest that 45%–65% of seismicity was driven by the slow slip event, 10%–35% by inter-earthquake interactions, and 10%–30% by fluids. Our model also provides constraints on the friction parameter and the pore pressure and suggests that this swarm behaved like an aftershock sequence but with the mainshock replaced by the slow slip event.

California↗

Conference explores mechanical involvement of fluids in faulting

A growing body of evidence suggests that fluids are intimately linked to a variety of faulting processes. These include the long-term structural and compositional evolution of fault zones; fault creep; and the nucleation, propagation, arrest, and recurrence of earthquake ruptures. Besides the widely recognized physical role of fluid pressures in controlling the strength of crustal fault zones, it is also apparent that fluids can exert mechanical influence through a variety of chemical effects. To address these issues, a “Red-Book” Conference on the Mechanical Effects of Fluids in Faulting was sponsored by the U.S. Geological Survey under the auspices of the National Earthquake Hazards Reduction Program at Fish Camp, Calif., from June 6–10, 1993. The coconvenors were Steve Hickman, Rick Sibson, and Ron Bruhn.

California↗

The October 17, 1989, Loma Prieta, California, Earthquake and its aftershocks: Geometry of the sequence from high-resolution locations

Hypocenters of the Loma Prieta sequence form a dipping zone that rises from the mainshock hypocenter and is parallel to the mainshock nodal plane. Most aftershocks cluster around the perimeter of the zone, surrounding a relatively aseismic center which approximates the region of mainshock rupture. At its southeastern end, the dipping aftershock zone warps into a vertical surface that corresponds to the San Andreas fault. In the central and northwestern parts of the zone at depths above ∼10 km, the aftershocks define numerous disjoint fault structures. The large component of reverse-slip observed in this event agrees with a simple model for slip on a dipping plane within a compressional fault bend. We do not believe that the Loma Prieta earthquake occurred on the Sargent fault. However, we are unable to conclude whether it ruptured the principal plate boundary fault or a less frequently active fault.

California↗

Earthquakes of the Holocene

Geologic studies of earthquakes involve mapping of coseismic surface faulting and secondary deformation from historical events, trenching and geomorphic analysis to define the timing and size of past earthquakes, and investigations of fault zone structure and geometry in both unconsolidated sediments and bedrock. This research is now being referred to as paleoseismology, seismic geology, and earthquake geology. Since the mid-1970s, it has led to some of the most exciting and important contributions to the understanding of earthquake behavior in space and time [Hanks, 1985; Allen, 1986]. The present report is the first of what is hoped will become a regularly contributed summary of the geologic aspects of the study of earthquakes.

Reviews of Geophysics↗

Analysing the 1811-1812 New Madrid earthquakes with recent instrumentally recorded aftershocks

Although dynamic stress changes associated with the passage of seismic waves are thought to trigger earthquakes at great distances, more than 60 per cent of all aftershocks appear to be triggered by static stress changes within two rupture lengths of a mainshock. The observed distribution of aftershocks may thus be used to infer details of mainshock rupture geometry. Aftershocks following large mid-continental earthquakes, where background stressing rates are low, are known to persist for centuries, and models based on rate-and-state friction laws provide theoretical support for this inference. Most past studies of the New Madrid earthquake sequence have indeed assumed ongoing microseismicity to be a continuing aftershock sequence. Here we use instrumentally recorded aftershock locations and models of elastic stress change to develop a kinematically consistent rupture scenario for three of the four largest earthquakes of the 1811-1812 New Madrid sequence. Our results suggest that these three events occurred on two contiguous faults, producing lobes of increased stress near fault intersections and end points, in areas where present-day microearthquakes have been hitherto interpreted as evidence of primary mainshock rupture. We infer that the remaining New Madrid mainshock may have occurred more than 200 km north of this region in the Wabash Valley of southern Indiana and Illinois-an area that contains abundant modern microseismicity, and where substantial liquefaction was documented by historic accounts. Our results suggest that future large midplate earthquake sequences may extend over a much broader region than previously suspected.

Nature↗

Effects of acoustic waves on stick-slip in granular media and implications for earthquakes

It remains unknown how the small strains induced by seismic waves can trigger earthquakes at large distances, in some cases thousands of kilometres from the triggering earthquake, with failure often occurring long after the waves have passed. Earthquake nucleation is usually observed to take place at depths of 10-20 km, and so static overburden should be large enough to inhibit triggering by seismic-wave stress perturbations. To understand the physics of dynamic triggering better, as well as the influence of dynamic stressing on earthquake recurrence, we have conducted laboratory studies of stick-slip in granular media with and without applied acoustic vibration. Glass beads were used to simulate granular fault zone material, sheared under constant normal stress, and subject to transient or continuous perturbation by acoustic waves. Here we show that small-magnitude failure events, corresponding to triggered aftershocks, occur when applied sound-wave amplitudes exceed several microstrain. These events are frequently delayed or occur as part of a cascade of small events. Vibrations also cause large slip events to be disrupted in time relative to those without wave perturbation. The effects are observed for many large-event cycles after vibrations cease, indicating a strain memory in the granular material. Dynamic stressing of tectonic faults may play a similar role in determining the complexity of earthquake recurrence. ??2007 Nature Publishing Group.

Nature↗

Strike-slip 23 January 2018 MW 7.9 Gulf of Alaska rare intraplate earthquake: Complex rupture of a fracture zone system

Large intraplate earthquakes in oceanic lithosphere are rare and usually related to regions of diffuse deformation within the oceanic plate. The 23 January 2018 M W 7.9 strike-slip Gulf of Alaska earthquake ruptured an oceanic fracture zone system offshore Kodiak Island. Bathymetric compilations show a muted topographic expression of the fracture zone due to the thick sediment that covers oceanic basement but the fracture zone system can be identified by offset N-S magnetic anomalies and E-W linear zones in the vertical gravity gradient. Back-projection from global seismic stations reveals that the initial rupture at first propagated from the epicenter to the north, likely rupturing along a weak zone parallel to the ocean crustal fabric. The rupture then changed direction to eastward directed with most energy emitted on Aka fracture zone resulting in an unusual multi-fault earthquake. Similarly, the aftershocks show complex behavior and are related to two different tectonic structures: (1) events along N-S trending oceanic fabric, which ruptured mainly strike-slip and additionally, in normal and oblique slip mechanisms and (2) strike-slip events along E-W oriented fracture zones. To explain the complex faulting behavior we adopt the classical stress and strain partitioning concept and propose a generalized model for large intra-oceanic strike-slip earthquakes of trench-oblique oriented fracture zones/ocean plate fabric near subduction zones. Taking the Kodiak asperity position of 1964 maximum afterslip and outer-rise Coulomb stress distribution into account, we propose that the unusual 2018 Gulf of Alaska moment release was stress transferred to the incoming oceanic plate from co- and post-processes of the nearby great 1964 M W 9.2 megathrust earthquake.

Alaska↗

Three-dimensional seismic structure and moment tensors of non-double-couple earthquakes at the Hengill-Grensdalur volcanic complex, Iceland

The volcanic and geothermal areas of Iceland are rich sources of non-double-couple (non-DC) earthquakes. A state-of-the-art digital seismometer network deployed at the Hengill–Grensdalur volcanic complex in 1991 recorded 4000 small earthquakes. We used the best recorded of these to determine 3-D V P and V P /V S structure tomographically and accurate earthquake moment tensors. The V P field is dominated by high seismic wave speed bodies interpreted as solidified intrusions. A widespread negative (−4 per cent) V P /V S anomaly in the upper 4km correlates with the geothermal field, but is too strong to be caused solely by the effect of temperature upon liquid water or the presence of vapour, and requires in addition mineralogical or lithological differences between the geothermal reservoir and its surroundings. These may be caused by geothermal alteration. Well-constrained moment tensors were obtained for 70 of the best-recorded events by applying linear programming methods to P - and S -wave polarities and amplitude ratios. About 25 per cent of the mechanisms are, within observational error, consistent with DC mechanisms consistent with shear faulting. The other 75 per cent have significantly non-DC mechanisms. Many have substantial explosive components, one has a substantial implosive component, and the deviatoric component of many is strongly non-DC. Many of the non-DC mechanisms are consistent, within observational error, with simultaneous tensile and shear faulting. However, the mechanisms occupy a continuum in source-type parameter space and probably at least one additional source process is occurring. This may be fluid flow into newly formed cracks, causing partial compensation of the volumetric component. Studying non-shear earthquakes such as these has great potential for improving our understanding of geothermal processes and earthquake source processes in general.

Geophysical Journal International↗

From coseismic offsets to fault-block mountains

In the Basin and Range extensional province of the western United States, coseismic offsets, under the influence of gravity, display predominantly subsidence of the basin side (fault hanging wall), with comparatively little or no uplift of the mountainside (fault footwall). A few decades later, geodetic measurements [GPS and interferometric synthetic aperture radar (InSAR)] show broad (∼100 km) aseismic uplift symmetrically spanning the fault zone. Finally, after millions of years and hundreds of fault offsets, the mountain blocks display large uplift and tilting over a breadth of only about 10 km. These sparse but robust observations pose a problem in that the coesismic uplifts of the footwall are small and inadequate to raise the mountain blocks. To address this paradox we develop finite-element models subjected to extensional and gravitational forces to study time-varying deformation associated with normal faulting. Stretching the model under gravity demonstrates that asymmetric slip via collapse of the hanging wall is a natural consequence of coseismic deformation. Focused flow in the upper mantle imposed by deformation of the lower crust localizes uplift, which is predicted to take place within one to two decades after each large earthquake. Thus, the best-preserved topographic signature of earthquakes is expected to occur early in the postseismic period.

Proceedings of the National Academy of Sciences of↗

Stress- and structure-controlled anisotropy in a region of complex faulting—Yuha Desert, California

We examine shear velocity anisotropy in the Yuha Desert, California using aftershocks of the 2010 M7.2 El Mayor-Cucapah earthquake. The Yuha Desert is underlain by a complex network of right- and left-lateral conjugate faults, some of which experienced triggered slip during the El Mayor-Cucapah earthquake. An automated method that implements multiple measurement windows and a range of bandpass filters is used to estimate the fast direction ( ϕ ) and delay time ( δt ) of the split shear waves. We find an average ϕ oriented approximately north–south suggesting it is primarily controlled by the regional maximum compressive stress direction. However, the spatial variability in ϕ reveals that the fault structures that underlie the Yuha Desert also influence the measured splitting parameters. We infer that the northeast- and northwest-oriented ϕ reflect shear fabric subparallel to the conjugate fault structures. We do not observe a simple correlation between δt and hypocentral distance. Instead, the observed spatial variation in δt suggests that near-source variation in anisotropic strength may be equal to or more important than effects local to the station. No temporal variation in splitting parameters is observed during the 70-day period following the main shock. In this region of complex faulting, we observe a spatially variable pattern of anisotropy that is both stress- and structure-controlled. This study suggests that shear fabric can form even along short, discontinuous fault strands with minimal offset.

California↗

Fault locking, block rotation and crustal deformation in the Pacific Northwest

We interpret Global Positioning System (GPS) measurements in the northwestern United States and adjacent parts of western Canada to describe relative motions of crustal blocks, locking on faults and permanent deformation associated with convergence between the Juan de Fuca and North American plates. To estimate angular velocities of the oceanic Juan de Fuca and Explorer plates and several continental crustal blocks, we invert the GPS velocities together with seafloor spreading rates, earthquake slip vector azimuths and fault slip azimuths and rates. We also determine the degree to which faults are either creeping aseismically or, alternatively, locked on the block-bounding faults. The Cascadia subduction thrust is locked mainly offshore, except in central Oregon, where locking extends inland. Most of Oregon and southwest Washington rotate clockwise relative to North America at rates of 0.4–1.0 ° Myr −1 . No shear or extension along the Cascades volcanic arc has occurred at the mm/yr level during the past decade, suggesting that the shear deformation extending northward from the Walker Lane and eastern California shear zone south of Oregon is largely accommodated by block rotation in Oregon. The general agreement of vertical axis rotation rates derived from GPS velocities with those estimated from palaeomagnetic declination anomalies suggests that the rotations have been relatively steady for 10–15 Ma. Additional permanent dextral shear is indicated within the Oregon Coast Range near the coast. Block rotations in the Pacific Northwest do not result in net westward flux of crustal material—the crust is simply spinning and not escaping. On Vancouver Island, where the convergence obliquity is less than in Oregon and Washington, the contractional strain at the coast is more aligned with Juan de Fuca—North America motion. GPS velocities are fit significantly better when Vancouver Island and the southern Coast Mountains move relative to North America in a block-like fashion. The relative motions of the Oregon, western Washington and Vancouver Island crustal blocks indicate that the rate of permanent shortening, the type that causes upper plate earthquakes, across the Puget Sound region is 4.4 ± 0.3 mm yr −1 . This shortening is likely distributed over several faults but GPS data alone cannot determine the partitioning of slip on them. The transition from predominantly shear deformation within the continent south of the Mendocino Triple Junction to predominantly block rotations north of it is similar to changes in tectonic style at other transitions from shear to subduction. This similarity suggests that crustal block rotations are enhanced in the vicinity of subduction zones possibly due to lower resisting stress.

Geophysical Journal International↗

Long-term dynamics of earthquake swarms in the Yellowstone caldera

The factors controlling the spatial distribution and temporal evolution of earthquake swarms in volcanic systems remain unclear. We leverage leading-edge deep learning algorithms and a detailed three-dimensional velocity model to construct a 15-year high-resolution earthquake catalog of the Yellowstone caldera region. More than half of the region’s earthquakes are clustered into swarm-like families characterized by episodes of hypocenter expansion and migration. Adjacent earthquake swarms, separated by long quiescent periods, are found to be a dominant feature. We suggest that these swarms are controlled by the interplay between slowly diffusing aqueous fluids and rapid episodic fluid injections, which may result from the breaking of permeability seals. Our analyses also indicate that clustered seismicity beneath the caldera occurs on relatively immature, rougher fault structures, compared to more planar faults outside. Our results provide additional context for understanding seismicity in hydrothermal systems, highlighting the key role played by long-term fluid diffusion processes in driving the occurrence of earthquake swarms.

Idaho, Montana, Wyoming↗

Block rotation in east-central Alaska: A framework for evaluating earthquake potential?

Geologic and seismic data reveal a set of parallel, active, strike-slip faults in east-central Alaska between the Denali and Tintina fault systems. The faults strike northeast to north-northeast, at a high angle to the bounding dextral fault systems, and exhibit sinistral slip. We hypothesize that this set of faults divides the crust into long blocks that are rotating clockwise in response to northerly compression resulting from Pacific–North American plate convergence. We suggest that these faults have produced most of the large historical earthquakes in east-central Alaska between the Alaska Range and the Yukon River.

Geology↗

Cascadia subduction tremor muted by crustal faults

Deep, episodic slow slip on the Cascadia subduction megathrust of western North America is accompanied by low-frequency tremor in a zone of high fluid pressure between 30 and 40 km depth. Tremor density (tremor epicenters per square kilometer) varies along strike, and lower tremor density statistically correlates with upper plate faults that accommodate northward motion and rotation of forearc blocks. Upper plate earthquakes occur to 35 km depth beneath the faults. We suggest that the faults extend to the overpressured megathrust, where they provide fracture pathways for fluid escape into the upper plate. This locally reduces megathrust fluid pressure and tremor occurrence beneath the faults. Damping of tremor and related slow slip caused by fluid escape could affect fault properties of the megathrust, possibly influencing the behavior of great earthquakes.

Geology↗