A stem taper equation for Eucalyptus camaldulensis in northeast Nigeria

Abstract

Eucalyptus (Eucalyptus camaldulensis Dehnh.) is an important exotic species in northern Nigeria commonly used for poles and timber. Sustainable management of this resource would require quantifying its volume. Stem taper equations are one of the main and most efficient methods for estimating stem volume to any merchantable limit of a species. There is currently no taper equation for Eucalyptus species in Nigeria. Therefore, this study developed taper equations for E. camaldulensis in northern Nigeria. Data for this study were obtained from a private plantation in Jalingo Local Government Area, Taraba State, Nigeria. 68 trees were felled and sectioned into 1-m bolt across the stem to a merchantable limit of 5 cm, which were used as the fitting dataset. An additional 22 trees were felled and used to validate the taper equations for stem volume estimation. Seven taper equations were initially fitted to the dataset using nonlinear least squares. The best taper equation was then refitted using a nonlinear mixed-effects approach and calibrated using diameters of one to five sections from the butt end. The taper equations were numerically integrated to obtain the stem volume, which was compared with empirical volume equations. The result shows that the Kozak (Can J For Res 27(5): 619–629. 10.1139/x97-011, 1997) equation, which included eight parameters, provided the best fit for predicting section diameters for under and over bark. The mixed-effects taper equation (NLME-TE) explained most stem diameter variations in the fitting dataset (pseudo-R²: 0.986–0.987; RMSE: 0.547–0.578 cm) without substantial residual trends. The validation showed that the prediction accuracy of the integrated NLME-TE improved as the number of sectional diameter measurements increased, with at least a 35% reduction in volume estimate error. For practical implementation, two calibration sectional diameter measurements taken from the butt end per tree are recommended. This approach would reduce measurement effort and cost while improving model performance.

Keywords

Eucalyptus plantation / Stem volume / Nonlinear least squares / Mixed-effects / Response calibration

References

[1] Abubakar B, Shamaki SB, Usman U, Bello AG, Dantani A, Gada ZY, Abubakar A. Diameter growth and mortality assessment of river red gum (Eucalyptus camaldulensis (Dehnh.)) stands in Afaka forest reserve, Kaduna-Nigeria. Int J Agric Plant Sci, 2019, 1(2): 24-29

[2] Adamec Z, Adolt R, Drápela K, Závodský J. Evaluation of different calibration approaches for merchantable volume predictions of Norway spruce using nonlinear mixed effects model. Forests, 2019, 10(12): 1104

[3] Adegbehin JO, Okojie JA, Nokoe S. Growth of Eucalyptus tereticornis and Pinus caribaea at different sites in Northern Nigeria. For Ecol Manage, 1988, 23(4261-272

[4] Alo AA, Onyekwelu JC, Akindele SO. Taper equations for Gmelina arborea in Omo forest reserve, southwestern Nigeria. J Appl Trop Agric, 2011, 16(1–2): 120-127

[5] Arias-Rodil M, Castedo-Dorado F, Cámara-Obregón A, Diéguez-Aranda U. Fitting and calibrating a multilevel mixed-effects stem taper model for maritime pine in NW Spain. PLoS ONE, 2015, 10(12 e0143521

[6] Bi HQ. Trigonometric variable-form taper equations for Australian eucalypts. For Sci, 2000, 46(3): 397-409

[7] Boczniewicz D, Mason EG, Morgenroth JA. Developing fully compatible taper and volume equations for all stem components of Eucalyptus globoidea Blakely trees in New Zealand. NZ J for Sci, 2022, 52(6): 1-14

[8] Bronisz K, Mehtätalo L. Mixed-effects generalized height–diameter model for young silver birch stands on post-agricultural lands. For Ecol Manage, 2020, 460 117901

[9] Bronisz K, Zasada M. Taper models for black locust in west Poland. Silva Fenn, 2020, 54(5 10351

[10] Burkhart HE, Tomé M. Modeling forest trees and stands, 2012, Dordrecht, Springer Science & Business Media458

[11] Cao QV, Wang J. Evaluation of methods for calibrating a tree taper equation. For Sci, 2015, 61(2): 213-219

[12] Corral-Rivas J, Vega-Nieva D, Rodríguez-Soalleiro R, López-Sánchez C, Wehenkel C, Vargas-Larreta B, Álvarez-González J, Ruiz-González A. Compatible system for predicting total and merchantable stem volume over and under bark, branch volume and whole-tree volume of pine species. Forests, 2017, 8(11): 417

[13] Da Silva LMS, Rodriguez LCE, Caixeta Filho JV, Bauch SC. Fitting a taper function to minimize the sum of absolute deviations. Sci Agric (Piracicaba, Braz), 2006, 63(5460-470

[14] Dávila-Molina DE, Sáenz-Romero C, Aguirre-Calderón OA, Lopez-Toledo L. Age contributes to volume estimation and form factor of Pinus pseudostrobus Lindley in commercial forest plantations from western Mexico. J Sustain for, 2022, 42(3): 336-351

[15] Demaerschalk JP. Converting volume equations to compatible taper equations. For Sci, 1972, 18(3): 241-245

[16] Demaerschalk JP, Kozak A. The whole-Bole system: a conditioned dual-equation system for precise prediction of tree profiles. Can J for Res, 1977, 7(3): 488-497

[17] de-Miguel S, Mehtätalo L, Shater Z, Kraid B, Pukkala T. Evaluating marginal and conditional predictions of taper models in the absence of calibration data. Can J for Res, 2012, 42(7): 1383-1394

[18] Ekhuemelo D, Onah G, Wuam L. Evaluation of the uses of Eucalyptus species in makurdi local government area of Benue state Nigeria. GSC Biol Pharm Sci, 2017, 1(1): 25-34

[19] Elzhov TV, Mullen KM, Spiess A, Bolker B (2023) minpack.lm: R Interface to the Levenberg-Marquardt Nonlinear Least-Squares Algorithm Found in MINPACK, Plus Support for Bounds. R package version 1.2–3. https://CRAN.R-project.org/package=minpack.lm

[20] Fonweban J. An evaluation of numerical integration of taper functions for volume estimation in Eucalyptus saligna stands. J Trop for Sci, 1999, 11(2): 410-419

[21] Hada S. On the taper of the Sugi (Cryptomeria japonica. D. DON.) bole by the form quotient. J Jpn for Soc, 1958, 40(9): 379-382

[22] Heiðarsson L, Pukkala T. Taper functions for lodgepole pine (Pinus contorta) and Siberian larch (Larix sibirica) in Iceland. Icel Agric Sci, 2011, 24(2011): 3-11

[23] Jaiyeoba IA. Soil rehabilitation through afforestation: evaluation of the performance of Eucalyptus and pine plantations in a Nigerian savanna environment. Land Degrad Dev, 2001, 12(2): 183-194

[24] Kozak A. A variable-exponent taper equation. Can J for Res, 1988, 18(11): 1363-1368

[25] Kozak A. Effects of multicollinearity and autocorrelation on the variable-exponent taper functions. Can J for Res, 1997, 27(5): 619-629

[26] Kozak A. My last words on taper equations. Fortschr Chron, 2004, 80(4): 507-515

[27] Kozak A, Munro DD, Smith JHG. Taper functions and their application in forest inventory. For Chron, 1969, 45(4): 278-283

[28] Kozak A, Smith JHG. Standards for evaluating taper estimating systems. Fortschr Chron, 1993, 69(4): 438-444

[29] Krisnawati H. A compatible estimation model of stem volume and taper for Acacia mangium Willd. plantations. Ina J for Res, 2016, 3(1): 49-64

[30] Kublin E, Augustin NH, Lappi J. A flexible regression model for diameter prediction. Eur J for Res, 2008, 127(5): 415-428

[31] Kurowska EE, Czerniak A, Garba ML. Afforestation of transformed savanna and resulting land cover change: a case study of Zaria (Nigeria). Sustainability, 2022, 14(3): 1160

[32] Lanssanova LR, et al.. Mixed-effect non-linear modelling for diameter estimation along the stem of Tectona grandis in mid-western Brazil. South A J Sci, 2019, 81(2167-217

[33] Lee WK, Seo JH, Son YM, Lee KH, von Gadow K. Modeling stem profiles for Pinus densiflora in Korea. For Ecol Manage, 2003, 172(1): 69-77

[34] Li DD, Guo HT, Jia WW, Wang F. Analysis of taper functions for Larix olgensis using mixed models and TLS. Forests, 2021, 12(2): 196

[35] Maraseni T, Cockfield G, Apan A. Estimation of taper rates and volume of smaller-sized logs in spotted gum saw timber plantations in Southeast Queensland, Australia. South Hemisph J, 2007, 69(3169-173

[36] Max T, Burkhart H. Segmented polynomial regression applied to taper equations. For Sci, 1976, 22: 283-289

[37] McTague JP, Weiskittel A. Evolution, history, and use of stem taper equations: a review of their development, application, and implementation. Can J for Res, 2021, 51(2): 210-235

[38] Mehtätalo L, Lappi J (2020) Biometry for forestry and environmental data: with examples in R (1st ed), Chapman and Hall/CRC, New York, p 426. https://doi.org/10.1201/9780429173462

[39] Morley T, Little K. Comparison of taper functions between two planted and coppiced eucalypt clonal hybrids, South Africa. New for, 2012, 43(2): 129-141

[40] Morrone S, Green PC. Regional differences in stem form between southern and northern red spruce (Picea rubens Sarg.) populations. Forestry, 2024, 97(5): 771-784

[41] Muhairwe CK. Tree form and taper variation over time for interior lodgepole pine. Can J for Res, 1994, 24(9): 1904-1913

[42] Muhairwe CK. Taper equations for Eucalyptus pilularis and Eucalyptus grandis for the north coast in New South Wales, Australia. For Ecol Manage, 1999, 113(2–3): 251-269

[43] Newberry JD, Burkhart HE. Variable-form stem profile models for loblolly pine. Can J for Res, 1986, 16(1): 109-114

[44] Newnham RM. Variable-form taper functions for four Alberta tree species. Can J for Res, 1992, 22(2): 210-223

[45] Ogana FN, Holmström E, Sharma RP, Langvall O, Nilsson U. Optimizing height measurement for the long-term forest experiments in Sweden. For Ecol Manage, 2023, 532 120843

[46] Ogana FN, Corral-Rivas S (2023) EMAR: Empirical Model Assessment. R package version 1.0.0. https://CRAN.R-project.org/package=EMAR.

[47] Oluwajuwon TV, Hossain MS, Ludjen AN, Hasert DP, Sitanggang M, Israel R, Offiah CJ, Alawode GL, Ogana FN. Form factors and volume models for Falcataria moluccana in smallholder plantations, Central Kalimantan, Indonesia. Aust, 2024, 87(3): 101-113

[48] Oluwajuwon TV, Ogbuka CE, Ogana FN, Hossain MS, Israel R, Lee DJ. Describing and modelling stem form of tropical tree species with form factor: a comprehensive review. Forests, 2025, 16(1): 29

[49] Ormerod DW. A simple Bole model. Fortschr Chron, 1973, 49(3136-138

[50] Otegbeye GO, Samarawira I. Growth and form of Eucalyptus camaldulensis Dehnh. provenances in northern Nigeria. For Ecol Manage, 1991, 42(3–4): 219-228

[51] Pinheiro JC, Bates DM. Mixed-effects models in S and S-PLUS, 2000, New York, Springer-Verlag528

[52] Pinheiro JC, Bates DM (2023) nlme: Linear and nonlinear mixed effects models. R package version 3.1–162. https://CRAN.R-project.org/package=nlme

[53] Poudel KP, Cao QV. Evaluation of methods to predict weibull parameters for characterizing diameter distributions. For Sci, 2013, 59(2): 243-252

[54] PTPU (1988) Agroforestry manual. Nigeria Second Forestry Project (Afforestation Prog.). Plateau State Project Unit (PTPU), Jos, pp 5–20

[55] Pukkala T, Hanssen K, Andreassen K. Stem taper and bark functions for Norway spruce in Norway. Silva Fenn, 2019

[56] R Core Team (2023) R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna, Austria

[57] Salekin S, Catalán CH, Boczniewicz D, Phiri D, Morgenroth J, Meason DF, Mason EG. Global tree taper modelling: a review of applications, methods, functions, and their parameters. Forests, 2021, 12(7): 913

[58] Saud P, Chapagain TR, Bhandari SK, Moser WK. Taper functions to predict the upper stem diameter of Chir pine (Pinus roxburghii) in the mid-hills of Nepal. Trees for People, 2024, 17 100627

[59] Subedi N, Sharma M, Parton J. Effects of sample size and tree selection criteria on the performance of taper equations. Scand J Res, 2011, 26(6555-567

[60] Tarp-Johansen MJ, Skovsgaard JP, Madsen SF, Johannsen VK, Skovgaard I. Compatible stem taper and stem volume functions for oak (Quercus robur L and Q petraea (Matt) Liebl) in Denmark. Ann for Sci, 1997, 54(7): 577-595

[61] Vonesh E, Chinchilli VM (1997) Linear and nonlinear models for the analysis of repeated measurements (1st ed), CRC Press, Boca Raton, p 560. https://doi.org/10.1201/9781482293272

[62] Warner AJ, Jamroenprucksa M, Puangchit L. Development and evaluation of teak (Tectona grandis L.f.) taper equations in northern Thailand. Agric Nat Resour, 2016, 50(5): 362-367

[63] Watanabe Y, Masunaga T, Fashola OO, Agboola A, Oviasuyi PK, Wakatsuki T. Eucalyptus camaldulensis and Pinus caribaea growth in relation to soil physico-chemical properties in plantation forests in Northern Nigeria. Soil Sci Plant Nutr, 2009, 55(1): 132-141

[64] Wimbush SH. Afforestation of restored tin-mining land in Nigeria. Commonw for Rev, 1963, 42(3): 255-262

[65] Zhang SS, Sun JJ, Duan AG, Zhang JG. Variable-exponent taper equation based on multilevel nonlinear mixed effect for Chinese fir in China. Forests, 2021, 12(2): 126

[66] Zheng C, Wang Y, Jia L, Mason EG, We S, Sun C, Duan J. Compatible taper-volume models of Quercus variabilis Blume forests in North China. Iforest, 2017, 10(3567-575

RIGHTS & PERMISSIONS

This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply