Fault healing promotes high-frequency earthquakes in laboratory experiments and on natural faults.

Faults strengthen or heal with time in stationary contact, and this healing may be an essential ingredient for the generation of earthquakes. In the laboratory, healing is thought to be the result of thermally activated mechanisms that weld together micrometre-sized asperity contacts on the fault surface, but the relationship between laboratory measures of fault healing and the seismically observable properties of earthquakes is at present not well defined. Here we report on laboratory experiments and seismological observations that show how the spectral properties of earthquakes vary as a function of fault healing time. In the laboratory, we find that increased healing causes a disproportionately large amount of high-frequency seismic radiation to be produced during fault rupture. We observe a similar connection between earthquake spectra and recurrence time for repeating earthquake sequences on natural faults. Healing rates depend on pressure, temperature and mineralogy, so the connection between seismicity and healing may help to explain recent observations of large megathrust earthquakes which indicate that energetic, high-frequency seismic radiation originates from locations that are distinct from the geodetically inferred locations of large-amplitude fault slip.

Tremor-tide correlations and near-lithostatic pore pressure on the deep San Andreas fault.

Since its initial discovery nearly a decade ago, non-volcanic tremor has provided information about a region of the Earth that was previously thought incapable of generating seismic radiation. A thorough explanation of the geologic process responsible for tremor generation has, however, yet to be determined. Owing to their location at the plate interface, temporal correlation with geodetically measured slow-slip events and dominant shear wave energy, tremor observations in southwest Japan have been interpreted as a superposition of many low-frequency earthquakes that represent slip on a fault surface. Fluids may also be fundamental to the failure process in subduction zone environments, as teleseismic and tidal modulation of tremor in Cascadia and Japan and high Poisson ratios in both source regions are indicative of pressurized pore fluids. Here we identify a robust correlation between extremely small, tidally induced shear stress parallel to the San Andreas fault and non-volcanic tremor activity near Parkfield, California. We suggest that this tremor represents shear failure on a critically stressed fault in the presence of near-lithostatic pore pressure. There are a number of similarities between tremor in subduction zone environments, such as Cascadia and Japan, and tremor on the deep San Andreas transform, suggesting that the results presented here may also be applicable in other tectonic settings.

Assessing Margin-Wide Rupture Behaviors Along the Cascadia Megathrust With 3-D Dynamic Rupture Simulations.

From California to British Columbia, the Pacific Northwest coast bears an omnipresent earthquake and tsunami hazard from the Cascadia subduction zone. Multiple lines of evidence suggests that magnitude eight and greater megathrust earthquakes have occurred - the most recent being 321 years ago (i.e., 1700 A.D.). Outstanding questions for the next great megathrust event include where it will initiate, what conditions are favorable for rupture to span the convergent margin, and how much slip may be expected. We develop the first 3-D fully dynamic rupture simulations for the Cascadia subduction zone that are driven by fault stress, strength and friction to address these questions. The initial dynamic stress drop distribution in our simulations is constrained by geodetic coupling models, with segment locations taken from geologic analyses. We document the sensitivity of nucleation location and stress drop to the final seismic moment and coseismic subsidence amplitudes. We find that the final earthquake size strongly depends on the amount of slip deficit in the central Cascadia region, which is inferred to be creeping interseismically, for a given initiation location in southern or northern Cascadia. Several simulations are also presented here that can closely approximate recorded coastal subsidence from the 1700 A.D. event without invoking localized high-stress asperities along the down-dip locked region of the megathrust. These results can be used to inform earthquake and tsunami hazards for not only Cascadia, but other subduction zones that have limited seismic observations but a wealth of geodetic inference.

Deep long-period earthquakes generated by second boiling beneath Mauna Kea volcano.

Deep long-period earthquakes (DLPs) are an enigmatic type of volcanic seismicity that sometimes precedes eruptions but mostly occurs at quiescent volcanoes. These earthquakes are depleted in high-frequency content and typically occur near the base of the crust. We observed a near-periodic, long-lived sequence of more than one million DLPs in the past 19 years beneath the dormant postshield Mauna Kea volcano in Hawai'i. We argue that this DLP sequence was caused by repeated pressurization of volatiles exsolved through crystallization of cooling magma stalled beneath the crust. This "second boiling" of magma is a well-known process but has not previously been linked to DLP activity. Our observations suggest that, rather than portending eruptions, global DLP activity may more commonly be indicative of stagnant, cooling magma.

Mega-earthquakes rupture flat megathrusts.

The 2004 Sumatra-Andaman and 2011 Tohoku-Oki earthquakes highlighted gaps in our understanding of mega-earthquake rupture processes and the factors controlling their global distribution: A fast convergence rate and young buoyant lithosphere are not required to produce mega-earthquakes. We calculated the curvature along the major subduction zones of the world, showing that mega-earthquakes preferentially rupture flat (low-curvature) interfaces. A simplified analytic model demonstrates that heterogeneity in shear strength increases with curvature. Shear strength on flat megathrusts is more homogeneous, and hence more likely to be exceeded simultaneously over large areas, than on highly curved faults.