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.


Title: Homogenization of fault frictional properties
SUMMARY The frictional properties of large faults are expected to vary in space. However, fault models often assume that properties are homogeneous, or nearly so. We investigate the conditions under which the details of variations may be neglected and properties homogenized. We do so by examining the behaviour of nonlinear solutions for unstably accelerating fault slip under frictional heterogeneity. We consider a rate- and state-dependent fault friction in which the characteristic wavelength for the property variations is a problem parameter. We find that homogenization is permissible only when that wavelength shows scale separation from an elasto-frictional length scale. However, fault models also often include property transitions that occur over distances comparable to the elasto-frictional length. We show that under such comparable variations, the dynamics of earthquake-nucleating instabilities is controlled by the properties’ spatial distribution.  more » « less
Award ID(s):
1653382
PAR ID:
10163009
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Geophysical Journal International
Volume:
219
Issue:
2
ISSN:
0956-540X
Page Range / eLocation ID:
1203 to 1211
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT Large earthquakes on strike-slip faults often rupture multiple fault segments by jumping over stepovers. Previous studies, based on field observations or numerical modeling with a homogeneous initial stress field, have suggested that stepovers more than ∼5  km wide would stop the propagation of rupture, but many exceptions have been observed in recent years. Here, we integrate a dynamic rupture model with a long-term fault stress model to explore the effects of background stress perturbation on rupture propagation across stepovers along strike-slip faults. Our long-term fault models simulate steady-state stress perturbation around stepovers. Considering such stress perturbation in dynamic rupture models leads to prediction of larger distance a dynamic rupture can jump over stepovers: over 15 km for a releasing stepover or 7 km for a restraining stepover, comparing with the 5 km limit in models with the same fault geometry and frictional property but assuming a homogeneous initial stress. The effect of steady-state stress perturbations is stronger in an overlapping stepover than in an underlapping stepover. The maximum jumping distance can reach 20 km in an overlapping releasing stepover with low-static frictional coefficients. These results are useful for estimating the maximum length of potential fault ruptures and assessing seismic hazard. 
    more » « less
  2. There is scientific and industrial interest in understanding how geologic faults respond to transient sources of fluid. Natural and artificial sources can elevate pore fluid pressure on the fault frictional interface, which may induce slip. We consider a simple boundary value problem to provide an elementary model of the physical process and to provide a benchmark for numerical solution procedures. We examine the slip of a fault that is an interface of two elastic half-spaces. Injection is modelled as a line source at constant pressure and fluid pressure is assumed to diffuse along the interface. The resulting problem is an integro-differential equation governing fault slip, which has a single dimensionless parameter. The expansion of slip is self-similar and the rupture front propagates at a factor $$\lambda$$ of the diffusive length scale $$\sqrt {\alpha t}$$ . We identify two asymptotic regimes corresponding to $$\lambda$$ being small or large and perform a perturbation expansion in each limit. For large $$\lambda$$ , in the regime of a so-called critically stressed fault, a boundary layer emerges on the diffusive length scale, which lags far behind the rupture front. We demonstrate higher-order matched asymptotics for the integro-differential equation, and in doing so, we derive a multipole expansion to capture successive orders of influence on the outer problem for fault slip for a driving force that is small relative to the crack dimensions. Asymptotic expansions are compared with accurate numerical solutions to the full problem, which are tabulated to high precision. 
    more » « less
  3. SUMMARY This paper examines the linear stability of sliding on faults embedded in a 2-D elastic medium that obey rate and state friction and have a finite length and/or are near a traction-free surface. Results are obtained using a numerical technique that allows for analysis of systems with geometrical complexity and heterogeneous material properties; however only systems with homogeneous frictional and material properties are examined. Some analytical results are also obtained for the special case of a fault that is parallel to a traction-free surface. For velocity-weakening faults with finite length, there is a critical fault length $$L^{*}$$ for unstable sliding that is analogous to the critical wavelength $$h^{*}$$ that is usually derived from infinite fault systems. Faults longer than $$L^{*}$$ are linearly unstable to perturbations of any length. On vertical strike-slip faults or faults in a full-space $$L^{*} \approx h^{*}/e$$, where e is Euler’s number. For dip-slip faults near a traction-free surface $$L^{*} \le h^{*}/e$$ and is a function of dip angle $$\beta$$, burial depth d of the fault’s up-dip edge and friction coefficient. In particular, $$L^{*}$$ is at least an order of magnitude smaller than $$h^{*}$$ on shallowly dipping ($$\beta < 10^\circ$$) faults that intersect the traction-free surface. Additionally, $$L^{*} \approx h^{*}/e$$ on dip-slip faults with burial depths $$d \ge h^{*}$$. For sliding systems that can be treated as a thin layer, such as landslides, glaciers or ice streams, $$L^{*} = h^{*}/2$$. Finally, conditions are established for unstable sliding on infinitely-long, velocity-strengthening faults that are parallel to a traction-free surface. 
    more » « less
  4. Injection-induced seismicity and aseismic slip often involve the reactivation of long-dormant faults, which may have extremely low permeability prior to slip. In contrast, most previous models of fluid-driven aseismic slip have assumed linear pressure diffusion in a fault zone of constant permeability and porosity. Slip occurs within a frictional shear crack whose edge can either lag or lead pressure diffusion, depending on the dimensionless stress-injection parameter that quantifies the prestress and injection conditions. Here, we extend this foundational work by accounting for permeability enhancement and dilatancy, assumed to occur instantaneously upon the onset of slip. The fault zone ahead of the crack is assumed to be impermeable, so fluid flow and pressure diffusion are confined to the interior, slipped part of the crack. The confinement of flow increases the pressurization rate and reduction of fault strength, facilitating crack growth even for severely understressed faults. Suctions from dilatancy slow crack growth, preventing propagation beyond the hydraulic diffusion length. Our new two-dimensional and three-dimensional solutions can facilitate the interpretation of induced seismicity data sets. They are especially relevant for faults in initially low permeability formations, such as shale layers serving as caprock seals for geologic carbon storage, or for hydraulic stimulation of geothermal reservoirs. This article is part of the theme issue ‘Induced seismicity in coupled subsurface systems’. 
    more » « less
  5. Abstract Frictional sliding along grain boundaries in brittle shear zones can result in the fragmentation of individual grains, which ultimately can impact slip dynamics. During deformation at small scales, stick–slip motion can occur between grains when existing force chains break due to grain rearrangement or failure, resulting in frictional sliding of granular material. The rearrangement of the grains leads to dilation of the granular package, reducing the shear stress and subsequently leading to slip. Here, we conduct physical experiments employing HydroOrbs, an elasto-plastic material, to investigate grain comminution in granular media under simple shear conditions. Our findings demonstrate that the degree of grain comminution is dependent on both the normal force and the size of the grains. Using the experimental setup, we benchmark Discrete Element Method (DEM) numerical models, which are capable of simulating the movement, rotation, and fracturing of elasto-plastic grains subjected to simple shear. The DEM models successfully replicate both grain comminution patterns and horizontal force fluctuations observed in our physical experiments. They show that increasing normal forces correlate with higher horizontal forces and more fractured grains. The ability of our DEM models to accurately reproduce experimental results opens up new avenues for investigating various parameter spaces that may not be accessible through traditional laboratory experiments, for example, in assessing how internal friction or cohesion affect deformation in granular systems. 
    more » « less