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

Subduction megathrust creep governed by pressure solution and frictional-viscous flow

Subduction megathrust slip speeds range from slow creep at plate convergence rates (centimetres per year) to seismic slip rates (metres per second) in the largest earthquakes on Earth. The deformation mechanisms controlling whether fast slip or slow creep occurs, however, remain unclear. Here, we present evidence that pressure solution creep (fluid-assisted stress driven mass transfer) is an important deformation mechanism in megathrust faults. We quantify megathrust strength using a laboratory-constrained microphysical model for fault friction, involving viscous pressure solution and frictional sliding. We find that at plate-boundary deformation rates, aseismic, frictional–viscous flow is the preferred deformation mechanism at temperatures above 100 °C. The model thus predicts aseismic creep at temperatures much cooler than the onset of crystal plasticity, unless a boundary condition changes. Within this model framework, earthquakes may nucleate when a local increase in strain rate triggers velocity-weakening slip, and we speculate that slip area and event magnitude increase with increasing spacing of strong, topographically derived irregularities in the subduction interface

Mid-crustal shear zone development under retrograde conditions: pressure–temperature–fluid constraints from the Kuckaus Mylonite Zone, Namibia

The Kuckaus Mylonite Zone (KMZ) forms part of the larger Marshall Rocks–Pofadder shear zone system, a 550 km-long, crustal-scale strike-slip shear zone system that is localized in high-grade granitoid gneisses and migmatites of the Namaqua Metamorphic Complex. Shearing along the KMZ occurred ca. 40 Ma after peak granulite-facies metamorphism during a discrete tectonic event and affected the granulites that had remained at depth since peak metamorphism. Isolated lenses of metamafic rocks within the shear zone allow the P–T–fluid conditions under which shearing occurred to be quantified. These lenses consist of an unsheared core that preserves relict granulite-facies textures and is mantled by a schistose collar and mylonitic envelope that formed during shearing. All three metamafic textural varieties contain the same amphibolite-facies mineral assemblage, from which calculated pseudosections constrain the P–T conditions of deformation at 2.7–4.2 kbar and 450–480 °C, indicating that deformation occurred at mid-crustal depths through predominantly viscous flow. Calculated T–MH2O diagrams show that the mineral assemblages were fluid saturated and that lithologies within the KMZ must have been rehydrated from an external source and retrogressed during shearing. Given that the KMZ is localized in strongly dehydrated granulites, the fluid must have been derived from an external source, with fluid flow allowed by local dilation and increased permeability within the shear zone. The absence of pervasive hydrothermal fractures or precipitates indicates that, even though the KMZ was fluid bearing, the fluid/rock ratio and fluid pressure remained low. In addition, the fluid could not have contributed to shear zone initiation, as an existing zone of enhanced permeability is required for fluid infiltration. We propose that, following initiation, fluid infiltration caused a positive feedback that allowed weakening and continued strain localization. Therefore, the main contribution of the fluid was to produce retrograde mineral phases and facilitate grain-size reduction. Features such as tectonic tremor, which are observed on active faults under similar conditions as described here, may not require high fluid pressure, but could be explained by reaction weakening under hydrostatic fluid pressure conditions

Mid-crustal shear zone development under retrograde conditions: pressure–temperature–fluid constraints from the Kuckaus Mylonite Zone, Namibia

The Kuckaus Mylonite Zone (KMZ) forms part of the larger Marshall Rocks–Pofadder shear zone system, a 550 km-long, crustal-scale strike-slip shear zone system that is localized in high-grade granitoid gneisses and migmatites of the Namaqua Metamorphic Complex. Shearing along the KMZ occurred ca. 40 Ma after peak granulite-facies metamorphism during a discrete tectonic event and affected the granulites that had remained at depth since peak metamorphism. Isolated lenses of metamafic rocks within the shear zone allow the P–T–fluid conditions under which shearing occurred to be quantified. These lenses consist of an unsheared core that preserves relict granulite-facies textures and is mantled by a schistose collar and mylonitic envelope that formed during shearing. All three metamafic textural varieties contain the same amphibolite-facies mineral assemblage, from which calculated pseudosections constrain the P–T conditions of deformation at 2.7–4.2 kbar and 450–480 °C, indicating that deformation occurred at mid-crustal depths through predominantly viscous flow. Calculated T–MH2O diagrams show that the mineral assemblages were fluid saturated and that lithologies within the KMZ must have been rehydrated from an external source and retrogressed during shearing. Given that the KMZ is localized in strongly dehydrated granulites, the fluid must have been derived from an external source, with fluid flow allowed by local dilation and increased permeability within the shear zone. The absence of pervasive hydrothermal fractures or precipitates indicates that, even though the KMZ was fluid bearing, the fluid/rock ratio and fluid pressure remained low. In addition, the fluid could not have contributed to shear zone initiation, as an existing zone of enhanced permeability is required for fluid infiltration. We propose that, following initiation, fluid infiltration caused a positive feedback that allowed weakening and continued strain localization. Therefore, the main contribution of the fluid was to produce retrograde mineral phases and facilitate grain-size reduction. Features such as tectonic tremor, which are observed on active faults under similar conditions as described here, may not require high fluid pressure, but could be explained by reaction weakening under hydrostatic fluid pressure conditions

Sources and effects of fluids in continental retrograde shear zones: Insights from the Kuckaus Mylonite Zone, Namibia

Midcrustal rocks in retrograde metamorphic settings are typically H2O-undersaturated and fluid-absent and have low permeability. Exhumed continental retrograde faults, nonetheless, show evidence for the operation of fluid-mediated weakening mechanisms during deformation at midcrustal conditions. To explore the origin and effects of fluids in retrograde faults, we study the Kuckaus Mylonite Zone (KMZ), an exhumed crustal-scale, strike-slip shear zone in the southern Namibian Namaqua Metamorphic Complex. The KMZ deformed quartzofeldspathic migmatised gneisses at midcrustal retrograde conditions (450-480°C, 270-420 MPa) in the Mesoproterozoic, 40 Ma after granulite facies peak metamorphism at 825°C and 550 MPa. The mylonites contain fully hydrated retrograde mineral assemblages, predominantly adjacent to anastomosing high-strain zones, providing evidence of local H2O saturation and fluid presence during deformation. Whole rock and quartz vein δ18O values suggest that at least some of the fluids were meteoric in origin. The rocks across the shear zone retain the effect of different protoliths, implying little effect of fluid-rock interaction on whole rock major element chemistry. Together with a general scarcity of quartz veins, this suggests that fluid/rock ratios remained low in the KMZ. However, even small amounts of H2O allowed reaction weakening and diffusion-precipitation, followed by growth and alignment of phyllosilicates. In the ultramylonites, a fine grain size in the presence of fluids allowed for grain size sensitive creep. We conclude that the influx of even small volumes of fluids into retrograde shear zones can induce drastic weakening by facilitating grain size sensitive creep and retrograde reactions. In retrograde settings, these reactions consume fluids, and therefore elevated fluid pressures will only be possible after considerable weakening has already occurred. Our findings imply that the range of seismic styles recently documented at active retrograde transform faults may not require high fluid pressures but could also arise from other local weakening mechanisms

San Andreas Fault tremor and retrograde metamorphism

[1] Tectonic tremor is an enigmatic low-frequency seismic phenomenon mainly observed in subduction zones, but also documented along the deep extension of the central San Andreas Fault. The physical mechanisms behind this unusual seismic event are not yet determined for any tectonic setting; however, low effective stress conditions arising from metamorphic fluid production are commonly inferred for subduction-related tremor. We investigate the petrologic conditions at which the San Andreas tectonic tremor is inferred to occur through calculations of the pressure – temperature – time evolution of stable mineral assemblages and their water content in the dominant lithologies of the Franciscan Complex. We find that tremor locations around Parkfield and Cholame are currently experiencing retrograde metamorphic conditions. Within the temperature-depth conditions of observed tremor activity, at approximately 500°C and 20 km depth, several mineralogical transitions may occur in cooling greywacke and mafic rocks, leading to localised, significant removal of free water and an associated volume decrease. This indicates that, contrary to subduction-related tremor, tremor on the San Andreas Fault is not linked to prograde, crustal metamorphic fluid production within the fault zone; rather it might be related to mantle-derived fluids from below the tremor zone, and/or fault zone weakening that occurs as phyllosilicates replace more competent and granular mineral phases

Mixed brittle and viscous strain localisation in pelagic sediments seaward of the Hikurangi margin, New Zealand

Calcareous‐pelagic input sediments are present at several subduction zones and deform differently to their siliciclastic counterparts. We investigate deformation in calcareous‐pelagic sediments drilled ~20 km seaward of the Hikurangi megathrust toe at Site U1520 during IODP Expeditions 372 and 375. Clusters of normal faults and subhorizontal stylolites in the sediments indicate both brittle faulting and viscous pressure solution operated at 150°C where frictional (possibly seismic) slip likely predominates. Plain Language Summary The type of sediments entering subduction zones will influence the way the plates in the subduction zone slide past one another. We looked at limestones in sediments drilled before they reach the subduction zone and found that because of the pressure they are under, they begin to crack and dissolve at very shallow depths. Most of the dissolution happens on thin layers where it concentrates clay by removing other, more soluble minerals. We compare how much vertical shortening we see in the sediments to a computer model. The model overestimates vertical shortening over the history of the sediment unless either high pressure fluids reduce the pressure felt by the sediments, or dissolution is governed by the largest sediment grains rather than their average size. Dissolving and cracking make the sediments weaker by concentrating soft materials such as clay. When these sediments enter the subduction zone, the two plates might slip past one another more easily on these weak regions, possibly during slow slip events

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

Quartz vein formation by local dehydration embrittlement along the deep, tremorgenic subduction thrust interface

Hydrothermal quartz veins are ubiquitous in exhumed accretionary complexes, including the Namibian Damara belt. Here, subduction-related deformation occurred at temperatures ≤550 °C, and vein geometry is consistent with plate interface shear, low effective normal stresses, and mixed-mode deformation. Quartz vein δ18O values relative to Standard Mean Ocean Water (SMOW) range from 9.4‰ to 17.9‰ (n = 30), consistent with precipitation from metamorphic fluids. A dominant subset of quartz veins away from long-lived high-strain zones and basaltic slivers have δ18O values in a smaller range of 14.9‰ ± 1‰, requiring precipitation from a fluid with δ18O of 12‰ ± 1‰ at 470–550 °C. This uniform fluid isotope value is consistent with progressive local breakdown of chlorite allowing extensive hydrofracture at temperatures typical of the plastic regime. In active subduction zones, brittle deformation within the plastic regime is inferred from observations of tectonic tremor, a noise-like seismic signal including overlapping low- and very low-frequency earthquakes, which occurs below the seismogenic zone. Both tremor and hydrothermal veins correlate with zones of inferred high fluid pressure, could represent a mixture of shear and dilatant failure, and may therefore be controlled by episodic hydrofracturing within a dominantly plastic and aseismic regime

Embrittlement within viscous shear zones across the base of the subduction thrust seismogenic zone

Geophysical observations indicate that patches of localised fracturing occur within otherwise viscous regions of subduction plate boundaries. These observations place uncertainty on the possible down-dip extent of the seismogenic zone, and as a result the maximum magnitude of subduction thrust earthquakes. However, the processes controlling where and how localised fracturing occurs within otherwise viscous shear zones are unclear. We examined three exposures of exhumed plate boundary on Kyushu, Japan, which contain subducted sediments and hydrated oceanic crust deformed at ∼300 to ∼500 °C. These exposures preserve subduction-related viscous deformation, which in two of the studied exposures has a mutually overprinting relationship with quartz veins, indicating localised cyclical embrittlement. Where observed, fractures are commonly near lithological contacts that form viscosity contrasts. Mineral equilibrium calculations for a metabasalt composition indicate that exposures showing cyclical embrittlement deformed at pressure-temperature conditions near dehydration reactions that consume prehnite and chlorite. In contrast, dominantly viscous deformation occurred at intervening pressure-temperature conditions. We infer that at conditions close to metamorphic dehydration reactions, only small stress perturbations are required for transient embrittlement, driven by localised dehydration reactions reducing effective stress, and/or locally increased shear stresses along rheological contrasts. Our results show that the protolith composition of the subducting oceanic lithosphere controls the locations and magnitudes of dehydration reactions, and the viscosity of metamorphosed oceanic crust. Therefore, compositional variations might drive substantial variations in slip style