Of Small Streams and Great Lakes: Integrating Tributaries to Understand the Ecology and Biogeochemistry of Lake Superior

Lake Superior receives inputs from approximately 2,800 tributaries that provide nutrients and dissolved organic matter (DOM) to the nearshore zone of this oligotrophic lake. Here, we review the magnitude and timing of tributary export and plume formation in Lake Superior, how these patterns and interactions may shift with global change, and how emerging technologies can be used to better characterize tributary–lake linkages. Peak tributary export occurs during snowmelt‐driven spring freshets, with additional pulses during rain‐driven storms. Instream processing and transformation of nitrogen, phosphorus, and dissolved organic carbon (DOC) can be rapid but varies seasonally in magnitude. Tributary plumes with elevated DOC concentration, higher turbidity, and distinct DOM character can be detected in the nearshore during times of high runoff, but plumes can be quickly transported and diluted by in‐lake currents and mixing. Understanding the variability in size and load of these tributary plumes, how they are transported within the lake, and how long they persist may be best addressed with environmental sensors and remote sensing using autonomous and unmanned vehicles. The connections between Lake Superior and its tributaries are vulnerable to climate change, and understanding and predicting future changes to these valuable freshwater resources will require a nuanced and detailed consideration of tributary inputs and interactions in time and space.     

River Basin Export Reduction Optimization Support Tool; a tool to screen options for reducing nutrient loads while minimizing cost

Excess loading of nitrogen and phosphorus to river networks causes environmental harm, but reducing loads from large river basins is difficult and expensive. We developed a new tool, the River Basin Export Reduction Optimization Support Tool (RBEROST) to identify the least-cost combinations of management practices that will reduce nutrient loading to target levels in downstream and mid-network waterbodies. We demonstrate the utility of the tool in a case study in the Upper Connecticut River Basin in New England, USA. The total project cost of optimized lowest-cost plans ranged from $18.0 million to $41.0 million per year over 15 years depending on user specifications. Plans include both point source and non-point source management practices, and most costs are associated with urban stormwater practices. Adding a 2% margin of safety to loading targets improved the estimated probability of success from 37.5% to 99%. The large spatial scale of RBEROST, and the consideration of both point and non-point source contributions of nutrients, make it well suited as an initial screening tool in watershed planning.     

Landscape Hydrogeology and its Influence on Patterns of Groundwater Flux and Nitrate Removal Efficiency in Riparian Buffers

This study uses data from 46 riparian sites to examine the influence of landscape hydrogeology on patterns of groundwater flux and the buffer width required for effective nitrate removal in humid temperate agricultural regions. There is a considerable imbalance in the research focus on different hydrogeologic settings. More than 40% of the buffers are located in landscapes with surficial sand aquifers, whereas few buffers have been studied in glacial till and weathered bedrock landscapes which cover large areas. Annual groundwater fluxes for 29 of these sites ranged from <20 L/m/day for buffers on flat sand plains and uplands with fine‐textured deposits to 50‐1,200 L/m/day for many sites with upland sand aquifers. Despite a similar range of water fluxes, buffers in gently to moderately sloping landscapes with <4 m depths of sand sediments reached a 90% removal efficiency within 30‐60 m while sites with >4 m depths required a 150‐200 m width. The width for 90% efficiency in buffers with loamy sand and sandy loam sediments also increased from 10‐20 m with <4 m sediment depths to 50‐100 m for >4 m depths. Limited data for buffers with fine‐textured sediments suggest that 90% of the nitrate flux was often depleted in a 10‐20 m width. Groundwater flux did not have a significant relationship with nitrate removal percent per meter buffer width because of the variation in efficiency that occurred in buffers with similar fluxes in different hydrogeologic settings.     

Connectivity of Streams and Wetlands to Downstream Waters: An Integrated Systems Framework

Interest in connectivity has increased in the aquatic sciences, partly because of its relevance to the Clean Water Act. This paper has two objectives: (1) provide a framework to understand hydrological, chemical, and biological connectivity, focusing on how headwater streams and wetlands connect to and contribute to rivers; and (2) briefly review methods to quantify hydrological and chemical connectivity. Streams and wetlands affect river structure and function by altering material and biological fluxes to the river; this depends on two factors: (1) functions within streams and wetlands that affect material fluxes; and (2) connectivity (or isolation) from streams and wetlands to rivers that allows (or prevents) material transport between systems. Connectivity can be described in terms of frequency, magnitude, duration, timing, and rate of change. It results from physical characteristics of a system, e.g., climate, soils, geology, topography, and the spatial distribution of aquatic components. Biological connectivity is also affected by traits and behavior of the biota. Connectivity can be altered by human impacts, often in complex ways. Because of variability in these factors, connectivity is not constant but varies over time and space. Connectivity can be quantified with field‐based methods, modeling, and remote sensing. Further studies using these methods are needed to classify and quantify connectivity of aquatic ecosystems and to understand how impacts affect connectivity.     

Relative Importance of Water‐Quality Stressors in Predicting Fish Community Responses in Midwestern Streams

Fish, habitat, and water chemistry data were collected from 98 streams in the midwestern United States, an area dominated by intense cultivation of row crops, in order to identify important water‐quality stressors to fish communities. We focused on 10 stressors including riparian disturbance, riparian vegetative cover, instream fish cover, streambed sedimentation, streamflow variability, total nitrogen, total phosphorus, minimum dissolved oxygen, pesticides, and bed sediment contaminants. Fish community response variables included a measure of observed/expected taxonomic completeness; species‐specific tolerances to nitrogen, phosphorus, dissolved oxygen, and water temperature; the percent of species classified as macrohabitat generalists; and an index of pesticide toxicity to fish. Multivariate analysis indicated that total nitrogen was the most important stressor, signifying that fish communities were responding to total nitrogen despite relatively high levels common to an agricultural setting. Individually, fish taxonomic completeness decreased with increasing streambed sedimentation, whereas fish community tolerance to total phosphorus increased with increasing streambed sedimentation, riparian disturbance, and total nitrogen. These findings underscore the importance of multiple biological response metrics to better understand the effects of water‐quality stressors on fish communities and highlight the complex relations between total phosphorus and fish communities.     

The Impact of Asynchronicity on Event‐Flow Estimation in Basin‐Scale Hydrologic Model Calibration1

Abstract: The calibration of basin‐scale hydrologic models consists of adjusting parameters such that simulated values closely match observed values. However, due to inevitable inaccuracies in models and model inputs, simulated response hydrographs for multiyear calibrations will not be perfectly synchronized with observed response hydrographs at the daily time step. An analytically derived formula suggests that when timing errors are significant, traditional calibration approaches may generally underestimate the total event‐flow volume. An event‐adaptive time series is developed and incorporated into the Nash‐Sutcliffe Efficiency objective function to diagnose the potential impact of event‐flow synchronization errors. Test sites are the 50 km² Subwatershed I of the Little River Experimental Watershed (LREWswI) in southeastern Georgia, and the 610 km² Little Washita River Experimental Watershed (LWREW) in southwestern Oklahoma, with the Soil and Water Assessment Tool used as the hydrologic model. Results suggest that simulated surface runoff generation is 55% less for LREWswI when the daily time series is used compared with when the event‐adaptive technique is used. Event‐flow generation may also be underestimated for LWREW, but to a lesser extent than it may be for LREWswI, due to a larger portion of the event flow being lateral flow.     

Predicting the Planform Configuration of the Braided Toklat River,Alaska, With a Suite of Rule‐Based Models1

Abstract: An ensemble of rule‐based models was constructed to assess possible future braided river planform configurations for the Toklat River in Denali National Park and Preserve, Alaska. This approach combined an analysis of large‐scale influences on stability with several reduced‐complexity models to produce the predictions at a practical level for managers concerned about the persistence of bank erosion while acknowledging the great uncertainty in any landscape prediction. First, a model of confluence angles reproduced observed angles of a major confluence, but showed limited susceptibility to a major rearrangement of the channel planform downstream. Second, a probabilistic map of channel locations was created with a two‐parameter channel avulsion model. The predicted channel belt location was concentrated in the same area as the current channel belt. Finally, a suite of valley‐scale channel and braid plain characteristics were extracted from a light detection and ranging (LiDAR)‐derived surface. The characteristics demonstrated large‐scale stabilizing topographic influences on channel planform. The combination of independent analyses increased confidence in the conclusion that the Toklat River braided planform is a dynamically stable system due to large and persistent valley‐scale influences, and that a range of avulsive perturbations are likely to result in a relatively unchanged planform configuration in the short term.     

Development of Regional Curves for Hydrologic Landscape Regions (HLR) in the Contiguous United States

Regional curves relate drainage area to the bankfull channel characteristics discharge, cross‐sectional area, width, and mean depth. These curves are used for a variety of purposes, including aiding in the field identification of bankfull elevation and in the natural channel design process. When developing regional curves, the degree to which landform, geology, climate, and vegetation influence stream systems within a single physiographic province may not be fully considered. This study examined the use of the U.S. Geological Survey's Hydrologic Landscape Regions (HLR), as well as data from 2,856 independent sites throughout the contiguous United States (U.S.), to develop a set of regional curves (bankfull discharge, cross‐sectional area, width, and mean depth) for (1) the contiguous U.S., (2) each of the 20 HLRs, (3) each of the eight physiographic divisions, (4) 22 of the 25 physiographic provinces, and (5) individual HLRs within the physiographic provinces. These regional curves were then compared to each other, as well as those from the literature. Regional curves developed for individual HLRs, physiographic divisions, and physiographic provinces tended to outperform the contiguous U.S. indicating increased stratification was beneficial. Further stratifying physiographic provinces by HLR markedly improved regional curve reliability. Use of HLR as a basis of regional curve development, rather than physiographic region alone, may allow for the development of more robust regional curves.     

Continental‐Scale River Flow Modeling of the Mississippi River Basin Using High‐Resolution NHDPlus Dataset

As a key component of the National Flood Interoperability Experiment (NFIE), this article presents the continental scale river flow modeling of the Mississippi River Basin (MRB), using high‐resolution river data from NHDPlus. The Routing Application for Parallel computatIon of Discharge (RAPID) was applied to the MRB with more than 1.2 million river reaches for a 10‐year study (2005‐2014). Runoff data from the Variable Infiltration Capacity (VIC) model was used as input to RAPID. This article investigates the effect of topography on RAPID performance, the differences between the VIC‐RAPID streamflow simulations in the HUC‐2 regions of the MRB, and the impact of major dams on the streamflow simulations. The model performance improved when initial parameter values, especially the Muskingum K parameter, were estimated by taking topography into account. The statistical summary indicates the RAPID model performs better in the Ohio and Tennessee Regions and the Upper and Lower Mississippi River Regions in comparison to the western part of the MRB, due to the better performance of the VIC model. The model accuracy also increases when lakes and reservoirs are considered in the modeling framework. In general, results show the VIC‐RAPID streamflow simulation is satisfactory at the continental scale of the MRB.