- Research
- Open access
- Published:

# Evaluating diameter distribution series of small-leaved lime (*Tilia cordata* Mill.) in forest stands

*Plant Methods*
**volume 17**, Article number: 45 (2021)

## Abstract

### Background

The paper provides studies on the structure of planted small-leaved lime (*Tilia cordata* Mill.) in the conditions of the Bashkir Cis-Urals. This work aimed to analyze their assortment by diameters and compile appropriate tables. This is the first study of lime for this region. The results of the study are based on data from 69 temporary sample plots. Stands are represented by trees of 11–79 years old, not affected by thinning. They belong to the I–III growth classes and the most common goutweed forests.

### Results

It was found, that small-leaved limes have a specific structure—coefficient of diameter variation in stands decreases with a higher average diameter, reaching 26 and 16% at 28 cm with thin and dense initial planting, respectively. The variability of tree diameters is related to the average stand diameter and is due to the initial density of the grown plantations. Correlations between the coefficients of asymmetry and excess with age and the average diameter was revealed.

### Conclusions

The verification of theoretically calculated frequencies of distribution series to empirically observed frequencies showed a discrepancy in 29 and 19% of the total number, respectively. For the remaining series of distributions, the scale and shape parameters of the Weibull function were modeled using the average, minimum and maximum diameters, standard deviation, coefficients of asymmetry and excess. This made it possible to develop stand tables for small-leaved limes depending on the average diameter.

## Introduction

The application of forest stand structure and development regularities is an integral part of forest management planning, decision-making, scientific research. It leads to higher forest quality and productivity, determines economic efficiency. The structure of stands affects the intensity of biogeocenotic processes, the efficiency of production and deposition of organic matter, the stability and biosphere functions of forest ecosystems. The stand structure is understood as a cumulative combination of examination indicators variation, the distribution of trees by these values, their relationship at certain age stages of the forest community development.

When studying the structure of stands, the analysis of trunk diameter variation is of the greatest interest. The diameter at breast height (d_{1,3}) is an important indicator of a tree. It is widely used in forest inventory to calculate some stand parameters such as stem volume and aboveground biomass that cannot be measured directly. The d_{1,3} distribution is an indicator of the forest structure [1], its stock, assortment composition, and so on [2,3,4].

There is continuing interest in studying the structure of stands, despite more than a century of history [5,6,7]. The stand structure is analyzed in terms of different tree species [1, 8, 9], composition [10,11,12], average age [4, 13], age structure [10, 14], growth site conditions [3, 15], density [4], factors that affect the tree number distribution by diameter class [3, 6, 16]. Bassil et al. [17] studied the temporal stability of the number of trees by diameter classes in different-aged Northern hardwoods that are subject to logging [17]. Many publications are aimed at searching for models and evaluating their effectiveness in describing the d_{1,3} distribution [1, 3, 18], including those subjected to thinning [13]. Thus, Safonov et al. [12] describe the functionality of the Forest-Fit software package and its features that simplify the estimation of common probability distributions for modeling d_{1,3} distributions, using the example analysis of mixed *Pinus ponderosa* forest compartments. Paradis and Lebel [19] propose a methodology that can be used to build a d_{1,3} distribution model for any combination of forest species and types in Quebec (Canada), using easily accessible data from the state program for the inventory of sample areas.

Though there is a significant amount of literature on distribution series of different tree species and natural conditions, studies on small-leaved limes (*Tilia cordata* Mill.) are not numerous [10, 15, 20].

Much attention is paid to determining the structure of planted stands [16, 21, 22], because even the same species, due to differences in the initial density, tree distribution, tree ages, the quality of planting material and different growing conditions, will form a variety of structures.

Pure and mixed artificial stands of small-leaved limes have been planted in the Bashkir Cis-Urals since the end of the 30 s of the last century. There have been no studies describing their structure and dynamics. Since forest stand structure regularities by estimation indicators are statistical in nature, there is a need to develop differentiated stand inventory standards [2]. Therefore, the goal of the work is to consider diameter distribution characteristics in planted stands of small-leaved lime and to compile appropriate tables for the conditions of the Bashkir Cis-Urals.

## Materials and methods

The paper studies the diameter structure of 11–79-year-old small-leaved lime plantations in 69 temporary sample plots located in the Bashkir Cis-Urals and belonging to the forest-steppe zone of the South Ural forest-steppe region of the European part of the Russian Federation. According to the S.F. Kurnaev’s fractional forest zoning, this region belongs to the broad-leaved forest zone of the forest-steppe subzone within the Russian plain [23].

The studied plots (except for three sample areas) are located on the territory of Ufa city and the Ufa municipal district, located at 54°70′ north longitude and 55°90′ east longitude 150 m above sea level (Fig. 1). The climate of the district is continental, rather humid. The average annual air temperature is 3.0 °C, the average temperature in January is − 14.5 °C, in July 19.5 °C with an absolute maximum of 40 °C and an absolute minimum of − 50 °C. The average annual precipitation is within 500–600 mm. It is about 350 mm during the growing season [24]. In these conditions, small-leaved limes develop according to the I–III growth classes and represent productive phytocenoses.

The data was collected from rectangular sample plots of 0.1 ha or more, depending on the stand density. Each plot is a homogeneous plantation. All trees > 3.9 cm were d_{1,3} estimated by one-centimeter diameter class in every plot. Tree heights were measured. The age is determined based on the data of mensurational descriptions and the list of forest plantations. The remaining stand estimation indicators were calculated based on the trees counted [2].

The studied plantations were of different densities, and were not thinned. They are grouped for an analysis as follows: thin—from 2 to 6 thousand pieces per ha (52 plots), medium density—from 6 to 11 thousand pieces per ha (10 plots), dense and very dense—from 11 to 20 thousand pieces per ha (7 plots). The main dendrometric characteristics of the selected plantings are shown in Table 1.

Basic statistics are calculated for all series of d_{1,3} distributions: sample size (n), arithmetic mean (X), standard deviation (S), coefficients of variation (C_{v}), asymmetry (A_{s}), excess (E_{x}), and others. By conducting a correlation analysis, the relationships between evaluation indicators and statistics of the d_{1,3} distribution series of trees were revealed, being approximated by the following functions:

where y_{i} is the minimum (*d*_{min}) and maximum (*d*_{max}), arithmetic mean (X) values of the diameters of the distribution series (in cm); b_{i} is the coefficients of models; D is the average inventory diameter of the stand, cm.

Tree diameter distribution series statistics (X, d_{min} and d_{max}, S, A_{s} and E_{x}) was compared by calculating the mean-square deviation of pairwise matched models (σ, %) [25]:

where a_{i} and C_{i} are pairwise compared data of distribution statistics calculated from models 2–4; n is the number of pairs compared, pcs.

To describe these distribution series, the Pearson family curves of I–VII types were used [10]. The number of trunks was equalized by one-centimeter diameter classes in Microsoft Excel. The correspondence of theoretically calculated frequencies of distribution series to empirically observed frequencies was estimated by the B. S. Iastremskii criterion [26].

Since there was no special software for leveling the curves of the Pearson family, a two-parameter Weibull distribution in the Statistica environment was used to calculate the expected values of d_{1,3} frequencies with an accuracy of up to a centimeter, expressed as a percentage of the total number of trunks of each sample plot. It has the following density function (for positive parameters γ and β):

where γ is the parameter of the distribution form; β is the distribution scale parameter; *e* is the Euler number (2.71…).

Estimation accuracy of the Weibull function was verified using the Kolmogorov–Smirnov test (K–S) and the Anderson–Darlig test (AD) with a probability of 5%. After evaluating the γ and β parameters of this function, a correlation and regression analysis was performed and multiple regression models were obtained to predict them depending on the distribution statistics (X, d_{min} and d_{max}, S, A_{s} and E_{x}).

## Results

The tree diameter is the main examination indicator for studying the structure of the stand. It varies from 1 to 47 cm by sample plots. The average arithmetic value of d_{1,3} for individual objects ranges from 2.6 cm to 25.4 cm. The coefficient of variation, which is an indicator of the homogeneous structure of the stand, varies from 5 to 56%. In three sample plots, the d_{1,3} distribution is characterized by left—sided asymmetry (from − 0.23 to − 0.01), while it is right-sided in the rest. it is insignificant in five plots (from 0.01 to 0.09). The coefficient of asymmetry in other compartments ranges from 0.11 to 1.93. D_{1,3} distributions are characterized by different density grades from − 0.97 to 5.61. Analysis of the asymmetry and excess coefficients for 34 and 13 sample plots, respectively, showed that they are beyond their twofold basic errors. This indicates the difference between the given series from the standard one. For the rest of the sample areas, the asymmetry and excess confidence estimation by the t-Student’s test, on the contrary, shows that there is no deviation of the distribution curves from the standard one (t_{0, 05} < 1.97).

To find out the regularities of the diameter distribution dynamics and its statistics, a correlation analysis was performed (Table 2).

The studied limes showed the following changes in the distribution series with tree aging. The relation of the asymmetry and excess coefficients with age is not expressed, but there is a tendency to decrease the asymmetry coefficient and increase the excess coefficient with increasing age. The coefficient of variation significantly decreases with the age of plantations. The standard deviation increases with aging and higher average diameters of thin and dense plantations. The value of the average tree rank for thin planting significantly correlates with the coefficient of asymmetry (r = − 0.624) and the minimum diameter of the tree in the distribution series (r = 0.717). The identified relationships made it possible to create models (1–3), presented in Table 3.

The data in Table 3 indicate the significance of the obtained equations (except for two cases) and the possibility of their application for forecasting.

Checking the statistics of the distribution series according to formula 4 using the developed models (1–3, Table 3) showed that the initial density of plantations (within the considered planting schemes) affects them: the degree of difference is significant—for *d*_{min} (27–101%), *d*_{max} (16–19%) and S (18–21%). A similar test for the coefficients of asymmetry and excess demonstrated several times greater differences: up to 535% for A_{s}, up to 1817% for E_{x}. Differences in the distribution series depending on the density of forest plantation can be visually estimated from Fig. 2, which illustrates aligned indicators of variability in tree diameters.

The given data indicate that diameters in thin plantations vary greatly compared to dense and very dense stands. Meanwhile, medium dense plantations have an intermediate position of variability. With average diameters of 8–14 cm, the variability is closer to dense and very dense stands, and with large diameters (22–28 cm)—to thin. With higher average stand diameter, the variability of trees in thickness decreases, reaching values of 26 and 16% at 28 cm with thin and dense initial planting, respectively.

The use of Pearson curves to approximate the tree distribution series by diameter showed that the type I function provides more satisfactory description for 46 of the 69 sample plots, types II and VII for 7 plots, type IV for 6 sample compartments, type VI for the rest. This indicates the heterogeneity of the tree distribution by their size with the age of stands. The verification of theoretically calculated frequencies of distribution series to empirically observed frequencies by the B.S. Iastremskii criterion proved to be relevant for 49 applied models (71%). Figure 3 presents the simulation results for six sample plots with the same initial planting density (thin, row spacing 3.0 m, plant spacing 0.7 m).

It shows that trees of the smaller diameter classes are more frequent in young stands and their number decreases with age. At the same time, as the forest grows, the frequency of large diameter trees increases, shifting the curve to the right. With a higher average diameter of the stand, the distribution curves stretch, they become less steep and less elongated upwards in the center, approach the abscissas axis more slowly when moving away from the middle, and become the closest to the standard distribution.

Due to the small number of samples of medium, dense and very dense plantations, insignificant correlations between the coefficients of asymmetry and excess with the average diameter of stands (Table 2), as well as the deviation of two models (Table 3), the following calculations were performed only for stands with thin planting density. The distribution series by one-centimeter diameter classes as a percentage were estimated by the two-parameter Weibull function. The model validity test evaluated by the K–S and AD criteria proved their adequacy for 42 sample plots, since the values of these tests were not significant (p > 0.05). This condition indicates that the observed and expected frequencies are not statistically different from each other, and therefore the null hypothesis cannot be rejected. On the other hand, for 10 sample areas (19%), significant values of the K–S and AD tests indicate inadequate estimates of diameter frequencies (p < 0.05).

A preliminary linear correlation analysis was performed to assess trends between the two parameters of the Weibull function and the distribution statistics (X, d_{min}, d_{max}, S, A_{s}, and E_{x}) (Table 4).

The data in Table 4 indicates a high chance of successfully projecting the probability distribution by restoring the coefficients of the Weibull function by multiple regression. The following models were offered:

Both models presented have high F scores (p < 0.0001) and explain 99 and 96% of the variability, respectively.

The coefficients obtained for the Weibull function were used to estimate the evolution of the diameter distribution in planted stands of small-leaved lime depending on the class of the average diameter (one-centimeter diameter class), the calculated values of which are given in Appendix.

## Discussion

Planted stands of small-leaved limes in distribution series by d_{1,3}, being the main forest inventory indicator, change the place of the average tree, the reduction numbers by rank, the minimum and maximum diameters, the concentration of the trunk number percentage in the diameter class, the standard deviation value, asymmetry and excess of the series, which is consistent with other studies [2, 4, 6].

The use of Pearson curves made it possible to effectively describe the d_{1,3} distribution series. This is consistent with the view of Bachioua [27], who considers that they can flexibly adapt to a set of well-known theoretical and practical distributions. The work of Shakil et al. [28] showed a better result of the proposed Pearson model than the gamma, lognormal and inverse Gaussian distributions.

The discrepancy between the theoretically calculated frequencies of the distribution series and the empirically observed frequencies of the Pearson curves and the two-parameter Weibull function (29 and 19% of the total number, respectively) indicates the complex nonlinear nature of the d_{1,3} distributions and modeling the probability-density function, which conforms to the research of Diamantopoulou et al. [29] and Pach and Podlaski [1].

Many authors report that the Weibull function is widely used in many forestry applications to model d_{1,3} distributions [3, 13, 22], because it has a simple cumulative function and is flexible in selecting distributions for different forms and degrees of asymmetry [11, 29, 30]. Moreover, the simple two-parameter form of its distribution often gives better results compared to the three-parameter form [3, 6, 30]. In addition, the two-parameter Weibull function was used in a large number of related studies due to the high correlation of its parameters (γ and β) with the stand characteristics [13, 17, 18]. Therefore, its application is suitable for describing the d_{1,3} distributions in this study.

The correlation analysis results obtained between the two parameters (γ and β) of the Weibull function and the characteristics of *d*_{min} and *d*_{max} are similar to those observed by Binoti et al. [13]. Except that there is no connection between *d*_{max} and γ. Linear models that link the parameters β and γ to the main distribution statistics show very high determination coefficients (99 and 96%, respectively). It is noted that their performance deteriorates as the average value of d_{1,3}, asymmetry and excess of the distribution series increases, which is consistent with studies by Lima et al. [3].

Diameter distribution models can be useful in describing and analyzing the stand structure, age distribution, growth and yields. They can be employed in assessing the stand stability and calculating a number of trees in each diameter class to plan forest management activities, reveal previous violations, predict successions of forests and land biomass reserves, etc. [3, 4, 12]. Tree diameter distribution modelling based on remote sensing by laser scanning also requires indirect estimation of these distributions using forecasting models [22]. Determining stand diameter distribution is costly due to measuring diameters for a large number of trees during inventory. Costs can be reduced by using diameter distribution models based on diameter classes, depending on the average diameter, the number of trunks and the stand basal area [22, 30]. However, the given models should be treated with caution, since these stand parameters are used in many different cutting strategies resulting in different diameter distributions [6]. It becomes particularly evident for forest crops of different densities, especially during the stand formation, since there is a different number of trunks and the stand basal area at the same average diameter. In practical application, such models, being a gross simplification of the reality, can be disaggregated into more detailed resolutions, providing the forest manager with more detailed information [6].

## Conclusions

The given paper analyzes the diameter distribution of small-leaved limes in forest plantations. Lime plantations were found to have certain features in the structure due to the even age, planting density and biological characteristics as a shade-tolerant specie. There is regularity in variation coefficient changes. It significantly decreases with the age of plantations. The experimental material revealed correlations between the coefficients of asymmetry and excess with the age and average diameter of small-leaved lime.

The analysis of the numerical characteristics and models of the structure of plantings by diameter gives grounds to consider the initial density of the grown plantations as the determining factor.

The analysis of tree distribution by diameter classes using the Pearson curves and the Weibull function provided the best processing of the experimental material. The simplicity of algebraic manipulation and the ability to take various forms of distribution curves make Pearson curves and the Weibull function a useful tool for forestry models.

In fact, this is the first study of this type for planted stands of small-leaved lime. The obtained data of the compiled distribution series of the total number of trunks by diameter class are recommended for describing and analyzing the structure of stands, for forest management (conducting various logging operations), development of forest inventory standards (drawing up commodity tables, guidelines for thinning), assessing sustainability, environmental factors, inventory, etc.

## Data availability

All data will be available on request.

## References

Pach M, Podlaski R. Tree diameter structural diversity in Central European forests with

*Abies alba*and*Fagus sylvatica*: managed versus unmanaged forest stands. Ecol Res. 2015;30(2):367–84.Verhunov PM. Forest taxation: a textbook. Yoshkar-Ola: Povolzhskii State Technological University Publishing; 2007.

Lima RAF, Batista JLF, Prado PI. Modeling tree diameter distributions in natural forests: an evaluation of 10 statistical models. For Sci. 2015;61(2):320–7.

Schmidt L, Machado S, Pelissari A, Silva G. Dynamics of eucalyptus diameter distribution in the State of Minas Gerais. Floresta e Ambiente. 2019;26(2):e20170156.

Meyer WH. Diameter distribution series in evenaged forest stands. Yale Sch For Bull. 1930;28:105.

Nord-Larsen T, Cao Q. A diameter distribution model for even-aged beech in Denmark. For Ecol Manag. 2006;231:218–25.

Zasada M. Evaluation of the double normal distribution for tree diameter distribution modeling. Silva Fenn. 2013;47(2):1–17.

Huang KL, Wang CH, Yen TM, Li LE. Evaluation of the applicability of the Weibull probability density function with different approaches for diameter distributions of plantations of four conifer species. Taiwan J For Sci. 2015;30(1):1–13.

Liu C, Lianjun Zh, Davis CJ, Solomon DS, Gove JH. A finite mixture model for characterizing the diameter distributions of mixed-species forest stands. For Sci. 2002;48:653–61.

Jaworski A, Kołodziej Z, Bartkowicz L. Structure and dynamics of stands of primeval character composed of the little-leaf linden (

*Tilia cordata*Mill.) in the “Las lipowy Obrożyska” reserve (southern Poland). J For Sci. 2005;51(7):283–304.Mulverhill C, Coops N, White J, Tompalski P, Marshall P, Bailey T. Enhancing the estimation of stem-size distributions for unimodal and bimodal stands in a boreal mixedwood forest with airborne laser scanning data. Forests. 2018;9(2):95.

Safonov VA, Ermakov VV, Degtyarev AP, Dogadkin NN. Prospects of biogeochemical method implementation in identifying rhenium anomalies. In: IOP conference series: earth and environmental science, vol. 421, No. 6. IOP Publishing; 2020. p. 062035.

Binoti DHB, Binoti MLM, Leite HG, Fardin L, Oliveira JC. Probability density functions for description of diameter distribution in thinned stands of

*Tectona grandis*. Cerne. 2012;18(2):185–96.Konashova SI, Sultanova RR, Khayretdinov AF, Gabdrakhimov KM, Konovalov VF, Rakhmatullin ZZ, Isyanyulova RR, Nasyrova ER, Gubydullin AF, Muftakhova SI. Forestry and ecological aspects of the broad-leaved forest formation. J Eng Appl Sci. 2018;13(11):8789–95.

Martynova M, Sultanova R, Odintsov G, Sazgutdinova R, Khanova E. Growth of

*Tilia cordata*Mill. in urban forests. SEEFOR. 2020;11(1):51–9.Souza L, Matias H, Binoti D, Leite H, Mendonça A, Fernandes da Silva G, Cruz J. Funções densidade de probabilidade para a descrição da distribuição de diâmetros em povoamentos florestais desbastados. Revista Brasileira de Biometria. 2016;34(3):421–34.

Bassil S, Nyland RD, Kern CC, Kenefic LS. Dynamics of the diameter distribution after selection cutting in uneven- and even-aged northern hardwood stands: a long-term evaluation. Can J For Res. 2019;49(12):1525–39.

Sghaier T, Cañellas I, Calama R, Sánchez-González M. Modelling diameter distribution of

*Tetraclinis articulata*in Tunisia using normal and Weibull distributions with parameters depending on stand variables. iForest. 2016;9(5):702–9.Paradis G. Lebel L. Diameter distribution models for Quebec, Canada. CIRRELT; 2017.

Sokolov PA. Identifying the structure of linden stands. For Sci. 1975;2:80–2.

Vendruscolo D, Cerqueira C, Madi JP, Silva R, Carvalho S, Medeiros R. Alterações temporais na distribuição diamétrica de teca em diferentes espaçamentos. Nativa. 2020;8(2):210–5.

Maltamo M, Mehtätalo L, Valbuena R, Vauhkonen J, Packalen P. Airborne laser scanning for tree diameter distribution modelling: a comparison of different modelling alternatives in a tropical single-species plantation. Forestry. 2018;91(1):121–31.

Kurnaev SF. Forest zoning in the USSR. Moscow: Nauka Publishing; 1973.

Bashkir encyclopedia. http://rus.bashenc.ru/. Accessed 23 Aug 2020.

Safonov VA, Danilova VN, Ermakov VV, Vorobyov VI. Mercury and methylmercury in surface waters of arid and humid regions, and the role of humic acids in mercury migration. Periodico Tche Quimica. 2019;16(31):892–902.

Iastremskii BS. Selected works. Moscow: Bukinist Publishing; 1964.

Bachioua L. On Pearson families of distributions and its applications. Afr J Math Comput Sci Res. 2013;6(5):108–17.

Shakil M, Kibria BMG, Singh J. A new family of distributions based on the generalized Pearson differential equation with some applications. Austrian J Stat. 2010;39(3):259–78.

Diamantopoulou M, Ozcelik R, Crecente-Campo F, Eler Ü. Estimation of Weibull function parameters for modelling tree diameter distribution using least squares and artificial neural networks methods. Biosyst Eng. 2015;133:33–45.

Pogoda P, Ochał W, Orzeł S. Modeling diameter distribution of black alder (

*Alnus glutinosa*(L.) Gaertn.) stands in Poland. Forests. 2019;10(5):c412.

## Acknowledgements

Not applicable.

## Funding

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

## Author information

### Authors and Affiliations

### Contributions

AG, ZR, MM, IF, and IM contributed equally to the experimentation. MM, IM, and ZR wrote and edited the article. AG, and MM equally designed and conducted the experiment. ZR, IF, studied scientific literature about the topic. All authors read and approved the final manuscript.

### Corresponding author

## Ethics declarations

### Ethics approval and consent to participate

All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.

### Consent for publication

Not applicable.

### Competing interests

The authors declare that they have no competing interests.

## Additional information

### Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Appendix

### Appendix

See Table 5.

## Rights and permissions

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

## About this article

### Cite this article

Gabdelkhakov, A., Rakhmatullin, Z., Martynova, M. *et al.* Evaluating diameter distribution series of small-leaved lime (*Tilia cordata* Mill.) in forest stands.
*Plant Methods* **17**, 45 (2021). https://doi.org/10.1186/s13007-021-00741-6

Received:

Accepted:

Published:

DOI: https://doi.org/10.1186/s13007-021-00741-6