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

USGS2009

Co-seismic ruptures of the 12 May 2008, M s 8.0 Wenchuan earthquake, Sichuan: East-west crustal shortening on oblique, parallel thrusts along the eastern edge of Tibet

The M s 8.0, Wenchuan earthquake, which devastated the mountainous western rim of the Sichuan basin in central China, produced a surface rupture over 200??km-long with oblique thrust/dextral slip and maximum scarp heights of ~ 10??m. It thus ranks as one of the world's largest continental mega-thrust events in the last 150??yrs. Field investigation shows clear surface breaks along two of the main branches of the NE-trending Longmen Shan thrust fault system. The principal rupture, on the NW-dipping Beichuan fault, displays nearly equal amounts of thrust and right-lateral slip. Basin-ward of this rupture, another continuous surface break is observed for over 70??km on the parallel, more shallowly NW-dipping Pengguan fault. Slip on this latter fault was pure thrusting, with a maximum scarp height of ~ 3.5??m. This is one of the very few reported instances of crustal-scale co-seismic slip partitioning on parallel thrusts. This out-of-sequence event, with distributed surface breaks on crustal mega-thrusts, highlights regional, ~ EW-directed, present day crustal shortening oblique to the Longmen Shan margin of Tibet. The long rupture and large offsets with strong horizontal shortening that characterize the Wenchuan earthquake herald a re-evaluation of tectonic models anticipating little or no active shortening of the upper crust along this edge of the plateau, and require a re-assessment of seismic hazard along potentially under-rated active faults across the densely populated western Sichuan basin and mountains. ?? 2009 Elsevier B.V.

Earth and Planetary Science Letters↗

USGS2025

Dendroseismological investigation of redwood trees along the North Coast section of the San Andreas Fault

Sequoia sempervirens (coast redwood) tree rings have the potential to annually resolve late-Holocene earthquakes on the northern San Andreas Fault based on direct (e.g., physical damage) and indirect (e.g., co-seismic environmental change) impacts, but scarcity of suitable samples and challenges crossdating this long-lived species have limited progress. More precise dating of the pre-1906 (penultimate) earthquake can improve hazard assessment and understanding of rupture segmentation. We target old trees (maximum >815 yr) along the North Coast section of the fault (increment cores via rope-climbing, 11 living trees; plunge cuts, 23 stumps) and employ complementary disturbance detection methods including radial-growth averaging (tree- and series-level), cataloging anatomical indicators (e.g., traumatic resin ducts, TRD), and dating structural components (e.g., reiterated trunks, leans). Multi-centennial ring-width chronologies at Fort Ross (1569−2023) and Gualala (1397−2023) promote continued study with incomplete crossdating limiting utilization of some series. Growth pulses (reductions, releases) and TRD dispersed across the record reflect dynamic environments that obfuscate detection of earthquake signals. The 1906 earthquake did not leave strong signatures on most trees, and when it did, within-tree response varied from normal presentation to discoloration, TRD, and missing rings. Synchrony of indicators at both locations identified 1678−1680 (6 of 15 trees) and 1698−1700 (8 of 16 trees) as the strongest disturbances among dated rings in the time range of the penultimate earthquake, peaking at 1698 (15.7 % of possible growth and anatomical indicators), but the triggering mechanisms for these events are unknown.

California↗

USGS2003

Paleoseismology of the Chelungpu Fault during the past 1900 years

The 1999 earthquake brought about 80-km-long surface ruptures along the Shihkang, Chelungpu, and Tajienshan Faults, central Taiwan. Several trenches have been excavated across the Chelungpu Fault of the middle segment. The surface ruptures display clear scarps ranging from 0.2 to 4 m high, showing a complex geomorphic pattern due to coseismic faulting and folding. In the study, measurement of the vertical offset or structural relief was taken with reference to the hanging wall beyond the trishear deformation zone. Therefore we suggest that, for the measurement of offset, we should disregard the trishear zone, and that structural relief on the hanging wall should be represented as a real vertical offset. The net slip is then calculated from the structural relief and dip angle of the thrust on a vertical plane along the slip direction. Through the excavation of a pineapple field across the Chelungpu Fault, we are able to provide evidence of at least four earthquake events for the past about 1900 years, including the 1999 earthquake. Furthermore, based on the radiocarbon dates and historical record, the timing of the penultimate event is bracketed to be between 430 and 150 years ago, and the average recurrence interval is less than 700 years. These data indicate that the average slip rate is about 8.7 mm/yr for the past 1900 years. ?? 2003 Published by Elsevier Ltd.

Quaternary International↗

USGS1999

Stress transfer by the 1988-1989 M=5.3 and 5.4 Lake Elsman foreshocks to the Loma Prieta fault: Unclamping at the site of peak mainshock slip

We study the stress transferred by the June 27, 1988, M = 5.3 and August 8, 1989, M = 5.4 Lake Elsman earthquakes, the largest events to strike within 15 km of the future Loma Prieta rupture zone during 74 years before the 1989 M = 6.9 Loma Prieta earthquake. We find that the first Lake Elsman event brought the rupture plane of the second event 0.3–1.6 bars (0.03–0.16 MPa) closer to Coulomb failure but that the Lake Elsman events did not bring the future Loma Prieta hypocentral zone closer to failure. Instead, the Lake Elsman earthquakes are calculated to have reduced the normal stress on (or “undamped”) the Loma Prieta rupture surface by 0.5–1.0 bar (0.05–0.10 MPa) at the site where the greatest slip subsequently occurred in the Loma Prieta earthquake. This association between the sites of peak unclamping and slip suggests that the Lake Elsman events did indeed influence the Loma Prieta rupture process. Unclamping the fault would have locally lowered the resistance to sliding. Such an effect could have been enhanced if the lowered normal stress permitted fluid infusion into the undamped part of the fault. Although less well recorded, the M L = 5.0 1964 and M L = 5.3 1967 Corralitos events struck within 10 km of the southwest end of the future Loma Prieta rupture. No similar relationship between the normal stress change and subsequent Loma Prieta slip is observed, although the high-slip patch southwest of the Loma Prieta epicenter corresponds roughly to the site of calculated Coulomb stress increase for a low coefficient of friction. The Lake Elsman-Loma Prieta result is similar to that for the 1987 M = 6.2 Elmore Ranch and M = 6.7 Superstition Hills earthquakes, suggesting that foreshocks might influence the distribution of mainshock slip rather than the site of mainshock nucleation.

Journal of Geophysical Research B: Solid Earth↗

USGS2003

Source model for the Mw 6.7, 23 October 2002, Nenana Mountain earthquake (Alaska) from InSAR

The 23 October 2002 Nenana Mountain Earthquake (M w ∼ 6.7) occurred on the Denali Fault (Alaska), to the west of the M w ∼ 7.9 Denali Earthquake that ruptured the same fault 11 days later. We used 6 interferograms, constructed using radar images from the Canadian Radarsat-1 and European ERS-2 satellites, to determine the coseismic surface deformation and a source model. Data were acquired on ascending and descending satellite passes, with incidence angles between 23 and 45 degrees, and time intervals of 72 days or less. Modeling the event as dislocations in an elastic half space suggests that there was nearly 0.9 m of right-lateral strike-slip motion at depth, on a near-vertical fault, and that the maximum slip in the top 4 km of crust was less than 0.2 m. The Nenana Mountain Earthquake increased the Coulomb stress at the future hypocenter of the 3 November 2002, Denali Earthquake by 30–60 kPa.

Alaska↗

USGS2003

High precision earthquake locations reveal seismogenic structure beneath Mammoth Mountain, California

In 1989, an unusual earthquake swarm occurred beneath Mammoth Mountain that was probably associated with magmatic intrusion. To improve our understanding of this swarm, we relocated Mammoth Mountain earthquakes using a double difference algorithm. Relocated hypocenters reveal that most earthquakes occurred on two structures, a near-vertical plane at 7–9 km depth that has been interpreted as an intruding dike, and a circular ring-like structure at ∼5.5 km depth, above the northern end of the inferred dike. Earthquakes on this newly discovered ring structure form a conical section that dips outward away from the aseismic interior. Fault-plane solutions indicate that in 1989 the seismicity ring was slipping as a ring-normal fault as the center of the mountain rose with respect to the surrounding crust. Seismicity migrated around the ring, away from the underlying dike at a rate of ∼0.4 km/month, suggesting that fluid movement triggered seismicity on the ring fault.

California↗

USGS2009

Rupture parameters of the 2003 Zemmouri (M w 6.8), Algeria, earthquake from joint inversion of interferometric synthetic aperture radar, coastal uplift, and GPS

We study the surface deformation associated with the 21 May 2003 (M w = 6.8) Zemmouri (Algeria) earthquake, the strongest seismic event felt in the Algiers region since 1716. The thrust earthquake mechanism and related surface deformation revealed an average 0.50 m coastal uplift along ??55-km-long coastline. We obtain coseismic interferograms using Envisat advanced synthetic aperture radar (ASAR) (IS2) and RADARSAT standard beam (ST4) data from both the ascending and descending orbits of Envisat satellite, whereas the RADARSAT data proved useful only in the descending mode. While the two RADARSAT interferograms cover the earthquake area, Envisat data cover only the western half of the rupture zone. Although the interferometric synthetic aperture radar (InSAR) coherence in the epicenter area is poor, deformation fringes are observed along the coast in different patches. In the Boumerdes area, the maximum coseismic deformation is indicated by the high gradient of fringes visible in all interferograms in agreement with field measurements (tape, differential GPS, leveling, and GPS). To constrain the earthquake rupture parameters, we model the interferograms and uplift measurements using elastic dislocations on triangular fault patches in an elastic and homogeneous half-space. We invert the coseismic slip using first, a planar surface and second, a curved fault, both constructed from triangular elements using Poly3Dinv program that uses a damped least square minimization. The best fit of InSAR, coastal uplift, and GPS data corresponds to a 65-km-long fault rupture dipping 40?? to 50?? SE, located at 8 to 13 km offshore with a change in strike west of Boumerdes from N60??-65?? to N95??-105??. The inferred rupture geometry at depth correlates well with the seismological results and may have critical implications for the seismic hazard assessment of the Algiers region. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗

USGS2012

Coseismic and postseismic stress rotations due to great subduction zone earthquakes

The three largest recent great subduction zone earthquakes (2011 M9.0 Tohoku, Japan; 2010 M8.8 Maule, Chile; and 2004 M9.2 Sumatra-Andaman) exhibit similar coseismic rotations of the principal stress axes. Prior to each mainshock, the maximum compressive stress axis was shallowly plunging, while immediately after the mainshock, both the maximum and minimum compressive stress axes plunge at ~45°. Dipping faults can be oriented for either reverse or normal faulting in this post-mainshock stress field, depending on their dip, explaining the observed normal-faulting aftershocks without requiring a complete reversal of the stress field. The significant stress rotations imply near-complete stress drop in the mainshocks, with >80% of the pre-mainshock stress relieved in the Tohoku and Maule earthquakes and in the northern part of the Sumatra-Andaman rupture. The southern part of the Sumatra-Andaman rupture relieved ~60% of the pre-mainshock stress. The stress axes rotated back rapidly in the months following the Tohoku and Maule mainshocks, and similarly in the southern part of the Sumatra-Andaman rupture. A rapid postseismic rotation is possible because the near-complete stress drop leaves very little “background” stress at the beginning of the postseismic reloading. In contrast, there has been little or no postseismic rotation in the northern part of the Sumatra-Andaman rupture over the 7 years since the mainshock. All M ≥8.0 subduction earthquakes since 1990 with an adequate number of pre- and post-mainshock events were evaluated, and not all show similar coseismic stress rotations. Deeper earthquakes exhibit smaller coseismic stress rotations, likely due to increasing deviatoric stress with depth.

Geophysical Research Letters↗

USGS2019

Earthquake stress drop and Arias Intensity

We determine earthquake stress drops directly from the Arias intensity database of NGA-West2. Arias intensity (Arias, 1970) is an engineering measure proportional to the integral of the absolute value of acceleration squared, over the significant duration of the signal. As such, it is closely related to root-mean-square acceleration, and can readily be connected to earthquake stress drop (Hanks and McGuire, 1981). Arias intensity records out to 100 km yield stable stress drops for moderate-to-large magnitude earthquakes, M6.5+; for smaller events ~M4.5 – 6.5, only closer-in records yield stable results. For the 116 events considered, stress drops are about 35% larger for Class 1 mainshocks than for traditional on-fault Class 2 aftershocks, and smaller for those aftershocks close to the main fault plane. Aftershock stress drops show large variability, however, implying that on average they re-rupture weakened patches, but can also rupture intact rock or high-stress asperities. We observe an increase of stress drop with earthquake depth similar to that of other studies but do not find any significant faulting mechanism dependence. The variability of the Arias intensity-based stress drop is lower than that of eGf-based stress drops from Baltay et al. (2010, 2011), and nearly on par with variability seen in ground-motion prediction equations. The Arias intensity stress drop is a novel and promising method to estimate stress drop without the need for path and site corrections, and yields further insight into the connection between source physics and ground-motion.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

Comment on “222Rn Premonitory signals for earthquakes?”

Radon and other terrestrial gases in groundwater and soil air have been studied in recent years in seismically active areas, especially in China, Japan, the former Soviet Union, and the United States, in search of changes that may be useful for earthquake prediction. Concentrations of radon, helim, hydrogen, mercury vapor, carbon dioxide, and several other gases have been found to be high along active faults. Faults may thus be paths of least resistance along which terrestrial gases generated or stored in the Earth can escape into the atmosphere. Gas-concentration changes with durations of hours to months have been observed before and after earthquakes at relatively few favorably situated sites, and their patterns have not always been repeated. Some of these sites are located as far as several hundred kilometers away from earthquake epicenters. These “sensitive” sites are generally situated along active faults, whether the seismogenic ones or not, especially at their intersections or bends, or at other structurally weak zones, possibly because of concentration of tectonic strain and abundance of pore fluids at such places (see review papers by King , 1986, 1990). Due to this selectivity and for two other important reasons, studies of terrestrial gases for earthquake prediction have not been viewed favorably by most geophysicists, especially those in the United States.

Eos Science News↗

USGS2012

Multicycle slip distribution along a laboratory fault

Slip distribution along a laboratory fault, which consists of eight spring-connected blocks that are elastically driven to slide on a frictional surface, has been examined for a “long” sequence of slip events to test the applicability of some conceptual models proposed recently in the literature. The distributions of large slip events are found to be quite variable and do not fit the uniform slip or characteristic earthquake models. The rupture initiation points are usually not near the corresponding maximum slip points, in contrast to observations by Thatcher (1990) and by Fukao and Kikuchi (1987) that earthquake hypocenters are commonly near corresponding regions of maximum slip in the fault planes. This contrast may suggest that either the present observations or theirs are not representative or the teleseismically determined hypocenters may not always be true rupture initiation points as usually assumed. Large slip events are also found to be a stress-roughening process. They are triggered by some small events after the stresses have been adjusted by some earlier small-to-moderate events to be near the critical levels at most locations along the fault. This suggests that earthquake prediction monitoring efforts should not be limited to a small region near an asperity but should be spread out to cover the entire fault segment in a seismic gap in order to detect the condition of simultaneous strain buildup.

Journal of Geophysical Research↗

USGS2012

Formation of left-lateral fractures within the Summit Ridge Shear Zone, 1989 Loma Prieta, California, Earthquake

The 1989 Loma Prieta, California, earthquake is characterized by the lack of major, throughgoing, coseismic, right-lateral faulting along strands of the San Andreas fault zone in the epicentral area. Instead, throughout the Summit Ridge area there are zones of tension cracks and left-lateral fracture zones oriented about N45°W, that is, roughly parallel to the San Andreas fault in this area. The left-lateral fractures zones are enigmatic because their left-lateral slip is opposite to the right-lateral sense of the relative motion between the Pacific and North American plates. We suggest that the enigmatic fractures can be understood if we assume that coseismic deformation was by right-lateral shear across a broad zone, about 0.5 km wide and 4 km long, beneath Summit Ridge. The shear zone has an orientation of about N70° W and is oblique to the San Andreas fault zone, which has a trend of N45° to 50°W, so that right-lateral shearing, along with some dilation, can account for the orientation of the fracture zones. Using an analog, kinematic model of the Summit Ridge shear zone and the observable geometric parameters (a shear zone about 0.5 km wide, an acute angle of 25° between the tension cracks and the shear-zone walls, left-lateral slip of 5 cm, and spacing of about 100 m in the tension cracks) we estimate that the amount of right-lateral shift across the Summit Ridge shear zone was on the order of 1.4 m. This is comparable to the amount of slip for coseismic faulting at depth, 1.9 to 2.4 m, as computed by several investigators. Thus contrary to most previous reports on the Loma Prieta earthquake, which assert that coseismic, right-lateral ground rupture was restricted to considerable (greater than 4 km) depths in the epicentral area, we find that nearly all the right-lateral offset is represented at the ground surface by the Summit Ridge shear zone.

Journal of Geophysical Research B: Solid Earth↗

USGS1998

Rupture characteristics of the three M ∼ 4.7 (1992-1994) Parkfield earthquakes

Slip on the San Andreas fault was determined for three M &sim; 4.7 earthquakes using a tomographic inverse system [ Beroza and Spudich , 1988] to invert seismic source time functions (STFs) from S waves. STFs were obtained by deconvolving mainshock accelerograms by those from collocated smaller earthquakes. Accelerograms were from the U.S. Geological Survey Parkfield Small Aperture Array (UPSAR) and from a distributed array of digital accelerometer stations at Parkfield. Eight or nine STFs are used in each of the three inversions. STFs are typically symmetrical pulses with a duration of about 0.3–0.5 s. In the inversion, mainshock rise time was set to 0.05 s, and we allowed the rupture time to vary slightly from a constant rupture velocity of approximately 0.85β. Rupture for all three events, which are located in or close to the Middle Mountain preparation zone or box (MMB), quickly reaches a local maximum in slip and then propagates outward to peaks, ridges, or plateaus in the slip distribution. Slip for the October 20, 1992, event (located just inside the southern edge of the MMB) propagates from an initial spike north and updip along a curving ridge for about 2 km. The initial spike continued to grow in the November 14, 1993, event (located north of the October 20, 1992, event just beneath the hypocenter of the 1966 Parkfield earthquake), which shows little directivity, although there is a smaller patch of slip updip and to the south. In contrast, rupture for the December 20, 1994, event (located just south of the October 20, 1992, event) propagated north and slightly updip, creating a rough plateau in slip a few kilometers wide on a side. Directivity for this event also is to the north. Directivity for all three events points in the approximate direction of the 1966 hypocenter. Small pulses, which comprise a coda, are found on the STFs for several seconds after the initial impulsive event. Several tests based on the assumption that the average of all STFs from UPSAR for each event is an estimate of the true slip at the source suggest that the codas in the STFs are S waves from a long-duration source rather than uncorrected site response. An initiation phase is found on the array average for the November 14, 1993, and December 20, 1994, events. These precursory phases are the result of a spike in slip at the hypocenter. A value of 2.4–4 mm is obtained for D c , the slip-weakening distance, by interpreting the initial spike as a critical patch. The few aftershocks for the October 20, 1992, event are distributed to the north and updip of the mainshock, but the November 14, 1993, event had a strong burst of aftershock activity that propagated to the north of its hypocenter at roughly the same depth. Aftershocks of the December 20, 1994, event are mostly updip. The November 14, 1993, event had the simplest slip distribution, appeared to be the most impulsive, and had the most active aftershock sequence and the greatest depth. If the eventual Parkfield earthquake initiates near the 1966 hypocenter, then the directivity of the three events studied here will have pointed to it. However, it is certainly possible that both the initiation of characteristic Parkfield shocks and the directivity of smaller events are controlled by fault properties on a larger scale such as by fault bends or jogs.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Seismicity in the twenty years preceding the Loma Prieta California Earthquake

Persistent seismicity occurred at a low rate during at least the twenty years before the Loma Prieta earthquake along the 60 km-long rupture zone. The depth distribution of this seismicity forms a broad “U”-shape that delineates the previously locked rupture zone. Relocations of seismicity during the ten years before the earthquake relative to the Loma Prieta aftershocks show that this “U”-shaped distribution can be partitioned lengthwise into activity on two adjacent subparallel structures: a vertical fault beneath the San Andreas fault trace and an eastward-dipping blind fault beneath the Sargent fault trace. The 11–18 km deep, southwest-dipping part of the Loma Prieta rupture was not active during at least the preceding ten years. The slip geometry of intersecting fault structures in this zone could contribute to both the preparation process and the complexity of the rupture.

California↗

USGS2012

Seismic slip, segmentation, and the Loma Prieta Earthquake

We have plotted the cumulative seismic slip projected onto a vertical plane for earthquakes occurring during the last 20 years along 210 km of the San Andreas fault that includes the section that moved in the Loma Prieta earthquake. These plots illustrate the differences in depth and character of the seismicity between the locked and creeping portions of the fault or fault zone and define the segment upon which the Loma Prieta earthquake occurred. Working by analogy from the relation between pre-main shock microseismicity and presumed main shock slip regions at Parkfield and Loma Prieta, we identify a segment on the San Francisco Peninsula where we believe the M 7 1838 earthquake occurred, and which we believe may have accumulated sufficient strain that rupture should be expected in the coming decades.

California↗

USGS1992

Creep, compaction and the weak rheology of major faults

Field and laboratory observations suggest that the porosity within fault zones varies over earthquake cycles so that fluid pressure is in long-term equilibrium with hydrostatic fluid pressure in the country rock. Between earthquakes, ductile creep compacts the fault zone, increasing fluid pressure, and finally allowing frictional failure at relatively low shear stress. Earthquake faulting restores porosity and decreases fluid pressure to below hydrostatic. This mechanism may explain why major faults, such as the San Andreas system, are weak.

Nature↗

USGS2008

Stress changes from the 2008 Wenchuan earthquake and increased hazard in the Sichuan basin

On 12 May 2008, the devastating magnitude 7.9 (Wenchuan) earthquake struck the eastern edge of the Tibetan plateau, collapsing buildings and killing thousands in major cities aligned along the western Sichuan basin in China. After such a large-magnitude earthquake, rearrangement of stresses in the crust commonly leads to subsequent damaging earthquakes. The mainshock of the 12 May earthquake ruptured with as much as 9 m of slip along the boundary between the Longmen Shan and Sichuan basin, and demonstrated the complex strike-slip and thrust motion that characterizes the region. The Sichuan basin and surroundings are also crossed by other active strike-slip and thrust faults. Here we present calculations of the coseismic stress changes that resulted from the 12 May event using models of those faults, and show that many indicate significant stress increases. Rapid mapping of such stress changes can help to locate fault sections with relatively higher odds of producing large aftershocks. ??2008 Macmillan Publishers Limited. All rights reserved.

Nature↗

USGS2003

Mantle fault zone beneath Kilauea Volcano, Hawaii

Relocations and focal mechanism analyses of deep earthquakes (???13 kilometers) at Kilauea volcano demonstrate that seismicity is focused on an active fault zone at 30-kilometer depth, with seaward slip on a low-angle plane, and other smaller, distinct fault zones. The earthquakes we have analyzed predominantly reflect tectonic faulting in the brittle lithosphere rather than magma movement associated with volcanic activity. The tectonic earthquakes may be induced on preexisting faults by stresses of magmatic origin, although background stresses from volcano loading and lithospheric flexure may also contribute.

Science↗