Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Using niche‐based modelling to assess the impact of climate change on tree functional diversity in Europe

Using niche‐based modelling to assess the impact of climate change on tree functional diversity... INTRODUCTION Rapid anthropogenic global climate change has already affected a range of species, communities and ecosystems ( Parmesan & Yohe, 2003 ). These impacts have taken the form of loss, expansion and relocation of habitats, changes in phenology and physiology ( Hughes, 2000 ; Walther ., 2002 ), coupled with possible repercussions on biotic interactions ( Hughes, 2000 ). Recent projections forecast that such a range of effects may be widespread in the future and affect both fauna and flora globally ( Peterson ., 2002 ; Thomas ., 2004 ; Thuiller ., 2005b ). Furthermore, there is abundant evidence from the Quaternary that plant species responses to past climate change have predominantly been expressed as large‐scale shifts in their geographical distributions ( Huntley, 1990 ; Webb, 1992 ). This expectation has been also supported by recent analyses showing that plants have greater degree of equilibrium with current climate than vertebrates with generally poorer dispersal abilities ( Araújo & Pearson, 2005 ). Nevertheless, the analogies form past to future species distributions must be made with caution because there is a possibility that non‐analogue situations might be created with current climate change. For example, levels of atmospheric CO 2 concentrations are greater now than in any period of the Pleistocene and these may cause changes in water use efficiency by plants, which may in turn influence competitive interactions ( Cowling & Sykes, 1999 ; Berry & Roderick, 2004 ). However, as rapid climate change ranks high among drivers of expected changes in biodiversity ( Sala ., 2000 ), it is useful to address this driver separately. Projections of the role of climate for future vegetation and habitat distributions are based on the idea of modelling stable current relationships between species distributions and environmental variables (predominantly climatic at large scales). These modelled relationships are then used to predict how species potential distributions may be altered in response to potential future climate change scenarios ( Sykes ., 1996 ; Cramer ., 2001 ; Bachelet ., 2003 ; Woodward & Lomas, 2004 ; Thuiller ., 2005a ). Recent studies have highlighted the potential ecosystem impacts of changes in functional diversity due to climate change ( Chapin III, 2003 ; Díaz & Cabido, 2001 ; Lavorel & Garnier, 2002 ). An important reason for concern about changes in species distributions and diversity is indeed that species can differ substantially in their effects on landscape structure and ecosystem processes ( Loreau ., 2001 ). The replacement of a species in a selected community by an ‘equivalent’ species ( Rosenfeld, 2002 ) belonging to the same functional group is likely to have a minimal impact on community and landscape structure or on ecosystem functioning as compared to the replacement by a species belonging to a different functional group ( Walker, 1992 ; Rosenfeld, 2002 ). For example, the replacement of a grass species by a woody deciduous broadleaved species would considerably modify carbon and nitrogen dynamics in a temperate ecosystem. Moreover organisms seldom affect a single process in isolation. Instead, effects often cascade through a broad range of ecosystem processes due to the tight linkages among carbon, nutrient and water cycles of ecosystems ( Tateno & Chapin III, 1997 ). Species also influence landscape structure and processes through their effects on the spread of materials, organisms or disturbance and on the response of landscapes to spatially heterogeneous perturbations. Large‐scale dynamic vegetation models have been developed in recent years to link dynamically changes in vegetation through time and space with the carbon and water cycles (e.g. Woodward ., 1998 ; Daly ., 2000 ; Sitch ., 2003 ). Given the number and variety of species at continental or global scales, appropriate levels of descriptions of plant function to simulate ecosystem dynamics and biogeochemistry have necessarily been coarse, using plant functional types rather than species ( Nemani & Running, 1996 ). Such plant functional types (PFT) do provide some information about plant diversity at a coarse level, but cannot be applied when information on individual plant species responses is required. To achieve this, there remains a fundamental need to establish bridges between these ecosystem models and detailed projections of changes in distributions for single species. Despite this, to date, few studies have attempted to extrapolate the likely impacts of climate change on functional diversity from single‐species projections. While various estimates of species‐extinction risk or species turnover are often taken as a surrogate for the sensitivity of species diversity to climate change, they do not allow modellers to project impacts on functional diversity with any degree of accuracy. Analyses that examine modifications to the functional diversity of ecosystems rather than focusing solely on extinction risks for single species would be needed to start assessing the functional impacts of species range shifts associated with climate change ( Peterson ., 2005 ). Here we focus on tree diversity within Europe (EU + Switzerland). Forests, woodlands and even individual trees are major elements in the European landscape, whether natural or managed. Forests occupy one‐third of the land area of Europe and have major functions for wood production, recreation, hunting and conservation ( Broekmeyer & Vos, 1993 ). Thus, trees are a major element in the functioning and the services provided by many European ecosystems. Tree species differ in the role they play in ecosystems. For example, broadleaved summer‐deciduous and needle‐leaved evergreen trees have very different effects on the structure and functioning of the whole ecosystem, including light and moisture regimes, or litter decomposition rates ( Broekmeyer & Vos, 1993 ). Hence, changing tree species ranges can affect ecosystem function but also the composition and diversity of understorey vegetation and associated biota ( Broekmeyer & Vos, 1993 ). In this paper we assess the potential impacts of a particular climate change scenario for 2080 on 112 tree species in Europe, and address the following questions: • What are the most susceptible functional types to climate change in Europe? • How will changes in species richness affect the tree functional diversity? • How will the spatial distribution of tree functional types in Europe be affected by climate warming as a result of changes in species distributions? METHODS Tree distributions and characteristics A total of 112 tree and tall shrubs species were selected for this study, encompassing all Gymnospermae families (Coniferales and Taxales), and a major part of the Angiospermae families (Salicales, Myricales, Juglandales, Fagales and Urticales). We did not include Rosales and Liliales as they were not available to us in a digital format. Distributional data for these taxa were obtained from the Atlas Florae Europaeae (AFE, Jalas & Suominen, 1972–96 ). Species data were originally located in 4419 UTM (Universal Transverse Mercator) 50 × 50 km grid cells ( Lahti & Lampinen, 1999 ), but we used only 2089 grid cells excluding most of the eastern European countries (except for the Baltic States) because of few records in these areas ( Williams ., 2000 ). Species were classified into broad functional types similar to those used by dynamic global vegetation models ( Sitch ., 2003 ). We retained four classes: deciduous broadleaved, evergreen broadleaved, evergreen coniferous and deciduous coniferous ( Tutin ., 1964–93; Bolòs ., 1990 ; Hodgson ., 1999 ; Peat & Fitter, 2002 ). Finally, because changes in functional diversity could be spatially structured, we also extracted the chorology (phytogeographical classes) of every species from Tutin (1964–93 ) and Bolòs . (1990, no. 1814). This allowed us to create subfunctional types based on the chorology of the species composing them (e.g. boreal evergreen coniferous vs. mediterranean evergreen coniferous). Bioclimatic data Seven bioclimatic variables interpolated for Europe in the context of the EU Framework 5 ATEAM project ( Schröter.2005 ) were used. All data were produced at a spatial resolution of 10′ for European grid cells based on the ATEAM geographical window and then aggregated to the Atlas Florae Europaeae 50 × 50 km grid. Bioclimate data were averaged for the period of 1961–90 and included mean annual temperature, mean temperature of the coldest month per year, mean annual precipitation sum, mean winter precipitation sum, mean summer precipitation sum, mean annual growing degree days and mean ratio of annual actual to potential evapotranspiration. Future projections for 2080 (averages for years 2070–99) were derived using one general circulation model experiment (HadCM3, Mitchell , 2004 ). This GCM projected mean climate under the A1Fi storyline ( IPCC, 2001 ), which describes a future world of very rapid economic growth, and with global population peaking by mid‐century and then declining, accompanied by the rapid introduction of new and more efficient technologies. This scenario is the one assuming greater intensity of energy use in the IPCC scenarios ( IPCC, 2001 ), and we chose it to address the extreme effects of climate change on functional diversity. For this exercise we ignored the impacts of changing land use on species range, although we acknowledge it can compound the effects of climate change on species ( Travis, 2003 ). However, at the spatial extent and resolution considered in this study, land use has been found to be largely overridden by climate ( Thuiller ., 2004a ) and should not affect the power of our analysis. Ecological niche modelling Models relating species distributions to the seven bioclimatic variables were fitted using the BIOMOD framework (see details in Thuiller, 2003, 2004 ) on a random sample of the initial data (70%). For each species, generalized linear models (GLM), generalized additive models (GAM), classification tree analysis (CTA) and artificial neural networks (ANN) were calibrated. Then each model for each species was evaluated on the remaining 30% of the initial data set using the values obtained either for Cohen's k ( Cohen, 1960 ) and for the area under the curve (AUC) of a receiver operating characteristic (ROC) plot of sensitivity against (1‐specificity) ( Swets, 1988 ). Sensitivity is defined as the proportion of true positives correctly predicted, whereas specificity is the proportion of true negatives correctly predicted ( Fielding & Bell, 1997 ). We transformed the probability of occurrence from models into presence–absence using either a threshold maximizing the percentage of presence and absence correctly predicted ( Pearce & Ferrier, 2000 ) or the probability threshold that maximized its performance as measured by Cohen's k ( Manel ., 2001 ). As we required only one projection per species for impact assessment into the future, we derived for each species eight projections (4 models × 2 threshold methods) and then, carried out an assessment of consensus using principal component analysis to explore central tendencies in the projections and select the niche‐based model representing the greatest commonality among projections ( Thuiller, 2004 ). There is an increasing evidence that model projections can be extremely variable and there remains a need to test the accuracy of models and to reduce uncertainties ( Thuiller, 2004 ; Thuiller ., 2004b ; Araújo ., 2005a ). One recent analysis has, however, provided the first test of the predictive accuracy of such models using observed range shifts of bird species and climate change in two periods of the recent past ( Araújo ., 2005b ). This work provides the first validation of niche‐based models under climate change and demonstrated how uncertainty can be reduced by selecting the most consensual projections, as done in the study presented here. We are therefore confident that this strategy provides a robust and defensible approach to the projections of species ranges for the purposes of conservation planning and biodiversity management. For each species, we estimated the number of pixels lost, gained (assuming unlimited migration), and unchanged, which allowed us to assess the degree of range change for each species. Current species richness (SR) and future SR assuming unlimited and no migration ability were estimated by summing the number of species present each pixel. Functional diversity estimation Concepts of functional diversity vary; according to Martinez (1997 ), functional diversity can be quantified by determining the nature and extent to which functional groups are represented in an ecological system. Functional diversity can also refer to the number of such groups in a community each of which contains one or more species ( Huston, 1994 ). Whatever the nature of the functional groups, it is generally accepted they will be fewer than the species under study ( Mooney, 1996 ). In this sense, functional ‘diversity’ is simply a measure of group richness rather than an estimate of evenness or dominance based on the abundance of individuals per group ( Petchey & Gaston, 2002 ). As with species diversity, it would seem reasonable to derive a parallel measure of functional diversity based on the abundance of individuals per functional type but without species‐weighting. While logically viable, this is likely to be limiting in practice as to record all individuals in a 50 × 50 km grid cell can be excessively time‐consuming at the European scale. Unsurprisingly, this information was not available to us. For these reasons, we explore the possibility of using species instead of individuals to serve as a ‘higher‐order ’ measure of abundance by deriving a species‐weighted, rather than a spatial or density‐driven, measure of functional diversity based on abundances of individuals. A species‐weighted form of functional diversity can therefore be defined as the diversity of functional types expressed as a function of the number of species per type. While the definition can be compared with that of Mooney (1996 ), this approach is more sensitive to evenness and dominance. We achieve this in the same way that species abundance is used to calculate species diversity but with the important difference that counts of species per functional type are used instead of counts of individuals per species. For this we apply the Simpson index of diversity as it takes into account both richness and evenness ( Magurran, 1988 ; Wilsey ., 2005 ): Index of Simpson: , where f i is the proportion of the PFT i (number of species) in a pixel. RESULTS The niche‐based models performed very well on average in the validation set (0.81 < AUC < 1, mean AUC = 0.96, 0.6 < k < 0.95, k = 0.74). Although not being an indisputable measure of performance ( Araújo ., 2005a ), good accuracy on the validation set increases our confidence that tree distributions might be governed by climatic factors at the European scale. Changes in functional richness The observed patterns of species range loss (assuming no dispersal ability) and range expansion (assuming unlimited migration) displayed distinct patterns according to the functional types and phytogeographical classes analysed ( Fig. 1 ). With regard to species range loss, two groups could be discerned. The first group, which encompasses boreal (both broadleaved deciduous and coniferous), Iberian (both broadleaved deciduous and evergreen), mountain (both broadleaved deciduous and deciduous and evergreen coniferous) and temperate (both broadleaved evergreen and coniferous) species, was predicted to lose on average > 60% of the current species ranges. The second group consisting of Mediterranean and temperate broadleaved deciduous species did not show consistent patterns. 1 Predicted percentage of tree species’ range change (a), expansion (b) and loss (c) by functional types and phytogeographical classes. Abbreviations: Bor, Boreal; Iber, Iberian (Splain and Portugal); Med, Mediterranean; Mount, Mountainous; Temp, temperate. Deci and Ever mean deciduous and evergreen, respectively. B and C mean broadleaved and coniferous, respectively. With regard to species expansions, two groups were identified. The first included only Mediterranean trees and showed greater range expansions (ranging from 40 to 140%) but with a large variation around the mean. The other group included all of the remaining species with small range expansions (< 50%) and low variability. Current and future (no and unlimited dispersal) species richness maps by functional type highlighted the northwards shift of deciduous communities ( Fig. 2 ). Assuming no dispersal, a remarkable decrease in the richness of deciduous species was shown in temperate regions and mountains. Evergreen species richness was also predicted to drop mainly to the south of Iberian Peninsula; however, this effect may not in fact be real as it could be due to the lack of information on the climatic preferences of the northern African ranges of these species ( Thuiller ., 2004c ). Finally, coniferous species showed contrasting patterns according to latitude. Mediterranean coniferous richness was predicted to strongly decrease in Central Spain, southern France and the Adriatic coast, supporting patterns revealed in Fig. 1 . Conversely, boreal coniferous richness was predicted to remain stable even in the absence of migration. Under the assumption of unlimited migration, there was a strong shift of deciduous richness from central Europe towards Scandinavia, compensated by a northward expansion of broadleaved evergreen trees. 2 Spatial patterns of current species richness (a), future species richness assuming null migration (b), and future species richness assuming unlimited migration (c) for each functional type. Respectively, broadleaved deciduous (1), broadleaved evergreen (2) and coniferous evergreen (3) trees. The grey scale legend corresponds to the number of species per class. Changes in functional diversity To assess the potential impact of climate change on functional diversity, we plotted the variation between future and current Simpson's diversity index ( Fig. 3 ). The severity of the change is much greater without migration. 3 Variability in tree functional diversity under null (a) and unlimited (b) migration assumptions. FD means functional diversity. Assuming no migration, the loss of deciduous species ( Figs 1 & 2 ) in mountain areas and on the Atlantic coast, in addition to the loss of richness among coniferous species in Scandinavia, made functional diversity higher in these two areas. This was realized by providing a greater evenness between deciduous and evergreen trees in southern Europe, and between deciduous and coniferous trees in northern Europe. Assuming unlimited migration, the expansion of evergreen species in central France compensated by a shift of deciduous species resulted in a net increase in functional diversity by providing an increase in the evenness of diversity in these areas currently dominated by deciduous broadleaved trees. Central western Europe was predicted to have a decrease of functional diversity because of a decrease in species richness of both deciduous and coniferous functional types. Relationships between variation in the index of diversity and variation in species richness, assuming unlimited dispersal, according to the biogeographical region showed that functional diversity was not systematically, nor always positively, related to species diversity ( Fig. 4 and Table 1 ). 4 Relationship between potential future change in tree species ‘richness’ change in tree functional diversity (change in Simpson's diversity index) for the different biogeographical regions analysed. Rho values represent Spearman rank correlations, which are all significant ( P < 0.05). 1 Spearman rank correlation between variation in index of diversity (δD) and current species (SR) richness, and between variation in index of diversity (δD) and variation in species richness (δSR) δD vs. SR δD vs. δSR No migration Unlimited migration No migration Unlimited migration Alpine north −0.01 (NS) 0.27 −0.55 0.49 Boreal 0.16 0.45 −0.44 0.36 Nemoral −0.23 −0.14 0.59 0.57 Atlantic −0.50 −0.66 0.37 0.39 Alpine south 0.22 0.19 0.32 0.28 Continental −0.42 −0.29 0.32 0.05 Atlantic central 0.29 0.15 −0.26 −0.58 Pannnonian −0.15 −0.33 0.30 0.26 Lusitanian 0.13 0.27 −0.29 −0.40 Med. mountain −0.14 −0.21 0.13 0.13 Med. north 0.33 0.47 −0.15 −0.40 Med. south 0.31 0.35 −0.27 −0.45 The Mediterranean, Central Atlantic coast and Lusitanian regions showed a negative correlation (increase in functional diversity but decrease in species richness). In the case of the Mediterranean and Lusitanian regions, the trend reflected a decrease in both broadleaved deciduous and Mediterranean coniferous. Their decrease was likely due to increasing drought, which should lead to more open vegetation and thereby more shrubs ( Gritti ., 2005 ). The decrease was strongest for broadleaved deciduous, therefore, the evenness of the tree communities decreases, leading to decreased functional diversity. For the central Atlantic region the main pattern was a loss in the strong dominance of broadleaved deciduous trees along with colonization by Mediterranean broadleaved evergreen trees. DISCUSSION Advantages and caveats of the approach We will not discuss the pros and cons of the ecological niche modelling approach as several papers have already highlighted their limitations ( Woodward & Beerling, 1997 ; Davis ., 1998 ; Hampe, 2004 ; Svenning & Skov, 2004 ; Guisan & Thuiller, 2005 ) and the underlying methodological uncertainties ( Thuiller, 2004 ; Araújo ., 2005b ). Nevertheless, one of the major weaknesses of current niche‐based models is that they cannot take into account the migration process of species as they track environmental change. Although palaeological records and simulations for the Quaternary have suggested a range of migration rates up to 1000 m year −1 ( Cain ., 1998 ; Clark, 1998 ), it is difficult to estimate the migration process for each simulated species ( Higgins ., 2003 ). Furthermore, without simulations of both the current highly fragmented landscape within Europe and the possible future changes in land use, time lags are unknown. In this paper, therefore, we contrasted two simple assumptions, no dispersal beyond a species present site, as opposed to unlimited dispersal to all sites ( Peterson ., 2002 ; Thuiller ., 2005b ). Reality for species lies within these boundaries. In this study, we deliberately focused only on functional and structural aspects of impacts of climate change using broadly defined tree functional types. Functional diversity, the extent of functional differences among the species in a community ( Tilman, 2001 ), is an important determinant of ecosystem processes. There is an ongoing debate on how to estimate functional diversity. A common measure of functional diversity is the number of functional groups represented by the species in a community ( Hooper, 1998 ; Tilman, 2001 ). Of the problems associated with assigning species to groups ( Lavorel ., 1997 ), perhaps the least tractable is the arbitrary scale at which differences between species qualify as functionally significant ( Simberloff & Dayan, 1991 ). There is no simple, satisfactory or standardized measure of functional diversity ( Díaz & Cabido, 2001 ). Recently, some new measures of functional diversity have been proposed based on functional distances among taxa ( Petchey & Gaston, 2002 ), similar to methods used to estimate phylogenetic diversity ( Faith, 2004 ; but see Clarke & Warwick, 1999 ). Further analyses would thus be required for using such estimates. Instead, we used Simpson's diversity index, which combines a measure of the number of functional types in one pixel and a measure of evenness ( Magurran, 1988 ). Impact of climate change on functional diversity Climate is one major driver of species, biomes and functional type distributions ( Woodward, 1987 ). In this study, we apply this premise to assess the potential impact of climate change on tree functional diversity. The major objective of our analyses was then to assess how species range shifts and variation in species richness could produce changes in functional assemblages. Our results clearly show that is a likely outcome. Based on phytogeographical classes, we showed that boreal broadleaved deciduous and coniferous species are likely to undergo large range reductions (no migration hypothesis) and remain unlikely to find new suitable habitats (unlimited migration hypothesis). Potential replacements may be by temperate broadleaved deciduous and temperate coniferous species. Such changes are unlikely to produce strong functional changes if they largely represent a replacement by species belonging to different phytogeographical class but within the same functional types that are already present. We, however, predicted an increase of functional diversity in these areas, resulting from an increased evenness of the tree communities. These changes in functional richness and diversity could then have significant structural and functional effects on ecosystems ( Cramer ., 1999 ). For example, the timing and frequency of phenological events such as leafing out or leaf fall influence the distribution of light both in the forest canopy and on the forest floor, thereby directly affecting ecosystem functioning ( Martens ., 2000 ). Soil carbon content and thus soil biogeochemistry depend on litter type and different timings for litter fall ( Gill & Burke, 1999 ). The observed increase in the share of deciduous species could also strongly modify ecosystems and especially their microclimate and surface properties ( Chapin III ., 2000 ). For instance, the albedo of the vegetation soil surface determines the proportion of incoming radiation that is absorbed, which in turn influences the quantity of heat that reaches the ground and the energy available to drive water loss from the ecosystem through transpiration and soil evapotranspiration. Forests dominated by conifers have a lower albedo (i.e. reflecting less energy back to the atmosphere) than do deciduous‐dominated forests because the low albedo of coniferous canopies and clumped distribution of leaves cause much of the reflected radiation in these forests to be absorbed by other leaves rather than being reflected directly back to space ( Chapin III, 2003 ). Finally, fire regimes, forest productivity and thus carbon storage are often species‐dependent and may be modified by changes in species and/or functional composition ( Jonhnstone & Chapin III, 2003 ). Changes in the ratio between coniferous and deciduous trees could influence the probability or severity of fires. Boreal conifers are more flammable than deciduous trees because of their large leaf and twig surface area, low moisture content and high resin content ( Johnson, 1992 ). Similar effects have been documented during the invasion of boreal forests in Alaska by black spruce in the middle Holocene, which caused an increase in fire frequency ( Lynch ., 2003 ). Conversely, dynamic models predict an increase in deciduous species and concurrent changes in fire regimes for these regions under scenarios of climate change ( Rupp ., 2000 ). In the southern part of temperate Europe, we predicted a replacement of temperate deciduous broadleaved species by Mediterranean trees (all functional types). These resulted in an increase in functional diversity, producing communities more functionally mixed. Greater phenological diversity may lead to a longer period of annual functioning and therefore greater stocks in carbon. The ranges of temperate tree species can often be related to the length of the growing season (GDD5) and especially to absolute minimum temperature ( Woodward, 1987 ). In particular, the limit where broadleaved evergreen species are replaced by deciduous species is mostly climatically determined ( Woodward, 1987 ). However, the northern limit of Mediterranean species in Europe appears to be controlled by a combination of growing season length (or cumulated heat) and minimum winter temperature ( Woodward, 1987 ). Climate change is likely to be reflected in these with increasing growing season length and decreasing minimum temperatures. Such changes are likely to lead to the northwards expansion of Mediterranean species, principally evergreen, but also deciduous, causing range reductions at least at the southern edge of temperate deciduous and coniferous species distributions. Our simulated patterns are confirmed by observations in Europe ( Walther ., 2001 ). For example, Ilex aquifolium (holly) has been used as a classical example of a climate‐limited distribution ( Larcher, 2003 ). A synchronized shift in both its distribution and climate was observed in southern Scandinavia ( Walther ., 2005 ), confirming the relevance of climatic envelopes to model and predict the distribution of temperate trees. Walther . (2001 ) also noted the capacity of other climate‐sensitive transition zones to respond to changes in environmental conditions and observed that in Switzerland, Japan and Georgia, indigenous and exotic evergreen broadleaved species are already expanding their ranges into deciduous forests habitats. This replacement of deciduous by evergreen broadleaved trees could have significant functional impacts such as a reduction of cycling for biogeochemistry and especially of decomposition and nitrogen mineralization. Litter from species typical of productive environments (e.g. deciduous trees) typically decomposes more rapidly than that species from less productive environments (e.g. evergreen trees) ( Cornelissen, 1996 ). The possible increase of shrubs in Mediterranean and Lusitanian regions as a result of the regression by broadleaved deciduous trees could increase the evenness and have profound impacts on the landscape pattern and the disturbance regimes. For instance, shrub‐dominated landscapes are more prone to invasion by grasses or herbs than forests, which could produce strong fire regime changes ( Grigulis ., 2005 ). Unfortunately, the static nature of our modelling, as well as the scale and resolution considered, does not allow us to derive more detailed conclusions about the effect of climate change on ecosystem processes. Even if the link between plant traits and ecosystems is now well recognized ( Chapin III, 2003 ), the lack of information about abundance of each species in a pixel precludes the projection of community structure, and thereby the possibility to link trait values with ecosystem function ( Garnier ., 2004 ). ACKNOWLEDGEMENTS This research was funded by the European Commission's FP5 ATEAM (Advanced Terrestrial Ecosystem Analysis and Modelling, EVK2‐CT‐ 2000–00075) and FP6 ALARM (Assessing LArge‐scale environmental Risks with tested Methods, GOCE‐CT‐2003‐506675) integrated projects. MBA thanks T. Lahti and R. Lampinen for providing a digital version of Atlas Florae Europaeae . We thank T.D. Mitchell of the Tyndall Centre for Climate Change Research for providing climate and scenario data, and S. Zaehle and M. Erhard at the Potsdam Institute for Climate Impact Research for aggregating climate data. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Diversity and Distributions Wiley

Using niche‐based modelling to assess the impact of climate change on tree functional diversity in Europe

Loading next page...
 
/lp/wiley/using-niche-based-modelling-to-assess-the-impact-of-climate-change-on-Qmux3gTTs7

References (96)

Publisher
Wiley
Copyright
Copyright © 2006 Wiley Subscription Services, Inc., A Wiley Company
ISSN
1366-9516
eISSN
1472-4642
DOI
10.1111/j.1366-9516.2006.00216.x
Publisher site
See Article on Publisher Site

Abstract

INTRODUCTION Rapid anthropogenic global climate change has already affected a range of species, communities and ecosystems ( Parmesan & Yohe, 2003 ). These impacts have taken the form of loss, expansion and relocation of habitats, changes in phenology and physiology ( Hughes, 2000 ; Walther ., 2002 ), coupled with possible repercussions on biotic interactions ( Hughes, 2000 ). Recent projections forecast that such a range of effects may be widespread in the future and affect both fauna and flora globally ( Peterson ., 2002 ; Thomas ., 2004 ; Thuiller ., 2005b ). Furthermore, there is abundant evidence from the Quaternary that plant species responses to past climate change have predominantly been expressed as large‐scale shifts in their geographical distributions ( Huntley, 1990 ; Webb, 1992 ). This expectation has been also supported by recent analyses showing that plants have greater degree of equilibrium with current climate than vertebrates with generally poorer dispersal abilities ( Araújo & Pearson, 2005 ). Nevertheless, the analogies form past to future species distributions must be made with caution because there is a possibility that non‐analogue situations might be created with current climate change. For example, levels of atmospheric CO 2 concentrations are greater now than in any period of the Pleistocene and these may cause changes in water use efficiency by plants, which may in turn influence competitive interactions ( Cowling & Sykes, 1999 ; Berry & Roderick, 2004 ). However, as rapid climate change ranks high among drivers of expected changes in biodiversity ( Sala ., 2000 ), it is useful to address this driver separately. Projections of the role of climate for future vegetation and habitat distributions are based on the idea of modelling stable current relationships between species distributions and environmental variables (predominantly climatic at large scales). These modelled relationships are then used to predict how species potential distributions may be altered in response to potential future climate change scenarios ( Sykes ., 1996 ; Cramer ., 2001 ; Bachelet ., 2003 ; Woodward & Lomas, 2004 ; Thuiller ., 2005a ). Recent studies have highlighted the potential ecosystem impacts of changes in functional diversity due to climate change ( Chapin III, 2003 ; Díaz & Cabido, 2001 ; Lavorel & Garnier, 2002 ). An important reason for concern about changes in species distributions and diversity is indeed that species can differ substantially in their effects on landscape structure and ecosystem processes ( Loreau ., 2001 ). The replacement of a species in a selected community by an ‘equivalent’ species ( Rosenfeld, 2002 ) belonging to the same functional group is likely to have a minimal impact on community and landscape structure or on ecosystem functioning as compared to the replacement by a species belonging to a different functional group ( Walker, 1992 ; Rosenfeld, 2002 ). For example, the replacement of a grass species by a woody deciduous broadleaved species would considerably modify carbon and nitrogen dynamics in a temperate ecosystem. Moreover organisms seldom affect a single process in isolation. Instead, effects often cascade through a broad range of ecosystem processes due to the tight linkages among carbon, nutrient and water cycles of ecosystems ( Tateno & Chapin III, 1997 ). Species also influence landscape structure and processes through their effects on the spread of materials, organisms or disturbance and on the response of landscapes to spatially heterogeneous perturbations. Large‐scale dynamic vegetation models have been developed in recent years to link dynamically changes in vegetation through time and space with the carbon and water cycles (e.g. Woodward ., 1998 ; Daly ., 2000 ; Sitch ., 2003 ). Given the number and variety of species at continental or global scales, appropriate levels of descriptions of plant function to simulate ecosystem dynamics and biogeochemistry have necessarily been coarse, using plant functional types rather than species ( Nemani & Running, 1996 ). Such plant functional types (PFT) do provide some information about plant diversity at a coarse level, but cannot be applied when information on individual plant species responses is required. To achieve this, there remains a fundamental need to establish bridges between these ecosystem models and detailed projections of changes in distributions for single species. Despite this, to date, few studies have attempted to extrapolate the likely impacts of climate change on functional diversity from single‐species projections. While various estimates of species‐extinction risk or species turnover are often taken as a surrogate for the sensitivity of species diversity to climate change, they do not allow modellers to project impacts on functional diversity with any degree of accuracy. Analyses that examine modifications to the functional diversity of ecosystems rather than focusing solely on extinction risks for single species would be needed to start assessing the functional impacts of species range shifts associated with climate change ( Peterson ., 2005 ). Here we focus on tree diversity within Europe (EU + Switzerland). Forests, woodlands and even individual trees are major elements in the European landscape, whether natural or managed. Forests occupy one‐third of the land area of Europe and have major functions for wood production, recreation, hunting and conservation ( Broekmeyer & Vos, 1993 ). Thus, trees are a major element in the functioning and the services provided by many European ecosystems. Tree species differ in the role they play in ecosystems. For example, broadleaved summer‐deciduous and needle‐leaved evergreen trees have very different effects on the structure and functioning of the whole ecosystem, including light and moisture regimes, or litter decomposition rates ( Broekmeyer & Vos, 1993 ). Hence, changing tree species ranges can affect ecosystem function but also the composition and diversity of understorey vegetation and associated biota ( Broekmeyer & Vos, 1993 ). In this paper we assess the potential impacts of a particular climate change scenario for 2080 on 112 tree species in Europe, and address the following questions: • What are the most susceptible functional types to climate change in Europe? • How will changes in species richness affect the tree functional diversity? • How will the spatial distribution of tree functional types in Europe be affected by climate warming as a result of changes in species distributions? METHODS Tree distributions and characteristics A total of 112 tree and tall shrubs species were selected for this study, encompassing all Gymnospermae families (Coniferales and Taxales), and a major part of the Angiospermae families (Salicales, Myricales, Juglandales, Fagales and Urticales). We did not include Rosales and Liliales as they were not available to us in a digital format. Distributional data for these taxa were obtained from the Atlas Florae Europaeae (AFE, Jalas & Suominen, 1972–96 ). Species data were originally located in 4419 UTM (Universal Transverse Mercator) 50 × 50 km grid cells ( Lahti & Lampinen, 1999 ), but we used only 2089 grid cells excluding most of the eastern European countries (except for the Baltic States) because of few records in these areas ( Williams ., 2000 ). Species were classified into broad functional types similar to those used by dynamic global vegetation models ( Sitch ., 2003 ). We retained four classes: deciduous broadleaved, evergreen broadleaved, evergreen coniferous and deciduous coniferous ( Tutin ., 1964–93; Bolòs ., 1990 ; Hodgson ., 1999 ; Peat & Fitter, 2002 ). Finally, because changes in functional diversity could be spatially structured, we also extracted the chorology (phytogeographical classes) of every species from Tutin (1964–93 ) and Bolòs . (1990, no. 1814). This allowed us to create subfunctional types based on the chorology of the species composing them (e.g. boreal evergreen coniferous vs. mediterranean evergreen coniferous). Bioclimatic data Seven bioclimatic variables interpolated for Europe in the context of the EU Framework 5 ATEAM project ( Schröter.2005 ) were used. All data were produced at a spatial resolution of 10′ for European grid cells based on the ATEAM geographical window and then aggregated to the Atlas Florae Europaeae 50 × 50 km grid. Bioclimate data were averaged for the period of 1961–90 and included mean annual temperature, mean temperature of the coldest month per year, mean annual precipitation sum, mean winter precipitation sum, mean summer precipitation sum, mean annual growing degree days and mean ratio of annual actual to potential evapotranspiration. Future projections for 2080 (averages for years 2070–99) were derived using one general circulation model experiment (HadCM3, Mitchell , 2004 ). This GCM projected mean climate under the A1Fi storyline ( IPCC, 2001 ), which describes a future world of very rapid economic growth, and with global population peaking by mid‐century and then declining, accompanied by the rapid introduction of new and more efficient technologies. This scenario is the one assuming greater intensity of energy use in the IPCC scenarios ( IPCC, 2001 ), and we chose it to address the extreme effects of climate change on functional diversity. For this exercise we ignored the impacts of changing land use on species range, although we acknowledge it can compound the effects of climate change on species ( Travis, 2003 ). However, at the spatial extent and resolution considered in this study, land use has been found to be largely overridden by climate ( Thuiller ., 2004a ) and should not affect the power of our analysis. Ecological niche modelling Models relating species distributions to the seven bioclimatic variables were fitted using the BIOMOD framework (see details in Thuiller, 2003, 2004 ) on a random sample of the initial data (70%). For each species, generalized linear models (GLM), generalized additive models (GAM), classification tree analysis (CTA) and artificial neural networks (ANN) were calibrated. Then each model for each species was evaluated on the remaining 30% of the initial data set using the values obtained either for Cohen's k ( Cohen, 1960 ) and for the area under the curve (AUC) of a receiver operating characteristic (ROC) plot of sensitivity against (1‐specificity) ( Swets, 1988 ). Sensitivity is defined as the proportion of true positives correctly predicted, whereas specificity is the proportion of true negatives correctly predicted ( Fielding & Bell, 1997 ). We transformed the probability of occurrence from models into presence–absence using either a threshold maximizing the percentage of presence and absence correctly predicted ( Pearce & Ferrier, 2000 ) or the probability threshold that maximized its performance as measured by Cohen's k ( Manel ., 2001 ). As we required only one projection per species for impact assessment into the future, we derived for each species eight projections (4 models × 2 threshold methods) and then, carried out an assessment of consensus using principal component analysis to explore central tendencies in the projections and select the niche‐based model representing the greatest commonality among projections ( Thuiller, 2004 ). There is an increasing evidence that model projections can be extremely variable and there remains a need to test the accuracy of models and to reduce uncertainties ( Thuiller, 2004 ; Thuiller ., 2004b ; Araújo ., 2005a ). One recent analysis has, however, provided the first test of the predictive accuracy of such models using observed range shifts of bird species and climate change in two periods of the recent past ( Araújo ., 2005b ). This work provides the first validation of niche‐based models under climate change and demonstrated how uncertainty can be reduced by selecting the most consensual projections, as done in the study presented here. We are therefore confident that this strategy provides a robust and defensible approach to the projections of species ranges for the purposes of conservation planning and biodiversity management. For each species, we estimated the number of pixels lost, gained (assuming unlimited migration), and unchanged, which allowed us to assess the degree of range change for each species. Current species richness (SR) and future SR assuming unlimited and no migration ability were estimated by summing the number of species present each pixel. Functional diversity estimation Concepts of functional diversity vary; according to Martinez (1997 ), functional diversity can be quantified by determining the nature and extent to which functional groups are represented in an ecological system. Functional diversity can also refer to the number of such groups in a community each of which contains one or more species ( Huston, 1994 ). Whatever the nature of the functional groups, it is generally accepted they will be fewer than the species under study ( Mooney, 1996 ). In this sense, functional ‘diversity’ is simply a measure of group richness rather than an estimate of evenness or dominance based on the abundance of individuals per group ( Petchey & Gaston, 2002 ). As with species diversity, it would seem reasonable to derive a parallel measure of functional diversity based on the abundance of individuals per functional type but without species‐weighting. While logically viable, this is likely to be limiting in practice as to record all individuals in a 50 × 50 km grid cell can be excessively time‐consuming at the European scale. Unsurprisingly, this information was not available to us. For these reasons, we explore the possibility of using species instead of individuals to serve as a ‘higher‐order ’ measure of abundance by deriving a species‐weighted, rather than a spatial or density‐driven, measure of functional diversity based on abundances of individuals. A species‐weighted form of functional diversity can therefore be defined as the diversity of functional types expressed as a function of the number of species per type. While the definition can be compared with that of Mooney (1996 ), this approach is more sensitive to evenness and dominance. We achieve this in the same way that species abundance is used to calculate species diversity but with the important difference that counts of species per functional type are used instead of counts of individuals per species. For this we apply the Simpson index of diversity as it takes into account both richness and evenness ( Magurran, 1988 ; Wilsey ., 2005 ): Index of Simpson: , where f i is the proportion of the PFT i (number of species) in a pixel. RESULTS The niche‐based models performed very well on average in the validation set (0.81 < AUC < 1, mean AUC = 0.96, 0.6 < k < 0.95, k = 0.74). Although not being an indisputable measure of performance ( Araújo ., 2005a ), good accuracy on the validation set increases our confidence that tree distributions might be governed by climatic factors at the European scale. Changes in functional richness The observed patterns of species range loss (assuming no dispersal ability) and range expansion (assuming unlimited migration) displayed distinct patterns according to the functional types and phytogeographical classes analysed ( Fig. 1 ). With regard to species range loss, two groups could be discerned. The first group, which encompasses boreal (both broadleaved deciduous and coniferous), Iberian (both broadleaved deciduous and evergreen), mountain (both broadleaved deciduous and deciduous and evergreen coniferous) and temperate (both broadleaved evergreen and coniferous) species, was predicted to lose on average > 60% of the current species ranges. The second group consisting of Mediterranean and temperate broadleaved deciduous species did not show consistent patterns. 1 Predicted percentage of tree species’ range change (a), expansion (b) and loss (c) by functional types and phytogeographical classes. Abbreviations: Bor, Boreal; Iber, Iberian (Splain and Portugal); Med, Mediterranean; Mount, Mountainous; Temp, temperate. Deci and Ever mean deciduous and evergreen, respectively. B and C mean broadleaved and coniferous, respectively. With regard to species expansions, two groups were identified. The first included only Mediterranean trees and showed greater range expansions (ranging from 40 to 140%) but with a large variation around the mean. The other group included all of the remaining species with small range expansions (< 50%) and low variability. Current and future (no and unlimited dispersal) species richness maps by functional type highlighted the northwards shift of deciduous communities ( Fig. 2 ). Assuming no dispersal, a remarkable decrease in the richness of deciduous species was shown in temperate regions and mountains. Evergreen species richness was also predicted to drop mainly to the south of Iberian Peninsula; however, this effect may not in fact be real as it could be due to the lack of information on the climatic preferences of the northern African ranges of these species ( Thuiller ., 2004c ). Finally, coniferous species showed contrasting patterns according to latitude. Mediterranean coniferous richness was predicted to strongly decrease in Central Spain, southern France and the Adriatic coast, supporting patterns revealed in Fig. 1 . Conversely, boreal coniferous richness was predicted to remain stable even in the absence of migration. Under the assumption of unlimited migration, there was a strong shift of deciduous richness from central Europe towards Scandinavia, compensated by a northward expansion of broadleaved evergreen trees. 2 Spatial patterns of current species richness (a), future species richness assuming null migration (b), and future species richness assuming unlimited migration (c) for each functional type. Respectively, broadleaved deciduous (1), broadleaved evergreen (2) and coniferous evergreen (3) trees. The grey scale legend corresponds to the number of species per class. Changes in functional diversity To assess the potential impact of climate change on functional diversity, we plotted the variation between future and current Simpson's diversity index ( Fig. 3 ). The severity of the change is much greater without migration. 3 Variability in tree functional diversity under null (a) and unlimited (b) migration assumptions. FD means functional diversity. Assuming no migration, the loss of deciduous species ( Figs 1 & 2 ) in mountain areas and on the Atlantic coast, in addition to the loss of richness among coniferous species in Scandinavia, made functional diversity higher in these two areas. This was realized by providing a greater evenness between deciduous and evergreen trees in southern Europe, and between deciduous and coniferous trees in northern Europe. Assuming unlimited migration, the expansion of evergreen species in central France compensated by a shift of deciduous species resulted in a net increase in functional diversity by providing an increase in the evenness of diversity in these areas currently dominated by deciduous broadleaved trees. Central western Europe was predicted to have a decrease of functional diversity because of a decrease in species richness of both deciduous and coniferous functional types. Relationships between variation in the index of diversity and variation in species richness, assuming unlimited dispersal, according to the biogeographical region showed that functional diversity was not systematically, nor always positively, related to species diversity ( Fig. 4 and Table 1 ). 4 Relationship between potential future change in tree species ‘richness’ change in tree functional diversity (change in Simpson's diversity index) for the different biogeographical regions analysed. Rho values represent Spearman rank correlations, which are all significant ( P < 0.05). 1 Spearman rank correlation between variation in index of diversity (δD) and current species (SR) richness, and between variation in index of diversity (δD) and variation in species richness (δSR) δD vs. SR δD vs. δSR No migration Unlimited migration No migration Unlimited migration Alpine north −0.01 (NS) 0.27 −0.55 0.49 Boreal 0.16 0.45 −0.44 0.36 Nemoral −0.23 −0.14 0.59 0.57 Atlantic −0.50 −0.66 0.37 0.39 Alpine south 0.22 0.19 0.32 0.28 Continental −0.42 −0.29 0.32 0.05 Atlantic central 0.29 0.15 −0.26 −0.58 Pannnonian −0.15 −0.33 0.30 0.26 Lusitanian 0.13 0.27 −0.29 −0.40 Med. mountain −0.14 −0.21 0.13 0.13 Med. north 0.33 0.47 −0.15 −0.40 Med. south 0.31 0.35 −0.27 −0.45 The Mediterranean, Central Atlantic coast and Lusitanian regions showed a negative correlation (increase in functional diversity but decrease in species richness). In the case of the Mediterranean and Lusitanian regions, the trend reflected a decrease in both broadleaved deciduous and Mediterranean coniferous. Their decrease was likely due to increasing drought, which should lead to more open vegetation and thereby more shrubs ( Gritti ., 2005 ). The decrease was strongest for broadleaved deciduous, therefore, the evenness of the tree communities decreases, leading to decreased functional diversity. For the central Atlantic region the main pattern was a loss in the strong dominance of broadleaved deciduous trees along with colonization by Mediterranean broadleaved evergreen trees. DISCUSSION Advantages and caveats of the approach We will not discuss the pros and cons of the ecological niche modelling approach as several papers have already highlighted their limitations ( Woodward & Beerling, 1997 ; Davis ., 1998 ; Hampe, 2004 ; Svenning & Skov, 2004 ; Guisan & Thuiller, 2005 ) and the underlying methodological uncertainties ( Thuiller, 2004 ; Araújo ., 2005b ). Nevertheless, one of the major weaknesses of current niche‐based models is that they cannot take into account the migration process of species as they track environmental change. Although palaeological records and simulations for the Quaternary have suggested a range of migration rates up to 1000 m year −1 ( Cain ., 1998 ; Clark, 1998 ), it is difficult to estimate the migration process for each simulated species ( Higgins ., 2003 ). Furthermore, without simulations of both the current highly fragmented landscape within Europe and the possible future changes in land use, time lags are unknown. In this paper, therefore, we contrasted two simple assumptions, no dispersal beyond a species present site, as opposed to unlimited dispersal to all sites ( Peterson ., 2002 ; Thuiller ., 2005b ). Reality for species lies within these boundaries. In this study, we deliberately focused only on functional and structural aspects of impacts of climate change using broadly defined tree functional types. Functional diversity, the extent of functional differences among the species in a community ( Tilman, 2001 ), is an important determinant of ecosystem processes. There is an ongoing debate on how to estimate functional diversity. A common measure of functional diversity is the number of functional groups represented by the species in a community ( Hooper, 1998 ; Tilman, 2001 ). Of the problems associated with assigning species to groups ( Lavorel ., 1997 ), perhaps the least tractable is the arbitrary scale at which differences between species qualify as functionally significant ( Simberloff & Dayan, 1991 ). There is no simple, satisfactory or standardized measure of functional diversity ( Díaz & Cabido, 2001 ). Recently, some new measures of functional diversity have been proposed based on functional distances among taxa ( Petchey & Gaston, 2002 ), similar to methods used to estimate phylogenetic diversity ( Faith, 2004 ; but see Clarke & Warwick, 1999 ). Further analyses would thus be required for using such estimates. Instead, we used Simpson's diversity index, which combines a measure of the number of functional types in one pixel and a measure of evenness ( Magurran, 1988 ). Impact of climate change on functional diversity Climate is one major driver of species, biomes and functional type distributions ( Woodward, 1987 ). In this study, we apply this premise to assess the potential impact of climate change on tree functional diversity. The major objective of our analyses was then to assess how species range shifts and variation in species richness could produce changes in functional assemblages. Our results clearly show that is a likely outcome. Based on phytogeographical classes, we showed that boreal broadleaved deciduous and coniferous species are likely to undergo large range reductions (no migration hypothesis) and remain unlikely to find new suitable habitats (unlimited migration hypothesis). Potential replacements may be by temperate broadleaved deciduous and temperate coniferous species. Such changes are unlikely to produce strong functional changes if they largely represent a replacement by species belonging to different phytogeographical class but within the same functional types that are already present. We, however, predicted an increase of functional diversity in these areas, resulting from an increased evenness of the tree communities. These changes in functional richness and diversity could then have significant structural and functional effects on ecosystems ( Cramer ., 1999 ). For example, the timing and frequency of phenological events such as leafing out or leaf fall influence the distribution of light both in the forest canopy and on the forest floor, thereby directly affecting ecosystem functioning ( Martens ., 2000 ). Soil carbon content and thus soil biogeochemistry depend on litter type and different timings for litter fall ( Gill & Burke, 1999 ). The observed increase in the share of deciduous species could also strongly modify ecosystems and especially their microclimate and surface properties ( Chapin III ., 2000 ). For instance, the albedo of the vegetation soil surface determines the proportion of incoming radiation that is absorbed, which in turn influences the quantity of heat that reaches the ground and the energy available to drive water loss from the ecosystem through transpiration and soil evapotranspiration. Forests dominated by conifers have a lower albedo (i.e. reflecting less energy back to the atmosphere) than do deciduous‐dominated forests because the low albedo of coniferous canopies and clumped distribution of leaves cause much of the reflected radiation in these forests to be absorbed by other leaves rather than being reflected directly back to space ( Chapin III, 2003 ). Finally, fire regimes, forest productivity and thus carbon storage are often species‐dependent and may be modified by changes in species and/or functional composition ( Jonhnstone & Chapin III, 2003 ). Changes in the ratio between coniferous and deciduous trees could influence the probability or severity of fires. Boreal conifers are more flammable than deciduous trees because of their large leaf and twig surface area, low moisture content and high resin content ( Johnson, 1992 ). Similar effects have been documented during the invasion of boreal forests in Alaska by black spruce in the middle Holocene, which caused an increase in fire frequency ( Lynch ., 2003 ). Conversely, dynamic models predict an increase in deciduous species and concurrent changes in fire regimes for these regions under scenarios of climate change ( Rupp ., 2000 ). In the southern part of temperate Europe, we predicted a replacement of temperate deciduous broadleaved species by Mediterranean trees (all functional types). These resulted in an increase in functional diversity, producing communities more functionally mixed. Greater phenological diversity may lead to a longer period of annual functioning and therefore greater stocks in carbon. The ranges of temperate tree species can often be related to the length of the growing season (GDD5) and especially to absolute minimum temperature ( Woodward, 1987 ). In particular, the limit where broadleaved evergreen species are replaced by deciduous species is mostly climatically determined ( Woodward, 1987 ). However, the northern limit of Mediterranean species in Europe appears to be controlled by a combination of growing season length (or cumulated heat) and minimum winter temperature ( Woodward, 1987 ). Climate change is likely to be reflected in these with increasing growing season length and decreasing minimum temperatures. Such changes are likely to lead to the northwards expansion of Mediterranean species, principally evergreen, but also deciduous, causing range reductions at least at the southern edge of temperate deciduous and coniferous species distributions. Our simulated patterns are confirmed by observations in Europe ( Walther ., 2001 ). For example, Ilex aquifolium (holly) has been used as a classical example of a climate‐limited distribution ( Larcher, 2003 ). A synchronized shift in both its distribution and climate was observed in southern Scandinavia ( Walther ., 2005 ), confirming the relevance of climatic envelopes to model and predict the distribution of temperate trees. Walther . (2001 ) also noted the capacity of other climate‐sensitive transition zones to respond to changes in environmental conditions and observed that in Switzerland, Japan and Georgia, indigenous and exotic evergreen broadleaved species are already expanding their ranges into deciduous forests habitats. This replacement of deciduous by evergreen broadleaved trees could have significant functional impacts such as a reduction of cycling for biogeochemistry and especially of decomposition and nitrogen mineralization. Litter from species typical of productive environments (e.g. deciduous trees) typically decomposes more rapidly than that species from less productive environments (e.g. evergreen trees) ( Cornelissen, 1996 ). The possible increase of shrubs in Mediterranean and Lusitanian regions as a result of the regression by broadleaved deciduous trees could increase the evenness and have profound impacts on the landscape pattern and the disturbance regimes. For instance, shrub‐dominated landscapes are more prone to invasion by grasses or herbs than forests, which could produce strong fire regime changes ( Grigulis ., 2005 ). Unfortunately, the static nature of our modelling, as well as the scale and resolution considered, does not allow us to derive more detailed conclusions about the effect of climate change on ecosystem processes. Even if the link between plant traits and ecosystems is now well recognized ( Chapin III, 2003 ), the lack of information about abundance of each species in a pixel precludes the projection of community structure, and thereby the possibility to link trait values with ecosystem function ( Garnier ., 2004 ). ACKNOWLEDGEMENTS This research was funded by the European Commission's FP5 ATEAM (Advanced Terrestrial Ecosystem Analysis and Modelling, EVK2‐CT‐ 2000–00075) and FP6 ALARM (Assessing LArge‐scale environmental Risks with tested Methods, GOCE‐CT‐2003‐506675) integrated projects. MBA thanks T. Lahti and R. Lampinen for providing a digital version of Atlas Florae Europaeae . We thank T.D. Mitchell of the Tyndall Centre for Climate Change Research for providing climate and scenario data, and S. Zaehle and M. Erhard at the Potsdam Institute for Climate Impact Research for aggregating climate data.

Journal

Diversity and DistributionsWiley

Published: Jan 1, 2006

There are no references for this article.