Geology topics

T.E. Tullis

Publications and source records attributed to T.E. Tullis.

7 recordsLinked to original sources

USGS2018

On the depth extent of co-seismic rupture

We investigate the implications of deformation experiments for the coseismic down‐dip extent of rupture in quasi‐dynamic, whole‐cycle earthquake models of a fault for which the depth of the transition between seismic and aseisimic fault slip depends on strain rate. The calculations use a dislocation fault model from Tse and Rice (1986) with a vertical strike‐slip orientation, mode III rupture, and variable along‐strike length. Our reference calculation is the original rheological representation of Tse and Rice with a strain‐rate‐independent transition. The primary calculations use two different representations of a strain‐rate‐dependent transition: (1) between rate‐weakening friction and dislocation creep and (2) between rate‐weakening and rate‐strengthening frictions. For both these cases, when fault strength is high (friction between 0.5 and 0.6) and the transition is sharp, coseismic slip extends a small distance (1–2 km) below the fixed temperature (depth) that is commonly used to define the rheological transition at the plate‐motion rate. Thus, coseismic slip occurs below the depth assumed in seismic hazard models using microseismicity or a chosen fixed‐temperature contour. Though significant coseismic slip occurs below the plate‐rate transition depth, the added moment is < 10 % "> < 10 % of the total. The deep extension is a region that is rheologically distinct; for example, deep coseismic slip can produce a stress increase rather than a stress drop. If friction is smaller, the deepening effect and its contribution to moment are larger. For all representations of the transition, average and surface slip increase with the along‐strike rupture length in a manner consistent with the limited data from natural observations. However, this property is not controlled by the assumed fault rheology; instead, it arises because the stiffness of the fault decreases weakly with fault length, an intrinsic and unrealistic property of the particular crustal scale fault model used.

Bulletin of the Seismological Society of America↗

USGS2008

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↗

USGS1998

Effects of slip, slip rate, and shear heating on the friction of granite

The stability of fault slip is sensitive to the way in which frictional strength responds to changes in slip rate and in particular to the effective velocity dependence of steady state friction Δμ ss /Δ ln V . This quantity can vary substantially with displacement, temperature and slip rate. To investigate the physical basis for this behavior and the possible influence of shear heating, we slid initially bare granite surfaces in unconfined rotary shear to displacements of hundreds of millimeters at normal stresses, σ n of 10 and 25 MPa and at room temperature. We imposed step changes in slip rate within the range 10 −2 to 10 3.5 μm/s and also monitored frictional heating with thermistors embedded in the granite. The transient response of μ to slip rate steps was fit to a rate- and state-dependent friction law using two state variables to estimate the values of several parameters in the constitutive law. The first 20 mm of slip shows rising friction and falling Δμ ss /Δ ln V ; further slip shows roughly constant friction, Δμ ss /Δ ln V and parameter values, suggesting that a steady state condition is reached on the fault surface. At V ≤ 10 μm/s, Δμ ss /Δ ln V = −0.004 ± 0.001. At higher rates the response is sensitive to normal stress: At σ n = 25 MPa granite shows a transition to effective velocity strengthening (Δμ ss /Δ ln V = 0.008 ± 0.004) at the highest slip rates tested. At 10 MPa granite shows a less dramatic change to Δμ ss /Δ ln V ≈ 0 at the highest rates. The maximum temperature measured in the granite is ∼60°C at 25 MPa and 10 3.5 μm/s. Temperatures are in general agreement with a numerical model of heat conduction which assumes spatially homogeneous frictional heating over the sliding surface. The simplest interpretation of our measurements of Δμ ss /Δ ln V is that the granite is inherently velocity weakening (∂μ ss /∂ ln V < 0) and temperature strengthening (∂μ ss /∂ T −1 < 0) at all velocities. At high slip rates the response of μ to changes in temperature from shear heating may outweigh the response to changing velocity, such that the net effect Δμ ss /Δ ln V > 0 mimics velocity strengthening. These results have implications for the frictional behavior of faults during earthquakes. High slip rates may cause a switch to effective velocity strengthening which could limit peak coseismic slip rate and stress drop. For fluid-saturated faults, strengthening by this mechanism may be partly or fully offset by weakening due to thermal pressurization of a poorly drained pore fluid.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

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↗

USGS1996

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↗

USGS1996

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↗

USGS1995

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↗