skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.
Attention:The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 7:00 AM ET to 7:30 AM ET on Friday, April 24 due to maintenance. We apologize for the inconvenience.


Title: Homogeneous and Membrane Pore Fluid Diffusion in Spring Block Simulations of Fault Slip With Rate and State Friction
Abstract Spring block, and sometimes continuum, models of the effects of the coupling of fluid flow and frictional slip often employ the membrane approximation. This approximation assumes that the fluid flux between the slipping layer and a remote pore fluid reservoir is proportional to the difference between the value of pore pressure in the reservoir and in the slipping zone. In contrast, Darcy's law states that the fluid flux is proportional to the gradient of the pore pressure. We analyze and compare these two formulations by asymptotic analysis of the fluid flow equations and numerical simulations of a spring–block model using rate and state friction. This analysis shows that membrane diffusion agrees with homogeneous diffusion in the limit of undrained conditions and for nearly drained conditions. For homogeneous diffusion and essentially undrained conditions, both the asymptotic analysis and numerical simulations indicate the formation of a boundary layer near the shear zone where gradients of pore pressure are large. Outside this layer the pore pressure rapidly approaches the drained solution. A linearized stability analysis derives the dependence of the non‐dimensional critical stiffness () on fluid diffusivity, dilatancy factor (), and shear zone thickness. The homogeneous and membrane diffusion models exhibit nearly identical in the limits of drained and undrained conditions, but differ between the two limits. In this intermediate range, numerical simulations show that the two models produce similar slip behavior for but significant differences in slip velocity and recurrence patterns for .  more » « less
Award ID(s):
2120374
PAR ID:
10673624
Author(s) / Creator(s):
 ;  
Publisher / Repository:
American Geophysical Union
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
130
Issue:
12
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Marine structures placed in the shallower seabed can experience pore water drainages with more complexity than those in onshore environments, particularly in coarse-grained soils where drainage is neither purely “drained” nor “undrained,” but Partially Drained (PD). However, current laboratory approaches for characterizing soil behavior are limited to modeling drainage conditions as fully drained or undrained. This paper presents results from a series of confined monotonic saturated simple shear tests under various drainage conditions on reconstituted medium dense to dense Monterey sand specimens to fill this knowledge gap. Although others have performed limited PD element-level tests under triaxial conditions, no documentation exists for tests using a simple monotonic shear configuration. To achieve PD, a special filter was fabricated and connected between the bottom of the specimen and the backpressure controller. The hydraulic filter comprises a series of needle valves to provide various hydraulic impedances. All simple shear tests in this paper were backpressure-saturated. Two different degrees of PD were considered and compared with fully drained and undrained conditions. Results show that the excess pore water pressure generation and measured volumetric changes in the PD tests are bounded between those measured from fully drained and undrained, proving the PD filter provided the hydraulic resistance to achieve PD condition. 
    more » « less
  2. ABSTRACT This study investigates linkages between volume change, pore fluid drainage, shear wave velocity, and temperature of soft clays using a thermal triaxial cell equipped with bender elements, a measurement approach that has not been explored widely in past thermo-mechanical studies. Two kaolinite specimens were consolidated mechanically to a normally consolidated state and then subjected to drained and undrained heating-cooling cycles, respectively. After cooling, the specimens were subjected to further mechanical consolidation to evaluate changes in apparent preconsolidation stress. Both specimens showed net contractive thermal strains after a heating-cooling cycle and overconsolidated behavior during mechanical compression immediately after cooling. The shear wave velocity increased during drained heating, but negligible changes were observed during drained cooling, indicating permanent hardening because of thermal consolidation during the heating-cooling cycle. The shear wave velocity decreased during undrained heating because of a reduction in effective stress associated with thermal pressurization of the pore fluid but subsequently increased when drainage was permitted at elevated temperature. The shear wave velocity increased slightly during undrained cooling but decreased when drainage was permitted at room temperature. Net increases in small-strain shear modulus of 17 and 11 % after heating-cooling cycles under drained and undrained (with drainage after reaching stable temperatures) conditions, respectively, provide further evidence to the potential of thermal soil improvement of normally consolidated clays. Transient changes in shear modulus also highlight the importance of considering drainage conditions and corresponding changes in effective stress state during heating-cooling cycles. 
    more » « less
  3. Key Points High pore fluid pressure stabilizes fault propagation in porous sandstone deformed under drained conditions Slow faulting was associated with pervasive microcracking and diffuse shear bands only in samples deformed sufficiently slow Pervasive subcritical cracking enables slow faulting at high pore fluid pressure under drained conditions at the sample scale 
    more » « less
  4. ABSTRACT Frictional heating during earthquake rupture raises the fault-zone fluid pressure, which affects dynamic rupture and seismic radiation. Here, we investigate two key parameters governing thermal pressurization of pore fluids – hydraulic diffusivity and shear-zone half-width – and their effects on earthquake rupture dynamics, kinematic source properties, and ground motions. We conduct 3D strike-slip dynamic rupture simulations assuming a rate-and-state dependent friction law with strong velocity weakening coupled to thermal-pressurization of pore fluids. Dynamic rupture evolution and ground shaking are densely evaluated across the fault and Earth’s surface to analyze the variations of rupture parameters (slip, peak slip rate, rupture speed, and rise time), correlations among rupture parameters, and variability of peak ground velocity. Our simulations reveal how variations in thermal-pressurization affect earthquake rupture properties. We find that the mean slip and rise time decrease with increasing hydraulic diffusivity, whereas mean rupture speed and peak slip-rate remain almost constant. Mean slip, peak slip-rate, and rupture speed decrease with increasing shear-zone half-width, whereas mean rise time increases. Shear-zone half-width distinctly affects the correlation between rupture parameters, especially for parameter pairs (slip, rupture speed), (peak slip-rate, rupture speed), and (rupture speed, rise time). Hydraulic diffusivity has negligible effects on these correlations. Variations in shear-zone half-width primarily impact rupture speed, which then may affect other rupture parameters. We find a negative correlation between slip and peak slip-rate, unlike simpler dynamic rupture models. Mean peak ground velocities decrease faster with increasing shear-zone half-width than with increasing hydraulic diffusivity, whereas ground-motion variability is similarly affected by both the parameters. Our results show that shear-zone half-width affects rupture dynamics, kinematic rupture properties, and ground shaking more strongly than hydraulic diffusivity. We interpret the importance of shear-zone half-width based on the characteristic time of diffusion. Our findings may inform pseudodynamic rupture generators and guide future studies on how to account for thermal-pressurization effects. 
    more » « less
  5. null (Ed.)
    Abstract Fault-zone fluids control effective normal stress and fault strength. While most earthquake models assume a fixed pore fluid pressure distribution, geologists have documented fault valving behavior, that is, cyclic changes in pressure and unsteady fluid migration along faults. Here we quantify fault valving through 2-D antiplane shear simulations of earthquake sequences on a strike-slip fault with rate-and-state friction, upward Darcy flow along a permeable fault zone, and permeability evolution. Fluid overpressure develops during the interseismic period, when healing/sealing reduces fault permeability, and is released after earthquakes enhance permeability. Coupling between fluid flow, permeability and pressure evolution, and slip produces fluid-driven aseismic slip near the base of the seismogenic zone and earthquake swarms within the seismogenic zone, as ascending fluids pressurize and weaken the fault. This model might explain observations of late interseismic fault unlocking, slow slip and creep transients, swarm seismicity, and rapid pressure/stress transmission in induced seismicity sequences. 
    more » « less