Geology Reports⌕ Search

USGS · 70003749

Comparisons of watershed sulfur budgets in southeast Canada and northeast US: New approaches and implications

Abstract

Most of eastern North America receives elevated levels of atmospheric deposition of sulfur (S) that result from anthropogenic SO 2 emissions from fossil fuel combustion. Atmospheric S deposition has acidified sensitive terrestrial and aquatic ecosystems in this region; however, deposition has been declining since the 1970s, resulting in some recovery in previously acidified aquatic ecosystems. Accurate watershed S mass balances help to evaluate the extent to which atmospheric S deposition is retained within ecosystems, and whether internal cycling sources and biogeochemical processes may be affecting the rate of recovery from decreasing S atmospheric loads. This study evaluated S mass balances for 15 sites with watersheds in southeastern Canada and northeastern US for the period 1985 to 2002. These 15 sites included nine in Canada (Turkey Lakes, ON; Harp Lake, ON; Plastic Lake, ON; Hermine, QC; Lake Laflamme, QC; Lake Clair, QC; Lake Tirasse, QC; Mersey, NS; Moosepit, NS) and six in the US (Arbutus Lake, NY; Biscuit Brook, NY; Sleepers River, VT; Hubbard Brook Experimental Forest, NH; Cone Pond, NH; Bear Brook Watershed, ME). Annual S wet deposition inputs were derived from measured bulk or wet-only deposition and stream export was obtained by combining drainage water fluxes with SO 4 2− concentrations. Dry deposition has the greatest uncertainty of any of the mass flux calculations necessary to develop accurate watershed balances, and here we developed a new method to calculate this quantity. We utilized historical information from both the US National Emissions Inventory and the US (CASTNET) and the Canadian (CAPMoN) dry deposition networks to develop a formulation that predicted SO 2 concentrations as a function of SO 2 emissions, latitude and longitude. The SO 2 concentrations were used to predict dry deposition using relationships between concentrations and deposition flux derived from the CASTNET or CAPMoN networks. For the year 2002, we compared the SO 2 concentrations and deposition predictions with the predictions of two continental-scale air quality models, the Community Multiscale Air Quality (CMAQ) model and A Unified Regional Air-quality Modeling System (AURAMS) that utilize complete inventories of emissions and chemical budgets. The results of this comparison indicated that the predictive relationship provides an accurate representation of SO 2 concentrations and S deposition for the region that is generally consistent with these models, and thus provides confidence that our approach could be used to develop accurate watershed S budgets for these 15 sites. Most watersheds showed large net losses of SO 4 2− on an annual basis, and the watershed mass balances were grouped into five categories based on the relative value of mean annual net losses or net gains. The net annual fluxes of SO 4 2− showed a strong relationship with hydrology; the largest net annual negative fluxes were associated with years of greatest precipitation amount and highest discharge. The important role of catchment hydrology on S budgets suggests implications for future predicted climate change as it affects patterns of precipitation and drought. The sensitivity of S budgets is likely to be greatest in watersheds with the greatest wetland area, which are particularly sensitive to drying and wetting cycles. A small number of the watersheds in this analysis were shown to have substantial S sources from mineral weathering, but most showed evidence of an internal source of SO 4 2− , which is likely from the mineralization of organic S stored from decades of increased S deposition. Mobilization of this internal S appears to contribute about 1–6 kg S ha −1 year −1 to stream fluxes at these sites and is affecting the rate and extent of recovery from acidification as S deposition rates have declined in recent years. This internal S source should be considered when developing critical deposition loads that will promote ecosystem recovery from acidification and the depletion of nutrient cations in the northeastern US and southeastern Canada.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 39.977120098439634° to 52.02545860348814° latitude; -87.1875° to -63.28125° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Myron J. Mitchell, Gary Lovett, Scott Bailey, Fred Beall, Doug Burns, Don Buso, Thomas A. Clair, Francois Courchesne, Louis Duchesne, Cathy Eimers, Ivan Fernandez, Daniel Houle, Dean S. Jeffries, Gene E. Likens, Michael D. Moran, Christopher Rogers, Donna Schwede, Jamie Shanley, Kathleen C. Weathers, Robert Vet. 2010-05-19. Comparisons of watershed sulfur budgets in southeast Canada and northeast US: New approaches and implications. https://doi.org/10.1007/s10533-010-9455-0

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Water table rise sustains carbon release from soils in wetland-dominated landscapes: An intact soil core study

Wetland-dominated landscapes influence carbon cycling through their potential to act as both carbon sinks and sources. Wetlands in low-relief landscapes have dynamic terrestrial-aquatic interfaces that change seasonally with variable surface water and groundwater levels. However, few studies have directly quantified dissolved organic matter (DOM) release and greenhouse gas (CO 2 , CH 4 ) fluxes from wetland soils along terrestrial-aquatic interfaces as they are seasonally re-saturated by groundwater. To estimate groundwater-mediated soil DOM and gas fluxes, we performed laboratory simulations of vertical groundwater rise on intact soil cores collected from four Delmarva bay wetlands located in the Mid-Atlantic United States. At each wetland, one core was collected from within the wetland basin and the other from the transitional zone at the basin edge. Cores were re-saturated with groundwater over 15 days and then kept fully saturated for an additional 25 days. Source groundwater, soil porewater, and exfiltrated surface water samples were collected and analyzed for pH, ORP, DOM concentration, and DOM optical indices. In both the wetland and transition zone cores, porewater DOM concentrations increased over the wet up and were sustained during prolonged saturation. Optical indices shifted from recently produced, microbial-like signatures towards aromatic, terrestrial-like signatures. Fluxes of CO 2 decreased as the duration of soil saturation increased and soil cores switched from CH 4 sinks to sources upon full soil core saturation. Results indicate that groundwater rise sustains carbon mobilization from soils in wetland-dominated landscapes, emphasizing the need to understand how climate-driven changes to groundwater dynamics may affect carbon fluxes along terrestrial-aquatic interfaces.

Delaware, Maryland↗

Seasonal synchronicity and multi-decadal stability of headwater biogeochemistry in the northern temperate zone

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.

Biogeochemistry↗

Rare earth element and phosphorus mobility depend on adsorption to Al-, Fe-, and Mn-oxyhydroxides in a headwater stream in Vermont, USA

We evaluated base cations, Al, Fe, Mn, REEs, DOC, anions, and P mobilization during three discharge events (E-1, E-2, E-3) at the headwater catchment W-9, Sleepers River, Vermont, USA. Peak discharge ranged from 3.696 (E-1) to 0.073 (E-3) mm h −1 . Eight samples from each event were speciated for total (unfiltered-acidified) and dissolved (0.45 µm filtered-acidified). During E-1, total Al, Fe, and Mn increased to maxima of 376 (1790X), 161 (194X), and 38 (45X) µmol L −1 , respectively. Concurrently, total La, Ce, Pr, and Nd increased to maxima of 87 (590X), 114 (671X), 15 (375X), and 53 (408X) nmol L −1 , respectively, greatly exceeding the dissolved fraction. Totals for Er > Yb > Tm > Lu reached comparable enrichments near or at maximum discharge. Discharge ranged from 0.41 to 0.94 mm h −1 during E-2, a snowmelt event. Total Al was comparatively stable; total Fe and Mn increased 15X and 79X, remaining less than total Al. Total La, Ce, Pr, and Nd peaked with total Fe and Mn. E-3, a late summer rain, resembled E-1 but had much lower maxima for all REEs, P, Al, Fe, and Mn. Particulate, total, and dissolved REEs and P correlated with DOC and with total Al, Fe, and Mn over all discharges, with molar Al > Fe > Mn (most samples) for E-1 and E-3. Maximum total and dissolved P declined from 15.4 and 0.89 (E-1) to 0.98 and 0.27 µmol L −1 (E-3), respectively. Particulate REEs correlated strongly (R 2 = 0.95–0.96) with Al, Fe, and Mn particulates eroded from the stream bed and continuously precipitated during and after high discharge of groundwater.

Vermont↗