skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Groove Generation and Coalescence on a Large‐Scale Laboratory Fault
Abstract Faults are the products of wear processes acting at a range of scales from nanometers to kilometers. Grooves produced by wear are a first‐order observable feature of preserved surfaces. However, their interpretation is limited by the complex geological histories of natural faults. Here we explore wear processes on faults by forensically examining a large‐scale controlled, laboratory fault which has a maximum offset between the sides of 42 mm and has been reset multiple times for a cumulative slip of approximately 140 mm. We find that on both sides of the fault scratches are formed with lengths that are longer than the maximum offset but less than the cumulative slip. The grooves are explained as a result of interaction with detached gouge rather than as toolmarks produced by an intact protrusion on one side of the fault. The density of grooves increases with normal stress. The experiment has a range of stress of 1–20 MPa and shows a density of 10 grooves/m/MPa in this range. This value is consistent with recent inferences of stress‐dependent earthquake fracture energy of 0.2 J/m2/MPa. At normal stresses above 20 MPa, the grooves are likely to coalesce into a corrugated surface that more closely resembles mature faults. Groove density therefore appears to be an attractive target for field studies aiming to determine the distribution of normal stress on faults. At low stresses the groove spacing can be measured and contrasted with areas where high stresses produce a corrugated surface.  more » « less
Award ID(s):
1763499
PAR ID:
10452393
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
AGU Advances
Volume:
1
Issue:
4
ISSN:
2576-604X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Few natural examples exist where climate’s influence on tectonics is clear. Based on a study of the Sangre de Cristo Mountains in southern Colorado, we argue that climate-driven changes in ice loads affected spatial and temporal slip patterns on the range-front normal fault. Relict glacial features enable the reconstruction of paleoglacier extents and show variable amounts of footwall ice coverage during the Last Glacial Maximum (LGM). Line load models indicate post-LGM ice melting reduced fault clamping stress by ~20–55 kPa at seismic depths. Flexural isostatic modeling shows several meters of footwall uplift due to ice unloading with spatial patterns and magnitudes consistent with post-LGM fault throw measured from offset Holocene and late Pleistocene alluvial fans. Post-LGM fault throw rates are at least a factor of five higher than middle and early Pleistocene rates. We infer that climate-modulated ice-load changes can pace fault clamping stress and slip patterns on range-bounding normal faults. 
    more » « less
  2. Abstract During rock friction experiments at large displacement, room temperature and humidity, and following a hold test, the fracture energy increases approximately as the square of the logarithm of hold duration. While it's been long known that failure strength increases with log hold time, here the slip weakening distance,dh, also increases. The weakening distance increase is large, hundreds of percent change over a few thousand seconds. The initial bare surface and simulated fault gouge experiments were conducted in rotary shear at 25 MPa normal stress, 21 MPa confining stress and at displacements greater than 100 mm. In contrast, initially bare surface experiments at 5 MPa normal stress, unconfined at displacements less than 10 mm show effectively no change indh. We attribute the difference to the presence of an appreciable shear zone that develops due to wear over significant displacements, confined at elevated normal stress. Prior published studies of sheared simulated fault gouge at short displacement show both acknowledged and unacknowledged increases indhthat may relate to our observations. Since natural faults have well‐developed shear zones, the observations have more direct relevance to earthquake nucleation than prior laboratory studies that use short displacement data and focus on frictional strength recovery alone. However, the physics underlying this increase in weakening distance are not known; candidates are compaction (Nakatani, 1998) and delocalization (Sleep et al., 2000). Additional caveats are that these are room temperature and humidity experiments, at a single normal stress that have not yet been reproduced in other laboratories. 
    more » « less
  3. Abstract Many low-angle normal faults (dip ≤30°) accommodate tens of kilometers of crustal extension, but their mechanics remain contentious. Most models for low-angle normal fault slip assume vertical maximum principal stress σ1, leading many authors to conclude that low-angle normal faults are poorly oriented in the stress field (≥60° from σ1) and weak (low friction). In contrast, models for low-angle normal fault formation in isotropic rocks typically assume Coulomb failure and require inclined σ1 (no misorientation). Here, a data-based, mechanical-tectonic model is presented for formation of the Whipple detachment fault, southeastern California. The model honors local and regional geologic and tectonic history and laboratory friction measurements. The Whipple detachment fault formed progressively in the brittle-plastic transition by linking of “minidetachments,” which are small-scale analogs (meters to kilometers in length) in the upper footwall. Minidetachments followed mylonitic anisotropy along planes of maximum shear stress (45° from the maximum principal stress), not Coulomb fractures. They evolved from mylonitic flow to cataclasis and frictional slip at 300–400 °C and ∼9.5 km depth, while fluid pressure fell from lithostatic to hydrostatic levels. Minidetachment friction was presumably high (0.6–0.85), based upon formation of quartzofeldspathic cataclasite and pseudotachylyte. Similar mechanics are inferred for both the minidetachments and the Whipple detachment fault, driven by high differential stress (∼150–160 MPa). A Mohr construction is presented with the fault dip as the main free parameter. Using “Byerlee friction” (0.6–0.85) on the minidetachments and the Whipple detachment fault, and internal friction (1.0–1.7) on newly formed Reidel shears, the initial fault dips are calculated at 16°–26°, with σ1 plunging ∼61°–71° northeast. Linked minidetachments probably were not well aligned, and slip on the evolving Whipple detachment fault probably contributed to fault smoothing, by off-fault fracturing and cataclasis, and to formation of the fault core and fractured damage zone. Stress rotation may have occurred only within the mylonitic shear zone, but asymmetric tectonic forces applied to the brittle crust probably caused gradual rotation of σ1 above it as a result of: (1) the upward force applied to the base of marginal North America by buoyant asthenosphere upwelling into an opening slab-free window and/or (2) basal, top-to-the-NE shear traction due to midcrustal mylonitic flow during tectonic exhumation of the Orocopia Schist. The mechanical-tectonic model probably applies directly to low-angle normal faults of the lower Colorado River extensional corridor, and aspects of the model (e.g., significance of anisotropy, stress rotation) likely apply to formation of other strong low-angle normal faults. 
    more » « less
  4. Abstract Finite-element models of neotectonics require transform faults to rupture seismically even where preseismic shear stresses are low, presumably by dynamic-weakening mechanisms. A long-standing objection is that, if a rupture initiated at an asperity with high static friction stresses, which then transitioned to low dynamic-weakening stresses, local stress drop would be near total and on the order of 80 MPa, which is 4×–40× greater than observed. But the 5 Mw ≥ 7.8 transform earthquakes since 2000 initially ruptured on the branch faults of small net slip (Stein and Bird, 2024). If the slip initiates on a branch fault with different slip physics and no dynamic weakening, this solves the stress-drop problem. We propose that most large shallow earthquakes are hybrid ruptures, which begin on branch faults of small slip with high shear stresses, and then continue propagating on a connected dynamically weakened fault of large slip, even where shear stresses are low. One prediction of this model is that most large shallow ruptures should be unilateral. We test this prediction against the 100 largest (m ≥ 6.49) shallow continental strike-slip earthquakes 1977–2022, using information from the Global Centroid Moment Tensor and International Seismological Centre catalogs. The differences in time and location between the epicenter and the epicentroid define a horizontal “migration” velocity vector for the evolving centroid of each rupture. Early aftershock locations are summarized by a five-parameter elliptical model. Using the geometric relations between these (and mapped traces of active faults) and guided by a symmetrical decision table, we classified 55 ruptures as apparently unilateral, 30 as bilateral, and 15 as ambiguous. Our finding that a majority (55%–70%) of these ruptures are unilateral permits the interpretation that a majority of ruptures are hybrids, both in terms of geometry (branch fault to transform) and in terms of the physics of their fault slip. 
    more » « less
  5. Abstract Many natural faults are believed to consist of velocity weakening (VW) patches surrounded by velocity strengthening (VS) sections. Numerical studies routinely employ this framework to study earthquake sequences including repeating earthquakes. In this laboratory study, we made a VW asperity, of lengthL, from a bare Poly(methyl methacrylate) PMMA frictional interface and coated the surrounding interface with Teflon to make VS fault sections. Behavior of this isolated asperity was studied as a function ofL(ranging from 100 to 400 mm) and the critical nucleation length, , which is inversely proportional to the applied normal stress (2–16 MPa). Consistent with recent numerical simulations, we observed aseismic slip for  < 2, periodic slip for 2 <  < 6, and non‐periodic slip for 10 < . Furthermore, we compared the experiments whereLwas contained by VS material to standard stick‐slip events whereLwas bounded by free surfaces (i.e.,L = the total sample length). The free surface case produced ∼10 times larger slip during stick‐slip events compared to the contained fault ruptures, even with identical . This disparity highlights how standard, complete‐rupture stick‐slip events differ from contained events expected in nature, due to both the free surface conditions and the heterogeneous normal stress along the fault near the free ends, as confirmed by Digital Image Correlation analysis. This study not only introduces the Teflon coating experimental technique for containing laboratory earthquake ruptures, but also highlights the utility of as a predictive parameter for earthquake behavior. 
    more » « less