Module 7: Ungulate Herds — Collective Biophysics of the Savanna

The African savanna supports the largest mammalian migration on Earth—1.3–1.5 million wildebeest sweeping through the Serengeti-Mara ecosystem in an annual cycle. This module decomposes ungulate herd biophysics into collective motion (Reynolds boids, selfish-herd theory), anti-predator statistics (confusion, dilution, stotting signals), grass-height niche differentiation (Bell 1971 grazing succession), species-specific physiology (eland heterothermy, oryx rete mirabile, kudu acacia-tannin defence), and interspecies coupling across the savanna food web.

1. Wildebeest (Connochaetes) — Migration and Demography

Two species of wildebeest occupy different biogeographical zones: the blue wildebeest (Connochaetes taurinus), abundant across East and Southern Africa, and the black wildebeest (C. gnou), endemic to the highveld grasslands of South Africa. The Serengeti-Mara blue wildebeest population is estimated at 1.3–1.5 million individuals—a global biomass peak for a single large-herbivore species.

Rain-food regulation (Mduma 1999)

Mduma, Sinclair & Hilborn (1999) showed through 30 years of census data that wildebeest population dynamics in the Serengeti are bottom-up regulated by dry-season rainfall. The regulation mechanism runs through grass biomass. Let \(R\) = dry-season rainfall, \(G\) = grass biomass, \(N\) = wildebeest density:

\[\frac{dN}{dt} = r\,N\,\frac{G - G_c}{G + G_h} - \mu_0 N\]

with \(G \propto R\), giving rainfall-driven population equilibrium.

The migration loop: 3000 km per year

The wildebeest herd traces an annual loop approximately 3000 km long, tracking greenness through the ecosystem. December through March: southern Serengeti short-grass plains (calving). April–May: northwest movement through the Western Corridor. June–July: northern crossings of the Mara River. August–October: Masai Mara grazing. November–December: southern return. Each leg is punctuated by river crossings that exact enormous predation and drowning tolls.

Calving synchrony — predator swamping

Approximately 80% of wildebeest calves are born within a three-week window in late January / early February. This extraordinary synchrony is a classic example of predator swamping: by producing far more neonates than the local predator guild can consume during the vulnerable window, the herd ensures per-capita predation mortality crashes. If we denote predator saturation rate \(P_{\max}\) and neonate flux \(B(t)\):

\[p_{\text{kill}}(t) = \frac{\min(P_{\max},\; B(t))}{B(t)}\]

When \(B(t) \gg P_{\max}\), per-capita kill probability \(\to 0\).

The Mara River crossing — hydrodynamics and crocodiles

The Mara River crossings are defining events. Subalusky et al. (2017) estimated that approximately 6200 wildebeest carcasses per yeardrown or are killed crossing, delivering an enormous pulse of carbon, nitrogen, and phosphorus into the river ecosystem—a continent-scale trophic cascade. Physics of the crossing:

\[F_{\text{drag}} = \tfrac{1}{2} \rho_w C_d A U^2, \quad \text{resisting: } \mu_{\text{foot}}(M g - F_{\text{buoy}})\]

With river current \(U \sim 2\) m/s and peak flows higher still, adult wildebeest (250 kg) are often swept off their feet. Narrow hoofprint area and powerful neck extension allow swimmers to keep muzzles above water; juveniles suffer highest mortality.

Wildebeest locomotion efficiency

Wildebeest achieve minimum cost-of-transport at approximately 60 km/h gallop. Their narrow hoof (contact patch \(\sim 20\,\text{cm}^2\)) concentrates load on dry firm ground. Bone stress at mid-stance approaches the Alexander-predicted safety margin of 2–4: the tibia and metatarsals are loaded near their fatigue limit, consistent with a species that covers 3000 km/yr.

2. Herd Dynamics — Collective Motion Physics

Ungulate herds exhibit the same three classes of self-organized behavior as starling flocks and fish schools: alignment, cohesion, and separation. Reynolds (1987) showed these three local rules alone produce realistic group motion. Krause & Ruxton (2002) applied fish-school physics to ungulates: each rule has a characteristic distance scale.

Reynolds boid model

\[\vec v_i(t+\Delta t) = \vec v_i + w_s \vec F_s + w_a \vec F_a + w_c \vec F_c + w_p \vec F_p\]

\(\vec F_s\): separation; \(\vec F_a\): alignment within neighborhood; \(\vec F_c\): cohesion to centroid; \(\vec F_p\): predator repulsion.

Confusion effect (Treisman 1975)

Predators attacking a herd experience a confusion effect: the visual system cannot lock onto a single target when surrounded by many near-identical moving objects. Attention is modeled as a fixed-capacity channel:

\[p_{\text{lock}} = \frac{1}{1 + k N_{\text{visible}}}\]

Dilution effect

For a predator that captures exactly one prey per attack, each individual in a group of \(N\) has capture probability \(1/N\):

\[p_{\text{capture}}(i) = \frac{1}{N}\]

Selfish herd (Hamilton 1971)

Hamilton’s selfish-herd theory states that each individual minimizes its domain of danger—the Voronoi region around it, within which it is the nearest prey to any predator entering. The optimal strategy is movement toward the centroid, producing spontaneous aggregation without explicit cohesion.

\[\text{Domain of danger}_i = |\text{Voronoi cell}(i)|\]

Migration as a traveling wave

The spatial advance of the wildebeest front across the Serengeti can be modeled as a Fisher–KPP traveling wave with invasion speed

\[c^* = 2 \sqrt{r D}\]

With \(r \sim 0.02\)/day (forage uptake rate) and \(D \sim 1\,\text{km}^2/\text{day}\) (diffusive mixing), the wave advances at \(\sim 0.3\) km/day—consistent with the observed seasonal front.

Reynolds Boid Rules

Separationavoid crowdingAlignmentmatch neighbors' headingCohesionsteer toward centroid

3. Zebra (Equus quagga, E. grevyi) — The Stripe Problem

Two zebra species occupy the African savanna. The plains zebra (Equus quagga) is abundant across East and Southern Africa, while Grevy’s zebra (E. grevyi) is the larger, narrower-striped northern species. Why do zebras have stripes? Four hypotheses have been tested:

Hypothesis 1: crypsis — largely disproven

How & Zanker (2014) measured zebra spectral reflectance and modeled lion visual system detection. At 10 m range under African daylight, lions readily detect zebras. Crypsis cannot account for stripes.

Hypothesis 2: social cohesion

Caro (2014) found weak evidence that striped coats aid kin and conspecific recognition in groups. Statistical signal exists but is small.

Hypothesis 3: tabanid-fly deterrence — LEADING

Caro et al. (2019) ran painted-horse experiments in which horses were dressed in black, white, and striped coats and blood-feeding tabanid flies were counted. Landings were reduced by approximately 80% on striped horses. The proposed mechanism is optical-flow disruption at the scale of fly vision: stripe-period matches the insect’s motion-detecting receptive field, producing aliasing during final approach.

\[\lambda_{\text{stripe}} \approx \frac{v_{\text{approach}}}{f_{\text{flicker}}} \approx \frac{1\,\text{m/s}}{50\,\text{Hz}} = 2\,\text{cm}\]

Close to the observed zebra stripe period (\(\sim 3\)–5 cm at adult body scale).

Hypothesis 4: thermoregulation — disproven but nuanced

Larison et al. (2021) tested the hypothesis that alternating black/white stripes set up micro-convection currents that cool the animal. High-resolution thermal imaging showed no meaningful cooling effect. However the temperature difference between black and white stripes (~15 °C at noon) does set up small-scale convection that may aid pelage turnover.

Polarized-reflectance optics

How & Zanker (2014) also examined polarization of reflected light. Tabanid flies preferentially orient to strongly polarized reflections from smooth dark surfaces (often water). Zebra stripes break up the polarization signature, reducing attractiveness:

\[I_{\text{pol}}(\theta) = \frac{I_\parallel - I_\perp}{I_\parallel + I_\perp}\]

4. Buffalo (Syncerus caffer) — Herd Voting and Disease

The African buffalo is a 600–900 kg bovid numbering approximately 800,000 in Africa. It has a fearsome reputation among human hunters (200+ deaths per year) and is a keystone mover of both vegetation structure and pathogen dynamics.

Group “voting” behavior (Prins 1996)

Prins (1996) demonstrated in Manyara buffalo that when the herd rises from a rest bout to graze, the travel direction is determined by a consensus vote. Adult females stand up, face in a chosen direction, and the herd travel-direction is the vector sum of these female orientations. Formally:

\[\vec d_{\text{herd}} = \frac{1}{N_{\text{vote}}} \sum_i \hat n_i\]

Only individuals in an “engaged” posture cast votes; a quorum threshold triggers movement.

Bovine tuberculosis dynamics

Michel et al. (2007) documented bovine tuberculosis (Mycobacterium bovis) in Kruger Park buffalo herds. The pathogen transmits through droplet aerosol during close contact, particularly at shared waterholes. The SIR-like model with herd-size dependence is:

\[R_0 = \frac{\beta N}{\gamma + \mu}\]

Large contact rate \(\beta\) and long infectious duration \(1/\gamma\) produce \(R_0 > 1\) in large herds, ensuring endemic persistence.

5. Bell (1971) Grazing Succession — Niche Differentiation by Body Size

One of the most elegant hypotheses in African ecology, due to Bell (1971), explains the coexistence of multiple ungulate species on shared grassland by body-size-differentiated grass-height niches. Large herbivores process coarse, tall grass profitably because of capacious rumen volume; small herbivores cannot, and require young, tender, protein-rich short grass.

\[\text{Retention time} \propto M^{0.25}, \quad \text{Intake rate} \propto M^{0.75}\]

Bigger animals can tolerate longer gut retention, extracting nutrition from lower-quality forage.

Succession order

  1. Zebra first: crops tall coarse grass stems; body mass ~300 kg.
  2. Wildebeest second: graze shortened intermediate-height grass; body mass ~250 kg.
  3. Thomson gazelle last: takes the remaining tender, short protein-rich regrowth; body mass ~25 kg.

Facilitation vs competition

In the Serengeti the three species show grazing facilitation: zebra clearing tall stems allows wildebeest to access the medium layer, and wildebeest grazing stimulates tender regrowth favorable for gazelle. The coupling is formalized by cross-diffusion in the grass height compartments (simulated below).

6. Thomson Gazelle, Eland, Kudu, Oryx — Specialist Physiology

Thomson’s gazelle — stotting

Eudorcas thomsonii (25 kg) produces the distinctive stotting (or pronking) behavior: stiff-legged vertical leaps to 60–80 cm even when pursued by predators. FitzGibbon (1988) showed stotting is an honest signal of condition to cheetahs: only well-conditioned gazelles can afford the energetic cost, and cheetahs preferentially abandon pursuit of strong-stotting individuals. The energy cost of a single stott is

\[E_{\text{stott}} \approx M g h + \tfrac{1}{2} M v_{\text{takeoff}}^2\]

For \(M = 25\) kg and \(h = 0.7\) m, \(E \approx 170\) J per leap—sustainable only with good body condition.

Eland — heterothermy

The eland (Taurotragus oryx) is the largest antelope, 500–900 kg, and extreme drought adaptation rests on heterothermy. Taylor (1969) showed elands store body heat during the day, allowing core temperature to rise from 39 to 42 °C, and dump that heat by non-evaporative cooling overnight. This strategy conserves water; eland have the lowest water turnover rate of any savanna bovid:

\[Q_{\text{stored}} = m c_p \Delta T_{\text{body}}, \quad m c_p \Delta T \approx 500 \cdot 3500 \cdot 3 = 5.2\,\text{MJ}\]

Equivalent to evaporating 2.1 L of water—but without drinking it.

Oryx — carotid rete mirabile

Oryx (Oryx gazella) possess a carotid rete mirabile—a heat exchanger in which warm arterial blood to the brain is cooled against cool venous blood returning from the evaporative nasal cavity. Brain temperature is held 2–3 °C below body core even as the body heats to 45 °C. Oryx also recycle urea into rumen microbes, recapturing nitrogen without excreting water. This combination (cross-referenced with savanna-adaptation biochemistry in climate-biodiversity M13) lets oryx survive indefinitely without drinking.

Kudu — tannin-induced mortality

In 1982 van Hoven documented a mass-mortality event of greater kudu (Tragelaphus strepsiceros) in South African game ranches. Prolonged browsing on fenced acacia triggered the acacias to raise leaf condensed-tannin concentrations, apparently via an airborne ethylene signal transmitted from heavily browsed trees. Tannin concentrations above 6% of dry weight are lethal to kudu, which lack a sufficient proline-rich salivary protein defence. The cascade illustrates plant-herbivore chemical warfare at the ecosystem scale.

7. Giraffe & Cross-Module Couplings

Although the giraffe is treated in detail in Module 3, its migratory ecology couples strongly to ungulate herd dynamics. Giraffes browse the canopy stratum and remove up to 15% of annual acacia leaf production, indirectly opening understory for grazers. Giraffe abundance affects grass height distributions through trophic top-down control on woody encroachment, which in turn modulates Bell grazing succession described in Section 5.

Migration as a traveling wave — redux

Treating the wildebeest migration front as a Fisher-KPP reaction-diffusion wave:

\[\frac{\partial u}{\partial t} = D \nabla^2 u + r u \left(1 - \frac{u}{K(\vec x, t)}\right)\]

The carrying capacity \(K(\vec x, t)\) pulses seasonally with rainfall, producing the observed annual loop rather than a stationary advance.

Coupled Lotka-Volterra food web

The savanna herbivore guild is coupled to lions (Module 5), hyenas, cheetahs, and leopards (Module 6) through a generalized Lotka-Volterra system. Stability is preserved because apex predators exert top-down control that prevents any single ungulate from monopolizing forage:

\[\frac{dN_i}{dt} = r_i N_i \left(1 - \frac{\sum_j a_{ij} N_j}{K_i}\right) - \sum_k \phi_{ik} P_k \frac{N_i}{\sum_l N_l + h}\]

8. Self-Organization & the Physics of Collective Motion

Ungulate herds, starling murmurations, and fish schools share a common statistical mechanics. At the level of coarse-grained variables, the herd field behaves like a self-propelled fluid described by the Toner-Tu continuum theory:

\[\partial_t \vec v + \lambda (\vec v \cdot \nabla) \vec v = -\nabla P + D \nabla^2 \vec v + \alpha \vec v - \beta |\vec v|^2 \vec v + \vec \eta\]

Toner-Tu (1995) gives phase transitions between disordered and aligned flocking states, with a true-long-range-order phase at finite temperature—a purely non-equilibrium effect.

Scale-free correlations in flocks

Cavagna et al. (2010) measured velocity correlations in starling flocks and showed that \(C_{\parallel}(r)\) decays as a scale-free power law with correlation length equal to the flock diameter. This means collective response to a perturbation (a hawk attack, an environmental cue) propagates instantly across the flock. The same phenomenon almost certainly occurs in wildebeest herds but has not been quantitatively verified at the one-million-animal scale.

Information percolation

How does a Mara crocodile alarm propagate across 10,000 wildebeest? Each animal couples to its nearest 6–7 neighbors (Ballerini et al. 2008 topological interaction). Let reaction latency \(\tau_r \sim 100\) ms and inter-animal spacing \(\ell \sim 2\) m. Information wave speed

\[v_{\text{info}} \approx \ell / \tau_r \approx 20\,\text{m/s}\]

Much faster than the wildebeest running speed—the herd can preemptively shift before any individual could outrun the predator.

Phase transition at critical density

Vicsek et al. (1995) showed that self-propelled agents with simple alignment rule undergo a disorder-order phase transition at a critical noise level. In herd biology this corresponds to the dichotomy between scattered grazing (disordered) and migration (ordered).

\[|\langle \vec v \rangle| \sim (\eta_c - \eta)^\beta, \qquad \beta \approx 0.45\]

Simulation 1: Boid Herd with Predator — Confusion & Dilution

A Reynolds boid model augmented with predator pursuit. Each trial runs 400 time-steps with the herd evading a predator that pursues its nearest prey. Herd size is swept over \(N \in \{5, 10, 20, 40, 80, 150\}\) and 30 replicates per size. The simulation reports predator success rate (testing the dilution effect), a bearing-angle confusion metric, and produces a final snapshot visualization.

Python
script.py164 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Simulation 2: Bell Grazing Succession — Three-Compartment Grass Dynamics

A compartmental grass-biomass model coupled to three grazer guilds: zebra (prefers tall), wildebeest (medium), Thomson gazelle (short). Rainfall drives growth; senescence and selective removal cascade the grass through height classes. Each grazer population follows a Ricker-logistic dynamic tied to available forage. The output shows the characteristic succession pattern predicted by Bell (1971) and observed in the Serengeti.

Python
script.py139 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Key References

• Mduma, S. A. R., Sinclair, A. R. E., & Hilborn, R. (1999). “Food regulates the Serengeti wildebeest.” Journal of Animal Ecology, 68, 1101–1122.

• Subalusky, A. L. et al. (2017). “Annual mass drownings of the Serengeti wildebeest migration influence nutrient cycling in the Mara River.” PNAS, 114, 7647–7652.

• Reynolds, C. W. (1987). “Flocks, herds and schools: a distributed behavioral model.” Computer Graphics, 21, 25–34.

• Krause, J. & Ruxton, G. D. (2002). Living in Groups. Oxford University Press.

• Treisman, M. (1975). “Predation and the evolution of gregariousness.” Animal Behaviour, 23, 779–800.

• Hamilton, W. D. (1971). “Geometry for the selfish herd.” Journal of Theoretical Biology, 31, 295–311.

• How, M. J. & Zanker, J. M. (2014). “Motion camouflage induced by zebra stripes.” Zoology, 117, 163–170.

• Caro, T. et al. (2019). “Benefits of zebra stripes: behaviour of tabanid flies around zebras and horses.” PLoS ONE, 14, e0210831.

• Larison, B. et al. (2021). “How the zebra got its stripes: a problem with too many solutions.” Royal Society Open Science, 2, 140452.

• Prins, H. H. T. (1996). Ecology and Behaviour of the African Buffalo. Chapman & Hall.

• Michel, A. L. et al. (2007). “Wildlife tuberculosis in South African conservation areas.” Veterinary Microbiology, 112, 91–100.

• Bell, R. H. V. (1971). “A grazing ecosystem in the Serengeti.” Scientific American, 225, 86–93.

• FitzGibbon, C. D. (1988). “Stotting in Thomson’s gazelles: an honest signal of condition.” Behavioral Ecology and Sociobiology, 23, 69–74.

• Taylor, C. R. (1969). “The eland and the oryx.” Scientific American, 220, 88–95.

• van Hoven, W. (1982). “Mortalities in kudu populations related to chemical defence in trees.” Journal of African Zoology, 105, 141–145.