Exploring dynamical complexity in diffusion driven predator–prey systems: Effect of toxin producing phytoplankton and spatial heterogeneities

In this paper, dynamical complexities in two reaction–diffusion (RD) model systems are explored. A spatial heterogeneity in the form of linear spatial gradient in the reproductive growth rate of the phytoplankton is incorporated in both the model systems. Extra mortality of the zooplankton due to toxin production by the phytoplankton is included in the second reaction diffusion model system. Effect of toxin production and spatial heterogeneity in the model systems are studied. Toxin production does not seem to have an appreciable effect on the asymptotic dynamics of the model systems. On the other hand, spatial heterogeneity does influence the dynamics. In particular, it increases the frequency of occurrence of chaos as evident from two dimensional parameter scans. Both these model systems display short term recurrent chaos [Rai V. Chaos in natural populations: edge or wedge? Ecol Complex 2004;1: 127–38] as they reside on ‘edges of chaos’ (EOC) [Rai V, Upadhyay RK. Evolving to the edge of chaos: chance or necessity? Chaos, Solitons and Fractals 2006;30:1074–87]. This suggests that the ecological systems have a tendency to evolve to EOC. The study corroborates the inferences drawn from an earlier study by Rai and Upadhyay [Rai V, Upadhyay RK. Evolving to the edge of chaos: chance or necessity? Chaos, Solitons and Fractals 2006;30:1074–87]. The system’s dynamics is largely unpredictable and admits bursts of short-term predictability.