Inherited Heterogeneities Can Control Viscous Subduction Zone Deformation of Carbonates at Seismogenic Depths

This work links mineral-scale deformation mechanisms with structural evolution during subduction, providing examples showing how grain-scale heterogeneities facilitated viscous creep in calcite at nominally seismogenic temperatures. Carbonates commonly enter subduction zones, either highly concentrated in irregularly distributed sediments or as more distributed precipitates in seafloor volcanics. We present shear zones, localized in calcite veins formed during shallow subduction of calcareous sediment and seafloor volcanics, with viscous shear strains of ≥5. Shear strain localized because secondary phases and chemical variations maintained fine grain sizes in calcite aggregates, activating relatively rapid grain size-sensitive and frictional-viscous creep at temperatures (260 ± 10°C), cooler than predicted from extrapolation of experimental data. Creep at increased strain rates may limit elastic strain accumulation during interseismic periods, reducing the likelihood of large megathrust earthquakes. As shown here for calcite, common inherited natural heterogeneities may induce weakening of viscous mechanisms in other rocks, or at larger scales in the lithosphere.

Rheology of Naturally Deformed Antigorite Serpentinite: Strain and Strain-Rate Dependence at Mantle-Wedge Conditions

Antigorite serpentinite is expected to occur in parts of subduction plate boundaries, and may suppress earthquake slip, but the dominant deformation mechanisms and resultant rheology of antigorite are unclear. An exhumed plate boundary shear zone exposed near Nagasaki, Japan, contains antigorite deformed at 474°C ± 30°C. Observations indicate that a foliation defined by (001) crystal facets developed during plate-boundary shear. Microstructures indicating grain-scale dissolution at high-stress interfaces and precipitation in low-stress regions suggest that dissolution-precipitation creep contributed to foliation development. Analysis of crystal orientations indicate a small contribution from dislocation activity. We suggest a frictional-viscous rheology for antigorite, where dissolution-precipitation produces a foliation defined by (001) crystal facets and acts to resolve strain incompatibilities, allowing for efficient face-to-face sliding between facets. This rheology can not only explain aseismic behavior at ambient plate boundary conditions, but also some of the contrasting behaviors shown by previous field and laboratory studies.

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.

Afterslip Moment Scaling and Variability From a Global Compilation of Estimates

Aseismic afterslip is postseismic fault sliding that may significantly redistribute crustal stresses and drive aftershock sequences. Afterslip is typically modeled through geodetic observations of surface deformation on a case-by-case basis, thus questions of how and why the afterslip moment varies between earthquakes remain largely unaddressed. We compile 148 afterslip studies following 53 M _w 6.0-9.1 earthquakes, and formally analyze a subset of 88 well-constrained kinematic models. Afterslip and coseismic moments scale near-linearly, with a median Spearman's rank correlation coefficient (CC) of 0.91 after bootstrapping (95% range: 0.89-0.93). We infer that afterslip area and average slip scale with coseismic moment as M o 2 / 3 and M o 1 / 3 , respectively. The ratio of afterslip to coseismic moment (M _rel ) varies from <1% to >300% (interquartile range: 9%-32%). M _rel weakly correlates with M _o (CC: -0.21, attributed to a publication bias), rupture aspect ratio (CC: -0.31), and fault slip rate (CC: 0.26, treated as a proxy for fault maturity), indicating that these factors affect afterslip. M _rel does not correlate with mainshock dip, rake, or depth. Given the power-law decay of afterslip, we expected studies that started earlier and spanned longer timescales to capture more afterslip, but M _rel does not correlate with observation start time or duration. Because M _rel estimates for a single earthquake can vary by an order of magnitude, we propose that modeling uncertainty currently presents a challenge for systematic afterslip analysis. Standardizing modeling practices may improve model comparability, and eventually allow for predictive afterslip models that account for mainshock and fault zone factors to be incorporated into aftershock hazard models.

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'.

Earthquake nucleation in the lower crust by local stress amplification

Deep intracontinental earthquakes are poorly understood, despite their potential to cause significant destruction. Although lower crustal strength is currently a topic of debate, dry lower continental crust may be strong under high-grade conditions. Such strength could enable earthquake slip at high differential stress within a predominantly viscous regime, but requires further documentation in nature. Here, we analyse geological observations of seismic structures in exhumed lower crustal rocks. A granulite facies shear zone network dissects an anorthosite intrusion in Lofoten, northern Norway, and separates relatively undeformed, microcracked blocks of anorthosite. In these blocks, pristine pseudotachylytes decorate fault sets that link adjacent or intersecting shear zones. These fossil seismogenic faults are rarely >15 m in length, yet record single-event displacements of tens of centimetres, a slip/length ratio that implies >1 GPa stress drops. These pseudotachylytes represent direct identification of earthquake nucleation as a transient consequence of ongoing, localised aseismic creep.

This study investigates deep intracontinental earthquakes. Based on field data from exhumed lower crustal pseudotachylytes and mylonites from Lofoten, northern Norway, the authors describe a novel model of earthquake nucleation in the lower crust as a transient consequence of ongoing localised aseismic creep.