JEQ Journal of Natural Resources and Life Sciences Education
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Published online 2 February 2006
Published in J Environ Qual 35:495-504 (2006)
DOI: 10.2134/jeq2005.0012
© 2006 American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via ISI Web of Science (2)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Panno, S. V.
Right arrow Articles by Hwang, H.-H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Panno, S. V.
Right arrow Articles by Hwang, H.-H.
Agricola
Right arrow Articles by Panno, S. V.
Right arrow Articles by Hwang, H.-H.
Related Collections
Right arrow Surface Water Quality
Right arrow Nutrients
Right arrow Nutrient Cycling
Right arrow Water Pollution
Right arrow Nitrogen

TECHNICAL REPORTS

Surface Water Quality

Isotopic Evidence of Nitrate Sources and Denitrification in the Mississippi River, Illinois

Samuel V. Pannoa,*, Keith C. Hackleya, Walton R. Kellyb and Hue-Hwa Hwanga

a Illinois State Geological Survey, Natural Resources Building, 615 E. Peabody Street, Champaign, IL 61820
b Illinois State Water Survey, 2204 Griffith Drive, Champaign, IL 61820-7495

* Corresponding author (panno{at}isgs.uiuc.edu)

Received for publication January 14, 2005.

    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Anthropogenic nitrate (NO3) within the Mississippi–Atchafalaya River basin and discharge to the Gulf of Mexico has been linked to serious environmental problems. The sources of this NO3 have been estimated by others using mass balance methods; however, there is considerable uncertainty in these estimates. Part of the uncertainty is the degree of denitrification that the NO3 has undergone. The isotopic composition of NO3 in the Mississippi River adjacent to Illinois and tile drain (subsurface drain) discharge in agricultural areas of east-central Illinois was examined using N and O isotopes to help identify the major sources of NO3 and assess the degree of denitrification in the samples. The isotopic evidence suggests that most of the NO3 in the river is primarily derived from synthetic fertilizers and soil organic N, which is consistent with published estimates of N inputs to the Mississippi River. The 1:2 relationship between {delta}18O and {delta}15N also indicate that, depending on sample location and season, NO3 in the river and tile drains has undergone significant denitrification, ranging from about 0 to 55%. The majority of the denitrification appears to have occurred before discharge into the Mississippi River.

Abbreviations: {delta}15N, isotopic composition of nitrogen • {delta}18O, isotopic composition of oxygen • {varepsilon}, isotope enrichment factor • MARB, Mississippi–Atchafalaya River Basin • OM, organic matter


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
EXCESS N, primarily as NO3, is a major pollutant of rivers, lakes, and estuaries and can adversely affect aquatic life (Parry, 1998). Nitrogen has been linked to the hypoxic zone that develops annually in the bottom waters of the Gulf of Mexico, threatening aquatic life and commercial fisheries (Goolsby et al., 1999; Rabalais et al., 2001, 2002). The N load in the Mississippi–Atchafalaya River Basin (MARB) doubled from the 1960s to the 1980s, and the mean NO3–N concentration in the Mississippi River also doubled during that period (Goolsby et al., 1999). This increase has been attributed primarily to the use of synthetic fertilizers in the Corn Belt states of the midwestern USA (Rabalais et al., 2001). Goolsby et al. (1999) estimated that ~60% of the N inputs into the MARB are from synthetic fertilizer and mineralized soil N. The remaining 40% comes from legumes (N2 fixation), manure, atmospheric deposition, and point sources such as effluent from sewage treatment facilities. Policymakers are using estimated contributions of these N sources, typically calculated by mass balance approaches, to develop strategies to reduce NO3 inputs to the Mississippi River in the hope of decreasing the hypoxic zone. There is, however, a great deal of uncertainty in these estimates and the link between fertilizer use and hypoxia in the Gulf of Mexico is poorly understood. When comparing the inputs and outputs of the N mass-balance calculations for various river systems, there is typically a significant deficit in the output portion, most of which is generally attributed to denitrification (David and Gentry, 2000).

The N isotopic ratio has been used to evaluate NO3 sources and geochemical processes affecting the concentration of NO3 (Létolle, 1980, Hübner, 1986, Heaton, 1986, Komor and Anderson, 1993). For example, NO3 originating from sewage and livestock effluent is typically more enriched in the heavier N isotope (15N) compared with N and NO3 derived from atmosphere deposition, fertilizer N, and soil organic N (Kendall, 1998). There is significant overlap, however, of {delta}15N values for NO3 derived from the latter three sources. When attempts have been made to distinguish NO3 derived from such sources using only {delta}15N values, many uncertainties arise because of the initial overlap in {delta}15N values and the isotopic changes that occur due to fractionation effects caused by geochemical processes such as denitrification (Kohl et al., 1971; Hauck et al., 1972). In recent years, researchers have included the measurement of O isotope ratios of NO3 as well as N isotopes, producing more definitive results. A number of researchers have used this dual-isotope technique to identify dominant sources of NO3 in natural waters and to determine whether or not denitrification has occurred (Heaton, 1986; Böttcher et al., 1990; Aravena et al., 1993; Wassenaar, 1995; Aravena and Roberston, 1998; Panno et al., 2001; Burns and Kendall, 2002; Chang et al., 2002).

Illinois is a major source of N to the Mississippi River. Most of the N from Illinois is agricultural in origin, but sewage is also significant, primarily from the Chicago area. David and Gentry (2000) estimated that sewage contributes about 16% of the total N to rivers in Illinois. Battaglin et al. (2001) used the dual-isotope approach to, in part, assess the degree of denitrification in the Mississippi River at several different sites from Illinois to Louisiana, but drew no conclusions as to the dominant sources of NO3. Chang et al. (2002) were more successful by comparing land use with isotopic results; they observed distinct, although overlapping, isotopic signatures. The objectives of our investigation were to identify the major sources contributing to the total dissolved NO3 in the upper Mississippi River using N and O isotopes from the NO3 ion together with published mass balance data and to assess the occurrence and degree of denitrification.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Sample Collection and Analysis
Water samples from the Mississippi River were collected quarterly from August 2001 through May 2002 near the intakes for local water utilities at Davenport, IA, Quincy, IL, and Chester, IL (Fig. 1). These utilities pump their water through pipes that either extend ~10 m from the shore or from a concrete sea wall that extends to a depth of 10 m. No major tributaries enter the Mississippi River between the Illinois–Wisconsin border and Davenport. Between Davenport and Quincy, the Rock, Iowa, and Des Moines rivers discharge into the Mississippi. Between Quincy and Chester, two of the largest tributaries in the Midwest, the Illinois and Missouri rivers, enter the Mississippi River. All of these tributaries drain land in the Corn Belt region of the midwestern USA. Because there were no major tributaries <60 km upstream of any of the sampling sites, we assumed that the concentrations of dissolved constituents were fairly well mixed and uniform across the river at each sampling site. This is supported by Meade et al. (1995), who compared single-point pumped samples collected from the center of the river with composite sampling across the width and depth of the Mississippi River and reported that the dissolved composition of the single pumped samples was comparable to the composite river water samples. To help test this assumption, during one sampling event we took samples at three locations across the river at Quincy and compared results with a sample taken at the municipal utility plant intake.


Figure 1
View larger version (31K):
[in this window]
[in a new window]
 
Fig. 1. Sampling locations. The light gray area on the U.S. map shows the extent of the Mississippi River watershed.

 
The dominant sources of N from the tile drains (subsurface drains) are synthetic fertilizers and soil organic matter (OM) from the cultivated fields. A large portion of agricultural land in the Midwest is tile drained, approximately 35% in Illinois and 25% in Iowa (Fausey et al., 1995), and these tiles contribute a large percentage of water and N to streams and rivers. Our tile water samples were collected seasonally from two tile drains in Champaign County in east-central Illinois (11 occasions between April 2000 and September 2004). Land drained at both locations is dominated by row-crop agriculture (>90%), and the tiles are typically emplaced at depths of 1 to 2 m. Nitrogen fertilizer is predominantly anhydrous NH4; it is highly unlikely that manure was used in any part of the tiled areas during our study because manure is applied to fields in close proximity to livestock-raising facilities, and there are none in the vicinity of the drained areas. Based on infrared images of the tile-drain areas, about 1.32 km2 of cropland is drained by Tile 1 and about 1.82 km2 is drained by Tile 2 (S. Walker, personal communication, 2005). At the time of sampling, the crops grown on the land drained by Tile 1 included both corn (Zea mays .) and soybean [Glycine max (L.) Merr.]. Site 2 was exclusively in soybean at the time of the first sampling (September 2003), followed by the application of N fertilizer in the fall of 2003 for corn in 2004 (the next three sampling times).

Measurements of temperature, pH, Eh (redox potential), and specific conductance (SpC) were taken in the field during sample collection. Water samples were collected using a peristaltic pump, and were filtered through in-line 0.45-µm filter capsules and stored in polyethylene bottles. The tile water samples were collected directly from the stream of water exiting the tile drains or in a clean 19-L bucket, and then pumped through an in-line filter. Filtered samples were analyzed for anions, dissolved organic carbon (DOC), NH4–N, dissolved Kjeldahl N (DKN), and NO3 isotopes. Samples collected for the analysis of {delta}15N and {delta}18O of the NO3 ion were acidified with HCl to a pH <2. All samples were transported in ice-filled coolers to the laboratory, and kept refrigerated at approximately 4°C until analysis.

Anion concentrations, including NO3–N, were determined using ion chromatography. The DOC was determined by difference (total carbon – inorganic carbon) using persulfate–oxidation and infrared detection (Greenburg et al., 1987). The DKN was determined by converting organic N to NH4+ by digestion with H2SO4; the resultant solution was analyzed for NH4–N using a titrimetric procedure (Bremner and Mulvaney, 1982). Ammonium was determined using semi-automated colorimetry (R.A. Cahill, unpublished data, 1985). Analytical precision for the above methods is as follows: ion chromatography (±0.02 mg/L for Cl; ±0.02 mg/L for NO3–N; ±0.02 mg/L for SO4; ±0.23 mg/L for alkalinity), DOC (±0.02 mg/L), DKN (±0.25 mg/L), NH4–N (±0.10 mg/L). These data are based on 2{sigma} of the repeated analysis of a standard (n = 10).

Isotopic analyses of N and O were conducted using published methods (Silva et al., 1994; 2000; Wassenaar, 1995) with some modifications (Hwang et al., 1999; Panno et al., 2001). The samples were first boiled under acidification to remove HCO3 and dissolved CO2. The dissolved organic matter (DOM) and SO42– were also removed to minimize contamination of the {delta}18O due to the O in the SO42– ions and the DOM, and to help eliminate anion interference during the ion-exchange step for NO3 extraction. The DOM was removed using a silicalite molecular sieve, and the SO42– was removed by precipitation as BaSO4. After removal of HCO3, SO42–, and DOM, the NO3 was extracted using an anion-exchange column packed with BioRad AG 1-X8 resin (BioRad, Hercules, CA). Nitrate collected on the anion-exchange column was eluted with HBr solution (Hwang et al., 1999) and converted to AgNO3 by adding AgO. The AgNO3 was precipitated by freeze-drying the sample in a vacuum system. The dried AgNO3 was converted to N2 and CO2 for {delta}15N and {delta}18O analysis, respectively, by combustion techniques (Hwang et al., 1999; Silva et al., 2000). Both gases were analyzed on an isotope ratio mass spectrometer. International isotope standards IAEA-N1, IAEA-N2, USGS25, and USGS26 were used for {delta}15N calibration. International standard IAEA-N3 was used for {delta}18O calibration. Reproducibility of duplicate analyses was ≤0.9{per thousand} for {delta}18O (averaging ±0.3{per thousand}), and ≤0.3{per thousand} for {delta}15N (averaging ±0.1{per thousand}).

Recent work by Révész and Böhlke (2002) showed that the methods used to determine {delta}18O values from the combustion of NO3 in sealed quartz glass tubes resulted in the interaction of the quartz O with the resultant CO2 from the NO3 combustion, altering the {delta}18O measured value of the NO3. The new high-temperature conversion technique discussed by Révész and Böhlke (2002) does not contain any quartz in the reaction chamber and results in a more accurate {delta}18O value for the NO3. Because most of the published {delta}18O data for NO3 up to this point were produced using quartz glass sealed-tube combustion techniques, however, we have used the sealed-tube combustion {delta}18O values for this investigation. We have analyzed three National Institute of Standards and Technology (NIST) NO3 standards recently prepared by the USGS for {delta}18O analyses and calculated a correction for the sealed-tube {delta}18O values to adjust for the quartz O influence so that the values can be compared with the new high-temperature conversion technique in the future (Table 1).


View this table:
[in this window]
[in a new window]
 
Table 1. Results from comparison of National Institute of Standards and Technology (NIST) standard {delta}18O values with our replicate values and adjusted values based on a least squares equation.

 
Stable Isotope Values for Sources of Nitrate
The stable isotopes of N and O are 14N (earth abundance of 99.64%), 15N (0.36%), 16O (99.763%), 17O (0.0375%), and 18O (0.1995%) (Hoefs, 1987). Stable isotope concentrations are expressed in the {delta} notation, which is parts per thousand difference in the ratio of the less abundant isotope to the most abundant isotope relative to the same ratio in a known standard. This may be represented by:

Formula 1[1]
where X is the isotope of interest (e.g., 18O or 15N), and R is the ratio of 18O/16O or 15N/14N, respectively. Reference standards used in Eq. [1] are the Vienna-Standard Mean Ocean Water for O and air for N. The isotopic enrichment factor ({varepsilon}) for a chemical reaction can be defined as:

Formula 2[2]
where {alpha} = Rp/Rs, Rp is the isotopic ratio of the reaction product, and Rs is the isotopic ratio of the reaction substrate (reactant). The Rayleigh equation describes the evolution of the isotopic composition of the residual substrate during certain reactions where the products are removed from chemical or isotopic equilibrium with the reactant. The following is a simplified version of the Rayleigh equation:

Formula 3[3]
where {delta}o is the initial composition of the substrate and f is the remaining fraction of the substrate (Clark and Fritz, 1997).

Typical ranges of {delta}15N and {delta}18O for various natural and anthropogenic sources of NO3 are defined by boxes in Fig. 2 (Kendall et al., 1995; Clark and Fritz, 1997; Mengis et al., 2001). Soil OM has a large range of {delta}15N values, although most soils have values between 2 and 10{per thousand} (Fogg et al., 1998; Kendall, 1998).


Figure 2
View larger version (25K):
[in this window]
[in a new window]
 
Fig. 2. Ranges of {delta}18O and {delta}15N values for potential NO3 sources and values measured in Mississippi River and tile drain samples. The domain of soil organic matter is shaded to show the range of greatest (darkest) to lowest (lightest) frequency. The domain of soil organic matter partially overlaps that of reduced N fertilizer, manure, and sewage with respect to {delta}15N. The {delta}18O of NO3 from these reduced N sources was calculated as described in the discussion. The large arrow represents a denitrification vector; as denitrification proceeds, the {delta}15N and {delta}18O values of the remaining NO3 progressively increase in the direction of the vector. The size of circles represents NO3–N concentration; shaded circles indicate tile drain samples. In all cases, analytical errors are smaller than the circle size. Synthetic NO3 Fertilizer refers to synthetic fertilizer applied as NO3 (e.g., the NO3 in NH4NO3) and Reduced N Fertilizer refers to NO3 produced by nitrification of applied NH3 (as NH3 or the NH4+ in NH4NO3) and urea.

 
The range of {delta}18O values for NO3 originating from the reduced N sources including soil OM and sewage was calculated as described below.

Most of the reduced N applied to fields is quickly oxidized to NO3, which is taken up by plants, leached from the soil zone, or denitrified to N2. A portion of the N from fertilizers can be lost via volatilization or be immobilized in the soil by uptake, storage, and recycling in the microbial biomass within the soil zone (Burkart and James, 1999; Mengis et al., 2001). Some of these processes can greatly affect the isotopic composition of the NO3 formed in the soil zone. For example, fertilizer N added as NO3 and taken up by plants or microbes, converted to organic N, and then returned to the soil and reoxidized to NO3 would no longer have the {delta}18O value of the initial applied NO3, but instead would have a {delta}18O value like that of other reduced N sources that were nitrified in the soil zone (Mengis et al., 2001). A sample so affected would be expected to plot in the Soil OM domain in Fig. 2. Thus, NO3 in ground or surface water whose original source was synthetic fertilizer may have different isotopic signatures depending on how long it remained in the soil–plant environment; i.e., fertilizer NO3 leached soon after application is isotopically different than NO3 from fertilizer taken up by plants and later returned to the soil.

Biologically mediated reactions such as denitrification cause the N and O isotopes to fractionate, increasing the {delta}15N and {delta}18O values of the remaining (undenitrified) NO3. The fractionation factors for {delta}15N and {delta}18O vary depending on the local conditions and rates of reaction; however, the ratio of the change in {delta}18O and {delta}15N during denitrification is typically close to 1:2 (Böttcher et al., 1990; Aravena and Robertson, 1998; Kendall, 1998; Mengis et al., 1999), as shown by the vector in Fig. 2.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
The NO3–N concentrations in the 2001–2002 Mississippi River samples ranged from 0.2 to 2.9 mg/L, with mean concentrations for Davenport, Quincy, and Chester being 1.2, 1.7, and 2.1 mg/L, respectively (Table 2). These concentrations are similar to mean values measured between 1980 and 1996 at stations in the Mississippi River near our sites (Goolsby et al., 1999). Nitrate-N concentrations in the tile drain samples ranged from 6.5 to 15.3 mg/L, with a mean value of 11.5 mg/L (Table 2).


View this table:
[in this window]
[in a new window]
 
Table 2. Water chemistry data.

 
Nitrate-N concentrations increased downstream in the Mississippi River at each sampling event except between Quincy and Chester in the spring (May 2002). At that time, the river was about 3.6 m above flood stage at Chester, due to flooding of the Illinois, Missouri, and Kaskaskia rivers; dilution from these rivers presumably caused the decrease in the NO3–N concentrations. The increasing NO3–N concentrations from north to south along the Mississippi River adjacent to Illinois are probably due to greater NO3–N concentrations in the tributaries draining watersheds in Iowa, Illinois, and Missouri. The mean NO3–N concentrations for the Rock, Iowa, Des Moines, and Illinois rivers from 1980 to 1996 were 3.5 to 5.0 mg/L (Goolsby et al., 1999). In contrast, the mean NO3–N concentration in the Missouri River during that period was 1.2 mg/L, about equal to the NO3–N in the Mississippi River at Davenport, IA. Thus, most of the increase in NO3–N concentration in the Mississippi River between Davenport and Chester is from northern and central Illinois and eastern Iowa.

The {delta}15N values of the NO3 in the water samples from the Mississippi River and the tile drains range was 4.8 to 16.4{per thousand} and the {delta}18O values range was 6.6 to 13.9{per thousand} (Table 2); the data plot within a narrow band on a {delta}15N–{delta}18O plot (Fig. 2). The samples having the greatest NO3–N concentrations tended to have the lowest {delta}15N and {delta}18O values. Plots of NO3–N concentrations vs. {delta}15N and {delta}18O values for the river samples show a significant negative correlation (r2 = 0.64 and 0.61, respectively), and the isotopic values tended to decrease downstream (Fig. 3). The isotopic values of NO3 were lowest and showed less variability among the samples collected during the winter and spring, when NO3–N concentrations were greatest. The tile drain samples had a similar isotopic range to the river samples, and negative correlations between NO3–N and the isotopic values, although the correlations were not as strong as that observed for the Mississippi River samples (Fig. 3). Nitrate-N concentrations also tended to be greater in the spring and winter for the tile drains.


Figure 3
View larger version (16K):
[in this window]
[in a new window]
 
Fig. 3. Nitrate-N concentrations vs. (a) {delta}15N and (b) {delta}18O of NO3 for Mississippi River and tile drain samples. Letters indicate sample locations (D = Davenport, Q = Quincy, C = Chester). The gray hexagon represents the average composition of the four river samples collected in April 2004 at Quincy. The month of sample collection is shown for tile samples. Lines are linear regressions for river samples (r2 = 0.64 for {delta}15N and 0.61 {delta}18O).

 
For the samples collected from the Mississippi River transect at Quincy and at the municipal water plant intake (Table 2) there were no significant differences, within analytical error, in NO3 isotopic compositions except for the sample collected near the west bank. The sample from near the west bank was collected in a shallow, slow-flowing part of the river at this site and obviously was not representative of the majority of the water, as indicated by the chemical analyses. It appears, however, that the sampling technique of using municipal water supply intakes provided NO3 isotopic compositions that were representative of the bulk of the river water. Because we conducted our transect sampling during a relatively high flow period, we cannot rule out the possibility that, during low flow periods, there may be more variability. It is likely, however, that the greatest variability would be found during high flow because of input from local runoff.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Sources of Nitrate to the Mississippi River
The isotopic data plot outside of the estimated domains representing the range of isotopic values expected for the major NO3 sources. The distribution of data could be explained by denitrification, mixing of different NO3 sources, or both. Plots of {delta}15N and {delta}18O vs. (NO3–N)–1 and ln(NO3–N), which have been used by some researchers to distinguish between mixing and fractionation processes, were inconclusive. Linear relationships (r2 > 0.5) were found for the river samples in both types of plots.

If mixing between atmospheric NO3 and NO3 from reduced sources was the dominant process affecting the isotopic values, the fact that the {delta}18O values of the samples are significantly greater than the values estimated for the reduced N sources would require a large amount of the NO3 to reach the river as precipitation, without passing through the soil zone. Using a {delta}18O value of 45{per thousand} for NO3 in precipitation (Kendall, 1998), we calculate this amount to be at least 10% of the total NO3, clearly an unrealistic value. Similarly, samples with relatively elevated {delta}15N values would require a large input of manure and sewage if mixing was the only process affecting isotopic values. More than 50% of the NO3 would have to come from manure and sewage sources for the {delta}15N value to exceed 12{per thousand}. If we expanded the {delta}18O range of the source domains to reach 14 or 16{per thousand}, as observed in some samples of NO3 in forest soils (Mayer et al., 2001; Burns and Kendall, 2002), and then used mixing to explain the isotopic results for the Mississippi River and tile drain samples, it would suggest that the samples with the most positive isotopic values were generated from nearly 100% sewage or manure NO3 inputs. These scenarios are also obviously unrealistic; for comparison, the Illinois River, which receives discharge from some of the largest wastewater treatment plants in the world, has been estimated to receive 21% of its N load from sewage based on mass balance approaches (David and Gentry, 2000). It is clear, however, that some of the N in the Mississippi River comes from sewage sources, especially downstream of the mouth of the Illinois River.

Because mixing of sources alone cannot account for the range of isotopic values for the samples, it appears that denitrification has played an important role. If the NO3 is primarily from one source, then following the denitrification vector back toward its intersection with a source domain would indicate the dominant source. Doing this with the Mississippi River and tile data suggests that the NO3 in the samples comes primarily from fertilizer N and soil OM (see dotted line in Fig. 2). The vector shown in Fig. 2 was passed through the center of the observed data. This vector suggests a {delta}18O value for the fertilizer and soil OM sources that is in the heavy part of these source domains for which the {delta}18O composition was calculated. These slightly greater {delta}18O values may be due to the contribution of NO3 from atmospheric and synthetic NO3 fertilizer sources. Approximately 15% of the synthetic N fertilizer applied to Illinois crop land is synthetic NO3 (USDA, 1999). To account for the possible input of higher {delta}18O values that are associated with synthetic NO3 fertilizers, typically 19 to 25{per thousand} (Mayer et al., 2001), we added a 15% contribution of NO3 with a {delta}18O value of 22{per thousand} (measured from a sample of NH4NO3 fertilizer obtained in southwestern Illinois) to estimate the range of values that would be expected from a combination of NO3 from nitrification of reduced N fertilizer and synthetic NO3 (15% NO3 + 85% Reduced Fertilizer in Fig. 2). However, we would normally not expect to observe the impact of the greater {delta}18O values from the synthetic NO3 because of the rapid uptake and transformations of NO3 to organic N within the soil zone followed by the eventual oxidation back to NO3, resetting the {delta}18O value (Mengis et al., 2001). We could expect that the isotopically heavier O of synthetic NO3 may impact the {delta}18O of the overall NO3 during periods of slow biological activity or rapid runoff of the NO3 from the fields to the streams, as observed by Mengis et al. (2001). Additionally, in some forested areas 18O-enriched atmospheric NO3 was observed to impact the {delta}18O of NO3 in the shallow soil water and stream runoff (Durka et al., 1994; Burns and Kendall, 2002), indicating that NO3 with isotopically heavy {delta}18O can, under certain conditions, affect the {delta}18O value of the NO3 in shallow soil water and runoff.

Because NO3 is not initially present in the reduced N sources (reduced fertilizer, soil OM, manure, and sewage), the {delta}18O values for NO3 that are produced by the oxidation of reduced N must be estimated. This was done by assuming that the biologically mediated nitrification process for most of the soils in the midwestern USA derives one-third of the O for the NO3 from atmospheric O and two-thirds from the surrounding water (Amberger and Schmidt, 1987; Böttcher et al., 1990; Kendall, 1998). Recent studies suggest that this is probably an oversimplification for predicting the {delta}18O of NO3 during nitrification (Mayer et al., 2001; Burns and Kendall, 2002). These studies found that, in some cases, the {delta}18O is greater than expected based on incubation experiments and on measurements of NO3 from forest soils, suggesting that perhaps the contribution of atmospheric O is sometimes greater than one-third. Mayer et al. (2001) suggested that in low-pH environments, which they observed in forest soils, a different bacterial process dominates the nitrification reaction and utilizes a greater amount of atmospheric O. Soil pHs in the forest studied by Mayer et al. (2001) were 3.2 to 3.3, whereas in Illinois, agricultural soil pHs are closer to 6.5 and range from 5.2 to 7.5 for the state (Illinois State Water Survey, 2005). There is also the possibility that the {delta}18O of the soil water may be enriched in 18O due to respiratory isotope fractionation or evaporative effects that may occur within the soil zone (Kendall, 1998; Burns and Kendall, 2002).

In agricultural settings, {delta}18O values for NO3 in shallow ground water samples are commonly quite low, with many of the larger values attributed to denitrification processes or fast throughput of NO3 fertilizers (Mengis et al., 2001; Aravena et al., 1993; Beaumont, 2003). The lower {delta}18O values in these studies are consistent with the two-thirds water, one-third atmospheric O contribution for nitrification. There is always a range of values, however, which probably reflects the variable {delta}18O of the soil water as well as the influence of some of the processes discussed above. It appears the environmental setting and climatic conditions affect the initial {delta}18O value of NO3 during nitrification of reduced N sources.

Because our samples are primarily from agricultural settings, we decided to use the two-thirds water and one third atmospheric contribution of O to estimate the values for the initial {delta}18O of NO3 for nitrification of reduced N sources. The range of {delta}18O of NO3 from the nitrification reaction will vary depending on the {delta}18O of the shallow ground water in contact with the nitrifying bacteria. Because the isotopic composition of shallow ground water reflects the average isotopic composition of precipitation for a geographic region, we used the average range of {delta}18OH2O for precipitation in Illinois and bordering states, which generally ranges from about –4 to –11{per thousand} (Coplen and Kendall, 2000; summarized in Clark and Fritz, 1997). Using these {delta}18OH2O values gives an expected range of {delta}18ONO3 of ~5.2 to 0.5{per thousand} for the reduced N sources of NO3 (Fig. 2).

Denitrification
While the dominant source(s) of the NO3 for the tile drains and the Mississippi River appears to be the same, the NO3–N concentrations were much lower for the river samples than those from the tile drains. The fact that the isotopic values for the tile and river samples cover a similar range (Fig. 2) is indirect evidence that in-stream denitrification was probably minimal in the Mississippi River. Battaglin et al. (2001) also concluded that there was little or no in-stream denitrification in the Mississippi River in their study of the MARB. The lack of evidence of denitrification in the Mississippi River was not unexpected, as N loss in streams declines rapidly with increasing channel size (Alexander et al., 2000). Dilution and biological uptake of NO3 are probably the dominant mechanisms for the lower NO3–N concentrations in the river samples, as Battaglin et al. (2001) suggested. The majority of the NO3 entering streams and rivers in the Corn Belt passes through the soil zone and into the ground water before discharging to surface water. Therefore, it seems likely that most of the denitrification, indicated by the isotopic data, probably occurs within soils and ground water rather than in surface water. This conclusion is supported by Beaumont (2003), who obtained a similar range of isotopic values for shallow ground water samples collected in an agricultural watershed in central Illinois; the majority clustered along the same vector as the tile and Mississippi River samples reported here.

We approximated the degree of denitrification that has probably occurred in the Upper Mississippi River basin and the tile samples collected for this study by applying typical isotopic enrichment factors ({varepsilon}) to the Rayleigh equation. Although Battaglin et al. (2001) and Chang et al. (2002) were not able to distinguish denitrification in their MARB samples, our seasonally collected, more localized water samples followed a well-defined denitrification pattern. Further, we were primarily interested in the overall dentrification exhibited by water in the Mississippi River in comparison with that of tile drain water. To calculate the degree of denitrification for these data, an initial isotopic composition is needed. We chose to use two initial {delta}15N and {delta}18O values. The first estimate uses the average isotopic composition of four samples from the tile drains collected in the spring, during April and May ({delta}15N = 6.1{per thousand} and {delta}18O = 7.2{per thousand}), as the starting isotopic composition. About one-half of the N fertilizer applied to cropland in the midwestern USA is applied in the spring (typically April–May). Precipitation is usually greatest in the spring in Illinois, resulting in more rapid runoff and probably a better representation of initial isotopic composition of recently formed NO3 in the soil zone; rainfall amounts in April and May were slightly below average in 2000 and above average in 2002 and 2004. However, this estimated initial isotopic composition may be conservatively high because these two tiles maintained discharge all year, indicating the influence of local ground water with a longer flow path than most tile drains. Thus, there could be some denitrification already recorded in these spring samples. Lower {delta}15N and {delta}18O values for NO3 samples from tile drains, lysimeters, and shallow wells have been observed in other agricultural settings (Beaumont, 2003; Mengis et al., 2001; Aravena et al., 1993). For example, the average {delta}15N and {delta}18O values for spring samples of tile drains that discharged seasonally in a nearby agricultural field in central Illinois were 4.7{per thousand} and 5.7{per thousand}, respectively (Beaumont, 2003).

The second set of initial {delta}15N and {delta}18O values for calculating the degree of denitrification was determined using theoretical contributions from the different major sources of NO3 as determined by mass-balance calculations for large river systems with emphasis on the Upper Mississippi River basin. Using published studies, we estimated that approximately 14% of the N load is from sewage and manure discharge, 5% from atmospheric deposition, 30% from synthetic fertilizers, and 51% from soil OM (including legumes and pasture) (Howarth et al., 1996; Burkart and James, 1999; Goolsby et al., 1999; David and Gentry, 2000). Using these percentages and the average isotopic compositions of the different major sources of NO3 (Table 3), the theoretical {delta}15N would be ~4.4 {per thousand} and the {delta}18O ~6.4 {per thousand} for the initial NO3 before denitrification effects for the Mississippi River data.


View this table:
[in this window]
[in a new window]
 
Table 3. Typical isotopic values for various sources of NO3 used in the source term model.

 
The isotopic enrichment factor for 15N during the microbial denitrification process can vary significantly, with {varepsilon} values ranging from –5 to nearly –40{per thousand} (Kendall, 1998). The {varepsilon} values observed in soils or laboratory experiments are generally larger than in ground water (Mariotti et al., 1988), except the very large enrichment factors observed in a desert environment (Vogel et al., 1981). Laboratory experiments showed that smaller {varepsilon} values are typically associated with very rapid denitrification rates while the larger values are associated with slow denitrification rates (Mariotti et al., 1988). Enrichment factors observed for studies in temperate environments range from about –4 to –16{per thousand} (Mariotti et al., 1988, Böttcher et al., 1990; Bates et al., 1998; Bates and Spalding, 1998). The environmental conditions, NO3 concentration, availability and concentration of electron donors and bacterial consortium, and reaction rate can result in a large range of {varepsilon} values for denitrification (Kendall, 1998). For this particular study we decided to use {varepsilon} values (–15.9{per thousand} for 15N and –8{per thousand} for 18O) observed by Böttcher et al. (1990) in an unconfined shallow ground water aquifer in an agricultural setting of a temperate climate zone.

Using these published {varepsilon} values and the "initial" isotopic compositions of NO3 estimated from the spring tile drain samples and the theoretical contributions from major NO3 sources, as discussed above, the amount of denitrification for the NO3 in the Mississippi River and the tile drain samples was calculated (Fig. 4). Because the drainage basin of the Upper Mississippi River is large and contains many different soil types and land uses, we feel that applying average or typical isotopic values for the various end-member NO3 sources and using typical {varepsilon} values is a reasonable approach for modeling the data. Assuming that the spread of the isotopic values of the NO3 in water samples collected from the Mississippi River is due primarily to denitrification, these calculated vectors indicate that 0 to 55% of the original NO3 has been denitrified, with most of the samples falling between 10 and 40%. The isotopic values of samples from the tile drains suggest that ~0 to 45% denitrification may have occurred in the water from the two cultivated fields sampled in central Illinois (Fig. 4). The sample from Tile 1 collected in August 2002 (with the largest {delta}15N and {delta}18O values) had the lowest flow and probably represented base-flow conditions (i.e., shallow ground water), suggesting a longer residence time for the NO3 in this sample compared with the other tile samples. We should keep in mind that if smaller {varepsilon} values were used for these calculations, greater degrees of denitrification would be estimated, while larger {varepsilon} values would result in smaller degrees of denitrification. For example, if {varepsilon} values of –10{per thousand} for 15N and –5{per thousand} for 18O were used, then the sample with the largest {delta}15N and {delta}18O values for the Mississippi River samples would be estimated to have gone through approximately 65% denitrification; if {varepsilon} values of –20{per thousand} were used, then the largest degree of denitrification would be approximately 40%. Even with all the uncertainties involved with these calculations, the results indicate that denitrification is an important mechanism for NO3 loss, and help to validate what others have postulated using mass-balance calculations (Howarth et al., 1996; Burkart and James, 1999; Goolsby et al., 1999; David and Gentry, 2000).


Figure 4
View larger version (23K):
[in this window]
[in a new window]
 
Fig. 4. Detailed portion from Fig. 2, showing the two denitrification lines calculated from the two estimated initial {delta}15N and {delta}18O values (gray hexagons). Open symbols represent river samples, closed symbols represent tile drain samples. Letters indicate river sample locations (D = Davenport, Q = Quincy, C = Chester). Tick marks on the vectors indicate an estimation of the extent of denitrification (percentage of original NO3 lost, using {varepsilon} values of –15.9 for {delta}15N and –8 for {delta}18O in Eq. [3] for NO3 having initial values of 4.4{per thousand} for {delta}15N and 6.4{per thousand} for {delta}18O [solid lines] and for NO3 having initial values of 6.1{per thousand} for {delta}15N and 7.2{per thousand} for {delta}18O [dashed lines]).

 
Seasonal Variations
The lower isotopic values of NO3 observed during the winter and spring, when NO3–N concentrations were greatest, correspond to fertilization activities. About one-half of the fertilization in the study area is in early to mid-spring in the midwestern USA (G.F. McIsaac, Univ. of Illinois, personal communication, 2004), a season when precipitation is commonly high. This would lead to relatively rapid movement of NO3 from the soil zone into tile drains, shallow ground water, and river systems. The other one-half of the N fertilizer is applied in Illinois after harvest in late fall. However, soil moisture levels are generally high in the winter and biological activity, including plant growth and denitrification, is inhibited during the cold weather that follows fall fertilizer applications. These conditions lead to less denitrification and increased leaching of NO3 from the soil zones during winter precipitation and thaws.

The isotopic data from the Mississippi River also suggest there was greater denitrification of the NO3 in upstream samples during the summer and fall seasons. The greater degrees of denitrification in the upstream samples may indicate relatively longer travel times in the more northern parts of the Upper Mississippi River basin than the southern parts during the summer and fall seasons. A greater percentage of the soils in Illinois and Iowa are drained with tiles compared with the soils in Wisconsin and Minnesota (Zucker and Brown, 1998). Tile drains provide a shorter and faster route for water to pass through the subsurface, limiting denitrification.


    CONCLUSIONS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
The greatest concentrations of NO3–N for our seasonally collected water samples in both the Mississippi River and tile drain samples were observed in the samples collected during the winter and spring. Although the number of our samples is limited, the timing of the greatest NO3 concentrations agrees with previous studies of the Mississippi and Illinois rivers and follows or coincides with the timing of the major fertilizer application periods for the midwestern USA. The isotopic values of the NO3 ions from the Mississippi River and tile drain samples were lowest during the winter and spring collection periods, when the NO3–N concentrations were greatest.

The isotopic composition of the NO3 from the Mississippi River and the tile drains overlap one another, suggesting that the dominant source of the NO3 for both data sets is similar. The {delta}15N and {delta}18O values for NO3 from the Mississippi River and tile drains follow a linear trend that suggests the isotopic compositions were significantly affected by denitrification. Mixing of different sources of NO3 did not sufficiently explain the range of isotopic data observed for the Mississippi River and tile drain samples. The trend and seasonal variation of the isotopic composition in both data sets indicated that the NO3 is primarily derived from agricultural fertilizers and soil organic N. Although the Mississippi River NO3 should have a greater variety of sources, such as manure and sewage inputs, the isotopic data are not consistent with a large fraction of the NO3 coming from either source, but small contributions from these two sources are likely. Thus, our isotope data support the findings of other investigators who used mass-balance approaches and concluded that elevated NO3–N concentrations in the Mississippi River are primarily due to crop-related agricultural N sources (Turner and Rabalais, 1994; Howarth et al., 1996; Burkart and James, 1999; Goolsby et al., 1999).

Using published {varepsilon} values and calculated starting {delta}15N and {delta}18O compositions for NO3, the estimated degree of denitrification varied from about 0 to nearly 55% for the Mississippi River samples and from about 0 to 45% for the tile drain samples, depending on sample location and season. The majority of the denitrification appears to have occurred before discharge into the Mississippi River, presumably in the soil zone and shallow ground water environments. To establish better estimates of the degree of denitrification that actually take place, additional research is needed to determine the best isotope enrichment factors and help establish the initial isotopic composition of the NO3 for the various sources in the agriculturally dominated midwestern USA.


    ACKNOWLEDGMENTS
 
This research was supported through a grant from the Illinois Groundwater Consortium at Southern Illinois University in Carbondale and by the Illinois Department of Natural Resources. We thank the Iowa-American Water Co., the Quincy Water Co., and the Chester Water Co. for their assistance in collecting river-water samples. We also thank E. Mehnert, B. Herzog, and J. Goodwin of the ISGS, and H. Wehrmann, E. Krug, T. Holm, and D. Winstanley of the ISWS for their critical reviews and insightful comments. Publication of this article has been authorized by the Chiefs of the Illinois State Geological and Water Surveys.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 





This Article
Right arrow Abstract Freely available <