Slip-rate-dependent friction as a universal mechanism for slow slip events

A growing body of observations worldwide has documented fault slip transients that radiate little or no seismic energy. The mechanisms that govern these slow slip events (SSEs) and their wide range of depths, slip rates, durations, stress drops and recurrence intervals remain poorly known. Here we show that slow slip can be explained by a transition from rate-weakening frictional sliding at low slip rates towards rate-neutral or rate-strengthening behaviour at higher slip rates, as has been observed experimentally. We use numerical simulations to illustrate that this rate-dependent transition quantitatively explains the experimental data for natural fault rocks representative of materials in the source regions of SSEs. With a standard constant-parameter rate-and-state friction law, SSEs arise only near the threshold for slip instability. The inclusion of velocity-dependent friction parameters substantially broadens the range of conditions for slow slip occurrence, and produces a wide range of event characteristics, which include stress drop, duration and recurrence, as observed in nature. Upscaled numerical simulations that incorporate parameters consistent with laboratory measurements can reproduce geodetic observations of repeating SSEs on tectonic faults. We conclude that slip-rate-dependent friction explains the ubiquitous occurrence of SSEs in a broad range of geological environments

Foreshock properties illuminate nucleation processes of slow and fast laboratory earthquakes

Laboratory experiments demonstrate that prior to fast laboratory earthquakes the fault begins to unlock and creep, causing foreshocks to coalesce in both space and time. This demonstrates that the evolution of foreshocks is closely connected to the fault slip velocity.Understanding the connection between seismic activity and the earthquake nucleation process is a fundamental goal in earthquake seismology with important implications for earthquake early warning systems and forecasting. We use high-resolution acoustic emission (AE) waveform measurements from laboratory stick-slip experiments that span a spectrum of slow to fast slip rates to probe spatiotemporal properties of laboratory foreshocks and nucleation processes. We measure waveform similarity and pairwise differential travel-times (DTT) between AEs throughout the seismic cycle. AEs broadcasted prior to slow labquakes have small DTT and high waveform similarity relative to fast labquakes. We show that during slow stick-slip, the fault never fully locks, and waveform similarity and pairwise differential travel times do not evolve throughout the seismic cycle. In contrast, fast laboratory earthquakes are preceded by a rapid increase in waveform similarity late in the seismic cycle and a reduction in differential travel times, indicating that AEs begin to coalesce as the fault slip velocity increases leading up to failure. These observations point to key differences in the nucleation process of slow and fast labquakes and suggest that the spatiotemporal evolution of laboratory foreshocks is linked to fault slip velocity

Frictional properties of low-angle normal fault gouges and implications for low-angle normal fault slip

The mechanics of slip on low-angle normal faults (LANFs) remain an enduring problem in structural geology and fault mechanics. In most cases, new faults should form rather than having slip occur on LANFs, assuming values of fault friction consistent with Byerlee's Law. We present results of laboratory measurements on the frictional properties of natural clay-rich gouges from low-angle normal faults (LANF) in the American Cordillera, from the Whipple Mts. Detachment, the Panamint range-front detachment, and the Waterman Hills detachment. These clay-rich gouges are dominated by neoformed clay minerals and are an integral part of fault zones in many LANFs, yet their frictional properties under in situ conditions remain relatively unknown. We conducted measurements under saturated and controlled pore pressure conditions at effective normal stresses ranging from 20 to 60 MPa (corresponding to depths of 0.9-2.9 km), on both powdered and intact wafers of fault rock. For the Whipple Mountains detachment, friction coefficient (μ) varies depending on clast content, with values ranging from 0.40 to 0.58 for clast-rich material, and 0.29-0.30 for clay-rich gouge. Samples from the Panamint range-front detachment were clay-rich, and exhibit friction values of 0.28 to 0.38, significantly lower than reported from previous studies on fault gouges tested under room humidity (nominally dry) conditions, including samples from the same exposure. Samples from the Waterman Hills detachment are slightly stronger, with μ ranging from 0.38 to 0.43. The neoformed gouge materials from all three localities exhibits velocity-strengthening frictional behavior under almost all of the experimental conditions we explored, with values of the friction rate parameter (a − b) ranging from −0.001 to +0.025. Clast-rich samples exhibited frictional healing (strength increases with hold time), whereas clay-rich samples do not. Our results indicate that where clay-rich neoformed gouges are present along LANFs, they provide a mechanically viable explanation for slip on faults with dips <20 • , requiring only moderate (P f < σ 3 ) overpressures and/or correcting for ∼5 • of footwall tilting. Furthermore, the low rates of frictional strength recovery and velocity-strengthening frictional behavior we observe provide an explanation for the lack of observed seismicity on these structures. We suggest that LANFs in the upper crust (depth <8 km) slip via a combination of a) reaction-weakening of initially high-angle fault zones by the formation of neoformed clay-rich gouges, and b) regional tectonic accommodation of rotating fault blocks

Determination of stress state in deep subsea formation by combination of hydraulic fracturing in situ test and core analysis: A case study in the IODP Expedition 319

[1] In situ test of hydraulic fracturing (HF) provides the only way to observe in situ stress magnitudes directly. The maximum and minimum horizontal stresses, SHmax and Shmin, are determined from critical borehole pressures, i.e., the reopening pressure Pr and the shut-in pressure Ps, etc, observed during the test. However, there is inevitably a discrepancy between actual and measured values of the critical pressures, and this discrepancy is very significant for Pr. For effective measurement of Pr, it is necessary for the fracturing system to have a sufficiently small compliance. A diagnostic procedure to evaluate whether the compliance of the employed fracturing system is appropriate for SHmax determination from Pr was developed. Furthermore, a new method for stress measurement not restricted by the system compliance and Pr is herein proposed. In this method, the magnitudes and orientations of SHmax and Shmin are determined from (i) the cross-sectional shape of a core sample and (ii) Ps obtained by the HF test performed near the core depth. These ideas were applied for stress measurement in a central region of the Kumano fore-arc basin at a water depth of 2054?m using a 1.6?km riser hole drilled in the Integrated Ocean Drilling Program (IODP) Expedition 319. As a result, the stress decoupling through a boundary at 1285?m below seafloor was detected. The boundary separates new upper layers and old lower ones with an age gap of ~1.8?Ma, which is possibly the accretionary prism. The stress state in the lower layers is consistent with that observed in the outer edge of accretionary prism

Does hydrologic circulation mask frictional heat on faults after large earthquakes?

Knowledge of frictional resistance along faults is important for understanding the mechanics of earthquakes and faulting. The clearest in situ measure of fault friction potentially comes from temperature measurements in boreholes crossing fault zones within a few years of rupture. However, large temperature signals from frictional heating on faults have not been observed. Unambiguously interpreting the coseismic frictional resistance from small thermal perturbations observed in borehole temperature profiles requires assessing the impact of other potentially confounding thermal processes. We address several issues associated with quantifying the temperature signal of frictional heating including transient fluid flow associated with the earthquake, thermal disturbance caused by borehole drilling, and heterogeneous thermal physical rock properties. Transient fluid flow is investigated using a two‐dimensional coupled fluid flow and heat transport model to evaluate the temperature field following an earthquake. Simulations for a range of realistic permeability, frictional heating, and pore pressure scenarios show that high permeabilities (>10−14 m2) are necessary for significant advection within the several years after an earthquake and suggest that transient fluid flow is unlikely to mask frictional heat anomalies. We illustrate how disturbances from circulating fluids during drilling diffuse quickly leaving a robust signature of frictional heating. Finally, we discuss the utility of repeated borehole temperature profiles for discriminating between different interpretations of thermal perturbations. Our results suggest that temperature anomalies from even low friction should be detectable at depths >1 km 1 to 2 years after a large earthquake and that interpretations of low friction from existing data are likely robust.Keywords: Frictional heat on faults, Hydrologic circulationKeywords: Frictional heat on faults, Hydrologic circulatio

Fluid budgets along the northern Hikurangi subduction margin, New Zealand: the effect of a subducting seamount on fluid pressure

We estimate fluid sources around a subducted seamount along the northern Hikurangi subduction margin of New Zealand, using thermomechanical numerical modelling informed by wedge structure and porosities from multichannel seismic data. Calculated fluid sources are input into an independent fluid-flow model to explore the key controls on overpressure generation to depths of 12 km. In the thermomechanical models, sediment transport through and beneath the wedge is calculated assuming a pressure-sensitive frictional rheology. The change in porosity, pressure and temperature with calculated rock advection is used to compute fluid release from compaction and dehydration. Our calculations yield more precise information about source locations in time and space than previous averaged estimates for the Hikurangi margin. The volume of fluid release in the wedge is smaller than previously estimated from margin-averaged calculations (∼14 m3 yr−1 m−1), and is exceeded by fluid release from underlying (subducting) sediment (∼16 m3 yr−1 m−1). Clay dehydration contributes only a small quantity of fluid by volume (∼2 m3 yr−1 m−1 from subducted sediment), but the integrated effect is still significant landward of the seamount. Fluid source terms are used to estimate fluid pressures around a subducting seamount in the fluid-flow models, using subducted sediment permeability derived from porosity, and testing two end-members for décollement permeability. Models in which the décollement acts as a fluid conduit predict only moderate fluid overpressure in the wedge and subducting sediment. However, if the subduction interface becomes impermeable with depth, significant fluid overpressure develops in subducting sediment landward of the seamount. The location of predicted fluid overpressure and associated dehydration reactions is consistent with the idea that short duration, shallow, slow slip events (SSEs) landward of the seamount are caused by anomalous fluid pressures; alternatively, it may result from frictional effects of changing clay content along the subduction interface

Produced Water Management and Beneficial Use

Large quantities of water are associated with the production of coalbed methane (CBM) in the Powder River Basin (PRB) of Wyoming. The chemistry of co-produced water often makes it unsuitable for subsequent uses such as irrigated agriculture. However, co-produced waters have substantial potential for a variety of beneficial uses. Achieving this potential requires the development of appropriate water management strategies. There are several unique characteristics of co-produced water that make development of such management strategies a challenge. The production of CBM water follows an inverse pattern compared to traditional wells. CBM wells need to maintain low reservoir pressures to promote gas production. This need renders the reinjection of co-produced waters counterproductive. The unique water chemistry of co-produced water can reduce soil permeability, making surface disposal difficult. Unlike traditional petroleum operations where co-produced water is an undesirable by-product, co-produced water in the PRB often is potable, making it a highly valued resource in arid western states. This research project developed and evaluated a number of water management options potentially available to CBM operators. These options, which focus on cost-effective and environmentally-sound practices, fall into five topic areas: Minimization of Produced Water, Surface Disposal, Beneficial Use, Disposal by Injection and Water Treatment. The research project was managed by the Colorado Energy Research Institute (CERI) at the Colorado School of Mines (CSM) and involved personnel located at CERI, CSM, Stanford University, Pennsylvania State University, the University of Wyoming, the Argonne National Laboratory, the Gas Technology Institute, the Montana Bureau of Mining and Geology and PVES Inc., a private firm

Operational Review of the First Wireline In Situ Stress Test in Scientific Ocean Drilling

Scientific ocean drilling’s first in situ stress measurement was made at Site C0009A during Integrated Ocean Drilling Program (IODP) Expedition 319 as part of Nankai Trough Seismogenic Zone Experiment (NanTroSEIZE) Stage 2. The Modular Formation Dynamics Tester (MDT, Schlumberger)wireline logging tool was deployed in riser Hole C0009A to measure in situ formation pore pressure, formation permeability (often reported as mobility=permeability/viscosity), and the least principal stress (S3) at several isolated depths (Saffer et al., 2009; Expedition 319 Scientists, 2010). The importance of in situ stress measurements is not only for scientific interests in active tectonic drilling, but also for geomechanical and well bore stability analyses. Certain in situ tools were not previously available for scientific ocean drilling due to the borehole diameter and open hole limits of riserless drilling. The riser-capable drillship, D/V Chikyu,now in service for IODP expeditions, allows all of the techniques available to estimate the magnitudes and orientations of 3-D stresses to be used. These techniques include downhole density logging for vertical stress, breakout and caliper log analyses for maximum horizontal stress, core-based anelastic strain recovery (ASR, used in the NanTroSEIZE expeditions in 2007–2008), and leak-off test (Lin et al., 2008) and minifrac/hydraulic fracturing (NanTroSEIZE Expedition319 in 2009). In this report, the whole operational planning process related to in situ measurements is reviewed, and lessons learned from Expedition 319 are summarized for efficient planning and testing in the future

Asymmetric brittle deformation at the Pāpaku Fault, Hikurangi Subduction Margin, NZ, IODP Expedition 375

Quantifying fault damage zones provides a window into stress distribution and rheology around faults. International Ocean Discovery Program (IODP) Expeditions 372/375 drilled an active thrust splay fault within the Hikurangi subduction margin. The fault, which is hosted in Pleistocene clastic sediments, is surrounded by brittle fractures and faults as well as ductile deformation features. We find that fracture density in the damage zone enveloping the fault is asymmetric, with the hanging wall showing greater overall fracture density and at greater distances from the fault than the footwall. Furthermore, the peak in fracture density occurs within an area of mesoscale folding and localized slip in the hanging wall rather than adjacent to the main fault zone. We attribute the asymmetry in damage to disparate deformation histories between the hanging wall and footwall, greater ductile deformation within the footwall, and/or dynamic stress asymmetry around a propagating rupture. Damage asymmetry is common at shallow depths in subduction zones and influences the mechanical and hydrological properties of the fault, such as channelized fluid flow and fault stability. Finally, we demonstrate that subduction zone faults show similar damage-displacement scaling as continental faults