Earthquake Source Properties from Pseudotachylite

Earthquake‐radiated motions contain information that can be interpreted as source displacement and therefore related to stress drop. Except in a few notable cases, these displacements cannot be easily related to the absolute stress level or the fault strength, or attributed to a particular physical mechanism. In contrast, paleoearthquakes recorded by exhumed pseudotachylite have a known dynamic mechanism whose properties constrain the coseismic fault strength. Pseudotachylite can be used to directly address a discrepancy between seismologically measured stress drops, which are typically a few MPa, and much larger dynamic stress drops expected from thermal weakening during slip at seismic speeds in crystalline rock (Mckenzie and Brune, 1972; Sibson, 1973; Lachenbruch, 1980; Mase and Smith, 1987; Rice, 2006), and as have been observed in laboratory experiments at high slip rates (Di Toro, Hirose, Nielsen, Pennacchioni, et al., 2006). This places pseudotachylite‐derived estimates of fault strength and inferred crustal stress within the context and bounds of naturally observed earthquake source parameters: apparent stress, stress drop, and overshoot, including consideration of fault‐surface roughness, off‐fault damage, fracture energy, and the strength excess. The analysis, which assumes stress drop is related to corner frequency as in the Madariaga (1976) source model, is restricted to earthquakes of the Gole Larghe fault zone in the Italian Alps, where the dynamic shear strength is well constrained by field and laboratory measurements. We find that radiated energy is similar to or exceeds the shear‐generated heat and that the maximum strength excess is ∼16  MPa. These events have inferred earthquake source parameters that are rare, for instance, a low percentage of the global earthquake population has stress drops as large, unless fracture energy is routinely greater than in existing models, pseudotachylite is not representative of the shear strength during the earthquake that generated it, or the strength excess is larger than we have allowed

Quantifying damage, saturation and anisotropy in cracked rocks by inverting elastic wave velocities

Crack damage results in a decrease of elastic wave velocities and in the development of anisotropy. Using non-interactive crack effective medium theory as a fundamental tool, we calculate dry and wet elastic properties of cracked rocks in terms of a crack density tensor, average crack aspect ratio and mean crack fabric orientation from the solid grains and fluid elastic properties. Using this same tool, we show that both the anisotropy and shear-wave splitting of elastic waves can be derived. Two simple crack distributions are considered for which the predicted anisotropy depends strongly on the saturation, reaching up to 60% in the dry case. Comparison with experimental data on two granites, a basalt and a marble, shows that the range of validity of the non-interactive effective medium theory model extends to a total crack density of approximately 0.5, considering symmetries up to orthorhombic. In the isotropic case, Kachanov's (1994) non-interactive effective medium model was used in order to invert elastic wave velocities and infer both crack density and aspect ratio evolutions. Inversions are stable and give coherent results in terms of crack density and aperture evolution. Crack density variations can be interpreted in terms of crack growth and/or changes of the crack surface contact areas as cracks are being closed or opened respectively. More importantly, the recovered evolution of aspect ratio shows an exponentially decreasing aspect ratio (and therefore aperture) with pressure, which has broader geophysical implications, in particular on fluid flow. The recovered evolution of aspect ratio is also consistent with current mechanical theories of crack closure. In the anisotropic cases—both transverse isotropic and orthorhombic symmetries were considered—anisotropy and saturation patterns were well reproduced by the modelling, and mean crack fabric orientations we recovered are consistent with in situ geophysical imaging. Our results point out that: (1) It is possible to predict damage, anisotropy and saturation in terms of a crack density tensor and mean crack aspect ratio and orientation; (2) using well constrained wave velocity data, it is possible to extrapolate the contemporaneous evolution of crack density, anisotropy and saturation using wave velocity inversion as a tool; 3) using such an inversion tool opens the door in linking elastic properties, variations to permeability