skip to main content


Title: Poloidal‐ and Toroidal‐Mode Mantle Flows Underneath the Cascadia Subduction Zone
Abstract

Several hypotheses have been proposed to explain intriguing circular shear wave splitting patterns in the Pacific Northwest, invoking either 2‐D entrained flows or 3‐D return flows. Here, we present some hitherto unidentified, depth‐dependent anisotropic signatures to reconcile different conceptual models. At depths shallower than 200 km, the fast propagation directions of seismic waves to the west of the Rocky Mountain are aligned sub‐parallel to the subduction direction of the Juan de Fuca and Gorda Plates. This pattern is consistent with previous onshore/offshore shear wave splitting measurements and indicates that 2‐D entrained flows dominate at shallower depths. From 300 to 500 km, two large‐scale return flows are revealed, one circulating around Nevada and Colorado and the other running around the edge of the descending Juan de Fuca slab. These observations suggest the development of toroidal‐mode mantle flows, driven by the fast rollback of the narrow, fragmented Juan de Fuca and Gorda slabs.

 
more » « less
Award ID(s):
1924282
NSF-PAR ID:
10386844
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
47
Issue:
14
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 and 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.

     
    more » « less
  2. 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
  3. Abstract

    A new azimuthal anisotropy model for the North American and Caribbean Plates, namely,, is constructed based on full waveform inversion and records from the USArray and other temporary/permanent networks deployed in the study region. A total of 180 earthquakes and 4,516 seismographic stations are employed in the inversion to simultaneously constrain radially and azimuthally anisotropic model parameters:,,, and, within the crust and mantle. Thirty‐two preconditioned conjugate gradient iterations have been utilized to minimize frequency‐dependent phase discrepancies between observed and predicted seismograms for three‐component short‐period (15–40 s) body waves and long‐period (25–100 s) surface waves. Modelexhibits complicated variations in anisotropic fabrics underneath the western and eastern United States, especially at depths shallower than 100 km. For instance, the fast axis orientations in modelsuggest the presence of trench‐perpendicular mantle flows underneath the Cascadia Subduction Zone and also follow the strikes of the Snake River Plain, the Ouachita Orogenic Front, and the Grenville and Appalachian Orogenic Belts. The amplitudes of azimuthal anisotropy reduce to around 1% at depths greater than 200 km, and the orientations are subparallel to the global plate motion directions to the east of the Rocky Mountain, except for large discrepancies in central and eastern Canada. At a depth of 700 km, the fast axes change along the trajectory of the Farallon slab underneath the Great Lakes region and Gulf of Mexico, which might indicate the development of 2‐D poloidal‐mode mantle flows perpendicular to the strike of the sinking slab within the uppermost lower mantle. Comparisons between modelwith a western U.S. model from ambient noise tomography and SKS splitting measurements demonstrate a relatively good agreement for the fast axis orientations, considering the usage of different data sets and imaging techniques. However, the absolute magnitude of azimuthal anisotropy in modelmight be underestimated, especially at greater depths, given the poor agreement on the amplitudes of predicted and observed SKS splitting times. At the current stage, the agreement among different azimuthal anisotropy models at global and continental scales is still poor even for the United States with a dense station coverage.

     
    more » « less
  4. Abstract

    Laboratory experiments and geodynamic simulations demonstrate that poloidal- and toroidal-mode mantle flows develop around subduction zones. Here, we use a new 3-D azimuthal anisotropy model constructed by full waveform inversion, to infer deep subduction-induced mantle flows underneath Middle America. At depths shallower than 150 km, poloidal-mode flow is perpendicular to the trajectory of the Middle American Trench. From 300 to 450 km depth, return flows surround the edges of the Rivera and Atlantic slabs, while escape flows are inferred through slab windows beneath Panama and central Mexico. Furthermore, at 700 km depth, the study region is dominated by the Farallon anomaly, with fast axes perpendicular to its strike, suggesting the development of lattice-preferred orientations by substantial stress. These observations provide depth-dependent seismic anisotropy for future mantle flow simulations, and call for further investigations about the deformation mechanisms and elasticity of minerals in the transition zone and uppermost lower mantle.

     
    more » « less
  5. SUMMARY

    The Earth's mantle transition zone (MTZ) plays a key role in the thermal and compositional interactions between the upper and lower mantle. Seismic anisotropy provides useful information about mantle deformation and dynamics across the MTZ. However, seismic anisotropy in the MTZ is difficult to constrain from surface wave or shear wave splitting measurements. Here, we investigate the sensitivity to anisotropy of a body wave method, SS precursors, through 3-D synthetic modelling and apply it to real data. Our study shows that the SS precursors can distinguish the anisotropy originating from three depths: shallow upper mantle (80–220 km), deep upper mantle above 410 km, and MTZ (410–660 km). Synthetic resolution tests indicate that SS precursors can resolve $\ge $3 per cent azimuthal anisotropy where data have an average signal-to-noise ratio (SNR = 7) and sufficient azimuthal coverage. To investigate regional sensitivity, we apply the stacking and inversion methods to two densely sampled areas: the Japan subduction zone and a central Pacific region around the Hawaiian hotspot. We find evidence for significant VS anisotropy (15.3 ± 9.2 per cent) with a trench-perpendicular fast direction (93° ± 5°) in the MTZ near the Japan subduction zone. We attribute the azimuthal anisotropy to the grain-scale shape-preferred orientation of basaltic materials induced by the shear deformation within the subducting slab beneath NE China. In the central Pacific study region, there is a non-detection of MTZ anisotropy, although modelling suggests the data coverage should allow us to resolve at least 3 per cent anisotropy. Therefore, the Hawaiian mantle plume has not produced detectable azimuthal anisotropy in the MTZ.

     
    more » « less