Macroecology and the Latitudinal Diversity Gradient

Macroecology and the latitudinal diversity gradient (LDG) investigate the statistical patterns and evolutionary mechanisms governing the distribution of life across continental to global spatial scales and evolutionary timescales. Coined in a 1989 Science article by James H. Brown and Brian A. Maurer, macroecology complements reductionist small-plot field experiments by taking a "top-down" statistical-mechanics approach to large ensembles of species—analyzing universal emergent laws in Species Abundance Distributions (SADs), Species–Area Relationships (SARs), body-size allometries, and geographic range geometries.

The most pervasive macroecological pattern on Earth is the Latitudinal Diversity Gradient (LDG): the systematic increase in taxonomic, phylogenetic, and functional biodiversity from the polar regions toward the equatorial tropics across nearly all terrestrial, freshwater, and marine clades (including vascular plants, insects, birds, mammals, amphibians, reef fishes, and planktonic foraminifera). First documented by pioneering naturalists Alexander von Humboldt (1807), Alfred Russel Wallace (1878), and Charles Darwin, the LDG reflects the coupled interplay of solar energy and water budgets, tropical biome area, geological climate stability, metabolic mutation rates, and cradles versus museums of speciation.

Fisher's Log-Series, Preston's Log-Normal, and the Veil Line

Across every natural community ever censused, a few species are hyper-abundant while the vast majority of species are rare, generating a hollow-curve distribution of individuals per species. Three canonical statistical models describe this Species Abundance Distribution (SAD):

1. Fisher's Log-Series (1943): Derived by Ronald A. Fisher, A. Steven Corbet, and C. B. Williams to describe Malayan butterfly and Rothamsted moth light-trap catches, the log-series predicts the number of species containing n individuals as S(n) = α xⁿ / n (where 0 < x < 1 is a sample-size constant and Fisher's α is a sample-size-independent diversity parameter given by S = α ln(1 + N/α)). In Fisher's log-series, the modal class is always singletons (n = 1).
2. Preston's Canonical Log-Normal and Veil Line (1948, 1962): Engineer and ecologist Frank W. Preston binned species abundances into base-2 logarithmic "octaves" (1–2, 2–4, 4–8, 8–16 individuals) and discovered a bell-shaped Gaussian normal curve: S(R) = S₀ exp(−*a*²R²), where R is octave distance from the modal octave S₀. Preston explained why small samples appear to follow Fisher's monotonic hollow curve: in incomplete censuses, the left tail of rare species (those with <1 expected individual in the sample) is hidden behind a sampling cutoff he termed the veil line. Doubling sampling effort shifts the veil line one octave to the left, unveiling the true interior mode.
3. Hubbell's Unified Neutral Theory of Biodiversity and Biogeography (UNTBB, 2001): Formulated by Stephen P. Hubbell, UNTBB makes the radical null assumption of ecological equivalence among all individuals of all trophically similar species (equal per-capita birth, death, dispersal m, and point speciation rate ν in a zero-sum community of size J_M). Driven solely by stochastic demographic drift ("ecological drift"), dispersal limitation, and speciation governed by the fundamental biodiversity number θ = 2*J*_Mν, UNTBB predicts a zero-sum multinomial (ZSM) abundance distribution that fits tropical forest plot censuses (such as Barro Colorado Island, Panama) with remarkable precision.

Rosenzweig's Triphasic SAR across Spatial Scales

Synthesizing global biogeographic datasets in Species Diversity in Space and Time (1995), Michael L. Rosenzweig demonstrated that when log S is plotted against log A from 1 m² quadrats up to the entire surface of Earth (5.1 × 10¹⁴ m²), the Species–Area Relationship (S = cA^z) follows a triphasic sigmoidal curve:

• Phase 1: Intra-Provincial / Local Scale (z ≈ 0.10–0.18): Sampling nested subplots within a single continuous biotic province adds habitat heterogeneity, while transient individuals immigrating from surrounding habitat (spatial mass / rescue effects) inflate richness in small plots and flatten the slope z.
• Phase 2: Archipelagic / Isolate Scale (z ≈ 0.25–0.35): Across true oceanic islands or isolated continental mountaintops, MacArthur–Wilson immigration–extinction dynamics without continuous rescue steepen the exponent around Preston's canonical z ≈ 0.26.
• Phase 3: Inter-Provincial / Evolutionary Scale (z ≈ 0.65–1.0): Comparing entire continents and biogeographic realms separated by oceanic and tectonic barriers where independent in situ cladogenesis replaces dispersal as the source of diversity yields a near-unity exponent.

Rapoport's Rule, Bergmann's Rule, and Janzen's Mountain Passes

Macroecologists have established several classical ecogeographical regularities across latitude and elevation:

• Rapoport's Rule and Janzen's "Mountain Passes are Higher in the Tropics" (1967): Formulated by Eduardo H. Rapoport (1975) and named by George C. Stevens (1989), Rapoport's rule observes that the latitudinal and elevational geographic range sizes of species are narrowest in the tropics and broadest at high temperate/boreal latitudes. Earlier, Daniel H. Janzen (1967) provided the physiological mechanism: because a tropical lowland organism experiences minimal seasonal temperature variation across the year, it evolves a narrow thermal tolerance breadth (stenothermy). Consequently, a 2,000-meter tropical mountain pass acts as an insurmountable physiological thermal barrier to gene flow—promoting allopatric speciation—whereas a temperate organism adapted to a 40 °C annual winter-to-summer temperature swing easily crosses a 2,000-meter temperate mountain range.
• Bergmann's Rule (1847) and Allen's Rule (1877): Within broadly distributed endothermic clades (mammals and birds), populations in colder high-latitude climates exhibit larger body mass (Bergmann's rule, reducing surface-area-to-volume ratio L²/L³) and shorter protruding extremities such as pinnae, bills, and limbs (Allen's rule)—as seen across the progression from the desert fennec fox (Vulpes zerda) to the red fox (V. vulpes) and Arctic fox (V. lagopus).

Geographic Area, Climate Stability, and Metabolic Kinetics

Over 30 hypotheses have been proposed to explain why a single hectare of Ecuadorian Amazonian rainforest at Yasuní harbors ~655 tree species—exceeding the native tree richness of all of North America north of Mexico. Modern synthesis groups these into four non-mutually-exclusive causal frameworks:

1. Geographic Area and the Mid-Domain Effect: John Terborgh (1973) and Michael L. Rosenzweig (1995) noted that because Earth is a sphere, the tropical belt (23.5° N to 23.5° S) encompasses the largest land and ocean surface area of any climatic zone, and is unique in Having North and South tropical zones contiguous across the equator. Larger biome area supports larger population sizes (N), lowering stochastic extinction rates (E) while expanding the geographic target for allopatric speciation.
2. Historical Time-for-Speciation and Tropical Niche Conservatism: Formulated by John J. Wiens and Michael J. Donoghue (2004), the Tropical Niche Conservatism Hypothesis notes that most major angiosperm, vertebrate, and insect clades originated during the warm Paleogene/Eocene greenhouse when tropical environments covered much of the globe. Because evolving physiological tolerance to freezing cellular water is evolutionarily difficult, relatively few lineages breached the temperate frost line, and Pleistocene glaciations repeatedly scoured high latitudes while tropical rainforest refugia accumulated lineages over tens of millions of years.
3. Energy–Water Productivity and Metabolic Kinetics: David J. Currie (1991) and Bradford A. Hawkins (2003) showed that actual evapotranspiration (AET) and solar radiation explain >70% of spatial variance in continental richness. Under the Metabolic Theory of Ecology (Shane D. Wright and James H. Brown), elevated tropical ambient temperatures accelerate ectothermic metabolic rates (e^−*E*/kT), shorten generation times, increase reactive oxygen species DNA mutation rates, and speed up molecular evolution and speciation ("evolutionary speed hypothesis").
4. Biotic Interactions (Dobzhansky–MacArthur–Janzen–Connell): Theodosius Dobzhansky (1950) proposed that while abiotic mortality dominates at the poles, year-round tropical warmth intensifies biotic selection—parasitism, predation, pollination specialization, and Janzen–Connell negative density dependence—driving finer niche partitioning and coevolutionary diversification.

Jablonski's "Out of the Tropics" Cradle-and-Museum Model

For decades, evolutionary biologists debated whether the tropics are a "cradle" of diversity (generating higher origination/speciation rates λ) or a "museum" of diversity (preserving older lineages through lower extinction rates μ). Analyzing a comprehensive paleontological database of Cenozoic marine bivalve genera across latitude over the past 11 million years, David Jablonski, Kaustuv Roy, and James W. Valentine (2006) proved that the tropics are both a cradle and a museum: the Out of the Tropics (OTT) model shows that most marine lineages originate in the tropics (high λ), persist there with low extinction (low μ), and simultaneously expand their geographic ranges poleward into extratropical latitudes without losing their tropical presence.

What clades are exceptions to the Latitudinal Diversity Gradient?

Although >85% of studied taxonomic groups peak in richness near the equator, several instructive inverse latitudinal diversity gradients exist. Ichneumonid parasitoid wasps (Ichneumonidae, studied by Daniel H. Janzen and Ian D. Gauld), aphids (Aphidoidea), sawflies (Symphyta), penguins (Spheniscidae), pinnipeds (seals and walruses), ectomycorrhizal fungi, and marine kelps (Laminariales) reach peak species richness in temperate or subpolar latitudes. In ichneumonids and aphids, high tropical host diversity and plant chemical defenses fragment host populations to low local densities (resource fragmentation hypothesis), making it difficult for host-specialist temperate guilds to locate hosts in tropical forests.

What is Rapoport's rescue effect and the Mid-Domain Effect?

Proposed by Robert K. Colwell and David C. Lees (2000), the Mid-Domain Effect (MDE) is a geometric null model demonstrating that if species' geographic ranges of varying sizes are shuffled randomly within a bounded domain (such as between the North and South Poles, or between sea level and a mountain summit) without any environmental gradients, range overlaps inevitably peak in the geometric center of the domain (the equator or mid-elevation), analogous to shaking pencils of different lengths inside a pencil box. While MDE contributes to broad-range overlap geometry, empirical residuals confirm that climate, energy, and evolutionary history remain the primary drivers of the LDG.

How does marine biodiversity differ in its latitudinal and longitudinal distribution from terrestrial biodiversity?

Whereas terrestrial biodiversity peaks across three equatorial continental blocks (the Neotropical Amazon/Andes, Afrotropical Congo, and Indo-Malayan rainforests), shallow coastal marine biodiversity exhibits both a strong latitudinal gradient and a pronounced longitudinal hotspot centered in the Coral Triangle (Indonesia, Philippines, Papua New Guinea, and Solomon Islands), where >600 reef-building scleractinian coral species and >2,000 reef fish species coexist. By contrast, pelagic open-ocean taxa (such as euphausiid krill, foraminifera, and tuna) often display a bimodal latitudinal gradient peaking at subtropical/temperate latitudes (20°–40° N and S) where oceanic thermocline dynamics balance nutrient upwelling and solar irradiation.

1. Janzen, D. H. (1967). "Why Mountain Passes are Higher in the Tropics." The American Naturalist, 101(919), 233–249. doi:10.1086/282487
2. Brown, J. H., & Maurer, B. A. (1989). "Macroecology: The Division of Food and Space Among Species on Continents." Science, 243(4895), 1145–1150. doi:10.1126/science.243.4895.1145
3. Hubbell, S. P. (2001). The Unified Neutral Theory of Biodiversity and Biogeography. Monographs in Population Biology, Vol. 32. Princeton University Press. https://press.princeton.edu/books/paperback/9780691021287/the-unified-neutral-theory-of-biodiversity-and-biogeography
4. Jablonski, D., Roy, K., & Valentine, J. W. (2006). "Out of the Tropics: Evolutionary Dynamics of the Latitudinal Diversity Gradient." Science, 314(5796), 102–106. doi:10.1126/science.1130880