David R. Shelly

Low-frequency earthquakes in Shikoku, Japan, and their relationship to episodic tremor and slip

Non-volcanic seismic tremor was discovered in the Nankai trough subduction zone in southwest Japan and subsequently identified in the Cascadia subduction zone. In both locations, tremor is observed to coincide temporally with large, slow slip events on the plate interface downdip of the seismogenic zone. The relationship between tremor and aseismic slip remains uncertain, however, largely owing to difficulty in constraining the source depth of tremor. In southwest Japan, a high quality borehole seismic network allows identification of coherent S-wave (and sometimes P-wave) arrivals within the tremor, whose sources are classified as low-frequency earthquakes. As low-frequency earthquakes comprise at least a portion of tremor, understanding their mechanism is critical to understanding tremor as a whole. Here, we provide strong evidence that these earthquakes occur on the plate interface, coincident with the inferred zone of slow slip. The locations and characteristics of these events suggest that they are generated by shear slip during otherwise aseismic transients, rather than by fluid flow. High pore-fluid pressure in the immediate vicinity, as implied by our estimates of seismic P- and S-wave speeds, may act to promote this transient mode of failure. Low-frequency earthquakes could potentially contribute to seismic hazard forecasting by providing a new means to monitor slow slip at depth.

Non-volcanic tremor and low-frequency earthquake swarms

Non-volcanic tremor is a weak, extended duration seismic signal observed episodically on some major faults, often in conjunction with slow slip events. Such tremor may hold the key to understanding fundamental processes at the deep roots of faults, and could signal times of accelerated slip and hence increased seismic hazard. The mechanism underlying the generation of tremor and its relationship to aseismic slip are, however, as yet unresolved. Here we demonstrate that tremor beneath Shikoku, Japan, can be explained as a swarm of small, low-frequency earthquakes, each of which occurs as shear faulting on the subduction-zone plate interface. This suggests that tremor and slow slip are different manifestations of a single process.

A scaling law for slow earthquakes.

Recently, a series of unusual earthquake phenomena have been discovered, including deep episodic tremor, low-frequency earthquakes, very-low-frequency earthquakes, slow slip events and silent earthquakes. Each of these has been demonstrated to arise from shear slip, just as do regular earthquakes, but with longer characteristic durations and radiating much less seismic energy. Here we show that these slow events follow a simple, unified scaling relationship that clearly differentiates their behaviour from that of regular earthquakes. We find that their seismic moment is proportional to the characteristic duration and their moment rate function is constant, with a spectral high-frequency decay of f-1. This scaling and spectral behaviour demonstrates that they can be thought of as different manifestations of the same phenomena and that they comprise a new earthquake category. The observed scale dependence of rupture velocity for these events can be explained by either a constant low-stress drop model or a diffusional constant-slip model. This new scaling law unifies a diverse class of slow seismic events and may lead to a better understanding of the plate subduction process and large earthquake generation.

Deep low-frequency earthquakes in tremor localize to the plate interface in multiple subduction zones

Deep tremor under Shikoku, Japan, consists primarily, and perhaps entirely, of swarms of low-frequency earthquakes (LFEs) that occur as shear slip on the plate interface. Although tremor is observed at other plate boundaries, the lack of cataloged low-frequency earthquakes has precluded a similar conclusion about tremor in those locales. We use a network autocorrelation approach to detect and locate LFEs within tremor recorded at three subduction zones characterized by different thermal structures and levels of interplate seismicity: southwest Japan, northern Cascadia, and Costa Rica. In each case we find that LFEs are the primary constituent of tremor and that they locate on the deep continuation of the plate boundary. This suggests that tremor in these regions shares a common mechanism and that temperature is not the primary control on such activity. Citation: Brown, J.R., G. C. Beroza, S. Ide, K. Ohta, D. R. Shelly, S. Y. Schwartz, W. Rabbel, M. Thorwart, and H. Kao (2009), Deep low-frequency earthquakes in tremor localize to the plate interface in multiple subduction zones, Geophys. Res. Lett., 36, L19306, doi:10.1029/2009GL040027.

Earthquakes triggered by silent slip events on Kīlauea volcano, Hawaii.

Slow-slip events, or ‘silent earthquakes’, have recently been discovered in a number of subduction zones including the Nankai trough in Japan, Cascadia, and Guerrero in Mexico, but the depths of these events have been difficult to determine from surface deformation measurements. Although it is assumed that these silent earthquakes are located along the plate megathrust, this has not been proved. Slow slip in some subduction zones is associated with non-volcanic tremor, but tremor is difficult to locate and may be distributed over a broad depth range. Except for some events on the San Andreas fault, slow-slip events have not yet been associated with high-frequency earthquakes, which are easily located. Here we report on swarms of high-frequency earthquakes that accompany otherwise silent slips on Kīlauea volcano, Hawaii. For the most energetic event, in January 2005, the slow slip began before the increase in seismicity. The temporal evolution of earthquakes is well explained by increased stressing caused by slow slip, implying that the earthquakes are triggered. The earthquakes, located at depths of 7–8 km, constrain the slow slip to be at comparable depths, because they must fall in zones of positive Coulomb stress change. Triggered earthquakes accompanying slow-slip events elsewhere might go undetected if background seismicity rates are low. Detection of such events would help constrain the depth of slow slip, and could lead to a method for quantifying the increased hazard during slow-slip events, because triggered events have the potential to grow into destructive earthquakes.

Migrating tremors illuminate complex deformation beneath the seismogenic San Andreas fault

The San Andreas fault is one of the most extensively studied faults in the world, yet its physical character and deformation mode beneath the relatively shallow earthquake-generating portion remain largely unconstrained. Tectonic ‘non-volcanic’ tremor, a recently discovered seismic signal probably generated by shear slip on the deep extension of some major faults, can provide new insight into the deep fate of such faults, including that of the San Andreas fault near Parkfield, California. Here I examine continuous seismic data from mid-2001 to 2008, identifying tremor and decomposing the signal into different families of activity based on the shape and timing of the waveforms at multiple stations. This approach allows differentiation between activities from nearby patches of the deep fault and begins to unveil rich and complex patterns of tremor occurrence. I find that tremor exhibits nearly continuous migration, with the most extensive episodes propagating more than 20 kilometres along fault strike at rates of 15–80 kilometres per hour. This suggests that the San Andreas fault remains a localized through-going structure, at least to the base of the crust, in this area. Tremor rates and recurrence behaviour changed markedly in the wake of the 2004 magnitude-6.0 Parkfield earthquake, but these changes were far from uniform within the tremor zone, probably reflecting heterogeneous fault properties and static and dynamic stresses decaying away from the rupture. The systematic recurrence of tremor demonstrated here suggests the potential to monitor detailed time-varying deformation on this portion of the deep San Andreas fault, deformation which unsteadily loads the shallower zone that last ruptured in the 1857 magnitude-7.9 Fort Tejon earthquake.

High-resolution subducting-slab structure beneath northern Honshu, Japan, revealed by double-difference tomography

The high-resolution seismic velocity structure of the subducting slab beneath northern Honshu, Japan, has been obtained by double-difference tomography, capitalizing on the existence of two planes of seismicity. The upper plane lies in the region with average to high Vp/Vs ratios (1.72–1.85), which may correspond to the transformations of metabasalt and metagabbro to blueschist. The lower plane is associated with very low Vp/Vs ratios (1.6–1.7), in sharp contrast with high Vp/Vs ratios (∼1.8–1.85) in the region between the two planes. These features may be explained by forsterite-enstatite-H2O formation from serpentine dehydration in the lower plane and partial hydration of the region between the two planes. Our results provide direct constraints for petrologic modeling of the subducting slab.

Non-volcanic Tremor: A Window into the Roots of Fault Zones

The recent discovery of non-volcanic tremor in Japan and the coincidence of tremor with slow-slip in Cascadia have made earth scientists reevaluate our models for the physical processes in subduction zones and on faults in general. Subduction zones have been studied very closely since the discovery of slow-slip and tremor. This has led to the discovery of a number of related phenomena including low frequency earthquakes and very low frequency earthquakes. All of these events fall into what some have called a new class of events that are governed under a different frictional regime than simple brittle failure. While this model is appealing to many, consensus as to exactly what process generates tremor has yet to be reached. Tremor and related events also provide a window into the deep roots of subduction zones, a poorly understood region that is largely devoid of seismicity. Given that such fundamental questions remain about non-volcanic tremor, slow-slip, and the region in which they occur, we expect that this will be a fruitful field for a long time to come.

Inferring fault rheology from low‐frequency earthquakes on the San Andreas

Families of recurring low-frequency earthquakes (LFEs) within nonvolcanic tremor (NVT) on the San Andreas fault in central California show strong sensitivity to shear stress induced by the daily tidal cycle. LFEs occur at all levels of the tidal shear stress and are in phase with the very small, ~400 Pa, stress amplitude. To quantitatively explain the correlation, we use a model from the existing literature that assumes the LFE sources are small, persistent regions that repeatedly fail during shear of a much larger scale, otherwise aseismically creeping fault zone. The LFE source patches see tectonic loading, creep of the surrounding fault which may be modulated by the tidal stress, and direct tidal loading. If the patches are small relative to the surrounding creeping fault then the stressing is dominated by fault creep, and if patch failure occurs at a threshold stress, then the resulting seismicity rate is proportional to the fault creep rate or fault zone strain rate. Using the seismicity rate as a proxy for strain rate and the tidal shear stress, we fit the data with possible fault rheologies that produce creep in laboratory experiments at temperatures of 400 to 600°C appropriate for the LFE source depth. The rheological properties of rock-forming minerals for dislocation creep and dislocation glide are not consistent with the observed fault creep because strong correlation between small stress perturbations and strain rate requires perturbation on the order of the ambient stress. The observed tidal modulation restricts ambient stress to be at most a few kilopascal, much lower than rock strength. A purely rate dependent friction is consistent with the observations only if the product of the friction rate dependence and effective normal stress is ~ 0.5 kPa. Extrapolating the friction rate strengthening dependence of phyllosilicates (talc) to depth would require the effective normal stress to be ~50 kPa, implying pore pressure is lithostatic. If the LFE source is on the order of tens of meters, as required by the model, rate-weakening friction rate dependence (e.g., olivine) at 400 to 600°C requires that the minimum effective pressure at the LFE source is ~ 2.5 MPa.

Fluid-faulting interactions: Fracture-mesh and fault-valve behavior in the February 2014 Mammoth Mountain, California, earthquake swarm

Faulting and fluid transport in the subsurface are highly coupled processes, which may manifest seismically as earthquake swarms. A swarm in February 2014 beneath densely monitored Mammoth Mountain, California, provides an opportunity to witness these interactions in high resolution. Toward this goal, we employ massive waveform-correlation-based event detection and relative relocation, which quadruples the swarm catalog to more than 6000 earthquakes and produces high-precision locations even for very small events. The swarms main seismic zone forms a distributed fracture mesh, with individual faults activated in short earthquake bursts. The largest event of the sequence, M 3.1, apparently acted as a fault valve and was followed by a distinct wave of earthquakes propagating ~1 km westward from the updip edge of rupture, 1–2 h later. Late in the swarm, multiple small, shallower subsidiary faults activated with pronounced hypocenter migration, suggesting that a broader fluid pressure pulse propagated through the subsurface.