Columbia River System Research Articles

Nested Benders’s decomposition of capacity-planning problems for electricity systems with hydroelectric and renewable generation

Nested Benders’s decomposition is an efficient means to solve large-scale optimization problems with a natural time sequence of decisions. This paper examines the use of the technique to decompose and solve efficiently capacity-expansion problems for electricity systems with hydroelectric and renewable generators. To this end we develop an archetypal planning model that captures key features of hydroelectric and renewable generators and apply it to a case study that is based on the Columbia River system in the northwestern United States of America. We apply standard network and within-year temporal simplifications to reduce the problem’s size. Nevertheless, the remaining problem is large-scale and we demonstrate the use of nested Benders’s decomposition to solve it. We explore refinements of the decomposition method which yield further performance improvements. Overall, we show that nested Benders’s decomposition yields good computational performance with minimal loss of model fidelity.

Irrigation intensification impacts sustainability of streamflow in the Western United States

Quantifying the interconnected impacts of climate change and irrigation on surface water flows is critical for the proactive management of our water resources and the ecosystem services they provide. Changes in streamflow across the Western U.S. have generally been attributed to an aridifying climate, but in many basins flows can also be highly impacted by irrigation. We developed a 35-year dataset consisting of streamflow, climate, irrigated area, and crop water use to quantify the effects of both climate change and irrigation water use on streamflow across 221 basins in the Colorado, Columbia, and Missouri River systems. We demonstrate that flows have been altered beyond observed climate-related changes and that many of these changes are attributable to irrigation. Further, our results indicate that increases in irrigation water use have occurred over much of the study area, a finding that contradicts government-reported irrigation statistics. Increases in crop consumption have enhanced fall and winter flows in some portions of the Upper Missouri and northern Columbia River basins, and have exacerbated climate change-induced flow declines in parts of the Colorado basin. We classify each basin’s water resources sustainability in terms of flow and irrigation trends and link irrigation-induced flow changes to irrigation infrastructure modernization and differences in basin physiographic setting. These results provide a basis for determining where modern irrigation systems benefit basin water supply, and where less efficient systems contribute to return flows and relieve ecological stress.

Simplifying capacity planning for electricity systems with hydroelectric and renewable generation

This work investigates approaches to simplify capacity planning for electricity systems with hydroelectric and renewable generators with three specific foci. First, we examine approaches to represent the efficiency of hydroelectric units. Next, we explore the effects of water-travel times and the representation of run-of-river units within cascaded hydroelectric systems. Third, we analyze the use of representative operating periods to capture electricity-system operations. We conduct these analyses using an archetypal planning models that is applied to the Columbia River system in the northwestern United States of America. We demonstrate that planning models can be simplified significantly, which improves model tractability with little loss of fidelity.

Theoretical analysis of green hydrogen from hydropower: A case study of the Northwest Columbia River system

Within the Pacific Northwest region of the United States, there is the unique opportunity to explore alternative energy management solutions of the Columbia River's multi-use hydropower system. As with various European hydropower systems that experience large variability in water runoff, but lack adequate reservoir storage capacity, the Columbia River System is a viable source for renewable hydrogen production. This paper studies the theoretical potential of green hydrogen production from excess hydropower energy from the Columbia River System. The potential surplus hydroelectric energy and hydrogen production potential from surplus energy (during March through July months) are estimated from 11 hydroelectric projects along with the Columbia River System. Results show that the system's total monthly average hydrogen production potential ranges from 2.22 × 106 to 8.96 × 106 kg H2 with the utilization of surplus energy over a historical 80 water year period (1928–2008). This study concludes that hydrogen production from spilled hydropower energy and its use in the transportation sector is a viable opportunity to lead the country towards a hydrogen economy.

Food Security: Impact of European Settlement and Infrastructure on Columbia River Salmon Migration

The 2000 km Columbia River is the longest river in the Pacific Northwest region of North America. It starts in British Columbia, Canada and flows through the states of Oregon and Washington before discharging into Pacific Ocean near Cape Disappointment. The mouth of the Columbia River is a deep water harbor and 180 km of the river can be accessed by navigation. Based on flow volume, the river is the 4th largest river in the United States. The headwaters and approximately 800 km of the Columbia River lie in Canada. Columbia Lake and the Columbia Wetlands are located in British Columbia and are the headwaters of the Columbia River. The lake has an elevation of 820 m above sea level and drains to the Pacific Ocean near Astoria, Oregon. The Ocean tides flow up river to Portland, Oregon. This paper highlights how the geological and landscape resources of the Columbia River and tributaries watershed contributed to the economic development of this historically rich region of North America. The Columbia River is one of the most biologically diverse freshwater systems in the United States. The Columbia River system with trails and cruise ships was designed to increase use of the Columbia River, to promote recreational tourism, and to create a generation of people who are willing to protect and provide environmental stewardship of the river basin resources.

Palaeogeographical reconstruction and hydrology of glacial Lake Purcell during MIS 2 and its potential impact on the Channeled Scabland, USA

Large, ice‐marginal lakes that were impounded by the maximally extended Cordilleran Ice Sheet (CIS) provided source waters for the extraordinarily large floods that formed the Channeled Scabland of Washington and Idaho, USA. However, flood flows that drained CIS meltwater and contributed to landscape evolution during later stages of deglaciation have hitherto been poorly investigated. This paper provides the first evidence for such a late deglacial floodwater source: glacial Lake Purcell (gLP). Sedimentary evidence records the northward extension of gLP from Idaho, USA into British Columbia, Canada and establishes its minimum palaeogeographical extent. Sedimentary evidence suggests that the deglacial Purcell Lobe was a capable ice dam that impounded large volumes of gLP water. A review of glacio‐isostatically affected lakes during CIS deglaciation suggests that gLP could have been subjected to tilts ranging from 0 to >1.25 m km−1. Sedimentary evidence suggests high lake plane tilts (⪆1.25 m km−1) are the most likely to have affected gLP. Using this, the palaeogeography and volume of gLP are modelled, revealing that ~116 km3 of water was susceptible to sudden drainage into the Channeled Scabland via the Columbia River system. This calculation is supported by sedimentary and geomorphic evidence compatible with energetic flood flows along the gLP drainage route and suggests gLP drained suddenly, causing significant landscape change.

Dam trout: Genetic variability in Oncorhynchus mykiss above and below barriers in three Columbia River systems prior to restoring migrational access.

Restoration of access to lost habitat for threatened and endangered fishes above currently impassable dams represents a major undertaking. Biological monitoring is critical to understand the dynamics and success of anadromous recolonization as, in the case of Oncorhynchus mykiss, anadromous steelhead populations are reconnected with their conspecific resident rainbow trout counterparts. We evaluate three river systems in the Lower Columbia River basin: the White Salmon, Sandy, and Lewis rivers that are in the process of removing and/or providing passage around existing human-made barriers in O. mykiss riverine habitat. In these instances, now isolated resident rainbow trout populations will be exposed to competition and/or genetic introgression with steelhead and vice versa. Our genetic analyses of 2,158 fish using 13 DNA microsatellite (mSAT) loci indicated that within each basin anadromous O. mykiss were genetically distinct from and significantly more diverse than their resident above-dam trout counterparts. Above long-standing natural impassable barriers, each of these watersheds also harbors unique rainbow trout gene pools with reduced levels of genetic diversity. Despite frequent releases of non-native steelhead and rainbow trout in each river, hatchery releases do not appear to have had a significant genetic effect on the population structure of O. mykiss in any of these watersheds. Simulation results suggest there is a high likelihood of identifying anadromous x resident individuals in the Lewis and White Salmon rivers, and slightly less so in the Sandy River. These genetic data are a prerequisite for informed monitoring, managing, and conserving the different life history forms during upstream recolonization when sympatry of life history forms of O. mykiss is restored.

A box model for representing estuarine physical processes in Earth system models

Appropriately treating riverine freshwater discharge into the oceans in Earth system models is a challenging problem. Commonly, the river runoff is discharged into the ocean models with zero salinity and arbitrarily distributed either horizontally or vertically over several grid cells. Those approaches entirely neglect estuarine physical processes that modify river inputs before they reach the open ocean. In order to realistically represent riverine freshwater inputs in Earth system models, a physically based Estuary Box Model (EBM) is developed to parameterize the mixing processes in estuaries. The EBM represents the estuary exchange circulation with a two-layer box structure. It takes as input the river volume flux from the land surface model and the subsurface salinity at the estuary mouth from the ocean model. It delivers the estuarine outflow salinity and net volume flux into and out of the estuary to the ocean model. An offline test of the EBM forced with observed conditions for the Columbia River system shows good agreement with observations of outflow salinity and high-resolution simulations of the exchange flow volume flux. To illustrate the practicality of use of the EBM in an Earth system model, the EBM is implemented for all coastal grid cells with river runoff in the Community Earth System Model (CESM). Compared to the standard version of CESM, which treats runoff as an augmentation to precipitation, the EBM increases sea surface salinity and reduces stratification near river mouths. The EBM also leads to significant regional and remote changes in CESM ocean surface salinities.

Natural and human impacts on centennial sediment accumulation patterns on the Umpqua River margin, Oregon

Abstract Quantifying patterns of sediment accumulation rates (SARs) over the past ~ 125 years on continental margins provides insight into diverse processes spanning the land–ocean boundary. In particular, temporal changes in the export of fluvial sediment can lead to changes in ocean margin SARs, whereas spatial patterns of accumulation reflect the net effect of wave-driven resuspension and current transport averaged over many years. To explore these issues we quantified SARs using 210 Pb geochronology at 73 stations on the shelf and upper slope off the Umpqua River, Oregon. Three types of 210 Pb profiles were observed: a well-defined roughly 10-cm-thick surface mixing layer (SML) underlain by a zone of logarithmic decrease on the distal mid to outer shelf and slope (type 1, n = 45); a deep (20–30-cm thick) SML underlain by a zone of logarithmic decrease in the sandy inner shelf (type 2, n = 8); and a composite profile with two distinct zones of logarithmic decrease below an ~ 10-cm-thick SML in a mid-shelf depocenter adjacent to the river mouth (type 3, n = 21). Type 3 profiles imply a 2–4-fold increase in the SAR that occurred, on average, in 1967 ± 13 y. Such an increase in SARs is consistent with the history of industrial logging in the Umpqua basin, which peaked in the two decades after World War II and coincided with a wet phase of the Pacific Decadal Oscillation (1944–1978) when average and peak river flows were elevated. Comparison of the Umpqua River shelf depocenter with the well-studied Columbia and Eel River systems reveals some important similarities and differences between these ‘marine-dispersal dominated’ systems, whereby forcing (e.g., river–ocean coherence) on the Umpqua is comparable to that of the Eel, whereas the spatial pattern is more similar to that of the Columbia. These seemingly paradoxical results can be reconciled by considering the relative significant wave height during periods of elevated sediment delivery.

Correlation of Tephra Marker Beds in Latest Pleistocene and Holocene Fill of the Submerged Lower Columbia River Valley, Washington and Oregon, U.S.A

Abstract Peterson, C.; Minor, R.; Gates, E.B.; Vanderburgh, S., and Carlisle, K., 2012. Correlation of tephra marker beds in latest Pleistocene and Holocene fill of the submerged Lower Columbia River Valley, Washington and Oregon, U.S.A. Seven tephra layers in the submerged fill of the Lower Columbia River Valley (LCRV) are radiocarbon dated (0.5–13 ka in age). They are correlated to reported eruption sources in the adjacent Cascade Range volcanic arc. The tephras provide a stratigraphic framework for the transgressive depositional filling of the tidal Columbia River system (∼200 km in length). The thickest and most continuous tephra marker beds in the LCRV are from the Mount Mazama set-O eruption at 7.7 ka. The eruption source for this tephra (1–10 m thick) is confirmed by glass geochemistry and radiocarbon dating of borehole samples from depths of −13 to −38 m (North American Vertical Datum of 1988 [NAVD88]). The Mazama set-O tephra conveniently divides the depositional sequences of the Early Holocene f...

Plunging inflows and the summer photic zone in reservoirs

Kinbasket and Revelstoke Reservoirs are part of the Columbia River system in eastern British Columbia, Canada. Hydroelectricity from these large reservoirs represents about 30% of the province's generation capacity. Of interest to water use planning is the effect of reservoir operation on pelagic productivity. We address one small part of this question, namely, the supply of nutrients carried by inflows that plunge below the photic zone during the summer. Using an idealized water balance for the photic zone, three cases are considered: (1) a shallow outlet, (2) a deep outlet, and (3) no outflow. For a shallow outlet, all inflow that plunges below the photic zone is upwelled into the photic zone on its way to the outlet. For a deep outlet, inflow that plunges below the photic zone will short circuit or pass directly to the outlet unless entrainment generates upwelling of the inflow into the photic zone. For a reservoir with no outflow, such as a reservoir that is filling, inflow that plunges below the photic zone remains at depth unless either entrainment or a bathymetric effect generates flow into the photic zone; nutrients are then released when the reservoir is drawn down, often in winter.

Population-specific escapement of Columbia River fall Chinook salmon: Tradeoffs among estimation techniques

In the multi-stock Columbia River system, managers estimate fall Chinook salmon, Oncorhynchus tshawytscha (Walbaum), escapements using various combinations of spawning ground surveys, harvest data and fish counts at dams and hatcheries. Our objectives were to improve upon the traditional methods, and to evaluate trade-offs among methods. Using data from radio-tagged (n=4421) and PIT-tagged (n=1950) adult salmon over eight years, we applied a mark-recapture method to estimate population-specific escapements, both aggregating data within year and stratifying them by week. Mark-recapture estimates differed between estimation techniques and from estimates generated using traditional methods. Stratifying data by week measured escapement estimate uncertainty more reasonably than aggregating data within year. Radiotelemetry provided better spatial resolution among populations for tributary spawners whereas PIT tags provided low-cost, easily replicated estimates using an existing detection system. Mark-recapture techniques had several advantages over current practices: quantifying uncertainty, transparent methods and reduced sensitivity to survey biases.

Human and climate impacts on Columbia River hydrology and salmonids

AbstractThe Columbia River had until recently the world's largest Chinook salmon runs. Restoration of the system's severely decimated runs requires understanding changes in the hydrologic variables (e.g. flow and sediment transport) important to salmonids. We describe here methods to distinguish the human and climate‐induced contributions to Columbia River hydrologic processes relevant to the crucial seaward spring migration of juveniles through the tidal river and estuary. Flow regulation has caused most of the decrease in peak flow and sediment transport; it has contributed to changes in spring freshet timing. Climate change has reduced peak and average flows and sediment transport, changing spring‐freshet timing by several weeks. Irrigation diversion has reduced the annual average flow as much as climate change. A better understanding of historical changes in hydrologic processes entailed in this paper tells us how management and climate have changed the Columbia River system over time. The separation of the climate and anthropogenic influences used here may assist in policy analyses and strategies aimed at restoration of the Columbia River endangered salmonids, and in management of other systems. Copyright © 2010 John Wiley & Sons, Ltd.

First Documented Spawning of White Sturgeon in the Lower Willamette River, Oregon

White sturgeon (Acipenser transmontanus) in the lower Columbia River system were previously known to spawn in only one area; within the 12 km of mainstem Columbia River downstream of Bonneville Dam. Our work provides the first documentation of white sturgeon spawning in the Willamette River, Oregon, a major tributary of the lower Columbia River. We used artificial substrates to sample for white sturgeon eggs downstream of Willamette Falls from 18 May to 20 May 2009 and collected a total of 22 fertilized eggs. Embryonic developmental stages ranged from pigmentation change (Stage 2) to early epithelial (Stage 8), corresponding to fertilization times of approximately 5–17 hours prior to collection. We estimated that spawning occurred between 2100 hours on 19 May and 0900 hours on 20 May 2009, and that a minimum of three independent spawning events took place during that time. Results suggest that the area immediately downstream of Willamette Falls may be important white sturgeon spawning habitat, and that the lower Willamette River is likely an additional source of production for the white sturgeon population in the lower Columbia River system.

Autonomous Acoustic Receiver Deployment and Mooring Techniques for Use in Large Rivers and Estuaries

AbstractAcoustic telemetry receivers are used across a range of aquatic habitats to study a diversity of aquatic species. The Juvenile Salmon Acoustic Telemetry System autonomous acoustic receiver system was deployed and moored in the Columbia River and its estuary. A high receiver loss rate during 2005 in the Columbia River estuary, an area with high water velocities and unstable substrates, prompted improvements to the receiver mooring system, and in 2006 the mooring system was redesigned. This change included elimination of surface buoys, a cable bridle, and an anchor tagline (for anchor recovery). The new mooring configuration, consisting of an acoustic receiver, acoustic release, and mooring line sections that were anchored to the riverbed, improved receiver recovery rates and crew safety. Additionally, a reward program was implemented to provide an incentive for people to return found receivers. The mooring design presented here performed well compared with previous acoustic receiver mooring methods used in the Columbia River system and should be useful for similar applications in large rivers and estuaries with high water velocities and shifting substrates.

A semi‐Lagrangian water temperature model for advection‐dominated river systems

This paper describes a one‐dimensional stream temperature model that is computationally efficient and highly scalable in both time and space. The model is developed within the framework of state space structure. The time‐dependent equations for the conservation of thermal energy in a flowing stream or river are solved using a mixed Eulerian‐Lagrangian, or semi‐Lagrangian, numerical scheme. Solutions are obtained by tracking individual water parcels along their flow characteristics and storing the simulated results at discrete points on a fixed grid. Computational efficiency and accuracy of the numerical scheme are demonstrated by comparison of model estimates with observations of stream temperatures from rivers in the Pacific Northwest as well as with results from a closed‐form solution of the energy equation. A preliminary analysis of the impact of climate changes on stream temperature in the Columbia River system illustrates the strengths of the semi‐Lagrangian method for addressing water quality issues of regional, national, and, ultimately, global scale. Further development of the semi‐Lagrangian method has the potential to improve the ability of water quality planners to perform uncertainty analysis, risk analysis, and forecasting for large, complex river systems.

The Portland Basin: A (big) river runs through it

Metropolitan Portland, Oregon, USA, lies within a small Neogene to Holocene basin in the forearc of the Cascadia subduction system. Although the basin owes its existence and structural development to its convergent-margin tectonic setting, the stratigraphic architecture of basin-fill deposits chiefly reflects its phys iographic position along the lower reaches of the continental-scale Columbia River system. As a result of this globally unique setting, the basin preserves a complex record of aggradation and incision in response to distant as well as local tectonic, volcanic, and climatic events. Voluminous flood basalts, continental and locally derived sediment and volcanic debris, and catastrophic flood deposits all accumulated in an area influenced by contemporaneous tectonic deformation and variations in regional and local base level.

Propagation of the myxozoan parasite Myxobolus cerebralis by different geographic and genetic populations of Tubifex tubifex: An Oregon perspective

Tubifex tubifex are obligate invertebrate hosts in the life cycle of Myxobolus cerebralis, the myxozoan parasite that causes whirling disease in salmonid fishes. This exotic parasite is established to varying degrees across Oregon’s Columbia River system (Pacific Northwest, USA) and characteristics of local T. tubifex populations likely play a role in the pattern of disease occurrence. To better understand these patterns, we collected T. tubifex from three Oregon river basins (Willamette, Deschutes, and Grande Ronde), determined their genotype (mitochondrial 16S rDNA lineage and RAPD genotype) and exposed 10 different populations to M. cerebralis in the laboratory . Four mt lineages were identified: I, III, V and VI. Lineage III was found in all river basins but dominated both central and eastern sites. The RAPD assay further divided these lineages into geographic sub-populations; no RAPD genotype was common to all basins. There was a significant difference in prevalence of infection and level of parasite production among the populations we exposed to M. cerebralis that was attributed to genotypic composition. Only lineage III worms released actinospores and only populations dominated by this lineage amplified the parasite. These populations had the lowest survival, however, the lineage dominant before exposure remained dominant despite the high prevalence of infection. The distribution and infection dynamics of susceptible T. tubifex throughout Oregon may contribute to the differences in M. cerebralis occurrence; our studies further support the influence of oligochaete genotypes on the manifestation of whirling disease in salmonid populations.

Bayesian analysis of mark‐recapture data with travel time‐dependent survival probabilities

AbstractThe authors extend the classical Cormack‐Jolly‐Seber mark‐recapture model to account for both temporal and spatial movement through a series of markers (e.g., dams). Survival rates are modeled as a function of (possibly) unobserved travel times. Because of the complex nature of the likelihood, they use a Bayesian approach based on the complete data likelihood, and integrate the posterior through Markov chain Monte Carlo methods. They test the model through simulations and apply it also to actual salmon data arising from the Columbia river system. The methodology was developed for use by the Pacific Ocean Shelf Tracking (POST) project.

Non‐direct homing behaviours by adult Chinook salmon in a large, multi‐stock river system

Two non‐direct homing behaviours, overshoot of natal tributaries and temporary non‐natal tributary use, were evaluated for 5150 radio‐tagged spring–summer Chinook salmon Oncorhynchus tshawytscha from 40 populations in the large Columbia River system. Over 7 years, 2–44% (mean = 15%) of individuals within each group temporarily entered presumed non‐natal tributaries. In addition, many Chinook salmon from lower river tributaries initially travelled 3 to >250 km upstream in the main‐stem river beyond confluences with presumed natal tributaries before returning to the natal sites (‘overshoot’). Both overshoot and temporary tributary use behaviours declined exponentially with increasing distance from the natal tributary. Non‐direct homing also increased later in the season as water temperatures rose and was associated with hatchery origin in some cases. The behaviours may reflect a mix of active searching for olfactory cues from natal sites, behavioural thermoregulation and orientation challenges in a large‐river migration corridor transformed by dams and reservoirs. While anadromous salmonid homing is generally accurate and precise, these results indicate that route finding can be non‐direct, potentially increasing energetic costs and harvest risks during migration.