Biodiversity-Stability Relationships
Insurance hypothesis, functional diversity metrics, metabolic scaling, species-energy theory, and biodiversity-ecosystem function experiments
Does biodiversity matter for ecosystem function? This fundamental question has driven decades of ecological research, from Tilman’s Cedar Creek grassland experiments to the metabolic theory of ecology. The answer is unequivocally yes — more diverse ecosystems are more productive, more stable, and more resilient to perturbation. The mechanisms involve biochemical complementarity, where different species exploit distinct chemical niches, and the portfolio effect, where variance is reduced through the averaging of independent biochemical responses to environmental fluctuation.
6.1 Insurance Hypothesis
Yachi & Loreau (1999) formalized the insurance hypothesis: biodiversity provides a buffer against environmental fluctuations because different species respond differently to the same perturbation. If species A’s photosynthetic biochemistry collapses under drought (e.g., C3 grasses with low water-use efficiency), species B (a C4 grass with efficient PEP carboxylase) compensates, maintaining total ecosystem productivity.
The key insight is that biochemical redundancy at the ecosystem level is not wasteful but essential. Each species carries a unique set of enzyme isoforms, stress-response pathways, and metabolic strategies. Under stable conditions, this redundancy appears unnecessary. Under perturbation, it becomes the difference between ecosystem collapse and persistence.
The Portfolio Effect: Formal Derivation
Consider an ecosystem with \(S\) species, each contributing to total ecosystem function \(F = \sum_{i=1}^{S} f_i\), where \(f_i\) is the contribution of species \(i\). The variance of total ecosystem function over time is:
Variance of Ecosystem Function
\[ \sigma^2_{\text{eco}} = \text{Var}\left(\sum_{i=1}^{S} f_i\right) = \sum_{i=1}^{S} \sigma_i^2 + \sum_{i \neq j} \text{cov}(f_i, f_j) \]
Now assume all species have equal variance \(\sigma^2\) and equal pairwise covariance \(\rho \sigma^2\), where \(\rho\) is the average correlation. There are \(S\) variance terms and \(S(S-1)\)covariance terms:
\[ \sigma^2_{\text{eco}} = S \sigma^2 + S(S-1) \rho \sigma^2 = S\sigma^2 \left[1 + (S-1)\rho\right] \]
The per-species variance of total function, normalized by \(S^2\)(since mean function scales as \(S\)):
Coefficient of Variation Squared
\[ \text{CV}^2 = \frac{\sigma^2_{\text{eco}}}{\mu_{\text{eco}}^2} = \frac{S\sigma^2[1 + (S-1)\rho]}{S^2 \mu^2} = \frac{\sigma^2}{S\mu^2}[1 + (S-1)\rho] \]
When species respond independently (\(\rho = 0\)), we get the classic portfolio result: \(\text{CV}^2 = \sigma^2 / (S\mu^2)\), which decreases as \(1/S\). When species are perfectly positively correlated (\(\rho = 1\)), CV is independent of \(S\) — no insurance benefit. The insurance hypothesis operates whenever \(\rho < 1\), meaning species have at least partially independent biochemical responses to environmental change.
Biochemical Basis of Response Diversity
What determines \(\rho\)? The correlation depends on the degree to which species share the same biochemical vulnerability. Consider drought stress:
- C3 grasses (e.g., Festuca): RuBisCO-only carbon fixation, high photorespiration under heat/drought, stomata close early
- C4 grasses (e.g., Andropogon): PEP carboxylase concentrates CO₂ around RuBisCO, maintaining photosynthesis under drought
- CAM plants (e.g., Opuntia): nighttime CO₂ fixation via malic acid storage, extreme drought tolerance
- Legumes (e.g., Trifolium): N₂ fixation via nitrogenase in root nodules, independent nitrogen supply
A community containing all four functional types has low \(\rho\) because drought decimates C3 grasses but barely affects CAM plants. The portfolio effect is maximized when species span the widest possible range of biochemical strategies.
Yachi-Loreau Insurance Result (1999)
\[ \text{Stability} = \frac{\mu_{\text{eco}}}{\sigma_{\text{eco}}} = \frac{S\mu}{\sigma\sqrt{S[1 + (S-1)\rho]}} \]
Stability increases with \(S\) whenever \(\rho < 1\), and increases fastest when species have independent (\(\rho \approx 0\)) or negatively correlated (\(\rho < 0\)) responses.
6.2 Functional Diversity Metrics
Taxonomic diversity (species counts) may not capture the biochemical diversity that matters for ecosystem function. Two communities with the same number of species can differ enormously in functional diversity — the range and distribution of biochemical traits such as photosynthetic pathway, nitrogen fixation capacity, secondary metabolite production, mycorrhizal association type, and leaf mass per area.
Rao’s Quadratic Entropy
Rao (1982) introduced a metric that combines species abundances with functional distances. Given \(S\) species with relative abundances \(p_1, \ldots, p_S\)and a matrix of functional distances \(d_{ij}\) between species\(i\) and \(j\):
Rao’s Quadratic Entropy
\[ Q = \sum_{i=1}^{S} \sum_{j=1}^{S} d_{ij} \cdot p_i \cdot p_j \]
where \(d_{ij}\) can be computed from trait dissimilarity (e.g., Gower distance across multiple biochemical traits) and \(\sum p_i = 1\).
This can be expanded as:
\[ Q = \sum_{i < j} 2 d_{ij} p_i p_j + \sum_{i=1}^{S} d_{ii} p_i^2 \]
Since \(d_{ii} = 0\), this simplifies to \(Q = 2\sum_{i < j} d_{ij} p_i p_j\).
Components of Functional Diversity
Functional diversity decomposes into three independent components (Villeger et al. 2008):
- Functional Richness (FRic): Volume of trait space occupied. Measured as the convex hull volume in multivariate trait space. A high FRic means the community spans a wide range of biochemical strategies (e.g., from shade-tolerant understory herbs to canopy trees with sun-adapted photosynthesis).
- Functional Evenness (FEve): How evenly abundances are distributed across trait space. If one biochemical strategy dominates, FEve is low even if FRic is high. Computed from the minimum spanning tree in trait space:\[ \text{FEve} = \frac{\sum_{b=1}^{S-1} \min\left(\text{PEW}_b, \frac{1}{S-1}\right) - \frac{1}{S-1}}{1 - \frac{1}{S-1}} \]where PEW\(_b\) is the partial weighted evenness of branch \(b\).
- Functional Divergence (FDiv): How abundances are distributed relative to the center of trait space. High FDiv indicates that the most abundant species have extreme trait values — the community is dominated by biochemical specialists rather than generalists.
Taxonomic vs. Functional Diversity Relationship
The relationship between taxonomic diversity (\(S\)) and functional diversity (\(Q\)) is not linear. When new species are functionally redundant with existing ones, adding species increases \(S\) without increasing \(Q\). Formally, adding species\(k\) increases function when:
Condition for Functional Gain
\[ \Delta Q = 2 p_k \sum_{i=1}^{S} d_{ik} p_i > 0 \]
This is always positive when \(d_{ik} > 0\) for at least one existing species — any species that is biochemically distinct from the current community adds functional diversity.
However, the marginal gain decreases as the community fills trait space. The relationship is typically saturating: \(Q \approx Q_{\max}(1 - e^{-\beta S})\), meaning functional diversity plateaus at moderate species richness. This has profound implications for conservation: preserving functionally distinct species (e.g., nitrogen fixers, mycorrhizal hosts) matters more than preserving taxonomically rich but functionally redundant assemblages.
6.3 Metabolic Theory of Ecology: Predictions
The Metabolic Theory of Ecology (Brown et al. 2004) unifies diverse ecological patterns through the fundamental biochemical constraint that all organisms must metabolize. The theory begins with the allometric scaling of metabolic rate and the Boltzmann temperature dependence of biochemical reaction rates.
Starting Point: Kleiber’s Law with Temperature
Individual metabolic rate \(B\) scales with body mass \(M\) and temperature \(T\):
Fundamental MTE Equation
\[ B = B_0 \cdot M^{3/4} \cdot \exp\!\left(-\frac{E_a}{k_B T}\right) \]
where \(B_0\) is a normalization constant, \(E_a \approx 0.65\) eV is the average activation energy of metabolism (reflecting rate-limiting steps in the citric acid cycle and electron transport chain), \(k_B\) is Boltzmann’s constant, and \(T\) is absolute temperature.
Derivation: Population Density
The total energy flux through a population must equal the energy supply rate \(R\)(resource supply per unit area). If there are \(N\) individuals per unit area, each consuming at rate \(B\):
\[ N \cdot B = R \quad \Rightarrow \quad N = \frac{R}{B_0 M^{3/4} \exp(-E_a/k_B T)} \]
At constant temperature and resource supply:
Population Density Scaling
\[ N \propto M^{-3/4} \]
This explains why elephants are rare and bacteria are abundant: smaller organisms have lower per-capita metabolic demands, so more individuals can be supported per unit resource.
Derivation: Generation Time
The energy required to produce an offspring is proportional to adult body mass: \(E_{\text{repro}} \propto M\). The rate of energy allocation to reproduction is proportional to metabolic rate \(B \propto M^{3/4}\). Therefore, generation time \(\tau\) is:
\[ \tau = \frac{E_{\text{repro}}}{B} \propto \frac{M}{M^{3/4}} = M^{1/4} \]
Bacteria divide in hours; elephants reproduce every ~5 years; trees live for centuries. All are described by \(\tau \propto M^{1/4}\).
Derivation: Species Richness
The MTE predicts species richness by combining the temperature dependence of metabolic rate with the species-area relationship. More energy (\(\propto \exp(-E_a/k_B T)\)) supports more individuals, which supports more species via the species-individuals relationship \(S \propto N^z\):
MTE Species Richness Prediction
\[ \ln S = \ln S_0 + \frac{E_a}{k_B}\left(\frac{1}{T_0} - \frac{1}{T}\right) + z \ln A \]
This predicts that \(\ln S\) should increase linearly with \(1/k_B T\), explaining the latitudinal diversity gradient (more species in warmer tropics) as a consequence of the Boltzmann factor governing biochemical reaction rates.
Allen et al. (2002) tested this prediction across trees, amphibians, and marine fish: plotting \(\ln S\) vs \(-1/k_B T\) yielded slopes of ~0.65 eV, matching the average activation energy of metabolism. This remarkable result suggests that the global distribution of biodiversity is fundamentally constrained by the temperature sensitivity of mitochondrial biochemistry.
6.4 Species-Energy Relationship
The species-energy hypothesis proposes that areas with greater energy availability (sunlight, net primary productivity) support more species. The mechanism is biochemical: more energy means more total biomass production, creating more biochemical niches through finer resource partitioning.
Wright’s Species-Energy Theory
Wright (1983) proposed that species richness depends not on area alone but on total available energy, the product of area and energy per unit area:
Species-Energy Relationship
\[ S = c \cdot E^z \quad \text{where} \quad E = A \cdot \text{NPP} \]
where \(S\) is species richness, \(E\) is total energy (area \(\times\) NPP),\(c\) is a taxon-specific constant, and \(z \approx 0.5\text{--}1.0\) depending on the organism group.
Taking logarithms:
\[ \ln S = \ln c + z \cdot \ln E = \ln c + z \cdot [\ln A + \ln(\text{NPP})] \]
This generalizes the classical species-area relationship (\(S = cA^z\)) by incorporating energy availability. The mechanism involves the more individuals hypothesis: more NPP supports more total individuals (\(J\)), and more individuals sample more species from the regional species pool:
\[ J \propto \text{NPP} \cdot A, \quad S \propto J^{z'} \quad \Rightarrow \quad S \propto (\text{NPP} \cdot A)^{z'} \]
Biochemical Niche Partitioning
Higher energy availability enables finer partitioning of the biochemical niche space. In high-NPP tropical forests, light gradients from canopy to forest floor create distinct biochemical environments:
- Canopy: High light, high photorespiration rates, thick leaves with high \(V_{\max}\) for RuBisCO
- Mid-story: Moderate light, enhanced chlorophyll b/a ratios for efficient light capture
- Understory: Deep shade specialists with high antenna complex size, low light compensation point
- Forest floor: Saprophytic and mycoheterotrophic species, no photosynthesis at all
Each stratum represents a distinct biochemical niche. In low-energy environments (tundra, desert), these strata compress or disappear, supporting fewer species. The species-energy relationship thus reflects the thermodynamic constraint that biochemical niche space expands with available energy.
6.5 Biodiversity-Ecosystem Function (BEF)
Tilman’s landmark Cedar Creek experiments (beginning in 1994) provided the definitive experimental evidence that biodiversity enhances ecosystem function. Plots planted with 1, 2, 4, 8, or 16 grassland species showed that species-rich plots had higher aboveground biomass, better nutrient retention, and lower year-to-year variability.
The BEF Saturating Curve
The relationship between species richness and ecosystem function is typically a saturating curve, well described by Michaelis-Menten kinetics:
BEF Saturating Relationship
\[ F(S) = F_{\max} \cdot \frac{S}{S + K_S} \]
where \(F_{\max}\) is maximum function, \(S\) is species richness, and \(K_S\) is the half-saturation constant (species number at which function reaches half its maximum).
Additive Partitioning: Complementarity vs. Selection Effect
Loreau & Hector (2001) developed a method to partition the net biodiversity effect into two components: the complementarity effect(niche partitioning + facilitation) and the selection effect(dominance by productive species). The derivation proceeds as follows.
Let \(M_i\) be the monoculture yield of species \(i\) and\(Y_i\) its yield in mixture. The expected yield if species do not interact is\(E_i = p_i \cdot M_i\), where \(p_i = 1/S\) is the planted proportion. Define the relative yield deviation:
\[ \Delta \text{RY}_i = \frac{Y_i}{M_i} - p_i = \text{RY}_{Oi} - \text{RY}_{Ei} \]
The net biodiversity effect — the difference between observed mixture yield and expected yield based on monocultures — is:
\[ \Delta Y = Y_{\text{obs}} - Y_{\text{exp}} = \sum_{i=1}^{S} Y_i - \sum_{i=1}^{S} p_i M_i = \sum_{i=1}^{S} \Delta \text{RY}_i \cdot M_i \]
Using the identity \(\sum a_i b_i = N\overline{a}\cdot\overline{b} + N\text{cov}(a, b)\), where \(N = S\):
Loreau-Hector Partition (2001)
\[ \Delta Y = \underbrace{S \cdot \overline{\Delta \text{RY}} \cdot \overline{M}}_{\text{Complementarity effect}} + \underbrace{S \cdot \text{cov}(\Delta \text{RY}, M)}_{\text{Selection effect}} \]
- Complementarity effect (\(S \cdot \overline{\Delta \text{RY}} \cdot \overline{M}\)): Positive when species, on average, yield more in mixture than expected from monocultures. This occurs through biochemical niche partitioning (e.g., deep-rooted vs. shallow-rooted species accessing different soil nutrient pools) or facilitation (e.g., legume N₂fixation benefiting neighboring grasses).
- Selection effect (\(S \cdot \text{cov}(\Delta \text{RY}, M)\)): Positive when species with high monoculture yields also gain relatively more in mixtures. This reflects competitive dominance by the most productive species, essentially a “sampling effect” — more species means a higher chance of including the most productive one.
Meta-analyses show that in most BEF experiments, the complementarity effect dominates after multiple years, while the selection effect is more important initially. This suggests that biochemical niche partitioning, not mere dominance, is the primary long-term mechanism linking biodiversity to ecosystem function.
6.6 BEF Relationship & Portfolio Effect
The diagram below illustrates two key concepts: (left) the saturating relationship between species richness and ecosystem productivity observed in BEF experiments, and (right) the portfolio effect showing how variance in ecosystem function decreases with species diversity.
Left panel: Cedar Creek-type BEF experiment data showing the saturating relationship between plant species richness and aboveground productivity. The curve follows \(F = F_{\max} \cdot S/(S + K_S)\). Most of the productivity gain occurs in the first 4–8 species due to complementarity in resource use.
Right panel: Portfolio effect — coefficient of variation (CV) of ecosystem function decreases with species richness. The rate of decrease depends on species response correlation (\(\rho\)). When species respond independently (\(\rho = 0\)), CV decreases as \(1/\sqrt{S}\). Higher correlation reduces the insurance benefit.
6.7 Computational Simulations
We explore the biodiversity-stability relationship through three simulations: (1) the portfolio effect demonstrating variance reduction with diversity, (2) MTE predictions of species richness across latitudes, and (3) the Loreau-Hector partition of complementarity vs. selection effects.
Portfolio Effect: Variance Reduction with Biodiversity
PythonClick Run to execute the Python code
Code will be executed with Python 3 on the server
References
- Yachi, S. & Loreau, M. (1999). Biodiversity and ecosystem productivity in a fluctuating environment: the insurance hypothesis. Proceedings of the National Academy of Sciences, 96(4), 1463–1468.
- Rao, C.R. (1982). Diversity and dissimilarity coefficients: a unified approach. Theoretical Population Biology, 21(1), 24–43.
- Villeger, S., Mason, N.W.H. & Mouillot, D. (2008). New multidimensional functional diversity indices for a multifaceted framework in functional ecology. Ecology, 89(8), 2290–2301.
- Brown, J.H., Gillooly, J.F., Allen, A.P., Savage, V.M. & West, G.B. (2004). Toward a metabolic theory of ecology. Ecology, 85(7), 1771–1789.
- Allen, A.P., Brown, J.H. & Gillooly, J.F. (2002). Global biodiversity, biochemical kinetics, and the energetic-equivalence rule. Science, 297(5586), 1545–1548.
- Wright, D.H. (1983). Species-energy theory: an extension of species-area theory. Oikos, 41(3), 496–506.
- Tilman, D., Wedin, D. & Knops, J. (1996). Productivity and sustainability influenced by biodiversity in grassland ecosystems. Nature, 379(6567), 718–720.
- Loreau, M. & Hector, A. (2001). Partitioning selection and complementarity in biodiversity experiments. Nature, 412(6842), 72–76.
- Tilman, D., Reich, P.B. & Knops, J.M.H. (2006). Biodiversity and ecosystem stability in a decade-long grassland experiment. Nature, 441(7093), 629–632.
- Doak, D.F., Bigger, D., Harding, E.K., Marvier, M.A., O’Malley, R.E. & Thomson, D. (1998). The statistical inevitability of stability-diversity relationships in community ecology. The American Naturalist, 151(3), 264–276.