Cellular Automata/Examples of Plankton and Fish Dynamics

From Wikibooks, open books for an open world
Jump to navigation Jump to search

Definition of prey-predator systems[edit | edit source]

Prey-predator systems are a well known biological example of excitable media. This example is based upon an article (see references) about plankton and fish dynamics. The dynamics are given as a PDE system. The purpose of this example is to observe pattern formation and other dynamic properties of prey-predator systems using both finite difference equations and cellular automata.

PDE describing prey-predator systems[edit | edit source]

We define the prey-predator system using the a system of general reaction-diffusion differential equations.

The first equation describes the prey population and the seccond equation describes the predator population . Nonlinear functions and describe local prey-predator dynamics, without them the system is decomposed into two separate diffusion equations.

We are interested in prey-predator dynamics of aquatic systems, where phytoplankton (small plants, food producers) is the prey and zooplankton (small animals) is the predator. Since we observe only slowly moving plankton the diffusion coefficient depends mainly on marine turbulence. The two nonlinear functions describing prey-predator dynamics are

describes local growth and mortality of the prey, describes the predation, is food utilization and is the mortality rate of the predator. The selection of functions , and parameters , depends on the observed population. Here a simple example is used, for a detailed description of the equations see the references.

is the maximal growth rate of the prey, is the carrying capacity of the prey population, and H is the half-saturation abundance of prey.

Algebraic analysis[edit | edit source]

The first step is to introduce dimensionless variables

and the new dimensionless parameters are

the PDE system with only dimensionless quantities is (tildes are omitted)

Local Dynamics[edit | edit source]

Local dynamics cycle (starting in (u*,v*))

The local dynamics of the system are defined by the nonlinear reaction part of equations

We first observe stationary conditions There are three stationary states

  • total extinction is a sedal point for
  • predator extinction is a sedal point if is in the biologically meaningful region and a stable point otherwise
  • and coexistence can be a sedal or stationary point

Although the reaction dynamics depend on three parameters is the most important. If

than the population is in a biologically meaningful region . If

than there is a stable limit cycle surrounding the steady state. There is a figure of such a cycle where the path begins near from the stationary point and the parameters are , and .

Numerical simulation using finite difference equations[edit | edit source]

A direct approach to simulate prey-predator dynamics in the plankton example is to discretize the reaction-diffusion differential equation, and to solve it numerically. The dynamics is first observed in an one dimensional space and than on two dimensions. In both cases using a zero-flux boundary .

The PDE discretization method used is FTCS, which is the simplest. It is capable of producing complex behavior we would like to observe, but it has stability and convergence problems discussed below.

One dimension[edit | edit source]

The simulation is first performed in an one dimensional space of elements. The discretization parameters are

The prey-predator system parameters are

For initial conditions the uniform stationary point is used

a gradual increase of predator population from the left to the right is added to ignite the spatial dynamics, the parameters are

There are two figures showing the system at times 0, 640, 2600. At a certain point in time and space the species distribution becomes chaotic. The chaotic part shows a fast growth.

Spatial distributions of populations in the 1D case at time t=640
Spatial distributions of populations in the 1D case at time t=2640


The system is not converging by decreasing the time difference . The source of the problem can be seen by observing the local dynamics of and inside the chaotic region from time to using two different time differences and . With decreasing the time difference the influence of diffusion grows over the influence of local dynamics. On the above figure the effect would be seen as a slower spreading of the chaotic region.

Local dynamics (u(t),v(t)) observed at x=1500 t=0-640 with Δt=0.20
Local dynamics (u(t),v(t)) observed at x=1500 t=0-640 with Δt=0.02


Two dimensions[edit | edit source]

The simulation is then performed in a two-dimensional space of elements. The discretization parameters are

The prey-predator system parameters are

For initial conditions the uniform stationary point is used

a gradual increase of predator population is added to ignite the spatial dynamics, the parameters are

2D dynamics are observed on time points . Spiral patterns can be seen arising from the initial gradient, at the end patterns brake and the system ends in chaos.

t=0
t=100
t=150
t=200
t=300
t=400
t=1000


Another problem arises because the FTCS method is not unconditionally stable. Some results may become negative, and the functions describing local dynamics diverge very fast in the negative quadrants.

References[edit | edit source]

  1. Alexander B. Medvinsky, Sergei V. Petrovskii, Irene A. Tikhonova, Horst Malchow, Bai-Lian Li, Spatiotemporal Complexity of Plankton and Fish Dynamics, SIAM Review, Volume 44, Number 3, pp. 311-370
  1. Marcus R. Garvie, Finite difference schemes for reaction-diffusion equations modeling predator-prey interactions in MATLAB PDF