High-resolution temperature models are needed to describe and predict the thermal regime of streams/rivers under current and future management. Statistical models (e.g., spatial stream network models – SSN) can cover large extents at fine resolutions. Mechanistic models (e.g., Soil and Water Assessment Tool – SWAT) offer refined temporal resolution and process-based drivers in predictions, but are time consuming to parameterize. Hybrid modeling (e.g., hybrid SSN-SWAT) approaches can leverage benefits of both types to improve temporal robustness of thermal regime model predictions across entire networks. In the Upper Rogue River Basin, OR, we compared performance of a NorWeST August mean temperature (AMT) SSN model with a hybrid SSN-SWAT AMT (Hybrid AMT) model. We extended the hybrid SSN-SWAT approach to the 7-day average daily maximum temperature metric (Hybrid 7DADM). The Hybrid AMT and Hybrid 7DADM models performed well, with similar model diagnostics between the Hybrid AMT and NorWeST AMT models, and with other SSN and mechanistic temperature models (used for temperature total maximum daily loads – TMDLs). The hybrid SSN-SWAT approach was able to incorporate the effect of dams on downstream hydrology and temperature regimes. The best-fitting hybrid SSN models incorporated upstream shade effects spanning distances of 1-5 kilometers or 1-3 hours travel time, highlighting the need to maintain upstream riparian zones to protect downstream thermal regimes.
Forest practices can exacerbate or moderate risks to coldwater habitat (CWH) related to climate variability, wildfire, and forest disease. We applied a high-resolution hybrid mechanistic/statistical stream network (SSN) model to evaluate effects of riparian management practices and future climate in the Upper Rogue, OR basin on flow and temperature. We compared recent (1990 – 2018) and 2040’s climate scenarios for temperature and precipitation, considering both average, minimum, and maximum based on the Palmer Drought Severity Index. Riparian management conditions included historical reference, current vegetation (CV), and vegetation managed via the Oregon Forest Practices (FPA) Act which prescribes minimum forest buffer widths of 100-110 ft for fish-bearing (F) and 35 ft for non fish-bearing (NF) streams. Under the 2040s scenario, flow will be lost by the end of August in 79% of F reaches and 93% of NF reaches in a dry year. Interannual variability of CWH extent will increase during the 2040s. In cool/wet years, core F CWH (< 16 0C) extent will decline by ~ 10 % while habitat > 18 0C will increase under the FPA scenario as compared to the current vegetation regime. For 2040s flowing streams under mean climate, the extent of NF reaches with <15 °C 7DADM will be reduced under the FPA to 2.2% of reaches as compared to current vegetation (4.5 % of reaches). Thus, reduced protections prescribed for NF streams will impact downstream CWH as well as headwater NF habitat.
Many cold water-dependent aquatic organisms are experiencing habitat and population declines from increasing water temperatures. Identifying mechanisms which drive local and regional stream thermal regimes facilitates restoration at ecologically relevant scales. Stream temperatures vary spatially and temporally both within and among river basins. We developed a modeling process to identify statistical relationships between drivers of stream temperature and covariates representing landscape, climate, and management-related processes. The modeling process was tested in three study areas of the Pacific Northwest United States during the growing season (May [start], August [warmest], September [end]). Across all months and study systems, covariates with the highest relative importance represented the physical landscape (elevation [1st], catchment area [3rd], main channel slope [5th]) and climate covariates (mean monthly air temperature [2nd] and discharge [4th]). Two management covariates (groundwater use [6th] and riparian shade [7th]) also had high relative importance. Across the growing season (for all basins), local reach slope had high relative importance in May, but transitioned to a regional main channel slope covariate in August and September. This modeling process identified regionally similar and locally unique relationships among drivers of stream temperature. High relative importance of management-related covariates suggested potential restoration actions for each system.
River temperatures are expected to increase this century harming species requiring cold‐water habitat unless restoration activities protect or improve habitat availability. Local shading by riparian vegetation can cool water temperatures, but uncertainty exists over the scaling of this local effect to larger spatial extents. We evaluate this issue using a regional spatial stream network temperature model with covariates representing shade effects to predict mean August stream temperatures across 78,195 km of tributaries flowing into the Columbia River in the northwestern United States. We evaluate nine scenarios predicting stream temperatures for three riparian shade conditions (current, restored, and no riparian vegetation) within three different climate periods (2000s, 2040s, and 2080s). Results suggest riparian shade restoration (2000s climate) could decrease mean August stream temperatures by 0.62°C across the study network. Under the same restored shade conditions, temperature predictions for tributaries at their confluence with the Columbia River range from 0.02 to 2.08°C cooler than under current shade conditions. The climate warming effect predicted for the 2040s and 2080s, however, is greater than the cooling effect from restoring riparian shade. Streams less than 10‐m bankfull width cooled more frequently with riparian shade restoration. In Oregon, the proportion of fish habitat for salmon and trout rearing and migration that meet temperature numeric water quality criteria could be increased by 20% under restored shade conditions although net habitat declines may still occur in the future. We conclude riparian vegetation restoration could partially mitigate future warming and help maintain cold‐water habitats that function as thermal refuges if implemented strategically.
Ecosystems in the Anthropocene face pressures from multiple, interacting forms of environmental change. These pressures, resulting from land use change, altered hydrologic regimes, and climate change, will likely change the synchrony of ecosystem processes as distinct components of ecosystems are impacted in different ways. However, discipline-specific definitions and ad hoc methods for identifying synchrony and asynchrony have limited broader synthesis of this concept among studies and across disciplines. Drawing on concepts from ecology, hydrology, geomorphology, and biogeochemistry, we offer a unifying definition of synchrony for ecosystem science and propose a classification framework for synchrony and asynchrony of ecosystem processes. This framework classifies the relationships among ecosystem processes according to five key aspects: (1) the focal variables or relationships representative of the ecosystem processes of interest, (2) the spatial and temporal domain of interest, (3) the structural attributes of drivers and focal processes, (4) consistency in the relationships over time, and (5) the degree of causality among focal processes. Using this classification framework, we identify and differentiate types of synchrony and asynchrony, thereby providing the basis for comparing among studies and across disciplines. We apply this classification framework to existing studies in the ecological, hydrologic, geomorphic, and biogeochemical literature and discuss potential analytical tools that can be used to quantify synchronous and asynchronous processes. Furthermore, we seek to promote understanding of how different types of synchrony or asynchrony may shift in response to ongoing environmental change by providing a universal definition and explicit types and drivers with this framework.
Under a warmer future climate, thermal refuges could facilitate the persistence of species relying on cold-water habitat. Often these refuges are small and easily missed or smoothed out by averaging in models. Thermal infrared (TIR) imagery can provide empirical water surface temperatures that capture these features at a high spatial resolution (<1 m) and over tens of kilometers. Our study examined how TIR data could be used along with spatial stream network (SSN) models to characterize thermal regimes spatially in the Middle Fork John Day (MFJD) River mainstem (Oregon, USA). We characterized thermal variation in seven TIR longitudinal temperature profiles along the MFJD mainstem and compared them with SSN model predictions of stream temperature (for the same time periods as the TIR profiles). TIR profiles identified reaches of the MFJD mainstem with consistently cooler temperatures across years that were not consistently captured by the SSN prediction models. SSN predictions along the mainstem identified ~80% of the 1-km reach scale temperature warming or cooling trends observed in the TIR profiles. We assessed whether landscape features (e.g., tributary junctions, valley confinement, geomorphic reach classifications) could explain the fine-scale thermal heterogeneity in the TIR profiles (after accounting for the reach-scale temperature variability predicted by the SSN model) by fitting SSN models using the TIR profile observation points. Only the distance to the nearest upstream tributary was identified as a statistically significant landscape feature for explaining some of the thermal variability in the TIR profile data. When combined, TIR data and SSN models provide a data-rich evaluation of stream temperature captured in TIR imagery and a spatially extensive prediction of the network thermal diversity from the outlet to the headwaters.
The success of species reintroductions can depend on a combination of environmental, demographic, and genetic factors. Although the importance of these factors in the success of reintroductions is well-accepted, they are typically evaluated independently, which can miss important interactions. For species that persist in metapopulations, movement through and interaction with the landscape is predicted to be a vital component of persistence. Simulation-based approaches are a promising technique for evaluating the independent and combined effects of these factors on the outcome of various reintroduction and associated management actions. We report results from a simulation study of bull trout (Salvelinus confluentus) reintroduction to three watersheds of the Pend Oreille River system in northeastern Washington State, USA. We used an individual-based, spatially explicit simulation model to evaluate how reintroduction strategies, life history variation, and riverscape structure (e.g., network topology) interact to influence the demographic and genetic characteristics of reintroduced bull trout populations in three watersheds. Simulation scenarios included a range of initial genetic stocks (informed by empirical bull trout genetic data), variation in migratory tendency and life history, and two landscape connectivity alternatives representing a connected network (isolation-by-distance) and a fragmented network (isolation-by-barrier, using the known existing barriers). A novel feature of these simulations was the ability to consider the interaction of both demographic and genetic (i.e., demogenetic) factors in riverscapes with implicit asymmetric movement probabilities across the barriers. We found that connectivity (presence or absence of barriers) had the largest effect on demographic and genetic outcomes over 200 yr, with a greater effect than both initial genetic diversity and life history variation. We also identified regions of the study system in which bull trout populations persisted across a wide range of demographic, life history, and environmental connectivity parameters. Finally, we found no evidence that initial neutral genetic diversity influenced genetic diversity and structure after 200 yr; instead, genetic drift due to stray rate and population isolation dominated and erased any initial differences in genetic diversity. Our results highlight the utility of spatially explicit demogenetic approaches in exploring and understanding population dynamics and their implications for management strategies-in fresh waters.
Existing publicly available digital elevation models (DEMs) provide global-scale data but are often not precise enough for studying processes that depend on small-scale topographic features in rivers. For example, slope breaks and knickpoints in rivers can be important in understanding tectonic processes, and riffle-pool structures are important drivers of riverine ecology. More precise data (e.g. lidar) are available in some areas, but their spatial extent limits large-scale research. The upcoming Surface Water and Ocean Topography (SWOT) satellite mission is planned to launch in 2021 and will provide measurements of elevation and inundation extent of surface waters between 78° north and south latitude on average twice every 21 days. We present a novel noise reduction method for multitemporal river water surface elevation profiles from SWOT that combines a truncated singular value decomposition and a slope-constrained least-squares estimator. We use simulated SWOT data of 85-145 km sections of the Po, Sacramento, and Tanana Rivers to show that 3-12 months of simulated SWOT data can produce elevation profiles with mean absolute errors of 5.38-12.55 cm at 100-200 m along-stream resolution. Mean absolute errors can be reduced further to 4-11 cm by averaging all observations. The average profiles have errors much lower than existing DEMs, allowing new advances in riverine research globally. We consider two case studies in geomorphology and ecology that highlight the scientific value of the more accurate in-river DEMs expected from SWOT. Simulated SWOT elevation profiles for the Po reveal convexities in the river longitudinal profile that are spatially coincident with the upward projection of blind thrust faults that are buried beneath the Po Plain at the northern termination of the Apennine Mountains. Meanwhile, simulated SWOT data for the Sacramento River reveals locally steep sections of the river profile that represent important habitat for benthic invertebrates at a spatial scale previously unrecognizable in large-scale digital elevation models presently available for this river.
Habitat fragmentation restricts the movement of individuals across a landscape. In terrestrial and aquatic systems, barriers to movement can modify population and community dynamics at local or regional scales. This study contrasted life history traits related to lifespan with habitat fragmentation to determine impacts on species population genetic structure in the Neuse River Basin, USA. For this, we simulated gene flow among evenly-spaced populations in a river network and tracked individual and population genetics for 200 years. The modeled scenarios represent a full cross between five life history strategies and four riverscapes representing varying degrees of fragmentation. The five life history strategies include species (based on freshwater mussels) with average lifespans ranging from 10 to 50 years and age at maturity from 2 to 6 years. The movement landscapes included a (1) panmictic, (2) stepping-stone landscape allowing movement to only neighboring populations during each dispersal event, (3) partially-fragmented landscape divided by dams currently in the network, and (4) fully-fragmented landscape. Results suggest species with shorter lifespans have higher population genetic structure in fragmented landscapes than species with longer lifespans. Furthermore, species with shorter lifespans in highly fragmented landscapes may be harboring genetic degradation or decline as allele fixation and loss. Although anthropogenic fragmentation of many river systems is only 100–200 years old, the simulation indicates that species can respond genetically in that period of time. Additionally, the time frame of the simulation suggests that genetic impacts of habitat fragmentation in some species present in the Neuse River Basin may not yet be manifesting and restoration activities could be successful.
Increases in river fragmentation globally threaten freshwater biodiversity. Rivers are fragmented by many agents, both natural and anthropogenic. We review the distribution and frequency of these major agents, along with their effects on connectivity and habitat quality. Most fragmentation research has focused on terrestrial habitats, but theories and generalizations developed in terrestrial habitats do not always apply well to river networks. For example, terrestrial habitats are usually conceptualized as two-dimensional, whereas rivers often are conceptualized as one-dimensional or dendritic. In addition, river flow often leads to highly asymmetric effects of barriers on habitat and permeability. New approaches tailored to river networks can be applied to describe the network-wide effects of multiple barriers on both connectivity and habitat quality. The net effects of anthropogenic fragmentation on freshwater biodiversity are likely underestimated, because of time lags in effects and the difficulty of generating a single, simple signal of fragmentation that applies to all aquatic species. We conclude by presenting a decision tree for managing freshwater fragmentation, as well as some research horizons for evaluating fragmented riverscapes.
Increases in river fragmentation globally threaten freshwater biodiversity. Rivers are fragmented by many agents, both natural and anthropogenic. We review the distribution and frequency of these major agents, along with their effects on connectivity and habitat quality. Most fragmentation research has focused on terrestrial habitats, but theories and generalizations developed in terrestrial habitats do not always apply well to river networks. For example, terrestrial habitats are usually conceptualized as two‐dimensional, whereas rivers often are conceptualized as one‐dimensional or dendritic. In addition, river flow often leads to highly asymmetric effects of barriers on habitat and permeability. New approaches tailored to river networks can be applied to describe the network‐wide effects of multiple barriers on both connectivity and habitat quality. The net effects of anthropogenic fragmentation on freshwater biodiversity are likely underestimated, because of time lags in effects and the difficulty of generating a single, simple signal of fragmentation that applies to all aquatic species. We conclude by presenting a decision tree for managing freshwater fragmentation, as well as some research horizons for evaluating fragmented riverscapes.
Differences among lake morphologies often explain variation in characteristics of lentic ecosystems. Although beaver ponds also vary in morphology, previous studies have not examined the effects of such variation on downstream ecosystems. This study evaluated downstream effects of multiple beaver ponds in the Colorado Rocky Mountains during one low and one high-flow year. Beaver pond morphology was described as the natural log transformed ratio of beaver dam height (which determines hydraulic head) to pond surface area and related to pond spillover phytoplankton and characteristics of the ecosystem downstream (nutrient concentrations, limiting nutrients, periphyton, benthic organic matter (BOM), and benthic invertebrate consumers). Nitrate concentration increased systematically downstream of beaver ponds, but only in the low flow year when groundwater influences predominated. Effects of beaver ponds on soluble reactive phosphorus concentration depended on pond morphology, increasing downstream of small ponds with high dams, but only during the low-flow year. In situ experiments showed that neither beaver activity nor pond morphology predicted periphyton-limiting nutrients downstream. Both periphyton biomass and BOM decreased downstream of small ponds with high dams but pond morphology did not predict abundance of invertebrate grazers or detritus-feeding consumers. While suspension feeding invertebrates increased downstream from small ponds with high dams, variation in chlorophyll a from water spilling over beaver dams did not follow a similar pattern. We conclude that the effects of beaver ponds on downstream nutrients, resources and consumers are rarely systematic, but instead depend on variation in pond morphology and on annual hydrologic variation.
P>1. The North American beaver has been studied as a model ecosystem engineer for many decades. Previous studies have documented physical, chemical and biological impacts attributed to beaver engineering in both aquatic and terrestrial environments. This study focused on the effects of ecosystem engineering by beavers on life histories of a common mayfly and on the potential consequences for mayfly populations.2. We studied 18 montane beaver ponds of varying size and shape in western Colorado near the Rocky Mountain Biological Laboratory. Our goal was to test whether variation in beaver pond morphology (pond size and shape) explains downstream changes in stream temperature, mayfly size and timing of emergence.3. Downstream water temperatures varied predictably with pond morphology, being colder downstream of high-head dams and warmer downstream of low-head dams. Pond morphology was also a significant predictor of variation in the size of mature female Baetis bicaudatus (the most abundant mayfly), with larger females emerging downstream of high-head dams and smaller females downstream of low-head dams. The size of male B. bicaudatus was not significantly related to pond morphology or stream temperature. There was no relationship between pond morphology and variation in the timing of emergence of Baetis (males or females) between upstream and downstream reaches.4. Our results have implications for the effects of beaver ponds on Baetis individual fitness because large Baetis females are more fecund. Therefore, predictable female size variation associated with beaver pond morphology makes it possible to model the effects of beaver activity on local contributions of Baetis to the regional pool of reproductive adults at the catchment scale. Additionally, predictable changes in the size of emerging mayflies may have important consequences for the magnitude of aquatic to terrestrial resource subsidies in beaver-modified systems.
THE ISSUE Ecological processes and patterns interact at various scales across landscapes. Spatially explicit consideration of pattern and process can better inform ecological questions. This Figure Set introduces students to the variability in scales of pattern and process and to the fundamentals of Landscape Ecology. STUDENT-ACTIVE APPROACHES Think-pair-share, sketching conceptual figures STUDENT ASSESSMENTS sharing sketches of conceptual figures with class, informal short answer synthesis and application questions ACKNOWLEDGEMENTS We would like to thank Dean Urban for his initial encouragement to develop this teaching module, as well as comments on previous drafts. We also thank two anonymous reviewers for their helpful feedback and suggestions to improve the Figure Set. In addition, we thank Justin Wright for allowing us to trial our teaching module with his undergraduate ecology class. OVERVIEW WHAT IS THE ECOLOGICAL ISSUE? Landscape Ecology seeks to understand the interactions between ecological patterns and processes across varying spatial scales (Turner 2005). In some landscapes, extreme gradients drive patterns through their effect on physical processes. For example, the large elevation gradients in mountainous landscapes influence temperature and moisture, thereby causing plant species to sort into communities along the gradient according to their ability to grow and compete at different temperatures and soil moistures. However, subtle gradients also occur frequently in nature and can similarly influence variation in the physical and ecological variables that influence landscape pattern. To appreciate variability in the spatial scales of patterns and processes, this introductory lesson in Landscape Ecology compares two landscapes with large and small elevation gradients. This figure set introduces the field of Landscape Ecology by asking students: How do ecological processes and landscape patterns interact? Three sections ask students to examine how three types of drivers influence landscape pattern: the physical template, biotic interactions, and disturbance regimes. Using the Sierra Nevada Mountains and the Florida Scrub as examples, students explore how landscape patterns and ecological processes interact. Class size: Small Time: Moderate (time requirements include reading to be completed before class, an activity designed to take one 50 minute class period, and possibly a short follow-up as homework) FIGURE SET BACKGROUND This figure set is intended for undergraduate students nearing the end of an introductory ecology course. The reading and exercises draw on ecological concepts students should already know; the synthesis and application of these concepts are discussed in new ways to give a spatially explicit context to-3-TIEE understanding ecology. A glossary is supplied …