Repulsive-attractive models for the impact of two predators on prey species varying in anti-predator response

Abstract

This study considers the dynamical interaction of two predatory carnivores (Lions (Panthera leo) and Spotted Hyaenas (Crocuta crocuta)) and three of their common prey (Buffalo (Syncerus caffer), Warthog (Phacochoerus africanus) and Kudu (Tragelaphus strepsiceros)). The dependence on spatial structure of species’ interaction stimulated the author to formulate reaction-diffusion models to explain the dynamics of predator-prey relationships in ecology. These models were used to predict and explain the effect of threshold populations, predator additional food and prey refuge on the general species’ dynamics. Vital parameters that model additional food to predators, prey refuge and population thresholds were given due attention in the analyses. The stability of a predator-prey model for an ecosystem faced with a prey out-flux which is analogous to and modelled as an Allee effect was investigated. The results highlight the bounds for the conversion efficiency of prey biomass to predator biomass (fertility gain) for which stability of the three species ecosystem model can be attained. Global stability analysis results showed that the prey (warthog) population density should exceed the sum of its carrying capacity and threshold value minus its equilibrium value i.e., W >(Kw + $) −W . This result shows that the warthog’s equilibrium population density is bounded above by population thresholds, i.e., W < (Kw+$). Besides showing the occurrence under parameter space of the so-called paradox of enrichment, early indicators of chaos can also be deduced. In addition, numerical results revealed stable oscillatory behaviour and stable spirals of the species as predator fertility rate, mortality rate and prey threshold were varied. The stabilising effect of prey refuge due to variations in predator fertility and proportion of prey in the refuge was studied. Formulation and analysis of a robust mathematical model for two predators having an overlapping dietary niche were also done. The Beddington-DeAngelis functional and numerical responses which are relevant in addressing the Principle of Competitive Exclusion as species interact were incorporated in the model. The stabilizing effect of additional food in relation to the relative diffusivity D, and wave number k, was investigated. Stability, dissipativity, permanence, persistence and periodicity of the model were studied using the routine and limit cycle perturbation methods. The periodic solutions (b 1 and b 3), which influence the dispersal rate (') of the interacting species, have been shown to be controlled by the wave number. For stability, and in order to overcome predator natural mortality, the nutritional value of predator additional food has been shown to be of high quality that can enhance predator fertility gain. The threshold relationships between various ecosystem parameters and the carrying capacity of the game park for the prey species were also deduced to ensure ecosystem persistence. Besides revealing irregular periodic travelling wave behaviour due to predator interference, numerical results also show oscillatory temporal dynamics resulting from additional food supplements combined with high predation rates