Laboratory studies of frictional sliding and the implications of precursory seismicity

The dynamic transition from slow to rapid sliding along a frictional interface is of interest to geophysicists, engineers and scientists alike. In our direct shear experiment, we simulated a pre-existing frictional fault similar to those occurring naturally in the Earth. The laboratory study reported here has incorporated appropriate sensors that can detect foreshock events on the fringe of a nucleation zone prior to a gross fault rupture (main shock).;During loading we observed foreshocks sequences as slip transitioned from slow to rapid sliding. These laboratory-induced foreshocks showed similar acoustic characteristics and spatio-temporal evolution as those detected in nature. Through direct observation (video camera), foreshocks were found to be the rapid, localized (millimeter length scale) failure of highly stresses asperities formed along the interface. The interface was created by the meshing of two rough polymethyl methacrylate (PMMA) bodies in a direct shear configuration. A carefully calibrated pressure sensitive film was used to map the contact junctions (asperities) throughout the interface at a range of applied normal loads Fn. Foreshocks were found to coalesce in a region of the fault that exhibited a more dense distribution of asperities (referred to as the seismogenic region).;Microscopy of the interface in the seismogenic region displayed a variety of surface roughness at various length scales. This may have been introduced from the surface preparation techniques use to create a mature interface. The mature interface consisted of 'flat-topped' asperity regions with separating sharp valleys. The 'flat-topped' sections spanned millimetric length scales and were considerably flatter (nanometric roughness) that the roughness exhibited at longer length scales (tens of millimeters). We believe that the smoother, 'flat-topped' sections were responsible for the individual asperity formation (determining their size and strength), whereas the longer length scale roughness influenced the asperity-asperity interaction during the nucleation phase. Asperities in the seismogenic region where shown to exist close enough to each other so that elastic communication (through the off-fault material) could not be neglected.;Prior to gross fault rupture (i.e. mainshock), we measured the propagation of a slow nucleating rupture into the relatively 'locked', seimsogenic region of the fault. Slow slip dynamics were captured using slip sensors placed along the fault that measured a non-uniform slip profile leading up to failure. We found that the propagation of the slow rupture into the locked region was dependent on the normal force Fn. Higher Fn was found to slow the propagation of shear rupture into the locked region. Within the relatively 'locked' region, a noticeable increase in size and a more compact spatial-temporal distribution of foreshocks were measured when Fn was increased.;In order to develop an understanding of the relationship between Fn and the resistance of the fault to slow rupture, a quasi-static finite element (FE) model was developed. The model used distributions of asperities measured directly from the pressure sensitive film in a small section of the interface where foreshocks coalesced; specifically, the region where the slowly propagating slip front encountered the more dense distribution of asperities. A single asperity was modeled and followed the Cattaneo partial slip asperity solution. As the shear force increased along the fault, the asperities in this model were able to accommodate tangential slip by entering a partial sliding regime; the central contact of the asperities remained adhered while sliding accumulated along its periphery. Partial slip on the asperity propagated inwards as the shear force was incrementally increased. A further increase in the shear force caused the asperity to enter a full sliding condition. Increasing confining loads caused increased stiffness and increased capacity to store potential shear strain energy -- a possible measure of the 'degree of coupling' between the fault surfaces. Physics from the numerical model followed the qualitative observations made using photometry of asperities along the interface, which visualized asperities in the 'locked' region -- larger asperities remained stuck throughout the loading cycle and the light transmitted through individual asperities decreased from the periphery as shear loads increased.;The numerical partial slip, quantified by the potential energy stored by the asperity, increased relative to the normal pressure p. Asperity-asperity interactions were modeled along the interface using a quasi-static analysis. Progression of slip into the asperity field was increasingly inhibited as the normal confining force Fn was increased. The computational model provided an explanation as to why an increased confining force Fn could result in an increased resistance to slow rupture as well as an increased potential for larger foreshocks within the resistive, relatively 'locked' section of a fault. This study lays the foundation for more innovative laboratory work that could potentially improve the phenomenological models currently used to estimate earthquake hazard. (Abstract shortened by UMI.).