skip to main content


Title: Subduction Zone Interface Structure Within the Southern M W 9.2 1964 Great Alaska Earthquake Asperity: Constraints From Receiver Functions Across a Spatially Dense Node Array
Abstract

We conduct a high‐resolution teleseismic receiver function investigation of the subducting plate interface within the Alaskan forearc beneath Kodiak Island using data collected as part of the Alaska Amphibious Community Seismic Experiment in 2019. The Kodiak node array consisted of 398 nodal geophones deployed at ∼200–300 m spacing on northeastern Kodiak Island within the southern asperity of the 1964 Mw9.2 Great Alaska earthquake. Receiver function images at frequencies of 1.2 and 2.4 Hz show a coherent, slightly dipping velocity increase at ∼30–40 km depth consistent with the expected slab Moho. In contrast to studies within the northern asperity of the 1964 rupture, we find no evidence for a prominent low‐velocity layer above the slab Moho thick enough to be resolved by upgoing P‐to‐S conversions. These results support evidence from seismicity and geodetic strain suggesting that the 1964 rupture connected northern (Kenai) and southern (Kodiak) asperities with different plate interface properties.

 
more » « less
Award ID(s):
1654568 1950328
NSF-PAR ID:
10372459
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
49
Issue:
15
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The top of the subducting plate is a thick, complex zone where the heterogeneous structure likely controls earthquake rupture processes. Imaging this heterogeneous channel typically involves active‐source methods with limited depth penetration, or low‐resolution teleseismic methods. To access short wavelengths at greater depth, we use high‐frequency P‐to‐S (PS, 1–15 Hz) mode‐converted arrivals from nearby earthquakes >50 km deep to image the plate interface at vertical scales <1 km. We use 37 broadband stations in southcentral Alaska between 2007 and 2008 at 10–15 km spacing, spanning the great 1964 earthquake rupture zone and adjacent deeper slow slip and tremor regions. The central 21 stations record high‐amplitude PS arrivals converting from the megathrust region, at depths corresponding to the top of a prominent low‐velocity zone (LVZ) in receiver function images. The PS/P amplitude ratio (APS_P) varies along strike and with depth of the conversion point but is independent of earthquake location and varies slowly between adjacent stations. APS_Pchanges with slab depth, indicating changes in lithology or fluid content of the plate interface, consistent with transitions in slip behavior from locked to slow slip. High APS_Pcannot be explained by a velocity step or a single low‐velocity, high Vp/Vs layer, but requires several alternating high and low‐velocity layers. These observations indicate that the LVZ is a highly heterogeneous channel at multiple scales, resembling a subduction channel or sheared zone of metasediment and altered crust as observed in many exhumed subduction zones.

     
    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

    Seamounts are found at many subduction zones and act as seafloor heterogeneities that affect slip behavior on megathrusts. At the Hikurangi subduction zone offshore the North Island, New Zealand, seamounts have been identified on the incoming Pacific plate and below the accretionary prism, but there is little concrete evidence for seamounts subducted beyond the present‐day coastline. Using a high‐resolution, adjoint tomography‐derived velocity model of the North Island, we identify two high‐velocity anomalies below the East Coast and an intraslab low‐velocity zone up‐dip of one of these anomalies. We interpret the high‐velocity anomalies as previously unidentified, deeply subducted seamounts, and the low‐velocity zone as fluid in the subducting slab. The seamounts are inferred to be 10–30 km wide and on the plate interface at 12–15 km depth. Resolution analysis using point spread functions confirms that these are well‐resolved features. The locations of the two seamounts coincide with bathymetric features whose geometries are consistent with those predicted from analog experiments and numerical simulations of seamount subduction. The spatial characteristics of seismicity and slow slip events near the inferred seamounts agree well with previous numerical modeling predictions of the effects of seamount subduction on megathrust stress and slip. Anomalous geophysical signatures, magnetic anomalies, and swarm seismicity have also been observed previously at one or both seamount locations. We propose that permanent fracturing of the northern Hikurangi upper plate by repeated seamount subduction may be responsible for the dichotomous slow slip behavior observed geodetically, and partly responsible for along‐strike variations in plate coupling on the Hikurangi subduction interface.

     
    more » « less
  4. Abstract

    The crust and upper mantle beneath the New England Appalachians exhibit a large offset of the Moho across the boundary between Laurentia and accreted terranes and several dipping discontinuities, which reflect Paleozoic or younger tectonic movements. We apply scattered wavefield migration to the SEISConn array deployed across northern Connecticut and obtain insights not previously available from receiver function studies. We resolve a doubled Moho at a previously imaged Moho offset, which may reflect westward thrusting of rifted Grenville crust. The migration image suggests laterally variable velocity contrasts across the Moho, perhaps reflecting mafic underplating during continental rifting. A west‐dipping feature in the lithospheric mantle is further constrained to have a slab‐like geometry, representing a relict slab subducted during an Appalachian orogenic event. Localized low seismic velocities in the upper mantle beneath the eastern portion of the array may indicate that the Northern Appalachian Anomaly extends relatively far to the south.

     
    more » « less
  5. Abstract

    The Patagonian slab window has been proposed to enhance the solid Earth response to ice mass load changes in the overlying Northern and Southern Patagonian Icefields (NPI and SPI, respectively). Here, we present the first regional seismic velocity model covering the entire north‐south extent of the slab window. A slow velocity anomaly in the uppermost mantle indicates warm mantle temperature, low viscosity, and possibly partial melt. Low velocities just below the Moho suggest that the lithospheric mantle has been thermally eroded over the youngest part of the slab window. The slowest part of the anomaly is north of 49°S, implying that the NPI and the northern SPI overlie lower viscosity mantle than the southern SPI. This comprehensive seismic mapping of the slab window provides key evidence supporting the previously hypothesized connection between post‐Little Ice Age anthropogenic ice mass loss and rapid geodetically observed glacial isostatic uplift (≥4 cm/yr).

     
    more » « less