Global methane emissions from rivers and streams.

Methane (CH4) is a potent greenhouse gas and its concentrations have tripled in the atmosphere since the industrial revolution. There is evidence that global warming has increased CH4 emissions from freshwater ecosystems1,2, providing positive feedback to the global climate. Yet for rivers and streams, the controls and the magnitude of CH4 emissions remain highly uncertain3,4. Here we report a spatially explicit global estimate of CH4 emissions from running waters, accounting for 27.9 (16.7-39.7) Tg CH4 per year and roughly equal in magnitude to those of other freshwater systems5,6. Riverine CH4 emissions are not strongly temperature dependent, with low average activation energy (EM = 0.14 eV) compared with that of lakes and wetlands (EM = 0.96 eV)1. By contrast, global patterns of emissions are characterized by large fluxes in high- and low-latitude settings as well as in human-dominated environments. These patterns are explained by edaphic and climate features that are linked to anoxia in and near fluvial habitats, including a high supply of organic matter and water saturation in hydrologically connected soils. Our results highlight the importance of land-water connections in regulating CH4 supply to running waters, which is vulnerable not only to direct human modifications but also to several climate change responses on land.

Long-range dependence and extreme values of precipitation, phosphorus load, and Cyanobacteria.

Extreme daily values of precipitation (1939-2021), discharge (1991-2021), phosphorus (P) load (1994-2021), and phycocyanin, a pigment of Cyanobacteria (June 1-September 15 of 2008-2021) are clustered as multi-day events for Lake Mendota, Wisconsin. Long-range dependence, or memory, is the shortest for precipitation and the longest for phycocyanin. Extremes are clustered for all variates and those of P load and phycocyanin are most strongly clustered. Extremes of P load are predictable from extremes of precipitation, and precipitation and P load are correlated with later concentrations of phycocyanin. However, time delays from 1 to 60 d were found between P load extremes and the next extreme phycocyanin event within the same year of observation. Although most of the lake's P enters in extreme events, blooms of Cyanobacteria may be sustained by recycling and food web processes.

Methane concentrations and fluxes in agricultural and preserved tropical headwater streams.

Tropical streams have been intensively impacted by agricultural activities. Among the most important agricultural activities in Brazil, sugarcane production represents a large impact for economic development and for environmental conditions. Permeating sugarcane fields, several headwater streams can be affected by sugarcane cultivation, in special, aquatic biogeochemical cycles because of the deforestation, fertilization, crop residues and higher temperatures in the tropics. In this study, we analyzed the effects of sugarcane cultivation on methane fluxes and concentrations, assuming that carbon cycles are influenced by agricultural activities in headwater streams. Our study aimed to (1) measure methane fluxes and concentrations in tropical streams located in Southeastern Brazil, (2) Analyze whether seasonal cycles influence methane fluxes and concentrations, (3) Evaluate the influence of sugarcane cultivation on methane fluxes and (4) Analyze the association between water chemistry in the methane concentrations in tropical streams. We found mean fluxes of CH4 of 0.280 mmol m-2 d-1, with higher fluxes during the summer and in streams draining preserved catchments. The average CH4 concentrations were 0.695 µmol L-1, with higher values during the summer and in streams draining preserved catchments. Methane concentrations in the studied streams was influenced by dissolved oxygen (negatively), dissolved organic carbon (negatively), water velocity (positively) and conductivity (negatively). Methane concentrations were significantly higher than concentrations found in Temperate Grasslands, Savannas & Shrublands and similar to concentrations found in other tropical biomes (excluding Tropical & Subtropical Moist Broadleaf Forests which receives large amounts of organic inputs). We conclude that sugarcane influence methane concentrations and fluxes in tropical streams by reducing the organic matter availability provided by the native vegetation in soil and water.

Unexpectedly minor nitrous oxide emissions from fluvial networks draining permafrost catchments of the East Qinghai-Tibet Plateau.

Streams and rivers emit substantial amounts of nitrous oxide (N2O) and are therefore an essential component of global nitrogen (N) cycle. Permafrost soils store a large reservoir of dormant N that, upon thawing, can enter fluvial networks and partly degrade to N2O, yet the role of waterborne release of N2O in permafrost regions is unclear. Here we report N2O concentrations and fluxes during different seasons between 2016 and 2018 in four watersheds on the East Qinghai-Tibet Plateau. Thawing permafrost soils are known to emit N2O at a high rate, but permafrost rivers draining the East Qinghai-Tibet Plateau behave as unexpectedly minor sources of atmospheric N2O. Such low N2O fluxes are associated with low riverine dissolved inorganic N (DIN) after terrestrial plant uptake, unfavorable conditions for N2O generation via denitrification, and low N2O yield due to a small ratio of nitrite reductase: nitrous oxide reductase in these rivers. We estimate fluvial N2O emissions of 0.432 - 0.463 Gg N2O-N yr-1 from permafrost landscapes on the entire Qinghai-Tibet Plateau, which is marginal (~0.15%) given their areal contribution to global streams and rivers (0.7%). However, we suggest that these permafrost-affected rivers can shift from minor sources to strong emitters in the warmer future, likely giving rise to the permafrost non-carbon feedback that intensifies warming.

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.

Outsized nutrient contributions from small tributaries to a Great Lake.

Excessive nitrogen (N) and phosphorus (P) loading is one of the greatest threats to aquatic ecosystems in the Anthropocene, causing eutrophication of rivers, lakes, and marine coastlines worldwide. For lakes across the United States, eutrophication is driven largely by nonpoint nutrient sources from tributaries that drain surrounding watersheds. Decades of monitoring and regulatory efforts have paid little attention to small tributaries of large water bodies, despite their ubiquity and potential local importance. We used a snapshot of nutrient inputs from nearly all tributaries of Lake Michigan-the world's fifth largest freshwater lake by volume-to determine how land cover and dams alter nutrient inputs across watershed sizes. Loads, concentrations, stoichiometry (N:P), and bioavailability (percentage dissolved inorganic nutrients) varied by orders of magnitude among tributaries, creating a mosaic of coastal nutrient inputs. The 6 largest of 235 tributaries accounted for ∼70% of the daily N and P delivered to Lake Michigan. However, small tributaries exhibited nutrient loads that were high for their size and biased toward dissolved inorganic forms. Higher bioavailability of nutrients from small watersheds suggests greater potential to fuel algal blooms in coastal areas, especially given the likelihood that their plumes become trapped and then overlap in the nearshore zone. Our findings reveal an underappreciated role that small streams may play in driving coastal eutrophication in large water bodies. Although they represent only a modest proportion of lake-wide loads, expanding nutrient management efforts to address smaller watersheds could reduce the ecological impacts of nutrient loading on valuable nearshore ecosystems.

Spatial patterns of enzymatic activity in large water bodies: Ship-borne measurements of beta-D-glucuronidase activity as a rapid indicator of microbial water quality.

This study used automated enzymatic activity measurements conducted from a mobile research vessel to detect the spatial variability of beta­d­glucuronidase (GLUC) activity in large freshwater bodies. The ship-borne observations provided the first high-resolution spatial data of GLUC activity in large water bodies as rapid indication of fecal pollution and were used to identify associations with hydrological conditions and land use. The utility of this novel approach for water quality screening was evaluated by surveys of the Columbia River, the Mississippi River and the Yahara Lakes, covering up to a 500 km river course and 50 km2 lake area. The ship-borne measurements of GLUC activity correlated with standard E. coli analyses (R2 = 0.71) and revealed the effects of (1) precipitation events and urban run-off on GLUC activity in surface waters, (2) localized point inlets of potential fecal pollution and (3) increasing GLUC signals along gradients of urbanization. We propose that this ship-borne water quality screening to be integrated into future water inventory programs as an initial or complementary tool (besides established fecal indicator parameters), due to its ability to provide near real-time spatial information on potential fecal contamination of large surface water resources and therefore being helpful to greatly reduce potential human health risks.

The metabolic regimes of 356 rivers in the United States.

A national-scale quantification of metabolic energy flow in streams and rivers can improve understanding of the temporal dynamics of in-stream activity, links between energy cycling and ecosystem services, and the effects of human activities on aquatic metabolism. The two dominant terms in aquatic metabolism, gross primary production (GPP) and aerobic respiration (ER), have recently become practical to estimate for many sites due to improved modeling approaches and the availability of requisite model inputs in public datasets. We assembled inputs from the U.S. Geological Survey and National Aeronautics and Space Administration for October 2007 to January 2017. We then ran models to estimate daily GPP, ER, and the gas exchange rate coefficient for 356 streams and rivers across the continental United States. We also gathered potential explanatory variables and spatial information for cross-referencing this dataset with other datasets of watershed characteristics. This dataset offers a first national assessment of many-day time series of metabolic rates for up to 9 years per site, with a total of 490,907 site-days of estimates.

Unexpected stasis in a changing world: Lake nutrient and chlorophyll trends since 1990.

The United States (U.S.) has faced major environmental changes in recent decades, including agricultural intensification and urban expansion, as well as changes in atmospheric deposition and climate-all of which may influence eutrophication of freshwaters. However, it is unclear whether or how water quality in lakes across diverse ecological settings has responded to environmental change. We quantified water quality trends in 2913 lakes using nutrient and chlorophyll (Chl) observations from the Lake Multi-Scaled Geospatial and Temporal Database of the Northeast U.S. (LAGOS-NE), a collection of preexisting lake data mostly from state agencies. LAGOS-NE was used to quantify whether lake water quality has changed from 1990 to 2013, and whether lake-specific or regional geophysical factors were related to the observed changes. We modeled change through time using hierarchical linear models for total nitrogen (TN), total phosphorus (TP), stoichiometry (TN:TP), and Chl. Both the slopes (percent change per year) and intercepts (value in 1990) were allowed to vary by lake and region. Across all lakes, TN declined at a rate of 1.1% year-1 , while TP, TN:TP, and Chl did not change. A minority (7%-16%) of individual lakes had changing nutrients, stoichiometry, or Chl. Of those lakes that changed, we found differences in the geospatial variables that were most related to the observed change in the response variables. For example, TN and TN:TP trends were related to region-level drivers associated with atmospheric deposition of N; TP trends were related to both lake and region-level drivers associated with climate and land use; and Chl trends were found in regions with high air temperature at the beginning of the study period. We conclude that despite large environmental change and management efforts over recent decades, water quality of lakes in the Midwest and Northeast U.S. has not overwhelmingly degraded or improved.

Lake nutrient stoichiometry is less predictable than nutrient concentrations at regional and sub-continental scales.

Production in many ecosystems is co-limited by multiple elements. While a known suite of drivers associated with nutrient sources, nutrient transport, and internal processing controls concentrations of phosphorus (P) and nitrogen (N) in lakes, much less is known about whether the drivers of single nutrient concentrations can also explain spatial or temporal variation in lake N:P stoichiometry. Predicting stoichiometry might be more complex than predicting concentrations of individual elements because some drivers have similar relationships with N and P, leading to a weak relationship with their ratio. Further, the dominant controls on elemental concentrations likely vary across regions, resulting in context dependent relationships between drivers, lake nutrients and their ratios. Here, we examine whether known drivers of N and P concentrations can explain variation in N:P stoichiometry, and whether explaining variation in stoichiometry differs across regions. We examined drivers of N:P in ~2,700 lakes at a sub-continental scale and two large regions nested within the sub-continental study area that have contrasting ecological context, including differences in the dominant type of land cover (agriculture vs. forest). At the sub-continental scale, lake nutrient concentrations were correlated with nutrient loading and lake internal processing, but stoichiometry was only weakly correlated to drivers of lake nutrients. At the regional scale, drivers that explained variation in nutrients and stoichiometry differed between regions. In the Midwestern U.S. region, dominated by agricultural land use, lake depth and the percentage of row crop agriculture were strong predictors of stoichiometry because only phosphorus was related to lake depth and only nitrogen was related to the percentage of row crop agriculture. In contrast, all drivers were related to N and P in similar ways in the Northeastern U.S. region, leading to weak relationships between drivers and stoichiometry. Our results suggest ecological context mediates controls on lake nutrients and stoichiometry. Predicting stoichiometry was generally more difficult than predicting nutrient concentrations, but human activity may decouple N and P, leading to better prediction of N:P stoichiometry in regions with high anthropogenic activity.

Ecology under lake ice.

Winter conditions are rapidly changing in temperate ecosystems, particularly for those that experience periods of snow and ice cover. Relatively little is known of winter ecology in these systems, due to a historical research focus on summer 'growing seasons'. We executed the first global quantitative synthesis on under-ice lake ecology, including 36 abiotic and biotic variables from 42 research groups and 101 lakes, examining seasonal differences and connections as well as how seasonal differences vary with geophysical factors. Plankton were more abundant under ice than expected; mean winter values were 43.2% of summer values for chlorophyll a, 15.8% of summer phytoplankton biovolume and 25.3% of summer zooplankton density. Dissolved nitrogen concentrations were typically higher during winter, and these differences were exaggerated in smaller lakes. Lake size also influenced winter-summer patterns for dissolved organic carbon (DOC), with higher winter DOC in smaller lakes. At coarse levels of taxonomic aggregation, phytoplankton and zooplankton community composition showed few systematic differences between seasons, although literature suggests that seasonal differences are frequently lake-specific, species-specific, or occur at the level of functional group. Within the subset of lakes that had longer time series, winter influenced the subsequent summer for some nutrient variables and zooplankton biomass.

Controls on methane concentrations and fluxes in streams draining human-dominated landscapes.

Streams and rivers are active processors of carbon, leading to significant emissions of CO2 and possibly CH4 to the atmosphere. Patterns and controls of CH4 in fluvial ecosystems remain relatively poorly understood. Furthermore, little is known regarding how major human impacts to fluvial ecosystems may be transforming their role as CH4 producers and emitters. Here, we examine the consequences of two distinct ecosystem changes as a result of human land use: increased nutrient loading (primarily as nitrate), and increased sediment loading and deposition of fine particles in the benthic zone. We did not find support for the hypothesis that enhanced nitrate loading down-regulates methane production via thermodynamic or toxic effects. We did find strong evidence that increased sedimentation and enhanced organic matter content of the benthos lead to greater methane production (diffusive + ebullitive flux) relative to pristine fluvial systems in northern Wisconsin (upper Midwest, USA). Overall, streams in a human-dominated landscape of southern Wisconsin were major regional sources of CH4 to the atmosphere, equivalent to ~20% of dairy cattle emissions, or ~50% of a landfill's annual emissions. We suggest that restoration of the benthic environment (reduced fine deposits) could lead to reduced CH4 emissions, while decreasing nutrient loading is likely to have limited impacts to this ecosystem process.

Building a multi-scaled geospatial temporal ecology database from disparate data sources: fostering open science and data reuse.

Although there are considerable site-based data for individual or groups of ecosystems, these datasets are widely scattered, have different data formats and conventions, and often have limited accessibility. At the broader scale, national datasets exist for a large number of geospatial features of land, water, and air that are needed to fully understand variation among these ecosystems. However, such datasets originate from different sources and have different spatial and temporal resolutions. By taking an open-science perspective and by combining site-based ecosystem datasets and national geospatial datasets, science gains the ability to ask important research questions related to grand environmental challenges that operate at broad scales. Documentation of such complicated database integration efforts, through peer-reviewed papers, is recommended to foster reproducibility and future use of the integrated database. Here, we describe the major steps, challenges, and considerations in building an integrated database of lake ecosystems, called LAGOS (LAke multi-scaled GeOSpatial and temporal database), that was developed at the sub-continental study extent of 17 US states (1,800,000 km(2)). LAGOS includes two modules: LAGOSGEO, with geospatial data on every lake with surface area larger than 4 ha in the study extent (~50,000 lakes), including climate, atmospheric deposition, land use/cover, hydrology, geology, and topography measured across a range of spatial and temporal extents; and LAGOSLIMNO, with lake water quality data compiled from ~100 individual datasets for a subset of lakes in the study extent (~10,000 lakes). Procedures for the integration of datasets included: creating a flexible database design; authoring and integrating metadata; documenting data provenance; quantifying spatial measures of geographic data; quality-controlling integrated and derived data; and extensively documenting the database. Our procedures make a large, complex, and integrated database reproducible and extensible, allowing users to ask new research questions with the existing database or through the addition of new data. The largest challenge of this task was the heterogeneity of the data, formats, and metadata. Many steps of data integration need manual input from experts in diverse fields, requiring close collaboration.

Ebullitive methane emissions from oxygenated wetland streams.

Stream and river carbon dioxide emissions are an important component of the global carbon cycle. Methane emissions from streams could also contribute to regional or global greenhouse gas cycling, but there are relatively few data regarding stream and river methane emissions. Furthermore, the available data do not typically include the ebullitive (bubble-mediated) pathway, instead focusing on emission of dissolved methane by diffusion or convection. Here, we show the importance of ebullitive methane emissions from small streams in the regional greenhouse gas balance of a lake and wetland-dominated landscape in temperate North America and identify the origin of the methane emitted from these well-oxygenated streams. Stream methane flux densities from this landscape tended to exceed those of nearby wetland diffusive fluxes as well as average global wetland ebullitive fluxes. Total stream ebullitive methane flux at the regional scale (103 Mg C yr(-1) ; over 6400 km(2) ) was of the same magnitude as diffusive methane flux previously documented at the same scale. Organic-rich stream sediments had the highest rates of bubble release and higher enrichment of methane in bubbles, but glacial sand sediments also exhibited high bubble emissions relative to other studied environments. Our results from a database of groundwater chemistry support the hypothesis that methane in bubbles is produced in anoxic near-stream sediment porewaters, and not in deeper, oxygenated groundwaters. Methane interacts with other key elemental cycles such as nitrogen, oxygen, and sulfur, which has implications for ecosystem changes such as drought and increased nutrient loading. Our results support the contention that streams, particularly those draining wetland landscapes of the northern hemisphere, are an important component of the global methane cycle.

Lake microbial communities are resilient after a whole-ecosystem disturbance.

Disturbances act as powerful structuring forces on ecosystems. To ask whether environmental microbial communities have capacity to recover after a large disturbance event, we conducted a whole-ecosystem manipulation, during which we imposed an intense disturbance on freshwater microbial communities by artificially mixing a temperate lake during peak summer thermal stratification. We employed environmental sensors and water chemistry analyses to evaluate the physical and chemical responses of the lake, and bar-coded 16S ribosomal RNA gene pyrosequencing and automated ribosomal intergenic spacer analysis (ARISA) to assess the bacterial community responses. The artificial mixing increased mean lake temperature from 14 to 20 °C for seven weeks after mixing ended, and exposed the microorganisms to very different environmental conditions, including increased hypolimnion oxygen and increased epilimnion carbon dioxide concentrations. Though overall ecosystem conditions remained altered (with hypolimnion temperatures elevated from 6 to 20 °C), bacterial communities returned to their pre-manipulation state as some environmental conditions, such as oxygen concentration, recovered. Recovery to pre-disturbance community composition and diversity was observed within 7 (epilimnion) and 11 (hypolimnion) days after mixing. Our results suggest that some microbial communities have capacity to recover after a major disturbance.

Fate of allochthonous dissolved organic carbon in lakes: a quantitative approach.

Inputs of dissolved organic carbon (DOC) to lakes derived from the surrounding landscape can be stored, mineralized or passed to downstream ecosystems. The balance among these OC fates depends on a suite of physical, chemical, and biological processes within the lake, as well as the degree of recalcintrance of the allochthonous DOC load. The relative importance of these processes has not been well quantified due to the complex nature of lakes, as well as challenges in scaling DOC degradation experiments under controlled conditions to the whole lake scale. We used a coupled hydrodynamic-water quality model to simulate broad ranges in lake area and DOC, two characteristics important to processing allochthonous carbon through their influences on lake temperature, mixing depth and hydrology. We calibrated the model to four lakes from the North Temperate Lakes Long Term Ecological Research site, and simulated an additional 12 'hypothetical' lakes to fill the gradients in lake size and DOC concentration. For each lake, we tested several mineralization rates (range: 0.001 d(-1) to 0.010 d(-1)) representative of the range found in the literature. We found that mineralization rates at the ecosystem scale were roughly half the values from laboratory experiments, due to relatively cool water temperatures and other lake-specific factors that influence water temperature and hydrologic residence time. Results from simulations indicated that the fate of allochthonous DOC was controlled primarily by the mineralization rate and the hydrologic residence time. Lakes with residence times <1 year exported approximately 60% of the DOC, whereas lakes with residence times >6 years mineralized approximately 60% of the DOC. DOC fate in lakes can be determined with a few relatively easily measured factors, such as lake morphometry, residence time, and temperature, assuming we know the recalcitrance of the DOC.

Changes in the dissolved nitrogen pool across land cover gradients in Wisconsin streams.

Increases in anthropogenic nitrogen fixation have resulted in wide-scale enrichment of aquatic ecosystems. Existing biogeochemical theory suggests that N enrichment is associated with increasing concentrations of nitrate; however, dissolved organic nitrogen (DON) is often a major component of the total dissolved nitrogen (TDN) pool in streams and rivers, and its concentration can be significantly elevated in human-influenced basins. We examined N concentrations during summer base flow conditions in 324 Wisconsin streams to determine whether DON was a significant component of TDN and how its relative contribution changed across a gradient of increasing human (agriculture and urban) land use for 84 of these sites. Total dissolved nitrogen varied from 0.09 to 20.74 mg/L, and although DON was significantly higher in human-dominated basins relative to forested and mixed-cover basins, its concentration increased relatively slowly in response to increasing human land cover. This limited response reflected a replacement of wetland-derived DON in low-N streams by anthropogenic sources in human-dominated sites, such that net changes in DON were small across the land use gradient. Nitrate-N increased exponentially in response to greater human land cover, and NH4-N and NO2-N were present at low levels. Nitrite-N exceeded NH4-N at 20% of sites and reached a maximum concentration of 0.10 mg/L. This examination suggests that basic mechanisms driving N losses from old-growth forests subject to N saturation also shape the summertime N pool in Wisconsin streams, in addition to other processes dictated by landscape context. The overwhelming role of human land use in determining the relative and absolute composition of the summertime N pool included (1) rapid increases in NO3-N, (2) limited changes in DON, and (3) the unexpected occurrence of NO2-N. High (>3 mg/L) TDN conditions dominated by NO3-N, regardless of landscape context or forms of N inputs, indicate a state of "N hypersaturation", which appears to be increasingly common in human-influenced streams and rivers. Many sites in agriculturally rich areas had NO2-N and NO3-N concentrations that, if sustained, are at chronically toxic levels for sensitive aquatic biota, suggesting that N enrichment now has local consequences for resident stream biota in addition to contributing to coastal eutrophication.