Thibault Candela

Mirror-like faults and power dissipation during earthquakes

Earthquakes occur along faults in response to plate tectonic movements, but paradoxically, are not widely recognized in the geological record, severely limiting our knowledge of earthquake physics and hampering accurate assessments of seismic hazard. Light-reflective (so-called mirror like) fault surfaces are widely observed geological features, especially in carbonate-bearing rocks of the shallow crust. Here we report on the occurrence of mirror-like fault surfaces cutting dolostone gouges in the Italian Alps. Using friction experiments, we demonstrate that the mirror-like surfaces develop only at seismic slip rates (∼1 m/s) and for applied normal stresses and sliding displacements consistent with those estimated on the natural faults. Under these experimental conditions, the frictional power density dissipated in the samples is comparable to that estimated for natural earthquakes (1–10 MW/m 2 ). Our results indicate that mirror-like surfaces in dolostone gouges are a signature of seismic faulting, and can be used to estimate power dissipation during ancient earthquake ruptures.

Stress Drop during Earthquakes: Effect of Fault Roughness Scaling

We propose that a controlling parameter of static stress drop during an earthquake is related to the scaling properties of the fault-surface topography. Using high resolution laser distance meters, we have accurately measured the roughness scaling properties of two fault surfaces in different geological settings (the French Alps and Nevada). The data show that fault-surface topography is scale dependent and may be accurately described by a self-affine geometry with a slight anisotropy characterized by two extreme roughness exponents ( H R ), H ||=0.6 in the direction of slip and H ⊥=0.8 perpendicular to slip. Disregarding plastic processes like rock fragmentation and focusing on elastic deformation of the topography, which is the dominant mode at large scales, the stress drop is proportional to the deformation, which is a spatial derivative of the slip. The evolution of stress-drop fluctuations on the fault plane can be derived directly from the self-affine property of the fault surface, with the length scale ( λ ) as std Δσ ( λ )∝ λ H R -1. Assuming no characteristic length scale in fault roughness and a rupture cascade model, we show that as the rupture grows, the average stress drop, and its variability should decrease with increasing source dimension. That is for the average stress drop Δσ ( r )∝ r H R -1, where r is the radius of a circular rupture. This result is a direct consequence of the elastic squeeze of fault asperities that induces the largest spatial fluctuations of the shear strength before and after the earthquake at local (small) scales with peculiar spatial correlations.

Constraints from fault roughness on the scale-dependent strength of rocks

Principal slip surfaces in faults have measurable roughness generated during slip. The roughness both records previous events and poses the boundary conditions for future rupture. Digital, high-precision roughness data are now available at the field scale (tens of centimeters to tens of meters) for at least 22 faults, and at the laboratory scale (millimeters to tens of centimeters) for a subset of these. We quantify the slip surface roughness by measuring the aspect ratio, which is the average asperity height divided by the profile length. Higher aspect ratios indicate rougher surfaces. From the field studies, two major trends have emerged: (1) fault surface roughness lies in a restricted range with aspect ratios in the slip-parallel direction of 0.07%–0.5% for profiles of 1 m length, and (2) fault surfaces are rougher at small scales than large ones. These features can both be interpreted as fingerprints of scale-dependent strength, which sets a limit to the aspect ratio of the surface. The measurements imply that shear strength scales with the observation scale, L , as L –0.4 . The new understanding of the physical controls on roughness allows generalization of the extant measurements of a wide array of faults.

The Burst‐Like Behavior of Aseismic Slip on a Rough Fault: The Creeping Section of the Haiyuan Fault, China

Recent observations suggesting the influence of creep on earthquakes nucleation and arrest are strong incentives to investigate the physical mechanisms controlling how active faults slip. We focus here on deriving generic characteristics of shallow creep along the Haiyuan fault, a major strike‐slip fault in China, by investigating the relationship between fault slip and geometry. We use optical images and time series of Synthetic Aperture Radar data to map the surface fault trace and the spatiotemporal distribution of surface slip along the creeping section of the Haiyuan fault. The fault trace roughness shows a power‐law behavior similar to that of the aseismic slip distribution, with a 0.8 roughness exponent, typical of a self‐affine regime. One possible interpretation is that fault geometry controls to some extent the distribution of aseismic slip, as it has been shown previously for coseismic slip along active faults. Creep is characterized by local fluctuations in rates that we define as creep bursts. The potency of creep bursts follows a power‐law behavior similar to the Gutenberg–Richter earthquake distribution, whereas the distribution of bursts velocity is non‐Gaussian, suggesting an avalanche‐like behavior of these slip events. Such similarities with earthquakes and lab experiments lead us to interpret the rich dynamics of creep bursts observed along the Haiyuan fault as resulting from long‐range elastic interactions within the heterogeneous Earth’s crust.

Frequency, pressure, and strain dependence of nonlinear elasticity in Berea Sandstone

Acoustoelasticity measurements in a sample of room dry Berea sandstone are conducted at various loading frequencies to explore the transition between the quasi‐static ( ) and dynamic (few kilohertz) nonlinear elastic response. We carry out these measurements at multiple confining pressures and perform a multivariate regression analysis to quantify the dependence of the harmonic content on strain amplitude, frequency, and pressure. The modulus softening (equivalent to the harmonic at 0f) increases by a factor 2–3 over 3 orders of magnitude increase in frequency. Harmonics at 2f, 4f, and 6f exhibit similar behaviors. In contrast, the harmonic at 1f appears frequency independent. This result corroborates previous studies showing that the nonlinear elasticity of rocks can be described with a minimum of two physical mechanisms. This study provides quantitative data that describes the rate dependency of nonlinear elasticity. These findings can be used to improve theories relating the macroscopic elastic response to microstructural features.

Flow rate dictates permeability enhancement during fluid pressure oscillations in laboratory experiments

Flow rate dictates permeability enhancement during fluid pressure oscillations in laboratory experiments Thibault Candela 1,2 , Emily E. Brodsky 2 , Chris Marone 1 , Derek Elsworth 1 1.Dept Geosciences, Penn State Univ, University Park, PA, United States. 2.Dept Earth Sciences, UC Santa Cruz, Santa Cruz, CA, United States.

The minimum scale of grooving on faults

At the field scale, nearly all fault surfaces contain grooves generated as one side of the fault slips past the other. Grooves are so common that they are one of the key indicators of principal slip surfaces. Here, we show that at sufficiently small scales, grooves do not exist on fault surfaces. A transition to isotropic roughness occurs at 4-500 mm. Although the scale of the transition can vary even between locales on a single fault, the aspect ratio of the roughness at the transition is well defined for a given fault. We interpret the transition between grooved and ungrooved scales as a transition in deformation mode of asperities on the slip surface. Grooves can form when a hard indenter slides past a softer surface. At small scales, the asperities appear to yield plastically and therefore do not generate grooves as hard indenters. The plastic yielding can be a consequence of the high shear strains required to deform the surfaces at small scales where the aspect ratio (roughness) is high. The transition to plastic yielding is predicted to occur at a specific aspect ratio for each fault, as observed. The new observation both shows a limit to one of the most commonly observed features of faults and suggests a change in the mode of failure of faults as a function of scale.

Dynamic acousto-elasticity in Berea sandstone: Influence of the strain rate

In comparison with standard nonlinear ultrasonic methods such as frequency mixing or resonance based measurements that allow one to extract average, bulk variations of modulus and attenuation versus strain level, dynamic acousto-elasticity (DAE) allows to obtain the elastic behavior over the entire dynamic cycle, detailing the full nonlinear behavior under tension and compression, including hysteresis and memory effects. To improve our understanding of these phenomena, this work aims at comparing static and dynamic acousto-elasticity to evaluate the influence of strain rate. To this purpose, we perform acousto-elasticity on a sample of Berea sandstone and a glass beads pack, oscillating them from 0.001 to 10 Hz. These results are then compared to DAE measurements made in the kHz range. We observe that the average decrease in modulus increases with frequency, meaning that conditioning effects are higher at high strain rate, when relaxation characteristic time is higher than the oscillation period. This res...