Search

The Japan Trench Fast Drilling Project (JFAST) endeavored to establish the stress state on the shallow subduction megathrust that slipped during the M9 Tohoku earthquake. Borehole breakout data from the drill hole can constrain both the orientation and magnitude of the principal stresses. Here we reanalyze those data to refine our understanding of the stress state on the fault. In particular, we (1) improve the identification of breakouts, (2) consider a fuller range of stress states consistent with the data, and (3) incorporate new and more robust laboratory constraints on rock strength. The original conclusion that the region is in a normal faulting regime after the earthquake is strengthened by the new analysis. The combined analysis suggests that the earthquake released sufficient elastic strain energy to reset the local stress field. [ABSTRACT FROM AUTHOR]

Earthquakes are the result of slip along faults and are due to the decrease of rock frictional strength (dynamic weakening) with increasing slip and slip rate. Friction experiments simulating the abrupt accelerations (>>10 m/s2), slip rates (∼1 m/s), and normal stresses (>>10 MPa) expected at the passage of the earthquake rupture along the front of fault patches, measured large fault dynamic weakening for slip rates larger than a critical velocity of 0.01-0.1 m/s. The dynamic weakening corresponds to a decrease of the friction coefficient (defined as the ratio of shear stress vs. normal stress) up to 40%-50% after few millimetres of slip (flash weakening), almost independently of rock type. The microstructural evolution of the sliding interfaces with slip may yield hints on the microphysical processes responsible for flash weakening. At the microscopic scale, the frictional strength results from the interaction of micro- to nano-scale surface irregularities (asperities) which deform during fault sliding. During flash weakening, the visco-plastic and brittle work on the asperities results in abrupt frictional heating (flash heating) and grain size reduction associated with mechano-chemical reactions (e.g., decarbonation in CO2-bearing minerals such as calcite and dolomite; dehydration in water-bearing minerals such as clays, serpentine, etc.) and phase transitions (e.g., flash melting in silicate-bearing rocks). However, flash weakening is also associated with grain size reduction down to the nanoscale. Using focused ion beam scanning and transmission electron microscopy, we studied the micro-physical mechanisms associated with flash heating and nanograin formation in carbonate-bearing fault rocks. Experiments were conducted on pre-cut Carrara marble (99.9% calcite) cylinders using a rotary shear apparatus at conditions relevant to seismic rupture propagation. Flash heating and weakening in calcite-bearing rocks is associated with a shock-like stress release due to the migration of fast-moving dislocations and the conversion of their kinetic energy into heat. From a review of the current natural and experimental observations we speculate that this mechanism tested for calcite-bearing rocks, is a general mechanism operating during flash weakening (e.g., also precursory to flash melting in the case of silicate-bearing rocks) for all fault rock types undergoing fast slip acceleration due to the passage of the seismic rupture front. [ABSTRACT FROM AUTHOR]

Mapping full-field displacement and strain changes on the Earth's surface following an earthquake is of paramount importance to enhance our understanding of earthquake mechanics. Currently, aerial and satellite images taken pre- and post-earthquake can be processed with sub-pixel correlation algorithms to infer the co-seismic ground deformations (e.g., [, ]). However, the interpretation of this data is not straightforward due to the inherent complexity of natural faults and deformation fields. To gain understanding into rupture mechanics and to help interpret complex rupture features occurring in nature, we develop a laboratory earthquake setup capable of reproducing displacement and strain maps similar to those obtained in the field, while maintaining enough simplicity so that clear conclusions can be drawn. Earthquakes are mimicked in the laboratory by dynamic rupture propagating along an inclined frictional interface formed by two Homalite plates under compression (e.g., []). In our study, the interface is partially glued, in order to confine the rupture before it reaches the ends of the specimen. The specimens are painted with a speckle pattern to provide the surface with characteristic features for image matching. Images of the specimens are taken before and after dynamic rupture with a 4 Megapixels resolution CCD camera. The digital images are analyzed with two software packages for sub-pixel correlation: VIC-2D (Correlated Solutions Inc.) and COSI-Corr []. Both VIC-2D and COSI-Corr are able to characterize the full-field static displacement of the experimentally produced dynamic shear ruptures. The correlation analysis performed with either software clearly shows (i) the relative displacement (slip) along the frictional interface, (ii) the rupture arrest on the glued boundaries, and (iii) the presence of wing cracks. The obtained displacement measurements are converted to strains, using non-local de-noising techniques; stresses are obtained by introducing Homalite's constitutive properties. This study is a first step towards using the digital image correlation method in combination with high-speed photography to capture the highly transient phenomena involved in dynamic rupture. [ABSTRACT FROM AUTHOR]

Earthquakes are the result of slip along faults and are due to the decrease of rock frictional strength (dynamic weakening) with increasing slip and slip rate. Friction experiments simulating the abrupt accelerations (>>10 m/s2), slip rates (~1 m/s), and normal stresses (>>10 MPa) expected at the passage of the earthquake rupture along the front of fault patches, measured large fault dynamic weakening for slip rates larger than a critical velocity of 0.01–0.1 m/s. The dynamic weakening corresponds to a decrease of the friction coefficient (defined as the ratio of shear stress vs. normal stress) up to 40%–50% after few millimetres of slip (flash weakening), almost independently of rock type. The microstructural evolution of the sliding interfaces with slip may yield hints on the microphysical processes responsible for flash weakening. At the microscopic scale, the frictional strength results from the interaction of micro- to nano-scale surface irregularities (asperities) which deform during fault sliding. During flash weakening, the visco-plastic and brittle work on the asperities results in abrupt frictional heating (flash heating) and grain size reduction associated with mechano-chemical reactions (e.g., decarbonation in CO2-bearing minerals such as calcite and dolomite; dehydration in water-bearing minerals such as clays, serpentine, etc.) and phase transitions (e.g., flash melting in silicate-bearing rocks). However, flash weakening is also associated with grain size reduction down to the nanoscale. Using focused ion beam scanning and transmission electron microscopy, we studied the micro-physical mechanisms associated with flash heating and nanograin formation in carbonate-bearing fault rocks. Experiments were conducted on pre-cut Carrara marble (99.9% calcite) cylinders using a rotary shear apparatus at conditions relevant to seismic rupture propagation. Flash heating and weakening in calcite-bearing rocks is associated with a shock-like stress release due to the migration of fast-moving dislocations and the conversion of their kinetic energy into heat. From a review of the current natural and experimental observations we speculate that this mechanism tested for calcite-bearing rocks, is a general mechanism operating during flash weakening (e.g., also precursory to flash melting in the case of silicate-bearing rocks) for all fault rock types undergoing fast slip acceleration due to the passage of the seismic rupture front.

Abstract: During fluid injection experiments at the geothermal site of Soultz-sous-Forêts (France), more than 114,000 induced seismic events with magnitudes between −2.0 and +2.9 were detected by a local downhole monitoring network. Of these, 35,039 events are sufficiently constrained to be located. Hypocenters align along a sub-vertical, planar structure with the apparent width being dominated by data scattering indicating that seismic activity predominantly occurs along a (pre-existing) larger scale fault structure. For this scenario, we present a numerical model to simulate hydraulic overpressures and induced seismicity during hydraulic injection. The numerical model is based on the physical processes of fluid pressure and stress diffusion with triggering of the induced seismicity being controlled by Coulomb friction. Even in its simplest form of a fault zone without any structural heterogeneity, the numerical model reproduces typical observations at Soultz-sous-Forêts, such as number and magnitude of induced events, hypocenter locations (including the Kaiser effect), occurrence of post-injection seismicity, and the largest magnitude event occurring several days after shut-in. [Copyright &y& Elsevier]

SUMMARY A model of a ductile fault zone is incorporated into a forward 3-D earthquake model to better constrain fault-zone hydraulics. The conceptual framework of the model fault zone was chosen such that two distinct parts are recognized. The fault core, characterized by a relatively low permeability, is composed of a coseismic fault surface embedded in a visco-elastic volume that can creep and compact. The fault core is surrounded by, and mostly sealed from, a high permeability damaged zone. The model fault properties correspond explicitly to those of the coseismic fault core. Porosity and pore pressure evolve to account for the viscous compaction of the fault core, while stresses evolve in response to the applied tectonic loading and to shear creep of the fault itself. A small diffusive leakage is allowed in and out of the fault zone. Coseismically, porosity is created to account for frictional dilatancy. We show in the case of a 3-D fault model with no in-plane flow and constant fluid compressibility, pore pressures do not drop to hydrostatic levels after a seismic rupture, leading to an overpressured weak fault. Since pore pressure plays a key role in the fault behaviour, we investigate coseismic hydraulic property changes. In the full 3-D model, pore pressures vary instantaneously by the poroelastic effect during the propagation of the rupture. Once the stress state stabilizes, pore pressures are incrementally redistributed in the failed patch. We show that the significant effect of pressure-dependent fluid compressibility in the no in-plane flow case becomes a secondary effect when the other spatial dimensions are considered because in-plane flow with a near-lithostatically pressured neighbourhood equilibrates at a pressure much higher than hydrostatic levels, forming persistent high-pressure fluid compartments. If the observed faults are not all overpressured and weak, other mechanisms, not included in this model, must be at work in nature, which need to be investigated. Significant leakage perpendicular to the fault strike (in the case of a young fault), or cracks hydraulically linking the fault core to the damaged zone (for a mature fault) are probable mechanisms for keeping the faults strong and might play a significant role in modulating fault pore pressures. Therefore, fault-normal hydraulic properties of fault zones should be a future focus of field and numerical experiments. [ABSTRACT FROM AUTHOR]

Earthquake mechanics may be determined by the geometry of a fault system. Slip on a fractal branching fault surface can explain: 1) regeneration of stress irregularities in an earthquake; 2) the concentration of stress drop in an earthquake into asperities; 3) starting and stopping of earthquake slip at fault junctions, and 4) self-similar scaling of earthquakes. Slip at fault junctions provides a natural realization of barrier and asperity models without appealing to variations of fault strength. Fault systems are observed to have a branching fractal structure, and slip may occur at many fault junctions in an earthquake. Consider the mechanics of slip at one fault junction. In order to avoid a stress singularity of order 1/r, an intersection of faults must be a triple junction and the Burgers vectors on the three fault segments at the junction must sum to zero. In other words, to lowest order the deformation consists of rigid block displacement, which ensures that the local stress due to the dislocations is zero. The elastic dislocation solution, however, ignores the fact that the configuration of the blocks changes at the scale of the displacement. A volume change occurs at the junction; either a void opens or intense local deformation is required to avoid material overlap. The volume change is proportional to the product of the slip increment and the total slip since the formation of the junction. Energy absorbed at the junction, equal to confining pressure times the volume change, is not large enongh to prevent slip at a new junction. The ratio of energy absorbed at a new junction to elastic energy released in an earthquake is no larger than P/µ where P is confining pressure and µ is the shear modulus. At a depth of 10 km this dimensionless ratio has th value P/µ= 0.01. As slip accumulates at a fault junction in a number of earthquakes, the fault segments are displaced such that they no longer meet at a single point. For this reason the volume increment for a given slip increment becomes larger. A juction with past accumulated slip ??0 is a strong barrier to earthquakes with maximum slip um < 2 (P/µ) u0 = u0/50. As slip continues to occur elsewhere in the fault system, a stress concentration will grow at the old junction. A fresh fracture may occur in the stress concentration, establishing a new triple junction, and allowing continuity of slip in the fault system. The fresh fracture could provide the instability needed to explain earthquakes. Perhaps a small fraction (on the order of P/µ) of the surface that slips in any earthquake is fresh fracture. Stress drop occurs only on this small fraction of the rupture surface, the asperities. Strain change in the asperities is on the order of P/µ. Therefore this model predicts average strais change in an earthquake to be on the order of (P/µ)2 = 0.0001, as is observed.