Plan
Comptes Rendus

Ecology/Écologie
Cedrus libani (A. Rich) distribution in Lebanon: Past, present and future
[Distribution de Cedrus libani (A. Rich) au Liban : passé, présent et futur]
Comptes Rendus. Biologies, Volume 333 (2010) no. 8, pp. 622-630.

Résumé

Long-term vegetation studies are needed to better predict the impact of future climate change on vegetation structure and distribution. According to the IPCC scenario, the Mediterranean region is expected to undergo significant climatic variability over the course of this century. Cedrus libani (A. Rich), in particular, is currently distributed in limited areas in the Eastern Mediterranean region, which are expected to be affected by such climate change. In order to predict the impact of future global warming, we have used fossil pollen data and model simulations. Palaeobotanical data show that C. libani has been affected by both climate change and human activities. Populations of C. libani survived in refugial zones when climatic conditions were less favourable and its range extended during periods of more suitable climate conditions. Simulations of its future geographical distribution for the year 2100 using a dynamic vegetation model show that only three areas from Mount Lebanon may allow its survival. These results extrapolated for cedar forests for the entire Eastern Mediterranean region show that forests in Syria are also threatened by future global warming. In southern Turkey, cedar forests seem to be less threatened. These results are expected to help in the long-term conservation of cedar forests in the Near East.

Métadonnées
Reçu le :
Accepté le :
Publié le :
DOI : 10.1016/j.crvi.2010.05.003
Mots clés : Cedar, Eastern Mediterranean, Pollen, Model simulations, Conservation

Lara Hajar 1, 2 ; Louis François 3 ; Carla Khater 4 ; Ihab Jomaa 4 ; Michel Déqué 5 ; Rachid Cheddadi 2

1 CNRS-université de Toulouse, LMTG, observatoire Midi-Pyrénées, 14, avenue Édouard-Belin, 31400 Toulouse, France
2 CNRS, université Montpellier-2, institut des sciences de l’évolution CC 061, place Eugène-Bataillon, 34095 Montpellier cedex 05, France
3 Unité de modélisation du climat et des cycles biogéochimiques, université de Liège, bâtiment B5c, 17, allée du 6-Août, 4000 Liège, Belgium
4 Center for Remote Sensing, National Council for Scientific Research in Lebanon. BP 11-8281, Bir Hassan, Beirut, Lebanon
5 Météo-France, CNRS/GAME, centre national de recherches météorologiques, 42, avenue Coriolis, Toulouse, France
@article{CRBIOL_2010__333_8_622_0,
     author = {Lara Hajar and Louis Fran\c{c}ois and Carla Khater and Ihab Jomaa and Michel D\'equ\'e and Rachid Cheddadi},
     title = {\protect\emph{Cedrus libani} {(A.} {Rich)} distribution in {Lebanon:} {Past,} present and future},
     journal = {Comptes Rendus. Biologies},
     pages = {622--630},
     publisher = {Elsevier},
     volume = {333},
     number = {8},
     year = {2010},
     doi = {10.1016/j.crvi.2010.05.003},
     language = {en},
}
TY  - JOUR
AU  - Lara Hajar
AU  - Louis François
AU  - Carla Khater
AU  - Ihab Jomaa
AU  - Michel Déqué
AU  - Rachid Cheddadi
TI  - Cedrus libani (A. Rich) distribution in Lebanon: Past, present and future
JO  - Comptes Rendus. Biologies
PY  - 2010
SP  - 622
EP  - 630
VL  - 333
IS  - 8
PB  - Elsevier
DO  - 10.1016/j.crvi.2010.05.003
LA  - en
ID  - CRBIOL_2010__333_8_622_0
ER  - 
%0 Journal Article
%A Lara Hajar
%A Louis François
%A Carla Khater
%A Ihab Jomaa
%A Michel Déqué
%A Rachid Cheddadi
%T Cedrus libani (A. Rich) distribution in Lebanon: Past, present and future
%J Comptes Rendus. Biologies
%D 2010
%P 622-630
%V 333
%N 8
%I Elsevier
%R 10.1016/j.crvi.2010.05.003
%G en
%F CRBIOL_2010__333_8_622_0
Lara Hajar; Louis François; Carla Khater; Ihab Jomaa; Michel Déqué; Rachid Cheddadi. Cedrus libani (A. Rich) distribution in Lebanon: Past, present and future. Comptes Rendus. Biologies, Volume 333 (2010) no. 8, pp. 622-630. doi : 10.1016/j.crvi.2010.05.003. https://comptes-rendus.academie-sciences.fr/biologies/articles/10.1016/j.crvi.2010.05.003/

Version originale du texte intégral

Abbreviations

C. libani

Cedrus libani

GCM

General Circulation Models

gC m-2 year-1

Gram of carbon per m2 per year

IPCC

Intergovernmental Panel On Climate Change

1 Introduction

The Mediterranean region is a hot spot of biodiversity [1] that has been under continuous anthropogenic pressure for millennia [2]. Today, it faces the consequences of ongoing climate change and human pressure. During the next century, the Mediterranean is expected to undergo stronger warming and water stress which might affect all ecosystems. Most of the forests and woodlands located in the southern and Eastern Mediterranean regions are threatened by decreased moisture availability that may be related to the global warming. Among these, cedar forests constitute an emblematic ecosystem, still relatively extensive but rapidly declining.

Previous research studies, mainly in Turkey and Lebanon, have aimed to: define the bioclimatic limits of C. libani [3]; evaluate its economic value and justify its conservation status [4]; determine the cause(s) of its current regression [5]; highlight and protect new potential areas [6] and define new methods to preserve the more resistant populations [5,6]. Taberlet and Cheddadi (2002) [7] underline the importance of past vegetation dynamics in conservation studies and suggest that lessons from the past might improve the long-term efficiency of our conservation efforts. These authors [7] also outline the fact that species respond independently to climate change and not identically from one cycle to another, which requires the study of individual species dynamics. Petit et al. (2008) [8] also highlight the value of palaeoclimate data, fossils records and genealogical analyses in helping to interpret simulations of future species distribution and understanding their adaptive response to climate change.

The history of C. libani is only partially known in the Eastern Mediterranean through palynological studies performed in Lebanon, Syria and Turkey [9–17]. These studies show that the range of C. libani forests has changed since the Late Glacial.

Currently, in Lebanon, cedar forests are restricted to small mountain areas also due to increasing human activities since the beginning of the Neolithic [17,18]. Consequently, there presently remain only 12 scattered C. libani forests in Lebanon. These forests are currently shelters for high diversity and home to flora species (Convention on International Trade in Endangered Species [CITES]), and harbor a range of endemic plant species [19–21]. In Lebanon, C. libani forests require a better management strategy. Today, despite strict conservation policies, the decline of the forest is still significant and ongoing [22,23].

Model simulations [24] show that the predicted future decrease of precipitation will lead to a high probability of desertification throughout the Mediterranean region. The increases of atmospheric CO2 and temperature, and the water stress predicted by several GCM outputs [25] within the next century, may lead to an expansion of xeric vegetation in the Eastern Mediterranean region. However, the resolution of climate data used as input by Cheddadi et al. [24] does not allow a precise analysis of future changes at the scale of Lebanon. Here, we downscaled a climate dataset to simulate the future potential distribution of C. libani using a dynamic vegetation model at the scale of Lebanon.

2 Regional setting

Lebanon is situated in the Eastern Mediterranean region between Israel, Syria and the Mediterranean Sea (Fig. 1). The landscape is shaped by four parallel structures: the littoral zone, Mount Lebanon, the Bekaa valley and the Anti-Lebanon Mountains (Fig. 1). Mountain elevations range from more than 3000 m.a.s.l. for Mount Lebanon to c. 2680 m.a.s.l. for the Anti-Lebanon Mountains. These two mountain ranges are separated by the Bekaa valley (elevation c. 800 m.a.s.l.).

Fig. 1

Cedrus libani geographical distribution in Turkey [62], Syria [63] and Lebanon [45] and topography of the Eastern Mediterranean region in meters.

A Mediterranean climate prevails in Lebanon with hot and dry summers [26]. Most of the precipitation falls during autumn and winter [27]. Precipitation falls mostly on the western part of Mount Lebanon because of a foehn effect, while the major part of the Bekaa valley and the northern part of the Anti-Lebanon Mountains are arid [28].

An important herbaceous and arboreal diversity is observed in Lebanon [28]. Forests and wooded lands are mainly located on the slopes of Mount Lebanon and the southern part of the Anti-Lebanon Mountains [28]. In summary, along the littoral zone and on the western slope of Mount Lebanon, several coniferous (cedar, cypress, fir, juniper and pine), deciduous (carob tree, European hop-hornbeam, manna ash, oaks, pistachio tree), semi-deciduous (Aleppo oak) and evergreen (Palestine oak) forests and wooded lands have developed (for details, see Abi-Saleh [28] and Ministry of Agriculture and FAO [29]). The eastern slopes of the Mount Lebanon range are characterised by drier vegetation units with the development of forests and wooded lands of conifers (cypress, junipers), evergreen (Palestine oak) and semi-deciduous (Aleppo oak) populations. The Bekaa valley is presently a cultivated area producing vegetables, fruit trees and industrial crops (e.g. beetroot, cereals, potato and grains). The northern part of the Anti-Lebanon mountain range is characterised by the development of herbaceous species and junipers. Finally, on the southern part of the Anti-Lebanon Mountains, semi-deciduous (Aleppo oak), evergreen (Palestine oak) and conifers (cypress, pine) populations are found.

3 Material and methods

3.1 Palaeodata

In this article, we review the palynological data available in the Eastern Mediterranean region (Fig. 2) and synthesise them to provide an overview of past C. libani distribution (Fig. 2). All ages used in this article have been calibrated in years before present (cal. years BP) (Fig. 2), whenever possible, using CALIB 5.0 program ([30], using atmospheric data from [31]) to allow comparison between the sites. When the pollen record was poorly dated, which prevents the application of a reliable age-depth model, we decided to provide only the calibrated radiocarbon dates (i.e. the palynological sequence from Ghab valley in Syria [9,14]).

Fig. 2

Synthesis of Cedrus libani percentages from Eastern Mediterranean region main pollen records. Scales are in years cal. BP when an age-depth model is available and otherwise in cm. Pollen sequences are from Lebanon [16,17], Syria [9,14] and Turkey [11,12].

3.2 Model simulations

Simulations have been performed for both modern and future geographical distribution of C. libani using the CARAIB model. CARAIB simulates carbon cycle in the continental biosphere [32,33], biomass growth and vegetation carbon storage [34,35]. It is composed of five modules [36], which are: the hydrological cycle, photosynthesis and stomatal regulation, carbon allocation and biomass growth, heterotrophic respiration and litter/soil carbon, and biogeography. CARAIB requires climate, soil and vegetation data as input [37].

CARAIB can simulate the geographical distribution of C. libani using its modern climate boundaries. For the climate data, a local climate dataset was interpolated from meteorological stations (50 for temperature and 90 for precipitation) in Lebanon, (1930–1970) [27] with a 2’ × 2’ resolution (Fig. 3). Kriging with the R statistical software and gstat library [38] was used to spatially interpolate climate variables. Monthly mean temperatures were interpolated using universal kriging over altitude, longitude and latitude (Fig. 3a). For precipitation, ordinary kriging provided more robust interpolated data (Fig. 3b). Residuals (Table 1), which permit the evaluation of discrepancies between observed and interpolated data, are low for temperature values and slightly higher for precipitation values. Other climate data (wind speed, relative humidity, temperature amplitude, and sunshine hours) were obtained using bilinear interpolation over the same grid [39]. To obtain daily climatic data as input, CARAIB was coupled with a stochastic generator [40]. The soil texture class is taken from a data set originally at 1 × 1 resolution [41].

Fig. 3

Elevation, instrumental and interpolated climate data in Lebanon. a. Instrumental and interpolated (universal kriging) January and August temperatures in Lebanon. b. Instrumental and interpolated (ordinary kriging) January and August precipitations in Lebanon.

Table 1

Residuals (absolute dissimilarity) temperatures (°C) and precipitations (mm) between interpolated data and measurements made in meteorological stations in Lebanon.

Jan. Feb. March April May June July Aug. Sep. Oct. Nov. Dec.
Temp.
Res. min −1.58 −1.64 −1.31 −2.55 −2.38 −2.28 −2.43 −2.65 −3.02 −2.56 −1.83 −1.52
Res. max 1.77 1.38 1.49 1.26 1.91 2.08 2.44 1.86 1.93 3.46 1.57 1.48
Prec.
Res. min −73 −58 −58 −36 −20 −1.5 −0.5 −0.5 −5.5 −23 −53 −67
Res. max 88 95 69 29 16 2.4 0.5 0.4 5.7 43 50 82

The climatic boundaries of C. libani (Table 2) used for its simulation are based on the largest populations from southern Turkey and western Lebanon (Fig. 1). The small patches of C. libani found in northern Turkey and western Syria (Fig. 1) were not used in this computation because of the low accuracy of the forests map. These climatic boundaries were computed from monthly climatic data [42] and present-day distribution using GRASS GIS (http://grass.itc.it/). These values differ from those of Khouzami and Nahal (1983) [3] where C. libani climatic boundaries were computed using 30 meteorological stations located in forests in Lebanon, Syria and Turkey; these differences can be attributed to the fewer meteorological stations in Khouzami and Nahal's study (1983) [3].

Table 2

Climatic boundaries used as input for Cedrus libani simulations.

Climatic boundaries Tmin Tmaxg Wmaxg Sw Gdd
Cedrus libani −16 4 0.5 0.03 1700

The output of the CARAIB model used to predict the distribution of the species is net primary productivity (npp). The npp corresponds to photosynthesis minus autotrophic respiration and allows an accurate representation of potential distribution of a species; the presence of the species being proportional. The presence of these species is proportional to the npp values. The calibration of the model is based on the comparison of the npp values and the current geographical distribution of C. libani (Fig. 4a, b). The climate scenarios from the ARPEGE-Climat GCM [43] were interpolated on a 2’ × 2’ grid [34].

Fig. 4

a. In black: current distribution of Cedrus libani in Lebanon [45]. In gray: simulated net primary productivity above 160 gC m-2 year-1 b. Simulated present-day distribution of Cedrus libani in Lebanon using CARAIB model.

3.3 Model calibration

To calibrate the model, simulations of the present-day distribution of C. libani were performed. In order to obtain the most accurate output which best reflects the present-day distribution of C. libani as possible, we computed and tested different quantiles of current climatic boundaries of C. libani. The most suitable predicted current distribution was simulated with the 0.05 and 0.95 quantiles (Fig. 4b, Table 2). The lowest npp value for which all C. libani forests in Lebanon are predicted is around 160 gC m-2 year-1 (Fig. 4a). Thus, we considered that all values of npp equal or higher than c. 160 gC m-2 year-1 correspond to a potential presence of C. libani (Fig. 4b). This threshold is coherent with that computed by Evrendilek et al. [44]. Based on this npp threshold, the CARAIB model predicts an expanded and continuous potential forest on Mount Lebanon (Fig. 4b). However, it should be noted that human activity has not been taken into account in the vegetation model and that the forest map of Lebanon [45] has a much finer spatial resolution than the climate dataset used here (2’ × 2’). In Lebanon, C. libani trees are gathered into small patches (sometimes a stand comprising just a few trees) which does not allow satellite images of medium resolution to detect patch edges. None of the available maps in Lebanon depict all the extant occurrences of C. libani spots altogether either because of the map scale or because of the resolution of the satellite image used to derive the map.

4 Past distribution – review

Palynological records from the Eastern Mediterranean region providing percentages of C. libani are from Lebanon [16,17], Syria [9,14] and Turkey [11,12]. In order to interpret these palynological data accurately, modern pollen data obtained in/near cedar forests [16] have been used. The latter data show that C. libani pollen grains are not transported over long distances and that percentages above 5% indicate its presence near the site [16]. Using these results and palynological records from the Eastern Mediterranean region, we attempted to interpret past distributions of cedars since the Late Glacial.

Only three sequences have samples dated from the Late Pleistocene: Aammiq [16], Ghab [9,14] and Karamik Batakliği [11] (Fig. 2). The Aammiq and Karamik Batakliği records have reliable age-depth models that can be used to reconstruct past vegetation dynamics. The Ghab record has no age-depth model because the dates are not reliable [9,14].

Reconstructing C. libani expansion and development since the Late Glacial with such few pollen data is difficult. Nevertheless, the available sequences show that C. libani percentages from the Aammiq (Lebanon) and Ghab (Syria) cores decrease during the Holocene in comparison to the Bølling/Allerød period. In the Karamik Batakliği core (Turkey), we observe an opposite trend.

These results suggest that cedar forests from Lebanon and Syria were better adapted to the climate prevailing near the Bekaa and the Ghab valleys during the Bølling/Allerød than during the Holocene (Fig. 2). Probably, in order to survive the warmer climate prevailing during the Holocene, C. libani forests migrated to more adapted areas at a higher elevations on Mount Lebanon and Jabal an-Nusayriya [9,14,16], where they occur nowadays (Fig. 1).

However, during the second part of the Holocene, human interference might have affected cedar forests [18] and exaggerated the limited expansion of these populations in Lebanon and Syria. Pollen records from the Holocene in these regions do not only reflect the impact of climate change on past vegetations. In spite of human interferences, cedar forests seem to have less developed during the Holocene in the southern part of the Eastern Mediterranean region than during the Late Pleistocene.

Conversely, in Turkey, the palynological records show rather a major development of cedar forests during the Holocene (Fig. 2) which may suggest that climate was more favourable than in Lebanon and Syria.

5 Present-day distribution

Today, C. libani forests are present naturally in Lebanon, Syria and Turkey between 800 and 2100 m a.s.l. [5,6] (Fig. 1). More than 90% of the total distribution area is located in Turkey [46]. In Lebanon there are 12 scattered forests situated on the western slope of Mount Lebanon [28] (Fig. 1). The Ministry of Environment and the Ministry of Agriculture in Lebanon have designated some of these cedars forests as protected areas [47]. Thus, three forests are natural forest reserves (Barouk, Ehden and Tannourine); Al Quammoua area is designated as a protected natural site and the following forests: Barouk, Tannourine, Jaje, Quammoua, Sir El Dinniyeh and Sweisse are now protected forests. Most of the cedar forests in Turkey are located in the Taurus mountains but there are scattered populations in other parts of Anatolia [48] (Fig. 1). In Syria, as a consequence of increasing human impact, C. libani populations are currently located in an extremely restricted area in the Jabal an-Nusayriya in the northwestern part of the country [6,49–51] (Fig. 1).

6 Future distribution

Three simulations with different climate scenarios from the IPCC [25] were performed using the ARPEGE-Climat climate dataset as input to the CARAIB model. The A2 scenario projects an increase in global CO2 emissions till the year 2100. The A1B scenario predicts an increase of global CO2 values until c. 2050 then a decrease towards the year 2100. Lastly, the B1 scenario is more optimistic than the two previous ones, with a slight increase in CO2 emissions until the year c. 2040 and a decrease by the year 2100. The simulations performed using these three scenarios predict higher npp values of C. libani on the most elevated regions on Mount Lebanon (Fig. 5). However, for each scenario, major geographical differences are observed. For scenario B1, the simulated area encompasses only the Bcharre, Ehden, Jaje and Tannourine/Haddath al Jibbe cedar forests and a small area on the southern part of the Anti-Lebanon mountain range. For scenario A1B, only Ehden and Bcharre forests are simulated. For the A2 scenario, all cedar forests, except for Bcharre forest, are simulated to expand into areas where a Juniperus excelsa-dominated forest-steppe vegetation unit is currently present.

Fig. 5

Simulations of future potential areas for Cedrus libani in Lebanon using different IPCC climate scenarios.

7 Discussion

7.1 Future of Cedrus libani in Lebanon

Palaeobotanical records in the Eastern Mediterranean region show that the distribution of C. libani forests changed since the Late Glacial (Fig. 2). In Lebanon, palynological records suggest that cedars were more adapted to the Bølling/Allerød climate, which is assumed to be colder than the climate prevailing during the mid-Holocene [52,53]. The same observations have already been made for Cedrus atlantica forests in North Africa, where cedar forests were more expanded during the Late Pleistocene than at present [54]. The currently reduced development of cedar forests in the southern Mediterranean region is also attributable to human interferences such as the use of wood for construction by a succession of civilizations [18]. Consequently, cedar forests in the southern Mediterranean region are threatened not only by the climate warming during the Holocene but also by different human activities.

As the current climate in the Eastern Mediterranean region is warmer and drier than during the Bølling/Allerød [52,53], C. libani forests are becoming even more threatened nowadays in Lebanon. Consequently, the modern populations are already surviving in refugium areas. For future conservation strategies, the identification of new areas with suitable climate is crucial as C. libani may be able to survive in these locations under drier climatic conditions predicted by the various GCMs [55]. Climate during the next century is expected to be even more arid in Lebanon which will undoubtedly affect the distribution of C. libani. The available network of fossil data does not allow the evaluation of its migration rate, but a migration to higher altitudes should be expected.

Three climate scenarios allowed us to simulate areas where cedars could potentially be preserved until the year 2100. These simulations show that only a few of the existing C. libani forests will be able to survive in the northern part of Mount Lebanon (Fig. 5). The three scenarios show that cedars may be preserved near their present-day geographical distribution at higher elevations (Fig. 5). However, it is important to consider the rate of climate change during the next century and the ongoing strong human impact versus the migrational capacity of the species. Propagation of cedar seeds is accomplished by wind, gravity or animals [56]. Of these, wind is the most important agent that can carry seeds over varying distances, depending on wind velocity and plant location [56]. This natural dispersal ability allows forests to expand from their borders at several distances, depending on abiotic (i.e. wind velocity, topography, climate) and biotic influences (i.e. plant competition, grazing, deforestation) [56]. Today in Lebanon, most of C. libani forests cannot expand any higher because they are already located at the summit of mountains. Only Bcharre, Ehden and Tannourine/Haddath al Jibbe have higher altitudes. Fady et al. [57] show that cedars from Tannourine/Haddath al Jibbe forest (Fig. 4) have a high genetic diversity. As suggested by Cheddadi et al. [54] for Cedrus atlantica, by combining model simulations and genetic data, it would be appropriate to rely on populations with high genetic diversity for eventual reforestation projects. Consequently, in order to preserve C. libani forests in Lebanon efficiently, conservation efforts should focus on protecting forests where cedars present high genetic diversity, such as the Tannourine/Haddath al Jibbe forest and where high elevations are available.

Importantly, the ARPEGE-Climat climate dataset has a resolution between 0.5° (30’) and 1° (60’) in the Eastern Mediterranean region. The downscaling of such GCM resolution will be necessary in order to produce more accurate model simulations that portray the potential distribution at local scales in mountainous areas.

7.2 Future of Cedrus libani in the Eastern Mediterranean region

Palynological studies show that C. libani dynamics depend upon a range of natural and human-induced forcing mechanisms in the entire Eastern Mediterranean region. Our simulations allowed us to identify suitable areas for C. libani. In Turkey, most cedar populations have not yet reached the highest elevation in the Taurus mountain and are currently well expanded (Fig. 1). In these areas, there exist higher elevations where cedars could migrate and adapt. In the northern part of Turkey, there are small remnant populations [58] which are under threat because of the low mountain elevations. In Syria, scattered C. libani populations have also already reached the highest elevations (Fig. 1). Cedar populations in Turkey and Syria have high genetic diversity [57,59]. Therefore, there is less threat of loss of genetic diversity if these populations are preserved. However, it should be noted that Bou Dagher-Kharrat et al. (2007) [59] reported that genetic erosion has occurred in some Cedrus populations in the Eastern Mediterranean. Thus, maintaining population dynamics and regeneration of C. libani populations should be a priority matter for conservation managers [60]. Our simulations allow the identification of areas where cedars could potentially be preserved. Future studies need to take into account migrational processes. A network of well-dated fossil data may help to estimate these migrational rates.

7.3 Implications for conservation strategies

Previous studies [7,8,60,61] have shown that increased knowledge and understanding of long-term plant dynamics studies is needed for better evaluation of the predicted impacts of future global warming. The integration of palaeoecology, model simulations and genetics is necessary for conservation issues [7,8]. However, very few studies dealing with conservation do include long-term palaeoecological data [4–6,23,58], thus many conservation strategies are ill-informed regarding long-term trajectories of individual species. The results of this study aim at evaluating the impact of future global warming on C. libani and proposing where long-term conservation efforts should be applied:

  • 1. The threats on Mediterranean species should be re-evaluated at different temporal scales, taking into account the long-term effect of future climate change based on palaeobotanical data and model simulations.
  • 2. In order to preserve a species or a population, we need to identify and protect potential distribution areas and migrational corridors through a long-term conservation policy.

Potential migration areas for the survival of C. libani populations from Bcharre, Ehden and Tannourine/Haddah-al-Jibbé have been identified but these are not currently protected areas. If human activities (such as grazing) are allowed to continue over the future decades, then this will be a serious threat to the preservation of C. libani populations.

8 Conclusion

Studies on the long-term dynamics of C. libani are necessary, as the current distribution of this species is restricted to mountainous areas. Model simulations show that C. libani could be preserved around its current distribution area or/and on more elevated areas over Mount Lebanon. Such prospective is coherent with the simulated velocity of temperature change at the global scale [65]. In south Turkey, potential areas have also been identified. In Syria, the only remaining C. libani population has already reached the highest elevations. More importantly, the conservation of C. libani will depend on its ability to migrate across a predominantly cultivated landscape. The available palaeobotanical data are not sufficient to evaluate its migration rates. In order to better interpret our model simulations, we have investigated the long-term dynamics of C. libani using palaeobotanical data.

Our study shows that the threats on Mediterranean species should be re-evaluated considering the long-term effect of future global warming and that highlighting future potential areas for the survival of C. libani (or others species) is necessary to protect these areas against future human activities.

Conflict of interest statement

No conflict of interest.

Acknowledgements

The authors thank Simon Brewer for his contribution and interesting discussions about geostatistical analysis and Serge Muller for interesting discussions, which helped to improve the manuscript. This work is supported by a Ph.D. grant from the National Council of Scientific Research of Lebanon (NCSR) and from the Belgium National Foundation for Scientific Research (F.R.S.-FNRS; Crédit aux chercheurs) as well as by a research program of the NCSR and the ECOCHANGE European program.

This is an LMTG contribution and ISEM contribution 2010-048.


Bibliographie

[1] N. Myers; R.A. Mittermeier; C.G. Mittermeier; G.A.B. Da Fonseca; J. Kent Biodiversity hot spots for conservation priorities, Nature, Volume 403 (2000), pp. 853-858

[2] J.-L. Vernet Man and vegetation in the Mediterranean area during the last 20,000 years (F di Castri; A.J. Hansen; M. Debussche, eds.), Biological Invasions in Europe and the Mediterranean Basin, Kluwer Academic Publishers, Dordrecht, Netherlands, 1990

[3] M. Khouzami; I. Nahal Les bioclimats du Cèdre du Liban (Cedrus libani A. Rich.) et leurs particularités dans son aire naturelle, Agricultural Sciences, Volume 5 (1983), pp. 39-62

[4] E.J. Sattout; S.N. Talhouk; P.D.S. Caligari Economic value of cedar relics in Lebanon: an application contingent valuation method for conservation, Ecological economics, Volume 61 (2007), pp. 315-322

[5] M. Boydak Regeneration of Lebanon cedar (Cedrus libani A. Rich.) on karstic lands in Turkey, Forest Ecology and Management, Volume 178 (2003), pp. 231-243

[6] S. Khuri; M.R. Shmoury; R. Baalbaki; M. Maunder; S.N. Talhouk Conservation of the Cedrus libani populations in Lebanon: history, current status and experimental application of somatic embryogenesis, Biodiversity and Conservation, Volume 9 (2000), pp. 1261-1273

[7] P. Taberlet; R. Cheddadi Quaternary Refugia and Persistence of Biodiversity, Science, Volume 297 (2002), pp. 2009-2010

[8] R.J. Petit; F.S. Hu; C.W. Dick Forests of the past: a window to future changes, Science, Volume 320 (2008), pp. 1450-1452

[9] J. Niklewski; W. Van Zeist A Late Quaternary pollen diagram from NW Syria, Acta Botanica Neerlandica, Volume 9 (1970), pp. 737-754

[10] S. Bottema A pollen diagram from the Syrian anti-Lebanon, Paléorient, Volume 3 ( 1975–1977 ), pp. 259-268

[11] W. van Zeist; H. Woldring; D. Stapert Late Quaternary vegetation and climate of South-western Turkey, Palaeohistoria, Volume 17 (1975), pp. 53-145

[12] W.J. Eastwood; N. Roberts; H.F. Lamb; J.C. Tibby Holocene environmental change in the southwest Turkey: a palaeoecological record of lake and catchment-related changes, Quaternary Science Reviews, Volume 18 (1999), pp. 671-695

[13] M. Vermoere; E. Smets Late Holocene Environmental change and the Record of Human Impact at Gravgaz near Sagalassos, southwest Turkey, Journal of Archaeological Science, Volume 27 (2000), pp. 571-595

[14] Y. Yasuda; H. Kitagawa; T. Nakagawa The earliest record of major anthropogenic deforestation in the Ghab Valley, northwest Syria: a palynological study, Quaternary International, Volume 73/74 (2000), pp. 127-136

[15] M. Vermoere; T.V. Thuyne; S. Six; L. Vanhecke; M. Waelkens; E. Paulissen; E. Smets Late Holocene local vegetation dynamics in the marsh of Gravgaz (southwest Turkey), Journal of Paleolimnology, Volume 27 (2002), pp. 429-451

[16] L. Hajar; C. Khater; R. Cheddadi Vegetation changes during the late Pleistocene and Holocene in Lebanon: a pollen record from the Bekaa Valley, The Holocene, Volume 18 (2008), pp. 1089-1099

[17] L. Hajar; M. Haïdar-Boustani; C. Khater; R. Cheddadi Environmental changes in Lebanon during the Holocene: man vs climate impacts, Journal of Arid Environments, Volume 74 (2010) no. 7, pp. 746-755

[18] H.C. Loffet Sur quelques espèces d’arbres de la zone syro-palestinienne et libanaise exportées vers l’Égypte Pharaonique, Archaeology and History in Lebanon, Volume 9 (2004), pp. 10-33

[19] M. Semaan M, R. Haber. Scientific Study Report on the Plant Diversity of Ehden Forest Natural reserve, Ministry of environment, Beirut, Lebanon, 1998.

[20] G Tohme, N.Z Jeze, H. Tohme, S.H. Bloquet, M. Karakireh, R. Jeze. Horsh Ehden Nature Reserve, National Council for Scientific Research and Ministry of Environment, Beirut, 1999.

[21] G Tohme, s.N. Talhouk, H. Tohme, S.H. Bloquet, M. Karakireh, R. Jeze. Arz Al-shouf Nature Reserve, National Council for Scientific Research and Ministry of Environment, Beirut, 1999.

[22] Jomaa I. Analyse diachronique de la fragmentation des forêts du Liban. Thesis, Université Paul Sabatier, Toulouse, 2008.

[23] E.J. Sattout; P.D.S. Caligari; S.N. Talhouk Perspectives for sustainable management of cedars forests in Lebanon: situation analysis and guidelines, Environment, Development and Sustainability, Volume 10 (2008), pp. 107-127

[24] R. Cheddadi; J. Guiot; D. Jolly The Mediterranean vegetation: what if the atmospheric CO2 increased?, Landscape Ecology, Volume 16 (2001), pp. 667-675

[25] IPCC. Climate change 2007: the AR4 synthesis report, 2007.

[26] M.C. Peel; B.L. Finlayson; T.A. McMahon Updated world map of the Köppen-Geiger climate classification, Hydrology and Earth System Sciences, Volume 11 (2007), pp. 1633-1644

[27] Service météorologique du Liban, Atlas climatique du Liban/République libanaise, Direction de l’aviation civile- Beyrouth: impr. Catholique, Tome. I (pluie, température, pression, nébulosité), 1977.

[28] Abi-Saleh B. Étude phytosociologique, phytodynamique et écologique des peuplements sylvatiques du Liban. Thesis, Fac. des Sciences et Techniques Saint-Jérôme, 1978.

[29] Ministry of Agriculture and FAO, Forest Map of Lebanon, 2005.

[30] M. Stuiver; P.J. Reimer Extended 14C database and revised CALIB radiocarbon calibration program, Radiocarbon, Volume 35 (1993), pp. 215-230

[31] P.J. Reimer; M.G.L. Baillie; E. Bard; A. Bayliss; J.W. Beck; C.J.H. Bertrand; P.G. Blackwell; C.E. Buck; G.S. Burr; K.B. Cutler; P.E. Damon; R.L. Edwards; R.G. Fairbanks; M. Friedrich; T.P. Guilderson; A.G. Hogg; K.A. Hughen; B. Kromer; G. McCormac; S. Manning; C.B. Ramsey; R.W. Reimer; S. Remmele; J.R. Southon; M. Stuiver; S. Talamo; F.W. Taylor; J. van der Plicht; C.E. Weyhenmeyer Intcal04 terrestrial radiocarbon age calibration 0-26 cal kyr BP, Radiocarbon, Volume 46 (2004), pp. 1029-1058

[32] P. Warnant; L. François; S. Strivay; J.-C. Gérard CARAIB: A global model of terrestrial biological productivity, Global Biogeochemical Cycles, Volume 8 (1994), pp. 255-270

[33] B. Nemry; L.M. François; P. Warnant; F. Robinet; J.-C. Gérard The seasonality of the CO2 exchange between the atmosphere and the land biosphere: a study with a global mechanistic model, Journal of Geophysical Research, Volume 101 (1996), pp. 7111-7125

[34] L.M. François; C. Delire; P. Warnant; G. Munhoven Modelling the glacial-interglacial changes in the continental biosphere, Global and Planetary Change, Volume 16–17 (1998), pp. 37-52

[35] L.M. François; Y. Goddéris; P. Warnant; G. Ramstein; N. de Noblet; S. Lorenz Carbon stocks and isotopic budgets of the terrestrial biosphere at mid-Holocene and last glacial maximum times, Chemical Geology, Volume 159 (1999), pp. 163-189

[36] L. François; M. Ghislain; D. Otto; A. Micheels Late Miocene vegetation reconstruction with the CARAIB model, Palaeogeography, Palaeoclimatology, Palaeoecology, Volume 238 (2006), pp. 302-320

[37] D. Otto; D. Rasse; J. Kaplan; P. Warnant; L. François Biospheric carbon stocks reconstructed at the Last Glacial Maximum: comparison between general circulation models using prescribed and computed sea surface temperatures, Global and Planetary Change, Volume 33 (2002), pp. 117-138

[38] E.J. Pebesma Multivariable geostatistics in S: the gstat package, Computers and Geosciences, Volume 30 (2004), pp. 683-691

[39] H. Akima A method of bivariate interpolation and smooth surface fitting for irregularly distributed data points, ACM Transactions on Mathematical Software, Volume 4 (1978), pp. 148-159

[40] B. Hubert; L. François; P. Warnant; D. Strivay Stochastic generation of meteorological variables and effects on global models of water and carbon cycles in vegetation and soils, Journal of Hydrology, Volume 212–213 (1998), pp. 318-334

[41] L. Zobler. A world soil file for global climate modelling, NASA Technical memorandum, Goddard Institute for Space Studies (GISS), New York, 1986.

[42] M. New; D. Lister; M. Hulme; I. Makin A hight-resolution data set of surface climate over global land areas, Climate Research, Volume 21 (2002), pp. 1-25

[43] M. Déqué; C. Dreveton; A. Braun; D. Cariolle The ARPEGE/IFS atmosphere model: a contribution to the French community climate modelling, Climate Dynamics, Volume 10 (1994), pp. 249-266

[44] F. Evrendilek; S. Berberoglu; S. Taskinsu-Meydan; E. Yilmaz Quantifying carbon budgets of conifer Mediterranean forest ecosystems, Turkey, Environmental Monitoring and Assessment, Volume 119 (2006), pp. 527-543

[45] R. Baltaxe. Projet de bonification intégrale de la montagne libanaise. Report on mapping the forests of Lebanon. Lebanon: FAO and Green Plan, 1965.

[46] M. Bariteau. Ferrandes Les Cèdres in: Gallais A, Bannerot H (Eds.), L’amélioration des plantes. INRA Paris, 1992. pp. 732–43.

[47] Association for Forests, Development and Conservation, State of Lebanon's Forests. Georges Mitri, 2007

[48] M. Sevim Lübnan sedirinin Türkiye’deki tabii yayılışı ve ekolojik şartları (Die natürliche verbreitung und standortsbedingungen der Libanonzeder in der Türkei), Orman Umum Müdürlüğü Yayını (Ankara) 143/24, 1955

[49] FAO, Databook on endangered tree and shrub species and provenances, in: FAO, Forestry paper 77, Food and Agriculture Organization of the United Nations, Rome, 1986.

[50] M. Khouzami. The Lebanese cedar forests, in: Proceedings of the First National Conference on the Cedar of Lebanon, Present and Future. American University of Beirut, Lebanon, 1994.

[51] S. Khuri, S.N. Talhouk, Cedar of Lebanon (Cedrus libani A.Rich.), in: A. Farjon, C.N. Page (compilers) Conifers. Status Survey and Conservation Action Plan. Gland, Switzerland and Cambridge, UK: IUCN/SSC Conifer Specialist Group, IUCN, 1999, pp. 108–11.

[52] M. Bar-Mathews; A. Ayalon Late Quaternary paleoclimate in the Eastern Mediterranean region from stable isotope analysis of speleothems at Soreq Cave, Israel, Quaternary Research, Volume 47 (1997), pp. 155-168

[53] S.A. Robinson; S. Black; B.W. Sellwood; P.J. Valdes A review of palaeoclimates and palaeoenvironments in the Levant and Eastern Mediterranean from 25,000 to 5000 years BP: setting the environmental background for the evolution of human civilisation, Quaternary Science Reviews, Volume 25 (2006), pp. 1517-1541

[54] R. Cheddadi; B. Fady; L. François; L. Hajar; J.P. Suc; K. Huang; M. Demarteau Glacial refugia and genetic diversity of Cedrus atlantica, Journal of Biogeography, Volume 36/7 (2009), pp. 1361-1371

[55] F. Médail; K. Diadema Glacial refugia influence plant diversity patterns in the Mediterranean Basin, Journal of Biogeography, Volume 13 (2009), pp. 1333-1345

[56] A.J. Fordham; L. Spraker Propagation manual of selected Gymnosperms, Arnoldia, Volume 37 (1977), pp. 1-88

[57] B. Fady; F. Lefèvre; G.G. Vendramin; A. Ambert; C. Régnier; M. Bariteau Genetic consequences of past climate and human impact on Eastern Mediterranean Cedrus libani forests. Implications for their conservation, Conservation Genetics, Volume 9 (2008), pp. 85-95

[58] M. Boydak, Ü. Eler, N. Pehlivan. Antalya-Elmalı Yöresi Sedirlerinin (Cedrus libani A. Rich.) Gençlestirilmesinde Denetimli Yakma ve Diger Bazı Faktörlerin Basarı Üzerine Etkileri, Batı Akdeniz Ormancılık Arastırma Müdürlügü, Teknik Rapor 2, 1996.

[59] M. Bou Dagher-Kharrat; S. Mariette; F. Lefèvre; B. Fady; G. Grenier-de March; C. Plomion; A. Savouré Geographical diversity and genetic relationships among Cedrus species estimated by AFLP, Tree Genetics and Genomes, Volume 3 (2007), pp. 275-285

[60] H.R. Delcourt; P.A. Delcourt; T. Webb Dynamic plant ecology: the spectrum of vegetational change in space and time, Quaternary Science Reviews, Volume 1 (1983), pp. 153-175

[61] J.C. Ritchie Current trends in studies of long-term plant community dynamics, New Phytologist, Volume 130 (1995), pp. 469-494

[62] I. Atalay General ecological properties of natural occurrence areas of cedar (Cedrus libani A. Rich) forest regioning of seed transfer of cedar in Turkey, Orman Genel Müdürlüğü yayini (Ankara), Volume 663 (1987)

[63] P. Quézel, F. Médail. Écologie et biogéographie des forêts du bassin méditerranéen, Paris, Elsevier, 2003.

[64] J.-M. Laurent; L. François; A. Bar-Hen; L. Bel; R. Cheddadi European bioclimatic affinity groups: data-model comparisons, Global and Planetary Change, Volume 61 (2008), pp. 28-40

[65] S.R. Loarie; P.B. Duffy; H. Hamilton; G.P. Asner; C.B. Field; D.D. Ackerly The velocity of climate change, Nature, Volume 462 (2009), pp. 1052-1055


Commentaires - Politique