Next Article in Journal
Mean Leaf Angles Affect Irrigation Efficiency and Physiological Responses of Tropical Species Seedling
Previous Article in Journal
Fine Root Production and Soil Available Nutrients in Rubber Monoculture versus Rubber–Flemingia macrophylla Agroforestry
Previous Article in Special Issue
Impact of Ice-Storms and Subsequent Salvage Logging on the Productivity of Cunninghamia lanceolata (Chinese Fir) Forests
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Influence of the Provenance and Spatial Structure on the Growth of European Silver Fir (Abies alba Mill.) of Autochthonous Origin in a Forest Plantation in the Białowieża Forest

1
Institute of Forest Sciences, Bialystok University of Technology, Wiejska 45E Street, 15-351 Bialystok, Poland
2
Faculty of Forestry, Belarusian State Technological University, Sverdlova Street 13a, 220006 Minsk, Belarus
3
Institute of Biological Sciences, Cardinal Stefan Wyszynski University in Warsaw, Woycickiego 1/3 Street, 01-938 Warsaw, Poland
4
Department of Biochemistry and Pharmacogenomics, Faculty of Pharmacy, Medical University of Warsaw, 02-097 Warsaw, Poland
5
Rocky Mountain Research Station, Forest Service, U.S. Department of Agriculture, 2500 S. Pine Knoll Dr., Flagstaff, AZ 86001, USA
6
Environmental Sciences, University of Guelph, Guelph, ON N1G 2W1, Canada
7
Forest Protection Department, Forest Research Institute, Braci Lesnej 3, 05-090 Sekocin Stary, Poland
*
Authors to whom correspondence should be addressed.
Submission received: 9 May 2022 / Accepted: 22 May 2022 / Published: 26 May 2022

Abstract

:
Currently, a change in the species composition of the Białowieża Forest, eastern Poland, is occurring. Parallel to the dying of spruce (Picea abies L. Karst.), caused by Ips typographus (L.), there is a spread of deciduous tree species, among which hornbeam (Carpinus betulus L.) is the most active. Against the background of climate change, it is of interest to consider the possibility of reintroducing relict species in areas occupied until recently by spruce. One of these is silver fir (Abies alba Mill.), which is biologically and ecologically very similar to spruce. The Tisovik tract (the Belarusian part of the Białowieża Forest) is the most northeastern and the only preserved refuge of autochthonous silver fir in the region. Therefore, it is the most suitable source of propagation of this species outside the mountainous part of its range. The target area of our study was a 26-year old artificial stand included in the information system of forest genetic resources conserved in Europe. It was created in 1996 in the Polish part of the Białowieża Forest in the Hajnówka Forest District (Forestry Wilczy Jar) from 10 half-sib families originating from seeds collected in the Tisovik tract. The goal of the study was to assess the influence of provenance factors and spatial structure on productivity and phenotypic variation as measured by diameter (DBH). The study’s initial results showed that the provenance factor’s influence at the pole wood stage manifests itself more clearly than spatial structure. Even within the framework of individual half-sibs, there was a clear differentiation of diameters between investigated trees. Interpretation of the obtained results allowed us to conclude that there is a “threshold value” of the distance between trees in a stand of a given age, which determines the point when its increase no longer contributes to an increase in DBH. The most promising phenotypes were selected for further research and practical actions to reproduce silver fir further.

1. Introduction

Since 2012, there has been extensive mortality of the Norway spruce (Picea abies L. Karst.) due to insect attacks [1,2]. The most harmful pest is the European spruce or eight-toothed bark beetle (Ips typographus L.)—a secondary pest that primarily affects weakened and depressed trees. A change in the hydrological regime due to global warming and acid rain facilitated the weakening of spruce. It should be borne in mind that the spruce in the Białowieża Forest in Poland (which means further Forest) is located in the south of the northeastern part of the range [1]. Therefore, it is especially susceptible to changes in hydrology. The bark beetle actively attacks even healthy and robust trees during mass reproduction. This negative phenomenon occurs in the Forest, in Poland, including areas with rich soil conditions that fully meet the ecological and biological requirements of A. alba [3,4,5]. This species is natural to the Forest but has survived only in a small refugium, in the Tisovik tract in the Belarusian part of the Forest [6,7]. It lies 500 km northeast of the edge of the continuous range of silver fir in Poland and 120 km from the nearest natural island habitat in the Jata Reserve [8]. At present, it is the rarest coniferous species in the Forest.
Taking into account the autochthonous nature of the silver fir from Tisovik, it is reasonable to assume that the A. alba does not pose a threat to biodiversity, but, on the contrary, enhances it. This is unlike, for example, a very aggressive invasive species, Quercus rubra (L.), which displaces Q. robur (L.) and thus poses a real threat to the stability of forest ecosystems [9,10].
It is interesting to study the possibility of including silver fir as a candidate for artificial reforestation in the Forest and adjacent forest districts. According to our previous research, silver fir has exhibited high growth rates and vigor in stands of different ages and origins [11,12]. These performance measures justify establishing a seed base of silver fir for targeted afforestation of those areas freed up by spruce mortality.
A half-sib family is a collection of plants obtained from seeds from free pollination from one tree. In this context, the study of the unique gene pool of autochthonous silver fir, presented in the form of half-sib families (what we understand as “provenance”) in the Hajnówka Forest District (FD), is necessary to identify the most promising genotypes. These promising genotypes could be preferred seed sources, especially for establishing seed plantations. Therefore, the studies of the first stage [11,12], continued in this research, were conducted using traditional forest mensuration and phenotypic analyses.
It is essential to note that despite its small size—approximately 20 trees at present [7,13,14]—the autochthonous population of silver fir in the Tisovik tract is similar in the degree of polymorphism and heterozygosity (based on isoenzymatic and DNA analyses) compared to the populations from the main distribution area in the Carpathians [15]. The list of publications on the genetics of fir in the Carpathians, without comparison with the population from Tisovik, is quite extensive [16,17,18]. Unfortunately, only this publication provides a genetic comparison of the population from Tisovik with populations from the main area. This small population shows no inbreeding effect [15], suggesting the value of a detailed genetic study of half-sib families, especially in selecting candidates for preferred seed trees. A study by Pawlaczyk et al. [19] showed an excess of heterozygotes on a half-sib plantation in the Hajnówka FD, which is the offspring of such a small population as the fir group in the Tisovik tract. The authors attribute this to a “bottleneck effect” and assume that pollination of the remnant firs in Tisovik could have been carried out, among other things, from two stands (artificial, of unknown geographical origin) in the Polish part of the Białowieża Forest. Based on genetic analysis, Mejnartowicz [20] became convinced that these two stands are not descendants of autochthonous fir from the Tisovik tract (even though the nearest stand is located 8.5 km in a straight line from this refuge). If the assumption of Pawlaczyk et al. [19] is correct, then not only the progeny of local cross-pollinated silver fir from Tisovik is represented in the half-sib plantation, but also individuals that arose as a result of gene flow together with pollen from these two stands of unknown origin, which are located in two forest subcompartments, 453Aa and 498 Ci of the Białowieża FD. We were unsuccessful in finding documents on the history of their appearance (between 1920 and 1939) in Polish archives, likely due to the destruction occurring during the Second World War. They might exist in some archives in Belarus, which are not available to us at the moment.
It would be logical to assume that there is a hereditary high degree of differentiation of DBH (diameter at breast height) between the families of the studied silver fir half-sibs and within families. Confirming or refuting this assumption was one of our research objectives, presented below. Note that at this stage, the genetic factor is conventionally represented as the membership of trees in one of the families of half-sibs (provenance). In addition, only the next step will be a study based on the evaluation of DNA polymorphism of the identified best specimens for this phenotype (DBH) trait.
Competitive interactions between trees occur at all stages of the forest stand dynamics, but they are the most pronounced at the second—stem—exclusion stage of forest formation [21]. The density and spacing of the trees determine the competitive influence experienced by any single tree. These interactions, in turn, influence the spatial structure of the stand, affecting both productivity and species diversity of forest communities. Optimising the spatial structure of forest stands as a function of age classes, and taking into account the biological characteristics of the breed is one of the most effective methods of increasing a stand’s productivity and stability. Thus, regulating the spatial distribution of trees is the most important tool for managing artificial forest ecosystems [22].
The development of stand growth modeling raises the question of how to simulate the initial states of these models and, in particular, how to simulate realistic spatial structures [23]. The parameters of the spatial structure of forest stands determine and characterize them quantitatively. However, only understanding how to change the spatial structure effectively is crucial for forest management [24]. Knowledge of the spatial structure of stands can be used to develop silvicultural regimes to increase the area of semi-natural forests established as even-aged plantations [25]. Growth models of stands considering their spatial characteristics can predict productivity and evaluate changes in forest fund parameters and are recommended for use in forestry [26].
Thus, the aim of this study was a comparative evaluation of the influence of the provenance factor and spatial structure on the growth of trees in a 22-year-old artificial stand of autochthonous silver fir from the Tisovik tract.
The objectives of this study were as follows:
  • For the provenance factor:
    • to evaluate the variability of the most important silvicultural parameter—DBH—and thus reveal the degree of possible differentiation both between the families of the studied silver fir half-sibs and within families;
    • to identify the most productive phenotypes, which are potential sources of vegetative and seed material for effective reforestation of the autochthonous population of silver fir and objects of further genetic studies;
  • For the spatial structure factor:
    • to determine the parameters of the spatial structure effectively affecting the DBH,
    • to determine the approaches of effective regulation of the spatial structure of the stand in the near term.

2. Materials and Methods

2.1. Object of Research

The bulk of the study was conducted in 2018 in a half-sib plantation represented by ten families of silver fir from the Tisovik tract, located in the Belarusian part of the Białowieża Forest. The plantation was established by Professor A. Korczyk in Hajnówka FD (subcompartment 416Ag (now 416Ad)) in rich soil conditions (mixed fresh forest) in 1996 with 4-year-old seedlings grown from seeds [12]. The half-sib number corresponds to the number of the tree in the refugium. Unfortunately, according to our information, some mother trees have died for various reasons. However, they are presented in the form of half-sib families at the research object. The arrangement of plants at planting in furrows created in a recent cutting is 1.3 × 1.0 m. A total of 1551 seedlings (7692 specimens/ha) were planted in 65 rows of 27 m length (north-south), numbered in ascending order from west to east. The half-sibs located in the first, oldest plot were examined. At the time of the study, it had been 22 years since the trees were planted, and their biological age was 26 years. The average diameter of the trees in the stand at the time of the study was 10.0 cm, and the average height was 9.4 m [12]. The other two plots, planted in 1998, are included in our future research plans under a scientific grant. All three plots are fenced to protect them from animals (Figure 1).
After conducting research in 2018, we did not let the experimental site out of sight. In 2020 and 2021, the collection of cones was carried out. Seeds isolated from them were sown in the nursery Budy of Hajnówka FD. The pictures (Figure 2 and Figure 3) taken from the lift used to collect the cones show that the trees are in excellent condition. In 2021, the diameter and height of the 9 tallest and thickest trees were also measured. In the same year, we made a preliminary comparative assessment of trees assigned for felling in 2015 but not harvested due to the change in the status of the object.
The crowns are pointed (positive growth trend) with dense, bright green needles. One can notice a pronounced difference in heights, which corresponds to a difference in diameters. However, the crowns of trees of the lower growth classes also look good. The aerial view shows no apparent openings in the canopy. Furthermore, 38 trees formed cones in 2021, and some showed repeated (after 2020) fruiting. As a rule, these are trees with the largest diameter. However, fruiting was also recorded on small-diameter trees. All these facts indicate that the condition of the trees on the site for three years after the study continues to be very good and is in no way a cause for concern.

2.2. Collection of Field Data and Its Standard Processing

Tree growth was assessed by analyzing the diameter at 1.3 m from the root collar or DBH, which represents it well and is easy to measure in dense stands (in contrast to height measurement). To determine the influence of the provenance factor for the ten half-sib families, we measured the diameters of all 701 trees to an accuracy of 0.1 cm. The representation of the families varied considerably: from 11 specimens in families 16 and 21 to 167 and 180 in 17 and 15, respectively. A relatively large number of seedlings were planted in sibs Nos. 2 and 5, respectively, 336 and 258 trees. A fairly close number of trees was initially represented by No. 3—75 pcs., No. 4—70 pcs., No. 11—66 pcs., and No. 12—38 pcs. [12].
At the time of the study, the survival of plants (as a percentage) for siblings relative to their initially planted number was No. 2—46.7, No. 3—41.3, No. 4—34.3, No. 5—28.3, No. 11—51.3, No. 12—34.2, No. 15—54.1, No. 16—61.1, No. 17—51.1, and No. 21—40.7. Despite a low survival rate, some sibs had started with such a large number of individuals that we still had many trees to evaluate. For example, sib No. 5 had the lowest survival rate—28.3. However, thanks to a high number of initial plantings—258, we have the opportunity to study 73 individuals. While sib No. 16 has the highest survival rate—61.1%, it is represented by only 11 individuals. A table with tree survival in different years of observation is presented in our previous article [12].
We wished to determine the influence of spatial (horizontal) structure on DBH in 126 randomly selected specimens (n). Accordingly, we measured the distances in four directions (north, south, west, and east) to nearest neighbors with a tape, and their average distance was calculated [12].
For our analyses, we utilized rank-based nonparametric methods to overcome distributional requirements that could not be confirmed by the data. We used the Kruskal–Wallis (K.W.) nonparametric analysis of variance to test whether our treatment samples originated from the same distribution [27]. Before employing the K.W. test, we used the Fligner–Killeen test to confirm the assumption of equal dispersion among all the treatment levels [28]. We used Dunn’s test to assess pairwise differences between treatment levels [29] and Holm’s step-down method to control family-wise error rates [30]. Because the K.W. test is not one of the medians, we report the appropriate Hodges–Lehman estimators of the treatment levels [31], consistent with the K.W. test and analogous in interpretation to the other location parameters such as the median. We allowed a Type I error rate of α = 0.05 for significance testing.
The relationship between DBH and distance to neighboring trees was modeled using a Gompertz (or Logistic) asymptotic function.
Analyses were performed using the R statistical environment version 4.0.5 [32] and R libraries PMCMRplus [33] and DescTools [34].
Identifying the influence of provenance and spatial structure factors on tree growth was based on comparing the statistical analyses of DBH (Table 1). The original data were clustered into three classes in both cases, but slightly different approaches were used to isolate them.

2.3. Methodological Approach for Determining the Influence of a Provenance Factor

When determining the influence of provenance as a factor, the data were grouped into three classes based on a formal clustering procedure. The classification was based on productivity, as expressed through the average DBH of each half-sib (Table 1) based on the measure of their proximity to the average value for the entire sample.
The number of clusters (classes)—3 was established subjectively or in simplification: the average DBH is small (Class 1), medium (Class 2), and large (Class 3). The number of clusters, in this case, is optimal, and this approach is only formally subjective as it does not reduce the quality of clustering. It is widely used in practical forestry, for example, when preliminary assessing the level of fruiting, the state of natural and artificial regeneration, determining the vitality of forest stands, determining the degree of their damage by animals, insects, diseases, etc. [35,36,37]. Additionally, the number of objects (half-sibs) in the sample, which had to be divided into disjoint clusters, was relatively small (10), and therefore it would be inappropriate to increase the number of clusters by more than three (i.e., their fragmentation).
The solution to the clustering problem was reduced by determining for each kind the deviation of its average DBH from the average value of this parameter for the entire population (10 half-sibs), equal to 10.0 cm [12]. Half-sibs whose average DBH decreased more than 5% from the average DBH were assigned to the first class, those deviating no more than 5% in either direction (increasing or decreasing)—to the second, and finally those whose DBH was greater than the average DBH by 5% or more—to the third (Table 2).
As you can see, the deviation of the DBH of the species presented in the central, second class, is insignificant and amounts, as we initially assumed, no more than 5% of the average diameter for the entire population (701 trees), while for the first and third classes these deviations (“−” and “+”) more than 10%, sometimes much more (Table 2). These data, however, are only preliminary evidence of the legitimacy of clustering, i.e., there is a need to confirm it using better statistics. At the same time, it will answer the question of the influence of provenance on DBH.

2.4. Methodological Approach for Determining the Influence of the Factor of Spatial Structure

In our argument concerning spatial structure, we work with the concept of “growth area” [38,39]. It corresponds to the “feeding area” [40,41] and conditionally implies the space available to a particular tree and not used by the neighboring trees.
When determining the influence of the spatial structure (defined as the average distance to the nearest four trees in the directions of the cardinal points), the data on DBH were also divided into three classes. That was done for better comparison with the provenance factor. We used the data obtained earlier [12] and presented it in Table 3. Moreover, the partitioning of a set of objects (126 randomly selected trees), in contrast to clustering by provenance factor, was more consistent with the understanding of automatic classification, which excluded any subjectivity.
The following mathematical manipulation was applied. The minimum value was subtracted from the maximum value of the mean distance between trees, and class boundaries were established by dividing the difference by three. It was not difficult to establish the middle of the class interval. Additionally, the actual average distance between trees in each class was determined (Table 4).
Note that the last two parameters, which have a similar meaning but are determined in different ways (empirically and based on statistical processing of actual data), practically do not differ. This testifies to the legitimacy of the methodological approach used in the cluster classification of objects in studying the significance of the role of the spatial structure. Next, the statistical parameters DBH were determined for each class (Table 4). Thus, dividing the sample into three classes and grouping each class tree with a similar spatial structure (average distance to neighbors) allowed us to identify its effect on productivity (DBH) and compare it with the influence of the provenance factor.

3. Results

3.1. The Influence of the Provenance Factor

Dunn’s test of three classes of DBH by provenance factor showed a statistically significant difference between the DBH of trees of all three classes (Table 5), even in comparable pairs of classes with a close DBH, for example, Class 1/Class 2 and Class 2/Class 3.
All three classes have a wide range of values. Class 2 is represented by the largest sample (371 trees), the mean value of the DBH, which is 10.08 cm (Table 6). However, the smaller samples, represented by Class 1 (150 trees) and Class 3 (180 trees), are still sufficient to obtain representative data [42].

3.2. The Influence of the Spatial Structure

The analyzed feature also depends on the spatial structure, which the Kruskal–Wallis test confirmed. However, Dunn’s test for pairwise differences between all pairs revealed that only Class 1 compared to Class 2 (Class 1/Class2) and Class1/Class3 differ significantly (Table 7).
Class 1 is represented by the largest number of trees (68 trees), with the smallest average DBH. Class 2 has the greatest range in DBH. Classes 2 and 3, despite the difference in mean DBH (Table 8), do not differ statistically (Table 7). The range of variation in class 3 is the smallest.
The influence of the spatial features is superimposed on the phenotypic manifestation of genotypes (average DBH) characterizing certain half-sibs. We performed a correlation analysis (Table 9) using the tree DBH and spatial distribution of the neighboring trees (Table 1 and Table 3).
The relationship between the analyzed indicators is weak if the correlation coefficient is less than ±0.3; from ±0.3 to ±0.7—the relationship is the moderate strength; if the correlation coefficient is greater than ±0.7—the relationship is strong [42]. In our case, the relationship between the average stem DBH and the average distance to the nearest trees is strongly positive. Moreover, with the distances in specific directions, it has a moderate positive strength (Table 9).

3.3. Comparison of the Influence of Provenance Factor and Spatial Structure

We examined the differences between the DBH classes and the breadth of the variation within each class (Figure 4 and Figure 5). Using these data, we can make a preliminary conclusion about the presence (or absence) of dependence of DBH variation on factors of provenance and spatial structure and compare their influence on DBH.

3.4. Most Productive Phenotypes

Nine trees ranging from 20.3 to 23.7 cm in DBH and 18.0 to 22.4 m were identified (Table 10). While predictions of near-term growth can be challenging, it seems to us that it is unlikely to slow down. If we compare the estimated features of these 29-year-old trees with the tabular ones [43], we will see that they correspond to 55–60-year-old natural stands in diameter and 50–60 years old in height, growing in conditions of the 1st class of stand quality.
Half (55.5%) of the most productive phenotypes belonged to half-sib No. 17. Interestingly, the energy of the best trees is enough not only for intensive growth but also for frequent fruiting—some of them have recorded the formation of cones for two consecutive years. However, among these nine trees, not one is in the 1st class in terms of the spatial structure factor (Table 4).

4. Discussion

We note that the established significant variability of average DBH of half-sibs (Table 1) occurs in the absence of differences in edapho-climatic conditions, characteristics of the forest cultivated area, age, and type of planting material, timing, method of establishment, and maintenance of the plantation. The provenance factor’s influence is reflected by the membership of trees in a particular half-sib family. However, within each half-sib, a certain amount of phenotypic variability occurs. This variability might be due, among other things, to possible cross-fertilization with the influence of paternal components, which may differ genetically. Moreover, as noted above (see section “Introduction,” [19]), trees not only from Tisovik but also from stands of unknown origin in the Polish part of the Forest may participate in pollination. The closest location (subcompartment 453Aa) is 8.5 km from the refugium (Figure 6). If this is the case, then given the growing aridity of the climate, an increase in genetic mixing is regarded as a positive factor. However, the opposite process is also possible—with pollen from Tisovik pollinating the macrostrobils on trees in sub-compartment 453Aa.
As shown above, the parameters of phenotypic variability differ in different half-siblings, which may evolve into actual genetic differences. This potential for future genetic variability suggests the need for maximum conservation of each tree—a potential source of heritable valuable economic traits—in thinning.
Note that the biological features of fir fully contribute to this outcome. On the one hand, due to the pronounced taproot system of fir, the growing space of individual trees is relatively small. On the other hand, this species is characterized by exceptional shade tolerance [44,45,46]. These results suggest that managers reduce harvest impacts of future silvicultural and increase genetic opportunities by minimizing felling intensity at the cleaning stage. For example, in 2015, the marked tree in Figure 7 was initially designated for felling because of its reduced growth compared to the neighboring trees. However, since this experimental site was assigned to the UNESCO protected zone II, where any felling is prohibited, this tree was preserved. This “reprieve” was borne out by its subsequent performance similar to neighboring trees by 2021.
In this regard, we also note that, according to Goncharenko [47], the fir in Tisovik “despite its isolation and small size, has not yet experienced the influence of inbreeding, and although only 20 trees were preserved, as a result of crosses (cross-pollination) between them, one can obtain 13,122 genetically different offspring”.
Nevertheless, the issue of the manifestation of inbreeding in the Tisovik tract should be approached more carefully and, perhaps, not so optimistically [48]. Because the population, represented by only about 20 trees [7,12,47], is very small, and they live in a very small area. In addition, perhaps it is the effect of inbreeding that can explain the lowest survival rate in family No. 5—28.3 (only 73 trees out of 258 survived) [12], as well as the significant differences in diameter both between sibs and within sibs, as can be seen from our data. Although, in the latter case, the effect of inbreeding, even if it takes place, is manifested against the background of the natural process of biological differentiation of trees (Kraft classes or social classes of trees position in canopy [48]).
Moreover, for an exceptionally shade-tolerant fir, it is characteristic that not only trees of IV, but even sometimes Va of Kraft classes are not doomed to fail. Most of them are not in a state of gradual mortality. On the contrary, their condition corresponds to a live plant, only with very slow growth rates. In addition, this is precisely why the statistical analysis showed a large standard deviation and, accordingly, the coefficient of variation of DBH (Table 1, Table 2 and Table 4). In natural stands, when the dominant trees approach the climax stage, it is the representatives of the lower tiers that begin to take over their role [44,45,46].
To obtain a genetically richer silver fir population originating from the seeds of the planned vegetative seed plantation created by grafting, we must first select the best phenotypes from all sibs.
Therefore, when creating it, you should not focus only on the top nine trees identified at the moment. According to Barzdajn [49], who was involved in the restitution of silver fir in Sudety, clones should number at least 100. This conclusion implies that phenotypic studies must be repeated periodically to identify the best trees (Table 10). In our opinion, a frequency of 5 years will be necessary. These remeasurements will make it possible to record the emergence of new candidates for superior trees since, as we noted earlier, there is a pronounced positive dynamics of tree growth [12] (Figure 2 and Figure 3). Additionally, two more plots are awaiting their turn at this object. The most critical condition for the success of future vegetative seed plantations is the maximum representation of all the best phenotypes from all half-sibs families.
The closest positive correlation of DBH occurs with the distance to the nearest neighboring tree in the southern direction (0.69) and the lowest (0.45) in the eastern direction (Table 9). These results suggest that the creation of conditions for an increase in illumination in the stand without excessive felling of trees (i.e., an increase in the distance to a neighboring tree in the southern direction with the highest intensity of solar radiation) has a positive effect on the growth of trees (an increase in their DBH). As it turns out, A. alba, which is a classically shade-tolerant species, responds positively to improving light conditions under these conditions, even at a young age.
The relationship between the DBH and the average distance to neighboring trees (Figure 8) suggests a certain “threshold value” of the average distance between trees in a stand of a given age. Beyond this threshold, further thinning will not lead to an increase in the productivity of the stand nor an improvement in its general condition. At the moment, this threshold indicator is 2.3 m (Figure 8). However, this distance is not permanent and will change with the age of the stand.
Thus, 2.3 m is a “threshold,” which we believe determines the end of the influence of spatial structure on DBH increase. This factor ceases to be limiting [50,51]. Even though the curve of the graph continues to move very slightly upward, from a silvicultural point of view, a further increase in the average distance will lead to an unnecessary decrease in the density of the stand, a reduction in its overall growth and volume, and underutilization of forest area by root systems of trees. Furthermore, an excessive increase in the sparseness of the stand negatively affects the light regime of the shade-tolerant fir. A simple calculation shows that with an increase in the average distance between trees from 2.3 m to 2.5 m (i.e., when the red cross in Figure 8 moves to the right), it is necessary to additionally remove about 300 trees from a hectare, which will lead to a decrease in stand density by 0.1 units.
The same cautious approach to tree cutting should be followed in the future. Of course, the parameters of the spatial structure will change, and this is the subject of future research. However, the principles of silvicultural care remain as follows: (a) only obvious outsiders are cut down (trees at the left-hand side of the distribution curve), and (b) tree stand care is carried out using more frequent, low-intensity entries rather than infrequent intensive thinnings.
In this case, the presence of individuals with an increased DBH and indicators of spatial structure exceeding the “threshold” can only be interpreted as a consequence of the influence of a provenance factor. As shown in Table 10, in seven out of nine trees (78%) with the largest DBH (candidates for plus trees), the average distance to neighbors exceeds 2.3 m. A genetic study of these trees, planned for the future, will allow a more substantiated statement on the validity of this assumption.
Following the methodological approach of Rogozin [41], we use half the average distance between trees (1.15 m) as the radius of a sufficient “feeding area” [40,41] or “growth area” [38,39]. As a result, we get that the growth area of each tree is 4.15 m2, and the density of the stand is 2410 pcs/ha (where it is 3186 pcs/ha at present). It is necessary to approach this indicator, as already mentioned, by selecting the weak intensity according to the grass-root method, removing especially the trees of the 2 and 4 cm diameter classes. In this way: (1) the biology of the species (shade tolerance) is taken into account, (2) the vertical structure is formed, and (3) the gene pool of A. alba is preserved.
In a situation where a pronounced differentiation is currently taking place in the forest stand (Table 1), the future of the dominant trees, especially advantage tree candidates, is of interest. Can their growth be improved? It was found that there is an unequal dependence of the change in DBH on the distance to the nearest tree in the direction of the different world sides (Figure 9).
The starting point for evaluating this dependence is the intersection of the two perpendiculars on the axis of this diagram. The vertical perpendicular corresponds to a DBH of 12 cm, representing the diameter class above the average stand DBH (9.97 cm, [12]). The horizontal line corresponds to the “threshold” of the spatial structure (see above).
The thinning required to increase the DBH, expressed as the distance to the nearest tree in the various directions, takes the following form (in decreasing order): north, south, west, and east. That is, by focusing on the removal of trees in the rows (north-south direction) and at the same time cutting less intensively in the direction between the rows (west-east), it is possible to influence the increase in DBH of the dominant trees. This opportunity arises from the trees’ initial very dense arrangement in the rows (1 m, and sometimes less, Table 3).
The complexity and ambiguity determining the influence of the spatial structure are illustrated by this analysis of the randomly located sample (Figure 10).
As can be seen, the planting lines are still visible. However, the initially regular (systematic) placement of trees according to the 1.3 × 1 m design becomes uneven. Biogroups of trees (highlighted in green) and gaps in the forest stand (highlighted in orange) begin to form. Note that from a bird’s eye view, the unevenness of the canopy of the crowns, due to their development, is not yet so apparent (Figure 2). What may be the cause of the uneven spatial structure of the stand? First of all, we note the possible accretion of fir roots [52] and the formation, thus, of microgroups, each of which has a common root system. Some plants died due to late spring frosts in the early years (see Korczyk [53]). Perhaps the negative influence of excessive solar radiation was also uneven in different parts of the plot. In addition, this non-uniformity in the distribution of plants persisted subsequently (we did not find any information in the literature about the addition of crops). After thinning, it may have even increased somewhat under the influence of the subjective factor (selection of trees for cleanings). However, the root causes of this phenomenon (uneven distribution of trees) are not yet known [41]. Therefore, we can only note the appearance of fundamentally new attributes of the morphogenesis of the stand: biogroups (microcenoses) and gaps in the stand.
From the analysis of the fragment of the spatial model (Figure 10), it follows that the trees with neighbors that are behind in growth or have a greater distance to them also have a larger DBH (Nos. 1, 4). Conversely, the DBH of the focus trees decreases with decreasing distance to neighboring trees or an increase in the neighboring trees’ DBH (Nos. 2, 3). These two situations are well explained by the weakening and strengthening of intraspecific competition. However, then, do we explain the same DBH of the trees Nos. 7–14 growing in the biogroup? Perhaps there is already the partnership effect (group effect) [41] (see Figure 11). In addition, in this context—perhaps there is no need to strive for an even distribution of trees [54]? In nature, A. alba is usually found in an admixture of the main tree species in groups because it has relatively heavy seeds [46]. In addition, as it turns out, it grows better with artificial reproduction in groups; in particular, it is highly resistant to biotic and abiotic factors [55].
And what influenced the insignificant DBH of the trees Nos. 5 and 6 growing in the gaps? After all, the competition factor is not the limiting influence at present. Perhaps they were suppressed before the late cleaning and still cannot rebuild? Or is this perhaps the result of genetic differentiation? Approximately the same questions arise in analyzing any other fragment of a half-sib plantation.
It must be admitted that the above analysis is objectively a simplified approach that does not take into account other complex influences that are challenging to determine and not consistently observable by managers and scientists (allelopathy, root competition, local specificity of the soil, the state of the plants before planting and late cleaning, etc.). According to Rogozin [41], the proportion of unknown factors that determine the size of trees in old pine stands can reach 60%.
When assessing the influence of the provenance factor on the diameter, it was shown that all three classes are statistically significantly different from each other in terms of the p-level of significance. In addition, there is a pronounced dependence of the DBH on the provenance factor (Figure 4).
A different situation occurs when assessing the influence of spatial structure. While classes 1 and 2, as well as 1 and 3, differ significantly from each other (p < 0.01), the difference between classes 2 and 3 is not significantly significant (p = 0.12). Nevertheless, we note that the tendency to increase in diameter from class 2 to class 3 is visually apparent (Figure 5).
In light of the above, there is a reason to once again turn to the conclusion illustrated in Figure 8. The distance between the trees (the density of the stand), which provides them with favorable conditions at a given age is 2.3 m, the middle of the class interval No. 2 (Table 4). The silvicultural interpretation of the statistically insignificant difference between classes 2 and 3 is the same: a borderline acceptable distance between trees—“threshold”—exists. We found the same result using two different methods, supporting our conclusions.

5. Conclusions

The study results give reason to believe that at this particular age, provenance has a more significant influence on the growth of the trees in DBH throughout the whole range of its variation than has the spatial structure. However, it was also found that the importance of the last factor should not be discounted, especially in the range of small and medium distances between trees (classes 1 and 2), where it can determine their parameters.
Free pollination (father unknown) of a particular mother tree in the Tisovik tract is the basis for the emergence of the half-sibs we studied. In this respect, each of the trees on the plantation of European silver fir, even within the same family, is genetically differentiated and may differ in characteristics to one degree or another. Therefore, to maximize future silvicultural and genetic opportunities, we should preserve as many trees as possible.
The nine best trees were identified, which are the prototypes of plus trees. Their seeds and shoots can already find a practical use for the early creation of a seed base of autochthonous silver fir. By starting this work now, in an experimental stand of about 30 years, we will accelerate the emergence of seed plantations by at least 50 years.
Family No. 17 is of particular interest for further genetic study. A significant number of specimens were identified in this half-sib, with DBH significantly exceeding the average of the plantation. However, this does not mean we should concentrate only on this half-sib. On the contrary, all semi-sibs should be maximally represented in the seed base in order to maintain or even enhance the genetic diversity of future autochthonous silver fir forests.
The actual planting density at the time of the study was 3186 plants/ha, while the threshold value determined by us is 2410 plants/ha. This particular hold can be used in practical forestry when thinning silver fir trees at about 25–30 years. The stand should be “brought” to the threshold value gradually, performing regular fellings of low intensity and removing only obvious outliers.
In light of the catastrophic situation with Norway spruce in Białowieża Forest, studies on various aspects of growing local silver fir as an alternative are of undoubted scientific and practical value.

Author Contributions

Conceptualization, A.M., T.O. and U.K.; methodology, A.M., U.K., D.B., J.N. and T.O.; software, U.K. and D.B.; validation, D.B., J.N., T.O., W.K.M. and T.H.; formal analysis, A.M., U.K., D.B., J.N. and T.O.; investigation, A.M., T.O. and J.N.; resources, A.M. and T.O.; writing—original draft preparation, A.M., U.K., T.O. and J.N.; editing, T.O., J.N., D.B., W.K.M. and T.H.; supervision, A.M., D.B., J.N. and T.O.; project administration, A.M. and T.O.; funding acquisition, A.M. All authors have read and agreed to the published version of the manuscript.

Funding

The study was carried out within the framework: WZ/WB-INL/2/2021 and financed from the science funds from the Ministry of Science and Higher Education in Poland.

Data Availability Statement

No publicly archived dataset.

Acknowledgments

The authors thank student Adam Kaliszuk for participation in fieldwork and primary data processing. Some material from the engineering thesis of Adam Kaliszuk “Analysis of the breeding status of a 23-year-old ancestral cultivation of silver fir (Abies alba Mill.) in the Hajnówka Forest District” at Bialystok University of Technology, Poland, was referenced in the article with full consent of the main author.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the study’s design, in the collection, analyses, or interpretation of data, in the writing of the manuscript, or in the decision to publish the results.

References

  1. Brzeziecki, B.; Hilszczański, J.; Kowalski, T.; Łakomy, P.; Małek, S.; Miścicki, S.; Modrzyński, J.; Sowa, J.M.; Starzyk, J.R. Problem masowego zamierania drzewostanów świerkowych w Leśnym Kompleksie Promocyjnym “Puszcza Białowieska”. Sylwan 2018, 162, 373–386. (In Polish) [Google Scholar]
  2. Nowakowska, J.A.; Hsiang, T.; Patynek, P.; Stereńczak, K.; Olejarski, I.; Oszako, T. Health assessment and genetic structure of monumental Norway spruce trees during a bark beetle (Ips typographus L.) outbreak in the Białowieża Forest District, Poland. Forests 2020, 11, 647. [Google Scholar] [CrossRef]
  3. Szafer, W. Jodła w Puszczy Białowieskiej [Fir in the Białowieza Primeval Forest]. Sylwan 1920, 38, 65–74. (In Polish) [Google Scholar]
  4. Paczoski, J. Lasy Białowiezy [Forests of Białowieża]; Państwowa Rada Ochrony Przyrody: Poznań, Poland, 1930. (In Polish)
  5. Karpiński, J.J. Białowieża; Drukarnia państwowa w Płocku: Warszawa, Poland, 1947. (In Polish)
  6. Korczyk, A.F.; Kawecka, A.; Martysevič, V.V.; Strelkov, A.Z. European fir (Abies alba Mill.) natural stand in Białowieża Primeval Forest. Pr. Inst. Badaw. Leśnictwa 1997, 837, 27–62. [Google Scholar]
  7. Ermochin, M.V.; Barsukova, T.L.; Knys, N.V.; Byśko, V.E.; Bernackij, D.I. Dinamika i sostojanie populjacii pichty beloj v uročišče Tisovik [Dynamics and state of the silver fir’s population in the Tisovik tract]. In Belovezhskaia Pushcha. Research Compendium of Scientific Articles; Buryj, A., Arnolbik, V., Cherkas, N., Bunevich, A., Eds.; Alternativa: Brest, Belarus, 2016; pp. 65–88. (In Russian) [Google Scholar]
  8. Korczyk, A.F. Ocena wartości genetycznej i hodowlanej naturalnych populacji jodły pospolitej (Abies alba Mill.) ze wschodniego zasięgu w Polsce [Assessment of the genetic and breeding value of natural populations of silver fir (Abies alba Mill.) from theeastern range in Poland]. Zesz. Nauk. Akad. Rol. Im. H. Kołłątaja W Krakowie 1999, 61, 155–170. (In Polish) [Google Scholar]
  9. Cherkas, E.V.; Morozov, O.V. Influence of the oak red (Quercus rubra) on components of wood plantings of the SNPE «Belovezhsky dense forest». In Problemy Lesowedenija I Lesowodstwa; INL NAN Belarusi: Gomel, Belarus, 2012; Volume 72, pp. 150–156. (In Russian) [Google Scholar]
  10. Cherkas, E.V.; Morozov, O.V. Population of a red oak (Quercus rubra L.) in the Bialowieza Forest. In Ekologija, Lesowodstwo i Lesoohotničje Hozjajstwo; Belarusian State Technological University: Minsk, Belarus, 2012; Volume 1, pp. 127–129. (In Russian) [Google Scholar]
  11. Marozau, A.; Mielcarek, M.; Krok, G.; Paluch, R.; Chiliński, K. European silver fir—An alternative for the dying Norway spruce in Białowieża Forest? Folia For. Pol. 2021, 63, 150–166. [Google Scholar] [CrossRef]
  12. Marozau, A.; Kotsan, U.; Kalishuk, A. Reintroduction of the European silver fir (Abies alba Mill.) in Białowieża Forest. Balt. For. 2021, 27, 121–131. [Google Scholar] [CrossRef]
  13. Strelkov, A.Z.; Romaniuk, I.G.; Dworak, L.E. Sostojanie i perspektiwy wosproizwodstwa i ochrany pichty beloj (Abies alba Mill.) w Belovežskoj pušče. In Sochranenie Biologičeskogo Raznoobrazija Lesow Belovežskoj Pušči [Conservation of the Biological Diversity of the Forests of Belovezhskaia Pushcha]; Lučkov, A.I., Tolkač, W.N., Parfenov, W.I., Sawickij, B.P., Januško, A.D., Eds.; Izdanie: Kameniuki-Minsk, Belarus, 1996; pp. 122–130. (In Russian) [Google Scholar]
  14. Parfenov, V.I.; Savčuk, S.S. Abies alba Mill. Red Book of the Republic of Belarus. Rare and Threatened Species of Wild Plants; Encyklapedyja im. P. Brouki; Belarusian State Technological University: Minsk, Belarus, 2015; pp. 25–26. (In Russian) [Google Scholar]
  15. Gončarenko, G.G.; Padutov, A.E.; Padutov, V.E.; Silin, A.E.; Bolsun, S.I. Opredelenie genetičeskogo raznoobrazija osnownych lesoobrazujuščich i redkich drewesnych vidov GNP «Belovežskaja pušča». In Sochranenie Biologičeskogo Raznoobrazija Lesow Belovežskoj Pušči [Conservation of the Biological Diversity of the Forests of Belovezhskaia Pushcha]; Lučkov, A.I., Tolkač, W.N., Parfenov, W.I., Sawickij, B.P., Januško, A.D., Eds.; Izdanie: Kameniuki-Minsk, Belarus, 1996; pp. 141–160. (In Russian) [Google Scholar]
  16. Longauer, R. Genetic variation of European silver fir (Abies alba Mill.) in the Western Carpathians. J. For. Sci. 2001, 47, 429–438. [Google Scholar]
  17. Kempf, M.; Zarek, M.; Paluch, J. The pattern of genetic variation, survival and growth in the Abies alba Mill. Population within the introgression zone of two refugial lineages in the carpathians. Forests 2020, 11, 849. [Google Scholar] [CrossRef]
  18. Mihai, G.; Alexandru, A.M.; Stoica, E.; Birsan, M.V. Intraspecific growth response to drought of Abies alba in the Southeastern Carpathians. Forests 2021, 12, 387. [Google Scholar] [CrossRef]
  19. Pawlaczyk, E.M.; Bączkiewicz, A.; Buczkowska, K.; Bobowicz, M.A. Genetic variation of silver fir progeny from Tisovik Reserve population determined via microsatellite and isozyme markers. Folia For. Pol. 2017, 59, 45–58. [Google Scholar] [CrossRef] [Green Version]
  20. Mejnartowicz, L. Cisovka—The relic population of Abies alba and its relationship to man-made silver-fir stands in Białowieża Primeval Forest. Acta Soc. Bot. Pol. 1996, 65, 319–328. [Google Scholar] [CrossRef] [Green Version]
  21. Oliver, C.D.; Larson, B.A. Forest Stand Dynamics, Update Edition; John Wiley & Sons (Copyright now held by CD Oliver & BC Larson): Hoboken, NY, USA, 1996. [Google Scholar]
  22. Rossi, S.; Morin, H.; Gionest, F.; Laprise, D. Spatially explicit structure of natural stands dominated by black spruce. Silva Fenn. 2013, 47, 1–14. [Google Scholar] [CrossRef] [Green Version]
  23. Goreaud, F.; Loussier, B.; Marie Ange, N.B.; Allain, R. Simulating Realistic Spatial Structure for Forest Stands: A Mimetic Point Process; Interdisciplinary Spatial Statistics Workshop: Paris, France, 2004. [Google Scholar]
  24. Wang, H.; Zhang, G.; Hui, G.; Li, Y.; Hu, Y.; Zhao, Z. The influence of sampling unit size and spatial arrangement patterns on neighborhood-based spatial structure analyses of forest stands. For. Syst. 2016, 25, 56. [Google Scholar] [CrossRef] [Green Version]
  25. Mason, W.L.; Connolly, T.; Pommerening, A.; Edwards, C. Spatial structure of semi-natural and plantation stands of Scots pine (Pinus sylvestris L.) in Northern Scotland. Forestry 2007, 80, 567–586. [Google Scholar] [CrossRef] [Green Version]
  26. Corral-Rivas, J.J. Models of Tree Growth and Spatial Structure for Multi-Species, Uneven-Aged Forests in Durango (Mexico)—Cuvillier Verlag. Available online: https://cuvillier.de/de/shop/publications/2111 (accessed on 23 September 2021).
  27. Kruskal, W.H.; Wallis, W.A. Use of ranks in one-criterion variance analysis. J. Am. Stat. Assoc. 2012, 47, 583–621. [Google Scholar] [CrossRef]
  28. Conover, W.J.; Johnson, M.E.; Johnson, M.M. A comparative study of tests for homogeneity of variances, with applications to the outer continental shelf bidding data. Technometrics 1981, 23, 351–361. [Google Scholar] [CrossRef]
  29. Dunn, O.J. Multiple comparisons using rank sums. Technometrics 1964, 6, 241–252. [Google Scholar] [CrossRef]
  30. Holm, S. A simple sequentially rejective multiple test procedure. Scand. J. Stat. 1979, 6, 65–70. [Google Scholar]
  31. Hodges, J.L.; Lehmann, E.L. Estimates of location based on rank tests. Ann. Math. Stat. 1963, 34, 598–611. [Google Scholar] [CrossRef]
  32. R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing; EEA: Vienna, Austria, 2021.
  33. Signorell, A.; Aho, K.; Alfons, A.; Anderegg, N.; Aragon, T.; Arachchige, C.; Arppe, A.; Baddeley, A.; Barton, K.; Bolker, B.; et al. Borchers DescTools: Tools for Descriptive Statistics; R Package Version 0.99.43; R Core Team: Vienna, Austria, 2021; pp. 1–648. [Google Scholar]
  34. Pohlert, T. PMCMRplus: Calculate Pairwise Multiple Comparisons of Mean Rank Sums Extended; R Package Version 1.9.2; R Core Team: Vienna, Austria, 2021; pp. 1–207. [Google Scholar]
  35. Haze, M. (Ed.) Zasady hodowli hodowli lasu. In Pincipes of Siviculture; Centrum Informacyjne Lasów Państwowych: Warszawa, Poland, 2012. (In Polish) [Google Scholar]
  36. Sierota, Z. Criteria and methods of estimation of trees and stands health state. Pr. Inst. Badaw. Leśnictwa 1998, 854, 28. (In Polish) [Google Scholar]
  37. Jaworski, A.; Paluch, J. Charakterystyka cech morfologicznych jodeł w drzewostanach o strukturze przerębowej Beskidów. Leśne Pr. Badaw. 2007, 3, 7–31. (In Polish) [Google Scholar]
  38. Nagimov, Z.J. Evaluation of methods for determining tree growth areas. Lesa Ural. I Chozajstwo W Nich 1999, 19, 82–98. (In Russian) [Google Scholar]
  39. Kotsan, V.V. Classification of trees on the basis of spatial structure in appointments thinning. Tr. BGTU 2015, 1, 24–27. (In Russian) [Google Scholar]
  40. Rogozin, M.V. Nutrion area of tree: Analysis of methods. Bull. Sci. Pract. 2018, 4, 34–37. (In Russian) [Google Scholar]
  41. Rogozin, M.V. The Structure of the Forest: Competition or Partnership? Perm State University: Perm, Russia, 2019; p. 222. (In Russian) [Google Scholar]
  42. Bühl, A.; Zöfel, P. SPSS: Iskusstvo Obrabotki Informatsii. Analiz Statisticheskikh Dannykh I Vosstanovlenie Skrytykh Zakonomernostey [SPSS: The Art of Information Processing. Statistical Data Analysis and Mining]; DiaSoftYuP: St Petersburg, Russia, 2005. (In Russian) [Google Scholar]
  43. Szymkiewicz, B. Tablice Zasobności I Przyrostu Drzewostanów Ważniejszych Gatunków Drzew Leśnych [Stands Abundance and Growth Tables of The Most Important Forest Tree Species]; PWRiL: Warszawa, Poland, 2001; p. 179. (In Polish) [Google Scholar]
  44. Jaworski, A.; Zarzycki, K. Ekologia. In Jodła Pospolita Abies Alba Mill. Nasze Drzewa Leśne; Białobok, S., Ed.; Państwowe Wydawn. Nauk: Warszawa-Pozanań, Poland, 1983; Volume 4, pp. 317–430. (In Polish) [Google Scholar]
  45. Bernadzki, E. Jodła pospolita. In Ekologia—Zagrożenia—Hodowla; Fonder, W., Ed.; PWRiL: Warszawa, Poland, 2008; p. 209. (In Polish) [Google Scholar]
  46. Jaworski, A. Hodowla Lasu. Charakterystyka Hodowlana Drzew I Krzewów Leśnych; PWRiL: Warszawa, Poland, 2019; Volume 3, p. 696. (In Polish) [Google Scholar]
  47. Gonczarenko, G.G. Analysis of Population Genetic Resources (on the Example of White Fir in Belarus); Gomelskij Gosudarstwennyj Universitet Imeni F. Skoriny: Gomel, Belarus, 2002; p. 112. ISBN 978-985-6234-78-4. (In Russian) [Google Scholar]
  48. Jaworski, A. Hodowla Lasu. Pielęgnowanie Lasu; PWRiL: Warszawa, Poland, 2020; Volume 2, p. 400. (In Polish) [Google Scholar]
  49. Bardzain, W. Strategia restytucji jody pospolitej (Abies alba Mill.) w Sudetach. Sylwan 2000, 2, 63–77. (In Polish) [Google Scholar]
  50. Zeide, B. Optimal stand density: A solution. Can. J. For. Res. 2004, 34, 846–854. [Google Scholar] [CrossRef]
  51. Assmann, E. Waldertragskunde. Organische Produktion, Struktur, Zuwachs und Ertrag Von Waldbeständen; BLV Verlagsgesellschaft: München, Germany, 1961; p. 490. [Google Scholar]
  52. Kobendza, R. O zarastaniu ściętych pni jodłowych. Acta Soc. Bot. Pol. 1932, 9, 123–153. [Google Scholar] [CrossRef] [Green Version]
  53. Korczyk, A.F. The ancestral conservative tillage of silver fir in the ‘Tisovik’ reserve of the Białowieża Primeval Forest. For. Res. Pap. 2015, 76, 153–167. [Google Scholar] [CrossRef] [Green Version]
  54. Tjabera, A.P. Issues of territorial distribution in pine stands. Lesn. Žurnal 1980, 5, 5–8. [Google Scholar]
  55. Będkowski, M.; Czacharowski, M.; Opanowska, J.; Patro, S.; Urbańczyk, E.; Bielak, K. Effect of planting method on the growth and silvicultural quality of silver fir (Abies alba Mill.) in the Bieszczady Mountains. Sylwan 2021, 165, 542–553. [Google Scholar] [CrossRef]
Figure 1. Plot No. 1 of a half-sib silver fir plantation (photo by M. Kaczyński).
Figure 1. Plot No. 1 of a half-sib silver fir plantation (photo by M. Kaczyński).
Forests 13 00831 g001
Figure 2. General view of plot No. 1 of the half-sib silver fir plantation from the lift when collecting cones in 2021 (photo by K. Wilamowski).
Figure 2. General view of plot No. 1 of the half-sib silver fir plantation from the lift when collecting cones in 2021 (photo by K. Wilamowski).
Forests 13 00831 g002
Figure 3. Fruiting of individual trees on plot No. 1 of the half-sib silver fir plantation in 2021 (photo by K. Wilamowski).
Figure 3. Fruiting of individual trees on plot No. 1 of the half-sib silver fir plantation in 2021 (photo by K. Wilamowski).
Forests 13 00831 g003
Figure 4. The relationship between tree DBH and the provenance factor. The lower and upper hinges correspond to the first and third quartiles (the 25th and 75th percentiles), and the centerline corresponds to the median. The upper whisker extends from the upper hinge to the largest value no further than 1.5*IQR from the hinge (where IQR is the interquartile range or distance between the first and third quartiles). Likewise, the lower whisker extends from the lower hinge to the smallest value no further than 1.5*IQR from the hinge. The different lower case letters (a, b, c) above the means of the treatments (Class 1, 2, 3) are significantly different at p = 0.05. Points on the plot were offset to avoid overplotting of identical values.
Figure 4. The relationship between tree DBH and the provenance factor. The lower and upper hinges correspond to the first and third quartiles (the 25th and 75th percentiles), and the centerline corresponds to the median. The upper whisker extends from the upper hinge to the largest value no further than 1.5*IQR from the hinge (where IQR is the interquartile range or distance between the first and third quartiles). Likewise, the lower whisker extends from the lower hinge to the smallest value no further than 1.5*IQR from the hinge. The different lower case letters (a, b, c) above the means of the treatments (Class 1, 2, 3) are significantly different at p = 0.05. Points on the plot were offset to avoid overplotting of identical values.
Forests 13 00831 g004
Figure 5. Dependence of the DBH on the spatial structure factor. The lower and upper hinges correspond to the first and third quartiles (the 25th and 75th percentiles), and the centerline corresponds to the median. The upper whisker extends from the upper hinge to the largest value no further than 1.5*IQR from the hinge (where IQR is the interquartile range or distance between the first and third quartiles). Likewise, the lower whisker extends from the lower hinge to the smallest value no further than 1.5*IQR from the hinge. The different lower case letters (a, b) above the means of the treatments (Class 1, 2, 3) are significantly different at p = 0.05. Points on the plot were offset to avoid overplotting of identical values.
Figure 5. Dependence of the DBH on the spatial structure factor. The lower and upper hinges correspond to the first and third quartiles (the 25th and 75th percentiles), and the centerline corresponds to the median. The upper whisker extends from the upper hinge to the largest value no further than 1.5*IQR from the hinge (where IQR is the interquartile range or distance between the first and third quartiles). Likewise, the lower whisker extends from the lower hinge to the smallest value no further than 1.5*IQR from the hinge. The different lower case letters (a, b) above the means of the treatments (Class 1, 2, 3) are significantly different at p = 0.05. Points on the plot were offset to avoid overplotting of identical values.
Forests 13 00831 g005
Figure 6. Localization of the silver fir refugium and the nearest fir stand in the Polish part of the Białowieża Forest (photo by W. Klimiuk).
Figure 6. Localization of the silver fir refugium and the nearest fir stand in the Polish part of the Białowieża Forest (photo by W. Klimiuk).
Forests 13 00831 g006
Figure 7. Silver fir, which was lagging in growth in 2015 (it was then marked in green, which means an initial assignment for harvest, but not felled). Compare this fir in 2021 to neighboring trees that were not assigned for felling 6 years ago (photo by A. Marozau).
Figure 7. Silver fir, which was lagging in growth in 2015 (it was then marked in green, which means an initial assignment for harvest, but not felled). Compare this fir in 2021 to neighboring trees that were not assigned for felling 6 years ago (photo by A. Marozau).
Forests 13 00831 g007
Figure 8. Dependence of the DBH on the average distance to neighboring trees. The line was fitted using a Logistic (or Gompertz) nonlinear relationship that allows for estimation of the asymptote of 2.3 m (lower 95% bound, upper 95% bound).
Figure 8. Dependence of the DBH on the average distance to neighboring trees. The line was fitted using a Logistic (or Gompertz) nonlinear relationship that allows for estimation of the asymptote of 2.3 m (lower 95% bound, upper 95% bound).
Forests 13 00831 g008
Figure 9. Dependence of DBH on the distance to neighboring trees according to the geographical orientation.
Figure 9. Dependence of DBH on the distance to neighboring trees according to the geographical orientation.
Forests 13 00831 g009
Figure 10. Arrangement of trees placement in an arbitrarily selected stand’s fragment (differences in the DBH are reflected by different sizes of circles). Drawing by A. Kaliszuk.
Figure 10. Arrangement of trees placement in an arbitrarily selected stand’s fragment (differences in the DBH are reflected by different sizes of circles). Drawing by A. Kaliszuk.
Forests 13 00831 g010
Figure 11. Biogroup of old firs in the Polish part of Białowieża FD (subdistrict Grudki, subcompartment 498 Ci) (photo by A. Marozau).
Figure 11. Biogroup of old firs in the Polish part of Białowieża FD (subdistrict Grudki, subcompartment 498 Ci) (photo by A. Marozau).
Forests 13 00831 g011
Table 1. Parameters of the DBH of the examined half-sibs.
Table 1. Parameters of the DBH of the examined half-sibs.
Half-Sib FamilyMean and Standard Deviation of the DBH, cmMedianLower Confidence IntervalUpper Confidence Interval
47.71 (3.36)7.756.29.3
38.46 (3.55)8.57.19.9
218.52 (4.68)8.54.9511.9
58.63 (3.50)8.57.659.35
168.64 (2.40)8.686.710.5
119.97 (3.28)9.88.711.05
1510.05 (4.15)9.89.210.45
210.13 (4.04)10.059.410.7
1211.25 (3.23)10.998.9513.45
1711.18 (4.52)11.110.411.8
Table 2. Data clustering in the context of the provenance factor.
Table 2. Data clustering in the context of the provenance factor.
ClassHalf-Sibs (Deviation from the Mean Value of the DBH for the Entire Sample in %) DBH
Mean and
Standard Deviation, cm
Coefficient of Variation,%
14 (−22.9), 3 (−15.4),
21 (−14.8), 5 (−13.7),
16 (−13.6)
8.44 (3.49)41
211(−0.3), 15 (+0.5),
2 (+1.3)
10.08 (4.02)40
312 (+12.5), 17 (+11.8)11.19 (4.43)37
Table 3. Basic statistical indicators are used to group data by class in determining the influence of the factor of spatial structure [12].
Table 3. Basic statistical indicators are used to group data by class in determining the influence of the factor of spatial structure [12].
Indicator Mean and Standard DeviationMinimum ValueMaximum ValueCoefficient of Variation, %
Ln a2.13 (1.09)0.776.1552
Ls b1.94 (0.88)0.696.3046
Le c1.89 (0.72)0.583.7538
Lw d1.81 (0.73)0.594.2641
Lav e1.93 (0.49)1.103.4226
a—distance to the nearest tree in a northerly direction (distance in a row), m; b—distance to the nearest tree in a southerly direction (distance in a row), m; c—distance to the nearest tree in the east (distance between rows), m; d—distance to the nearest tree in the west direction (distance between rows), m; e—average distance to the nearest trees, m.
Table 4. Data clustering in the context of the spatial structure factor.
Table 4. Data clustering in the context of the spatial structure factor.
ClassDistanceDBH
Class BoundariesThe Middle of the Class IntervalMean and Standard Deviation, mCoefficient of Variation, %Mean and Standard Deviation, cmCoefficient of Variation, %
11.10–1.871.491.59 (0.19)129.6 (4.1)43
21.88–2.652.262.18 (0.20)913.4 (4.9)37
32.66–3.423.043.07 (0.21)716.9 (4.2)25
Table 5. The p-value result of pairwise comparisons of classes by provenance (for α = 0.05).
Table 5. The p-value result of pairwise comparisons of classes by provenance (for α = 0.05).
Comparison Groupsp-Value
Class 1/Class 20.00015
Class 1/Class 30.0000000072
Class 2/Class 30.00227
Table 6. The volume of classes of the provenance factor and the amplitude of indicators of DBH in them.
Table 6. The volume of classes of the provenance factor and the amplitude of indicators of DBH in them.
ClassN, pcsMode, cmHodges-Lehman Estimator Lower Confidence IntervalUpper Confidence Interval
1150multi *8.47.818.99
237199.99.4910.35
31801111.110.4511.75
* distribution having several modes.
Table 7. The p-value result of pairwise comparisons of classes by spatial structure (for α = 0.05).
Table 7. The p-value result of pairwise comparisons of classes by spatial structure (for α = 0.05).
Comparison Groupsp-Value
Class 1/Class 20.000338039927
Class 1/Class 30.000000019719
Class 2/Class 30.12146648694
Table 8. The volume of classes of the spatial structure factor and the amplitude of indicators of DBH in them.
Table 8. The volume of classes of the spatial structure factor and the amplitude of indicators of DBH in them.
ClassN, pcsMode, cmHodges-Lehman EstimatorLower Confidence IntervalUpper Confidence Interval
1688.59.38.2910.35
2452313.411.7514.9
313multi *17.014.5519.65
* distribution having several modes.
Table 9. The relationship between DBH and the distance to the nearest trees in the directions of the cardinal points.
Table 9. The relationship between DBH and the distance to the nearest trees in the directions of the cardinal points.
Indicators *LnLsLeLwLav
Coefficient of correlation+0.65+0.69+0.45+0.67+0.72
* See Table 3.
Table 10. Characteristic of the best trees (2021).
Table 10. Characteristic of the best trees (2021).
RowTree NumberHalf-SibDBH, cmTree Height, mAverage Distance to the Nearest Trees, m
1420221.518.03.05
21201522.620.12.13
22161722.019.21.89
2824220.321.23.14
39171523.720.92.76
42151721.020.53.85
5091720.621.82.50
50131723.722.43.24
50161720.420.32.93
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Marozau, A.; Kotszan, U.; Nowakowska, J.; Berezovska, D.; Moser, W.K.; Hsiang, T.; Oszako, T. The Influence of the Provenance and Spatial Structure on the Growth of European Silver Fir (Abies alba Mill.) of Autochthonous Origin in a Forest Plantation in the Białowieża Forest. Forests 2022, 13, 831. https://0-doi-org.brum.beds.ac.uk/10.3390/f13060831

AMA Style

Marozau A, Kotszan U, Nowakowska J, Berezovska D, Moser WK, Hsiang T, Oszako T. The Influence of the Provenance and Spatial Structure on the Growth of European Silver Fir (Abies alba Mill.) of Autochthonous Origin in a Forest Plantation in the Białowieża Forest. Forests. 2022; 13(6):831. https://0-doi-org.brum.beds.ac.uk/10.3390/f13060831

Chicago/Turabian Style

Marozau, Aleh, Uladzimir Kotszan, Justyna Nowakowska, Daria Berezovska, Warren Keith Moser, Tom Hsiang, and Tomasz Oszako. 2022. "The Influence of the Provenance and Spatial Structure on the Growth of European Silver Fir (Abies alba Mill.) of Autochthonous Origin in a Forest Plantation in the Białowieża Forest" Forests 13, no. 6: 831. https://0-doi-org.brum.beds.ac.uk/10.3390/f13060831

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