Geology topics

Research about Cascadia subduction zone

Source-linked reports with geographic coverage including Cascadia subduction zone.

17 recordsLinked to original sources

USGS2026

Widespread abyssal turbidites record megathrust earthquake-triggered landslides and coseismic deformation in the Cascadia subduction zone

Abyssal marine turbidites provide some of the longest and most spatially extensive records of subduction zone earthquake recurrence globally; however, correlation of these deposits over long distances and interpretation of synchronous emplacement requires both an understanding of the turbidite generating systems and precise dating. Here, we present an integrated suite of high-resolution bathymetry, subbottom profiles, and sediment cores from combined autonomous underwater vehicle, remotely operated vehicle, and ship-based studies at a key paleoseismic site in the southern Cascadia subduction zone. We demonstrate how widespread, earthquake-triggered landslides on the lower slope deposit discrete, proximal mass transport deposits (MTDs) that grade offshore into complex, interfingered abyssal turbidites, which correspond to records of megathrust earthquake history. We propose accretion and oversteepening of thrust folds on the lower slope both preconditions the slope to fail and provides a perpetual source of unstable material to fail during every earthquake cycle. Furthermore, we suggest the periodic and pervasive landsliding indicates coseismic deformation of the outer accretionary wedge during megathrust rupture.

California, Oregon↗

USGS2025

Structural controls on splay fault rupture dynamics during Cascadia megathrust earthquakes

Great subduction earthquakes ( M w ≥ 8.0) can generate devastating tsunamis by rapidly displacing the seafloor and overlying water column. These potentially tsunamigenic seafloor offsets result from coseismic fault slip and deformation beneath or within the accretionary wedge. The mechanics of these shallow rupture phenomena and their dependence on subduction zone properties remain unresolved, partly due to the sparsity of offshore observations of shallow megathrust earthquake deformation. Here, we analyze how offshore structure influences shallow rupture mechanics and slip partitioning using 3D dynamic earthquake simulations of the Cascadia subduction zone (CSZ) megathrust with and without variably dipping seaward- or landward-vergent splay faults in the wedge that sole into the megathrust. Resulting tradeoffs between splay and megathrust slip reveal structural controls on rupture partitioning, with greater splay slip leading to less shallow megathrust slip updip. Gently dipping and seaward-vergent splays host more slip than those with steeper, landward-vergent splays. To isolate the underlying mechanisms, we compare models with Andersonian and plunging principal stresses. Results suggest distinct static and dynamic processes control the dip- and vergence-dependence of splay rupture: static (mis)alignment relative to far-field tectonic loading favors slip on more optimally oriented, shallowly dipping splay faults. In contrast, dynamic stress interactions of an updip-propagating megathrust rupture front with the free surface and potential branch faults favor forward branching onto seaward-vergent splays and inhibit backward branching onto landward-vergent splays. Resulting seafloor displacements suggest splay fault structure may influence coseismic tsunami source processes, highlighting the importance of dynamically viable rupture scenarios in subduction hazard assessments.

Cascadia subduction zone↗

USGS2025

Using subducting plate motion to constrain Cascadia slab geometry and interface strength

Subduction zones are home to multiple geohazards driven by the evolution of the regional tectonics, including earthquakes, volcanic eruptions and landslides. Past evolution builds the present-day structure of the margin, while the present-day configuration of the system determines the state-of-stress in which individual hazardous events manifest. Regional simulations of subduction zones provide a tool to synthesize the tectonic history of a region and investigate how geologic features lead to variations in the state of stress across the subduction system. However, it is challenging to design regional models that provide a force-balance that is consistent with the large-scale motion of surrounding tectonic plates while also not over-constraining the solution. Here, we present new models for the Cascadia subduction zone that meet these criteria and demonstrate how the motion of the subducting Juan de Fuca plate can be used to determine the along-strike variations in the viscous (long-term) coupling across the plate boundary. All successful models require lower viscous coupling in the northern section of the trench compared to the central and southern sections. However, due to uncertainties in the geometry of the Cascadia slab, we find that there is a trade-off between along-strike variation in viscous coupling and slab shape. Better constraints on the slab shape, and/or use of other observations are needed to resolve this trade-off. The approach presented here provides a framework for further exploring how geologic features in the overriding plate and the properties of the plate boundary region affect the state-of-stress across this and other subduction zones.

Cascadia subduction zone↗

USGS2025

No evidence for an active margin-spanning megasplay fault at the Cascadia Subduction Zone

It has been previously proposed that a megasplay fault within the Cascadia accretionary wedge, spanning from offshore Vancouver Island to Oregon, has the potential to slip during a future Cascadia subduction zone earthquake. This hypothetical fault has major implications for tsunami size and arrival times and is included in disaster-planning scenarios currently in use in the region. This hypothesis is evaluated in this study using CASIE21 deep-penetrating and U.S. Geological Survey high-resolution seismic reflection profiles. We map changes in wedge structural style and seismic character to identify the inner-outer wedge transition zone where a megasplay fault has been previously hypothesized to exist and evaluate evidence for active faulting within this zone. Our results indicate that there is not an active, through-going megasplay fault in Cascadia, but instead, the structure and activity of faulting at the inner-outer wedge transition zone is highly variable and segmented along strike, consistent with the segmentation of other physical and mechanical properties in Cascadia. Wedge sedimentation, plate dip, and subducting topography are proposed to play a major role in controlling megasplay fault development and evolution. Incorporating updated megasplay fault location, geometry, and activity into modeling of Cascadia earthquakes and tsunamis could help better constrain associated hazards.

British Columbia, Oregon, Washington↗

USGS2024

Evidence for low effective stress within the crust of the subducted Gorda plate from the 2022 December Mw 6.4 Ferndale earthquake sequence

Stress levels on and adjacent to megathrust faults at seismogenic depths remain a key but difficult to constrain parameter for assessing seismic hazard in subduction zones. Although strong ground motions have been observed to be generated from distinct, high-stress regions on the downdip end of the megathrust rupture areas in many great earthquakes, we lack direct constraints on the stress level in the lower seismogenic portion of the Cascadia megathrust. On 2022 December 20, a M w 6.4 strike-slip earthquake occurred near Ferndale, California in southern Cascadia and likely ruptured the Gorda slab crust in the lower seismogenic portion, providing an opportunity to assess the stress level in this region. Here, we relocate the Ferndale mainshock and the first two weeks of aftershocks using a high-resolution 3-D velocity model and estimate rupture dimensions, directivity, and stress drop for several M w 4-5 aftershocks and recent earthquakes. The aftershocks define a strike-slip fault in the slab crust striking ENE, consistent with the mainshock focal mechanism. The orientation of this fault is about 45 ° off the ideally oriented fault plane given the stress state in the slab. The aftershock zone is extensive and broad in the forward direction of the mainshock rupture but still constrained within the volume of high Vp/Vs within the slab crust. Our stress drop estimates are generally lower for M w 4-5 earthquakes located within the slab crust compared to those a few km deeper in the slab mantle. Combined, our results support a relatively low effective stress level in the vicinity of the megathrust in the lower portion of the seismogenic zone in southern Cascadia, likely due to elevated fluid pressures. Consequently, the ground motion in the onshore region above this low-stress seismogenic portion in southern Cascadia may not be as intense as that observed during great earthquakes in other subduction zones.

California↗

USGS2022

Mass wasting along the Cascadia subduction zone: Implications for abyssal turbidite sources and the earthquake record

The only submarine records of large (>M w 7) prehistoric earthquakes along the Cascadia subduction zone are derived from sequences of deep sea turbidites interpreted to represent synchronous, shaking-induced failures along the continental slope. However, the spatial correlation of these deposits along the margin is complicated and the chronological constraints involve significant uncertainties, raising questions about how these deposits were generated. Here we present the most comprehensive spatial database of seafloor failure scarps across the Cascadia margin to date, with more than 8700 features mapped. We observe pervasive mass wasting of the steep lower slope along the ∼800 km length of the margin. This portion of the margin is heavily deformed, uplifted and oversteepened, such that it may be optimally preconditioned to fail during intense shaking events. Our results suggest disintegration of the steep lower slope is likely the primary source of turbidity flows triggered by earthquake shaking in Cascadia. This more proximal source for seismoturbidite generation may be the principal process along subduction margins with steep lower slope topography and may explain why seismoturbidites are found in areas containing few submarine canyons. Furthermore, the extensive mass wasting present on the steep lower slope of central Cascadia, a critical 250 km long section where no records currently exist, suggests additional turbidite records may be found along this portion of the margin by selecting optimal sample sites based on failure patterns. The strong clustering of seafloor failures along the lower slope also suggests the most earthquake intense shaking may be spatially restricted to a narrow zone along the deformation front, which could have important implications for hazard predictions across the region.

California, Oregon, Washington↗

USGS2021

Cataloging tectonic tremor energy radiation in the Cascadia subduction zone

For the past ∼12 years the Pacific Northwest Seismic Network has been automatically detecting and locating tectonic tremor across the Cascadia subduction zone, resulting in a catalog of more than 500,000 tremor epicenters to date, which has served as a valuable resource for tremor and slip research. This manuscript presents an updated methodology for routine tremor detection in Cascadia and a new catalog of over 180,000 tremor epicenters including amplitudes detected along the subduction zone margin from 2017 to 2021. The events are detected via cross-correlation of continuous vertical envelope data of 128 stations from northern California to northern Vancouver Island. The modified approach results in less scatter and a 55% increase in detected epicenters than previously observed, as well as a newly identified tremor source offset updip from the main tremor and slip region at the southern edge of the subduction zone. Radiated seismic energy in the 1.5–5 Hz band is used to assign epicenters an energy magnitude ( M eL ), which is calibrated to the M L of local earthquakes. Southern Cascadia is most active, but the highest tremor energy rates occur in northern Cascadia. Tremor in central Cascadia is systematically weaker and less frequent. Individual epicenter magnitudes range from ∼0.5–2 and spatiotemporally cluster into 1,060 swarms with cumulative M eL ranging from ∼0.8 to 3.7. The swarms reflect underlying slow slip events and occur with an earthquake-like energy distribution with a b value ∼1. Tremor epicenters, however, follow a tapered Gutenberg-Richter distribution with high b values, suggesting individual tremor bursts and their constituent low-frequency earthquakes are fault-dimension limited.

California, Oregon, Washington↗

USGS2021

Imaging the next Cascadia earthquake: Optimal design for a seafloor GNSS- A network

The Cascadia subduction zone in the Pacific Northwest of the United States of America capable of producing magnitude ∼9 earthquakes, likely often accompanied by tsunamis. An outstanding question in this region is the degree and spatial extent of interseismic strain accumulation on the subduction megathrust. Seafloor geodetic methods combining GNSS and underwater acoustic ranging (GNSS-A) are capable of imaging this strain accumulation on the offshore portion of the subduction zone and therefore anticipating the potential size and rupture pattern of a future earthquake. However, the high cost of seafloor geodesy means that only a limited number of stations may be deployed and monitored. To facilitate expansion of current geodetic networks offshore, we develop a quantitative recommendation of optimal locations for future seafloor geodetic observations, based on the amount of new information provided by that observation. The optimal network depends on the problem that one is trying to solve with those observations (mapping subduction locking rates, coupling rates, constraining total moment rate, etc.), and on a number of modelling and data uncertainty assumptions. In particular, data uncertainty assumptions will change over time, as more position observations reduce velocity uncertainties. We find that near-trench observations on the megathrust hangingwall, distributed along-strike, consistently provide significant reduction in differential entropy over a large suite of assumptions, and that a well-placed seafloor observation can provide up to ∼30 times the information gain of the most optimal onshore observation.

British Columbia, California, Oregon, Washington↗

USGS2021

Commentary: The role of geodetic algorithms for earthquake early warning in Cascadia

The ShakeAlert earthquake early warning (EEW) system issues public alerts in California and will soon extend to Oregon and Washington. The Cascadia subduction zone presents significant new challenges and opportunities for EEW. Initial publications suggested that EEW algorithms based on Global Navigation Satellite System (GNSS) data could provide improved warning for intraslab events and dramatically improved warning for offshore megathrust events, both of which contribute significantly to hazard in Cascadia. We find that some expectations in these publications were unrealistic, and we demonstrate that in general geodetic algorithms would not produce timely warnings for intraslab events nor warning times of two minutes or more for severe shaking from megathrust earthquakes. Nonetheless, lessons from recent earthquakes in Japan and California, for which alerts from seismic algorithms suffered from magnitude saturation and high data latencies, demonstrate the urgent need for rigorous testing of geodetic EEW as a potential complement to seismic EEW.

California, Oregon, Washington↗

USGS2020

A long-term geothermal observatory across subseafloor gas hydrates, IODP Hole U1364A, Cascadia accretionary prism

We report 4 years of temperature profiles collected from May 2014 to May 2018 in Integrated Ocean Drilling Program Hole U1364A in the frontal accretionary prism of the Cascadia subduction zone. The temperature data extend to depths of nearly 300 m below seafloor (mbsf), spanning the gas hydrate stability zone at the location and a clear bottom-simulating reflector (BSR) at ∼230 mbsf. When the hole was drilled in 2010, a pressure-monitoring Advanced CORK (ACORK) observatory was installed, sealed at the bottom by a bridge plug and cement below 302 mbsf. In May 2014, a temperature profile was collected by lowering a probe down the hole from the ROV ROPOS. From July 2016 through May 2018, temperature data were collected during a nearly two-year deployment of a 24-thermistor cable installed to 268 m below seafloor (mbsf). The cable and a seismic-tilt instrument package also deployed in 2016 were connected to the Ocean Networks Canada (ONC) NEPTUNE cabled observatory in June of 2017, after which the thermistor temperatures were logged by Ocean Networks Canada at one-minute intervals until failure of the main ethernet switch in the integrated seafloor control unit in May 2018. The thermistor array had been designed with concentrated vertical spacing around the bottom-simulating reflector and two pressure-monitoring screens at 203 and 244 mbsf, with wider thermistor spacing elsewhere to document the geothermal state up to seafloor. The 4 years of data show a generally linear temperature gradient of 0.055°C/m consistent with a heat flux of 61–64 mW/m 2 . The data show no indications of thermal transients. A slight departure from a linear gradient provides an approximate limit of ∼10 −10 m/s for any possible slow upward advection of pore fluids. In-situ temperatures are ∼15.8°C at the BSR position, consistent with methane hydrate stability at that depth and pressure.

Cascadia subduction zone↗

USGS2020

Systematic characterization of morphotectonic variability along the Cascadia convergent margin: Implications for shallow megathrust behavior and tsunami hazards

Studies of recent destructive megathrust earthquakes and tsunamis along subduction margins in Japan, Sumatra, and Chile have linked forearc morphology and structure to megathrust behavior. This connection is based on the idea that spatial variations in the frictional behavior of the megathrust influence the tectono-morphological evolution of the upper plate. Here we present a comprehensive examination of the tectonic geomorphology, outer wedge taper, and structural vergence along the marine forearc of the Cascadia subduction zone (offshore northwestern North America). The goal is to better understand geologic controls on outer wedge strength and segmentation at spatial scales equivalent to rupture lengths of large earthquakes (≥M 6.7), and to examine potential linkages with shallow megathrust behavior. We use cross-margin profiles, spaced 25 km apart, to characterize along-strike variation in outer wedge width, steepness, and structural vergence (measured between the toe and the outer arc high). The width of the outer wedge varies between 17 and 93 km, and the steepness ranges from 0.9° to 6.5°. Hierarchical cluster analysis of outer wedge width and steepness reveals four distinct regions that also display unique patterns of structural vergence and shape of the wedge: Vancouver Island, British Columbia, Canada (average width, linear wedge, seaward and mixed vergence); Washington, USA (higher width, concave wedge, landward and mixed vergence); northern and central Oregon, USA (average width, linear and convex wedge, mixed and seaward vergence); and southern Oregon and northern California, USA (lower width, convex wedge, seaward and mixed vergence). Variability in outer wedge morphology and structure is broadly associated with along-strike megathrust segmentation inferred from differences in oceanic asthenospheric velocities, patterns of episodic tremor and slow slip, GPS models of plate locking, and the distribution of seismicity near the plate interface. In more detail, our results appear to delineate the extent, geometry, and lithology of dynamic and static backstops along the margin. Varying backstop configurations along the Cascadia margin are interpreted to represent material-strength contrasts within the wedge that appear to regulate the along- and across-strike taper and structural vergence in the outer wedge. We argue that the morphotectonic variability in the outer wedge may reflect spatial variations in shallow megathrust behavior occurring over roughly the last few million years. Comparing outer wedge taper along the Cascadia margin to a global compilation suggests that observations in the global catalog are not accurately representing the range of heterogeneity within individual margins and highlights the need for detailed margin-wide morphotectonic analyses of subduction zones worldwide.

British Columbia, California, Oregon, Washington↗

USGS2020

Identifying the greatest earthquakes of the past 2000 years at the Nehalem River Estuary, Northern Oregon Coast, USA

We infer a history of three great megathrust earthquakes during the past 2000 years at the Nehalem River estuary based on the lateral extent of sharp (≤3 mm) peat-mud stratigraphic contacts in cores and outcrops, coseismic subsidence as interpreted from fossil diatom assemblages and reconstructed with foraminiferal assemblages using a Bayesian transfer function, and regional correlation of 14 C-modeled ages for the times of subsidence. A subsidence contact from 1700 CE (contact A), sometimes overlain by tsunami-deposited sand, can be traced over distances of 7 km. Contacts B and D, which record subsidence during two earlier megathrust earthquakes, are much less extensive but are traced across a 700-m by 270-m tidal marsh. Although some other Cascadia studies report evidence for an earthquake between contacts B and D, our lack of extensive evidence for such an earthquake may result from the complexities of preserving identifiable evidence of it in the rapidly shifting shoreline environments of the lower river and bay. Ages (95% intervals) and subsidence for contacts are: A, 1700 CE (1.1 ± 0.5 m); B, 942–764 cal a BP (0.7 ± 0.4 m and 1.0 m ± 0.4 m); and D, 1568–1361 cal a BP (1.0 m ± 0.4 m). Comparisons of contact subsidence and the degree of overlap of their modeled ages with ages for other Cascadia sites are consistent with megathrust ruptures many hundreds of kilometers long. But these data cannot conclusively distinguish among different types or lengths of ruptures recorded by the three great earthquake contacts at the Nehalem River estuary.

California, Oregon, Washington, British Columbia↗

USGS2017

Sediment gravity flows triggered by remotely generated earthquake waves

Recent great earthquakes and tsunamis around the world have heightened awareness of the inevitability of similar events occurring within the Cascadia Subduction Zone of the Pacific Northwest. We analyzed seafloor temperature, pressure, and seismic signals, and video stills of sediment-enveloped instruments recorded during the 2011–2015 Cascadia Initiative experiment, and seafloor morphology. Our results led us to suggest that thick accretionary prism sediments amplified and extended seismic wave durations from the 11 April 2012 M w 8.6 Indian Ocean earthquake, located more than 13,500 km away. These waves triggered a sequence of small slope failures on the Cascadia margin that led to sediment gravity flows culminating in turbidity currents. Previous studies have related the triggering of sediment-laden gravity flows and turbidite deposition to local earthquakes, but this is the first study in which the originating seismic event is extremely distant (> 10,000 km). The possibility of remotely triggered slope failures that generate sediment-laden gravity flows should be considered in inferences of recurrence intervals of past great Cascadia earthquakes from turbidite sequences. Future similar studies may provide new understanding of submarine slope failures and turbidity currents and the hazards they pose to seafloor infrastructure and tsunami generation in regions both with and without local earthquakes.

Cascadia subduction zone↗

USGS2014

Rethinking turbidite paleoseismology along the Cascadia subduction zone

A stratigraphic synthesis of dozens of deep-sea cores, most of them overlooked in recent decades, provides new insights into deep-sea turbidites as guides to earthquake and tsunami hazards along the Cascadia subduction zone, which extends 1100 km along the Pacific coast of North America. The synthesis shows greater variability in Holocene stratigraphy and facies off the Washington coast than was recognized a quarter century ago in a confluence test for seismic triggering of sediment gravity flows. That test compared counts of Holocene turbidites upstream and downstream of a deep-sea channel junction. Similarity in the turbidite counts among seven core sites provided evidence that turbidity currents from different submarine canyons usually reached the junction around the same time, as expected of widespread seismic triggering. The fuller synthesis, however, shows distinct differences between tributaries, and these differences suggest sediment routing for which the confluence test was not designed. The synthesis also bears on recent estimates of Cascadia earthquake magnitudes and recurrence intervals. The magnitude estimates hinge on stratigraphic correlations that discount variability in turbidite facies. The recurrence estimates require turbidites to represent megathrust earthquakes more dependably than they do along a flow path where turbidite frequency appears limited less by seismic shaking than by sediment supply. These concerns underscore the complexity of extracting earthquake history from deep-sea turbidites at Cascadia.

Cascadia subduction zone↗

USGS2001

Geophysics: Silent slip on the Cascadia subduction interface

G eodetic satellites of the Global Positioning System (GPS) now permit continuous recording of surface motions around earthquake faults and volcanoes with millimeter precision. Data from distinct points on Earth's surface can be combined to infer the locations of the sources of deformation at depth in the crust. These data are radically altering our understanding of earthquake processes, justifying ambitious new sensor arrays to image active deformation sources in Earth's crust. In recent years, large continuous GPS arrays have been deployed in Japan and southern California, and smaller networks have been installed in other seismically active regions. On page 1525 of this issue ( 1 ), Dragert et al. provide a glimpse of the kinds of insights we can expect from these arrays.

Cascadia subduction zone↗

USGS2000

Models of downdip frictional coupling for the Cascadia Megathrust

We have developed models of downdip frictional coupling along two transects across the Cascadia subduction zone in northern Washington and central Oregon. The models involve complicated downdip coupling profiles. Although not unique, our models closely predict available GPS shortening rates and vertical uplift data. We are able to explain relatively low horizontal shortening rates along the Washington coast and small vertical uplift rates in central Oregon. Our models depart from previous models by inclusion of a deeply coupled region assumed to be related to mafic rocks in both the upper and lower plates of the thrust.

Oregon, Washington↗

USGS1998

New seismic images of the cascadia subduction zone from cruise SO 108-ORWELL

In April and May 1996, a geophysical study of the Cascadia continental margin off Oregon and Washington was conducted aboard the German R/V Sonne. This cooperative experiment by GEOMAR and the USGS acquired wide-angle reflection and refraction seismic data, using ocean-bottom seismometers (OBS) and hydrophones (OBH), and multichannel seismic reflection (MCS) data. The main goal of this experiment was to investigate the internal structure and associated earthquake hazard of the Cascadia subduction zone and to image the downgoing plate. Coincident MCS and wide-angle profiles along two tracks are presented here. The plate boundary has been imaged precisely beneath the wide accretionary wedge close to shore at c13km depth. Thus, the downgoing plate dips more shallowly than previously assumed. The dip of the plate changes from 2?? to 4?? at the eastern boundary of the wedge on the northern profile, whereas approximately 3km of sediment is entering the subduction zone. On the southern profile, where the incoming sedimentary section is about 2.2km thick, the plate dips about 0.5?? to 1.5?? near the deformation front and increases to 3.5?? further landwards. On both profiles, the deformation of the accretionary wedge has produced six ridges on the seafloor, three of which represent active faulting, as indicated by growth folding. The ridges are bordered by landward verging faults which reach as deep as the top of the oceanic basement. Thus, the entire incoming sediment package is being accreted. At least two phases of accretion are evident, and the rocks of the older accretionary phase(s) forms the backstop for the younger phase, which started around 1.5 Ma ago. This documents that the 30 to 50km wide frontal part of the accretionary wedge, which is characterized by landward vergent thrusts, is a Pleistocene feature which was formed in response to the high input of sediment building the fans during glacial periods. Velocities increase quite rapidly within the wedge, both landward and downward. At the toe of the deformation front, velocities are higher than 4.0 km/s, indicating extensive dewatering of deep, oceanic sediment. Further landward, considerable velocity variation is found, which indicates major breaks throughout the accretionary history.

British Columbia, California, Idaho, Nevada, Orego↗