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: The Impact of Rheology on the Transition From Stick‐Slip to Creep in a Semibrittle Analog
Abstract Faults can release energy via a variety of different slip mechanisms ranging from steady creep to fast and destructive earthquakes. Tying the rheology of the crust to various slip dynamics is important for our understanding of plate tectonics and earthquake generation. Here we propose that the interplay of fractures and viscous flow leads to a spectrum between stick‐slip and creep. We use an elasto‐visco‐plastic rock analog (Carbopol U‐21) where we vary the yield stress to investigate its impact on slip dynamics in shear experiments. The experiments are performed using a simple shear apparatus, which provides distributed shear across the entire width of the experiment and allows in situ observations of deformation. We record force and displacement during deformation and use time lapse photography to document fracture development. A low yield stress (25 Pa) leads to creep dynamics in the absence of fractures. An intermediate yield stress (144 Pa) leads to the development and interaction of opening (mode I) and shear (mode II) fractures. This interaction leads to a spectrum in slip dynamics ranging from creep to stick‐slip. A high yield stress (357 Pa) results in the development of many mode I fractures and a deformation signal dominated by stick‐slip. These results show that bulk yield stress, fracture formation, and slip dynamics are closely linked and can lead to a continuum between creep and stick‐slip. We suggest that rheology should be considered as an additional mechanism to explain the broad range of slip dynamics in natural faults.  more » « less
Award ID(s):
1547492
PAR ID:
10453442
Author(s) / Creator(s):
 ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
124
Issue:
3
ISSN:
2169-9313
Page Range / eLocation ID:
p. 3144-3154
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    We observed and modeled the elastic, inelastic and time-dependent viscous properties of damaged Berea Sandstone samples to investigate the impact of damage on the rheological properties of rocks. Cylindrical Berea Sandstone plugs were prepared both parallel and perpendicular to bedding. We impacted the samples with Split Hopkinson Pressure Bar to pervasively fracture the specimens at different strain rates. Longitudinal mode-I fractures are dominant in specimens impacted at relatively low strain rates (about 130 /s), whereas shear fractures also form in specimens deformed at high strain rates (up to 250 /s). The damaged rocks were subjected to multiple steps of differential stress loading and hold stages under 15 MPa confining pressure. A key observation is that higher damaged specimens showed greater axial and volumetric creep strain deformation during loading and hold stages. Poisson ratio also increase with increasing damage. We modeled the volumetric strain of the sandstone specimens using a Perzyna viscoplasticity law that employs the Modified Cam Clay model as the yield criterion (Haghighat et al. 2020). We deduced that fractured rocks undergo substantial bulk time-dependent deformation due to volumetric compaction and fracture closure. Damage increase results in decrease of the effective viscosity of the material. 
    more » « less
  2. Abstract The rheology of the crust and mantle and the interaction of viscoelastic flow with seismic/aseismic slip on faults control the state of stress in the lithosphere over multiple seismic cycles. The rheological behavior of rocks is well constrained in a laboratory setting, but thein situproperties of the lithosphere and its lateral variations remain poorly known. Here, we access the lower‐crustal rheology in Southern California by exploiting 8 years of geodetic postseismic deformation following the 2010 El Mayor‐Cucapah earthquake. The data illuminate viscoelastic flow in the lower crust with lateral variations of effective viscosity correlated with the geological province. We show that a Burgers assembly with dashpots following a nonlinear constitutive law can approximate the temporal evolution of stress and strain rate, indicating the activation of nonlinear transient creep before steady‐state dislocation creep. The transient and background viscosities in the lower crust of the Salton Trough are on the order of ~1018and ~1019 Pa s, respectively, about an order of magnitude lower than those in the surrounding regions. We highlight the importance of transient creep, nonlinear flow laws, and lateral variations of rheological properties to capture the entire history of postseismic relaxation following the El Mayor‐Cucapah earthquake. 
    more » « less
  3. Abstract Tectonic faults fail through a spectrum of slip modes, ranging from slow aseismic creep to rapid slip during earthquakes. Understanding the seismic radiation emitted during these slip modes is key for advancing earthquake science and earthquake hazard assessment. In this work, we use laboratory friction experiments instrumented with ultrasonic sensors to document the seismic radiation properties of slow and fast laboratory earthquakes. Stick‐slip experiments were conducted at a constant loading rate of 8 μm/s and the normal stress was systematically increased from 7 to 15 MPa. We produced a full spectrum of slip modes by modulating the loading stiffness in tandem with the fault zone normal stress. Acoustic emission data were recorded continuously at 5 MHz. We demonstrate that the full continuum of slip modes radiate measurable high‐frequency energy between 100 and 500 kHz, including the slowest events that have peak fault slip rates <100 μm/s. The peak amplitude of the high‐frequency time‐domain signals scales systematically with fault slip velocity. Stable sliding experiments further support the connection between fault slip rate and high‐frequency radiation. Experiments demonstrate that the origin of the high‐frequency energy is fundamentally linked to changes in fault slip rate, shear strain, and breaking of contact junctions within the fault gouge. Our results suggest that having measurements close to the fault zone may be key for documenting seismic radiation properties and fully understanding the connection between different slip modes. 
    more » « less
  4. Abstract Imaging tectonic creep along active faults is critical for measuring strain accumulation and ultimately understanding the physical processes that guide creep and the potential for seismicity. We image tectonic deformation along the central creeping section of the San Andreas Fault at the Dry Lake Valley paleoseismic site (36.468°N, 121.055°W) using three data sets with varying spatial and temporal scales: (1) an Interferometric Synthetic Aperture Radar (InSAR) velocity field with an ~100‐km footprint produced from Sentinel‐1 satellite imagery, (2) light detection and ranging (lidar) and structure‐from‐motion 3‐D topographic differencing that resolves a decade of deformation over a 1‐km aperture, and (3) surface fractures that formed over the 3‐ to 4‐m wide fault zone during a drought from late 2012 to 2014. The InSAR velocity map shows that shallow deformation is localized to the San Andreas Fault. We demonstrate a novel approach for differencing airborne lidar and structure‐from‐motion topography that facilitates resolving deformation along and adjacent to the San Andreas Fault. The 40‐m resolution topographic differencing resolves a 2.5 ± 0.2 cm/yr slip rate localized to the fault. The opening‐mode fractures accommodate cm/yr of fault slip. A 90% ± 30% of the 1‐km aperture deformation is accommodated over the several meter‐wide surface trace of the San Andreas Fault. The extension direction inferred from the opening‐mode fractures and topographic differencing is 40°–48° from the local trend of the San Andreas Fault. The localization of deformation likely reflects the well‐oriented and mature fault. 
    more » « less
  5. Abstract Localized frictional sliding on faults in the continental crust transitions at depth to distributed deformation in viscous shear zones. This brittle‐ductile transition (BDT), and/or the transition from velocity‐weakening (VW) to velocity‐strengthening (VS) friction, are controlled by the lithospheric thermal structure and composition. Here, we investigate these transitions, and their effect on the depth extent of earthquakes, using 2D antiplane shear simulations of a strike‐slip fault with rate‐and‐state friction. The off‐fault material is viscoelastic, with temperature‐dependent dislocation creep. We solve the heat equation for temperature, accounting for frictional and viscous shear heating that creates a thermal anomaly relative to the ambient geotherm which reduces viscosity and facilitates viscous flow. We explore several geotherms and effective normal stress distributions (by changing pore pressure), quantifying the thermal anomaly, seismic and aseismic slip, and the transition from frictional sliding to viscous flow. The thermal anomaly can reach several hundred degrees below the seismogenic zone in models with hydrostatic pressure but is smaller for higher pressure (and these high‐pressure models are most consistent with San Andreas Fault heat flow constraints). Shear heating raises the BDT, sometimes to where it limits rupture depth rather than the frictional VW‐to‐VS transition. Our thermomechanical modeling framework can be used to evaluate lithospheric rheology and thermal models through predictions of earthquake ruptures, postseismic and interseismic crustal deformation, heat flow, and the geological structures that reflect the complex deformation beneath faults. 
    more » « less