Earthquake rupture dynamics frozen in exhumed ancient faults

Multiscale off-fault brecciation records coseismic energy budget of principal fault zone

Breccia and pulverized rock are typical textures in off-fault damage adjacent to a main seismogenic zone. Previously, by estimating the energy required to advance the rupture in this zone using particle size distribution at sub-millimeter/micrometer scales, we could constrain the energy budget during coseismic events. However, whether microscopic estimation is sufficient to capture surface energy fragmentation during an earthquake and the effect of measurement scale variation on calculation of co-seismic energy partitioning remained unclear. Here, we investigated the mechanism of coseismic off-fault damage based on field and microstructural observations of a well-exposed breccia body in Ichinokawa, Japan. We used in situ clast measurements coupled with thin-section analysis of breccia clasts to estimate the energy budget of the damage zone adjacent to the principal slip zone of the Median Tectonic Line (MTL). The total surface energy density and corresponding surface energy per unit fault for a width of ~ 500 m of the dynamical damage zone were estimated. The moment magnitude estimated based on surface energy was 5.8-8.3 Mw. In Ichinokawa, off-fault fragmentation is initiated by coseismic activity and is followed by fluid activity. Under dynamic fragmentation conditions, the scale is important to calculate the surface energy.

On-fault earthquake energy density partitioning from shocked garnet in an exhumed seismic midcrustal fault

The energy released during an earthquake is mostly dissipated in the fault zone and subordinately as radiated seismic waves. The on-fault energy budget is partitioned into frictional heat, generation of new grain surface by microfracturing, and crystal-lattice distortion associated with dislocation defects. The relative contribution of these components is debated and difficult to assess, but this energy partitioning strongly influences earthquake mechanics. We use high-resolution scanning-electron-microscopy techniques, especially to analyze shocked garnet in a fault wall-rock, to provide the first estimate of all three energy components for a seismic fault patch exhumed from midcrustal conditions. Fault single-jerk seismicity is recorded by the presence of pristine quenched frictional melt. The estimated value of energy per unit fault surface is ~13 megajoules per square meter for heat, which is predominant with respect to both surface energy (up to 0.29 megajoules per square meter) and energy associated with crystal lattice distortion (0.02 megajoules per square meter).

Coseismic fault sealing and fluid pressurization during earthquakes

Earthquakes occur because faults weaken with increasing slip and slip rate. Thermal pressurization (TP) of trapped pore fluids is deemed to be a widespread coseismic fault weakening mechanism. Yet, due to technical challenges, experimental evidence of TP is limited. Here, by exploiting a novel experimental configuration, we simulate seismic slip pulses (slip rate 2.0 m/s) on dolerite-built faults under pore fluid pressures up to 25 MPa. We measure transient sharp weakening, down to almost zero friction and concurrent with a spike in pore fluid pressure, which interrupts the exponential-decay slip weakening. The interpretation of mechanical and microstructural data plus numerical modeling suggests that wear and local melting processes in experimental faults generate ultra-fine materials to seal the pressurized pore water, causing transient TP spikes. Our work suggests that, with wear-induced sealing, TP may also occur in relatively permeable faults and could be quite common in nature.

Dynamic earthquake rupture in the lower crust

Earthquakes in the continental crust commonly occur in the upper 15 to 20 km. Recent studies demonstrate that earthquakes also occur in the lower crust of collision zones and play a key role in metamorphic processes that modify its physical properties. However, details of the failure process and sequence of events that lead to seismic slip in the lower crust remain uncertain. Here, we present observations of a fault zone from the Bergen Arcs, western Norway, which constrain the deformation processes of lower crustal earthquakes. We show that seismic slip and associated melting are preceded by fracturing, asymmetric fragmentation, and comminution of the wall rock caused by a dynamically propagating rupture. The succession of deformation processes reported here emphasize brittle failure mechanisms in a portion of the crust that until recently was assumed to be characterized by ductile deformation.

Heating, weakening and shear localization in earthquake rupture

Field and borehole observations of active earthquake fault zones show that shear is often localized to principal deforming zones of order 0.1-10 mm width. This paper addresses how frictional heating in rapid slip weakens faults dramatically, relative to their static frictional strength, and promotes such intense localization. Pronounced weakening occurs even on dry rock-on-rock surfaces, due to flash heating effects, at slip rates above approximately 0.1 m s^-1 (earthquake slip rates are typically of the order of 1 m s^-1). But weakening in rapid shear is also predicted theoretically in thick fault gouge in the presence of fluids (whether native ground fluids or volatiles such as H₂O or CO₂ released by thermal decomposition reactions), and the predicted localizations are compatible with such narrow shear zones as have been observed. The underlying concepts show how fault zone materials with high static friction coefficients, approximately 0.6-0.8, can undergo strongly localized shear at effective dynamic friction coefficients of the order of 0.1, thus fitting observational constraints, e.g. of earthquakes producing negligible surface heat outflow and, for shallow events, only rarely creating extensive melt. The results to be summarized include those of collaborative research published with Nicolas Brantut (University College London), Eric Dunham (Stanford University), Nadia Lapusta (Caltech), Hiroyuki Noda (JAMSTEC, Japan), John D. Platt (Carnegie Institution for Science, now at *gramLabs), Alan Rempel (Oregon State University) and John W. Rudnicki (Northwestern University).This article is part of the themed issue 'Faulting, friction and weakening: from slow to fast motion'.

State-variable friction for the Burridge-Knopoff model

This work shows the relationship of the state variable rock-friction law proposed by Dieterich to the Carlson and Langer friction law commonly used in the Burridge-Knopoff (BK) model of earthquakes. Further to this, the Dieterich law is modified to allow slip rates of zero magnitude yielding a three parameter friction law that is included in the BK system. Dynamic phases of small scale and large scale events are found with a transition surface in the parameter space. Near this transition surface the event size distribution follows a power law with an exponent that varies as the transition is approached contrasting with the invariant exponent observed using the Carlson and Langer friction. This variability of the power-law exponent is consistent with the range of exponents measured in real earthquake systems and is more selective than the range observed in the Olami-Feder-Christensen model.