Abyssal marine turbidites provide some of the longest and most spatially extensive records of subduction zone earthquake recurrence globally; however, correlation of these deposits over long distances and interpretation of synchronous emplacement requires both an understanding of the turbidite generating systems and precise dating. Here, we present an integrated suite of high-resolution bathymetry, subbottom profiles, and sediment cores from combined autonomous underwater vehicle, remotely operated vehicle, and ship-based studies at a key paleoseismic site in the southern Cascadia subduction zone. We demonstrate how widespread, earthquake-triggered landslides on the lower slope deposit discrete, proximal mass transport deposits (MTDs) that grade offshore into complex, interfingered abyssal turbidites, which correspond to records of megathrust earthquake history. We propose accretion and oversteepening of thrust folds on the lower slope both preconditions the slope to fail and provides a perpetual source of unstable material to fail during every earthquake cycle. Furthermore, we suggest the periodic and pervasive landsliding indicates coseismic deformation of the outer accretionary wedge during megathrust rupture.

Ozette Lake is an ~100-m-deep coastal lake located along the outer coast of the Olympic Peninsula (Washington, USA); it is situated above the locked portion of the northern Cascadia megathrust but also relatively isolated from active crustal faults and intraslab earthquakes. Here we present a suite of geophysical and geological evidence for earthquake-triggered mass transport deposits (MTDs) and related turbidite deposition in Ozette Lake since ca. 14 ka. Comprehensive high-resolution bathymetry data, seismic reflection profiles, and sediment cores are used to characterize the post-glacial stratigraphic framework and examine paleoseismic evidence in the lacustrine sediments. Stacked sequences of MTDs along the steep eastern flanks of the lake appear to grade basin-ward from thick, chaotic, blocky masses to thin, parallel-bedded turbidite beds. The discrete turbidite event layers are separated by fine-grained (silt and clay) lake sedimentation. The event layers are observed throughout the lake, but the physical characteristics of the deposits vary considerably depending on proximity to primary depocenters, steep slopes, and subaqueous deltas. A total of 30–34 event deposits are observed in the post-glacial record. Radiometric dating was used to reconstruct a detailed sedimentation history over the last ~5.5 k.y., develop an age model, and estimate the recurrence (365–405 yr) for the most recent 12 event layers. Based on sedimentological characteristics, temporal overlap with other regional paleoseismic chronologies, and recurrence estimates, at least 10 of the dated event layers appear to be sourced from slope failures triggered by intense shaking during megathrust ruptures; the recurrence interval for these 10 events is 440–560 yr. Thus, Ozette Lake contains one of the longest and most robust geological records of repeated shaking along the northern Cascadia subduction zone.

Intraslab earthquakes do not produce primary paleoseismic evidence at the Earth’s surface, making efforts to develop an event chronology challenging. However, the strong ground motion from intraslab events may initiate gravity-driven turbidity flows in subaqueous basins; the resulting deposits (turbidites) can provide a paleoseismic proxy if the conditions that initiate these flows are known. To better constrain the initiating conditions, we use two recent intraslab earthquakes in southcentral Alaska, the Mw 7.1 November 30, 2018, Anchorage and the Mw 7.1 January 24, 2016, Iniskin earthquakes, as calibration events. Through a multi-lake investigation spanning a range of shaking intensities and based on a combined geological and geophysical dataset, we document the occurrence, or absence, of earthquake-generated turbidity flows from these two earthquakes. The 2018 and 2016 earthquakes are recorded by centimeter-scale turbidites that can be differentiated from climatically generated deposits, as well as other seismic sources (i.e., the 1964 Alaska megathrust earthquake) based on deposit thickness, sedimentological properties, and deposit age. We show that a Modified Mercalli Intensity (MMI) of ~V-V1/2 is the minimum shaking intensity required to generate localized sediment remobilization from deltaic slopes, and a MMI of ~V1/2 is required to produce a deposit of sufficient thickness that a seismic origin can be confidently assigned. Deltaic slopes are the major source of remobilized sediment that record the 2018 and 2016 events, however sediment from non-tributary sourced basin slopes may become remobilized in steep-sloped, high sedimentation areas, and under elevated shaking intensity. The documentation of seismically generated deposits in quick succession (~2 years) with diagnostic features that can be assigned to the seismic source highlights the utility of using recent earthquakes as calibration events to investigate the subaqueous response to strong ground motion.

Previous margin-wide studies of methane seep distribution along the Cascadia Subduction Zone indicate peaks in seep density within the landward limit of the of gas hydrate stability zone (GHSZ; ≤500 m depth), suggesting a link between current ocean warming, acceleration of hydrate dissociated, and methane emissions. This inferred connection, however, may not account for regional geologic and/or structural complexities driving methane seepage. Expanding upon an existing seep database by adding new seep data, we conducted statistical and spatial analyses to determine margin-wide distribution trends and offer a tectonic framework for understanding the tendency toward non-normality and spatial clustering. We then highlight the role of local-scale drivers of seep formation in addition to the first-order tectonic framework, using systematic geologic/geomorphic characterization of seep emission sites in southern Cascadia and case studies using meta-attribute analysis of seismic reflection data. Seep distribution along the margin is non-random, but instead of clustering along the 500-m isobath, regions of high seep density occur in canyons and topographic highs. New findings from this study conclude that co-location of the outer arc high (OAH) and the landward limit of the GHSZ may explain high concentrations of seeps where deformation is the greatest and hydrates are unstable. Detailed analysis of the spatial relationships between seep sites and geologic-geomorphic features in southern Cascadia reveal a link between seeps and anticlines, with 52% of the seeps found in association with anticlines, 36% found at faults, 16% associated with canyons, and 11% at seafloor failure scarps. Given that a majority of anticlines are located along or seaward of the OAH in the actively deforming outer wedge, we suggest that the location of the OAH is a primary structural control on seep distribution. This scenario is supported by neural network analysis of multichannel seismic data revealing zones of probable fluid migration along vertical pipes, faults, and chimneys in the vicinity of active seep sites on anticlines. Determining linkages between seeps and submarine tectonic geomorphology is a crucial first step for understanding and forecasting the distribution of methane seepage, but also a necessity for evaluating causal relationships between ocean warming and gas hydrate stability.

•Mass wasting is pervasive along the steep lower slope in Cascadia.•Earthquake triggered turbidity flows may be primarily sourced from the lower slope.•Fine-grained turbidites in Cascadia may represent smaller magnitude earthquakes.•The most intense earthquake shaking may occur closest to the deformation front.

High resolution seafloor mapping shows extraordinary evidence that massive (>300 m thick) icebergs once drifted >5,000 km south along the eastern United States, with >700 iceberg scours now identified south of Cape Hatteras. Here we report on sediment cores collected from several buried scours that show multiple plow marks align with Heinrich Event 3 (H3), ~31,000 years ago. Numerical glacial iceberg simulations indicate that the transport of icebergs to these sites occurs during massive, but short-lived, periods of elevated meltwater discharge. Transport of icebergs to the subtropics, away from deep water formation sites, may explain why H3 was associated with only a modest increase in ice-rafting across the subpolar North Atlantic, and implies a complex relationship between freshwater forcing and climate change. Stratigraphy from subbottom data across the scour marks shows there are additional features that are both older and younger, and may align with other periods of elevated meltwater discharge.

The lack of historic seismicity along the Cascadia subduction zone (CSZ) fuels the uncertainty of the megathrust rupture process and resulting hazards in the US and Canadian Pacific Northwest. Diverse types of paleoseismic proxies uniquely record different aspects of the earthquake rupture process. Thus, margin-wide integration of paleoseismic proxy records leverages these unique transcripts to provide a richer understanding of past earthquake variability in time and space. In a series of virtual workshops, the USGS Powell Center research group on CSZ earthquake science and hazards have compiled paleoseismic records of coseismic coastal land-level change, tsunami inundation, and ground shaking proxies such as landslides, liquefaction, and turbidites. To consistently integrate different datasets, we constructed ranking schemes to quasi-quantitatively evaluate the quality of geochronologic age control and certainty of earthquake evidence at 49 coastal land level change sites, 81 tsunami sites, 76 terrestrial paleoseismic shaking sites, and 13 marine paleoseismic shaking sites. At each site, we evaluated how well geochronologic samples and results (14C dates, stratigraphy, geological context etc.) constrain the timing of candidate earthquake events. Similarly, we applied an evidence quality ranking scheme to assess how certain the observations of coastal land-level change, tsunami inundation, and strong ground motion were the result of a CSZ megathrust earthquake source. We also focused on historic records of non-local events, such as the 1964 Alaska earthquake and the 1992 Cape Mendocino earthquake, to identify the characteristics of distal subduction zone earthquakes and small local events, respectively, in the paleoseismic record. In the future we plan to identify the likely rupture characteristics of major CSZ earthquakes during the Holocene and to synthesize paleoseismic, geophysical, and instrumental datasets to identify the existence and persistence of rupture boundaries along the subduction zone.

Marine turbidite records have been used to infer palaeoseismicity and estimate recurrence intervals for large (>M(w)7) earthquakes along the Cascadia Subduction Zone. Conventional models propose that upper slope failures are funneled into submarine canyons and develop into turbidity flows that are routed down-canyon to deep-water channel and fan systems. However, the sources and pathways of these turbidity flows are poorly constrained, leading to uncertainties in the connections between ground shaking, slope failure and deep-water turbidites. We examine the spatial distribution of submarine landslides along the southern Cascadia margin to identify source regions for slope failures that may have developed into turbidity flows. Using multibeam bathymetry, sparker multichannel seismic and chirp sub-bottom data, we observe relatively few canyon head slope failures and limited evidence of large landslides on the upper and middle slope. Most of the submarine canyons are draped with sediment infill in the upper reaches and do not appear to be active sediment conduits during the recent sea-level highstand. In contrast, there is evidence of extensive mass wasting of the lower slope and non-channelized downslope flows. Contrary to previous studies, we propose that failures along the lower slope are the primary sources for deep-sea seismoturbidites in southern Cascadia.

Abstract The Cape Fear Slide is one of the largest (>25 000 km3) submarine slope failure complexes on the US Atlantic margin. Here we use a combination of new high-resolution multichannel seismic data (MCS) from the National Science Foundation Geodynamic Processes at Rifting and Subducting Margins (NSF GeoPRISMS) Community Seismic Experiment and legacy industry MCS to derive detailed stratigraphy of this slide and constrain the conditions that lead to slope instability. Limited outer-shelf and upper-slope accommodation space during the Neogene, combined with lowstand fluvial inputs and northwards Gulf Stream sediment transport, appears to have contributed to thick Miocene and Pliocene deposits that onlapped the lower slope. This resulted in burial of an upper-slope bypass zone developed from earlier erosional truncation of Paleogene strata. These deposits created a broad ramp that allowed accumulation of thick Quaternary strata across a low-gradient (<3.5°) upper slope. Upslope of one of the larger headwalls, undulating Quaternary strata appear to downlap onto a buried failure plane. Many of the nested headwalls of the upper-slope portion of slide complex are underlain by deformed strata, which may be the result of fluid migration associated with localized subsidence from salt migration. These new data and observations suggest that antecedent margin physiography, sediment loading and substrate fluid flow were key factors in preconditioning the Cape Fear slope for failure.

Linear scour depressions observed along the upper slope offshore of South Carolina are among the largest and most extensive iceberg ploughmarks found at subtropical latitudes. These features record the drifting of large, tabular icebergs south from the Laurentide Ice Sheet (LIS), along the entire length of US Atlantic margin during the last deglaciation (Hill et al. 2008). Over 400 individual, kilometre-length-scale furrows are observed in swath bathymetry along the upper slope (170–380 m water depth) of the southern US Atlantic margin (Fig. 1). The furrows are present along the entire continental slope from Cape Hatteras ( c. 35° N) to the Florida Straits ( c. 24.5° N), but are most prominent and well documented along a shoaling platform offshore of South Carolina ( c. 32.5° N; Hill & Condron 2014). The furrows in this region are orientated predominantly WSW, sub-parallel to regional bathymetric contours (Fig. 1c). While some larger furrows are up to 400 m wide and 20 m deep, the scoured troughs are generally 10–100 m wide and <10 m deep (Fig. 1a, e). The furrows range from linear to arcuate with numerous cross-cutting tracks, particularly in deeper water where the ploughmarks are more densely spaced. Several of the more prominent linear furrows are traceable for more than 30 km along the margin. Many of …