Time Scale Separation

Imagine two situations (prey and predator) with the logistic growth :

  1. Prey, changes fast:
dNdt=r(1−N/K)N−βNP
  1. Predator, changes very slowly:
dPdt=γβNP−δP

K is prey carrying capacity
β is capture rate of searching predators
δ is predator death rate
γ is conversion rate - if predator eats one prey, it turns it into γ offpsring
r is prey growth rate

Since the number of predators changes slowly, we can sort of consider P as a constant in formula 1. This gives us the prey "quasi-equilibrium" N^2=r−βPrK, if r>βP.

(using same method as Harvested Logistic Model)

For the predator function, if the prey can be assumed to always be close to its quasi-equilibrium, we can model the predator change over this slower time (years for example):

dPdt=[γβr−βPrK−δ]PdPdt=[(γβK−δ)−γβ2KrP]PdPdt=(γβK−δ)[1−Prγβ2K(γβK−δ)]P

This last one looks like the logistic function again. So we know that we have the intrinsic growth rate γβK−δ and the carrying capacity rγβ2K(γβK−δ)

This shows that the trivial equilibrium P^1=0 is stable if δ>γβK, meaning the predator population will always trend to zero if this is met.

The predator is only viable if γβK>δ - they need to be able to catch enough prey, produce enough offspring and prey needs a large enough carrying capacity to exceed their death rate.

Predator population model

left=-0.5
top=5
---
g = 1
B = 1
K = 3
d = 1.2
R = 0.8

y=(gBK-d) (1-x/((R/gB^2K)gBK-d))x