Although prevalent for the late Holocene, relative sea level (RSL) constraints during and immediately after the Last Glacial Maximum (LGM) are sparse. This scarcity of data is particularly pronounced along mid-latitude shelves such as central California, which lack post LGM RSL constraints older than 12 ka. In this study we collected 7 sediment cores and high-resolution seismic data from Estero Bay to constrain RSLs across the central California shelf between similar to 9 and similar to 16 ka. We reconstructed these RSLs using two sea-level indicators found within our sediment cores: the wave ravinement shell hash burial surface (WRSHBS) and the sedimentary contact between offshore mud facies and ripple cross-laminated sands. To determine the indicative meaning of these two sea-level indicators, we examined the relationship between the local wave regime, modern bathymetric profiles, and the depth of preservation of each sea-level indicator. After correcting for tectonic uplift, we estimated sea levels in central California to have been similar to 39 +/- 7.5 and 49 +/- 7.5 m below present sea level between 9 and 12 ka, in agreement with previous RSL reconstructions along this coast. Between 13.8 and 15.9 ka, we estimate sea levels to have reached similar to 86 +/- 8-99 +/- 8 m below present sea level. Our findings offer a Late Pleistocene RSL reconstruction for central California and develop new methodologies for estimating past RSLs on similar mid-latitude shelves.

The distribution and intensity of fault damage zones provides insight into fault activity and its relationship to fluid flow in the crust. Presently, measures of the in‐situ distribution of fault damage remain limited and along‐strike studies are rare. This study focuses on an offshore section Palos Verdes Fault damage zone that spans 28 km, near Los Angeles, California. To investigate the previously unresolved shallow (∼400 m below the seafloor) fault damage zone we use densely spaced (∼500 m line separation) newly collected sparker multichannel seismic lines and sub‐bottom profiles. The combination of high‐resolution acquisition methods and specialized seismic processing workflows provide improved imaging of shallow faulting. We apply a multi‐trace similarity technique to identify discontinuities in the seismic data that may be attributed to faults and fractures. This fault detection approach reveals diverse fault damage patterns on adjacent seismic profiles. However, a discernible damage zone pattern emerges by stacking multiple damage detection profiles along strike. We find that peak damage identified in this way corresponds to the active main fault strand, confirmed in this study, and thus the technique may be useful for identifying active fault strands elsewhere. Additionally, we observe that the variable width of the damage zone along strike is controlled by fault obliquity. Furthermore, our observations reveal a correlation between fault damage and seafloor fluid seeps visible in the water column, suggesting that damage plays a role in controlling fluid flow around the fault.

Abstract It is critical to characterize submarine landslide hazards near dense coastal populations, especially in areas with active faults, which can trigger slope failure, subsequent tsunamis, and damage seabed infrastructure during earthquake shaking. Offshore southern California, numerous marine geophysical surveys have been conducted over the past decade, and high‐resolution bathymetric and subsurface data now cover about 60 percent of the total region between Point Conception and the United States‐Mexico border from the California coast out to the base of Patton Escarpment ∼200 km offshore. In a comprehensive compilation and interpretive mapping effort, we find evidence of seafloor failure throughout offshore southern California with nearly 1,500 submarine landslide‐related features, including 63 discrete slide deposits with debris and >1,400 slide‐related scarps. In our analysis, we highlight new mapping of submarine landslides in Catalina Basin, the Del Mar slide, the San Gabriel slide complex, and the 232 km2 San Nicolas slide, the largest area of any known submarine landslide mass offshore southern California. Analysis of the spatial distribution of submarine landslide features suggests that most mapped slide features are located relatively near coastal sediment sources, particularly during sea‐level lowstand conditions, which underscores the importance of sediment supply and sediment accumulation on low‐gradient slopes as failure preconditioning processes. Tectonically driven uplift at shelf edges and along basin flanks is another key preconditioning factor, and our results also suggest that earthquakes along active faults trigger mass wasting, especially for repeated, small‐scale failures on tectonically steepened slopes.

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.

AbstractSeismic reflection data collected within the US Exclusive Economic Zone (EEZ) represent a scientific resource with inestimable value. The US Geological Survey (USGS), in partnership with the Bureau of Ocean Energy Management (BOEM), provide unrestricted public access to seismic data acquired since 1975 through the National Archive of Marine Seismic Surveys (NAMSS) portal. This map‐searchable web portal provides access to over 1,100 two‐ and three‐dimensional (2D, 3D) marine seismic data sets collected by the USGS and by the offshore resource industry. The industry data sets are released by BOEM under The Code of Federal Regulations 30 CFR Parts 551 and 580, allowing for the release of seismic reflection data after the completed 25‐year proprietary term from the issuance of the Geological and Geophysical (G&G) permit. The web portal provides tools to filter and identify surveys within the US EEZ and navigation and shapefiles can be downloaded with each survey. Seismic data sets downloaded from NAMSS have been instrumental in numerous scientific investigations along the US EEZ. These include studies focused on offshore marine geohazards, such as fault systems and submarine landslides, as well as gas hydrates and other geologic processes. Acquisition of such data sets, especially pneumatic‐sourced surveys, is becoming progressively more difficult due to environmental concerns and associated costs, thus preservation and public access to such high‐quality data sets is critical to maintaining the ability to conduct offshore seismic reflection scientific research along the US EEZ.

The Cross-Hosgri slope is a bathymetric lineament that crosses the main strand of the Hosgri fault offshore Point Estero, central California. Recently collected chirp seismic reflection profiles and sediment cores provide the basis for a reassessment of Cross-Hosgri slope origin and the lateral slip rate of the Hosgri fault based on offset of the lower slope break of the Cross-Hosgri slope. The Cross-Hosgri slope is comprised of two distinct stratigraphic units. The lower unit (unit 1) overlies the post-Last Glacial Maximum transgressive erosion surface and is interpreted as a Younger Dryas (ca. 12.85-11.65 ka) shoreface deposit based on radiocarbon and optically stimulated luminescence (OSL) ages, Bayesian age modeling, seismic facies, sediment texture, sediment infauna, and heavy mineral component. The shoreface was abandoned and partly eroded during rapid sea-level rise from ca. 11.5 to 7 ka. Unit 2 consists of fine sand and silt deposited in a midshelf environment when the rate of sea-level rise slowed between ca. 7 ka and the present. Although unit 2 provides a thin, relatively uniform cover over the lower slope break of the older shoreface, this feature still represents a valuable piercing point, providing a Hosgri fault slip rate of 2.6 +/- 0.8 mm/yr. Full-waveform processing of chirp data resulted in significantly higher resolution in coarser-grained strata, which are typically difficult to interpret with more traditional envelope processing. Our novel combination of offshore radiocarbon and OSL dating is the first application to offshore paleoseismic studies, and our results indicate the utility of this approach for future marine neotectonic investigations.

File S1 consists of additional sediment core and chirp figures. File S2 provides additional details and figures on the optical stimulated luminescence sampling and dating. File S3 provides details on the Bayesian age modeling conducted.

Subaqueous clinoforms are an important shelf feature. Their origins are typically associated with subaerial deltas but recent work has identified similar features in settings without a significant fluvial source (Mitchell, 2012; Patruno et al., 2015; Patruno and Helland-Hansen, 2018). These other studies have shown that such subaqueous clinoforms, also known as infralittoral prograding wedges (IPWs), are created largely by wave-induced processes. This study uses geophysical, sedimentological, and radiocarbon data to determine the sedimentary character-istics and genesis of a shore-parallel subaqueous clinoform developed far from any significant river on the central California continental shelf; a sedimentary shelf buildup known locally as the Cross Hosgri Slope. Sediment cores through the clinoformal sedimentary unit reveal that it is composed of beds with an erosive base, overlain by a thin (-15 cm) coarsening upward sequence of shelly fine sands transitioning to a fining upward sequence marked by alternating parallel and ripple cross laminated very fine sands. The sedimentary succession is often capped by fine silts that are commonly interbedded with thin very fine sand beds. Radiocarbon dating of shells collected just above the erosive base indicate the subaqueous clinoform initiated progradation-7 ka, nucleating on an older Younger Dryas relict shoreface. We suggest the Cross Hosgri Slope was created by winter-storm waves mobilizing sands in water depths up to-70 m that transitioned into wave-supported gravity flows. The wave-supported gravity flows traveled downslope to water depths of up to-80 m, corresponding to the foot of the subaqueous clinoform. They did not travel beyond this depth as wave influence at these depths is negli-gible and the shelf gradient is insufficient to maintain movement of the load alone. Our work suggests that wave -supported gravity flows can entrain very fine sands and silts and build subaqueous clinoforms, even in the absence of a significant river source. Furthermore, we provide a facies model for sandy wave-supported gravity flow deposits.

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.

First posted November 2, 2022 For additional information, contact: Director, Woods Hole Coastal and Marine Science CenterU.S. Geological Survey384 Woods Hole RoadQuissett CampusWoods Hole, MA 02543–1598 In 2015, U.S. Geological Survey scientists in collaboration with scientists from other institutions began a study of the Queen Charlotte fault—the first systematic study of the fault in more than three decades. The primary goal of the study was to gain a better understanding of the earthquake, tsunami, and underwater-landslide hazards throughout southeastern Alaska, as well as gather data to develop geologic models that can be applied to similar plate boundaries around the globe, such as the San Andreas fault system in southern California, the Alpine fault in New Zealand, and the North Anatolian fault in Turkey. A bathymetric terrain model was compiled from six different multibeam surveys of the previously unmapped Queen Charlotte fault offshore of southeastern Alaska and Haida Gwaii archipelago.

Abstract High‐resolution bathymetry and three‐dimensional seismic data along the Cocos Ridge reveal a 245 km2 field of ∼1–4 km in diameter seafloor depressions. The seafloor depressions are part of a two‐tiered honeycomb pattern. The lower‐tier depressions have steep faults that truncate strata with chaotic internal reflections consistent with sediment collapse into the depression. These extend into a lens shaped interval just above igneous basement. Overlying these depressions is a second broader set with rough seafloor morphology with gently dipping boundaries defined by pinch‐out stratigraphic patterns. Drilling results indicate that the lens‐shaped zones that host the deeper depressions represent anomalous regions of high porosity, low velocity, and low density within calcareous rich sediment. Analysis of nannofossils from IODP Site U1414 suggests the collapse structures formed during the late Miocene, whereas the younger shallower depressions likely formed between the early Pliocene and the Pliocene‐Pleistocene boundary. Geochemical and petrological analysis at Site U1414 suggests that hydrothermal circulation during the late Miocene led to carbonate dissolution and collapse. Following collapse, focused fluid‐flow and bottom current scouring resulted in formation of the overlying set of depressions and a honeycomb seafloor morphology. Similar sets of depressions along the Carnegie Ridge to the south support the hypothesis that two‐tiered depressions formed in response to processes that occurred broadly across the Panama Basin between the late Miocene and the Pliocene‐Pleistocene transition. Geochemical results at Site U1414, combined with geophysical data, suggest this two‐tiered system of depressions currently guides ongoing fluid outflow.

Fault damage zones provide a window into the non-elastic processes of an earthquake. Geological and seismic tomography methods have been unable to measure damage zones at depth with sufficient spatial sampling to evaluate the relative influence of depth, distance, and lithological variations. Here, we identify and analyze the damage zone of the Palos Verdes Fault offshore southern California using two 3D seismic reflection datasets. We apply a novel algorithm to identify discontinuities attributed to faults and fractures in large seismic volumes and examine the spatial distribution of fault damage in sedimentary rock surrounding the Palos Verdes Fault. Our results show that damage through fracturing is most concentrated around mapped faults and decays exponentially to a distance of ∼2 km, where fracturing reaches a clearly defined and relatively undamaged background for all examined depths and lithologies (450 m to 2.2 km). This decrease in fracturing with distance from the central fault strand exhibits similar functional form to outcrop studies. However, here we extend analysis to distances seldom accessible (∼10 km lateral distance). Separating the data by geologic units we find that the damage decay and background level differs for each unit, with the older and deeper units having higher levels of background fracturing and shallower exponential decays of fracturing with distance from the fault. Surprisingly, these differences in damage decay and background level trade-off result in a consistent damage zone width regardless of lithology or depth. We find that the damage zone has similar decay trends on both sides of the fault. When examining the damage zone at shorter (4 km vs 17 km) along strike distances, the damage zone has a more complex decay trend and at least two strands are resolvable.

ABSTRACT The Seattle fault zone (SFZ) is a north-directed thrust fault system that underlies the greater Seattle metropolitan area. Evidence of past land level changes, landslides, liquefaction, and a local tsunami indicate that this 70-km-long fault system can host up to M 7–7.5 earthquakes. Both the geometry and earthquake recurrence of the SFZ are debated and surveys of the shallow subsurface have not yet been incorporated into deeper crustal-scale structural interpretations, especially where the SFZ cuts across marine portions of the Puget Lowland. Here we use a new high-resolution marine seismic reflection dataset to image fault-related deformation in Quaternary sediments and Tertiary bedrock throughout Puget Sound and Lake Washington. We use this perspective of shallow geology as a link between existing crustal-scale geophysical insights into fault geometry at depth and paleoseismological observations of faulting at the surface and propose a refined structural model for the SFZ. We interpret that our new seismic reflection data in the Rich Passage area of Puget Sound images evidence of an inactive, south-dipping strand of the SFZ, which is overprinted by Quaternary folding and slip along north-dipping backthrusts within the hanging wall of a blind, south-dipping fault located 6 km farther north. To explain these results, we propose that the SFZ is a normal sequence fault propagation fold that has stepped northward through time, and we show the plausibility of this model through trishear forward modeling. Growth strata and faulting imaged in Quaternary sediments in Lake Washington and Rich Passage are consistent with the spatial distribution of folding and backthrusting that occurred during an M 7–7.5 earthquake in A.D. 900–930, corroborating existing evidence that the SFZ has been active throughout the Quaternary.

First posted July 30, 2021 For additional information, contact: Pacific Coastal and Marine Science CenterU.S. Geological SurveyPacific Coastal and Marine Science Center2885 Mission St.Santa Cruz, CA 95060 The California Continental Borderland is the complex continental margin in southern California that extends from Point Conception southward into northern Baja California (Fisher and others, 2009). This colored shaded-relief bathymetry map of the northern continental borderland in southern California was generated primarily from multibeam-echosounder data collected by the University of Washington in 2016, the Ocean Exploration Trust-Nautilus Exploration Program in 2015–17, and the National Oceanic and Atmospheric Administration in 2017. These datasets were processed in part by the U.S. Geological Survey. Additional smaller amounts of publicly available multibeam-bathymetry data collected by other federal and local agencies, academic institutions, and private firms were also incorporated into this map. Since the production of this map, other multibeam-bathymetry data have been collected in this region.

Abstract The Santa Barbara Channel, offshore California, contains several submarine landslides and ample evidence for incipient failure. This region hosts active thrust and reverse faults that accommodate several mm/yr of convergence, yet the relationships between tectonic deformation and slope failure remain unclear. We present 3‐D and 2‐D multichannel seismic reflection (MCS) data sets, multibeam bathymetry, and chronostratigraphic constraints to investigate the controls on slope failure. Splay faulting along the North Channel Deformation Trend (NCDT) coincides with a distinct zone of compressional uplift and onlapping of steeply dipping Quaternary strata. The NCDT is spatially correlated with seafloor fissures, and 3‐D seismic analyses reveal an intricate system of en echelon reverse faults that offset sediments younger than ~25 ka. Localized uplift zones are located between faults, one of which underlies the Gaviota landslide headscarp. We observe a direct relationship between slope failure and along‐strike variations in the tectonostratigraphic framework. Based on geophysical properties at Ocean Drilling Program (ODP) Site 893, we predict a trend in compaction and porosity reduction in the basin that drives pore fluids up‐dip, toward the zone of onlap above the NCDT, thus reducing slope stability. This interplay between tectonic, sedimentary, and fluid‐flow processes along the NCDT has created a confluence of preconditioning factors, with Gaviota and Goleta landslides being distinguished from the surrounding slopes by their position above the NCDT. The distribution of seafloor fissures suggests sections of the slope remain unstable and are prone to future landsliding. These results provide insights into the processes and 3‐D feedbacks that lead to slope instability along other convergent margins.

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.

Megathrust roughness and structural complexity are thought to be controls on earthquake slip at subduction zones because they result in heterogeneity in shear strength and resolved stress. However, because active megathrust faults are difficult to observe, the causes and scales of complexity are largely unknown. Here we measured the in situ properties of the megathrust of the Middle America subduction zone in a three-dimensional seismic reflection volume to determine how fault properties vary. We quantify spatial variability in the megathrust roughness, overburden and rock physical properties. Heterogeneity in the megathrust roughness exists at length scales of a few kilometres because the megathrust is dissected by active lower-plate normal faults, which offset the megathrust and renewed fault roughness. Spatial variations in the rock physical properties at the plate interface are characterized by correlation length scales of hundreds of metres. Frontal prism taper, historical seismicity and the variation in earthquake stress drop values local to the megathrust are all affected by the heterogeneity at these length scales. Both geometric and rheological complexities may therefore control the mechanical behaviour of the subduction plate interface, which includes earthquake rupture characteristics. Geometric and rheological complexities may control the mechanical behaviour of megathrusts, according to an analysis of the heterogeneity in roughness and rock properties of the Middle America megathrust from 3D seismic reflection data.