Patterns of Biomass and Carbon Allocation across Chronosequence of Chir Pine ( Pinus roxburghii ) Forest in Pakistan : Inventory-Based Estimate

*e quantitative relationship between carbon sequestration potential and stand ages of Pine (Pinus roxburghii) forest is not documented in Pakistan. Using field inventory data, this study underlines the patterns of biomass and carbon allocation across a chronosequence of Chir Pine forest. Based on the uniform shelterwood silvicultural management system, the forest was classified into three stand age classes representing the young stand (<50 years), mature stand (50–75 years), and overmature stand (> 75 years). *e results showed an increasing trend in living tree biomass carbon with stand age. However, soil carbon showed gradually decreasing trend from young to overmature stand. Similarly, deadwood, litter, and understory biomass carbon showed an increase pattern of changes. Altogether, the results highlighted that the mean carbon values of all components varied between 90.3 t·C·ha in the young stand and 309.5 t·C·ha in the overmature stand. Furthermore, our results confirm that the current management operations affect the forest floor and soil carbon. *erefore, we suggest that different protection measures should be considered during management operations to enhance soil and forest floor carbon.


Introduction
Biomass and carbon storage in any forest ecosystem play a multiple and significant role in the global carbon cycle [1].Forest ecosystems have more potential of storing carbon than another terrestrial ecosystem [2].In a forest, carbon is stored in different pools, but tree and soil are the main components that store more carbon than other components [2][3][4][5].In order to enhance the global carbon sequestration and to increase carbon-absorbing capacity, reforestation, afforestation, and the preservation of existing forest have been suggested as an effective way for mitigation of elevated CO 2 concentrations [6][7][8].e contribution of forest carbon storage in total carbon of terrestrial ecosystem is significantly important.Generally, forest stand development has a strong relationship with carbon pools of forest ecosystem because tree growth rate of forest ecosystem mostly varies with stand age [8,9].Trees increase their biomass as well as enhance growth by absorbing CO 2 from the atmosphere and sequester it in different tree components [10].As biomass increases, the growth of tree species also increases [11].Carbon sequestration of a forest depends on forest type, forest stand age, and dominant tree species [11].Previous studies showed that forest stand age may have a significant effect on the carbon stock among different ecosystem components [12].
Stand age is an important factor affecting storage of carbon in a forest in various components such as in a tree, understory, deadwood, and soil [6,13,14] and a strong predictor of the structure and function of a forest ecosystem which may affect the carbon density among carbon pools [15].Furthermore, various carbon pools of a forest respond differently to stand age [16].erefore, it is important to understand the relationship of stand age with different carbon components and sequestration rate [17].A number of studies underline the effect of stand age on carbon allocation among different carbon banks [6,7,18].e biomass carbon of living trees increases with an increase in stand age [19,20].However, this increasing trend followed a sigmoidal trend, as young trees carbon grow rapidly up to a certain age and then gradually decline [21].
In Pakistan, the researchers mainly focused on the measurement of forest carbon in the subtropical forest [22], temperate forest [23][24][25], and in the planted forest [26].However, no up-to-date studies have been carried out on the carbon allocation with respect to stand age.In addition, biomass carbon measurement is based on already available inventory data in which the minimum diameter is selected from 10 cm to 16 cm in various forests types of Pakistan, which may underestimate the available carbon storage and sequestration rate [22].In order to assess the carbon allocation among different carbon pools in relation to stand age, the present study was conducted in subtropical Chir forest.Pinus roxburghii is the dominant species of subtropical Chir zone in Pakistan.It is widely distributed in the Murree Hill including Dir, Swat, and Azad Kashmir.Large plantations have been raised in Mansehra and Abbottabad areas of KPK.It is mostly found in the form of monospecific community.
is species is easily cultivated and naturally regenerated in the northern area of Punjab and KPK province and considered as fast-growing and valuable species for forestation and reforestation of denuded areas of Pakistan [27].
Subtropical Chir Pine forest is managed under shelterwood management system.Under this system, the area of the forest is divided into different blocks on the base of age carrying with different management operations.Cutting and felling operations are concentrated to overmature blocks, thinning operation are carried out in mature block, and cleaning operations are in the young block.
ese management operations in various blocks may affect the carbon allocation in different carbon pools particularly in the litter, deadwood, and soil.us, complete accurate information is needed regarding biomass and carbon storage in different age stand of the Chir with respect to the current management operations.
erefore, we designed this study with overall objectives of assessing biomass and carbon stock of P. roxburghii forest ecosystem on basis of stand age and outlining the effect of management practices on carbon stock in respective stand age classes.

Study Area.
e study was carried out in Murree Hill of Pakistan.e latitude and longitude of the study area range from 33 °47′15″ to 33 °54′ 47″ N and from 73 °16′ 54″ to 73 °29′ 18″ E. e elevations from sea level range from 939 to 1873 m.
e average precipitation varies from 500 to 1200 mm while temperature ranges from −5 °C in winter to 40 °C in summer in the study area.e rocks are sedimentary in origin and comprise of shales, sandstones, limestone, and marls.e soil is loamy in nature with a variable proportion of sand, silt, and clay [27].e study site is the natural zone of Chir forest that is managed under shelterwood system.e dominant tree species of the area is Chir Pine (PR).e major associated tree species are Pinus wallichiana (kail), Quercus incana (rhin), and Pyrus pashia (batangi).e associated understory flora includes Myrsine africana (khukhal), Berberis lycium spp.(sumblu), Carissa spinarum (granda), Dodonaea viscosa (sanatha), Adhatoda vasica (Bahekar), Capparis decidua (karir), and cannabis sativa (Bang) (Figure 1).

Description of Shelterwood System and Research Design.
Based on the uniform shelterwood silvicultural system, the forest area is divided into four periodic blocks (PBI, PBII, PBIII, and PBIV) on the base of rotation age (100 years) and regeneration period (25 years), thus forming four age stands that are (1-25, 26-50, 51-75, and 76-100 years).Under this system, the first in the mature crop (PBIV), a regeneration felling, is carried out in which all the trees are not felled, but prescribed numbers of trees (seed bearers) are retained over the area for seed production.Mostly, 10-20 trees•ha −1 are retained for the seed production.When the regeneration is established, then the final felling is carried out to remove the seed bearers for the development of the new crop.e regeneration and final felling are completed in 25 years.In PBIII, the thinning operations are carried out with the object to enhance the growth and to prepare the crop for felling.
inning is mostly concentrated in the lower diameter classes.In PBII, and PBI thinning and cleaning operations are carried out.A theoretical sketch of the system is given in Table 1.
Because of the same management operations in PBI and PBII, we placed these two periodic blocks into young age class (1-50 years).Similarly, PBIII represents the mature stand (50-75 years) and PBIV represents overmature stand (>75 years).In each block, sample plots were selected randomly for data collection.e size of each sample plot size was 20 m × 30 m.In each plot, tree diameter (cm) and height (m) along with stem density (ha −1 ) were measured.In each sample plot of the respective block, subplots of 1 m 2 were laid out for soil samples collection.Similarly, samples were also collected from these subplots for the litter, deadwood, and understory carbon analysis.

2
Advances in Meteorology e BEF and R values were obtained from the literature [22,26,[30][31][32]. e whole relations are as follows: where BA � basal area (m 2 •ha −1 ), H � total tree height (m), and F � form factor. e value of FF was sourced from the previous literature [26,29]: For the biomass estimation of herbs, the vegetation was destructively harvested.Similarly, the deadwood and litter were also collected within 1 m 2 subplots following [22][23][24][25][26]33]. e fresh weight of all samples was recorded, and then one Kg subsample was transferred to the laboratory and dried at 72 °C for 48 hours for biomass measurement and then extrapolated to the contemporary weight of the whole plot [34][35][36]: where WS � weight of wet field sample in gram (g), WDS � weight of dry subsample in gram (g), WFS � weight of fresh dry subsample in Gram (g), and A � area in hectare (ha).

Carbon Stocks Assessment.
e carbon stocks values in trees, herbs, and deadwood were calculated from biomass using a conversion factor having value 0.5 (Equation ( 5)) [22,23,[34][35][36][37]: e carbon stocks values in litter were calculated from biomass using a conversion factor having value 0.47 (Equation ( 6)) following IPCC (2006) guidelines: 2.4.Soil Organic Carbon.For the estimation of soil carbon stock, soil samples were taken with soil auger at the depth of 0-15 cm and 16-30 cm in each plot.e weight of each sample was measured in the field, and the samples were put in labeled bags and were brought to the laboratory for further analysis.e soil bulk density (g•cm 3 ) of each sample was measured from the weight of soil sample and volume of the core.Soil carbon content was measured following the procedure of Walkley and Black method [38], and soil carbon in t•ha −1 was calculated from the following equation [22,23,39]: where SBK � soil bulk density (g/cm 3 ), SOC � soil organic content (%), and SHT � soil horizon thickness (cm).

Statistical Analysis.
Statistical software (sigma plot 12.5) was used for the manipulation of statistical data.Regression models were developed for the relationship between diameter (cm) and tree height (m), basal area (m 2 •ha −1 ), and volume (m 3 •ha −1 ), biomass (t•ha −1 ) and carbon stocks (t•ha −1 ).One-way ANOVA (all pairwise comparison) was applied to test the significance of the means values.

Growing Stock and Biomass.
e results revealed that the average stem density ranged from 636 ± 93.7 (trees•ha −1 ) in the young stand to 147 ± 56.7 (trees ha −1 ) in the overmature stand.A statistical higher stem density was found in the young stand while lower stem density was found in the overmature stand.
e basal area in young, mature, and overmature stands was 15.51 ± 1.6 (m 2 •ha −1 ), 35.57± 18.0 (m 2 •ha −1 ), and 45.91 ± 15.1 (m 2 •ha −1 ), respectively.Statistically, no significant variation was found in the mature and overmature stand in the basal area; however, young stand showed significantly lowest basal area.e highest total volume was found in the overmature stand, while lowest volume was recorded in the young stand (Table 2).Further regression analysis between basal area and volume showed a significant relationship among all stands (Table 3).e characteristics of growing stock are figured in Table 2.  Advances in Meteorology e total tree biomass distribution of Chir Pine forest of the young, mature, and overmature stand was recorded at 80.0 ± 12.4 (t•ha −1 ), 343.1 ± 167.6 (t•ha −1 ), and 529.5 ± 176.8 (t•ha −1 ), respectively, demonstrating a rapid increase of biomass from young stand age trees to overmature stand age trees (Table 4).e stem made the largest contribution to total tree biomass.e results of ANOVA showed that the mean value of stem biomass and total tree biomass significantly different among the age classes (Table 4).Moreover, the regression analysis between basal area and biomass revealed strong positive correlation (Table 3, Figure 2).In addition to total tree biomass, we also measured the value of biomass in understory vegetation, deadwood, and litter, and no significant variation was observed among all age classes in their values (Table 4).

Carbon Stocks.
e results show the values of a living tree, understory vegetation, litter, deadwood, and soil carbon of the respective age stands (Table 5).It can be seen from the table that the average living tree carbon varied between 40.0 ± 6.2 (t•ha −1 ) and 264.5 ± 88.4 (t•ha −1 ), respectively.Among the all age stands, significantly larger value was recorded in the overmature stand, whereas smaller value was recorded in the young age stand.e mean value of carbon in understory vegetation was found the maximum in the mature stand.Similarly, deadwood and litter carbon were recorded higher in the overmature stand and lower in the mature stand.However, no significant variation was found in carbon values of understory vegetation, deadwood, and litter in all age stands.Soil organic carbon analysis showed that young stand holds the maximum soil carbon and the overmature stand holds minimum carbon.
e results showed a decreasing trend from young to overmature stand in soil carbon (see Figure 3).

Discussion
e present findings showed that stem density of Chir Pine in all age stands decreases with increasing diameter.In Chir Pine forest, not only natural thinning is common, but stem removal under the shelterwood system is also practiced which are mostly concentrated in the mature and overmature age class [22].
ose areas of the forest where Superscripts in each column show significant differences at α � 0.1 and p ≤ 0.0001.
, and H � height (m).[6,19,40,41], confirming that stand age is the important variable in biomass carbon estimation.However, the mean values of understory vegetation (USV) and deadwood carbon showed contrasting pattern against the living tree biomass carbon.e USV biomass carbon was recorded higher in the mature stand as compared to young stand and overmature stand. is can be attributed to management operations in the forest.It has been discussed earlier that Chir forest in Pakistan is managed under shelterwood system.Under this system, the young crops received the cleaning, and clearing operations hence resulted in lower biomass carbon values.Similarly, in overmature stand heavy thinning, regeneration felling and final felling operations are carried out that may affect the biomass carbon of USV.In contrast, the mature block receives mostly the thinning and selected felling operations.
is creates canopy gaps, thus providing a favorable condition for the growth of USV.Consequently, it gives the larger value of biomass carbon.
e biomass carbon of deadwood was found higher in the overmature block.is higher value of biomass carbon in overmature block reflects the effect of management operations.In the mature block, due to felling operation, the accumulation of more wood materials on forest floor resulted in higher deadwood carbon.
e previous results also reported that biomass carbon for conifer chronosequence increases with stand age [42][43][44], and some researches showed that the carbon stock of forest floor and understory vegetation decreased with increasing stand age [8,19] because of high susceptibility to disturbances, variation, and decomposition rate in stand age [6,8,45].e results of our soil carbon analysis showed decreasing trend with respect to age. is decreasing trend in soil may be attributed to the soil disturbance during the felling operations in the overmature stand and decrease in the natural thinning with the stand age.
e soil organic carbon has many limitations that whether or not soil organic carbon stock could change with respect to stand ages [46,47].Some previous studies showed that there is no significant increase in soil organic carbon stock with stand age [46,17], while some other previous studies showed increasing of soil organic carbon with increasing stand age [45,47,48].e difference may due to different factors such as climatic factor, soil factor, and forest type [6], but in case of the present study, the lower soil carbon with respect to increasing stand age is the result of management operations such as thinning and felling.During these operations, the soil becomes disturbed and the top layer gets exposed to wind and rainfall that wash out the surface soil, reducing the soil organic matter and soil carbon.However, generally, soil organic carbon increases with stand age in conifer forest, due to the accumulation of organic matter in older stages [14,40].
In conclusion, our results depicted that the biomass carbon of living tree (LT) increases with stand age.However, the value of soil carbon, deadwood, and USV carbon did not follow the same trend.Altogether, our current findings figure out that the mean carbon values of all components varied between 90.3 t•C•ha −1 in young, 218.1 t•C•ha −1 in mature, and 309.5 t•C•ha −1 in the overmature stand.ese results indicated the significant increasing trend in carbon stocks with stand age.ese present results highlighted the importance of age in assessing carbon measurement of a forest.Furthermore, our results showed the effects of management operations on the forest floor and soil carbon.
e lower soil carbon value in the overmature stand is the result of soil disturbances during management operations.

Limitations and Uncertainties
is study was conducted using the field inventory method to explore patterns of biomass and carbon allocation across Chronosequence of Chir Pine (Pinus roxburghii) forest in Pakistan.We calculate the stem biomass from stem volume  and wood density and then converted stem biomass to total tree (above ground) biomass using a fixed generic BEF value of 1.51 following Haripriya [30], IPCC [31], and Nizami [22].For measuring the belowground biomass carbon, we used a root to shoot ration (R) of 0.19 following Rana et al. [32].
As the BEF and R value may vary with stand age, the present estimates of carbon allocation along the Chronosequence may contain some uncertainties, which might induce errors in carbon stock estimation along the stand age.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Table 1 :
eoretical sketch of the system.

Table 3 :
Relationship type, equation, and R 2 value of the Pinus roxburghii stands.

Table 4 :
Biomass of Chir Pine forest.

Table 5 :
Carbon stocks of Chir Pine forest.