USGS Science⌕ Search

Geology topics

Annett B. Sullivan

Publications and source records attributed to Annett B. Sullivan.

At least 19 recordsLinked to original sources
USGS2025

Review of the Lake Washington Ship Canal and Ballard Locks model, Seattle, Washington, 2014–20

Executive Summary The Hiram M. Chittenden (Ballard) Locks and Lake Washington Ship Canal connect freshwater Lake Washington and saline Shilshole Bay of Puget Sound in Seattle, Washington. The locks and canal allow for ships to traverse this reach. Anadromous salmonids also migrate through, transitioning between saline and freshwater environments, and making use of a fish ladder at the locks when traveling upstream. WEST Consultants, Inc., constructed a two-dimensional hydrodynamic and water-quality model (CE-QUAL-W2) simulating flow, water temperature, and salinity for the Ballard Locks and the Lake Washington Ship Canal. An initial model was built for calendar years 2014–15, and the model was updated using a more recent and modern dataset for calendar years 2016–20. The U.S. Army Corps of Engineers requested that the U.S. Geological Survey review this model and its documentation to evaluate the technical aspects of its development and calibration. Findings from this review include the following: Overall, the Lake Washington Ship Canal CE-QUAL-W2 model was well-documented and constructed largely following typical model-development methods. The Lake Washington Ship Canal model was built with CE-QUAL-W2 model version 4.5, compiled and released by Portland State University in April 2021. CE-QUAL-W2 updates and improvements are regularly released with bug fixes and new features, so any model updates would benefit from the use of the most-recent software release. The model grid that represents the Lake Washington Ship Canal bathymetry was 9.2 kilometers (5.7 miles) long, matching the expected length of the waterway. The deepest model segments were near sampling site LLLW (Large Locks site) near the locks. Lake Union is reported to constitute most of the volume of the Lake Washington Ship Canal and is depicted as such in the model grid. The model includes several water outflows at Ballard Locks, including the large and small locks, a saltwater drain, a spillway, smolt flumes, and a fish ladder. Flows from the spillway, smolt flumes, and fish ladder were combined into one structure outflow in the model and assigned one withdrawal elevation from the Lake Washington Ship Canal. The smolt flume and spillway withdraw from the same elevation, but the fish ladder flow withdraws from a higher elevation in Lake Washington Ship Canal, and that flow could be separated into its own withdrawal. The model input files were created using the Coordinated Universal Time standard instead of the more typical choice of using local standard time. This is not incorrect, but sub-daily results would need to be converted to local time for science-communication purposes. The meteorological dataset had some unexpected anomalies, such as a baseline shift in the wind-speed dataset. Other nearby meteorological datasets could be used instead or used to correct the current meteorological inputs. The upstream boundary was configured with water-temperature data from a continuous monitor buoy in Lake Washington. The boundary salinity was set at 0 parts per thousand for the duration of the model simulation. A more realistic estimate of salinity at the upstream boundary could be constructed using data from the same buoy. Saline inflow at the downstream boundary of the Lake Washington Ship Canal model through lock exchanges at the large lock was included as a tributary in the model. Salinity and temperature inputs in this tributary at the large locks were set as constant values for the entire simulation. Saline inflow through the small lock was not included in the model because few data were available, and the input was likely to be small because of the smaller surface area and volume of the small lock relative to the large lock. The model did not include any flow, water temperature, or salinity inputs to the Lake Washington Ship Canal other than at the locks and at the upstream boundary. Any point sources, small tributaries, or stormwater inputs were omitted from the model. It is unclear whether this is a substantial omission relative to model results. Most model parameters were set as defaults or to reasonable values. However, the value of the WINDH parameter, the height of the wind speed measurement, was different than the height of the meteorological site. Compared to measured data, the model simulated water-surface elevations and water temperatures with reasonable accuracy. Differences in the modeled and measured salinities revealed some opportunities to improve the simulation of salinity, both baseline salinity and the salinity maxima in summer and autumn.

Washington↗
USGS2023

Modeling the water-quality effects to the Klamath River from recirculation in drains and canals, Oregon and California, 2006–15

The potential recirculation of Klamath Strait Drain (hereafter called by its local name, “Klamath Straits Drain”) water into Ady Canal to reduce the drain discharge of high nutrient loads into the Klamath River was assessed by the U.S. Geological Survey for the Bureau of Reclamation. To study the feasibility of recirculation, this investigation evaluated three recirculation scenarios over a 10-year period from 2006 to 2015, as a series of 1-year model simulations. A combination of two existing hydrodynamic, water-temperature, and water-quality models (CE-QUAL-W2) were used, including (1) the Link-Keno reach of the Klamath River, using Klamath Straits Drain as a tributary and for calendar years 2006–11, and (2) the same Link-Keno model used for calendar years 2012–15 in combination with an independent Klamath Straits Drain model from 2012 to 2015. Model simulations using the water-quality models were configured for the base case conditions and three different sets of recirculation scenarios: the maximum year-round recirculation without limits (scenario 1), limited year-round recirculation fixed by the current pipe flow configuration (scenario 2), and limited seasonal recirculation (May–September) also fixed by the current pipe flow configuration (scenario 3). In the base case, estimates of annual average daily total nitrogen loads and daily total phosphorus loads exported to the Klamath River from the Klamath Straits Drain were as much as 3,060 and 457 pounds per day (lbs/day), respectively. Currently (2023), the Total Maximum Daily Loads allocations for the Klamath Straits Drain are 21 and 268 lbs/day for total phosphorus and total nitrogen, respectively, so these maximum estimates exceed the current Total Maximum Daily Loads by greater than an order of magnitude. With scenario 1, load reductions occurred year-round for all constituents evaluated (total nitrogen, total phosphorus, 5-day biochemical oxygen demand [BOD5], 5-day carbonaceous biochemical oxygen demand) for the Klamath Straits Drain discharging to the Klamath River. Scenario 2 also had large reductions in total nitrogen, total phosphorus, and BOD5 loads. Substantial reductions did occur for scenario 3 but were constrained to only the active recirculation period from May through September. Despite the restricted period, the average reductions in the annual average daily load for total phosphorus and total nitrogen were 32.1 percent and 26.5 percent, respectively. The Ady Canal diverts high nutrient loads from the Klamath River, so the loading tradeoffs to the Klamath River between no recirculation and the recirculation scenarios were calculated. On an annual basis, the overall net balance between the Klamath Straits Drain and Ady Canal resulted in more total nitrogen and total phosphorus load reductions to the Klamath River for the three recirculation scenarios than the base case, for most years. In contrast, the net balance for BOD5 loads was higher to the Klamath River for the three recirculation scenarios than the base case, for most years. With the recirculation scenarios, the optimal recirculation periods to benefit Ady Canal, Klamath River, and Klamath Straits Drain did not always coincide. Recirculation would be most effective at reducing loads toward the Klamath Straits Drain Total Maximum Daily Load allocations in the spring (March–May) of each year. However, recirculation during these months would also increase salinity in the Ady Canal. In summer, recirculation would reduce Klamath Straits Drain loads toward the Total Maximum Daily Load allocations, though recirculation could decrease Klamath River water quality mostly because of decreased withdrawals of Klamath River water by the Ady Canal. Scenario 3 avoided recirculation into Ady Canal in the early spring months when salinity concerns would be the highest, while still decreasing nutrient loads exported from the Klamath Straits Drain to the Klamath River in the summer months.

California, Oregon↗
USGS2023

Modeling flow and water quality in reservoir and river reaches of the Mahoning River Basin, Ohio

Executive Summary The U.S. Army Corps of Engineers (USACE) is considering changes to the management of water surface elevation in four lakes in the Mahoning River Basin. These changes would affect the timing and amounts of water released to the Mahoning River and could affect the water quality of those releases. To provide information on possible water-quality effects from these operational changes, flow and water-quality models were constructed for Berlin Lake, Lake Milton, Michael J Kirwan Reservoir, Mosquito Creek Lake, Mosquito Creek, and the Mahoning River from the dams downstream to Lowellville, Ohio. The models were calibrated for two calendar years each, with model years selected depending on the availability of water-quality data. Models were developed with CE-QUAL-W2 version 4.2 (Wells, S.A., 2020, CE-QUAL-W2—A two-dimensional, laterally averaged, hydrodynamic and water quality model [version 4.2]: Portland State University, variously paged), a two-dimensional, laterally averaged hydrodynamic and water-quality model. Modeled constituents included flow, velocity, ice cover, water temperature, total dissolved solids (TDS), sulfate, chloride, inorganic suspended sediment, nitrate, ammonia, total Kjeldahl nitrogen, orthophosphate, total phosphorus, dissolved and particulate organic matter, algae, and dissolved oxygen. Iron was included for the lake models, but not the river. A whole-basin model, with the four lake models and river model, was used to run model scenarios to examine the effects of altered lake water surface elevations on flow and water quality in the lakes, the lake outflows, and the Mahoning River. The initial whole-basin model, with calendar year 2013 hydrology and measured or typical water quality, was designated as scenario 0. Mahoning River flows for calendar year 2013 were close to a 20-year median flow. Four additional scenarios were constructed based on reservoir operations model (RES-SIM) model water surface elevations for the four lakes as provided by USACE. Scenario 1 was the RES-SIM base case, scenario 2 kept Berlin Lake water surface elevations higher in summer, scenario 3 allowed 25 percent of summer flood storage to extend the guide curve, and scenario 4 allowed more flexibility in lake management by removing any downstream Mahoning River minimum flow requirements. The Mahoning River model was not changed in any scenarios but received altered flows from the lakes. Significant findings from this study include the following: In two of the four lakes (Berlin and Mosquito Creek Lakes), development of lake model grids using recent bathymetric surveys suggests that sedimentation in these lakes has occurred since they were constructed, altering volume-elevation curves. Tests of model parameter sensitivity showed that modeled water temperature, TDS, and dissolved oxygen were relatively insensitive to model parameter values. Modeled chlorophyll a , a measure of algal concentration, was most sensitive to parameter values; nitrate and total phosphorus concentrations were affected by a few of the parameters tested. As a group, the lake model results were more sensitive to model parameter values compared to the Mahoning River model. Data gaps were identified for inflows, both for water quantity and water quality, that could be filled through future sampling programs. Ample data were available from within the waterbodies for model calibration. The model simulated the general spatial and temporal patterns of water temperature, TDS, chloride, sulfate, nutrients, suspended sediment, organic matter, chlorophyll a , and dissolved oxygen in the lakes and Mahoning River. From late spring to autumn in the years modeled (2006, 2013, 2017–19 depending on the lake), all lakes developed thermal stratification and periods of anoxia in bottom waters. Stratification was most stable in Michael J Kirwan Reservoir and least stable in Mosquito Creek Lake. The stratification and anoxia in Berlin Lake, Lake Milton, and Mosquito Creek Lake could be interrupted by high-flow inputs moving through those lakes. The model predicted the release of ammonia and iron during anoxic periods in the lake hypolimnions. Concentrations of TDS, nitrate, orthophosphate, and total phosphorus increased in the Mahoning River down to Lowellville, the end of the river model, in the years modeled. These concentrations were greater than those in upstream lake releases. Chloride and sulfate concentrations were underpredicted in the Mahoning River, suggesting the presence of unreported chloride and sulfate inputs to the river, at least in the years modeled. Model scenario 4 kept water surface elevations the highest in all lakes in the April to mid-December period, compared to scenarios 1–3. Model scenario 2 kept water surface elevations in Berlin Lake higher in summer and late autumn, compared to scenarios 1 and 3, but to satisfy downstream minimum flow requirements, water surface elevations in the other lakes had periods of lower water surface elevation. As a group, scenarios 1–3 had largely similar effects on flow and water surface elevation in the Mahoning River because the lake releases in those scenarios still met downstream Mahoning River flow targets. Modeling the removal of downstream flow targets, scenario 4 had periods of lower flow in the Mahoning River from April to mid-September as water was held in the lakes, and periods of higher Mahoning River flow from mid-September through November as the lakes were drawn down to prepare for winter flood-risk management. In the four scenarios, all the lakes and lake outflows had generally similar seasonal cycles of water quality, though some differences were predicted. For instance, higher concentrations of iron and ammonia in the Lake Milton hypolimnion were modeled during a period of both low inflows from Berlin Lake and low outflows at Lake Milton dam. It is possible that those changes could be minimized by maintaining more flow or pulses of higher flow through the lake. Compared to the scenario 1 base case, changes to Mahoning River water quality were relatively minor for scenarios 2 and 3, which maintained downstream flows but shifted the flow source among the upstream lakes. The largest changes in Mahoning River water quality were predicted between Leavittsburg and Lowellville for scenario 4. The periods of lower lake outflows between April and mid-September led to correspondingly higher concentrations of TDS, orthophosphate, total phosphorus, and nitrate in the river, compared to the base case scenario 1. Conversely, the overall greater lake outflows from mid-September through November in scenario 4 led to periods of lower concentrations of TDS and nutrients in that portion of the river, at that time of year.

Ohio↗
USGS2021

Modeling water temperature response to dam operations and water management in Green Peter and Foster Lakes and the South Santiam River, Oregon

Significant Findings Green Peter and Foster Dams have altered natural seasonal temperature patterns in the South and Middle Santiam Rivers of the Willamette River Basin in northwestern Oregon. Cold-water releases from Green Peter Dam, upstream of Foster Lake, contribute to the cool-water conditions at Foster Dam. In summer, unseasonably cold water typically is discharged from Foster Dam into the Foster Dam fish ladder, which may be one factor contributing to the low numbers of upstream migrating Chinook salmon ( Oncorhynchus tshawytscha ) that enter the fish ladder. The U.S. Army Corps of Engineers is leading efforts to improve conditions for Chinook salmon upstream and downstream of these dams by considering structural alterations to Foster Dam and by exploring changes to the way the dams are operated. The U.S. Geological Survey assisted the U.S. Army Corps of Engineers by using previously calibrated numerical models of flow and water quality for Green Peter and Foster Lakes and for the South Santiam River downstream of Foster Dam. These two-dimensional hydrodynamic and water-quality (CE-QUAL-W2) models were used to test scenarios of altered dam operations and alternate water-management strategies. Results of these scenarios provide information and insights into how the mixing and thermal characteristics of the lakes are affected by dam operations, how the mixing and timing of upstream source waters reaching Foster Dam are affected by dam operations, how river and fish-ladder temperature targets might be achieved, and how quickly (or slowly) such changes in the lakes and downstream river reaches occur, relative to typical unmodified operations at Green Peter and Foster Dams.

Oregon↗
USGS2020

Evaluation of restoration alternatives using hydraulic models of lake outflow at Wapato Lake National Wildlife Refuge, northwestern Oregon

Wapato Lake National Wildlife Refuge near the city of Gaston in northwestern Oregon was established in 2013, and planning is underway to restore a more natural lake and wetland system after more than 100 years of agricultural activity on the lakebed. Several water-management and restoration alternatives are under consideration, one of which involves opening and reconnecting Wapato Lake’s outlet to allow flow in and out of the lake to Wapato Creek and downstream to the Tualatin River. The effects of this and other alternatives are being evaluated, partly through a detailed examination of the lake’s water budget. The water budget for the lake during 2011–13 was quantified by the U.S. Geological Survey in partnership with U.S. Fish and Wildlife Service and others. Results were incorporated in a spreadsheet-based Water Management Scenario Tool (WMST) for Wapato Lake, which predicts the effects of various management actions on daily lake level and potential habitat areas for waterfowl or other target species. Incorporating the effects of a hypothetical open outlet between the lake and the downstream river network in the WMST was accomplished by using a hydraulic model to simulate the flow-exchange rate between Wapato Lake and Wapato Creek over a wide range of lake levels and downstream river conditions. A Hydraulic Engineering Center-River Analysis System (HEC-RAS) one-dimensional unsteady flow model was constructed and calibrated for Wapato Creek and part of the Tualatin River using data from October 2011 to April 2013, and then was used to simulate daily lake/creek exchange flows in water years 1992–2014 under hypothetically constant lake levels. Results were used to populate a table of lake/creek flow-exchange rates for use in the WMST; a dynamic link between the WMST and HEC-RAS was unrealistic because it would require hundreds of calls to HEC-RAS and result in long run times for a single water-year’s WMST calculations with daily time steps. Predictions of daily outlet flows from the WMST were checked against HEC-RAS simulated flows under daily varying lake levels to ensure that the timing and magnitude of lake/creek exchange flows used by the WMST were consistent with those of the hydraulic model. Two scenarios were tested with a hypothetical open lake outlet to show how the WMST could be used to inform restoration planning—one scenario used a year-round open lake outlet, and the other scenario closed that outlet for part of the high-water winter season. Results showed that flows in and out of a year-round open lake outlet would dominate the lake’s water budget and produce water depths during winter and through mid-summer that might be too deep to support waterbird species that require shallow water. Closing the lake outlet during large winter storms and high-water conditions in the downstream river network would isolate the lake from surrounding rivers, keep the lake level lower, and retain substantially more shallow-water areas. Because of the ease with which management alternatives can be evaluated, a water-budget spreadsheet tool such as the WMST has been a valuable part of an analysis of potential water-management and restoration alternatives for Wapato Lake National Wildlife Refuge.

Oregon↗
USGS2020

Modeling a 2- and 4-foot drawdown in the Link River to Keno Dam reach of the upper Klamath River, south-central Oregon

Executive Summary The most upstream, pooled reach of the Klamath River in south-central Oregon, from Link River mouth to Keno Dam (Link-Keno), has a water-surface elevation that remains relatively constant throughout the year. Two model scenarios, using an existing two-dimensional hydrodynamic and water-quality model (CE-QUAL-W2), were constructed to examine the effects of lowering the water-surface elevation by 2 and 4 feet (ft) (0.61 and 1.2 meters) throughout an entire calendar year to mimic some of the potential effects of removal or modification of Keno Dam. Model results for these drawdown scenarios were analyzed for changes in velocity, travel time, water temperature, total dissolved solids, inorganic suspended sediment, nutrients, organic matter, chlorophyll a , and dissolved oxygen, compared to the base-case model. The model used in this study had been previously calibrated with the presence of aquatic plants (macrophytes). However, most model analyses were completed for model runs where macrophytes were “turned off” because the species, abundance, and distribution of macrophytes in a lowered-water scenario were all highly uncertain. For comparison, a few model scenario runs were completed with macrophytes enabled within the model. Findings from this study include the following: Modeled water velocity increased and travel time decreased substantially throughout the reach with the 2- and 4-ft drawdown scenarios, with travel time roughly halved in the 4-ft scenario under unchanged flow conditions. For many water-quality constituents in the drawdown scenarios, the model showed little to no change in outflow concentrations at Keno Dam compared to the base case, which represents conditions based on year 2007 water-level elevations. However, chlorophyll a and particulate organic carbon concentrations increased in the Keno Dam outflow in summer with the drawdowns. This seemed to be related to the increased water velocity. As suspended algae and other particulate organic matter such as dead algal cells and detritus moved farther downstream with the faster velocity, there was relatively less deposition of that material in the most upstream model segments and relatively more deposition in downstream segments of the Link-Keno reach, especially in summer and autumn, the period with the largest algae bloom. Despite minor differences in the Keno Dam outflow, modeled dissolved oxygen concentration and water temperature showed some changes in the Link-Keno reach compared to the base case, with differences as large as 2 milligrams per liter and 2 degrees Celsius, respectively, computed as segment averages and daily averages. With macrophytes enabled, the model outflow in the drawdown scenarios showed a decrease in dissolved oxygen concentration in mid-summer that seemed to be related to decreased macrophyte abundance and primary production in those scenarios. The with-macrophyte model results have greater uncertainty, but the results indicated that macrophytes are likely to be an important factor in the dissolved oxygen budget under reduced water-surface elevations.

Oregon↗
USGS2018

Modeling hydrodynamics, water temperature, and water quality in Klamath Straits Drain, Oregon and California, 2012–15

Executive Summary Located southwest of Klamath Falls, Oregon, Klamath Straits Drain is a 10.1-mile-long canal that conveys water uphill and northward through the use of pumps before discharging to the Klamath River. Klamath Straits Drain traverses an area that historically encompassed Lower Klamath Lake. Currently, the Drain receives water from farmland and from parts of the Lower Klamath Lake National Wildlife Refuge. To support water-quality improvement in Klamath Straits Drain, a hydrodynamic and water-temperature model was constructed and calibrated for calendar years 2012–15 with the two-dimensional model CE-QUAL-W2 (version 4.0). Water quality was calibrated for a subset of that time, from April 1, 2012 to March 31, 2015. Flows in calendar year 2012 were within the normal range, while calendar years 2013–15 were dry years. Significant findings from this study include: In the years studied, only limited flow entered Klamath Straits Drain at the upstream Headworks (KSDH) site. Most flow entered the Drain between KSDH and the E-EE pumps near Township Road through several irrigation channels and ditches. Few data were available to describe the quality of this water for the period of study. The E-EE and F-FF pumps along Klamath Straits Drain mainly operated automatically to keep water levels relatively steady. Ten-minute flow data at streamgage 11509340, downstream of the F-FF pumps, showed high-frequency on/off switching of the F-FF pumps. Combined with daily mean flow data from the F-FF pumps, the downstream 10-minute flow data allowed estimation of 10-minute pumping rates for the F-FF pumps. Paper pump charts showed the existence of short-term variability at the E-EE pumps; however, daily pump data were used at the E-EE pump location in the model. Water temperature in Klamath Straits Drain varied from less than 5 degrees Celsius (°C) (with occasional ice cover in December–January) to greater than 20 °C in May–September. In the years studied, specific conductance was typically 250–850 microsiemens per centimeter, higher than Klamath River specific conductance (typically 100–200 microsiemens per centimeter). Increased chlorophyll a in autumn and winter, along with supersaturated oxygen concentrations, indicated algal blooms in the Drain at that time of year. The blooms were most likely diatoms, based on the timing of blooms sampled elsewhere. Total nitrogen concentration was as much as 5.5 mg/L, with most in dissolved organic and particulate forms, and lower amounts in ammonia and nitrate+nitrite. Total phosphorus concentrations were distributed between orthophophorus (at a median concentration of 0.15 mg/L) and organic and particulate forms (at a median concentration of 0.13 mg/L). Most of the organic carbon in the Klamath Straits Drain was in dissolved rather than particulate form. Newly collected water-quality data for April 1, 2012–March 31, 2015 helped provide the impetus for this modeling study. However, a lack of some data still hindered the construction and calibration of this model. The model would benefit from additional data to describe water-quality boundary conditions, water-quality calibration data upstream of the F-FF pumps, short-term E-EE pump operations, and channel bathymetry in the reach between Highway 97 and the confluence with the Klamath River. Klamath River water mixed upstream into the Klamath Straits Drain, up to the Klamath Straits Drain F-FF pumps at Highway 97, when the F-FF pumps were not operating for periods of hours to days. The F-FF pumps were off for many days during this study, especially during dry years. The boundary between Klamath Straits Drain and the Klamath River was best modeled with an external head condition, which allows exchange of water between the river and the drain in both directions, upstream and downstream. Currently there is a flow gage, water-quality monitor, and a water-quality sampling site located downstream of the F-FF pumps, in the reach where Klamath Straits Drain water can mix with Klamath River water. To sample solely Klamath Straits Drain water, water samples would need to be collected only when the F-FF pumps are actively pumping. Alternately, the sampling location could be moved upstream of the pumps. Interpretation and use of historical water-quality data at the Klamath Straits Drain at Highway 97 site should be done in conjunction with information on pump activity to help inform whether mixing with Klamath River water may have occurred. Total 2014 (a dry year) phosphorus loads from the Drain to the Klamath River were lower and closer to total maximum daily load (TMDL) allocations, as compared to 2013, a year with greater flow and pumping. Modeled travel time through the Klamath Straits Drain, from Headworks to its confluence with the Klamath River, ranged from approximately 24 hours at high flow to 16 days or more, depending on how many days the pumps were turned off. The longer travel times are sufficient for important water-quality transformations, such as algal growth and organic-matter decomposition. This newly constructed model of the Klamath Straits Drain simulates flow, water levels, water temperature, and water quality with acceptable accuracy but with certain data limitations. This model should prove useful in evaluating potential strategies for flow and water-quality management and restoration.

California, Oregon↗
USGS2017

Characteristics of dissolved organic matter in the Upper Klamath River, Lost River, and Klamath Straits Drain, Oregon and California

Concentrations of particulate organic carbon (POC) and dissolved organic carbon (DOC), which together comprise total organic carbon, were measured in this reconnaissance study at sampling sites in the Upper Klamath River, Lost River, and Klamath Straits Drain in 2013–16. Optical absorbance and fluorescence properties of dissolved organic matter (DOM), which contains DOC, also were analyzed. Parallel factor analysis was used to decompose the optical fluorescence data into five key components for all samples. Principal component analysis (PCA) was used to investigate differences in DOM source and processing among sites. At all sites in this study, average DOC concentrations were higher than average POC concentrations. The highest DOC concentrations were at sites in the Klamath Straits Drain and at Pump Plant D. Evaluation of optical properties indicated that Klamath Straits Drain DOM had a refractory, terrestrial source, likely extracted from the interaction of this water with wetland peats and irrigated soils. Pump Plant D DOM exhibited more labile characteristics, which could, for instance, indicate contributions from algal or microbial exudates. The samples from Klamath River also had more microbial or algal derived material, as indicated by PCA analysis of the optical properties. Most sites, except Pump Plant D, showed a linear relation between fluorescent dissolved organic matter (fDOM) and DOC concentration, indicating these measurements are highly correlated (R 2 =0.84), and thus a continuous fDOM probe could be used to estimate DOC loads from these sites.

California, Oregon↗
USGS2016

Modeling water quality, temperature, and flow in Link River, south-central Oregon

The 2.1-km (1.3-mi) Link River connects Upper Klamath Lake to the Klamath River in south-central Oregon. A CE-QUAL-W2 flow and water-quality model of Link River was developed to provide a connection between an existing model of the upper Klamath River and any existing or future models of Upper Klamath Lake. Water-quality sampling at six locations in Link River was done during 2013–15 to support model development and to provide a better understanding of instream biogeochemical processes. The short reach and high velocities in Link River resulted in fast travel times and limited water-quality transformations, except for dissolved oxygen. Reaeration through the reach, especially at the falls in Link River, was particularly important in moderating dissolved oxygen concentrations that at times entered the reach at Link River Dam with marked supersaturation or subsaturation. This reaeration resulted in concentrations closer to saturation downstream at the mouth of Link River.

Oregon↗
USGS2014

Water-quality modeling of Klamath Straits Drain recirculation, a Klamath River wetland, and 2011 conditions for the Link River to Keno Dam reach of the Klamath River, Oregon

The upper Klamath River and adjacent Lost River are interconnected basins in south-central Oregon and northern California. Both basins have impaired water quality with Total Maximum Daily Loads (TMDLs) in progress or approved. In cooperation with the Bureau of Reclamation, the U.S. Geological Survey (USGS) and Watercourse Engineering, Inc., have conducted modeling and research to inform management of these basins for multiple purposes, including agriculture, endangered species protection, wildlife refuges, and adjacent and downstream water users. A water-quality and hydrodynamic model (CE-QUAL-W2) of the Link River to Keno Dam reach of the Klamath River for 2006–09 is one of the tools used in this work. The model can simulate stage, flow, water velocity, ice cover, water temperature, specific conductance, suspended sediment, nutrients, organic matter in bed sediment and the water column, three algal groups, three macrophyte groups, dissolved oxygen, and pH. This report documents two model scenarios and a test of the existing model applied to year 2011, which had exceptional water quality. The first scenario examined the water-quality effects of recirculating Klamath Straits Drain flows into the Ady Canal, to conserve water and to decrease flows from the Klamath Straits Drain to the Klamath River. The second scenario explicitly incorporated a 2.73×10 6 m 2 (675 acre) off-channel connected wetland into the CE-QUAL-W2 framework, with the wetland operating from May 1 through October 31. The wetland represented a managed treatment feature to decrease organic matter loads and process nutrients. Finally, the summer of 2011 showed substantially higher dissolved-oxygen concentrations in the Link-Keno reach than in other recent years, so the Link-Keno model (originally developed for 2006–09) was run with 2011 data as a test of model parameters and rates and to develop insights regarding the reasons for the improved water-quality conditions.

Oregon↗
USGS2013

Review of revised Klamath River Total Maximum Daily Load models from Link River Dam to Keno Dam, Oregon

Flow and water-quality models are being used to support the development of Total Maximum Daily Load (TMDL) plans for the Klamath River downstream of Upper Klamath Lake (UKL) in south-central Oregon. For riverine reaches, the RMA-2 and RMA-11 models were used, whereas the CE-QUAL-W2 model was used to simulate pooled reaches. The U.S. Geological Survey (USGS) was asked to review the most upstream of these models, from Link River Dam at the outlet of UKL downstream through the first pooled reach of the Klamath River from Lake Ewauna to Keno Dam. Previous versions of these models were reviewed in 2009 by USGS. Since that time, important revisions were made to correct several problems and address other issues. This review documents an assessment of the revised models, with emphasis on the model revisions and any remaining issues. The primary focus of this review is the 19.7-mile Lake Ewauna to Keno Dam reach of the Klamath River that was simulated with the CE-QUAL-W2 model. Water spends far more time in the Lake Ewauna to Keno Dam reach than in the 1-mile Link River reach that connects UKL to the Klamath River, and most of the critical reactions affecting water quality upstream of Keno Dam occur in that pooled reach. This model review includes assessments of years 2000 and 2002 current conditions scenarios, which were used to calibrate the model, as well as a natural conditions scenario that was used as the reference condition for the TMDL and was based on the 2000 flow conditions. The natural conditions scenario included the removal of Keno Dam, restoration of the Keno reef (a shallow spot that was removed when the dam was built), removal of all point-source inputs, and derivation of upstream boundary water-quality inputs from a previously developed UKL TMDL model. This review examined the details of the models, including model algorithms, parameter values, and boundary conditions; the review did not assess the draft Klamath River TMDL or the TMDL allocations. Attention to the details of a model is one of the best ways to identify potential problems, correct them if possible, and begin to assess the magnitude of potential model errors and uncertainty. Model users need to determine the level of acceptable uncertainty associated with their objectives, identify all sources of potential uncertainty (model uncertainty, data uncertainty, etc.), and assess their approach and results accordingly. In the draft Klamath River TMDL, the Oregon Department of Environmental Quality identified the upstream boundary conditions as the largest source of uncertainty for both the current and natural conditions scenarios, not the model algorithms or choice of model parameters. We agree that the upstream boundary conditions are one of the largest, if not the largest, source of model uncertainty; therefore, the derivation of upstream boundary conditions may be more important to the TMDL than some other model-related issues identified in this review. The revised models contain a number of changes, some of which were done to solve small problems and are largely inconsequential to model results, but others of which are important and affect model predictions of instream concentrations. A consistent version of the model is now applied to all scenarios, and an error in the source code was corrected that had inadvertently discarded 20 percent of the incoming solar radiation in the original model. The baseline light-extinction coefficient for water was decreased and set to a consistent and defensible value across all models of reservoir reaches. Inconsistencies among the values of certain parameters in the original models, such as the ammonia nitrification rate and the decomposition rates of organic matter, have been eliminated, although the reasoning behind the final selections was not documented. The dependence of the rate of sediment oxygen demand (SOD) on temperature was modified such that the SOD rate was substantially decreased at temperatures less than 20°C, causing the model to predict higher dissolved oxygen (DO) concentrations in spring, autumn, and winter. Although that change to the temperature dependence function was done to make the function more similar to the model’s default, this change was not accompanied by any documentation of recalibration or sensitivity exercises. The maximum SOD rate for the 2002 current conditions scenario was decreased from 3.0 grams per square meter per day (g/m 2 /d) in the original model to 2.0 g/m 2 /d in the revised model, a considerable adjustment that appears to have been needed to offset effects of a change to another variable (O2LIM) that would have resulted in a substantial increase in the effective SOD rate for 2002. A 50-percent decrease in the SOD rate over a 2-year period, however, is not likely to be mirrored by field measurements, so this change may be compensating for some process that is not represented correctly in the DO budget for the current conditions scenarios. Several important changes were made to the natural conditions scenario. First, the elevation of the Keno reef was corrected; the elevation specified in the original model was 1 foot too high, which affected the volume of the pooled reach and the travel time through it. The most important changes to this scenario were to the upstream boundary inputs of organic matter and algae, which affect incoming fluxes of nitrogen and phosphorus. Algal biomass inputs were increased by approximately 60 percent during summer because of a change in the way those inputs were derived from results of the UKL TMDL model. Non-algal organic matter inputs were decreased, particularly in summer to correct a problem attributed to double-counting of phosphorus in the original inputs. The distribution of non-algal organic matter was changed from 20 percent dissolved in the original model to 90 percent dissolved in the revised model in response to review comments and published data. The overall sum of algal biomass and non-living organic matter was decreased, which resulted in lower inputs of total phosphorus and nitrogen. Total phosphorus inputs were less than 0.03 mg/L, and although the inputs were derived from selected results of the UKL TMDL model, these concentrations seem too low to be representative of a historically eutrophic system surrounded by extensive wetlands, peat soils, and a groundwater system high in phosphorus. The draft TMDL states that the upstream boundary conditions are the greatest source of uncertainty, greater than any uncertainty associated with the models. Efforts to improve existing models of algal growth and nutrient cycling in UKL, therefore, would provide a substantial benefit to downstream modeling efforts on the Klamath River. Although many improvements were made in revising the Klamath River TMDL models, some issues and uncertainties remain. Several errors in the model source code remain, but do not affect model results for this application as long as certain options and rates are not changed; future users of these models should be aware of these issues. Although the distribution of dissolved and particulate organic matter was modified for the natural conditions scenario, that distribution was not changed for the current conditions scenarios. Recent data on that distribution and the likely rates of organic matter decomposition could be used to improve these models in the future. Nitrate predictions at Keno (Highway 66) still are too high for the current conditions scenarios; future efforts should re-evaluate the model’s denitrification rates and the release rate of ammonia from anoxic sediments. Possibly the most important of the remaining issues are tied to the two-state (healthy/unhealthy) hypothesis for the algae population that was coded into the model. Some of the rates and conversion functions could be refined to make them more acceptable; currently, the published literature does not support the concept of moderately low dissolved-oxygen concentrations as a stressor of algae in the ranges used by the model. More research is needed before these algorithms can be truly tested. The algorithms currently appear to help the model fit the patterns in the available data, and that is useful and perhaps sufficient for some purposes, but those algorithms are not truly predictive or reliable for certain purposes until they can be tested through well-designed experiments and research. In summary, the TMDL models used to simulate Link and Klamath Rivers from Link River Dam to Keno Dam were revised to fix several problems and address various issues. The resulting models are an improvement over those that were reviewed by USGS in 2009, and represent a useful advance in the simulation of a complex system that is difficult to model. However, several issues remain that cause increased uncertainty in the model results. Depending on the objectives of the modeling, now or in the future, these remaining issues could be more or less important. For the Klamath River TMDL, the upstream boundary conditions may be a larger source of uncertainty than the concerns with model algorithms and model parameters identified in this review. Efforts to re-evaluate the available models of algal growth and nutrient cycling in UKL would be highly beneficial to downstream modeling efforts in the Klamath River.

Oregon;California↗
USGS2013

Development of CE-QUAL-W2 models for the Middle Fork Willamette and South Santiam Rivers, Oregon

Hydrodynamic (CE-QUAL-W2) models of Hills Creek Lake (HCL), Lookout Point Lake (LOP), and Dexter Lake (DEX) on the Middle Fork Willamette River (MFWR), and models of Green Peter Lake and Foster Lake on the South Santiam River systems in western Oregon were updated and recalibrated for a wide range of flow and meteorological conditions. These CE-QUAL-W2 models originally were developed by West Consultants, Inc., for the U.S. Army Corps of Engineers. This study by the U.S. Geological Survey included a reassessment of the models’ calibration in more recent years—2002, 2006, 2008, and 2011—categorized respectively as low, normal, high, and extremely high flow calendar years. These years incorporated current dam-operation practices and more available data than the time period used in the original calibration. Modeled water temperatures downstream of both HCL and LOP-DEX on the MFWR were within an average of 0.68 degree Celsius (°C) of measured values; modeled temperatures downstream of Foster Dam on the South Santiam River were within an average of 0.65°C of measured values. A new CE-QUAL-W2 model was developed and calibrated for the riverine MFWR reach between Hills Creek Dam and the head of LOP, allowing an evaluation of the flow and temperature conditions in the entire MFWR system from HCL to Dexter Dam. The complex bathymetry and long residence time of HCL, combined with the relatively deep location of the power and regulating outlet structures at Hills Creek Dam, led to a HCL model that was highly sensitive to several outlet and geometric parameters related to dam structures (STR TOP, STR BOT, STR WIDTH). Release temperatures from HCL were important and often persisted downstream as they were incorporated in the MFWR model and the LOP-DEX model (downstream of MFWR). The models tended to underpredict the measured temperature of water releases from Dexter Dam during the late-September-through-December drawdown period in 2002, and again (to a lesser extent) in 2011, but simulations were much more accurate in 2006 and 2008. This episodic model bias may have been a result of hot, dry conditions; lower lake elevations; and earlier drawdown at both HCL and LOP in 2002. These dry conditions in 2002 may have contradicted assumptions inherent in the estimation of certain model inputs, such as unmeasured inflows and water temperatures, which may respond differently during dry years than during normal and wet years. This report documents the development and calibration of new and revised flow and water-temperature models for riverine and reservoir reaches in the Middle Fork Willamette River and South Santiam River systems. Methods and model parameter values were established for the accurate simulation of flows and temperatures in these systems under current conditions. By extension, these models should be able to accurately simulate flows and temperatures under potential future conditions in which dam operations and dam outlet structures may be changed as part of a strategy to improve habitat, fish passage, and temperature conditions for endangered fish.

Oregon↗
USGS2013

Macrophyte and pH buffering updates to the Klamath River water-quality model upstream of Keno Dam, Oregon

A hydrodynamic, water temperature, and water-quality model of the Link River to Keno Dam reach of the upper Klamath River was updated to account for macrophytes and enhanced pH buffering from dissolved organic matter, ammonia, and orthophosphorus. Macrophytes had been observed in this reach by field personnel, so macrophyte field data were collected in summer and fall (June-October) 2011 to provide a dataset to guide the inclusion of macrophytes in the model. Three types of macrophytes were most common: pondweed (Potamogeton species), coontail (Ceratophyllum demersum), and common waterweed (Elodea canadensis). Pondweed was found throughout the Link River to Keno Dam reach in early summer with densities declining by mid-summer and fall. Coontail and common waterweed were more common in the lower reach near Keno Dam and were at highest density in summer. All species were most dense in shallow water (less than 2 meters deep) near shore. The highest estimated dry weight biomass for any sample during the study was 202 grams per square meter for coontail in August. Guided by field results, three macrophyte groups were incorporated into the CE-QUAL-W2 model for calendar years 2006-09. The CE-QUAL-W2 model code was adjusted to allow the user to initialize macrophyte populations spatially across the model grid. The default CE-QUAL-W2 model includes pH buffering by carbonates, but does not include pH buffering by organic matter, ammonia, or orthophosphorus. These three constituents, especially dissolved organic matter, are present in the upper Klamath River at concentrations that provide substantial pH buffering capacity. In this study, CE-QUAL-W2 was updated to include this enhanced buffering capacity in the simulation of pH. Acid dissociation constants for ammonium and phosphoric acid were taken from the literature. For dissolved organic matter, the number of organic acid groups and each group's acid dissociation constant (K a ) and site density (moles of sites per mole of carbon) were derived by fitting a theoretical buffering response to measured upper Klamath River alkalinity titration curves. The organic matter buffering in the Klamath River was modeled with two monoprotic organic acids: carboxylic acids with a mean pK a of 5.584 and site density of 0.1925, and phenolic organic acids with a mean pK a of 9.594 and site density of 0.6466. Total inorganic carbon concentrations in the model boundary inputs were recalculated based on the new buffering equations. CE-QUAL-W2 was also adjusted to allow the simulation of nonconservative alkalinity caused by nitrification, denitrification, photosynthesis, and respiration. The Klamath River model was recalibrated after the macrophyte and pH buffering updates producing improved predictions for pH, dissolved oxygen, and particulate carbon.

Oregon↗
USGS2013

Modeling the Water - Quality Effects of Changes to the Klamath River Upstream of Keno Dam, Oregon

The Link River to Keno Dam (Link-Keno) reach of the Klamath River, Oregon, generally has periods of water-quality impairment during summer, including low dissolved oxygen, elevated concentrations of ammonia and algae, and high pH. Efforts are underway to improve water quality in this reach through a Total Maximum Daily Load (TMDL) program and other management and operational actions. To assist in planning, a hydrodynamic and water-quality model was used in this study to provide insight about how various actions could affect water quality in the reach. These model scenarios used a previously developed and calibrated CE-QUAL-W2 model of the Link-Keno reach developed by the U.S. Geological Survey (USGS), Watercourse Engineering Inc., and the Bureau of Reclamation for calendar years 2006-09 (referred to as the "USGS model" in this report). Another model of the same river reach was previously developed by Tetra Tech, Inc. and the Oregon Department of Environmental Quality for years 2000 and 2002 and was used in the TMDL process; that model is referred to as the "TMDL model" in this report. This report includes scenarios that (1) assess the effect of TMDL allocations on water quality, (2) provide insight on certain aspects of the TMDL model, (3) assess various methods to improve water quality in this reach, and (4) examine possible water-quality effects of a future warmer climate. Results presented in this report for the first 5 scenarios supersede or augment those that were previously published (scenarios 1 and 2 in Sullivan and others [2011], 3 through 5 in Sullivan and others [2012]); those previous results are still valid, but the results for those scenarios in this report are more current.

Oregon↗
USGS2012

Dissolved oxygen analysis, TMDL model comparison, and particulate matter shunting—Preliminary results from three model scenarios for the Klamath River upstream of Keno Dam, Oregon

Efforts are underway to identify actions that would improve water quality in the Link River to Keno Dam reach of the Upper Klamath River in south-central Oregon. To provide further insight into water-quality improvement options, three scenarios were developed, run, and analyzed using previously calibrated CE-QUAL-W2 hydrodynamic and water-quality models. Additional scenarios are under development as part of this ongoing study. Most of these scenarios evaluate changes relative to a "current conditions" model, but in some cases a "natural conditions" model was used that simulated the reach without the effect of point and nonpoint sources and set Upper Klamath Lake at its Total Maximum Daily Load (TMDL) targets. These scenarios were simulated using a model developed by the U.S. Geological Survey (USGS) and Watercourse Engineering, Inc. for the years 2006–09, referred to here as the "USGS model." Another model of the reach was developed by Tetra Tech, Inc. for years 2000 and 2002 to support the Klamath River TMDL process; that model is referred to here as the "TMDL model." The three scenarios described in this report included (1) an analysis of whether this reach of the Upper Klamath River would be in compliance with dissolved oxygen standards if sources met TMDL allocations, (2) an application of more recent datasets to the TMDL model with comparison to results from the USGS model, and (3) an examination of the effect on dissolved oxygen in the Klamath River if particulate material were stopped from entering Klamath Project diversion canals. Updates and modifications to the USGS model are in progress, so in the future these scenarios will be reanalyzed with the updated model and the interim results presented here will be superseded. Significant findings from this phase of the investigation include: * The TMDL analysis used depth-averaged dissolved oxygen concentrations from model output for comparison with dissolved oxygen standards. The Oregon dissolved oxygen standards do not specify whether the numeric criteria are based on depth-averaged dissolved oxygen concentration; this was an interpretation of the standards rule by the Oregon Department of Environmental Quality (ODEQ). In this study, both depth-averaged and volume-averaged dissolved oxygen concentrations were calculated from model output. Results showed that modeled depth-averaged concentrations typically were lower than volume-averaged dissolved oxygen concentrations because depth-averaging gives a higher weight to small volume areas near the channel bottom that often have lower dissolved oxygen concentrations. Results from model scenarios in this study are reported using volume-averaged dissolved oxygen concentrations. * Under all scenarios analyzed, violations of the dissolved oxygen standard occurred most often in summer. Of the three dissolved oxygen criteria that must be met, the 30-day standard was violated most frequently. Under the base case (current conditions), fewer violations occurred in the upstream part of the reach. More violations occurred in the down-stream direction, due in part to oxygen demand from the decay of algae and organic matter from Link River and other inflows. * A condition in which Upper Klamath Lake and its Link River outflow achieved Upper Klamath Lake TMDL water-quality targets was most effective in reducing the number of violations of the dissolved oxygen standard in the Link River to Keno Dam reach of the Klamath River. The condition in which point and nonpoint sources within the Link River to Keno Dam reach met Klamath River TMDL allocations had no effect on dissolved oxygen compliance in some locations and a small effect in others under current conditions. On the other hand, meeting TMDL allocations for nonpoint and point sources was predicted to be important in meeting dissolved oxygen criteria when Upper Klamath Lake and Link River also met Upper Klamath TMDL water-quality targets. * The location of greatest dissolved oxygen improvement from nutrient and organic matter reductions was downstream from point and nonpoint source inflows because time and distance are required for decay to occur and for oxygen demand to be exerted. * After assessing compliance with dissolved oxygen standards at all 102 model segments in the Link River to Keno Dam reach, it was determined that the seven locations used by ODEQ appear to be a representative subset of the reach for dissolved oxygen analysis. * The USGS and TMDL models were qualitatively compared by running both models for the 2006–09 period but preserving the essential characteristics of each, such as organic matter partitioning, bathymetric representation, and parameter rates. The analysis revealed that some constituents were not greatly affected by the differing algorithms, rates, and assumptions in the two models. Conversely, other constituents, especially organic matter, were simulated differently by the two models. Organic matter in this river system is best represented by a mixture of relatively labile particulate material and a substantial concentration of refractory dissolved material. In addition, the use of a first-order sediment oxygen demand, as in the USGS model, helps to capture the seasonal and dynamic effect of settled organic and algal material. * Simulation of shunting (diverting) particulate material away from the intake of four Klamath Project diversion canals, so that the material stayed in the river and out of the Project area, caused higher concentrations of particulate material to occur in the river. In all cases modeled, the increase in in-river particulate material also produced decreased dissolved oxygen concentrations and an increase in the number of days when dissolved oxygen standards were violated. * If particulate material were shunted back into the river at the Klamath Project diversion canals, less organic matter and nutrients would be taken into the Klamath Project area and the Lost River basin, resulting in return flows to the Klamath River via Lost River Diversion Channel that may have reduced nutrient concentrations. Model scenarios bracketing potential end-member nutrient concentrations showed that the composition of the return flows had little to no effect on dissolved oxygen compliance under simulated conditions.

Oregon↗
USGS2012

Simulating potential structural and operational changes for Detroit Dam on the North Santiam River, Oregon, for downstream temperature management

Detroit Dam was constructed in 1953 on the North Santiam River in western Oregon and resulted in the formation of Detroit Lake. With a full-pool storage volume of 455,100 acre-feet and a dam height of 463 feet, Detroit Lake is one of the largest and most important reservoirs in the Willamette River basin in terms of power generation, recreation, and water storage and releases. The U.S. Army Corps of Engineers operates Detroit Dam as part of a system of 13 reservoirs in the Willamette Project to meet multiple goals, which include flood-damage protection, power generation, downstream navigation, recreation, and irrigation. A distinct cycle in water temperature occurs in Detroit Lake as spring and summer heating through solar radiation creates a warm layer of water near the surface and isolates cold water below. Controlling the temperature of releases from Detroit Dam, therefore, is highly dependent on the location, characteristics, and usage of the dam's outlet structures. Prior to operational changes in 2007, Detroit Dam had a well-documented effect on downstream water temperature that was problematic for endangered salmonid fish species, releasing water that was too cold in midsummer and too warm in autumn. This unnatural seasonal temperature pattern caused problems in the timing of fish migration, spawning, and emergence. In this study, an existing calibrated 2-dimensional hydrodynamic water-quality model [CE-QUAL-W2] of Detroit Lake was used to determine how changes in dam operation or changes to the structural release points of Detroit Dam might affect downstream water temperatures under a range of historical hydrologic and meteorological conditions. The results from a subset of the Detroit Lake model scenarios then were used as forcing conditions for downstream CE-QUAL-W2 models of Big Cliff Reservoir (the small reregulating reservoir just downstream of Detroit Dam) and the North Santiam and Santiam Rivers. Many combinations of environmental, operational, and structural options were explored with the model scenarios. Multiple downstream temperature targets were used along with three sets of environmental forcing conditions representing cool/wet, normal, and hot/dry conditions. Five structural options at Detroit Dam were modeled, including the use of existing outlets, one hypothetical variable-elevation outlet such as a sliding gate, a hypothetical combination of a floating outlet and a fixed-elevation outlet, and a hypothetical combination of a floating outlet and a sliding gate. Finally, 14 sets of operational guidelines for Detroit Dam were explored to gain an understanding of the effects of imposing different downstream minimum streamflows, imposing minimum outflow rules to specific outlets, and managing the level of the lake with different timelines through the year. Selected subsets of these combinations of operational and structural scenarios were run through the downstream models of Big Cliff Reservoir and the North Santiam and Santiam Rivers to explore how hypothetical changes at Detroit Dam might provide improved temperatures for endangered salmonids downstream of the Detroit-Big Cliff Dam complex. Conclusions that can be drawn from these model scenarios include: *The water-temperature targets set by the U.S. Army Corps of Engineers for releases from Detroit Dam can be met through a combination of new dam outlets or a delayed drawdown of the lake in autumn. *Spring and summer dam operations greatly affect the available release temperatures and operational flexibility later in the autumn. Releasing warm water during midsummer tends to keep more cool water available for release in autumn. *The ability to meet downstream temperature targets during spring depends on the characteristics of the available outlets. Under existing conditions, although warm water sometimes is present at the lake surface in spring and early summer, such water may not be available for release if the lake level is either well below or well above the spillway crest. *Managing lake releases to meet downstream temperature targets depends on having outlet structures that can access both (warm) lake surface water and (cold) deeper lake water throughout the year. The existing outlets at Detroit Dam do not allow near-surface waters to be released during times when the lake surface level is below the spillway (spring and autumn). *Using the existing outlets at Detroit Dam, lake level management is important to the water temperature of releases because it controls the availability and depth of water at the spillway. When lake level is lowered below the spillway crest in late summer, the loss of access to warm water at the lake surface can result in abrupt changes to release temperatures. *Because the power-generation intakes (penstocks) are 166 feet below the full-pool lake level, imposing minimum power production requirements at Detroit Dam limits the amount of warm surface water that can be expelled from the lake in midsummer, thereby postponing and amplifying warm outflows from Detroit Lake into the autumn spawning season. *Likewise, imposing minimum power production requirements at Detroit Dam in autumn can limit the amount of cool hypolimnetic water that is released from the lake, thereby limiting cool outflows from Detroit Lake during the autumn spawning season. *Model simulations indicate that a delayed drawdown of Detroit Lake in autumn would result in better control over release temperatures in the immediate downstream vicinity of Big Cliff Dam, but the reduced outflows necessary to retain more water in the lake in late summer are more susceptible to rapid heating downstream. *Compared to the existing outlets at Detroit Dam, floating or sliding-gate outlet structures can provide greater control over release temperatures because they provide better access to warm water at the lake surface and cooler water at depth. These conclusions can be grouped into several common themes. First, optimal and flexible management and achievement of downstream temperature goals requires that releases of warm water near the surface of the lake and cold water below the thermocline are both possible with the available dam outlets during spring, summer, and autumn. This constraint can be met to some extent with existing outlets, but only if access to the spillway is extended into autumn by keeping the lake level higher than called for by the current rule curve (the typical target water-surface elevation throughout the year). If new outlets are considered, a variable-elevation outlet such as a sliding gate structure, or a floating outlet in combination with a fixed-elevation outlet at sufficient depth to access cold water, is likely to work well in terms of accessing a range of water temperatures and achieving downstream temperature targets. Furthermore, model results indicate that it is important to release warm water from near the lake surface during midsummer. If not released downstream, the warm water will build up at the top of the lake as a result of solar energy inputs and the thermocline will deepen, potentially causing warm water to reach the depth of deeper fixed-elevation outlets in autumn, particularly when the lake level is drawn down to make room for flood storage. Delaying the drawdown in autumn can help to keep the thermocline above such outlets and preserve access to cold water. Although it is important to generate hydropower at Detroit Dam, minimum power-production requirements limit the ability of dam operators to meet downstream temperature targets with existing outlet structures. The location of the power penstocks below the thermocline in spring and most of summer causes the release of more cool water during summer than is optimal. Reducing the power-production constraint allows the temperature target to be met more frequently, but at the cost of less power generation. Finally, running the Detroit Dam, Big Cliff Dam, and North Santiam and Santiam River models in series allows dam operators to evaluate how different operational strategies or combinations of new dam outlets might affect downstream temperatures for many miles of critical endangered salmonid habitat. Temperatures can change quickly in these downstream reaches as the river exchanges heat with its surroundings, and heating or cooling of 6 degrees Celsius is not unusual in the 40–50 miles downstream of Big Cliff Dam. The results published in this report supersede preliminary results published in U.S. Geological Survey Open-File Report 2011-1268 (Buccola and Rounds, 2011). Those preliminary results are still valid, but the results in this report are more current and comprehensive.

Oregon↗
USGS2011

Modeling hydrodynamics, water temperature, and water quality in the Klamath River upstream of Keno Dam, Oregon, 2006-09

A hydrodynamic, water temperature, and water-quality model was constructed for a 20-mile reach of the Klamath River downstream of Upper Klamath Lake, from Link River to Keno Dam, for calendar years 2006-09. The two-dimensional, laterally averaged model CE-QUAL-W2 was used to simulate water velocity, ice cover, water temperature, specific conductance, dissolved and suspended solids, dissolved oxygen, total nitrogen, ammonia, nitrate, total phosphorus, orthophosphate, dissolved and particulate organic matter, and three algal groups. The Link-Keno model successfully simulated the most important spatial and temporal patterns in the measured data for this 4-year time period. The model calibration process provided critical insights into water-quality processes and the nature of those inputs and processes that drive water quality in this reach. The model was used not only to reproduce and better understand water-quality conditions that occurred in 2006-09, but also to test several load-reduction scenarios that have implications for future water-resources management in the river basin. The model construction and calibration process provided results concerning water quality and transport in the Link-Keno reach of the Klamath River, ranging from interesting circulation patterns in the Lake Ewauna area to the nature and importance of organic matter and algae. These insights and results include: * Modeled segment-average water velocities ranged from near 0.0 to 3.0 ft/s in 2006 through 2009. Travel time through the model reach was about 4 days at 2,000 ft 3 /s and 12 days at 700 ft3/s flow. Flow direction was aligned with the upstream-downstream channel axis for most of the Link-Keno reach, except for Lake Ewauna. Wind effects were pronounced at Lake Ewauna during low-flow conditions, often with circulation in the form of a gyre that rotated in a clockwise direction when winds were towards the southeast and in a counterclockwise direction when winds were towards the northwest. * Water temperatures ranged from near freezing in winter to near 30 degrees C at some locations and periods in summer; seasonal water temperature patterns were similar at the inflow and outflow. Although vertical temperature stratification was not present at most times and locations, weak stratification could persist for periods up to 1-2 weeks, especially in the downstream parts of the reach. Thermal stratification was important in controlling vertical variations in water quality. * The specific conductance, and thus density, of tributaries within the reach usually was higher than that of the river itself, so that inflows tended to sink below the river surface. This was especially notable for inflows from the Klamath Straits Drain, which tended to sink to the bottom of the Klamath River at its confluence and not mix vertically for several miles downstream. * The model was able to capture most of the seasonal changes in the algal population by modeling that population with three algal groups: blue-green algae, diatoms, and other algae. The blooms of blue-green algae, consisting mostly of Aphanizomenon flos aquae that entered from Upper Klamath Lake, were dominant, dwarfing the populations of the other two algae groups in summer. A large part of the blue-green algae population that entered this reach from upstream tended to settle out, die, and decompose, especially in the upper part of the Link-Keno reach. Diatoms reached a maximum in spring and other algae in midsummer. * Organic matter, occurring in both dissolved and particulate forms, was critical to the water quality of this reach of the Klamath River, and was strongly tied to nutrient and dissolved-oxygen dynamics. Dissolved and particulate organic matter were subdivided into labile (quickly decaying) and refractory (slowing decaying) groups for modeling purposes. The particulate matter in summer, consisting largely of dead blue-green algae, decayed quickly. Consequently, this particulate matt

Oregon↗
USGS2010

Controls on biochemical oxygen demand in the upper Klamath River, Oregon

A series of 30-day biochemical oxygen demand (BOD) experiments were conducted on water column samples from a reach of the upper Klamath River that experiences hypoxia and anoxia in summer. Samples were incubated with added nitrification inhibitor to measure carbonaceous BOD (CBOD), untreated to measure total BOD, which included demand from nitrogenous BOD (NBOD), and coarse-filtered to examine the effect of removing large particulate matter. All BOD data were fit well with a two-group model, so named because it considered contributions from both labile and refractory pools of carbon: BOD t = a 1 (1 − e − a 0 t ) + a 2 t . Site-average labile first-order decay rates a 0 ranged from 0.15 to 0.22/day for CBOD and 0.11 to 0.29/day for BOD. Site-average values of refractory zero-order decay rates a 2 ranged from 0.13 to 0.25 mg/L/day for CBOD and 0.01 to 0.45 mg/L/day for BOD; the zero-order CBOD decay rate increased from early- to mid-summer. Values of ultimate CBOD for the labile component a 1 ranged from 5.5 to 28.8 mg/L for CBOD, and 7.6 to 30.8 mg/L for BOD. Two upstream sites had higher CBOD compared to those downstream. Maximum measured total BOD 5 and BOD 30 during the study were 26.5 and 55.4 mg/L; minimums were 4.2 and 13.6 mg/L. For most samples, the oxygen demand from the three components considered here were: labile CBOD > NBOD > refractory CBOD, though the relative importance of refractory CBOD to oxygen demand increased over time. Coarse-filtering reduced CBOD for samples with high particulate carbon and high biovolumes of Aphanizomenon flos-aquae . There was a strong positive correlation between BOD, CBOD, and the labile component of CBOD to particulate C and N, with weaker positive correlation to field pH, field dissolved oxygen, and total N. The refractory component of CBOD was not correlated to particulate matter, instead showing weak but statistically significant correlation to dissolved organic carbon, UV absorbance at 254 nm, and total N. Particulate organic matter, especially the alga A. flos-aquae , is an important component of oxygen demand in this reach of the Klamath River, though refractory dissolved organic matter would continue to exert an oxygen demand over longer time periods and as water travels downstream.

Oregon↗