We investigate the physics of laboratory earthquake precursors in a biaxial shear configuration. We conduct laboratory experiments at room temperature and humidity in which we shear layers of glass ...beads under applied normal loads of 2–8 MPa and with shearing rates of 5–10 µm/s. We show that above ~ 3 MPa load, acoustic emission (AE), and shear microfailure (microslip) precursors exhibit an exponential increase in rate of occurrence, culminating in stick‐slip failure. Precursors take place where the material is in a critical state—still modestly dilating, yet while the macroscopic frictional strength is no longer increasing.
Key Points
Sheared granular media exhibits precursor events to stick‐slip failure
Rate of precursor activity increases exponentially prior to stick‐slip failure
Acoustic emission and microslip precursors are correlated in time
We present a unified analysis of physical properties of cataclastic fault rocks collected from surface exposures of the central Alpine Fault at Gaunt Creek and Waikukupa River, New Zealand. Friction ...experiments on fault gouge and intact samples of cataclasite were conducted at 30–33 MPa effective normal stress (σn′) using a double‐direct shear configuration and controlled pore fluid pressure in a true triaxial pressure vessel. Samples from a scarp outcrop on the southwest bank of Gaunt Creek display (1) an increase in fault normal permeability (k = 7.45 × 10−20 m2 to k = 1.15 × 10−16 m2), (2) a transition from frictionally weak (μ = 0.44) fault gouge to frictionally strong (μ = 0.50–0.55) cataclasite, (3) a change in friction rate dependence (a‐b) from solely velocity strengthening, to velocity strengthening and weakening, and (4) an increase in the rate of frictional healing with increasing distance from the footwall fluvioglacial gravels contact. At Gaunt Creek, alteration of the primary clay minerals chlorite and illite/muscovite to smectite, kaolinite, and goethite accompanies an increase in friction coefficient (μ = 0.31 to μ = 0.44) and fault‐perpendicular permeability (k = 3.10 × 10−20 m2 to k = 7.45 × 10−20 m2). Comminution of frictionally strong (μ = 0.51–0.57) cataclasites forms weaker (μ = 0.31–0.50) foliated cataclasites and fault gouges with behaviors associated with aseismic creep at low strain rates. Combined with published evidence of large magnitude (Mw ∼ 8) surface ruptures on the Alpine Fault, petrological observations indicate that shear failure involved frictional sliding within previously formed, velocity‐strengthening fault gouge.
Key Points
Cataclasites and fault gouges have distinct physical properties
Quantitative XRD results reveal in situ and along‐strike changes in mineralogy
Fault gouges have repeatedly accommodated coseismic slip
Frictional Mechanics of Slow Earthquakes Leeman, J. R.; Marone, C.; Saffer, D. M.
Journal of geophysical research. Solid earth,
September 2018, 2018-09-00, 20180901, Letnik:
123, Številka:
9
Journal Article
Recenzirano
Odprti dostop
Tectonic faults slip in a wide range of modes that span from slow slip events to dynamic rupture. A growing body of observations document this spectrum of failure modes in many geologic settings. ...However, the physical mechanisms that dictate slow slip are not understood. Here we investigate the mechanics of slow slip using carefully controlled laboratory experiments that demonstrate a complete spectrum of slip modes: Laboratory stick‐slip event durations span from seconds to milliseconds, representing the equivalent of failure events that span the range from slow to dynamic earthquakes. The rheological critical stiffness kc is the primary control on the mode of slip, but higher‐order effects including velocity dependence of the frictional rate parameter and critical slip distance also play an important role. We also find that quasi‐dynamic instability results from negligible stress drop near the stability boundary, in tandem with negative feedback during slip acceleration rooted in the rate dependence of kc. Our work shows that a broad spectrum of slip behaviors can arise from a common frictional mechanism modulated by fault zone rheology and elastic properties.
Key Points
Laboratory experiments illuminate a friction‐based explanation for the mechanics underlying the spectrum of fault slip
Fault slip mode is dictated by the ratio of the loading stiffness to a critical rheological stiffness and varies with slip velocity
Slow slip can be explained by velocity dependence of kc and rupture velocity by potential stress drop near the stability boundary
Sand‐shale mélanges from the Kodiak accretionary complex and Shimanto belt of Japan record deformation during underthrusting along a paleosubduction interface in the range 150 to 350 °C. We use ...observations from these mélanges to construct a simple kinetic model that estimates the maximum time required to seal a single fracture as a measure of the rate of fault zone healing. Crack sealing involves diffusive redistribution of Si from mudstones with scaly fabric to undersaturated fluid‐filled cracks in sandstone blocks. Two driving forces are considered for the chemical potential gradient that drives crack sealing: (1) a transient drop in fluid pressure ∆Pf, and (2) a difference in mean stress between scaly slip surfaces in mudstones and cracks in stronger sandstone blocks. Sealing times are more sensitive to mean stress than ∆Pf, with up to four orders of magnitude faster sealing. Sealing durations are dependent on crack spacing, silica diffusion kinetics, and magnitude of the strength contrast between block and matrix, each of which is loosely constrained for conditions relevant to the seismogenic zone. We apply the model to three active subduction zones and find that sealing rates are fastest along Cascadia and several orders of magnitude slower for a given depth along Nicaragua and Tohoku slab‐top geotherms. The model provides (1) a framework for geochemical processes that influence subduction mechanics via crack sealing and shear fabric development and (2) demonstration that kinetically driven mass redistribution during the interseismic period is a plausible mechanism for creating asperities along smooth, sediment‐dominated convergent margins.
Plain Language Summary
Geophysical monitoring of active subduction zones has revealed plate boundary slip behaviors such as creep, slow slip events, and earthquakes that vary spatially and temporally for different plate boundaries and downdip along a given boundary. Fault rocks exposed on land from paleosubduction plate boundaries provide a record of the deformation processes that likely occur during slip along active boundaries, so we review the characteristics of these ancient rocks to develop insight into slip behavior in subduction zones. We find that plate boundary deformation in these cases occurs within a wide fault zone through processes that involve redistribution of silica from shearing mudstones to cracking sandstone blocks. We use a geochemical model to calculate how long it would take to seal a crack by this process and conclude that cracks seal at rates that could influence the earthquake cycle, with rates of crack healing dependent on the temperature structure and the depth where slip occurs. Our results suggest that processes of frictional failure and geochemical healing in downgoing sediments may influence the slip stability along the subduction interface.
Key Points
Microstructural observations of subduction mélange record dissolution along a scaly fabric and diffusion to cracks in sandstone blocks
A silica kinetics model depicts crack healing that is highly dependent on temperature and slab‐top geotherm
Kinetic healing of cracks could modulate fault zone strength during the interseismic period
We explore relationships among bulk modulus, crack density, and permeability through repetitive loading of Marcellus shale. Cumulative cyclic stressing (22–26 MPa with confinement of 24 MPa) is ...applied at a frequency of 0.05 Hz over 100,000 cycles. Changes in acoustic velocities are used to follow changes in dynamic bulk modulus, Poisson ratio, and crack density and to correlate these with bedding-parallel measurements of methane permeability. The shale is represented as an orthotropic elastic medium containing a dominant, noninteracting fracture set separated by thin laminae. An effective continuum model links permeability evolution to the evolution of the bulk modulus and crack density. Bulk modulus is linearly related to crack density by a scaling parameter representing rock fabric and fracture geometry. The Poisson ratio and bulk modulus of the intact, uncracked shale are approximated from our data. We propose a method for tracking permeability evolution of finely laminated shales using acoustic waves.
•Some shales are orthotropic materials with a dominant, noninteracting fracture set oriented parallel to bedding.•When a shale exhibits a dominant fracture set, permeability evolution can be tracked with changing acoustic velocities.•Changes in bulk modulus and Poisson ratio are inversely related to changes in permeability andcrack density parameter.
Plain Language Summary: Shales are increasingly important rocks for oil and gas production and carbon sequestration. Due to their tight and often disconnected internal flow paths, it is often difficult to track changes in permeability—how easily a fluid flows through a rock—that arise within the shale itself. This study shows that for some shales sound waves that can originate on the surface or within bore holes can be used to track changes in permeability. Being able to monitor permeability during liquid or gas production/injection will provide useful information to scientists and engineers for a variety of applications.
SUMMARY
Serpentinites are polymineralic rocks distributed almost ubiquitously across the globe in active tectonic regions. Magnetite-rich serpentinites are found in the low-strain domains of ...serpentinite shear zones, which act as potential sites of nucleation of unstable slip. To assess the potential of earthquake nucleation in these materials, we investigate the link between mechanical properties and fabric of these rocks through a suite of laboratory shear experiments. Our experiments were done at room temperature and cover a range of normal stress and slip velocity from 25 to 100 MPa and 0.3 to 300 µm s−1, respectively. We show that magnetite-rich serpentinites are ideal materials since they display strong sensitivity to the loading rate and are susceptible to nucleation of unstable slip, especially at low forcing slip velocities. We also aim at the integration of mechanical and microstructural results to describe the underlying mechanisms that produce the macroscopic behaviour. We show that mineralogical composition and mineral structure dictates the coexistence of two deformation mechanisms leading to stable and unstable slip. The weakness of phyllosilicates allows for creep during the interseismic phase of the laboratory seismic cycle while favouring the restoration of a load-bearing granular framework, responsible of the nucleation of unstable events. During dynamic slip, fault zone shear fabric determines the mode of slip, producing either asymmetric or Gaussian slip time functions for either fast or slow events. We report rate/state friction parameters and integrate our mechanical data with microstructural observations to shed light on the mechanisms dictating the complexity of laboratory earthquakes. We show that mineralogical and fabric heterogeneities control fault slip behaviour.
To improve our understanding of the complex coupling between circulating fluids and the development of crack damage, we performed flow‐through tests on samples of Etna basalt and Westerly granite ...that were cyclically loaded by deviatoric stresses. The basalt was naturally microfractured, while the relatively crack‐free Westerly granite was thermally pretreated to 500°C and 800°C to generate microcrack damage. Samples were repeatedly loaded and then unloaded under deviatoric stress paths and ultimately to failure. Permeability and water volume content were measured throughout the loading history together with the differential stress. Permeability decreases at low differential stresses and increases at intermediate differential stresses up to a steady value at failure. We use water volume content as a proxy for fluid storage and show that both permeability and storage evolve with damage and evolution of crack density. We use crack models to represent the evolution of permeability as a function of loading state and are able to independently link it to the observed evolution of deformability, used as an independent measure of crack density.
Key Points
Identification of 3 regimes of permeability approaching failure
Identification of the role of cooling rates on crack geometry and connectivity
Prediction of the coupled permeability and microcrack network evolution
We present results from a comprehensive laboratory study of the frictional strength and constitutive properties for all three active strands of the San Andreas Fault penetrated in the San Andreas ...Observatory at Depth (SAFOD). The SAFOD borehole penetrated the Southwest Deforming Zone (SDZ), the Central Deforming Zone (CDZ), both of which are actively creeping, and the Northeast Boundary Fault (NBF). Our results include measurements of the frictional properties of cuttings and core samples recovered at depths of ~2.7 km. We find that materials from the two actively creeping faults exhibit low frictional strengths (μ = ~0.1), velocity‐strengthening friction behavior, and near‐zero or negative rates of frictional healing. Our experimental data set shows that the center of the CDZ is the weakest section of the San Andreas Fault, with μ = ~0.10. Fault weakness is highly localized and likely caused by abundant magnesium‐rich clays. In contrast, serpentine from within the SDZ, and wall rock of both the SDZ and CDZ, exhibits velocity‐weakening friction behavior and positive healing rates, consistent with nearby repeating microearthquakes. Finally, we document higher friction coefficients (μ > 0.4) and complex rate‐dependent behavior for samples recovered across the NBF. In total, our data provide an integrated view of fault behavior for the three active fault strands encountered at SAFOD and offer a consistent explanation for observations of creep and microearthquakes along weak fault zones within a strong crust.
Key Points
Comprehensive laboratory study of faults at SAFOD
Creeping faults are frictionally weak, slide stably, and exhibit no healing
Results consistent with observed creep and microearthquakes
Interseismic recovery of fault strength (healing) following earthquake failure is a fundamental requirement of the seismic cycle and likely plays a key role in determining the stability and slip ...behavior of tectonic faults. We report on laboratory measurements of time‐ and slip‐dependent frictional strengthening for natural and synthetic gouges to evaluate the role of mineralogy in frictional strengthening. We performed slide‐hold‐slide (SHS) shearing experiments on nine natural fault gouges and eight synthetic gouges at conditions of 20 MPa normal stress, 100% relative humidity (RH), large shear strain (~15), and room temperature. Phyllosilicate‐rich rocks show the lowest rates of frictional strengthening. Samples rich in quartz and feldspar exhibit intermediate rates of frictional strengthening, and calcite‐rich gouges show the largest values. Our results show that (1) the rates of frictional strengthening and creep relaxation scale with frictional strength, (2) phyllosilicate‐rich fault gouges have low strength and healing characteristics that promote stable, aseismic creep, (3) most natural fault gouges exhibit intermediate rates of frictional strengthening, consistent with a broad range of fault slip behaviors, and (4) calcite‐rich fault rocks show the highest rates of frictional strengthening, low values of dilation upon reshear, and high frictional strengths, all of which would promote seismogenic behavior.
Key Points
Frictional restrengthening is a fundamental requirement of the earthquake cycle
Performed carefully controlled experiments on 17 natural and synthetic gouges
Gouge mineralogy exerts significant control on healing behavior and mode of slip
Shale gas reservoirs like coalbed methane (CBM) reservoirs are promising targets for geological sequestration of carbon dioxide (CO2). However, the evolution of permeability in shale reservoirs on ...injection of CO2 is poorly understood unlike CBM reservoirs. In this study, we report measurements of permeability evolution in shales infiltrated separately by nonsorbing (He) and sorbing (CO2) gases under varying gas pressures and confining stresses. Experiments are completed on Pennsylvanian shales containing both natural and artificial fractures under nonpropped and propped conditions. We use the models for permeability evolution in coal (Journal of Petroleum Science and Engineering, Under Revision) to codify the permeability evolution observed in the shale samples. It is observed that for a naturally fractured shale, the He permeability increases by approximately 15% as effective stress is reduced by increasing the gas pressure from 1 MPa to 6 MPa at constant confining stress of 10 MPa. Conversely, the CO2 permeability reduces by a factor of two under similar conditions. A second core is split with a fine saw to create a smooth artificial fracture and the permeabilities are measured for both nonpropped and propped fractures. The He permeability of a propped artificial fracture is approximately 2‐ to 3fold that of the nonpropped fracture. The He permeability increases with gas pressure under constant confining stress for both nonpropped and propped cases. However, the CO2 permeability of the propped fracture decreases by between one‐half to one‐third as the gas pressure increases from 1 to 4 MPa at constant confining stress. Interestingly, the CO2 permeability of nonpropped fracture increases with gas pressure at constant confining stress. The permeability evolution of nonpropped and propped artificial fractures in shale is found to be similar to those observed in coals but the extent of permeability reduction by swelling is much lower in shale due to its lower organic content. Optical profilometry is used to quantify the surface roughness. The changes in surface roughness indicate significant influence of proppant indentation on fracture surface in the shale sample. The trends of permeability evolution on injection of CO2 in coals and shales are found analogous; therefore, the permeability evolution models previously developed for coals are adopted to explain the permeability evolution in shales.