Elevation and latitude drives structure and tree species composition in Andean forests: Results from a large-scale plot network

Our knowledge about the structure and function of Andean forests at regional scales remains limited. Current initiatives to study forests over continental or global scales still have important geographical gaps, particularly in regions such as the tropical and subtropical Andes. In this study, we assessed patterns of structure and tree species diversity along ~ 4000 km of latitude and ~ 4000 m of elevation range in Andean forests. We used the Andean Forest Network (Red de Bosques Andinos, https://redbosques.condesan.org/) database which, at present, includes 491 forest plots (totaling 156.3 ha, ranging from 0.01 to 6 ha) representing a total of 86,964 identified tree stems ≥ 10 cm diameter at breast height belonging to 2341 identified species, 584 genera and 133 botanical families. Tree stem density and basal area increases with elevation while species richness decreases. Stem density and species richness both decrease with latitude. Subtropical forests have distinct tree species composition compared to those in the tropical region. In addition, floristic similarity of subtropical plots is between 13 to 16% while similarity between tropical forest plots is between 3% to 9%. Overall, plots ~ 0.5-ha or larger may be preferred for describing patterns at regional scales in order to avoid plot size effects. We highlight the need to promote collaboration and capacity building among researchers in the Andean region (i.e., South-South cooperation) in order to generate and synthesize information at regional scale.


Introduction
The Tropical Andes have been identified as one of the most important hotspots of global biodiversity [1,2]. Indeed, this region is one of the most diverse terrestrial hotspots on Earth [1]. There are approximately 45,000 identified vascular plant species occurring in the Andes, 20,000 of which are endemic to the region [1]. There are also likely to be many more species (up to 35%) that have not yet been described [3]. Moreover around 60 million people depend directly or indirectly on ecosystem services provided by Andean forests, such as water and the regulation of regional climates [4,5,6,7,8].
Andean forests face a high risk of degradation as a result of climate change [9,10] and land use change [11,12] due to human population growth [7] and migration [12,13,14,15], and a combination of these factors [15]. Changes in Andean forest cover includes both forest regrowth mostly frequent above 2000 m asl and deforestation that often concentrates below 1000 m asl [12]. Some Andean countries have already lost~50-60% (~550.500 km 2 total) of their cloud forests [16] due to historical and ongoing land use changes, linked to the expansion of the agricultural frontier, increase in extension of cattle ranching areas for pastures, and other economic activities such as mining, road constructions among others [17,18,19]. Furthermore, as a result of warming, some forests and grasslands are experiencing an apparent shift in plant species composition (i.e. increasing relative abundances of species from lower elevations) [10,20]. Moreover, model projections for biomes [21] and vascular plant species [22] predict upward shifts and the loss of optimal conditions in their lower and mid ranges. These changes may compromise the persistence of Andean forest ecosystems and may reduce the provision of benefits for human populations.
Recently, Mathez-Stiefel et al. [23] developed a research agenda for Andean forests landscapes where they highlighted critical knowledge gaps on the ecology of Andean mountain ecosystems, and the effect of socio-environmental changes over forest functions and dynamics. This includes the need to understand general patterns of tree diversity along elevational and latitudinal gradients, as well as forest structure and forest dynamics across environmental gradients at continental scales throughout integrated and collaborative research initiatives [23].
Tree or stem density and basal area may increase with elevation in extra Andean tropical montane forests [24,25,26,27,28]. However, in Andean forests, our understanding of structural parameters is limited to a handful of locations that have reported differing trends [29,30,31]. For example, in the tropical section of the Andes, basal area may increase with altitude until mid-elevations [32], and then decrease [30] or remain stable [29]. In subtropical Andes, basal area tended to increase with elevation [31]. Baez et al. [33] made an effort to describe regional patterns of Andean forests structure and dynamics along elevation and latitudinal gradients. These authors found that in tropical mountain forests, basal area increased at lower elevation but did not change along the elevation gradient. They also found that in subtropical Andes, changes in basal area appear to be influenced by land use history together with environmental variation [33].
Plant species richness and diversity are thought to decline away from the equator and towards higher latitudes and elevations [34,35]. However, elevational patterns of species richness have been studied much more intensively than latitudinal ones [34]. For example, a hump-shaped diversity pattern has been reported for various plant groups along elevational gradients in the Andes, including vascular epiphytes in Bolivia [36], non-vascular bryophytes in Colombia [37], and ferns in Ecuador and Bolivia [35]. Tree species richness (as well as tree genera and family richness) also shows the hump-shaped pattern, reaching a maximum at intermediate elevations (i.e. 1000 m asl) in subtropical Andes of Argentina [38], while decreasing above 1800 m asl in tropical Andes of Ecuador [30]. Still, regional-scale studies of tree community structure and composition along broader environmental gradients in the Andes are lacking.
This study is one of the first attempts to describe and characterize regional patterns of forest structure and diversity in Andean forest communities across~4000 km of latitude and~4000 m of elevation gradient. In this sense, this study incorporates a higher dataset than Baez et al. [33] (i.e. 8 times more) as well as diversity patterns. Specifically, for structural patterns we assessed: stem density and basal area, and for diversity patterns we assessed: species richness, tree species composition and floristic similarity along these elevation and latitudinal gradients. To accomplish this, we used the Andean Forest Network, AFN (Red de Bosques Andinos, https://redbosques.condesan.org/) database that provides data on forests throughout the tropical and subtropical Andes. Our main hypothesis is that the climatic gradient generated by elevation and latitude model the structural and diversity patterns of Andean forests, with plot size playing a role in explaining these trends. Specifically, we hypothesis that: (1) Tree communities located near the equator will have higher species richness, higher stem density and higher basal area than those located further away from the equator; and (2) Tree communities from lower elevations have a higher species richness, lower stem density and lower basal area than those located at higher elevations.

Andes origin and vegetation
The Andes Mountains originated through a major uplift of the Central Andes during the Paleogene (65 to 34 Ma) and subsequent plates collisions intensified mountain building in the Northern Andes (23 Ma) [39]. However, the Andes reached their modern elevation during the late mid Miocene (~12 Ma) and early Pliocene (~4.5 Ma) [39,40]. The formation of the Andes mountains had an enormous effect on regional climate regulation and was crucial for the evolution of landscapes and ecosystems in South America, including current biodiversity patterns of Amazon ecosystems [39]. For example, the uplift of the cordillera changed the Amazonian landscape by re-configuring the patterns of drainage and creating a vast influx of sediments into the basin [39]. Recent Pleistocene glaciations (110,000 to 10,000 b.p.) were also important in determining Andean forest biodiversity and species distributions [41]. For example, by acting as a refuge for biodiversity and by allowing immigration of holarctic species (e.g. Alnus), not only due to changes in temperature and CO 2 worldwide, but also due to variations in precipitation over the Amazon which had direct influence in Andean forests [42].
The Andean mountain forest ecosystem zonation is expressed in changes of forest architecture (i.e. decreases in tree stature and stem diameter, trends in stem deformation, hard, thick and smaller leaves; Fig 1) and tree community composition [16,29,30,32,33,43] as elevation increases. The tropical and subtropical slopes vary in floristic composition and forest structure, from premontane forests (also called sub-andean forests) at lower elevations (from 500/800 to 700/1500 m asl) toward lower montane forests (1500/2700 m asl), upper montane forests (2700/2800 to 3300 m asl), and finally the upper treeline forests at higher elevations (3300/ 3500 -exceptionally 4000 m asl (Fig 1) [44]. Elevation ranges and vegetation zonation vary according to latitude, environmental humidity and topography [44]. For example, in the Andes of Colombia and Ecuador, moist inter-Andean valley slopes present similar gradations in vegetation zones. Further, vegetation zonation is related to changes in the frequency and persistence of cloud formation, which affects solar radiation, air temperature, and water regimes. Cloud formation varies from low cloud cover in premontane forest towards more persistent cloudiness in upper montane forest [16,45,46].

Andean Forest Network
Background and standardized protocol of measurement. The AFN was created in 2012, being the first Network at this regional scale, and brings together scientists and decision makers interested in research, management and conservation of Andean forests. The main objective of the AFN is to generate knowledge about the ecology of the Andean forests throughout the collaborative work of its members, exchange, systematization and synthesis of information, development of research protocols, strengthening of research capabilities, and the articulation with decision making processes in the region. Specifically, the Andean Forest Network aims to: 1) promote the establishment and maintenance of long-term research sites through the installation and monitoring of permanent vegetation plots; 2) enhance the ability to understand the structure, composition, dynamics and functioning of Andean forests in a global context of environmental change; 3) promote collaboration and capacity building among researchers and technicians in the Andean region (i.e. South-South cooperation). For general information about the Network see https://redbosques.condesan.org/. Members of the Andean Forest Network have detailed forest tree data (i.e. species identity, diameter, spatial coordinates) derived from research plots (both permanent and non-permanent). The Network prioritizes collaborative efforts amongst network members, encouraging them to share their data under confidence and transparent concepts, with previous consent and approval of the principal researchers and/or institutions [i.e. 10,33,47,48,49].
An important contribution of the Andean Forest Network has been the development of protocols for the establishment of permanent monitoring plots of plant diversity adapted to the conditions prevalent in the Andean forests [47,48] which provide methodological approaches to study forest structure and dynamics of tree communities. Osinaga Acosta et al. [47], also includes an approach to consider and measure lianas, palms, ferns and herbs. For trees (with diameter �10 cm at breast height, DBH), the protocol encourages the establishment of 1-ha permanent forest plots subdivided in 20 × 20 m quadrants, where all individuals are identified to the lowest taxonomic level, marked with a numbered tags, and mapped in a x-y coordinate system to 1-m resolution. Botanical samples are collected and incorporated into local institutional herbaria. Plots are periodically recensused during which time mortality and recruitment events are recorded and the DBH of marked trees is remeasured. Recensus intervals vary among research objectives and available resources, although a five-year interval is suggested in order to maintain consistent monitoring over time [47].
Forest plot distribution, general information and climate data. At present, the Andean Forest Network comprises 491 forest plots spanning a latitudinal range from tropical Venezuela (8.63˚N) to subtropical Argentina (26.77˚S), a longitudinal range from the Pacific to the Atlantic versants (-80.14 to -63.84˚W), and an elevation range from 41 to 3980 (mean ± SD = 1901 ± 903) m asl (Fig 2). Forest plots setup were based on the combination of several factors: (1) accessibility, (2) sustainable over time to ensure long term monitoring for some plots, and (3) homogeneity in terms of forest type, topography, and disturbances in order to be representative of particular forest types. The majority of the plots are located in mature forests but the AFN also includes several older secondary forests plots, i.e. > 30 years old [50] as a result of human activities and the associated land use changes [19,51,52,53].  Table 1, S1 Table). Considering DBH cut off points, 47% of forest plots (n = 229) measured trees �10 cm, 37% (n = 184) measured trees �5 cm, 12% of forest plots (n = 60) measured trees �2.5 cm, and 4% (n = 18) measured trees �1cm (S1 Table). All relevant data upon which all the presented results in this manuscript are based on is included in the Supporting Information files; number of stems, basal area (both extrapolated to 1-ha), and species richness for the 491 plots in S1 Table while the abundance of species per country in S1 Appendix.
We extracted the following climatic data for each of the 491 plots from the CHELSA (Climatologies at high resolution for the earth's land surface areas) dataset at 30 arc sec (~1 km) resolution (period 1970-2013) [54]: estimated ranges in total annual rainfall varied from 608 to 5059 mm yr -1 while estimated ranges in mean annual temperature varied from 5.9 to 27˚C ( Table 1, S1 Table). Estimated mean temperature for the three coldest months varied between -1 and 24˚C, while the estimated mean temperature for the three warmest months varied from 11 to 32˚C. Predicted seasonal thermal amplitude (mean temperature from the warmest quarter -mean temperature from the coldest quarter) varied from 1˚C in Ecuador, Colombia and Venezuela to 10˚C in Argentina, showing intermediates values in Peru (2˚C) and Bolivia (4˚C), depicting a clear latitudinal trend in thermal seasonality such that the thermal amplitude increases with distance from the equator.

Data analyses
To analyze patterns in stem density, basal area and species richness, we considered those plots � 0.1-ha (n = 236) where stems �10 cm DBH were measured. Due to high variability and dispersion when extrapolating basal area and stem density at a larger spatial scale we discarded plots of 0.01-ha, 0.04-ha, 0.05-ha and 0.08-ha (total n = 255). For those plots � 0.1-ha, we estimated stem density and basal area at 1-ha for comparison while species richness was evaluated on a per plot basis. For calculations of stem density and basal area we used all stems �10 cm DBH (n = 97,054) while for species richness we used stems �10 cm DBH identified to species level (n = 86,964). We performed generalized linear models (GLM) [55] and used elevation, absolute latitude, and plot size as explanatory variables. We included the quadratic terms of elevation and latitude in the models in order to describe unimodal trends. We used a Quasi-Poisson distribution for stem density and species richness since Poisson distribution showed overdispersion in both variables, and a Gaussian distribution with log-link function for basal area [55]. We calculated the percentage of variance explained with the following equation VE = (Null Deviance-Residual Deviance)/Null Deviance x 100.
To analyze patterns of species composition and floristic similarity, we used all data available from all plots (n = 491) and from all species fully identified (2341) with stems �10 cm DBH. To describe species composition, we applied a Non-Metric Multidimensional Scaling approach (NMDS) [56] based on a Bray-Curtis distance matrix [57] calculated from species abundance (stems per plots) between pairs of plots. We used a two-dimensional configuration because the final stress, an index of agreement between the distances in the graph configuration and the distances in the Bray-Curtis matrix, was 9.2 (most ecological community data sets have solutions with stress between 10 and 20; [58]. To explore the relationships between plots based on their species composition along elevation and latitudinal gradients we used Spearman's correlation coefficients [59] between the scores in the axes of the NMDS and elevation and latitude. To describe floristic similarity, we estimated Bray Curtis distance based on species abundance data [60]. Bray Curtis distances were calculated among pairs of plots within latitudinal bands. For this, we created eight latitudinal bands every 5˚along the latitudinal gradient (from 8˚N to 27˚S) ( Table 2). Floristic similarity (i.e. beta diversity) varies between 0 and 100%, i.e. values close to 0 imply less shared species. To explore the relationship between Bray Curtis distance with the geographical distance (km) and elevational difference (m asl) between pairs of plots we performed a linear mixed-effects model (LMM) [61] using latitudinal bands as a random factor. Geographical distances and elevational differences were obtained with Euclidean distances among pairs of plots within each latitudinal band, considering geographical coordinates and elevation, respectively. We calculated the proportion of variance explained within and among latitudinal bands, i.e., marginal and conditional R 2 , respectively, according to the definitions given by Nakagawa & Schielzeth [62]. All analyses were performed in R 3.4.3 [63], using AER to test overdispersion in GLM, vegan for ordination analysis, and lme4 for LMM.

Results
The 491 plots were distributed along a 4138 km latitudinal gradient, encompassing a wide range of climates (Fig 3). Elevation ranges decreased from tropical to subtropical latitudes depending on the treeline location ( Fig 3A). Mean annual air temperature and rainfall were Table 2 positively correlated (r = 0.46, p < 0.001) ( Fig 3B). As expected, mean annual temperature decreased with elevation (r = -0.95, p < 0.001) but did not vary with latitude (r = -0.06, p = 0.31) (Fig 3C). However, seasonal thermal amplitude increased with distance from the equator (r = -0.93, p < 0.001). Rainfall tended to decrease with elevation (r = -0.38, p < 0.001) and decreased towards higher latitudes (r = 0.45, p < 0.001) (Fig 3D). We found that stem density tended to increase with elevation but decreased above 3000 m asl (S1 Appendix, Fig 4A). Stem density showed a unimodal trend with latitude, peaking at 10-15 degrees from the equator (Fig 4B). Elevation and latitude (and their quadratics terms) explained 43% of the variation in stem density, while plot size did not have a significant effect on stem density (S1 Appendix). Basal area tended to increase with elevation but peaked around 2000 m asl (Fig 4C). Basal area did not relate with latitude ( Fig 4D). Plot size had a significant effect on basal area, small plots of 0.1-ha tended to present higher basal area than 1-ha plots (S1 Appendix). Combined, elevation, latitude and plot size explained 9.5% of the variation in basal area. Species richness showed a unimodal trend with elevation, peaking at 1000 m asl ( Fig 4E) and decreasing with latitude being higher in the tropics (Fig 4F). Plot size also had a significant effect on species richness such that richness increased with the area sampled. Elevation, latitude and plot size together explained 70% of the variation in species richness. Taking into account the effect of plot size, we fitted a model only with 1-ha plots (n = 109) that showed, on average, higher species richness in tropical forest plots (Intercept = 68 tree species ha -1 , above 15˚S) than in subtropical forest plots of Argentina (Intercept = 27 tree species ha -1 , between 22˚-27˚S) (Fig 4G, S1 Appendix).

. Summary of total species identified per latitudinal band, followed by mean Bray-Curtis distance (max value in brackets), mean geographical distance (max value in brackets), number of plots and mean plot size in hectares (range values in brackets) across the latitudinal gradient of the Andean Forest
Considering all plots (n = 491), we registered a total of 97,054 tree stems' � 10 cm DBH, of which 90% (86,964) were identified to species level belonging to 2341 species, 584 genera and 133 botanical families (S2 Appendix), while 1095 remained unidentified as morphospecies of which 1019 remain to genera, 41 to family level, and 35 were not determined. Considering identified species (n = 2341), the richest genera were Miconia (105 species), followed by Inga The study plots tended to group along two NMDS dimensions based on their species composition (Fig 5). Elevation and latitude correlated well with both NMDS axes (Axis 1: elevation Bray-Curtis distance (i.e., floristic similarity) across all Andean plots was on average 4% (the lower the value the less similarity among plots). Latitudinal bands between 10˚N and 15S (i.e. tropical montane forest plots) showed mean values between 3 to 9% in floristic similarity while in latitudinal bands below 20˚S (i.e. subtropical montane forest plots), floristic similarity presented mean values of 13 and 16% (Table 2). These latitudinal bands explained 30.7% of floristic similarity variation among latitudinal bands (i.e., tropical vs. subtropical forest plots) while geographical and elevational distance between plots explained only 2% of floristic similarity decay within latitudinal bands (Fig 6A and 6B; S3 Appendix). Together, geographical distance and elevational difference between pairs of plots within and among latitudinal bands explained 32.8%.

Discussion and conclusions
The 491 forest plots considered in this study covered broad range of latitudes, elevations and environmental conditions. As a general pattern, temperature and rainfall decreased with elevation, implying that forest plots located at lower elevations were subject to warmer and wetter conditions. While this pattern is well established for temperature, it was different from what Urrutia et al. [5] reported for the Andes where rainfall did not co-vary linearly with elevation. In the latitudinal gradient, rainfall increased towards the equator while mean temperature did not show a clear pattern but seasonal thermal amplitude was higher with increasing latitude (i.e. further from the equator) [64].
Both stem density and basal area were related with the gradients covered by our dataset. However, patterns of stem density were more consistent with elevation and latitude than for basal area. Forest plots with higher stem density (and basal area) were located at higher elevations (> 2000 m asl) and towards the equator. Stem density tended to peaked at 10-15˚latitude (e.g. plots of Colombia). The observed increase in the mean values of structural variables with elevation has not been frequently reported; on the contrary, some authors reported the opposite pattern of decreasing density [32], basal area or productivity with elevation [30,33]. Probably, there are other variables masked in our study, such as the latitudinal gradient (because the highest densest plots are in the tropics) or maybe some local climate associated to topographic conditions (i.e. wetter or drier slopes, soil properties). In addition land use change history may influence tree density as well in the case of secondary forests, i.e. if associated with tree invasive species [65]. Furthermore, the use of small plots (i.e. 0.1-ha) which represent about 12% of plots in AFN, tended to overestimate basal area in comparison to larger plots (mean basal area per hectare 33.9± 1.9 for 0.1-ha plots; and 26.5± 0.7 for 1-ha plots). Moreover, the basal area of these small plots presented the opposite pattern to the general trend, i.e. decreasing basal area with elevation, being consistent with other studies that use small plots [30,33]. A stronger edge effect and differences in disturbance intensities (disturbed vs. mature forests) as well as majestic forest effect could be expected in these small plots. In addition, smaller plots (e.g. between 0.01 and 0.05-ha), which were discarded from our main analyses, showed a very large variability and dispersion when scaled to 1-ha (e.g. 0.01 ha plot with 1 stem would extrapolated to 100 stems ha -1 , and in a 0.04-ha plot with 35 trees extrapolated value of basal area would be as high as 131.2 m 2 ha -1 ). Overall, plots~0.5-ha or larger may be preferred for describing patterns at regional scales.
The greater species richness at lower latitudes compared to the higher latitudes has been widely studied. Explanations for this pattern include the Janzen and Connell conspecific density dependent hypothesis which maintain diversity in plant communities by reducing survival rates of conspecific seedlings located close to reproductive adults or in areas of high conspecific density [66], the narrow environmental niche [67], the recruitment limitation hypothesis [68]. Also, the tropical niche conservatism hypothesis may explain the observed pattern where species diversity and biogeographic processes are driven by historical climate where the stable, wet, warm and a-seasonal tropical climate promote high diversity [69,70,71] and narrower niches (i.e. the latitudinal gradient hypothesis) [64,72].
Subtropical forest plots (i.e., plots in Argentina) clearly segregated from the rest of the tropical forest plots based on species composition, although the elevation gradient was also important, particularly in tropical plots. In terms of floristic similarity (i.e., beta diversity or species turnover) tropical forest plots showed between 3% to 9% of species shared implying high species turnover among their plots within latitudinal bands, while subtropical plots shared more species, between 13% and 16%, implying lower species turnover among plots associated with low diversity of tree species. The high species turnover found in the tropics was reported by several authors and attributed to the combination of climatic variables [73], to multiple coexistence mechanisms [74], and to the high biotic interaction found at these lower latitudes [75] as well as high habitat heterogeneity driven by ample environmental gradients [43,76].
As well, tropical forest plots separated by larger geographical distances (km) may be less similar due the presence of intermontane valleys and the influence of different floristic regions. This pattern (i.e. less floristic similarity or high species turnover found in the tropics) has also been reported for Andean high alpine vegetation [77]. As the variation explained among latitudinal bands in this regional context was very prominent, geographical and elevational distances between forest plots seemed to be relatively less important. Considering the elevation gradient of the Andes, our study suggests a mid-elevation peak in species diversity for both tropics and subtropics which respond to the mid domain effect hypotheses [78] related to mountain topographic constraints as a result of low temperatures although the mechanisms involved are not clear [36]. Our findings may imply that forests located in tropical and towards mid elevations sections of the Andes have the potential to accumulate more biomass, and consequently sequestrate more carbon, than forests located away from the equator or towards lower elevations. Conservation efforts may be particularly important within this forests.
The creation of the Andean Forest Network has served as an important platform for communication where researchers were able to discuss the use and interpretation of forest plot information in a context of global environmental change. The Andean Forest Network expects to continue expanding its collaborative work across the region adding forest plots along the Andes. Some of the key issues that emerged from its functioning includes: 1) strengthening the institutional arrangements within the Andean Forest Network, and between it and other actors interested in conducting research in Andean forests; 2) implementing mechanisms to promote the long-term sustainability of the Network and the monitoring activities of its members. This will require identifying funding opportunities within the region and abroad, from diverse sources such as climate funding, grants for field research, public funding for science, among others. Also, the compilation and management of the large database includes the need of: 1) improving species identifications and, when possible, standardizing morphospecies along plots; 2) developing a "decentralized" database or an online platform to share data quickly and easily with the least possible errors; 3) strengthening the on-site collection of climate data from meteorological stations; 4) filling current geographical gaps (for example, Bolivian montane dry forest or central Peruvian yungas) and encouraging institutions/researchers to contribute with their forests plot data to the Network, increasing the area covered per country and promoting re-measurement of the installed plots; 5) promoting the implementation of other modules available in protocols (e.g. carbon stocks, species functional traits) [47] and other useful complementary measurements (e.g. dendrometer bands) in all the permanent forest plots of the Network; 6) detecting situations that are locally important; 7) promoting comparative research for monitoring the impacts of anthropogenic activities across the region, including the dynamics of invasive species (i.e. the case of Ligustrum lucidum in Argentina) [79], or primary/secondary succession as a result of volcanism, forest fires, hurricanes, or land use change (abandonment of agricultural and livestock lands); 8) exploring in detail the association between changes in permanent plots and changes in remotely sensed descriptors of functioning (e.g., NDVI) [80], promoting joint research to develop high resolution models of climate change for the Andean region, and descriptions of land use change [12].
Overall, this study is one of the first attempts to integrate forest structure and composition in the Andes showing clear patterns along its elevation and latitudinal gradients. In this sense, the AFP Network constitutes an important initiative to fill geographical gaps in regions such as the tropical and subtropical Andes. .
Supporting information S1  plots, the Instituto Nacional de Parques in Venezuela for facilitating establishment of permanent plots in the Sierra Nevada National Park. We thank Andrea Izquierdo, Javier Foguet and Silvia Pacheco for their contribution with climatic data and collaboration with the map of Fig 2.