The Lotka–Volterra Model
~11 min read · Tier toggle changes depth (T1 plain · T2 standard · T3 rigour)
You should be comfortable with basic differential calculus, specifically solving simple first-order differential equations, and have a general understanding of population growth models such as exponential and logistic growth.
Why this matters
The Lotka–Volterra model provides the foundational mathematical framework for understanding how predator and prey populations interact over time.
In ecology, interactions between species are rarely isolated; they form complex networks that drive community structure and stability. While previous lectures established that predation reduces prey abundance and that predators depend on prey for energy, we now need a quantitative tool to predict the resulting dynamics. The Lotka–Volterra model is the simplest mathematical representation of this interaction. It demonstrates that even a system with only two species and simple assumptions can produce complex, oscillatory behaviour. This insight is crucial because it challenges the intuitive notion that predator and prey populations should settle into a static equilibrium. Instead, it reveals that biological systems often exhibit dynamic stability, where populations fluctuate around a mean value. Understanding this model is essential for interpreting real-world data on population cycles, managing fisheries, and controlling pest populations, as it highlights the inherent lag and feedback mechanisms that drive ecological change.
The Intuitive Logic of the Model
To understand the Lotka–Volterra model, we must first strip away the complexity of the real world and focus on the core assumptions that make the mathematics tractable. The model assumes a closed system where no individuals enter or leave the population, and where the only factor affecting birth and death rates is the interaction between the two species. For the prey population, denoted as \( N \), we assume that in the absence of predators, they grow exponentially at a rate \( r \). This is a simplification, as real prey populations are usually limited by resources, but it allows us to isolate the effect of predation. The predation term is assumed to be proportional to the encounter rate between predators and prey. This is the "mass action" assumption, meaning that if you double the number of predators or the number of prey, the number of predation events doubles. Consequently, the loss of prey is proportional to the product of the two populations, \( NP \). For the predator population, denoted as \( P \), we assume that in the absence of prey, they die off exponentially at a rate \( m \). Their growth is driven entirely by the consumption of prey. A fraction \( b \) of the energy gained from eating prey is converted into new predator biomass. Therefore, the predator growth rate is also proportional to the product \( NP \). The key insight from these assumptions is that the system is tightly coupled. An increase in prey leads to an increase in predators, which then leads to a decrease in prey, which leads to a decrease in predators, and the cycle repeats. This feedback loop is the engine of the oscillations. It is important to note that this model does not include intraspecific competition (density dependence) within the prey population, which is why the oscillations are neutral and do not dampen or explode over time. This is a major limitation, but it provides the baseline against which more complex models are compared. The model predicts that the populations will cycle indefinitely around a stable equilibrium point, provided the parameters remain constant. This equilibrium is not a state of zero growth, but a state where the birth and death rates for both species balance out perfectly. The timing of these cycles is determined by the interaction rates and the intrinsic growth and death rates of the species. This intuitive understanding sets the stage for the formal mathematical description, where we can see exactly how these assumptions translate into the differential equations that govern the system's behaviour.
Formulating the Differential Equations
The mathematical formulation of the Lotka–Volterra model consists of two coupled first-order ordinary differential equations. The first equation describes the rate of change of the prey population \( N \) with respect to time \( t \). It is given by \( \frac{dN}{dt}=rN-aNP \). Here, \( r \) is the intrinsic growth rate of the prey, representing the per capita birth rate minus the natural death rate in the absence of predators. The term \( aNP \) represents the loss of prey due to predation. The parameter \( a \) is the attack rate, which quantifies the efficiency of predators in finding and consuming prey. The second equation describes the rate of change of the predator population \( P \). It is given by \( \frac{dP}{dt}=baNP-mP \). In this equation, \( m \) is the natural mortality rate of the predators. The term \( baNP \) represents the growth of the predator population due to predation. The parameter \( b \) is the conversion efficiency, representing the fraction of consumed prey biomass that is converted into new predator biomass. The product \( ba \) thus represents the net efficiency of energy transfer from prey to predators. These two equations are coupled because the rate of change of one population depends on the size of the other. To find the equilibrium points, we set both derivatives to zero. For the prey equation, \( rN-aNP=0 \) implies \( N(r-aP)=0 \). This gives two solutions: \( N=0 \) (extinction of prey) or \( P=r/a \). For the predator equation, \( baNP-mP=0 \) implies \( P(baN-m)=0 \). This gives two solutions: \( P=0 \) (extinction of predators) or \( N=m/(ba) \). The non-trivial equilibrium point, where both species coexist, is therefore \( N^* = \frac{m}{ba} \) and \( P^* = \frac{r}{a} \). This result is significant because it shows that the equilibrium prey density is determined by the predator's mortality rate and conversion efficiency, while the equilibrium predator density is determined by the prey's growth rate and the attack rate. This mathematical structure reveals the delicate balance required for coexistence. If the parameters shift, the equilibrium point moves, potentially leading to different dynamic outcomes. The derivation of these equilibrium points is a standard exercise in mathematical ecology and provides a benchmark for evaluating the impact of environmental changes on predator-prey systems.
Analysis of Population Cycles
While the equilibrium point represents a state of zero net growth, the Lotka–Volterra model predicts that populations will not remain at this point but will oscillate around it. To understand the nature of these oscillations, we can analyse the system's behaviour near the equilibrium. By linearising the equations around the equilibrium point \( (N^*, P^*) \), we can determine the stability of the system. The Jacobian matrix of the system, evaluated at the equilibrium, has eigenvalues that are purely imaginary. This indicates that the equilibrium is a centre, meaning that the trajectories in the phase plane are closed orbits. These closed orbits correspond to periodic cycles of predator and prey abundance. For small cycles close to the equilibrium the period is approximately \( \frac{2\pi}{\sqrt{rm}} \); larger cycles take somewhat longer. That the cycles neither grow nor shrink is a direct consequence of the lack of density dependence in the prey population. In real ecosystems, density dependence usually causes the oscillations to dampen or grow, but in this idealised model, the energy is conserved in a sense, leading to perpetual cycles. The phase relationship between the two populations is also important. The prey population peaks first, followed by the predator population. This lag is a critical feature of predator-prey dynamics and is often observed in empirical data. The amplitude of the cycles depends on the initial conditions of the system. If the initial populations are far from the equilibrium, the cycles will have a larger amplitude. However, because the orbits are closed, the system will always return to its initial state after one full cycle. This neutral stability is a key feature of the model and highlights the sensitivity of the system to initial conditions. Understanding these cycles is crucial for interpreting long-term population data and for predicting the impact of perturbations such as harvesting or environmental change. The model provides a theoretical framework for understanding why predator and prey populations often exhibit correlated fluctuations, with the predator population lagging behind the prey population. This lag is a direct result of the time required for predators to respond to changes in prey abundance through reproduction and mortality.
Limitations and Extensions
Despite its elegance and simplicity, the Lotka–Volterra model has several significant limitations that must be considered when applying it to real-world ecological systems. The most critical limitation is the assumption of exponential prey growth in the absence of predators. In reality, prey populations are almost always limited by resources, leading to logistic growth. When logistic growth is incorporated into the model, the dynamics change significantly. The equilibrium point becomes stable, and the oscillations dampen over time, eventually settling at the equilibrium. This is a more realistic representation of many natural systems, where populations do not oscillate indefinitely but rather fluctuate around a stable mean. Another limitation is the assumption of a single predator and a single prey species. In reality, ecosystems are complex, with multiple predator and prey species interacting in food webs. The presence of multiple species can stabilise or destabilise the system, depending on the nature of the interactions. For example, if predators can switch between prey species, this can reduce the amplitude of oscillations and increase stability. Additionally, the model assumes that the environment is constant, with no seasonal variation or stochastic events. In reality, environmental fluctuations can drive population dynamics, leading to irregular cycles or even extinction. The model also assumes that the attack rate and conversion efficiency are constant, which is not always the case. Predators may change their hunting behaviour in response to prey density, and conversion efficiency may vary with the quality of the prey. Despite these limitations, the Lotka–Volterra model remains a valuable tool for understanding the basic principles of predator-prey dynamics. It provides a baseline for more complex models and helps ecologists identify the key factors that drive population cycles. By understanding the assumptions and limitations of the model, ecologists can better interpret empirical data and develop more realistic models of ecosystem dynamics. The model also serves as a pedagogical tool, introducing students to the concepts of coupled differential equations, equilibrium analysis, and stability. It highlights the power of mathematical modelling in ecology and the importance of making explicit assumptions when constructing models. In summary, while the Lotka–Volterra model is a simplification of reality, it captures the essential features of predator-prey interactions and provides a foundation for more advanced ecological theory.
The concept
Predator-prey dynamics
Coupled cycles of predator and prey abundance.
A common misconception is that the Lotka–Volterra model predicts stable, damped oscillations that eventually settle at the equilibrium point. In reality, the model predicts neutral, undamped oscillations that persist indefinitely. This is because the model lacks density dependence in the prey population, which is the primary mechanism for damping oscillations in real ecosystems. Confusing neutral stability with asymptotic stability can lead to incorrect predictions about long-term population behaviour. It is crucial to recognise that the Lotka–Volterra model is a special case of a broader class of predator-prey models, and its predictions should be interpreted with caution when applied to real-world systems. C
Worked example
Problem. Consider a predator-prey system with the following parameters: prey intrinsic growth rate \( r = 0.5 \ \mathrm{yr^{-1}} \), predator natural mortality rate \( m = 0.2 \ \mathrm{yr^{-1}} \), attack rate \( a = 0.1 \ \mathrm{yr^{-1}} \), and conversion efficiency \( b = 0.5 \). Calculate the equilibrium densities of the prey and predator populations, and estimate the period of the population cycles.
- Identify the equilibrium conditions by setting the derivatives to zero: \( \frac{dN}{dt}=0 \) and \( \frac{dP}{dt}=0 \).
- Solve for the equilibrium predator density \( P^* \) using the prey equation: \( P^* = \frac{r}{a} \).
- Substitute the given values into the equation for \( P^* \): \( P^* = \frac{0.5}{0.1} = 5 \).
- Solve for the equilibrium prey density \( N^* \) using the predator equation: \( N^* = \frac{m}{ba} \).
- Substitute the given values into the equation for \( N^* \): \( N^* = \frac{0.2}{0.5 \times 0.1} = \frac{0.2}{0.05} = 4 \).
- Estimate the period of the cycles using the formula \( T \approx \frac{2\pi}{\sqrt{rm}} \).
- Substitute the given values into the period formula: \( T \approx \frac{2\pi}{\sqrt{0.5 \times 0.2}} = \frac{2\pi}{\sqrt{0.1}} \approx \frac{6.28}{0.316} \approx 19.9 \ \mathrm{yr} \).
Units check. The equilibrium densities are dimensionless in this context, as they represent relative abundances, and the period is in years, consistent with the time units of the rate parameters.
Three quick tests
- What is the primary assumption about prey growth in the Lotka–Volterra model?
Answer
Prey grow exponentially in the absence of predators. - How does the equilibrium predator density depend on the prey intrinsic growth rate?
Answer
It is directly proportional to the prey intrinsic growth rate. - What is the nature of the oscillations predicted by the Lotka–Volterra model?
Answer
They are neutral and undamped, persisting indefinitely.
Back to BU-105 →
Every lecture and every concept in this unit, with prerequisites and threads.
Lecture notes checked line by line by the editor; every number in the worked example was recomputed.