A generalized algebraic difference approach allows an improved estimation of aboveground biomass dynamics of Cunninghamia lanceolata and Castanopsis sclerophylla forests

Abstract

Key message

A generalized algebraic difference approach (GADA) developed in this study improved the estimation of aboveground biomass dynamics of Cunninghamia lanceolata (Lamb.) Hook and Castanopsis sclerophylla (Lindl.) Schott forests. This could significantly improve the fieldwork efficiency for dynamic biomass estimation without repeated measurements.

Context

The estimation of biomass growth dynamics and stocks is a fundamental requirement for evaluating both the capability and potential of forest carbon sequestration. However, the biomass dynamics of Cunninghamia lanceolata and Castanopsis sclerophylla using the generalized algebraic difference approach (GADA) model has not been made to date.

Aims

This study aimed to quantify aboveground biomass (AGB, including stem, branch and leaf biomass) dynamics and AGB increment in C. lanceolata and C. sclerophylla forests by combining a GADA for diameter prediction with allometric biomass models.

Methods

A total of 12 plots for a C. lanceolata plantation and 11 plots for a C. sclerophylla forest were selected randomly from a 100 m × 100 m systematic grid placed over the study area. GADA model was developed based on tree ring data for each stand.

Results

GADA models performed well for diameter prediction and successfully predicted AGB dynamics for both stands. The mean AGB of the C. lanceolata stand ranged from 69.4 ± 7.7 Mg ha⁻¹ in 2010 to 102.5 ± 11.4 Mg ha⁻¹ in 2013, compared to 136.9 ± 7.0 Mg ha⁻¹ in 2010 to 154.8 ± 8.0 Mg ha⁻¹ in 2013 for C. sclerophylla. The stem was the main component of AGB stocks and production. Significantly higher production efficiency (stem production/leaf area index) and AGB increment was observed for C. lancolata compared to C. sclerophylla.

Conclusion

Dynamic GADA models could overcome the limitations posed by within-stand competition and limited biometric data, can be applied to study AGB dynamics and AGB increment, and contribute to improving our understanding of net primary production and carbon sequestration dynamics in forest ecosystems.

1 Introduction

Forest biomass plays an important role in the mitigation of climate change caused by rising atmospheric CO₂ concentrations because around 50% of tree biomass is carbon (Bonan 2008; IPCC 2007). The estimation of biomass growth dynamics and stocks has become a fundamental requirement for evaluating both the capability and potential of forest ecosystems to sequester carbon (Gower et al. 1997). It represents the net carbon input from the atmosphere to vegetation and it has received great attention during the past few decades (Fang et al. 2003). Due to the practical difficulties in belowground biomass estimation, there has been a focus on aboveground biomass (AGB), and its dynamics have become an important topic in estimating the effects of afforestation, deforestation and the role of improved forest management on the global carbon balance (Nunes et al. 2013; Zianis and Mencuccini 2005). AGB is calculated as the sum of stem, branch and leaf biomass, which can be calculated by allometric models using easily measured variables, such as tree diameter at breast height (dbh) and height (Albert et al. 2014; Cai et al. 2013). However, due to the diversity of climatic zones, site conditions and tree species, different allometric models generate great discrepancies in biomass estimates (Neumann et al. 2016). Therefore, localized- and species-specific allometric models are encouraged for biomass estimates. Besides the allometric models, biogeochemical-mechanistic approaches, such as BIOME BGC (Pietsch et al. 2005), CLM (Lawrence et al. 2011) and C-FIX (Veroustraete et al. 2002), are effective tools to estimate regional or even global biomass.

AGB increment represents a major part of net primary production (NPP) and is a key element required to understand ecosystem processes (Foster et al. 2014). It can be measured as the difference in AGB between two points in time (Foster et al. 2014). Repeated measurements of inventory plots are commonly used to estimate AGB increment or to calibrate forest growth models that can be used to predict AGB increment (Landsberg and Gower 1997; Nunes et al. 2013). However, repeated measurements require long-term planning and resources. As an alternative to long-term growth data from inventory plots and growth models, tree ring analysis by coring trees or felling a sample of trees in new inventory plots offers a retrospective approach to study forest biomass dynamics as long as past mortality, thinning or other disturbances are known (Foster et al. 2014; Liu et al. 2012; Woolley et al. 2015). Tree rings could provide reliable data on historical growth rates (Brienen and Zuidema 2006). Tree age and the long-term growth rate reflected in the tree ring width could be directly applied to estimate the development of biomass using allometric models. Tree ring analysis has successfully been used to estimate stand level AGB dynamics and carbon sequestration trends (Liu et al. 2012; Mund et al. 2002). For instance, the annual increment of woody biomass estimated by coring individual trees ranged from 1.3 to 4.5 Mg ha⁻¹ year⁻¹ (equaling 0.65–2.25 Mg C ha⁻¹ year⁻¹) in the main European forest types (Babst et al. 2014). Although the stand biomass dynamics have been successfully estimated using tree ring analysis of the standard trees (representing the mean diameter and height) in our previous published results (Tang et al. 2015), that study did not address the biomass dynamics of each individual tree, which is necessary for studying the stand development and competition. Alternatively, age-related growth models based on tree ring data from a sample of trees could be used to improve the efficiency of AGB increment estimation. Although there are many studies on age-related diameter growth models using a generalized algebraic difference approach (GADA) in different forest types (Gea-Izquierdo et al. 2008; Palahí et al. 2008; Tang 2015; Tomé et al. 2006; Wu et al. 2001), these GADA models are mainly used to predict diameter growth based on repeated measurements, with few studies focusing on AGB and AGB increment (Ogawa 2012). To our knowledge, there are so far no studies about predicting AGB increment using GADA models derived from tree ring analysis and individual-tree allometric biomass equations in Cunninghamia lanceolata (Lamb.) Hook and Castanopsis sclerophylla (Lindl.) Schott forests.

C. lanceolata is one of the most popular plantation tree species in subtropical China with relatively high timber quality and fast growth (Zhao et al. 2009). In the past few decades, the area of C. lanceolata has increased rapidly. According to the seventh national forest inventory of China, it covers an area of 8.54 million ha, which corresponds to 21% of the total plantation area in China (Jia et al. 2009). C. lanceolata plantations are usually managed in a clear cutting system with a relatively short rotation period of 25 years. Information about AGB increment with a sufficient accuracy over the rotation period is increasingly important to evaluate the carbon sequestration potential in the C. lanceolata forests and carbon market (Zhang et al. 2004), as well as the supply of biomass for bio-energy, timber industry and bio-based products.

C. sclerophylla is another important evergreen and shade-tolerant tree species. It is widely distributed across subtropical areas of China from 200 to 1000 m above sea level, and it is one of the dominant and native species in evergreen broadleaved forests (An et al. 2001). C. sclerophylla plays a significant role in water and soil conservation (Li et al. 2011). It is also an economically important species because its seeds are used to produce tofu, which is a staple in the diet of many local residents. However, during the establishment of the C. lanceolata plantations for climate mitigation and timber production, many C. sclerophylla forests have been cleared.

The objectives of this study were to (1) develop GADA models to predict diameter at breast height (dbh, 1.3 m) at any age for each individual trees in C. lanceolata and C. scleraphylla forests; (2) estimate AGB dynamics and inter-annual variability in AGB increment using GADA models and allometric biomass models; (3) compare the production efficiency of the two species. The results will provide a better understanding of AGB dynamics of these two important species in subtropical China and the application of GADA models in AGB dynamic estimation.

2 Materials and methods

2.1 Study area

The study was conducted in the Jitan township of Shitai county (29°59′-30°24′ N, 117°12′-117°59′ E), in the Southern part of Anhui province, China (Fig. 1). It is a mountainous area with a forest cover of about 80%, an elevation range of 50 to 1000 m with an average elevation of 260 m. The soils mainly originated from limestones, moorstones and mudstones (Forest Bureau of Shitai County 2004). The region has a mid-subtropical, humid, mountainous climate with distinct seasonality (Geng and Wang 2011). The recorded annual average temperature is 16 °C but can vary from −13.2 to 40.9 °C (Lu 2010). The mean precipitation is 1668 mm with high inter-annual variability, with about 70% of annual precipitation occurring from April to September (Geng and Wang 2011). The average annual sunshine duration is 1704 h and evaporative capacity is 1256 mm (Lu 2010).

Fig. 1

Location of Shitai County in the southern part of Anhui Province, South-East China

2.2 Plot design

One pure Cunninghamia lanceolata stand and one pure Castanopsis sclerophylla stand were selected for our study. Both stands had similar stand origins, originating from clear cuts of secondary forests. The area was 22 ha for C. lanceolata and 10 ha for C. scerophylla. In the first 3 years after planting C. lanceolata, understory and grass were removed manually to improve the survival rate of seedlings. Poorly formed, dead and dying trees were removed from both stands in 2004. A thinning from below with an intensity of 10% (volume removal), or 13% (trees per ha removal), was conducted in the two stands in 2010 (unpublished document from local Forest Bureau). In each stand, a systematic base grid of 100 m × 100 m was established, from which six locations were chosen randomly to establish permanent plots while an additional temporary plot was installed 25 m north of each permanent plot for destructive biomass sampling. Altogether, 12 plots were installed for C. lanceolata and 11 plots for C. sclerophylla. A forest inventory was conducted from March to May 2013, during which, dbh of all trees, height of one dominant, co-dominant and suppressed tree, tree horizontal distance, azimuth, plot slope, elevation and coordinates were measured and recorded. Tree dominance was determined based on the dbh, height and crown layer distributions. The crown of dominant trees extends above the general layer of the stand and intercepts direct sunlight across the top and upper branches. The diameter is usually amongst the largest in the stand. The crown of co-dominant trees lies within the main layer. Tree diameter lies in the middle range in the stand. The crown of suppressed trees lies entirely below the main canopy. Tree diameter is amongst the smallest in the stand (Tang et al. 2015). The plot design was a nested circular plot with different horizontal radii for different tree diameter classes (Tang 2015). An inner smaller plot with a 6-m radius was used to include trees with a dbh between 10 and 20 cm, while trees with dbh larger than 20 cm where measured up to a horizontal distance of 10 m (Tang et al. 2015).

2.3 Biomass analysis

All trees in sample plots were numbered and assigned to three dominance classes—dominant, co-dominant and suppressed trees. A sample tree representing the average dbh and height in each dominance class was destructively sampled for AGB (including stems, branches and leaves), leaf area determination and stem analysis in the temporary plots. In total, 18 C. lanceolata and 15 C. sclerophylla trees were felled, which corresponded to six and five trees for each dominance class for C. lanceolate and C. sclerophylla, respectively. The dbh and tree height of sample trees of C. lanceolata ranged from 8.2 to 27.3 cm and 8.7 to 16.6 m, respectively, and were 12.9 to 28.4 cm and 8.1 to 15.9 m for the dbh and height of C. sclerophylla. Before felling, the dbh was measured to the nearest 0.1 cm with a diameter tape and the North direction was marked. After the trees were felled, the total tree height and the heights and diameters of all branches were measured.

Live branches (and additional dead branches, if applicable) from each tree were sampled to develop allometric equations for predicting branch and leaf biomass from branch diameter, and hence total tree branch biomass and leaf biomass. From the sampled branches, a subsample of leaves was randomly selected to determine the specific leaf area (SLA), which is the ratio of leaf area and leaf biomass (Guisasola 2014). For C. sclerophylla, a random selection of approximately 20 leaves per branch was used to obtain SLA. Due to the very high number of small leaves of C. lanceolata, a subsample per branch was selected for SLA determination, as described in detail in Guisasola (2014). Leaves were scanned on an A3 scanner (Epson Expression 1000XL). The leaf area of each branch was calculated by multiplying its leaf mass with the SLA for that branch. Allometric branch mass and leaf area equations were then developed (Guisasola 2014) and used to predict the mass and leaf area of all branches, which was then summed up to the tree and the stand levels.

2.4 Stem analysis

Stems were cut into different sections at heights of 0.3, 1.3, 3 and 2-m intervals thereafter up to the treetop. Cross-sectional stem discs (about 5 cm thick) were cut at each of these heights, as well as one disc at ground level. If there was a branch at the height of the disc, a replacement disc was collected 5 cm above or below. In the end, 213 C. lanceolata discs and 119 C. sclerophylla discs were collected. To fit the GADA model for dbh predictions, 1 disc at breast height (1.3 m) of each felled tree was used, resulting in 18 C. lanceolata discs and 15 C. sclerophylla discs.

All branches, leaves and discs were dried at 100 °C to constant weight. The mean fresh to oven-dried weight ratios of the discs of a given tree were used to convert the whole stem fresh weight into dry biomass. The biomass of each individual tree was the sum of all biomass compartments. After drying, the discs were processed according to standard techniques applied in dendrochronology: the discs were sanded with progressively finer grades of sandpaper (up to 240 grits) to improve the visibility of year rings (Stokes 1996). The polished discs were rewetted and scanned at 1200 dpi resolution (Epson Expression 1000XL) and visually cross-dated by Lignovision and TSAPwin (LINTAB, RINNTECH, Inc. Freiburg) software. All the discs were cross-dated by matching wide and narrow rings to assign the same calendar year in four directions (East, West, North, and South) (Tang 2015). The annual diameter increments of each tree were calculated by averaging the cardinal radii and multiplying this value by two. The bark thickness was also measured. The tree diameter calculated from cross-dating was not equal to that measured in the field using a diameter tape due to the irregular growth rate of boles in different directions. Therefore diameters needed to be calibrated according to the procedures proposed by Liu et al. (2012). Finally, annual diameter with bark was estimated according to following the processes (Bush and Brand 2008):

$$ \mathrm{BR}=\frac{B_{\mathrm{T}}}{D_{\mathrm{OB}}}=\frac{B_{\mathrm{T}}}{\left({B}_{\mathrm{T}}+{D}_{\mathrm{IB}}\right)} $$

(1)

$$ {D}_{\mathrm{OB}}=\frac{D_{\mathrm{IB}}}{1-\mathrm{BR}} $$

(2)

Where BR is the bark ratio, B _T is the measured bark thickness from tree ring analysis, and D _IB and D _OB are the diameter inside and outside bark. The D _OB was calculated by dividing D _IB by 1-BR. This conversion was required because we developed the allometric models using D _OB as the input variable.

2.5 Development of GADA models

The GADA model was developed based on the model of Bertalanffy (1957), whose integral form can be expressed as follows:

$$ Y={a}_1{\left(1- \exp \left(-{a}_2\cdot t\right)\right)}^{a_3} $$

(3)

Where Y is the target variable; a ₁ is an asymptote or limiting value; a ₂ is a rate parameter, and a ₃ is an initial pattern parameter; t is the tree age (years). To derive a model with both polymorphism and asymptote from Eq. 3, more than one parameter has to be site-specific (Álvarez-González et al. 2010; Cieszewski and Bailey 2000). Therefore, the asymptote parameter a ₁ and the shape parameter a ₃ are assumed to be dependent on the site, and a site parameter (X) was included into the model (in Eq. 4). To achieve such a derivation, the base equation is re-parameterized into a form that is more suitable for the manipulation of these two parameters (using $\exp \left({\mathrm{a}}_1^{\prime}\right)$ instead of a ₁, and taking the natural logarithm of the function) as follows (Álvarez-González et al. 2010):

$$ Y= \exp \left({a}_1^{\hbox{'}}\right){\left(1- \exp \left(-{a}_2\cdot t\right)\right)}^{a_3} $$

(4)

$$ \ln Y={a}_1^{\hbox{'}}+{a}_3 \ln \left(1- \exp \left(-{a}_2\cdot t\right)\right) $$

(5)

These site-dependent parameters are conditioned to be consistently proportional to each other’s inverse over the site productivity by the following definitions:

$$ {a}_1^{\hbox{'}}= X\kern0.5em \mathrm{and}\kern0.5em {a}_3={b}_2+\frac{b_3}{X},\kern0.5em \mathrm{while}\kern0.5em {a}_2={b}_1 $$

Eq. 3 can be expressed as follows:

$$ \ln Y= X+{b}_2+\frac{b_3}{X} \ln \left(1- \exp \left(-{b}_1\cdot t\right)\right) $$

(6)

The solution of X involves solving roots of a quadratic equation and a selection of the suitable root to substitute into the dynamic equation. The selection process may depend on the equation parameters that in return depend on the data and the domain of the applicable ages (Cieszewski and Bailey 2000). With the initial condition value Y ₀ and t ₀, the solution of Eq. 6 is expressed as follows:

$$ {X}_0=\frac{1}{2}\Big(\left( \ln {Y}_0-{b}_2{L}_0\right)\pm \sqrt{{\left( \ln {Y}_0-{b}_2{L}_0\right)}^2-4{b}_3{L}_0\Big)} $$

(7)

where L ₀ =  ln (1 −  exp (−b ₁ t ₀)).

In this study, t ₀ is the tree age in the year 2013 derived from the stem analysis, while Y ₀ is the dbh in the year 2013.

According to the tree growth characteristics, the root most likely to be real and positive is selected. Thus, after substituting it into Eq. 6, the GADA formulation applied in our study is as follows:

$$ Y={Y}_0{\left[\frac{1- \exp \left(-{b}_1 t\right)}{1- \exp \left(-{b}_1{t}_0\right)}\right]}^{b_2+\frac{b_3}{X_0}} $$

(8)

Although this model has commonly been applied to calculate site index (Alvarez-González et al. 2010; Sharma et al. 2011) and basal area (Alvarez-González et al. 2010; Barrio-Anta et al. 2006), some authors have successfully applied this model for diameter prediction (Gea-Izquierdo et al. 2008). GADA formulations using different base models such as the log-logistic or the Hossfeld models can also be analysed (Cieszewski 2002); however, the GADA model (Eq. 8) was used because comparing the performance of different GADA models is not the focus of this study and Eq. 8 could achieve our objectives.

2.6 Model fitting

All the possible growth intervals were used to fit the GADA model; therefore, autocorrelated residuals are expected because the database contains multiple observations for each individual tree. This, however, violates the assumption of independent error terms. Therefore, we corrected the autocorrelation using a modified continuous autoregressive error structure (mCAR(x)). To account for first-order autocorrelation, a mCAR(1) model form can be used which expands the error terms according to Eq. 9 (Álvarez-González et al. 2010):

$$ {e}_{i j}={d}_i{\rho}_1^{h_{i j}-{h}_{i j-1}}{e}_{i j-1}+{\varepsilon}_{i j} $$

(9)

where e _ij is the jth ordinary residual on the ith individual (i.e. the difference between the observed and the estimated diameters of the tree i at height measurement j), d _i = 1 for j > 1 and it is zero for j = 1, ρ ₁ is the first-order continuous autoregressive parameter to be estimated, and h _ij -h _ij-1 is the distance separating the jth from the jth-1 observations within each tree, h _ij  > h _ij-1. ε _ij is an independent normal random variable with a mean value of zero. To test for the presence of autocorrelation, we used the Durbin-Watson test. The modified mCAR(1) error structure was programmed in the MODEL procedure of SAS/ETS® (SAS Institute Inc 2007), which allows for dynamic updating of the residuals.

Heteroscedasticity could lead to non-minimum variance parameter estimates and unreliable predictor intervals. To avoid this problem, the error variance was assumed to be a power function of the predicted diameter (Huang et al. 2000). The weighting factor (W _i ) was:

$$ {\mathrm{W}}_i=\frac{1}{\mathrm{Pred}.{dt}_i^k} $$

(10)

where k is a constant; Pred.dt _i is the predicted values from the fitted model. Generally, k = −2, −3/2, −1, −1/2, 1/2, 1, 3/2, 2. Since the predicted diameters are initially unknown, weighting is an iterative process. In our study, k = 3/2 showed the best results based on plots of studentized residuals against the predicted diameter.

2.7 Estimation of AGB and AGB increment

Allometric models are widely used to estimate tree biomass, either for the biomass of different biomass compartments or the tree as a whole (González-García et al. 2013; Xiang et al. 2011). The tree-level allometric models for each biomass compartment are given in Supplement Table S1. The stem, branch and leaf biomass and leaf area of all trees was calculated in the inventory plots and scaled to the reference area of 1 ha by considering an expansion factor (32 for 10-m radius plots and 88 for 6-m radius plots) for the year 2013. The total AGB of 2013 was calculated as the sum of stem, branch and leaf biomass. Following the same procedure, biomass compartment, total AGB and leaf area of 2012, 2011 and 2010 were calculated using these allometric models and the diameter of each individual tree estimated using Eq. 8. Consequently, AGB increment was calculated as the sum of the biomass increment of stems, branches and leaves of each tree per inventory plot. AGB increment of plants with dbh < 10 cm was ignored because it is assumed to contribute less than 5% to the total AGB increment for these forests and is difficult to measure in the field (Yang et al. 2010). Since a thinning from below took place in 2010 without precise data acquisition (extracted volume, tree heights and dbh), we only estimated AGB and AGB increment for the period from 2010 to 2013. We also assumed that there was no mortality after the thinning in 2010, which is consistent with the absence of any dead trees in both stands.

2.8 Data analysis

Linear regression was applied to examine the relationships between stand basal area and AGB or AGB increment. One-way ANOVA was used to test for significant differences in stand characteristics, biomass compartment, production efficiency (stem biomass production/leaf area index), AGB and AGB increment. Statistical analyses were performed using R (R Core Team 2014).

3 Results

3.1 Stand characteristics

The stand characteristics of the two studied stands are shown in Table 1. Although there was a large difference in age between C. lanceolata and C. sclerophylla stands, no significant difference was found for mean stand height, basal area or stand volume (all p > 0.05). However, leaf area index (LAI) of the C. lanceolata stand was 5.4 ± 0.6 m² m⁻², (mean ± standard error) which was significantly lower than that of the C. sclerophylla stand (6.9 ± 0.4; p = 0.049). Mean stand biomass followed the reverse trend such that the C. lanceolata stand was two times of that in the C. sclerophylla stand.

Table 1 Structural characteristics of C. lanceolata and C. sclerophylla stands in 2013 (mean ± standard error)

3.2 Performance of the GADA model

At first, the models were fitted without considering any autocorrelation. There was a strong linear trend in the residuals as a function of lag1 residuals (one lag residual from previous observations of tree diameter) for both stands (Fig. 2a–d). This indicated a strong autocorrelation among the multiple observations in each tree. Therefore, a first-order continuous time autoregressive error structure was derived. By considering this autocorrelation, although the trend in residuals as a function of lag1 residuals did not disappear completely, the trend was clearly reduced (Fig. 2e–h). The only purpose of eliminating autocorrelation was to improve the statistical interpretation of the model, and it has no practical use unless repeated measurements were conducted in the sample plot (Álvarez-González et al. 2010; Diéguez-Aranda et al. 2006). All parameter estimates were significant at the 0.001 level (Table 2). The GADA model performed very well when predicting tree diameters for both stands. We observed an adjusted R ² of 0.94 for the C. lanceolata stand and an adjust R ² of 0.96 for the C. sclerophylla stand, with a RMSE of 1.810 and 1.308 cm for both stands, respectively.

Table 2 Estimated regression coefficients and goodness-of-fit statistics of Eq. 8

3.3 AGB and AGB increment

The total AGB of C. lanceolata was 69.4 ± 7.7 Mg ha⁻¹ in 2010, increasing to 80.0 ± 8.9, 91.0 ± 10.1 and 102.5 ± 11.4 Mg ha⁻¹ in 2011, 2012 and 2013, respectively (Table 3). These values were significantly lower than those of the C. sclerophylla stand (p < 0.05). Stem biomass contributed 75%–76% to the total AGB for both stands. Branch and leaf biomass contributed 14% and 11%, respectively, to the total AGB for the C. lanceolata stand, compared to 21% and 3%, respectively, for the C. sclerophylla stand. Although the stem and branch biomass of the C. lanceolata stand was significantly lower than that of the C. sclerophylla stand, the leaf biomass of C. lanceolata was about twice as high as that of C. sclerophylla. Total AGB was strongly correlated with stand basal area for both stands (Fig. 3).

Table 3 Stem, branch, leaf and total AGB (Mg ha⁻¹) of the C. lanceolata and C. sclerophylla stands from 2010 to 2013, Mean ± standard error (percentage of total AGB in each biomass compartment , n = 12 for C. lanceolata and n = 11 for C. sclerophylla)

Fig. 3

Relationships between aboveground biomass and stand basal area for C. lanceolata and C. sclerophylla forests. For both linear regressions, a significance level of p < 0.001 was achieved

AGB increment showed a reverse trend of the AGB such that AGB increment of C. lanceolata was significantly higher (p < 0.05) than that of C. sclerophylla for the three-year period (Fig. 4). AGB increment of C. lanceolata was 10.6 ± 1.2 for 2010–2011, 11.0 ± 1.2 for 2011–2012 and 11.5 ± 1.3 Mg ha⁻¹ year⁻¹ for 2012–2013, respectively, compared to 5.8 ± 0.3, 5.9 ± 0.3 and 6.3 ± 0.4 Mg ha⁻¹ year⁻¹ for the C. sclerophylla stand. Stem biomass increment was the main contributor to AGB increment. As expected, leaf biomass increment contributed least to AGB increment with a value of 8% for C. lanceolata and 1.5% for C. sclerophylla. There was no significant inter-annual variability in compartment-wise biomass increment and total AGB increment for the two stands. Again, strong relationships between AGB increment and stand basal area were found with values of the adjusted R ² of 0.96 for C. lanceolata and 0.69 for C. sclerophylla (Fig. 5).

Fig. 4

Inter-annual biomass increment of stem, branch and leaf from 2011 to 2013 for C. lanceolata (a) and C. sclerophylla (b) forests (n = 12 for C. lanceolata and n = 11 for C. sclerophylla). The error bars mean the standard error

Fig. 5

Relationships between AGB increment stand basal area

To further compare the AGB increment of the two stands, AGB increment at the age of 16 years was estimated (Fig. 6) using the GADA model and the allometric biomass equations. AGB increments of stem, branch and leaf of C. lanceolata at age 16 years were significantly higher than those of C. sclerophylla (p < 0.01). Total AGB increment was 11.5 ± 1.3 Mg ha⁻¹ year⁻¹ for C. lanceolata and 1.5 ± 0.1 Mg ha⁻¹ year⁻¹ for C. sclerophylla.

Fig. 6

Comparison of AGB increment at age of 16 years for C. lanceolata and C. sclerophylla forests

The production efficiency was defined as the ratio of stem production and leaf area index because the carbon sequestration in stem biomass is much larger than that of other compartments (Nunes et al. 2013). As expected, stem production efficiency of C. lanceolata, ranging from 1.8 to 2.0, was significantly higher than that of the C. sclerophylla stand for the 3 years observed (p < 0.001) (about 0.8, Fig. 7), indicating that C. lanceolata is a fast-growing and efficient species. This can be confirmed by the mean stand biomass (Table 1) and AGB increment (Fig. 4).

Fig. 7

Production efficiency (stem production/leaf are index) for C. lanceolata and C. sclerophylla forests. Different letters above the error bars indicate significant differences between the tree species (p < 0.05) (n = 12 for C. lanceolata and n = 11 for C. sclerophylla). The error bars mean the standard error

4 Discussion

4.1 Development of GADA model

In our study, GADA models were developed based on data from tree ring analyses to predict dbh at any age for each individual tree. The GADA models were flexible enough to predict diameter increment at any time intervals, such as two- or five-year intervals, which allowed us to reconstruct historical diameters. Using tree-individual allometric biomass models to estimate tree compartment biomass combined with the GADA model, AGB dynamics and annual AGB increment were successfully estimated from 2010 to 2013. Such estimates of AGB increment are urgently needed to reduce uncertainties in ecosystem production and associated feedback to climate change (Frank et al. 2010). The dynamic GADA models could overcome the limitations of within-stand competition and limited biometric data and improve our understanding of net primary production in these forest ecosystems. In addition, the methodology described in the present paper is that the dbh and biomass increments can be assessed much quicker than with a permanent experiment where dbh is re-measured in several points in time. This methodology may be of interest for developing growth models in remote areas (or in other conditions) where the establishment of permanent plots is unaffordable. Another advantage of our approach is that the complete dbh growth pattern with age is reconstructed, something which may be unaffordable with permanent plots which are being re-measured regularly. The final advantage of our method is that the growing conditions of a given stand can be modelled quickly. The lag between the growing conditions and obtaining the data in permanent plots is so large that this is not a suitable method for developing models for the actual climatic conditions. To develop GADA models, there are two basic approaches used to gain the tree ring data: collecting increment cores (Babst et al. 2014; Foster et al. 2014; Groenendijk et al. 2014) or collecting complete discs after felling the trees (Brienen and Zuidema 2006; Groenendijk et al. 2014; Mbow et al. 2013). Collecting the increment cores is quicker and less destructive than sampling complete discs. However, caution is needed with tree ring analyses from increment cores because the difficulty in separating true year rings and false rings can lead to a significant bias. This is especially important for C. lanceolata where false rings are frequent and there is a high probability of incorrectly classifying false rings when using increment cores. Another problem that arises from increment cores is associated with no-pith cores, because stems are not always round when growing on a slope, due to competition and other stressors, even though there are many approaches to predict the missing rings (Brienen and Zuidema 2006; Groenendijk et al. 2014). Therefore, as for many other species (Liu et al. 2012; Maingi 2006), complete discs are preferred.

However, because the quality of the fit does not necessarily reflect the quality of the prediction, assessment of their validity is often needed to ensure that the predictions represent the most likely outcome in the real world (Yang et al. 2004). The only method that can be regarded as “true” validation involves the use of a new independent data set (Pretzsch et al. 2002; Yang et al. 2004), but the scarcity of such data forces the use of alternative approaches. The common method of splitting the data set in two portions does not provide additional information (Huang et al. 2003), and according to Myers (1990) and Hirsch (1991), the final estimation of the model parameters should come from the entire data set, because the estimates obtained with this approach will be more precise than those obtained with the model fitted from only one portion of the data. Moreover, other techniques such as double cross-validation or statistical testing provide very limited information about the predictive ability of the models (Kozak and Kozak 2003; Yang et al. 2004). Therefore, we decided to defer model validation until a new data set is available for assessing the quality of the predictions.

4.2 Leaf area index

LAI is strongly related to carbon uptake, transpiration, leaf respiration, light variability and stand growth and structure (Chen et al. 1997; Moser et al. 2007; Olivas et al. 2013), and therefore affects many biological and physical processes in plant canopies (Chen et al. 1997). The average LAI in 2013 was 5.4 ± 0.6 m² m⁻² for the C. lanceolata stand, which is close to the world average leaf area index (5.55 ± 4.14 m² m⁻²) based on a meta-analysis of 809 field plots in a literature review by Luo et al. (2002). Other authors have reported lower leaf area indices for this species, for example, 4.3 m² m⁻² in Qianyanzhou, subtropical China (Li et al. 2007). The leaf area index of the C. sclerophylla stand (6.9 m² m⁻²) was significantly higher than for the C. lanceolata stand (p = 0.049). This is consistent with Luo et al. (2002) who report that the LAI in China is commonly >6 m² m⁻² for tropical rainforest, subtropical evergreen broadleaved forest and temperate mixed forest.

4.3 Comparisons of AGB

As expected, AGB increased with age for both stands, which agrees with other studies (Chen et al. 2013; Chen 1998). The AGB of C. sclerophylla was significantly higher than that of C. lanceolata. Similar results were found for Castanopsis kawakamii forests that had a significantly higher AGB compared to nearby C. lanceolata forests (Yang et al. 2006). This was mainly due to the differences in wood density of the two species. Although the stand volume was similar (Table 1), the wood density was 0.40 kg m⁻³ for C. lanceolata and 0.64 kg m⁻³ for C. sclerophylla according to our measured data. Difference in ages was another important reason for the differences in AGB. It is generally recognized that AGB increases with stand age for young forests (Chen et al. 2013; Foster et al. 2014). In our study, the C. lanceolata stand was 17 years old, while the C. sclerophylla stand was 52 years old (Table 1).

AGB of the C. lanceolata stand was 102.5 ± 11.4 Mg ha⁻¹ in 2013 (Table 3), which was comparable to a C. lanceolata stand with similar stand age in Huitong, Hunan Province (106 Mg ha⁻¹) (Xiao et al. 2007). But it was higher than inventory-based national average biomass (80 Mg ha⁻¹) of C. lanceolata (Guo et al. 2010) and a C. lanceolata stand in Guangxi Province (86 Mg ha⁻¹) (Bing et al. 2009). Other authors observed higher values with an average of 155 Mg ha⁻¹ in a 17-year C. lanceolata stand in Tongdao County (Chen 1998) and 160 Mg ha⁻¹ in Nanping City (Chen et al. 2013). Similarly, AGB of C. sclerophylla was 154.8 ± 8.0 Mg ha⁻¹ in 2013, which lies in the range of reported biomass of the main forest types in China (48.8–215.8 Mg ha⁻¹) (Guo et al. 2010). It is higher than that of natural Castanopsis carlesii–Schima superba mixed stands in the Xiaokeng region of the Nanling Mountains, in subtropical China (115.2 Mg ha⁻¹) (Xie et al. 2013) but lower than the mean AGB in Castanopsis, Cinnamomum and Schima superba forests in China (256.4 Mg ha⁻¹) (Jiang et al. 1999). These differences in AGB are mainly related to inherent variations in growth rates, stand age, stand density and starting dbh (Houghton 2005; Jandl et al. 2007; Pan et al. 2004). For example, the AGB of C. lanceolata with a stand age of 40 years was twice as high as a 21-year-old stand (Chen et al. 2013). However, even with a similar age, site and climate conditions, the AGB of C. lanceolata in our study was lower than that reported by Chen et al. (2013) due to the high stand density of the latter study (2800–3875 trees·ha⁻¹ versus 1553 trees·ha⁻¹ in our study), indicating how stand density can be used to increase stand biomass (Jandl et al. 2007).

Stem mass was the dominant AGB compartment for every stand and contributed 76% to total AGB on average. This is consistent with C. lanceolata forests in Taiwan (80%, Yen and Lee 2011), the Huitong national research station of forest ecosystems (78%, Wang et al. 2013) and the northwest Fujian Province in subtropical China (72%, Huang et al. 2013). This value is also very similar to the mean ratio of stem biomass to AGB in forest ecosystems across China (79%, Cheng et al. 2007).

4.4 Annual AGB increment

AGB increment of the C. lanceolata stand was 10.6 ± 1.2 for 2010–2011, 11.0 ± 1.2 for 2011–2012 and 11.5 ± 1.3 Mg ha⁻¹ year⁻¹ for 2012–2013. These values were significantly higher than those of the C. sclerophylla stand (Fig. 5). This can be explained by the tree physiology, structural characteristics of individual trees and the stand age (McMurtrie et al. 1994; Nunes et al. 2013). The young age in C. lanceolata could potentially lead to a higher AGB increment because it is generally accepted that the biomass accumulation reaches its maximum relatively early in the life of a stand, and thereafter declines with forest age (Acker et al. 2002; Ryan et al. 1997).

Because there were large age difference between the two stands (17 years and 52 years), the comparison of AGB increment at the same age was more meaningful. For simplicity, AGB increment at age 16 years was estimated and AGB increment of C. lanceolata was significantly higher than that of C. sclerophylla (Fig. 6). The results indicated that the carbon sequestration rate of C. lanceolata was significantly higher than that of C. sclerophylla. It should be noted that the AGB increment of the trees thinned in 2010, and tree mortality during stand development for both stands was not included because such data was not available in this study, as mentioned above. However, thinning from below and tree mortality often occur more among small trees in these stands (Zhang et al. 2015); therefore, it is believed that these trees only account for a small part of the total AGB increment and their exclusion will not have significantly influenced the predicted AGB increment trend. Moreover, these results should be considered with caution due to the large simulation period used to estimate the AGB increment of C. sclerophylla (36 years). By considering only the dbh increment, a constant height/dbh relationship is assumed, an assumption which is reasonable for the specific conditions of this study, as all the field plots were established in a relatively homogeneous area.

In this study, the AGB increment of C. lanceolata was higher than the average net primary production of main forest types in China (3.5–6.7 Mg ha⁻¹ year⁻¹) (Fang et al. 2003). The AGB increment of the C. sclerophylla stand ranged from 5.8 ± 0.3 Mg ha⁻¹ year⁻¹ for 2010–2011 to 6.3 ± 0.4 Mg ha⁻¹ year⁻¹ for 2010–2011. It was comparable to the AGB increment of spruce, secondary birch and secondary conifer-broadleaved forests (about 6 Mg ha⁻¹ year⁻¹) at ages from 30 to 50 years in Sichuan province, subtropical China (Zhang et al. 2012), and higher than the mixed Schima superba–Castanopsis carlesii forests in subtropical China (3.9 Mg ha⁻¹ year⁻¹) (Yang et al. 2010). However, the value was lower than the national average AGB increment of the evergreen broadleaved forests in China (Fang et al. 2003).

5 Conclusions

Understanding biomass dynamics is a fundamental requirement for evaluating the carbon sequestration capabilities and potentials of forest ecosystems. In this study, AGB dynamics and AGB increment from 2010 to 2013 of C. lanceolata and C. sclerophylla were successfully estimated by using GADA models and allometic biomass models. The stand AGB of C. sclerophylla was significantly higher than that of C. lanceolata, and stem contributed 76% to the total AGB in both stand types, highlighting the importance of the stem production in carbon sequestration of forest ecosystems.

The high determination coefficients of the GADA models for both stands indicated that GADA models are a powerful tool to predict the diameter dynamics. Although repeated measurements are most accurate for increment estimation, GADA model developed using tree ring data were far more efficient and easy to operate. These dynamic models can overcome the limitations of within-stand competition and limited biometric data and can be applied to study biomass dynamics. If mortality data are available, the GADA models could contribute to the prediction of long-term biomass dynamics without repeated measurements in the field and improve our understanding of the net primary production in forest ecosystems.