Phosphorus transport in agricultural subsurface drainage: a review

Denitrification in wetlands: A review towards a quantification at global scale

Research to understand the nitrogen cycle has been thriving. The production of reactive nitrogen by humans exceeds the removal capacity through denitrification of any natural ecosystem. The surplus of reactive nitrogen is also a significant pollutant that can shift biological diversity and distribution, promotes eutrophication in aquatic ecosystems, and affects human health. Denitrification is the microbial respiration in anoxic conditions and is the main process that removes definitively nitrates from the ecosystem by returning of reactive nitrogen (Nr) to the atmosphere as N₂ and N₂O emissions. This process occurs in the oceans, aquatic ecosystems and temporary flooded terrestrial ecosystems. Wetlands ecosystems are rich in organic matter and they have regular anoxic soil conditions ideal for denitrification to occur. In the current paper, we provide a meta-analysis that aims at exploring how research around global nitrogen, denitrification and wetlands had evolved in the last fifty years. Back in the time, wetland ecosystems were seen as non-exploitable elements of the landscape, and now they are being integrated as providers of ecosystem services. A significant improvement of molecular biology techniques and genetic extraction have made the denitrification process fully understood allowing constructed wetlands to be more efficient and popular. Yet, large uncertainties remain concerning the dynamic quantification of the global denitrification capacity of natural wetland ecosystems. The contribution of the current investigation is to provide a way forward for reducing these uncertainties by the integration of satellite-based Earth Observation (EO) technology with parsimonious physical based models.

Cultivar-Dependent Responses in Plant Growth, Leaf Physiology, Phosphorus Use Efficiency, and Tuber Quality of Potatoes Under Limited Phosphorus Availability Conditions

The limited availability of phosphorus (P) in soils causes a major constraint in the productivity of potatoes, which requires increased knowledge of plant adaptation responses in this condition. In this study, six potato cultivars, namely, Agria, Lady Claire, Milva, Lilly, Sieglinde, and Verdi, were assessed for their responses on plant growth, leaf physiology, P use efficiency (PUE), and tuber quality with three P levels (P_low, P_med, and P_high). The results reveal a significant variation in the cultivars in response to different P availabilities. P-efficient cultivars, Agria, Milva, and Lilly, possessed substantial plant biomass, tuber yield, and high P uptake efficiency (PUpE) under low P supply conditions. The P-inefficient cultivars, Lady Claire, Sieglinde, and Verdi, could not produce tubers under P deprivation conditions, as well as the ability to efficiently uptake P under low-level conditions, but they were efficient in P uptake under high soil P conditions. Improved PUpE is important for plant tolerance with limited P availability, which results in the efficient use of the applied P. At the leaf level, increased accumulations of nitrate, sulfate, sucrose, and proline are necessary for a plant to acclimate to P deficiency-induced stress and to mobilize leaf inorganic phosphate to increase internal PUE and photosynthesis. The reduction in plant biomass and tuber yield under P-deficient conditions could be caused by reduced CO₂ assimilation. Furthermore, P deficiency significantly reduced tuber yield, dry matter, and starch concentration in Agria, Milva, and Lilly. However, contents of tuber protein, sugars, and minerals, as well as antioxidant capacity, were enhanced under these conditions in these cultivars. These results highlight the important traits contributing to potato plant tolerance under P-deficient conditions and indicate an opportunity to improve the P efficiency and tuber quality of potatoes under deficient conditions using more efficient cultivars. Future research to evaluate molecular mechanisms related to P and sucrose translocation, and minimize tuber yield reduction under limited P availability conditions is necessary.

High flow event induced the subsurface transport of particulate phosphorus and its speciation in agricultural tile drainage system

Subsurface storm flow of phosphorus (P), including particulate P, has been recently discussed as an important P transport path in contrast to typical surface runoff events. However, P speciation, and P concentration during storm events has not been extensively investigated; therefore, its contribution to the water quality is not clearly understood. In this study, the physicochemical properties of particulate P in tile water samples during a high flow event were investigated in Midwestern agricultural lands using wet chemical methods, 31P Nuclear Magnetic Resonance spectroscopy and P K-edge X-ray absorptions near edge structure spectroscopy. In slightly alkaline pH tile water, total P was ranging from ∼0.06 to 0.22 mg L^-1, which is significantly greater than dissolved reactive P (DRP) (∼0.02-0.08 mg L^-1). The tile water contains P enriched particulate matters (∼200-660 mg L^-1). Total P in the colloidal fraction was from 1013 to 2270 mg kg^-1. Phosphate and organic P species, especially monoesters, are sorbed in soil colloids like calcite, and iron oxides, and colloids are effective carriers of P in the subsurface transport process during storm events. The results of this study show that storm events can accelerate the subsurface transport of P with soil particles in addition to DRP.

Surface runoff and nutrient dynamics in cover crop-soybean systems in the Upper Midwest

Relay-cropping of the novel oilseeds winter camelina (Camelina sativa L.) and pennycress (Thlaspi arvense L.) with short-season crops such as soybean [Glycine max (L.) Merr.] can provide economic and environmental incentives for adopting winter cover crop practices in the U.S. Upper Midwest. However, their ability to reduce nutrient loss in surface runoff is unknown. Accordingly, surface runoff and quality were evaluated during three seasonal phases (cover, intercrop, and soybean) over 2 yr in four cover crop-soybean treatments (pennycress, winter camelina, forage radish [Raphanus sativus L.], and winter rye [Secale cereale L.]) compared with no-till and chisel-till fallow treatments. Runoff was collected with Gerlach troughs and assessed for concentrations and loads of NO3 - -N, total mineral N, soluble reactive P (SRP), and total suspended solids (TSS). Cumulative runoff and nutrient loads were greater during the winter cover phase because of increased snow melt and freeze-thaw released nutrients from living vegetation. In contrast, cumulative TSS was greater during intercrop and soybean phases due to high-intensity rainfall events with an open soybean canopy. Average TSS loads during the intercrop phase were reduced by 75% in pennycress compared with fallow and radish treatments. During the soybean phase, average TSS, total mineral N, and SRP loads were generally elevated in cover crop treatments compared with no-till. Overwintering cover crops may contribute to mobility of nutrients solubilized from living or decomposing vegetation; however, this was balanced by their potential to reduce runoff and TSS during high-intensity spring rains.

Phosphorus transport in riverbed sediments and related adsorption and desorption characteristics in the Beiyun River, China

Riverbed sediments are the interface layer in riverine ecosystems connecting the overlying medium of water and the vadose zone. The transport behavior of phosphorus (P), which has been recognized as the primary cause of freshwater eutrophication, in riverbed sediments remains unclear. Understanding the impact of riverbed sediments on P transport is a necessary prerequisite for the development of appropriate strategies to reduce potential groundwater pollution. In this study, riverbed sediments were collected from the upstream, midstream, and downstream sections of the Beiyun River, China, and packed into vertical soil columns to perform leaching experiments to quantify P transport characteristics. In addition, the impact mechanisms were further explored by conducting laboratory batch tests of P adsorption and desorption. The results demonstrated that approximately 80% of P can be adsorbed by riverbed sediments in soil column leaching experiment, and a tailing phenomenon was observed in its desorption. The hydraulic conductivity properties of riverbed sediments were evaluated by the advection-dispersion equation, showing a gradually decreasing adsorption capacity for P from upstream to downstream sections, which was supported by the results obtained from adsorption-desorption thermodynamic and kinetic batch tests. The estimated annual leaching masses of P increased from 60.72 g/(m² a) in the upstream section to 132.31 g/(m² a) in the downstream section. The role of riverbed sediments as a source or sink of P is possibly determined by their coarse sand particles content, and the mean equilibrium P concentration (EPC0). The competitive relationship between P and other forms of nutrients also has an important influence on its source-sink role. These findings suggest that the prevention of the potential P leaching is most needed in the downstream sections of Beiyun River, and corresponding control strategies should be developed to avoid groundwater pollution.

Riparian buffer effectiveness as a function of buffer design and input loads

Although many agricultural watersheds rely heavily on riparian buffer adoption to meet water quality goals, design and management constraints in current policies create adoption barriers. Based on focus group feedback, we developed a flexible buffer design paradigm that varies buffer width, vegetation, and harvesting. Sixteen years of daily-scale nutrient and sediment loads simulated with the Soil and Water Assessment Tool (SWAT) were coupled to the three-zone Riparian Ecosystem Management Model (REMM) to compare the effectiveness of traditional, policy-based buffer designs with designs that are more flexible and integrate features important to local farmers. Buffer designs included (i) 10 m grass, (ii) 15 m grass, (iii) 15 m deciduous trees, (iv) 30 m grass and trees, (v) 30 m grass and trees with trees harvested every 3 yr, and (vi) 30 m grass and trees with grass harvested every year. Allowing harvesting in one zone of the buffer vegetation (either trees or grasses) minimally affected water quality, with annual average percent reductions differing by <5% (p > .05; 76-78% for total nitrogen [TN], 51-55% for total phosphorus [TP], and 68% for sediment). Under the highest input loading conditions, buffers with lower removal efficiencies removed more total mass than did buffers with high removal efficiencies. Thus, by focusing on mass reduction in addition to percent reduction, watershed-wide buffer implementation may be better targeted to TN, TP, and sediment reduced. These findings have important implications for informing flexible buffer design policies and enhanced placement of buffers in watersheds impaired by nutrient and sediment.

The Shallow and Deep Hypothesis: Subsurface Vertical Chemical Contrasts Shape Nitrate Export Patterns from Different Land Uses

Eutrophication has threatened water resources worldwide, yet mechanistic understanding on controls of nutrient export remains elusive. This work tests the shallow and deep hypothesis: subsurface vertical chemical contrasts regulate nitrate export patterns under different land use conditions. We synthesized data from 228 watersheds and used reactive transport modeling (500 simulations) under broad land use, climate, and geology conditions. Data synthesis indicated that human perturbation has amplified chemical contrasts in shallow water (e.g., soil water) versus deep waters (e.g., groundwater), inducing primarily flushing patterns (concentrations increase with streamflow) in agriculture lands and dilution patterns (concentrations decrease with streamflow) in urban watersheds. Results revealed a quantitative relationship between export patterns and shallow-versus-deep concentration contrasts, underscoring the often-overlooked role of nutrient distribution over depth. Results challenge the commonly held perception that legacy stores in agricultural lands induce chemostasis where concentrations vary negligibly with streamflow. They suggest that nitrate concentrations from agricultural lands will escalate during large hydrological events, which can exacerbate nutrient export problems as flooding events intensify in the future climate.

Effects of cropping system and fertilization regime on soil phosphorous are mediated by rhizosphere-microbial processes in a semi-arid agroecosystem

In semi-arid regions, soil phosphorus (P) dynamics in cereal-legume intercropping are not yet fully elucidated, particularly in relation to integrated application of fertilizers. To this aim, we investigate the effects of different fertilizers on various P fractions in relation to the rhizosphere-microbial processes in a cowpea/maize intercropping system. Field experiments were conducted during two consecutive years (2016-2017) in a split-plot design by establishing cowpea/maize alone or intercropped onto the main plot, while the sub-plot was treated with four types of fertilization, i.e. no fertilizer addition (control), organic amendment (compost), mineral fertilizers (NPK) and multi-nutrient enriched compost (NPKEC). Our results showed that NPKEC fertilizer increased NaHCO₃-P_i by 69% in maize, 62% in cowpea and 93% in intercropped plots compared to control plots. Similarly, a significant increase in the NaHCO₃-P_o fraction was also recorded with NPKEC treatment in all cropping systems. In case of moderately labile P, NPKEC fertilizer caused the highest increase of NaOH-P_o (12.87 ± 0.50 mg P kg^-1 soil) and NaOH-P_i (22.29 ± 0.83 mg P kg^-1 soil) fractions in intercropped plots. Except for intercropping, NPK application caused an increase in the non-available P fraction (HCl-P_i), while the use of NPKEC decreased the HCl-P_i concentration in all cropping systems, suggesting stronger merits both for intercropping and NPKEC. Surprisingly, maize exhibited substantially higher phosphatases activity compared to cowpea in monoculture amended with compost, implying distinct crop strategies for adaptation under low P conditions. Based on the multi-factor analysis, the close association of NaHCO₃-P with P solubilizing bacteria, root carboxylates and pH indicated that rhizosphere processes are the strongest predictors of immediately available P. Since alkaline phosphatase (ALP) is a P-degrading enzyme of microbial origin, rhizosphere related ALP association may have originated from root-associated microflora promoting P mobilization. Furthermore, the strong association of microbial biomass P (MBP) and acid phosphates (ACP) with NaOH-P fraction indicated moderately available P cycle in soil was mainly driven by microbial-related processes. Factor analysis map and two-way ANOVA confirmed that fertilization regime had a stronger effect on all tested variables compared to cropping system. Altogether, our results suggest that a combination of microbial-rhizosphere processes controls the dynamics of P fertility in semi-arid soils. In the broader context of improving soil P fertility, it is highly recommended the use of environmentally sustainable sources of fertilizer, such as NPKEC, which can enhance the competitive performance of legume-cereal intercropping under semi-arid agroecosystems.

Sources and subsurface transport of dissolved reactive phosphorus in a semiarid, no-till catchment with complex topography

The subsurface transport of dissolved reactive phosphorus (DRP) from artificially drained agricultural fields can impair water quality, especially in no-till fields. The distribution of soil P in the wheat (Triticum aestivum L.)-dominated Palouse region in the inland U.S. Pacific Northwest varies greatly due to its steep and complex topography, and a legacy (∼130 yr) of excessive soil erosion and deposition processes. The primary goal of this research was to better understand the magnitude and temporal dynamics of DRP export from an artificial drain line and the variability of subsurface DRP leaching within a long-term, no-till field. Dissolved reactive P in drain line effluent was monitored across three water years. Large intact soil cores were extracted at contrasting field locations (toe and top slope positions) to measure DRP leachate concentration and relative P sorption. Drain line DRP concentration was predominantly >0.05 mg L^-1 and often exceeded 0.1 mg L^-1 during winter and early spring. Mean leachate DRP levels were significantly higher in toe slope cores than in top slope cores (0.11 and 0.02 mg L^-1 , respectively). Saturated hydraulic conductivity varied widely across cores and was not correlated with leachate DRP concentration. All soil cores exhibited high P sorption potential, even under conditions of preferential flow. These findings suggest that much of the DRP transport in these landscapes is derived from P hotspots located in toe slope positions. Application of soil P fertilizer amounts in variable rates that account for spatial variability in P transport may minimize P enrichment and subsequent leaching in these locations.

Modeling of phosphorus loss from field to watershed: A review

Phosphorus (P) losses from nonpoint sources into surface water resources through surface runoff and tile drainage play a significant role in eutrophication. Accordingly, the number of studies involving the modeling of agricultural P losses, the uncertainties of such models, and the best management practices (BMPs) supported by the modeling of hypothetical P loss reduction scenarios has increased significantly around the world. Many improvements have been made to these models: separate manure P pools, variable source areas allowing the determination of critical source areas of P loss, analyses of modeling uncertainties, and understanding of legacy P. However, several elements are still missing or have yet to be sufficiently addressed: the incorporation of preferential flow into models, the modification of P sorption-desorption processes considering recent research data (e.g., pedotransfer functions for labile, active, or stable P, along with P sorption coefficients), BMP parameterization, and scale-up issues, as well as stakeholder-scientist and experimentalist-modeler interactions. The accuracy of P loss modeling can be improved by (a) incorporating dynamic P sorption-desorption processes and new P subroutines for direct P loss from manure, fertilizer, and dung, (b) modeling preferential flow, connectivity between field and adjacent water bodies, and P in-stream processes, (c) including an assessment of model uncertainty, (d) integrating field and watershed models for BMP calibration and scaling field results up to larger areas, and (e) building a holistic interaction between stakeholders, experimentalists, and modelers.

Closed depressions and soil phosphorus influence subsurface phosphorus losses in a tile-drained field in Illinois

Artificial subsurface (tile) drainage systems can convey phosphorus (P) from agricultural fields to surface waters; however, controls of subsurface dissolved reactive P (DRP) losses at the sub-field scale are not fully understood. We characterized subsurface DRP loads and flow-weighted mean concentration (FWMC) from January 2015 through September 2017 to determine seasonal (growing vs. non-growing) patterns from 36 individually monitored plots across a farm under a corn (Zea mays L.) and soybean [Glycine max (L.) Merr.] rotation in east-central Illinois. Using linear mixed models, we investigated the effects of soil test P (STP), depression depth, and their interaction with precipitation and P fertilization on subsurface DRP losses. Dissolved reactive P loads in drainage tiles increased with precipitation and were greatest during the non-growing season (NGS) in 2016 and 2017. Annual subsurface DRP loads were positively related to STP, and during the NGS, there was a positive relationship between depression depth quantified at the plot-scale and subsurface DRP loads and FWMC. Along a depression-depth gradient, piecewise regression displayed a threshold at a depth of 0.38 m at which STP increased, indicating soil P accumulation in deeper closed depressions. Our study highlights the need to identify areas with the greatest risk of subsurface P losses to implement sub-field scale nutrient management practices.

Impacts of preferential flow and agroecosystem management on subsurface particulate phosphorus loadings in tile-drained landscapes

Recent research on tile-drainage has placed emphasis on dissolved reactive phosphorus (DRP) delivery and transport pathways but less emphasis on particulate P (PP), resulting in its exclusion from agricultural water management models. In this study, we quantified the fluxes, mechanisms, and factors driving PP delivery into tiles through statistical analysis of a long-term hydrologic and water quality dataset. The dataset includes 5 yr of surface and tile discharge, total P (TP), DRP, total nitrogen (TN), and dissolved inorganic N concentrations from two edge-of-field study sites with contrasting soil and management practices. Hydrograph recession techniques were coupled with multiple linear regression for understanding hydrologic flow pathways, and empirical mode decomposition (EMD) time-series analysis was used to determine the significance of PP seasonality processes and the effect of management practices. The analysis of hydrologic flow pathways demonstrated that quickflow contributed 66 and 36% of subsurface discharge in the clay and loam sites, respectively. Phosphorus loading analysis showed that macropore flow plays a significant role in PP delivery to subsurface P loading and that PP significantly contributed to TP and DRP delivery; however, greater PP loadings were observed at the clay site despite greater subsurface discharge and soil test P levels at the loam site. Furthermore, PP delivery was significantly affected by environmental conditions and management practices. We highlight the efficacy of hydrograph recession analysis for identifying macropore and diffuse drainage, of P/N ratios to characterize sediment delivery mechanisms in tiles, and of EMD to detect management impacts on TP and DRP at the field scale.

Performance of a Ditch-Style Phosphorus Removal Structure for Treating Agricultural Drainage Water with Aluminum-Treated Steel Slag

Several structural, treatment, and management approaches exist to minimize phosphorus (P) transport from agricultural landscapes (e.g., cover cropping and conservation tillage). However, many of these practices are designed to minimize particulate P transport and are not as effective in controlling dissolved P (DP) losses. Phosphorus removal structures employ a P sorption material (PSM) to trap DP from flowing water. These structures have been successful in treating surface runoff by utilizing aluminum (Al)-treated steel slag, but subsurface tile drainage has never been tested with this material. The goal of this study was to evaluate the performance and economics of a ditch-style P removal structure using Al-treated steel slag for treating agricultural subsurface drainage discharge. The structure treated subsurface drainage water from a 4.5 ha agricultural field with elevated soil test P levels. Overall, the structure removed approximately 27% and 50% of all DP and total P (TP) entering the structure, respectively (i.e., 2.4 and 9.4 kg DP and TP removal). After an initial period of strong DP removal, the discrete DP removal became highly variable. Flow-through analysis of slag samples showed that the slag used to construct the structure was coarser and less sorptive compared to the slag samples collected prior to construction that were used to design and size the structure. Results of this study highlight the importance of correctly designing the P removal structures using representative PSMs.

Variations in foliar carbon:nitrogen and nitrogen:phosphorus ratios under global change: a meta-analysis of experimental field studies

Foliar-level stoichiometry plays an important role in ecosystem elemental cycling. Shifts in foliar ratios of carbon to nitrogen (C:N) and nitrogen to phosphorus (N:P) in response to global change can therefore have a large impact upon ecosystem function. We conducted a meta-analysis with 2,236 paired observations from 123 published studies to investigate the responses of foliar C:N and N:P ratios to experimental global change treatments, i.e. warming, increased precipitation, drought, N addition and elevated carbon dioxide concentration (eCO2), in field conditions. Foliar C:N and N:P ratios were neither affected by warming nor by increased precipitation. Foliar C:N ratio increased with drought and eCO2, and decreased with N addition. Foliar N:P ratios declined with eCO2, and increased under drought and N addition. Our results suggested the responses of the C:N ratio to global change were mainly related to shifts in foliar [N], whereas changes in the N:P ratio were related to the responses of both [N] and [P]. Moreover, the response magnitude of foliar N:P ratio decreased with treatment duration under increased precipitation, N addition and eCO2. Our findings are important for our understanding of plant nutrient dynamic and modeling of nutrient biogeochemistry under global change.

Nutrient loss from floodplain soil with controlled subsurface drainage under forage production

Expansion of subsurface drainage into forage production may have a deleterious effect on surface waters due to increased nitrogen and phosphorus loading. The impact of controlled subsurface drainage (CD) on nitrogen and phosphorus loss compared with free subsurface drainage (FD) in tile drainage water has been explored to a lesser extent from forage production systems. This study quantifies the effects of CD and FD on average seasonal concentrations and cumulative loads of the total suspended solids (TSS), nitrate nitrogen (NO₃ -N), and dissolved reactive phosphorus (DRP) in subsurface drainage water from a poorly drained floodplain soil in a cereal rye (Secale cereale L.)-sorghum [Sorghum bicolor (L.) Moench] rotation with rotational cattle grazing. During all crop seasons of sorghum production (2010-2013), CD had 6.03-9.63 mg L^-1 less NO₃ -N than FD. Mean DRP concentration was significantly higher for CD than for FD during all seasons except for sorghum in 2012-2013. Average cumulative discharge was 38 and 314 m³ ha^-1 less for CD than for FD during sorghum and cereal rye growing seasons, respectively. Controlled drainage had 0.68-6.14 kg ha^-1 lower cumulative NO₃ -N loads than FD. The DRP loads were dependent on discharge. During sorghum growing seasons, TSS and DRP loads were 79-90% lower in CD compared with FD. The ability to reduce drainage water flow from tiles and subsequent nitrogen and phosphorus loading with CD compared with FD in a floodplain soil indicates that CD can be effective best management practice for forage production systems.

Drivers of excess phosphorus and stream sediments in a nested agricultural catchment during base and stormflow conditions

A variety of landscape and hydrological characteristics influence nutrient concentrations and suspended sediments in freshwater systems, yet the combined influence of these characteristics within nested agricultural catchments is still poorly understood, particularly across varying flow states. To tease apart potential drivers at within-catchment scales, it is necessary to sample at a spatiotemporal resolution that captures how landscape drivers change with time. The overall objective of this study was to evaluate the relative influence of landscape and hydrological characteristics at sub-catchment scales in relation to total P (TP), soluble reactive P (SRP), the ratio of SRP and TP (SRP/TP), and total suspended solids (TSS) across varying flow conditions. Synoptic surveys were conducted at 13 longitudinal sampling sites under a variety of flow conditions (n = 14) between 2016 and 2017 in the Innisfil Creek watershed, southern Ontario. The surveys were grouped into baseflow and stormflow conditions, and partial least squares regression (PLSR) was used to characterize the relationships between catchment characteristics, median concentrations of P, and TSS. Soil texture (i.e., clay dominated), winter wheat (Triticum aestivum L.), and constructed drain density had the largest influences on stormflow SRP and SRP/TP ratios, but measures of soil erosion, like the Bank Erosion Hazard Index and sinuosity, had the largest influence on stormflow TSS. During baseflow periods, these landscape characteristics were not informative, and they were difficult to tie to in-stream conditions. Overall, our PLSR models indicated that buried tile drainage was a major source of SRP in Innisfil Creek, whereas bank erosion was a dominant source of TSS.

Valorization of farm pond biomass as fertilizer for reducing basin-scale phosphorus losses

Long-term fertilizer phosphorus (P) inputs are causing phosphorous saturation of agricultural soils globally. The saturation is spreading to the edge-of-the-farm stormwater detention systems (SDSs) from where the legacy P is potentially being released to downstream surface waters. We use site-specific and literature data for P-saturated SDSs, to develop and evaluate the biogeochemical and economic feasibility of a P recycling program that targets both low (LIC, sugarcane) and high intensity cropping (HIC, fresh-produce) systems within a watershed. The focus is to close the P cycle loop to rejuvenate P sink function of SDSs. It involves harvesting and composting the SDS's biomass and it's on-farm use as an organic fertilizer for crops. Results showed that harvesting-composting can conservatively increase the P retention from 50% to 77% for HIC and almost complete treatment for LIC. Beyond potentially increasing yield and improving soil health, compost use can further increase in-field retention of P (and water). Additional costs incurred in harvesting and composting can be offset by the economic value of compost and the reduction in State's expenditure on regional P treatment systems. Treatment costs were $26/kg of P for HIC and $42/kg for LIC, 10 times less than the current state expenditure of $355-$909/kg P using constructed wetlands. We propose an incentivized, payment for services (PS) program, where producers are paid for P recycling. The PS program considers the intensity of cropping systems and their location along the drainage network from headwaters to the outlet, to achieve basin-scale P load reduction. The LIC SDSs recover regional P by passing the public water through them while recycling is implemented at the HIC. The estimated basin-scale P retention with harvest-compost approach was 854 metric tons, 5 times the P that entered the Everglades Protection Area in 2018, at 88%-93% less cost than the State treatment systems.

Temperature and freezing effects on phosphorus release from soils to overlying floodwater under flooded-anaerobic conditions

Increased phosphorus (P) availability under flooded, anaerobic conditions may accelerate P loss from soils to water bodies. Existing knowledge on P release to floodwater from flooded soils is limited to summer conditions and/or room temperatures. Spring snowmelt runoff, which occurs under cold temperatures with frequent freeze-thaw events, is the dominant mode of P loss from agricultural lands to water bodies in the Canadian Prairies. This research examined the effects of temperature on P dynamics under flooded conditions in a laboratory study using five agricultural soils from Manitoba, Canada. The treatments were (a) freezing for 1 wk at -20 °C, thawing and flooding at 4 ± 1 °C (frozen, cold); (b) flooding unfrozen soil at 4 ± 1 °C (unfrozen, cold); and (c) flooding unfrozen soil at 20 ± 2 °C (warm). Pore water and surface water were collected weekly over 8 wk and analyzed for dissolved reactive phosphorus (DRP), pH, calcium, magnesium, iron (Fe), and manganese (Mn). Soils under warm flooding showed enhanced P release with significantly higher DRP concentrations in pore and surface floodwater compared with cold flooding of frozen and unfrozen soils. The development of anaerobic conditions was slow under cold flooding with only a slight decrease in Eh, whereas under warm flooding Eh declined sharply, favoring reductive dissolution reactions releasing P, Fe, and Mn. Pore water and floodwater DRP concentrations were similar between frozen and unfrozen soil under cold flooding, suggesting that one freeze-thaw event prior to flooding had minimal effect on P release under simulated snowmelt conditions.

Nutrient control in water bodies: A systems approach

Nutrient pollution is considered a wicked problem because of its many significant economic, social, and environmental impacts that are caused by multiple pollutants originating from a variety of sources and pathways that exist across different temporal and spatial scales. Further adding to the difficulty in managing nutrient pollution is that it is a global, rural, and urban problem. A systems approach can improve nutrient management by incorporating technological, environmental, and societal considerations. This approach can consider valuation of monetized and nonmonetized co-benefits and the inherent consequences that make up a nutrient management program. In this introduction to a special collection of papers on nutrient pollution, we describe several systems frameworks that can be used to support nutrient management and evaluation of system performance as it relates to impacts, then highlight several attributes and barriers of nutrient management that point to the need for a systems framework, and conclude with thoughts on implementing systems approaches to nutrient management with effective community engagement and use of new technologies. This special collection presents results from a USEPA Science to Achieve Results (STAR) initiative to advance solutions to nutrient pollution through innovative and sustainable research and demonstration projects for nutrient management based on a systems approach. These studies evaluate several promising nutrient control technologies for stormwater or domestic wastewater, investigate the effects of agricultural conservation practices and stream restoration strategies on nutrient loads, and discuss several challenges and opportunities-social, policy, institutional, and financial considerations-that can accelerate adoption of reliable technologies to achieve system-level outcomes.

Impact factors and mechanisms of dissolved reactive phosphorus (DRP) losses from agricultural fields: A review and synthesis study in the Lake Erie basin

Dissolved Reactive Phosphorus (DRP) losses from agricultural fields promote algae growth in water bodies, and may increase the risk of Harmful Algal Blooms (HABs). Using existing data from the Lake Erie Basin, we applied multiple regression analysis to better understand the impacts of both site-specific conditions (e.g., soil types/properties) and management practices (e.g., Agricultural Conservation Practices [ACP]) on annual DRP losses in subsurface and surface runoff. Results showed that soil properties significantly impact DRP losses. Greater DRP losses were associated with increased soil pH and Soil Test Phosphorus (STP). By contrast, soil organic matter (SOM) was inversely correlated with DRP losses. Soil clay content was also inversely correlated with DRP subsurface losses, but had no impact on DRP surface losses. The ACPs evaluated had varied effectiveness on DRP loss reduction. Cropping systems involving soybean could reduce DRP subsurface losses, whereas winter cover crops could cause unintended DRP subsurface losses. Cropping systems involving soybean and cover crops, however, had no impact on DRP surface losses. In addition, no-till and conservation tillage also enhanced DRP losses compared to conventional tillage, particularly for soils with high SOM and/or high clay content. Precipitation amount and fertilizer application rate significantly increased DRP surface losses as expected. Fertilizer application rate, however, had no impact on DRP subsurface losses. The impact of precipitation amount on DRP subsurface losses depends on STP levels. Precipitation amount significantly increases DRP subsurface losses when STP is higher (>41 mg kg-1 in this analysis). The optimal STP level for crop growth is 30 to 50 mg kg-1. Results from this study help us to better understand DRP losses and the effectiveness of ACPs for controlling them. We suggest taking soil surveys and soil tests into consideration when designing and/or implementing ACPs to manage DRP losses. Furthermore, the method we used for this study could be applied to other agricultural regions to investigate impacts of site-specific conditions and management practices on water quality.

Spatiotemporal variations of nitrogen and phosphorus in a clay plain hydrological system in the Great Lakes Basin

Nutrient imbalance in groundwater and surface water resources can have severe implications on human and aquatic life, including contamination of drinking water sources and the degradation of ecosystems. A field-based watershed-scale study was completed to investigate nutrient dynamics and hydrologic processes in an agriculturally-dominant clay plain system within the Great Lakes Basin. Spatial and temporal variations of nitrogen and phosphorus were examined by sampling groundwater and surface water regularly over a period of one year (June 2017 to July 2018) for nutrients including nitrate, soluble reactive phosphorus, total dissolved phosphorus and total reactive phosphorus. Nitrate transport from surrounding agricultural land to surface water was intensified with an increase in precipitation events in spring and early winter and phosphorus transport to surface water was increased during freeze-thaw cycles in the winter. The results are pertinent to the improvement of current nutrient and water management policies in clay plain systems where nutrient imbalances in surface water are a concern.

Impact of P inputs on source-sink P dynamics of sediment along an agricultural ditch network

Phosphorus (P) loss from intensive dairy farms is a pressure on water quality in agricultural catchments. At farm scale, P sources can enter in-field drains and open ditches, resulting in transfer along ditch networks and delivery into nearby streams. Open ditches could be a potential location for P mitigation if the right location was identified, depending on P sources entering the ditch and the source-sink dynamics at the sediment-water interface. The objective of this study was to identify the right location along a ditch to mitigate P losses on an intensive dairy farm. High spatial resolution grab samples for water quality, along with sediment and bankside samples, were collected along an open ditch network to characterise the P dynamics within the ditch. Phosphorus inputs to the ditch adversely affected water quality, and a step change in P concentrations (increase in mean dissolved reactive phosphorus (DRP) from 0.054 to 0.228 mg L^-1) midway along the section of the ditch sampled, signalled the influence of a point source entering the ditch. Phosphorus inputs altered sediment P sorption properties as P accumulated along the length of the ditch. Accumulation of bankside and sediment labile extractable P, Mehlich 3 P (M3P) (from 13 to 97 mg kg^-1) resulted in a decrease in P binding energies (k) to < 1 L mg^-1 at downstream points and raised the equilibrium P concentrations (EPC₀) from 0.07 to 4.61 mg L^-1 along the ditch. The increase in EPC₀ was in line with increasing dissolved and total P in water, demonstrating the role of sediment downstream in this ditch as a secondary source of P to water. Implementation of intervention measures are needed to both mitigate P loss and remediate sediment to restore the sink properties. In-ditch measures need to account for a physicochemical lag time before improvements in water quality will be observed.

Development of phosphorus sorption capacity-based environmental indices for tile-drained systems

The persistent environmental relevance of phosphorus (P) and P sorption capacity (PSC) on P loss to surface waters has led to proposals for its inclusion in soil fertility and environmental management programs. As fertility and environmental management decisions are made on a routine basis, the use of laborious P sorption isotherms to quantify PSC is not feasible. Alternatively, pedotransfer functions (pedoTFs) estimate PSC from routinely assessed soil chemical properties. Our objective was to examine the possibility of developing a suitable pedoTF for estimating PSC and to evaluate subsequent PSC-based indices (P saturation ratio [PSR] and soil P storage capacity [SPSC]) using data from an in-field laboratory where tile drain effluent is monitored daily. Phosphorus sorption capacity was well predicted by a pedoTF derived from soil aluminum and organic matter (R² = .60). Segmented-line relationships between PSR and soluble P were observed in both desorption assays (R² = .69) and drainflows (R² = .66) with apparent PSR thresholds in close agreement at 0.21 and 0.24, respectively. Negative SPSC values exhibited linear relationships with increasing soluble P concentrations in both desorption assays and drainflows (R² = .52 and R²  = .53 respectively), whereas positive SPSC values were associated with low SP concentrations. Therefore, PSC-based indices determined using pedoTFs could estimate the potential for subsurface soluble P losses. Also, we determined that both index thresholds coincided with the critical soil-test P level for agronomic P sufficiency (22 mg kg^-1 Mehlich-3 P) suggesting that the agronomic threshold could serve as an environmental P threshold.

Ranking connectivity risk for phosphorus loss along agricultural drainage ditches

Agricultural drainage systems comprising both in-field pipe drains and surface ditches are typically installed to remove excess water from agricultural land. These drainage networks can provide connectivity between phosphorus (P) sources and surface waters thereby increasing the risk of P loss to rivers and streams. The objective of this study was to derive a farm-scale drainage ranking that categorises drainage ditches in terms of P loss risk based on connectivity and physic-chemical characteristics. Ten pilot farms were selected to characterise drainage networks through ground survey and, sediment and water sampling. Five drainage ditch categories were derived based on landscape setting and connectivity. Each category recorded soluble and reactive P concentrations above environmental water quality standards. To assess the risk of surface ditches as a connectivity vector between agricultural P and surface waters ditches were ranked in order of P loss risk by integrating landscape position and sediment P chemistry. Elevated sediment P with high equilibrium P concentration (EPCo) were associated with ditches connected to farm yards, and in sediment sampled at ditch outlets, suggesting P deposition over time indicative of a legacy P source. The greatest risk of P loss was attributed to ditches connecting farm yards to streams, and ditches that connected the drainage network to surface waters, or Outlets. These results rank connectivity risk for P loss along agricultural drainage ditches for farm level risk assessment to target P loss mitigation measures to the appropriate locations.

Impacts of Tile Drainage on Phosphorus Losses from Edge-of-Field Plots in the Lake Champlain Basin of New York

Quantifying the influence of tile drainage on phosphorus (P) transport risk is important where eutrophication is a concern. The objective of this study was to compare P exports from tile-drained (TD) and undrained (UD) edge-of-field plots in northern New York. Four plots (46 by 23 m) were established with tile drainage and surface runoff collection during 2012–2013. Grass sod was terminated in fall 2013 and corn (Zea mays L.) for silage was grown in 2014 and 2015. Runoff, total phosphorus (TP), soluble reactive phosphorus (SRP), and total suspended solids (TSS) exports were measured from April 2014 through June 2015. Mean total runoff was 396% greater for TD, however, surface runoff for TD was reduced by 84% compared to UD. There was no difference in mean cumulative TP export, while SRP and TSS exports were 55% and 158% greater for UD, respectively. A three day rain/snowmelt event resulted in 61% and 84% of cumulative SRP exports for TD and UD, respectively, with over 100% greater TP, SRP and TSS exports for UD. Results indicate that tile drainage substantially reduced surface runoff, TSS and SRP exports while having no impact on TP exports, suggesting tile drains may not increase the overall P export risk.

Edge-of-Field Technologies for Phosphorus Retention from Agricultural Drainage Discharge

Agriculture is often responsible for the eutrophication of surface waters due to the loss of phosphorus—a normally limiting nutrient in freshwater ecosystems. Tile-drained agricultural catchments tend to increase this problem by accelerating the transport of phosphorus through subsurface drains both in dissolved (reactive and organic phosphorus) and particulate (particle-bound phosphorus) forms. The reduction of excess phosphorus loads from agricultural catchments prior to reaching downstream surface waters is therefore necessary. Edge-of-field technologies have been investigated, developed and implemented in areas with excess phosphorus losses to receive and treat the drainage discharge, when measures at the farm-scale are not able to sufficiently reduce the loads. The implementation of these technologies shall base on the phosphorus dynamics of specific catchments (e.g., phosphorus load and dominant phosphorus form) in order to ensure that local retention goals are met. Widely accepted technologies include constructed wetlands, restored wetlands, vegetated buffer strips and filter materials. These have demonstrated a large variability in the retention of phosphorus, and results from the literature can help targeting specific catchment conditions with suitable technologies. This review provides a comprehensive analysis of the currently used edge-of-field technologies for phosphorus retention in tile-drained catchments, with great focus on performance, application and limitations.

Short-term impacts of biochar, tillage practices, and irrigation systems on nitrate and phosphorus concentrations in subsurface drainage water

Leaching of nitrogen (N) and phosphorus (P) from agricultural lands can cause serious environmental problems such as eutrophication. The objective of this study was to investigate the impacts of biochar application, tillage practices, and irrigation systems on nitrate and dissolved phosphorus (DP) concentrations in subsurface drainage water and grain yield of winter wheat using a strip-split plot design with 3 replications. Irrigation at three different levels (flood (I_fl), furrow (I_fu), and sprinkler (I_s) systems) considered as main factor, tillage at two levels (reduced tillage (T_r) and conventional systems (T_c)) as subplot factor, and bagasse biochar at two levels (without biochar (B₀) and 20 ton ha^-1 biochar (B₁)) as sub-subplot factor. Polyvinyl chloride (PVC) standpipes were used in each sub-subplot to collect leachate water at 100-cm depth. The results indicated that irrigation had significant effects on yield, collected water volume (CWV), nitrate, and DP concentrations (P < 0.01). Interaction of tillage and irrigation was significant for grain yield (P < 0.05). Biochar application only caused a significant decrease in nitrate concentration under sprinkler irrigation (P < 0.05), while no significant impact was observed under flood and furrow irrigation systems. Under sprinkler irrigation, the total nitrate collected in the PVC standpipes decreased by 37.51 and 34.29% compared with flood and furrow irrigations, respectively. Biochar application reduced the total nitrate collected by 16.84%, while difference among tillage treatments was negligible (4.51%). The total DP collected under sprinkler irrigation was lower in comparison with flood and furrow irrigations by 42.24 and 38.76%, respectively. Biochar application reduced the total DP collected by 10.84%, while reduced tillage increased the total DP collected by 8.90% compared with the conventional tillage.

Synergistic effects of anionic polyacrylamide and gypsum to control phosphorus losses from biogas slurry applied soils

Excessive application of phosphorus (P) fertilizer to the agriculture soil can lead to the P loss to the aquatic ecosystem. Three soils from different land use (tea, paddy, and vegetable soils) were treated with biogas slurry (BS) at rate of 79.8 kg ha^-1, anionic polyacrylamide (PAM) at rate of 12.5, 25 and 50 kg ha^-1, Gypsum (Gy) at rate of 5t ha^-1, and the control (Ck) without treatment, to control P losses in a laboratory batch experiment. The soil solution was analyzed for total P (TP), dissolved reactive P (DRP) and total dissolved P (TDP). Particulate P (PP) and dissolved unreactive P (DURP) were calculated by subtracting DRP from TP and TDP, respectively. The DRP was within a range of 68-98% of TP whereas PP varied in a range of 2-32% of TP and DURP was within a range of 0.1-19% of TP, depending on soil types and different treatments. The addition of PAM at a high rate and Gy together on the soil after BS application reduced the release of TP by 54, 46, and 51% and that of DRP by 61, 49, and 53% for tea, paddy and vegetable soils, respectively. However, the application of BS with PAM alone promoted the release of DRP by 77, 86, 70% for tea, paddy, and vegetable soils, respectively. This study suggests that the synergic application of PAM and Gy on soils after BS addition can be a good strategy to reduce P losses.

Direct Exports of Phosphorus from Fertilizers Applied to Grazed Pastures

Since its discovery in 1669, phosphorus (P) in the form of fertilizer has become an essential input for many agroecosystems. By introducing a concentrated P source, fertilizers increase short-term P export potential soon after their application and longer-term export potential by increasing soil fertility (legacy P). The 4R concept was developed to help mitigate P exports from the fertilizers that sustain agricultural productivity. This review investigates the factors affecting P exports soon after the application of mineral fertilizers to pasture-based grazing systems and studies quantifying its potential impact in different systems, with an emphasis on Australasia. Initially, P fertilizers and reactions that might affect their short-term P export potential are reviewed, along with P transport pathways, the forms of P exported from grazing systems, factors affecting P mobilization into water, and studies demonstrating the possible short-term effects of fertilizer application on P exports. Using that foundation, we review studies quantifying the short-term impact of fertilizer application in different regions; they show that under poor management, recently applied fertilizer can contribute a considerable proportion (30-80%) of total farm P exports in drainage, but when fertilizer is well-managed, that figure is expected to be <10%. We then use three model systems of varying hydrology that are common to Australasia to demonstrate the principles for selecting fertilizers that are likely to minimize P exports soon after their application.

The Latitudes, Attitudes, and Platitudes of Watershed Phosphorus Management in North America

Phosphorus (P) plays a crucial role in agriculture as a primary fertilizer nutrient-and as a cause of the eutrophication of surface waters. Despite decades of efforts to keep P on agricultural fields and reduce losses to waterways, frequent algal blooms persist, triggering not only ecological disruption but also economic, social, and political consequences. We investigate historical and persistent factors affecting agricultural P mitigation in a transect of major watersheds across North America: Lake Winnipeg, Lake Erie, the Chesapeake Bay, and Lake Okeechobee/Everglades. These water bodies span 26 degrees of latitude, from the cold climate of central Canada to the subtropics of the southeastern United States. These water bodies and their associated watersheds have tracked trajectories of P mitigation that manifest remarkable similarities, and all have faced challenges in the application of science to agricultural management that continue to this day. An evolution of knowledge and experience in watershed P mitigation calls into question uniform solutions as well as efforts to transfer strategies from other arenas. As a result, there is a need to admit to shortcomings of past approaches, plotting a future for watershed P mitigation that accepts the sometimes two-sided nature of Hennig Brandt's "Devil's Element."

Increasing the Effectiveness and Adoption of Agricultural Phosphorus Management Strategies to Minimize Water Quality Impairment

Phosphorus (P) is essential for optimum agricultural production, but it also causes water quality degradation when lost through erosion (sediment-attached P), runoff (soluble reactive P; SRP), or leaching (sediment-attached P or SRP). Implementation of conservation practices (CP) affects P at the source (avoiding), during transport (controlling), or at the water resource edge (trapping). Trade-offs often occur with CP implementation. For instance, multiple researchers have shown that conservation tillage reduces total P by over 50%, while increasing SRP by upward of 40%. Conservation tillage may increase water quality degradation as SRP is more bioavailable than is particulate P. Conservation practices must be implemented as a system of practices to increase redundancy and to address all loss pathways, such as P management with conservation tillage and a riparian buffer. Further, planning and adoption must be at a watershed scale to ensure practices are placed in critical source areas, thereby providing the most treatment for the least price. Farmers must be involved in watershed planning, which should include financial backstopping and educational outreach. It is imperative that CPs be used more effectively to reduce and retard off-site P losses. New and innovative CPs are needed to improve control of P leaching, address legacy stores of soil test P, and mitigate increased P losses expected with climate change. Without immediate changes to CP implementation, P losses will increase due to climate change, with a concomitant degradation of water quality. These changes must be made at a watershed scale and in an intentional and transparent manner.

Phosphorus and the Chesapeake Bay: Lingering Issues and Emerging Concerns for Agriculture

Hennig Brandt's discovery of phosphorus (P) occurred during the early European colonization of the Chesapeake Bay region. Today, P, an essential nutrient on land and water alike, is one of the principal threats to the health of the bay. Despite widespread implementation of best management practices across the Chesapeake Bay watershed following the implementation in 2010 of a total maximum daily load (TMDL) to improve the health of the bay, P load reductions across the bay's 166,000-km watershed have been uneven, and dissolved P loads have increased in a number of the bay's tributaries. As the midpoint of the 15-yr TMDL process has now passed, some of the more stubborn sources of P must now be tackled. For nonpoint agricultural sources, strategies that not only address particulate P but also mitigate dissolved P losses are essential. Lingering concerns include legacy P stored in soils and reservoir sediments, mitigation of P in artificial drainage and stormwater from hotspots and converted farmland, manure management and animal heavy use areas, and critical source areas of P in agricultural landscapes. While opportunities exist to curtail transport of all forms of P, greater attention is required toward adapting P management to new hydrologic regimes and transport pathways imposed by climate change.

Evaluating Hydrologic Response in Tile-Drained Landscapes: Implications for Phosphorus Transport

Phosphorus (P) loss in agricultural discharge has typically been associated with surface runoff; however, tile drains have been identified as a key P pathway due to preferential transport. Identifying when and where these pathways are active may establish high-risk periods and regions that are vulnerable for P loss. A synthesis of high-frequency, runoff data from eight cropped fields across the Great Lakes region of North America over a 3-yr period showed that both surface and tile flow occurred year-round, although tile flow occurred more frequently. The relative timing of surface and tile flow activation was classified into four response types to infer runoff-generation processes. Response types were found to vary with season and soil texture. In most events across all sites, tile responses preceded surface flow, whereas the occurrence of surface flow prior to tile flow was uncommon. The simultaneous activation of pathways, indicating rapid connectivity through the vadose zone, was seldom observed at the loam sites but occurred at clay sites during spring and summer. Surface flow at the loam sites was often generated as saturation-excess, a phenomenon rarely observed on the clay sites. Contrary to expectations, significant differences in P loads in tiles were not apparent under the different response types. This may be due to the frequency of the water quality sampling or may indicate that factors other than surface-tile hydrologic connectivity drive tile P concentrations. This work provides new insight into spatial and temporal differences in runoff mechanisms in tile-drained landscapes.

Study on Characteristics of Nitrogen and Phosphorus Loss under an Improved Subsurface Drainage

Agricultural drainage plays an effective role in preventing waterlogging and salinity disasters and also is the main transport pathway for agricultural non-point source pollutants into rivers and lakes. Hence, the water quality of agricultural drainage should be a point of focus. In this paper, nitrogen and phosphorus loss under improved subsurface drainage with different filter materials (gravel, layered sand-gravel, mixed sand-gravel, straw) were studied by a three-year field experiment (2016–2018) compared with the conventional subsurface drainage. The pH values, total nitrogen, ammonia nitrogen, nitrate nitrogen, total phosphorus and soluble reactive phosphate were considered. The results showed that the nitrogen and phosphorus concentrations of drain outflow under improved subsurface drainage with gravel filter were larger than that with layered sand-gravel filter and mixed sand-gravel filter. The improved subsurface drainages with layered sand-gravel filter and mixed sand-gravel filter had an effect on reducing the ammonia nitrogen, total phosphorus and soluble reactive phosphate concentrations of the outflow. Meanwhile, the characteristics of nitrogen and phosphorus loss under the improved subsurface drainage with straw filter were different from that with layered sand-gravel filter and mixed sand-gravel filter. For the improved subsurface drainage with layered sand-gravel filter outflow, the ammonia nitrogen, total phosphorus, and soluble reactive phosphate concentrations were about 13%–78%, 38%–63%, 40%–68% less, and total nitrogen, nitrate nitrogen concentrations were 24%–80%,18%–96% more than that under conventional subsurface drainage. Meanwhile, for the improved subsurface drainage with straw filter outflow, compared with conventional subsurface drainage outflow, the percentage changes of the total nitrogen, nitrate nitrogen, ammonia nitrogen, total phosphorus and the soluble reactive phosphate concentrations were about −76%–62%, −77%–78%, −152%–−274%, −103%–−400% and −221%–−291%, respectively. Additionally, in the outflow of all subsurface drainage patterns, there were much higher total nitrogen and nitrate nitrogen concentrations which should be focused on and the agricultural water management should be adopted.

Long-term Variation in Agricultural Edge-of-Field Phosphorus Transport during Snowmelt, Rain, and Mixed Runoff Events

Runoff of nutrients and erosion of soil from agricultural lands affect soil fertility and are important nonpoint contributors of P and N to surface and ground waters, yet studies of edge-of-field nutrient transport from snowmelt or rainfall runoff on frozen ground are limited. The objective of this study was to quantify the temporal and spatial variation in edge-of-field snowmelt, rain, and mixed (rain on snow) runoff events for sediment and P loadings in five agricultural subwatersheds over a 12-yr period. Edge-of-field runoff events from five subwatersheds at Pioneer Farm near Platteville, WI, ranging in size from approximately 4 to 30 ha were sampled using automated samplers from 2002 through 2014 to determine sediment and P yields (mass loads). Mean dissolved reactive P (DRP) runoff concentrations for each event type (rain = 1.24 mg L, snow = 1.90 mg L, mix = 2.23 mg L) were above total P (TP) water quality guidelines for surface waters. The percentages of TP that was DRP for snow, mixed, and rain events were 74, 84, and 39%, respectively. Although variation in total annual P yield in edge-of-field runoff was noted between years and among sites within a given year, when aggregated over the study period, the subwatersheds showed similar transport characteristics with respect to DRP and TP yield. This study highlights the importance of examining long-term datasets in quantifying annual yields and understanding the timing of DRP and TP transport for developing best management practices and improving model accuracy in cold weather agricultural systems.

Winter Phosphorus Release from Cover Crops and Linkages with Runoff Chemistry

Cover crops (CC) have both agronomic and environmental benefits but also have the potential to increase losses of dissolved reactive P after freeze-thaw cycles (FTC). This field study, conducted over one nongrowing season (NGS) in Ontario, Canada, characterized water-extractable P (WEP) content in different CC species and compared observed changes in plant WEP content with changes in P content in soil, surface runoff, and shallow groundwater (5-25 cm). Five plots (0.4 ha) of cereal rye ( L.), oilseed radish ( L. var. Metzg Stokes), oat ( L.), and hairy vetch ( Roth) were established after winter wheat ( L.) harvest. Throughout the NGS (October-April), CC shoot tissues and surface soil were routinely sampled for WEP analyses, and groundwater and runoff water samples were collected after rain and snowmelt. Responses to FTC varied among CC species, with P released from frost-intolerant species but not frost-tolerant species. Although CC released P, the top 5 cm of soil contained greater WEP than plants at all times, and the changing WEP content in CC over the NGS was not reflected in soil or water P concentrations. These results suggest that the degree of frost exposure should be considered in the selection of CC species in cold regions; however, in temperate regions with snow cover that insulates the soil surface from heavy frost, P release from vegetation may not lead to increased P loss in runoff.

Impacts of Cover Crops and Crop Residues on Phosphorus Losses in Cold Climates: A Review

The use of cover crops and crop residues is a common strategy to mitigate sediment and nutrient losses from land to water. In cold climates, elevated dissolved P losses can occur associated with freeze-thaw of plant materials. Here, we review the impacts of cover crops and crop residues on dissolved P and total P loss in cold climates across ∼41 studies, exploring linkages between water-extractable P (WEP) in plant materials and P loss in surface runoff and subsurface drainage. Water-extractable P concentrations are influenced by plant type and freezing regimes. For example, WEP was greater in brassica cover crops than in non-brassicas, and increased with repeated freeze-thaw cycles. However, total P losses in surface runoff and subsurface drainage from cropped fields under cold climates were much lower than plant WEP, owing to retention of 45 to >99% of released P by soil. In cold climatic regions, cover crops and crop residues generally prevented soil erosion and loss of particle-bound P during nongrowing seasons in erodible landscapes but tended to elevate dissolved P loss in nonerodible soils. Their impact on total P loss was inconsistent across studies and complicated by soil, climate, and management factors. More research is needed to understand interactions between these factors and plant type that influence P loss, and to improve the assessment of crop contributions to P loss in field settings in cold climates. Further, tradeoffs between P loss and the control of sediment loss and N leaching by plants should be acknowledged.

Phosphorus Release from Unamended and Gypsum- or Biochar-Amended Soils under Simulated Snowmelt and Summer Flooding Conditions

Prolonged flooding changes the oxidation-reduction status of soils, often enhancing P release to overlying floodwater. We studied P release from unamended, gypsum-amended, and biochar-amended soils under simulated snowmelt flooding (previously frozen, cold flooding at +4°C) and summer flooding (unfrozen, warm flooding at +22°C) using two soils, Fyala clay (FYL-Cl) and Neuenberg sandy loam (NBG-SL), from Manitoba, Canada. Amended and unamended soils were packed into vessels and flooded under cold and warm temperatures in the laboratory. Pore water and floodwater samples were taken weekly for 6 wk after flooding (WAF) and thereafter biweekly for 10 WAF and analyzed for dissolved reactive P (DRP), pH, and cation concentrations. The NBG-SL showed a significantly higher DRP concentration in pore water and floodwater despite its low Olsen P content. Redox potential (Eh) decreased slowly under cold versus warm flooding; hence, redox-induced P release was substantially lower under cold flooding. Gypsum amendment significantly decreased the floodwater DRP concentrations in NBG-SL by 38 and 35% under cold and warm flooding, respectively, but had no significant effect in FYL-Cl, which had low DRP concentrations (<1.2 mg L) throughout the flooding period. Biochar amendment significantly increased floodwater DRP concentrations by 27 to 68% in FYL-Cl under cold and warm flooding, respectively, but had no significant effect in NBG-SL. The results indicate substantially less P release under cold than under warm flooding. Gypsum was effective in reducing floodwater DRP concentrations only at high DRP concentrations; thus, the effectiveness was greater under warm than under cold flooding conditions.

Modeling and Mitigating Phosphorus Losses from a Tile-Drained and Manured Field Using RZWQM2-P

Prediction of P losses from manured agricultural fields through surface runoff and tile drainage is necessary to mitigate widespread eutrophication in water bodies. However, present water quality models are weak in predicting P losses, particularly in tile-drained and manure-applied cropland. We developed a field-scale P management model, the Root Zone Water Quality Model version 2-Phosphorus (RZWQM2-P), whose accuracy in simulating P losses from manure applied agricultural field is yet to be tested. The objectives of this study were (i) to assess the accuracy of this new model in simulating dissolved reactive phosphorus (DRP) and particulate phosphorus (PP) losses in surface runoff and tile drainage from a manure amended field, and (ii) to identify best management practices to mitigate manure P losses including water table control, manure application timing, and spreading methods by the use of model simulation. The model was evaluated against data collected from a liquid cattle manure applied field with maize ( L.)-soybean [ (L.) Merr.] rotation in Ontario, Canada. The results revealed that the RZWQM2-P model satisfactorily simulated DRP and PP losses through both surface runoff and tile drainage (Nash-Sutcliffe efficiency > 0.50, percentage bias within ±25%, and index of agreement > 0.75). Compared with conventional management practices, manure injection reduced the P losses by 18%, whereas controlled drainage and winter manure application increased P losses by 13 and 23%, respectively. The RZWQM2-P is a promising tool for P management in manured and subsurface drained agricultural field. The injection of manure rather than controlled drainage is an effective management practice to mitigate P losses from a subsurface-drained field.

Springs drive downstream nitrate export from artificially-drained agricultural headwater catchments

Excessive nutrient loading from small agricultural headwaters can substantially degrade downstream water quality and ecological conditions. But, our understanding of the scales and locations to implement nutrient attenuation tools within these catchments is poor. To help inform farm- and catchment-scale management, we quantified nitrate export in nine one-kilometre-long lowland agricultural headwaters fed by tile and open tributary drains in a region with high groundwater nitrate (<1 to >15 mg L^-1 NO₃-N) over four years. Across-catchment differences in upstream spring water nitrate concentrations predicted differences in annual nitrate loads at catchment outlets (range <1-72 megagrams NO₃-N 365 d^-1), and nitrate loads were higher in wet seasons and wet years, reflecting strong groundwater influences. Partitioning the sources of variability in catchment nitrate fluxes revealed that ~60% of variation was accounted for by a combination of fluxes from up-stream springs and contributions from tile and open tributary drains (46% and 15%, respectively), with ~40% of unexplained residual variation likely due to groundwater upwellings. Although tile and open tributary drains contributed comparatively less to catchment loads (tile drains: <0.01 and up to 50 kg NO₃-N d^-1; open drains: <5 kg and up to 100 kg NO₃-N d^-1), mitigation targeted at these localised, farm-scale sources will contribute to decreasing downstream nitrate fluxes. However, high nitrate loads from groundwater mean current NO₃-N waterway management and rehabilitation practices targeting waterway stock exclusion by fencing alone will be insufficient to reduce annual NO₃-N export. Moreover, managing catchment nutrient fluxes will need to acknowledge contributions from groundwater as well as farm-scale losses from land. Overall, our results highlight how nutrient fluxes in spring-fed waterways can be highly dynamic, dominated more by groundwater than local run-off, and point to the scales and locations where nitrate attenuation tools should be implemented.

Phosphorus Leaching from Soil Profiles in Agricultural and Forest Lands Measured by a Cascade Extraction Method

The risk of P leaching from topsoil based on the change-point estimated via a split-line model between Olsen P and leachable P extracted by 0.01 M CaCl has been reported. However, little information is available for the assessment of P leaching from soil profiles. In this study, samples were collected at three depth profiles (0-20 cm, topsoil; 20-40 cm, subsoil; 40-60 cm, third-layer soil) at each of 74 sites under agriculture and forest in an agroforestry area. A cascade extraction method was proposed to determine the leachable P in the subsoil, extracted by the topsoil extraction solution; a similar extracted process was followed in the third-layer soil, and in the topsoil, it was still extracted by 0.01 M CaCl. A positive linear correlation was found between the subsoil leachable P extracted by the topsoil extraction solution and the accumulated P obtained from the subsoil leached by topsoil leachates, and so on. Therefore, the cascade extraction method for determining leachable P from topsoils and underlying soils could be valuable for predicting the potential of P leaching from soil profiles. Approximately 81, 73, and 73% of the agricultural sampling sites were at or above the change-points for the soil depths of 0 to 20, 20 to 40, and 40 to 60 cm (30.4, 32.9, and 18.2 mg kg respectively); these values were higher than those for the forest site, implying a high risk of P leaching from agricultural soil profiles in the study area.

Emulating natural wetlands oxygen conditions for the removal of N and P in agricultural wastewaters

The aim of this research was to evaluate a constructed wetland system (CW) operated under aerobic-anoxic-aerobic conditions to remove C, N and P from water with high concentrations of the last two nutrients. A series of three CW were operated continuously for 190 days. An aerobic vertical CW was used in the first and third stages and an anoxic horizontal CW was used in the second stage. The total nitrogen (TN) removal efficiency was 70 ± 1.5%. Similar removal efficiency behavior was observed in others nitrogen compounds, where a removal of 85 ± 1.5% for NO₃^--N and 97 ± 2.2% for NH₃⁺N were achieved. The combination of different oxygen conditions enhanced oxidation of nitrates and the assimilation of ammonium by vegetation. On the other hand, 54 ± 6.5% total phosphorus (TP) was removed in the entire system, which is higher than the reported in several investigations, including mechanized and controlled systems such as activated sludge. The phosphorous removal efficiency was attributed to the adequate design and configuration of CW, which facilitated dissolved oxygen (DO) conditions required for phosphorus capture. Despite in this investigation the CW was not designed for an optimal removal of organic matter the removal efficiency of this parameter was 64 ± 7.5%. The successful results suggest that the combination of aerobic-anoxic-aerobic stages is a technically suitable option for the treatment of agricultural wastewater with high content of N and P.

Hydrograph separation of subsurface tile discharge

Baseflow is an important component of streamflow and watershed hydrologic budgets, yet quantifying the baseflow fraction of tile drainage has rarely been reported. In this study, we used two common hydrograph separation methods (local minimum method, recursive digital filter) to separate the discharge hydrographs from three drainage district tiles located in Iowa. Based on data collected from 2009 to 2013, annual baseflow ranged from 116 to 162 mm and comprised approximately 60% of the annual discharge. Baseflow was greatest during June (average of 34% of annual baseflow) and the March through August period produced 86% of the total annual baseflow. We found that the two methods of hydrograph separation produced similar results but the digital filter method was less erratic in estimating baseflow fraction. Study results can be used to better quantify hydrologic pathways in tiled landscapes and improve the design, implementation, and evaluation of nutrient reduction strategies.

Harnessing Soil Microbes to Improve Plant Phosphate Efficiency in Cropping Systems

Phosphorus is an essential macronutrient required for plant growth and development. It is central to many biological processes, including nucleic acid synthesis, respiration, and enzymatic activity. However, the strong adsorption of phosphorus by minerals in the soil decreases its availability to plants, thus reducing the productivity of agricultural and forestry ecosystems. This has resulted in a complete dependence on non-renewable chemical fertilizers that are environmentally damaging. Alternative strategies must be identified and implemented to help crops acquire phosphorus more sustainably. In this review, we highlight recent advances in our understanding and utilization of soil microbes to both solubilize inorganic phosphate from insoluble forms and allocate it directly to crop plants. Specifically, we focus on arbuscular mycorrhizal fungi, ectomycorrhizal fungi, and phosphate-solubilizing bacteria. Each of these play a major role in natural and agroecosystems, and their use as bioinoculants is an increasing trend in agricultural practices.

Development of PLEAD: A Database Containing Event-based Runoff Phosphorus Loadings from Agricultural Fields

Computer models are commonly used for predicting risks of runoff P loss from agricultural fields by enabling simulation of various management practices and climatic scenarios. For P loss models to be useful tools, however, they must accurately predict P loss for a wide range of climatic, physiographic, and land management conditions. A complicating factor in developing and evaluating P loss models is the relative scarcity of available measured field data that adequately capture P losses before and after implementing management practices in a variety of physiographic settings. Here, we describe the development of the P Loss in runoff Events from Agricultural fields Database (PLEAD)-a compilation of event-based, field-scale dissolved and/or total P loss runoff loadings from agricultural fields collected at various research sites located in the US Heartland and southern United States. The database also includes runoff and erosion rates; soil-test P; tillage practices; planting and harvesting rates and practices; fertilizer application rate, method, and timing; manure application rate, method, and timing; and livestock grazing density and timing. In total, >1800 individual runoff events-ranging in duration from 0.4 to 97 h-have been included in the database. Event runoff P losses ranged from <0.05 to 1.3 and 3.0 kg P ha for dissolved and total P, respectively. The data contained in this database have been used in multiple research studies to address important modeling questions relevant to P management planning. We provide these data to encourage additional studies by other researchers. The PLEAD database is available at .

Point sources and agricultural practices control spatial-temporal patterns of orthophosphate in tributaries to Chesapeake Bay

Orthophosphate (PO₄) is the most bioavailable form of phosphorus (P). Excess PO₄ may cause harmful algal blooms in aquatic ecosystems. A major restoration effort is underway for Chesapeake Bay (CB) to reduce P, nitrogen, and sediment loading to CB. Although PO₄ cycling and delivery to streams has been characterized in small-scale studies, regional drivers of PO₄ patterns remain poorly understood because most water quality trend assessment focus on total P. Moreover, these trend assessments are usually at an annual timestep. To address this research gap, we analyzed PO₄ patterns over a 9-year period at 53 monitoring stations across the CB watershed to: 1) characterize the role of PO₄ in total P fluxes and trends; 2) describe spatial and temporal patterns of PO₄ concentrations across seasons and streamflow; and 3) explore factors explaining these patterns. Agricultural watersheds exported the most total P compared with watersheds under different land uses (e.g., urban or forest), with PO₄ comprising up to 50% of those exports. Although PO₄ exports are declining at many sites, some agricultural regions are experiencing increasing trends at a rate sufficient to drive total P trends. Regression modeling results suggest that point source load reductions are likely responsible for decreasing PO₄ concentrations observed at many sites. Watersheds with more Conservation Reserve Program enrollment had lower summer PO₄ concentrations, highlighting the effectiveness of this practice. Manure inputs strongly predicted PO₄ concentrations at high flows across all seasons. Both manure applications and conservation tillage were correlated with changes in PO₄ concentrations at high flow, suggesting these activities could contribute to increasing PO₄ concentrations. This study highlights the effectiveness of point source control for reducing PO₄ exports and underscores the need for management strategies to target sources, practices, and landscape factors determining PO₄ loss from soils where manure inputs remain high.

The implications of weather, nutrient prices, and other factors on nutrient concentrations in agricultural watersheds

This paper examines how nutrient prices, weather, and other factors influenced P outputs in agricultural watersheds using a detailed daily dataset of water quality observations over a 40-year period. Because policies have focused differentially on soluble P through federal permitting programs for point sources and sediments through federal subsidies for conservation, we examine sediment, particulate P and soluble P separately. A novel element of this study is the inclusion of farm fertilizer and output (i.e., corn) prices, which affect agricultural sources of P in these watersheds. We do not find that sediment concentrations are influenced by P prices, but sediment has trended downward, and is seasonally lower in all months except February and March in the Maumee. In contrast, we find that soluble P concentrations are heavily influenced by P prices. They trended downward through 1995, but upwards since. While concerns about fall and winter P application have emerged, we do not find evidence that the distribution of soluble P concentrations shifted towards winter over time. Weather accounts for about 50% of the higher soluble P loadings in 1996-2011, but higher P prices in 2005-2011 lowered P concentrations relative to what they would have been. Other factors account for the remaining 50% of the increase in soluble P concentrations in 1996-2011.

Gypsum Amendment Reduces Redox-Induced Phosphorous Release from Freshly Manured, Flooded Soils to Floodwater

The effectiveness of gypsum in reducing runoff P losses from soils and the mechanisms responsible are well documented; however, gypsum amendment effects in reducing redox-induced P losses from flooded soils are less researched and documented. We examined the effect of gypsum amendment on P release from freshly manured soils to pore water and floodwater with continuous flooding for 56 d in the laboratory. Three soils (Pembina, Denham, and Dencross series) collected from Manitoba, Canada, were preincubated with liquid swine manure. Each preincubated manured soil was packed into vessels with or without recycled wallboard gypsum in triplicates and flooded for 56 d, during which pore water and floodwater were sampled weekly and analyzed for pH and dissolved reactive P (DRP), Ca, Mg, Fe, and Mn concentrations. Change in soil redox potential (Eh) with flooding was also monitored. Wallboard gypsum amendment significantly decreased the pore water and surface floodwater DRP concentrations in all three soils for most days after flooding (DAF). The Dencross soil, which had Olsen P about fivefold greater than the other soils, showed the greatest magnitude decrease in DRP concentration with gypsum amendment, by 1.27 mg L on 49 DAF and 0.99 mg L on 21 DAF for pore water and floodwater, respectively. Gypsum amendment (i) delayed the Eh reduction with flooding beyond +200 mV, (ii) decreased pore water pH, and (iii) increased concentrations of Ca, Mg, and Mn in pore water favoring precipitation of P, all of which may have directly or indirectly reduced the P release from flooded soils to overlying floodwater.

Surface and Subsurface Phosphorus Discharge from a Clay Soil in a Nine-Year Study Comparing No-Till and Plowing

No-till as a water protection measure is highly efficient in controlling erosion and particulate P (PP) loss but tends to increase dissolved reactive P (DRP) concentrations in runoff water. In a 9-yr field study on a clay soil in Southwest Finland, the effects of no-till and autumn plowing on surface runoff and subsurface drainage water quality were compared. The site had a 2% slope and was under spring cereal cropping, with approximately replacement fertilizer P rates. Vertical stratification of soil-test P that had developed during a preceding 6-yr grass ley was undone by plowing but continued to develop under no-till. During the 9-yr study period, no-till soil had 27% lower cumulative total P losses than plowed soil (10.0 vs. 13.7 kg total P ha). Concentrations and losses of PP were clearly lower under no-till than under plowing (5.6 vs. 12.3 kg PP ha), but DRP loss showed the opposite trend (4.3 vs. 1.4 kg DRP ha). There was an increasing trend in subsurface drainflow DRP concentration under no-till, possibly because of development of a conductive pore structure from soil surface to drain depth. The potential benefit of no-till in water protection depends on how much of the PP transported to water is transformed into a bioavailable form and used by aquatic organisms. The beneficial effect of no-till in controlling P-induced eutrophication at the study site would only be realized if the bioavailable share of PP exceeds 43%. Otherwise, no-till would not be an efficient eutrophication control measure at this site.