Description: Using older and in part fl awed data, Ruff (1989) suggested that thick sediment entering the subduction zone (SZ) smooths and strengthens the trench-parallel distribution of interplate coupling. This circumstance was conjectured to favor rupture continuation and the generation of high-magnitude (≥Mw8.0) interplate thrust (IPT) earthquakes. Using larger and more accurate compilations of sediment thickness and instrumental (1899 to January 2013) and pre-instrumental era (1700–1898) IPTs (n = 176 and 12, respectively), we tested if a compelling relation existed between where IPT earthquakes ≥Mw7.5 occurred and where thick (≥1.0 km) versus thin (≤1.0 km) sedimentary sections entered the SZ. Based on the new compilations, a statistically supported statement (see Summary and Conclusions) can be made that high-magnitude earthquakes are most prone to nucleate at well-sedimented SZs. For example, despite the 7500 km shorter global length of thicksediment trenches, they account for ~53% of instrumental era IPTs ≥Mw8.0, ~75% ≥Mw8.5, and 100% ≥Mw9.1. No megathrusts 〉Mw9.0 ruptured at thin-sediment trenches, whereas three occurred at thick-sediment trenches (1960 Chile Mw9.5, 1964 Alaska Mw9.2, and 2004 Sumatra Mw9.2). However, large Mw8.0–9.0 IPTs commonly (n = 23) nucleated at thin-sediment trenches. These earthquakes are associated with the subduction of low-relief ocean floor and where the debris of subduction erosion thickens the plate-separating subduction channel. The combination of low bathymetric relief and subduction erosion is inferred to also produce a smooth trench-parallel distribution of coupling posited to favor the characteristic lengthy rupturing of highmagnitude IPT earthquakes. In these areas subduction of a weak sedimentary sequence further enables rupture continuation

Description: Large intraplate earthquakes in oceanic lithosphere are rare and usually related to regions of diffuse deformation within the oceanic plate. The 23 January 2018 MW 7.9 strike-slip Gulf of Alaska earthquake ruptured an oceanic fracture zone system offshore Kodiak Island. Bathymetric compilations show a muted topographic expression of the fracture zone due to the thick sediment that covers oceanic basement but the fracture zone system can be identified by offset N-S magnetic anomalies and E-W linear zones in the vertical gravity gradient. Back-projection from global seismic stations reveals that the initial rupture at first propagated from the epicenter to the north, likely rupturing along a weak zone parallel to the ocean crustal fabric. The rupture then changed direction to eastward directed with most energy emitted on Aka fracture zone resulting in an unusual multi-fault earthquake. Similarly, the aftershocks show complex behavior and are related to two different tectonic structures: (1) events along N-S trending oceanic fabric, which ruptured mainly strike-slip and additionally, in normal and oblique slip mechanisms and (2) strike-slip events along E-W oriented fracture zones. To explain the complex faulting behavior we adopt the classical stress and strain partitioning concept and propose a generalized model for large intra-oceanic strike-slip earthquakes of trench-oblique oriented fracture zones/ocean plate fabric near subduction zones. Taking the Kodiak asperity position of 1964 maximum afterslip and outer-rise Coulomb stress distribution into account, we propose that the unusual 2018 Gulf of Alaska moment release was stress transferred to the incoming oceanic plate from co- and post-processes of the nearby great 1964 MW 9.2 megathrust earthquake.

Description: Some of the largest earthquakes worldwide, including the 1964 9.2 Mw megathrust earthquake, occurred in Alaskan subduction zones. To better understand rupture processes and their mechanisms, we relate seafloor morphology from multibeam and regional bathymetric compilations with sub-seafloor images and seismic P-wave velocity structures. We re-processed legacy multichannel seismic (MCS) data including shot- and intra-shotgather interpolation, multiple removal and Kirchhoff depth migration. These images even reveal the shallow structure of the subducting oceanic crust. Traveltime tomography of a coincident vintage (1994) wide angle dataset reveals the P-wave velocity distribution as well as the deep structure of the subducting plate to the ocean crust Moho. The subducting oceanic crust morphology is rough and partly hidden by a thick sediment cover that reaches ~3 km depth at the trench axis. Bathymetry shows two major contrasting upper plate morphologies: the shallow dipping lower slope consists of trench-parallel ridges that form the accreted prism whereas the steep rough middle and upper slopes are composed of competent older rock. Thrust faults are distributed across the entire slope, some of which connect with the subducted plate interface. A subtle change in seafloor gradient from the lower to the middle slope coincides with a thrust fault zone marking the boundary between the margin framework and the frontal prism. It corresponds to the most prominent lateral increase in seismic P-wave velocities, ~25 km landward of the trench axis. Major thrusts in several MCS-lines are correlated with bathymetric data, showing their 〉 100 km lateral extent, which might also be tsunamigenic paths of earthquake rupture from the seismogenic zone to the seafloor.

Description: In 1946, megathrust seismicity along the Unimak segment of the Alaska subduction zone generated the largest ever recorded Alaska/Aleutian tsunami. The tsunami severely damaged Pacific islands and coastal areas from Alaska to Antarctica. It is the charter member of “tsunami” earthquakes that produce outsized far-field tsunamis for the recorded magnitude. Its source mechanisms were unconstrained by observations because geophysical data for the Unimak segment were sparse and of low resolution. Reprocessing of legacy geophysical data reveals a deep water, high-angle reverse or splay thrust fault zone that leads megathrust slip upward to the mid-slope terrace seafloor rather than along the plate boundary toward the trench axis. Splay fault uplift elevates the outer mid-slope terrace and its inner area subsides. Multibeam bathymetry along the splay fault zone shows recent but undated seafloor disruption. The structural configuration of the nearby Semidi segment is similar to that of the Unimak segment, portending generation of a future large tsunami directed toward the US West coast.

Description: Some of the largest earthquakes worldwide, including the 1964 9.2 Mw megathrust earthquake, occurred in Alaskan subduction zones. To better understand rupture propagation, with a focus on barriers that limit rupture, i.e. the southwest end of the Albatross segment at the boundary with the 1938 Semidi Mw 8.3 earthquake segment, we relate multibeam seafloor morphology with sub-seafloor images and seismic P-wave velocity structure. We re-processed legacy multichannel seismic (MCS) data including shot- and intra-shotgather interpolation, multiple removal and Kirchhoff depth migration and/or MCS traveltime tomography. These images even reveal the shallow structure of the subducting oceanic crust. Traveltime tomography of a coincident vintage wide-angle dataset reveals the P-wave velocity and the deep structure of the subducting plate to depths of the ocean crust Moho. The subducting oceanic crust morphology is rough and partly hidden by thick sediment cover that reaches ~3 km in the Albatross and ~1.5-2 km thickness in the Semidi segment at the trench axis. In both segments, bathymetry shows two major contrasting upper plate morphologies: trench-parallel ridges form the accreted prism of the shallow dipping lower slope whereas the steep rough middle and upper slopes are composed of competent older rock. Thrust faults are distributed across the entire slope, some of which connect with the subducted plate interface. A subtle change in seafloor gradient from the lower to the middle slope coincides with a splay fault zone (SFZ) marking the boundary between the margin framework and the frontal prism. This SFZ corresponds to the most prominent lateral increase in seismic P-wave velocities, ~25 km landward of the trench axis. Major differences in the Albatross/Semidi segments are 1.) origin of subducting sediment (Surveyor vs. Kodiak fan) 2.) geometry of the subducting plate (gentle vs. steep dip) beneath the lower and middle slopes.

Description: Some of the largest earthquakes worldwide, including the 1964 9.2 Mw megathrust earthquake, occurred in Alaskan subduction zones. To better understand rupture processes and their mechanisms, we relate seafloor morphology from multibeam and regional bathymetric compilations with sub-seafloor images and seismic P-wave velocity structures. We re-processed legacy multichannel seismic (MCS) data including shot- and intra-shotgather interpolation, multiple removal and Kirchhoff depth migration. These images even reveal the shallow structure of the subducting oceanic crust. Traveltime tomography of a coincident vintage (1994) wide angle dataset reveals the P-wave velocity distribution as well as the deep structure of the subducting plate to the ocean crust Moho. The subducting oceanic crust morphology is rough and partly hidden by a thick sediment cover that reaches ~3 km depth at the trench axis. Bathymetry shows two major contrasting upper plate morphologies: the shallow dipping lower slope consists of trench-parallel ridges that form the accreted prism whereas the steep rough middle and upper slopes are composed of competent older rock. Thrust faults are distributed across the entire slope, some of which connect with the subducted plate interface. A subtle change in seafloor gradient from the lower to the middle slope coincides with a splay fault zone marking the boundary between the margin framework and the frontal prism. It corresponds to the most prominent lateral increase in seismic P-wave velocities, ~25 km landward of the trench axis. Major thrusts in several MCS-lines are correlated with bathymetric data, showing their 〉 100 km lateral extent, which might also be tsunamigenic paths of earthquake rupture from the seismogenic zone to the seafloor. The splay fault zone has been recognized as a potential tsunamigenic structure in the 1938 and 1946 earthquake rupture areas.

Description: The potential for a major earthquake in the Shumagin seismic gap, and the tsunami it could generate, was reported in 1971. However, while potentially tsunamigenic splay faults in the adjacent Unimak and Semidi earthquake segments are known, such features along the Shumagin segment were undocumented until recently. To investigate margin structure and search for splay faults, we reprocessed six legacy seismic records and also processed seismic data acquired by RV Langseth during the ALEUT project (cf. Bécel et al., 2017). All records show splay faults separating the frontal prism from the margin framework. A ridge uplifted by the splay fault hanging wall extends along the entire segment. At the plate interface, the splay fault cuts across subducted sediment strata in some images, whereas in others, the plate interface sediment cuts across the fault. Splay fault zones are commonly associated with subducting lower-plate relief. Along the upper slope, beneath a sediment cover, major normal faults dipping landward and seaward border a ridge of basement rock. The faults displace a regional unconformity that elsewhere received Oligocene–Miocene sediment. Low seafloor scarps above some normal faults indicate recent tectonism. The buried ridge is a continuation of the Unimak Ridge structure that extends NE of the Unimak/Shumagin segment boundary. Some geological characteristics of the Shumagin segment differ from those of other Alaskan earthquake segments, but a causal link to the proposed Shumagin creeping seismic behavior is equivocal.