Module 5: Migration & Phenology
Climate change is reshuffling the geography and timing of life. Species are shifting poleward and upslope, phenological events are advancing asymmetrically, and novel species combinations are emerging. This module derives range shift velocities, the fitness cost of phenological mismatch, invasion front dynamics, and community dissimilarity metrics that quantify the emergence of no-analog ecosystems.
1. Range Shifts & Climate Velocity
A landmark meta-analysis by Chen et al. (2011) documented that species are shifting their ranges poleward at a median rate of 16.9 km per decade and upslope at 11 m per decade. However, these shifts vary enormously across taxa, with marine fish moving fastest (~72 km/decade) and amphibians slowest (~3.5 km/decade).
Climate Velocity
The concept of climate velocity quantifies how fast organisms must move to stay within their thermal envelope. It is the ratio of the temporal temperature gradient to the spatial temperature gradient:
$$v_{\text{climate}} = \frac{\partial T/\partial t}{|\nabla_{\text{spatial}} T|}$$
Current rate: $\partial T/\partial t \approx 0.02$°C/yr; spatial gradient: $|\nabla T| \approx 0.005$°C/km in midlatitudes
This gives a current latitudinal climate velocity of approximately 4 km/yr (40 km/decade). In mountainous terrain, the altitudinal lapse rate ($\Gamma \approx 5.5$°C/km) means organisms need only move ~4 m/yr upslope to track the same temperature change. Mountains thus provide short-distance climate refugia.
Dispersal Limitation & Climatic Debt
When climate velocity exceeds a species’ dispersal capacity, a climatic debt accumulates: the species occupies conditions warmer than its optimum and has not yet colonized suitable habitat at higher latitudes. The tracking ratio quantifies this:
$$\text{Tracking ratio} = \frac{v_{\text{species}}}{v_{\text{climate}}}$$
Ratio < 1 indicates climatic debt; trees (ratio ~ 0.3) accumulate the largest debt
The Mountain-Top Extinction Trap
Species adapted to high-altitude or high-latitude conditions face a geometric constraint: as they shift upslope, the available area decreases (mountains narrow toward the summit). This “escalator to extinction” (Freeman et al., 2018) means that even species that can track climate velocity may face extinction through habitat compression:
$$A(z) \propto e^{-z/z_0} \quad \Longrightarrow \quad \frac{dA}{dz} < 0$$
Area decreases exponentially with altitude, so upslope shifts reduce habitat even when tracking succeeds
Range Shifts: Poleward & Upslope Movement
2. Phenological Mismatch
Phenology—the timing of seasonal biological events—is advancing across the globe. Spring events (budburst, emergence, breeding) are occurring approximately 2.3 days earlier per decade (Parmesan & Yohe, 2003). However, different trophic levels advance at different rates, creating mismatches between consumers and their resources.
The Fitness Cost of Mismatch
When a consumer’s peak demand (e.g., chick-rearing) no longer coincides with peak resource availability (e.g., caterpillar biomass), fitness declines. The relationship follows a Gaussian function:
$$W = W_{\max} \exp\!\left(-\frac{(\Delta t)^2}{2\sigma^2}\right)$$
$\Delta t$ = temporal mismatch (days); $\sigma$ = tolerance window (5–8 days for insectivores)
The Great Tit–Caterpillar System
The best-documented case of phenological mismatch involves the great tit (Parus major) and its caterpillar prey in European oak forests (Both et al., 2006; Visser et al., 2006). The caterpillar peak has advanced ~3 days/decade (directly tracking spring temperatures), while great tit egg-laying has advanced only ~1.5 days/decade (partially constrained by photoperiod cues). The growing mismatch reduces nestling survival and fledgling weight.
Long-Distance Migrants: The Pied Flycatcher
Long-distance migrants face a more severe problem. The pied flycatcher (Ficedula hypoleuca) overwinters in West Africa and times its spring migration using photoperiod—a cue that does not change with climate. Both et al. (2006) documented population declines of up to 90% in areas where the caterpillar peak had advanced most, because the flycatchers arrived too late to exploit the food pulse:
$$\frac{dN}{dt} = r_{\max}\left(2W(\Delta t) - 1\right) \cdot N$$
Population growth becomes negative when $W < 0.5$, i.e., when mismatch exceeds $\sigma\sqrt{2\ln 2} \approx 1.18\sigma$
Phenological Mismatch: Trophic Timing Gaps
3. Migratory Disruption
Arctic-breeding shorebirds exemplify the challenge of migratory disruption. Species like the red knot (Calidris canutus) time their arrival at Arctic breeding grounds to coincide with the peak emergence of invertebrate prey. However, their departure from tropical wintering grounds is cued by photoperiod, while Arctic food peaks are shifting with temperature.
Optimal Migration Timing
The optimal departure date can be derived from an energy budget model. The net energy gain at the breeding ground is:
$$\frac{dE}{dt} = I(t) - C(t)$$
$I(t)$ = food intake rate (peaks with invertebrate emergence), $C(t)$ = metabolic cost (higher in cold conditions)
The bird must accumulate sufficient energy $E_{\text{breed}}$ for breeding within the narrow Arctic season. Optimal arrival maximizes the integral:
$$E_{\text{net}} = \int_{t_{\text{arrive}}}^{t_{\text{depart}}} [I(t) - C(t)]\, dt$$
As the food peak shifts earlier while arrival remains fixed, the overlap between $I(t)$ and the bird’s presence window decreases, reducing $E_{\text{net}}$ and breeding success. van Gils et al. (2016) showed that red knots arriving late relative to the food peak produced smaller offspring with shorter bills, reducing their survival by 30%.
4. Invasive Species Facilitation
Climate change facilitates biological invasions by removing the thermal barriers that previously limited species’ ranges. As temperatures warm, invasive species from lower latitudes can expand poleward into territory previously too cold for establishment.
Invasion Front Velocity (Fisher-KPP)
The speed of an invasion front is predicted by the Fisher-Kolmogorov-Petrovsky-Piskunov (KPP) equation, which combines population growth with spatial diffusion:
$$\frac{\partial n}{\partial t} = rn\left(1 - \frac{n}{K}\right) + D\frac{\partial^2 n}{\partial x^2}$$
The minimum invasion front speed is:
$$v = 2\sqrt{rD}$$
$r$ = intrinsic growth rate (temperature-dependent), $D$ = diffusion coefficient (dispersal capacity)
As warming increases local temperature toward the invader’s optimum, $r$ increases and the invasion front accelerates.
Case Studies
Lionfish (Pterois volitans) invaded the western Atlantic from the Indo-Pacific, spreading along the US East Coast at rates exceeding 10 km/yr. Warming waters have extended their viable range northward past Cape Hatteras (Whitfield et al., 2014). They reduce native reef fish recruitment by up to 79%.
Asian tiger mosquito (Aedes albopictus) has expanded across southern Europe, facilitated by warming winters that no longer kill overwintering eggs. By 2050, climate models project suitable habitat extending to Scandinavia (Kraemer et al., 2015), bringing dengue and chikungunya risk to 500+ million additional people.
5. Community Reshuffling & No-Analog Ecosystems
Because species respond individualistically to climate change—each tracking its own thermal niche at its own dispersal rate—existing species assemblages are being disassembled and novel combinations are forming. Williams et al. (2007) termed these no-analog communities: species combinations with no historical precedent.
Bray-Curtis Dissimilarity
Community change can be quantified using the Bray-Curtis dissimilarity index between a community at time $t$ and its historical baseline:
$$BC_{ij} = 1 - \frac{2\sum_k \min(n_{ik}, n_{jk})}{\sum_k n_{ik} + \sum_k n_{jk}}$$
$BC = 0$: identical communities; $BC = 1$: completely different species; $BC > 0.5$: no-analog threshold
The Bray-Curtis index increases approximately as $BC \approx 1 - e^{-\beta \Delta T}$, where $\beta$ depends on the sensitivity of the biome. Tropical communities are most sensitive ($\beta \approx 0.3$) because species there are already near their upper thermal limits and small temperature changes cause disproportionate turnover.
At +3°C warming, Williams et al. estimated that 12–39% of Earth’s land surface would host no-analog communities. These novel ecosystems present challenges for conservation because historical baselines no longer apply—we cannot “restore” an ecosystem to a state that never existed.
Simulation: Range Shift Velocities
Poleward and altitudinal range shifts across taxa, climate velocity vs dispersal capacity, elevational range squeeze, and climatic debt tracking ratios.
Click Run to execute the Python code
Code will be executed with Python 3 on the server
Simulation: Phenological Mismatch & Population Dynamics
Phenological shift timelines, Gaussian fitness cost model, population trajectories for partially-tracking (great tit) vs poorly-tracking (pied flycatcher) species, and sensitivity to advancement rates.
Click Run to execute the Python code
Code will be executed with Python 3 on the server
Simulation: Invasion Fronts & Community Dissimilarity
Fisher-KPP invasion velocities for key invasive species, cumulative front spread projections, Bray-Curtis community dissimilarity under warming, and no-analog community area projections.
Click Run to execute the Python code
Code will be executed with Python 3 on the server
Key References
• Chen, I.-C. et al. (2011). “Rapid range shifts of species associated with high levels of climate warming.” Science, 333, 1024–1026.
• Parmesan, C. & Yohe, G. (2003). “A globally coherent fingerprint of climate change impacts across natural systems.” Nature, 421, 37–42.
• Both, C. et al. (2006). “Climate change and population declines in a long-distance migratory bird.” Nature, 441, 81–83.
• Visser, M. E. et al. (2006). “Shifts in caterpillar biomass phenology due to climate change and its impact on the breeding biology of an insectivorous bird.” Oecologia, 147, 164–172.
• Freeman, B. G. et al. (2018). “Climate change causes upslope shifts and mountaintop extirpations in a tropical bird community.” Proceedings of the National Academy of Sciences, 115, 11982–11987.
• van Gils, J. A. et al. (2016). “Body shrinkage due to Arctic warming reduces red knot fitness in tropical wintering range.” Science, 352, 819–821.
• Williams, J. W. et al. (2007). “Projected distributions of novel and disappearing climates by 2100 AD.” Proceedings of the National Academy of Sciences, 104, 5738–5742.
• Whitfield, P. E. et al. (2014). “Native fish community structure and Indo-Pacific lionfish Pterois volitans densities along a depth-temperature gradient in Onslow Bay, North Carolina, USA.” Marine Ecology Progress Series, 509, 241–254.
• Kraemer, M. U. G. et al. (2015). “The global distribution of the arbovirus vectors Aedes aegypti and Ae. albopictus.” eLife, 4, e08347.
• Loarie, S. R. et al. (2009). “The velocity of climate change.” Nature, 462, 1052–1055.