Lotka-Volterra Simulator: Predator-Prey Population Oscillations

simulator intermediate ~10 min
Loading simulation...
T ≈ 5.1 yr — stable predator-prey oscillation

With default parameters, the system exhibits stable oscillations with a period of about 5.1 years. Prey equilibrium is 20 individuals and predator equilibrium is 10 individuals.

Formula

dx/dt = αx − βxy (prey growth minus predation)
dy/dt = δxy − γy (predator growth minus mortality)
T ≈ 2π / √(αγ) (approximate oscillation period)

The Dance of Predator and Prey

In 1925, Alfred Lotka proposed a mathematical model for oscillating chemical reactions that Vito Volterra independently rediscovered in 1926 to explain fish population cycles in the Adriatic Sea. Their coupled differential equations — dx/dt = αx − βxy for prey and dy/dt = δxy − γy for predators — became the foundation of mathematical ecology. The model reveals how two species locked in a consumption relationship generate perpetual, self-sustaining oscillations.

Phase Space and Conservation

Unlike most nonlinear systems, the Lotka-Volterra equations possess a conserved quantity H that constrains trajectories to closed orbits in the prey-predator phase plane. Each orbit corresponds to a different initial condition, and populations trace ellipse-like loops around the equilibrium point. This conservation law means the system never settles to a steady state and never diverges — a mathematically elegant but biologically idealized property.

Ecological Oscillations in Nature

The most celebrated real-world example is the Canadian lynx-snowshoe hare cycle, documented through Hudson Bay Company fur records spanning over a century. Population peaks recur approximately every 10 years, with predator peaks lagging prey peaks by 1-2 years — qualitatively matching Lotka-Volterra predictions. Similar oscillations appear in wolf-moose systems on Isle Royale and in microbial predator-prey communities in laboratory chemostats.

Beyond the Classic Model

Modern extensions incorporate carrying capacity (logistic prey growth), functional responses (Holling type II/III), spatial diffusion, stochastic noise, and multiple interacting species. The Rosenzweig-MacArthur model adds prey self-limitation and predator saturation, producing limit cycles and the paradox of enrichment. Despite its simplicity, the original Lotka-Volterra framework remains the starting point for understanding all predator-prey dynamics.

FAQ

What is the Lotka-Volterra model?

The Lotka-Volterra equations are a pair of first-order nonlinear differential equations describing predator-prey interactions. Proposed independently by Alfred Lotka (1925) and Vito Volterra (1926), they predict perpetual oscillations where prey and predator populations cycle out of phase.

Why do populations oscillate in the Lotka-Volterra model?

Oscillations arise from a negative feedback loop: abundant prey fuels predator growth, but increased predators deplete prey, causing predator starvation and decline, which lets prey recover. This cycle repeats indefinitely in the idealized model.

Is the Lotka-Volterra model realistic?

The basic model captures qualitative oscillatory behavior observed in real systems (e.g., lynx-hare cycles), but it omits carrying capacity, stochastic effects, spatial structure, and multi-species interactions. Extensions like the Rosenzweig-MacArthur model address some of these limitations.

What is the conserved quantity in Lotka-Volterra?

The system has a conserved quantity H = δx + βy − γ ln(x) − α ln(y), which remains constant along trajectories. This makes the orbits closed curves in phase space, preventing the system from reaching a stable equilibrium or diverging.

Sources

Embed

<iframe src="https://homo-deus.com/lab/mathematical-biology/lotka-volterra/embed" width="100%" height="400" frameborder="0"></iframe>
View source on GitHub