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.

Preventing Instabilities and Inducing Controlled, Slow-Slip in Frictionally Unstable Systems

We propose a theory for preventing instabilities and inducing controlled, slow-slip in frictionally unstable systems, such as the Generalized-Burridge-Knopoff (GBK) model and seismic fault models. We exploit the dependence of friction on pressure and use it as a backdoor for altering the dynamics of the underlying dynamical system. We use the mathematical Theory of Control and, for the first time, we manage to (a) stabilize and restrict chaos in this kind of systems, (b) guarantee slow frictional dissipation and (c) tune the system toward desirable global asymptotic equilibria of lower energy. Our control approach is robust and does not require exact knowledge of the frictional or elastic behavior of the system. Numerical examples of control are given for a Burridge-Knopoff system and a strike-slip fault model obeying rate-and-state friction. GBK models are known to present Self-Organized Critical (SOC) behavior. Therefore, the presented methodology shows an additional example of SOC Control. Even though further developments are necessary before any practical application, we expect our methodology to inspire earthquake mitigation strategies regarding anthropogenic and/or natural seismicity.

Strain localization in planar shear of granular media: the role of porosity and boundary conditions

Shear strain localization into shear bands is associated with velocity weakening instabilities and earthquakes. Here, we simulate steady-state plane-shear flow of numerical granular material (gouge), confined between parallel surfaces. Both constant shear stress and constant strain-rate boundary conditions are tested, and the two types of boundary conditions are found to yield distinct velocity profiles and friction laws. The inertial number, I, exerts the largest control on the layers' behavior, but additional dependencies of friction on normal stress and thickness of the layer are observed under constant stress boundary condition. We find that shear-band localization, which is present in the quasistatic regime ([Formula: see text]) in rate-controlled shear, is absent under stress-controlled loading. In the latter case, flow ceases when macroscopic friction coefficient approaches the quasistatic friction value. The inertial regime that occurs at higher inertial numbers ([Formula: see text]) is associated with distributed shear, and friction and porosity that increase with shear rate (rate-strengthening regime). The finding that shear under constant stress boundary condition produces the inertial, distributed shear but never quasistatic, localized deformation is rationalized based on low fluctuations of shear forces in granular contacts for stress-controlled loading. By examining porosity within and outside a shear band, we also provide a mechanical reason why the transition between quasistatic and inertial shear coincides with the transition between localized and distributed strain.

Nanoscale insight into the relation between pressure solution of calcite and interfacial friction

Pressure solution of carbonate-based rocks participates in many geophysical and geochemical processes, but fundamental knowledge of the interfacial processes is still lacking. By concurrently pressing and sliding two single calcite crystals past each other, the pressure solution rate and the friction force between the crystals were concurrently measured in calcium-carbonate saturated water with an extended surface forces apparatus. These studies reveal that both a decrease and an increase in frictional strength can originate from the pressure-solution of calcite single crystals. By conducting nanoscale force measurements with an atomic force microscope, ion specific effects were unveiled at the level of a single asperity. Pressure solution is promoted when the interfacial water layers of calcite remain undisturbed under stress (e.g. with Ca(II)) and the dissolved ions and water lubricate the interface - a phenomenon called pressure-solution facilitated slip. The mechanically induced disruption of the hydration layers of the calcite surface (e.g. with Mg(II) and low Ni(II) concentration) correlates with the more fluid-like and lubricious behavior of the confined fluid in the absence of pressure solution. Charge neutralization of the calcite surface leads to an abrupt change of calcite's hydration layers, which promotes pressure-solution facilitated slip. This work advances the fundamental understanding of physicochemical interactions occurring at confined surfaces of stressed calcite.

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.

Excitation of San Andreas tremors by thermal instabilities below the seismogenic zone

The relative motion of tectonic plates is accommodated at boundary faults through slow and fast ruptures that encompass a wide range of source properties. Near the Parkfield segment of the San Andreas fault, low-frequency earthquakes and slow-slip events take place deeper than most seismicity, at temperature conditions typically associated with stable sliding. However, laboratory experiments indicate that the strength of granitic gouge decreases with increasing temperature above 350°C, providing a possible mechanism for weakening if temperature is to vary dynamically. Here, we argue that recurring low-frequency earthquakes and slow-slip transients at these depths may arise because of shear heating and the temperature dependence of frictional resistance. Recurring thermal instabilities can explain the recurrence pattern of the mid-crustal low-frequency earthquakes and their correlative slip distribution. Shear heating associated with slow slip is sufficient to generate pseudotachylyte veins in host rocks even when fault slip is dominantly aseismic.

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.

The role of water in fault lubrication

The friction between two adjacent tectonic plates under shear loading may dictate seismic activities. To advance the understanding of mechanisms underlying fault strength, we investigate the frictional characteristics of calcite in an aqueous environment. By conducting single-asperity friction experiments using an atomic force microscope, here we show three pathways of energy dissipation with increasing contact stresses: viscous shear of a lubricious solution film at low normal stresses; shear-promoted thermally activated slip, similar to dry friction but influenced by the hydrated ions localized at the interface; and pressure-solution facilitated slip at sufficiently high stresses and slow sliding velocities, which leads to a prominent decrease in friction. It is also shown that the composition of the aqueous solution affects the frictional response. We use this nanoscale evidence to scrutinize the role of brines on fault behavior and argue that pressure solution provides a weakening mechanism of the fault strength at the level of single-asperity contacts.

Temperature dependence of ice-on-rock friction at realistic glacier conditions

Using a new biaxial friction apparatus, we conducted experiments of ice-on-rock friction in order to better understand basal sliding of glaciers and ice streams. A series of velocity-stepping and slide-hold-slide tests were conducted to measure friction and healing at temperatures between -20°C and melting. Experimental conditions in this study are comparable to subglacial temperatures, sliding rates and effective pressures of Antarctic ice streams and other glaciers, with load-point velocities ranging from 0.5 to 100 µm s-1 and normal stress σn = 100 kPa. In this range of conditions, temperature dependences of both steady-state friction and frictional healing are considerable. The friction increases linearly with decreasing temperature (temperature weakening) from μ = 0.52 at -20°C to μ = 0.02 at melting. Frictional healing increases and velocity dependence shifts from velocity-strengthening to velocity-weakening behaviour with decreasing temperature. Our results indicate that the strength and stability of glaciers and ice streams may change considerably over the range of temperatures typically found at the ice-bed interface.This article is part of the themed issue 'Microdynamics of ice'.

Large-displacement, hydrothermal frictional properties of DFDP-1 fault rocks, Alpine Fault, New Zealand: Implications for deep rupture propagation

The Alpine Fault, New Zealand, is a major plate-bounding fault that accommodates 65-75% of the total relative motion between the Australian and Pacific plates. Here we present data on the hydrothermal frictional properties of Alpine Fault rocks that surround the principal slip zones (PSZ) of the Alpine Fault and those comprising the PSZ itself. The samples were retrieved from relatively shallow depths during phase 1 of the Deep Fault Drilling Project (DFDP-1) at Gaunt Creek. Simulated fault gouges were sheared at temperatures of 25, 150, 300, 450, and 600°C in order to determine the friction coefficient as well as the velocity dependence of friction. Friction remains more or less constant with changes in temperature, but a transition from velocity-strengthening behavior to velocity-weakening behavior occurs at a temperature of T = 150°C. The transition depends on the absolute value of sliding velocity as well as temperature, with the velocity-weakening region restricted to higher velocity for higher temperatures. Friction was substantially lower for low-velocity shearing (V < 0.3 µm/s) at 600°C, but no transition to normal stress independence was observed. In the framework of rate-and-state friction, earthquake nucleation is most likely at an intermediate temperature of T = 300°C. The velocity-strengthening nature of the Alpine Fault rocks at higher temperatures may pose a barrier for rupture propagation to deeper levels, limiting the possible depth extent of large earthquakes. Our results highlight the importance of strain rate in controlling frictional behavior under conditions spanning the classical brittle-plastic transition for quartzofeldspathic compositions.