Aquatic carbon cycling in the conterminous United States and implications for terrestrial carbon accounting.

Inland water ecosystems dynamically process, transport, and sequester carbon. However, the transport of carbon through aquatic environments has not been quantitatively integrated in the context of terrestrial ecosystems. Here, we present the first integrated assessment, to our knowledge, of freshwater carbon fluxes for the conterminous United States, where 106 (range: 71-149) teragrams of carbon per year (TgC⋅y(-1)) is exported downstream or emitted to the atmosphere and sedimentation stores 21 (range: 9-65) TgC⋅y(-1) in lakes and reservoirs. We show that there is significant regional variation in aquatic carbon flux, but verify that emission across stream and river surfaces represents the dominant flux at 69 (range: 36-110) TgC⋅y(-1) or 65% of the total aquatic carbon flux for the conterminous United States. Comparing our results with the output of a suite of terrestrial biosphere models (TBMs), we suggest that within the current modeling framework, calculations of net ecosystem production (NEP) defined as terrestrial only may be overestimated by as much as 27%. However, the internal production and mineralization of carbon in freshwaters remain to be quantified and would reduce the effect of including aquatic carbon fluxes within calculations of terrestrial NEP. Reconciliation of carbon mass-flux interactions between terrestrial and aquatic carbon sources and sinks will require significant additional research and modeling capacity.

Inland waters and their role in the carbon cycle of Alaska.

The magnitude of Alaska (AK) inland waters carbon (C) fluxes is likely to change in the future due to amplified climate warming impacts on the hydrology and biogeochemical processes in high latitude regions. Although current estimates of major aquatic C fluxes represent an essential baseline against which future change can be compared, a comprehensive assessment for AK has not yet been completed. To address this gap, we combined available data sets and applied consistent methodologies to estimate river lateral C export to the coast, river and lake carbon dioxide (CO2 ) and methane (CH4 ) emissions, and C burial in lakes for the six major hydrologic regions in the state. Estimated total aquatic C flux for AK was 41 Tg C/yr. Major components of this total flux, in Tg C/yr, were 18 for river lateral export, 17 for river CO2 emissions, and 8 for lake CO2 emissions. Lake C burial offset these fluxes by 2 Tg C/yr. River and lake CH4 emissions were 0.03 and 0.10 Tg C/yr, respectively. The Southeast and South central regions had the highest temperature, precipitation, terrestrial net primary productivity (NPP), and C yields (fluxes normalized to land area) were 77 and 42 g C·m-2 ·yr-1 , respectively. Lake CO2 emissions represented over half of the total aquatic flux from the Southwest (37 g C·m-2 ·yr-1 ). The North Slope, Northwest, and Yukon regions had lesser yields (11, 15, and 17 g C·m2 ·yr-1 ), but these estimates may be the most vulnerable to future climate change, because of the heightened sensitivity of arctic and boreal ecosystems to intensified warming. Total aquatic C yield for AK was 27 g C·m-2 ·yr-1 , which represented 16% of the estimated terrestrial NPP. Freshwater ecosystems represent a significant conduit for C loss, and a more comprehensive view of land-water-atmosphere interactions is necessary to predict future climate change impacts on the Alaskan ecosystem C balance.

Effects of Backpacker Use, Pack Stock Trail Use, and Pack Stock Grazing on Water-Quality Indicators, Including Nutrients, E. coli, Hormones, and Pharmaceuticals, in Yosemite National Park, USA.

We investigated how visitor-use affects water quality in wilderness in Yosemite National Park. During the summers of 2012-2014, we collected and analyzed surface-water samples for water-quality indicators, including fecal indicator bacteria Escherichia coli, nutrients (nitrogen, phosphorus, carbon), suspended sediment concentration, pharmaceuticals, and hormones. Samples were collected upstream and downstream from different types of visitor use at weekly to biweekly intervals and during summer storms. We conducted a park-wide synoptic sampling campaign during summer 2014, and sampled upstream and downstream from meadows to evaluate the mitigating effect of meadows on water quality. At pack stock stream crossings, Escherichia coli concentrations were greater downstream from crossings than upstream (median downstream increase in Escherichia coli of three colony forming units 100 mL-1), with the greatest increases occurring during storms (median downstream increase in Escherichia coli of 32 CFU 100 mL-1). At backpacker use sites, hormones, and pharmaceuticals (e.g., insect repellent) were detected at downstream sites, and Escherichia coli concentrations were greater at downstream sites (median downstream increase in Escherichia coli of 1 CFU 100 mL-1). Differences in water quality downstream vs. upstream from meadows grazed by pack stock were not detectable for most water-quality indicators, however, Escherichia coli concentrations decreased downstream, suggesting entrapment and die-off of fecal indicator bacteria in meadows. Our results indicate that under current-use levels pack stock trail use and backpacker use are associated with detectable, but relatively minor, effects on water quality, which are most pronounced during storms.

Biogeochemistry of beetle-killed forests: explaining a weak nitrate response.

A current pine beetle infestation has caused extensive mortality of lodgepole pine (Pinus contorta) in forests of Colorado and Wyoming; it is part of an unprecedented multispecies beetle outbreak extending from Mexico to Canada. In United States and European watersheds, where atmospheric deposition of inorganic N is moderate to low (<10 kg⋅ha⋅y), disturbance of forests by timber harvest or violent storms causes an increase in stream nitrate concentration that typically is close to 400% of predisturbance concentrations. In contrast, no significant increase in streamwater nitrate concentrations has occurred following extensive tree mortality caused by the mountain pine beetle in Colorado. A model of nitrate release from Colorado watersheds calibrated with field data indicates that stimulation of nitrate uptake by vegetation components unaffected by beetles accounts for significant nitrate retention in beetle-infested watersheds. The combination of low atmospheric N deposition (<10 kg⋅ha⋅y), tree mortality spread over multiple years, and high compensatory capacity associated with undisturbed residual vegetation and soils explains the ability of these beetle-infested watersheds to retain nitrate despite catastrophic mortality of the dominant canopy tree species.

Organic Carbon Burial in Lakes and Reservoirs of the Conterminous United States.

Organic carbon (OC) burial in lacustrine sediments represents an important sink in the global carbon cycle; however, large-scale OC burial rates are poorly constrained, primarily because of the sparseness of available data sets. Here we present an analysis of OC burial rates in water bodies of the conterminous U.S. (CONUS) that takes advantage of recently developed national-scale data sets on reservoir sedimentation rates, sediment OC concentrations, lake OC burial rates, and water body distributions. We relate these data to basin characteristics and land use in a geostatistical analysis to develop an empirical model of OC burial in water bodies of the CONUS. Our results indicate that CONUS water bodies sequester 20.8 (95% CI: 9.4-65.8) Tg C yr(-1), and spatial patterns in OC burial are strongly influenced by water body type, size, and abundance; land use; and soil and vegetation characteristics in surrounding areas. Carbon burial is greatest in the central and southeastern regions of the CONUS, where cultivation and an abundance of small water bodies enhance accumulation of sediment and OC in aquatic environments.

Links between N deposition and nitrate export from a high-elevation watershed in the Colorado Front Range.

Long-term patterns of stream nitrate export and atmospheric N deposition were evaluated over three decades in Loch Vale, a high-elevation watershed in the Colorado Front Range. Stream nitrate concentrations increased in the early 1990 s, peaked in the mid-2000s, and have since declined by over 40%, coincident with trends in nitrogen oxide emissions over the past decade. Similarities in the timing and magnitude of N deposition provide evidence that stream chemistry is responding to changes in atmospheric deposition. The response to deposition was complicated by a drought in the early 2000s that enhanced N export for several years. Other possible explanations, including forest disturbance, snow depth, or permafrost melting, could not explain patterns in N export. Our results show that stream chemistry responds rapidly to changes in N deposition in high-elevation watersheds, similar to the response observed to changes in sulfur deposition.

Effects of stock use and backpackers on water quality in wilderness in Sequoia and Kings Canyon National Parks, USA.

During 2010-2011, a study was conducted in Sequoia and Kings Canyon National Parks (SEKI) to evaluate the influence of pack animals (stock) and backpackers on water quality in wilderness lakes and streams. The study had three main components: (1) a synoptic survey of water quality in wilderness areas of the parks, (2) paired water quality sampling above and below several areas with differing types and amounts of visitor use, and (3) intensive monitoring at six sites to document temporal variations in water quality. Data from the synoptic water quality survey indicated that wilderness lakes and streams are dilute and have low nutrient and Escherichia coli concentrations. The synoptic survey sites were categorized as minimal use, backpacker-use, or mixed use (stock and backpackers), depending on the most prevalent type of use upstream from the sampling locations. Sites with mixed use tended to have higher concentrations of most constituents (including E. coli) than those categorized as minimal-use (P ≤ 0.05); concentrations at backpacker-use sites were intermediate. Data from paired-site sampling indicated that E. coli, total coliform, and particulate phosphorus concentrations were greater in streams downstream from mixed-use areas than upstream from those areas (P ≤ 0.05). Paired-site data also indicated few statistically significant differences in nutrient, E. coli, or total coliform concentrations in streams upstream and downstream from backpacker-use areas. The intensive-monitoring data indicated that nutrient and E. coli concentrations normally were low, except during storms, when notable increases in concentrations of E. coli, nutrients, dissolved organic carbon, and turbidity occurred. In summary, results from this study indicate that water quality in SEKI wilderness generally is good, except during storms; and visitor use appears to have a small, but statistically significant influence on stream water quality.

Assessing possible visitor-use impacts on water quality in Yosemite National Park, California.

There is concern that visitor-use associated activities, such as bathing, dish washing, wastewater production, and stock animal use near lakes and streams, could cause degradation of water quality in Yosemite National Park. A study was conducted during 2004-2007 to assess patterns in nutrient and Escherichia coli (E. coli) concentrations in the Merced and Tuolumne Rivers and characterize natural background concentrations of nutrients in the park. Results indicated that nutrient and E. coli concentrations were low, even compared to other undeveloped sites in the United States. A multiple linear regression approach was used to model natural background concentrations of nutrients, with basin characteristics as explanatory variables. Modeled nitrogen concentrations increased with elevation, and modeled phosphorus concentrations increased with basin size. Observed concentrations (±uncertainty) were compared to modeled concentrations (±uncertainty) to identify sites that might be impacted by point sources of nutrients, as indicated by large model residuals. Statistically significant differences in observed and modeled concentrations were observed at only a few locations, indicating that most sites were representative of natural background conditions. The empirical modeling approach used in this study can be used to estimate natural background conditions at any point along a study reach in areas minimally impacted by development, and may be useful for setting water-quality standards in many national parks.

Preferential elution of ionic solutes in melting snowpacks: Improving process understanding through field observations and modeling in the Rocky Mountains.

The preferential elution of ions from melting snowpacks is a complex problem that has been linked to temporary acidification of water bodies. However, the understanding of these processes in snowpacks around the world, including the polar regions that are experiencing unprecedented warming and melting, remains limited despite being instrumental in supporting climate change adaptation. In this study, data collected from a snowmelt lysimeter and snowpits at meadow and forest-gap sites in a high elevation watershed in Colorado were combined with the PULSE multi-phase snowpack chemistry model to investigate the controls of meltwater chemistry and preferential elution. The snowdepth at the meadow site was 64% of that at the forest-gap site, and the snowmelt rate was greater there (meadow snowpit) due to higher solar irradiance. Cations such as Ca2+ and NH4+ were deposited mostly within the upper layers of both the meadow and forest-gap snowpacks, and acid anions such as NO3- and SO42- were more evenly distributed. The snow ion concentrations were generally greater at the forest-gap snowpit, except for NH4+, which indicates that wind erosion of wet and dry deposited ions from the meadow may have reduced concentrations of residual snow. Furthermore, at the forest-gap site, snow interception and scavenging processes such as sublimation, ventilation, and throughfall led to particular ion enrichment of Ca2+, Mg2+, K+, Cl-, SO42- and NO3-. Model simulations and observations highlight that preferential elution is enhanced by low snowmelt rates, with the model indicating that this is due to lower dilution rates and increased contact time and area between the percolating meltwater and the snow. Results suggest that low snowmelt rates can cause multiple early meltwater ionic pulses for ions subject to lower ion exclusion. Ion exclusion rates at the grain-size level have been estimated for the first time.

Comparison of total mercury and methylmercury cycling at five sites using the small watershed approach.

The small watershed approach is well-suited but underutilized in mercury research. We applied the small watershed approach to investigate total mercury (THg) and methylmercury (MeHg) dynamics in streamwater at the five diverse forested headwater catchments of the US Geological Survey Water, Energy, and Biogeochemical Budgets (WEBB) program. At all sites, baseflow THg was generally less than 1ng L(-1) and MeHg was less than 0.2ng L(-1). THg and MeHg concentrations increased with streamflow, so export was primarily episodic. At three sites, THg and MeHg concentration and export were dominated by the particulate fraction in association with POC at high flows, with maximum THg (MeHg) concentrations of 94 (2.56)ng L(-1) at Sleepers River, Vermont; 112 (0.75)ng L(-1) at Rio Icacos, Puerto Rico; and 55 (0.80)ng L(-1) at Panola Mt., Georgia. Filtered (<0.7microm) THg increased more modestly with flow in association with the hydrophobic acid fraction (HPOA) of DOC, with maximum filtered THg concentrations near 5ng L(-1) at both Sleepers and Icacos. At Andrews Creek, Colorado, THg export was also episodic but was dominated by filtered THg, as POC concentrations were low. MeHg typically tracked THg so that each site had a fairly constant MeHg/THg ratio, which ranged from near zero at Andrews to 15% at the low-relief, groundwater-dominated Allequash Creek, Wisconsin. Allequash was the only site with filtered MeHg consistently above detection, and the filtered fraction dominated both THg and MeHg. Relative to inputs in wet deposition, watershed retention of THg (minus any subsequent volatilization) was 96.6% at Allequash, 60% at Sleepers, and 83% at Andrews. Icacos had a net export of THg, possibly due to historic gold mining or frequent disturbance from landslides. Quantification and interpretation of Hg dynamics was facilitated by the small watershed approach with emphasis on event sampling.

Organic carbon burial in global lakes and reservoirs.

Burial in sediments removes organic carbon (OC) from the short-term biosphere-atmosphere carbon (C) cycle, and therefore prevents greenhouse gas production in natural systems. Although OC burial in lakes and reservoirs is faster than in the ocean, the magnitude of inland water OC burial is not well constrained. Here we generate the first global-scale and regionally resolved estimate of modern OC burial in lakes and reservoirs, deriving from a comprehensive compilation of literature data. We coupled statistical models to inland water area inventories to estimate a yearly OC burial of 0.15 (range, 0.06-0.25) Pg C, of which ~40% is stored in reservoirs. Relatively higher OC burial rates are predicted for warm and dry regions. While we report lower burial than previously estimated, lake and reservoir OC burial corresponded to ~20% of their C emissions, making them an important C sink that is likely to increase with eutrophication and river damming.

Mapping critical loads of nitrogen deposition for aquatic ecosystems in the Rocky Mountains, USA.

Spatially explicit estimates of critical loads of nitrogen (N) deposition (CL(Ndep)) for nutrient enrichment in aquatic ecosystems were developed for the Rocky Mountains, USA, using a geostatistical approach. The lowest CL(Ndep) estimates (<1.5 ± 1 kg N ha(-1) yr(-1)) occurred in high-elevation basins with steep slopes, sparse vegetation, and abundance of exposed bedrock and talus. These areas often correspond with areas of high N deposition (>3 kg N ha(-1) yr(-1)), resulting in CL(Ndep) exceedances ≥ 1.5 ± 1 kg N ha(-1) yr(-1). CL(Ndep) and CL(Ndep) exceedances exhibit substantial spatial variability related to basin characteristics and are highly sensitive to the NO(3)(-) threshold at which ecological effects are thought to occur. Based on an NO(3)(-) threshold of 0.5 µmol L(-1), N deposition exceeds CL(Ndep) in 21 ± 8% of the study area; thus, broad areas of the Rocky Mountains may be impacted by excess N deposition, with greatest impacts at high elevations.

Melting Alpine glaciers enrich high-elevation lakes with reactive nitrogen.

Alpine glaciers have receded substantially over the last century in many regions of the world. Resulting changes in glacial runoff not only affect the hydrological cycle, but can also alter the physical (i.e., turbidity from glacial flour) and biogeochemical properties of downstream ecosystems. Here we compare nutrient concentrations, transparency gradients, algal biomass, and fossil diatom species richness in two sets of high-elevation lakes: those fed by snowpack melt alone (SF lakes) and those fed by both glacial and snowpack meltwaters (GSF lakes). We found that nitrate (NO(3)(-)) concentrations in the GSF lakes were 1-2 orders of magnitude higher than in SF lakes. Although nitrogen (N) limitation is common in alpine lakes, algal biomass was lower in highly N-enriched GSF lakes than in the N-poor SF lakes. Contrary to expectations, GSF lakes were more transparent than SF lakes to ultraviolet and equally transparent to photosynthetically active radiation. Sediment diatom assemblages had lower taxonomic richness in the GSF lakes, a feature that has persisted over the last century. Our results demonstrate that the presence of glaciers on alpine watersheds more strongly influences NO(3)(-)concentrations in high-elevation lake ecosystems than any other geomorphic or biogeographic characteristic.