Sherri L. Johnson

We examined the relationship between winter snow accumulation and summer stream flow and temperature in the Pacific Northwest. Several papers have hypothesized that climate change will reduce snow accumulation, thereby reducing summer stream flows, and leading to higher stream temperatures, especially in the late summer. Long-term data from gaged watersheds located in the H. J. Andrews Experimental Forest allowed us to test this hypothesis. Over 20 years, maximum snowpacks have been highly variable, with annual maximum snowpacks ranging from 103 to 875 mm, but summers are typically dry and often rain free. Thus, snowmelt would appear to have a high potential to influence summer streamflow. Our analyses did not support the hypothesized relationships. While the size of snowpack explained variability in spring stream temperatures and streamflow, this relationship did not persist into the summer. Summer stream temperatures were poorly explained by spring precipitation and snowmelt metrics, and were, alternatively, best explained by streamflow, precipitation, and air temperature metrics across all watersheds and elevations. Our findings suggest for Mediterranean climates in the Pacific Northwest, snowpack variability has a limited relationship to summer stream temperatures and late summer streamflow. Instead, late-summer stream conditions appeared to be most influenced by multiple factors, including precipitation after snowmelt has ceased.

Paired‐catchment studies are widely used to examine the effects of land management practices (“treatments”) on hydrologic processes. Catchments are matched and a pretreatment calibration regression is used to identify the hydrological relationship between the reference and treated catchments. This method assumes that the calibration regression represents the actual relationship between the catchments (assumption of representativeness) and that the relationship will remain stable over time (assumption of stability). Errors are assumed to be small and similar between reference and treated catchments. Thus, observed differences between the catchments following treatment are assumed to result from that treatment alone. However, calibration periods are often short and it is impossible to know if the calibration period is representative. Further, because the study is unreplicated, it is impossible to determine if stability is maintained. Consequently, it is difficult to determine a minimum detectable effect sizes (MDES) below which estimates of changes in streamflow are statistically uncertain. Here, we use bootstrapped sampling from reference‐by‐reference (RxR) comparisons in a paired‐catchment study design to evaluate the MDES. We generate frequency distributions of the potential changes in flow—changes that cannot be caused by treatment effects. From these, we estimate bootstrapped ±95% confidence intervals encompassing the non‐treatment effects which we use as the MDES. We apply this method to long‐term paired‐catchment studies and reexamine changes in both annual water yields and late summer low flows at the HJA Experimental Forest. This bootstrapping method is widely transferable to any long‐term paired catchment study sites where multiple reference catchments exist.

Riparian zones are a critical terrestrial‐aquatic ecotone. They play important roles in ecosystems including (1) harboring biodiversity, (2) influencing light and carbon fluxes to aquatic food webs, (3) maintaining water quality and streamflow, (4) enhancing aquatic habitat, (5) influencing greenhouse gas production, and (6) sequestering carbon. Defining what qualifies as a riparian zone is a first step to delineation. Many definitions of riparian boundaries focus on static attributes or a subset of potential functions without recognizing that they are spatially continuous, temporally dynamic, and multi‐dimensional. We emphasize that definitions should consider multiple ecological and biogeochemical functions and physical gradients, and explore how this approach influences spatial characterization of riparian zones. One or more of the following properties can guide riparian delineation: (1) distinct species, elevated biodiversity, or species with specific adaptations to flooding and inundation near streams relative to nearby upland areas; (2) unique vegetation structure directly influencing irradiance or organic material inputs to aquatic ecosystems; (3) hydrologic and geomorphic features or processes maintaining floodplains; (4) hydric soil properties that differ from the uplands; and/or (5) elevated retention of dissolved and suspended materials relative to adjacent uplands. Considering these properties for an operational and dynamic definition of riparian zones recognizes that riparian boundaries vary in space (e.g., variation of riparian corridor widths within or among watersheds) and time (e.g., responses to hydrological variance and climate change). Inclusive definitions addressing multiple riparian functions could facilitate attainment of research and management goals by linking properties of interest to specific outcomes.

Temporal patterns in chemistry of headwater streams reflect responses of water and elemental cycles to perturbations occurring at local to global scales. We evaluated multi-scale temporal patterns in up to 32 y of monthly observations of stream chemistry (ammonium, calcium, dissolved organic carbon, nitrate, total dissolved phosphorus, and sulfate) in 22 reference catchments within the northern temperate zone of North America. Multivariate autoregressive state-space (MARSS) models were applied to quantify patterns at multi-decadal, seasonal, and shorter intervals during a period that encompassed warming climate, seasonal changes in precipitation, and regional declines in atmospheric deposition. Significant long-term trends in solute concentrations within a subset of the catchments were consistent with recovery from atmospheric deposition (e.g., calcium, nitrate, sulfate) and increased precipitation (e.g., dissolved organic carbon). Lack of evidence for multi-decadal trends in most catchments suggests resilience of northern temperate ecosystems or that subtle net effects of simultaneous changes in climate and disturbance regimes do not result in directional trends. Synchronous seasonal oscillations of solute concentrations occurred across many catchments, reflecting shared climate and biotic drivers of seasonality within the northern temperate zone. Despite shared patterns among catchments at a seasonal scale, multi-scale temporal patterns were statistically distinct among even adjacent headwater catchments, implying that local attributes of headwater catchments modify the signals imparted by atmospheric phenomena and regional disturbances. To effectively characterize hydrologic and biogeochemical responses to changing climate and disturbance regimes, catchment monitoring programs could include multiple streams with contributing areas that encompass regional heterogeneity in vegetation, topography, and elevation. Overall, detection of long-term patterns and trends requires monitoring multiple catchments at a frequency that captures periodic variation (e.g., seasonality) and a duration encompassing the perturbations of interest.

Accurate mapping of headwater streams and their flow status has important implications for understanding and managing water resources and land uses. However, accurate information is rare, especially in rugged, forested terrain. We developed a streamflow permanence classification model for forested lands in western Oregon using the latest light detection and ranging‐derived hydrography published in the National Hydrography Dataset. Models were trained using 2,518 flow/no flow field observations collected in late summer 2019–2021 across headwaters of 129 sub‐watersheds. The final model, the Western Oregon WeT DRy model, used Random Forest and 13 environmental covariates for classifying every 5‐m stream sub‐reach across 426 sub‐watersheds. The most important covariates were annual precipitation and drainage area. Model output included probabilities of late summer surface flow presence and were subsequently categorized into three streamflow permanence classes—Wet, Dry, and Ambiguous. Ambiguous denoted model probabilities and associated prediction intervals that extended over the 50% classification threshold between wet and dry. Model accuracy was 0.83 for sub‐watersheds that contained training data and decreased to 0.67 for sub‐watersheds that did not have observations of late summer surface flow. The model identified where predictions extrapolated beyond the domain characterized by the training data. The combination of spatially continuous estimates of late summer streamflow status along with uncertainty and extrapolation estimates provide critical information for strategic project planning and designing additional field data collection.

The body size of aquatic vertebrates is declining across populations and ecosystems worldwide owing to warmer water temperature and changing streamflow. In freshwaters, the effects of stream network position and density‐dependent factors on body size are less understood. We used an extensive dataset spanning 41 stream sites over 7 years to evaluate how density‐dependent and density‐independent factors influence the size of two top predators in small watersheds, Coastal Cutthroat Trout Oncorhynchus clarkii clarkii and Coastal Giant Salamanders Dicamptodon tenebrosus. We tested three hypotheses of body‐size variation for trout and salamanders, including intraspecific density dependence, interspecific density dependence, and resource availability, using empirical observations in hierarchical linear mixed models in a model‐selection framework. In our best‐supported models, the strongest predictors of size were conspecific negative density dependence, as expected, suggesting greater intraspecific interactions probably owing to conspecific individuals having similar requirements. We reveal a biogeographic pattern in which body size peaks in middle stream‐network positions and plateaus or declines at lower and upper locations, proposing that stream network position also plays a role in determining body size in small watersheds. Salamander density also has a quadratic effect on adult trout size, with salamanders having a greater overall effect on the body size of both species than trout, suggesting that salamanders might be more dominant than trout in some interactions. Collectively, we found that biotic interactions, mainly conspecific but also interspecific, and stream‐network position affect trout and salamander body sizes in small watersheds.

Large‐scale wildfires are becoming increasingly common in the wet forests of the Pacific Northwest (USA), with predicted increases in fire prevalence under future climate scenarios. Wildfires can alter streamflow response to precipitation and mobilize water quality constituents, which pose a risk to aquatic ecosystems and downstream drinking water treatment. Research often focuses on the impacts of high‐severity wildfires, with stream biogeochemical responses to low‐ and mixed‐severity fires often understudied, particularly during seasonal shifts in hydrologic connectivity between hillslopes and streams. We studied the impacts of the 2020 Holiday Farm Fire at the HJ Andrews Experimental Forest where rare pre‐fire stream discharge and chemistry data allowed us to evaluate the influence of mixed‐severity fire on stream water quantity and quality. Our research design focused on two well‐studied watersheds with low and low‐moderate burn severity where we examined long‐term data (pre‐ and post‐fire), and instantaneous grab samples collected during four rain events occurring immediately following wildfire and a prolonged dry summer. We analysed the impact of these rain events, which represent the transition from low‐to‐high hydrologic connectivity of the subsurface to the stream, on stream discharge and chemistry behaviour. Long‐term data revealed total annual flows and mean flows remained fairly consistent post‐fire, while small increases in baseflow were observed in the low‐moderately burned watershed. Stream water concentrations of nitrate, phosphate and sulfate significantly increased following fire, with variance in concentration increasing with fire severity. Our end member mixing models suggested that during rain events, the watershed with low‐moderate severity fire had greater streamflow inputs from soil water and groundwater during times of low connectivity compared to the watershed with low severity fire. Finally, differences in fire severity impacts on concentration‐discharge relationships of biogenic solutes were most expressed under low catchment connectivity conditions. Our study provides insights into post‐wildfire impacts to stream water quality, with the goal of informing future research on stream chemistry responses to low, moderate and mixed severity wildfire.

Disturbances to forested watersheds often result in increases of nutrients and light to nearby streams. Such changes are generally expected to produce a shift to a more autotrophic aquatic ecosystem, with measurable increases in algae, and associated implications for food webs and fisheries. Although this paradigm is widely established, results from our 10-year study (2007–2016) in 12 headwater streams and four sites downstream in the Trask River Watershed (Oregon, USA), did not concur. In 2012, one watershed was thinned, three were clearcut harvested with variable buffers and three with uniform riparian buffers. After harvest, light to the stream surface significantly increased at the three watersheds with variable buffers while dissolved inorganic nitrogen (DIN) significantly increased in all of the clearcut harvested streams. Despite the increase in DIN and light, algal standing stocks and chlorophyllaconcentrations did not significantly increase. The common assumption of increased autotrophic responses in stream food webs following increases of nitrogen and light was not supported here. We postulate the co-limitation of nutrients, driven by low phosphorus concentrations, which unlike DIN did not increase post-harvest, and the characteristics of the algal community, which were dominated by low light adapted diatoms rather than green algae, contributed to our findings of no responses for standing stocks of epilithic algae or concentrations of chlorophylla. The inclusion of multiple statistical analyses provided more certainty around our findings. This study documents responses to current forest practices and provides cautionary information for management and restoration activities aiming to increase fish abundance and standing stocks by opening riparian canopies and adding nutrients.