Introduction
The contamination of soil and groundwater with nitrate from agricultural operations is a global water quality issue that has been extensively documented (Power and Schepers, 1989; Spalding and Exner, 1993; Rodvang and Simpkins, 2001; Galloway et al., 2008; Zirkle et al., 2016; Arauzo, 2017; Ascott et al., 2017). Leaching of nitrate () from animal waste or fertilisers can result in groundwater concentrations that exceed drinking water guidelines and pose human health risks (Fan and Steinberg, 1996; Gulis et al., 2002; Yang et al., 2007). The discharge of high- groundwater, runoff, or drainage can contaminate streams and lakes, resulting in eutrophication and ecosystem decline (Deutsch et al., 2006; Kaushal et al., 2011). In saturated groundwater systems with low oxygen concentrations, elevated can be naturally attenuated by microbial denitrification (Wassenaar, 1995; Robertson et al., 1996; Smith et al., 1996; Tesoriero et al., 2000; Singleton et al., 2007). Concentrations of will also decrease along groundwater flow paths due to attenuation via dilution by hydrodynamic dispersion (referred to hereafter as mixing). Because of these natural attenuation mechanisms, infiltration to groundwater may be preferable to off-site drainage and runoff of nitrate-rich waters. Many agricultural operations are undertaken on fertile soils associated with glacial sediments (Spalding and Exner, 1993; Ernstsen et al., 2015; Zirkle et al., 2016). Understanding the sources and fate of agriculturally derived nitrate in groundwater systems in glacial sediments is therefore critical for managing impacts of human food production on the environment.
Identification of the sources and fate of at agricultural operations can be challenging because of spatial and temporal variations in sources (e.g. earthen manure storage, temporary manure piles, or fertiliser) and heterogeneity in hydrogeologic systems (Spalding and Exner, 1993; Rodvang et al., 2004; Showers et al., 2008; Kohn et al., 2016). These spatial and temporal variations can result in complex subsurface solute distributions that are difficult to interpret using classical transect studies or numerical groundwater models (Green et al., 2010; Baily et al., 2011).
Groundwater containing significant agriculturally derived also typically has elevated chloride () concentrations (Saffigna and Keeney, 1977; Rodvang et al., 2004; Menció et al., 2016). Decreasing (or ) ratios have been used to define denitrification based on the assumption that is reactive while is non-reactive (conservative), such that denitrification results in a decrease in the ratio (Kimble et al., 1972; Weil et al., 1990; Liu et al., 2006; McCallum et al., 2008). However, ratios can also change in response to mixing of groundwater with different ratios or when groundwater sampling traverses hydraulically disconnected formations (Bourke et al., 2015b). If ratios vary among potential sources and the ratio at the point of entry to the groundwater system can be reconstructed, this information could be used to show that anthropogenic at different locations within an aquifer is derived from the same or different sources.
The stable isotopes of ( and ) provide an alternative approach to characterising the source and fate of in groundwater systems. In agricultural areas, multiple sources of are common and could include precipitation, soil , inorganic fertiliser, manure, and septic waste (Komor and Anderson, 1993; Liu et al., 2006; Pastén-Zapata et al., 2014; Clague et al., 2015; Xu et al., 2015). While source identification is theoretically possible using and (particularly with a dual-isotope approach), in practice this can be difficult due to geologic heterogeneity, overlapping source values, and the complexity of biologically mediated reactions (Aravena et al., 1993; Wassenaar, 1995; Mengis et al., 2001; Choi et al., 2003; Granger et al., 2008; Vavilin and Rytov, 2015; Xu et al., 2015).
attenuation by denitrification in groundwater systems can be identified based on the characteristic enrichment of and . Numerous studies have made qualitative assessments that identified denitrification in groundwater using the stable isotope approach (Böttcher et al., 1990; Wassenaar, 1995; Singleton et al., 2007; Baily et al., 2011; Clague et al., 2015; Xu et al., 2015). Recently published papers have also used stable isotopic values of and water as the basis for mixing models in agricultural settings (Ji et al., 2017; Lentz and Lehersch, 2019). Isotopic fractionation effects can also allow for quantitative assessment of the proportion of substrate that has undergone a given reaction, if enrichment factors and source values are known; as in the case of evaporative loss of water, for example (Dogramaci et al., 2012). To date, there have been very few attempts to quantify denitrification using dual-isotope enrichment, largely due to uncertainty in source values and enrichment factors (Böttcher et al., 1990, Xue et al., 2009).
The only published calculations of the fraction of remaining after denitrification the that we are aware of assumed a constant enrichment factor and the same isotopic source values across the field site (Otero et al., 2009). However, the enrichment factor will vary across a field site in response to reaction rates (Kendall and Aravena, 2000), and isotopic values of even the same type of source (e.g. manure) can vary substantially (Xue et al., 2009).
If the variation in source values and enrichment factors can be characterised from measured data then these uncertainties can be accounted for using a Monte Carlo approach (Joerin et al., 2002; Bourke et al., 2015a; Ji et al., 2017), thereby extending the application of the dual-isotope technique to allow for a robust quantitative assessment of denitrification in agricultural settings.
A synthesised analysis of stable isotopes of with additional ionic tracers can further improve the assessment of attenuation mechanisms and sources of NOin agricultural settings (Showers et al., 2008; Vitòria et al., 2008; Xue et al., 2009; Xu et al., 2015; Ji et al., 2017). We hypothesise that if the amount of denitrification can be quantified based on and , then this estimate of the fraction of removed through denitrification can be used to constrain a mixing model based on and concentrations. This novel approach allows for the ratio of at the point of entry to the groundwater system to be reconstructed from measured and concentrations (see Sect. 2.4). Where the ratio varies between sources, this ratio can then be used to assess the source of the in groundwater (e.g. temporary manure piles or feeding pens). These data can also then be used to estimate the initial concentrations of and at the point of entry to the groundwater system and quantify attenuation by mixing.
Figure 1Map of study sites CFO1 and CFO4, showing locations of groundwater monitoring wells, core collection, earthen manure storages (EMSs), dairy and feedlot pens, manure piles, and irrigated land. Blue rectangle indicates extent of CFO1 inset.
[Figure omitted. See PDF]
In this study, we present the application of this approach at two confined feeding operations (CFOs) in Alberta, Canada, with differing lithologies and durations of operation (Fig. 1). Concentrations of and nitrogen species (N species) and the stable isotopes of were measured in groundwater samples collected from monitoring wells and continuous soil cores, as well as manure filtrate at both sites. These data were interpreted to (1) assess the extent of agriculturally derived in groundwater, (2) identify sources and initial concentrations of at the point of entry to the groundwater system, and (3) assess mixing and denitrification as attenuation mechanisms at these sites.
Materials and methods
Experimental sites
This study was conducted using data from two of the five sites investigated by Alberta Agriculture and Forestry during an assessment of the impacts of livestock manure on groundwater quality (Lorenz et al., 2014). To the best of our knowledge (including discussions with farm operators) fertilisers have not been applied at either of these sites. As such, manure waste from livestock is assumed to be the sole source of agricultural nitrogen (N) and elevated concentrations in groundwater at these sites.
The first study site (CFO1) is located 25 km northeast of Lethbridge, Alberta (Fig. 1). Agricultural operations at this site were initiated with the construction of a dairy in 1928, which has the capacity for 150 dairy cattle. A feedlot for beef cattle was added in 1960s along with an earthen manure storage (EMS) facility for storing liquid dairy manure (approx. 4 m deep) and a catch basin that receives surface water runoff. This feedlot was expanded in the 1980s to the 2000-head capacity it was at the time of this study. There is also a dugout (or slough, a shallow wetland) on-site that receives local runoff and an irrigation drainage canal at the southern boundary of the property.
The second study site (CFO4) is located approximately 30 km north of Red Deer, Alberta, and 300 km north of CFO1. This dairy and associated EMS (approx. 6 m deep) were constructed in 1995 and the facility had 350 head of dairy cattle at the time of the study. Runoff will drain either to the small dugout in the northwest of the site, or the natural drainage features (ephemeral ponds or a creek approx. 1.5 km east).
Sampling and instrumentation
Groundwater monitoring wells
Groundwater samples were collected from water table wells and piezometers (hereafter both are referred to as wells) installed at both sites (Table 1). At CFO1, groundwater samples were collected from six individual water table wells (DMW1, DMW2, DMW3, DMW4, DMW5, DMW6) and eight sets of nested wells with one well screened at the water table and one well screened 20 m below ground (BG) (DP10-2 and DP10-1, DMW10 and DP11-10b, DMW11 and DP11-11b, DMW12 and DP11-12b, DMW13 and DP11-13b, DMW14 and DP11-14b, DMW15 and DP11-15b, and DMW16 and DP11-16b). Wells DP10-2 and DP10-1 were located directly adjacent to the EMS on the hydraulically down-gradient side. At CFO4, groundwater samples were collected from eight water table wells (BC1, BC2, BC3, BC4, BC5, BMW1, BMW3, BMW7) and four sets of nested wells, with wells screened across the water table and at 15 m BG. Two of these nests were located adjacent to the EMS (BMW2 and BP10-15e, BMW4 and BP10-15w) and two were hydraulically down-gradient of the EMS (BMW5 and BP5-15, BMW6 and BP6-15).
Table 1Details of groundwater monitoring wells and continuous core collection at CFO1 and CFO4 (all screens installed at bottom of the well).
Site | Well/core | Type | Lateral | Ground | Total | Screen | Lithology of | (m s) |
---|---|---|---|---|---|---|---|---|
hole ID | distance | elevation | depth (m | length | screened interval | |||
from | (m a.s.l.) | below | (m) | |||||
EMS (m) | ground) | |||||||
CFO1 | DMW1 | WTW | 60 | 869.7 | 5.0 | 4.0 | Sand | |
DMW2 | WTW | 10 | 867.2 | 6.0 | 4.0 | Sand | ||
DMW3 | WTW | 2 | 867.5 | 3.7 | 2.0 | Sand | ||
DMW4 | WTW | 160 | 4.2 | 4 | Sand | |||
DMW5 | WTW | 270 | 866.4 | 6.8 | 4.0 | Clayey sand | ||
DMW6 | WTW | 310 | 6.7 | 4 | ||||
DP10-1 | Piezo | 2 | 867.8 | 18.6 | 0.5 | Clay | ||
DP10-2 | Piezo | 2 | 867.9 | 8.0 | 1.5 | Sand | ||
DMW10 | WTW | 340 | 868.0 | 7.2 | 3.0 | Clay | ||
DP11-10b | Piezo | 340 | 868.0 | 20 | 0.5 | Clay | ||
DMW11 | WTW | 470 | 864.8 | 7.0 | 3.0 | Sand and clay | ||
DP11-11b | Piezo | 470 | 20 | 0.5 | Clay | |||
DMW12 | WTW | 50 | 867.6 | 7.0 | 3.0 | Sand and clay | ||
DP11-12b | Piezo | 50 | 867.6 | 20.1 | 1.0 | Clay | ||
DMW13 | WTW | 35 | 867.1 | 7.0 | 3.0 | Sand | ||
DP11-13b | Piezo core | 35 | 867.1 | 20.0 | 0.5 | Clay | ||
DMW14 | WTW | 105 | 865.7 | 7.0 | 3.0 | Clay | ||
DP11-14b | Piezo | 105 | 865.7 | 20.0 | 0.5 | Sand | ||
DMW15 | WTW | 185 | 7.0 | 3 | Clay | |||
DP11-15b | Piezo | 185 | 20.0 | 0.5 | Clay | |||
DMW16 | WTW | 320 | 866.0 | 6.0 | 3.0 | Sand and clay | – | |
DP11-16b | Piezo | 320 | 20.0 | 0.5 | Clay | |||
DC15-20 | Core | 76 | 15 | |||||
DC15-21 | Core | 45 | 10.5 | |||||
DC15-22 | Core | 22 | 12 | |||||
DC15-23 | Core | 9 | 15 | |||||
CFO4 | BC1 | WTW | 110 | 857.0 | 6.9 | 3.1 | Clay and sandstone | |
BC2 | WTW | 365 | 859.4 | 7.0 | 3.1 | Clay and sandstone | ||
BC3 | WTW | 145 | 858.6 | 6.8 | 3.1 | Clay and sandstone | ||
BC4 | WTW | 95 | 858.8 | 5.9 | 3.0 | Clay and sandstone | ||
BC5 | WTW | 105 | 859.5 | 7.5 | 4.5 | Clay and sandstone | ||
BMW1 | WTW | 4 | 858.6 | 7.1 | 3.1 | Clay and sandstone | ||
BMW2 | WTW | 3 | 857.9 | 7.5 | 4.5 | Clay and sandstone | ||
BMW3 | WTW | 8 | 858.6 | 6.0 | 3.0 | Clay and sandstone | ||
BMW4 | WTW | 14 | 858.0 | 7.5 | 4.8 | Clay and sandstone | ||
BMW5 | WTW | 60 | 858.0 | 7.5 | 4.5 | Clay and sandstone | ||
BP5-15 | Piezo | 60 | 858.1 | 15.3 | 1.5 | Sandstone | ||
BMW6 | WTW | 150 | 856.9 | 7.5 | 4.5 | Clay and sandstone | ||
BP6-15 | Piezo | 150 | 856.8 | 15.2 | 1.5 | Sandstone | ||
BMW7 | WTW | 140 | 856.7 | 7.5 | 4.5 | Clay and sandstone | ||
BP10-15e | Piezo | 4 | 858.2 | 14.9 | 1.5 | Sandstone | ||
BP10-15w | Piezo | 10 | 858.0 | 15.0 | 1.5 | Sandstone |
WTW: water table well, Piezo: piezometer, Core: continuous core. EMS: earthen manure storage.
Groundwater samples were collected for ion analysis ( and N species) quarterly between April 2010 and August 2015. All water samples were collected using a bailer after purging (1–3 casing volumes) and stored at C prior to analysis. Samples for and were collected from wells at CFO1 on 1 January and 1 May 2013. Samples for and at CFO4 were collected on 27 October 2014. Wells were purged prior to sample collection (1–3 casing volumes), and samples filtered into high-density polyethylene (HDPE) bottles in the field and frozen until analysis.
Hydraulic heads in monitoring wells were determined using manual measurements (approximately monthly, 2010–2015). Hydraulic head response tests were conducted on the majority of the wells at the sites to determine hydraulic conductivity () of the formation media surrounding the intake zone. These tests were either a slug test (water level decline after water addition) or bail test (water level recovery after water removal) depending on the location of the water level within the well at the time of testing. was determined from the hydraulic head responses using the method of Hvorslev (1951).
Continuous core
A continuous core was collected at CFO1 immediately adjacent to well DP11-13b on 1 May 2013 (Fig. 1). Additional core samples were collected from 1 to 5 June 2015 along a transect hydraulically down-gradient of the southeastern side of the EMS at CFO1, where hydrochemistry data suggested leakage from the EMS (see Sect. 3). During this 2015 drilling campaign, core samples were collected at four locations (DC15-20, DC15-21, DC15-22, DC15–23) to depths of up to 15 m below surface and distances of up to 100 m from the EMS between wells DMW3 and DP11–14.
Continuous core samples were retrieved using a hollow stem auger (1.5 m core lengths) with 0.3 m sub-samples collected at approximately 1 m intervals ensuring that visually consistent lithology could be sampled. Core samples for were stored in Ziploc bags and kept cool until analysis. Core samples for N-species analysis were stored in Ziploc bags filled with an atmosphere of argon (99.9 % Ar) to minimise oxidation and kept cool until analysis. Subsamples of each core (250–300 g) were placed under 50 MPa pressure in a Carver Auto Series NE mechanical press with a 0.5 m filter placed at the base of the squeezing chamber, which was placed within an Ar atmosphere to minimise oxidation. A syringe was attached to the base of the apparatus and 15 mL of filtered pore water were collected for analyses within 3.5 to 6.0 h (Hendry et al., 2013).
Liquid manure storages
Samples of liquid manure slurry were collected directly from the EMS at both sites and the catch basin (containing local runoff from the feedlot) at CFO1 using a pipe and plunger apparatus to sample from approximately 0.5 m below the surface. The slurry collected was subsequently filtered (0.45 m) to separate the liquid and solid components. The water filtered from samples collected from the EMS or catch basin is hereafter referred to as manure filtrate.
Laboratory analysis
Groundwater samples from wells were analysed by Alberta Agriculture and Forestry (Lethbridge, Alberta). Concentrations of were determined using potentiometric titration of , with a detection limit of 5.0 mg L and accuracy of 5 % (APHA 4500-). Concentrations of as N (), as N (), and as N () were measured by air-segmented continuous-flow analysis (APHA 4500-NH3 G, APHA 4500-NO3 F). Total nitrogen (TN) was determined by high temperature catalytic combustion and chemiluminescence detection using a Shimadzu TOC-V with attached TN unit (ASTM D8083-16). Total organic nitrogen (TON) was calculated by subtracting , , and from TN. Bicarbonate () was analysed by titration (APHA 2320 B). Dissolved organic carbon (DOC) was analysed by a combustion infrared method (APHA 5310 B) using a Shimadzu TOC-V system. Manure filtrate was analysed by ALS (Saskatoon, Saskatchewan) using similar methods for (APHA 4110 B), TN (RMMA A3769 3.3), as N (APHA 4500-NO3-F), (APHA 4500-NH3 D), (APHA 2320), and DOC (APHA 5310 B).
Pore-water samples squeezed from the continuous core were analysed at the University of Saskatchewan (Saskatoon, Canada) for , , and using a Dionex IC25 ion chromatograph (IC) coupled to a Dionex As50 autosampler (EPA Method 300.1, accuracy and precision of 5.0 %) (Hautman and Munch, 1997). Ammonia as N () was measured by Exova laboratories using the automated phenate method (APHA Standard 4500-NH3 G, detection limit of 0.025 mg L, accuracy of 2 % of the measured concentration, and a precision of 5 % of the measured concentration).
and in groundwater samples (from wells and pore water from the continuous core) and manure filtrate were measured at the University of Calgary (Calgary, Alberta) using the denitrifier method (Sigman et al., 2001) with an accuracy and precision of 0.3 ‰ for and 0.7 ‰ for . Groundwater samples collected for isotope analysis in January 2013 were also analysed for by the University of Calgary (denitrifier technique, Delta XL).
Modelling approach
Quantification of denitrification based on and
Nitrate in groundwater that has undergone denitrification is commonly reported as being identified by enrichment of and with a slope of about 0.5 on a cross-plot (Clark and Fritz, 1997). However, published studies of denitrification in groundwater report slopes of up to 0.77 (Mengis et al., 1999; Fukada et al., 2003; Singleton et al., 2007). The relationship between isotopic enrichment of and and the fraction of remaining during denitrification can be described by a Rayleigh equation: where is the initial isotope ratio (relative to the standard) of the ( or ), is the isotopic ratio when fraction of remains, and is the kinetic fractionation factor () (Böttcher et al., 1990; Clark and Fritz, 1997; Otero et al., 2009; Xue et al., 2009). Kinetic fraction effects are commonly also expressed as the enrichment factor, . In the case of a constant enrichment factor, can be calculated from measured (or ), if the initial () is known; The fraction of removed from groundwater through denitrification is then given by (). The concentration of that would have been measured if mixing was the only attenuation mechanism () can also be calculated by dividing the measured concentration by .
A subset of 20 samples with isotopic values of indicative of denitrification were identified, and for each of these samples (mean and standard deviation) was calculated from Eq. (2) using a Monte Carlo approach with 500 realizations. The distribution of values was defined based on measured data. If the initial is known, for () can be determined from the slope of the linear regression line on a plot of vs. (Böttcher et al., 1990). If the initial and are not known, as is the case here, can be determined from the slope of the regression line on a plot of () vs. , which will be the same as on a plot of vs. . In situ variations in temperature and reaction rates may affect the enrichment factor (Kendall and Aravena, 2000) and this was accounted for by allowing for variation in within the Monte Carlo analysis. The enrichment factor for () was calculated by multiplying the by a linear coefficient of proportionality determined for each CFO from the slope of the denitrification trend on an isotope cross-plot (see Sect. 3.2).
For each realisation, initial isotopic values ( and ) were determined by Excel Solver such that the difference between calculated from and was minimised ( % difference). The ranges of and were limited based on measured data and literature values (see Sect. 3.2). This approach neglects the effect of mixing of groundwater with differing isotopic values and is valid if the concentration of in the source is much greater than background concentrations such that the isotopic composition of is dominated by the agriculturally derived end-member.
Quantification of mixing and initial concentrations of and
A binary mixing model that also accounts for decreasing concentrations in response to denitrification was used to quantify attenuation by mixing and estimate the initial concentrations of and . The measured concentration of was assumed to be a function of two end-members mixing, described by where Cl is the measured concentration of in the groundwater sample, Cl is the concentration of at the initial point of entry of the agriculturally derived to the groundwater system, Cl is the concentration of in the background ambient groundwater, and is the fraction of water in the sample from the source of agriculturally derived (and ) remaining in the mixture.
The concentration of was also assumed to be a function of two end-members mixing but with an additional coefficient, (the fraction of remaining after denitrification), applied to account for denitrification. The measured concentration was thus described by where is the concentration of measured in the groundwater sample, is the concentration of in the source of agriculturally derived at the initial point of entry to the groundwater system, and is the concentration of in the background ambient groundwater. This mixing calculation was only conducted on samples for which dominated total-N ( 10 % of ) so that nitrification of could be neglected.
If Cl is much greater than Cl and is much greater than , then is insensitive to background concentrations and these terms can be neglected (see Sect. 4.2 for further discussion of this assumption). In this case, Eqs. (3) and (4) reduce to Solving Eq. (6) for and substituting into Eq. (5) yields Thus, for each groundwater sample, the ratio of at the initial point of entry of the agriculturally derived to the groundwater system can be simply calculated using measured concentrations, and estimated from isotope data. This provides a relatively simple method to identify agriculturally derived from different sources (e.g. EMS vs. manure piles) if they have different ratios. Estimated Cl and are reported as the mid-range value with uncertainty described by the minimum and maximum values. These initial concentrations are at the water table for top-down inputs, or at the saturated point of contact between the EMS and the aquifer for leakage from the EMS. This analysis assumes that a sampled water parcel consists of water with agriculturally derived that entered the aquifer from one source at one point in time and space and has since mixed with natural ambient groundwater. Any produced during nitrification after the anthropogenic source water enters the aquifer is implicitly included in . The error in was assumed to be dominated by error in the estimated , with the measurement error in and considered negligible.
The initial concentrations of the agriculturally derived source ( and Cl) were estimated by simultaneously solving Eqs. (5) and (6) using Excel Solver (GRG nonlinear). The absolute minimum values of and Cl were defined by measured concentrations (e.g. if Cl, ). Maximum values of and Cl were defined based on measured concentrations of and in groundwater and manure filtrate ( 150 mg L and 1300 mg L; see Sect. 3.2). These maximum values of and Cl correspond to the minimum . The value of was assumed to be the mean estimated from isotopes using Eq. (2), and was required to be within 1 standard deviation of the estimate from Eq. (7).
The resulting estimates of are reported as the mid-range, with uncertainty described by the minimum and maximum values. Larger values of indicate less mixing (a shorter path for advection–dispersion) and suggest a source close to the well. Smaller values of indicate extensive mixing (a longer path for advection–dispersion) and suggest a source further away from the well. The relative contributions of mixing and denitrification to attenuation at each site were evaluated by comparing and for each sample. This analysis was conducted using isotope values from the samples collected on 1 May 2013 at CFO1, which were combined with the and data from 6 June 2013. At CFO4, results from stable isotopes collected on 27 October 2014 were combined with and data collected on 7 October 2014.
Results
Site hydrogeology
CFO1
The geology at CFO1 consists of clay and clay–till interspersed with sand layers of varying thickness to the maximum depth of investigation (20 m BG, bedrock not encountered). Hydraulic conductivities () calculated from slug tests on wells ranged from to m s () for sand, to m s () for clay–till, and to m s () for clay. Depth to the water table throughout the study site ranged from 0.5 m at DMW14 to 3.8 m at DMW11. Seasonal water table variations were about 0.5 m with no obvious change in the annual average during the 6-year measurement period. Water table elevation was highest at DMW10 and DMW1 on the west side of the site and lowest at DMW11 on the northeast side of the site (see Supplement). Measured heads indicate groundwater flow from the vicinity of the EMS to the northeast and southeast. Mean horizontal hydraulic gradients at the water table ranged from to m m. Vertical gradients were predominantly downward in the upper 20 m of the profile (mean gradients ranging from to 0.18 m m), with the exception of DMW11 where the vertical gradient was upward (mean gradient m m). Using the geometric mean for the sand ( m s) and a lateral head gradient of m m yields a specific discharge (Darcy flux, ) of 2.2 m yr. Assuming an effective porosity of 0.3 (Rodvang et al., 1998), the average linear velocity () is 7.4 m yr. This suggests that, in the absence of attenuation by mixing or denitrification, agriculturally derived could have been transported through the groundwater system by advection about 400 m since 1960 and 630 m since 1930.
CFO4
The geology at CFO4 consists of about 5 m of clay (with minor till) underlain by sandstone, to the maximum depth investigated (20 m BG). Hydraulic conductivities measured using slug tests on wells were to m s () for the clay and sandstone (many shallow wells were screened across the clay–till and into the sandstone) and to m s () for the sandstone. The depth to water table ranged from 1.0 to 3.4 m, increasing from west to east across the study site. Seasonal water table variations were on the order of 1.5 m with water table declines on the order of 0.3 m yr. The horizontal hydraulic gradient was consistently from west to east, with a mean gradient at the water table of m m between BC2 and BMW2 and m m between BMW2 and BMW7. Vertical hydraulic gradients were to m m downward. Using the geometric mean for the site ( m s) and a lateral head gradient of m m yields a of 0.4 m yr. Assuming an effective porosity of 0.3 yields a of 1.3 m yr. These values suggest that, in the absence of attenuation by mixing or denitrification, anthropogenic could have been transported through the groundwater systems about 10 m by advection between 1995 and the time of sampling.
Values and evolution of stable isotopes of nitrate
The range of isotopic values of in groundwater was similar at both sites (Fig. 2). At CFO1, ranged from to 20.1 ‰ and from to 61.0 ‰. At CFO4, ranged from to 31.6 ‰ and from to 70.5 ‰. The isotopic values of in groundwater are commonly assumed to be derived from a mix of a one-third atmospheric-derived oxygen ( ‰) and two-thirds water-derived oxygen (Xue et al., 2009). Given the average for both sites ( ‰; see Supplement), a one-third atmospheric two-thirds groundwater mix would result in a of ‰. Manure filtrate from the EMS at CFO1 had ranging from 0.4 to 5.0 ‰ and ranging from 7.1 to 19.0 ‰. A curve showing the co-evolution of (mixing of atmospheric with groundwater-derived ) and (Rayleigh distillation, ) during nitrification is shown in Fig. 2. Isotopic values in DMW3, where direct leakage from the EMS was evident, are consistent with partial nitrification following this trend of isotopic evolution ( of ‰ and of 7.8 ‰).
Figure 2(a) Cross-plot of stable isotopes of nitrate at CFO1 and CFO4 showing hypothetical nitrification trend, boundary of manure-sourced values and linear enrichment trends associated with denitrification. (b) Enrichment of during denitrification (only samples within source region and with evidence of denitrification are shown); dashed lines represent 1 SD of enrichment factor ().
[Figure omitted. See PDF]
At both sites, co-enrichment of and characteristic of denitrification was evident in some samples (slopes of 0.42 and 0.72 in Fig. 2a). At CFO1, this includes samples from DP10-2, DMW5, DMW11, DMW12, DP11-12b, and DMW13 (and associated core) and some pore water from cores DC15-22 and DC15-23. These samples had concentrations of 0.6 to 23.7 mg L, ranging from 4.8 to 20.6 ‰, and ranging from 22.9 to 61.3 ‰. At CFO4, samples exhibiting evidence of denitrification were from BMW2, BMW5, BMW6, BMW7, and BC4. These samples had concentrations ranging from 0.4 to 35.1 mg L, ranging from 1.6 to 22.1 ‰, and ranging from 20.9 to 70.1 ‰. Although the isotopic values of DMW5 suggest enrichment by denitrification, the data plot away from the rest of the CFO1 data and close to the denitrification trend at CFO4 (Fig. 2), suggesting these samples were affected by some other process (possibly mixing or nitrification); therefore, was not calculated. Also, well DMW3, which clearly receives leakage from the EMS, did not contain substantial and so was not calculated.
In the Monte Carlo analysis the potential range of original isotopic values of the source prior to denitrification ( and ) varied from 5 to 27 ‰ for and from to 7 ‰ for based on isotopic values measured during this study (Fig. 2a). These values are consistent with literature values for manure-sourced , which report ranging from 5 to 25 ‰ and ranging from to 5 ‰ (Wassenaar, 1995; Wassenaar et al., 2006; Singleton et al., 2007; McCallum et al., 2008; Baily et al., 2011). was defined by a normal distribution with a mean of ‰ and standard deviation of 2.5 ‰ (Fig. 2b). At CFO1, the coefficient of proportionality between the enrichment factor of and was described by a normal distribution with mean of 0.72 and standard deviation of 0.05. At CFO4, the coefficient of proportionality was also described by a normal distribution with a mean of 0.42 and standard deviation of 0.035 (see Fig. 2a). These enrichment factors are consistent with values from denitrification studies that report ranging from to ‰ and ranging from to ‰ (Vogel et al., 1981; Mariotti et al., 1988; Böttcher et al., 1990; Spalding and Parrott, 1994; Mengis et al., 1999; Pauwels et al., 2000; Otero et al., 2009).
Distribution and sources of agricultural nitrate in groundwater
At both sites TN concentrations in filtrate from the EMS and catch basin were generally an order of magnitude larger than concentrations in groundwater (Table 2). The one exception is well DMW3 at CFO1, which intercepted direct leakage from the EMS (see Sect. 3.3.1 for further discussion of this well). The dominant form of N differed between manure filtrate and groundwater. In the EMS filtrate, N was predominately organic N (TON up to 71 %) or (up to 90 %), with 0.1 % of TN. In the catch basin at CFO1 TON was % of TN. In groundwater TN concentrations ranged from to 84.6 mg L, and this N was predominantly (again, with the exception of DMW3).
Table 2Range of measured concentrations of TN, , (), and TON at each study site. At CFO1 results from monitoring well DMW3 are presented separately because values in this well differed substantially from all other wells.
Site | N pool | TN | TON | ||
---|---|---|---|---|---|
(mg L) | (mg L) | (mg L) | (mg L) | ||
CFO1 | EMS | 550–1820 | 275–747 | –0.4 | 73–1301 |
Catch basin | 200–1440 | 2.5–7.3 | 196–1437 | ||
DMW3 | 278–548 | 219–479 | –50 | 31.3–73.9 | |
Other monitoring wells | –33.4 | –2.9 | –31.4 | –3.7 | |
CF04 | EMS | 1000–1240 | 724–747 | 0.25–0.29 | 275–492 |
Monitoring wells | –84.6 | –0.23 | –80.4 | –13.9 |
of 50 mg L in DMW3 consisted of 12.6 mg L as and 37.4 mg L as . max in groundwater was measured in core ( 66.4 mg L, 67.8 mg L). Range across three replicates was measured on 25 August 2011.
CFO1
Agriculturally derived was generally restricted to the upper 20 m (or less) at CFO1 ( 0.2 mg L and 57 mg L in seven wells screened at 20 m). The one exception was DP11-12b, which had up to 4.1 mg L of . The southeast portion of the site also does not appear to have been significantly contaminated by agriculturally derived , with concentrations 1 mg L in five water table wells (DMW4, DMW6, DMW14, DMW15, DMW16). In DMW6, and TN concentrations were elevated (see Supplement) but concentrations were mg L. Collectively, these data suggest the catch basin is not a significant source of to the groundwater at this site.
Figure 3Temporal variations in (a) , (b) , and (c) at CFO1. Only wells with 10 mg L are shown.
[Figure omitted. See PDF]
Leakage of manure slurry from the EMS at CFO1 is clearly indicated by the data from DMW3, which feature the highest concentrations of TN in groundwater (up to 548 mg L) and elevated , , and DOC in concentrations similar to EMS manure filtrate (see Supplement). Nevertheless, concentrations in this well were consistently low ( mg L, ). The potential for nitrification in the vicinity of this well is indicated by production ( mg L, ). However, the data demonstrate that only a small proportion of the in DMW3 ( mg L, ) could have been converted to within the subsurface ( in groundwater 66 mg L). Further work is required to assess the importance of cation exchange as an attenuation mechanism for direct leakage from the EMS at this site.
Figure 4(a) Estimated ratios (mean and SD) in water table wells with evidence of denitrification at CFO1, plotted with distance from earthen manure storage (EMS), where dashed lines are the upper and lower bounds of DP10-2 (EMS source) and labelled values are maximum measured (mg L. (b) Estimated concentrations of and at CFO1 (mid-range, error bars are max and min values).
[Figure omitted. See PDF]
Contamination by agricultural that exceeds the drinking water guidelines ( 10 mg L) was observed in four wells (DMW1, DMW11, DMW13, and DP10-2) and in the continuous core (DC15-23) (Fig. 3). DMW2 and DMW12 also had concentrations that were elevated but did not exceed the drinking water guideline ( mg L). Given the evidence of partial nitrification in DMW3 (and low concentrations), the ratio of contamination from the EMS was assumed to be best represented by DP10-2, which is located directly down-gradient of the EMS. Data for this well indicate values of predominantly ranging from 0.1 to 0.3 with Cl estimated at (Fig. 4).
The maximum concentration in groundwater at CFO1 (66.4 mg L) was measured in core sample DC15–23 (clay at 2 m b.g.l., 7 m hydraulically down-gradient of DMW3). Pore water extracted from the unsaturated zone (sand) at the top of this core profile contained 865 mg L of and had a ratio of 1.04, consistent with the ratio of 0.95 in the core sample. Given this consistency, and that concentrations in the well immediately up-gradient were low (DMW3), the in this core sample was most likely introduced into the groundwater system by vertical infiltration or diffusion from above. In contrast, elevated (up to 21.1 mg L) within the sand between 6 and 12 m depth in this core had ratios consistent with an EMS source (0.07 to 0.31). Stable isotope values in pore water from this sand layer do not indicate substantial denitrification ( 5.9 ‰, 16.7 ‰), suggesting these ratios will be similar to the initial ratios at the point of entry to the groundwater system.
In DMW13 (33 m down-gradient from DP10-2) the ratio of Cl was , similar to the ratio in DC15-23 at 2 m (0.95), which is interpreted as reflecting a top-down source. The in DMW13 is therefore unlikely to be sourced solely from leakage from the EMS, and could be sourced from the adjacent dairy pens or a temporary manure pile that was observed adjacent to this well during core collection in 2015 (or a combination of EMS and top-down sources).
In DMW12 the Cl ratio was not inconsistent with an EMS source, but the hydraulic gradient between DMW2 and DMW12 is negligible, indicating a lack of driving force for advective transport from the EMS towards DMW12. This is also the case for well DMW1, which is up-gradient of the EMS but had elevated concentrations (, ). The source of nitrate in these wells is therefore unlikely to be related to leakage from the EMS, but alternative sources (i.e. nearby temporary manure piles) are not known.
Well DMW11, 470 m from the EMS, had consistently low ratios (), similar to DP10-2, but estimates of Cl were 3 times higher than Cl for DP10-2 (Fig. 4b). and Cl estimated for DMW11 were consistent with measured values in that well, indicating a local top-down source. Well DMW11 is located hydraulically down-gradient of feedlot pens and adjacent to a solid manure storage area, in a local topographic low. Elevated in this well is therefore interpreted to be from surface runoff and top-down infiltration, rather than lateral advection from the EMS.
CFO4
At CFO4, measured data indicate that effects from agricultural operations on concentrations in groundwater are restricted to the upper 15 m of the subsurface. concentrations in wells screened at 15 m depth were mg L, with the exception of one sample from BP10-15w (May 2012) with 4.3 mg L of . Water table wells in the west and north of the study site (BC1, BC2, and BC3) also indicate negligible impacts of agricultural operations, with 10 mg L and 0.1 mg L.
Figure 5Temporal variations in (a) , (b) , and (c) at CFO4. Only wells with 10 mg L are shown.
[Figure omitted. See PDF]
Concentrations of 10 mg L were measured in three water table wells (BMW2, BMW3, BMW4) adjacent to the EMS, indicating that they have been impacted by the EMS (Fig. 5). Of these, BMW2 had much higher concentrations ( mg L, in BMW2 compared to mg L in BMW3 and mg L in BMW4), and therefore lower ratios (). concentrations in BMW2 were consistent with concentrations in the EMS suggesting direct leakage, while stable isotopes of and initial concentrations ( mg L) indicate substantial denitrification (Table 3, Fig. 6). The Cl ratio in BMW2 is consistent with measured in BMW4, which therefore likely reflects leakage from the EMS without denitrification (consistent with stable isotope of values of ).
Figure 6(a) Estimated ratios (mean and SD) in water table wells with evidence of denitrification at CFO4, plotted with distance from earthen manure storage (EMS), where dashed lines are upper and lower bounds of BMW2 (EMS source) and values are maximum measured (mg L). (b) Estimated concentrations of and at CFO4 (mid-range, error bars are max and min values).
[Figure omitted. See PDF]
Given that the estimated subsurface travel distance during operations at this site is 10 m, agriculturally derived in other wells not immediately adjacent to the EMS is unlikely to be related to leakage from the EMS. Wells BMW5 and BMW7 are 60 and 140 m hydraulically down-gradient from the EMS, respectively. Cl ratios in these wells were not inconsistent with BMW2 (i.e. the range of values overlap), but given the distance from the EMS the source of in these wells is most likely the adjacent dairy pens. Concentrations of 10 mg L were also measured in BC4, which is located 95 m hydraulically up-gradient of the EMS. The ratio of Cl at BC4 was the highest at CFO4 (0.6) and did not overlap with BMW2. The in this well is interpreted to have been sourced from an adjacent manure pile, which was observed during the study.
Mechanisms of attenuation of agriculturally derived
Attenuation of agriculturally derived in groundwater is dominated by denitrification at both CFO1 and CFO4, with estimates of consistently higher than estimates of (Tables 3 and S10, Fig. 7). Calculated values indicate that where denitrification was identified, at least half of the present at the initial point of entry to the groundwater system has been removed by this attenuation mechanism. Comparison of (the concentration of that would be measured if mixing was the only attenuation mechanism) with measured concentrations (which reflect attenuation by both mixing and denitrification) suggests that the sample from 20 m depth (DP11-12b) is the only sample that would be below the drinking water guideline if mixing was the only attenuation mechanism (Fig. 8).
Table 3Calculated and based on measured and concentrations and stable isotope values of .
Study | Sample ID | ||||||
---|---|---|---|---|---|---|---|
area | (mg L) | (mg L) | (‰) | (‰) | (mean | (mid- | |
SD) | range) | ||||||
CFO1 | DP11-13_4.3m | 28.5 | 7.0 | 30.3 | 9.8 | 0.58 | |
DP11-13_5.2m | 25.0 | 7.8 | 31.0 | 10.8 | 0.58 | ||
DP11-13_7m | 72.3 | 12.0 | 31.6 | 10.2 | 0.65 | ||
DP11-13_7.9m | 70.8 | 9.1 | 36.4 | 14.0 | 0.68 | ||
DP11-13_8.8m | 81.7 | 10.9 | 29.6 | 9.9 | 0.63 | ||
DC15-22_10m | 73.0 | 11.0 | 26.1 | 7.4 | 0.63 | ||
DP10-2 | 74.5 | 11.8 | 24.2 | 4.8 | 0.63 | ||
DMW11 | 436.1 | 17.1 | 33.3 | 10.9 | 0.83 | ||
DMW12 | 78.0 | 2.57 | 29.8 | 14.3 | 0.54 | ||
DMW13 | 56.7 | 23.7 | 23.0 | 6.8 | 0.65 | ||
DP11-12b | 95.7 | 0.6 | 35.9 | 17.0 | 0.54 | ||
CFO4 | BC4 | 163.1 | 35.1 | 30.6 | 1.6 | 0.82 | |
BMW2 | 595.6 | 16.5 | 41.6 | 8.3 | 0.92 | ||
BMW5 | 131.2 | 12.9 | 28.9 | 6.5 | 0.63 | ||
BMW6 | 156.0 | 0.4 | 70.5 | 22.1 | 0.56 | ||
BMW7 | 134.7 | 11.6 | 34.0 | 5.9 | 0.68 |
Central depth of core samples, , indicated as SampleID_xm. Maximum is 1 for all samples, which implies no mixing.
Figure 7Relative contributions to attenuation by mixing and denitrification, as indicated by estimated and at (a) CFO1 and (b) CFO4, for groundwater samples with denitrification indicated by stable isotope values of .
[Figure omitted. See PDF]
At both sites, the stable isotope values of indicate that denitrification proceeds within metres of the source. At CFO1, calculated in well DP10-2 (2 m from the EMS) is ; at CFO4, in well BMW2 (3 m from the EMS) is . Denitrification also substantially attenuated concentrations in wells where the source is not the EMS but instead is adjacent solid manure piles (e.g. DMW11 at CFO1, BC4 at CFO4). In BMW6 at CFO4, denitrification completely attenuated the agriculturally derived . This well had negligible ( mg L, ) and the lowest of 0.01. Measured DOC in this well was consistent with other wells at both sites ( mg L, ), suggesting DOC depletion does not limit denitrification at these CFOs.
Figure 8Measured concentrations of (blue circles – attenuation by mixing and denitrification) and (red triangles – attenuation by mixing only) vs. mid-range estimate of at (a) CFO1 and (b) CFO4. Dashed lines are drinking water guideline (10 mg L of ).
[Figure omitted. See PDF]
Discussion
Implications for on-farm waste management
Agriculturally derived at these two sites with varying lithology was generally restricted to depths 20 m, consistent with previous studies at CFOs (Robertson et al., 1996; Rodvang and Simpkins, 2001; Rodvang et al., 2004; Kohn et al., 2016). Attenuation of agriculturally derived in groundwater was a spatially varying combination of mixing and denitrification, with denitrification playing a greater role than mixing at both sites. In the samples for which could be determined, denitrification reduced concentrations by at least half and, in some cases, back to background concentrations. Given that the range of source isotopic composition was allowed to vary to its maximum justifiable extent, these quantitative estimates of denitrification based on stable isotopes of are likely to be conservative. Redox conditions within the groundwater system were not able to be determined in this study due to the sampling method used to collect groundwater from wells screened across low- formations (well bailed dry then sample collected after water level recovery). However, denitrification appears to proceed within metres of the source, suggesting relatively short subsurface residence times are required and that redox conditions close to the water table are conducive to denitrification reactions (Critchley et al., 2014; Clague et al., 2015).
The substantial role of denitrification within the saturated glacial sediments at these study sites indicates the potential for significant attenuation of agriculturally derived by denitrification in similar groundwater systems across the North American interior and Europe (Ernstsen et al., 2015; Zirkle et al., 2016). Denitrification in the unsaturated zone is limited by low water contents and oxic conditions, resulting in substantial stores of in vadose zones (Turkeltaub et al., 2016; Ascott et al., 2017). in water that is removed rapidly from the site is also unlikely to be substantially attenuated by denitrification due to oxic conditions and rapid transit times (Ernstsen et al., 2015). Therefore, water management focussed on reducing the effects of contamination in similar hydrogeological settings to this study should aim to maximise infiltration into the saturated zone where concentrations can be naturally attenuated, provided that local groundwater is not used for potable water supply.
At both sites there is evidence of elevated due to leakage from the EMS, but the impact appears to be limited to within metres of the EMS. This suggests that saturation within the clay lining of the EMS has limited the development of extensive secondary porosity that would allow rapid water percolation (Baram et al., 2012). Infiltration of -rich water that has passed through temporary solid manure piles and dairy pens has resulted in groundwater concentrations as high as those associated with leakage from the EMS (e.g. DMW11, BC4). At CFO4, this is in spite of the presence of clay at the surface, reflecting secondary porosity in the upper part of the profile that has led to hydraulic conductivities comparable to sand. This is consistent with the findings of Showers et al. (2008), who investigated sources of at an urbanised dairy farm in North Carolina, USA. Construction of EMS facilities in Alberta has been regulated under the Agriculture Operation Practices Act since 2002, which requires them to be lined with clay to minimise leakage (Lorenz et al., 2014). On-farm waste management should increasingly focus on minimising temporary manure piles that are in direct contact with the soil to reduce contamination associated with dairy farms and feedlots.
Critique of this approach and applicability at other sites
At both sites, leakage from the EMS had Cl of between 0.1 and 0.4, but this alone was not diagnostic of the source. The sources of manure-derived (manure piles vs. EMS) are distinguishable based on Cl ratios, provided there is also an understanding of the history of each site, local hydrogeology, and potential sources. Calculated and generally decreased with increasing subsurface residence time and distance from source, providing additional evidence for source attribution. For example, at CFO4, well BMW2, which is adjacent to the EMS, had the highest (0.92), indicating the least attenuation of by mixing and consistent with the EMS being the source of to this well. Temporal variability in Cl for each source could not be determined based on the snapshot isotope sampling conducted, but this could be investigated by measuring isotopes in conjunction with and at multiple times.
Figure 9Effect of neglecting background concentrations ( or ) in the mixing model on calculated over the range of values in this study.
[Figure omitted. See PDF]
Calculation of Cl assumed that background concentrations could be neglected in the mixing model. At these study sites, background concentrations are likely to be mg L for and mg L for . Estimated values were at least 20 times background concentrations, and over 100 times background concentrations in some wells. The estimated Cl values were at least 3 times as high as the background concentrations at CFO1 and at least 10 times as high as the background concentrations at CFO4. The error introduced by neglecting background concentrations was assessed by comparing calculated with and without background concentrations included, using the full range of values in this study (Fig. 9). Neglecting background concentrations results in overestimation of (i.e. underestimation of the amount of attenuation mixing) with the largest errors occurring when measured concentrations are close to background concentrations. For the maximum difference of 0.13 is in the mid-range of values. For , the difference is consistently with the largest errors at the lowest values of . The uncertainty in is primarily related to uncertainty in the initial concentrations (Cl and ), which depends on measured and . The largest uncertainties in and Cl correspond to the lowest measured concentrations (i.e. furthest from the upper limit), with less uncertainty at higher measured concentrations as they approach the maximum values.
Although applicable at these sites, this approach may not be valid at other sites if additional sources of in groundwater (e.g. fertiliser or nitrification) are significant, or if concentrations in groundwater are naturally elevated (Hendry et al., 1984). The combination of the approach outlined here with measurement of groundwater age indicators would allow for better constraints on groundwater flow velocities and determination of denitrification rates (Böhlke and Denver, 1995; Katz et al., 2004; McMahon et al., 2004; Clague et al., 2015).
Comparison with isotopic values of in previous studies
Nitrate isotope values in groundwater at the two CFOs studied were generally consistent with previous studies reporting denitrification of manure-derived at dairy farms (Wassenaar, 1995; Wassenaar et al., 2006; Singleton et al., 2007; McCallum et al., 2008; Baily et al., 2011). However, the isotopic values of in the manure filtrate from the EMS at CFO1 were not consistent with values for manure-sourced reported in other groundwater studies (Wassenaar, 1995; Wassenaar et al., 2006; Singleton et al., 2007; McCallum et al., 2008; Baily et al., 2011). This is likely to be because nitrification within the EMS was negligible ( 0.7 mg L), such that the isotopic values of in the manure filtrate reflect volatilisation of and partial nitrification within the EMS. values may also have been affected by evaporative enrichment of the being incorporated into (Showers et al., 2008).
A number of groundwater samples collected during this study had relatively enriched ( ‰) with depleted ( ‰). Some of these isotopic values are within the range previously reported for derived from inorganic fertiliser ( from to 3 ‰ and from to 25 ‰), with the depending on whether the is from or in the fertiliser (Mengis et al., 2001; Wassenaar et al., 2006; Xue et al., 2009). To the best of our knowledge, however, no inorganic fertilisers have been applied at these study sites. Another potential source is derived from soil organic N, but this should have values of 0 to 10 ‰ and values of to 15 ‰ (Durka et al., 1994; Mayer et al., 2001; Mengis et al., 2001; Xue et al., 2009; Baily et al., 2011). Incomplete nitrification of can result in lower than the manure source (Choi et al., 2003), but as there was no measurable in these samples this is also unlikely. These isotope values may reflect the influence of from precipitation, which usually has values ranging from to 5 ‰ for and 40 to 60 ‰ for and has been reported to dominate isotope values of groundwater under forested landscapes (Durka et al., 1994). Alternatively, they may be affected by microbial immobilisation and subsequent mineralisation and nitrification, which can mask the source in aquifers with long residence times (Mengis et al., 2001; Rivett et al., 2008).
Conclusions
A mixing model constrained by quantitative estimates of denitrification from isotopes substantially improved our understanding of nitrate contamination at these sites. This novel approach has the potential to be widely applied as a tool for monitoring and assessment of groundwater in complex agricultural settings. concentrations in excess of the drinking water guideline were measured at both sites, with sources including manure piles, pens, and the EMS. Even though these sites are dominated by clay-rich glacial sediments, the input of to groundwater from temporary manure piles and pens resulted in concentrations comparable to (or greater than) leakage from the EMS. This is attributed to the development of secondary porosity within unsaturated clays.
Nitrate attenuation at both sites is dominated by denitrification, which is evident even in wells directly adjacent to the source. In the wells for which denitrification was identified, concentrations of agriculturally derived had been reduced by at least half and, in some wells, completely. In the absence of denitrification all but one of these wells would have had concentrations above the drinking water guideline.
These results indicate that infiltration to groundwater systems in glacial sediments where can be naturally attenuated is likely to be preferable to off-farm export via runoff or drainage networks, provided that local groundwater is not a potable water source. On-farm management of manure waste at similar operations should increasingly focus on limiting manure piles that are in direct contact with the soil to limit contamination of groundwater.
Data availability
Alberta Agriculture and Forestry are the custodians of the data used in this paper.
The supplement related to this article is available online at:
Author contributions
Investigation was carried out by SB, MI and JM with assistance from staff at AAF, NRCB and USask. SB developed the data analysis methodology and prepared the paper. All co-authors contributed to supervision, conceptualization, review and editing.
Competing interests
The authors declare that they have no conflict of interest.
Acknowledgements
This research was supported by Alberta Agriculture and Forestry (AAF) and the Natural Resources Conservation Board (NRCB), who provided assistance with field work and laboratory analysis. Funding was also provided by a Natural Sciences and Engineering Research Council of Canada (NSERC) Industrial Research Chair (IRC) (184573) awarded to M. Jim Hendry. The authors thank Barry Olson at AAF for reviewing the paper. Our thanks also to the local producers, whose cooperation made this research possible. And finally, our thanks to Sebastien Lamontagne, Huaiwei Sun and anonymous reviewers for their valuable comments and suggestions during the review process. Edited by: Bill X. Hu Reviewed by: Huaiwei Sun and three anonymous referees
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2019. This work is published under https://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Leaching of nitrate (
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details
1 Department of Geological Sciences, University of Saskatchewan, Saskatchewan, SK, S7N 5C9, Canada; School of Earth Sciences, University of Western Australia, Crawley, WA, 6009, Australia
2 Natural Resources Conservation Board, Calgary, AB, T2P 0R4, Canada
3 Alberta Agriculture and Forestry, Irrigation and Farm Water Branch, Edmonton, AB, T6H 5T6, Canada
4 Department of Geological Sciences, University of Saskatchewan, Saskatchewan, SK, S7N 5C9, Canada