Aaron G. Wech

Slip rate and tremor genesis in Cascadia

At many plate boundaries, conditions in the transition zone between seismogenic and stable slip produce slow earthquakes. In the Cascadia subduction zone, these events are consistently observed as slow, aseismic slip on the plate interface accompanied by persistent tectonic tremor. However, not all slow slip at other plate boundaries coincides spatially and temporally with tremor, leaving the physics of tremor genesis poorly understood. Here we analyze seismic, geodetic, and strainmeter data in Cascadia to observe for the first time a large, tremor-generating slow earthquake change from tremor-genic to silent and back again. The tremor falls silent at reduced slip speeds when the migrating slip front pauses as it loads the stronger adjacent fault segment to failure. The finding suggests that rheology and slip-speed-regulated stressing rate control tremor genesis, and the same section of fault can slip both with and without detectable tremor, limiting tremors use as a proxy for slip.

Cascadia subducting plate fluids channelled to fore-arc mantle corner: ETS and silica deposition

In this study we first summarize the constraints that on the Cascadia subduction thrust, there is a 70 km gap downdip between the megathrust seismogenic zone and the Episodic Tremor and Slip (ETS) that lies further landward; there is not a continuous transition from unstable to conditionally stable sliding. Seismic rupture occurs mainly offshore for this hot subduction zone. ETS lies onshore. We then suggest what does control the downdip position of ETS. We conclude that fluids from dehydration of the downgoing plate, focused to rise above the fore-arc mantle corner, are responsible for ETS. There is a remarkable correspondence between the position of ETS and this corner along the whole margin. Hydrated mineral assemblages in the subducting oceanic crust and uppermost mantle are dehydrated with downdip increasing temperature, and seismic tomography data indicate that these fluids have strongly serpentinized the overlying fore-arc mantle. Laboratory data show that such fore-arc mantle serpentinite has low permeability and likely blocks vertical expulsion and restricts flow updip within the underlying permeable oceanic crust and subduction shear zone. At the fore-arc mantle corner these fluids are released upward into the more permeable overlying fore-arc crust. An indication of this fluid flux comes from low Poissons Ratios (and Vp/Vs) found above the corner that may be explained by a concentration of quartz which has exceptionally low Poissons Ratio. The rising fluids should be silica saturated and precipitate quartz with decreasing temperature and pressure as they rise above the corner.

Reconsidering earthquake scaling

The relationship (scaling) between scalar moment, M0, and duration, T, potentially provides key constraints on the physics governing fault slip. The prevailing interpretation of M0-T observations proposes different scaling for fast (earthquakes) and slow (mostly aseismic) slip populations and thus fundamentally different driving mechanisms. We show that a single model of slip events within bounded slip zones may explain nearly all fast and slow slip M0-T observations, and both slip populations have a change in scaling, where the slip area growth changes from 2-D when too small to sense the boundaries to 1-D when large enough to be bounded. We present new fast and slow slip M0-T observations that sample the change in scaling in each population, which are consistent with our interpretation. We suggest that a continuous but bimodal distribution of slip modes exists and M0-T observations alone may not imply a fundamental difference between fast and slow slip.

Ambient tremors in a collisional orogenic belt

Deep-seated tectonic tremors have been regarded as an observation tied to interconnected fluids at depth, which have been well documented in worldwide subduction zones and transform faults but not in a collisional mountain belt. In this study we explore the general features of collisional tremors in Taiwan and discuss the possible generation mechanism. In the 4 year data, we find 231 ambient tremor episodes with durations ranging from 5 to 30 min. In addition to a coseismic slip-induced stress change from nearby major earthquake, increased tremor rate is also highly correlated with the active, normal faulting earthquake swarms at the shallower depth. Both the tremor and earthquake swarm activities are confined in a small, area where the high attenuation, high thermal anomaly, the boundary between high and low resistivity, and localized veins on the surfaces distributed, suggesting the involvement of fluids from metamorphic dehydration within the orogen.

Cascadia subduction tremor muted by crustal faults

Deep, episodic slow slip on the Cascadia subduction megathrust of western North America is accompanied by low-frequency tremor in a zone of high fluid pressure between 30 and 40 km depth. Tremor density (tremor epicenters per square kilometer) varies along strike, and lower tremor density statistically correlates with upper plate faults that accommodate northward motion and rotation of forearc blocks. Upper plate earthquakes occur to 35 km depth beneath the faults. We suggest that the faults extend to the overpressured megathrust, where they provide fracture pathways for fluid escape into the upper plate. This locally reduces megathrust fluid pressure and tremor occurrence beneath the faults. Damping of tremor and related slow slip caused by fluid escape could affect fault properties of the megathrust, possibly influencing the behavior of great earthquakes.

Extending Alaska’s plate boundary: Tectonic tremor generated by Yakutat subduction

The tectonics of the eastern end of the Alaska-Aleutian subduction zone are complicated by the inclusion of the Yakutat microplate, which is colliding into and subducting beneath continental North America at near-Pacific-plate rates. The interaction among these plates at depth is not well understood, and further east, even less is known about the plate boundary or the source of Wrangell volcanism. The drop-off in Wadati-Benioff zone (WBZ) seismicity could signal the end of the plate boundary, the start of aseismic subduction, or a tear in the downgoing plate. Further compounding the issue is the possible presence of the Wrangell slab, which is faintly outlined by an anemic, eastward-dipping WBZ beneath the Wrangell volcanoes. In this study, I performed a search for tectonic tremor to map slow, plate-boundary slip in south-central Alaska. I identified ∼11,000 tremor epicenters, which continue 85 km east of the inferred Pacific plate edge marked by WBZ seismicity. The tremor zone coincides with the edges of the downgoing Yakutat terrane, and tremors transition from periodic to continuous behavior as they near the aseismic Wrangell slab. I interpret tremor to mark slow, semicontinuous slip occurring at the interface between the Yakutat and North America plates. The slow slip region lengthens the megathrust interface beyond the WBZ and may provide evidence for a connection between the Yakutat slab and the aseismic Wrangell slab.

Temporal variation of tectonic tremor activity in southern Taiwan around the 2010 ML6.4 Jiashian earthquake

Deep tectonic tremor, which is extremely sensitive to small stress variations, could be used to monitor fault-zone processes during large earthquake cycles and aseismic processes before large earthquakes. In this study, we develop an algorithm for the automatic detection and location of tectonic tremor beneath the southern Central Range of Taiwan and examine the spatio-temporal relationship between tremor and the 4 March 2010 ML6.4 Jiashian earthquake, located about 20 km from active tremor sources. We find that tremor in this region has a relatively short duration, short recurrence time, and no consistent correlation with surface GPS data. We find a short-term increase in the tremor rate 19 days before the Jiashian mainshock, and around the time when the tremor rate began to rise, one GPS station recorded a flip in its direction of motion. We hypothesize that tremor is driven by a slow-slip event that preceded the occurrence of the shallower Jiashian mainshock, even though the inferred slip is too small to be observed by all GPS stations. Our study shows that tectonic tremor may reflect stress variation during the pre-nucleation process of a nearby earthquake.

Plateau subduction, intraslab seismicity, and the Denali (Alaska) volcanic gap

Tectonic tremors in Alaska (USA) are associated with subduction of the Yakutat plateau, but their origins are unclear due to lack of depth constraints. We have processed tremor recordings to extract low-frequency earthquakes (LFEs), and generated a set of six LFE waveform templates via iterative network matched filtering and stacking. The timing of impulsive P (compressional) wave and S (shear) wave arrivals on template waveforms places LFEs at 40–58 km depth, near the upper envelope of intraslab seismicity and immediately updip of increased levels of intraslab seismicity. S waves at near-epicentral distances display polarities consistent with shear slip on the plate boundary. We compare characteristics of LFEs, seismicity, and tectonic structures in central Alaska with those in warm subduction zones, and propose a new model for the region’s unusual intraslab seismicity and the enigmatic Denali volcanic gap (i.e., an area of no volcanism where expected). We argue that fluids in the Yakutat plate are confined to its upper crust, and that shallow subduction leads to hydromechanical conditions at the slab interface in central Alaska akin to those in warm subduction zones where similar LFEs and tremor occur. These conditions lead to fluid expulsion at shallow depths, explaining strike-parallel alignment of tremor occurrence with the Denali volcanic gap. Moreover, the lack of double seismic zone and restriction of deep intraslab seismicity to a persistent low-velocity zone are simple consequences of anhydrous conditions prevailing in the lower crust and upper mantle of the Yakutat plate.

Linking magma transport structures at Kīlauea volcano

Identifying magma pathways is important for understanding and interpreting volcanic signals. At Kīlauea volcano, seismicity illuminates subsurface plumbing, but the broad spectrum of seismic phenomena hampers event identification. Discrete, long-period (LP) events dominate the shallow (5–10 km) plumbing, and deep (40+ km) tremor has been observed offshore. However, our inability to routinely identify these events limits their utility in tracking ascending magma. Using envelope cross-correlation, we systematically catalog non-earthquake seismicity between 2008 and 2014. We find that the LPs and deep tremor are spatially distinct, separated by the 15–25 km deep, horizontal mantle fault zone (MFZ). Our search corroborates previous observations, but we find broader band (0.5–20 Hz) tremor comprising collocated earthquakes and reinterpret the deep tremor as earthquake swarms in a volume surrounding and responding to magma intruding from the mantle plume beneath the MFZ. We propose that the overlying MFZ promotes lateral magma transport, linking this deep intrusion with Kīlaueas shallow magma plumbing.

Short-Term Forecasting and Detection of Explosions During the 2016–2017 Eruption of Bogoslof Volcano, Alaska

We describe a multidisciplinary approach to forecast, rapidly detect, and characterize explosive events during the 2016–2017 eruption of Bogoslof volcano, a back-arc shallow submarine volcano in Alaska’s Aleutian arc. The eruptive sequence began in December 2016 and included over 60 discrete explosive events. Because the volcano has no local monitoring stations, we used distant stations on the nearest volcanoes, Okmok (54 km) and Makushin (72 km), combined with regional infrasound sensors and lightning detection from the Worldwide Lightning Location Network (WWLLN). Monitoring of activity used a combination of scheduled checks combined with automated alarms. Alarms triggered on real-time data included real-time seismic amplitude measurement (RSAM); infrasound from several arrays, the closest being on Okmok; and lightning strokes detected from WWLLN within a 20-km radius of the volcano. During periods of unrest, a multidisciplinary response team of four people fulfilled specific roles to evaluate geophysical and remote-sensing data, event-specific ash-cloud dispersion modeling, interagency coordination, and development and distribution of formalized warning products. Using this approach, for events that produced ash clouds ≥7.5 km above sea level, AVO called emergency response partners 15 minutes, and issued written notices 30 minutes, after event onset (mean times). Factors that affect timeliness of written warnings include event size and number of data streams available; bigger events and more data both decrease uncertainty and allow for faster warnings. In remote areas where airborne ash is the primary hazard, the approach used at Bogoslof is an effective strategy for hazard mitigation.