Is complex fault zone behaviour a reflection of rheological heterogeneity?

Fault slip speeds range from steady plate boundary creep through to earthquake slip. Geological descriptions of faults range from localized displacement on one or more discrete planes, through to distributed shearing flow in tabular zones of finite thickness, indicating a large range of possible strain rates in natural faults. We review geological observations and analyse numerical models of two-phase shear zones to discuss the degree and distribution of fault zone heterogeneity and effects on active fault slip style. There must be certain conditions that produce earthquakes, creep and slip at intermediate velocities. Because intermediate slip styles occur over large ranges in temperature, the controlling conditions must be effects of fault properties and/or other dynamic variables. We suggest that the ratio of bulk driving stress to frictional yield strength, and viscosity contrasts within the fault zone, are critical factors. While earthquake nucleation requires the frictional yield to be reached, steady viscous flow requires conditions far from the frictional yield. Intermediate slip speeds may arise when driving stress is sufficient to nucleate local frictional failure by stress amplification, or local frictional yield is lowered by fluid pressure, but such failure is spatially limited by surrounding shear zone stress heterogeneity. This article is part of a discussion meeting issue 'Understanding earthquakes using the geological record'.

No Significant Effect of Coulomb Stress on the Gutenberg-Richter Law after the Landers Earthquake

Coulomb-stress theory has been used for years in seismology to understand how earthquakes trigger each other. Whenever an earthquake occurs, the stress field changes, and places with positive increases are brought closer to failure. Earthquake models that relate earthquake rates and Coulomb stress after a main event, such as the rate-and-state model, assume that the magnitude distribution of earthquakes is not affected by the change in the Coulomb stress. By using different slip models, we calculate the change in Coulomb stress in the fault plane for every aftershock after the Landers event (California, USA, 1992, moment magnitude 7.3). Applying several statistical analyses to test whether the distribution of magnitudes is sensitive to the sign of the Coulomb-stress increase, we are not able to find any significant effect. Further, whereas the events with a positive increase of the stress are characterized by a much larger proportion of strike-slip events in comparison with the seismicity previous to the mainshock, the events happening despite a decrease in Coulomb stress show no relevant differences in focal-mechanism distribution with respect to previous seismicity.

Seismological asperities from the point of view of dynamic rupture modeling: the 2007 Mw6.6 Chuetsu-Oki, Japan, earthquake

We study the ground motion simulations based on three finite-source models for the 2007 Mw6.6 Niigata Chuetsu-oki, Japan, earthquake in order to discuss the performance of the input ground motion estimations for the near-field seismic hazard analysis. The three models include a kinematic source inverted from the regional accelerations, a dynamic source on a planar fault with three asperities inferred from the very-near-field ground motion particle motions, and another dynamic source model with conjugate fault segments. The ground motions are calculated for an available 3D geological model using a finite-difference method. For the comparison, we apply a goodness-of-fit score to the ground motion parameters at different stations, including the nearest one that is almost directly above the ruptured fault segments. The dynamic rupture models show good performance. We find that seismologically inferred earthquake asperities on a single fault plane can be expressed with two conjugate segments. The rupture transfer from one segment to another can generate a significant radiation; this could be interpreted as an asperity projected onto a single fault plane. This example illustrates the importance of the fault geometry that has to be taken into account when estimating the very-near-field ground motion.

Earthquake rupture dynamics frozen in exhumed ancient faults

Most of our knowledge about co-seismic rupture propagation is derived from inversion and interpretation of strong-ground-motion seismograms, laboratory experiments on rock and rock-analogue material, or inferred from theoretical and numerical elastodynamic models. However, additional information on dynamic rupture processes can be provided by direct observation of faults exhumed at the Earth's surface. Pseudotachylytes (solidified friction-induced melts) are the most certain fault-rock indicator of seismicity on ancient faults. Here we show how the asymmetry in distribution and the orientation of pseudotachylyte-filled secondary fractures around an exhumed fault can be used to reconstruct the earthquake rupture directivity, rupture velocity and fracture energy, by comparison with the theoretical dynamic stress field computed around propagating fractures. In particular, the studied natural network of pseudotachylytes is consistent with a dominant propagation direction during repeated seismic events and subsonic rupture propagation close to the Rayleigh wave velocity.

A supershear transition mechanism for cracks

Seismic data indicate that fault ruptures follow complicated paths with variable velocity because of inhomogeneities in initial stress or fracture energy. We report a phenomenon unique to three-dimensional cracks: Locally stronger fault sections, rather than slowing ruptures, drive them forward at velocities exceeding the shear wave speed. This supershear mechanism differentiates barrier and asperity models of fault heterogeneity, which previously have been regarded as indistinguishable. High strength barriers concentrate energy, producing potentially destructive pulses of strong ground motion.