Articles | Volume 18, issue 4
Biogeosciences, 18, 1525–1541, 2021
Biogeosciences, 18, 1525–1541, 2021

Research article 03 Mar 2021

Research article | 03 Mar 2021

Factors controlling the productivity of tropical Andean forests: climate and soil are more important than tree diversity

Factors controlling the productivity of tropical Andean forests: climate and soil are more important than tree diversity
Jürgen Homeier1,2 and Christoph Leuschner1,2 Jürgen Homeier and Christoph Leuschner
  • 1Plant Ecology and Ecosystems Research, University of Göttingen, Untere Karspüle 2, 37073 Göttingen, Germany
  • 2Centre for Biodiversity and Sustainable Land Use, University of Göttingen, Untere Karspüle 2, 37073 Göttingen, Germany

Correspondence: Jürgen Homeier (


Theory predicts positive effects of species richness on the productivity of plant communities through complementary resource use and facilitative interactions between species. Results from manipulative experiments with tropical tree species indicate a positive diversity–productivity relationship (DPR), but the existing evidence from natural forests is scarce and contradictory. We studied forest aboveground productivity in more than 80 humid tropical montane old-growth forests in two highly diverse Andean regions with large geological and topographic heterogeneity and related productivity to tree diversity and climatic, edaphic and stand structural factors with a likely influence on productivity. Main determinants of wood production in the perhumid study regions were elevation (as a proxy for temperature), soil nutrient (N, P and base cation) availability and forest structural parameters (wood specific gravity, aboveground biomass). Tree diversity had only a small positive influence on productivity, even though tree species numbers varied largely (6–27 species per 0.04 ha). We conclude that the productivity of highly diverse Neotropical montane forests is primarily controlled by thermal and edaphic factors and stand structural properties, while tree diversity is of minor importance.

1 Introduction

Research into the diversity–productivity relationship (DPR) has recently shifted its focus on forests because of their importance for the global carbon cycle and the high biodiversity, especially of tropical forests (Mori et al., 2017; Fei et al., 2018). Results from mixed species plantations with tropical trees indicate a positive DPR in the majority of cases (Sapijanskas et al., 2014; Huang et al., 2018; Schnabel et al., 2019), suggesting complementary resource use and facilitative interactions between different tree species (McIntire and Fajardo, 2014; Chen et al., 2016) The results of these experiments may provide valuable insights into the mechanisms of diversity effects on productivity (Tuck et al., 2016; Schnabel et al., 2019), but they are difficult to extrapolate to tropical old-growth forests due to the young age and cohort structure of these tree plantations. Observational studies on the linkage between tree diversity and productivity in natural or near-natural tropical forests are usually burdened with additional confounding factors and thus are more difficult to interpret than experiments, but they have the advantage of considering natural systems. The few existing observational studies on the DPR in tropical forests have produced contradictory evidence. On Mt. Kinabalu, Malaysia, aboveground net primary production (NPP) was higher in more diverse tropical montane forest plots (Aiba et al., 2005), but no effects of species richness on productivity were found in logged tropical lowland forests in Guyana (van der Sande et al., 2018). Chisholm et al. (2013) analyzed the tree diversity–productivity relation in 25 permanent forest plots from temperate to tropical regions and found a positive diversity effect when plot size was small (0.04 ha), but the effect disappeared at larger spatial scales (0.25 and 1 ha). In a global analysis of forest productivity data, Liang et al. (2016) found that positive DPRs predominate, but negative relations do also exist, though they are restricted to a few of the tested data sets.

The existing observational data from tropical forests suggest that (i) tree diversity may enhance productivity under certain conditions, and (ii) a positive DPR might be more prominent in small plots. To date, it is not clear under which abiotic conditions and stand structural settings a positive DPR does emerge and how important tree diversity is compared to abiotic and other biotic determinants of productivity.

The search for abiotic determinants of tropical-forest productivity has mainly concentrated on soil nutrient availability, solar irradiance and precipitation, revealing clear dependencies in studies along environmental gradients (Schuur, 2003; Girardin et al., 2014; Malhi et al., 2017). Many tropical forests, especially those in the lowlands, seem to be P-limited (Aragão et al., 2009; Quesada et al., 2012), while N limitation may be more influential in tropical montane forests (TMFs) (Tanner et al., 1998; Moser et al., 2008; Homeier et al., 2012; Fisher et al., 2013). Fertilization experiments and studies on foliar nutrient levels across environmental gradients suggest that the NPP of tropical forests tends to increase with soil fertility, notably N and P availability, but locally also with base cation supply (Kitayama and Aiba, 2002; Cleveland et al., 2011; Malhi et al., 2004; Homeier et al., 2010; Unger et al., 2012; Fisher et al., 2013; Banin et al., 2014; Hofhansel et al., 2015). Bruijnzeel and Veneklaas (1998) assumed that the productivity of TMFs is primarily limited by cloudiness, i.e., low solar radiation, which is supported by a modeling study (Fyllas et al., 2017). Precipitation should play a minor role in TMFs (Zimmermann et al., 2010). Less clear is the role of temperature. The elevational temperature decrease can act on plants in multiple ways: directly through a limitation of photosynthetic activity and of carbon, water and nutrient cycling in the plant, and it may influence stem, leaf and root morphology. The impairment of nutrient supply by low soil temperatures is an indirect effect. Yet, plants can adapt to lower temperatures, thereby partly reducing these limitations (Lambers and Oliveira, 2019). The majority of studies along elevation transects found an NPP decline with decreasing temperature in tropical mountains (Kitayama and Aiba, 2002; Girardin et al., 2010; Cleveland et al., 2011; Leuschner et al., 2013; van de Weg et al., 2014). However, an analysis of a global database failed to show a temperature dependence of NPP in tropical forests (Luyssaert et al., 2007). This fits to modeling results suggesting that trait variation with elevational species turnover in TMFs may partly capture the effect of temperature on tree metabolism and growth (Fyllas et al., 2017).

Biotic controls related to stand structural properties, notably stem density, basal area and wood density, and the functional composition of the community can also affect forest productivity. Some studies found that forest NPP increased with standing aboveground biomass (AGB) (Keeling and Phillips, 2007; Pan et al., 2013; Lohbeck et al., 2015), but the relation does not seem to be very close since woody tissue residence time, i.e., the ratio between AGB and aboveground productivity, varies considerably with soil fertility and climate in tropical moist forests (Malhi, 2012; Quesada et al., 2012; Chisholm et al., 2013; Malhi et al., 2017). Variation in leaf area index and leaf functional traits such as LMA (leaf mass per area) and foliar N and P content should be major determinants of net primary productivity (van de Weg et al., 2014; Wittich et al., 2012; Finegan et al., 2015) as they influence the photosynthetic capacity of the canopy and radiation interception. Indeed, the modeling study of Fyllas et al. (2017) indicated for tropical forests along an Andean elevation gradient that the influence of LMA, foliar N and P content, and wood density of the species was a more important determinant of NPP than temperature, which is in line with the findings of Luyssaert et al. (2007). From the partly contradicting results we conclude that a multi-factorial approach is needed to understand the determinants of tropical-forest productivity. Yet, most of the existing studies have investigated only one or two of the mentioned factors.

Tropical montane forests cover all humid tropical mountains between  1000 and 3000 m a.s.l., with the widest distribution in the Andes and Central American mountains and more restricted occurrence in the mountains of the Paleotropis. They differ from tropical lowland forest in terms of lower canopy heights, higher stem densities and lower AGB. TMFs play important roles in the provision of water for human populations and as stores of carbon, and tropical mountains are known as hotspots of biodiversity (Ashton, 2003; Bruijnzeel et al., 2010; Spracklen and Righelato, 2014; Fahey et al., 2015; Rahbeck et al., 2019). TMFs occur along steep thermal and edaphic gradients, making them attractive for the study of abiotic controls of tropical-forest productivity and the role of tree diversity. In this study in two Andean elevation transects, we combine biomass, productivity and tree diversity data with comprehensive soil chemical data to identify the most important abiotic and biotic drivers of forest productivity in Neotropical montane forests (Fig. 1). In permanent plots in old-growth forests, we measured coarse-wood production (WP) and fine-litter production (only one of the transects) and related these components of aboveground primary production to the availability of the five macronutrients N, P, Ca, K and Mg in the soil; to elevation as a proxy for temperature; and to stand structural properties and tree diversity in the plots. We chose a plot size of 0.04 ha (20 m × 20 m) in order to meet plot homogeneity requirements in the rugged terrain. Moreover, assumed diversity effects are more likely to emerge at small spatial scales (Chisholm et al., 2013). This allowed us to investigate a larger number of plots across extended environmental gradients (2000 m of elevation distance or  10 K variation in mean annual temperature; large bedrock, topographic and soil variation) and to encounter considerable variation in tree diversity (6 to 27 species per 0.04 ha). The analyses were restricted to old-growth forests with closed canopies in order to exclude variation in stand structure and species composition related to successional processes. We employed structural equation modeling for exploring mutual interrelationships between the likely abiotic and biotic drivers of productivity.

Figure 1Conceptual model of plausible interaction pathways between environment (elevation, soil) and forest structure as drivers of productivity. Arrows drawn from the two boxes “stand properties” and “soil properties” indicate that pathways from the respective variables (forest stand properties: wood specific gravity, stem density, leaf area index (LAI); soil properties: first two principal component analysis (PCA) axes) were included in the models.


Based on the assumption that precipitation is not a growth-limiting factor in these humid to perhumid study regions, we hypothesized in accordance with the existing literature that (i) temperature is the principal productivity-determining factor in the study regions; (ii) soil nutrient availability is another – but secondary – influential abiotic factor acting in TMFs mainly through N availability; and (iii) tree diversity is a less important driver of productivity than abiotic factors.

2 Methods

2.1 Study transects

The study was conducted in permanent plots in old-growth montane forests established along two elevation transects on the eastern slope of the Ecuadorian Andes in the provinces Napo (Napo transect), Loja and Zamora-Chinchipe (Loja transect) (Fig. A1). The chosen plot area was 20 m × 20 m (400 m2), a size that has previously been used in surveys of tropical-forest diversity, because the area is large enough to give representative information on stand structure and species composition, while the patch size in tropical montane old-growth forests is small enough to guarantee sufficient stand structural and edaphic homogeneity on the plot level. All plots were placed in forest patches without signs of recent disturbance and with sufficient distance to special microhabitats such as ravines and ridges. By selecting small-sized plots solely in old-growth portions of the forest, we attempted to reduce the variation in structural and biological parameters resulting from the mosaic of different successional stages typically existing in natural forests.

The climate in both transect regions is perhumid throughout the year, with a mean annual precipitation > 2000 mm at all elevations and lack of a distinct dry season (Bendix et al., 2008; Salazar et al., 2015). All study sites are characterized by aridity index values > 0.74 after Trabucco and Zomer (2009).

Soil properties vary with elevation and differ between the two transects. The Loja transect is characterized by relatively nutrient-poor soils, mostly on metamorphic schists and sandstones, with a slightly better nutrient supply at lower elevations and in valleys and less favorable supply of N and P on upper slopes and at higher elevations (Wolf et al., 2011; Werner and Homeier, 2015). In contrast, the Napo transect has more fertile soils, which developed on volcanic deposits or limestone. In this transect, the N mineralization rate decreases with elevation, while plant-available P increases (Unger et al., 2010). The availability of the five plant macronutrients N, P, Ca, K and Mg was analyzed in the soil of all study plots together with soil acidity in order to cover the physiologically most meaningful soil chemical properties. For characterizing N availability, incubation experiments to determine the net release of NH4+ and NO3- through mineralization and nitrification (Nmin) were conducted in the topsoil. In addition, the bulk soil C/N ratio was determined in the topsoil for characterizing the decomposability of soil organic matter and thus gives an independent measure of potential N supply (Pastor et al., 1984). The increase in organic-layer depth with elevation in both transects, accompanied by greater C/N ratios and lower pH, is related to higher soil-organic-matter content and indicates in our study areas a decrease in nutrient availability as a result of reduced organic-matter turnover at lower temperatures. P availability was estimated as resin-exchangeable P, i.e., the exchange of phosphate ions by anion exchange resins, which may give a minimum estimate of plant-available P (Pav). The plant availability of Ca, K and Mg (CaKMgex) and also of potentially toxic Al3+ was quantified by salt exchange (0.2 N BaCl2 solution), applying a standard protocol for the chemical analysis of forest soils (for analytical details see Unger et al., 2010; Wolf et al., 2011). Soil pH was measured in soil suspended in 1 M KCl. Since most fine roots are concentrated in the topsoil (Soethe et al., 2006), where the vegetation effect on the soil is also largest, the statistical analyses focus on the chemistry of the 0–10 cm layer and its importance for the vegetation (see Supplement S2 for the soil chemical properties of the plots).

2.2 Plot inventory and forest productivity

The Loja transect consists of 54 plots that are equally spread over three chosen elevation levels (18 plots each at  1000,  2000 and  3000 m a.s.l.) (see Supplement Sect. S1 for the stand properties of the plots in both transects). From the Napo transect, we included data from 66 plots that were located in a comparable elevation range (960–3100 m a.s.l.). In both transects, all trees with a diameter at breast height (dbh)  10 cm in the 400 m2 plot were recorded with their dbh and identified according to species. From unknown tree species, specimens were collected for later identification at Ecuadorian herbaria (HUTPL, LOJA, QCA, QCNE) or by international specialists for difficult groups.

Data on AGB, annual AGB increment (coarse-wood production, WP) and WSG (wood specific gravity) of all trees  10 cm dbh in the plots were available from earlier studies (Loja transect: Wolf et al., 2011; Leuschner et al., 2013; Napo transect: Unger et al., 2012; Kessler et al., 2014). AGB and WP of all trees  10 cm dbh in a plot were estimated as the sum of the biomass of all tree individuals and the corresponding increment of these trees in a second stem diameter inventory after 1–5 years. For quantifying AGB, we applied the allometric equation of Chave et al. (2005) for tropical wet forests. WP was calculated as the difference between the two AGB estimates on an annual basis. For the Loja transect we additionally calculated net aboveground productivity (NPPa) as the sum of WP and annual fine-litter production that was recorded in litter traps (data available from Wallis et al., 2019).

Leaf area index (LAI) was estimated by synchronous measurement with two LAI 2000 plant canopy analyzers (LI-COR Inc., Lincoln, NE, USA) that were operated in the remote mode during periods of overcast sky. Simultaneous readings were taken below the canopy at 2 m height above the ground and in nearby open areas (“above-canopy” reading) to record incoming radiation. Data for the Napo transect were taken from Unger et al. (2013). For the Loja transect, the LAI measurements were conducted in 2011, adopting the methods described in detail in Unger et al. (2013).

2.3 Data analysis

Tree diversity was calculated according to the individual-based rarefaction method (Gotelli and Colwell, 2001) as the number of species (stems with dbh  10 cm) expected in a random sample of 14 trees in a 400 m2 plot (14 being the smallest number of tree individuals recorded in our plots). Linear regression analyses were applied to identify significant relationships between AGB, WP, tree diversity, LAI, stem density, WSG and soil properties as dependent variables and elevation (as a proxy for temperature) as an independent variable.

We used principal component analysis (PCA) to reduce the number of soil variables and to ensure that the subsequent analyses were not affected by the problem of multi-collinearity. We conducted two PCAs using the R package “pcaMethods”, one with the merged data from both transects and another one for the Loja transect. The following 10 soil variables were included in the PCAs: organic-layer depth; mineral soil pH (KCl); exchangeable K, Mg, Ca and Al content; resin-exchangeable P content; C/N ratio; and topsoil N mineralization and nitrification rate. The plot scores on the first two resulting principal components (PC 1 and PC 2) were included in the subsequent analyses.

Structural equation modeling was used for identifying the direct influence of abiotic factors (elevation, soil principal components) on forest stand properties and tree diversity and to assess the direct and indirect influences of these factors on productivity. For each initial structural equation model (SEM), we included the two principal soil components (PC1 and PC2) from the respective PCA and the factors elevation, WSG and tree diversity to predict AGB and stand productivity (WP or NPPa). WSG was selected from the stand properties as it showed a stronger correlation to stand productivity than the other measured variables (stem density and LAI).

Based on the existing knowledge about the dependency of productivity and AGB on environmental and stand properties and assumed diversity effects on productivity, we developed an initial model of interaction paths between the environmental parameters (elevation, soil PC 1 and PC 2), WSG and tree diversity, and the target variables (AGB and WP or NPPa) (Fig. 1). We then fitted two different SEMs, one for the complete set of 83 plots from both transects (where all information was available: 54 Loja plots and 29 Napo plots) and another one for the Loja transect (54 plots), where NPPa (and not only WP) data are available.

IBM SPSS AMOS 24 software (Arbuckle, 2016) was used to fit the normalized data to the hypothesized path model and to determine path coefficients using the maximum likelihood method. We assessed the goodness of model fit using the χ2 value, the associated p value, Akaike information criterion (AIC), RMSEA (root mean square error of approximation) and CFI (comparative fit index). Since we used SEM in an explorative way, the original model has been subject to modification. We iteratively removed insignificant paths to test whether incorporation of those paths in the model significantly increased the χ2 value and CFI and reduced the AIC and the RMSEA of the model. For all dependent variables, we calculated R2 values that indicate the proportion of variance explained by the model.

3 Results

3.1 Patterns of tree diversity and productivity

The two transects differed considerably with respect to floristic composition. The most important tree families (in order of abundance) were Lauraceae, Melastomataceae, Cunoniaceae, Moraceae and Clusiaceae in the Loja transect and Euphorbiaceae, Fabaceae, Lauraceae, Moraceae and Rubiaceae in the Napo transect. The number of stems per plot varied between 14 and 61 and the number of tree species between 6 and 27 (Supplement Sect. S1). Tree species richness showed no significant decline with elevation (Fig. 2a). AGB decreased on average by 50–90 Mg ha−1 per kilometer elevation in both transects (Fig. 2b). Biomass was between 200 Mg ha−1 (lower elevations) and 100 Mg ha−1 (upper elevations) higher in the Napo transect on more fertile soil than in the Loja transect.

Figure 2Variation in tree diversity (a), aboveground biomass (b), wood production (c), LAI (d), wood specific gravity (e) and stem density (f) along the two studied forest transects (red dots: Loja transect; blue dots: Napo transect). The number of plots available for a given parameter and transect varies as not all parameters were measured in all plots (a: Loja 54 plots, Napo 64 plots; b: 54, 66; c: 54, 29; d: 53, 60; e: 54, 66; f: 54, 66). Data presented in (b), (c), (e) and (f) were compiled from Unger et al. (2012), Leuschner et al. (2013) and Kessler et al. (2014); data from the Napo transect in (d) are from Unger et al. (2013).


The AGB decrease was linked to a decline in WP between 1000 and 3000 m by about 1.3–1.5 Mg ha−1 yr−1 per kilometer elevation, with a generally lower productivity (by about 1.0–1.5 Mg ha−1 yr−1) in the Loja transect compared to the Napo transect (Fig. 2c). LAI decreased with elevation in both transects by about 1 m2 m−2 per kilometer elevation, with slightly higher LAI values in the Loja transect (Fig. 2d). WSG showed no systematic variation with elevation (Fig. 2e). The expected stem density increase with elevation was only observed in the Loja transect (Fig. 2f).

3.2 Identification of principal soil factors

In the two data sets (both transects merged, Loja transect only), the first two PCA axes explained 63 % and 61 %, respectively, of the variation in soil factors across all plots. In the merged data set (Figs. 3 and B1 and Table B1), the first axis represented mainly the basic cations Mg, Ca and K; soil pH; and organic-layer depth; the second axis stood for N supply and Pav. In the PCA of the Loja transect (Fig. B2 and Table B2), soil parameters were represented similarly, but Pav had a higher loading on the first axis. In general, the soils of the Napo transect were characterized by higher fertility, as indicated by on average markedly higher topsoil net N mineralization rates and higher P availability in the mineral topsoil (see Supplement Sect. S2).

Figure 3Variation in soil properties along the two forest transects (red dots: Loja transect; blue dots: Napo transect). Shown are the (a) first and (b) second principal components of the soil PCA (see Fig. B1 and Table B1).


3.3 Factors determining wood production in the merged data set

In the merged data set (83 plots; Fig. 4, Table 1), aboveground tree biomass (AGB) was positively influenced by soil fertility (standardized total effects of PC1 and PC2 of 0.30 and 0.45, respectively) and tree diversity (0.23). The three most important determinants of WP were elevation (negative influence, standardized direct effect of 0.58), soil PC2 (0.29) and WSG (0.25).

The direct negative effect of elevation on WP was prominent, and this influence was enhanced through indirect relationships via negative effects of elevation on soil properties and tree diversity (standardized indirect effect of 0.18). However, soil factors were also important drivers of wood production, either directly (in case of PC2, which was related to N mineralization; 0.29) or indirectly through a soil N effect (PC2) on AGB and diversity (together 0.10) and effects of base cation availability and pH (PC1) on WSG and AGB. Tree diversity itself had only a relatively weak indirect positive effect on WP through its influence on AGB (standardized total effect of 0.05).

3.4 Factors determining NPPa in the Loja transect

The model for aboveground net primary production (i.e., WP plus fine-litter production), built only with the Loja transect data, identified elevation and WSG as the main direct determinants of NPPa (both effects negative; standardized direct effects of 0.67 and 0.39, respectively; Fig. 5, Table 2). Soil properties were also influential (mainly PC1 with N availability) but only indirectly through a negative effect on WSG and a positive effect on diversity (standardized total effect of 0.30). Tree diversity had a significant positive but relatively weak direct effect on NPPa (0.18), which was enhanced by an indirect influence via AGB (standardized total effect of 0.27).

4 Discussion

In a global perspective, forest productivity increases with temperature in cooler regions and with precipitation in drier regions (Reich and Bolstad, 2001; Schuur, 2003). In the absence of temperature and precipitation measurements and soil moisture data from the plots, we used elevation as a proxy for the thermal conditions. We further assumed that soil moisture rarely constrains forest productivity in both study regions as precipitation totals range well above 2000 mm yr−1 throughout the Napo and Loja transects (Bendix et al., 2008; Salazar et al., 2015), and soil moisture is usually high in Andean TMFs (Zimmermann et al., 2010). Our assumption is supported by the highly significant negative influence of elevation on WP, suggesting that the temperature decrease with elevation was physiologically much more important than a possible water shortage effect on productivity. Both the SEMs and correlation analyses show that tropical-forest productivity largely decreases (by 1.3–1.5 Mg ha−1 yr−1 km−1) with a concurrent decrease in mean annual temperature from  20 to  10 C in the two transects, confirming earlier observations (Aiba et al., 2005; Cleveland et al., 2011). Reduced temperatures may negatively affect tree growth through a variety of direct and indirect influences on carbon gain, carbon allocation and nutrient acquisition. Photosynthesis measurements revealed community means of light-saturated net photosynthesis rates (Amax) that were 36 % and 18 % lower at 3000 m compared to the community means at 2000 and 1000 m, respectively (Wittich et al., 2012). Reduced canopy carbon gain, in combination with the elevational LAI decrease by 1 m2 m−2 km−1, might largely explain the productivity decline from 1000 to 3000 m in the Ecuadorian Andes (Leuschner et al., 2013). In apparent contradiction, the tree-based carbon balance model of Fyllas et al. (2017) suggests for the Peruvian Andes that the apparent temperature effect on photosynthesis and productivity in the 3300 m elevation gradient is manifested through leaf trait variation associated with the elevational species turnover rather than through a temperature dependency of photosynthesis itself. In fact, light-saturated net photosynthesis (Amax) did not change between 250 and 3500 m a.s.l. in this transect (Fyllas et al., 2017). These results contradict a pantropical analysis (ca. 170 tree species, 18 sites) showing that Amax decreases on average by 1.3 µmol m−2 s−1 per kilometer elevation in tropical mountains (Wittich et al., 2012). We assume that local differences in soil fertility and climatic factors such as cloudiness or very high precipitation, causing temporal soil anoxia, may have a large influence on elevational patterns in foliar nutrient content and photosynthetic capacity. Moreover, the Amax study of Wittich et al. (2012) with 40 Andean tree species from 21 families points at a considerable influence of phylogeny on photosynthetic capacity. Species turnover between elevations and also between different locations might therefore explain some of the differences in photosynthetic carbon gain and productivity independently of temperature differences.

Low-temperature effects on nutrient supply and acquisition processes, which enhance nutrient deficiency, could well be additional factors causing productivity reductions and increased belowground C allocation at higher elevations (Graefe et al., 2008; Moser et al., 2011; Leuschner et al., 2013). Indeed, pronounced leaf trait changes were observed in the Loja transect from larger, thinner and N-richer leaves at 1000 m to smaller, thicker and N-poorer leaves at 3000 m (Homeier et al., unpublished), suggesting that temperature may act on the forest carbon cycle partly through N and P availability as an environmental filter, sorting tree species according to their leaf traits, from more acquisitive to more conservative, at higher elevations.

Figure 4Final model for wood production (Loja transect and Napo transect: 83 plots). Structural equation model (chi-square = 8.4; 9 degrees of freedom; p= 0.50; AIC = 60.4; RMSEA <0.01; CFI = 1.00) with standardized path coefficients. The size of the arrows is proportional to the strength of the paths; their significance is indicated by asterisks (* P<0.05, ** P<0.01, *** P<0.001). Blue arrows indicate positive and red negative estimates. AGB: aboveground biomass; WSG: wood specific gravity; PC1, PC2: first two axes of the soil PCA (Fig. B1 and Table B1).


Table 1Standardized direct, indirect and total effects of elevation, various stand structural parameters and soil chemical parameters on coarse-wood production (WP) and aboveground biomass (AGB) according to the SEM analysis in Fig. 4. Standardized path coefficients are shown.

Download Print Version | Download XLSX

Temporal water limitation, if it were relevant in the two regions, should be more prominent at lower elevations with lower rainfall, which should lead to a positive and not a negative elevation effect on productivity. We cannot exclude, however, that high soil moisture in combination with temporal hypoxia in the rhizosphere of the high-elevation stands is contributing to the productivity decline toward the timberline (Moser et al., 2011).

Other climatic factors with a potential impact on productivity such as solar radiation, air humidity and leaf wetness were not considered in our analysis due to lack of suitable local data. Solar radiation typically reaches a minimum at montane or upper montane elevation due to cloud immersion, which can reduce carbon gain (Malhi et al., 2017). Including these factors in the analysis may lead to somewhat different conclusions on the key determinants of forest productivity (Finegan et al., 2015; Fyllas et al., 2017; Malhi et al., 2017).

Our study in two regions of the northern Andes is the first worldwide to analyze the dependence of tropical-montane-forest production on the availability of all five quantitatively most important plant nutrients (N, P, Ca, K, Mg). P availability has been reported to limit or co-limit forest productivity at many tropical lowland sites (Tanner et al., 1998; Paoli and Curran, 2007; Cleveland et al., 2011; Wright et al., 2011; Dalling et al., 2016b), whereas N seems to limit tropical-forest productivity more often at higher elevations (Tanner et al., 1998; Benner et al., 2010; Wolf et al., 2011; Dalling et al., 2016a). The prominent direct and indirect effects of PC1 and PC2 on productivity suggest that all three nutrient components (N, P and base cations) play roles in the limitation of forest productivity in the study regions.

Figure 5Final model for NPPa (Loja transect: 54 plots). Structural equation model (chi-square = 9.9; 10 degrees of freedom; p= 0.45; AIC = 59.9; RMSEA <0.01; CFI = 0.99) with standardized path coefficients. The size of the arrows is proportional to the strength of the paths; their significance is indicated by asterisks (* P<0.05, ** P<0.01, *** P<0.001). Blue arrows indicate positive and red negative estimates. AGB: aboveground biomass; WSG: wood specific gravity; PC1, PC2: first two axes of the soil PCA (see Fig. B2 and Table B2).


Table 2Standardized direct, indirect and total effects of elevation, various stand structural parameters and soil chemical parameters on NPPa and aboveground biomass (AGB) (Loja transect) according to the SEM analysis in Fig. 5. Standardized path coefficients are shown.

Download Print Version | Download XLSX

A high net N mineralization rate positively influenced WP in the merged data set, but the indirect effect through AGB dominated. The large difference in N mineralization rates between the two transects, with high values in the predominantly volcanic soils in the Napo region but low values in the non-volcanic soils of the Loja region, probably also explains the considerably higher AGB and consequently WP in the Napo transect as compared to the Loja transect (Fig. 2). The trees in the less fertile Loja transect most likely allocate a higher proportion of their photosynthates to root production and to root symbionts (Vicca et al., 2012; Doughty et al., 2018), which must reduce wood production. The availability of base cations and soil acidity influenced productivity in both models only indirectly, mainly through an effect on WSG (higher wood density on acidic, base-poor soils).

The PCAs also showed that the three main nutrient categories (P, N, base cations) varied more or less independently from each other across the plots. Moreover, their relative importance varied between the Loja and Napo data sets, suggesting that soil nutrient limitation of forest productivity depends more on region, bedrock type and soil age than on elevation. The different paths from soil nutrients to productivity and the weak correlations among Pav, Nmin and base cations across our data set imply that no single nutrient element can serve as a good indicator of overall nutrient availability or soil fertility in the study regions. This is also valid for pH, which showed the tightest relation to CaKMgex (Loja transect) or organic-layer depth (merged transects) but weaker or no association to Pav, Nmin and soil C/N ratio. Soil fertility assessments in Andean montane forests should consider N, P and also basic cations when forest productivity shall be related to soil fertility. The widely held assumption that net primary production in tropical lowland forests is primarily limited by low P availability and TMF productivity predominantly by low N availability (Vitousek and Sanford, 1986; Tanner et al., 1998; Fisher et al., 2013) may thus require adaptation.

Various studies have identified WSG as a good predictor of the diameter increment of tropical trees, and high WSG has been linked to low aboveground productivity in tropical lowland forests (Malhi et al., 2004; Poorter et al., 2008; Finegan et al., 2015). In the Ecuadorian Andes, the negative effect of average WSG on productivity was as strong as the positive biomass (AGB) influence, suggesting that negative edaphic effects on productivity, notably high soil acidity and low availability of basic cations, are exerted in part indirectly through an increase in wood specific gravity (Unger et al., 2012). A recent study in Southeast Asian tropical forests showed that WSG has a direct effect neither on biomass production (Kotowska et al., 2021) nor on tree water consumption even though harder wood tends to be associated with lower growth rates (Muller-Landau, 2004; Hoeber et al., 2014). Yet, WSG is known to be associated with most structural and functional wood properties (Chave et al., 2009), and it may also be related to anatomical attributes such as pit membrane characteristics that influence sap flux density, which itself is positively related to productivity (Kotowska et al., 2021). In addition, our data support a causal link from higher N availability through lower WSG to higher productivity, which is suggested by N fertilization experiments in conifers (Jozsa and Brix, 1989) and a negative correlation between wood density and soil fertility in continent-wide wood property analyses (Chave et al., 2009). N fertilization experiments in Andean TMFs support the positive effect of N on productivity (Homeier et al., 2012; Fisher et al., 2013), but they were not designed to prove a WSG-mediated mechanism.

Unexpected is the insignificant effect of LAI on WP in both transects, visible in the exclusion of leaf area from the models. This is difficult to reconcile with the fact that canopy carbon gain depends on the amount of intercepted light and that the LAI–productivity relationship is usually tight in forests (Reich, 2012). A possible explanation could be that LAI measurements with the LAI2000 canopy analyzer are known to lead to both underestimation where leaf clumping is dominant and overestimation where stem density is high. This may have blurred LAI differences between stands differing in productivity.

Tree diversity apparently has only a small influence on productivity in the studied tropical montane forests even though tree species numbers varied largely across the plots (6–27 species per 0.04 ha), and plot size was small so that the effect should be more visible. A positive diversity effect apparently was confined to the component fine-litter production (in the Loja transect) as the effect became more visible when NPPa (which includes fine-litter production) was considered and not wood production alone (Fig. 4). A possible explanation of the weak diversity effect is that the species richness in our plots may have been too high to show a positive DPR as diversity effects on productivity generally are more pronounced at lower species numbers (Paquette and Messier, 2011; Tilman et al., 1997). It would be interesting to see how tree species richness translates into functional diversity in our two transects and if functional diversity is already saturating at intermediate levels of tree diversity (as a result of increasing functional redundancy). That would explain the weak diversity effect. But functional trait data are only available for a small fraction of the recorded tree species.

Our results support our third hypothesis and fit the controversial picture that emerges from the existing observational studies on the DPR in tropical forests. Chisholm et al. (2013) systematically explored the role of diversity for productivity in natural tropical forests and found only a small absolute effect: when controlling for stem density effects in the 0.04 ha plots, a doubling of species richness corresponded in their study to only a 5 % increase in WP on average. At a larger scale (0.25, 1 ha), results were mixed, and negative relationships became more common. This points at a negligible role of tree diversity for productivity in tropical-forest landscapes, in line with our findings. These results from tropical forests match findings from more than 100 000 forest plots in the United States, where the predominant DPR was concave and negative in humid climates and primarily non-significant in harsh climates (Fei et al., 2018). Similar to the DPR, a significant positive diversity–aboveground biomass relationship was detectable in tropical lowland forests only in small plots (0.04 ha), while it disappeared at a larger scale (1 ha plots) (Sullivan et al., 2017). Moreover this study showed that the diversity effect on AGB was small: doubling species richness at 0.04 ha increased biomass only by 6.9 %.

An additional important result of our study is the high variation in stand properties and WP across the plots of an elevation level (Fig. 2). Environmental variation related to topography can be a strong driver of forest structure and tree species composition, and it influences productivity and nutrient cycling through changes in traits, as has been shown for several tropical lowland forests (Fortunel et al., 2018; Jucker et al., 2018) and also for the TMF of the Loja transect (Werner and Homeier, 2015; Pierick et al., 2021).

Our analysis of multiple drivers of productivity further shows that a diversity effect, even if significant, is only of minor importance when other productivity-influencing factors are also included in the analysis. In our case, diversity ranked behind effects of elevation (temperature), soil nutrient availability and wood specific gravity. For understanding the controls of forest productivity, it is more important to study the influences of environmental factors and stand structural and functional properties in more detail.

5 Conclusions

From the analysis of more than 80 old-growth forest plots in two highly diverse Andean regions with large geological and topographic heterogeneity, we conclude that the main determinants of aboveground productivity in tropical montane forests are elevation (primarily as a proxy for temperature), nutrient supply and stand structural properties, given non-limiting moisture conditions. Nutrient availability seems to affect productivity not only directly but also indirectly through stand properties, while tree diversity has only a small influence on productivity. Tree species turnover and associated changes in functional traits are far more important for productivity than changes in species diversity. Functional biodiversity research has to proceed from the search for the significance of diversity effects, which has dominated biodiversity–ecosystem functioning (BEF) research in the past, to an assessment of the importance of diversity effects for productivity and other ecosystem processes in order to better understand the functional role of diversity. A deeper understanding of foliar and root trait variation with elevation and among the species of a community will help to achieve a more mechanistic understanding of tropical-forest productivity change along environmental gradients. Our findings are of practical relevance as the analysis addresses spatial scales relevant for old-growth forest conservation and management.

Appendix A: Study sites and plots

Figure A1Map showing the location of the two study sites in Ecuador.

Appendix B: PCAs of soil properties

Figure B1Principal component analysis of the soil properties of 83 study plots (Loja transect and Napo transect); PC1 and PC2 explained 37.4 % and 26.1 % of variance, respectively.


Table B1Factor loadings of the variables in the PCA (Fig. B1).

Download Print Version | Download XLSX

Figure B2Principal component analysis of the soil properties of 54 study plots of the Loja transect; PC1 and PC2 explained 46.9 % and 14.6 % of variance, respectively.


Table B2Factor loadings of the variables in the PCA (Fig. B2).

Download Print Version | Download XLSX

Data availability

All data sets used are included in the Supplement.


Supplement S1 contains a table with the stand parameters of the study plots. Supplement S2 contains a table with the soil parameters of the study plots. The supplement related to this article is available online at:

Author contributions

JH and CL designed the study, JH performed the research and analyzed the data, and JH and CL wrote the manuscript.

Competing interests

The authors declare that they have no conflict of interest.


We thank the universities in Loja (UTPL, UNL) and in Quito (PUCE) for continuous support during our field studies. We thank the Ministerio de Ambiente del Ecuador for granting the research permits. We further acknowledge the assistance of Nixon Cumbicus, Miguel-Angel Chinchero, Jaime Peña, Roman Link, Malte Unger and Katrin Mikolajewski during fieldwork and lab work. We thank the two anonymous reviewers for their constructive criticism of an earlier draft of the paper.

Financial support

This research has been supported by the Deutsche Forschungsgemeinschaft (grant nos. Ho3296/2, Ho3296/4 and Le762/10) and the Bundesministerium für Bildung und Forschung (grant no. Pro Benefit).

This open-access publication was funded
by the University of Göttingen.

Review statement

This paper was edited by Sara Vicca and reviewed by two anonymous referees.


Aragão, L. E. O. C., Malhi, Y., Metcalfe, D. B., Silva-Espejo, J. E., Jiménez, E., Navarrete, D., Almeida, S., Costa, A. C. L., Salinas, N., Phillips, O. L., Anderson, L. O., Alvarez, E., Baker, T. R., Goncalvez, P. H., Huamán-Ovalle, J., Mamani-Solórzano, M., Meir, P., Monteagudo, A., Patiño, S., Peñuela, M. C., Prieto, A., Quesada, C. A., Rozas-Dávila, A., Rudas, A., Silva Jr., J. A., and Vásquez, R.: Above- and below-ground net primary productivity across ten Amazonian forests on contrasting soils, Biogeosciences, 6, 2759–2778,, 2009. 

Arbuckle, J. S.: IBM SPSS Amos 24 User's Guide, IBM, Chicago, USA, 2016. 

Aiba, S., Takyu, M., and Kitayama, K.: Dynamics, productivity and species richness of tropical rainforests along elevational and edaphic gradients on Mount Kinabalu, Borneo, Ecol. Res., 20, 279–286, 2005. 

Ashton, P. S.: Floristic zonation of tree communities on wet tropical mountains revisited, Perspect. Plant Ecol., 6, 87–104,, 2003. 

Banin, L., Lewis, S. L., Lopez-Gonzalez, G., Baker, T. R., Quesada, C. A., Chao, K. J., Burslem, D. F. R. P., Nilus, R., Salim, K. A., Keeling, H. C., Tan, S., Davies, S. J., Mendoza, A. M., Vásquez, R., Lloyd, J., Neill, D. A., Pitman, N., and Phillips, O. L.: Tropical forest wood production: a cross-continental comparison, J. Ecol., 104, 1025–1037,, 2014. 

Bendix, J., Rollenbeck, R., Fabian, P., Emck, P., Richter, M., and Beck, E.: Climate variability, in: Gradients in a tropical mountain ecosystem of Ecuador, edited by: Beck, E., Kottke, I., Bendix, J., Makeschin, F., and Mosandl, R., Springer, Berlin and Heidelberg, Germany, 281–290, 2008. 

Benner, J., Vitousek, P. M., and Ostertag, R.: Nutrient cycling and nutrient limitation in tropical montane cloud forests, in: Tropical montane cloud forests: science for conservation and management, edited by: Bruijnzeel, L. A., Scatena, F. N., and Hamilton, L., Cambidge University Press, Cambidge, UK, 90–100, 2010. 

Bruijnzeel, L. A. and Veneklaas, E. J.: Climatic conditions and tropical montane forest productivity: the fog has not lifted yet, Ecology, 79, 3–9, 1998. 

Bruijnzeel, L. A., Scatena, F. N., and Hamilton, L. S.: Tropical montane cloud forests: science for conservation and management, Cambridge University Press, Cambridge, UK, 2010. 

Chave, J., Andalo, C., Brown, S., Cairns, M. A., Chambers, J. Q., Eamus, D., Fölster, H., Fromard, F., Higuchi, N., Kira, T., Lescure, J. P., Nelson, B. W., Ogawa, H., Puig, H., Riéra, B., and Yamakura, T.: Tree allometry and improved estimation of carbon stocks and balance in tropical forests, Oecologia, 145, 87–99,, 2005. 

Chave, J., Coomes, D., Jansen, S., Lewis, S. L., Swenson, N. G., and Zanne, A. E.: Towards a worldwide wood economics spectrum, Ecol. Lett., 12, 351–366,, 2009. 

Chen, Y., Wright, S. J., Muller-Landau, H. C., Hubbell, S. P., Wang, Y., and Yu, S.: Positive effects of neighborhood complementarity on tree growth in a Neotropical forest, Ecology, 97, 776–785,, 2016. 

Chisholm, R. A., Muller-Landau, H. C., Abdul Rahman, K., Bebber, D. P., Bin, Y., Bohlman, S. A., Bourg, N. A., Brinks, J., Bunyavejchewin, S., Butt, N., Cao, H., Cao, M., Cárdenas, D., Chang, L.-W., Chiang, J.-M., Chuyong, G., Condit, R., Dattaraja, H. S., Davies, S., Duque, A., Fletcher, C., Gunatilleke, N., Gunatilleke, S., Hao, Z., Harrison, R. D., Howe, R., Hsieh, C.-F., Hubbell, S. P., Itoh, A., Kenfack, D., Kiratiprayoon, S., Larson, A. J., Lian, J., Lin, D., Liu, H., Lutz, J. A., Ma, K., Malhi, Y., McMahon, S., McShea, W., Meegaskumbura, M., Mohd. Razman, S., Morecroft, M. D., Nytch, C. J., Oliveira, A., Parker, G. G., Pulla, S., Punchi-Manage, R., Romero-Saltos, H., Sang, W., Schurman, J., Su, S.-H., Sukumar, R., Sun, I.-F., Suresh, H. S., Tan, S., Thomas, D., Thomas, S., Thompson, J., Valencia, R., Wolf, A., Yap, S., Ye, W., Yuan, Z., and Zimmerman, J. K.: Scale-dependent relationships between tree species richness and ecosystem function in forests, J. Ecol., 101, 1214–1224,, 2013. 

Cleveland, C. C., Townsend, A. R., Taylor, P., Alvarez-Clare, S., Bustamante, M. M., Chuyong, G., Dobrowski, S. Z., Grierson, P., Harms, K. E., Houlton, B. Z., Marklein, A., Parton, W., Porder, S., Reed, S. C., Sierra, C. A., Silver, W. L., Tanner, E. V. J., and Wieder, W. R.: Relationships among net primary productivity, nutrients and climate in tropical rain forest: a pan-tropical analysis, Ecol. Lett., 14, 939–947,, 2011. 

Dalling, J. W., Heineman, K., Lopez, O. R., Wright, S. J., and Turner, B. L.: Nutrient availability in tropical rain forests: the paradigm of phosphorus limitation, in: Tropical Tree Physiology, edited by: Goldstein, G. and Santiago, L. S., Springer, Cham, UK, 261–273, 2016a. 

Dalling, J. W., Heineman, K., González, G., and Ostertag, R.: Geographic, environmental and biotic sources of variation in the nutrient relations of tropical montane forests, J. Trop. Ecol., 32, 368–383,, 2016b. 

Doughty, C. E., Goldsmith, G. R., Raab, N., Girardin, C. A. J., Farfan-Amezquita, F., Huaraca-Huasco, W., Silva-Espejo, J. E., Araujo-Murakami, A., da Costa, A. C. L., Rocha, W., Galbraith, D., Meir, P., Metcalfe, D. B., and Malhi, Y.: What controls variation in carbon use efficiency among Amazonian tropical forests?, Biotropica, 50, 16–25,, 2018. 

Fahey, T. J., Sherman, R. E., and Tanner, E.: Tropical montane cloud forest: environmental drivers of vegetation structure and ecosystem function, J. Trop. Ecol., 32 355–367,, 2015. 

Fei, S., Jo, I., Guo, Q., Wardle, D. A., Fang, J., Chen, A., Oswalt, C. M., and Brockerhoff, E. G.: Impacts of climate on the biodiversity-productivity relationship in natural forests, Nat. Commun., 9, 1–7,, 2018. 

Finegan, B., Peña-Claros, M., de Oliveira, A., Ascarrunz, N., Bret-Harte, M. S., Carreño-Rocabado, G., Casanoves, F., Díaz, S., Eguiguren Velepucha, P., Fernandez, F., Licona, J. C., Lorenzo, L., Salgado Negret, B., Vaz, M., and Poorter, L.: Does functional trait diversity predict above-ground biomass and productivity of tropical forests? Testing three alternative hypotheses, J. Ecol., 103, 191–201,, 2015. 

Fisher, J. B., Malhi, Y., Torres, I. C., Metcalfe, D. B., van de Weg, M. J., Meir, P., Silva-Espejo, J. E., and Huaraca Huasco, W.: Nutrient limitation in rainforests and cloud forests along a 3000-m elevation gradient in the Peruvian Andes, Oecologia, 172, 889–902,, 2013. 

Fortunel, C., Lasky, J. R., Uriarte, M., Valencia, R., Wright, S. J., Garwood, N. C., and Kraft, N. J. B.: Topography and neighborhood crowding can interact to shape species growth and distribution in a diverse Amazonian forest, Ecology, 99, 2272–2283,, 2018. 

Fyllas, N. M., Bentley, L. P., Shenkin, A., Asner, G. P., Atkin, O. K., Díaz, S., Enquist, B. J., Farfan-Rios, W., Gloor, E., Guerrieri, R., Huasco, W. H., Ishida, Y., Martin, R. E., Meir, P., Phillips, O., Salinas, N., Silman, M., Weerasinghe, L. K., Zaragoza-Castells, J., and Malhi, Y.: Solar radiation and functional traits explain the decline of forest primary productivity along a tropical elevation gradient, Ecol. Lett., 20, 730–740,, 2017. 

Girardin, C. A. J., Malhi, Y., Aragão, L. E. O. C., Mamani, M., Huaraca Huasco, W., Durand, L., Feeley, K. J., Rapp, J., Silva-Espejo, J. E., Silman, M., Salinas, N., and Whittaker, R. J.: Net primary productivity allocation and cycling of carbonalong a tropical forest elevational transect in the Peruvian An-des, Glob. Change Biol., 16, 3176–3192,, 2010. 

Girardin, C. A. J., Malhi, Y., Feeley, K. J., Rapp, J. M., Silman, M. R., Meir, P., Huaraca Huarasco, W., Salinas, N., Mamani, M., Silva-Espejo, J. E., Garcia Cabrera, K., Farfan Rios, W., Metcalfe, D. B., Doughty, C. E., and Aragão, L.: Seasonality of above-ground net primary productivity along an Andean altitudinal transect in Peru, J. Trop. Ecol., 30, 503–519,, 2014. 

Gotelli, N. J. and Colwell, R. K.: Quantifying biodiversity: procedures and pitfalls in the measurement and comparison of species richness, Ecol. Lett., 4, 379–391,, 2001. 

Graefe, S., Hertel, D., and Leuschner, C.: Estimating Fine Root Turnover in Tropical Forests along an Elevational Transect using Minirhizotrons, Biotropica, 40, 536–542,, 2008. 

Hoeber, S., Leuschner, C., Köhler, L., Arias-Aguilar, D., and Schuldt, B.: The importance of hydraulic conductivity and wood density to growth performance in eight tree species from a tropical semi-dry climate, Forest Ecol. Manag., 300, 126–136,, 2014. 

Hofhansl, F., Schnecker, J., Singer, G., and Wanek, W.: New insights into mechanisms driving carbon allocation in tropical forests, New Phytol., 205, 137–146,, 2015. 

Homeier, J., Breckle, S.-W., Günter, S., Rollenbeck, R. T., and Leuschner, C.: Tree Diversity, Forest Structure and Productivity along Altitudinal and Topographical Gradients in a Species-Rich Ecuadorian Montane Rain Forest, Biotropica, 42, 140–148,, 2010. 

Homeier, J., Hertel, D., Camenzind, T., Cumbicus, N., Maraun, M., Martinson, G. O., Poma, L. N., Rillig, M. C., Sandmann, D., Scheu, S., Veldkamp, E., Wilcke, W., Wullaert, H., and Leuschner, C.: Tropical Andean forests are highly susceptible to nutrient inputs – Rapid effects of experimental N and P addition to an Ecuadorian montane forest, PLOS ONE, 7, e47128,, 2012. 

Huang, Y., Chen, Y., Castro-Izaguirre, N., Baruffol, M., Brezzi, M., Lang, A., Li, Y., Härdtle, W., von Oheimb, G., Yang, X., Liu, X., Pei, K., Both, S., Yang, B., Eichenberg, D., Assmann, T., Bauhus, J., Behrens, T., Buscot, F., Chen, X.-Y., Chesters, D., Ding, B.-Y., Durka, W., Erfmeier, A., Fang, J., Fischer, M., Guo, L.-D., Guo, D., Gutknecht, J. L. M., He, J.-S., He, C.-L., Hector, A., Hönig, L., Hu, R.-Y., Klein, A.-M., Kühn, P., Liang, Y., Li, S., Michalski, S., Scherer-Lorenzen, M., Schmidt, K., Scholten, T., Schuldt, A., Shi, X., Tan, M.-Z., Tang, Z., Trogisch, S., Wang, Z., Welk, E., Wirth, C., Wubet, T., Xiang, W., Yu, M., Yu, X.-D., Zhang, J., Zhang, S., Zhang, N., Zhou, H.-Z., Zhu, C.-D., Zhu, L., Bruelheide, H., Ma, K., Niklaus, P. A., and Schmid, B.: Impacts of species richness on productivity in a large-scale subtropical forest experiment, Science, 362, 80–83,, 2018. 

Jozsa, L. A. and Brix, H.: The effects of fertilization and thinning on wood quality of a 24-year-old Douglas-fir stand, Can. J. Forest Res., 19, 1137–1145,, 1989. 

Jucker, T., Bongalov, B., Burslem, D. F. R. P., Nilus, R., Dalponte, M., Lewis, S. L., Phillips, O. L., Qie, L., and Coomes, D. A.: Topography shapes the structure, composition and function of tropical forest landscapes, Ecol. Lett., 21, 989–1000,, 2018. 

Keeling, H. C. and Phillips, O. L.: The global relationship between forest productivity and biomass, Global Ecol. Biogeogr., 16, 618–631,, 2007. 

Kessler, M., Salazar, L., Homeier, J., and Kluge, J.: Species richness – productivity relationships of tropical terrestrial ferns at regional and local scales, J. Ecol., 102, 1623–1633,, 2014. 

Kitayama, K. and Aiba, S.-I.: Ecosystem structure and productivity of tropical rain forests along altitudinal gradients with contrasting soil phosphorus pools on Mount Kinabalu, Borneo, J. Ecol., 90, 37–51,, 2002. 

Kotowska, M., Link, R., Röll, A., Hertel, D., Hölscher, D., Waite, P.-A., Moser, G., Tjoa, A., Leuschner, C., and Schuldt, B.: Effects of wood hydraulic properties on water use and productivity of tropical rainforest trees, Front. For. Glob. Change, 3, 598759,, 2021. 

Lambers, H. and Oliveira, R. S.: Plant Physiological Ecology, 3rd edition, Springer Nature, Cham, UK, 2019. 

Leuschner, C., Zach, A., Moser, G., Homeier, J., Graefe, S., Hertel, D., Wittich, B., Soethe, N., Iost, S., Röderstein, M., Horna, V., and Wolf, K.: The carbon balance of tropical mountain forests along an altitudinal transect, in: Ecosystem services, biodiversity and environmental change in a tropical mountain ecosystem of South Ecuador, edited by: Bendix, J., Beck, E., Bräuning, A., Makeschin, F., Mosandl, R., Scheu, S., and Wilcke, W., Springer, Berlin and Heidelberg, Germany, 117–139,, 2013. 

Liang, J., Crowther, T. W., Picard, N., Wiser, S., Zhou, M., Alberti, G., Schulze, E.-D., Mcguire, A. D., Bozzato, F., Pretzsch, H., Paquette, A., Hérault, B., Scherer-Lorenzen, M., Barrett, C. B., Glick, H. B., Hengeveld, G. M., Nabuurs, G.-J., Pfautsch, S., Viana, H., Vibrans, A. C., Ammer, C., Schall, P., Verbyla, D., Tchebakova, N., Fischer, M., Watson, J. V., Chen, H. Y. H., Lei, X., Schelhaas, M.-J., Lu, H., Gianelle, D., Parfenova, E. I., Salas, C., Lee, E., Lee, B., Kim, H. S., Bruelheide, H., Coomes, D. A., Piotto, D., Sunderland, T., Schmid, B., Gourlet-Fleury, S., Sonké, B., Tavani, R., Zhu, J., Brandl, S., Baraloto, C., Frizzera, L., Ba, R., Oleksyn, J., Peri, P. L., Gonmadje, C., Marthy, W., Brien, T. O., Martin, E. H., Marshall, A. R., Rovero, F., Bitariho, R., Niklaus, P. A., Alvarez-Loayza, P., Chamuya, N., Valencia, R., Mortier, F., Wortel, V., Engone-Obiang, N. L., Ferreira, L. V., Odeke, D. E., Vasquez, R. M., Lewis, S. L., and Reich, P. B.: Positive biodiversity-productivity relationship predominant in global forests, Science, 354, aaf8957,, 2016. 

Lohbeck, M., Poorter, L., Martínez-Ramos, M., and Bongers, F.: Biomass is the main driver of changes in ecosystem process rates during tropical forest succession, Ecology, 96, 1242–1252,, 2015. 

Luyssaert, S., Inglima, I., Jung, M., Richardson, A. D., Reichstein, M., Papale, D., Piao, S. L., Schulze, E.-D., Wingate, L., Matteuci, G., Aragao, L., Aubinet, M., Beer, C., Bernhofer, C., Black, K. G., Bonal, D., Bonnefond, J.-M., Chambers, J., Ciais, P., Cook, B., Davis, K. J., Dolman, A. J., Gielen, B., Goulden, M., Grace, J., Granier, A., Grelle, A., Griffis, T., Grünwald, T., Guidolotti, G., Hanson, P. J., Harding, R., Hollinger, D. Y., Hutyra, L. R., Kolari, P., Kruijt, B., Kuzsch, W., Lagergren, F., Laurila, T., Law, B. E., Le Maire, G., Lindroth, A., Loustau, D., Malhi, Y., Mateus, J., Migliavacca, M., Misson, L., Montagnani, L., Moncrieff, J., Moors, E., Munger, J. W., Nikinmaa, E., Ollinger, S. V., Pita, G., Rebmann, C., Roupsard, O., Saigusa, N., Sanz, M. J., Seufert, G., Sierra, C., Smith, M.-L., Tang, J., Valentini, R., Vesala, T., and Janssens, I. A.: CO2 balance of boreal, temperate, and tropical forests derived from a global database, Glob. Change Biol, 13, 2509–2537,, 2007. 

Malhi, Y.: The productivity, metabolism and carbon cycle of tropical forest vegetation, J. Ecol., 100, 65–75,, 2012. 

Malhi, Y., Baker, T. R., Phillips, O. L., Almeida, S., Alvarez, E., Arroyo, L., Chave, J., Czimczik, C. I., Fiore, A. D., Higuchi, N., Killeen, T. J., Laurance, S. G., Laurance, W. F., Lewis, S. L., Montoya, L. M. M., Monteagudo, A., Neill, D. A., Vargas, P. N., Patiño, S., Pitman, N. C., Quesada, C. A., Salomão, R., Silva, J. N. M., Lezama, A. T., Martínez, R. V., Terborgh, J., Vinceti, B., and Lloyd, J.: The above-ground coarse wood productivity of 104 Neotropical forest plots, Glob. Change Biol., 10, 563–591,, 2004. 

Malhi, Y., Girardin, C. A. J., Goldsmith, G. R., Doughty, C. E., Salinas, N., Metcalfe, D. B., Huaraca Huasco, W., Silva-Espejo, J. E., del Aguilla-Pasquell, J., Farfán Amézquita, F., Aragão, L. E. O. C., Guerrieri, R., Ishida, F. Y., Bahar, N. H. A., Farfan-Rios, W., Phillips, O. L., Meir, P., and Silman, M.: The variation of productivity and its allocation along a tropical elevation gradient: a whole carbon budget perspective, New Phytol., 214, 1019–1032,, 2017. 

McIntire, E. J. B. and Fajardo, A.: Facilitation as a ubiquitous driver of biodiversity, New Phytol., 201, 403–416,, 2014. 

Mori, A. S., Lertzman, K. P., and Gustafsson, L.: Biodiversity and ecosystem services in forest ecosystems: a research agenda for applied forest ecology, J. Appl. Ecol., 54, 12–27,, 2017. 

Moser, G., Röderstein, M., Soethe, N., Hertel, D., and Leuschner, C.: Altitudinal changes in stand structure and biomass allocation of tropical mountain forests in relation to microclimate and soil chemistry, in: Gradients in a Tropical Mountain Ecosystem of Ecuador, edited by: Beck, E., Bendix, J., Kottke, I., Makeschin, F., and Mosandl, R., Springer, Berlin and Heidelberg, Germany, 229–242,, 2008. 

Moser, G., Leuschner, C., Hertel, D., Graefe, S., Soethe, N., and Iost, S.: Elevation effects on the carbon budget of tropical mountain forests (S Ecuador): the role of the belowground compartment, Glob. Change Biol., 17, 2211–2226,, 2011. 

Muller-Landau, H. C.: Interspecific and inter-site variation in wood specific gravity of tropical trees, Biotropica, 36, 20–32,, 2004. 

Pan, Y., Birdsey, R. A., Phillips, O. L., and Jackson, R. B.: The structure, distribution, and biomass of the world's forests, Annu. Rev. Ecol. Evol. S., 44, 593–622,, 2013. 

Paoli, G. D. and Curran, L. M.: Soil nutrients limit fine litter production and tree growth in mature lowland forest of southwestern Borneo, Ecosystems, 10, 503–518,, 2007. 

Paquette, A. and Messier, C.: The effect of biodiversity on tree productivity: from temperate to boreal forests, Global Ecol. Biogeogr., 20, 170–180,, 2011. 

Pastor, J., Aber, J. D., McClaugherty, C. A., and Melillo, J. M.: Aboveground Production and N and P Cycling Along a Nitrogen Mineralization Gradient on Blackhawk Island, Wisconsin, Ecology, 65, 256–268,, 1984. 

Pierick, K., Leuschner, C., and Homeier, J.: Topography as a factor driving small-scale variation in tree fine root traits and root functional diversity in a species-rich tropical montane forest, New Phytol.,, 2021. 

Poorter, L., Wright, S. J., Paz, H., Ackerly, D. D., Condit, R., Ibarra-Manríquez, G., Harms, K. E., Licona, J. C., Martínez-Ramos, M., Mazer, S. J., Muller-Landau, H. C., Peña-Claros, M., Webb, C. O., and Wright, I. J.: Are functional traits good predictors of demographic rates? Evidence from five neotropical forests, Ecology, 89, 1908–1920,, 2008. 

Quesada, C. A., Phillips, O. L., Schwarz, M., Czimczik, C. I., Baker, T. R., Patiño, S., Fyllas, N. M., Hodnett, M. G., Herrera, R., Almeida, S., Alvarez Dávila, E., Arneth, A., Arroyo, L., Chao, K. J., Dezzeo, N., Erwin, T., di Fiore, A., Higuchi, N., Honorio Coronado, E., Jimenez, E. M., Killeen, T., Lezama, A. T., Lloyd, G., López-González, G., Luizão, F. J., Malhi, Y., Monteagudo, A., Neill, D. A., Núñez Vargas, P., Paiva, R., Peacock, J., Peñuela, M. C., Peña Cruz, A., Pitman, N., Priante Filho, N., Prieto, A., Ramírez, H., Rudas, A., Salomão, R., Santos, A. J. B., Schmerler, J., Silva, N., Silveira, M., Vásquez, R., Vieira, I., Terborgh, J., and Lloyd, J.: Basin-wide variations in Amazon forest structure and function are mediated by both soils and climate, Biogeosciences, 9, 2203–2246,, 2012. 

Rahbek, C., Borregaard, M. K., Colwell, R. K., Dalsgaard, B., Holt, B. G., Morueta-Holme, N., Nogues-Bravo, D., Whittaker, R. J., and Fjeldså, J.: Humboldt's enigma: What causes global patterns of mountain biodiversity?, Science, 365, 1108–1113,, 2019. 

Reich, P. B.: Key canopy traits drive forest productivity, Proc. R. Soc. Ser. B, 279, 2128–2134,, 2012. 

Reich, P. B. and Bolstad, P.: Productivity of evergreen and deciduous temperate forests, in: Terrestrial global productivity, edited by: Roy, J., Saugier, B., and Mooney, H. A., Academic Press, San Diego, USA, 245–283, 2001. 

Salazar, L., Homeier, J., Kessler, M., Abrahamczyk, S., Lehnert, M., Krömer, T., and Kluge, J.: Diversity patterns of ferns along elevational gradients in Andean tropical forests, Plant Ecol. Divers., 8, 13–24,, 2015. 

Sapijanskas, J., Paquette, A., Potvin, C., Kunert, N., and Loreau, M.: Tropical tree diversity enhances light capture through crown plasticity and spatial and temporal niche differences, Ecology, 95, 2479–2492,, 2014. 

Schnabel, F., Schwarz, J. A., Dănescu, A., Fichtner, A., Nock, C. A., Bauhus, J., and Potvin, C.: Drivers of productivity and its temporal stability in a tropical tree diversity experiment, Glob. Change Biol., 25, 4257–4272,, 2019. 

Schuur, E. A.: Productivity and global climate revisited: the sensitivity of tropical forest growth to precipitation, Ecology, 84, 1165–1170,[1165:PAGCRT]2.0.CO;2, 2003. 

Soethe, N., Lehmann, J., and Engels, C.: The vertical pattern of rooting and nutrient uptake at different altitudes of a south Ecuadorian montane forest, Plant Soil, 286, 287–299,, 2006. 

Spracklen, D. V. and Righelato, R.: Tropical montane forests are a larger than expected global carbon store, Biogeosciences, 11, 2741–2754,, 2014. 

Sullivan, M. J. P., Talbot, J., Lewis, S. L., Phillips, O. L., Qie, L., Begne, S. K., Chave, J., Cuni-Sanchez, A., Hubau, W., Lopez-Gonzalez, G., Miles, L., Monteagudo-Mendoza, A., Sonké, B., Sunderland, T., ter Steege, H., White, L. J. T., Affum-Baffoe, K., Aiba, S., de Almeida, E. C., de Oliveira, E. A., Alvarez-Loayza, P., Dávila, E. Á., Andrade, A., Aragão, L. E. O. C., Ashton, P. C. G. A. A., Baker, T. R., Balinga, M., Banin, L. F., Baraloto, C., Bastin, J.-F., Berry, N., Bogaert, J., Bonal, D., Bongers, F., Brienen, R., Camargo, J. L. C., Cerón, C., Moscoso, V. C., Chezeaux, E., Clark, C. J., Pacheco, Á. C., Comiskey, J. A., Valverde, F. C., Coronado, E. N. H., Dargie, G., Davies, S. J., De Canniere, C. K. M. N. D., Doucet, J.-L., Erwin, T. L., Espejo, J. S., Ewango, C. E. N., Fauset, S., Feldpausch, T. R., Herrera, R., Gilpin, M., Gloor, E., Hall, J. S., Harris, D. J., Hart, T. B., Kartawinata, K., Kho, L. K., Kitayama, K., Laurance, S. G. W., Laurance, W. F., Leal, M. E., Lovejoy, T., Lovett, J. C., Lukasu, F. M., Makana, J.-R., Malhi, Y., Maracahipes, L., Marimon, B. S., Junior, B. H. M., Marshall, A. R., Morandi, P. S., Mukendi, J. T., Mukinzi, J., Nilus, R., Vargas, P. N., Camacho, N. C. P., Pardo, G., Peña-Claros, M., Pétronelli, P., Pickavance, G. C., Poulsen, A. D., Poulsen, J. R., Primack, R. B., Priyadi, H., Quesada, C. A., Reitsma, J., Réjou-Méchain, M., Restrepo, Z., Rutishauser, E., Salim, K. A., Salomão, R. P., Samsoedin, I., Sheil, D., Sierra, R., Silveira, M., Slik, J. W. F., Steel, L., Taedoumg, H., Tan, S., Terborgh, J. W., Thomas, S. C., Toledo, M., Umunay, P., Valenzuela Gamarra, L., Vieira, I. C. G., Vos, V. A., Wang, O., Willcock, S., and Zemagho, L.: Diversity and carbon storage across the tropical forest biome, Sci. Rep.-UK, 7, 39102,, 2017. 

Tanner, E. V. J., Vitousek, P. M., and Cuevas, E.: Experimental investigation of nutrient limitation of forest growth on wet tropical mountains, Ecology, 79, 10–23[0010:EIONLO]2.0.CO;2, 1998. 

Tilman, D., Lehman, C. L., and Thomson, K. T.: Plant diversity and ecosystem productivity: theoretical considerations, P. Natl. Acad. Sci. USA, 94, 1857–1861,, 1997. 

Trabucco, A. and Zomer, R. J.: Global Aridity Index (Global-Aridity) and Global Potential Evapo-Transpiration (Global-PET) Geospatial Database, CGIAR Consortium for Spatial Information, CGIAR-CSI GeoPortal, available at: (last access: 31 August 2018), 2009. 

Tuck, S. L., O'Brien, M. J. O., Philipson, C. D., Saner, P., Tanadini, M., Dzulkifli, D., Godfray, H. C. J., Godoong, E., Nilus, R., Ong, R. C., Schmid, B., Sinun, W., Snaddon, J. L., Snoep, M., Tangki, H., Tay, J., Ulok, P., Wai, Y. S., Weilenmann, M., Reynolds, G., and Hector, A.: The value of biodiversity for the functioning of tropical forests: insurance effects during the first decade of the Sabah biodiversity experiment, Proc. R. Soc. B, 283, 20161451,, 2016. 

Unger, M., Leuschner, C., and Homeier, J.: Variability of indices of macronutrient availability in soils at different spatial scales along an elevation transect in tropical moist forests (NE Ecuador), Plant Soil, 336, 443–458,, 2010. 

Unger, M., Homeier, J., and Leuschner, C.: Effects of soil chemistry on tropical forest biomass and productivity at different elevations in the equatorial Andes, Oecologia, 170, 263–274,, 2012. 

Unger, M., Homeier, J., and Leuschner, C.: Relationships among leaf area index, below-canopy light availability and tree diversity along a transect from tropical lowland to montane forests in NE Ecuador, Trop. Ecol., 54, 33–45, 2013. 

van der Sande, M. T., Arets, E. J., Peña-Claros, M., Hoosbeek, M. R., Cáceres-Siani, Y., van der Hout, P., and Poorter, L.: Soil fertility and species traits, but not diversity, drive productivity and biomass stocks in a Guyanese tropical rainforest, Funct. Ecol., 32, 461–474,, 2018. 

van de Weg, M. J., Meier, P., Williams, M., Girardin, C., Malhi, Y., Silva-Espejo, J., and Grace, J.: Gross primary productivity of a high elevation tropical montane cloud forest, Ecosystems, 17, 751–764,, 2014. 

Vicca, S., Luyssaert, S., Peñuelas, J., Campioli, M., Chapin III, F. S., Ciais, P., Heinemeyer, A., Högberg, P., Kutsch, W. L., Law, B. E., Malhi, Y., Papale, D., Piao, S. L., Reichstein, M., Schulze, E. D., and Janssens, I. A.: Fertile forests produce biomass more efficiently, Ecol. Lett., 15, 520–526,, 2012.  

Vitousek, P. M. and Sanford, R. L.: Nutrient cycling in moist tropi-cal forest, Annu. Rev. Ecol. Syst., 17, 137–167, 1986. 

Wallis, C. I., Homeier, J., Peña, J., Brandl, R., Farwig, N., and Bendix, J.: Modeling tropical montane forest biomass, productivity and canopy traits with multispectral remote sensing data, Remote Sens. Environ., 225, 77–92,, 2019. 

Werner, F. A. and Homeier, J.: Is tropical montane forest heterogeneity promoted by a resource-driven feedback cycle? Evidence from nutrient relations, herbivory and litter decomposition along a topographical gradient, Funct. Ecol., 29, 430–440,, 2015. 

Wittich, B., Horna, V., Homeier, J., and Leuschner, C.: Altitudinal change in the photosynthetic capacity of tropical trees: a case study from Ecuador and a pantropical literature analysis, Ecosystems, 15, 958–973,, 2012. 

Wolf, K., Veldkamp, E., Homeier, J., and Martinson, G. O.: Nitrogen availability links forest productivity, soil nitrous oxide and nitric oxide fluxes of a tropical montane forest in southern Ecuador, Global Biogeochem. Cy., 25, GB4009,, 2011. 

Wright, S. J., Yavitt, J. B., Wurzburger, N., Turner, B. L., Tanner, E. V. J., Sayer, E. J., Santiago, L. S., Kaspari, M., Hedin, L. O., Harms, K. E., Garcia, M. N., and Corre, M. D.:Potassium, phosphorus, or nitrogen limit root allocation, tree growth, or litter production in a lowland tropical forest, Ecology, 92, 1616–1625,, 2011. 

Zimmermann, M., Meir, P., Bird, M. I., Malhi, Y., and Ccahuana, A. J. Q.: Temporal variation and climate dependence of soil respiration and its components along a 3000 m altitudinal tropical forest gradient, Global Biogeochem. Cy., 24, GB4012,, 2010. 

Short summary
We studied aboveground productivity in humid tropical montane old-growth forests in two highly diverse Andean regions with large geological and topographic heterogeneity and related productivity to tree diversity and climatic, edaphic and stand structural factors. From our results we conclude that the productivity of highly diverse Neotropical montane forests is primarily controlled by thermal and edaphic factors and stand structural properties, while tree diversity is of minor importance.
Final-revised paper