Next Article in Journal
Phylogenetics and Biogeography of Lilium ledebourii from the Hyrcanian Forest
Next Article in Special Issue
Small Mammal Diversity in Response to Land Transformation and Seasonal Variation in South Africa
Previous Article in Journal / Special Issue
Pronounced Seasonal Diet Diversity Expansion of Golden Eagles (Aquila chrysaetos) in Northern Greece during the Non-Breeding Season: The Role of Tortoises
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Biological Interaction as a Possible Ultimate Driver in the Local Extinction of Cedrus atlantica in the Iberian Peninsula

1
Department of Botany, University of Granada, C.U. Fuentenueva, 18071 Granada, Spain
2
Department of Botany, Ecology and Plant Physiology, University of Cordoba, C.U. Rabanales, 14014 Cordoba, Spain
*
Author to whom correspondence should be addressed.
Submission received: 16 December 2021 / Revised: 9 February 2022 / Accepted: 12 February 2022 / Published: 15 February 2022
(This article belongs to the Special Issue 2021 Feature Papers by Diversity’s Editorial Board Members)

Abstract

:
The presence of Cedrus atlantica on the European continent, including, especially, the determination of the time of its disappearance from the Iberian Peninsula, is one of the most controversial issues in recent decades regarding the successive extinction of conifers in the Western Mediterranean. This work propounds the possibility that C. atlantica and Pinus nigra could have co-habited in the past, mutually excluding each other in the areas with suitable conditions for both species, where, ultimately, the one that was the most competitive would have remained. The niche overlap in the two-dimensional ecological space was analyzed. In addition, the potential distribution of both species in the Western Mediterranean today and two past periods (Last Glacial Maximum and Mid-Holocene) was modeled to identify their common geographic area of distribution. The species showed very well differentiated niches and a distribution of their habitats virtually segregated by continents since the Mid-Holocene (P. nigra in Europe and C. atlantica in Africa), which responds to differences in climatic affinities. However, the contact of the bordering areas of their distributions in the Baetic mountain range suggests that C. atlantica could have maintained its presence in the Iberian Peninsula until recent times. P. nigra would have displace it in later stages due to its greater prevalence on the continent, so it would have had greater opportunities to occupy the available space.

1. Introduction

The presence of Cedrus atlantica (Endl.) Carrière on the European continent and the time of its disappearance from the Iberian Peninsula are issues that have raised controversy in recent decades regarding the successive extinction of conifers in the Western Mediterranean [1,2,3]. The presence of Cedrus in the Iberian Peninsula has been referenced from the Eocene to the Holocene [4] and widely in Europe since the Miocene [5]. However, the only fossil macro-remain that confirms its autochthonous origin corresponds to leaves of C. blombley found on the island of Rhodes dated to the Middle Pleistocene (500 ka) [6]. The absence of macro-remains (charred wood) of Cedrus on the European continent led to the statement that the presence of cedar pollen in the Iberian fossil record responded to its airborne uptake by winds from Africa [2,5,7,8]. However, the constancy of the pollen type in numerous deposits [1,3,7,9,10,11,12], its proportion and the low dispersal power of grains [1,13,14]; the Cedrus–Abies assemblage in the records [10]; and the low correlation between the presence of pollen and the origin of the winds from Africa [15,16] question this possibility and support the hypothesis that Cedrus lived in southern Europe during the Holocene as a component of mixed conifers woods or as small copses [3].
The distribution dynamics of C. atlantica in Europe and Africa have been conditioned by climate, especially by the warmer temperatures in winter and the dryness of summer [17,18]. Thus, the warmest periods of the Quaternary, such as the last interglacial (200–140 ka) or the Holocene (11 ka-present), determined the migration of its populations towards the south and towards higher altitudes, in a progressive reduction in its distribution [17]. In the Iberian Peninsula, its last populations had to take refuge in the southern mountains of the Baetics as small, isolated nuclei [3,15]. In this context, it has been proposed that the definitive extinction of Cedrus in Europe during the Holocene could be motivated by human action [16,19,20], or even by competition with other species [3,16]. In this sense, an examination of the pollen diagrams of the available deposits on both sides of the Strait of Gibraltar (Baetic Mountain Range and Rif) revealed inverse patterns in the relative abundance of the Pinus and Cedrus pollen types; while, in the Baetics, Pinus dominates and Cedrus is scarce, in the Rif, Pinus appears meager where Cedrus is the dominant type [1,16,17,21,22,23,24,25,26]. If the ecological affinities of the species are not mutually exclusive and in the absence of barriers to dispersal, this negative association between species suggests the possibility that species of both genera have interacted in the past, mutually excluding each other in areas where they shared a niche and where the one that was more competitive would have ultimately remained [27,28].
The objective of this study is to identify the interspecific interaction through a niche overlap analysis [29,30]. Pinus nigra Arnold was selected among the mountain conifers distributed in the Western Mediterranean due to its ecological affinity with C. atlantica [31,32,33,34]. Likewise, the distribution of both species in the Iberian Peninsula was explored and the location of the possible common area was identified, both currently and in two moments in the past, Last Glacial Maximum (LGM; 21 ka) and Holocene Climate Optimum (HCO; 6 ka). The variations in the distribution of both species would make it possible to understand their respective historical geographical dynamics.
The use of Species Distribution Models (SDMs) has been a complement to the traditional techniques used in Paleobiogeography that can help identify knowledge gaps derived from inherent geographic, temporal and taxonomic biases [35,36,37,38,39]. In fact, SDMs have been previously applied in the projections of different climatic conditions in the past with C. atlantica [14,40,41] and P. nigra [31,42,43].

2. Materials and Methods

To identify overlap in the potential geographic distribution of C. atlantica and P. nigra, Species Distribution Models (SDM) were used. These models were projected into the past (LGM, 21 ka; HCO, 6 ka). On the other hand, the relationship of the niches in the two-dimensional ecological space was revealed by a principal coordinate analysis (PCA), from which the equivalence and similarity of the niches of both species was evaluated.

2.1. Study Area

The study area was framed in the Western Mediterranean, where C. atlantica may have reached its maximum extension during the last glaciation (Würm, 21 ka) [5]. It covers the Iberian and Italic peninsulas, the Western Alps and the Dinarides, in Europe; the Atlas Mountains and Rif, in Africa; and Corsica, Sardinia, Sicily and the Balearic Islands.
For the purposes of the different analyses, the study area was defined by the coordinates lat 29° N–46° N; lon 12° W–19° E. For the niche overlap, we defined the surrounding environment of each species as the geographic space delimited by the extreme coordinates of the observations used (Figure 1).

2.2. Species

Cedrus atlantica is currently distributed in mountainous areas of the Rif and Middle Atlas in Morocco and the Tell Atlas and Aures Mountains in Algeria [44], between 1300 and 2600 m of altitude, where annual rainfall ranges between 500 and 2000 mm and the minimum temperatures of the coldest month range from −1 to −8 °C [45]. For the analysis of this work, observations made directly in the field, in Morocco, were used; while in Algeria, the recognition of the species was conducted by photo interpretation (Figure 1).
Pinus nigra is a species of medium and high mountain (1500–2500 m of altitude) of circum-Mediterranean distribution, which extends through southern Europe, from the Iberian Peninsula to Anatolia, and reaches the Caucasus, with occasional presence in North Africa, in the Rif and Tell Atlas. The subspecies P. nigra subsp. salzmannii (Dunal) Franco, the one contemplated in this study, is limited to the Western Mediterranean (Cévennes, Pyrenees, Iberian and Central Systems and Baetic mountain range in Europe; Rif and Tell Atlas in Africa). It inhabits rocky terrain where it presents greater resistance to drought than other species of medium- and high-mountain pine [18,32,46]. For the analysis of this work, 257 localities were randomly selected, at the working resolution (30 arc seconds), from the database of the Global Biodiversity Information Facility (GBIF; www.gbif.net). Redundant observations at working resolution were eliminated (Figure 1).

2.3. Variables

The niche predictors were selected from the 19 bioclimatic variables of the 1960–1990 series for the present, from the WorldClim 1.4 portal, (www.worldclim.com, accessed on 1 February 2013). Preliminarily, a survey of the variables that best fit the distribution models of C. atlantica was carried out with MaxEnt [47], implemented in the ‘dismo’ R package [48]. Those variables with the highest percentage of contribution and importance of permutation were selected alternatively when the model was adjusted [49], discarding the rest of the highly correlated variables (values higher than 0.7) in each choice (Table A1).
Finally, the following four predictor variables were selected:
  • Minimum temperature of the coldest month (bio6);
  • Temperature annual range (bio7);
  • Mean temperature of the wettest quarter (bio8);
  • Precipitation of the driest month (bio14).
To project the niche model in the geographical space, the climatic layers of WorldClim 1.4 were used, at a resolution of 30 arc seconds for current conditions and HCO (6 ka) and at a resolution of 2.5 arc minutes for LGM (21 ka). MIROC-ESM was chosen among the possible general circulation models [50], because it presented the best fit in its evaluation for the two moments in time in the past for C. atlantica (Table A2). In this evaluation, the available evidence of Cedrus atlantica in the Western Mediterranean from the fossil record for the two periods referred was used. In this review, the sites of the European Pollen Database [51] were considered, as well as compilations of works in the Iberian Peninsula [22] and in North Africa [14,17,21,23,24,25,26].

2.4. Analysis

The niche of each species was modeled and projected in the geographical space using the MaxEnt algorithm, included in the ‘dismo’ R’ package [48], a machine learning-type algorithm that only requires presence records and a reference environment (background) to calibrate the models [52,53,54]. The potential habitat of each species was expressed through the binary transformation of suitability (suitable/unsuitable habitat for the species) applying, as threshold, the maximum of the sum (sensitivity + specificity). The overlap of the potential habitat of both species indicates the favorable areas for both species, thus where interaction between them is possible. There are numerous antecedents in the application of the SDM to the past and, especially, in the Iberian Peninsula [31,35,55].
The validation of the models was carried out by contrasting the value of the observed area under the curve (AUC) with the null model AUC [56]. The null model was built with the distribution of AUC values from 100 SDM iterations built with n randomly selected points (257 for P. nigra and 211 for C. atlantica) from the background. The null model AUC represents the 0.95 quantile in this distribution and the model is significantly valid (p < 0.05) if the observed AUC value is greater than the null model AUC.
On the other hand, the niche overlap in the two-dimensional environmental space was evaluated using the ‘ecospat’ R’ package [57]. This is defined by the two main axes of a PCA built from the values of the four variables selected in the surrounding environment of C. atlantica and P. nigra. (Figure 2). This two-dimensional environmental space was gridded (100 × 100 cells) and the density surfaces for the observations were projected on this grid to represent the niches [29]. The relative positions of the niches of both species allowed us to infer possible interspecific relationships.
Niche overlap was measured with Schoener’s D index (ranging from 0, no overlap, to 1, total overlap). The niche overlap was contrasted using the equivalence and similarity tests [29,30]. In each of them, the observed D-value is compared with a null distribution of D-values obtained as the result of the iteration of 100 superposition operations between pairs of sets of random observations.
The equivalence test is aimed at identifying identical niches in the whole surrounding area. In this test, each of the two sets of observations of the null distribution is built by random selection among the whole surrounding area of the two species. The niches of the two species are admitted to be equivalent if the observed overlap value is significantly higher than the 0.95 quantile of the null distribution. That is, the two niches are more similar to each other than any two random sets of the whole surrounding area.
The similarity test assesses the similarity of one niche to another in relation to the surrounding environment in which it is found. The surrounding environment of each species is delimited by the extreme coordinates of the observations (Figure 1). In this test, the total of observations from which the pairs are drawn at random, for the null distribution, was constructed with the observations of one species and an extract from the surrounding environment of the opposite species (in a number equal to that of their observations). Thus, in this test, two asymmetric analyzes are carried out in which the observations of each species are compared with the surrounding environment of the species being confronted. Therefore, two species with environmental affinities would have higher overlap values than expected when comparing the niche of one species with the surrounding environment of the opposite. The niches of the two species are significantly similar if the D-value exceeds the 0.95 quantile of the null distribution [29,30].

3. Results

3.1. Geographical Distribution

The four selected environmental variables defined the fundamental niche of both species with a high level of prediction of the optimal habitat for the current conditions. The validation of the models showed that the observed AUC values greatly exceeded those marked for a null distribution [56], with a significance of 95% (Table 1; Figure A1).
According to the projections made in MaxEnt, the potential habitats of C. atlantica and P. nigra are mostly geographically segregated from each other at present (Figure 3c). Currently, the optimal habitat for C. atlantica is distributed mainly in North Africa, from the High Atlas in the immediate vicinity of the Atlantic Ocean to the Tell Atlas and Aures Mountains, via the Rif and the Middle Atlas. In Europe, this species finds its optimal habitat in the Iberian Peninsula, Baetic Mountain Range, and in the Central System, extending to the Leon Mountains. The predicted distribution of P. nigra is framed in the European Mediterranean basin, from the Iberian and Central Systems and the Cantabrian Mountains, in the Iberian Peninsula; to the Apennines and western Alps, through the Pyrenees and Massif Central. Marginally, it is also represented in the Aures Mountains, Rif and Middle Atlas on the African continent. At present, C. atlantica and P. nigra would share their distribution in the Baetic Mountain Range and Central System in Europe and in the Middle Atlas and Aures Mountains in Africa (Figure 3c).
Compared with current conditions, the projection to the LGM (21 ka) shows the descent in altitude of the most favorable conditions for C. atlantica, where it would share the habitat with P. nigra (Figure 3a). This was revealed by a greater extension of the potential habitat of the species throughout the mountain systems of Africa and the prevalence of the plateaus in the Iberian Peninsula, as well as their presence in the mountain systems near the Tyrrhenian. In the Iberian Peninsula, the common areas would extend through the Iberian System and the inland areas close to the Cantabrian Mountains. In Africa, the common habitat would be linked to the High, Middle and Saharan Atlas and the Aures Mountains (Figure 3a).
The warmer conditions of the HCO (6 ka) with respect to the LGM would have determined the ascent in altitude of the optimal conditions for both species (Figure 3b). This would have resulted in a drastic reduction in the potential area of C. atlantica and its isolation in a number of disjoint areas, largely coinciding with the current distribution of the species. The presence of the habitat of each of the species would have been virtually restricted to one continent, C. atlantica to Africa and P. nigra to Europe, with negligible presence on the opposite continent, where both species would have shared limits of habitat distribution marginally.

3.2. Niche Overlap

The analysis of niche overlap in the two-dimensional ecological space reflects the overlap with a D-value close to 0.37. According to the strict nature of the test, the equivalence test result was that both species niches were non-significantly equivalent (they did not match completely). Contrary to the equivalence test, both similarity tests indicated that both species niches were significantly similar to each other considering their different surrounding environment (Table 2; Figure A2).
The arrangement of the density grids in the two-dimensional ecological space shows a relative displacement between the niches of both species along the second axis of the PCA. (Figure 4). This could be interpreted as a coincidence in the affinity of the species for the variables that determine the first axis of the PCA (precipitation of the driest month and mean temperature of the wettest quarter; Figure 2), while differences appear in the niche determined by the two remaining variables, minimum temperature of the coldest month, whose lowest values are better tolerated by C. atlantica (from −6 to −2 °C, compared to −3–2 °C, where P. nigra lives); and the temperature annual range, with C. atlantica better tolerating higher ranges in temperature (from 29 to 36 °C, compared to 23–26 °C in which P. nigra lives; Figure A3 and Figure A4).

4. Discussion

According to our results, C. atlantica currently extends the distribution of its potential habitat in the Iberian Peninsula even though the species does not currently inhabit Europe. These results are consistent with other models made for C. atlantica [17,41]. During the LGM (21 ka), the extension of the habitat would occupy most of the center of the Iberian Peninsula, dominating the plateaus and mountainous areas of the Iberian System, the interior areas near the Cantabrian Mountains and most of the Baetic Mountain Range. The presence of Cedrus in the fossil records could support these projections to the past [11,58,59]. The increase in temperature during the HCO (6 ka) determined the increase in altitude of both species. Populations of C. atlantica were isolated in the mountain systems, which served as a refuge, from where the species spread to the present day. This is consistent with the dynamics of cedars in North Africa during the Holocene found in fossil records [1,15,17,21,23,24,26]. Therefore, the dynamics of the cedars in Europe have been linked to the climate since the last ice age, so their disappearance on the continent could have been due to a greater extent to a progressive reduction in their potential habitat [17,18].
This would not be incompatible with the existence of biotic relationships as a determining agent of local extinction in the redoubts present in the Baetic Mountain Range (Postigo-Mijarra et al., 2010). Bearing the possibility that C. atlantica and P. nigra interacted in the Iberian Peninsula in the past, occupying similar environmental spaces, they could mutually have excluded each other in mountainous areas where the more competitive species would have remained. The starting hypothesis arises from the fact that the distribution of each species is restricted to separate continents and that the Cedrus and Pinus spp. show some negative temporal association in the pollen profiles of the fossil records from both sides of the Mediterranean [1,16,17,21,22,23,24,25,26]. However, both species show a similar ecological behavior (they form forests in the high Mediterranean mountains), which raises the question of whether and to what extent they share a niche. If so, the negative association of their presence (geographical and temporal) could be associated with biological interaction [27,28].
The results obtained indicate that the distribution of the potential habitat of C. atlantica and P. nigra overlapped both in Africa and in Europe during the Pleistocene, being able to maintain contact between their populations in their distribution limits. In addition, the niche overlap tests indicate that both species habitats are significantly similar to each other. In the geographical space, the results obtained indicate that the distribution of the potential habitat of C. atlantica and P. nigra overlapped both in Africa and in Europe during the Pleistocene, being able to maintain contact between their populations in their distribution limits. Warmer temperatures during the Holocene caused both species, C. atlantica and P. nigra, to migrate upwards [17]. In this situation, both species reduced their shared habitat and, even so, it is possible that their populations maintained contact. Thus, they could have been growing in the same stands as mixed coniferous forests; Cedrus could have formed part of pine forests as a companion species, or formed small separate copses in the same territory that P. nigra would have inhabited. Going up in altitude, the smaller distribution area of Cedrus in the Iberian Peninsula with respect to Pinus would make the species more sensitive to certain stochastic processes that would lead to habitat fragmentation and, in the extreme case, to its local disappearance, despite the existence of available habitat.
According to our results, the spatial segregation between species is much greater than that of its niches. Both species share their niches and both species have also shared habitat on the same continent since the Pleistocene. Nevertheless, their current distributions have drifted towards complete spatial segregation. It cannot be ruled out that this fact is due to the interaction between both species. In fact, the geographic segregation may have affected the results of the analysis by introducing a bias in the observations used to calibrate the distribution models. If there had been any observations of C. atlantica currently in Europe, the resulting model would likely have been different, with a much larger potential habitat distribution and greater overlap with the P. nigra niche.
In conclusion, C. atlantica and P. nigra show differentiated niches based on greater tolerance to extreme cold and continentality of Cedrus. Currently, the distribution of both species, segregated on different continents—P. nigra in Europe and C. atlantica in Africa—is virtually responding to dissociation from their optimal habitats. However, both species were able to share habitat in their distribution limits on the European continent until the end of the Pleistocene. The increase in temperatures of the Holocene Optimum Climate must have forced the rise of both species in altitude, causing the isolation of their populations in mountain refuges and the consequent loss of connectivity that has lasted until today. Although C. atlantica has continuously maintained its ideal habitat in the Baetic Mountains to date, it cannot be ruled out that the higher prevalence of P. nigra on mainland Europe has eventually displaced Cedrus.

Author Contributions

Conceptualization, A.G.-H., D.N.-L. and J.P.; methodology and formal analysis, A.G.-H. and D.N.-L.; writing—original draft preparation, A.G.-H.; writing—review and editing, D.N.-L., F.A.-S. and J.P.; project administration and funding acquisition, F.A.-S. All authors have read and agreed to the published version of the manuscript.

Funding

This research study was funded by (i) the Spanish government, State R&D Program Oriented to the Challenges of the Society, MED-REFUGIA Research Project (RTI2018-101714-B-I00); (ii) Andalusian Plan for Research, Development and Innovation, OROMEDREFUGIA Research Project (P18-RT-4963); (iii) ERDF Operational Programme in Andalusia (EU regional programme), RELIC-FLORA 2 Research Project (B-RNM-404-UGR18); and (iv) State Program for the Promotion of Scientific Research and Excellence Technique, PALEOPINSAPO Research Project (CSO2017-83576-P). The APC was funded by (i) the Spanish government, State R&D Program Oriented to the Challenges of the Society, MED-REFUGIA Research Project (RTI2018-101714-B-I00); and (ii) Andalusian Plan for Research, Development and Innovation, OROMEDREFUGIA Research Project (P18-RT-4963).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used in this study are available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A

Table A1. Correlation between environmental variables (Pearson). TCOL, min temperature of coldest month; TRAN, temperature annual range; TWET, temperature of wettest quarter; PDRM, precipitation of the driest month.
Table A1. Correlation between environmental variables (Pearson). TCOL, min temperature of coldest month; TRAN, temperature annual range; TWET, temperature of wettest quarter; PDRM, precipitation of the driest month.
TCOLTRANTWETPDRM
TCOL1
TRAN0.0551
TWET0.4640.4341
PDRM−0.687−0.504−0.4361
Table A2. Evaluation of the SDM of Cedrus atlantica projected into the past: observed and null model AUC.
Table A2. Evaluation of the SDM of Cedrus atlantica projected into the past: observed and null model AUC.
HCOLGM
CCSM4MIROC-ESMCCSM4MIROC-ESM
Obs. AUC0.73960.78300.50110.5425
Null AUC0.75840.73980.76160.7769
The observed AUC (area under the curve) values are only significantly higher (p < 0.05) than the null model AUC values in the models projected to the HCO with the MIROC-ESM general circulation model.
Figure A1. Observed and null model AUC in MaxEnt. The null model AUC is the 0.95 quantile in a random distribution. The observed AUC values (red diamond) were significantly (p < 0.05) higher than the AUC of the null models.
Figure A1. Observed and null model AUC in MaxEnt. The null model AUC is the 0.95 quantile in a random distribution. The observed AUC values (red diamond) were significantly (p < 0.05) higher than the AUC of the null models.
Diversity 14 00136 g0a1
Figure A2. Niche overlap test. Histograms show the observed D-value (red diamond) versus the null distribution in the equivalence and similarity tests.
Figure A2. Niche overlap test. Histograms show the observed D-value (red diamond) versus the null distribution in the equivalence and similarity tests.
Diversity 14 00136 g0a2
Figure A3. Density curves of the predictor variables for Pinus nigra, Cedrus atlantica and the background in the study area (see Figure 1).
Figure A3. Density curves of the predictor variables for Pinus nigra, Cedrus atlantica and the background in the study area (see Figure 1).
Diversity 14 00136 g0a3
Figure A4. Response curves of the four environmental variables in the models in MaxEnt. The differences in the pattern of the curves between the Pinus nigra and Cedrus atlantica models are consistent with the divergences observed in the distribution of their habitats. TCOL, min temperature of coldest month, TRAN, temperature annual range, TWET, temperature of wettest quarter, PDRM, precipitation of driest month.
Figure A4. Response curves of the four environmental variables in the models in MaxEnt. The differences in the pattern of the curves between the Pinus nigra and Cedrus atlantica models are consistent with the divergences observed in the distribution of their habitats. TCOL, min temperature of coldest month, TRAN, temperature annual range, TWET, temperature of wettest quarter, PDRM, precipitation of driest month.
Diversity 14 00136 g0a4

References

  1. Abel-Schaad, D.; Iriarte, E.; López-Sáez, J.A.; Pérez-Díaz, S.; Sabariego Ruiz, S.; Cheddadi, R.; Alba-Sánchez, F. Are Cedrus atlantica forests in the Rif Mountains of Morocco heading towards local extinction? Holocene 2018, 28, 1023–1037. [Google Scholar] [CrossRef] [Green Version]
  2. Magri, D.; Di Rita, F.; Aranbarri, J.; Fletcher, W.; González-Sampériz, P. Quaternary disappearance of tree taxa from Southern Europe: Timing and trends. Quat. Sci. Rev. 2017, 163, 23–55. [Google Scholar] [CrossRef] [Green Version]
  3. Postigo-Mijarra, J.M.; Morla, C.; Barrón, E.; Morales-Molino, C.; García, S. Patterns of extinction and persistence of Arctotertiary flora in Iberia during the Quaternary. Rev. Palaeobot. Palynol. 2010, 162, 416–426. [Google Scholar] [CrossRef]
  4. Postigo-Mijarra, J.M.; Barrón, E.; Gómez Manzaneque, F.; Morla, C. Floristic changes in the Iberian Peninsula and Balearic Islands (South-West Europe) during the Cenozoic. J. Biogeogr. 2009, 36, 2025–2043. [Google Scholar] [CrossRef] [Green Version]
  5. Magri, D. Quaternary History of Cedrus in Southern Europe. Ann. Bot. 2012, 2, 57–66. [Google Scholar] [CrossRef]
  6. Boyd, A. Relict Conifers from the Mid-Pleistocene of Rhodes, Greece. Hist. Biol. 2009, 21, 1–15. [Google Scholar] [CrossRef]
  7. Jiménez-Moreno, G.; Anderson, R.S.; Ramos-Román, M.J.; Camuera, J.; Mesa-Fernández, J.M.; García-Alix, A.; Jiménez-Espejo, F.J.; Carrión, J.S.; López-Avilés, A. The Holocene Cedrus pollen record from Sierra Nevada (S Spain), a proxy for climate change in N Africa. Quat. Sci. Rev. 2020, 242, 106468. [Google Scholar] [CrossRef]
  8. Magri, D.; Parra, I. Late Quaternary Western Mediterranean pollen records and African winds. Earth Planet. Sci. Lett. 2002, 200, 401–408. [Google Scholar] [CrossRef]
  9. Aranbarri, J.; González-Sampériz, P.; Valero-Garcés, B.; Moreno, A.; Gil-Romera, G.; Sevilla-Callejo, M.; García-Prieto, E.; Di Rita, F.; Mata, M.P.; Morellón, M.; et al. Rapid climatic changes and resilient vegetation during the Lateglacial and Holocene in a continental region of South-Western Europe. Glob. Planet. Chang. 2014, 114, 50–65. [Google Scholar] [CrossRef] [Green Version]
  10. González-Sampériz, P.; Leroy, S.A.G.; Carrión, J.S.; Fernández, S.; García-Antón, M.; Gil-García, M.J.; Uzquiano, P.; Valero-Garcés, B.; Figueiral, I. Steppes, savannahs, forests and phytodiversity reservoirs during the Pleistocene in the Iberian Peninsula. Rev. Palaeobot. Palynol. 2010, 162, 427–457. [Google Scholar] [CrossRef] [Green Version]
  11. González-Sampériz, P.; Gil-Romera, G.; García-Prieto, E.; Aranbarri, J.; Moreno, A.; Morellón, M.; Sevilla-Callejo, M.; Leunda, M.; Santos, L.; Franco-Múgica, F.; et al. Strong continentality and effective moisture drove unforeseen vegetation dynamics since the Last Interglacial at inland Mediterranean areas: The Villarquemado sequence in NE Iberia. Quat. Sci. Rev. 2020, 242, 106425. [Google Scholar] [CrossRef]
  12. Manzano, S.; Carrión, J.S.; López-Merino, L.; Ochando, J.; Munuera, M.; Fernández, S.; González-Sampériz, P. Early to Mid-Holocene spatiotemporal vegetation changes and tsunami impact in a paradigmatic coastal transitional system (Doñana National Park, Southwestern Europe). Glob. Planet. Chang. 2018, 161, 66–81. [Google Scholar] [CrossRef] [Green Version]
  13. Bell, B.A.; Fletcher, W.J. Modern Surface Pollen assemblages from the Middle and High Atlas, Morocco: Insights into pollen representation and transport. Grana 2016, 55, 286–301. [Google Scholar] [CrossRef] [Green Version]
  14. Cheddadi, R.; Henrot, A.-J.; François, L.; Boyer, F.; Bush, M.; Carré, M.; Coissac, E.; De Oliveira, P.E.; Ficetola, F.; Hambuckers, A.; et al. Microrefugia, climate change, and conservation of Cedrus atlantica in the Rif Mountains, Morocco. Front. Ecol. Evol. 2017, 5, 114. [Google Scholar] [CrossRef] [Green Version]
  15. Alba-Sánchez, F.; Abel-Schaad, D.; López-Sáez, J.A.; Sabariego Ruiz, S.; Pérez-Díaz, S.; González-Hernández, A. Paleobiogeografía de Abies spp. y Cedrus atlantica en el Mediterráneo Occidental (Península Ibérica y Marruecos). Rev. Ecosist. 2018, 27, 26–37. [Google Scholar] [CrossRef] [Green Version]
  16. Alba-Sánchez, F.; Abel-Schaad, D.; López-Sáez, J.A.; Sabariego-Ruiz, S.; Pérez-Díaz, S.; Luelmo-Lautenschlaeger, R.; Garrido-García, J.A. Early anthropogenic change in Western Mediterranean mountains (Sierra Nevada, SE Spain). Anthropocene 2021, 33, 100278. [Google Scholar] [CrossRef]
  17. Cheddadi, R.; Fady, B.; François, L.; Hajar, L.; Suc, J.-P.; Huang, K.; Demarteau, M.; Vendramin, G.G.; Ortu, E. Putative glacial refugia of Cedrus Atlantica deduced from Quaternary pollen records and modern genetic diversity. J. Biogeogr. 2009, 36, 1361–1371. [Google Scholar] [CrossRef]
  18. Linares, J.C.; Tíscar, P.A.; Camarero, J.J.; Taïqui, L.; Viñegla, B.; Seco, J.I.; Merino, J.; Carreira, J.A. Tree growth decline on relict Western-Mediterranean mountain forest: Causes and impacts. In Forest Decline: Causes and Impacts; Environmental Science, Engineering and Technology; Jenkins, J.A., Ed.; Nova Science Publishers: Hauppauge, NY, USA, 2011; pp. 91–110. ISBN 978-1-61470-002-9. [Google Scholar]
  19. Carrión, J.S.; Sánchez-Gómez, P.; Mota, J.F.; Yll, R.; Chaín, C. Holocene vegetation dynamics, fire and grazing in the Sierra de Gádor, Southern Spain. Holocene 2003, 13, 839–849. [Google Scholar] [CrossRef] [Green Version]
  20. Carrión, J.S.; Fuentes, N.; González-Sampériz, P.; Sánchez Quirante, L.; Finlayson, J.C.; Fernández, S.; Andrade, A. Holocene environmental change in a montane region of Southern Europe with a long history of human settlement. Quat. Sci. Rev. 2007, 26, 1455–1475. [Google Scholar] [CrossRef]
  21. Campbell, J.F.E.; Fletcher, W.J.; Joannin, S.; Hughes, P.D.; Rhanem, M.; Zielhofer, C. Environmental drivers of Holocene forest development in the Middle Atlas, Morocco. Front. Ecol. Evol. 2017, 5, 113. [Google Scholar] [CrossRef] [Green Version]
  22. Carrión, J.S. Paleoflora y Paleovegetación de la Península Ibérica e Islas Baleares; Universidad de Murcia, Secretariado PU: Murcia, Spain, 2012; ISBN 978-84-615-9026-1. [Google Scholar]
  23. Muller, S.D.; Rhazi, L.; Andrieux, B.; Bottollier-Curtet, M.; Fauquette, S.; Saber, E.-R.; Rifai, N.; Daoud-Bouattour, A. Vegetation History of the Western Rif Mountains (NW Morocco): Origin, Late-Holocene Dynamics and Human Impact. Veg. Hist. Archaeobotany 2015, 24, 487–501. [Google Scholar] [CrossRef]
  24. Reille, M. Analyse pollinique de sédiments postglaciaires dans le Moyen Atlas et le Haut Atlas Marocains: Premiers résultats. Ecol. Mediterr. 1976, 2, 153–170. [Google Scholar] [CrossRef]
  25. Reille, M. Contribution pollenanalytique à l’histoire holocène de la végétation des montagnes du Rif (Maroc septentrional). Suppl. Bull. AFEQ 1977, 50, 53–76. [Google Scholar]
  26. Stambouli-Essassi, S.; Roche, E.; Bouzid, S. Evolution de la végétation et du climat dans le Nord-ouest de la Tunisie au cours des 40 derniers millénaires. Geo-EcoTrop 2007, 31, 171–214. [Google Scholar]
  27. Blois, J.L.; Gotelli, N.J.; Behrensmeyer, A.K.; Faith, J.T.; Lyons, S.K.; Williams, J.W.; Amatangelo, K.L.; Bercovici, A.; Du, A.; Eronen, J.T.; et al. A Framework for evaluating the influence of climate, dispersal limitation, and biotic interactions using fossil pollen associations across the Late Quaternary. Ecography 2014, 37, 1095–1108. [Google Scholar] [CrossRef]
  28. D’Amen, M.; Mod, H.K.; Gotelli, N.J.; Guisan, A. Disentangling Biotic interactions, environmental filters, and dispersal limitation as drivers of species co-occurrence. Ecography 2018, 41, 1233–1244. [Google Scholar] [CrossRef] [Green Version]
  29. Broennimann, O.; Fitzpatrick, M.C.; Pearman, P.B.; Petitpierre, B.; Pellissier, L.; Yoccoz, N.G.; Thuiller, W.; Fortin, M.-J.; Randin, C.; Zimmermann, N.E.; et al. Measuring ecological niche overlap from occurrence and spatial environmental data. Glob. Ecol. Biogeogr. 2012, 21, 481–497. [Google Scholar] [CrossRef] [Green Version]
  30. Warren, D.L.; Glor, R.E.; Turelli, M. Environmental niche equivalency versus conservatism: Quantitative approaches to niche evolution. Evolution 2008, 62, 2868–2883. [Google Scholar] [CrossRef] [PubMed]
  31. Benito Garzón, M.; Sánchez de Dios, R.; Sáinz Ollero, H. Predictive modelling of tree species distributions on the Iberian Peninsula during the Last Glacial Maximum and Mid-Holocene. Ecography 2007, 30, 120–134. [Google Scholar] [CrossRef]
  32. López-Tirado, J.; Hidalgo, P.J. A high resolution predictive model for relict trees in the Mediterranean-mountain forests (Pinus sylvestris L., P. nigra Arnold and Abies pinsapo Boiss.) from the South of Spain: A Reliable management tool for reforestation. For. Ecol. Manag. 2014, 330, 105–114. [Google Scholar] [CrossRef]
  33. Navarro-Cerrillo, R.M.; Manzanedo, R.D.; Rodriguez-Vallejo, C.; Gazol, A.; Palacios-Rodríguez, G.; Camarero, J.J. Competition modulates the response of growth to climate in pure and mixed Abies pinsapo subsp. maroccana forests in Northern Morocco. For. Ecol. Manag. 2020, 459, 117847. [Google Scholar] [CrossRef]
  34. Valladares, A. Abetales de Abies pinsapo Boiss. In Bases Ecológicas Preliminares para la Conservación de los Tipos de Hábitat de Interés Comunitario en España; Ministerio Medio Ambiente y Medio Rural y Marino: Madrid, Spain, 2009; Volume 9520, ISBN 978-84-491-0911-9. [Google Scholar]
  35. Alba-Sánchez, F.; López-Sáez, J.A.; Benito de Pando, B.M.; Linares, J.C.; Nieto-Lugilde, D.; López-Merino, L. Past and present potential distribution of the Iberian Abies species: A phytogeographic approach using fossil pollen data and species distribution models. Divers. Distrib. 2010, 16, 214–228. [Google Scholar] [CrossRef] [Green Version]
  36. Blois, J.L. Recent advances in using species distributional models to understand past distributions. Front. Biogeogr. 2012, 3, 123–124. [Google Scholar] [CrossRef] [Green Version]
  37. Brewer, S.; Jackson, S.T.; Williams, J.W. Paleoecoinformatics: Applying geohistorical data to ecological questions. Trends Ecol. Evol. 2012, 27, 104–112. [Google Scholar] [CrossRef]
  38. Svenning, J.-C.; Fløjgaard, C.; Marske, K.A.; Nógues-Bravo, D.; Normand, S. Applications of species distribution modeling to paleobiology. Quat. Sci. Rev. 2011, 30, 2930–2947. [Google Scholar] [CrossRef]
  39. Varela, S.; Lobo, J.M.; Hortal, J. Using species distribution models in paleobiogeography: A matter of data, predictors and concepts. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2011, 310, 451–463. [Google Scholar] [CrossRef]
  40. Bouahmed, A.; Vessella, F.; Schirone, B.; Krouchi, F.; Derridj, A. Modeling Cedrus atlantica potential distribution in North Africa across time: New putative glacial refugia and future range shifts under climate change. Reg. Environ. Chang. 2019, 19, 1667–1682. [Google Scholar] [CrossRef]
  41. Demarteau, M.; Francois, L.; Cheddadi, R.; Roche, E. Réponses de Cedrus atlantica aux changements climatiques passés et futurs. Geo-Eco-Trop 2007, 31, 105–146. [Google Scholar]
  42. Benito Garzón, M.; Sánchez de Dios, R.; Sáinz Ollero, H. The evolution of the Pinus sylvestris L. area in the Iberian Peninsula from the Last Glacial Maximum to 2100 under climate change. Holocene 2008, 18, 705–714. [Google Scholar] [CrossRef]
  43. Di Pasquale, G.; Saracino, A.; Bosso, L.; Russo, D.; Moroni, A.; Bonanomi, G.; Allevato, E. Coastal pine-oak glacial refugia in the Mediterranean Basin: A biogeographic approach based on charcoal analysis and spatial modelling. Forests 2020, 11, 673. [Google Scholar] [CrossRef]
  44. Terrab, A.; Hampe, A.; Lepais, O.; Talavera, S.; Vela, E.; Stuessy, T.F. Phylogeography of North African Atlas cedar (Cedrus atlantica, Pinaceae): Combined molecular and fossil data reveal a complex Quaternary history. Am. J. Bot. 2008, 95, 1262–1269. [Google Scholar] [CrossRef] [PubMed]
  45. Thomas, P. Cedrus atlantica. The IUCN Red List of Threatened Species. e.T42303A2970716. 2013. Available online: https://www.iucnredlist.org/species/42303/2970716 (accessed on 27 November 2014).
  46. Farjon, A. Pinus nigra. The IUCN Red List of Threatened Species. e.T42386A2976817. 2013. Available online: https://www.iucnredlist.org/species/42386/2976817 (accessed on 13 January 2020).
  47. Phillips, S.J.; Dudık, M.; Shapire, R. Maxent Software for Species Habitat Modeling. 2013. Available online: http://biodiversityinformatics.amnh.org/open_source/maxent/ (accessed on 25 January 2013).
  48. Hijmans, R.J.; Phillips, S.J.; Leathwick, J.; Elith, J. Package ‘Dismo’. Species Distribution Modeling. R Package Version 0.7-17. 2014. Available online: https://rspatial.org/raster/sdm/ (accessed on 10 September 2014).
  49. Phillips, S.J. A Brief Tutorial on Maxent 2017. Available online: https://pdf4pro.com/view/a-brief-tutorial-on-maxent-american-museum-of-natural-5a321d.html (accessed on 2 April 2021).
  50. Hijmans, R.J.; Cameron, S.E.; Parra, J.L.; Jones, P.G.; Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 2005, 25, 1965–1978. [Google Scholar] [CrossRef]
  51. EPD European Pollen Database. Available online: http://www.europeanpollendatabase.net/data/ (accessed on 26 April 2012).
  52. Booth, T.H.; Nix, H.A.; Busby, J.R.; Hutchinson, M.F. Bioclim: The first species distribution modelling package, its early applications and relevance to most current MaxEnt studies. Divers. Distrib. 2014, 20, 1–9. [Google Scholar] [CrossRef]
  53. Elith, J.; Phillips, S.J.; Hastie, T.; Dudík, M.; Chee, Y.E.; Yates, C.J. A statistical explanation of MaxEnt for ecologists. Divers. Distrib. 2011, 17, 43–57. [Google Scholar] [CrossRef]
  54. Merow, C.; Smith, M.J.; Silander, J.A. A practical guide to MaxEnt for modeling species’ distributions: What it does, and why inputs and settings matter. Ecography 2013, 36, 1058–1069. [Google Scholar] [CrossRef]
  55. Rodríguez-Sánchez, F.; Hampe, A.; Jordano, P.; Arroyo, J. Past tree range dynamics in the Iberian Peninsula inferred through phylogeography and palaeodistribution modelling: A review. Rev. Palaeobot. Palynol. 2010, 162, 507–521. [Google Scholar] [CrossRef]
  56. Raes, N.; ter Steege, H. A null-model for significance testing of presence-only species distribution models. Ecography 2007, 30, 727–736. [Google Scholar] [CrossRef]
  57. Di Cola, V.; Broennimann, O.; Petitpierre, B.; Breiner, F.T.; D’Amen, M.; Randin, C.; Engler, R.; Pottier, J.; Pio, D.; Dubuis, A.; et al. Ecospat: An R package to support spatial analyses and modeling of species niches and distributions. Ecography 2017, 40, 774–787. [Google Scholar] [CrossRef]
  58. Camuera, J.; Jiménez-Moreno, G.; Ramos-Román, M.J.; García-Alix, A.; Toney, J.L.; Anderson, R.S.; Jiménez-Espejo, F.; Bright, J.; Webster, C.; Yanes, Y.; et al. Vegetation and climate changes during the last two glacial-interglacial cycles in the Western Mediterranean: A new long pollen record from Padul (Southern Iberian Peninsula). Quat. Sci. Rev. 2019, 205, 86–105. [Google Scholar] [CrossRef]
  59. Pons, A.; Reille, M. The Holocene- and Upper Pleistocene pollen record from Padul (Granada, Spain): A new study. Palaeogeogr. Palaeoclimatol. Palaeoecol. 1988, 66, 243–263. [Google Scholar] [CrossRef]
Figure 1. Study area (outlined in red) showing the distribution of the observations for Pinus nigra (blue, n = 257) and Cedrus atlantica (red, n = 211). The surrounding environment of each species, marked by dotted lines, were applied in the niche overlap analysis.
Figure 1. Study area (outlined in red) showing the distribution of the observations for Pinus nigra (blue, n = 257) and Cedrus atlantica (red, n = 211). The surrounding environment of each species, marked by dotted lines, were applied in the niche overlap analysis.
Diversity 14 00136 g001
Figure 2. Two-dimensional environmental space, determined by the two main axes of a PCA from four descriptor variables. TCOL, minimum temperature of the coldest month; TRAN, temperature annual range; TWET, mean temperature of the wettest quarter; PDRM, precipitation of the driest month.
Figure 2. Two-dimensional environmental space, determined by the two main axes of a PCA from four descriptor variables. TCOL, minimum temperature of the coldest month; TRAN, temperature annual range; TWET, mean temperature of the wettest quarter; PDRM, precipitation of the driest month.
Diversity 14 00136 g002
Figure 3. Potential habitat for Pinus nigra, Cedrus atlantica and their overlap in three time periods: (a) Last Glacial Maximum; (b) Holocene Climatic Optimum; (c) present. These were built from the species distribution model (MaxEnt) with environmental variables from the past according to the MIROC-ESM General Circulation Models.
Figure 3. Potential habitat for Pinus nigra, Cedrus atlantica and their overlap in three time periods: (a) Last Glacial Maximum; (b) Holocene Climatic Optimum; (c) present. These were built from the species distribution model (MaxEnt) with environmental variables from the past according to the MIROC-ESM General Circulation Models.
Diversity 14 00136 g003aDiversity 14 00136 g003bDiversity 14 00136 g003c
Figure 4. Niche overlap in environmental space. The solid line delimits the surrounding environment and the discontinuous line 50% of its area. The arrows show the relative displacement between the centroids of the niches in the overlap analysis. The overlap diagram shows the niches of Pinus nigra and Cedrus atlantica (polygons in light and dark green, respectively). The reddish polygon indicates the common niche (overlap).
Figure 4. Niche overlap in environmental space. The solid line delimits the surrounding environment and the discontinuous line 50% of its area. The arrows show the relative displacement between the centroids of the niches in the overlap analysis. The overlap diagram shows the niches of Pinus nigra and Cedrus atlantica (polygons in light and dark green, respectively). The reddish polygon indicates the common niche (overlap).
Diversity 14 00136 g004
Table 1. Evaluation of the distribution models: observed and null model AUC.
Table 1. Evaluation of the distribution models: observed and null model AUC.
SPECIESnAUCNull Model AUC
Pinus nigra2570.9650.591
Cedrus atlantica2110.9950.604
The observed AUC (area under the curve) values were significantly higher (p < 0.05) than the AUC values of the null model; n is the number of presences used in the calibration of the species distribution model and of each iteration in the null model.
Table 2. Niche overlap of Pinus nigra and Cedrus atlantica and equivalence and similarity tests. Similarity implies two asymmetric analyses (niche→surrounding environment; 1 and 2 represent P. nigra and C. atlantica, respectively).
Table 2. Niche overlap of Pinus nigra and Cedrus atlantica and equivalence and similarity tests. Similarity implies two asymmetric analyses (niche→surrounding environment; 1 and 2 represent P. nigra and C. atlantica, respectively).
SPECIESOverlap (D)EquivalenceSimilarity 1→2Similarity 2→1
P. nigra–C. atlantica0.365nsSimilar *Similar **
D, Schoener Overlap Index. * Significant at p < 0.05; ** Significant at p < 0.01; ns, the niches of the species are not significantly equivalent to each other.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

González-Hernández, A.; Nieto-Lugilde, D.; Alba-Sánchez, F.; Peñas, J. Biological Interaction as a Possible Ultimate Driver in the Local Extinction of Cedrus atlantica in the Iberian Peninsula. Diversity 2022, 14, 136. https://0-doi-org.brum.beds.ac.uk/10.3390/d14020136

AMA Style

González-Hernández A, Nieto-Lugilde D, Alba-Sánchez F, Peñas J. Biological Interaction as a Possible Ultimate Driver in the Local Extinction of Cedrus atlantica in the Iberian Peninsula. Diversity. 2022; 14(2):136. https://0-doi-org.brum.beds.ac.uk/10.3390/d14020136

Chicago/Turabian Style

González-Hernández, Antonio, Diego Nieto-Lugilde, Francisca Alba-Sánchez, and Julio Peñas. 2022. "Biological Interaction as a Possible Ultimate Driver in the Local Extinction of Cedrus atlantica in the Iberian Peninsula" Diversity 14, no. 2: 136. https://0-doi-org.brum.beds.ac.uk/10.3390/d14020136

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop