Slow Slip Transients Research Articles

Parallel dynamics of slow slips and fluid-induced seismic swarms

Earthquake swarms may be driven by fluids, through hydraulic injections or natural fluid circulation, but also by slow and aseismic slip transients. Understanding the driving factors for these prolific sequences and how they can potentially develop into larger ruptures remains a challenge. A notable and almost ubiquitous feature of swarms is their hypocenters migration, which occurrence is closely related to the processes driving the observed seismicity, in a similar way as seismicity accompanies slow-slip events at subduction zones. Here, we analyze global data on migrating sequences, and identify scaling laws for migration velocity, moment and duration measured on natural and injection-induced swarms, foreshock sequences, and slow slip events. We highlight two different behaviors among these sequences: one linked to slow slips, with elevated migration velocities and moments, and the other related to fluid-induced processes, featuring lower velocities and moments. These results provide metrics for distinguishing between the drivers of earthquake swarms, fluid or slow-slip related, and prompt a reevaluation of scaling laws of fault slip transients, especially for swarms.

Fault rock heterogeneity can produce fault weakness and reduce fault stability

Geological heterogeneity is abundant in crustal fault zones; however, its role in controlling the mechanical behaviour of faults is poorly constrained. Here, we present laboratory friction experiments on laterally heterogeneous faults, with patches of strong, rate-weakening quartz gouge and weak, rate-strengthening clay gouge. The experiments show that the heterogeneity leads to a significant reduction in strength and frictional stability in comparison to compositionally identical faults with homogeneously mixed gouges. We identify a combination of weakening effects, including smearing of the weak clay; differential compaction of the two gouges redistributing normal stress; and shear localization producing stress concentrations in the strong quartz patches. The results demonstrate that geological heterogeneity and its evolution can have pronounced effects on fault strength and stability and, by extension, on the occurrence of slow-slip transients versus earthquake ruptures and the characteristics of the resulting events, and should be further studied in lab experiments and earthquake source modelling.

Cataclastic and crystal-plastic deformation in shallow mantle-wedge serpentinite controlled by cyclic changes in pore fluid pressures

Episodic tremor and slip (ETS) events in the forearc mantle wedge of a warm subduction zone may reflect mixed brittle–ductile deformation of serpentinite in association with high pore fluid pressures. To understand deformation mechanisms and processes occurring in the hydrated mantle wedge, we examined in the Sanbagawa metamorphic belt, SW Japan, an antigorite serpentinite shear zone derived from mantle wedge that was formed under pressure and temperature conditions that correspond to the ETS regions. The serpentinite underwent multiple extensional (mode I) and extensional–shear (mode I–II) failure events at supralithostatic pore fluid pressures (Pf). Such failure events led to drops in Pf (several MPa) and formation of a distributed ‘fault–fracture mesh’. Antigorite precipitation in the fracture openings contributed to an increase in Pf until the failure condition was reached again, and thereby antigorite kinetics controlled the recurrence interval of seismic events. We also suggest that under the low-Pf conditions that facilitate intracrystalline plasticity rather than cataclasis, the newly precipitated antigorite aggregates (localized along shear bands) deform by dislocation creep at a high strain rate and high shear stress, resulting in the transient, accelerated viscous creep that may characterize slow slip transients.

The generation of large earthquakes

Despite decades of observational, laboratory and theoretical studies, the processes leading to large earthquake generation remain enigmatic. However, recent observations provide new promising perspectives that advance knowledge. Here, we review data on the initiation processes of large earthquakes and show that they are multiscale and diverse, involving localization of deformation, fault heterogeneities and variable local loading rate effects. Analyses of seismic and geodetic data reveal evidence for regional weakening by earthquake-induced rock damage and progressive localization of deformation around the eventual rupture zones a few years before some large earthquakes. The final phase of deformation localization includes, depending on conditions, a mixture of slow slip transients and foreshocks at multiple spatial and temporal scales. The evolution of slip on large, localized faults shows a step-like increase that might reflect stress loading by previous failures, which can produce larger dynamic slip, in contrast to the smooth acceleration expected for a growing aseismic nucleation phase. We propose an integrated model to explain the diversity of large earthquake generation from progressive volumetric deformation to localized slip, which motivates future near-fault seismic and geodetic studies with dense sensor networks and improved analysis techniques that can resolve multiscale processes. The processes leading to large earthquakes remain enigmatic. Using detailed seismic and geodetic data, this Review examines how tectonic deformation and evolving fault behaviour initiate large earthquakes, and proposes an integrated model accounting for the diversity of observations.

Excitation of San Andreas tremors by thermal instabilities below the seismogenic zone.

The relative motion of tectonic plates is accommodated at boundary faults through slow and fast ruptures that encompass a wide range of source properties. Near the Parkfield segment of the San Andreas fault, low-frequency earthquakes and slow-slip events take place deeper than most seismicity, at temperature conditions typically associated with stable sliding. However, laboratory experiments indicate that the strength of granitic gouge decreases with increasing temperature above 350°C, providing a possible mechanism for weakening if temperature is to vary dynamically. Here, we argue that recurring low-frequency earthquakes and slow-slip transients at these depths may arise because of shear heating and the temperature dependence of frictional resistance. Recurring thermal instabilities can explain the recurrence pattern of the mid-crustal low-frequency earthquakes and their correlative slip distribution. Shear heating associated with slow slip is sufficient to generate pseudotachylyte veins in host rocks even when fault slip is dominantly aseismic.

Daily measurement of slow slip from low-frequency earthquakes is consistent with ordinary earthquake scaling.

Slow slip transients on faults can last from seconds to months and stitch together the earthquake cycle. However, no single geophysical instrument is able to observe the full range of slow slip because of bandwidth limitations. Here, we connect seismic and geodetic data from the Mexican subduction zone to explore an instrumental blind spot. We establish a calibration of the daily median amplitude of the seismically recorded low-frequency earthquakes to the daily geodetically recorded moment rate of previously established slow slip events. This calibration allows us to use the precise evolution of low-frequency earthquake activity to quantitatively measure the moment of smaller, subdaily slip events that are unresolvable by geodesy alone. The resulting inferred slow slip moments scale with duration and inter-event time like ordinary earthquakes. These new quantifications help connect slow and fast events in a broad spectrum of transient slip and suggest that slow slip events behave much like ordinary earthquakes.

Strength of tremor patches along deep transition zone of a megathrust

Deep low frequency tremors are indicators of slow slip transients in the brittle-ductile transition zone along subducting plates. Investigation of comprehensive tremor activities is therefore an important issue for understanding the seismic/aseismic characteristics in transition zones. Here, we focus on the radiated energy from tremors to reveal the along-strike heterogeneity in the strength of tremor patches. Based on a tremor catalog that more accurately evaluates radiated energy, we examine the spatio-temporal activity of tremors accompanied by slow slip events (SSEs) in western Shikoku, southwestern Japan. The new finding of this study is that the energy radiated from tremors is positively correlated with the speed of tremor migration front and the slip rate along the plate interface during a SSE. This can be qualitatively explained by a stress diffusion model, which consists of along-strike heterogeneities in the effective strength of tremor patches embedded in a ductile shear zone. This effective strength heterogeneity is supported by a lateral variation in the stress drop of a SSE; it is consistent with the fluid pressure distribution along the plate boundary fault and the tidal sensitivity of tremors. Accurate evaluation of tremor activities, especially the radiated energy, can be used to infer the spatial distribution of the strength of tremor patches in transition zones worldwide.

Imbricated slip rate processes during slow slip transients imaged by low-frequency earthquakes

Low Frequency Earthquakes (LFEs) often occur in conjunction with transient strain episodes, or Slow Slip Events (SSEs), in subduction zones. Their focal mechanism and location consistent with shear failure on the plate interface argue for a model where LFEs are discrete dynamic ruptures in an otherwise slowly slipping interface. SSEs are mostly observed by surface geodetic instruments with limited resolution and it is likely that only the largest ones are detected. The time synchronization of LFEs and SSEs suggests that we could use the recorded LFEs to constrain the evolution of SSEs, and notably of the geodetically-undetected small ones. However, inferring slow slip rate from the temporal evolution of LFE activity is complicated by the strong temporal clustering of LFEs. Here we apply dedicated statistical tools to retrieve the temporal evolution of SSE slip rates from the time history of LFE occurrences in two subduction zones, Mexico and Cascadia, and in the deep portion of the San Andreas fault at Parkfield. We find temporal characteristics of LFEs that are similar across these three different regions. The longer term episodic slip transients present in these datasets show a slip rate decay with time after the passage of the SSE front possibly as t−1/4. They are composed of multiple short term transients with steeper slip rate decay as t−α with α between 1.4 and 2. We also find that the maximum slip rate of SSEs has a continuous distribution. Our results indicate that creeping faults host intermittent deformation at various scales resulting from the imbricated occurrence of numerous slow slip events of various amplitudes.

Slow Slip Transients Before the 2011 Tohoku-Oki Earthquake

I review a spatiotemporal evolution of slow-slip transients on the plate interface of the subducting Pacific plate that happened in and around the mainshock rupture area prior to the 2011 Tohoku-Oki earthquake. Based on foreshock activity before the mainshock, two sequences of slow-slip transients were identified by earthquake migrations toward the initiation point of the mainshock rupture. These two slow-slip transients were also detected by geodetic measurement. The second sequence of slow-slip transients, which involved large slip rates, may have caused significant stress loading onto the hypocenter of the mainshock and prompted the initiation of unstable dynamic rupture. In addition, decadal slip-behavior on the plate interface revealed by geodetic measurement and small repeating earthquakes show that slow-slip transients occurred in the down-dip and up-dip portions of the mainshock rupture area. These slow-slip transients imply the reduction of coupling between the subducting and overlying plates, that could be interpreted as the late stage of mega-thrust earthquake cycle, although this notion remains conjectural.

Seismic and aseismic slip pulses driven by thermal pressurization of pore fluid

There are several lines of evidence that suggest that thermal pressurization (TP) of pore fluid within a low‐permeability fault core may play the key role in the development of earthquake slip. To elucidate effects of TP on spontaneous fault slip, we consider solutions for a steadily propagating slip pulse on a fault with a constant sliding friction, the level of which may reflect other thermally‐activated processes at the rupture front (such as the flash heating on asperities). Upon arrival of the pulse front, essentially undrained‐adiabatic TP takes place during the initial slip acceleration from the locked state with a corresponding reduction of the fault strength. With passage of time, the diminishing rate of heating (due to the reduced fault strength) and increasing rate of hydrothermal diffusion from the shear zone offset TP and result in partial recovery of the strength, slip deceleration and eventual locking and healing of the slip. We show that the rupture speedvr decreases with thickness hof the principal shear zone. For lab‐constrained values of fault‐gouge parameters, the TP‐pulse solution predicts seismic (vr∼ km/s) slip on a millimeter‐to‐cm thin principal shear zone; and aseismic slip withvr ∼ 10 km/day and slip rates 1–2 orders above the plate rate on a relatively thick (h∼ 1 m) shear zone. These and other predictions of the TP‐pulse model are consistent with the independent sets of observational constraints for large crustal and subduction interplate earthquakes, and slow slip transients (North Cascadia), respectively. Locking of the slip soon after the diffusive transport of the heat and pore fluid becomes efficient significantly limits the maximum co‐seismic temperature rise to values well below previous theoretical estimates. As a result, the onset of macroscopic melting and some of thermal decomposition reactions, recently suggested to explain strong co‐seismic fault weakening, are precluded over much of the seismogenic zone.

Propagation of Slow Slip Leading Up to the 2011 M w 9.0 Tohoku-Oki Earthquake

Many large earthquakes are preceded by one or more foreshocks, but it is unclear how these foreshocks relate to the nucleation process of the mainshock. On the basis of an earthquake catalog created using a waveform correlation technique, we identified two distinct sequences of foreshocks migrating at rates of 2 to 10 kilometers per day along the trench axis toward the epicenter of the 2011 moment magnitude (M(w)) 9.0 Tohoku-Oki earthquake in Japan. The time history of quasi-static slip along the plate interface, based on small repeating earthquakes that were part of the migrating seismicity, suggests that two sequences involved slow-slip transients propagating toward the initial rupture point. The second sequence, which involved large slip rates, may have caused substantial stress loading, prompting the unstable dynamic rupture of the mainshock.

Slow slip predictions based on granite and gabbro friction data compared to GPS measurements in northern Cascadia

For episodic slow slip transients in subduction zones, a large uncertainty in comparing surface deformations predicted by forward modeling based on rate and state friction to GPS measurements lies in our limited knowledge of the frictional properties and fluid pore pressure along shallow subduction faults. In this study, we apply the laboratory rate and state friction data of granite and gabbro gouges under hydrothermal conditions to a Cascadia‐like 2‐D model to produce spontaneous aseismic transients, and we compare the resulting intertransient and transient surface deformations to GPS observations along the northern Cascadia margin. An inferred region along dip of elevated fluid pressure is constrained by seismological observations where available and by thermal and petrological models for the Cascadia and SW Japan subduction zones. For the assumed friction parameters a and a − b profiles, we search the model parameter space, by varying the level of effective normal stress , characteristic slip distance L in the source areas of transients, and the fault width under , to identify simulation cases that produce transients of total aseismic slip and recurrence interval similar to the observed 20–30 mm and 1–2 years, respectively, in northern Cascadia. Using a simple planar fault geometry and extrapolating the 2‐D fault slip to a 3‐D distribution, we find that the friction data for gabbro gouge, a better representation of the seafloor, fit GPS observations of transient deformation in northern Cascadia much better than do the granite data, which, for lack of a suitable alternative, have been the basis for most previous modeling.

Slow slip transients along the Oaxaca subduction segment from 1993 to 2007

SUMMARY We use data from 12 continuous GPS stations in southern Mexico, including eight new stations, to better characterize transient slip episodes along the Mexican subduction zone. Continuous GPS recording in Oaxaca that began 14 years ago, constituting the longest continuous GPS record in southern Mexico, defines nine distinct episodes of transient slip from 1993 to 2007, including previously unreported transient slip episodes in early 1995 and 2006. All transient slip episodes recorded in Oaxaca City were also recorded at a GPS site ∼400 km away in the state of Guerrero after measurements began there in early 1997, demonstrating that transient slip affects widespread areas of southern Mexico. Well-recorded transients in 2004 and 2006 appear to have originated along the Oaxaca trench segment, hundreds of kilometres from the potentially hazardous Guerrero seismic gap. During the 2004 slip transient, displacements of more than 20 mm were recorded at stations in the Oaxaca GPS array, a factor of 2‐3 larger than at sites in Guerrero. Elastic half-space modelling of the 2004 displacements at sites throughout southern Mexico indicates that most slip was focused beneath Oaxaca in deeper areas of the subduction interface immediately downdip from previous megathrust rupture zones, but that slip also may have extended significantly updip into a seismic gap along the Oaxaca coast. The best-fitting model is also able to explain transient vertical offsets of up to 30 mm observed during this event. Shortly before the onset of transient slip in 2004, a pair of moderate (M w 5.1, 5.5) earthquakes ruptured the downdip end of the seismogenic zone immediately east of the Oaxaca seismic gap, suggesting a possible relationship between the two. In 2006, transient slip began 1‐3 months earlier in Oaxaca than in other areas, offering the clearest evidence to date for the existence of one or more source regions for transient slip outside the Guerrero seismic gap, where large amplitude transient slip originated in 2002. The still-limited data indicate that transient slip that originates elsewhere in Mexico may trigger aseismic slip in Guerrero, and further indicate that transient slip beneath Oaxaca is limited to areas of the subduction interface that surround the rupture zones of previous large shallow-thrust earthquakes.

Foreshock sequences and short-term earthquake predictability on East Pacific Rise transform faults

East Pacific Rise transform faults are characterized by high slip rates (more than ten centimetres a year), predominantly aseismic slip and maximum earthquake magnitudes of about 6.5. Using recordings from a hydroacoustic array deployed by the National Oceanic and Atmospheric Administration, we show here that East Pacific Rise transform faults also have a low number of aftershocks and high foreshock rates compared to continental strike-slip faults. The high ratio of foreshocks to aftershocks implies that such transform-fault seismicity cannot be explained by seismic triggering models in which there is no fundamental distinction between foreshocks, mainshocks and aftershocks. The foreshock sequences on East Pacific Rise transform faults can be used to predict (retrospectively) earthquakes of magnitude 5.4 or greater, in narrow spatial and temporal windows and with a high probability gain. The predictability of such transform earthquakes is consistent with a model in which slow slip transients trigger earthquakes, enrich their low-frequency radiation and accommodate much of the aseismic plate motion.

Imaging of aseismic fault slip transients recorded by dense geodetic networks

SUMMARY The large-scale continuous GPS networks that have emerged in the last decade have documented a wide-range of aseismic deformation transients that resulted from physical processes such as aseismic fault slip and magma intrusion. In particular, a new class of (M ∼7) slow earthquakes with durations ranging from days to months have been observed with GPS arrays located above the downdip portion of subduction zone thrust interfaces. Interpretation of the displacement time-series resulting from these events is not straightforward owing to the contaminating effects of multiple contributing signals such as fault-slip, local benchmark motion, seasonal effects, and reference frame errors. We have developed a time-dependent inversion algorithm based on the extended Kalman filter which can separate the various signals and allow the space‐ time evolution of these slow-slip transients to be studied in detail. We applied the inversion algorithm to the 1999 Cascadia slow earthquake. This event had two primary episodes of moment-release separated by a two week period in which relatively little moment-release occurred. The Cascadia event and other slow earthquakes share numerous similarities with both ordinary earthquakes and afterslip transients suggesting that they may represent fault slip under a conditionally stable regime.