Biomass and Carbon Stocks Estimation of Lowland Dipterocarp, Riparian and Hill Dipterocarp Forests in Pahang National Park, Malaysia Biomass and Carbon Stocks Estimation of Lowland Dipterocarp, Riparian and Hill Dipterocarp Forests in Pahang National Park, Malaysia

Forest biomass estimation at local or global scale is very crucial and served as an important indicator for monitoring and estimating the forest carbon ecosystem especially in the context of climate change. Pahang National Park (PNP) is considered as a primary forest, and therefore, it is expected that more carbon can be absorbed and stored by forest biomass. Despite the multifunctional roles of forest biomass, lack of research had been done with regard to the extent of above-ground biomass (AGB) and below-ground biomass (BGB) in lowland dipterocarp (LDF), riparian (RF) and hill dipterocarp forests (HDF). Therefore, this study was conducted to provide an estimation of the AGB, BGB and carbon stocks with respect to different localities in PNP. A total of 60 plots were randomly set up and each forest type contains 20 plots measuring 20 × 20 m. The diameter at breast height (DBH) and height (H) were used to calculate the AGB and BGB, and the carbon conversion coefficient of 0.50 was used to calculate the carbon stocks. Based on the results, the estimation of biomass within LDF, RF and HDF not greatly varies between different species with the mean total tree biomass (TTB) values of 415.11, 323.33 and 579.05 t/ha, respectively. The estimation of carbon storage demonstrated that HDF attained the highest carbon stocks in TTB with the value of 289.52 t/ha. The information from this study is expected to provide baseline information and an understanding on the role of trees in the natural forest in sequestrating carbon.


Introduction
Tree biomass is a product from photosynthesis as a result of carbon sequestration by tree. A tree can absorb approximately 23 kg of carbon per year. Indeed, a tree can increase biomass as an effect of tree growth and loss biomass through mortality that is due to natural death or logging. Tree biomass can be divided into above (AGB) and below ground biomass (BGB) in which the AGB includes the stem, leaves and branch biomass whereas the BGB is the biomass of tree roots. Each component of the AGB varies in biomass density. The estimation of biomass is significantly important for the environment which is a critical aspect of studies of carbon stocks and the effects of carbon sequestration on the global carbon balance. In recent years, carbon dioxide (CO 2 ) has received much attention from the world because its concentration in the atmosphere has risen to approximately 30% above natural background levels [1]. The need for biomass and carbon stocks estimation is critical and can be measured using destructive or non-destructive sampling method. That is why a field inventory is conducted where the measurement of tree diameter is recorded to estimate the biomass of tree and later the carbon stocks. According to Brown [2], for closed forest such as Pahang National Park (PNP), a minimum diameter of tree to be measured is greater than or equal to 10 cm. However, for open or secondary forest, a smaller minimum diameter should be chosen [2,3].
Most of the researches focus on the estimation of the AGB rather than BGB because the process to estimate the AGB is easier and less complicated as compared to BGB. In addition, the above ground tree components are the largest contributor of biomass from the total tree biomass (TTB) whereas the BGB only constitutes a small portion of the TTB. Lajuni and Latiff [4] reported the BGB value in their study plots at Khao Chong forest was one tenth of the AGB. Besides, a study conducted by Mohamad [5] in Kenaboi Forest Reserve, Negeri Sembilan found that root biomass in his study plots was six times smaller than the AGB with values of 463.81 and 73.57 t/ha for AGB and BGB, respectively.
Tree biomass and carbon stocks also varied in accordance to forest types and geographical regions. As such, forest biomass and carbon stocks in tropical forest are higher than temperate forest. This might be due to the different in tree species and climatic condition between both forests. Furthermore, in any forest types, tree biomass and carbon stocks in primary forest are higher than secondary forest. Secondary forest is a forest that has been logged or naturally disturbed whereas primary forest is a forest that has never been logged and free from anthropogenic disturbance. In this case, PNP is considered as a primary forest since anthropogenic activities such as logging have never occurred in this forest. Therefore, it is expected that more carbon can be stored by forest biomass in PNP.
Despite the multi-functional roles of forest biomass, lack of research had been conducted with regards to the extent of AGB, BGB and carbon stocks in lowland dipterocarp (LDF), riparian (RF) and hill dipterocarp forests (HDF) in PNP. In addition, information on biomass estimation and carbon stocks from tree inventory data is currently unavailable for protected forest of PNP. Therefore, this study was conducted to provide the estimation of the AGB, BGB and carbon stocks with respect to different localities in PNP. Considering the fact that biomass represents the role of tree as a key indicator of carbon source and sink, the information from this study is expected to provide baseline information and an understanding on the role of trees in sequestrating carbon. This study aims to estimate the AGB and BGB as well as the TTB of LDF, RF, and HDF in PNP. This study also aims to estimate the carbon stocks of LDF, RF, and HDF in PNP, and to investigate the interaction between forest and five similar family and five similar species in the study areas.

Study area and field data collection
This study was conducted in PNP in the state of Pahang. PNP has a tropical climate with an annual rainfall of about 2.260 mm and rich in forest vegetation such as trees, climbers, shrubs, epiphytes and palms. Average temperature throughout the year ranges from 20 to 35°C with more than 80% humidity [6]. There are differences in the soil series in the LDF, RF and HDF mainly due to the variations of parent material between localities [7].
This study was conducted in three types of forests of LDF, RF and HDF of PNP and the location of study area are shown in Figure 1. The description for each location is summarized in Table 1. A total of 60 plots were set up in which each forest contains 20 plots measuring at 20 × 20 m (0.04 ha). Study plots for LDF were located in Kuala Keniam while plots for RF were scattered; 10 plots were located along Keniam River while another 10 plots were located along Tembeling River near to Kampung Pagi. As for the HDF, data collection was conducted in the Teresek Hill at an elevation around 330 m above sea level.
For the field measurement, diameter at breast height (DBH) tape was used to measure the diameter of sampled trees with DBH ≥ 10 cm which is 1.3 m up from the ground [8]. In the case of big buttressed stems, the tree height was measured just above the upper end of plank buttress [9]. Each tree was permanently tagged using laminated label. Tree height was measured using a clinometer, a device that can be used to measure the slope to points on a tree, which can subsequently be used to determine the tree height. The sampled trees were identified to species level and for unknown species, the botanical specimens (e.g. leaves, flower or fruit) were collected for species identification at herbarium laboratories of Universiti Kebangsaan Malaysia (UKM) and Forest Research Institute Malaysia (FRIM).

Data analysis of tree biomass
Throughout this study, AGB was estimated using Kato et al.'s function [9] (Eqs. (1)-(4)) while BGB using a function from Niyama et al. [10] (Eq. (5)). According to reference [10], the total root biomass is the summation of coarse and fine roots in which fine root is defined as root with diameter less than 5 mm. The TTB is the summation of AGB and BGB. From the values of measured DBH and tree height; the dry mass of stem, branch and leaves of sample trees were estimated. The equations used to estimate these components are as follows: Where, M s, M B , M L were denoted as dry mass of stem, branch and leaves in kg, respectively (Eqs. (1)-(4)). The AGB was computed from the summation of these components as in Eq. (4). The biomass functions developed by Kato et al. [9] can be applied irrespective of tree species since these equations were developed without taking regards of tree species in the study area of Pasoh Forest Reserve [9,11]. In this study, the dry mass of the tree biomass components was presented in t/ha. The dry mass (kg) for each component was converted into tonne by dividing the values with 1000 then divided with 0.04 ha which is the size of each plot. For an estimation of carbon storage, the biomass value was divided with 0.8 ha which is the total size for each study area. The carbon storage in the forests was calculated in accordance to method from Brown [2] whereby 50% of the biomass in the forest is assumed as carbon.

Statistical analysis
Means of AGB, BGB and TTB between LDF, RF and HDF of PNP were obtained and analyzed using 3 × 5 factorial two-way ANOVA. The PROC GLM was applied in Statistical Analysis Software (SAS) version 9.3 to study the interaction between forests and five similar family and species based on the highest AGB in LDF, RF and HDF of PNP. The normality of the dataset is test using frequency distribution or histogram. Based on the analysis, the data distribution is normal and the statistical tests are considered as parametric tests.

The total AGB, BGB and TTB for different types of forests
The total AGB, BGB and TTB for lowland dipterocarp, riparian and HDF are shown in Table 2.
From the Table 2, it appears that HDF recorded the highest AGB, BGB and TTB among study areas. This is because HDF consists of higher trees (n = 579) and number of trees with DBH of more than 80 cm was higher than the other two forests (14 trees/ha) ( Table 3). Furthermore, dominant family in HDF based on basal area was Dipterocarpaceae with tree count of 58 from 579 trees ( Table 3). These dipterocarp trees have diameter ranges from 10.8 to 103.5 cm. RF recorded the lowest AGB, BGB and TTB among the three forests as it recorded contains less number of trees (n = 285) and most of trees in RF have smaller diameter. Big-sized trees in RF with diameter more than 80 cm was lower than LDF and HDF (3 trees/ha) thus less contributed to tree biomass of RF.
As comparison with the previous studies, Cairns et al. [12] presented the AGB for 195 sampled trees with diameter of more than 10 cm in dry forest of Mexico's Yucatan Peninsula with value of 191.5 t/ha. Hikmat [13] conducted a study in three virgin jungle reserves in Mata Ayer, Bukit Bauk and Gunung Pulai each in 2 ha plot. A total of 2341, 2702 and 2070 trees with diameter greater than 5 cm were enumerated in Mata Ayer, Bukit Bauk and Gunung Pulai, respectively. From this study, he found that the AGB of each forest was 402.6, 551 and 320.57 t/ha, respectively. The BGB in Hikmat's [13] study was computed following method from [14] in which the root biomass was estimated to be one tenth of the AGB. In this case, the BGB values in Mata Ayer, Bukit Bauk and Gunung Pulai were 40.26, 55.12 and 32.06 t/ ha, respectively. The summation of AGB and BGB in the three study areas resulted in total tree biomass of 415.11, 323.33 and 579.05 t/ha. A study at Bangi Permanent Forest Reserve by Lajuni and Latiff [4] revealed that the AGB in 1 ha study plot was 362.13 t/ha derived from 1018 trees of more than 5 cm diameter. Most of trees in their study were distributed in class 5.0-14.90 cm (65.71%) causing the biomass value to be quite low than other studies. Table 4 shows results from the analysis of AGB, BGB and TTB (t/ha) of lowland dipterocarp, riparian and HDF of PNP. Values presented in Table 4 are mean values of AGB, BGB and TTB per plot. Result from ANOVA revealed that HDF recorded significantly higher mean of AGB, BGB and TTB than LDF and RF with the values of 499.97, 85.27 and 585.25 t/ha, respectively (p ≤ 0.05). This is because HDF comprises the highest number of tree and basal area compared to LDF and RF. Family Dipterocarpaceae contributed 10% from the total individuals in HDF. Mostly, dipterocarp trees in this forest especially Shorea curtisii have tree height ranges from 30 to 45 m and form the emergent layer of the forest. Even though Dipterocarpaceae was not the highest in term of tree density in the forest, they contributed the most in basal area with value of 13.91 m 2 /ha as these trees have larger diameter and height as compared to the other family. This value was the highest compared to LDF and RF. Therefore, this contributed to the higher values of AGB, BGB and TTB in HDF. Generally, basal area indicates the cross section of tree stem at breast height. Therefore, this value can be assumed as proportional to the stem biomass of a tree which also indicates the productivity of a forest [4]. This was supported by a result from Proctor and Newberry [15] in their study in four types of lowland forest in Gunung Mulu. They reported that TTB values in each forest types were in accordance to the value of mean basal area.

The analysis of mean of AGB, BGB and TTB between forests
As for LDF, family Euphorbiaceae recorded the highest density (90 trees/ha), more than the highest family in HDF. However, the basal areas contributed only 3.10 m 2 /ha, considerably lower than family Dipterocarpaceae from HDF. Euphorbiaceae is known as a pioneer species and commonly have small diameter at the range of 10 to 30 cm in this forest. RF on the other hand, recorded the lowest number of trees compared to the other two forests (285 trees). Family Meliaceae recorded 26% (75 trees) from the total of 285 trees in RF and mostly composed of small trees with diameter 10 to 30 cm and seldom can exceed more than 40 cm, thus causing the tree biomass in RF to be lower than LDF and HDF. The AGB value for tropical, temperate and boreal forests were 330-370 t/ha, 155-170 t/ha and 50-60 t/ha, respectively. The heterogeneity in the AGB values between forests was attributed to the climatic factors that affected the soil nutrients in the forest. In this case, due to the seasonality and temperature of boreal forest, nutrient availability is limited by slow decomposition in cold and water-freeze soil. Tropical forest on the other hand, even though has all year warm temperature but have poor soil nutrient and water availability as a result from high soil porosity and heavily leach soil. In general, higher tree biomass is expected on fertile soil simply because there are more resources available for tree growth. According to Laurance et al. [19] a high fraction of forest biomass could be associated with the most fertile soils as well as the tree size. Castilho et al. [20] claimed that texture was strongly associated with the variation in AGB value in their study area at Amazon Forest rather than soil nutrients. Soil texture influences the soil moisture, nutrient availability and nutrient cycling as well. Figure 2 shows the above ground, below ground and total tree biomass in LDF, RF and HDF of PNP, respectively. Based on Figure 2, the total tree biomass in the study areas were not uniformly increased according to diameter class. HDF attained the highest total tree biomass for most diameter class except for diameter class 40.0-69.9 cm. RF achieved the lowest total tree biomass except for diameter class 50.0-69.9 cm whereas LDF only obtained the highest total tree biomass for diameter class 40.0-49.9 cm. With respect to Figure 2, lowland dipterocarp, riparian and HDF acquire highest biomass for diameter class of more than 70 cm with biomass value of 83.39, 70.58 and 202.72 t/ha, respectively. The biomass value for class >70 cm dominated 34, 38 and 35% of the total tree biomass in lowland dipterocarp, riparian and HDF, respectively. This indicates that tree diameter is a deciding factor for producing high biomass value in a forest.

Hill dipterocarp forest (n = 20)
Above   As comparison between diameter class, sample trees at class 10.0-19.9 cm recorded the highest number of trees which is 256, 157 and 344 trees in lowland dipterocarp, riparian and HDF, respectively. However, this diameter class recorded lower biomass even though the number of trees was high. Diameter class >70.0 cm recorded the highest biomass though the number of trees was lower which are 5, 6 and 14 sample trees in lowland dipterocarp, riparian and HDF, respectively. Higher biomass value in diameter class >70.0 cm in HDF was due to large trees from family Dipterocarpaceae that constitute 11 trees of the total 14 trees from this diameter class.
A comparison with other studies indicated a similar result whereby a larger diameter class achieved a higher biomass in the study area. For example, a study from Kusin [21] [16] revealed that diameter class of 30.0-39.9 cm recorded the highest TTB in Air Hitam Forest Reserve. The TTB value of diameter class 30.0-39.9 cm was 232.73 t/ ha whereas for diameter class >70 cm was 151.54 t/ha. The differences of TTB values between diameter classes in their study were due to the different in tree density. Furthermore, the TTB value in their study for five compartments was higher than LDF and RF from the present study because higher number of trees at diameter > 70 cm (15 trees) compared to this study (five trees). A study by [18] at tropical seasonal rainforest in Xishuangbanna, China found that TTB value for diameter class >70 cm was 115.01 t/ha. This value was higher than LDF and RF but lower than HDF in this study. This might be attributed to the different forest type and environmental factor that cause the biomass to be higher.

A comparison of similar tree families between forests
In factorial ANOVA experiment, forest and families are considered as two types of treatments. In each treatment, forest for example, consist of three levels; lowland dipterocarp, riparian and HDF while family has five levels; Anacardiaceae, Burseraceae, Dipterocarpaceae, Euphorbiaceae and Leguminosae. Therefore, in this study, the factorial design is 3 × 5 factorial. The non-significant interaction between forest types and tree families is shown graphically in Figure 3 which indicated by parallel line trend of mean of AGB distribution among families in each forest (P > 0.05). This indicates that the five families in the forest types in this study respond similarly towards the forest types.
From   Since Euphorbiaceae was a fast-growing pioneer species, therefore most of the biomass in their study was recorded by species in this family.
Zani and Suratman [24] attained a similar result in which there was no significant different detected in the mean of AGB between five transect lines (20 × 100 m) in LDF of Kuala Keniam at PNP. In another study, Rayachhetry et al. [25] observed a similar result in a study to quantify the dry weight of the above ground components of Melaleuca quinquenervia trees in three different localities (dry, seasonally flooded, and permanently flooded) at southern Florida. From their study, the effects of locality on the above ground components (total wood, trunk, branch, leaf, seed capsule, and seed) were found to be no significant.

A comparison of similar tree species between forests
The result of comparison of five similar species between forests namely Canarium littorale (Burseraceae), Elateriospermum tapos (Euphorbiaceae), Ochanostachys amentacea (Olacaceae), Pimelodendron griffithianum (Euphorbiaceae) and Shorea leprosula (Dipterocarpaceae) based on AGB in lowland dipterocarp, riparian and HDF of PNP is presented in Figure 4. Table 6 shows mean of AGB for five similar species in lowland dipterocarp, riparian and HDF of PNP. From analysis of variance, there are statistically significant differences in the mean of AGB values both in the forest types and species (p ≤ 0.05). This result indicated there were significant main effects of forests types and species on the AGB value suggesting that the AGB was influenced by the forest types and families. Based on Figure 4, there was significant interaction between forest and species (p ≤ 0.05). This indicates that there is a variation in the AGB value among species. That is to say, species behaves differently in different forest types.
Based on Table 6, Shorea leprosula from family Dipterocarpaceae appeared as the species with the higher AGB value among the five-similar species in lowland dipterocarp, riparian and HDF with the AGB of 0.69, 7.16 and 2.64 t/ha, respectively. Even though the presence of Shorea leprosula in each forest was not the highest, Shorea leprosula managed to attain higher AGB value due to the large DBH of sample trees. The presence of big trees with diameter of more than 80 cm contributed to the higher AGB values especially in the HDF. The mean AGB of Shorea leprosula in RF was higher than LDF and HDF because this forest consists of only two sample trees of Shorea leprosula therefore, mean AGB value per tree was higher. In fact, both sample trees of Shorea leprosula in this forest have large diameter of 69.0 and 69.1 cm that caused of higher AGB values for Shorea leprosula in RF. The higher AGB of Shorea leprosula in the RF is anticipated because the species is most common in lowland forest. That is to say, Shorea leprosula found in RF might be located at the continuum between lowland and RF. That is why only two trees of Shorea leprosula with large diameter were found in the RF plots.
Among the five species, Elateriospermum tapos recorded the highest number of trees in LDF and HDF but greatly lower in RF. Based on this result, it might suggest that Elateriospermum tapos grows abundantly in LDF and HDF rather than RF. However, mean AGB in RF was significantly higher than HDF (p ≤ 0.05) but not significant to LDF. This is due to similar reason as stated in the case of Shorea leprosula.
The AGB values and tree density varies among five similar species between lowland dipterocarp, riparian and HDF might be due to the environmental factors in the study areas (e.g., soil nutrient, topography, water, light). Each species adapts and respond differently to the limiting factors in the area Shono et al. [26]. For example, Elateriospermum tapos favors forest soil that is dry and less preferable on soil that often wet. This might be the reason the tree density of Elateriospermum tapos in RF was lower than the other two forests.
These five-similar species that can be found in all forests in this study were due to the adaptability of these species to the environmental factors in the areas. For example, Shorea leprosula is a dipterocarp species that can easily adapt to full sunlight and fast growing once the seeds have been germinated whereas Canarium littorale was capable to survive in full sunlight and water stress [26]. According to Pereira Da Silva et al. [27], many factors influence the tree growth in the forest. Usually, tropical tree species exhibits different behavior under different environmental conditions regardless of species or families. Macdicken and Brewbaker [28] agreed with this finding in which they found a significant different between site location and species interactions which indicate different environmental requirements for each species. In support to these findings, Brackand and Wood [29] provided a fact that tree growth was influenced by the environmental factors in the forest. Factors such as climatic, soil, topographic and competition combine to create a site. Therefore, the biomass value in a forest is indirectly affected by these factors because tree biomass value depends on the tree diameter.

Carbon stocks
Global climate change is the current major threat to the earth. Due to the rapid deforestation and land clearing and conversion that have been actively taking place since 1850 [3] the emission of CO 2 keeps increasing. Referring to the report from National Research Council [30], these activities contribute 17% from the total CO 2 released to the atmosphere. However, it was reported that forests can remove twice the amount that is lost by deforestation. It was estimated that the total carbon pool in the forest ecosystems approximately 1150 Gt, of which 14% in temperate forests, 37% in tropical forests and 49% is in the boreal forests [3].  Table 7).
The carbon storage in HDF was the highest due to the higher biomass in this forest. This is because the tree density in HDF was higher compared to the other two forests types (n = 579).
As comparison to other study, Hikmat [13] found nearly the same result in three virgin jungle reserves (VJR) in Peninsular Malaysia. Carbon storage in Mata Ayer VJR, Bukit Bauk VJR and Gunung Pulai VJR recorded 221.43, 303.16 and 176.33 t/ha, respectively. In another study, [16] estimated carbon storage in Air Hitam Forest Reserve was 89.57 t/ha. This value was considerably lower than the present study because Air Hitam Forest Reserve was recovering from the past disturbances. Therefore, most of the sample trees were composed of small diameter trees with average diameter of 24.0 cm.
The differences of estimated carbon storage among tropical forests might be due to some limiting factors such as species composition, soil fertility, disturbance history, successional stage and climate Kang et al. [31]. The AGB in the secondary forest would not be the same as the primary forest. Primary forest contains old-growth and large trees since this forest have never been disturbed whereas secondary forest that had been logged or naturally disturbed contains trees with smaller diameter. Therefore, the tree biomass in secondary forest is less than the primary forest. This was supported by Kang et al. [31] who conducted a study to quantify carbon stocks in primary and secondary forests of Bukit Timah Nature Reserve in Singapore. The result from their study revealed that primary forest obtained higher carbon stock than secondary forest with value of 337 and 274 t/ha, respectively. The values in their study were lower than LDF and HDF but higher than RF from this study.

Conclusions
In this chapter, the AGB, BGB and TTB of lowland dipterocarp, riparian and HDF have been estimated. Analysis of AGB, BGB and TTB between forests showed that means of AGB, BGB and TTB values in HDF were significantly higher than LDF and riparian (p ≤ 0.05). The distribution of AGB, BGB and TTB according to diameter class revealed higher AGB, BGB and TTB values in >70 cm class for all forests. HDF was highest in most diameter class except for 40.0-69.9 cm. LDF obtained highest biomass in 40.0-49.9 cm whereas RF for 50.0-69.9 cm. There was no significant interaction between lowland dipterocarp, riparian and HDF and five similar families (i.e. Anacardiaceae, Burseraceae, Dipterocarpaceae, Euphorbiaceae and Leguminosae) with (p > 0.05). However, the interaction between lowland dipterocarp, riparian and HDF and five similar species (i.e. Canarium littorale (Burseraceae), Elateriospermum tapos (Euphorbiaceae), Ochanostachys amentacea (Olacaceae), Pimelodendron griffithianum (Euphorbiaceae) and Shorea leprosula (Dipterocarpaceae) was significant at (p ≤ 0.05). The estimation of carbon storage in the study areas demonstrated HDF attained the highest carbon value in above ground, below ground and total tree with value of 246.88, 42.64 and 289.52 t/ha, respectively.