Turing instability in pioneer/climax species interactions

Systems of pioneer and climax species are used to model interactions of species whose reproductive capacity is sensitive to population density in their shared ecosystem. Intraspecies interaction coefficients can be adjusted so that spatially homogeneous solutions are stable to small perturbations. In a reaction-diffusion pioneer/climax model we will determine the critical value of the diffusion rate of the climax species, below which the equilibrium solution is unstable to non-homogeneous perturbations. For diffusion rates smaller than this critical value, an equilibrium solution remains stable to spatially homogeneous perturbations but is unstable to non-homogeneous perturbations. A Turing (diffusional) bifurcation leads to the formation of spatial patterns in species' densities. Forcing, interpreted as stocking or harvesting of the species, can reverse the bifurcation and establish equilibrium solutions which are stable to small perturbations. The implicit function theorem is used to determine whether stocking or harvesting of one of the species in the model is the appropriate remedy for diffusional instability. The use of stocking or harvesting by a natural resource manager thus influences the long-term dynamics and spatial distribution of species in a pioneer/climax ecosystem.