Old, abandoned feedlots may serve as a source of nutrients that can degrade groundwater and downstream water quality. We characterized the distribution and concentration of nutrients at the Crookston Cattle Company feedlot (northwest Minnesota, USA), 15 years after it ended operations in 1999. Groundwater nitrate concentration decreased from 55 mg/L (as nitrogen) in 2003 to less than 5 mg/L since 2007. Results from stable isotope analysis, with δ15N and δ18O in groundwater nitrate ranging up to +44 and +30‰, respectively, suggest denitrification as the cause, rather than either nitrate transport from the site or dilution. Phosphorus, with soil B-horizon concentrations as much as 112 and averaging 24 mg/kg, is sequestered by carbonate-rich glacial sediments and, serendipitously, an iron-rich sand deposit formed millennia ago by wave action along the shore of glacial Lake Agassiz. Map analysis indicates roughly 20,000 kg of P in excess of background concentration remains in soil at the 15 ha site. Evidence suggests that the former feedlot has not affected water quality significantly in an agricultural ditch that drains the feedlot and its vicinity. Rather than originating from the feedlot, small increases of total phosphorus observed in the downstream ditch likely result from release of phosphorus from nearby recently restored wetlands. More consequential than elevated nutrient concentrations to the future reclamation of this and similar sites is the persistence of robust non-native species. Our results suggest that before development, feedlot sites should be evaluated for their phosphorus sequestration and denitrification potential, thus mitigating the potential for later off-site transport of nutrients.

Excessive agricultural nutrients – nitrogen (N) and phosphorus (P) – in soils and waterways are a global problem (Carpenter et al. 1998; Smith 2003), with abundant quantities generated in feedlots. In the U.S. alone, cattle feedlots generate an estimated 18 million metric tons of nutrient-rich manure per year (Sweeten 2000), which affects runoff water quality. On average, an adult cow excretes 55 kg N and 18 kg P per year (Goolsby et al. 1999), although these nutrient amounts can vary significantly (e.g., Davis et al. 2002). Thus, for even small feedlots, the nutrients generated can easily total in the tens of thousands of kilograms. While much N is lost to the atmosphere as ammonia (Hristov et al. 2011) and some of the N and P is spread and broadcast beyond the confines of the feedlot (Everett and Vickery 2005), a large portion may be transported vertically and then laterally after leaching into soils that underlie the facility (García et al. 2012; Vaillant et al. 2009; Chang and Entz 1996). In this report, we suggest that feedlots may create “hot spots” of N and P that can serve as a potentially long-term source of nutrients, which can lead to downstream eutrophication, harm nearby wetlands, and degrade natural plant communities.

The fate of nutrients in feedlots is controlled by biogeochemical processes that depend significantly on hydrogeological factors. Nitrate from agricultural sources, including manure and feedlots, has resulted in the widespread degradation of water resources (i.e., Power and Schepers 1989), with limnic and marine zones effected by runoff from agricultural watersheds, such as the Gulf of Mexico near the outlet of the Mississippi River (i.e., Rabalais et al. 2002). In response to this major environmental problem, much research to-date has focused on the transport and transformation of nitrate that occurs by leaching and denitrification (i.e., Korom 1992; Schipper et al. 2010). Dual analyses of oxygen and nitrogen isotopes in nitrate (reviewed by Kendall 1998), in particular, reveal the source and fate of nitrate in a variety of environments, including those that host livestock (Fenech et al. 2012). In many cases, denitrification processes can be tracked by the dual isotope approach (Böttcher et al. 1990; Fukada et al. 2003).

In contrast to nitrate, phosphorus is typically the limiting nutrient in freshwater eutrophication (Schindler 1981; Vollenweider and Kerekes 1982). P is not lost to the atmosphere and is instead sequestered by precipitation of mineral phases and adsorption onto Fe-hydroxides, clays, carbonates, and organic matter. In some circumstances, P can be mobilized and lost to nearby streams and wetlands if the soils hold excessive P concentrations, contain abundant organic matter that may cause chemical reduction, or are deep and permeable (Eghball et al. 1996; Sims et al. 1998; Eghball 2003). Recent research indicates subsurface P transport is not negligible (King et al. 2015). Furthermore, large volumes of solid manure may lead to saturation of P retention capacity, which will enhance the capacity for P transport (Ciapparelli et al. 2016).

The fate of P in soil is determined by a variety of competing processes. Depending on the pH and oxidation reduction conditions, P can be precipitated as Ca-, Fe-, or Al-phases, or can be adsorbed onto organic or clay particles (Djodjic 2001; Sekhon 2002). Phosphorus adsorption, however, competes with fixation of organic acids generated from the manure (Chardon and Schoumans 2002) and the loss of adsorption capacity by the formation of Fe and Al complexes (Bolster and Sistani 2009). Once soil adsorption capacity is reached, which may occur in soils influenced strongly by organic matter (Ciapparelli et al. 2016), P can be mobilized (Campbell and Racz 1975; Eghball 2003; Glaesner et al. 2012) and transported to the water table where low oxygen and concomitant reduction may increase P solubility. Thus, coupled with the subsurface movement of organic matter, the geological characteristics of soils and near-surface sediments – texture, composition, and the oxidation-reduction conditions in the subsurface – likely control the transport and concentration of P in and around feedlots.

Coupled with excessive input of nutrients, site operations adjacent to pens and large numbers of stock lead to disturbance of native vegetation communities and establishment of invasive exotic species. Tognetti and Chaneton (2015) showed how mowing in pastures (and often done adjacent to feedlots) combined with N addition led to rapid invasion by short-lived exotic forbs and later by exotic perennial grasses. Also sensitive to disturbance are wetlands and riparian zones that may lie near or adjacent to disturbed areas, including feedlots. For example, alien and invasive species such as Phragmites australis and Lythrum salicaria are clear indicators of wetland disturbance along Connecticut (USA) roadways (Moore et al. 1999). In Illinois (USA) restored wetlands, non-native invasion and dominance were found to be controlled in larger measure by abiotic factors, such as nitrogen availability and disturbance, rather than propagule strength and resiliency, which enhanced non-native richness (Matthews et al. 2009). In other wetlands and riparian areas, disturbance combined with increasing nutrients favor invasion of monospecific, fast-growing, robust, and morphologically plastic species such as Phalaris (Kercher et al. 2007) and Phragmites (Minchinton and Bertness 2003).

Lastly, feedlot status-of-use (Vaillant et al. 2009) and age of the site (Zhu et al. 2004) are important considerations related to possible environmental degradation and contamination, with older and abandoned sites posing a greater risk because of diminished management and typically fewer environmental protection regulations that were active during operation. Mitigating “legacy” P derived from agricultural sources presents challenges because the traditional conservation approach, reducing soil erosion and transport to surface water, may not provide a sufficiently robust management strategy (Kleinman et al. 2011). Under optimum hydrological conditions, even low concentrations of legacy P in cropped fields, which are likely to have lower concentrations of P than feedlots, can contribute significant loads of P to runoff. Tillage management, field amendments, and engineered runoff filters can be used to mitigate P release and transport, but these alternatives often demand costly, focused efforts that are rarely used in agricultural operations (Kleinman et al. 2011).

In our report, we show to what degree a single, abandoned feedlot influenced and continues to affect its surrounding environment. For this case study, we investigate the spatial and temporal distribution of soil and groundwater N and P in and near an abandoned feedlot in north-central USA, which operated between 1970 and 1999. First, we describe the landscape and current distribution of nutrient concentrations at the former feedlot site. Next, in the context of soil composition and potential biogeochemical processes, we suggest reasons for the observed variability of nutrients across the site and evaluate the degree of transport into adjacent wetlands and surface-water channels. Finally, our report describes more generally the likely biogeochemical processes that determine the mobility of feedlot nutrients, along with the ecological consequences.

Site description of the abandoned cattle feedlot. Our case example, the former Crookston Cattle Company feedlot (47°43.6’ N, 96°19.2’ W), lies within the US Fish and Wildlife Service’s Glacial Ridge National Wildlife Refuge near Crookston, northwestern Minnesota, USA. (Figure 1, insert). The feedlot contained as many as 2,500 cattle at one time during its 30-year operation, although the number of livestock was generally much less, with an average of 500 head for the period of operation, including the times when the livestock was on the range. Similar numbers of sheep occupied the lot during the last few years of the feedlot’s operation (1995–2000). The pens on the northeast end of the feedlot were used more than those toward the southwest (Figure 1) (Terrill Bradford, University of Minnesota-Crookston, personal communication).

Figure 1

Map of the feedlot site. Insert shows the location of Glacial Ridge National Wildlife Refuge in northwest Minnesota, USA. Dots and triangles indicate the locations of soil and groundwater samples (Table S1), respectively, with narrow gray boxes showing former feedlot pens. The open circle symbol locates US Geological Survey well G-8. Overlay symbol represent current wetland areas. X- and Y-scales are in UTM zone 14 coordinates (NAD83, in meters). DOI: https://doi.org/10.1525/elementa.271.f1

Figure 1

Map of the feedlot site. Insert shows the location of Glacial Ridge National Wildlife Refuge in northwest Minnesota, USA. Dots and triangles indicate the locations of soil and groundwater samples (Table S1), respectively, with narrow gray boxes showing former feedlot pens. The open circle symbol locates US Geological Survey well G-8. Overlay symbol represent current wetland areas. X- and Y-scales are in UTM zone 14 coordinates (NAD83, in meters). DOI: https://doi.org/10.1525/elementa.271.f1

The feedlot was established on a gentle, low-relief sandy ridge. Formed roughly 9,900 years ago, the ridge developed along a prominent strandline on the eastern margin of glacial Lake Agassiz (Wright, 1972). This location was selected for development of the livestock facility because of enhanced air circulation on the ridge and ample natural subsurface drainage. The former feedlot pens and wastewater ponds cover about 15 ha and are drained by a ditch system, which conveys surface runoff into the Red Lake River, a tributary to the Red River of the North and downstream Lake Winnipeg. In recent decades, Lake Winnipeg has become increasingly eutrophic, due largely to agricultural activity and the regional expansion and increasing efficiency of drainage within the Red River portion of its watershed (e.g., McCullough et al. 2012; Wassenaar and Rao 2012).

The Natural Resources Conservation Service (NRCS) Digital Soil Survey Geographic Database (SSURGO) digital maps show that the soil series underlying the site consist of excessively well-drained Sandberg – Radium complex loamy sand along the highest axis of the ridge, Radium loamy sand or Syrene sandy loam along the flanks, and Strathcona fine sandy loam in adjacent wetlands. These soils are Mollisols with a combined O-A-B horizon that ranges from about 0.4 m in wetlands to 0.8 m in more elevated sites (NRCS 2015).

Historical imagery (Figure 2) shows that partial drainage of wetlands in the vicinity of the feedlot occurred before 1939, with legal records suggesting excavation of the main ditch to the east and south in the 1910s. Because of the sandy, excessively well-drained soils that underlie the feedlot, this area was never cultivated but used for grazing and forage before the feedlot was established in the 1970s. In contrast, effective drainage of the former wetland basins toward the east and south allowed for cultivation, which began in the 1990s. As part of the Glacial Ridge Project (Gerla et al. 2012), many of these ditches were filled in 2005 and formerly drained areas, where accessible, were planted in native perennial grasses. Neither the National Wetlands Inventory (US Fish and Wildlife Service 2015) nor SSURGO hydric soils (NRCS 2015) clearly delineate wetlands near the feedlot. In addition, only components of some of the soil series to the northeast and south are hydric; not the entire area shown actually functions as wetland. Thus, wetland-upland boundaries shown near the feedlot are generalized (Figure 2).

Figure 2

Historical imagery. Aerial images of the feedlot taken in 1939, 1954, 2003, and 2015 and surrounding wetlands. Most of the area shown within the image was originally converted from pasture and forage to cultivated crops between 1970 and 1995, and then restored to its original, more natural land cover after 2005. Hydric soil distribution (SSURGO digital data, NRCS 2015) is shown with the bright overlay. Note that areas to the northeast and southeast of the feedlot are shown as underlain by hydric soils, but these soils occur as only a component in the dominantly upland soils in these areas. The former feedlot pens, which were constructed in about 1970, are delineated by bold lines in the 2003 and 2015 images. DOI: https://doi.org/10.1525/elementa.271.f2

Figure 2

Historical imagery. Aerial images of the feedlot taken in 1939, 1954, 2003, and 2015 and surrounding wetlands. Most of the area shown within the image was originally converted from pasture and forage to cultivated crops between 1970 and 1995, and then restored to its original, more natural land cover after 2005. Hydric soil distribution (SSURGO digital data, NRCS 2015) is shown with the bright overlay. Note that areas to the northeast and southeast of the feedlot are shown as underlain by hydric soils, but these soils occur as only a component in the dominantly upland soils in these areas. The former feedlot pens, which were constructed in about 1970, are delineated by bold lines in the 2003 and 2015 images. DOI: https://doi.org/10.1525/elementa.271.f2

The wetland basin north of the feedlot was never entirely drained (Figure 2) but received feedlot runoff and waste that flowed through an unlined sewage lagoon excavated in the 1970s. Prior to that, the wetland was likely used for pasture and grazing (Figure 2, upper images), followed by no management and non-use while the feedlot was in operation, and then finally outlets were filled at least partially during the Glacial Ridge Project. This attempt at limited restoration led to a recent condition where open water occurs seasonally (Figure 2, lower right). Because of the wetter than long-term average conditions that have affected the region since 1993 (Winter and Rosenberry 1998), the presence of open water in the basin unlikely existed except recently, at least during historical times.

Groundwater and soil investigation. Groundwater monitoring began in 2002, with extensive sampling and analysis of groundwater and soil carried out in 2013 using methods described in detail by Gbolo (2014). After laying out an irregular grid with transects oriented perpendicular to the trend of the beach ridge (Figure 1 and Table S1), soil samples were collected at 63 sites, prepared, and analyzed for texture, pH, electrical conductivity, total nitrogen, nitrate, ammonium, available phosphorus (Olsen et al. 1954), total carbon, organic carbon, calcium, and iron. Samples were collected from both the combined O- and A-horizon and soil B-horizon (Table S2 and S3), so depths varied depending on the landscape position and soil series. Plant associations were mapped by combining field locations of boundaries determined with GPS together with visual identification of vegetation communities on recent Farm Service Agency aerial imagery.

Groundwater was sampled from 16 shallow water table wells at the site (Figure 1 and Table S4), which ranged in depth from 2.4–3.2 meters and had either 0.3 or 0.6 m sand-packed, mill-slot PVC, 5 cm (2 inch) diameter screens. The wells were thoroughly purged before sampling, and samples were transported on ice to the University of North Dakota Environmental Analytical Research Laboratory (UND EARL), with analyses completed within the time limit recommended by the US Environmental Protection Agency (1983). Groundwater samples were obtained from these wells with screens that were placed either across or immediately below the mean level of the water table, always within 0 to 1 m of the tension saturated zone. An additional 21 more single and nested wells, along with four open pits fitted with a standpipe were used to monitor water levels estimate lateral and vertical hydraulic gradients at the site (Gbolo 2014). Lastly, monitoring data from one monitoring well installed by the U.S. Geological Survey (G08-R 149N44W17ABAD 0000620668) (Figure 1, open circle) supplemented groundwater geochemistry presented in this report.

Surface water investigation. Surface water was sampled within wetlands adjacent to the feedlot on the northwest and southeast, and along Polk County Judicial Ditch 66, the nearest surface-water channel, which lies roughly 500 m east of the feedlot (Figure 3). Much of the water that is conveyed by the ditch is groundwater discharged from a large sand and gravel production pit that lies 800 m east of the feedlot, maintaining at least 0.006 cubic meters per second during the driest times and as much as 0.2 cubic meters per second during spring runoff. Three sites were sampled at various times between 2007 and 2013: near the outlet of the aggregate production pit, downstream about 800 m (nearest the feedlot), and roughly 1,700 m farther downstream. Along with the feedlot and adjacent wetlands, this ditch also receives runoff from surrounding former cropland, which is now reconstructed prairie and wetland. The surface-water analytes and laboratory techniques used were the same as those for groundwater analysis.

Figure 3

Vegetation communities in 2013 and groundwater flow direction. Arrows indicate the general direction of groundwater flow based on water-table elevations monitored during the course of the study. Dotted lines show nearby ditches (Judicial Ditch 66 lies on the far east of the map and flows north), with the former feedlot pens outlined and highlighted. DOI: https://doi.org/10.1525/elementa.271.f3

Figure 3

Vegetation communities in 2013 and groundwater flow direction. Arrows indicate the general direction of groundwater flow based on water-table elevations monitored during the course of the study. Dotted lines show nearby ditches (Judicial Ditch 66 lies on the far east of the map and flows north), with the former feedlot pens outlined and highlighted. DOI: https://doi.org/10.1525/elementa.271.f3

Chemical analysis. Chemical analyses of the soil extracts and water samples were performed at the UND EARL. A Shimadzu TOC analyzer (Model VCSH) was used to analyze for total carbon, inorganic carbon, and total organic carbon. The filtered (0.45 μm) unpreserved samples were used for nitrate-nitrogen, nitrite-nitrogen, and soluble or dissolved reactive phosphorus analyses. Nitrate-nitrogen and nitrite-nitrogen were measured using ion chromatography (Dionex DX-120). The filtered (0.45 μm) and acid-preserved samples were used for ammonium-nitrogen analysis. Soluble reactive phosphorus, total phosphorus, and ammonium-nitrogen were measured using a HACH DR/2010 spectrophotometer. Soluble reactive phosphorus (SRP), total phosphorus (TP), and ammonium-nitrogen (NH4-N) analyses used ammonia nesslerization, acid-hydrolysable, and ascorbic acid methods, respectively. Calcium and iron were analyzed using the EARL’s inductively coupled argon plasma-atomic emission spectrometer.

Nitrogen and oxygen stable isotopes in groundwater nitrate were collected between 2009 and 2013 and analyzed by mass spectrophotometry at the University of Nebraska – Lincoln’s Water Sciences Center using the laboratory’s VG Optima dual inlet isotope ratio mass spectrometer. Water samples for isotope analysis were collected from wells using a disposable bailer, filtered (0.45 µm), frozen at < –18°C, and shipped by overnight delivery to the laboratory in a well-insulated package.

Statistical analysis and mapping. A square grid (x = y = 11.6 m) was created from the scattered soil and groundwater analytical data using a multi-quadric radial basis function algorithm (Carlson and Foley 1991a), with an R-squared smoothing factor (Carlson and Foley 1991b) selected to highlight the trends in analyte concentration across the feedlot and its vicinity. To create a grid, the radial basis function approach uses a summation of quadric surface regression equations with coefficients determined by values of known points across the mapped area. This method generates a grid that filters local irregularities while maintaining fully the trend of the variable (e.g., Hardy 1971). The output from the gridding process was also used to estimate the volume of soil and groundwater that remains affected by feedlot-related nutrients. This was done by applying the extended Simpson’s 3/8 Rule, whereby parabolas are used to approximate z-values for each part of the curving surface. Volumes were then estimated by integrating z between the underlying clay till and the ground surface (for soil) or water table (for groundwater).

Groundwater flow. Monitoring from 2010–2014 showed that the water-table level beneath the feedlot varied within a range of 0.5 to 1 m, the largest difference occurring in the coarsest and thickest sediments, and diminishing towards the wetlands. Clay till underlies the site at greater depth, forming an essentially impervious lower boundary at a depth ranging from 0.5 m along the margins to 9 m along the northeast-southwest trending axis of the sandy ridge (Gbolo and Gerla 2013). These observations revealed that a broad, low-relief groundwater mound or crest, maintained by shallow groundwater recharge, occurs within the glacial Lake Agassiz beach sediments that underlie the feedlot.

Groundwater flow, based on water-level monitoring (Gbolo 2014), follows a simple pattern of recharge along the northeast to southwest-trending sandy ridge, accompanied by lateral movement and discharge toward wetlands on the northwest and southeast (Figure 3). This relatively uncomplicated and constrained flow system, underlain by clay till, provided conditions conducive to assess nutrient transport.

Ecohydrological setting. Plant community maps (Figure 3) show that the feedlot is flanked on the northwest and southeast by partially drained wetlands characterized by generally monospecific stands of reed canary grass (Phalaris arundinacea), hybrid cattail (Typha sp.), and common reed (Phragmites australis), with more upland sites covered by grassland vegetation dominated by a largely non-native bluegrass (Poa sp.), quackgrass (Elymus repens), field thistle (Cirsium discolor), and smooth brome (Bromus inermis), with some original (not restored) big blue stem (Andropogon gerardii) in less disturbed areas, including the areas that underlie the former feedlot animal pens.

These generally invasive plant communities tend to correspond closely to the distribution of soil series, with Typha, Phalaris, and Phragmites occurring in areas underlain by Strathcona fine-sandy loam and the well-drained Radium and Syrene series hosting Poa, Elymus, Cirsium, and Bromus. Regionally, less disturbed well-drained soils have somewhat lower concentrations of organic matter and nutrients when compared to soils on beach-ridge flanks and wetland margins, but the pattern is reversed at the former feedlot. For example, mean nitrogen and phosphorus content of the sandy soils on the ridge is at least double or as much as an order of magnitude more, respectively, even with only small differences in the organic matter content (Tables 1 and 2).

Table 1

O- and A-horizon soil chemical composition related to plant communitiesa. DOI: https://doi.org/10.1525/elementa.271.t1

NO3PSOMpHEC

 
All Communities (n = 63) (mg/kg) (mg/kg) (%)  (µS/cm) 
minimum, maximum 0.5, 26.5 2, 121 0.1, 44.7 7.5, 10.0 22, 2260 
mean and standard deviation 5.2 ± 6.0 38 ± 43 12.8 ± 9.4 8.5 ± 0.6 400 ± 469 
Typha (n = 8)      
minimum, maximum 0.5, 8.0 2, 8 5.2, 44.7 7.8, 9.8 153, 2240 
mean and standard deviation 3.8 ± 3.1 5 ± 2 18.8 ± 11.4 8.5 ± 0.7 790 ± 693 
Phalaris (n = 15)      
minimum, maximum 2.0, 26.5 3, 54 12.0, 33.1 7.6, 9.0 162, 2260 
mean and standard deviation 12.0 ± 8.0 10 ± 13 23.3 ± 7.9 8.3 ± 0.4 613 ± 523 
Bromus and Andropogon (n = 21)      
minimum, maximum 0.5, 10.5 3, 104 0.8, 21.0 7.5, 9.4 22, 1567 
mean and standard deviation 2.8 ± 2.8 25 ± 31 8.1 ± 5.1 8.5 ± 0.6 242 ± 338 
Poa, Cirsium, Elymus (n = 19)      
minimum, maximum 0.5, 12.0 3, 121 0.1, 13.4 7.5, 10.0 30, 1142 
mean and standard deviation 3.0 ± 3.2 91 ± 30 7.3 ± 3.5 8.5 ± 0.7 243 ± 265 
NO3PSOMpHEC

 
All Communities (n = 63) (mg/kg) (mg/kg) (%)  (µS/cm) 
minimum, maximum 0.5, 26.5 2, 121 0.1, 44.7 7.5, 10.0 22, 2260 
mean and standard deviation 5.2 ± 6.0 38 ± 43 12.8 ± 9.4 8.5 ± 0.6 400 ± 469 
Typha (n = 8)      
minimum, maximum 0.5, 8.0 2, 8 5.2, 44.7 7.8, 9.8 153, 2240 
mean and standard deviation 3.8 ± 3.1 5 ± 2 18.8 ± 11.4 8.5 ± 0.7 790 ± 693 
Phalaris (n = 15)      
minimum, maximum 2.0, 26.5 3, 54 12.0, 33.1 7.6, 9.0 162, 2260 
mean and standard deviation 12.0 ± 8.0 10 ± 13 23.3 ± 7.9 8.3 ± 0.4 613 ± 523 
Bromus and Andropogon (n = 21)      
minimum, maximum 0.5, 10.5 3, 104 0.8, 21.0 7.5, 9.4 22, 1567 
mean and standard deviation 2.8 ± 2.8 25 ± 31 8.1 ± 5.1 8.5 ± 0.6 242 ± 338 
Poa, Cirsium, Elymus (n = 19)      
minimum, maximum 0.5, 12.0 3, 121 0.1, 13.4 7.5, 10.0 30, 1142 
mean and standard deviation 3.0 ± 3.2 91 ± 30 7.3 ± 3.5 8.5 ± 0.7 243 ± 265 

a NO3 – N, nitrate as nitrogen; SOM, soil organic matter; EC, electrical conductivity; n, number of samples.

Table 2

B-horizon soil chemical composition related to plant communitiesa. DOI: https://doi.org/10.1525/elementa.271.t2

NO3PFeCaSOMpHEC

 
All Communities (n = 60) (mg/kg) (mg/kg) (mg/kg) (mg/kg) (%)  (µS/cm) 
minimum, maximum 0.5, 20.0 1, 112 4.2, 97.5 2720, 11540 0.2, 10.2 7.7, 9.6 30, 1308 
mean, sample standard deviation 2.5 ± 3.9 24 ± 35 26.4 ± 25.5 6843 ± 2619 3.7 ± 2.4 8.6 ± 0.5 327 ± 294 
Typha (n = 8)        
minimum, maximum 0.5, 3.0 1,4 6.9, 61.5 6060, 10720 0.2, 5.7 8.2, 9.1 103, 1308 
mean, sample standard deviation 1.4 ± 0.9 2 ± 1 18.2 ± 18.0 8770 ± 1506 3.6 ± 1.7 8.6 ± 0.3 448 ± 384 
Phalaris (n = 15)        
minimum, maximum 0.5, 11.0 1, 12 5.2, 97.5 4320, 11540 1.8, 6.6 7.8, 9.4 34, 1284 
mean, sample standard deviation 2.6 ± 2.6 4 ± 3 20.7 ± 24.6 8839 ± 2283 3.3 ± 1.6 8.5 ± 0.5 418 ± 397 
Bromus and Andropogon (n = 21)        
minimum, maximum 0.5, 3.0 2, 83 4.2, 92.5 2720, 10940 1.0, 10.2 7.7, 9.6 30, 664 
mean, sample standard deviation 1.2 ± 0.7 13 ± 18 19.0 ± 20.7 5904 ± 2495 3.7 ± 2.6 8.7 ± 0.6 222 ± 200 
Poa, Cirsium, Elymus (n = 16)        
minimum, maximum 0.5, 20.0 2, 112 5.4, 93.0 3020, 8580 0.3, 10.0 7.7, 9.4 61, 617 
mean, sample standard deviation 4.6 ± 6.6 70 ± 37 45.5 ± 27.3 5240 ± 1680 4.0 ± 3.2 8.5 ± 0.5 317 ± 195 
NO3PFeCaSOMpHEC

 
All Communities (n = 60) (mg/kg) (mg/kg) (mg/kg) (mg/kg) (%)  (µS/cm) 
minimum, maximum 0.5, 20.0 1, 112 4.2, 97.5 2720, 11540 0.2, 10.2 7.7, 9.6 30, 1308 
mean, sample standard deviation 2.5 ± 3.9 24 ± 35 26.4 ± 25.5 6843 ± 2619 3.7 ± 2.4 8.6 ± 0.5 327 ± 294 
Typha (n = 8)        
minimum, maximum 0.5, 3.0 1,4 6.9, 61.5 6060, 10720 0.2, 5.7 8.2, 9.1 103, 1308 
mean, sample standard deviation 1.4 ± 0.9 2 ± 1 18.2 ± 18.0 8770 ± 1506 3.6 ± 1.7 8.6 ± 0.3 448 ± 384 
Phalaris (n = 15)        
minimum, maximum 0.5, 11.0 1, 12 5.2, 97.5 4320, 11540 1.8, 6.6 7.8, 9.4 34, 1284 
mean, sample standard deviation 2.6 ± 2.6 4 ± 3 20.7 ± 24.6 8839 ± 2283 3.3 ± 1.6 8.5 ± 0.5 418 ± 397 
Bromus and Andropogon (n = 21)        
minimum, maximum 0.5, 3.0 2, 83 4.2, 92.5 2720, 10940 1.0, 10.2 7.7, 9.6 30, 664 
mean, sample standard deviation 1.2 ± 0.7 13 ± 18 19.0 ± 20.7 5904 ± 2495 3.7 ± 2.6 8.7 ± 0.6 222 ± 200 
Poa, Cirsium, Elymus (n = 16)        
minimum, maximum 0.5, 20.0 2, 112 5.4, 93.0 3020, 8580 0.3, 10.0 7.7, 9.4 61, 617 
mean, sample standard deviation 4.6 ± 6.6 70 ± 37 45.5 ± 27.3 5240 ± 1680 4.0 ± 3.2 8.5 ± 0.5 317 ± 195 

a NO3 – N, nitrate as nitrogen; P, total phosphorus; Fe, total iron; Ca, calcium; SOM, soil organic matter; EC, electrical conductivity; n, number of samples.

Soil and groundwater nitrogen. Soil nitrate in the shallow O- and A-horizon soils (Table 1 and Table S2) has its greatest concentrations along the margins of the feedlot pens (in excess of 10 mg/kg), with the largest concentrations associated with the waste lagoon southeast of the site (Figure 4). In contrast, soil nitrogen in the B-horizon (Table 2 and Table S3) is more strongly concentrated along the southeastern margin of the feedlot. Nitrogen distribution in B-horizon soils is similarly distributed spatially to shallow groundwater nitrate concentrations and situated almost directly below the feedlot pens, suggesting that B-horizon soil and nitrate have a similar geochemical control on concentrations. Nitrite and ammonium concentrations in water samples were all near or below the detection limit (Table S4).

Figure 4

Nitrogen in soil and groundwater. Concentration of total nitrogen (mg/kg) in soil O-A and B-horizon at the feedlot, indicated in the solid colors. Contours lines show the mean concentration of nitrate (mg/L, as nitrogen) measured during the monitoring period. Dots show the locations of soil samples, with narrow gray boxes indicating former feedlot pens. Small triangles are shallow groundwater monitoring wells. The X- and Y-scales are in UTM zone 14 coordinates (NAD83, in meters). DOI: https://doi.org/10.1525/elementa.271.f4

Figure 4

Nitrogen in soil and groundwater. Concentration of total nitrogen (mg/kg) in soil O-A and B-horizon at the feedlot, indicated in the solid colors. Contours lines show the mean concentration of nitrate (mg/L, as nitrogen) measured during the monitoring period. Dots show the locations of soil samples, with narrow gray boxes indicating former feedlot pens. Small triangles are shallow groundwater monitoring wells. The X- and Y-scales are in UTM zone 14 coordinates (NAD83, in meters). DOI: https://doi.org/10.1525/elementa.271.f4

Nitrogen isotopes. A large range in oxygen and nitrogen isotopes in nitrate was recorded for the feedlot and surrounding crop land, ranging from 0 to +44‰ for δ15N and 0 to +53‰ for δ18O (Figure 5 and Table S5). In comparison to nitrate, only minor concentrations of ammonium-N and nitrite-N were detected in soils and groundwater at the feedlot site, with mean concentrations for both roughly one order of magnitude less than nitrate-N (Table S4).

Figure 5

Nitrogen and oxygen isotope ratios from nitrate in groundwater. Nitrogen and oxygen isotope ratios (‰) in shallow groundwater nitrate measured for the feedlot and crop land soils within a two-km radius of the feedlot. Boxes show the typical fields for various sources of nitrate (adapted from Kendall 1998 and Mengis et al. 2001). The insert shows the trend of isotopic differentiation during various biogeochemical processes. DOI: https://doi.org/10.1525/elementa.271.f5

Figure 5

Nitrogen and oxygen isotope ratios from nitrate in groundwater. Nitrogen and oxygen isotope ratios (‰) in shallow groundwater nitrate measured for the feedlot and crop land soils within a two-km radius of the feedlot. Boxes show the typical fields for various sources of nitrate (adapted from Kendall 1998 and Mengis et al. 2001). The insert shows the trend of isotopic differentiation during various biogeochemical processes. DOI: https://doi.org/10.1525/elementa.271.f5

Soil and groundwater phosphorus. As expected in the coarse textured, aerobic soils along the beach ridge, soil and groundwater phosphorus concentrations are largest directly below the pens at the feedlot: 112 mg/kg and <0.1 – 0.35 mg/L, respectively. Results suggest minimal transport of P from the feedlot (Figure 6) when compared to the background concentrations, which are soil P at roughly 3 mg/kg and total dissolved P < 0.1 mg/L, inferred from values near the fringe of the feedlot and ungrazed/undisturbed areas elsewhere within the Glacial Ridge National Wildlife Refuge (Gbolo 2014). Note that the two “hotspots” are in roughly the same position as elevated nitrate in the soil B-horizon (Figure 4). The pattern and composition of P in the O-A soil horizon (Table S2) is similar to the B-horizon (Table S3) and therefore not shown or described in detail.

Figure 6

Phosphorus in soil and groundwater. Concentration of total phosphorus (mg/kg) in soil B-horizon at the feedlot is indicated in the solid colors. Contours lines show the mean concentration of total phosphorus (mg/L) recorded during the monitoring period. Dots show the locations of soil samples, with narrow gray boxes indicating former feedlot pens. Small triangles are shallow groundwater monitoring wells. The X- and Y-scales are in UTM zone 14 coordinates (NAD83, in meters). DOI: https://doi.org/10.1525/elementa.271.f6

Figure 6

Phosphorus in soil and groundwater. Concentration of total phosphorus (mg/kg) in soil B-horizon at the feedlot is indicated in the solid colors. Contours lines show the mean concentration of total phosphorus (mg/L) recorded during the monitoring period. Dots show the locations of soil samples, with narrow gray boxes indicating former feedlot pens. Small triangles are shallow groundwater monitoring wells. The X- and Y-scales are in UTM zone 14 coordinates (NAD83, in meters). DOI: https://doi.org/10.1525/elementa.271.f6

Calcium, iron, and organic carbon concentrations in soil. As noted in the introduction, the transport and sequestration of P and nitrate are likely controlled by calcium, iron, and organic carbon. These were analyzed in the soils sampled at the feedlot (Tables S2 and S3). Results showed that calcium ranges from about 3,500 to 10,000 mg/kg, with a spatial distribution that is strongly related to the sandy beach ridge that underlies the former feedlot (Figure 7). Greatest concentrations lie along the lower flanks of the ridge, and appear largest in the proximity of the feedlot lagoons.

Figure 7

Calcium and total iron in soil. Concentration of total calcium (mg/kg) (top) and total extractable iron (mg/kg) (bottom) in soil B-horizon at the feedlot, indicated in the solid colors separated by contours. Dots show the locations of soil samples, with narrow gray boxes indicating former feedlot pens. The X- and Y-scales are in UTM zone 14 coordinates (NAD83, in meters). DOI: https://doi.org/10.1525/elementa.271.f7

Figure 7

Calcium and total iron in soil. Concentration of total calcium (mg/kg) (top) and total extractable iron (mg/kg) (bottom) in soil B-horizon at the feedlot, indicated in the solid colors separated by contours. Dots show the locations of soil samples, with narrow gray boxes indicating former feedlot pens. The X- and Y-scales are in UTM zone 14 coordinates (NAD83, in meters). DOI: https://doi.org/10.1525/elementa.271.f7

In contrast, the largest concentrations of Fe (about 65 mg/kg) generally coincide with the southeastern flank of the ridge, and fall to less than 20 mg/kg in the wetlands that lie both to the northwest and southeast (Figure 7). Organic carbon content in the A horizon revealed a pattern that is similar to calcium (Figure 8), with the lowest soil concentration (<5%) occurring beneath the feedlot (within the beach ridge) and increasing to as much as 30% toward wetlands, and averaging 15%.

Figure 8

Organic carbon in soil. Percent concentration of soil organic carbon (OC) soil O and A-horizon at the feedlot, indicated in the solid colors separated by contours. Dots show the locations of soil samples, with narrow gray boxes indicating former feedlot pens. The X- and Y-scales are in UTM zone 14 coordinates (NAD83, in meters). DOI: https://doi.org/10.1525/elementa.271.f8

Figure 8

Organic carbon in soil. Percent concentration of soil organic carbon (OC) soil O and A-horizon at the feedlot, indicated in the solid colors separated by contours. Dots show the locations of soil samples, with narrow gray boxes indicating former feedlot pens. The X- and Y-scales are in UTM zone 14 coordinates (NAD83, in meters). DOI: https://doi.org/10.1525/elementa.271.f8

Nitrate and phosphorus in surface water. Nitrate and phosphorus concentration in surface water drained from the vicinity of the feedlot varied widely during the sampling and analysis from July 2007 through June 2013. This was also the time during which the surface water drain was retrofitted to a more natural pattern (Gerla et al. 2012), with excavation and re-grading contemporaneous with wetland restoration downstream in June 2010.

Nitrate concentration in surface waters remained well below the U.S. EPA drinking water standard of <10 mg/L nitrate-N for the upstream and downstream stations for all sample events (Figure 9). Generally, the nitrate varied between not detected (<0.1 mg/L) and 1 mg/L, with no obvious difference between the two stations, although a concentration of 2.2 mg/L was reported at the upstream – aggregate pit outlet sample site on one occasion (Figure 9).

Figure 9

Downstream water quality. Measured concentrations of nitrate (as N, mg/L) (left) and total phosphorus (mg/L) (right) at the aggregate pit (top), which lies 1.5 km upstream from the feedlot, and at the point where Judicial Ditch 66 flows under US Highway 2, approximately 4 km downstream from the feedlot. Annotations show when the pit was inactive in mid 2008 (“no pit excavation”) and when a major ditch retrofit was constructed. DOI: https://doi.org/10.1525/elementa.271.f9

Figure 9

Downstream water quality. Measured concentrations of nitrate (as N, mg/L) (left) and total phosphorus (mg/L) (right) at the aggregate pit (top), which lies 1.5 km upstream from the feedlot, and at the point where Judicial Ditch 66 flows under US Highway 2, approximately 4 km downstream from the feedlot. Annotations show when the pit was inactive in mid 2008 (“no pit excavation”) and when a major ditch retrofit was constructed. DOI: https://doi.org/10.1525/elementa.271.f9

Total phosphorus showed generally very low concentrations throughout the monitoring period 2007–2011 (Figure 9). Concentrations had greater variability at the upstream station in when compared to the lower sample location. All but four of the analyses revealed concentrations at or below 0.1 mg/L, with the other four ranging up to 0.23 mg/L at the upper sampling point in July 2007 (Figure 9). A surface-water concentration of 0.1 mg/L total P (100 ppb) has been adopted by many states as the water quality threshold.

An increasing trend in total P concentration occurred at the downstream sample site but not at the upstream monitoring station, with this latter site showing a wide range of concentration. Although not investigated further, this sharply varying concentration may correspond to periods of aggregate extraction or the presence and absence of large flocks of migratory birds in flooded aggregate pits. In contrast, the apparent rising trend in phosphorus concentration at the downstream site from 2009–2011 may have resulted from soil and channel disturbance during the ditch retrofit and then later with the restoration of a large wetland basin immediately upstream from the sample station.

Land use changes on plant communities. Aerial imagery from 1939, 1954, 2003, and 2013 shows significant changes in land use and land cover between 1954 and 2003. Prior to at least 1954, the upland areas at the feedlot were used for forage (Figure 3). Changes of ownership of the tract in the 1970s and 1980s began the conversion of pasture and forage to crops (Gerla et al. 2012). After abandonment of the feedlot, especially in the 1990s, areas to the southeast and northeast adjacent to the feedlot were drained, tilled, and cultivated with wheat, soybean, corn rotation. During the Glacial Ridge project, ditches were partially filled and the formerly drained areas were seeded with native perennial grass species, starting in 2005. Previous wetland areas that were not drained and the area encompassing the feedlot pens, however, were not reclaimed. Currently, most if not all of these areas are characterized by plant communities with dominantly non-native, generally invasive species (Figure 2). In particular, the ubiquitous hybrid Typha and Phalaris arundinacea stands in and around wetlands indicate that restoration to native species with moderate diversity will be difficult and likely impossible (Zedler and Kercher 2005), even with a return to former hydrological patterns. Diminishing nitrate and overall sequestration of phosphorus will likely do little to mitigate the disturbance related to the feedlot, with present conditions providing a niche for invasive species (e.g., James et al. 2010; Levine et al. 2003).

Oxygen and nitrogen isotopes in groundwater nitrate. Stable isotopes of nitrogen and oxygen in groundwater nitrate suggest a variety of nitrate sources at the site (Figure 5), especially manure together with nitrate- and ammonium-based fertilizers. Although it is unlikely that nitrogen fertilizers were ever directly applied to soil at the feedlot, cultivation and fertilizing equipment was stored at the northeastern end of the site. In addition, handling and storage of forage used for cattle may have leached naturally and released nutrients into soils. Fractionation revealed by the unusually large range of δ15N from about 0‰ to more than +40‰, suggests strongly that past and on-going active denitrification occurs at the feedlot (e.g., Kendall 1998; Kendall et al. 2007; Komor and Anderson 1993). Although data are not shown temporally, the largest δ15N ratios (>+30‰) were recorded for oldest samples collected (2010), and have isotopic signatures that are significantly different from nearby crop land. Three conditions must be met for denitrification to occur: 1) the presence of denitrifying bacteria, 2) sufficiently reduced, anaerobic conditions, and 3) availability of electron donors (Firestone 1982; Korom 1992). These conditions are likely met at the feedlot, with the flux of organic carbon through the groundwater flow system serving to create sufficiently reducing conditions and providing an electron donor for denitrification reactions.

The US Geological Survey has been monitoring groundwater conditions since 2003 (Cowdery et al. 2008) at a water-table well located immediately northeast of the feedlot. Data show a marked change in nitrate concentration (Figure 10), with a major drop occurring in 2005–2006, perhaps triggered by unusually dry conditions in 2006. Nitrate near the water table would have been transported to deeper levels where anaerobic conditions prevail during normal conditions; this translocation of nitrate followed by recharge and a return to normal levels of the water table may have accelerated denitrification and loss of nitrate from the groundwater system. Alternatively, recent work indicates that large changes in nitrate loss to streams can be triggered by a strong change in prevailing precipitation, or “weather whiplash” (Loecke et al. 2017). Dry springs following wet autumns in the upper Mississippi River basin have been shown to decrease nitrate concentrations in streams and rivers. Conversely, a wet spring following a dry autumn results in elevated spring nitrate concentrations (Loecke et al. 2017). Although unusually dry conditions during 2006 correlate with the sudden drop in groundwater nitrate at the feedlot, concentrations did not recover even partially when normal weather conditions and water-table levels returned (Figure 10), nor is there evidence for transport of nitrate to nearby drains (Figure 9). When considered together with stable isotope trends (Figure 5), denitrification seems to provide a better explanation for the loss of nitrate.

Figure 10

Temporal change in shallow groundwater nitrogen concentration. Concentration of nitrate + nitrate (mg/L as N) from US Geological Survey well G-8, which has been monitored since 2003. (Data are from the National Water Information Service, accessed June 2015, http://waterdata.usgs.gov/nwis, site number: 474346096185501, site name: G08-R 149N44W17ABAD 0000620668). DOI: https://doi.org/10.1525/elementa.271.f10

Figure 10

Temporal change in shallow groundwater nitrogen concentration. Concentration of nitrate + nitrate (mg/L as N) from US Geological Survey well G-8, which has been monitored since 2003. (Data are from the National Water Information Service, accessed June 2015, http://waterdata.usgs.gov/nwis, site number: 474346096185501, site name: G08-R 149N44W17ABAD 0000620668). DOI: https://doi.org/10.1525/elementa.271.f10

Phosphorus sequestration. Phosphorus tends to be sequestered by oxidizing conditions coupled with (1) formation of ferric iron and calcium-phosphate minerals, (2) adsorption onto clay, organic matter, and iron hydroxide minerals, as previously noted, and (3) assimilated and recycled by vegetation (Kaila 1949; Cole et al. 1977; Damon et al. 2014). Although the greatest potential for sequestration lies at the margins of the beach ridge and feedlot, much P remains beneath the feedlot (Figure 6). This suggests effective isolation of P by calcium, iron, and possibly wetland organic matter. These processes are likely to occur within the soil B-horizon, characterized by an accumulation of soluble or suspended organic material, clay, iron, and aluminum transported from shallower soil horizons (Table 2). The water table tends to rise and fall through this soil horizon at the feedlot, resulting in local sorption and desorption of soil P.

Elevated concentrations of calcium occur in wetlands directly adjacent to the sandy ridge (Figure 7, top). The beach ridge sediments are reworked from the coarse fraction of the Red Lake Falls Formation, which directly underlies the feedlot and contains a large proportion of carbonate rock fragments (Harris et al. 1974). One explanation for the calcium-concentration distribution may be that groundwater recharge becomes increasingly calcareous as infiltration occurs in the carbonate-rich sediments of the ridge. As this carbonate-saturated/supersaturated groundwater is discharged in adjacent wetland margins, carbonate precipitates as carbon dioxide is lost from the water (Stumm and Morgan 1996) thereby increasing the concentration of calcium in these wetland-fringe soils.

Total extractable iron shows a peculiar pattern of elevated concentration in the soil B-horizon (Figure 7, bottom) that is enriched along the southeastern margin of the feedlot and continues at elevated concentrations up to 60 mg/kg toward the southwest. One explanation may be the occurrence of heavy sand placers along the back-beach shore of glacial Lake Agassiz. Beach wave action can concentrate relatively inert minerals with high specific gravity, commonly magnetite and hematite. These heavy sand placers are well known and have been mapped along a wide range of shore line environments (e.g., Stanaway 2012; Komar 2007), although in the glacial Lake Agassiz region placers are likely small and of little value. For example, at several locations along the east coast of Australia, heavy sands accumulated in a “back-barrier wash-over facies, at the rear of swash-aligned barriers in coastal embayments” (Roy 1999), which describes features similar to the glacial Lake Agassiz-age beach sediments at the feedlot. These small iron-rich zones were likely deposited as fine magnetite and hematite sand retained in the back-barrier area during storms on glacial Lake Agassiz, later weathering to form an iron-rich zone in the soil B-horizon.

Large concentrations of both iron and calcium, which developed naturally along the margins of the beach ridge (Figure 7), suggest that P from the feedlot was and continues to be effectively isolated from down-gradient receptors. This occurs because of potentially strong sequestration by Fe- and Ca-phosphate mineral precipitation, with transport further mitigated by organic matter adsorption in wetlands surrounding the feedlot.

Phosphorus remaining at the feedlot. The total mass of phosphorus remaining at the feedlot can be estimated roughly from the dimensions of the groundwater and soil P anomaly. For groundwater, the average porosity of the saturated sediments above the underlying clay (0.15) was multiplied by the sediment volume (2.7 × 105 m3) and the average concentration of P that exceeds the background value (0.1 mg/L). Results showed that only about 4 kg of P occur dissolved in sediment pores. A similar estimate can be carried out for the P sequestered within the soils. In this case, the volume of the sand and silt that lies above the clay till (3.6 × 105 m3) is multiplied by both the dry bulk density (1,400 kg m3) and the average concentration of P above background (35 mg/kg), giving a result of 1.8 × 104 kg of P. Therefore, only a very small fraction of the total P at the site occurs in groundwater (0.02%).

Unfortunately, no documentation exists on how many cattle were sheltered at the feedlot during its 30 years of existence, other than an average of about 500 head during the years of the operation, nor is there any record of the amount of manure removed from the site. Regardless, much of the manure would have been moved off-site and spread on nearby crop fields. Assuming that 500 cattle were present in the lot, on average, along with a yearly P deposition rate of 18 kg P per cow (Goolsby et al. 1999), then perhaps 2.7 × 105 kg of P would have been generated, which is roughly 15 times more than presently sequestered at the site.

An additional consideration on the fate of phosphorus relates to the oxidation-reduction dynamics once all of the nitrate has been denitrified. Carbon continues to be transported to groundwater from decaying and leaching of vegetation at the surface, which will lead to increasingly reduced conditions in the shallow groundwater. The sequence of microbial-mediated oxidation-reduction processes proceeds from NO3 / N2 to MnO2(mineral) / Mn2+ to Fe(OH)3 (mineral) / Fe2+ oxidation-reduction couples (e.g., Stumm and Morgan 1996). The oxidant with the highest thermodynamic potential is used initially, followed by the one with the next lower potential. Because only small concentrations of nitrite and manganese occur at the site, microbial reduction and dissolution of iron minerals which sequester P may occur once nitrate is completely used. As a result, this process may trigger a later mobilization of phosphorus.

Nutrient transport and downstream water quality. In addition to the focused study at the feedlot, the water quality within nearby Judicial Ditch 66 (Figure 2) was tracked during and prior to the geochemical characterization and monitoring at the abandoned feedlot. As observed elsewhere, a rise in phosphorus transport in runoff often occurs following restoration of wetlands (Aldous et al. 2005, 2007). In the case of Judicial Ditch 66, an increase of total P had been observed in ditch discharge flowing northward out of the Glacial Ridge National Wildlife Refuge following wetland restoration in the basin, although nitrate had decreased during the same period (Figure 9). Because of the strong sequestering of P at the feedlot, the observed transport of P from the ditch system was instead a likely result of the chemical reduction of soils and sediments in the recently restored wetlands. Newly saturated soils will lose much of their capacity for ferric adsorption, thereby releasing P to surface waters. This process is also consistent with the observed reduction in nitrate concentration (Figure 9), as nitrate is readily denitrified in reducing conditions (e.g., Korom 1992) brought on by soil saturation. Monitoring of Judicial Ditch 66 both up- and downstream from the feedlot suggested that aggregate mining and the disturbance from prairie restoration may account for the small downstream increase in P, rather than transport from the abandoned feedlot. Furthermore, restoration and a concomitant return to anaerobic conditions has likely helped mitigate excessive nitrate throughout the system.

Old, abandoned feedlots may serve as a source of nutrients that can degrade groundwater and downstream water quality. At the Crookston Cattle Company feedlot, which operated from 1970 through 1999 and abandoned in 2000, some nutrients remain and a former waste lagoon continues to serve as a source of nitrate. N isotopes and temporal monitoring, however, suggest that most nitrate has been denitrified. Phosphorus remains strongly fixed beneath the feedlot and the site is surrounded by soils that have even greater sequestration capacity. Diminishing organic matter and dissolved organic carbon occurred after cattle were removed from the feedlot, thereby mitigating continued transport of organic P. This indicates that P will likely be retained to a greater extent at the site, provided the land use/land cover do not change significantly on the site’s calcareous soils.

In the case of the Crookston Cattle Company’s landscape setting, results indicate that the partitioning of nutrients was similar between the soil B-horizon and groundwater, suggesting that site characterization may not have required costly installation of groundwater monitoring equipment. Instead, soil sampling and routine analysis would suffice to characterize the transport and storage of P. This suggestion, however, might only apply to sites with a similar shallow water table and deep B-horizon soils.

Results indicate that although N and P concentrations are not excessive in adjacent wetlands and do not appear to have been transported into downstream channels. Nonetheless, robust invasive vegetation is well established in the nutrient-rich soils surrounding the feedlot and would be difficult to reclaim to more natural conditions. Disturbance, even decades earlier, provided an avenue for establishment of non-native weeds, which was later and continues to be augmented by abundant nutrient resources. Ditch closure and return to a more temporally and spatially consistent, natural hydrological pattern following the surrounding site restoration did little to mitigate weeds, and therefore does not alone assure or provide a trajectory toward recovery.

The ability to assess thoroughly the outcome of the nutrients released to the surrounding environment is limited greatly by our lack of feedlot records, including the number of animal units held in the pens at various times, the nutrient concentration of excrement, and the approximate proportion of manure retained and removed from the site. Nonetheless, learning from the legacy nutrients at Crookston Cattle Company feedlot suggests that on-site investigation might be carried out to characterize the suitability of potential future feedlot sites. Locations could be selected that provide the best conditions for denitrification, which requires the presence of anoxic soils and labile organic carbon or other electron donors. Soils that are well drained and contain large concentrations of organic matter, calcium, and ferric iron, similar to those at the Crookston Cattle Company feedlot, would likely be those most favorable for mitigating nitrate transport and enhancing phosphorus sequestration.

All data used are presented in this published report or its accompanying supplemental files.

The authors wish to thank the U.S. Fish & Wildlife Service Plains & Prairie Potholes Landscape Conservation Cooperative (PPP-LCC, administered through the Red Lake Watershed District) and the North Dakota Water Resources Research Institute for financial support. We also appreciate The Nature Conservancy’s permission to access and monitor the former feedlot site. Several anonymous reviewers, along with Ariane Peralta, East Carolina University, especially, provided thorough reviews and suggestions for improving the manuscript.

Most funding for this project came from the U.S. Fish & Wildlife Service’s Plains & Prairie Potholes Landscape Conservation Cooperative agreement 301817J119M4, administered through the Red Lake Watershed District. Additional funding came from the North Dakota Water Resources Research Institute and the Harold Hamm School of Geology and Geological Engineering.

The authors have no competing interests to declare.

  • Contributed to conception and design: PJG, PG

  • Contributed to acquisition of data: PJG, PG

  • Contributed to analysis and interpretation of data: PJG, PG

  • Drafted and/or revised the article: PJG, PG

  • Approved the submitted version for publication: PJG

1
Aldous
,
AR
,
Craft
,
CB
,
Stevens
,
CJ
,
Barry
,
MJ
and
Bach
,
LB
2007
Soil phosphorus release from a restoration wetland, upper Klamath Lake, Oregon
.
Wetlands
27
(
4
):
1025
1035
. DOI:
2
Aldous
,
AR
,
McCormick
,
P
,
Ferguson
,
C
,
Graham
,
S
and
Craft
,
CB
2005
Hydrologic regime controls soil phosphorus fluxes in restoration and undisturbed wetlands
.
Restoration Ecology
13
(
2
):
341
347
. DOI:
3
Bolster
,
CH
and
Sistani
,
KR
2009
Sorption of phosphorus from swine, dairy, and poultry manures
.
Communications in Soil Science and Plant Analysis
40
:
1106
1123
. DOI:
4
Böttcher
,
J
,
Strebel
,
O
,
Voerkelius
,
S
and
Schmidt
,
H-L
1990
Using isotope fractionation of nitrate-nitrogen and nitrate-oxygen for evaluation of microbial denitrification in a sandy aquifer
.
Journal of Hydrology
114
(
3–4
):
413
424
. DOI:
5
Campbell
,
LB
and
Racz
,
GJ
1975
Organic and inorganic P content, movement, and mineralization of P in soil beneath a feedlot
.
Canadian Journal of Soil Science
55
:
457
466
. DOI:
6
Carlson
,
RE
and
Foley
,
TA
1991a
Radial Basis Interpolation Methods on Track Data, Lawrence Livermore National Laboratory
, UCRL-JC-1074238.
7
Carlson
,
RE
and
Foley
,
TA
1991b
The parameter R2 in multiquadric interpolation
.
Computers Math. Applic
.
21
(
9
):
29
42
. DOI:
8
Carpenter
,
SR
,
Caraco
,
NF
,
Correll
,
DL
,
Howarth
,
RW
,
Sharpley
,
AN
and
Smith
,
VH
1998
Nonpoint pollution of surface waters with phosphorus and nitrogen
.
Ecological Applications
8
(
3
):
559
568
. DOI:
9
Chang
,
C
and
Entz
,
T
1996
Nitrate leaching losses under repeated cattle feedlot manure applications in southern Alberta
.
Journal of Environmental Quality
25
(
1
):
145
153
. DOI:
10
Chardon
,
WJ
and
Schoumans
,
OF
2002
Solubilization of phosphorus: concepts and process description of chemical mechanisms. In:
Chardon
,
WJ
and
Schoumans
,
OF
(eds.),
Phosphorus losses from agricultural soils: processes at the field scale
,
42
52
.
Wageningen, The Netherlands
:
COST Action 832 Alterra
.
11
Ciapparelli
,
IC
,
de Iorio
,
AF
and
García
,
AR
2016
Phosphorus downward movement in soil highly charged with cattle manure
.
Environmental Earth Sciences
75
(
7
), art. no. 568. DOI:
12
Cole
,
CV
,
Innis
,
GS
and
Stewart
,
JWB
1977
Simulation of phosphorus cycling in semiarid grasslands
.
Ecology
58
:
2
15
. DOI:
13
Cowdery
,
TK
,
Lorenz
,
DL
and
Arntson
,
AD
2008
Hydrology Prior to Wetland and Prairie Restoration in and around the Glacial Ridge National Wildlife Refuge, Northwestern Minnesota, 2002–2005
vi
:
68
.
Washington DC
:
US Geological Survey Scientific Investigations Report
2007
5200
.
14
Damon
,
PM
,
Bowden
,
B
,
Rose
,
T
and
Rengel
,
Z
2014
Crop residue contributions to phosphorus pools in agricultural soils: A review
.
Soil Biology and Biochemistry
74
:
127
137
. DOI:
15
Davis
,
JG
,
Iversen
,
KV
and
Vigil
,
MF
2002
Nutrient variability in manures: Implications for sampling and regional database creation
.
Journal of Soil and Water Conservation
57
(
6
):
473
478
.
16
Djodjic
,
F
2001
Displacement of phosphorus in structured soils [dissertation].
Uppsala, Sweden
:
Swedish University of Agricultural Sciences
.
Acta Universitatis Agriculturae Sueciae, Agraria
283
:
30
. ISSN: 1401-6249. ISBN: 91-576-5826-9.
17
Eghball
,
B
2003
Leaching of phosphorus fractions following manure or compost application
.
Communications in Soil Science and Plant Analysis
34
(
19–20
):
2803
2815
. DOI:
18
Eghball
,
B
,
Binford
,
GD
and
Baltensperger
,
DD
1996
Phosphorus movement and adsorption in a soil receiving long-term manure and fertilizer application
.
J. Environmental Quality
25
:
1339
1343
. DOI:
19
Everett
,
LA
and
Vickery
,
J
2005
Land application of manure: Minnesota livestock producers’ practices and educational needs
.
Journal of Extension
43
(
3
):
11
, art. no. 3RIB8.
20
Fenech
,
C
,
Rock
,
L
,
Nolan
,
K
,
Tobin
,
J
and
Morrissey
,
A
2012
The potential for a suite of isotope and chemical markers to differentiate sources of nitrate contamination: A review
.
Water Research
46
(
7
):
2023
2041
. DOI:
21
Firestone
,
MK
1982
Biological denitrification. In:
Stevenson
,
FJ
(ed.),
Nitrogen in Agricultural Soils
.
Madison Wisconsin, USA
:
American Society of Agronomy, Agronomy Monograph
22
:
289
326
. DOI:
22
Fukada
,
T
,
Hiscock
,
KM
,
Dennis
,
PF
and
Grischek
,
T
2003
A dual isotope approach to identify denitrification in groundwater at a river-bank infiltration site
.
Water Research
37
(
13
):
3070
3078
. DOI:
23
García
,
AR
,
Maisonnave
,
R
,
Massobrio
,
MJ
and
Fabrizio de Iorio
,
AR
2012
Field-scale evaluation of water fluxes and manure solution leaching in feedlot pen soils
.
Journal of Environmental Quality
41
(
5
):
1591
1599
. DOI:
24
Gbolo
,
P
2014
Quantifying nutrient cycling and fate within an abandoned feedlot and adjacent wetlands [dissertation].
Grand Forks, North Dakota, USA
:
University of North Dakota, Harold Hamm School of Geology and Geological Engineering
.
25
Gbolo
,
P
and
Gerla
,
PJ
2013
Statistical analysis to characterize the transport of nutrients in groundwater near an abandoned feedlot
.
Hydrol. Earth Syst. Sci
.
17
:
4897
4906
. DOI:
26
Gerla
,
PJ
,
Cornett
,
MW
,
Ekstein
,
JD
and
Ahlering
,
MA
2012
Talking big: Lessons learned from a 9000 hectare restoration in the northern tallgrass prairie
.
Sustainability
4
(
11
):
3066
3087
. DOI:
27
Glaesner
,
N
,
Donner
,
E
,
Magid
,
J
,
Rubaek
,
GH
,
Zhang
,
H
and
Lombi
,
E
2012
Characterization of leached phosphorus from soil, manure, and manure-amended soil by physical and chemical fractionation and diffusive gradients in thin films (DGT)
.
Environmental Science & Technology
46
:
10564
10571
. DOI:
28
Goolsby
,
DA
,
Battaglin
,
WA
,
Lawrence
,
GB
,
Artz
,
RS
,
Aulenbach
,
BT
,
Hooper
,
RP
,
Keeney
,
DR
and
Stensland
,
GJ
1999
Flux and Sources of Nutrients in the Mississippi-Atchafalaya River Basin: Topic 3 Report for the Integrated Assessment on Hypoxia in the Gulf of Mexico.
NOAA Coastal Ocean Program Decision Analysis Series No. 17
.
Silver Spring
,
MD NOAA Coastal Ocean Program
.
29
Hardy
,
RL
1971
Multiquadric equations of topography and other irregular surfaces
.
J. Geophysical Research
76
(
8
):
1905
1915
. DOI:
30
Harris
,
KL
,
Moran
,
SR
and
Clayton
,
L
1974
Late Quaternary Stratigraphic Nomenclature, Red River Valley, North Dakota and Minnesota.
North Dakota Geological Survey Miscellaneous Series 52
.
Bismarck, North Dakota
,
47
.
31
Hristov
,
AN
,
Hanigan
,
M
,
Cole
,
A
,
Todd
,
R
,
McAllister
,
TA
,
Ndegwa
,
PM
and
Rotz
,
A
2011
Review: Ammonia emissions from dairy farms and beef feedlots
.
Canadian Journal of Animal Science
91
(
1
):
1
35
. DOI:
32
James
,
JJ
,
Ziegenhagen
,
L
and
Aanderud
,
ZT
2010
Exploitation of nutrient-rich soil patches by invasive annual and native perennial grasses
.
Invasive Plant Science and Management
3
(
2
):
169
177
. DOI:
33
Kaila
,
A
1949
Biological absorption of phosphorus
.
Soil Science
68
(
4
):
279
290
. DOI:
34
Kendall
,
C
1998
Chapter 16 – Tracing nitrogen sources and cycling in catchment. In:
Kendall
,
C
and
McDonnell
,
JJ
(eds.),
Isotope Tracers in Catchment Hydrology
,
519
576
.
Amsterdam, Netherlands
:
Elsevier
, ISBN-13:978-0444501554. DOI:
35
Kendall
,
C
,
Elliott
,
EM
and
Wankel
,
SD
2007
Chapter 12 – Tracing anthropogenic inputs of nitrogen to ecosystems. In:
Michener
,
RH
and
Lajtha
,
K
(eds.),
Stable Isotopes in Ecology and Environmental Science
(2nd edition),
375
449
.
Oxford, England
:
Blackwell Publishing
, ISBN-13: 978-1405126809. DOI:
36
Kercher
,
SM
,
Herr-Turoff
,
A
and
Zedler
,
JB
2007
Understanding invasion as a process: The case of Phalaris arundinacea in wet prairies
.
Biological Invasions
9
(
6
):
657
665
. DOI:
37
King
,
KW
,
Williams
,
MR
,
Macrae
,
ML
,
Fausey
,
NR
,
Frankenberger
,
J
,
Smith
,
DR
,
Kleinman
,
PJA
and
Brown
,
LC
2015
Phosphorus transport in agricultural subsurface drainage: A review
.
Journal of Environmental Quality
44
(
2
):
467
485
. DOI:
38
Kleinman
,
PJA
,
Sharpley
,
AN
,
Buda
,
AR
,
McDowell
,
RW
and
Allen
,
AL
2011
Soil controls of phosphorus in runoff: Management barriers and opportunities
.
Canadian Journal of Soil Science
91
(
3
):
329
338
. DOI:
39
Komar
,
PD
2007
Chapter 1 – The Entrainment, Transport and Sorting of Heavy Minerals by Waves and Currents. In:
Mange
,
MA
and
Wright
,
DT
(eds.),
Heavy Minerals in Use – Developments in Sedimentology
58
.
Amsterdam
:
Elsevier
. ASIN: B01A0C5DJU. DOI:
40
Komor
,
SC
and
Anderson
,
HW
Jr.
,
1993
Nitrogen isotopes as indicators of nitrate sources in Minnesota sand-plain aquifers
.
Ground Water
31
(
2
):
260
270
. DOI:
41
Korom
,
SF
1992
Natural denitrification in the saturated zone: a review
.
Water Resources Research
28
(
6
):
1657
1668
. DOI:
42
Levine
,
JM
,
Vilà
,
M
,
D’Antonio
,
CM
,
Dukes
,
JS
,
Grigulis
,
K
and
Lavorel
,
S
2003
Mechanisms underlying the impacts of exotic plant invasions
.
Proceedings of the Royal Society B: Biological Sciences
270
(
1517
):
775
781
. DOI:
43
Loecke
,
TD
,
Burgin
,
AJ
,
Riveros-Iregui
,
DA
,
Ward
,
AS
,
Thomas
,
SA
,
Davis
,
CA
and
Clair
,
MAS
2017
Weather whiplash in agricultural regions drives deterioration of water quality
.
Biogeochemistry
133
(
1
):
7
15
. DOI:
44
Matthews
,
JW
,
Peralta
,
AL
,
Soni
,
A
,
Baldwin
,
P
,
Kent
,
AD
and
Endress
,
AG
2009
Local and landscape correlates of non-native species invasion in restored wetlands
.
Ecography
32
(
6
):
1031
1039
. DOI:
45
McCullough
,
GK
,
Page
,
SJ
,
Hesslein
,
RH
,
Stainton
,
MP
,
Kling
,
HJ
,
Salki
,
AG
and
Barber
,
DG
2012
Hydrological forcing of a recent trophic surge in Lake Winnipeg
.
Journal of Great Lakes Research
38
(
SUPPL. 3
):
95
105
. DOI:
46
Mengis
,
M
,
Walther
,
U
,
Bernasconi
,
SM
and
Wehrli
,
B
2001
Limitations of using δ18O for the source identification of nitrate in agricultural soils
.
Environ. Sci. Technol
.
35
(
2001
):
1840
1844
. DOI:
47
Minchinton
,
TE
and
Bertness
,
MD
2003
Disturbance-mediated competition and the spread of Phragmites australis in a coastal marsh
.
Ecological Applications
13
(
5
):
1400
1416
. DOI:
48
Moore
,
HH
,
Niering
,
WA
,
Marsicano
,
LJ
and
Dowdell
,
M
1999
Vegetation change in created emergent wetlands (1988–1996) in Connecticut (USA)
.
Wetlands Ecology and Management
7
(
4
):
177
191
. DOI:
49
Natural Resources Conservation Service (NRCS)
2015
Soil Survey Geographic (SSURGO) Database.
Washington DC
:
US Department of Agriculture, Natural Resources Conservation Service
. Accessed at: http://soildatamart.nrcs.usda.gov/ on 6 June 2015.
50
Olsen
,
SR
,
Cole
,
CV
,
Watanabe
,
FS
and
Dean
,
LA
1954
Estimation of available phosphorus in soils by extraction with sodium bicarbonate.
Washington DC
:
US Department of Agriculture Circular
939
.
51
Power
,
JF
and
Schepers
,
JS
1989
Nitrate contamination of groundwater in North America
.
Agriculture, Ecosystems and Environment
26
(
3–4
):
165
187
. DOI:
52
Rabalais
,
NN
,
Turner
,
RE
and
Wiseman
,
WJ
, Jr.
2002
Gulf of Mexico hypoxia, a.k.a. “The dead zone”
.
Annual Review of Ecology and Systematics
33
:
235
263
. DOI:
53
Roy
,
PS
1999
Heavy mineral beach placers in southeastern Australia: Their nature and genesis
.
Economic Geology
94
(
4
):
567
588
. DOI:
54
Schindler
,
DW
1981
Studies of eutrophication in lakes and their relevance to the estuarine environment. In:
Neilson
,
BJ
and
Cronin
,
LE
(eds.),
Contemporary Issues in Science and Society – Estuaries and Nutrients
71
82
.
Humana Press
, DOI:
55
Schipper
,
LA
,
Robertson
,
WD
,
Gold
,
AJ
,
Jaynes
,
DB
and
Cameron
,
SC
2010
Denitrifying bioreactors-An approach for reducing nitrate loads to receiving waters
.
Ecological Engineering
36
(
11
):
1532
1543
. DOI:
56
Sekhon
,
BS
2002
Modeling of soil phosphorus sorption and control of phosphorus pollution with acid mine drainage floc [dissertation].
Morgantown, West Virginia, USA
:
West Virginia University, Davis College of Agriculture, Forestry, and Consumer Sciences Dissertation
,
224
.
57
Sims
,
JT
,
Simard
,
RR
and
Joern
,
BC
1998
Phosphorus loss in agricultural drainage: Historical perspective and current research
.
Journal of Environmental Quality
27
(
2
):
277
293
. DOI:
58
Smith
,
VH
2003
Eutrophication of freshwater and coastal marine ecosystems: A global problem
.
Environmental Science and Pollution Research
10
(
2
):
126
139
. DOI:
59
Stanaway
,
KJ
2012
Ten placer deposit models from five sedimentary environments
.
Transactions of the Institutions of Mining and Metallurgy, Section B: Applied Earth Science
121
(
1
):
43
51
. DOI:
60
Stumm
,
W
and
Morgan
,
JJ
1996
Aquatic Chemistry
, Third Edition.
New York
:
John Wiley and Sons
. ASIN: B010WIN2RC.
61
Sweeten
,
JM
2000
Manure and wastewater management for cattle feedlots
.
Reviews of Environmental Contamination and Toxicology
167
:
121
153
. DOI:
62
Tognetti
,
PM
and
Chaneton
,
EJ
2015
Community disassembly and invasion of remnant native grasslands under fluctuating resource supply
.
Journal of Applied Ecology
52
(
1
):
119
128
. DOI:
63
US Environmental Protection Agency
1983
Sample Preservation. In:
Methods for Chemical Analysis of Water and Wastes
.
xv
xx
. EPA-600/4-79-020
Cincinnati, Ohio, USA
.
64
US Fish and Wildlife Service
2015
National Wetlands Inventory.
Washington DC
:
Department of the Interior. US Fish and Wildlife Service
. Accessed at: http://www.fws.gov/wetlands/ on 9 June 2015.
65
Vaillant
,
GC
,
Pierzynski
,
GM
,
Ham
,
JM
and
DeRouchey
,
J
2009
Nutrient accumulation below cattle feedlot pens in Kansas
.
Journal of Environmental Quality
38
(
3
):
909
918
. DOI:
66
Vollenweider
,
RA
and
Kerekes
,
J
1982
Eutrophication of Waters. Monitoring, Assessment and Control
.
Organization for Economic Co-Operation and Development (OECD)
.
Paris, France
,
156
.
67
Wassenaar
,
LI
and
Rao
,
YR
2012
Lake Winnipeg: The forgotten great lake
.
Journal of Great Lakes Research
38
(
SUPPL. 3
):
1
5
. DOI:
68
Winter
,
TC
and
Rosenberry
,
DO
1998
Hydrology of prairie pothole wetlands during drought and deluge: A 17-year study of the Cottonwood Lake wetland complex in North Dakota in the perspective of longer term measured and proxy hydrological records
.
Climatic Change
40
(
2
):
189
209
. DOI:
69
Wright
,
HE
, Jr.
1972
Quaternary History of Minnesota. In:
Sims
,
PK
and
Morey
,
GB
(eds.),
Geology of Minnesota
515
548
.
Minneapolis, Minnesota, USA
:
Minnesota Geological Survey Geological Survey
.
70
Zedler
,
JB
and
Kercher
,
S
2005
Wetland resources: Status, trends, ecosystem services, and restorability
.
Annual Review of Environment and Resources
30
:
39
74
. DOI:
71
Zhu
,
J
,
Luo
,
A
,
Zhou
,
Y
,
Ndegwa
,
PM
and
Schmidt
,
DR
2004
The age effect of dairy feedlots on manure nutrient seepage in loam soils
.
Biosystems Engineering
89
(
2
):
223
229
. DOI:
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International License (CC-BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. See http://creativecommons.org/licenses/by/4.0/.