Role of void space geometry in permeability evolution in crustal rocks at elevated pressure

[1] A key consequence of the presence of void space within rock is their significant influence upon fluid transport properties. In this study, we measure changes in elastic wave velocities (P and S) contemporaneously with changes in permeability and porosity at elevated pressure for three rock types with widely different void space geometries: a high‐porosity sandstone (Bentheim), a tight sandstone (Crab Orchard), and a microcracked granodiorite (Takidani). Laboratory data are then used with the permeability models of Guéguen and Dienes and Kozeny‐Carman to investigate the characteristics that different void space geometries impart to measured permeabilities. Using the Kachanov effective medium theory, elastic wave velocities are inverted, permitting the recovery of crack density evolution with increasing effective pressure. The crack densities are then used as input to the microcrack permeability model of Guéguen and Dienes. The classic Kozeny‐Carman approach of Walsh and Brace is also applied to the measured permeability data via a least squares fit in order to extract tortuosity data. We successfully predict the evolution of permeability with increasing effective pressure, as directly measured in experiments, and report the contrast between permeability changes observed in rock where microcracks or equant pores dominate the microstructure. Additionally, we show how these properties are affected by anisotropy of the rock types via the measured anisotropic fabrics in each rock. The combined experimental and modeling results illustrate the importance of understanding the details of how rock microstructure changes in response to an external stimulus in predicting the simultaneous evolution of different rock physical properties

Dimensional Reduction without Extra Continuous Dimensions

We describe a novel approach to dimensional reduction in classical field theory. Inspired by ideas from noncommutative geometry, we introduce extended algebras of differential forms over space-time, generalized exterior derivatives and generalized connections associated with the "geometry" of space-times with discrete extra dimensions. We apply our formalism to theories of gauge- and gravitational fields and find natural geometrical origins for an axion- and a dilaton field, as well as a Higgs field.Comment: 23 page

Dynamic weakening and amorphization in serpentinite during laboratory earthquakes

The mechanical properties of serpentinites are key factors in our understanding of the dynamics of earthquake ruptures in subduction zones, especially intermediate-depth earthquakes. Here, we performed shear rupture experiments on natural antigorite serpentinite, which showed that friction reaches near-zero values during spontaneous dynamic rupture propagation. Rapid coseismic slip (>1 m/s), although it occurs over short distances (<1 mm), induces significant overheating of microscale asperities along the sliding surface, sufficient to produce surface amorphization and likely some melting. Antigorite dehydration occurs in the fault walls, which leaves a partially amorphized material. The water generated potentially contributes to the production of a low-viscosity pressurized melt, explaining the near-zero dynamic friction levels observed in some events. The rapid and dramatic dynamic weakening in serpentinite might be a key process facilitating the propagation of earthquakes at intermediate depths in subduction zones

Modeling the permeability evolution of microcracked rocks from elastic wave velocity inversion at elevated isostatic pressure

[1] A key consequence of the presence of microcracks within rock is their significant influence upon elastic anisotropy and transport properties. Here two rock types (a basalt and a granite) with contrasting microstructures, dominated by microcracks, have been investigated using an advanced experimental arrangement capable of measuring porosity, P wave velocity, S wave velocity, and permeability contemporaneously at effective pressures up to 100 MPa. Using the Kachanov (1994) noninteractive effective medium theory, the measured elastic wave velocities are inverted using a least squares fit, permitting the recovery of the evolution of crack density and aspect ratio with increasing isostatic pressure. Overall, the agreement between measured and predicted velocities is good, with average error less than 0.05 km/s. At larger scales and above the percolation threshold, macroscopic fluid flow also depends on the crack density and aspect ratio. Using the permeability model of Guéguen and Dienes (1989) and the crack density and aspect ratio recovered from the elastic wave velocity inversion, we successfully predict the evolution of permeability with pressure for direct comparison with the laboratory measurements. We also calculate the evolution of the crack porosity with increasing isostatic pressure, on the basis of the calculated crack density, and compare this directly with the experimentally measured porosity. These combined experimental and modeling results illustrate the importance of understanding the details of how rock microstructures change in response to an external stimulus when predicting the simultaneous evolution of rock physical properties

Modeling the permeability evolution of microcracked rocks from elastic wave velocity inversion at elevated isostatic pressure

[1] A key consequence of the presence of microcracks within rock is their significant influence upon elastic anisotropy and transport properties. Here two rock types (a basalt and a granite) with contrasting microstructures, dominated by microcracks, have been investigated using an advanced experimental arrangement capable of measuring porosity, P wave velocity, S wave velocity, and permeability contemporaneously at effective pressures up to 100 MPa. Using the Kachanov (1994) noninteractive effective medium theory, the measured elastic wave velocities are inverted using a least squares fit, permitting the recovery of the evolution of crack density and aspect ratio with increasing isostatic pressure. Overall, the agreement between measured and predicted velocities is good, with average error less than 0.05 km/s. At larger scales and above the percolation threshold, macroscopic fluid flow also depends on the crack density and aspect ratio. Using the permeability model of Guéguen and Dienes (1989) and the crack density and aspect ratio recovered from the elastic wave velocity inversion, we successfully predict the evolution of permeability with pressure for direct comparison with the laboratory measurements. We also calculate the evolution of the crack porosity with increasing isostatic pressure, on the basis of the calculated crack density, and compare this directly with the experimentally measured porosity. These combined experimental and modeling results illustrate the importance of understanding the details of how rock microstructures change in response to an external stimulus when predicting the simultaneous evolution of rock physical properties

Nucleation of Laboratory Earthquakes: Quantitative Analysis and Scalings

In this study we use the precursory acoustic emission (AE) activity during the nucleation of stick-slip instability as a proxy to investigate foreshock occurrence prior to natural earthquakes. We report on three stick-slip experiments performed on cylindrical samples of Indian metagabbro under upper crustal stress conditions (30–60 MPa). AEs were continuously recorded by eight calibrated acoustic sensors during the experiments. Seismological parameters (moment magnitude, corner frequency and stress-drop) of the detected AEs (−8.8 ≤ Mw ≤ −7) follow the scaling law between moment magnitude and corner frequency that characterizes natural earthquakes. AE activity always increases toward failure and is driven by along fault slip velocity. The stacked AE foreshock sequences follow an inverse Omori type law, with a characteristic Omori time c inversely proportional to normal stress. AEs moment magnitudes increase toward failure, as manifested by a decrease in b-value from ∼1 to ∼0.5 at the end of the nucleation process. During nucleation, foreshocks migrate toward the mainshock epicenter location, and stabilize at a distance from the latter compatible with the predicted Rate-and-State nucleation size. Importantly, the nucleation characteristic timescale also scales inversely with applied normal stress and the expected nucleation size. Finally, we infer that foreshocks are the byproducts of the nucleation phase which is an almost fully aseismic process. Nevertheless, the seismic/aseismic energy release ratio continuously increases during nucleation, highlighting that, the nucleation process starts as a fully aseismic process, and evolves toward a cascading process at the onset of dynamic rupture

Frictional Heating Processes and Energy Budget During Laboratory Earthquakes

International audienceDuring 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

Laboratory micro-seismic signature of shear faulting and fault slip in shale

This article reports the results of a triaxial deformation experiment conducted on a transversely isotropic shale specimen. This specimen was instrumented with ultrasonic transducers to monitor the evolution of the micro-seismic activity induced by shear faulting (triaxial failure) and subsequent fault slip at two different rates. The strain data demonstrate the anisotropy of the mechanical (quasi-static) compliance of the shale; the P-wave velocity data demonstrate the anisotropy of the elastic (dynamic) compliance of the shale. The spatio-temporal evolution of the micro-seismic activity suggests the development of two distinct but overlapping shear faults, a feature similar to relay ramps observed in large-scale structural geology. The shear faulting of the shale specimen appears quasi-aseismic, at least in the 0.5 MHz range of sensitivity of the ultrasonic transducers used in the experiment. Concomitantly, the rate of micro-seismic activity is strongly correlated with the imposed slip rate and the evolution of the axial stress. The moment tensor inversion of the focal mechanism of the high quality micro-seismic events recorded suggests a transition from a non-shear dominated to a shear dominated micro-seismic activity when the rock evolves from initial failure to larger and faster slip along the fault. The frictional behaviour of the shear faults highlights the possible interactions between small asperities and slow slip of a velocity-strengthening fault, which could be considered as a realistic experimental analogue of natural observations of non-volcanic tremors and (very) low-frequency earthquakes triggered by slow slip events

Fracture and damage localization in volcanic edifice rocks from El Hierro, Stromboli and Tenerife

© 2018 The Author(s). We present elastic wave velocity and strength data from a suite of three volcanic rocks taken from the volcanic edifices of El Hierro and Tenerife (Canary Islands, Spain), and Stromboli (Aeolian Islands, Italy). These rocks span a range of porosity and are taken from volcanoes that suffer from edifice instability. We measure elastic wave velocities at known incident angles to the generated through-going fault as a function of imposed strain, and examine the effect of the damage zone on P-wave velocity. Such data are important as field measurements of elastic wave tomography are key tools for understanding volcanic regions, yet hidden fractures are likely to have a significant effect on elastic wave velocity. We then use elastic wave velocity evolution to calculate concomitant crack density evolution which ranges from 0 to 0.17: highest values were correlated to the damage zone in rocks with the highest initial porosity