Abstract

Estimation of forest carbon storage can be of great significance to the research on the productivity of terrestrial ecosystem, carbon cycle, and global warming. China has more than 54 million hm2 barren hills and waste land suitable for forestation, which provides a great potential for developing carbon sink forestry by means of forestation. This research analyzed the volume increments, volume densities, and carbon contents of 15 analytical samples of five main plantation tree species in North China, including Pinus tabulaeformis (A), Robinia pseudoacacia (B), Populus euramericana (C), Larix olgenisis (D), and Larix kaempferi (E). Results showed that carbon storage dynamic process can be expressed as follows: the ages of quantitative maturity of each tree species are 67a, 40a, 30a, 48a, 49a, respectively; the average wood densities of each tree species at different age classes are 550.93 kg/m3, 629.25 kg/m3, 404.56 kg/m3, 592.33 kg/m3, and 544.11 kg/m3,t. The average carbon contents of each tree species at different age classes are 51.48%, 46.88%, 47.81%, 46.76%, and 47.24%. It showed a significant difference between the above tree species through variance test. The maximum values of average carbon storage are 70a, 40a, 30a, 48.7a, and 49.2a, respectively. The corresponding average carbon storages are A 2.527 kg, B 3,794 kg, C 2.781 kg, D 2.996 kg, and E 3,322 kg, in a descending order of . This research, through experiment on four tree species with clear growth rings and one tree species with unclear growth rings, verified the scientific character and the scope of application of the carbon storage dynamic analysis method, providing a new method for the measurement and analysis of forest carbon storage.

1. Introduction

Large-scale devastation of natural vegetation including forests and grasslands has impaired the carbon sequestration capacity of the biosphere. Moreover the release of inactive carbon from consumption and burning of fossil fuels (coal and petroleum) has led to a considerable increase of CO2 in the atmosphere [15]. To prevent the global warming, it is legally imperative to cut down the emission of greenhouse gases [6]. The international community is making various efforts to control the emission of greenhouse gases. As forests maintain over 86% of the terrestrial carbon stock on earth, an accurate estimate of forest carbon storage will be of great significance to the research on the productivity of terrestrial ecosystem, carbon cycle, and global warming, and so forth. Researches on the forest carbon storage and carbon flux have attracted extensive attentions [6].

Over the past 20 years, plenty of researches focusing on the carbon storage, carbon density, and carbon sink function of forest ecosystem have been made by many countries [7, 8]. Using the volume-derived method, a number of scholars have estimated the changes in carbon storage of global, national, and regional forest vegetation on a time-scale basis [911]. Also, China’s scholars have estimated and calculated the carbon storage in forest vegetation planted in some of China’s provinces or regions [1216]. Although the means and methods for the research are being improved and innovated, their differences in method, time, and space scale make it difficult to apply the research findings and prediction more widely or internationally [17, 18].

As a country with scarce forest resources, China has more than 54 million hm2 barren hills and waste land suitable for forestation, which provides a great potential for developing carbon sink forestry by means of forestation. Since 2004, the State Forestry Administration has launched pilot projects of forestry carbon sinks in 6 provinces (autonomous regions) including Inner Mongolia. The rapid development of carbon sink forestry presses for a method for carbon storage measurement, and in the matter of plantation and operation of carbon sink forests, measurement and evaluation methods that can reflect the change of forest carbon storage are required so as to scientifically determine the tree species and operating cycle for carbon sink forestation.

2. Study Area

In this research, we selected 5 primary forestation tree species planted in Liaoning: Pinus tabulaeformis, Larix olgenisis, Larix kaempferi, Robinia pseudoacacia, and Populus euramericana. Pinus tabulaeformis and Robinia pseudoacacia samples were taken from Wendao Forest Farm of Fushun County (41°45′~41°53′ N, 124°1′~124°8′ E). This region features its temperate continental monsoon climate with a mean annual precipitation of 800~900 mm, a maximum temperature of 35.8°C, a minimum temperature of −34°C, a mean temperature of 7.8°C, and a frost-free period of 130–140 d. The forest land is characterized by its brown forest soil with a soil thickness of 30~55 cm and a humus horizon of 10 cm. It is a low mountainous region with an altitude of 200~300 m and a gradient of 15~25°. The most common tree species in this region are Larix olgenisis, Larix kaempferi, Pinus tabulaeformis, Pinus koraiensis, Pinus sylvestnis, Robinia pseudoacacia, and other man-made forests. In the natural forests are mainly planted with Quercus liaotungensis [19]. The Populus euramericana samples were taken from Xinmin Mechanical Forest Farm (41°48′~42°16′ N, 122°26′~122°45′ E), which is located west of Liaohe River and belongs to the sand-drift region north of Liaohe River Plain with an average altitude of 29 m. This region is low lying and flat, featured by its temperate continental monsoon climate with a mean annual temperature of 7.8°C, an extreme high temperature of 36.5°C, an extreme low temperature of −31.5°C, a mean annual precipitation of 605.7 mm, and a forest-free period of 150–160 d. The underground water level is 1~4 m deep, and most of the soil here is carbonate meadow soil formed on the alluvium parent materials out of river inundations. The main tree species here include Populus euramericana, Populus simonii, Populus pseudo-simonii, Salix matsudana and Ulmus pumila [19]. The Larix olgenisis and Larix kaempferi samples were taken from Dabiangou Forest Farm (41°59~42°04′ N, 124°59~125°10′ E) of Qingyuan Manchu Autonomous County. Belonging to Longgang branch range of Changbai Mountain with an average altitude of 500~600 m, this region is situated at a mid-low mountain land, featured by its continental monsoon climate, cold and dry in winter, hot and rainy in summer, with a mean annual temperature of 5°C, a minimum temperature of −37.7°C, a frost-free period of 120~125 d, and an annual rainfall of 700~800 mm. The forest soil here is mostly brown forest soil. Man-made forests are mainly planted with such conifers as Larix kaempferi, Larix olgenisis, Pinus koraiensis and Pinus tabulaeformis, and in the natural secondary forests are Fraxinus mandshurica, Fraxinus rhynchophylla, Quercus mongolica, and Quercus Liaotungensis [20].

3. Methods

3.1. Theory

The carbon storage dynamic analysis is a method to explore the change in carbon storage speed of forest trees at different age classes during growth. This method bases on the close correlation between growing stock (stem volume) and tree biomass [12] and the widely verified Theoretical Volume Equation with high accuracy and fitness that is of universal significance to the measurement of the stem volume of conifers [21]. Since the annual carbon storage in forests is a variable attached to the tree volume growth percentage [22], thus we can just devote ourselves to the measurement of carbon storage in stem biomass in place of the measurement of stem volume increment and subdivide carbon storage into such indexes as volume increment, wood density, and carbon content. This method aims to explore the carbon storage in forest trees and its change law at different stages (each age gradation) [20] through subdividing the carbon storage process into such indexes as the volume increment, wood density and carbon content of forest trees at different stages, measuring directly with instruments, and then a chain of calculations and analyses of these data. For man-made forests that feature their straight trunks and less branches, we can measure and evaluate the result of carbon storage in trunks of different tree species during growth in a direct-viewing and handy way and further calculate the tree biomass by use of such a characteristic as high positive correlation between growing stock and tree biomass, calculate the tree biomass based on the law of change in carbon storage on the time series, and ultimately realize the measurement and evaluation of carbon storage in forest trees.

3.2. Experimental Design and Measuring Method

Plots with the highest stand age and moderate conditions for tree growth were selected from five of main tree species planting areas as the standard woodlands, including one Pinus tabulaeformis and Robinia pseudoacacia sample plot at Wendao Forest Farm of Fushun County, one Populus euramericana sample plot at Xinmin Mechanical Forest Farm, one Larix olgenisis sample plot and one Larix kaempferi sample plot at Dabiangou Forest Farm of Qingyuan Manchu Autonomous county, each having a size of 20 m 30 m. The experiment was carried out by selecting as samples Robinia pseudoacacia, Pinus tabulaeformis, Larix olgenisis, and Larix kaempferi that have clear growth rings and Populus euramevicana that has unclear growth rings, in order to verify the feasibility of this experimental method when it is used for the tree species with different growth ring structures. The stand ages of Pinus tabulaeformis, Robinia pseudoacacia, Populus euramevicana, Larix olgenisis, and Larix kaempferi were 85a, 45a, 28a, 45a, and 45a respectively; their average diameters at breast height were 24.6 cm, 26.5 cm, 24.8 cm, 23.4 cm, and 25.8 cm respectively; their average heights were 18.3 m, 20.4 m, 24.8 m, 25.2 m, and 26.1 m respectively; their stand densities were 672 N/hm2, 716 N/hm2, 640 N/hm2, 655 N/hm2, and 623 N/hm2, respectively. Single tree measurement was made in each sample plot to calculate the average diameter at breast height, and three standard trees were selected according to their average diameter at breast height, from each sample plot for a stem analysis. Belt transects were taken from each disc to measure the wood density, and block samples were taken on an age gradation basis from each belt transect to measure the carbon content. All the data obtained through measurements and analyses were the average values of the three analytic tree samples, as shown in Table 1.

3.2.1. Growth Increment Measurement

The stem analysis method, a method to explore the growth process of forest trees, was adopted to measure the growth increment. Cut each trunk into sections, and on each cross section measure the diameter increments at every age gradation based on the growth ring widths, and on the profile calculate the tree height increment at different ages based on the cross-section height and the difference of growth rings between adjacent profiles, thus the volumes and form factors [23] can be calculated. Cut down the sample trees at a height of ground diameter, and cut the trunks into 2~4 cm thick discs at intervals of 0 m, 1 m, 1.3 m, 3 m, 5 m, and so forth for a diameter measurement. The age gradation of Pinus tabulaeformis, Robinia pseudoacacia, Larix olgenisis and Larix kaempferi was 5a, and that of Populus euramevicana was 2a. Use LINTAB5 Growth ring Measuring Apparatus to measure the growth increment of each disc at different age classes along eastwest and northsouth diameter lines (to the accuracy of 0.01 mm), take their average value as the age gradation diameter value, then calculate the current annual increment and periodic average increment of stem volume at different age classes, to express the growth speed of forest trees at a certain moment (one year) and their average growth speed at a certain stage (stand age).

3.2.2. Wood Density Measurement

Wood density means the mass of wood per unit volume. Wood density is related with growth ring width, and the latter is related directly with carbon storage, so the wood ring density plays an important role in analyzing the change of carbon storage [24]. In general, there are four types of wood density [25] according to the difference of water content in wood, and the basic density of wood (hereinafter called wood density) was adopted in this paper. Across the pith, cut the discs into 1 cm wide and high belt transects one by one in the south-north direction, and then cut the belt transects into block samples on an age gradation length basis, which should be sealed with market-sold plastic membrane and marked with numbers for later use. Put all the block samples into DHG-9053A electrothermal blowing dry box for drying. Select 2~3 clock samples for a trial weighing for the first time, carry out the second trail weighing in 2 hours, and repeat the experiment several times. At last, carry out the third trial weighing. When the difference of measured value between the third and the second weighing is within an order of magnitude of , the samples will be judged to be fully dry: where is the basic density of wood, kgm−3; is the mass of fully dried wood, kg; is the volume of green wood.

3.2.3. Carbon Content Measurement

Sample from the numbered block samples with different age classes, of which Pinus tabulaeformis and Populus euramevicana were sampled mechanically from every two age classes, but for other tree species, from every age gradation. Put the carbon content samples into FZ102 mini plant disintegrator for crushing (to a diameter of 100 meshes), and take 10 mg of the crushed samples to pack up and mark with number, and at last, put the samples into a Vario ELIII elemental analyzer (made by Germany’s Elementar Company) to measure the carbon contents. Burn the samples at 1150°C in a combustion chamber by the dynamic combustion method, with the help of an oxygen-adding device so that the samples could oxidize and decompose in the flow of high-concentration oxygen and helium. Since the samples were crushed as fine as 100 meshes and only a little amount was taken, the samples could be burnt completely. The combustion product was reduced to N2 and SO2 in a reduction tube that was provided with a copper wire and heated up to 850°C. SO2 and H2O in the mixed gas were removed one by one, and CO2 was absorbed when passing through a U-shape tube. Heat the U-shaped tube to release CO2, which simultaneously was detected by TCD, thus C content could be calculated [26].

3.3. Data Calculation

The data sorting and analysis were realized under Excel 2003 through variance analysis and multiple comparison and calculation by use of SPSS17 software. Corresponding wood density, carbon content and volume increment were used to calculate the carbon storage at different stages. Such two indexes as current annual carbon storage and average carbon storage were adopted to calculate and analyze the carbon storage dynamic change [12, 21], calculated through the following equations: where is the current annual carbon storage, kg; is the average carbon storage, kg; is the carbon content, %; is the wood density, kg.m−3; is the current annual increment of volume, m3; : average increment of volume, m3.

4. Results

4.1. Analyses of the Growth Process of Stem Volume

The diameter data of each age classes was input into the analysis software called Forstat 2.0 Wood of Statistics. By using of the stem analysis function of the statistic software, the average increments and current annual increments of volume of all the experimental tree species were calculated, as shown in Table 2.

The current annual increment of Pinus tabulaeformis keeps increasing in a uniform speed basically before 40a, but speeds up after 40a and reaches the peak at 55~60a, then drops down gradually. The average increment of Pinus tabulaeformis increases slowly before 25a, but speeds up gradually after that, and tends to be steady with the increase in age. The curves of the current annual increment and the average increment intersect at 65a. The current annual increment of Robinia pseudoacacia begins increasing rapidly from 10a, slows down gradually from 25a, and begins dropping down year by year since reaching the peak at 35a. Its average increment keeps increasing in a uniform speed after 10a and reaches the peak at 45a. The curves of the current annual increment and the average increment intersect at 40a. The current annual increment of Populus euramevicana increases rapidly before 10a, then begins slowing down after that until 20a when its growth speeds up, reaches the peak at 24a, and drops down gradually after that. Its average increment increases rapidly before 10a, and slows down after 10a until 20a when its average increment continues to increase. No intersect between the current annual increment and the average increment is found yet, but it is estimated through the tendency chart that they will intersect at 30a. The current annual increment of Larix olgensis increases slowly before 10a, begins speeding up after 10a, drops down a little at 30a, reaches the peak at 35a, and slows down gradually at 35a. Its average increment increases slowly before 10a and increases significantly after 10a until 45a. No intersect between the current annual increment and the average increment is found yet, but it is estimated through the tendency chart that they will intersect at 49a. The current annual increment of Larix kaempferi increases slowly before 10a, speeds up gradually after that until 25a when an inflection point appears, as is just like Larix olgenisis, and reaches the peak at 25a. Its average increment increases rapidly before 20a and tends to be steady after 30a. It can be estimated that the curves of the current annual increment and average increment of this tree species tend to intersect.

4.2. Analysis of Wood Density of Trunks at Different Stages

Measurement and calculation of wood density were performed using of formula (1). The wood density of the five tree species varies regularly with the increase in age (Table 2). The wood density of Pinus tabulaeformis and Robinia pseudoacacia tends to increase slowly all the way since they begin growing until 35a when reaching the peak. After 35a, the wood density of Pinus tabulaeformis drops down slowly, but that of Robinia pseudoacacia drops down rapidly; that of Populus euramevicana increases in a quick speed before 4a, but slowly at 4a~18a and reaches the peak at 18a, and tends to drop down after 18a till 28a; that of Larix olgenisis and Larix kaempferi tend to increase since they begin growing, reach the peak at 30a, and drop down gradually and tend to be steady after 30a. This demonstrates that the wood density of the same tree species fluctuates greatly at different growth stages.

The average values (standard deviation) of wood density of the above tree species are as follows: Robinia pseudoacacia 629.25 kg/m3 (43.86), Larix olgenisis 592.33 kg/m3 (73.41), Pinus tabulaeformis 550.93 kg/m3 (42.70), Larix kaempferi 544.11 kg/m3 (56.23), and Populus euramevicana 404.56 kg/m3 (14.84), showing an extremely remarkable difference between the five tree species through variance test. Through multiple comparison by LSD method, the difference is extremely remarkable between Pinus tabulaeformis, Robinia pseudoacacia, and Populus euramevicana, between Robinia pseudoacacia, Larix kaempferi, and Populus euramevicana, and between Populus euramevicana, Larix olgenisis, and Larix kaempferi; the difference is remarkable between Larix olgenisis, Pinus tabulaeformis , and Larix kaempferi . The wood density analysis results demonstrate that wood densities between different tree species are great.

4.3. Analysis of Carbon Content in Trunks at Different Stages

Carbon contents of these five tree species at different age classes were measured as shown in Table 4. According to the table, Pinus tabulaeformis has higher carbon content at 10a and 20a, and reaches the peak at 20a up to 60.52%, minimum 46.14% at 40a. In general, the carbon content of Pinus tabulaeformis fluctuates greatly during its whole growth process, but tends to be steady after 40a; the carbon content of Robinia pseudoacacia is high at 5a~20a, reaches the peak at 5a, up to 56.96%, and drops down rapidly after 25a, minimum 42.58% at 25a, having a great fluctuation during the whole growth process, but tends to be steady after 25a; changing slightly during its whole growth process, the carbon content of Populus euramevicana increases slowly at 4a~16a and reaches the peak at 16a, up to 48.64%, minimum 46.64% at 28a; the maximum carbon content of Larix olgenisis is 48.76%, the minimum is 48.00%, having a great fluctuation. Its carbon content is very high at 5a~10a, then drops down gradually, to the minimum at 30a, but increases again gradually with the increase in age, then keeps steady basically with the growth of forest trees; the maximum carbon content of Larix kaempferi is 48.00%, the minimum 46.80%, having a slight fluctuation. Its carbon content is high at 5a, slows down after that, and drops down to the minimum at 30a, then increases slowly after 30a. Larix olgenisis and Larix kaempferi have a more consistent carbon content than the other tree species.

On the basis of analysis, the average carbon contents of Pinus tabulaeformis, Robinia pseudoacacia, Larix olgenisis, Larix kaempferi, and Populus euramevicana were worked out to be 51.48% (5.94%), 46.88% (4.41%), 46.76% (0.84%), 47.24% (0.43%), and 47.81% (0.82%), respectively. Through variance test, the carbon contents between different tree species are remarkably different . The result of multiple comparisons by LSD method is that there is an extremely remarkable difference between Pinus tabulaeformis, Robinia pseudoacacia and Larix olgenisis; there is a remarkable difference between Pinus tabulaeformis, Larix kaempferi, and Populus euramevicana. The carbon content research results demonstrate that there is a great difference in the carbon content of the five tree species, and the carbon content fluctuates at different growth stages. The research also demonstrates that carbon storage dynamic analysis method can be used to measure the true carbon content of different tree species and the same tree species at different growth stages in a more accurate way.

4.4. Calculation of Carbon Storage in Trunks at Different Stages

The carbon storage was worked out through analyzing the increment, density, and carbon content of stem volume; therefore, the differences in increment, wood density and carbon content of different tree species would surely lead to the difference in carbon storage. The carbon storage can be divided into current annual carbon storage and average carbon storage (Figure 1). The current annual carbon storage and average carbon storage of Pinus tabulaeformis, Robinia pseudoacacia and Populus euramevicana show a change trend basically alike at different age classes—increasing with age and slowing down when reaching the peak. The current annual carbon storage of Pinus tabulaeformis reaches the peak at 60a, and its average carbon storage reaches the peak at 70a; the current annual carbon storage of Robinia pseudoacacia reaches the peak at 35a, and its average carbon storage reaches the peak at 40a; the current annual carbon storage of Populus euramevicana reaches the peak at 24a, but no maximum of its average carbon storage appears yet. The current carbon storage and average carbon storage of Larix olgenisis intersect at 48.7a; the current carbon storage and average carbon storage of Larix kaempferi intersect at 49.2a.

According to the forest maturity theory, the carbon storage maturity age refers to the time when the average carbon storage of a tree reaches the peak [19, 22]. Based on the method to determine the maturity age, we worked out the carbon storage maturity ages of the previous tree species as follows (average carbon storage): Pinus tabulaeformis 70a (2.5267 kg); Robinia pseudoacacia 40a (3.7939 kg); Populus euramevicana 30a (2.7810 kg); Larix olgenisis 48.7a (2.996 kg); Larix kaempferi 49.2a (3.322 kg), in a descending order of Robinia pseudoacacia Larix kaempferi Larix olgenisis) Populus euramevicana Pinus tabulaeformis.

In order to forecast and evaluate the average carbon storage of the above tree species at different growth stages, a regression analysis was made for the calculation results of their carbon storage, by use of the following regression formula:

Pinus tabulaeformis:

Robinia pseudoacacia:

Populus euramevicana:

Larix olgenisis:

Larix kaempferi: where “” represents the carbon contents of the above tree species, and “” represents the stand age.

5. Conclusion and Discussion

Scientific verification of the carbon storage dynamic analysis method: the carbon storage dynamic analysis method realizes the scientific measurement on the long-time storage of carbon in stem volume of the same tree species. The dynamic data about carbon storage during the whole growth process obtained through this method can evaluate the true change in carbon storage in a more accurate way than the carbon content coefficient method that is commonly used. This method provides a theoretical basis for how to select plantation tree species for small-scale sites and how to determine the operation cycle. Favorable verifications were obtained through experiment on four types of broadleaf trees and conifers with clear growth rings and a type of broadleaf tree with unclear growth rings; simultaneously, the scope of application of this method was verified, that is, this method is applied to not only the tree species with clear growth rings, but those with unclear growth rings.

Analyses of the volume increment of different tree species: the increments of the five tree species are greatly different on account of such factors as tree species characteristics and site conditions, and so forth, but change regularly in the form of quadratic or cubic equation curve at different growth stages. The mature ages of Pinus tabulaeformis, Robinia pseudoacacia, Robinia pseudoacacia, Larix olgenisis, and Larix kaempferi obtained from this experiment are 67a, 40a, 30a, 48a, and 49a, respectively, consistent with most research results achieved before. This demonstrates that the selected experimental plots are typical ones, which make the subsequent wood density and carbon content analysis results more reliable.

Analyses of wood density of the above tree species: there is a remarkable difference between the above five tree species in wood density, which shows a change tendency similar to a conic. Pinus tabulaeformis shows a tendency to increase slowly at its initial growth stage, reaches the peak at 35a, drops down gradually after that, and tends to be steady; the wood density of Robinia pseudoacacia increases slowly before 30a, but drops down in a quick speed after 30a; the wood density of Populus euramevicana increases in a quick speed before 4a, but slowly at 4a~18a, and tends to decrease after 18a; the wood density of Larix olgenisis and Larix kaempferi, which belong to the same tree species, changes in a very similar way, but the former is far higher than the later.

Analyses of carbon content of the above tree species: among the above five experimental tree species, Pinus tabulaeformis has the highest carbon content, over 50%, but for the other tree species, lower than 50%. According to the variance analysis results, the difference in carbon content of Pinus tabulaeformis, Robinia pseudoacacia, and Larix olgenisis at different stages is extremely remarkable; for Populus euramevicana and Larix kaempferi, the difference is remarkable. The carbon contents of these tree species are similar to precious researches [27]. The dynamic analysis results demonstrate that the carbon content of Pinus tabulaeformis reaches 60.52% at 15a~20a, drops down to the minimum 46.14% at 40a, with a severe fluctuation, but tends to be steady after 40a; the carbon content maximum of Robinia pseudoacacia appears at 20a, up to 56.96%, drops down to the minimum 42.58% at 25a, with a severe fluctuation, but tends to be steady after 25a; the carbon content of Populus euramevicana increases slowly all the way before 16a, but drops down slowly after 16a, with a maximum of 48.64% and minimum of 46.64%; the change in carbon contents of Larix olgenisis and Larix kaempferi is steady and similar, but Larix kaempferi has a higher carbon content than Larix olgenisi, though the difference is not remarkable. This is just opposite to the wood density of the above tree species.

Carbon storage maturity of the above tree species: carbon storage maturity is a term put forward by reference to the forest quantitative maturity theory. According to the quantitative maturity theory, the average carbon storage reaches the peak when the curves of current annual carbon storage and average carbon storage intersect, that is, the carbon storage maturity age is reached, so felling trees at such a time can realize an optimal utilization of land [19, 24, 28]. This research demonstrates that the average carbon storages of Pinus tabulaeformis, Robinia pseudoacacia, Populus euramevicana, Larix olgenisis, and Larix kaempferi are 2.527 kg, 3.794 kg, 2.781 kg, 2.996 kg, and 3.322 kg, appearing at their 70a, 40a, 30a, 48.7a, and 49.2a, respectively, in a descending order of Robinia pseudoacacia >  Larix kaempferi>  Larix olgenisis>  Populus euramevicana>  Pinus tabulaeformis. The carbon storage regression formulas for the above tree species obtained from this paper may also be applied to the prediction and evaluation of carbon storage.

Knowledge on the carbon storage in trunks that was obtained through this research: first, the carbon storage of the above tree species calculated by this method relies on the comprehensive effect of such three factors as carbon content, wood basic density, and the average increment of volume. This experiment demonstrates that some tree species, for example, Pinus tabulaeformis, have lower carbon content, though their carbon content is very high. Also, the maximum values of average carbon storage of different tree species appear at different age classes, which should be given full consideration in the process of forestation and operation, especially in the carbon sink function management of the existing forests; in addition, though the carbon contents of the five tree species are similar to those of constructive species of main forest types in North China [27], which demonstrates that there is a slight difference in the carbon content for the same tree species at different regions, however, whether there is regional difference in carbon storage for the same tree species? What is the effect of various operation means and site conditions on the carbon storage of each tree species? There is lack of researches yet on these aspects nowadays, for which we must make more efforts.

Acknowledgment

This study was funded by the National Natural Science Fund Projects of China, no. 30972364.