The dynamics of accretionary prisms and the processes that take place along subduction interfaces are controlled, in part, by the porosity and fluid overpressure of both the forearc wedge and the sediments transported to the system by the subducting plate. The Hikurangi Margin, located offshore the North Island of New Zealand, is a particularly relevant area to investigate the interplay between the consolidation state of incoming plate sediments, dewatering and fluid flow in the accretionary wedge and observed geodetic coupling and megathrust slip behaviour along the plate interface. In its short geographic extent, the margin hosts a diversity of properties that impact subduction processes and that transition from north to south. Its southernmost limit is characterized by frontal accretion, thick sediment subduction, the absence of seafloor roughness, strong interseismic coupling and deep slow slip events. Here we use seafloor magnetotelluric (MT) and controlled-source electromagnetic (CSEM) data collected along a profile through the southern Hikurangi Margin to image the electrical resistivity of the forearc and incoming plate. Resistive anomalies in the shallow forearc likely indicate the presence of gas hydrates, and we relate deeper forerarc resistors to thrust faulting imaged in colocated seismic reflection data. Because MT and CSEM data are highly sensitive to fluid phases in the pore spaces of seafloor sediments and oceanic crust, we convert resistivity to porosity to obtain a representation of fluid distribution along the profile. We show that porosity predicted by the resistivity data can be well fit by an exponential sediment compaction model. By removing this compaction trend from the porosity model, we are able to evaluate the second-order, lateral changes in porosity, an approach that can be applied to EM data sets from other sedimentary basins. Using this porosity anomaly model, we examine the consolidation state of the incoming plate and accretionary wedge sediments. A decrease in porosity observed in the sediments approaching the trench suggests that a protothrust zone is developing ∼25 km seaward of the frontal thrust. Our data also imply that sediments deeper in the accretionary wedge are slightly underconsolidated, which may indicate incomplete drainage and elevated fluid overpressures of the deep wedge.
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- NSF-PAR ID:
- 10424063
- 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
-
more » « less
SUMMARY -
more » « less
Abstract P ‐to‐S ‐converted waves observed in controlled‐source multicomponent ocean bottom seismometer (OBS) records were used to derive theVp /Vs structure of Cascadia Basin sediments. We usedP ‐to‐S waves converted at the basement to derive an empirical function describing the averageVp /Vs of Cascadia sediments as a function of sediment thickness. We derived one‐dimensional intervalVp /Vs functions from semblance velocity analysis ofS ‐converted intrasediment and basement reflections, which we used to define an empiricalVp /Vs versus burial depth compaction trend. We find that seaward from the Cascadia deformation front,Vp /Vs structure 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 /Vs are 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 /Vs functions of Cascadia Basin sediments. -
null (Ed.)International Ocean Discovery Program (IODP) Expedition 358 was carried out from October 2018 through March 2019 on the D/V Chikyuin an attempt to reach a plate boundary fault zone at seismogenic depths for the first time in scientific ocean drilling. The goal was to extend Hole C0002P from ~2900 to ~5200 meters below seafloor (mbsf) and cross the seismically interpreted main décollement fault zone with logging while drilling (LWD), downhole stress measurements, cuttings sampling, mud gas sampling, and partial coring by drilling a sidetrack to create a new hole (C0002Q). Although drilling reached 3262.5 mbsf, the deepest to date in all of scientific ocean drilling, the effort to drill to and sample the target—the megathrust fault zone—was not successful. Operational challenges in establishing sidetrack holes and advancing them at reasonable rates of penetration limited the new cased hole interval to less than 60 m total at a depth shallower than the previously established casing depth of 2922 mbsf. Combined, the cuttings, logs, and ~60 cm of recovered core from sidetrack Holes C0002Q–C0002T revealed hemipelagic sediments and fine silty turbidites consistent in lithology and physical properties with those recovered in the same depth interval at the same site during Integrated Ocean Drilling Program Expedition 348. Cuttings revealed evidence of only weakly deformed rock, with relatively common calcite veins but few other structural indicators. Because no downhole leak-off tests were made and very little borehole imaging was performed, no further insight into the tectonic context was acquired. After riser drilling at Site C0002 was terminated, drilling at alternate contingency Sites C0024 and C0025 was carried out. Site C0024 targeted the frontal thrust region to sample and log hanging wall rocks and the shallow portion of the décollement zone, and Site C0025 accessed sediments in the Kumano fore-arc basin. At Site C0024, a dedicated logging hole was drilled and a very complete suite of logs were acquired from 0 to 869 mbsf. Preliminary interpretation of log response and images suggests the frontal thrust zone was encountered from about 813 mbsf to the base of the hole, with a zone of notably low resistivity and steep bedding from 850 mbsf to the bottom of the hole. Core samples revealed lithologic units interpreted to be hemipelagic and turbiditic basin fill, trench fill, and Shikoku Basin sediments and encountered deformation potentially associated with a back thrust imaged in seismic reflection data. However, coring had to be terminated at about 620 mbsf, well short of the frontal thrust zone. Site C0025 recovered fore-arc basin sediments underlain by those interpreted to have been deposited in a trench-slope basin setting; no clear transition into older, inner accretionary wedge material was identified during preliminary analysis. Coring from 400 to 571 mbsf yielded datable material and possible evidence for diapiric intrusion of sediments.more » « less
-
null (Ed.)International Ocean Discovery Program (IODP) Expedition 358 was carried out from October 2018 through March 2019 on the D/V Chikyu in an attempt to reach a plate boundary fault zone at seismogenic depths for the first time in scientific ocean drilling. The goal was to extend Hole C0002P from ~2900 to ~5200 meters below seafloor (mbsf) and cross the seismically interpreted main décollement fault zone with logging while drilling, downhole stress measurements, cuttings sampling, mud gas sampling, and partial coring by drilling a sidetrack to create a new hole (C0002Q). Although drilling reached 3262.5 mbsf, the deepest to date in all of scientific ocean drilling, the effort to drill to and sample the target—the megathrust fault zone—was not successful. Operational challenges in establishing sidetrack holes and advancing them at reasonable rates of penetration limited the new cased hole interval to less than 60 m total at a depth shallower than the previously established casing depth of 2922 mbsf. Combined, the cuttings, logs, and ~60 cm of recovered core from sidetrack Holes C0002Q–C0002T revealed hemipelagic sediments and fine silty turbidites consistent in lithology and physical properties with those recovered in the same depth interval at the same site during Integrated Ocean Drilling Program Expedition 348. Cuttings revealed evidence of only weakly deformed rock, with relatively common calcite veins but few other structural indicators. Because no downhole leak-off tests were made and very little borehole imaging was performed, no further insight into the tectonic context was acquired. After riser drilling at Site C0002 was terminated, drilling at alternate contingency Sites C0024 and C0025 was carried out. Site C0024 targeted the frontal thrust region to sample and log hanging wall rocks and the shallow portion of the décollement zone, and Site C0025 accessed sediments in the Kumano fore-arc basin. At Site C0024, a dedicated logging hole was drilled and a very complete suite of logs were acquired from 0 to 869 mbsf. Preliminary interpretation of log response and images suggests the frontal thrust zone was encountered from about 813 mbsf to the base of the hole, with a zone of notably low resistivity and steep bedding from 850 mbsf to the bottom of the hole. Core samples revealed lithologic units interpreted to be hemipelagic and turbiditic basin fill, trench fill, and Shikoku Basin sediments and encountered deformation potentially associated with a back thrust imaged in seismic reflection data. However, coring had to be terminated at about 620 mbsf, well short of the frontal thrust zone. Site C0025 recovered fore-arc basin sediments underlain by those interpreted to have been deposited in a trench-slope basin setting; no clear transition into older, inner accretionary wedge material was identified during the preliminary analysis. Coring from 400 to 571 mbsf yielded datable material and possible evidence for diapiric intrusion of sediments.more » « less
-
more » « less
Abstract High‐resolution three‐dimensional discrete element method (DEM) simulations of sandbox‐scale models of accretionary wedges suggest thrusts follow a variety of propagation processes and orientations depending on a number of factors. These include the stage of development of the wedge (precritical vs. critical), basal friction, and type of thrust (forward vs. backward‐vergent). In terms of propagation processes, two clear mechanisms are identified. The first involves propagation from the décollement to the wedge top, similar to the standard model of thrust propagation seen in many kinematic models, and in the second, thrusts grow downward from an initial nucleation point just below the top surface of the wedge as well as upward from the décollement joining in the middle. In terms of orientation, forward‐vergent thrusts initially form at Roscoe (
θ R = 45° − Ψ/2) or Arthur orientations (θ A = 45° − (ϕ + Ψ)/4), and over greater shortening, rotate into Coulomb orientations (θ C = 45° −ϕ /2). To arrive at these results, a wide array of continuum parameters and fields were extracted from the DEM simulations, including stress, strain, strain rate, kinetic energy, Mohr‐Coulomb parameters, and proximity to yielding using the Drucker‐Prager criterion to visualize thrust nucleation and propagation. Lastly, the advantages and disadvantages of these continuum proxies for discerning failure in the granular assembly are considered, and the spatial and temporal relationship between proximity to yielding and strain localization (both pre‐peak and subsequent persistent shear banding) in the granular model of an accretionary wedge is explored.
