Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Large-Scale Patterns of Turnover and Basal Area Change in Andean Forests

  • Selene Báez ,

    Contributed equally to this work with: Selene Báez, Agustina Malizia, Julieta Carilla, Álvaro Duque, Jürgen Homeier, Oliver L. Phillips, Kenneth J. Feeley

    selene.baez@condesan.org

    Affiliations Consorcio para el Desarrollo Sostenible de la Ecoregión Andina (CONDESAN), Quito, Ecuador, Universidad Técnica Particular de Loja, Loja, Ecuador

  • Agustina Malizia ,

    Contributed equally to this work with: Selene Báez, Agustina Malizia, Julieta Carilla, Álvaro Duque, Jürgen Homeier, Oliver L. Phillips, Kenneth J. Feeley

    Affiliations Consejo Nacional de Ciencias de Investigaciones Científicas y Técnicas (CONICET), Buenos Aires, Argentina, Instituto de Ecología Regional (IER), Universidad Nacional de Tucumán, Tucumán, Argentina

  • Julieta Carilla ,

    Contributed equally to this work with: Selene Báez, Agustina Malizia, Julieta Carilla, Álvaro Duque, Jürgen Homeier, Oliver L. Phillips, Kenneth J. Feeley

    Affiliations Consejo Nacional de Ciencias de Investigaciones Científicas y Técnicas (CONICET), Buenos Aires, Argentina, Instituto de Ecología Regional (IER), Universidad Nacional de Tucumán, Tucumán, Argentina

  • Cecilia Blundo,

    Affiliations Consejo Nacional de Ciencias de Investigaciones Científicas y Técnicas (CONICET), Buenos Aires, Argentina, Instituto de Ecología Regional (IER), Universidad Nacional de Tucumán, Tucumán, Argentina

  • Manuel Aguilar ,

    ‡ These authors also contributed equally to this work.

    Affiliation Facultad de Ciencias Forestales, Universidad Nacional Agraria La Molina, Lima, Perú

  • Nikolay Aguirre ,

    ‡ These authors also contributed equally to this work.

    Affiliation Universidad Nacional de Loja, Loja, Ecuador

  • Zhofre Aquirre ,

    ‡ These authors also contributed equally to this work.

    Affiliation Universidad Nacional de Loja, Loja, Ecuador

  • Esteban Álvarez ,

    ‡ These authors also contributed equally to this work.

    Affiliation Laboratorio de Servicios Ecosistémicos y Cambio Climático, Jardín Botánico de Medellín, Medellín, Colombia

  • Francisco Cuesta ,

    ‡ These authors also contributed equally to this work.

    Affiliation Consorcio para el Desarrollo Sostenible de la Ecoregión Andina (CONDESAN), Quito, Ecuador

  • Álvaro Duque ,

    Contributed equally to this work with: Selene Báez, Agustina Malizia, Julieta Carilla, Álvaro Duque, Jürgen Homeier, Oliver L. Phillips, Kenneth J. Feeley

    Affiliation Departamento de Ciencias Forestales, Universidad Nacional de Colombia, Medellín, Colombia

  • William Farfán-Ríos ,

    ‡ These authors also contributed equally to this work.

    Affiliation Department of Biology, Wake Forest University, Winston-Salem, North Carolina, United States of America

  • Karina García-Cabrera ,

    ‡ These authors also contributed equally to this work.

    Affiliation Department of Biology, Wake Forest University, Winston-Salem, North Carolina, United States of America

  • Ricardo Grau ,

    ‡ These authors also contributed equally to this work.

    Affiliation Consejo Nacional de Ciencias de Investigaciones Científicas y Técnicas (CONICET), Buenos Aires, Argentina

  • Jürgen Homeier ,

    Contributed equally to this work with: Selene Báez, Agustina Malizia, Julieta Carilla, Álvaro Duque, Jürgen Homeier, Oliver L. Phillips, Kenneth J. Feeley

    Affiliation Plant Ecology, University of Göttingen, Göttingen, Germany

  • Reynaldo Linares-Palomino ,

    ‡ These authors also contributed equally to this work.

    Affiliation Facultad de Ciencias Forestales, Universidad Nacional Agraria La Molina, Lima, Perú

  • Lucio R. Malizia ,

    ‡ These authors also contributed equally to this work.

    Affiliations Facultad de Ciencias Agrarias, Universidad Nacional de Jujuy, Jujuy, Argentina, Fundación ProYungas, Jujuy, Argentina

  • Omar Melo Cruz ,

    ‡ These authors also contributed equally to this work.

    Affiliation Grupo de investigación en Biodiversidad y Dinámica de Ecosistemas Tropicales, Universidad del Tolima, Bogotá, Colombia

  • Oriana Osinaga ,

    ‡ These authors also contributed equally to this work.

    Affiliation Instituto de Ecología Regional (IER), Universidad Nacional de Tucumán, Tucumán, Argentina

  • Oliver L. Phillips ,

    Contributed equally to this work with: Selene Báez, Agustina Malizia, Julieta Carilla, Álvaro Duque, Jürgen Homeier, Oliver L. Phillips, Kenneth J. Feeley

    Affiliation School of Geography, University of Leeds, Leeds, United Kingdom

  • Carlos Reynel ,

    ‡ These authors also contributed equally to this work.

    Affiliation Facultad de Ciencias Forestales, Universidad Nacional Agraria La Molina, Lima, Perú

  • Miles R. Silman ,

    ‡ These authors also contributed equally to this work.

    Affiliation Department of Biology, Wake Forest University, Winston-Salem, North Carolina, United States of America

  •  [ ... ],
  • Kenneth J. Feeley

    Contributed equally to this work with: Selene Báez, Agustina Malizia, Julieta Carilla, Álvaro Duque, Jürgen Homeier, Oliver L. Phillips, Kenneth J. Feeley

    Affiliation International Center for Tropical Botany, Department of Biological Sciences, Florida International University, Miami, Florida, United States of America

  • [ view all ]
  • [ view less ]

Abstract

General patterns of forest dynamics and productivity in the Andes Mountains are poorly characterized. Here we present the first large-scale study of Andean forest dynamics using a set of 63 permanent forest plots assembled over the past two decades. In the North-Central Andes tree turnover (mortality and recruitment) and tree growth declined with increasing elevation and decreasing temperature. In addition, basal area increased in Lower Montane Moist Forests but did not change in Higher Montane Humid Forests. However, at higher elevations the lack of net basal area change and excess of mortality over recruitment suggests negative environmental impacts. In North-Western Argentina, forest dynamics appear to be influenced by land use history in addition to environmental variation. Taken together, our results indicate that combinations of abiotic and biotic factors that vary across elevation gradients are important determinants of tree turnover and productivity in the Andes. More extensive and longer-term monitoring and analyses of forest dynamics in permanent plots will be necessary to understand how demographic processes and woody biomass are responding to changing environmental conditions along elevation gradients through this century.

Introduction

Our understanding of the responses of tropical forests to environmental factors is still limited [13]. Indeed, general patterns of forest dynamics and productivity in the Andes Mountains in particular remain poorly characterized due to the scarcity of studies, as well as the complexity of environmental variation in these topographically complex systems [46]. This lack of knowledge is troubling since the Andes are among the most important areas for biological conservation in the world [7] and deliver valuable environmental services (e.g., water provision and carbon storage) to large human populations [8, 9]. In addition, the ecological functioning of the Andes is tightly linked to that of the Amazonian rainforests lying to their eastern flank, which is another area critical for biological conservation. The lack of knowledge about the influence of climate on dynamics and productivity of tropical forests is of special concern in the face of on-going and future climate change [10].

Tree turnover and primary productivity are controlled to a large extent by environmental factors. In mountain ecosystems, tree turnover (the mortality and recruitment of stems) and net primary productivity (gross production of fixed carbon) both tend to decrease with elevation and lower environmental temperatures [1115]. The influence of precipitation on tree turnover and productivity is less clear than that of temperature. In general, forest productivity increases with precipitation but the relationship is not linear [16], and high precipitation may have negative effects on plant production in some tropical mountain forests [17]. There is increasing evidence that temperature and precipitation patterns are shifting in the Andean region [1820], and that Andean forests are already responding to environmental changes. For example, trees in the Peruvian Andes have shifted their distributions upslope over the past decades potentially in response to rising temperatures [21].

The main objective of this paper is to describe how rates of tree turnover, tree growth and basal area net change, a simple proxy for woody productivity, are affected by environmental and geographic variation in the Andes. Thus far, studies on this subject have only been conducted over relatively small geographic scales [12, 13, 15, 22, 23]. Here we used field data collected from permanent forest monitoring plots distributed widely from Colombia to northern Argentina. Using data from these plots, we investigated how tree turnover and growth rates vary across gradients of latitude, elevation, temperature, and precipitation within two regions: the North-Central Andes and North-Western Argentina. Moreover, we evaluated changes of basal area as a function of environmental variation in 32 permanent plots in the Northern-Central Andes. This is the first study to characterize and synthesize patterns of forest dynamics across a large spatial scale in Andean forests.

Materials and Methods

This study was conducted using data collected from 63 permanent forest monitoring plots located across a latitudinal gradient spanning approximately 4000 km (Fig 1) in Colombia, Ecuador, Peru and Argentina. The plots were distributed across a large range of elevations, from 57 to 3940 m above sea level (Fig 1 and S1 Table). Fifty-two study plots were located in the Andean cordillera above 500 m (1812 m ± 121), while 11 were in the Amazonian and Caribbean lowlands adjacent to the cordillera.

thumbnail
Fig 1. Distribution of the 63 permanent forest plots used in this study.

North-central Andes plots were only used for analyses of tree turnover.

https://doi.org/10.1371/journal.pone.0126594.g001

The permanent plots had a uniform distribution across the elevation gradient. The mean elevation of the North-Central Andean plots was 1801 m (range = 3883 m, median 1980 m). Only the permanent plots in southern Peru were placed along a large elevation gradient, from 1800 to 3450 m above sea level (Fig 1). The other sites with a high concentration of permanent plots did not have strong altitudinal variation. On average, Argentinian plots were located at 886 m above sea level (range = 1285, median = 898).

The plots varied in size from 0.4 to 1 ha (S1 Table) and were established by independent teams of researchers for different purposes; however, all plots share a common core set of installation and census methods stated in or comparable to the RAINFOR protocol [24]. Plots were established in mature forests with no, or only minor, signs of recent human influence, and also avoiding natural landscape-scale disturbance processes such as landslides. In each of the study plots, all trees ≥ 10 cm of diameter at breast height (DBH) were identified and measured for diameter at least twice. First censuses were conducted between 1996 and 2009 (mean year of first census = 2003.1±0.34), and plot re-censuses took place between 2000 and 2012 (mean year of last census = 2008.3 ±0.31). On average, plots were thus resampled 5.2 years after their establishment (minimum and maximum resampling periods = 1.75 and 10.16 years, respectively). We used DBH measurements to estimate the total basal area of individuals (cross sectional area of stems at breast height; m2 ha-1) in each census.

There was a negative relationship between altitude and temperature in our data set corresponding to lapse rates of 4.9°C and 5.1°C for every 1 km increase in elevation in the North-Central Andes (y = 27.44–0.0049 * elevation (m), P<0.001) and in the Argentinean sites (y = 23.8–0.0051 * elevation (m), P<0.001). Only in the North-Central Andes precipitation rates were somewhat negatively related with elevation (y = 2971.3–0.2 * elevation (m), P>0.05), as in the Argentinian sites there was no relationship between these two variables (y = 1109–0.01 * elevation (m), P>0.05).

We calculated demographic rates for each plot using the following equations: Annual mortality rates were estimated as m = 1-[1- (N0—N1)/ N0]1/t * 100, annual recruitment rates (trees reaching 10 cm DBH): r = 1-(1- Nr—N1)1/t * 100; where N0 is the number of individuals alive in the initial census, N1 is the number of individuals surviving to the second census, Nr is the number of individuals recruited between censuses, t corresponds to time in years. Turnover was estimated as the mean of tree mortality and recruitment [25]. Basal area net change was estimated as BAn = ((BA1-BA0)/ BA0)/t * 100 and relative tree growth rate was calculated as g = ((BA1-BAr)-(BA0-BAd) / (BA0- BAr))/t * 100 BA0 and BA1 correspond to basal area values at the initial and last census, respectively; BAr is the basal area of trees recruited, and BAd, is the basal area of dead trees in the last census, t is time in years [26]. Tree growth was estimated as changes of the basal area of trees surviving through the census period. To correct for census interval effects in estimates of turnover for plots in the North-Central Andes we used the formula ʎcorr = ʎ t0.08, where ʎcorr is turnover rate standardized to a 1-year census interval, ʎ is raw turnover rate, ant t is the length of census interval in years [27]. Demographic rates for each of the individual study plots are presented in the S2 and S3 Tables.

We analysed the relationships between our measures of forest demography and several environmental and geographic variables (Table 1). Our analyses included the following information for each plot: latitude, elevation, total annual precipitation, minimum temperature (i.e., average monthly minimum temperatures), and a pluvio-thermic index (Iod2). Iod2 = Pp/Tp where Pp is the total precipitation in mm, and Tp is the sum of the monthly temperature higher than 0°C in tenths of °C for the two consecutive driest months of the year [28]. Pluvio-thermic (or ombrothermic) indices express the relationship between precipitation and temperature at a given site and thus are used to classify vegetation in categories that range from dry to hyper-humid (dry sites have low scores and humid sites have high scores), providing a way to express the outcome in one index of the large variations in rates of both moisture supply and moisture loss that exist in the geographically highly variable climates of our study region. Latitudinal location, minimum annual temperatures, and the pluvio-thermic index were included in the statistical models described below to account for extremes of environmental variation, which were partly related to the geographic distribution of our forest plots. Minimum annual temperature, total annual precipitation, and the variables needed to calculate the pluvio-thermic index were obtained from the WorldClim extrapolated climate database (www.worldclim.org) [29]. For some plots, total annual precipitation rates were recorded at nearby meteorological stations. Data analyses were conducted using the statistical package JMP [30]. The ethics statement is presented in the S1 File.

thumbnail
Table 1. Environmental and demographic forest variables in the North-Central Andes, and North-Western Argentina.

https://doi.org/10.1371/journal.pone.0126594.t001

Data Analyses

We used subsets of our data base to explore aspects of forest demography in areas with different environmental or geographical features. Hence, responses of tree turnover, growth, and basal area net change to environmental variation were examined separately for the 45 North-Central Andean plots and for the 18 North-Western Argentina plots. In addition, patterns of forest demography were examined for 32 North-Central Andean plots, considering two forest groups, as explained below.

We conducted Principal Components Analysis (PCA) to characterize the environmental heterogeneity recorded in our 63 forest monitoring plots in fewer and uncorrelated variables. The PCA included data on elevation, latitude, minimum temperature, total annual precipitation, and the pluvio-thermic index (Table 2). We conducted linear regression analyses to examine the relationships between the first two PCA axis and tree turnover, tree growth rates, and basal area net change. Analyses were conducted separately for two groups of plots, those located in the North-Central Andes (Colombia to Peru), and plots in North-Western Argentina.

thumbnail
Table 2. Loadings and variation explained in two Principal Components Analysis using five environmental variables recorded in 63 permanent plots, and for a subset of 32 permanent plots located in the North-Central Andes.

https://doi.org/10.1371/journal.pone.0126594.t002

A second PCA included the same environmental variables considered above, but only for 32 northern plots (Fig 1). The scores of each plot in the first two PCA axes were used to identify two groups of plots with contrasting environmental conditions (Figs 1 and 2, Tables 2 and 3). The group of Higher Montane Humid Forests (HMHF) included permanent plots located at higher elevation, with lower minimum temperature, receiving high precipitation rates, and higher pluvio-thermic scores compared to Lower Montane Moist Forests (LMMF). Latitude also differed between forest groups; HMHFs were located toward the south (lower latitude) and LMMFs occurred at northern locations (higher latitude). Each forest group had 16 permanent plots. We used ANOVA to compare environmental, geographic features, and forest demography between both groups (Table 3).

thumbnail
Fig 2. Classification of 32 North-Central Andean permanent plots.

Classification according to the plot’s scores in the first two factors of a Principal Component Analysis.

https://doi.org/10.1371/journal.pone.0126594.g002

thumbnail
Table 3. Descriptive statistics of the geographic, environmental, and demographic variables of Higher Montane Humid Forests (HMHF) and Lower Montane Moist Forests (LMMF).

https://doi.org/10.1371/journal.pone.0126594.t003

Results

The PCA, including 63 permanent plots, produced a first PCA axis (PCA 1) that captured 45% of the environmental variation and was positively correlated with latitude, total annual precipitation, and high pluvio-thermic scores. The second PCA axis (PCA 2) described 33% of the environmental variation and was positively correlated with elevation and negatively correlated with minimum annual temperature. Combined, both axes explained 78% of the environmental variation (Table 2).

North-Central Andes

Our results showed that forest demography in this region was affected by the environmental variation captured in PCA 2. Rates of tree turnover, total and relative growth, and basal area net change were negatively related to PCA 2, indicating that these demographic variables decreased with increasing elevation and declining minimum temperatures (Fig 3, Table 4, S4 Table). None of our demographic variables were related to PCA 1.

thumbnail
Fig 3. Linear regressions of forest demographic variables as a function of the first two factors of a Principal Component Analysis.

https://doi.org/10.1371/journal.pone.0126594.g003

thumbnail
Table 4. Summary of the results of Linear Regression analyses using PCA axes as predictive variables of forest demography for plots in North-Central, and North-Western Argentina.

https://doi.org/10.1371/journal.pone.0126594.t004

North-Western Argentina

We found a negative relationship between PCA 1 and basal area net change, as it decreased with increasing latitude, annual precipitation and pluvio-thermic scores (Fig 4, Table 4, S4 Table). We also found that in this region tree growth was positively related to PCA 2, thus it increased with elevation and with minimum environmental temperatures. Neither tree turnover nor relative growth was related to our PCA axes (S1 Fig).

thumbnail
Fig 4. Linear regressions of forest demographic variables in North-Western Argentina.

Basal area net change and individual tree growth and as a function of the first two factors of a Principal Component Analysis.

https://doi.org/10.1371/journal.pone.0126594.g004

North-Central Andes: Higher Montane Humid and Lower Montane Moist Forests

The first PCA axis captured 48% of the environmental variation, and it was positively correlated to low elevation, high minimum temperatures, and high latitude (Table 2, Fig 2). The second PCA axis explained 32% of the variation and was positively correlated with high annual precipitation rates and high pluvio-thermic scores. Therefore, the first two axes captured 80% of the environmental variation.

Our analyses showed that LMMFs increased in basal area over the study period (mean = 0.84 ± 0.26%), whereas HMHFs did not have statistically significant changes in net basal area change (mean = -0.19 ± 0.25%), and that these differences were significantly different (Table 3). In addition, LMMFs also had significantly higher relative tree growth rates than HMHFs, although this pattern did not hold for total tree growth rates (Table 3). In addition, tree recruitment was faster in LMMFs than at HMHFs. Mortality rates did not differ between forest groups (Table 3). Within forest groups, tree mortality was significantly higher than tree recruitment in HMHFs (F1,15 = 3.76, P<0.001), but these demographic rates did not differ in LMMFs (F1,15 = 1.13, P = 0.275).

Discussion

This study shows that the rates of tree turnover and basal area net change of Andean forests are strongly related to large-scale variation in geographic and environmental factors. Moreover, we found contrasting patterns of basal area change in forests with different environmental features. Tree turnover and tree growth decreased with elevation and with lower minimum temperatures in the North-Central Andean forests (Fig 3). The same tendency was observed for basal area net change. Thus, woody growth declined with elevation at this latitudinal zone. This pattern has been found several times in temperate regions [11, 31], and at smaller spatial scales in the Andes [12, 13, 15, 23, 32]; however, this is the first time this elevational trend is documented across a large spatial scale in the tropical Andes.

Given the strong connection between elevation and temperature, the decrease in tree growth at high elevations is possibly due to effects of low temperature on plant metabolism [33] and its indirect effect on nutrient cycling [23, 32]. Other abiotic factors that reduce plant growth often vary with elevation and may contribute to explain this trend. For example, deceasing soil fertility with elevation reduces productivity at higher altitudes in certain Andean forests [23, 32], and at other mountain areas [31]. In the same way, increasing cloudiness at higher elevations decreases direct solar radiation and reduces productivity [14, 34].

Biotic factors can also increase tree turnover at lower elevations. For example high rates of tree mortality at low elevations may be due to increased aboveground competition among tree species, stronger top-down regulation exerted by herbivores, disease [35], and parasitism [31]. In addition, tree communities at lower elevations generally have higher relative abundances of tree species that favour growth over defence strategies, which results in high recruitment rates [35].

It is interesting to note that tree turnover, tree growth rates, and basal area net change were not associated with PCA 1, which described a gradient of low to high precipitation, pluvio-thermicity, and latitude (Fig 2). The lack of demographic responses to these environmental factors is likely related to methodological issues, rather than a lack of their effect on tree turnover and basal area net change, as it has been demonstrated in previous studies [17, 36, 37]. Our dataset is certainly an incomplete sample of the precipitation gradient in the Andes. In addition, we used modeled rather than field climatic variables, which may be problematic due to the complexity of modeling climatic variables in mountain regions [38], where ‘horizontal precipitation’ and cloud cover may affect the amount of water available for plant growth [34]. Using long-term data taken within narrower latitudinal zones in the Andes may be a first step to understand how precipitation affects forest demography in this region.

The patterns of forest demography in North-Western Argentina differed from those of the North-Central Andes. Here basal area net change decreased with higher precipitation, pluvio-thermicity, and latitude (PCA 1) and tree growth increased with elevation and with lower minimum temperatures (PCA 2) (Table 2, Fig 3). These unexpected patterns are explained by the abiotic factors considered and ecological processes related to forest recovery. The results were driven to a large extent by the more humid Tucuman plots located at higher elevations (1000–1700 m), and at lower latitude (Fig 4). Tucuman plots had higher pluvio-thermic index scores (but not higher precipitation rates) compared to the other two Argentinian sites, Salta and Jujuy. These plots had a 50% increase in stem density and 6% in basal area between 1992 and 2007, which has been attributed to forest recovery after the removal of livestock (i.e., cattle and horses) [39] once the area was declared a natural reserve in 1973 [40]. Hence, higher tree growth rates at higher elevation sites may also be due to the recovery of populations of trees palatable to herbivores, which may favor growth over persistence trade-offs.

Taken together, our results indicate that combinations of abiotic and biotic factors that vary across elevation gradients are important determinants of tree turnover and productivity in the Andes. It is noteworthy that this pattern emerges despite the large variation in species composition (among many other variables) associated to the geographical variation considered in the study. It is likely that the effect of our environmental variables on forest demography vary in strength and even in sign at narrower latitudinal zones, that are represented as different countries or sites in this study (Figs 3 and 4). In addition, in the North-Western Argentina, where we considered a narrower latitudinal zone, these environmental controls of tree turnover and productivity were modified by past land-use and ecological recovery.

North-Central Andes: Higher Montane Humid and Lower Montane Moist Forests

Our results overall indicate a balance between stem recruitment and mortality at lower elevation moist forests (LMMFs), while all have increased basal area. In contrast, at higher elevations the plots showed an excess of mortality over tree recruitment, and no significant change in basal area (Table 3). The increase in basal area at lower elevations could be caused by a number of reasons. It is for example consistent with evidence that Amazonian and other lowland forests have experienced net biomass increases in recent decades [1, 4143], which has been attributed to increased nutrient deposition, atmospheric CO2, and shifts on precipitation regimes [44], but other factors may also be in play in the Andes. The lack of basal area net change at higher elevations and excess of mortality over recruitment suggests negative environmental impacts.

Globally, as well as in tropical forests, drought and higher environmental temperatures affect tree mortality [45, 46], and reduce aboveground productivity [47]. Between 2000–2012, a period of sampling of most of our permanent plots, some areas in the Northern and Central Andes experienced lower precipitation rates than those recorded in previous decades [47]. Thus, higher environmental temperatures could be interacting with lower precipitation rates to affect tree performance and mortality in HMHFs. Lower precipitation rates or short drought periods can enhance tree mortality even in forests that are not considered to be water-limited [4850]. Hydraulic stress can kill trees through xylem embolism or carbon starvation and by enhancing negative effects of insect attacks and disease [46].

High tree mortality rates reflect various biological responses to environmental change in our study area. For instance, upward migration of Andean trees, possibly linked to higher environmental temperatures [21] may not be a process in which mortality and recruitment rates occur at the same time. Periods of high mortality and low recruitment of species sensitive to higher environmental temperatures may take place. During such periods, low recruitment rates may be magnified by lagged migration of species adapted to higher environmental temperatures. More detailed, longer-term analyses of forest dynamics in permanent plots are clearly necessary to understand how demographic processes and woody biomass are responding to changing environmental conditions along elevation gradients.

Taken together our results suggest a high sensitivity of Andean forest dynamics to climate. Perhaps more importantly, these results suggest that Andean forests are responding to the changes in climate that are already occurring in this region and that are predicted to occur in the future. Further, Andean montane forests seem to be responding to their changing environment in different ways than lowland tropical forests. In the Andes, mean annual temperatures increased 0.34°C decade-1 over the period from 1974–1998 [51], a rate much higher than that observed for the Neotropical region as a whole [52]. Moreover, temperature increases have been more pronounced at higher elevations in the Andes [53], which may cause stronger biological responses in these areas. These findings provide a baseline for future studies exploring the responses of Andean tree species to environmental variation and to future climate scenarios. This on-the-ground monitoring network is still at an early stage. Detailed analyses of the existing information in the context of forest responses to directional climate change, and further field monitoring will help to understand if the small changes observed so far are part of a larger, cyclical pattern, perhaps induced by decadal climate fluctuations.

Supporting Information

S1 Table. Location, geographical and environmental features, and establishment information of the permanent plots used in this study.

Country codes: Argentina (ARG), Colombia (COL), Ecuador (ECU), Peru (PER). Authors codes: AD: A. Duque, AM: A. Malizia, CB: C. Blundo, EA: E. Álvarez, JC: J. Carilla, JH: J. Homeier, KJF: K.J. Feeley, LM: L. Malizia, NA: N. Aguirre, OO: O. Osinaga, SB: S. Báez, ZA: Z. Aquirre, WF: W. Farfán, RLP: R. Linares-Palomino.

https://doi.org/10.1371/journal.pone.0126594.s001

(DOCX)

S2 Table. Demographic rates, and dominant or most common species of the 45 permanent plots located in the North-Central Andes.

“Forest group” indicates the classification of the forest plot in one of the groups: HMHFs = Higher Montane Humid Forest, LMMFs = Lower Montane Moist Forests. Country codes: COL = Colomobia, ECU = Ecuador, PER = Peru. Upperscripts on the Plot code indicate the methodological protocol used. The scientific nomenclature was updated according to the Global Biodiversity Information Facility databases (GBIF; www.gbif.org).

https://doi.org/10.1371/journal.pone.0126594.s002

(DOCX)

S3 Table. Demographic rates, and dominant or most common species of the 18 permanent plots located in South-Western Argentina.

Scientific nomenclature after the Global Biodiversity Information Facility databases (GBIF; www.gbif.org).

https://doi.org/10.1371/journal.pone.0126594.s003

(DOCX)

S4 Table. Results of the Linear regression analyses using PCA factors as predictors of forest demography in North-Central Andean, and North-Western Argentina.

https://doi.org/10.1371/journal.pone.0126594.s004

(DOCX)

S1 Fig. Graphs of the Linear regression analyses using PCA factors as predictors of forest demography in North-Western Argentina.

https://doi.org/10.1371/journal.pone.0126594.s006

(TIFF)

Acknowledgments

The comments and suggestions of Paul V.A. Fine, Kristina Anderson-Teixeira, and an anonymous reviewer greatly improved this study. EAD thanks the Cooperation Agreement Colciencias-JBMED (No. 393 of 2012) ECOPARQUE Los Besotes Colombia for permission to work in the protected area and to Oswaldo Velasquez for forest dynamics data of Farallones. MRS, KJF, KGC and WF thank the Andes Biodiversity and Ecosystems Research Group (http://www.andesconservation.org/), SERNANP Peru and the personnel of Manu National Park for the assistance with logistics and permission to work in the protected area, and the Amazon Conservation Association (ACA) for the logistic support. RLP, CR and MA thank APRODES (Asociación Peruana para la Promoción del Desarrollo Sostenible). SB thanks Wouter Buytaert for support to obtain forest dynamics data in Southern Ecuador.

Author Contributions

Conceived and designed the experiments: SB AM JC AD JH KJF. Performed the experiments: SB AM JC CB MA NA ZA EA FC AD WFR KGC RG JH RLP LM OMC OO OP CR MRS KJF. Analyzed the data: SB JH AM JC CB. Wrote the paper: SB KJF JH OP AM JC RG.

References

  1. 1. Lewis SL, Phillips OL, Baker TR, Lloyd J, Malhi Y, Almeida S, et al. Concerted changes in tropical forest structure and dynamics: evidence from 50 South American long-term plots. Philosophical Transactions of the Royal Society of London Series B: Biological Sciences. 2004;359(1443):421–36. pmid:15212094
  2. 2. Clark DA. Sources or sinks? The responses of tropical forests to current and future climate and atmospheric composition. Philosophical Transactions of the Royal Society of London Series B: Biological Sciences. 2004;359(1443):477–91. pmid:15212097
  3. 3. Wright SJ. Tropical forests in a changing environment. Trends in Ecology & Evolution. 2005;20:553–60.
  4. 4. Pitman N, Widmer J, Jerkins CN, Stocks G, Seales L, Paniagua F, et al. Volume and Geographical Distribution of Ecological Research in the Andes and the Amazon, 1995–2008. Tropical Conservation Science. 2011;4:64–81.
  5. 5. Killeen TJ, Douglas M, Consiglio T, Jørgensen PM, Mejia J. Dry spots and wet spots in the Andean hotspot. Journal of Biogeography. 2007;34(8):1357–73.
  6. 6. Buytaert W, Vuille M, Dewulf A, Urrutia R, Karmalkar A, Célleri R. Uncertainties in climate change projections and regional downscaling in the tropical Andes: implications for water resources management. Hydrol Earth Syst Sci. 2010;14(7):1247–58.
  7. 7. Myers N, Mittermeier RA, Mittermeier CG, Fonseca GAd, Kent J. Biodiversity hotspots for conservation priorities. Nature. 2000;403:843–58. pmid:10706267
  8. 8. Cincotta RP, Wisnewski J, Engelman R. Human population in the biodiversity hotspots. Nature. 2000;404(6781):990–2. pmid:10801126
  9. 9. Spracklen DV, Righelato R. Tropical montane forests are a larger than expected global carbon store. Biogeosciences. 2014;11(10):2741–54.
  10. 10. Silva CE, Kellner JR, Clark DB, Clark DA. Response of an old-growth tropical rainforest to transient high temperature and drought. Global Change Biology. 2013;19(11):3423–34. pmid:23824759
  11. 11. Stephenson NL, van Mantgem PJ. Forest turnover rates follow global and regional patterns of productivity. Ecology Letters. 2005;8(5):524–31. pmid:21352456
  12. 12. Girardin CAJ, Malhi YM, Aragao LEOC, Mamani M, Huaraca Huasco W, Durand L, et al. Net primary productivity allocation and cycling of carbon along a tropical forest elevational transect in the Peruvian Andes. Global Change Biology. 2010.
  13. 13. Homeier J, Breckle S-W, Günter S, Rollenbeck RT, Leuschner C. Tree Diversity, Forest Structure and Productivity along Altitudinal and Topographical Gradients in a Species-Rich Ecuadorian Montane Rain Forest. Biotropica. 2010;42(2):140–8.
  14. 14. Huasco WH, Girardin CAJ, Doughty CE, Metcalfe DB, Baca LD, Silva-Espejo JE, et al. Seasonal production, allocation and cycling of carbon in two mid-elevation tropical montane forest plots in the Peruvian Andes. Plant Ecology & Diversity. 2013;7(1–2):125–42.
  15. 15. Leuschner C, Zach A, Moser G, Homeier J, Graefe S, Hertel D, et al. The carbon balance of tropical mountain forests along an altitudinal transect, Southern Ecuador. In: Bendix J, Beck E, Bräuning A, Makeschin F, Mosandl R, Scheu S, et al., editors. Ecosystem services, biodiversity and environmental change in a tropical mountain ecosystem of South Ecuador Ecological Studies 221. Berlin: Springer; 2013. p. 117–39.
  16. 16. Anderson-Teixeira KJ, Miller AD, Mohan JE, Hudiburg TW, Duval BD, DeLucia EH. Altered dynamics of forest recovery under a changing climate. Global Change Biology. 2013;19(7):2001–21. pmid:23529980
  17. 17. Schuur EAG, Matson PA. Net primary productivity and nutrient cycling across a mesic to wet precipitation gradient in Hawaiian montane forest. Oecologia. 2001;128(3):431–42. pmid:24549913
  18. 18. Ruiz D, Moreno HA, Gutiérrez ME, Zapata PA. Changing climate and endangered high mountain ecosystems in Colombia. Science of The Total Environment. 2008;398(1–3):122–32. http://dx.doi.org/10.1016/j.scitotenv.2008.02.038. pmid:18452974
  19. 19. Mora DE, Campozano L, Cisneros F, Wyseure G, Willems P. Climate changes of hydrometeorological and hydrological extremes in the Paute basin, Ecuadorean Andes. Hydrol Earth Syst Sci. 2014;18(2):631–48.
  20. 20. Urrutia R, Vuille M. Climate change projections for the tropical Andes using a regional climate model: Temperature and precipitation simulations for the end of the 21st century. Journal of Geophysical Research: Atmospheres. 2009;114(D2):D02108.
  21. 21. Feeley KJ, Silman MR, Bush M, Farfan W, Garcia Cabrera K, Malhi Y, et al. Upslope migration of Andean trees. Journal of Biogeography. 2011;38:783–91.
  22. 22. Delaney M, Brown S, Lugo AE, Torres Lezama A, Bello Quintero N. The distribution of organic carbon in major components of forests located in five life zones of Venezuela. Journal of Tropical Ecology. 1997;13:697–708.
  23. 23. Unger M, Homeier J, Leuschner C. Effects of soil chemistry on tropical forest biomass and productivity at different elevations in the equatorial Andes. Oecologia. 2012;170:263–74. pmid:22410639
  24. 24. Phillips O, Baker T. Manual de Campo para la Remedición y Establecimiento de Parcelas. RAINFOR. Project for the Advance of Networked Science in Amazonia. Sixth frame-work Programme (2002–2006): RAINFOR. Project for the Advance of Networked Science in Amazonia. Sixth Framework Programme (2002–2006); 2002.
  25. 25. Bellingham PJ, Sparrow AD. Multi-stemmed trees in montane rain forests: their frequency and demography in relation to elevation, soil nutrients and disturbance. Journal of Ecology. 2009;97(3):472–83.
  26. 26. Peña MA, Duque A. Patterns of stocks of aboveground tree biomass, dynamics, and their determinants in secondary Andean forests. Forest Ecology and Management. 2013;302(0):54–61. http://dx.doi.org/10.1016/j.foreco.2013.03.025.
  27. 27. Lewis SL, Phillips OL, Sheil D, Vinceti B, Baker TR, Brown S, et al. Tropical forest tree mortality, recruitment and turnover rates: calculation, interpretation and comparison when census intervals vary. Journal of Ecology. 2004;92(6):929–44.
  28. 28. Rivas-Martínez S, Rivas-Saenz S. Worldwide Bioclimatic Classification System, 1996–2009 Spain: Phytosociological Research Center; 2008. Available from: http://www.globalbioclimatics.org/book/bioc/global_bioclimatics-2008_00.htm.
  29. 29. Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A. Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology. 2005;25:1965–78.
  30. 30. Inc SI. JMP. 9 ed. Cary, NC,1989–2007.
  31. 31. Sundqvist MK, Sanders NJ, Wardle DA. Community and ecosystem responses to elevational gradients: processes, mechanisms, and insights for global change. Annual Review of Ecology, Evolution, and Systematics. 2013;44(1):261–80.
  32. 32. Wolf K, Veldkamp E, Homeier J, Martinson GO. Nitrogen availability links forest productivity, soil nitrous oxide and nitric oxide fluxes of a tropical montane forest in southern Ecuador. Global Biogeochemical Cycles. 2011;25(4):GB4009.
  33. 33. Brown JH, Gillooly JF, Allen AP, Savage VM, West GB. Toward a metabolic theory of ecology. Ecology. 2004;85(7):1771–89.
  34. 34. Bruijnzeel LA, Veneklaas EJ. Climatic conditions and tropical monante forest productivity: the fog has not lifted yet. Ecology. 1998;79(1):3–9.
  35. 35. Stephenson NL, van Mantgem PJ, Bunn AG, Bruner H, Harmon ME, O'Connell KB, et al. Causes and implications of the correlation between forest productivity and tree mortality rates. Ecological Monographs. 2011;81(4):527–55.
  36. 36. Cleveland CC, Townsend AR, Taylor P, Alvarez-Clare S, Bustamante MMC, Chuyong G, et al. Relationships among net primary productivity, nutrients and climate in tropical rain forest: a pan-tropical analysis. Ecology Letters. 2011;14(9):939–47. pmid:21749602
  37. 37. Raich JW, Russell AE, Kitayama K, Parton WJ, Vitousek PM. Temperature influences carbon accumulation in moist tropical forests. Ecology. 2006;87(1):76–87. pmid:16634298
  38. 38. Daly C, Halbleib M, Smith JI, Gibson WP, Doggett MK, Taylor GH, et al. Physiographically sensitive mapping of climatological temperature and precipitation across the conterminous United States. International Journal of Climatology. 2008;28(15):2031–64.
  39. 39. Malizia A, Easdale T, Grau HR. Rapid structural and compositional change in an old-growth subtropical forest: Using plant traits to identify probable drivers. PlosOne. 2013. https://doi.org/10.1371/journal.pone.0073546
  40. 40. Grau HR, Hernández ME, Gutierrez J, Gasparri NI, Casavecchia C, Flores-Ivaldi EE, et al. A Peri-urban neotropical forest transition and its consequences for environmental services. Ecology and Society. 2008;13. http://www.ecologyandsociety.org/vol13/iss1/art35/.
  41. 41. Phillips OL, Baker TR, Arroyo L, Higuchi N, Killeen TJ, Laurance WF, et al. Pattern and process in Amazon tree turnover, 1976–2001. Philosophical Transactions of the Royal Society of London Series B: Biological Sciences. 2004;359(1443):381–407. pmid:15212092
  42. 42. Chave J, Condit R, Muller-Landau HC, Thomas SC, Ashton PS, Bunyavejchewin S, et al. Assessing Evidence for a Pervasive Alteration in Tropical Tree Communities. PLoS Biol. 2008;6(3):e45. pmid:18318600
  43. 43. Phillips OL, Malhi Y, Higuchi N, Laurance WF, Núñez PV, Vásquez RM, et al. Changes in the Carbon Balance of Tropical Forests: Evidence from Long-Term Plots. Science. 1998;282(5388):439–42. pmid:9774263
  44. 44. Lewis SL, Malhi Y, Phillips OL. Fingerprinting the impacts of global change on tropical forests. Philosophical Transactions of the Royal Society of London Series B: Biological Sciences. 2004;359(1443):437–62. pmid:15212095
  45. 45. Clark DB, Clark DA, Oberbauer SF. Annual wood production in a tropical rain forest in NE Costa Rica linked to climatic variation but not to increasing CO2. Global Change Biology. 2009;16(2):747–59.
  46. 46. Allen CD, Macalady AK, Chenchouni H, Bachelet D, McDowell N, Vennetier M, et al. A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. Forest Ecology and Management. 2010;259(4):660–84. http://dx.doi.org/10.1016/j.foreco.2009.09.001.
  47. 47. Hilker T, Lyapustin AI, Tucker CJ, Hall FG, Myneni RB, Wang Y, et al. Vegetation dynamics and rainfall sensitivity of the Amazon. Proceedings of the National Academy of Sciences. 2014. https://doi.org/10.1073/pnas.1404870111
  48. 48. Condit R, Hubbell SP, Foster RB. Mortality Rates of 205 Neotropical Tree and Shrub Species and the Impact of a Severe Drought. Ecological Monographs. 1995;65(4):419–39.
  49. 49. Chazdon RL, Redondo Brenes A, Vilchez Alvarado B. Effects of climate and stand age on annual tree dynamics in tropical secondary growth rain forests. Ecology. 2005;86(7):1808–15.
  50. 50. Phillips OL, van der Heijden G, Lewis SL, López-González G, Aragão LEOC, Lloyd J, et al. Drought–mortality relationships for tropical forests. New Phytologist. 2010;187(3):631–46. pmid:20659252
  51. 51. Vuille M, Bradley RS, Werner M, Keimig F. 20th century climate change in the tropical Andes: observations and model results. Climate Change. 2003;59:75–99.
  52. 52. Malhi Y. The carbon balance of tropical forest regions, 1990–2005. Current Opinion in Environmental Sustainability. 2010;2(4):237–44. http://dx.doi.org/10.1016/j.cosust.2010.08.002.
  53. 53. Vuille M, Bradley R. Mean temperature trends and their vertical structure in the tropical Andes. Geophysical Research Letters. 2000;27:3885–8.