Velocity-dependent slip weakening by the combined operation of pressure solution and foliation development

Cathodoluminescence as a tracing technique for quartz precipitation in low velocity shear experiments

Two simulated gouges (a pure quartz and a quartz-muscovite mixture) were experimentally deformed in a ring shear apparatus at a constant low velocity under hydrothermal conditions favourable for dissolution-precipitation processes. Microstructural analysis using scanning electron microscope cathodoluminescence imaging and cathodoluminescence spectroscopy combined with chemical analysis showed that quartz dissolution and precipitation occurred in both experiments. The starting materials and deformation conditions were chosen so that dissolution-precipitation microstructures could be unambiguously identified from their cathodoluminescence signal. Precipitated quartz was observed as blue luminescent fracture fills and overgrowths with increased Al content relative to the original quartz. In the pure quartz gouge, most of the shear deformation was localized on a boundary-parallel slip surface. Sealing of fractures in a pulverized zone directly adjacent to the slip surface may have helped keeping the deformation localized. In the quartz-muscovite mixture, some evidence was observed of shear-accommodating precipitation of quartz in strain shadows, but predominantly in fractures, elongating the original grains. Precipitation of quartz in fractures implies that the length scale of diffusive mass transfer in frictional-viscous flow is shorter than the length of the quartz domains. Additionally, fracturing might play a more important role than generally assumed. Our results show that cathodoluminescence, especially combined with chemical analysis, is a powerful tool in microstructural analyses of experimentally deformed quartz-bearing material and visualizing quartz precipitation.

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.

Surface microstructures developed on polished quartz crystals embedded in wet quartz sand compacted under hydrothermal conditions

Intergranular pressure solution plays a key role as a deformation mechanism during diagenesis and in fault sealing and healing. Here, we present microstructural observations following experiments conducted on quartz aggregates under conditions known to favor pressure solution. We conducted two long term experiments in which a quartz crystal with polished faces of known crystallographic orientation was embedded in a matrix of randomly oriented quartz sand grains. For about two months an effective axial stress of 15 MPa was applied in one experiment, and an effective confining pressure of 28 MPa in the second. Loading occurred at 350 °C in the presence of a silica-saturated aqueous solution. In the first experiment, quartz sand grains in contact with polished quartz prism (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline10{1 }0$$\end{document}1¯010) faces became ubiquitously truncated against these faces, without indenting or pitting them. By contrast, numerous sand-grain-shaped pits formed in polished pyramidal (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$17\overline{6 }3$$\end{document}176¯3) and (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{4 }134$$\end{document}4¯134) crystal faces in the second experiment. In addition, four-leaved and (in some cases) three-leafed clover-shaped zones of precipitation formed on these prism faces, in a consistent orientation and pattern around individual pits. The microstructures observed in both experiments were interpreted as evidence for the operation of intergranular pressure solution. The dependence of the observed indentation/truncation microstructures on crystal face orientation can be explained by crystallographic control of stress-induced quartz dissolution kinetics, in line with previously published experimental and petrographic data, or possibly by an effect of contact orientation on the stress-induced driving force for pressure solution. This should be investigated in future experiments, providing data and microstructures which enable further mechanism-based analysis of deformation by pressure solution and the effect of crystallographic control on its kinetics in quartz-rich sands and sandstones.

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.

Nanocrystalline Principal Slip Zones and Their Role in Controlling Crustal Fault Rheology

Principal slip zones (PSZs) are narrow (<10 cm) bands of localized shear deformation that occur in the cores of upper-crustal fault zones where they accommodate the bulk of fault displacement. Natural and experimentally-formed PSZs consistently show the presence of nanocrystallites in the <100 nm size range. Despite the presumed importance of such nanocrystalline (NC) fault rock in controlling fault mechanical behavior, their prevalence and potential role in controlling natural earthquake cycles remains insufficiently investigated. In this contribution, we summarize the physical properties of NC materials that may have a profound effect on fault rheology, and we review the structural characteristics of NC PSZs observed in natural faults and in experiments. Numerous literature reports show that such zones form in a wide range of faulted rock types, under a wide range of conditions pertaining to seismic and a-seismic upper-crustal fault slip, and frequently show an internal crystallographic preferred orientation (CPO) and partial amorphization, as well as forming glossy or “mirror-like” slip surfaces. Given the widespread occurrence of NC PSZs in upper-crustal faults, we suggest that they are of general significance. Specifically, the generally high rates of (diffusion) creep in NC fault rock may play a key role in controlling the depth limits to the seismogenic zone.

Ripplocations provide a new mechanism for the deformation of phyllosilicates in the lithosphere

Deformation in Earth's lithosphere is localised in narrow, high-strain zones. Phyllosilicates, strongly anisotropic layered minerals, are abundant in these rocks, where they accommodate much of the strain and play a significant role in inhibiting or triggering earthquakes. Until now it was understood that phyllosilicates could deform only by dislocation glide along layers and could not accommodate large strains without cracking and dilation. Here we show that a new class of atomic-scale defects, known as ripplocations, explain the development of layer-normal strain without brittle damage. We use high-resolution transmission electron microscopy (TEM) to resolve nano-scale bending characteristic of ripplocations in the phyllosilicate mineral biotite. We demonstrate that conjugate delamination arrays are the result of elastic strain energy release due to the accumulation of layer-normal strain in ripplocations. This work provides the missing mechanism necessary to understand phyllosilicate deformation, with important rheological implications for phyllosilicate bearing seismogenic faults and subduction zones.

Velocity-dependent slip weakening by the combined operation of pressure solution and foliation development

Phyllosilicate-bearing faults are characterized by an anastomosing foliation with intervening hard clasts and are believed to be long-term weak structures. Here, I present results of sliding experiments on gouges of 80 wt% quartz and 20 wt% muscovite, sheared under hydrothermal conditions at constant velocity. The results show that significant strengthening occurs over a narrow range of sliding velocities (0.03–1 μm/s). At the lowest velocity investigated, weakness is achieved after a considerable sliding distance of over 20 mm with friction reaching a value of 0.3. Microstructural observations and the application of existing models point to the operation of frictional-viscous flow (FVF), through the serial operation of frictional sliding over a weak foliation and pressure solution of intervening clasts, resulting in low frictional strength and pronounced velocity-strengthening. At higher velocities, grain size reduction becomes dominant in a localized zone, which results in disruption of the foliation and the cessation of the FVF mechanism. In natural settings, earthquakes originating elsewhere on the fault would be rapidly arrested when encountering a foliated part of the fault deforming via FVF. Furthermore, pulses of elevated slip velocity would lead to grain size reduction which would destroy the foliation and cause a long-term strengthening of the fault.