Start Submission

Reading: Legacies of stream channel modification revealed using General Land Office surveys, with imp...

Download

A- A+
Alt. Display

Research Article

Legacies of stream channel modification revealed using General Land Office surveys, with implications for water temperature and aquatic life

Authors:

Seth M. White ,

Department of Fishery Science, Columbia River Inter-Tribal Fish Commission, Portland, Oregon, US
X close

Casey Justice,

Department of Fishery Science, Columbia River Inter-Tribal Fish Commission, Portland, Oregon, US
X close

Denise A. Kelsey,

Department of Fishery Science, Columbia River Inter-Tribal Fish Commission, Portland, Oregon, US
X close

Dale A. McCullough,

Department of Fishery Science, Columbia River Inter-Tribal Fish Commission, Portland, Oregon, US
X close

Tyanna Smith

School of the Environment, Portland State University, Portland, Oregon, US
X close

Abstract

Land use legacies can have a discernible influence in present-day watersheds and should be accounted for when designing conservation strategies for riverine aquatic life. We describe the environmental history of three watersheds within the Grande Ronde subbasin of the Columbia River using General Land Office survey field notes from the 19th century. In the two watersheds severely impacted by Euro-American land use, stream channel widths—a metric representing habitat simplification—increased from an average historical width of 16.8 m to an average present width of 20.8 m in large streams; 4.3 m to 5.5 m in small, confined or partly confined streams; and 3.5 m to 6.5 m in small, laterally unconfined steams. Conversely, we did not detect significant change in stream widths in an adjacent, wilderness stream with minimal human impact. Using a mechanistic water temperature model and restoration scenarios based on the historical condition, we predicted that stream restoration in the impacted watersheds could notably decrease average water temperatures—especially when channel narrowing is coupled with riparian restoration—up to a 6.6°C reduction in the upper Grande Ronde River and 3.0°C in Catherine Creek. These reductions in water temperature translated to substantial changes in the percentage of stream network habitable to salmon and steelhead migration (from 29% in the present condition to 79% in the fully restored scenario) and to core juvenile rearing (from 13% in the present condition to 36% in the fully restored scenario). We conclude that land use legacies leave an important footprint on the present landscape and are critical for understanding historic habitat-forming processes as a necessary first step towards restoration.

Knowledge Domain: Ecology
How to Cite: White, S.M., Justice, C., Kelsey, D.A., McCullough, D.A. and Smith, T., 2017. Legacies of stream channel modification revealed using General Land Office surveys, with implications for water temperature and aquatic life. Elem Sci Anth, 5, p.3. DOI: http://doi.org/10.1525/elementa.192
466
Views
145
Downloads
Altmetric
 Published on 21 Feb 2017
 Accepted on 22 Nov 2016            Submitted on 17 Jun 2016
Domain Editor-in-Chief: Don Zak; University of Michigan, Michigan, United States
Guest Editor: Marcelo Ardon-Sayao; North Carolina State University, North Carolina, United States


Introduction

A major challenge of the Anthropocene—the period in which human activity is the dominant influence on climate and the environment—is to solve the interrelated problems leading to irreversible damage to planetary life support systems. These intertwined problems include human population growth, overconsumption, land use, climate change, and subsequent extinctions to biodiversity and elimination of ecosystem services (Foley et al., 2005; Barnosky et al., 2016). A common approach to forecasting the effects of human activity on the environment is through modeling scenarios of land use change and climate conditions, revealing various possible futures that can be embraced, avoided, or mitigated (Moss et al., 2010; Jantz et al., 2015; Isaak et al., 2016). In order to accomplish this, a comprehensive understanding of past human activities is needed, especially when past actions propagate a legacy extending to the present (Foster et al., 2003). The integrity of rivers and streams is especially vulnerable to human activities because hydrology and water temperature are strongly influenced by climatic effects (Dittmer, 2013) and the landscapes over which they flow (Hynes, 1975; Fausch et al., 2002; Allan, 2004). Streams and rivers provide important ecosystem services including clean and abundant water supply that are difficult to valuate but nonetheless essential (Arthington et al., 2010). Degradation of riverine ecosystems represents an important loss in terms of aquatic biodiversity (Dudgeon et al., 2006) and to people that depend upon rivers for food and other cultural values (Close et al., 2002).

Modifications to river ecosystems in Europe, U.S., and other locations across the globe have been well documented. The European subcontinent has experienced land use change—specifically urbanization—since 700 B.C. (Antrop, 2004); these patterns have been manifested in several ways, but primarily as landscape fragmentation (Jaeger et al., 2011) and river channelization (Jurajda, 1995) with negative consequences to aquatic life (Muxika et al., 2007; Haidvogl et al., 2015). Tropical and subtropical rivers have experienced intense localized impacts from human settlement over several decades (Webb, 1992), with more ominous threats looming from increased interest in building large dams (Bergkamp et al., 2000; White et al., 2012). Perhaps because of a longer period of human settlement as compared to the U.S. West coupled with extensive historical documentation, land use legacies in the U.S. have been broadly implicated in degradation of streams and rivers east of the Continental Divide (Harding et al., 1998; Scott, 2006; Wenger et al., 2008; Gardiner et al., 2009; Walter and Merritss, 2008; Maloney and Weller, 2011; Einheuser et al., 2013). Land use impacts in watersheds of the U.S. West are also pervasive (McIntosh et al., 1994; Robbins and Wolf, 1994; Wallin et al., 1994), yet legacies have been described for streams and rivers in the region less frequently (see however Sedell and Froggatt, 1984; McIntosh et al., 2000; White and Rahel, 2008).

Land use has been implicated as a leading cause of river channel simplification with subsequent consequences to aquatic life. Channel widening is one common form of simplification (Figure 1) and occurs through various human-related causes including increased flooding after removal of native hillslope vegetation, riparian vegetation, and large woody debris (Knox, 1977; Faustini and Jones, 2003); eroding banks and sediment deposition in the stream channel (Beschta, 1983; Simon and Rinaldi, 2006; Allan, 2004); and is linked to various upstream land use such as timber harvest, road networks, and livestock grazing (Dose and Roper, 1994; Ralph et al., 1994; Knapp et al., 1998; Kondolf et al., 2002). The capacity for river width adjustment is also strongly reliant on geomorphic setting of the river channel. Broader alluvial channels with less confinement by hillslopes are more sensitive to changes in channel morphology (Montgomery and Buffington, 1997; Thorne, 1998; Faustini et al., 2009). In addition to decreasing habitat complexity important for rearing salmonids and other aquatic biota, channel widening increases surface water area and the capacity for solar radiation to reach the stream, thereby increasing water temperature (Poole and Berman, 2001) which can have negative physiological and behavioral impacts on organisms adapted to cold temperatures (Margesin and Schinner, 1999; Dell et al., 2014).

Figure 1 

Illustration of channel complexity in a wilderness versus modified stream. Stream channel complexity contrasted in two watersheds including the Minam River—a designated wilderness area representing historical conditions—as seen from (A) aerial photography and (B) ground photography in 2013, and a simplified and unnaturally widened stream channel in the heavily impacted upper Grande Ronde River as seen from (C) aerial photography and (D) ground photography in 2015. Aerial images courtesy of Google Earth; ground images courtesy of the Columbia Habitat Monitoring Program. DOI: https://doi.org/10.1525/elementa.192.f1

In this study, we describe the implications of watershed history to changes in average stream conditions from the 19th century to present. To accomplish this, we made use of General Land Office (GLO) surveys within three watersheds of the Columbia River basin, each encompassing distinct spring Chinook Salmon Oncorhynchus tshawytscha and steelhead O. mykiss populations (NMFS, 2016). The GLO surveys are a unique resource of historical information prior to major Euro-American impacts, and have been used by historical ecologists to describe past conditions and set restoration targets (Egan, 2005). The GLO surveys were intended to provide information on the quality of conditions for rangeland, agriculture, and forestry for prospective land claims. However several studies have used GLO information for other purposes, primarily for describing historical riparian vegetation communities (Johnson, 1994; Galat et al., 1998; Hulse et al., 2000; McAllister, 2008; Dilts et al., 2012). Less frequently, GLO surveys have been employed to represent in-channel characteristics of streams and rivers, such as wood recruitment (Sedell and Froggatt, 1984; Collins et al., 2002) or river channel morphology (McDowell, 2000; Collins et al., 2003; Hereford and Betancourt, 2009).

To define a historical baseline for fish habitat, we described changes to stream channel widths since the late 1800s, with expectations that the magnitude of change would be greater in areas with more intense ranching, logging, agriculture, and other forms of land use. Therefore, the specific objectives of this study were to: (1) evaluate overall patterns of stream channel widening since the 19th century, with reference to the geomorphic context where modification has been most severe and where restoration efforts may have the greatest physical capacity for improvement; and (2) simulate the effects of stream channel widening and riparian vegetation on water temperature using a mechanistic stream temperature model, with implications for aquatic life.

Methods

Study area

This study was conducted in tributaries of the Grande Ronde River originating in the Blue Mountains and Wallowa Mountains of Northeast Oregon, United States, and flowing 334 km to its confluence with the Snake River and eventually the Columbia River. Focal watersheds include two tributaries heavily impacted by anthropogenic land use—the upper Grande Ronde River above the town of La Grande (draining 1,896 km2) and Catherine Creek (1,051 km2)—and one least-impacted watershed in the Eagle Cap Wilderness—the Minam River (618 km²) (Figure 2A). Topography is characterized by rugged mountains in the headwaters (2,269 m) and a broad, low gradient valley at the confluence of the upper Grande Ronde River and Catherine Creek (820 m). The climate is characterized by cold, moist winters and warm, dry summers with average temperatures near La Grande averaging –0.42°C in January and 21°C in July. Average annual precipitation ranges from 36 cm in the valleys to 152 cm in the mountains, with most of the precipitation in the mountains falling as winter snow. Due to the lower elevation of the Blue Mountains relative to the Wallowas, snowmelt generally occurs earlier in its tributaries, often resulting in very low stream flows during summer (Kelly and White, 2016).

Figure 2 

Study area, stream classification, and historical changes to channel widths in three focal watersheds. Location of study watershed in northeast Oregon including (A) major salmon-bearing tributaries and the stream classification described in the methods and (B) values of channel change estimates where historical General Land Office surveys intersected with contemporary Aquatic Inventory Program surveys. Focal watersheds include the upper Grande Ronde River, Catherine Creek, and Minam River. The upper Grande Ronde River and Catherine Creek have significantly modified stream conditions from over a century of intensive land use. The Minam River is in the Eagle Cap Wilderness area and most approximates historical reference conditions. DOI: https://doi.org/10.1525/elementa.192.f2

Watershed conditions in the upper Grande Ronde River and Catherine Creek, like many watersheds in the U.S. West, have experienced degraded ecological health caused by cumulative influence of past human activity since the 19th century (McIntosh et al., 1994; Wissmar et al., 1994). In the upper Grande Ronde River and Catherine Creek, intensive land use impacts occurred beginning in 1812 (Table 1) with beaver (Castor canadensis) trapping and proceeded through the late 1980s with activities as diverse as draining marshland and diking river sections for agriculture in the Grande Ronde Valley; logging hillslopes and riparian areas with associated splash damming and, later, road building; cattle and sheep grazing on public and private land; dredge mining; and damming (Gildemeister, 1998). These land use practices have been implicated in decades-long trends in simplification of stream habitat through loss of deep pools across the Columbia River basin, whereas streams in wilderness or roadless areas retained their more complex nature relative to streams in managed watersheds (McIntosh et al., 2000).

Table 1

Land use history of the upper Grande Ronde River and Catherine Creek from 1812 to 1989a. DOI: https://doi.org/10.1525/elementa.192.t1

Year Event

1812 Fur trader Robert Stuart observes beaver as common; War of 1812 intensifies beaver trapping
1850 Donation Land Claim Law enacted, encouraging settlement by emigrants
1855 Treaty between U.S. and upper Columbia River Indian Tribes exchanging ceded lands for reservations and reserving rights to traditional hunting, fishing, and gathering
1861–1862 First land claim in Grande Ronde Valley; first sawmill and salmon-blocking dam built on Grande Ronde River
1865–1869 Water-powered flourmill established on Catherine Creek; railroad route laid out across Grande Ronde Valley
1870 Construction begins on State Ditch and Catherine Creek ditch draining lakes and swamplands in the Grande Ronde Valley
1890 Grande Ronde Lumber Company acquires timberland and begins constructing splash dams on Grande Ronde River and tributaries
1890s–1900 Development of railroad network in Grande Ronde tributaries; estimated 50 sawmills in watershed, annual timber export estimated at 32.5 million board feet
1934 Taylor Grazing Act leads to decline of livestock grazing on public lands
1939 Mine dredging begins in Grande Ronde River
1946 Establishment of Union Co. Soil and Water Conservation District leads to substantial land leveling, ditching, and stream channeling projects
1984–1985 Log jams blasted on Catherine Creek to alleviate flooding; Army Corps of Engineers clears willow and cottonwood from riparian zones
1989 Recognition that peak flows shifting as much as 30 days earlier based on 1904–1989 record, partly attributed to land use in watershed

aAbbreviated from Gildemeister, 1998.

We selected the Minam River as a reference watershed based on its less intensive land use history and its near proximity to the other focal watersheds. The Minam River flows through the Wallowa-Whitman National Forest and Eagle Cap Wilderness Area at an average elevation of 1251 m in our study extent. The Eagle Cap was established as a primitive area in 1930, designated as wilderness in 1940, and registered in the National Wilderness Preservation System in 1964. In 1988, the Minam River was registered as a Wild and Scenic River from its headwaters at Minam Lake, 62.8 river km downstream to Cougar Creek. The protected status of the Minam River provides a stark contrast to the intense present and historical agricultural, grazing, and logging use in the upper Grande Ronde and Catherine Creek basins, making the Minam River a good candidate stream to represent reference conditions.

However, intrinsic physical conditions unrelated to land use are somewhat different among the basins. While most of the Grande Ronde River and lower sections of Catherine Creek flow out of Miocene and younger volcanic and sedimentary rocks (with some higher elevations in the Grande Ronde River flowing from older, pre-Cenozoic sedimentary and volcanic rock), the Minam River and upper sections of Catherine Creek flow out of Oligocene and lower Miocene volcanic and sedimentary rock (Walker, 1990). Both the Minam River and upper sections of Catherine Creek have notable, U-shaped valleys carved by glaciers, and wet climate due to the orographic effect. The geomorphology of reaches in the Minam River is most similar to that of Catherine Creek, but with a smaller proportion of tributaries with low gradient and with a larger proportion of valleys constrained by hillslope walls. The upper Grande Ronde River is unique compared to the other two watersheds having a higher proportion of tributaries with low gradient and unconstrained valleys. Study sites in the Minam River were selected within a range of intrinsic watershed characteristics (elevation, upstream watershed area, cumulative precipitation, valley width index, etc.) that most corresponded with the impacted watersheds.

Spring Chinook Salmon populations in the upper Grande Ronde River and Catherine Creek were listed as threatened under the Endangered Species Act in 1992. Population declines over the past century were due in part to severely degraded habitat conditions resulting from the aforementioned anthropogenic disturbances. Specifically, stream temperature, streamflow, fine sediment, habitat diversity, and quantity of key habitats such as large pools in these basins have been identified as key limiting factors for recovery of fish populations (Nowak and Kuchenbecker, 2004).

Assessing historical stream conditions

Historical estimates of channel width were based on GLO surveys conducted in the mid- to late-1800s within the area of study in the Grande Ronde subbasin. The GLO established the Public Land Survey System of townships and ranges in 1812 (White, 1983), with subsequent modifications to methods in Oregon that were applied in other U.S. States (Principle Clerk of Surveys, 1855). Public lands were apportioned into townships 9.7 km (6 mi) on a side, and townships were further divided into 36 sections, each 1.6 km (1 mi) on a side. The survey involved GLO surveyors walking the section lines for each section in a township. In addition to recording the general character of vegetation, soil, and rangeland conditions, surveyors recorded the location and bank-to-bank channel width (active channel width) of any streams or rivers crossed. We accessed GLO field notes on the U.S. Bureau of Land Management’s Official Land Records Site (BLM, 2016) and translated the handwritten notes into spatial data in a geographic information system (GIS) (ESRI, 2011). Estimates of channel width were converted from chains and links (1 chain = 100 links) to meters (1 link = 0.20 m).

Assessing contemporary stream conditions

Contemporary estimates of channel width were based on Oregon Department of Fish and Wildlife’s Aquatic Inventories Project (AIP) (Moore et al., 2008). The AIP survey is a rapid assessment of common fish habitat characteristics collected in a spatially continuous fashion across the stream network. Two AIP surveyors walked smaller streams or canoed unwadeable sections and recorded the characteristics and location of channel units (i.e., pools, riffles, and glides) with a hand-held global positioning system (GPS) with accuracy 5–7 m. Data from the 1990s were used as the baseline for present conditions, except where surveys were conducted outside the low flow period (ordinal date 200–300; Kelly and White, 2016). When surveys did not match those criteria we used surveys from years 2000 or 2010 that fell within the low flow period. In our study, channel width was used as a proxy of width:depth ratio—a metric strongly tied to integrity of stream channels (e.g., Beschta and Platts 1986; Myers and Swanson 1996) and commonly used in fish-habitat models (Fausch et al., 1988)—because historical estimates of water depths were not available. Constraining the use of all survey data to only the low flow period presumably provided consistency in discharge over the years that would allow change in width to be a valid surrogate for change in width:depth ratio.

Whereas AIP surveyors recorded wetted width at every channel unit, active channel width was recorded only at every 10th channel unit and at tributary junctions. We therefore developed a linear relationship between co-occurrences of measured wetted and active channel width using data from the Columbia Habitat Monitoring Program (CHaMP, 2016) collected during low-flow periods in 2011–2015 within the study basins. CHaMP is an intensive stream habitat survey focusing on detailed, reach-scale stream channel geomorphology; a much larger sample size was available from this program (n = 131) than from AIP (n = 23) and encompassed a broader range of stream sizes. This provided confidence in extrapolating spatially extensive channel widths comparable to active width as recorded in historical GLO surveys:

WP=2.12 + 1.22  ×  WW,

where WP is the predicted present active channel width and WW is the present wetted width (m) measured by AIP crews (n = 131, R2 = 0.89, p < 0.001) (Figure 3). In estimating channel change over time (described below), locations where the historical active channel width was equal to or smaller than the model intercept (i.e., ≤ 2.12 m) were removed from the analysis to avoid potential upward bias in estimated stream widths for small channels in the present time period.

Figure 3 

Relationship between active versus wetted stream width. Linear relationship between active (WP) versus wetted (WW) stream channel widths from contemporary Aquatic Inventory Program stream surveys (WA = 2.12 + 1.22 × WW, n = 131, R2 = 0.89, p < 0.001). Shaded area represents 95% confidence band. DOI: https://doi.org/10.1525/elementa.192.f3

Estimating historical changes to stream channel width

We calculated the percentage change in channel width from the historical to present periods based on GLO and AIP estimates of active channel width:

Δ  W=WPWHWH  ×  100,

where ΔW is percentage change in channel width, WP is present channel width, and WH is historical channel width (m). We then evaluated the magnitude of change since the historical period according to watershed identity and a geomorphic valley setting classification. Watershed identity was defined by the extent of the following spring Chinook Salmon populations: Catherine Creek Chinook (CCC), upper Grande Ronde Chinook (UGC), and Minam River Chinook (MRC). The classification system consisted of dividing the stream network into small and large streams using an 8-m bankfull width criterion based on the work of Beechie and Imaki (2014). Next, the stream network was further divided into three different valley types based on valley confinement (laterally unconfined, partly confined, and confined) following the methodology described in the River Styles Framework (Brierley and Fryirs, 2005). Based on exploratory analysis of fish habitat conditions among stream types, we simplified the classification into three classes for our study: large streams (LS), small/partly confined and confined streams (SC), and small/laterally unconfined streams (SU) (Figure 2A). The effect of watershed identity on magnitude of channel change was tested using one-way analysis of variance (ANOVA). One-way ANOVA was also used to test the effect of valley setting on magnitude of channel change, but only for sites in the impacted watersheds (CCC and UGC) because all locations in the Minam River where estimates of channel change existed were in large stream types. Model assumptions were confirmed as valid by visually evaluating residuals versus fits, normal Q-Q plots, scale-location plots, residuals versus leverage, and histograms of residuals. Tukey’s HSD test was used for post-hoc evaluations of individual group differences.

Water temperature modeling and assessing importance to aquatic life

To evaluate the contribution towards cooling stream temperatures from restoring channel widths to their historical state in combination with revegetation of riparian zones, we employed a mechanistic water temperature model called Heat Source (Boyd and Kasper, 2003). The model integrates stream channel geometry, hydrology, climatic conditions, and riparian vegetation cover and height to simulate stream temperature and effective shade at 100-m intervals throughout the stream network. The model was calibrated for a 10-week period between 10 July and 20 September 2010. This period was chosen to best represent present conditions for summer base-flow conditions when water temperatures are typically highest and salmonids are consequently at risk.

Model parameters from present conditions were used as a baseline for evaluating restoration scenarios including (1) restoring stream channels to their historical widths, (2) restoring riparian vegetation to its potential natural state, and (3) a combination of channel width and riparian restoration. Channel width scenarios were developed by assigning the average value of channel width change by geomorphic stream classification across the modeled extent. Potential natural vegetation scenarios were developed by estimating the potential height and canopy cover of trees and shrubs in the riparian zone under natural historical conditions using a detailed map of present vegetation and potential natural vegetation (PNV) along the entire extent of the Chinook-bearing portion of the upper Grande Ronde River and Catherine Creek watersheds (Wells et al., 2015). Additional details of water temperature model development, riparian scenario development, and application of additional restoration and climate change scenarios are discussed in McCullough et al. (2016) and Justice et al. (2017).

We summarized the water temperature predictions for each model scenario by calculating the maximum 7-day running average of the daily maximum temperature (MWMT) in degrees Celsius. This metric has been demonstrated as an important metric for various salmon life stages in Pacific Northwest streams (EPA, 2003). The median MWMT was compared for both upper Grande Ronde River and Catherine Creek (including major tributaries) for present conditions and all restoration scenarios. To visualize spatial patterns in how restoration scenarios affected water temperature, we mapped MWMT across the upper Grande Ronde River and Catherine Creek for present conditions and under the channel width restoration scenario. Lastly, we reported the percentage of stream length in the river networks having water temperatures below critical Pacific salmon (Oncorhynchus spp.) and steelhead (O. mykiss) water temperature thresholds for both watersheds combined for present conditions and restoration scenarios Thresholds corresponded to 16°C for adult holding and core juvenile rearing and 20°C for fish migration (EPA, 2003).

Results

Changes to channel widths since historical period

A total of 193 intersections between GLO and AIP surveys throughout the three watersheds allowed for the estimation of historical and present stream widths (Figure 2B). There was a significant effect of watershed identity on the percentage change in channel width as determined by one-way ANOVA (F[2, 190] = 10.71, p < 0.001) (Figure 4A). Post-hoc Tukey comparisons revealed that percentage change in channel width in Catherine Creek (CCC) was greater than in the Minam River (MRC) by 63.8 (p-adj < 0.001, 95% CI = 30.5–97.1). Percentage change in channel width in the upper Grande Ronde River (UGC) was greater than in the Minam River by 37.9 (p-adj = 0.02, 95% CI = 4.4–71.4). Percentage change in channel width in the upper Grande Ronde River was less than Catherine Creek by 25.9 (p-adj = 0.03, 95% CI = –50.1–1.7). Because the Minam River was more homogenous regarding stream classification type (i.e., having a majority of reaches in the large stream category, some reaches in the small/partly confined and confined category, and no reaches in the small/laterally unconfined category), we additionally tested for the effect of watershed identity on percentage change in channel width using only sites in the large stream category to ensure that the differences were not a function of disparity in stream classification types. Again, we noted a significant effect of watershed identity on the percentage change in channel width as determined by one-way ANOVA (F[2, 167] = 8.2, p < 0.001). A similar pattern of differences in percentage change in channel widths among impacted versus wilderness watersheds was confirmed using post-hoc Tukey comparisons (p-adj[CCC-MRC] > 0.001, p-adj[UGC-MRC] = 0.01), whereas channel change between the two impacted streams did not significantly differ (p-adj[UGC-CCC] = 0.30). These results provided justification for evaluating percentage change in channel width as a function of stream classification for upper Grande Ronde and Catherine Creek watersheds combined, and separately from the Minam River.

Figure 4 

Percentage change from historical stream widths by watershed and stream type. Boxplots of percentage change in channel widths from historical (1880s) to present (1990s and later) by (A) Chinook population (all locations combined, n = 193) and (B) Stream type (impacted watersheds only, n = 174). Channel width change in Catherine Creek (CCC) and upper Grande Ronde River (UGC) was significantly different than in the Minam River (MRC), but not from each other. All stream classification groups—large streams (LS); small, partly confined or confined streams (SC); and small, laterally unconfined streams (SU)—were significantly different from one another. DOI: https://doi.org/10.1525/elementa.192.f4

A total of 164 intersections between GLO and AIP surveys throughout Catherine Creek and the upper Grande Ronde River allowed for the estimation of historical and present stream widths in these two watersheds where intensive land use had occurred (Table 2). A majority of the locations available for comparison were in large streams (LS) (n = 141) with fewer in the small, partly confined and confined streams (SC) (n = 15) and small, laterally unconfined streams (SU) (n = 8). However this roughly matched the proportion of stream types by river kilometer within the Heat Source model extent: LS (66.7%), SC (20.7%), and SU (12.6%). Stream channel widths increased from an average historical width of 16.8 m to an average present width of 20.8 m in large streams; 4.3 m to 5.5 m in small, confined or partly confined streams; and 3.5 m to 6.5 m in small, laterally unconfined steams. There was a significant effect of stream classification type on the percentage change in channel width as determined by one-way ANOVA (F[2, 161] = 4.5, p = 0.01) (Figure 4B; Table 2). Post-hoc Tukey comparisons revealed that percentage change in channel width in small/laterally unconfined sites (SU) was significantly greater than in large stream sites by 69.8 (p-adj = 0.02, 95% CI = 11.2–128.5). Percentage change in channel width in small/laterally unconfined sites was significantly greater than in small/partly confined and confined sites by 84.5 (p-adj = 0.01, 95% CI = 13.9–155.2). Percentage change in channel width in small/partly confined and confined sites (SC) was not significantly different than in large stream sites (LS) (p-adj = 0.71).

Table 2

Mean stream channel widths, increases from historical dimensions, and percentage change from historical condition in impacted watersheds. DOI: https://doi.org/10.1525/elementa.192.t2

Stream type Sample size (n) Mean historical channel width (m) ±SE Mean present channel width (m) ±SE Mean increase (m) ±SE Mean change (%) ±SE

Large streams (LS) 141 16.8 ± 0.8 20.8 ± 0.8 4.0 ± 0.6 45.9 ± 5.7
Small/partly confined & confined (SC) 15 4.3 ± 0.2 5.5 ± 0.4 1.2 ± 0.5 31.2 ± 13.0
Small/laterally unconfined (SU) 8 3.5 ± 0.5 6.5 ± 0.7 3.0 ± 1.1 115.8 ± 37.6

Water temperature simulations and implications for aquatic life

Mean percentage change in channel width by stream classification (Table 2) was applied to channel widths in the Heat Source water temperature model to yield estimates for water temperature under the restoration scenarios (Table 3). Water temperatures under present conditions were substantially higher in the upper Grande Ronde River (median MWMT = 24.4°C) compared with Catherine Creek (median MWMT = 18.3°C). The predicted change in median MWMT relative to the present condition for the restored channel width scenario was –2.2°C in the upper Grande Ronde River, compared with –0.6°C for Catherine Creek. The predicted change in median MWMT relative to the present condition for the restored potential natural vegetation (PNV) was –5.5°C in the upper Grande Ronde River, compared with –2.4°C for Catherine Creek. The combined PNV and channel width restoration was estimated to change median water temperatures by –6.6°C in the upper Grande Ronde River compared with –3.0°C in Catherine Creek.

Table 3

Simulated water temperature under current conditions and restoration scenarios. DOI: https://doi.org/10.1525/elementa.192.t3

Mainstem Median MWMTa (°C)
Current Widthb PNVc PNV_Widthd

Upper Grande Ronde River (UGC) 24.4 22.2 18.9 17.8
Catherine Creek (CCC) 18.3 17.7 15.9 15.3

aMaximum weekly maximum water temperature.

bRestored potential channel width scenario.

cRestored potential natural vegetation scenario.

dCombination of vegetation and channel width restoration scenarios.

The longitudinal profile of present (i.e., 2010) MWMT in the mainstem upper Grande Ronde River (Figure 5A) showed rapidly increasing water temperatures from about 12°C near its headwaters to approximately 25°C just downstream of the Sheep Creek confluence. At that point, the river enters a canyon with considerable topographic shade and higher tree cover, and consequently, river temperatures declined moderately to approximately 22°C near the mouth of Fly Creek. From that point downstream to its confluence with Catherine Creek, the river temperature increased gradually with some relatively minor cooling effects at tributary junctions to a maximum of approximately 29°C. In Catherine Creek, water temperature starts out at approximately 16°C at the confluence of North and South Fork Catherine Creek (with upstream MWMT < 12°C), with gradual warming to approximately 22°C in the Grande Ronde Valley.

Figure 5 

Simulated water temperature for present conditions and restoration scenarios. Map of simulated maximum weekly maximum water temperatures (MWMT) in the upper Grande Ronde River and Catherine Creek under (A) present conditions (Current) and (B) the restored channel width scenario (Width). See Figure 2 for stream names. DOI: https://doi.org/10.1525/elementa.192.f5

In the channel width restoration scenario, we noted patterns of decreasing MWMT throughout the two impacted watersheds after mapping model results (Figure 5B). For example, in the upper Grande Ronde River, narrowing channel widths in headwater tributaries (upper mainstem Grande Ronde and Sheep Creek) cooled water with beneficial effects extending downstream to the mouth of Fly Creek. From that point downstream, MWMT increased gradually to temperatures matching those of the present condition. In Catherine Creek, restored channel widths led to cooler water temperatures extending from the confluence of North and South Fork Catherine Creek downstream, through the river’s canyon and alluvial reaches, and into the Grande Ronde Valley. However MWMT values at the model’s downstream extent in the valley were similar between present condition and the restored channel width scenario.

For both upper Grande Ronde River and Catherine Creek watersheds combined, the percentage stream length with MWMT below critical salmon and steelhead thresholds increased with restoration scenarios as compared to present conditions (Figure 6). Percentage stream length having MWMT below migration and holding/rearing thresholds was approximately 29% and 13% for present conditions (respectively), 39% and 17% for restoration of channel widths, 67% and 32% for restoration of potential natural vegetation, and 79% and 36% for the combination of channel width and vegetation restoration.

Figure 6 

Percentage stream length below biological water temperature thresholds for model scenarios. Estimated percentage of stream length below critical salmon and steelhead thresholds for maximum weekly maximum water temperatures (MWMT) (EPA, 2003) in the upper Grande Ronde River and Catherine Creek watersheds combined. Model scenarios represent current conditions (Current), restored channel width (Width), restored potential natural vegetation (PNV), and the combination of vegetation and channel width restoration (Width_PNV). DOI: https://doi.org/10.1525/elementa.192.f6

Discussion

Stream channel widening since the 19th century

Using estimates of channel width from 1880s General Land Office surveys as compared to present day estimates, we noted significant widening of stream channels in the watersheds impacted by anthropogenic land use (upper Grande Ronde River and Catherine Creek) as opposed to a wilderness stream with less human disturbance (Minam River), where channel widening was absent or minimal. Channel widening has been described as one response of unstable, alluvial stream channels to watershed land use via eroding stream banks and channel incision (Thorne, 1998; Simon and Rinaldi, 2006). The upper Grande Ronde River and Catherine Creek—alongside other watersheds in the American West—have been subjected to intensive land use activities including removal of beaver, logging and associated road building, railroad and road encroachments, diking, ditching, dredging, sheep ranching, and cattle grazing (Robbins and Wolf, 1994; Gildemeister, 1998; McIntosh et al., 2000). This cocktail of land use has likely contributed to stream channel simplification in numerous ways, including channel widening as noted in our study, but also loss of large pools (McIntosh et al., 2000) that are important refugia for spawning and rearing fish (Torgersen et al., 2006). Channel widening in the wilderness stream (Minam River) was absent or minimal, lending evidence to the premise that multiple forms of land use are linked to stream channel simplification.

In the upper Grande Ronde River and Catherine Creek, an important driver of stream channel widening is likely increased flooding after widespread modifications to riparian and hillslope vegetation. Historical accounts of settlers and newspaper articles describe drastic increases of the magnitude and timing of flooding corresponding to increasing land use, especially forest harvest and associated activities (Gildemeister, 1998) (Table 1). In New England watersheds, changes to channel morphology from increased flooding were attributed to a large-scale natural disturbance (hurricane) combined with reduced interception of precipitation by vegetation, evaporation from leaf surfaces, and transpiration of moisture from the forest canopy (Foster et al., 2004). Subsequent increases in flood magnitude can lead to channel widening by scouring the stream channel with sediment, bedload material, and coarse woody debris (Ralph et al., 1994). Splash damming—the practice of staging logs in a dammed pond with a sudden release for transporting logs to downstream mill sites—also occurred historically in the upper Grande Ronde River (Gildemeister, 1998) and has been implicated in simplifying stream channels elsewhere in the Pacific Northwest (Miller, 2010). Stream channels in the upper Grande Ronde River may have also been intentionally widened to promote better drainage and decrease the impacts of localized flooding (Gildemeister, 1998).

Cattle grazing is widespread in the upper Grande Ronde River and Catherine Creek watersheds, an activity implicated in increasing channel widths through bank erosion and deposition of fine sediments over the stream channel in other watersheds (Dose and Roper, 1994; Ralph et al., 1994; Knapp et al., 1998; Kondolf et al., 2002). Unnaturally high populations of native ungulates, specifically elk (Cervus canadensis), may also have similar effects due to lack of native predators (i.e., wolves, Canis lupus) that would otherwise reduce ungulate populations (Beschta and Ripple, 2006) or cause behavioral shifts in foraging, reducing grazing impacts to streams (Ripple and Beschta, 2004).

Understanding the geomorphic context of channel widening may help inform the kinds of restoration activities that will halt or reverse these trends. For example, researchers in the nearby Middle Fork John Day River concluded that the potential for adjusting the channel planform to desired conditions was limited by natural planform, but that adding large woody debris may overcome this impediment (McDowell, 2000). In our study, channel widening as a proportion of the original width was more predominant in smaller channels that were laterally unconfined by hillslopes. This was not surprising, given that laterally unconfined stream channels are zones of deposition and net sediment accumulation (Brierley and Fryirs, 2005); any anthropomorphic changes in the upstream watershed would register downstream in these unconfined reaches (Allan, 2004). However, small channels that were either partly confined or confined by hillslopes exhibited less overall proportional change in channel widths, indicating that width adjustments were indeed more prevalent in laterally unconfined channels. Streams in alluvial valleys—corresponding to the laterally unconfined stream classification in our study—have a higher potential for channel geometry responses to changes in sediment supply and discharge (Montgomery and Buffington, 1997) and are known to respond to increases in discharge and sediment by becoming wider and shallower (Ralph et al., 1994). Channel widening in alluvial streams occurs through various pathways including bank erosion without incision; retreat of outer banks when toe scouring exceeds the rate of advance of the opposite bank; or in braided channels, by erosion from flows deflected around advancing bars (Thorne, 1998).

A significant body of literature indicates that channel widening via various pathways can be arrested or reversed through restoration activities. Restoration or protection of riparian vegetation enhances root strength, contributing to increased stability of streambanks (Simon and Collison, 2002), especially in loosely-packed alluvial deposits that are characteristic of unconfined channels (Montgomery and Buffington, 1997). Rooting of riparian vegetation on bars and streambanks can help channels narrow by trapping sediment and removing soils that would otherwise remain suspended in the water column (Boon and Raven, 2012). Riparian vegetation and large woody debris additionally provide sources of roughness that can reduce erosion during high streamflows (Gregory et al., 1991; Ralph et al., 1994). Reduction in the intensity of land use that reduces hillslope or riparian vegetation (e.g., forest harvest and cattle grazing) decreases erosion and sedimentation (Allan, 2004) and increases water storage capacity of soils, thereby reducing the potential for unnaturally high peak flows (Poff et al., 1997). Fortunately, these and other restoration activities meant to address channel simplification are already being initiated in the upper Grande Ronde River and Catherine Creek basins (Booth et al., 2016). However, it remains to be seen whether the extent and intensity of stream channel restoration, along with changes to upstream land use, are sufficient to meet biological objectives for the ESA-listed fish populations (Simon, 2016).

Implications to water temperature and aquatic life

Restoration scenarios yielded positive results in terms of reducing water temperature. Riparian vegetation restoration yielded the greatest benefit alone of any one single approach, but the combination of riparian and channel width restoration yielded the greatest benefit overall. If intensive and widespread restoration actions were successfully implemented in the upper Grande Ronde River, temperatures could be reduced below the present temperature by as much as 6.6°C in the upper Grande Ronde River and 3.0°C in Catherine Creek. The pattern of increasing water temperature with increasing channel width has been documented in other modeling analyses investigating land use impacts (LeBlanc et al., 1997; ODEQ, 2009; Butcher et al., 2010). Simulations of channel narrowing yielded a small cooling benefit as compared to restored vegetation in another study in the upper Grande Ronde River (ODEQ, 2000); however, justification for baseline channel widths in that study was unstated. In the nearby John Day River, a modeling analysis demonstrated that a 30% reduction in channel width yielded an approximately equivalent reduction in water temperature compared with vegetation restoration in the upper 100 km of the river, and a substantially greater reduction compared with vegetation restoration in the lower 325 km (Butcher et al., 2010); justification for baseline channel widths for that study came from unreferenced “basin literature” indicating historical channel widths were 5–50% narrower. Restoration scenarios in our study did not account for climate change and hyporheic exchange, factors that are also important determinants of stream temperature (Poole and Berman, 2001) but were beyond the scope of our analysis. A fruitful restoration program with the goal of reducing water temperature would address riparian shade, channel morphology (width, average depth, sinuosity, bed roughness, etc.), groundwater-hyporheic-surface water connectivity, and upstream sources of sedimentation. The potential benefits of these actions should additionally be evaluated in the context of climate change (Wu et al., 2012).

In our study, all permutations of riparian and channel width restoration scenarios increased the amount of stream length below critical biological thresholds (EPA, 2003). The greatest single benefit was from riparian vegetation; however restoring channel width alone provided an increase from 29% to 39% of stream length with water temperatures less than the 20°C migration threshold and from 13% to 17% of stream length less than the 16°C adult holding and juvenile rearing thresholds for salmon and steelhead. Not surprisingly, the greatest benefit was from the combined riparian and channel width scenarios: the percentage of usable habitat increased to 79% and 36% of stream length for migration and holding/rearing, respectively. This finding underlines the importance of applying a combination of restoration efforts in a comprehensive program accounting for short-term and long-term benefits from habitat actions (Roni et al., 2002, 2008), especially considering that riparian restoration can take decades to yield improvements (Hasselquist et al., 2015).

We used a simplified approach to assessing the biological importance of water temperature by merely tallying the stream kilometers that could be utilized by fish according to published thresholds (EPA, 2003). Determining the true benefit to fish from reduction in water temperature would include accounting for other local factors linked to riparian restoration that also improve fish habitat, such as large wood delivery or pool development (Fausch et al., 1988), food availability and growth (Weber et al., 2014), fish carrying capacity (Lobón-Cerviá, 2008), fish behavioral responses (White et al., 2014), physiological responses of fish (Feldhaus et al., 2010), and other spatial factors such as the juxtaposition of habitats meeting requirements for multiple life stages (Jackson et al., 2001; White and Rahel, 2008). Fish response to temperature regimes may also be highly dependent on spatially discontinuous coldwater refuges (Ebersole et al., 2003) that are not captured in coarse-grained stream temperature models. These factors would ideally be used in a life cycle modeling framework accounting for survival bottlenecks in multiple life history stages (e.g., Scheuerell et al., 2006). However, a broad assessment of water temperature across the stream network was helpful for documenting the existing and potential template over which more nuanced factors affecting fish distribution can play out. The concept of thresholds implies that above certain thermal limits, physiological performance is severely limited enough to preempt colonization or success in warm-water sections of the stream network (EPA, 2003; McCullough, 2010). Therefore, the simulated increases in available stream length from restoration (Figure 6) should be considered the potential gain in available habitat; whether or not fish occupy or thrive there will depend on additional factors.

Using General Land Office surveys for constructing watershed histories

In combination with other historical data sources, General Land Office surveys can provide valuable information on watershed conditions prior to major Euro-American settlement and land use impacts (McAllister, 2008). These pre-impact descriptions can provide insights into how watershed conditions have changed over time, what are the major drivers, and how fish and other aquatic life may still be responding to land use legacies (Harding et al., 1998). Whereas numerous studies have employed GLO surveys to provide information on historical vegetation, to our knowledge, fewer studies have used GLO surveys to describe historical modifications to stream channel morphology across an entire stream network. Fitzpatrick and Knox (2000) used GLO surveys to describe historical versus present-day channel widths and sediment conditions in an assessment of flooding effects on channel geomorphology in North Fish Creek, Wisconsin. Beckham (1995) documented historical stream channel widths along the mainstem upper Grande Ronde River, Oregon, using GLO surveys, but the data were used for descriptive purposes and were not compared to contemporary estimates. Graf (1981) used a time series of GLO plat maps (along with other data sources) to assess potential zones of hazardous channel migration in the Gila River, Arizona. McDowell (2000) used GLO records from 1881 to describe original channel location (along with riparian vegetation) in an assessment of anthropogenic versus natural drivers for channel change in the Middle Fork John Day River, Oregon. Collins et al. (2003) used GLO surveys and other sources of information to reconstruct the historical riverine landscapes—including previous channel locations and riparian vegetation—of Puget lowlands, Washington. Each of these studies provided invaluable insights on reconstructing river channel and floodplain characteristics to help inform historical patterns of riverine habitat development, a necessary first step towards restoration (Ebersole and Liss, 1997).

Conclusions

Historical ecology involves using multiple information sources as lines of evidence, often in a manner inconsistent with the original purpose of the data collection; however, if caution is taken we can begin to discern helpful insights regarding the character of the changing landscape (Fuller et al., 2004). Setting target conditions for restoration typically involves inferring conditions from nearby, undisturbed reference areas and using statistical models to extrapolate the expected, unimpacted conditions from within the existing range of anthropogenic disturbance in a watershed (Pollock et al., 2012). However, when historical information is available, it can provide a more realistic estimate of baseline conditions and range of natural variability (Motzkin and Foster, 2004). Reconstructing historical conditions does not necessarily imply a target for restoration. However, understanding the manner in which watersheds have been altered can improve our understanding of how and why conditions have changed (Pedroli et al., 2002). Furthermore, an understanding of the past may help us to avoid future mistakes. Using GLO surveys in comparison to contemporary stream surveys, we found that watersheds impacted by human land use had widened stream channels, especially in smaller streams in less-confined valleys. Restoration activities meant to return channels to their historical dimensions may have a greater physical impact in small, unconfined channels because those channels were most impaired by land use and have a greater capacity for geomorphic change. However, we emphasize that ability to affect change in a geomorphic context is only one of several criteria for prioritizing restoration actions; other important factors for planning restoration include biological benefits associated with actions and the social, economic, and overall land use objectives that set the context for restoration (Beechie et al., 2008). Restoration scenarios that included both restoration of riparian vegetation and stream channel narrowing projected reduced water temperature and increased length of the stream network habitable by salmonids. GLO surveys as a source of historical information can be valuable in describing broadscale watershed histories. Perhaps the most important benefit of reconstructing watershed histories in the present study—and in general—is the ability to use historical baselines to shed light on the legacy of processes constraining the abundance, productivity, and spatial distribution of aquatic life.

Data accessibility statement

The following publicly available datasets were used for analyses:

Acknowledgements

Thanks to Ariane Peralta and Marcelo Ardón of East Carolina University and other participants of the “Ghosts of land-use past” session at the 2014 Joint Aquatic Sciences Meeting in Portland, Oregon, that was the impetus for this paper. We additionally thank David Graves at Columbia River Inter-Tribal Fish Commission for assistance with geographical information systems during this project, Jazzmine Allen of Roosevelt High School for assistance collecting General Land Office data during her internship with CRITFC, and Edwin Sedell at Oregon Department of Fish & Wildlife for providing Aquatic Inventories Project stream survey data. Timothy Beechie and Hiroo Imaki at National Oceanic at Atmospheric Administration provided a draft protocol for transferring GLO records into GIS. The multi-agency Grande Ronde and Catherine Creek Atlas restoration planning group in La Grande, Oregon, provided invaluable feedback on a preliminary workshop presentation of this analysis.

Funding Information

This work was funded by Bonneville Power Administration (BPA) Project No. 2009-004-00. Contributions from TS were supported by the Ronald E. McNair Scholars Program, Portland State University.

Competing Interests

The authors have no competing interests to declare.

Author Contributions

  • Contributed to conception and design: SMW, CJ, DAK, DAM
  • Contributed to acquisition of data: SMW, CJ, DAM, DAK, TS
  • Contributed to analysis and interpretation of data: SMW, CJ, DAK
  • Drafted and/or revised the article: SMW, CJ, DAK, DAM
  • Approved the submitted version for publication: SMW, CJ, DAK, DAM, TS

References

  1. Allan JD 2004. Landscapes and riverscapes: The influence of land use on stream ecosystems. Annu Rev Ecol Syst 35: 257–284, DOI: http://dx.doi.org/10.1146/annurev.ecolsys.35.120202.110122 

  2. Antrop M 2004. Landscape change and the urbanization process in Europe. Landscape Urban Plan 67(1–4): 9–26, DOI: http://dx.doi.org/10.1016/S0169-2046(03)00026-4 

  3. Arthington AH, Naiman RJ, McClain ME and Nilsson C 2010. Preserving the biodiversity and ecological services of rivers: New challenges and research opportunities. Freshwater Biol 55(1): 1–16, DOI: http://dx.doi.org/10.1111/j.1365-2427.2009.02340.x 

  4. Barnosky AD, Ehrlich PR and Hadly EA 2016. Avoiding collapse: Grand challenges for science and society to solve by 2050. Elem Sci Anth 4: 94.DOI: http://dx.doi.org/10.12952/journal.elementa.000094 

  5. Beckham SD 1995. Grande Ronde River, Oregon: River widths, vegetative environment, and conditions shaping its condition, Imbler vicinity to headwaters In: Walla Walla, WA: Eastside Ecosystem Management Project.  

  6. Beechie T and Imaki H 2014. Predicting natural channel patterns based on landscape and geomorphic controls in the Columbia River Basin, USA. Water Resour Res 50(1): 39–57, DOI: http://dx.doi.org/10.1002/2013WR013629 

  7. Beechie T, Pess G, Roni P and Giannico G 2008. Setting river restoration priorities: A review of approaches and a general protocol for identifying and prioritizing actions. N Am J Fish Man 28: 891–905, DOI: http://dx.doi.org/10.1577/M06-174.1 

  8. Bergkamp G, McCartney M, Dugan P, McNeely J and Acreman M 2000. Dams, ecosystem functions and environmental restoration In: WCD Thematic Review. Capetown, South Africa: World Commission on Dams (WCD).  

  9. Beschta RL 1983. Long-term changes in channel widths of the Kowai River, Torlesse Range, New Zealand. J Hydro: New Zealand 22(2): 112–122.  

  10. Beschta RL and Platts WS 1986. Significance and function of morphological features of small streams. Water Resour Bull 22(3): 369–379, DOI: http://dx.doi.org/10.1111/j.1752-1688.1986.tb01891.x 

  11. Beschta RL and Ripple WJ 2006. River channel dynamics following extirpation of wolves in northwestern Yellowstone National Park, USA. Earth Surf Proc Land 31(12): 1525–1539, DOI: http://dx.doi.org/10.1002/esp.1362 

  12. BLM (US Bureau of Land Management) 2016. The Official Land Records Site [dataset].  https://www.glorecords.blm.gov/default.aspx 

  13. Boon, PJ and Raven, PJ eds. 2012. River Conservation and Management. 1st ed. Oxford, UK: Wiley-Blackwell.  

  14. Booth D, Scholz J, Beechie T and Ralph S 2016. Integrating limiting-factors analysis with process-based restoration to improve recovery of endangered salmonids in the Pacific Northwest, USA. Water 8(5): 174.DOI: http://dx.doi.org/10.3390/w8050174 

  15. Boyd M and Kasper B 2003. Analytical methods for dynamic open channel heat and mass transfer: Methodology for the Heat Source model version 7.0 In: Portland, OR: Watershed Sciences, Inc..  

  16. Brierley GJ and Fryirs KA 2005. Geomorphology and River Management: Applications of the River Styles Framework. Oxford: Blackwell.  

  17. Butcher D, Crown J, Brannan K, Kishida K and Hubler S 2010. John Day River basin total maximum daily load (TMDL) and water quality management plan (WQMP) In: Portland, OR: Oregon Department of Environmental Quality. DEQ 10-WQ-025. 

  18. CHaMP (Columbia Habitat Monitoring Program) 2016. Scientific protocol for salmonid habitat surveys within the Columbia Habitat Monitoring Program In: Prepared by CHaMP for the Bonneville Power Administration.  

  19. Close DA, Fitzpatrick MS and Li HW 2002. The ecological and cultural importance of a species at risk of extinction, Pacific lamprey. Fisheries 27(7): 19–25, DOI: http://dx.doi.org/10.1577/1548-8446(2002)027<0019:TEACIO>2.0.CO;2 

  20. Collins BD, Montgomery DR and Haas AD 2002. Historical changes in the distribution and functions of large wood in Puget lowland rivers. Can J Fish Aquat Sci 59(1): 66–76, DOI: http://dx.doi.org/10.1139/f01-199 

  21. Collins BD, Montgomery DR and Sheikh AJ 2003. Reconstructing the historical riverine landscape of the Puget Lowland In: Montgomery, DR, Bolton, S, Booth, DB and Wall, L eds.  Restoration of Puget Sound Rivers. Seattle, WA: University of Washington Press.  

  22. Dell AI, Pawar S and Savage VM 2014. Temperature dependence of trophic interactions are driven by asymmetry of species responses and foraging strategy. J Anim Ecol 83(1): 70–84, DOI: http://dx.doi.org/10.1111/1365-2656.12081 

  23. Dilts TE, Weisberg PJ, Yang J, Olson TJ, Turner PL and Condon LA 2012. Using historical General Land Office survey notes to quantify the effects of irrigated agriculture on land cover change in an arid lands watershed. Ann Assoc Am Geogr 102(3): 531–548, DOI: http://dx.doi.org/10.1080/00045608.2011.641479 

  24. Dittmer K 2013. Changing streamflow on Columbia basin tribal lands—climate change and salmon. Climate Change 120(3): 627–641, DOI: http://dx.doi.org/10.1007/s10584-013-0745-0 

  25. Dose JJ and Roper BB 1994. Long-term changes in low-flow channel widths within the South Umpqua Watershed, Oregon. J Am Water Resour As 30(6): 993–1000, DOI: http://dx.doi.org/10.1111/j.1752-1688.1994.tb03347.x 

  26. Dudgeon D Arthington AH Gessner MO Kawabata ZI Knowler DJ et al. 2006. Freshwater biodiversity: Importance, threats, status and conservation challenges. Biol Rev 81(2): 163–182, DOI: http://dx.doi.org/10.1017/S1464793105006950 

  27. Ebersole JL and Liss WJ 1997. Restoration of stream habitats in the western United States: Restoration as reexpression of habitat capacity. Environ Manage 21(1): 1–14, DOI: http://dx.doi.org/10.1007/s002679900001 

  28. Ebersole JL, Liss WJ and Frissell CA 2003. Thermal heterogeneity, steam channel morphology, and salmonid abundance in northeastern Oregon streams. Can J Fish Aquat Sci 60(10): 1266–1280, DOI: http://dx.doi.org/10.1139/f03-107 

  29. Egan D 2005. The Historical Ecology Handbook: A Restorationist’s Guide to Reference Ecosystems. Washington, DC: Island Press.  

  30. Einheuser MD, Nejadhashemi AP, Wang L, Sowa SP and Woznicki SA 2013. Linking biological integrity and watershed models to assess the impacts of historical land use and climate changes on stream health. Environ Manage 51(6): 1147–1163, DOI: http://dx.doi.org/10.1007/s00267-013-0043-7 

  31. EPA 2003. Environmental Protection Agency Region 10 guidance for Pacific Northwest state and tribal temperature water quality standards In: Seattle, WA: Environmental Protection Agency, Region 10 Office of Water. EPA 910-B-03-002. 

  32. ESRI 2011. ArcGIS Desktop: Release 11 In: Redlands, CA: Environmental Systems Research Institute.  

  33. Fausch KD, Hawkes CL and Parsons MG 1988. Models that predict standing crop of stream fish from habitat variables: 1950–85 In: General Technical Report PNW-GTR-213. Portland, OR: U.S. Department of Agriculture, Forest Service, Pacific Northwest Research Station.  

  34. Fausch KD, Torgersen CE, Baxter CV and Li HW 2002. Landscapes to riverscapes: Bridging the gap between research and conservation of stream fishes. BioScience 52(6): 483–498, DOI: http://dx.doi.org/10.1641/0006-3568(2002)052[0483:LTRBTG]2.0.CO;2 

  35. Faustini JM and Jones JA 2003. Influence of large woody debris on channel morphology and dynamics in steep, boulder-rich mountain streams, Western Cascades, Oregon. Geomorphology 51(1–3): 187.DOI: http://dx.doi.org/10.1016/S0169-555X(02)00336-7 

  36. Faustini JM, Kaufmann PR and Herlihy AT 2009. Downstream variation in bankfull width of wadeable streams across the conterminous United States. Geomorphology 108(3/4): 292–311, DOI: http://dx.doi.org/10.1016/j.geomorph.2009.02.005 

  37. Feldhaus JW, Heppell SA, Li HW and Mesa MG 2010. A physiological approach to quantifying thermal habitat quality for redband rainbow trout (Oncorhynchus mykiss gairdneri) in the South Fork John Day River, Oregon. Environ Biol Fish 87(4): 277–290, DOI: http://dx.doi.org/10.1007/s10641-010-9580-6 

  38. Fitzpatrick FA and Knox JC 2000. Spatial and temporal sensitivity of hydrogeomorphic response and recovery to deforestation, agriculture, and floods. Phys Geogr 21(2): 89–108, DOI: http://dx.doi.org/10.1080/02723646.2000.10642701 

  39. Foley JA DeFries R Asner GP Barford C Bonan G et al. 2005. Global consequences of land use. Science 309(5734): 570–574, DOI: http://dx.doi.org/10.1126/science.1111772 

  40. Foster D Cooper-Ellis S Plotkin AB Carlton G Bowden R et al. 2004. Simulating a catostrophic hurricane In: Foster, D and Aber, J eds.  Forests in Time: The Environmental Consequences of 1,000 Years of Change in New England. New Haven & London: Yale University Press.  

  41. Foster D Swanson F Aber J Burke I Brokaw N et al. 2003. The importance of land-use legacies to ecology and conservation. BioScience 53(1): 77–88, DOI: http://dx.doi.org/10.1641/0006-3568(2003)053[0077:TIOLUL]2.0.CO;2 

  42. Fuller J, Foster D, Motzkin G, McLachlan J and Barry S 2004. Broadscale forest response to land use and climate change In: Foster, D and Aber, J eds.  Forests in Time: The Environmental Consequences of 1,000 Years of Change in New England. New Haven & London: Yale University Press.  

  43. Galat DL Fredrickson LH Humburg DD Bataille KJ Bodie JR et al. 1998. Flooding to restore connectivity of regulated, large-river wetlands. BioScience 48(9): 721–733, DOI: http://dx.doi.org/10.2307/1313335 

  44. Gardiner EP Sutherland AB Bixby RJ Scott MC Meyer JL et al. 2009. Linking stream and landscape trajectories in the southern Appalachians. Environ Monit Assess 156(1): 17–36, DOI: http://dx.doi.org/10.1007/s10661-008-0460-x 

  45. Gildemeister J 1998. Watershed history: Middle & upper Grande Ronde River subbasins, northeast Oregon In: La Grande, OR: Oregon Department of Environmental Quality, U.S. Environmental Protection Agency, and the Confederated Tribes of the Umatilla Indian Reservation.  

  46. Graf WL 1981. Channel instability in a braided, sand bed river. Water Resour Res 17(4): 1087–1094, DOI: http://dx.doi.org/10.1029/WR017i004p01087 

  47. Gregory SV, Swanson FJ, McKee SW and Cummins KW 1991. An ecosystem perspective of riparian zones: Focus on links between land and water. Bioscience 41(8): 540–551, DOI: http://dx.doi.org/10.2307/1311607 

  48. Haidvogl G, Hoffmann R, Pont D, Jungwirth M and Winiwarter V 2015. Historical ecology of riverine fish in Europe. Aquat Sci 77(3): 315–324, DOI: http://dx.doi.org/10.1007/s00027-015-0400-0 

  49. Harding JS, Benfield EF, Bolstad PV, Helfman GS and Jones EDB 1998. Stream biodiversity: The ghost of land use past. P Natl Acad Sci 95(25): 14843–14847, DOI: http://dx.doi.org/10.1073/pnas.95.25.14843 

  50. Hasselquist EM Nilsson C Hjältén J Jørgensen D Lind L et al. 2015. Time for recovery of riparian plants in restored northern Swedish streams: A chronosequence study. Ecol Appl 25(5): 1373–1389, DOI: http://dx.doi.org/10.1890/14-1102.1 

  51. Hereford R and Betancourt JL 2009. Historic geomorphology of the San Pedro River: Archival and physical evidence In: Stromberg, J and Tellman, B eds.  Ecology and Conservation of Desert Riparian Ecosystems: The San Pedro River Example. Tucson, AZ: University of Arizona Press.  

  52. Hulse D, Eilers J, Freemark K, Hummon C and White D 2000. Planning alternative future landscapes in Oregon: Evaluating effects on water quality and biodiversity. Landscape J 19(1–2): 1–19, DOI: http://dx.doi.org/10.3368/lj.19.1-2.1 

  53. Hynes HBN 1975. The stream and its valley. Proceedings of the International Association of Theoretical and Applied Limnology. 19(1): 1–15.  

  54. Isaak DJ Young MK Luce CH Hostetler SW Wenger SJ et al. 2016. Slow climate velocities of mountain streams portend their role as refugia for cold-water biodiversity. P Natl Acad Sci 113(16): 4374–4379, DOI: http://dx.doi.org/10.1073/pnas.1522429113 

  55. Jackson DA, Peres-Neto PR and Olden JD 2001. What controls who is where in freshwater fish communities: The roles of biotic, abiotic, and spatial factors. Can J Fish Aquat Sci 58(1): 157–170, DOI: http://dx.doi.org/10.1139/f00-239 

  56. Jaeger JAG, Soukup T, Madriñán LF, Schwick C and Kienast F 2011. Landscape fragmentation in Europe In: Copenhagen, Denmark: European Environment Agency. EEA Report No 2/2011. 

  57. Jantz SM Barker B Brooks TM Chini LP Huang Q et al. 2015. Future habitat loss and extinctions driven by land-use change in biodiversity hotspots under four scenarios of climate-change mitigation. Conserv Biol 29(4): 1122–1131, DOI: http://dx.doi.org/10.1111/cobi.12549 

  58. Johnson WC 1994. Woodland expansions in the Platte River, Nebraska: Patterns and causes. Ecol Monogr 64(1): 45–84, DOI: http://dx.doi.org/10.2307/2937055 

  59. Jurajda P 1995. Effect of channelization and regulation on fish recruitment in a floodplain river. Regul River 10(2–4): 207–15, DOI: http://dx.doi.org/10.1002/rrr.3450100215 

  60. Justice C, White SM, McCullough DM, Graves DS and Blanchard MR 2017. Can stream and riparian restoration offset climate change impacts to salmon populations?. J Environ Manage 188: 212–217, DOI: http://dx.doi.org/10.1016/j.jenvman.2016.12.005 

  61. Kelly VJ and White SM 2016. A method for characterizing late-season low-flow regime in the upper Grand Ronde River Basin, Oregon. USGS Numbered Series 2016–5041 In: Reston, VA: U.S. Geological Survey. http://pubs.er.usgs.gov/publication/sir20165041 

  62. Knapp RA, Vredenburg VT and Mathews KR 1998. Effects of stream channel morphology on Golden Trout spawning habitat and recruitment. Ecol Appl 8(4): 1104–1117, DOI: http://dx.doi.org/10.1890/1051-0761(1998)008[1104:EOSCMO]2.0.CO;2 

  63. Knox JC 1977. Human impacts on Wisconsin stream channels. Ann Assoc Am Geogr 67(3): 323–342, DOI: http://dx.doi.org/10.1111/j.1467-8306.1977.tb01145.x 

  64. Kondolf GM, Piégay H and Landon N 2002. Channel response to increased and decreased bedload supply from land use change: Contrasts between two catchments. Geomorphology 45(1–2): 35–51, DOI: http://dx.doi.org/10.1016/S0169-555X(01)00188-X 

  65. LeBlanc RT, Brown RD and FitzGibbon JE 1997. Modeling the effects of land use change on the water temperature in unregulated urban streams. J Environ Manage 49(4): 445–469, DOI: http://dx.doi.org/10.1006/jema.1996.0106 

  66. Lobón-Cerviá J 2008. Habitat quality enhances spatial variation in the self-thinning patterns of stream-resident brown trout (Salmo trutta). Can J Fish Aquat Sci 65(9): 2006–2015, DOI: http://dx.doi.org/10.1139/F08-105 

  67. Maloney KO and Weller DE 2011. Anthropogenic disturbance and streams: Land use and land-use change affect stream ecosystems via multiple pathways. Freshwater Biol 56(3): 611–626, DOI: http://dx.doi.org/10.1111/j.1365-2427.2010.02522.x 

  68. Margesin, R and Schinner, F eds. 1999. Cold-Adapted Organisms. Berlin, Heidelberg: Springer Berlin Heidelberg.  

  69. McAllister LS 2008. Reconstructing historical riparian conditions of two river basins in eastern Oregon, USA. Environ Manage 42(3): 412–425, DOI: http://dx.doi.org/10.1007/s00267-008-9127-1 

  70. McCullough DA 2010. Are coldwater fish populations of the United States actually being protected by temperature standards?. Freshwater Rev 3(2): 147–199, DOI: http://dx.doi.org/10.1608/FRJ-3.2.4 

  71. McCullough DA White SM Justice C Blanchard M Lessard R et al. 2016. Assessing the status and trends of spring Chinook habitat in the upper Grande Ronde River and Catherine Creek. Annual Report to Bonneville Power Administration In: Portland, OR: Columbia River Inter-Tribal Fish Commission.  

  72. McDowell PF 2000. Human impacts and river channel adjustment, northeastern Oregon: Implications for restoration. American Water Resources Association, Annual International Summer Specialty Conference. August 27–30, 2000, Portland, Oregon 

  73. McIntosh BA Sedell JR Smith JE Wissmar RC Clarke SE et al. 1994. Management history of eastside ecosystems: Changes in fish habitat over 50 years, 1935 to 1992 In: General Technical Report PNW-GTR-321. Portland, OR: U.S. Department of Agriculture, Forest Service, Pacific Northwest Research Station.  

  74. McIntosh BA, Sedell JR, Thurow RF, Clarke SE and Chandler GL 2000. Historical changes in pool habitats in the Columbia River basin. Ecol Appl 10(5): 1478–1496, DOI: http://dx.doi.org/10.1890/1051-0761(2000)010[1478:HCIPHI]2.0.CO;2 

  75. Miller RR 2010. Is the past present? Historical splash-dam mapping and stream disturbance detection in the Oregon Coastal Province [M.S. thesis]. Corvallis: Oregon State University, Department of Fisheries and Wildlife. http://hdl.handle.net/1957/18998 

  76. Montgomery DR and Buffington JM 1997. Channel-reach morphology in mountain drainage basins. Geol Soc Am Bull 109(5): 596–611, DOI: http://dx.doi.org/10.1130/0016-7606(1997)109<0596:CRMIMD>2.3.CO;2 

  77. Moore K, Jones KK, Dambacher J and Stein C 2008. Methods for stream habitat surveys: Aquatic Inventories Project In: Corvallis, OR: Oregon Department of Fish and Wildlife, Conservation and Recovery Program. http://oregonstate.edu/dept/ODFW/freshwater/inventory/pdffiles/habmethod.pdf 

  78. Moss RH Edmonds JA Hibbard KA Manning MR Rose SK et al. 2010. The next generation of scenarios for climate change research and assessment. Nature 463(7282): 747–756, DOI: http://dx.doi.org/10.1038/nature08823 

  79. Motzkin G and Foster D 2004. Insights for ecology and conservation In: Foster, D and Aber, J eds.  Forests in Time: The Environmental Consequences of 1,000 Years of Change in New England. New Haven & London: Yale University Press.  

  80. Muxika I, Borja Á and Bald J 2007. Using historical data, expert judgment and multivariate analysis in assessing reference conditions and benthic ecological status, according to the European Water Framework Directive. Mar Pollut Bull 55(1–6): 16–29, DOI: http://dx.doi.org/10.1016/j.marpolbul.2006.05.025 

  81. Myers TJ and Swanson S 1996. Long-term aquatic habitat restoration: Mahogany Creek, Nevada, as a case study. Water Resour Bull 32(2): 16–29, DOI: http://dx.doi.org/10.1111/j.1752-1688.1996.tb03448.x 

  82. NMFS 2016. 5-Year review: Summary & evaluation of Snake River Sockeye, Spring-Chinook Salmon, Fall-Run Chinook, and Steelhead In: Portland, OR: National Marine Fisheries Service.  

  83. Nowak MC and Kuchenbecker L 2004. Grande Ronde subbasin plan In: Northwest Power and Conservation Council.  

  84. ODEQ 2000. Upper Grande Ronde River subbasin: Total maximum daily load (TMDL) In: Oregon Department of Environmental Quality.  

  85. ODEQ 2009. Water monitoring and assessment mode of operations manual (MOMs) Version 3.2 In: Hillsboro, Oregon: Oregon Department of Environmental Quality, Laboratory and Environmental Assessment Division. DEQ03-LAB-0036-SOP. 

  86. Pedroli B, de Blust G, van Looy K and van Rooij S 2002. Setting targets in strategies for river restoration. Landscape Ecol 17(1): 5–18, DOI: http://dx.doi.org/10.1023/A:1015221425315 

  87. Poff LN Allan DJ Bain MB Karr JR Prestegaard KL et al. 1997. The natural flow regime: A paradigm for river conservation and restoration. BioScience 41(11): 769–784, DOI: http://dx.doi.org/10.2307/1313099 

  88. Pollock MM, Beechie TJ and Imaki H 2012. Using reference conditions in ecosystem restoration: An example for riparian conifer forests in the Pacific Northwest. Ecosphere 3(11): 1–23, DOI: http://dx.doi.org/10.1890/ES12-00175.1 

  89. Poole GC and Berman CH 2001. An ecological perspective on instream temperature: Natural heat dynamics and mechanisms of human-caused thermal degradation. Environ Manage 27(6): 787–802, DOI: http://dx.doi.org/10.1007/s002670010188 

  90. Principle Clerk of Surveys 1855. Instructions to the Surveyors General of Public Lands of the United States for Those Surveying Districts Established in and Since the Year 1850. Containing Also a Manual of Instructions to Regulate the Field Operations of Deputy Surveyors In: Washington: A.O.P. Nicholson, Public Printer.  

  91. Ralph SC, Poole GC, Conquest LL and Naiman RJ 1994. Stream channel morphology and woody debris in logged and unlogged basins of western Washington. Can J Fish Aquat Sci 51(1): 37–51, DOI: http://dx.doi.org/10.1139/f94-006 

  92. Ripple WJ and Beschta RL 2004. Wolves and the ecology of fear: Can predation risk structure ecosystems?. BioScience 54(8): 755.DOI: http://dx.doi.org/10.1641/0006-3568(2004)054[0755:WATEOF]2.0.CO;2 

  93. Robbins WG and Wolf DW 1994. Landscape and the intermontane Northwest: An environmental history In: General Technical Report PNW-GTR-319. Portland, Oregon: Pacific Northwest Research Station, U.S. Department of Agriculture.  

  94. Roni P Beechie TJ Bilby RE Leonetti FE Pollock MM et al. 2002. A review of stream restoration techniques and a hierarchical strategy for prioritizing restoration in Pacific Northwest watersheds. N Am J Fish Man 22(1): 1–20, DOI: http://dx.doi.org/10.1577/1548-8675(2002)022<0001:AROSRT>2.0.CO;2 

  95. Roni P, Hanson K and Beechie TJ 2008. Global review of the physical and biological effectiveness of stream habitat rehabilitation techniques. N Am J Fish Man 28(3): 856–890, DOI: http://dx.doi.org/10.1577/M06-169.1 

  96. Scheuerell MD Hilborn R Ruckelshaus MH Bartz KK Lagueux KM et al. 2006. The Shiraz model: A tool for incorporating anthropogenic effects and fish–habitat relationships in conservation planning. Can J Fish Aquat Sci 63(7): 1596–1607, DOI: http://dx.doi.org/10.1139/f06-056 

  97. Scott MC 2006. Winners and losers among stream fishes in relation to land use legacies and urban development in the southeastern US. Biol Conserv 127(3): 301–309, DOI: http://dx.doi.org/10.1016/j.biocon.2005.07.020 

  98. Sedell JR and Froggatt JL 1984. Importance of streamside forests to large rivers: The isolation of the Willamette River, Oregon, USA, from its floodplain by snagging and streamside forest removal. Verh Internat Verein Limnol 22: 1828–1834.  

  99. Simon A and Collison AJC 2002. Quantifying the mechanical and hydrologic effects of riparian vegetation on streambank stability. Earth Surf Proc Land 27(5): 527–546, DOI: http://dx.doi.org/10.1002/esp.325 

  100. Simon A and Rinaldi M 2006. Disturbance, stream incision, and channel evolution: The roles of excess transport capacity and boundary materials in controlling channel response. Geomorphology 79(3–4): 361–383, DOI: http://dx.doi.org/10.1016/j.geomorph.2006.06.037 

  101. Simon MH 2016. National Wildlife Federation v. National Marine Fisheries Service 3: 01-cv-00640-SI 1 In: U.S. District Court for the District of Oregon.  

  102. Thorne CR 1998. River width adjustment: Processes and mechanisms. J Hydraul Eng 124(9): 881–902, DOI: http://dx.doi.org/10.1061/(ASCE)0733-9429(1998)124:9(881) 

  103. Torgersen CE, Baxter CV, Li HW and McIntosh BA 2006. Landscape influences on longitudinal patterns of river fishes: Spatially continuous analysis of fish-habitat relationships. Am Fish S S 48: 473–492.  

  104. Walker GW 1990. Cenozoic geology of the Blue Mountains region In: Washington, DC: U.S. Geological Survey. Professional Paper USGS-1437. 

  105. Wallin DO, Swanson FJ and Marks B 1994. Landscape pattern response to changes in pattern generation rules: Land-use legacies in forestry. Ecol Appl 4(3): 569–580, DOI: http://dx.doi.org/10.2307/1941958 

  106. Walter RC and Merritts DJ 2008. Natural streams and the legacy of water-powered mills. Science 319(5861): 299–304, DOI: http://dx.doi.org/10.1126/science.1151716 

  107. Webb JLA 1992. Ecological and economic change along the middle reaches of the Gambia River, 1945–1985. African Affairs 91(365): 543–565.  

  108. Weber N, Bouwes N and Jordan CE 2014. Estimation of salmonid habitat growth potential through measurements of invertebrate food abundance and temperature. Can J Fish Aquat Sci 71(8): 1158–1170, DOI: http://dx.doi.org/10.1139/cjfas-2013-0390 

  109. Wells AF, Crow E and Blaha R 2015. Riparian vegetation mapping in the Grande Ronde watershed, Oregon: Monitoring and validation of spring Chinook habitat recovery and population viability In: Anchorage, AK: ABR, Inc.-Environmental Research & Services. (prepared for the Columbia River Inter-Tribal Fish Commission). 

  110. Wenger SJ, Peterson JT, Freeman MC, Freeman BJ and Homans DD 2008. Stream fish occurrence in response to impervious cover, historic land use, and hydrogeomorphic factors. Can J Fish Aquat Sci 65(7): 1250–1264, DOI: http://dx.doi.org/10.1139/F08-046 

  111. White CA 1983. A history of the rectangular survey system In: U.S. Department of the Interior, Bureau of Land Management.  

  112. White SM, Giannico GR and Li HW 2014. A ‘behaviorscape’ perspective on stream fish ecology and conservation: Linking fish behavior to riverscapes. WIRES: Water 1(4): 385–400, DOI: http://dx.doi.org/10.1002/wat2.1033 

  113. White SM, Ondračková M and Reichard M 2012. Hydrologic connectivity affects fish assemblage structure, diversity, and ecological traits in the unregulated Gambia River, West Africa. Biotropica 44(4): 521–530, DOI: http://dx.doi.org/10.1111/j.1744-7429.2011.00840.x 

  114. White SM and Rahel FJ 2008. Complementation of habitats for Bonneville Cutthroat Trout in watersheds influenced by beavers, livestock, and drought. T Am Fish Soc 137(3): 881–894, DOI: http://dx.doi.org/10.1577/T06-207.1 

  115. Wissmar RC Smith JE McIntosh BA Li HW Reeves GH et al. 1994. Ecological health of river basins in forested regions of Eastern Washington and Oregon In: Portland, OR: U.S. Department of Agriculture, Forest Service, Pacific Northwest Research Station, DOI: http://dx.doi.org/10.2737/PNW-GTR-326 PNW-GTR-326. 

  116. Wu H Kimball JS Elsner MM Mantua N Adler RF et al. 2012. Projected climate change impacts on the hydrology and temperature of Pacific Northwest rivers. Water Resour Res 48(11): 1–23, DOI: http://dx.doi.org/10.1029/2012WR012082 

comments powered by Disqus