Numerical Response

The numerical response is the quantitative change in the population density of a predator, herbivore, or parasitoid in response to variations in the population density of its prey, host, or food resource within an ecosystem (Solomon, 1949; Holling, 1959; Beddington, Free, & Lawton, 1976). Coined by British ecologist Maurice E. Solomon in 1949 alongside the functional response and expanded by Crawford S. Holling (1959) and John R. Beddington, C. A. Free, and John H. Lawton (1976), the numerical response explains how local consumer abundance (P) tracks resource abundance (N) through two distinct ecological mechanisms operating on different spatial and temporal scales: the demographic (reproductive) numerical response and the aggregative (behavioral) numerical response.

In trophic and population ecology, neither the functional response (how many prey an individual predator eats per unit time, f(N)) nor the numerical response (how many predators are present per unit area at a given prey density, P(N)) can independently determine whether a predator guild will regulate a prey population. Instead, the total predatory response—the product of per-capita consumption and predator density—governs whether natural enemies stabilize prey equilibria, drive delayed population cycles across the biosphere, or generate indirect apparent competition among co-occurring prey species.

Ecologists partition the numerical response into two fundamental processes that operate across contrasting ecological timescales:

1. The Demographic (Reproductive and Survival) Numerical Response

The demographic numerical response (or intergenerational numerical response) is the change in predator population size driven by shifts in per-capita birth rates (fecundity), juvenile developmental survival, and adult starvation mortality as a consequence of per-capita prey ingestion f(N).

In continuous-time consumer–resource differential equations—such as the Lotka–Volterra equations and the Rosenzweig–MacArthur model—the per-capita population growth rate of the predator (1/P · dP/dt) is modeled as a function g(N) of the energy and macronutrients assimilated from consumed prey:

• Predator Population Rate of Change: dP / dt = P · g(N) = P [ e · f(N) − m ]
• With Type II Functional Response: dP / dt = P [ e · aN / (1 + ahN) − m ]

where f(N) is the predator's per-capita functional response (governed by attack rate a and handling time h), e is the trophic conversion efficiency converting consumed prey biomass into new predator offspring, and m is the per-capita metabolic maintenance and background mortality rate of the predator.

Setting dP / dt = 0 yields the zero-net-growth isocline (ZNGI) of the predator—the minimum threshold prey density N^∗ required for predator births to balance predator deaths:

• Minimum Sustaining Prey Threshold (N^∗): N^∗ = m / [ a (e − mh) ]

When prey density falls below this threshold (N < N^∗), per-capita energy intake cannot cover basal metabolic maintenance; female predators cease reproduction, neonate starvation rises, and the demographic numerical response is negative (dP/dt < 0). When prey density exceeds N^∗, surplus assimilated energy fuels clutch size expansion and higher juvenile survival, yielding a positive demographic numerical response that plateaus at a physiological maximum intrinsic growth rate r_P,max = (e / h) − m once individual feeding rates saturate (f(N) → 1/h). Because gestation, egg development, and maturation require weeks to years, the demographic numerical response exhibits an inherent generational time lag (τ) relative to prey fluctuations.

2. The Aggregative (Behavioral Migration) Numerical Response

The aggregative numerical response (or within-generation numerical response) is the rapid spatial redistribution of mobile predators, parasitoids, or herbivores that congregate in habitat patches of high prey density and emigrate from patches of low prey density (Hassell & May, 1974; Readshaw, 1973). Unlike the demographic response, the aggregative response occurs without any change in total regional predator abundance, operating within hours or days as foraging animals follow the ideal free distribution across a heterogeneous landscape.

When foraging predators sample multiple habitat patches differing in prey density N_i, individuals prolong their patch residence time (via area-restricted search and increased turning angles after prey capture) inside high-density prey patches while departing rapidly from depleted patches. In the Hassell–May (1974) aggregation model and the stabilized Nicholson–Bailey model, spatial aggregation of searching parasitoids or predators toward high-density host patches concentrates mortality in outbreak epicenters while leaving low-density patches as partial spatial refugia.

In 1959, C. S. Holling synthesized Solomon's two concepts by defining the total response (E(N)) as the total number of prey killed per unit area per unit time by the entire predator population at prey density N:

• Total Predatory Response: E(N) = f(N) × P(N)
• Percentage (Proportional) Prey Mortality: M(N) = E(N) / N = [ f(N) / N ] × P(N)

For predators to regulate a prey population around a stable equilibrium N_eq, the proportional prey mortality M(N) must be positively density-dependent (dM(N) / dN > 0) across low-to-moderate prey densities:

1. Complementary Stabilization When Functional Response Is Type II: Although a Type II hyperbolic functional response alone produces inversely density-dependent per-capita predation (f(N) / N declines with N), a steep, unlagged aggregative or demographic numerical response (P(N) rising rapidly with N) allows the product [ f(N) / N ] × P(N) to increase with N across low-to-moderate prey densities, generating a dome-shaped percentage mortality curve M(N).
2. The Escape Threshold ("Predator Pit"): Because both f(N) (limited by handling time 1/h) and P(N) (limited by territoriality, nesting sites, or maximum fecundity) saturate at high prey densities, percentage mortality M(N) peaks at an intermediate prey density N_peak and declines at higher densities. If an insect pest—such as the eastern spruce budworm (Choristoneura fumiferana) in boreal forests—is pushed above an escape threshold N_escape by favorable weather or propagule pressure, it outstrips the combined functional and numerical ceiling of avian and mammalian predators, erupting into a landscape-scale outbreak.

1. Resident Generalists vs. Nomadic Specialists in Fennoscandian Vole Cycles

Across Fennoscandian boreal forests, Erkki Korpimäki and Kai Norrdahl (1991) contrasted how mammalian and avian predators respond to 3-to-4-year population cycles of Microtus and Myodes voles:

• Resident Mammalian Specialists (Least Weasels, Mustela nivalis): Weasels have high reproductive output but cannot disperse long distances. Their numerical response is purely demographic and lags 6 to 9 months behind vole peaks, driving delayed density-dependent mortality that deepens the crash phase of the multi-annual vole cycle.
• Nomadic Avian Predators (Kestrels Falco tinnunculus and Short-Eared Owls Asio flammeus): Because these raptors fly hundreds of kilometers to settle and breed wherever spring vole densities are high, their local numerical response is aggregative and unlagged, tracking vole density linearly without a generational delay.

2. Biological Control and Intraguild Interference

In agricultural biological control, selecting an effective natural enemy requires matching both its functional response and its numerical response to the target pest. Classic successes—such as the vedalia ladybird beetle (Rodolia cardinalis) controlling cottony cushion scale (Icerya purchasi) in California citrus groves—depend on enemies with short generation times whose demographic and aggregative numerical responses outpace pest reproduction, provided that intraguild predation among competing biocontrol agents does not suppress the specialist consumer.

1. Solomon, M. E. (1949). "The natural control of animal populations." Journal of Animal Ecology, 18(1), 1–35. https://doi.org/10.2307/1578
2. Holling, C. S. (1959). "The components of predation as revealed by a study of small-mammal predation of the European pine sawfly." The Canadian Entomologist, 91(5), 293–320. https://doi.org/10.4039/Ent91293-5
3. Beddington, J. R., Free, C. A., & Lawton, J. H. (1976). "Concepts of stability and resilience in predator-prey models." Journal of Animal Ecology, 45(3), 791–816. https://doi.org/10.2307/3581
4. Korpimäki, E., & Norrdahl, K. (1991). "Numerical and functional responses of kestrels, short-eared owls, and long-eared owls to vole densities." Ecology, 72(3), 814–826. https://doi.org/10.2307/1940584