skip to main content

Title: Seismic anisotropy, dominant slip systems and phase transitions in the lowermost mantle
SUMMARY The presence of seismic anisotropy at the base of the Earth's mantle is well established, but there is no consensus on the deformation mechanisms in lower mantle minerals that could explain it. Strong anisotropy in magnesium post-perovskite (pPv) has been invoked, but different studies disagree on the dominant slip systems at play. Here, we aim to further constrain this by implementing the most recent results from atomistic models and high-pressure deformation experiments, coupled with a realistic composition and a 3-D geodynamic model, to compare the resulting deformation-induced anisotropy with seismic observations of the lowermost mantle. We account for forward and reverse phase transitions from bridgmanite (Pv) to pPv. We find that pPv with either dominant (001) or (010) slip can both explain the seismically observed anisotropy in colder regions where downwellings turn to horizontal flow, but only a model with dominant (001) slip matches seismic observations at the root of hotter large-scale upwellings. Allowing for partial melt does not change these conclusions, while it significantly increases the strength of anisotropy and reduces shear and compressional velocities at the base of upwellings.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Geophysical Journal International
Page Range / eLocation ID:
1665 to 1681
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. In recent years there have been several attempts to make the link between mineral properties and seismic anisotropy in the D’’ region but have yet to reach consensus with regards to the dynamics in lower mantle minerals that could give rise to the observed seismic anisotropy. Here, we aim to provide further constraints on the observed long wavelength shear velocity patterns seen in seismic tomography studies. We introduce a forward model of deformation in a subducting slab as it impacts the core mantle boundary (D’’ layer) and proceeds to upwelling at the edge of a simulated LLSVP. By implementing the most recent results from atomistic modeling and high-pressure deformation experiments coupled with a 3-dimensional geodynamic model, we compare the microstructural evolution of an aggregate with a pyrolytic composition to the macroscopically observed seismic anisotropy of the lowermost mantle. We account for topotaxial relations in the forward and reverse phase transitions of MgSiO3-perovskite (Pv) to post-perovskite (pPv) within the slab as well as explore the effects introduced by partial melting near the CMB. Comparisons in the two leading candidate deformation mechanisms in the post-perovskite phase, (001) and (010), are compared. In this study we find that the reverse transition (pPv to Pv) occurs at a depth which is ~ 150 km deeper than that of the forward transition due to increasing temperature near the CMB providing a varying topography of the D’’ discontinuity. Our model also produces good fits with the isotropic velocities of PREM for the bulk lower mantle. When coupled with temperature and pressure dependent forward and reverse phase transitions, a pPv system with dominant (001) slip provides good correlation with the currently observed VSH fast horizontal (~ 1 – 6%) in D’’ and with VSV consistently fast in upwelling areas. Azimuthal variations along the streamline are also investigated showing a symmetry lower than that of the assumed VTI in D’’ introduced by ‘rolling’ effects near the slab’s edge. The addition of 1% partial melting at the CMB is shown to increase S and P wave anisotropy beneath the slab at the base of upwelling with up to ~2.5 & 4.0% P and S wave reductions respectively compared to the global reference. 
    more » « less
  2. Abstract

    The exact mechanism for lowermost mantle seismic anisotropy remains unknown; however, work on the elasticity and deformation of lower mantle materials has constrained a few possible options. The most probable minerals producing anisotropy are bridgmanite, postperovskite, and ferropericlase. While there is an extensive literature on the elasticity and deformation of lower mantle minerals, we create a comprehensive uniform database ofD″ anisotropy scenarios. In order to characterize a range of the possible fabrics forD″ anisotropy, we carry out VPSC (visco‐plastic self‐consistent modeling) to predict textures for each proposed mineral and dominant slip system. We numerically deform each mineral under different geometrical scenarios: simple shear, pure shear, and extension. By using published single crystal elasticity values, we produce a library of 336 candidate elastic tensors. We used the elastic tensor library to revisit previously publishedD″‐associated seismic anisotropy studies for crossing raypaths (Siberia, North America, the Afar region of Africa, and Australia). While we cannot identify a single, unique mechanism that explains all of these data sets, we find that postperovskite (dominant slip on [100](010) or [100](001)) and periclase (dominant slip on {100}<011>) provide the best fit to the observations and suggest reasonable shear directions for each region of interest. Bridgmanite generally provides a poor fit to the observations; however, we cannot completely rule out any particular model. As the number of anisotropy observations forD″ increases, this elastic tensor library will be helpful for observational seismologists in identifying possible mechanisms of anisotropy and shear directions at in the lowermost mantle.

    more » « less
  3. Abstract

    Observations of seismic anisotropy just above the core‐mantle boundary are a powerful way to understand the dynamics of the lowermost mantle. Here we present models of seismic anisotropy in the lowermost mantle beneath Siberia based on new and previously published body wave observations. We compile a set of measurements based on a variety of data types, including the differential splitting of SK(K)S‐PKS and S‐ScS phases and polarities of PdP and SdS reflections off the D″ discontinuity. We carry out ray theoretical forward modeling of these data sets in combination, using a novel approach that allows for tighter constraints on the geometry of seismic anisotropy than would be possible using a single data type. Observations for each seismic phase alone only provide limited information; by combining different data types into one modeling approach, we can constrain D″ seismic anisotropy more tightly. We test a range of plausible mechanisms for seismic anisotropy, including a variety of candidate minerals (bridgmanite, post‐perovskite, and ferropericlase), dominant slip systems, and orientations, and consider both single‐crystal elasticity and elastic tensors based on polycrystalline texture modeling. In general, we find that post‐perovskite (with either a [100](010) or [100](001) dominant slip system) or bridgmanite models provide the best fit to the observations. The best‐fitting models suggest a dominant shear direction to the southwest or northeast direction at the base of the mantle. This is consistent with flow directed toward the Perm anomaly, potentially driven by slab remnants impinging on the core‐mantle boundary.

    more » « less
  4. Abstract

    Seismic anisotropy arises in the upper mantle due to the alignment of olivine crystal lattices and is often used to interpret mantle flow direction. Experiments on the evolution of olivine crystal‐preferred orientation (CPO) have found that the texture that develops is dependent on many factors, including water content, differential stress, preexisting CPO, and deformation kinematics. To evaluate the role of these factors in naturally deformed samples, we present microstructural transects across three shear zones in the Josephine Peridotite. Samples from these shear zones exhibit a mixture of A‐type textures, which have been associated with dry conditions and primary activation of the olivine [100](010) slip system, and of E‐type textures, which have been associated with wetter conditions and primary activation of the [100](001) slip system. CPOs with characteristics of both A‐type and E‐type textures are also present. CPO type does not evolve systematically as a function of either strain or water content. We used a micromechanical model to evaluate the roles of preexisting texture and kinematics on olivine CPO evolution. We find that the preexisting texture controls CPO evolution at strains up to 5 during simple shear. Kinematics involving a combination of simple shear and pure shear can explain the olivine CPOs at higher strain. Hence, preexisting CPOs and deformation kinematics should be considered in the interpretation of CPOs measured in naturally deformed rocks and of large‐scale patterns in upper‐mantle seismic anisotropy.

    more » « less
  5. Abstract

    Seismic anisotropy produced by aligned olivine in oceanic lithosphere offers a window into mid‐ocean ridge (MOR) dynamics. Yet, interpreting anisotropy in the context of grain‐scale deformation processes and strain observed in laboratory experiments and natural olivine samples has proven challenging due to incomplete seismological constraints and length scale differences spanning orders of magnitude. To bridge this observational gap, we estimate an in situ elastic tensor for oceanic lithosphere using co‐located compressional‐ and shear‐wavespeed anisotropy observations at the NoMelt experiment located on ∼70 Ma seafloor. The elastic model for the upper 7 km of the mantle, NoMelt_SPani7, is characterized by a fast azimuth parallel to the fossil‐spreading direction, consistent with corner‐flow deformation fabric. We compare this model with a database of 123 petrofabrics from the literature to infer olivine crystallographic orientations and shear strain accumulated within the lithosphere. Direct comparison to olivine deformation experiments indicates strain accumulation of 250%–400% in the shallow mantle. We find evidence for D‐type olivine lattice‐preferred orientation (LPO) with fast [100] parallel to the shear direction and girdled [010] and [001] crystallographic axes perpendicular to shear. D‐type LPO implies similar amounts of slip on the (010)[100] and (001)[100] easy slip systems during MOR spreading; we hypothesize that grain‐boundary sliding during dislocation creep relaxes strain compatibility, allowing D‐type LPO to develop in the shallow lithosphere. Deformation dominated by dislocation‐accommodated grain‐boundary sliding (disGBS) has implications for in situ stress and grain size during MOR spreading and implies grain‐size dependent deformation, in contrast to pure dislocation creep.

    more » « less