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

USGS2025

Ultralong, supershear rupture of the 2025 Mw 7.7 Mandalay earthquake reveals unaccounted risk

The 28 March 2025 moment magnitude ( M w ) 7.7 earthquake in Mandalay, Burma (Myanmar), ruptured 475 kilometers of the Sagaing Fault, which was more than twice the length predicted by magnitude scaling relationships. Kinematic slip models and observation of a Rayleigh Mach wave that passed through parts of Thailand confirmed that rupture occurred at supershear velocities of greater than 5 kilometers per second. The anomalous length exposed a vast population to violent near-fault shaking. The Mandalay earthquake is a modern analog for the M w 7.9 1906 San Francisco earthquake, another atypically long and fast rupture. Probabilistic seismic hazard analyses use scaling relations that do not account for such long ruptures at moderate magnitudes. This limitation, in conjunction with a likely increased population and infrastructure exposure for atypically long ruptures, contributes to a potential mischaracterization of seismic risk.

Mandalay↗

USGS2007

Long-term perspectives on giant earthquakes and tsunamis at subduction zones

Histories of earthquakes and tsunamis, inferred from geological evidence, aid in anticipating future catastrophes. This natural warning system now influences building codes and tsunami planning in the United States, Canada, and Japan, particularly where geology demonstrates the past occurrence of earthquakes and tsunamis larger than those known from written and instrumental records. Under favorable circumstances, paleoseismology can thus provide long-term advisories of unusually large tsunamis. The extraordinary Indian Ocean tsunami of 2004 resulted from a fault rupture more than 1000 km in length that included and dwarfed fault patches that had broken historically during lesser shocks. Such variation in rupture mode, known from written history at a few subduction zones, is also characteristic of earthquake histories inferred from geology on the Pacific Rim. Copyright ?? 2007 by Annual Reviews. All rights reserved.

Annual Review of Earth and Planetary Sciences↗

USGS2011

A paleoseismic study along the central Denali Fault, Chistochina Glacier area, south-central Alaska

In the Chistochina Glacier area of south-central Alaska, the active trace of the Denali fault is well defined by prominent tectonic geomorphology, including scarps, grabens, and mole tracks associated with the 2002 Mw=7.9 Denali fault earthquake. Interpretation of a trench excavated across the 2002 rupture trace places a constraint on the timing of the penultimate earthquake to after 550 to 660 yr before 2002, consis- tent with other paleoseismic studies along the central Denali fault. Field measurement of offset moraine deposits and use of existing geochronologic data is the basis to estimate a minimum latest Pleistocene slip rate of 8–14 mm/yr.

Alaska↗

USGS2002

Primary surface rupture associated with the Mw 7.1 16 October 1999 Hector Mine earthquake, San Bernardino County, California

The Mw 7.1 Hector Mine earthquake occurred within the Mojave Desert portion of the eastern California shear zone and was accompanied by 48 km of dextral surface rupture. Complex northward rupture began on two branches of the Lavic Lake fault in the northern Bullion Mountains and also propagated southward onto the Bullion fault. Lesser amounts of rupture occurred across two right steps to the south. Surface rupture was mapped using postearthquake, 1:10,000-scale aerial photography. Field mapping provided additional detail and more than 400 fault-rupture observations; of these, approximately 300 measurements were used to characterize the slip distribution. En echelon surface rupture predominated in areas of thick alluvium, whereas in the bedrock areas, rupture was more continuous and focused within a narrower zone. Measured dextral offsets were relatively symmetrical about the epicentral region, with a maximum displacement of 5.25 ?? 0.85 m. Vertical slip was a secondary component and was variable, with minor west-side-down displacements predominat.ing in the Bullion Mountains. Field and aerial photographic evidence indicates that most of the faults that ruptured in 1999 had had prior late-Quaternary displacement, although only limited sections of the rupture show evidence for prior Holocene displacement.

Bulletin of the Seismological Society of America↗

USGS2009

Array deployment to observe rotational and translational ground motions along the Meishan Fault, Taiwan: A progress report

Measurements in the near field of earthquakes in Japan ( Takeo, 1998 ) and in Taiwan ( Huang et al. , 2006 ; Liu et al. , 2009 ) indicate that rotational ground motions are many times larger than expected from the classical elasticity theory. The Central Weather Bureau (CWB) deployed four instrumentation sets on the campus of the National Chung-Cheng University, Chia-yi, Taiwan, in order to study in detail earthquakes in the near field for both seismology and earthquake engineering. Both rotational and translational ground motions are being monitored along the active Meishan fault, where a major earthquake occurred in 1906, more than a century ago. The deployed instruments are: (1) a 32-element seismic array in the free field, (2) a 32-element seismic array in a building, (3) a six-channel unit with a broadband seismometer and an accelerometer, and (4) a six-channel unit with an accelerometer and an external rotational velocity sensor.

Bulletin of the Seismological Society of America↗

USGS2014

Modeling the effects of source and path heterogeneity on ground motions of great earthquakes on the Cascadia Subduction Zone Using 3D simulations

We ran finite‐difference earthquake simulations for great subduction zone earthquakes in Cascadia to model the effects of source and path heterogeneity for the purpose of improving strong‐motion predictions. We developed a rupture model for large subduction zone earthquakes based on a k −2 slip spectrum and scale‐dependent rise times by representing the slip distribution as the sum of normal modes of a vibrating membrane. Finite source and path effects were important in determining the distribution of strong motions through the locations of the hypocenter, subevents, and crustal structures like sedimentary basins. Some regions in Cascadia appear to be at greater risk than others during an event due to the geometry of the Cascadia fault zone relative to the coast and populated regions. The southern Oregon coast appears to have increased risk because it is closer to the locked zone of the Cascadia fault than other coastal areas and is also in the path of directivity amplification from any rupture propagating north to south in that part of the subduction zone, and the basins in the Puget Sound area are efficiently amplified by both north and south propagating ruptures off the coast of western Washington. We find that the median spectral accelerations at 5 s period from the simulations are similar to that of the Zhao et al. (2006) ground‐motion prediction equation, although our simulations predict higher amplitudes near the region of greatest slip and in the sedimentary basins, such as the Seattle basin.

Bulletin of the Seismological Society of America↗

USGS2016

A note on adding viscoelasticity to earthquake simulators

Here, I describe how time‐dependent quasi‐static stress transfer can be implemented in an earthquake simulator code that is used to generate long synthetic seismicity catalogs. Most existing seismicity simulators use precomputed static stress interaction coefficients to rapidly implement static stress transfer in fault networks with typically tens of thousands of fault patches. The extension to quasi‐static deformation, which accounts for viscoelasticity of Earth’s ductile lower crust and mantle, involves the precomputation of additional interaction coefficients that represent time‐dependent stress transfer among the model fault patches, combined with defining and evolving additional state variables that track this stress transfer. The new approach is illustrated with application to a California‐wide synthetic fault network.

Bulletin of the Seismological Society of America↗

USGS2024

Prediction of regional broadband strong ground motions using a teleseismic source model of the 18 April 2014 Mw 7.3 Papanoa, Mexico, earthquake

To estimate predicted ground motion from a teleseismic slip model, we use a low‐ and high‐frequency hybrid method to simulate the regional, strong ground motions observed following the 18 April 2014 moment magnitude ( ⁠M w ⁠ ) 7.3 Papanoa, Mexico, earthquake. To generate the regional ground motion at low frequencies (<1 Hz), a teleseismically derived, finite‐fault, kinematic model is used to define the earthquake source, taking into account slip‐model variations identified with a parameter sampling approach that considers possible errors in the fault geometry, the hypocenter depth, and the rupture velocity. A 3D crustal model is used to calculate the low‐frequency ground motions using a finite‐element calculation that includes topography and considers variations in the source model to estimate the uncertainty in the calculations. High frequencies (>1 Hz) are added using a 1D full‐wave propagation code that estimates uncertainties by considering multiple random distributions of slip with different spatial correlation lengths. The synthetic, broadband (0.05–10.0 Hz) ground motions are obtained by combining the low‐ and high‐frequency portions match filtered at 1 Hz. These synthetic ground motions are compared with the regional observations using velocity records, peak ground acceleration, and medians of the orientation‐independent response spectra of the horizontal components (RotD50) calculated at periods of 0.2, 0.3, 0.5, 1.0, 2.0, 3.0, 5.0, 7.5, and 10.0 s. The results indicate that ground motions estimated at these periods using our hybrid approach based primarily on a teleseismically derived source model are comparable to the values observed for the 2014 Papanoa earthquake at regional distances. The approach could be used to estimate strong‐motion spectral levels expected for regions with limited local and regional recordings and could also fill in magnitude or distance gaps in ground‐motion prediction relations utilized in the assessment of seismic hazard.

Papanoa↗

USGS2024

Static and dynamic strain in the 1886 Charleston, South Carolina, earthquake

During the 1886 Mw 7.3 Charleston, South Carolina, earthquake, three railroads emanating from the city were exposed to severe shaking. Expansion joints in segmented railroad tracks are designed to allow railroad infrastructure to withstand a few parts in 10,000 of thermoelastic strain. We show that, in 1886, transient contractions exceeding this limiting value buckled rails, and transient extensions pulled rails apart. Calculated values for dynamic strain in the meizoseismal region are in reasonable agreement with those anticipated from the relation between strain and moment magnitude proposed by Barbour et al. (2021) and exceed estimated tectonic strain released by the earthquake by an order of magnitude. Almost all of the documented disturbances of railroad lines, including evidence for shortening of the rails, can thus be ascribed to the effects of dynamic strain changes, not static strain. Little or no damage to railroads was reported outside the estimated 10 −4 dynamic strain contour. The correspondence between 10 −3 and 2×10 −4 contours of dynamic strain and Mercalli intensity 9 and 8, anticipated from the dependence of each quantity on peak ground velocity, suggests it may be possible to use railroad damage to quantitatively estimate shaking intensity. At one location, near Rantowles, ≈20 km west of Charleston, a photograph of buckled track taken one day after the earthquake has been cited as evidence for shallow dextral slip and has long focused a search for a causal fault in this region. Photogrammetric analysis reveals that the buckle was caused by transient contraction of <10 cm with no dextral offset. Our results further weaken the evidence for faulting in the swamps and forests south of the Ashley River in 1886, hitherto motivated by the photograph and limited macroseismic evidence for high‐intensity shaking.

South Carolina↗

USGS2026

The U.S. Geological Survey 2025 Puerto Rico and U.S. Virgin Islands time-independent earthquake rupture forecast

We present the 2025 U.S. Geological Survey Puerto Rico and U.S. Virgin Islands (PRVI) time‐independent earthquake rupture forecast (ERF), developed for the 2025 update to the National Seismic Hazard Model (NSHM) for PRVI. The updated ERF improves upon a prior model from 2003, including an expanded fault inventory with slip‐rate estimates, updated seismicity catalogs, and refined subduction zone geometries and deformation models. It applies the fault‐system inversion methodology to solve for rates of ruptures on modeled faults, adapted from the 2023 NSHM (NSHM23) for the western United States, including the first application of the inversion to model rates on a U.S. subduction interface. Off‐fault and intraslab seismicity are constrained by observed seismicity and use updated methods developed for NSHM23. Uncertainties in model components are substantial, and the ERF represents epistemic uncertainties through a comprehensive logic tree consisting of 1.7 billion logic‐tree branches combined across all sources.

Puerto Rico, U.S. Virgin Islands↗

USGS2026

Constraining source and path effects of large magnitude earthquakes using ground motion simulations

The purpose of this study is to use ground‐motion simulations to investigate ways in which source and path effects for large‐magnitude earthquakes can be represented in nonergodic ground‐motion models (GMMs). To achieve this, we designed a ground‐motion study in the San Francisco Bay Area that includes earthquakes with a broad range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of kinematic simulations (magnitude 4–7), we then develop a nonergodic GMM with the simulated data. We find that trends in the within‐site residuals are affected significantly by the earthquake radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within‐site residuals among multiple source realizations. Finally, after removing the source effects from the within‐site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site as the large events, whereas the second approach considers all small events on the fault plane. The results indicate that it is difficult to satisfactorily approximate the path effects of large events with those of small events using either approach, at least in the case of simulations.

California↗

USGS2020

Airborne lidar and electro-optical imagery along surface ruptures of the 2019 Ridgecrest earthquake sequence, Southern California

Surface rupture from the 2019 Ridgecrest earthquake sequence, initially associated with the M 6.4 foreshock, occurred on July 4 on a ~17 km long, northeast-southwest oriented, left-lateral zone of faulting. Following the M 7.1 mainshock on July 5 (local time), extensive northwest-southeast-oriented, right-lateral faulting was then also mapped along a ~50 km long zone of faults, including sub-parallel splays in several areas. The largest slip was observed in the epicentral area, and crossing the dry lakebed of China Lake to the southeast. Surface fault rupture mapping by a large team, reported elsewhere, was used to guide the airborne data acquisition reported here. Rapid rupture mapping allowed for accurate and efficient flight line planning for the high-resolution lidar and aerial photography. Flight line planning trade-offs were considered to allocate the medium (25 pulses per square meter, or ppsm) and high resolution (80 ppsm) lidar data collection polygons. The National Center for Airborne Laser Mapping (NCALM) acquired the airborne imagery with a Titan multispectral lidar system and DiMAC aerial digital camera, and USGS acquired GPS ground control data. This effort required extensive coordination with the Navy as much of the airborne data acquisition occurred within their restricted airspace at the China Lake Ranges.

Calilfornia↗

USGS2025

Deterministic physics-based earthquake sequence simulators match empirical ground-motion models and enable extrapolation to data poor regimes: Application to multifault multimechanism ruptures

We use the deterministic earthquake simulator RSQSim to generate complex sequences of ruptures on fault systems used for hazard assessment. We show that the source motions combined with a wave propagation code create surface ground motions that fall within the range of epistemic uncertainties for the Next Generation Attenuation‐West2 set of empirical models. We show the model is well calibrated where there are good data constraints, and has good correspondence in regions with fewer data constraints. We show magnitude, distance, and mechanism dependence all arising naturally from the same underlying friction. The deterministic physics‐based approach provides an opportunity for better understanding the physical origins of ground motions. For example, we find that reduced stress drops in shallow layers relative to constant stress drop with depth lead to peak ground velocities in the near field that better match empirical models. The simulators may also provide better extrapolations into regimes that are poorly empirically constrained by data because physics, rather than surface shaking data parameterizations, is underlying the extrapolations. Having shown the model is credible, we apply it to a problem where observations are lacking. We examine the case of crustal faults above a shallow subduction interface seen to break coseismically in simulations of the New Zealand fault system. These types of events were left out of consideration in the most recent New Zealand national seismic hazard model due to the modeling complexity and lack of observational data to constrain ground‐motion models (GMMs). Here, we show that in the model, by breaking up the coseismic crustal and interface rupturing fault motions into two separate subevents, and then recombining the resulting ground‐motion measures in a square‐root‐of‐sum‐of‐squares incoherent manner, we reproduce well the ground‐motion measures from the full event rupture. This provides a new method for extrapolating GMMs to more complex multifault ruptures.

Seismological Research Letters↗

USGS2023

Rapid surface rupture mapping from satellite data: The 2023 Kahramanmaraş, Turkey (Türkiye), earthquake sequence

The 6 February 2023 Kahramanmaraş, Turkey (Türkiye), earthquake sequence produced > 500 km of surface rupture primarily on the left‐lateral East Anatolian (~345 km) and Çardak (~175 km) faults. Constraining the length and magnitude of surface displacement on the causative faults is critical for loss estimates, recovery efforts, rapid identification of impacted infrastructure, and fault displacement hazard analysis. To support these efforts, we rapidly mapped the surface rupture from satellite data with support from remote sensing and field teams, and released the results to the public in near‐real time. Detailed surface rupture mapping commenced on 7 February and continued as high‐resolution (< 1.0 m/pixel) optical images from WorldView satellites (2023 Maxar) became available. We interpreted the initial simplified rupture trace from subpixel offset fields derived from Advanced Land Observation Satellite2 and Sentinel‐1A synthetic aperture radar image pairs available on 8 and 10 February, respectively. The mapping was released publicly on 10 February, with frequent updates, and published in final form four months postearthquake ( Reitman, Briggs, et al. , 2023 ). This publicly available, rapid mapping helped guide fieldwork and constrained U.S. Geological Survey finite‐fault and loss estimate models, as well as stress change estimates and dynamic rupture models.

Kahramanmaraş↗

USGS1988

Inversion for slip distribution using teleseismic P waveforms: North Palm Springs, Borah Peak, and Michoacan earthquakes

We have inverted the teleseismic P waveforms recorded by stations of the Global Digital Seismograph Network for the 8 July 1986 North Palm Springs, California, the 28 October 1983 Borah Peak, Idaho, and the 19 September 1985 Michoacan, Mexico, earthquakes to recover the distribution of slip on each of the faults using a point-by-point inversion method with smoothing and positivity constraints. In the inversion procedure, a fault plane with fixed strike and dip is placed in the region of the earthquake hypocenter and divided into a large number of subfaults. Rupture is assumed to propagate at a constant velocity away from the hypocenter, and synthetic ground motions for pure strike-slip and dip-slip dislocations are calculated at the teleseismic stations for each subfault. The observed seismograms are then inverted to obtain the distribution of strike-slip and dip-slip displacement for the earthquake.

California↗

USGS2005

Preliminary report on the 28 September 2004, M 6.0 Parkfield, California earthquake

The M w 6.0 Parkfield earthquake struck central California at 17:15:14 UTC on 28 September 2004. The epicenter was located 11 km southeast of the rural town of Parkfield, adjacent to Gold Hill and on the San Andreas Fault ( Figure 1 ). The California Integrated Seismic Network (CISN) reported that the hypocenter was located at 35.819°N, 120.364°W at a depth of 8.8 km. From the distribution of aftershocks and from models of seismograms, strain changes, and geodetic displacements from the earthquake, it appears that the rupture propagated to the northwest along the San Andreas Fault from its hypocenter beneath Gold Hill to Middle Mountain. Because of the earthquake's moderate size and the low population density, only minimal damage was reported, but strong ground motions of approximately 1 g were recorded at a few isolated points (Shakal et al. , 2005, this issue).

California↗

USGS1993

The San Andreas Fault

The presence of the San Andreas fault was brought dramatically to world attention on April 18, 1906, when sudden displacement along the fault produced the great San Francisco earthquake and fire. This earthquake, however, was but one of many that have resulted from episodic displacement along the fault throughout its life of about 15-20 million years.

California↗

USGS1967

Effects of the Truckee, California, earthquake of September 12, 1966

The Truckee, Calif., earthquake of September 12, 1966, had a magnitude of 5.4 on the Richter scale, as reported by the U.S. Coast and Geodetic Survey. It was felt over an area extending from San Francisco eastward to Salt Lake City and from Bakersfield northward beyond Chico, Calif. Ground breakage due to seismic shaking occurred in unconsolidated alluvium along a zone about 10 miles long that trends about N. 30° E. and extends from a point west of Prosser Reservoir, 4 miles northeast of Truckee, to Hoke Valley. The concentration of ground breakage along this zone suggests that it may be related to a subsurface northeastward-trending fault rather than to buried extensions of the northwestward-trending faults that dominate the regional geologic structure. Other effects of the earthquake include the following: (1) Damage to Prosser and Boca Dams, (2) minor damage to several bridges on Interstate Highway 80, (3) several earthquake-generated slumps or landslides on Interstate Highway 80 between Boca Dam and the California-Nevada boundary, (4) damage to the caretaker's house at Boca Dam and a lumber shed at Loyalton, (5) damage by falling rocks to the Farad powerhouse, to the wooden flume which supplies water to the powerhouse, and to the Southern Pacific Lines in Truckee Canyon, (6) minor damage to water wells and water distribution systems, particularly at Truckee, and (7) rejuvenation of springs throughout the Truckee-Russel Valley area.

Arizona, California, Idaho, Nevada, Oregon, Utah↗