Biomass functions and carbon estimation

From AWF-Wiki
(Difference between revisions)
Jump to: navigation, search
(Constructing biomass functions)
 
(8 intermediate revisions by 3 users not shown)
Line 1: Line 1:
{{Content Tree|HEADER=Forest mensuration lecture notes|NAME=Forest mensuration}}
+
{{Ficontent}}  
__TOC__
+
The [[stem volume]] is obviously closely related to stem biomass and also to tree biomass. And the same holds for [[carbon content]]. Therefore, the basic principles of [[volume functions]] do also hold for biomass functions and for carbon functions.
 
+
 
+
==General observations==
+
+
The [[stem volume]] is obviously closely related to stem biomass and also to tree biomass. And the same holds for [[carbon content]]. Therefore, the basic principles of volume functions do also hold for biomass functions and for carbon functions.
+
 
   
 
   
 
[[Forest Definition|Forests]] play a significant role in the climate system. They are important [[carbon sinks]] and sources and the assessment of their carbon budgets has received much attention in recent years  
 
[[Forest Definition|Forests]] play a significant role in the climate system. They are important [[carbon sinks]] and sources and the assessment of their carbon budgets has received much attention in recent years  
Line 14: Line 9:
 
[[File:2.8.5-fig42.png|right|thumb|400px|'''Figure 1''' Relationship  between dbh, tree height and aboveground biomass on a metric (original)  scale (left) and after logarithmic transformation of the variables  (right) (Kleinn 2007<ref name="kleinn2007">Kleinn, C. 2007.  Lecture Notes for the Teaching Module Forest Inventory. Department of  Forest Inventory and Remote Sensing. Faculty of Forest Science and  Forest Ecology, Georg-August-Universität Göttingen. 164  S.</ref>).]]
 
[[File:2.8.5-fig42.png|right|thumb|400px|'''Figure 1''' Relationship  between dbh, tree height and aboveground biomass on a metric (original)  scale (left) and after logarithmic transformation of the variables  (right) (Kleinn 2007<ref name="kleinn2007">Kleinn, C. 2007.  Lecture Notes for the Teaching Module Forest Inventory. Department of  Forest Inventory and Remote Sensing. Faculty of Forest Science and  Forest Ecology, Georg-August-Universität Göttingen. 164  S.</ref>).]]
  
For the construction of biomass functions, a data set with observed biomass and [[Diameter at breast height|DBH]] is required, for example. However, individual tree biomass is a variable that cannot be directly observed so that also the observation of individual tree biomass is usually sample based. Destructive sampling of the different biomass compartments of an individual tree is the standard approach, where normally the main compartments leaf/needles, branches, stem and roots are treated separately for practical reasons as well as for differences in the wood density between these different functional parts.
+
For the construction of biomass functions, a data set with observed biomass and [[Diameter at breast height|DBH]] or DBH and [[tree height]] is required, for example. However, individual tree biomass is a variable that cannot be measured directly so that also the observation of individual tree biomass is usually sample based. Therefore any estimation of biomass is not affected by a measurement error alone, but carries further errors as result of modeling.
While the whole compartment is weighted for its fresh mass in the field, only samples from each compartment is taken to lab for oven drying (at about ~105°C until weight is stable). The relation between fresh and dry mass is then used as ''expansion factor'' for the whole compartment.  
+
Destructive sampling of the different biomass compartments of an individual tree is the standard approach, where normally the main compartments leaf/needles, branches, stem and roots are treated separately for practical reasons as well as for differences in the wood density between these different functional parts.
 +
While the whole compartment is weighted for its fresh mass in the field, only samples from each compartment are taken to lab for oven drying (at about ~105°C until weight is stable). The relation between fresh and dry mass is then used as ''conversion factor'' for the whole compartment.  
  
 
The dry weights of compartments can be correlated separately or summed up to the overall aboveground biomass. In cases that also the root system was sampled, the belowground part can be incorporated to calculate the total biomass. In the context of carbon budgeting on large scale or national level, the approved methodologies of the IPCC (Intergovernmental Panel on Climate Change) for the calculation of greenhouse gas balances further propose so called [[Biomass Expansion Factors]] (BEFs) for biomass quantification. BEFs are derived as relationship between [[stem volume]] and (aboveground) biomass as:
 
The dry weights of compartments can be correlated separately or summed up to the overall aboveground biomass. In cases that also the root system was sampled, the belowground part can be incorporated to calculate the total biomass. In the context of carbon budgeting on large scale or national level, the approved methodologies of the IPCC (Intergovernmental Panel on Climate Change) for the calculation of greenhouse gas balances further propose so called [[Biomass Expansion Factors]] (BEFs) for biomass quantification. BEFs are derived as relationship between [[stem volume]] and (aboveground) biomass as:
Line 41: Line 37:
  
  
{|border="1" cellspacing="0"
+
{|  
|colspan="2"|'''Table 1 Series of three conversion factors whose combination allows converting<br>stem volume to total biomass C-content (from Karjalainen and Kellomäki 1996<ref name="karjalainen_kellomäki1996">
+
|'''Table 1:''' Series of three conversion factors whose combination allows converting<br>stem volume to total biomass C-content (from Karjalainen and Kellomäki 1996<ref name="karjalainen_kellomäki1996">
Karjalainen T and S Kellomäki. 1996. Greenhouse gas inventory for land use changes and forestry in Finland based on international guidelines. Mitigation Adapt. Strategies Global Clim. 1:51-71</ref>).'''
+
Karjalainen T and S Kellomäki. 1996. Greenhouse gas inventory for land use changes and forestry in Finland based on international guidelines. Mitigation Adapt. Strategies Global Clim. 1:51-71</ref>).
|-
+
{| class="wikitable"
|colspan="2" align="center"|<math>cf=ef*x*dw*x*cc\,</math>
+
! colspan="2" |<math>cf=ef*x*dw*x*cc</math>
 
|-
 
|-
 
|<math>cf\,</math>
 
|<math>cf\,</math>
Line 58: Line 54:
 
|<math>cc\,</math>
 
|<math>cc\,</math>
 
|C-content
 
|C-content
 +
|}
 
|}
 
|}
  
 
[[File:2.8.5-fig43.png|right|thumb|400px|'''Figure 2''' Allometric functions for different tree compartments and total biomass for a Norway spruce dataset (Kleinn 2007<ref name="kleinn2007">Kleinn, C. 2007.  Lecture Notes for the Teaching Module Forest Inventory. Department of  Forest Inventory and Remote Sensing. Faculty of Forest Science and  Forest Ecology, Georg-August-Universität Göttingen. 164  S.</ref>).]]
 
 
   
 
   
 
The expansion factor approach can also be applied for tree carbon estimation, because wood carbon content is fairly constant. By combining three conversion factors, one arrives then a conversion factor that converts stem volume (which is a standard result of [[Forest inventory|forest inventories]]) into total biomass carbon content.
 
The expansion factor approach can also be applied for tree carbon estimation, because wood carbon content is fairly constant. By combining three conversion factors, one arrives then a conversion factor that converts stem volume (which is a standard result of [[Forest inventory|forest inventories]]) into total biomass carbon content.
  
The standard methodology of single-tree biomass estimation is fitting parametric regression models just as with [[volume functions]]; frequently, however, the data base for biomass functions is relatively small, simply because of the high cost of obtaining biomass data. Numerous models have been derived, most of which expressed as allometric functions of easily observable variables such as diameter in breast height (''dbh'') and [[tree height]]. Typically, these models are specific to tree species and site conditions of the particular study. Because growing conditions as well as the particular stand history are influencing the parameter estimates based on chronosequences, the derived models are restrictedly suitable for estimations on a larger scale. An early mathematical description of scaling relations in organisms was presented by Huxley (1924<ref name="huxley1924">Huxley, J. S. 1924. Constant differential growth-ratios and their significance. Nature 114; 895.</ref>, 1932<ref name="huxley1932">Huxley, J. S. 1932. Problems of Relative Growth. Methuen & Co., Ltd London.</ref>).
+
==Allometric models==
 +
The standard methodology of single-tree biomass estimation is fitting parametric regression models just as with [[volume functions]]; frequently, however, the data base for biomass functions is relatively small, simply because of the high cost of obtaining biomass data. Numerous models have been derived, most of which expressed as allometric models of easily observable variables such as diameter in breast height (''dbh'') and [[tree height]]. Typically, these models are specific to tree species and site conditions of the particular study. Because growing conditions as well as the particular stand history are influencing the parameter estimates based on chronosequences, the derived models are restrictedly suitable for estimations on a larger scale. An early mathematical description of scaling relations in organisms was presented by Huxley (1924<ref name="huxley1924">Huxley, J. S. 1924. Constant differential growth-ratios and their significance. Nature 114; 895.</ref>, 1932<ref name="huxley1932">Huxley, J. S. 1932. Problems of Relative Growth. Methuen & Co., Ltd London.</ref>).
 +
 
 +
[[File:2.8.5-fig43.png|right|thumb|400px|'''Figure 2''' Allometric  functions for different tree compartments and total biomass for a Norway  spruce dataset (Kleinn 2007<ref name="kleinn2007">Kleinn, C.  2007.  Lecture Notes for the Teaching Module Forest Inventory.  Department of  Forest Inventory and Remote Sensing. Faculty of Forest  Science and  Forest Ecology, Georg-August-Universität Göttingen. 164  S.</ref>).]]
  
 
Until today, the allometric equation <math>agb=aD^b</math> is the most commonly used regression model in empirical biomass studies, where <math>agb</math> stands for aboveground biomass, <math>a</math> for the integration factor and <math>b</math> for the scaling factor (Pretzsch, 2001<ref name="pretzsch2001">Pretzsch, H. 2001. Modellierung des Waldwachstums. Parey Buchverlag Berlin, 341 S.</ref>). The scaling factor <math>b</math> describes the proportion between two relative growth rates where <math>dM/M=b*dD/D</math>. Allometry in this sense is defined as study of the relative growth of '''a part''' of an organism in relation to the growth of the '''whole'''. For regression analysis the equation is transformed to the linear logarithmic form:
 
Until today, the allometric equation <math>agb=aD^b</math> is the most commonly used regression model in empirical biomass studies, where <math>agb</math> stands for aboveground biomass, <math>a</math> for the integration factor and <math>b</math> for the scaling factor (Pretzsch, 2001<ref name="pretzsch2001">Pretzsch, H. 2001. Modellierung des Waldwachstums. Parey Buchverlag Berlin, 341 S.</ref>). The scaling factor <math>b</math> describes the proportion between two relative growth rates where <math>dM/M=b*dD/D</math>. Allometry in this sense is defined as study of the relative growth of '''a part''' of an organism in relation to the growth of the '''whole'''. For regression analysis the equation is transformed to the linear logarithmic form:
 
  
 
:<math>agb=a*D^b\Rightarrow\mbox{ln }agb=\mbox{ln }a+b*\mbox{ln }D\,</math>
 
:<math>agb=a*D^b\Rightarrow\mbox{ln }agb=\mbox{ln }a+b*\mbox{ln }D\,</math>
 
  
 
This transformation is necessary to achieve a homogeneous distribution of variances and allows the application of a simple [[linear regression]]. As back-transformation causes estimation bias, corrections are required. Figure 1 shows the relationship between ''dbh'', tree height and aboveground biomass for a spruce dataset on a metric scale (left) and after logarithmic transformation (right).
 
This transformation is necessary to achieve a homogeneous distribution of variances and allows the application of a simple [[linear regression]]. As back-transformation causes estimation bias, corrections are required. Figure 1 shows the relationship between ''dbh'', tree height and aboveground biomass for a spruce dataset on a metric scale (left) and after logarithmic transformation (right).
Line 77: Line 73:
 
Figure 2 shows allometric functions for different tree compartments and the corresponding total biomass in a Norway spruce dataset. In most cases the coefficient of determination of biomass functions is very high (>0.9) and the consideration of tree height as additional explanatory variable is not able to enhance the model performance remarkable. This is also due to the fact that ''dbh'' and tree height are oftentimes highly correlated on the spatial limited study sites. As consequence tree height does not deliver much explanatory power.
 
Figure 2 shows allometric functions for different tree compartments and the corresponding total biomass in a Norway spruce dataset. In most cases the coefficient of determination of biomass functions is very high (>0.9) and the consideration of tree height as additional explanatory variable is not able to enhance the model performance remarkable. This is also due to the fact that ''dbh'' and tree height are oftentimes highly correlated on the spatial limited study sites. As consequence tree height does not deliver much explanatory power.
  
 
+
{|
{|border="1" cellspacing="0"
+
|colspan="6"|'''Table  2:''' Biomass regression equations for estimating biomass of tropical trees (from Brown 1996<ref name="brown1996">Brown S. 1996. Estimating Biomass and Biomass Change of Tropical Forests: a Primer. FAO Forestry Paper 134.</ref>).
|colspan="6"|'''Table  2 Biomass regression equations for estimating biomass of tropical trees (from Brown 1996<ref name="brown1996">Brown S. 1996. Estimating Biomass and Biomass Change of Tropical Forests: a Primer. FAO Forestry Paper 134.</ref>).'''
+
|-
+
|colspan="6" align="center"|<math>Y</math> = biomass per tree in kg, <math>D</math> = ''dbh'' in cm, <math>BA</math> = basal area in cm²
+
 
|-
 
|-
|align="center"|'''Equation Nr.'''
+
|colspan="6" |<math>Y</math> = biomass per tree in kg, <math>D</math> = ''dbh'' in cm, <math>BA</math> = basal area in cm²
|align="center"|'''Climatic Zone'''
+
{| class="wikitable"
|align="center"|'''Equation'''
+
!align="center"|'''Equation Nr.'''
|align="center"|'''Range in dbh'''
+
!align="center"|'''Climatic Zone'''
|align="center"|'''Number of trees'''
+
!align="center"|'''Equation'''
|align="center"|'''Adjusted r²'''
+
!align="center"|'''Range in dbh'''
 +
!align="center"|'''Number of trees'''
 +
!align="center"|'''Adjusted r²'''
 
|-  
 
|-  
 
|3.2.1
 
|3.2.1

Latest revision as of 15:44, 26 October 2013

The stem volume is obviously closely related to stem biomass and also to tree biomass. And the same holds for carbon content. Therefore, the basic principles of volume functions do also hold for biomass functions and for carbon functions.

Forests play a significant role in the climate system. They are important carbon sinks and sources and the assessment of their carbon budgets has received much attention in recent years (Apps and Price, 1996[1]; IPCC 2000, 2001[2]). Losses or gains of carbon may be due to deforestation or afforestation and these activities are explicitly included into Article 3.3 of the Kyoto Protocol as “accountable activities” in the national commitments to reduce net greenhouse gas emissions (UNFCCC, 1997). In contrast, “additional human-induced activities” related to forest management of existing forests in Article 3.4 of the Kyoto Protocol are less obvious with respect to their impact on the global carbon budget (Mund and Schulze, 2006[3]). All Parties to the UNFCCC must report periodically on forest carbon stock changes in their National Communications. There is also rising interest in establishing potentials for biomass from forests as a source of sustainable energy. Forest carbon emissions and sequestration are also relevant in establishment and management of forest plantations, silviculture, logging, forest management, deforestation, forest degradation, and desertification.

[edit] Constructing biomass functions

Figure 1 Relationship between dbh, tree height and aboveground biomass on a metric (original) scale (left) and after logarithmic transformation of the variables (right) (Kleinn 2007[4]).

For the construction of biomass functions, a data set with observed biomass and DBH or DBH and tree height is required, for example. However, individual tree biomass is a variable that cannot be measured directly so that also the observation of individual tree biomass is usually sample based. Therefore any estimation of biomass is not affected by a measurement error alone, but carries further errors as result of modeling. Destructive sampling of the different biomass compartments of an individual tree is the standard approach, where normally the main compartments leaf/needles, branches, stem and roots are treated separately for practical reasons as well as for differences in the wood density between these different functional parts. While the whole compartment is weighted for its fresh mass in the field, only samples from each compartment are taken to lab for oven drying (at about ~105°C until weight is stable). The relation between fresh and dry mass is then used as conversion factor for the whole compartment.

The dry weights of compartments can be correlated separately or summed up to the overall aboveground biomass. In cases that also the root system was sampled, the belowground part can be incorporated to calculate the total biomass. In the context of carbon budgeting on large scale or national level, the approved methodologies of the IPCC (Intergovernmental Panel on Climate Change) for the calculation of greenhouse gas balances further propose so called Biomass Expansion Factors (BEFs) for biomass quantification. BEFs are derived as relationship between stem volume and (aboveground) biomass as:

\[BEF_i=\frac{W_i}{V_{stem}}\,\],

where:

\(W_i\,\) = dry weight (Kg)
\(V\,\) = Volume (m³)
\(i\,\) = total tree or compartment.

BEFs are simple “biomass functions”. Because the relative shares of different tree compartments are changing over the tree age, BEFs are not stable over a trees lifetime. New approaches therefore try to derive age-depending BEFs. Even if criticised concerning their uncertainty and inflexibility BEFs represents the only feasible pathway to carbon assessment for many developing countries, where forest inventories and reliable, recent data on growing stock volumes and increment are rudimentary or missing.


Table 1: Series of three conversion factors whose combination allows converting
stem volume to total biomass C-content (from Karjalainen and Kellomäki 1996[5]).
\(cf=ef*x*dw*x*cc\)
\(cf\,\) Conversion factor from stem volume to total biomass C-content
\(ef\,\) Expansion factor from stem volume to total tree biomass
\(dw\,\) Conversion factor to dry matter
\(cc\,\) C-content


The expansion factor approach can also be applied for tree carbon estimation, because wood carbon content is fairly constant. By combining three conversion factors, one arrives then a conversion factor that converts stem volume (which is a standard result of forest inventories) into total biomass carbon content.

[edit] Allometric models

The standard methodology of single-tree biomass estimation is fitting parametric regression models just as with volume functions; frequently, however, the data base for biomass functions is relatively small, simply because of the high cost of obtaining biomass data. Numerous models have been derived, most of which expressed as allometric models of easily observable variables such as diameter in breast height (dbh) and tree height. Typically, these models are specific to tree species and site conditions of the particular study. Because growing conditions as well as the particular stand history are influencing the parameter estimates based on chronosequences, the derived models are restrictedly suitable for estimations on a larger scale. An early mathematical description of scaling relations in organisms was presented by Huxley (1924[6], 1932[7]).

Figure 2 Allometric functions for different tree compartments and total biomass for a Norway spruce dataset (Kleinn 2007[4]).

Until today, the allometric equation \(agb=aD^b\) is the most commonly used regression model in empirical biomass studies, where \(agb\) stands for aboveground biomass, \(a\) for the integration factor and \(b\) for the scaling factor (Pretzsch, 2001[8]). The scaling factor \(b\) describes the proportion between two relative growth rates where \(dM/M=b*dD/D\). Allometry in this sense is defined as study of the relative growth of a part of an organism in relation to the growth of the whole. For regression analysis the equation is transformed to the linear logarithmic form:

\[agb=a*D^b\Rightarrow\mbox{ln }agb=\mbox{ln }a+b*\mbox{ln }D\,\]

This transformation is necessary to achieve a homogeneous distribution of variances and allows the application of a simple linear regression. As back-transformation causes estimation bias, corrections are required. Figure 1 shows the relationship between dbh, tree height and aboveground biomass for a spruce dataset on a metric scale (left) and after logarithmic transformation (right).

Figure 2 shows allometric functions for different tree compartments and the corresponding total biomass in a Norway spruce dataset. In most cases the coefficient of determination of biomass functions is very high (>0.9) and the consideration of tree height as additional explanatory variable is not able to enhance the model performance remarkable. This is also due to the fact that dbh and tree height are oftentimes highly correlated on the spatial limited study sites. As consequence tree height does not deliver much explanatory power.

Table 2: Biomass regression equations for estimating biomass of tropical trees (from Brown 1996[9]).
\(Y\) = biomass per tree in kg, \(D\) = dbh in cm, \(BA\) = basal area in cm²
Equation Nr. Climatic Zone Equation Range in dbh Number of trees Adjusted r²
3.2.1 DRY (a) \(Y=\mbox{exp}\left\{-1.996+2.32*\mbox{ln}(D)\right\}\,\) 5 - 40 28 0.89
3.2.2 \(Y=10^\left\{-0.535+\mbox{log}_{10}(BA)\right\}\,\) 3 - 30 191 0.94
3.2.3 MOIST (b) \(Y=42.69-12.800(D)+1.242(D^2)\,\) 5 - 148 170 0.84
3.2.4 \(Y=\mbox{exp}\left\{-2.134+2.530*\mbox{ln}(D)\right\}\,\) 0.97
3.2.5 WET (c) \(Y=21.297-6.953(D)+0.749(D^2)\,\) 4 - 112 169 0.92


None of the regression equations should be used for estimating the biomass of trees whose diameter greatly exceeds the range of the original data.

  • a Equation 3.2.1 revised from Brown et al. (1989[10]) for dry forest in India, and Eq. 3.2.2 from Martinez-Yrizar et al. (1992[11]) for dry forest in Mexico (original equation based on BA). For dry zones with rainfall less than 900 mm/a use euqation 3.2.2 and for dry zones with rainfall > 900mm/a use equation 3.2.1. "exp" means "e to the power of".
  • b Both equations are based on the same data base; A.J.R. Gillespie, pers. comm. based on a revision of equation in Brown et. al. (1989[10]).
  • c From Brown and Iverson (1992[12]).

There are also tables and functions which do give predictions of the different biomass compartments like stem wood, stem bark, large branches, fine branches, leaves, roots.

[edit] References

  1. Apps, M. J. and D.T. Price (eds.) 1996. Forest Ecosystems, Forest Management and the Global Carbon Cycle, NATO ASI Series 1: Global Environmental Change. 292 p.
  2. IPCC 2001. Third Assessment Report: Climate 2001 (TAR). 304 p.
  3. Mund, M. and Schulze, E.-D. 2006. Impacts of forest management on the carbon budget of European beech (Fagus sylvatica) forests. Allgemeine Forst- und Jagdzeitung 177, 47-63.
  4. 4.0 4.1 Kleinn, C. 2007. Lecture Notes for the Teaching Module Forest Inventory. Department of Forest Inventory and Remote Sensing. Faculty of Forest Science and Forest Ecology, Georg-August-Universität Göttingen. 164 S.
  5. Karjalainen T and S Kellomäki. 1996. Greenhouse gas inventory for land use changes and forestry in Finland based on international guidelines. Mitigation Adapt. Strategies Global Clim. 1:51-71
  6. Huxley, J. S. 1924. Constant differential growth-ratios and their significance. Nature 114; 895.
  7. Huxley, J. S. 1932. Problems of Relative Growth. Methuen & Co., Ltd London.
  8. Pretzsch, H. 2001. Modellierung des Waldwachstums. Parey Buchverlag Berlin, 341 S.
  9. Brown S. 1996. Estimating Biomass and Biomass Change of Tropical Forests: a Primer. FAO Forestry Paper 134.
  10. 10.0 10.1 Brown, S., A.J.R. Gillespie and A.E. Lugo. 1989. Biomass estimation methods fortropical forests with applications to forest inventory data. Forest Science, 35:881-902.
  11. Martinez-Yrizar, A. J. Sarukhan, A. Perez-Jimenez, E. Rincon, J. M. Maass, A. Solis-Magallanes, and L. Cervantes, 1992. Above-ground phytomass of a tropical deciduous forest on the coast of Jalisco, Mexico. Journal of Tropical Ecology 8:87-96.
  12. Brown, S. and L. R. Iverson 1992. Biomass estimates for tropical forests. World Resources Review 4:366-384.

Personal tools
Namespaces

Variants
Actions
Navigation
Development
Toolbox
Print/export