Skip to main content
Advertisement
Browse Subject Areas
?

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

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Wood Specific Gravity Variations and Biomass of Central African Tree Species: The Simple Choice of the Outer Wood

  • Jean-François Bastin ,

    bastin.jf@gmail.com

    Affiliations Landscape Ecology and Plant Production Systems Unit, Université libre de Bruxelles, CP264-2, B-1050 Bruxelles, Belgium, BIOSE Department, Gembloux Agro-Bio Tech, Université de Liège, B-5030 Gembloux, Belgium, Ecole Régionale post-universitaire d’Aménagement et de gestion Intégrés des Forêts et Territoires tropicaux, Kinshasa, DR Congo

  • Adeline Fayolle,

    Affiliation BIOSE Department, Gembloux Agro-Bio Tech, Université de Liège, B-5030 Gembloux, Belgium

  • Yegor Tarelkin,

    Affiliation Landscape Ecology and Plant Production Systems Unit, Université libre de Bruxelles, CP264-2, B-1050 Bruxelles, Belgium

  • Jan Van den Bulcke,

    Affiliation UGCT, University Ghent Centre for X-ray Tomography, Proeftuinstraat 86, 9000 Ghent, Belgium

  • Thales de Haulleville,

    Affiliations BIOSE Department, Gembloux Agro-Bio Tech, Université de Liège, B-5030 Gembloux, Belgium, Laboratory for Wood Biology and Xylarium, Royal Museum for Central Africa, Tervuren, Belgium

  • Frederic Mortier,

    Affiliation UPR BSEF, CIRAD, Campus International de Baillarguet, F-34398 Montpellier, France

  • Hans Beeckman,

    Affiliation Laboratory for Wood Biology and Xylarium, Royal Museum for Central Africa, Tervuren, Belgium

  • Joris Van Acker,

    Affiliation UGCT, University Ghent Centre for X-ray Tomography, Proeftuinstraat 86, 9000 Ghent, Belgium

  • Adeline Serckx,

    Affiliations Ecole Régionale post-universitaire d’Aménagement et de gestion Intégrés des Forêts et Territoires tropicaux, Kinshasa, DR Congo, Behavioural Biology Unit, University of Liege, Liege, Belgium, Conservation Biology Unit, Royal Belgian Institute of Natural Sciences, Brussels, Belgium

  • Jan Bogaert,

    Affiliation BIOSE Department, Gembloux Agro-Bio Tech, Université de Liège, B-5030 Gembloux, Belgium

  • Charles De Cannière

    Affiliation Landscape Ecology and Plant Production Systems Unit, Université libre de Bruxelles, CP264-2, B-1050 Bruxelles, Belgium

Abstract

Context

Wood specific gravity is a key element in tropical forest ecology. It integrates many aspects of tree mechanical properties and functioning and is an important predictor of tree biomass. Wood specific gravity varies widely among and within species and also within individual trees. Notably, contrasted patterns of radial variation of wood specific gravity have been demonstrated and related to regeneration guilds (light demanding vs. shade-bearing). However, although being repeatedly invoked as a potential source of error when estimating the biomass of trees, both intraspecific and radial variations remain little studied. In this study we characterized detailed pith-to-bark wood specific gravity profiles among contrasted species prominently contributing to the biomass of the forest, i.e., the dominant species, and we quantified the consequences of such variations on the biomass.

Methods

Radial profiles of wood density at 8% moisture content were compiled for 14 dominant species in the Democratic Republic of Congo, adapting a unique 3D X-ray scanning technique at very high spatial resolution on core samples. Mean wood density estimates were validated by water displacement measurements. Wood density profiles were converted to wood specific gravity and linear mixed models were used to decompose the radial variance. Potential errors in biomass estimation were assessed by comparing the biomass estimated from the wood specific gravity measured from pith-to-bark profiles, from global repositories, and from partial information (outer wood or inner wood).

Results

Wood specific gravity profiles from pith-to-bark presented positive, neutral and negative trends. Positive trends mainly characterized light-demanding species, increasing up to 1.8 g.cm-3 per meter for Piptadeniastrum africanum, and negative trends characterized shade-bearing species, decreasing up to 1 g.cm-3 per meter for Strombosia pustulata. The linear mixed model showed the greater part of wood specific gravity variance was explained by species only (45%) followed by a redundant part between species and regeneration guilds (36%). Despite substantial variation in wood specific gravity profiles among species and regeneration guilds, we found that values from the outer wood were strongly correlated to values from the whole profile, without any significant bias. In addition, we found that wood specific gravity from the DRYAD global repository may strongly differ depending on the species (up to 40% for Dialium pachyphyllum).

Main Conclusion

Therefore, when estimating forest biomass in specific sites, we recommend the systematic collection of outer wood samples on dominant species. This should prevent the main errors in biomass estimations resulting from wood specific gravity and allow for the collection of new information to explore the intraspecific variation of mechanical properties of trees.

Introduction

Recent development of international programs aiming to reduce emissions from deforestation and forest degradation in the tropics (i.e., the REDD+;[1]) drew the attention of a wide scientific community to wood density. Indeed, to calculate the carbon budget of a forest, the biomass of all the trees composing the forest is estimated through the measurement of several structural parameters (i.e., the diameter at breast height, the wood density and the total height) in allometric models[2,3]. While the main source of error in biomass estimation remains the choice of the allometric model[25], potential errors on wood density measurement cannot be overlooked as it is recognized to be the second best predictor of the biomass of a tree[2,3,6]. Wood density, or more conventionally the wood specific gravity (WSG), i.e., wood oven-dry mass divided by its green volume, as used in most studies developing biomass allometric models[2,3], ranges from 0.1 to 1.5 g.cm-³ among tropical trees[7,8]. This variation is well conserved across the world’s tropical regions[9,10]. However, most species exhibit values close to the mean[9], oscillating between 0.56 and 0.63 g.cm-³[8]. The wood specific gravity is strongly conserved across phylogenies[8,11], so genus averages are often used for biomass estimation when the species are not identified in the field[12]. The variation in wood specific gravity has been shown to be greater among than within species[13,14], with substantial differences between light-demanding and shade-tolerant species[8,15]. Wood specific gravity integrates many aspects of wood mechanical properties[8,16] and is consequently often used as a proxy to understand the stature and functioning of tropical tree species[8,11,17]. The fast growth of light-demanding species in the early stages is often associated with the production of soft wood with low density[18,19]. Such a trade-off between growth and density is sometimes assumed to provide a competitive advantage but reduce tree longevity[20]. Conversely, many shade-tolerant species are believed to invest in denser wood and grow more slowly, but to persist longer in the understory. This higher wood specific gravity provides greater resistance to physical damage and potential pathogens[21,22], and reduces the investment in metabolism and conduction[23,24]. As a result, the mean wood specific gravity of a forest stand is often related to the successional status, with a lower wood specific gravity in secondary than in old-growth forests reported in central and southern America[25,26], in south-eastern Asia[27], and more recently in Central Africa[28].

The wood specific gravity also presents substantial variations within a tree, both along its vertical[10,29] and its radial profile[18,26,30]. Exceptional radial variations have been reported for large tropical pioneer species where the ratio between the outer and the inner wood can reach 4.3-fold[20]. This radial increase of wood density from pith-to-bark is often related to the variation of the tree growth rate during ontogeny[9]. As light-demanding species are particularly sensitive to light exposure, it is assumed that they need to invest in a denser wood to persist in a mature forest, i.e., when the competition for light is stronger[18,20]. By contrast, shade-tolerant species tend to show a decreasing trend through ontogeny, possibly due to gradual shift from the shaded understory to increased light exposure, and the consequent increase in metabolic activity[19]. Non-linear patterns (i.e., U-shaped) have also been described[30], reflecting the flexibility of growth in response to changing environmental conditions[18]. Indeed, several events of growth suppression and release are usually experienced across the lifetime of tropical trees[31].

Most studies investigating forest biomass variations globally[32], regionally[33,34] and locally[35,36] use average wood specific gravity values at the species or the genus level extracted from global repositories such as DRYAD[7]. However, it has been shown locally that the use of such repositories[7,8] can lead to an overestimation of the wood specific gravity of approximately 16% for the species community[37]. Studies focused on biomass often neglect both “within genera” and “within species” variations in wood specific gravity, the latter potentially depending on tree size and mechanical constraints, or on environmental conditions. In contrast, studies describing the variability of wood specific gravity, and stressing the potential consequences of such variations on the biomass of the tree[26,30], barely consider biomass estimations in the framework of the carbon budget of the forest. Neglecting the variability of wood specific gravity might be essentially problematic for species that prominently contribute to forest biomass; i.e., species reaching large dimensions[32,38], or those that are frequent and/or locally ‘dominant’[33,39].

In the present study we aimed (i) to extract pith-to-bark wood specific gravity profiles of the 120 wood cores collected on 14 ‘dominant’ tree species, (ii) to characterize the variance partition of wood specific gravity profiles using linear mixed models and (iii) to quantify errors and biases in biomass estimations when wood specific gravity is not fully measured or only extracted from global repositories. To our knowledge, this is the first study addressing the problem of wood specific gravity radial variation in the framework of biomass estimations, i.e., with a particular attention on both intercepting detailed radial profiles and selecting species contributing mainly to the biomass of the forest.

Material and Methods

Data collection

The study area was located at the southern edge of the Congo Basin, north of the Bateke plateau, in the Bandundu province of the Democratic Republic of the Congo (WGS 1984; 2°29’35” S, 16°30’5.5” E). Mean annual temperature and mean annual rainfall are 25°C and 1500 mm.year-1, respectively, with a long dry season occurring between June and August and a short dry season in February[40]. The land cover is characterized as a forest-savanna mosaic[36,41] and forest species composition is identified as typical for Moist Central African forests[42]. Based on 26 1-ha plots sampled in 2011 and 2012 we estimated variations in above-ground biomass (AGB) between 27 and 460 Mg.ha-1 from young secondary to old-growth forests[36] (Table A in S2 File). The biomass of each tree with a diameter at breast height (dbh, measured at 130 cm or 50 cm above any buttresses) greater than or equal to 10 cm was estimated from the measurement of the dbh, total height and WSG using the pantropical equation developed for moist forests[2]. Measured trees were identified up to the species-level in the field, and samples were deposited in the Herbarium and botanical African library of the Université libre de Bruxelles (references Bastin-Serckx, #1–474). We selected 14 species from 14 genera and 8 families (Table 1) foremost based on their contribution to forest biomass but also considering contrasted regeneration guilds (7 non-pioneer light demanding, 5 shade-tolerant and 2 swamp species) as categorized by Hawthorne[43]. Non-pioneer light demanding species will be referred to as light-demanding species. The 14 species account for 50.7% of the total AGB estimated from inventory data[36] (Table B in S2 File). According to the global wood density repository registered on DRYAD [7,8] (registered under the following DOI: http://dx.doi.org/10.5061/dryad.234/1), average values for the 14 species range between 0.28 and 0.95 g.cm-3, with a mean of 0.65 g.cm-3.

thumbnail
Table 1. Characteristics of the 14 study species including botanical family, regeneration guild sensu Hawthorne[43], number of trees sampled (n), WSG, DBH, sample quantity (n) of the trees sampled, as well as frequency (% of plot presence), density (# of stems per hectare) and aboveground biomass (AGB) in the Malebo study site, in the Democratic Republic of Congo[36].

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

A total of 200 individual trees were sampled in the field. Because the pith was not noticeable on 80 of the 200 cores scanned (S1 File), these were excluded from further analyses. Remaining samples covered a wide range of diameters (from 6.5 to 77 cm, i.e., pith-to-bark profile from 3.25 to 38.5 cm), with 3 to 14 sampled trees per species (Table 1). Wood cores were extracted manually using Häglof Pressler augers of 30 or 60 cm length and 5 mm diameter. No specific permissions were required for the extraction of wood cores because these constitute only superficial and non-destructive vegetative samples from non-endangered species. We measured the wood density at 8% moisture content (WD8%) along the radial profile of each wood core, i.e., from pith-to-bark, using an X-ray CT scanner built at the Ghent University Centre for X-ray Tomography (UGCT; http://www.ugct.ugent.be). Cores were scanned using a closed microfocus X-ray tube to obtain a profile resolution of 50 μm[44] and extracted as wood density profiles using the Fiji software 1.6.0_24. Core snapshots reconstructed with a cross section view allowed for the identification of the exact location of the pith (Figs A and B in S2 File). For each core, a 1-D microdensitometric profile from pith-to-bark was extracted and smoothed at 1-mm resolution using a moving average window (Fig 1).

thumbnail
Fig 1. Variation in wood density measured at 8% of moisture content (g.cm-³) along the distance to the pith (cm) for the 14 species investigated in Malebo, the Democratic Republic of the Congo.

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

Data analysis

We first used the linear mixed framework to identify variations in the radial profile[45]. Regeneration guild and distance to the pith were included as fixed effects whereas species and individual within species were treated as random effects (see Eq 1). This modeling approach accounts for dependencies between measurements on the same species and individuals. As variance components for random slopes cannot be easily integrated with other variance components[46], we only considered random intercepts. (equation 1) where WDlkij is the wood density (g.cm-³) for the ith 1-mm interval in the jth core of the kth species and the lth regeneration guild; dij is the radial distance from the ith interval to the pith of the jth core (in mm); β0 and β1 are the fixed effects; α0 and α1 are the random intercepts; and εijkl is the residual variance within individual. Model diagnostics, i.e., variance constancy, homogeneity and normality, were assessed graphically[47]. Model selection was undertaken using likelihood ratio tests against several reduced models (see Table 2 [48]), completed by Akaike criterion and Kenward-Roger’s analysis of variance for fixed effects (Table C in S2 File).

thumbnail
Table 2. Coefficients from Gaussian linear mixed models predicting the wood density along the radial profile for the 14 species investigated in Malebo, the Democratic Republic of Congo.

Coefficient estimates are provided for the fixed effects (at 95% of confidence interval). The contribution to the total variance is estimated for the random effect through the maximum of likelihood.

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

Then, we calculated the mean WSG of each tree accounting for the increasing surface of wood along the profile[49], hereafter referenced as the weighted-WSG. We used the conversion relationship developed by Sallenave[50] to convert WD8% in WSG: (equation 2) where D is the wood density measured at moisture content M; d is the correction of D when the moisture content varies by 1%; v is the variation of the wood volume when the moisture content varies by 1%; and S is the fiber saturation point. Average d, v and S for each species were obtained from literature data on wood physical and mechanical properties[5052]. Species not recorded were converted to WSG using the coefficients from the relationship between WSG and WD8% fitted on recorded species. This relationship is very stable and is comparable to previous work performed at a larger scale[14] with a slope of 0.860 and an intercept close to 0 (n = 94; r² = 0.987; P<0.001).

We tested the validity of the WSG derived from the scans by comparing the results with WSG measured from the water displacement method[13] for a random subset of 25 wood cores belonging to 10 of the 14 study species (Figure C in S2 File). Samples were oven-dried at 105°C for 48 to 72 h. Dry mass and dry volume displacement were measured on a 1 mg precision scale and converted in WSG using Eq 2. We found a systematic bias (~ 9%; R² = 0.99; Figure C in S2 File), and consequently applied a systematic correction on the 120 wood cores[44]. This correction has been performed using the parameters of the linear relationship between WSG measured from 3D X-rays and WSG measured by water displacement method (Figure C in S2 File).

We used paired t-test to determine any significant differences between the mean WSG of the 14 species between weighted and the global datasets, and we used linear regression analysis to identify potential biases.

The potential error in WSG and AGB resulting from ignoring WSG radial variations was assessed by comparing the weighted-WSG, used as the reference, with the WSG from partial samples[26], i.e., the inner-WSG (2 cm from the pith) and the outer-WSG (2 cm under the bark). We then calculated the slope between the inner-WSG and the outer-WSG and estimated the influence on the potential error in WSG at the individual, species, and regeneration guild level. These slopes indicate the importance of radial variations[53]. Slopes were scaled-up to meters (multiplied by 100) to simplify interpretation, such that they express the variation of WSG in g.cm-3 for a radial length of 1 meter[26].

All statistical analyses were performed within the R environment (R Development Team) using the lmer() function of the lme4 and lmerTest packages[54].

Results

Pith-to-bark profiles

We first examined the major variations in WSG. On Fig 1, we can see positive trends on 3 of the light-demanding species (Piptadeniastrum africanum, Uapaca guineensis and Pentaclethra eetveldeana) and negative trends on 2 light-demanding species (Klainedoxa gabonensis, Entandrophragma angolense), and 3 shade-bearer (Gilbertiodendron dewevrei, Polyalthia suaveolens and Strombosia pustulata). However, as only 3 samples were recorded on K.gabonensis, such trends need to be interpreted with caution. The rest of the species examined presented no obvious global trends.

Linear mixed models and variance partition

The best model, regarding AIC and LRT, accounted for regeneration guilds and nested effect between species and individuals (Table 2). This model explained 92% of the total variability in wood density. Species alone explained 45%, with an additional 36% explained both by species and regeneration guilds. Individuals alone explained 11%. The studied species showed contrasting mean WSG values, ranging from 0.35 for the light-wooded Macaranga staudtii to 0.93 g.cm-³ for the dense-wooded Klainedoxa gabonensis (Fig 2A). When aggregating species by regeneration guild (Fig 2B), we found that swamp species had the lowest mean WSG (~0.35 g.cm-³) while shade-tolerant species had WSG significantly above the average (~0.71 g.cm-³). Light-demanding species did not show any particular pattern (~0.63 g.cm-³) with values oscillating around the average of the 14 species (~0.65 g.cm-³).

thumbnail
Fig 2. Boxplots of the WSG and of the slope calculated between the outer-WSG and the inner-WSG for the 14 species (A) and the corresponding regeneration guilds (B) investigated in Malebo, the Democratic Republic of the Congo.

The color corresponds to the regeneration guild with swamp, light-demanding and shade-tolerant species respectively colored in light grey, grey and dark grey, respectively.

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

We then analyzed radial variations in WSG based on the comparison of the weighted-WSG along the radial profile and the WSG from DRYAD. We showed that average weighted-WSG for the 14 species of the 120 trees investigated was not significantly different from the average WSG calculated from DRYAD (t = -0.1539; df = 24.635; P = 0.88; Figure D in S2 File). However, we observed a substantial bias (intercept = 0.17; slope = 0.74), with the underestimation of WSG for light-wooded species and an overestimation for dense-wooded species when using WSG from DRYAD.

Errors and biases in biomass estimations

Finally, we quantified the errors in AGB estimations when inner-WSG, outer-WSG or WSG from DRYAD are used as a proxy of WSG (Fig 3), using as reference the AGB estimated with the weighted-WSG. In addition, we analyzed the relationship between the error in AGB, the WSG (from inner, outer and DRYAD) and the slope between the inner and the outer WSG for each single tree sampled. Both the use of inner-WSG and WSG from DRYAD revealed substantial errors (Fig 3A, 3C and 3D). The use of the inner-WSG presented an error strongly related to the slope (Fig 3D), with an underestimation of WSG for light-wooded species and an overestimation for dense-wooded species (Fig 3A). Interestingly, using the outer-WSG yielded much less error with no significant dependence on tree WSG or slope (Fig 3B and 3E). In the end, regarding all the 120 trees together, we found the error consequently averages out, whatever the proxy of WSG used (Fig 3F).

thumbnail
Fig 3. Relative error of the estimation of the AGB for each sampled tree in Malebo, the Democratic Republic of the Congo.

Relative errors were calculated using weighted-WSG as reference and using the inner-WSG (A,D), the outer-WSG (B,E) and the global-WSG (C) as estimators. The dependence of the relative error was tested against the absolute value of WSG and against the slope. The size of the dots is proportional to tree diameter. The final boxplots summarize the distribution of the errors according to each estimator (F).

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

Regarding species and guilds (Fig 4), we found that the AGB estimated from inner-WSG tended to be overestimated for shade-tolerant species while no general trend was identified for light-demanding species. In particular, the WSG and AGB of E. angolense and K. gabonensis tended to be overestimated while P. africanum and U. guineensis tended to be underestimated. Interestingly, no specific or general pattern was observed when using the outer-WSG. When using global-WSG, no general trends were found but the error at species level was considerable. Specifically, the use of global-WSG values for Dialium pachyphyllum led to 40% overestimation of the AGB, followed by Pycnanthus angolensis (-26%), E. angolense (-23%), and Plagiostyles africana (+20%).

thumbnail
Fig 4. Boxplots of the relative errors calculated for the 14 species (A) and the three regeneration guilds (B) sampled in Malebo, the Democratic Republic of the Congo.

Relative errors were calculated using weighted-WSG as reference and using the inner-WSG, the outer-WSG and the global-WSG as proxies. Swamp, light-demanding and shade-tolerant species are respectively colored in light grey, grey and dark grey, respectively.

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

Discussion

The first aim of our study was to examine the variation in WSG among and within Central African tree species, including radial variation within individuals. Based on an extensive analysis conducted at local scale, we found that species explained most wood density radial variance whereas individuals explained only a minor part. These results are consistent with previous observations in Panama[30] and Australia[55]. In addition, we found that the share of variance generally explained by species only (>70%) can actually be partitioned into 45% explained only by species, 36% shared with regeneration guilds, and 11% explained only by individuals. This observation highlights the strong relationship between wood density and regeneration guild and supports the idea of an integrated wood economic spectrum [8,17].

Shade-tolerant species systematically showed high WSG values with a decreasing trend along the radial profile (negative slopes, Fig 4). This observation supports assumptions about slow growth and investment in dense wood by shade-tolerant species, in particular during the early stages of their life[21,22]. The decrease of wood density from pith-to-bark is often attributed to increased light exposure and improved growth conditions in the later stages[19] with possible access to the canopy for some trees during ontogeny[31].

Wood density profiles with an increasing trend from pith-to-bark and a low WSG were not systematic for light-demanding species. Only three species out of seven showed an increasing profile: P.africanum, U.guineensis and Pentaclethra eetveldeana. In the study area, U.guineensis is the dominant species in young secondary forests colonizing the savannas[36] and P.eetveldeana and P.africanum both exhibit a distribution of diameters with a low relative abundance of small trees indicative of a regeneration shortage (Figure E in S2 File). These three species are unlikely to persist in old-growth forests, supporting the assumption that strictly light-demanding species present increasing trends of WSG from pith-to-bark[18,20]. The other light-demanding species showed a reverse-J shaped relationship between tree size and frequency (Figure E in S2 File), characteristic of species well established in a mature forest and tolerant to shade[56]. This confirms that the regeneration guild of a species may vary between sites[43] and should therefore be attributed with caution.

The second aim of our study was to assess the potential error in AGB estimations resulting from radial variations in WSG and from intraspecific variations due to differences between observed and global-WSG values. Most importantly, we found radial variations have minor consequences for AGB estimations as the outer-WSG is strongly correlated to the weighted-WSG at the individual (RMSE of 3.7% and error centered on 0), the species, and the regeneration guild level. This result is mainly explained by the 3-D geometrical properties of trees: the outer wood occupying more volume than the inner wood in a tree and being consequently much more representative of the weighted-WSG. For instance, for a tree with a dbh of 30 cm and assuming a perfect circular shape, the volume of outer wood (2 cm under the bark; 91 cm²) is 30-fold larger than the volume of inner wood (2 cm from the pith; 3 cm²). To evaluate the potential generalization of this result, we compared outer-WSG with weighted-WSG of published datasets from Costa Rica[26]. We found a strong correlation between the two (intercept = -0.03; slope = 1.05; Figure F in S2 File), which reinforce that the outer-WSG can safely be used as an estimator of the weighted-WSG, regardless of the location. Consequently, intersecting the full radial profile of WSG is not required to properly estimate the AGB of a tree. However, as tree growth is a dynamic process[57], further consideration should be given to the variations of the carbon stored within a forest through time. In that context, radial profiles of WSG acquired at a very high spatial resolution constitute a great opportunity to improve our understanding of the relationship between tree growth, tree dimension and tree anatomical properties[58].

We also found that the use of WSG from global repositories can strongly affect the estimation of AGB depending on the species concerned (up to 40% for Dialium pachyphyllum). Because radial variations do not explain such bias in AGB estimations, we suggest that this may result from intraspecific variation for widely distributed species occurring under various climatic regimes. WSG from global repositories must consequently be used with caution when estimating the AGB of these species. Consequently, to minimize the potential error in AGB estimations related to the WSG, we recommend the extraction of a superficial sample of wood (2 cm under the bark). Considering that a small proportion of tropical tree species disproportionately contributes to the regional stem abundance[59] and biomass[33], we particularly recommend the systematic collection of outer wood samples for these ‘biomass dominant’ species. This should offer new opportunities to study the intraspecific variation of WSG along bioclimatic gradients.

Supporting Information

S1 File. Dataset of the detailed pith-to-bark profiles.

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

(TXT)

S2 File. Supporting informations.

Detailed wood density profile illustration with X-ray (Figure A), Pith localization (Figure B), Xray and water displacement correlation (Figure C), Dryad vs. observed values (Figure D), Species local diametric structures (Figure E), Outer WSG vs. inner WSG wood as proxy of weighted WSG in Costa Rica (Figure F), Biomass inventories metadata (Table A), Species contribution to total biomass (Table B), Kenward-Rogers approximation (Table C).

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

(DOC)

Acknowledgments

We are grateful to the FRIA (FNRS) and the ERAIFT for support and funding, to the WWF Belgium and DRC, to MbouMonTour NGOs for logistical support. We also would like to thank J.Lejoly, L.Pauwels, J-.L.Doucet and J-.F.Gillet (ULg–Gembloux, Belgium) for advice and aid in botanical identifications, to J.Hébert and P.Thibaut (ULg–Gembloux, Belgium) for general advice and water displacement analysis on the samples and to Nicholas Rowe (CNRS, AMAP, France), Olivier Hardy and Marius Gilbert (ULB, Belgium) for their insightful comments. Finally, we are very grateful to the local community chiefs from Nkala and Mpelu villages, and particularly, to P.Eyiano, N.Mongama, Leboncoeur and Bovik for fieldwork assistance.

Author Contributions

Conceived and designed the experiments: J-FB AF JVdB JVA JB CDC. Performed the experiments: J-FB YT JVdB. Analyzed the data: J-FB AF YT FM. Contributed reagents/materials/analysis tools: J-FB JVdB JVA. Wrote the paper: J-FB AF YT JVdB TdH HB AS JB CDC.

References

  1. 1. Agrawal A, Nepstad D, Chhatre A. Reducing Emissions from Deforestation and Forest Degradation. Annu Rev Environ Resour. Annual Reviews; 2011;36: 373–396.
  2. 2. Chave J, Andalo C, Brown S, Cairns M a, Chambers JQ, Eamus D, et al. Tree allometry and improved estimation of carbon stocks and balance in tropical forests. Oecologia. 2005;145: 87–99. pmid:15971085
  3. 3. Chave J, Réjou-Méchain M, Búrquez A, Chidumayo E, Colgan MS, Delitti WBC, et al. Improved allometric models to estimate the aboveground biomass of tropical trees. Glob Chang Biol. 2014;20: 3177–3190. pmid:24817483
  4. 4. Chave J, Condit R, Aguilar S, Hernandez A, Lao S, Perez R. Error propagation and scaling for tropical forest biomass estimates. Philos Trans R Soc Lond B Biol Sci. 2004;359: 409–20. pmid:15212093
  5. 5. Molto Q, Rossi V, Blanc L. Error propagation in biomass estimation in tropical forests. Freckleton R, editor. Methods Ecol Evol. 2013;4: 175–183.
  6. 6. Vieilledent G, Vaudry R, Andriamanohisoa SFD, Rakotonarivo OS, Randrianasolo HZ, Razafindrabe HN, et al. A universal approach to estimate biomass and carbon stock in tropical forests using generic allometric models. Ecol Appl. Ecological Society of America; 2012;22: 572–83. Available: http://www.ncbi.nlm.nih.gov/pubmed/22611855 pmid:22611855
  7. 7. Zanne AE, Lopez-Gonzalez G, Coomes DA, Ilic J, Jansen S, Lewis SL, et al. Global wood density database. Dryad. Dryad. Identifier; 2009;235: 33. Available: http://hdl.handle.net/10255/dryad.235.
  8. 8. Chave J, Coomes D, Jansen S, Lewis SL, Swenson NG, Zanne AE. Towards a worldwide wood economics spectrum. Ecol Lett. 2009;12: 351–66. pmid:19243406
  9. 9. Détienne P, Chanson B. L’éventail de la densité du bois des feuillus : Comparaison entre différentes régions du monde. Bois forêts des Trop. CIRAD—Forêt; 1996; 19–30. Available: http://cat.inist.fr/?aModele=afficheN&cpsidt=2472307
  10. 10. Turner IM. The Ecology of Trees in the Tropical Rain Forest. Cambridge: Cambridge University Press; 2001.
  11. 11. Swenson NG, Enquist BJ. Ecological and evolutionary determinants of a key plant functional trait: wood density and its community-wide variation across latitude and elevation. Am J Bot. 2007;94: 451–9. pmid:21636415
  12. 12. Slik JWF. Estimating species-specific wood density from the genus average in Indonesian trees. J Trop Ecol. 2006;22: 481.
  13. 13. Maniatis D, Saint André L, Temmerman M, Malhi Y, Beeckman H. The potential of using xylarium wood samples for wood density calculations: a comparison of approaches for volume measurement. iForest—Biogeosciences For. 2011;4: 150–159.
  14. 14. Chave J, Muller-Landau HC, Baker TR, Easdale T a, ter Steege H, Webb CO. Regional and phylogenetic variation of wood density across 2456 Neotropical tree species. Ecol Appl. 2006;16: 2356–67. Available: http://www.ncbi.nlm.nih.gov/pubmed/17205910 pmid:17205910
  15. 15. King DA, Davies SJ, Supardi MNN, Tan S. Tree growth is related to light interception and wood density in two mixed dipterocarp forests of Malaysia. Funct Ecol. 2005;19: 445–453.
  16. 16. Panshin AJ, de Zeeuw C. textbook of wood technology. McGraw-Hill Book Co.; 1980.
  17. 17. Reich PB. The world-wide “fast-slow” plant economics spectrum: a traits manifesto. Cornelissen H, editor. J Ecol. 2014;102: 275–301.
  18. 18. Nock CA, Geihofer D, Grabner M, Baker PJ, Bunyavejchewin S, Hietz P. Wood density and its radial variation in six canopy tree species differing in shade-tolerance in western Thailand. Ann Bot. 2009;104: 297–306. pmid:19454592
  19. 19. Woodcock DW, Shier a. D. Wood specific gravity and its radial variations: the many ways to make a tree. Trees. 2002;16: 437–443.
  20. 20. Wiemann M, Williamson G. Extreme Radial Changes in Wood Specific Gravity in Some Tropical Pioneers. Wood Fiber Sci. 1988;20: 344–349. Available: http://swst.metapress.com/content/QR6U1843MR344J38
  21. 21. Chao K-J, Phillips OL, Gloor E, Monteagudo A, Torres-Lezama A, Martínez RV. Growth and wood density predict tree mortality in Amazon forests. J Ecol. 2008;96: 281–292.
  22. 22. King DA, Davies SJ, Tan S, Noor NSM. The role of wood density and stem support costs in the growth and mortality of tropical trees. J Ecol. 2006;94: 670–680.
  23. 23. McCulloh KA, Meinzer FC, Sperry JS, Lachenbruch B, Voelker SL, Woodruff DR, et al. Comparative hydraulic architecture of tropical tree species representing a range of successional stages and wood density. Oecologia. 2011;167: 27–37. pmid:21445684
  24. 24. Fournier M, Dlouhá J, Jaouen G, Almeras T. Integrative biomechanics for tree ecology: beyond wood density and strength. J Exp Bot. 2013;64: 4793–815. pmid:24014867
  25. 25. Urquiza-Haas T, Dolman PM, Peres CA. Regional scale variation in forest structure and biomass in the Yucatan Peninsula, Mexico: Effects of forest disturbance. For Ecol Manage. 2007;247: 80–90.
  26. 26. Plourde BT, Boukili VK, Chazdon RL. Radial changes in wood specific gravity of tropical trees: inter- and intraspecific variation during secondary succession. Anten N, editor. Funct Ecol. 2015;29: 111–120.
  27. 27. Slik JWF, Aiba S-I, Brearley FQ, Cannon CH, Forshed O, Kitayama K, et al. Environmental correlates of tree biomass, basal area, wood specific gravity and stem density gradients in Borneo’s tropical forests. Glob Ecol Biogeogr. 2010;19: 50–60.
  28. 28. Fayolle A, Picard N, Doucet J-L, Swaine M, Bayol N, Bénédet F, et al. A new insight in the structure, composition and functioning of central African moist forests. For Ecol Manage. 2014;329: 195–205.
  29. 29. Wassenberg M, Chiu H-S, Guo W, Spiecker H. Analysis of wood density profiles of tree stems: incorporating vertical variations to optimize wood sampling strategies for density and biomass estimations. Trees. 2014;29: 51–561.
  30. 30. Osazuwa-Peters OL, Wright SJ, Zanne AE. Radial variation in wood specific gravity of tropical tree species differing in growth-mortality strategies. Am J Bot. 2014;101: 803–811. pmid:24793318
  31. 31. Brienen RJW, Zuidema PA. Lifetime growth patterns and ages of Bolivian rain forest trees obtained by tree ring analysis. J Ecol. 2006;94: 481–493.
  32. 32. Slik JWF, Paoli G, McGuire K, Amaral I, Barroso J, Bastian M, et al. Large trees drive forest aboveground biomass variation in moist lowland forests across the tropics. Glob Ecol Biogeogr. 2013;22: 1261–1271.
  33. 33. Bastin J-F, Barbier N, Réjou-Méchain M, Fayolle A, Gourlet-Fleury S, Maniatis D, et al. Seeing Central African forests through their largest trees. Sci Rep. Nature Publishing Group; 2015;5: 13156. pmid:26279193
  34. 34. Gourlet-Fleury S, Rossi V, Rejou-Mechain M, Freycon V, Fayolle A, Saint-André L, et al. Environmental filtering of dense-wooded species controls above-ground biomass stored in African moist forests. J Ecol. 2011;99: 981–990.
  35. 35. Gourlet-Fleury S, Mortier F, Fayolle A, Baya F, Ouédraogo D, Bénédet F, et al. Tropical forest recovery from logging: a 24 year silvicultural experiment from Central Africa. Philos Trans R Soc Lond B Biol Sci. 2013;368: 20120302. pmid:23878332
  36. 36. Bastin J-F, Barbier N, Couteron P, Adams B, Shapiro A, Bogaert J, et al. Aboveground biomass mapping of African forest mosaics using canopy texture analysis: toward a regional approach. Ecol Appl. 2014;24: 1984–2001. Available: http://www.esajournals.org/doi/abs/10.1890/13-1574.1
  37. 37. Ramananantoandro T, Rafidimanantsoa HP, Ramanakoto MF. Forest aboveground biomass estimates in a tropical rainforest in Madagascar: new insights from the use of wood specific gravity data. J For Res. 2015;26: 47–55.
  38. 38. Chave J, Riera B, Dubois M-A, Riéra B. Estimation of biomass in a neotropical forest of French Guiana : spatial and temporal variability. J Trop Ecol. CUP; 2001;17: 79–96.
  39. 39. Fauset S, Johnson O., Gloor M, Baker TRTR, Monteagudo A, Brienen RJW, et al. Hyperdominance in Amazonian forest carbon cycling. Nat Commun. Nature Publishing Group; 2015;6: 6857. pmid:25919449
  40. 40. Vancutsem C, Pekel J, Kibambe L, Blaes X, de Waseige C, Defourny F. République démocratique du Congo—occupation du sol. Carte Géographique. Presses Universitaire de Louvain. Bruxelles, Belgique. 2006.
  41. 41. Serckx A, Huynen M-C, Bastin J-F, Hambuckers A, Beudels-Jamar RC, Vimond M, et al. Nest grouping patterns of bonobos (Pan paniscus) in relation to fruit availability in a forest-savannah mosaic. PLoS One. 2014;9: e93742. pmid:24695637
  42. 42. Fayolle A, Swaine MD, Bastin J-F, Bourland N, Comiskey JA, Dauby G, et al. Patterns of tree species composition across tropical African forests. J Biogeogr. 2014;41: 2320–2331.
  43. 43. Hawthorne WD. Ecological Profiles of Ghanaian Forest Trees [Internet]. Oxford Forestry Institute, Department of Plant Sciences, University of Oxford; 1995. Available: http://books.google.be/books/about/Ecological_Profiles_of_Ghanaian_Forest_T.html?id=qYvKAAAACAAJ&pgis=1
  44. 44. De Ridder M, Van den Bulcke J, Vansteenkiste D, Van Loo D, Dierick M, Masschaele B, et al. High-resolution proxies for wood density variations in Terminalia superba. Ann Bot. 2011;107: 293–302. pmid:21131386
  45. 45. Zuur A, Ieno E, Walker N, Saveliev A, Smith G. Mixed effects models and extensions in ecology with R [Internet]. Survival Analysis, Edition …. Springer Science & Business Media; 2009. Available: http://link.springer.com/content/pdf/10.1007/978-1-4757-2739-5.pdf
  46. 46. Nakagawa S, Schielzeth H. A general and simple method for obtaining R 2 from generalized linear mixed-effects models. O’Hara RB, editor. Methods Ecol Evol. 2013;4: 133–142.
  47. 47. Pinheiro J, Bates D. Mixed-Effects Models in S and S-PLUS. Springer-V. New York; 2000.
  48. 48. Pinheiro J, Bates D. Mixed effects models in S and S-PLUS [Internet]. Springer-V. New York; 2000. Available: http://books.google.com/books?hl=en&lr=&id=RFDe_BKxvRIC&oi=fnd&pg=PR7&dq=Mixed+effect+models+in+S+and+S-PLUS&ots=mOsYED2TK5&sig=kxgG9b-1AAQnWTvjP6-S705WFHc
  49. 49. Williamson GB, Wiemann MCM. Measuring wood specific gravity correctly. Am J Bot. 2010;97: 519–524. pmid:21622413
  50. 50. Sallenave P. Propriétés Physiques et Mécaniques des Bois Tropicaux. Deuxième Supplément. CTFT. Nogent sur Marne, France.; 1971.
  51. 51. Sallenave P. Propriétés Physiques et Mécaniques des Bois Tropicaux de l’Union Française. CTFT. Nogent sur Marn, France; 1955.
  52. 52. Sallenave P. Propriétés Physiques et Mécaniques des Bois Tropicaux. Premier Supplément. CTFT. Nogent sur Marne, France.; 1964.
  53. 53. Nogueira EM, Fearnside PM, Nelson BW. Normalization of wood density in biomass estimates of Amazon forests. For Ecol Manage. 2008;256: 990–996.
  54. 54. Bates D. Fitting linear mixed models in R. R News. 2005;5: 27–30.
  55. 55. Onoda Y, Richards AE, Westoby M. The relationship between stem biomechanics and wood density is modified by rainfall in 32 Australian woody plant species. New Phytol. 2010;185: 493–501. pmid:19925557
  56. 56. Morse DR, Lawton JH, Dodson MM, Williamson MH. Fractal dimension of vegetation and the distribution of arthropod body lengths. Nature. 1985;314: 731–733.
  57. 57. Stephenson NL, Das AJ, Condit R, Russo SE, Baker PJ, Beckman NG, et al. Rate of tree carbon accumulation increases continuously with tree size. Nature. Nature Publishing Group; 2014; pmid:24429523
  58. 58. Van den Bulcke J, Wernersson ELG, Dierick M, Van Loo D, Masschaele B, Brabant L, et al. 3D tree-ring analysis using helical X-ray tomography. Dendrochronologia. 2014;32: 39–46.
  59. 59. Ter Steege H, Pitman NCA, Sabatier D, Baraloto C, Salomão RP, Guevara JE, et al. Hyperdominance in the Amazonian tree flora. Science (80-). 2013;342: 1243092. Available: http://www.sciencemag.org/content/342/6156/1243092.short pmid:24136971