Sediment phosphorus buffering in streams at baseflow: A meta-analysis

Multi-decadal impacts of effluent loading on phosphorus sorption capacity in a restored wetland

Natural wetlands are widely used and cost-effective systems for the passive remediation of phosphorus (P)-rich surface waters from various effluent sources. Yet the long-term biogeochemical impacts of effluent loading on wetland P retention capacity is unclear. Here, we had a unique opportunity to document the spatio-temporal evolution of sediment P sorption over a ∼25-year period of constant municipal and industrial effluent loading, as part of a wetland restoration and wastewater treatment strategy in one of the largest restored wetlands in Canada. Sediment P sorption experiments across Frank Lake's three basins revealed a wide spatial variation in sorption capacity, closely linked to sediment geochemistry gradients (Ca, Fe, and Mn). Relative to a similar study ∼25 years prior, P sorption capacity has become exhausted near the effluent inlet, but remarkably, remains elevated throughout the rest of the wetland. Compared to other prairie wetlands and global aquatic ecosystems, Frank Lake has a greater capacity overall to retain P through sediment sorption. Given the paucity of long-term (multi-decade) data on wetland response to effluent loading, we provide key insights into the dynamics of wetland P cycling in human-dominated watersheds.

Linking sediment geochemistry with catchment processes, internal phosphorus loading and lake water quality

Research in the field of sediment geochemistry suggests potential linkages between catchment processes (land use), internal phosphorus (P) loading and lake water quality, but evidence is still poorly quantified due to a limited amount of data. Here we address the issues based on a comprehensive data set from 27 lakes in southern Finland. Specifically, we aimed at: 1) elucidating factors behind spatial variations in sediment geochemistry; 2) assessing the impact of diagenetic transformation on sediment P regeneration across lakes based on the changes in the vertical distribution of sediment components; 3) exploring the role of the sediment P forms in internal P loading (IL), and 4) determining the impact of IL on lake water quality. The relationship between sediment P concentration and field area percentage (FA%) was statistically significant in (mainly eutrophic) lakes with catchments that included more than 10 % of fields. We found that sediment iron-bound P (Fe-P) increased with increasing FA%, which agrees with the high expected losses from the cultivated areas. Additionally, populated areas increased the pool of sediment Fe-P. Internal P loading was significantly positively related to both sediment Fe-P and sediment organic P (Org-P). However, Org-P was not significant (as the third predictor) in models that had a trophic state variable as the first predictor and Fe-P as the second predictor. Further, the vertical profiles of sediment components indicated a role of diagenetic transformations in the long-term sediment P release, especially in lakes with deeper maximum depth and longer water residence time. Finally, IL was significantly positively correlated to water quality variables including phytoplankton biomass, its proportion of cyanobacteria, chlorophyll a concentration and trophic state index. Our findings suggest that reduction of P losses from the field and populated areas will decrease internal P loads and increase water quality through a reduced pool of Fe-P.

Phosphorus lability across diverse agricultural contexts with legacy sources

The buffering of phosphorus (P) in the landscape delays management outcomes for water quality. If stored in labile form (readily exchangeable and bioavailable), P may readily pollute waters. We studied labile P and its intensity for >600 soils and sediments across seven study locations in the United States. Stocks of labile P were large enough to sustain high P losses for decades, indicating the transport-limited regime typical of legacy P. Sediments were commonly more P-sorptive than nearby soils. Soils in the top 5 cm had 1.3-3.0 times more labile P than soils at 5-15 cm. Stratification in soil test P and total P was, however, less consistent. As P exchange via sorption processes follows the difference in intensities between soil/sediment surface and solution, we built a model for the equilibrium phosphate concentration at net zero sorption (EPC0) as a function of labile P (quantity) and buffer capacity. Despite widely varying properties across sites, the model generalized well for all soils and sediments: EPC0 increased sharply with more labile P and to greater degree when buffer capacity was low or sorption sites were likely more saturated. This quantity-intensity-capacity relationship is central to the P transport models we rely on today. Our data inform the improvement of such P models, which will be necessary to predict the impacts of legacy P. Further, this work reaffirms the position of labile P as a key focus for environmental P management-a view Dr. Sharpley developed in the 1980s with fewer data and resources.

Trade-offs between nitrogen and phosphorus removal with floodplain remediation in agricultural streams

To improve water quality and reduce instream erosion, floodplain remediation along agricultural streams can provide multiple ecosystem services through biogeochemical and fluvial processes. During floodplain inundation, longer water residence time and periodic anoxic conditions can lead to increased nitrogen (N) removal through denitrification but also mobilization of phosphorus (P), impeding overall water quality improvements. To investigate the capacity for N and P processing in remediated streams, we measured potential denitrification and nitrous oxide production and yields together with potential P desorption and P fractions in floodplain and stream sediments in ten catchments in Sweden. Sediment P desorption was measured as equilibrium P concentration, using P isotherm incubations. Denitrification rates were measured with the acetylene inhibition method. Sediment nutrient process rates were combined with hydrochemical monitoring along remediated streams and their paired upstream control reaches of trapezoidal shape to determine the impact of floodplains on water quality. The correlation between floodplain denitrification rates and P desorption (r = 0.53, p = 0.02) revealed a trade-off between soluble reactive P (SRP) and nitrate removal, driven by stream water connectivity to floodplains. Nitrous oxide production was not affected by differences in P processing, but nitrous oxide yields decreased with higher denitrification and P desorption. The release of SRP from floodplains (0.03 ± 0.41 mg P kg-1 day-1) was significantly lower than from trapezoidal stream banks (0.38 ± 0.37 mg P kg-1 day-1), predicted by long-term SRP concentrations in stream water and floodplain inundation frequency. The overall impact of SRP release from floodplains on stream SRP concentrations in remediated reaches was limited. However, the remediated reaches showing increased stream SRP concentrations were also frequently inundated and had higher labile P content and coarse soil texture in floodplain sediments. To fully realize the potential for water quality improvements with constructed floodplains in agricultural streams, the promotion of denitrification through increased inundation should be balanced against the risk of P release from sediments, particularly in streams with high SRP inputs.

Evaluating the longevity of in-stream phosphorus legacies: A downstream cascade of recovery following point source remediation

In-stream phosphorus (P) legacies cause lags between upstream remediation and downstream load reductions. However, the length of these lags is largely unknown, especially for long stream distances. As a result, lag time estimates at the large-watershed scale have been abstract and sometimes understated. Here, we leverage a large area watershed model with newly improved in-stream P simulation (SWAT+P.R&R) to evaluate the magnitude, longevity, and spatial cascade of legacy P remobilization in a U.S. corn belt watershed. Our results illustrate the "spiraling recovery" of P loads after a hypothetical point source remediation, where locations further downstream take longer to recover to baseline load levels. At the watershed outlet, in-stream legacy P contributions are equivalent to 30% of the baseline average annual P loads for three years after remediation. In-stream legacies do not approach exhaustion (95% remobilized) until at least 9 years after remediation. In hypothetical weather scenarios beginning with dry years, legacy contributions persist even longer. These findings (1) suggest that in-stream legacies could impact P loads for years to decades in large river basins, (2) support explicit accounting for spatial scale in future studies of in-stream legacies, and (3) provide concerning implications for water quality recovery in large river basins.

Biogeochemical Processes and Microbial Dynamics Governing Phosphorus Retention and Release in Sediments: A Case Study in Lower Great Lakes Headwaters

The ability of headwater bed and suspended sediments to mitigate non-point agricultural phosphorus (P) loads to the lower Great Lakes is recognized, but the specific biogeochemical processes promoting sediment P retention or internal P release remain poorly understood. To elucidate these mechanisms, three headwater segments located within priority watersheds of Southern Ontario, Canada, were sampled through the growing season of 2018-2020. The study employed equilibrium P assays along with novel assessments of legacy watershed nutrients, nitrogen (N) concentrations, sediment redox, and microbial community composition. 20-year data revealed elevated total P (TP) and total Nitrogen (TN) at an inorganic fertilizer and manure fertilizer-impacted site, respectively. Overall, sampled sites acted as P sinks; however, agricultural sediments exhibited significantly lower buffering capacity compared to a reference forested watershed. Collection of fine suspended sediment (<63 µm) through time-integrated sampling showed the suspended load at the inorganic-fertilized site was saturated with P, indicating a greater potential for P release into surface waters compared to bed sediments. Through vertical microsensor profiling and DNA sequencing of the sediment microbial community, site-specific factors associated with a distinct P-source event were identified. These included rapid depletion of dissolved oxygen (DO) across the sediment water interface (SWI), as well as the presence of nitrate-reducing bacterial and ammonia-oxidizing archaeal (AOA) genera. This research provides valuable insights into the dynamics of P in headwaters, shedding light on P retention and release. Understanding these processes is crucial for effective management strategies aimed at mitigating P pollution to the lower Great Lakes.

The spatial and temporal dynamics of sediment phosphorus attenuation and release in impacted stream catchments

Sediments can attenuate phosphorus (P) from overlying water and reduce trophic status in zero and first order ditches and streams. These features can be considered as intermediate mitigation features between P mobilised from land, and onward delivery to river systems, if the risk of chemical P release from sediments is minimal. However, risk assessments are rarely based on temporal scale dynamics and especially at fine scale in both sediment and water column environments. In this study, in eutrophic stream catchments, bed sediments were tested fortnightly and spatially over one year for EPC0 (to derive phosphate exchange potential-PEP) and for P across a spectrum from labile to recalcitrant fractions. At the same time stream discharge and P concentrations were measured synchronously at high frequency and resolved to 1-hour intervals and indicated high water quality pressures at all flow rates. PEP indicated spatial and temporal changes most likely caused by periods of source disconnection/reconnection and sediment mobilisation during storm events, moving from periods of high attenuation potential to near saturation. Despite these spatial and temporal changes, PEP did not indicate much potential for chemical P release from the sediments (distributing mostly below or close to zero). However, this may be a misleading risk assessment by itself as physical P release, especially of the labile bicarbonate-dithionite (B-D) P fraction of sediments, was a more dominant process mobilised during storm events reducing by up to 84 % during a succession of summer storm events. The total P and total reactive P loads monitored leaving the catchments were coincident with these changes. The specific downstream trophic effects of this episodic P release will need to be assessed in terms of its bioavailability, in combination with other more noted diffuse and point P source processes.

The influence of turbulence on sediment phosphorus sorption

The effect of mean flow velocity on phosphorus (P) partitioning between water and sediment has received much attention in recent decades. However, the impact of turbulence on the efficiency and capability of sediment adsorbing and desorbing dissolved inorganic phosphorus (DIP) is still unclear. A series of contrasting experiments on the sediment sorption and desorption of DIP with the flow turbulence kinetic energy (TKE) ranging from 1.95 to 2.93 pa have been conducted. It was found that the adsorbed P onto unit mass of sediment increases with the increase in TKE. It is because an increase in TKE results in a rise in the effective adsorption capacity of sediment (b_m) by 20-30% during the adsorption process. The b_m shows the maximum rise from 0.18 to 0.25 mg/g when TKE increases from 1.95 to 2.93 pa with a fixed sediment concentration of 0.5 g/L. To account for the direct effect of TKE on P adsorption, the Langmuir model is modified by introducing a newly defined coefficient (f_A-TKE). The f_A-TKE shows a good linear relationship with TKE. Comparison between the modified model and the classic model shows that the amount of adsorbed P could be overestimated by over 50% if the direct effect of turbulence intensity is ignored. The experimental data show that the increase in TKE also enhances the desorption process, with the degree of P desorption (D_des) increased by 44%. The relation between D_des and TKE can be well represented using a logarithmic function to quantify the direct effect of turbulence intensity on desorption of P.

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.

Periphyton Phosphorus Uptake in Response to Dynamic Concentrations in Streams: Assimilation and Changes to Intracellular Speciation

Effective modeling and management of phosphorus (P) losses from landscapes to receiving waterbodies requires an adequate understanding of P retention and remobilization along the terrestrial-aquatic continuum. Within aquatic ecosystems, the stream periphyton can transiently store bioavailable P through uptake and incorporation into biomass during subscouring and baseflow conditions. However, the capacity of stream periphyton to respond to dynamic P concentrations, which are ubiquitous in streams, is largely unknown. Our study used artificial streams to impose short periods (48 h) of high SRP concentration on stream periphyton acclimated to P scarcity. We examined periphyton P content and speciation through nuclear magnetic resonance spectroscopy to elucidate the intracellular storage and transformation of P taken up across a gradient of transiently elevated SRP availabilities. Our study demonstrates that the stream periphyton not only takes up significant quantities of P following a 48-h high P pulse but also sustains supplemental growth over extended periods of time (10 days), following the reestablishment of P scarcity by efficiently assimilating P stored as polyphosphates into functional biomass (i.e., phospho-monoesters and phospho-diesters). Although P uptake and intracellular storage approached an upper limit across the experimentally imposed SRP pulse gradient, our findings demonstrate the previously underappreciated extent to which the periphyton can modulate the timing and magnitude of P delivery from streams. Further elucidating these intricacies in the transient storage potential of periphyton highlights opportunities to enhance the predictive capacity of watershed nutrient models and potentially improve watershed P management.

Watershed- and reach-scale drivers of phosphorus retention and release by streambed sediment in a western Lake Erie watershed during summer

Reducing phosphorus (P) concentrations in aquatic ecosystems, is necessary to improve water quality and reduce the occurrence of harmful cyanobacterial algal blooms. Managing P reduction requires information on the role rivers play in P transport from land to downstream water bodies, but we have a poor understanding of when and where river systems are P sources or sinks. During the summers of 2019 and 2021, we sampled streambed sediment at 78 sites throughout the Maumee River network (a major source of P loads to Lake Erie) focusing on the zero equilibrium P concentration (EPC₀), the soluble reactive phosphorus (SRP) concentration at which sediment neither sorbs nor desorbs P. We used structural equation modeling to identify direct and indirect drivers of EPC₀. Stream sediment was a P sink at 40 % and 67 % of sites in 2019 and 2021, respectively. During both years, spatial variation in EPC₀ was shaped by stream water SRP concentrations, sediment P saturation, and sediment physicochemical characteristics. In turn, SRP concentrations and sediment P saturation (PSR) were influenced by agricultural land use and stream size. Effect of stream size differed among years with stream size having a greater effect on SRP in 2019 and on PSR in 2021. Streambed sediment is currently a net P sink across the sites sampled in the Maumee River network during summer, but sediment at these locations, especially sites in headwater streams, may become a P source if stream water SRP concentrations decrease. Our results improve the understanding of watershed- and reach-scale controls on EPC₀ but also indicate the need for further research on how changes in SRP concentration as a result of conservation management implementation influences the role of streambed sediment in P transport to Lake Erie.

River phosphorus cycling during high flow may constrain Lake Erie cyanobacteria blooms

Cyanobacterial harmful blooms have been increasing worldwide, due in part to excessive phosphorus (P) losses from agriculture-dominated watersheds. Unfortunately, cyanobacteria bloom management is often complicated by uncertainty associated with river P cycling. River P cycling mediates P exports during low flow but has been assumed to be unimportant during high flows. Thus, we examined interactions between dissolved reactive phosphorus (DRP) and suspended sediment P during high flows in the Maumee River network, focusing on March-June Maumee River DRP exports, which fuel recurring cyanobacteria blooms in Lake Erie. We estimate that during 2003-2019 March to June high flow events, P sorption reduced DRP exports by an average of 13-27%, depending upon the colloidal-P:DRP ratio, decreasing the bioavailability of P exports, and potentially constraining cyanobacteria blooms by 13-40%. Phosphorus sorption was likely lower during 2003-2019 than 1975-2002 due to reductions in suspended sediment loads, associated with soil-erosion-minimizing agricultural practices. This unintended outcome of erosion management has likely decreased P sorption, increased DRP exports to Lake Erie, and subsequent cyanobacteria blooms. In other watersheds, DRP-sediment P interactions during high flow could have a positive or negative effect on DRP exports; therefore, P management should consider riverine P cycles, particularly during high flow events, to avoid undermining expensive P mitigation efforts.

Applying the nutrient transfer continuum framework to phosphorus and nitrogen losses from livestock farmyards to watercourses

Farmyards are commonly conceptualized as point sources of nutrient pollution nested within the wider agricultural landscape. However, within farmyards there are individual sources and delivery pathways, each of which is affected by a range of management practices and infrastructure. Rainfall mobilizes these nutrients, which may then be delivered to a receptor or to the wider drainage network. As such, the nutrient transfer continuum (NTC), which has been established as a framework to understand and mitigate nutrient loss at a landscape scale, can be similarly applied to disentangle the stages of nutrient transfer from farmyards. The NTC differentiates nutrient transfer into source, mobilization, delivery, and impact stages. This differentiation allows targeting of mitigation measures and evaluation of costs and benefits. This review paper applies the NTC template to farmyard nitrogen and phosphorus transport to conceptualize causative factors and to identify mitigation options.