Next Article in Journal
Knowledge Mapping on Nepal’s Protected Areas Using CiteSpace and VOSviewer
Next Article in Special Issue
Joining Policy Forums Together to Develop Ki-no-Eki, a Community Currency System for Forest Management in Japan: Dynamics of Policy Communication Networks
Previous Article in Journal
Geographical Data and Metadata on Land Administration in Spain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Species Composition Affects the Accuracy of Stand-Level Biomass Models in Hemiboreal Forests

Latvian State Forest Research Institute ‘Silava’, Rigas Str. 111, LV-2169 Salaspils, Latvia
*
Author to whom correspondence should be addressed.
Land 2022, 11(7), 1108; https://doi.org/10.3390/land11071108
Submission received: 28 June 2022 / Revised: 14 July 2022 / Accepted: 18 July 2022 / Published: 19 July 2022

Abstract

:
Various tree species contribute differently to total biomass stock, making the development of species-specific stand-level equations critical for better estimation of forest biomass and quantification of carbon stocks. Previously derived dry weight biomass models did not assess the effect of dominant species composition according to stand growing stock. Growing stock definitions and forest species composition differ by country, justifying the need for national stand-level biomass equations. We explored the relationship between growing stock volume and stand biomass density of above- and below-ground components in six common forest categories in Latvia using plot-level data from the National Forest Inventory from 2016 to 2020. Additionally, we explored model dependence on region, forest type, and species composition index. Models that considered growing stock and dominant species composition index performed better than models with growing stock as the only variable, especially for heterogeneous deciduous forests with greater species diversity. The elaborated models are a useful alternative to individual-level assessment for estimating forest biomass stocks in circumstances where individual tree data are not available.

1. Introduction

The increasing attention on the role of the forestry sector in reducing the amount of CO2 in the atmosphere and storing carbon in forest biomass requires continuously updated information on the condition of forests and their development [1,2,3]. Total forest biomass is an important factor in carbon balance modeling [4,5,6,7]. National Forest Inventory (NFI) data are the main source of information for international programs and statistics, such as the Forest Resource Assessment Program for the Food and Agriculture Organization (FAO) and the national greenhouse gas inventory report for the Land Use, Land-Use Change, and Forestry sector under the United Nations Framework Convention on Climate Change. The importance of NFI data in estimating forest biomass and carbon stocks is widely accepted and recognized [3,8,9].
Allometric relationships between above- and below-ground biomass components and stand characteristics are highly variable by species, stand growing stock [10,11,12], stand age [9,10,11,12,13], aridity [10], site [14,15], and other factors [16,17]. Considering this, the extrapolation of existing biomass models can lead to divergent biomass estimates, justifying the need for national and site-specific equations [18,19].
Forest stand biomass calculation can be carried out using allometric equations with individual tree or stand attributes, biomass expansion factors, and various remote sensing techniques, which are the key tools to improve forest biomass estimation on a broader scale [20,21]. Using stand attributes is more common in practical forestry since the plot data are not always available to forest owners and companies who want to assess the contribution of their forests to the total carbon stock. In ecological studies, such parameters as biomass are not easily measured directly; therefore, national biomass equations are created. Regionally, equations based on extensive and representative material for the prediction of above- and below-ground biomass components of individual trees of Scots pine, Norway spruce, birch, European aspen, common alder, and grey alder have recently been adapted in Latvia [22,23]. The obtained equations are the only ones available for larger-scale estimation of above- and below-ground biomass components across the region between boreal and temperate forests. However, stand-level biomass models are still not available.
The territory of Latvia is located in the European hemiboreal forest zone [24], where coniferous species are mixed with fast-growing broad-leaved species [25]. According to 2021 NFI data, the most widespread and commercially valuable tree species in Latvia are Scots pine (Pinus sylvestris L.), birch (mainly silver birch (Betula pendula Roth)), and Norway spruce (Picea abies [L.] Karst), which contribute 32.7%, 23.6%, and 20.0% of the total growing stock volume, respectively. Species that account for the remaining 21.3% of the growing stock are European aspen (Populus tremula L., 9.3%), common alder (Alnus glutinosa (L.) Gaertn., 6.3%), and grey alder (Alnus incana (L.) Moench., 5.8%).
In Latvia and the European hemiboreal forest zone as a whole, knowledge of factors influencing stand-level forest biomass appraisals based on inventories is limited. There is a need for reliable algorithms to obtain accurate estimates of forest biomass if individual tree data are unavailable. The objective of this study is to develop stand-level biomass models for six common species-based forest categories in Latvia’s forests. We hypothesize that derived biomass appraisals for above- and below-ground stand components are influenced by the factors of region, forest type, and species composition. Specifically, our objectives were to (I) investigate the relationship between species-based forest category biomass and stand characteristics, (II) study the dependence of the stand-level models on region, forest type, and species composition index, and (III) assess the suitability of the developed equations for national forest biomass estimates in Latvia.

2. Materials and Methods

In Latvia, information about forest resources and their dynamics has been obtained via the NFI since 2003. In the present study, data from the Latvian NFI from 2016 to 2020 were used to derive stand-level biomass models for forest categories dominated by tree species common to Latvia. The dominance of a tree species was determined as a proportion of the growing stock in the sample plot in relation to other tree species. The NFI uses a 4 × 4 km grid across Latvia with 4 permanent 500 m2 sample plots at each grid point. In total, 16,157 plots are spread throughout the territory of Latvia, and each individual plot is measured once every 5 years [26]. According to the data (NFI 2016–2020), forest land covers 3.24 million ha (52%) of Latvia, with a total growing stock of 682.3 million m3.
The stem volume of a tree in the Latvian NFI is defined as the volume of stem wood over bark above the stump, including the top of a tree. The total volume of growing stock includes the individual stem volumes of all living trees with a minimum diameter at breast height (DBH) of 2.1 cm. The calculation of growing stock is carried out using species-specific models developed by Liepa [27]. Above- and below-ground dry weight biomass components of Scots pine, Norway spruce, birch, European aspen, common alder, and grey alder were computed through species-specific allometric models [22,23]. Total tree biomass was calculated at the plot level by summing the tree biomass of the individuals within the plot. These models predict stem volume and tree component biomass using DBH and tree height (H) as explanatory variables.
If an individual NFI plot contains more than one land use category, it is divided into smaller units referred to as sectors. Because sector size varied, the growing stock and biomass were converted to a per hectare basis. For the development of stand-level biomass models, only the plots or sectors with an area of at least 400 m2 were selected, assuming they contained an adequate distribution of trees. This was done to avoid unusual observations, such as sample plots at the side of a river or field that could lead to undesirable outliers. We used NFI data for the land categories of forest land and afforested agricultural land. Ultimately, data from 6530 NFI plots met the selection criteria and were used for stand-level biomass assessment (Table 1). In Latvian forests, two tree-like birch species—silver birch (Betula pendula Roth) and downy birch (Betula pubescens Ehrh.)—are not distinguished in the forest inventory registers; therefore, we used the general term birch and did not distinguish between them.
We determined the relationships between stand age, DBH, H, stand basal area (G), growing stock (GS), and stand biomass components (total above-ground biomass (AGB), total below-ground biomass (BGB), stem biomass (SB), and branch biomass (BB) using Spearman’s correlation.
The single variable nonlinear mixed-effects function “nlmer” from the lme4 package in R statistical software [28] was fitted to the stand biomass data to analyze the sources of variation in biomass components. In this model (Equation (1), forest stand AGB, BGB, SB, and BB estimates were predicted as a function of stand variables as a fixed effect, while the region, forest type, and species composition division were accounted for separately as random effects:
Y r s c = a X r s c b 1 + u r + ε r s c ,
where Yrsc is the observed stand biomass component of the target species (t), X is the stand fixed effect variable, a and b1 are parameters, and u r is the stand random effect variable assigned to the b1 parameter. The subscripts r, s, and c refer to the random effect (region, forest type, or species composition index), stand, and biomass components, respectively.
The random effects were incorporated into this model to isolate possible biomass differences between administrative districts of Latvia and tree species since the forest type characterizes the species composition and productivity of the forest stand. The variance in random factor or standard deviation of the random variable indicates how much biomass variability is present between individual forest stands across all treatments. According to NFI data, Latvia is divided into 33 districts. Based on forest typology, the forests of Latvia are divided into dry site type forests or forests on wet mineral soils, wet peaty soils, drained mineral soils, or drained peaty soils [29]. A total of 23 forest types in Latvia are distinguished. The species composition index (CI) was calculated as a proportion of the standing volume of species of interest from the standing volume of a stand expressed in tenth parts (numeric).
Stand AGB, BGB, SB, and BB per hectare were fitted using Equation (2):
Y s c = a G S s b 1 + ε s c ,
where Ysc is the observed stand biomass component of the target species (t), GS is the stand growing stock volume of the target species (m3), and a and b1 are parameters. The subscripts s and c refer to the stand and biomass components, respectively.
The estimated random effects variances (associated with b1 parameter in Equation (1) linked with region and forest type were close to 0, indicating that the level of between-group variability was not sufficient to justify the incorporation of these variables into further models. To evaluate how the species composition index influenced the relationships between biomass components and growing stock for each species, this variable was included in Equation (3) as:
Y s c = a G S s b 1 C I s b 2 + ε s c ,
where Ysc is the observed stand biomass component of the target species (t), GS is the stand growing stock volume of the target species (m3), CI is the species composition index, and a, b1, and b2 are parameters. The subscripts s and c refer to the stand and biomass components, respectively.
We used the standard “nls” function from the stats package in R statistical software to determine the nonlinear (weighted) least squares estimates of the parameters in nonlinear Equations (2) and (3). Among the common goodness-of-fit measures, the root mean square error (RMSE), mean average percent error (MAPE), adjusted coefficient of determination (adj. R2), and residual studies were used to evaluate the accuracy of model estimates. The models were compared on the basis of Akaike’s Information Criterion (AIC). The estimated AIC value indicates the likelihood that a model is correct with a higher penalty for extra parameters in the models [30].
Finally, assessing the performance of Equations (2) and (3) at the national level, the stand-level biomass estimates were compared with the NFI reference data. The results were evaluated by comparing relative differences in relation to the NFI biomass estimates.

3. Results

Spearman’s correlation analysis revealed a significant (p < 0.001) positive correlation (trend) between AGB, BGB, SB, BB, and all examined forest stand characteristics obtained from the NFI plot-level data (Table 2). Overall, the stand growing stock, followed by the basal area, correlated best with the various stand biomass components, whereas stand age and diameter at breast height were less informative, especially for Scots pine-dominated stands. Among the biomass components, the weakest correlation (0.53–0.96) was found for stand characteristics and the BB, which is the component typically associated with the most biomass variation. Considering the relatively high correlation coefficient (0.82–0.99) with all studied biomass components, the growing stock was included as the only variable to avoid collinearity between predictor parameters in the equations.
Three random effects (region, forest type, and species composition index) were included in Equation (1) to isolate possible biomass differences. The estimated standard deviation of the random variable (associated with b1 parameter in Equation (1)) linked with region and forest type in each model was 0, indicating that there was no variability in forest stand biomass across regions and forest types in Latvia. In contrast, the effect of species composition index on the volume–biomass relationship varied according to the biomass components of the different tree species, with BB the component with the greatest influence on the standard deviation of the random variable (Table S1 Supplementary Material).
Figure 1 shows an example of how much variability exists within individual forest stand biomass across all composition index values of the dominant species modeled by Equation (1) for Scots pine forests. There is a close positive relationship between biomass components and stand growing stock for all forest categories, with mixed stands producing the greatest stand biomass for a given growing stock in Scots pine, European aspen, grey alder, and common alder forests. In birch forests, the opposite was true in that the higher the species composition index (CI) value (pure stands), the greater the stand biomass for a given growing stock.
Separate allometric equations were derived for species-based forest categories and stand-level biomass components. The parameter estimates, including the goodness-of-fit statistics of the AGB, BGB, SB, and BB for the basic model (Equation (2)), are shown in Table 3 All parameter estimates were significant (p < 0.001) and generally explained between 64% and 99% of the observed stand biomass variability for all forest categories and biomass components. There were no clear trends in model residuals (Figure S1 Supplementary Material) that could reveal systematic discrepancies between the predicted values and the data set.
The equations of stand-level dry biomass in six forest categories were tested using one- and two-variable models. In Equation (3), CI and GS were significant predictors for all studied species, improving the basic model in most cases in terms of AIC, RMSE, MAPE, and adjR2 (Table 4). The elimination of CI only improved the model fit to the data set for two derived biomass functions (the birch BGB component and common alder SB component), resulting in lower AIC scores. The difference between the two AIC values was −1 for both biomass components. In general, the RMSE and MAPE estimates based on Equation (2) were slightly larger than those based on Equation (3), indicating poorer model performance for Equation (2).
Figure S2 (Supplementary Material) illustrates examples of the variation of the AGB and BGB observed in the NFI plots and their modeled relationship with GS alone (Equation (2)) and GS and CI (Equation (3)) after plot AGB or BGB were estimated using the derived equations for the Scots pine forest category.
A comparison of Equations (2) and (3) for judging the fit of the model to observed biomass data were analyzed for each biomass component and for the total forest stand biomass, as shown in Figure 2. Statistical tests of intercept (a = 0) and slope (b = 1) did not detect a significant deviation from the ideal model (1:1 line) for all equations fitted to the dataset.
The modeling results confirmed that the tendency of the relationship between total stand dry biomass and growing stock was similar in European aspen, grey, and common alder forests and that total biomass was slightly lower than for the other studied forest categories for a given growing stock (Figure 3). In general, if the growing stock was below 300 m3 ha−1, the Norway spruce dominant forests had greater total biomass than the other forests. Compared to the other studied forest categories, the total biomass of birch-dominated forests was greater with the increase in the stand growing stock above 300 m3 ha−1.
Assessing total stand biomass using Equations (2) and (3) instead of currently used individual tree biomass equations [22,23] resulted in small differences in the estimates at the national level of 0.17% and 0.14%, respectively (Table S2 Supplementary Material). At the level of studied forest categories, the differences seem acceptable and smaller for Equation (3) (0.01–0.21%) than Equation (2) (0.09–0.32%). These results indicate that Equations (2) and (3) can be used to transform growing stock volume data into national biomass estimates by forest category.

4. Discussion

At the national level, NFI data are the most practical means of modeling forest stand biomass and estimating carbon stocks, as forest data are usually collected from all populations of interest in a statistically sound and verified manner. The elaboration of models was intended to contribute to the knowledge of the relationship between stand growing stock, a variable often measured in forestry, and the forest stand biomass of different components. Our study meets a need to improve forest stand biomass predictions by making it possible to calculate the species composition index according to stand growing stock data of different forest categories. This study was limited to NFI sample plots established in Latvia on forest land according to the FAO definition [31].
Frequently, total forest stand or above- and below-ground biomass components are estimated using easily measurable stand characteristics obtained via NFI plot-level data or field investigations [3,10,11,16]. We observed a strong and significant relationship between studied stand components (AGB, BGB, SB, and BB) and stand characteristics (stand age, DBH, H, G, and GS). Based on Spearman’s rank correlation coefficients linking forest stand variables and biomass of a particular stand component (Table 2), the biomass of every forest category was estimated more or less accurately using these stand variables in allometric equations. However, among the studied stand variables, the closest correlation, especially for stand AGB and BGB portions, was found with growing stock. It is not surprising that the growing stock volume is the most commonly used variable converted into estimates of above- and below-ground woody biomass, as the volume of growing stock is reported as one of the most important forest characteristics monitored by NFIs to quantify wood resources [8]. Methods to estimate growing stock differ by country, with growing stock definitions and calculation methods differing due to country-specific conditions, justifying the need for national stand-level equations.
The mean growing stock, according to the Latvian NFI data (from 2016 to 2020), was 210.4 m3 ha−1, varying by forest category from 123.6 m3 ha−1 in gray alder forests to 265.8 m3 ha−1 in Scots pine forests (Table 1). Similarly, the forest stand biomass density (AGB + BGB) differed among forest categories, where the estimated stand biomass in grey alder forests (78.2 t ha−1) was more than two times less than in Scots pine forests (178.9 t ha−1). The stand biomass density per unit of forest area in Latvia was approximately 147.9 t ha−1 (AGB 117.2 t ha−1 and BGB 30.7 t ha−1), which was more, on average, than in boreal forests but less than in temperate forests [6,32].
In this study, stand growing stock alone and in combination with species composition index as variables were examined in Equations (2) and (3), respectively, to calculate the stand biomass of different components. AGB, BGB, SB, and BB by forest category were better estimated using Equation (3) than the simplest model (Equation (2)). The values of all the model parameters were biologically consistent. Although the simpler model met the requirements of heteroscedastic model residuals (Figure S1 Supplementary Material) and fitted the observed data well (Figure 2), its accuracy in terms of RMSE, MAPE, and adjR2 was lower compared to the model including species composition index as a variable. Equation (2) had a higher probability of being a superior function only for the prediction of the birch BGB component and common alder SB component, resulting in lower AIC scores than Equation (3). The difference between the two AIC values was −1 for both biomass components. According to Motulsky and Arthur [30], with such a small difference between two AIC scores, we cannot be sure that the best model is more appropriate; the data are simply ambiguous.
The biomass stock of the six studied forest categories accounted for 97.8% of the total forest biomass in Latvia, indicating that the study material represents the bulk of all Latvian forests. More than half (55%) of the biomass stock was accumulated in coniferous forests, and the Scots pine forests alone, as the dominant forest category, contributed 34.1% to the total biomass stock. The majority of the forest biomass stock in Latvia can be characterized by a few tree species, in contrast to temperate forests, where there is a slightly higher tree species diversity [24,25].
The forest category and composition of the dominant species have an impact on stand biomass density, indicating that forests with similar growing stock can have different structures and component biomass. The dominant species in each forest category were the essential contributors to the accumulation of forest stand biomass (Figure 4). Compared to birch, aspen, and alder forests that were heterogeneous with mixed species biomass, coniferous forests were more homogeneous in terms of species composition. The contributions of the dominant species to total biomass for pine and spruce forests were 74.6% and 77.3%, respectively, while in the range of 54.4% to 67.4% in deciduous forests. The greater the admixture proportion of other species, the greater the differences in biomass density are expected at the same growing stock.
Numerous studies have revealed differences in stand-level biomass allocation patterns [19,33,34]. The contribution of the various biomass components (branches, stems, and roots) to the total biomass depends on the forest category and stand age, but most studies indicate that the largest part of the total tree or stand biomass corresponds to stems. From the studied species, the birch has the highest stem wood basic density, while density for coniferous species—Scots pine and Norway spruce—and deciduous species—aspen, grey alder, and common alder—is lower [35,36,37], indicating that birch stem biomass will be greater at the same stem volume. Therefore, the admixture of birch in stands dominated by other tree species enhances total biomass yield at a given growing stock. The admixture of Norway spruce in forest stands (Figure 4) also increases total biomass yield, mainly due to the biomass of the branches [19,38,39], the share of which is higher in the total biomass of spruce than in other studied species (negative b2 parameter in Table 4).
In this study, we did not aim to improve the Latvian NFI to report biomass and C stocks, but rather, we intended to present a methodology for calculating stand-level biomass components in cases where individual tree measurements are not available. Applying our models (Equations (2) and (3)), it is possible to calculate the biomass and carbon stocks [40] for forest categories according to the dominant species common in Latvia and hemiboreal forests. We found that forests with the same growing stock can have different amounts of biomass due to their dominant and admixture species. Overall, both Equations (2) and (3) are appropriate for estimating stand biomass components for all studied forest categories. However, unless information on species-specific growing stock is available, we recommend using Equation (3) with stand growing stock and composition index according to dominant species as variables. The elaborated models can be used for calculating the national scale biomass of forest categories and carbon stocks in Latvia regardless of the region and forest type.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/land11071108/s1, Table S1: the estimated allometric equations incorporating species composition index as random effect and fitting statistics for Equation (1); Figure S1: an example of how Pearson residuals are analyzed for the Equations (2) and (3) by RMSE for above- and below-ground biomass components of Scots pine forests category; Figure S2: an example of the relationship between stand level growing stock and above- and below-ground dry biomass (AGB and BGB, respectively) for the Scots pine forest category. The predicted values were estimated with the AGB and BGB functions presented in Table 3 and Table 4; Table S2: National level forest stand biomass (AGB + BGB) estimated by the NFI (reference) and Equations (2) and (3). The difference, %, is expressed as the absolute percent difference between the NFI reference value and the estimates obtained from equations presented in Table 3 and Table 4.

Author Contributions

Conceptualization, J.L. and A.L.; methodology, J.L. and S.K.; validation, J.L., S.K. and K.L.; formal analysis, K.L.; investigation, J.L. and S.K.; resources, J.L. and A.L.; data curation, J.L.; writing—original draft preparation, J.L.; writing—review and editing, K.L. and S.K.; visualization, J.L.; supervision, A.L.; project administration, S.K.; funding acquisition, A.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Regional Development Fund project, grant number 1.1.1.1/19/A/064.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

J.L.’s contribution was supported by the European Regional Development Fund, support for post-doctoral studies in Latvia “Reducing uncertainty in the calculation of forest stand biomass and carbon stock in Latvia” (No.: 1.1.1.2/VIAA/4/20/687). S.K.’s contribution was supported by the European Regional Development Fund project “Development of greenhouse gas emission factors and decision support tools for management of peatlands after peat extraction” (1.1.1.1/19/A/064). A.L.’s contribution was supported by the Joint Stock Company “Latvia’s state forests” research program “Carbon turnover in forest lands”.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Neumann, M.; Moreno, A.; Mues, V.; Härkönen, S.; Mura, M.; Bouriaud, O.; Lang, M.; Achten, W.M.J.; Thivolle-Cazat, A.; Bronisz, K.; et al. Comparison of Carbon Estimation Methods for European Forests. For. Ecol. Manag. 2016, 361, 397–420. [Google Scholar] [CrossRef]
  2. Martin, A.R.; Doraisami, M.; Thomas, S.C. Global Patterns in Wood Carbon Concentration across the World’s Trees and Forests. Nat. Geosci. 2018, 11, 915–920. [Google Scholar] [CrossRef]
  3. Lu, J.; Feng, Z.; Zhu, Y. Estimation of Forest Biomass and Carbon Storage in China Based on Forest Resources Inventory Data. Forests 2019, 10, 650. [Google Scholar] [CrossRef] [Green Version]
  4. Pan, Y.; Birdsey, R.; Fang, J.; Houghton, R.; Kauppi, P.E.; Kurz, W.; Phillips, O.L.; Shvidenko, A.; Lewis, S.L.; Canadell, J.G.; et al. A Large and Persistent Carbon Sink in the World’s Forests. Science 2011, 333, 988–993. [Google Scholar] [CrossRef] [Green Version]
  5. Ni, Y.; Eskeland, G.S.; Giske, J.; Hansen, J.P. The Global Potential for Carbon Capture and Storage from Forestry. Carbon Balance Manag. 2016, 11, 8. [Google Scholar] [CrossRef] [Green Version]
  6. Pan, Y.; Birdsey, R.A.; Phillips, O.L.; Jackson, R.B. The Structure, Distribution, and Biomass of the World’s Forests. Annu. Rev. Ecol. Evol. Syst. 2013, 44, 593–622. [Google Scholar] [CrossRef] [Green Version]
  7. Bārdule, A.; Petaja, G.; Butlers, A.; Purviņa, D.; Lazdiņš, A. Estimation of Litter Input in Hemiboreal Forests with Drained Organic Soils for Improvement of GHG Inventories. Balt. For. 2021, 29, 232–246. [Google Scholar]
  8. Gschwantner, T.; Alberdi, I.; Bauwens, S.; Bender, S.; Borota, D.; Bosela, M.; Bouriaud, O.; Breidenbach, J.; Donis, J.; Fischer, C.; et al. Growing Stock Monitoring by European National Forest Inventories: Historical Origins, Current Methods and Harmonisation. For. Ecol. Manag. 2022, 505, 119868. [Google Scholar] [CrossRef]
  9. Repo, A.; Rajala, T.; Henttonen, H.M.; Lehtonen, A.; Peltoniemi, M.; Heikkinen, J. Age-Dependence of Stand Biomass in Managed Boreal Forests Based on the Finnish National Forest Inventory Data. For. Ecol. Manag. 2021, 498, 119507. [Google Scholar] [CrossRef]
  10. Castedo-Dorado, F.; Gómez-García, E.; Diéguez-Aranda, U.; Barrio-Anta, M.; Crecente-Campo, F. Aboveground Stand-Level Biomass Estimation: A Comparison of Two Methods for Major Forest Species in Northwest Spain. Ann. For. Sci. 2012, 69, 735–746. [Google Scholar] [CrossRef] [Green Version]
  11. Di Cosmo, L.; Gasparini, P.; Tabacchi, G. A National-Scale, Stand-Level Model to Predict Total above-Ground Tree Biomass from Growing Stock Volume. For. Ecol. Manag. 2016, 361, 269–276. [Google Scholar] [CrossRef]
  12. Fang, J.-Y.; Wang, G.G.; Liu, G.-H.; Xu, S.-L. Forest Biomass of China: An Estimate Based on the Biomass-Volume Relationship. Ecol. Appl. 1998, 8, 1084. [Google Scholar]
  13. Lehtonen, A.; Mäkipää, R.; Heikkinen, J.; Sievänen, R.; Liski, J. Biomass Expansion Factors (BEFs) for Scots Pine, Norway Spruce and Birch According to Stand Age for Boreal Forests. For. Ecol. Manag. 2004, 188, 211–224. [Google Scholar] [CrossRef]
  14. Jagodziński, A.M.; Zasada, M.; Bronisz, K.; Bronisz, A.; Bijak, S. Biomass Conversion and Expansion Factors for a Chronosequence of Young Naturally Regenerated Silver Birch (Betula Pendula Roth) Stands Growing on Post-Agricultural Sites. For. Ecol. Manag. 2017, 384, 208–220. [Google Scholar] [CrossRef]
  15. Konôpka, B.; Pajtík, J.; Šebeň, V. Biomass Functions and Expansion Factors for Young Trees of European Ash and Sycamore Maple in the Inner Western Carpathians. Austrian J. For. Sci. 2015, 132, 1–26. [Google Scholar]
  16. Khan, M.N.I.; Islam, M.R.; Rahman, A.; Azad, M.S.; Mollick, A.S.; Kamruzzaman, M.; Sadath, M.N.; Feroz, S.M.; Rakkibu, M.G.; Knohl, A. Allometric Relationships of Stand Level Carbon Stocks to Basal Area, Tree Height and Wood Density of Nine Tree Species in Bangladesh. Glob. Ecol. Conserv. 2020, 22, e01025. [Google Scholar] [CrossRef]
  17. Forrester, D.I.; Tachauer, I.H.H.; Annighoefer, P.; Barbeito, I.; Pretzsch, H.; Ruiz-Peinado, R.; Stark, H.; Vacchiano, G.; Zlatanov, T.; Chakraborty, T.; et al. Generalized Biomass and Leaf Area Allometric Equations for European Tree Species Incorporating Stand Structure, Tree Age and Climate. For. Ecol. Manag. 2017, 396, 160–175. [Google Scholar] [CrossRef]
  18. Smith, A.; Granhus, A.; Astrup, R. Functions for Estimating Belowground and Whole Tree Biomass of Birch in Norway. Scand. J. For. Res. 2016, 31, 568–582. [Google Scholar] [CrossRef]
  19. Wirth, C.; Schumacher, J.; Schulze, E.D. Generic Biomass Functions for Norway Spruce in Central Europe-a Meta-Analysis Approach toward Prediction and Uncertainty Estimation. Tree Physiol. 2004, 24, 121–139. [Google Scholar] [CrossRef] [Green Version]
  20. Jalkanen, A.; Makipaa, R.; Stahl, G.; Lehtonen, A.; Petersson, H. Estimation of the Biomass Stock of Trees in Sweden: Comparison of Biomass Equations and Age-Dependent Biomass Expansion Factors. Ann. For. Sci. 2005, 62, 845–851. [Google Scholar] [CrossRef] [Green Version]
  21. Shi, L. Methods of Estimating Forest Biomass: A Review; Tumuluru, S.L.E.-J.S., Ed.; IntechOpen: Rijeka, Croatia, 2017. [Google Scholar]
  22. Liepiņš, J.; Liepiņš, K.; Lazdiņš, A. Equations for Estimating the Above- and Belowground Biomass of Grey Alder (Alnus Incana (L.) Moench.) and Common Alder (Alnus Glutinosa L.) in Latvia. Scand. J. For. Res. 2021, 36, 389–400. [Google Scholar] [CrossRef]
  23. Liepiņš, J.; Lazdiņš, A.; Liepiņš, K. Equations for Estimating Above- and Belowground Biomass of Norway Spruce, Scots Pine, Birch Spp. and European Aspen in Latvia. Scand. J. For. Res. 2018, 33, 58–70. [Google Scholar] [CrossRef]
  24. Ahti, T.; Hämet-Ahti, L.; Jalas, J. Vegetation Zones and Their Sections in Northwestern Europe. Ann. Bot. Fenn. 1968, 3, 169–211. [Google Scholar]
  25. EEA. European Forest Types. Categories and Types for Sustainable Forest Management Reporting and Policy; EEA: Copenhagen, Denmark, 2007; Volume 3. [Google Scholar]
  26. Jansons, J.; Līcīte, I. Latvia. In National Forest Inventories: Pathways for Common Reporting; Erkki, T., Thomas, G., Mark Lawrence, R.E.M., Eds.; Springer: Dordrecht, The Netherlands, 2010; pp. 341–350. [Google Scholar]
  27. Liepa, I. Pieauguma Mācība; LLU: Jelgava, Latvia, 1996. [Google Scholar]
  28. Bates, D.; Mächler, M.; Bolker, B.; Walker, S. Fitting Linear Mixed-Effects Models Using Lme4. J. Stat. Softw. 2015, 67, 51. [Google Scholar] [CrossRef]
  29. Zālītis, P.; Jansons, J. Latvian Forest Typology and Its Origin; Saule: Salaspils, Latvia, 2013. [Google Scholar]
  30. Motulsky, H.; Arthur, C. Fitting Models to Biological Data Using Linear and Nonlinear Regression. A Practical Guide to Curve Fitting; GraphPad Software: San Diego, CA, USA, 2003. [Google Scholar]
  31. FAO. Global Forest Resources Assessment 2020: Main Report; FAO: Rome, Italy, 2020. [Google Scholar]
  32. Avitabile, V.; Camia, A. An Assessment of Forest Biomass Maps in Europe Using Harmonized National Statistics and Inventory Plots. For. Ecol. Manag. 2018, 409, 489–498. [Google Scholar] [CrossRef]
  33. Pajtik, J.; Konopka, B.; Lukac, M. Individual Biomass Factors for Beech, Oak and Pine in Slovakia: A Comparative Study in Young Naturally Regenerated Stands. Trees Struct. Funct. 2011, 25, 277–288. [Google Scholar] [CrossRef]
  34. Petersson, H.; Holm, S.; Ståhl, G.; Alger, D.; Fridman, J.; Lehtonen, A.; Lundström, A.; Mäkipää, R. Individual Tree Biomass Equations or Biomass Expansion Factors for Assessment of Carbon Stock Changes in Living Biomass-A Comparative Study. For. Ecol. Manag. 2012, 270, 78–84. [Google Scholar] [CrossRef] [Green Version]
  35. Rytter, L.; Ingerslev, M.; Kilpeläinen, A.; Torssonen, P.; Lazdina, D.; Löf, M.; Madsen, P.; Muiste, P.; Stener, L.-G. Increased Forest Biomass Production in the Nordic and Baltic Countries-A Review on Current and Future Opportunities. Silva Fenn. 2016, 50, 1660. [Google Scholar] [CrossRef] [Green Version]
  36. Liepiņš, J.; Ivanovs, J.; Lazdiņš, A.; Jansons, J.; Liepiņš, K. Mapping of Basic Density within European Aspen Stems in Latvia. Silva Fenn. 2017, 51, 7798. [Google Scholar] [CrossRef] [Green Version]
  37. Liepiņš, J.; Liepiņš, K. Mean Basic Density and Its Axial Variation in Scots Pine, Norway Spruce and Birch Stems. Res. Rural Dev. 2017, 1, 21–27. [Google Scholar]
  38. Konôpka, B.; Pajtík, J.; Marušák, R. Canopy Closure Altered Biomass Allocation in Young Spruce Stand. J. For. Sci. 2015, 61, 62–71. [Google Scholar]
  39. Runģis, D.E.; Lībiete, Z.; Korica, A.; Katrevičs, J.; Jansons, Ā.; Veinberga, I.; Jansons, J. Genetic Diversity and Differentiation of Even-Aged Norway Spruce Stands in Latvia. Balt. For. 2019, 25, 45–51. [Google Scholar] [CrossRef]
  40. Bārdule, A.; Liepiņš, J.; Liepiņš, K.; Stola, J.; Butlers, A.; Lazdiņš, A. Variation in Carbon Content among the Major Tree Species in Hemiboreal Forests in Latvia. Forests 2021, 12, 1292. [Google Scholar] [CrossRef]
Figure 1. Above- and below-ground dry biomass variance, when converted to estimates, from random effects in a mixed model (Equation (1) by composition index (CI) groups and stand growing stock for the Scots pine forest category.
Figure 1. Above- and below-ground dry biomass variance, when converted to estimates, from random effects in a mixed model (Equation (1) by composition index (CI) groups and stand growing stock for the Scots pine forest category.
Land 11 01108 g001
Figure 2. The relationship between observed and predicted total (above- and below-ground components) forest stand biomass. The predicted biomass was estimated with Equations (2) and (3) presented in Table 3 and Table 4.
Figure 2. The relationship between observed and predicted total (above- and below-ground components) forest stand biomass. The predicted biomass was estimated with Equations (2) and (3) presented in Table 3 and Table 4.
Land 11 01108 g002aLand 11 01108 g002b
Figure 3. The observed and modeled (Equation (2)) relationship between stand level growing stock and total stand biomass (AGB + BGB) in Scots pine, Norway spruce, birch, European aspen, grey alder, and common alder-dominated forests.
Figure 3. The observed and modeled (Equation (2)) relationship between stand level growing stock and total stand biomass (AGB + BGB) in Scots pine, Norway spruce, birch, European aspen, grey alder, and common alder-dominated forests.
Land 11 01108 g003
Figure 4. Contribution of tree species to the total forest stand biomass for main species-based forest categories in Latvia.
Figure 4. Contribution of tree species to the total forest stand biomass for main species-based forest categories in Latvia.
Land 11 01108 g004
Table 1. Descriptive statistics of data used to develop stand-level biomass models.
Table 1. Descriptive statistics of data used to develop stand-level biomass models.
Species-Based Forest CategoryNFI PlotsGS, m3 ha−1AGB, t ha−1BGB, t ha−1SB, t ha−1BB, t ha−1
Average (Min–Max)Average (Min–Max)Average (Min–Max)Average (Min–Max)Average (Min–Max)
Scots pine1838265.8
(0.1–958.5)
143.7
(0.1–456.4)
35.2
(0.1–111.3)
112.4
(0.1–403.6)
28
(0.1–124.2)
Norway spruce1231224.4
(0.1–832.2)
128
(0.2–398.8)
35.5
(0.1–112.6)
89.9
(0.1–321.1)
34.2
(0.2–96.2)
Birch1800182.9
(0.1–776.4)
105.5
(0.1–427.5)
28.6
(0.1–114)
84.3
(0.1–352.1)
21
(0.1–91.9)
European aspen513231.1
(0.1–1015.2)
114.6
(0.1–483.9)
28
(0.1–100.7)
90.8
(0.1–404.2)
21.3
(0.1–80.8)
Grey alder582123.6
(0.1–494.3)
61.7
(0.1–248.8)
16.5
(0.1–63.3)
50.4
(0.1–215.2)
11.8
(0.1–75.4)
Common alder396210.3
(0.1–1035.6)
108.9
(0.1–494.8)
28.1
(0.1–88)
91.7
(0.1–484.6)
17.2
(0.1–53.8)
Other170173.4
(0.1–731.8)
90.9
(0.1–319.1)
30.3
(0.1–182.2)
68.1
(0.1–238.8)
32.4
(0.2–136.0)
Table 2. Spearman’s rank correlation coefficients for species-based forest category total above-ground biomass (AGB), below-ground biomass (BGB), stem biomass (SB), branch biomass (BB), and stand characteristics. All correlations were significant at a level of α = 0.001.
Table 2. Spearman’s rank correlation coefficients for species-based forest category total above-ground biomass (AGB), below-ground biomass (BGB), stem biomass (SB), branch biomass (BB), and stand characteristics. All correlations were significant at a level of α = 0.001.
Forest CategoryVariableAgeDiameter at
Breast Height
HeightBasal AreaGrowing
Stock
Scots pineAGB0.560.750.860.940.99
BGB0.550.760.860.930.99
SB0.550.750.890.930.99
BB0.530.640.620.820.82
Norway spruceAGB0.750.770.860.960.99
BGB0.750.780.850.960.99
SB0.790.800.890.940.99
BB0.620.660.700.920.90
BirchAGB0.820.830.910.950.99
BGB0.830.850.900.940.99
SB0.810.820.910.950.99
BB0.830.840.840.860.93
European aspenAGB0.930.900.940.960.99
BGB0.920.900.930.950.99
SB0.920.890.940.960.99
BB0.920.910.910.910.96
Grey alderAGB0.860.850.920.960.99
BGB0.830.810.880.970.99
SB0.860.850.920.970.99
BB0.860.850.850.880.92
Common alderAGB0.830.780.890.960.99
BGB0.790.730.840.960.98
SB0.820.770.900.960.99
BB0.770.770.740.810.85
Table 3. The estimated allometric equations and fit statistics for Equation (2).
Table 3. The estimated allometric equations and fit statistics for Equation (2).
Forest CategoryComponent *Parameter Values ± Standard ErrorsAICRMSEMAPEadjR2
ab1
Scots pineAGB1.0360.0160.8890.00312,837.58.26.00.992
BGB0.2480.0060.8930.0049473.83.28.00.981
SB0.3750.0081.0210.00313,022.58.39.00.989
BB1.6850.1110.5170.01112,837.57.929.70.703
Norway spruceAGB1.4280.0310.8400.0048791.98.69.480.991
BGB0.5530.0180.7820.0066766.53.817.40.976
SB0.2930.0061.0540.0037709.95.58.40.994
BB2.8950.1670.4800.0108505.27.626.70.845
BirchAGB0.7870.0110.9450.00211,816.16.410.30.995
BGB0.3220.0090.8710.0059537.33.419.90.977
SB0.5220.0100.9780.00312,135.77.010.20.990
BB0.5030.0380.7340.01312,283.07.335.90.802
European aspenAGB0.6440.0260.9570.0063850.210.318.60.992
BGB0.3540.0230.8210.0112984.84.421.80.955
SB0.4890.0240.9640.0083823.710.016.10.988
BB0.3960.0520.7580.0213473.37.139.80.878
Grey alderAGB0.5020.0200.9990.0073754.86.117.00.989
BGB0.3510.0250.8160.0132966.33.115.40.971
SB0.3340.0101.0370.0053161.13.616.80.994
BB0.2990.0570.7820.0353769.66.145.80.709
Common alderAGB0.7010.0240.9470.0062663.46.99.70.993
BGB0.6750.0500.7150.0132247.74.124.00.952
SB0.3220.0101.0530.0052411.75.011.10.996
BB1.0810.2180.5430.0352741.87.758.40.643
* The allometric equation of species-based forest category biomass components: B i o m a s s ( k g ) = a   G r o w i n g   s t o c k b 1 . Total above-ground biomass (AGB), below-ground biomass (BGB), stem biomass (SB), and branch biomass (BB).
Table 4. The estimated allometric equations and fit statistics for Equation (3).
Table 4. The estimated allometric equations and fit statistics for Equation (3).
Forest CategoryComponent *Parameter Values ± Standard ErrorsAICRMSEMAPEadjR2
ab1b2
Scots pineAGB1.1870.0220.8820.003−0.0480.00412,802.27.96.00.993
BGB0.3920.0090.8700.003−0.1610.0058578.12.57.00.988
SB0.3440.0091.0250.0030.0300.00512,986.78.38.90.989
BB4.3300.3230.4770.010−0.3520.01712,477.57.230.10.756
Norway spruceAGB1.3640.0350.8410.0040.0190.0068783.78.59.40.991
BGB0.4770.0190.7860.0060.0620.0106727.13.717.10.977
SB0.3560.0081.0480.003−0.0780.0057476.35.07.80.995
BB1.5830.1170.4910.0090.2680.0218352.57.225.10.863
BirchAGB0.6770.0120.9560.0020.0490.00411,640.56.110.30.995
BGB0.3140.0110.8730.0050.0090.0079537.93.419.70.977
SB0.3390.0071.0100.0030.1410.00411,147.95.38.90.994
BB1.1320.1030.6790.013−0.2760.02012,114.87.036.90.820
European aspenAGB0.7100.0250.9710.006−0.1040.0083692.38.818.20.996
BGB0.4750.0230.8480.008−0.2610.0112587.53.019.80.987
SB0.5120.0250.9710.008−0.0500.0113803.79.816.20.988
BB0.6340.0800.7800.020−0.3440.0293351.76.332.50.904
Grey alderAGB0.6930.0260.9860.006−0.1310.0073512.94.916.00.994
BGB0.6410.0420.7980.010−0.2620.0152718.72.515.40.969
SB0.3550.0121.0350.005−0.0250.0073150.13.616.60.994
BB1.1750.2230.7380.030−0.5880.0433611.85.343.20.778
Common alderAGB0.7480.0210.9720.005−0.1110.0072460.65.49.70.993
BGB0.8110.0560.7460.012−0.1940.0182146.43.616.40.965
SB0.3230.0101.0540.005−0.0070.0082412.25.011.30.996
BB2.2260.4410.6050.034−0.5800.0502639.96.739.00.725
* The allometric equation of species-based forest category biomass components: B i o m a s s ( k g ) = a   G r o w i n g   s t o c k b 1 C o m p o s i t i o n   i n d e x b 2 . Total above-ground biomass (AGB), below-ground biomass (BGB), stem biomass (SB), and branch biomass (BB).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liepiņš, J.; Lazdiņš, A.; Kalēja, S.; Liepiņš, K. Species Composition Affects the Accuracy of Stand-Level Biomass Models in Hemiboreal Forests. Land 2022, 11, 1108. https://doi.org/10.3390/land11071108

AMA Style

Liepiņš J, Lazdiņš A, Kalēja S, Liepiņš K. Species Composition Affects the Accuracy of Stand-Level Biomass Models in Hemiboreal Forests. Land. 2022; 11(7):1108. https://doi.org/10.3390/land11071108

Chicago/Turabian Style

Liepiņš, Jānis, Andis Lazdiņš, Santa Kalēja, and Kaspars Liepiņš. 2022. "Species Composition Affects the Accuracy of Stand-Level Biomass Models in Hemiboreal Forests" Land 11, no. 7: 1108. https://doi.org/10.3390/land11071108

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop