Estimations of tree annual biomass increments are used by a variety of studies related to forest productivity or carbon fluxes. Biomass increment estimations can be easily obtained from diameter surveys or historical diameter reconstructions based on tree rings' records. However, the biomass models rely on the assumption that wood density is constant. Converting volume increment into biomass also requires assumptions about the wood density. Wood density has been largely reported to vary both in time and between trees. In Norway spruce, wood density is known to increase with decreasing ring width. This could lead to underestimating the biomass or carbon deposition in bad years. The variations between trees of wood density have never been discussed but could also contribute to deviations. A modelling approach could attenuate these effects but will also generate errors.
Here a model of wood density variations in Norway spruce, and an allometric model of volume growth were developed. We accounted for variations in wood density both between years and between trees, based on specific measurements. We compared the effects of neglecting each variation source on the estimations of annual biomass increment. We also assessed the errors of the biomass increment predictions at tree level, and of the annual productivity at plot level.
Our results showed a partial compensation of the decrease in ring width in
bad years by the increase in wood density. The underestimation of the
biomass increment in those years reached 15 %. The errors related to the
use of an allometric model of volume growth were modest, around
Predicting trees' biomass increment is a key step in quantifying and understanding forest productivity. Considerable efforts have been spent in evaluating forest productivity and carbon sink strength (Ciais et al., 2008). While productivity has long referred to volume growth, amply used in the forest management and displayed in yield tables, the focus has recently switched to biomass, for its relationships with energy or carbon storage. Field-based estimations of biomass growth have a wide variety of applications, from forestry to carbon fluxes estimation, for example in comparison against eddy covariance (Barford et al., 2001; Rocha et al., 2006; Gough et al., 2008; Curtis et al., 2002; Ilvesniemi et al., 2009). Considerable efforts have been spent in estimating annual forest productivity in relation to climate fluctuations and forests' carbon sink strength (Richardson et al., 2010; Wu et al., 2014). The importance of having both annual resolution and high spatial coverage has been illustrated by numerous studies (e.g. Reichstein et al., 2007; Ciais et al., 2005; Beer et al., 2010). Several methods are used to estimate forest productivity and carbon sink: eddy covariance, modelling, or field-based estimations such as inventories or tree-ring studies. Tree-ring-based studies have the advantage of offering a large spatial covering, a potentially long timescale and also an annual resolution. They are therefore often used to produce reference annual biomass production estimations, to compare with other methods (Beck et al., 2011; Babst et al., 2014a) or to provide complementary information (Babst et al., 2013). However, several issues are associated with the use of tree-ring-based estimations and the estimation of their error remains a critical issue that is poorly documented (Nickless et al., 2011).
In the reconstruction of the annual productivity or of the above-ground carbon uptake from field-based studies, one limiting element is the estimation of the wood density variations (Babst et al., 2014a). Indeed, volume increment time series can be produced by a variety of methods, such as the reconstruction of the diameter growth, based on tree rings (Wirth et al., 2004; Rocha et al., 2006) or inventory reconstruction (Ohtsuka et al., 2007), but none of these methods provide information on the variation of wood density. Converting volume into biomass requires an estimation of the wood density, which is most likely based on literature and therefore related neither to site conditions, nor to trees' growth rate, as for example in Vilàet al. (2013). In the same manner, biomass equations implicitly rely on the use of an average and constant wood density despite the many evidences of substantial wood density variations. In both cases, wood density is considered constant in time, and equal between trees.
Wood density has however been acknowledged as a highly variable characteristic, and several major sources of annual density variations have been identified. Very high precision in the description of the wood density variations with new techniques (e.g. SilviScan; Evans, 1994) is possible but not widely available, while other techniques based on X-ray are rather time-consuming and thus not applied to forest productivity studies. Within-tree variations occur at distinct timescales (Jyske et al., 2008). Over medium or long timescales, annual wood density was proved to be related to ring age or to tree diameter, with higher values close to the pith in many species (Schweingruber, 1988). At an interannual scale, wood density variations can be substantial. There were several reports that (annual) ring density decreases with increasing ring width, for instance in Norway spruce (Bergqvist, 1998; Dutilleul et al., 1998; Lundgren, 2004; Bouriaud et al., 2005; Franceschini et al., 2010, 2013). Wood density was also proved to vary between trees (Wilhelmsson et al., 2002; Guilley et al., 2004), a fact which is never accounted for in studies using diameter surveys to produce biomass increment estimations.
The variations in wood density between trees and between years could compensate the variations in annual volume increment, or at least soften them. Recent studies have brought evidence of such compensation, proving that neglecting annual wood density fluctuations could lead to substantial errors or bias in estimating the biomass (Molto et al., 2013; Babst et al., 2014a). The errors generated by neglecting the variations in wood density have been considered small compared to those resulting from the volume increment estimation, but to our knowledge, such assumptions have never been tested and the consequences not documented.
To be properly quantified, the consequences of neglecting wood density fluctuations between years and between trees had to be tested using an integrated approach, whereby the errors of the density model are propagated and combined with those of the model for volume growth. Such chains can be decomposed, and the impact of each step was studied by modelling the steps into a single Markov chain Monte Carlo (MCMC) process (e.g. Molto et al., 2013). Analytical solutions to estimate the biomass estimation error, based e.g. on Taylor expansion can sometimes be determined, depending on the model's complexity. But the errors of biomass increment, obtained by subtracting subsequent estimations, are anyhow less predictable and particularly challenging at the plot level, when summing tree-level estimations (Nickless et al., 2011). The MCMC approach therefore appears to be the most suitable to estimate the biomass increment, where such estimations and the propagation of the errors from one model to another are done without assumptions.
Our study aimed to quantify the impact of density variations, both between years and between trees, on the estimations of annual biomass increment in Norway spruce (Picea abies), and compare it with the impact of volume increment estimation errors. The objectives were as follows: (i) to quantify and model the influence of annual radial growth variations on wood density, (ii) to quantify the consequences of annual and between-tree variations of wood density on biomass increment estimations, and (iii) to compare the errors related to wood density estimations to those of volume increment.
All samples analysed for this study were taken from the Wetzstein site near
the village of Lehesten in Thuringia, central Germany
(50
This study combines data from three successive samplings realized in this site: (i) stem analysis performed to quantify the relationship between breast-height radial growth and stem volume increment. This was achieved in connection with a biomass harvest of the four stands (see below); (ii) wood density measurements taken for selected harvest trees to establish a relation between ring-width and wood density variations; and (iii) a dendrochronological analysis of interannual growth variation of many trees carried out using micro-cores for scaling up to the plot scale. The volume increment and wood density measurements are used exclusively to develop models, while the micro-core sampling is used as an application to quantify and compare the errors of each model in this representative case study.
The stem volume increment model was fit based on a stem analysis realized on
22 trees – seven samples in the footprint stand W72 and five in each of the
additional stands (W15, W33, W116). Trees were selected to represent
seven/five DBH (diameter at breast height) classes, defined based on the
population of all inventoried trees (W15:
For the annual wood density (WD) measurements, wood discs were sampled at
breast height from trees, representing the lowest, the central, and the
highest diameter class in each of the four stands. This yielded a total of
12 sample trees, again representing the size and age range of Norway spruce
trees at the site. Two 1–2 cm wide slices from opposite radii were sawn from
the wood discs, for which wood density was measured by X-ray densitometry in
the densitometric laboratory of Krasnoyarsk, Russia (Walesch Electronics,
Switzerland), using the standard procedure described by Schweingruber (1988).
Longitudinal strips with a constant thickness of 1.2 mm were sawn, air
dried, and exposed to X-ray radiation for 1 h on a Kodak TL film using
standard exposure conditions: acceleration tension of 8.5 kV, flux intensity
of 15.0 mA, and a distance to the source of 3.5 m. Annual wood density (WD, kg m
The volume increment and WD models were applied together on an independent
set of trees sampled in 13 randomly placed inventory plots inside the
footprint stand W72. The plots were established within the context of the
project FORCAST (Rey and Jarvis, 2006). Thirty-one to 62 trees per plots (551
in total) with diameter varying from 8 to 51 cm (thus well within the range
of the sample trees) were sampled for historical diameter reconstruction
based on micro-cores. The micro-cores enabled the reconstruction of the past
growth over the last 10 years only, since these short cores are
Models of WD or annual volume increment were fit using both maximum likelihood methods and the MCMC approach. The structure of the two models was first determined using likelihood fits before being implemented in a Bayesian MCMC framework using WinBUGS 1.4 (Spiegelhalter et al., 2003), based on exactly the same data sets, using non-informative flat priors. The maximum-likelihood estimations were realized using the “nlme” package (version 3.1-102, Pinheiro et al., 2011) of R (R version 3.0.1, R Development Team, 2014).
Following recent publications on Norway spruce wood density (Franceschini et al., 2010, 2013), the diameter and the ring cambial age (as counted from the pith) were used as independent variables. The selection of the model was based both on the AIC (Akaike information criterion) and the examination of the residual distribution. Fixed and random tree-level effects were considered. The principle of parsimony was also followed in the model-building process, and random effect parameters were considered only if improvements were observed based on the likelihood ratio test.
Several candidate models were tested, as follows:
The annual volume increment was modelled as a non-linear function of ring
width and tree diameter, based on the annual estimations of volume growth
resulting from the detailed stem analysis. The model reflects the fact that,
for a given ring width, volume increment depends strongly on the current
size of the tree, here its diameter, mostly for geometrical reasons. The
taper was therefore not supposed to be constant in time, and the trends in
tree growth with age were directly absorbed in the model since the volume
increments resulted directly from the stem analysis measurements, not from
using models. Another specificity of this model was the specification of a
variance function in order to cope with the heteroscedasticity in the
errors. The resulting model is given in Eq. (4) and includes random
coefficients for the exponent
The residual
The micro-core data set was used as a concrete case study for estimating the consequences of wood density variations and comparing the errors resulting from the wood density and from the volume increment model. Both models were fit based on their specific data sets within the MCMC framework (code provided in Supplement S2); then the estimated parameters and the variance terms were applied to compute the biomass increment of the micro-cores trees, which represents an external set. The models were therefore fit using the same structure as that used in the likelihood method, the estimated parameters being further used to produce estimations of WD or annual volume increment on the micro-core trees. Having both the fitting and the application run in a single MCMC loop enables the propagation of the errors of each model.
The tree-level biomass increment estimations were the product of the WD and
the volume increment, then summed up to obtain stand level per hectare biomass
estimations, also based on the plot size. But according to the way the errors
could be accounted for, four different scenarios were distinguished:
In the baseline scenario (sc1), a constant wood density was used, set to be equal
to the average observed value across the data set (475 kg m Only the fixed part of the parameters In the second scenario (sc2), the annual wood density was held constant but the
volume increment included both the random tree-level variation and the
residual error. For tree where In the third scenario (sc3), the biomass increment was defined as the product
of the parametric estimations of both the WD and the annual volume
increment: here only the fixed part of the models was used to produce both
the WD and the volume increment estimations, while not accounting for random
effects or residual variance. This represents the most common and probable
use of such models, when no data are available for a calibration. Only the fixed part of the parameters was used. In the last scenario (sc4), a full error propagation was conducted: the random
and the residual errors of both the WD and the volume increment models were
used to produce the biomass increment estimation. having
Thus, four different biomass increment estimations were produced, according
to the density estimation and the error propagation, and their difference was
summed at plot level. In all the scenarios, volume increment was estimated
based on measured ring-width series and the historical diameter of the
trees.
The MCMC process generated posterior distributions of the model parameter
estimates, with their associated errors, and the estimations of the variance
of the random effects, based on the Metropolis–Hastings algorithm over
10
The (annual) ring wood density (WD) varied from 287 to 787 kg m
Variations in WD were mostly related to ring width with a linear correlation
of
Relation between annual wood density and cambial age (left) or ring width (right) at tree level. Two trees with very distinct average wood density were highlighted (dark grey/black).
The selection of the WD model resulted from the comparison of several models
based on independent variables such as ring width, cambial age, and diameter.
The models offered very comparable results (Table 1) although model 2 had a
greater root-mean-square error (RMSE) and bias. Using cambial age or diameter
as a second independent variable did not lead to significant differences in
the fit according to the likelihood ratio test (LRT). Nevertheless, models
differed in the ease of the convergence or in the sensitivity to initial
parameters provided. The exponent parameters
Fit statistics and parameters for the wood density models.
WD: (annual) wood density, RW: (annual) ring width,
The volume increment model was fit as a function of diameter and ring width,
with fixed and random tree-level effects, to a set of 22 trees. The
intercept was kept free after testing its significance using the LRT by
comparing models with intercept held constant or forced to 0. It appeared
that a free intercept increases the likelihood, while the estimated value of
the intercept was very realistic. The use of a weight function (constant
plus power) was also amply confirmed by the LRT (L.ratio
Observed and fitted annual volume increment model and standardized residuals of the volume increment model fit.
Provided that there was an overall decrease in wood density with increasing
ring width, a compensation of ring-width annual variability by wood density
was also probable. The ring-width series showed peak years of growth (e.g.
1967, 1989) or depressions (1976, 1983). In these years, the radial growth
was much more affected than the wood density, as suggested by the deviations
relative to the mean value calculated over the entire series length. The
deviations peaked in 1967 at
Comparison of the fixed parameters estimated for the wood density and the volume models, obtained by maximum likelihood and MCMC. Standard deviations are provided in brackets.
WD: (annual) wood density; RW: (annual) ring width; DBH: (annual)
breast-height diameter;
Left, comparison of biomass increment estimations for Norway
spruce trees growing in Wetzstein, based on constant density hypothesis vs.
actual wood density measurements; right, time course of the average ratio of
biomass increment estimations (actual over constant density) and time course
of the detrended mean ring width (spline smoothing, for illustration
purposes). The
The two models presented and fitted above were introduced in the Bayesian framework, with the same structure exactly and on the same data, and further re-fitted using the MCMC method. A comparison of the parameters estimated by both methods is presented in Table 2. Expectedly, the parameters were not exactly the same but very close, and the correlation between the predictions was very high.
When applied on the independent application set, the estimated wood density
varied from 278 to 541 kg m
Comparable results were obtained with the volume increment model: the contribution of the random effects and the inclusion of the residual variance inflated the prediction interval substantially (Fig. 4). Nevertheless, the relative prediction interval was substantially lower than that of the wood density: typically less than 40 % of the predicted value, compared to 60 % for WD.
Left, variations of the MCMC annual predictions and prediction intervals (95 %) of wood density and volume increment for one given tree, randomly chosen while accounting for different error sources: regression only/regression and random effects/regression, random effects and residual variance; right, distribution density of the relative prediction interval (expressed in percentage of the prediction) for all trees included in the simulation, according to the error sources included.
The annual variations of the predicted biomass increment resulting from considering a dynamic wood density were always smaller than predictions based on a constant density (Fig. 5). The prediction uncertainty was considerably higher when accounting for random effects on either the WD or the volume increment. The full error propagation (sc4) had a relative prediction uncertainty up to 60 % of the predicted value on average, occasionally reaching or overcoming 100 %. Constant density predictions had logically the lowest uncertainties (less than 10 %) since they only included the error from the volume increment estimation. Wood density had the greatest contribution to the prediction uncertainty, and mainly through the between-tree variations. The parametric estimation (sc3) had a prediction interval 4 times lower than the full error propagation prediction (sc4), showing an underestimation of the error made by considering the uncertainty related to the regression coefficients only.
At plot level, which is the aggregation of the tree-level predictions and
errors, the prediction errors tended to compensate each other since the
relative prediction intervals of the annual biomass production were smaller
than at tree-level (Fig. 6). Thus the interval of biomass production
estimates varied from
Annual biomass increment (posterior MCMC distribution) for one
given tree chosen as representative with its associated prediction error for
scenario 1 and 2
The variation between years in the prediction error was also very low (Fig. 6) despite contrasted ring widths. The error of the predictions based on regression errors only (sc1 and sc3) did not vary with increasing number of trees in the plot (Fig. 6). In contrast, the predictions error decreased slightly with increasing number of trees for the scenarios that used a (tree-level) random-effect term (sc2 and sc4).
Wood density was found to decrease when ring width increased, in agreement
with previous studies on Norway spruce (e.g. Lindström,
1996; Dutilleul et al., 1998). Despite the seemingly high correlation
between ring width and WD, the decrease of WD was not enough to compensate
the increase in ring width but contributed to attenuate its effects. The
order of magnitude of the WD variability was not – and by far – as large as
that of ring width. Hence, it is logical to find a moderate compensation
between radial growth and wood density variations even in extreme years such
as 1976: 15 % at plot level. Nevertheless, when the focus is put on key
years, such as years of climatic extremes, the measurements of WD is
necessary to avoid a systematic underestimation of the biomass increment or
carbon uptake. Climate is indeed probably the most important driver of WD
variations with influences at both inter- and intra-annual time steps (e.g.
Gindl et al., 2000; Bouriaud et al., 2005). These results are consistent with
those reported in Babst et al. (2014a) showing that accounting for the
variations in WD strongly improved the match between the tree-ring-based
above-ground wood biomass increment estimations and the seasonal CO
Plot-level annual relative prediction intervals of the biomass
increment as a function of the number of trees sampled in the plot for the
four scenarios
A constant value of wood density, such as that implicitly used in biomass equations, can generate systematic deviations because it only has a few chances to be equal to the mean density of the trees to which the model is applied. Even if a site-specific WD value is used, neglecting the radial increment of WD (i.e. the age-related trend) will also lead to an underestimation of the biomass increment. This source of error can unfortunately not be compensated by a larger sampling since it affects all the trees simultaneously. This has consequences not only for the annual productivity estimations but also for periodical productivity assessments, such as those conducted on permanent sample plots over a 5- or 10-year period.
Compensation of increased growth rate by a decrease in wood density was documented for Norway spruce but over a long timescale (Bontemps et al., 2013). The trends in radial growth and in WD reported for many species could lead to such deviations between the actual WD and the modelled or implicit WD. In this context, a local calibration would reduce such errors but cannot solve the problem of the variations between years and between trees.
The anticorrelation between ring width and wood density seems to be a general feature in Norway spruce according to the literature (e.g. Lindström, 1996; Dutilleul et al., 1998) but the phenomenon is not limited to this species (Babst et al., 2014a). The attenuation therefore probably occurs at a large scale. The between-tree variability in the relationship has also been reported in several studies and probably is a widespread feature with potentially large consequences on the error of annual biomass increment predictions, as demonstrated by this study. The fact that the trees used to assess both the wood density variations and to model the volume increment came from the same site as those used for the error estimations has ruled out the issues of using locally inappropriate models. Additional errors should be considered in practice when using models that may not be locally valid.
The inventory-based or tree-ring-based estimations of annual biomass production or carbon uptake are often used for comparisons against other methods such as remote sensing, vegetation models or eddy covariance (Beck et al., 2011; Bunn et al., 2013; Babst et al., 2014a). To be conclusive, the benchmarking however supposes that prediction errors are known or can be estimated. High prediction errors would invalidate the biometric approaches but the errors are not always accounted for. Analytical solutions are indeed not always available to estimate the errors of the allometric models, and their estimation remains very complex or based on assumptions. In the case of the biomass increment, the error results from the combination of several models, and the estimation is even more challenging. The use of the MCMC framework here avoids the cumbersome analytical approximations for prediction variances (e.g. Wutzler et al., 2008).
The prediction interval at plot level was on average between 20 and 40 % of the predicted biomass increment value. The uncertainty related to the regression parameters was only about 10 % for both models. Reduced variance may be inherent to the use of local trees and the Bayesian modelling method (Zapata-Cuartas et al., 2012) but these values are similar to those found by Nickless et al. (2011) for biomass estimations following a parametric approach – as opposed to the MCMC method used here. Unlike our results, this study did not however include the random tree-level variations, which appeared to be quite an important source of uncertainty. Indeed, accounting for random tree-level variations in the relation between wood density and ring width increased the prediction interval of the tree-level biomass increment drastically (i.e. decreased the prediction confidence), by a factor of 5. Further errors related to the residual non-explained variance, were, in comparison, very small. Consequently, the prediction interval of the biomass annual increment at plot level increased twofold by accounting for the random-tree effects. Hence, the contribution of WD to the prediction error of the biomass increment was much larger than that of the volume increment model.
The tree-level prediction error (in percentage of the prediction value) was found to be greater the error at plot level. Thus, compensation occurred at plot level when summing up trees predictions. We speculate that this compensation happens because the variations are centred by construction around zero and have both negative and positive values. This explains also why the mean prediction values were always unaffected by accounting for random effects. Hence, neglecting random effects affected the prediction interval more than the predictions themselves.
The relation between wood density, ring width, and cambial age were proven to fluctuate between trees sampled within a same stand for many species: oak (Guilley et al., 2004; Bergès et al., 2008), common beech (Bouriaud et al., 2004), Norway spruce (Mäkinen et al., 2002; Jaakola et al., 2005; Franceschini et al., 2010). For a given radial growth rate, the trees are building more or less biomass and so storing more or less carbon, according to the density of the wood.
This fluctuation is considered random because it cannot be attributed to a measurable factor. Random tree-level variations were nevertheless reported as a major source of wood density variation in a population (Zhang et al., 1994; Guilley et al., 2004; Bouriaud et al., 2004; Jaakola et al., 2005). It is often hypothesized to be related to genetics, although this is not proven. Provenance studies have brought some insight on it (Hylen, 1999; Rozenberg et al., 2004), but much of the determinism remains unknown. Other factors, such as crown development (Lindström, 1996), could also be invoked to explain this variation source in wood density.
The changes in silvicultural practices, whereby the focus is put on targeted individuals, further stress the importance of errors in tree-level estimations of biomass and biomass increments. The tree-level variations were the largest error source and showed that the inter-tree variations can be seen as a limitation to the tree-level biomass prediction. Despite the many evidences of tree-level random effects, this variation source was largely ignored. Our study proved that the between-tree variations in the relation between ring width and wood density – although within the same species – contributed the most to the uncertainty in the biomass increment predictions. The variations are hypothesized to follow a normal distribution (Lindström and Bates, 1990). Thus, at plot level, compensation is likely to occur. But this situation may not be true for all samples, and certain designs could generate additional biases in the biomass production estimations. In this study, all the trees in a plot were sampled. Other samples, for instance the selection of the biggest trees in a plot as classically done in dendrochronology, could lead to serious deviations as it could involve sampling faster-growing trees. Apart from the bias in productivity caused by a sample focusing on faster-growing trees (Nehrbass-Ahles et al., 2014), the productivity at stand level would probably generate an overestimation related to a decreased wood density as trees producing larger rings would be sampled. Another issue in using the tree-ring parameters (width and density) to produce annual productivity estimations is the presence of autocorrelation or carry-over effects in the series, which are reflected in the derived productivity estimations but are generally not observed in the measured or modelled carbon fluxes (Babst et al., 2014a, b; Rammig et al., 2015).
Apart from the climate, the two most-used variables used in modelling annual WD variations are ring width and ring (cambial) age. The relation between WD and radial growth was strong in our study and probably dominant in Norway spruce but may not be so for other species. In beech, for example, the relation between ring width and WD was shown to be weak (Bouriaud et al., 2004) and there was no clear trend in WD related to the age either. Several studies reported a lack of significant correlations between ring width and WD for Norway spruce (e.g. Dutilleul et al., 1998). The relative stability in annual WD values is not calling for a correction of the biomass increment in such situation. It is probable that variations in WD would affect the estimation of biomass increment in species for which a relationship with ring width was already observed such oaks (Zhang et al., 1993; Bergès et al., 2008) or larch (Karlman et al., 2005). The contribution to the error in the prediction of biomass production is however likely to be important.
Conversely to ring width, ring age was found to be only slightly influent on the annual wood density in Norway spruce. Ring age is often considered in density models for representing the age trend or for the variations observed near the pith – the juvenile versus mature wood transition (e.g. Franceschini et al., 2010). WD in Norway spruce has been shown to present an age-dependent trend from pith to bark (Dutilleul et al., 1998; Hylen, 1999; Mäkinen et al., 2002), apart from the juvenile wood effect. In our study, the juvenile effect was not included for simplicity (series were pruned to exclude the first 3 years) but also because rings near pith anyway are often missing when working with increment cores. Part of the age effect can be absorbed by the irregular ring-width variations exhibited by trees growing in stands where thinnings induce successive episodes of growth surge.
Wood density should not be mistaken for stem-specific gravity (Williamson et al., 2010). Bark has a different mass to volume ratio than wood. The contribution of bark to the annual increment is however negligible. The calculated approximation confirms that the variations in specific gravity are proportional to that of wood density. Variations in ring width and WD at upper stem positions were however documented for different species (Bouriaud et al., 2005; Repola, 2006; van der Maaten-Theunissen and Bouriaud, 2012). These variations were mostly in the sense of a lesser reduction in growth of upper stem parts during years of limited growth. Altogether with the WD density effect, these effects show that the reaction of trees to unfavourable climate conditions is exacerbated or overestimated by the breast-height radial growth.
Annual variations in wood density were proved to compensate partially (up to 15 %) to the variations in radial growth. Ignoring the relation between ring width and wood density would result in an underestimation of the biomass production in bad years. The use of allometric equations generated estimations with large prediction intervals at tree level, up to 60 %, but the prediction errors at plot level compensated for each other. Most of the error in the prediction of a tree's annual biomass increment comes from the great between-tree variability in wood density. Plot-level errors were found to range between 10 and 20 % only. This study validates the approach based on historical diameter records for estimating tree annual biomass increment and stand annual biomass production, but a local calibration of the allometric models reduces the prediction errors considerably.
This work was supported by a grant of the Romanian National Authority for Scientific Research, CNCS-UEFISCDI, project number PN-II-ID-PCE-2011-3-0781. C. Wirth is thankful to Sebastian Weist and Ulrich Pruschitzki for help in the field and Anja Kahl for dendrochronological analysis. C. Wirth acknowledges the support of the Max-Planck Society. A. V. Kirdyanov acknowledges EU COST Action FP1106 STReESS. Ernst-Detlef Schulze and Reiner Zimmermann established the Wetzstein site. Edited by: J. Schöngart