Phosphorus retention and internal loading in the Bay of Quinte, Lake Ontario, using diagenetic modelling

Internal phosphorus (P) loading significantly contributes to hysteresis in ecosystem response to nutrient remediation, but the dynamics of sediment P transformations are often poorly characterized. Here, we applied a reaction-transport diagenetic model to investigate sediment P dynamics in the Bay of Quinte, a polymictic, spatially complex embayment of Lake Ontario, (Canada). We quantified spatial and temporal variability of sediment P binding forms and estimated P diffusive fluxes and sediment P retention in different parts of the bay. Our model supports the notion that diagenetic recycling of redox sensitive and organic bound P forms drive sediment P release. In the recent years, summer sediment P diffusive fluxes varied in the range of 3.2-3.6 mg P m^-2 d^-1 in the upper bay compared to 1.5 mg P m^-2 d^-1 in the middle-lower bay. Meanwhile sediment P retention ranged between 71% and 75% in the upper and middle-lower bay, respectively. The reconstruction of temporal trends of internal P loading in the past century, suggests that against the backdrop of reduced external P inputs, sediment P exerts growing control over the lake nutrient budget. Higher sediment P diffusive fluxes since mid-20th century with particular increase in the past 20 years in the shallower upper basins, emphasize limited sediment P retention potential and suggest prolonged ecosystem recovery, highlighting the importance of ongoing P control measures.