Mixed Species Allometric Models for Estimating above-Ground Liana Biomass in Tropical Primary and Secondary Forests , Ghana

The study developed allometricmodels for estimating liana stemand total above-ground (TAGB) biomass in primary and secondary forests in the Asenanyo Forest Reserve, Ghana. Liana biomass was determined for 50 individuals for each forest using destructive sampling. Various predictors involving liana diameter and length were run against liana biomass in regression analysis, and R, RMSE, and Furnival’s index of fit (FI) were used for model comparison. The equations comprised models fitted to untransformed and log-transformed data. Forest type had a significant influence (P < 0.05) on liana allometric models in the current study, resulting in the development of forest-type-specific equations. There were significant and strong linear relationships between liana biomass and the predictors in both forests (R > 0.970). Liana diameter was a better predictor of biomass than liana length. Generally, the models which were based on log-transformed data showed better fit (higher FI values) than those fitted to untransformed data. Comparison of the site specific models in the current study with previously published models indicated that the models of the current study differed from the previous ones.This indicates the need for forest specific equations to be used for accurate determination of above-ground liana biomass.


Introduction
Lianas are important structural component of tropical forests [1].They perform a number of ecological functions which help to sustain tropical forest ecosystems [2].Lianas add substantially to plant assemblages in tropical forests in terms of number of species [3] and stem density [4].Apart from contributing directly to species diversity in tropical forests, lianas also play a number of roles which contribute to maintain diversity of other organism [5].Due to relatively high dominance of lianas in tropical forests, they also contribute a lot to forest biomass, especially in heavy liana infested forests.Specifically, lianas can accumulate as high as 30% of total above-ground biomass in tropical forest ecosystems [2].Comparatively, lianas devote much less biomass for stem support than tress, and therefore, they are able to allocate more biomass for their growth compared with trees [2,5].
Consequently, lianas have higher biomass growth than trees [6].Lianas allocate more biomass for leaf production than the amount of biomass allocated to their stems.Because lianas allocate less biomass to their stems they produce less dense wood compared to trees [7].
Estimating tropical forest biomass and determining its dynamics are important aspects of tropical forest ecology.These are more important in areas where changes in composition and structure of forests are apparent.Although changes in biomass levels in secondary forests may occur often due to persistent human disturbance, primary forest biomass may also undergo changes due to other factors which bring about changes in forest composition, structure, diversity, and productivity [8,9].In view of this, it is necessary for biomass studies to be conducted in both tropical primary and secondary forests.However, at the moment, most biomass assessment studies have been limited to secondary forests, neglecting tropical primary forests.Available data on the amounts of biomass within a forest is not only important in determining the amounts of carbon stored by that forest [10] but also essential in assessing the productivity, structure, and conditions of that forest [11].In addition, biomass data on secondary forests could be useful in explaining the effects of deforestation and carbon sequestration on the global carbon balance [10].As indicated above, lianas already contribute substantially to total above-ground biomass in heavily lianainfested tropical forests (cf.[4]).As human disturbance continues to increase in tropical forests, especially in developing countries, lianas would most likely continue to increase in abundance, which could ultimately lead to an increase in the amount of biomass they store.In spite of the ability of lianas to store high amounts of biomass in tropical forest ecosystems, they have been ignored in most biomass assessment studies in the world [12].It appears that many experts in the area are oblivious of the significant contribution lianas can make to forest biomass.This could explain why most recent calls about the need to increase the number of biomass and carbon stock assessments of tropical forests have all been centered on trees only [13].
The easiest and most practical way to determine liana biomass in tropical forests is to employ allometric equations.The use of allometric equations in estimating plant biomass does not only avoid forest destruction but also allows for the estimation of large forest areas.This is because data on the estimators of plant biomass in allometric models can easily be obtained for large areas of forests within a relatively short period of time without harvesting them.Nonetheless, allometric models can sometimes yield biomass estimates that do not reflect the biomass content of a forest [14].Due to the apparent lack of knowledge about the importance of lianas in storing biomass, and their neglect in forest biomass assessments, only a few liana allometric equations have been developed for them.The limited number of liana allometric equations is also partly due to the difficulty in accessing the whole length of lianas from trees [15].In spite of this challenge, conscious efforts should be made to develop many liana allometric equations that can be used to accurately estimate the increasing liana biomass in tropical forests [2].This would ensure that biomass estimates of tropical forests are a true reflection of the actual biomass contents of forests.The availability of plant part allometric equations enable biomass allocation of plant parts to be determined and assessed.Knowledge of biomass allocation changes in different plant parts can be used to assess changes in plant structure and biogeochemical cycles in tropical forests (cf.[16]).In spite of this, only one study [17] has developed allometric equations for liana leaves and total above-ground part.The present study therefore sought out to develop allometric equations for different liana parts.
Ghana has some of the most complex and biodiversityrich tropical forests in which lianas feature prominently [21,22].Nevertheless, there is no allometric equation for lianas in Ghana and also in the whole of Africa.The current study developed allometric models for the estimation of aboveground biomass of lianas in primary and secondary forests within the Asenanyo Forest Reserve, Ghana.Allometric equations in the present study were developed for both stem and total above-ground components of lianas.Total aboveground biomass of lianas is usually made up of stem and leaf (shoot) biomass components.Therefore, development of stem and total above-ground biomass equations in the current study makes it possible to also estimate liana leaf biomass.This would enable changes in relative contribution of various liana parts to total above-ground biomass to be assessed from time to time.This information together with other similar ones from other plant life forms, such as trees, can be useful in assessing forest ecosystem dynamics which can help in developing forest management strategies.The average elevation of the forest reserve is 162 m a.s.l.Daily temperatures range from 20 to 32.9 ∘ C, and the average annual rainfall is 1856 mm.Relative humidity for the area is high (91%).

Field Sampling and Biomass Measurements.
A total of 22 and 24 liana species were destructively sampled from the primary and secondary forest, respectively (Table 1), for development of mixed species allometric equations (from August to December 2012).Sampling was purposely conducted in the rainy season when leaf biomass was highest [17].Lianas were harvested from different habitats (flatlands, slopes, and undulating lands) which were comparable in the primary and secondary forests.Because in each forest type sampling occurred in flatlands, slopes, and undulating lands which were comparable with those in the other forest type, habitat biased sampling was avoided.The sampling within different habitats in each forest type was to ensure that the models produced could be a reflection of habitat variations in the forests.Liana diameter (at 1.3 m from the rooting base) was measured before harvesting the individuals, whereas their length was measured after they were harvested.A total number of 100 liana individuals (primary forest: 50 individuals; secondary forest: 50 individuals) were harvested in the study (Table 1), and their allometric characteristics are indicated in Table 2.Both single and multiple stem liana individuals were included in the study.Liana leaves and stems were separated from each other, and they were sun dried to constant weights over different periods of time (leaves: 3 weeks; stems: 3 months).The constant dry weights of the species were recorded as their biomass.The FI was used to compare models that had different response variables while the RMSE and  2 were used for models with the same response variables.The RMSE and  2 could not be used to compare models with different response variables because they have the potential of producing misleading results for that purpose [23,24].For this reason, the Furnival's index [25] was used to compare untransformed and log-transformed models.The index was computed as follows:

Data
where   () is the derivative of the dependant variable with respect to biomass, MSE is the mean square error of the fitted equation, and the square bracket ([⋅]) is the geometric mean.Comparatively, models with lower FI and RMSE values have better goodness of fit.On the other hand, the higher the  2 value of a model the better its goodness of fit.Due to downward bias which usually occurs when log biomass are back transformed to arithmetic units [26], a correction factor (CF) indicated below [27] was calculated for the models, which could be used to correct them. Consider where SEE is the standard error of the estimate.Linear regression analyses were conducted to determine the relationships between liana biomass and the response variables (diameter, [diameter] 2 , length, log 10 [diameter], or log 10 [diameter] 2 ) in the case of both untransformed and log-transformed data.Analysis of covariance (ANCOVA) was conducted to examine possible differences in regression slopes of models between the forest types.Forest type (primary and secondary) was used as the main factor, whereas the response variables (diameter, [diameter] 2 ,    In the current study, the overall best total above-ground biomass models were determined for the forest types according to the Furnial's index of fit.These were compared with previous total above-ground liana biomass models indicated in Table 3 [4,12,[17][18][19][20], using paired t-tests.The equations were applied to data sets collected from the forests from which the current allometric equations were developed.The data are comprised of 92 (from 32 species) and 100 (from 32 species) liana individuals in the primary and secondary forests, respectively, (diameter range; primary forest: 2-10.7 cm and secondary forest: 2 to 14 cm).Some of the data pairs were transformed (square root and log transformations) to meet t-test assumptions.
Regression analyses and t-test were performed using GenStat software (VSN International Ltd., Hemel Hempstead, UK) whereas ANCOVA was run with Minitab 15 Software (Minitab Inc.).All analyses were conducted at a significance level of 5%.

Results and Discussion
Same set of allometric equations can be developed for use in both primary and secondary forests in areas where forest type does not influence the equations significantly [15,17].However, in this study, there were significant differences in the slopes of regression models between the primary and secondary forest types (ANCOVA;  < 0.05).This suggested forest-specific differences in regressions between the primary and secondary forests, resulting in the development of forestspecific models for the primary and secondary forests.
In each forest, a total of four different models were developed for the estimation of liana stem and total aboveground biomass (Tables 4, 5, 6, and 7).These models were developed based on untransformed (models 1-2, 5-6, 9-10, and 13-14) and log-transformed (3-4, 7-8, 11-12, and 15-16) data.There were strong and significant linear relationships between liana biomass and the various predictors in all the models developed (Tables 4-7; Figures 1 and 2).Although liana diameter and length were good predictors of liana biomass in the models fitted to data on arithmetic scale, liana diameter ( 2 = 0.994-0.996;RMSE = 0.414-0.484)was slightly a better predictor of liana biomass than liana length ( 2 = 0.972-0.981;RMSE = 0.821-1.100) in the current study.The use of liana allometric equations that use length as a predictor of biomass has a practical challenge.
Measuring liana length on the field is impossible unless they are harvested.Therefore, the allometric models of the current study which use liana diameter as a predictor of biomass are recommended for use in liana biomass determination.Logarithmic transformation of data resulted in increased homogeneity of variance compared with data on arithmetic scale (Figures 1 and 2).The decline in homogeneity of variance of the log-transformed models is consistent with previous studies [13,28].
On the basis of Furnival's index of fit, the log-transformed allometric equations performed better than the equations fitted to data on original arithmetic scale for both liana stem and total above-ground biomass models.The higher goodness of fit of the log-transformed models provides evidence in support of another study which indicated that log-transformed models show better fit than models of untransformed data [29].Based on the FI, models 3 and 4 in the primary forest and models 11 and 12 in the secondary forest were the best allometric models for liana stem biomass.
In the same way, models 7 and 8 in the primary forest and models 15 and 16 in the secondary forest were the best models for total above-ground liana biomass.Since the log-transformed models that use diameter are simpler than the models that employ [diameter] 2 , models 3, and 11 are recommended for liana stem biomass estimation in primary and secondary forests, respectively, in line with the Ockham's Razer principle [30].The principle suggests the selection of simplest models if there are many models that have equal goodness of fit.Likewise, models 7 and 15 are recommended for estimating total above-ground biomass of lianas in primary and secondary forests, respectively.The correction factor values determined for the various logtransformed allometric equations (Table 8) were very low (primary forest: 0.3-0.35%;secondary forest: 0.13-0.14%),indicating that downward bias of the models was negligible.A similar finding was reported by Addo-Fordjour and Rahmad [15].One of the best total above-ground models in the primary forest (model 7) was compared with some previously published models by applying them on a data set taken from the same primary forest.The results indicated that on the whole, most of the previous equations underestimated total above-ground biomass of lianas in this study (Table 9).The mean total above-ground biomass per liana species estimated by five of the previous equations were significantly lower (Table 9;  < 0.05) than that predicted by the current study.Nonetheless, one of the previous allometric equations estimated a mean total above-ground biomass that was similar ( = 0.502) to that of the current allometric equation.Another allometric equation selected as one of the best total above-ground model in the secondary forest of this study (model 15) was compared with the previously mentioned equations in the same manner, using a data set collected from the secondary forest.All the previous equations highly over-predicted total above-ground biomass of lianas in the secondary forest (Table 10).The mean total above-ground biomass per species estimated by all the previous equations differed significantly from that estimated by the current model (Table 10;  < 0.05).The equation of Schnitzer et al. [4] which is considered as somewhat general (because it was developed from large data sets from four different countries) also estimated total above-ground biomass that differed considerably from those of the current models in both the primary and secondary forests.Similarly, in a study in Malaysia, a universal equation of Brown [31] highly overestimated total above-ground biomass of trees by as much as 100% compared with the equation developed in that study [32].The wide disparities between the estimation powers of the site specific and general allometric equations of plants observed in this study as well as that of Kenzo et al. [32] indicate that general allometric equations may not provide accurate biomass estimation in specific forests.The differences in estimation of total above-ground biomass between the previous and the current models may be due to differences in liana wood density in the various forests from which the equations were developed [32,33].The comparison of the previous equations to the current ones has revealed the need for site specific models to be encouraged for accurate determination of liana biomass in tropical forests.However, where site specific allometric equations are not available, care must be taken when choosing allometric equations for forests.As much as possible, equations from the same region or continent should be given preference to equations from different continents [15].

Conclusion
The current study developed allometric relationships between liana biomass, diameter, [diameter] 2 , and length for the estimation of stem and total above-ground biomass.Forest type had a significant influence on liana allometric models in the current study, resulting in the development of forest-type-specific equations.Models were developed on data fitted to log-transformed and untransformed data.
In both forest types, log-transformed data fitted better compared to untransformed data.Comparison of the site specific models in the current study with previously published models indicated that the models of the current study differed from the previous ones.

Table 1 :
Number of individuals of liana species in the primary and secondary forests used for the allometric equations.
The study was conducted in the Asenanyo Forest Reserve in the Nkwawie District, Ghana (06 ∘ 26  23  N, 02 ∘ 06  28 2.1.Study Area. W).The forest reserve lies within the moist semideciduous forest zone in Ghana and has both primary and secondary forests.The secondary forest has undergone selective and illegal logging as well as silvicultural treatments in the past, and the relics of these human activities are still evident.The dominant species in the forest are Celtis mildbraedii, Triplochiton scleroxylon, Albizia zygia, and Cedrella odorata.

Table 2 :
Summary of allometric properties of liana individuals in the two forest types (Asenanyo Forest Reserve) used for the study.

Table 3 :
Six previously published allometric equations (total aboveground liana biomass) used in comparing the current allometric equation.All the equations are based on the model form; total aboveground biomass = exp[ +  ln(Diameter)].
and log 10 [diameter] 2 as estimators of liana biomass to obtain models that best fit the data.A series of models were developed using original untransformed and log-transformed (log 10 ) data.In all cases, only models that complied with regression assumptions (homogeneity of variance, linearity, normality, and nonautocorrelation) and showed high goodness of fit ( 2 > 0.97) were retained.
Analyses.Data analyses involved data exploration and model fitting using liana diameter, [diameter] 2 , length, log 10 [diameter], Homogeneity of variance and linearity of data were assessed using residual plots while autocorrelation, and normality were verified with Durbin-Watson statistics and probability plots, respectively.The Furnival's index (FI), root mean square error (RMSE), and coefficient of determination (adjusted  2 ) were used for model selection and comparison.

Table 4 :
Allometric equations of mixed species for estimating liana stem biomass (kg) in the primary forest.

Table 5 :
Allometric equations of mixed species for estimating total above-ground biomass (TAGB) (kg) in the primary forest.

Table 6 :
Allometric equations of mixed species for estimating liana stem biomass (kg) in the secondary forest.

Table 7 :
Allometric equations of mixed species for estimating total above-ground biomass (TAGB) (kg) in the secondary forest.

Table 8 :
Correction factor (CF) of logarithmic models for stem and total above-ground biomass of lianas.

Table 9 :
Comparison of mean estimated total above-ground biomass (per species) between model 7 and previous models (see Table3) with -test in the primary forest.The value in parenthesis represents the mean biomass for the current equation, whereas those outside the parenthesis represent the means of the previous equations.

Table 10 :
Comparison of mean estimated total above-ground biomass (per species) between model 15 and previous models (see Table3) with -test in the secondary forest.The value in parenthesis represents the mean biomass for the current equation, whereas those outside the parenthesis represent the means of the previous equations.