Search

International audience; Two predator-prey model formulations are studied: for the classical Rosenzweig-MacArthur (RM) model and the Mass Balance (MB) chemostat model. When the growth and loss rate of the predator is much smaller than that of the prey these models are slow-fast systems leading mathematically to singular perturbation problem. In contradiction to the RM-model, the resource for the prey are modelled explicitly in the MB-model but this comes with additional parameters. These parameter values are chosen such that the two models become easy to compare. In both models a transcritical bifurcation, a threshold above which invasion of predator into prey-only system occurs, and the Hopf bifurcation where the interior equilibrium becomes unstable leading to a stable limit cycle. The fast-slow limit cycles are called relaxation oscillations which for increasing differences in time scales leads to the well known degenerated trajectories being concatenations of slow parts of the trajectory and fast parts of the trajectory. In the fast-slow version of the RM-model a canard explosion of the stable limit cycles occurs in the oscillatory region of the parameter space. To our knowledge this type of dynamics has not been observed for the RM-model and not even for more complex ecosystem models. When a bifurcation parameter crosses the Hopf bifurcation point the amplitude of the emerging stable limit cycles increases. However, depending of the perturbation parameter the shape of this limit cycle changes abruptly from one consisting of two concatenated slow and fast episodes with small amplitude of the limit cycle, to a shape with large amplitude of which the shape is similar to the relaxation oscillation, the well known degenerated phase trajectories consisting of four episodes (concatenation of two slow and two fast). The canard explosion point is accurately predicted by using an extended asymptotic expansion technique in the perturbation and bifurcation parameter simultaneously where the small amplitude stable limit cycles exist. The predicted dynamics of the MB-model is in a large part of the parameter space similar to that of the RM-model. However, the fast-slow version of MB-model does not predict a canard explosion phenomenon.

We consider a predator-prey model in a two-patch environment and assume that migration between patches is faster than prey growth, predator mortality and predator-prey interactions. Prey (resp. predator) migration rates are considered to be predator (resp. prey) density-dependent. Prey leave a patch at a migration rate proportional to the local predator density. Predators leave a patch at a migration rate inversely proportional to local prey population density. Taking advantage of the two different time scales, we use aggregation methods to obtain a reduced (aggregated) model governing the total prey and predator densities. First, we show that for a large class of density-dependent migration rules for predators and prey there exists a unique and stable equilibrium for migration. Second, a numerical bifurcation analysis is presented. We show that bifurcation diagrams obtained from the complete and aggregated models are consistent with each other for reasonable values of the ratio between the two time scales, fast for migration and slow for local demography. Our results show that, under some particular conditions, the density dependence of migrations can generate a limit cycle. Also a co-dim two Bautin bifurcation point is observed in some range of migration parameters and this implies that bistability of an equilibrium and limit cycle is possible. © 2007.

The aim of this work is to study the influence of patch selection on the dynamics of a system describing the interactions between two populations, generically called `population N' and `population P'. Our model may be applied to prey–predator systems as well as to certain host–parasite or parasitoid systems. A situation in which population P affects the spatial distribution of population N is considered. We deal with a heterogeneous environment composed of two spatial patches: population P lives only in patch 1, while individuals belonging to population N migrate between patch 1 and patch 2, which may be a refuge. Therefore they are divided into two patch sub-populations and can migrate according to different migration laws. We make the assumption that the patch change is fast, whereas the growth and interaction processes are slower. We take advantage of the two time scales to perform aggregation methods in order to obtain a global model describing the time evolution of the total populations, at a slow time scale. At first, a migration law which is independent on population P density is considered. In this case the global model is equivalent to the local one, and under certain conditions, population P always gets extinct. Then, the same model, but in which individuals belonging to population N leave patch 1 proportionally to population P density, is studied. This particular behavioral choice leads to a dynamically richer global system, which favors stability and population coexistence. Finally, we study a third example corresponding to the addition of an aggregative behavior of population N on patch 1. This leads to a more complicated situation in which, according to initial conditions, the global system is described by two different aggregated models. Under certain conditions on parameters a stable limit cycle occurs, leading to periodic variations of the total population densities, as well as of the local densities on the spatial patches.