Andrew S. Gendaszek

The Enloe Dam was built in 1920 for hydropower generation and impounds a steep-banked, narrow reach of the Similkameen River in north-central Washington. During the subsequent century, sediment from the Similkameen River watershed, which includes historical mining operations, accumulated within Enloe Dam’s impoundment. Enloe Dam ceased hydropower production in 1958 and is currently under consideration for removal, but there are concerns that the remobilization of sediments may have harmful impacts on downstream water quality. To complement previously published analyses of heavy metal concentrations within sediments and assess the total volume that may be transported following dam removal, this report presents estimates of the volume of sediment trapped behind Enloe Dam which was measured in 2020. The volume of sediment was estimated by comparing a bathymetric survey, collected using an acoustic Doppler continuous profiler, and a survey of the bedrock-sediment interface, collected using a continuous resistivity profiler.The study area spanned a 2.6-kilometer reach of the Similkameen River upstream from Enloe Dam. The volume of impounded sediment was calculated in 2020 by subtracting the elevation of the bedrock-sediment interface measured using the continuous resistivity profiler from a bathymetric surface measured by the acoustic Doppler current profiler. In 2020, the estimated volume of sediment impounded by Enloe Dam was 2.17±0.04 million cubic meters (Mm3) compared to 1.37 Mm3 measured in 1971. This equates to a deposition rate of approximately 16,300 cubic meters per year from 1971 to 2020. Continuous resistivity profiles revealed that bedrock within Enloe Dam’s impoundment was deepest (about 20 m) in the southern, downstream end of the profiles and shallowest (less than 5 m) in the northern upstream end of the profiles.

Water temperature is a primary control on the occurrence and distribution of fish and other ectothermic aquatic species. In the Pacific Northwest, cold-water species such as Pacific salmon (Oncorhynchus spp.) and bull trout (Salvelinus confluentus) have specific temperature requirements during different life stages that must be met to ensure the viability of their populations. Rivers draining Mount Rainier in western Washington, including the White River along its northern flank, support a number of cold-water fish populations, but the spatial distribution of water temperatures, particularly during late-summer baseflow during August and September, and the climatic, hydrologic, and physical processes regulating it are not well constrained. Spatial stream network (SSN) models, which are generalized linear models that incorporate streamwise spatial autocovariance structures, were fit to mean and 7-day average daily maximum water temperature for August and September for the White River Basin. The SSN models were calibrated using water temperature measurements collected in 2010 through 2020. The extent of the models included the White River and its tributaries upstream from its confluence with Silver Creek in Mount Rainier National Park, Washington. SSN models incorporated covariates hypothesized to represent the climatic, hydrologic, and physical processes that influence water temperature. SSN models were fit to the measured data and compared to generalized linear models that lacked spatial autocovariance structures. Statistically significant covariates within the best-fit models included the proportion of ice cover and forest cover within the basin, mean August air temperature, the proportion of consolidated geologic units, and snow-water equivalent. Statistical models that included spatial autocovariance structures had better predictive performance than those that did not. Additionally, models of mean August and September water temperature had better predictive performance than those of 7-day average daily maximum temperature in August and September. Predictions of the spatial distribution of water temperature were similar between August and September with a general warming in the downstream part of the mainstem White River compared to cooler water temperatures in the high-elevation headwater streams. The proportion of ice cover emerged as an inversely related significant covariate to both mean August and September water temperature because streams that receive glacial meltwater are colder than non-glaciated streams. Water temperatures of the upper White River increased downstream and are attributed to warming of water temperature from accumulated solar radiation and inflow of non-glaciated tributaries. Estimated water temperatures for the upper White River model are 3–4 degrees Celsius (°C) warmer for tributaries, but 1–2 °C cooler for the mainstem compared to the regional-scale model. Differences between the upper White River SSN model and the regional-scale NorWeST model are attributed to the fact that the upper White River SSN included water temperature observations specific to the upper White River, whereas water temperature observations from lower elevation streams and downstream from the Mount Rainer National Park boundary were used in the regional scale model.

More than 1 million people live within the active model area (AMA) in the southeastern part of the lowlands surrounding Puget Sound, or Puget Lowland, Washington, and groundwater is the source for approximately one-half of their public, domestic, and irrigation water demands. The 887-square-mile AMA, located in King and Pierce Counties, represents the area of analysis for the conceptual hydrogeologic framework and numerical groundwater-flow models within the study area and includes the Puyallup River and Chambers-Clover Creek watersheds. To assess the potential hydrologic and anthropogenic impacts to groundwater and the connected surface-water resources, conceptual and numerical groundwater-flow models of groundwater flow were developed by the U.S. Geological Survey Washington Water Science Center in close cooperation with 18 water-resource agencies and stakeholders.This report presents information used to characterize the groundwater-flow system and the development of a numerical model in the AMA. Included are descriptions of the geology and conceptual hydrogeologic framework, groundwater levels and flow directions, groundwater recharge and discharge, numerical groundwater-flow model construction and results, and model limitations. The study area encompasses the western part of Pierce County and the southwestern part of King County, Washington. The study area extends south to the Nisqually River, southwest to Tanwax Creek, northeast to the Green River, and north through the valley near Auburn and adjacent uplands. It is bounded on the east by foothills of the Cascade Range, and on the northwest by Puget Sound.

Accurate mapping of streams that maintain surface flow during annual baseflow periods in mountain headwater streams is important for informing water availability for human consumption and is a fundamental determinant of in‐channel conditions for stream‐dwelling organisms. Yet accurate mapping that captures local spatial variability and associated local controls on surface flow presence is limited. An empirical random‐forest model was developed to predict streamflow permanence (late summer surface‐flow presence) for Mount Rainier National Park and the surrounding mountainous area in western Washington, USA. This model was developed to improve upon the existing multi‐state, regional‐scale probability of stream permanence developed for the greater Pacific Northwest Region (PROSPER_PNW). The model was trained on 544 wet/dry observations collected during the late summer, baseflow period from 2018 to 2020 using the crowd‐source mobile application, FLOwPER. Final model accuracy was 0.74 with drainage area and covariates describing geology, topography, and land cover as top predictors of streamflow permanence compared to coarser resolution climatic covariates. The prevalence of static covariates over climatic covariates as top ranked important covariates highlights the importance of scale when evaluating controls on streamflow permanence. Cross validation of the model indicates that streamflow permanence probabilities from this model is an improvement over the regional‐scale PROSPER_PNW model demonstrating the utility of relatively simple, crowd‐sourced data to address water resource needs, and that determination of important predictors of streamflow permanence is influenced by the spatial and temporal resolution of analysis.

In streams supporting Pacific salmon (Oncorhynchus spp.) within the southern Puget Lowland, high water temperatures during late summer are a primary water-quality concern. The metabolic rates of fish and other ectothermic (in other words, cold-blooded) species are regulated by water temperature; salmon and other cold-water fish have specific thermal tolerances outside of which they are susceptible to infection, disease, increased predation, and decreased reproductive success. Mill and Gosnell Creeks, which collectively drain a 30-square mile area of the Puget Lowland in Mason County, Washington, support several species of anadromous salmonids. Whereas previous studies documented relatively cool water temperatures in Gosnell Creek, which drains the watershed upstream from Lake Isabella, water temperatures in Mill Creek, which heads at the outlet of Lake Isabella, regularly exceed thermal tolerances for cold-water fish. The occurrence and distribution of cold-water anomalies in less-than-ambient water temperatures in Mill Creek, however, have not been assessed. In this report, we present spatially and temporally continuous measurements of near-streambed water temperature measured using fiber-optic distributed temperature sensing for three reaches of Mill Creek during August–September 2020 when the water temperatures of streams in western Washington were near their annual maximum. Water temperature was collected every hour and averaged spatially over 1.015-meter sections of the fiber-optic cable deployed at the streambed of Mill Creek. The lengths of the fiber-optic cables deployed in Reaches A, B, and C were 883, 270, and 1,014 meters, respectively. Daily maximum water temperature and daily temperature variability, as measured by standard deviation of water temperature during the deployment, progressively decreased downstream as distance from Lake Isabella increased. However, no abrupt decreases in daily maximum or standard deviation of water temperature were detected in longitudinal temperature profiles of any of the three reaches. Collectively, these results suggest that warm water discharged from Lake Isabella was progressively buffered downstream as it equilibrated with downstream heat fluxes mediated by physical processes including riparian shading and diffuse groundwater input. Although parts of the surveyed reaches associated with deep pools were cooler than other locations, no large (less than 2 °C) water-temperature anomalies characteristic of discrete sources of cold groundwater or surface-water inputs were measured in any of the three surveyed reaches.

Warm water temperatures in the lower Yakima River in central Washington are key limitations to the restoration of Pacific salmon (Onchorhynchus spp.) populations within the Yakima River Basin. Identification of the location and magnitude of cold-water anomalies, which are cooler than ambient river temperatures during summer months, and the processes that create and maintain them is needed to inform salmon restoration efforts within the Yakima River Basin. Longitudinal thermal profiles of nine reaches in the lower Yakima River were surveyed at ambient river velocity during summer 2018 when surface-water temperatures were near their annual maximum and the difference between surface-water and groundwater temperatures was greatest. The profiles were compared to previously published profiles of the same reaches measured in 2001, 2002, 2008, and 2009, and analyzed in the context of hydrologic, geomorphic, and hydrogeologic conditions that may create and maintain cold-water inputs to the river. Cold-water anomalies that departed from expected diurnal increases in water temperature were measured in all nine study reaches and were attributed to diffuse groundwater discharge through the streambed, discrete groundwater discharge at seeps and springs, and cold-water tributaries entering the river. Some cold-water anomalies were measured during repeated surveys in different years, whereas other cold-water anomalies did not persist across surveys. Additionally, some discrete cold-water anomalies were confined to one side of the channel, but others associated with diffuse groundwater discharge were present across the channel for several river miles. Hydrogeologic conditions including the extent and thickness of aquifers connected to the Yakima River, geomorphic conditions including channel gradient, channel geometry, and floodplain extent, and the location of tributaries, irrigation returns, and other surface-water inputs created the large-scale conditions that facilitate the formation and maintenance of cold-water anomalies. Finer-scale geomorphic features such as side channels, gravel-bar alcoves, deep pools, and other locations, where colder water collected and remained relatively unmixed with upstream surface water, were also important factors in the occurrence and distribution of cold-water anomalies. These hydrogeologic and geomorphic conditions, coupled with the alteration of the Yakima River’s hydrologic regime to support irrigation within the Yakima Valley, contributed to the surveyed distribution of cold-water anomalies within the river.

First posted December 16, 2021 Director, Washington Water Science Center U.S. Geological Survey 934 Broadway, Suite 300 Tacoma, Washington 98402 Sockeye salmon (Oncorhynchus nerka) spawn at beaches along Lake Ozette’s shoreline and within its tributary streams including Umbrella Creek and Big River in western Washington. The tributary-spawning aggregate of the Lake Ozette sockeye salmon population has been increasing from very low abundance through hatchery supplementation, but the beach-spawning aggregate has decreased from the early 20th century resulting in an Endangered Species Act listing of the Lake Ozette sockeye salmon population as “Threatened” in 1999. Among several factors inhibiting the recovery of beach spawning sockeye salmon, the quality of spawning habitat in beaches and low dissolved-oxygen concentrations in beach gravels during incubation were identified as important limitations on the recovery of this population. Proliferation of primarily native near-shore vegetation during the 20th century, as a result of alterations in the lake hydroperiod, accompanied by fine-grained sediment deposition was hypothesized as a potential cause of low rates of water circulation and dissolved-oxygen concentrations in beach spawning gravels. The potential for shoreline vegetation removal to restore spawning habitat function, including dissolved-oxygen concentrations, was evaluated in this report by measuring continuous dissolved-oxygen concentrations with data-logging dissolved-oxygen sensors, by measuring particle-size distribution of beach sediment, and by estimating groundwater/surface-water exchange using vertical sediment temperature profiles at three shoreline sampling areas. These sampling areas included an area of current spawning devoid of shoreline vegetation, an adjacent vegetated area, and an adjacent treatment area where a 3-meter-wide swath of existing above-ground vegetation was removed in 2018 prior to the study. Substrate particle-size distributions, dissolved-oxygen concentrations, sediment temperatures, and groundwater/surface-water exchange were compared among the three shoreline sampling areas. Median grain size (D₅₀) of sediment varied at sampling stations from medium sand fine to coarse gravel. The coarsest sediment generally occurred in the current spawning area that was devoid of vegetation; whereas the vegetated shoreline and the shoreline where above-ground vegetation was removed were characterized by finer sediment. Removal of above-ground vegetation resulted in increased D₅₀ at the most shoreward station at the treatment sampling area from 8.2 millimeters in 2018 to 21.6 millimeters in 2019 but other changes in substrate particle-size distribution in the treatment area were negligible. Increased grain size from 2018 to 2019 at this site suggests that while higher wave energy was allowed to mobilize sediment in the backshore area of the treatment area and winnow fine sediment during the winter, residual root structure in the treatment area may have limited the ability of wave energy to mobilize sediment after removal. During the November 2018 to March 2019 incubation period for sockeye salmon, dissolved-oxygen concentrations at the depth of sockeye salmon egg pockets (15–25 centimeters) within all three shoreline sampling areas were less than 1 milligram per liter throughout the deployment time (October 2018—May 2019) and below the threshold to sustain sockeye salmon embryo development (3 milligrams per liter). In addition, the similarity of dissolved-oxygen concentrations among all three shoreline sampling areas indicates that above-ground vegetation removal did not increase subsurface dissolved-oxygen concentrations. Groundwater/surface-water exchange measured from sediment temperature profiles were variable both within and across shoreline sampling areas. At the most lakeward stations, groundwater discharge to the lake ranged from 0.25 to 0.007 meter per day and was highest at the control station and lowest at the vegetated station. However, in general, the differences in groundwater/surface-water exchange across the three shoreline sampling areas were negligible. Collectively, these results suggest the process of removing above-ground vegetation had little effect on subsurface dissolved-oxygen concentrations and groundwater/surface-water exchange during the study period, but limitations of the study design, including retention of below-surface root cohesion after above-ground vegetation removal, too narrow of a band of vegetation removal, and a limited duration of the monitoring period, may have prevented wave energy from winnowing fine-grained sediment along the shoreline and altering subsurface dissolved-oxygen concentrations during the study period. Response of the substrate particle-size distribution, groundwater/surface-water exchange, and subsurface dissolved-oxygen concentrations to shoreline vegetation removal that includes root-zone removal over a larger extent and longer periods than the 7-month study period from October 2018 to May 2019, however, remain unknown and warrant further investigation.

The autumn and winter flood season of western Washington coincides with the incubation period of many Pacific salmon (Onchorhynchus spp.) populations. During this period, salmon embryos incubating within gravel nests called “redds” are vulnerable to mobilization of surrounding sediment during floods. As overlying sediment is transported downstream, the vertical position of the streambed can be lowered, a process termed streambed scour; thus developing salmon embryos may be destroyed resulting in decreasing egg-to-fry survival rates. The Sauk River, which drains a 1,900 km2 (733.5 mi2) area of the central Cascade Range of Washington State, provides spawning and rearing habitat for several species of Pacific salmon including Chinook salmon (O. tshawytscha), which were listed as threatened under the Endangered Species Act (ESA) in 1999. In order to assess the hydrologic conditions when streambed scour and concomitant geomorphic changes occur, accelerometer scour monitors (ASMs), which record the time when streambed scour lowers the streambed to the level of salmon egg pockets, were deployed in two geomorphically different reaches of the Sauk River to monitor scour during water year 2018. Nineteen ASMs were deployed in an upstream reach, which was largely confined by valley walls with vegetated, stable banks and low channel-migration rates near the confluence of the Sauk and White Chuck Rivers. Twelve additional ASMs were deployed in a downstream reach within an unconfined valley with unvegetated, unstable banks and high channel-migration rates between the town of Darrington and the confluence of the Sauk and Suiattle Rivers. During the ASM deployment, discharge measured at the U.S. Geological Survey (USGS) streamgage Sauk River above White Chuck River, near Darrington, Washington (12186000), peaked at 479 m3/s (16,900 ft3/s) with an estimated 0.18 probability of annual exceedance (5.7-year recurrence interval). During the flood season, large-scale geomorphic changes, including channel migration and bar deposition, were measured at the downstream reach, but only minimal geomorphic changes were measured at the upstream reach. ASMs deployed at the downstream reach were not recovered after the flood season and total scour depth was presumed to have exceeded ASM anchor depth. At the upstream reach, 7 of the 19 deployed ASMs were recovered after the flood season and all recovered ASMs recorded scour at discharges that equaled or exceeded 204 m3/s (7,210 ft3/s). The remaining 12 ASMs deployed at the upstream reach were not recovered and total scour depth was presumed to have exceeded ASM anchor depth. Collectively, this analysis enhances the ability of fisheries managers to forecast egg-to-fry survival rates of salmonids by determining the hydrologic conditions at which scour at the level of salmon redds initiates.

The seasonal and inter-annual variability of flow presence and water temperature within headwater streams of the Great Basin of the western United States limit the occurrence and distribution of coldwater fish and other aquatic species. To evaluate changes in flow presence and water temperature during seasonal dry periods, we developed spatial stream network (SSN) models from remotely sensed land-cover and climatic data that account for autocovariance within stream networks to predict the May to August flow presence and water temperature between 2015 and 2017 in two arid watersheds within the Great Basin: Willow and Whitehorse Creeks in southeastern Oregon and Willow and Rock Creeks in northern Nevada. The inclusion of spatial autocovariance structures improved the predictive performance of the May water temperature model when the stream networks were most connected, but only marginally improved the August water temperature model when the stream networks were most fragmented. As stream network fragmentation increased from the spring to the summer, the SSN models revealed a shift in the scale of processes affecting flow presence and water temperature from watershed-scale processes like snowmelt during high-runoff seasons to local processes like groundwater discharge during sustained seasonal dry periods.