Beverton-Holt Model

The Beverton-Holt Model is a foundational discrete-time, density-dependent mathematical model of population dynamics and fisheries stock-recruitment developed in 1957 by British fisheries scientists Raymond J. H. Beverton (1922–1995) and Sidney J. Holt (1926–2019) in their landmark monograph On the Dynamics of Exploited Fish Populations. Originally formulated to predict the number of juvenile recruits entering commercial North Sea demersal fisheries (such as European plaice Pleuronectes platessa and haddock Melanogrammus aeglefinus) as a function of adult spawning stock biomass, the Beverton-Holt equation serves throughout theoretical ecology as the canonical discrete-time analogue of the continuous-time Pierre-François Verhulst logistic differential equation.

Unlike the overcompensatory Ricker Model—which generates humped recruitment curves, multi-year oscillations, and deterministic chaos at high intrinsic growth rates—the Beverton-Holt model describes compensatory contest competition: as parental population density or initial larval cohort density increases toward infinity, the number of surviving offspring increases monotonically toward a finite horizontal asymptote representing the environmental carrying capacity (K). Because its functional form is mathematically equivalent to a Michaelis-Menten or Monod rectangular hyperbola, the Beverton-Holt map exhibits globally stable, non-oscillatory equilibrium behavior across all positive parameter values.

Depending on whether the model is applied to annual plant and insect generations in theoretical ecology or to spawning biomass and juvenile recruitment in fisheries stock assessment, the Beverton-Holt difference equation is written in three algebraically equivalent forms:

1. Discrete-Time Population Ecology Parameterization (R₀, M, and K)

Let N_t denote the population density of adults (or seeds) in generation t, and let N_t+1 denote the population density in generation t + 1. The standard ecological Beverton-Holt map is:

N_t+1 = ( R₀ · N_t ) / ( 1 + N_t / M )

where R₀ (R₀ > 0) is the maximum per-capita net reproductive rate (geometric finite rate of increase per generation) realized at vanishingly low population density (N_t → 0), and M (M > 0) is the half-saturation threshold density at which intraspecific competition reduces per-capita survival to half of its maximum value (N_t+1 / N_t = 0.5 R₀). Setting N_t+1 = N_t = K for a growing population (R₀ > 1) yields the non-trivial equilibrium carrying capacity (K):

K = M · (R₀ − 1) ,     so that     N_t+1 = ( R₀ · N_t ) / [ 1 + ( (R₀ − 1) / K ) N_t ]

2. Classical Fisheries Stock-Recruitment Parameterization (α and β)

In fisheries stock assessment (Beverton & Holt, 1957; Hilborn & Walters, 1992), let S represent the spawning stock biomass (or total eggs spawned) and R represent the number of recruits surviving to enter the fishery:

R(S) = ( α · S ) / ( 1 + β · S )

Here α = dR / dS|_S=0 is the initial slope at the origin (maximum recruits produced per unit spawner biomass at low stock density), β measures the intensity of density-dependent mortality among pre-recruits, and R_max = α / β is the asymptotic maximum recruitment ceiling as S → ∞.

3. Steepness Parameterization (h) in Stock Assessment

To compare resilience across fish stocks of differing sizes, Mace and Doonan (1988) and Myers et al. (1999) introduced the dimensionless steepness parameter (h), defined as the fraction of unfished equilibrium recruitment (R₀) produced when spawning stock biomass is reduced to 20% of its unfished virgin level (S = 0.20 S₀):

R(S) = ( 4 · h · R₀ · S ) / [ S₀ (1 − h) + S (5h − 1) ]

Under the Beverton-Holt model, steepness h ranges between 0.20 (linear density independence with zero compensatory buffer) and 1.0 (complete recruitment compensation). Demersal flatfishes and gadoids typically exhibit high Beverton-Holt steepness (h ≈ 0.75 to 0.90), whereas slow-growing K-selected elasmobranchs and rockfishes (Sebastes) exhibit low steepness (h ≈ 0.35 to 0.55).

Continuous Within-Cohort Mortality Derivation

Beverton and Holt (1957) derived their equation mechanistically from a continuous-time differential equation governing linear density-dependent larval mortality across the juvenile nursery window (from hatching t = 0 to recruitment age t = T). Suppose a cohort of newly hatched larvae starts at initial abundance N(0) = E · S and suffers both a density-independent mortality rate μ₁ and a density-dependent mortality rate μ₂ · N(t) proportional to instantaneous cohort density N(t):

dN(t) / dt = − [ μ₁ + μ₂ · N(t) ] N(t)

Solving this Bernoulli differential equation via the reciprocal substitution u(t) = 1 / N(t) and integrating from t = 0 to t = T gives the exact Beverton-Holt hyperbola R = (α · S) / (1 + β · S).

Nicholson's Contest Competition vs. Scramble Competition

In the ecological taxonomy of intraspecific competition introduced by A. J. Nicholson (1954) and formalized by Michael P. Hassell (1975), the Beverton-Holt model represents pure contest competition (with Hassell exponent b = 1):

• Under contest competition, the nursery habitat is partitioned into a fixed number of discrete refuges or defended territories (such as coral reef crevices for damselfish recruits or stream riffle territories for salmonid parr, compare the ideal free distribution).
• Winning individuals secure a territory and obtain the full resource quota needed to mature, whereas excess individuals excluded from territories perish or emigrate.
• Consequently, regardless of how many larvae arrive, surviving recruitment saturates smoothly at the ceiling of available territories (R_max = α / β).

Linear Reciprocal Transform and Eigenvalue Stability

Taking the reciprocal of both sides of the Beverton-Holt map yields a linear difference equation in inverse population density (1 / N_t+1 − 1 / K = R₀⁻¹ [1 / N_t − 1 / K]), whose closed-form analytical solution for any generation t is:

N_t = K / [ 1 + ( (K − N₀) / N₀ ) · R₀^−t ]

This trajectory matches the continuous Verhulst logistic growth curve evaluated at integer time steps with r = ln(R₀). Evaluating the derivative at the positive equilibrium N* = K yields the stability eigenvalue λ = f′(K) = 1 / R₀. Because 0 < λ < 1 for all R₀ > 1, perturbations from carrying capacity decay monotonically without oscillations or chaos.

Multispecies Annual Plant Coexistence and Allee Extensions

In community ecology, the multispecies Beverton-Holt model (Leslie-Gower discrete competition model, applied by Peter Chesson and Jonathan M. Levine to annual plant communities) incorporates seed-bank survival (s_i), germination (g_i), intrinsic seed yield (λ_i), and competition coefficients (α_ij) to extend Lotka-Volterra equations and resource partitioning to discrete generations:

N_i,t+1 = [ (1 − g_i)s_i + ( g_i λ_i ) / ( 1 + ∑_j α_ij g_j N_j,t ) ] N_i,t

Furthermore, generalizing the stock-recruitment exponent to R(S) = α S^c / (1 + β S^c) with c > 1 incorporates a depensatory Allee effect at low spawning biomass, modeling delayed recovery in collapsed marine stocks.

1. Beverton, R. J. H., & Holt, S. J. (1957). On the Dynamics of Exploited Fish Populations (Fishery Investigations Series II, Vol. XIX). Her Majesty's Stationery Office, London. https://doi.org/10.1007/978-94-011-2106-4
2. Hassell, M. P. (1975). Density-dependence in single-species populations. Journal of Animal Ecology, 44(1), 283–295. https://doi.org/10.2307/3863
3. Myers, R. A., Bowen, K. G., & Barrowman, N. J. (1999). Maximum reproductive rate of fish at low population sizes. Canadian Journal of Fisheries and Aquatic Sciences, 56(12), 2404–2419. https://doi.org/10.1139/f99-201
4. Cushing, J. M., & Henson, S. M. (2001). Global dynamics of some periodically forced, monotone difference equations. Journal of Difference Equations and Applications, 7(6), 859–872. https://doi.org/10.1080/10236190108808308
5. Levine, J. M., & HilleRisLambers, J. (2009). The importance of niches for the maintenance of species diversity. Nature, 461(7261), 254–257. https://doi.org/10.1038/nature08251