High-frequency NO − 3 isotope ( δ 15 N , δ 18 O ) patterns in groundwater recharge reveal that short-term changes in land use and precipitation influence nitrate contamination trends

Poultry manure is the primary cause of nitrate (NO−3 ) exceedances in the transboundary Abbotsford–Sumas aquifer (ASA; Canada–USA) based on synoptic surveys two decades apart, but questions remained about seasonal and spatial aspects of agricultural nitrate fluxes to the aquifer to help better focus remediation efforts. We conducted over 700 monthly δ15N and δ18O of nitrate assays, focusing on shallow groundwater (< 5 years old) over a 5-year period to gain new insight on spatio-temporal sources and controls of groundwater nitrate contamination. NO−3 concentrations in these wells ranged from 1.3 to 99 mg N L−1 (n= 1041) with a mean of 16.2±0.4 mg N L−1. The high-frequency 15N and 18O isotope data allowed us to identify three distinctive NO−3 source patterns: (i) primarily from synthetic fertilizer, (ii) dynamic changes in nitrate due to changes in land use, and (iii) from a mix of poultry manure and fertilizer. A key finding was that the source(s) of nitrate in recharge could be quickly influenced by short-term near-field management practices and stochastic precipitation events, which ultimately impact long-term nitrate contamination trends. Overall, the isotope data affirmed a subtle decadal-scale shift in agricultural practices from manure increasingly towards fertilizer nitrate sources; nevertheless, poultry-derived N remains a predominant source of nitrate contamination. Because the aquifer does not generally support denitrification, remediation of the Abbotsford–Sumas aquifer is possible only if agricultural N sources are seriously curtailed, a difficult proposition due to longstanding high-value intensive poultry and raspberry and blueberry operations over the aquifer.

Abstract. Poultry manure is the primary cause of nitrate (NO − 3 ) exceedances in the transboundary Abbotsford-Sumas aquifer (ASA; Canada-USA) based on synoptic surveys two decades apart, but questions remained about seasonal and spatial aspects of agricultural nitrate fluxes to the aquifer to help better focus remediation efforts. We conducted over 700 monthly δ 15 N and δ 18 O of nitrate assays, focusing on shallow groundwater (< 5 years old) over a 5-year period to gain new insight on spatio-temporal sources and controls of groundwater nitrate contamination. NO − 3 concentrations in these wells ranged from 1.3 to 99 mg N L −1 (n = 1041) with a mean of 16.2±0.4 mg N L −1 . The high-frequency 15 N and 18 O isotope data allowed us to identify three distinctive NO − 3 source patterns: (i) primarily from synthetic fertilizer, (ii) dynamic changes in nitrate due to changes in land use, and (iii) from a mix of poultry manure and fertilizer.
A key finding was that the source(s) of nitrate in recharge could be quickly influenced by short-term near-field management practices and stochastic precipitation events, which ultimately impact long-term nitrate contamination trends. Overall, the isotope data affirmed a subtle decadal-scale shift in agricultural practices from manure increasingly towards fertilizer nitrate sources; nevertheless, poultry-derived N remains a predominant source of nitrate contamination. Because the aquifer does not generally support denitrification, remediation of the Abbotsford-Sumas aquifer is possible only if agricultural N sources are seriously curtailed, a difficult proposition due to longstanding high-value inten-sive poultry and raspberry and blueberry operations over the aquifer.

Introduction
The global widespread use and over-application of synthetic and manure N nutrients in agriculture have caused widespread groundwater nitrate (NO − 3 ) contamination in surficial aquifers around the world (Haslauer et al., 2005;Hamilton and Helsel, 1995;Spalding and Exner, 1993). Furthermore, with global trends towards increased agricultural intensification, threats to groundwater quality are correspondingly heightened (Vorosmarty et al., 2000;Böhlke, 2002). In agricultural settings, elevated groundwater NO − In the phreatic transboundary Abbotsford-Sumas aquifer (Canada-USA, Fig. 1), long-term nitrate contamination trends and isotopic studies have been conducted over several decades. The isotopic apportionment of NO − 3 sources in the aquifer was based on two summertime synoptic isotopic studies spanning a full decade, which revealed that poultry manure was the predominant source of NO − 3 , with localized long-term shifts towards inorganic fertilizer sources (Wassenaar, 1995;Wassenaar et al., 2006) due to changes in agricultural practices (Zebarth et al., 2015). These agricultural changes include decadal shifts from beef and dairy feed production to lower N uptake berry production, changes in irrigation methods from sprinklers to low-pressure drip fertigation, planting of cover crop for berry alley management, and a slow transition from raspberries to blueberries. One critique of the previous synoptic isotope studies was that sampling (and hence interpretations) was prejudiced to summer snapshots and thereby could be biased, especially for the numerous shallow and highly responsive water table wells spanning the aquifer and winter-biased recharge (Environment Canada, 2014). The seasonal dynamics of NO − 3 sources and fluxes and the potential for isotopic changes due to soil and unsaturated zone NO − 3 cycling were not fully evaluated and need to be considered to improve surface nutrient applications and agricultural management practices.
To address this knowledge gap, we conducted highfrequency (monthly) NO − 3 concentration and isotope analysis of the ASA over a 5-year period, with a focus on shallow wells having groundwater residence times of < 5 years as determined by 3 H-He age dating. Our aim was to determine whether high-frequency (monthly) nitrate and isotope (δ 15 N, δ 18 O) sampling improved upon previous interpretations of N sources and processes and whether important seasonal changes in the proportion of NO − 3 sources recharging to groundwater were overlooked by synoptic snapshots. Our goal was to gain improved insight on the spatio-temporal sources and controls of groundwater-nitrate dynamics and thereby to help better inform agricultural nutrient management practices and potential NO − 3 remediation efforts in the aquifer.

Study area and hydrogeologic setting
The Abbotsford-Sumas aquifer is a shallow phreatic transboundary aquifer located in southwestern British Columbia, Canada, and northwestern Washington state, USA (Fig. 1). The ASA is the most intensively studied nitratecontaminated aquifer in Canada (Zebarth et al., 1998(Zebarth et al., , 2015 and covers an area of about 200 km 2 , with approximately 40 % of the total surface area in Canada (Cox and Kahle, 1999). Our study area encompassed approximately 40 km 2 on the Canadian side of the aquifer, between the Abbotsford M. Suchy et al.: High-frequency NO − 3 isotope patterns 4269 Figure 1. Location of the Abbotsford-Sumas aquifer (ASA), southwestern BC, Canada, and northwestern Washington state, USA, along with simplified agricultural land use and sampling locations with groundwater mean residence times (MRT) of < 5 years. Arrows show the approximate groundwater flow direction. +7.3 ± 1.2 Loo et al. (2017) International Airport and the Canada-USA border (Fig. 1). Land use on the aquifer is predominantly commercial raspberry and blueberry production, mixed with intensive commercial poultry barn operations ( Fig. 1) and is < 5 % rural residential; (BC Ministry of Agriculture, unpublished data). The aquifer is typically 10-25 m thick but reaches 70 m thickness towards the southeast portion (Cox and Kahle, 1999). The aquifer comprises coarse glacio-fluvial sand and gravel with minor till and clayey silt lenses (Armstrong et al., 1965), with glacio-marine clays confining the aquifer below (Halstead, 1965). The high sand and gravel content results in a high transmittance of water, with mean hydraulic conductivities (K) of 1.6 × 10 −3 m s −1 (Chesnaux and Allen, 2007) to 9.5×10 −4 m s −1 (Cox and Kahle, 1999). The thin organicpoor soils (0-70 cm) are medium-textured aeolian deposits, moderately well to well drained, and are classified as Orthic Humo-Ferric Podzols (Luttmerding, 1980 (Zebarth et al., 2015), and water table depths typically vary widely between 2 and 20 m below surface depending on the location and season. Annual water table fluctuations average ∼ 3.6 m (Scibek and Allen, 2006). The overall flow direction in the aquifer is south (Fig. 1), southeast, and southwest at linear velocities of up to 450 m yr −1 (Liebscher et al., 1992;Cox and Kahle, 1999). The aquifer is highly vulnerable to surface-derived NO − 3 and other contamination because of (i) intensive agricultural activity; (ii) highly permeable soil, coarse sand, and gravel lithology; and (iii) high precipitation amounts in the fall and winter when nutrient uptake by crops is lowest and the NO − 3 leaching potential is greatest (Kohut et al., 1989;Liebscher et al., 1992). Elevated groundwater-nitrate concentrations exceeding drinking water guidelines are observed since the 1970s (Zebarth et al., 2015). Mitchell et al. (2003) and others (Wassenaar et al., 1995) showed vertical stratification of nitrate in groundwater was linked to agricultural practices, with the greatest nitrate concentrations (> 20 mg N L −1 ) occurring in shallow water table regions (< 10 m b.g.l., metres below ground level), while average groundwater-nitrate concentrations (< 10 mg N L −1 ) in deep wells (> 10 m screen depth below average static water level) were lower and relatively stable over time. Based on Environment and Climate Change Canada (ECCC) monitoring, the highest seasonal and temporal variations in NO − 3 are found in wells screened near the water table. Both seasonal and long-term temporal variations in groundwater-nitrate over decadal timeframes are well documented (Liebscher et al., 1992;Graham et al., 2015). The aquifer has little widespread intrinsic capacity to sustain microbial denitrification (self-remediation) because of largely aerobic conditions and the low organic content of the aquifer materials (Wassenaar, 1995), but it can occur in localized pockets around riparian discharge zones.

Sample collection and analysis
Monthly groundwater samples (n = 56 per well) were collected from 19 selected monitoring wells from September 2008 to March 2013. These wells were selected based on the following criteria: (1) groundwater having a < 5-year residence time based on 3 H-He age dating (Wassenaar et al., 2006), (2) representative spatial coverage within the monitoring network, and (3) wells representing aerobic groundwater where denitrification does not occur (Tesoriero, 2000;Wassenaar et al., 2006). These criteria helped to ensure that high-frequency nitrate and isotopic patterns stem from shortterm nitrate responses unaffected by historical or subsurface biogeochemical processes or mixing with deeper water and could therefore be more explicitly linked to contemporary landscape and agricultural activities and practices happening roughly within a 5-year timeframe.
Static water level measurements were taken prior to pumping and were reported in metres above mean sea level (m a.s.l.). Groundwater was sampled from the wells using a Grundfos ® stainless-steel submersible pump, Teflon ® lined LDPE (low-density polyethylene) tubing, and stainless-steel fittings and valves. Well water was pumped through a flowthrough cell housing a calibrated YSI ® multi-probe sonde (temperature; pH; specific conductance; oxidation reduction potential, ORP; and dissolved oxygen, DO). General chemistry and NO − 3 isotope water samples were collected after at least three well volumes were purged and the YSI ® field parameters were stabilized. All bottles were rinsed three times with sample water prior to filling. Water samples for major ion and nutrient concentrations were taken in 1 L LDPE bottles, filtered through 0.45 µm cellulose acetate membrane filters, stored at 5 • C, and analyzed within 5 days for nitrate using standard ion chromatography techniques. Nitrate concentrations were determined at the Pacific and Yukon Laboratory for Environmental Testing in North Vancouver, BC, Canada. Nitrate results are reported as mg N L −1 with a minimum detectable limit of 0.02 mg L −1 .
Nitrate and chloride concentrations were log-transformed prior to analysis to ensure normal distributions and were evaluated using principal component analysis (PCA) and factor analysis. Statistical analyses (at the 95 % confidence level), including multivariate time series analyses, were conducted using the Kruskal-Wallis methods for determining seasonality, log-normal transformations, Mann-Kendall trend analyses, and Gaussian mixture and Bayesian clustering models using WQHydro ® , ProUCL 5 ® , and XLSTAT ® (Lettenmaier, 1988;Thas et al., 1998). Seasonal Mann-Kendall trend analyses were deemed inappropriate for evaluating nitrate concentration seasonality as the repeating periods were correlated to precipitation patterns instead of calendar month and because peak nitrate concentration timings varied from year to year, resulting in a determination of nonseasonality. . These nitrate exceedances were consistent with previous observations of high nitrate concentrations in shallow wells in the aquifer (Hii et al., 1999). Previous studies reported NO − 3 concentrations exceeding the MAC in 58 %, 69 %, and 59 % of wells (Wassenaar, 1995;Zebarth et al., 1998;Wassenaar et al., 2006), respectively. The current study only had ∼ 50 % monitoring wells in common with previous investigations because early studies also sampled deeper monitoring wells containing older groundwater.
A time-series analysis showed that overall NO − 3 concentrations steadily increased from 14.1 to 18.4 mg L −1 in the targeted shallow wells over our 5-year study period, which contrasted with long-term declines observed for a wider depth variety of wells in the Canadian portion of the aquifer (Zebarth et al., 2015). Graham et al. (2015) identified several key drivers causing the short-term (intra-and inter-year) nitrate trends (increases or declines) that contrasted with the long-term (inter-decadal) declines. These key drivers were primarily stochastic rainfall patterns (wet vs. dry years, with first major rainfalls occurring in October) and short-term land use change factors. The overall increasing nitrate trend in the 19 wells could be attributed to the marked increases in NO − 3 concentrations in three of the wells occurring in the second half of our study, without which the average NO − 3 concentration in the remaining 16 wells (15.5 mg L −1 ) did not change. These nitrate increases were attributed to (i) clearing of an adjacent woodlot, and the subsequent application of large quantities of poultry manure as a soil amendment up-gradient of PC-25 and PC-35 in 2011; and (ii) a raspberry field up-gradient of US-02 that underwent a renovation cycle (described in Zebarth et al., 2015), which likely also included soil organic amendments. Wells 94Q-14, PA-25, and PA-35 did not exceed the nitrate MAC because these sites were located up-gradient of the most intensive agricultural production areas.
Almost half the 19 shallow monitoring wells (47 %) showed NO − 3 seasonality, with maximum concentrations usually occurring in the springtime. Nitrate accumulates in the soil and root zones over the summer, and a large proportion of nitrate flushing to the water table happens with the first major recharge events in the fall and winter rainy season (Kowalenko, 2000). Subsequent recharge typically has lower nitrate concentrations as the availability of dissolved soil nitrate drops. Previous evidence of NO − 3 flushing in the fall is shown by Wassenaar (1995) and Zebarth et al. (1998), when precipitation, recharge rates, and soil-NO − 3 are at their peak. Coupled with vadose zone infiltration lag times of several months, which do not vary significantly between sites (Herod et al., 2015), accordingly peak NO − 3 concentrations reaching the water table are usually observed in the springtime.
All wells have aerobic groundwater (mean = 8.9 mg O 2 L −1 ; Table S1 in the Supplement); however, two sites (ABB-03 and US-02) showed short intervals of lowered DO levels (< 1 mg L −1 ) in the winter months, coinciding with higher water tables. Chloride levels were on average 8.7 ± 3.0 mg L −1 . At six sites (91-10, 91-15, PA-25, PA-35, US-02, and US-05), NO − 3 and Cl concentrations exhibited a covariance (Pearson's R correlation coefficients > 0.5), suggesting similar sources. Three sites  exhibited an offset between NO − 3 and Cl; however, the Cl peaks usually lagged behind NO − 3 peaks by 1-3 months, which is inconsistent with the conservative properties of Cl − . This was also observed by Malekani (2012) and may be caused by differences in timing of NO − 3 and Cl availability in combination with growing season leaching under drip irrigation or variable transit times though the unsaturated zone. The remaining sites exhibited limited seasonal nitrate and chloride variability or showed no correlation between these variables, suggesting more mixing within local groundwater and possibly longer transit times in the unsaturated zone.

Nitrate N and O isotopes
Overall, the mean (±SE) groundwater nitrate δ 15 N value from the 19 study wells was +7.9 ± 0.1 ‰ (n = 717), which was consistent with positive δ 15 N values of local poultry manure solid sources (Wassenaar, 1995;Loo et al., 2017) as summarized in Table 1. Higher δ 15 N source values were observed in nitrate extracted from soil below former manure stockpiles and ranging up to +14 ‰, due to ammonia volatilization (Wassenaar, 1995). Mean nitrate δ 18 O was −1.7±0.1 ‰ (n = 717), which was typical of values derived during the nitrification of manure or synthetic fertilizers (Xue et al., 2009). Previously measured groundwater δ 18 O H 2 O in the aquifer ranged between −10 ‰ and −12 ‰ and, coupled with O derived from air (+23.5 ‰), the current nitrate δ 18 O values were comparable with earlier δ 18 O values (Wassenaar, 1995) of −1.0 ± 0.3 ‰ (n = 16) and +0.5 ± 0.8 ‰ (n = 40) a decade later (Wassenaar et al., 2006).  To further assess sources and seasonality of groundwater nitrate from these 19 shallow wells, the results were evaluated using nitrate concentrations and isotopic compositions. A Keeling plot of 1/NO − 3 vs. δ 15 N (Fig. 2a), supported by a Gaussian mixing model, suggests three main nitrate groupings with the following proportions and interpretation: (i) a historical mixing (47 % , with a mixed source of manure or soil and N fertilizer (6 %). A Bayesian VVV clustering model (which considers the volume, shape, and orientation of data points) using δ 15 N and NO − 3 suggested five possible groupings (Fig. 2b), with means shown in Table 2. These findings altogether suggest that field-scale agricultural management practices upgradient of the monitoring wells resulted in four quantifiably distinctive nitrate isotopic clusters (Table 3 -source grouping). The first two Bayesian groups were amalgamated as they suggested the same isotopic source, but with partial isotopic enrichment.
Another clustering approach, based on δ 15 N trends and seasonality in the 19 wells over the course of the study was also evaluated. In this case, sites were separated into four clusters (Table 3 -

Nitrate isotopic variations
Considering the Bayesian and Gaussian clustering approaches altogether, we separated the nitrate and isotope data into four distinctive groups (Fig. 3) based on their isotopic values (three primary groups and one subgroup), both in relation to each other and to well-known NO − 3 sources. Group 1a groundwater was impacted by synthetic fertilizer and/or organic N and showed little isotopic variability, while group 1b was similar but impacted by clear short-term spikes in δ 15 N and NO − 3 . Group 2 was dominated by poultry manure with some influence of 15 N-depleted sources, while group 3 was dominated solely by poultry manure N.
The four wells categorized into group 1a, with δ 15 N values of +3 ‰ to +8 ‰ representing 21 % of the 19 sites (PA-25, PA-35, 91-07, and US-04), had a mean δ 15 N value of +5.0 ‰. The N isotope distribution of these samples suggests they were dominated by synthetic fertilizers and natural (background) soil N sources (δ 15 N of −1.0 ‰ and +4.0 ‰, respectively). Loo et al. (2017) reported that weighted δ 15 N of fertilizer treatment leachate in the ASA is +3.2 ± 2.3 ‰. Sampling wells in this group did not exhibit large seasonal swings in NO − 3 concentration or δ 15 N values, although strong seasonality was found for NO − 3 in wells PA-25 and PA-35. These isotope data suggest a combination of annual synthetic fertilizer applications with occasional poultry manure application as a soil amendment, which is a common agricultural practice in this area, particularly with blueberry crops.
The group 1b wells were distinctive because the mean nitrate δ 15 N value was more depleted than poultry manure (+6.7 ‰), but spanned a wider δ 15 N range from +2 ‰ to +16 ‰ and represented only two wells (PC-25 and PC-35). In addition, both exhibited nitrate 18 O enrichment, coupled with increasing δ 15 N values (Fig. 3a) and NO − 3 concentrations. Well PC-25 was likely subjected to localized and temporal root zone denitrification since some δ 18 O values increased above +5 ‰; however, groundwater DO values were never below 8.8 mg L −1 , suggesting groundwater microbial denitrification processes were unlikely in the aquifer. The positive δ 15 N values coupled with elevated NO − 3 (Fig. 3b) concentrations were more likely the result of soil amendment practices whereby poultry manure is applied to fields during crop replacement cycles to augment soil carbon and nitrogen content (Zebarth et al., 2015). As previously indicated, this site may also have been affected by recent adjacent woodlot clearing and poultry manure application following planting of a new blueberry crop in 2011-2012. If the elevated δ 15 N values after January 2012 are omitted from these two wells, the mean δ 15 N value drops to +4.2 ‰, which corresponds to group 1a. Furthermore, most of the group 1a and 1b wells fall along the same groundwater flow path (Fig. 1).
Wells categorized as group 2 had a mean δ 15 N of +7.8 ‰, which corresponded to both leachate from manure-treated fields (+7.3 ± 1.2 ‰; Loo et al., 2017) and poultry manure in general (∼ +7.8 ‰). More 15 N-depleted samples were likely influenced by synthetic fertilizers or residual soil N, while 15 N-enriched samples represented temporal soil zone denitrification. Group 2 wells include 91-03, 91-15, 94Q-14, ABB-02, ABB-03, ABB-05, . Wells in this group were in the majority, representing 53 % of the sites, and as with group 1 did not exhibit large seasonal or inter-annual swings in NO − 3 concentrations or their δ 15 N values, other than NO − 3 concentrations and δ 15 N values that were elevated compared to group 1. Based on these results, it appeared that poultry manure applications or excess residual soil N from historical poultry manure applications influenced the nitrate contamination level in these wells.
The group 3 wells (91-10, US-02, and US-05) had a mean δ 15 N value of +12.6 ‰ and a δ 15 N range of +9 ‰ to +16 ‰,  (Wassenaar, 1995;Loo et al., 2017). (b) δ 15 N vs. nitrate Bayesian clustering model suggests five distinct groupings. which was more 15 N enriched than local poultry manure or manure leachates (Table 1). These 15 N-enriched results likely resulted from ammonia volatilization of the source poultry manure and temporal soil zone denitrification. Ammonia volatilization occurs in poultry manure piles and during field application. The mineralized residual ammonium can have δ 15 N values up to +25 ‰ but is dependent on pH, temperature, humidity, and other environmental factors (Kendall, 1998). Group 3 sites are all located down-gradient of current and former poultry barns or known locations of on-field poultry storage piles, which was shown by Wassenaar (1995) to result in 15 N-enriched values in soil N from +7.5 ‰ to +13.6 ‰ that are transported to the aquifer.

5-year isotopic trends
The 19 monitoring wells were evaluated based on their nitrate δ 15 N and δ 18 O isotopic trends over the 5-year study period (Table 4). The trend evaluation was conducted using Mann-Kendall (monthly data) and Seasonal Kendall (bimonthly data) non-parametric tests for detection of upward or downward trends in a time series at the p > 0.05 level of significance. For individual wells, if there was insufficient evidence to detect a trend, individual well results were grouped as being stable or variable, depending on whether the δ 15 N standard deviation was < or > 1.0 ‰, respectively. Wells exhibiting seasonality were identified as group B. The analysis showed no statistically significant temporal trend in δ 15 N during the study period; however, if results from the three nitrate "spiking" sites (US-02, PC-25, and PC-35) were removed, a statistically significant δ 15 N depletion trend was observed over the 5-year period, with δ 15 N declining from +9.1 ‰ to +7.4 ‰. This finding corresponded to the previously reported finding of a decadal-scale nitrate 15 N depletion trend in the aquifer, which was attributed to a long-term shift from manure to fertilizer use (Wassenaar et al., 2006). Four wells (91-15, ABB-02, ABB-05, and FT5-12) were classified into trend group A, where analyses did not support a significant upward or downward δ 15 N trend and ±SD ≤ 1.0 ‰ (Fig. 5a). All four wells (21 %) were from distribution group 2, where δ 15 N values were +6 ‰ to +10 ‰. Interestingly, all group A sites exhibited appreciable NO −  Table 4. Results summary with 3 H / 3 He groundwater ages in years (Wassenaar et al., 2006); average water column height (metres; midscreen depth below average static water level); isotopic distribution and trend groupings; NO 3 -N, δ 15 N, and δ 18 O values (mean, standard deviation, and confidence intervals (α = 0.05)). NA means not available.  (Table S1) were synchronous, suggesting the same source. In fact, the up-gradient field of this well had undergone a renovation cycle in the preceding months, where old raspberry plants were removed followed by a poultry manure soil amendment prior to replanting. Sites 91-10 and US-05 showed similar δ 15 N and NO − 3 fluctuations, albeit smaller in magnitude, with corresponding increases in chloride and elevated DO concentrations. Sites 91-10 and US-05 are close to each another (< 200 m apart) along a similar groundwater flow path, suggesting these variations are linked. No other sites in this group were spatially proximal. Sites PB-20, PA-25, and PA-35 exhibited varying degree of coupled δ 15 N and NO − 3 seasonality, suggesting nitrate leaching was the primary driver of NO − 3 variability. For PA-25, increasing NO − 3 concentrations with 15 N enrichment (although variable in degree) were systematically observed each winter, suggesting nitrate mobilization occurred during peak winter rainfall periods.
Six sites were identified as trend group C, with increasing δ 15 N trends . These sites were evenly distributed between distribution groups 1a (3), 1b (2), and 2 (1), suggesting one driver controlling local NO − 3 concentrations and δ 15 N values. Enriching 15 N trends (often along a flow path) are usually associated with progressive microbial denitrification; however, all sites had high DO aerobic concentrations (> 8 mg L −1 ). Sites PC-25 and PC-35, which exhibited some degree of coupled δ 15 N and δ 18 O enrichment at a 2 : 1 ratio, also showed increasing NO − 3 concentrations, suggesting heavy loading of poultry manure. Prior to the marked increase in NO − 3 and δ 15 N in the spring of 2012, PC-25 and PC-35 exhibited a significant, albeit gradual, increasing 15 N trend (Fig. 5c). This revealed a second subtle driver -the increased precipitation that occurred between 2008 and 2011 and its effect on groundwater nitrate concentrations, as shown by Graham et al. (2015). Wells 91-03, FT5-25, and US-04 did not undergo any upgradient crop replacement or soil amendments and exhibited various degrees of NO − 3 and δ 15 N seasonality, further strengthening the precipitation link as a potential driver. The δ 15 N trend could be linked to the enhanced mobilization and infiltration of 15 N-depleted organic soil N.
Group D exhibited a 15 N depletion trend (Fig. 5d); consisted of monitoring wells 94Q-14, PB-35, and ABB-03; and had a negative δ 15 N shift of 1 ‰-3 ‰ and δ 15 N values between +6 ‰ and +10 ‰ (group 2). Well 94Q-14 showed δ 15 N seasonality, but not in NO be attributed to the longer-term change in nitrogen sources used from poultry manure to synthetic fertilizers. Lastly, ABB-03 showed no significant trend in NO − 3 concentrations or in δ 18 O of nitrate; however, δ 15 N and δ 18 O were correlated, while δ 15 N and NO − 3 were inversely correlated. Furthermore, ABB-03 exhibited 3-4-month intervals of lowerthan-usual DO conditions ( ≥ 4 mg L −1 compared to other seasons; Table S1), which corresponds to periods of 15 N enrichment and decreasing NO 3 , suggesting localized denitrification, that were repeatable to various degrees on a seasonal basis but were most prominent in 2011. These findings suggest localized and temporally limited denitrification may be occurring in the soil root zone in some areas, contributing to 15 N enrichment and variability of NO − 3 concentrations. Site ABB-03 was not near Fishtrap Creek (Fig. 1), which Tesoriero (2000) and Wassenaar et al. (2006) identified as a localized riparian denitrification hot spot. Enrichment in 15 N of nitrate at these sites appeared to be from temporal drivers that could be overlooked in one-time synoptic sampling (Wassenaar, 1995).

Conclusions and outlook
This study represents an unprecedented high-frequency 5year seasonal spatio-temporal study of the water table well with over 700 nitrate isotopic assays, revealing the dynamics of nitrate recharging the transboundary Abbotsford-Sumas aquifer. The high (monthly) temporal frequency of nitrate and isotopic data aimed to address concerns that infrequent nitrate isotopic or concentration synoptic samplings of shallow groundwater overlook important factors of seasonality that may be key drivers of nitrate sources and fluxes to shallow aquifers. Indeed, our study revealed new important scientific information not previously seen in the synoptic surveys that will help managers better tackle nutrient management strategies to help reduce groundwater pollution.
Overall, and unsurprisingly, we found the predominant perennial source of nitrate to the aquifer at all spatiotemporal scales within the 5-year intensive sampling period was animal waste (poultry) sources, which was already known for decades. Nitrate concentrations in young (< 5 years old) and shallow groundwater were persistently high in nitrate, ranging from 1.3 to 99 mg-N L −1 , with a mean of 16.2 mg-N L −1 , and well in exceedance of the Canadian drinking water MAC of 10 mg-N L −1 for 76 % of the wells. The study also verified a postulated and subtle decadal-scale shift towards 15 N-depleted nitrate sources, likely reflecting systematic changes in agriculture practices from the early days of indiscriminate manure disposal towards more targeted use of synthetic fertilizers or from changes in crop types and associated nutrient practices, as evidenced by the mean δ 15 N value for nitrate of +7.9 ± 3.0 ‰ compared to +10.2 ± 4.0 ‰ in the 1990s. Synthetic fertilizer and soil N represent a comparatively higher N loading in the central portions of the ASA but are flanked on both sides by higher poultry-manure-dominated N loadings (Fig. 4). The high nitrate concentrations in contemporary recharging groundwater and lack of targeted nitrogen reduction measures means widespread nitrate contamination of the aquifer is likely to persist into the foreseeable future, and our data affirm little evidence for persistent or widespread natural attenuation of nitrate by subsurface denitrification processes at any time of the year. Nitrate remediation of the aquifer will only be possible if agricultural N sources are dramatically reduced or eliminated.
In some wells we found that localized agricultural practices (i.e. nearby organic N soil amendment) had a nearly immediate multi-year negative impact, mainly exhibited by marked increases in poultry-derived N, and lasting for several years and across seasons. This common practice at field renovation and post-woodlot clearing resulted in spatial clustering and differing short-term trends for water table nitrate and its isotopes across the aquifer (Figs. 4 and 6), further revealing that infiltrating NO − 3 and its isotopic composition can change quickly in direct response to contemporary nearfield practices. Conversely, this suggests N source cutoff as M. Suchy et al.: High-frequency NO − 3 isotope patterns 4277 Figure 6. Spatial distribution of δ 15 N trend groupings, along with agricultural land use.
a groundwater remediation effort could be similarly as effective. Despite 53 % of shallow wells showing no isotopic trends, 47 % showed an 15 N enrichment or depletion trend in nitrate, and about half of the wells exhibited nitrate seasonality in NO − 3 concentrations and/or δ 15 N values controlled by temporal infiltration of residual mineralized N or weak, short-term denitrification.
Due to the rapid shift in NO − 3 and isotopic values of recharging groundwater immediately following field renovation and soil amendment practices, this study reinforces the importance of designing and conducting appropriate spatiotemporal nitrate sampling to reduce the risk of misinterpreting nitrate concentrations and its isotopic data through the more common practice of occasional synoptic surveys. The dynamics of nitrate in younger (< 5 years old) water table wells, however, also imply it would be prudent to monitor deeper, older groundwater which smoothes out short-term fluctuations and hence record longer-term and aquifer-wide trends.
For the ASA agricultural area specifically, measuring the impact of changes in nutrient management practices associated with the switch from raspberry to blueberry crops or field renovation is required to determine its impacts on groundwater nitrate dynamics. Decisions on future groundwater aquifer nitrate management need to take into consideration (i) permanent or cyclical changes in the planned crop types and (ii) the associated nutrient management practices involved with them. Subtle shifts in nitrate in the ASA may be unexpectedly influenced by the recent replacement of raspberries with blueberries, which may to be less reliant of cyclical poultry manure soil amendments.
Data availability. All data used in this study are available in the Supplement.
Author contributions. LIW and BZ prepared the funding proposal. LIW, GG, and MS established the site-selection criteria and selected sites, while MS undertook field activities. MS and LIW analyzed data statistically and conducted the data interpretation. MS prepared the manuscript with the help of LIW, and all co-authors were asked to review the manuscript.