Seismological evidence for a multifault network at the subduction interface

Subduction zones generate the largest earthquakes on Earth, yet their detailed structure, and its influence on seismic and aseismic slip, remains poorly understood. Geological studies of fossil subduction zones characterize the seismogenic interface as a 100 m-1 km thick zone1-3 in which deformation occurs mostly on metres-thick faults1,3-6. Conversely, seismological studies, with their larger spatial coverage and temporal resolution but lower spatial resolution, often image the seismogenic interface as a kilometres-wide band of seismicity7. Thus, how and when these metre-scale structures are active at the seismic-cycle timescale, and what influence they have on deformation is not known. Here we detect these metres-thick faults with seismicity and show their influence on afterslip propagation. Using a local three-dimensional velocity model and dense observations of more than 1,500 double-difference relocated earthquakes in Ecuador, we obtain an exceptionally detailed image of seismicity, showing that earthquakes occur sometimes on a single plane and sometimes on several metres-thick simultaneously active subparallel planes within the plate interface zone. This geometrical complexity affects afterslip propagation, demonstrating the influence of fault continuity and structure on slip at the seismogenic interface. Our findings can therefore help to create more realistic models of earthquake rupture, aseismic slip and earthquake hazard in subduction zones.

Heterogeneous Subgreenschist Deformation in an Exhumed Sediment-Poor Mélange

Many described subduction complexes (or mélanges) exhumed from seismogenic depths comprise thick, turbidite-dominated sequences with deformed zones containing clasts or boudins of more competent sandstone and/or basalt. In contrast, many active subduction zones have a relatively small thickness of sedimentary inputs (<2 km), turbidite sequences are commonly accreted rather than subducted, and the role of pelagic sediments and basalt (lavas and hyaloclastites) in the deforming zone near the plate interface at <20 km depth is poorly understood. Field investigation of Neoproterozoic oceanic sequences accreted in the Gwna Complex, Anglesey, UK, reveals repeated lenticular slices of variably sampled ocean plate stratigraphy (OPS) bounded by thin mélange-bearing shear zones. Mélange matrix material is derived from adjacent OPS lithologies and is either dominantly illitic, likely derived from altered siliciclastic sediment, or chloritic, likely derived from altered volcanics. In the illitic mélange, mutually cross-cutting phyllosilicate foliation and variably deformed chlorite-quartz-calcite veins suggest ductile creep was cyclically punctuated by transient, localized fluid pulses. Chlorite thermometry indicates the veins formed at 260 ± 10°C. In the chloritic mélange, recrystallized through-going calcite veins are deformed to shear strains of 4-5 within a foliated chlorite matrix, suggesting calcite veins in subducting volcanics may localize deformation in the seismogenic zone. Shear stress-strain rate curves constructed using existing empirical relationships in a simplified shear zone geometry predict that slip velocities varied depending on pore fluid pressure; models predict slow slip velocities preferentially by frictional sliding in chlorite, at pore fluid pressures greater than hydrostatic but less than lithostatic.

Is complex fault zone behaviour a reflection of rheological heterogeneity?

Fault slip speeds range from steady plate boundary creep through to earthquake slip. Geological descriptions of faults range from localized displacement on one or more discrete planes, through to distributed shearing flow in tabular zones of finite thickness, indicating a large range of possible strain rates in natural faults. We review geological observations and analyse numerical models of two-phase shear zones to discuss the degree and distribution of fault zone heterogeneity and effects on active fault slip style. There must be certain conditions that produce earthquakes, creep and slip at intermediate velocities. Because intermediate slip styles occur over large ranges in temperature, the controlling conditions must be effects of fault properties and/or other dynamic variables. We suggest that the ratio of bulk driving stress to frictional yield strength, and viscosity contrasts within the fault zone, are critical factors. While earthquake nucleation requires the frictional yield to be reached, steady viscous flow requires conditions far from the frictional yield. Intermediate slip speeds may arise when driving stress is sufficient to nucleate local frictional failure by stress amplification, or local frictional yield is lowered by fluid pressure, but such failure is spatially limited by surrounding shear zone stress heterogeneity. This article is part of a discussion meeting issue 'Understanding earthquakes using the geological record'.

Flow-to-Friction Transition in Simulated Calcite Gouge: Experiments and Microphysical Modeling

A (micro)physical understanding of the transition from frictional sliding to plastic or viscous flow has long been a challenge for earthquake cycle modeling. We have conducted ring-shear deformation experiments on layers of simulated calcite fault gouge under conditions close to the frictional-to-viscous transition previously established in this material. Constant velocity (v) and v-stepping tests were performed, at 550°C, employing slip rates covering almost 6 orders of magnitude (0.001-300 μm/s). Steady-state sliding transitioned from (strong) v-strengthening, flow-like behavior to v-weakening, frictional behavior, at an apparent "critical" velocity (v cr ) of ~0.1 μm/s. Velocity-stepping tests using v < v cr showed "semi-brittle" flow behavior, characterized by high stress sensitivity ("n-value") and a transient response resembling classical frictional deformation. For v ≥ v cr , gouge deformation is localized in a boundary shear band, while for v < v cr , the gouge is well-compacted, displaying a progressively homogeneous structure as the slip rate decreases. Using mechanical data and post-mortem microstructural observations as a basis, we deduced the controlling shear deformation mechanisms and quantitatively reproduced the steady-state shear strength-velocity profile using an existing micromechanical model. The same model also reproduces the observed transient responses to v-steps within both the flow-like and frictional deformation regimes. We suggest that the flow-to-friction transition strongly relies on fault (micro)structure and constitutes a net opening of transient microporosity with increasing shear strain rate at v < v cr , under normal stress-dependent or "semi-brittle" flow conditions. Our findings shed new insights into the microphysics of earthquake rupture nucleation and dynamic propagation in the brittle-to-ductile transition zone.