Geology topics

Nicholas M. Beeler

Publications and source records attributed to Nicholas M. Beeler.

At least 19 recordsLinked to original sources

USGS2025

Rapid fault healing from cementation controls the dynamics of deep slow slip and tremor

Despite its status as one of the most important discoveries in geophysics, the physical mechanism(s) responsible for slow slip events (SSEs) are not well understood. Here, we synthesize observations of deep SSEs in the Cascadia Subduction Zone and argue that rapid, cohesive fault strengthening may control the dynamics of deep SSEs. Cohesive strength is frequently ignored in constitutive laws used to describe fault rheology in numerical simulations of earthquakes and SSEs alike. To demonstrate its importance, we perform and analyze a suite of petrological experiments that simulate fault healing under representative pressure and temperature conditions. We show that significant cohesive strength recovery caused by dissolution-precipitation processes occurs on timescales of just a few hours. Together, our experimental and observational results support the idea that cohesion is a key component of fault strength under SSE conditions and highlight the need for its inclusion in both future experiments and numerical models of fault slip.

British Columbia, Washington↗

USGS2025

To heal or not to heal?: 1. The effect of pore fluid pressure on the frictional healing behavior of lithologies in Oklahoma

The competition between fault healing (i.e., re-strengthening) and fault loading determines the timing and magnitude of fault failure within the seismic cycle. Repeating earthquakes can give observational estimates of fault healing rates, however, it is difficult to link laboratory studies of frictional healing and observed healing rates from repeating earthquakes in part because of uncertainty in lithology at depth. Due to well-constrained and relatively simple geology, earthquakes in Oklahoma can be linked to the granitic basement rock and to the Arbuckle Group, which is primarily composed of dolomite at earthquake depths. Here, we conduct friction experiments to measure healing rates of the two earthquake-bearing lithologies at confining pressures representative of earthquake depths and pore pressures ranging from 0% to 80% of the confining pressure. We measure frictional healing by executing slide-hold-slide tests with hold times ranging from 3 s to 3000 s. The friction experiments on the Troy Granite indicate that pore fluid pressure does not greatly affect healing rate. On the other hand, the dolomite of the Arbuckle Group exhibits decreased healing with increased pore fluid pressure, with weakening at the highest pore pressure. We hypothesize that this is due to an increase in dissolution of dolomite at high pore pressures/low effective normal stress. These healing rates are used in the companion paper to understand the moment-recurrence time behavior of repeating earthquakes in Prague, Oklahoma. This work has implications for possible enhanced dissolution and weakening behavior of the Arbuckle Group during wastewater injection activities.

Oklahoma↗

USGS2025

Earthquake stress drop, source spectra from laboratory experiments, and the apparent unsteadiness of frictional slip

Over a range of co-seismic total energy release, high sample rate slip, stress, and slip velocity measurements were made during frictional failure to determine seismic source properties. The near-field measurements were used to estimate stress drops using spectral approaches that are analogous to the standard seismological moment-corner determined Δσ and acceleration amplitude Δσ arms stress drops. Rapid frictional sliding produces slip velocity amplitude spectra that vary with frequency -1 so long as the energy dissipated by friction approaches the total energy release. Under these conditions Δσ and Δσ arms are similar in size. Since these are nearly direct measurements of fault slip velocity and because they produce flat acceleration spectra between the corner and an apparent fmax, friction in these experiments seems to produce dynamic acceleration and deceleration at all resolved frequencies. The physical origin of this unsteady motion is not known exactly but likely reflects rapid dynamic variations in fault strength or applied stress. On the other hand, when the total energy stored prior to failure is artificially made somewhat larger relative to that which can be dissipated by friction the velocity spectra vary with frequency -2 immediately below the corner. At these conditions Δσ and Δσ arms diverge. The excess stored energy is preferentially partitioned into radiated energy upon failure. Collectively the experiments imply that to produce unsteady, white noise accelerations that are observed for natural earthquakes, requires that at least 95% of the energy released does not escape the source region to be radiated to the far-field.

Bulletin of the Seismological Society of America↗

USGS2025

Probing lower-crustal fault properties with frequency-dependent tidal tremor triggering

The way seismicity responds to periodic stress perturbations offers crucial insights into the processes that can trigger an earthquake. Laboratory and theoretical analyses have shown that the period of imposed forcing and source properties affect the sensitivity to triggering, but frequency-dependent triggering of tectonic faults is poorly understood. The rate of low-frequency earthquakes (LFEs) near Parkfield, California has been found to be strongly correlated with solid earth tides. Tidal forcing acts over multiple frequencies, and the sensitivity to tidal triggering of LFEs therefore provides a unique opportunity to probe the physics of earthquake triggering and underlying fault properties. Here, we find that the response of LFEs to solid earth tides at diurnal and semi-diurnal frequencies is highly variable but spatially coherent along the San Andreas Fault. Using rate-state friction modeling, we find that the variation of the amplitude of tidal modulation is mainly affected by the spatial variation of the background effective normal stress, whereas the spatially varied tidal modulation at diurnal and semi-diurnal frequencies is mainly affected by the fault frictional property. The spatial variations of the LFE response reveal a heterogeneous lower-crustal geologic structure and complex physical faulting processes below the rupture zone of eventual great earthquakes.

California↗

USGS2023

Time-dependent weakening of granite at hydrothermal conditions

The evolution of a fault's frictional strength during the interseismic period is a critical component of the earthquake cycle, yet there have been relatively few studies that examine the time-dependent evolution of strength at conditions representative of seismogenic depths. Using a simulated fault in Westerly granite, we examined how frictional strength evolves under hydrothermal conditions up to 250°C during slide-hold-slide experiments. At temperatures ≤100°C, frictional strength generally increases with hold duration but, at 200 and 250°C, an initial increase in strength transitions to rapid time-dependent weakening for holds longer than 14 hr. Forward modeling of long hold periods at 250°C using the rate and state friction constitutive equations requires a second, strongly negative, state variable with a long evolution distance. This implies that significant hydrothermal alteration is occurring at 250°C, consistent with microstructural observations of dissolution and secondary mineral precipitation.

Geophysical Research Letters↗

USGS2023

Strength recovery in quartzite is controlled by changes in friction in experiments at hydrothermal conditions up to 200°C

The rate of fault zone restrengthening between earthquakes can be influenced by both frictional and cohesive healing processes. Friction is dependent on effective normal stress while cohesion is independent of normal stress, potentially explaining—in part—the lack of depth dependence of earthquake stress drops. Although amenable to laboratory testing, few studies have systematically addressed the normal stress dependence of restrengthening rate. This is partially due to difficulty in separating relative contributions of friction and cohesion in recovery of fault strength. We present results from a series of slide-hold-slide tests on thin layers (≤10 𝜇m) of ultrafine quartz gouge that develop during shearing of initially bare-surface quartzite. Tests were conducted at 10 MPa constant pore pressure, 20–200 MPa constant effective normal stress, and temperatures of 22°–200°C. Restrengthening, defined as the difference between peak shear stress measured after resumption of sliding and steady-state sliding shear stress, increases with the log of hold duration. The 200°C healing rate, 0.014 per e-fold increase in time, is comparable to that determined from seismological observations along the Calaveras Fault, California. Construction of Mohr-Coulomb failure envelopes shows that changes in cohesion are small (<1 MPa) and independent of hold durations to 10 5 s, indicating that the increased strength is due to changes in the friction coefficient. These experimental results are inconsistent with the hypothesis that cohesive healing explains the depth independence of earthquake stress drop, but higher temperatures, longer time-scales, and more complex mineralogy could facilitate cohesive healing in natural fault systems.

Journal of Geophysical Research - Solid Earth↗

USGS2023

On the scale-dependence of fault surface roughness

Defining roughness as the ratio of height to length, the standard approach to characterize amplitudes of single fault, joint and fracture surfaces is to measure average height as a function of profile length. Empirically, this roughness depends strongly on scale. The ratio is approximately 0.01 at a few mm but 10× smaller at a few tens of meters. Surfaces are rougher at small scales. However, these conclusions are metric-dependent. If instead height is averaged over wavelength, roughness is nearly Brown spatial noise, having almost scale-independent apparent surface height to wavelength ratio. The small deviation from scale-independence is of the opposite sense than found using the standard metric; surfaces are slightly rougher at long wavelengths. Some natural surfaces may be Brownian within the measurement uncertainties. These contradictions are curiosities of surfaces that have Hurst exponents between 0.5 and 1, as natural fault surfaces do. The wavelength-based analysis of roughness and how it changes with scale are straight-forward; a normalized Fourier transform approximately preserves amplitude and its scale dependence in the wavelength domain. Among the conclusions from reconsideration of scale dependence are that the scale dependence is weak and much smaller than that of other fault and shear zone properties. Background and aftershock seismicity, jogs and step-overs indicate strong localization (smoothing) with slip and scale. The lack of strong scale dependence to surface roughness suggests it is not the dominant control on brittle shear zone evolution.

JGR Solid Earth↗

USGS2023

Effect of thermal and mechanical processes on hydraulic transmissivity evolution

Fracture healing is a critical component of enhanced geothermal systems, the earthquake cycle, and induced seismicity. Accordingly, there is significant interest in understanding the process of healing and its effects on fluid transport. The creation, reactivation, and sustainability of fracture networks depend on complex coupling among thermal, hydraulic, mechanical, and chemical processes. We use laboratory slide-hold-slide experiments, at temperatures from 22 to 200 ˚C, to examine effects of fracture reactivation and quasi-static loading on the evolution of fluid transport properties of simulated fractures in Westerly granite. At all temperatures, the in-plane hydraulic transmissivity consistently decays during hold periods resulting in an overall reduction in transmissivity. During the first three to fifteen hours of an experiment, transmissivity decreases rapidly due to the generation of wear products, development of a sliding surface, and compaction of the resulting gouge. Once the sliding surface has developed, the long-term transmissivity decay rate at 22 and 100 ˚C is significantly lower than the transmissivity decay rate during the initial 3 to 15 hours of the experiment. However, at 200 °C, the decay of hydraulic transmissivity remains high throughout the experiment. The long-term decay of hydraulic transmissivity can be fitted with a power law model with more rapid reduction of hydraulic transmissivity at higher temperature. Periods of sliding on the fracture surface result in transient increases in the transmissivity, due to shear dilation, as is expected for Coulomb materials. These transients are superimposed on the long-term decay. When sliding ceases and a new hold period commences, there is a rapid reduction in transmissivity and return to the long-term rate of transmissivity decay. The rate of decay of the transmissivity transients is inversely proportional to temperature, in contrast to the long-term decay and the expected behavior for processes like subcritical crack growth and indentation creep. The higher decay rates that are observed during the initial 3-15 hours of the tests and following sliding, are associated with times that the porosity of the gouge is expected to be high. The difference in decay rates suggests that when the gouge is driven far from equilibrium by active shearing, densification may be dominated by a different mechanism from long-term compaction.

Conference Paper↗

USGS2022

Brittle faulting at elevated temperature and vanishing effective stress

If brittle fault strength depends only on friction, slip instability is discouraged at low effective normal stress, σ . Stress drop and the critical stiffness necessary for unstable sliding both vanish with σ ; small earthquakes cannot occur. Very low σ is inferred in the source region of low-frequency earthquakes (LFEs) on the San Andreas fault (SAF). Moreover, if pore pressure, p , is undrained at low σ , then instabilities are prevented at all scales. This is due to dilatant strengthening which arises due to a dependence of porosity on strain rate. Dilatant strengthening is σ -independent and dominates at low σ . Undrained p is inferred over time scales of less than a few days for the SAF LFEs. Based on experiments that measure rapid contact overgrowth between 350 and 530°C at very low σ , fault failure controlled by time-dependent cementation is invoked as an explanation for the SAF LFEs. Because this “cohesion” is σ -independent, stress drops can occur at σ = 0. If in addition cohesion exceeds any dilatant strengthening during slip, cohesion dominates strength at low σ . Dilatancy measured in prior faulting and shear experiments indicate that at all stress levels steady-state porosity depends on σ in addition to strain rate. Moreover, porosity at low σ depends elastically on the confining and differential stresses. A model with these additional pore pressure effects, friction, and time-dependent cohesion, applied to the SAF LFEs produces stress drops, slip speeds, and durations that are consistent with the observations, when the shear-induced dilatancy is not extreme.

JGR Solid Earth↗

USGS2022

The evolution of rock friction is more sensitive to slip than elapsed time, even at near-zero slip rates

For many decades, frictional strength increase at low slip rates has been ascribed to time-dependent contact-area growth across the sliding interface. As a result, phenomenological models that correctly predict contact-area growth, as observed in laboratory experiments, have also been widely assumed to be appropriate descriptors of frictional strength evolution. We present experiments that impose more than 5-orders-of-magnitude slip-rate reductions on granite to show that frictional strength evolution in these rocks unequivocally refutes such models. Instead, the data suggest that, even at subnanometric slip rates, frictional strength dominantly evolves with accrued slip. This remarkable slip-sensitivity of friction requires changes of intrinsic strength of the interface with slip that are absent from popular conceptual models of friction at the microscopic contact scale.

Proceedings of the National Academy of Sciences↗

USGS2022

Apparent age dependence of the fault weakening distance in rock friction

During rock friction experiments at large displacement, room temperature and humidity, and following a hold test, the fracture energy increases approximately as the square of the logarithm of hold duration. While it's been long known that failure strength increases with log hold time, here the slip weakening distance, d h , also increases. The weakening distance increase is large, hundreds of percent change over a few thousand seconds. The initial bare surface and simulated fault gouge experiments were conducted in rotary shear at 25 MPa normal stress, 21 MPa confining stress and at displacements greater than 100 mm. In contrast, initially bare surface experiments at 5 MPa normal stress, unconfined at displacements less than 10 mm show effectively no change in d h . We attribute the difference to the presence of an appreciable shear zone that develops due to wear over significant displacements, confined at elevated normal stress. Prior published studies of sheared simulated fault gouge at short displacement show both acknowledged and unacknowledged increases in d h that may relate to our observations. Since natural faults have well-developed shear zones, the observations have more direct relevance to earthquake nucleation than prior laboratory studies that use short displacement data and focus on frictional strength recovery alone. However, the physics underlying this increase in weakening distance are not known; candidates are compaction (Nakatani, 1998) and delocalization (Sleep et al., 2000). Additional caveats are that these are room temperature and humidity experiments, at a single normal stress that have not yet been reproduced in other laboratories.

Journal of Geophysical Research↗

USGS2021

Characterizing fault roughness—Are faults rougher at long or short wavelengths?

Changes in fault roughness with scale, “scaling,” is the topic of this report; changes are considered using a general power law relation between some measure of surface height, H , and another of length, L , H=kL n , where k is a constant and n is an exponent that characterizes the scaling. Extensive profile measurements of natural fault surfaces show that the ratio of average surface height to profile length decreases with scale. Average height is defined using the root mean squared height, Rq. For this analysis, fault surfaces are smoother at long wavelengths (have smaller average height to profile length ratios) than they are at shorter wavelengths. These and other statistical properties of natural fault surfaces hold for more than five orders of magnitude, a huge range from tens of micrometers to 10 meters. However, a different roughness metric, the average height (amplitude) that is specifically associated with a wavelength shows the opposite sense of scaling. The ratio of average amplitude to wavelength increases with wavelength. Thus, the same fault surface can be deemed rougher at long wavelength, or smoother, depending on the chosen metric. This apparent contradiction is a curiosity of the statistics of rough surfaces that have scaling exponents that relate profile length to Rq between 0.5 and 1, as most natural faults do. To add context, the implied roughness scaling for reference synthetic surfaces is determined. These span the natural range of scaling exponents and have moderate to strong point to point amplitude correlation. The potential payoff of expanded descriptions of natural fault roughness and of reference surfaces are improved constraints on physical mechanisms that generate and modify roughness during shear.

Open-File Report↗

USGS2021

Teleseismic waves reveal anisotropic poroelastic response of wastewater disposal reservoir

Connecting earthquake nucleation in basement rock to fluid injection in basal, sedimentary reservoirs, depends heavily on choices related to the poroelastic properties of the fluid-rock system, thermo-chemical effects notwithstanding. Direct constraints on these parameters outside of laboratory settings are rare, and it is commonly assumed that the rock layers are isotropic. With the Arbuckle wastewater disposal reservoir in Osage County, Oklahoma, high-frequency formation pressure changes and collocated broadband ground velocities measured during the passing of large teleseismic waves show a poroelastic response of the reservoir that is both azimuthally variable and anisotropic; this includes evidence of static shifts in pressure that presumably relate to changes in local permeability. The azimuthal dependence in both the static response and shear coupling appears related to tectonic stress and strain indicators such as the orientations of the maximum horizontal stress and faults and fractures. Using dynamic strains from a nearby borehole strainmeter, we show that the ratio of shear to volumetric strain coupling is ∼ 0.41 "> ∼ 0.41 ∼0.41 which implies a mean Skempton's coefficient of A = 0.24 "> A = 0.24 A=0.24 over the plausible range of the undrained Poisson's ratio. Since these observations are made at relatively low confining pressure and differential stress, we suggest that the hydraulically conductive fracture network is a primary control on the coupling between pore pressure diffusion and elastic stresses in response to natural or anthropogenic sources.

Earth and Planetary Physics↗

USGS2021

Evolution of fluid transmissivity and strength recovery of shear fractures under hydrothermal conditions

Geothermal systems rely on the presence of long-lived and high-volume, permeable fracture systems. The creation, reactivation, and sustainability of these systems depend on complex coupling among thermal, hydraulic, mechanical, and chemical (THMC) processes occurring in geothermal reservoirs. In part due to a paucity of experimental data, the evolution of fractures at geothermal conditions in response to THMC processes is poorly understood, particularly during the process of shear. We present preliminary results of triaxial slide-hold-slide experiments, with hold periods ranging in duration from 103 s to 106 s, to constrain rates and mechanisms of healing and sealing. Experiments were conducted on simulated fault gouge composed of Westerly granite and on bare surfaces of Westerly granite. The tests were run at temperatures of 22˚ and 200˚C with confining and average pore pressures of 30 MPa and 10 MPa, respectively. We used an axial displacement rate of 0.1 μm/s during sliding periods. Deionized water was continuously flowed along the simulated fracture so we could determine in-plane transmissivity during the tests. In gouge and bare surface experiments conducted at 200˚C, we observe significant decreases in transmissivity over the course of the experiments. For the hydrothermal gouge experiment we measured an order of magnitude net reduction in transmissivity from 1.73x10-18 to 0.17x10-18 m3, over the course of 220 hours while in the room temperature gouge experiment transmissivity only decreased by 0.35x10-18 m3 over the same amount of time. In the experiments, we observe an up to 16% recovery in transmissivity during sliding periods. At room temperature the friction data showed limited fault re-strengthening with time; healing rates are on the order of 0.1 MPa/decade. A similar healing rate was observed at 200˚C in the gouge but we observe an increase in the healing rate, to 0.75 MPa/decade, for a bare surface experiment at 200˚C. The differences in the healing rate of the gouge and bare surface experiments suggest that the generation of fine particles by grinding down of asperities on the bare surface promote quartz dissolution and reprecipitation at elevated temperatures. Further work is needed to test this possibility and provide better constraints on factors influencing the evolution of fluid transport properties and strength of shear fractures at geothermal conditions.

Conference Paper↗

USGS2021

Strength recovery and sealing under hydrothermal conditions

While there is significant evidence for healing in natural faults, geothermal reservoirs, and lab experiments, the thermal, hydraulic, mechanical, and chemical interactions that influence healing are poorly understood. We present preliminary results of triaxial slide-hold-slide experiments to constrain rates and mechanisms of healing. Experiments were conducted on gouge composed of Westerly granite and on bare surfaces of Westerly granite and Eureka quartzite. Tests were run at 22, 100, and 200˚C. In some experiments, we also determined the in-plane fluid transmissivity. In bare surface experiments we observe that restrengthening depends on both time and temperature. At 200˚C the simulated fractures restrengthen at a rate of ∆ µ /∆log(t hold ) = 0.009/decade while at 22˚C the healing rate is ~ 0.002/decade. In the gouge experiments restrengthening appears to be independent of temperature. This may be related to the heterogenous mineral composition and thickness of the gouge layer which could allow shearing to be accommodated in unhealed zones. In the experiments, an overall reduction in fluid transmissivity is observed but sliding periods are often associated with increases in the fluid transmissivity. The transmissivity reduction tends to be greater at 200˚C relative to room temperature. Our preliminary results suggest that multiple healing mechanisms are operating under hydrothermal conditions.

Conference Paper↗

USGS2020

Direct evidence for fluid pressure, dilatancy, and compaction affecting slip in isolated faults

Earthquake instability occurs as a result of strength loss during sliding on a fault. It has been known for over 50 years that fault compaction or dilatancy may cause significant weakening or strengthening by dramatically changing the fluid pressure trapped in faults. Despite this fundamental importance, we have no real understanding of the exact conditions that lead to compaction or dilation during nucleation or rupture. To date, no direct measurements of pore pressure changes during slip in hydraulically isolated faults have been reported. We show direct examples of fluid pressure variations during nucleation and rupture using a miniature pressure transducer embedded in an experimental fault. We demonstrate that fluids are not only significant in controlling fault behavior, but can provide the dominant mechanism controlling fault stability. The effect of fluid pressure changes can exceed frictional variations predicted by rate- and state-dependent friction laws, exerting fundamental controls on earthquake rupture initiation.

Geophysical Research Letters↗

USGS2020

Near-fault velocity spectra from laboratory failures and their relation to natural ground motion

We compared near-fault velocity spectra recorded during laboratory experiments to that of natural earthquakes. We fractured crystalline rock samples at room temperature and intermediate confining pressure (50 MPa). Subsequent slip events were generated on the fracture surfaces under higher confinement (300 MPa). Velocity spectra from rock fracture resemble the inverse frequency (1/ f ) decay of natural earthquake velocity. This spectrum can be attributed to fault creation via seismic fracturing over a wide range of spatial scales. In contrast, subsequent slips on the rough fracture surfaces are depleted in high frequency energy and falloff approximately as 1/ f 2 . The 1/ f 2 spectrum is more consistent with a slider-block model obeying static-kinetic friction than a natural earthquake. The depleted high frequency content precludes the rough fault experiments from being directly analogous to natural sources. The suppression of high frequencies may have resulted from two possible factors: (1) the presence of a well-developed shear zone and coseismic damping of the fault motion by dissipation within it or, in our favored interpretation, (2) a smaller amount of energy dissipated by shearing relative to the total energy release at elevated confining pressure. In context of the latter explanation, a unifying concept that applies to these experiments, earthquakes, ground motion, and models of complex radiated motion is that high frequency radiated energy is relatively enhanced when total energy release is nearly balanced within the source region by dissipative processes. This near-critical energy release condition can be accessed at low normal stress in laboratory experiments.

Journal of Geophysical Research↗

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↗