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: Numerical exploration of factors that control stress amplification in Earth's lithosphere
SUMMARY The mechanical heterogeneity of Earth's lithosphere leads to significant amplification of stresses across spatial scales ranging from mineral grains to tectonic plates. These stress amplifications play a key role in mechanical and chemical processes within the rock that affect bulk rock strength. Identifying the most effective causes of stress amplification is critical for understanding processes such as strain localization and fluid transport at scales ranging from microshear zones to tectonic plate boundaries. However, studies quantifying and predicting stress heterogeneities and amplifications are limited. We used numerical modelling of two-phase isotropic viscous systems to explore the factors influencing and controlling stress amplification and the potential magnitude of stress amplification in viscous regimes. We found the most geologically relevant amplification factors to be weak-phase spacing, rheological contrast and loading type. Our results indicate that stress amplification can reach a factor of ∼9 under specific conditions, but most of our experiments suggest amplifications at or below a factor of 2. Pressure differences across the model domains generally do not exceed ∼55 MPa, but some are as high as ∼110 MPa. The stress and pressure amplifications resulting from our analyses are large enough to drive a variety of geologically important processes such as failure and strain localization, as well as transient permeability and fluid migration.  more » « less
Award ID(s):
2150831
PAR ID:
10472833
Author(s) / Creator(s):
; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Geophysical Journal International
Volume:
235
Issue:
3
ISSN:
0956-540X
Format(s):
Medium: X Size: p. 2927-2947
Size(s):
p. 2927-2947
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Critical processes including seismic faulting, reservoir compartmentalization, and borehole failure involve high‐pressure mechanical behavior and strain localization of sedimentary rocks such as sandstone. Sand is often used as a model material to study the mechanical behavior of poorly lithified sandstone. Recent studies exploring the multi‐scale mechanics of sand have characterized the brittle, low‐pressure regime of behavior; however, limited work has provided insights into the ductile, high‐pressure regime of behavior viain‐situmeasurements. Critical features of the ductile regime, including grain breakage, grain micromechanics, and volumetric strain behavior therefore remain under‐explored. Here, we use a new high‐pressure triaxial apparatus within‐situx‐ray tomography to provide new insights into deformation banding, grain breakage, and grain micromechanics in Ottawa sand subjected to triaxial compression under confining pressures between 10 and 45 MPa. We observed strain‐hardening at pressures above 15 MPa and strain‐neutral responses at pressures below 15 MPa. Compacting shear bands and grain breakage were observed at all pressures with no significant variation due to grain size, except for minor increases in breakage in less‐rounded sands. Grain breakage emerged at stress levels lower than the assumed yield threshold and more intense breakage was associated with thinner deformation bands. Contact sliding at inter‐grain contacts demonstrated a bifurcation into a bimodal distribution, with intense sliding within deformation bands and reduced but non‐negligible sliding outside of deformation bands, suggesting that off‐band zones remain mechanically active during strain hardening. 
    more » « less
  2. Abstract Carbonate sediments play a prominent role on the global geological stage as they store more than $$60\%$$ 60 % of world’s oil and $$40\%$$ 40 % of world’s gas reserves. Prediction of the deformation and failure of porous carbonates is, therefore, essential to minimise reservoir compaction, fault reactivation, or wellbore instability. This relies on our understanding of the mechanisms underlying the observed inelastic response to fluid injection or deviatoric stress perturbations. Understanding the impact of deformation/failure on the hydraulic properties of the rock is also essential as injection/production rates will be affected. In this work, we present new experimental results from triaxial deformation experiments carried out to elucidate the behaviour of a porous limestone reservoir analogue (Savonnières limestone). Drained triaxial and isotropic compression tests were conducted at five different confining pressures in dry and water-saturated conditions. Stress–strain data and X-ray tomography images of the rock indicate two distinct types of deformation and failure regimes: at low confinement (10 MPa) brittle failure in the form of dilatant shear banding was dominant; whereas at higher confinement compaction bands orthogonal to the maximum principal stress formed. In addition to the pore pressure effect, the presence of water in the pore space significantly weakened the rock, thereby shrinking the yield envelope compared to the dry conditions, and shifted the brittle–ductile transition to lower effective confining pressures (from 35 MPa to 29 MPa). Finally, permeability measurements during deformation show a reduction of an order of magnitude in the ductile regime due to the formation of the compaction bands. These results highlight the importance of considering the role of the saturating fluid in the brittle–ductile response of porous rocks and elucidate some of the microstructural processes taking place during this transition. 
    more » « less
  3. Abstract Dilatant hardening is one proposed mechanism that causes slow earthquakes along faults. Previous experiments and models show that dilatant hardening can stabilize fault rupture and slip in several lithologies. However, few studies have systematically measured the mechanical behavior across the transition from dynamic to slow rupture or considered how the associated damage varies. To constrain the processes and scales of dilatant hardening, we conducted triaxial compression experiments on cores of Crab Orchard sandstone and structural analyses using micro‐computed tomography imaging and petrographic analysis. Experiments were conducted at an effective confining pressure of ∼10 MPa, while varying confining pressure (10–130 MPa) and pore fluid pressure (1–120 MPa). Above 15 MPa pore fluid pressure, dilatant hardening slows the rate of fault rupture and slip and deformation becomes more distributed amongst multiple faults as microfracturing increases. The resulting increase in fracture energy has the potential to control fault slip behavior. 
    more » « less
  4. SUMMARY The occurrence of plate tectonics on Earth is rooted in the physics of lithospheric ductile weakening and shear-localization. The pervasiveness of mylonites at lithospheric shear zones is a key piece of evidence that localization correlates with reduction in mineral grain size. Most lithospheric mylonites are polymineralic and the interaction between mineral phases, such as olivine and pyroxene, especially through Zener pinning, impedes normal grain growth while possibly enhancing grain damage, both of which facilitate grain size reduction and weakening, as evident in lab experiments and field observations. The efficacy of pinning, however, relies on the mineral phases being mixed and dispersed at the grain scale, where well-mixed states lead to greater mylonitization. To model grain mixing between different phases at the continuum scale, we previously developed a theory treating grain-scale processes as diffusion between phases, but driven by imposed compressive stresses acting on the boundary between phases. Here we present a new model for shearing rock that combines our theory for diffusive grain mixing, 2-D non-Newtonian flow and two-phase grain damage. The model geometry is designed specifically for comparison to torsional shear-deformation experiments. Deformation is either forced by constant velocity or constant stress boundary conditions. As the layer is deformed, mixing zones between different mineralogical units undergo enhanced grain size reduction and weakening, especially at high strains. For constant velocity boundary experiments, stress drops towards an initial piezometric plateau by a strain of around 4; this is also typical of monophase experiments for which this initial plateau is the final steady state stress. However, polyphase experiments can undergo a second large stress drop at strains of 10–20, and which is associated with enhanced phase mixing and resultant grain size reduction and weakening. Model calculations for polyphase media with grain mixing and damage capture the experimental behaviour when damage to the interface between phases is moderately slower or less efficient than damage to the grain boundaries. Other factors such as distribution and bulk fraction of the secondary phase, as well as grain-mixing diffusivity also influence the timing of the second stress drop. For constant stress boundary conditions, the strain rate increases during weakening and localization. For a monophase medium, there is theoretically one increase in strain rate to a piezometric steady state. But for the polyphase model, the strain rate undergoes a second abrupt increase, the timing for which is again controlled by interface damage and grain mixing. The evolution of heterogeneity through mixing and deformation, and that of grain size distributions also compare well to experimental observations. In total, the comparison of theory to deformation experiments provides a framework for guiding future experiments, scaling microstructural physics to geodynamic applications and demonstrates the importance of grain mixing and damage for the formation of plate tectonic boundaries. 
    more » « less
  5. Abstract Elevated pore fluid pressure is proposed to contribute to slow earthquakes along shallow subduction plate boundaries. However, the processes that create high fluid pressure, disequilibrium compaction and dehydration reactions, lead to different effective stress paths in fault rocks. These paths are predicted by granular mechanics frameworks to lead to different strengths and deformation modes, yet granular mechanics do not predict their effects on fault stability. To evaluate the role of fluid overpressure on shallow megathrust strength and slip behavior, we conducted triaxial shear experiments on chlorite and celadonite rich gouge layers. Experiments were conducted at constant temperature (130 and 100°C), confining pressure (130 and 140 MPa), and pore fluid pressures (between 10 and 120 MPa). Fluid overpressure due to disequilibrium compaction was simulated by increasing confining and pore fluid pressure synchronously without exceeding the target effective pressure, whereas overpressure due to dehydration reactions was simulated by first loading the sample to a target isotropic effective pressure and then increasing pore fluid pressure to reduce the effective pressure. We find that the effects of fluid pressure and stress path on the mechanical behavior of the chlorite and celadonite gouges can generally be described using the critical state soil mechanics (CSSM) framework. However, path effects are more pronounced and persist to greater displacements in chlorite because its microstructure is more influenced by stress path. Due to its effects on microstructure, the stress path also imparts greater control on the rate‐dependence of chlorite strength, which is not predicted by CSSM. 
    more » « less