Multiscale modeling of upper mantle plasticity: From single-crystal rheology to multiphase aggregate deformation

We report a first application of an improved second-order (SO) viscoplastic self-consistent model for multiphase aggregates, applied to an olivine + diopside aggregate as analogue for a dry upper mantle peridotite deformed at 10 15 s 1 shear strain rate along a 20-Ma ocean geotherm. Beside known dislocation slip systems, this SO-model version accounts for an isotropic relaxation mechanism representing ‘diffusionrelated’ creep mechanisms in olivine. Slip-system critical resolved shear stress (CRSS) are evaluated in both phases – as functions of P, T, oxygen fugacity (fO2) and strain rate – from previously reported experimental data obtained on single crystals and first-principle calculations coupled with the Peierls–Nabarro model for crystal plasticity; and the isotropic-mechanism dependence on T and P matches that of Si selfdiffusion in olivine, while its relative activity is constrained by reported data. The model reproduces well the olivine and diopside lattice preferred orientations (LPO) produced experimentally and observed in naturally deformed rocks, as well as observed sensitivities of multiphase aggregate strength to the volume fraction of the hard phase (here diopside). It shows a significant weakening of olivine LPO with increasing depth, which results from the combined effects of the P-induced [100]/[001] dislocation-slip transition and the increasing activity with T of ‘diffusion-related’ creep. This work thus provides a first quantification of the respective effects of [100]/[001] slip transition and diffusion creep on the olivine LPO weakening inducing the seismic anisotropy attenuation observed in the upper mantle

Microstructures et rhéologie du manteau terrestre supérieur déduites d'une approche multi-échelle

The strongly anisotropic rheology of olivine polycrystals, associated to their microstructure, constitutes a key feature affecting the dynamics of the Earth's upper mantle. High pressure deformation experiments carried out on olivine single crystals under synchrotron radiation, together with estimations of lattice friction based on first-principle calculations, show a transition from easy [100] to easy [001] slips as pressure and temperature (thus depth) increases. We input these data at the slip system level into the second-order extension of the self-consistent scheme to assess microstructure evolution along a typical flow pattern beneath an oceanic spreading cente

Textures in deforming forsterite aggregates up to 8 GPa and 1673 K

We report results from axisymmetric deformation experiments carried out on forsterite aggregates in the deformation-DIA apparatus, at upper mantle pressures and temperatures (3.1–8.1 GPa, 1373–1673 K). We quantified the resulting lattice preferred orientations (LPO) and compare experimental observations with results from micromechanical modeling (viscoplastic second-order self-consistent model—SO). Up to 6 GPa (~185-km depth in the Earth), we observe a marked LPO consistent with a dominant slip in the (010) plane with one observation of a dominant [100] direction, suggesting that [100](010) slip system was strongly activated. At higher pressures (deeper depth), the LPO becomes less marked and more complex with no evidence of a dominant slip system, which we attribute to the activation of several concurrent slip systems. These results are consistent with the pressure-induced transition in the dominant slip system previously reported for olivine and forsterite. They are also consistent with the decrease in the seismic anisotropy amplitude observed in the Earth’s mantle at depth greater than ~200 km

Earth's oldest mantle fabrics indicate Eoarchaean subduction

The extension of subduction processes into the Eoarchaean era (4.0-3.6 Ga) is controversial. The oldest reported terrestrial olivine, from two dunite lenses within the ~3,720 Ma Isua supracrustal belt in Greenland, record a shape-preferred orientation of olivine crystals defining a weak foliation and a well-defined lattice-preferred orientation (LPO). [001] parallel to the maximum finite elongation direction and (010) perpendicular to the foliation plane define a B-type LPO. In the modern Earth such fabrics are associated with deformation of mantle rocks in the hanging wall of subduction systems; an interpretation supported by experiments. Here we show that the presence of B-type fabrics in the studied Isua dunites is consistent with a mantle origin and a supra-subduction mantle wedge setting, the latter supported by compositional data from nearby mafic rocks. Our results provide independent microstructural data consistent with the operation of Eoarchaean subduction and indicate that microstructural analyses of ancient ultramafic rocks provide a valuable record of Archaean geodynamics

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

Thermal weakening friction during seismic slip experiments and models with heat sources and sinks

Experiments that systematically explore rock friction under crustal earthquake conditions reveal that faults undergo abrupt dynamic weakening. Processes related to heating and weakening of fault surfaces have been invoked to explain pronounced velocity weakening. Both contact asperity temperature Ta and background temperature T of the slip zone evolve significantly during high-velocity slip due to heat sources (frictional work), heat sinks (e.g., latent heat of decomposition processes), and diffusion. Using carefully calibrated High-Velocity Rotary Friction experiments, we test the compatibility of thermal weakening models: (1) a model of friction based only on T in an extremely simplified, Arrhenius-like thermal dependence; (2) a flash heating model which accounts for the evolution of both V and T; (3) same but including heat sinks in the thermal balance; and (4) same but including the thermal dependence of diffusivity and heat capacity. All models reflect the experimental results but model (1) results in unrealistically low temperatures and model (2) reproduces the restrengthening phase only by modifying the parameters for each experimental condition. The presence of dissipative heat sinks in stage (3) significantly affects T and reflects on the friction, allowing a better joint fit of the initial weakening and final strength recovery across a range of experiments. Temperature is significantly altered by thermal dependence of (4). However, similar results can be obtained by (3) and (4) by adjusting the energy sinks. To compute temperature in this type of problem, we compare the efficiency of three different numerical approximations (finite difference, wavenumber summation, and discrete integral)