Modeling Predator-Prey Dynamics

Simulating Predator-Prey Dynamics With The Lotka-Volterra Model

Decoding Biology Shorts is a new weekly newsletter sharing tips, tricks, and lessons in bioinformatics. Enjoyed this piece? Show your support by tapping the ❤️ in the header above. Your small gesture goes a long way in helping me understand what resonates with you and in growing this newsletter. Thank you!

---

Modeling Predator-Prey Dynamics

The Lotka-Volterra (LV) model, often called the predator-prey model, is a cornerstone of mathematical ecology. It provides a powerful framework for understanding the dynamic interactions between two species within an ecosystem: one acting as a predator and the other as prey. This model is particularly useful for illustrating the principles of population dynamics and offers a foundation for studying more complex biological systems. For example, using the LV model we can better understand ecosystem stability and species coexistence, design conservation strategies by predicting the impact of introducing or removing species, and build intuition for or complex population models involving multiple species or environmental factors.

Simulating Predator-Prey Dynamics With The Lotka-Volterra Model

The Lotka-Volterra model uses a system of differential equations to describe the cyclical population dynamics of two interacting species. These equations capture how the populations of predators and prey influence each other over time. For example, the prey population grows exponentially in the absence of predators but is diminished by predation. Additionally, the predator population depends on the availability of prey for sustenance and declines in the absence of sufficient prey.

We can represent this system of predator-prey relationships with the following system of equations where dD/dt and dL/dt are the change in prey and predator populations over time respectively, b is the birth rate of prey, p is the predation rate, r is the predator reproduction rate per prey consumed, d is the death rate of predators, D is the prey population, and L is the predator population, as demonstrated below:

\(\frac{dD}{dt} = bD -pDL \)

\(\frac{dL}{dt}=rDL-dL\)

These equations predict cyclical patterns of population booms and crashes, often observed in nature.

Below is a Python implementation of the Lotka-Volterra model. This code simulates and visualizes the population dynamics of a predator (lions) and prey (deer):

Note: You can find the code for the Lotka-Volterra model on GitHub here.

import numpy as np
from scipy.integrate import odeint
import matplotlib.pyplot as plt

# Model parameters
b = 0.5  # Birth rate of prey
p = 0.3  # Predation rate
r = 0.5  # Predator reproduction rate per prey consumed
d = 0.1  # Death rate of predators

D0 = 20  # Initial prey population
L0 = 20  # Initial predator population

t = np.linspace(0, 1000, 1000)  # Time range for simulation

def deriv(y, t, b, p, r, d):
    D, L = y
    dDdt = b * D - p * D * L
    dLdt = r * D * L - d * L
    return dDdt, dLdt

# Initial conditions
y0 = D0, L0

# Solve differential equations
ret = odeint(deriv, y0, t, args=(b, p, r, d))
D, L = ret.T

# Plot the results
fig = plt.figure(facecolor='w')
ax = fig.add_subplot(111, facecolor='#dddddd', axisbelow=True)
ax.plot(t, D, 'b', alpha=0.5, lw=2, label='Deers')
ax.plot(t, L, 'r', alpha=0.5, lw=2, label='Lions')
ax.set_xlabel('Time (days)')
ax.set_ylabel('Population (thousands)')
ax.set_ylim(0, 20)
ax.grid(which='major', c='w', lw=2, ls='-')
legend = ax.legend()
legend.get_frame().set_alpha(0.5)
plt.show()

Which produces the following output:

The code above can be modified and repurposed for various purposes. For example, we can change the initial conditions D0 and L0 to explore how different starting populations of prey and predators impact the populations dynamics, we can expand the system of equations to include additional predator and prey species, allowing us to simulate more complex ecosystems, or we can incorporate external factors and introduce terms to represent seasonal effects, habitat changes, or human interventions.

Closing Thoughts

Working with the Lotka-Volterra model allows us to visualize the interdependence of species populations, understand the sensitivity of ecosystems to parameter changes, and build intuition for real-world ecological applications, such as managing invasive species or conserving endangered ones. By experimenting with the provided code, we also gain hands-on experience with differential equations and computational modeling—key skills in bioinformatics and systems biology.

---