Skip to main content
Have a personal or library account? Click to login
Height–diameter relationship of plantation-grown juvenile black locust trees is differentiated according to their growth rate, which is positively affected by spacing Cover

Height–diameter relationship of plantation-grown juvenile black locust trees is differentiated according to their growth rate, which is positively affected by spacing

Open Access
|Dec 2022

Figures & Tables

Table 1

Description of the experimental data used to derive the height–diameter relationships for Robinia pseudoacacia L.

Sample plotParental genotypeGeographic regionPlot size (m2)Growing space (m2)Age (years)Number of measured treesd (cm)h (m)
PP 1southern185.002.56483.6–7.14.7–6.3
PP 2southern173.252.52513.1–6.04.7–6.5
PP 3southern202.502.53503.8–4.66.4–6.5
PP 4southern165.003.54671.5–4.62.9–4.7
PP 5Roszin Varganorthern304.004.55642.3–6.72.8–7.7
PP 6Karaisennorthern306.044.55522.4–7.63.9–9.2
PP 7Tsarevetsnorthern300.004.55702.2–9.14.2–9.2

[i] Abbreviations: d – breast height diameter of the tree (cm), h – total tree height (m). Data range (minimum–maximum) of d and h is shown.

Figure 1

Mixed-effects height–diameter model estimated by plots

Figure 2

Generalised height–diameter model estimated by spacings

Table 2

Local height–diameter models and goodness-of-fit statistics

ModelMean Bias (m)Aggregated MSE (m)*90th percentile of ARBi*
abbreviationformula
M1h = 1.3 + a1d0.15156.90680.2184
M2h = 1.3 + a1d + a2d2–0.00243.84960.1367
M3h = a0da1–0.00314.09650.1176
M4h = 1.3 + exp[a0 + a1/(d + 1)]–0.00343.86160.1100
M5h = 1.3 + a0 exp(a1/d)–0.00053.77900.1077
M6h = 1.3 + a0da1–0.00534.17140.1191
M7h = 1.3 + a0d/(a1 + d)0.12216.92140.2412

Abbreviations: h – tree height (m), d – diameter at breast height (cm), a0, a1, a2 – model parameters.

* Formulae: mean bias – i=1PMbiasPM ; Aggregated MSE – mean squared errors aggregated from all plots – i=1PM(biasi2+vari) ; 90th percentile of the absolute values of the relative biases ARBi, where biasi is the mean error per plot biasi=j=1N(hj^hj)N , vari is the error variance per plot, vari=j=1N(hj^hj)2N , ARBi is the absolute value of the relative bias per plot ARBi=|j=1N(hj^hj)hjN| , PM – total number of plots, N – number of the measured tree heights per plot, hj – observed tree height, hj^ – predicted tree height.

Table 3

Comparison of the mixed-effects height–diameter models by Likelihood Ratio Test

ModelMixed-effects parametersdflogLikTestLikelihood Ratiop-value
abbreviationlocalisation level
1GR.PP.01sample plot within geographic regionb0, b19–413.990
2PP.01sample plotb0, b16–414.3241 vs. 20.6690.881
3PP.0sample plotb04–419.8082 vs. 310.9670.004
4GR.PP.0sample plot within geographic regionb05–419.4073 vs. 40.8000.371
5PP.1sample plotb14–430.3274 vs. 521.839<0.001
6GR.PP.1sample plot within geographic regionb15–430.1305 vs. 60.3930.531
7GR.01geographic regionb0, b16–438.4716 vs. 716.682<0.001
8GR.0geographic regionb04–440.7487 vs. 84.5530.103
9GR.1geographic regionb14–445.830
10simple localno localisation3–453.0599 vs 1014.457<0.001

[i] Abbreviations: df – degrees of freedom; logLik – log-likelihood.

Table 4

Regression estimates and goodness-of-fit statistics of the mixed-effects and the generalised deterministic height–diameter models

Mixed-effects model: h = 1.3 + (u + b0)eb1
ME 0.773Model parametersModel variances
b0b1 σu2 σ2
MARE% 10.322Estimate8.302–2.6310.6390.555
Bias 0.010SE0.3950.1460.3630.042
Generalised deterministic model: h = 1.3 + (a0 + a1Spacing)eb1/d
ME 0.724Variance functionModel parameters
(θ+dη)2a0a1b1
MARE% 11.634Variance function parametersEstimate7.3220.422–2.918
Bias 0.025ΘHSE0.4170.0870.123
37.2431.701

[i] Abbreviations: ME – model efficiency, ME=1(hi,hi^)2(hi,h¯)2 , MARE% is the average of the absolute values of the relative errors ARE%, ARE%=|hi,hi^|hi100 Bias is the average error, Bias=(hi,hi^)n , where hi, hi^ are experimental and predicted height values of the i-th measurement and represents the mean observed height value; d – diameter at breast height (cm), h – total tree height (m), Spacing – growing space per plant (m2), a0, a1, b0, b1 – fixed model parameters, θ, η – variance function parameters, u – random parameter component, σu2 – variance of the random component, σ2 – residual variance of the model, SE – standard error.

Table 5

Calibration of the random component of the mixed-effects model with 1-tree height–diameter measurement: accuracy estimates by plots

PlotuRelative errorsMARE%Bias
Perc0Perc25Perc50Perc75Perc100
PP 1–0.789–0.719–0.150–0.077–0.0030.16611.829–0.363
PP 2–0.561–0.579–0.133–0.0370.0160.20511.001–0.233
PP 30.369–0.100–0.0180.0170.0760.1535.4230.163
PP 4–0.981–0.350–0.182–0.0670.0630.39114.612–0.228
PP 5–0.125–0.405–0.116–0.0110.0860.19810.536–0.050
PP 61.118–0.390–0.0470.0860.1670.24213.0040.506
PP 70.968–0.673–0.0250.0680.1180.23812.6610.403

[i] Abbreviations: u – estimated random component of the rate parameter b0; Perc0, Perc25, Perc50, Perc75, Perc100 – 0th, 25th, 50th, 75th and 100th percentiles of the relative errors estimated as: hihi^hi ; MARE% is the average of the absolute values of the relative errors ARE%, ARE%=|hihi^|hi100 ; Bias is the average error, Bias=(hihi^)n , where hi, hi^ are experimental and predicted height values of the i-th measurement and represents the mean observed height value.

DOI: https://doi.org/10.2478/ffp-2022-0019 | Journal eISSN: 2199-5907 | Journal ISSN: 0071-6677
Language: English
Page range: 195 - 205
Submitted on: Aug 4, 2022
Accepted on: Sep 13, 2022
Published on: Dec 10, 2022
Published by: Forest Research Institute
In partnership with: Paradigm Publishing Services
Related subjects:

© 2022 Tatiana Stankova, Proletka Dimitrova, Veselka Gyuleva, Penka Stefanova, published by Forest Research Institute
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.