skip to main content


Title: Constraints on the Depth, Thickness, and Strength of the G Discontinuity in the Central Pacific From S Receiver Functions
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 that 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.

 
more » « less
PAR ID:
10359825
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
126
Issue:
4
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. SUMMARY

    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 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.

     
    more » « less
  2. Abstract

    This study probes the lithosphere‐asthenosphere system beneath 155 Ma Pacific seafloor using teleseismic S‐to‐p receiver functions at the Pacific Lithosphere Anisotropy and Thickness Experiment project ocean‐bottom‐seismometers. Within the lithosphere, a significant velocity decrease at 33–50 km depth is observed. This mid‐lithospheric discontinuity is consistent with the velocity contrast between the background mantle and thin, trapped layers of crystallized partial melt, in the form of either dolomite or garnet granulite. These melts possibly originated from deeper asthenospheric melting beneath the flanks of spreading centers, and were transported within the cooling lithosphere. A positive velocity increase of 3%–6% is observed at 130–155 km depth and is consistent with the base of a layer with partial melt in the asthenosphere. A shear velocity decrease associated with the lithosphere‐asthenosphere boundary at 95–115 km depth is permitted by the data, but is not required.

     
    more » « less
  3. The asthenosphere plays a fundamental role in present-day plate tectonics as its low viscosity controls how convection in the mantle below it is expressed at the Earth’s surface above. The origin of the asthenosphere, including the role of partial melting in reducing its viscosity and facilitating deformation, remains unclear. Here we analysed receiver-function data from globally distributed seismic stations to image the lower reaches of the asthenospheric low-seismic-velocity zone. We present globally widespread evidence for a positive seismic-velocity gradient at depths of ~150 km, which represents the base of a particularly low-velocity zone within the asthenosphere. This boundary is most commonly detected in regions with elevated upper-mantle temperatures and is best modelled as the base of a partially molten layer. The presence of the boundary showed no correlation with radial seismic anisotropy, which represents accumulated mantle strain, indicating that the inferred partial melt has no substantial effect on the large-scale viscosity of the asthenosphere. These results imply the presence of a globally extensive, partially molten zone embedded within the asthenosphere, but that low asthenospheric viscosity is controlled primarily by gradual pressure and temperature variations with depth. 
    more » « less
  4. Abstract

    Seismic anisotropy beneath eastern North America, as expressed in shear wave splitting observations, has been attributed to plate motion‐parallel shear in the asthenosphere, resulting in fast axes aligned with the plate motion. However, deviations of fast axes from plate motion directions are observed near major tectonic boundaries of the Appalachians, indicating contributions from lithospheric anisotropy associated with past tectonic processes. In this study, we conduct anisotropic receiver function (RF) analysis using data from a dense seismic array traversing the New England Appalachians in Connecticut to examine anisotropic layers in the crust and upper mantle and correlate them with past tectonic processes as well as present‐day mantle flow. We use the harmonic decomposition method to separate directionally‐dependent variations of RFs and focus on features with the same harmonic signals observed across multiple stations. Within the crust, there are multiple features that may be correlated with stratification in the Hartford Basin, faults in the Taconic thrust belt, shear zones formed during Salinic/Acadian terrane accretion events, and orogen‐parallel crustal flow in the Acadian orogenic plateau. We apply a Bayesian inversion method to obtain quantitative constraints on the direction and strength of intra‐crustal anisotropy beneath the Hartford Basin. In the upper mantle, we identify a fossil shear zone possibly formed during oblique subduction of Rheic Ocean lithosphere. We also find evidence for a plate motion‐parallel flow zone in the asthenosphere that is likely disturbed by mantle upwelling near the southern margin of the Northern Appalachian Anomaly in the eastern part of the study area.

     
    more » « less
  5. Seismic waves with different propagation and oscillation directions can exhibit different velocities when going through a medium with some directional properties; this phenomenon is called seismic anisotropy. Seismic anisotropy observed beneath eastern North America is often attributed to present-day flow in the upper mantle. The mantle flow causes shear waves oscillating in the direction of flow (e.g., in the direction of North America plate motion) to travel faster than those that travel in other directions. However, this pattern does not hold true for some regions along the Appalachian orogen, suggesting that past tectonic events can result in long-lived, ‘frozen-in’ anisotropy in the lithosphere, which modifies the predicted anisotropic behavior beneath these regions. In this study, we investigate sources of seismic anisotropy beneath southern New England using a method based on directionally dependent variations of P-wave to S-wave conversions at interfaces with contrasts in anisotropy. This method can separate signals caused by different anisotropic features and constrain the depth distribution of anisotropy. Within the crust there are multiple features that may be correlated with stratification in the Hartford Basin, faults in the Taconic thrust belt, shear zones formed during Salinic/Acadian terrane accretion events, and orogen-parallel crustal flow in the Acadian orogenic plateau. We apply a Bayesian inversion method to obtain quantitative constraints on the direction and strength of intra-crustal anisotropy beneath the Hartford Basin. In the upper mantle, we identify a fossil shear zone possibly formed during oblique subduction of Rheic Ocean lithosphere. We also find evidence for a plate motion-parallel flow zone in the asthenosphere that is likely disturbed by mantle upwelling near the southern margin of the Northern Appalachian Anomaly in the eastern part of the study area. 
    more » « less