Module 10: Temperate & Boreal Biomes
The temperate and boreal biomes span Europe, North America, and Siberia and together cover over a third of the terrestrial surface. They are undergoing climate-driven range reorganisation: treelines advance, southern range edges contract, boreal fire regimes intensify, and tightly co-evolved trophic chains (oak → winter moth → pied flycatcher) slip out of phase. This module derives the ecological theory behind latitudinal and altitudinal range shifts, phenological mismatch, and fire-climate feedbacks, and closes with two quantitative projections: an oak/beech species-distribution model and a simulated pied-flycatcher mismatch trajectory.
1. Latitudinal & Altitudinal Range Shifts
Two meta-analyses anchor our quantitative view of recent range shifts. Parmesan & Yohe (2003) and Parmesan (2006) pooled 1700+ species and reported an average poleward shift of 6.1 km per decade and an upslope shift of 6.1 m/decade. Chen et al. (2011, Science) revised these to 16.9 km/decade poleward and 11.0 m/decade upslopewith faster shifts in regions of stronger warming.
The theoretical expectation for the latitudinal shift \(\Delta y\) needed to conserve climatic niche is:
\[\Delta y = \frac{\Delta T}{\Gamma_{\text{lat}}}, \quad \Gamma_{\text{lat}} \approx 0.6\;\text{°C / 100 km in mid-latitudes}\]
A regional warming of \(+2\)°C demands a \(\sim 330\) km poleward shift. Few temperate trees achieve more than 0.5 km/yr, highlighting a “climatic debt”.
Treeline Advance
Harsch et al. (2009, Ecology Letters) meta-analysed 166 sites globally and found 52% of treelines advancing, 47% static, 1% retreating. Advance was strongly associated with winter warming (not summer). The Alps alpine treeline has risen ~24 m/decade; Scandinavia mountain birch (Betula pubescens ssp. czerepanovii) treeline has climbed ~35–100 m over the past century. The diffusion-reaction model for treeline propagation is:
\[\frac{\partial N}{\partial t} = D\nabla^2 N + r\,N\bigl(1 - N/K(T,z)\bigr)\]
The Fisher wave speed is \(c = 2\sqrt{rD}\); for boreal conifers \(r\sim 0.05\)/yr, \(D\sim 0.01\) km²/yr gives \(c \sim 45\) m/yr—consistent with observed treeline velocities.
Biome Transitions: Holocene vs. Anthropocene
2. Oak Expansion vs. Beech Contraction
Svenning & Skov (2008, Ecology Letters) used palaeodata and range-filling metrics to show that European trees today occupy only 38–76% of their climatically suitable area: species are still catching up to the mid-Holocene optimum. Under business-as-usual warming, oaks (Quercus robur, Q. petraea) can expand northward by ∼500 km via pollen dispersal and acorn caching (jays are the dominant vector, with mean dispersal kernels 1–5 km per event), while beech (Fagus sylvatica) is constrained southward by increasing summer drought stress.
Cheaib et al. (2012) modelled French forests under A1B and showed beech declining by 30–50% in the Mediterranean margin by 2050. The key mechanism is \(c_i/c_a\)stomatal down-regulation that fails once vapour-pressure deficit exceeds 2 kPa for sustained periods. Beech's shallow root architecture makes it drought-vulnerable compared to deeply-rooted Turkey oak (Q. cerris) and pubescent oak (Q. pubescens) replacing it.
Niche Model
For a species, the niche is a function of the \(k\)-dimensional climate vector \(\mathbf{x}\). A Gaussian niche model is:
\[P(\text{present}\mid\mathbf{x}) = \exp\!\Bigl(-\tfrac12 (\mathbf{x}-\boldsymbol{\mu})^\top \Sigma^{-1}(\mathbf{x}-\boldsymbol{\mu})\Bigr)\]
Fitting via a logistic GLM with quadratic BIO1 and BIO12 terms yields predictions consistent with the Gaussian form. The simulation in this module uses this approach on CHELSA-like bioclim fields.
3. Boreal Fire Regimes & the Black-Carbon Feedback
Boreal Eurasia and North America have entered a new fire regime. The 2023 Canadian wildfire season burned 18.5 Mha—more than double the previous record—releasing approximately 4 Pg CO&sub2;e. Siberian “zombie fires” that overwinter in peat and reignite in spring have been documented in the Yakutian republic; Witze (2020) described their persistence via smouldering combustion at depth.
Soja et al. (2007, Global and Planetary Change) proposed a fire–climate feedback loop: warmer spring temperatures lengthen the fire weather season, black carbon deposited on the Arctic snowpack lowers albedo and accelerates melt, and permafrost degradation exposes peat carbon to combustion. The fire weather index (FWI) follows Canadian Forest Service formalism combining temperature, relative humidity, wind and precipitation:
\[\text{FWI} = f(\text{ISI},\text{BUI}), \quad \text{ISI}\propto U\,e^{-0.05\cdot\text{FFMC}^{-1}}\]
FFMC = fine fuel moisture code, ISI = initial spread index, BUI = build-up index.
Empirically, burn area scales exponentially with summer temperature anomaly (Balch et al., 2017): \(A_{\text{burn}}\propto e^{\beta\Delta T}\) with \(\beta \approx 0.6\)/°C in western North America. Under SSP5-8.5, boreal burn area is projected to double by 2100.
Black Carbon Deposition
Black carbon (soot) deposited on sea ice and snow reduces surface albedo by \(\Delta\alpha \approx -0.01\) to \(-0.08\) depending on concentration. This translates to a radiative forcing at the surface of:
\[\Delta F_{\text{BC}} = S_0(1-A_{\text{cloud}})\bigl|\Delta\alpha\bigr| \approx 2\text{--}15 \;\text{W/m}^2\]
on a local basis. Globally averaged, boreal fires contribute 0.05–0.25 W/m² (Bond et al., 2013), enough to matter in the Arctic radiative budget.
4. Moose, Caribou, Beetles
Moose (Alces alces) Heat Stress
Moose are large-bodied heat-retaining ungulates with a thermal neutral zone bounded above by ~14 °C in summer and ~−5 °C in winter (Renecker & Hudson, 1986). Above these thresholds, individuals experience increased respiration rate, reduced feeding time, and immunosuppression. The heat-stress hours per day:
\[H = \int_0^{24}\mathbb{1}\!\bigl[T(t) > T_{\text{crit}}\bigr]\,dt\]
Moose also carry increasing winter tick (Dermacentor albipictus) burdens with warmer autumns. Murray et al. (2006) documented “ghost moose” (hair loss from tick grooming) in New England, with calf mortality exceeding 70% in years with late snowfall. Southern range contractions are documented in Minnesota (50% decline since 2006), Nova Scotia, and New Hampshire.
Caribou & Plant-Phenology Mismatch
Post & Forchhammer (2008, Phil. Trans. R. Soc. B) showed that caribou (Rangifer tarandus) calf survival in West Greenland dropped as snowmelt advanced faster than the female migration schedule adjusted. Calves are born when energy demand peaks but now miss the peak of high-protein forage (graminoids, willow buds). The mismatch index \(\Delta\tau\) (days between calving and vegetation green-up) rose from ~−6 to +10 days between 1993 and 2007; recruitment fell proportionally.
Mountain Pine Beetle (Dendroctonus ponderosae)
The mountain pine beetle has expanded into previously uninhabitable high-elevation whitebark pine and into the boreal jack-pine zone via the Peace River corridor (Carroll et al., 2004; Cullingham et al., 2011). The limiting factor—winter temperatures below −40 °C that kill overwintering larvae—has retreated northward and upslope. Beetle generation time compresses from 2 years to 1 year when summer degree-days > 833 °C·d. The resulting outbreaks in British Columbia killed ~18 Mha of lodgepole pine by 2015, turning the biome from a net carbon sink to a net carbon source for over a decade (Kurz et al., 2008, Nature).
5. European Passerines: Phenological Mismatch
Visser & Both (2005, Proc. R. Soc. B) formalised the concept of phenological mismatch: a trophic level responds to warming with a different sensitivity from the one below, so demand and supply of food decouple. The canonical three-trophic system is:
- Oak bud-burst (Quercus robur): sensitivity ~−5 days / +1 °C spring T.
- Winter-moth caterpillar peak(Operophtera brumata): ~−5 days / +1 °C, tracking budburst.
- Pied flycatcher chick demand(Ficedula hypoleuca): ~−1 day / +1 °C (arrival in Europe constrained by Sahel rainfall, photoperiod, and intrinsic circannual clock).
The fitness penalty from mismatch follows a Gaussian around the caterpillar peak (\(\sigma_{\text{match}}\sim 7\) days):
\[W(\Delta) = W_{\max}\exp\!\Bigl(-\tfrac{\Delta^2}{2\sigma_{\text{match}}^2}\Bigr), \quad \Delta = t_{\text{chick}} - t_{\text{cat}}\]
Both et al. (2006, Nature) showed that pied flycatcher populations in mismatch-prone Dutch woodlands had declined by 90% since 1987, while those in mismatch-buffered habitats (early oaks + early arrival) held steady. The simulation in this module reproduces this dynamic.
Other Passerines
The common nightingale (Luscinia megarhynchos) breeding range has contracted northward in the UK but arrival dates are only modestly shifted. The willow warbler (Phylloscopus trochilus) shows declining productivity in southern England linked to shifting vegetation structure. The Arctic warbler (P. borealis) is extending its breeding range northward through Siberia by ∼80 km/decade (Saino et al., 2011).
Woodpecker Habitat
The three-toed woodpecker (Picoides tridactylus), an old-growth boreal specialist, is losing habitat to beetle-killed pine that no longer offers a sustained bark-beetle prey base after outbreak collapse. The white-backed woodpecker (Dendrocopos leucotos) in Scandinavian deciduous forests has declined >60% since 1975 (Swedish Species Information Centre, 2020).
6. Soil Carbon & Mycorrhizal Shifts (ECM vs AM)
Temperate and boreal soils store ~40% of terrestrial carbon. Their fate under warming is mediated by the mycorrhizal type of the dominant trees. Ectomycorrhizal (ECM) associations dominate in beech, oak, pine and spruce; arbuscular mycorrhizal (AM) associations dominate in maple, ash and cherry. Steidinger et al. (2019, Nature) mapped the global ECM/AM distribution and showed ECM forests store 1.7× more soil organic carbon than AM forests because ECM fungi mine organic-N directly, slowing litter decomposition.
With warming, AM trees expand northward and uphill faster, so the biome-integrated ECM:AM ratio is projected to fall. Extrapolating Steidinger's maps, a transition from ECM to AM dominance on 10% of the boreal zone would release an estimated 20–35 Pg C to the atmosphere.
Priming, Nitrogen, and the Microbial Carbon Pump
Warmer soils do not simply respire faster; the mix of microbial carbon-use efficiency (CUE), extracellular enzyme activity, and N mineralisation reshapes soil carbon stocks non-monotonically. The microbial carbon pump concept (Liang et al., 2017, Nature Microbiology) posits that microbially-derived compounds are major contributors to stable soil organic matter through the ex-vivo modification and in-vivo turnover pathways. Rising temperatures increase turnover but may also reduce CUE below ∼0.25, tilting balance toward CO&sub2; release.
Mass-balance for a single soil C pool \(C\):
\[\frac{dC}{dt} = I(T) - k(T,N)\,C, \quad k(T,N) = k_0\,Q_{10}^{\,(T-T_{\text{ref}})/10}\,(1 + \eta N)\]
Input \(I(T)\) rises with NPP; decomposition rate \(k\) rises with \(T\) and available \(N\). Nitrogen deposition in central Europe (10–30 kg N / ha / yr) can accelerate decomposition by 20–50%.
Decomposition Temperature Sensitivity
The Arrhenius / Q10 formalism for soil organic matter decay:
\[R(T) = R_{\text{ref}}\cdot Q_{10}^{(T-T_{\text{ref}})/10}, \quad Q_{10}\approx 2\text{--}3\]
Boreal soils typically have \(Q_{10}\approx 2.5\) at low temperatures, so a warming of 4 °C nearly doubles respiration (Davidson & Janssens, 2006). Whether this is matched by productivity gains (greening feedback) remains the central uncertainty in the boreal carbon balance.
7. Migratory Passerines & Sahel Connection
Afro-Palaearctic long-distance migrants face a double constraint: they must time arrival in Europe to the emerging spring food peak, yet their departure schedule from sub-Saharan Africa is fixed by photoperiod and Sahel rainfall anomalies (Both & Visser, 2001; Saino et al., 2011). The common nightingale, garden warbler, sedge warbler and willow warbler all suffer productivity penalties because their internal clock responds more slowly than European spring phenology.
The relative migration-phenology advance \(\Delta t_{\text{arr}}\) can be decomposed into a plastic response to en-route conditions and a genetic shift via selection on photoperiod thresholds (Pulido & Berthold 2010). Plastic response dominates on decadal timescales:
\[\Delta t_{\text{arr}} = h^2_{\text{mig}}\,S\,\Delta T + \beta_{\text{weather}}\,\Delta T_{\text{route}}\]
\(h^2_{\text{mig}}\sim 0.4\) narrow-sense heritability of arrival date; \(S\sim -0.3\) days per +1°C selection gradient; \(\beta_{\text{weather}}\sim -0.6\) days/°C en-route.
Sahel Rainfall Connectivity
Pre-breeding Sahel rainfall in West Africa determines insect food availability and departure body-condition for willow warblers, common whitethroats (Sylvia communis) and sand martins (Riparia riparia). The Sahel's 1970–80s drought and subsequent recovery explain a large fraction of European migratory-bird population trends (Sanderson et al., 2006). Future Sahel rainfall under SSP scenarios is a key but uncertain CMIP6 projection.
Arctic Warbler Range Expansion
The Arctic warbler (Phylloscopus borealis), which winters in Southeast Asia and breeds across Siberia and into western Alaska, is extending its breeding range northward by ~80 km/decade and upslope in the Altai by ~30 m/decade (Sokolov, 2006). Its arrival date has advanced by 5–7 days since 1970, but breeding productivity has declined at the southern range edge due to drought and habitat conversion.
Woodpecker Habitat & Deadwood
Old-growth boreal woodpecker species depend on deadwood for both foraging substrate (bark-boring beetle larvae) and nest cavities. Fennoscandian managed forests retain <2–5 m³/ha of deadwood against 40–100 m³/ha in reference old-growth (Siitonen, 2001). The white-backed woodpecker needs >20 m³/ha of coarse deciduous deadwood and is now Critically Endangered in Sweden. The three-toed woodpecker experiences a boom after beetle outbreaks and a bust as snags fall. Managing for continuity of dead-standing wood is the key conservation intervention.
Simulation 1: Oak & Beech SDM under CHELSA SSPs
We fit a logistic-GLM species-distribution model with quadratic BIO1/BIO12 terms to synthesised European presence-absence data, then project oak and beech suitability maps for present, SSP2-4.5 and SSP5-8.5 in 2080, reporting latitudinal range-centroid shifts and probability-mass retention.
Click Run to execute the Python code
Code will be executed with Python 3 on the server
8. Mismatch Quantification
For a predator (chick demand) and resource (caterpillar peak) with respective climate sensitivities \(b_{\text{pred}}\) and \(b_{\text{res}}\) days per degree, the mismatch under warming \(\Delta T\) is:
\[\Delta(\Delta T) = (b_{\text{res}} - b_{\text{pred}})\,\Delta T + (t_{\text{pred},0} - t_{\text{res},0})\]
For pied flycatcher vs. caterpillar, \(b_{\text{pred}}\approx 1\), \(b_{\text{res}}\approx 5\) days/°C, so +2 °C warming adds +8 days to the mismatch. Converting to a fitness-cost curve (Gaussian with \(\sigma\sim 7\) d) gives projected fledgling output. This is exactly the machinery of Simulation 2.
Demographic Consequences
For a shortlived songbird with adult survival \(s_a\approx 0.5\), fecundity of 4 fledglings per pair is replacement. Below ~2 per pair, populations decline deterministically. Our simulation captures this threshold.
Simulation 2: Pied Flycatcher Mismatch 1980-2075
We simulate caterpillar peak date, pied flycatcher arrival date, and chick-demand peak date under SSP2-4.5 and SSP5-8.5 warming, and translate the mismatch to fledglings per pair via a Gaussian fitness kernel. Output uses np.trapezoid for time-integrated productivity.
Click Run to execute the Python code
Code will be executed with Python 3 on the server
9. Conservation & Management Implications
Four management levers have demonstrated traction in temperate and boreal conservation under climate change:
- Assisted migration: Trans-locating provenances of oak and Douglas fir from warmer source regions to match projected future climates. Practised operationally in British Columbia (Aitken & Bemmels, 2016) and France's REINFFORCE network.
- Habitat heterogeneity: Retaining micro-refugia (north-facing slopes, riparian corridors, forest canopies > 25 m tall) that buffer understorey temperatures by 1–4 °C (De Frenne et al., 2019). Tree canopies act as living air-conditioners for herbaceous diversity.
- Fire management: Prescribed low-intensity burns combined with fire-break networks reduce future crown-fire risk. Indigenous burning traditions in Ontario, Alberta and Yukon are being revived (Christianson, 2015).
- Deadwood retention: Retaining 20–60 m³/ha of standing and fallen deadwood restores saproxylic beetle, fungal, and cavity-nesting bird communities (Gustafsson et al., 2012).
The European Union's Nature Restoration Law (2024) sets binding targets to restore 20% of degraded ecosystems by 2030 and all degraded ecosystems by 2050. For temperate and boreal biomes this translates to roughly 50 Mha of reforestation and re-wilding, with explicit climate-adaptive provenance selection.
Key References
• Aitken, S. N. & Bemmels, J. B. (2016). “Time to get moving: assisted gene flow of forest trees.” Evolutionary Applications, 9, 271–290.
• De Frenne, P. et al. (2019). “Global buffering of temperatures under forest canopies.” Nature Ecology & Evolution, 3, 744–749.
• Gustafsson, L. et al. (2012). “Retention forestry to maintain multifunctional forests.” BioScience, 62, 633–645.
• Siitonen, J. (2001). “Forest management, coarse woody debris and saproxylic organisms.” Ecological Bulletins, 49, 11–41.
• Sokolov, L. V. (2006). “Effect of global warming on the timing of migration and breeding of passerines in the 20th century.” Entomological Review, 86, S59–S81.
• Pulido, F. & Berthold, P. (2010). “Current selection for lower migratory activity will drive the evolution of residency in a migratory bird population.” PNAS, 107, 7341–7346.
• Sanderson, F. J. et al. (2006). “Long-term population declines in Afro-Palaearctic migrant birds.” Biological Conservation, 131, 93–105.
• Christianson, A. (2015). “Social science research on Indigenous wildfire management in the 21st century.” International Journal of Wildland Fire, 24, 190–200.
• Parmesan, C. & Yohe, G. (2003). “A globally coherent fingerprint of climate change impacts across natural systems.” Nature, 421, 37–42.
• Parmesan, C. (2006). “Ecological and evolutionary responses to recent climate change.” Annual Review of Ecology, Evolution, and Systematics, 37, 637–669.
• Chen, I.-C. et al. (2011). “Rapid range shifts of species associated with high levels of climate warming.” Science, 333, 1024–1026.
• Harsch, M. A. et al. (2009). “Are treelines advancing? A global meta-analysis of treeline response to climate warming.” Ecology Letters, 12, 1040–1049.
• Svenning, J.-C. & Skov, F. (2008). “Limited filling of the potential range in European tree species.” Ecology Letters, 7, 565–573.
• Cheaib, A. et al. (2012). “Climate change impacts on tree ranges: model intercomparison facilitates understanding and quantification of uncertainty.” Ecology Letters, 15, 533–544.
• Thuiller, W. et al. (2005). “Climate change threats to plant diversity in Europe.” Proc. Natl. Acad. Sci. USA, 102, 8245–8250.
• Soja, A. J. et al. (2007). “Climate-induced boreal forest change: Predictions versus current observations.” Global and Planetary Change, 56, 274–296.
• Balch, J. K. et al. (2017). “Human-started wildfires expand the fire niche across the United States.” Proc. Natl. Acad. Sci. USA, 114, 2946–2951.
• Bond, T. C. et al. (2013). “Bounding the role of black carbon in the climate system.” J. Geophys. Res. Atmos., 118, 5380–5552.
• Witze, A. (2020). “The Arctic is burning like never before.” Nature, 585, 336–337.
• Renecker, L. A. & Hudson, R. J. (1986). “Seasonal energy expenditures and thermoregulatory responses of moose.” Canadian Journal of Zoology, 64, 322–327.
• Murray, D. L. et al. (2006). “Pathogens, nutritional deficiency, and climate influences on a declining moose population.” Wildlife Monographs, 166, 1–30.
• Post, E. & Forchhammer, M. C. (2008). “Climate change reduces reproductive success of an Arctic herbivore through trophic mismatch.” Phil. Trans. R. Soc. B, 363, 2367–2373.
• Carroll, A. L. et al. (2004). “Effects of climate change on range expansion by the mountain pine beetle in British Columbia.” Mountain Pine Beetle Symposium Proceedings, 223–232.
• Cullingham, C. I. et al. (2011). “Mountain pine beetle host-range expansion threatens the boreal forest.” Molecular Ecology, 20, 2157–2171.
• Kurz, W. A. et al. (2008). “Mountain pine beetle and forest carbon feedback to climate change.” Nature, 452, 987–990.
• Visser, M. E. & Both, C. (2005). “Shifts in phenology due to global climate change: the need for a yardstick.” Proc. R. Soc. B, 272, 2561–2569.
• Both, C. et al. (2006). “Climate change and population declines in a long-distance migratory bird.” Nature, 441, 81–83.
• Reed, T. E. et al. (2013). “Phenological mismatch strongly affects individual fitness but not population demography in a woodland passerine.” Journal of Animal Ecology, 82, 131–144.
• Saino, N. et al. (2011). “Climate warming, ecological mismatch at arrival and population decline in migratory birds.” Proc. R. Soc. B, 278, 835–842.
• Steidinger, B. S. et al. (2019). “Climatic controls of decomposition drive the global biogeography of forest-tree symbioses.” Nature, 569, 404–408.
• Davidson, E. A. & Janssens, I. A. (2006). “Temperature sensitivity of soil carbon decomposition and feedbacks to climate change.” Nature, 440, 165–173.
• Karger, D. N. et al. (2017). “Climatologies at high resolution for the earth's land surface areas (CHELSA).” Scientific Data, 4, 170122.