Dynamics of a Producer–Grazer Model Incorporating the Effects of Phosphorus Loading on Grazer’s Growth

Phosphorus is an essential element for all life forms, and it is also a limiting nutrient in many aquatic ecosystems. To keep track of the mismatch between the grazer’s phosphorus requirement and producer phosphorus content, stoichiometric models have been developed to explicitly incorporate food quality and food quantity. Most stoichiometric models have suggested that the grazer dynamics heavily depends on the producer phosphorus content when the producer has insufficient nutrient content [low phosphorus (P):carbon (C) ratio]. However, recent laboratory experiments have shown that the grazer dynamics are also affected by excess producer nutrient content (extremely high P:C ratio). This phenomenon is known as the “stoichiometric knife edge.” While the Peace et al. (Bull Math Biol 76(9):2175–2197, 2014) model has captured this phenomenon, it does not explicitly track P loading of the aquatic environment. Here, we extend the Peace et al. (2014) model by mechanistically deriving and tracking P loading in order to investigate the growth response of the grazer to the producer of varying P:C ratios. We analyze the dynamics of the system such as boundedness and positivity of the solutions, existence and stability conditions of boundary equilibria. Bifurcation diagram and simulations show that our model behaves qualitatively similar to the Peace et al. (2014) model. The model shows that the fate of the grazer population can be very sensitive to P loading. Furthermore, the structure of our model can easily be extended to incorporate seasonal P loading.