Abstract

Total tree height (H) and diameter at beast height (D) are important independent variables in predicting volume, biomass, and other forest stand attributes. However, unlike D measurement, which is easy to measure with high accuracy, H measurement is laborious. This study, therefore, developed H-D relationships for ten different forest types in Tanzania Mainland. Extents in which climate and forest stand variables explain the variation in H-D allometry were also assessed. A total of 31782 sample trees covering miombo woodlands, humid montane, lowland forests, bushlands, grasslands, mangroves, cultivated land, wetlands forests, and pines and Eucalyptus species plantations were used for model development. The H estimating model without climate and forest stand variables referred herein as “base model” was first developed followed by “generalized model” which included climate and stand variables. All the data were fitted using nonlinear mixed effect modelling approach. Results indicated that generalized H estimating models had better fit than the base models. We therefore confirm a significant contribution of climate and forest structure variables in improving H-D allometry. Among the forest structure variables, basal area (BA) was far more important explanatory variable than other variables. In addition, it was found that the mean tree H tends to increase with the increase of mean precipitation (PRA). We therefore conclude that forest specific generalized H model is to be applied when predicting H. When forest type information is not available, generalized regional model may be applied. Base model may be applied when forest stand or climate information are missing.

1. Introduction

Diameter at breast height (D) and tree height (H) are among the important single tree parameters that describe forest structure and other forest characteristics [13]. Most tree variables, which are difficult to measure directly such as tree biomass and volume, are highly correlated with D and H [1, 4]. Consequently, single tree biomass and volume are estimated by means of allometric models which use D and H as predictor variables. On the other hand, these parameters are also highly correlated. A good example is the relationships between H and D, D and crown diameter, and D and stem height [47].

Often, H has been excluded and only D or a combination of D and wood basic density is often applied to estimate biomass [810]. This assumes that the ratio of H-D, stem taper, and crown mass fraction is constant from smaller to larger trees. However, this assumption may result in systematic error for biomass and other stand level variables estimates since the relationship between H and D is not linear [3]. This is because in tropics H-D relationship varies within the same species with different sizes, stand density, species composition, and site conditions [1, 2]. In addition, while D increases at a more or less constant rate until the tree is lost to mortality, H normally reaches an upper limit asymptotically [11]. Consequently, various H patterns affect other tree variable such as individual tree biomass and volume [9, 12, 13]. A good example is the reported reduction of allometric model standard error and relative root mean square error from 357 kg to 228 kg and from 47% to 33% when applying biomass and volume allometric equation for miombo woodlands of Tanzania Mainland, respectively, with D as sole predictor and with both D and H as predictor variables [4, 14].

Conventionally, when carrying out forest inventory all trees in sample plots are measured for D and few sample trees are measured for H [15]. The measured sample trees are used to develop simple H-D equations and applied to estimate H of unmeasured trees. This approach is often used because measuring H of all trees is laborious and time consuming especially in the tropical forests. The amount of work involved is aggravated by situations encountered in the field such as steep terrain, wide tree crown, and obstruction of tree tops by other trees crowns and leaning trees [16, 17]. Advantage of this approach is that it takes consideration of site conditions, e.g., climate and structural variation at plot level. However, often, the number of sample trees is inadequate to capture the H-D allometry variations which may jeopardize the precision and accuracy of estimated H and other associated forest stand parameters such as volume and biomass (e.g., [8]). Such H estimating models which utilize only D have been applied for decades [1, 1821]. Studies have reported that H-D allometry is regulated by other factors, particularly abiotic factors such as soil nutrients, depth, and availability of water and biotic factors such as those describing the interactions of trees, e.g., competition for soil and light resources [2, 2224]. These factors have been reported to trigger tree species specific phenotypic responses which affect allocation of biomass in tree parts. Irrespective of the reported significant effects of abiotic factors in H-D, there are a limited number of studies particularly in the tropical forests of Tanzania which have attempted to report the effect of abiotic factors on the tree H-D relationship. In this study we consider that in addition to D, H is also modulated by forest structure and climate variable.

At local scale where the climate gradient is relatively negligible, precipitation and temperature at spatial scale may not have effect on H. At large scale where climate variation is substantial there is indication that H-D allometry also varies significantly among sites [2]. It has been shown that environmental resources especially water affect carbon allocation in tree compartments which may eventually affect H-D allometry. For example, Ledo et al. [22] have reported root-shoot ratio to decrease with the increase in water stress. It has been further reported that maximum attained H corresponds to amount of precipitation while for D irregular pattern has been observed [25, 26]. In Tanzania Mainland where in most cases altitude coincides with moisture availability and precipitation, trees found in miombo woodlands at an altitude ranging from 300 m to 1100 m, on average, are shorter than trees in humid montane forest found at an altitude above 1000 m [27, 28]. Despite such remarkable and obvious effect of climate and other associated abiotic factor, to date, there is no study which has attempted to report evidence suggesting the effect of climate on both H and D simultaneously in the miombo woodlands and other vegetation types in Tanzania. On the other hand, forest stand variables which describe density, e.g., basal area per unit area (BA; m2 ha−1), may vary considerably within a given site compared to precipitation. The effect of forest stands variables on H-D allometry had been reported elsewhere [29]. It is well known that trees of the same size in diameter growing in a dense stand are taller compared to tree growing in less dense stand [29] although Zahabu et al. [30] reported taller trees in widely spaced than in closely spaced teak. Since these studies investigated and reported the effect of forest stand variables to only one dimension of the tree, i.e., H, there is no adequate evidence on their effect on H-D allometry.

Robust H estimating models developed from adequate sample trees and which consider other factors that modify H-D allometry are prerequisite for explaining majority of H variations. Few H estimating models have been developed in Tanzania Mainland but with some limitations; e.g., they covered single species [31] or utilized limited number of sample trees [32, 33]. It is until recent where H estimating models for miombo woodlands (1,392 sample trees), humid montane (751), and lowland forest (314) were developed for few selected sites [14, 28]. However, like other previously developed models, they do not incorporate climate variables. Due to heterogeneity of H-D allometry associated with differences in climate and topography in Tanzania, application of these models beyond their range may not be appropriate. In addition, Tanzania has made an effort to develop accurate volume and biomass allometric equations where those utilizing H in addition to D have better performance in terms of mean prediction error [4, 34, 35]. These efforts will be wasted if H is not accurately estimated and therefore risk tree volume and biomass estimates.

This study utilized large data set generated from the first National Forest Inventory (NFI) in Tanzania popularly known as National Forest Resources Monitoring and Assessments (NAFORMA) with the aim of developing H-D relationship and assessing the contribution of climate and forest structure variables to the H-D allometry. The developed models will mainly serve the following purposes, i.e., increase precision of forest biomass and volume which utilize H as predictor in addition to other variables, e.g., D and wood basic density and reduce forest inventory costs associated with H measurements.

2. Methodology

2.1. Site Description

Data used in this study were obtained from a network of plots established by NAFORMA in Tanzania Mainland (Figure 1). Data were collected over the entire Tanzania Mainland covering all primary vegetation types. The country has a diverse climate with mean annual rainfall ranging from below 400 mm to over 2000 mm per annum. The rainfall for large parts of the country is bimodal with short rains from October to December and long rains from March to May. The biodiversity of forests in Tanzania Mainland is high consisting over 10,000 plant species, hundreds of which are nationally endemic. Of the plant species, 305 are identified as threatened in the IUCN Red List, with 276 species classified as endangered (IUCN 2013). The main forest types include miombo woodlands, humid montane forest, lowland forest, bushland, grassland, mangroves, wetlands, cultivated areas, and plantations [36]. Some of these vegetation types are described briefly in subsequent subsections.

2.1.1. Woodlands and Bushlands

The woodlands constitute the largest vegetation type in Tanzania Mainland, occupying about 44.7 million hectares which is equivalent to 93% of the entire forest area [36]. The woodland has three subtypes: Closed (>40% canopy cover), Open (10–40% canopy cover), and Woodland with scattered cropland. The height ranges between 5–20 m although occasionally being taller than 20 m.

2.1.2. Humid Montane and Lowland Forest

The humid montane and lowland forest occupy 995,000 and 1,656,500 ha, respectively, in Tanzania Mainland which is altogether about 5.5% of the countries forest area. These forests are characterized by high richness of flora and fauna and have high catchment values [43].

2.1.3. Grassland

Grassland is another vegetation type possessing marked variety, with four subtypes [44]. For the most part, this type occurs as its subtypes in combination with either a limited Wooded or Bushed component, or with scattered subsistence cultivation. The grassland subtypes are Wooded grassland, Bushed grassland, Grassland with scattered cropland, and Open grassland.

2.1.4. Cultivated Land

The cultivated land is a type of land with four vegetation subtypes [44]. The physiognomy varies widely in accordance with the significance of the tree and crop component associated with each unit. The agroforestry systems which contain permanent tree crops (timber and fruit) that are mixed with permanent and annual agricultural crops (yam, beans, banana, coffee, etc.) are recognized as one vegetation subtype. The tree crops (Grevillea, Albizia, Cordia, Citrus, Acrocarpus, etc.) which form the upper canopy act as shade to the lower canopy crops (banana, coffee, beans). The other cultivated land subtype is where the woody crops are mixed in varying proportions of fruit tree species such as mango, coconut, citrus, and cashew.

2.1.5. Mangroves

Mangroves are forests found in the tropical and subtropical coastlines between 30° south and north of equator [45]. Mangroves comprise trees and shrubs which grow in or adjacent to the intertidal zone. They are found at primarily along sheltered shorelines where freshwater (rainfall or river flow) dilutes the ocean. At the eastern coast of Africa, 14 mangrove species are growing naturally, and 10 among these are found in Tanzania Mainland. Avicennia marina (Forssk.) Vierh, Sonneratia alba J. Smith, and Rhizophora mucronata Lam. are the three most dominant mangrove species in Tanzania [35].

2.1.6. Pine and Eucalyptus Species

Pine and Eucalyptus species are among the planted exotic tree species in Tanzania Mainland. Pine species in Tanzania Mainland include Pinus patula, P. elliottii, and P. caribaea. These species are the dominant species in most of the government and private plantations with about 78% of the total area planted and the remaining 22% is shared among hardwoods and other softwood species. In Tanzania Mainland, the area under Eucalyptus species is estimated to be 25,000 ha [46] of which 4,665 ha. The planted Eucalyptus species include E. saligna, E. grandis, E. camaldulensis, E. globulus, E. viminalis, E. citriodora, E. regnans, and E. microtheca [46].

2.2. Sampling Design and Field Measurements

Data used in this study were collected during the first ever NFI in Tanzania carried out in year between 2010 and 2014. The double sampling for stratification approach was applied. The sampling was designed based on a simulation study described by Tomppo et al. [47]. The first-phase sample consisted of clusters of plots on a 5×5 km grid. The first-phase clusters were stratified based on predicted growing stock, time consumption for cluster measurements, and slope of the terrain. Altogether, the first-phase clusters that contain 6 to 10 plots were assigned to 18 predefined strata. The second-phase samples were systematically selected from the first-phase sample, with different sampling intensities in each of the 18 strata following an optimal allocation procedure with cost functions tailored for each stratum (see. [47]). Greater sampling intensity was allocated to strata with large predicted growing stock and smaller sampling intensity to strata with small-predicted growing stock. Only the clusters selected during the second phase of sampling were measured in the field. The distance between field plots within a cluster was 250 m, while the distance between clusters varies from 5 km to 45 km. For each cluster plots spaced 250 m apart varied from 6 to 10 depending on estimated difficulty to access the cluster (Figure 2). Accordingly, a total of 3,420 clusters with 32,660 plots were adopted. However, during the fieldwork, a total of 3,219 clusters which are equivalent to 30,773 plots were actually measured. Some clusters were not measured for various reasons, mainly difficult terrain and remoteness. Circular fixed area nested plot designs were applied. The nested plots included 1, 5, 10, and 15 m radius concentric plots.

Measurements of D for all the trees within the plot were in accordance with the tree size and distance from plot centre where larger trees where measured in larger plot and vice versa [47]. For all the trees measured in the plot, species names were recorded both botanical and local names. Every fifth tree in the cluster was selected as sample tree and measured for H using Suunto hypsometers. We refer the readers to MNRT [36] for more details on the field measurements.

2.3. Data Preparation

Following data cleaning to remove outliers in the NAFORMA data, a total of 31782 observations were obtained for analysis. During preliminary data analysis scatterplots revealed relationships between H and D varying with vegetation types. Accordingly, for modelling purpose, the entire data were grouped into vegetation types in order to account for variations in H-D allometry that might be attributed by the tree allometry associated with the specific vegetation type. The relationship between response and predictor variables was assessed by means of scatter plots. The H-D plots indicated nonlinear relationship (Figure 3). Summary statistics for variable under consideration for each forest type is presented in Table 1.

2.3.1. Climate Data and Forest Stand Parameters

Climate data of 2.5 min resolution were downloaded from http://worldclim.org/version2 covering a period between 1970 and 2000 [48]. Climate data parameters of interest were mean annual precipitation (PRA), number of months with precipitation below 100 mm (dm) and mean annual temperature (TEA). Candidate stand forest parameter variables tested were stand level competition variables, i.e., basal area (BA, m2 ha−1), a distance-independent tree-level competition index (CI, a ratio of D to the plot-level quadratic mean diameter (e.g., [49]). Summary statistics of climate and stand forest variables are shown in Table 1.

2.4. Modelling
2.4.1. Height Estimating Base Models

Several models describing H-D relationships have been reported in the past [2, 19, 28]. No uniformly best function has been recommended, suggesting that H-D allometry differs significantly among forests (e.g., [21, 50]). We first developed H estimating models which had only as sole explanatory variable herein referred as “base model”. This was necessary so as to make comparison of base models with other H estimating models which include climate and forest stand variables herein referred as “generalized models”. In addition, H base models are required for estimating H when forest structure and climate variables are not available. Furthermore, we fitted H estimating models to all data set combined herein referred as “regional models”. We fitted six common H estimating model forms, i.e., three with two parameters and other three with three parameters (Table 2).

2.4.2. Generalized Tree Height Models

For complex models with several independent variables, i.e., forest structure and climate variables (generalized model), simple log-log model has been preferred partly because it is easy to converge when estimating model parameters [2]. Limitation associated with log-log model include problem of back transformation and lack of flexibility to set maximum H a tree must not exceed [21, 50]. It is through this background that in this study we decided to fit nonlinear models which agree with known biological H-D allometry. We included climate and forest structure variables into six nonlinear H model forms (Table 2) and refitted the models. It was apparent that majority of the fitted model forms failed to converge. Weibull model was the most successful in terms of convergence and fit and therefore was selected to fit generalized models (see [51]). The formulation of three parameters Weibull model is presented in equation (1)where , , and are model parameters to be estimated and is unexplained variation in H. Parameter describes the maximum H the model can predict; describes the rate of change of H with change in D and describes the shape of the curve.

The NAFORMA data represent a hierarchical structure where the field plots were nested within the clusters. In this case, nonlinear mixed effect modelling approach was considered to be an ideal for developing predictive models that will account for dependence of the plots within the clusters. In addition, to account for variation (i.e., heteroscedasticity due to cluster and D) not accounted by the random effects, we also included power variance function structure, i.e., varPower, implemented in the nlme packages of R software for nonlinear and linear models [52]. Preliminary model fitting showed that model with plots nested in cluster was not able to converge probably due to few observations [21] in some vegetation types. Model with only cluster as random variable converged successful. Therefore, the fixed effects (equation (1)) are common to all subjects, while random effect parameters are specific to each subject (in this case cluster) [52], (Zuur et al. 2009). At the beginning, we allowed all model parameters to vary between clusters through inclusion of random effects. The model fitted data well by only allowing random effect on parameter . With the inclusion of random effects, the following equation (2) emerged: where expresses the difference in parameter of cluster i from mean value obtained from (1) or typical cluster.

The next step was to ascertain the effect of forest structure and climate parameters on H-D allometry. We refer this model as “generalized model”. Before fitting generalized H estimation models, we explored relationship of climate and forest structure variables with H in each vegetation type. Preliminary results show that at least one of the dependent variables (BA, CI, PRA, and dm) is highly correlated with H of the studied forest types (Figure 4). We consider forest structure and climate variables to mostly affect asymptotic parameter, i.e., . The parameter was further expanded as shown in (3). The generalized model variables, i.e., PRA, CI, dm and BA were allowed to enter into the model linearly for easy evaluation of their contribution (equation (3)). where are model parameters.

2.5. Model Fitting and Performance Evaluation

Height-diameter models were fitted by considering the contribution of the random effect associated with clusters. In addition, to account for variation (i.e., heteroscedasticity due to clusters and D) not accounted by the random effects, we included power variance function, i.e., varPower, implemented in the nlme packages of R software for nonlinear regression equations [52, 53]. The models were fitted using Maximum Likelihood procedure. To ensure that our modelling strategy has accounted for heteroscedasticity, models were assessed using residuals plots. We assessed the contribution of forest stand and climate variables by looking at the significance of their respective parameters. Variables with nonsignificant parameter were considered irrelevant to H-D allometry.

To present generalized H estimation models for practical application, we selected those with low Akaike Information Criterion (AIC) and evaluated the best first three models. Percentage Root Mean Square Error (SE%) was also assessed. We also assessed the models by means of residuals plots where models with more balanced residuals were selected.

In addition, preliminary findings (Figures 3 and 4) show that climate variables specifically PRA varies among forest types which consequently may affect H-D allometry. Variation of H among forest types of Tanzania Mainland has been reported also by other scholars [27, 54]. Therefore, assuming that the variation in PRA within the same forest type may not necessarily influence H-D allometry significantly, it is expected that other factors such as competition expressed in terms of CI and BA may have the role of modifying H-D allometry. With large data set, collected from NAFORMA, it is expected that large variation associated with species, site condition, and stand structure difference has been adequately covered. Since the regional models have some of these potential variables, we therefore tested mean prediction error (PE %) of the regional models with and without climate and/or stand forest variables to ascertain which one may adequately predict H of all forest types. The mean prediction error was calculated using the following equation: where and H are predicted and observed total tree height, respectively and n is the number of sample trees.

3. Results

3.1. Height Base Model

Performance of base models is shown in Table 3. The SE% were found to range from 25.4% to 41.3% across forest types. Forests with highest SE% were found to be bushlands followed by lowland and montane forests. Forests with low SE% were pines (26.5%) and Eucalyptus (25.5%) plantations and mangroves forest (28.1%) for natural forests. In terms of performance, based on AIC, Model had a good fit in all forest types except for wetland forests. Model (1) was the second in performance. Model (1), i.e., Ratkowsky [41], was selected for wetland forests while for other forests Model , Weibull [42], was selected. Parameter estimates and mode expressions are provided in Table 4.

3.2. Height Estimating Models with Climate and Forest Stand Variables
3.2.1. Regional Generalized Height Models

The findings show that all parameters of the fitted regional models were significant; i.e., all fitted variables had significant contribution in explaining H (Table 5). The magnitude of importance of independent variables decreased in the following order, BA, PRA, dm, and CI. Based on AIC, the model which combines all the independent variables was found to have best fit. Positive parameters for PRA and BA suggest that as precipitation and stand density increases, the mean H also increases (e.g., Figure 5). Trend of mean PRA versus mean H is shown in Figure 4.

3.2.2. Generalized Height Models for Lowland and Montane

For both forest types, all parameters estimates were significant when forest structure and climate variables were included in the model one by one except for dm in the montane forest (Table 5). The magnitude of importance of independent variables for lowland and montane forests decreased in the following order: BA > PRA> CI and BA > CI > PRA> dm, respectively. In many cases combinations which include CI, CI were found to be nonsignificant. For example, for lowland, a model with a combination of BA and CI and PRA and CI consistently yielded nonsignificant parameters for CI. For the lowland, the five best models in decreasing order of model fit had the following combination of variables, i.e., BA, dm, PRA> BA, PRA> dm, PRA > BA, dm > BA while for the montane they were BA, CI > BA, PRA > BA > CI, PRA> CI.

3.2.3. Generalized Height Models for Miombo Woodlands

All variable combinations yielded significant model parameters. Model with all variables had the best fit. The importance of stand and climate variables decreased in following order: BA > PRA> CI > dm. The models fit decreased in the following order of combination: BA, CI, dm, PRA> BA, CI, PRA> BA, dm, PRA> BA, PRA> BA, CI, dm.

3.2.4. Generalized Height Models for Mangroves and Wetlands

Climate variables were not found to influence H-D allometry in both wetlands and mangroves forests. Except for mangroves where BA was a significant variable, forest stands variables did not improve the H estimating model fit in wetlands. This implies that H base models were found to describe H-D allometry adequately for wetland forests (Table 5). In mangroves, model with BA included had better fit than base model.

3.2.5. Generalized Height Models for Bushland, Cultivated, and Grasslands

The PRA in combination with dm and/or CI did not improve model fit for bushlands. Important variables in decreasing order were BA, dm, and PRA. In cultivated lands, BA, dm, and PRA were found to be important in that order in describing H-D allometry while CI was found to be nonsignificant. In grasslands, CI and dm did not add value to the model while BA was the most significant variable followed by PRA. The combination with good fit in decreasing order in the bushlands was found to be BA, CI, dm > BA, PRA > BA, dm > BA, for the grassland, BA, dm, PRA> BA, PRA> BA > dm, PRA > PRA, and for cultivated lands, BA, dm > BA > CI, dm.

3.2.6. Generalized Height Models for Pines and Eucalyptus Species

The most important variables describing variation in H were found to be forest stand variables, i.e., BA and CI for both plantation tree species of Eucalyptus and pines. Climate variables, i.e., PRA and dm, were found to be nonsignificant. When PRA or dm are combined with BA, the climate variables were found to be important in describing H-D allometry for Eucalyptus species.

3.3. Application of the Developed Models

Although H estimating models with climate and forest structure variables proved to be superior, except for wetlands forest (failure to converge), H base models for all forest types are also presented for practical reasons, i.e., to be applied when climate and forest stand variables are not available (Table 4; e.g., [14, 28]). Of the fitted generalized H estimating models (Table 5), we selected at least two best models for each forest type and regional models for further evaluation except for mangrove forest (one model selected) and wetlands (none was selected). We prioritize and select simple and practical models instead of complex models with large number of variables.

The selected models were refitted by allowing stand and climate variables to enter into the model nonlinearly. This is because H increases to asymptote with the increase of climate and stand variables. However, not all the models were able to converge, e.g., bushlands, grasslands, mangroves, and Eucalyptus species, while also for some, we maintained variables getting into the model linearly due to poor model fit. Models in which climate or forest stand variables were not significant are not presented in Table 6. It was evident that BA was an important variable for every vegetation except wetlands. Addition of PRA improved model fit for reginal model, lowland, miombo, and grassland forests (Table 6). The CI improved model fit for montane, miombo, bushland, and pine forests. The PRA, however, was not important to wetlands, mangroves, montane forest, and pine species when BA enters into the model nonlinearly. The number of dry month (dm) improved the model fit for only cultivated land and bushlands.

At this end, the best regional models were found to be that which consist of sole BA and a combination of BA and PRA. Base model had the poorest performance. The residual box plot and mean prediction error (PE %) for the three regional models are shown in Figure 6 and Table 7, respectively. The performance of the model with sole BA and that of BA and PRA when evaluated against forest types did not differ significantly. The prediction error ranged between -1.26% and -21.68%. The two models produced relatively lower prediction error to the bushlands, lowlands, mangroves, montane, and miombo woodlands. The two models consistently produced similar pattern of residuals; i.e., both were biased to wetlands and grassland trees (Figure 6). The PE% for the basic model were found to be relatively high to all forest types except for lowland and mangrove forests (Table 7). This was also evident when evaluating the model by means of residuals (Figure 6).

4. Discussion

4.1. Height Base Models

This study aimed at developing H estimating models for different forest types in Tanzania Mainland using the National Forest Inventory data collected through NAFORMA program. We also assessed how climate and forest stand variables modify H-D allometry. This is the first study with H estimating models developed from sample trees covering the entire country of Tanzania. As such the models developed from this study are robust enough to account for different tree allometry variations which might be attributed by differences in forest types as well as climatic conditions. Several H estimating model forms were tested and evaluated but generally Weibull model turned out to be the best and was able to smoothly converge when climate and forest stand variables are included (Tables 5 and 6). In addition, the fact that the model form has asymptotic characteristics complies with biological features of H growth. This characteristic may explain why Weibull model has been reported to be more appropriate for prediction of above ground biomass and reduces errors for trees of extreme D compared to nonasymptote which tend to overestimate biomass of extreme trees [50]. With such a big data set used in this study where the confidence that most of H-D variation has been captured is high, e.g., good presentation of trees which have reached maximum H growth, the asymptotic Weibull model may be more appropriate than the log-log model.

The Weibull model was found to be superior in all forest types except in wetland. The performance of Weibull models has been widely reported and used to describe tree allometric relationships, especially H-D relationships reported by other scholars [19, 50, 5557]. The SE% for the base model was found to be relatively higher than that of generalized models (Tables 3 and 6). This may be explained by large variations in H which could not be explained by D alone. However, the error margin corresponds to the reported findings elsewhere [19, 28, 57]. Mangroves and plantation tree species, i.e., pine and Eucalyptus, had relatively lower SE% compared to other forest types. Mangroves have rather fewer woody tree species [35] which may translate to low H variation. This also applies to pines and Eucalyptus trees species which have more uniform allometry in plantation.

4.2. Height-Diameter Allometry and Climate

The findings show mean H of all forest types which were affected by PRA and dm in a regular pattern (Figure 2) and show high coefficient of variation of both PRA and dm (Table 1); PRA and dm were also found to affect the H-D allometry significantly. The affected forest types include all data combined, bushlands, cultivated land, grasslands, miombo woodlands, and lowlands (Table 1). The observed trend may be explained by the fact that the reported forest types are widely distributed in the country and therefore exposed to different PRA and patterns (Table 1; [36]). Taking an example of lowland forests, where two subclasses exist, i.e., semideciduous normally found inland receiving relative low rainfall compared to evergreen forests in the coastal regions [58] with more or less comparable tree species composition, it was found that trees in coastal regions were taller than those found inland [10]. Phenotypic plasticity of trees responses to climate and other environmental conditions such as altitude and soil characteristics may explain these differences [59]. This pattern is also supported by Vizcaíno-Palomar et al. [13] who studied the effect of climate and population origin to the H-D allometry and found that some species modify H-D allometry (e.g., Pinus halepensis and Pinus pinaster) when growing in environment with varying PRA while others maintain the H-D allometry (e.g., Pinus sylvestris and P. nigra). However, pines and Eucalyptus species were not affected by PRA and . Dominant pine species grown in Tanzania Mainland at large scale is Pinus patula in southern and northern highlands [60]. These areas experience similar climatic conditions which may explain why this study did not find H-D allometry variation associated with PRA. On the other hand, Eucalyptus species is widely distributed in both dry and wet areas where we would expect to find the effect of PRA and dm on their H-D allometry. However, the effect of PRA and dm were insignificant suggesting inconsistent H-D allometry response to climatic condition [61, 62].

Furthermore, when comparing mean H of forest types, the findings show that mean H increases with the increase of PRA (Figure 5) and decreases with the increase of dm with few exceptions, i.e., grasslands (Figure 4). Generally, these findings are in line with those reported by Burgess et al. [27] who found that maximum H decreases with decreasing PRA in this order, i.e., montane, lowland, closed woodlands, and open woodlands. The contradicting results with grassland maybe due to the fact that the vegetation type is mainly found in the southern highlands which experience high PRA and also in patches woodlands with moderate PRA. In the grasslands trees are normally scattered and therefore it is likely that they are not struggling for light resource for vertical growth and therefore invest more on branches (horizontal growth; [63, 64]). This is further justified by nonsignificant parameter for CI as shown in Table 5 since the trees experience no competition. On the other hand, although wetlands and mangroves forests receive moderate PRA, it was expected that trees would be taller than those found in other forests since these areas are wetter throughout the year. This implies that there are other factors other than PRA such as genetics, soils fertility and depth, winds, and humidity which also modify H-D allometry which were not captured in our data set [2, 59]. Trend observed in this study provides vital information on tree allometry response to climatic condition. Climate fluctuations are even more urgent because climate changes are expected to cause stronger climatic variability, with increased frequency of extreme dry years which will affect distribution, allometry of tree species and consequently forest biomass [6568].

4.3. Height-Diameter Allometry and Forest Structure Variables

The findings from this study show that BA in most cases (except for wetlands forests) have a strong influence on H-D allometry, i.e., increasing of H with increasing BA. This pattern suggests that H-D allometry varies within forest type, i.e., among plots. This implies that when excluding BA as a predictor, different H-D equation would be required for each plot [69]. There are contradicting findings reported on the effect of stand density to H (e.g., [70]). Other scholars have claimed that the vertical and horizontal dimensional change of a tree is linked with the availability of resources: solar radiation, water, and nutrients. Consequently, site quality has been conventionally determined by dominant H which is considered to be less sensitive to stand density such as BA [71]. Dunning [72] and Baker [73] have shown that H may be affected by stand density. However, the inconsistence of the reported findings has been shown to be attributed by size of plots and variable responses of different tree species to stand density [70]. Data collected in larger plots have shown to be more explanatory to this phenomenon than those collected in small plots [70]. Due to the consistence of our findings, we rule out that the plot size applied in this study was adequate to explain H-stand density relationships. In addition, while few tree species are shade tolerant, light demanders in a tightly packed setup may be forced to shoot upward in pursuit of sunlight [74, 75]. On the other hand, trees growing in a densely packed stems are not susceptible to wind damage due to sheltering effect and therefore they are likely to invest in H than D since mechanical support may not be the priority [2, 76].

While increase in BA aggravates competition for resources and affects tree growth, the magnitude of the impact may be experienced differently from one tree to another depending on their sizes [77, 78]. This is when the tree status in terms of CI becomes important, thus explaining how a tree is affected by neighboring trees. It is expected that tree with large CI to have “proportionally” large H, i.e., larger trees have more advantages in resources capturing than the small trees [22]. However, this interaction favors D growth rather than H. For example, Mugasha et al. [79] reported diameter increment to be positively linked to CI. This imply that while small trees have narrow ground for resources since adjacently they are outmatched with larger trees [80], they capitalize on H growth in pursuit for sunlight [75]. Therefore, it is not surprising that in this study CI was inversely proportional to H in all forest types except for wetland where CI was not significant but maintained similar trend. Insignificant CI parameter for wetland forest may be associated with the limited number of observations used in modelling, i.e., 97 observations (Table 1).

4.4. Application of Height-Diameter Models

Prediction of H is critical step when estimating other forest parameters such as forest biomass and volume (e.g., [1]). The accuracy of the estimates will depend on the prediction power of H estimating models. Comparison between the fitted base and generalized H estimating model for H predictions of forest types indicated that for montane, mangrove, and wetlands, parameters for PRA were not significant suggesting that when BA enters into the model nonlinearly, the model explains most of H and PRA remain redundant. We dropped models having insignificant parameter for PRA. It was expected that the pattern displayed by PRA to the lowland and humid montane forest to be the same as shown in Table 5 (positive and significant PRA parameters). However, while PRA was important for lowland forests in modifying H-D allometry, it was not the case for humid montane forest. This may be explained by wide distribution scale of lowland forests to various precipitation patterns [58]. The fact that humid montane forests occur in wetter areas, i.e., leeward of eastern, southern and northern mountains (elevated areas and therefore low evapotranspiration; [43]), may explain why PRA were found insignificant to H-D allometry.

As expected climate and forest structure variables were highly significant when fitting generalized regional H estimating model since the modelling data covers large spatial variability in terms of climate and forest structure which consequently affects H-D allometry (e.g., [2]). However, the magnitude of their contribution in explaining H-D allometry varied. The residuals of generalized regional model with BA included were far more balanced than base model. Addition of PRA on top of BA did not improve the residual plots. This is further indicated by PE% values where generalized regional model with BA had PE% below 12% except for grassland and wetland forest. These findings indicate that generalized regional model can adequately predict H of trees found in the studied forests with exception of few when climate and forest stand variables are available [2]. For improved accuracy, however, forest specific generalized H estimating models are recommended. Forest specific H estimating model is recommended when climate and stand variables information is not available.

5. Conclusion

The H-D allometry was developed with comprehensive data set collected from NFI covering eight main forest types in Tanzania Mainland. The findings indicated that spatial variability in H-D allometry in Tanzania Mainland is significant. The mean tree height tends to increase with the increase of mean precipitation (PRA). Total tree height increases with the increase in BA. The BA in addition to D was far more important explanatory variable than other climate and forest stand variables. We therefore confirm a significant contribution of climate and forest structure variables in modifying H-D allometry. Consequently, generalized H models turned out to be the best compared to base models. It is therefore recommended the forest specific generalized H model to be applied when predicting H. When forest type information is not available, generalized regional model may be applied. Height base model may be applied when forest stand or climate information are missing.

Data Availability

The National Forest Inventory data also known as NAFORMA used to support the findings of this study were supplied with agreement with three institutions, i.e., National Carbon Monitoring Centre (NCMC), Tanzania Forest Service (TFS), and Tanzania Forestry Research Institute (TAFORI) who are custodian of the said data set. The data set cannot be made freely available. Requests for access to these data should be made through the following address: National Carbon Monitoring Centre, College of Forestry, Wildlife and Tourism, Sokoine University of Agriculture, P.O. Box 3009, Morogoro, Tanzania. E-Mail: NCMCTanzania@sua.ac.tz.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors wish to thank Tanzania Forest Service Agency (TFS) in the Ministry of Natural Resources and Tourism (MNRT) for availing the tree height-diameter data from the National Forest Assessments and Monitoring (NAFORMA) database to the INFORES project supported by Ministry of Foreign Affairs of Finland [grant number: 28235767 and ID: UHA2013-008402] which supported the process of data analysis and finally publishing the paper.