Ultra-thin clay layers facilitate seismic slip in carbonate faults

Fault rock heterogeneity can produce fault weakness and reduce fault stability

Geological heterogeneity is abundant in crustal fault zones; however, its role in controlling the mechanical behaviour of faults is poorly constrained. Here, we present laboratory friction experiments on laterally heterogeneous faults, with patches of strong, rate-weakening quartz gouge and weak, rate-strengthening clay gouge. The experiments show that the heterogeneity leads to a significant reduction in strength and frictional stability in comparison to compositionally identical faults with homogeneously mixed gouges. We identify a combination of weakening effects, including smearing of the weak clay; differential compaction of the two gouges redistributing normal stress; and shear localization producing stress concentrations in the strong quartz patches. The results demonstrate that geological heterogeneity and its evolution can have pronounced effects on fault strength and stability and, by extension, on the occurrence of slow-slip transients versus earthquake ruptures and the characteristics of the resulting events, and should be further studied in lab experiments and earthquake source modelling.

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.

Molecular origin of sliding friction and flash heating in rock and heterogeneous materials

It is generally believed that earthquakes occur when faults weaken with increasing slip rates. An important factor contributing to this phenomenon is the faults’ dynamic friction, which may be reduced during earthquakes with high slip rates, a process known as slip-rate weakening. It has been hypothesized that the weakening phenomenon during fault slip may be activated by thermal pressurization of pores’ fluid and flash heating, a microscopic phenomenon in which heat is generated at asperity contacts due to high shear slip rates. Due to low thermal conductivity of rock, the heat generated at the contact points or surfaces cannot diffuse fast enough, thus concentrating at the contacts, increasing the local contact temperature, and reducing its frictional shear strength. We report the results of what we believe to be the first molecular scale study of the decay of the interfacial friction force in rock, observed in experiemntal studies and attributed to flash heating. The magnitude of the reduction in the shear stress and the local friction coefficients have been computed over a wide range of shear velocities V. The molecular simulations indicate that as the interfacial temperature increases, bonds between the atoms begin to break, giving rise to molecular-scale fracture that eventually produces the flash heating effect. The frequency of flash heating events increases with increasing sliding velocity, leaving increasingly shorter times for the material to relax, hence contributing to the increased interfacial temperature. If the material is thin, the heat quickly diffuses away from the interface, resulting in sharp decrease in the temperature immediately after flash heating. The rate of heat transfer is reduced significantly with increasing thickness, keeping most of the heat close to the interface and producing weakened material. The weakening behavior is demonstrated by computing the stress–strain diagram. For small strain rates there the frictional stress is essentially independent of the materials’ thickness. As the strain rate increases, however, the dependence becomes stronger. Specifically, the stress–strain diagrams at lower velocities V manifest a pronounced strength decrease over small distances, whereas they exhibit progressive increase in the shear stress at higher V, which is reminiscent of a transition from ductile behavior at high velocities to brittle response at low velocities.

Diagenetic and shear-induced transitions of frictional strength of carbon-bearing faults and their implications for earthquake rupture dynamics in subduction zones

Subduction-related diagenetic reactions affect fault strength and are thus important for understanding earthquake rupture dynamics in subduction zones. Carbonaceous material (CM) is found worldwide in active plate-boundary and intracontinental faults, yet the effect of its transformation on frictional strength and rupture dynamics remains unknown. We conducted high-velocity friction experiments together with organochemical analyses on CM in the form of lignite, bituminous coal, anthracite and graphite. Results clearly show that an increase in CM maturity and crystallinity leads to a decrease in the peak friction coefficient (from 0.5 to 0.2). We also infer that friction applied to low-grade CM increases its maturity, but friction applied to high-grade CM reduces its maturity. These findings suggest that both diagenetic and shear-induced transformations of CM strongly affect the frictional strength of CM-bearing faults, potentially affecting the depth-dependences of frictional strength and rupture dynamics on plate-subduction faults.

Fault Dynamics of the 1999 Chi-Chi earthquake: clues from nanometric geochemical analysis of fault gouges

The 1999 Chi-Chi, Taiwan earthquake (Mw 7.6) produced surface ruptures for about 90 km along the north-south trending Chelungpu fault, with surface displacements of up to 12 m. Based on the combination of nanoscopic investigation and geochemistry analysis of core samples from a 450 m long inclined borehole drilled through the slip zone, we suggest the dynamical processes that likely occurred in the northern portion of the Chelungpu fault during the faulting. Our analysis revealed that the frictional heating could have reached 1200 °C, which would cause most of the siderite in the fault gouge to evaporate, resulting in a large amount of nano-size siderite grains with a mean diameter 20 nm. These nano grains could have acted as a mechanical lubricant to reduce the dynamic frictional resistance during sliding, giving rise to the large but smooth type of slipping seen in the north.