Introduction
Release of nitrogen (N) to the atmosphere has been increasing globally due to greater emissions of N from human activities (
Vitousek et al. 1997;
Galloway et al. 2003,
2004). These N emissions react to form N species that are deposited onto ecosystems hundreds of kilometres from the source (
Asman et al. 1998). Nitrogen deposition is one component in the calculation of critical loads for acid sensitive lakes and soils (
Schulze et al. 1989;
Driscoll et al. 2001;
Jeffries et al. 2003). Although N is often a limiting element for forest growth, the loss of base cations from soils as a result of oxidized N loading can compromise forest response to N fertilization. At elevated levels of N deposition, forest productivity and biomass decrease, whereas emissions of the potent greenhouse gas, nitrous oxide, increase (
Aber et al. 1989;
Matson et al. 2002;
Venterea et al. 2003). Global N deposition modeling indicates that 50–80% of N deposition falls on natural (non-agricultural) areas and N accumulation is driving changes in terrestrial plant diversity (
Dentener et al. 2006;
Bobbink et al. 2010).
There are few studies on atmospheric NO
3 − isotopes in Canada (
Mayer et al. 2001;
Spoelstra 2004;
Spoelstra et al. 2001,
2004,
2007,
2010). Most published δ
15N−NO
3 − data are from heavily impacted systems in Europe, the eastern USA, northern Europe, and China (e.g.,
Hastings et al. 2003;
Elliott et al. 2007,
2009;
Zhang et al. 2008;
Fang et al. 2011;
Koszelnik and Gruca-Rokosz 2013;
Korth et al. 2014;
Yang et al. 2014;
Guerrieri et al. 2015;
Yu et al. 2016). In general, fewer studies with NH
4 + isotopes in atmospheric deposition are available (e.g.,
Garten 1992;
Zhang et al. 2008). Recently,
Holtgrieve et al. (2011) suggested that by more than doubling the reactive N in the biosphere and increasing atmospheric CO
2 concentration, the δ
15N of lake sediment has been altered. Even in remote areas, the deposition of N species and their associated δ
15N values have changed (
Dentener et al. 2006;
Knapp et al. 2008,
2010) and lake sediments have recorded a decrease in δ
15N over the last century (
Wolfe et al. 2003;
Hobbs et al. 2010;
Holtgrieve et al. 2011).
Here, we examine trends in N deposition and the stable isotopic values of N species at a site remote from urban and agricultural activities. The IISD–Experimental Lakes Area (ELA) in northwestern Ontario, Canada, has been the site of long-term whole-ecosystem research for over 45 years (
Blanchfield et al. 2009). Nitrogen deposition at the ELA is currently less than one third that of the most impacted sites in eastern North America and Europe (
Watmough et al. 2005;
Pardo et al. 2011). Our objectives in this paper are to (1) assess annual and open water season changes in the speciation of atmospheric N deposition at the ELA over 44 years of record; (2) characterize the seasonal variability in δ
15N−NO
3 −, δ
18O–NO
3 −, δ
15N−NH
4 +, and δ
15N–total nitrogen (TN) of atmospheric deposition at the ELA; and (3) compare these values with dissolved organic matter (DOM), particulate organic matter (POM), and zooplankton at the base of ELA food webs.
Methods
The ELA is situated on the Canadian Shield in the boreal forest of northwestern Ontario, Canada (93°41′W 49°41′N;
Fig. 1). The ELA climate is continental with cold winters and warm summers. Average annual precipitation is 707 mm/year (1970–2013). Climate and deposition measurements, including chemistry, have been made in the Lake 239 catchment since June 1969. In the early 1970s, N deposition at the ELA was among the lowest recorded in North America; however, this is no longer the case (
Parker et al. 2009). Bulk deposition was collected during the ice-free season at the Lake 239 island from 1969 to 1983 before being moved to the nearby Lake 240 island. There is no evidence indicating that the change in location affected the data. Precipitation has also been collected at the ELA Meteorological Station since 1970. Detailed methods are outlined by
Beaty (1981) and
Linsey et al. (1987). Briefly, in summer, 0.5 m × 0.5 m plexiglass bulk deposition collectors for wet and dry deposition screened to 100 μm were emptied after enough water was collected for analyses, filtered, and stored cold until chemical analyses were performed within a few days. In winter, snow from the collector was pushed with a clean paddle into a plastic bag and kept frozen until analysis.
Bulk deposition samples for isotopic analyses were collected at the meteorological site in a collector from 6 June 2012 to 21 March 2013. Samples from a bulk collector in the Lake 302 catchment, about 3 km away, collected from April to September 1996, were archived frozen (as per
Lamontagne and Schiff 1999) and analysed for this study.
Stream and lake survey samples of DOM, POM, and zooplankton were collected from the epilimnion of eight lakes (114, 239, 303, 304, 737, 626, 658, and 979) and five streams (water only) at the ELA (Lake 114 inflow, Lake 114 cliff inflow, Lake 240 inflow from Lake 470, Lake 302 upland 8, and Lake 979 east inflow) in mid-July 2010 (except Lake 114 inflow samples, which were collected during the summers of 2000–2002). Water samples were filtered with pre-combusted QMA filters (nominal pore size of 0.8 μm). Filtrate was acidified with concentrated HCl to pH 4 before being freeze-dried. Zooplankton were collected by towing a 50 cm diameter net with a 150 μm mesh behind a boat. Samples were washed and freeze-dried. Both freeze-dried material and filters were analysed for δ15N on a Carlo Erba 1108 elemental analyser coupled to a Thermo Finnigan Delta+ continuous flow isotope-ratio mass spectrometer (EA-IRMS).
ELA deposition samples were analysed by standard methods (
Stainton et al. 1977). NO
2 − was analysed colorimetrically by the azo dye method on a Technicon Autoanalyzer. NO
3 − was analysed colorimetrically following reduction to NO
2 − via a copper-cadmium couple by the azo dye method on a Technicon Autoanalyzer. NH
4 + was analysed colorimetrically by the indophenol blue method with nitroprusside catalyst on a Technicon Autoanalyzer. Total dissolved nitrogen (TDN) was analysed by photo-oxidation of alkaline samples, and the subsequent NO
3 − was reduced to NH
4 + by zinc in acid, then analysed colorimetrically as per NH
4 +. Dissolved organic nitrogen (DON) was calculated as the difference between TDN and the sum of NO
3 −, NO
2 −, and NH
4 +. Prior to 1989, NO
2 − was not reported separately from NO
3 − and the values reported here for NO
3 − include NO
2 −. After 1989, we used only the NO
3 − values. NO
2 − values in bulk collected deposition are very small, averaging 1.4% of total NO
3 − + NO
2 − from 1990 to 2013. NO
2 − is included in the TN and TDN analyses in the entire data record. Suspended N was collected by filtration on a glass fibre filter and analysed on an elemental analyser.
Bulk deposition samples in 2010–2012 were collected at the ELA meteorological site in plastic bags in about two week intervals. Bags were sealed, kept cold, screened to 150 μm, transferred to Nalgene bottles, and frozen within 2 d. Samples were kept frozen until analysis. Bulk deposition samples at the Lake 302 uplands were collected in 1996 on a storm event basis and were treated similarly.
For isotopic analysis, NO
3 − + NO
2 − concentrations were determined using a Westco SmartChem 200 discrete analyzer with a method based on
USEPA 353.2 Revision 2.0 (1993). NO
3 − was reduced to NO
2 − by passage of the sample through a tubular copperized cadmium reactor from which NO
2 − was treated with sulphanilamide and N-(naphthyl)-ethylenediamine dihydrochloride to form a dye measured colorimetrically at 550 nm. Precision was ±0.02 mgN/L. NH
4 + concentrations were analysed manually using spectrophotometer on unfiltered samples. NH
4 + was reacted with alkaline phenol and then hypochlorite-forming indophenol blue, which was intensified by adding sodium nitroprusside before analysis at 600 nm. Precision was ±0.08 mgN/L.
For δ
15N and δ
18O–NO
3 − analysis, NO
3 − + NO
2 − was chemically reduced to N
2O (
McIlvin and Altabet 2005). Briefly, NO
3 − was reduced to NO
2 − with cadmium and then NO
3 − + NO
2 − was reduced to N
2O with NaN
3. N
2O was analysed with a VG IsoPrime continuous flow isotope-ratio mass spectrometer with a VG TraceGas pre-concentrator. Analyses were performed in duplicate. Precisions were ±0.2‰ for δ
15N−NO
3 − and ±0.5‰ for δ
18O–NO
3 −.
For δ
15N−NH
4 + analysis, NH
4 + was chemically converted to N
2O (
Zhang et al. 2007). Briefly, NH
4 + was oxidized to NO
2 − with BrO
− and then reduced to N
2O with NaN
3. The resulting N
2O was analysed as above in duplicate. Precision was ±0.3‰.
For δ15N–TN analysis, sample volumes were reduced by evaporation from about 2 L to about 100 mL. Samples were freeze-dried and packed into tin cups for analysis on the EA-IRMS as above. Precision was ± 0.3‰.
Concentrations and precipitation depths were combined to yield areal deposition rates on an annual basis and for the open water season of late-April to the end of September (ice-off to thermal destratification, d 120–273). There are a few external inputs of N to these small headwater lakes other than atmospheric deposition once they are thermally stratified and after snowmelt. Long-term trends were assessed with the non-parametric Mann–Kendall test (
McLeod 2011) in R (
R Core Team 2016).