Competition–Colonization Trade-Off
The Competition–Colonization Trade-Off is the spatial metacommunity model developed by Robert Levins, David Culver, Alan Hastings, and David Tilman proving that an unlimited number of species can stably coexist on a single limiting habitat resource if superior local competitors have lower colonization rates, while predicting that habitat destruction selectively extinguishes the best competitor first.
The Competition–Colonization Trade-Off is a foundational mathematical model in spatial ecology, metacommunity theory, and conservation biology that explains how an arbitrary number of competing species can stably coexist within a physically homogeneous habitat by partitioning space dynamically through an inverse trade-off between local competitive dominance and regional dispersal (colonization) ability. First introduced in two-species patch-occupancy form by Robert Levins and David Culver in 1971, made rigorous with competitive displacement by Alan Hastings in 1980, and generalized to n-species communities and fragmented landscapes by David Tilman in 1994 (Ecology) and David Tilman, Robert M. May, Clarence L. Lehman, and Martin A. Nowak in 1994 (Nature), the model assumes a hierarchical competitive ranking (species 1 > species 2 > ... > species n) in which propagules of a superior competitor i instantaneously displace established adults of any inferior competitor j > i upon arriving at the same microsite. Despite this asymmetric local dominance and the absence of resource heterogeneity, inferior competitors persist regionally as "fugitive species" by colonizing empty patches vacated by mortality before superior competitors arrive. In 1994, Tilman and colleagues used this framework to uncover the extinction debt: when permanent habitat destruction (D) reduces the fraction of habitable patches across a landscape, the best competitor (species 1, with the lowest colonization rate c1 and highest spatial threshold 1 − m1/c1) is deterministically driven to extinction first.
Historical Development: From Fugitive Species to Spatial Patch Dynamics
Hutchinson's Fugitive Species and the Levins Metapopulation Framework
In 1951, the ecologist G. Evelyn Hutchinson introduced the concept of a "fugitive species"—an organism that is excluded in direct local competition by a dominant competitor, yet survives in nature by dispersing rapidly among newly disturbed or vacant habitat patches before superior competitors arrive. In 1969, Richard Levins formulated the single-species metapopulation differential equation governing the fraction of habitat patches p(t) ∈ [0, 1] occupied by a species across an infinite network of identical microsites:
dp / dt = c p(1 − p) − m p
where c > 0 is the per-patch colonization rate, 1 − p is the fraction of empty patches available for colonization, and m > 0 is the local patch extinction (mortality) rate. Setting dp/dt = 0 yields the single-species equilibrium occupancy p̂ = 1 − (m/c), which is positive and globally stable whenever c > m (compare Incidence Function Model).
Levins–Culver (1971), Hastings (1980), and Tilman (1994)
In 1971, Levins and Culver extended the metapopulation equation to two competing species, and in 1980 Alan Hastings formulated the hierarchical displacement model in which a superior competitor (species 1) can colonize patches occupied by an inferior competitor (species 2) and displace it, whereas species 2 can only colonize empty patches.
In 1994, David Tilman generalized Hastings's two-species equations to an n-species community competing for space in Minnesota cedar-savanna grasslands—where late-successional perennial bunchgrasses (Schizachyrium scoparium, Andropogon gerardii) invest heavily in dense root biomass that draws soil nitrogen down to low levels (Tilman's R* rule) at the expense of seed production, while weedy grasses (Agrostis scabra, Poa pratensis) allocate biomass to prolific small seeds at the expense of root competitive ability.
Mathematical Formulation of the n-Species Hierarchical Model
Hierarchical Competitive Ranking and Dynamical Equations
Consider a landscape composed of identical single-individual microsites (or local patches), in which a fraction D ∈ [0, 1) of sites has been permanently destroyed, leaving 1 − D habitable sites. The community consists of n species numbered in order of decreasing competitive ability:
Species 1 > Species 2 > Species 3 > ⋯ > Species i > ⋯ > Species n
This hierarchical ranking obeys three biological rules:
- Preemption Against Inferior Competitors: A propagule of inferior species i cannot establish in any site currently occupied by a superior competitor j < i (nor in a site already occupied by conspecifics of species i).
- Instantaneous Competitive Displacement: Whenever a propagule of superior species j < i arrives at a site occupied by inferior species i, species j displaces species i and takes over the site.
- Autonomous Dynamics of Higher Ranks: Because species 1 through i − 1 are unaffected by inferior species i, i + 1, ..., n, the dynamics of species 1 are identical to a single-species metapopulation, and the dynamics of species i depend only on species 1 through i.
Let pi(t) denote the fraction of total landscape sites occupied by species i at time t, ci its colonization rate, and mi its local mortality rate. The rate of change of pi is governed by Tilman's (1994) n-species patch-occupancy differential equation:
dpi/dt = ci pi (1 − D − ∑j=1i pj) − mi pi − ∑j=1i−1 cj pj pi
Each of the three terms on the right-hand side has a direct mechanistic meaning:
- Colonization Gain (ci pi (1 − D − ∑j=1i pj)): Species i successfully colonizes both truly empty habitable sites and sites currently occupied by inferior competitors (j > i), because inferior competitors offer no resistance to invasion by species i. Thus, the total fraction of sites open to colonization by species i is 1 − D − ∑j=1i pj.
- Intrinsic Mortality Loss (−mi pi): Occupied sites of species i become vacant through density-independent adult death or disturbance at rate mi.
- Competitive Displacement Loss (−∑j=1i−1 cj pj pi): Occupied sites of species i are invaded and usurped by arriving propagules of all superior competitors j ∈ {1, ..., i − 1} at combined rate ∑j=1i−1 cj pj.
Equilibrium Analysis and Multispecies Coexistence Thresholds
Equilibrium of the Best Competitor (Species 1)
For the top competitor (i = 1), the displacement summation is empty, reducing the equation to dp1/dt = c1p1(1 − D − p1) − m1p1. Setting dp1/dt = 0 yields the equilibrium occupancy p̂1 of the superior competitor:
p̂1 = 1 − D − m1/c1
As long as m1 > 0, even the best competitor in the community never occupies 100% of the habitable landscape (p̂1 < 1 − D). At equilibrium, species 1 always leaves an uncolonized vacant fraction of habitable space equal to m1/c1 because finite colonization rate c1 cannot instantaneously refill sites vacated by mortality m1.
Invasion Threshold and Equilibrium of Inferior Competitors
This persistent uncolonized gap m1/c1 provides the spatial niche opportunity for inferior competitor species 2. Substituting i = 2 and p1 = p̂1 = 1 − D − m1/c1 gives the per-capita invasion growth rate of species 2 when rare (p2 → 0):
(1 / p2)(dp2/dt)|p2→0 = c2(m1 / c1) − m2 − c1p̂1 > 0
Solving for the minimum colonization rate c2 required for species 2 to invade and stably coexist with species 1 yields the Tilman–Hastings threshold:
c2 > c1 (c1(1 − D) − m1 + m2) / m1
When species 1 and species 2 share the same mortality rate (m1 = m2 = m) in an intact landscape (D = 0), this condition simplifies to the quadratic colonization trade-off law:
c2 > c12/m (with c1 > m)
By mathematical induction, Tilman (1994) proved that for any community with uniform mortality mi = m and D = 0, the i-th ranked species has positive equilibrium occupancy p̂i > 0 if and only if its colonization rate satisfies ci > m / (1 − ∑j=1i−1 p̂j)2. Because the sum of equilibrium occupancies of any finite number of species remains less than 1, arbitrarily many species can stably coexist on a single limiting spatial resource provided inferior competitors have sufficiently superior colonization rates (cn > ... > c2 > c1).
| Species Rank & Ecological Strategy | Competitive Rank (i) | Root Allocation & Local Dominance (R*) | Seed Output & Colonization Rate (ci) | Minimum Threshold for Persistence (D = 0, mi = m) | Vulnerability to Habitat Destruction (D) |
|---|---|---|---|---|---|
| Apex Climax Competitor (Species 1) | 1st (Displaces all j > 1) | Highest root biomass; lowest soil resource R* | Lowest seed dispersal (c1 > m) | c1 > m / (1 − D) | Highest vulnerability: First species driven extinct at D1 = 1 − m1/c1 |
| Secondary Sub-Dominant (Species 2) | 2nd (Displaced by 1; displaces j > 2) | Intermediate-high root investment | Intermediate colonization (c2 > c12/m) | c2 > c12/m | Second in extinction sequence once species 1 declines |
| Mid-Successional Fugitive (Species i) | i-th in hierarchy | Moderate root investment | High colonization (ci ≫ c1) | ci > m / (1 − ∑j<i p̂j)2 | Temporarily increases when superior competitors decline, then goes extinct at higher D |
| Weedy Ruderal Colonizer (Species n) | n-th (Displaced by all j < n) | Minimal root allocation; highest R* | Highest seed production & dispersal (cn ≫ m) | Persists in the residual empty fraction 1 − ∑j<n p̂j | Lowest vulnerability: Last surviving species as D → 1 − m/cn |
Habitat Destruction and the Extinction Debt (Tilman et al., 1994)
Why the Best Competitor Goes Extinct First
In a landmark 1994 paper in Nature, David Tilman, Robert M. May, Clarence L. Lehman, and Martin A. Nowak analyzed the response of the n-species competition–colonization community to permanent habitat destruction D (compare habitat amount hypothesis). Inspecting the equilibrium occupancy of the best competitor (species 1),
p̂1(D) = 1 − D − m1/c1
reveals a counterintuitive conservation law:
- Because species 1 sits at the top of the competitive hierarchy, it invests in interference or resource-uptake structures at the expense of dispersal fecundity, giving it the lowest colonization rate c1 in the community and the highest spatial threshold 1 − m1/c1.
- Consequently, species 1 goes deterministically extinct (p̂1 → 0) as soon as permanent habitat destruction reaches the threshold D1 = 1 − m1/c1 = p̂1(0). For example, if a dominant tree or bunchgrass occupies 20% of an intact landscape (p̂1(0) = 0.20 because m1/c1 = 0.80), destroying 20% of the habitat (D = 0.20) drives the top competitor to regional extinction even though 80% of the landscape remains intact.
- As habitat destruction D increases progressively, species are eliminated in order from best competitor to worst competitor (species 1 first, then species 2, ..., leaving weedy fugitive species n last).
Time-Lagged Deterministic Extinction ("Extinction Debt") and Empirical Tests
Tilman et al. (1994) further demonstrated that when habitat destruction breaches 1 − mi/ci, patch-occupancy dynamics decay slowly over 50 to 400 generations before the final disappearance of the doomed competitive dominants. Tilman et al. named this generational overhang the extinction debt, showing that fragmented reserves can continue losing climax biodiversity for centuries after habitat loss halts. Subsequent empirical and theoretical studies (Levine & Rees, 2002; Calcagno et al., 2006) tested the seed size–number trade-off across vascular plants and integrated partial preemption with the spatial insurance hypothesis.
Frequently Asked Questions
How does the Competition–Colonization Trade-Off allow unlimited species to coexist on a single resource?
In a local patch, the superior competitor (species 1) always wins and displaces any inferior competitor. However, because local mortality (m1 > 0) continually kills individuals and colonization rate (c1 < ∞) is finite, species 1 can never occupy more than a fraction 1 − m1/c1 of the landscape at equilibrium. This permanently vacant fraction of patches (m1/c1) acts as an open spatial niche that an inferior competitor (species 2) can exploit provided its colonization rate satisfies c2 > c12/m. Repeating this mechanism down a hierarchy (cn > ⋯ > c2 > c1) allows an arbitrary number of species to coexist regionally.
Why does habitat destruction drive the best competitor extinct first?
Because the top competitor (species 1) invests energy in root or structural dominance rather than high seed dispersal, it has the lowest colonization rate c1 and occupies p̂1 = 1 − D − m1/c1. As habitat destruction D increases, p̂1 declines linearly until reaching zero at D = 1 − m1/c1. Paradoxically, as species 1 declines toward extinction, inferior weedy competitors are released from competitive displacement and can temporarily increase in abundance before going extinct at higher levels of D.
What is the "extinction debt" discovered by Tilman, May, Lehman, and Nowak (1994)?
The extinction debt is the deterministic future extinction of species caused by past habitat destruction (D ≥ 1 − m1/c1) that has not yet been realized because patch-occupancy dynamics decay slowly near the extinction threshold. Even after all further habitat destruction is halted, poor-colonizing dominant competitors whose minimum habitat threshold has been breached will continue declining to extinction over 50 to 400 generations.
References
- Levins, R., & Culver, D. (1971). "Regional coexistence of species and competition between rare species." Proceedings of the National Academy of Sciences, 68(6), 1246–1248. https://doi.org/10.1073/pnas.68.6.1246
- Hastings, A. (1980). "Disturbance, coexistence, history, and competition for space." Theoretical Population Biology, 18(3), 363–373. https://doi.org/10.1016/0040-5809(80)90059-390059-3)
- Tilman, D. (1994). "Competition and biodiversity in spatially structured habitats." Ecology, 75(1), 2–16. https://doi.org/10.2307/1939377
- Tilman, D., May, R. M., Lehman, C. L., & Nowak, M. A. (1994). "Habitat destruction and the extinction debt." Nature, 371(6492), 65–66. https://doi.org/10.1038/371065a0
- Calcagno, V., Mouquet, N., Jarne, P., & David, P. (2006). "Coexistence in a metacommunity: The competition–colonization trade-off is not dead." Ecology Letters, 9(8), 897–907. https://doi.org/10.1111/j.1461-0248.2006.00930.x