Observations of seismic anisotropy can provide direct constraints on the character of mantle flow in subduction zones, critical for our broader understanding of subduction dynamics. Here we present over 750 new SKS splitting measurements in the vicinity of Mount St. Helens in the Cascadia subduction zone using a combination of stations from the iMUSH broadband array and Cascades Volcano Observatory network. This provides the highest density of splitting measurements yet available in Cascadia, acting as a focused “telescope” for seismic anisotropy in the subduction zone. We retrieve spatially consistent splitting parameters (mean fast direction
Seismic anisotropy provides insight into past episodes of lithospheric deformation and the orientations of strain in the underlying asthenosphere. The Greenland mantle has played host to a rich history of tectonic processes, including multiple orogenies and plume‐lithosphere interactions. This study presents new measurements of SKS splitting that reveal strong variations in fast polarization direction with back‐azimuth that are consistent across Greenland, including at stations where splitting measurements have not previously been reported. We compared observed fast polarization directions to the predictions of two‐layer models with olivine‐orthopyroxene anisotropy. The family of models that provides acceptable misfits at 95% confidence indicates an upper layer olivine a‐axis azimuth of 222–236°, a lower layer olivine a‐axis azimuth of 114–130°, and non‐zero a‐axis plunges. These models are consistent with an upper layer of lithospheric anisotropy due to Proterozoic and Archean orogenic fabrics, and a lower layer of anisotropy corresponding to either asthenospheric flow aligned approximately parallel to the direction of absolute plate motion and plunging due to lithospheric topography or dipping lithospheric structures created by episodes of paleo‐convergence.
more » « less- NSF-PAR ID:
- 10450750
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geochemistry, Geophysics, Geosystems
- Volume:
- 22
- Issue:
- 5
- ISSN:
- 1525-2027
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract Φ : 74°, mean delay time∂t : 1.0 s) with the azimuthal occurrence of nulls in agreement with the fast direction of splitting. When averaged across the array, a 90° periodicity in splitting parameters as a function of back azimuth is revealed, which has not been recovered previously with single‐station observations. The periodicity is characterized by a sawtooth pattern inΦ with a clearly defined 45° trend. We present new equations that reproduce this behavior based upon known systematic errors when calculating shear wave splitting from data with realistic seismic noise. The corrected results suggest a single layer of anisotropy with an ENE‐WSW fast axis parallel to the motion of the subducting Juan de Fuca plate; in agreement with predictions for entrained subslab mantle flow. The splitting pattern is consistent with that seen throughout Cascadia, suggesting that entrainment of the underlying asthenosphere with the subducting slab is coherent and widespread. -
more » « less
Abstract Shear‐wave splitting observations of SKS and SKKS phases have been used widely to map azimuthal anisotropy, as caused by the occurrence of olivine, to constrain the dominant directions of upper mantle deformation. While SK(K)S splitting measurements at individual seismic stations are often averaged before interpretation, it is useful to consider additional information, for example, based on the variation of splitting parameters with azimuth due to the non‐vertical incidence of core‐phases. These constraints in theory enable a differentiation between various types of olivine and may allow us to infer otherwise poorly known upper mantle parameters such as stress, temperature, and water content. In this study, we predict the azimuthal variation of splitting parameters for A‐, C‐, and E‐type olivine fabrics and match them with observations from the High Lava Plains, Northwestern Basin and Range, and Western Yellowstone Snake River Plain in the Pacific Northwest US. This helps to constrain the amount of water in the upper mantle in the back‐arc of the Cascadia subduction zone, known for its consistent E‐W oriented seismic anisotropy, and particularly large splitting delay times. Our investigation renders the C‐type olivine mechanism improbable for this location; A‐ and E‐type fabrics match the observations, although differentiating between them is difficult. However, the agreement of the amplitude of backazimuthal variation of the fast orientation, plus the potential to explain large splitting delay times, suggest the occurrence of E‐type olivine, and thus the likely presence of a moderately hydrated upper mantle beneath Cascadia's back‐arc.
-
SUMMARY The vast majority of teleseismic XKS (including SKS, SKKS and PKS) shear wave splitting studies interpret the observed splitting parameters (fast orientation and splitting time) based on the assumption of a spatially invariant anisotropy structure in the vicinity of a recording station. For such anisotropy structures the observed splitting parameters are either independent of the arriving azimuth of the seismic ray paths if the medium traversed by the ray paths can be represented by a single layer of anisotropy with a horizontal axis of symmetry (i.e. simple anisotropy), or demonstrate a periodic variation with respect to the arriving azimuth for a more complicated structure of anisotropy (e.g. multiple layers with a horizontal axis of symmetry, or a single layer with a dipping axis). When a recording station is located near the boundary of two or more regions with different anisotropy characteristics, the observed splitting parameters are dependent on the location of the ray piercing points. Such a piercing-point dependence is clearly observed using a total of 360 pairs of XKS splitting parameters at three stations situated near the northeastern edge of the Sichuan Basin in central China. For a given station, the fast orientations differ as much as 90°, and the azimuthal variation of the fast orientations lacks a 90° or 180° periodicity which is expected for double-layered or dipping axis anisotropy. The observed splitting parameters from the three stations are spatially most consistent when they are projected at a depth of ∼250 km, and can be explained by shear strain associated with the absolute plate motion and mantle flow deflected by the cone-shaped lithospheric root of the Sichuan Basin.more » « less
-
more » « less
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 (V S H ) faster than vertical velocities (V S V ). Azimuthal anisotropy is strong in the mantle, with 4.5–6% 2θ variation inV S V with fast propagation parallel to the fossil‐spreading direction (FSD), and 2–2.5% 4θ variation inV S H with 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.V S V remains 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 models 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. -
more » « less
SUMMARY Determinations of seismic anisotropy, or the dependence of seismic wave velocities on the polarization or propagation direction of the wave, can allow for inferences on the style of deformation and the patterns of flow in the Earth’s interior. While it is relatively straightforward to resolve seismic anisotropy in the uppermost mantle directly beneath a seismic station, measurements of deep mantle anisotropy are more challenging. This is due in large part to the fact that measurements of anisotropy in the deep mantle are typically blurred by the potential influence of upper mantle and/or crustal anisotropy beneath a seismic station. Several shear wave splitting techniques are commonly used that attempt resolve seismic anisotropy in deep mantle by considering the presence of multiple anisotropic layers along a raypath. Examples include source-side S-wave splitting, which is used to characterize anisotropy in the deep upper mantle and mantle transition zone beneath subduction zones, and differential S-ScS and differential SKS-SKKS splitting, which are used to study anisotropy in the D″ layer at the base of the mantle. Each of these methods has a series of assumptions built into them that allow for the consideration of multiple regions of anisotropy. In this work, we systematically assess the accuracy of these assumptions. To do this, we conduct global wavefield modelling using the spectral element solver AxiSEM3D. We compute synthetic seismograms for earth models that include seismic anisotropy at the periods relevant for shear wave splitting measurements (down to 5 s). We apply shear wave splitting algorithms to our synthetic seismograms and analyse whether the assumptions that underpin common measurement techniques are adequate, and whether these techniques can correctly resolve the anisotropy incorporated in our models. Our simulations reveal some inaccuracies and limitations of reliability in various methods. Specifically, explicit corrections for upper mantle anisotropy, which are often used in source-side direct S splitting and S-ScS differential splitting, are typically reliable for the fast polarization direction ϕ but not always for the time lag δt, and their accuracy depends on the details of the upper mantle elastic tensor. We find that several of the assumptions that underpin the S-ScS differential splitting technique are inaccurate under certain conditions, and we suggest modifications to traditional S-ScS differential splitting approaches that lead to improved reliability. We investigate the reliability of differential SKS-SKKS splitting intensity measurements as an indicator for lowermost mantle anisotropy and find that the assumptions built into the splitting intensity formula can break down for strong splitting cases. We suggest some guidelines to ensure the accuracy of SKS-SKKS splitting intensity comparisons that are often used to infer lowermost mantle anisotropy. Finally, we suggest a new strategy to detect lowermost mantle anisotropy which does not rely on explicit upper mantle corrections and use this method to analyse the lowermost mantle beneath east Asia.
