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: A coupled model for phase mixing, grain damage and shear localization in the lithosphere: comparison to lab experiments
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
Award ID(s):
1853155
PAR ID:
10464846
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Geophysical Journal International
Volume:
232
Issue:
3
ISSN:
0956-540X
Page Range / eLocation ID:
2205 to 2230
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Oceanic transform faults represent abundant yet relatively unexplored components of the hydrologic cycle in the mantle lithosphere. Current models limit fluid circulation to 600 °C, the thermal limit of earthquakes recorded by teleseismic surveys. However, recent ocean‐bottom seismic surveys have located earthquakes at depths corresponding to >1000 °C in modeled thermal structure. To constrain the depth extent of brittle deformation and fluid infiltration, we analyzed peridotite mylonites dredged from the Shaka Transform Fault, Southwest Indian Ridge. Samples range from high strain mylonites that preserve ductile microstructures to lower strain mylonites that are fractured and overprinted by hydrothermal alteration. Microstructural analysis of the high strain samples reveals brittle deformation of pyroxene concomitant with ductile deformation of olivine and growth of amphibole. Porphyroclasts preserve healed fractures filled with fluid inclusions, implying repeated episodes of fracture, fluid infiltration, and healing. The association of hydration features with brittle structures points to seawater, rather than melt, as the fluid source. Textural analysis indicates that strain localization was initiated by grain boundary pinning and that olivine grain size was reduced to ~1 μm in the presence of amphibole. Comparing the amphibole stability field to thermometry estimates for the limit of recrystallization suggests that fluid flow extended to ~650–850 °C. Our results indicate that the hydrologic cycle extends past the brittle‐ductile transition and promotes strain localization via hydrolytic weakening and hydration reactions. We propose that seawater infiltration on oceanic transform faults is driven by the seismic cycle and represents a first order control on the rheology of the oceanic lithosphere. 
    more » « less
  2. Abstract Quartz deformation fabrics reflect stress and strain conditions in mylonites, and their interpretation has become a mainstay of kinematic and structural analysis. Quantification of grain size and shape and interpretation of textures reflecting deformation mechanisms can provide estimates of flow stress, strain rate, kinematic vorticity, and deformation temperatures. Empirical calibration and determination of quartz flow laws is based on laboratory experiments of pure samples; however, pure quartzite mylonites are relatively uncommon. In particular, phyllosilicates may localize and partition strain that can inhibit or enhance different deformation mechanisms. Experimental results demonstrate that even minor phyllosilicate content (<15 vol%) can dramatically alter the strain behavior of quartz; however, few field studies have demonstrated these effects in a natural setting. To investigate the role of phyllosilicates on quartz strain fabrics, we quantify phyllosilicate content and distribution in quartzite mylonites from the Miocene Raft River detachment shear zone (NW Utah, USA). We use microstructural analysis and electron backscatter diffraction to quantify quartz deformation fabrics and muscovite spatial distribution, and X-ray computed tomography to quantify muscovite content in samples with varying amounts of muscovite collected across the detachment shear zone. Phyllosilicate content has a direct control on quartz deformation mechanisms, and application of piezometers and flow laws based on quartz deformation fabrics yield strain rates and flow stresses that vary by up to two orders of magnitude across our samples. These findings have important implications for the application of flow laws in quartzite mylonites and strain localization mechanisms in mid-crustal shear zones. 
    more » « less
  3. How subduction—the sinking of cold lithospheric plates into the mantle—is initiated is one of the key mysteries in understanding why Earth has plate tectonics. One of the favored locations for subduction triggering is at passive margins, where sea floor abuts continental margins. Such passive margin collapse is problematic because the strength of the old, cold ocean lithosphere should prohibit it from bending under its own weight and sinking into the mantle. Some means of mechanical weakening of the passive margin are therefore necessary. Spontaneous and accumulated grain damage can allow for considerable lithospheric weakening and facilitate passive margin collapse. Grain damage is enhanced where mixing between mineral phases in lithospheric rocks occurs. Such mixing is driven both by compositional gradients associated with petrological heterogeneity and by the state of stress in the lithosphere. With lateral compressive stress imposed by ridge push in an opening ocean basin, bands of mixing and weakening can develop, become vertically oriented, and occupy a large portion of lithosphere after about 100 million y. These bands lead to anisotropic viscosity in the lithosphere that is strong to lateral forcing but weak to bending and sinking, thereby greatly facilitating passive margin collapse. 
    more » « less
  4. Abstract To understand the effects of secondary minerals on changes in the mechanical properties of upper mantle rocks due to phase mixing, we conducted high‐strain torsion experiments on aggregates of iron‐rich olivine + orthopyroxene (opx) with opx volume fractions offopx = 0.15, 0.26, and 0.35. For samples with larger amounts of opx,fopx = 0.26 and 0.35, the value of the stress exponent decreases with increasing strain fromn ≈ 3 for γ ≲ 5 ton ≈ 2 for 5 ≲ γ ≲ 25, indicating that the deformation mechanism changes as strain increases. In contrast, for samples withfopx = 0.15, the stress exponent is constant atn ≈ 3.3 for 1 ≲ γ ≲ 25, suggesting that no change in deformation mechanism occurs with increasing strain for samples with smaller amounts of opx. The microstructures of samples with larger amounts of opx provide insight into the change in deformation mechanism derived from the mechanical data. Elongated grains align subparallel to the shear direction for samples of all three compositions deformed to lower strains. However, strain weakening with grain size reduction and the formation of a thoroughly mixed, fine‐grained texture only develops in samples withfopx = 0.26 and 0.35 deformed to higher strains of γ ≳ 16. These mechanical and associated microstructural properties imply that rheological weakening due to phase mixing only occurs in the samples with largerfopx, which is an important constraint for understanding strain localization in the upper mantle of Earth. 
    more » « less
  5. Abstract To study the microstructural evolution of polymineralic rocks, we performed deformation experiments on two‐phase aggregates of olivine (Ol) + ferropericlase (Per) with periclase fractions (fPer) between 0.1 and 0.8. Additionally, single‐phase samples of both Ol and Per were deformed under the same experimental conditions to facilitate comparison of the microstructures in two‐phase and single‐phase materials. Each sample was deformed in torsion atT = 1523 K,P = 300 MPa at a constant strain rate up to a final shear strain of γ = 6 to 7. Microstructural developments, analyzed via electron backscatter diffraction (EBSD), indicate differences in both grain size and crystalline texture between single‐ and two‐phase samples. During deformation, grain size approximately doubled in our single‐phase samples of Ol and Per but remained unchanged or decreased in two‐phase samples. Zener‐pinning relationships fit to the mean grain sizes in each phase for samples with 0.1 ≤ fPer≤ 0.5 and for those with 0.8 ≥ fPer ≥ 0.5 demonstrate that the grain size of the primary phase is controlled by phase‐boundary pinning. Crystallographic preferred orientations, determined for both phases from EBSD data, are significantly weaker in the two‐phase materials than in the single‐phase materials. 
    more » « less