Pressure sensitivity of olivine slip systems and seismic anisotropy of Earth's upper mantle

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.

Fluid-assisted grain size reduction leads to strain localization in oceanic transform faults

Oceanic Transform Faults are major plate boundaries representing the most seismogenic part of the mid ocean ridge system. Nonetheless, their structure and deformation mechanisms at depth are largely unknown due to rare exposures of deep sections. Here we study the mineral fabric of deformed mantle peridotites - ultramafic mylonites - collected from the transpressive Atobá ridge, along the northern fault of the St. Paul transform system in the Equatorial Atlantic Ocean. We show that, at pressure and temperature conditions of the lower oceanic lithosphere, the dominant deformation mechanism is fluid-assisted dissolution-precipitation creep. Grain size reduction during deformation is enhanced by dissolution of coarser pyroxene grains in presence of fluid and contextual precipitation of small interstitial ones, leading to strain localization at lower stresses than dislocation creep. This mechanism potentially represents the dominant weakening factor in the oceanic lithosphere and a main driver for the onset and maintenance of oceanic transform faults.

Glacial isostatic adjustment: physical models and observational constraints

By far the most prescient insights into the interior structure of the planet have been provided on the basis of elastic wave seismology. Analysis of the travel times of shear or compression wave phases excited by individual earthquakes, or through analysis of the elastic gravitational free oscillations that individual earthquakes of sufficiently large magnitude may excite, has been the central focus of Earth physics research for more than a century. Unfortunately, data provide no information that is directly relevant to understanding the solid state 'flow' of the polycrystalline outer 'mantle' shell of the planet that is involved in the thermally driven convective circulation that is responsible for powering the 'drift' of the continents and which controls the rate of planetary cooling on long timescales. For this reason, there has been an increasing focus on the understanding of physical phenomenology that is unambiguously associated with mantle flow processes that are distinct from those directly associated with the convective circulation itself. This paper reviews the past many decades of work that has been invested in understanding the most important of such processes, namely that which has come to be referred to as 'glacial isostatic adjustment' (GIA). This process concerns the response of the planet to the loading and unloading of the high latitude continents by the massive accumulations of glacial ice that have occurred with almost metronomic regularity over the most recent million years of Earth history. Forced by the impact of gravitationaln-body effects on the geometry of Earth's orbit around the Sun through the impact upon the terrestrial regime of received solar insolation, these surface mass loads on the continents have left indelible records of their occurrence in the 'Earth system' consisting of the oceans, continents, and the great polar ice sheets on Greenland and Antarctica themselves. Although this ice-age phenomenology has been clearly recognized since early in the last century, it was for over 50 years considered to be no more than an interesting curiosity, the understanding of which remained on the periphery of the theoretical physics of the Earth. This was the case in part because no globally applicable theory was available that could be applied to rigorously interpret the observations. Equally important to understanding the scientific lethargy that held back the understanding of this phenomenon involving mantle flow processes was the lack of appreciation of the wide range of observations that were in fact related to GIA physics. This paper is devoted to a review of the global theories of the GIA process that have since been developed as a means of interpreting the extensive variety of observations that are now recognized as being involved in the response of the planet to the loading and unloading of its surface by glacial ice. The paper will also provide examples of the further analyses of Earth physics and climate related processes that applications of the modern theoretical structures have enabled.

Using Multigrain Crystallography to Explore the Microstructural Evolution of the α-Olivine to γ-Ringwoodite Transformation and ε-Mg2SiO4 at High Pressure and Temperature

The introduction of multigrain crystallography (MGC) applied in a laser-heated diamond anvil cell (LH-DAC) using synchrotron X-rays has provided a new path to investigate the microstructural evolution of materials at extreme conditions, allowing for simultaneous investigations of phase identification, strain state determination, and orientation relations across phase transitions in a single experiment. Here, we applied this method to a sample of San Carlos olivine beginning at ambient conditions and through the α-olivine → γ-ringwoodite phase transition. At ambient temperatures, by measuring the evolution of individual Bragg reflections, olivine shows profuse angular streaking consistent with the onset of yielding at a measured stress of ~1.5 GPa, considerably lower than previously reported, which may have implications for mantle evolution. Furthermore, γ-ringwoodite phase was found to nucleate as micron to sub-micron grains imbedded with small amounts of a secondary phase at 15 GPa and 1000 °C. Using MGC, we were able to extract and refine individual crystallites of the secondary unknown phase where it was found to have a structure consistent with the ε-phase previously described in chondritic meteorites.

Stress-induced amorphization triggers deformation in the lithospheric mantle

The mechanical properties of olivine-rich rocks are key to determining the mechanical coupling between Earth's lithosphere and asthenosphere. In crystalline materials, the motion of crystal defects is fundamental to plastic flow^1-4. However, because the main constituent of olivine-rich rocks does not have enough slip systems, additional deformation mechanisms are needed to satisfy strain conditions. Experimental studies have suggested a non-Newtonian, grain-size-sensitive mechanism in olivine involving grain-boundary sliding^5,6. However, very few microstructural investigations have been conducted on grain-boundary sliding, and there is no consensus on whether a single or multiple physical mechanisms are at play. Most importantly, there are no theoretical frameworks for incorporating the mechanics of grain boundaries in polycrystalline plasticity models. Here we identify a mechanism for deformation at grain boundaries in olivine-rich rocks. We show that, in forsterite, amorphization takes place at grain boundaries under stress and that the onset of ductility of olivine-rich rocks is due to the activation of grain-boundary mobility in these amorphous layers. This mechanism could trigger plastic processes in the deep Earth, where high-stress conditions are encountered (for example, at the brittle-plastic transition). Our proposed mechanism is especially relevant at the lithosphere-asthenosphere boundary, where olivine reaches the glass transition temperature, triggering a decrease in its viscosity and thus promoting grain-boundary sliding.

Exceptional increase in the creep life of magnesium rare-earth alloys due to localized bond stiffening

Several recent papers report spectacular, and unexpected, order of magnitude improvement in creep life of alloys upon adding small amounts of elements like zinc. This microalloying effect raises fundamental questions regarding creep deformation mechanisms. Here, using atomic-scale characterization and first principles calculations, we attribute the 600% increase in creep life in a prototypical Mg-rare earth (RE)-Zn alloy to multiple mechanisms caused by RE-Zn bonding-stabilization of a large volume fraction of strengthening precipitates on slip planes, increase in vacancy diffusion barrier, reduction in activated cross-slip, and enhancement of covalent character and bond strength around Zn solutes along the c-axis of Mg. We report that increased vacancy diffusion barrier, which correlates with the observed 25% increase in interplanar bond stiffness, primarily enhances the high-temperature creep life. Thus, we demonstrate that an approach of local, randomized tailoring of bond stiffness via microalloying enhances creep performance of alloys.

Seismic anisotropy evidence for dehydration embrittlement triggering intermediate-depth earthquakes

It has been proposed that dehydration embrittlement of hydrous materials can trigger intermediate-depth earthquakes and form a double seismic zone in a subducting slab. Seismic anisotropy may provide a possible insight into intermediate-depth intraslab seismicity, because anisotropic properties of minerals change with varying water distribution, temperature and pressure. Here we present a high-resolution model of P-wave radial anisotropy tomography of the Japan subduction zone down to ~400 km depth, which is obtained using a large number of arrival-time data of local earthquakes and teleseismic events. Our results reveal a close correlation between the pattern of intermediate-depth seismicity and anisotropic structures. The seismicity occurs in portions of the Pacific and Philippine Sea slabs where positive radial anisotropy (i.e., horizontal velocity being faster than vertical one) dominates due to dehydration, whereas the inferred anhydrous parts of the slabs are found to be aseismic where negative radial anisotropy (i.e., vertical velocity being faster than horizontal one) dominates. Our anisotropic results suggest that intermediate-depth earthquakes in Japan could be triggered by dehydration embrittlement of hydrous minerals in the subducting slabs.

In situ analysis of texture development from sinusoidal stress at high pressure and temperature

Here, we present a new experimental protocol to investigate the relationship between texture, plastic strain, and the mechanisms of plastic deformation at high pressure and temperature. The method utilizes synchrotron X-ray radiation as the probing tool, coupled with a large-volume high pressure deformation device (D-DIA). The intensity of X-ray diffraction peaks within the spectrum of the sample is used for sampling texture development in situ. The unique feature of this study is given by the sinusoidal variation of the intensity when a sinusoidal strain is applied to the sample. For a sample of magnesium oxide at elevated pressure and temperature, we demonstrate observations that are consistent with elasto-plastic models for texture development and for diffraction-peak measurements of apparent stress. The sinusoidal strain magnitude was 3%.

Messengers from the deep: Fossil wadsleyite-chromite microstructures from the Mantle Transition Zone

Investigations of the Mantle Transition Zone (MTZ; 410-660 km deep) by deformation experiments and geophysical methods suggest that the MTZ has distinct rheological properties, but their exact cause is still unclear due to the lack of natural samples. Here we present the first direct evidence for crystal-plastic deformation by dislocation creep in the MTZ using a chromitite from the Luobusa peridotite (E. Tibet). Chromite grains show exsolution of diopside and SiO2, suggesting previous equilibration in the MTZ. Electron backscattered diffraction (EBSD) analysis reveals that olivine grains co-existing with exsolved phases inside chromite grains and occurring on chromite grain boundaries have a single pronounced crystallographic preferred orientation (CPO). This suggests that olivine preserves the CPO of a high-pressure polymorph (wadsleyite) before the high-pressure polymorph of chromite began to invert and exsolve. Chromite also shows a significant CPO. Thus, the fine-grained high-pressure phases were deformed by dislocation creep in the MTZ. Grain growth in inverted chromite produced an equilibrated microstructure during exhumation to the surface, masking at first sight its MTZ deformation history. These unique observations provide a window into the deep Earth, and constraints for interpreting geophysical signals and their geodynamic implications in a geologically robust context.

Seismic signature of the Alpine indentation, evidence from the Eastern Alps

The type of collision between the European and the Adriatic plates in the easternmost Alps is one of the most interesting questions regarding the Alpine evolution. Tectonic processes such as compression, escape and uplift are interconnected and shape this area. We can understand these ongoing processes better, if we look for signs of the deformation within the Earth's deep crust of the region. By collecting records from permanent and temporary seismic networks, we assemble a receiver function dataset, and analyze it with the aim of giving new insights on the structure of the lower crust and of the shallow portion of the upper mantle, which are inaccessible to direct observation. Imaging is accomplished by performing common conversion depth stacks along three profiles that crosscut the Eastern Alpine orogen, and allow isolating features consistently persistent in the area. The study shows a moderately flat Moho underlying a seismically anisotropic middle-lower crust from the Southern Alps to the Austroalpine nappes. The spatial progression of anisotropic axes reflects the orientation of the relative motion and of the stress field detected at the surface. These observations suggest that distributed deformation is due to the effect of the Alpine indentation. In the shallow upper mantle right below the Moho interface, a further anisotropic layer is recognized, extended from the Bohemian Massif to the Northern Calcareous Alps.

Energy dissipation of materials at high pressure and high temperature

We report an experimental method to study the anelastic properties of materials at high pressure and high temperature. The multianvil high pressure deformation device, used to apply a cyclic loading force onto the sample, can reach 15 GPa and 2000 K. A synchrotron x-ray radiation source provides time resolved images of the sample and reference material. The images yield stress and strain as a function of time; stresses are derived from the reference material, and strains from the sample. This method has been tested by applying a sinusoidal stress at megahertz to hertz frequency on a San Carlos olivine specimen at 5 GPa and up to 2000 K. Strain as small as 10(-5) can be resolved. We have obtained experimental results which exhibit resolvable attenuation factor (Q(-1)) and shear modulus (M) at deep Earth conditions. These results are in quantitative agreement with previously reported lower pressure data and suggest that temperature and grain size have dominating effect on these properties.

Atomic scale modelling of the cores of dislocations in complex materials part 2: applications

In an accompanying article, we have described a methodology for the simulation of dislocations in structurally complex materials. We illustrate the applicability of this method through studies of screw dislocations in a structurally simple ionic ceramic (MgO), a molecular ionic mineral (forsterite, Mg2SiO4), a semi-ionic zeolite (siliceous zeolite A) and a covalent molecular crystalline material (the pharmaceutical, orthorhombic paracetamol-II). We focus on the extent of relaxation and the structure of the dislocation cores and comment on similarities and points of disparity between these materials. It is found that the magnitude of the relaxation varies from material to material and does not simply correlate with the magnitude of the principal elastic constants in an easily predictable fashion, or with the size of the cohesive lattice energy or length of the Burgers vector, which emphasises the need to model the non-linear forces and atomic structure of the core.