Reﬁning Species Traits in a Dynamic Vegetation Model to Project the Impacts of Climate Change on Tropical Trees in Central Africa

: African tropical ecosystems and the services they provide to human society suffer from an increasing combined pressure of land use and climate change. How individual tropical tree species respond to climate change remains relatively unknown. In this study, we reﬁned the species characterization in the CARAIB (CARbon Assimilation In the Biosphere) dynamic vegetation model by replacing plant functional type morpho-physiological traits by species-speciﬁc traits. We focus on 12 tropical tree species selected for their importance in both the plant community and human society. We used CARAIB to simulate the current species net primary productivity (NPP), biomass and potential distribution and their changes in the future. Our results indicate that the use of species-speciﬁc traits does not necessarily result in an increase of predicted current NPPs. The model projections for the end of the century highlight the large uncertainties in the future of African tropical species. Projected changes in species distribution vary greatly with the general circulation model (GCM) and, to a lesser extent, with the concentration pathway. The question about long-term plant response to increasing CO 2 concentrations also leads to contrasting results. In absence of fertilization effect, species are exposed to climate change and might lose 25% of their current distribution under RCP8.5 (12.5% under RCP4.5), considering all the species and climatic scenarios. The vegetation model projects a mean biomass loss of − 21.2% under RCP4.5 and − 34.5% under RCP8.5. Potential range expansions, unpredictable due to migration limitations, are too limited for offsetting range contraction. By contrast, if the long-term species response to increasing [CO 2 ] is positive, the range reduction is limited to 5%. However, despite a mean biomass increase of 12.2%, a positive CO 2 feedback might not prevent tree dieback. Our analysis conﬁrms that species will respond differently to new climatic and atmospheric conditions, which may induce new competition dynamics in the ecosystem and affect ecosystem services.

with so particular worth showed specific responses to past environmental changes [43]. Moreover, with narrower bioclimatic spectrum than PFTs, individual species may prove to be more vulnerable to climate change. The difficulty of DVM application on individual species resides in the importance of specific morpho-physiological traits for optimal model performance [39]; yet, many traits of African tropical tree species remain unknown [44].
In this contribution, we address two questions. First, what is the relevance of using morpho-physiological traits specific to species instead of traits characterizing PFTs? Secondly, what could be the magnitude of tree species responses to climate change in equatorial Africa? We particularly consider changes in species net primary productivity, biomass and distribution, important for evaluating the services provided by plants in the future. With the CARAIB DVM [45][46][47][48][49][50][51], we focus on 12 species which are keystone plants for the ecosystem and for people as goods and services. To take the range of possibilities the future may hold into account, we used several climate projections produced by climate models of the Coupled Model Intercomparison Project-Phase 5 (CMIP5; [52]) to 2100 time horizon under two Representative Concentration Pathways, RCP4.5 and RCP8.5 [53].

Studied Species
Equatorial Africa is dominated by tropical rainforests with broadleaved evergreen or raingreen trees, the deciduous character of the forest increasing with the distance from the equator. A total of 12 tree and palm species were selected according to their great importance for local communities in central Africa, especially for the populations of the WWF (World Wide Fund for Nature) Lake Tele-Lake Tumba Landscape (RD Congo), the largest area of swamp and flooded forest in Africa situated at the lowest latitude zone of the Congo Basin (Table 1; species spatial distributions in Supplementary S1). The uses of their ecosystem services have recently been investigated in this region by the multidisciplinary BIOSERF (Sustainability of tropical forest biodiversity and services under climate and human pressure) project [54]. The selection criteria were: a relatively good knowledge of the species distribution (climatic envelop), the number of species-specific morpho-physiological traits available and the quantity and quality of available data on provisioning ecosystem services, covering a wide range of services. Species occurrences were derived from the RAINBIO (African RAIN forest community dynamics: implications for tropical BIOdiversity conservation and climate change mitigation) database [55], compiling distribution data on tropical African plants (e.g., the occurrences available via the Global Biodiversity Information Facility portal, GBIF) and from Fayolle et al. [56]. The species morpho-physiological trait data used in this study mainly come from the TRY database [57] and from the COBIMFO (Congo Basin integrated monitoring for Forest carbon mitigation and biodiversity) project [58,59] (Table 2). Since they are very diverse across the African continent, the selected ES should be interpreted as potential services. It is consistent with the recommendations of Kumar [60] to distinguish between the ability of an ecosystem to provide a service and the actual use of that service. Some of these services have international economic value, such as the palm oil of Elaeis guineensis and the timber of Lophira alata and Pericopsis elata, but most of them are only appraised at local level. We added three functional types of grasses to represent the lower strata in our model, including herbaceous C3 (wet and dry) and C4 types.

Climate Data and Future Scenarios
The climate database of the Climatic Research Unit (CRU; [61]) was used as reference for historical climate conditions. For future climate, the outputs of 33 general circulation models (GCMs) were acquired from the Coupled Model Intercomparison Project Phase 5 (CMIP5) database [52] of the Intergovernmental Panel on Climate Change (IPCC). We analyzed the capacity of those models to reproduce the annual and seasonal climate over Africa in terms of mean, interannual variability and trend of temperature and precipitation over the reference 1980-1999 period by comparing them to the CRU data [50]. We selected the five most performing models, here presented by decreasing rank: CNRM-CM5 [62], MPI-ESM-LR [63], CESM1-BGC [64], CanESM2 [65] and MPI-ESM-MR [63] (see Supplementary S2). This list of five models was completed by four GCMs producing extreme results in terms of future minimum and maximum precipitation and temperature changes in order to illustrate the range of potential climate change (see Supplementary S3). As such, we obtained under the RCP 8.5 scenario: CSIRO-Mk3. 6.0 [66] for the largest reduction in precipitation (ranked 22nd); CESM1-BGC for the largest increase in precipitation (already among the five best GCMs); CanESM2 for the largest temperature rise (already among the five best GCMs) and INMCM4 [67] for the lowest temperature rise (as all increase; ranked 25th). For each GCM, climatic anomalies between any given year and the reference 1980-1999 period were statistically downscaled to a 0.5 • spatial resolution with a 2D linear interpolation and then combined to 0.5 • CRU data for the reference period, using a procedure similar to that proposed by François et al. [46]. Table 1. The selected species and the ecosystem services they provide as reported by the BIOSERF (Sustainability of tropical forest biodiversity and services under climate and human pressure) project [54] supplemented with data from the PROTA (Plant Resources of Tropical Africa) database (www.prota4u.org). The services are food, fodder, medicinal uses (med), ritual or social uses (rit/soc), environmental uses (envi), poison, construction (constr), materials (mat) and fuel.

Species (Family)
Ecosystem Services food fodder med rit/soc envi poison constr mat fuel Elaeis guineensis (Arecaceae) Musanga cecropioides (Urticaceae) Pericopsis elata (Fabaceae) x x x 9 Pterocarpus soyauxii (Fabaceae) x x x x x 10 Pycnanthus angolensis (Myristicaceae) Symphonia globulifera (Clusiaceae) x x x x x 12 Uapaca guineensis (Euphorbiaceae) To account for the range of possibilities the future may hold, these climate models operated under four different scenarios of greenhouse gas forcing [53] established by the IPCC. Here, we excluded the unlikely strong mitigation scenario (RCP2.6) and concentrated on more extreme scenarios: RCP4.5 (538 ppmv of CO 2 in 2100) and RCP8.5 (936 ppmv of CO 2 in 2100), respectively representing stabilization without overshoot and rise of external forcing.

The CARAIB Dynamic Vegetation Model
CARAIB (CARBon Assimilation in the Biosphere) is a process-based model originally developed to understand the role of vegetation in the global carbon cycle [68,69]. It has also been used to reconstruct potential vegetation distribution under past climates (e.g., [46,70]) and to project land ecosystem evolution in the future [47], including in tropical regions [49,50]. It calculates the carbon and water transfers between the atmosphere, the vegetation and the soil. Depending on the studied area, vegetation is represented by plant functional types (PFTs) or by a set of selected species that are ecologically or economically important. Each species corresponds to a pool of carbon biomass (wood, leaves) per unit area calculated on every grid cell of the simulated domain. It is not an individual-based model. More detailed descriptions of CARAIB modules can be found in previous studies [47,48,71,72]. Raghunathan et al. [49,50] showed evidence of the rather good accuracy of the model to simulate distribution of species in the tropical Brazilian Atlantic Forest.
The environmental inputs of the model are climate data first: the monthly fields of the mean surface air temperature, diurnal temperature range, precipitation, wind speed, relative sunshine hours and air relative humidity; the daily values being determined by a stochastic generator [73]. Other environmental input data include soil characteristics, such as texture and color, the average altitude of the grid cell, and the atmospheric CO 2 mixing ratio (which is assumed to be constant on the simulated domain, but varies from year to year).
The outputs of the model are, among others, carbon fluxes and stocks, soil hydrology and surface energy balance. Those of main interest for this study focusing on provisioning services by plants are the calculation of species net primary productivity (NPP expressed per unit area, g C m −2 yr −1 ) and species biomass (g C m −2 ) by taking interspecific competition for space, light and water into account [47,71]. Besides the impacts of temperature, precipitation or sunshine on plant productivity that are easy to infer, relative humidity and wind speed influence respectively the stomatal regulation and the aerodynamics at the top of the canopy. If air relative humidity is low, the plant will tend to close its stomata to prevent water loss. The stronger the wind, the lower aerodynamic resistance and this allows more efficient carbon and water transfers between soil/plant and atmosphere. Increases in evapotranspiration will tend to decrease soil water and thus plant NPP while increased carbon fluxes will tend to increase NPP.
The model also needs input parameters characterizing the simulated plant species. Whereas previous studies used morpho-physiological traits of plant functional types (PFT; [72]), we looked for characteristics of individual species in an attempt to refine the simulations (Table 2). However, plant trait data available for African species remain scarce as already pointed out by Feeley and Silman [44]. In absence of specific traits, PFT parameters are kept [70]. The selected species belong to two different PFTs: broadleaved raingreen tropical trees and broadleaved evergreen tropical trees. The best-characterized parameters at the species level are the tree height, the foliage carbon/nitrogen ratio (C:N) and the specific leaf area (SLA). The C:N and the SLA modulate the V cmax according to the parameterization from Walker et al. [74]. C:N also impacts autotrophic respiration. The intraspecific variations (standard deviation) of C:N and SLA inferred from available values are presented in Table 2. They are respectively lower and higher than the interspecific variations (5.4 g g −1 for foliage C:N and 0.008 m 2 g −1 C for SLA). The climatic tolerances (minimum temperature, T mins , minimum soil water, SW mins ) governing mortality under stress conditions and the threshold values (minimum yearly sum of growing degree-days above 5 • C GDD5 ming , and, for drought deciduous species, maximum soil water of the driest month, SW maxg ) controlling germination are other species-specific parameters (Table 3). They are inferred from prescribed percentiles in their observed climate distribution, according to the method described by Francois et al. [72] and Laurent et al. [75]. Table 2. The selected species and the morpho-physiological traits for which some values were refined. Grey cells are species-specific values while white cells indicate initial plant functional type (PFT) values. The trait variations (standard deviation) for the foliage carbon/nitrogen ratio (C:N) and the specific leaf area (SLA) are given in addition to the mean value when the number of data is large enough. When species-specific trait data were found, the previous PFT values (standard in CARAIB 1 ) are noted between brackets.   4 Species-specific values from the PROTA database (www.prota4u.org). 5 Average values from the TRY database [57] and its individual datasets [79][80][81][82][83][84][85][86][87][88]. 6 Values from the COBIMFO project [58,59]. 7 Values for Cola sp. Table 3. Species climatic tolerances to stress (minimum temperature, T mins , and, minimum soil water, SW mins ) and requirements for germination (minimum yearly sum of growing degree-days above 5 • C GDD5 ming , and, for drought deciduous species, maximum soil water of the driest month, SW maxg ). Soil water thresholds refer to available soil water in relative units, i.e., in terms of the variable (W-WP)/(FC-WP) where W, WP and FC are respectively the soil water content, the wilting point and the field capacity in mm. In order to evaluate the simulated current distribution of the studied species (corresponding to the reference 1980-1999 period in our study) and the change in distribution between present and future, it is necessary to determine a threshold below which the species is considered to be absent [50]. We used NPP instead of biomass thresholds, because uncertainty is larger for the second due to its dependence on an additional parameter (i.e., tree mortality rate). In analogy with the ENM method of setting a threshold under which one considers the probability of presence to be too low, we set the presence threshold at 100 g C m −2 yr −1 (close to the threshold used in Raghunathan et al. [50]). Anyway, we preliminary observed that varying the presence threshold in the range of 100-200 g C m −2 yr −1 had little impact on the simulated species distributions because NPP sharply decreased below these values, once the climatic conditions are no more suitable. We also chose to maintain the same presence threshold for both PFT-and species-specific morpho-physiological parameters for the sake of comparison. The model validation was conducted using CARAIB outputs obtained with the CRU historical dataset using both PFT and mean species-specific traits. To inspect the ability of the model to simulate the observed species distribution, the actual observations of individual tree species occurrence were superimposed over the simulated species NPP and we calculated the sensitivity (Se, proportion of presence data points correctly predicted by the model). To quantify the uncertainty linked to intraspecific trait variations, two additional simulations were performed replacing respectively the mean species foliage C:N ratio and SLA values by an extreme value (mean ± standard deviation). The gross primary productivity (GPP) from a PFT-based simulation was compared to data at biome level of Beer et al. [89].

Species
An interesting asset of DVM is the potential to test several assumptions concerning the vegetation response to environmental changes, and specifically the CO 2 fertilization effect (as opposed to niche models which lack a direct mechanical connection between the physiology of the species and climate change), known as the increase of leaf photosynthesis with increasing atmospheric CO 2 mixing ratio. However, its long-term effect is uncertain given the interference of other processes on plant growth forcing them to reallocate available resources. Given these uncertainties, CARAIB was driven by climate scenarios under three CO 2 configurations:

1.
A constant atmospheric CO 2 mixing ratio of 330 ppmv ([330] configuration) in the hypothesis that there will be no fertilization effect due to nutrient limitation. Maintaining CO 2 concentrations constant in the model is indeed the only way to test this hypothesis since nutrient limitation in the model is only induced by the species C:N ratio which is a constant input; 2.
An increasing atmospheric CO 2 mixing ratio ([var] configuration) according to the projections of RCP4.5 and RCP8.5 scenarios which respectively reach 538 ppmv and 936 ppmv in 2100. This, then, is consistent with the hypothesis that there will be a fertilizing effect and no nutrient limitation; 3.
An increasing atmospheric CO 2 mixing ratio equal to those of the [var] configuration but accompanied by a downregulation of photosynthetic activity ([down] configuration) in an attempt to take the influence of decreasing nitrogen availability on the CO 2 fertilizing effect into account. Whereas the maximum carboxylation rate (V cmax ) for the first two is calculated as a function of foliage C/N and SLA [74], the third CO 2 configuration consists of a calibration of the photosynthesis based on empirical values of V cmax and maximum rate of electron transport (J max ) measured in free-air CO 2 enrichment (FACE) experiments following the meta-analysis of Ainsworth and Rogers [90]. These experiments were conducted mostly in temperate ecosystems. Their results are thus not directly transferable to the tropical species studied here. Nonetheless, we used them to construct an intermediate scenario of CO 2 fertilization, as a sensitivity test, since such experiments have essentially not been conducted up to now in the tropics [91].
Changes in NPP and biomass projected for the future under each climatic scenario and CO 2 configuration were calculated comparing CARAIB outputs averaged over the 2080-2099 period to the reference 1980-1999 values from the respective scenario.

Refinement of Morpho-Physiological Traits
Substitution of PFT morpho-physiological traits by species-specific ones has an important impact on species NPP with significant differences (p < 0.01) in mean simulated NPP (average over the distribution) for eight species (Figure 1). Both traits linked to water (height, root depth) and carbon (C:N ratio, SLA) uptake were adjusted. When only the tree height is adapted and set higher than the PFT one (Lophira alata, Nauclea diderrichii, Pericopsis elata, Uapaca guineensis), change in species productivity is not significant, except for Lophira alata with a specific height much higher than the PFT value of 20 m. Importantly, increase of tree height leads to a decrease of the mean simulated NPP of at least 50 g C m −2 yr −1 . Turbulence intensity, increasing with tree height, induces an increase in the evapotranspiration rate, which results in an increase of plant water stress and in a decrease of plant NPP, despite a concomitant strengthened assimilation (see Section 2.3). A reduction of root depth has little effect on plant productivity (Elaies guineensis). The most pronounced mean NPP decrease is observed for Ceiba pentendra whose decrease of 150 g C m −2 yr −1 is a result of a much higher height together with a lower root depth and a SLA increase. The most important increases (Cola acuminata, Musanga cecropioides) correspond to a large decrease of foliage C:N when there is a rather limited increase of tree height. The negative impact of increasing SLA is supported by the trait adjustment of Musanga cecropioides and Pycnanthus angolensis. With same slight height increase and big C:N ratio decrease, the NPP increase is less important for Pycnanthus angolensis with a SLA of 0.036 than for Musanga cecropioides with a SLA of 0.026. The large impact of modifying SLA and foliage C:N parameters is a logical consequence of their direct link with autotrophic respiration and V cmax .
impact on species NPP with significant differences (p < 0.01) in mean simulated NPP (average over the distribution) for eight species (Figure 1). Both traits linked to water (height, root depth) and carbon (C:N ratio, SLA) uptake were adjusted. When only the tree height is adapted and set higher than the PFT one (Lophira alata, Nauclea diderrichii, Pericopsis elata, Uapaca guineensis), change in species productivity is not significant, except for Lophira alata with a specific height much higher than the PFT value of 20 m. Importantly, increase of tree height leads to a decrease of the mean simulated NPP of at least 50 g C m −2 yr −1 . Turbulence intensity, increasing with tree height, induces an increase in the evapotranspiration rate, which results in an increase of plant water stress and in a decrease of plant NPP, despite a concomitant strengthened assimilation (see Section 2.3). A reduction of root depth has little effect on plant productivity (Elaies guineensis). The most pronounced mean NPP decrease is observed for Ceiba pentendra whose decrease of 150 g C m −2 yr −1 is a result of a much higher height together with a lower root depth and a SLA increase. The most important increases (Cola acuminata, Musanga cecropioides) correspond to a large decrease of foliage C:N when there is a rather limited increase of tree height. The negative impact of increasing SLA is supported by the trait adjustment of Musanga cecropioides and Pycnanthus angolensis. With same slight height increase and big C:N ratio decrease, the NPP increase is less important for Pycnanthus angolensis with a SLA of 0.036 than for Musanga cecropioides with a SLA of 0.026. The large impact of modifying SLA and foliage C:N parameters is a logical consequence of their direct link with autotrophic respiration and Vcmax. Figure 1. Bar chart of the mean simulated net primary productivity (NPP, g C m −2 yr −1 ) over the observed points of species occurrence, obtained with CRU (Climatic Research Unit) historical climatic data for the reference 1980-1999 period, using respectively PFT-(hatched bars) and species-specific (white bars) morpho-physiological traits (see Table 2). An asterisk behind species name indicates a significant difference (p-value of the Wilcoxon paired t-test < 0.01) between NPP values. Species name between brackets indicates that only height trait has been adjusted.
Trait adjustment impacts the frequency distribution of NPP values over the species spatial distribution. In the case of Pycnanthus angolensis (Figure 2), the frequency of pixels with a NPP below the threshold of species presence (100 g C m −2 yr −1 ) decreases. Above 100 g C m −2 yr −1 , the number of grid cells with low NPP values rises consequently but the frequency in high NPP classes also increases. To deepen this analysis, the impact of parameter adjustment on water use efficiency has Figure 1. Bar chart of the mean simulated net primary productivity (NPP, g C m −2 yr −1 ) over the observed points of species occurrence, obtained with CRU (Climatic Research Unit) historical climatic data for the reference 1980-1999 period, using respectively PFT-(hatched bars) and species-specific (white bars) morpho-physiological traits (see Table 2). An asterisk behind species name indicates a significant difference (p-value of the Wilcoxon paired t-test < 0.01) between NPP values. Species name between brackets indicates that only height trait has been adjusted.
Trait adjustment impacts the frequency distribution of NPP values over the species spatial distribution. In the case of Pycnanthus angolensis (Figure 2), the frequency of pixels with a NPP below the threshold of species presence (100 g C m −2 yr −1 ) decreases. Above 100 g C m −2 yr −1 , the number of grid cells with low NPP values rises consequently but the frequency in high NPP classes also increases. To deepen this analysis, the impact of parameter adjustment on water use efficiency has been examined ( Figure 2). It shows that the NPP increase is a direct result of enhanced water use efficiency due to reduced transpiration water loss. However, as already discussed above, trait refinement does not always lead to an increase in species productivity.
To quantify the uncertainty linked to intraspecific trait variations, two additional simulations were performed replacing respectively the mean species foliage C:N ratio and SLA values by an extreme value (mean ± standard deviation). It results in differences in the mean simulated NPP over the been examined (Figure 2). It shows that the NPP increase is a direct result of enhanced water use efficiency due to reduced transpiration water loss. However, as already discussed above, trait refinement does not always lead to an increase in species productivity.
To quantify the uncertainty linked to intraspecific trait variations, two additional simulations were performed replacing respectively the mean species foliage C:N ratio and SLA values by an extreme value (mean ± standard deviation). It results in differences in the mean simulated NPP over the observed species occurrences, sometimes significant (up to 90 g C m −2 yr −1 for Pycnanthus angolensis with an extreme C:N ratio and up to 170 g C m −2 yr −1 for Pterocarpus soyauxii with an extreme SLA).

Figure 2.
Net primary productivity and water use efficiency of the species Pycnanthus angolensis. (a) Frequency distribution histogram of simulated NPP (g C m −2 yr −1 ), using respectively PFT-(hatched bars) and species-specific (white bars) morpho-physiological traits (see Table 2). Dotted line indicates the 100 g C m −2 yr −1 NPP threshold. Annual mean water use efficiency (g C m −2 mm −1 ), using (b) PFTand (c) species-specific traits. Simulations were obtained with CRU historical climatic data for the reference 1980-1999 period.

Modelling Validation
The model projects relatively well the distribution of the 12 tree species with sensitivity values higher than 90% (Table 4), regardless of adjustment of morpho-physiological parameters. Simulated NPP maps with superimposed species occurrences are available for all the species in the Supplementary S1 and the example of Pycnanthus angolensis is given in Figure 3. We note that all the simulated species have a rather similar pattern of distribution. For Pycnanthus angolensis, we observe an extension of the species distribution when species-specific traits are used, encompassing some observed species occurrences (at the northernmost and southernmost latitudes) which were not considered in the species distribution using PFT traits. Yet more importantly, the increase in high NPP values occurs in regions where species presence is attested by observations. (a) Frequency distribution histogram of simulated NPP (g C m −2 yr −1 ), using respectively PFT-(hatched bars) and species-specific (white bars) morpho-physiological traits (see Table 2). Dotted line indicates the 100 g C m −2 yr −1 NPP threshold. Annual mean water use efficiency (g C m −2 mm −1 ), using (b) PFT-and (c) species-specific traits. Simulations were obtained with CRU historical climatic data for the reference 1980-1999 period.

Modelling Validation
The model projects relatively well the distribution of the 12 tree species with sensitivity values higher than 90% (Table 4), regardless of adjustment of morpho-physiological parameters. Simulated NPP maps with superimposed species occurrences are available for all the species in the Supplementary S1 and the example of Pycnanthus angolensis is given in Figure 3. We note that all the simulated species have a rather similar pattern of distribution. For Pycnanthus angolensis, we observe an extension of the species distribution when species-specific traits are used, encompassing some observed species occurrences (at the northernmost and southernmost latitudes) which were not considered in the species distribution using PFT traits. Yet more importantly, the increase in high NPP values occurs in regions where species presence is attested by observations. Simulated NPP and biomass values are in the range of the values for tropical trees [92]. According to Whittaker [93], annual NPP varies between 800 g C m −2 yr −1 in the tropical seasonal forests and 1100 g C m −2 yr −1 in the tropical rain forests. In a seasonal tropical forest of Mexico, wood and root increments are as low as 285 g m −2 yr −1 of dry matter, i.e., 143 g C m −2 yr −1 [94]. Our GPP results obtained from a PFT-based simulation are also coherent with the ecosystem GPP presented in Beer et al. [89]. For the studied region extended from 16 • S to 16 • N in latitude and from 17 • W to 40 • E in longitude, the correlation coefficient (r) is equal to 0.92 and the root-mean-square error (RMSE) is equal to 510.89 g C m −2 yr −1 . Due to the lack of data for the individual species, we were unable to perform an evaluation of productivity at the species level. The fact that the projected species distributions do not cover all species occurrences could suggest that the model underestimates NPP values, but it could also be a consequence of the calibration of the germination and stress thresholds for the plant species. Table 4. Observed species occurrences, species distribution area (10 6 km 2 ) delimited by the 100 g C m −2 yr −1 NPP threshold and sensitivity (proportion of species occurrences correctly predicted by the model) for both simulations using respectively PFT-and species-specific morpho-physiological traits. Simulations were obtained with CRU historical climatic data for the reference 1980-1999 period.  species-specific traits. The grey color corresponds to simulated NPP values below the 100 g C m −2 yr −1 threshold and the black circles to occurrence data [55,56]. Simulations were obtained with CRU historical climatic data for the reference 1980-1999 period. Table 4. Observed species occurrences, species distribution area (10 6 km 2 ) delimited by the 100 g C m −2 yr −1 NPP threshold and sensitivity (proportion of species occurrences correctly predicted by the model) for both simulations using respectively PFT-and species-specific morpho-physiological traits. Simulations were obtained with CRU historical climatic data for the reference 1980-1999 period. Simulated NPP and biomass values are in the range of the values for tropical trees [92]. species-specific traits. The grey color corresponds to simulated NPP values below the 100 g C m −2 yr −1 threshold and the black circles to occurrence data [55,56]. Simulations were obtained with CRU historical climatic data for the reference 1980-1999 period.

Species
To give a perspective to the model projections for the future, we calculated the change in simulated carbon stock in African tropical forests over historic period for the three CO 2 configurations. The biomass change is calculated in terms of annual change between 1968 and 2007 with the idea of a comparison with field plot biomass studies [24]. A biome classification was derived from PFT distributions (same simulation as for GPP validation). Tropical rainforest and tropical deciduous forest

Climate Scenarios
For central Africa, the climate models project under the RCP8.5 scenario a rise between 3 • C (INMCM4) and 7 • C (CanESM2) of the annual mean temperature between 2080-2099 and 1980-1999 periods. The CNRM-CM5 model, reproducing at best historical climate conditions, projects a rise of maximum 4 • C. For precipitation, the general trend is far less clear, with differences varying between −40% (CSIRO-Mk-3.6.0) and +50% (CESM1-BGC). The CNRM-CM5 model projects an increased rainfall for most grid cells up to +20%, whereas the other models project reduced rainfall in West Africa and an increase in East Africa. Under the RCP4.5 scenario, the same GCM trends are observed though they are limited to only half of the magnitude observed under the RCP8.5 scenario (temperature and precipitation anomalies for all the GCMs and both RCP scenarios in Supplementary S3).

Changes in Soil Water Content
Future soil water availability simulated by CARAIB varies considerably with the climate scenario ( Figure 4) and RCP scenario ( Figure 5). Where, with the GCM ensemble mean, the vegetation model simulates important decreases of soil water levels with respect to the reference period ( Figure 5), it projects rather minor changes with the best-ranked model (CNRM-CM5, Figure 4). Under the scenario projecting the largest precipitation reduction (CSIRO-Mk3.6.0), the model simulates significant soil water reduction in most areas. In case of the largest precipitation increase (CESM1-BGC model), soil water slightly increases in large parts of tropical Africa. Changes in soil water availability are most contrasted for the model projecting largest temperature increase (CanESM2). The GCM ensemble mean clearly shows the difference in soil water induced by RCP scenarios ( Figure 5): when considering a rise of external forcing (RCP8.5), soil water decreases with a mean factor of 0.80, while that decrease is limited to 0.90 when considering a stabilization without overshoot (RCP4.5). extrapolated to the continent, Lewis et al. [24] estimated that carbon storage in the biomass of African tropical forests increased by 90.4 g C m −2 yr −1 from 1968 to 2007 (95% confidence interval: 38.7-116.2 g C m −2 yr −1 ). CARAIB simulated changes are lower but the trends of both configurations including a fertilization effect fall in the confidence interval of Lewis et al. [24].

Climate Scenarios
For central Africa, the climate models project under the RCP8.5 scenario a rise between 3 °C (INMCM4) and 7 °C (CanESM2) of the annual mean temperature between 2080-2099 and 1980-1999 periods. The CNRM-CM5 model, reproducing at best historical climate conditions, projects a rise of maximum 4 °C. For precipitation, the general trend is far less clear, with differences varying between −40% (CSIRO-Mk-3.6.0) and +50% (CESM1-BGC). The CNRM-CM5 model projects an increased rainfall for most grid cells up to +20%, whereas the other models project reduced rainfall in West Africa and an increase in East Africa. Under the RCP4.5 scenario, the same GCM trends are observed though they are limited to only half of the magnitude observed under the RCP8.5 scenario (temperature and precipitation anomalies for all the GCMs and both RCP scenarios in Supplementary S3).

Changes in Soil Water Content
Future soil water availability simulated by CARAIB varies considerably with the climate scenario ( Figure 4) and RCP scenario ( Figure 5). Where, with the GCM ensemble mean, the vegetation model simulates important decreases of soil water levels with respect to the reference period ( Figure 5), it projects rather minor changes with the best-ranked model (CNRM-CM5, Figure  4). Under the scenario projecting the largest precipitation reduction (CSIRO-Mk3.6.0), the model simulates significant soil water reduction in most areas. In case of the largest precipitation increase (CESM1-BGC model), soil water slightly increases in large parts of tropical Africa. Changes in soil water availability are most contrasted for the model projecting largest temperature increase (CanESM2). The GCM ensemble mean clearly shows the difference in soil water induced by RCP scenarios ( Figure 5): when considering a rise of external forcing (RCP8.5), soil water decreases with a mean factor of 0.80, while that decrease is limited to 0.90 when considering a stabilization without overshoot (RCP4.5).
The impacts of the CO2 atmospheric concentrations through RCP4.5 and RCP8.5 climate scenarios and of the CO2 configuration prescribing the species response to increasing CO2 were first analyzed averaging the vegetation model results obtained with the different climate model projections ( Figure 6). The RCP8.5 scenario with a constant [330] (no CO2 fertilization effect) is the most pessimistic scenario for the loss of species range. Considering the species ensemble, the model projects a range contraction of 23.6% but the species may lose up to 30% of their current distribution (Musanga cecropioides). Assuming a CO2 fertilization effect, the contraction remains under 5% of the current distribution. Regarding the expansion of species beyond their current distribution (assuming no dispersal limitation), the scenario assuming a CO2 fertilization under RCP8.5 is the most optimistic (expansion of 55.6% considering the species ensemble). Expansion is at least 20% (up to 85% for Guibourtia demeusei) while it does not exceed 15% under a constant [330]. Reduction and extension of distribution are both less important under RCP4.5 (about the half). The simulations with a variable CO2 mixing ratio but a downregulation of photosynthetic activity do not perform as expected, as the difference with simulations with a variable CO2 mixing ratio and no downregulation is very limited. The differences lie in slightly lower projected expansion beyond the current distribution in the configuration with CO2 downregulation: on average, species range may "only" increase by 28 RCP4.5 and RCP8.5 scenarios. For area contraction, the differences between previous scenarios, respectively 22.6%, 11.5%, 11.1% and 0.5%, are less large. This suggests that the uncertainty in potential future distribution linked to plant response to CO2 is, for most of the species, larger than the one related to the representative concentration pathway (considering here only two RCPs).
The impacts of the CO 2 atmospheric concentrations through RCP4.5 and RCP8.5 climate scenarios and of the CO 2 configuration prescribing the species response to increasing CO 2 were first analyzed averaging the vegetation model results obtained with the different climate model projections ( Figure 6). The RCP8.5 scenario with a constant [330] (no CO 2 fertilization effect) is the most pessimistic scenario for the loss of species range. Considering the species ensemble, the model projects a range contraction of 23.6% but the species may lose up to 30% of their current distribution (Musanga cecropioides). Assuming a CO 2 fertilization effect, the contraction remains under 5% of the current distribution. Regarding the expansion of species beyond their current distribution (assuming no dispersal limitation), the scenario assuming a CO 2 fertilization under RCP8.5 is the most optimistic (expansion of 55.6% considering the species ensemble). Expansion is at least 20% (up to 85% for Guibourtia demeusei) while it does not exceed 15% under a constant [330]. Reduction and extension of distribution are both less important under RCP4.5 (about the half). The simulations with a variable CO 2 mixing ratio but a downregulation of photosynthetic activity do not perform as expected, as the difference with simulations with a variable CO 2 mixing ratio and no downregulation is very limited. The differences lie in slightly lower projected expansion beyond the current distribution in the configuration with CO 2 downregulation: on average, species range may "only" increase by 28 22.6%, 11.5%, 11.1% and 0.5%, are less large. This suggests that the uncertainty in potential future distribution linked to plant response to CO 2 is, for most of the species, larger than the one related to the representative concentration pathway (considering here only two RCPs).
The divergence in future species distribution due to climate model differences is also important.   For almost all the CO 2 scenarios and species, it is with the CNRM-CM5 (the best ranked model) and CESM1-BGC (largest increase in precipitation) scenarios that current species range is the best preserved (and even expanded) while it is with the CanESM2 (largest temperature rise) and MPI-ESM-LR (2nd ranked model) that the species are projected to be the most endangered.

Changes in Biomass
The future expansion of the species beyond their current distribution being hypothetical due to dispersal limitation, we restricted our assessment of change in biomass stock to current distribution. Our results display the uncertainty in the evolution of the tropical forest biomass linked to plant response to CO 2 (Figures 7 and 8). In absence of CO 2 fertilization, the species mean biomass is projected to decrease, whatever the climatic scenario (Supplementary S4). If we consider the climate model ensemble (Figure 7), the biomass change projected for the end of the century is significant for all the species. On average, species might lose 21.2% and 34.5% in biomass under RCP4.5 and RCP.8.5, with at least −17% (Elaeis guineensis) under RCP4.5 to −46% (Lophira alata) and −47% (Pericopsis elata) under RCP8.5. For these latter two species, biomass reductions even reach respectively −76% and −79% under RCP8.5 CanESM2 scenario but are less dramatic (only −11%) under RCP4.5 CESM1-BGC and CNRM-CM5 (Supplementary S4). Assuming a CO 2 fertilization effect, biomasses of the species ensemble increase by 10.1% and 12.2% under RCP4.5 and RCP8.5 (up to 28% for Elaeis guineensis with the climate ensemble and more than 40% with RCP8.5 INMCM4). The exceptions are Guibourtia demeusei and Pericopsis elata whose biomass remains unchanged or slightly decreases in the CO 2 configuration including a downregulation. It is important to note that, despite the CO 2 stimulation of growth, biomass drops (more than 40% with RCP8.5 CanESM2 and 20% with RCP4.5 MPI-ESM-LR) are still projected if we look at individual climatic scenarios. The analysis of the species biomass changes over time ( Figure 8) shows that, despite a CO 2 stimulation, the biomass of most species might stop increasing or even decrease after 2070, especially under RCP8.5 scenario. Under this scenario, the model projects very contrasted biomass changes between the species (as already noticed in Figure 7). preserved (and even expanded) while it is with the CanESM2 (largest temperature rise) and MPI-ESM-LR (2nd ranked model) that the species are projected to be the most endangered.

Changes in Biomass
The future expansion of the species beyond their current distribution being hypothetical due to dispersal limitation, we restricted our assessment of change in biomass stock to current distribution. Our results display the uncertainty in the evolution of the tropical forest biomass linked to plant response to CO2 (Figures 7 and 8). In absence of CO2 fertilization, the species mean biomass is projected to decrease, whatever the climatic scenario (Supplementary S4). If we consider the climate model ensemble (Figure 7), the biomass change projected for the end of the century is significant for all the species. On average, species might lose 21.2% and 34.5% in biomass under RCP4.5 and RCP.8.5, with at least −17% (Elaeis guineensis) under RCP4.5 to −46% (Lophira alata) and −47% (Pericopsis elata) under RCP8.5. For these latter two species, biomass reductions even reach respectively −76% and −79% under RCP8.5 CanESM2 scenario but are less dramatic (only −11%) under RCP4.5 CESM1-BGC and CNRM-CM5 (Supplementary S4). Assuming a CO2 fertilization effect, biomasses of the species ensemble increase by 10.1% and 12.2% under RCP4.5 and RCP8.5 (up to 28% for Elaeis guineensis with the climate ensemble and more than 40% with RCP8.5 INMCM4). The exceptions are Guibourtia demeusei and Pericopsis elata whose biomass remains unchanged or slightly decreases in the CO2 configuration including a downregulation. It is important to note that, despite the CO2 stimulation of growth, biomass drops (more than 40% with RCP8.5 CanESM2 and 20% with RCP4.5 MPI-ESM-LR) are still projected if we look at individual climatic scenarios. The analysis of the species biomass changes over time ( Figure 8) shows that, despite a CO2 stimulation, the biomass of most species might stop increasing or even decrease after 2070, especially under RCP8.5 scenario. Under this scenario, the model projects very contrasted biomass changes between the species (as already noticed in Figure 7).   We also examined the spatial distribution of species biomass changes (examples in Figure 9, maps of all species in Supplementary S5). Changes projected for the 2080-2099 period are evaluated in the light of the temporal variability (standard deviation, σ) of species biomass over the 1980-1999 period. The magnitude of change is expressed as a multiple of the mean biomass standard deviation (average σ over the species distribution; [48]). Unsurprisingly, contraction and expansion of species range go together respectively with biomass collapse (especially for Ceiba pentendra, Lophira alata and Pericopsis elata) and gain. Species disappearance is only projected on the distribution edges. Expansions appear in the mountainous areas in East Africa, around the African Great Lakes and the Ethiopian mountains, where climate conditions become suitable for germination. In absence of fertilization effect, our results show that climatic change might lead to an important forest dieback for most species, even in the preserved part of their distribution, with decrease in biomass two to three times the standard deviation for the reference period. When assuming the existence of a fertilization effect, species biomass substantially increases over most of the current distribution. Some biomass declines are however projected at the outer limits as well as in the core of the distribution. We also examined the spatial distribution of species biomass changes (examples in Figure 9, maps of all species in Supplementary S5). Changes projected for the 2080-2099 period are evaluated in the light of the temporal variability (standard deviation, σ) of species biomass over the 1980-1999 period. The magnitude of change is expressed as a multiple of the mean biomass standard deviation (average σ over the species distribution; [48]). Unsurprisingly, contraction and expansion of species range go together respectively with biomass collapse (especially for Ceiba pentendra, Lophira alata and Pericopsis elata) and gain. Species disappearance is only projected on the distribution edges. Expansions appear in the mountainous areas in East Africa, around the African Great Lakes and the Ethiopian mountains, where climate conditions become suitable for germination. In absence of fertilization effect, our results show that climatic change might lead to an important forest dieback for most species, even in the preserved part of their distribution, with decrease in biomass two to three times the standard deviation for the reference period. When assuming the existence of a fertilization effect, species biomass substantially increases over most of the current distribution. Some biomass declines are however projected at the outer limits as well as in the core of the distribution.

Discussion
We presented here an analysis of climate change impacts on tropical tree species in central Africa using the CARAIB dynamic vegetation model. For this purpose, we refined species parameterization from PFT-to species-specific traits we could find in morpho-physiological trait databases or collect from previous projects, for 12 species (trees and palms) known to be widely used by local populations in central Africa. Our results clearly show the sensitivity of species productivity to these traits and their intraspecific variations with a direction in NPP change depending on the new combination of traits. For instance, the productivity of Pycnanthus angolensis increases while the one of Ceiba pentandra decreases. This indicates that the trait acclimation does not necessarily results in NPP increase, which is rather in accord with Poorter et al. [95]. They concluded that traits at individual level are poor predictors of growth performances while traits strongly vary within species with environmental conditions. Though we did not introduce the traits separately, the analysis of the 12 species modifications (Table 2) allowed to disentangle the effects of each trait on species physiology. Following the formulation of Walker et al. [74], the maximum carboxylation rate (Vcmax) is a function of foliage C:N ratio and specific leaf area (SLA). A decrease of foliage C:N induces both an augmentation of photosynthetic activity and a slighter augmentation of leaf respiration driving at an increased species NPP. A decrease of the specific leaf area corresponds to an increase of the leaf thickness (if leaf density is constant) meaning an increase in the number of chloroplasts. This induces a rise of the maximum carboxylation rate (Vcmax) but the consequently more important leaf biomass (higher leaf thickness) also causes higher respiration costs. The

Discussion
We presented here an analysis of climate change impacts on tropical tree species in central Africa using the CARAIB dynamic vegetation model. For this purpose, we refined species parameterization from PFT-to species-specific traits we could find in morpho-physiological trait databases or collect from previous projects, for 12 species (trees and palms) known to be widely used by local populations in central Africa. Our results clearly show the sensitivity of species productivity to these traits and their intraspecific variations with a direction in NPP change depending on the new combination of traits. For instance, the productivity of Pycnanthus angolensis increases while the one of Ceiba pentandra decreases. This indicates that the trait acclimation does not necessarily results in NPP increase, which is rather in accord with Poorter et al. [95]. They concluded that traits at individual level are poor predictors of growth performances while traits strongly vary within species with environmental conditions. Though we did not introduce the traits separately, the analysis of the 12 species modifications (Table 2) allowed to disentangle the effects of each trait on species physiology. Following the formulation of Walker et al. [74], the maximum carboxylation rate (V cmax ) is a function of foliage C:N ratio and specific leaf area (SLA). A decrease of foliage C:N induces both an augmentation of photosynthetic activity and a slighter augmentation of leaf respiration driving at an increased species NPP. A decrease of the specific leaf area corresponds to an increase of the leaf thickness (if leaf density is constant) meaning an increase in the number of chloroplasts. This induces a rise of the maximum carboxylation rate (V cmax ) but the consequently more important leaf biomass (higher leaf thickness) also causes higher respiration costs. The resulting effect on species productivity is less direct. In humid conditions like in equatorial Africa, a reduced SLA results in improved photosynthetic activity and higher species productivities but, in drier climates, photosynthesis may be limited by stomatal closure.
The vegetation model projects relatively similar distribution patterns for all the species and, despite high sensitivity values, seems to overestimate some of them. A first reason might be an underestimation of actual species ranges, not fully covered by observations. Knowledge of plant species distribution in tropical Africa is improving but it remains limited [55]. However, there are other explanations for these discrepancies. The model projects the potential climatic species distribution ignoring notably past species migration and dispersal that have shaped-and maybe limited-their current pattern [96]. The interspecific competition occurring in the modelling between the 12 studied species (as well as with the herbaceous strata) does not integrate all the aspects of competition between the plants and does not take the other biotic interactions like parasitism and herbivory into account, which might affect species NPP [48]. The rather coarse 0.5 • spatial resolution of the simulations and of the gridded historical climate dataset [61] derived from monthly observations at meteorological stations, still very few in Africa, are other sources of uncertainties. The limited field data [44,97] make also difficult the validation of the simulated productivity and carbon stocks at the species level. We showed that the mean grid cell GPPs simulated by CARAIB are comparable to the observation-based estimate of ecosystem GPP presented in Beer et al. [89]. We also examined the change in carbon storage in biomass projected by CARAIB for recent period (between 1968 and 2007) in the light of the analysis of Lewis et al. [24] based on long-term monitoring plots in intact (no human disturbances like deforestation) African tropical forests. The model simulates an increase in tropical forest biomass as observed in long-term plots but only if plants are assumed to respond positively to increasing atmospheric carbon dioxide concentrations (CO 2 [var] and [down] configurations). However, the projected increase is clearly less important than reported by Lewis et al. [24], even if it falls in their confidence interval. Thus, our results suggest the existence of a fertilization effect. But, they also imply that the increase of atmospheric CO 2 cannot explain alone the increase in forest biomass in recent years. Recovery of long-lived trees from past anthropogenic or natural disturbances, a component difficult to take into account in modelling but also hardly monitored, may be another important explanation for the observed biomass increment. With a remote sensing approach, Liu et al. [98] estimated for 1993-2012 only a slight carbon increase over the intact tropical forests. Combined with a considerable deforestation (not included in our simulations), the net effect corresponds to a decline of 0.01 PgC yr −1 for the African tropical forests and more important decreases of 0.16 and 0.09 PgC yr −1 respectively for American and Southeast Asian tropical forests (0.26 PgC yr −1 loss of above-ground biomass carbon for global tropical forests).
The potential CO 2 effect in the long term is however uncertain given the interference of other processes on growth, forcing plants to reallocate available resources. Smith et al. [99] highlighted the current lack of long-term responses of vegetation in models while modelling studies can especially broaden our knowledge of these processes [100]. Among existing models that incorporate dynamic algorithms, there are many differences between these algorithms and between the types of responses considered. Some describe acclimation (mid-term responses, taking days to occur) of respiration to temperature and other acclimations of photosynthetic activity to temperature or atmospheric CO 2 mixing ratio. At present, no model integrates these three types of response simultaneously. In this study, we attempted to implement in our model an empirical algorithm based on FACE (Free-air CO 2 Enrichment) experimental data [90] to simulate an acclimation of photosynthesis to increased atmospheric CO 2 mixing ratio. To date, FACE measurements are the most appropriate data despite their limitations: experiments conducted at the PFT level in temperate or desert zones (without tropical species) with a duration limited to a few years and disputed causal mechanisms [99]. Therefore, our results with a downregulated CO 2 response must be considered as a first attempt and interpreted with caution (see below). Further research is recommended to implement photosynthesis downregulation in vegetation models, for instance by integration of a kind of optimization provided by the plant when photosynthetic activity tends to increase. Optimizing allocation and utilization of nutrients allow the plant to reach a compromise between the advantage of having a high Rubisco activity and the cost of maintaining it by respiration (e.g., [101]).
In our results, the differences between the climatic scenarios from seven GCMs under two RCPs and three CO 2 configurations highlight the large uncertainties in the future of African tropical tree species. They are first related to uncertainties in the climatic scenarios [102]. The GCMs selected for our study vary greatly in their projections of future precipitation although several of them project increases in precipitation in eastern Africa and decreases in western Africa likewise in James et al. [103]. Soil water content anomalies (increasing factor) projected by the model show that the temperature rise emphasizes soil water deficits, even in areas where climatic models project unchanged or increasing rainfall. Projected changes in species distribution vary greatly in function of the GCM climatic anomalies for the 21st century (uncertainties of 29.4% in the lost fraction and 41.8% in the area expansion beyond current distribution) and, to a lesser extent, of the concentration pathway, considering only two RCPs (differences of 11.1% between the [330] RCP4.5 and RCP8.5 scenarios for the area contraction and 24.3% between the [var] RCP4.5 and RCP8.5 scenarios for the range expansion). The question about long-term plant response to increasing CO 2 concentrations also leads to important differences in simulated future distributions. Our simulated range contraction and expansion for the species ensemble vary respectively of 22.6% and 46.8% between the [330] and [var] RCP8.5 scenarios. It is confirmed by Rammig et al. [104] in a study evaluating the risk of dieback in the Amazonian forest. The availability of the nutrients also raises questions as in Cusack et al. [27] and Huntingford et al. [105] studies which concluded that the largest uncertainties are associated with physiological and nutrient cycling responses of plants.
In absence of the fertilization effect, in our simulations (CO 2 [330] configuration), species are exposed to climate change and might lose under RCP8.5 25% (12.5% under RCP4.5) of their distribution on average (mean value for the ensemble of both species and climatic scenarios). The most vulnerable species seems to be Musanga cecropioides (30% of lost area) followed by Cola acuminata and Guirboutia demeusei (27%). Projected range expansions are also too limited to offset range contraction and moreover highly unpredictable. By contrast, if the long-term species response to increasing atmospheric [CO 2 ] is positive (CO 2 [var] configuration), the simulated future expansion might largely compensate for a range reduction limited to 5%. But, as already mentioned, the extent of expansion beyond current distribution is extremely dependent on the scenario, e.g., from 20% to 85% for Guibourtia demeusii. The simulations with CO 2 [down] configuration exhibit some downregulation, especially apparent in the area expansion less important than projected under CO 2 [var] configuration, but the effects of this acclimation are more limited than expected. In these two configurations, the current distribution of the 12 studied species has a good chance to be preserved. These results are consistent with those produced by other mechanistic models for tropical areas (e.g., [106]). According to Lloyd et al. [107], tropical forests are not "dangerously close" to their thermal optimum since the new temperature and atmospheric CO 2 mixing ratio conditions may produce higher growth rates [108]. Nevertheless, our findings also highlight that, if species range might be totally or partly preserved, there are large uncertainties on future biomass stocks. In absence of a long-term CO 2 fertilization effect, CARAIB projects significant dieback (i.e., amplitude of biomass loss two to three times the biomass variability calculated for the reference period) over all the conserved distribution. Considering the species ensemble, it means an average biomass loss of −21.2% under RCP4.5 and −34.5% under RCP8.5. Despite a mean biomass increase of 12.2%, assuming a positive feedback of increasing CO 2 does not prevent tree dieback both at the outer limits and in the core of the preserved distributions. This shows the competing effects that climate change and CO 2 fertilization might have [109,110].
This zoom on 12 individual species confirms that individual species respond differently to new climatic and atmospheric conditions, which may induce a new dynamics of competition in the ecosystem. For example, whereas the model projects under RCP8.5 CO 2 [var] configuration that mean biomass might increase by as much as 28% for Elaeis guineensis, it simulates unchanged biomass stock for Guibourtia demeusei and Pericopsis elata. Under RCP8.5 CO 2 [330] configuration, the three deciduous species (Ceiba pentandra, Lophira alata and Pericopsis elata) of our set are projected to experience the strongest dieback (around 50% of biomass loss over the current distribution). The impact differs also significantly between smaller study areas. For instance, over the WWF Lake Tele-Lake Tumba Landscape in DR Congo, a forest-savanna mosaic with strong spatial variation in forest biomass [111], the currently present species would not profit off the new conditions; they would sustain but not amplify current production rates. On the contrary, previously (nearly) absent species would appear and amplify their production rates exponentially. Such a species turnover would not be a favorable situation since those species not necessarily have the same role in the ecosystem nor provide the same services. Today, countryside way of life still partly relies on conventional plant uses.
Finally, we want to notice that the traits values used in this study are measured under current climate conditions and that, under influence of climate change, acclimation of these traits may occur. For example, Santiago et al. [112] highlight that rooting depth could be a major determinant of how tropical trees will fare during extreme droughts. Assumptions on the physiological responses of species to new climatic conditions thus remain necessary [49]. These results must also be analyzed keeping in mind the absence of anthropogenic actions and species dispersal/migration limitations. Synergistic interactions and feedbacks of clearing and forest fragmentation certainly have major implications for the vulnerability of forests to climate change. While the precise sensitivity is not known, human activities appear to increase the climate sensitivity of forest areas that are not directly affected by these activities, causing feedbacks on global and regional climates [30,102]. At the same time, Scheiter et al. [113] showed that appropriate ecosystem management can mitigate climate change impacts on vegetation and delay or avoid undesired vegetation shifts. We thus emphasize the potential gain of studies that examine human and climate pressures simultaneously (e.g., [114]). The development of approaches that integrate two spatial and dynamic models, one simulating vegetation and the other simulating land use by human agents (multi-agent model), should be encouraged to expand our knowledge of the complex interactions between human beings and ecosystems. The call for the integration of perspectives of natural sciences with those of social sciences [115,116] is then obvious. The question also arises whether the species would be able to migrate in order to reach newly suitable potential areas. In the WWF Lake Tele-Lake Tumba Landscape, the transport of zoochoric species seeds is clearly limited by hunting [117] and more widespread fires [102]. Other human activities such as slash-and-burn agriculture, road construction or habitat fragmentation may also represent barriers or challenges to migration, though the impact of anthropic actions in determining the ability of plant populations to disperse remains poorly understood [118].

Conclusions
Our study with the CARAIB dynamic vegetation model shows the specific response of tree species to climate change, urging for impact assessments carried out at the species level in parallel with ecosystem analyses. In Africa more than anywhere else, the lack of observations for species inventory and measurements of morpho-physiological traits limits thus our knowledge of species physiology and distribution. CARAIB results highlight the large uncertainties in the climate projections for the end of the century and in the long-term physiological responses to CO 2 concentration increase. The model projects contrasted future species distributions and functioning assuming or not a stimulation of plant growth by CO 2 , with NPP and biomass decreases (more or less pronounced) regardless of the long-term fertilization hypothesis.