Abstract:The evolution of the 4-variate probability distribution of the diameter at the breast height, total height, crown base height, and crown width against the age in a forest stand is of great interest to forest management and the evaluation of forest resources. This paper focuses on the Vasicek type 4-variate fixed effect stochastic differential equation (SDE) to quantify the dynamic of tree size components distribution against the age. The new derived 4-variate probability density function and its marginal univariate, bivariate, trivariate, and conditional univariate distributions are applied for the modeling of stand attributes such as the mean diameter, height, crown base height, crown width, volume, and slenderness. All parameters were estimated by the maximum likelihood procedure using a dataset of 1630 Scots pine trees (12 stands). The results were validated using a dataset of 699 Scots pine trees (five stands). A newly developed 4-variate simultaneous system of SDEs incorporated covariance structure driving changes in tree size components and improved predictions in one tree size component given the other tree size components in the system.