Dislocation interactions in olivine control postseismic creep of the upper mantle.

Changes in stress applied to mantle rocks, such as those imposed by earthquakes, commonly induce a period of transient creep, which is often modelled based on stress transfer among slip systems due to grain interactions. However, recent experiments have demonstrated that the accumulation of stresses among dislocations is the dominant cause of strain hardening in olivine at temperatures ≤600 °C, raising the question of whether the same process contributes to transient creep at higher temperatures. Here, we demonstrate that olivine samples deformed at 25 °C or 1150-1250 °C both preserve stress heterogeneities of ~1 GPa that are imparted by dislocations and have correlation lengths of ~1 μm. The similar stress distributions formed at these different temperatures indicate that accumulation of stresses among dislocations also provides a contribution to transient creep at high temperatures. The results motivate a new generation of models that capture these intragranular processes and may refine predictions of evolving mantle viscosity over the earthquake cycle

Seismic Anisotropy of Temperate Ice in Polar Ice Sheets

We present a series of simple shear numerical simulations of dynamic recrystallization of two‐phase nonlinear viscous materials that represent temperate ice. First, we investigate the effect of the presence of water on the resulting microstructures and, second, how water influences on P wave (Vp) and fast S wave (Vs) velocities. Regardless the water percentage, all simulations evolve from a random fabric to a vertical single maximum. For a purely solid aggregate, the highest Vp quickly aligns with the maximum c‐axis orientation. At the same time, the maximum c‐axis development reduces Vs in this orientation. When water is present, the developed maximum c‐axis orientation is less intense, which results in lower Vp and Vs. At high percentage of water, Vp does not align with the maximum c‐axis orientation. If the bulk modulus of ice is assumed for the water phase (i.e., implying that water is at high pressure), we find a remarkable decrease of Vs while Vp remains close to the value for purely solid ice. These results suggest that the decrease in Vs observed at the base of the ice sheets could be explained by the presence of water at elevated pressure, which would reside in isolated pockets at grain triple junctions. Under these conditions water would not favor sliding between ice grains. However, if we consider that deformation dominates over recrystallization, water pockets get continuously stretched, allowing water films to be located at grain boundaries. This configuration would modify and even overprint the maximum c‐axis‐dependent orientation and the magnitude of seismic anisotropy

Dislocation interactions in olivine control postseismic creep of the upper mantle

Abstract: Changes in stress applied to mantle rocks, such as those imposed by earthquakes, commonly induce a period of transient creep, which is often modelled based on stress transfer among slip systems due to grain interactions. However, recent experiments have demonstrated that the accumulation of stresses among dislocations is the dominant cause of strain hardening in olivine at temperatures ≤600 °C, raising the question of whether the same process contributes to transient creep at higher temperatures. Here, we demonstrate that olivine samples deformed at 25 °C or 1150–1250 °C both preserve stress heterogeneities of ~1 GPa that are imparted by dislocations and have correlation lengths of ~1 μm. The similar stress distributions formed at these different temperatures indicate that accumulation of stresses among dislocations also provides a contribution to transient creep at high temperatures. The results motivate a new generation of models that capture these intragranular processes and may refine predictions of evolving mantle viscosity over the earthquake cycle

Analytical Parametrization of Self-Consistent Polycrystal Mechanics: Fast Calculation of Upper Mantle Anisotropy

Progressive deformation of upper mantle rocks via dislocation creep causes their constituent crystals to take on a non-random orientation distribution (crystallographic preferred orientation or CPO) whose observable signatures include shear-wave splitting and azimuthal dependence of surface wave speeds. Comparison of these signatures with mantle flow models thus allows mantle dynamics to be unraveled on global and regional scales. However, existing self-consistent models of CPO evolution are computationally expensive when used in 3-D and/or time-dependent convection models. Here we propose a new method, called ANPAR, which is based on an analytical parameterisation of the crystallographic spin predicted by the second-order (SO) self-consistent theory. Our parameterisation runs approximately 2-6x10^4 times faster than the SO model and fits its predictions for CPO and crystallographic spin with a variance reduction > 99%. We illustrate the ANPAR model predictions for the deformation of olivine with three dominant slip systems, (010)[100], (001)[100] and (010)[001], for three uniform deformations (uniaxial compression, pure shear, simple shear) and for a corner-flow model of a spreading mid-ocean ridge

Can changes in deformation regimes be inferred from crystallographic preferred orientations in polar ice?

Creep due to ice flow is generally thought to be the main cause for the formation of crystallographic preferred orientations (CPOs) in polycrystalline anisotropic ice. However, linking the development of CPOs to the ice flow history requires a proper understanding of the ice aggregate's microstructural response to flow transitions. In this contribution the influence of ice deformation history on the CPO development is investigated by means of full-field numerical simulations at the microscale. We simulate the CPO evolution of polycrystalline ice under combinations of two consecutive deformation events up to high strain, using the code VPFFT (visco-plastic fast Fourier transform algorithm) within ELLE. A volume of ice is first deformed under coaxial boundary conditions, which results in a CPO. The sample is then subjected to different boundary conditions (coaxial or non-coaxial) in order to observe how the deformation regime switch impacts the CPO. The model results indicate that the second flow event tends to destroy the first, inherited fabric with a range of transitional fabrics. However, the transition is slow when crystallographic axes are critically oriented with respect to the second imposed regime. Therefore, interpretations of past deformation events from observed CPOs must be carried out with caution, particularly in areas with complex deformation histories

Can changes in deformation regimes be inferred from crystallographic preferred orientations in polar ice?

Creep due to ice flow is generally thought to be the main cause for the formation of crystallographic preferred orientations (CPOs) in polycrystalline anisotropic ice. However, linking the development of CPOs to the ice flow history requires a proper understanding of the ice aggregate's microstructural response to flow transitions. In this contribution the influence of ice deformation history on the CPO development is investigated by means of full-field numerical simulations at the microscale. We simulate the CPO evolution of polycrystalline ice under combinations of two consecutive deformation events up to high strain, using the code VPFFT (visco-plastic fast Fourier transform algorithm) within ELLE. A volume of ice is first deformed under coaxial boundary conditions, which results in a CPO. The sample is then subjected to different boundary conditions (coaxial or non-coaxial) in order to observe how the deformation regime switch impacts the CPO. The model results indicate that the second flow event tends to destroy the first, inherited fabric with a range of transitional fabrics. However, the transition is slow when crystallographic axes are critically oriented with respect to the second imposed regime. Therefore, interpretations of past deformation events from observed CPOs must be carried out with caution, particularly in areas with complex deformation histories