Light and flow regimes regulate the metabolism of rivers.

Mean annual temperature and mean annual precipitation drive much of the variation in productivity across Earth's terrestrial ecosystems but do not explain variation in gross primary productivity (GPP) or ecosystem respiration (ER) in flowing waters. We document substantial variation in the magnitude and seasonality of GPP and ER across 222 US rivers. In contrast to their terrestrial counterparts, most river ecosystems respire far more carbon than they fix and have less pronounced and consistent seasonality in their metabolic rates. We find that variation in annual solar energy inputs and stability of flows are the primary drivers of GPP and ER across rivers. A classification schema based on these drivers advances river science and informs management.

Shifting stoichiometry: Long-term trends in stream-dissolved organic matter reveal altered C:N ratios due to history of atmospheric acid deposition.

Dissolved organic carbon (DOC) and nitrogen (DON) are important energy and nutrient sources for aquatic ecosystems. In many northern temperate, freshwater systems DOC has increased in the past 50 years. Less is known about how changes in DOC may vary across latitudes, and whether changes in DON track those of DOC. Here, we present long-term DOC and DON data from 74 streams distributed across seven sites in biomes ranging from the tropics to northern boreal forests with varying histories of atmospheric acid deposition. For each stream, we examined the temporal trends of DOC and DON concentrations and DOC:DON molar ratios. While some sites displayed consistent positive or negative trends in stream DOC and DON concentrations, changes in direction or magnitude were inconsistent at regional or local scales. DON trends did not always track those of DOC, though DOC:DON ratios increased over time for ~30% of streams. Our results indicate that the dissolved organic matter (DOM) pool is experiencing fundamental changes due to the recovery from atmospheric acid deposition. Changes in DOC:DON stoichiometry point to a shifting energy-nutrient balance in many aquatic ecosystems. Sustained changes in the character of DOM can have major implications for stream metabolism, biogeochemical processes, food webs, and drinking water quality (including disinfection by-products). Understanding regional and global variation in DOC and DON concentrations is important for developing realistic models and watershed management protocols to effectively target mitigation efforts aimed at bringing DOM flux and nutrient enrichment under control.

Vegetation zones as indicators of denitrification potential in salt marshes.

Salt marsh vegetation zones shift in response to large-scale environmental changes such as sea-level rise (SLR) and restoration activities, but it is unclear if they are good indicators of soil nitrogen removal. Our goal was to characterize the relationship between denitrification potential and salt marsh vegetation zones in tidally restored and tidally unrestricted coastal marshes, and to use vegetation zones to extrapolate how SLR may influence high marsh denitrification at the landscape scale. We conducted denitrification enzyme activity assays on sediment collected from three vegetation zones expected to shift in distribution due to SLR and tidal flow restoration across 20 salt marshes in Connecticut, USA (n = 60 sampling plots) during the summer of 2017. We found lower denitrification potential in short-form Spartina alterniflora zones (mean, 95% CI: 4, 3-6 mg N h-1  m-2 ) than in S. patens (25, 15-36 mg N h-1  m-2 ) and Phragmites australis (56, 16-96 mg N h-1  m-2 ) zones. Vegetation zone was the single best predictor and explained 52% of the variation in denitrification potential; incorporating restoration status and soil characteristics (soil salinity, moisture, and ammonium) did not improve model fit. Because denitrification potential did not differ between tidally restored and unrestricted marshes, we suggest landscape-scale changes in denitrification after tidal restoration are likely to be associated with shifts in vegetation, rather than differences driven by restoration status. Sea-level-rise-induced hydrologic changes are widely observed to shift high marsh dominated by S. patens to short-form S. alterniflora. To explore the implications of this shift in dominant high marsh vegetation, we paired our measured mean denitrification potential rates with projections of high marsh loss from SLR. We found that, under low and medium SLR scenarios, predicted losses of denitrification potential due to replacement of S. patens by short-form S. alterniflora were substantially larger than losses due to reduced high marsh land area alone. Our results suggest that changes in vegetation zones can serve as landscape-scale predictors of the response of denitrification rates to rapid changes occurring in salt marshes.

Do waterbody classifications predict water quality?

Many states classify waterbodies according to groups of designated uses, which suggests that classifications may be correlated with water quality. The primary assessments of water quality in the United States (the Biennial Integrated Water Quality Reports) do not consider classification, so the relationship between classification and water quality is untested. Additionally, water quality has been shown to be influenced by watershed land use; however, land use is not typically part of waterbody classification systems. To determine the relationships between waterbody classification, water quality, watershed land cover, and forest fragmentation, we analyzed existing water quality data for the State of Connecticut from the United States Geological Survey and the Connecticut Department of Energy and Environmental Protection and land cover data from the National Land Cover Dataset. Connecticut uses a unique classification system that includes separation of drinking water sources (Class AA) and waterbodies receiving waste water discharges (Class B). Using a comparison of multiple means, we found that Class B waters had higher levels of nitrogen, solids, chloride, sodium, dissolved copper, total iron, and dissolved manganese than Class AA waters. Watersheds upstream of Class B segments had less forest cover, more development and more impervious cover than watersheds upstream of Class AA segments. Class A sites had some similarities in water quality and land cover with Class AA sites and some with Class B sites. The subset of Class B waterbodies with "Class AA-like" water quality also had "Class AA-like" land cover. Based on this and a multiple linear regression analysis, we found that water quality is more closely related to watershed land cover and forest fragmentation than to waterbody classification. Our results suggest that watershed land cover likely is a better proxy for water quality than waterbody classification.

Thermodynamic constraints on the utility of ecological stoichiometry for explaining global biogeochemical patterns.

Carbon and nitrogen cycles are coupled through both stoichiometric requirements for microbial biomass and dissimilatory metabolic processes in which microbes catalyse reduction-oxidation reactions. Here, we integrate stoichiometric theory and thermodynamic principles to explain the commonly observed trade-off between high nitrate and high organic carbon concentrations, and the even stronger trade-off between high nitrate and high ammonium concentrations, across a wide range of aquatic ecosystems. Our results suggest these relationships are the emergent properties of both microbial biomass stoichiometry and the availability of terminal electron acceptors. Because elements with multiple oxidation states (i.e. nitrogen, manganese, iron and sulphur) serve as both nutrients and sources of chemical energy in reduced environments, both assimilative demand and dissimilatory uses determine their concentrations across broad spatial gradients. Conceptual and quantitative models that integrate rather than independently examine thermodynamic, stoichiometric and evolutionary controls on biogeochemical cycling are essential for understanding local to global biogeochemical patterns.

Stream denitrification across biomes and its response to anthropogenic nitrate loading.

Anthropogenic addition of bioavailable nitrogen to the biosphere is increasing and terrestrial ecosystems are becoming increasingly nitrogen-saturated, causing more bioavailable nitrogen to enter groundwater and surface waters. Large-scale nitrogen budgets show that an average of about 20-25 per cent of the nitrogen added to the biosphere is exported from rivers to the ocean or inland basins, indicating that substantial sinks for nitrogen must exist in the landscape. Streams and rivers may themselves be important sinks for bioavailable nitrogen owing to their hydrological connections with terrestrial systems, high rates of biological activity, and streambed sediment environments that favour microbial denitrification. Here we present data from nitrogen stable isotope tracer experiments across 72 streams and 8 regions representing several biomes. We show that total biotic uptake and denitrification of nitrate increase with stream nitrate concentration, but that the efficiency of biotic uptake and denitrification declines as concentration increases, reducing the proportion of in-stream nitrate that is removed from transport. Our data suggest that the total uptake of nitrate is related to ecosystem photosynthesis and that denitrification is related to ecosystem respiration. In addition, we use a stream network model to demonstrate that excess nitrate in streams elicits a disproportionate increase in the fraction of nitrate that is exported to receiving waters and reduces the relative role of small versus large streams as nitrate sinks.

Nitrous oxide emission from denitrification in stream and river networks.

Nitrous oxide (N(2)O) is a potent greenhouse gas that contributes to climate change and stratospheric ozone destruction. Anthropogenic nitrogen (N) loading to river networks is a potentially important source of N(2)O via microbial denitrification that converts N to N(2)O and dinitrogen (N(2)). The fraction of denitrified N that escapes as N(2)O rather than N(2) (i.e., the N(2)O yield) is an important determinant of how much N(2)O is produced by river networks, but little is known about the N(2)O yield in flowing waters. Here, we present the results of whole-stream (15)N-tracer additions conducted in 72 headwater streams draining multiple land-use types across the United States. We found that stream denitrification produces N(2)O at rates that increase with stream water nitrate (NO(3)(-)) concentrations, but that <1% of denitrified N is converted to N(2)O. Unlike some previous studies, we found no relationship between the N(2)O yield and stream water NO(3)(-). We suggest that increased stream NO(3)(-) loading stimulates denitrification and concomitant N(2)O production, but does not increase the N(2)O yield. In our study, most streams were sources of N(2)O to the atmosphere and the highest emission rates were observed in streams draining urban basins. Using a global river network model, we estimate that microbial N transformations (e.g., denitrification and nitrification) convert at least 0.68 Tg·y(-1) of anthropogenic N inputs to N(2)O in river networks, equivalent to 10% of the global anthropogenic N(2)O emission rate. This estimate of stream and river N(2)O emissions is three times greater than estimated by the Intergovernmental Panel on Climate Change.

Temperature dependence of leaf breakdown in streams differs between organismal groups and leaf species.

Increased temperatures are altering rates of organic matter (OM) breakdown in stream ecosystems with implications for carbon (C) cycling in the face of global change. The metabolic theory of ecology (MTE) provides a framework for predicting temperature effects on OM breakdown, but differences in the temperature dependence of breakdown driven by different organismal groups (i.e., microorganisms vs. invertebrate detritivores) and litter species remain unresolved. Over two years, we conducted 12 60-day leaf litterbag incubations in 20 headwater streams in the southern Appalachian Mountains (USA). We compared temperature dependence (as activation energy, Ea) between microbial and detritivore-mediated breakdown, and between a highly recalcitrant (Rhododendron maximum) and a relatively labile (Acer rubrum) leaf species. Detritivore-mediated breakdown had a higher Ea than microbial breakdown for both leaf species (Rhododendron: 1.48 > 0.56 eV; Acer: 0.97 > 0.29 eV), and Rhododendron breakdown had a higher Ea than Acer breakdown for both organismal groups. Similarly, the Ea of total (coarse-mesh) Rhododendron breakdown was higher than the Ea of total Acer breakdown (0.89 > 0.52 eV). These effects for total breakdown were large, implying that the number of days to 95% mass loss would decline by 40% for Rhododendron and 26% for Acer between 12°C (our mean temperature value) and 16°C (+4°C, reflecting projected increases in global surface temperature due to climate change). Despite patterns in Ea, overall breakdown rates were higher for microbes than detritivores, and for Acer than Rhododendron over most of our temperature gradient. Additionally, the Ea for a subset of the microbial breakdown data declined from 0.40 to 0.22 eV when fungal biomass was included as a model predictor, highlighting the key role of fungi in determining the temperature dependence of litter breakdown. Our results imply that, as streams warm, routing of leaf litter C to detritivore-mediated fates will increase faster than predicted by previous studies and MTE, especially for labile litter. As temperatures rise, earlier depletion of autumn-shed, labile leaf litter combined with rapid breakdown rates of recalcitrant litter could exacerbate seasonal resource limitation and alter carbon storage and transport dynamics in temperate headwater stream networks.

How many mountains can we mine? Assessing the regional degradation of Central Appalachian rivers by surface coal mining.

Surface coal mining is the dominant form of land cover change in Central Appalachia, yet the extent to which surface coal mine runoff is polluting regional rivers is currently unknown. We mapped surface mining from 1976 to 2005 for a 19,581 km(2) area of southern West Virginia and linked these maps with water quality and biological data for 223 streams. The extent of surface mining within catchments is highly correlated with the ionic strength and sulfate concentrations of receiving streams. Generalized additive models were used to estimate the amount of watershed mining, stream ionic strength, or sulfate concentrations beyond which biological impairment (based on state biocriteria) is likely. We find this threshold is reached once surface coal mines occupy >5.4% of their contributing watershed area, ionic strength exceeds 308 µS cm(-1), or sulfate concentrations exceed 50 mg L(-1). Significant losses of many intolerant macroinvertebrate taxa occur when as little as 2.2% of contributing catchments are mined. As of 2005, 5% of the land area of southern WV was converted to surface mines, 6% of regional streams were buried in valley fills, and 22% of the regional stream network length drained watersheds with >5.4% of their surface area converted to mines.

Superlinear scaling of riverine biogeochemical function with watershed size.

River networks regulate carbon and nutrient exchange between continents, atmosphere, and oceans. However, contributions of riverine processing are poorly constrained at continental scales. Scaling relationships of cumulative biogeochemical function with watershed size (allometric scaling) provide an approach for quantifying the contributions of fluvial networks in the Earth system. Here we show that allometric scaling of cumulative riverine function with watershed area ranges from linear to superlinear, with scaling exponents constrained by network shape, hydrological conditions, and biogeochemical process rates. Allometric scaling is superlinear for processes that are largely independent of substrate concentration (e.g., gross primary production) due to superlinear scaling of river network surface area with watershed area. Allometric scaling for typically substrate-limited processes (e.g., denitrification) is linear in river networks with high biogeochemical activity or low river discharge but becomes increasingly superlinear under lower biogeochemical activity or high discharge, conditions that are widely prevalent in river networks. The frequent occurrence of superlinear scaling indicates that biogeochemical activity in large rivers contributes disproportionately to the function of river networks in the Earth system.

Continental-scale analysis of shallow and deep groundwater contributions to streams.

Groundwater discharge generates streamflow and influences stream thermal regimes. However, the water quality and thermal buffering capacity of groundwater depends on the aquifer source-depth. Here, we pair multi-year air and stream temperature signals to categorize 1729 sites across the continental United States as having major dam influence, shallow or deep groundwater signatures, or lack of pronounced groundwater (atmospheric) signatures. Approximately 40% of non-dam stream sites have substantial groundwater contributions as indicated by characteristic paired air and stream temperature signal metrics. Streams with shallow groundwater signatures account for half of all groundwater signature sites and show reduced baseflow and a higher proportion of warming trends compared to sites with deep groundwater signatures. These findings align with theory that shallow groundwater is more vulnerable to temperature increase and depletion. Streams with atmospheric signatures tend to drain watersheds with low slope and greater human disturbance, indicating reduced stream-groundwater connectivity in populated valley settings.