skip to main content

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):
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Geophysical Journal International
Page Range / eLocation ID:
2205 to 2230
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. Abstract

    Ultramylonites—intensely deformed rocks with fine grain sizes and well‐mixed mineral phases—are thought to be a key component of Earth‐like plate tectonics, because coupled phase mixing and grain boundary pinning enable rocks to deform by grain‐size‐sensitive, self‐softening creep mechanisms over long geologic timescales. In isoviscous two‐phase composites, “geometric” phase mixing occurs via the sequential formation, attenuation (stretching), and disaggregation of compositional layering. However, the effects of viscosity contrast on the mechanisms and timescales for geometric mixing are poorly understood. Here, we describe a series of high‐strain torsion experiments on nonisoviscous calcite‐fluorite composites (viscosity contrast,ηca/ηfl≈ 200) at 500°C, 0.75 GPa confining pressure, and 10−6–10−4 s−1shear strain rate. At low to intermediate shear strains (γ ≤ 10), polycrystalline domains of the individual phases become sheared and form compositional layering. As layering develops, strain localizes into the weaker phase, fluorite. Strain partitioning impedes mixing by reducing the rate at which the stronger (calcite) layers deform, attenuate, and disaggregate. Even at very large shear strains (γ ≥ 50), grain‐scale mixing is limited, and thick compositional layers are preserved. Our experiments (1) demonstrate that viscosity contrasts impede mechanical phase mixing and (2) highlight the relative inefficiency of mechanical mixing. Nevertheless, by employing laboratory flow laws, we show that “ideal” conditions for mechanical phase mixing may be found in the wet middle to lower continental crust and in the dry mantle lithosphere, where quartz‐feldspar and olivine‐pyroxene viscosity contrasts are minimized, respectively.

    more » « less
  4. 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
  5. Constraints on the state of stress in the lithosphere are fundamental to understanding a breadth of geological phenomena. Paleo-stresses are generally estimated using microstructural elements for which there are experimentally calibrated relationships with applied stress, with an emphasis on recrystallised grain-size piezometers. However, it is often difficult to clearly distinguish newly recrystallised grains from the relict matrix. Furthermore, these grain-size piezometers are only applicable to rocks consisting of a single mineral. An alternative proxy for paleo-stress in polymineralic rocks is the average subgrain size. Unfortunately, estimates of subgrain size differ significantly among different measurement methods, and therefore, piezometers must be individually calibrated for the method used. Existing subgrain-size piezometers are based on calibrations using optical or transmission electron microscopy. We use electron backscatter diffraction (EBSD), a common method of subgrain-boundary characterisation, to calibrate subgrain-size piezometers for both olivine and quartz. To test the application of our olivine subgrain-size piezometer to polymineralic rocks, we deformed synthetic mixtures of olivine and orthopyroxene. Experiments were conducted using a Deformation-DIA apparatus at beamline 6BM-B Advanced Photon Source, Argonne National Laboratory. These experiments offer the unique possibility of simultaneously deforming the sample and measuring the average stresses within each phase using X-ray diffraction, before applying subgrain-size piezometry to the recovered samples. The results provide tests of (1) the manner in which stress is partitioned between phases, (2) whether the stresses measured in each phase by X-ray diffraction are comparable to those estimated by subgrain-size piezometry, and (3) whether stresses from subgrain piezometry can be used to estimate the macroscopic average applied stress. Stresses estimated from X-ray diffraction agree well with those made from subgrain-size piezometry in both monomineralic and polymineralic samples. In harzburgites, average stresses are similar in both phases and indicate that in this system, subgrain-size piezometric measurements from a single phase can be used to estimate the bulk stress. 
    more » « less