Network analysis reveals multiscale controls on streamwater chemistry

Watershed features shape spatial patterns of fish tissue mercury in a boreal river network

Freshwater mercury (Hg) contamination is a widespread environmental concern but how proximate sources and downstream transport shape Hg spatial patterns in riverine food webs is poorly understood. We measured total Hg (THg) in slimy sculpin (Cottus cognatus) across the Kuskokwim River, a large boreal river in western Alaska and home to subsistence fishing communities which rely on fish for primary nutrition. We used spatial stream network models (SSNMs) to quantify watershed and instream conditions influencing sculpin THg. Spatial covariates for local watershed geology and slope accounted for 55 % of observed variation in sculpin THg and evidence for downstream transport of Hg in sculpins was weak. Empirical semivariograms indicated these spatial covariates accounted for most spatial autocorrelation in observed THg. Watershed geology and slope explained up to 70 % of sculpin THg variation when SSNMs accounted for instream spatial dependence. Our results provide network-wide predictions for fish tissue THg based largely on publicly available geospatial data and open-source software for SSNMs, and demonstrate how these emerging models can be used to understand contaminant behavior in spatially complex aquatic ecosystems.

Hydrologic Connectivity Regulates Riverine N2O Sources and Dynamics

Indirect nitrous oxide (N2O) emissions from streams and rivers are a poorly constrained term in the global N2O budget. Current models of riverine N2O emissions place a strong focus on denitrification in groundwater and riverine environments as a dominant source of riverine N2O, but do not explicitly consider direct N2O input from terrestrial ecosystems. Here, we combine N2O isotope measurements and spatial stream network modeling to show that terrestrial-aquatic interactions, driven by changing hydrologic connectivity, control the sources and dynamics of riverine N2O in a mesoscale river network within the U.S. Corn Belt. We find that N2O produced from nitrification constituted a substantial fraction (i.e., >30%) of riverine N2O across the entire river network. The delivery of soil-produced N2O to streams was identified as a key mechanism for the high nitrification contribution and potentially accounted for more than 40% of the total riverine emission. This revealed large terrestrial N2O input implies an important climate-N2O feedback mechanism that may enhance riverine N2O emissions under a wetter and warmer climate. Inadequate representation of hydrologic connectivity in observations and modeling of riverine N2O emissions may result in significant underestimations.

Outbreeding reduces survival during metamorphosis in a headwater stream salamander

Assessing direct fitness effects of individual genetic diversity is challenging due to the intensive and long-term data needed to quantify survival and reproduction in the wild. But resolving these effects is necessary to determine how inbreeding and outbreeding influence eco-evolutionary processes. We used 8 years of capture-recapture data and single nucleotide polymorphism genotypes for 1906 individuals to test for effects of individual heterozygosity on stage-specific survival probabilities in the salamander Gyrinophilus porphyriticus. The life cycle of G. porphyriticus includes an aquatic larval stage followed by metamorphosis into a semi-aquatic adult stage. In our study populations, the larval stage lasts 6-10 years, metamorphosis takes several months, and lifespan can reach 20 years. Previous studies showed that metamorphosis is a sensitive life stage, leading us to predict that fitness effects of individual heterozygosity would occur during metamorphosis. Consistent with this prediction, monthly probability of survival during metamorphosis declined with multi-locus heterozygosity (MLH), from 0.38 at the lowest MLH (0.10) to 0.06 at the highest MLH (0.38), a reduction of 84%. Body condition of larvae also declined significantly with increasing MLH. These relationships were consistent in the three study streams. With evidence of localised inbreeding within streams, these results suggest that outbreeding disrupts adaptations in pre-metamorphic and metamorphic individuals to environmental gradients along streams, adding to evidence that headwater streams are hotspots of microgeographic adaptation. Our results also underscore the importance of incorporating life history in analyses of the fitness effects of individual genetic diversity and suggest that metamorphosis and similar discrete life stage transitions may be critical periods of viability selection.

Stage-Specific Demographic Effects of Hydrologic Variation in a Stream Salamander

AbstractWe lack a strong understanding of how organisms with complex life histories respond to climate variation. Many stream-associated species have multistage life histories that are likely to influence the demographic consequences of floods and droughts. However, tracking stage-specific demographic responses requires high-resolution, long-term data that are rare. We used 8 years of capture-recapture data for the headwater stream salamander Gyrinophilus porphyriticus to quantify the effects of flooding and drying magnitude on stage-specific vital rates and population growth. Drying reduced larval recruitment but increased the probability of metamorphosis (i.e., adult recruitment). Flooding reduced adult recruitment but had no effect on larval recruitment. Larval and adult survival declined with flooding but were unaffected by drying. Annual population growth rates (λ) declined with flooding and drying. Lambda also declined over the study period (2012-2021), although mean λ was 1.0 over this period. Our results indicate that G. porphyriticus populations are resilient to hydrologic variation because of compensatory effects on recruitment of larvae versus adults (i.e., reproduction vs. metamorphosis). Complex life cycles may enable this resilience to climate variation by creating opportunities for compensatory demographic responses across stages. However, more frequent and intense hydrologic variation in the latter half of this study contributed to a decline in λ over time, suggesting that increasing environmental variability poses a threat even when demographic compensation occurs.

A multiscale analysis of the spatially heterogeneous relationships between non-point source pollution-related processes and their main drivers in Chaohu Lake watershed, China

A better understanding of the relationships between non-point source (NPS) pollution-related processes and their drivers will help to develop scientific watershed management measures. Although various studies have explored the drivers' impact on NPS pollution-related processes, quantitative knowledge of the properties within these relationships is still needed. This study uses the Integrated Valuation of Ecosystem Services and Trade-offs (InVEST) model to produce three related processes of NPS pollution, quick flow (QF), nitrogen export (NE), and sediment export (SE), in the upstream watershed of Chaohu Lake, China. The spatial distributions of QF, NE, and SE and their responses to multiple natural-socioeconomic drivers at nine spatial scales (1 km2, 10 km2, 20 km2, 30 km2, 50 km2, 75 km2, 100 km2, 200 km2, and town) were compared. The results showed that the spatial scale has little impact on the spatial distributions of NPS pollution-related processes. Across the nine scales, the socioeconomic drivers related to agricultural activities, area proportions of cultivated land (cultivated) and paddy field (paddy), have dominant impacts on NE, while the topographical drivers, the connectivity index (IC) and slope, have dominant impacts on both SE and QF. The magnitudes of single and paired natural-socioeconomic drivers' impacts on NPS pollution-related processes increase logarithmically or linearly with increasing spatial scale, but they tend to reach a stable threshold at a certain coarse scale. Our results emphasized the necessity and importance of embracing spatial scale effects in watershed water environmental management.

Differences in Corticosterone Release Rates of Larval Spring Salamanders (Gyrinophilus porphyriticus) in Response to Native Fish Presence

Invasive fish predators are an important factor causing amphibian declines and may have direct and indirect effects on amphibian survival. For example, early non-lethal exposure to these stressors may reduce survival in later life stages, especially in biphasic species. In amphibians, the glucocorticoid hormone corticosterone is released by the hypothalamo-pituitary-interrenal axis (HPI), as an adaptive physiological response to environmental stressors. The corticosterone response (baseline and response to acute stressors) is highly flexible and context dependent, and this variation can allow individuals to alter their phenotype and behavior with environmental changes, ultimately increasing survival. We sampled larvae of the spring salamander (Gyrinophilus porphyriticus) from two streams that each contained predatory brook trout (Slavelinus fontinalis) in the lower reaches and no predatory brook trout in the upper reaches. We measured baseline and stress-induced corticosterone release rates of larvae from the lower and upper reaches using a non-invasive water-borne hormone assay. We hypothesized that corticosterone release rates would differ between larvae from fish-present reaches and larvae from fish-free reaches. We found that baseline and stressor-induced corticosterone release rates were downregulated in larvae from reaches with fish predators. These results indicate that individuals from reaches with predatory trout are responding to fish predators by downregulating corticosterone while maintaining an active HPI axis. This may allow larvae more time to grow before metamorphosing, while also allowing them to physiologically respond to novel stressors. However, prolonged downregulation of corticosterone release rates can impact growth in post-metamorphic individuals.

Vulnerable Waters are Essential to Watershed Resilience

Watershed resilience is the ability of a watershed to maintain its characteristic system state while concurrently resisting, adapting to, and reorganizing after hydrological (for example, drought, flooding) or biogeochemical (for example, excessive nutrient) disturbances. Vulnerable waters include non-floodplain wetlands and headwater streams, abundant watershed components representing the most distal extent of the freshwater aquatic network. Vulnerable waters are hydrologically dynamic and biogeochemically reactive aquatic systems, storing, processing, and releasing water and entrained (that is, dissolved and particulate) materials along expanding and contracting aquatic networks. The hydrological and biogeochemical functions emerging from these processes affect the magnitude, frequency, timing, duration, storage, and rate of change of material and energy fluxes among watershed components and to downstream waters, thereby maintaining watershed states and imparting watershed resilience. We present here a conceptual framework for understanding how vulnerable waters confer watershed resilience. We demonstrate how individual and cumulative vulnerable-water modifications (for example, reduced extent, altered connectivity) affect watershed-scale hydrological and biogeochemical disturbance response and recovery, which decreases watershed resilience and can trigger transitions across thresholds to alternative watershed states (for example, states conducive to increased flood frequency or nutrient concentrations). We subsequently describe how resilient watersheds require spatial heterogeneity and temporal variability in hydrological and biogeochemical interactions between terrestrial systems and down-gradient waters, which necessitates attention to the conservation and restoration of vulnerable waters and their downstream connectivity gradients. To conclude, we provide actionable principles for resilient watersheds and articulate research needs to further watershed resilience science and vulnerable-water management.

Emergent dual scaling of riverine biodiversity

A prevailing paradigm suggests that species richness increases with area in a decelerating way. This ubiquitous power law scaling, the species-area relationship, has formed the foundation of many conservation strategies. In spatially complex ecosystems, however, the area may not be the sole dimension to scale biodiversity patterns because the scale-invariant complexity of fractal ecosystem structure may drive ecological dynamics in space. Here, we use theory and analysis of extensive fish community data from two distinct geographic regions to show that riverine biodiversity follows a robust scaling law along the two orthogonal dimensions of ecosystem size and complexity (i.e., the dual scaling law). In river networks, the recurrent merging of various tributaries forms fractal branching systems, where the prevalence of branching (ecosystem complexity) represents a macroscale control of the ecosystem's habitat heterogeneity. In the meantime, ecosystem size dictates metacommunity size and total habitat diversity, two factors regulating biodiversity in nature. Our theory predicted that, regardless of simulated species' traits, larger and more branched "complex" networks support greater species richness due to increased space and environmental heterogeneity. The relationships were linear on logarithmic axes, indicating power law scaling by ecosystem size and complexity. In support of this theoretical prediction, the power laws have consistently emerged in riverine fish communities across the study regions (Hokkaido Island in Japan and the midwestern United States) despite hosting different fauna with distinct evolutionary histories. The emergence of dual scaling law may be a pervasive property of branching networks with important implications for biodiversity conservation.

Riverscape approaches in practice: perspectives and applications

Landscape perspectives in riverine ecology have been undertaken increasingly in the last 30 years, leading aquatic ecologists to develop a diverse set of approaches for conceptualizing, mapping and understanding 'riverscapes'. Spatiotemporally explicit perspectives of rivers and their biota nested within the socio-ecological landscape now provide guiding principles and approaches in inland fisheries and watershed management. During the last two decades, scientific literature on riverscapes has increased rapidly, indicating that the term and associated approaches are serving an important purpose in freshwater science and management. We trace the origins and theoretical foundations of riverscape perspectives and approaches and examine trends in the published literature to assess the state of the science and demonstrate how they are being applied to address recent challenges in the management of riverine ecosystems. We focus on approaches for studying and visualizing rivers and streams with remote sensing, modelling and sampling designs that enable pattern detection as seen from above (e.g. river channel, floodplain, and riparian areas) but also into the water itself (e.g. aquatic organisms and the aqueous environment). Key concepts from landscape ecology that are central to riverscape approaches are heterogeneity, scale (resolution, extent and scope) and connectivity (structural and functional), which underpin spatial and temporal aspects of study design, data collection and analysis. Mapping of physical and biological characteristics of rivers and floodplains with high-resolution, spatially intensive techniques improves understanding of the causes and ecological consequences of spatial patterns at multiple scales. This information is crucial for managing river ecosystems, especially for the successful implementation of conservation, restoration and monitoring programs. Recent advances in remote sensing, field-sampling approaches and geospatial technology are making it increasingly feasible to collect high-resolution data over larger scales in space and time. We highlight challenges and opportunities and discuss future avenues of research with emerging tools that can potentially help to overcome obstacles to collecting, analysing and displaying these data. This synthesis is intended to help researchers and resource managers understand and apply these concepts and approaches to address real-world problems in freshwater management.

Citizen science reveals unexpected solute patterns in semiarid river networks

Human modification of water and nutrient flows has resulted in widespread degradation of aquatic ecosystems. The resulting global water crisis causes millions of deaths and trillions of USD in economic damages annually. Semiarid regions have been disproportionately affected because of high relative water demand and pollution. Many proven water management strategies are not fully implemented, partially because of a lack of public engagement with freshwater ecosystems. In this context, we organized a large citizen science initiative to quantify nutrient status and cultivate connection in the semiarid watershed of Utah Lake (USA). Working with community members, we collected samples from ~200 locations throughout the 7,640 km2 watershed on a single day in the spring, summer, and fall of 2018. We calculated ecohydrological metrics for nutrients, major ions, and carbon. For most solutes, concentration and leverage (influence on flux) were highest in lowland reaches draining directly to the lake, coincident with urban and agricultural sources. Solute sources were relatively persistent through time for most parameters despite substantial hydrological variation. Carbon, nitrogen, and phosphorus species showed critical source area behavior, with 10-17% of the sites accounting for most of the flux. Unlike temperate watersheds, where spatial variability often decreases with watershed size, longitudinal variability showed an hourglass shape: high variability among headwaters, low variability in mid-order reaches, and high variability in tailwaters. This unexpected pattern was attributable to the distribution of human activity and hydrological complexity associated with return flows, losing river reaches, and diversions in the tailwaters. We conclude that participatory science has great potential to reveal ecohydrological patterns and rehabilitate individual and community relationships with local ecosystems. In this way, such projects represent an opportunity to both understand and improve water quality in diverse socioecological contexts.

Emergent spatial patterns of competing benthic and pelagic algae in a river network: A parsimonious basin-scale modeling analysis

Algae, as primary producers in riverine ecosystems, are found in two distinct habitats: benthic and pelagic algae typically prevalent in shallow/small and deep/large streams, respectively. Over an entire river continuum, spatiotemporal patterns of the two algal communities reflect specificity in habitat preference determined by geomorphic structure, hydroclimatic controls, and spatiotemporal heterogeneity in nutrient loads from point- and diffuse-sources. By representing these complex interactions between geomorphic, hydrologic, geochemical, and ecological processes, we present here a new river-network-scale dynamic model (CⁿANDY) for pelagic (A) and benthic (B) algae competing for energy and one limiting nutrient (phosphorus, P). We used the urbanized Weser River Basin in Germany (7th-order; ~8.4 million population; ~46 K km²) as a case study and analyzed simulations for equilibrium mass and concentrations under steady median river discharge. We also examined P, A, and B spatial patterns in four sub-basins. We found an emerging pattern characterized by scaling of P and A concentrations over stream-order ω, whereas B concentration was described by three distinct phases. Furthermore, an abrupt algal regime shift occurred in intermediate streams from B dominance in ω≤3 to exclusive A presence in ω≥6. Modeled and long-term basin-scale monitored dissolved P concentrations matched well for ω>4, and with overlapping ranges in ω<3. Power-spectral analyses for the equilibrium P, A, and B mass distributions along hydrological flow paths showed stronger clustering compared to geomorphological attributes, and longer spatial autocorrelation distance for A compared to B. We discuss the implications of our findings for advancing hydro-ecological concepts, guiding monitoring, informing management of water quality, restoring aquatic habitat, and extending CⁿANDY model to other river basins.

Tundra wildfire triggers sustained lateral nutrient loss in Alaskan Arctic

Climate change is creating widespread ecosystem disturbance across the permafrost zone, including a rapid increase in the extent and severity of tundra wildfire. The expansion of this previously rare disturbance has unknown consequences for lateral nutrient flux from terrestrial to aquatic environments. Lateral loss of nutrients could reduce carbon uptake and slow recovery of already nutrient-limited tundra ecosystems. To investigate the effects of tundra wildfire on lateral nutrient export, we analyzed water chemistry in and around the 10-year-old  Anaktuvuk River fire scar in northern Alaska. We collected water samples from 21 burned and 21 unburned watersheds during snowmelt, at peak growing season, and after plant senescence in 2017 and 2018. After a decade of ecosystem recovery, aboveground biomass had recovered in burned watersheds, but overall carbon and nitrogen remained ~20% lower, and the active layer remained ~10% deeper. Despite lower organic matter stocks, dissolved organic nutrients were substantially elevated in burned watersheds, with higher flow-weighted concentrations of organic carbon (25% higher), organic nitrogen (59% higher), organic phosphorus (65% higher), and organic sulfur (47% higher). Geochemical proxies indicated greater interaction with mineral soils in watersheds with surface subsidence, but optical analysis and isotopes suggested that recent plant growth, not mineral soil, was the main source of organic nutrients in burned watersheds. Burned and unburned watersheds had similar δ¹⁵ N-NO₃ ^- , indicating that exported nitrogen was of preburn origin (i.e., not recently fixed). Lateral nitrogen flux from burned watersheds was 2- to 10-fold higher than rates of background nitrogen fixation and atmospheric deposition estimated in this area. These findings indicate that wildfire in Arctic tundra can destabilize nitrogen, phosphorus, and sulfur previously stored in permafrost via plant uptake and leaching. This plant-mediated nutrient loss could exacerbate terrestrial nutrient limitation after disturbance or serve as an important nutrient release mechanism during succession.

Forest management impacts on stream integrity at varying intensities and spatial scales: Do biological effects accumulate spatially?

The effects of forest harvesting on headwaters are quite well understood, yet our understanding of whether impacts accumulate or dissipate downstream is limited. To address this, we investigated whether several biotic indicators changed from smaller to larger downstream sites (n = 6) within three basins that had intensive, extensive or minimal forest management in New Brunswick (Canada). Biofilm biomass and grazer abundance significantly increased from upstream to downstream, whereas organic matter decomposition and the autotrophic index of biofilms decreased. However, some spatial trends differed among basins and indicated either cumulative (macroinvertebrate abundance, predator density, sculpin GSI) or dissipative (autotrophic index, cotton decomposition) effects downstream, potentially explained by sediment and nutrient dynamics related to harvesting. No such among-basin differences were observed for leaf decomposition, biofilm biomass, macroinvertebrate richness or sculpin condition. Additionally, results suggest that some of the same biological impacts of forestry observed in small headwaters also occurred in larger systems. Although the intensive and extensive basins had lower macroinvertebrate diversity, there were no other signs of biological impairment, suggesting that, overall, current best management practices protect biological integrity downstream despite abiotic effects.

Spatiotemporal dynamics of water sources in a mountain river basin inferred through δ2H and δ18O of water

In mountainous river basins of the Pacific Northwest, climate models predict that winter warming will result in increased precipitation falling as rain and decreased snowpack. A detailed understanding of the spatial and temporal dynamics of water sources across river networks will help illuminate climate change impacts on river flow regimes. Because the stable isotopic composition of precipitation varies geographically, variation in surface water isotope ratios indicates the volume-weighted integration of upstream source water. We measured the stable isotope ratios of surface water samples collected in the Snoqualmie River basin in western Washington over June and September 2017 and the 2018 water year. We used ordinary least squares regression and geostatistical Spatial Stream Network models to relate surface water isotope ratios to mean watershed elevation (MWE) across seasons. Geologic and discharge data was integrated with water isotopes to create a conceptual model of streamflow generation for the Snoqualmie River. We found that surface water stable isotope ratios were lowest in the spring and highest in the dry, Mediterranean summer, but related strongly to MWE throughout the year. Low isotope ratios in spring reflect the input of snowmelt into high elevation tributaries. High summer isotope ratios suggest that groundwater is sourced from low elevation areas and recharged by winter precipitation. Overall, our results suggest that baseflow in the Snoqualmie River may be resilient to predicted warming and subsequent changes to snowpack in the Pacific Northwest.

Forest management impacts on stream integrity at varying intensities and spatial scales: Do abiotic effects accumulate spatially?

Though effects of forest harvesting on small streams are well documented, little is known about the cumulative effects in downstream systems. The hierarchical nature and longitudinal connectivity of river networks make them fundamentally cumulative, but lateral and vertical connectivity and instream processes can dissipate the downstream transport of water and materials. To elucidate such effects, we investigated how a suite of abiotic indicators changed from small streams to larger downstream sites (n = 6) within three basins ranging in forest management intensity (intensive, extensive, minimal) in New Brunswick (Canada) in the summer and fall of 2017 and 2018. Inorganic sediments, the inorganic/organic ratios and water temperatures significantly increased longitudinally, whereas nutrients and the fluorescence index of dissolved organic carbon (DOC; indication of terrestrial source) decreased. However, some longitudinal trends differed across basins and indicated downstream cumulative (inorganic sediments, the inorganic/organic ratios and to a lesser extent DOC concentration and humification) as well as dissipative (temperatures, nutrients, organic sediments) effects of forest management. Overall, we found that the effects previously reported for small streams with managed forests also occur at downstream sites and suggest investigating whether different management practices can be used within the extensive basin to reduce these cumulative effects.

Variation in stream network relationships and geospatial predictions of watershed conductivity

Secondary salinization, the increase of anthropogenically-derived salts in freshwaters, threatens freshwater biota and ecosystems, drinking water supplies, and infrastructure. The various anthropogenic sources of salts and their locations in a watershed may result in secondary salinization of river and stream networks through multiple inputs. We developed a watershed predictive assessment to investigate the degree to which topology, land-cover, and land-use covariates affect stream specific conductivity (SC), a measure of salinity. We used spatial stream network models to predict SC throughout an Appalachian stream network in a watershed affected by surface coal mining. During high-discharge conditions, 8 to 44% of stream km in the watershed exceeded the SC benchmark of 300 μS/cm, which is meant to be protective of aquatic life in the Central Appalachian ecoregion. During low-discharge conditions, 96 to 100% of stream km exceeded the benchmark. The 2 different discharge conditions altered the spatial dependency of SC among the stream monitoring sites. During most low discharges, SC was a function of upstream-to-downstream network distances, or flow-connected distances, among the sites. Flow-connected distances are indicative of upstream dependencies affecting stream SC. During high discharge, SC was related to both flow-connected distances and flow-unconnected distances (i.e., distances between sites on different branches of the network). Flow-unconnected distances are indicative of processes on adjacent branches and their catchments affecting stream SC. With sites distributed from headwaters to the watershed outlet, the extent of impacts from secondary salinization could be better spatially predicted and assessed with spatial stream network models than with models assuming spatial independence. Importantly, the assessment also recognized the multi-scale spatial relationships that can occur between the landscape and stream network.

Assessing Contamination of Stream Networks near Shale Gas Development Using a New Geospatial Tool

Chemical spills in streams can impact ecosystem or human health. Typically, the public learns of spills from reports from industry, media, or government rather than monitoring data. For example, ∼1300 spills (76 ≥ 400 gallons or ∼1500 L) were reported from 2007 to 2014 by the regulator for natural gas wellpads in the Marcellus shale region of Pennsylvania (U.S.), a region of extensive drilling and hydraulic fracturing. Only one such incident of stream contamination in Pennsylvania has been documented with water quality data in peer-reviewed literature. This could indicate that spills (1) were small or contained on wellpads, (2) were diluted, biodegraded, or obscured by other contaminants, (3) were not detected because of sparse monitoring, or (4) were not detected because of the difficulties of inspecting data for complex stream networks. As a first step in addressing the last problem, we developed a geospatial-analysis tool, GeoNet, that analyzes stream networks to detect statistically significant changes between background and potentially impacted sites. GeoNet was used on data in the Water Quality Portal for the Pennsylvania Marcellus region. With the most stringent statistical tests, GeoNet detected 0.2% to 2% of the known contamination incidents (Na ± Cl) in streams. With denser sensor networks, tools like GeoNet could allow real-time detection of polluting events.

Elevation and spatial structure explain most surface-water isotopic variation across five Pacific Coast basins

The stable isotope ratios of stream water can be used to trace water sources within river basins; however, drivers of variation in water isotopic spatial patterns across basins must be understood before ecologically relevant and isotopically distinct water sources can be identified and this tool efficiently applied. We measured the isotope ratios of surface-water samples collected during summer low-flow across five basins in Washington and southeast Alaska (Snoqualmie, Green, Skagit, and Wenatchee Rivers, and Cowee Creek) and compared models (isoscapes) describing the spatial variation in surface-water isotope ratios across a range of hydraulic and climatic conditions. We found strong correlations between mean watershed (MWE) elevation and surface-water isotopic ratios on the windward west side of the Cascades and in Alaska, explaining 48-90% of variation in δ¹⁸O values. Conversely, in the Wenatchee basin, located leeward of the Cascade Range, MWE alone had no predicative power. The elevation relationship and predictive isoscapes varied between basins, even those adjacent to each other. Applying spatial stream network models (SSNMs) to the Snoqualmie and Wenatchee Rivers, we found incorporating Euclidean and flow-connected spatial autocovariance improved explanatory power. SSNMs improved the accuracy of river water isoscapes in all cases; however, their utility was greater for the Wenatchee basin, where covariates explained only a small proportion of total variation. Our study provides insights into why basinscale surface-water isoscapes may vary even in adjacent basins and the importance of incorporating spatial autocorrelation in isoscapes. For determining source water contributions to downstream waters, our results indicate that surface water isoscapes should be developed for each basin of interest.

Long-term survival probability, not current habitat quality, predicts dispersal distance in a stream salamander

Dispersal evolves as an adaptive mechanism to optimize individual fitness across the landscape. Specifically, dispersal represents a mechanism to escape fitness costs resulting from changes in environmental conditions. Decades of empirical work suggest that individuals use local habitat cues to make movement decisions, but theory predicts that dispersal can also evolve as a fixed trait, independent of local conditions, in environments characterized by a history of stochastic spatiotemporal variation. Until now, however, both conditional and fixed models of dispersal evolution have primarily been evaluated using emigration data (stay vs. leave), and not dispersal distances: a more comprehensive measure of dispersal. Our goal was to test whether conditional or fixed models of dispersal evolution predict variation in dispersal distance in the stream salamander Gyrinophilus porphyriticus. We quantified variation in habitat conditions using measures of salamander performance from 4 yr of spatially explicit, capture-mark-recapture (CMR) data across three headwater streams in the Hubbard Brook Experimental Forest in central New Hampshire, USA. We used body condition as an index of local habitat quality that individuals may use to make dispersal decisions, and survival probability estimated from multistate CMR models as an index of mortality risk resulting from the long-term history of environmental variation. We found that dispersal distances increased with declining survival probability, indicating that salamanders disperse further in risky environments. Dispersal distances were unrelated to spatial variation in body condition, suggesting that salamanders do not base dispersal distance decisions on local habitat quality. Our study provides the first empirical support for fixed models of dispersal evolution, which predict that dispersal evolves in response to a history of spatiotemporal environmental variation, rather than as a conditional response to current habitat conditions. More broadly, this study underscores the value of assessing alternative scales of environmental variation to gain a more complete and balanced understanding of dispersal evolution.

The Landscape Ecology of Rivers: from Patch-Based to Spatial Network Analyses

Purpose of Review

We synthesize recent methodological and conceptual advances in the field of riverscape ecology, emphasizing areas of synergy with current research in landscape ecology.

Recent Findings

Recent advances in riverscape ecology highlight the need for spatially explicit examinations of how network structure influences ecological pattern and process, instead of the simple linear (upstream-downstream) view. Developments in GIS, remote sensing, and computer technologies already offer powerful tools for the application of patch- and gradient-based models for characterizing abiotic and biotic heterogeneity across a range of spatial and temporal scales. Along with graph-based analyses and spatial statistical stream network models (i.e., geostatistical modelling), these approaches offer improved capabilities for quantifying spatial and temporal heterogeneity and connectivity relationships, thereby allowing for rigorous and high-resolution analyses of pattern, process, and scale relationships.

Summary

Spatially explicit network approaches are able to quantify and predict biogeochemical, hydromorphological, and ecological patterns and processes more precisely than models based on longitudinal or lateral riverine gradients alone. Currently, local habitat characteristics appear to be more important than spatial effects in determining population and community dynamics, but this conclusion may change with direct quantification of the movement of materials, energy, and organisms along channels and across ecosystem boundaries—a key to improving riverscape ecology. Coupling spatially explicit riverscape models with optimization approaches will improve land protection and water management efforts, and help to resolve the land sharing vs. land sparing debate.

Drivers and spatial structure of abiotic and biotic properties of lakes, wetlands, and streams at the national scale

Broad-scale studies have improved our ability to make predictions about how freshwater biotic and abiotic properties will respond to changes in climate and land use intensification. Further, fine-scaled studies of lakes, wetlands, or streams have documented the important role of hydrologic connections for understanding many freshwater biotic and abiotic processes. However, lakes, wetlands, and streams are typically studied in isolation of one another at both fine and broad scales. Therefore, it is not known whether these three freshwater types (lakes, wetlands, and streams) respond similarly to ecosystem and watershed drivers nor how they may respond to future global stresses. In this study, we asked, do lake, wetland, and stream biotic and abiotic properties respond to similar ecosystem and watershed drivers and have similar spatial structure at the national scale? We answered this question with three U.S. conterminous data sets of freshwater ecosystems. We used random forest (RF) analysis to quantify the multi-scaled drivers related to variation in nutrients and biota in lakes, wetlands, and streams simultaneously; we used semivariogram analysis to quantify the spatial structure of biotic and abiotic properties and to infer possible mechanisms controlling the ecosystem properties of these freshwater types. We found that abiotic properties responded to similar drivers, had large ranges of spatial autocorrelation, and exhibited multi-scale spatial structure, regardless of freshwater type. However, the dominant drivers of variation in biotic properties depended on freshwater type and had smaller ranges of spatial autocorrelation. Our study is the first to document that drivers and spatial structure differ more between biotic and abiotic variables than across freshwater types, suggesting that some properties of freshwater ecosystems may respond similarly to future global changes.

Hydrologic variability contributes to reduced survival through metamorphosis in a stream salamander

Changes in the amount, intensity, and timing of precipitation are increasing hydrologic variability in many regions, but we have little understanding of how these changes are affecting freshwater species. Stream-breeding amphibians-a diverse group in North America-may be particularly sensitive to hydrologic variability during aquatic larval and metamorphic stages. Here, we tested the prediction that hydrologic variability in streams decreases survival through metamorphosis in the salamander Gyrinophilus porphyriticus, reducing recruitment to the adult stage. Using a 20-y dataset from Merrill Brook, a stream in northern New Hampshire, we show that abundance of G. porphyriticus adults has declined by ∼50% since 1999, but there has been no trend in larval abundance. We then tested whether hydrologic variability during summers influences survival through metamorphosis, using capture-mark-recapture data from Merrill Brook (1999 to 2004) and from 4 streams in the Hubbard Brook Experimental Forest (2012 to 2014), also in New Hampshire. At both sites, survival through metamorphosis declined with increasing variability of stream discharge. These results suggest that hydrologic variability reduces the demographic resilience and adaptive capacity of G. porphyriticus populations by decreasing recruitment of breeding adults. They also provide insight on how increasing hydrologic variability is affecting freshwater species, and on the broader effects of environmental variability on species with vulnerable metamorphic stages.

A distance-performance trade-off in the phenotypic basis of dispersal

Across taxa, individuals vary in how far they disperse, with most individuals staying close to their origin and fewer dispersing long distances. Costs associated with dispersal (e.g., energy, risk) are widely believed to trade off with benefits (e.g., reduced competition, increased reproductive success) to influence dispersal propensity. However, this framework has not been applied to understand variation in dispersal distance, which is instead generally attributed to extrinsic environmental factors. We alternatively hypothesized that variation in dispersal distances results from trade-offs associated with other aspects of locomotor performance. We tested this hypothesis in the stream salamander Gyrinophilus porphyriticus and found that salamanders that dispersed farther in the field had longer forelimbs but swam at slower velocities under experimental conditions. The reduced swimming performance of long-distance dispersers likely results from drag imposed by longer forelimbs. Longer forelimbs may facilitate moving longer distances, but the proximate costs associated with reduced swimming performance may help to explain the rarity of long-distance dispersal. The historical focus on environmental drivers of dispersal distances misses the importance of individual traits and associated trade-offs among traits affecting locomotion.

Spatial and temporal variation of ecosystem properties at macroscales

Although spatial and temporal variation in ecological properties has been well-studied, crucial knowledge gaps remain for studies conducted at macroscales and for ecosystem properties related to material and energy. We test four propositions of spatial and temporal variation in ecosystem properties within a macroscale (1000 km's) extent. We fit Bayesian hierarchical models to thousands of observations from over two decades to quantify four components of variation - spatial (local and regional) and temporal (local and coherent); and to model their drivers. We found strong support for three propositions: (1) spatial variation at local and regional scales are large and roughly equal, (2) annual temporal variation is mostly local rather than coherent, and, (3) spatial variation exceeds temporal variation. Our findings imply that predicting ecosystem responses to environmental changes at macroscales requires consideration of the dominant spatial signals at both local and regional scales that may overwhelm temporal signals.

Advancing the use of molecular methods for routine freshwater macroinvertebrate biomonitoring – the need for calibration experiments

Over the last decade, steady advancements have been made in the use of DNA-based methods for detection of species in a wide range of ecosystems. This progress has culminated in molecular monitoring methods being employed for the detection of several species for enforceable management purposes of endangered, invasive, and illegally harvested species worldwide. However, the routine application of DNA-based methods to monitor whole communities (typically a metabarcoding approach) in order to assess the status of ecosystems continues to be limited. In aquatic ecosystems, the limited use is particularly true for macroinvertebrate communities. As part of the DNAqua-Net consortium, a structured discussion was initiated with the aim to identify potential molecular methods for freshwater macroinvertebrate community assessment and identify important knowledge gaps for their routine application. We focus on three complementary DNA sources that can be metabarcoded: 1) DNA from homogenised samples (bulk DNA), 2) DNA extracted from sample preservative (fixative DNA), and 3) environmental DNA (eDNA) from water or sediment. We provide a brief overview of metabarcoding macroinvertebrate communities from each DNA source and identify challenges for their application to routine monitoring. To advance the utilisation of DNA-based monitoring for macroinvertebrates, we propose an experimental design template for a series of methodological calibration tests. The template compares sources of DNA with the goal of identifying the effects of molecular processing steps on precision and accuracy. Furthermore, the same samples will be morphologically analysed, which will enable the benchmarking of molecular to traditional processing approaches. In doing so we hope to highlight pathways for the development of DNA-based methods for the monitoring of freshwater macroinvertebrates.

A spatial stream-network approach assists in managing the remnant genetic diversity of riparian forests

Quantifying the genetic diversity of riparian trees is essential to understand their chances to survive hydroclimatic alterations and to maintain their role as foundation species modulating fluvial ecosystem processes. However, the application of suitable models that account for the specific dendritic structure of hydrographic networks is still incipient in the literature. We investigate the roles of ecological and spatial factors in driving the genetic diversity of Salix salviifolia, an Iberian endemic riparian tree, across the species latitudinal range. We applied spatial stream-network models that aptly integrate dendritic features (topology, directionality) to quantify the impacts of multiple scale factors in determining genetic diversity. Based on the drift hypothesis, we expect that genetic diversity accumulates downstream in riparian ecosystems, but life history traits (e.g. dispersal patterns) and abiotic or anthropogenic factors (e.g. drought events or hydrological alteration) might alter expected patterns. Hydrological factors explained the downstream accumulation of genetic diversity at the intermediate scale that was likely mediated by hydrochory. The models also suggested upstream gene flow within basins that likely occurred through anemophilous and entomophilous pollen and seed dispersal. Higher thermicity and summer drought were related to higher population inbreeding and individual homozygosity, respectively, suggesting that increased aridity might disrupt the connectivity and mating patterns among and within riparian populations.

Roots Mediate the Effects of Snowpack Decline on Soil Bacteria, Fungi, and Nitrogen Cycling in a Northern Hardwood Forest

Rising winter air temperature will reduce snow depth and duration over the next century in northern hardwood forests. Reductions in snow depth may affect soil bacteria and fungi directly, but also affect soil microbes indirectly through effects of snowpack loss on plant roots. We incubated root exclusion and root ingrowth cores across a winter climate-elevation gradient in a northern hardwood forest for 29 months to identify direct (i.e., winter snow-mediated) and indirect (i.e., root-mediated) effects of winter snowpack decline on soil bacterial and fungal communities, as well as on potential nitrification and net N mineralization rates. Both winter snowpack decline and root exclusion increased bacterial richness and phylogenetic diversity. Variation in bacterial community composition was best explained by differences in winter snow depth or soil frost across elevation. Root ingrowth had a positive effect on the relative abundance of several bacterial taxonomic orders (e.g., Acidobacterales and Actinomycetales). Nominally saprotrophic (e.g., Saccharomycetales and Mucorales) or mycorrhizal (e.g., Helotiales, Russalales, Thelephorales) fungal taxonomic orders were also affected by both root ingrowth and snow depth variation. However, when grouped together, the relative abundance of saprotrophic fungi, arbuscular mycorrhizal fungi, and ectomycorrhizal fungi were not affected by root ingrowth or snow depth, suggesting that traits in addition to trophic mode will mediate fungal community responses to snowpack decline in northern hardwood forests. Potential soil nitrification rates were positively related to ammonia-oxidizing bacteria and archaea abundance (e.g., Nitrospirales, Nitrosomondales, Nitrosphaerales). Rates of N mineralization were positively and negatively correlated with ectomycorrhizal and saprotrophic fungi, respectively, and these relationships were mediated by root exclusion. The results from this study suggest that a declining winter snowpack and its effect on plant roots each have direct effects on the diversity and abundance of soil bacteria and fungal communities that interact to determine rates of soil N cycling in northern hardwood forests.

Spatial Convergence in Major Dissolved Ion Concentrations and Implications of Headwater Mining for Downstream Water Quality

Spatial patterns in major dissolved solute concentrations were examined to better understand impact of surface coal mining in headwaters on downstream water chemistry. Sixty sites were sampled seasonally from 2012 to 2014 in an eastern Kentucky watershed. Watershed areas (WA) ranged from 1.6 to 400.5 km² and were mostly forested (58%-95%), but some drained as much as 31% surface mining. Measures of total dissolved solutes and most component ions were positively correlated with mining. Analytes showed strong convergent spatial patterns with high variability in headwaters (<15 km² WA) that stabilized downstream (WA > 75 km²), indicating hydrologic mixing primarily controls downstream values. Mean headwater solute concentrations were a good predictor of downstream values, with % differences ranging from 0.55% (Na⁺) to 28.78% (Mg²⁺). In a mined scenario where all headwaters had impacts, downstream solute concentrations roughly doubled. Alternatively, if mining impacts to headwaters were minimized, downstream solute concentrations better approximated the 300 μS/cm conductivity criterion deemed protective of aquatic life. Temporal variability also had convergent spatial patterns and mined streams were less variable due to unnaturally stable hydrology. The highly conserved nature of dissolved solutes from mining activities and lack of viable treatment options suggest forested, unmined watersheds would provide dilution that would be protective of downstream aquatic life.

A patchy continuum? Stream processes show varied responses to patch- and continuum-based analyses

Many conceptual syntheses in ecology and evolution are undergirded by either a patch- or continuum-based model. Examples include gradualism and punctuated equilibrium in evolution, and edge effects and the theory of island biogeography in ecology. In this study, we sought to determine how patch- or continuum-based analyses could explain variation in concentrations of stream macronutrients and system metabolism, represented by measures of productivity and respiration rates, at the watershed scale across the Kanawha River Basin, USA. Using Strahler stream order (SSO; continuum) and functional process zone (FPZ; patch) as factors, we produced statistical models for each variable and compared model performance using likelihood ratio tests. Only one nutrient (i.e., PO 4 3 - ) responded better to patch-based analysis. Both models were significantly better than a null model for ecosystem respiration; however, neither outperformed the other. Importantly, in most cases, a combination model, including both SSO and FPZ, best described observed variation in the system. Our findings suggest that several patch- and continuum-based processes may simultaneously influence the concentration of macronutrients and system metabolism. Nutrient spiral- ing along a continuum and the patch mosaic of land cover may both alter macronutrients, for example. Similarly, increases in temperature and discharge associated with increasing SSO, as well as the differences in light availability and channel morphology associated with different FPZs, may influence system metabolism. For these reasons, we recommend a combination of patch- and continuum-based analyses when modeling, analyzing, and interpreting patterns in stream ecosystem parameters.

PHYSICAL AND CHEMICAL CONNECTIVITY OF STREAMS AND RIPARIAN WETLANDS TO DOWNSTREAM WATERS: A SYNTHESIS

Streams, riparian areas, floodplains, alluvial aquifers and downstream waters (e.g., large rivers, lakes, oceans) are interconnected by longitudinal, lateral, and vertical fluxes of water, other materials and energy. Collectively, these interconnected waters are called fluvial hydrosystems. Physical and chemical connectivity within fluvial hydrosystems is created by the transport of nonliving materials (e.g., water, sediment, nutrients, contaminants) which either do or do not chemically change (chemical and physical connections, respectively). A substantial body of evidence unequivocally demonstrates physical and chemical connectivity between streams and riparian wetlands and downstream waters. Streams and riparian wetlands are structurally connected to downstream waters through the network of continuous channels and floodplain form that make these systems physically contiguous, and the very existence of these structures provides strong geomorphologic evidence for connectivity. Functional connections between streams and riparian wetlands and their downstream waters vary geographically and over time, based on proximity, relative size, environmental setting, material disparity, and intervening units. Because of the complexity and dynamic nature of connections among fluvial hydrosystem units, a complete accounting of the physical and chemical connections and their consequences to downstream waters should aggregate over multiple years to decades.

Using synoptic tracer surveys to assess runoff sources in an Andean headwater catchment in central Chile

Headwater catchments in the Andes provide critical sources of water for downstream areas with large agricultural communities dependent upon irrigation. Data from such remote headwater catchments are sparse, and there is limited understanding of their hydrological function to guide sustainable water management. Here, we present the findings of repeat synoptic tracer surveys as rapid appraisal tools to understand dominant hydrological flow paths in the semi-arid Rio Grande basin, a 572-km² headwater tributary of the 11,696-km² Limarí basin in central Chile. Stable isotopes in stream water show a typical altitudinal effect, with downstream enrichment in δ²H and δ¹⁸O ratios. Seasonal signals are displayed in the isotopic composition of the springtime melting season water line with a steeper gradient, whilst evaporative effects are represented by lower seasonal gradients for autumn and summer. Concentrations of solutes indexed by electrical conductivity indicate that there are limited contributions of deeper mineralised groundwater to streamflow and that weathering rates vary in the different sub-catchments. Although simplistic, the insights gained from the study could be used to inform the structure and parameterisation of rainfall runoff models to provide seasonal discharge predictions as an evidence base for decision making in local water management.

IMPROVING PREDICTIVE MODELS OF IN-STREAM PHOSPHORUS CONCENTRATION BASED ON NATIONALLY-AVAILABLE SPATIAL DATA COVERAGES

Spatial data are playing an increasingly important role in watershed science and management. Large investments have been made by government agencies to provide nationally-available spatial databases; however, their relevance and suitability for local watershed applications is largely unscrutinized. We investigated how goodness of fit and predictive accuracy of total phosphorus (TP) concentration models developed from nationally-available spatial data could be improved by including local watershed-specific data in the East Fork of the Little Miami River, Ohio, a 1290 km² watershed. We also determined whether a spatial stream network (SSN) modeling approach improved on multiple linear regression (nonspatial) models. Goodness of fit and predictive accuracy were highest for the SSN model that included local covariates, and lowest for the nonspatial model developed from national data. Septic systems and point source TP loads were significant covariates in the local models. These local data not only improved the models but enabled a more explicit interpretation of the processes affecting TP concentrations than more generic national covariates. The results suggest that SSN modeling greatly improves prediction and should be applied when using national covariates. Including local covariates further increases the accuracy of TP predictions throughout the studied watershed; such variables should be included in future national databases, particularly the locations of septic systems.

Role of surface and subsurface processes in scaling N2O emissions along riverine networks

Riverine environments, such as streams and rivers, have been reported as sources of the potent greenhouse gas nitrous oxide ([Formula: see text]) to the atmosphere mainly via microbially mediated denitrification. Our limited understanding of the relative roles of the near-surface streambed sediment (hyporheic zone), benthic, and water column zones in controlling [Formula: see text] production precludes predictions of [Formula: see text] emissions along riverine networks. Here, we analyze [Formula: see text] emissions from streams and rivers worldwide of different sizes, morphology, land cover, biomes, and climatic conditions. We show that the primary source of [Formula: see text] emissions varies with stream and river size and shifts from the hyporheic-benthic zone in headwater streams to the benthic-water column zone in rivers. This analysis reveals that [Formula: see text] production is bounded between two [Formula: see text] emission potentials: the upper [Formula: see text] emission potential results from production within the benthic-hyporheic zone, and the lower [Formula: see text] emission potential reflects the production within the benthic-water column zone. By understanding the scaling nature of [Formula: see text] production along riverine networks, our framework facilitates predictions of riverine [Formula: see text] emissions globally using widely accessible chemical and hydromorphological datasets and thus, quantifies the effect of human activity and natural processes on [Formula: see text] production.

Water sources and mixing in riparian wetlands revealed by tracers and geospatial analysis

Mixing of waters within riparian zones has been identified as an important influence on runoff generation and water quality. Improved understanding of the controls on the spatial and temporal variability of water sources and how they mix in riparian zones is therefore of both fundamental and applied interest. In this study, we have combined topographic indices derived from a high-resolution Digital Elevation Model (DEM) with repeated spatially high-resolution synoptic sampling of multiple tracers to investigate such dynamics of source water mixing. We use geostatistics to estimate concentrations of three different tracers (deuterium, alkalinity, and dissolved organic carbon) across an extended riparian zone in a headwater catchment in NE Scotland, to identify spatial and temporal influences on mixing of source waters. The various biogeochemical tracers and stable isotopes helped constrain the sources of runoff and their temporal dynamics. Results show that spatial variability in all three tracers was evident in all sampling campaigns, but more pronounced in warmer dryer periods. The extent of mixing areas within the riparian area reflected strong hydroclimatic controls and showed large degrees of expansion and contraction that was not strongly related to topographic indices. The integrated approach of using multiple tracers, geospatial statistics, and topographic analysis allowed us to classify three main riparian source areas and mixing zones. This study underlines the importance of the riparian zones for mixing soil water and groundwater and introduces a novel approach how this mixing can be quantified and the effect on the downstream chemistry be assessed.

Validation and comparison of geostatistical and spline models for spatial stream networks

Scientists need appropriate spatial-statistical models to account for the unique features of stream network data. Recent advances provide a growing methodological toolbox for modelling these data, but general-purpose statistical software has only recently emerged, with little information about when to use different approaches. We implemented a simulation study to evaluate and validate geostatistical models that use continuous distances, and penalised spline models that use a finite discrete approximation for stream networks. Data were simulated from the geostatistical model, with performance measured by empirical prediction and fixed effects estimation. We found that both models were comparable in terms of squared error, with a slight advantage for the geostatistical models. Generally, both methods were unbiased and had valid confidence intervals. The most marked differences were found for confidence intervals on fixed-effect parameter estimates, where, for small sample sizes, the spline models underestimated variance. However, the penalised spline models were always more computationally efficient, which may be important for real-time prediction and estimation. Thus, decisions about which method to use must be influenced by the size and format of the data set, in addition to the characteristics of the environmental process and the modelling goals. ©2015 The Authors. Environmetrics published by John Wiley & Sons, Ltd.