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 289 records · Page 16Linked to original sources

USGS2017

3D ground‐motion simulations of Mw 7 earthquakes on the Salt Lake City segment of the Wasatch fault zone: Variability of long‐period (T≥1 s) ground motions and sensitivity to kinematic rupture parameters

We examine the variability of long‐period ( T ≥1 s) earthquake ground motions from 3D simulations of M w 7 earthquakes on the Salt Lake City segment of the Wasatch fault zone, Utah, from a set of 96 rupture models with varying slip distributions, rupture speeds, slip velocities, and hypocenter locations. Earthquake ruptures were prescribed on a 3D fault representation that satisfies geologic constraints and maintained distinct strands for the Warm Springs and for the East Bench and Cottonwood faults. Response spectral accelerations (SA; 1.5–10 s; 5% damping) were measured, and average distance scaling was well fit by a simple functional form that depends on the near‐source intensity level SA 0 ( T ) and a corner distance R c :SA( R , T )=SA 0 ( T )(1+( R / R c )) −1 . Period‐dependent hanging‐wall effects manifested and increased the ground motions by factors of about 2–3, though the effects appeared partially attributable to differences in shallow site response for sites on the hanging wall and footwall of the fault. Comparisons with modern ground‐motion prediction equations (GMPEs) found that the simulated ground motions were generally consistent, except within deep sedimentary basins, where simulated ground motions were greatly underpredicted. Ground‐motion variability exhibited strong lateral variations and, at some sites, exceeded the ground‐motion variability indicated by GMPEs. The effects on the ground motions of changing the values of the five kinematic rupture parameters can largely be explained by three predominant factors: distance to high‐slip subevents, dynamic stress drop, and changes in the contributions from directivity. These results emphasize the need for further characterization of the underlying distributions and covariances of the kinematic rupture parameters used in 3D ground‐motion simulations employed in probabilistic seismic‐hazard analyses.

Utah↗

USGS2021

Rupture passing probabilities at fault bends and steps, with application to rupture length probabilities for earthquake early warning

Earthquake early warning (EEW) systems can quickly identify the beginning of a significant earthquake rupture, but the first seconds of seismic data have not been found to predict the final rupture length. We present two approaches for estimating probabilities of rupture length given the rupture initiation from an EEW system. In the first approach, bends and steps on the fault are interpreted as physical mechanisms for rupture arrest. Arrest probability relations are developed from empirical observations and depend on bend angle and step size. Probability of arrest compounds serially with increasing rupture length as bends or steps are encountered. In the second approach, time‐independent rates among ruptures from the Uniform California Earthquake Rupture Forecast, Version 3 (UCERF3), are interpreted to apply to the time‐dependent condition in which rupture grows from a known starting point. Length probabilities from a Gutenberg–Richter magnitude–frequency relation provide a reference of comparison. We illustrate the new approach using the discretized fault model for California developed for UCERF3. For the case of rupture initiating on the southeast end of the San Andreas fault we find the geometric complexity of the Mill Creek section impedes most ruptures, and only ∼5% are predicted to reach to San Bernardino on the eastern edge of the greater Los Angeles region. Conditional probabilities of length can be precompiled in this manner for any initiation point on the fault system and thus are of potential value in seismic hazard and EEW applications.

California↗

USGS2026

Efficient physics‐informed ground‐motion simulations with reduced‐order models: CyberShake implications and high‐resolution site terms for southern San Andreas fault earthquakes

Recent advances in Probabilistic Seismic Hazard Analysis (PSHA) leverage physics‐based ground‐motion simulations to estimate seismic hazard, such as the CyberShake project. However, computational costs quickly escalate when performing PSHA for numerous faults or sites and can become prohibitively expensive. To reduce computational demands, CyberShake uses reciprocity and interpolates physics‐informed corrections from simulations conducted at fewer locations, but the accuracy of these interpolations remains poorly quantified. To quantify the interpolation accuracy, we derive high‐resolution, frequency‐dependent site terms for southern California and compare them with interpolated site terms using the CyberShake approach. We accomplish this by performing a set of earthquake point‐source simulations distributed along the nonplanar fault geometry for the southern San Andreas fault (SSAF) extending from Bombay Beach to Lake Hughes. Using SeisSol, we simulate three minutes of viscoelastic seismic wave propagation for these sources and store the horizontal‐component Green’s functions for 480,000 sites. We then use a scientific machine learning approach based on interpolated proper orthogonal decomposition to construct an accurate reduced‐order model of the Green’s functions to efficiently predict effective amplitude spectra (EAS) for finite‐source rupture models of SSAF earthquakes. Using minimum curvature interpolation with tension, as used in CyberShake, we compare the interpolated site terms against our high‐resolution site terms. We identify local discrepancies with EAS differing by up to a factor of approximately three. Furthermore, we identify locations where unexpectedly high or low ground motions are missed when using the interpolated dataset for these earthquakes. We estimate that our approach may be used within CyberShake to reduce the time‐to‐solution by a factor of 336 for the entire earthquake rupture forecast. Our analysis of physics‐based site terms provides more insight into the seismic hazard due to SSAF ruptures and guides future developments by combining high‐performance computing and reduced‐order modeling techniques for PSHA.

California↗

USGS2026

3D broadband ground motion simulations of Cascadia megathrust earthquakes: Full-margin scenarios incorporating variable magnitudes and frontal thrust splay faulting

Great megathrust earthquakes (∼ M w 8–9) that rupture along the Cascadia Subduction Zone pose significant shaking and tsunami hazards to the U.S. Pacific Northwest. Although these large earthquakes are documented in the geologic record, none have been recorded with modern instrumentation, complicating regional seismic hazard assessments. To improve constraints on seismic hazard estimates for future earthquakes, we develop 3D broadband (0–10 Hz) ground‐motion simulations for 36 full‐margin megathrust earthquake scenarios through the Cascadia Coastlines and Peoples Hazards Research Hub. These simulations expand on the M9 Project ( Frankel et al. , 2018 ; Wirth, Frankel, et al. , 2018 ) by incorporating a broader range of rupture parameters—variable magnitudes, splay faulting, and geodetically informed slip distributions. Benchmarking these simulations against the Next Generation Attenuation Subduction (NGA‐Sub) ground‐motion models shows reasonable fits for non‐basin sites, although simulations predict stronger shaking at short (0–2 s) and long (6–10 s) periods. NGA‐Sub basin terms improve fits to synthetic shaking within the Seattle Basin but underpredict amplitudes by a factor of 2 at periods of 2–5 s. We also evaluate how rupture characteristics influence regional inland shaking variability and basin response in the Puget Lowlands. Larger magnitude scenarios generally cause stronger inland shaking due to farther inland extents of slip and larger strong‐motion‐generating fault patches but show little difference in basin amplification within the Puget Lowlands compared to smaller magnitude events. Along‐dip hypocenter location strongly influences overall inland shaking, with offshore (updip) hypocenters generating on average ∼45% higher shaking at long periods (>6 s) compared to onshore (downdip) hypocenters due to rupture directivity effects. Offshore splay faults minimally affect inland shaking but alter offshore deformation important for coupled ground‐shaking and tsunami scenarios. Finally, we show that rupture directivity increases amplification within the Seattle Basin through the generation and interaction of higher order mode surface waves, an effect of 3D basin geometry that cannot be captured using 1D approximations.

British Columbia, California, Oregon, Washington↗

USGS2025

Potential for continental scientific drilling to inform fault mechanics and earthquake science

Our understanding of fault mechanics and earthquake processes remains limited, largely due to minimal direct observations near active faults at seismogenic depths. This lack of data restricts our ability to accurately assess and mitigate both natural and human-induced seismic hazards. However, recent advancements in drilling capabilities and downhole sensing technologies offer an opportunity: the ability to observe the physical conditions within a volume near active fault zones. In this contribution, we highlight how scientific drilling can provide access to the near-fault environment, enabling measurements of the stress, temperature, fluid pressure, and rock properties at depths where ruptures initiate, propagate, and arrest. These observations are essential to refine models of earthquake nucleation and dynamic rupture, bridging gaps between laboratory experiments, numerical simulations, and surface observations. These insights can advance fundamental understanding in earthquake science but also support the development of more effective seismic hazard assessments and risk mitigation strategies.

Seismica↗

USGS2010

Predicted liquefaction in the greater Oakland area and northern Santa Clara Valley during a repeat of the 1868 Hayward Fault (M6.7-7.0) earthquake

Probabilities of surface manifestations of liquefaction due to a repeat of the 1868 ( M 6.7-7.0) earthquake on the southern segment of the Hayward Fault were calculated for two areas along the margin of San Francisco Bay, California: greater Oakland and the northern Santa Clara Valley. Liquefaction is predicted to be more common in the greater Oakland area than in the northern Santa Clara Valley owing to the presence of 57 km 2 of susceptible sandy artificial fill. Most of the fills were placed into San Francisco Bay during the first half of the 20th century to build military bases, port facilities, and shoreline communities like Alameda and Bay Farm Island. Probabilities of liquefaction in the area underlain by this sandy artificial fill range from 0.2 to ~0.5 for a M 7.0 earthquake, and decrease to 0.1 to ~0.4 for a M 6.7 earthquake. In the greater Oakland area, liquefaction probabilities generally are less than 0.05 for Holocene alluvial fan deposits, which underlie most of the remaining flat-lying urban area. In the northern Santa Clara Valley for a M 7.0 earthquake on the Hayward Fault and an assumed water-table depth of 1.5 m (the historically shallowest water level), liquefaction probabilities range from 0.1 to 0.2 along Coyote and Guadalupe Creeks, but are less than 0.05 elsewhere. For a M 6.7 earthquake, probabilities are greater than 0.1 along Coyote Creek but decrease along Guadalupe Creek to less than 0.1. Areas with high probabilities in the Santa Clara Valley are underlain by young Holocene levee deposits along major drainages where liquefaction and lateral spreading occurred during large earthquakes in 1868 and 1906.

California↗

USGS1987

Fault failure with moderate earthquakes

High resolution strain and tilt recordings were made in the near-field of, and prior to, the May 1983 Coalinga earthquake ( M L = 6.7, Δ = 51 km), the August 4, 1985, Kettleman Hills earthquake ( M L = 5.5, Δ = 34 km), the April 1984 Morgan Hill earthquake ( M L = 6.1, Δ = 55 km), the November 1984 Round Valley earthquake ( M L = 5.8, Δ = 54 km), the January 14, 1978, Izu, Japan earthquake ( M L = 7.0, Δ = 28 km), and several other smaller magnitude earthquakes. These recordings were made with near-surface instruments (resolution 10 −8 ), with borehole dilatometers (resolution 10 −10 ) and a 3-component borehole strainmeter (resolution 10 −9 ). While observed coseismic offsets are generally in good agreement with expectations from elastic dislocation theory, and while post-seismic deformation continued, in some cases, with a moment comparable to that of the main shock, preseismic strain or tilt perturbations from hours to seconds (or less) before the main shock are not apparent above the present resolution. Precursory slip for these events, if any occurred, must have had a moment less than a few percent of that of the main event. To the extent that these records reflect general fault behavior, the strong constraint on the size and amount of slip triggering major rupture makes prediction of the onset times and final magnitudes of the rupture zones a difficult task unless the instruments are fortuitously installed near the rupture initiation point. These data are best explained by an inhomogeneous failure model for which various areas of the fault plane have either different stress-slip constitutive laws or spatially varying constitutive parameters. Other work on seismic waveform analysis and synthetic waveforms indicates that the rupturing process is inhomogeneous and controlled by points of higher strength. These models indicate that rupture initiation occurs at smaller regions of higher strength which, when broken, allow runaway catastrophic failure.

California↗

USGS1996

Role of stress triggering in earthquake migration on the North Anatolian fault

Ten M???6.7 earthquakes ruptured 1,000 km of the North Anatolian fault (Turkey) during 1939-92, providing an unsurpassed opportunity to study how one large shock sets up the next. Calculations of the change in Coulomb failure stress reveal that 9 out of 10 ruptures were brought closer to failure by the preceding shocks, typically by 5 bars, equivalent to 20 years of secular stressing. We translate the calculated stress changes into earthquake probabilities using an earthquake-nucleation constitutive relation, which includes both permanent and transient stress effects. For the typical 10-year period between triggering and subsequent rupturing shocks in the Anatolia sequence, the stress changes yield an average three-fold gain in the ensuing earthquake probability. Stress is now calculated to be high at several isolated sites along the fault. During the next 30 years, we estimate a 15% probability of a M???6.7 earthquake east of the major eastern center of Erzincan, and a 12% probability for a large event south of the major western port city of Izmit. Such stress-based probability calculations may thus be useful to assess and update earthquake hazards elsewhere. ?? 1997 Elsevier Science Ltd.

Physics and Chemistry of the Earth↗

USGS2002

Response of the San Andreas fault to the 1983 Coalinga-Nuñez earthquakes: an application of interaction-based probabilities for Parkfield

The Parkfield-Cholame section of the San Andreas fault, site of an unfulfilled earthquake forecast in 1985, is the best monitored section of the world's most closely watched fault. In 1983, the M = 6.5 Coalinga and M = 6.0 Nuñez events struck 25 km northeast of Parkfield. Seismicity rates climbed for 18 months along the creeping section of the San Andreas north of Parkfield and dropped for 6 years along the locked section to the south. Right-lateral creep also slowed or reversed from Parkfield south. Here we calculate that the Coalinga sequence increased the shear and Coulomb stress on the creeping section, causing the rate of small shocks to rise until the added stress was shed by additional slip. However, the 1983 events decreased the shear and Coulomb stress on the Parkfield segment, causing surface creep and seismicity rates to drop. We use these observations to cast the likelihood of a Parkfield earthquake into an interaction-based probability, which includes both the renewal of stress following the 1966 Parkfield earthquake and the stress transfer from the 1983 Coalinga events. We calculate that the 1983 shocks dropped the 10-year probability of a M &sim; 6 Parkfield earthquake by 22% (from 54 ± 22% to 42 ± 23%) and that the probability did not recover until about 1991, when seismicity and creep resumed. Our analysis may thus explain why the Parkfield earthquake did not strike in the 1980s, but not why it was absent in the 1990s. We calculate a 58 ± 17% probability of a M &sim; 6 Parkfield earthquake during 2001–2011.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

Imaging the San Andreas fault with explosion and earthquake sources

Mounting evidence suggests that fault zone heterogeneity may play a crucial role in the localization of rupture in earthquakes [ Aki , 1995]. The heterogeneity can take several forms: spatial variations in physical properties (elastic properties, pore fluid pressure, etc.) or complexity in the fault surface (bends, offsets, etc.). High-resolution, three-dimensional models of the P and S wave velocity ( V p and V s ) structure and accurate hypocenters provide tools for studying the effects of fault zone heterogeneity on rupture localization. A dense, passive array of 48 seismic instruments was deployed south of Hollister, California, from mid-November 1994 to late May 1995. The array obtained local-earthquake seismograms for high-resolution, three-dimensional imaging of the ( V p and V s ) structure in the San Andreas fault zone. In mid-May, anactive seismic experiment was also carried out. Figure 1 shows a 7-km wide, 3- to 4-km deep zone of very low V p imaged during the experiment. The zone is bounded on the southwest by the San Andreas fault and the adjacent Northern Gabilan Range. The San Andreas fault appears to be approximately vertical to at least 5 km depth (the limit of good model resolution), contrary to an existing model based on gravity data.

California↗

USGS1998

Evidence of shallow fault zone strengthening after the 1992 M7.5 Landers, California, earthquake

Repeated seismic surveys of the Landers, California, fault zone that ruptured in the magnitude (M) 7.5 earthquake of 1992 reveal an increase in seismic velocity with time. P, S, and fault zone trapped waves were excited by near-surface explosions in two locations in 1994 and 1996, and were recorded on two linear, three-component seismic arrays deployed across the Johnson Valley fault trace. The travel times of P and S waves for identical shot-receiver pairs decreased by 0.5 to 1.5 percent from 1994 to 1996, with the larger changes at stations located within the fault zone. These observations indicate that the shallow Johnson Valley fault is strengthening after the main shock, most likely because of closure of cracks that were opened by the 1992 earthquake. The increase in velocity is consistent with the prevalence of dry over wet cracks and with a reduction in the apparent crack density near the fault zone by approximately 1.0 percent from 1994 to 1996.

Science↗

USGS2002

A record of large earthquakes on the southern Hayward fault for the past 500 years

The Hayward fault, a major branch of the right-lateral San Andreas fault system, traverses the densely populated eastern San Francisco Bay region, California. We conducted a paleoseismic investigation to better understand the Hayward fault's past earthquake behavior. The site is near the south end of Tyson's Lagoon, a sag pond formed in a right step of the fault in Fremont. Because the Hayward fault creeps at the surface, we identified paleoseismic events using features that we judge to be unique to ground ruptures or the result of strong ground motion, such as the presence of fault-scarp colluvial deposits and liquefaction. We correlate the most recent event evidence (E1) to the historical 1868 M 6.9 earthquake that caused liquefaction in the pond and recognize three additional paleoruptures since A.D. 1470 ± 110 yr. Event ages were estimated by chronological modeling, which incorporated historical and stratigraphic information and radiocarbon and pollen data. Modeled, mean age and 95-percentile ranges of the three earlier events are A.D. 1730 (1650-1790) yr (E2), A.D. 1630 (1530-1740) yr (E3), and A.D. 1470 (1360-1580) (E4). The ages of these paleoearthquakes yield a mean recurrence of 130 ± 40 yr. Although the mean recurrence is well determined for the period A.D. 1470-1868, individual intervals are less well determined: E1-E2, 140 +80/-70 yr; E2-E3, 100 +90/-100 yr; and E3-E4, 150 +130/-110 yr.

California↗

USGS2002

Structure and earthquake offsets on the San Andreas fault at the Wrightwood, California paleoseismic site

Thirty-eight trenches and natural exposures across the San Andreas fault zone, four km northwest of Wrightwood, California, reveal the structure and stratigraphy of this paleoseismic site. A 25-m-thick stratigraphic section of late Holocene peat, debris flows, and fluvial sediments have accumulated syntectonically and are deformed in a complicated flower, or flake, structure within the fault zone. Whereas slip is everywhere parallel to the trend of the San Andreas fault, layering within the fan-shaped section guides the surface faulting, and transtensive portions of the structure are connected to transpressive portions by bedding-parallel detachment faults. The structure has grown through time, widely distributing the evidence for individual earthquakes. The broad zone of deformation and rapid sedimentation permit the characterization of an unusually large number of prehistoric earthquakes. The presence of at least four geomorphically similar sites within 10 km along the fault zone suggests that this kind of site is common and other sites should be exploited.

California↗

USGS2003

Maximum slip in earthquake fault zones, apparent stress, and stick-slip friction

The maximum slip , observed or inferred, for a small patch within the larger fault zone of an earthquake is a remarkably well-constrained function of the seismic moment. A large set of maximum slips, mostly derived from slip models of major earthquakes, indicate that this parameter increases according to the cube root of the seismic moment. Consistent with this finding, neither the average slip rate for the patches of maximum slip nor the apparent stresses of earthquakes show any systematic dependence on seismic moment. Maximum average slip rates are several meters per second independent of moment and, for earthquakes in continental crustal settings, the apparent stress is limited to about 10 MPa. Results from stick - slip friction experiments in the laboratory, combined with information about the state of stress in the crust, can be used to predict, quite closely, the maximum slips and maximum average slip rates within the fault zones of major earthquakes as well as their apparent stresses. These findings suggest that stick - slip friction events observed in the laboratory and earthquakes in continental settings, even with large magnitudes, have similar rupture mechanisms.

Bulletin of the Seismological Society of America↗

USGS1990

Source inversion of the 1988 Upland, California, earthquake: Determination of a fault plane for a small event

We examined short-period P waves to investigate if waveform data could be used to determine which of two nodal planes was the actual fault plane for a small (ML 4.6) earthquake near Upland, California. We removed path and site complications by choosing a small aftershock (ML 2.7) as an empirical Green function. The main shock P waves were deconvolved by using the empirical Green function to produce simple far-field displacement pulses. We used a least-squares method to invert these pulses for the slip distribution on a finite fault. Both nodal planes (strike 125°, dip 85° and strike 221°, dip 40°) of the first-motion focal mechanism were tested at various rupture velocities. The southwest trending fault plane consistently gave better fitting solutions than the southeast-trending plane. We determined a moment of 4.2 × 1022 dyne-cm. The rupture velocity, and thus the source area could not be well resolved, but if we assume a reasonable rupture velocity of 0.87 times the shear wave velocity, we obtain a source area of 0.97 km2 and a stress drop of 38 bars. Choice of a southwest-trending fault plane is consistent with the trend of the nearby portion of the Transverse Ranges frontal fault zone and indicates left-lateral motion. This method provides a way to determine the fault plane for small earthquakes that have no surface rupture and no obvious trend in aftershock locations.

California↗

USGS1994

Shear zones formed along long, straight traces of fault zones during the 28 June 1992 Landers, California, earthquake

Surface rupturing during the 28 June 1992 Landers, California, earthquake, east of Los Angeles, accommodated right-lateral offsets up to about 6 m along segments of distinct, en-echelon fault zones with a total length of 80 km. The offsets were accommodated generally not by faults—distinct slip surfaces—but rather by shear zones , tabular bands of localized shearing. Along simple stretches of fault zones at Landers the rupture is characterized by telescoping of shear zones and intensification of shearing: broad shear zones of mild shearing, containing narrow shear zones of more intense shearing, containing even narrower shear zones of very intense shearing, which may contain a fault. Thus the ground ruptured across broad belts of shearing with clearly defined, subparallel walls, oriented NW. Each broad belt consists of a broad zone of mild shearing, extending across its entire width (50 to 200 m), and much narrower (a few meters wide) shear zones that accommodate most of the offset of the belt and are portrayed by en-echelon tension cracks. In response to right-lateral shearing, the slices of ground bounded by the tension cracks rotated in a clockwise sense, producing left-lateral shearing, and the slices were forced against the walls of the shear zone, producing thrusting. Even narrower shear zones formed within the narrow shear zones. Although these probably are guides to right-lateral fault segments below, the surface rupturing during the earthquake is characterized not by faulting, but by the formation of shear zones at various scales.

California↗

USGS2004

Landslides triggered by the 2002 Denali fault, Alaska, earthquake and the inferred nature of the strong shaking

The 2002 M7.9 Denali fault, Alaska, earthquake triggered thousands of landslides, primarily rock falls and rock slides, that ranged in volume from rock falls of a few cubic meters to rock avalanches having volumes as great as 15 ?? 106 m3. The pattern of landsliding was unusual; the number of slides was less than expected for an earthquake of this magnitude, and the landslides were concentrated in a narrow zone 30-km wide that straddled the fault rupture over its entire 300-km length. The large rock avalanches all clustered along the western third of the rupture zone where acceleration levels and ground-shaking frequencies are thought to have been the highest. Inferences about near-field strong shaking characteristics drawn from the interpretation of the landslide distribution are consistent with results of recent inversion modeling that indicate high-frequency energy generation was greatest in the western part of the fault rupture zone and decreased markedly to the east. ?? 2004, Earthquake Engineering Research Institute.

Earthquake Spectra↗

USGS2011

The ShakeOut scenario: A hypothetical Mw7.8 earthquake on the Southern San Andreas Fault

In 2008, an earthquake-planning scenario document was released by the U.S. Geological Survey (USGS) and California Geological Survey that hypothesizes the occurrence and effects of a M w 7.8 earthquake on the southern San Andreas Fault. It was created by more than 300 scientists and engineers. Fault offsets reach 13 m and up to 8 m at lifeline crossings. Physics-based modeling was used to generate maps of shaking intensity, with peak ground velocities of 3 m/sec near the fault and exceeding 0.5 m/sec over 10,000 km 2 . A custom HAZUS ®MH analysis and 18 special studies were performed to characterize the effects of the earthquake on the built environment. The scenario posits 1,800 deaths and 53,000 injuries requiring emergency room care. Approximately 1,600 fires are ignited, resulting in the destruction of 200 million square feet of the building stock, the equivalent of 133,000 single-family homes. Fire contributes $ 87 billion in property and business interruption loss, out of the total $ 191 billion in economic loss, with most of the rest coming from shake-related building and content damage ( $ 46 billion) and business interruption loss from water outages ( $ 24 billion). Emergency response activities are depicted in detail, in an innovative grid showing activities versus time, a new format introduced in this study.

Calfornia↗