Alternative measures of trait-niche relationships: A test on dispersal traits in saproxylic beetles.

Functional trait approaches are common in ecology, but a lack of clear hypotheses on how traits relate to environmental gradients (i.e., trait-niche relationships) often makes uncovering mechanisms difficult. Furthermore, measures of community functional structure differ in their implications, yet inferences are seldom compared among metrics. Community-weighted mean trait values (CWMs), a common measure, are largely driven by the most common species and thus do not reflect community-wide trait-niche relationships per se. Alternatively, trait-niche relationships can be estimated across a larger group of species using hierarchical joint species distribution models (JSDMs), quantified by a parameter Γ. We investigated how inferences about trait-niche relationships are affected by the choice of metric. Using deadwood-dependent (saproxylic) beetles in fragmented Finnish forests, we followed a protocol for investigating trait-niche relationships by (1) identifying environmental filters (climate, forest age, and deadwood volume), (2) relating these to an ecological function (dispersal ability), and (3) identifying traits related to this function (wing morphology). We tested 18 hypothesized dispersal relationships using both CWM and Γ estimates across these environmental gradients. CWMs were more likely than Γ to show support for trait-niche relationships. Up to 13% of species' realized niches were explained by dispersal traits, but the directions of effects were consistent with fewer than 11%-39% of our 18 trait-niche hypotheses (depending on the metric used). This highlights the difficulty in connecting morphological traits and ecological functions in insects, despite the clear conceptual link between landscape connectivity and flight-related traits. Caution is thus warranted in hypothesis development, particularly where apparent trait-function links are less clear. Inferences differ when CWMs versus Γ estimates are used, necessitating the choice of a metric that reflects study questions. CWMs help explain the effects of environmental gradients on community trait composition, whereas the effects of traits on species' niches are better estimated using hierarchical JSDMs.

High-resolution 3D forest structure explains ecomorphological trait variation in assemblages of saproxylic beetles.

Climate, topography and the 3D structure of forests are major drivers affecting local species communities. However, little is known about how the specific functional traits of saproxylic (wood-living) beetles, involved in the recycling of wood, might be affected by those environmental characteristics.Here, we combine ecological and morphological traits available for saproxylic beetles and airborne laser scanning (ALS) data in Bayesian trait-based joint species distribution models to study how traits drive the distributions of more than 230 species in temperate forests of Europe.We found that elevation (as a proxy for temperature and precipitation) and the proportion of conifers played important roles in species occurrences while variables related to habitat heterogeneity and forest complexity were less relevant. Furthermore, we showed that local communities were shaped by environmental variation primarily through their ecological traits whereas morphological traits were involved only marginally. As predicted, ecological traits influenced species' responses to forest structure, and to other environmental variation, with canopy niche, wood decay niche and host preference as the most important ecological traits. Conversely, no links between morphological traits and environmental characteristics were observed. Both models, however, revealed strong phylogenetic signal in species' response to environmental characteristics.These findings imply that alterations of climate and tree species composition have the potential to alter saproxylic beetle communities in temperate forests. Additionally, ecological traits help explain species' responses to environmental characteristics and thus should prove useful in predicting their responses to future change. It remains challenging, however, to link simple morphological traits to species' complex ecological niches. Read the free Plain Language Summary for this article on the Journal blog.

Future supply of boreal forest ecosystem services is driven by management rather than by climate change.

Forests provide a wide variety of ecosystem services (ES) to society. The boreal biome is experiencing the highest rates of warming on the planet and increasing demand for forest products. To foresee how to maximize the adaptation of boreal forests to future warmer conditions and growing demands of forest products, we need a better understanding of the relative importance of forest management and climate change on the supply of ecosystem services. Here, using Finland as a boreal forest case study, we assessed the potential supply of a wide range of ES (timber, bilberry, cowberry, mushrooms, carbon storage, scenic beauty, species habitat availability and deadwood) given seven management regimes and four climate change scenarios. We used the forest simulator SIMO to project forest dynamics for 100 years into the future (2016-2116) and estimate the potential supply of each service using published models. Then, we tested the relative importance of management and climate change as drivers of the future supply of these services using generalized linear mixed models. Our results show that the effects of management on the future supply of these ES were, on average, 11 times higher than the effects of climate change across all services, but greatly differed among them (from 0.53 to 24 times higher for timber and cowberry, respectively). Notably, the importance of these drivers substantially differed among biogeographical zones within the boreal biome. The effects of climate change were 1.6 times higher in northern Finland than in southern Finland, whereas the effects of management were the opposite-they were three times higher in the south compared to the north. We conclude that new guidelines for adapting forests to global change should account for regional differences and the variation in the effects of climate change and management on different forest ES.

More future synergies and less trade-offs between forest ecosystem services with natural climate solutions instead of bioeconomy solutions.

To reach the Paris Agreement, societies need to increase the global terrestrial carbon sink. There are many climate change mitigation solutions (CCMS) for forests, including increasing bioenergy, bioeconomy, and protection. Bioenergy and bioeconomy solutions use climate-smart, intensive management to generate high quantities of bioenergy and bioproducts. Protection of (semi-)natural forests is a major component of "natural climate solution" (NCS) since forests store carbon in standing biomass and soil. Furthermore, protected forests provide more habitat for biodiversity and non-wood ecosystem services (ES). We investigated the impacts of different CCMS and climate scenarios, jointly or in isolation, on future wood ES, non-wood ES, and regulating ES for a major wood provider for the international market. Specifically, we projected future ES given by three CCMS scenarios for Sweden 2020-2100. In the long term, fulfilling the increasing wood demand through bioenergy and bioeconomy solutions will decrease ES multifunctionality, but the increased stand age and wood stocks induced by rising greenhouse gas (GHG) concentrations will partially offset these negative effects. Adopting bioenergy and bioeconomy solutions will have a greater negative impact on ES supply than adopting NCS. Bioenergy or bioeconomy solutions, as well as increasing GHG emissions, will reduce synergies and increase trade-offs in ES. NCS, by contrast, increases the supply of multiple ES in synergy, even transforming current ES trade-offs into future synergies. Moreover, NCS can be considered an adaptation measure to offset negative climate change effects on the future supplies of non-wood ES. In boreal countries around the world, forestry strategies that integrate NCS more deeply are crucial to ensure a synergistic supply of multiple ES.

Occupancy versus colonization-extinction models for projecting population trends at different spatial scales.

Understanding spatiotemporal population trends and their drivers is a key aim in population ecology. We further need to be able to predict how the dynamics and sizes of populations are affected in the long term by changing landscapes and climate. However, predictions of future population trends are sensitive to a range of modeling assumptions. Deadwood-dependent fungi are an excellent system for testing the performance of different predictive models of sessile species as these species have different rarity and spatial population dynamics, the populations are structured at different spatial scales, and they utilize distinct substrates. We tested how the projected large-scale occupancies of species with differing landscape-scale occupancies are affected over the coming century by different modeling assumptions. We compared projections based on occupancy models against colonization-extinction models, conducting the modeling at alternative spatial scales and using fine- or coarse-resolution deadwood data. We also tested effects of key explanatory variables on species occurrence and colonization-extinction dynamics. The hierarchical Bayesian models applied were fitted to an extensive repeated survey of deadwood and fungi at 174 patches. We projected higher occurrence probabilities and more positive trends using the occupancy models compared to the colonization-extinction models, with greater difference for the species with lower occupancy, colonization rate, and colonization:extinction ratio than for the species with higher estimates of these statistics. The magnitude of future increase in occupancy depended strongly on the spatial modeling scale and resource resolution. We encourage using colonization-extinction models over occupancy models, modeling the process at the finest resource-unit resolution that is utilizable by the species, and conducting projections for the same spatial scale and resource resolution at which the model fitting is conducted. Further, the models applied should include key variables driving the metapopulation dynamics, such as the availability of suitable resource units, habitat quality, and spatial connectivity.

Levels of forest ecosystem services depend on specific mixtures of commercial tree species.

Global and local ecosystem change resulting in diversity loss has motivated efforts to understand relationships between species diversity and ecosystem services. However, it is unclear how such a general understanding can inform policies for the management of ecosystem services in production systems, because these systems are primarily used for food or fibre, and are rarely managed for the conservation of species diversity. Here, using data from a nationwide forest inventory covering an area of 230,000 km2, we show that relative abundances of commercial tree species in mixed stands strongly influence the potential to provide ecosystem services. The mixes provided higher levels of ecosystem services compared to respective plant monocultures (overyielding or transgressive overyielding) in 35% of the investigated cases, and lower (underyielding) in 9% of the cases. We further show that relative abundances, not just species richness per se, of specific tree-species mixtures affect the potential of forests to provide multiple ecosystem services, which is crucial information for policy and sustainable forest management.

Discovery of long-distance gamete dispersal in a lichen-forming ascomycete.

Accurate estimates of gamete and offspring dispersal range are required for the understanding and prediction of spatial population dynamics and species persistence. Little is known about gamete dispersal in fungi, especially in lichen-forming ascomycetes. Here, we estimate the dispersal functions of clonal propagules, gametes and ascospores of the epiphytic lichen Lobaria pulmonaria. We use hierarchical Bayesian parentage analysis, which integrates genetic and ecological information from multiannual colonization and dispersal source data collected in a large, old-growth forest landscape. The effective dispersal range of gametes is several hundred metres to kilometres from potential paternal individuals. By contrast, clonal propagules disperse only tens of metres, and ascospores disperse over several thousand metres. Our study reveals the dispersal distances of individual reproductive units; clonal propagules, gametes and ascospores, which is of great importance for a thorough understanding of the spatial dynamics of ascomycetes. Sexual reproduction occurs between distant individuals. However, whereas gametes and ascospores disperse over long distances, the overall rate of colonization of trees is low. Hence, establishment is the limiting factor for the colonization of new host trees by the lichen in old-growth landscapes.

Impacts of forestry on boreal forests: An ecosystem services perspective.

Forests are widely recognized as major providers of ecosystem services, including timber, other forest products, recreation, regulation of water, soil and air quality, and climate change mitigation. Extensive tracts of boreal forests are actively managed for timber production, but actions aimed at increasing timber yields also affect other forest functions and services. Here, we present an overview of the environmental impacts of forest management from the perspective of ecosystem services. We show how prevailing forestry practices may have substantial but diverse effects on the various ecosystem services provided by boreal forests. Several aspects of these processes remain poorly known and warrant a greater role in future studies, including the role of community structure. Conflicts among different interests related to boreal forests are most likely to occur, but the concept of ecosystem services may provide a useful framework for identifying and resolving these conflicts.

Forest management could counteract distribution retractions forced by climate change.

Climate change is expected to drive the distribution retraction of northern species. However, particularly in regions with a history of intensive exploitation, changes in habitat management could facilitate distribution expansions counter to expectations under climate change. Here, we test the potential for future forest management to facilitate the southward expansion of an old-forest species from the boreal region into the boreo-nemoral region, contrary to expectations under climate change. We used an ensemble of species distribution models based on citizen science data to project the response of Phellinus ferrugineofuscus, a red-listed old-growth indicator, wood-decaying fungus, to six forest management and climate change scenarios. We projected change in habitat suitability across the boreal and boreo-nemoral regions of Sweden for the period 2020-2100. Scenarios varied in the proportion of forest set aside from production, the level of timber extraction, and the magnitude of climate change. Habitat suitabilities for the study species were projected to show larger relative increases over time in the boreo-nemoral region compared to the boreal region, under all scenarios. By 2100, mean suitabilities in set-aside forest in the boreo-nemoral region were similar to the suitabilities projected for set-aside forest in the boreal region in 2020, suggesting that occurrence in the boreo-nemoral region could be increased. However, across all scenarios, consistently higher projected suitabilities in set-aside forest in the boreal region indicated that the boreal region remained the species stronghold. Furthermore, negative effects of climate change were evident in the boreal region, and projections suggested that climatic changes may eventually counteract the positive effects of forest management in the boreo-nemoral region. Our results suggest that the current rarity of this old-growth indicator species in the boreo-nemoral region may be due to the history of intensive forestry. Forest management therefore has the potential to compensate for the negative effects of climate change. However, increased occurrence at the southern range edge would depend on the dispersal and colonization ability of the species. An increase in the amount of set-aside forest across both the boreal and boreo-nemoral regions is therefore likely to be required to prevent the decline of old-forest species under climate change.

Evaluating citizen science data for forecasting species responses to national forest management.

The extensive spatial and temporal coverage of many citizen science datasets (CSD) makes them appealing for use in species distribution modeling and forecasting. However, a frequent limitation is the inability to validate results. Here, we aim to assess the reliability of CSD for forecasting species occurrence in response to national forest management projections (representing 160,366 km2) by comparison against forecasts from a model based on systematically collected colonization-extinction data. We fitted species distribution models using citizen science observations of an old-forest indicator fungus Phellinus ferrugineofuscus. We applied five modeling approaches (generalized linear model, Poisson process model, Bayesian occupancy model, and two MaxEnt models). Models were used to forecast changes in occurrence in response to national forest management for 2020-2110. Forecasts of species occurrence from models based on CSD were congruent with forecasts made using the colonization-extinction model based on systematically collected data, although different modeling methods indicated different levels of change. All models projected increased occurrence in set-aside forest from 2020 to 2110: the projected increase varied between 125% and 195% among models based on CSD, in comparison with an increase of 129% according to the colonization-extinction model. All but one model based on CSD projected a decline in production forest, which varied between 11% and 49%, compared to a decline of 41% using the colonization-extinction model. All models thus highlighted the importance of protected old forest for P. ferrugineofuscus persistence. We conclude that models based on CSD can reproduce forecasts from models based on systematically collected colonization-extinction data and so lead to the same forest management conclusions. Our results show that the use of a suite of models allows CSD to be reliably applied to land management and conservation decision making, demonstrating that widely available CSD can be a valuable forecasting resource.

Local epiphyte establishment and future metapopulation dynamics in landscapes with different spatiotemporal properties.

Understanding the relative importance of different ecological processes on the metapopulation dynamics of species is the basis for accurately forecasting metapopulation size in fragmented landscapes. Successful local colonization depends on both species dispersal range and how local habitat conditions affect establishment success. Moreover, there is limited understanding of the effects of different spatiotemporal landscape properties on future metapopulation size. We investigate which factors drive the future metapopulation size of the epiphytic model lichen species Lobaria pulmonaria in a managed forest landscape. First, we test the importance of dispersal and local conditions on the colonization-extinction dynamics of the species using Bayesian state-space modelling of a large-scale data set collected over a 10-yr period. Second, we test the importance of dispersal and establishment limitation in explaining establishment probability and subsequent local population growth, based on a 10-yr propagule sowing experiment. Third, we test how future metapopulation size is affected by different metapopulation and spatiotemporal landscape dynamics, using simulations with the metapopulation models fitted to the empirical data. The colonization probability increased with tree inclination and connectivity, with a mean dispersal distance of 97 m (95% credible intervals, 5-530 m). Local extinctions were mainly deterministic set by tree mortality, but also by tree cutting by forestry. No experimental establishments took place on clearcuts, and in closed forest the establishment probability was higher on trees growing on moist than on dry-mesic soils. The subsequent local population growth rate increased with increasing bark roughness. The simulations showed that the restricted dispersal range estimated (compared to non-restricted dispersal range), and short tree rotation length (65 yr instead of 120) had approximately the same negative effects on future metapopulation size, while regeneration of trees creating a random tree pattern instead of an aggregated one had only some negative effect. However, using the colonization rate obtained with the experimentally added diaspores led to a considerable increase in metapopulation size, making the dispersal limitation of the species clear. The future metapopulation size is thus set by the number of host trees located in shady conditions, not isolated from occupied trees, and by the rotation length of these host trees.

Dynamic anthropogenic edge effects on the distribution and diversity of fungi in fragmented old-growth forests.

Diversity patterns and dynamics at forest edges are not well understood. We disentangle the relative importance of edge-effect variables on spatio-temporal patterns in species richness and occupancy of deadwood-dwelling fungi in fragmented old-growth forests. We related richness and log occupancy by 10 old-growth forest indicator fungi and by two common fungi to log conditions in natural and anthropogenic edge habitats of 31 old-growth Picea abies forest stands in central Sweden. We compared edge-to-interior gradients (100 m) to the forest interior (beyond 100 m), and we analyzed stand-level changes after 10 yr. Both richness and occupancy of logs by indicator species was negatively related to adjacent young clear-cut edges, but this effect decreased with increasing clear-cut age. The occupancy of logs by indicator species also increased with increasing distance to the natural edges. In contrast, the occupancy of logs by common species was positively related or unrelated to distance to clear-cut edges regardless of the edge age, and this was partly explained by fungal specificity to substrate quality. Stand-level mean richness and mean occupancy of logs did not change for indicator or common species over a decade. By illustrating the importance of spatial and temporal dimensions of edge effects, we extend the general understanding of the distribution and diversity of substrate-confined fungi in fragmented old-growth forests. Our results highlight the importance of longer forest rotation times adjacent to small protected areas and forest set-asides, where it may take more than 50 yr for indicator species richness levels to recover to occupancy levels observed in the forest interior. Also, non-simultaneous clear-cutting of surrounding productive forests in a way that reduces the edge effect over time (i.e., dynamic buffers) may increase the effective core area of small forest set-asides and improve their performance on protecting species of special concern for conservation.

Green Infrastructure Design Based on Spatial Conservation Prioritization and Modeling of Biodiversity Features and Ecosystem Services.

There is high-level political support for the use of green infrastructure (GI) across Europe, to maintain viable populations and to provide ecosystem services (ES). Even though GI is inherently a spatial concept, the modern tools for spatial planning have not been recognized, such as in the recent European Environment Agency (EEA) report. We outline a toolbox of methods useful for GI design that explicitly accounts for biodiversity and ES. Data on species occurrence, habitats, and environmental variables are increasingly available via open-access internet platforms. Such data can be synthesized by statistical species distribution modeling, producing maps of biodiversity features. These, together with maps of ES, can form the basis for GI design. We argue that spatial conservation prioritization (SCP) methods are effective tools for GI design, as the overall SCP goal is cost-effective allocation of conservation efforts. Corridors are currently promoted by the EEA as the means for implementing GI design, but they typically target the needs of only a subset of the regional species pool. SCP methods would help to ensure that GI provides a balanced solution for the requirements of many biodiversity features (e.g., species, habitat types) and ES simultaneously in a cost-effective manner. Such tools are necessary to make GI into an operational concept for combating biodiversity loss and promoting ES.

Evidence-based knowledge versus negotiated indicators for assessment of ecological sustainability: the Swedish Forest Stewardship Council standard as a case study.

Assessing ecological sustainability involves monitoring of indicators and comparison of their states with performance targets that are deemed sustainable. First, a normative model was developed centered on evidence-based knowledge about (a) forest composition, structure, and function at multiple scales, and (b) performance targets derived by quantifying the habitat amount in naturally dynamic forests, and as required for presence of populations of specialized focal species. Second, we compared the Forest Stewardship Council (FSC) certification standards' ecological indicators from 1998 and 2010 in Sweden to the normative model using a Specific, Measurable, Accurate, Realistic, and Timebound (SMART) indicator approach. Indicator variables and targets for riparian and aquatic ecosystems were clearly under-represented compared to terrestrial ones. FSC's ecological indicators expanded over time from composition and structure towards function, and from finer to coarser spatial scales. However, SMART indicators were few. Moreover, they poorly reflected quantitative evidence-based knowledge, a consequence of the fact that forest certification mirrors the outcome of a complex social negotiation process.

Higher levels of multiple ecosystem services are found in forests with more tree species.

Forests are of major importance to human society, contributing several crucial ecosystem services. Biodiversity is suggested to positively influence multiple services but evidence from natural systems at scales relevant to management is scarce. Here, across a scale of 400,000 km(2), we report that tree species richness in production forests shows positive to positively hump-shaped relationships with multiple ecosystem services. These include production of tree biomass, soil carbon storage, berry production and game production potential. For example, biomass production was approximately 50% greater with five than with one tree species. In addition, we show positive relationships between tree species richness and proxies for other biodiversity components. Importantly, no single tree species was able to promote all services, and some services were negatively correlated to each other. Management of production forests will therefore benefit from considering multiple tree species to sustain the full range of benefits that the society obtains from forests.

Estimates of connectivity reveal non-equilibrium epiphyte occurrence patterns almost 180 years after habitat decline.

Habitat loss is a major cause of species decline and extinction. Immediately after habitat loss, species occurrences are not in equilibrium with the new landscape and more closely reflect the previous landscape structure. Species with slow colonisation-extinction dynamics may display long time-lags before reaching a new equilibrium. We investigated the importance of connectivity to current and historical dispersal sources with the aim of explaining the occurrence pattern of epiphytic lichens with different traits among 104 old oaks. We used oak survey data collected from 1830 and 2009 for a Swedish landscape where oak densities declined drastically shortly after 1830. We fitted a commonly used connectivity measure and estimated the confidence interval for the spatial scale parameter. Small differences in the spatial scale parameter resulted in large differences in model fit. Connectivity to trees in 1830 better explained the occurrence of three of the four species compared to the connectivity in 2009. The explanatory power of the historical landscape structure was highest for the species with traits that may result in a low colonisation rate--both a narrow niche (here few suitable trees) and large dispersal propagules. The results suggest that oak-dependent epiphytic lichens have not reached equilibrium with the spatial landscape structure 180 years after the drastic decline in habitat. For the long-term persistence of epiphytes associated with old trees, conservation efforts should focus on (1) protecting and restoring stands where specialised species with large dispersal propagules (i.e. with low colonisation rates) occur today and (2) promoting tree regeneration in their near vicinity.

Metapopulation dynamics and future persistence of epiphytic cyanolichens in a European boreal forest ecosystem.

1.: One approach to biodiversity conservation is to set aside small woodland key habitats (WKHs) in intensively managed landscapes. The aim is to support species, such as epiphytes, which often depend on old trees and are negatively affected by intensive forestry. However, it is not known whether the number of host trees within these areas can sustain species in the long term. 2.: We studied metapopulation dynamics and assessed the future persistence of epiphytes assuming host tree numbers similar to those observed in large north European WKHs. The study species were seven cyanolichens confined to Populus tremula in the boreal study area. Colonizations and extinctions were recorded in 2008 on trees that had been surveyed 13 years earlier. We applied generalized (non)linear models to test the importance of environmental conditions, facilitation and spatial connectivity on the metapopulation dynamics. We also simulated the effects of tree numbers and tree fall rates on future species persistence. 3.: Metapopulation dynamics were explained by tree quality, size or tree fall. In one species, colonizations increased with increasing connectivity, and in a second species it increased if other lichens sharing the photobiont with the focal species were present, suggesting facilitation. Both stochastic extinctions from standing trees and deterministic extinctions caused by tree fall should be accounted for in projecting epiphyte metapopulation dynamics. 4.: One to three infrequent, sexually dispersed study species face a significant extinction risk within 50 years, especially in areas with low tree numbers. 5.: Synthesis and applications. During the coming decades, infrequent, sexually dispersed, epiphytic lichens are likely to be lost from small woodland habitat set asides in intensively managed landscapes. Local extinction will be a consequence of low colonization rates and tree fall. Low colonization rates can be prevented by retaining large trees on which lichen species colonization rates are the highest and by assuring a high density of occupied trees. The negative effect of tree fall should be compensated for by assuring continuous availability of old trees. This can be achieved by decreasing the populations of large browsers, or by retaining trees with high conservation value during management operations.

Epiphyte metapopulation dynamics are explained by species traits, connectivity, and patch dynamics.

The colonization-extinction dynamics of many species are affected by the dynamics of their patches. For increasing our understanding of the metapopulation dynamics of sessile species confined to dynamic patches, we fitted a Bayesian incidence function model extended for dynamic landscapes to snapshot data on five epiphytic lichens among 2083 mapped oaks (dynamic patches). We estimate the age at which trees become suitable patches for different species, which defines their niche breadth (number of suitable trees). We show that the colonization rates were generally low, but increased with increasing connectivity in accordance with metapopulation theory. The rates were related to species traits, and we show, for the first time, that they are higher for species with wide niches and small dispersal propagules than for species with narrow niches or large propagules. We also show frequent long-distance dispersal in epiphytes by quantifying the relative importance of local dispersal and background deposition of dispersal propagules. Local stochastic extinctions from intact trees were negligible in all study species, and thus, the extinction rate is set by the rate of patch destruction (tree fall). These findings mean that epiphyte metapopulations may have slow colonization-extinction dynamics that are explained by connectivity, species traits, and patch dynamics.

Disentangling effects of uncertainties on population projections: climate change impact on an epixylic bryophyte.

Assessment of future ecosystem risks should account for the relevant uncertainty sources. This means accounting for the joint effects of climate variables and using modelling techniques that allow proper treatment of uncertainties. We investigate the influence of three of the IPCC's scenarios of greenhouse gas emissions (special report on emission scenarios (SRES)) on projections of the future abundance of a bryophyte model species. We also compare the relative importance of uncertainty sources on the population projections. The whole chain global climate model (GCM)-regional climate model-population dynamics model is addressed. The uncertainty depends on both natural- and model-related sources, in particular on GCM uncertainty. Ignoring the uncertainties gives an unwarranted impression of confidence in the results. The most likely population development of the bryophyte Buxbaumia viridis towards the end of this century is negative: even with a low-emission scenario, there is more than a 65 per cent risk for the population to be halved. The conclusion of a population decline is valid for all SRES scenarios investigated. Uncertainties are no longer an obstacle, but a mandatory aspect to include in the viability analysis of populations.