Weighted Burgers Vector analysis of orientation fields from high-angular resolution electron backscatter diffraction

The Weighted Burgers Vector (WBV) method can extract information about dislocation types and densities present in distorted crystalline materials from electron backscatter diffraction (EBSD) maps, using no assumptions about which slip systems might be present. Furthermore, high-angular resolution EBSD (HR-EBSD) uses a cross-correlation procedure to increase the angular precision of EBSD measurements by an order of magnitude compared to conventional EBSD. However, the WBV technique has not previously been applied to HR-EBSD data and therefore it remains unclear as to which low-angle substructures can be reliably characterised by WBV analysis of conventional EBSD data and which require additional HR-EBSD processing. To establish some practical examples that can be used to guide future data-acquisition strategies, we compare the output of the WBV method when applied to conventional EBSD data and HR-EBSD data collected from the most common minerals in Earth's lower crust (plagioclase feldspar) and upper mantle (olivine). The results demonstrate that HR-EBSD and WBV processing are complementary techniques. The increase in angular precision achieved with HR-EBSD processing allows low-angle (on the order of 0.1°) structures, which are obscured by noise in conventional EBSD data, to be analyzed quantitatively using the WBV method. Combining the WBV and HR-EBSD methods increases the precision of calculated WBV directions, which is essential when using information about active slip systems to infer likely deformation mechanisms from naturally deformed microstructures. This increase in precision is particularly important for low-symmetry crystals, such as plagioclase, that have a wide range of available slip systems that vary in relative activity with changing pressure, temperature and differential stress. Because WBV directions are calculated using no assumptions about which slip systems may be present, combining this technique with HR-EBSD to refine the precision of lattice orientation gradients is ideal for investigating complex natural materials with unknown deformation histories.

Elastic wave propagation in anisotropic polycrystals: inferring physical properties of glacier ice

An optimization problem is proposed for inferring physical properties of polycrystals given ultrasonic (elastic) wave velocity measurements, made across multiple sample orientations. The feasibility of the method is demonstrated by inferring both the effective grain elastic parameters and the grain c -axis orientation distribution function (ODF) of ice-core samples from Priestley glacier, Antarctica. The method relies on expanding the ODF in terms of a spherical harmonic series, which allows for a non-parametric estimation of the sample ODF. Moreover, any linear combination of the Voigt (strain) and Reuss (stress) homogenization scheme is allowed, although for glacier ice, the exact choice is found to matter little for bulk elastic behaviour, and thus for inferred physical properties, too. Finally, the accuracy of the inferred grain elastic parameters is discussed, including the well-posedness and shortcomings of the inverse problem, relevant for future adoptions in glaciology, geology and elsewhere.

Explaining the dependence of M-site diffusion in forsterite on silica activity: a density functional theory approach

Experimentally, silica activity (aSiO₂) has been shown to have an effect on Mg diffusion in forsterite, but no fully satisfactory mechanism has yet been proposed. We calculated the effects of aSiO₂ and aluminium content (the main contaminant in some recent experimental studies), and their co-effect, on Mg diffusion in forsterite, using thermodynamic minimisations of defect formation energies [calculated using density functional theory (DFT)] and a Monte-Carlo diffusion model. These two variables, in isolation, do not appreciably change the defect concentrations of forsterite and thus do not affect the diffusivity of Mg. However, when elevated together, they cause large increases in the Mg vacancy content and thus can increase the Mg diffusivity by one to six orders of magnitude depending on temperature, with little pressure dependence. This effect is largely independent of Al₂O₃ concentration above ~ 1 wt. ppm, and thus, for all practical purposes, should occur wherever forsterite is in the presence of enstatite. It is also largely dependent upon configurational entropy and is thus highly sensitive to the chemistry of the crystal. A low concentration of structurally bound hydroxyl groups at low temperatures (1000 K) suppresses this effect in pure forsterite, but it is likely robust in the presence of water either when alternative water sinks (such as Ti or Fe) are present, or at high temperatures (> 1500 K). This effect is also robust in the presence of ferrous iron (or other substitutional Mg defects) at all temperatures. Fe₂O₃ can operate like Al₂O₃ in this reaction and should enhance its effect. These findings explain the experimentally observed dependency of Mg diffusion of aSiO₂, and elucidate how chemical activity variations in both experiments and natural settings could affect not only the diffusivity of Mg in forsterite, but of olivine-hosted cations in general.

Laboratory investigation of mechanisms for phase mixing in olivine + ferropericlase aggregates

To investigate the role of grain boundary pinning and the mechanisms by which phase mixing occurs during deformation of polymineralic rocks, we conducted high-strain torsion experiments on samples consisting of olivine plus 30 vol% ferropericlase. Experiments were performed in a gas-medium deformation apparatus at 1524 K and 300 MPa. Samples were deformed to outer radius shear strains of up to γ(R) = 14.1. The value of the stress exponent and the small grain sizes of our samples indicate that our two-phase material deformed by dislocation-accommodated grain boundary sliding. In samples deformed to 1 < γ < 7, elongated clusters of ferropericlase grains form thin layers in the olivine matrix, and small grains of ferropericlase appear at olivine grain boundaries and three- and four-grain junctions. By γ ≈ 14, a well-distributed mixture of small ferropericlase grains among the olivine grains developed. Microstructures exhibit similarities to both mechanical and chemical models proposed to describe the processes leading to phase mixing. Our results provide evidence for grain size reduction during phase mixing that results in a grain size significantly smaller than the value predicted by the single-phase recrystallization piezometer for olivine. Thus, phase mixing provides a mechanism for the persistent weakening of rocks that is important for developing and maintaining shear zones necessary for plate tectonics.This article is part of a discussion meeting issue 'Earth dynamics and the development of plate tectonics'.

Redox-influenced seismic properties of upper-mantle olivine

Lateral variations of seismic wave speeds and attenuation (dissipation of strain energy) in the Earth's upper mantle have the potential to map key characteristics such as temperature, major-element composition, melt fraction and water content. The inversion of these data into meaningful representations of physical properties requires a robust understanding of the micromechanical processes that affect the propagation of seismic waves. Structurally bound water (hydroxyl) is believed to affect seismic properties but this has yet to be experimentally quantified. Here we present a comprehensive low-frequency forced-oscillation assessment of the seismic properties of olivine as a function of water content within the under-saturated regime that is relevant to the Earth's interior. Our results demonstrate that wave speeds and attenuation are in fact strikingly insensitive to water content. Rather, the redox conditions imposed by the choice of metal sleeving, and the associated defect chemistry, appear to have a substantial influence on the seismic properties. These findings suggest that elevated water contents are not responsible for low-velocity or high-attenuation structures in the upper mantle. Instead, the high attenuation observed in hydrous and oxidized regions of the upper mantle (such as above subduction zones) may reflect the prevailing oxygen fugacity. In addition, these data provide no support for the hypothesis whereby a sharp lithosphere-asthenosphere boundary is explained by enhanced grain boundary sliding in the presence of water.

Size effects resolve discrepancies in 40 years of work on low-temperature plasticity in olivine

The strength of olivine at low temperatures and high stresses in Earth's lithospheric mantle exerts a critical control on many geodynamic processes, including lithospheric flexure and the formation of plate boundaries. Unfortunately, laboratory-derived values of the strength of olivine at lithospheric conditions are highly variable and significantly disagree with those inferred from geophysical observations. We demonstrate via nanoindentation that the strength of olivine depends on the length scale of deformation, with experiments on smaller volumes of material exhibiting larger yield stresses. This "size effect" resolves discrepancies among previous measurements of olivine strength using other techniques. It also corroborates the most recent flow law for olivine, which proposes a much weaker lithospheric mantle than previously estimated, thus bringing experimental measurements into closer alignment with geophysical constraints. Further implications include an increased difficulty of activating plasticity in cold, fine-grained shear zones and an impact on the evolution of fault surface roughness due to the size-dependent deformation of nanometer- to micrometer-sized asperities.

Inferring nonlinear mantle rheology from the shape of the Hawaiian swell

The convective circulation generated within the Earth's mantle by buoyancy forces of thermal and compositional origin is intimately controlled by the rheology of the rocks that compose it. These can deform either by the diffusion of point defects (diffusion creep, with a linear relationship between strain rate and stress) or by the movement of intracrystalline dislocations (nonlinear dislocation creep). However, there is still no reliable map showing where in the mantle each of these mechanisms is dominant, and so it is important to identify regions where the operative mechanism can be inferred directly from surface geophysical observations. Here we identify a new observable quantity--the rate of downstream decay of the anomalous seafloor topography (swell) produced by a mantle plume--which depends only on the value of the exponent in the strain rate versus stress relationship that defines the difference between diffusion and dislocation creep. Comparison of the Hawaiian swell topography with the predictions of a simple fluid mechanical model shows that the swell shape is poorly explained by diffusion creep, and requires a dislocation creep rheology. The rheology predicted by the model is reasonably consistent with laboratory deformation data for both olivine and clinopyroxene, suggesting that the source of Hawaiian lavas could contain either or both of these components.

Rheology of the upper mantle: a synthesis

Rheological properties of the upper mantle of the Earth play an important role in the dynamics of the lithosphere and asthenosphere. However, such fundamental issues as the dominant mechanisms of flow have not been well resolved. A synthesis of laboratory studies and geophysical and geological observations shows that transitions between diffusion and dislocation creep likely occur in the Earth's upper mantle. The hot and shallow upper mantle flows by dislocation creep, whereas cold and shallow or deep upper mantle may flow by diffusion creep. When the stress increases, grain size is reduced and the upper mantle near the transition between these two regimes is weakened. Consequently, deformation is localized and the upper mantle is decoupled mechanically near these depths.

Dynamic recrystallization and strain softening of olivine aggregates in the laboratory and the lithosphere

The effects of dynamic recrystallization on the deformation mechanisms and rheology of olivine aggregates in the laboratory and the lithosphere are reviewed in this paper. The low-strain rheology of olivine is well documented; however, deformation in the lithosphere often involves large strains. Large strain experiments show that recrystallization can result in both hardening and softening during deformation. Moderate strain softening in experimental shear and torsion can be explained by the operation of dislocation-accommodated grain boundary sliding in bands of fine recrystallized grains.

Data on the temperature dependence of recrystallized grain size are needed to extrapolate the effects of dynamic recrystallization to the lithosphere. Theories of dynamic recrystallization suggest that grain size is strongly stress dependent and moderately temperature dependent. A re-analysis of experimental grain size data indicates that the recrystallized grain size is temperature independent for olivine aggregates with low water content (<300 ppm H/Si).

Rheological regime maps have been constructed for the lithospheric mantle. The maps suggest that grain size sensitive power law creep, involving both grain boundary sliding and dislocation creep, will produce strong strain softening, greater than found so far in experimental studies, in dry and wet lithosphere shear zones.

Mantle Phase Changes and Deep-Earthquake Faulting in Subducting Lithosphere

Inclined zones of earthquakes are the primary expression of lithosphere subduction. A distinct deep population of subduction-zone earthquakes occurs at depths of 350 to 690 kilometers. At those depths ordinary brittle fracture and frictional sliding, the faulting processes of shallow earthquakes, are not expected. A fresh understanding of these deep earthquakes comes from developments in several areas of experimental and theoretical geophysics, including the discovery and characterization of transformational faulting, a shear instability connected with localized phase transformations under nonhydrostatic stress. These developments support the hypothesis that deep earthquakes represent transformational faulting in a wedge of olivine-rich peridotite that is likely to persist metastably in coldest plate interiors to depths as great as 690 km. Predictions based on this deep structure of mantle phase changes are consistent with the global depth distribution of deep earthquakes, the maximum depths of earthquakes in individual subductions zones, and key source characteristics of deep events.