Low-severity fires in the boreal region: reproductive implications for black spruce stands in between stand-replacing fire events.

BACKGROUND AND AIMS: Stand-replacing crown fires are the most prevalent type of fire regime in boreal forests in North America. However, a substantial proportion of low-severity fires are found within fire perimeters. Here we aimed to investigate the effects of low-severity fires on the reproductive potential and seedling recruitment in boreal forests stands in between stand-replacing fire events. METHODS: We recorded site and tree characteristics from 149 trees within twelve sites dominated by mature black spruce [Picea mariana (Mill.) B.S.P.] trees in the Northwest Territories, Canada. The presence of fire-scarred trees supported classification of sites as unburned or affected by low-severity fires in recent history. We used non-parametric tests to evaluate differences in site conditions between unburned and low-severity sites, and mixed effect models to evaluate differences in tree age, size, and reproductive traits among unburned trees and trees from low-severity sites. KEY RESULTS: Results showed significantly higher density of dead black spruce trees in low-severity sites, and marginally significant higher presence of permafrost. Trees from low-severity fire sites were significantly older, exhibited significantly lower tree growth, and showed a tendency towards a higher probability of cone presence and percentage of open cones compared to trees from unburned sites. Surviving fire-scarred trees affected by more recent low-severity fires showed a tendency towards higher probability of cone presence and cone production. Density of black spruce seedlings significantly increased with recent low-severity fires. CONCLUSIONS: Trees in low-severity sites appeared to have escaped mortality from up to three fires, as indicated by fire scar records and their older ages. Shallow permafrost at low-severity sites may cause lower flammability, allowing areas to act as fire refugia. Low-severity surface fires temporarily enhanced the reproductive capacity of surviving trees and the density of seedlings, likely as a stress response to the fire event.

Mycorrhizal feedbacks influence global forest structure and diversity.

One mechanism proposed to explain high species diversity in tropical systems is strong negative conspecific density dependence (CDD), which reduces recruitment of juveniles in proximity to conspecific adult plants. Although evidence shows that plant-specific soil pathogens can drive negative CDD, trees also form key mutualisms with mycorrhizal fungi, which may counteract these effects. Across 43 large-scale forest plots worldwide, we tested whether ectomycorrhizal tree species exhibit weaker negative CDD than arbuscular mycorrhizal tree species. We further tested for conmycorrhizal density dependence (CMDD) to test for benefit from shared mutualists. We found that the strength of CDD varies systematically with mycorrhizal type, with ectomycorrhizal tree species exhibiting higher sapling densities with increasing adult densities than arbuscular mycorrhizal tree species. Moreover, we found evidence of positive CMDD for tree species of both mycorrhizal types. Collectively, these findings indicate that mycorrhizal interactions likely play a foundational role in global forest diversity patterns and structure.

Permafrost thaw induces short-term increase in vegetation productivity in northwestern Canada.

Over the past several decades, various trends in vegetation productivity, from increases to decreases, have been observed throughout Arctic-Boreal ecosystems. While some of this variation can be explained by recent climate warming and increased disturbance, very little is known about the impacts of permafrost thaw on productivity across diverse vegetation communities. Active layer thickness data from 135 permafrost monitoring sites along a 10° latitudinal transect of the Northwest Territories, Canada, paired with a Landsat time series of normalized difference vegetation index from 1984 to 2019, were used to quantify the impacts of changing permafrost conditions on vegetation productivity. We found that active layer thickness contributed to the observed variation in vegetation productivity in recent decades in the northwestern Arctic-Boreal, with the highest rates of greening occurring at sites where the near-surface permafrost recently had thawed. However, the greening associated with permafrost thaw was not sustained after prolonged periods of thaw and appeared to diminish after the thaw front extended outside the plants' rooting zone. Highest rates of greening were found at the mid-transect sites, between 62.4° N and 65.2° N, suggesting that more southernly sites may have already surpassed the period of beneficial permafrost thaw, while more northern sites may have yet to reach a level of thaw that supports enhanced vegetation productivity. These results indicate that the response of vegetation productivity to permafrost thaw is highly dependent on the extent of active layer thickening and that increases in productivity may not continue in the coming decades.

Material Legacies and Environmental Constraints Underlie Fire Resilience of a Dominant Boreal Forest Type.

Resilience of plant communities to disturbance is supported by multiple mechanisms, including ecological legacies affecting propagule availability, species' environmental tolerances, and biotic interactions. Understanding the relative importance of these mechanisms for plant community resilience supports predictions of where and how resilience will be altered with disturbance. We tested mechanisms underlying resilience of forests dominated by black spruce (Picea mariana) to fire disturbance across a heterogeneous forest landscape in the Northwest Territories, Canada. We combined surveys of naturally regenerating seedlings at 219 burned plots with experimental manipulations of ecological legacies via seed addition of four tree species and vertebrate exclosures to limit granivory and herbivory at 30 plots varying in moisture and fire severity. Black spruce recovery was greatest where it dominated pre-fire, at wet sites with deep residual soil organic layers, and fire conditions of low soil or canopy combustion and longer return intervals. Experimental addition of seed indicated all species were seed-limited, emphasizing the importance of propagule legacies. Black spruce and birch (Betula papyrifera) recruitment were enhanced with vertebrate exclusion. Our combination of observational and experimental studies demonstrates black spruce is vulnerable to effects of increased fire activity that erode ecological legacies. Moreover, black spruce relies on wet areas with deep soil organic layers where other species are less competitive. However, other species can colonize these areas if enough seed is available or soil moisture is altered by climate change. Testing mechanisms underlying species' resilience to disturbance aids predictions of where vegetation will transform with effects of climate change. Supplementary Information: The online version contains supplementary material available at 10.1007/s10021-022-00772-7.

Understanding among-lake variability of mercury concentrations in Northern Pike (Esox lucius): A whole-ecosystem study in subarctic lakes.

Mercury concentrations ([Hg]) in fish reflect complex biogeochemical and ecological interactions that occur at a range of spatial and biological scales. Elucidating these interactions is crucial to understanding and predicting fish [Hg], particularly at northern latitudes, where environmental perturbations are having profound effects on land-water-animal interactions, and where fish are a critical subsistence food source. Using data from eleven subarctic lakes that span an area of ~60,000 km2 in the Dehcho Region of Northwest Territories (Canada), we investigated how trophic ecology and growth rates of fish, lake water chemistry, and catchment characteristics interact to affect [Hg] in Northern Pike (Esox lucius), a predatory fish of widespread subsistence and commercial importance. Results from linear regression and piecewise structural equation models showed that 83% of among-lake variability in Northern Pike [Hg] was explained by fish growth rates (negative) and concentrations of methyl Hg ([MeHg]) in benthic invertebrates (positive). These variables were in turn influenced by concentrations of dissolved organic carbon, MeHg (water), and total Hg (sediment) in lakes, which were ultimately driven by catchment characteristics. Lakes in relatively larger catchments and with more temperate/subpolar needleleaf and mixed forests had higher [Hg] in Northern Pike. Our results provide a plausible mechanistic understanding of how interacting processes at scales ranging from whole catchments to individual organisms influence fish [Hg], and give insight into factors that could be considered for prioritizing lakes for monitoring in subarctic regions.

Joint effects of climate, tree size, and year on annual tree growth derived from tree-ring records of ten globally distributed forests.

Tree rings provide an invaluable long-term record for understanding how climate and other drivers shape tree growth and forest productivity. However, conventional tree-ring analysis methods were not designed to simultaneously test effects of climate, tree size, and other drivers on individual growth. This has limited the potential to test ecologically relevant hypotheses on tree growth sensitivity to environmental drivers and their interactions with tree size. Here, we develop and apply a new method to simultaneously model nonlinear effects of primary climate drivers, reconstructed tree diameter at breast height (DBH), and calendar year in generalized least squares models that account for the temporal autocorrelation inherent to each individual tree's growth. We analyze data from 3811 trees representing 40 species at 10 globally distributed sites, showing that precipitation, temperature, DBH, and calendar year have additively, and often interactively, influenced annual growth over the past 120 years. Growth responses were predominantly positive to precipitation (usually over ≥3-month seasonal windows) and negative to temperature (usually maximum temperature, over ≤3-month seasonal windows), with concave-down responses in 63% of relationships. Climate sensitivity commonly varied with DBH (45% of cases tested), with larger trees usually more sensitive. Trends in ring width at small DBH were linked to the light environment under which trees established, but basal area or biomass increments consistently reached maxima at intermediate DBH. Accounting for climate and DBH, growth rate declined over time for 92% of species in secondary or disturbed stands, whereas growth trends were mixed in older forests. These trends were largely attributable to stand dynamics as cohorts and stands age, which remain challenging to disentangle from global change drivers. By providing a parsimonious approach for characterizing multiple interacting drivers of tree growth, our method reveals a more complete picture of the factors influencing growth than has previously been possible.

Increasing fire and the decline of fire adapted black spruce in the boreal forest.

Intensifying wildfire activity and climate change can drive rapid forest compositional shifts. In boreal North America, black spruce shapes forest flammability and depends on fire for regeneration. This relationship has helped black spruce maintain its dominance through much of the Holocene. However, with climate change and more frequent and severe fires, shifts away from black spruce dominance to broadleaf or pine species are emerging, with implications for ecosystem functions including carbon sequestration, water and energy fluxes, and wildlife habitat. Here, we predict that such reductions in black spruce after fire may already be widespread given current trends in climate and fire. To test this, we synthesize data from 1,538 field sites across boreal North America to evaluate compositional changes in tree species following 58 recent fires (1989 to 2014). While black spruce was resilient following most fires (62%), loss of resilience was common, and spruce regeneration failed completely in 18% of 1,140 black spruce sites. In contrast, postfire regeneration never failed in forests dominated by jack pine, which also possesses an aerial seed bank, or broad-leaved trees. More complete combustion of the soil organic layer, which often occurs in better-drained landscape positions and in dryer duff, promoted compositional changes throughout boreal North America. Forests in western North America, however, were more vulnerable to change due to greater long-term climate moisture deficits. While we find considerable remaining resilience in black spruce forests, predicted increases in climate moisture deficits and fire activity will erode this resilience, pushing the system toward a tipping point that has not been crossed in several thousand years.

Permafrost condition determines plant community composition and community-level foliar functional traits in a boreal peatland.

Boreal peatlands are critical ecosystems globally because they house 30%-40% of terrestrial carbon (C), much of which is stored in permafrost soil vulnerable to climate warming-induced thaw. Permafrost thaw leads to thickening of the active (seasonally thawed) layer and alters nutrient and light availability. These physical changes may influence community-level plant functional traits through intraspecific trait variation and/or species turnover. As permafrost thaw is expected to cause an efflux of carbon dioxide (CO2) and methane (CH4) from the soil to the atmosphere, it is important to understand thaw-induced changes in plant community productivity to evaluate whether these changes may offset some of the anticipated increases in C emissions. To this end, we collected vascular plant community composition and foliar functional trait data along gradients in aboveground tree biomass and active layer thickness (ALT) in a rapidly thawing boreal peatland, with the expectation that changes in above- and belowground conditions are indicative of altered resource availability. We aimed to determine whether community-level traits vary across these gradients, and whether these changes are dominated by intraspecific trait variation, species turnover, or both. Our results highlight that variability in community-level traits was largely attributable to species turnover and that both community composition and traits were predominantly driven by ALT. Specifically, thicker active layers associated with permafrost-free peatlands (i.e., bogs and fens) shifted community composition from slower-growing evergreen shrubs to faster-growing graminoids and forbs with a corresponding shift toward more productive trait values. The results from this rapidly thawing peatland suggest that continued warming-induced permafrost thaw and thermokarst development alter plant community composition and community-level traits and thus ecosystem productivity. Increased productivity may help to mitigate anticipated CO2 efflux from thawing permafrost, at least in the short term, though this response may be swamped by increase CH4 release.

Arbuscular mycorrhizal trees influence the latitudinal beta-diversity gradient of tree communities in forests worldwide.

Arbuscular mycorrhizal (AM) and ectomycorrhizal (EcM) associations are critical for host-tree performance. However, how mycorrhizal associations correlate with the latitudinal tree beta-diversity remains untested. Using a global dataset of 45 forest plots representing 2,804,270 trees across 3840 species, we test how AM and EcM trees contribute to total beta-diversity and its components (turnover and nestedness) of all trees. We find AM rather than EcM trees predominantly contribute to decreasing total beta-diversity and turnover and increasing nestedness with increasing latitude, probably because wide distributions of EcM trees do not generate strong compositional differences among localities. Environmental variables, especially temperature and precipitation, are strongly correlated with beta-diversity patterns for both AM trees and all trees rather than EcM trees. Results support our hypotheses that latitudinal beta-diversity patterns and environmental effects on these patterns are highly dependent on mycorrhizal types. Our findings highlight the importance of AM-dominated forests for conserving global forest biodiversity.

Unexpected greening in a boreal permafrost peatland undergoing forest loss is partially attributable to tree species turnover.

Time series of vegetation indices derived from satellite imagery are useful in measuring vegetation response to climate warming in remote northern regions. These indices show that productivity is generally declining in the boreal forest, but it is unclear which components of boreal vegetation are driving these trends. We aimed to compare trends in the normalized difference vegetation index (NDVI) to forest growth and demographic data taken from a 10 ha mapped plot located in a spruce-dominated boreal peatland. We used microcores to quantify recent growth trends and tree census data to characterize mortality and recruitment rates of the three dominant tree species. We then compared spatial patterns in growth and demography to patterns in Landsat-derived maximum NDVI trends (1984-2019) in 78 pixels that fell within the plot. We found that NDVI trends were predominantly positive (i.e., "greening") in spite of the ongoing loss of black spruce (the dominant species; 80% of stems) from the plot. The magnitude of these trends correlated positively with black spruce growth trends, but was also governed to a large extent by tree mortality and recruitment. Greening trends were weaker (lower slope) in areas with high larch mortality, and high turnover of spruce and birch, but stronger (higher slope) in areas with high larch recruitment. Larch dominance is currently low (~11% of stems), but it is increasing in abundance as permafrost thaw progresses and will likely have a substantial influence on future NDVI trends. Our results emphasize that NDVI trends in boreal peatlands can be positive even when the forest as a whole is in decline, and that the magnitude of trends can be strongly influenced by the demographics of uncommon species.

Seasonal thaw and landscape position determine foliar functional traits and whole-plant water use in tall shrubs on the low arctic tundra.

Climate warming is driving tundra shrub expansion with implications for ecosystem function and regional climate. Understanding associations between shrub ecophysiological function, distribution and environment is necessary for predicting consequences of expansion. We evaluated the role of topographic gradients on upland shrub productivity to understand potential constraints on shrub expansion. At a low arctic tundra site near Inuvik, Northwest Territories, Canada, we measured sap flow, stem water potential and productivity-related functional traits in green alder, and environmental predictors (water and nutrient availability and seasonal thaw depth) across a toposequence in alder patches. Seasonal thaw reduced stem sap flow whereas topographic position predicted stem water potential and productivity-related functional traits. Upslope shrubs were more water-limited than those downslope. Shrubs in drainage channels had traits associated with greater productivity than those on the tops of slopes. The effect of thaw depth on sap flow has implications for seasonal water-use patterns and warming impacts on tundra ecohydrology. Topographic variation in functional traits corresponds with observed spatial patterns of tundra shrub expansion along floodplains and concave hillslopes rather than in upland areas. Green alder is expanding rapidly across the low arctic tundra in northwestern North America; thus, anticipating the implications of its expansion is essential for predicting tundra function.

Increasing wildfires threaten historic carbon sink of boreal forest soils.

Boreal forest fires emit large amounts of carbon into the atmosphere primarily through the combustion of soil organic matter1-3. During each fire, a portion of this soil beneath the burned layer can escape combustion, leading to a net accumulation of carbon in forests over multiple fire events4. Climate warming and drying has led to more severe and frequent forest fires5-7, which threaten to shift the carbon balance of the boreal ecosystem from net accumulation to net loss1, resulting in a positive climate feedback8. This feedback will occur if organic-soil carbon that escaped burning in previous fires, termed 'legacy carbon', combusts. Here we use soil radiocarbon dating to quantitatively assess legacy carbon loss in the 2014 wildfires in the Northwest Territories of Canada2. We found no evidence for the combustion of legacy carbon in forests that were older than the historic fire-return interval of northwestern boreal forests9. In forests that were in dry landscapes and less than 60 years old at the time of the fire, legacy carbon that had escaped burning in the previous fire cycle was combusted. We estimate that 0.34 million hectares of young forests (<60 years) that burned in the 2014 fires could have experienced legacy carbon combustion. This implies a shift to a domain of carbon cycling in which these forests become a net source-instead of a sink-of carbon to the atmosphere over consecutive fires. As boreal wildfires continue to increase in size, frequency and intensity7, the area of young forests that experience legacy carbon combustion will probably increase and have a key role in shifting the boreal carbon balance.

Wildfire severity reduces richness and alters composition of soil fungal communities in boreal forests of western Canada.

Wildfire is the dominant disturbance in boreal forests and fire activity is increasing in these regions. Soil fungal communities are important for plant growth and nutrient cycling postfire but there is little understanding of how fires impact fungal communities across landscapes, fire severity gradients, and stand types in boreal forests. Understanding relationships between fungal community composition, particularly mycorrhizas, and understory plant composition is therefore important in predicting how future fire regimes may affect vegetation. We used an extreme wildfire event in boreal forests of Canada's Northwest Territories to test drivers of fungal communities and assess relationships with plant communities. We sampled soils from 39 plots 1 year after fire and 8 unburned plots. High-throughput sequencing (MiSeq, ITS) revealed 2,034 fungal operational taxonomic units. We found soil pH and fire severity (proportion soil organic layer combusted), and interactions between these drivers were important for fungal community structure (composition, richness, diversity, functional groups). Where fire severity was low, samples with low pH had higher total fungal, mycorrhizal, and saprotroph richness compared to where severity was high. Increased fire severity caused declines in richness of total fungi, mycorrhizas, and saprotrophs, and declines in diversity of total fungi and mycorrhizas. The importance of stand age (a surrogate for fire return interval) for fungal composition suggests we could detect long-term successional patterns even after fire. Mycorrhizal and plant community composition, richness, and diversity were weakly but significantly correlated. These weak relationships and the distribution of fungi across plots suggest that the underlying driver of fungal community structure is pH, which is modified by fire severity. This study shows the importance of edaphic factors in determining fungal community structure at large scales, but suggests these patterns are mediated by interactions between fire and forest stand composition.

Quantification of uncertainties in conifer sap flow measured with the thermal dissipation method.

Trees play a key role in the global hydrological cycle and measurements performed with the thermal dissipation method (TDM) have been crucial in providing whole-tree water-use estimates. Yet, different data processing to calculate whole-tree water use encapsulates uncertainties that have not been systematically assessed. We quantified uncertainties in conifer sap flux density (Fd ) and stand water use caused by commonly applied methods for deriving zero-flow conditions, dampening and sensor calibration. Their contribution has been assessed using a stem segment calibration experiment and 4 yr of TDM measurements in Picea abies and Larix decidua growing in contrasting environments. Uncertainties were then projected on TDM data from different conifers across the northern hemisphere. Commonly applied methods mostly underestimated absolute Fd . Lacking a site- and species-specific calibrations reduced our stand water-use measurements by 37% and induced uncertainty in northern hemisphere Fd . Additionally, although the interdaily variability was maintained, disregarding dampening and/or applying zero-flow conditions that ignored night-time water use reduced the correlation between environment and Fd . The presented ensemble of calibration curves and proposed dampening correction, together with the systematic quantification of data-processing uncertainties, provide crucial steps in improving whole-tree water-use estimates across spatial and temporal scales.

Boreal tree hydrodynamics: asynchronous, diverging, yet complementary.

Water stress has been identified as a key mechanism of the contemporary increase in tree mortality rates in northwestern North America. However, a detailed analysis of boreal tree hydrodynamics and their interspecific differences is still lacking. Here we examine the hydraulic behaviour of co-occurring larch (Larix laricina) and black spruce (Picea mariana), two characteristic boreal tree species, near the southern limit of the boreal ecozone in central Canada. Sap flux density (Js), concurrently recorded stem radius fluctuations and meteorological conditions are used to quantify tree hydraulic functioning and to scrutinize tree water-use strategies. Our analysis revealed asynchrony in the diel hydrodynamics of the two species with the initial rise in Js occurring 2 h earlier in larch than in black spruce. Interspecific differences in larch and black spruce crown architecture explained the observed asynchrony in their hydraulic functioning. Furthermore, the two species exhibited diverging stomatal regulation strategies with larch and black spruce employing relatively isohydric and anisohydric behaviour, respectively. Such asynchronous and diverging tree-level hydrodynamics provide new insights into the ecosystem-level complementarity in tree form and function, with implications for understanding boreal forests' water and carbon dynamics and their resilience to environmental stress.

Cross-scale controls on carbon emissions from boreal forest megafires.

Climate warming and drying is associated with increased wildfire disturbance and the emergence of megafires in North American boreal forests. Changes to the fire regime are expected to strongly increase combustion emissions of carbon (C) which could alter regional C balance and positively feedback to climate warming. In order to accurately estimate C emissions and thereby better predict future climate feedbacks, there is a need to understand the major sources of heterogeneity that impact C emissions at different scales. Here, we examined 211 field plots in boreal forests dominated by black spruce (Picea mariana) or jack pine (Pinus banksiana) of the Northwest Territories (NWT), Canada after an unprecedentedly large area burned in 2014. We assessed both aboveground and soil organic layer (SOL) combustion, with the goal of determining the major drivers in total C emissions, as well as to develop a high spatial resolution model to scale emissions in a relatively understudied region of the boreal forest. On average, 3.35 kg C m-2 was combusted and almost 90% of this was from SOL combustion. Our results indicate that black spruce stands located at landscape positions with intermediate drainage contribute the most to C emissions. Indices associated with fire weather and date of burn did not impact emissions, which we attribute to the extreme fire weather over a short period of time. Using these results, we estimated a total of 94.3 Tg C emitted from 2.85 Mha of burned area across the entire 2014 NWT fire complex, which offsets almost 50% of mean annual net ecosystem production in terrestrial ecosystems of Canada. Our study also highlights the need for fine-scale estimates of burned area that represent small water bodies and regionally specific calibrations of combustion that account for spatial heterogeneity in order to accurately model emissions at the continental scale.

Assessing local adaptation vs. plasticity under different resource conditions in seedlings of a dominant boreal tree species.

Under changing climate conditions, understanding local adaptation of plants is crucial to predicting the resilience of ecosystems. We selected black spruce (Picea mariana), the most dominant tree species in the North American boreal forest, in order to evaluate local adaptation vs. plasticity across regions experiencing some of the most extreme climate warming globally. Seeds from three provenances across the latitudinal extent of this species in northwestern Canada were planted in a common garden study in growth chambers. Two levels of two resource conditions were applied (low/high nutrient and ambient/elevated CO2) in a fully factorial design and we measured physiological traits, allocational traits, growth and survival. We found significant differences in height, root length and biomass among populations, with southern populations producing the largest seedlings. However, we did not detect meaningful significant differences among nutrient or CO2 treatments in any traits measured, and there were no consistent population-level differences in physiological traits or allocation patterns. We found that there was greater mortality after simulated winter in the high nutrient treatment, which may reflect an important shift in seedling growth strategies under increased resource availability. Our study provides important insight into how this dominant boreal tree species might respond to the changing climate conditions predicted in this region.

The world's biomes and primary production as a triple tragedy of the commons foraging game played among plants.

Plants appear to produce an excess of leaves, stems and roots beyond what would provide the most efficient harvest of available resources. One way to understand this overproduction of tissues is that excess tissue production provides a competitive advantage. Game theoretic models predict overproduction of all tissues compared with non-game theoretic models because they explicitly account for this indirect competitive benefit. Here, we present a simple game theoretic model of plants simultaneously competing to harvest carbon and nitrogen. In the model, a plant's fitness is influenced by its own leaf, stem and root production, and the tissue production of others, which produces a triple tragedy of the commons. Our model predicts (i) absolute net primary production when compared with two independent global datasets; (ii) the allocation relationships to leaf, stem and root tissues in one dataset; (iii) the global distribution of biome types and the plant functional types found within each biome; and (iv) ecosystem responses to nitrogen or carbon fertilization. Our game theoretic approach removes the need to define allocation or vegetation type a priori but instead lets these emerge from the model as evolutionarily stable strategies. We believe this to be the simplest possible model that can describe plant production.

Decreased competitive interactions drive a reverse species richness latitudinal gradient in subarctic forests.

The tendency for species richness to decrease toward the poles is one of the best-characterized patterns in biogeography. The mechanisms behind this pattern have received much attention, yet very few studies have investigated very high-latitude communities. Here, using data from 134 permanent sample plots from 60 degrees to 68 degrees N, we show that boreal forest plant communities in northwestern Canada increase in richness toward the poles, despite a strong increase in climatic harshness. We hypothesized three possible explanations for this pattern: (1) historical biogeography, (2) reduced competition for light at high latitudes (biotic interactions), and (3) changes in soil characteristics with latitude. We used multidimensional scaling to investigate the community composition at each site and found no clustering of communities by latitude, suggesting that historical biogeography was not constraining site diversity. We then investigated the mechanisms behind this gradient using both abiotic (climate and soil) and biotic (tree stand characteristics) variables in a multiple factor analysis. We found that the best predictor of species richness is an environmental gradient that describes an inverse relationship between temperature and tree-stand density, suggesting that reduced competition for light due to reduced tree growth at low temperatures at higher latitudes allows greater species richness. This study shows that low energy availability and climatic harshness may not be limiting species richness toward the poles, rather, abiotic effects act instead on the strength of biotic interactions.