Geology Reports⌕ Search

Geology topics

Research about Washington

Source-linked reports with geographic coverage including Washington.

At least 253 records · Page 14Linked to original sources

Hybridization between historically allopatric Chinook salmon populations in the White Salmon River, WA

Chinook Salmon spawning in the White Salmon River consist of members of three historically distinct populations: spring Chinook Salmon, Tule fall Chinook Salmon and Upriver Bright (URB) fall Chinook Salmon. Previous work examined juveniles captured in 2006-2008 and reported hybridization between introduced URBs, and the native threatened Tules. Recent increases in nearby hatchery URB release numbers raised the question of whether hybridization rates were increasing. We estimated hybrid frequencies among juveniles collected in the lower White Salmon River between 2016 and 2019. We also evaluated the frequencies at which non-target fish and hybrids were incorporated into the broodstocks of adjacent hatcheries. We observed that frequencies of hybrids in juvenile samples from the White Salmon River were greater in 2017-2019 (17-32%) than they had been in 2006-2008 (4-15%), but that a few (2/9) comparisons exhibited overlapping confidence intervals, suggesting that the rate has increased over time, but also that more sampling is needed to understand the importance of year-to-year variation. Further, differences in the habitat following dam removal and in the sampling sites complicated interpretation of our results. Examination of broodstocks of nearby hatcheries revealed low rates (< 0.5%) of incorporation of non-target populations and higher rates (< 9.0%) of incorporation of hybrids into those broodstocks. The relative compositions of all hatchery and natural-origin collections were similar: most individuals were one of the two parental stocks, a very small fraction were F1 hybrids, and a larger minority fraction were back-crosses. This pattern, in the context of hybridization which we know has been happening for several generations, is consistent with a hypothesis of selection against hybrids in which F1 hybrids are less fit than backcross hybrids.

Washington↗

Freezing in a warming climate: Marked declines of a subnivean hibernator after a snow drought

Recent snow droughts associated with unusually warm winters are predicted to increase in frequency and affect species dependent upon snowpack for winter survival. Changes in populations of some cold‐adapted species have been attributed to heat stress or indirect effects on habitat from unusually warm summers, but little is known about the importance of winter weather to population dynamics and how responses to snow drought vary among sympatric species. We evaluated changes in abundance of hoary marmots ( Marmota caligata ) over a period that included a year of record‐low snowpack to identify mechanisms associated with weather and snowpack. To consider interspecies comparisons, our analysis used the same a priori model set as a concurrent study that evaluated responses of American pikas ( Ochotona princeps ) to weather and snowpack in the same study area of North Cascades National Park, Washington, USA. We hypothesized that marmot abundance reflected mechanisms related to heat stress, cold stress, cold exposure without an insulating snowpack, snowpack duration, atmospheric moisture, growing‐season precipitation, or select combinations of these mechanisms. Changes in marmot abundances included a 74% decline from 2007 to 2016 and were best explained by an interaction of chronic dryness with exposure to acute cold without snowpack in winter. Physiological stress during hibernation from exposure to cold, dry air appeared to be the most likely mechanism of change in marmot abundance. Alternative mechanisms associated with changes to winter weather, including early emergence from hibernation or altered vegetation dynamics, had less support. A post hoc assessment of vegetative phenology and productivity did not support vegetation dynamics as a primary driver of marmot abundance across years. Although marmot and pika abundances were explained by strikingly similar models over periods of many years, details of the mechanisms involved likely differ between species because pika abundances increased in areas where marmots declined. Such differences may lead to diverging geographic distributions of these species as global change continues.

Washington↗

A probabilistic assessment of tephra-fall hazards at Hanford, Washington, from a future eruption of Mount St. Helens

Hanford, Washington (USA) is the construction site of a multi-billion-dollar high-level nuclear waste treatment facility. This site lies 200 kilometers (km) east of Mount St. Helens (MSH), the most active volcano in the contiguous United States. Tephra from a future MSH eruption could pose a hazard to the air intake and filtration systems at this plant. In this report, we present a probabilistic estimate of the amount of tephra that could fall, and the concentrations of airborne ash that could occur at the Hanford Site during a future eruption. Mount St. Helens has produced four large explosive eruptions in approximately the past 500 years, suggesting that its annual probability of eruption ( P 1 ) is roughly 4/500=0.008. Assuming that a large eruption occurs, we calculate the probability ( P 3|1 ) of a given fall deposit thickness or airborne concentration at Hanford by running about 10,000 simulations of ash-producing eruptions using the atmospheric transport model Ash3d. In each simulation, we calculate the pattern of tephra dispersal, deposit thickness at Hanford, and airborne ash concentration at ground level. As input for each simulation, we choose meteorological conditions from a randomly chosen time in the historical record between 1980 and 2010, using data from the European Centre for Medium-Range Weather Forecasting (ECMWF) Reanalysis (ERA) Interim model. The volume (dense-rock equivalent) of each simulated eruption is randomly chosen from a uniform probability distribution on a log scale from the range of magma volumes (0.008–2.3 cubic kilometers [km 3 ]) estimated for late Holocene eruptions at MSH. Plume heights and durations of each eruption are chosen using empirical correlations between volume, height, and eruption rate, which account for the fact that larger eruptions have higher plumes and last longer. We construct summary tables of final deposit thickness ( T ), maximum ground-level airborne concentration ( C max ), and average ground-level airborne concentration ( C avg ) during tephra-fall for each run. Each table is sorted and ranked by decreasing value of T , C max , or C avg . Conditional probabilities ( P 3|1 ) are derived by dividing rank by n+1, where n is the total number of successful runs. For example, a deposit thickness of 5.10 centimeters (cm) from run 446 is ranked 123 of 9,785 successful runs, yielding P 3|1 =123/9,786=0.01257. Its annual probability is P = P 1 · P 3|1 =0.008×0.01257=0.000101. By interpolation, the deposit thickness ( T 10k ) having an annual probability of 1 in 10,000 ( P = 0.0001) is 5.11 cm. Analogous concentration values are C max,10k =3,819 and C avg,10k =1,513 milligrams per cubic meter (mg/m 3 ), respectively. Independent calculations using the known mass accumulation rate of the deposit (=0.001–0.006 kilograms per square meter per second [kg/m 2 /s]), aggregate fall velocities ( u =0.3–0.8 meters per second [m/s]), and the simple formula , yield similar results, although highly variable fall velocities add significant uncertainty. This formula implies that deposit accumulation rates of millimeters (mm) to greater than 1 cm per hour, which are not uncommon during heavy ash fall, are associated with airborne concentrations of 10 2 –10 3 milligrams per cubic meter (mg/m 3 ). These concentrations are much higher than published measurements (10 -3 –10 1 mg/m 3 ), which record only suspended particles sampled in sheltered areas. During heavy ashfall, most fine ash falls as aggregates. Whether such aggregates will be ingested into air ducts will depend on the aggregate size and fall rate, the fragility of the aggregates, the air duct geometry, intake velocity, and other factors.

Washington↗

Improving paleoseismic earthquake magnitude estimates with rupture length information: Application to the Puget Lowland, Washington State, U.S.A.

Both earthquake displacement and rupture length correlate with magnitude, and therefore observations of each from past earthquakes can be used to estimate the magnitude of those earthquakes in the absence of instrumental records. We extend the Bayesian inversion method of Biasi and Weldon (2006), which estimates paleoearthquake magnitude from displacement observations, to incorporate both rupture length and surface displacement measurements into the magnitude inversion. We then use this method on 27 late Pleistocene to Holocene paleoearthquakes in the Puget Lowland region of Washington. Observations of (typically vertical) fault separation per event range from 0.6 to 7 m, implying net displacement per event of up to 10 ± 4 m for the largest event. Rupture lengths are estimated to vary between the smallest contiguous mapped scarps to the full extent of the faults mapped from geology and geophysical observations. Although a few of the ruptures may be longer than 150 km, the ruptures have a median of 53 km, indicating that earthquakes in the Puget Lowland have relatively high displacement to length ratios. By considering both datasets, we find that all events were between M 6.3 and 7.5, generally consistent with the expected seismicity from the USGS National Seismic Hazard Map for the region. The simultaneous use of both length and displacement data in the magnitude inversion decreases both the estimated earthquake magnitudes and the uncertainty. The magnitude reduction in particular is due to the relatively short rupture lengths possible for Puget Lowland faults. This implies a decrease in the seismic hazard (relative to a displacement-only assessment) to a highly populated and rapidly urbanizing region.

Washington↗

From forests to fish: Mercury in mountain lake food webs influenced by factors at multiple scales

Mountain lakes, while seemingly pristine, have been subjected to historical fish stocking practices and exposure to atmospherically deposited contaminants like mercury. Mercury bioaccumulation in these ecosystems varies widely due to strong environmental gradients, and there are complex, hierarchical factors that affect mercury transport and loading, methylmercury production, and food web biomagnification. We sought to assess how representative variables associated with watershed, lake, and food web‐scale processes—specifically, catchment tree cover, lake benthic primary production, and fish diet, respectively—are associated with mercury concentrations in mountain lake fish. Mean fish mercury concentrations varied threefold between lakes, with nearshore tree cover and fish diet accounting for the most variance in fish mercury. Tree cover was likely positively correlated to fish Hg due to its contributions to local deposition and its effect on lake biogeochemistry. Fish with benthic diets tended to have higher mercury concentrations, illustrating that food web processes are an important consideration when investigating drivers of contaminant bioaccumulation. Our results suggest that both landscape and ecological factors are determinants of fish mercury bioaccumulation, and thus variables at multiple scales should be considered when managing mountain lake food webs for mercury exposure risk.

Washington↗

Fish assemblages in eelgrass beds of Bellingham Bay, Washington, Northern Puget Sound, 2019

Puget Sound is a critical part of the Pacific Northwest, both culturally and economically. Eelgrass beds are an important feature of Puget Sound and are known to influence fish assemblages. As part of a larger site-characterization effort, and to gain a better understanding of the fish assemblages in Bellingham Bay, Washington, four eelgrass beds ( Zostera marina ) along the shoreline were surveyed. Fish were captured from 24 through 26 September 2019 by using three beach-seine hauls per eelgrass bed. In total, 12 hauls yielded 2,135 fish that comprised 20 species from 14 families. Shiner perch ( Cymatogaster aggregata ) accounted for 52 percent of the total catch. The other common species included three-spine stickleback ( Gasterosteus aculeatus ), bay pipefish ( Syngnathus leptorhynchus ), saddleback gunnel ( Pholis ornata ), Pacific staghorn sculpin ( Leptocottus armatus ), and Pacific sand lance ( Ammodytes personatus ). Total catch and species richness were highest at the two locations closest to the urban center of Bellingham; however, species diversity and evenness were highest at the two eelgrass beds farthest from the city center. Descriptions of fish assemblages in eelgrass beds are expected to be useful in the development of future process-based investigations by study partners and will focus on the movements of sediments and contaminants and their influence on biota.

Washington↗

Post-glacial Mw 7.0-7.5 earthquakes on the North Olympic fault zone, Washington

Holocene crustal faulting in the northern Olympic Peninsula of Washington State manifests in a zone of west‐northwest‐striking crustal faults herein named the North Olympic fault zone, which extends for &#x223C; 80 &#x2009;&#x2009; km "> ∼ 80 km ∼80 km along strike and includes the Lake Creek–Boundary Creek fault to the east and the Sadie Creek fault and newly discovered scarps to the west. This study focuses on the Sadie Creek fault, which extends for &gt; 14 &#x2009;&#x2009; km "> > 14 km >14 km west‐northwest from Lake Crescent. Airborne light detection and ranging (lidar) imagery reveals the trace of the Sadie Creek fault and offset postglacial landforms showing a history of Holocene surface‐rupturing earthquakes dominated by dextral displacement along a steeply dipping fault zone. Paleoseismic trenches at two sites on the Sadie Creek fault reveal till and outwash overlain by progressively buried forest and wetland soils developed on scarp‐derived colluvial wedges. Trench exposures of complex faulting with subhorizontal slickenlines indicate dextral displacement with lesser dip slip. Correlation of broadly constrained time intervals for earthquakes at the Sadie Creek sites and those to the east along the Lake Creek–Boundary Creek fault is consistent with rupture of much of the length of the North Olympic fault zone three to four times: at about 11, 7, 3, and 1 ka, with a shorter rupture at about 8.5 ka. Dated ruptures from trenches only partially coincide with coseismic landslides and megaturbidites in Lake Crescent, indicating that some earthquakes did not trigger megaturbidites, and some turbidites were unrelated to local fault rupture. Landform mapping suggests single‐event dextral displacement of 4 &#xB1; 1 &#x2009;&#x2009; m "> 4 ± 1 m 4±1 m on the Sadie Creek fault. Inferred maximum rupture length and single‐event slip imply earthquake magnitudes M w "> M w Mw 7.0–7.5. Dextral slip rates of 1.3 &#x2013; 2.3 &#x2009;&#x2009; mm / yr "> 1.3 – 2.3 mm / yr 1.3–2.3 mm/yr and the &#x223C; 11 , 000 &#x2009;&#x2009; yr "> ∼ 11 , 000 yr ∼11,000 yr slip history suggest that the North Olympic fault zone is a prominent contributor to permanent strain in the northern Cascadia fore‐arc.

Washington↗

Variable-density groundwater flow and contaminant transport, Operable Unit 1, Naval Base Kitsap, Keyport, Washington

Chlorinated volatile organic compounds (CVOCs) have migrated to groundwater beneath a former 9-acre landfill at Operable Unit 1 (OU-1) on Naval Base Kitsap, which was active from the 1930s through 1973 on the Keyport Peninsula, in Kitsap County, Washington. Biodegradation of CVOCs at OU-1 limits the mass of dissolved-phase CVOCs in groundwater that discharges to surface water, but contaminant concentrations up to 630 milligrams per liter persist in localized areas, likely from the dissolution of residual, non-aqueous phase liquids. Variable-density groundwater-flow and contaminant-transport models were developed using the SEAWAT-Version 4 computer program to simulate the direction and rate of groundwater flow in a 5.9 square-mile (mi 2 ) - area surrounding the Keyport Peninsula, to estimate the CVOC mass in groundwater and the rate of mass loading, and to assess possible remedial activities at OU-1. The study area is underlain by Quaternary deposits consisting of alternating glacial and interglacial sediments ranging from 500 to 1,500 feet (ft) thick. A hydrogeologic model delineated a sequence of 10 units including a relatively thin package (less than 100 ft) of recent sediments (Vashon Stade and younger) beneath the Keyport Peninsula that are underlain by the much thicker (more than 300 ft) Clover Park Aquitard, which overlies a confined, sea-level aquifer.

Washington↗

Landslide monitoring and runout hazard assessment by integrating multi-source remote sensing and numerical models: An application to the Gold Basin landslide complex, northern Washington

The landslide complex at Gold Basin, Washington, has been drawing considerable attention after a catastrophic runout of the nearby landslide at Oso, Washington, in 2014. To evaluate potential threats of the Gold Basin landslide to the campground down the slope, remote sensing and numerical modeling were integrated to monitor recent landslide activity and simulate hypothetical runout scenarios. Bare-earth LiDAR DEM (digital elevation model) differencing, InSAR (Interferometric Synthetic Aperture Radar), and offset tracking of SAR images reveal that localized collapses at the headscarps have been the primary type of landslide activity at Gold Basin from 2005 to 2019, and currently no signs indicative of movement of a large centralized block or a deep-seated main body were detected. The maximum horizontal deformation rate is 5 m/year occurring primarily from headscarp recession of the middle lobe, and the annual landsliding volume of the whole landslide complex averages 1.03 × 10 5 m 3 . From three-dimensional limit equilibrium analysis of generalized terrace structures, the maximum landslide volume is estimated as 2.0 × 10 6 m 3 . Simulations of hypothetical runout scenarios were carried out using the depth-averaged two-phase model D-claw with above-obtained landslide geometry constraints. The simulation results demonstrate that debris flows with volume less than 10 5 m 3 only pose limited threats to the campground, while volumes over 10 6 m 3 could cause severe damages. Consequently, the estimated maximum landslide volume of 2.0 × 10 6 m 3 suggests a potential risk to the campground nearby. Adaption of our methodology could prove useful for evaluating other similar landslides globally for hazards prevention and mitigation.

Washington↗

Shifting food web structure during dam removal—Disturbance and recovery during a major restoration action

We measured food availability and diet composition of juvenile salmonids over multiple years and seasons before and during the world’s largest dam removal on the Elwha River, Washington State. We conducted these measurements over three sediment-impacted sections (the estuary and two sections of the river downstream of each dam) and compared these to data collected from mainstem tributaries not directly affected by the massive amount of sediment released from the reservoirs. We found that sediment impacts from dam removal significantly reduced invertebrate prey availability, but juvenile salmon adjusted their foraging so that the amount of energy in diets was similar before and during dam removal. This general pattern was seen in both river and estuary habitats, although the mechanisms driving the change and the response differed between habitats. In the estuary, the dietary shifts were related to changes in invertebrate assemblages following a hydrological transition from brackish to freshwater caused by sediment deposition at the river’s mouth. The loss of brackish invertebrate species caused fish to increase piscivory and rely on new prey sources such as plankton. In the river, energy provided to fish by Ephemeroptera, Plecoptera, and Trichoptera taxa before dam removal was replaced first by terrestrial invertebrates, and then by sediment-tolerant taxa such as Chironomidae. The results of our study are consistent with many others that have shown sharp declines in invertebrate density during dam removal. Our study further shows how those changes can move through the food web and affect fish diet composition, selectivity, and energy availability. As we move further along the dam removal response trajectory, we hypothesize that food web complexity will continue to increase as annual sediment load now approaches natural background levels, anadromous fish have recolonized the majority of the watershed between and above the former dams, and revegetation and microhabitats continue to develop in the estuary.

Washington↗

Disentangling the effects of multiple fires on spatially interspersed sagebrush (Artemisia spp.) communities

Questions Relative to a landscape with a mosaic of two sagebrush community types and increasing fire frequency, we asked: (a) do vegetation characteristics vary significantly with number of times burned for each sagebrush community; (b) how do vegetation responses to different fire frequencies compare between the two sagebrush communities? Location Columbia Plateau Ecoregion, Washington, USA. Methods We sampled vegetation across a landscape that burned three times over a 10‐year period in two sagebrush community types that are interspersed on unique land forms: big sagebrush ( Artemisia tridentata ) communities that occur on small “mounds” and scabland sagebrush ( A. rigida ) communities that occur on surrounding “flats.” Spatially overlapping fires permitted a balanced sampling design to assess unburned and once‐, twice‐, and thrice‐burned locations for each land form/community type. We utilized a suite of statistical analyses to determine differences among plant functional groups and biomass among unburned/burned strata by land form and compared results between land forms. Results Big sagebrush and scabland sagebrush communities responded uniquely to multiple fires, due to different fuel loadings, fire severities, succession and invasion dynamics. Big sagebrush experienced nearly complete shrub loss and conversion from exotic‐invaded shrubland to exotic annual grassland after only one fire. In contrast, scabland sagebrush retained a minor shrub component and higher relative cover of native herbaceous species, even after three fires. Both communities retained cover of native perennial grasses, including shallow‐ and deep‐rooted species, likely reflecting decreasing fire intensity with number of times burned. Conclusions Despite different community‐level responses, increasing fire frequency is transforming the entire landscape to a non‐native/native grassland mix. Quantifying unique ecosystem responses to altered wildfire regimes is critical to understanding the relative resilience of communities to disturbance and their resistance to exotic species invasion (and community type conversion). Management actions may help to maintain spatial heterogeneity of ecosystems and fire‐tolerant native species.

Washington↗

Coarse sediment dynamics in a large glaciated river system: Holocene history and storage dynamics dictate contemporary climate sensitivity

The gravel-bedded White River drains a 1279 km 2 basin in Washington State, with lowlands sculpted by continental glaciation and headwaters on an actively glaciated stratovolcano. Chronic aggradation along an alluvial fan near the river’s mouth has progressively reduced flood conveyance. In order to better understand how forecasted climate change may influence coarse sediment delivery and aggradation rates in this lowland depositional setting, we assessed the contemporary delivery and routing of coarse sediment through the watershed; this assessment was based on a rich set of topographic, sedimentologic, and hydrologic data from the past century, with a focus on repeat high-resolution topographic surveys from the past decade. We found that most of the lower river’s contemporary bed-load flux originates from persistent erosion of alluvial deposits in the lower watershed. This erosion is a response to a drop in local base level caused by a major avulsion across the fan in 1906 and then augmented by subsequent dredging. The 1906 avulsion and modern disequilibrium valley profiles reflect landscape conditioning by continental glaciation and a massive mid-Holocene lahar. In the proglacial headwaters, infrequent large sediment pulses have accomplished most of the observed coarse sediment export, with exported material blanketing downstream valley floors; during typical floods, transported bed material is largely sourced from erosion of these valley floor deposits. Throughout the watershed, we observe decadal-scale coarse sediment dynamics strongly related to the filling or emptying of valley-scale sediment storage over 10 2 −10 4 yr time scales, often in response to major disturbances that either emplace large deposits or influence their redistribution. Paraglacial responses in large watersheds are suggested to be inherently complicated and punctuated as a result of internal landform interactions and stochastic/threshold-dependent events. We argue, in combination, that Holocene disturbance, storage dynamics, and human flow modification make coarse sediment fluxes in the lower White River relatively insensitive to decadal climate variability. Results highlight the degree to which river sensitivity to contemporary disturbance, climatic or otherwise, may be contingent on local and idiosyncratic watershed histories, underscoring the need to unpack those histories while demonstrating the utility of watershed-scale high-resolution topography toward that end.

Washington↗

Contaminants in fish and shellfish in the Stillaguamish River and Port Susan marine areas, Washington

The greater Port Susan area of Central Puget Sound, Washington, is home to some of the Stillaguamish Tribe’s fishing, hunting, and gathering areas since time immemorial. It is also a popular sport and commercial fishing area for the public. Large shellfish beds lie in the Port Susan and Stillaguamish estuary and several Pacific salmon species return to the Stillaguamish River and Tulalip fishery every year. Clams and salmon are a local and consumable resource for Tribal members and the public. This review largely confirms existing recommendations from the Washington State Department of Health regarding clam and salmon human consumption advisories.

Washington↗

Geomorphic map of western Whatcom County, Washington

Western Whatcom County has a rich history of glaciation, sea-level change, fluvial erosion and deposition, landsliding, nearby volcanic activity, and human landscape modification. This lidar-derived geomorphic map interprets this history from the form and position of the Earth’s surface. The geomorphic record is broken into nine phases, beginning with the peak of the Vashon stade of the Fraser glaciation of Armstrong and others (1965) (phase 1), shortly after 16,000 years ago. The Cordilleran ice sheet was ≥1.6 km thick in the Bellingham area. Glacial lineations on high ground demonstrate that ice flow was from north to south. Storage of water in ice sheets at this time resulted in global sea level ~120 m lower than at present. The weight of the ice sheet depressed the land so that local relative sea level was at least 150 m higher than at present. As the ice sheet melted and thinned, it floated, broke up, and was replaced by salt water. The margin of the ice sheet—or at least its grounding line—retreated to the northeast of the map area during or before phase 2. Marine deposition, currents, and waves smoothed earlier-formed surfaces in the western part of the map area. Global sea level rose (because of melting of continental ice sheets), but the Fraser Lowland rose even faster (due to glacio-isostatic rebound following the loss of ice-sheet load), and thus local relative sea level fell. The Cordilleran ice sheet readvanced during the Sumas stade of Armstrong and others (1965). Oldest Sumas moraines formed when relative sea level at Bellingham was ~55 m (phase 3). Younger moraines formed when relative sea level at Bellingham was ~25 m (phase 4). The amount of Sumas ice retreat and readvance between these times is unknown. Younger Sumas events are marked by local moraines, progressive isostatic rebound and lowering of relative sea level, and changes in the flow of ice-marginal water. During phase 5, the southeast margin of the ice sheet advanced, perhaps because capture of ice-marginal drainage by the Samish River (east and south of the map area) meant the ice sheet was no longer trimmed by high-discharge flow along Squalicum channel. Farther west and north, the ice margin retreated between phases 4 and 5. Phases 6 through 9 may mark stillstands during further ice retreat. There were glacial outburst floods (jökulhlaups) during phases 7 and 8, and perhaps during phase 5. When Sumas ice left the area, perhaps about 11,500 years ago, the Nooksack River appears to have discharged northeast through Sumas Valley to the Fraser River. Details of the switch to its modern course are speculative, but archaeological and sediment-supply arguments suggest that the modern Nooksack River delta south of Ferndale formed within the past 5,000 years. The foothills of the North Cascades are decorated with abundant post-glacial deep-seated landslides. Anomalously high late Holocene beaches are found at Birch Bay, Neptune Beach, perhaps at Maple Beach on the east side of Point Roberts, and perhaps at the northwest corner of the Lummi Peninsula. These beaches may have been uplifted by earthquakes that did not rupture the surface. The low-relief landscape shaped by the Cordilleran ice sheet, along with fluvial infilling of low areas, resulted in abundant wetland, at least 70 percent of which has been diked and (or) drained to control flooding and facilitate farming.

Washington↗

Latest Pleistocene to present geology of Mount Baker Volcano, northern Cascade Range, Washington

Mount Baker is the prominent andesitic stratocone that forms the youngest volcanic center in the Mount Baker volcanic field. Its heavily glaciated cone, rising to 3,286 meters, is an international landmark, dominating the skyline of Vancouver, British Columbia, even though the volcano is located 25 kilometers south of the international border. Mount Baker caught the attention of scientists and the public alike in 1975–76 during a period of increased steaming, thermal output, and near-vent lithic tephra falls. Although a magmatic eruption did not ensue, it awoke the populace to the possibility of renewed volcanic activity in the Cascade Range (the first since the 1914–17 eruptions of Lassen Peak, Calif.)—a possibility fulfilled just five short years later with the 1980 eruption of Mount St. Helens in southwest Washington. The 1980 Mount St. Helens eruption, with its dramatic edifice collapse, extraordinary pyroclastic density current, and catastrophic lahars, invigorated the scientific community into studying these then little-known processes. It also highlighted the need to better understand eruptive histories at other Cascade Range volcanoes in order to prepare for future eruptions. The 1975 unrest also spawned one of the earliest volcano hazard assessments in the Cascade Range, which recognized the rich history of postglacial events at Mount Baker and identified the risk posed by volcanic mudflows, or lahars. The focus of this study is to more fully describe the late-glacial to present surficial geology, to better constrain the timing of events (including 19th-century floods), and to dovetail this history with Hildreth and others’ (2003) bedrock study.

Washington↗

Development of a two-stage life cycle model for Oncorhynchus kisutch (coho salmon) in the upper Cowlitz River Basin, Washington

Recovery of salmon populations in the upper Cowlitz River Basin depends on trap-and-haul efforts owing to impassable dams. Therefore, successful recovery depends on the collection of out-migrating juvenile salmon at Cowlitz Falls Dam (CFD) for transport below downstream dams, as well as the collection of adults for transport upstream from the dams. Tacoma Power began downstream fish collection efforts at CFD in the mid-1990s and has been working consistently since then to improve collection efficiency to support self-sustaining salmon and steelhead ( Onchorhynchus spp.) populations in the upper Cowlitz River Basin. Although much work has focused on estimating fish collection efficiency (FCE), there has been relatively little focus on modeling population dynamics to understand how fish collection efficiency and other factors drive production of both juvenile and adult salmon over their life cycle. As a first step towards understanding the factors affecting population dynamics of Oncorhynchus kisutch (coho salmon) in the upper Cowlitz River Basin, we developed a statistical life cycle model using adult escapement and age structure data, juvenile collection data, and juvenile fish collection efficiency estimates. The goal of the statistical life cycle model is to estimate annual production and survival during two critical life-stage transitions: the freshwater production from escapement of adults upstream from CFD to collection of juveniles at CFD, and the juvenile-to-adult survival from the time of collection at the dam to the return of adults. To structure the life cycle model, we used the Ricker stock-recruitment model to estimate juvenile production from the number of parent spawners. This approach allowed us to account for density dependence at high spawner abundances while estimating annual productivity, defined as the number of juveniles produced per spawner at low spawner abundance. We then expressed productivity as a function two key variables affecting the number of juveniles collected and transported at CFD: (1) annual FCE, and (2) the annual number of days that spill occurred at CFD from September 1 to April 30. Our key findings were as follows: FCE was the primary factor affecting productivity of coho salmon upstream from CFD because FCE affects the number of juveniles that survive to continue downstream migration; Juvenile-to-adult return (JAR) rates were relatively high considering that harvest was included in the estimate, averaging about 3.6 percent and ranging as high as 9.1 percent, suggesting that adult coho salmon may be able to return to CFD at sustainable population sizes; and Much variation in the estimates of juvenile fish production upriver of CFD was unexplained even after adult escapement and FCE were accounted for, suggesting that the model may be improved by exploring different covariates and model structures for juvenile production as well as JAR rates. Additionally, by including FCE in the model, we estimated that the median pre-collection productivity, defined as the number of juveniles produced per spawner when FCE=1, was 108.4 juveniles per spawner. Because this two-stage life cycle model partitions factors that affect fish production in river compared to the ocean environment and fish life stages, the model estimates should help inform fishery managers about the overall role that fish collection at CFD may have on the recovery and sustainability of coho salmon populations.

Washington↗

Fish and habitat assessment in Rock Creek, Klickitat County, southeastern Washington, 2018

Executive Summary Native steelhead (anadromous form of rainbow trout [ Oncorhynchus mykiss ]) and bridgelip sucker ( Catostomus columbianus ) were historically used by the Kah-miltpah (Rock Creek) Band for sustenance, trade, and traditional practices in Rock Creek, a tributary to the Columbia River in southeastern Washington State. Rock Creek flows south to the Columbia River at river kilometer (rkm) 368 and is an intermittent stream of great significance to the Yakama Nation and to the Kah-miltpah Band in particular. Concern over declines in the abundance of these fish in Rock Creek prompted a research and monitoring program to better understand habitat conditions, population status, and limiting factors. In addition to steelhead and bridgelip sucker, coho salmon ( Oncorhynchus kisutch ) and resident rainbow trout are also present and monitored. Rainbow trout and steelhead will be collectively referred to as O. mykiss . Streamflow is a limiting habitat factor in this system, but steelhead and coho salmon still successfully return to spawn, rear, outmigrate, and survive over summer in many of the isolated pools that provide important refuge for juvenile rearing. We completed a habitat survey during autumn 2018 to assess the perennial pools during low-flow conditions. In Rock Creek, the overall percentage of habitat recorded as dry was 41, non-pool wet was 42, and pool was 17. The number of pools ( n =93) recorded was less than during previous years’ survey efforts (2015–17). The percentage of non-pool wet habitat was generally higher in 2018 than in previous years. This is a likely result of habitat reaches, which in the past, were considered pools but have become shallower and smaller and are now categorized as non-pool wet habitat. However, the fewer habitat reaches categorized as pools in 2018 now have an average length, area, and depth that are generally greater than in past years. In Walaluuks Creek, the percentage of habitat recorded as dry was 53, non-pool wet was 40, and pool was 7. The percentage of pool habitat was the lowest of all years surveyed since 2015. Fish sampling occurred during autumn after habitat surveys were completed from October 1 to November 9. Fish species distribution, relative abundance, length-frequency distribution, and pool fish density were determined using backpack electrofishing in stratified, systematically selected pools. During fish sampling, 855 O. mykiss were handled and 662 were tagged with a passive integrated transponder (PIT) tag, and 718 coho salmon were handled and 567 were PIT tagged. A total of 536 bridgelip suckers and largescale suckers ( Catostomus macrocheilus [ n =6]) were handled and 294 were PIT tagged. In Rock Creek, pool abundance estimates were calculated for six pools for both O. mykiss age classes (age 0 and age 1 or older [age 1+]) and one additional pool for age 1+. For pools where age-0 O. mykiss were present, the average pool population abundance was 0.144 ( n =6; range: 0.052–0.208) fish per square meter. For age-1+ O. mykiss , the average pool population abundance was 0.045 ( n =7; range: 0.002–0.179) fish per square meter. For age-0 O. mykiss in Walaluuks Creek, the average pool abundance was 0.207 fish per square meter ( n =7; range: 0.038–0.416), and for age-1+ fish, the average pool abundance was 0.382 fish per square meter ( n =6; range: 0.009–0.761). In Rock Creek, coho salmon were more abundant than O. mykiss in pools except for three pools upstream from rkm 20. The average pool abundance for coho salmon was 0.256 fish per square meter ( n =8; range: 0.019–0.756) in Rock Creek pools. In Walaluuks Creek, coho salmon were captured in four pools and were not captured in the upstream pools sampled. The average pool abundance for coho salmon in the four lower pools was 0.488 fish per square meter ( n =4; range: 0.417–0.548). Bridgelip suckers were captured in all pools in Rock Creek except the pool sampled at rkm 21.8. The average pool abundance for bridgelip suckers was 0.552 fish per square meter ( n =7; range: 0.015–1.554) in Rock Creek. Bridgelip suckers were captured in three downstream pools in Walaluuks Creek, and the average abundance was 0.044 fish per square meter ( n =3; range: 0.024–0.085). Overwinter and reach survival probabilities were estimated for O. mykiss and coho salmon using a Cormack-Jolly-Seber modeling approach. The best fit survival model for the O. mykiss and coho salmon was a reach only model. The upstream reach includes overwinter survival probability because fish are tagged and released in autumn and primarily migrate the following spring. During 2018, coho salmon (0.568, standard error [SE]=0.027) had a significantly higher probability of overwinter survival than O. mykiss (0.276, SE=0.019). The reach survival probability was higher for O. mykiss than coho salmon in the downstream migratory reaches. Survival was not modeled for bridgelip suckers. For bridgelip suckers, 147 were detected of 294, that were PIT tagged and released in the Rock Creek subbasin (50.0 percent). Information provided in this report increases our understanding of the status and trends of these populations. It further documents how intermittent streams can support salmonid populations. It also provides insight into potential management and restoration actions that could be beneficial and timing and allocation of resources. Ongoing monitoring work of this population will inform progress towards Rock Creek species recovery goals and contribution to recovery goals for the steelhead Middle Columbia River Distinct Population Segment.

Washington↗

Streambed scour of salmon (Oncorhynchus spp.) and steelhead (Oncorhynchus mykiss) redds in the South Fork Tolt River, King County, Washington

Prior to emergence as fry, salmonid embryos incubating within gravel nests called “redds” are vulnerable to substrate mobilization and lowering of the streambed, a process termed “streambed scour,” during floods. Water managers regulating discharge in salmonid-bearing rivers need information about the magnitude of discharge during which the scour of substrate surrounding salmonid redds occurs. The time when scour occurs, however, is difficult to measure and usually poorly constrained. The South Fork Tolt River in western Washington supplies the City of Seattle with hydroelectric power and about 40 percent of its municipal water needs, while providing spawning habitat for two salmonid species listed under the Endangered Species Act: Chinook salmon ( Oncorhynchus tshawytscha ) and steelhead trout ( O. mykiss ). The U.S. Geological Survey, in cooperation with Seattle City Light and Seattle Public Utilities, began a study in 2015 using accelerometer scour monitors (ASM) to characterize the timing of and hydrologic conditions associated with streambed scour at the depth of incubating salmonid embryos in the South Fork Tolt River. Prior to this study, operational thresholds for peak discharge on the South Fork Tolt River were 350 cubic feet per second (cfs) in the upper part of the river and 550 cfs in the lower part of the river as measured at USGS streamgages 12148000 and 12148300, respectively. These thresholds were developed from the peak discharge associated with observations of the flattening of redd structure and not from direct measurement of scour at the depth of egg pockets within redds. Accelerometer scour monitors were deployed at the level of salmonid egg pockets in spawning habitat of the South Fork Tolt River to record the temporal pattern of streambed scour at the depth of incubating salmon eggs during fall and winter flood seasons of water years (WY) 2016 and 2017. Thirteen of 48 ASMs deployed during the WY 2016 flood season recorded scour attributed to high streamflow when discharge measured at USGS streamgage 12148300 (the lower river streamgage used as an index gage) was between 969 and 1,360 cfs. Local discharge at individual scour sites varied depending on the timing of tributary inputs and downstream transport of water. During the subsequent flood season in WY 2017, peak discharge at the index gage reached 809 cfs. None of the 38 ASMs deployed recorded scour attributed to streamflow alone, although 10 ASMs recorded localized bed movement attributed to spawning activity of fish. Most scour at the depth of redds measured during WY 2016 occurred at or before peak flood discharge consistent with previous redd scour studies. The lack of scour measured in WY 2017 when peak discharge (809 cfs) was less than the minimum discharge when scour occurred in WY 2016 (969 cfs) suggests minimal to no scour of egg pockets in salmonid redds when discharge is less than 809 cfs.

Washington↗