Abstract The Alaska Peninsula has a long history of plate subduction with along‐arc variations in volcanic eruption styles and geochemistry. However, the sub‐arc melting processes that feed these volcanoes are unclear. The Alaska slab morphology below 200 km depth remains debated due to limited seismic data and thus low tomography resolution in this region. Here we utilize the newly available regional and teleseismic data to build 3‐D high‐resolutionVPandVSmodels to 660 km depth. We find that the high‐velocity Pacific Plate subducts to the bottom of the mantle transition zone (MTZ) with complex deformation and gaps. In the southwest, we observe a wide gap in the high‐velocity slab at 200–500 km depths. Toward the northeast, the slab becomes more continuous extending to the MTZ with a few holes below 200 km. We interpret these gaps as a slab tear that coincides with the subducted ancient Kula‐Pacific Ridge. We also invert for 3‐DVPandVP/VSmodels to 200 km depth with higher resolution and find strong along‐strike changes in slab dehydration and sub‐arc melting, indicated by lowVPand highVP/VSanomalies. Slab dehydration and sub‐arc melting are most extensive below the Pavlof and Shumagin segments in the southwest, becoming limited below the Chignik and Chirikof segments in the northeast, and extensive again beneath the Kodiak segment further to the northeast. We propose that the variations of slab hydration at the outer rise significantly influence slab dehydration at greater depths and further control sub‐arc melting beneath the Alaska Peninsula.
more »
« less
Along‐Strike Variations of Alaska Subduction Zone Structure and Hydration Determined From Amphibious Seismic Data
Abstract We develop a 3‐D isotropic shear velocity model for the Alaska subduction zone using data from seafloor and land‐based seismographs to investigate along‐strike variations in structure. By applying ambient noise and teleseismic Helmholtz tomography, we derive Rayleigh wave group and phase velocity dispersion maps, then invert them for shear velocity structure using a Bayesian Monte Carlo algorithm. For land‐based stations, we perform a joint inversion of receiver functions and dispersion curves. The forearc crust is relatively thick (35–42 km) and has reduced lower crustal velocities beneath the Kodiak and Semidi segments, which may promote higher seismic coupling. Bristol Bay Basin crust is relatively thin and has a high‐velocity lower layer, suggesting a dense mafic lower crust emplaced by the rifting processes. The incoming plate shows low uppermost mantle velocities, indicating serpentinization. This hydration is more pronounced in the Shumagin segment, with greater velocity reduction extending to 18 ± 3 km depth, compared to the Semidi segment, showing smaller reductions extending to 14 ± 3 km depth. Our estimates of percent serpentinization from VSreduction and VP/VSare larger than those determined using VPreduction in prior studies, likely due to water in cracks affecting VSmore than VP. Revised estimates of serpentinization show that more water subducts than previous studies, and that twice as much mantle water is subducted in the Shumagin segment compared to the Semidi segment. Together with estimates from other subduction zones, the results indicate a wide variation in subducted mantle water between different subduction segments.
more »
« less
- PAR ID:
- 10538864
- Publisher / Repository:
- American Geophysical Union
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 129
- Issue:
- 3
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract A shallow sub‐seafloor seismic model that includes well‐determined seismic velocities and clarifies sediment‐crust discontinuities is needed to characterize the physical properties of marine sediments and the oceanic crust and to serve as a reference for deeper seismic modeling endeavors. This study estimates the seismic structure of marine sediments and the shallow oceanic crust of the Alaska‐Aleutian subduction zone at the Alaska Peninsula, using data from the Alaska Amphibious Community Seismic Experiment (AACSE). We measure seafloor compliance and Ps converted wave delays from AACSE ocean‐bottom seismometers (OBS) and seafloor pressure data and interpret these measurements using a joint Bayesian Monte Carlo inversion to produce a sub‐seafloor S‐wave velocity model beneath each available OBS station. The sediment thickness across the array varies considerably, ranging from about 50 m to 2.80 km, with the thickest sediment located in the continental slope. Lithological composition plays an important role in shaping the seismic properties of seafloor sediment. Deep‐sea deposits on the incoming plate, which contain biogenic materials, tend to have reduced S‐wave velocities, contrasting with the clay‐rich sediments in the shallow continental shelf and continental slope. A difference in S‐wave velocities is observed for upper oceanic crust formed at fast‐rate (Shumagin) and intermediate‐rate (Semidi) spreading centers. The reduced S‐wave velocities in the Semidi crust may be caused by increased faulting and possible lithological variations, related to a previous period of intermediate‐rate spreading.more » « less
-
In subduction zones, along‐strike and downdip variations in megathrust slip behavior are linked to changes in the properties of the subducting and overriding plates. Although marine geophysical methods provide insights into subduction zone structures, most surveys consist of sparse 2D profiles, limiting our understanding of first‐order controls. Here, we use active‐source seismic data to derive a 3D crustal‐scale P‐wave velocity model of the Alaska Peninsula subduction zone that encompasses both plates and spans the Semidi segment and SW Kodiak asperity. Our results reveal modest variations within the incoming plate, attributed to a series of fracture zones, seamounts and their associated basement swell, collectively contributing to plate hydration. Basement swell appears to modulate the distribution and type of sediment entering the trench, likely impacting observed variations in slip behavior. The overriding plate exhibits significant heterogeneity. The updip limit and width of the dynamic backstop are similar between the SW Kodiak asperity and eastern Semidi segment but differ significantly from the Western Semidi segment. These distinctions likely account for differences in earthquake rupture patterns and interseismic coupling among these segments. Additionally, high velocities in the mid‐lower forearc crust coincide with the location of the megathrust slip during the Mw 8.2 2021 Chignik event. We interpret these velocities as intracrustal intrusions that contributed to the deep rupture of the 2021 event. Our findings suggest that the contrasting structural and material properties of both the incoming and overriding plates influence the spatially complex and semi‐persistent segmentation of the megathrust offshore the Alaska Peninsula.more » « less
-
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
-
Abstract The structure of the mantle transition zone (MTZ) in subduction zones is essential for understanding subduction dynamics in the deep mantle and its surface responses. We constructed the P (Vp) and SH velocity (Vs) structure images of the MTZ beneath Northeast Asia based on two‐dimensional (2‐D) triplicated waveform modeling. In the upper MTZ, a normalVpbut 2.5% lowVslayer compared with IASP91 are required by the triplication data. In the lower MTZ, our results show a relatively higher‐velocity layer (+2%Vpand −0.5%Vscompared to IASP91) with a thickness of ~140 km and length of ~1,200 km atop the 660‐km discontinuity. Taking this anomaly as the stagnant slab and considering the plate convergence rate of 7–10 cm/year in the western Pacific region during the late Cenozoic, we deduced that the stagnant slab has a subduction age of less than 30 Ma. This suggests that the observed stagnancy of the slab in the MTZ beneath Northeast Asia may have occurred no earlier than the Early Oligocene. From the constraints derived individually onVpandVsstructures, highVp/Vsratios are obtained for the entire MTZ beneath Northeast Asia, which may imply a water‐rich and/or carbonated environment. Within the overall higher‐velocity stagnant slab, a low‐velocity anomaly was further detected, with a width of ~150 km,VpandVsreductions of 1% and 3% relative to IASP91. Such a gap may have provided a passage for hot deep mantle materials to penetrate through the thick slab and feed the Changbaishan volcano.more » « less