Dennis Castillo-Figueroa
Andrés González-Melo
Juan M. Posada- Biology Department, Faculty of Natural Sciences, Universidad del Rosario, Bogota, Colombia
Wood density (WD) is a key functional trait related to ecological strategies and ecosystem carbon dynamics. Despite its importance, there is a considerable lack of information on WD in tropical Andean forests, particularly regarding its relationship with forest succession and ecosystem carbon cycling. Here, we quantified WD in 86 upper Andean tree and shrub species in central Colombia, with the aim of determining how WD changes with forest succession and how it is related to productivity. We hypothesized that WD will increase with succession because early successional forests will be colonized by acquisitive species, which typically have low WD, while the shaded understory of older forests should favor higher WD. We measured WD in 481 individuals from 27 shrub and 59 tree species, and quantified aboveground biomass (AGB), canopy height, net primary production (NPP) and species composition and abundance in 14, 400-m2, permanent plots. Mean WD was 0.513 ± 0.114 (g/cm3), with a range between 0.068 and 0.718 (g/cm3). Shrubs had, on average, higher WD (0.552 ± 0.095 g/cm3) than trees (0.488 ± 0.104 g/cm3). Community weighted mean WD (CWMwd) decreased with succession (measured as mean canopy height, AGB, and basal area); CWMwd also decreased with aboveground NPP and stem growth. In contrast, the percentage of NPP attributed to litter and the percent of shrubs in plots increased with CWMwd. Thus, our hypothesis was not supported because early successional forests had higher CWMwd than late successional forests. This was related to a high proportion of shrubs (with high WD) early in succession, which could be a consequence of: 1) a low seed availability of trees due to intense land use in the landscape and/or 2) harsh abiotic conditions early in succession that filter out trees. Forest with high CWMwd had a high %NPP attributed to litter because they were dominated by shrubs, which gain little biomass in their trunks. Our findings highlight the links between WD, succession and carbon cycling (biomass and productivity) in this biodiversity hotspot. Thus, WD is an important trait that can be used to understand upper Andean forest recovery and improve forest restoration and management practices.
Introduction
Wood density (WD) is a key functional trait for understanding tree functioning and community ecology (Chave et al., 2005; Chave et al., 2009; Oliveira et al., 2019; González-Melo, 2021). It is also important in ecosystem studies (Padilha and Marco Júnior, 2018), that require estimates of aboveground biomass, productivity, and carbon emissions (Fearnside, 1997; Muller-Landau, 2004; Chave et al., 2009). In an eco-physiological context, high WD has been related to lower water transport capacity in trunks and branches (Stratton et al., 2000; Scholz et al., 2007; Meinzer et al., 2008), but higher drought tolerance (Markesteijn et al., 2011). WD has also an inverse relationship with growth (Muller-Landau, 2004; Osunkoya et al., 2007; Kraft et al., 2010), mortality rates (Chave et al., 2009), and recruitment rates (Suzuki, 1999; Wright et al., 2003; King et al., 2006). Likewise, some studies have indicated that WD plays an important role in ecosystem processes as it is related to aboveground biomass (AGB), the rate of biomass accumulation (Baker et al., 2004; Chave et al., 2005; Stegen et al., 2009; Oliveira et al., 2019) and wood decomposition (Aryal et al., 2022). AGB is also related to forest successional status (i.e., age; Brown and Lugo, 1990; Garnier et al., 2004; Garnier et al., 2016; Chazdon et al., 2016; Poorter et al., 2021a) and successional changes in WD vary between ecosystems (Poorter et al., 2021b). For instance, WD increases with forest age in lowland tropical rain forests due to a shade tolerance filter later in succession (Poorter et al., 2021b). Another aspect to consider when studying WD is life form. In particular, WD can differ considerably between shrubs and trees (Martínez-Cabrera et al., 2009) because these functional groups have different height, and can vary in their resource needs (i.e., water, light or nutrients; Martínez-Cabrera et al., 2009).
Several studies related to WD have made significant contributions to both regional (Chave et al., 2006; Oliveira et al., 2019; Báez et al., 2022) and global trait databases (Chave et al., 2009; Zanne et al., 2009; Henry et al., 2013; Kattge et al., 2020). In the tropics, the bulk of research on WD has been carried out in lowland ecosystems like the Amazonian rainforest (Nogueira et al., 2008), the Brazilian Cerrado (Roitman et al., 2018), Panamanian lowland forest (Muller-Landau, 2004; Wright et al., 2010) and the Atlantic forest (Padilha and Marco Júnior, 2018), among others. Nevertheless, there is a considerable gap in our knowledge of WD variations in tropical montane forests, particularly in geographical regions above 2500 meters elevation (Álvarez-Dávila et al., 2017).
Tropical mountain forests play pivotal roles in biodiversity conservation, and water and carbon cycles (Galbraith et al., 2014; Duque et al., 2021). In particular, the upper Andes Mountains of Northern South America (2600-3200m) represent a global hotspot of biological diversity due to their significant of endemism and significant of transformation (Hernández-Camacho et al., 1992; Myers et al., 2000; Myster, 2021). However, Upper Andean Tropical Forests (UATF; above ca. 2500 m elevation) are among the least known ecosystems (Etter and van Wyngaarden, 2000), especially in relation to species functional traits, ecological succession, and ecosystem processes (Castillo-Figueroa, 2021). These forests have experienced extensive habitat loss due to large-scale deforestation and conversion of lands to agriculture and urban areas (Etter and van Wyngaarden, 2000; Etter and Villa, 2000; Etter et al., 2008; Etter et al., 2021), that has resulted in a mosaic of forests with different successional status, ranging from early-secondary to small and isolated old-growth forest patches (Rubiano et al., 2017; Calbi et al., 2020; Calbi et al., 2021; Hurtado-M et al., 2021).
Trait-based ecology predicts that WD should increase during succession, because in the first successional stages forests are commonly colonized by species with a high resource acquisition strategy, with rapid growth rates but low WD, which allow them to take advantage of high light and nutrient availability (Chave et al., 2009; Chazdon, 2014; Garnier et al., 2016; but see Poorter et al., 2021b). In contrast, old-growth forests are usually colonized by species with conservative strategies, with low growth rate and high WD that allow them to tolerate understory conditions (Chao et al., 2008; Chave et al., 2009; Poorter et al., 2021b). Yet, evidence supporting this prediction comes mainly from tropical lowland rainforests, and it is still unclear whether this pattern holds true for UATF. In this sense, the relationship of WD with forest successional status (measured by structural attributes such as AGB, basal area, and canopy height) would represent a step forward in our understanding of UATF function and regeneration.
In this study, we measured WD for 86 species along a successional gradient in UATF. We compared WD between life forms and assessed how forest structural attributes and net primary production (NPP) are related to community-weighted mean wood density (CWMwd) along a successional gradient. We addressed the following questions: (1) What is the range of variation of WD in species from UATF in comparison to regional and global databases? (2) Does WD differ between shrubs and trees? (3) How is CWMwd related to forest successional status (using basal area, AGB, and canopy height as proxies of succession) and to NPP along a successional gradient? Given that the studied UATF are exposed to mesic conditions, our main hypothesis was that CWMwd will increase with successional status because high resource availability early in succession will favor acquisitive species, that will be gradually replaced by species with more conservative traits.
Materials and methods
Study area
This study was conducted in the Eastern Colombian Andes (Figure 1). The Andean region of Colombia covers 24.5% of the country (Etter and van Wyngaarden, 2000), it is the most populated region, with almost the 70% of the population, and harbors the economic core of the country (Departamento Administrativo Nacional de Estadística (DANE), 2018). As such, the Andes of Colombia have been heavily transformed by human activities. The study was carried out in the Cundiboyacense Highplain, which includes the Capital city of Bogotá and its surrounding municipalities. In these peri-urban sites, cattle raising, agriculture and mining are the dominant economic activities (Montañez et al., 1994; Antonio-Fragala and Obregón-Neira, 2011). Average annual atmospheric temperature in the area is 14 °C, and mean annual precipitation varies from 600 mm in the center of the mountain range to 1200 mm in the western part (Clerici et al., 2016). Precipitation typically displays a bimodal pattern with two rainy seasons: one occurring from April to June, and another between September and November (Mendoza and Etter, 2002).
Figure 1 Study area. Triangles show the locations of the upper Andean tropical forests plots where sampling took place. Llanitos is a site within the municipality of Guatavita, where some species were sampled in 50x50 m plots.
We established 14 20x20 m plots and two 50x50 m plots in secondary and mature Andean Mountain forests at elevations ranging from 2685 to 3140 m. These plots are part of a larger network consisting of 36 20x20 m plots and eight 50x50 m plots. Forests were classified as secondary or mature by local experts, based on their structural attributes (basal area, tree height, tree density) and species composition. We mostly focused on the 14 20x20 m plots (see below), but we also sampled WD in some species in the 50x50m plots (Llanitos). The plots are on private properties/reserves, distributed in three municipalities (Tabio, Guasca and Guatavita) and one site (Torca) located close to the boundaries of the Capital District of Bogotá (Figure 1). In the 20x20 m plots, we tagged and identified all individuals with a diameter at ‘ankle’ height (DAH, 5 above the soil surface; Hurtado-M et al., 2021) larger or equal to 5 cm. We used this approach to include shrubs that usually branched near the ground or to include individuals with twisted stems. We also recorded diameter at breast height (DBH) to estimate biomass using allometric equations (see below). In the 50x50 m plots, we only sampled individuals with a DBH larger or equal to 10 cm. In the 20x20 m plots, the five dominant plant families were: Ericaceae, Melastomataceae, Cunoniaceae, Primulaceae, and Asteraceae, representing 57.7% of individuals. The five dominant species were: Weinmannia tomentosa Linnaeus filius 1782, Cavendishia bracteata Hoerold 1909, Miconia squamulosa Triana 1872, Myrcianthes leucoxyla McVaugh 1963, and Myrsine latifolia Spreng. 1824, with 46.16% of individuals. In total, 79 shrub and tree species were identified in the permanent plots. We also included seven common species found around the plots to complete a dataset of 86 plant species. Plant nomenclature was based on the Tropicos® Database from the Missouri Botanical Garden (https://www.tropicos.org/) and the World Flora Online Database (http://www.worldfloraonline.org).
Forest structural attributes and successional gradient
In each of the 14 20x20 m plots (Table 1) we quantified the percentage of shrubs (% shrubs), stem basal area (m2) based on DAH, canopy height (m) and AGB (Mg C ha-1); we assumed that average C content was 0.5 (Petersson et al., 2012). Species were classified as shrubs if they had permanent wood, ramified from their basal part or below the soil surface, had no distinctive trunk, and their height was lower than 5 m (Orsham, 1989; Rundel, 1991; Mahecha et al., 2012). To quantify AGB we measured DBH for all individuals with a DAH higher or equal to 5 cm. AGB per plot was estimated as the average of two allometric equations that are specifically formulated for Andean montane forests by (1) Sierra et al. (2007) and (2) Pérez and Díaz (2010), respectively:
Table 1 Characteristics of the 14 permanent plots in UATF, in Colombia.
Considering that AGB is a good indicator of forest recovery (Brown and Lugo, 1990; Garnier et al., 2004; Lohbeck et al., 2015a; Garnier et al., 2016), we used it as a proxy of ecological succession altogether with canopy height (h, m), and basal area calculated from DAH (m2). Canopy height was calculated as the weighted average of the upper quartile of all tree height measurements.
To measure net primary production (NPP; Mg C ha-1 y-1), we collected litterfall in 0.25 m2 litter traps in the 14 20x20 m plots (10 per plot, 140 in total). The trap frames were made of ½ inch PVC tubes and covered with a 1 mm nylon mesh screen; traps were elevated ca. 1 m above ground. Litter was collected between May 2016 and April 2017 every 15-20 days to minimize loss by decomposition. Litter samples were oven-dried at 60 °C until constant weight and then weighed with an analytical scale. Branches with a diameter above 1 cm were excluded, and branches with a diameter lower than 1 cm, but wider than the trap were cut to the length of the traps. We also measured stem growth by taking DBH measures between 2013 and 2016. NPP was calculated as the sum of annual litter production and stem biomass production; dead trees and new recruits were treated following the methodology of Clark et al. (2001).
Tree or shrub heights were measured with a clinometer or a telescopic measuring rod. Hemispherical photographs were taken using a Canon Inc. EW-77 digital camera with a 180° Fisheye Zoom Lens (EF 8-15mm 1:4, Japan). Nine hemispherical photographs were taken per plot and processed with Gap Light Analyzer 2.0 (Frazer et al., 1999) to obtain LAI and canopy openness. Plot LAI and canopy openness are reported as average values per plot (Table 1).
WD measurements
In 2019, we measured WD for all species found in the permanent plots, plus seven woody species sampled in the surrounding areas (Figure 1) following standardized methods (Pérez-Harguindeguy et al., 2013). Wood samples were collected from stems using a 5-mm increment borer. Samples were sapwood taken at 1.3 m for both trees and shrubs. Some shrubs had no defined stem at 1.3 m and sampling took place below that height, but never lower than 1.1 m. An average of five individuals were sampled per species. Yet, for dominant species we collected up to 12 individuals, while for some rare species we sampled between one and three individuals. For each wood sample, we measured green volume with the water-displacement method, and then oven-dried it at 105 °C, during 72 h, to obtain wood dry mass (Chave et al., 2005). WD (g/cm-3) was calculated as the ratio between dry mass and green volume (Williamson and Wiemann, 2010; Oliveira et al., 2019).
Data analyses
To characterize interspecific variation in WD for our study area, we plotted a kernel density function of our data, and compared it to global (Zanne et al., 2009) and regional density functions, including one in the Neotropics (Chave et al., 2006) and one in a lowland tropical forest of French Guiana from the BRIDGE project (Baraloto et al., 2010). We also explored variation in WD between families via one-way Kruskal-Wallis test. To compare WD between life forms (trees vs shrubs), we performed a nonparametric Mann Whitney test for pairwise groups. This test was also used for pairwise comparisons of WD between UATF and each of the three databases.
To examine whether forest attributes are related to WD at the community level, we calculated CWMwd for each of the 14 plots using species mean WD values weighted by species biomass following Garnier et al. (2004). To compare CWMwd to forest structural and productivity variables we used reduced major axis regression (RMA) that models residual variation of both axes (Pearson, 1901; Warton et al., 2012). We regressed CWMwd against AGB, basal area, canopy height, aboveground NPP, productivity of stems (i.e., annual biomass gain of stems), litterfall production and % shrub individuals in plots. All variables were log10 transformed before analyses. Statistical analyses were performed in Rwizard 4.3 (Guisande et al., 2014) and PAST 4.1 (Hammer et al., 2001).
Results
WD variation in UATF
In total, we measured WD in 481 individuals that belonged to 86 species, 59 genera and 39 families (Table 2). Species mean WD was 0.513 ± 0.114 g/cm3 and ranged from 0.068 (Vasconcellea pubescens with individuals as low as 0.065 g/cm3) to 0.718 (Ulex europeus with individuals up to 0.879 g/cm3) (Figure 2A). WD varied between families (H=220.75; df=38; p<0.0001), being Caricaceae, Betulaceae, Winteraceae, Clethraceae and Araliaceae the families with the lowest mean WD values (0.068-0.399 g/cm3), while Gentianaceae, Clusiaceae, Celastraceae, Rosaceae, and Fabaceae showed the highest values (0.619-0.718 g/cm3) (Figure 2B). The distribution of WD from UATF had a lower mean, a narrower range, and a smaller standard deviation than the global and Neotropical WD distributions (Figure 3). UATF also had a lower mean that a lowland tropical forest in French Guiana, but similar standard deviation and range. The overlap with the other datasets was: 66.59% (global), 59.18% (Neotropical) and 50.83% (lowland tropical forests from French Guiana). Significant differences were found in the pairwise comparisons between WD of UATF and the global (U= 222593; p<0.0001), Neotropical (U=55450; p<0.0001) and lowland tropical forests from French Guiana (U=6604; p<0.0001).
Table 2 List of species and their corresponding wood density in upper Andean tropical forests.
Figure 2 (A) Frequency histograms of WD by species in upper Andean tropical forests (B) Box plots of WD by families; the horizontal line is the median, the box is 50% of data, while the whisker extends to the furthest data point in each wing that is within 1.5 times the Interquartile range. Family acronyms are found in Table 2.
Figure 3 Density plots comparing wood density in upper Andean tropical forests, a lowland tropical rainforest in French Guiana (Baraloto et al., 2010), the Neotropics (Chave et al., 2006), and a global database (Zanne et al., 2009). Mean, median, standard deviation and range are shown on the upper right side.
Life form and WD variation
Shrubs comprised 27 species whereas there were 59 species of trees (Table 2). Trees had an average WD of 0.488 ± 0.104 g/cm3 whilst for shrubs it was 0.552 ± 0.095 g/cm3. Shrubs had significantly higher mean values of WD than trees (W=1095; p=0.005; n=86; Figure 4A), although the two distributions overlapped considerably. In addition, the proportion of shrubs per plot was positively related to CWMwd (R2 = 0.33; p=0.03; CI =[0.127, 0.273]; n=14; Figure 4B).
Figure 4 (A) Beanplots comparing shrubs (27 species) and trees (59 species) from upper Andean tropical forests plots. Each horizontal line represents the mean value of a species. (B) Reduced major axis regression between the % of individuals that are shrubs in each individual plot and community weighted mean (CWM) of wood density per plot. Values on the right corner of the plot are the coefficient of determination (R2) and slope (m). Grey circles correspond to secondary forests and black circles to mature forests plots.
WD, AGB and community attributes
AGB biomass ranged from 12.1 Mg C ha-1 to 110 Mg C ha-1, canopy openness between 5.85 to 17.89, and LAI between 2.52 and 4.38 (Table 1). On average, old-growth forests had an AGB of 83.16 Mg C ha-1 ( ± 18.87), a canopy openness of 7.83 (± 2.52), and a LAI of 3.32 ( ± 0.47), while secondary forests presented an AGB of 28.37 Mg C ha-1 ( ± 11.63), a canopy openness of 10.15 (± 4.53), and a LAI of 3.18 ( ± 0.67). AGB was significantly higher in mature forests than secondary ones (U= 28; p=0.017), but canopy openness (U=22; p= 0.17) and LAI (U=15; p=0.86) did not show any differences between successional stages.
We found significant negative RMA relationships between CWMwd and three variables related to forest successional status: basal area (R2=0.52; p=0.003; CI=[-0.232,-0.063]; n=14; Figure 5A), AGB (R2=0.54; p=0.002; CI=[-0.168,-0.075]; n=14; Figure 5B), and canopy height (R2=0.48; p=0.006; CI=[-0.449,-0.198]; n=14; Figure 5C). Regarding ecosystem productivity, we found an inverse relationship between CWMwd and aboveground NPP (R2=0.33; p=0.031; CI=[-6.508, -1.849]; n=14; Figure 5F) as well as CWMwd and stem growth (R2=0.34, p=0.029; CI=[-15.942,-2.770]; n=14; Figure 5D), but a positive one between CWMwd and %NPP attributed to litter (R2=0.35; p=0.024; CI=[0.973, 3.364]; n=14; Figure 5E).
Figure 5 Reduced major axis regressions of community weighted mean (CWM) of wood density vs. basal area (m2; A), aboveground biomass (AGB; Mg C ha-1; B), and forest canopy height (average of the upper quartile of tree height; C); each point correspond to a plot. The bottom panels represent CWM of wood density vs. aboveground net primary productivity (Total NPP; Mg C ha-1 y-1; D), stem biomass growth (NPP stems; Mg C ha-1 yr-1; E), and the % of aboveground NPP attributed to leaf litter (% NPP to litter; F). Values on the right corner of each graph are the coefficient of determination (R2) and slope (m). All variables are log10 transformed. Grey circles correspond to secondary forests and black circles to mature forests plots.
Discussion
WD in upper Andean tropical forests
Wood traits have been under-sampled above 2500 m in the Andean region, and previous studies have reported WD values for only a few species in UATF (e.g., Chave et al., 2006; Vasquez-Valderrama and Solorza-Bejarano, 2018; López-Camacho et al., 2020; González-Melo, 2021). Here, we presented a comprehensive dataset of WD for 53 additional plant species for this ecoregion. We found that species mean WD for UATF (0.513 ± 0.114 g/cm3) was significantly lower than that reported in a global WD database (0.616 g/cm3), in Neotropical forests (0.647 g/cm3), as well as for lowland tropical wet forests in French Guyana (0.647 g/cm3; Figure 3). The mean value that we reported falls in the light-weight wood category according to Riesco et al. (2019). Other studies have found that WD decreases with elevation and decreasing temperature (Williamson, 1984; Chave et al., 2009), which agrees with our results (Figure 3; Table 2). Lower WD at high altitude can have several non-exclusive explanations. First, in tropical mountain forests of the Ecuadorian Andes, Moser et al. (2011) reported a noticeable reduction in NPP above 2000 m and suggested that wood production can be environmentally constrained at these higher altitudes. At high elevation, soil nitrogen availability and atmospheric temperatures are lower, which can lead to a reduction in leaf area, and consequently a decrease in carbon gain (Moser et al., 2011; Homeier and Leuschner, 2020). Similarly, Leuschner et al (2007; 2013). found a remarkable increase in the allocation of NPP in roots at high elevations, and a lower investment in tree trunk biomass and AGB, suggesting nutrient limitation for tree growth in cold environments. Yet, estimations of forest carbon budget along an altitude gradient in Perú showed that the fraction of NPP allocated to stems and roots do not change with elevation, despite lower gross and net primary production at higher than at lower elevation (Malhi et al., 2017).
We suggest that, in addition, a lower WD in UATF can be related to a more open canopy. The understory of lowland tropical rain forests acts as an important filter that favor species with conservative ecological strategies with high WD (e.g., Kitajima, 1994; Alvarez-Clare and Kitajima, 2007; Poorter et al., 2021b). In contrast, vertical stratification in UATF is simplified because of lower height of emergent trees, lower tree stem diameter, and lower leaf area (Myster, 2021), which may increase the light transmittance to the understory. This agrees with our results as LAI in UATF was lower (mean=3.26, Table 1) compared to lowland forests (LAI=5.1; e.g., Moser et al., 2007), suggesting that selection for shade tolerance (and high WD) could be less pronounced in UATF. Moreover, lower canopy height as well as lower diameter of tree stems at higher elevation (Myster, 2021) likely reduces the needs to invest in higher mechanical strength (and hence higher WD; Chave et al., 2009). Lastly, although WD exhibits an important intra-specific variation in response to environmental conditions (Patiño et al., 2009), some studies have shown that this trait is highly influenced by phylogeny (Chave et al., 2006; Chave et al., 2009). Thus, low WD in Andean forests can also be influenced by environmental filters favoring certain taxonomic groups with lower WD.
Ulex europaeus was the only exotic species in our plots and had the highest WD of all the 86 species (Table 2). Note, however, that excluding it from our analyses did not alter the overall results, due to its low biomass in the plot where it was present (data not shown). U. europaeus is considered an aggressive invasive species worldwide (Lowe et al., 2000), including the Highlands of Colombia that have seen an important expansion in its range since the 1950’s (Vargas et al., 2009). Its invasive success has been related to high reproductive and survival rates, and high resource acquisition capacity (Beltrán and Cataño, 2014). This species comes from the Mediterranean coast and its high WD is likely related to its shrub life-form (see below) and its evolutionary origin from an arid biome. Yet, it will be interesting to determine if its high WD also contributes to its invasive potential in UATF.
Life forms
Shrubs showed higher WD than trees (Figure 4), as has been reported in previous studies (Jacobsen et al., 2007; Wheeler et al., 2007). This difference could be partially explained by habitat preferences and life-history variations between shrubs and trees. In general, shrubs are common in low-resource environments, such as dry ecosystems or shaded forest understories, characterized by light, water, and/or nutrient constraints (Petit and Hampe, 2006; Martínez-Cabrera et al., 2011). Moreover, trees and shrubs may differ in their wood structures because of differences in hydraulic demands. For example, since shrubs transport water to shorter distances compared to trees, allocation of stem cross-sectional area to vessels is expected to be lower in shrubs than trees. Shrubs may, therefore, allocate a higher percentage of stem volume to fibers (Martínez-Cabrera et al., 2009), via higher fiber wall to lumen ratios and/or higher fiber fractions (Ziemińska et al., 2015).
Relationship to successional gradient and forest structural attributes
Our hypothesis that WD will increase during succession was not supported, and our results were in contrast with the pattern found on lowland tropical rain forests (Chazdon, 2014; Lohbeck et al., 2015b; Poorter et al., 2019; Poorter et al., 2021b). WD is expected to reflect a trade-off between growth and survival (Poorter et al., 2008; Chave et al., 2009). That is, species that invest in dense wood have higher survival in the understory due to an increase protection against herbivory and mechanical damage, while species that produce lighter wood tend to grow faster, via higher hydraulic conductivity and low stem construction costs, which allow them to take advantage of higher resource availability early in succession (King et al., 2006; Chave et al., 2009; Chua and Potts, 2018). Here, we found that early successional forests (i.e., with low AGB, basal area and height) had higher CWMwd than late successional forests (i.e., with high AGB, basal area and height). Interestingly, this result is similar to the pattern found in tropical dry forests where species with drought tolerance traits (i.e., high WD) are dominant in early succession as a response of water limitations and are gradually replaced by species with more acquisitive strategies during forest recovery (Lohbeck et al., 2015b; Poorter et al., 2019). In our study, higher CWMwd is largely related to a high abundance of shrubs (with higher WD than trees, Figure 4B) early in succession (Figure 6).
Figure 6 Illustration of changes in aboveground biomass and the percentage of shrubs along the successional gradient. Reduced major axis regression between the percentage of shrubs and aboveground biomass (AGB) in plots (p=0.0007; R2=0.63; n=14). Grey circles correspond to secondary forests and black circles to mature forests plots.
The high proportion of shrubs during the first stages of forest recovery can have several explanations. First, seed availability of tree species can be limited due to intense and prolonged land use in the region, which likely resulted in a disproportionate cutting of trees over shrubs. Second, shrubs tend to produce more seeds per square meter of canopy per year than trees (Moles et al., 2004; Moles et al., 2005), which favor their colonization in the landscape, especially in human-modified ecosystems. Thirdly, more harsh abiotic conditions may have filter out trees in the early stages of succession. A study conducted in the same plots showed that tree seedling composition was highly sensitive to edaphic conditions including pH, available phosphorous and calcium, which can be altered by anthropogenic pressures (Hurtado-M et al., 2021). Considering adverse conditions for seedlings in early-successional mountain forests, such as nutrient-poor soils (Vitousek et al., 1988; Tanner et al., 1998), harsh winds (Sugden, 1986), periodic drought stress (Van Steenis, 1972), and frost, it is likely that environmental filtering at these first stages can favor the establishment of shrubs with denser woods and higher stress tolerance, instead of trees with lighter woods. In line with this, some studies have found that shrubs tend to present not only higher WD (Larjavaara and Muller-Landau, 2012; Santiago et al., 2018), but also thicker roots than trees (Fortunel et al., 2012; Valverde-Barrantes et al., 2017), which can favor their establishment under less favorable abiotic conditions of drought and low fertility that can be present early in succession in human-disturbed UATF.
We also found that the percentage of aboveground NPP attributed to litter increased with CWMwd (Figure 5F). This result can again be attributed to the dominance of shrubs in early successional forests (Figure 6). Shrubs gained little biomass in their trunks because they attain their maximum height after a few years, so most NPP was related to litter production and less to stem growth (Figures 5E, F). Conversely, in old-growth forests the presence of shrubs was lower, and many trees in the plots had not attained their maximum height, resulting in proportionally more stem growth, which is reflected in both the increase of total NPP and the fraction of NPP attributed to stem biomass gain (Figures 5D, E). These results are important in the context of carbon cycling and forest succession as some studies have found that WD is inversely related to carbon gains (Hubau et al., 2020). Thus, old-growth forests from UATF with lower CWMwd can produce more biomass and act as significant carbon sinks, as has been suggested recently (Duque et al., 2021).
The large gap of information about functional traits of species from UATF call the attention for future research efforts. WD is particularly relevant in trait-based ecology given its link to species ecological strategies, community organization and biogeochemical carbon cycling (Evans et al., 2000; Chave et al., 2005; Van Gelder et al., 2006; Chave et al., 2009). The study of local or global factors that influence WD could help better predict the biomass stored in Andean forests. This is especially important under a global change scenario, where the impact of human-dominated landscapes and climate change in Andean ecosystems keep growing, but is still poorly understood. In the Colombian Andes, less than 25% of the original forest coverage is still present (Armenteras et al., 2003; Etter et al., 2021). Therefore, to ensure adequate biomass recovery of secondary forests, it is imperative to fill knowledge gaps on the major factors that can influence carbon storage such as WD.
Conclusions
We measured WD in 86 species from UATF, which represent, for 53 of them, their first known report. These data could be incorporated into regional and global databases, thus contributing to help filling an important gap for these ecosystems (Báez et al., 2022). UATF had a lower average WD than the global, neotropical or lowland French Guyana forest. This low WD can be related to lower productivity, higher allocation to roots, lower tree stature, and/or a weaker “shade filter” due to higher light availability in the understory of UATF. We also found that disturbance, dispersal, and establishment filters appear to have favored shrubs with higher WD in early successional forests. In contrast, old-growth forests had a lower CWMwd due to a decrease in the abundance of shrubs. The decrease in aboveground NPP with WD can also be partially attributed to change in the abundance of shrubs, which showed reduced annual biomass gain in their trunks and a large fraction of NPP being related to litter production. Overall, we found that WD is a key functional trait related to plant community assembly and carbon cycling along successional gradients. Further research should determine how WD is related to other components of the carbon cycle and to different functional traits. That will deepen our understanding of the evolution of ecological strategies in UATF species and help us make better predictions about the benefit of the restoration of these forests on ecosystem services.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author.
Author contributions
DC-F: Writing – original draft, Writing – review & editing, Conceptualization, Formal Analysis, Investigation, Methodology, Data curation. AG-M: Conceptualization, Data curation, Investigation, Methodology, Writing – original draft. JP: Conceptualization, Investigation, Methodology, Writing – original draft, Formal Analysis, Funding acquisition, Project administration, Resources, Supervision, Validation, Visualization, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was part of the project “Estudio de dinámicas socio-ecológicas ante escenarios de cambio climático en bosques secundarios periurbanos Altoandinos” lead by JP and funded by MINCIENCIAS (No. FP44842-046-2017). Funding for the publication of this article came from the Dirección de Investigación e Innovación at University of el Rosario and from a Research Incentive Fund from University of el Rosario assigned to JP.
Acknowledgments
The permanent plots were originally set up by Natalia Norden, Ana Belén Hurtado, Carolina Álvarez and JP. We thank Carolina Álvarez, Sergio Venegas-Carrillo, María Alejandra Molina, Nicole Morales, and Luisa Vargas for their invaluable support in collecting and processing the data. We are grateful to the owners, residents and managers of the sites for their generosity in allowing us to conduct research in their properties; particularly, we thank Gonzalo Martínez, Martha Giraldo, Estefania Cabo, Gustavo Lobo-Guerrero, Carlos Castillo and the staff of Encenillo Biological Reserve (Fundación Natura) for facilitating access to the study sites. We are also grateful to Christopher Baraloto for sharing dataset of French Guiana from BRIDGE project and the Global Wood density database group. Finally, we want to thank to Laura Garzón-Salamanca for her help in the preparation of the map in Figure 1.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Alvarez-Clare, S., Kitajima, K. (2007). Physical defence traits enhance seedling survival of neotropical tree species. Funct. Ecol. 21 (6), 1044–1054. doi: 10.1111/j.1365-2435.2007.01320.x
Álvarez-Dávila, E., Cayuela, L., González-Caro, S., Aldana, A. M., Stevenson, P. R., Phillips, O., et al. (2017). Forest biomass density across large climate gradients in northern South America is related to water availability but not with temperature. PloS One 12 (3), p.e0171072. doi: 10.1371/journal.pone.0171072
Antonio-Fragala, F., Obregón-Neira, N. (2011). Recharge estimation in aquifers of the bogota savannah. Ing. Univ. 15, 145–169. doi: 10.11144/Javeriana.iyu15-1.erma
Armenteras, D., Gast, F., Villareal, H. (2003). Andean forest fragmentation and the representativeness of protected natural areas in the eastern Andes, Colombia. Biol. Conserv. 113, 245–256. doi: 10.1016/S0006-3207(02)00359-2
Aryal, D., De Jong, B. H. J., Gaona, S. O., Mendoza Vega, J., Olguín, L. E., López-Cruz, S. (2022). Fine wood decomposition rates decline with the age of tropical successional forests in southern Mexico: implications to ecosystem carbon storage. Ecosystems 25, 661–677. doi: 10.1007/s10021-021-00678-w
Báez, S., Cayuela, L., Macía, M. J., Álvarez-Dávila, E., Apaza-Quevedo, A., Arnelas, I., et al. (2022). FunAndes - A functional trait database of Andean plants. Sci. Data 9 (1), 511. doi: 10.1038/s41597-022-01626-6
Baker, T. R., Phillips, O. L., Malhi, Y., Almeida, S., Arroyo, L., Di Fiore, A., et al. (2004). Variation in wood density determines spatial patterns in Amazonian forest biomass. Glob. Chang. Biol. 10, 545–562. doi: 10.1111/j.1365-2486.2004.00751.x
Baraloto, C., Timothy Paine, C. E., Poorter, L., Beauchene, J., Bonal, D., Domenach, A. M., et al. (2010). Decoupled leaf and stem economics in rain forest trees. Ecol. Lett. 13 (11), 1338–1347. doi: 10.1111/j.1461-0248.2010.01517.x
Beltrán, H., Cataño, J. I. (2014). Caracterización de invasiones de Ulex europaeus L. de diferentes edades como herramienta para la restauración ecológica de bosques altoandinos, Colombia. Biota. Colombiana. 15 (2), 3–26.
Brown, S., Lugo, A. (1990). Tropical secondary forests. J.Trop. Ecol. 6 (1), 1–32. doi: 10.1017/S0266467400003989
Calbi, M., Clerici, N., Borsch, T., Brokamp, G. (2020). Reconstructing long term high andean forest dynamics using historical aerial imagery: A case study in Colombia. Forests 11 (88), 788. doi: 10.3390/f11080788
Calbi, M., Fajardo-Gutiérrez, F., Posada, J. M., Lücking, R., Brokamp, G., Borsch, T. (2021). Seeing the wood despite the trees: Exploring human disturbance impact on plant diversity, community structure, and standing biomass in fragmented high Andean forests. Ecol Evol. 11(5), 2110–2172. doi: 10.1002/ece3.7182
Castillo-Figueroa, D. (2021). Carbon cycle in tropical upland ecosystems: a global review. Web Ecol. 21, 109–136. doi: 10.5194/we-21-109-2021
Chao, K.-J., Phillips, O. L., Gloor, E., Monteagudo, A., Torres-Lezama, A., Martínez, R. V. (2008). Growth and wood density predict tree mortality in Amazon forests. J. Ecol. 96, 281–292. doi: 10.1111/j.1365-2745.2007.01343.x
Chave, J., Andalo, C., Brown, S., Cairns, M. A., Chambers, J. Q., Eamus, J. Q., et al. (2005). Tree allometry and improved estimation of carbon stocks and balance in tropical forests. Oecologia 145, 87–99. doi: 10.1007/s00442-005-0100-x
Chave, J., Coomes, D., Jansen, S., Lewis, S. L., Swenson, N. G., Zanne, A. E. (2009). Towards a worldwide wood economics spectrum. Ecol. Lett. 12, 351–366. doi: 10.1111/j.1461-0248.2009.01285.x
Chave, J., Muller-Landau, H. C., Baker, T. R., Easdale, T. A., ter Steege, H., Web, C. O. (2006). Regional and phylogenetic variation of wood density across 2456 neotropical tree species. Ecol. Appl. 16, 2356–2367. doi: 10.1890/1051-0761(2006)016[2356:rapvow]2.0.co;2
Chazdon, R. ,. L. (2014). Second growth: The Promise of Tropical Forest Regeneration in an Age of Deforestation (Chicago: University of Chicago Press).
Chazdon, R. L., Broadbent, E. N., Rozendaal, D. M., Bongers, F., Zambrano, A. M., Aide, T. M., et al. (2016). Carbon sequestration potential of second-growth forest regeneration in the Latin American tropics. Sci. Adv. 2 (5), e1501639. doi: 10.1126/sciadv.1501639
Chua, S. C., Potts, M. D. (2018). The role of plant functional traits in understanding forest recovery in wet tropical secondary forests. Sci. Total. Environ. 642, 1252–1262. doi: 10.1016/j.scitotenv.2018.05.397
Clark, D. A., Brown, S., Kicklighter, D. W., Chambers, J. Q., Thomlinson, J. N., Holland, E. A. (2001). Net primary production in tropical forests: an evaluation and synthesis of existing field data. Ecol. Appl. 11, 371–384. doi: 10.1890/1051-0761(2001)011[0371:NPPITF]2.0.CO;2
Clerici, N., Rubiano, K., Abd-Elrahman, A., Posada, J., Escobedo, F. J. (2016). Estimating aboveground biomass and carbon stocks in Periurban Andean secondary forests using very high resolution imagery. Forests 7 (7), 138. doi: 10.3390/f7070138
Departamento Administrativo Nacional de Estadística (DANE) (2018) Censo Nacional de Población y Vivienda. Available at: https://www.dane.gov.co/index.php/estadisticas-por-tema/demografia-y-poblacion/censo-nacional-de-poblacion-y-vivenda-2018/donde-estamos (Accessed December 24, 2021).
Duque, A., Peña, M. A., Cuesta, F., González-Caro, S., Kennedy, P., Phillips, O. L., et al. (2021). Mature Andean forests as globally important carbon sinks and future carbon refuges. Nat. Commun. 12, 2138. doi: 10.1038/s41467-021-22459-8
Etter, A., Andrade, A., Saavedra, K., Amaya, P., Cortés, J., Arévalo, P. (2021). Ecosistemas Colombianos. Amenazas y riesgos (Bogotá: Pontificia Universidad Javeriana).
Etter, A., McAlpine, C., Possingham, H. (2008). A historical analysis of the spatial and temporal drivers of landscape change in Colombia since 1500. Ann. Assoc. Am. Geogr. 98 (1), 2–23. doi: 10.1080/00045600701733911
Etter, A., van Wyngaarden, W. (2000). Patterns of landscape transformations in Colombia, whit emphasis in the Andean region. Ambio 29, 432–439. doi: 10.1579/0044-7447-29.7.432
Etter, A., Villa, A. (2000). Andean forest and farming systems in part of the Easter Cordillera (Colombia). Mt. Res. Dev. 20 (3), 236–245. doi: 10.1659/0276-4741(2000)020[0236:AFAFSI]2.0.CO;2
Evans, J., Senft, J., Green, D. (2000). Juvenile wood effect in red alder: analysis of physical and mechanical data to delineate juvenile and mature wood zones. For. Prod. J. 50 (7-8), 75–87.
Fearnside, P. M. (1997). Wood density for estimating forest biomass in Brazilian Amazonia. For. Ecol. Manage. 90, 59–87. doi: 10.1016/S0378-1127(96)03840-6
Fortunel, C., Fine, P. V. A., Baraloto, C. (2012). Leaf, stem and root tissue strategies across 758 Neotropical tree species. Funct. Ecol. 26, 1153–1161. doi: 10.1111/j.1365-2435.2012.02020.x
Frazer, G. W., Canham, C. D., Lertzman, K. P. (1999). Gap Light analyser, version 2.0. Bull. Ecol. Soc Amer. 81, 191–197. doi: 10.2307/20168436
Galbraith, D., Malhi, Y., Aragão, L., Baker, T. (2014). The ecosystem dynamics of Amazonian and Andean forests. Plant Ecol. Divers. 7 (1–2), 1–6. doi: 10.1080/17550874.2013.826744
Garnier, E., Cortez, J., Billès, G., Navas, M.-L., Roumet, C., Debussche, M., et al. (2004). Plant functional markers capture ecosystem properties during secondary succession. Ecology 85, 2630–2637. doi: 10.1890/03-0799
Garnier, E., Navas, M. L., Grigulis, K. (2016). Plant functional diversity (NewYork: Oxford University Press).
González-Melo, A. (2021). Radial variations in wood density, and their implications for above-ground biomass estimations, in a tropical high-andean forest. Dendrobiology 86, 19–28. doi: 10.12657/denbio.086.003
Guisande, C., Heine, J., González-DaCosta, J., García-Roselló, E. (2014) RWizard software. Available at: http://www.ipez.es/rwizard (Accessed January 14, 2021).
Hammer, O., Harper, D., Ryan, P. D. (2001). PAST: Paleontological Statistics software package for education and data analysis. Palaentol. Electron. 4 (1), 1–9.
Henry, M., Bombelli, A., Trotta, C., Alessandrini, A., Birigazzi, L., Sola, G., et al. (2013). GlobAllomeTree: International platform for tree allometric equations to support volume, biomass and carbon assessment. iForest 6, 326–330. doi: 10.3832/ifor0901-006
Hernández-Camacho, J., Ortiz-Quijano, R., Walshburger, T., Hurtado-Guerra, A. (1992). “Estado de la biodiversidad en Colombia,” in La diversidad biológica de Iberoamérica. Ed. Halffter, G. (México: INECOL), 41–239.
Homeier, J., Leuschner, C. (2020). Factors controlling the productivity of tropical Andean forests: Climate and soil are more important than tree diversity. Biogeosciences 18, 1525–1541. doi: 10.5194/bg-18-1525-2021
Hubau, W., Lewis, S. L., Phillips, O. L., Affum-Baffoe, K., Beeckman, H., Cuní-Sanchez, A., et al. (2020). Asynchronous carbon sink saturation in African and Amazonian tropical forests. Nature 579, 80–87. doi: 10.1038/s41586-020-2035-0
Hurtado-M, A. B., Echeverry-Galvis, M. Á., Salgado-Negret, B., Muñoz, J. C., Posada, J. M., Norden, N. (2021). Little trace of floristic homogenization in peri-urban Andean secondary forests despite high anthropogenic transformation. J. Ecol. 109 (3), 1468–1478. doi: 10.1111/1365-2745.13570
Jacobsen, A. L., Pratt, R. B., Ewers, F. W., Davis, S. D. (2007). Cavitation resistance among 26 chaparral species of southern California. Ecol. Monogr. 77, 99–115. doi: 10.1890/05-1879
Kattge, J., Bönisch, G., Díaz, S., Lavorel, S., Prentice, I. C., Leadley, P., et al. (2020). TRY plant trait database – enhanced coverage and open access. Global Change Biol. 26, 119–188. doi: 10.1111/gcb.14904
King, D. A., Davies, S. J., Tan, S., Noor, N. S. M. D. (2006). The role of wood density and stem support costs in the growth and mortality of tropical trees. J. Ecol. 94, 670–680. doi: 10.1111/j.1365-2745.2006.01112.x
Kitajima, K. (1994). Relative importance of photosynthetic traits and allocation patterns as correlates of seedling shade tolerance of 13 tropical trees. Oecologia 98 (3–4), 419–428. doi: 10.1007/bf00324232
Kraft, N. J. B., Metz, M. R., Condit, R. S., Chave, J. (2010). The relationship between wood density and mortality in a global tropical forest dataset. New Phytol. 188, 1124–1136. doi: 10.1111/j.1469-8137.2010.03444.x
Larjavaara, M., Muller-Landau, H. C. (2012). Still rethinking the value of high wood density. Am. J. Bot. 99 (1), 165–168. doi: 10.3732/ajb.1100324
Leuschner, C., Moser, G., Bertsch, C., Röderstein, M., Hertel, D. (2007). Large altitudinal increase in tree root/shoot ratio in tropical mountain forests of Ecuador. Basic. Appl. Ecol. 8, 219–230. doi: 10.1016/j.baae.2006.02.004
Leuschner, C., Zach, A., Moser, G., Homeier, J., Graefe, S., Hertel, D., et al. (2013). “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. Eds. Bendix, J., Beck, E., Bräuning, A., Makeschin, F., Mosandl, R., Scheu, S. (Berlin: Springer), 117–139.
Lohbeck, M., Lebrija-Trejos, E., Martínez-Ramos, M., Meave, J. A., Poorter, L., Bongers, F. (2015b). Functional trait strategies of trees in dry and wet tropical forests are similar but differ in their consequences for succession. PloS One 10 (4), e0123741. doi: 10.1371/journal.pone.0123741
Lohbeck, M., Poorter, L., Martínez-Ramos, M., Bongers, F. (2015a). Biomass is the main driver of changes in ecosystem process rates during tropical forest succession. Ecology 96, 1242–1252. doi: 10.1890/14-0472.1
López-Camacho, R., Quintero-Gómez, A., Amado-Ariza, S. (2020). Rasgos funcionales de la madera de tres bosques en Colombia: Bosque Seco, Andino y Alto-Andino. Ciênc. Florest. 30 (3), 856–872. doi: 10.5902/1980509839184
Lowe, S., Browne, M., Boudjelas, S., Poorter, M. (2000). 100 of the world’s worst invasive alien species A selection from the global invasive species database (Auckland: Hollands Printing).
Mahecha, G. E., Ovalle, A., Camelo, D., Rozo, A., Barrero, D. (2012). Vegetación del territorio CAR. 2nd ed. (Bogotá: Corporación Autónoma Regional de Cundinamarca, CAR).
Malhi, Y., Girardin, C. A., Goldsmith, G. R., Doughty, C. E., Salinas, N., Metcalfe, D. B., et al. (2017). The variation of productivity and its allocation along a tropical elevation gradient: a whole carbon budget perspective. New Phytol. 214, 1019–1032. doi: 10.1111/nph.14189
Markesteijn, L., Poorter, L., Bongers, F., Paz, H., Sack, L. (2011). Hydraulics and life history of tropical dry forest tree species: coordination of species’ drought and shade tolerance. New Phytol. 191, 480–495. doi: 10.1111/j.1469-8137.2011.03708.x
Martínez-Cabrera, H. I., Jones, C. S., Espino, S., Schenk, H. J. (2009). Wood anatomy and wood density in shrubs: responses to varying aridity along transcontinental transects. Am. J. Bot. 96, 1388–1398. doi: 10.3732/ajb.0800237
Martínez-Cabrera, H. I., Schenk, H. J., Cevallos-Ferriz, R. S., Jones, C. S. (2011). Integration of vessel traits, wood density, and height in angiosperm shrubs and trees. Am. J. Bot. 98, 915–922. doi: 10.3732/ajb.1000335
Meinzer, F. C., Campanello, P. I., Domec, J. C., Gatti, M. G., Goldstein, G., Villalobbos-Vega, R., et al. (2008). Constraints on physiological function associated with branch architecture and wood density in tropical forest trees. Tree Physiol. 28, 1609–1617. doi: 10.1093/treephys/28.11.1609
Mendoza, J. E., Etter, A. (2002). Multitemporal analysis, (1940–1996) of land cover changes in the southwestern Bogotá Highplain. Landsc. Urban. Plan. 59, 147–158. doi: 10.1016/S0169-2046(02)00012-9
Moles, A. T., Ackerly, D. D., Webb, C. O., Tweddle, J. C., Dickie, J. B., Pitman, A. J., et al. (2005). Factors that shape seed mass evolution. Proc. Natl. Acad. Sci. U.S.A. 102 (30), 10540–10544. doi: 10.1073/pnas.0501473102
Moles, A. T., Falster, D. S., Leishman, M. R., Westoby, M. (2004). Small-seeded species produce more seeds per square metre of canopy per year, but not per individual per lifetime. J. Ecol. 92 (3), 384–396. doi: 10.1111/j.0022-0477.2004.00880.x
Montañez, G., Arcila, O., Pacheco, J. C. (1994). Hacia Dónde va la Sabana de Bogotá: Modernización, Conflicto, Ambiente y Sociedad (Bogotá: Universidad Nacional de Colombia).
Moser, G., Hertel, D., Leuschner, C. (2007). Altitudinal change in LAI and stand leaf biomass in tropical montane forests: A transect study in Ecuador and a pan-tropical meta-analysis. Ecosystems 10, 924–935. doi: 10.1007/s10021-007-9063-6
Moser, G., Leuschner, C., Hertel, D., Graefe, S., Soethe, N., Iost, S. (2011). Elevation effects on the carbon budget of tropical mountain forests (S Ecuador): the role of the belowground compartment. Glob. Chang. Biol. 17 (6), 2211–2226. doi: 10.1111/j.1365-2486.2010.02367.x
Muller-Landau, H. C. (2004). Interspecific and inter-site variation in wood specific gravity of tropical trees. Biotropica 36, 20–32. doi: 10.1111/j.1744-7429.2004.tb00292.x
Myers, N., Mittermeier, R. A., Mittermeier, C. G., da Fonseca, G. A. B., Kent, J. (2000). Biodiversity hotspots for conservation priorities. Nature 403, 853–885. doi: 10.1038/35002501
Myster, R. W. (2021). “Introduction,” in The andean cloud forest. Ed. Myster, R. W. (Cham: Springer), 1–23.
Nogueira, E. M., Fearnside, P. M., Nelson, B. W. (2008). Normalization of wood density in biomass estimates of Amazon forests. For. Ecol. Manage. 256 (5), 990–996. doi: 10.1016/j.foreco.2008.06.001
Oliveira, L. Z., Uller, H. F., Klitzke, A. R., Eleotério, J. R., Vibrans, A. C. (2019). Towards the fulfillment of a knowledge gap: Wood densities for species of the subtropical Atlantic Forest. Data 4 (3), 104. doi: 10.3390/data4030104
Orsham, G. (1989). “Shrubs as a growth form,” in Biology and utilization of shrubs. Ed. McKell, C. M. (New York: Academic Press), 249–265.
Osunkoya, O. O., Sheng, T. K., Mahud, N., Damit, N. (2007). Variation in wood density, wood water content, stem growth and mortality among twenty-seven tree species in a tropical rainforest on Borneo Island. Austral Ecol. 32, 191–201. doi: 10.1111/j.1442-9993.2007.01678.x
Padilha, D. L., Marco Júnior, P. A. (2018). Gap in the woods: Wood density knowledge as impediment to develop sustainable use in Atlantic Forest. For. Ecol. Manage. 424, 448–457. doi: 10.1016/j.foreco.2018.05.012
Patiño, S., Lloyd, J., Paiva, R., Baker, T., Quesada, C., Mercado, L., et al. (2009). Branch xylem density variations across the Amazon Basin. Biogeosciences 6, 545–568. doi: 10.5194/bg-6-545-2009
Pearson, K. (1901). On lines and planes of closest fit to systems of points in space. Lond. Edinb. Dublin. Philos. mag. J. Sci. 2, 559–572. doi: 10.1080/14786440109462720
Pérez, M. C., Díaz, J. J. (2010). Estimación del Carbono Contenido en la Biomasa Forestal Aérea de dos Bosques Andinos en los Departamentos de Santander y Cundinamarca (Bogotá, Colombia: Universidad Distrital Francisco José de Caldas).
Pérez-Harguindeguy, N., Garnier, E., Lavorel, S., Poorter, H., Jaureguiberry, P., Cornwell, W. K., et al. (2013). New handbook for standardised measurement of plant functional traits worldwide. Aust. J. Bot. 61, 167–234. doi: 10.1071/BT12225
Petersson, H., Holm, S., Ståhl, G., Alger, D., Fridman, J., Lehtonen, A., et al. (2012). Individual tree biomass equations or biomass expansion factors for assessment of carbon stock changes in living biomass—a comparative study. For. Ecol. Manag. 270, 78–84. doi: 10.1016/j.foreco.2012.01.004
Petit, R. J., Hampe, A. (2006). Some evolutionary consequences of being a tree. Annu. Rev. Ecol. Evol. Syst. 37, 187–214. doi: 10.1146/annurev.ecolsys.37.091305.110215
Poorter, L., Craven, D., Jakovac, C. C., van der Sande, M. T., Amissah, L., Bongers, F., et al. (2021a). Multidimensional tropical forest recovery. Science 374 (6573), 1370–1376. doi: 10.1126/science.abh3629
Poorter, L., Rozendaal, D., Bongers, F., Almeida, J. S., Álvarez, F. S., Andrade, J. L., et al. (2021b). Functional recovery of secondary tropical forests. Proc. Natl. Acad. Sci. U.S.A. 118 (49), e2003405118. doi: 10.1073/pnas.2003405118
Poorter, L., Rozendaal, D. M. A., Bongers, F., Almeida-Cortes, J. S., Almayda Zambrano, A. M., Álvarez, F. S., et al. (2019). Wet and dry tropical forests show opposite successional pathways in wood density but converge over time. Nat. Ecol. Evol. 3, 928–934. doi: 10.1038/s41559-019-0882-6
Poorter, L., Wright, S. J., Paz, H., Ackerly, D. D., Condit, R., Ibarra-Manríquez, G., et al. (2008). Are functional traits good predictors of demographic rates? Evidence from five neotropical forests. Ecology 89, 1908–1920. doi: 10.1890/07-0207.1
Riesco, G., Encinas, J. I., de Paula, J. E. (2019). Densidad de la madera de 59 especies del orden Sapindales procedentes de bosques naturales brasileños. Madera. y Bosques. 25 (2), e2521817. doi: 10.21829/myb.2019.2521817
Roitman, I., Bustamante, M. M. C., Haidar, R. F., Shimbo, J. Z., Abdala, G. C., Eiten, G., et al. (2018). Optimizing biomass estimates of savanna woodland at different spatial scales in the Brazilian Cerrado: Re-evaluating allometric equations and environmental influences. PloS One 13 (8), e0196742. doi: 10.1371/journal.pone.0196742
Rubiano, K., Clerici, N., Norden, N., Etter, A. (2017). Secondary forest and shrubland dynamics in a highly transformed landscape in the Northern Andes of Colombia, (1985–2015). Forests 8, 1–17. doi: 10.3390/f8060216
Rundel, P. W. (1991). “Shrub life forms,” in Responses of plants to multiple stresses. Eds. Mooney, H. A., Winner, W. E., Pell, E. J., Chu, E. (New York: Academic Press), 345–370.
Santiago, L., De Guzman, M., Baraloto, C., Vogenberg, J. E., Brodie, M., Hérault, B., et al. (2018). Coordination and trade-ofs among hydraulic safety, efciency and drought avoidance traits in Amazonian rainforest canopy tree species. New Phytol. 218, 1015–1024. doi: 10.1111/nph.15058
Scholz, F. G., Bucci, S. J., Goldstein, G., Meinzer, F. C., Franco, A. C., Miralles-Wilhelm, F. (2007). Biophysical properties and functional significance of stem water storage tissues in Neotropical savanna trees. Plant Cell Environ. 30, 236–248. doi: 10.1111/j.1365-3040.2006.01623.x
Sierra, C. A., del Valle, J. I., Orrego, A., Moreno, F. H., Harmon, M. E., Zapata, M., et al. (2007). Total carbon stocks in a tropical forest landscape of the Porce region Colombia. For. Ecol. Manage. 243, 299–309. doi: 10.1016/j.foreco.2007.03.026
Stegen, J. C., Swenson, N. G., Valencia, R., Enquist, B. J., Thompson, J. (2009). Above-ground forest biomass is not consistently related to wood density in tropical forests. Glob. Ecol. Biogeogr. 18, 617–625. doi: 10.1111/j.1466-8238.2009.00471.x
Stratton, L., Goldstein, G., Meinzer, F. C. (2000). Stem water storage and efficiency of water transport: their functional significance in a Hawaiian dry forest. Plant Cell Environ. 23m, 99–106. doi: 10.1046/j.1365-3040.2000.00533.x
Sugden, A. M. (1986). The Montane vegetation and flora of Margarita Island, Venezuela. J. Arnold. Arbor. 67, 187–232. doi: 10.5962/p.185936
Suzuki, E. (1999). Diversity in specific gravity and water content of wood among Bornean tropical rain forest trees. Ecol. Res. 14, 211–224. doi: 10.1046/j.1440-1703.1999.143301.x
Tanner, E. V. J., Vitousek, P. M., Cuevas, E. (1998). Experimental investigation of nutrient limitation of forest growth on wet tropical mountains. Ecology 79, 10–22. doi: 10.1890/0012-9658(1998)079[0010:EIONLO]2.0.CO;2
Valverde-Barrantes, O. J., Freschet, G. T., Roumet, C., Blackwood, C. B. (2017). A worldview of root traits: the influence of ancestry, growth form, climate and mycorrhizal association on the functional trait variation of fine-root tissues in seed plants. New Phytol. 215 (4), 1562–1573. doi: 10.1111/nph.14571
Van Gelder, H. A., Poorter, L., Sterck, F. J. (2006). Wood mechanics, allometry, and life-history variation in a tropical rain forest tree community. New Phytol. 171, 367–378. doi: 10.1111/j.1469-8137.2006.01757.x
Vargas, O., Léon, O., Díaz, A. (2009). Restauración ecológica en zonas invadidas por retamo espinoso y plantaciones forestales de especies exóticas (Bogotá: Universidad Nacional de Colombia).
Vasquez-Valderrama, M. Y., Solorza-Bejarano, J. (2018). Agrupación funcional de especies vegetales para la restauración ecológica de ecosistemas de montaña, Bogotá, Colombia. Colombia. forestal. 21 (1), 5–17. doi: 10.14483/2256201X.11730
Vitousek, P. M., Matson, P. A., Turner, D. R. (1988). Elevational and age gradients in hawaiian montane rainforest: foliar and soil nutrients. Oecologia 77, 565–570. doi: 10.1007/BF00377275
Warton, D. I., Duursma, R. A., Falster, D. S., Taskinen, S. (2012). Smart 3-an R package for estimation and inference about allometric lines. Methods Ecol. Evol. 3, 257–259. doi: 10.1111/j.2041-210X.2011.00153.x
Wheeler, E. A., Baas, P., Rodgers, S. (2007). Variations in wood anatomy: a global analysis based on the Insidewood database. IAWA. J. 28, 229–258. doi: 10.1163/22941932-90001638
Williamson, G. B. (1984). Gradients of wood gravity in trees. Bull. Torrey. Bot. Club. 111, 51–55. doi: 10.2307/2996210
Williamson, G. B., Wiemann, M. C. (2010). Measuring Wood specific gravity … correctly. Am. J. Bot. 97 (3), 519–524. doi: 10.3732/ajb.0900243
Wright, S. J., Kitajima, K., Kraft, N. J., Reich, P. B., Wright, I. J., Bunker, D. E., et al. (2010). Functional traits and the growth-mortality trade-off in tropical trees. Ecology 91 (12), 3664–3674. doi: 10.1890/09-2335.1
Wright, S. J., Muller-Landau, H. C., Condit, R., Hubbel, S. P. (2003). Gap-dependent recruitment, realized vital rates, and size distributions of tropical trees. Ecology 84, 3174–3185. doi: 10.1890/02-0038
Zanne, A. E., Lopez-Gonzalez, G., Coomes, D. A., Ilic, J., Jansen, S., Lewis Simon, L., et al. (2009). Data from: Towards a worldwide wood economics spectrum [Data]. Dryad Digital Repository. doi: 10.5061/dryad.234
Keywords: montane forest, carbon cycle, Colombian Andes, functional traits, community weighted mean
Citation: Castillo-Figueroa D, González-Melo A and Posada JM (2023) Wood density is related to aboveground biomass and productivity along a successional gradient in upper Andean tropical forests. Front. Plant Sci. 14:1276424. doi: 10.3389/fpls.2023.1276424
Received: 11 August 2023; Accepted: 18 October 2023;
Published: 09 November 2023.
Edited by:
Barbara Vento, CONICET Mendoza, ArgentinaReviewed by:
Elivane Salete Capellesso, Federal University of Paraná, BrazilJoice Klipel, Federal University of Rio Grande do Sul, Brazil
Copyright © 2023 Castillo-Figueroa, González-Melo and Posada. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Juan M. Posada, juan.posada@urosario.edu.co; Dennis Castillo-Figueroa, dennis.castillo@urosario.edu.co