Testing the suitability of the Nepal-Somers stand table projection method for Eucalyptus grandis plantations in South Africa

This study presents an application of a stand table projection method for varying stand densities of Eucalyptus grandis in South Africa. The projections of diameter at breast height (dbh) class frequencies are compatible with future values of basal area and stem number per hectare. The application requires an individual tree dbh growth function, a current stand table, and estimates of future basal area, and tree survival. The resulting stand table is adjusted by an algorithm that guarantees that the future stand table is compatible either with observed or predicted values of basal area and stem number per hectare. The results, based on a relatively small experimental dataset, are satisfactory. However, further work based on a larger dataset is required to evaluate the general suitability of the method for stand table projection of Eucalyptus grandis in South Africa, especially regarding variable site conditions. The particular contribution of this study is the fact that we were able to demonstrate the suitability of the Nepal and Somers stand table projection method for a unique dataset from a very densitysensitive Correlated Curve Trend spacing trial. We also made an attempt to present a comprehensive example of a projection to facilitate practical application of this method in the KwaZulu-Natal coastal plain of South Africa.     

Characterising the diameter distribution of Sal plantations by comparing normal,lognormal and Weibull distributions at Tilagarh Eco-park,Bangladesh

For many years foresters have been using statistical probability density functions to describe and characterise stand structure. Predicting the current and future yields of a stand is essential for successful stand and timber management. Implicit prediction of current yield is accomplished by using diameter distribution methods. All diameter distribution yield systems predict the number of trees per unit area by diameter class. In this study, the normal, lognormal and the three-parameter Weibull probability density function were compared to characterise the diameter distributions of Sal (Shorea robusta) plantations grown at Tilagarh Eco-park, Bangladesh. Data from 70 plots, established in three plantations, were used for this study. The Weibull parameters were estimated by the maximum likelihood and moments estimator methods. A one-sample Kolmogorov–Smirnov test was used for the goodness of fit for all models. The Kolmogorov–Smirnov test results showed that both lognormal and Weibull distributions were suitable to characterise the diameter distributions of Sal plantations in the study area and may be applicable for other Sal forests in Bangladesh.     

Modeling diameter distribution of Austrian black pine (<Emphasis Type="Italic">Pinus nigra</Emphasis> Arn.) plantations: a comparison of the Weibull frequency distribution function and percentile-based projection methods

The main purpose of the present investigation is to examine and compare three methods for diameter distribution modeling in terms of their fitness to predict from stand level variables the diameter distributions of even-aged Austrian black pine (Pinus nigra Arn.) plantations in Bulgaria. The percentile-based projection method involving empirical probability density function based on 12 percentiles was the first method tested. A new modified approach based on the first method was proposed as the second alternative. The third method was the 2-parameter Weibull functional model in which parameters were recovered from the first and the second raw moments and the second central moment of the empirical distributions. The Kolmogorov–Smirnov test was applied to compare the experimental distributions with the predicted ones, and estimation of the error indices was employed to evaluate the total absolute deviation of the predicted numbers from the actual ones by diameter class. The two-parameter Weibull function proved superior to the examined alternative percentile-based projection methods and the newly proposed percentile method, without a driver percentile showed improved precision compared to the classical percentile method (with a driver percentile). The parameters of the Weibull frequency distribution function can be easily recovered from the stand quadratic mean diameter. Consequently, this diameter distribution model could be incorporated as a sub-model for stand horizontal structure characterization within the Stand Density Management Diagram modeling framework.     

Comparing normal,lognormal and Weibull distributions for fitting diameter data from Akashmoni plantations in the north-eastern region of Bangladesh

Statistical probability density functions are widely used to model tree diameter distributions and to describe stand structure. The objective of this study was to compare the performance of normal, logarithmic-normal and three-parameter Weibull distributions for fitting diameter data from Akashmoni (Acacia auriculiformis A. Cunn. ex Benth.) plantations grown in the north-eastern region of Bangladesh. Data from 96 plots, established in 24 plantations in north-eastern Bangladesh and ranging in age from 1 to 6 years, were used for this study. The parameters of the Weibull distribution were calculated using maximum likelihood estimation (MLE) and moment estimation (ME) methods. The goodness of fit of normal, lognormal, Weibull MLE and Weibull ME were tested using one-sample Kolmogorov-Smirnov (KS) tests. The KS test results showed that both lognormal and Weibull distributions were equally effective for describing the diameter distributions of these Akashmoni plantations grown in the north-eastern region of Bangladesh.     

Modeling diameter distributions of Gmelina arborea plantation in Omo Forest Reserve,Nigeria with Johnson’s SB

Diameter distribution modeling is an essential tool for obtaining reliable information on the structure, growth, and yield of forests. However, this tool has not been duly utilized in many forest plantations, especially in Nigeria. Thus, this article presents the effectiveness of using the Johnson S_B function for modeling diameter distributions of a Gmelina arborea plantation in Omo Forest Reserve, Nigeria. Fifty plots of 0.04 ha were randomly selected across three age series. All trees with DBH ≥ 5.0 cm in the selected plots were measured. Three fitting methods for the S_B distribution were compared: conditional maximum likelihood (CML), moments (MOM), and Knoebel and Burkhart (KB) methods. The parameter recovery model (PRM) was used to recover the parameters of the distribution. The assessment was based on Kolmogorov-Smirnov (K-S), mean square error (MSE), mean absolute error (MAE), and bias. The results showed that the underlying diameter distribution of the Gmelina arborea followed the Johnson’s S_B distribution fitted with MOM. The number of rejections by the K-S test was zero for MOM and CML. The mean K-S value for MOM was statistically different from CML and KB methods. The results obtained after recovering the parameters of the S_B distribution were comparable with the real distributions of the Gmelina arborea stand.     

A comparison of diameter distribution models for Khaya ivorensis A.Chev. plantations in Brazil

The purpose of this study was to compare Beta, Gamma, Johnson's S_B and Weibull functions fitted by different methods for describing the horizontal structure of Khaya ivorensis (African mahogany) plantations in Brazil. The database comprised 128 plots from six plantations at varying ages. The function fits were compared using the Kolmogoroff–Smirnoff test, mean bias and mean absolute error for the number of trees and basal area per hectare per diameter class. Johnson's S_B outperformed the other functions, although all functions provided an adequate fit. The best methods were method of moments and maximum likelihood fitted using 25% of the minimum observed diameter as the location parameter for the Johnson's S_B function. The errors were greater in diameter classes with higher frequencies. Location and scale parameters were highly correlated with mean diameter and age for the Weibull and Johnson's S_B functions, respectively, which is convenient for diameter prediction. Gamma's scale parameter had medium correlation with age. Beta's parameters had low correlation with stand attributes assessed.     

Estimating Stand Tables from Aerial Attributes: a Comparison of Parametric Prediction and Most Similar Neighbour Methods

Parameter prediction and the most similar neighbour (MSN) approaches were compared to estimate stand tables from aerial information. The study was based on 50 stands in the south-eastern interior of British Columbia, Canada. In the parametric prediction approach, stand tables were estimated from aerial attributes and three percentile points (16.7, 63 and 97%) of the diameter distribution. In the MSN analysis, stand tables were estimated from the MSN stand that was selected using 13 ground and 22 aerial variables. The accuracy of these approaches was evaluated by comparing the observed and estimated species composition, stand tables and volume per hectare. While the parametric prediction approach is easier and flexible to apply, the MSN approach provided reasonable projections, lower bias and lower root mean square error.     

Predicting and calibrating diameter distributions of Eucalyptus grandis (Hill) Maiden plantations in Zimbabwe

Diameter distribution models for even-aged Eucalyptus grandis plantations in Zimbabwe were developed using the two-parameter Weibull function. The analysis was based on data from Correlated Curve Trend (CCT) experiments replicated on four different sites. Parameters of the Weibull distribution were predicted using stand characteristics as regressors. Two sets of parameter models were estimated: a set with and one without stand basal area as a predictor. Stand variables such as dominant height, age, site index and number of stems were used in both sets. The models were further calibrated to result in a given number of stems and stand basal area simultaneously. The usability of constructed models was tested both in prediction of yield in a stand inventory situation and in simulation of growth in connection with different growth models. The results indicated that models not including stand basal area produce considerably less precise stand volume estimates compared to models including also stand basal area. Calibration improved the accuracy of diameter distribution models. In growth simulation diameter distribution models can be connected both to single tree growth models and to stand projection models. The usability of calibration in growth simulation depends on the accuracy of the prediction of stand characteristics.     

Estimating percentile-based diameter distributions in uneven-sized Norway spruce stands using airborne laser scanner data

Abstract

A model for prediction of stand basal area and diameters at 10 percentiles of a basal area distribution was estimated from small-footprint laser scanner data from primeval conifer forest using partial least squares regression. The regression explained 44–80% and 67% of the variability of the 10 percentiles and stand basal area, respectively. The predicted percentiles, scaled by the predicted stand basal area, were used to compute diameter distributions. A cross-validation showed that the mean differences between the predicted and observed number of stems by diameter class were non-significant (p>0.05) for 22 of 29 diameter classes. Moreover, plot volume was calculated from the predicted diameter distribution and cross-validation revealed a non-significant deviation between predicted and observed volume of ?3.3% (of observed volume). An independent validation showed non-significant mean differences for 20 of 21 diameter classes for data corresponding to the model calibration data. Plot volumes calculated from the predicted diameter distributions deviated from observed volume by ?4.4%. The model reproduced diameter distributions corresponding to the model calibration data (uneven-sized forest) well. However, the model is not flexible enough to reproduce normal and uniform diameter distributions. Volume estimates derived from predicted diameter distributions were generally well determined, irrespective of the observed distribution.     

Growth and yield prediction models for Acacia mangium grown in the plantations of the central region of Bangladesh

Acacia mangium is a very fast growing species belonging to the family Leguminosae that has been introduced in the plantations in Bangladesh for its faster growth and wide range of adaptability. The present study aimed at development of growth and yield prediction models for the species using simultaneous equation method. Models were selected for the species to estimate stand dominant height, stand diameter, stand basal area per hectare and total volume yield per hectare. Paired t-test, 45-degree line test, percent absolute deviation and biological principle of stand development were used for the validation of chosen models. The results suggest that the models derived were statistically and biologically acceptable and could be satisfactorily used for stands of Acacia mangium of ages 4–7 yrs based on a base age of 6 yr.     

Regional mixed-effects height–diameter models for loblolly pine (<Emphasis Type="Italic">Pinus taeda</Emphasis> L.) plantations

A height–diameter mixed-effects model was developed for loblolly pine (Pinus taeda L.) plantations in the southeastern US. Data were obtained from a region-wide thinning study established by the Loblolly Pine Growth and Yield Research Cooperative at Virginia Tech. The height–diameter model was based on an allometric function, which was linearized to include both fixed- and random-effects parameters. A test of regional-specific fixed-effects parameters indicated that separate equations were needed to estimate total tree heights in the Piedmont and Coastal Plain physiographic regions. The effect of sample size on the ability to estimate random-effects parameters in a new plot was analyzed. For both regions, an increase in the number of sample trees decreased the bias when the equation was applied to independent data. This investigation showed that the use of a calibrated response using one sample tree per plot makes the inclusion of additional predictor variables (e.g., stand density) unnecessary. A numerical example demonstrates the methodology used to predict random effects parameters, and thus, to estimate plot specific height–diameter relationships.     

Growth and structure development of semi-natural larch-spruce-fir (Larix olgensis–Picea jezoensis–Abies nephrolepis) forests in northeast China: 12-year results after thinning

Analyzing and understanding the structure and growth dynamics of semi-natural plantations is useful for their management. Since 1987, 16 plots with 4 treatments (CT: control; LT: light thinning; MT: medium thinning; and HT: heavy thinning) by 0, 20, 30 and 40% of basal area removal, respectively, and four replications were established in semi-natural larch-spruce-fir forests in northeast China. The structure and growth dynamics of semi-natural larch-spruce-fir stands and the effects of thinning on the growth, structure and diversity were examined. A mixed model repeated measures analysis of variance (RMANOVA) was used to test the effects of treatment and time.

Results showed that differences in periodic annual increment (PAI) of stand basal area and volume and the individual diameter and volume among treatments changed over time in a complex statistical interaction. Thinning, however, had a significant effect on growth at tree and stand levels 12 years after thinning while the PAI of the diameter, basal area and volume was positively correlated with thinning intensity. No significant differences were found in the total stand yield among treatments. Composition of tree species group (larch, other conifers and deciduous trees) during monitoring years did not change significantly. Moreover, no significant differences were observed in tree species and size diversity among treatments in the years following thinning. Both thinning and control plots had similar understory plant diversity after the 12 year period. Univariate point pattern analysis revealed that clumped and random distributions were dominant for tree species groups in this study. The current species composition and regeneration dynamics within these semi-natural plantations suggested a development towards mixed coniferous and broad-leaved forests. Management implications for the transformation from larch plantations towards mixed broad-leaved Korean forests with a more diverse structure, the climax vegetation in this region, were discussed.     

Study on tracheid variation ofPicea koraiensis (nakai) in natural stand and plantation

Tracheid characteristics ofPicea koraiensis from natural stand in Liangshui area and plantation in Mao'ershan area were studied. The results of variance analysis showed that the tracheid length, diameter, and thickness of tracheid wall were significant differences between the growth rings. All those characteristics tend to increase from heart-center to bark. In natural stand, it has positive correlation between the characteristics. In plantation, the tracheid length has positive correlation with the trachied diameter, while both the trachied length and the trachied diameter has negative correlation with the thickness of trachied wall. The tracheid length and diameter have no significant difference in growth between natural stand and plantations. The thickness of tracheid wall from the plantation is a little thicker than that from the natural stand before 15 years. Foundation item: This paper was part of State 9th Five-Year Plan project of “Early Mensuration ofPicea koraiensis (nakai)”. Biography: REN Xu-qin (1973-), female, graduate student in Northeast Forestry University, Harbin 150040, P. R. China. Responsible editor: Song Funan     

Modelling stand basal area growth for radiata pine plantations in Northwestern Spain using the GADA

A stand basal area growth system for radiata pine (Pinus radiata D. Don) plantations in Galicia (Northwestern Spain) was developed from data corresponding to 247 plots measured between one and five times. Six dynamic equations were considered for analysis and both numerical and graphical methods were used to compare alternative models. The equation that best described the data was a dynamic equation derived from the Korf growth function by the generalized algebraic difference approach (GADA) and by considering two parameters as site-specific. This equation was fitted in one stage by the base-age-invariant dummy variables method. The system also incorporated an equation for predicting initial stand basal area, expressed as a function of stand age, site index, and the number of trees per hectare. This information can be used to establish the starting point for the projection equation when no inventory data are available. The effect of thinning on stand basal area growth was also analyzed and the results showed that the same projection equation can be used to obtain reliable predictions of unit-area basal area development in thinned and unthinned stands.     

Allometry and evaluation of in situ optical LAI determination in Scots pine: a case study in Belgium

We evaluated several optical methods for in situ estimation of leaf area index (LAI) in a Belgian Scots pine (Pinus sylvestris L.) stand. The results obtained were compared with LAI determined from allometric relationships established in the same stand. We found high correlations between branch cross-sectional area, diameter at breast height (DBH) and basal area as dependent variables, and leaf mass, needle area and crown projection as independent variables. We then compared LAI estimated by allometry with LAI determined by three optical methods (LAI-2000, TRAC and digital hemispherical photography) both before and after corrections for blue light scattering, clumping and non-leafy material. Estimates of stand LAI of Scots pine ranged from 1.52 for hemispherical photography to 3.57 for the allometric estimate based on DBH. There was no significant difference (alpha = 0.01) between the allometric LAI estimates and the optical LAI values corrected for blue light scattering, clumping and interception by non-leafy material. However, we observed high sensitivity of the optical LAI estimates to the various conversion factors, particularly to the clumping factor, indicating the need for caution when correcting LAI measured by optical methods.     

Comparison of nearest neighbour approaches for small area estimation of tree species-specific forest inventory attributes in central Europe using airborne laser scanner data

The three nonparametric k nearest neighbour (kNN) approaches, most similar neighbour inference (MSN), random forests (RF) and random forests based on conditional inference trees (CF) were compared for spatial predictions of standing timber volume with respect to tree species compositions and for predictions of stem number distributions over diameter classes. Various metrics derived from airborne laser scanning (ALS) data and the characteristics of tree species composition obtained from coarse stand level ground surveys were applied as auxiliary variables. Due to the results of iterative variable selections, only the ALS data proved to be a relevant predictor variable set. The three applied NN approaches were tested in terms of bias and root mean squared difference (RMSD) at the plot level and standard errors at the stand level. Spatial correlations were considered in the statistical models. While CF and MSN performed almost similarly well, large biases were observed for RF. The obtained results suggest that biases in the RF predictions were caused by inherent problems of the RF approach. Maps for Norway spruce and European beech timber volume were exemplarily created. The RMSD values of CF at the plot level for total volume and the species-specific volumes for European beech, Norway spruce, European silver fir and Douglas fir were 32.8, 80.5, 99.0, 137.0 and 261.1%. These RMSD values were smaller than the standard deviation, although Douglas fir volume did not belong to the actual response variables. All three non-parametric approaches were also capable of predicting diameter distributions. The standard errors of the nearest neighbour predictions on the stand level were generally smaller than the standard error of the sample plot inventory. In addition, the employed model-based approach allowed kNN predictions of means and standard errors for stands without sample plots.     

A hybrid model for intensively managed Douglas-fir plantations in the Pacific Northwest,USA

Recent advances in traditional forest growth models have been achieved by linking growth predictions to key ecophysiological processes in a hybrid approach that combines the strengths of both empirical and process-based models. A hybrid model was constructed for intensively managed Douglas-fir plantations in the Pacific Northwest, USA, by embedding components representing fundamental physiological processes and detailed tree allometrics into an empirical growth model for projecting individual tree and stand development. The simulated processes operated at a variety of scales ranging from individual branches to trees and stands. The canopy structure submodel improved predictions of leaf area index at the stand level when compared to allometric and other empirical approaches (reducing mean square error by 30–42%). In addition, the hybrid model achieved accuracy in short-term volume growth prediction comparable to an empirical model. Biases in 4-year stand growth predictions from the hybrid model were similar to those from the empirical model under thinning, fertilization, and the combination of these treatments; however, volume growth predictions in unmanaged plantations averaged approximately 36% less bias. These improvements were attributed to detailed information on crown structure (i.e. size, location, and foliage mass of primary branches), simple representation of key physiological processes, and improved site characterization. Soil moisture, temperature, and nitrogen mineralization predicted by the hybrid model also agreed closely with observed values from several previous studies. Overall, the model framework will be helpful for future analyses as it can lend insight into the influence of weather and site edaphic factors on growth, help identify mechanisms of response to silvicultural treatments, and facilitate the design of sound management regimes for Douglas-fir plantations across the Pacific Northwest region.     

Yield capacity of planted Picea abies in northern Sweden

Norway spruce (Picea abies (L.) Karst.) is the dominant species in the older forests of central northern Sweden. However, spruce has seldom been planted in the area, partly because existing tools for site classification have indicated a low yield capacity for the species. The aim of this study was to examine the yield capacity of spruce on the basis of existing plantations. In total, 91 operational and experimental plantations in the age interval 27–46 yrs were sampled. Stands were located between 62° and 65° N at altitudes 130–620 m a.s.l. Site index was estimated by height growth and site‐factor equations previously developed from old‐growth stand data. Height developments in the plantations indicate that site index for these stands is on average 4.6 m higher than predicted by site‐factor equations. The differences between the two methods are larger on poor sites than on rich sites. No systematic deviations of top height development from the site index curves could be detected on remeasured sample plots. Existing growth models were applied on measured stand data to predict future growth. Calculated mean annual increments were on average 20% lower when site index was predicted by site factors instead of height and age. The bias means that the yield capacity of planted spruce in northern Sweden has been underestimated by about 35%.     

Comparison of beta,Johnson’s SB,Weibull and truncated Weibull functions for modeling the diameter distribution of forest stands in Catalonia (north-east of Spain)

The purpose of this study was to compare the beta, Johnson SB, Weibull and truncated Weibull functions in describing the diameter distributions of forest stands in Catalonia. The data consisted of permanent sample plots from the Spanish National Forest Inventory in Catalonia. The empirical data represent left-truncated distributions, as the smallest diameter measured in the field was 7.5 cm. A total of 1,242 plots were used to fit the functions and analyze their performance. The distribution functions were fitted to the diameter distributions of the number of stems (DD_N) and stand basal area (DD_G). The performance of the candidate functions was compared by means of their bias and RMSE for different diameter sums measuring the difference between the empirical and fitted distributions. The leftmost part (from 0 to 7.5 cm) of the non-truncated functions was ignored in this analysis. Sensitivity analyses were conducted to check whether the results depended on the number of trees measured in the stand, or the main species of the stand. The truncated Weibull function for the diameter distribution of stand basal area appeared to be in all cases the most accurate and consistent function. Generally, functions describing the distribution of stand basal area performed better than functions that described the distribution of the number of trees. Of the basal area distributions, beta and Johnson’s SB were the second best and nearly equally good with each other. The order of precision of the tested functions was: truncated Weibull for DD_G, truncated Weibull for DD_N, Johnson’s SB for DD_G, beta for DD_G, beta for DD_N and Weibull for DD_G, Weibull for DD_N, and Johnson’s SB for DD_N.     

桂西南马尾松人工林生长对不同强度采伐的动态响应

[目的]研究不同强度采伐下马尾松的生长动态,筛选适宜的采伐强度,为马尾松人工林近自然经营提供技术支撑。[方法]2007年10月在14年生马尾松人工林(保存密度1 100株·hm~(-2))内进行采伐试验,设置4个采伐强度,即保留密度分别为225、300、375、450株·hm~(-2),以不采伐为对照;其后,自2008年开始连续8 a,每2 a测定1次马尾松的胸径、树高、枝下高和冠面积等生长指标,并计算单株材积和林分蓄积量,应用方差分析和Duncan多重比较分析生长指标对不同采伐强度的动态响应。[结果]表明:采伐强度显著影响林分生长,其中,林分平均胸径、单株材积、冠面积的年均增长量随保留密度增大而减小,但均显著高于对照(P0.05)。采伐后第1 3年,马尾松冠面积增长量显著高于采伐后期,胸径则在采伐后第3 5年最高,而不同采伐强度对林分树高生长影响不明显。保留密度显著影响林分枝下高和蓄积量的动态变化,其年均增长量随密度增大而递增。5个处理间林分蓄积年均增长量的差异随林龄的增大而逐渐缩小。[结论]马尾松人工林生长对不同强度采伐的动态响应以树冠最敏感,冠面积首先陡然增大,进而引起胸径的快速生长。树高和枝下高在采伐后年均增量变化相对平稳。4个采伐强度均显著促进单株材积生长,而仅保留密度为225株·hm~(-2)的采伐对林分蓄积增长量影响显著。综合比较林分的单株材积和林分蓄积连年增长量,建议在桂西南15年生马尾松人工林近自然经营中宜选择300株·hm~(-2)的保留密度进行采伐。