skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on July 4, 2026

Title: Effects of incoming polygonal fault systems on subduction zone and slow slip behavior
The physical properties of subduction inputs profoundly influence megathrust slip behavior. Seismic data reveal extensive polygonal fault systems (PFSs) in the input sequences of the Hikurangi Margin and Nankai Trough. The mechanical and hydrological effects of these incoming PFSs on subduction zones are potentially substantial. Here, we investigate their effects following transport into the accretionary wedge by integrating discrete-element modeling with three-dimensional seismic interpretation. We find that the typical dips of the incoming PFSs overlap with modeled dips prone to reactivation and confirm that subducting PFSs can be reactivated and gradually evolve into major thrust faults. Comparisons with electromagnetic data indicate that PFSs may provide conduits for fluid leakage along the plate interface, coincide with disrupted strata and decreased shear stress, and enhance geometric and stress heterogeneity along the megathrust. These suggest that PFSs may play a previously unrecognized role in contributing to shallow slow earthquake phenomena in subduction zones.  more » « less
Award ID(s):
2123255
PAR ID:
10653686
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  
Publisher / Repository:
AAAS
Date Published:
Journal Name:
Science Advances
Volume:
11
Issue:
27
ISSN:
2375-2548
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Oceanic plate seamounts are believed to play an important role in megathrust rupture at subduction zones, although consistent relationships between subducting seamounts and plate interface seismicity patterns are not found. While most studies focus on impacts linked to their topography, seamounts are also sites of heterogeneity in incoming plate sediments that may contribute to megathrust properties. Here, we characterize incoming plate sediments along the Cascadia subduction zone using new high‐resolution seismic images and compressional wave (Vp) models from the CASIE21 multi‐channel‐seismic experiment. Nine fully‐to‐partially buried seamounts are identified seaward of the deformation front within a region of thick Plio‐Pleistocene sediment where the Juan de Fuca plate is bending into the subduction zone. Anomalously highVpsediment blankets two seamounts offshore Washington‐Central Oregon, with wavespeeds reaching 36% and 20% higher than adjacent sediment. Fluid seepage and temperatures warm enough for smectite diagenesis extending to shallow depths are inferred from heat flow studies and we attributeVpanomalies to sediment cementation linked primarily to smectite dehydration. Signatures of fluid seepage above seamounts are also identified offshore Vancouver Island, but anomalously lowVpsediment below distinct reverse polarity reflections are found, indicating trapped fluids, and cooler basement temperatures are inferred. Landward of one seamount, a zone of enhanced sediment compaction is found, consistent with the predicted stress modulating effects of seamount subduction. These new findings of variations in sediment diagenesis and strength around seamounts prior to subduction may contribute to the diverse megathrust frictional properties and seismicity patterns evident at subducting seamounts. 
    more » « less
  2. To explore controls on megathrust behavior and its connection with forearc deformation, we studied the Andreanof segment of the Aleutian Subduction Zone (offshore Alaska, USA), which has a simple geological history as a relatively young intra-oceanic subduction zone. Here, the forearc shows greater uplift and compression in the strongly coupled Adak region compared to the weakly coupled Atka region. Using multichannel seismic reflection data, we found that the incoming plate in both regions exhibits similar characteristics along the segment, suggesting that its properties do not account for the varying megathrust behavior and forearc deformation. Instead, differences between the Atka and Adak regions in the thickness of the methane hydrate stability zone, as marked by a bottom-simulating reflector, suggest more heat advection, and thus dewatering, in the Adak region, where the more developed fault network may enable fluid drainage, thereby lowering pore pressure at the megathrust and promoting coupling. Higher coupling allows for seismic and stress cycling that would sustain forearc permeability by faulting. Our results suggest a feedback between deformation and coupling that may be active or latent in other more complex subduction zones but in concert with or masked by other factors. 
    more » « less
  3. In subduction zones, along‐strike and downdip variations in megathrust slip behavior are linked to changes in the properties of the subducting and overriding plates. Although marine geophysical methods provide insights into subduction zone structures, most surveys consist of sparse 2D profiles, limiting our understanding of first‐order controls. Here, we use active‐source seismic data to derive a 3D crustal‐scale P‐wave velocity model of the Alaska Peninsula subduction zone that encompasses both plates and spans the Semidi segment and SW Kodiak asperity. Our results reveal modest variations within the incoming plate, attributed to a series of fracture zones, seamounts and their associated basement swell, collectively contributing to plate hydration. Basement swell appears to modulate the distribution and type of sediment entering the trench, likely impacting observed variations in slip behavior. The overriding plate exhibits significant heterogeneity. The updip limit and width of the dynamic backstop are similar between the SW Kodiak asperity and eastern Semidi segment but differ significantly from the Western Semidi segment. These distinctions likely account for differences in earthquake rupture patterns and interseismic coupling among these segments. Additionally, high velocities in the mid‐lower forearc crust coincide with the location of the megathrust slip during the Mw 8.2 2021 Chignik event. We interpret these velocities as intracrustal intrusions that contributed to the deep rupture of the 2021 event. Our findings suggest that the contrasting structural and material properties of both the incoming and overriding plates influence the spatially complex and semi‐persistent segmentation of the megathrust offshore the Alaska Peninsula. 
    more » « less
  4. Abstract From California to British Columbia, the Pacific Northwest coast bears an omnipresent earthquake and tsunami hazard from the Cascadia subduction zone. Multiple lines of evidence suggests that magnitude eight and greater megathrust earthquakes have occurred ‐ the most recent being 321 years ago (i.e., 1700 A.D.). Outstanding questions for the next great megathrust event include where it will initiate, what conditions are favorable for rupture to span the convergent margin, and how much slip may be expected. We develop the first 3‐D fully dynamic rupture simulations for the Cascadia subduction zone that are driven by fault stress, strength and friction to address these questions. The initial dynamic stress drop distribution in our simulations is constrained by geodetic coupling models, with segment locations taken from geologic analyses. We document the sensitivity of nucleation location and stress drop to the final seismic moment and coseismic subsidence amplitudes. We find that the final earthquake size strongly depends on the amount of slip deficit in the central Cascadia region, which is inferred to be creeping interseismically, for a given initiation location in southern or northern Cascadia. Several simulations are also presented here that can closely approximate recorded coastal subsidence from the 1700 A.D. event without invoking localized high‐stress asperities along the down‐dip locked region of the megathrust. These results can be used to inform earthquake and tsunami hazards for not only Cascadia, but other subduction zones that have limited seismic observations but a wealth of geodetic inference. 
    more » « less
  5. Abstract The dip angles of slabs are among the clearest characteristics of subduction zones, but the factors that control them remain obscure. Here, slab dip angles and subduction parameters, including subduction duration, the nature of the overriding plate, slab age, and convergence rate, are determined for 153 transects along subduction zones for the present day. We present a comprehensive tabulation of subduction duration based on isotopic ages of arc initiation and stratigraphic, structural, plate tectonic and seismic indicators of subduction initiation. We present two ages for subduction zones, a long‐term age and a reinitiation age. Using cross correlation and multivariate regression, we find that (1) subduction duration is the primary parameter controlling slab dips with slabs tending to have shallower dips at subduction zones that have been in existence longer; (2) the long‐term age of subduction duration better explains variation of shallow dip than reinitiation age; (3) overriding plate nature could influence shallow dip angle, where slabs below continents tend to have shallower dips; (4) slab age contributes to slab dip, with younger slabs having steeper shallow dips; and (5) the relations between slab dip and subduction parameters are depth dependent, where the ability of subduction duration and overriding plate nature to explain observed variation decreases with depth. The analysis emphasizes the importance of subduction history and the long‐term regional state of a subduction zone in determining slab dip and is consistent with mechanical models of subduction. 
    more » « less