Subscriber access provided by UNIV OF NEW ENGLAND ARMIDALE
Article
An isotope investigation of nitrate in soils and agricultural drains of the lower Yakima Valley, Washington Dallin P. Jensen, and Carey Gazis ACS Earth Space Chem., Just Accepted Manuscript • DOI: 10.1021/ acsearthspacechem.8b00086 • Publication Date (Web): 29 Oct 2018 Downloaded from http://pubs.acs.org on November 1, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Earth and Space Chemistry
1
An isotope investigation of nitrate in soils and agricultural drains of the lower
2
Yakima Valley, Washington
3
Dallin P. Jensen
4
Carey Gazis
5
[email protected],
[email protected], Department of Geological Sciences, 312 Science
6
II, Central Washington University, Ellensburg, WA 98926
7
Keywords: Stable isotope, oxygen-17, nitrate in soils, atmospheric nitrate, semi-arid, , caliche
8
Abstract
9
Nitrate in the groundwater of the lower Yakima Valley, Washington frequently
10
exceeds the EPA maximum contaminant level standard for potable water (10 mg/L),
11
impacting communities with disadvantaged socio-economic status. Nitrogen and
12
oxygen isotopic signatures were determined for nitrate in soil leachates and
13
irrigation return flow. Isotope signatures for nitrate from soil leachate had
14
significant overlap with both the isotope signatures of naturally occurring soil
15
nitrate at the nearby Hanford Site, Washington and of groundwater nitrate
16
attributed to manure and fertilizer application in a local EPA study. A mass balance
17
calculation using Δ17O data suggests that there is a consistent ~9% atmospheric
18
contribution to nitrate in soil accumulations below caliche layers at several
19
locations. This agrees with other research on the atmospheric contribution to
20
naturally occurring soil nitrates in areas with similar mean annual precipitation
21
values. We argue that this consistent ~9% atmospheric component indicates that
22
soil nitrate at depth is dominated by naturally occurring, biologically fixed nitrate
23
across multiple sites. We suggest the flushing of naturally occurring soil nitrate to
ACS Paragon Plus Environment
ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 2 of 22
24
groundwater during land use conversion to irrigated agriculture may represent a
25
previously overlooked significant nitrate input to aquifers in this region.
26 27
Introduction The global nitrogen cycle
28
has been greatly perturbed
29
by human activities since the
30
advent of the industrial
31
revolution. Resultant
32
environmental impacts have
33
resulted in great scientific
34
effort in understanding
35
nitrogen sources at many
36
spatial scales, with stable
37
isotopes frequently
38
employed as one of the few
39
successful approaches for
40
understanding nitrate
41
source and fate1. The
42
aquifers of the lower Yakima Valley, Washington (Figure 1) have been known since
43
2002 to contain water with nitrate in excess of the Environmental Protection
44
Agency’s (EPA) Maximum Contaminant Level (MCL) for drinking water of 10 mg/L2.
45
In 2008, the Yakima Herald Republic highlighted that nitrate-contaminated domestic
46
wells are disproportionately impacting a disadvantaged community of Latino
ACS Paragon Plus Environment
Page 3 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Earth and Space Chemistry
47
farmworkers3. This spurred the EPA to investigate and largely attribute
48
groundwater contamination to local dairy farms. This attribution was based on
49
agricultural chemicals in both dairy manure lagoons and downgradient wells,3 the
50
assumption that anthropogenic inputs dominate nitrate sources to groundwater,
51
and the use of nitrate stable isotope analysis which has a well-established scientific
52
backing4. In addition, the EPA study found nitrate in a few disparately distributed
53
wells to be anomalously enriched in 18O suggesting potential atmospheric nitrate
54
inputs to some groundwater. They suggested that in this dry environment,
55
atmospherically deposited nitrate may accumulate in caliche horizons in soils along
56
with carbonates and other minerals3. This study was undertaken to explore this
57
possibility, incorporating Δ17ONO3- data for the first time in this region, as it has been
58
successfully employed in other landscapes for nitrate source attribution to constrain
59
the significance of atmospheric nitrate5,6. A separate stable isotope study at the
60
Hanford Site (80 km away) concluded that vadose soil nitrate, naturally occurring at
61
depths of greater than 20 meters, was flushed into groundwater after the discharge
62
of wastewater onto the ground surface. Our study further investigates the role of
63
soils and atmospheric inputs as a potential source of nitrate to groundwater in the
64
lower Yakima Valley by analyzing stable isotopes of nitrate from soil leachates and
65
irrigation return flow in the region.
66
Site Description
67
The lower Yakima Valley naturally experiences limited natural groundwater
68
recharge (2.5 mg/L and the agricultural drain
111
water samples were analyzed at the University of Pittsburgh’s Regional Stable
112
Isotope Laboratory to obtain δ17ONO3― , δ18ONO3― , and δ15NNO3― values (±0.1‰) using
ACS Paragon Plus Environment
ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
113
the Sigman – Casciotti bacterial denitrifier method13, and a continuous flow GV
114
Instruments IsoPrimeTM stable isotope ratio mass spectrometer. The remaining
115
leachates with nitrate concentrations between 1.0 mg/L and 2.5 mg/L were
116
analyzed using the same method at the Washington State University Stable Isotope
117
Core Laboratory for δ18ONO3― and δ15NNO3― 14 (±0.05‰) using a Gas Bench II and a
118
ThermoFinnigan Delta V Ratio Mass Spectrometer. All oxygen isotope data
119
presented are calculated with respect to Vienna Standard Mean Ocean Water
120
(VSMOW) and reported in per mil. The relationship between δ18ONO3― and δ17ONO3―
121
values during mass dependent fractionation was used to calculate Δ17ONO3-15 where:
122
Δ17O =δ17O -0.52 x δ18O (1)
123
Positive Δ17O values were used to determine the fraction of nitrate from
124
atmospheric sources (fatm) which can be estimated based on mass balance
125
considerations that assume a mixture of atmospheric nitrate and nitrate along the
126
terrestrial fractionation line using the equation16,17:
127
𝑓𝑎𝑡𝑚 =
𝛥17𝑂 +23.4‰
(2)
128
Equation 2 uses the assumption that Δ17O values will fall under a two end-member
129
mixing line where non-atmospheric nitrate will have a value of 0.0‰ and
130
atmospheric nitrate will have a value of 23.4‰.
131
Results and Discussion
ACS Paragon Plus Environment
Page 6 of 22
Page 7 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Earth and Space Chemistry
132
Nitrate Concentrations. Analytical data for nine soil sites and two irrigation return
133
flow sites are presented in Supplementary Information Table S1. For the nine soil
134
sampling locations, nitrate in soil leachate was typically detected at levels above 1
135
mg/L for at least some portion of the soil profile. Of the nine soils sites sampled,
136
three were found to have leachate nitrate concentrations above 5 mg/L and selected
137
for detailed analysis. With respect to soil pore water, the nitrate concentrations
138
measured in this study have a relative, and not an absolute relationship. It is likely
139
that actual nitrate concentrations in soil pore waters are higher because of the 10:1
140
mixing ratio between DI water and soil, and the short (10 minute) period of time in
141
which leaching occurred. Downward percolating water interacts with the soil
142
column in a natural environment for much longer time periods and the ratio of
143
water to soil is usually less than 1:1.
144
Two of these sampling sites were from road cuts, while the third was
145
obtained from a moist silt loam likely wetted by irrigation water. All of these sites
146
allowed unusually easy access for sampling throughout the soil profile. In the two
147
roadcut soils, the highest nitrate concentrations lie within the caliche horizon. This
148
spatial pattern suggests a sampling bias between road cuts and soil cores due to
149
difficulty in penetrating nitrate rich caliche layers with a soil corer. IN3 (apple
150
orchard) and IN1 (shrub steppe roadcut) had relatively low nitrate concentrations,
151
however IN4 (shrub steppe eight meters south of orchard) had nitrate
152
concentrations as high as 9.49 mg/L. This may be due to migration of irrigation
153
water downgradient from the orchard to IN4. This is likely because IN4 is two
ACS Paragon Plus Environment
ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
154
meters lower elevation than the orchard, implying a significant hydrologic gradient
155
for irrigation water.
156
Notably, soil leachates from the two orchard locations (IN3 and I1) had the
157
lowest concentrations of nitrate out of any produced during this study, including
158
natural settings, despite known fertilizer application. It is probable that much of the
159
naturally occurring or added nitrate once present in these soils has been used by
160
plants or conveyed to groundwater by irrigation.
161
δ18O and δ15N nitrate data.δ18O and δ15N nitrate data collected are plotted in
162
Figure 2, along with typical source ranges4, data from an EPA groundwater study3,
163
and the range of values found in nitrate from soil extractions at the Hanford Site18
164
(~80 km east). The EPA groundwater study largely attributed well water nitrate
165
contamination to dairy manure with several outliers indicative of a significant
166
atmospheric nitrate contribution. The values determined for soil samples in this
167
study overlapped the range of values found to occur in groundwater in the EPA
168
study3. However, overlap in typical nitrate source values make it difficult to
ACS Paragon Plus Environment
Page 8 of 22
Page 9 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
169
ACS Earth and Space Chemistry
distinguish between a natural soil
170 171 172 173 174 175 176 177 178 179 180
FIGURE 2. Chart of δ18O vs δ15N for soil leachate and agricultural drain samples. Typical nitrate isotope source ranges after Kendall et al. (2007), shaded region of natural soil pore water values for a study at the Hanford Site18 (80 km east), groundwater values plotted for an EPA study3, soil leachate from an undergraduate thesis19 and soil leachate values from this study. Two agricultural drain samples from the Marion Drain, and the Sulfur Creek Wasteway are also plotted. Arrows signify typical alteration of isotope signatures from bacterial denitrification from arbitrarily selected δ18O and δ15N values. Drains contain agricultural run-off from a large portion of the study area.
181
nitrate versus a mixture of manure and ammonium fertilizer without prior land use
182
information or other chemical tracers. Farmers commonly apply a mixture of
183
ammonium fertilizer and synthetic nitrate fertilizer20 which will tend to result in
184
higher δ18O values and lower δ15N values when mixed with either naturally
185
occurring nitrate or manure in soils. During denitrification reactions, the remaining
ACS Paragon Plus Environment
ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
186
nitrate will move along 1:1 and 2:1 trajectories on at δ18O - δ15N plot4 (Figure 2).
187
δ15N values for ammonium fertilizer may additionally experience enrichments of up
188
to 15‰ during ammonium devolatilization4, this fertilizer may then be converted to
189
nitrate.
190
Δ17O nitrate values. Δ17O nitrate values of soil leachates and agricultural drains
191
ranged from -1.2 to +2.2‰ (Figure 3). It is unknown why Δ17O values may deviate
192
negatively from the mass dependent δ 17O versus δ18O line. However, significantly
193
negative Δ17O values occurred exclusively in irrigation return flow and irrigation
194
influenced soils (IN4). Other negative values have been reported before for biogenic
195
soil nitrate in a forested catchment21.
196 197 198 199 200
FIGURE 3. Plot of 𝛅𝟏𝟖𝐎𝐍𝐎𝟑― versus 𝛅𝟏𝟕𝐎𝐍𝐎𝟑― for soils and agricultural drains sampled in the lower Yakima Valley, the trend for caliche containing soils at depth, and the Terrestrial Fractionation Line (TFL) for mass dependent fractionation.
ACS Paragon Plus Environment
Page 10 of 22
Page 11 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Earth and Space Chemistry
201 202 203 204 205 206
Figure 4. Trends in nitrate isotope ratios with interpretations (a, d, g), nitrate and carbonate concentrations (b, e, h), and percent of nitrate with atmospheric origin (c, f, i) versus depth for locations RC2, RC1 and IN4. Caliche was observed in the field where carbonate content was found to exceed five percent.
207
Interpretation.Soil samples taken from RC2 presented complex trends in δ15NNO3―
208
and δ18ONO3― (Figure 4a). The shallowest sample analyzed (15-30 cm) had values
209
typical of biologic nitrate, potentially due to biotic processing in the root zone
210
(Figure 4c).. The 30-45 cm soil interval showed anomalously high values of δ15NNO3―
211
(+12.6‰) and δ18ONO3― (+12.9‰) (Figure 4a). These high values may be explained
ACS Paragon Plus Environment
ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
212
through denitrification fractionation of biologic soil nitrate4. Δ17O values indicated a
213
small atmospheric contribution of nitrate to these soil depth intervals. The depth
214
intervals 45–60 cm and 60–75 cm yielded δ15NNO3― (+0.2 to 4.6‰) and δ18ONO3―
215
(+10.8 to 12.9‰) values suggestive of a mix of biologic and atmospheric nitrate or
216
synthetic fertilizer. This depth interval also has relatively low leachate nitrate
217
concentrations between 1.5 and 1.0 mg/L. Nitrate concentrations then increase
218
through the depth intervals 75-90 cm, 90–105 cm, and 105–120 cm, from 1.2 to 3.0
219
mg/L, while δ15NNO3― (+1.0 to 6.4‰) and δ18ONO3― (+5.3 to 6.3‰) more strongly
220
reflect biologic nitrate. Δ17O values were near zero from 15–45 cm but were 1.8–
221
2.3‰ from 90–120 cm indicating an 8 to 10% atmospheric contribution which
222
reached these depths without biologic remediation (Figure 4a). RC1 soil leachate from the shallowest depth interval (0 cm to 15 cm) has δ
223 224
15
NNO3― (8.0‰) and δ18ONO3― (5.7‰) values consistent with biologic nitrate, with a
225
nitrate concentrations of 3.2 mg/L (Figures 4c and 4d). The depth interval from 15
226
to 30 cm shows δ15NNO3― (1.8 ‰) and δ18ONO3― (7.9‰) values which may be from a
227
mixture of biologic nitrate and nitrification of ammonium in precipitation. Leachate
228
from 30 to 45 cm yielded slightly elevated δ15NNO3― (9.0‰) and δ18ONO3― (8.2‰)
229
values interpreted to be from biologic fixation and then denitrification. Δ17O values
230
for 0 to 45 cm indicated a negligible atmospheric contribution (10‰), but δ
255
15
NNO3― values (0.2 to 4.6 ‰) were too low to suggest significant denitrification had
ACS Paragon Plus Environment
ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
256
occurred to alter these isotope signatures. The Yakima Valley is known to have high
257
concentrations of ammonium nitrate pollution in winter22, which may be deposited
258
in the snow pack and subsequently flushed to the subsurface during spring snow
259
melt. Over long periods of time this process may have the potential to significantly
260
increase nitrate accumulations in the subsurface.
261
All samples from soil pit IN4 (Figures 4g-i), except the sample taken from 75–
262
90 cm, exhibited δ15NNO3- values higher than the ammonium chemical fertilizer
263
range of -10 to +44, and within the natural soil range observed at the Hanford site
264
(+3‰ to +8‰)18. It is particularly challenging to uniquely determine nitrate
265
sources in an area with abundant natural soil nitrate18 and known usage of both
266
chemical ammonia and manure fertilizer. However, the land leaser at site IN stated
267
that only chemical ammonia fertilizer has been used on the site since at least 2004.
268
Δ17O values were slightly negative, indicating no significant atmospheric nitrate was
269
present. Considering that the natural soil nitrate in RC1 and RC2 has an
270
atmospheric component of between 7 and 10 percent, the absence of this
271
atmospheric component in IN4 is interpreted as largely reflecting agricultural
272
inputs. However, higher biological cycling from increased moisture availability via
273
throughflow from a nearby irrigated orchard may be an alternate explanation for
274
the lack of positive Δ17O values. The δ15NNO3― and δ18ONO3― values for IN4 are similar
275
to those for RC1 and RC2, demonstrating the difficulty associated with
276
distinguishing natural soil nitrate from anthropogenic inputs in this area using only
277
those isotopes.
ACS Paragon Plus Environment
Page 14 of 22
Page 15 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
278
ACS Earth and Space Chemistry
The Marion Drain, and the Sulfur Creek Wasteway samples (agricultural
279
return flow) yielded δ15NNO3― values of 7.1‰ and 10.1‰ respectively and δ18ONO3―
280
values of –4.2‰ and –1.9‰ respectively (Figure 2). These values are consistent
281
with a complex mixture of nitrate produced from nitrification of ammonium in
282
fertilizer and manure, natural soil nitrate, and nitrate fertilizer. Δ17ONO3- values were
283
–0.9‰ and –1.4‰, indicating no atmospheric contribution (Figure 3).
284
We conclude that natural soil nitrate represents a potential source of nitrate
285
in groundwater upon flushing during irrigation which may lead to elevated nitrate
286
concentrations. We interpret our soil leachate data at RC1 and RC2 to reflect
287
naturally occurring soil nitrate with a potential contribution of elevated
288
atmospheric ammonium nitrate in winter time. IN4 is interpreted to have nitrate
289
from commercial fertilizer, biologic fixation and denitrification within the soil.
290
Mechanism for nitrate accumulations in high desert soils. A potential
291
mechanism for large quantities of nitrate to be transported into soils is millennial-
292
scale atmospheric deposition, followed by partial biologic processing. A recent
293
study23, has shown that nitrogen cycling in soils depends on mean annual
294
precipitation (MAP) with Δ17O values that are increasingly shifted away from the
295
atmospheric value with increasing MAP due to increased biological nitrogen cycling.
296
Wang et al. (2016) developed an empirical equation to describe this shift:
297
MAP ln ( ) 253.8 = Δ17O (3) ―0.12
ACS Paragon Plus Environment
ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
298
Based on the MAP in the lower Yakima Valley (~190 mm/yr), naturally occurring
299
soil nitrate would be expected to exhibit a Δ17O value of approximately +2.5‰. This
300
is similar to Δ17O values observed at depth at sample sites RC2 and RC1 of +1.8 to
301
+2.3‰ (Figures 2 and 4).
302
The deposition of atmospheric nitrate and the cycling of nitrogen in soils
303
may work in conjunction with the formation of biological soil crusts, symbiotic
304
communities of fungi, cyanobacteria, bryophytes, algae and lichens. Studies in the
305
cold deserts of the Colorado Plateau, southwest Utah, the Mohave Desert, and the
306
Sonora Desert have investigated nitrogen cycling in these communities using micro
307
sensors, acetylene reduction assays to measure N2 fixation rates, acetylene
308
inhibition assays to measure denitrification rates, and measurements of ammonium
309
oxidation rates24–26. The results indicate that biological soil crusts fix an order of
310
magnitude more nitrate than is denitrified, due to a lack of denitrification activity,
311
leading to a flux of nitrate to the soil below during percolation events24–26. These
312
biological soil crusts have been documented to cover between 15 and 20 percent of
313
the ground surface of the Yakima Military Training Ground to the north of the study
314
area27,28. Therefore, it is possible a similar process has occurred prior to agriculture
315
in the lower Yakima Valley.
316
Over the long term, accumulated biologic and atmospheric nitrate are likely
317
flushed to the subsurface during snow melt which dominate moisture availability in
318
this basin. Downward percolating nitrate-rich soil water may be impeded by lower
319
permeability caliche deposits until the hot dry summer months when
ACS Paragon Plus Environment
Page 16 of 22
Page 17 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Earth and Space Chemistry
320
evapotranspiration rates,and matric potential increase. Similar to other arid
321
environments, as evaporation removes water from the soil column, nitrate and
322
other solutes accumulate29 Alternatively, moisture during the brief wet period may
323
penetrate the caliche layer and trace nitrate may become redistributed and then
324
evapo-concentrated into the caliche horizon during the dry summer months, as has
325
been suggest for nitrate accumulations in the Mojave Desert30. This nitrate may not
326
be remediated due to a lack of denitrifying bacteria at depth leading to incomplete
327
nitrogen cycling, similar to in other high desert environments25. This is supported
328
by EPA groundwater data with δ18ONO3― enrichments plotting near a mixing line
329
between biological soil nitrate and atmospheric nitrate, with little apparent δ15NNO3―
330
, enrichment associated with denitrification (Figure 2).
331
More recently, as this area experienced extensive land use conversion,
332
enhancing groundwater recharge by over an order of magnitude7, it is likely that soil
333
nitrate was flushed to the groundwater. This mechanism is supported by the nitrate
334
accumulations documented in this study, δ18ONO3― , enrichment in a subset of area
335
groundwaters3, the occurrence of similar processes in other semi-arid areas which
336
have experienced land-use conversion to irrigated agriculture29, and the lack of an
337
atmospheric Δ17ONO3― observed in irrigated agriculture influenced soils.
338
Implications. Studies use elevated δ18ONO3― , characteristic of atmospheric nitrate, in
339
groundwater to help assess if natural soil nitrates in caliche may represent a
340
significant input to nitrate contaminated water3,4. However, our study has found
341
natural nitrate in soil with the potential to impact groundwater lacking this well-
ACS Paragon Plus Environment
ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
342
known elevated δ18ONO3― signature. Studies into nitrate contamination of
343
groundwater in this, and other semi-arid regions, should therefore be careful to
344
avoid assigning contamination entirely to a mixture of agricultural fertilizer and
345
manure based purely on the absence of δ18ONO3― enrichment. Other studies have also
346
found nitrate-containing soils can be a significant input in an agricultural setting
347
upon land use conversion to irrigated agriculture29.
348
As IN4 and the agricultural return drains all exhibited negative Δ17O values of
349
up to –1.4‰ (Figure 2), the use of 0‰ for non-atmospheric nitrate in a two end-
350
member mixing model may underestimate the abundance of atmospheric nitrate in
351
soils and water. It is possible atmospheric nitrate was present in these agriculturally
352
impacted samples but had its characteristic positive Δ17O anomaly obscured by
353
mixing with nitrate sources with the observed negative values. If our most negative
354
observed value of–1.4‰ is used instead of 0‰ as an end-member in our mixing
355
model, atmospheric contributions to soil nitrate in sample RC2 90 to 105 cm are as
356
high as 15%, and atmospheric nitrate is present in all soil samples for which Δ17O
357
data was collected.
358
Acknowledgements
359
Funding for this study was provided by the Geological Society of America under the
360
2016 Graduate Student Research Grants program and by a 2016 graduate
361
fellowship from the Washington State Chapter of the American Water Resources
362
Association. The authors thank Ben Harlow at Washington State University for an
363
introduction to the bacterial nitrate method, and Katie Redling of the University of
ACS Paragon Plus Environment
Page 18 of 22
Page 19 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Earth and Space Chemistry
364
Pittsburgh for Δ17ONO3– sample analysis. We thank Anne Johansen and Ashleen
365
Reddy for assistance obtaining chemical analysis with the ion chromatograph. We
366
also thank the Washington State Department of Natural Resources for granting land
367
access.
368
Supporting Information
369
Table S1. Major anion concentrations, Δ17ONO3― , δ18ONO3― and δ15NNO3― values for soil
370
samples, soil leachates, and water samples. Organic content, moisture content, and
371
carbonate are given for soil samples.
372
Literature Cited
373 374 375 376
(1)
Xia, X.; Zhang, S.; Zhang, L.; Wang, G.; Zhang, L.; Wang, J.; Li, Z. The cycle of nitrogen in river systems: sources, transformation, and flux. Environ. Sci. Process. Impacts 2018, 6 (June 2018), 857–990.
377 378
(2)
Sell, R.; Knutson, L. Quality of Groundwater in Private Wells in the lower Yakima Valley; Washington State Department of Ecology; Yakima, WA, 2002.
379 380
(3)
U.S. Environmental Protection Agency. Relation Between Nitrate in Water Wells and Potential Sources in the Lower Yakima; 2013.
381 382 383
(4)
Kendall, C.; Wankel, S. D. Tracing anthropogenic inputs of nitrogen to ecosystems. In Stable Isotopes in Ecology and Environmental Science; Michener, R. ., Lajtha, K., Eds.; Blackwell Publishing: Malden, MA, 2007; pp 375–449.
384 385 386
(5)
Dejwakh, N. R.; Meixner, T.; Michalski, G.; Mcintosh, J. Using 17O to Investigate Nitrate Sources and Sinks in a Semi-Arid Groundwater System. Environ. Sci. Technol. 2012, 46, 745–751.
387 388
(6)
Liu, T.; Wang, F.; Michalski, G.; Xia, X.; Liu, S. Using 15N, 17O and 18O To Determine Nitrate Sources in the Yellow River, China. Environmental 2013, 47, 13412–13421.
389 390 391 392
(7)
Vacarro, J. J.; Jones, M. a.; Ely, D. M.; Keys, M. E.; Olsen, T. D.; Welch, W. B.; Cox, S. E. Hydrogeologic Framework of the Yakima River Basin Aquifer System, Washington: U.S. Geological Survey Scientific Investigations Report 2009 – 5152; U.S. Geological Survey: Reston, Virginia. 2009, 106.
ACS Paragon Plus Environment
ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
393 394
(8)
Lenfesty, C. D.; Reedy, T. E. Soil Survey of Yakima County Area, Washington; U.S. Department of Agriculture: Washington D.C.,1985.
395 396 397
(9)
Taylor, S. A.; Gazis, C. A. A geochemical study of the impact of irrigation and aquifer lithology on groundwater in the Upper Yakima River Basin, Washington, USA. Environ. Earth Sci. 2014, 72 (5), 1569–1587.
398 399 400
(10) Bachmann, M. Particle tracking for selected groundwater wells in the lower Yakima River Basin, Washington: U.S. Geological Survey Scientific Investigations Report 20155149; U.S. Geological Survey: Yakima, 2015.
401 402
(11) Pfaff, J. D. Method 300.0 Determination of Inorganic Anions by Ion Chromatography; U.S. Environmental Protection Agency: Cincinatti, Ohio.
403 404
(12) Horvath, B.; Opara-Nadi, O.; Beese, F. A Simple Method for Measuring the Carbonate Content of Soils. Soil Sci. Soc. Am. J. 2005, 69, 1066–1068.
405 406 407
(13) Sigman, D. M.; Casciotti, K. L.; Andreani, M.; Barford, C.; Galanter, M.; Böhlke, J. K. A bacterial method for the nitrogen isotopic analysis of nitrate in seawater and freshwater. Anal. Chem. 2001, 73 (17), 4145–4153.
408 409
(14) Mcilvin, M. R.; Casciotti, K. L. Technical Updates to the Bacterial Method for Nitrate Isotopic Analyses. Anal. Chem. 2011, 83, 1850–1856.
410 411
(15) Thiemens, M. H. Mass-independent isotope effects in planetary atmospheres and the early solar system. Science 1999, 283 (5400), 341–345.
412 413 414
(16) Michalski, G.; Meixner, T.; Fenn, M.; Hernandez, L.; Sirulnik, A.; Allen, E.; Thiemens, M. Tracing Atmospheric Nitrate Deposition in a Complex Semiarid Ecosystem Using Δ17O. Environ. Sci. Technol. 2004, 38 (7), 2175–2181.
415 416 417
(17) Costa, a. W.; Michalski, G.; Schauer, a. J.; Alexander, B.; Steig, E. J.; Shepson, P. B. Analysis of atmospheric inputs of nitrate to a temperate forest ecosystem from ∆17O isotope ratio measurements. Geophys. Res. Lett. 2011, 38 (15), 1–6.
418 419 420 421
(18) Singleton, M. J.; Woods, K. N.; Conrad, M. E.; DePaolo, D. J.; Dresel, P. E. Tracking sources of unsaturated zone and groundwater nitrate contamination using nitrogen and oxygen stable isotopes at the Hanford Site, Washington. Environ. Sci. Technol. 2005, 39 (10), 3563–3570.
422 423
(19) Cavanaugh, G. Source(s) of High Concentrations of Atmospheric Nitrates in the Lower Yakima River Basin, Washington, Central Washington University, 2013.
424 425
(20) Mattson, N.; Leatherwood, R. Nitrogen : All Forms Are Not Equal. Greenhouse Management and Production. 2009, pp 18–20, 22–23.
426 427 428
(21) Rose, L. A.; Elliott, E. M.; Adams, M. B. Triple Nitrate Isotopes Indicate Differing Nitrate Source Contributions to Streams Across a Nitrogen Saturation Gradient. Ecosystems 2015, 18 (7), 1209–1223.
ACS Paragon Plus Environment
Page 20 of 22
Page 21 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Earth and Space Chemistry
429 430 431
(22) VanReken, T.; Jobson, T.; Lamb, B.; Liu, H.; Kaspari, S. The Yakima Air Wintertime Nitrate Study ( YAWNS ); Laboratory for Environmental Research at Washington State University: Pullman, 2014.
432 433 434
(23) Wang, F.; Ge, W.; Luo, H.; Seo, J.-H.; Michalski, G. Oxygen-17 anomaly in soil nitrate: A new precipitation proxy for desert landscapes. Earth Planet. Sci. Lett. 2016, 438 (3), 103–111.
435 436 437
(24) Strauss, S. L.; Day, T. A.; Garcia-Pichel, F. Nitrogen cycling in desert biological soil crusts across biogeographic regions in the Southwestern United States. Biogeochemistry 2012, 108 (1–3), 171–182.
438 439 440
(25) Johnson, S. L.; Neuer, S.; Garcia-Pichel, F. Export of nitrogenous compounds due to incomplete cycling within biological soil crusts of arid lands. Environ. Microbiol. 2007, 9 (3), 680–689.
441 442 443
(26) Warren, S. D. Role of biological soil crusts in desert hydrology and geomorphology: Implications for military training operations. GSA Rev. Eng. Geol. 2014, 22 (16), 177– 186.
444 445
(27) Tyler, K. J. Biological Soil Crusts: Analysis of Monitoring Techniques at the Yakima Training Center, Washington, Central Washington University, 2006.
446 447
(28) Barringer, D. . Impacts caused by military vehicular tracking on the soils and vegetation of the Yakima Training Center, Colorado State University, 1997.
448 449 450
(29) Scanlon, B. R.; Reedy, R. C.; Stonestrom, D. A.; Prudic, D. E.; Dennehy, K. F. Impact of land use and land cover change on groundwater recharge and quality in the southwestern US. Glob. Chang. Biol. 2005, 11 (10), 1577–1593.
451 452 453
(30) Lybrand, R. A.; Michalski, G.; Graham, R. C.; Parker, D. R. The geochemical associations of nitrate and naturally formed perchlorate in the Mojave Desert, California, USA. Geochim. Cosmochim. Acta 2013, 104, 136–147.
454
ACS Paragon Plus Environment
ACS Earth and Space Chemistry
δ17O (‰)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 18 17 17δ O versus δ O for 1812 19 10 20 IN4 21 8 RC2 22 RC1 6 Drains 23 24 4 25 2 26 27 0 28 29 -2 30 -4 31 -10 -5 0 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nitrate in lower Yakima Valley Soils e ich Cal
5
th dep t a +2.3‰ ils g so n i TFC tain
con
10
δ O (‰) 18
15
20
25
ACS Paragon Plus Environment
Page 22 of 22