An official website of the United States government
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
The Evolution of Pore Pressure, Stress, and Physical Properties During Sediment Accretion at Subduction Zones
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
- Award ID(s):
- 2041496
- PAR ID:
- 10424062
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 128
- Issue:
- 6
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Abstract P‐to‐S‐converted waves observed in controlled‐source multicomponent ocean bottom seismometer (OBS) records were used to derive theVp/Vsstructure of Cascadia Basin sediments. We usedP‐to‐Swaves converted at the basement to derive an empirical function describing the averageVp/Vsof Cascadia sediments as a function of sediment thickness. We derived one‐dimensional intervalVp/Vsfunctions from semblance velocity analysis ofS‐converted intrasediment and basement reflections, which we used to define an empiricalVp/Vsversus burial depth compaction trend. We find that seaward from the Cascadia deformation front,Vp/Vsstructure offshore northern Oregon and Washington shows little variability along strike, while the structure of incoming sediments offshore central Oregon is more heterogeneous and includes intermediate‐to‐deep sediment layers of anomalously elevatedVp/Vs. These zones with elevatedVp/Vsare likely due to elevated pore fluid pressures, although layers of high sand content intercalated within a more clayey sedimentary sequence, and/or a higher content of coarser‐grained clay minerals relative to finer‐grained smectite could be contributing factors. We find that the proto‐décollement offshore central Oregon develops within the incoming sediments at a low‐permeability boundary that traps fluids in a stratigraphic level where fluid overpressure exceeds 50% of the differential pressure between the hydrostatic pressure and the lithostatic pressure. Incoming sediments with the highest estimated fluid overpressures occur offshore central Oregon where deformation of the accretionary prism is seaward vergent. Conversely, landward vergence offshore northern Oregon and Washington correlates with more moderate pore pressures and laterally homogeneousVp/Vsfunctions of Cascadia Basin sediments.more » « less
-
Abstract Elevated pore fluid pressure is proposed to contribute to slow earthquakes along shallow subduction plate boundaries. However, the processes that create high fluid pressure, disequilibrium compaction and dehydration reactions, lead to different effective stress paths in fault rocks. These paths are predicted by granular mechanics frameworks to lead to different strengths and deformation modes, yet granular mechanics do not predict their effects on fault stability. To evaluate the role of fluid overpressure on shallow megathrust strength and slip behavior, we conducted triaxial shear experiments on chlorite and celadonite rich gouge layers. Experiments were conducted at constant temperature (130 and 100°C), confining pressure (130 and 140 MPa), and pore fluid pressures (between 10 and 120 MPa). Fluid overpressure due to disequilibrium compaction was simulated by increasing confining and pore fluid pressure synchronously without exceeding the target effective pressure, whereas overpressure due to dehydration reactions was simulated by first loading the sample to a target isotropic effective pressure and then increasing pore fluid pressure to reduce the effective pressure. We find that the effects of fluid pressure and stress path on the mechanical behavior of the chlorite and celadonite gouges can generally be described using the critical state soil mechanics (CSSM) framework. However, path effects are more pronounced and persist to greater displacements in chlorite because its microstructure is more influenced by stress path. Due to its effects on microstructure, the stress path also imparts greater control on the rate‐dependence of chlorite strength, which is not predicted by CSSM.more » « less
-
Because splay faults branch at a steep dip angle from the plate-boundary décollement in an accretionary wedge, their coseismic displacement can potentially result in larger tsunamis with distinct characteristics compared to megathrust-only fault ruptures, posing an enhanced hazard to coastal communities. Elsewhere, there is evidence of coseismic slip on splay faults during many of the largest subduction zone earthquakes, but our understanding of potentially active splay faults and their hazards at the Cascadia subduction zone remains limited. To identify the most recently active splay faults at Cascadia, we conduct stratigraphic and structural interpretations of near-surface deformation in the outer accretionary wedge for the ~400 km along-strike length of the landward vergence zone. We analyze recently acquired high-frequency sparker seismic data and crustal-scale multi-channel seismic data to examine the record of deformation in shallow slope basins and the upper ~1 km of the surrounding accreted sediments and to investigate linkages to deeper décollement structure. We present a new fault map for widest, most completely locked portion of Cascadia from 45 to 48°N latitude, which documents the distribution of faults that show clear evidence of recent late Quaternary activity. We find widespread evidence for active splay faulting up to 30 km landward of the deformation front, in what we define as the active domain, and diminished fault activity landward outside of this zone. The abundance of surface-deforming splay faults in the active outer wedge domain suggests Cascadia megathrust events may commonly host distributed shallow rupture on multiple splay faults located within 30 km of the deformation front.more » « less
