RESUMO
Nutrient accumulation in man-made reservoirs has been documented worldwide. Therefore, quantifying phosphorus loading and understanding its temporal dynamics in reservoirs is mandatory for sustainable water management. In this study, the Vollenweider's complete-mix phosphorus budget model was adapted to account for high water level variations, which are a common feature in tropical reservoirs, and for internal phosphorus loads. First- and zero-order kinetics were adopted to simulate phosphorus settling and release from the sediment layer, respectively, considering variable area of phosphorus release according to the height of the anoxic layer. The modeling approach was applied for a 52-months period to a 31-years-old reservoir located in the semiarid region of Brazil with 7.7 hm3 storage capacity. The simulations were supported by hydrological, meteorological and water quality data, as well as analyses of phosphorus partitioning of the reservoir bed sediment. The external phosphorus load was estimated from a relationship adjusted between inflow and phosphorus concentration, revealing an u-shaped pattern. Sedimentary phosphorus linked to iron and aluminum (PFeAl) increased over time and along the reservoir. Such measurements were used to estimate the internal phosphorus load, i.e., the yield from the bed sediments to the water column. The adaptations proposed to the model's structure improved its capacity to simulate phosphorus concentration in the water column, from "not satisfactory" to "good". We estimate that the internal phosphorus load currently accounts for 44% of the total load. It prevailed during the wet period, when reservoir stratification and hypolimnetic hypoxia were more notable, resulting in higher phosphorus concentration in the water column due to the combined effects of internal and external loadings. However, if the reservoir were 70 years older, the internal load would reach 83% of the total and the reservoir would become a source instead of a sink of phosphorus.