Insect Decline & Conservation Biophysics
Quantifying the crisis: biomass decline, neonicotinoid toxicology, light pollution, and habitat fragmentation
8.1 The Krefeld Study: Quantifying Insect Decline
In 2017, Hallmann et al. published what became the most widely cited study in insect conservation: a 27-year dataset from the Krefeld Entomological Society in Germany, showing a 75% decline in flying insect biomass between 1989 and 2016. The data came from standardized malaise traps deployed across 63 protected nature reserves.
Exponential Decline Model
The decline follows an approximately exponential trajectory:
\(B(t) = B_0 \cdot e^{-\lambda t}\)
To find \(\lambda\): after 27 years, \(B(27) = 0.25 B_0\), so:
\(0.25 = e^{-\lambda \cdot 27} \quad \Rightarrow \quad \lambda = \frac{\ln 4}{27} \approx 0.051 \text{ yr}^{-1}\)
The halving time is:
\(t_{1/2} = \frac{\ln 2}{\lambda} = \frac{0.693}{0.051} \approx 13.5 \text{ years}\)
Every ~14 years, half of the remaining insect biomass is lost
Why This Matters
The Krefeld data was alarming because (1) the traps were in protected areas, not farmland; (2) the decline occurred regardless of habitat type, weather, or land-use changes in the surrounding landscape; and (3) subsequent studies in Puerto Rico (Lister & Garcia 2018), Denmark, and the Netherlands confirmed similar trends. This is not a local phenomenon — it appears to be global.
8.2 Neonicotinoid Biochemistry
Neonicotinoids (imidacloprid, clothianidin, thiamethoxam) are the world's most widely used insecticides, accounting for ~25% of the global insecticide market. They are systemic: applied as seed coatings, they are taken up by the plant and distributed to all tissues, including pollen and nectar.
Mechanism: nAChR Agonism
Neonicotinoids act as agonists at nicotinic acetylcholine receptors (nAChRs)in the insect central nervous system. They bind the same site as acetylcholine but are not degraded by acetylcholinesterase, causing sustained neuronal excitation followed by receptor desensitization, paralysis, and death.
Insect vs Mammalian Selectivity
The ~1000x selectivity for insect over mammalian nAChRs arises from a single amino acid difference at the binding loop:
\(\text{Insect nAChR: } \textbf{Tyr} \text{ at loop D} \quad \Rightarrow \quad K_d \approx 1 \text{ nM (high affinity)}\)
\(\text{Mammal nAChR: } \textbf{Arg} \text{ at loop D} \quad \Rightarrow \quad K_d \approx 1{,}000 \text{ nM (low affinity)}\)
The cationic Arg residue in mammalian receptors electrostatically repels the partially positive neonicotinoid nitrogen, while the neutral Tyr in insect receptors allows tight binding. This is why neonicotinoids were initially considered “safe” — but the sublethal effects on bee behavior were overlooked.
Dose-Response: LD\(_{50}\) Comparison
The dose-response relationship follows a sigmoidal (Hill equation) curve:
\(\text{Mortality}(D) = \frac{D^n}{D^n + \text{LD}_{50}^n}\)
where \(n\) = Hill coefficient (~2.5), \(D\) = dose (ng/bee)
| Compound | LD\(_{50}\) (ng/bee) | Rat oral LD\(_{50}\) (mg/kg) | Selectivity ratio |
|---|---|---|---|
| Imidacloprid | 3.7 | 450 | ~1,200x |
| Clothianidin | 2.5 | >5,000 | >2,000x |
| Thiamethoxam | 5.0 | 1,563 | ~300x |
8.3 Light Pollution: ALAN and Insect Navigation
Artificial Light at Night (ALAN) has increased globally by ~2% per year over the past century. For nocturnal insects, which constitute the majority of insect species, this represents a catastrophic disruption to navigation, reproduction, and predator-prey dynamics.
Moth Attraction: The Logarithmic Spiral
Moths navigate by maintaining a constant angle \(\alpha\) to distant light sources (the Moon, stars). When the light source is nearby (artificial), maintaining this constant angle produces a logarithmic spiral that draws the moth inexorably toward the light:
\(r(\theta) = r_0 \cdot e^{\theta \cdot \cot\alpha}\)
where \(r_0\) = initial distance, \(\alpha\) = angle between flight path and radial direction, \(\theta\) = cumulative angle traversed
For \(\alpha < 90^\circ\), \(\cot\alpha > 0\) and the moth spirals outward (normal celestial navigation). For \(\alpha > 90^\circ\) (attempted celestial navigation toward a nearby light),\(\cot\alpha < 0\) and the moth spirals inward toward the light, becoming trapped.
Broader Impacts of ALAN
- Circadian disruption: melatonin suppression, altered oviposition timing
- Predation amplification: congregations around lights create “ecological traps” for bats and spiders
- Pollination disruption: nocturnal pollinators (moths) fail to visit flowers, reducing seed set by up to 62% (Knop et al. 2017)
- Population sinks: street lights attract and kill billions of insects nightly, removing them from the breeding population
8.4 Habitat Fragmentation: Island Biogeography Applied
MacArthur and Wilson's Theory of Island Biogeography (1967), originally developed for oceanic islands, provides the mathematical framework for understanding how habitat fragmentation affects insect diversity.
Species-Area Relationship
\(S = cA^z\)
where \(S\) = number of species, \(A\) = area, \(c\) = taxon-specific constant, \(z \approx 0.15\text{-}0.35\)
The exponent \(z\) determines sensitivity to area loss. For mainland habitat fragments,\(z \approx 0.15\text{-}0.25\); for isolated fragments, \(z \approx 0.25\text{-}0.35\)(approaching island values). The species loss from halving the area is:
\(\frac{S_{A/2}}{S_A} = \left(\frac{A/2}{A}\right)^z = 2^{-z}\)
For \(z = 0.25\): losing 50% of habitat area eliminates ~16% of species
Extinction Debt
Extinction debt is the delayed loss of species following habitat fragmentation. Species may persist for decades after their habitat is reduced below the minimum viable area, but are committed to extinction — they are “living dead” species. The relaxation time to equilibrium follows:
\(S(t) = S_{\text{eq}} + (S_0 - S_{\text{eq}}) \cdot e^{-t/\tau}\)
where \(\tau\) = relaxation time (decades to centuries for insects with small populations)
Minimum Viable Population
From demographic stochasticity, the probability of extinction for a small population of size \(N\) is:
\(P_{\text{ext}}(t) \approx 1 - e^{-t/\tau_N}, \quad \text{where } \tau_N \propto e^{cN}\)
Mean time to extinction grows exponentially with population size — small populations are doomed by chance alone
8.5 The Pollinator Crisis: Allee Effects and Economic Consequences
Approximately 75% of global food crop species depend on animal pollinators, with an estimated economic value of $235-577 billion USD per year (IPBES 2016). The relationship between pollinator density and crop yield is not linear — it follows anAllee effect (Hill function):
\(f(P) = f_{\max} \cdot \frac{P^n}{P^n + K^n}\)
where \(P\) = pollinator density (visits/plant/day), \(K\) = half-saturation constant, \(n\) = Hill coefficient (cooperativity)
The Hill coefficient \(n > 1\) creates a sharp threshold: fruit set is near zero below a critical pollinator density and saturates above it. This means moderate pollinator declines may have disproportionately large effects on crop yields — the classic Allee effect.
Crop-Specific Vulnerability
Different crops have different \(K\) values reflecting their pollinator dependency:
- Almonds: \(K \approx 12\) visits/day — 100% pollinator-dependent, no self-pollination
- Apples: \(K \approx 8\) — high dependency, some self-incompatible cultivars
- Blueberries: \(K \approx 6\) — buzz-pollination specialists required
- Tomatoes: \(K \approx 3\) — can self-pollinate but yield improves with buzz pollinators
Simulation: Insect Decline Models
This simulation models (1) insect biomass decline with future projections, (2) neonicotinoid dose-response curves for three compounds, (3) moth flight paths as logarithmic spirals toward artificial lights, and (4) species-area relationships for habitat fragments with different isolation levels.
Click Run to execute the Python code
Code will be executed with Python 3 on the server
Advanced: Pollination Deficit & Multi-Stressor Synergy
This simulation visualizes (1) the pollination Allee effect across crops with different pollinator dependencies, and (2) how multiple simultaneous stressors interact synergistically to produce population declines greater than the sum of individual effects.
Click Run to execute the Python code
Code will be executed with Python 3 on the server
References
- Hallmann, C.A. et al. (2017). More than 75 percent decline over 27 years in total flying insect biomass in protected areas. PLoS ONE, 12(10), e0185809.
- Lister, B.C. & Garcia, A. (2018). Climate-driven declines in arthropod abundance restructure a rainforest food web. Proceedings of the National Academy of Sciences, 115(44), E10397-E10406.
- Sánchez-Bayo, F. & Wyckhuys, K.A.G. (2019). Worldwide decline of the entomofauna: a review of its drivers. Biological Conservation, 232, 8-27.
- Jeschke, P. et al. (2011). Overview of the status and global strategy for neonicotinoids. Journal of Agricultural and Food Chemistry, 59(7), 2897-2908.
- Tomizawa, M. & Casida, J.E. (2005). Neonicotinoid insecticide toxicology: mechanisms of selective action. Annual Review of Pharmacology and Toxicology, 45, 247-268.
- Woodcock, B.A. et al. (2017). Country-specific effects of neonicotinoid pesticides on honey bees and wild bees. Science, 356(6345), 1393-1395.
- Owens, A.C.S. et al. (2020). Light pollution is a driver of insect declines. Biological Conservation, 241, 108259.
- Knop, E. et al. (2017). Artificial light at night as a new threat to pollination. Nature, 548(7666), 206-209.
- MacArthur, R.H. & Wilson, E.O. (1967). The Theory of Island Biogeography. Princeton University Press.
- Hanski, I. (2011). Habitat loss, the dynamics of biodiversity, and a perspective on conservation. Ambio, 40(3), 248-255.
- IPBES (2016). Assessment Report on Pollinators, Pollination and Food Production. Secretariat of IPBES, Bonn.
- Klein, A.-M. et al. (2007). Importance of pollinators in changing landscapes for world crops. Proceedings of the Royal Society B, 274(1608), 303-313.