Changes of the potential distribution area of French Mediterranean forests under global warming

. This work aims at understanding future spatial and temporal distributions of tree species in the Mediterranean region of France under various climates. We focused on two different species ( Pinus Halepensis and Quercus Ilex ) and compared their growth under the IPCC-B2 climate scenario in order to quantify signiﬁcant changes between present and future. The inﬂuence of environmental factors such as atmospheric CO 2 increase and topography on the tree growth has also been quantiﬁed. We modeled species growth with the help of a process-based model (MAIDEN), previously calibrated over measured ecophysiological and dendrochronological series with a Bayesian scheme. The model was fed with the ARPEGE – MeteoFrance climate model, combined with an explicit increase in CO 2 atmospheric concentration. The main output of the model gives the carbon allocation in boles and thus tree production. Our results show that the MAIDEN model is correctly able to simulate pine and oak production in space and time, after detailed calibration and validation stages. Yet, these simulations, mainly based on climate, are indicative and not predictive. The comparison of simulated growth


Introduction
Forest productivity and distribution changes in the coming decades are a major concern of ecological studies.Differential impacts of global climatic changes or increasing atmospheric CO 2 concentration, as well as complex vegetation feedbacks on climate makes it hard to estimate the vegetation patterns in future (Prentice and Webb, 1998;Sitch et al., 2003).It is particularly true in Mediterranean region which has a limited area but with a high diversity of environments, vegetation and fauna (Joffre et al., 1999).Forests in this region are crucial for preserving this biodiversity and provide essential ecosystem services, such as soil protection, water conservation and climate regulation (Eamus et al., 2005).Circulation models predict a significant warming and decrease in precipitation in the Mediterranean basin (Gibelin and Deque, 2003).Giorgi (2006) identified the Mediterranean basin as the most prominent hot-spot of climate change over the world with more than 20% decline of precipitation over April to September.Precipitation and drought risk changes can trigger a positive feedback on climate change by decreasing CO 2 sequestration in ecosystems (Ciais et al., 2005).
Assessment of climatic change impacts must be supported by accurate and reliable models of climate-vegetation growth relationships (Misson, 2004;Misson et al., 2004;Rathgeber et al., 2003).These models must be carefully calibrated and validated on data at different time and space scales.At large spatial scales, these data may come from field-based studies on natural gradients (Rathgeber et al., 2005).To Published by Copernicus Publications on behalf of the European Geosciences Union.
include the time variability at these gradient, dendrochronological series are valuable as an exploratory tool for identifying relevant climatic variables and periods when tree growth is responding to climate (Tessier, 1989) or as a support to modelling (Misson, 2004;Rathgeber et al., 2005Rathgeber et al., , 2000)).At smaller space and time scale, field data (Rathgeber et al, 2005) and/or coupled micro-meteorological and biogeochemical data measured at fine resolution in forest stations (Misson, 2004) are essential to understand eco-physiological processes of the water and carbon cycles in the ecosystem.
Modelling approaches are diverse.Some of them have focused either on ecophysiological processes (Anselmi et al., 2004) or on spatial patterns by the use of statistical relationships (Luoto et al., 2005;Thuiller et al., 2004).We intend here to simultaneously investigate both issues.We explore the geographical distribution of two common Mediterranean species: Aleppo pine (Pinus halepensis M.) and evergreen Holm oak (Quercus ilex L.) by using an ecophysiological model (MAIDEN, (Misson et al., 2004) that has been calibrated independently (Gaucherel et al., 2008) for both species using a Bayesian approach.This model is driven here, for the South-East France region, by climate simulations for the 21st century, at the daily time step using the Global Climate Model ARPEGE for the IPCC-B2 scenario (Deque et al., 1994).In addition, the fertilization effect of CO 2 is assessed by comparing simulations with constant and with variable CO 2 concentration.
Our main assumptions for this work are common to other similar works (Yu et al., 2002): forests are mature, they do not genetically adapt to climatic changes, their colonization process is quick enough relatively to climate changes to neglect their population dynamics.Moreover, we only focus on potential distribution areas, neglecting human impacts (Debussche et al., 1999).The objectives of this work are threefold: i) to estimate the tendency of both Mediterranean species growth for the 21st century and to compare them to the past (20th century) in order to detect significant changes; ii) to quantify the influence of environmental factors such as CO 2 increase on the tree growth according to the topography, the other external factors, such as soil composition, being considered as fixed and uniform over the entire region; iii) to quantitatively compare both species behavior regarding to climate/topography/ CO 2 concentration changes, thus helping to extrapolate our results.

Dendrochronological and ecophysiological data
We use dendrochronological data from 21 Aleppo pine (P.halepensis) stands (Nicault, 1999) and one typical Holm oak (Q.ilex) stand (unpublished data sampled in the Puéchabon site (Rambal et al., 2003)).All sites are located on calcareous soil in PACA region (Long.0-8 • N×Lat.41.5-46.5 • N), southeast France.The pine tree-ring series are based on the annual earlywood width (E w ), the latewood width (L w ), the earlywood density (E d ) and the latewood density (L d ), which are combined into a biomass index, W (t) =E w (t) .E d (t) +L w (t) .L d (t), taking into account for both stem radial growth and wood density (Rathgeber et al., 2000).The W (t) curve of oak is calculated on the ring widths only as density measurements are not available.Ring width series are indexed using a digital filter as it is usually done in dendroclimatology (Cook and Kairiukstis, 1990).The oak tree-ring chronology was built by sampling 15 stem disks from 15 trees, among which only seven were kept.Because of the eccentricity of the stem, only the longest radius was measured in each disk.The tree-ring width chronologies on these 15 disks were interdated 2 by 2 in order to give the exact year for each tree ring.Ring with very high and very low width help us interdating the chronologies, as well as rings showing some cell disruption because of very low winter temperature that damaged the cambium (1963, 1985 and 1987).Tree-ring chronologies have been averaged into a 32-year regional chronology  in order to avoid over-parameterization and to perform a rigorous calibration.These 32 years are the common years for all sites.The 38year oak series  is more local as it is based on a single site.Both mean series will be used to calibrate some parameters of the MAIDEN tree-growth model.Botanical, ecological and topographical factors were also recorded at each stand (Rathgeber et al., 2005) and used to help in calibrating several physiological processes.Two years, 1956 and1985, are characterized by extremely low winter temperature for the region (much below the freezing level), which induced partial damaging of the cambium and had strong impacts on the growth during three years.As these processes are not integrated in the model, years 1956-1958 and 1985-1987 have been removed from the dataset.
In addition, we used the averaged biomass index of each 21 stands over the first 50 years (R50) of each site as an index of the unbiased fertility of the stand.R50 is computed by summing W(t) of each tree over the first 50 years and then averaging them over the stand.Although this index could still be influenced by stand density, it offers a synthetic idea of the spatial distribution of pine production (Rathgeber et al., 2005).At two sites, one in P. halepensis and one in Q. ilex, tree transpirations were measured continuously at 1 / 2 h time step for one year using the heat dissipation method (Granier, 1987).Xylem sap flux density (H 2 O h −1 dm −2 sapwood) was monitored on 6 trees in the two plots.The sampled trees were representative of the plots' mean basal area tree circumference.The sap flow sensors consist in a pair of probes, 2 cm long and 0.2 cm in diameter each, inserted in a radial orientation separated vertically by 12 cm behind the cambial zone.The temperature difference between the heated and reference probes ( T ) was recorded and sap flux density was calculated according to Granier (1985).Up-scaling from tree to stand was based on the estimation of the stand sapwood area.Biogeosciences, 5, 1493Biogeosciences, 5, -1504Biogeosciences, 5, , 2008 www.biogeosciences.net/5/1493/2008/

Ecophysiological model
The process-based tree-growth model MAIDEN is extensively described in Misson (2004); Misson et al. (2004); Vincke et al. (2005).The model calculates processes such as photosynthesis, stomatal conductance and carbon allocation.The water balance is computed at the ecosystem level, including canopy water interception, transpiration, soil evaporation, soil water transfer, drainage and runoff.MAIDEN separates daily net primary production (NPP) between carbon reservoirs (leaf, bole, root and storage) according to phenological phase-dependent rules.These phases are (1) winter: no activity, (2) spring: leaf and root expansion, (3) summer: bolewood production, (4) early falls: carbohydratereserve accumulation, (5) late fall: leaf and root senescence.
An original modeling procedure of carbon storage and mobilization was developed to reproduce the temporal autocorrelation structure of tree-ring series (Guiot, 1986).The annual increment of bole carbon reservoir at stand level is the modeled variable that will be compared with the observed dendrochronological tree-ring series.To be applied to a given species at a given site, the model needs several input variables such as altitude, latitude, maximum absolute Leaf Area Index (LAI), specific leaf area, initial bole biomass, soil thickness and soil textural classes, which can be obtained from site measurements.Moreover, it also needs eleven internal parameters that can be tuned to fit at best available ecophysiological and dendrochronological data, as explained in Gaucherel et al. (2008).Climatic driving variables are daily minimum and maximum temperature and precipitations.MAIDEN uses the MT-CLIM algorithm (Running et al., 1987) to estimate daily solar radiation and vapor pressure deficit from daily temperature and precipitations.

Climatic scenarios
Observed daily temperatures and precipitation series are obtained from the METEO-France meteorological station of Aix-en-Provence.Climatic simulations for the 20th and 21th centuries are obtained from the global climate model ARPEGE-IFS (hereafter APG (Gibelin and Deque, 2003)) driven by the IPCC-B2 scenario radiative forcing (including greenhouse gases, CO 2 , CH 4 , N 2 O, CFC, water vapor, ozone, and five types of aerosols).Doubling of atmospheric CO 2 concentration occurs towards the end of the 21st century with 610 ppm, while the 20th century concentrations are the observed values.Description of the model is given in Deque et al. (1994) and Gibelin and Deque (2003).The APG model, for which the time step is 30 min, provides daily maximum and minimum air temperatures as well as daily precipitations for the 1960-2099 period.The grid point used has a 0.5 • ×0.5 • cell size that we have interpolated as a function of latitude and longitude using a two-dimensional (2-D) bicubic technique.
We used the Climate Research Unit (CRU) gridded dataset (New et al., 1999) over the 1960-1998 control period to estimate possible biases of the APG simulation and correct it during the 21st century.This surface climatology has also a 0.5 • ×0.5 • resolution and has been interpolated for visualisation comfort in the same way than the APG fields.Although having the same resolution, the CRU grid does not perfectly superimpose with the APG grid: to obtain comparable data over the whole region, temperature and precipitation fields have been systematically corrected by associating the nearest grid points between the two grids (thus without any interpolation).The maximum and minimum temperature and precipitation fields simulated by the APG model are finally compared with the 1960-1998 CRU climatology over South East France (New et al., 1999).In addition, the regional topography is available with gridpoints ranging from 0 to 2200 m elevation in Pyrenees and Alps mountains.

Ecophysiological calibration and simulation
Bayesian techniques such as Monte Carlo Markov Chains (MCMC) offer the opportunity to estimate the (a posteriori) distribution of the calibration parameters, starting from their predefined (a priori) statistical distributions (Andrieu et al., 2001;Gelman et al., 1995).We basically proceeded in two main stages (details on the calibration are given in Gaucherel et al., 2008).The a priori distributions are assumed to be uniform with a very large range.First, the range of six parameters controlling the ecophysiological processes has been narrowed by using the MCMC algorithm and the transpiration observations available for both species.This method provided us with an a posteriori probability distribution and then a 90% confidence interval for each parameter.The limits of this interval have been used to define the a priori limits used in the second stage in which the eleven parameters (the six ones calibrated in the first stage, plus the five remaining ones related to carbon allocation processes) have been calibrated on the observed dendrochronological time series.For each calibration, adapted convergence tests proposed by the BUGS project (Spiegelhalter et al., 1993) have been implemented.We performed detailed elasticity tests of the model to estimate possible misspecifications of model's components, as well as validation on independent meteorological data to detect possible over-parameterization (Gaucherel et al., 2008).Finally, the model fit (determination coefficient r 2 ) between observations and simulation modes, after five draws of 5000 steps each, converges and summarizes the modeling efficiency.
The eleven parameters are species dependent and valid for the whole region.topography, elevation and coordinates (for radiative process estimation).We choose averaged regional soil textures for the four soil layers (being 10, 20, 20 and 10 cm deep respectively): percentages of clay (25, 30, 35, 35 resp.) and sand (10, 15, 20, 20 resp.) (Misson et al., 2004).Ecophysiological simulations are made for mature stands such that the only biomass gain between years is caused by the bole carbon allocation.The CO 2 atmospheric concentration has an influence on the climate and is taken into account by the climate models, but it has also a direct effect in the tree-growth model, called also fertilization effect because it stimulates the photosynthesis and reduces the water loss through the stomata.
3.2 Regionalization APG simulations cannot be directly used as inputs in the MAIDEN model because they are affected by some systematic biases.They have been corrected with the CRU data at every gridpoint and for every month of the studied period, by using the temporal discrepancies between both data sets over the common period .APG temperatures (resp.precipitations) are corrected by adding to daily values the difference (resp.ratio) between APG and CRU means of the corresponding month.No topographical correction has been applied on APG variable values as both data sets present similar topography.

Simulation validation
R50 index and simulated bole increment for pine species will also be compared in order to assess our model ability to simulate spatial variability of tree-growth, after interpolation using a 2-D bicubic interpolation for southeast France.The assessment of the climate impact on tree-growth is based on the comparison between two periods of 38 years each: 1960-1997 (hereafter called 20th century period) and 2062-2099 (21st century period).The direct impact of CO 2 is based on the comparison on the 21st century with and without taking into account for the CO 2 increase at the input of the MAIDEN model.Runs without CO 2 increase means that the 1995 CO 2 concentration (360 ppm) has been used for the whole 21st century.The vulnerability of both species is assessed by comparing the maps of their simulations for both periods and the time series of their averaged bole increments.

Ecophysiological and dendrochronological calibration
The differences (mean absolute values in mm/day) between observed and simulated transpiration over year 2004 is low and reach 4% for pine and 6% for oak (not shown).This discrepancy does not seem to be due to the parameter values as both observation and simulation curves present the same seasonal variations, but rather due to a higher soil depth inducing a lower simulation of the soil water stress.The bole C allocation of the pine is around 650±50 g/mm 2 and around 250±40 g/mm 2 for the oak.
The pine MCMC calibration of MAIDEN model with PACA dendrochronological series shows a modal fit at r 2 =0.37±0.04 (Fig. 1).This correlation is lower than the one computed in previous works (Misson, 2004), because model has been improved and because the (Bayesian) calibration stage is different and more rigorous (Gaucherel et al., 2008)   In these circumstances, the model seems to show too much interannual variability and to be unable to solve for low frequency growth cycle.Nevertheless, differences between calibrated bole increment simulations and observations are quite homogeneous, showing the ability of the model to simulate a large variability of growth (Fig. 1), while the oak growth is better fitted than the pine growth.It can be due to the negative trend of the oak growth observed since the beginning of the year 70's.The fact that the model is able to simulate this trend could confirm its climatic origin.As additional and independent validation, gridded pine R50 values and simulated bole increment are compared and show a correct qualitative agreement (Fig. 2).Maxima are located at the centre of both maps, even if R50 have two maxima and simulation maximum is slightly shifted westward by about 30 km.

Climate regionalization
Temperature and precipitation data for CRU and APG climatology are rather similar over the 38 common years in average (Table 1), but the APG simulations are more variable from month to month.Both temperature data sets show a lapse rate of about −0.0065 K/m.Spatial temperature and precipitation simulated distribution (not shown) are also more scattered, probably due to the fact that each APG gridpoint is a punctual value rather than a cell-grid average.The mean annual deviation (resp.ratio) between APG and CRU datasets over this control period is +3.18±1.68K for mean temperature (resp.0.82±0.47 for rainfall) and for the whole studied region.Note that these averaged values are not exactly equal with differences between lines of Table 1, as calculations before averaging stage are made separately for each grid points and each month.).The modern limit of its potential distribution is indicated in black dashed line.(b) Map of the mean annual growth variation of pine in southeast France (in %), between the twentieth and the twenty first centuries without direct effect of CO 2 (i.e. with atmospheric CO 2 being fixed to 360 ppm).The 21 pine stand area (Fig. 2) is located on the bole increment map.
APG simulations for the 21st century are corrected using monthly and spatially dependent deviation (or ratio) coefficients (Table 2).Temperature (Fig. 3a) and precipitation (Fig. 3b) values of the 2062-2099 period are significantly higher than those of the 1960-1997 period.Mean annual deviations of daily maximum and minimum temperatures are +2.62±0.66K and +2.27±0.54K, respectively, while mean annual precipitation will increase by a factor 1.75±0.98.The temperature diurnal amplitude (max-min) will also decrease with elevation.Precipitation patterns between both periods are very different at high elevations (almost a factor three) and similar at low elevations (Fig. 3b).
Table 2. Correction coefficients (up) of APG parameters for the twenty-first century based on monthly and spatially dependent deviation (for temperatures) or ratio (for precipitations) coefficients.All coefficients have been computed over the

Pine growth
Pine growth simulated for the 20th century with a CO 2 concentration increasing from 317 ppm (in 1960) to 360 ppm (1995), as expected, shows less productive zones at higher elevations where temperature is less favorable (Fig. 4a).For the 21st century, there will be a negligible mean decrease of about −2.72±4.57%compared to the first period.This hides a large spatial variability with increase up to +8% at high elevations (Fig. 4b).If we take into account the direct effect of atmospheric CO 2 (rising up to 612 ppm at the end of the 21st century), growth pattern will be very different (Fig. 5a): a mean increase of about +22.24±4.14%,with maxima at about +30% at intermediate elevations.Moreover, in relative values, low elevations are the most favored by the direct fertilizing effect (Fig. 5b).The average gain due to CO 2 is about +26.01±3.68%with a significant relative productivity loss of about 0.8% per 100 m elevation (up to the highest grid-point at elevation 2200 m).This difference between pine production simulation with and without direct CO 2 effect is linear in time and of about 0.34% per year.With CO 2 direct effect too, pine production is maximum at intermediate elevation, where the climate is a transition between warm and dry plains and cool and wet mountains (Fig. 5a).Pine growth seems to be optimum nowadays for elevation comprised below 800 m (Fig. 6a), while this optimum is shifted up to 1100 m for the 21st century (Fig. 6b), whatever is the CO 2 factor used.The twenty-year window moving average curves (Fig. 7a) show that pine growth reaches a maximum for present if we keep CO 2 constant during the 21st century, while it continues to increase if CO 2 increases according to the IPCC B2 scenario.The productivity increase seems to be constant, or weakly slowing down, at the end of 21st century (Fig. 7a).This increase over 100 years reaches 30% in average for the region, with a maximum at intermediate elevation.The present distribution of both species only covers a small fraction of these potential areas (approximately 15% for pine and 25% for oak) (Figs.5a and 8a, dashed lines).

Oak growth
Simulations of oak growth for 20th century without CO 2 direct effect show a weak spatial variability and much lower (two times) productivity than pine (Fig. 8a).Oak mean growth simulation during the 21st century period decreases by about −13.8±4.8% compared to the control period, but reaches a value of −25% in some low elevation zones (Fig. 8b).The differences between both periods are much more important for oak than for pine.If we include the direct effect of atmospheric CO 2 (Fig. 9a), mean oak growth increases by about +24.1±5.6% and reaches a value of about +43% at high elevations.The spatial distribution is quite homogenous regarding to the elevation factor.Oak does not show any optimum with elevation (Fig. 6c, d).The gain due to CO 2 direct effect is optimum at low elevations (Fig. 9b).
As highlighted by the twenty-year window moving average curve, the trend of regional oak growth with and without CO 2 direct effect (Fig. 7b), has the same general behavior than for pine: a progressive productivity increase throughout the 21st century with CO 2 direct effect versus a decrease with constant CO 2 .This increase over 100 years reaches 50% in average for the region, and even more at low elevation.Tree growth evolution for constant CO 2 concentration reaches a maximum around year 2005 for both species (Fig. 7).The simplified drought index of Fig. 6c (temperature/precipitation ratio) smoothed with a twenty year moving window confirms that the period around 2005 is simulated as the most humid of the analyzed period.Afterwards, drought increases considerably, inducing an important growth deficit if CO 2 direct effect does not compensate it.

Climate and growth modelling
The distribution of precipitation corrections (computed between CRU and APG data sets for the twentieth reference period) is highly heterogeneous in space (Fig. 3b).This is possibly due to difficulties to precisely simulate frequent and local raining events in Mediterranean climate (West et al., 1986).The precipitation changes are only significantly positive above 500 m, for the benefit to the tree growth.This is obvious for Alps and in a much lower extent for Pyrenees mountains.For temperature, Alps and Pyreneans do not behave in the same way: temperature anomalies will decrease with elevation in the Alps, and increase in the Pyreneans (not shown), amplifying so the constraints.The MAIDEN model, after calibration, has shown a correct ability to simulate tree growth of pine and oak species under various environmental conditions (Fig. 1).Some previous works succeeded in such modeling using remote sensing (LAI) calibration (Anselmi et al., 2004).Our approach is, in a way, different and original, as we have calibrated a complex ecophysiological model, using dendrochronological time series (Gaucherel et al., 2008).The simulations are generally realistic from 1960 to 1997.But some short periods (such as 1985-1987) are poorly simulated, mainly due to factors not included in the model, as cold winters that could have damaged the cambium cells of tree bole.The model calibration is based on the mean regional tree-ring chronology, and an independent validation has been designed on the spatial variability of the tree-growth.The spatial distribution of the pine R50 indices has been successfully compared to the spatial distribution of the simulations.Tree growth simulations are only regional trends and are still too broad to take into account for detailed abiotic factors such as soil characteristics.Nevertheless, spatial and temporal trends have been clearly identified.If we do not take into account for CO 2 fertilization effect, we find that tree growth has reached a maximum around year 2005 for both species (Fig. 7): this should not be taken as a precise date for the optimum of these species; it rather shows that the climate changes will not beneficiate to these species in the future in terms of ecophysiological processes at large scales in average.
Some assumptions made in this work have to be discussed.Simulated distribution areas are potential distributions, which do not taken into account any history of species colonization nor human influences (Debussche et al., 1999).The present distribution of both species, only covering a small fraction of their potential areas (Figs.5a and 8a, dashed lines), will likely extend or shift during the coming decades, as it has been simulated for example for broadleaf forests in China by similar approaches (Yu et al., 2002).For clarity too, we have chosen not hiding on maps high elevations at which tree species would not be able to grow due to freeze.Furthermore, another hypothesis behind our results is that there is no population dynamics included in the model, such as metapopulation processes (Rétho et al., 2008), and that species colonization is more rapid than climate dynamics.In case of low colonization rate, the driving factor of the potential distribution is no longer the climate, but the dispersion rate and its modulation by environmental factors (Hanski and Ovaskainen, 2000).Our simulations, mainly based on climate, are then indicative and not predictive.

Environmental factors
Topography and climate effects: the pine species production shows a significant optimum for intermediate elevations (around 650 m, Fig. 6a), indicating that pine does not accommodate with low temperatures.The oak species production increases with elevation without any optimum, thus directly beneficiating of precipitation increase in mountainous areas (Fig. 6c).Several processes that are not taken into account in our simulations are able to stop the oak progression at higher elevations: competition with other species and embolism of main vessels after winter colds.Our model simulates a shift of the pine production optimum at about 950 m (Fig. 6b), due to 2.5 K temperature increase.This is particularly true if temperature during the warm part of the year is the main control of forest limit in temperate zones (Jobbagy and Jackson, 2000).Climate is changing, with frost frequencies potentially changing too.It is probable that frost will be rarer and will favour colonization at higher elevations too.With this respect, we kept as much information as possible and we did not masked distribution areas computed (Figs. 5 and 9).Finally, both species will be favored during the twenty first century and will likely colonize higher elevations than occu-pied today, but for different reasons: pine species will climb up because it does not yet reached its productivity optimum, oak species because it probably has no optimum.
CO 2 concentration: without any direct fertilization effect of atmospheric CO 2 , both species show weak changes in their distribution area: −3 and −14% for pine and oak respectively in average for the region (Figs.4b and 8b).At least is it hard in our simulations to maintain species in plains (particularly for oak species).With the direct effect of CO 2 , both species have a significant increase in productivity: +26 and +43% for pine and oak respectively (Figs. 5b,9b).This fertilizing effect even, on the basis of our data and modeling, has a stronger impact than climate change.Higher CO 2 concentration allows the tree to close its stomata, leading to a better efficiency in the water use, even if the water budget decreases, such as after year 2030.We are confident about the simulation of this forcing because CO 2 fertilizing effect has been calibrated and validated by using the observed atmospheric CO 2 concentrations for the 20th century.In addition, the direct effect of CO 2 tends to reinforce the productivity/topography relationships and in particular the oak (r 2 shifting from 0.38 to 0.54).Finally, the direct effect of CO 2 will modify deeply the distribution areas of both species.Both species will certainly maintain in plains at low elevations, mainly thanks to the CO 2 fertilization, but this is more heterogeneous in space for oak species (Figs.5b, 9b).

Species comparison (Oak -Pine)
For pine and oak species, productivity increases by about +22 and +26% in the twenty first century under direct effect of CO 2 .Both productivities will increase at a rate of about 0.35% per year (Fig. 7).Nevertheless, model simulates a climatic optimum elevation at about 1000 m for pine while such optimum is not modeled for oak.This difference is probably due to the way each species beneficiate of temperature and precipitation.These species are representative of two functional types extensively found in the Mediterranean: the oak is an evergreen sclerophyllous often found in secondary successions, and the pine is drought-adapted conifer and a colonizing species.They have a very different physiological response to drought, which lead to contrasted ecosystem functioning: P. halepensis is typically drought-avoiding water-saving, whereas Q. ilex is more tolerant to precipitation variability (Ferrio et al., 2003;Martinez-Ferri et al., 2000;Methy et al., 1997).P. halepensis develops very well in dry environment, but might be less competitive in higher altitude where precipitation is high and temperature is low.Note here that our ecophysiological model does not limit tree growth at high elevations because freeze induced embolism is an effect not taken into account.
Our results also have limitations linked to both data and models.Our measurements still do not cover optimum periods of time.The vegetation models have some lacunae which must be filled up progressively in future.Finally, our Biogeosciences, 5, 1493Biogeosciences, 5, -1504Biogeosciences, 5, , 2008 www.biogeosciences.net/5/1493/2008/vulnerability studies are based on a single scenario (IPCC-B2) of a single climatic model.They therefore do not have value of prediction but simply of indication.It will be necessary, to complete this approach, to use simulation ensembles from several climate models to deal with the climate evolution under probabilistic forms.

Conclusion
Finally, considering their present distribution areas being mainly located in the central part of PACA region and at low elevations, pine and oak species have a good chance, according to our data and modeling, to extend their habitat.This general trend is mainly (but not only) caused by the CO 2 atmospheric increase.Pine should colonize eastward and southward (along to Alps and Pyrenees slopes), while oak should colonize eastward, in the Alps.While our model does not simulate species competition, nor embolism, the pine productivity optimum with elevation should favor this species in a first stage and then leave the place for oak, in particular at higher (above 1000 m) elevations.

Fig. 1 .
Fig. 1.Indexed tree-ring series (dotted line) and bole carbon allocation simulated by MAIDEN model after MCMC calibration (plain line) (a) Regional Aleppo pine (Pinus Halepensis) dendrochronology averaged over 21 series and (b) evergreen oak series at Puechabon site.Both simulations are based on a calibration stage obtained after the convergence of five draws (5000 steps each).Arrows show excluded years (1985-1986).
. The oak MCMC calibration of MAIDEN model with the Puéchabon dendrochronological series shows a modal fit at r 2 =0.50±0.06.The agreement between simulated and observed oak chronologies was good from 1942 to 1990 and Biogeosciences, 5, 1493-1504, 2008 www.biogeosciences.net/5/1493/2008/ degraded after, probably under the influence of a biological factor such as the canopy closure.The model overestimated growth in 1996 and 1997, maybe because trees might have shed some leaves during the 1995 extremely dry year.Drought and the reconstruction of the canopy might have necessitated carbohydrates reserves through processes that are not modelled.Conversely, the model seems having overestimated drought effects in 2003 and 2005.Tree-ring width during the 1999-2005 showed a decreasing trend probably due the accumulation of stress such as increasing drought in the nineties, high temperature in 2003 and insect defoliation in 2004 and 2005.

Fig. 2 .
Fig. 2. Map of the simulated bole increments of pine (a) and of the R50 productivity index (b) in PACA southeast France (inset) for the twentieth century.Coastline is highlighted and the 21 pine stands are superimposed on the R50 productivity index map.

Table 1 .
Monthly values of the three climate parameters used (mean maximum and minimum daily temperatures in • K and monthly amount of precipitations in mm) for the common periodbetween CRU (Climatic Research Unit, University of East Anglia) and APG (ARPEGE model Météo-France) data.Corrected APG values using CRU values are noted APG/CRU.Monthly means (for temperatures) and sums (precipitations) are highlighted in bold with their annual means and standard deviations.
control period using the whole PACA grid.Mean monthly variations of the corrected climate parameters (down) between the twenty-first and twentieth centuries (noted 21/20) are presented.Annual means and standard deviations of monthly values follow.