Lina Caballero-Villalobos
Francisco Fajardo-Gutiérrez
Mariasole Calbi
Gustavo A. Silva-Arias- 1Professorship for Population Genetics, Department of Life Science Systems, Technical University of Munich, Freising, Germany
- 2Jardín Botánico de Bogotá José Celestino Mutis, Bogotá, Colombia
- 3Botanischer Garten und Botanisches Museum Berlin, Institut für Biologie, Freie Universität Berlin, Berlin, Germany
It is predicted that climate change will strongly affect plant distributions in high elevation “sky islands” of tropical Andes. Polylepis forests are a dominant element of the treeline throughout the Andes Cordillera in South America. However, little is known about the climatic factors underlying the current distribution of Polylepis trees and the possible effect of global climate change. The species Polylepis quadrijuga is endemic to the Colombian Eastern Cordillera, where it plays a fundamental ecological role in high-altitude páramo-forest ecotones. We sought to evaluate the potential distribution of P. quadrijuga under future climate change scenarios using ensemble modeling approaches. We conducted a comprehensive assessment of future climatic projections deriving from 12 different general circulation models (GCMs), four Representative Concentration Pathways (R) emissions scenarios, and two different time frames (2041–2060 and 2061–2080). Additionally, based on the future projections, we evaluate the effectiveness of the National System of Protected Natural Areas of Colombia (SINAP) and Páramo Complexes of Colombia (PCC) in protecting P. quadrijuga woodlands. Here, we compiled a comprehensive set of observations of P. quadrijuga and study them in connection with climatic and topographic variables to identify environmental predictors of the species distribution, possible habitat differentiation throughout the geographic distribution of the species, and predict the effect of different climate change scenarios on the future distribution of P. quadrijuga. Our results predict a dramatic loss of suitable habitat due to climate change on this key tropical Andean treeline species. The ensemble Habitat Suitability Modeling (HSM) shows differences in suitable scores among north and south regions of the species distribution consistent with differences in topographic features throughout the available habitat of P. quadrijuga. Future projections of the HSM predicted the Páramo complex “Sumapaz-Cruz Verde” as a major area for the long-term conservation of P. quadrijuga because it provides a wide range of suitable habitats for the different evaluated climate change scenarios. We provide the first set of priority areas to perform both in situ and ex situ conservation efforts based on suitable habitat projections.
Introduction
The Andean region is a South American biodiversity hotspot, home to an impressively high diversity and endemism (Myers et al., 2000; Antonelli and Sanmartín, 2011). Regrettably, tropical Andean mountain species are classified as highly vulnerable to ongoing climate change, which may lead to strong shifts in distribution ranges, populations decline, and local extinction (Mittermeier et al., 2011; Mountain Research Initiative EDW Working Group, 2015). Although future distribution predictions for tropical Andean species are often challenging and carry a high uncertainty, due to the complex topography and climate of the region, they often point toward an extreme sensitivity to climate change (Lawler et al., 2009; Tovar et al., 2013; Tejedor Garavito et al., 2015; Peyre et al., 2020). Specifically, it is hypothesized that the main factors impacting future species distribution in the Andes will be the increase of temperature, changes in water fluxes, and the upward shift of the lower snowfall limit (Vuille et al., 2008, 2015; Herzog et al., 2011). Furthermore, several studies have already predicted an increase in threats and local extinction risks by the second half of the century in the Andes region (Ramirez-Villegas et al., 2014; Peyre et al., 2020; Valencia et al., 2020).
Polylepis forests are unique montane ecosystems, naturally occurring at high elevations in the Andes (Hoch and Körner, 2005; Navarro et al., 2005). These forests are generally present as sparse patches separated by vast areas of Páramo or puna grasslands (Fjeldså, 2002; Jameson and Ramsay, 2007) dominated mainly by Polylepis spp., that are endemic to South American mountain ranges, and tend to occur along the upper treeline (Simpson, 1979, 1986; Fjeldså et al., 1996). Polylepis treeline forests are the source of essential cultural, economic, and ecosystem services, providing pivotal ecological functions such as absorbing moisture from the clouds and releasing water in springs and rivers while being optimal habitats for other threatened and endemic species (Kessler, 2006; Donald et al., 2010). These forests require special conservation efforts (Purcell and Brelsford, 2004; Pinos, 2020) as they are threatened by deforestation, fires, farming, timber extraction, and land-use change (Fjeldså et al., 1996; Etter and Villa, 2000; Kessler, 2002; Hensen et al., 2011).
Several studies already support the importance of environmental variables such as elevation, precipitation, and mean annual temperature in shaping the current distribution patterns of different Polylepis species (Zutta et al., 2012; Boza Espinoza et al., 2019). Constraints related to water availability and soil water content during the dry season, hence precipitation and low-temperature extremes define the central eco-physiologic adaptations of Polylepis species (Velez et al., 1998; Morales et al., 2004; Ramos et al., 2013). Moreover, the decline of Polylepis diversity and abundance is one of the hypothesized effects of global or regional climate change patterns on Andean ecosystems (Gareca et al., 2010; Anderson et al., 2011). In particular, recent studies predicted severe reductions of suitable habitat for Polylepis tarapacana, P. incana, and P. reticulata before the end of the twenty-first century (Ramirez-Villegas et al., 2014; Cuyckens et al., 2016). It is known that closed forest stands can maintain their own microclimate and that their irregular cohort establishment depends on the occurrence of favorable (wet) years (Rada et al., 2011). The possible importance of moisture conditions for Polylepis spp. has already been reported for northern Chile (Rundel et al., 2003). Also, it has been postulated that any increase in temperature will further restrict Polylepis australis distribution in central Argentinian Andes to the uppermost elevations (Marcora et al., 2008), therefore underlining a high sensitivity to high temperatures for the genus. Therefore, understanding the consequences of climate change on Polylepis woodlands, as one of the most sensitive Andean biomes (Kessler et al., 2014), should be a priority for conservation biology in the Andean region, as necessary means to guide conservation policies and management strategies in the coming decades (Gosling et al., 2009; Sousa-Silva et al., 2014; Cuyckens et al., 2016). Also, it has to be taken into account that plausible niche specialization processes taking place in the different regions of the species distribution may lead to a reduced model ability to predict suitable areas under the assumption of niche equilibrium (Guisan and Thuiller, 2005; Parmesan, 2006; Hargreaves and Eckert, 2014).
In the specific case of the Colombian Andes, it is known that Polylepis quadrijuga Bitter was one of the first native species to be exploited (Rodríguez et al., 2018) and that continuous expansion of agricultural and livestock frontiers has led to increased fragmentation and isolation of P. quadrijuga woodlands. Thus, it is important to expand the understanding of the main ecological drivers of P. quadrijuga distribution, and to estimate the extent of its current distribution, while evaluating the effects of ongoing climate change on the species. Here, we selected this treeline species as a case study to address the influence of climate change on the distribution of key highland species in the northern Andean region. We apply Habitat Suitability Modeling (HSM) to infer the climate-change vulnerability of P. quadrijuga, one of the recorded species of Polylepis in Colombia (Boza Espinoza et al., 2019). Moreover, we specifically aim to (i) identify the key climate variables shaping the current geographical distribution of P. quadrijuga, (ii) estimate present and future suitable areas and, (iii) identify climatically stable areas as potential sites for the implementation of protection efforts. Additionally, we evaluate to what extent the current protected areas network and Páramo complexes of Colombia enclose the predicted future suitable areas under climate change for P. quadrijuga, and if these can be sufficient for the in situ conservation of this species. Finally, we discuss possible conservation strategies and effective management actions that can ensure the persistence of this Colombian iconic Andean treeline species, addressing the specific protected areas that could host suitable habitat for P. quadrijuga under current and future climate change conditions.
Materials and Methods
Study Species and Area
Polylepis quadrijuga (Rosaceae), locally known as “colorado” or “sietecueros” is an evergreen tree of up to 10 m height with compound leaves and a twisted trunk with continuously peeling reddish bark (Simpson, 1986). This species is easily recognized, its taxonomy resolved, and its distribution does not overlap with other Polylepis species (Kessler and Schmidt-Lebuhn, 2006; Segovia-Salcedo et al., 2018). P. quadrijuga is wind-pollinated and the seeds are dispersed by the wind during the dry season (Simpson, 1979; van Schaik et al., 1993; Schmidt-Lebuhn et al., 2007). Particularly, P. quadrijuga is endemic to the Colombian Eastern Cordillera slopes and rocky areas that are usually immersed within the Páramo grassland matrix (Simpson, 1986; Schmidt-Lebuhn et al., 2006; Rangel-Ch and Arellano-P, 2010; Pérez-Escobar et al., 2018; Peyre et al., 2018). Also, P. quadrijuga occurs often above the timberline and forms dense stands at elevations as high as 4,000 m, usually along streams (Simpson, 1979).
The Northern Andes present a high physiographic and topographic complexity, are considered one of the most diverse areas in the world (Myers et al., 2000; Jenkins et al., 2013; Jung et al., 2020), and are characterized by a very moist climate, low thermal seasonality, and marked diurnal temperature variations (Buytaert et al., 2006). Regionally, the climate of the Colombian Eastern Cordillera ranges from humid on the eastern flank to sub-humid or slightly dry on the western flank. The precipitation regime is bimodal, with pronounced dry and wet seasons, and annual rainfall ranging from 600 to 1,300 mm/y (Etter and Villa, 2000).
Polylepis species possess physiological adaptations to low temperatures and physiological drought that allow them to grow under colder conditions than the global mean for high elevation treeline species (Fjeldså et al., 1996; Kessler et al., 2014). Key physiological traits such as high foliar concentration of flavonoids and seasonal phenology could play an important role in the adaptation of this species to the severe environmental conditions of the Andean highlands (Velez et al., 1998). Despite these appealing features, Polylepis has been little studied. Moreover, as a result of the long history of deforestation in the high Andean forests of Colombia, the remaining patches of Polylepis forests are spread out over the open modified landscape (Etter and Villa, 2000; Rangel-Ch and Arellano-P, 2010), and several studies highlight that the relict patches of Polylepis forests in Colombia should be considered as a priority for conservation programs.
Species Occurrence Data
We compiled occurrence data for P. quadrijuga from three different sources: (i) the Global Biodiversity Information Facility (GBIF), initially 104 georeferenced records, where we also checked the nomenclature for synonyms (Borja-Acosta, 2017, 2019; Raz and Agudelo, 2019; García et al., 2020; Grant and Niezgoda, 2020; Marín and Moreno, 2020; Rodríguez Bolaños, 2020; Lozano Bernal, 2021; Magill et al., 2021); (ii) curated citizen science projects observations (iNaturalist, 49 records), which included geographical coordinates and a photographic record, allowing an accurate identification check (Ueda, 2021); and (iii) direct virtual exploration inside the main Páramo complexes of the Eastern Cordillera, using the Google Street View Imagery tool (Google Maps, 23 records). We surveyed all the rural roads with available street view imagery, crossing the Páramo districts “Boyacá,” “Cundinamarca,” and “Santanderes” (Sarmiento et al., 2013). Through direct inspection of hundreds of 3D georeferenced panoramas, we certified the occurrence of P. quadrijuga when we observed multi-trunk trees or shrubs showing the combination of vegetative characters as the distinctive reddish bark peeling in thin papery layers and compound leaves with grayish-green color1. No other group of plants could be mistakable with Polylepis in the Páramo landscape (Simpson, 1979), and P. quadrijuga is the only species of the genus recorded from the Eastern Cordillera of Colombia (Rangel-Ch and Arellano-P, 2010; Fajardo-Gutiérrez et al., 2018). Therefore we considered this methodological approach adapted from Hardion et al. (2016) to be adequate in this study.
Only records with global positioning system (GPS) coordinates and detailed localities were used, and then visually inspected for obvious geocoding errors. We included all the occurrences covering the entire known geographic distribution of the species (Figure 1 and Supplementary Table 1) except the ones belonging to the Páramo of Frontino, located in the Western Cordillera of Colombia in the Antioquia department, which were removed, as these are now believed to form a separate species, and are currently under taxonomic reevaluation (Boza and Kessler, pers. comm.). We removed records with uncertain locations and retained occurrences with a minimum distance of 1 km from each other to avoid spatial pseudo-replication, yielding a final set of 99 records: 37 from GBIF (41%), 40 from iNaturalist (40%) and 22 from Street View Imagery (22%) (Supplementary Table 1). We also evaluated potential geographical sampling bias in our species occurrence dataset due to differences in accessibility. For that, we assessed the biasing effect of roads on the sampling intensity using the Bayesian approach implemented in the R package sampbias (Zizka et al., 2020). We calculated the bias effect on the three different sources of observed records (GBIF, iNaturalist, and Google Street View Imagery) and ten simulated datasets of same size as our complete dataset under complete random sampling and other 10 highly biased simulated datasets limiting the random sampling to maximum 1 km from roadways within the species distribution range (both geographical and altitudinal). The biasing effects show that differences in the observations among the three data sources and the simulated full random sampling are much smaller compared to the random sampling around the roads (Supplementary Figure 1).
Figure 1. Habitat and distribution of Polylepis quadrijuga woodlands in the Colombian Eastern Cordillera. Left panel: habitat and architecture of the P. quadrijuga (top) and morphology of the plant in situ (bottom; Photos by F. Fajardo). Central panel: Projection of the P. quadrijuga ensemble Habitat Suitability Model to current climatic conditions (showing suitability values > 0.554). Continuous color scale represents the habitat suitability scores. The pink shading denotes the spatial extent of the Páramo complexes along the Colombian Eastern Cordillera sensu IAvH (2012). The protected areas of the Colombian National System of Protected Areas (SINAP) that overlap with Páramo regions are represented by the diagonal blue-lined areas. Right panel: broad geographical context of the main map in South America (top) and Colombia political boundaries (bottom) including the Andes Cordilleras (gray) and Eastern Cordillera’s Paramo complex areas (pink).
Environmental Data Sets and Abiotic Predictor Variables
We gathered a total of 23 biologically relevant climatic variables: 19 bioclimatic variables derived from monthly temperature and precipitation climatology for the years 1979–2013 available at 30 arc sec resolution (1 km at the equator) were retrieved from the CHELSA dataset (Karger et al., 2017)2; Slope and Aspect variables derived from a SRTM Digital Elevation Model using the raster package (Hijmans, 2015); and two topographic variables (Terrain roughness index and Topographic wetness index) at same resolution to the bioclimatic variables were obtained from the ENVIREM data set (Title and Bemmels, 2018). The bioclimatic variables have shown a high importance for defining the environmental niche of species of Polylepis (Zutta et al., 2012; Boza Espinoza et al., 2019) and are directly affected by climate change. The CHELSA climatic dataset has been proven to provide better results for topographically complex montane regions (Bobrowski and Schickhoff, 2017). Moreover, the topography-related variables that we incorporated, are secondary geographic predictors that are used to characterize spatial topographic diversity, terrain complexity and soil moisture patterns (Lin et al., 2006), and are very sensitive to local differences in elevation, being able to capture the presence of steep gradients and high spatial heterogeneity (Różycka et al., 2017).
We cropped and aligned the explanatory environmental layers to an area encompassing 0°–12.5° N and 79.5°–69.1° W, using the raster package (Hijmans, 2015) in R. This geographic extent includes the full known extant range of P. quadrijuga. In order to avoid multicollinearity-related noise, we performed variable selection using the function corselect of the R package fuzzySim (Barbosa, 2015). This function calculates Pairwise Pearson correlations among all explanatory variables and then, among each pair of highly correlated variables (correlation coefficient > 0.75), it retains the variable with the highest individual response in a preliminary binomial Generalized Linear model. The final explanatory data sets included 11 variables, all showing variation inflation factors below 10. Correlation values among the selected variables were calculated and plotted with the function ecospat.cor.plot of the package ecospat (Di Cola et al., 2017; Supplementary Figure 2).
Habitat Suitability Models
To perform the HSM we implemented an ensemble modeling framework (Araujo and New, 2007) using the BIOMOD package (Thuiller et al., 2009) in R version 4.0.3 (R Core Team, 2020). This approach combines models obtained by several algorithms to account for the uncertainty associated with each modeling technique (Thuiller et al., 2009, 2020). We selected seven modeling algorithms: (1) generalized linear models—GLM (McCullagh and Nelder, 1983); (2) generalized additive models—GAM (Hastie and Tibshirani, 1990); (3) generalized boosting models—GBM (Death, 2007); (4) classification tree analysis—CTA (Breiman et al., 1993); (5) multiple adaptive regression splines—MARS (Leathwick et al., 2005); (6) random forests—RF (Breiman, 2001); (7) artificial neural networks—ANN (Ripley, 1996); and (8) maximum entropy—MAXENT; (Phillips et al., 2006). We performed individual parameter tuning for each algorithm using the BIOMOD_tuning function.
Following Barve et al. (2011), we sampled the pseudo-absences within the accessible area for the species. We outlined the area including the Colombian Eastern Cordillera above 1,500 m of elevation and outside a 5 km buffer around the presence records (Supplementary Figure 3). Within this area, we randomly selected 10,000 pseudo-absences. Prevalence was maintained at 0.5, allowing the sum of presence-weights to be equal to the sum of pseudo-absence weights in the model calibration process (Barbet-Massin et al., 2012). We ran ten cross-validation replicates where presence records were split into sets of 75% for training and 25% for testing the models. For each cross-validation replicate we performed five runs with different pseudo-absence sets, therefore completing a total of 400 models. To assess the predictive performance of the models, we used both the threshold-independent area under the receiver operating characteristic curve statistic—AUC (Phillips et al., 2006) and the threshold-dependent true skill statistic—TSS (Allouche et al., 2006). Consensus models were obtained using a TSS-weighted average method to account for the predictive power of each model. Models with low predictive power (TSS < 0.7) were discarded.
As an additional characterization of the multi-dimensional occupied environmental space of P. quadrijuga, we extracted selected predictor values in correspondence of our species occurrences and performed a principal components analysis (PCA), with the function prcomp in R (R Core Team, 2020).
Evaluation of Future Range and Elevation Shifts and Representation in the Protected Areas System
We assessed the possible impacts of climate change on the distribution of P. quadrijuga through the projection of the HSMs to two time periods (2041–2060 and 2061–2080). For future projection we retrieved the projected bioclimatic variables from the CHELSA dataset, according to four greenhouse gas concentration scenarios, including low (RCP 2.6), moderate (RCP 4.5), high (RCP 6.0), and strong (RCP 8.5) increases in global radiative forcing (IPCC, 2013). These scenarios express the intensity of climatic changes by representative concentration pathways, using predicted values for radiative forcing for 2100, in comparison to the preindustrial era (Harris et al., 2014). The RCP2.6 scenario assumes CO2 concentrations of 450 ppm in 2100, that global mean temperatures will increase 0.2–1.8°C and has no strict reference in previous (Fourth Report) IPCC guidelines. RCP4.5 assumes 650 ppm CO2 and 1.0–2.6°C in 2100 and refers to the previous guideline scenario B1. RCP6.0 assumes 850 ppm CO2 and 1.3–3.2°C in 2100 and refers to the previous guideline scenario B2. RCP8.5 assumes 1350 ppm CO2 and 2.6–4.8°C in 2100 and refers to the A1F1 scenario of previous IPCC guidelines (Harris et al., 2014).
In addition, to consider the inherent uncertainty derived from each different general circulation model, we performed a comprehensive assessment including twelve different CMIP5 General Circulation Models (GCMs) [CESM1(CAM5), CSIRO-Mk3-6-0, FIO-ESM, GFDL-ESM2G, GFDL-ESM2M, GISS-E2-H, IPSL-CM5A-MR, MIROC-ESM-CHEM, MIROC5, MRI-CGCM3, NorESM1-M, and CCSM4]. For the resulting future climate projections, we generated consensus distribution maps for each time period and RCP scenarios, outlining the suitable areas predicted by at least three binary-transformed ensemble projections for the different GCMs. Future bioclimatic variables for all the GCMs, RCPs and time periods were retrieved from the CHELSA database. The topographic variables were kept invariant for all future projections.
We intersected the current and future HSMs projections to estimate the sensitivity to climate change at the species level, and assess the distribution range stability for all future RCPs scenarios and time periods. To this aim, we calculated the percentage of predicted area that remains stable, as well as range expansions or contractions in relation to the current projected area. Values of less than 100% indicate a decrease in environmental suitability over time (i.e., a high vulnerability of the species to climate change), whereas values greater than 100% indicate an increased extent of suitable habitat (i.e., a low vulnerability of the species to climate change).
In order to identify the specific Protected Areas (PAs) that could be relevant for the protection of P. quadrijuga in the future, we measured the overlapping area between the different projections of suitable habitat and the areas of 10 out of 36 Páramo complexes of Colombia (PCC), belonging to the Páramo districts “Boyacá,” “Cundinamarca,” and “Santanderes” (Sarmiento et al., 2013). Colombia has defined 17 Páramo districts and 36 Páramo complexes, by grouping natural areas that are divided by geographical features or human transformed landscapes that constitute special management and protection zones (Sarmiento et al., 2013). We also calculated the overlap between the suitable habitat projections and 36 protected areas of the National System of Protected Areas of Colombia (SINAP, 2020), including five National Natural Parks3. We used the PCC and the SINAP maps and calculated the portion of the suitable area for P. quadrijuga within these protected areas under the current conditions and future scenarios (RCP 2.6, RCP 4.5, RCP 6.0, and RCP 8.5), for the periods 2041–2060 and 2061–2080.
Results
Model Performance and Current Suitable Area Distribution
The projection of the ensemble-HSM to current climate conditions predicted a suitable area of 4444.9 km2, retrieving the highest and spatially continuous values of suitability throughout the southern range of the species distribution (SCV and CHG Table 1 and Figure 1), and less suitable and more fragmented areas in the northern range. The predicted potential current distribution is consistent with the known distribution of this species. Ten PCCs along the Colombian Eastern Cordillera contained suitable areas for P. quadrijuga (Table 1 and Figure 1). Less than a half of the predicted suitable area (46.1%, 1813.4 km2) was covered by the protected areas network of the SINAP, but most of the suitable areas (95.8%, 3765.3 km2) were located inside of the PCC boundaries, with only 4.2% was located outside of the PCC protection (163.4 km2) (Figure 1 and Table 1). The algorithms returning higher model performance metrics were GBM, GLM, Maxent, and RF and, the HSM returned overall high model evaluation metrics, with average AUC = 0.96 ± 0.04 and TSS = 0.9 ± 0.07 (Supplementary Figure 4), underscoring a good model performance in discriminating the species’ fundamental climatic niche.
Table 1. Number of occurrences of P. quadrijuga, and resulting suitable areas for each Páramo Complex of Colombia (PCC; Sarmiento et al., 2013).
Importance of Environmental Predictors
The PCA identified several important climatic variables that highlight differences in niche breadth through the species distribution. The analysis revealed differences in niche occupancy with a strong geographic signal represented by the different Páramos complexes (Figures 1, 2). The first axis of PCA (35% of the total variation) synthesized a distribution gradient between the eastern and western flanks of the Eastern Cordillera, with the variables Mean Diurnal Range and Precipitation of Warmest Quarter showing the highest loadings, as well as Temperature Seasonality and Topographic wetness index, but with an opposite direction (Figure 2). The second axis (ca. 21% of the variation) with Precipitation of Driest Quarter showing the highest loading, highlighted a north-south gradient, with the exception of the localities from the SCV Páramo that had a wide niche breadth along this axis (Figure 2).
Figure 2. Scatter plot of the principal components analysis (PCA) based on environmental space values occupied by Polylepis quadrijuga from Colombian Eastern Cordillera. Each different color dots represent the whole set of occurrences and the polygons highlight the occurrences superimposed to each Páramo complex. Niche area corresponds to the area of the hull relative to the largest polygon (P. quadrijuga = 1).
We found that Mean Temperature of Warmest Quarter (Bio 10) was the variable with consistently higher individual influence (mean = 77% ± 11.3) for all modeling techniques applied. Isothermality (Bio 3, mean = 51% ± 19.1) and Precipitation of Driest Quarter (Bio 17, mean = 45% ± 25.4) had also relatively high importance values (Supplementary Figure 5). Topographic variables appeared to be less informative in predicting the species occurrences but showed more relevance for differentiating the different geographical groups within the species (Figure 2 and Supplementary Figure 5).
Predicted Suitable Distribution of P. quadrijuga Under Climate Change
The potential geographic range projected for P. quadrijuga under different climate change scenarios presented strong reductions in suitable areas (Figure 3), and also in areas encompassed by PAs (Table 2). Our results predicted a substantial decline of the suitable habitat under all climate change scenarios by 2040s, the fragmentation of currently suitable areas, and the loss of most of the suitable habitat patches in the north of the Eastern Cordillera (Páramo districts “Boyacá” and “Santanderes”). In fact, the predicted future distribution of the suitable habitat for P. quadrijuga shows that only the population of SCV, would be able to persist under all the four different RCP scenarios (Figure 3), while for all the other areas, a strong reduction and fragmentation (or even disappearance) of suitable areas is predicted. Moreover, the HSMs predicted a potential loss of suitable area especially in the northern part of the Eastern Cordillera, even for the most optimistic scenario (RCP 2.6 for 2041–2060; Figure 3).
Figure 3. Projection of the Habitat Suitability Model (HMS) of P. quadrijuga to future climate conditions derived from twelve general circulation models (GCMs) for the periods (2041–2060 and 2061–2080) under four Representative Concentration Pathway (26, 45, 60, and 85). The discrete color scale represents the number of different projections derived from the 12GCMs that support the presence of the species for future climate scenarios.
Table 2. Change of suitable area for the projections of the ensemble Habitat Suitability Model of Polylepis quadrijuga to four climatic change scenarios (RCP) and two time periods.
For all RCPs and for both the evaluated time periods, a decrease of more than 80% of the suitable habitat for P. quadrijuga is evident, and underscores a general scenario of dramatic range reduction and loss. The effect of global climate change on the potential distribution of P. quadrijuga is particularly drastic under the most pessimistic scenario (RCP 8.5), leaving only small relictual patches of the current suitable area in the northern species range (SNC and PSB Páramo complexes; Figure 1), for the 2061–2080 period. Interestingly, for the 2041–2060 period, the HSM predicted a higher loss of suitable area under the RCP 4.5 scenario than under all the other RCPs, followed by the RCP 8.5 (Figure 3). The potential loss of climatically suitable areas ranges from 76.4% (RCP 2.6—2041–2060) to 98.5% (RCP 4.5 2041–2060), even for the more climatically stable regions (i.e., SCV Páramo complex; Figure 1).
Additionally, the ensemble projection of the HSMs to future climate conditions predicted a fast contraction of the altitudinal range around medium-high elevation values, together with more fragmentation of suitable areas and even local population extinction in several of the northernmost localities (Figures 3, 4).
Figure 4. Altitudinal range distribution of the projected suitable regions under current and future climate scenarios of P. quadrijuga for each of the three Páramo complex districts (Table 1). Santanderes (ALM, JSB), Boyacá (TBM, RRB, PSB, GLR, SNC, IGM), and Cundinamarca (SCV, SHG, GUE).
Distribution of P. quadrijuga in the Protected Areas Network
Within the predicted current potential distribution of P. quadrijuga (4444.9 km2), we identified about 2036.9 km2 (45.8%) overlaps SINAP network areas. Moreover, 4168.8 km2 (93.8% of the current suitable area) overlays Páramo complexes areas (Table 1). From the 99 presence records, 41 occur inside the protected areas of the SINAP. The Natural Parks Sumapaz and Pan de Azúcar el Consuelo are the protected areas with the highest number of records (16 and 11 occurrence points, respectively; Supplementary Table 2). Among the National Natural Parks (PNN) we found species records in the PNN Chingaza, PNN Sierra Nevada del Cocuy, and PNN Sumapaz. With regard to Regional Natural Parks (PNR), P. quadrijuga occurrences were recorded inside the PNR Pan de Azúcar el Consuelo, PNR Mutiscua Pamplona, and PNR Páramo de Santurbán.
The future projections of P. quadrijuga predicted a strong decline of protected areas covering current highly to moderately suitable habitats. The most dramatic shift of suitable areas was observed for the period 2041–2060, compared to the current distribution, affecting the northernmost portion of the species distribution. Although the species is of high conservation interest, the known occurrences located inside protected areas in Colombia are limited to only 41% of the reported occurrences (Supplementary Table 2).
Discussion
This study reinforces the growing evidence of significant range shifts and drastic reduction of high Andean plant species driven by global climate change. The predicted potential future distribution of the Colombian Eastern Cordillera highland endemic Polylepis quadrijuga, highlights a strong decline in its suitable habitat. Thus, the implementation of management and conservation strategies becomes necessary as a basis for ensuring P. quadrijuga persistence and the preservation of the entire ecosystems that support this species. To this aim, our study identifies the areas of highest vulnerability for P. quadrijuga under future climatic conditions and provides useful information about stable suitable areas, serving as a baseline for decision-makers to implement effective conservation and management strategies to preserve Polylepis forest in Colombia. Lastly, we underline the importance of Protected Natural Areas in conservation studies and planning in the higher Andes, and their relevance for mountain species whose distribution is often naturally fragmented and more susceptible to habitat loss.
Current Suitable Habitat and Climatic Predictors for the Distribution of P. quadrijuga
We present the most detailed model to date regarding the potential distribution of P. quadrijuga, which shows four main areas with high habitat suitability values (Figure 1). Compared to a previous assessment (Fajardo-Gutiérrez et al., 2018), we included an updated dataset of herbarium records and verified records from citizen science projects and street imagery sources, enhancing the model ability to characterize the species niche (Wisz et al., 2008). In addition, we applied a multiple algorithm ensemble modeling approach using a climatic dataset suited for montane regions, overall increasing the robustness in the identification of the climatic predictors of P. quadrijuga distribution and the spatial prediction of suitable areas.
The HSM predicts large suitable areas for P. quadrijuga in the western slopes of the Eastern Cordillera of Colombia, specifically in the SCV and CHG Páramo complexes, located in the southern range of the species, but lower and more fragmented suitable areas in the northern Páramos (JSB and ALM; Figure 1). In fact, our occurrences dataset shows a northward decrease in records, which may suggest a sampling bias problem (Table 1). However, we interpret this observed pattern of decreasing suitability toward the north as a natural response to the change of available niches driven by three interrelated factors: slope, topographic complexity, and soil water content. SCV and CHG Páramos host larger and topographically less complex valleys, allowing the establishment of bigger populations of P. quadrijuga. At higher latitudes, the suitable climatic belt for the species presents a more complex topography, with steep slopes and lower water availability. This restricts the species establishment to smaller forests on selected-facing slopes that provide a sheltering effect and enable the species to avoid long-term freezing, which limits tree growth, as well as seedling dehydration and photoinhibition caused by strong solar radiation (Toivonen et al., 2018).
Interestingly, the PCA grouped the populations occurring in the Páramos complexes of the larger Eastern Cordillera highland plateau (i.e., CHG and SCV) along similar values of the variables Topographic wetness index, Bio 4 and Bio 6, and separated them from the populations of the northeastern corner of the Eastern Cordillera (i.e., GLR, JSB, PSB, SNC), that were instead more related with variables such as Slope and Bio 3 (PC1, Figure 2), showing a temperature and topographic gradient. The PC2 highlighted a precipitation gradient spanning from the CHG Páramo (with high values of Bio 14, Bio 16, and Bio 19), to PSB, JSB, and SNC, that showed low values of those variables, and that in other words classify as atmospherically dry Páramos. Two slightly different trends were observed regarding the environmental niche of P. quadrijuga, one for the 35 records found in GLR, which forms part of the Boyacá Páramo district, and another for the 42 records of the SCV Páramo complex in the Cundinamarca district (Table 1 and Figure 2), where the variables Bio 2 and Bio 4 play an important explicative role in the PC1. The GLR Páramo complex had a higher temperature diurnal range (Bio 2) compared to SCV, and SCV had a higher temperature seasonality (Bio 4) compared to GLR (Figure 2).
An alternative explanation of the lower amount of occurrences of P. quadrijuga in the northernmost part of the species distribution (Boyacá and Santanderes Páramo districts) is that these areas, where environmental conditions are suitable to host high andean forest and Páramo communities, had undergone severe deforestation through the past decades, thus resulting in the impoverishment of the local flora. In fact, 22–47% of the land belonging to the ALM, GLR, GUE, IGM, JSB, and PSB Páramo complexes is highly transformed, and as of 2011, these hosted complessively 168 km2 of legal mining titles (Sarmiento et al., 2013). Oppositely, SCV and CHG Páramo complexes are also protected areas of the SINAP, hosting only 10 and 6% of transformed area, respectively, and 3556.4 km2 of natural ecosystems, mostly Páramos and high andean forests (Sarmiento et al., 2013). Also, the northern east corner of the Eastern Cordillera displays a remarkably high topographic complexity, due to the presence of numerous high peaks surrounding the deep Chicamocha inter-Andean canyon (Peyre et al., 2019; Valencia et al., 2020). As a consequence, we found that these Páramo complexes (SNC, GLR, ALM, and PSB) host smaller suitable areas for P. quadrijuga, as already pointed out for the Espeletia complex (Valencia et al., 2020).
Polylepis trees are exposed to the harsh climatic conditions that characterize highland areas. Constraints related to water availability, hence precipitation, soil water content during the dry season, and low-temperature extremes define the main eco-physiologic adaptations of Polylepis species (Velez et al., 1998; Morales et al., 2004; Ramos et al., 2013). Consistently, our results showed that temperature-related variables targeting extreme diurnal/annual ranges are among the better climatic predictors of P. quadrijuga occurrence (Supplementary Figure 5). These results are in agreement with previously published studies focusing on other Polylepis species showing that relatively warm diurnal and freezing night temperatures are expected to limit the distribution in Polylepis species (Purcell and Brelsford, 2004; Toivonen et al., 2014; Cuyckens et al., 2016; Zutta and Rundel, 2017). Field observations reports support that P. quadrijuga reproduces sexually inside the forest fragments, while in edge zones, clonal reproduction becomes the primary natural propagation strategy. This change relates to the scarce success of seedling establishment in ecotonal zones, where the open vegetation does not protect nor create favorable conditions for seedling survival.
Predicted Shifts and Putative Future Climatically Stable Areas for P. quadrijuga
According to our HSM future predictions, climate change can reduce the suitable area for P. quadrijuga by up to 80% for the most optimistic scenario (Figure 3 and Table 2). Notably, according to our HSMs, we expect P. quadrijuga to become severely threatened in the medium-term (2041–2060) and locally extinct in large part of its range under the pessimistic scenarios RCP 6.0 and 8.5 in the long-term (2061–2080). This issue is a concerning situation considering that high elevation tree species are often unable to colonize new habitats in higher (hence colder) areas (Feeley et al., 2011; Duque et al., 2015; Fadrique et al., 2018).
Several studies predict a higher sensitivity to climate change of highland species and mainly relate it to upward distribution shifts along the altitudinal ranges, thus leading to a contraction of the available suitable area (Thuiller et al., 2005; Morueta-Holme et al., 2015). In the case of P. quadrijuga, our results predict significant reductions in the projected area and an altitudinal range shrinkage to the uppermost edge as a possible consequence of climate change (Figures 3, 4). Our projections show that only 11–15% of the suitable habitat of P. quadrijuga will be able to persist under the different 2041–2060 RCP scenarios and that only 8–17% will persist under the 2061–2080 RCP scenarios. These results are compatible with the projection of a severe reduction of P. tarapacana in the South American Altiplano (Cuyckens et al., 2016) and with a climate vulnerability assessment of the Espeletia complex, which shows that several Páramos of the Eastern Cordillera are highly vulnerable to a changing climate (Valencia et al., 2020). Overall, our results agree with the prediction that ongoing climate change could compromise the long−term persistence of the genus Polylepis (Cuyckens et al., 2016; Pinos, 2020).
Warming temperatures and changes in precipitation patterns are causing many high Andean biomes to retreat upslope and suffer changes in their composition as higher-altitude environments undergo faster temperature changes than lower-altitude environments (Costion et al., 2015; Mountain Research Initiative EDW Working Group, 2015). In particular, for some species that will need to track the isotherms throughout the Tropical Andes, a general upslope migration of 600 m is predicted to happen by the year 2100 (Anderson et al., 2011). However, forests’ range shift dynamics in response to climate change are a complex process, and a possible Andean treeline upward displacement to areas occupied by Páramos can be limited by many factors, including high radiation levels, freezing events, and human disturbances (Rehm and Feeley, 2015a,b).
Even if P. quadrijuga forests’ extent will not decrease, their distribution is unlikely to increase at higher elevations. For example, some Polylepis woodlands from the Cordillera de Vilcanota (Perú) can be highly stable (Jameson and Ramsay, 2007; Rehm and Feeley, 2013, 2015b). However, since our result predicted an altitudinal range shrink under climate change conditions, the species reproduction needs to be considered because treeline stands may show slow and species-specific responses (Camarero and Gutiérrez, 2004; Dullinger et al., 2004). Interestingly, Toivonen et al. (2018) found that the topographic heterogeneity of the high-Andes could help P. subsericans forests to persist regionally, facilitating these treeline species by providing dispersal and establishment microsites, achievable by local relocation.
Importance of Protected Areas Network to the Conservation of P. quadrijuga
We observed that PAs covered almost 75% of all remaining P. quadrijuga forests in Colombia. However, in all the examined PAs, the total suitable area for this species tends to be reduced under future climate change scenarios, especially toward the northern Páramos complex (Figure 3). In this context, our results serve as a baseline to adequately protect Polylepis woodlands, while highlighting relevant areas to expand the current PAs network to include suitable areas for P. quadrijuga (e.g., SCV, GLR, SNC, and GUE; Figure 1). The future HSMs projections define some important climatically stable areas where P. quadrijuga could be able to persist. However, these stable areas are located principally inside the SCV and CHG Páramos, where the model predicts a lower decrease in suitable areas. Nonetheless, remaining Polylepis patches might also undergo further fragmentation and occupy sub-optimal conditions at the periphery of their future ranges. A severe concern underlines the areas located at the northeast edge of the Eastern Cordillera, where a dramatic decrease or even a complete loss of suitable habitat is predicted (Table 2).
Our study reiterated the value and importance of an adequate representation and protection of biological entities in the Páramo complexes network in Colombia. According to our predictions, the SCV area supports stable conditions within the intervened local landscape despite the regional biotic and abiotic change. Interestingly, Zutta and Rundel (2017) predicted refugial areas along the Andes for many Polylepis-associated species indicating that temporally stable areas will be crucial for conservation during anthropogenic-induced climate change, and that their protection and management may reflect the persistence of Polylepis woodlands in the future. Correspondingly, our findings suggest that conservation efforts targeting P. quadrijuga should be directed first toward the large National Natural Parks (Sumapaz, Chingaza, Sierra Nevada del Cocuy; Supplementary Table 2).
According to our future projections, the current PAs network will protect ∼36–48% of the extent of the predicted P. quadrijuga suitable area. Additionally, the few areas that could maintain part of its future suitable habitat (CE-CM-SCV) are currently strongly threatened by human activities. For example, only 142,112 of the 333,420 ha of the Sumapaz Páramo are currently protected as part of the Sumapaz Natural National Park (IAvH, 2012), and fracking projects currently threaten the area. Thus, for P. quadrijuga, habitat protection and connectivity restoration should become a central objective of conservation strategies. In this sense, management plans of already established protected areas such as SCV and CHG Páramos, which our study reports as climatically suitable, should receive special attention to preserve Polylepis woodlands. We recommend additional protection efforts in the north of the species distribution, where areas under protection are currently scarce.
Here we pinpoint that several areas within SCV, CHG, and SNC Páramos could be considered as stable suitable areas for P. quadrijuga, since our models predicted their persistence under some of the evaluated future time frames and for most of the evaluated RCPs. Interestingly, these Páramos coincide with identified high diversity and Pleistocene refuge areas (Cleef, 2008; Zutta and Rundel, 2017). These Páramos were already identified as future persistence areas for different highland plant species (Ramirez-Villegas et al., 2014; Valencia et al., 2020). Thus, we assume that our results and their conservation implications can be of high relevance for protecting the highland Andean biodiversity in the future. However, it is very worrying that Páramo ecosystems have long been at risk due to agriculture, mining and that now they are susceptible to climate change. For example, for the CHG complex, it was found that increasing temperatures and changing precipitation regimes will reduce the 39–52% of the currently suitable area for the Páramo ecosystem during the dry season, and 13–34% during the wet season (Cresso et al., 2020).
We suggest that five of the PAs that showed at least a minimal area of persistence in some of the future scenarios and time frames, should be initially targeted for management: “Parque Nacional Natural Chingaza,” “Parque Nacional Natural Sumapaz,” “Parque Nacional Natural Sierra Nevada del Cocuy,” “Parque Nacional Natural Pisba,” “Pan de Azúcar—El Consuelo,” and “Páramo de Guantivá y La Rusia” regional parks. Prioritizing and expanding small protected areas in the correspondent mountain ranges and PCCs will be essential for protecting this species in situ and ex situ, and pivotal to safeguard it from further habitat loss. Interestingly, our model projections to future scenarios highlighted several areas of climatic stability or even suitable habitat expansion in many protected areas that currently host no records of P. quadrijuga (e.g., Pisba and Páramo de las Oseras; Supplementary Table 2), or small private reserves (e.g., La Reserva, Villarica; Supplementary Table 2). These can be effective for ex situ conservation measures for the scattered mature P. quadrijuga individuals and seedlings, and can be particularly relevant to carry out seed germination or seedling transplant experiments in the field.
Implications and Recommendations for the Conservation of Polylepis Woodlands Under Climate Change
Currently, few sustainable management strategies have been proposed for Polylepis woodlands, and where implemented, they do not match the widespread management practices in Colombian rural areas (Gareca et al., 2010; Vásquez et al., 2014; Boza Espinoza et al., 2019). Since many Polylepis woodlands are located in densely populated areas and are under considerable pressure from human disturbance, it is urgent to identify hotspots of vulnerability and resilience and define management interventions. We also emphasize the need to maximize protection efforts targeted around treeline ecotones biodiversity in the Andes, as these are essential providers of ecosystem services and home to a diverse array of unique wildlife (Fjeldså, 1993; Arnal et al., 2014; Sevillano-Ríos et al., 2018).
We consider P. quadrijuga as an appropriate umbrella species to monitor the upward shift of treeline Andean ecosystems. However, it is essential to gather and consider more information about its interactions with other species for the selection of appropriate sites for protection. For example, high habitat specialization with strong confinement inside Polylepis forests was observed for some specialist bird species (Fjeldså, 1993; Cahill and Matthysen, 2007; Gareca et al., 2010; Meneses and Herrera, 2013). In agreement with Sevillano-Ríos and Rodewald (2017), we suggest protecting large remnants at lower elevations and maintaining all relictual patches of Polylepis forests irrespective of size become two key strategies constituting the cornerstones of future conservation efforts. Also, protected areas in the higher Andes must incorporate the study of potential species distributions and species responses to local climate changes.
Although it is known that the presence of a species inside PAs offers no guarantees for effective species long-term conservation (Hoffmann et al., 2019), protection efforts are the most successful tools to maximize landscape connectivity of Polylepis forest. Simultaneously, identifying areas of historical and potential future refugia or climatically stable areas is important to direct the focus of conservation efforts (Keppel et al., 2015). Moreover, studies predicting species’ future distribution considering the idiosyncratic species responses (i.e., phenotypic acclimation, rapid evolution, ecophysiological studies of selected species, and local adaptation patterns) are necessary to set the methodological baseline for the management of Polylepis woodlands inside PAs. Some examples include Ramos et al. (2013), showing that the physiological performance of P. quadrijuga was significantly affected by landscape fragmentation, and Velez et al. (1998) supporting that P. quadrijuga trees develop a high protection strategy against UV-B radiation, which tends to increase with altitude. However, more detailed studies are still needed to clarify relevant aspects of species biology and deliver comprehensive eco−evolutionary data in the context of global climate change.
Toivonen et al. (2018) suggested that a better understanding of the environmental niche preferences of Polylepis species is crucial to identifying potential climatic resilience areas for the Polylepis forests. Accordingly, we show that HSM could be crucial for possible conservation strategies to address changing environmental conditions. However, the incorporation of species−specific dispersal constraints and the influence of local topographic and edaphic factors are essential additions to projections of plant species distribution models in the future (Rull and Vegas-Vilarrúbia, 2017). Moreover, since treeline species responses to climate changes are species-specific, climate, land-use, fragmentation, and soil degradation are all critical factors determining whether an upwards migration is possible. For this reason, intensive fieldwork might be relevant to monitor plant altitudinal migration and could increase the knowledge of the potential responses of threatened treeline species in the higher Andes.
Under climate change conditions, plants face various barriers to their dispersal that may hamper gene-flow and recruitment, leading to extinction. Although Polylepis species are wind-pollinated, palynological and population genetics analyses show that pollen dispersal is limited to relatively short distances (Schmidt-Lebuhn et al., 2007), underlining that forest fragmentation can lead to an increase of genetic drift, inbreeding, and a consequent loss of adaptive capacity in local populations. On average, for the northern Andean region, is expected an increase in annual precipitation by values as high as 300 mm/year (Buytaert et al., 2011). This issue is crucial for biological processes such as Polylepis species dispersion. Also, the dispersion ability of tropical treeline species could be limited by the complex topography of the Andes and the harsh microclimatic conditions near the soil surface, which may be one of the reasons limiting seedlings recruitment (Wesche et al., 2008). We still have much to learn about the genetic mechanisms behind climate change adaptations of high-Andean treeline species, and the recognition of regional patterns of genetic diversity of Polylepis populations could offer the opportunity not only to preserve ecosystems with high biological diversity but also to understand and protect evolutionary processes from climate change (i.e., Schierenbeck, 2017).
In this work, we incorporated data deriving from citizen science projects to increase the number of occurrences of P. quadrijuga. We observed that integrating data from alternative databases could substantially increase sample size and could open up new opportunities to identify declining populations earlier to evaluate the effectiveness of policy interventions. We remark that citizen science has the potential to generate valuable biological observations suitable for species distribution modeling (Tiago et al., 2017; Roy-Dufresne et al., 2019; Petrovan et al., 2020) and ecological restoration projects (Edwards et al., 2018; Callaghan et al., 2019). Citizen science projects also provide the users with informal learning experiences that can contribute to the knowledge on ecological traits, seed collection, plant germination, and propagation of species of the higher Andes, which are often neglected in the scientific literature. However, our results should be interpreted with caution, given that the distribution of utilized species records is heavily biased toward the PAs and the central part of Colombian Eastern Cordillera. Probably, due also to the fact that patches of Polylepis are challenging to identify using current air photographs and satellite images (Aucca and Ramsay, 2005), most of our occurrences of P. quadrijuga are located in PAs situated and along roads, whereas large areas in the northeast and southeast of the Eastern Cordillera remain poorly sampled or not sampled at all.
Our model is subjected to limitations, since we did not account for physiological processes or any information about the species’ ecology and natural history (e.g., interspecific competition, migration rate), and recent human activity and disturbance (e.g., land use change, fires). Moreover, we could not incorporate data on microclimatic conditions into the modeling, due to the lack of available high resolution data for our study area. Nevertheless, an improved understanding of P. quadrijuga distribution and the possible future response to climate change represents a precondition for more precise modeling efforts of potential range shifts of treeline species in the Andes. Our results provide a better comprehension of the role of environmental constraints in shaping treeline species potential ranges. Specifically, we highlight a “Change in the climatically suitable area” and a “Change in average climatic suitability in already occupied cells” for P. quadrijuga, which have been defined as two quantifying components of risk for woody species under climate change (Ohlemüller et al., 2006). For this reason, based on our results, we suggest implementing an evaluation criterion to categorize the conservation status of P. quadrijuga, targeting a key treeline forming species on the northern Andes, a region that is still being heavily transformed and fragmented by agricultural and mining activities. Immediate implementation of in situ conservation programs could be crucial to ensure the persistence of the species. Monitoring the upward distributional shift would help us to detect the predicted decline in Polylepis populations over time. Lastly, the creation of a germplasm bank from specimens from all the different localities would be fundamental.
Conclusion
According to our future projections, the threat status of P. quadrijuga is alarming and needs attention due to the predicted considerable reduction in climatically suitable habitat for the species. The detailed information provided here could represent a valuable tool to guide the future establishment of new and efficient conservation efforts. Moreover, our findings can help restore Polylepis forests and preserve them through in situ and ex situ strategies, planned according to the predicted future stable climatic areas. Medium- and long-term conservation strategies for P. quadrijuga woodlands should be focused on actions to prevent deforestation of the extant patches and to implement ecological restoration measures in the Cruz Verde—Sumapaz and Páramo complex, which is the only predicted future suitable area for six out of eight of the future scenarios employed here. Our results can be used to define high-priority areas in the Andean region for conservation and management plans for Polylepis forest and can offer a baseline for similar analyses of other endangered and threatened treeline species in the region. Finally, we reinforce the need to strengthen management strategies inside those protected areas that comprise the remaining patches of P. quadrijuga forest in Colombia.
Data Availability Statement
The original occurrence data used in the study is included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
GS-A and LC-V designed the study and conducted the formal analysis. FF-G and MC performed data curation. GS-A, FF-G, and MC prepared the figures, tables, and Supplementary Data. LC-V and MC wrote the first draft of the manuscript. All authors interpreted the results, contributed substantially to revisions and worked together to define the structure and content of the manuscript, and gave final approval for publication.
Funding
The authors acknowledge the Technical University of Munich Publishing Fund for covering the costs for the publication of this article.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.2021.661550/full#supplementary-material
Footnotes
- ^ https://goo.gl/maps/AFXnvdGFjrtxPgkJ6
- ^ http://chelsa-climate.org/bioclim/
- ^ https://runap.parquesnacionales.gov.co
References
Allouche, O., Tsoar, A., and Kadmon, R. (2006). Assessing the accuracy of species distribution models: prevalence, kappa and the true skill statistic (TSS). J. Appl. Ecol. 43, 1223–1232. doi: 10.1111/j.1365-2664.2006.01214.x
Anderson, E. P., Marengo, J., Villalba, R., Halloy, S., Young, B., Cordero, D., et al. (2011). “Consequences of climate change for ecosystems and ecosystem services in the tropical Andes,” in Climate Change and Biodiversity in the Tropical Andes, eds S. K. Herzog, R. Martínez, P. M. Jørgensen, and H. Tiessen (São José dos Campos, Brazil: Inter-American Institute for Global Change Research (IAI) and Scientific Committee on Problems of the Environment).
Antonelli, A., and Sanmartín, I. (2011). Why are there so many plant species in the Neotropics? Taxon 60, 403–414. doi: 10.1002/tax.602010
Araujo, M., and New, M. (2007). Ensemble forecasting of species distributions. Trends Ecol. Evol. 22, 42–47. doi: 10.1016/j.tree.2006.09.010
Arnal, H., Sampson, A., Navarro, G., Palomino, W., Ferreira, W., Romoleroux, K., et al. (2014). Mapa Pan Andino de Bosques de Polylepis prioritarios para Conservación. American Bird Conservancy. USA: The Plains.
Aucca, C., and Ramsay, P. M. (2005). Management of biodiversity and land use in Southern Peru: ECOAN’s activities to help conserve Polylepis woodlands. Mountain Res. Dev. 25, 287–289. doi: 10.1659/0276-47412005025[0287:MOBALU]2.0.CO;2
Barbet-Massin, M., Jiguet, F., Albert, C. H., and Thuiller, W. (2012). Selecting pseudo-absences for species distribution models: how, where and how many?: how to use pseudo-absences in niche modelling? Methods Ecol. Evol. 3, 327–338. doi: 10.1111/j.2041-210X.2011.00172.x
Barbosa, A. M. (2015). fuzzySim: Applying fuzzy logic to binary similarity indices in ecology. Methods Ecol. Evol. 6, 853–858. doi: 10.1111/2041-210X.12372
Barve, N., Barve, V., Jiménez-Valverde, A., Lira-Noriega, A., Maher, S. P., Peterson, A. T., et al. (2011). The crucial role of the accessible area in ecological niche modeling and species distribution modeling. Ecol. Model. 222, 1810–1819. doi: 10.1016/j.ecolmodel.2011.02.011
Bobrowski, M., and Schickhoff, U. (2017). Why input matters: Selection of climate data sets for modelling the potential distribution of a treeline species in the Himalayan region. Ecol. Model. 359, 92–102. doi: 10.1016/j.ecolmodel.2017.05.021
Borja-Acosta, K. G. (2019). Plantas de la expedición en el páramo de Pisba, Boyacá - Proyecto Boyacá BIO. Instituto de Investigación de Recursos Biológicos Alexander von Humboldt. Occurrence dataset. Available online at: https://doi.org/10.15472/cnykdj accessed via GBIF.org on 2021-03-14 (accessed March 14, 2021).
Borja-Acosta, K. (2017). Colección Herbario Federico Medem Bogotá - FMB. Version 32.0. Instituto de Investigación de Recursos Biológicos Alexander von Humboldt. Occurrence dataset. Available online at : https://doi.org/10.15472/ighftu accessed via GBIF.org on 2021-03-14 (accessed March 14, 2021).
Boza Espinoza, T. E., Quispe-Melgar, H. R., and Kessler, M. (2019). Taxonomic reevaluation of the Polylepis sericea complex (Rosaceae), with the description of a new species. Syst. Bot. 44, 324–334. doi: 10.1600/036364419X15562052252225
Breiman, L., Friedman, J., Stone, C. J., and Olshen, R. A. (1993). Classification and regression trees. The Wadsworth statistics/probability series. New York, NY: Chapman & Hall.
Buytaert, W., Cuesta-Camacho, F., and Tobón, C. (2011). Potential impacts of climate change on the environmental services of humid tropical alpine regions. Glob. Ecol. Biogeogr. 20, 19–33. doi: 10.1111/j.1466-8238.2010.00585.x
Buytaert, W., Celleri, R., De Bievre, B., Wyseure, G., Deckers, J., and Hofstede, R. (2006). Human impact on the hydrology of the Andean Páramos. Earth-Sci. Rev. 79, 53–72. doi: 10.1016/j.earscirev.2006.06.002
Cahill, J. R. A., and Matthysen, E. (2007). Habitat use by two specialist birds in high-Andean Polylepis forests. Biol. Conserv. 140, 62–69. doi: 10.1016/j.biocon.2007.07.022
Camarero, J. J., and Gutiérrez, E. (2004). Pace and pattern of recent treeline dynamics: response of ecotones to climatic variability in the Spanish Pyrenees. Clim. Change. 63, 181–200.
Cleef, A. M. (2008). “Humid cloud super Páramo probably acts as a plant diversity centre and as a cool refuge: the case of Nevado de Sumapaz, Colombia,” in La Cordillera Oriental Colombiana, Transecto Sumapaz, eds T. van der Hammen, J. O. Rangel, and A. M. Cleef (Berlin-Stuttgart: J. Cramer).
Callaghan, C. T., Major, R. E., Lyons, M. B., Martin, J. M., Wilshire, J. H., Kingsford, R. T., et al. (2019). Using citizen science data to define and track restoration targets in urban areas. J. Appl. Ecol. 2019, 1998–2006. doi: 10.1111/1365-2664.13421
Costion, C. M., Simpson, L., Pert, P. L., Carlsen, M. M., Kress, W. J., and Crayn, D. (2015). Will tropical mountaintop plant species survive climate change? Identifying key knowledge gaps using species distribution modeling in Australia. Biol. Conserv. 191, 322–330. doi: 10.1016/j.biocon.2015.07.022
Cresso, M., Clerici, N., Sanchez, A., and Jaramillo, F. (2020). Future climate change renders unsuitable conditions for Páramo ecosystems in Colombia. Sustainability 12:8373. doi: 10.3390/su12208373
Cuyckens, G. A. E., Christie, D. A., Domic, A. I., Malizia, L. R., and Renison, D. (2016). Climate change and the distribution and conservation of the world’s highest elevation woodlands in the South American Altiplano. Glob. Environ. Change. 137, 79–87. doi: 10.1016/j.gloplacha.2015.12.010
Death, G. (2007). Boosted trees for ecological modeling and prediction. Ecology 88, 243–251. doi: 10.1890/0012-9658200788[243:BTFEMA]2.0.CO;2
Di Cola, V., Broennimann, O., Petitpierre, B., Breiner, F. T., D’Amen, M., Randin, C., et al. (2017). ecospat: an R package to support spatial analyses and modeling of species niches and distributions. Ecography 40, 774–787. doi: 10.1111/ecog.02671
Donald, P., Collar, N., Marsden, S., and Pain, D. (2010). Facing extinction: the world’s rarest birds and the race to save them. London: Y & AD Poyser.
Dullinger, S., Dirnböck, T., and Grabherr, G. (2004). Modelling climate change-driven treeline shifts: relative effects of temperature increase, dispersal and invasibility. J. Ecol. 92, 241–252. doi: 10.1111/j.0022-0477.2004.00872.x
Duque, A., Stevenson, P. R., and Feeley, K. J. (2015). Thermophilization of adult and juvenile tree communities in the northern tropical Andes. Proc. Natl. Acad. Sci. U.S.A. 112, 10744–10749. doi: 10.1073/pnas.1506570112
Edwards, P. M., Shaloum, G., and Bedell, D. (2018). A unique role for citizen science in ecological restoration: a case study in streams: Citizen science and ecological restoration. Restor. Ecol. 26, 29–35. doi: 10.1111/rec.12622
Etter, A., and Villa, L. A. (2000). Andean forests and farming systems in part of the Eastern Cordillera (Colombia). Mt. Res. Dev. 20, 236–245. doi: 10.1659/0276-47412000020[0236:AFAFSI]2.0.CO;2
Fadrique, B., Báez, S., Duque, Á, Malizia, A., Blundo, C., Carilla, J., et al. (2018). Widespread but heterogeneous responses of Andean forests to climate change. Nature 564, 207–212. doi: 10.1038/s41586-018-0715-9
Fajardo-Gutiérrez, F., Infante-Betancour, J., and Cabrera-Amaya, D. M. (2018). Modelización de la distribución potencial del género Polylepis en Colombia y consideraciones para su conservación. Ecol. Austral. 28, 202–215. doi: 10.25260/EA.18.28.1.1.585
Feeley, K. J., Silman, M. R., Bush, M. B., Farfan, W., Cabrera, K. G., Malhi, Y., et al. (2011). Upslope migration of Andean trees: Andean trees migrate upslope. J. Biogeogr. 38, 783–791. doi: 10.1111/j.1365-2699.2010.02444.x
Fjeldså, J. (2002). Polylepis forests–vestiges of a vanishing ecosystem in the Andes. Ecotropica 8, 111–123.
Fjeldså, J. (1993). The avifauna of the Polylepis woodlands of the Andean highlands: the efficiency of basing conservation priorities on patterns of endemism. Bird Conserv. Int. 3, 37–55. doi: 10.1017/S0959270900000770
Fjeldså, J., Kessler, M., Engblom, G., and Driesch, P. (1996). Conserving the biological diversity of Polylepis woodlands of the highland of Peru and Bolivia: A contribution to sustainable natural resource management in the Andes. Copenhagen: Nordeco.
Gareca, E. E., Hermy, M., Fjeldså, J., and Honnay, O. (2010). Polylepis woodland remnants as biodiversity islands in the Bolivian high Andes. Biodiv. Conserv. 19, 3327–3346. doi: 10.1007/s10531-010-9895-9
García, N., León-Gómez, M., Jácome-Reyes, J. H., Herrera-Collazos, E. E., and Rodríguez-Morales, M. A. (2020). Herbario Pontificia Universidad Javeriana. Version 4.4. Pontificia Universidad Javeriana. Occurrence. Available online at: https://doi.org/10.15472/oyj27o (accessed March 14, 2021).
Gosling, W. D., Hanselman, J. A., Knox, C., Valencia, B. G., and Bush, M. B. (2009). Long-term drivers of change in Polylepis woodland distribution in the central Andes. J. Veg. Sci. 20, 1041–1052. doi: 10.1111/j.1654-1103.2009.01102.x
Grant, S., and Niezgoda, C. (2020). Field Museum of Natural History (Botany) Seed Plant Collection. Version 11.12. Field Museum. Available online at: https://doi.org/10.15468/nxnqzf (accessed March 14, 2021).
Guisan, A., and Thuiller, W. (2005). Predicting species distribution: offering more than simple habitat models. Ecol. Lett. 8, 993–1009. doi: 10.1111/j.1461-0248.2005.00792.x
Hardion, L., Leriche, A., Schwoertzig, E., and Millon, A. (2016). Species distribution 2.0: An accurate time- and cost-effective method of prospection using street view imagery. PLoS One 11:e0146899. doi: 10.1371/journal.pone.0146899
Hargreaves, A. L., and Eckert, C. G. (2014). Evolution of dispersal and mating systems along geographic gradients: implications for shifting ranges. Funct. Ecol. 28, 5–21. doi: 10.1111/1365-2435.12170
Harris, R. M. B., Grose, M. R., Lee, G., Bindoff, N. L., Porfirio, L. L., and Fox-Hughes, P. (2014). Climate projections for ecologists. WIREs Clim. Chang. 5, 621–637. doi: 10.1002/wcc.291
Hastie, T. J., and Tibshirani, R. J. (1990). Generalized additive models, monographs on statistics & Applied Probability. London: Chapman & Hall/CRC.
Hensen, I., Cierjacks, A., Hirsch, H., Kessler, M., Romoleroux, K., Renison, D., et al. (2011). Historic and recent fragmentation coupled with altitude affect the genetic population structure of one of the world’s highest tropical tree line species. Glob. Ecol. Biogeogr. 21, 455–464. doi: 10.1111/j.1466-8238.2011.00691.x
Herzog, S. K., Martínez, R., Jørgensen, P. M., and Tiessen, H. (eds). (2011). Climate Change and Biodiversity in the Tropical Andes, Inter-American Institute for Global Change Research (IAI) and Scientific Committee on Problems of the Environment (SCOPE), 348. doi: 10.13140/2.1.3718.4969
Hoffmann, S., Irl, S. D. H., and Beierkuhnlein, C. (2019). Predicted climate shifts within terrestrial protected areas worldwide. Nat Commun 10, 4787. doi: 10.1038/s41467-019-12603-w
Hoch, G., and Körner, C. (2005). Growth, demography and carbon relations of Polylepis trees at the world’s highest treeline. Funct Ecol. 2005, 941–951. doi: 10.1111/j.1365-2435.2005.01040.x
IAvH (2012). Cartografía de Páramos de Colombia Esc. 1:100.000. Proyecto: Actualización del Atlas de Páramos de Colombia. D.C. Colombia: Convenio Interadministrativo de Asociación 11-103, Instituto de Investigación de Recursos Biológicos Alexander von Humboldt y Ministerio de Ambiente y Desarrollo Sostenible. Bogotá.
IPCC (2013). Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, MA: Cambridge University Press.
Jameson, J. S., and Ramsay, P. M. (2007). Changes in high-altitude Polylepis forest cover and quality in the Cordillera de Vilcanota, Perú, 1956–2005. Biol. Conserv. 138, 38–46. doi: 10.1016/j.biocon.2007.04.008
Jenkins, C. N., Pimm, S. L., and Joppa, L. N. (2013). Global patterns of terrestrial vertebrate diversity and conservation. Proc. Natl. Acad. Sci. U.S.A. 110, 2602–2610. doi: 10.1073/pnas.1302251110
Jung, M., Arnell, A., de Lamo, X., García-Rangel, S., Lewis, M., Mark, J., et al. (2020). Areas of global importance for terrestrial biodiversity, carbon, and water (preprint). bioRxiv 2020:021444. doi: 10.1101/2020.04.16.021444
Karger, D. N., Conrad, O., Böhner, J., Kawohl, T., Kreft, H., Soria-Auza, R. W., et al. (2017). Climatologies at high resolution for the earth’s land surface areas. Sci. Data 4:170122. doi: 10.1038/sdata.2017.122
Keppel, G., Mokany, K., Wardell-Johnson, G. W., Phillips, B. L., Welbergen, J. A., and Reside, A. E. (2015). The capacity of refugia for conservation planning under climate change. Front. Ecol. Environ. 13:106–112. doi: 10.1890/140055
Kessler, M. (2006). “Bosques de Polylepis,” in Botánica Económica de los Andes Centrales, eds R. Mónica Moraes, B. Øllgaard, L. P. Kvist, F. Borchsenius, and H. Balslev (La Paz: Universidad Mayor de San Andrés).
Kessler, M., and Schmidt-Lebuhn, A. N. (2006). Taxonomical and distributional notes on Polylepis (Rosaceae). Organ. Diver. Evol. 6, 67–69. doi: 10.1016/j.ode.2005.04.001
Kessler, M., Toivonen, J. M., Sylvester, S. P., Kluge, J., and Hertel, D. (2014). Elevational patterns of Polylepis tree height (Rosaceae) in the high Andes of Peru: Role of human impact and climatic conditions. Front. Plant Sci. 5:194. doi: 10.3389/fpls.2014.00194
Lawler, J. J., Shafer, S. L., White, D., Kareiva, P., Maurer, E. P., Blaustein, A. R., et al. (2009). Projected climate-induced faunal change in the Western Hemisphere. Ecology 90, 588–597. doi: 10.1890/08-0823.1
Leathwick, J. R., Rowe, D., Richardson, J., Elith, J., and Hastie, T. (2005). Using multivariate adaptive regression splines to predict the distributions of New Zealand’s freshwater diadromous fish. Freshw. Biol. 50, 2034–2052. doi: 10.1111/j.1365-2427.2005.01448.x
Lin, H. S., Kogelmann, W., Walker, C., and Bruns, M. A. (2006). Soil moisture patterns in a forested catchment: A hydropedological perspective. Geoderma 131, 345–368. doi: 10.1016/j.geoderma.2005.03.013
Lozano Bernal, M. F. (2021). Herbario Museo de La Salle Bogotá (BOG). Version 3.5. Universidad de la Salle. Available online at: https://doi.org/10.15472/ppzpea (accessed March 14, 2021).
Magill, B., Solomon, J., and Stimmel, H. (2021). Tropicos Specimen Data. Missouri Botanical Garden. Available online at: https://doi.org/10.15468/hja69f (accessed March 14, 2021).
Marcora, P., Hensen, I., Renison, D., Seltmann, P., and Wesche, K. (2008). The performance of Polylepis australis trees along their entire altitudinal range: implications of climate change for their conservation. Divers. Distrib. 14, 630–636. doi: 10.1111/j.1472-4642.2007.00455.x
Marín, C., and Moreno, D. (2020). Colección general de plantas vasculares del Herbario JBB. Version 2.2. Jardín Botánico de Bogotá “José Celestino Mutis”. Available online at: https://doi.org/10.15472/mpvbpd (accessed March 14, 2021).
McCullagh, P., and Nelder, J. A. (1983). Generalized linear models, monographs on statistics and applied probability. London: Chapman and Hall/CRC.
Meneses, L., and Herrera, M. (2013). Estudio preliminar de la avifauna asociada a parches de Polylepis quadrijuga (Rosaceae) del Páramo de la Rusia, Duitama (Boyacá-Colombia). Luna Azul 36, 40–54.
Mittermeier, R. A., Turner, W. R., Larsen, F. W., Brooks, T. M., and Gascon, C. (2011). “Global biodiversity conservation: The critical role of hotspots,” in Biodiversity Hotspots, eds F. Zachos and J. Habel (Berlin: Springer), doi: 10.1007/978-3-642-20992-5_1
Morales, M. S., Villalba, R., Grau, H. R., and Paolini, L. (2004). Rainfall-controlled tree growth in high-elevation subtropical treelines. Ecology 85, 3080–3089. doi: 10.1890/04-0139
Morueta-Holme, N., Engemann, K., Sandoval-Acuña, P., Jonas, J. D., Segnitz, R. M., and Svenning, J.-C. (2015). Strong upslope shifts in Chimborazo’s vegetation over two centuries since Humboldt. Proc. Nat. Acad. Sci. 41, 12741–12745. doi: 10.1073/pnas.1509938112
Mountain Research Initiative EDW Working Group (2015). Elevation-dependent warming in mountain regions of the world. Nat. Clim. Chang. 5, 424–430. doi: 10.1038/nclimate2563
Myers, N., Mittermeier, R. A., Mittermeier, C. G., da Fonseca, G. A. B., and Kent, J. (2000). Biodiversity hotspots for conservation priorities. Nature 403, 853–858. doi: 10.1038/35002501
Navarro, G., Molina, J. A., and De la Barra, N. (2005). Classification of the high-Andean Polylepis forests in Bolivia. Plant Ecol. 176, 113–130. doi: 10.1007/s11258-004-0025-1
Ohlemüller, R., Gritti, E. S., Sykes, M. T., and Thomas, C. D. (2006). Quantifying components of risk for European woody species under climate change. Glob. Chang. Biol. 12, 1788–1799. doi: 10.1111/j.1365-2486.2006.01231.x
Parmesan, C. (2006). Ecological and evolutionary responses to recent climate change. Annu. Rev. Ecol. Evol. Syst. 37, 637–669. doi: 10.1146/annurev.ecolsys.37.091305.110100
Pérez-Escobar, O. A., Cámara-Leret, R., Antonelli, A., Bateman, R., Bellot, S., Chomicki, G., et al. (2018). Mining threatens Colombian ecosystems. Science 359, 1475–1475. doi: 10.1126/science.aat4849
Petrovan, S. O., Vale, C. G., and Sillero, N. (2020). Using citizen science in road surveys for large-scale amphibian monitoring: are biased data representative for species distribution? Biodiv. Conserv. 29, 1767–1781. doi: 10.1007/s10531-020-01956-0
Peyre, G., Balslev, H., and Font, X. (2018). Phytoregionalisation of the Andean Páramo. PeerJ 6:e4786. doi: 10.7717/peerj.4786
Peyre, G., Balslev, H., Font, X., and Tello, J. S. (2019). Fine-scale plant richness mapping of the andean páramo according to macroclimate. Front. Ecol. Evol. 7:12. doi: 10.3389/ffgc.2020.00012
Peyre, G., Lenoir, J., Karger, D. N., Gomez, M., Gonzalez, A., Broennimann, O., et al. (2020). The fate of Páramo plant assemblages in the sky islands of the northern Andes. J. Veg. Sci. 31, 967–980. doi: 10.1111/jvs.12898
Phillips, S. J., Anderson, R. P., and Schapire, R. E. (2006). Maximum entropy modeling of species geographic distributions. Ecol. Model. 190, 231–259. doi: 10.1016/j.ecolmodel.2005.03.026
Pinos, J. (2020). Challenges and conservation implications of Polylepis woodlands in the Andean region: Defining actions for sustainable management. Hacquetia 19, 143–153.
Purcell, J., and Brelsford, A. (2004). Reassessing the causes of decline of Polylepis, a tropical subalpine forest. Ecotropica 10, 155–158.
Rada, F., García-Núñez, C., and Rangel, S. (2011). Microclimate and regeneration patterns of Polylepis sericea in a treeline forest of the Venezuelan Andes. Ecotrópicos 2011, 113–122.
R Core Team (2020). R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing.
Ramirez-Villegas, J., Cuesta, F., Devenish, C., Peralvo, M., Jarvis, A., and Arnillas, C. A. (2014). Using species distributions models for designing conservation strategies of Tropical Andean biodiversity under climate change. J. Nat. Conserv. 22, 391–404. doi: 10.1016/j.jnc.2014.03.007
Ramos, C., Buitrago, S. P., Pulido, K. L., and Vanegas, L. J. (2013). Variabilidad ambiental y respuestas fisiológicas de Polylepis quadrijuga (Rosaceae) en un ambiente fragmentado en el Páramo de la Rusia (Colombia). Rev. Biol. Trop. 61, 351–361.
Rangel-Ch, J. O., and Arellano-P, H. (2010). “Bosques de Polylepis: un tipo de vegetación condenado a la extinción,” in Colombia Diversidad Biótica X. Cambio Global (Natural) y Climático (Antrópico) en el Páramo colombiano, ed. J. O. Rangel-Ch (Bogotá: Instituto de Ciencias Naturales-Facultad de Ciencias-Universidad Nacional de Colombia).
Raz, L., and Agudelo, H. (2019). Herbario Nacional Colombiano (COL). Version 13.12. Universidad Nacional de Colombia. Available online at: https://doi.org/10.15472/ea8sek (accessed May 03, 2020).
Rehm, E. M., and Feeley, K. J. (2015a). Freezing temperatures as a limit to forest recruitment above tropical Andean treelines. Ecology 96, 1856–1865. doi: 10.1890/14-1992.1
Rehm, E. M., and Feeley, K. J. (2015b). The inability of tropical cloud forest species to invade grasslands above treeline during climate change: potential explanations and consequences. Ecography 38, 1167–1175. doi: 10.1111/ecog.01050
Rehm, E. M., and Feeley, K. J. (2013). Forest patches and the upward migration of timberline in the southern Peruvian Andes. For. Ecol. Manage. 305, 204–211. doi: 10.1016/j.foreco.2013.05.041
Ripley, B. D. (1996). Pattern recognition and neural networks. New York, NY: Cambridge University Press.
Rodríguez Bolaños, A. (2020). Rasgos funcionales de plantas, anfibios y aves en el complejo de páramos Cruz Verde - Sumapaz. Instituto de Investigación de Recursos Biológicos Alexander von Humboldt. Available online at: https://doi.org/10.15472/dqelbs (accessed March 14, 2021).
Rodríguez, M. A., Angueyra, A., Cleef, A. M., and Van Andel, T. (2018). Ethnobotany of the Sierra Nevada del Cocuy-Güicán: climate change and conservation strategies in the Colombian Andes. J. Ethnobiol. Ethnom. 14, 1–12. doi: 10.1186/s13002-018-0227-6
Roy-Dufresne, E., Saltré, F., Cooke, B. D., Mellin, C., Mutze, G., Cox, T., et al. (2019). Modeling the distribution of a wide-ranging invasive species using the sampling efforts of expert and citizen scientists. Ecol. Evol. 9, 11053–11063. doi: 10.1002/ece3.5609
Różycka, M., Migoń, P., and Michniewicz, A. (2017). Topographic Wetness Index and Terrain Ruggedness Index in geomorphic characterization of landslide terrains, on examples from the Sudetes, SW Poland. Zeitschrift für Geomorphol. 61, 61–80. doi: 10.1127/zfg_suppl/2016/0328
Rull, V., and Vegas-Vilarrúbia, T. (2017). Potential responses of vascular plants from the pristine “Lost World” of the Neotropical Guayana highlands to global warming: Review and new perspectives. Front. Plant Sci. 8:81. doi: 10.3389/fpls.2017.00081
Rundel, P. W., Gibson, A. C., Midgley, G. S., Wand, S. J. E., Palma, B., Kleier, C., et al. (2003). Ecological and ecophysiological patterns in a pre-altiplano shrubland of the Andean Cordillera in northern Chile. Plant Ecol. 169, 179–193. doi: 10.1023/A:1026075721045
Sarmiento, C., Cadena, C., Sarmiento, M., Zapata, J., and León, O. (2013). Aportes a la conservación estratégica de los Páramos de Colombia: Actualización de la cartografía de los complejos de Páramo a escala 1:100.000. Bogotá, D.C: Instituto de Investigación de Recursos Biológicos Alexander von Humboldt.
Schmidt-Lebuhn, A. N., Kessler, M., and Kumar, M. (2006). Promiscuity in the Andes: Species relationships in Polylepis (Rosaceae, Sanguisorbeae) based on AFLP and morphology. Syst. Bot. 31, 547–559. doi: 10.1600/036364406778388629
Schmidt-Lebuhn, A., Seltmann, P., and Kessler, M. (2007). Consequences of the pollination system on genetic structure and patterns of species distribution in the Andean genus Polylepis (Rosaceae): a comparative study. Plant Syst. Evol. 266, 91–103. doi: 10.1007/s00606-007-0543-0
Schierenbeck, K. A. (2017). Population-level genetic variation and climate change in a biodiversity hotspot. Annals Bot. 119, 215–228. doi: 10.1093/aob/mcw214
Segovia-Salcedo, M. C., Domic, A., Boza, T. E., and Kessler, M. (2018). Situación taxonómica de las especies del género de Polylepis. Implicancias para los estudios ecológicos, la conservación y la restauración de sus bosques. Ecol. Austral. 28, 188–201. doi: 10.25260/EA.18.28.1.1.527
Sevillano-Ríos, C. S., and Rodewald, A. D. (2017). Avian community structure and habitat use of Polylepis forests along an elevation gradient. PeerJ. 5:e3220. doi: 10.7717/peerj.3220
Sevillano-Ríos, C. S., Rodewald, A. D., and Morales, L. V. (2018). Ecología y conservación de las aves asociadas con Polylepis: ¿qué sabemos de esta comunidad cada vez más vulnerable?. Ecología Austral 28, 216–228. doi: 10.25260/EA.18.28.1.1.519
SINAP (2020). Registro Único Nacional de Áreas Protegidas - RNUP. Available online at: http://www.parquesnacionales.gov.co/portal/es/sistema-nacional-de-areas-protegidas-sinap/registro-unico-nacional-de-areas-protegidas/ (accessed May 03, 2020).
Simpson, B. B. (1979). A revision of the genus Polylepis (Rosaceae: Sanguisorbeae). Smithson. Contrib. Bot. 1979, 1–62. doi: 10.5479/si.0081024X.43.1
Simpson, B. B. (1986). “Speciation and specialization of Polylepis in the Andes,” in High altitude tropical biogeography, eds F. Vuilleumier and M. Monasterio (New York, NY: Oxford University Press).
Sousa-Silva, R., Alves, P., Honrado, J., and Lomba, A. (2014). Improving the assessment and reporting on rare and endangered species through species distribution models. Glob. Ecol. Conserv. 2, 226–237. doi: 10.1016/j.gecco.2014.09.011
Tejedor Garavito, N., Newton, A. C., Golicher, D., and Oldfield, S. (2015). The relative impact of climate change on the extinction risk of tree species in the montane tropical andes. PLoS One 10:e0131388. doi: 10.1371/journal.pone.0131388
Thuiller, W., Georges, D., Engler, R., and Breiner, F. (2020). BIOMOD2: Ensemble platform for species distribution modeling. R package version 3.4.6.
Thuiller, W., Lafourcade, B., Engler, R., and Araújo, M. B. (2009). BIOMOD - a platform for ensemble forecasting of species distributions. Ecography 32, 369–373. doi: 10.1111/j.1600-0587.2008.05742.x
Thuiller, W., Lavorel, S., and Araujo, M. (2005). Niche properties and geographical extent as predictors of species sensitivity to climate change. Glob. Ecol. Biogeo. 14, 347–357. doi: 10.1111/j.1466-822X.2005.00162.x
Tiago, P., Pereira, H. M., and Capinha, C. (2017). Using citizen science data to estimate climatic niches and species distributions. Basic Appl. Ecol. 20, 75–85. doi: 10.1016/j.baae.2017.04.001
Title, P. O., and Bemmels, J. B. (2018). ENVIREM: An expanded set of bioclimatic and topographic variables increases flexibility and improves performance of ecological niche modeling. Ecography 41, 291–307. doi: 10.1111/ecog.02880
Toivonen, J. M., Gonzales-Inca, C. A., Bader, M. Y., Ruokolainen, K., and Kessler, M. (2018). Elevational shifts in the topographic position of Polylepis forest stands in the Andes of Southern Peru. Forests 9:7. doi: 10.3390/f9010007
Toivonen, J. M., Horna, V., Kessler, M., Ruokolainen, K., and Hertel, D. (2014). Interspecific variation in functional traits in relation to species climatic niche optima in Andean Polylepis (Rosaceae) tree species: evidence for climatic adaptations. Funct. Plant Biol. 41, 301–312. doi: 10.1071/FP13210
Tovar, C., Arnillas, C. A., Cuesta, F., and Buytaert, W. (2013). Diverging responses of tropical andean biomes under future climate conditions. PLoS one 8:e63634. doi: 10.1371/journal.pone.0063634
Ueda, K. (2021). iNaturalist Research-grade Observations. iNaturalist.org. Available online at: https://doi.org/10.15468/ab3s5x (accessed March 14, 2021).
Valencia, J. B., Mesa, J., León, J. G., Madriñán, S., and Cortés, A. J. (2020). Climate vulnerability assessment of the Espeletia complex on Páramo sky islands in the northern Andes. Front. Ecol. Evol. 8:309. doi: 10.3389/fevo.2020.565708
van Schaik, C. P., Terborgh, J. W., and Wright, S. J. (1993). The phenology of tropical forests: Adaptive significance and consequences for primary consumers. Annu. Rev. Ecol. Syst. 24, 353–377. doi: 10.1146/annurev.es.24.110193.002033
Vásquez, E., Ladd, B., and Borchard, N. (2014). Carbon storage in a high-altitude Polylepis woodland in the Peruvian Andes. Alp. Bot. 124, 71–75. doi: 10.1007/s00035-014-0126-y
Velez, V., Cavelier, J., and Devia, B. (1998). Ecological traits of the tropical treeline species Polylepis quadrijuga (Rosaceae) in the Andes of Colombia. J. Trop. Ecol. 14, 771–787. doi: 10.1017/S026646749800056X
Vuille, M., Franquist, E., Garreaud, R., Lavado Casimiro, W. S., and Cáceres, B. (2015). Impact of the global warming hiatus on Andean temperature: Global warming hiatus in the Andes. J. Geo Res: Atmospheres 120, 3745–3757. doi: 10.1002/2015JD023126
Vuille, M., Francou, B., Wagnon, P., Juen, I., Kaser, G., Mark, B., et al. (2008). Climate change and tropical Andean glaciers: Past, present and future. Earth-Sc. Rev. 89, 79–96. doi: 10.1016/j.earscirev.2008.04.002
Wesche, K., Cierjacks, A., Assefa, Y., Wagner, S., Fetene, M., and Hensen, I. (2008). Recruitment of trees at tropical alpine treelines: Erica in Africa versus Polylepis in South America. Plant Ecol. Div. 1, 35–46. doi: 10.1080/17550870802262166
Wisz, M. S., Hijmans, R. J., Li, J., Peterson, A. T., Graham, C. H., Guisan, A., et al. (2008). Effects of sample size on the performance of species distribution models. Div. Distrib. 14, 763–773. doi: 10.1111/j.1472-4642.2008.00482.x
Zizka, A., Antonelli, A., and Silvestro, D. (2020). sampbias, a method for quantifying geographic sampling biases in species distribution data. Ecography 44, 25–32. doi: 10.1111/ecog.05102
Zutta, B. R., and Rundel, P. W. (2017). Modeled shifts in Polylepis species ranges in the andes from the last glacial maximum to the present. Forests 8:232. doi: 10.3390/f8070232
Keywords: climate change, citizen science data, habitat suitability modeling, high-Andean forests, Páramo, street view imagery
Citation: Caballero-Villalobos L, Fajardo-Gutiérrez F, Calbi M and Silva-Arias GA (2021) Climate Change Can Drive a Significant Loss of Suitable Habitat for Polylepis quadrijuga, a Treeline Species in the Sky Islands of the Northern Andes. Front. Ecol. Evol. 9:661550. doi: 10.3389/fevo.2021.661550
Received: 31 January 2021; Accepted: 22 April 2021;
Published: 09 June 2021.
Edited by:
Giuseppe Fenu, University of Cagliari, ItalyReviewed by:
Kenneth J. Feeley, University of Miami, United StatesEstrella Alfaro Saiz, Universidad de León, Spain
Copyright © 2021 Caballero-Villalobos, Fajardo-Gutiérrez, Calbi and Silva-Arias. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Gustavo A. Silva-Arias, gustavo.silva@tum.de