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
6.8 Keystone Species & Biochemical Cascades
Robert Paine (1966) introduced the keystone species concept after demonstrating that removing the predatory starfish Pisaster ochraceus from intertidal communities caused a dramatic collapse in species diversity. The key insight is that certain species exert disproportionate control over ecosystem biochemistry relative to their biomass. Removing a keystone species triggers biochemical cascades that propagate through multiple trophic levels.
The Sea Otter Cascade: A Biochemical Chain Reaction
The North Pacific kelp forest system provides a paradigmatic example of how keystone species removal cascades through ecosystem biochemistry:
- Step 1: Sea otter removal (fur trade hunting) eliminates top predation on sea urchins
- Step 2: Sea urchin populations explode (10–100\(\times\) increase), overgrazing kelp holdfasts
- Step 3: Kelp forest destruction eliminates \(\sim\)1,500 g C m\(^{-2}\) yr\(^{-1}\) of NPP
- Step 4: Loss of dissolved organic carbon (DOC) from kelp exudates reduces microbial production by 60–80%
- Step 5: Reduced microbial loop weakens the biological pump, decreasing CO₂ sequestration
- Step 6: Net CO₂ release to atmosphere — a single trophic perturbation becomes a biogeochemical perturbation
Quantifying Trophic Cascade Strength
The strength of a trophic cascade depends on the product of interaction strengths between adjacent trophic levels. Consider a linear food chain with \(n\) trophic levels, where \(\alpha_{ij}\) is the per-capita interaction strength (effect of species\(j\) on species \(i\)). When the top predator is removed, causing a change \(\Delta B_{\text{top}}\) in top predator biomass, the resulting change at trophic level \(i\) is:
Trophic Cascade Propagation
\[ |\Delta B_i| = \prod_{j=i}^{n-1} |\alpha_{j,j+1}| \times |\Delta B_{\text{top}}| \]
The cascade strength is the product of all pairwise interaction strengths along the chain. If any single interaction is weak (\(|\alpha_{ij}| \ll 1\)), the cascade attenuates. Keystone species maintain strong interactions (\(|\alpha| > 1\)) that amplify perturbations rather than damping them.
Paine’s keystone concept can be quantified using the community importance index (CI):
Community Importance Index
\[ \text{CI}_k = \frac{S_{\text{total}} - S_{-k}}{S_{\text{total}}} \cdot \frac{1}{p_k} \]
where \(S_{\text{total}}\) is species richness with species \(k\) present,\(S_{-k}\) is richness after removing species \(k\), and \(p_k\) is the proportional abundance of species \(k\). A keystone species has CI \(\gg 1\), meaning its effect on diversity is disproportionate to its abundance.
In alternating trophic cascades, even-numbered levels from the top experience release(population increase) while odd-numbered levels experience suppression. The sign alternates: \(\text{sign}(\Delta B_i) = (-1)^{n-i} \cdot \text{sign}(\Delta B_{\text{top}})\). This explains why otter removal (decrease at level 4) causes urchin increase (level 3) and kelp decrease (level 2).
6.9 Stoichiometric Constraints on Food Webs
Ecological stoichiometry (Sterner & Elser 2002) examines how the balance of chemical elements — primarily carbon (C), nitrogen (N), and phosphorus (P) — shapes ecological interactions. The fundamental insight is that organisms maintain relatively fixed elemental ratios (stoichiometric homeostasis), but their food sources vary enormously in composition, creating biochemical constraints on growth and trophic transfer.
Growth Rate Limitation by Elemental Mismatch
When a consumer with fixed C:N:P requirements feeds on a resource with different elemental ratios, the most limiting element determines growth rate. Following Liebig’s law of the minimum extended to multiple elements:
Multi-Element Growth Limitation
\[ \mu = \mu_{\max} \times \min\!\left(1,\; \frac{(C\text{:}N)_{\text{food}}}{(C\text{:}N)_{\text{consumer}}},\; \frac{(C\text{:}P)_{\text{food}}}{(C\text{:}P)_{\text{consumer}}}\right) \]
When the food has a higher C:N ratio than the consumer requires (more C relative to N), growth is N-limited. The consumer must void excess C (as CO₂ via respiration or as DOC excretion) to maintain its internal stoichiometry.
The threshold elemental ratio (TER) defines the critical food composition where limitation switches from one element to another:
\[ \text{TER}_{C:N} = \frac{(C\text{:}N)_{\text{consumer}}}{\text{GGE}_C / \text{GGE}_N} \]
where GGE is the gross growth efficiency for each element. When the food C:N exceeds TER, the consumer is N-limited; below TER, it is C-limited.
Stoichiometric Mismatch Reduces Transfer Efficiency
The mismatch between producer and consumer stoichiometry has profound consequences for trophic transfer efficiency. Consider the Daphnia-algae system:
- Daphnia body composition: C:N:P \(\approx\) 100:18:1 (P-rich, reflecting high RNA content)
- Nutrient-replete algae: C:N:P \(\approx\) 106:16:1 (Redfield ratio) — good match, high growth efficiency
- P-depleted algae: C:N:P \(\approx\) 500:30:1 — severe P limitation for Daphnia, \(\sim\)80% of ingested C must be voided
The gross growth efficiency (GGE) for carbon drops from ~40% when food is stoichiometrically matched to \(<\)10% under severe mismatch. This wasted carbon is released as CO₂, effectively shunting energy from the food web into respiratory losses.
The Growth Rate Hypothesis
Elser et al. (2003) proposed the growth rate hypothesis: fast-growing organisms require more phosphorus because rapid growth demands high rates of protein synthesis, which requires abundant ribosomes, and ribosomes are \(\sim\)85% ribosomal RNA (rRNA). RNA is ~9% phosphorus by mass. The connection to an “RNA world”:
Growth Rate–RNA–Phosphorus Link
\[ \mu \propto [\text{RNA}] \propto [\text{ribosomes}] \propto [P_{\text{body}}] \]
Quantitatively: body %P \(\approx 0.7 + 0.3 \times \mu\) (where \(\mu\) is specific growth rate in d\(^{-1}\)). This predicts that body C:P ratio decreases as growth rate increases — confirmed across bacteria, algae, zooplankton, and insects.
The growth rate hypothesis has far-reaching implications: ecosystems with low P availability (oligotrophic lakes, ancient soils) should favor slow-growing, high C:P species, while P-rich ecosystems (eutrophic systems, young volcanic soils) can support fast-growing, P-demanding species. This creates a stoichiometric filter on community composition that operates independently of traditional competitive mechanisms.
6.10 Resilience, Resistance & Alternative Stable States
Resilience is the magnitude of perturbation a system can absorb before shifting to an alternative stable state. Resistance is the degree to which a system remains unchanged during a perturbation. These are distinct properties: a system can be highly resistant (difficult to perturb) but have low resilience (if perturbed beyond a threshold, it collapses irreversibly).
Holling’s Adaptive Cycle
C.S. Holling (1973) proposed that ecosystems cycle through four phases, each characterized by distinct biochemical dynamics:
- Exploitation (r): Rapid colonization and growth. Pioneer species with fast biochemistry (high metabolic rates, simple secondary chemistry) capture resources. Nutrients flow from abiotic pools into biomass.
- Conservation (K): Slow accumulation and storage. Mature species with complex biochemistry (lignin, tannins, allelopathic compounds) dominate. Nutrients are locked in biomass and recalcitrant organic matter. System is rigid and brittle.
- Release (\(\Omega\)): Rapid collapse triggered by disturbance (fire, disease, drought). Stored biomass is rapidly mineralized, releasing nutrients as pulses of NH₄\(^+\), PO₄\(^{3-}\), and DOC.
- Reorganization (\(\alpha\)): Novel recombination of resources. Microbial communities process the nutrient pulse; mycorrhizal networks may persist as “biological memory.” The system can reorganize into the same or an alternative state.
Ball-in-Cup: Alternative Stable States
The ball-in-cup analogy visualizes resilience as the depth of a basin of attraction. The system state (ball) rests in a potential well (cup). Small perturbations are absorbed as the ball rolls back. If the perturbation exceeds the basin depth, the ball rolls into an alternative basin — an irreversible regime shift.
Catastrophe Theory: Fold Bifurcation
The mathematical framework for alternative stable states comes from catastrophe theory. A simple model exhibiting bistability is the fold (cusp) bifurcation:
Fold Bifurcation Model
\[ \frac{dx}{dt} = r + hx - x^3 \]
where \(x\) is the system state (e.g., water clarity), \(r\) is a slow driver (e.g., nutrient loading rate), and \(h\) is a symmetry-breaking parameter. The cubic term creates the possibility of multiple stable equilibria.
Setting \(dx/dt = 0\) and analyzing stability: the system has either one or three equilibria depending on the parameter values. When \(r\) is slowly varied (increasing nutrient load), the system tracks one equilibrium until a critical threshold \(r_{\text{crit}}\)where that equilibrium disappears and the system “catastrophically” jumps to the alternative state:
Critical Transition Point
\[ r_{\text{crit}} = \pm \frac{2}{3\sqrt{3}} h^{3/2} \]
The system exhibits hysteresis: the forward transition (clear \(\to\) turbid) occurs at a different \(r\) value than the reverse (turbid \(\to\) clear). This means that simply reducing nutrient loading to pre-eutrophication levels is insufficient to restore the clear-water state.
Lake Eutrophication: A Classic Bistable System
Shallow lakes exemplify alternative stable states driven by biochemical feedbacks:
- Clear-water state: Submerged macrophytes stabilize sediments (reducing resuspension), release allelopathic compounds that inhibit phytoplankton, provide refuge for Daphnia grazers, and oxygenate sediments (preventing P release). Positive feedback maintains clarity.
- Turbid state: Phytoplankton dominance shades macrophytes (light limitation), wind-driven resuspension increases turbidity further, anoxic sediments release P (Einsele-Mortimer mechanism: Fe(III) \(\to\) Fe(II) releases bound P), and Daphnia populations crash (no macrophyte refuge from fish predation). Positive feedback maintains turbidity.
The critical phosphorus loading that triggers the transition from clear to turbid is typically 25–100 \(\mu\)g P/L (total P), but the reverse transition requires reducing TP below 25–50 \(\mu\)g P/L — the hysteresis gap makes restoration much harder than prevention. This has been confirmed empirically in hundreds of lakes worldwide (Scheffer et al. 2001).
6.11 Advanced Simulations: Cascades, Stoichiometry & Resilience
Three simulations exploring the new topics: (1) a 4-level trophic cascade model showing biomass responses to keystone predator removal, (2) stoichiometric growth limitation across C:N:P space, and (3) a resilience landscape with bifurcation diagram for the lake eutrophication system.
Trophic Cascades, Stoichiometric Limitation & Resilience Landscapes
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.
- Paine, R.T. (1966). Food web complexity and species diversity. The American Naturalist, 100(910), 65–75.
- Estes, J.A., Terborgh, J., Brashares, J.S., Power, M.E., Berger, J., Bond, W.J., ... & Wardle, D.A. (2011). Trophic downgrading of planet Earth. Science, 333(6040), 301–306.
- Sterner, R.W. & Elser, J.J. (2002). Ecological Stoichiometry: The Biology of Elements from Molecules to the Biosphere. Princeton University Press.
- Elser, J.J., Acharya, K., Kyle, M., Cotner, J., Makino, W., Markow, T., ... & Sterner, R.W. (2003). Growth rate–stoichiometry couplings in diverse biota. Ecology Letters, 6(10), 936–943.
- Holling, C.S. (1973). Resilience and stability of ecological systems. Annual Review of Ecology and Systematics, 4(1), 1–23.
- Scheffer, M., Carpenter, S., Foley, J.A., Folke, C. & Walker, B. (2001). Catastrophic shifts in ecosystems. Nature, 413(6856), 591–596.