Geology ReportsSearch

Geology topics

N.M. Beeler

Publications and source records attributed to N.M. Beeler.

15 recordsLinked to original sources

Effective stress, friction and deep crustal faulting

Studies of crustal faulting and rock friction invariably assume the effective normal stress that determines fault shear resistance during frictional sliding is the applied normal stress minus the pore pressure. Here we propose an expression for the effective stress coefficient α f at temperatures and stresses near the brittle-ductile transition (BDT) that depends on the percentage of solid-solid contact area across the fault. α f varies with depth and is only near 1 when the yield strength of asperity contacts greatly exceeds the applied normal stress. For a vertical strike-slip quartz fault zone at hydrostatic pore pressure and assuming 1 mm and 1 km shear zone widths for friction and ductile shear, respectively, the BDT is at ~13 km. α f near 1 is restricted to depths where the shear zone is narrow. Below the BDT α f = 0 is due to a dramatically decreased strain rate. Under these circumstances friction cannot be reactivated below the BDT by increasing the pore pressure alone and requires localization. If pore pressure increases and the fault localizes back to 1 mm, then brittle behavior can occur to a depth of around 35 km. The interdependencies among effective stress, contact-scale strain rate, and pore pressure allow estimates of the conditions necessary for deep low-frequency seismicity seen on the San Andreas near Parkfield and in some subduction zones. Among the implications are that shear in the region separating shallow earthquakes and deep low-frequency seismicity is distributed and that the deeper zone involves both elevated pore fluid pressure and localization.

Journal of Geophysical Research B: Solid Earth

Constructing constitutive relationships for seismic and aseismic fault slip

For the purpose of modeling natural fault slip, a useful result from an experimental fault mechanics study would be a physically-based constitutive relation that well characterizes all the relevant observations. This report describes an approach for constructing such equations. Where possible the construction intends to identify or, at least, attribute physical processes and contact scale physics to the observations such that the resulting relations can be extrapolated in conditions and scale between the laboratory and the Earth. The approach is developed as an alternative but is based on Ruina (1983) and is illustrated initially by constructing a couple of relations from that study. In addition, two example constitutive relationships are constructed; these describe laboratory observations not well-modeled by Ruina's equations: the unexpected shear-induced weakening of silica-rich rocks at high slip speed (Goldsby and Tullis, 2002) and fault strength in the brittle ductile transition zone (Shimamoto, 1986). The examples, provided as illustration, may also be useful for quantitative modeling.

Pure and Applied Geophysics

Constitutive relationships and physical basis of fault strength due to flash heating

We develop a model of fault strength loss resulting from phase change at asperity contacts due to flash heating that considers a distribution of contact sizes and nonsteady state evolution of fault strength with displacement. Laboratory faulting experiments conducted at high sliding velocities, which show dramatic strength reduction below the threshold for bulk melting, are well fit by the model. The predicted slip speed for the onset of weakening is in the range of 0.05 to 2 m/s, qualitatively consistent with the limited published observations. For this model, earthquake stress drops and effective shear fracture energy should be linearly pressure-dependent, whereas the onset speed may be pressure-independent or weakly pressure-dependent. On the basis of the theory, flash weakening is expected to produce large dynamic stress drops, small effective shear fracture energy, and undershoot. Estimates of the threshold slip speed, stress drop, and fracture energy are uncertain due to poor knowledge of the average ontact dimension, shear zone thickness and gouge particle size at seismogenic depths. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth

Attempting to bridge the gap between laboratory and seismic estimates of fracture energy

To investigate the behavior of the fracture energy associated with expanding the rupture zone of an earthquake, we have used the results of a large-scale, biaxial stick-slip friction experiment to set the parameters of an equivalent dynamic rupture model. This model is determined by matching the fault slip, the static stress drop and the apparent stress. After confirming that the fracture energy associated with this model earthquake is in reasonable agreement with corresponding laboratory values, we can use it to determine fracture energies for earthquakes as functions of stress drop, rupture velocity and fault slip. If we take account of the state of stress at seismogenic depths, the model extrapolation to larger fault slips yields fracture energies that agree with independent estimates by others based on dynamic rupture models for large earthquakes. For fixed stress drop and rupture speed, the fracture energy scales linearly with fault slip.

Geophysical Research Letters

Stress-induced, time-dependent fracture closure at hydrothermal conditions

Time-dependent closure of fractures in quartz was measured in situ at 22–530°C temperature and 0.1-150 MPa water pressure. Unlike previous crack healing and rock permeability studies, in this study, fracture aperture is monitored directly and continuously using a windowed pressure vessel, a long-working-distance microscope, and reflected-light interferometry. Thus the fracture volume and geometry can be measured as a function of time, temperature, and water pressure. Relatively uniform closure occurs rapidly at temperatures and pressures where quartz becomes significantly soluble in water. During closure the aperture is reduced by as much as 80% in a few hours. We infer that this closure results from the dissolution of small particles or asperities that prop the fracture open. The driving force for closure via dissolution of the prop is the sum of three chemical potential terms: (1) the dissolution potential, proportional to the logarithm of the degree of undersaturation of the solution; (2) the coarsening potential, proportional to the radius of curvature of the prop; and (3) the pressure solution potential, proportional to the effective normal stress at the contact between propping particles and the fracture wall. Our observations suggest that closure is controlled by a pressure solution-like process. The aperture of dilatant fractures and microcracks in the Earth that are similar to those in our experiments, such as ones generated from thermal stressing or brittle failure during earthquake rupture and slip, will decrease rapidly with time, especially if the macroscopic stress is nonhydrostatic.

Journal of Geophysical Research B: Solid Earth

On the expected relationships among apparent stress, static stress drop, effective shear fracture energy, and efficiency

We consider expected relationships between apparent stress τ a and static stress drop Δτ s using a standard energy balance and find τ a = Δτ s (0.5 - ξ), where ξ is stress overshoot. A simple implementation of this balance is to assume overshoot is constant; then apparent stress should vary linearly with stress drop, consistent with spectral theories (Brune, 1970) and dynamic crack models (Madariaga, 1976). Normalizing this expression by the static stress drop defines an efficiency η sw = τ s a/Δτ s as follows from Savage and Wood (1971). We use this measure of efficiency to analyze data from one of a number of observational studies that find apparent stress to increase with seismic moment, namely earthquakes recorded in the Cajon Pass borehole by Abercrombie (1995). Increases in apparent stress with event size could reflect an increase in seismic efficiency; however, η sw for the Cajon earthquakes shows no such increase and is approximately constant over the entire moment range. Thus, apparent stress and stress drop co-vary, as expected from the energy balance at constant overshoot. The median value of η sw for the Cajon earthquakes is four times lower than η sw for laboratory events. Thus, these Cajon-recorded earthquakes have relatively low and approximately constant efficiency. As the energy balance requires η sw = 0.5 - ξ, overshoot can be estimated directly from the Savage-Wood efficiency; overshoot is positive for Cajon Pass earthquakes. Variations in apparent stress with seismic moment for these earthquakes result primarily from systematic variations in static stress drop with seismic moment and do not require a relative decrease in sliding resistance with increasing event size (dynamic weakening). Based on the comparison of field and lab determinations of the Savage-Wood efficiency, we suggest the criterion η sw > 0.3 as a test for dynamic weakening in excess of that seen in the lab.

Bulletin of the Seismological Society of America

Earthquake stress drop and laboratory-inferred interseismic strength recovery

We determine the scaling relationships between earthquake stress drop and recurrence interval t r that are implied by laboratory-measured fault strength. We assume that repeating earthquakes can be simulated by stick-slip sliding using a spring and slider block model. Simulations with static/kinetic strength, time-dependent strength, and rate- and state-variable-dependent strength indicate that the relationship between loading velocity and recurrence interval can be adequately described by the power law V L ∝ t r n where n≈−1. Deviations from n=−1 arise from second order effects on strength, with n>−1 corresponding to apparent time-dependent strengthening and n<−1 corresponding to weakening. Simulations with rate and state-variable equations show that dynamic shear stress drop Δτ d scales with recurrence as d Δτ d / d ln t r ≤σ e ( b-a ), where σ e is the effective normal stress, μ=τ/σ e , and ( a-b )= d μ ss / d ln V is the steady-state slip rate dependence of strength. In addition, accounting for seismic energy radiation, we suggest that the static shear stress drop Δτ s scales as d Δτ s / d ln t r ≤σ e (1 +ζ)( b-a ), where ζ is the fractional overshoot. The variation of Δτ s with ln t r for earthquake stress drops is somewhat larger than implied by room temperature laboratory values of ζ and b-a . However, the uncertainty associated with the seismic data is large and the discrepancy between the seismic observations and the rate of strengthening predicted by room temperature experiments is less than an order of magnitude.

Journal of Geophysical Research B: Solid Earth

Stress drop with constant, scale independent seismic efficiency and overshoot

To model dissipated and radiated energy during earthquake stress drop, I calculate dynamic fault slip using a single degree of freedom spring-slider block and a laboratory-based static/kinetic fault strength relation with a dynamic stress drop proportional to effective normal stress. The model is scaled to earthquake size assuming a circular rupture; stiffness varies inversely with rupture radius, and rupture duration is proportional to radius. Calculated seismic efficiency, the ratio of radiated to total energy expended during stress drop, is in good agreement with laboratory and field observations. Predicted overshoot, a measure of how much the static stress drop exceeds the dynamic stress drop, is higher than previously published laboratory and seismic observations and fully elasto-dynamic calculations. Seismic efficiency and overshoot are constant, independent of normal stress and scale. Calculated variation of apparent stress with seismic moment resembles the observational constraints of McGarr [1999].

Geophysical Research Letters

Pore fluid pressure, apparent friction, and Coulomb failure

Many recent studies of stress-triggered seismicity rely on a fault failure model with a single free parameter, the apparent coefficient of friction, presumed to be a material constant with possible values 0 ≤ μ′ ≤ 1. These studies may present a misleading view of fault strength and the role of pore fluid pressure in earthquake failure. The parameter μ′ is intended to incorporate the effects of both friction and pore pressure, but is a material constant only if changes in pore fluid pressure induced by changes in stress are proportional to the normal stress change across the potential failure plane. Although specific models of fault zones permit such a relation, neither is it known that fault zones within the Earth behave this way, nor is this behavior expected in all cases. In contrast, for an isotropic homogeneous poroelastic model the pore pressure changes are proportional to changes in mean stress, μ′ is not a material constant, and −∞ ≤ μ′ ≤ +∞. Analysis of the change in Coulomb failure stress for tectonically loaded reverse and strike-slip faults shows considerable differences between these two pore pressure models, suggesting that such models might be distinguished from one another using observations of triggered seismicity (e.g., aftershocks). We conclude that using the constant apparent friction model exclusively in studies of Coulomb failure stress is unwise and could lead to significant errors in estimated stress change and seismic hazard.

Journal of Geophysical Research B: Solid Earth

Earthquake triggering by transient and static deformations

Observational evidence for both static and transient near-field and far-field triggered seismicity are explained in terms of a frictional instability model, based on a single degree of freedom spring-slider system and rate- and state-dependent frictional constitutive equations. In this study a triggered earthquake is one whose failure time has been advanced by Δ t (clock advance) due to a stress perturbation. Triggering stress perturbations considered include square-wave transients and step functions, analogous to seismic waves and coseismic static stress changes, respectively. Perturbations are superimposed on a constant background stressing rate which represents the tectonic stressing rate. The normal stress is assumed to be constant. Approximate, closed-form solutions of the rate-and-state equations are derived for these triggering and background loads, building on the work of Dieterich [1992, 1994]. These solutions can be used to simulate the effects of static and transient stresses as a function of amplitude, onset time t 0 , and in the case of square waves, duration. The accuracies of the approximate closed-form solutions are also evaluated with respect to the full numerical solution and t 0 . The approximate solutions underpredict the full solutions, although the difference decreases as t 0 approaches the end of the earthquake cycle. The relationship between Δ t and t 0 differs for transient and static loads: a static stress step imposed late in the cycle causes less clock advance than an equal step imposed earlier, whereas a later applied transient causes greater clock advance than an equal one imposed earlier. For equal Δ t , transient amplitudes must be greater than static loads by factors of several tens to hundreds depending on t 0 . We show that the rate-and-state model requires that the total slip at failure is a constant, regardless of the loading history. Thus a static load applied early in the cycle, or a transient applied at any time, reduces the stress at the initiation of failure, whereas static loads that are applied sufficiently late raise it. Rate-and-state friction predictions differ markedly from those based on Coulomb failure stress changes (ΔCFS) in which Δ t equals the amplitude of the static stress change divided by the background stressing rate. The ΔCFS model assumes a stress failure threshold, while the rate-and-state equations require a slip failure threshold. The complete rate-and-state equations predict larger Δ t than the ΔCFS model does for static stress steps at small t 0 , and smaller Δ t than the ΔCFS model for stress steps at large t 0 . The ΔCFS model predicts nonzero Δ t only for transient loads that raise the stress to failure stress levels during the transient. In contrast, the rate-and-state model predicts nonzero Δ t for smaller loads, and triggered failure may occur well after the transient is finished. We consider heuristically the effects of triggering on a population of faults, as these effects might be evident in seismicity data. Triggering is manifest as an initial increase in seismicity rate that may be followed by a quiescence or by a return to the background rate. Available seismicity data are insufficient to discriminate whether triggered earthquakes are “new” or clock advanced. However, if triggering indeed results from advancing the failure time of inevitable earthquakes, then our modeling suggests that a quiescence always follows transient triggering and that the duration of increased seismicity also cannot exceed the duration of a triggering transient load. Quiescence follows static triggering only if the population of available faults is finite.

Journal of Geophysical Research B: Solid Earth

The roles of time and displacement in velocity-dependent volumetric strain of fault zones

The relationship between measured friction μ A and volumetric strain during frictional sliding was determined using a rate and state variable dependent friction constitutive equation, a common work balance relating friction and volume change, and two types of experimental faults: initially bare surfaces of Westerly granite and rock surfaces separated by a 1 mm layer of less than 90 μm Westerly granite gouge. The constitutive equation is the sum of a constant term representing the nominal resistance to sliding and two smaller terms: a rate dependent term representing the shear viscosity of the fault surface (direct effect), and a term which represents variations in the area of contact (evolution effect). The work balance relationship requires that μ A differs from the frictional resistance that leads to shear heating by the derivative of fault normal displacement with respect shear displacement, d δ n / d δ s . An implication of this relationship is that the rate dependence of d δ n / d δ s contributes to the rate dependence of μ A . Experiments show changes in sliding velocity lead to changes in both fault strength and volume. Analysis of data with the rate and state equations combined with the work balance relationship preclude the conventional interpretation of the direct effect in the rate and state variable constitutive equations. Consideration of a model bare surface fault consisting of an undeformable indentor sliding on a deformable surface reveals a serious flaw in the work balance relationship if volume change is time-dependent. For the model, at zero slip rate indentation creep under the normal load leads to time-dependent strengthening of the fault surface but, according to the work balance relationship, no work is done because compaction or dilatancy can only be induced by shearing. Additional tests on initially bare surfaces and gouges show that fault normal strain in experiments is time-dependent, consistent with the model. This time-dependent fault normal strain, which is not accounted for in the work balance relationship, explains the inconsistency between the constitutive equations and the work balance. For initially bare surface faults, all rate dependence of volume change is due to time dependence. Similar results are found for gouge. We conclude that μ A reflects the frictional resistance that results in shear heating, and no correction needs to be made for the volume changes. The result that time-dependent volume changes do not contribute to μ A is a general result and extends beyond these experiments, the simple indentor model and particular constitutive equations used to illustrate the principle.

Journal of Geophysical Research B: Solid Earth

Transient triggering of near and distant earthquakes

We demonstrate qualitatively that frictional instability theory provides a context for understanding how earthquakes may be triggered by transient loads associated with seismic waves from near and distance earthquakes. We assume that earthquake triggering is a stick-slip process and test two hypotheses about the effect of transients on the timing of instabilities using a simple spring-slider model and a rate- and state-dependent friction constitutive law. A critical triggering threshold is implicit in such a model formulation.

Bulletin of the Seismological Society of America

Self-healing slip pulses in dynamic rupture models due to velocity-dependent strength

Seismological observations of short slip duration on faults (short rise time on seismograms) during earthquakes are not consistent with conventional crack models of dynamic rupture and fault slip. In these models, the leading edge of rupture stops only when a strong region is encountered, and slip at an interior point ceases only when waves from the stopped edge of slip propagate back to that point. In contrast, some seismological evidence suggests that the duration of slip is too short for waves to propagate from the nearest edge of the ruptured surface, perhaps even if the distance used is an asperity size instead of the entire rupture dimension. What controls slip duration, if not dimensions of the fault or of asperities?

Bulletin of the Seismological Society of America

Frictional behavior of large displacement experimental faults

The coefficient of friction and velocity dependence of friction of initially bare surfaces and 1-mm-thick simulated fault gouges (<90 μm) of Westerly granite were determined as a function of displacement to >400 mm at 25°C and 25 MPa normal stress. Steady state negative friction velocity dependence and a steady state fault zone microstructure are achieved after ∼18 mm displacement, and an approximately constant strength is reached after a few tens of millimeters of sliding on initially bare surfaces. Simulated fault gouges show a large but systematic variation of friction, velocity dependence of friction, dilatancy, and degree of localization with displacement. At short displacement (<10 mm), simulated gouge is strong, velocity strengthening and changes in sliding velocity are accompanied by relatively large changes in dilatancy rate. With continued displacement, simulated gouges become progressively weaker and less velocity strengthening, the velocity dependence of dilatancy rate decreases, and deformation becomes localized into a narrow basal shear which at its most localized is observed to be velocity weakening. With subsequent displacement, the fault restrengthens, returns to velocity strengthening, or to velocity neutral, the velocity dependence of dilatancy rate becomes larger, and deformation becomes distributed. Correlation of friction, velocity dependence of friction and of dilatancy rate, and degree of localization at all displacements in simulated gouge suggest that all quantities are interrelated. The observations do not distinguish the independent variables but suggest that the degree of localization is controlled by the fault strength, not by the friction velocity dependence. The friction velocity dependence and velocity dependence of dilatancy rate can be used as qualitative measures of the degree of localization in simulated gouge, in agreement with previous studies. Theory equating the friction velocity dependence of simulated gouge to the sum of the friction velocity dependence of bare surfaces and the velocity dependence of dilatancy rate of simulated gouge fails to quantitatively account for the experimental observations.

Journal of Geophysical Research B: Solid Earth

Implications of Coulomb plasticity for the velocity dependence of experimental faults

Simulated fault gouges often deform more stably than initially bare surfaces of the same composition. It is important to understand why the sliding stability is enhanced because the presence of gouge on natural faults may have the same effect as seen in experiments, and thus explain the absence of earthquakes at shallow depths. Gouge stabilization in experiments has been attributed to positive contributions to velocity dependence within gouge layers from either dilation (Maroneet al., 1990) or grain fracture (Biegelet al., 1989). In this study we test the hypothesis that some aspects of gouge and initially bare surface velocity dependence are identical by measuring the time-dependent constitutive parameter b. An important result follows however from stress analysis: if both sample configurations are frictional in the Mohr-Coulomb sense, each configuration is required to deform on planes of distinctly different orientation. The measured strength and velocity dependence will reflect this geometric difference. Our observed values of b for simulated granite and quartz gouge are two to two and a half times smaller than b for initially bare surfaces. This difference is completely accounted for if gouge is represented as a cohesionless-Coulomb plastic material. The analysis demonstrates the following points: 1) gouge deformation is fully consistent with Coulomb plasticity, 2) observed gouge velocity dependence is a function of observed strength and 3) the constitutive parameter b is the same for both bare surfaces and gouge. Furthermore, the results suggest that there is no time-dependent strengthening associated with stabilizing effects in gouge. These observations provide a framework for understanding how slip on initially bare surfaces and gouge deformation are related. ?? 1995 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH