Search USGS⌕ Search

Geology topics

Research about California, Oregon, Washington

Source-linked reports with geographic coverage including California, Oregon, Washington.

At least 19 recordsLinked to original sources

A geographic dataset of rocky reefs habitat areas of particular concern for the United States West Coast

The United States National Marine Fisheries Service determines “essential fish habitat (EFH)” for federally managed species in coordination with regional fishery management councils, considers adverse effects to those habitats, and provides information to further habitat conservation and enhancement. Identifying discrete subsets of EFH as “habitat areas of particular concern (HAPC)” can help focus conservation, management, and research efforts. In 2006, the Pacific Fishery Management Council designated rocky reefs along the United States (U.S.) West Coast as HAPCs for groundfishes because of their ecological significance, sensitivity to human impacts, and relative rarity. To better understand where rocky reefs occur, we (1) located rocky reef areas that were not included in the 2006 rocky reef dataset, and (2) incorporated best available data into a refined geographic dataset that enables visualization. Our update shows that rocky reefs are distributed throughout the U.S. West Coast continental margin, are patchier than previously known, and comprise 8% of the extent of all data inputs. This updated dataset will inform resource management decisions in coastal and marine environments.

California, Oregon, Washington↗

Listening for extinction: Range-wide occupancy highlights critical risks and priorities for northern spotted owls

The northern spotted owl ( Strix occidentalis caurina ) continues to decline across its range, threatened by habitat disturbances and invasive barred owls ( Strix varia ). We conducted a range-wide spotted owl occupancy assessment using 2.1 million h of passive acoustic recordings from 4081 survey stations within 1027 randomly selected 5-km 2 hexagons across six physiographic regions. Applying multistate occupancy models, we estimated probabilities of landscape use (by at least one northern spotted owl) and pair occupancy (occupied by both a male and a female), while accounting for imperfect and sex-biased detection probability. Detection probability of pairs was low in most regions, primarily due to undetected females, but increased with survey effort. If management decisions can only be made when there is a high probability of detecting both sexes when they are present, then longer surveys are one approach to support the decisions. Landscape use and pair occupancy were positively associated with old-growth forest structure and topographic diversity. We also found a pronounced latitudinal gradient in spotted owl landscape use rates with estimated mean landscape use of 0.17 (SE = 0.03) and 0.16 (SE = 0.02) in the two most northern regions, and 0.54 (SE = 0.04) and 0.52 (SE = 0.05) in the two most southern regions. Pair occupancy rates were higher in southern regions (0.29, SE = 0.03, 0.43, SE = 0.05) and very low in northern regions (0.04, SE = 0.01, 0.03, SE = 0.01). In contrast to spotted owls, barred owls were detected in 90% of northern hexagons and 50% of southern hexagons, and barred owl calling intensity was over eight times higher than spotted owls. These results highlight extinction risks in several regions and suggest that management actions, such as spatially targeted barred owl control and protection of structurally complex habitats, have a narrow temporal window in which they can be effective in lowering the probability of regional extirpations. We generated predictive maps of landscape use and pair occupancy, which could help guide regional conservation strategies. This study demonstrates how passive acoustic monitoring and machine learning can be used for broad-scale ecological monitoring, identifying extinction risk thresholds, supporting adaptive management, and improving conservation outcomes for wide-ranging and elusive species that can reliably be detected through unique vocalizations.

California, Oregon, Washington↗

The Cascadia Region Earthquake Science Center (CRESCENT): Advancing understanding of Cascadia's earthquake hazards

The Cascadia Subduction Zone (CSZ) is capable of producing large megathrust earthquakes and tsunamis, posing significant risks to the Pacific Northwest. This paper examines the seismic characteristics of the CSZ, including its offshore coupling, sparse seismicity, and a complex paleoseismic record. While recent advances in geophysical methods and data collection have improved our understanding of subduction zone processes, important uncertainties remain regarding rupture dynamics, ground shaking, and tsunami impacts. The Cascadia Region Earthquake Science Center (CRESCENT) has been established to address these challenges through coordinated research, education, and community engagement. By integrating scientific studies, fostering collaboration, and leveraging new technologies, CRESCENT seeks to improve hazard assessment and resilience in the region. This paper outlines the primary challenges and opportunities in advancing subduction zone science and enhancing societal preparedness for future CSZ events.

California, Oregon, Washington↗

Spatially concentrating logging could mitigate climate-magnified fragmentation risks to a globally endangered bird

1. Rising timber demand is transforming forest structure globally, profoundly affecting biodiversity and climate resilience. Logging-driven fragmentation is potentially a major driver of biodiversity loss in production landscapes, yet its interactions with escalating climate stressors remain poorly understood. 2. We combine two decades of Landsat-derived habitat metrics with 29,000 surveys of the marbled murrelet ( Brachyramphus marmoratus )—an iconic Pacific Northwest old-forest specialist seabird affecting management of >10 million hectares. Controlling for habitat amount and detection probability, increasing landscape-scale forest edge amount sharply reduces murrelet occupancy, with impacts worsening under unfavourable climate-driven ocean conditions. 3. Comparing alternative landscape-scale timber harvest strategies, spatially concentrated logging consistently supports higher murrelet populations than fragmented approaches producing equivalent wood volumes, with benefits amplified under adverse ocean conditions. However, historical harvesting policies in the Pacific Northwest have instead driven severe habitat fragmentation, which we show is eroding the value of core set-aside forests on federal and conservation lands and ultimately rendering murrelets more vulnerable to climate change. 4. Synthesis and applications : We map key opportunities to boost populations by reducing edginess around remaining nesting habitat and investigate these opportunities' spatial distribution across land ownership and timber productivity gradients. Concentrating logging could be critical for mitigating fragmentation and climate threats for murrelets and potentially other forest-dependent species amid rising timber demand.

California, Oregon, Washington↗

ENSO and PDO drive shoreline position anomalies in the U.S. Pacific Northwest

Sandy beaches act as buffers against various coastal hazards but are vulnerable to episodic (seasonal) and chronic (interannual) erosion. Understanding the variation of shoreline position, a key metric in coastal morphology, over a spectrum of time scales is therefore crucial in assessing hazard vulnerability. Long-standing research has investigated the role of El Niño-Southern Oscillation (ENSO), the dominant mode of climate variability in the Pacific Basin, in seasonal shoreline variability. Yet, ENSO’s chronic influence—and that of another Pacific climate mode, the Pacific Decadal Oscillation (PDO)—on shoreline anomalies remains poorly understood. Here, we examine the variability of sandy beaches in the US Pacific Northwest, a ∼750 km long coastal region on the US West Coast. We leverage 40 years (1984–2024) of shoreline data from publicly available Earth-observing (Landsat) satellite imagery at a high spatial resolution (>10,000 shore-normal transects at 50-m alongshore spacing) and employ Convergent Cross Mapping (CCM), a methodology for inferring causality in dynamical systems. We discover that strong El Niño years are signified by erosion (75.1% of transects), and strong La Niña years exhibit accretional behavior (73.4% of transects). Furthermore, we establish, for the first time, that both ENSO and PDO exert a statistically significant control on interannual shoreline variability, particularly on the alongshore component (in 95 and 100% of littoral cells, respectively), with water level fluctuations playing a critical role. This effort advances our understanding of the seasonal-to-interannual interactions between Pacific Basin climate variability and the PNW’s coastal morphodynamics, with implications for sediment management and coastal adaptation.

California, Oregon, Washington↗

Fostering landscape resilience and species conservation in frequent fire landscapes

Restoration of dry forest landscapes of the Pacific Northwest often involves the removal of some closed-canopy forest that has expanded due to fire exclusion, but this can conflict with species conservation when species require closed-canopy forest as part of their habitat. Fire refugia—areas within fire perimeters that burn less severely or frequently—might provide a solution to this conflict. We used a species-specific approach to understanding fire refugia for northern spotted owls to model the post-fire persistence of nesting and roosting (NR) forest in the frequent fire landscapes of the eastern Washington Cascades and the Klamath province and parts of the southern Cascades. Persistent NR forest occurred in topographically sheltered locations and lower slope positions whereas less persistent NR forest occurred in topographically exposed locations. These patterns of persistence reflect historical vegetation patterns when complex interactions among topography, vegetation and weather maintained resilient landscapes. These vegetation patterns may also resemble higher quality heterogenous habitat for spotted owls. Our models can help managers achieve the dual objectives of maintaining viable populations of native species, including those that rely on closed-canopy forests, while also restoring resilient landscapes. In an era of rapid ecological change, we suggest that northern spotted owl conservation and fire-adapted landscape restoration are not opposing goals. By identifying the biophysical settings that promote persistent NR forest and acknowledging the role of heterogeneous landscapes for maintaining resilient ecosystems and high quality spotted owl habitat, the dual objectives of supporting long-term resilience for native species and ecosystems can be achieved.

California, Oregon, Washington↗

The impact of 3D structure on coseismic coastal land-level change and tsunami generation in the Cascadia Subduction Zone

Estimates of coseismic vertical displacements from past and potential future subduction zone earthquakes provide critical constraints on regional seismic and tsunami hazard. Many studies use elastic homogeneous half-space models to calculate vertical displacements from a specified earthquake slip distribution, neglecting complexities of 3D structure in subduction zone settings. Here, we use 3D ground motion simulations of potential Cascadia Subduction Zone earthquakes to investigate the impact of realistic 3D Earth structure on estimated vertical displacements. We show that offshore uplift increases when including 3D structure due to high fault slip in low-rigidity accretionary wedge sediments and coastal subsidence decreases due to higher rigidity material at depth. These larger offshore uplifts cause increased tsunami maximum wave heights, and coastal subsidence at paleoseismic sites decreases on average by ∼17 cm, or ∼60%, which can have important implications for the amplitude and extent of slip in earthquake scenarios developed based on paleoseismic data constraints.

California, Oregon, Washington↗

Fiber-optic sensing for earthquake hazards research, monitoring and early warning

The use of fiber‐optic sensing systems in seismology has exploded in the past decade. Despite an ever‐growing library of ground‐breaking studies, questions remain about the potential of fiber‐optic sensing technologies as tools for advancing if not revolutionizing earthquake‐hazards‐related research, monitoring, and early warning systems. A working group convened to explore these topics; we comprehensively examined the application of fiber optics in various aspects of earthquake hazards, encompassing earthquake source processes, crustal imaging, data archiving, and technological challenges. There is great potential for fiber‐optic systems to advance earthquake monitoring and understanding, but to fully unlock their capabilities requires continued progress in key areas of research and development, including instrument testing and validation, increased dynamic range for applications focused on larger earthquakes, and continued improvement in subsurface and source imaging methods. A key current stumbling block results from the lack of clear data archiving requirements, and we propose an initial strategy that balances data volume requirements with preserving key data for a broad range of future studies. In addition, we demonstrate the potential for fiber‐optic sensing to impact monitoring efforts by documenting the data completeness in a number of long‐term experiments. Finally, we outline the features of a instrument testing facility that would enable progress toward reliable and standardized distributed acoustic sensing data. Overcoming these current obstacles would facilitate progress in fiber‐optic sensing and unlock its potential application to a broad range of earthquake hazard problems.

California, Oregon, Washington↗

Refinements to the attenuated propagation of local earthquake shaking (APPLES) ground-motion-based earthquake early warning algorithm

We refined the Attenuated ProPagation of Local Earthquake Shaking (APPLES) ground-motion-based earthquake early warning (EEW) approach, and directly compare APPLES performance with that of the source-characterization-based U.S. ShakeAlert EEW system for a suite of historical earthquakes in the U.S. West Coast and Japan. APPLES is an extension of the Propagation of Local Undamped Motion (PLUM) algorithm in which observed shaking intensity at seismic stations is used to forward-predict intensity distributions to surrounding areas using an attenuation model derived from an intensity prediction equation. We test new configuration options within APPLES, such as using the second highest estimated ground motion rather than the maximum, to better match median ground-motion observations and reduce alerts for small magnitude earthquakes, both of which are key alerting priorities within ShakeAlert. We evaluate these configurations alongside ShakeAlert by comparing the ground-motion estimation accuracy and available warning times relative to station observations and ShakeMap distributions. Our preferred APPLES configuration produces accurate ground-motion estimates and corresponds better with median observations compared to ShakeAlert’s estimates. This preferred configuration substantially reduces alert issuance for M < 5.0 earthquakes compared to the previous APPLES configuration, and alert-release criteria can further restrict alerts to primarily M ≥ 5.5 earthquakes without requiring magnitude estimation. Prioritizing matching median-observed ground motions may reduce APPLES warning times compared to configurations that were tuned to avoid missed alerts (such as those that use the maximum estimated ground motions), which can lead to shorter warning times compared to ShakeAlert for the same alert threshold. However, station-based warning time assessments demonstrate that APPLES can outperform ShakeAlert for high target thresholds. APPLES is a simple, independent EEW approach that may improve the robustness of EEW for the West Coast of the U.S.

California, Oregon, Washington↗

Earthquake probabilities and hazards in the U.S. Pacific Northwest

Earthquakes and their cascading consequences pose a significant threat to the people, environment, infrastructure, and economy of the U.S. Pacific Northwest. The Pacific Northwest is susceptible to three types of earthquakes: deep (intraslab) earthquakes, subduction zone (megathrust) earthquakes, and shallow crustal earthquakes. For each of these earthquake types, earth scientists can use a variety of methods to estimate the probability of occurrence for future events, which constrains seismic hazard and informs building codes. The timing of past earthquakes indicates that there is an 85-percent chance of a magnitude 6.5 or greater deep earthquake in the Puget Sound region; a 10-15-percent chance of an approximately magnitude 9 earthquake on the Cascadia Subduction Zone; and a 17-percent chance of a magnitude 6.5 or greater crustal fault earthquake in the Puget Sound region in the next 50 years. Individuals and communities can take simple steps to prepare for and reduce the impact of future earthquakes.

California, Oregon, Washington↗

Origin and evolution of mafic volcanism associated with 3 m.y. of andesite production at the Goat Rocks volcanic cluster, southern Washington Cascade Range

More than 3 m.y. of mafic volcanism near the Goat Rocks volcanic cluster in the southern Washington Cascade Range, USA, lends insight into the evolution of basalts and the subarc mantle at a long-lived, major arc volcanic locus. We contribute field observations, 40 Ar/ 39 Ar dates, paleomagnetic directions, and bulk rock and mineral compositions to characterize nine mafic units that erupted in association with the Goat Rocks volcanic cluster. The time frame of mafic volcanism, ca. 3.6 Ma to 60 ka, encompasses the lifespan of the central volcanic cluster (3.1 Ma to 115 ka), with a lull from ca. 2.7 Ma to 1.4 Ma. A climactic period of voluminous mafic activity and far-traveled lava flows, including construction of the Hogback Mountain shield volcano, coincided with voluminous andesite eruptions from the central volcanic cluster. The basaltic rocks in the Goat Rocks area are calc-alkaline to barely tholeiitic and have high field strength element depletion relative to large-ion lithophile elements characteristic of calc-alkaline basalts (CAB) of the Cascade volcanic arc. Unlike at neighboring andesitic volcanic centers (Mounts Adams, St. Helens, and Rainier), no other mafic end members such as high-aluminum olivine tholeiite (HAOT) or intraplate-type basalt (IPB) are present at or near the Goat Rocks volcanic cluster, although some of the calc-alkaline basalts in this study have IPB-like affinities. The Goat Rocks mafic units exhibit two main temporal trends in composition: (1) the most primitive basalts erupted earlier, compared to less primitive and more evolved compositions later, and (2) high field strength element concentrations are higher in the younger basalt units relative to the oldest two. In contrast to these temporal trends, the mafic units define two compositional groups that recur through time, a low-Sr and a high-Sr group, each with distinct trace element and Sr and Nd isotope ratios. Although radiogenic isotope ratios are generally aligned with High Cascades CAB and HAOT, some extend toward IPB of Mount Adams and Simcoe Mountains volcanic field. Olivine-dominated crystal fractionation at shallow pressure from a small range of parent magma compositions accounts for much of the variation among the basalts and basaltic andesites. A high-pressure fractionation model is plausible for only one of the youngest basalt units (basalt of Walupt Lake volcano). Mafic recharge and crustal assimilation accounts for the incompatible-element enriched composition of basaltic andesites erupted during construction of the largest andesitic centers, further supporting sustained basalt mass flux and thermal energy driving andesite genesis. We model the most primitive members of the Goat Rocks mafic units as partial melts of successively less depleted mantle in time. Variable degrees of fluxing with fluids and melts from subduction explain the distinction between high-Sr and low-Sr groups. We propose that mantle metasomatism by ancestral subduction and fluid-flux melting is heterogeneously distributed through the local subarc mantle and played a greater role in the genesis of the high-Sr basalt group. The limited range of primitive basalt types around the Goat Rocks volcanic cluster contrasts with the much greater diversity of basalts throughout the southern Washington to northern Oregon Cascade arc. On the other hand, the central volcanic cluster encompasses nearly the entire diversity observed at neighboring composite volcanoes. In the case of the Goat Rocks area at least, and perhaps attributable to the entire region, this means that the genesis of diverse intermediate magmas is independent from and does not require vastly different parental basalt compositions.

California, Oregon, Washington↗

Rupture continuity through intermittent pauses in Cascadia slow slip events

Cascadia slow slip events (SSEs) are often envisioned as smooth, continuous ruptures, progressively activating tremor asperities as they propagate. Macroscopically, geodetic inversions and spatiotemporal maps of tremor epicenters show steady, uniform migration. In detail tremor is more chaotic and discontinuous. Larger long-term SSEs observed in daily geodetic solutions are inferred to exhibit intermittent pauses that reflect temporary re-locking of the fault, but this temporal resolution limits tests for similar re-locking on shorter timescales. We use temporal measurements of the areal growth and radiated energy of tremor clusters to investigate SSE intermittence. We find that ruptures mirror tremor pauses. Areal growth rate, however, does not reset, and removing the pauses results in smoother and more similar growth measurements among all SSEs. The rupture similarity occurs regardless of size or location and hints at an underlying uniformity and lack of predeterminism in eventual SSE size. Epicentral uncertainty precludes quantifying early rupture stages, but for larger events areal growth follows a power-law and slows with increasing size. Temporal correlations in tremor energy with inferred SSE propagation velocities and tremor rates suggest its use as a proxy for slip velocity. We find that tremor energy is tidally modulated at daily and sub-daily frequencies, and this modulation is continuous through pauses, suggesting a memory of slip state is sustained through them. We argue these pauses reflect unsteady propagation of the slip front, marked by rapid re- and un-locking, and excluding them removes rupture complexity to reveal a diffusive-like slip process and underlying universality in growth.

California, Oregon, Washington↗

Organic matter availability and the production of methane-derived dissolved organic carbon at methane seeps

Methane seeps located along continental margins and slopes export sedimentary methane and dissolved organic carbon (DOC) into the ocean. The flux of these reduced carbon molecules from the seafloor into the ocean impacts ocean chemistry and supports deep-sea life. While significant effort has been made to understand how the anaerobic oxidation of methane (AOM) regulates the release of methane from the seafloor, little is known about the production of DOC in association with AOM or its flux and fate in the ocean. We hypothesize a mechanism for methane incorporation into DOC at seeps and investigate the relationship between sediment total organic carbon (TOC) availability and the incorporation of methane-derived carbon into DOC at four methane seep regions along the Cascadia margin, with a range of microbial and thermogenic methane sources. At sites with <2.0 wt.% TOC (Hydrate Ridge and Bullseye Vent), up to 60%–80% of carbon in DOC is methane-carbon, much more than sites with >2.0 wt.% TOC (Astoria Canyon and Barkley Canyon). We attribute the greater methane contribution at the more TOC-limited sites to a greater role of AOM in the carbon cycle, whereas at the organic matter-rich sites, microbial competition for sulfate as an electron acceptor for organic matter decomposition limits AOM and hence the transfer of carbon from methane to DOC. We estimate that the global diffusive flux of methane-derived DOC from the seafloor is 0.07–10.1 Tg C/yr, contributing to the stock of DOC present in the deep ocean and/or fueling the deep-sea microbial loop.

California, Oregon, Washington↗

Increased flood exposure in the Pacific Northwest following earthquake-driven subsidence and sea-level rise

Climate-driven sea-level rise is increasing the frequency of coastal flooding worldwide, exacerbated locally by factors like land subsidence from groundwater and resource extraction. However, a process rarely considered in future sea-level rise scenarios is sudden (over minutes) land subsidence associated with great (>M8) earthquakes, which can exceed 1 m. Along the Washington, Oregon, and northern California coasts, the next great Cascadia subduction zone earthquake could cause up to 2 m of sudden coastal subsidence, dramatically raising sea level, expanding floodplains, and increasing the flood risk to local communities. Here, we quantify the potential expansion of the 1% floodplain (i.e., the area with an annual flood risk of 1%) under low (~0.5 m), medium (~1 m), and high (~2 m) earthquake-driven subsidence scenarios at 24 Cascadia estuaries. If a great earthquake occurred today, floodplains could expand by 90 km 2 (low), 160 km 2 (medium), or 300 km 2 (high subsidence), more than doubling the flooding exposure of residents, structures, and roads under the high subsidence scenario. By 2100, when climate-driven sea-level rise will compound the hazard, a great earthquake could expand floodplains by 170 km 2 (low), 240 km 2 (medium), or 370 km 2 (high subsidence), more than tripling the flooding exposure of residents, structures, and roads under the high subsidence scenario compared to the 2023 floodplain. Our findings can support decision-makers and coastal communities along the Cascadia subduction zone as they prepare for compound hazards from the earthquake cycle and climate-driven sea-level rise and provide critical insights for tectonically active coastlines globally.

California, Oregon, Washington↗

3D viscoelastic models of slip-deficit rate along the Cascadia subduction zone

Interseismic deformation in the Pacific Northwest is constrained by the horizontal crustal velocity field derived from the Global Positioning System (GPS) in addition to vertical rates derived from GPS, leveling, and tide gauge measurements. Such measurements were folded in to deformation models of fault slip rates as part of the 2023 National Seismic Hazard Model (NSHM) update. Here I build upon one of the contributing models, the viscoelastic earthquake-cycle model of Pollitz [2022]. This model permits inclusion of effects of time-dependent viscoelastic relaxation within earthquake cycles (i.e., ‘ghost transients’) and laterally variable elastic and/or ductile material properties. I lever-age these capabilities to incorporate the Cascadia megathrust into Western U.S.-wide deformation models in which crustal fault slip rates are estimated simultaneously with slip deficit rates along the interplate boundary between the descending Juan de Fuca plate and North American plate. This effort includes construction of a margin-wide model of viscoelastic structure founded on the Slab 2.0 model and probes different models of the ductile properties of the surrounding oceanic asthenosphere, continental lower crust, and mantle asthenosphere. This results in new estimates of the distribution of slip deficit rate along the ∼ 1000 km long margin, highlights the importance of correcting for glacial-isostatic adjustment effects, and permits assessment of sensitivity of results to assumed ductile properties.

California, Oregon, Washington↗

Long-lived partial melt beneath Cascade volcanoes

Quantitative estimates of magma storage are fundamental to evaluating volcanic dynamics and hazards. Yet our understanding of subvolcanic magmatic plumbing systems and their variability remains limited. There is ongoing debate regarding the ephemerality of shallow magma storage and its volume relative to eruptive output, and so whether an upper-crustal magma body could be a sign of imminent eruption. Here we present seismic imaging of subvolcanic magmatic systems along the Cascade Range arc from systematically modelling the three-dimensional scattered wavefield of teleseismic body waves. This reveals compelling evidence of low-seismic-velocity bodies indicative of partial melt between 5 and 15 km depth beneath most Cascade Range volcanoes. The magma reservoirs beneath these volcanoes vary in depth, size and complexity, but upper-crustal magma bodies are widespread, irrespective of the eruptive flux or time since the last eruption of the associated volcano. This indicates that large volumes of melts can persist at shallow depth throughout eruption cycles beneath large volcanoes.

California, Oregon, Washington↗

Strong shaking from past Cascadia Subduction Zone earthquakes encoded in coastal landforms

Strong earthquakes along subduction zones are often devastating events, but sparse records along some tectonic margins limit our understanding of seismic hazards. Constraining shaking intensities is critical, especially in subduction zones with infrequent but large-magnitude earthquakes like the Cascadia Subduction Zone (CSZ), where the lack of recorded ground motions has led to uncertainty in the severity and potential impacts of future earthquakes. Here we fill this observational gap with a novel inventory of quantitative estimates of past shaking intensities from geotechnical modeling of coastal landforms. One hundred fifty-four deep-seated landslides and 65 fragile geologic features constrain minimum and maximum peak ground accelerations, respectively. These estimates are broadly consistent with model predictions of M9 ruptures, suggesting strong shaking of 0.4–0.8 g during past CSZ earthquakes. Local discrepancies between our geologic shaking constraints and earthquake simulations may inform past rupture behavior, leading to better predictions of shaking intensity for future earthquakes.

California, Oregon, Washington↗

Factors structuring macrofaunal communities of hydrocarbon seeps along the Cascadia margin

Cold seeps, fueled by hydrocarbon-based chemosynthesis, support unique benthic communities that can vary across small spatial scales influenced by local geochemistry. We examined the community structure and function of macrofauna in cold seeps along the Cascadia margin on the edge of gas hydrate stability—a dynamic environment that may fluctuate as seafloor temperatures change. These communities were assessed in relation to their sediment environment to provide baseline data for these previously uninvestigated seeps and decipher the most significant variables structuring them. Specific environmental drivers investigated include sediment organic carbon, total nitrogen, C:N ratios, stable isotopes (δ 13 C and δ 15 N), redox potential, grain size, and porefluid chemistry. Macrofaunal community composition varied across locations separating into twelve distinct clusters. Several co-located cores clustered separately, indicating high heterogeneity in these communities at small spatial scales. Significant drivers of macrofauna communities included clay and sand content, C:N, δ 15 N values, organic carbon content, and depth. Functional trait composition was influenced by similar drivers, including depth, δ 15 N values, C:N, and sand content. Our results indicate similar environmental variables structure macrofaunal community composition and function across small- and large-scale gradients, contributing to our understanding of heterogeneity in local and regional factors that shape seep macrofaunal communities. This is the first investigation of macrofauna at these recently discovered seeps and provides baseline data on the hydrate stability zone for future comparisons, advancing our knowledge of broad-scale trends in seep macrofauna ecology.

California, Oregon, Washington↗