skip to main content


Title: Porosity Evolution in Rate and State Friction
Abstract

This letter compares the predictions of two expressions proposed for the porosity evolution in the context of rate and state friction. One (Segall & Rice, 1995,https://doi.org/10.1029/95jb02403) depends only on the sliding velocity; the other (Sleep, 1995,https://doi.org/10.1029/94jb03340) depends only on the state variable. Simulations of both are similar for velocity stepping and slide‐hold‐slide experiments. They differ significantly for normal effective stress jumps at constant sliding velocity. Segall and Rice (1995,https://doi.org/10.1029/95jb02403) predicts no change in the porosity; Sleep (1995,https://doi.org/10.1029/94jb03340) does. Simulation with a spring‐block model indicates that the magnitude of rapid slip events is essentially the same for the two formulations. Variations of porosity and induced pore pressure near rapid slip events are similar and consistent with experimental observations. Predicted porosity variations during slow slip intervals and the time at which rapid slip events occur are significantly different. The simulation indicates that changes in friction stress due to pore pressure changes exceed those due to rate and state effects.

 
more » « less
Award ID(s):
2120374
NSF-PAR ID:
10380746
Author(s) / Creator(s):
 
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
49
Issue:
22
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Rate and state frictional parameters are typically determined using two types of experimental protocols: velocity steps and slide‐hold‐slide events. Here we take a new approach by examining the frictional response to controlled, harmonic oscillations in load point velocity. We present a Matlab graphical user interface software package, called RSFitOSC, that allows users to easily determine frictional parameters by fitting oscillation events using the rate and state friction equations. We apply our new methods to a set of ice‐rock friction experiments conducted over a temperature range of −16.4°C to −2°C, and described in a companion paper: McCarthy et al. (2021,https://doi.org/10.1002/essoar.10509831.110.1002/essoar.10509831.1). Values of the frictional stability parameter (ab) determined from oscillations reveal dominantly velocity‐weakening behavior across the entire range of experimental conditions. However, values of (a–b) determined from velocity steps in the same experiments yield velocity‐strengthening behavior. We also show that the elastic stiffness of the ice‐rock system depends on the temperature, and is unlikely to be explained by changes in the elastic properties of ice. Load point velocity oscillations induce oscillations in applied shear stress. Many natural fault systems exhibit slip behaviors that depend on harmonic oscillations in applied tidal stresses. Our new method provides a way to study how frictional properties directly depend on parameters relevant to tidal forcing, and how oscillatory loading must be considered when extracting friction parameters.

     
    more » « less
  2. Abstract

    Earthquake focal mechanisms, determined with P‐wave polarities and S/P amplitude ratios, are primary data for analyzing fault zone geometry, sense of slip, and the crustal stress field. Solving for the focal mechanisms of small earthquakes is often challenging because phase arrivals and first‐motion polarities are hard to be separated from noise. To overcome this challenge, we implement convolutional‐neural‐network algorithms (Ross, Meier, & Hauksson, 2018, Ross, Meier, Hauksson, & Heaton, 2018,https://doi.org/10.1029/2017jb015251,https://doi.org/10.1785/0120180080) to detect additional phases and polarities. Using both existing and these new data, we build a high‐quality focal mechanism catalog of 297,478 events that occurred from 1981 to 2021 in southern California with the HASH method of Hardebeck and Shearer (2002),https://doi.org/10.1785/0120010200, Hardebeck and Shearer (2003),https://doi.org/10.1785/0120020236. The new focal mechanism catalog is overall consistent with the standard catalog (Yang et al., 2012,https://doi.org/10.1785/0120110311) but includes 40% more focal mechanisms, and is more consistent with moment tensor solutions derived using waveform‐fitting methods. We apply the new catalog to identify changes in focal mechanism properties caused by the occurrences of large mainshocks such as the 2010Mw7.2 El Mayor‐Cucapah and 2019Mw7.1 Ridgecrest earthquakes. Such changes may be associated with co‐seismic stress drops, post‐seismic deformation processes, and static stress changes on a regional scale. The new high‐resolution catalog will contribute to improved understanding of the crustal stress field, earthquake triggering mechanisms, fault zone geometry, and sense of slip on the faults in southern California.

     
    more » « less
  3. Abstract

    Wave breaking induced bubbles contribute a significant part of air‐sea gas fluxes. Recent modeling of the sea state dependent CO2flux found that bubbles contribute up to ∼40% of the total CO2air‐sea fluxes (Reichl & Deike, 2020,https://doi.org/10.1029/2020gl087267). In this study, we implement the sea state dependent bubble gas transfer formulation of Deike and Melville (2018,https://doi.org/10.1029/2018gl078758) into a spectral wave model (WAVEWATCH III) incorporating the spectral modeling of the wave breaking distribution from Romero (2019,https://doi.org/10.1029/2019gl083408). We evaluate the accuracy of the sea state dependent gas transfer parameterization against available measurements of CO2gas transfer velocity from 9 data sets (11 research cruises, see Yang et al. (2022,https://doi.org/10.3389/fmars.2022.826421)). The sea state dependent parameterization for CO2gas transfer velocity is consistent with observations, while the traditional wind‐only parameterization used in most global models slightly underestimates the observations of gas transfer velocity. We produce a climatology of the sea state dependent gas transfer velocity using reanalysis wind and wave data spanning 1980–2017. The climatology shows that the enhanced gas transfer velocity occurs frequently in regions with developed sea states (with strong wave breaking and high significant wave height). The present study provides a general sea state dependent parameterization for gas transfer, which can be implemented in global coupled models.

     
    more » « less
  4. Abstract

    We report on the mountain top observation of three terrestrial gamma‐ray flashes (TGFs) that occurred during the summer storm season of 2021. To our knowledge, these are the first TGFs observed in a mountaintop environment and the first published European TGFs observed from the ground. A gamma‐ray sensitive detector was located at the base of the Säntis Tower in Switzerland and observed three unique TGF events with coincident radio sferic data characteristic of TGFs seen from space. We will show an example of a “slow pulse” radio signature (Cummer et al., 2011,https://doi.org/10.1029/2011GL048099; Lu et al., 2011,https://doi.org/10.1029/2010JA016141; Pu et al., 2019,https://doi.org/10.1029/2019GL082743; Pu et al., 2020,https://doi.org/10.1029/2020GL089427), a −EIP (Lyu et al., 2016,https://doi.org/10.1002/2016GL070154; Lyu et al., 2021,https://doi.org/10.1029/2021GL093627; Wada et al., 2020,https://doi.org/10.1029/2019JD031730), and a double peak TGF associated with an extraordinarily powerful and complicated positive‐polarity sferic, where each TGF peak is possibly preceded by a short burst of stepped leader emission.

     
    more » « less
  5. Abstract

    The mechanical, physical, and frictional properties of incoming materials play an important role in subduction zone structure and slip behavior because these properties influence the strength of the accretionary wedge and megathrust plate boundary faults. Incoming sediment sections often show an increase in compressional wave speed (Vp) and a decrease in porosity with depth due to consolidation. These relations allow seismic‐velocity models to be used to elucidate properties and conditions at depth. However, variations in these properties are controlled by lithology and composition as well as cementation and diagenesis. We present an analysis of shipboard measurements of Vpand porosity on incoming sediment cores from International Ocean Discovery Program (IODP) expeditions at the Hikurangi Margin, Nankai Trough, Aleutian Trench, Middle America Trench, and Sunda Trench. Porosity for these samples ranges from 5% to 85% and Vpranges from 1.5 to 6 km/s. Vp‐porosity relations developed by Erikson & Jarrad (1998),https://doi.org/10.1029/98JB02128and Hoffman & Tobin (2004)https://10.2973/odp.proc.sr.190196.355.2004, with a critical porosity of ∼30%, can represent carbonate‐poor (<50 wt% CaCO3), mainly hemipelagic, incoming sediment regardless of the margin. But these relations tend to underestimate porosity in incoming sediments with carbonate content greater than 50 wt%, which appear to have a critical porosity of between 45% and 50%. This discrepancy will lead to inaccuracy in estimates of fluid budget and overpressure in subduction zones. The velocity‐porosity relation in carbonate sediments is non‐unique due to the complexity that results from the greater susceptibility of carbonate rocks to diagenetic processes.

     
    more » « less