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 June 8, 2026

Title: Impact of Spontaneous Generation of Splay Faults on Stress Heterogeneity in Accreting Sediments
ABSTRACT:We model the progressive generation and propagation of faults in accreting sediments and study how this process affects stress and deformation in the accretionary wedge. We develop large-scale evolutionary drained geomechanical models. We simulate sediment behavior using a porous elasto-plastic constitutive formulation that incorporates the effect of both mean and shear stress to compaction. We integrate a FEM-DEM framework that enables faulting when the deviatoric plastic strain exceeds a specified threshold over a specified length. New faults are introduced to the mesh as contact surfaces with a friction angle lower than that of the intact sediment. We find that the weaker faults lead to a decrease in sediment differential stress in zones that extend several hundred meters away from the fault. This introduces a significant heterogeneity in both stress and strength within the accretionary wedge. The heterogeneity propagates seaward as the wedge evolves and new faults are generated. We also explore the effect of fault frictional strength through parametric analyses. We find that weaker faults result in more extensive areas of low differential stress and delay the generation of the next set of faults at the toe of the wedge. Our results offer a significant improvement over previous models of continuum wedge sediments that predict Coulomb failure throughout the wedge. Our study provides insights into the state of stress in faulted accretionary wedges, highlighting the spatial heterogeneity of stress and its potential impact on factors influencing seismic cycles in subduction zones.  more » « less
Award ID(s):
2041496
PAR ID:
10647992
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
ARMA
Date Published:
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT:We study the influence of an upper-plate fault on the stress state of accreting sediments under large-scale deformation. We develop drained evolutionary geomechanical models using the Finite Element program Elfen. We simulate sediments as porous-elastoplastic material, and we model the fault as a pre-existing contact surface with a varying frictional strength that is lower than the intact sediment. The weaker fault results in a decrease in sediment differential stress near and especially seaward of the fault. A significant section of the wedge is affected by this stress variation. In contrast, the stress ratio is that of Coulomb failure further away from the fault. We also show that the maximum principal stress the sediments can support decreases with decreasing fault strength. This study offers a significant improvement over previous models of continuum wedge sediments that predict Coulomb failure throughout the wedge. Our results improve our understanding of near-fault stress state, hence improving our understanding of seismic hazards in subduction zones and providing practical insights for reservoir quality and the design of safe and economic well trajectories. 1 INTRODUCTIONAccretionary wedges are geological structures that develop at convergent plate boundaries, particularly at subduction zones. They form as a result of offscraping sediments from the subducting oceanic plate, which then accumulates on the leading edge of the overriding plate (e.g., Moore et al., 2011; Buiter et al., 2016; Gao et al., 2018).Understanding the state of stress in accretionary wedges is crucial for the study of earthquake mechanics in subduction zones (e.g., Brodsky et al., 2017; Huffman and Saffer, 2016; Suppe, 2007). It provides insights for earthquake occurrence, large or slow-slip events (e.g., Kodaira et al., 2012; Tobin et al., 2022; Liu and Rice, 2007) and tsunami potential (e.g., Dean et al., 2010; Riedel et al., 2016). In addition, estimates of mean and differential stress, as well as porosity evolution can help improve the assessment of reservoir quality and the design of safe and economic well trajectories (e.g., Morley et al., 2011).The accumulation of sediments in an accretionary wedge resembles a bulldozer gathering snow as it moves. Hence, wedge sediments are often assumed to be at compressional failure (e.g., Davis et al., 1983; Dahlen et al., 1984; Flemings and Saffer, 2018). Several field observations of faulting, folding and lateral compression as seen in seismic images and drilling measurements (e.g., Flemings and Saffer, 2018; Henry et al., 2003; Moore et al., 1990; Westbrook et al., 1988) provide evidence that the accretionary wedge is at a state of failure in compression. 
    more » « less
  2. Abstract We study stress, pressure, and rock properties in evolving accretionary wedges using analytical formulations and geomechanical models. The evolution of the stress state from that imposed by uniaxial burial seaward of the trench to Coulomb failure within the wedge generates overpressure and drives compaction above the décollement. Changes in both mean and shear stress generate overpressure and shear‐induced pressures play a particularly important role in the trench area. In the transition zone between uniaxial burial and Coulomb failure, shear‐induced overpressures increase more than overburden and are higher than footwall pressures. This rapid increase in overpressure reduces the effective normal stress and weakens the plate interface along a zone that onsets ahead of the trench and persists well into the subduction zone. It also drives dewatering at the trench, which enables compaction of the hanging‐wall sediments and a porosity offset at the décollement. Within the accretionary wedge, sediments are at Coulomb failure and the pore pressure response is proportional to changes in mean stress. Low permeability and high convergence rates promote overpressure generation in the wedge, which limits sediment strength. Our results may provide a hydromechanical explanation for a wide range of observed behaviors, including the development of protothrust zones, widespread occurrence of shallow slow earthquake phenomena, and the propagation of large shallow coseismic slip. 
    more » « less
  3. 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
  4. Abstract Geological heterogeneity is abundant in crustal fault zones; however, its role in controlling the mechanical behaviour of faults is poorly constrained. Here, we present laboratory friction experiments on laterally heterogeneous faults, with patches of strong, rate-weakening quartz gouge and weak, rate-strengthening clay gouge. The experiments show that the heterogeneity leads to a significant reduction in strength and frictional stability in comparison to compositionally identical faults with homogeneously mixed gouges. We identify a combination of weakening effects, including smearing of the weak clay; differential compaction of the two gouges redistributing normal stress; and shear localization producing stress concentrations in the strong quartz patches. The results demonstrate that geological heterogeneity and its evolution can have pronounced effects on fault strength and stability and, by extension, on the occurrence of slow-slip transients versus earthquake ruptures and the characteristics of the resulting events, and should be further studied in lab experiments and earthquake source modelling. 
    more » « less
  5. Abstract We numerically investigate the role of plastic strain accumulation on the mechanical response of a planar strike‐slip fault. Our models show that fault‐zone strength significantly impacts the ensuing sequence of earthquakes. Weaker fault zones accumulating more plastic strain promote more complexity in the seismicity pattern through aperiodic earthquake occurrences and intermittent episodes of rupture and arrest. However, if the fault zone strength is high enough, the overall earthquake sequence is characterized by periodic fault‐spanning events. We find that both the fault normal stress and the fault geometric profile evolve throughout the earthquake sequence, suggesting a self‐roughening mechanism. Despite the significant impact of plasticity on the fault response, the width of the plastically deforming region in the fault zone is small compared to the fault length. Our results suggest a rich behavior in dynamically evolving fault zones and support the need for further high‐resolution studies of the highly non‐linear near‐fault region. 
    more » « less