Next Article in Journal
Combining Global Remote Sensing Products with Hydrological Modeling to Measure the Impact of Tropical Forest Loss on Water-Based Ecosystem Services
Next Article in Special Issue
Modeling Dynamics of Structural Components of Forest Stands Based on Trivariate Stochastic Differential Equation
Previous Article in Journal
Transcriptome Sequencing of Different Avocado Ecotypes: de novo Transcriptome Assembly, Annotation, Identification and Validation of EST-SSR Markers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling Diameter Distribution of Black Alder (Alnus glutinosa (L.) Gaertn.) Stands in Poland

Department of Biometry and Forest Productivity, Faculty of Forestry, University of Agriculture in Krakow, Al. 29 Listopada 46, 31-425 Krakow, Poland
*
Author to whom correspondence should be addressed.
Forests 2019, 10(5), 412; https://doi.org/10.3390/f10050412
Submission received: 9 April 2019 / Revised: 7 May 2019 / Accepted: 10 May 2019 / Published: 13 May 2019
(This article belongs to the Special Issue Modeling of Forest Structure and Dynamics)

Abstract

:
We present diameter distribution models for black alder (Alnus glutinosa (L.) Gaertn.) derived from diameter measurements made at breast height in 844 circular sample plots set in 163 managed stands located in south-eastern Poland. A total of 22,530 trees were measured. Stand age ranged from six to 89 years. The model formulation was based on the two-parameter Weibull function and a non-parametric percentile-based method. Weibull function parameters were recovered from the first raw and second central moments estimated using the stand quadratic mean diameter. The same stand characteristic was used to predict values of 12 percentiles in the percentile-based method. The model performance was assessed using the k-fold cross-validation method. The goodness-of-fit statistics include the Kolmogorov–Smirnov statistic, mean error, root mean squared error, and two variants of the error index introduced by Reynolds. The percentile model developed, accurately predicted diameter distributions in 88.4% of black alder stands, as compared to 81.9% for the Weibull model (Kolmogorov–Smirnov test). Alternative statistical metrics assessing goodness-of-fit to empirical distributions suggested that the non-parametric percentile model was superior to the parametric Weibull model, especially in stands older than 20 years. In younger stands, the two models were accurate only in 57% of the cases, and did not differ significantly with respect to goodness-of-fit measures.

1. Introduction

The natural range of the black alder (Alnus glutinosa (L.) Gaertn.) extends through most of Europe, from central Scandinavia to the southern coast of the Mediterranean Sea [1,2]. In Poland, the black alder is widespread in the lowlands, especially in the northern and central parts of the country; in the south it is less common and generally does not occur above 500 m above sea level [3]. The black alder is an economically important species in Poland, covering 522,000 ha (5.7%) with a growing stock estimated at 137 million m3 (263 m3·ha−1), which accounts for 5.3% of the country’s overall growing stock [4]. In wet soils with high groundwater levels, it mostly forms pure, even-aged, one-storey stands. On mesic sites, in mixed deciduous forests, the black alder may be used as a valuable admixture for improving timber production or for biocenotic considerations. Due to its robust root system and symbiosis with bacteria of genus Frankia, which produce nitrogen-fixing nodules on fine roots, it actively transforms its habitats [3]. In addition, as a result of open mineral metabolism, black alder leaf fall contributes 30–130 kg·ha−1·year−1 of nitrogen to the soil [1], and increases the availability of phosphorus, magnesium, manganese, and zinc [1]. These ecological properties of the black alder make it suitable for reclamation of post-industrial sites [5,6]. In addition to providing valuable timber, black alder stands for unique niches, offering habitats to riverine and riparian flora and fauna.
The implementation of economic and non-economic forest functions requires efficient forest management planning that provides reliable and detailed information about forest resources. The optimum level of detail required for economic decision-making includes data about the distribution of stand characteristics by tree-size classes [7]. Diameter distribution defined as the number of trees within a certain diameter range at breast height (DBH) [8,9], is crucial for the determination of stand volume and value [8], biodiversity [10], stability [11], development stage [12], and for predicting growth dynamics [13].
Accurate determination of diameter distribution entails substantial costs associated with measuring DBH for a large number of trees while conducting an inventory [14,15]. In practice, in order to reduce survey costs, only selected stand features are measured; these include quadratic mean diameter (QMD), mean tree height, basal area (BA), etc. Diameter distribution models can be used to estimate the number of trees falling into different DBH classes. The application of diameter distribution models in conjunction with stand characteristics generates results convergent with actual DBH measurements [16].
Diameter distribution modeling may also improve the utility of existing growth models. While complex individual tree-level models, which account for interactions between trees and their environment, accurately describe tree growth processes, they are impracticable due to the prohibitive costs of acquiring the necessary input data [17]. Diameter distribution models can be used with individual tree-level or taper models and can minimize the measurement work needed for their implementation [16,17,18]. The combination of diameter distribution models and individual tree-level or stand models gives rise to an integrated system that not only enhances the accuracy of stand characteristics assessment, but also ensures a convergence of estimates made at different levels of prediction [19].
Diameter distribution models are typically constructed using various theoretical distributions, including normal [20,21], log-normal [20,21], double normal [21,22], Weibull [8,23,24,25,26,27,28,29,30,31], Johnson SB [32,33,34,35], and beta [16,18]. From among the above, the Weibull function is most often used for modeling unimodal diameter distributions, as it has a simple form of the cumulative function [7,8,36] and due to its flexibility in fitting distributions of various shapes and degrees of skewness and ease of parameter estimation. The parameters of the probability density function may be estimated using two methods—the parameter prediction method (PPM), in which they are predicted directly from stand characteristics, and the parameter recovery method (PRM), in which they are computed from the moments or percentiles of diameter distribution predicted from stand characteristics [37].
An alternative approach to the widely used parametric distributions is the percentile-based method, in which diameter distribution is defined by 12 percentiles estimated directly from stand characteristics [38]. The method has been applied in practice by Siipilehto [39], Kangas and Maltamo [40,41], Mehtätalo [14], Stankova and Diéguez-Aranda [42], and Stankova and Zlatanov [43], and has been deemed appropriate for describing multimodal distributions and those with irregular probability density functions [38,43].
The objective of the study was to develop a diameter distribution model for black alder stands using a two-parameter Weibull function and a non-parametric percentile-based approach. The predictive power of the resulting models was evaluated using a random sample of stands from the western part of the Sandomierz Basin in southeastern Poland.

2. Materials and Methods

2.1. Data Collection

The study material consisted of DBH measurements from 844 circular sample plots in 163 managed black alder stands, aged between 6 and 89 years, and located in southeastern Poland. Regularly spaced circular sample plots, numbering about 2 to 10, with area ranging from 100 to 1000 m2 (depending on stand age) were established in each stand. In each plot, height (with an accuracy of 0.1 m) and DBH (twice, at right angles, with an accuracy of 0.1 cm, with the arithmetic mean of the two measurements adopted as the actual DBH) were measured for all living alders. A total of 22,530 trees were measured. An empirical diameter distribution for a given stand was obtained by combining the measurements from all sampled plots in that stand [21].
The collected data were used to calculate the basic stand variables widely used for the construction of diameter distribution models, that is, quadratic mean diameter (QMD; cm), stand density (SD; trees·ha−1), Reineke stand density index (SDI; trees·ha−1), Lorey’s mean height (HL; m), and basal area (BA; m2·ha−1). Furthermore, site index (SI; m) was estimated by means of a formula developed for alder stands in Poland, assuming a base age of 100 years [44].
The parameters of diameter distribution models were developed for the entire set of 163 empirical diameter distributions. The predictive accuracy and precision of the models were evaluated by the k-fold cross-validation method [45,46,47]. To that end, the set of observations was divided at random into 10 parts; nine of them constituted a training dataset. The remaining one was a testing dataset. The procedure was repeated 10 times, so that each group of stands was used as a testing dataset once. The basic statistics for the studied stands and trees are given in Table 1.

2.2. Diameter Distribution Models

The diameter distribution of the studied black alder stands was described using a two-parameter Weibull distribution and a non-parametric percentile-based method.
The probability density function f(x) of a two-parameter Weibull distribution is given as [8]:
f ( x ) = ( c b ) ( x b ) c 1 exp { ( x b ) c }   ; x 0 ;   b ,   c > 0 ,
with its cumulative distribution function F(x) being:
F ( x ) = 1 exp { ( x b ) c } ;   x 0 ;   b ,   c > 0 ,
were x is DBH (cm) and b and c are the scale and shape parameters, respectively.
Among the many ways of estimating the parameters of a Weibull distribution, the moment method is deemed superior due to its simplicity and accuracy [48,49,50]. Moreover, if parameters are calculated based on the first raw moment (m1) and the second central moment (m2), the results are unbiased estimations of these moments from the sample which guarantees agreement of statistics between the empirical and fitted distributions. The parameters b and c of the Weibull distribution were estimated based on the first raw moment and the second central moment according to the formulas:
m 1 = b Γ ( 1 + 1 c ) ,
m 2 = b 2 ( Γ ( 1 + 2 c ) Γ 2 ( 1 + 1 c ) ) ,
where Γ ( ) is the gamma function.
Transforming Equation (3) in terms of parameter b and substituting the result into Equation (4) gives:
b = m 1 Γ ( 1 + 1 c ) ,
m 2 = m 1 2 Γ 2 ( 1 + 1 c ) ( Γ ( 1 + 2 c ) Γ 2 ( 1 + 1 c ) ) ,
When the first raw moment and the second central moment are known, Equation (6) depends only on parameter c and may be solved iteratively after being transformed to:
m 1 2 Γ 2 ( 1 + 1 c ) ( Γ ( 1 + 2 c ) Γ 2 ( 1 + 1 c ) ) m 2 = 0 ,
Parameter b may be then calculated directly from Equation (5).
In the PRM, moments, m1 and m2 are estimated from the stand characteristics recorded in the inventory. The first raw moment m1, being the (arithmetic) mean diameter, may be expressed as a linear function of the QMD [42,43]:
m 1 = a + b Q M D ,
where a and b are the equation coefficients.
The second central moment m2, which is the variance of diameters, is calculated according to the formula [42,43]:
m 2 = Q M D 2 m 1 2 ,
The other approach to constructing a diameter distribution model is the percentile method proposed by Borders et al. [38], which defines the PDF by means of 12 percentiles, that is, 0th, 5th, 15th... 85th, 95th, and 100th. These percentiles are determined using empirical regression equations based on the actual stand characteristics. In the original version of the method, the percentile predicted with the greatest confidence, based on preliminary regression equations, is selected as a “driver percentile”. This work uses a version of the percentile method modified by Stankova and Zlatanov [43], without a driver percentile, in which the 12 percentiles are determined solely from stand characteristics. Preliminary analysis showed, that QMD was the best performing stand feature to predict values of the required percentiles. The parameters of the system of equations describing the empirical relationship between the 12 percentiles and the QMD of the stand were determined using seemingly unrelated regressions:
{ P 0 = a 0 + b 0 Q M D P 5 = a 5 + b 5 Q M D P 15 = a 15 + b 15 Q M D . . . P 100 = a 100 + b 100 Q M D ,
where P0, P5, P15, ..., P100 are percentiles of the empirical diameter distribution (cm), and a0, a5, a15, ..., a100 and b0, b5, b15, ..., b100 are equation coefficients.

2.3. Model Evaluation and Validation

At the stage of model parameterisation, goodness-of-fit to empirical data was evaluated using the adjusted coefficient of determination (R2adj), the mean squared error (MSE), and the significance levels of regression models (F-test), and their coefficients (t-test).
The predictive performance of the models developed was estimated using the k-fold cross-validation method. The fit of the models to data from the actual stands (testing dataset) was evaluated by comparing empirical F and predicted F ^ CDFs. The Kolmogorov–Smirnov statistic Dn was calculated to measure the maximum absolute distance between the observed and predicted CDFs [7,36,43]:
D n = sup d | F ( d i ) F ^ ( d i ) | ,
The average distances between the compared CDFs, expressed by means of the mean error (ME) and the root mean squared error (RMSE), were calculated for each stand from the following formulas:
M E = i = 1 n ( F ^ ( d i ) F ( d i ) ) n ,
R M S E = i = 1 n ( F ^ ( d i ) F ( d i ) ) 2 n ,
where di is the DBH of i-th tree in the stand and n is the number of trees.
In the case of diameter distribution modeling, it is important to not only determine differences between empirical and predicted CDFs, but also to find in which part of the distribution those differences occur. Practical implications of an error of one tree will substantially differ depending on whether that error arose in a small or large DBH class. A statistic which accounts for this additional criterion is the error index e, proposed by Reynolds et al. [36], which is the weighted sum of absolute differences between predicted and observed trees per hectare in all DBH classes [51]:
e = N k = 1 m | d L k d U k w ( d ) d F ^ ( d ) d L k d U k w ( d ) d F ( d ) | ,
where N is the number of trees per unit of area (trees·ha−1), w(d) is the weight assigned to a tree with a DBH of d, d L k and d U k are the lower and upper limits of k-th DBH class respectively, and m is the number of DBH classes. If the weight w(d) corresponds to tree volume, then the expression N d L k d U k w ( d ) d F ^ ( d ) is the predicted volume per hectare in the k-th DBH class, while N d L k d U k w ( d ) d F ( d ) is the observed volume per hectare in that DBH class. If the effect of weights is eliminated, that is w(d) = 1, the above expressions correspond to the predicted and observed number of trees per hectare for the k-th DBH class [36]. This study used two variants of the error index, eN and eG, as calculated from the formulas below:
e N = k = 1 m | N ^ k N k | ,
e G = k = 1 m | G ^ k G k | ,
where N ^ k and N k correspond to the predicted and observed number of trees per hectare for the k-th DBH class, G ^ k and G k are the predicted and observed BAs (m2·ha−1) for the k-th DBH class.
In the case of the Weibull distribution, the predicted number of trees in a given DBH class N ^ k was calculated as the product of the overall number of trees per hectare (N) determined on the basis of field measurement and the probability pk of the presence of trees in the k-th DBH class calculated from the difference in the CDF between the upper d U k and lower d L k DBH limits:
N ^ k = N p k ,
p k = F ^ ( d U k ) F ^ ( d L k ) ,
In the case of the percentile method, the predicted number of trees in a given DBH class N ^ k was determined from the following formula [38]:
N ^ k = { P i d L k P i P i 1 ( t i t i 1 ) + ( t j t i ) + ( t j + 1 t j ) d U k P j P j + 1 P j } N ,
in which:
P i 1 d L k < P i ,
P j < d U k P j + 1 ,
where: Pi, Pj are values of the r-th percentile, ti, tj are the order of the r-th percentile expressed as (ti = r/100), where {i – 1, i, j, j + 1} r; r{0, 5, 15, 25, 35, 45, 55, 65, 75, 85, 95, 100}.
The predicted BA G ^ k in the k-th DBH class was calculated from the formula:
G ^ k = N ^ k π d M k 2 40000 ,
where d M k is the middle DBH of the k-th DBH class.
Models were finally compared based on the differences revealed by the goodness-of-fit statistics for the 163 stands comprising the testing dataset in the cross-validation step. As the distributions of differences for all the analyzed statistics were not normal, they were compared using the Wilcoxon signed rank–test. At a significance level of α = 0.05, the null hypothesis (H0: Median(differences) = 0) was tested against the alternative hypothesis (HA: Median(differences) ≠ 0).
All calculations and analyses were performed in the R statistical software [52], using the following packages: base, caret [53], stats, and systemfit [54].

3. Results

3.1. Parametrization of the Diameter Distribution Models

The value of the first raw moment m1, required to determine the scale parameter b and the shape parameter c of the two-parameter Weibull distribution (Equation (1)), was calculated from the empirical formula:
m 1 = 0.1695 + 0.9847 Q M D ,
The relationship between the first raw moment m1 and the QMD was statistically significant (F(1,161) = 554072.2; p < 0.0001), with the values of the constant term (t = −5.36; p < 0.0001) and slope (t = 744.36; p < 0.0001) being significantly different from 0 (Figure 1a). The adjusted determination coefficient was 0.9997. After substituting the above equation (Equation (23)) into the empirical formula (Equation (9)) for the second central moment m2, the latter assumes the following form:
m 2 = Q M D 2 ( 0.1695 + 0.9847 Q M D ) 2 ,
with its graphical interpretation shown in Figure 1b.
Relationships between the 12 percentiles and the QMD in the stands studied are given in Figure 2. The fit of the developed system of seemingly unrelated equations assessed by McElroy’s R2 is 0.9994, with R2adj ranging from 0.850 to 0.998 for individual equations (Table 2). The lowest goodness-of-fit, as expressed by the adjusted determination coefficient and mean squared error (MSE), was found for the extreme percentiles P0 (R2adj = 0.850; MSE = 4.70) and P100 (R2adj = 0.934; MSE = 13.27). Goodness-of-fit increased for percentiles closer to the center with the best fit obtained for P55 (R2adj = 0.998; MSE = 0.15).

3.2. Evaluation and Comparison of the Distribution Models

The Kolmogorov–Smirnov test used for comparing empirical and predicted distributions, revealed that the percentile and the Weibull models accurately reflected diameter distributions in 88.4% and 81.9% of black alder stands, respectively, at a significance level of α = 0.05. The share of cases in which the null hypothesis was not rejected depended on stand age, and in consecutive 20-year age classes amounted to 57.1%, 95.1%, 95.3%, and 90.2% for the percentile model and 57.1%, 77.1%, 90.2%, and 90.2% for the Weibull model. Both models poorly predicted diameter distributions for the youngest stands (1–20 years), with the proportion of poor fits significantly declining with stand age. In the second and third age classes, the percentile model was clearly superior to the Weibull model. Figure 3 shows the examples of best and worst fitted distributions in the various age classes.
The mean Kolmogorov–Smirnov statistic Dn was 0.0786 and 0.0919 for the overall set of studied stands, for the percentile and Weibull models, respectively (Table 3). The Wilcoxon signed-rank test revealed that the two models differed significantly in terms of that criterion (Z = 7.097; p < 0.05; Table 4). In both models, the highest Dn statistic was found for stands in the first age class (Table 3), without significant differences between the models (Table 4). In subsequent age classes, the mean Dn statistic decreased, especially for the percentile model (Table 3, Figure 4), with the differences between the models being significant (Wilcoxon test, Table 4).
The mean difference (ME) between predicted and empirical CDFs was −0.0067 for the percentile model and −0.0207 for the Weibull model (Table 3). The observed absolute average differences were significantly smaller for the percentile models in all stand age classes (Table 3 and Table 4; Figure 4). RMSE differences between the predicted and empirical CDFs were 0.0372 and 0.0471 for the percentile and Weibull models, respectively (Table 3; Figure 4), and were statistically significant (Z = 8.936; p < 0.05; Table 4). For both models, RMSE was the highest in the first age class and generally decreased with stand age (Table 3; Figure 4). In the youngest stands (1–20 years), RMSE differences between the two models were not significant (Z = 1.512; p > 0.05; Table 4), in contrast to all other age classes, where statistical significance was obtained (Table 4).
The distribution of the error index eN was similar for both diameter distribution models (Figure 5). The index ranged from 83.0 to 1677.4 trees·ha−1 for the Weibull model and from 78.8 to 1997.3 trees·ha−1 for the percentile model, with the mean values being 315.0 trees·ha−1 and 314.0 trees·ha−1, respectively (Table 5). Despite these similar results, the Wilcoxon signed-rank test for all stands revealed significant differences between the median values of eN for the two models (Z = 4.499; p < 0.05; Table 6). Significant differences between the Weibull and percentile models were found in stands older than 20 years, whereas in the youngest group of stands (1–20 years), the two models revealed similar eN (Table 6; Figure 5). The same pattern was found for the error index eG, which was calculated using the BAs of the various DBH classes as weights. Also in this case, the mean error index for the percentile model was lower than for the Weibull model, both, for all stands as a whole and for individual age classes (Table 5; Figure 5). Statistically significant differences were again found for stands older than 20 years, while in the youngest ones the null hypothesis could not be rejected (Wilcoxon signed-rank test; Table 6).

4. Discussion

The study led to the development of two diameter distribution models, which in conjunction with the previously proposed generalized height-diameter model for black alder [55] may serve as a useful tool supporting the management of alder stands. In the earlier models, the frequency of trees in the various DBH classes is predicted only from the QMD, which is a commonly used measure in forestry [56]. Borders and Patterson [26] emphasize that parsimony is a desirable feature of empirical models. A similar view is expressed by the authors, the more so as the constructed models of diameter distribution will be used mainly for practical applications rather than in scientific research. That is why, a simpler two-parameter form of the Weibull distribution was applied, especially, as it often gives better results as compared to the three-parameter form [18], and excludes the problem associated with the risk of obtaining a negative scale parameter as stated by Lei [57], and also by the present authors in a previous paper [21]. To prevent a negative value for this parameter, it is often arbitrarily fixed as half of the smallest DBH in a given stand [19,28,30,58,59]. On the other hand, the act of adopting an arbitrary value for a parameter describing the PDF remains controversial and theoretically unjustified.
The other approach to modeling diameter distributions was the non-parametric percentile-based method [38], in which it is not necessary to adopt a priori any theoretical function form, and so it can be used for modeling distributions of any shape. Maltamo et al. [18] pointed out, that there are no mathematical, and especially, no biological grounds for assuming any theoretical distribution for tree diameters. Furthermore, according to Zhang et al. [51] and Palahí et al. [16], the final evaluation of a given density function often depends on the applied method of estimation of its coefficient. In the case of the method proposed by Borders et al. [38], some objections may be raised as to the selection of a given number of percentiles and the assumption as to the existence of a relationship between percentile values and stand characteristic(s). Thus, for the sake of parsimony, the present study adopted 12 percentiles describing the distribution function, but without the application of a driver percentile, especially as Stankova and Zlatanov [43] reported that the simplified method was more accurate.
The main challenge in constructing diameter distribution models for a stand, group of stands, or another set of trees is to obtain the best possible approximation of their observed distributions using the basic stand characteristics available from inventory data. The developed percentile model accurately predicted diameter distributions in 88.4% of black alder stands as compared to 81.9% for the Weibull model (Kolmogorov–Smirnov test). Similar models developed for black pine plantations in Bulgaria led to 97% and 98% agreement for the percentile and Weibull models, respectively, but it should be noted that those results were obtained with the teaching dataset, for which the models were constructed [43]. A comparison of the proposed models using other statistics evaluating goodness-of-fit to empirical distributions indicated that the non-parametric percentile model was superior to the parametric Weibull model, especially in stands older than 20 years. In the youngest stands, the two models were accurate only in 57% of the cases and did not differ significantly with respect to goodness-of-fit. These findings differ from the models developed by Stankova and Zlatanov [43] for Austrian black pine plantations, where the two-parameter Weibull distribution was markedly better than the percentile model, especially in the youngest stands. The two-parameter Weibull distribution was similarly effective in characterizing DBH for young stands in a study by Bullock and Burkhart [60].
Both models were assessed by means of different statistics including ME, RMSE and error indexes eN and eG. While the first two rely only on differences in CDFs, the third requires information about the number of trees per hectare. In our analysis, we obtain this characteristics based on field measurements on sample plots of known area. There are also other ways like the application of existing whole stand models [34] or using airborne laser scanning data [15]. It seems that, regardless of the method used to determine the number of trees per hectare, the final conclusions regarding the comparison of our two models should not be affected. However, it should be expected that the range and variability of eN and eG will be larger compared to our results, especially when the stand density is predicted based on whole stand models.
According to Nord-Larsen and Cao [17], the use of the basic stand characteristics to estimate the frequency of trees in various DBH classes may be insufficient. This objection, in particular, pertains to the fact that stands treated according to different thinning regimes may reveal similar characteristics, while being completely different in terms of diameter distribution. Our study results, especially the values of goodness-of-fit statistics and their variation, confirm that while the distribution models generated accurate estimates for some stands, in other cases the predicted results were significantly different from the observed values (Figure 3). Such a situation may be caused not only by thinning, but also by natural factors affecting diameter distribution, such as heterogeneous site conditions and irregular tree density caused by wind or snow damages. Von Gadow et al. [61] noted that a given diameter distribution may correspond to many different spatial tree distributions, and so one should incorporate information about tree location in the process of diameter distribution modeling [62]. While a greater number of factors included in a model may certainly improve prediction quality, it may also decrease the application potential of the resulting solution, thus running contrary to the parsimony principle.
While some of the latest works devoted to diameter distributions of stands have applied remote sensing tools, such as airborne laser scanning [15], at the current stage of technological development, it is still necessary to estimate those distributions indirectly by means of prediction models [15]. It seems that the constructed models may provide valuable assistance to cutting-edge technologies in modeling diameter distributions for black alder stands. Nevertheless, the primary application of the developed models is the projection of future diameter distributions to optimize management processes in black alder stands [63]. This task may be facilitated by inclusion of predicted changes in QMD over time into elaborated distribution models.

5. Conclusions

Two diameter distribution models were developed for even-aged black alder stands in Poland to estimate the frequency of trees in the various DBH classes exclusively on the basis of the QMD. The percentile-based method was found to predict the diameter distribution very accurately and was superior to the Weibull model, but its application should be limited to stands older than 20 years. In younger stands, diameter distribution was not predicted well either by the percentile or Weibull models.

Author Contributions

Data analysis: P.P. and W.O.; writing original draft and revising: P.P.; supervision, review and modification: W.O. and S.O.

Funding

The study have been financed by the Ministry of Science and Higher Education (MNiSW, Poland) as statutory funds DS-3418/ZBiPL/2018 awarded to the Department of Biometry and Forest Productivity, University of Agriculture in Krakow.

Acknowledgments

We thank the University of Agriculture in Krakow for providing empirical data collected in the 2P06L 0-5626 (KBN 2004–2005) and N N309 160,638 (NCN 2010–2012) projects that were financed by The Ministry of Science and Higher Education, Poland.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Claessens, H.; Oosterbaan, A.; Rondeux, J.; Savill, P. A review of the characteristics of black alder (Alnus glutinosa (L.) Gaertn.) and their implications for silvicultural practices. Forestry 2010, 83, 163–175. [Google Scholar] [CrossRef]
  2. Kajba, D.; Gračan, J. EUFORGEN Technical Guidelines for Genetic Conservation and Use for Black Alder (Alnus glutinosa); International Plant Genetic Resources Institute: Rome, Italy, 2003; p. 4. [Google Scholar]
  3. Pancer-Kotejowa, E.; Zarzycki, K. Zarys ekologii. In Nasze Drzewa Leśne Olsze; PWN: Warszawa-Poznań, Polska, 1980; pp. 229–249. [Google Scholar]
  4. Główny Urząd Statystyczny. Leśnictwo. In Rocznik Statystyczny Rzeczypospolitej Polskiej; Rozkrut, D., Ed.; Zakład Wydawnictw Statystycznych: Warszawa, Polska, 2018; pp. 493–497. [Google Scholar]
  5. Pietrzykowski, M.; Woś, B.; Chodak, M.; Sroka, K.; Pająk, M.; Wanic, T.; Krzaklewski, W. Effects of alders (Alnus sp.) used for reclamation of lignite combustion wastes. J. Am. Soc. Min. Reclam. 2018, 7, 30–55. [Google Scholar] [CrossRef]
  6. Sroka, K.; Chodak, M.; Klimek, B.; Pietrzykowski, M. Effect of black alder (Alnus glutinosa) admixture to Scots pine (Pinus sylvestris) plantations on chemical and microbial properties of sandy mine soils. Appl. Soil Ecol. 2018, 124, 62–68. [Google Scholar] [CrossRef]
  7. Liu, C.; Zhang, S.Y.; Lei, Y.; Newton, P.F.; Zhang, L. Evaluation of three methods for predicting diameter distributions of black spruce (Picea mariana) plantations in central Canada. Can. J. For. Res. 2004, 34, 2424–2432. [Google Scholar] [CrossRef]
  8. Bailey, R.L.; Dell, T.R. Quantifying diameter distribution with the Weibull function. For. Sci. 1973, 19, 97–104. [Google Scholar]
  9. Kangas, A.; Maltamo, M. Calibrating predicted diameter distribution with additional information. For. Sci. 2000, 46, 390–396. [Google Scholar]
  10. Maltamo, M.; Kangas, A.; Uuttera, J.; Torniainen, T.; Saramäki, J. Comparison of percentile based prediction methods and the Weibull distribution in describing the diameter distribution of heterogeneous Scots pine stands. For. Ecol. Manag. 2000, 133, 263–274. [Google Scholar] [CrossRef]
  11. Calzado Carretero, A.; Torres Alvarez, E. Modelling diameter distributions of Quercus suber L. stands in “Los Alcornocales” Natural Park (Cádiz-Málaga, Spain) by using the two parameter Weibull functions. For. Syst. 2013, 22, 15–24. [Google Scholar] [CrossRef]
  12. de Lima, R.A.F.; Batista, J.L.F.; Prado, P.I. Modeling tree diameter distributions in natural forests: An evaluation of 10 statistical models. For. Sci. 2015, 61, 320–327. [Google Scholar] [CrossRef]
  13. Gove, J.H.; Patil, G.P. Modeling the basal area-size distribution of forest stands: A compatible approach. For. Sci. 1998, 44, 285–297. [Google Scholar]
  14. Mehtätalo, L. An algorithm for ensuring compatibility between estimated percentiles of diameter distribution and measured stand variables. For. Sci. 2004, 50, 20–32. [Google Scholar]
  15. 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, 121–131. [Google Scholar] [CrossRef]
  16. Palahí, M.; Pukkala, T.; Blasco, E.; Trasobares, A. 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). Eur J. For. Res. 2007, 126, 563–571. [Google Scholar] [CrossRef]
  17. Nord-Larsen, T.; Cao, Q.V. A diameter distribution model for even-aged beech in Denmark. For. Ecol. Manag. 2006, 231, 218–225. [Google Scholar] [CrossRef]
  18. Maltamo, M.; Puumalainen, J.; Päivinen, R. Comparison of beta and Weibull functions for modelling basal area diameter distribution in stands of Pinus sylvestris and Picea abies. Scand. J. For. Res. 1995, 10, 284–295. [Google Scholar] [CrossRef]
  19. Zhang, X.; Lei, Y. A linkage among whole-stand model, individual-tree model and diameter-distribution model. J. For. Sci. 2010, 56, 600–608. [Google Scholar] [CrossRef] [Green Version]
  20. Ochał, W.; Pająk, M.; Pietrzykowski, M. Struktura grubości wybranych drzewostanów sosnowych wzrastających na zrekultywowanych dla leśnictwa obiektach pogórniczych. [Diameter structure of selected pine stands growing on post-mining sites reclaimed for forestry]. Sylwan 2010, 154, 323–332. [Google Scholar]
  21. Orzeł, S.; Pogoda, P.; Ochał, W. Ocena przydatności wybranych funkcji do modelowania rozkładu pierśnic w drzewostanach olszy czarnej (Alnus glutinosa (L.) Gaertn.). [Evaluation of usefulness of selected functions for modeling distribution of breast height diameter in black alder stands (Alnus glutinosa (L.) Gaertn.)]. Sylwan 2017, 161, 101–113. [Google Scholar]
  22. Zasada, M. Evaluation of the double normal distribution for tree diameter distribution modeling. Silva. Fenn. 2013, 47, 1–17. [Google Scholar] [CrossRef]
  23. Ek, A.R.; Issos, J.N.; Bailey, R.L. Solving for Weibull diameter distribution parameters to obtain specified mean diameters. For. Sci. 1975, 21, 290–292. [Google Scholar]
  24. Shiver, B.D. Sample sizes and estimation methods for the Weibull distribution for unthinned slash pine plantation diameter distributions. For. Sci. 1988, 34, 809–814. [Google Scholar]
  25. Kilkki, P.; Maltamo, M.; Mykkänen, R.; Päivinen, R. Use of the Weibull function in estimating the basal area dbh-distribution. Silva. Fenn. 1989, 23, 311–318. [Google Scholar] [CrossRef]
  26. Borders, B.E.; Patterson, W.D. Projecting stand tables: A comparison of the Weibull diameter distribution method, a percentile-based projection method, and a basal area growth projection method. For. Sci. 1990, 36, 413–424. [Google Scholar]
  27. Knowe, S.A.; Ahrens, G.R.; DeBell, D.S. Comparison of diameter-distribution-prediction, stand-table-projection and individual-tree-growth modelling approaches for young red alder plantations. For. Ecol. Manag. 1997, 98, 49–60. [Google Scholar] [CrossRef]
  28. Cao, Q.V. Predicting parameters of a Weibull function for modeling diameter distribution. For. Sci. 2004, 50, 682–685. [Google Scholar]
  29. Newton, P.F.; Lei, Y.; Zhang, S.Y. A parameter recovery model for estimating black spruce diameter distributions within the context of a stand density management diagram. For. Chron. 2004, 80, 349–358. [Google Scholar] [CrossRef] [Green Version]
  30. Ercanli, I.; Bolat, F.; Kahriman, A. Comparing Parameter Recovery Methods for Diameter Distribution Models of Oriental Spruce (Picea orientalis (L.) Link.) and Scotch Pine (Pinus sylvestris L.) mixed stands located Trabzon and Giresun Forest Regional Directorate; Oral presentation. International Caucasian Forestry Symposium: Artvin, Turkey, 24–26 October 2013. [Google Scholar]
  31. Siipilehto, J.; Mehtätalo, L. Parameter recovery vs. parameter prediction for the Weibull distribution validated for Scots pine stands in Finland. Silva. Fenn. 2013, 47. [Google Scholar] [CrossRef] [Green Version]
  32. Mønness, E. Diameter Distributions and Height Curves in Even-Aged Stands of Pinus sylvestris L; Reports of the Norwegian Forest Research Institute: Ås, Norway, 1982; Volume 36, pp. 1–40. [Google Scholar]
  33. Fonseca, T.F.; Marques, C.P.; Parresol, B.R. Describing Maritime Pine diameter distribution with Johnson’s SB distribution using a new all-parameter recovery approach. For. Sci. 2009, 55, 367–373. [Google Scholar]
  34. Mateus, A.; Tomé, M. Modelling the diameter distribution of eucalyptus plantations with Johnson’s SB probability density function: Parameters recovery from a compatible system of equations to predict stand variables. Ann. For. Sci. 2011, 68, 325–335. [Google Scholar] [CrossRef]
  35. Siipilehto, J. Methods and Applications for Improving Parameter Prediction Models for Stand Structures in Finland. Ph.D Thesis, University of Helsinki, Helsinki, Finland, September 2011; p. 56. [Google Scholar]
  36. Reynolds, M.R.; Burk, T.E.; Huang, W.C. Goodness-of-fit tests and model selection procedures for diameter distribution models. For. Sci. 1988, 34, 373–399. [Google Scholar]
  37. Hyink, D.M.; Moser, J.W., Jr. A generalized framework for projecting forest yield and stand structure using diameter distributions. For. Sci. 1983, 29, 85–95. [Google Scholar]
  38. Borders, B.E.; Souter, R.A.; Bailey, R.L.; Ware, K.D. Percentile-based distributions characterize forest stand tables. For. Sci. 1987, 33, 570–576. [Google Scholar]
  39. Siipilehto, J. Improving the accuracy of predicted basal-area diameter distribution in advanced stands by determining stem number. Silva. Fenn. 1999, 33, 281–301. [Google Scholar] [CrossRef]
  40. Kangas, A.; Maltamo, M. Percentile based basal area diameter distribution models for Scots Pine, Norway Spruce and Birch Species. Silva. Fenn. 2000, 34, 371–380. [Google Scholar] [CrossRef]
  41. Kangas, A.; Maltamo, M. Performance of percentile based diameter distribution prediction and Weibull method in independent data sets. Silva. Fenn. 2000, 34, 381–398. [Google Scholar] [CrossRef]
  42. Stankova, T.; Diéguez-Aranda, U. Diameter distribution model for Scots Pine plantations in Bulgaria. For. Ideas 2010, 16, 155–162. [Google Scholar]
  43. Stankova, T.V.; Zlatanov, T.M. Modeling diameter distribution of Austrian black pine (Pinus nigra Arn.) plantations: A comparison of the Weibull frequency distribution function and percentile-based projection methods. Eur J. For. Res. 2010, 129, 1169–1179. [Google Scholar] [CrossRef]
  44. Socha, J.; Ochał, W. Dynamic site index model and trends in changes of site productivity for Alnus glutinosa (L.) Gaertn. in southern Poland. Dendrobiology 2017, 77, 45–57. [Google Scholar] [CrossRef]
  45. Pandey, R.; Dhall, S.P.; Kumar, R. Cross-validation–a tool for forest modelers. Indian For. 1999, 125, 1224–1227. [Google Scholar]
  46. Refaeilzadeh, P.; Tang, L.; Liu, H. Cross-Validation. In Encyclopedia of Database Systems; Liu, L., Özsu, M.T., Eds.; Springer: Boston, MA, USA, 2009; p. 6. [Google Scholar]
  47. Arlot, S.; Celisse, A. A survey of cross-validation procedures for model selection. Stat. Surv. 2010, 4, 40–79. [Google Scholar] [CrossRef] [Green Version]
  48. Merganič, J.; Sterba, H. Characterisation of diameter distribution using the Weibull function: Method of moments. Eur J. For. Res. 2006, 125, 427–439. [Google Scholar] [CrossRef]
  49. Al-Fawzan, M. Methods for Estimating the Parameters of Weibull Distribution; King Abdulaziz City for Science and Technology: Riyadh, Saudi Arabia, 2000. [Google Scholar]
  50. Nanang, D. Suitability of the Normal, Log-normal and Weibull distributions for fitting diameter distributions of neem plantations in Northern Ghana. For. Ecol. Manag. 1998, 103, 1–7. [Google Scholar] [CrossRef]
  51. Zhang, L.; Packard, K.C.; Liu, C. A comparison of estimation methods for fitting Weibull and Johnson’s SB distributions to mixed spruce-fir stands in northeastern North America. Can. J. For. Res. 2003, 33, 1340–1347. [Google Scholar] [CrossRef]
  52. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  53. Kuhn, M.; Wing, J.; Weston, S.; Williams, A.; Keefer, C.; Engelhardt, A.; Cooper, T.; Mayer, Z.; Kenkel, B.; Benesty, M.; et al. Caret: Classification and Regression Training. R package version 6.0-81. Available online: https://CRAN.R-project.org/package=caret (accessed on 10 February 2019).
  54. Henningsen, A.; Hamann, J.D. Systemfit: A Package for Estimating Systems of Simultaneous Equations in R. J. Stat. Softw. 2007, 23, 1–40. [Google Scholar] [CrossRef]
  55. Orzeł, S.; Pogoda, P.; Ochał, W. Stała krzywa wysokości dla olszy czarnej (Alnus glutinosa (L.) Gaertn.) z zachodniej części Kotliny Sandomierskiej. [Generalized height-diameter model for black alder (Alnus glutinosa (L.) Gaertn.) in the western part of the Sandomierz Basin]. Sylwan 2014, 158, 840–849. [Google Scholar]
  56. Curtis, R.O.; Marshall, D.D. Why quadratic mean diameter? West. J. Appl. For. 2000, 15, 137–139. [Google Scholar]
  57. Lei, Y. Evaluation of three methods for estimating the Weibull distribution parameters of Chinese pine (Pinus tabulaeformis). J. For. Sci. 2008, 54, 566–571. [Google Scholar] [CrossRef]
  58. Gorgoso, J.J.; Rojo, A.; Camara-Obregon, A.; Dieguez-Aranda, U. A comparison of estimation methods for fitting Weibull, Johnson’s SB and beta functions to Pinus pinaster, Pinus radiata and Pinus sylvestris stands in northwest Spain. For. Syst. 2012, 21, 446–459. [Google Scholar] [CrossRef]
  59. Poudel, K.P.; Cao, Q.V. Evaluation of methods to predict Weibull parameters for characterizing diameter distributions. For. Sci. 2013, 59, 243–252. [Google Scholar] [CrossRef]
  60. Bullock, B.P.; Burkhart, H.E. Juvenile diameter distributions of loblolly pine characterized by the two-parameter Weibull function. New For. 2005, 29, 233–244. [Google Scholar] [CrossRef]
  61. Von Gadow, K.; Zhang, C.Y.; Wehenkel, C.; Pommerening, A.; Corral-Rivas, J.; Korol, M.; Myklush, S.; Hui, G.Y.; Kiviste, A.; Zhao, X.H. Forest structure and diversity. In Continuous Cover Forestry; Pukkala, T., Von Gadow, K., Eds.; Springer: Dordrecht, The Netherlands, 2012; pp. 29–83. [Google Scholar]
  62. García, O. What is a diameter distribution? In Proceedings of the Symposium on Integrated Forest Management Information Systems, Tsukuba, Japan, 13–18 October 1991; Minowa, M., Tsuyuki, S., Eds.; Japan Society of Forest Planning Press: Tsukuba, Japan, 1992; pp. 11–29. [Google Scholar]
  63. Goelz, J.C.G. Forest growth and yield models viewed from different perspective. Math. Modeling Sci. Comput. 2002, 13, 177–189. [Google Scholar]
Figure 1. Relationship between (a) the first raw moment m1, (b) the second central moment m2 and the quadratic mean diameter (QMD) of black alder stands.
Figure 1. Relationship between (a) the first raw moment m1, (b) the second central moment m2 and the quadratic mean diameter (QMD) of black alder stands.
Forests 10 00412 g001
Figure 2. Relationship between percentiles (P0, P5, P15, ..., P100) and the quadratic mean diameter (QMD) of black alder stands.
Figure 2. Relationship between percentiles (P0, P5, P15, ..., P100) and the quadratic mean diameter (QMD) of black alder stands.
Forests 10 00412 g002
Figure 3. Examples of the best (ad) and worst (eh) fit of developed percentile (solid line) and Weibull (dotted line) models for test data in subsequent age classes (from top to bottom). Dn(W) and Dn(P) are the Kolmogorov–Smirnov statistics for Weibull and percentile models, respectively.
Figure 3. Examples of the best (ad) and worst (eh) fit of developed percentile (solid line) and Weibull (dotted line) models for test data in subsequent age classes (from top to bottom). Dn(W) and Dn(P) are the Kolmogorov–Smirnov statistics for Weibull and percentile models, respectively.
Forests 10 00412 g003
Figure 4. Distributions of goodness-of-fit statistics: (a,b) Kolmogorov–Smirnov Dn statistic, (c,d) mean error (ME), (e,f) root mean squared error (RMSE) for the fitted Weibull and percentile models in black alder stands.
Figure 4. Distributions of goodness-of-fit statistics: (a,b) Kolmogorov–Smirnov Dn statistic, (c,d) mean error (ME), (e,f) root mean squared error (RMSE) for the fitted Weibull and percentile models in black alder stands.
Forests 10 00412 g004
Figure 5. Distributions of error indexes eN (a,b), and eG (c,d) for the fitted Weibull and percentile models in black alder stands.
Figure 5. Distributions of error indexes eN (a,b), and eG (c,d) for the fitted Weibull and percentile models in black alder stands.
Forests 10 00412 g005
Table 1. Descriptive statistics for stand and tree variables.
Table 1. Descriptive statistics for stand and tree variables.
StatisticsStand VariableTree Variable
QMDAgeSDSDIHLSIBADBHTHBA
Mean21.9944102052820.4932.425.2120.318.90.04120
StDev9.44238851296.372.49.1110.66.80.03951
Min.4.4762232126.2826.55.481.02.80.00008
Max.42.9589436087332.7540.249.7967.036.60.35257
QMD (cm), quadratic mean diameter; SD (stems·ha−1), stand density; SDI (stems·ha−1), Raineke stand density index; HL (m), Lorey’s mean height; SI (m), site index (base age 100 years), BA (m2·ha−1), basal area; DBH (cm), diameter at breast height; TH (m), tree height; StDev, standard deviation.
Table 2. Parameters and goodness-of-fit statistics for equations describing relationships between percentiles and the quadratic mean diameter (QMD) of black alder stands.
Table 2. Parameters and goodness-of-fit statistics for equations describing relationships between percentiles and the quadratic mean diameter (QMD) of black alder stands.
PercentilesParameters a0–100Parameters b0–100R2adjMSE
P0 = a0 + b0 QMD0.0000.534 ***0.8504.70
P5 = a5 + b5 QMD−0.663 ***0.690 ***0.9492.34
P15 = a15 + b15 QMD−1.039 ***0.810 ***0.9850.89
P25 = a25 + b25 QMD−0.929 ***0.872 ***0.9920.56
P35 = a35 + b35 QMD−0.714 ***0.918 ***0.9960.32
P45 = a45 + b45 QMD−0.451 ***0.958 ***0.9970.21
P55 = a55 + b55 QMD−0.155 *0.996 ***0.9980.15
P65 = a65 + b65 QMD0.0001.048 ***0.9980.17
P75 = a75 + b75 QMD0.284 ***1.101 ***0.9980.23
P85 = a85 + b85 QMD0.681 ***1.164 ***0.9960.53
P95 = a95 + b95 QMD1.186 ***1.280 ***0.9881.76
P100 = a100 + b100 QMD2.574 ***1.468 ***0.93413.27
Significance levels: *** p < 0.001, * p < 0.05; R2adj, adjusted coefficient of determination; MSE, mean squared error.
Table 3. Goodness-of-fit statistics for Weibull and percentile diameter distribution models for different stand age classes.
Table 3. Goodness-of-fit statistics for Weibull and percentile diameter distribution models for different stand age classes.
Distribution ModelsAge Classes of StandsDn StatisticMERMSE
MeanSDMeanSDMeanSD
Weibull1–20 years0.09770.0308−0.01470.01750.05050.0180
21–40 years0.09370.0299−0.02480.01620.04910.0181
41–60 years0.08810.0240−0.02060.00990.04440.0133
>60 years0.09000.0277−0.02090.01250.04560.0163
All stands0.09190.0279−0.02070.01420.04710.0164
Percentile1–20 years0.09780.0368−0.00180.02140.04730.0209
21–40 years0.07850.0259−0.01220.01880.03840.0151
41–60 years0.07260.0282−0.00630.01250.03310.0143
>60 years0.07170.0316−0.00510.01550.03360.0177
All stands0.07860.0315−0.00670.01710.03720.0175
ME, mean error; RMSE, root mean squared error; SD, standard deviation.
Table 4. Results of the Wilcoxon signed-rank test for comparison of goodness-of-fit statistics for the Weibull and percentile diameter distribution models.
Table 4. Results of the Wilcoxon signed-rank test for comparison of goodness-of-fit statistics for the Weibull and percentile diameter distribution models.
Age Classes of StandsDn StatisticMERMSE
Z Statisticp-ValueZ Statisticp-ValueZ Statisticp-Value
1–20 years0.2570.79714.7410.00001.5120.1310
21–40 years4.1300.00005.7110.00005.2040.0000
41–60 years4.6220.00005.8410.00005.4010.0000
>60 years4.6560.00005.8410.00005.1080.0000
All stands7.0970.000011.0700.00008.9360.0000
ME, mean error; RMSE, root mean squared error.
Table 5. Error indexes eN and eG for the Weibull and percentile diameter distribution models for different stand age classes.
Table 5. Error indexes eN and eG for the Weibull and percentile diameter distribution models for different stand age classes.
Distribution ModelsAge Classes of StandseNeG
MeanSDMeanSD
Weibull1–20 years675.7378.64.102.17
21–40 years352.3144.87.743.72
41–60 years207.470.310.243.67
>60 years146.534.613.043.42
All stands315.0261.49.224.58
Percentile1–20 years749.1483.73.921.59
21–40 years322.2132.56.863.08
41–60 years194.366.09.503.55
>60 years136.029.611.742.91
All stands314.0309.48.394.05
Table 6. Results of the Wilcoxon signed-rank test for comparison of error indexes eN and eG for Weibull and percentile diameter distribution models.
Table 6. Results of the Wilcoxon signed-rank test for comparison of error indexes eN and eG for Weibull and percentile diameter distribution models.
Age Classes of StandseNeG
Z Statisticp-ValueZ Statisticp-Value
1–20 years1.1000.27120.5660.5717
21–40 years3.3450.00084.3710.0000
41–60 years3.3470.00083.9790.0001
>60 years4.3060.00005.0170.0000
All stands4.4990.00007.5570.0000

Share and Cite

MDPI and ACS Style

Pogoda, P.; Ochał, W.; Orzeł, S. Modeling Diameter Distribution of Black Alder (Alnus glutinosa (L.) Gaertn.) Stands in Poland. Forests 2019, 10, 412. https://doi.org/10.3390/f10050412

AMA Style

Pogoda P, Ochał W, Orzeł S. Modeling Diameter Distribution of Black Alder (Alnus glutinosa (L.) Gaertn.) Stands in Poland. Forests. 2019; 10(5):412. https://doi.org/10.3390/f10050412

Chicago/Turabian Style

Pogoda, Piotr, Wojciech Ochał, and Stanisław Orzeł. 2019. "Modeling Diameter Distribution of Black Alder (Alnus glutinosa (L.) Gaertn.) Stands in Poland" Forests 10, no. 5: 412. https://doi.org/10.3390/f10050412

APA Style

Pogoda, P., Ochał, W., & Orzeł, S. (2019). Modeling Diameter Distribution of Black Alder (Alnus glutinosa (L.) Gaertn.) Stands in Poland. Forests, 10(5), 412. https://doi.org/10.3390/f10050412

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop