skip to main content


Title: Anisotropy-revealed change in hydration along the Alaska subduction zone
Abstract Megathrust earthquake behavior in subduction zones is controlled by a variety of factors including the hydration state of the subducting slab. Increased hydration reduces the occur-rence of great, damaging earthquakes by diminishing the strength of the material along the interface between tectonic plates. Understanding variations in hydration in subductions zones is necessary for properly assessing the overall hazard posed by each region. Fortunately, seismic anisotropy is strongly dependent upon hydration of the subducting crust and litho-sphere. I present shear-wave splitting measurements that illuminate changes in anisotropy, and therefore hydration, of the subducting Pacific plate beneath the Alaska subduction zone (northern Pacific Ocean). Variations in shear-wave splitting directly correlate to changes in the behavior of great, megathrust earthquakes. My measurements show that the Shumagin seismic gap is characterized by a hydrated subducting slab, explaining the long-term lack of great earthquakes. Observations in the immediately adjacent Semidi segment, which experiences great events regularly, indicate a far less hydrated slab. These results are driven by the preferential alignment of paleo-spreading fabrics of the Pacific plate. Where fabrics are more closely aligned with the orientation of the trench, outer-rise faulting and plate hydration is enhanced. These results highlight the importance of changes in preexisting slab structures and subsequent hydration in the production of great, damaging earthquakes.  more » « less
Award ID(s):
1949210
NSF-PAR ID:
10293874
Author(s) / Creator(s):
Date Published:
Journal Name:
Geology
Volume:
49
Issue:
9
ISSN:
0091-7613
Page Range / eLocation ID:
1122 to 1125
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Shear‐wave splitting observations can provide insight into mantle flow, due to the link between the deformation of mantle rocks and their direction‐dependent seismic wave velocities. We identify anisotropy in the Cook Inlet segment of the Alaska subduction zone by analyzing splitting parameters of S waves from local intraslab earthquakes between 50 and 200 km depths, recorded from 2015–2017 and emphasizing stations from the Southern Alaska Lithosphere and Mantle Observation Network experiment. We classify 678 high‐quality local shear‐wave splitting observations into four regions, from northwest to southeast: (L1b) splitting measurements parallel to Pacific plate motion, (L1a) arc‐perpendicular splitting pattern, (L2) sharp transition to arc‐parallel splitting, and (L3) splitting parallel to Pacific plate motion. Forward modeling of splitting from various mantle fabrics shows that no one simple model fully explains the observed splitting patterns. An A‐type olivine fabric with fast direction dipping 45° to the northwest (300°)—aligned with the dipping slab—predicts fast directions that fit L1a observations well, but not L2. The inability of the forward model fabrics to fit all the observed splitting patterns suggests that the anisotropy variations are not due to variable ray angles, but require distinct differences in the anisotropy regime below the arc, forearc, and subducting plate.

     
    more » « less
  2. Subduction of the very young (<15 Myr old) oceanic lithosphere of the Nazca plate in central to southern Colombia is observationally related to an unusually high and unusually variable amount of intermediate (>50 km) depth seismicity. From 2010 through 2019 89% of central and southern Colombia’s 11,466 intermediate depth events occurred between 3.5°N and 5.5°N, highlighting these unusual characteristics of the young slab. In addition, morphologic complexity and possible tears characterize the Nazca slab in Colombia and complicate mantle flow in the region. Prior SKS-phase shear-wave splitting results indicate sub-slab anisotropy is dominated by plate motion parallel-to-subparallel orientations in the region, suggesting the young slab has entrained a relatively thick portion of the sub-slab mantle. These observations suggest the subduction of young lithosphere has significant effects on both the overlying and underlying asthenosphere in the Colombia subduction zone. Here we use more than 10 years of data to calculate receiver functions for the Red Sismológica Nacional de Colombia’s network of broadband seismometers. These receiver functions allow us to tie these prior observations of the Colombia subduction zone to distinct, structural features of the slab. We find that the region of high seismicity corresponds to a low seismic velocity feature along the top of the subducting plate between 3.5°N and 5.5°N that is not present to the south. Moderately elevated P-wave velocity to S-wave velocity ratios are also observed within the slab in the north. This feature likely represents hydrated slab mantle and/or uneclogitized oceanic crust extending to a deeper depth in the north of the region which may provide fluids to drive slab seismicity. We further find evidence for a thick layer of material along the slab’s lithosphere-asthenosphere boundary characterized by spatially variable anisotropy. This feature likely represents entrained asthenosphere at the base of the plate sheared by both the overlying plate and complex flow related to proposed slab tears just north and south of the study region. These observations highlight how structural observations provide key contextual constraints on short-term (seismogenic) and long-term (anisotropic fabric) dynamic processes in the Colombia subduction zone. Plain-language Summary The Nazca oceanic plate is very young (<15 million years old) where it is pulled or subducted beneath the South America plate in central and southern Colombia. Earthquakes occurring in the subducted Nazca plate at depths greater than 50 km are nearly 9x more common in central Colombia than in southern Colombia. The subducted Nazca plate also has a complex shape in this region and may have been torn both in northern Colombia and to the south near the Colombia-Ecuador border. The slow flow of mantle rock beneath the subducted plate is believed to be affected by this and earlier studies have inferred this flow is mostly in the same direction as the subducting plate's motion. We have used 10+ years of data to calculate receiver functions, which can detect changes in the velocity of seismic waves at the top and bottom of the subducted plate to investigate these features. We found that the Nazca plate is either hydrated or has rocks with lower seismic velocities at its top in the central part of Colombia where earthquakes are common. We also find that a thick layer of mantle rock at the base of the subducted plate has been sheared. 
    more » « less
  3. 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
  4. Abstract

    Seismic anisotropy measurements show that upper mantle hydration at the Middle America Trench (MAT) is limited to serpentinization and/or water in fault zones, rather than distributed uniformly. Subduction of hydrated oceanic lithosphere recycles water back into the deep mantle, drives arc volcanism, and affects seismicity at subduction zones. Constraining the extent of upper mantle hydration is an important part of understanding many fundamental processes on Earth. Substantially reduced seismic velocities in tomography suggest that outer rise plate‐bending faults provide a pathway for seawater to rehydrate the slab mantle just prior to subduction. Estimates of outer‐rise hydration based on tomograms vary significantly, with some large enough to imply that, globally, subduction has consumed more than two oceans worth of water during the Phanerozoic. We found that, while the mean upper mantle wavespeed is reduced at the MAT outer rise, the amplitude and orientation of inherited anisotropy are preserved at depths >1 km below the Moho. At shallower depths, relict anisotropy is replaced by slowing in the fault‐normal direction. These observations are incompatible with pervasive hydration but consistent with models of wave propagation through serpentinized fault zones that thin to <100‐m in width at depths >1 km below Moho. Confining hydration to fault zones reduces water storage estimates for the MAT upper mantle from ∼3.5 wt% to <0.9 wt% H20. Since the intermediate thermal structure in the ∼24 Myr‐old MAT slab favors serpentinization, limited hydration suggests that fault mechanics are the limiting factor, not temperatures. Subducting mantle may be similarly dry globally.

     
    more » « less
  5. SUMMARY The Ecuadorian forearc is a complex region of accreted terranes with a history of large megathrust earthquakes. Most recently, a Mw 7.8 megathrust earthquake ruptured the plate boundary offshore of Pedernales, Ecuador on 16 April 2016. Following this event, an international collaboration arranged by the Instituto Geofisico at the Escuela Politécnica Nacional mobilized a rapid deployment of 65 seismic instruments along the Ecuadorian forearc. We combine this new seismic data set with 14 permanent stations from the Ecuadorian national network to better understand how variations in crustal structure relate to regional seismic hazards along the margin. Here, we present receiver function adaptive common conversion point stacks and a shear velocity model derived from the joint inversion of receiver functions and surface wave dispersion data obtained through ambient noise cross-correlations for the upper 50 km of the forearc. Beneath the forearc crust, we observe an eastward dipping slow velocity anomaly we interpret as subducting oceanic crust, which shallows near the projected centre of the subducting Carnegie Ridge. We also observe a strong shallow positive conversion in the Ecuadorian forearc near the Borbon Basin indicating a major discontinuity at a depth of ∼7 km. This conversion is not ubiquitous and may be the top of the accreted terranes. We also observe significant north–south changes in shear wave velocity. The velocity changes indicate variations in the accreted terranes and may indicate an increased amount of hydration beneath the Manabí Basin. This change in structure also correlates geographically with the southern rupture limit of multiple high magnitude megathrust earthquakes. The earthquake record along the Ecuadorian trench shows that no event with a Mw >7.4 has ruptured south of ∼0.5°S in southern Ecuador or northern Peru. Our observations, along with previous studies, suggest that variations in the forearc crustal structure and subducting oceanic crust may influance the occurrence and spatial distribution of high magnitude seismicity in the region. 
    more » « less