A model is developed to study the phosphorus dynamics in aquatic sediments and to conduct dynamic predictions of phosphorus release across a sediment-water interface. The model focuses on the sediment active layer below the sediment-water interface and is based on primary mechanisms regulating phosphorus behavior in sediments, including effective diffusion, bioturbation mixing and burial processes (transport), organic decomposition, sorption kinetic processes and non-linear partitioning (mobilization). The effects of environmental factors such as dissolved oxygen and temperature are taken into account. The model is solved by numerical integration. The primary difference from models in the literature is that the model directly describes the dynamic behavior of dissolved, particulate exchangeable ortho-phosphorus and organic phosphorus in sediments, and incorporates dynamic sorption and non-linear partitioning processes. These improve model mechanisms and allow regulation of phosphorus flux through the sediment reservoir that acts as both a source and sink of phosphorus.