skip to main content


Title: Imminent seismic velocity changes in earthquake cycle simulations
Earthquake prediction is the holy grail of seismology. Many previous studies have searched for robust precursory signals to inform us of imminent earthquakes, the most significant of which are seen in laboratory experiments as temporal changes in pressure and shear wave velocities during the seismic cycle. Similar changes are seen in natural faults and the surrounding structurally complex network of fractures with nested hierarchy of localized deformation, referred to as fault damage zone. However, little is known whether such temporal changes in material properties contains any precursory signals for imminent earthquakes.Conversely, the effect of precursory velocity changes on the seismic cycle is not well understood. By imposing shear wave velocity changes in fault damage zones, we investigate the effects of these precursors on multiple stages of the seismic cycle, including nucleation, coseismic, postseismic, and interseismic stages. We perform 2D fully dynamic earthquake cycle simulations with a fault-parallel damage zone for strike-slip fault systems with antiplane geometry. The fault is governed by rate-state-dependent friction laws, and the fault damage zone material is considered elastic. Our preliminary results show that the temporal onset of shear wave velocity drop causes a reduction in earthquake recurrence intervals over the seismic cycle. Furthermore, a dynamic earthquake rupture within the seismic cycle terminates much faster and abruptly in models with precursory velocity changes. We will also discuss how the precursory velocity changes affect the fault-slip behavior, including fast-slip, slow-slip, and aseismic creep, for different amplitudes of shear wave velocity changes at different compliance contrast of the fault damage zones. Our results highlight the importance of short and long-term monitoring of fault zone structures for better assessment of regional seismic hazard.  more » « less
Award ID(s):
1943742
NSF-PAR ID:
10299110
Author(s) / Creator(s):
;
Date Published:
Journal Name:
AGU Fall Meeting 2021
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Predicting the onset and timing of fault failure is one of the ultimate goals of seismology. However, our current understanding of the earthquake preparation and nucleation process is limited. One direction towards understanding this process is looking at precursory signals preceding large earthquakes. Previous laboratory experiments have studied robust precursory signals, observed as temporal changes in pressure and shear wave velocities during the seismic cycle. The effects of such precursory velocity changes on the seismic cycle are not well understood. We use numerical models to simulate fully-dynamic earthquake cycles in 2D strike-slip fault systems with antiplane geometry, surrounded by a narrow fault-parallel damage zone. By imposing shear wave velocity changes inside fault damage zones, we investigate the effects of these precursors on multiple stages of the seismic cycle, including nucleation, coseismic, postseismic, and interseismic stages. Our modeling results show a wide spectrum of fault-slip behaviors including fast earthquakes, slow-slip events, and variable creep. One primary effect of the imposed velocity precursor is the facilitation of the otherwise slow-slip event to grow into a fully dynamic earthquake. Furthermore, the onset time of these precursors have significant effects on the nucleation phase of the earthquakes, and earlier onset of precursors causes the earthquakes to nucleate earlier with a smaller nucleation size. Our results highlight the importance of short and long-term monitoring of fault zone structures for better assessment of regional seismic hazard. 
    more » « less
  2. null (Ed.)
    The fault damage zone is a well-known structure of localized deformation around faults. Its material properties evolve over earthquake cycles due to coseismic damage accumulation and interseismic healing. We will present fully dynamic earthquake cycle simulations to show how the styles of earthquake nucleation and rupture propagation change as fault zone material properties vary temporally. First, we will focus on the influence of fault zone structural maturity quantified by near-fault seismic wave velocities in simulations. The simulations show that immature fault zones promote small and moderate subsurface earthquakes with irregular recurrence intervals, whereas mature fault zones host pulse-like earthquake rupture that can propagate to the surface, extend throughout the seismogenic zone, and occur at regular intervals. The interseismic healing in immature fault zones plays a key role in allowing the development of aseismic slip episodes including slow-slip events and creep, which can propagate into the seismogenic zone, and thus limit the sizes of subsequent earthquakes by releasing fault stress. In the second part, we will discuss how the precursory changes of seismic wave velocities of fault damage zones may affect earthquake nucleation process. Both laboratory experiments and seismic observations show that the abrupt earthquake failure can be preceded by accelerated fault deformation and the accompanying velocity reduction of near-fault rocks. We will use earthquake cycle simulations to systematically test the effects of timing and amplitudes of such precursory velocity changes. Our simulations will provide new insights into the interplay between fault zone structure and earthquake nucleation process, which can be used to guide future real-time monitoring of major fault zones. 
    more » « less
  3. 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
  4. Abstract

    Mature strike‐slip faults are usually surrounded by a narrow zone of damaged rocks characterized by low seismic wave velocities. Observations of earthquakes along such faults indicate that seismicity is highly concentrated within this fault damage zone. However, the long‐term influence of the fault damage zone on complete earthquake cycles, that is, years to centuries, is not well understood. We simulate aseismic slip and dynamic earthquake rupture on a vertical strike‐slip fault surrounded by a fault damage zone for a thousand‐year timescale using fault zone material properties and geometries motivated by observations along major strike‐slip faults. The fault damage zone is approximated asan elastic layer with lower shear wave velocity than the surrounding rock. We find that dynamic wave reflections, whose characteristics are strongly dependent on the width and the rigidity contrast of the fault damage zone, have a prominent effect on the stressing history of the fault. The presence of elastic damage can partially explain the variability in the earthquake sizes and hypocenter locations along a single fault, which vary with fault damage zone depth, width and rigidity contrast from the host rock. The depth extent of the fault damage zone has a pronounced effect on the earthquake hypocenter locations, and shallower fault damage zones favor shallower hypocenters with a bimodal distribution of seismicity along depth. Our findings also suggest significant effects on the hypocenter distribution when the fault damage zone penetrates to the nucleation sites of earthquakes, likely being influenced by both lithological (material) and rheological (frictional) boundaries.

     
    more » « less
  5. Abstract

    Tectonic faults fail in a continuum of modes from slow earthquakes to elastodynamic rupture. Precursory variations in elastic wavespeed and amplitude, interpreted as indicators of imminent failure, have been observed in limited natural settings and lab experiments where they are thought to arise from contact rejuvenation and microcracking within and around the fault zone. However, the physical mechanisms and connections to fault creep are poorly understood. Here we vary loading stiffness during frictional shear to generate a range of slip modes and measure fault zone properties using transmitted elastic waves. We find that elastic wave amplitudes show clear changes before fault failure. The temporal onset of amplitude reduction scales with lab earthquake magnitude and the magnitude of this reduction varies with fault slip. Our data provide clear evidence of precursors to lab earthquakes and suggest that continuous seismic monitoring could be useful for assessing fault state and seismic hazard potential.

     
    more » « less