skip to main content

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


  • NSF Public Access
  • Search Results
  • Inferring damage state and evolution with increasing stress using direct and coda wave velocity measurements in faulted and intact granite samples
Title: Inferring damage state and evolution with increasing stress using direct and coda wave velocity measurements in faulted and intact granite samples
SUMMARY

A better understanding of damage accumulation before dynamic failure events in geological material is essential to improve seismic hazard assessment. Previous research has demonstrated the sensitivity of seismic velocities to variations in crack geometry, with established evidence indicating that initial crack closure induces rapid changes in velocity. Our study extends these findings by investigating velocity changes by applying coda wave interferometry (CWI). We use an array of 16 piezoceramic transducers to send and record ultrasonic pulses and to determine changes in seismic velocity on intact and faulted Westerly granite samples. Velocity changes are determined from CWI and direct phase arrivals. This study consists of three sets of experiments designed to characterize variations in seismic velocity under various initial and boundary conditions. The first set of experiments tracks velocity changes during hydrostatic compression from 2 and 191 MPa in intact Westerly granite samples. The second set of experiments focuses on saw-cut samples with different roughness and examines the effects of confining pressure increase from 2 to 120 MPa. The dynamic formation of a fracture and the preceding damage accumulation is the focus of the third type of experiment, during which we fractured an initially intact rock sample by increasing the differential stress up to 780 MPa while keeping the sample confined at 75 MPa. The tests show that: (i) The velocity change for rough saw cut samples suggests that the changes in bulk material properties have a more pronounced influence than fault surface apertures or roughness. (ii) Seismic velocities demonstrate higher sensitivity to damage accumulation under increasing differential stress than macroscopic measurements. Axial stress measured by an external load cell deviates from linearity around two-third through the experiment at a stress level of 290 MPa higher than during the initial drop in seismic velocities. (iii) Direct waves exhibit strong anisotropy with increasing differential stress and accumulating damage before rock fracture. Coda waves, on the other hand, effectively average over elastic wave propagation for both fast and slow directions, and the resulting velocity estimates show little evidence for anisotropy. The results demonstrate the sensitivity of seismic velocity to damage evolution at various boundary conditions and progressive microcrack generation with long lead times before dynamic fracture.

  more » « less
Award ID(s):
2142489
PAR ID:
10471926
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Geophysical Journal International
Volume:
235
Issue:
3
ISSN:
0956-540X
Format(s):
Medium: X Size: p. 2846-2861
Size(s):
p. 2846-2861
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ( , Geophysical Journal International)
    null (Ed.)
    SUMMARY Earthquake ruptures are generally considered to be cracks that propagate as fracture or frictional slip on pre-existing faults. Crack models have been used to describe the spatial distribution of fault offset and the associated static stress changes along a fault, and have implications for friction evolution and the underlying physics of rupture processes. However, field measurements that could help refine idealized crack models are rare. Here, we describe large-scale laboratory earthquake experiments, where all rupture processes were contained within a 3-m long saw-cut granite fault, and we propose an analytical crack model that fits our measurements. Similar to natural earthquakes, laboratory measurements show coseismic slip that gradually tapers near the rupture tips. Measured stress changes show roughly constant stress drop in the centre of the ruptured region, a maximum stress increase near the rupture tips and a smooth transition in between, in a region we describe as the earthquake arrest zone. The proposed model generalizes the widely used elliptical crack model by adding gradually tapered slip at the ends of the rupture. Different from the cohesive zone described by fracture mechanics, we propose that the transition in stress changes and the corresponding linear taper observed in the earthquake arrest zone are the result of rupture termination conditions primarily controlled by the initial stress distribution. It is the heterogeneous initial stress distribution that controls the arrest of laboratory earthquakes, and the features of static stress changes. We also performed dynamic rupture simulations that confirm how arrest conditions can affect slip taper and static stress changes. If applicable to larger natural earthquakes, this distinction between an earthquake arrest zone (that depends on stress conditions) and a cohesive zone (that depends primarily on strength evolution) has important implications for how seismic observations of earthquake fracture energy should be interpreted. 
    more » « less
  2. ; ; ( , Geophysical Research Letters)
    Abstract

    Earthquakes rarely occur in isolation but rather as complex sequences of fore, main and aftershocks. Assessing the associated seismic hazard requires a holistic view of event interactions. We conduct frictional sliding experiments on faulted Westerly Granite samples at mid‐crustal stresses to investigate fault damage and roughness effects on aftershock generation. Abrupt laboratory fault slip is followed by periods of extended stress relaxation and aftershocks. Large roughness promotes less co‐seismic slip and high aftershock activity whereas smooth faults promote high co‐seismic slip with few aftershocks. Conditions close to slip instability generate lab‐quake sequences that exhibit similar statistical distributions to natural earthquakes. Aftershock productivity in the lab is linearly related to the residual strain energy on the fault which, in turn, is controlled by the level of surface heterogeneity. We conclude that roughness and damage govern slip stability and seismic energy partitioning between fore, main and aftershocks in lab and nature.

     
    more » « less
  3. ; ( , Journal of Geophysical Research: Solid Earth)
    Abstract

    Brittle fracture propagation in rocks is a complex process due to significant grain‐scale heterogeneity and evolving stress states under dynamic loading conditions. In this work, we use digital image correlation and linear elastic fracture mechanics to make instantaneous measurements of the opening (mode I) and in plane shear (mode II) components of the stress intensity field during dynamic mixed mode crack initiation and propagation in crystalline and granular rocks. Both rock types display some similar fracture behaviors as observed in engineered materials, including rate dependent fracture initiation toughness and a direct relationship between propagation toughness and crack velocity; however, measured propagation toughness is higher than quasi‐static values at crack velocities well below the branching velocity in both rocks. Additionally, due to grain scale controls on the fracture process, mixed mode crack propagation is fundamentally different between these two rock types. Mixed mode propagation is energetically more favorable than pure opening mode propagation in sandstone, while the opposite is true in granite. Furthermore, following initiation, propagation in granite occurs so as to minimize the mode II contribution, irrespective of the initiation conditions, while fractures in sandstone maintain a non‐negligible mode II contribution during propagation across the sample.

     
    more » « less
  4. ; ( , Geophysical Research Letters)
    Abstract

    Large strike‐slip faults experience numerous earthquakes during which transient tensile and compressive mean normal stress perturbations travel along opposing sides of the fault. Research exploring dynamic rock fracture through multiple earthquake cycles has focused predominantly on transient compressive loading, but little is known about off‐fault damage development due to successive tensile loading. We investigate damage accumulation by transient tensile loading over multiple earthquake cycles using a modified sample configuration for uniaxial compressive loading apparatuses consisting of a Westerly granite rock disk bonded to two lead disks. We show that fracture density increases during each successive loading cycle, and pulverized rock can be produced under tension at strain rates as low as 10−3s−1. Therefore, pulverized rock can form at low strain rates, and its texture and extent may be controlled by the size of the coseismic tensile stress perturbation and the number of slip events on the fault.

     
    more » « less
  5. ; ( , 2023 SSA Annual Meeting)
    The temporal variation of elastic property of the bulk material surrounding the fault is considered an important contribution to the observed co-seismic velocity reduction and interseismic healing. Paglialunga et al. [2021] found that as fault normal stress increases, co-seismic velocity reduction becomes larger because more cracks reopen with higher stress drops. Larger normal stress can lead to smaller nucleation size and contribute to larger co-seismic slip. By contrast, with larger co-seismic velocity reduction and interseismic healing, more slow slip events can propagate in the seismogenic zone [Thakur and Huang, 2021], because the temporal velocity change related to fault zone damage modulates earthquake nucleation. Hence, fault normal stress and temporal damage zone structure evolution have opposite influences on the spatial distribution and recurrence intervals of earthquakes. We conducted 2-D anti-plane fully-dynamic seismic cycle simulations and explored the effects of fault normal stress on seismic cycle when there is coseismic damage and interseismic healing in the fault damage zone. The normal stress is in a range of 40-70 MPa and the co-seismic rigidity reduction is in a range of 5-8%. We find larger normal stress results in larger co-seismic slip and fewer slow slip events, while more co-seismic velocity reduction and interseismic healing leads to more partial ruptures as well as slow slip events. With the increase of both normal stress and seismic velocity change, more regular earthquakes occur and slow slip events gradually disappear. For the selected parameter space, the influence of seismic velocity change is not as significant as the effect of normal stress. However, fault zone maturity or the initial rigidity of fault damage zones should also affect the competitive relationship between normal stress and seismic velocity change, and we will characterize earthquakes and slow-slip events in immature and mature fault damage zones when both on-fault normal stress and off-fault seismic velocity vary over earthquake cycles. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email