The Azimuth Project
Logistic equation

Contents

Idea

The logistic equation is a simple model of population growth in conditions where there are limited resources. When the population is low it grows in an approximately exponential way. Then, as the effects of limited resources become important, the growth slows, and approaches a limiting value, the equilibrium population or carrying capacity.

Details

The logistic growth model is

dxdt=rK(Kx)x.{d x \over d t} = {r\over K} (K-x) x \, .

Here xx is the population, which is a function of time tt. KK is the equilibrium population, and rr is the growth rate.

Note that in the limit KK \to \infty, we get the simpler model:

dxdt=rx {d x \over d t} = r x

describing exponential population growth:

x(t)=x 0e rt. x(t) = x_0 e^{r t} \, .

When KK is finite and positive, the logistic model describes population growth that is approximately exponential when the population is much less than KK, but levels off as the population approaches KK. If the population is larger than KK, it will decrease. Every positive solution has

lim t+x(t)=K. lim_{t \to +\infty} x(t) = K \, .

The logistic model can be normalised by rescaling the units of population and time. Define y:=x/Ky\colon = x/K and s:=rts\colon = rt. The result is

dyds=y(1y).{d y \over d s} = y(1-y) \, .

It is easy to find the explicit solution of the logistic equation, since it is a first-order separable differential equation.

The growing solutions are all time-translated versions of the logistic function (See Wikipedia)

y(s)=e s1+e s=11+e s y(s) = {e^s\over 1 + e^s} = {1\over 1 + e^{-s}}

which looks like this:

This function goes from 00 to 11 as tt goes from -\infty to ++\infty. All other growing solutions have the same limiting behavior and are time-translated versions of this one. After rescaling back to the original variables, we have

x(t)=K1+e r(tt 0). x(t) = {K\over 1 + e^{-r(t-t_0)}} \, .

There are also decreasing solutions where x>1x \gt 1 and solutions (irrelevant to population biology) where x<0x \lt 0 decreases explosively to -\infty.

Chaos in the logistic model

Discrete-time version

When the logistic equation is discretised it displays chaos. It is, in fact, the canonical ground for the study of the period-doubling cascade.

Assume that the Euler formula is used to discretise the logistic growth model, that is,

ΔxΔt=rK(Kx)x.{\Delta x \over \Delta t} = {r\over K} (K-x) x \, .

Assume t n+1t n=Δtt_{n+1}-t_n = \Delta t and x n=x(t n)x_n = x(t_n) for all nn. Then,

x n+1=(1+rΔt)x nrΔtKx n 2.x_{n+1} = (1+r\Delta t)x_n - {r\Delta t\over K} x_n^2 \, .

This can be normalised by letting λ=1+rΔt\lambda = 1+r\Delta t:

x n+1=λx n(111/λKx n)x_{n+1} = \lambda x_n\Bigl(1 - {1 - 1/\lambda\over K} x_n\Bigr)

and measuring xx in units of K11/λ{K\over 1-1/\lambda}, that is, x=Ky11/λx = {Ky\over 1-1/\lambda} and one obtains the so-called Logistic Map.

y n+1=λy n(1y n).y_{n+1} = \lambda y_n(1 - y_n)\, .

It is known (Wikipedia, again) that, depending on the value of λ\lambda, this equation has a stable fixed point, multiperiodic stationary states, or displays chaotic behaviour. Note that λ\lambda depends only on the growth rate and the time step, and not on the carrying capacity.

When the logistic map is interpreted as a numerical approximation to the continuous logistic growth model, the time step Δt\Delta t is largely arbitrary and can be chosen to be less than 1/r1/r, so that 1<λ<21 \lt \lambda \lt 2 ensuring numerical stability.

However, it is possible to interpret the continuous logistic model as an approximation to a fundamental discrete model. In that case, Δt\Delta t is not arbitrary (representing, for instance, the periodicity with which a species reproduces - usually annually) and neither is rr, in which case λ=1+rΔt\lambda = 1 + r\Delta t cannot be tuned. If the population reproduces too fast (high values of lambda), then the sizes of successive generations may be chaotic. If λ>3\lambda \gt 3 the stationary state fluctuates among at least two states, and if λ>3.57\lambda \gt 3.57 (approximately) the dynamics is chaotic.

Continuous version

Here chaos arises as a consequence of delayed feedback.

Rewrite the normalised logistic growth model as

dxx=(1x)dt {d x\over x} = (1 - x)d t

and add a delay to the right-hand side, obtaining

dx(t)x(t)=[1x(tΔ)]dt {d x(t)\over x(t)} = \bigl[1 - x(t-\Delta)\bigr]d t

This introduces the possibility of overshoot?, as when x(t)x(t) reaches the carrying capacity 11 the earlier value x(tΔ)x(t-\Delta) is still away from the carrying capacity and so the system continues in the same direction as it is moving in when it touches 1. In fact, the population trend reverses a time Δ\Delta after reaching the carrying capacity. Depending on the speed with which the x=1x=1 level is crossed this may result in growing oscillations as when the constant solution to the discrete version becomes unstable at λ=3\lambda = 3.

To see this, write x=e yx = e^y:

dy(t)dt=1e y(tΔ). {d y(t)\over d t} = 1 - e^{y(t-\Delta)}.

Assuming now that y0y\approx 0,

dy(t)dty(tΔ). {d y(t)\over d t} \approx - y(t-\Delta).

An analysis of this equation is carried out in the article on delayed feedback.

Complementary approaches

More details on the populations age structure has been added by several research groups. Also the Individual based ecology approach has bee used both for complementary and comparative research. See Law,Murrell and Dieckmann.

References

How great an effect does self-generated spatial structure have on logistic population growth? Results are described from an individual based model (IBM) with spatially localised dispersal and competition, and from a deterministic approximation to the IBM describing the dynamics of the first and second spatial moments. The dynamical system incorporates a novel closure that gives a close approximation to the IBM in the presence of strong spatial structure. Population growth given by the spatial logistic equation can differ greatly from that of the non-spatial logistic model. Numerical simulations show that populations may grow more slowly or more rapidly than would be expected from the non-spatial model, and may reach their maximum rate of increase at densities other than half of the carrying capacity. Populations can achieve asymptotic densities substantially greater than or less than the carrying capacity of the non-spatial logistic model,

and can even tend towards extinction. These properties of the spatial logistic equation are caused by local dispersal and competition which affect spatial structure, which in turn affects population growth. Accounting for these local spatial processes brings the theory of single-species population growth a step closer to the growth of real spatially structured populations.

See also Allee effect for a generalization of the logistic equation where the per capita population growth rate becomes smaller for very small populations.

category: ecology