Ecology & Life Systems — Food Webs, Trophic Cascades and Population Dynamics

Why Mathematical Biology Works

Living systems look impossibly complex from the outside. A temperate forest contains thousands of interacting species across dozens of trophic levels, all embedded in fluctuating abiotic conditions. Yet remarkably simple mathematical models — systems of coupled ordinary differential equations — capture the core qualitative behaviour: the oscillating predator–prey cycles, the cascade that follows a keystone extinction, the convergence to a steady state when growth and mortality balance.

The reason is that ecology, like physics, is governed by conservation laws and rate processes. Energy enters as primary production, flows up trophic levels at roughly 10% efficiency per step, and eventually dissipates as heat. The mathematics describing these flows is the same mathematics used for chemical kinetics and fluid dynamics. This structural similarity is what makes ecological modelling powerful — and why simulations built on the same equations can reveal genuine ecological insight.

Predator–Prey Cycles: The Lotka-Volterra Equations

The foundational model of population ecology was published independently by Alfred Lotka in 1925 and Vito Volterra in 1926. Volterra was motivated by a puzzle in the Adriatic fishery data: during World War I, when fishing ceased, sharks increased and prey fish decreased — the opposite of what most people expected. His equations explained why.

The Lotka-Volterra model predicts neutrally stable oscillations: the populations cycle around a fixed point without converging or diverging. Real systems show damped oscillations because of carrying capacity and density-dependent effects, but the basic cycle — prey abundance drives predator growth, predator abundance drives prey decline, prey decline drives predator decline, which allows prey recovery — is robustly observed in nature. Classic data comes from the Hudson's Bay Company fur trading records of lynx and snowshoe hare populations from 1845 to 1935.

Trophic Cascades: When Wolves Change Rivers

A trophic cascade is an indirect effect that propagates through the food web from top predators downward (top-down) or from primary productivity upward (bottom-up). The most famous documented example is the reintroduction of grey wolves to Yellowstone National Park in 1995. The wolves reduced elk numbers and, more importantly, changed elk behaviour — elk stopped overgrazing stream banks for fear of ambush. Riparian vegetation recovered, riverbanks stabilised through root systems, and river channels narrowed and deepened. Beavers returned. Songbirds returned. The wolves, by changing a single top-predator–herbivore interaction, changed the physical geography of the park.

The mathematics of trophic cascades involves three-level Lotka-Volterra systems: plant–herbivore–carnivore. The equilibrium of the plant level can be controlled by the top predator through indirect suppression of the herbivore, a phenomenon ecologists call "mesopredator release" when a top predator is removed.

Circadian Rhythms: Biological Clocks

Almost all living organisms — from cyanobacteria to humans — maintain an internal clock that tracks the ~24-hour day-night cycle. These circadian rhythms regulate gene expression, hormone secretion, body temperature, feeding behaviour, and sleep–wake cycles. Disruption of circadian rhythms (shift work, jet lag, artificial light at night) correlates with elevated risk of metabolic syndrome, cancer, and cardiovascular disease.

The molecular mechanism was elucidated by Hall, Rosbash, and Young (Nobel Prize in Physiology or Medicine, 2017). The core is a transcription–translation feedback loop: the CLOCK/BMAL1 protein complex activates transcription of the PER and CRY genes; PER and CRY proteins accumulate, inhibit CLOCK/BMAL1, reduce their own transcription, degrade, and the cycle repeats. The mathematical description is a Goodwin oscillator — a negative feedback loop with a time delay that sustains oscillations.

The Goodwin oscillator is a prototypical negative-feedback oscillator: a gene product represses its own production through a chain of intermediates. Sustained oscillations require a Hill coefficient n ≥ 9 — a highly cooperative nonlinear repression. The simulation lets you explore how changing degradation rates shifts the period, and how light-pulse phase-response curves (PRCs) describe entrainment.

Brownian Motion and Diffusion in Living Systems

In 1827, botanist Robert Brown observed pollen grains moving erratically in water. In 1905, Einstein derived the diffusion coefficient D = k_B T / (6πηr) from first principles using the kinetic theory of heat — and showed that the mean-squared displacement grows linearly with time: ⟨r²⟩ = 6Dt for 3D diffusion, 4Dt for 2D, 2Dt for 1D. This provided the first quantitative evidence for molecular reality.

In living cells, Brownian motion is not just a curiosity — it is the primary transport mechanism at the nanoscale. Molecules smaller than ~5 nm diffuse rapidly enough to explore the entire cell volume in seconds. Larger complexes and organelles require active transport by molecular motors (kinesin, dynein) along cytoskeletal tracks. The crossover length between passive diffusion and active transport is set by the balance of thermal energy k_B T and the energy available per ATP molecule (~80 pN·nm).

Ecological Networks: From Simple to Complex

Real ecosystems are not simple predator–prey pairs. The food web of a temperate lake may contain hundreds of species connected by thousands of trophic links. Ecologists characterise these networks by their connectance C (fraction of possible links that exist), average chain length, and proportion of top, intermediate, and basal species.

Robert May's 1972 paradox showed mathematically that random networks become less stable (in the linear-stability sense) as they become larger and more connected — contradicting the then-prevailing intuition that diversity promotes stability. The resolution is that real food webs are not random: they have structure (weak links dominate, strong interactions are rare) that allows large, complex networks to be dynamically stable.

Energy Flow and Trophic Efficiency

The "10% rule" of ecology states that roughly 10% of the energy at one trophic level is available at the next. A primary producer captures solar radiation; a herbivore eating that plant captures about 10% as biomass; a carnivore eating the herbivore captures about 10% of that, meaning only ~1% of the original solar energy. This is why food chains are rarely longer than 5–6 links: beyond that, there is simply too little energy to support a viable predator population.

The reason for the inefficiency is thermodynamic: most of the energy in food is dissipated as heat during respiration, lost in faeces, or used for non-growth metabolic processes (maintenance costs). Ecological efficiency = assimilation efficiency × net production efficiency × consumption efficiency. For warm-blooded vertebrates, maintenance costs are very high (endothermy is expensive), reducing ecological efficiency to 1–5%. For ectotherms (insects, fish), efficiencies can reach 15–20%.

Algorithms & Methods