Mechanistically derived spatially heterogeneous producer-grazer model subject to stoichiometric constraints

Known stoichiometric models of a two species producer-grazer ecosystem have either neglected spatial dynamics or failed to track free phosphorus in the media. In this paper we present a spatially heterogeneous model that tracks phosphorus content in the producer and free phosphorus in the media. We simulate our model numerically under various environmental conditions. Multiple equilibria, with bistability and deterministic extinction of the grazer, are possible here. In conditions that had been previously studied without tracking free phosphorus we find cases where qualitatively different behavior is observed. In particular under certain environmental conditions previous models predict stable equilibria where our model predicts stable limit cycles near the surface. Oscillatory dynamics can have consequences on the population densities, which may spend some time at low values throughout the cycles where they are in danger of stochastic extinction.