Search USGS⌕ Search

Geology topics

David A. Lockner

Publications and source records attributed to David A. Lockner.

At least 37 records · Page 2Linked to original sources

Data Report: Permeability, porosity, and frictional strength of core samples from IODP Expedition 366 in the Mariana forearc

Core samples from the International Ocean Discovery Program (IODP) Expedition 366 were tested in the laboratory to determine permeability, porosity, density, and frictional strength and their relation to mineralogy as part of an effort to understand hydro-mechanical processes at convergent plate margins. Seven samples were tested from a depth range of 19.6 to 197.9 m below the sea floor. The samples were derived from three serpentinite mud volcanoes in the Mariana forearc region, formed where slab-derived fluids and materials ascend along faults. The physical characteristics mirror compositional differences between predominantly serpentine-rich and saponite-rich samples. Permeability values ranged from 10-17 to 10-19 m2, low enough to facilitate the formation of high fluid pressures, which have been observed in the Mariana and other subduction megathrust environments. Porosities ranged from 0.37 to 0.51 and densities from 1.66 to 2.01 gm/cc. Serpentine-rich samples have coefficients of friction of 0.2 to 0.4, consistent with crustal serpentinite from a variety of fault zones, whereas saponite-rich samples have friction values below 0.2, consistent with saponite fault gouge from the San Andreas Fault Zone at Depth (SAFOD) drillhole in California.

Book chapter↗

KG²B, a collaborative benchmarking exercise for estimating the permeability of the Grimsel granodiorite - Part 1: measurements, pressure dependence and pore-fluid effects

Measuring the permeability of tight rocks remains a challenging task. In addition to the traditional sources of errors that affect more permeable formations (e.g. sample selection, non-representative specimens, disturbance introduced during sample acquisition and preparation), tight rocks can be particularly prone to solid–fluid interactions and thus more sensitive to the methods, procedures and techniques used to measure permeability. To address this problem, it is desirable to collect, for a single material, measurements obtained by different methods and pore-fluids. For that purpose a collaborative benchmarking exercise involving 24 laboratories was organized for measuring the permeability of a single low permeability material, the Grimsel granodiorite, at a common effective confining pressure (5 MPa). The objectives of the benchmark were: (i) to compare the results for a given method, (ii) to compare the results between different methods, (iii) to analyze the accuracy of each method, (iv) to study the influence of experimental conditions (especially the nature of pore fluid), (v) to discuss the relevance of indirect methods and models and finally (vi) to suggest good practice for low permeability measurements. In total 39 measurements were collected that allowed us to discuss the influence of (i) pore-fluid, (ii) measurement method, (iii) sample size and (iv) pressure sensitivity. Discarding some outliers from the bulk data set (4 out of 39) an average permeability of 1.11 × 10 −18 m² with a standard deviation of 0.57 × 10 −18 m² was obtained. The most striking result was the large difference in permeability for gas measurements compared to liquid measurements. Regardless of the method used, gas permeability was higher than liquid permeability by a factor approximately 2 ( k gas = 1.28 × 10 −18 m² compared to k liquid = 0.65 × 10 −18 m²). Possible explanations are that (i) liquid permeability was underestimated due to fluid-rock interactions (ii) gas permeability was overestimated due to insufficient correction for gas slippage and/or (iii) gases and liquids do not probe exactly the same porous networks. The analysis of Knudsen numbers shows that the gas permeability measurements were performed in conditions for which the Klinkenberg correction is sufficient. Smaller samples had a larger scatter of permeability values, suggesting that their volume were below the Representative Elementary Volume. The pressure dependence of permeability was studied by some of the participating teams in the range 1–30 MPa and could be fitted to an exponential law k = k o .exp(– γP eff ) with γ = 0.093 MPa −1 . Good practice rules for measuring permeability in tight materials are also provided.

Geophysical Journal International↗

KG²B, a collaborative benchmarking exercise for estimating the permeability of the Grimsel granodiorite - Part 2: modeling, microstructures and complementary data

Measuring and modelling the permeability of tight rocks remains a challenging task. In addition to the traditional sources of errors that affect more permeable formations (e.g. sample selection, non-representative specimens, disturbance introduced during sample acquisition and preparation), tight rocks can be particularly prone to solid–fluid interactions and thus more sensitive to the methods, procedures and techniques used to measure permeability. To address this problem, it is desirable to collect, for a single material, measurements obtained by different methods and pore fluids. For that purpose, a benchmarking exercise involving 24 laboratories was organized for measuring and modelling the permeability of a single low-permeability material, the Grimsel granodiorite. The objectives of the benchmark were: (i) to compare the results for a given method, (ii) to compare the results between different methods, (iii) to analyse the accuracy of each method, (iv) to study the influence of experimental conditions (especially the nature of pore fluid), (v) to discuss the relevance of indirect methods and models and finally (vi) to suggest good practice for low-permeability measurements. To complement the data set of permeability measurements presented in a companion paper, we focus here on (i) quantitative analysis of microstructures and pore size distribution, (ii) permeability modelling and (iii) complementary measurements of permeability anisotropy and poroelastic parameters. Broad ion beam—scanning electron microscopy, micro-computerized tomography, mercury injection capillary pressure (MICP) and nuclear magnetic resonance (NMR) methods were used to characterize the microstructures and provided the input parameters for permeability modelling. Several models were used: (i) basic statistical models, (ii) 3-D pore network and effective medium models, (iii) percolation model using MICP data and (iv) free-fluid model using NMR data. The models were generally successful in predicting the actual range of measured permeability. Statistical models overestimate the permeability because they do not adequately account for the heterogeneity of the crack network. Pore network and effective medium models provide additional constraints on crack parameters such as aspect ratio, aperture, density and connectivity. MICP and advanced microscopy techniques are very useful tools providing important input data for permeability estimation. Permeability measured—orthogonal to foliation is lower that—parallel to foliation. Combining the experimental and modelling results provide a unique and rich data set.

Geophysical Journal International↗

Interaction between hydraulic fracture and a preexisting fracture under triaxial stress conditions

Enhanced reservoir connectivity generally requires maximizing the intersection between hydraulic fracture (HF) and preexisting underground natural fractures (NF), while having the hydraulic fracture cross the natural fractures (and not arrest). We have studied the interaction between a hydraulic fracture and a polished saw-cut fault. The experiments include a hydraulic fracture initiating from a pressurized axial borehole (using water) that approaches a dry fault that is inclined at an angle θ with respect to the borehole axis. The experiments are conducted on Poly(methyl) Meta Acrylate (PMMA) and Solnhofen limestone, a finely grained (<5 μm grain), low permeability (<10 nD) carbonate. The confining pressure in all experiments is 5 MPa, while the differential stress (1-80 MPa) and approach angle, θ (30, 45, 60, 90°) are experimental variables. During the hydraulic fracture, acoustic emissions (AE), slip velocity, slip magnitude, stress drop and pore pressure are recorded at a 5 MHz sampling rate. A Doppler laser vibrometer measures piston velocity outside the pressure vessel to infer fault slip duration and a strain gauge adjacent to the saw-cut provides a near-field measure of axial stress. For PMMA, the coefficient of friction was 0.30 and sliding was unstable (stick-slip). The approaching HF in PMMA created a tensile fracture detected by AE transducers ~100 μs before the significant stick-slip event (45% stress drop and slip velocity of ~60 mm/s) and was arrested by the fault at all fault orientations and differential stresses, even at 90° fault orientation and 80 MPa differential stress. For Solnhofen limestone, we observed stable sliding at a coefficient of friction of 0.12. In contrast to PMMA, the HF in Solnhofen consistently crossed to the other side of the fault. When the HF crossed the fault, it produced a small stress drop (<10%) and slip velocity of only 0.5 mm/s. Theoretical models by Blanton (1986) and Renshaw and Pollard (1995) predict that HF will be arrested for Solnhofen limestone and cross PMMA 90° fault at 80 MPa differential stress. Although the exact cause for the discrepancy between experiments and the theory is not known, one feature present in the experiments but not considered in the models, is the diffusion of fluid driven by the fault slip. Thus, the formation of a "fluid-filled patch" on the fault surface as it is intersected by the HF may substantially impact the crossing/arrest behavior. The approach angle and differential stress also influence the HF initiation azimuth and breakdown pressure. In most cases, the HF initiation azimuth was normal to the fault strike. These observations suggest that the presence of natural fractures could result in rotation of hydraulic fractures to be more normal to their strike and a subsequent change in the downhole pressure recordings. The latter could be used as a diagnostic tool for predicting this interaction.

Conference Paper↗

Shear failure of a granite pin traversing a sawcut fault

Fault heterogeneities such as bumps, bends, and stepovers are commonly observed on natural faults, but are challenging to recreate under controlled laboratory conditions. We study deformation and microseismicity of a 76 mm-diameter Westerly granite cylinder with a sawcut fault with known frictional properties. An idealized asperity is added by emplacing a precision-ground 21 mm-diameter solid granite dowel that crosses the center of the fault at right angles. This intact granite ‘pin’ provides a strength contrast that resists fault slip. Upon loading to 80 MPa in a triaxial machine, we first observed a M -4 slip event that ruptured the sawcut fault, slipped 40 µm, but was halted by the granite pin. With continued loading, the pin failed in a swarm of thousands of M -6 to M -8 events known as acoustic emissions (AEs). Once the pin was fractured to a critical point, it permitted complete rupture events (M -3) on the sawcut fault (stick-slip instabilities). Subsequent slip events were preceded by clusters of foreshock-like AEs, all located on the fault plane, and the spatial extent of the foreshock clusters is consistent with our estimate of a critical nucleation dimension h*. We also identified an aseismic zone on the fault plane surrounding the fractured rock pin. A post-mortem analysis of the sample showed a thick gouge layer where the pin intersected the fault, suggesting that dilatancy of this gouge propped open the fault and prevented microseismic events in its vicinity. Recorded microseismicity separates into three categories: slip on the sawcut fault, fracture of the intact rock pin, and off-fault seismicity associated with pin-related rock joints. We found that pin fracture events were exclusively implosive (anticrack) even though the shear process zone was overall dilatant. This shows how aseismic effects can lead to unexpected seismic manifestations of certain faulting processes.

International Journal of Rock Mechanics and Mining↗

Frictional properties and 3-D stress analysis of the southern Alpine Fault, New Zealand

New Zealand's Alpine Fault (AF) ruptures quasi-periodically in large-magnitude earthquakes. Paleoseismological evidence suggests that about half of all recognized AF earthquakes terminated at the boundary between the Central and South Westland sections of the fault. There, fault geometry and the polarity of uplift change. The South Westland AF exhibits oblique-normal fault motion on a structure oriented 052°/82°SE that, for at least 35 km along strike, contains saponite-rich principal slip zone gouges. New hydrothermal friction experiments reveal that the saponite fault gouge is frictionally weak, exhibiting friction coefficients between μ = 0.12 and μ = 0.16 for a range of temperatures ( T = 25–210 °C) and effective normal stresses ( σ n ' = 31.2–93.6 MPa). The saponite gouge is rate-strengthening in all velocity steps performed at velocities between 0.01 and 3.0 μm/s, behavior conducive to aseismic creep. A three-dimensional stress analysis shows that the South Westland AF is favorably oriented with respect to the regional stress field for slip within the frictionally weak saponite fault gouge. Geometrically, the fault is severely misoriented for slip in any fault-forming materials with friction coefficients exceeding μ ∼0.5. The combination of weak gouges prone to aseismic creep, strong asperities , and low resolved shear stress may impede earthquake rupture propagation along the South Westland Alpine Fault.

Journal of Structural Geology↗

Role of fault gouge during Interaction between hydraulic fracture and a preexisting fracture

Enhanced reservoir connectivity generally requires maximizing the intersection between hydraulic fracture (HF) and preexisting underground natural fractures (NF), while having the hydraulic fracture continue to propagate across the natural fractures. Observations of downhole core samples suggest that these natural fractures are in fact veins filled with minerals such as calcite (Mighani et al., 2016). We study this interaction during the approach of a hydraulic fracture to a smooth saw-cut fracture under triaxial stress conditions. The specimen is Solnhofen limestone, a fine-grained (<5 µm grain), low permeability (<10 nD) carbonate. The differential stress (1-20 MPa) and inclination of the fault which determines the approach angle, θ (30, 60°) are the experimental variables. We conduct the experiments on both bare surface and gouge-filled fault surfaces. The gouge is a 1 mm thick crushed powder of Solnhofen limestone with <106 µm grain size. During the hydraulic fracture, acoustic emissions (AE), inferred slip velocity, axial stress and pore pressure are recorded at a 5 MHz sampling rate. The hydraulic fracture was able to cross the bare surface fault with small induced fault slip. The fault gouge increased the coefficient of friction significantly from 0.12 (bare, polished surface) to > 0.44 (gouged layer). However, the gouge-filled fault arrested the hydraulic fracture and generated a slip event with different characteristics: 1- The stress drop was larger while the generated AE signals had lower magnitude. 2- Slip velocity recorded by the vibrometer was of the same order of magnitude for the bare and gouge-filled faults, but the slip duration increased from 29 µsec for bare surface to ~2.5 msec (~90 times longer rise time) for the gouge-filled fault. The experiments suggest that the gouge-filled fault can accommodate much larger displacement while promoting slow slip on the fault which is harder to detect as AE signals. The observed long duration slip events are similar to the field observations of the long period and long duration (LPLD) events during the stimulation of clay-rich shale formations (Zoback et al., 2012). While the intrinsic low strength, high ductility, and unfavorably oriented natural fractures in shale formations are expected to reduce the occurrence of induced seismicity, our experiments suggest an additional mechanism for the observed LPLD events, i.e. the role of fault gouge. They also suggest that the microseismic detection techniques may under-predict the stimulated volume as the activation of natural gouge-filled fractures may proceed aseismically.

Conference Paper↗

Conversion of wet glass to melt at lower seismogenic zone conditions: Implications for pseudotachylyte creep

Coseismic frictional melting and the production of quenched glass called pseudotachylyte is a recurring process during earthquakes. To investigate how glassy materials affect the postseismic strength and stability of faults, obsidian gouges were sheared under dry and wet conditions from 200°C to 300°C at ~150 MPa effective normal stress. Dry glass exhibited a brittle rheology at all conditions tested, exhibiting friction values and microstructures consistent with siliciclastic materials. Likewise, wet glass at 200°C exhibited a brittle rheology. In contrast, wet gouges at 300°C transitioned from brittle sliding to linear‐viscous (Newtonian) flow at strain rates <3 × 10 −4 s −1 , indicating melt‐like behavior. The viscosity ranged from 2 × 10 11 to 7.8 × 10 11 Pa‐s. Microstructures show that viscous gouges were fully welded with rod‐shaped microlites rotated into the flow direction. Fourier transform infrared spectroscopy along with electron backscatter imaging demonstrate that hydration of the glass by diffusion of pore water was the dominant process reducing the viscosity and promoting viscous flow. As much as 5 wt % water diffused into the glass. These results may provide insight into postseismic‐slip behaviors and challenge some interpretations of fault kinematics based on studies assuming that pseudotachylyte formation and flow is solely coseismic.

Geophysical Research Letters↗

Physical properties of sidewall cores from Decatur, Illinois

To better assess the reservoir conditions influencing the induced seismicity hazard near a carbon dioxide sequestration demonstration site in Decatur, Ill., core samples from three deep drill holes were tested to determine a suite of physical properties including bulk density, porosity, permeability, Young’s modulus, Poisson’s ratio, and failure strength. Representative samples of the shale cap rock, the sandstone reservoir, and the Precambrian basement were selected for comparison. Physical properties were strongly dependent on lithology. Bulk density was inversely related to porosity, with the cap rock and basement samples being both least porous (<3 percent) and densest (~2.6 grams per cubic centimeter [g/cc]). Permeability was highest in the reservoir sandstones (10-15 to 10-18 meters squared [m2]) relative to the cap rock and basement rocks (<10-21 m2). Young’s modulus was distinctly higher in the basement rocks (45 to 80 gigapascal [GPa]) compared to the cap rock and sandstones (19 to 57 GPa). Poisson’s ratio for the sandstones varied widely (0.14 to 0.27), but the highest values were similar to the cap rock and basement rocks (0.24 to 0.28). These physical properties reflect the layered structure of the reservoir and adjacent rocks at the Decatur site. However, within the sandstone there is a great deal of lithologic variety, accounting for the large range in physical parameters for this geologic unit.

Illinois↗

Earthquake source properties from instrumented laboratory stick-slip

Stick-slip experiments were performed to determine the influence of the testing apparatus on source properties, develop methods to relate stick-slip to natural earthquakes and examine the hypothesis of McGarr [2012] that the product of stiffness, k , and slip duration, Δ t , is scale-independent and the same order as for earthquakes. The experiments use the double-direct shear geometry, Sierra White granite at 2 MPa normal stress and a remote slip rate of 0.2 µm/sec. To determine apparatus effects, disc springs were added to the loading column to vary k . Duration, slip, slip rate, and stress drop decrease with increasing k , consistent with a spring-block slider model. However, neither for the data nor model is k Δ t constant; this results from varying stiffness at fixed scale. In contrast, additional analysis of laboratory stick-slip studies from a range of standard testing apparatuses is consistent with McGarr's hypothesis. k Δ t is scale-independent, similar to that of earthquakes, equivalent to the ratio of static stress drop to average slip velocity, and similar to the ratio of shear modulus to wavespeed of rock. These properties result from conducting experiments over a range of sample sizes, using rock samples with the same elastic properties as the Earth, and scale-independent design practices.

Book chapter↗

Frictional strength of wet and dry montmorillonite

Montmorillonite is a common mineral in fault zones, and its low strength relative to other common gouge minerals is important in many models of fault rheology. However, the coefficient of friction, μ , varies with degree of saturation and is not well constrained in the literature due to the difficulty of establishing fully drained or fully dried states in the laboratory. We measured μ of both saturated and oven-dried montmorillonite at normal stresses up to 700 MPa. Care was taken to shear saturated samples slowly enough to avoid pore fluid overpressure. For saturated samples, μ increased from 0.10 to 0.28 with applied effective normal stress, while for dry samples μ decreased from 0.78 to 0.45. The steady state rate dependence of friction, ( a − b ), was positive, promoting stable sliding. The wide disparity in reported frictional strengths can be attributed to experimental procedures that promote differing degrees of partial saturation or overpressured pore fluid conditions.

Journal of Geophysical Research↗

Buried shallow fault slip from the South Napa earthquake revealed by near-field geodesy

Earthquake-related fault slip in the upper hundreds of meters of Earth’s surface has remained largely unstudied because of challenges measuring deformation in the near field of a fault rupture. We analyze centimeter-scale accuracy mobile laser scanning (MLS) data of deformed vine rows within ±300 m of the principal surface expression of the M (magnitude) 6.0 2014 South Napa earthquake. Rather than assuming surface displacement equivalence to fault slip, we invert the near-field data with a model that allows for, but does not require, the fault to be buried below the surface. The inversion maps the position on a preexisting fault plane of a slip front that terminates ~3 to 25 m below the surface coseismically and within a few hours postseismically. The lack of surface-breaching fault slip is verified by two trenches. We estimate near-surface slip ranging from ~0.5 to 1.25 m. Surface displacement can underestimate fault slip by as much as 30%. This implies that similar biases could be present in short-term geologic slip rates used in seismic hazard analyses. Along strike and downdip, we find deficits in slip: The along-strike deficit is erased after ~1 month by afterslip. We find no evidence of off-fault deformation and conclude that the downdip shallow slip deficit for this event is likely an artifact. As near-field geodetic data rapidly proliferate and will become commonplace, we suggest that analyses of near-surface fault rupture should also use more sophisticated mechanical models and subsurface geomechanical tests.

California↗

Pseudotachylyte increases the post-slip strength of faults

Solidified frictional melts, or pseudotachylytes, are observed in exhumed faults from across the seismogenic zone. These unique fault rocks, and many experimental studies, suggest that frictional melting can be an important process during earthquakes. However, it remains unknown how melting affects the post-slip strength of the fault and why many exhumed faults do not contain pseudotachylyte. Analyses of triaxial stick-slip events on Westerly Granite (Rhode Island, USA) sawcuts at confining pressures from 50 to 400 MPa show evidence for frictional heating, including some events energetic enough to generate surface melt. Total and partial stress drops were observed with slip as high as 6.5 mm. We find that in dry samples following melt-producing stick slip, the shear failure strength increased as much as 50 MPa, while wet samples had <10 MPa strengthening. Microstructural analysis indicates that the strengthening is caused by welding of the slip surface during melt quenching, suggesting that natural pseudotachylytes may also strengthen faults after earthquakes. These results predict that natural pseudotachylyte will inhibit slip reactivation and possibly generate stress heterogeneities along faults. Wet samples do not exhibit melt welding, possibly because of thermal pressurization of water reducing frictional heating during slip.

Geology↗

Gallery of melt textures developed in Westerly Granite during high-pressure triaxial friction experiments

Introduction Melting occurred during stick-slip faulting of granite blocks sheared at room-dry, room-temperature conditions in a triaxial apparatus at 200–400 megapascals (MPa) confining pressure. Petrographic examinations of melt textures focused largely on the 400-MPa run products. This report presents an overview of the petrographic data collected on those samples, followed by brief descriptions of annotated versions of all the images. Scanning electron microscope (SEM) images of the starting materials and the three examined 400-MPa samples are presented in this report. Secondary-electron (SE) and backscattered-electron (BSE) imaging techniques were used on different samples. The SE images look down on the sawcut surfaces, yielding topographic and three-dimensional textural information. The BSE imaging was done on samples cut to provide cross-sectional views of the glass-filled shear band (or zone) that developed along the sawcut. Brightness in the BSE images increases with increasing mean atomic number of the material. Additional chemical information about the quenched melt and adjoining minerals was obtained using the energy dispersive system of the SEM during BSE examinations. However, the very narrow shear-band thicknesses and common occurrence of very fine lamellar compositional layering limited the usefulness of this technique for estimating melt chemistry.

Open-File Report↗

Calibrated acoustic emission system records M -3.5 to M -8 events generated on a saw-cut granite sample

Acoustic emission (AE) analyses have been used for decades for rock mechanics testing, but because AE systems are not typically calibrated, the absolute sizes of dynamic microcrack growth and other physical processes responsible for the generation of AEs are poorly constrained. We describe a calibration technique for the AE recording system as a whole (transducers + amplifiers + digitizers + sample + loading frame) that uses the impact of a 4.76-mm free-falling steel ball bearing as a reference source. We demonstrate the technique on a 76-mm diameter cylinder of westerly granite loaded in a triaxial deformation apparatus at 40 MPa confining pressure. The ball bearing is dropped inside a cavity within the sample while inside the pressure vessel. We compare this reference source to conventional AEs generated during loading of a saw-cut fault in a second granite sample. All located AEs occur on the saw-cut surface and have moment magnitudes ranging from M −5.7 down to at least M −8. Dynamic events rupturing the entire simulated fault surface (stick–slip events) have measurable stress drop and macroscopic slip and radiate seismic waves similar to those from a M −3.5 earthquake. The largest AE events that do not rupture the entire fault are M −5.7. For these events, we also estimate the corner frequency (200–300 kHz), and we assume the Brune model to estimate source dimensions of 4–6 mm. These AE sources are larger than the 0.2 mm grain size and smaller than the 76 × 152 mm fault surface.

Rock Mechanics and Rock Engineering↗

Evolution of fracture permeability of ultramafic rocks undergoing serpentinization at hydrothermal conditions: An experimental study

We performed flow-through laboratory experiments on five cylindrically cored samples of ultramafic rocks, in which we generated a well-mated through-going tensile fracture, to investigate evolution of fracture permeability during serpentinization. The samples were tested in a triaxial loading machine at a confining pressure of 50 MPa, pore pressure of 20 MPa, and temperature of 260&deg;C, simulating a depth of 2 km under hydrostatic conditions. A pore pressure difference of up to 2 MPa was imposed across the ends of the sample. Fracture permeability decreased by 1&ndash;2 orders of magnitude during the 200&ndash;330 h experiments. Electron microprobe and SEM data indicated the formation of needle-shaped crystals of serpentine composition along the walls of the fracture, and chemical analyses of sampled pore fluids were consistent with dissolution of ferro-magnesian minerals. By comparing the difference between fracture permeability and matrix permeability measured on intact samples of the same rock types, we concluded that the contribution of the low matrix permeability to flow is negligible and essentially all of the flow is focused in the tensile fracture. The experimental results suggest that the fracture network in long-lived hydrothermal circulation systems can be sealed rapidly as a result of mineral precipitation, and generation of new permeability resulting from a combination of tectonic and crystallization-induced stresses is required to maintain fluid circulation.

Geochemistry, Geophysics, Geosystems↗

Hydrothermal frictional strengths of rock and mineral samples relevant to the creeping section of the San Andreas Fault

We compare frictional strengths in the temperature range 25–250 °C of fault gouge from SAFOD (CDZ and SDZ) with quartzofeldspathic wall rocks typical of the central creeping section of the San Andreas Fault (Great Valley sequence and Franciscan Complex). The Great Valley and Franciscan samples have coefficients of friction, μ > 0.35 at all experimental conditions. Strength is unchanged between 25° and 150 °C, but μ increases at higher temperatures, exceeding 0.50 at 250 °C. Both samples are velocity strengthening at room temperature but show velocity-weakening behavior beginning at 150 °C and stick-slip motion at 250 °C. These rocks, therefore, have the potential for unstable seismic slip at depth. The CDZ gouge, with a high saponite content, is weak (μ = 0.09–0.17) and velocity strengthening in all experiments, and μ decreases at temperatures above 150 °C. Behavior of the SDZ is intermediate between the CDZ and wall rocks: μ < 0.2 and does not vary with temperature. Although saponite is probably not stable at depths greater than ∼3 km, substitution of the frictionally similar minerals talc and Mg-rich chlorite for saponite at higher temperatures could potentially extend the range of low strength and stable slip down to the base of the seismogenic zone.

California↗