Next Article in Journal
Estimate Forest Aboveground Biomass of Mountain by ICESat-2/ATLAS Data Interacting Cokriging
Next Article in Special Issue
Evaluation of the Community Land Model-Simulated Specific Leaf Area with Observations over China: Impacts on Modeled Gross Primary Productivity
Previous Article in Journal
Effect of Type of Forest Growth Conditions and Climate Elements on the Dynamics of Radial Growth in English Oak (Quercus robur L.) of Early and Late Phenological Forms
Previous Article in Special Issue
Feasibility of Agarwood Cultivation in Indonesia: Dynamic System Modeling Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling Number of Trees per Hectare Dynamics for Uneven-Aged, Mixed-Species Stands Using the Copula Approach

by
Petras Rupšys
1,2,* and
Edmundas Petrauskas
2
1
Faculty of Informatics, Vytautas Magnus University, LT 44404 Kaunas, Lithuania
2
Agriculture Academy, Vytautas Magnus University, LT 53361 Kaunas, Lithuania
*
Author to whom correspondence should be addressed.
Forests 2023, 14(1), 12; https://doi.org/10.3390/f14010012
Submission received: 19 November 2022 / Revised: 7 December 2022 / Accepted: 15 December 2022 / Published: 21 December 2022
(This article belongs to the Special Issue Modelling Forest Ecosystems)

Abstract

:
For the monitoring and management of forest resources, the main index is the stand volume, which is determined on the basis of the tree diameter, height, and number of trees per hectare of three-dimensional distribution. The development of trees in the forest stand is dynamic and is driven by random phenomena. In this study, the tree diameter, the potentially available area, and the height are described by the mixed-effect parameters of the Gompertz-type diffusion process. A normal copula function is used to connect a three-dimensional distribution to its one-dimensional margins. The newly developed model was illustrated using empirical data from 53 permanent experimental plots (measured for seven cycles), which were characterized as follows: pine forests (Pinus sylvestris), 63.8%; spruce (Picea abies), 30.2%; silver birch (Betula pendula Roth and Betula pubescens Ehrh.), 5.8%; and others, 0.2%. An analysis of the tree diameter and height of growth, including current and mean increments and inflection points, is presented. The models for the change in the number of trees per hectare with age are presented on the basis of the probabilistic density functions of the solutions of stochastic differential equations and the copula function. The dynamics of the number of trees per hectare are visualized graphically, and the goodness of fit of the newly developed models is evaluated using standard statistical measures.

1. Introduction

1.1. Background

Most often, tree species in the forest develop from generation to generation. The progressive changes in the size, shape, and function of the trees in a stand during their lifetimes, which transform the stand’s genetic potential into a living mature stand, define the dynamism of the stand and are referred to as stand growth. The development of trees in a forest stand can be understood as the interrelationship of the different tree sizes, which allows a particular tree to overcome the chaotic behavior of nature and proves that the stand exists as an integrated system, where changes in one tree size ultimately cause a corresponding change in the size of another. Trees in a forest stand do not live in isolation, as they are shaped in accordance with their interactions with their environment. The growth equation mathematically combines the two dependent variables of age and size to explain the continuous process of the two hidden forces of body expansion and size limitation [1]. As the size of a tree is determined by several components (for example, diameter, height, crown width, etc.), there exists a specific relationship between these size components during the entire growth process. Therefore, it is appropriate to expand the interpretation of the growth process of a particular size component by additionally introducing other size components and clarifying the extent of the dependencies between the response and the explanatory variables [2,3].

1.2. Research Motivation

Complex forest systems are characterized by important data features, such as dynamism, uncertainty, stochasticity, deviation of variables from the normal distribution, multicollinearity, and false correlations, which complicate the modeling and prediction of these systems. An approach is required that, under the conditions of the described data features, is able to adequately identify the processes in such systems, ensuring a high level of accuracy and efficiency. To assess the level of productivity of forest stands and to formalize their dynamics, it is appropriate to determine the link between the main components of the stand size, namely the number of trees per hectare, the average diameter of the trees at breast height (or average square diameter), and the average height of the trees using univariate or multivariate diffusion process models. The growth processes of forest stands evolve under the influence of random noise, which is mostly expressed by Brownian motion. Traditionally, special types of stochastic differential equations have been used to model random phenomena in forestry on the basis of the statistical properties of the observed data or on the basis of the researcher’s experience in stochastic calculus since the stochastic differential equations governing random phenomena are generally unknown. In stochastic, stand-level growth models, the transition probability density functions were pioneered by Suzuki [4], Sloboda [5], Garcia [6], Tanaka [7], and Rennolls [8] to describe the movement of individual trees with respect to the diameter or height classes. The parametric identification of the developed stochastic differential equations [4,5,6,7,8] required unrealistic information storage and computing resources that were difficult to satisfy at that time. Naturally, this methodology was not applied in practice due to its theoretical complexity for foresters and the large volume of computer calculations required. Stochastic differential equations can simultaneously capture the known sigmoidal or exponential deterministic dynamics of underlying variables of interest, such as tree height [9], tree diameter [10], tree potentially available area [11], tree crown width, crown base height [12,13], and the others [14,15], while enabling the researcher to capture the unknown random dynamics in a stochastic setting. In the long run, as age approaches infinity, the (unconditional) stationary approximations hold for the Gompertz- [16] and Vasicek-type [17] stochastic growth models.

1.3. Research Objectives

In this paper, we propose a copula technique [18] that provides methodological support for multi-dimensional data analysis with the aim of uncovering the hidden dependencies between size components and the responses to changing behavior in time and identifying the non-obvious regularities in forestry data using machine learning techniques. Copula functions have been used in forestry applications since 2000 [19], as have the references therein. The applications of multivariate diffusion processes in forestry have shown that the multivariate diffusion process is suitable for modeling the multivariate distribution of tree size components, but the implementation of the maximum likelihood method for finding the parameter estimates of such a model requires a lot of computer calculation time and memory. Copula models are based on Sklar’s theorem [20], which states that every multivariate cumulative distribution function can be expressed as a copula function depending on the univariate marginal distributions that capture the structure of the dependencies between the marginal components.
Previously published works have presented the cases of two-dimensional [21], three-dimensional [22], and four-dimensional [23] stochastic differential equations that relate singular stochastic differential equations using the covariance matrix between the individual variables, thus revealing the dependency structure. The growth model defined by the multi-dimensional stochastic differential equations had drawbacks in that, firstly, the discrete variable measurements for each observed variable need to have the same amount of measurements, and, secondly, the large number of unknown parameters characterizing the multi-dimensional system necessitates computer calculations that require a very long time before reaching an appropriate level of convergence. In this study, we focused on a trivariate normal copula that is defined by a trivariate normal probability density function with a corresponding correlation matrix and with marginal distributions defined by asymmetric lognormal distributions. Theoretical studies of tree growth have confirmed that tree size distribution is inverted J-shaped, with many small trees and few larger trees due to asymmetric competition [24].
The main objective of this study is to improve the estimates of the growth and yield capacity in Lithuanian stands of uneven-aged, mixed-species forests and to reveal the main moments of stand age, which are of special importance in the growth process. The functional relationships between the number of trees per hectare were modeled using stochastic differential equations and copula methods. In this paper, we propose a copula technique that provides methodological support for multivariate data analysis using stochastic differential equations, making it possible to reveal the degree of hidden dependencies and to identify implicit patterns in the forest growth data using machine learning methods. The proposed technique makes it possible, firstly, to carry out a systematic analysis of the growth of uneven-aged, mixed-species forest stands, taking into account the uncertainty of the internal and external environment, based on quantitative statistical analysis tools supported by univariate and multivariate density functions. Secondly, it builds a theoretically acceptable model of the observed forest stand and creates machine learning algorithms that enable the implementation of the modeled object in practice. Thirdly, it can serve as the basis of a decision-support system for the management of uneven-aged, mixed-species forest stands. Fourthly, the random effect scenario mode accounts for the unobserved ecological factors such as species composition, temperature, solar radiation, precipitation, nutrient condition, and much more. We illustrate the potential of this new framework using the example of longitudinal measurements (up to seven times) in the Kazlų Rūda district in Lithuania on the basis of the measurements of 53 stands, recording the tree’s age, diameter, height (approximately every fifth tree), and the coordinates of the tree’s position.

2. Materials and Methods

In this section, we introduce a stochastic four-parameter Gompertz-type model that describes the dynamics of the tree diameter at breast height (henceforth: diameter), the potentially available area, and the height over time (age). In the subsequent analysis, the potentially available area of a tree will be expressed by the corresponding area of the polygon of the Voronoi diagram since the coordinates of the position of each tree in the plot were recorded during the measurements (see Figure 1, rectangular plot 50 m × 100 m, 0.5 ha, 751 trees).

2.1. Stochastic Differential Equation Framework

The mixed-effect parameters of the Gompertz-type one-dimensional stochastic differential equations for the tree diameter, X d i ( t ) , potentially available area, X p i ( t ) , and height, X h i ( t ) , take the following forms, respectively:
d X d i ( t ) = ( ( α d + ϕ d i ) β d l n ( X d i ( t ) ɣ d ) ) ( X d i ( t ) ɣ d ) d t + σ d ( X d i ( t ) ɣ d ) d W d i ( t ) , P ( X d i ( t 0 ) = d 0 ) = 1 , i = 1 , , M ,
d X p i ( t ) = ( ( α p + ϕ p i ) β p l n ( X p i ( t ) ɣ p ) ) ( X p i ( t ) ɣ p ) d t + σ p ( X p i ( t ) ɣ p ) d W p i ( t ) , P ( X p i ( t 0 ) = δ ) = 1 , i = 1 , ,   M ,
d X h i ( t ) = ( ( α h + ϕ h i ) β h l n ( X h i ( t ) ɣ h ) ) ( X h i ( t ) ɣ h ) d t + σ h ( X h i ( t ) ɣ h ) d W h i ( t ) , P ( X d i ( t 0 ) = h 0 ) = 1 ,   i = 1 , ,   M ,
where M is the number of observed plots; α d , α p , α h and βd, βp, βh are the birth and death rate fixed-effect parameters α j   > 0, β j > 0, j = d, p, h); σ d , σ p , σ h are the volatility fixed-effect parameters; ɣ d , ɣ p , ɣ h are the threshold fixed-effect parameters; δ is the fixed-effect parameter; the random effects ϕ d i , ϕ p i , ϕ h i are the independent and normally distributed random variables with zero mean and constant variances τ j 2 , respectively, ϕ j i ~ N ( 0 ; τ j 2 ) , j = d, p, h; and d W j i ( t ) , j = d, p, h, represent the Brownian motion increments, which are considered to be independent across all the stands. The unknown fixed-effect parameters, θ = { α j , β j , ɣ j , σ j , τ j , δ ; j = d ,   p , h }, must be estimated.
Transforming Equations (1)–(3) by e β j t l n ( X j i ( t ) ) i = 1, …, M, j = d, p, h, and using Itô’s lemma [25], it can be determined that ( X j i ( t ) | X j i ( t 0 ) = x j 0 ) conforms to a lognormal distribution L N 1 ( μ j i ( t ) ; v j ( t ) ) , with the mean μ j i ( t ) , variance v j ( t ) and probability density function f j i ( x j , t | θ , ϕ j i ) in the following form:
μ j i ( t ) = l n ( x j 0 ɣ j ) e β j ( t t 0 ) + α j + ϕ j i β j ( 1 e β j ( t t 0 ) ) ,
v j ( t ) = 1 e 2 β j ( t t 0 ) 2 β j σ j ,
f j i ( x j , t | θ , ϕ j i ) = 1 2 π v j ( t ) ( x ɣ j ) exp ( ( l n ( x ɣ j ) μ j i ( t ) ) 2 2 v j ( t ) ) .
The newly derived probability density functions prove themselves for the quantitative assessment of tree and stand growth using mathematical–statistical methods and form a useful basis for research in forestry. In the analysis based on stochastic growth processes, any repeated surveys provide data for the fitting of stochastic differential Equations (1)–(3), and the various attributes and characteristics of growth, such as current annual increment, etc., are calculated from the fitted probability density functions rather than directly from the observed data.
In our study, the dynamics of the state (diameter, potentially available area and height) of a tree’s statistical growth characteristics at time t, e.g., the mean, m j i ( t ) ; median, m e j i ( t ) ; mode, m o j i ( t ) ; qth quantile (0 < q < 1), q m j i ( t ) ; and variance, w j i ( t ) , i = 1, …, M, j = d, p, h, using the properties of the lognormal distribution, take the following expressions:
m j i ( t ) = ɣ j + exp ( μ j i ( t ) + 1 2 v j ( t ) ) ,
m e j i ( t ) = ɣ j + exp ( μ j i ( t ) ) ,
m o j i ( t ) = ɣ j + exp ( μ j i ( t ) v j ( t ) ) ,
q m j i ( t , q ) = ɣ j + exp ( μ j i ( t ) + v j ( t ) Φ q 1 ( 0 ; 1 ) ) ,
w j i ( t ) = exp ( 2 μ j i ( t ) ( t ) + v j ( t ) ) ( exp ( v ( t ) ) 1 ) .
Using a three-dimensional normal distribution, we define a three-dimensional normal copula cumulative distribution function C that would relate all three tree state variables (diameter, potentially available space, height) by defining the corresponding three-dimensional copula-type distribution (see Appendix A), as follows:
C ( u 1 , u 2 , u 3 ; t ) = H t ( F 1 1 ( u 1 , t ) , F 2 1 ( u 2 , t ) , F 3 1 ( u 3 , t ) ) , F 1 ( x 1 , t ) = P ( X 1 ( t ) x 1 ) ,   F 2 ( x 2 , t ) = P ( X 2 ( t ) x 2 ) ,   F 3 ( x 3 , t ) = P ( X 3 ( t ) x 3 ) .  
This three-dimensional distribution function enables us to formalize the relationships of one of the tree state variables with the others. Using the integration operation and the copula-type conditional probability density functions f j | k i ( x j , t | x k ) , i = 1, …, M, j = d, p, h; kj, and f j | k , l i ( x j , t | x k , x l ) ,  i = 1, …, M, j = d, p, h; k, lj (see Appendix A), we are able to express the relationships of the mean trend of one of the state variables with the others, including the dependences of the state variables, in the following form:
m j | k i ( x k , t ) = ɣ j + x j f j | k i ( x j , t | x k ) d x j .
m j | k , l i ( x k , x l , t ) = ɣ j + x j f j | k , l i ( x j , t | x k , x l ) d x j .
An analogy with the solution of mechanical problems makes it possible to define the current annual and mean increments and their confidence intervals, relative increment, growth inflection point, etc., as follows, i = 1, …, M, j = d, p, h; k, lj:
The current annual increment:
c m j i ( t ) = m j i ( t ) t   or   c m j | k i ( x k , t ) = m j | k i ( x k , t ) t ,   c m j | k , l i ( x k , x l , t ) = m j | k , l i ( x k , x l , t ) t ,
The mean annual increment:
m m j i ( t ) = m j i ( t ) t   or   m m j | k i ( x k , t ) = m j | k i ( x k , t ) t ,   m m j | k , l i ( x k , x l , t ) = m j | k , l i ( x k , x l , t ) t ,
The 95% confidence limits for the current annual increment:
l c q m j i ( t , 0.025 ) = q m j i ( t , 0.025 ) t ,   u c q m j i ( t , 0.975 ) = q m j i ( t , 0.975 ) t ,
The 95% confidence limits for the mean annual increment:
l m q m j i ( t , 0.025 ) = q m j i ( t , 0.025 ) t ,   u m q m j i ( t , 0.975 ) = q m j i ( t , 0.975 ) t ,
The relative increment:
r m j i ( t ) = 1 m j i ( t ) m j i ( t ) t   or   r m j | k i ( x k , t ) = 1 m j | k i ( x k , t ) m j | k i ( x k , t ) t ,   r m j | k , l i ( x k , x l , t ) = 1 m j | k , l i ( x k , x l , t ) m j | k , l i ( x k , x l , t ) t ,
The growth inflection point, a, is obtained by solving the following equation:
2 m j i ( a ) t 2 = 0   or   2 m j | k i ( x k , a ) t 2 = 0 ,   2 m j | k , l i ( x k , x l , a ) t 2 = 0 ,
where 2 m j i ( t ) t 2 < 0 and 2 m j i ( t ) t 2 > 0 are satisfied on the increasing and decreasing sides of point a.

2.2. Study Area and Data

Pine (Pinus sylvestris), spruce (Picea abies), and silver birch (Betula pendula Roth and Betula pubescens Ehrh.) tree stands dominate Lithuanian forests (Lithuanian Statistical Yearbook of Forestry, 2009) and grow at the Arenosol and Podzol forest sites. All the data were collected between 1983 and 2020 across the municipality of Kazlų Rūda in Lithuania. The mean temperatures vary from −16.4 °C in winter to +22° in summer. Precipitation is distributed throughout the year, although it is predominantly in the summer, averaging approximately 680 mm a year. From 1983 to 1987, 53 permanent experimental plots were installed in the forests of the Kazlų Rūda region [11]. Each sample plot consisted of about 0.16–0.72 ha and was remeasured several times, from 1 to 6, at 2- to 37-year intervals. In terms of the regeneration regime, all the plots vary between being naturally regenerated and artificially regenerated and are located in pure or mixed stands. The distribution of the permanent experimental plots can be characterized as follows: pine forests (Pinus sylvestris), 63.8%; spruce (Picea abies), 30.2%; silver birch (Betula pendula Roth and Betula pubescens Ehrh.), 5.8%; and others, 0.2 percent. The age of the ith tree (ranging from all trees to the 10th) in the first measurement cycle was recorded by counting its growth rings in the growth core (for even-aged stands, from the documented records), and the ages of the remaining trees were obtained using the arithmetic mean. The position accuracy of the plane coordinates was 1 dcm, and the diameter measurements were performed by rounding to the nearest 1 mm. The height of approximately every 5th tree was measured, and the height measurements were performed with an accuracy of approximately 1 dcm. Given the two-level measurements, two different datasets were used to obtain the parameter estimates: the first (48 plots; 39,437 mixed-species trees) was used to estimate the fixed-effect parameters for the tree diameter and the potentially available area; see Equations (1) and (2). The second (48 plots; 8604 mixed-species trees) was used to estimate the fixed-effect parameters for the tree height (Equation (3)) to calibrate the random effects (Equation (A2)) and to estimate the correlation matrix of the three-dimensional normal copula by maximizing the pseudo maximum likelihood function defined by Equation (A22). The validation set consisted of 15 plots with 5- and 15-year remeasurement cycles (a total of 2329 tree measurements), in which age, diameter, potentially available area, and height were determined for each measured tree. First of all, a set of validation data was created using only plots with three remeasurement cycles where the time difference between the first and second remeasurement cycles is 5 years and that between the first and third remeasurement cycles is 15 years. If more remeasurement cycles were performed in these validation plots, then the remaining remeasurement cycles of the plot were directed to the evaluation set. The summary of the measurements is presented in Table 1.

3. Results and Discussion

3.1. Parameter Estimates

The fixed-effect parameters of the joint copula-type three-dimensional growth model defined by the cumulative distribution function (A4) are estimated using the pseudo maximum likelihood procedure. In the first step, the fixed-effect parameter vector Ɵ was estimated separately for diameter, the potentially available area, and the height by maximizing an approximated maximum log-likelihood function and using discrete measurements of the tree diameter (x1), the potentially available area (x2), and the height (x3) { ( x 11 i , x 21 i , x 31 i ) , ( x 12 i , x 22 i , x 32 i ) , , ( x 1 n i i , x 2 n i i , x 3 n i i ) } at discrete times (ages) { t 1 i , t 2 i , , t n i i } [9]. The measurements of all the observed trees were used to estimate the parameters of the stochastic differential Equations (1) and (2), which determine the growth of the tree diameter and potentially available area, while the parameters of the tree height (Equation (3)) were estimated using a smaller part of the trees that had height measurements (see Table 1). The correlation matrix P of the three-dimensional normal copula-type probability density function defined by (A12) was estimated by maximizing the pseudo maximum likelihood function (A22). The parameter estimates of the stochastic differential Equations (1)–(3) and their standard errors and the parameters of the dependencies for the copula-type probability density function are presented in Table 2 and Table 3, respectively. All the parameter estimates are significant (p < 0.05).

3.2. Analysis of Tree Growth

The foundational aspects of the stochastic differential equation framework of tree growth in a given forest stand can be viewed through the diffusion processes, establishing the one-dimensional probability density functions of the corresponding tree size component and their dependency structure. The functions defined in Equation (7) represent the mean cumulative growth in diameter, the potentially available area, and the height achieved by a tree of a given species at any age in the particular forest stand with growth affecting factors such as artificial thinning, site conditions, changes in temperature, precipitation, or atmospheric CO2, which are included in the form of random effects. Environmental change has a decisive effect on the growth of the diameter, height, and density (potentially available area), and all these variables also affect each other in the complex form of feedback.
The functions defined in Equations (13) and (14) represent the mean cumulative growth in diameter, the potentially available area, and the height reached by trees of a given species, depending on the values of age and other tree size components in a given forest stand. In this study, we examined in detail the interaction of the components of tree size using the fixed-effect parameter mode (random effects set to zero) (shown in Figure 2), representing three tree species: pine, spruce, and birch. As shown in Figure 2(p1,p2,s1,s2,b1,b2), tree diameter growth is marginally affected by tree potentially available area, but in contrast, diameter growth is significantly affected by height changes. Changes in the height of all tree species are strongly affected by changes in diameter (see Figure 2(p5,s5,b5)), but the potentially available area has a small influence on changes in height (see Figure 2(p6,s6,b6)).
The forestry literature has established that the growth curves of fast-growing forests do not have an inflection point (apparently due to sustained growth) and express a convex shape [26]. Equation (7) has a sigmoidal or convex form and asymptotes. Unfortunately, the growth equation (Equation (7)) does not necessarily have an inflection point, and due to accelerated growth, the inflection point disappears and the growth process takes the form of a convex curve. Forest statisticians have frequently used sigmoidal functions such as logistic, Gompertz, Richards, etc. [27,28]. Graphically, the inflection point can be represented by plotting the second derivative of the diameter, potentially available area, and height growth (Equation (7) (growth acceleration)). The point at which the acceleration curve crosses the abscissa axis indicates the inflection point. Figure 3 illustrates that the growth acceleration curves of the tree diameter, potentially occupied area, and height cross the x-axis when both the fixed- and the mixed-effect modes are used, and the random effects are calibrated using Equation (A2) for all the plots in the validation dataset (among them, 15 plots have pine species trees, 9 have spruce species trees, and 13 have birch species trees).
The forms of the acceleration versus the time equations for both the fixed- and the mixed-effect modes in Figure 3 are in perfect agreement. From the Figure 3(pf1,pm3), it can be seen that the acceleration in the growth of the diameter and height of the pine trees is negative everywhere, meaning that the growth curve is convex and that an inflection point does not exist. The remaining two Figure 3(pm1,pf1), represent the growth acceleration curves of diameter (mixed effect) and height (fixed effect), respectively, of the pine species, showing that the inflection point occurs at a young age and not necessarily in all plots. It would be reasonable to assert that the soils and other environmental conditions in the Kazlų Rūda district are the most suitable for the growth of pine trees. Figure 3(sf1,sf3,sm1,sm3), for spruce trees, show that the diameter and height growth curves have a single inflection point. Additionally, the diameter and height growth acceleration curves for the spruce trees show that there exist two points at which the third derivative of the growth curves is equal to zero; these are referred to as the maximum acceleration point and maximum deceleration point. The maximum acceleration point occurs at a very young age. The diameter growth acceleration curves for the birch trees in Figure 3(pf1,pm1) show that there is a single inflection point, and maximum acceleration and deceleration points. For the birch trees, the height growth acceleration curves presented in Figure 3(pf3,pm3) show that there is no tipping point, and the height growth curve is influenced only by the point of maximum deceleration. The growth acceleration curves of the potentially available area for all the species of trees shown in Figure 3(pf2,sf2,bf2,pm2,sm2,bm2) show that there is an inflection point, and the height growth curve is influenced by the points of maximum acceleration and deceleration.
The relationships between the current and mean annual increments of the tree diameter and the height against the age of the forest stand are illustrated in Figure 4. As seen in Figure 4, the age of a stand exerts a strong influence on the current and mean annual diameter and height increments. From Figure 4, it can be seen that the peaks of the mean annual diameter and height increments are reached even later than those of the current annual increments. The highest current annual increment in tree diameter is reached at about 20–30 years for pine trees (see Figure 4(pf1,pm1)), 30–40 years for spruce trees (see Figure 4(sf1,sm1)), and 10–20 years for birch trees (see Figure 4(bf1,bm1)). The highest current annual increment in tree height occurs until the age of 10 years for pine trees (see Figure 4(pf2,pm2)) 30–40 years for spruce trees (see Figure 4(sf2,sm2)), and until the age of 10 years for birch trees (see Figure 4(bf2,bm2)), respectively. The current annual diameter increment becomes equal to the mean annual increment at 30–50 years of age for pine trees (see Figure 4(pf1,pm1)), 50–70 years of age for spruce trees (see Figure 4(sf1,sm1)), and 20–30 years of age for birch trees (see Figure 4(bf1,bm1)). The increments in diameter and height for birch trees show the highest values for the current annual increment at 1.2 cm (see Figure 4(bf1,bm1)) and 89 cm (see Figure 4(bf2,bm2)), respectively, and at the peak ages of 20 years and 10 years, respectively.
The age at which mean annual increment reaches its peak can be obtained by modifying the relative increment (Equation (19)) by multiplying it by time, t; then, the peak point can be visualized by the point at which this modified curve crosses the line y = 1. Figure 5 visualizes the curve of the relative increment defined by Equation (19) multiplied by time, t, together with the signal line y = 1 for the mixed-effect mode. As can be seen from Figure 5, trees of the pine and birch species reach their peaks of mean annual increment in diameter and height at almost the same age in all the plots in the validation set (see Figure 4(bm2) and Figure 5(pm1,bm1,pm2)). Thus, for the pine trees, the peak of mean annual increment in diameter is reached at approximately 50 years (see Figure 4(pm1)), and in height at approximately 40 years (see Figure 4(pm2)), while for the birch trees, the peak of the mean annual increment in diameter is reached at approximately 27 years (see Figure 4(bm1)), and in height at approximately 30 years (see Figure 4(bm2)). For the spruce trees, the peak of mean annual increment in diameter and height in the investigated plots had values in the intervals of 60–70 and 50–65 (see Figure 4(sm1,sm2)), respectively. The differences in the peaks of the mean annual increment between the plots can be explained by the fact that the conditions of the growing region are unfavorable for the growth of spruce trees, as demonstrated previously in Figure 2 and Figure 4. It should be noted that trees of all species reach their peak mean annual increment in height earlier than in diameter.
It makes sense to link the current annual increment in the tree diameter or height growth (see Equation (15)) to the growth of the diameter or height (see Equations (7), (13) and (14)). Figure 6 shows the linkage of the tree diameter or height with the current increment in the diameter and height growth. Figure 6(p1,s1,b1) show that the current increment in the pine tree diameter reaches its growth peak when the tree diameter reaches 10 cm, and the diameters of the spruce and birch trees reach their current increment peaks when the tree diameter is in the range of 5–15 cm. Likewise, from Figure 6(p2,s2,b2), it can be observed that the current height increments of the pine and spruce trees reach their peak when the tree height reaches 5 m, and the height of the birch tree reaches its current height increment peak when the height of the tree is in the range of 5–15 m.
In this study, using the stochastic differential Equations (1)–(3) and the three-dimensional copula-type distribution function, the three curves for diameter, the potentially available area, and the height were newly derived; for the first curve, the independent variable was only time (see Equation (7)); for the second curve, the independent variables were time and another tree size variable (see Equation (13)); finally, for the third curve, the independent variables were time and the other two remaining tree size variables (see Equation (14)). The goodness of fit of the three presented curves discussed above for modeling the tree diameter, potentially available area, and height is evaluated using statistical measures, and the results are presented in Table 4 using a validation dataset. The statistical measures presented in Table 4 show that the proposed models best corresponded to the observed validation dataset of the pine trees in uneven-aged mixed-species stands. The predictions of the pine tree heights showed stable and high statistical measures; for example, the mean bias did not exceed 0.02 m (−1.03%); the absolute mean bias did not exceed 1.36 m (7.52%); the root mean square error did not exceed 1.85 m (9.57%); and the coefficient of determination reached to as high as 88.2%. The pine tree diameter predictions showed slightly lower statistical measures, and the tree potentially occupied area predictions showed a lower performance on the statistical measures compared to the tree height predictions. It should be noted that the predictions of the diameter, height, and potentially available area for spruce trees generated the lowest statistical measures among the tree species. Additionally, it can be stated that the statistical measures of the predictions of the potentially available area were the lowest compared to the predictions of the height and diameter, and the additional variable diameter or height increased the statistical measures slightly. Furthermore, with the introduction of an additional independent variable, the tree potentially available area, into the tree height or diameter curve, the increase in the statistical measures is negligible, and in contrast, the inclusion of the tree diameter and height in the relationship of the potentially available area does not show a significant effect.
B = 1 n i = 1 n ( y i y i ) ,   ( % B = 1 n i = 1 n y i y i y i 100 ) ,
A B = 1 n i = 1 n | y i y i | ,   ( % B = 1 n i = 1 n | y i y i y i | 100 ) ,
R M S E = 1 n i = 1 n ( y i y i ) 2 ,   ( % R M S E = 1 n i = 1 n ( y i y i y i ) 2 100 ) ,
R 2 = 1 i = 1 n ( y i y i ) 2 i = 1 n ( y i y ¯ ) 2 .
Table 5 summarizes the statistical measures of the stand mean diameter, potentially available area, and height curves (7), (13), and (14). All of the mean diameter, potentially available area, and height models of a stand used, as defined by Equations (7), (13), and (14), produced predictions that showed a significantly higher level of statistical measures compared to the accuracy of individual-tree predictions assessed by the same statistical measures and presented in Table 4. The best statistical measures were demonstrated by model (14), defining the mean height of the pine trees in the stand as a function of the average age of the stand, the average diameter of the pine trees in the stand, and the average potentially available area of the pine trees in the stand. The predictions for the mean height of the pine trees in the stand demonstrated that the statistical measure for mean bias did not exceed −0.17 m (−0.86%), absolute mean bias did not exceed 0.39 m (1.83%), root mean square error did not exceed 0.58 m (2.72%), and the coefficient of determination reached as high as 98.2%.

3.3. Evolution of the Number of Trees

Environmental resources can only satisfy a limited number of trees in a stand; therefore, the number of trees gradually decreases with increasing tree diameter and height, depending on the intensity of interspecific competition. To model the number of trees per hectare during the development of a stand, the link between the number of trees per hectare and the mean potentially available area, diameter, and height was determined using the mean trend Equations (7), (13), and (14), as follows:
N i ( t | t i n , x 2   i n , θ ^ ,   ϕ p i ^ ) = ɣ p + k s i 10 , 000 . x p f p i ( x p , t | t i n , x 2   i n , θ ^ ,     ϕ p i ^ ) d x p ,   i = 1 ,   2 ,   ,   K ,
N i ( t ,   x d | t i n , x d   i n , x p   i n , θ ^ ,   ϕ d i ^ ,   ϕ p i ^ , ρ d p ^ ) = ɣ p + k s i 10 , 000 . x p p i f p | d i ( x p , t | x d , t i n , x d   i n , x p   i n , θ ^ ,     ϕ d i ^ ,     ϕ p i ^ , ρ d p ^ ) d x p ,   i = 1 ,   2 ,   ,   K ,
N i ( t ,   x h | t i n , x p   i n , x h   i n , θ ^ ,   ϕ p i ^ ,   ϕ h i ^ , ρ p h ^ ) = ɣ p + k s i 10 , 000 . x p p i f p | d i ( x p , t | x h , t i n , x p   i n , x h   i n , θ ^ ,     ϕ p i ^ ,     ϕ h i ^ , ρ p h ^ ) d x p ,   i = 1 ,   2 ,   ,   K ,
N i ( t , x d ,   x h | t i n , x d   i n , x p   i n , x h   i n θ ^ ,   ϕ d i ^ ,   ϕ p i ^ , ϕ h i ^ , P ^ ) = ɣ p + k s i 10 , 000 . x p p i f p | d h i ( x p , t | x d , x h , t i n , x d   i n , x p   i n , x h   i n , θ ^ , ϕ d i ^ ,   ϕ p i ^ ,     ϕ h i ^ , P ^ ) d x p , i = 1 , 2 , , K ,
where K is the number of observed plots in the validation dataset; k s l is the part of the stand occupied by the tree species concerned ( k a l l l = 1 , 0 k p i n e l 1 , 0 k s p r u c e l 1 , 0 k b i r c h l 1 ); the estimates of the fixed-effect parameters θ ^ are taken from Table 2; the estimates of the dependency parameters   P ^   , ρ d p ^ , and ρ p h ^ are taken from Table 3; the random effects ϕ d i ^ ,   ϕ p i ^ ,   ϕ h i ^ are calibrated by Equation (A2) using the validation dataset; for the prediction mode, t i n = 4 , x d   i n = 0.1 , x p   i n =   δ ^   , x h   i n = 0.4 (where   δ ^ was taken from Table 2); and for the forecast mode, t i n , x d   i n , x p   i n   , x h   i n are the average values of the first cycle measurements of age, diameter, potentially available area, and height. The evolution of the number of trees per hectare of a particular stand against age for all species and pine, spruce, and birch trees and four different stands is shown in Figure 7. The results of the statistical measures of the number of trees per hectare predictions are presented in Table 6. The statistical measures of the number of trees per hectare for the forecast mode are presented in Table 7. Figure 7 shows that the model curve of the number of trees per hectare defined by Equation (25) corresponds well to the measurement data from the validation dataset. The best fitness of the predictions of the number of trees per hectare was confirmed by all the statistical measures for the trees of the pine species, while the aforementioned statistical measures were the worst for the trees of the spruce species. Table 6 and Table 7 demonstrate the higher goodness-of-fit statistical measures compared to those achieved using the generalized algebraic difference approach, which is the one used traditionally [29]. The statistical measures computed in Table 6 and Table 7 for the number of trees per hectare models, with the additional explanatory variables for diameter and height, indicated a slight improvement in prediction accuracy, with a lower performance of the first model, as expected. Therefore, as demonstrated in Table 6 and Table 7, the maximum impact on the tree density dynamics demonstrated the tree potentially available area and age.
The number of trees per hectare varies with the mean diameter, height, and volume of the stand. Stands of small mean diameter, height, or volume have a large number of trees, while stands of large mean diameter, height, or volume have relatively few [30]. The relationship between the mean diameter, height, and volume (increasing over age) and the number of live trees per hectare (decreasing over age) can been described by means of a “limiting relationship” [31]. The intensity of competition between trees varies with the initial number of trees per hectare, as shown in Figure 7. The process of self-thinning (decay) in the different stands shown in Figure 7 indicates that the stands formed with different initial densities do not thin to the same final number of trees per hectare (see Figure 7a–d). The process of the current rate of the decay of the number of trees per hectare, expressed as the ratio of the current annual decay (Equation (15)) to the number of trees per hectare (Equation (25)), is shown in Figure 8. From Figure 8a–d, we can state that the decay process drops to 1% for the pine species trees at approximately 100 years; for the spruce species trees at 80 years; and finally at 70 years for the birch species trees. It should be noted that the soil characteristics of this region are not favorable for spruce trees.

4. Conclusions

Forest system data are characterized by important features such as dynamism, uncertainty, stochasticity, deviation of variables from the normal distribution, multicollinearity, and false correlations, which complicate the modeling and prediction of forest systems. Generally, forestry researchers treat probabilistic modeling as the mathematical study of an observed random variable, such as tree diameter, height, etc. The theory of stochastic processes enables the description of the evolution of random variables with respect to time, which plays an important role in the study of various growth processes in a forest stand. Stochastic processes, in a sense, could be referred to as the dynamic part of probability theory. Our previous works used the multi-dimensional stochastic process paradigm, where for parameter estimates a lot of computing resources are consumed. Linking one-dimensional stochastic processes according to the copula function reduced the volume of computer calculations and made it possible to analyze the growth processes of various trees, stands and their connections in the form of nonlinear equations.

Author Contributions

Conceptualization, P.R. and E.P.; methodology, P.R. and E.P.; software, P.R.; validation, P.R. and E.P.; formal analysis, P.R.; resources, E.P.; data curation, E.P.; writing—original draft preparation, P.R.; writing—review and editing, P.R. and E.P.; visualization, P.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Original data presented in the study are included in the main text, and further inquiries can be directed to the corresponding author.

Acknowledgments

We thank the Edmundas Petrauskas research group for their help with sampling in the field. The authors would like to express their appreciation for the support from the Lithuanian Association of Impartial Timber Scalers. Thanks are also due to anonymous reviewers for their constructive criticism.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Let F j i ^ ( x j , t ) , j = d, p, h (diameter, potentially available area, height), and i = 1, …, K (where K is the number of the plots from the validation dataset) denote the cumulative distribution function of the Gompertz-type diffusion processes, defined by Equations (1)–(3), which take the following forms:
F j i ^ ( x j , t ) = ɣ j ^ + f j i ^ ( x j , t ) d x j .
It should be noted that in Equation (A1), the probability density functions f j i ^ ( x j , t )  j = d, p, h; i = 1, …, K incorporate the values of the fixed-effect parameter θ ^ estimated using the approximated maximum likelihood procedure (see Table 2), and the random effects for each validation plot are calibrated in the form [32,33]:
ϕ j i ^ = a r g m a x ( ϕ j i ) ( j = 1 m i l n ( f j i ( x j , t | θ ^ , ϕ j i ) ) + k = 1 4 l n ( p ( ϕ j i | σ ^ j 2 ) ) ) ,
where mi is the number of observed trees in the ith validation plot, p ( ϕ j i | τ ^ j 2 ) is the normal probability density function with 0 mean and variance τ j 2 , and, finally,
f j i ^ ( x j , t ) = f j i ( x j , t | θ ^ , ϕ j i ^ ) .
The three-dimensional normal copula cumulative distribution functions i = 1, …, K with correlation matrix P are defined by
C 3 i ( u 1 , u 2 , u 3 ; P , t ) = Φ 3 ( Φ 1 ( u 1 ) , Φ 1 ( u 2 ) , Φ 1 ( u 3 ) ; P ) ,
where
u 1 = F d i ^ ( x d , t ) ,   u 2 = F p i ^ ( x p , t ) ,   and   u 3 = F h i ^ ( x h , t )
Φ 3 ( x d , x p , x h ; P ) = x d x p x h φ 3 ( z 1 , z 2 , z 3 ; P ) d z 1 d z 2 d z 3 ,
φ 3 ( x d , x p , x h ; P ) = e w 2 ( ρ 12 2 + ρ 13 2 ± ρ 23 2 2 ρ 12 ρ 13 ρ 23 ) 2 π 3 2 1 ( ρ 12 2 + ρ 13 2 + ρ 23 2 ) + 2 ρ 12 ρ 13 ρ 23 ,
w = x d 2 ( ρ 23 2 1 ) + x p 2 ( ρ 13 2 1 ) + x h 2 ( ρ 12 2 1 ) + 2 ( x d x p ( ρ 12 ρ 13 ρ 23 ) + x d x h ( ρ 13 ρ 12 ρ 23 ) + x p x h ( ρ 23 ρ 12 ρ 13 ) ) ,
Φ ( x ) = x φ ( z ) d z ,
φ ( x ) = 1 2 π e x 2 2 ,
P = ( 1 ρ 12 ρ 13 ρ 12 1 ρ 23 ρ 13 ρ 23 1 ) .
The three-dimensional normal copula probability density functions c 3 i ( u 1 , u 2 , u 3 ; P , t )  i = 1, …, K take the following forms:
c 3 i ( u 1 , u 2 , u 3 ; P , t ) = 3 u 1 u 2 u 3 C 3 i ( u 1 , u 2 , u 3 ; P , t ) = φ 3 ( Φ 1 ( u 1 ) , Φ 1 ( u 2 ) , Φ 1 ( u 3 ) ; P ) φ ( Φ 1 ( u 1 ) ) φ ( Φ 1 ( u 1 ) ) φ ( Φ 1 ( u 1 ) ) = 1 1 ( ρ 12 2 + ρ 13 2 + ρ 23 2 ) + 2 ρ 12 ρ 13 ρ 23 e w 2 ( ρ 12 2 + ρ 13 2 ± ρ 23 2 2 ρ 12 ρ 13 ρ 23 1 ) , w = x 1 2 ( 2 ρ 12 ρ 13 ρ 23 ρ 12 2 ρ 13 2 ) + x 2 2 ( 2 ρ 12 ρ 13 ρ 23 ρ 12 2 ρ 23 2 ) + x 3 2 ( 2 ρ 12 ρ 13 ρ 23 ρ 13 2 ρ 23 2 ) + 2 ( x 1 x 2 ( ρ 12 ρ 13 ρ 23 ) + x 1 x 3 ( ρ 13 ρ 12 ρ 23 ) + x 2 x 3 ( ρ 23 ρ 12 ρ 13 ) ) ,
where x 1 = Φ 1 ( u 1 ) , x 2 = Φ 1 ( u 2 ) , and x 3 = Φ 1 ( u 3 ) .
The joint three-dimensional copula-type probability density functions f 3 c i ( x d , x p , x h , t )   i = 1, …, K are given as
f 3 c i ( x d , x p , x h : P , t ) = c 3 i ( F d i ^ ( x d , t ) , F p i ^ ( x p , t ) , F h i ^ ( x h , t ) ; P , t ) f d i ^ ( x d , t ) f p i ^ ( x p , t ) f h i ^ ( x h , t ) .
The two-dimensional normal copula cumulative distribution functions C 2 i ( u j , u k ; ρ j k , t ) , jk; j, k { d , p , h } ; i = 1, …, K with correlation coefficient ρ j k are defined analogously:
C 2 i ( u j , u k ; ρ j k , t ) = Φ 2 ( Φ 1 ( u j ) , Φ 1 ( u k ) ; ρ j k )
where
u j = F j i ^ ( x j , t ) ,   and   u k = F k i ^ ( x k , t ) ,
Φ 2 ( x j , x k , P ) = x j x k φ 2 ( z 1 , z 2 ; ρ j k ) d z 1 d z 2 ,
φ 2 ( x j , x k ; ρ j k ) = 1 2 π 1 ρ j k 2 e x j 2 2 ρ j k x j x k + x k 2 2 ( 1 ρ j k 2 ) .
The two-dimensional normal copula probability density functions c 2 i ( u j , u k ; ρ j k , t ) , jk; j, k { d , p , h } ; i = 1, …, K take the following form:
c 2 i ( u j , u k ; ρ j k , t ) = 2 u j u k C 2 i ( u j , u k ; ρ j k , t ) = φ 2 ( Φ 1 ( u j ) , Φ 1 ( u k ) ; ρ j k ) φ ( Φ 1 ( u j ) ) φ ( Φ 1 ( u k ) ) .
The joint two-dimensional copula-type probability density functions f 2 c i ( x j , x k ; ρ j k , t ) , jk; j, k { d , p , h } ; i = 1, …, K are given as
f 2 c i ( x j , x k ; ρ j k , t ) = c 2 i ( u j , u k ; ρ j k , t ) f j i ^ ( x j , t ) f k i ^ ( x k , t ) .
Therefore, the conditional probability density functions of X l i ( t ) , i = 1, …, K; l = 1, 2, 3 at a given ( X j ( t ) = x j , X k ( t ) = x k ) , j, k ≠ l, are defined as
f d | p , h i ( x d , t | x p , x h ) = f 3 c i ( x d , x p , x h ; P ^ , t ) f 2 c i ( x p , x h ; ρ p h ^ , t ) , f p | d , h i ( x p , t | x d , x h ) = f 3 c i ( x d , x p , x h ; P ^ , t ) f 2 c i ( x d , x h ; ρ d h ^ , t ) , f h | d , p i ( x h , t | x d , x p ) = f 3 c i ( x d , x p , x h ; P ^ , t ) f 2 c i ( x d , x p ; ρ d p ^ , t ) .
The conditional probability density functions of X j i ( t ) , i = 1, …, K; j = d, p, h at a given ( X k ( t ) = x k ) , kj, are defined as:
f d | p ( x d , t | x p ) = f 2 c i ( x d , x p ; ρ d p ^ , t ) f p i ^ ( x p , t ) , f d | h ( x d , t | x h ) = f 2 c i ( x d , x h ; ρ d h ^ , t ) f h i ^ ( x h , t ) , f p | d ( x p , t | x d ) = f 2 c i ( x d , x p ; ρ d p ^ , t ) f d i ^ ( x d , t ) , f p | h ( x p , t | x h ) = f 2 c i ( x p , x h ; ρ p h ^ , t ) f h i ^ ( x h , t ) , f h | d ( x h , t | x d ) = f 2 c i ( x d , x h ; ρ d h ^ , t ) f d i ^ ( x d , t ) , f h | p ( x h , t | x p ) = f 2 c i ( x p , x h ; ρ p h ^ , t ) f p i ^ ( x p , t ) .
The pseudo maximum log-likelihood function for the copula-type probability density function c 3 i is represented as
L L ( ρ 12 , ρ 13 , ρ 23 ) = i = 1 M j = 1 n i l n ( c 3 i ( Φ 1 ( F d i ^ ( x 1 j i , t j i ) ) , Φ 1 ( F p i ^ ( x 1 j i , t j i ) ) , Φ 1 ( F h i ^ ( x 1 j i , t j i ) ) | ( ρ 12 , ρ 13 , ρ 23 ) , t j i ) )

References

  1. Zeide, B. Analysis of Growth Equations. For. Sci. 1993, 39, 594–616. [Google Scholar] [CrossRef]
  2. Petrauskas, E.; Rupšys, P.; Narmontas, M.; Aleinikovas, M.; Beniušienė, L.; Šilinskas, B. Stochastic Models to Qualify Stem Tapers. Algorithms 2020, 13, 94. [Google Scholar] [CrossRef] [Green Version]
  3. Rupšys, P.; Petrauskas, E. Symmetric and Asymmetric Diffusions through Age-Varying Mixed-Species Stand Parameters. Symmetry 2021, 13, 1457. [Google Scholar] [CrossRef]
  4. Suzuki, T. Forest transition as a stochastic process (I). J. Jpn. For. Sci. 1966, 48, 436–439. [Google Scholar]
  5. Sloboda, B. Kolmogorow–Suzuki und die stochastische Differentialgleichung als Beschreibungsmittel der Bestandesevolution. Mitt Forstl Bundes Vers. Wien 1977, 120, 71–82. [Google Scholar]
  6. Garcia, O. Modelling stand development with stochastic differential equations. In Mensuration for Management Planning of Exotic Forest Plantations; Elliot, D.A., Ed.; New Zealand Forest Service, Forest Research Institute Symposium: Rotorua, New Zealand, 1979; 20, pp. 315–333. [Google Scholar]
  7. Tanaka, K. A stochastic model of height growth in an even-aged pure forest stand—why is the coefficient of variation of the height distribution smaller than that of the diameter distribution. J. Jpn. For. Soc. 1988, 70, 20–29. [Google Scholar]
  8. Rennolls, K. Forest height growth modelling. For. Ecol. Manag. 1995, 71, 217–225. [Google Scholar] [CrossRef]
  9. Rupšys, P. New insights into tree height distribution based on mixed-effects univariate diffusion processes. PLoS ONE 2016, 11, e0168507. [Google Scholar] [CrossRef]
  10. Narmontas, M.; Rupšys, P.; Petrauskas, E. Construction of Reducible Stochastic Differential Equation Systems for Tree Height–Diameter Connections. Mathematics 2020, 8, 1363. [Google Scholar] [CrossRef]
  11. Rupšys, P.; Petrauskas, E. Analysis of Longitudinal Forest Data on Individual-Tree and Whole-Stand Attributes Using a Stochastic Differential Equation Model. Forests 2022, 13, 425. [Google Scholar] [CrossRef]
  12. Rupšys, P. Stochastic Mixed-Effects Parameters Bertalanffy Process, with Applications to Tree Crown Width Modeling. Math. Probl. Eng. 2015, 2015, 375270. [Google Scholar] [CrossRef] [Green Version]
  13. Rupšys, P. Understanding the Evolution of Tree Size Diversity within the Multivariate nonsymmetrical Diffusion Process and Information Measures. Mathematics 2019, 7, 761. [Google Scholar] [CrossRef] [Green Version]
  14. Petrauskas, E.; Bartkevičius, E.; Rupšys, P.; Memgaudas, R. The use of stochastic differential equations to describe stem taper and volume. Baltic For. 2013, 19, 43–151. [Google Scholar]
  15. Narmontas, M.; Rupšys, P.; Petrauskas, E. Models for Tree Taper Form: The Gompertz and Vasicek Diffusion Processes Framework. Symmetry 2020, 12, 80. [Google Scholar] [CrossRef]
  16. Rupšys, P. Generalized fixed-effects and mixed-effects parameters height–diameter models with diffusion processes. Int. J. Biomath. 2015, 8, 1550060. [Google Scholar] [CrossRef]
  17. Rupšys, P. Modeling Dynamics of Structural Components of Forest Stands Based on Trivariate Stochastic Differential Equation. Forests 2019, 10, 506. [Google Scholar] [CrossRef] [Green Version]
  18. Rupšys, P.; Petrauskas, E. On the Construction of Growth Models via Symmetric Copulas and Stochastic Differential Equations. Symmetry 2022, 14, 2127. [Google Scholar] [CrossRef]
  19. Wang, M.; Upadhyay, A.; Zhang, L. Trivariate distribution modeling of tree diameter, height, and volume. For. Sci. 2010, 56, 290–300. [Google Scholar]
  20. Sklar, M. Fonctions de repartition an dimensions et leurs marges. Publ. Inst. Statist. Univ. Paris 1959, 8, 229–231. [Google Scholar]
  21. Rupšys, P.; Petrauskas, E. Evolution of Bivariate Tree Diameter and Height Distribution via Stand Age: Von Bertalanffy Bivariate Diffusion Process Approach. J. For. Res. 2019, 24, 16–26. [Google Scholar] [CrossRef]
  22. Rupšys, P.; Narmontas, M.; Petrauskas, E. A Multivariate Hybrid Stochastic Differential Equation Model for Whole-Stand Dynamics. Mathematics 2020, 8, 2230. [Google Scholar] [CrossRef]
  23. Rupšys, P.; Petrauskas, E. A Linkage among Tree Diameter, Height, Crown Base Height, and Crown Width 4-variate Distribution and Their Growth Models: A 4-variate Diffusion Process Approach. Forests 2017, 8, 479. [Google Scholar] [CrossRef] [Green Version]
  24. Ishihara, M.I.; Konno, Y.; Umeki, K.; Ohno, Y.; Kikuzawa, K. A new model for size-dependent tree growth in forests. PLoS ONE 2016, 11, e0152219. [Google Scholar] [CrossRef] [Green Version]
  25. Itô, K. Stochastic integral. Proc. Imp. Acad. 1944, 20, 519–524. [Google Scholar] [CrossRef]
  26. Yuancai, L.; Marques, C.; Macedo, F. Comparison of Schnute’s and Bertalanffy-Richards’ growth functions. For. Ecol. Manag. 1997, 96, 283–288. [Google Scholar] [CrossRef]
  27. Monti, C.A.U.; Oliveira, R.M.; Roise, J.P.; Scolforo, H.F.; Gomide, L.R. Hybrid Method for Fitting Nonlinear Height–Diameter Functions. Forests 2022, 13, 1783. [Google Scholar] [CrossRef]
  28. Rodrigo, M.; Zulkarnaen, D. Mathematical Models for Population Growth with Variable Carrying Capacity: Analytical Solutions. AppliedMath 2022, 2, 466–479. [Google Scholar] [CrossRef]
  29. Gómez-García, E.; Crecente-Campo, F.; Tobin, B.; Hawkins, M.; Nieuwenhuis, M.; Diéguez-Aranda, U. A dynamic volume and biomass growth model system for even-aged downy birch stands in south-western Europe. Forestry 2014, 87, 165–176. [Google Scholar] [CrossRef]
  30. Clutter, J.L.; Bennett, F.A. Diameter Distributions in Old-Field Slash Pine Plantations; Georgia Forest Research Council: Atlanta, GA, USA, 1965; Volume 13, p. 9. [Google Scholar]
  31. Reineke, L.H. Perfecting a stand-density index for evenaged forests. J. Agric. Res. 1933, 46, 627–638. [Google Scholar]
  32. Rupšys, P. The use of copulas to practical estimation of multivariate stochastic differential equation mixed effects models. AIP Conf. Proc. 2015, 1684, 080011. [Google Scholar]
  33. Rupšys, P. Univariate and Bivariate Diffusion Models: Computational Aspects and Applications to Forestry. In Stochastic Differential Equations: Basics and Applications; Tony, G.D., Ed.; Nova Science Publisher’s: New York, NY, USA, 2018; pp. 1–77. [Google Scholar]
Figure 1. Tree potentially available area for randomly selected plot: 1983rd year cycle of measurement (mean age, 50.38 years); red—Scots pine trees; green—Norway spruce trees; yellow—birch trees; circles—tree position.
Figure 1. Tree potentially available area for randomly selected plot: 1983rd year cycle of measurement (mean age, 50.38 years); red—Scots pine trees; green—Norway spruce trees; yellow—birch trees; circles—tree position.
Forests 14 00012 g001
Figure 2. Relationships between components of two tree sizes at a particular age: solid black line—20 years; solid blue line—60 years; solid red line—120 years; (p1p6) pine trees; (s1s6) spruce trees; (b1b6) birch trees; (p1,s1,b1) diameter growth against potentially available area; (p2,s2,b2) diameter growth against height; (p3,s3,b3) potentially available area growth against diameter; (p4,s4,b4) potentially available area growth against height; (p5,s5,b5) height growth against diameter; (p6,s6,b6) height growth against potentially available area.
Figure 2. Relationships between components of two tree sizes at a particular age: solid black line—20 years; solid blue line—60 years; solid red line—120 years; (p1p6) pine trees; (s1s6) spruce trees; (b1b6) birch trees; (p1,s1,b1) diameter growth against potentially available area; (p2,s2,b2) diameter growth against height; (p3,s3,b3) potentially available area growth against diameter; (p4,s4,b4) potentially available area growth against height; (p5,s5,b5) height growth against diameter; (p6,s6,b6) height growth against potentially available area.
Forests 14 00012 g002
Figure 3. Acceleration of tree diameter, potentially available area, and height growth for both the fixed- and the mixed-effect modes via age: (pf1pf3,sf1sf3,bf1bf3) fixed-effect mode; (pm1pm3,sm1sm3,bm1bm3) mixed-effect mode; (pf1pf3,pm1pm3) pine trees; (sf1sf3,sm1sm3) spruce trees; (bf1bf3,bm1bm3) birch trees; first column—diameter; second column—potentially available area; third column—height.
Figure 3. Acceleration of tree diameter, potentially available area, and height growth for both the fixed- and the mixed-effect modes via age: (pf1pf3,sf1sf3,bf1bf3) fixed-effect mode; (pm1pm3,sm1sm3,bm1bm3) mixed-effect mode; (pf1pf3,pm1pm3) pine trees; (sf1sf3,sm1sm3) spruce trees; (bf1bf3,bm1bm3) birch trees; first column—diameter; second column—potentially available area; third column—height.
Forests 14 00012 g003
Figure 4. Current and mean annual increments of tree diameter and height growth for both the fixed- and the mixed-effect modes via age: (pf1,pf2,sf1,sf2,bf1,bf2) fixed-effect mode; (pm1,pm2,sm1,sm2,bm1,bm2) mixed-effect mode; (pf1,pf2,pm1,pm2) pine trees; (sf1,sf2,sm1,sm2) spruce trees; (bf1,bf2,bm1,bm2) birch trees; first column—diameter; second column—height; black—first stand; blue—second stand; red—third stand; solid line—current annual increment; dotted line—mean annual increments.
Figure 4. Current and mean annual increments of tree diameter and height growth for both the fixed- and the mixed-effect modes via age: (pf1,pf2,sf1,sf2,bf1,bf2) fixed-effect mode; (pm1,pm2,sm1,sm2,bm1,bm2) mixed-effect mode; (pf1,pf2,pm1,pm2) pine trees; (sf1,sf2,sm1,sm2) spruce trees; (bf1,bf2,bm1,bm2) birch trees; first column—diameter; second column—height; black—first stand; blue—second stand; red—third stand; solid line—current annual increment; dotted line—mean annual increments.
Forests 14 00012 g004aForests 14 00012 g004b
Figure 5. Relative increments multiplied by time of tree diameter, and height growth along with the signal line y = 1: (pm1,pm2) pine trees; (sm1,sm2) spruce trees; (bm1,bm2) birch trees; first column—diameter; second column—height; black—first stand; blue—second stand; red—third stand; dashed line—signal line; black—first stand; blue—second stand; red—third stand.
Figure 5. Relative increments multiplied by time of tree diameter, and height growth along with the signal line y = 1: (pm1,pm2) pine trees; (sm1,sm2) spruce trees; (bm1,bm2) birch trees; first column—diameter; second column—height; black—first stand; blue—second stand; red—third stand; dashed line—signal line; black—first stand; blue—second stand; red—third stand.
Forests 14 00012 g005
Figure 6. Linkage of current increments of tree diameter and height growth with mean diameter and height: (p1,p2) pine trees; (s1,s2) spruce trees; (b1,b2) birch trees; first column—diameter; second column—height; black—first stand; blue—second stand; red—third stand.
Figure 6. Linkage of current increments of tree diameter and height growth with mean diameter and height: (p1,p2) pine trees; (s1,s2) spruce trees; (b1,b2) birch trees; first column—diameter; second column—height; black—first stand; blue—second stand; red—third stand.
Forests 14 00012 g006
Figure 7. Evolution of the number of trees per hectare for three randomly selected stands: black—first stand; blue—second stand; red—third stand; (a) all trees; (b) pine trees ( k p i n e 1 = 0.97 in black, k p i n e 2 = 0.25 in blue; k p i n e 3 = 0.73 in red); (c) spruce trees ( k s p r u c e 1 = 0.05 in black, k s p r u c e 2 = 0.72 in blue; k s p r u c e 3 = 0.31 in red); (d) birch tree ( k b i r c h 1 = 0.22 in black, k b i r c h 2 = 0.16 in blue; k b i r c h 3 = 0.05 in red).
Figure 7. Evolution of the number of trees per hectare for three randomly selected stands: black—first stand; blue—second stand; red—third stand; (a) all trees; (b) pine trees ( k p i n e 1 = 0.97 in black, k p i n e 2 = 0.25 in blue; k p i n e 3 = 0.73 in red); (c) spruce trees ( k s p r u c e 1 = 0.05 in black, k s p r u c e 2 = 0.72 in blue; k s p r u c e 3 = 0.31 in red); (d) birch tree ( k b i r c h 1 = 0.22 in black, k b i r c h 2 = 0.16 in blue; k b i r c h 3 = 0.05 in red).
Forests 14 00012 g007
Figure 8. Evolution of the relative annual decay process of the number of trees per hectare for three randomly selected stands: black—first stand; blue—second stand; red—third stand; (a) all trees; (b) pine trees ( k p i n e 1 = 0.97 in black, k p i n e 2 = 0.25 in blue; k p i n e 3 = 0.73 in red); (c) spruce trees ( k s p r u c e 1 = 0.05 in black, k s p r u c e 2 = 0.72 in blue; k s p r u c e 3 = 0.31 in red); (d) birch tree ( k b i r c h 1 = 0.22 in black, k b i r c h 2 = 0.16 in blue; k b i r c h 3 = 0.05 in red).
Figure 8. Evolution of the relative annual decay process of the number of trees per hectare for three randomly selected stands: black—first stand; blue—second stand; red—third stand; (a) all trees; (b) pine trees ( k p i n e 1 = 0.97 in black, k p i n e 2 = 0.25 in blue; k p i n e 3 = 0.73 in red); (c) spruce trees ( k s p r u c e 1 = 0.05 in black, k s p r u c e 2 = 0.72 in blue; k s p r u c e 3 = 0.31 in red); (d) birch tree ( k b i r c h 1 = 0.22 in black, k b i r c h 2 = 0.16 in blue; k b i r c h 3 = 0.05 in red).
Forests 14 00012 g008
Table 1. Tree age, diameter, potentially available area, and height summary statistics for model estimation and validation datasets *.
Table 1. Tree age, diameter, potentially available area, and height summary statistics for model estimation and validation datasets *.
SpeciesDataNumber of TreesMinMaxMeanSt. Dev.Number of TreesMinMaxMeanSt. Dev.
EstimationValidation
Pinet (year)24,17612.0211.056.4926.75153125.0119.051.8820.73
d (cm)24,1760.161.019.3010.3215315.5052.8019.607.53
p (m2)24,1760.09124.1910.508.9815311.2046.819.986.48
h (m)53460.2037.9017.299.1215316.2037.5019.385.39
Sprucet (year)13,36012.0207.064.4225.276647.0101.049.8916.90
d (cm)13,3600.2072.2012.938.706643.6044.8010.935.58
p (m2)13,3600.11160.2410.168.946640.6131.157.375.24
h (m)28430.5038.012.568.526641.031.0011.845.34
Bircht (year)176112.0127.8257,4923.2213325.075.043.2410.32
d (cm)17610.9050.015.179.421335.1035.1016.126.59
p (m2)17610.33173.8210.378.871331.5132.557.594.87
h (m)3880.5031.9014.858.431337.8031.018.695.19
Allt (year)39,43712.0211.059.2526.3623297.0119.050.8119.35
d (cm)39,4370.172.2016.9510.2223293.652.816.937.98
p (m2)39,4370.09173.8210.378.9523290.6146.819.106.19
h (m)86040.2038.0015.629.1623291.037.5017.196.34
* t—age, d—diameter, p—potentially available area (the area of the polygon formed by the Voronoi diagram [3]), h—tree height.
Table 2. Parameter estimates for mixed-effect mode of the stochastic differential Equations (1)–(3).
Table 2. Parameter estimates for mixed-effect mode of the stochastic differential Equations (1)–(3).
SpeciesAΒɣσδτj
Diameter
Pine0.08780.0169−23.50980.0006-0.0162
Spruce0.09260.0291−1.55240.0102-0.0146
Birch0.27250.1060−0.21370.0337-0.0679
All0.08500.0226−7.11080.0042-0.0069
Potentially available area
Pine0.05380.0157−1.84350.00711.54780.0076
Spruce0.06960.0230−0.76300.01601.74670.0125
Birch0.06690.0207−2.33950.00881.63750.0099
All0.06170.0186−1.32600.01021.61510.0094
Height
Pine0.09030.0213−36.74860.0001-0.0021
Spruce0.09140.0264−2.17430.0062-0.0085
Birch0.20250.0577−12.03700.0032-0.0147
All0.08270.0213−13.34590.0013-0.0043
Table 3. Parameter estimates for copula probability density function (19).
Table 3. Parameter estimates for copula probability density function (19).
Species𝝆12𝝆13𝝆23
Pine0.14760.69640.0446
Spruce0.20830.87860.1512
Birch0.16120.78540.1261
All0.15130.85280.0877
Table 4. Results of statistical measures * and p-values of Student’s t-test for the individual-tree predictions using the mixed-effect mode.
Table 4. Results of statistical measures * and p-values of Student’s t-test for the individual-tree predictions using the mixed-effect mode.
Curve
(Variables)
B
(%)
AB
(%)
RMSE
(%)
R2T
p-Value
Curve
(Variables)
B
(%)
AB
(%)
RMSE
(%)
R2T
p-Value
Pine Tree DiameterPine Tree Height
Equation (7)
(t)
−0.0523
(−6.76)
3.7676
(21.31)
4.8276
(24.63)
0.58910.6716Equation (7)
(t)
0.0094
(−2.28)
2.0417
(11.65)
2.7259
(14.06)
0.74410.8921
Equation (13)
(t, p)
−0.0525
(−6.72)
3.7081
(20.99)
4.8127
(24.56)
0.59160.6693Equation (13)
(t, d)
0.0201
(−1.07)
1.3732
(7.59)
1.8684
(9.64)
0.87980.6733
Equation (13)
(t, h)
−0.0586
(−3.08)
2.4951
(13.16)
3.3515
(17.10)
0.80190.4941Equation (13)
(t, p)
0.0081
(−2.30)
2.0428
(11.67)
2.7290
(14.08)
0.74350.9040
Equation (13)
(t, p, h)
−0.0598
(−3.07)
2.4439
(12.84)
3.3119
(16.90)
0.80660.4866Equation (13)
(t, d, p)
0.0214
(−1.03)
1.3647
(7.52)
1.8553
(9.57)
0.88150.6508
Pine Tree Potentially Available AreaSpruce Tree Diameter
Equation (7)
(t)
−0.1330
(−27.75)
3.6345
(49.64)
4.8685
(48.79)
0.43580.2852Equation (7)
(t)
−0.7109
(−14.84)
3.1111
(32.20)
4.3715
(39.97)
0.37020.0001
Equation (13)
(t, d)
−0.1270
(−27.22)
3.6019
(49.05)
4.8270
(48.37)
0.44540.3034Equation (13)
(t, p)
−0.5898
(−13.70)
2.8473
(29.85)
4.0189
(36.75)
0.47020.0002
Equation (13)
(t, h)
−0.0864
(−36.59)
3.8832
(60.09)
5.9846
(63.67)
0.44450.1800Equation (13)
(t, h)
−0.2101
(−3.12)
1.4505
(14.29)
2.1234
(19.42)
0.85380.0110
Equation (13)
(t, d, h)
−0.0775
(−35.72)
3.8191
(59.12)
5.9044
(62.82)
0.45930.2229Equation (13)
(t, p, h)
−0.2032
(−3.05)
1.3917
(13.77)
2.0059
(18.34)
0.86950.0092
Spruce Tree HeightSpruce Tree Potentially Available Area
Equation (7)
(t)
−0.4554
(−13.65)
3.0009
(31.30)
4.0033
(33.82)
0.43090.0035Equation (7)
(t)
−0.0902
(−50.87)
3.6998
(5.0414)
5.0114
(68.39)
0.07370.6448
Equation (13)
(t, d)
0.0029
(−4.52)
1.4286
(15.20)
1.8768
(15.85)
0.87650.9682Equation (13)
(t, d)
−0.0370
(−47.16)
3.5327
(73.93)
4.7696
(64.70)
0.17110.8416
Equation (13)
(t, p)
−0.3581
(−13.11)
2.8746
(30.40)
3.8185
(32.26)
0.48440.0160Equation (13)
(t, h)
−0.0558
(−49.12)
3.6247
(76.41)
4.9000
(66.47)
0.12510.7691
Equation (13)
(t, d, p)
0.0024
(−4.49)
1.4145
(15.07)
1.8656
(15.76)
0.87800.9732Equation (13)
(t, d, h)
−0.0366
(−46.64)
3.5100
(73.27)
4.7422
(64.33)
0.18060.8425
Birch Tree DiameterBirch Tree Height
Equation (7)
(t)
−0.1782
(−13.61)
3.91.9
(31.24)
4.8552
(30.11)
0.45680.6739Equation (7)
(t)
−0.0788
(−5.45)
2.7530
(17.44)
3.5389
(18.93)
0.52930.7984
Equation (13)
(t, p)
−0.1535
(−13.27)
3.8440
(30.73)
4.8083
(29.79)
0.46830.7140Equation (13)
(t, d)
−0.1668
(−1.93)
1.3747
(7.92)
1.7571
(9.40)
0.88300.2772
Equation (13)
(t, h)
0.0969
(−3.92)
2.1346
(14.93)
2.7308
(16.93)
0.82820.6839Equation (13)
(t, p)
−0.0735
(−5.42)
2.7269
(17.33)
3.5634
(19.06)
0.52280.8129
Equation (13)
(t, p, h)
0.1013
(−3.87)
2.1098
(14.76)
2.7082
(16.79)
0.83100.6679Equation (13)
(t, d, p)
−0.1668
(−1.93)
1.3745
(7.92)
1.7569
(9.40)
0.88300.2771
* Statistical measures: the mean bias, B (the percentage mean bias, % B ); the absolute mean bias, AB (the percentage absolute mean bias, % A B ); the root mean square error, RMSE (the percentage root mean square error, % R M S E ); and the coefficient of determination, R2, defined by Equations (21)–(24).
Table 5. Results of statistical measures and p-values of Student’s t-test for whole-stand predictions using the mixed-effect model.
Table 5. Results of statistical measures and p-values of Student’s t-test for whole-stand predictions using the mixed-effect model.
Curve
(Variables)
B
(%)
AB
(%)
RMSE
(%)
R2T
p-Value
Curve
(Variables)
B
(%)
AB
(%)
RMSE
(%)
R2T
p-Value
Pine Tree DiameterPine Tree Height
Equation (7)
(t)
0.0137
(0.02)
0.9450
(4.23)
1.4187
(6.38)
0.93180.9489Equation (7)
(t)
−0.1262
(−0.52)
0.5926
(2.76)
0.8655
(4.08)
0.96200.3385
Equation (13)
(t, p)
−0.1319
(−0.63)
0.9479
(4.18)
1.3652
(6.14)
0.93630.5247Equation (13)
(t, d)
−0.2203
(−1.02)
0.3980
(1.85)
0.5834
(2.75)
0.98070.0160
Equation (13)
(t, h)
0.2499
(1.00)
0.5783
(2.61)
0.8827
(3.97)
0.97150.0669Equation (13)
(t, p)
−0.1513
(−0.64)
0.5889
(2.73)
0.8621
(4.07)
0.96200.2505
Equation (13)
(t, p, h)
0.1306
(0.46)
0.5367
(2.47)
0.8347
(3.75)
0. 97580.3051Equation (13)
(t, d, p)
−0.1879
(−0.86)
0.3918
(1.83)
0.5776
(2.72)
0.98170.0364
Pine Tree Potentially Available AreaSpruce Tree Diameter
Equation (7)
(t)
−0.1371
(−1.43)
0.8767
(8.01)
1.2111
(11.03)
0.92650.4567Equation (7)
(t)
−0.7202
(−5.74)
1.1156
(9.56)
1.1587
(9.20)
0.88370.0047
Equation (13)
(t, d)
−0.1622
(−1.62)
0.8050
(7.43)
1.1211
(10.21)
0.93650.3423Equation (13)
(t, p)
−0.9502
(−7.81)
1.3568
(11.19)
1.2785
(10.15)
0.84150.0010
Equation (13)
(t, h)
−0.1294
(−1.36)
0.8629
(7.91)
1.1905
(10.85)
0.92900.4745Equation (13)
(t, h)
−0.0306
(−0.14)
0.8023
(6.40)
1.0270
(8.15)
0.93400.8713
Equation (13)
(t, d, h)
−0.1906
(−35.72)
0.7973
(7.32)
1.1201
(10.20)
0.93610.2650Equation (13)
(t, p, h)
−0.1389
(−1.02)
0.8442(6.61)1.0505
(8.34)
0.92990.5145
Spruce Tree HeightSpruce Tree Potentially Available Area
Equation (7)
(t)
−0.6456
(−5.18)
1.0927
(8.50)
1.4656
(11.56)
0.78570.0371Equation (7)
(t)
0.1726
(−1.81)
1.3916
(13.67)
1.9892
(19.75)
0.87890.6689
Equation (13)
(t, d)
−0.3298
(−3.01)
0.7273
(5.79)
0.9914
(7.82)
0.90890.1086Equation (13)
(t, d)
0.1667
(−2.07)
1.3726
(13.52)
1.9766
(19.63)
0.88030.6767
Equation (13)
(t, p)
−0.8052
(−9.57)
1.1729
(9.19)
1.5147
(11.94)
0.75410.0135Equation (13)
(t, h)
0.1931
(−1.77)
1.4267
(14.02)
2.0112
(19.98)
0.87580.6353
Equation (13)
(t, d, p)
−0.2877
(−2.66)
0.7246
(5.73)
0.9878
(7.79)
0.91160.1577Equation (13)
(t, d, h)
0.1511
(−2.24)
1.3364
(13.15)
1.9744
(19.61)
0.88070.7051
Birch Tree DiameterBirch Tree Height
Equation (7)
(t)
0.3437
(1.72)
1.1631
(6.86)
1.4458
(8.03)
0.89580.2556Equation (7)
(t)
0.1826
(0.61)
0.8385
(4.53)
0.9437
(4.90)
0.94120.3526
Equation (13)
(t, p)
0.0591
(0.03)
1.2562
(7.27)
1.4959
(8.31)
0.89430.8481Equation (13)
(t, d)
−0.5721
(−3.38)
0.7000
(3.94)
0.6981
(3.63)
0.94760.0004
Equation (13)
(t, h)
0.5554
(3.24)
0.9046
(5.14)
1.0379
(5.77)
0.93440.0151Equation (13)
(t, p)
0.0041
(−0.42)
0.8570
(4.67)
0.9856
(5.12)
0.93820.9836
Equation (13)
(t, p, h)
0.4326
(2.54)
0.8625
(4.87)
1.0415
(5.78)
0.94000.0530Equation (13)
(t, d, p)
−0.5783
(−3.37)
0.6992
(3.94)
0.6980
(3.63)
0.94770.0001
Table 6. Results of statistical measures and p-values of Student’s t-test for predictions of the number of trees per hectare using the mixed-effect mode.
Table 6. Results of statistical measures and p-values of Student’s t-test for predictions of the number of trees per hectare using the mixed-effect mode.
Curve
(Variables)
B
(%)
AB
(%)
RMSE
(%)
R2T
p-Value
Curve
(Variables)
B
(%)
AB
(%)
RMSE
(%)
R2T
p-Value
Number of All Trees per haNumber of Pine Trees per ha
Equation (25)
(t)
25.21
(2.48)
83.19
(7.78)
103.01
(9.09)
0.94230.1116Equation (25)
(t)
25.40
(4.97)
65.95
(9.80)
84.40
(11.80)
0.94680.0521
Equation (26)
(t, d)
37.91
(3.42)
80.13
(7.48)
94.71
(8.36)
0.94660.0110Equation (26)
(t, d)
30.50
(5.50)
66.09
(9.81)
80.97
(11.32)
0.94870.0163
Equation (27)
(t, h)
29.56
(2.75)
81.77
(7.65)
99.30
(8.76)
0.94500.0545Equation (27)
(t, h)
26.66
(5.07)
66.04
(9.82)
83.15
(11.62)
0.94770.0391
Equation (28)
(t, d, h)
44.49
(3.99)
81.81
(7.59)
94.53
(8.34)
0. 94400.0032Equation (28)
(t, d, h)
32.89
(5.88)
67.03
(9.94)
81.76
(11.43)
0.94680.0106
Number of Spruce Trees per haNumber of Birch Trees per ha
Equation (25)
(t)
−25.09
(−7.19)
69.48
(14.45)
108.44
(19.03)
0.93670.2581Equation (25)
(t)
1.50
(−5.47)
17.22
(18.27)
22.27
(15.80)
0.96840.7437
Equation (26)
(t, d)
−22.97
(−6.69)
61.56
(13.31)
94.47
(16.58)
0.95170.2353Equation (26)
(t, d)
4.46
(−3.58)
18.06
(17.94)
23.59
(16.74)
0.96340.3635
Equation (27)
(t, h)
−25.57
(−7.16)
65.93
(13.95)
99.81
(17.52)
0.94580.2118Equation (27)
(t, h)
2.17
(−4.89)
17.10
(18.0)
22.25
(15.79)
0.96830.6369
Equation (28)
(t, d, h)
−21.39
(−6.49)
62.38
(13.35)
95.09
(16.69)
0.95150.2714Equation (28)
(t, d, h)
4.47
(−3.57)
18.07
(17.94)
23.60
(16.75)
0.96340.3624
Table 7. Results of statistical measures and p-values of Student’s t-test for forecasts of the number of trees per hectare using the mixed-effect mode and Equation (25).
Table 7. Results of statistical measures and p-values of Student’s t-test for forecasts of the number of trees per hectare using the mixed-effect mode and Equation (25).
Curve
(Variables)
B
(%)
AB
(%)
RMSE
(%)
R2T
p-Value
Curve
(Variables)
B
(%)
AB
(%)
RMSE
(%)
R2T
p-Value
5-Year Forecast Period15-Year Forecast Period
Equation (25)
(t)
25.21
(2.48)
83.19
(7.78)
103.01
(9.09)
0.94230.1116Equation (25)
(t)
25.40
(4.97)
65.95
(9.80)
84.40
(11.80)
0.94680.0521
Equation (26)
(t, d)
37.91
(3.42)
80.13
(7.48)
94.71
(8.36)
0.94660.0110Equation (26)
(t, d)
30.50
(5.50)
66.09
(9.81)
80.97
(11.32)
0.94870.0163
Equation (27)
(t, h)
29.56
(2.75)
81.77
(7.65)
99.30
(8.76)
0.94500.0545Equation (27)
(t, h)
26.66
(5.07)
66.04
(9.82)
83.15
(11.62)
0.94770.0391
Equation (28)
(t, d, h)
44.49
(3.99)
81.81
(7.59)
94.53
(8.34)
0. 94400.0032Equation (28)
(t, d, h)
32.89
(5.88)
67.03
(9.94)
81.76
(11.43)
0.94680.0106
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rupšys, P.; Petrauskas, E. Modeling Number of Trees per Hectare Dynamics for Uneven-Aged, Mixed-Species Stands Using the Copula Approach. Forests 2023, 14, 12. https://doi.org/10.3390/f14010012

AMA Style

Rupšys P, Petrauskas E. Modeling Number of Trees per Hectare Dynamics for Uneven-Aged, Mixed-Species Stands Using the Copula Approach. Forests. 2023; 14(1):12. https://doi.org/10.3390/f14010012

Chicago/Turabian Style

Rupšys, Petras, and Edmundas Petrauskas. 2023. "Modeling Number of Trees per Hectare Dynamics for Uneven-Aged, Mixed-Species Stands Using the Copula Approach" Forests 14, no. 1: 12. https://doi.org/10.3390/f14010012

APA Style

Rupšys, P., & Petrauskas, E. (2023). Modeling Number of Trees per Hectare Dynamics for Uneven-Aged, Mixed-Species Stands Using the Copula Approach. Forests, 14(1), 12. https://doi.org/10.3390/f14010012

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