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

USGS2006

Using Modified Mercalli Intensities to estimate acceleration response spectra for the 1906 San Francisco earthquake

We derive and test relations between the Modified Mercalli Intensity (MMI) and the pseudo-acceleration response spectra at 1.0 and 0.3 s - SA(1.0 s) and SA(0.3 s) - in order to map response spectral ordinates for the 1906 San Francisco earthquake. Recent analyses of intensity have shown that MMI ??? 6 correlates both with peak ground velocity and with response spectra for periods from 0.5 to 3.0 s. We use these recent results to derive a linear relation between MMI and log SA(1.0 s), and we refine this relation by comparing the SA(1.0 s) estimated from Boatwright and Bundock's (2005) MMI map for the 1906 earthquake to the SA(1.0 s) calculated from recordings of the 1989 Loma Prieta earthquake. South of San Jose, the intensity distributions for the 1906 and 1989 earthquakes are remarkably similar, despite the difference in magnitude and rupture extent between the two events. We use recent strong motion regressions to derive a relation between SA(1.0 s) and SA(0.3 s) for a M7.8 strike-slip earthquake that depends on soil type, acceleration level, and source distance. We test this relation by comparing SA(0.3 s) estimated for the 1906 earthquake to SA(0.3 s) calculated from recordings of both the 1989 Loma Prieta and 1994 Northridge earthquakes, as functions of distance from the fault. ?? 2006, Earthquake Engineering Research Institute.

Earthquake Spectra↗

USGS1993

The October 12, 1992, Dahshur, Egypt, Earthquake

Cairo and northeastern Egypt experienced a rare, damaging earthquake on October 12, 1992. The earthquake, which measured 5.9 on the Richter magnitude scale, was centered near the village of Dahshur, about 18 km south of Cairo. The computed hypocentral depth of the earthquake, about 25 km, is consistent with the fact that fault rupture associated with the earthquake did not reach the surface. Despite its relatively moderate size, the earthquake caused many casualties and heavy damage. These losses included more than 500 fatalities, more than 6,500 injuries, and about 8,300 damaged or destroyed buildings. The Foreign Broadcast Service estimated monetary losses directly attributable to the earthquake at $300 million. We were part of an international reconnaissance team that investigated the Dahsur earthquake. This article summarizes our findings and points out how even a relatively moderate sized earthquake can cause widespread damage and a large number of casualities.

Earthquakes & Volcanoes (USGS)↗

USGS2014

High-frequency imaging of elastic contrast and contact area with implications for naturally observed changes in fault properties

During localized slip of a laboratory fault we simultaneously measure the contact area and the dynamic fault normal elastic stiffness. One objective is to determine conditions where stiffness may be used to infer changes in area of contact during sliding on nontransparent fault surfaces. Slip speeds between 0.01 and 10 µm/s and normal stresses between 1 and 2.5 MPa were imposed during velocity step, normal stress step, and slide-hold-slide tests. Stiffness and contact area have a linear interdependence during rate stepping tests and during the hold portion of slide-hold-slide tests. So long as linearity holds, measured fault stiffness can be used on nontransparent materials to infer changes in contact area. However, there are conditions where relations between contact area and stiffness are nonlinear and nonunique. A second objective is to make comparisons between the laboratory- and field-measured changes in fault properties. Time-dependent changes in fault zone normal stiffness made in stress relaxation tests imply postseismic wave speed changes on the order of 0.3% to 0.8% per year in the two or more years following an earthquake; these are smaller than postseismic increases seen within natural damage zones. Based on scaling of the experimental observations, natural postseismic fault normal contraction could be accommodated within a few decimeter wide fault core. Changes in the stiffness of laboratory shear zones exceed 10% per decade and might be detectable in the field postseismically.

Journal of Geophysical Research↗

USGS2005

The Olmsted fault zone, southernmost Illinois: A key to understanding seismic hazard in the northern new Madrid seismic zone

Geological deformation in the northern New Madrid seismic zone, near Olmsted, Illinois (USA), is analyzed using integrated compressional-wave (P) and horizontally polarized-wave (SH) seismic reflection and regional and dedicated borehole information. Seismic hazards are of special concern because of strategic facilities (e.g., lock and dam sites and chemical plants on the Ohio River near its confluence with the Mississippi River) and because of alluvial soils subject to high amplification of earthquake shock. We use an integrated approach starting with lower resolution, but deeper penetration, P-wave reflection profiles to identify displacement of Paleozoic bedrock. Higher resolution, but shallower penetration, SH-wave images show deformation that has propagated upward from bedrock faults into Pleistocene loess. We have mapped an intricate zone more than 8 km wide of high-angle faults in Mississippi embayment sediments localized over Paleozoic bedrock faults that trend north to northeast, parallel to the Ohio River. These faults align with the pattern of epicenters in the New Madrid seismic zone. Normal and reverse offsets along with positive flower structures imply a component of strike-slip; the current stress regime favors right-lateral slip on northeast-trending faults. The largest fault, the Olmsted fault, underwent principal displacement near the end of the Cretaceous Period 65 to 70 million years ago. Strata of this age (dated via fossil pollen) thicken greatly on the downthrown side of the Olmsted fault into a locally subsiding basin. Small offsets of Tertiary and Quaternary strata are evident on high-resolution SH-wave seismic profiles. Our results imply recent reactivation and possible future seismic activity in a critical area of the New Madrid seismic zone. This integrated approach provides a strategy for evaluating shallow seismic hazard-related targets for engineering concerns. ?? 2005 Elsevier B.V. All rights reserved.

Engineering Geology↗

USGS2016

Active tectonics within the NW and SE extensions of the Pambak-Sevan-Syunik fault: Implications for the present geodynamics of Armenia

This study analyzes the active tectonics within the northwestern and southeastern extensions of the Pambak-Sevan-Syunik fault (PSSF), a major right-lateral strike-slip fault cutting through Armenia. Quantifying the deformations in terms of geometry, kinematics, slip rates and earthquake activity, using cosmogenic 3 He, OSL/IRSL and radiocarbon dating techniques, reveal different behaviors between the two regions. Within the northwestern extension, in the region of Amasia, the PSSF bends to the west and splits into two main WNW–ESE trending reverse faults defining a compressional pop-up structure. We estimate an uplift rate and a shortening rate of 0.5 ± 0.1 mm/y and 1.4 ± 0.6 mm/y, respectively. This suggests that most of the &sim;2 mm/y right lateral movement of the PSSF seems to be absorbed within the Amasia pop-structure. Within the southeastern extension, the PSSF shows signs of dying out within the Tsghuk Volcano region at the southernmost tip of the Syunik graben. There, the tectonic activity is characterized by a very slow NS trending normal faulting associated with a slight right-lateral movement. Slip rates analyses (i.e. vertical slip rate, EW stretching rate at 90° to the fault, and right-lateral slip rate of &sim;0.2 mm/y, &sim;0.1 mm/y and &sim;0.05 mm/y, respectively) lead to the conclusion that the right lateral movement observed further north along the PSSF is mainly transferred within other active faults further west within the Karabagh (Hagari fault or other structures further northwestwards). Comparing our slip rates with those estimated from GPS data suggests that most of the deformation is localized and seismic, at least within the Tsghuk region. The geometrical and kinematic pattern observed within the two terminations of the PSSF suggests that the fault and its surrounding crustal blocks are presently rotating anticlockwise, as also observed within the GPS velocity field. This is consistent with the recent kinematic models proposed for the Caucasus-Kura-South Caspian region and brings a new insight into the present geodynamics of Armenia.

Quaternary International↗

USGS1998

The influence of the San Gregorio fault on the morphology of Monterey Canyon

A side-scan sonar survey was conducted of Monterey Canyon and the San Gregorio fault zone, off shore of Monterey Bay. The acoustic character and morphology of the sonar images, enhanced by SeaBeam bathymetry, show the path of the San Gregorio fault zone across the shelf, upper slope, and Monterey Canyon. High backscatter linear features a few kilometers long and 100 to 200 m wide delineate the sea-floor expression of the fault zone on the shelf. Previous studies have shown that brachiopod pavements and carbonate crusts are the source of the lineations backscatter. In Monterey Canyon, the fault zone occurs where the path of the canyon makes a sharp bend from WNW to SSW (1800 m). Here, the fault is marked by NW-SE-trending, high reflectivity lineations that cross the canyon floor between 1850 m and 1900 m. The lineations can be traced to ridges on the northwestern canyon wall where they have ~ 15 m of relief. Above the low-relief ridges, bowl-shaped features have been excavated on the canyon wall contributing to the widening of the canyon. We suggest that shear along the San Gregorio fault has led to the formation of the low-relief ridges near the canyon wall and that carbonate crusts, as along the shelf, may be the source of the high backscatter features on the canyon floor. The path of the fault zone across the upper slope is marked by elongated tributary canyons with high backscatter floors and 'U'-shaped cross-sectional profiles. Linear features and stepped scarps suggestive of recent crustal movement and mass-wasting, occur on the walls and floors of these canyons. Three magnitude-4 earthquakes have occurred within the last 30 years in the vicinity of the canyons that may have contributed to the observed features. As shown by others, motion along the fault zone has juxtaposed diverse lithologies that outcrop on the canyon walls. Gully morphology and the canyon's drainage patterns have been influenced by the substrate into which the gullies have formed.

Marine Geology↗

USGS2000

Recent tectonics of the Blanco Ridge, eastern blanco transform fault zone

Bathymetric, hydro-acoustic, seismic, submersible, and gravity data are used to investigate the active tectonics of the eastern Blanco Transform Fault Zone (BTFZ). The eastern BTFZ is dominated by the ∼150 km long transform-parallel Blanco Ridge (BR) which is a right-lateral strike-slip fault bordered to the east and west by the Gorda and Cascadia Depressions. Acoustic locations, fault-parameter information, and slip vector estimates of 43 earthquakes ( M w ≥3.8) that occurred along the eastern BTFZ over the last 5 years reveal that the Blanco Ridge is a high-angle right-lateral strike-slip fault, with a small component of dip-slip motion, where the Juan de Fuca plate is the hanging wall relative to the Pacific plate. Furthermore, the Cascadia and Gorda basins are undergoing normal faulting with extension predominantly oblique to the transform trend. Seafloor submersible observations agree with previous hypotheses that the active transform fault trace is the elongate basin that runs the length of the BR summit. Brecciated and undeformed basalt, diabase, and gabbro samples were collected at the four submersible survey sites along the Blanco Ridge. These petrologic samples indicate the Blanco Ridge is composed of an ocean crustal sequence that has been uplifted and highly fractured. The petrologic samples also appear to show an increase in elevation of the crustal section from east to west along the Blanco Ridge, with gabbros exposed at a shallower depth farther west along the southern (Pacific plate side) BR ridge flank. Further supporting evidence for BR uplift exists in the seismic reflection profiles across the BR showing uplift of turbidite sequences along the north and south ridge base, and gravity and magnetics profiles that indicate possible basement uplift and a low-density zone centered on the ridge's Pacific plate side. The BR formation mechanism preferred here is first, uplift achieved partially through strike-slip motion (with a small dip-slip component). Second, seawater penetration along the fault into the lower crust upper mantle, which then enhanced formation and intrusion of a mantle-derived serpentinized-peridotite diapir into the shallow ocean crust, causing further uplift along the fault.

Marine Geophysical Research↗

USGS2023

Rift basins and intraplate earthquakes: New high-resolution aeromagnetic data provide insights into buried structures of the Charleston, South Carolina seismic zone

The delineation of faults that pose seismic risk in intraplate seismic zones and the mapping of features associated with failed rift basins can help our understanding of links between the two. We use new high-resolution aeromagnetic data, previous borehole sample information, and reprocessed seismic reflection profiles to image subsurface structures and evaluate recent fault activity within the Charleston seismic zone, the associated Mesozoic South Georgia rift basin, and surrounds. The new aeromagnetic data provide an unprecedented view of buried basement structures. NE- and NW-trending lineaments of various lengths throughout the survey area are interpreted as Paleozoic orogenic structures and Mesozoic dikes, respectively. Within the rift basin, 15- to 20-km long ESE-trending lineaments are associated with faults in pre-Cretaceous strata of the reflection data and are interpreted as Mesozoic rift structures. Various intersections and terminations of interpreted faults suggest rift-related reactivation of Paleozoic faults and corresponding inheritance for Mesozoic structures. The reflection data show that several Paleozoic and Mesozoic faults are associated with deformation in Cretaceous and younger sediments, suggesting reactivation in the more recent passive margin setting. Two of these faults, one NE-striking and one ESE-striking, are coincident with surficial landforms, suggesting Quaternary slip; the ESE-striking fault is also well-aligned with a plan-view offset in modern seismicity. A favorable orientation for reverse motion on ESE-striking Mesozoic faults, a possible sub-basin, and potentially weakened lithosphere are failed rift basin features that may influence intraplate seismicity within the Charleston seismic zone.

South Carolina↗

USGS2011

Earthquake prediction lessons from Parkfield Experiment

In 1985, the National Earthquake Prediction Evaluation Council (NEPEC), and later its California counterpart (CEPEC), approved an earthquake prediction for the first time. This was the Bakun and Lindh [1985] forecast, with 95% confidence, that an M = 5.5 to 6 earthquake would rupture along the San Andreas fault near the town of Parkfield before 1993. This forecast was an important factor in selecting the Parkfield area for a major geophysical and geological monitoring study. The experiments installed there are aimed at producing a better understanding of the earthquake generation process and earthquake effects, gaining experience with real-time monitoring of a variety of geophysical data, and, if possible, making a short-term prediction of the future earthquake. Twenty-eight different observational networks now exist in the Parkfield area. The forecast also generated interest that led researchers to reexamine past Parkfield earthquakes with current knowledge and to develop new models of the earthquake cycle.

California↗

USGS2012

Seismomagnetic effect generated by the October, 1989, ML, 7.1 Loma Prieta, California, Earthquake

A differentially connected array of proton magnetometers operated within the epicentral region of the October 18, 1989, M L 7.1 Loma Prieta earthquake for 12 years from 1974 to 1986. The closest magnetometer station was located 7.3 km from the epicenter of the earthquake and within 3 km of the site where anomalous ULF magnetic noise measurements were observed. Following the earthquake, the magnetometers were reinstalled with sensors replaced in the original undisturbed sensor holders. Comparison of pre-1986 total intensity magnetic field data with data obtained during the months following the earthquake indicate local offsets of about 1 nT may have been generated at stations nearest the epicenter. Tests on other continuous differenced data from 1983 to present indicate that offsets determined could be biased by as much as 0.7 nT. The offsets can be approximately fit with a simple seismomagnetic model of the earthquake for which 1.9 m of right lateral and 1.3 m of dip slip (southwest side up) occurred on a fault patch between 6 km and 18 km deep and 45 km long. The total rock magnetization is assumed to be 1.5 A/m. Since the offset has persisted following the earthquake, an alternate explanation in terms of electokinetic effects is unlikely even though transient ground water flow occurred following the earthquake. Comparison of pre-1986 and similar post-seismic total magnetic field noise does not indicate any change caused by aliasing of ULF (0.01 Hz-10 Hz) magnetic noise in the vicinity of the Loma Prieta earthquake.

California↗

USGS2015

San Andreas tremor cascades define deep fault zone complexity

Weak seismic vibrations - tectonic tremor - can be used to delineate some plate boundary faults. Tremor on the deep San Andreas Fault, located at the boundary between the Pacific and North American plates, is thought to be a passive indicator of slow fault slip. San Andreas Fault tremor migrates at up to 30 m s - 1 , but the processes regulating tremor migration are unclear. Here I use a 12-year catalogue of more than 850,000 low-frequency earthquakes to systematically analyse the high-speed migration of tremor along the San Andreas Fault. I find that tremor migrates most effectively through regions of greatest tremor production and does not propagate through regions with gaps in tremor production. I interpret the rapid tremor migration as a self-regulating cascade of seismic ruptures along the fault, which implies that tremor may be an active, rather than passive participant in the slip propagation. I also identify an isolated group of tremor sources that are offset eastwards beneath the San Andreas Fault, possibly indicative of the interface between the Monterey Microplate, a hypothesized remnant of the subducted Farallon Plate, and the North American Plate. These observations illustrate a possible link between the central San Andreas Fault and tremor-producing subduction zones.

California↗

USGS2020

Mechanics of near-field deformation during co- and post-seismic shallow fault slip

Poor knowledge of how faults slip and distribute deformation in the shallow crust hinders efforts to mitigate hazards where faults increasingly intersect with the expanding global population at Earth’s surface. Here we analyze two study sites along the 2014 M 6.0 South Napa, California, earthquake rupture, each dominated by either co- or post-seismic shallow fault slip. We combine mobile laser scanning (MLS), active-source seismic tomography, and finite element modeling to investigate how deformation rate and mechanical properties of the shallow crust affect fault behavior. Despite four orders-of-magnitude difference in the rupture velocities, MLS-derived shear strain fields are remarkably similar at the two sites and suggest deceleration of the co-seismic rupture near Earth’s surface. Constrained by the MLS and seismic data, finite element models indicate shallow faulting is more sensitive to lithologic layering and plastic yielding than to the presence of fault compliant zones (i.e., regions surrounding faults with reduced stiffness). Although both elastic and elastoplastic models can reproduce the observed surface displacement fields within the uncertainty of MLS data, elastoplastic models likely provide the most reliable representations of subsurface fault behavior, as they produce geologically reasonable stress states and are consistent with field, geodetic, and seismological observations.

California↗

USGS1994

Late Quaternary normal faulting of the Hat Creek basalt, northern California

The Hat Creek fault is a major, young, north-striking, normal fault along the western boundary of extensional Basin and Range deformation in the Lessen region of northeastern California. Volcanic rocks of Quaternary and late Pliocene age are displaced a total of >500 m down to the west along west-facing, en echelon scarps now retreated to ∼35° slopes. Fresh, young scarps as much as 30 m high cut the Hat Creek Basalt (erupted between 15 and ∼40 ka) a few tens of meters west of the retreated scarps. Prior to the late 1980s, these young scarps were interpreted as lava slump scarps formed as the Hat Creek Basalt ponded against the older fault scarps and then drained away to the northwest. Numerous pieces of geologic evidence, however, show that the young scarps formed after the Hat Creek Basalt solidified and cooled and are true fault features formed by the youngest displacements of the Hat Creek fault. Structural details are remarkably well preserved along the series of left-stepping scarps cutting the Hat Creek Basalt. Near the central parts of individual segments, the fault is displayed as a single, vertical scarp. Near the ends of the segments, the scarp decreases in height and becomes a monodinal flexure on which the recent dip separation has been taken up by small-scale offset along columnar cooling joints in the basalt These monodinal flexures commonly rotate into east-west monoclines that join adjacent north-south segments. Displacement of outwash gravel overlying the Hat Creek Basalt shows that vertical separation on the Hat Creek fault has averaged ∼1.3 mm yr −1 for the past 15,000 yr. The Hat Creek fault thus represents a potential earthquake hazard, despite the low level and diffuse nature of modern seismidty in the region.

Geological Society of America Bulletin↗

USGS1992

Lithofacies analysis of colluvial sediments - an aid in interpreting the recent history of Quaternary normal faults in the Basin and Range Province, western United States

Inferring the frequency and magnitude of past earthquakes from the stratigraphy in exposures of normal-faulted sediments is difficult because colluvial lithofacies assemblages adjacent to faults are complex. Similarities in facies assemblages adjacent to young fault scarps in arid to semiarid areas, such as the Basin and Range province, allow lithofacies to be grouped into two genetic architectural elements: debris and wash elements. Upper and lower facies associations can commonly be recognized within each element. A lithofacies code scheme, similar to those used in the analysis of fluvial and glacial lithofacies sequences, provides a concise way of illustrating lithofacies relations in fault exposures. The source lithology of colluvial lithofacies is shown in the code, and soil-horizon symbols can be included. The architecture of lithofacies assemblages near fault scarps in semiarid areas is explained by a model of colluvial sedimentation in response to a single surface faulting event. Analysis of lithofacies assemblages exposed in three trenches across normal faults in the eastern Basin and Range shows how the model can be used to interpret fault histories. Similar facies analysis methods may be useful in interpreting colluvial sequences formed by non-tectonic processes.

Journal of Sedimentary Petrology↗

USGS2005

Teleseismic body waves from dynamically rupturing shallow thrust faults: Are they opaque for surface-reflected phases?

We investigate whether a shallow-dipping thrust fault is prone to waveslip interactions via surface-reflected waves affecting the dynamic slip. If so, can these interactions create faults that are opaque to radiated energy? Furthermore, in this case of a shallow-dipping thrust fault, can incorrectly assuming a transparent fault while using dislocation theory lead to underestimates of seismic moment? Slip time histories are generated in three-dimensional dynamic rupture simulations while allowing for varying degrees of wave-slip interaction controlled by fault-friction models. Based on the slip time histories, P and SH seismograms are calculated for stations at teleseismic distances. The overburdening pressure caused by gravity eliminates mode I opening except at the tip of the fault near the surface; hence, mode I opening has no effect on the teleseismic signal. Normalizing by a Haskell-like traditional kinematic rupture, we find teleseismic peak-to-peak displacement amplitudes are approximately 1.0 for both P and SH waves, except for the unrealistic case of zero sliding friction. Zero sliding friction has peak-to-peak amplitudes of 1.6 for P and 2.0 for SH waves; the fault slip oscillates about its equilibrium value, resulting in a large nonzero (0.08 Hz) spectral peak not seen in other ruptures. These results indicate wave-slip interactions associated with surface-reflected phases in real earthquakes should have little to no effect on teleseismic motions. Thus, Haskell-like kinematic dislocation theory (transparent fault conditions) can be safety used to simulate teleseismic waveforms in the Earth.

Bulletin of the Seismological Society of America↗

USGS2015

Long Valley caldera and the UCERF depiction of Sierra Nevada range-front faults

Long Valley caldera lies within a left-stepping offset in the north-northwest-striking Sierra Nevada range-front normal faults with the Hilton Creek fault to the south and Hartley Springs fault to the north. Both Uniform California Earthquake Rupture Forecast (UCERF) 2 and its update, UCERF3, depict slip on these major range-front normal faults as extending well into the caldera, with significant normal slip on overlapping, subparallel segments separated by &sim;10  km. This depiction is countered by (1) geologic evidence that normal faulting within the caldera consists of a series of graben structures associated with postcaldera magmatism (intrusion and tumescence) and not systematic down-to-the-east displacements consistent with distributed range-front faulting and (2) the lack of kinematic evidence for an evolving, postcaldera relay ramp structure between overlapping strands of the two range-front normal faults. The modifications to the UCERF depiction described here reduce the predicted shaking intensity within the caldera, and they are in accord with the tectonic influence that underlapped offset range-front faults have on seismicity patterns within the caldera associated with ongoing volcanic unrest.

Bulletin of the Seismological Society of America↗

USGS2026

Simulations of earthquakes in California with multiple velocity models

Estimating seismic hazard using 3D earthquake simulations of moderate‐to‐large magnitude events is becoming a reality. In this study, we further this endeavor by simulating seismic wavefields for five earthquakes in California, United States, and evaluating their pseudospectral acceleration at periods of 2, 5, and 10 s. We compare results from the recently developed U.S. Geological Survey (USGS) National Crustal Model with three well‐established and validated regional community velocity models (CVMs): two CVMs for Southern California (CVMS4.26.M01 and CVMH15.1.0) and the USGS San Francisco Bay region 3D CVM. A minimum shear‐wave velocity of 500 m/s is imposed in all four of the velocity models. We use finite‐difference calculations to simulate the earthquakes and specify a grid and timestep capable of resolving periods down to 1 s. For our simulations, we find that for each earthquake, the available velocity models produce broadly similar amplitudes and patterns of earthquake ground motions. Use of the 3D velocity models generally perform better than ergodic ground‐motion models, and variability in the results between earthquakes may indicate a combination of factors including uncertainty in earthquake magnitude and source parameterization, signal‐to‐noise ratio in the observed time series, and the lack of fault damage zones, topography, and other features in the computational seismic velocity models. Of particular note, we find that for most cases, averaging the ground motions for all available velocity models outperforms the ground‐motion results from any single velocity model. Our results indicate that the 3D models considered in this study can be sufficiently independent and have enough predictive power that their combination could be used to improve estimates of earthquake ground motions.

California↗

USGS2026

Intensity prediction equations for intraslab earthquakes in the Cascadia Subduction Zone

In the U.S. Pacific Northwest (PNW), intraslab earthquakes located within the subducting Juan de Fuca slab are responsible for most of the region’s moderate‐magnitude seismicity. Historically, these earthquakes have caused significant damage, such as the 2001 M 6.8 Nisqually, 1965 M 6.5 Seattle, and 1949 M 6.8 Olympia events. Intensity prediction equations (IPEs) aim to relate source‐to‐site distance and earthquake magnitude to the resultant felt shaking intensity at a particular location, and help to estimate hazard from earthquakes. We find that a commonly used IPE developed by Atkinson et al. (2014) does not accurately describe shaking intensity from most intraslab events in the PNW, likely because it was developed using primarily crustal fault earthquakes. In this study, we develop two new IPEs that provide an improved fit to observed intensities from Cascadia subduction zone intraslab earthquakes. To develop our IPEs, we compute linear mixed‐effects regressions using community decimal intensity data from the U.S. Geological Survey’s ‘Did You Feel It?’ program from 36 intraslab events that occurred in British Columbia, Canada, Washington, Oregon, and northern California from 2000 to 2025. We also investigate the relationship between the new intraslab IPEs and select parameters related to near‐surface conditions and deeper geologic structure that may be related to ground‐motion intensities, but find no correlations. These new IPEs provide a better fit to observations of intraslab earthquakes in the PNW than Atkinson et al. (2014) , particularly within 300 km of the epicenter. They also provide a better fit to the 2018 M 7.1 Anchorage, Alaska, intraslab earthquake, suggesting that they may have broader applicability to global subduction zones.

California, Oregon, Washington↗