Allocation of Wood Density in European Oak ( Quercus robur L.) Trees Grown under a Canopy of Scots Pine

: The allocation of wood density in trees depends on many factors, but mainly on the tree species. A great number of studies have analysed wood density in dominant or codominant trees, but only a few have focused on trees grown under canopy. We examined the basic wood density and fresh wood density of natural origin oaks growing under canopy of artiﬁcially planted Scots pine. The major purpose of the work was to study the wood density allocation in different parts of the trees such as the trunk and branches. From a total of 80 oaks we selected eight model trees and measured biometric features of their trunks and crowns. Wood samples from different parts of the trunk and crown were collected after the trees were felled. We observed signiﬁcant differences between the average basic wood density (595 kg · m − 3 ) and the average fresh wood density (1031 kg · m − 3 ). The central part of the trunk and heartwood shows much higher density than the outer part of the trunk and sapwood, which corresponds to the model of ring-porous trees. Both types of wood density (basic and fresh) were also higher in the trunk than in the branches. The wood density of the branches differed between two zones: A1, which was closer to the trunk and had higher density; and A2, which was farther away from the trunk with lower density. Wood density shows positive correlation with crown length but not with crown width, which was more connected with diameter at breast height. We found lower value of slenderness than the value reported by other authors in oaks planted without canopy. The allocation of wood density in trees is associated with the potential mechanical load.

There are several types of wood density, which are defined differently in different scientific works. Mostly, wood density is defined as density of oven-dry wood density [6,[12][13][14], basic density or relative density [6,10], green density or fresh density (FD) [15], the density at 12% moisture content, according to various technical standards [16,17], or density calculated for unconventional levels of moisture content [18]. Some researchers have also used relative density based on various types of interpolation [19][20][21][22][23][24][25]. A completely novel approach is represented by X-ray densitometry [26][27][28]. Each type of wood density shows results on different scales, so comparing the results of various works is difficult and sometimes impossible. Predicting other wood properties based on various types of wood density is also complicated [29]. Many researchers have used another type of densityspecific gravity, which is used as an equivalent to wood density. Williamson [30] reported on the correct definitions of specific gravity. The work of Gao [24,25] is also helpful for understanding types of wood density. In this work, we used two types of wood density as described below.
Because of its ring-porous structure, the wood density of oak is often determined by methods based on X-ray densitometry and computed tomography [31][32][33][34][35][36][37], resistance drilling [3] and other methods that calculate the density relatively [38]. The results of such measurements allow comparison among different populations [5,39,40] and explanations of reactions in the trees [31]. From the point of view of forestry and industrial practices, they are less useful. Wood density calculated using technical standards and with detailed methodological description is more useful for the wood industry and forestry [14,[41][42][43][44].
Most of the research relating to oak wood density focuses on dominant or codominant trees [37,[45][46][47], with clone studies occasionally conducted on seed orchards [48]. Studies on trees growing in different biosocial positions are less common [3] and there have been few studies on trees growing under canopy [49]. The general focus on dominant trees comes from the analysis of tree rings [46], which are typically collected from dominant trees.
For many years after the Second World War, the Scots pine planting model dominated in Poland as a commercial basis for the wood industry. Pine was also planted in locations where other, primarily deciduous, species occurred. These monoculture stands were sensitive to many factors, such as gradations, fires and strong winds. For the last 30 years, foresters have rebuilt many coniferous forests into mixed or broadleaved forests. Our experimental research area represents a mixture of pine and oak, which was established in an unusual way. Scots pine (115 years old) was planted on sites that prefer oaks and other deciduous trees. The under canopy oak forest was not planted with all trees being of natural origin, growing under the process of natural succession.
The principal aim of this study was to record wood density allocation in different parts of oak trees. As an additional goal, we wanted to provide estimated correlations between wood density and the biometric features of the trees. Biometric characteristics in combination with wood density can help predict tree behaviour after pine removal. We hypothesised that: (H1) the radial variability in wood density would follow the typical model for ring-porous trees, with the wood density decreasing from pith to bark; (H2) we expected the wood density to depend on the height and lower wood density in the branches due to juvenile character of the wood; and (H3) we also assumed that the wood density would be related to the light-dependent biometric features of trees.

Study Area, Tree Selection and Wood Sampling Procedure
The study area was in the central part of Poland, at latitude − 52 • 09 N, 17 • 27 E and an altitude 78 m a.s. The research was carried out in a mixed forest with 115-year-old Scots pine and younger oak under the canopy. The Scots pine was planted as a monoculture, and all the trees are of the same age. The oaks, however, occur as a naturally established forest, and are aged between 60 and 80 years. The overstory basal area was 28.9 m 2 /ha. Site index in the Polish classification was III on a six-grade scale, meaning average conditions. The forest is treated as water-protective, and the last thinning operations were carried out more than 20 years ago. A sample plot with 80 oaks was selected. The diameter of breast height (DBH) and height (H) of all the trees were measured. Measured traits were used to calculate slenderness as the ratio between H and DBH. Based on distribution and height curve, we selected eight as representative of the sample plot. The selected trees had regular healthy crowns. We measured the following biometric features of these model trees: DBH, H, length of living crown, width of crown and the direction of north (Table 1). The selected trees were felled and the wood samples were collected according to the following methods. Four 2 cm thick discs were taken from each tree from: (1)-base of the trunk (at 0.2 m); (2) the level of the DBH (1.3 m); (3)-base of the crown (relative height); (4) the middle of the crown length; and (5) six discs from branches based on the scheme described in Figure 1. The first disc (from zone A1) was taken 10 cm from the trunk, the second (from zone A2) from the middle of the crown width. All the samples were collected in July 2019.  The selected trees were felled and the wood samples were collected according to the following methods. Four 2 cm thick discs were taken from each tree from: (1)-base of the trunk (at 0.2 m); (2) the level of the DBH (1.3 m); (3)-base of the crown (relative height); (4) the middle of the crown length; and (5) six discs from branches based on the scheme described in Figure 1. The first disc (from zone A1) was taken 10 cm from the trunk, the second (from zone A2) from the middle of the crown width. All the samples were collected in July 2019.

Measurement of Wood Density
The mass of fresh (green) wood was determined on site for complete discs (from the trunk and branches) after removing the bark. The dry wood mass was measured in two ways: from discs (but only from the branches, because they were too small to split into pieces) and from small samples after cutting (only from the trunk). Small cubed samples were cut from the discs in radial rows, oriented in the east-west and north-south directions ( Figure 1). Not all the diameters of the samples were measured because the volume calculation was performed using a water displacement method, but they ranged between 1.4 cm and 1.6 cm in every direction. The number of samples from each radius depended on the diameter of the disk and ranged between four and eight ( Figure 1). We collected 645 small samples in total. The samples were then dried in a laboratory dryer at 105 °C until they reached constant weight. The mass of each sample was measured by using two laboratory scales to provide accuracy of 0.1 g (for discs) and 0.01 g (for the small samples). Scales were also used to measure the volume of wood above the fibre saturation point by using the water displacement method [50]. For the volume measurement, wood samples

Measurement of Wood Density
The mass of fresh (green) wood was determined on site for complete discs (from the trunk and branches) after removing the bark. The dry wood mass was measured in two ways: from discs (but only from the branches, because they were too small to split into pieces) and from small samples after cutting (only from the trunk). Small cubed samples were cut from the discs in radial rows, oriented in the east-west and north-south directions ( Figure 1). Not all the diameters of the samples were measured because the volume calculation was performed using a water displacement method, but they ranged between 1.4 cm and 1.6 cm in every direction. The number of samples from each radius depended on the diameter of the disk and ranged between four and eight ( Figure 1). We collected 645 small samples in total. The samples were then dried in a laboratory dryer at 105 • C until they reached constant weight. The mass of each sample was measured by using two laboratory scales to provide accuracy of 0.1 g (for discs) and 0.01 g (for the small samples). Scales were also used to measure the volume of wood above the fibre saturation point by using the water displacement method [50]. For the volume measurement, wood samples were immersed in water for 72 h. BD and FD were measured using the following formulae: where: BD = basic density (kg·m −3 ), M 0% = mass of sample at 0% of moisture content and V max = volume of sample at maximum swelling.
where: FD = fresh density (kg·m −3 ), M = mass of sample as fresh wood (the moisture content refers to wood removed from just a fallen tree) and V max = volume of sample as fresh wood (as above).

Statistical Analysis
Statistical analyses were applied in two ways-to compare different groups of wood density and to correlate wood density with the biometric parameters. The basic statistics were plotted as graphs using a standard boxplot (frame = Q25 to Q75, centre line = median). We also used a violin chart to show the density of the values. In order to verify the statistical differences between the groups, we used nonparametric tests because of lack of a normal distribution in the wood density. The Mann-Whitney test was used for two groups, and the Kruskal-Wallis test for more than two groups. All the tests were performed at a level α = 0.05 using Statistica software. In order to estimate the relationships between the wood density and biometric features, we built a correlation matrix based on Spearman's rank coefficient.

Wood Density of the Trunk
The mean BD value of the trunks was 595 kg·m −3 with a wide variation in range. The maximum value was 880 kg·m −3 while the minimum was only 318 kg·m −3 . The mean FD was much higher, amounting to 1031 kg·m −3 . The lower quartile (Q25) for the FD also had a high value (983 kg·m −3 ; Table 2). The difference between BD and FD oscillated around 40%. The average values for the branches were lower for both the BD at 573 kg·m −3 and the FD at 1007 kg·m −3 ( Table 2). All of the differences between the trunk and branches were significant at a level of α = 0.05. The radial allocation of the BD presented as typical distribution for ring-porous wood, where the wood density was higher close to the pith, decreasing towards the bark. In the central part of the stem, most samples had values above 600 kg·m −3 . The median in the subsequent zones of trunk decreased in line with the cambial age, down to 500 kg·m −3 close to the bark (Figure 2a). A similar pattern was apparent when the trunk was divided into heartwood and sapwood, with the mean value of the heartwood density being 624 kg·m −3 , much higher than the sapwood density (538 kg·m −3 ). Sapwood represented a very narrow part of the crosssection but dominated in the uppermost discs. The broad width of the heartwood also resulted in higher values in most samples than in the sapwood, and not only the average value. The violin charts and boxplots showed that most samples were focused on the median in the heartwood rather than in the sapwood, which was more variable (Figure 2b). A comparison of wood densities in relation to compass directions (azimuth) showed no differences.
In the axial direction, the highest wood density occurred at the base of the trunk, where mean values equalled 632 kg·m −3 and, in most samples, exceeded the threshold of 600 kg·m −3 . The lowest values were observed at 2 m high, the wood density increasing again at the top of the tree. However, the differences between the sample heights (except the basal level) were not significant ( Figure 3). The boxes containing the values between Q25 and Q75 have a narrow range, with a tendency to expand at the top of the trunk. The rest of the sample population has a wide range. A similar pattern was apparent when the trunk was divided into heartwood and sapwood, with the mean value of the heartwood density being 624 kg·m −3 , much higher than the sapwood density (538 kg·m −3 ). Sapwood represented a very narrow part of the cross-section but dominated in the uppermost discs. The broad width of the heartwood also resulted in higher values in most samples than in the sapwood, and not only the average value. The violin charts and boxplots showed that most samples were focused on the median in the heartwood rather than in the sapwood, which was more variable (Figure 2b). A comparison of wood densities in relation to compass directions (azimuth) showed no differences.
In the axial direction, the highest wood density occurred at the base of the trunk, where mean values equalled 632 kg·m −3 and, in most samples, exceeded the threshold of 600 kg·m −3 . The lowest values were observed at 2 m high, the wood density increasing again at the top of the tree. However, the differences between the sample heights (except the basal level) were not significant ( Figure 3). The boxes containing the values between Q25 and Q75 have a narrow range, with a tendency to expand at the top of the trunk. The rest of the sample population has a wide range. 12, x FOR PEER REVIEW 6 of 12 The first hypothesis (H1) was confirmed by the results, supporting the model described in most publications. The decrease in wood density from pith to bark in ring-porous species was described a long time ago [41,51]; however, current research has suggested that this relationship is more complicated and depends on various factors. Bahmani et al. [44] reported a decrease in wood density for Persian oak, but only for small diameters (10-30 cm). For larger diameters, the relationship is not so clear. A decrease in density has been observed in European oak and is more obvious in the upper part of the trunk [31]. A greater difference between the central and outer parts of the trunk has been noted previously in classical silvicultural regimes more than in dynamic silvicultural regimes. Decreasing wood density with cambial age has also been observed in cork oak [35], the authors proving that wood density decreases from pith to bark most rapidly up to the 15th ring. The opposite correlation was found, surprisingly, during the testing of ovendry wood density in Quercus. robur and Q. rubra [14]. The researchers also reported high variability in the wood density in the outermost part of the trunk. In a long-term series of annual rings taken from Q. petraea in France, a reduction in wood density was observed throughout the period of the study. In periods of both 50 and 100 years, the wood density decreased, while stability occurred only in small time intervals [33]. Decreasing wood density may be the result of the narrowing of annual rings. This has been a source of discussion among some scientists [47,51] who observed that narrowing rings correspond with cambial age. Lei et al. [51] also found an increase in vessel diameter with cambial age. The difficulty in determining how basic features affect wood density arises from the complicated relationships between wood elements such as vessels, fibres, wood rays, cell walls and various forms of parenchyma [52]. The correlation between the number of annual rings and wood density has not been clearly explained. In our tests, we found a very weak relationship among the tested wood samples (Figure 4), with correlation being negative, but not significant, and the slope of the lines being the result of the presence of a few very narrow rings in two of the trees. The first hypothesis (H1) was confirmed by the results, supporting the model described in most publications. The decrease in wood density from pith to bark in ring-porous species was described a long time ago [41,51]; however, current research has suggested that this relationship is more complicated and depends on various factors. Bahmani et al. [44] reported a decrease in wood density for Persian oak, but only for small diameters (10-30 cm). For larger diameters, the relationship is not so clear. A decrease in density has been observed in European oak and is more obvious in the upper part of the trunk [31]. A greater difference between the central and outer parts of the trunk has been noted previously in classical silvicultural regimes more than in dynamic silvicultural regimes. Decreasing wood density with cambial age has also been observed in cork oak [35], the authors proving that wood density decreases from pith to bark most rapidly up to the 15th ring. The opposite correlation was found, surprisingly, during the testing of oven-dry wood density in Quercus. robur and Q. rubra [14]. The researchers also reported high variability in the wood density in the outermost part of the trunk. In a long-term series of annual rings taken from Q. petraea in France, a reduction in wood density was observed throughout the period of the study. In periods of both 50 and 100 years, the wood density decreased, while stability occurred only in small time intervals [33]. Decreasing wood density may be the result of the narrowing of annual rings. This has been a source of discussion among some scientists [47,51] who observed that narrowing rings correspond with cambial age. Lei et al. [51] also found an increase in vessel diameter with cambial age. The difficulty in determining how basic features affect wood density arises from the complicated relationships between wood elements such as vessels, fibres, wood rays, cell walls and various forms of parenchyma [52]. The correlation between the number of annual rings and wood density has not been clearly explained. In our tests, we found a very weak relationship among the tested wood samples (Figure 4), with correlation being negative, but not significant, and the slope of the lines being the result of the presence of a few very narrow rings in two of the trees. Other authors have been sceptical about there being a clear connection between these two features and between the proportion of latewood and ring width. They have described these relationships as low intra-ring variability with slight differences between the densities of earlywood and latewood--known as the low ring heterogeneity index [35]. The correlation between ring width and mean density has also been recorded as being weak or insignificant in sessile oak in northern France, with the authors suggesting that early wood width varies between sample plots under the influence of various site conditions [33]. Some researchers have reported a unique point of view that involves species dependences. A very weak correlation between proportion of latewood and ring width has been observed in Q. rubra, but this occurs much more strongly in Q. robur [14]. The negative correlation between ring width and latewood density has been described by Mosedale et al. [48] for a clone of Q. robur and Q. petraea in Germany in a seed orchard.
A decrease in wood density with cambial age lowers the density of sapwood, as can be seen from our results (Figure 2b). A similar difference between heartwood and sapwood has been reported by Kiaei [42], in relation to oven-dry wood density in Q. castaneaefolia.
The second hypothesis (H2) was only partially confirmed in relation to the lower part of the trunk having higher wood density (Figure 3). We did not observe a decrease in wood density moving up the trunk; the data were variable. Probably the only part of the trunk that bears the most load reacts by increasing its density. King et al. [53], in studies on over 20 tropical species, indicated a strong correspondence between wood density and stem biomass, and a poor relationship between wood density and stem volume, which may suggest that the parts of the trunk carrying the most load may be denser.

Wood Density of the Branches
The wood density in the branches differed significantly in both the studied zones (A1 and A2). In zone A1 closest to the trunk, the branches had higher densities (both BD and FD) than those farther from the trunk (zone A2) ( Figure 5). The differences were not great, but were distinct and significant. The average wood density for zone A1 equalled 595 (BD) and 1040 (FD) kg·m −3 , and for zone A2, respectively, 551 (BD) and 974 (FD) kg·m −3 . The Other authors have been sceptical about there being a clear connection between these two features and between the proportion of latewood and ring width. They have described these relationships as low intra-ring variability with slight differences between the densities of earlywood and latewood--known as the low ring heterogeneity index [35]. The correlation between ring width and mean density has also been recorded as being weak or insignificant in sessile oak in northern France, with the authors suggesting that early wood width varies between sample plots under the influence of various site conditions [33]. Some researchers have reported a unique point of view that involves species dependences. A very weak correlation between proportion of latewood and ring width has been observed in Q. rubra, but this occurs much more strongly in Q. robur [14]. The negative correlation between ring width and latewood density has been described by Mosedale et al. [48] for a clone of Q. robur and Q. petraea in Germany in a seed orchard.
A decrease in wood density with cambial age lowers the density of sapwood, as can be seen from our results (Figure 2b). A similar difference between heartwood and sapwood has been reported by Kiaei [42], in relation to oven-dry wood density in Q. castaneaefolia.
The second hypothesis (H2) was only partially confirmed in relation to the lower part of the trunk having higher wood density (Figure 3). We did not observe a decrease in wood density moving up the trunk; the data were variable. Probably the only part of the trunk that bears the most load reacts by increasing its density. King et al. [53], in studies on over 20 tropical species, indicated a strong correspondence between wood density and stem biomass, and a poor relationship between wood density and stem volume, which may suggest that the parts of the trunk carrying the most load may be denser.

Wood Density of the Branches
The wood density in the branches differed significantly in both the studied zones (A1 and A2). In zone A1 closest to the trunk, the branches had higher densities (both BD and FD) than those farther from the trunk (zone A2) ( Figure 5). The differences were not great, but were distinct and significant. The average wood density for zone A1 equalled 595 (BD) and 1040 (FD) kg·m −3 , and for zone A2, respectively, 551 (BD) and 974 (FD) kg·m −3 .
The violin displays had more focused values for BD and more scattered values for FD ( Figure 5). violin displays had more focused values for BD and more scattered values for FD ( Figure  5).

Figure 5.
Comparison between the basic density (BD) and fresh density (FD) of branches at different distances from the trunk (zones A1 and A2).
Wood density was much lower in the branches, and decreased with distance from the trunk (Figure 5), which may result from a biomechanical relationship in the tree. Yang et al. [54] studied competition in Liaodong oak, describing the reaction of biomass allocation to changes in the tree architecture depending on specific conditions. They found that the biomass was reduced in the branches rather than the stem when the diameter and length of the branches were reduced, but not their number.

Biometric Tree Relationships
A matrix heat map based on Spearman ranks highlighted a few interesting aspects, especially concerning the relationships between the crown parameters and wood density ( Figure 6). Both densities, BD and FD, were related to length of crown (0.7 and 0.8), but only a weak correlation was observed between wood density and width of crown. Crown length was highly correlated with slenderness (0.8), which could be a key feature for explaining the growth under the canopy. Slenderness also showed a moderate correlation with wood density. The second element worth noting was the high correlation between width of crown and area of sapwood (0.8). Wood density was much lower in the branches, and decreased with distance from the trunk (Figure 5), which may result from a biomechanical relationship in the tree. Yang et al. [54] studied competition in Liaodong oak, describing the reaction of biomass allocation to changes in the tree architecture depending on specific conditions. They found that the biomass was reduced in the branches rather than the stem when the diameter and length of the branches were reduced, but not their number.

Biometric Tree Relationships
A matrix heat map based on Spearman ranks highlighted a few interesting aspects, especially concerning the relationships between the crown parameters and wood density ( Figure 6). Both densities, BD and FD, were related to length of crown (0.7 and 0.8), but only a weak correlation was observed between wood density and width of crown. Crown length was highly correlated with slenderness (0.8), which could be a key feature for explaining the growth under the canopy. Slenderness also showed a moderate correlation with wood density. The second element worth noting was the high correlation between width of crown and area of sapwood (0.8).
Heartwood area was highly related to DBH. The two parameters that showed the most relationships in the matrix table-DBH and wood density-did not correlate with each other. Likewise, the tree H did not correlate with wood density. Heartwood area was highly related to DBH. The two parameters that showed the most relationships in the matrix table-DBH and wood density-did not correlate with each other. Likewise, the tree H did not correlate with wood density.
The correlation heat map resulting from the biometric features of the trees also provided some pointers to potential future work. There was a high correlation between heartwood area and DBH ( Figure 6). A comparable correlation between heartwood width and stem diameter at DBH was described by Sousa et al. [47], who studied dependencies in Q. faginea. Heartwood width also depends on site conditions, as the authors emphasised [47]. Our results show a lack of correlation between wood density and DBH, which supports the findings of other authors, who have pointed out that the relationship between DBH and wood density in Q. robur is not clear [55]. Different results were provided by Kranjc et al. [3], who observed that pubescent oak diameter has a positive effect on wood density. In a specific study area, the competition between trees, potential gaps in the main pine forest and other relationships that exist in the living forest environment must be taken into account. Additionally, we calculated the slenderness of 80 trees (Figure 7) in order to compare that with a study carried out in the same region by Kaźmierczak et al. [56] (Table  3), who showed the variability in slenderness in oak growth in a basal overstory with full access to light. We obtained an average value of 0.81 for all the trees. All the oaks were grown in a mixed-aged forest, so calculations for 10-year intervals were not possible. Therefore, we present a histogram for the total area. Most trees had values below 1.0, with the slenderness being lower than reported by the abovementioned authors. It is likely that the shadow of the pine trees affected oak growth, causing lower values. The correlation heat map resulting from the biometric features of the trees also provided some pointers to potential future work. There was a high correlation between heartwood area and DBH ( Figure 6). A comparable correlation between heartwood width and stem diameter at DBH was described by Sousa et al. [47], who studied dependencies in Q. faginea. Heartwood width also depends on site conditions, as the authors emphasised [47]. Our results show a lack of correlation between wood density and DBH, which supports the findings of other authors, who have pointed out that the relationship between DBH and wood density in Q. robur is not clear [55]. Different results were provided by Kranjc et al. [3], who observed that pubescent oak diameter has a positive effect on wood density. In a specific study area, the competition between trees, potential gaps in the main pine forest and other relationships that exist in the living forest environment must be taken into account. Additionally, we calculated the slenderness of 80 trees (Figure 7) in order to compare that with a study carried out in the same region by Kaźmierczak et al. [56] ( Table 3), who showed the variability in slenderness in oak growth in a basal overstory with full access to light. We obtained an average value of 0.81 for all the trees. All the oaks were grown in a mixed-aged forest, so calculations for 10-year intervals were not possible. Therefore, we present a histogram for the total area. Most trees had values below 1.0, with the slenderness being lower than reported by the abovementioned authors. It is likely that the shadow of the pine trees affected oak growth, causing lower values.

Conclusions
The basic wood density accounted for only 60% of the fresh wood density, so these two parameters represent two different scales. A visible decrease in density in the radial direction was observed. The central part of the trunk and heartwood had much higher density than the outer part and the sapwood. Both types of wood density (BD and FD) were also higher in the trunk than in the branches. The wood density of the branches differed between the two zones (i.e., closer to versus farther from the trunk), which could signal a mechanical relationship in the tree. The other biometric features of the tree elements correlated with different strength, with wood density showing a positive correlation with crown length, but not with crown width, which corresponds better with DBH. We found lower slenderness than the values reported by other authors in oaks planted without canopy. On one hand, low slenderness predicts better resistance to wind after removing the canopy, but on the other hand, the oaks had grown for many years in uncomfortable conditions, which may have caused problems, for example, related to health. These observations may be used to direct future research in a wider field.

Conclusions
The basic wood density accounted for only 60% of the fresh wood density, so these two parameters represent two different scales. A visible decrease in density in the radial direction was observed. The central part of the trunk and heartwood had much higher density than the outer part and the sapwood. Both types of wood density (BD and FD) were also higher in the trunk than in the branches. The wood density of the branches differed between the two zones (i.e., closer to versus farther from the trunk), which could signal a mechanical relationship in the tree. The other biometric features of the tree elements correlated with different strength, with wood density showing a positive correlation with crown length, but not with crown width, which corresponds better with DBH. We found lower slenderness than the values reported by other authors in oaks planted without canopy. On one hand, low slenderness predicts better resistance to wind after removing the canopy, but on the other hand, the oaks had grown for many years in uncomfortable conditions, which may have caused problems, for example, related to health. These observations may be used to direct future research in a wider field.

Data Availability Statement:
The study did not report any data.