skip to main content

Title: Three-dimensional seismic anisotropy in the Pacific upper mantle from inversion of a surface-wave dispersion data set

We present a new, 3-D model of seismic velocity and anisotropy in the Pacific upper mantle, PAC13E. We invert a data set of single-station surface-wave phase-anomaly measurements sensitive only to Pacific structure for the full set of 13 anisotropic parameters that describe surface-wave anisotropy. Realistic scaling relationships for surface-wave azimuthal anisotropy are calculated from petrological information about the oceanic upper mantle and are used to help constrain the model. The strong age dependence in the oceanic velocities associated with plate cooling is also used as a priori information to constrain the model. We find strong radial anisotropy with vSH > vSV in the upper mantle; the signal peaks at depths of 100–160 km. We observe an age dependence in the depth of peak anisotropy and the thickness of the anisotropic layer, which both increase with seafloor age, but see little age dependence in the depth to the top of the radially anisotropic layer. We also find strong azimuthal anisotropy, which typically peaks in the asthenosphere. The azimuthal anisotropy at asthenospheric depths aligns better with absolute-plate-motion directions while the anisotropy within the lithosphere aligns better with palaeospreading directions. The relative strengths of radial and azimuthal anisotropy are consistent with A-type olivine more » fabric. Our findings are generally consistent with an explanation in which corner flow at the ridge leads to the development and freezing-in of anisotropy in the lithosphere, and shear between the lithosphere and underlying asthenosphere leads to anisotropy beneath the plate. We also observe large regions within the Pacific basin where the orientation of anisotropy and the absolute-plate-motion direction differ; this disagreement suggests the presence of shear in the asthenosphere that is not aligned with absolute-plate-motion directions. Azimuthal-anisotropy orientation rotates with depth; the depth of the maximum vertical gradient in the fast-axis orientation tends to be age dependent and agrees well with a thermally controlled lithosphere–asthenosphere boundary. We observe that azimuthal-anisotropy strength at shallow depths depends on half-spreading rate, with higher spreading rates associated with stronger anisotropy. Our model implies that corner flow is more efficient at aligning olivine to form lattice-preferred orientation anisotropy fabrics in the asthenosphere when the spreading rate at the ridge is higher.

« less
; ;
Publication Date:
Journal Name:
Geophysical Journal International
Page Range or eLocation-ID:
p. 355-383
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Lithospheric seismic anisotropy illuminates mid‐ocean ridge dynamics and the thermal evolution of oceanic plates. We utilize short‐period (5–7.5 s) ambient‐noise surface waves and 15‐ to 150‐s Rayleigh waves measured across the NoMelt ocean‐bottom array to invert for the complete radial and azimuthal anisotropy in the upper ∼35 km of ∼70‐Ma Pacific lithospheric mantle, and azimuthal anisotropy through the underlying asthenosphere. Strong azimuthal variations in Rayleigh‐ and Love‐wave velocity are observed, including the first clearly measured Love‐wave 2θand 4θvariations. Inversion of averaged dispersion requires radial anisotropy in the shallow mantle (2‐3%) and the lower crust (4‐5%), with horizontal velocities (VSH) faster than vertical velocities (VSV). Azimuthal anisotropy is strong in the mantle, with 4.5–6% 2θvariation inVSVwith fast propagation parallel to the fossil‐spreading direction (FSD), and 2–2.5% 4θvariation inVSHwith a fast direction 45° from FSD. The relative behavior of 2θ, 4θ, and radial anisotropy in the mantle are consistent with ophiolite petrofabrics, linking outcrop and surface‐wave length scales.VSVremains fast parallel to FSD to ∼80 km depth where the direction changes, suggesting spreading‐dominated deformation at the ridge. The transition at ∼80 km perhaps marks the dehydration boundary and base of the lithosphere. Azimuthal anisotropy strength increases from the Moho to ∼30 km depth, consistent with flow modelsmore »of passive upwelling at the ridge. Strong azimuthal anisotropy suggests extremely coherent olivine fabric. Weaker radial anisotropy implies slightly nonhorizontal fabric or the presence of alternative (so‐called E‐type) peridotite fabric. Presence of radial anisotropy in the crust suggests subhorizontal layering and/or shearing during crustal accretion.

    « less
  2. Abstract

    The breakup of supercontinent Pangea occurred ∼200 Ma forming the Eastern North American Margin (ENAM). Yet, the precise timing and mechanics of breakup and onset of seafloor spreading remain poorly constrained. We investigate the relict lithosphere offshore eastern North America using ambient‐noise Rayleigh‐wave phase velocity (12–32 s) and azimuthal anisotropy (17–32 s) at the ENAM Community Seismic Experiment (CSE). Incorporating previous constraints on crustal structure, we construct a shear velocity model for the crust and upper ∼60 km of the mantle beneath the ENAM‐CSE. A low‐velocity lid (VSof 4.4–4.55 km/s) is revealed in the upper 15–20 km of the mantle that extends ∼200 km from the margin, terminating at the Blake Spur Magnetic Anomaly (BSMA). East of the BSMA, velocities are fast (>4.6 km/s) and characteristic of typical oceanic mantle lithosphere. We interpret the low‐velocity lid as stretched continental mantle lithosphere embedded with up to ∼15% retained gabbro. This implies that the BSMA marks successful breakup and onset of seafloor spreading ∼170 Ma, consistent with ENAM‐CSE active‐source studies that argue for breakup ∼25 Myr later than previously thought. We observe margin‐parallel Rayleigh‐wave azimuthal anisotropy (2%–4% peak‐to‐peak) in the lithosphere that approximately correlates with absolute plate motion (APM) at the time of spreading. We hypothesize that lithosphere formed duringmore »ultra‐slow seafloor spreading records APM‐modified olivine fabric rather than spreading‐parallel fabric typical of higher spreading rates. This work highlights the importance of present‐day passive margins for improving understanding of the fundamental rift‐to‐drift transition.

    « less
  3. Abstract

    The mantle transition zone (MTZ) of Earth is demarcated by solid‐to‐solid phase changes of the mineral olivine that produce seismic discontinuities at 410 and 660‐km depths. Mineral physics experiments predict that wadsleyite can have strong single‐crystal anisotropy at the pressure and temperature conditions of the MTZ. Thus, significant seismic anisotropy is possible in the upper MTZ where lattice‐preferred orientation of wadsleyite is produced by mantle flow. Here, we use a body wave method, SS precursors, to study the topography change and seismic anisotropy near the MTZ discontinuities. We stack the data to explore the azimuthal dependence of travel‐times and amplitudes of SS precursors and constrain the azimuthal anisotropy in the MTZ. Beneath the central Pacific, we find evidence for ~4% anisotropy with a SE fast direction in the upper mantle and no significant anisotropy in the MTZ. In subduction zones, we observe ~4% anisotropy with a trench‐parallel fast direction in the upper mantle and ~3% anisotropy with a trench‐perpendicular fast direction in the MTZ. The transition of fast directions indicates that the lattice‐preferred orientation of wadsleyite induced by MTZ flow is organized separately from the flow in the upper mantle. Global azimuthal stacking reveals ~1% azimuthal anisotropy in themore »upper mantle but negligible anisotropy (<1%) in the MTZ. Finally, we correct for the upper mantle and MTZ anisotropy structures to obtain a new MTZ topography model. The anisotropy correction produces±3 km difference and therefore has minor overall effects on global MTZ topography.

    « less
  4. Abstract

    The relative motion of the lithosphere with respect to the asthenosphere implies the existence of a boundary zone that accommodates shear between the rigid plates and flowing mantle. This shear zone is typically referred to as the lithosphere‐asthenosphere boundary (LAB). The width of this zone and the mechanisms accommodating shear across it have important implications for coupling between mantle convection and surface plate motion. Seismic observations have provided evidence for several physical mechanisms that might help enable relative plate motion, but how these mechanisms each contribute to the overall accommodation of shear remains unclear. Here we present receiver function constraints on the discontinuity structure of the oceanic upper mantle at the NoMelt site in the central Pacific, where local constraints on shear velocity, anisotropy, conductivity, and attenuation down to ∼300 km depth provide a comprehensive picture of upper mantle structure. We image a seismic discontinuity with a Vsv decrease of 4.5% or more over a 0–20 km thick gradient layer centered at a depth of ∼65 km. We associate this feature with the Gutenberg discontinuity (G), and interpret our observation of G as resulting from strain localization across a dehydration boundary based on the good agreement between the discontinuity depth and thatmore »of the dry solidus. Transitions in Vsv, azimuthal anisotropy, conductivity, and attenuation observed at roughly similar depths suggest that the G discontinuity represents a region of localized strain within a broader zone accommodating shear between the lithosphere and asthenosphere.

    « less
  5. Abstract

    We use surface wave measurements to reveal anisotropy as a function of depth within the Juan de Fuca and Gorda plate system. Using a two‐plane wave method, we measure phase velocity and azimuthal anisotropy of fundamental mode Rayleigh waves, solving for anisotropic shear velocity. These surface wave measurements are jointly inverted with constraints fromSKSsplitting studies using a Markov chain approach. We show that the two data sets are consistent and present inversions that offer new constraints on the vertical distribution of strain beneath the plates and the processes at spreading centers. Anisotropy of the Juan de Fuca plate interior is strongest (~2.4%) in the low‐velocity zone between ~40‐ to 90‐km depth, with ENE direction driven by relative shear between plate motion and mantle return flow from the Cascadia subduction zone. In disagreement withPnmeasurements, weak (~1.1%) lithospheric anisotropy in Juan de Fuca is highly oblique to the expected ridge‐perpendicular direction, perhaps connoting complex intralithospheric fabrics associated with melt or off‐axis downwelling. In the Gorda microplate, strong shallow anisotropy (~1.9%) is consistent withPninversions and aligned with spreading and may be enhanced by edge‐driven internal strain. Weak anisotropy with ambiguous orientation in the low‐velocity zone can be explained by Gorda's youth andmore »modest motion relative to the Pacific. Deeper (≥90 km) fabric appears controlled by regional flow fields modulated by the Farallon slab edge: anisotropy is strong (~1.8%) beneath Gorda, but absent beneath the Juan de Fuca, which is in the strain shadow of the slab.

    « less