Frictional Heating Processes and Energy Budget During Laboratory Earthquakes

During an earthquake, part of the released elastic strain energy is dissipated within the slip zone by frictional and fracturing processes, the rest being radiated away via elastic waves. While frictional heating plays a key role in the energy budget of earthquakes, it could not be resolved by seismological data up to now. Here we investigate the dynamics of laboratory earthquakes by measuring frictional heat dissipated during the propagation of shear instabilities at stress conditions typical of seismogenic depths. We estimate the complete energy budget of earthquake rupture and demonstrate that the radiation efficiency increases with thermal-frictional weakening. Using carbon properties and Raman spectroscopy, we map spatial heat heterogeneities on the fault surface. We show that an increase in fault strength corresponds to a transition from a weak fault with multiple strong asperities and little overall radiation, to a highly radiative fault behaving as a single strong asperity. Plain Language Summary In nature, earthquakes occur when the stress accumulated in a medium is released by frictional sliding on faults. The stress released is dissipated into fracture and heat energy or radiated through seismic waves. The seismic efficiency of an earthquake is a measure of the fraction of the energy that is radiated away into the host medium. Because faults are at inaccessible depths, we reproduce earthquakes in the laboratory under natural in situ conditions to understand the physical processes leading to dynamic rupture. We estimate the first complete energy budget of an earthquake and show that increasing heat dissipation on the fault increases the radiation efficiency. We develop a novel method to illuminate areas of the fault that get excessively heated up. We finally introduce the concept of spontaneously developing heat asperities, playing a major role in the radiation of seismic waves during an earthquake.

On the scale dependence in the dynamics of frictional rupture: Constant fracture energy versus size-dependent breakdown work

Fabian Barras, Mathias Alexandre David Lebihain, Federica Paglialunga, François Xavier Thibault Passelègue, Marie Estelle Solange Violay

Potential energy stored during the inter-seismic period by tectonic loading around faults is released during earthquakes as radiated energy, frictional dissipation and fracture energy. The latter is of first importance since it is expected to control the nucleation, the propagation and the arrest of the seismic rupture. On one side, the seismological fracture energy estimated for natural earthquakes (commonly called breakdown work) ranges between 1 J/m2and tens of MJ/m2for the largest events, and shows a clear slip dependence. On the other side, recent experimental studies highlighted that, concerning rupture experiments, fracture energy is a material property (energy required to break the fault interface) independently of the size of the event, i.e. of the seismic slip. To reconcile these contradictory observations and definitions, we performed stick-slip experiments, as analog for earthquakes, in a bi-axial shear configuration. We estimated fracture energy through both Linear Elastic Fracture Mechanics (LEFM) and a Cohesive Zone Model (CZM) and through the integration of the near-fault stress-slip evolution. We show that, at the scale of our experiments, fault weakening is divided into a near-tip weakening, corresponding to an energy of few J/m2, consistent with the one estimated through LEFM and CZM, and a long-tailed weakening corresponding to a larger energy not localized at the rupture tip, increasing with slip. Through numerical simulations, we demonstrate that only near-tip weakening controls the rupture initiation and that long-tailed weakening can enhance slip during rupture propagation and allow the rupture to overcome stress heterogeneity along the fault. We conclude that the origin of the seismological estimates of breakdown work could be related to the energy dissipated in the long-tailed weakening rather than to the one dissipated near the tip.

2022

Can Precursory Moment Release Scale With Earthquake Magnitude? A View From the Laboratory

Mateo Alejandro Acosta, François Xavier Thibault Passelègue, Marie Estelle Solange Violay

Today, earthquake precursors remain debated. While precursory slow slip is an important feature of earthquake nucleation, foreshock sequences are not always observed, and their temporal evolution remains poorly constrained. We report on laboratory earthquakes conducted under upper-crustal stress and fluid pressure conditions. The dynamics of precursors (slip, seismicity, and fault coupling) prior to the mainshock are dramatically affected by slight changes in fault conditions (fluid pressure and slip history). A relationship between precursory moment release and mainshock magnitude is systematically observed, independent of fault conditions. Based on nucleation theory, we derive a semiempirical scaling relationship which explains this trend for laboratory earthquakes. Several natural observations of earthquakes ranging from similar to M-w 6.0-9.0, where precursory moment release could be estimated, seem to follow the extrapolation of the laboratory-derived scaling law. Notwithstanding spatiotemporal complexity in natural seismicity, some moderate to large earthquake magnitudes might be estimated through integrated seismological and geodetic measurements of both seismic and aseismic slips during nucleation. Plain Language Summary Understanding the preparation phase that precedes earthquake ruptures (nucleation) is crucial for seismic hazard assessment because it might provide information on the impending earthquake. Here, we show that the temporal evolution of laboratory earthquake precursors (precursory slow slip, precursory seismicity, and fault coupling) is of little use in assessing an impending earthquake's size. Nevertheless, independent of fault slip behavior (seismic or aseismic) and environmental conditions (stress state and fluid pressure and slip history), the amount of moment released during the preparation phase scales with the earthquake's magnitude. This property is demonstrated by laboratory observations and earthquake nucleation theory and seems compatible with several natural observations of earthquakes ranging from M-w 6.0 to M-w 9.0. As a consequence, if earthquakes exhibit a preparation phase, it could be possible that this phase is larger or longer for higher magnitude earthquakes and consequently, more likely to be detectable.

2019

Energy budget of laboratory earthquakes: a comparison between linear elastic fracture mechanics approach and experimental approach

Federica Paglialunga, François Xavier Thibault Passelègue, Marie Estelle Solange Violay

Earthquakes correspond to a sudden release of elastic energy stored during inter-seismic period by tectonic loading around fault. The earthquake energy budget consists of four non-independent terms: the energy release rate (by unit crack length), the fracture energy, the heat term and finally the radiated energy. These terms depend on the rupture and sliding velocities, the amount of slip and the stress drop. Because of the impossibility to access to stress and strain conditions at depth, the earthquake energy budget cannot be fully constrained from seismological data, limiting our understanding of its influence on rupture propagation. To address this issue, we conducted stick-slip experiments with large samples in a biaxial configuration apparatus. By imposing constant normal load and increasing shear load, seismic events were produced on a 20 cm long fault, for which the energy budget was estimated using different methods. Fracture energy was estimated by recording the strain field around the crack tip through high frequency (2 MHz) strain gage rosettes array and comparing it to the theoretical LEFM strain field predictions obtained for same conditions (i.e. rupture velocity, distance from the fault). Fracture energies were then inverted and found to range in between 1 and 10 J/m2. At the same time the energy partitioning was estimated through stress-slip evolution during rupture. The fracture energies obtained from this method are almost one order of magnitude larger than the ones inverted from LEFM and range in between 1 and 90 J/m2. Moreover, the energy partitioning shows the radiated energy ranging between 80 and 300 J/m2 and finally the heat/thermal energy as the largest fraction of the energy partitioning with values ranging from 200 to 2500 J/m2. Our preliminary results highlight the importance of understanding the contribution of heat energy in frictional processes, since this term cannot be estimated from seismological data.

2020