skip to main content


Title: A systematic investigation of piercing-point-dependent seismic azimuthal anisotropy
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
Award ID(s):
1919789
NSF-PAR ID:
10308578
Author(s) / Creator(s):
 ;  ;  ;  ;  
Date Published:
Journal Name:
Geophysical Journal International
Volume:
227
Issue:
3
ISSN:
0956-540X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    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Φ: 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
  2. Abstract

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

    Seismic azimuthal anisotropy beneath Australia is investigated using splitting of the teleseismic PKS, SKKS, and SKS phases to delineate asthenospheric flow and lithospheric deformation beneath one of the oldest and fast‐moving continents on Earth. In total 511 pairs of high‐quality splitting parameters were observed at 116 seismic stations. Unlike other stable continental areas in Africa, East Asia, and North America, where spatially consistent splitting parameters dominate, the fast orientations and splitting times observed in Australia show a complex pattern, with a slightly smaller than normal average splitting time of 0.85 ± 0.33 s. On the North Australian Craton, the fast orientations are mostly N‐S, which is parallel to the absolute plate motion (APM) direction in the hotspot frame. Those observed in the South Australian Craton are mostly NE‐SW and E‐W, which are perpendicular to the maximum lithospheric horizontal shortening direction. In east Australia, the observed azimuthal anisotropy can be attributed to either APM induced simple shear or lithospheric fabric parallel to the strike of the orogenic belts. The observed spatial variations of the seismic azimuthal anisotropy, when combined with results from depth estimation utilizing the spatial coherency of the splitting parameters and seismic tomography studies, suggest that the azimuthal anisotropy in Australia can mostly be related to simple shear in the rheologically transition layer between the lithosphere and asthenosphere. Non‐APM parallel anisotropy is attributable to modulations of the mantle flow system by undulations of the bottom of the lithosphere, with a spatially variable degree of contribution from lithospheric fabric.

     
    more » « less
  4. Many researchers have used the birefringence of P‑to‑S converted waves from the Moho discontinuity to constrain the anisotropy of Earth’s crust. However, this practice ignores the substantial influence that anisotropy has on the initial amplitude of the converted wave, which adds to the splitting acquired during its propagation from Moho to the seismometer. We find that large variations in Ps birefringence estimates with back-azimuth occur theoretically in the presence of P‑wave anisotropy, which normally accompanies S‑wave anisotropy. The variations are largest for crustal anisotropy with a tilted axis of symmetry, a geometry that is often neglected in birefringence interpretations, but is commonly found in Earth’s crust. We simulated globally-distributed P‑coda datasets for 36 distinct 4‑layer crustal models with combinations of elliptical shear anisotropy or compressional anisotropy, and also incorporated the higher-order anisotropic Backus parameter C. We tested both horizontal and tilted symmetry-axis geometries and tested the birefringence tradeoff associated with Ps converted phases at the top and bottom of a thin high‑ or low‑velocity basal layer. We computed composite receiver functions (RFs) with harmonic regression over back azimuth, using multipletaper correlation with moveout corrections for the epicentral distances of 471 events, to simulate a realistic data set. We estimate Ps birefringence from the radial and transverse RFs, a strategy that is similar to previous studies. We find that Ps splitting can be a useful indicator of bulk crustal anisotropy only under restricted circumstance, either in media with no compressional anisotropy, or if the symmetry axis is horizontal throughout. In other, more-realistic cases, the inferred fast polarization of Ps birefringence estimated from synthetic RFs tends either to drift with back-azimuth, form weak penalty-function minima, or return splitting times that depend on the thickness of an anisotropic layer, rather than the birefringence accumulated within it.  

     
    more » « less
  5. Abstract

    This study represents the first campaign‐style teleseismic shear wave splitting (SWS) investigation of central Myanmar, an area that is tectonically controlled by the oblique subduction of the Indian Plate underneath the Eurasian Plate. The resulting 678 well‐defined and 247 null SWS measurements obtained from recently deployed 71 broadband seismic stations show that the Indo‐Burma Ranges (IBR) possess mostly N‐S fast orientations that are parallel to the trend of the depth contours of the subducted slab. Relative to the global average of 1.0 s, extremely large splitting times with station‐averaged values ranging from 1.28 to 2.79 s and an area‐averaged value of 2.09 ± 0.55 s are observed in the IBR. In contrast, the Central Basin (CB) and the Shan Plateau (SP) are characterized by slightly larger than normal splitting times. The fast orientations observed in the CB are mostly NE‐SW in the northern part of the study area, N‐S in the central part, and NW‐SE in the southern part. The fast orientations change from nearly N‐S along the N‐S oriented Sagaing Fault, to NW‐SE in the central and eastern portions of the SP. These observations, together with SWS measurements using local S events, crustal anisotropy measurements using P‐to‐S receiver functions, and the estimated depth of the source of anisotropy using the spatial coherency of the splitting parameters, suggest the presence of a trench‐parallel sub‐slab flow system driven by slab rollback, a trench‐perpendicular corner flow, and a trench‐parallel flow possibly entering the mantle wedge through a slab window or gap.

     
    more » « less