Slow Fault Slip Research Articles

Refined Holocene Slip Rate for the Western and Central Segments of the Garlock Fault: Record of Alternating Millennial-Scale Periods of Fast and Slow Fault Slip

Refined Holocene Slip Rate for the Western and Central Segments of the Garlock Fault: Record of Alternating Millennial-Scale Periods of Fast and Slow Fault Slip

Slow Slip as an Indicator of Fault Stress Criticality

AbstractFault regions inferred to be slowly slipping are interpreted to accommodate much of tectonic plate motion aseismically and potentially serve as barriers to earthquake rupture. Here, we build on prior work using simulations of earthquake sequences with enhanced dynamic fault weakening to show how fault regions that exhibit decades of steady creep or transient slow‐slip events can be driven to dynamically fail by incoming earthquake ruptures. Following substantial earthquake slip, such regions can be under‐stressed and locked for centuries prior to slowly slipping again. Our simulations illustrate that slow fault slip indicates that a region is sufficiently loaded to be failing about its quasi‐static strength. Hence, if a fault region is susceptible to failing dynamically, then observations of slow slip could serve as an indication that the region is critically stressed and ready to fail in a future earthquake, posing a qualitatively different interpretation of slow slip for seismic hazard.

Comparison of statistical low-frequency earthquake activity models

Slow earthquakes are slow fault slip events. Quantifying and monitoring slow earthquake activity characteristics are important, because they may change before large earthquakes occur. Statistical seismicity models are useful for quantifying seismicity characteristics. However, no standard statistical model exists for slow earthquake activity. This study used a high-quality catalog of low-frequency earthquakes (LFEs), a type of slow earthquake, in the Nankai subduction zone from April 2004 to August 2015 and conducted the first comparison of existing statistical LFE activity models to determine which model better describes LFE activity. Based on this comparison, this study proposes a new hybrid model that incorporates existing model features. The new model considers the LFE activity history in a manner similar to the epidemic-type aftershock sequence (ETAS) model and represents the LFE aftershock rate (subsequent LFE occurrence rate) with a small number of model parameters, as in the Omori–Utsu aftershock law for regular earthquakes. The results show that the proposed model outperforms other existing models. However, the new model cannot reproduce a feature of LFE activity: the sudden cessation of intense LFE bursts. This is because the new model superimposes multiple aftershock activities and predicts extremely high seismicity rates during and after the LFE bursts. I suggest that reproducing and successfully predicting the sudden cessation of intense LFE bursts is critical for the further improvement of statistical LFE activity models. In addition, the empirical equations formulated in this study for the LFE aftershock rates may be useful for future statistical and physical modeling of LFE activity.Graphical

The effect of temperature on injection-induced shear slip of laboratory faults in sandstone

Fluid injection into subsurface reservoirs may cause existing faults/fractures to slip seismically. To study the effect of temperature on injection-induced fault slip, at a constant confining pressure of 10 MPa, we performed a series of injection-induced shear slip experiments on critically stressed sandstone samples containing saw-cut fractures (laboratory-simulated faults) under varying fluid pressurization rates (0.1 and 0.5 MPa/min, respectively) and temperatures (25, 80, and 140 °C, respectively). At 25 °C, slow fault slip events with a peak slip velocity of about 0.13 μm/s were observed on a tested sample in response to a low fluid pressurization rate of 0.1 MPa/min. In contrast, fluid injection with a high pressurization rate of 0.5 MPa/min caused fault slip events with a peak slip rate up to about 0.38 μm/s. In response to a given fluid pressurization rate, several episodes of slip events with a higher slip velocity were induced at an elevated temperature of 140 °C, indicating an appreciable weakening effect at elevated temperatures. We also experimentally constrained the rate-and-state frictional (RSF) parameters at varying effective normal stresses and temperatures by performing velocity-stepping tests. The obtained RSF parameters demonstrate that for a relatively high normal stress, increasing temperature tends to destabilize fault slip. Post-mortem microstructural observations reveal that elevated temperatures promote the generation of abundant fine-grained gouge particles associated with injection-induced shear slip. Our experiments highlight that injection-induced fault slip is affected by temperature-related wear production over the fault surface.

Imaging the Spatiotemporal Evolution of Plate Coupling With Interferometric Radar (InSAR) in the Hikurangi Subduction Zone

AbstractThe coupling at the interface between tectonic plates is a key geophysical parameter to capture the frictional locking across plate boundaries and provides a means to estimate where tectonic strain is accumulating through time. Here, we use both interferometric radar (InSAR) and Global Navigation Satellite System (GNSS) data to investigate the plate coupling of the Hikurangi subduction zone beneath the North Island of New Zealand, where multiple slow slip cycles are superimposed on the long‐term loading. We estimate the plate coupling across the subduction zone over three multi‐year observational periods targeting different stages of the slow slip cycle. Our results highlight the importance of the observational time period when interpreting coupling maps, emphasizing the temporal variability of plate coupling. Leveraging multiple geodetic data sets, we demonstrate how InSAR provides powerful constraints on the spatial resolution of both plate coupling and slow fault slip, even in a region where a dense GNSS network exists.

Subdaily Slow Fault Slip Dynamics Captured by Low‐Frequency Earthquakes

AbstractGeodetic positioning is the geophysical record of reference for slow slip events, but typical daily solutions limit studies of the evolution of slow slip to its long‐term dynamics. Accompanying seismic low‐frequency earthquakes located precisely in time and space provide an opportunity to image slow slip dynamics at subdaily time scales. Here we show that a high‐resolution time history of low‐frequency earthquake fault slip alone can reproduce the geodetic record of slow slip that we observe to be dominated by subdaily fault slip dynamics. However, a simple linear model cannot accommodate the complex dynamics present throughout the slow slip cycle, and an analysis of different phases of the slow slip cycle shows that the ratio of geodetic to seismic fault slip varies as a function of time. This suggests that the low‐frequency earthquake source region saturates as slow slip grows in moment and area. We propose that rheological heterogeneities at the plate boundary associated with low‐frequency earthquakes do not play a significant role in the slow slip rupture process, thus implying that their activity is incidental to the driving aseismic slip.

Measurements of Wave‐Induced Attenuation in Saturated Metapelite and the Band‐Limitation of Low‐Frequency Earthquakes

AbstractThe most common explanation for the depletion of high frequency waves that defines low‐frequency earthquakes (LFEs) and very low‐frequency earthquakes (VLFEs) is that fault rupture and slip are slower than typical earthquakes. However, it is difficult to rule out the possibility that the high frequency waves are produced during slip, but attenuated near the LFE source. One reason this hypothesis has been poorly tested is that there are no measurements of attenuation on the relevant rocks. We present the results of forced oscillation experiments that measure the frequency‐dependent attenuation of a chlorite‐rich metapelitic schist, a lithology found along the subduction plate boundary where LFEs and VLFEs have been documented. Experiments were run on dry and water‐saturated schist at effective pressures of 2–10 MPa and at frequencies of 2 × 10−5–30 Hz. We find that pore fluids and low effective pressure result in the attenuation of high frequencies. The frequency‐dependent attenuation is consistent with the concomitant operation of two wave‐induced fluid flow mechanisms, squirt flow, and patchy saturation. When the effects of these mechanisms are extrapolated to geologic conditions using rock physics models, our results predict that attenuation is capable of completely diminishing the frequencies depleted in LFEs and VLFEs. Therefore, LFEs and VLFEs may not necessarily record slow fault slip, but possibly the presence of high fluid pressure.

Implications of slow fault slip for hydraulic-fracturing-induced seismicity

Since the discovery of slow-slip phenomena, scientific understanding of the behavior of active fault systems has been transformed significantly. It is now recognized that tectonic fault systems are characterized by a spectrum of slip behavior, from “regular” (stick-slip) earthquakes that radiate elastic wave energy and occur on a timescale of seconds, to slow-slip events with durations ranging from minutes to years. More recently, slow-slip phenomena have been observed and modeled in association with injection-induced seismicity. This includes evidence for predominantly slow fault slip during injection that triggered dynamic rupture elsewhere on a fault. In the case of hydraulic fracturing, slow-slip behavior is consistent with the frictional characteristics of faults in clay-rich rocks. A change in pore pressure or slip rate can cause a fault to transition from slow to unstable slip. Through real-time monitoring of slip-slip processes and, potentially, the development of operational adjustments to reduce the hazard of damaging ground motions, a better understanding of slow-slip processes could contribute to improved risk mitigation for induced earthquakes. Effective tools for direct observation of slow-slip processes include tiltmeters, strainmeters, global navigation satellite systems, interferometric synthetic aperture radar, and distributed fiber-optic sensing.

Identification of Low‐Frequency Earthquakes on the San Andreas Fault With Deep Learning

AbstractLow‐frequency earthquakes are a seismic manifestation of slow fault slip. Their emergent onsets, low amplitudes, and unique frequency characteristics make these events difficult to detect in continuous seismic data. Here, we train a convolutional neural network to detect low‐frequency earthquakes near Parkfield, CA using the catalog of Shelly (2017), https://doi.org/10.1002/2017jb014047 as training data. We explore how varying model size and targets influence the performance of the resulting network. Our preferred network has a peak accuracy of 85% and can reliably pick low‐frequency earthquake (LFE) S‐wave arrival times on single station records. We demonstrate the abilities of the network using data from permanent and temporary stations near Parkfield, and show that it detects new LFEs that are not part of the Shelly (2017), https://doi.org/10.1002/2017jb014047 catalog. Overall, machine‐learning approaches show great promise for identifying additional low‐frequency earthquake sources. The technique is fast, generalizable, and does not require sources to repeat.

Stress Orientations in the Nankai Trough Constrained Using Seismic and Aseismic Slip

AbstractThe strength of subduction thrust faults is key to understanding seismogenesis at the provenance of Earth's largest earthquakes. Earthquake focal mechanisms are routinely inverted to constrain the stress state at seismogenic depths. However, on some megathrusts, deformation is accommodated by both earthquakes and types of slow fault slip. We employ focal mechanisms of short‐term slow slip events (SSEs), a type of slow fault slip, and earthquakes in a regional stress inversion to investigate the stress state of the Nankai Trough megathrust and interpret the results in the context of regional tectonics Previous studies using earthquake‐only stress inversions found principal stress orientations in this region that are incompatible with thrust faulting on the megathrust. When both SSEs and earthquakes are considered, the stress state of the central and eastern Nankai Trough megathrust is well oriented for thrust faulting. Our results suggest that slow fault slip source regions may appear to have misoriented stress fields if slow fault slip constitutes a substantial proportion of fault slip and the stress field is not well constrained by earthquakes. In the SSE region, we find that faults are well oriented for failure, suggesting they have strengths similar to their surroundings. Combined with low Vp/Vs ratios and sensitivity to small stress changes, our results imply that the megathrust and surroundings operate at low deviatoric stresses in the SSE source region. Further, we show that the coefficient of friction for areas hosting SSEs is low (μ=0.19–0.50), implying frictionally weak materials in the SSE source region.

Slow Slip Events in New Zealand

Continuously operating global positioning system sites in the North Island of New Zealand have revealed a diverse range of slow motion earthquakes on the Hikurangi subduction zone. These slow slip events (SSEs) exhibit diverse characteristics, from shallow (<15 km), short (<1 month), frequent (every 1–2 years) events in the northern part of the subduction zone to deep (>30 km), long (>1 year), less frequent (approximately every 5 years) SSEs in the southern part of the subduction zone. Hikurangi SSEs show intriguing relationships to interseismic coupling, seismicity, and tectonic tremor, and they exhibit a diversity of interactions with large, regional earthquakes. Due to the marked along-strike variations in Hikurangi SSE characteristics, which coincide with changes in physical characteristics of the subduction margin, the Hikurangi subduction zone presents a globally unique natural laboratory to resolve outstanding questions regarding the origin of episodic, slow fault slip behavior. ▪ New Zealand's Hikurangi subduction zone hosts slow slip events with a diverse range of depth, size, duration, and recurrence characteristics. ▪ Hikurangi slow slip events show intriguing relationships with seismicity ranging from small earthquakes and tremor to larger earthquakes. ▪ Slow slip events play a major role in the accommodation of plate motion at the Hikurangi subduction zone. ▪ Many aspects of the Hikurangi subduction zone make it an ideal natural laboratory to resolve the physical processes controlling slow slip.

Implications of basement rock alteration in the Nankai Trough, Japan for subduction megathrust slip behavior

The Nankai Trough is an exceptionally well-studied convergent margin known to host damaging megathrust earthquakes as well as various forms of slow fault slip. Outcrop studies of exhumed analogues for the modern subduction thrust, as well as seismic reflection images of the active subduction zone, suggest that the basaltic basement participates in plate-boundary shearing at depths of a few km and greater. We obtained altered basaltic upper basement from the modern Nankai Trough seaward of the trench, recovered by drilling on IODP Expedition 333. We performed laboratory friction experiments on bare surfaces and gouge powders of this material, in addition to and in combination with other materials for comparison. The altered basalt exhibits predominantly velocity-strengthening frictional behavior, indicating a tendency for stable slip. For bare surface experiments, few instances of velocity-weakening friction occur and are restricted to slip rates of <10−6 m/s. Thin sections and XRD analyses of the starting material indicate that the velocity-strengthening behavior is likely associated with the presence of clay minerals, mostly Mg-smectite, which coat the rims of larger, stronger grains. Based our friction data, we suggest that in addition to creep the altered Nankai basalts may also allow the possibility of slow slip events for two reasons: (1) velocity-weakening friction may be expected at low slip rates, but velocity-strengthening at higher rates will damp any potential slip instabilities; and (2) based on a critical stiffness criterion for accelerating (stable) slip, SSEs are capable of nucleating in the altered Nankai basalt despite velocity-strengthening friction. In light of these observations, and given the depth at which altered basement appears to be incorporated along the plate interface, we raise the possibility that some shallow slow slip events in the Nankai accretionary prism could originate from shearing in altered basalt.

Fluid overpressure from chemical reactions in serpentinite within the source region of deep episodic tremor

Slow fault slip includes a range of transient phenomena that occur over timescales longer than those of standard earthquakes. Slow slip events are often closely associated with swarms of tectonic tremor. Deep episodic tremor and slip close to the slab–mantle interface in subduction zones has been linked to high fluid pressures produced by dehydration of the subducting slab at greater depths. The slab–mantle interface is a fundamental chemical boundary, where mantle rocks are sheared and mixed with oceanic slab lithologies in a highly reactive environment to form serpentinite. Here we present field and microstructural observations from the plate boundary-scale crustal Livingstone Fault in New Zealand that suggest chemical reactions involving serpentinite can promote rock hardening and generate in situ fluid overpressures. We infer that these processes collectively can result in hydrofracturing and a transition from distributed creep to localized brittle failure and faulting. Serpentinite-related reactions occur over a wide range of pressure and temperature conditions that overlap with those in many forearc mantle wedges. We conclude that the release of fluids derived from such reactions may be an additional and widespread mechanism to generate high fluid pressure patches and brittle failure in the source region of deep tremor along the slab–mantle interface. Chemical reactions between slab and mantle rocks may lead to brittle failure where deep episodic tremor occurs in subduction zones, according to field and microstructural observations of a shear zone in New Zealand.

Transient compression fault slip detected within andesitic rocks of the Casma Group, Lima, Peru

A fault slip within the Nana tunnel Lima, Peru has been monitored since 2012. The data are recorded using an optical-mechanical 3D extensometer, capable of providing very precise long-term three-dimensional measurements of relative displacement across discontinuities. The Nana tunnel has an extremely stable environment and cannot possibly be affected by gravitationally-induced mass movements. The host rock of the tunnel is an aphanitic hornblende-bearing basaltic trachyandesite. Several fault and fracture zones detected in the tunnel represent the major discontinuities of the broader surrounding of the monitored site. The recorded fault slip on the NNW-SSE and E-W striking fracture and fault, with inclinations of 76° to WSW and 78° to N respectively, appoint to compressional event (discontinuity contraction) with a maximum horizontal compression stress axis oriented approximately WSW to ENE corresponding to the direction of the Nazca and South America plates convergence. This event resulted in an aseismic slow fault slip between July 2012 and May 2013. The anticipated compression orientation matches the previously published in-situ stress measurements and fault plane solutions, as well as GPS measurements of the movements of the corresponding part of the Peruvian coast. Nevertheless, the presented monitoring results reflect only short-term fault slip dynamics and need to be considered with caution, even though they correspond to the overall tectonic activity driven by continental subduction.

The Role of Near‐Fault Relief Elements in Creating and Maintaining a Strike‐Slip Landscape

AbstractStrike‐slip landscapes are often associated with a suite of characteristic geomorphic features that provide primary evidence for interpreting fault slip histories. Here we explore the role of shutter ridges, areas of relief advected laterally along faults, in generating two classic strike‐slip processes: progressive lateral offset of channels and stream capture. Landscape models and comparative analysis of the Marlborough Fault System, NZ, show that the length of channel offsets observable in a landscape is primarily controlled by the length of shutter ridges. In our simple landscape model, this scale is controlled by the drainage spacing, and therefore by the geometry of the mountain range. In a more complex landscape, this scale may be controlled by lithologic or structural contrasts. We also find that shutter ridge relief inhibits stream capture, especially at slow fault slip rates relative to hillslope erosion rates. In this case, lateral drainage advection enables streams to “outrun” capture.

The geophysics, geology and mechanics of slow fault slip

Modern geodetic and seismologic observations describe the behavior of fault slip over a vast range of spatial and temporal scales. Slip at sub-seismogenic speeds is evident from top to bottom of lithospheric faults and plays an important role throughout the earthquake cycle. Where earthquakes and tremor accompany slow slip, they help illuminate the spatiotemporal evolution of fault slip. Geophysical subsurface imaging and geologic field studies provide information about suitable environments of slow slip. In particular, exhumed fault and shear zones from various depths reveal the importance of multiple deformation processes and fault-zone structures. Most geologic examples feature frictionally weak and velocity-strengthening materials, well-developed mineral fabrics, and abundant veining indicative of near-lithostatic fluid pressure. To produce transient slow slip events and tremor, in addition to the presence of high-pressure fluids a heterogeneous fault-zone structure, composition, and/or metamorphic assemblage may be needed. Laboratory and computational models suggest that velocity-weakening slip patches smaller than a critical dimension needed for earthquake nucleation will also fail in slow slip events. Changes in fluid pressure or slip rate can cause a fault to transition between stable and unstable fault slip behavior. Future interdisciplinary investigations of slow fault slip, directly integrating geophysical, geological and modeling investigations, will further improve our understanding of the dynamics of slow slip and aid in providing more accurate earthquake hazard characterizations.

Kinematics of the 2015 San Ramon, California earthquake swarm: Implications for fault zone structure and driving mechanisms

Earthquake swarms represent a sudden increase in seismicity that may indicate a heterogeneous fault-zone, the involvement of crustal fluids and/or slow fault slip. Swarms sometimes precede major earthquake ruptures. An earthquake swarm occurred in October 2015 near San Ramon, California in an extensional right step-over region between the northern Calaveras Fault and the Concord–Mt. Diablo fault zone, which has hosted ten major swarms since 1970. The 2015 San Ramon swarm is examined here from 11 October through 18 November using template matching analysis. The relocated seismicity catalog contains ∼4000 events with magnitudes between −0.2<Md<3.6. The swarm illuminated three sub-parallel, southwest striking and northwest dipping fault segments of km-scale dimension and thickness of up to 200 m. The segments contain coexisting populations of different focal-mechanisms, suggesting a complex fault zone structure with several sets of en échelon fault orientations. The migration of events along the three planar structures indicates a complex fluid and faulting interaction processes. We searched for correlations between seismic activity and tidal stresses and found some suggestive features, but nothing that we can be confident is statistically significant.

Large-scale dynamic triggering of shallow slow slip enhanced by overlying sedimentary wedge

Seismic waves can trigger further fault slip. Analysis of seismic and geodetic data shows that seismic waves from the 2016 Kaikōura, New Zealand earthquake were amplified by subduction zone sediments, triggering slow fault slip up to 600 km away.

Coseismic slip propagation on the Tohoku plate boundary fault facilitated by slip‐dependent weakening during slow fault slip

AbstractSlow slip phenomena, including earthquake afterslip and discrete slow slip events (SSEs), may have a major effect on megathrust earthquakes in subduction zones; however, the nature of this effect is still incompletely understood. Here we report that in friction experiments using samples from the plate boundary fault zone in the 2011 Tohoku‐Oki earthquake, an increase in sliding velocity can induce a change from steady state frictional strength or slip‐strengthening friction to slip‐weakening frictional behavior. Specifically, significant slip weakening is caused by velocity steps with initial slip rates consistent with afterslip observed after the largest Tohoku earthquake foreshock on 9 March 2011. This suggests that the foreshock afterslip, which is slightly faster than discrete SSEs, is favorable for fault weakening that may have contributed to the large coseismic slip during the main shock on the shallow plate boundary fault during the Tohoku‐Oki earthquake.

Slow fault propagation in serpentinite under conditions of high pore fluid pressure

The rupture, localization, and slip of faults in serpentinite were studied under varying pore fluid pressure conditions to understand deformation mechanisms potentially responsible for slow slip in fault zones. Experiments were conducted at a constant effective confining pressure of 10 MPa and under pore fluid pressures from 0 to 120 MPa and at temperatures from 23 to 110 °C. With no fluid pressure, faulting occurs rapidly and audibly, and the duration of failure increases monotonically with increasing fluid pressure and temperature. Although non-dilatant during initial strain hardening, the serpentinite dilates during strain weakening concomitant with fault rupture and slip. Non-dilatant strain hardening occurs by microcracking along serpentine basal planes and grain boundaries and rarely in mode I orientations, consistent with previous studies. Dilatant fault rupture produces a network of transgranular shear fractures in conjugate orientations, generally with one dominant fracture. Structural observations show that as fluid pressure increases, the number of transgranular fractures increases. We propose that when faulting occurs over a distributed zone rather than a pre-existing principal slip surface, dilatant hardening causes deformation to migrate. This process causes an increase in slip weakening distance and fracture energy at elevated fluid pressures that can lead to more stable failure. Further, thermally-activated processes caused deformation at propagating crack tips, which also increases the slip weakening distance and the effective fracture energy with increasing temperature. Given the geologic settings for slow slip, our results indicate that high fluid pressure, distributed deformation, and thermally-activated processes may all contribute to slow fault rupture and slip.