Predator-Prey Simulator: Lotka-Volterra Population Dynamics

simulator intermediate ~10 min
Loading simulation...
Oscillation period ≈ 15.7 time units, prey peak ≈ 44, predator peak ≈ 18

With default parameters (α=1.1, β=0.4, γ=0.4, δ=0.1), the Lotka-Volterra system produces stable oscillations with a period of approximately 15.7 time units. Prey and predator populations cycle out of phase — the predator peak lags behind the prey peak, a quarter-cycle delay that is a hallmark of coupled oscillatory dynamics.

Formula

dx/dt = αx − βxy
dy/dt = δxy − γy
Equilibrium: x* = γ/δ, y* = α/β
Conserved quantity: V = δx − γ·ln(x) + βy − α·ln(y)

The Lotka-Volterra Predator-Prey Model

In 1925, Alfred Lotka proposed a mathematical model for oscillating chemical reactions that was independently rediscovered by Vito Volterra in 1926 to explain fluctuations in Adriatic fish catches. The resulting Lotka-Volterra equations became one of the most influential models in mathematical ecology, providing the theoretical foundation for understanding predator-prey dynamics.

The Equations

The model consists of two coupled differential equations:

dx/dt = αx − βxy — Prey grow exponentially at rate α in the absence of predators. The term βxy represents predation: encounters between predators and prey are proportional to the product of their populations (the mass-action principle).

dy/dt = δxy − γy — Predators die exponentially at rate γ without food. The term δxy represents predator reproduction fueled by successful predation, where δ reflects the efficiency of converting consumed prey into new predators.

Perpetual Oscillations

The most striking feature of the Lotka-Volterra model is its prediction of perpetual, undamped oscillations. The system possesses a conserved quantity (analogous to energy in physics): V = δx − γ·ln(x) + βy − α·ln(y), which remains constant along any trajectory. This conservation law ensures that orbits in the phase plane are closed curves — the populations cycle forever without damping or amplification.

The equilibrium point (x* = γ/δ, y* = α/β) sits at the center of these orbits. Trajectories closer to the equilibrium have smaller amplitude oscillations; those further away swing more wildly. Initial conditions determine which orbit the system follows, but all orbits share the same period — a unique property of the linearized system near equilibrium.

Ecological Insights

The quarter-cycle lag between prey and predator peaks is a key prediction confirmed in many real ecosystems. The classic example is the Hudson's Bay Company fur trade records showing coupled oscillations of snowshoe hare and Canada lynx populations over more than a century. While real ecosystems include complications absent from the basic model — carrying capacity, functional responses, spatial structure — the Lotka-Volterra framework remains the starting point for understanding all predator-prey interactions.

Explore the phase portrait on the right side of the visualization: notice how the closed orbit traces out the perpetual cycle. The shape and size of the orbit change with parameters — increasing predation efficiency (δ) or prey birth rate (α) dramatically alters the dynamics.

FAQ

What is the Lotka-Volterra model?

The Lotka-Volterra equations are a pair of coupled first-order nonlinear differential equations that model the dynamics of two interacting species — a predator and its prey. Independently proposed by Alfred Lotka (1925) and Vito Volterra (1926), the model predicts perpetual oscillations in both populations.

What are the Lotka-Volterra equations?

dx/dt = αx − βxy (prey growth minus predation) and dy/dt = δxy − γy (predator growth from predation minus natural death), where x is the prey population, y is the predator population, and α, β, γ, δ are positive parameters.

Why do predator and prey populations oscillate?

Oscillations arise from a feedback loop: abundant prey fuels predator growth, but increasing predators reduce prey, which then causes predator decline from starvation, allowing prey to recover. This negative feedback with time delays creates sustained cycles.

What is a phase portrait in ecology?

A phase portrait plots predator population (y) against prey population (x), removing time as an explicit variable. In the Lotka-Volterra model, trajectories form closed orbits around the equilibrium point (x* = γ/δ, y* = α/β), with the orbit shape depending on initial conditions.

Sources

Embed

<iframe src="https://homo-deus.com/lab/evolution/predator-prey/embed" width="100%" height="400" frameborder="0"></iframe>
View source on GitHub