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 predictable chaos of slow earthquakes
Slow earthquakes, like regular earthquakes, result from unstable frictional slip. They produce little slip and can therefore repeat frequently. We assess their predictability using the slip history of the Cascadia subduction between 2007 and 2017, during which slow earthquakes have repeatedly ruptured multiple segments. We characterize the system dynamics using embedding theory and extreme value theory. The analysis reveals a low-dimensional (<5) nonlinear chaotic system rather than a stochastic system. We calculate properties of the underlying attractor like its correlation and instantaneous dimension, instantaneous persistence, and metric entropy. We infer that the system has a predictability horizon of the order of days weeks. For the better resolved segments, the onset of large slip events can be correctly forecasted by high values of the instantaneous dimension. Longer-term deterministic prediction seems intrinsically impossible. Regular earthquakes might similarly be predictable but with a limited predictable horizon of the order of their durations.  more » « less
Award ID(s):
1821853
PAR ID:
10291993
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Science Advances
Volume:
6
Issue:
27
ISSN:
2375-2548
Page Range / eLocation ID:
eaaz5548
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We study the mechanical response of two‐dimensional vertical strike‐slip fault to coseismic damage evolution and interseismic healing of fault damage zones by simulating fully dynamic earthquake cycles. Our models show that fault zone structure evolution during the seismic cycle can have pronounced effects on mechanical behavior of locked and creeping fault segments. Immature fault damage zone models exhibit small and moderate subsurface earthquakes with irregular recurrence intervals and abundance of slow‐slip events during the interseismic period. In contrast, mature fault damage zone models host pulse‐like earthquake ruptures that can propagate to the surface and extend throughout the seismogenic zone, resulting in large stress drop, characteristic rupture extents, and regular recurrence intervals. Our results suggest that interseismic healing and coseismic damage accumulation in fault zones can explain the observed differences of earthquake behaviors between mature and immature fault zones and indicate a link between regional seismic hazard and fault structural maturity. 
    more » « less
  2. Recent seismic and geodetic observations indicate that interseismic creep rate varies in both time and space. The spatial extent of creep pinpoints locked asperities, while its temporary accelerations, known as slow-slip events, may trigger earthquakes. Although the conditions promoting fault creep are well-studied, the mechanisms for initiating episodic slow-slip events are enigmatic. Here we investigate surface deformation measured by radar interferometry along the central San Andreas Fault between 2003 and 2010 to constrain the temporal evolution of creep. We show that slow-slip events are ensembles of localized creep bursts that aseismically rupture isolated fault compartments. Using a rate-and-state friction model, we show that effective normal stress is temporally variable on the fault, and support this using seismic observations. We propose that compaction-driven elevated pore fluid pressure in the hydraulically isolated fault zone and subsequent frictional dilation cause the observed slow-slip episodes. We further suggest that the 2004 Mw 6 Parkfield earthquake might have been triggered by a slow-slip event, which increased the Coulomb failure stress by up to 0.45 bar per year. This implies that while creeping segments are suggested to act as seismic rupture barriers, slow-slip events on these zones might promote seismicity on adjacent locked segments. 
    more » « less
  3. Abstract Megathrust earthquakes and their associated tsunamis cause some of the worst natural disasters. In addition to earthquakes, a wide range of slip behaviors are present at subduction zones, including slow earthquakes that span multiple orders of spatial and temporal scales. Understanding these events may shed light on the stress or strength conditions of the megathrust fault. Out of all types of slow earthquakes, very low frequency earthquakes (VLFEs) are most enigmatic because they are difficult to detect reliably, and the physical nature of VLFEs are poorly understood. Here we show three VLFEs in Cascadia that were dynamically triggered by a 2009 Mw 6.9 Canal de Ballenas earthquake in the Gulf of California. The VLFEs likely locate in between the seismogenic zone and the Cascadia episodic tremor and slip (ETS) zone, including one event with a moment magnitude of 5.7. This is the largest VLFE reported to date, causing clear geodetic signals. Our results show that the Cascadia megathrust fault might slip rapidly at some spots in this gap zone, and such a permissible slip behavior has direct seismic hazard implications for coastal communities and perhaps further inland. Further, the observed seismic sources may represent a new class of slip events, whose characteristics do not fit current understandings of slow or regular earthquakes. 
    more » « less
  4. Abstract. Fault slip is a complex natural phenomenon involving multiple spatiotemporal scales from seconds to days to weeks. To understand the physical and chemical processes responsible for the full fault slip spectrum, a multidisciplinary approach is highly recommended. The Near Fault Observatories (NFOs) aim at providing high-precision and spatiotemporally dense multidisciplinary near-fault data, enabling the generation of new original observations and innovative scientific products. The Alto Tiberina Near Fault Observatory is a permanent monitoring infrastructure established around the Alto Tiberina fault (ATF), a 60 km long low-angle normal fault (mean dip 20°), located along a sector of the Northern Apennines (central Italy) undergoing an extension at a rate of about 3 mm yr−1. The presence of repeating earthquakes on the ATF and a steep gradient in crustal velocities measured across the ATF by GNSS stations suggest large and deep (5–12 km) portions of the ATF undergoing aseismic creep. Both laboratory and theoretical studies indicate that any given patch of a fault can creep, nucleate slow earthquakes, and host large earthquakes, as also documented in nature for certain ruptures (e.g., Iquique in 2014, Tōhoku in 2011, and Parkfield in 2004). Nonetheless, how a fault patch switches from one mode of slip to another, as well as the interaction between creep, slow slip, and regular earthquakes, is still poorly documented by near-field observation. With the strainmeter array along the Alto Tiberina fault system (STAR) project, we build a series of six geophysical observatory sites consisting of 80–160 m deep vertical boreholes instrumented with strainmeters and seismometers as well as meteorological and GNSS antennas and additional seismometers at the surface. By covering the portions of the ATF that exhibits repeated earthquakes at shallow depth (above 4 km) with these new observatory sites, we aim to collect unique open-access data to answer fundamental questions about the relationship between creep, slow slip, dynamic earthquake rupture, and tectonic faulting. 
    more » « less
  5. 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