Phosphorus transport in agricultural subsurface drainage: a review

Lanthanum-loaded geopolymer for phosphate removal from agricultural runoff

The removal of phosphate from agricultural runoff is of great importance to mitigate an overabundance of nutrients discharging into receiving water bodies, which are susceptible to eutrophication. In this study, a La-loaded geopolymer was produced by adding metakaolin within an alkaline medium, consisting of sodium silicate (Na2SiO3) and sodium hydroxide (NaOH). The ratio of Na2SiO3:NaOH within the geopolymer slurry was manipulated to evaluate its effect on phosphate adsorption capacity. The 1.54 ratio yielded the highest adsorption capacity of 33.65 mg g-1. However, due to structural strength, safety, and economic considerations, the 2.0 ratio was used for the isotherm and kinetic adsorption testing. The La-loaded geopolymer exhibited higher phosphate removal in batch adsorption experiments at an acidic level (pH 4) and elevated temperature (40 °C). Moreover, ionic strength (3.5-20 mM) had a negligible impact on phosphate removal, indicating inner-sphere complexation as the main mechanism of adsorption. Additionally, bicarbonate and humic acid increased phosphate removal, whereas sulfate slightly decreased adsorption capacity. The La-loaded geopolymer was further evaluated using a synthetic agricultural solution, which yielded a Type III adsorption isotherm, demonstrating unrestricted multilayer phosphate adsorption. Results from this study verified the La-loaded geopolymer is a promising adsorbent for phosphate removal from agricultural runoff and underlined the importance of matrix effects on treatment performance.

A regional examination of the footprint of agriculture and urban cover on stream water quality

Freshwater systems in cold regions, including the Laurentian Great Lakes, are threatened by both eutrophication and salinization, due to excess nitrogen (N), phosphorus (P) and chloride (Cl-) delivered in agricultural and urban runoff. However, identifying the relative contribution of urban vs. agricultural development to water quality impairment is challenging in watersheds with mixed land cover, which typify most developed regions. In this study, a self-organizing map (SOM) analysis was used to evaluate the contributions of various forms of land cover to water quality impairment in southern Ontario, a population-dense, yet highly agricultural region in the Laurentian Great Lakes basin where urban expansion and agricultural intensification have been associated with continued water quality impairment. Watersheds were classified into eight spatial clusters, representing four categories of agriculture, one urban, one natural, and two mixed land use clusters. All four agricultural clusters had high nitrate-N concentrations, but levels were especially high in watersheds with extensive corn and soybean cultivation, where exceedances of the 3 mg L-1 water quality objective dramatically increased above a threshold of ‍∼30 % watershed row crop cover. Maximum P concentrations also occurred in the most heavily tile-drained cash crop watersheds, but associations between P and land use were not as clear as for N. The most urbanized watersheds had the highest Cl- concentrations and expansions in urban area were mostly at the expense of surrounding agricultural land cover, which may drive intensification of remaining agricultural lands. Expansions in tile-drained corn and soybean area, often at the expense of mixed, lower intensity agriculture are not unique to this area and suggest that river nitrate-N levels will continue to increase in the future. The SOM approach provides a powerful means of simplifying heterogeneous land cover characteristics that can be associated with water quality patterns and identify problem areas to target management.

Perspectives on total phosphorus response in rivers: Examining the influence of rainfall extremes and post-dry rainfall

Eutrophication is a significant environmental problem caused by nutrient loads from both point and non-point sources. Weather variables, particularly precipitation, affect the concentration of nutrients in water bodies, particularly those from non-point sources, in two contrasting ways. Heavy precipitation causes surface runoff which transports pollutants to rivers and increases nutrient concentration. Conversely, increased river flow can dilute the concentration, lowering it. This study investigates the impact of extreme precipitation, prolonged precipitation, and precipitation after a dry period on the total phosphorus concentration in the Moehne and Erft rivers in Germany, given the projected increase in frequency of extreme precipitation events and long drought periods due to climate change. The study comprises two parts: selecting extreme weather days from 2001 to 2021 and comparing observed Total Phosphorus concentrations with estimated concentrations derived from Generalized Additive Models and linear regression based on the discharge-concentration relationship. Changes in river TP concentration in response to continuous precipitation and precipitation after a dry period were also studied. Our results showed that during wet extreme and post-dry period rainfall events, TP concentration consistently surpassed expected values, underscoring the profound influence of intense rainfall on nutrient mobilization. However, we observed the impact of continuous rainfall to be non-unidirectional. Our work is distinguished by three key innovations: 1) addressing limitations in studying the effects of extreme weather on water quality due to limited temporal resolution, 2) incorporating both linear and non-linear modeling approaches for discharge-concentration relationships, and 3) performing a comprehensive analysis of temporal and spatial patterns of Total Phosphorus concentrations in response to varying rainfall patterns.

Mapping agricultural tile drainage in the US Midwest using explainable random forest machine learning and satellite imagery

There has been an increase in tile drained area across the US Midwest and other regions worldwide due to agricultural expansion, intensification, and climate variability. Despite this growth, spatially explicit tile drainage maps remain scarce, which limits the accuracy of hydrologic modeling and implementation of nutrient reduction strategies. Here, we developed a machine-learning model to provide a Spatially Explicit Estimate of Tile Drainage (SEETileDrain) across the US Midwest in 2017 at a 30-m resolution. This model used 31 satellite-derived and environmental features after removing less important and highly correlated features. It was trained with 60,938 tile and non-tile ground truth points within the Google Earth Engine cloud-computing platform. We also used multiple feature importance metrics and Accumulated Local Effects to interpret the machine learning model. The results show that our model achieved good accuracy, with 96 % of points classified correctly and an F1 score of 0.90. When tile drainage area is aggregated to the county scale, it agreed well (r2 = 0.69) with the reported area from the Ag Census. We found that Land Surface Temperature (LST) along with climate- and soil-related features were the most important factors for classification. The top-ranked feature is the median summer nighttime LST, followed by median summer soil moisture percent. This study demonstrates the potential of applying satellite remote sensing to map spatially explicit agricultural tile drainage across large regions. The results should be useful for land use change monitoring and hydrologic and nutrient models, including those designed to achieve cost-effective agricultural water and nutrient management strategies. The algorithms developed here should also be applicable for other remote sensing mapping applications.

Ironstone and red mud barriers to reduce subsurface movement of soil phosphorus

Loss of phosphorus in seepage may contribute to eutrophication of downstream water bodies. This study examined the potential use of pedogenic ironstone and untreated red mud (bauxite refining residue) as P sorbents in a permeable reactive barrier (PRB) to mitigate such loss. Effects of ironstone and red mud on P sorption (batch), transport (columns), saturated hydraulic conductivity (KS), and growth of common bermudagrass (Cynodon dactylon; greenhouse) were examined. Both materials had sorption maxima of ∼30 mmol P kg-1 or about five times that of a P-enriched sandy soil; however, sorption by red mud greatly increased with decreasing pH. Transport of P through columns of ironstone and red mud (diluted with nonreactive sand) was similar and slower compared to soil + sand. However, when red mud was mixed with soil, increased sorption at lower pH resulted in greater P retention compared to ironstone + soil (76% vs. 13%). Although addition of ironstone to soil up to 20% did not reduce KS, red mud at even 5% did. Soil amendment with red mud increased bermudagrass growth and P uptake. Given long-term neutralization of red mud in an acidic soil and increased P sorption, it may be suitable in a PRB if incorporated at a low rate and/or co-incorporated with a coarser material.

Soil Microbiomes and their Arsenic Functional Genes in Chronically High-Arsenic Contaminated Soils

Microbial arsenic transformations play essential roles in controlling pollution and ameliorating risk. This study combined high-throughput sequencing and PCR-based approaches targeting both the 16 S rRNA and arsenic functional genes to investigate the temporal and spatial dynamics of the soil microbiomes impacted by high arsenic contamination (9.13 to 911.88 mg/kg) and to investigate the diversity and abundance of arsenic functional genes in soils influenced by an arsenic gradient. The results showed that the soil microbiomes were relatively consistent and mainly composed of Actinobacteria (uncultured Gaiellales and an unknown_67 - 14 bacterium), Proteobacteria, Firmicutes (particularly, Bacillus), Chloroflexi, and Acidobacteria (unknown_Subgroup_6). Although a range of arsenic functional genes (e.g., arsM, arsC, arrA, and aioA) were identified by shotgun metagenomics, only the arsM gene was detected by the PCR-based method. The relative abundance of the arsM gene accounted for 0.20%-1.57% of the total microbial abundance. Combining all analyses, arsenic methylation mediated by the arsM gene was proposed to be a key process involved in the arsenic biogeochemical cycle and mitigation of arsenic toxicity. This study advances our knowledge about arsenic mechanisms over the long-term in highly contaminated soils.

New phosphorus losses via tile drainage depend on fertilizer form, placement, and timing

Agricultural phosphorus (P) losses are harmful to water quality, but knowledge gaps about the importance of fertilizer management practices on new (recently applied) sources of P may limit P loss mitigation efforts. Weighted regression models applied to subsurface tile drainage water quality data enabled estimating the new P losses associated with 155 P applications in Ohio and Indiana, USA. Daily discharge and dissolved reactive P (DRP) and total P (TP) loads were used to detect increases in P loss following each application which was considered new P. The magnitude of new P losses was small relative to fertilizer application rates, averaging 79.3 g DRP ha-1 and 96.1 g TP ha-1 , or <3% of P applied. The eight largest new P losses surpassed 330 g DRP ha-1 or 575 g TP ha-1 . New P loss mitigation strategies should focus on broadcast liquid manure applications; on average, manure applications caused greater new P losses than inorganic fertilizers, and surface broadcast applications were associated with greater new P losses than injected or incorporated applications. Late fall applications risked having large new P losses applications. On an annual basis, new P contributed an average of 14% of DRP and 5% of TP losses from tile drains, which is much less than previous studies that included surface runoff, suggesting that tile drainage is relatively buffered with regard to new P losses. Therefore old (preexisting soil P) P sources dominated tile drain P losses, and P loss reduction efforts will need to address this source.

A review of the development and implementation of the critical source area concept: A reflection of Andrew Sharpley's role in improving water quality

Critical source areas (CSAs) are small areas of a field, farm, or catchment that account for most contaminant loss by having both a high contaminant availability and transport potential. Most work on CSAs has focused on phosphorus (P), largely through the work in the 1990s initiated by Dr. Sharpley and colleagues who recognized the value in targeting mitigation efforts. The CSA concept has been readily grasped by scientists, farmers, and policymakers across the globe. However, experiences and success have been mixed, often caused by the variation in where and how CSAs are defined. For instance, analysis of studies from 1990 to 2023 shows that the proportion of the annual contaminant load coming from a CSA decreases from field to farm to catchment scale. This finding is consistent with increased buffering of CSAs and greater contribution of other sources with scale, or variation in the definition of CSAs. We therefore argue that the best application of CSAs to target mitigation actions should be at small areas that truly account for most contaminant loss. This article sheds light on the development and utilization of CSAs, paying tribute to Dr. Sharpley's remarkable contributions to the improvement of water quality, and reflecting upon where the CSA concept has succeeded or not in reducing contaminant (largely P) loss.

The importance of consensus science to managing phosphorus in the environment: SERA-17 and the legacy of Andrew Sharpley

Phosphorus (P) loss from agricultural systems to surface waters, and ultimately, eutrophication, presents a wicked problem requiring transdisciplinary solutions. The mission of SERA-17 (Southern Extension and Research Advisory Information Exchange Group-17) has been to address this problem by developing "Innovative Solutions to Minimize Phosphorus Losses from Agriculture." Over the course of his career, Dr. Andrew Sharpley demonstrated a rare ability to collaboratively achieve consensus around issues related to the science and management of P. The SERA-17 organization served as the central community of experts and stakeholders where that consensus was built and applied. The consensus-based approach, demonstrated by Sharpley and at the core of the SERA-17 organization, was routinely applied to key areas of P science to produce applied outcomes that have been readily adopted: advance foundational science to resolve knowledge gaps and to promote innovation; promote consistency in methods to facilitate comprehensive investigations and conclusions across a diversity of systems; engage diverse stakeholders to prioritize research, and ultimately, ensure that outcomes reflect a plurality of perspectives; and deliver pragmatic solutions that reflect the best information available at a particular time. We review the history of SERA-17 in delivering new science and management recommendations for P, with an eye to elucidating Sharpley's role and legacy in this process.

Citizen scientists filling knowledge gaps of phosphate pollution dynamics in rural areas

In situ monitoring is fundamental to manage eutrophication in rivers and streams. However, in recent decades, the frequency and spatial coverage of regulatory monitoring have often been reduced due to funding and infrastructure limitations. This reduction has made it impossible to provide adequate coverage for most water bodies. In this study, trained citizen scientists filled spatial and temporal gaps in agency monitoring across a major catchment in rural England. By integrating data from citizen scientists, regulatory agencies, and the local water company, it was possible to demonstrate the opportunities for hypothesis-based citizen scientist monitoring to identify continuous and event-driven sources of phosphate pollution. Local citizen scientists effectively covered important spatial gaps, investigating river conditions both upstream and downstream of suspected pollution point sources, improving the identification of their temporal dynamics. When combined with long-term monitoring data from regulatory agencies, it became possible to identify areas within the catchment that exhibited increased phosphate concentrations during periods of low river discharge (summer). Inter-annual trends and anomaly detection suggested that continuous pollution sources dominated over event-driven sources in many sub-basins, allowing for the prioritisation of mitigation actions. This study highlights the opportunity for citizen scientists to fill gaps in regulatory monitoring efforts and contribute to the improved management of eutrophication in rural catchments.

The role of agricultural drainage, storm-events, and natural filtration on the biogeochemical cycling capacity of aquatic and sediment environments in Lake Erie's drainage basin

Lake Erie is the most at risk of the Great Lakes for degraded water quality due to non-point source pollution caused by agricultural activities in the lake's watershed. The extent and temporal patterns of nutrient loading from these agricultural activities is influenced by the timing of agronomic events, precipitation events, and water flow through areas of natural filtration within the watershed. Downstream impacts of these nutrient loading events may be moderated by the co-loading of functionally relevant biogeochemical cycling microbial communities from agricultural soils. This study quantified loading patterns of these communities from tile drain sources, assessed whether functional communities from agricultural sources influenced downstream microbial functionality, and investigated how distance from agricultural sources, storm events, and areas of natural filtration altered nutrient cycling and nutrient fluxes in aquatic and sediment environments. Water and sediment samples were collected in the Wigle Creek watershed in Ontario, from tile drains through to Lake Erie, from May to November 2021, and microbial nitrogen (N) and phosphorous (P) cycling capacity (quantitative PCR), and nutrient levels were evaluated. Results showed that N and P functional groups were co-loaded with nutrients, with increased loading occurring during storm events and during agricultural activities including fertilization and harvest. Overall functional capacity in the aquatic environment decreased with distance from the agricultural sources and as water transited through natural filtration areas. In contrast, the sediment environment was more resilient to both agricultural disturbances and abiotic factors. This study expands our understanding of when and where different stages of N and P cycling occurs in agriculturally impacted watersheds, and identifies both seasons and regions to target with nutrient mitigation strategies.

A review of ecosystem services from edge-of-field practices in tile-drained agricultural systems in the United States Corn Belt Region

Edge-of-field management practices that reduce nutrient pollution from tile drainage while contributing habitat and other ecosystem services are needed to enhance agricultural systems in the US Corn Belt Region. In this review, we identified edge-of-field and catchment scale agricultural conservation practices for intercepting and treating tile drainage. The reviewed conservation practices were (1) controlled drainage, also known as drainage water management (USDA-NRCS Code 554); (2) drainage water recycling (USDA-NRCS Code 447); (3) denitrifying bioreactors (USDA-NRCS Code 605); (4) saturated buffers (USDA-NRCS Code 604); and (5) constructed or restored wetlands designed for water quality improvement (USDA-NRCS Code 656) herein referred to as water quality wetlands. We examined 119 studies that had information on one or more of the following ecosystem services: water retention, water quality improvement (e.g., nitrate, phosphate, sediment, or pesticide retention), wetland habitat (for birds, aquatic organisms, and pollinators), crop yield improvement, and other benefits (e.g., recreation, education, aesthetic appreciation, greenhouse gas retention). We found the five edge-of-field practices were all effective at removing nitrate with varying degrees of other potential benefits and disservices (e.g., greenhouse gas production). Drainage water recycling and water quality wetlands have the potential to provide the most co-benefits as they provide surface water systems for capturing surface flows in addition to tile drainage while also potentially providing habitat and recreation opportunities. However, the following research needs are identified: 1) the disservices and benefits associated with drainage water recycling have not been adequately evaluated; 2) surface flow dynamics are understudied across all reviewed management practices; 3) a complete accounting of phosphorus species and flow pathways for all management practices is needed; 4) field evaluations of the habitat benefit of all management practices are needed; and 5) greenhouse gas dynamics are understudied across all management practices. While all management practices are expected to reduce nitrate loads, addressing these knowledge gaps will help inform holistic management decisions for diverse stakeholders across the US Corn Belt.

Long-term phosphorus removal by Ca and Fe-rich drainage filter materials under variable flow and inlet concentrations

Phosphorus (P) losses from tile-drained agricultural fields may degrade surface water quality by accelerating eutrophication. Among the different edge-of-field technologies, compact filter systems using different filter materials have been identified as potentially effective solutions for removing P from drainage water before discharge downstream. This study investigated the long-term (>696 days) P removal efficiency of 5 different filter materials in a column setup, using artificial drainage water (pH 6). Filter materials included two iron-based granulates (calcinated diatomaceous earth (CDE), ferric hydroxide granules (CFH)), and three calcium-based granulates (seashells, limestone, calcinated silicate/calcium oxide (Filtralite-P)). Experiments were performed under variable flow rates (0.037 and 1.52 L h-1; hydraulic retention time of 26-43 min and 18-30 h) and inlet P concentrations (0.14 and 0.7 mg L-1). An overall analysis revealed that the Fe-based materials achieved higher P retention than Ca-based materials. In particular, CFH was capable of retaining 99 and 98 % of the high and low inlet P concentrations, respectively. Conversely, limestone retained only 25 % of the high P load. CDE performed moderately well, independently of the inlet P concentration. Filtralite-P and Seashells performed well at high inlet P concentration but relatively poorly at low P concentration. The sensitivity of filter material P removal efficiency to variations in P loading was generally lowest for CFH and highest for limestone.

Trade-offs in nutrient and sediment losses in tile drainage from no-till versus conventional conservation-till cropping systems

Nutrient and soil loss from agricultural areas impairs surface water quality globally. In the Great Lakes region, increases in the frequency and magnitude of harmful and nuisance algal blooms in freshwater lakes have been linked to elevated phosphorus (P) losses from agricultural fields, some of which are transported via tile drainage. This study examined whether concentrations and loads of P fractions, total suspended sediments (TSS), nitrate (NO3 - ), and ammonium (NH4 + ) in tile drainage in a clay soil differed between a continuous no-till system combining cover crops and surface broadcast fertilizer (no-till cover crop [NTCC]), and a more conventional tillage system with shallow tillage, fertilizer incorporation and limited use of cover crops (conventional conservation-till, CT). Both sites had modest soil fertility levels. Year-round, high-frequency observations of tile drainage flow and chemistry are described over 4 full water years and related to management practices on the associated fields. There were similar water yields in tile drainage between the two systems; however, losses of TSS, particulate P (PP), and NO3 - were consistently greater from the CT site, which received larger quantities of fertilizer. In contrast, dissolved reactive P (DRP) losses were considerably greater from the NTCC site, offsetting the lower PP losses, such that there was little difference in TP losses between sites. Approximately 60% of the DRP losses from the NTCC site over the 4 years were associated with incidental losses following surface application of fertilizer in fall. This study provides insight into trade-offs in controlling losses of different nutrient fractions using different management systems.

Linking transport pathways and phosphorus distribution in a loamy soil: a case study from a North-Eastern German Stagnosol

Heterogeneous flow pathways through the soil determine the transport of dissolved and particle-bound nutritional elements like phosphorus (P) to ground and surface waters. This study was designed to understand the spatial patterns of P in agriculturally used soils and the mechanisms causing P accumulation and depletion at the centimetre scale. We conducted dye tracer experiments using Brilliant Blue on a loamy Stagnosol in North-Eastern-Germany. The plant-available P was analysed using double lactate extraction (DL-P). The plant-available P content of the topsoil was significantly higher than that of the subsoil in all three replicates (p < 0.001). The topsoil's stained areas showed significantly higher P contents than unstained areas (p < 0.05), while the opposite was found for the subsoil. The P content varied enormously across all observed soil profiles (4 to 112 mg P kg-1 soil) and different categories of flow patterns (matrix flow, flow fingers, macropore flow, and no visible transport pathways). The P contents of these transport pathways differed significantly and followed the order: Pmatrix flow > Pfinger flow > Pno visible transport pathways > Pmacropore flow. We conclude that P tends to accumulate along flow pathways in the topsoil in the observed fertilized and tilled mineral soil. In contrast, in the subsoil at a generally lower P level, P is depleted from the prominent macroporous flow domains.

Contribution of non-point source pollution that migrated with underground runoff process based on the SWAT model and a digital filter algorithm

Non-point source (NPS) pollution has always been the focus of research worldwide, and understanding the migration process is the basis for effective control of NPS pollution. In this study, the SWAT model and digital filtering algorithm were combined to explore the contribution of NPS pollution that migrated with underground runoff (UR) process to the Xiangxi River watershed. The results showed that the surface runoff (SR) was the main migration process of NPS pollution, while the contribution of NPS pollution that migrated with the UR process only accounted for 30.9%. With the decrease in annual precipitation among the three selected hydrological years, the proportion of NPS pollution that migrated with the UR process for TN decreased, whereas the proportion for TP increased. The contribution of NPS pollution migrated with UR process varied remarkably during different months. Although the maximum total load and the load of NPS pollution that migrated with the UR process for TN and TP all appeared in the wet season, due to the hysteresis effect, the load of NPS pollution that migrated with the UR process for TP appeared 1 month later than the total load of NPS pollution. With an increase in precipitation from the dry season to the wet season, the proportion of NPS pollution that migrated with the UR process for TN and TP decreased gradually, and the degree of decrease in NPS pollution that migrated with the UR process for TP was more evident than that for TN. Besides, being affected by topography, land use, and other factors, the proportion of NPS pollution that migrated with the UR process for TN decreased from 80% in upstream areas to 9% in downstream areas, while that for TP reached a maximum of 20% in downstream areas. Based on the research results, the contribution of soil and groundwater cumulative nitrogen and phosphorus should be considered, and different managements and control measures for different migration routes should be adopted in controlling pollution.

Challenges in linking soil health to edge-of-field water quality across the Great Lakes basin

To better understand agricultural nutrient losses, we evaluated relationships between management (e.g., manure and tillage), soil health measurements, and resulting edge-of-field (EOF) surface water quality. This work was conducted before or early into conservation implementation at 14 Great Lakes Restoration Initiative EOF sites spanning Wisconsin, Michigan, Indiana, Ohio, and New York. Analyses of site characteristics (hydroclimate, management, catchment properties) along with 3 yr of soil health measurements (chemical, biological, and physical properties) showed EOF-nutrient export depended on both site and soil properties. A pattern emerged whereby sites not receiving manure and sites with manure defined opposite ends of several gradients for soil and water data. Sites receiving manure had increased microbial activity, organic matter (3.2 vs. 2.7%), and soil test phosphorus (P) (2.8 times more) relative to sites without manure. Suspended sediments (SS), total P (TP), and total nitrogen (TN) in EOF surface runoff varied over three to five orders. Multivariate analysis among sites showed covariant linkages between soil nutrients, soil C, microbial properties, and nutrients in EOF water. There were positive univariate relationships between water-extractable soil P and annual EOF-water concentrations and yields of orthophosphate, TP, TN, and SS (p < .01). Some soil physical properties (e.g., bulk density and infiltration) also covaried among sites but were not consistently related to runoff index or water yield variables. Given the observed among-site variability, we were not able to isolate desirable soil health signals on EOF surface water quality.

Phosphorus sources, forms, and abundance as a function of streamflow and field conditions in a Maumee River tributary, 2016-2019

Total phosphorus (TP), dissolved P (DP), and suspended sediment (SS) were sampled in Black Creek, Indiana, monthly during base flow and for 100 storm events during water years 2016-2019, enabling analysis of how each of these varied as a function of streamflow and field conditions at nested edge-of-field sites. Particulate P was normalized for SS (P_SS  = [TP - DP]/SS). Streamflow events were differentiated by maximum TP concentrations co-occurring with maximum SS (SED) or DP (SOL). The combination of new precipitation and high antecedent soil-water storage during months when fields were exposed coincided with higher streamflow that drove SED events. These SED events carried more SS, including sediment eroded from streambanks that added sediment P but also may have provided for sorption of DP. During SOL events, DP was higher and contributed approximately half of TP; SS was lower. These SOL events had higher P_SS , more similar to that in base flow as well as composited samples of overland flow and tile-drain discharge from fields. Base-flow samples had significantly higher P_SS concentrations than most event samples, with ≤25 times enrichment relative to soil P concentrations in fine-grained source material. Combining base-flow and event samples showed that P_SS integrates SS, DP, and streamflow. Addition of new suspended sediment during events may provide for sorption of DP during and after events and storage in the system, delaying delivery of this P to Lake Erie relative to what would be expected for the dissolved form but adding to the legacy P stored in the stream system.

Investigating relationships between climate controls and nutrient flux in surface waters, sediments, and subsurface pathways in an agricultural clay catchment of the Great Lakes Basin

Water quality within agricultural catchments is governed by management practices and climate conditions that control the transport, storage, and exchange of nutrients between components of the hydrologic cycle. This study aims to improve knowledge of nitrogen (N) and phosphorus (P) transport in low permeability agricultural watersheds by considering spatial and temporal trends of surface water nutrient concentrations in relation to hydroclimatic drivers, sediment quality, shallow hyporheic exchange, groundwater quality, and tile drain discharge over a 14-month field study in a clay hydrosystem of the Lake Huron basin, one of the five Great Lakes. Results found that events of varying magnitude and intensity enhanced nutrient release from overland flow and subsurface pathways. Tile drain discharge was found to be a consistent and elevated source of P and N to surface waters when flowing, mobilizing both diffuse nutrients from fertilizer application and legacy stores in the vadose zone. Surface water quality was temporally variable at the seasonal and event scale. Targeted sampling following fertilization periods, snowmelt, and moderate precipitation events revealed catchment wide elevated nutrient concentrations, emphasizing the need for targeted sampling regimes. Controls other than discharge magnitude and overland flow were found to contribute to peak nutrient concentrations, including internal nitrate loading, soil-snowmelt interaction, catchment wetness, and freeze thaw cycles. Sediments were found to store P in calcium minerals and have a high P storage capacity. Instream mechanisms such as sediment P fixation and hyporheic exchange may play a role in mediating surface water quality, but currently have no discernable benefit to year-round surface water nutrient concentrations. Best management practices need to focus on reducing sources of agricultural nutrients (e.g., field phosphorus concentrations and tile drain discharge loading) at the watershed scale to reduce nutrient concentrations and export in flashy clay catchments.

Phosphorus and iron-oxide transport from a hydrologically isolated grassland hillslope

Dissolved reactive phosphorus (DRP) loss from agricultural soils can negatively affect water quality. Shallow subsurface pathways can dominate P losses in grassland soils, especially in wetter months when waterlogging is common. This study investigated the processes controlling intra- and inter-event and seasonal DRP losses from poorly drained permanent grassland hillslope plots. Temporal flow related water samples were taken from surface runoff and subsurface (in-field pipe) discharge, analysed, and related to the likelihood of anaerobic conditions and redoximorphic species including nitrate (NO₃^-) over time. Subsurface drainage accounted for 89% of total losses. Simple linear regression and correlation matrices showed positive relationships between DRP and iron and soil moisture deficit; and negative relationships between these three factors and NO₃^- concentrations in drainage. These data indicate that waterlogging and low NO₃^- concentrations control the release of P in drainage, potentially via reductive dissolution. The relationship between DRP and metal release was less obvious in surface runoff, as nutrients gathered from P-rich topsoil camoflaged redox reactions. The data suggest a threshold in NO₃^- concentrations that could exacerbate P losses, even in low P soils. Knowledge of how nutrients interact with soil drainage throughout the year can be used to better time soil N and P inputs via, for example, fertiliser or grazing to avoid to excessive P loss that could harm water quality.

Effectiveness of denitrification bioreactors with woodchips, corn stover, and phosphate-sorbing media for simultaneous removal of drainage water N and P in a corn-soybean system

Millions of acres of farmland in the midwestern United States (US) are artificially drained, and this contributes to the export of nitrogen (N) and phosphorus (P) from agricultural land to surface water. Using a 36-acre tile-drained farm field, effects of P-sorbing media in combination with a denitrifying bioreactor system constructed with woodchips (WC) and corn stover (CS) on reducing nutrient export in drainage water were tested for 3 cropping years (2018-2020). The field was divided into three subfields as replicates. In each subfield, the drainage water was divided and separately channeled into three bioreactors, each of which contains one of the three different substrates: WC, CS, and CS-WC (1:1 v/v mixture of CS and WC), randomly assigned. The outlet of each compartment contained a 2.25 L flow-through chamber filled with activated iron (Fe) filings as P-sorbing material. Both WC and CS bioreactors were effective in removing drainage NO₃ ^- with a 77% (WC), 86% (CS), and 89% (CS-WC) reduction in mean NO₃ ^- -N concentration. For the three cropping years, the WC bioreactor reduced the total drainage inorganic N (NO₃ ^- -N + NH₄ ⁺ -N) load by 72%, but the CS bioreactor increased the total inorganic N load in the drainage water due to the substantial release of NH₄ ⁺ with the decomposition of CS. The breakdown of CS also increased drainage P. The NH₄ ⁺ and P release decreased with the decrease in the proportion of CS; thus, not more than 10% of CS is recommended for blending with WC to enhance the performance of a bioreactor. The P-sorbing Fe filing media reduced the P loads in drainage by an average of 19% during the 2-year study.

Resolving new and old phosphorus source contributions to subsurface tile drainage with weighted regressions on discharge and season

Agricultural losses of dissolved reactive phosphorus (DRP) emanate from both historic P applications (i.e., "old P") and recently applied fertilizer (i.e., "new P"). Understanding the relative contributions of these sources is important for mitigating DRP losses from agriculture. This study provides a proof-of-concept for resolving new P vs. old P source contributions to DRP losses in subsurface tile drainage using edge-of-field water quality data and management records from eight fields in Ohio. Weighted regressions on discharge and season (WRDS) were fitted using data from periods without P fertilizer applications and then used to predict DRP losses in tile drainage during new P loss risk periods (default length, 90 d) after fertilizer applications. Differences between observed and predicted DRP concentrations during the new P loss risk period were attributed to the new P source. Remaining losses were attributed to the old soil P source. The WRDS model performance was modest (modified Kling-Gupta efficiency ranged from -0.074 to 0.484). New P sources contributed between 0 and 17% of overall DRP losses (average, 7%), with old soil P contributing 83-100%. Individual P fertilizer applications were associated with new DRP losses up to 192 g P ha^-1 . Increasing the length of the risk period for new P losses up to 180 d after fertilizer application marginally increased the estimated contribution of the new P source. The WRDS-based analysis provides a novel approach for resolving the contributions of new and old sources to edge-of-field DRP losses.

Analysis of reactive phosphorus treatment by filter materials at the edge of tile-drained agricultural catchments: A global view of the current status and challenges

Phosphorus losses from agriculture have long generated concern due to the ecological impact on surface waters. Here tile-drained agricultural catchments are a critical source for concentrating and transporting phosphorus bioavailable forms or dissolved reactive phosphorus (DRP). Hence, edge-of-field technologies have been introduced to reduce DRP loads. Filter systems have received special attention due to their targeted approach using a permeable filter material (FM) rich in DRP sorbents. This review explores the performance and applicability of FMs in the aforementioned context because of the growing number of studies. An overall analysis revealed that sorption is preferable to precipitation for DRP retention at the edge-of-field, and that FM pH and particle size affect sorption properties and subsequently DRP retention and lifetime. Thus, FMs with predominant amounts of iron and/or aluminium can be recommended. Such materials generally have an appreciable availability of DRP binding sites, strong bonds with DRP and short reaction times, as well as low desorption, which lead to good operation. On the other hand, FMs with predominant amounts of calcium and/or magnesium are restricted to catchments with favourable conditions unless they have optimal reactivity for DRP. The review also found that hydraulic retention time plays a key role in the performance and applicability of FMs, especially in those dependent on precipitation reactions. Therefore, it is crucial that FMs are designed, constructed and managed according to the catchment conditions-including normally varying flow rates and DRP concentrations-in order to ensure successful operation. This reflects in long-term, high and steady net DRP retention along with low costs, thus improving the FM cost-effectiveness, besides discharging non-harmful effluents to aquatic ecosystems.

The seasonality of nitrate and phosphorus leaching from manure and chemical fertilizer added to a chernozemic soil in Canada

Identifying seasons sensitive to nutrient losses could help farmers and policymakers to formulate effective nutrient loss reduction strategies. This long-term study monitored water percolation as well as nitrate (NO₃ -N) and total phosphorus (TP) leaching from liquid swine manure and chemical fertilizer applied to intact core lysimeters in a sandy loam soil in Manitoba, Canada. Water percolation, NO₃ -N, and TP leaching were monitored from 2005 to 2016. Chemical fertilizer showed greater average annual mean water percolation (p = .01), annual flow-weighted mean concentration (FWMC) of NO₃ -N (22 mg L^-1 ; p < .001), and annual NO₃ -N leaching (36 kg N ha^-1 ; p = .002) compared with the manure treatment (FWMC NO₃ -N, 15 mg L^-1 ; NO₃ -N leaching load, 22 kg N ha^-1 ). Average annual mean TP loss did not differ between treatments (p = .86). Spring (April-June) was the most sensitive season, when >75% of annual percolation, >80% of annual NO₃ -N, and >68% of annual TP leaching losses occurred from both manure and chemical fertilizer. Annual NO₃ -N and TP leaching increased exponentially with cumulative winter and spring precipitation (control, r²  = .69; manure, r²  = .79; chemical fertilizer, r²  = .63) and decreased with winter and spring air temperatures. The largest spring NO₃ -N and TP leaching losses were observed in 2013, which followed the dry year of 2012, indicating the potential for nutrient flushing. The findings emphasize the need for environmentally sound N and P management strategies in cold North American regions underlain by coarse-textured soils, particularly during the spring season.

Assessing the concept of control points for dissolved reactive phosphorus losses in subsurface drainage

Agricultural phosphorus (P) loss, which is highly variable in space and time, has been studied using the hot spot/hot moment concept, but increasing the rigor of these assessments through a relatively newer "ecosystem control point" framework may help better target management practices that provide a disproportionate water quality benefit. Sixteen relatively large (0.85 ha) subsurface drainage plots in Illinois were used as individual observational units to assess dissolved reactive P (DRP) concentrations and losses within a given field over four study years. Three plot-months were identified as DRP control points (one export and two transport control points), where each plot-month contributed >10% of the annual DRP load from the field. These control points occurred on separate plots and in both the growing and nongrowing seasons but were likely related to agronomic P applications. Elevated soil test P, especially near a historic farmstead, and soil clay content were spatial drivers of P loss across the field. The nongrowing season was hypothesized to be the most significant period of P loss, but this was only documented in two of the four study years. A cereal rye (Secale cereale L.) cover crop did not significantly reduce DRP loss in any year, but there was also no evidence of increased drainage P losses due to freezing and thawing of the cover crop biomass. This work confirmed annual subsurface drainage DRP losses were agronomically small (<3% of P application rate), although the range of DRP concentrations relative to eutrophication criteria still demonstrated a potential for negative environmental impact. The control point concept may provide a new lens to view drainage DRP losses, but this framework should be refined through additional within-field studies because mechanisms of P export at this field were more nuanced than just the presence of tile drainage (i.e., a transport control point).

Updates to the Annual P Loss Estimator (APLE) model

The Annual P Loss Estimator (APLE) is a spreadsheet-based model developed for predicting annual field-scale P loss in surface runoff and changes in soil test P. This empirically based model was designed for use by those without significant modeling experience. However, a significant limitation with the model is that it does not calculate runoff. Moreover, APLE is deterministic and thus predicts a single value for a given set of inputs, thereby ignoring any uncertainties associated with model inputs. Here, we describe modifications to APLE that allow users to estimate runoff using the Curve Number method. Using Monte Carlo simulations, the updated version of APLE also provides users the ability to account for model input uncertainties in estimating model prediction errors. We provide examples of using the revised version of APLE (ver. 3.0) for calculating P loss from two fields in Mississippi over a 4-yr period and calculating the change in Mehlich-3 P concentrations over a 9-yr period at three locations in Maryland following cessation of P application. Both examples demonstrate that incorporating estimates of uncertainties in both measured data and model predictions provides modelers with a more realistic understanding of the model's performance.

Contribution of nonpoint source pollution from baseflow of a typical agriculture-intensive basin in northern China

Baseflow is an essential component of total surface runoff that is widely considered one of the most influential factors regarding water quality via nonpoint source (NPS) pollution. Previously, many researchers and policy makers have directed their efforts toward surface runoff pollution, largely ignoring nutrient delivery via baseflow. Taking a typical agriculture-intensive basin of northern China as an example, this study explored the spatiotemporal characteristics of baseflow and pollution load in relation to NPS pollution. Baseflow was quantified using digital filtering techniques, and the results together with observed pollution data were used to validate a physically based hydrological model, i.e., the Soil and Water Assessment Tool. Then, the spatial and temporal distribution characteristics of NPS and baseflow pollution were investigated using the modeling results. Results indicated that baseflow contribution to total runoff accounted for more than 70% during the studied years (2016-2018), and 84.15% of the basin area showed non-point source pollution dominated by baseflow pollution; both baseflow and its pollution load were greater in the nonflood seasons (spring, autumn, and winter) than in the flood season (summer); the spatial distribution of baseflow total nitrogen and total phosphorus pollution intensity showed higher values in the east and lower values in the west; the scaling effects of baseflow and its pollution load was found with increasing basin area. The results of our study highlighted the necessity for management of pollution load via baseflow in the river basin and provided reference information for improvement of NPS pollution management in other similar basins.

Sensitivity of the Transport of Plastic Nanoparticles to Typical Phosphates Associated with Ionic Strength and Solution pH

The influence of phosphates on the transport of plastic particles in porous media is environmentally relevant due to their ubiquitous coexistence in the subsurface environment. This study investigated the transport of plastic nanoparticles (PNPs) via column experiments, paired with Derjaguin–Landau–Verwey–Overbeek calculations and numerical simulations. The trends of PNP transport vary with increasing concentrations of NaH₂PO₄ and Na₂HPO₄ due to the coupled effects of increased electrostatic repulsion, the competition for retention sites, and the compression of the double layer. Higher pH tends to increase PNP transport due to the enhanced deprotonation of surfaces. The release of retained PNPs under reduced IS and increased pH is limited because most of the PNPs were irreversibly captured in deep primary minima. The presence of physicochemical heterogeneities on solid surfaces can reduce PNP transport and increase the sensitivity of the transport to IS. Furthermore, variations in the hydrogen bonding when the two phosphates act as proton donors will result in different influences on PNP transport at the same IS. This study highlights the sensitivity of PNP transport to phosphates associated with the solution chemistries (e.g., IS and pH) and is helpful for better understanding the fate of PNPs and other colloidal contaminants in the subsurface environment.

Influence of climate, topography, and soil type on soil extractable phosphorus in croplands of northern glacial-derived landscapes

Delineating the relative solubility of soil phosphorus (P) in agricultural landscapes is essential to predicting potential P mobilization in the landscape and can improve nutrient management strategies. This study describes spatial patterns of soil extractable P (easily, moderately, and poorly soluble P) in agricultural landscapes of the Red River basin and the southern Great Lakes region. Surface soils (0-30 cm) and select deeper cores (0-90 cm) were collected from 10 cropped fields ranging in terrain (near-level to hummocky), soil texture (clay to loam), composition (calcareous to noncalcareous), and climate across these differing glacial landscapes. Poorly soluble P dominated (up to 91%) total extractable P in the surface soils at eight sites. No differences in the relative solubilities of soil extractable P with microtopography were apparent in landscapes without defined surface depressions. In contrast, in landscapes with pronounced surface depressions, increased easily soluble P (Sol-P), and decreased soil P sorption capacity were found in soil in wetter, low-slope zones relative to drier upslope locations. The Sol-P pool was most important to soil P retention (up to 28%) within the surface depressions of the Red River basin and at sites with low-carbonate soils in the southern Lake Erie watershed (up to 28%), representing areas at elevated risk of soil P remobilization. This study demonstrates interrelationships among soil extractable P pools, soil development, and soil moisture regimes in agricultural glacial landscapes and provides insight into identifying potential areas for soil P remobilization and associated P availability to crops and runoff.

Effect of water sampling strategies on the uncertainty of phosphorus load estimation in subsurface drainage discharge

Accurate phosphorus (P) load estimation in subsurface drainage water is critical to assess the field-scale efficacy of conservation practices. The HydroCycle-PO₄ instrument measures real-time total reactive P (TRP) concentration without the need for sample filtration, thereby enabling comparative evaluation of different sampling strategies. The main objective of this study was to evaluate the effects of water sampling strategies on the uncertainty of P load estimation. Hourly TRP concentration and hourly drainage discharge measurements formed the reference P load dataset. Four hypothetical water sampling strategies were evaluated: (a) time-proportional discrete sampling, (b) time-proportional composite sampling, (c) flow-proportional discrete sampling, and (d) flow-proportional composite sampling. All sampling strategies underestimated TRP load compared with the reference dataset. Total reactive P load underestimation changed from 0.2 to 51% as time-proportional discrete sampling intervals increased from 3 h to 14 d. Total reactive P load underestimation changed from 12 to 43% as the time-proportional compositing scenario increased from 1 to 7 d, each with one aliquot per day. In the case of flow-proportional discrete sampling scenario, the lowest (0.6%) and the highest (-5.1%) uncertainties were observed when 1- and 5-mm flow intervals were used. The relative error based on the results provided by the flow-proportional composite sampling ranged from 0.2% when using 1-mm flow interval to -6.7% when using 5-mm flow interval. In conclusion, the flow-proportional sampling strategies provided a more accurate estimate of cumulative P load with fewer number of samples because a greater portion of samples were taken at higher flow rates compared with time-proportional sampling strategies.

Simulation of Subsurface Drainage in the Sugarcane Crop under Different Spacing and Drain Depths

Agricultural land drainage is an instrument for growing production and a tool for the conservation of land resources. The performance of land drainage systems is thus critical for achieving sustainable agricultural production Recently, many types of software have been developed in this field for modeling and simulating the performance of these systems. SISDRENA is a simulation model of the performance of underground drainage systems. The main objectives of this paper are to simulate different combination of depths and spaces between drains and to analyze their impact on potential sugarcane productivity in the western plains of Venezuela using a land drainage system model. Therefore, three climatic scenarios were defined by annual precipitation: dry years (25% below average), normal (mean) and humid (75% above average). The scenarios were implemented in three different soil types: sandy loam, loam and silt loam, with a hydraulic conductivity of 0.19, 0.26 and 0.04 m day−1, respectively. The simulation of the yield related to soil deficit (YRD) and water stress (YRW) indicated that the highest yields were reached for the larger spacing between drains and the high conductivity hydraulic of soils. In relation to the average relative productivity (YT), it was shown that in soils with a greater water retention capacity there is an inversely proportional relationship between the spacing between drains and the productivity. We concluded that in order to reach the maximum sugarcane yield, the effect of hydraulic conductivity is more important than the changes in the precipitation pattern.

Use of iron-coated sand for removing soluble phosphorus from drainage water

Mitigation measures are needed for reducing chronic dissolved phosphorus (P) losses from agricultural soils with a legacy of excessive P inputs to surface waters. Since pipe drains are an important pathway for P transport from agricultural soils to surface waters in flat areas, removing P from drainage water can be an effective measure. During a 4.5 year-field experiment, we tested the performance of a pipe drain enveloped with Fe-coated sand for removing soluble P from drainage water. Iron-coated sand is a by-product of the drinking water industry and has a high ability to bind P. The P concentration in the effluent from the enveloped pipe drain remained at a very low level over the entire monitoring period, with a removal percentage amounting to 93% for total P. During the field experiment, the enveloped pipe drain was below the groundwater level for a prolonged time. Nevertheless, no reduction of Fe(III) in the Fe-coated sand occurred during the first two years, most likely due to preferential reduction of Mn oxides present in the coatings of the sand particles, as reflected in elevated effluent Mn concentrations. Thereafter, reductive dissolution of Fe oxides in the coatings caused a gradual increase in the Fe concentration in the enveloped pipe drain effluent over time. Concomitantly, the dissolved Mn concentration decreased, most probably due to the depletion in easily accessible Mn oxides in the Fe-coated sand. The Fe in the Fe-coated sand was identified as silicate-containing ferrihydrite (Fh). The submerged conditions of the enveloped pipe drain neither affected the stability of Fh in the Fe-coated sand nor the ability of this measure to capture P from drainage water. Enveloping pipe drains with Fe-coated sand is an effective method for reducing dissolved P inputs from agricultural soils to surface waters and holds great promise for implementation in practice.

Conditioning of Feed Material Prior to Feeding: Approaches for a Sustainable Phosphorus Utilization

A circular phosphorus (P) bioeconomy is not only worthwhile for conserving limited mineral P reservoirs, but also for minimizing negative environmental impacts caused by human-made alterations. Although P is an essential nutrient, most of the P in concentrates based on cereals, legumes and oilseed byproducts is organically bound to phytate. The latter cannot be efficiently utilized by monogastric animals and is therefore diluted into the environment through the manure pathway. This review examines various strategies for improved P utilization in animals and reflects the respective limitations. The strategies considered include feeding of debranned feedstuffs, pre-germinated feed, co-feeding of phytase and feeding material with high native phytase activity. All these approaches contribute to an improved P bioavailability. However, about half of the organic P content continues to be excreted and therefore remains unused by the animals. Nevertheless, technologies for an efficient utilization of P from cereal-based feed already exist; however, these are not industrially established. Conditioning feed material prior to feeding fosters P-reduced feed; meanwhile, P bound to phytate can be recovered. Based on known techniques for P separation and solubilisation from cereal products and phytate conversion, potential designs for feed material conditioning processes are proposed and evaluated.

External Phosphorus Loading in New Lakes

The water quality of lakes is highly dependent on external phosphorus (P) loading. The vast external loadings from sewage and other wastewater discharge that European lakes have historically received have been dramatically reduced today by improved wastewater treatment. Gaining knowledge of the catchment characteristics that influence external P-loading should enable predictions of the achievable water quality of lakes. In this study, we tested this proposition. Data from 90 new Danish lakes show no apparent relationship between the mean summer P-concentration and the size or land use of the catchments. The external P-loading and resulting annual P-concentration were further investigated on a representative subset of 12 of the new lakes, using six methods. Three of the methods used empirical estimates of P-transport from catchments, based on the national average P-transport, runoff-dependent P-transport, and crop-dependent P-transport, and the other three methods used different empirical models tested on the lakes. External P-loading was reliably predicted by several of the methods. The predictions of the annual P-concentration were highly dependent on the inclusion of annual runoff. However, the predicted P-concentrations were generally overestimated, most pronounced for the nutrient-poor and most recently established lakes. In these lakes, internal P-loading was found to be the most important factor in predicting achievable water quality.

SWAT model application for evaluating agricultural conservation practice effectiveness in reducing phosphorous loss from the Western Lake Erie Basin

Lake Erie is threatened by eutrophication and harmful algal blooms due to excess nutrient loading from agricultural sources. Agricultural conservation practices (ACPs) have been developed and implemented to reduce nutrient losses but estimating ACP effectiveness is challenging. The Soil and Water Assessment Tool (SWAT) has been used to investigate ACP effectiveness for water quality improvement. Many SWAT applications have been developed by different investigators to evaluate ACP effectiveness for reducing nutrient, particularly phosphorus (P), loading in the agriculturally-dominated Western Lake Erie Basin (WLEB). Our objective is to establish what has been achieved by past modeling research and make suggestions for future applications and improvements. We synthesized the findings of 28 SWAT modeling studies within the WLEB. Models generally performed satisfactorily against accepted criteria for streamflow and sediment, but performance for P loads, like soluble reactive P, was mostly "unsatisfactory". The "unsatisfactory" performance maybe due to imperfections and idealizations in model formulations and/or parameterization. Thus, simulations of P transport and transformation processes need improvement. In addition, model parameter selection is the key part of model set-up. Most SWAT modeling studies used default values during initial set-up, then performed calibration and validation. It was found that the calibrated P related parameter values varied widely across different studies, even within the same watershed with some values unrealistic for the study areas. The phenomena of different combinations of model parameters producing similar outputs indicates equifinality. Equifinality in the baseline model may impact results when ACPs are incorporated. Furthermore, the unrealistic values used in ACP assessment undermine the credibility of ACP effectiveness. Future model applications should try to re-examine the calibrated P parameters and make sure they are realistic for the study area as well as reduce equifinality by constraining the model with characterization of watershed conditions, better understanding of hydrologic processes, and parameter values based on real-world observations. In summary, future model applications should focus on improving P transport and transformation processes, using measured watershed characteristics for parameterization, and improving reflections of climate change, which could result in more accurate assessments of ACP effectiveness to meet targeted goals.

Desorption Kinetics of Legacy Soil Phosphorus: Implications for Non-Point Transport and Plant Uptake

Soil phosphorus (P) solubility and kinetics partly control dissolved P losses to surface water and uptake by plants. While previous studies have focused on batch techniques for measuring soil P desorption kinetics, flow-through techniques are more realistic because they simulate P removal from the system, akin to runoff, leaching, and plant uptake. The objectives were to measure soil P desorption by a flow-through technique at two flow rates and several batch methods, and utilize both for understanding how flow rate impacts the thermodynamics and kinetics of soil P desorption. Desorption obeyed first-order kinetics in two different phases: an initial rapid desorption phase followed by a gradual release. Desorption was limited by equilibrium and the kinetics of physical processes as demonstrated by an interruption test. Dilution-promoted desorption occurred with increasing cumulative volume, which increased desorption rate via first-order kinetics. The batch tests that simulated cumulative solution volume and time of flow-through were similar to the flow-through results; however, the batch methods overestimated the desorption rates due to less limitations to diffusion. Fast flow rates desorbed less P, but at a greater speed than slow flow rates. The differences were due to contact time, cumulative time, and solution volume, which ultimately controlled the potential for chemical reactions to be realized through physical processes. The interaction between these processes will control the quantity and rate of desorption that buffer P in non-point drainage losses and plant uptake.

Phosphorus availability and leaching losses in annual and perennial cropping systems in an upper US Midwest landscape

Excessive phosphorus (P) applications to croplands can contribute to eutrophication of surface waters through surface runoff and subsurface (leaching) losses. We analyzed leaching losses of total dissolved P (TDP) from no-till corn, hybrid poplar (Populus nigra X P. maximowiczii), switchgrass (Panicum virgatum), miscanthus (Miscanthus giganteus), native grasses, and restored prairie, all planted in 2008 on former cropland in Michigan, USA. All crops except corn (13 kg P ha⁻¹ year⁻¹) were grown without P fertilization. Biomass was harvested at the end of each growing season except for poplar. Soil water at 1.2 m depth was sampled weekly to biweekly for TDP determination during March–November 2009–2016 using tension lysimeters. Soil test P (0–25 cm depth) was measured every autumn. Soil water TDP concentrations were usually below levels where eutrophication of surface waters is frequently observed (> 0.02 mg L⁻¹) but often higher than in deep groundwater or nearby streams and lakes. Rates of P leaching, estimated from measured concentrations and modeled drainage, did not differ statistically among cropping systems across years; 7-year cropping system means ranged from 0.035 to 0.072 kg P ha⁻¹ year⁻¹ with large interannual variation. Leached P was positively related to STP, which decreased over the 7 years in all systems. These results indicate that both P-fertilized and unfertilized cropping systems may leach legacy P from past cropland management.

Assessment of water quality in ponds across the rural, peri-urban, and urban gradient

The Prairie Pothole Region is one of the most wetland rich areas of the world and has experienced intense disturbance from increased agricultural demands and urban sprawl. This study assessed ponds across the urban gradient for the first time in the region to determine the impacts of urbanization on water quality. Thirty ponds (ten rural, ten peri-urban, and ten urban) were randomly selected and compared based on land use type and the impervious to pervious surface ratio within 1.6 km of each pond. Water quality samples were taken monthly in 2015 and 2016, across 3 and 6 months respectively. Assessment included chemical and physical parameters, which were compared spatially across the gradient and temporally between sampling periods. Results indicate disturbance from urbanization negatively impacts water quality. Spatially across the gradient, rural pond water quality was significantly different from both peri-urban and urban ponds, whereas peri-urban and urban pond water quality was not significantly different. Temporally, differences between water quality parameters and sampling periods indicate that surrounding land use, land cover, and precipitation influence parameter concentrations across the urbanization gradient. Information from this study is useful to water professionals dealing with urban development and sprawl that continue to impact water and natural habitat.

Effects of biochar application and irrigation rate on the soil phosphorus leaching risk of fluvisol profiles in open vegetable fields

Biochar application was reported to influence soil phosphorus (P) leaching, but the reports are conflicting, and could be related to soil depth and water management. A field trial of a Wild Cabbage-Chinese Cabbage rotation was used to investigate the effect of biochar application and irrigation volume on P leaching risk in fluvisol soil profiles (0-20 cm, 20-50 cm, 50-100 cm) in the Chaobai River basin. The experiment included two biochar levels [0 (-BC), 30 t/hm2 (+BC)], and two irrigation levels [conventional irrigation (CI) and water-saving irrigation (WSI)]. The irrigation rate of WSI was 80% of CI. The results demonstrated that there was no significant difference in soil leachable P in the soil profiles under the two irrigation volumes, while biochar application tended to increase soil leachable P in the top layer soil (0-20 cm) and subsurface layer soil (20-50 cm) irrespective of the irrigation rate. The average value of the P leaching "change point" in the soil profiles with +BC was significantly higher than that with -BC (0-20 cm: 35.52 mg kg-1 vs. 25.86 mg kg-1; 20-50 cm: 27.61 mg kg-1 vs. 20.02 mg kg-1). Additionally, the P leaching risk was observed in all top layer soil (0-20 cm) irrespective of irrigation rate and biochar application, and the P leaching risk in the subsurface layer (20-50 cm) with +BC was lower than that with -BC, especially under WSI. Therefore, it is recommended that biochar application combined with water-saving irrigation could be used as a measure for controlling soil phosphorus leaching under open field vegetable rotation in the alluvial soil of Chaobai River basin.

Numerical Representation of Groundwater-Surface Water Exchange and the Effect on Streamflow Contribution Estimates

The effects of streams and drainage representation in 3D numerical catchment scale models on estimated streamflow contribution were investigated. MODFLOW-USG was used to represent complex geology and a stream network with two different conceptualizations—one with equal cell discretization in the entire model domain and another with refined cell discretization along stream reaches. Both models were calibrated against a large data set including hydraulic heads and streamflow measurements. Though the optimized hydraulic parameters and statistical performance of both model conceptualizations were comparable, their estimated streamflow contribution differed substantially. In the conceptualization with equal cell discretization, the drainage contribution to the streamflow was 13% compared to 41% in the conceptualization with refined cell discretization. The increase in drainage contribution to streamflow was attributed to the increase in drainage area in proximity to the stream reaches arising from the refined discretization. e.g., the cell refinement along stream reaches reduced the area occupied by stream cells allowing for increased drain area adjacent to the stream reaches. As such, an increase in drainage area equivalent to 7% yielded a 146% increase in drainage contribution to streamflow. In-stream field measurements of groundwater-surface water exchange fluxes that were qualitatively compared to calculated fluxes from the models indicated that estimates from the refined model discretization were more representative. Hence, the results of this study accentuate the importance of being able to represent stream and drain flow contribution correctly, that is, to achieve representative exchange fluxes that are crucial in simulating groundwater–surface water exchange of both flow and solute transport in catchment scale modeling. To that end, the in-stream measurements of exchange fluxes showed the potential to serve as a proxy to numerically estimate drainage contribution that is not readily available at the catchment scale.

A generic approach to evaluate costs and effectiveness of agricultural Beneficial Management Practices to improve water quality management

Nutrient export from agricultural areas is among the main contributors to water pollution in various watersheds. Agricultural Beneficial Management Practices (BMPs) are commonly used to reduce excessive nutrient runoff and improve water quality. The successful uptake of BMPs not only depends on their effectiveness but also on their costs of implementation. This study conducts a set of cost-effectiveness analyses to help stakeholders identify their preferred combinations of BMPs in the Qu'Appelle River Basin, a typical watershed in the Canadian Prairies. The considered BMPs are related to cattle and cropping farms and are initially selected by agricultural producers in this region. The analyses use a water quality model to estimate the impact of implementing BMPs on nutrient export, and the cost estimation model to approximate the cost of implementing BMPs at tributary and watershed scales. Our results show that BMPs' effectiveness, total costs of implementation and costs per kilogram of nutrient abatement vary between tributaries. However, wetland conservation is among the optimal practices to improve water quality across the watershed. It is also found that the rates of BMP adoption by stakeholders can influence the effectiveness of practices in a large watershed scale, which highlights the importance of stakeholder engagement in water quality management. This type of analyses can help stakeholders choose single or a combination of BMPs according to their available budget and acceptable levels of reduction in nutrients.

Characterizing bioavailable phosphorus concentrations in an agricultural stream during hydrologic and streambed disturbances

In freshwater ecosystems, phosphorus (P) is often considered a growth-limiting nutrient. The use of fertilizers on agricultural fields has led to runoff-driven increases in P availability in streams, and the subsequent eutrophication of downstream ecosystems. Isolated storms and periodic streambed dredging are examples of two common disturbances that contribute dissolved and particulate P to agricultural streams, which can be quantified as soluble reactive P (SRP) using the molybdate-blue method on filtered water samples, or total P (TP) measured using digestions on unfiltered water reflecting all forms of P. While SRP is often considered an approximation of bioavailable P (BAP), research has shown that this is not always the case. Current methods used to estimate BAP do not account for the role of biology (e.g., NaOH extractions) or require specialized platforms (e.g., algal bioassays). Here, in addition to routine analysis of SRP and TP, we used a novel yeast-based bioassay with unfiltered sample water to estimate BAP concentrations during two storms (top 80% and > 95% flow quantiles), and downstream of a reach where management-associated dredging disturbed the streambed. We found that the BAP concentrations were often greater than SRP, suggesting that SRP is not fully representative of P bioavailability. The SRP concentrations were similarly elevated during the two storms, but remained consistently low during streambed disturbance. In contrast, turbidity and TP were elevated during all events. The BAP concentrations were significantly related to turbidity during all disturbance events, but with TP only during storms. The novel yeast assay suggests that BAP export can exceed SRP, particularly when streams are not in equilibrium, such as the rising limb of storms or during active dredging.

One size does not fit all: Toward regional conservation practice guidance to reduce phosphorus loss risk in the Lake Erie watershed

Agricultural phosphorus (P) losses to surface water bodies remain a global eutrophication concern, despite the application of conservation practices on farm fields. Although it is generally agreed upon that the use of multiple conservation practices ("stacking") will lead to greater improvements to water quality, this may not be cost effective to farmers, reducing the likelihood of adoption. At present, wholesale recommendations of conservation practices are given; however, the application of specific conservation practices in certain environments (e.g., no-till with surface application, cover crops) may not be effective and can even lead to unintended consequences. In this paper, we present the Lake Erie watershed as a case study. The Lake Erie watershed contains regions with unique physical geographies that include differences in climate, soil, topography, and land use, which have implications for both P transport from agricultural fields and the efficacy of conservation practices in mitigating P losses. We define major regions within the Lake Erie watershed where common strategies for conservation practice implementation are appropriate, and we propose a five-step plan for bringing regionally tailored, adaptive, and cost-conscious conservation practice into watershed planning. Although this paper is specific to the Lake Erie watershed, our framework can be transferred across broader geographic regions to provide guidance for watershed planning.

Mapping of Agricultural Subsurface Drainage Systems Using Unmanned Aerial Vehicle Imagery and Ground Penetrating Radar

Agricultural subsurface drainage systems are commonly installed on farmland to remove the excess water from poorly drained soils. Conventional methods for drainage mapping such as tile probes and trenching equipment are laborious, cause pipe damage, and are often inefficient to apply at large spatial scales. Knowledge of locations of an existing drainage network is crucial to understand the increased leaching and offsite release of drainage discharge and to retrofit the new drain lines within the existing drainage system. Recent technological developments in non-destructive techniques might provide a potential alternative solution. The objective of this study was to determine the suitability of unmanned aerial vehicle (UAV) imagery collected using three different cameras (visible-color, multispectral, and thermal infrared) and ground penetrating radar (GPR) for subsurface drainage mapping. Both the techniques are complementary in terms of their usage, applicability, and the properties they measure and were applied at four different sites in the Midwest USA. At Site-1, both the UAV imagery and GPR were equally successful across the entire field, while at Site-2, the UAV imagery was successful in one section of the field, and GPR proved to be useful in the other section where the UAV imagery failed to capture the drainage pipes’ location. At Site-3, less to no success was observed in finding the drain lines using UAV imagery captured on bare ground conditions, whereas good success was achieved using GPR. Conversely, at Site-4, the UAV imagery was successful and GPR failed to capture the drainage pipes’ location. Although UAV imagery seems to be an attractive solution for mapping agricultural subsurface drainage systems as it is cost-effective and can cover large field areas, the results suggest the usefulness of GPR to complement the former as both a mapping and validation technique. Hence, this case study compares and contrasts the suitability of both the methods, provides guidance on the optimal survey timing, and recommends their combined usage given both the technologies are available to deploy for drainage mapping purposes.

Strong adsorption of phosphate by amorphous lanthanum carbonate nano-adsorbents

Phosphorus removal is a crucial aspect of controlling water pollution and eutrophication. In this study, the preparation of lanthanum carbonate (LC) nano-adsorbents for the efficient removal of phosphate (P) from water and wastewater was investigated. Results from XRD, SEM and Zeta potential analyses revealed that addition of magnesium ions and adjustment of the reaction temperature could control the morphology and microstructure of LC. Effects of initial pH, adsorbent dosage, contact time, and the water matrix on P adsorption were investigated. Batch adsorption experiments revealed that LC showed strong performance on P removal over a wide pH range (3.0 to 11.0). The kinetic data followed a pseudo-second-order model, and equilibrium data were well fitted by the Langmuir model with a maximum adsorption capacity of 112.9 mg P/g. Adsorption thermodynamics showed that the adsorption process was exothermic and spontaneous. Results of a monolayer model for single adsorption indicated that P could completely interact with two or more functional groups from the LC surface. In the presence of competing ions (F^-, Cl^-, SO₄^2-, NO₃^-, and HCO₃^-), LC maintained high selectivity for phosphate. For a real effluent, the P concentration was efficiently reduced from 3.2 mg P/L to below 0.5 mg P/L at a dose of 0.5 g/L LC. All the results suggested that LC can serve as a promising adsorbent for P removal in a wide range of pH, and thus could meet the stricter discharge regulations from actual wastewater.

Phosphorus Transport along the Cropland–Riparian–Stream Continuum in Cold Climate Agroecosystems: A Review

Phosphorus (P) loss from cropland to ground and surface waters is a global concern. In cold climates (CCs), freeze–thaw cycles, snowmelt runoff events, and seasonally wet soils increase P loss potential while limiting P removal effectiveness of riparian buffer zones (RBZs) and other practices. While RBZs can help reduce particulate P transfer to streams, attenuation of dissolved P forms is more challenging. Moreover, P transport studies often focus on either cropland or RBZs exclusively rather than spanning the natural cropland–RBZ–stream gradient, defined here as the cropland–RBZ–stream continuum. Watershed P transport models and agronomic P site indices are commonly used to identify critical source areas; however, RBZ effects on P transport are usually not included. In addition, the coarse resolution of watershed P models may not capture finer-scale soil factors affecting P mobilization. It is clear that site microtopography and hydrology are closely linked and important drivers of P release and transport in overland flow. Combining light detection and ranging (LiDAR) based digital elevation models with P site indices and process-based models show promise for mapping and modeling P transport risk in cropland-RBZ areas; however, a better mechanistic understanding of processes controlling mobile P species across regions is needed. Broader predictive approaches integrating soil hydro-biogeochemical processes with real-time hydroclimatic data and risk assessment tools also hold promise for improving P transport risk assessment in CCs.

Critical Review of Polyphosphate and Polyphosphate Accumulating Organisms for Agricultural Water Quality Management

Despite ongoing management efforts, phosphorus (P) loading from agricultural landscapes continues to impair water quality. Wastewater treatment research has enhanced our knowledge of microbial mechanisms influencing P cycling, especially regarding microbes known as polyphosphate accumulating organisms (PAOs) that store P as polyphosphate (polyP) under oxic conditions and release P under anoxic conditions. However, there is limited application of PAO research to reduce agricultural P loading and improve water quality. Herein, we conducted a meta-analysis to identify articles in Web of Science on polyP and its use by PAOs across five disciplines (i.e., wastewater treatment, terrestrial, freshwater, marine, and agriculture). We also summarized research that provides preliminary support for PAO-mediated P cycling in natural habitats. Terrestrial, freshwater, marine, and agriculture disciplines had fewer polyP and PAO articles compared to wastewater treatment, with agriculture consistently having the least. Most meta-analysis articles did not overlap disciplines. We found preliminary support for PAOs in natural habitats and identified several knowledge gaps and research opportunities. There is an urgent need for interdisciplinary research linking PAOs, polyP, and oxygen availability with existing knowledge of P forms and cycling mechanisms in natural and agricultural environments to improve agricultural P management strategies and achieve water quality goals.

Cyanobacterial bloom phenology in Saginaw Bay from MODIS and a comparative look with western Lake Erie

Saginaw Bay and western Lake Erie basin (WLEB) are eutrophic catchments in the Laurentian Great Lakes that experience annual, summer-time cyanobacterial blooms. Both basins share many features including similar size, shallow depths, and equivalent-sized watersheds. They are geographically close and both basins derive a preponderance of their nutrient supply from a single river. Despite these similarities, the bloom phenology in each basin is quite different. The blooms in Saginaw Bay occur at the same time and place and at the same moderate severity level each year. The WLEB, in contrast, exhibits far greater interannual variability in the timing, location, and severity of the bloom than Saginaw Bay, consistent with greater and more variable phosphorus inputs. Saginaw Bay has bloom biomass that corresponds to relatively mild blooms in WLEB, and also has equivalent phosphorus loads. This result suggests that if inputs of P into the WLEB were reduced to similarly sized loads as Saginaw Bay the most severe blooms would be abated. Above 500 t P input, which occur in WLEB, blooms increase non-linearly indicating any reduction in P-input at the highest inputs levels currently occurring in the WLEB, would yield disproportionately large reductions in cyanobacterial bloom intensity. As the maximum phosphorus loads in Saginaw Bay lie just below this inflection point, shifts in the Saginaw Bay watershed toward greater agriculture uses and less wetlands may substantially increase the risk of more intense cyanobacterial blooms than presently occur.

Using nutrient transport data to characterize and identify the presence of surface inlets in regions with subsurface drainage

Surface inlets route ponded surface water into subsurface drainage networks and are prevalent throughout North America. Despite serving as a nutrient loss pathway, contributing to downstream water quality degradation, surface inlets are thought to be underreported in drainage studies within the literature. Previous studies have demonstrated the footprint that surface inlets have on nutrient transport and drainage effluent but are site specific and focused on individual events. Moreover, although their ubiquitous presence is assumed, no regional surface inlet database exists. To this end, a structured review was undertaken with two goals. First, the MANAGE Drain Load database, consisting of nearly 1,500 site-years of drainage and nutrient data, was analyzed to determine distinctions between areas with and without surface inlets. The median annual total phosphorus (TP) load was greater at site-years with surface inlets (0.40 kg ha^-1 ) than site-years without (0.21 kg ha^-1 ). The opposite emerged for dissolved nitrogen (DN) loads as site-years with surface inlet had a smaller median annual load (3.3 kg ha^-1 ) than site-years without (23.0 kg ha^-1 ). This relationship is attributed to immobile TP being transported primarily through overland flow and routed to subsurface drains via surface inlets and to relatively more mobile DN being subsurface driven, bypassed in settings with surface inlets. No statistical differences were found in annual drainage or ratios of particulate P to TP between site-years with and without surface inlets. Second, a logistic regression model was developed that predicts the presence of surface inlets within MANAGE. Eighteen percent of site-years and 21% of sites were predicted to have surface inlets.