Reducing Nitrogen Export from the Corn Belt to the Gulf of Mexico: Agricultural Strategies for Remediating Hypoxia

SPAtially Referenced Regression on Watershed models developed for the Upper Midwest were used to help evaluate the nitrogen‐load reductions likely to be achieved by a variety of agricultural conservation practices in the Upper Mississippi‐Ohio River Basin (UMORB) and to compare these reductions to the 45% nitrogen‐load reduction proposed to remediate hypoxia in the Gulf of Mexico (GoM). Our results indicate that nitrogen‐management practices (improved fertilizer management and cover crops) fall short of achieving this goal, even if adopted on all cropland in the region. The goal of a 45% decrease in loads to the GoM can only be achieved through the coupling of nitrogen‐management practices with innovative nitrogen‐removal practices such as tile‐drainage treatment wetlands, drainage–ditch enhancements, stream‐channel restoration, and floodplain reconnection. Combining nitrogen‐management practices with nitrogen‐removal practices can dramatically reduce nutrient export from agricultural landscapes while minimizing impacts to agricultural production. With this approach, it may be possible to meet the 45% nutrient reduction goal while converting less than 1% of cropland in the UMORB to nitrogen‐removal practices. Conservationists, policy makers, and agricultural producers seeking a workable strategy to reduce nitrogen export from the Corn Belt will need to consider a combination of nitrogen‐management practices at the field scale and diverse nitrogen‐removal practices at the landscape scale.     

Nutrient Sources and Transport in the Missouri River Basin,with Emphasis on the Effects of Irrigation and Reservoirs1

Brown, Juliane B., Lori A. Sprague, and Jean A. Dupree, 2011. Nutrient Sources and Transport in the Missouri River Basin, With Emphasis on the Effects of Irrigation and Reservoirs. Journal of the American Water Resources Association (JAWRA) 47(5):1034‐1060. DOI: 10.1111/j.1752‐1688.2011.00584.x Abstract: SPAtially Referenced Regressions On Watershed attributes (SPARROW) models were used to relate instream nutrient loads to sources and factors influencing the transport of nutrients in the Missouri River Basin. Agricultural inputs from fertilizer and manure were the largest nutrient sources throughout a large part of the basin, although atmospheric and urban inputs were important sources in some areas. Sediment mobilized from stream channels was a source of phosphorus in medium and larger streams. Irrigation on agricultural land was estimated to decrease the nitrogen load reaching the Mississippi River by as much as 17%, likely as a result of increased anoxia and denitrification in the soil zone. Approximately 16% of the nitrogen load and 33% of the phosphorus load that would have otherwise reached the Mississippi River was retained in reservoirs and lakes throughout the basin. Nearly half of the total attenuation occurred in the eight largest water bodies. Unlike the other major tributary basins, nearly the entire instream nutrient load leaving the outlet of the Platte and Kansas River subbasins reached the Mississippi River. Most of the larger reservoirs and lakes in the Platte River subbasin are upstream of the major sources, whereas in the Kansas River subbasin, most of the source inputs are in the southeast part of the subbasin where characteristics of the area and proximity to the Missouri River facilitate delivery of nutrients to the Mississippi River.     

INTEGRATING PHOSPHORUS AND NITROGEN DECISION MANAGEMENT AT WATERSHED SCALES1

ABSTRACT: The persistence of water quality problems has directed attention towards the reduction of agricultural nonpoint sources of phosphorus (P) and nitrogen (N). We assessed the practical impact of three management scenarios to reduce P and N losses from a mixed land use watershed in central Pennsylvania, USA. Using Scenario 1 (an agronomic soil P threshold of 100 mg Mehlich‐3 P kg^‐1, above which no crop response is expected), 81 percent of our watershed would receive no P as fertilizer or manure. Under Scenario 2 (an environmental soil P threshold of 195 mg Mehlich‐3 P kg^‐1, above which the loss of P in surface runoff and subsurface drainage increases greatly), restricts future P inputs in only 51 percent of the watershed. Finally, using scenario 3 (P and N indices that account for likely source and transport risks), 25 percent of the watershed was at high risk or greater of P loss, while 60 percent of the watershed was classified as of high risk of nitrate (NO₃) leaching. Areas at risk of P loss were near the stream channel, while areas at risk of NO₃ leaching were near the boundaries of the watershed, where freely draining soils and high manure and fertilizer N applications coincide. Remedial measures to minimize P export should focus on critical source areas, while remedial measures to reduce N losses should be source based, concentrating on more efficient use of N by crops.     

Stream Nitrogen Sources Apportionment and Pollution Control Scheme Development in an Agricultural Watershed in Eastern China

A modeling system that couples a land-use-based export coefficient model, a stream nutrient transport equation, and Bayesian statistics was developed for stream nitrogen source apportionment. It divides a watershed into several sub-catchments, and then considers the major land-use categories as stream nitrogen sources in each sub-catchment. The runoff depth and stream water depth are considered as the major factors influencing delivery of nitrogen from land to downstream stream node within each sub-catchment. The nitrogen sources and delivery processes are lumped into several constant parameters that were calibrated using Bayesian statistics from commonly available stream monitoring and land-use datasets. This modeling system was successfully applied to total nitrogen (TN) pollution control scheme development for the ChangLe River watershed containing six sub-catchments and four land-use categories. The temporal (across months and years) and spatial (across sub-catchments and land-use categories) variability of nonpoint source (NPS) TN export to stream channels and delivery to the watershed outlet were assessed. After adjustment for in-stream TN retention, the time periods and watershed areas with disproportionately high-TN contributions to the stream were identified. Aimed at a target stream TN level of 2 mg L^?1, a quantitative TN pollution control scheme was further developed to determine which sub-catchments, which land-use categories in a sub-catchment, which time periods, and how large of NPS TN export reduction were required. This modeling system provides a powerful tool for stream nitrogen source apportionment and pollution control scheme development at the watershed scale and has only limited data requirements.     

NUTRIENT MASS BALANCE FOR THE ALBEMARLE-PAMLICO DRAINAGE BASIN,NORTH CAROLINA AND VIRGINIA, 19901

ABSTRACT: A 1990 nitrogen and phosphorus mass balance calculated for eight National Stream Quality Accounting Network (NASQAN) basins in the Albemarle-Pamlico Drainage Basin indicated the importance of agricultural nonpoint sources of nitrogen and phosphorus and watershed nitrogen retention and processing capabilities. Basin total nitrogen and phosphorus input estimates were calculated for atmospheric deposition (which averaged 27 percent of total nitrogen inputs and 22 percent of total phosphorus inputs); crop fertilizer (27 and 25 percent); animal-waste (22 and 50 percent, respectively); point sources (3 percent each of total nitrogen and total phosphorus inputs); and biological nitrogen fixation (21 percent of total nitrogen inputs). Highest in-stream nitrogen and phosphorus loads were measured in predominantly agricultural drainage areas. Intermediate loads were observed in mixed agricultural/urban drainage areas; the lowest loads were measured in mixed agricultural/forested drainage areas. The difference between the sum of the nutrient input categories and the sum of the in-stream nutrient loads and crop-harvest nutrient removal was assigned to a residual category for the basin. The residual category averaged 51 percent of total nitrogen inputs and 54 percent of total phosphorus inputs.     

Spatial Optimization of Six Conservation Practices Using Swat in Tile‐Drained Agricultural Watersheds

Targeting of agricultural conservation practices to the most effective locations in a watershed can promote wise use of conservation funds to protect surface waters from agricultural nonpoint source pollution. A spatial optimization procedure using the Soil and Water Assessment Tool was used to target six widely used conservation practices, namely no‐tillage, cereal rye cover crops (CC), filter strips (FS), grassed waterways (GW), created wetlands, and restored prairie habitats, in two west‐central Indiana watersheds. These watersheds were small, fairly flat, extensively agricultural, and heavily subsurface tile‐drained. The targeting approach was also used to evaluate the model's representation of conservation practices in cost and water quality improvement, defined as export of total nitrogen, total phosphorus, and sediment from cropped fields. FS, GW, and habitats were the most effective at improving water quality, while CC and wetlands made the greatest water quality improvement in lands with multiple existing conservation practices. Spatial optimization resulted in similar cost‐environmental benefit tradeoff curves for each watershed, with the greatest possible water quality improvement being a reduction in total pollutant loads by approximately 60%, with nitrogen reduced by 20‐30%, phosphorus by 70%, and sediment by 80‐90%.     

Source and Delivery of Nutrients to Receiving Waters in the Northeastern and Mid‐Atlantic Regions of the United States1

Moore, Richard B., Craig M. Johnston, Richard A. Smith, and Bryan Milstead, 2011. Source and Delivery of Nutrients to Receiving Waters in the Northeastern and Mid‐Atlantic Regions of the United States. Journal of the American Water Resources Association (JAWRA) 47(5):965‐990. DOI: 10.1111/j.1752‐1688.2011.00582.x Abstract: This study investigates nutrient sources and transport to receiving waters, in order to provide spatially detailed information to aid water‐resources managers concerned with eutrophication and nutrient management strategies. SPAtially Referenced Regressions On Watershed attributes (SPARROW) nutrient models were developed for the Northeastern and Mid‐Atlantic (NE US) regions of the United States to represent source conditions for the year 2002. The model developed to examine the source and delivery of nitrogen to the estuaries of nine large rivers along the NE US Seaboard indicated that agricultural sources contribute the largest percentage (37%) of the total nitrogen load delivered to the estuaries. Point sources account for 28% while atmospheric deposition accounts for 20%. A second SPARROW model was used to examine the sources and delivery of phosphorus to lakes and reservoirs throughout the NE US. The greatest attenuation of phosphorus occurred in lakes that were large relative to the size of their watershed. Model results show that, within the NE US, aquatic decay of nutrients is quite limited on an annual basis and that we especially cannot rely on natural attenuation to remove nutrients within the larger rivers nor within lakes with large watersheds relative to the size of the lake.     

A Comparative Assessment of Runoff Nitrogen from Turf,Forest, Meadow,and Mixed Landuse Watersheds

Landscaping paradigms that encourage high‐input, intensively managed and mono‐culture turf/lawn landscapes have raised concerns about water quality. We conducted a watershed‐scale assessment of landscaping practices that included turf, urban, forest, native meadow, and mixed landuse watersheds with a professional golf course and a parking lot. The turf site was moderately managed and had lower fertilizer inputs than those typically used by homeowners and golf courses. Stream water sampling was performed during base flow and storm events. Highest nitrate and total nitrogen concentrations in runoff were observed for the mixed watershed draining the golf course. In contrast, concentrations in base flow from the turf watershed were lower than expected and were comparable to those measured in the surrounding meadow and forest sites. Total nitrogen concentrations from the turf site increased sharply during the first storms following fertilization, suggesting that despite optimal management there exists a risk for nutrient runoff following fertilization. Overall, this study suggests that turf or lawns, when managed properly, pose minimal water quality risk to surface waters. Rate, timing of application, and the type of fertilizer appear to be the key factors affecting water quality. Better education of homeowners and landscaping professionals with regard to these factors may be a cost‐effective strategy to reduce nonpoint source pollution.     

Regional scale variability in sediment and nutrient delivery from small agricultural watersheds

Although many studies have pointed out the various controlling factors of sediment and nutrient delivery on a plot or watershed scale, little is known on the spatial variability of sediment and nutrient delivery on a regional scale. This study was conducted to reveal regional variations in sediment-associated nutrient delivery in central Belgium. Sediment deposited in 13 small retention ponds was sampled and analyzed for total phosphorus (TP), K, Mg, and Ca content. The TP content of the sediment deposits varied from 510 to 2001 mg P per kg sediment. Nutrients are predominantly fixed on the very fine sediment fraction (<16 microm), which is the reason why the nutrient trap efficiency of the ponds is only a fraction of the sediment trap efficiency. Average nutrient trap efficiency of the studied ponds varies between 4 and 31%, whereas sediment trap efficiency varies between 10 and 72%. For watersheds ranging from 7 to 4873 ha, sediment yield ranged between 1.2 and 20.6 Mg ha(-1) yr(-1), whereas TP export varied from 1.8 to 39.7 kg ha(-1) yr(-1). The observed spatial variability in nutrient losses is primarily attributed to regional variations in erosion and sediment yield values and to a far lesser degree to the spatial variations in fertilizer application. Redistribution of manure in the framework of an agricultural policy may increase the rate of nutrient delivery by ways of erosion and sediment transport.     

Quantitative Assessment of Agricultural Runoff and Soil Erosion Using Mathematical Modeling: Applications in the Mediterranean Region

Three mathematical models, the runoff curve number equation, the universal soil loss equation, and the mass response functions, were evaluated for predicting nonpoint source nutrient loading from agricultural watersheds of the Mediterranean region. These methodologies were applied to a catchment, the gulf of Gera Basin, that is a typical terrestrial ecosystem of the islands of the Aegean archipelago. The calibration of the model parameters was based on data from experimental plots from which edge-of-field losses of sediment, water runoff, and nutrients were measured. Special emphasis was given to the transport of dissolved and solid-phase nutrients from their sources in the farmers' fields to the outlet of the watershed in order to estimate respective attenuation rates. It was found that nonpoint nutrient loading due to surface losses was high during winter, the contribution being between 50% and 80% of the total annual nutrient losses from the terrestrial ecosystem. The good fit between simulated and experimental data supports the view that these modeling procedures should be considered as reliable and effective methodological tools in Mediterranean areas for evaluating potential control measures, such as management practices for soil and water conservation and changes in land uses, aimed at diminishing soil loss and nutrient delivery to surface waters. Furthermore, the modifications of the general mathematical formulations and the experimental values of the model parameters provided by the study can be used in further application of these methodologies in watersheds with similar characteristics.     

Toward Explaining Nitrogen and Phosphorus Trends in Chesapeake Bay Tributaries, 1992–2012

Understanding trends in stream chemistry is critical to watershed management, and often complicated by multiple contaminant sources and landscape conditions changing over varying time scales. We adapted spatially referenced regression (SPARROW) to infer causes of recent nutrient trends in Chesapeake Bay tributaries by relating observed fluxes during 1992, 2002, and 2012 to contemporary inputs and watershed conditions. The annual flow‐normalized nitrogen flux to the bay from its watershed declined by 14% to 127,000 Mg (metric tons) between 1992 and 2012, due primarily (more than 80% of the decline) to reduced point sources. The remainder of the decline was due to reduced atmospheric deposition (13%) and urban nonpoint sources. Agricultural inputs, which contribute most nitrogen to the bay, changed little, although trends in the average nitrogen yield (flux per unit area) from cropland and pasture to streams in some settings suggest possible effects of evolving nutrient applications or other land management practices. Point sources of phosphorus to local streams declined by half between 1992 and 2012, while nonpoint inputs were relatively unchanged. Annual phosphorus delivery to the bay increased by 9% to 9,570 Mg between 1992 and 2012, however, due mainly to reduced retention in the Susquehanna River at Conowingo Reservoir.     

Modeling riverine nitrate export from an East-Central Illinois watershed using SWAT

Reliable water quality models are needed to forecast the water quality consequences of different agricultural nutrient management scenarios. In this study, the Soil and Water Assessment Tool (SWAT), version 2000, was applied to simulate streamflow, riverine nitrate (NO(3)) export, crop yield, and watershed nitrogen (N) budgets in the upper Embarras River (UER) watershed in east-central Illinois, which has extensive maize-soybean cultivation, large N fertilizer input, and extensive tile drainage. During the calibration (1994-2002) and validation (1985-1993) periods, SWAT simulated monthly and annual stream flows with Nash-Sutcliffe coefficients (E) ranging from 0.67 to 0.94 and R(2) from 0.75 to 0.95. For monthly and annual NO(3) loads, E ranged from -0.16 to 0.45 and R(2) from 0.36 to 0.74. Annual maize and soybean yields were simulated with relative errors ranging from -10 to 6%. The model was then used to predict the changes in NO(3) output with N fertilizer application rates 10 to 50% lower than original application rates in UER. The calibrated SWAT predicted a 10 to 43% decrease in NO(3) export from UER and a 6 to 38% reduction in maize yield in response to the reduction in N fertilizer. The SWAT model markedly overestimated NO(3) export during major wet periods. Moreover, SWAT estimated soybean N fixation rates considerably greater than literature values, and some simulated changes in the N cycle in response to fertilizer reduction seemed to be unrealistic. Improving these aspects of SWAT could lead to more reliable predictions in the water quality outcomes of nutrient management practices in tile-drained watersheds.     

Application of Tabu Search Algorithm With a Coupled AnnAGNPS‐CCHE1D Model to Optimize Agricultural Land Use1

Abstract: A principal contributor to soil erosion and nonpoint source pollution, agricultural activities have a major influence on the environmental quality of a watershed. Impact of agricultural activities on the quality of water resources can be minimized by implementing suitable agriculture land‐use types. Currently, land uses are designed (location, type, and operational schedule) based on field study results, and do not involve a science‐based approach to ensure their efficiency under particular regional, climatic, geological, and economical conditions. At present, there is a real need for new methodologies that can optimize the selection, design, and operation of agricultural land uses at the watershed scale by taking into account environmental, technical, and economical considerations, based on realistic simulations of watershed response. In this respect, the present study proposes a new approach, which integrates computational modeling of watershed processes, fluvial processes in the drainage network, and modern heuristic optimization techniques to design cost effective land‐use plans. The watershed model AnnAGNPS and the channel network model CCHE1D are linked together to simulate the sediment and pollutant transport processes. Based on the computational results, a multi‐objective function is set up to minimize soil losses, nutrient yields, and total associated costs, while the production profits from agriculture are maximized. The selected iterative optimization algorithm uses adaptive Tabu Search heuristic to flip (switching from one alternative to another) land‐change variables. USDA’s Goodwin Creek experimental watershed, located in Northern Mississippi, is used to demonstrate the capabilities of the proposed approach. The results show that the optimized land‐use design with BMPs using an integrated approach at the watershed level can provide efficient and cost‐effective conservation of the environmental quality by taking into account both productivity and profitability.     

Applying Statistical Causal Analyses to Agricultural Conservation: A Case Study Examining P Loss Impacts

Estimating the effect of agricultural conservation practices on reducing nutrient loss using observational data can be confounded by factors such as differing crop types and management practices. As we may not have the full knowledge of these confounding factors, conventional statistical meta‐analysis methods can be misleading. We discuss the use of two statistical causal analysis methods for quantifying the effects of water and soil conservation practices in reducing P loss from agricultural fields. With the propensity score method, a subset of data was used to form a treatment group and a control group with similar distributions of confounding factors. With the multilevel modeling method, data were stratified based on important confounding factors, and the conservation practice effect was evaluated for each stratum. Both methods resulted in similar estimates of the conservation practice effect (total P load reduction avg. ~70%). In addition, both methods show evidence of conservation practices reducing the incremental increase in total P export per unit increase in fertilizer application. These results are presented as examples of the types of outcomes provided by statistical causal analyses, not to provide definitive estimates of P loss reduction. The enhanced meta‐analysis methods presented within are applicable for improved assessment of agricultural practices and their effects and can be used for providing realistic parameter values for watershed‐scale modeling.     

WATER QUALITY OF AN INTENSWE AGRICULTURAL WATERSHED IN QUEBEC1

ABSTRACT: . Under a watershed based approach being examined by the Quebec Ministry of Agriculture to accelerate the adoption of conservation practices, a study on the impacts of agricultural practices on the St. Esprit watershed was initiated in the fall of 1993. The water quality of this 26 km² intensive agricultural watershed was studied over an 18 month period. Water samples taken at the outlet of the watershed were analyzed for nitrate, phosphate, suspended sediment, and atrazine. Water quality data were analyzed to establish seasonal trends in pollutant concentration and load in the watercourse. Spring snowmelt was identified as a significant period of pollutant material export. All pollutant materials displayed seasonal variability in the export process. Peak pollutant concentrations were associated with high flow events. Mean observed pollutant concentrations did not exceed drinking water quality standards.     

Application of SPARROW Modeling to Understanding Contaminant Fate and Transport from Uplands to Streams

Understanding spatial variability in contaminant fate and transport is critical to efficient regional water‐quality restoration. An approach to capitalize on previously calibrated spatially referenced regression (SPARROW) models to improve the understanding of contaminant fate and transport was developed and applied to the case of nitrogen in the 166,000 km² Chesapeake Bay watershed. A continuous function of four hydrogeologic, soil, and other landscape properties significant (α = 0.10) to nitrogen transport from uplands to streams was evaluated and compared among each of the more than 80,000 individual catchments (mean area, 2.1 km²) in the watershed. Budgets (including inputs, losses or net change in storage in uplands and stream corridors, and delivery to tidal waters) were also estimated for nitrogen applied to these catchments from selected upland sources. Most (81%) of such inputs are removed, retained, or otherwise processed in uplands rather than transported to surface waters. Combining SPARROW results with previous budget estimates suggests 55% of this processing is attributable to denitrification, 23% to crop or timber harvest, and 6% to volatilization. Remaining upland inputs represent a net annual increase in landscape storage in soils or biomass exceeding 10 kg per hectare in some areas. Such insights are important for planning watershed restoration and for improving future watershed models.     

Evaluating agricultural best management practices in tile-drained subwatersheds of the Mackinaw River, Illinois

Best management practices (BMPs) are widely promoted in agricultural watersheds as a means of improving water quality and ameliorating altered hydrology. We used a paired watershed approach to evaluate whether focused outreach could increase BMP implementation rates and whether BMPs could induce watershed-scale (4000 ha) changes in nutrients, suspended sediment concentrations, or hydrology in an agricultural watershed in central Illinois. Land use was >90% row crop agriculture with extensive subsurface tile drainage. Outreach successfully increased BMP implementation rates for grassed waterways, stream buffers, and strip-tillage within the treatment watershed, which are designed to reduce surface runoff and soil erosion. No significant changes in nitrate-nitrogen (NO-N), total phosphorus (TP), dissolved reactive phosphorus, total suspended sediment (TSS), or hydrology were observed after implementation of these BMPs over 7 yr of monitoring. Annual NO-N export (39-299 Mg) in the two watersheds was equally exported during baseflow and stormflow. Mean annual TP export was similar between the watersheds (3.8 Mg) and was greater for TSS in the treatment (1626 ± 497 Mg) than in the reference (940 ± 327 Mg) watershed. Export of TP and TSS was primarily due to stormflow (>85%). Results suggest that the BMPs established during this study were not adequate to override nutrient export from subsurface drainage tiles. Conservation planning in tile-drained agricultural watersheds will require a combination of surface-water BMPs and conservation practices that intercept and retain subsurface agricultural runoff. Our study emphasizes the need to measure conservation outcomes and not just implementation rates of conservation practices.     

Weighting Nitrogen and Phosphorus Pixel Pollutant Loads to Represent Runoff and Buffering Likelihoods

Watershed models often estimate annual nitrogen (N) or phosphorus (P) pollutant loads in rural areas with export coefficient (EC) (kg/ha/yr) values based on land cover, and in urban areas as the product of spatially uniform event mean concentration (EMC) (mg/L) values and runoff volume. Actual N and P nonpoint source (NPS) pollutant loading has more spatial complexity due to watershed variation in runoff likelihood and buffering likelihood along surface and subsurface pathways, which can be represented in a contributing area dispersal area (CADA) NPS model. This research develops a CADA NPS model to simulate how watershed properties of elevation, land cover, and soils upslope and downslope of each watershed pixel influence nutrient loading. The model uses both surface and subsurface runoff indices (RI), and surface and subsurface buffer indices (BI), to quantify the runoff and buffering likelihood for each watershed pixel, and generate maps of weighted EC and EMC values that identify NPS pollutant loading hotspots. The research illustrates how CADA NPS model maps and pixel loading values are sensitive to the spatial resolution and accuracy of elevation and land cover data, and model predictions can represent the lower and upper bounds of NPS loading. The model provides managers with a tool to rapidly visualize, rank, and investigate likely areas of high nutrient export.     

Quantifying Legacy Phosphorus Using a Mass Balance Approach and Uncertainty Analysis

Classic agricultural‐conservation practices may not address decades of phosphorus (P) accumulation, known as legacy P. Identifying and quantifying legacy P sources are necessary to identify the most cost‐efficient conservation practices. A method was developed to identify and quantify legacy P at the watershed scale using a mass‐balance approach and uncertainty analysis. The method was applied to two nutrient‐rich watersheds in northeast Oklahoma and northwest Arkansas. Each P import and export to and from the two watersheds was identified and quantified using a probability distribution and uncertainty analysis. The P retained in the soils, reservoirs, and stream systems were estimated from 1925 to 2015. Over 8.5 and 6.1 kg/ha/year of P were added to the Illinois River and Eucha‐Spavinaw watersheds with 53 and 55% from poultry production, respectively. Other major historical sources were attributed to human population and commercial fertilizer. Though currently the net addition of P in the watersheds is small due to the export of approximately 90% of the poultry litter, historically only 14‐19% of all P imported to the Illinois River and Eucha‐Spavinaw watersheds was removed via the reservoir spillways, poultry litter, and food exports. The majority of the retained P is located in the soil, 3.6‐5.8 kg/ha/year, and stream systems, 0.01‐3.0 per ha/year.