Ultralow friction of carbonate faults caused by thermal decomposition

Interfacial processes underlying the temperature-dependence of friction and wear of calcite single crystals

HYPOTHESIS

This study posits that thermal effects play a substantial role in influencing interfacial processes on calcite, and consequently impacting its mechanochemical properties.

EXPERIMENTS

This work interrogates the temperature-dependence of friction and wear at nanoscale contacts with calcite single crystals at low air humidity (≤ 3-10 % RH) by AFM.

FINDINGS

Three logarithmic regimes for the velocity-dependence of friction are identified. BelowTc ∼ 70 °C, where friction increases with T, there is a transition from velocity-weakening (W1) to velocity-strengthening friction (S1). AboveTc ∼ 70 °C, where friction decreases with T, a second velocity-strengthening friction regime (S0) precedes velocity-weakening friction (W1). The low humidity is sufficient to induce atomic scale changes of the calcite cleavage plane due to dissolution-reprecipitation, and more so at higher temperature and 10 % RH. Meanwhile, the surface softens above Tc -likely owing to lattice dilation, hydration and amorphization. These interfacial changes influence the wear mechanism, which transitions from pit formation to plowing with increase in temperature. Furthermore, the softening of the surface justifies the appearance of the second velocity-strengthening friction regime (S0). These findings advance our understanding of the influence of temperature on the interfacial and mechanochemical processes involving calcite, with implications in natural processes and industrial manufacturing.

Effect of the Fracturing Degree of the Source Rock on Rock Avalanche River-Blocking Behavior Based on the Coupled Eulerian-Lagrangian Technique

In this study, the effect of the fracturing degree of the source rock on rock avalanche river-blocking behavior was investigated. The study included the analysis of mass movement behavior, impulse wave behavior, and the formation of landslide dams. The study included a series of simulations of rock avalanche river-blocking based on the coupled Eulerian-Lagrangian (CEL) technique. Prior to the simulation, a water column collapse model was applied to validate the use of the CEL technique on fluid-structure interaction, and to calibrate the material parameters. The source rock in the rock avalanche simulation was cut by different groups of structural planes, with the number of 0 × 0 × 0, 1 × 1 × 1, 4 × 4 × 4, 9 × 9 × 9, 14 × 14 × 14, 19 × 19 × 19 in each dimension, respectively, to represent different fracturing degrees, on the premise of the same volume and shape of the source rock. The simulation results showed that the sliding mass exhibited structure stabilization, such that the structure of the sliding mass gradually stabilized to a steady status over time, in the mass movement process. The structure stabilization made the center of the sliding mass constantly decrease, and provided a higher speed of movement for the rock avalanches with higher fracturing degrees of the source rock. As for the impulse wave behavior, with the increase in the fracturing degree of the source rock, the maximum kinetic energy of the water decreased, and the maximum height and propagation speed of the impulse waves decreased, which indicated that the maximum height and the propagation speed of the impulse waves were positively correlated with the maximum kinetic energy of the water. In regard to the formation of the landslide dams, when the fracturing degree of the source rock was low, the shape of the landslide dam was very different. With the increase of the fracturing degree of the source rock, the shapes of the landslide dams stabilized, and varied slightly after the fracturing degree of the source rock reached a threshold value.

Fast and Localized Temperature Measurements During Simulated Earthquakes in Carbonate Rocks

The understanding of earthquake physics is hindered by the poor knowledge of fault strength and temperature evolution during seismic slip. Experiments reproducing seismic velocity (∼1 m/s) allow us to measure both the evolution of fault strength and the associated temperature increase due to frictional heating. However, temperature measurements were performed with techniques having insufficient spatial and temporal resolution. Here we conduct high velocity friction experiments on Carrara marble rock samples sheared at 20 MPa normal stress, velocity of 0.3 and 6 m/s, and 20 m of total displacement. We measured the temperature evolution of the fault surface at the acquisition rate of 1 kHz and over a spatial resolution of ∼40 µm with an optical fiber conveying the infrared radiation to a two-color pyrometer. Temperatures up to 1,250°C and low coseismic fault shear strength are compatible with the activation of grain size dependent viscous creep.

Microphysical Modeling of Carbonate Fault Friction at Slip Rates Spanning the Full Seismic Cycle

Laboratory studies suggest that seismogenic rupture on faults in carbonate terrains can be explained by a transition from high friction, at low sliding velocities (V), to low friction due to rapid dynamic weakening as seismic slip velocities are approached. However, consensus on the controlling physical processes is lacking. We previously proposed a microphysically based model (the "Chen-Niemeijer-Spiers" [CNS] model) that accounts for the (rate-and-state) frictional behavior of carbonate fault gouges seen at low velocities characteristic of rupture nucleation. In the present study, we extend the CNS model to high velocities (1 mm/s ≤ V ≤ 10 m/s) by introducing multiple grain-scale deformation mechanisms activated by frictional heating. As velocity and hence temperature increase, the model predicts a continuous transition in dominant deformation mechanisms, from frictional granular flow with partial accommodation by plasticity at low velocities and temperatures, to grain boundary sliding with increasing accommodation by solid-state diffusion at high velocities and temperatures. Assuming that slip occurs in a localized shear band, within which grain size decreases with increasing velocity, the model results capture the main mechanical trends seen in high-velocity friction experiments on room-dry calcite-rich rocks, including steady-state and transient aspects, with reasonable quantitative agreement and without the need to invoke thermal decomposition or fluid pressurization effects. The extended CNS model covers the full spectrum of slip velocities from earthquake nucleation to seismic slip rates. Since it is based on realistic fault structure, measurable microstructural state variables, and established deformation mechanisms, it may offer an improved basis for extrapolating lab-derived friction data to natural fault conditions.

Nanoparticles Observed in a Shear Fracture of Dolomite and a Probable Formation Mechanism

Nanoparticles have been extensively found in brittle faults or ductile shear zones, and their formation is closely related to shear movement along the fault plane. However, the formation mechanisms of these nanoparticles are not yet clear. In this study, dolomite samples were triaxially compressed, at a confining pressure of 200-300 MPa, a temperature between 27 °C and 900 °C and a strain rate of approximately 10^-5s^-1, with a Paterson designed gas medium high-temperature and high-pressure deformation apparatus (HTPDA). Samples deformed at room temperature were characterized by universal microcracks and undulatory extinctions in some grains; when at a temperature between 300 °C and 500 °C, well-developed mechanical twins dominated the microstructure, while at a temperature ≥800 °C, displacements of twin lamellae along a cleavage and a well-developed fracture zone could be seen. Nanoparticles of different shapes were discovered on the slip surfaces of a shear fracture or in microcracks by field emission scanning electron microscopy (FESEM). Nanoparticles on deformed samples under low differential stress were usually of sporadic spherical shapes and uneven distribution; while deformed samples under high differential stress had more dense distributions that were identified. Moreover, grain-overlap and nanofine granulation could be recognized in high strain samples. Based on a mechanical data analysis and microstructural observations, it was suggested that the initial formation of nanoparticles was macroscopically determined by the differential stress subjected to the host rocks, and had nothing to do with temperature; whereas the aggregation morphology of the nanoparticles was related to the temperature during the formation and evolution processes of the nanoparticles.

Weak phases production and heat generation control fault friction during seismic slip

The triggering and magnitude of earthquakes is determined by the friction evolution along faults. Experimental results have revealed a drastic decrease of the friction coefficient for velocities close to the maximum seismic one, independently of the material studied. Due to the extreme loading conditions during seismic slip, many competing physical phenomena occur (like mineral decomposition, nanoparticle lubrication, melting among others) that are typically thermal in origin and are changing the nature of the material. Here we show that a large set of experimental data for different rocks can be described by such thermally-activated mechanisms, combined with the production of weak phases. By taking into account the energy balance of all processes during fault movement, we present a framework that reconciles the data, and is capable of explaining the frictional behavior of faults, across the full range of slip velocities (10^-9 to 10 m/s).

The role of aseismic slip in hydraulic fracturing-induced seismicity

Models for hydraulic fracturing-induced earthquakes in shales typically ascribe fault activation to elevated pore pressure or increased shear stress; however, these mechanisms are incompatible with experiments and rate-state frictional models, which predict stable sliding (aseismic slip) on faults that penetrate rocks with high clay or total organic carbon. Recent studies further indicate that the earthquakes tend to nucleate over relatively short injection time scales and sufficiently far from the injection zone that triggering by either poroelastic stress changes or pore pressure diffusion is unlikely. Here, we invoke an alternative model based on recent laboratory and in situ experiments, wherein distal, unstable regions of a fault are progressively loaded by aseismic slip on proximal, stable regions stimulated by hydraulic fracturing. This model predicts that dynamic rupture initiates when the creep front impinges on a fault region where rock composition favors dynamic and slip rate weakening behavior.

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.

Strain localization in dry sheared granular materials: A compactivity-based approach

Shear banding is widely observed in natural fault zones as well as in laboratory experiments on granular materials. Understanding the dynamics of strain localization under different loading conditions is essential for quantifying strength evolution of fault gouge and energy partitioning during earthquakes and characterizing rheological transitions and fault zone structure changes. To that end, we develop a physics-based continuum model for strain localization in sheared granular materials. The grain-scale dynamics is described by the shear transformation zone (STZ) theory, a nonequilibrium statistical thermodynamic framework for viscoplastic deformation in amorphous materials. Using a finite strain computational framework, we investigate the initiation and growth of complex shear bands under a variety of loading conditions and identify implications for strength evolution and the ductile to brittle transition. Our numerical results show similar localization patterns to field and laboratory observations and suggest that shear zones show more ductile response at higher confining pressures, lower dilatancy, and loose initial conditions. Lower pressures, higher dilatancy, and dense initial conditions favor a brittle response and larger strength drops. These findings shed light on a range of mechanisms for strength evolution in dry sheared granular materials and provide a critical input to physics-based multiscale models of fault zone instabilities.

Carbonaceous Materials in the Fault Zone of the Longmenshan Fault Belt: 2. Characterization of Fault Gouge from Deep Drilling and Implications for Fault Maturity

In recent works on the determination of graphitization of carbonaceous materials (CM) within the principal slip zone (PSZ) of the Longmenshan fault (China), we demonstrated that the formation of graphite, resulted from strain and frictional heating, could be evidence of past seismic slip. Here we utilize Raman Spectroscopy of CM (RSCM) on the CM-bearing gouges in the fault zone of the Longmenshan fault belt, at the borehole depth of 760 m (FZ760) from the Wenchuan earthquake Fault Scientific Drilling project-1 (WFSD-1), to quantitatively characterize CM and further retrieve ancient fault deformation information in the active fault. RSCM shows that graphitization of CM is intense in the fault core with respect to the damage zone, with the graphitized carbon resembling those observed on experimentally formed graphite that was frictionally generated. Importantly, compared to the recognized active fault zone of the Longmenshan fault, the RSCM of measured CM-rich gouge shows a higher degree of graphitization, likely derived from high-temperature-perturbation faulting events. It implies that FZ760 accommodated numerous single-event displacement and/or at higher normal stresses and/or in the absence of pore fluid and/or along a more localized slip surface(s). Because graphite is a well-known lubricant, we surmise that the presence of the higher degree graphitized CM within FZ760 will reduce the fault strength and inefficiently accumulate tectonic stress during the seismic cycle at the current depth, and further infer a plausible mechanism for fault propagation at the borehole depth of 590 m during the Mw 7.9 Wenchuan earthquake.

Surface Forces Apparatus Measurements of Interactions between Rough and Reactive Calcite Surfaces

nm-Range forces acting between calcite surfaces in water affect macroscopic properties of carbonate rocks and calcite-based granular materials and are significantly influenced by calcite surface recrystallization. We suggest that the repulsive mechanical effects related to nm-scale surface recrystallization of calcite in water could be partially responsible for the observed decrease of cohesion in calcitic rocks saturated with water. Using the surface forces apparatus, we simultaneously followed the calcite reactivity and measured the forces in water in two surface configurations: between two rough calcite surfaces (CC) and between rough calcite and a smooth mica surface (CM). We used nm-scale rough, polycrystalline calcite films prepared by atomic layer deposition. We measured only repulsive forces in CC in CaCO₃-saturated water, which was related to roughness and possibly to repulsive hydration effects. Adhesive or repulsive forces were measured in CM in CaCO₃-saturated water depending on calcite roughness, and the adhesion was likely enhanced by electrostatic effects. The pull-off adhesive force in CM became stronger with time, and this increase was correlated with a decrease of roughness at contacts, the parameter which could be estimated from the measured force-distance curves. That suggested a progressive increase of real contact areas between the surfaces, caused by gradual pressure-driven deformation of calcite surface asperities during repeated loading-unloading cycles. Reactivity of calcite was affected by mass transport across nm- to μm-thick gaps between the surfaces. Major roughening was observed only for the smoothest calcite films, where gaps between two opposing surfaces were nm-thick over μm-sized areas and led to force of crystallization that could overcome confining pressures of the order of MPa. Any substantial roughening of calcite caused a significant increase of the repulsive mechanical force contribution.

Slow-slip events in semi-brittle serpentinite fault zones

Slow-slip events are earthquake-like events only with much lower slip rates. While peak coseismic velocities can reach tens of meters per second, slow-slip is on the order of 10^−7±2 m/s and may last for days to weeks. Under the rate-and-state model of fault friction, slow-slip is produced only when the asperity size is commensurate with the critical nucleation size, a function of frictional properties. However, it is unlikely that all subduction zones embody the same frictional properties. In addition to friction, plastic flow of antigorite-rich serpentinite may significantly influence the dynamics of fault slip near the mantle wedge corner. Here, we show that the range of frictional parameters that generate slow slip is widened in the presence of a serpentinized layer along the subduction plate interface. We observe increased stability and damping of fast ruptures in a semi-brittle fault zone governed by both brittle and viscoelastic constitutive response. The rate of viscous serpentinite flow, governed by dislocation creep, is enhanced by high ambient temperatures. When effective viscosity is taken to be dynamic, long-term slow slip events spontaneously emerge. Integration of rheology, thermal effects, and other microphysical processes with rate-and-state friction may yield further insight into the phenomenology of slow slip.

Experimental constraints on dynamic fragmentation as a dissipative process during seismic slip

Various fault damage fabrics, from gouge in the principal slip zone to fragmented and pulverized rocks in the fault damage zone, have been attributed to brittle deformation at high strain rates during earthquake rupture. Past experimental work has shown that there exists a critical threshold in stress-strain rate space through which rock failure transitions from failure along a few discrete fracture planes to intense fragmentation. We present new experimental results on Arkansas Novaculite (AN) and Westerly Granite (WG) in which we quantify fracture surface area produced by dynamic fragmentation under uniaxial compressive loading and examine the controls of pre-existing mineral anisotropy on dissipative processes at the microscale. Tests on AN produced substantially greater new fracture surface area (approx. 6.0 m² g^-1) than those on WG (0.07 m² g^-1). Estimates of the portion of energy dissipated into brittle fracture were significant for WG (approx. 5%), but appeared substantial in AN (10% to as much as 40%). The results have important implications for the partitioning of dissipated energy under extreme loading conditions expected during earthquakes and the scaling of high-speed laboratory rock mechanics experiments to natural fault zones.This article is part of the themed issue 'Faulting, friction and weakening: from slow to fast motion'.

From slow to fast faulting: recent challenges in earthquake fault mechanics

Faults-thin zones of highly localized shear deformation in the Earth-accommodate strain on a momentous range of dimensions (millimetres to hundreds of kilometres for major plate boundaries) and of time intervals (from fractions of seconds during earthquake slip, to years of slow, aseismic slip and millions of years of intermittent activity). Traditionally, brittle faults have been distinguished from shear zones which deform by crystal plasticity (e.g. mylonites). However such brittle/plastic distinction becomes blurred when considering (i) deep earthquakes that happen under conditions of pressure and temperature where minerals are clearly in the plastic deformation regime (a clue for seismologists over several decades) and (ii) the extreme dynamic stress drop occurring during seismic slip acceleration on faults, requiring efficient weakening mechanisms. High strain rates (more than 10⁴ s^-1) are accommodated within paper-thin layers (principal slip zone), where co-seismic frictional heating triggers non-brittle weakening mechanisms. In addition, (iii) pervasive off-fault damage is observed, introducing energy sinks which are not accounted for by traditional frictional models. These observations challenge our traditional understanding of friction (rate-and-state laws), anelastic deformation (creep and flow of crystalline materials) and the scientific consensus on fault operation.This article is part of the themed issue 'Faulting, friction and weakening: from slow to fast motion'.

Ultra-thin clay layers facilitate seismic slip in carbonate faults

Many earthquakes propagate up to the Earth’s surface producing surface ruptures. Seismic slip propagation is facilitated by along-fault low dynamic frictional resistance, which is controlled by a number of physico-chemical lubrication mechanisms. In particular, rotary shear experiments conducted at seismic slip rates (1 ms⁻¹) show that phyllosilicates can facilitate co-seismic slip along faults during earthquakes. This evidence is crucial for hazard assessment along oceanic subduction zones, where pelagic clays participate in seismic slip propagation. Conversely, the reason why, in continental domains, co-seismic slip along faults can propagate up to the Earth’s surface is still poorly understood. We document the occurrence of micrometer-thick phyllosilicate-bearing layers along a carbonate-hosted seismogenic extensional fault in the central Apennines, Italy. Using friction experiments, we demonstrate that, at seismic slip rates (1 ms⁻¹), similar calcite gouges with pre-existing phyllosilicate-bearing (clay content ≤3 wt.%) micro-layers weaken faster than calcite gouges or mixed calcite-phyllosilicate gouges. We thus propose that, within calcite gouge, ultra-low clay content (≤3 wt.%) localized along micrometer-thick layers can facilitate seismic slip propagation during earthquakes in continental domains, possibly enhancing surface displacement.

Preservation of amorphous ultrafine material: A proposed proxy for slip during recent earthquakes on active faults

The criteria for designating an "Active Fault" not only are important for understanding regional tectonics, but also are a paramount issue for assessing the earthquake risk of faults that are near important structures such as nuclear power plants. Here we propose a proxy, based on the preservation of amorphous ultrafine particles, to assess fault activity within the last millennium. X-ray diffraction data and electron microscope observations of samples from an active fault demonstrated the preservation of large amounts of amorphous ultrafine particles in two slip zones that last ruptured in 1596 and 1999, respectively. A chemical kinetic evaluation of the dissolution process indicated that such particles could survive for centuries, which is consistent with the observations. Thus, preservation of amorphous ultrafine particles in a fault may be valuable for assessing the fault's latest activity, aiding efforts to evaluate faults that may damage critical facilities in tectonically active zones.

Near-trench slip potential of megaquakes evaluated from fault properties and conditions

Near-trench slip during large megathrust earthquakes (megaquakes) is an important factor in the generation of destructive tsunamis. We proposed a new approach to assessing the near-trench slip potential quantitatively by integrating laboratory-derived properties of fault materials and simulations of fault weakening and rupture propagation. Although the permeability of the sandy Nankai Trough materials are higher than that of the clayey materials from the Japan Trench, dynamic weakening by thermally pressurized fluid is greater at the Nankai Trough owing to higher friction, although initially overpressured fluid at the Nankai Trough restrains the fault weakening. Dynamic rupture simulations reproduced the large slip near the trench observed in the 2011 Tohoku-oki earthquake and predicted the possibility of a large slip of over 30 m for the impending megaquake at the Nankai Trough. Our integrative approach is applicable globally to subduction zones as a novel tool for the prediction of extreme tsunami-producing near-trench slip.

G: Fracture energy, friction and dissipation in earthquakes

Recent estimates of fracture energy G^' in earthquakes show a power-law dependence with slip u which can be summarized as G^' ∝ u^a where a is a positive real slightly larger than one. For cracks with sliding friction, fracture energy can be equated to G_f : the post-failure integral of the dynamic weakening curve. If the dominant dissipative process in earthquakes is friction, G^' and G_f should be comparable and show a similar scaling with slip. We test this hypothesis by analyzing experiments performed on various cohesive and non-cohesive rock types, under wet and dry conditions, with imposed deformation typical of seismic slip (normal stress of tens of MPa, target slip velocity > 1 m/s and fast accelerations ≈ 6.5 m/s²). The resulting fracture energy G_f is similar to the seismological estimates, with G_f and G^' being comparable over most of the slip range. However, G_f appears to saturate after several meters of slip, while in most of the reported earthquake sequences, G^' appears to increase further and surpasses G_f at large magnitudes. We analyze several possible causes of such discrepancy, in particular, additional off-fault damage in large natural earthquakes.

Scale dependence of rock friction at high work rate

Determination of the frictional properties of rocks is crucial for an understanding of earthquake mechanics, because most earthquakes are caused by frictional sliding along faults. Prior studies using rotary shear apparatus revealed a marked decrease in frictional strength, which can cause a large stress drop and strong shaking, with increasing slip rate and increasing work rate. (The mechanical work rate per unit area equals the product of the shear stress and the slip rate.) However, those important findings were obtained in experiments using rock specimens with dimensions of only several centimetres, which are much smaller than the dimensions of a natural fault (of the order of 1,000 metres). Here we use a large-scale biaxial friction apparatus with metre-sized rock specimens to investigate scale-dependent rock friction. The experiments show that rock friction in metre-sized rock specimens starts to decrease at a work rate that is one order of magnitude smaller than that in centimetre-sized rock specimens. Mechanical, visual and material observations suggest that slip-evolved stress heterogeneity on the fault accounts for the difference. On the basis of these observations, we propose that stress-concentrated areas exist in which frictional slip produces more wear materials (gouge) than in areas outside, resulting in further stress concentrations at these areas. Shear stress on the fault is primarily sustained by stress-concentrated areas that undergo a high work rate, so those areas should weaken rapidly and cause the macroscopic frictional strength to decrease abruptly. To verify this idea, we conducted numerical simulations assuming that local friction follows the frictional properties observed on centimetre-sized rock specimens. The simulations reproduced the macroscopic frictional properties observed on the metre-sized rock specimens. Given that localized stress concentrations commonly occur naturally, our results suggest that a natural fault may lose its strength faster than would be expected from the properties estimated from centimetre-sized rock samples.

Fast-moving dislocations trigger flash weakening in carbonate-bearing faults during earthquakes

Rupture fronts can cause fault displacement, reaching speeds up to several ms⁻¹ within a few milliseconds, at any distance away from the earthquake nucleation area. In the case of silicate-bearing rocks the abrupt slip acceleration results in melting at asperity contacts causing a large reduction in fault frictional strength (i.e., flash weakening). Flash weakening is also observed in experiments performed in carbonate-bearing rocks but evidence for melting is lacking. To unravel the micro-physical mechanisms associated with flash weakening in carbonates, experiments were conducted on pre-cut Carrara marble cylinders using a rotary shear apparatus at conditions relevant to earthquakes propagation. In the first 5 mm of slip the shear stress was reduced up to 30% and CO₂ was released. Focused ion beam, scanning and transmission electron microscopy investigations of the slipping zones reveal the presence of calcite nanograins and amorphous carbon. We interpret the CO₂ release, the formation of nanograins and amorphous carbon to be the result of a shock-like stress release associated with the migration of fast-moving dislocations. Amorphous carbon, given its low friction coefficient, is responsible for flash weakening and promotes the propagation of the seismic rupture in carbonate-bearing fault patches.

Fault rheology beyond frictional melting

During earthquakes, comminution and frictional heating both contribute to the dissipation of stored energy. With sufficient dissipative heating, melting processes can ensue, yielding the production of frictional melts or "pseudotachylytes." It is commonly assumed that the Newtonian viscosities of such melts control subsequent fault slip resistance. Rock melts, however, are viscoelastic bodies, and, at high strain rates, they exhibit evidence of a glass transition. Here, we present the results of high-velocity friction experiments on a well-characterized melt that demonstrate how slip in melt-bearing faults can be governed by brittle fragmentation phenomena encountered at the glass transition. Slip analysis using models that incorporate viscoelastic responses indicates that even in the presence of melt, slip persists in the solid state until sufficient heat is generated to reduce the viscosity and allow remobilization in the liquid state. Where a rock is present next to the melt, we note that wear of the crystalline wall rock by liquid fragmentation and agglutination also contributes to the brittle component of these experimentally generated pseudotachylytes. We conclude that in the case of pseudotachylyte generation during an earthquake, slip even beyond the onset of frictional melting is not controlled merely by viscosity but rather by an interplay of viscoelastic forces around the glass transition, which involves a response in the brittle/solid regime of these rock melts. We warn of the inadequacy of simple Newtonian viscous analyses and call for the application of more realistic rheological interpretation of pseudotachylyte-bearing fault systems in the evaluation and prediction of their slip dynamics.

Dynamic weakening of serpentinite gouges and bare surfaces at seismic slip rates

To investigate differences in the frictional behavior between initially bare rock surfaces of serpentinite and powdered serpentinite ("gouge") at subseismic to seismic slip rates, we conducted single-velocity step and multiple-velocity step friction experiments on an antigorite-rich and lizardite-rich serpentinite at slip rates (V) from 0.003 m/s to 6.5 m/s, sliding displacements up to 1.6 m, and normal stresses (σn ) up to 22 MPa for gouge and 97 MPa for bare surfaces. Nominal steady state friction values (μnss) in gouge at V = 1 m/s are larger than in bare surfaces for all σn tested and demonstrate a strong σn dependence; μnss decreased from 0.51 at 4.0 MPa to 0.39 at 22.4 MPa. Conversely, μnss values for bare surfaces remained ∼0.1 with increasing σn and V. Additionally, the velocity at the onset of frictional weakening and the amount of slip prior to weakening were orders of magnitude larger in gouge than in bare surfaces. Extrapolation of the normal stress dependence for μnss suggests that the behavior of antigorite gouge approaches that of bare surfaces at σn ≥ 60 MPa. X-ray diffraction revealed dehydration reaction products in samples that frictionally weakened. Microstructural analysis revealed highly localized slip zones with melt-like textures in some cases gouge experiments and in all bare surfaces experiments for V ≥ 1 m/s. One-dimensional thermal modeling indicates that flash heating causes frictional weakening in both bare surfaces and gouge. Friction values for gouge decrease at higher velocities and after longer displacements than bare surfaces because strain is more distributed.

KEY POINTS

Gouge friction approaches that of bare surfaces at high normal stressDehydration reactions and bulk melting in serpentinite in < 1 m of slipFlash heating causes dynamic frictional weakening in gouge and bare surfaces.