International audienceReactivation of structures inherited from previous collisional or rifting events, especially lithospheric-scale faults, is a major feature of plate tectonics. Its expression ranges from continental break-up along ancient collisional belts(1,2) to linear arrays of intraplate magmatism and seismicity(3,4). Here we use multiscale numerical models to show that this reactivation can result from an anisotropic mechanical behaviour of the lithospheric mantle due to an inherited preferred orientation of olivine crystals. We explicitly consider an evolving anisotropic viscosity controlled by the orientation of olivine crystals in the mantle. We find that strain is localized in domains where shear stresses on the inherited mantle fabric are high, and that this leads to shearing parallel to the inherited fabric. During rifting, structural reactivation induced by anisotropy results in oblique extension, followed by either normal extension or failure. Our results suggest that anisotropic viscosity in the lithospheric mantle controls the location and orientation of intraplate deformation zones that may evolve into new plate boundaries, and causes long-lived lithospheric-scale wrench faults, contributing to the toroidal component of plate motions on Earth

S U M M A R YThe Earth's non-hydrostatic gravity field, or geoid, provides a first-order constraint on mantle density structure and dynamics. Geodynamic models for the geoid have proliferated since the advent of seismic mapping of mantle heterogeneity structure (tomography) because the geoid offers perhaps the best-measured independent constraint on mantle density heterogeneity. However, dynamic geoid models involve a number of questionable physical assumptions and uncertainties whose effects need to be evaluated before geodynamic inferences based upon the geoid can be considered sound. Troubling issues include the appropriate surface boundary conditions (free-slip, no-slip, plates?) and parametrization of radial viscosity variations (how many layers can be resolved?), in addition to lateral viscosity variations, possible chemical layering of the mantle, phase transitions, etc. There are also uncertainties in the density heterogeneity models used as 'input' to dynamic geoid models, most of which are derived from seismic tomography and require weakly constrained, empirical conversion factors to go from seismic velocity variations to density variations. Here we address several of the most straightforward problems inherent in geoid modelling, namely the issues of viscosity structure resolution, uncertainties in appropriate boundary conditions, and differences among mantle heterogeneity models. A robust feature of all models is a lower-mantle viscosity at least a factor of 30 greater than that of the upper mantle, but there is little resolution with regard to finer details such as lithospheric or uppermost mantle ('low-viscosity zone') viscosity. Ironically, free-slip boundary conditions result in the best fits to the geoid in all cases, but all boundary conditions exhibit predictable trade-offs with the uppermostmantle viscosity. Models with a single viscosity layer representing the lower mantle yield similar dynamic topography estimates of the order of 700-1000 m in amplitude, regardless of the finer details of upper-mantle viscosity structure, boundary conditions or input heterogeneity models. Comparing mantle heterogeneity models based on two independent seismological determinations (Harvard and Berkeley models) and on the history of subduction, we find that these models are virtually indistinguishable regarding inferences of mantle viscosity structure and amplitude of dynamic topography, and in terms of the effects of different boundary conditions. Uncertainties concerning which type of boundary condition is appropriate are much more important than which mantle heterogeneity model is chosen. Given other uncertainties in modelling the geoid, particularly the strong effects due to lateral viscosity variations for intermediate (< 10 000 km) wavelengths, we conclude that the class of dynamic geoid models explored so far cannot reliably elucidate the details of upper-mantle viscosity structure.

[1] We use forward models based on recent high-pressure experimental data on mantle minerals to predict the seismic anisotropy produced by plastic strain of orthorhombic wadsleyite, the dominant mineral in the upper transition zone. These models predict a weak seismic anisotropy for a polycrystal of pyrolitic composition (60% wadsleyite, 40% garnet) at transition zone conditions: $2% for P and $1% for S waves for a shear strain of 1. Both P and S wave anisotropy patterns show an orthorhombic symmetry. P waves propagate faster at low angle to the shear direction and slower at high angle to the shear plane. S wave anisotropy is characterized by faster propagation of waves polarized at low angle to the shear direction. Horizontal shearing results therefore in higher velocities for horizontally propagating P waves (PH ) and horizontally polarized S waves (SH ), as well as in weak azimuthal variation of SV and SH velocities. On the other hand, vertical flow leads to higher velocities for vertically propagating P waves (PV ) and vertically polarized S waves (SV) and to a weak azimuthal variation of SV velocity but to a roughly constant SH velocity. Analysis of global observations of seismic anisotropy in the transition zone in the light of these models supports dominant horizontal flow in the uppermost transition zone, in agreement with predictions of geodynamical models that explicitly introduce phase transitions.