Nucleation of Stick-Slip Instability Within a Large-Scale Experimental Fault: Effects of Stress Heterogeneities Due to Loading and Gouge Layer Compaction.

Geodetic observations and large-scale laboratory experiments show that seismic instability is preceded by slow slip within a finite nucleation zone. In laboratory experiments rupture nucleation is studied mostly using bare (rock) interfaces, whereas upper crustal faults are typically filled with gouge. To investigate effects of gouge on rupture nucleation, we performed a biaxial shearing experiment on a 350 mm long saw-cut fault filled with gypsum gouge, at room temperature and a minimum horizontal stress σ 2 = 0.3-5 MPa. The gouge layer was sandwiched between polymethylmethacrylate (PMMA) plates For reference also a fault without gouge was deformed. Strain gauges and Digital Image Correlation were used to monitor the deformation field along the fault zone margins. Stick-slip behavior occurred on both the gouge-filled fault and the PMMA fault. Nucleation of instability on the PMMA fault persistently occurred from one location 2/3 to 3/4 along the fault adjacent to a slow slip zone at the fault end, but nucleation on the gouge-filled fault was more variable, nucleating at the ends and/or at approximately 2/3 along the fault, with precursory slip occurring over a large fraction of the fault. Nucleation correlated to regions of high average fault stress ratio τ/σ n , which was more variable for the gouge-filled fault due to small length scale variations in normal stress caused by heterogeneous gouge compaction. Rupture velocities and slip rates were lower for the gouge-filled fault than for the bare PMMA fault. Stick-slip persisted when σ 2 was lowered and the nucleation zone length increased, expanding from the center to the sample ends before transitioning into instability.

Investigating Compaction by Intergranular Pressure Solution Using the Discrete Element Method.

Intergranular pressure solution creep is an important deformation mechanism in the Earth's crust. The phenomenon has been frequently studied and several analytical models have been proposed that describe its constitutive behavior. These models require assumptions regarding the geometry of the aggregate and the grain size distribution in order to solve for the contact stresses and often neglect shear tractions. Furthermore, analytical models tend to overestimate experimental compaction rates at low porosities, an observation for which the underlying mechanisms remain to be elucidated. Here we present a conceptually simple, 3-D discrete element method (DEM) approach for simulating intergranular pressure solution creep that explicitly models individual grains, relaxing many of the assumptions that are required by analytical models. The DEM model is validated against experiments by direct comparison of macroscopic sample compaction rates. Furthermore, the sensitivity of the overall DEM compaction rate to the grain size and applied stress is tested. The effects of the interparticle friction and of a distributed grain size on macroscopic strain rates are subsequently investigated. Overall, we find that the DEM model is capable of reproducing realistic compaction behavior, and that the strain rates produced by the model are in good agreement with uniaxial compaction experiments. Characteristic features, such as the dependence of the strain rate on grain size and applied stress, as predicted by analytical models, are also observed in the simulations. DEM results show that interparticle friction and a distributed grain size affect the compaction rates by less than half an order of magnitude.

Misorientation distributions in hot deformed NaCl using electron backscattered diffraction.

EBSD orientation mapping has been used to derive subgrain boundary misorientation distributions in a series of hot deformed and etched NaCl samples. The main objective of this study has been to examine the influence of data processing, noise caused by angular resolution limits and step size on the subgrain misorientation distributions in hot deformed NaCl. Processing of non-indexed EBSD patterns increased the average misorientations in etched NaCl. Noise contributed significantly to low angle misorientation peaks for step sizes less than the minimum subgrain size. Orientation data collected using a step size larger than the average subgrain size cumulated misorientations across individual subgrains and effectively measured an orientation gradient between steps. Orientation gradient distributions were not influenced by noise. Average misorientation values calculated from large step data correlated well with average misorientation from small step size data, Average misorientations showed a power law relationship with strain. Three types of substructures were identified using scanning electron microscopy and EBSD mapping, equiaxed subgrains, long subgrain boundaries and a core-mantle subgrain arrangement.