Homogenization of Temperate Mixed Deciduous Forests in Białowie ˙za Forest: Similar Communities Are Becoming More Similar

: Many studies show the signiﬁcant impact of direct and indirect human activity on the functioning of terrestrial ecosystems, including forests. The increase in the number of invasive species, changes caused by climate change, or eutrophication of habitats resulting from air pollution can irrevocably a ﬀ ect biodiversity, species composition, or species interactions. Many of these e ﬀ ects cannot be seen in commercial forests due to the signiﬁcant impact of direct human use of the forest and the high degree of transformation of forest ecosystems. In this work, we ask: how have forest communities changed over the past 70 years? What was the reason for these changes? To answer the above questions, we conducted research on repeated observations in the core area of the Białowie˙za National Park, which is characterized by one of the highest degrees of naturalness in Europe, where ecological processes have occurred without direct human intervention since the last glaciation. Studies have shown directional changes in species composition and biotic homogenization of three forest communities. Directional changes were found to be associated with both eutrophication of habitats as well as with changes in humidity and temperature. However, the observed changes in species composition were opposite to the hypotheses based on the observed global change. In contrast, changes in the species composition of the stand and the ability to shade and bu ﬀ er the temperature and humidity under the canopy caused changes in the species composition of forest communities. In the mixed deciduous forest, homogenization occurred along with the simultaneous change of species composition of forest communities. This was caused by an increase in fertility caused by increased nitrogen deposition and changes in environmental conditions prevailing under the canopy of trees, which, however, were caused by changes in the species composition of the stand.


Introduction
Global change is the main factor causing changes in plant communities [1,2]. Its most important components are (1) climate change, and especially increase in mean annual temperature, changes in the amount, character and seasonality of precipitation, and the resulting increase in frequency of droughts [3][4][5]; (2) an increase in the atmospheric deposition of nitrogen, which leads to eutrophication of habitats and changes in the proportion of nitrogen-demanding species to non-demanding species or even to a decrease in species diversity [6][7][8]; (3) biological invasions resulting in the displacement of native species and modification of local environmental conditions [9][10][11][12]; and (4) fragmentation and loss of habitats caused by land use change and overexploitation [13]. The effects of each of these factors, both separately and cumulative at the same time, are often difficult to predict and model [14], and can lead to many serious

Study Area
Our study was conducted in Białowieża Forest (BF), which is the transboundary forest complex (52.74 N; 23.87 E) with an area of approximately 1450 km 2 -600 km 2 belonging to Poland and 850 km 2 to Belarus. BF is one of the best preserved temperate forests in Europe, with continuous tree cover and a course of ecological processes not significantly modulated by man during the last 12 thousands years [34][35][36]. Substantial fragments of BF's ecosystems have preserved the natural character, expressed in diverse stand structure, species composition, and large amounts of deadwood [37,38], which is an effect of the late introduction of modern silviculture, including large scale wood extraction, which started there as late as the First World War [39]. The most common type of forest community in BF is oak-lime-hornbeam forest of the Tilio-Carpinetum type, which covers over 60% of the area [40]. This type of community develops on various types of eutrophic soils: typical eutrophic brown soils, typical clay-illuvial soils, and stagnogleyic clay-illuvial soils [41], with these English soil names based on [42]. In effect, species composition in different patches of vegetation is diverse, and consequently this plant community has been typically divided into three subcommunities [43], although it has also been divided up to seven subcommunities [37,40].
Since 1921, the best preserved part of BF has been strictly protected (i.e., no intervention approach), which recently has become the core area of the Białowieża National Park (BNP), where human activity is limited to tourism on a few restricted paths, for scientific exploration, and for maintenance of a few roads and touristic trails. Such a strict approach to nature conservation allows the study of long-term changes and dynamics of forest ecosystems without the need to take into account direct human interference in their course. The boundaries of the core area of BNP, in its original shape, were delimited by the rivers Hwoźna (northern border) and Narewka (western border), the state border between Poland and Belarus (eastern border), and the open habitats of the Białowieża village glade (southern border; Figure 1).
Forests 2020, 11, x FOR PEER REVIEW 3 of 14 eutrophic brown soils, typical clay-illuvial soils, and stagnogleyic clay-illuvial soils [41], with these English soil names based on [42]. In effect, species composition in different patches of vegetation is diverse, and consequently this plant community has been typically divided into three subcommunities [43], although it has also been divided up to seven subcommunities [37,40]. Since 1921, the best preserved part of BF has been strictly protected (i.e., no intervention approach), which recently has become the core area of the Białowieża National Park (BNP), where human activity is limited to tourism on a few restricted paths, for scientific exploration, and for maintenance of a few roads and touristic trails. Such a strict approach to nature conservation allows the study of long-term changes and dynamics of forest ecosystems without the need to take into account direct human interference in their course. The boundaries of the core area of BNP, in its original shape, were delimited by the rivers Hwoźna (northern border) and Narewka (western border), the state border between Poland and Belarus (eastern border), and the open habitats of the Białowieża village glade (southern border; Figure 1). Most of the core area (98%) is covered with forest communities, and 53% of forests are oak-limehornbeam forests, which we chose for our study. Matuszkiewicz, in the baseline vegetation survey carried out in 1949, distinguished three types of oak-lime-hornbeam forest subcommunities [43]: Querceto-Carpinetum caricetosum pilosae (QCCaricetosum) (n = 25), Querceto-Carpinetum typicum (QCTypicum) (n = 22), and Querceto-Carpinetum stachyetosum silvaticae (QCStachyetosum) (n = 29). These plant communities in 1949 clearly differed by species composition (Figure 2) and co-occurred in the following gradient of soil fertility: QCStachyetosum occupied the most fertile and humid habitats, whereas the poorer and drier habitats of this type were covered by QCTypicum and QCCaricetosum. The dominance of European hornbeam (Carpinus betulus L.) and Norway maple (Acer platanoides L.), with an admixture of common ash (Fraxinus excelsior L.) in the stands, distinguished QCStachyetosum forest from the other two types. The herb layer in this type of forest was inhomogeneous, often Most of the core area (98%) is covered with forest communities, and 53% of forests are oak-lime-hornbeam forests, which we chose for our study. Matuszkiewicz, in the baseline vegetation survey carried out in 1949, distinguished three types of oak-lime-hornbeam forest subcommunities [43]: Querceto-Carpinetum caricetosum pilosae (QCCaricetosum) (n = 25), Querceto-Carpinetum typicum (QCTypicum) (n = 22), and Querceto-Carpinetum stachyetosum silvaticae (QCStachyetosum) (n = 29). These plant communities in 1949 clearly differed by species composition ( Figure 2) and co-occurred in the following gradient of soil fertility: QCStachyetosum occupied the most fertile and humid habitats, whereas the poorer and drier habitats of this type were covered by QCTypicum and QCCaricetosum. The dominance of European hornbeam (Carpinus betulus L.) and Norway maple (Acer platanoides L.), with an admixture of common ash (Fraxinus excelsior L.) in the stands, distinguished QCStachyetosum forest from the other two types. The herb layer in this type of forest was inhomogeneous, often combining species typical to QCTypicum or ash-alder riparian forest Fraxineto-Alnetum, with large abundances of Asarum europaeum L. and Stellaria nemorum L. As indicator species of QCStachyetosum, Matuszkiewicz pointed out Impatiens noli-tangere L., Festuca gigantea L. Vill., and Carex remota L. QCTypicum was a community in which the hornbeam dominated in the tree stand, creating strong shading of the forest floor. The herb layer was homogenous, with high cover of Galium odoratum L., Galeobdolon luteum Huds, Stellaria holostea L., and Oxalis acetosella L. Querceto-Carpinetum caricetosum pilosae was a community very similar in structure to QCTypicum, but QCCaricetosum occupied the least moist habitats among the oak-lime-hornbeam forests. QCCaricetosum was a very homogeneous community with a small number of enclaves enabling the development of plants with various requirements for the habitat [43]. Despite the great importance of hornbeam in the stand layer, the presence of pedunculate oak (Quercus robur L.) and Norway maple as an admixture species of great importance made QCCaricetosum the most open type of oak-lime-hornbeam forest in Białowieża National Park, with the lowest crown cover. The understory can be considered grassy with a high proportion of Carex pilosa Scop. Matuszkiewicz mentioned Anemone nemorosa L., Daphne mezereum L., Carex digitata L., and Sorbus aucuparia L. as indicator species [43].  The indicator species analysis, based on the relationship between the species occurrence and abundance values from the set of sampled sites and the classification of the same sites into site groups, revealed a substantial decrease in the number of indicator species. Many plant species lost their diagnostic value-their number decreased from 25 to 7 in QCStachyetosum, from 9 to 2 in QCCaricetosum, and from 1 to 0 in QCTypicum. Only a few species gained the status of indicator species-three species in QCStachyetosum and one species in QCCaricetosum (Table 2).

Data Sampling
In 2018, we resurveyed 76 semi-permanent vegetation plots sensu Kapfer et al. [44] with an area of 100 m 2 (10 m × 10 m), representing three types of the studied oak-lime-hornbeam forest, sampled originally by the team of Professor Matuszkiewicz in 1949. The plots were located by Professor Matuszkiewicz arbitrarily on the territory of the Strict Reserve of BNP, "in the homogeneous vegetation patches, representative of the plant community" [43]. Both vegetation surveys (baseline and resurvey) were carried out with the use of the six-step Braun-Blanquet scale for estimating plant species cover [45]. Three forest layers were defined, following the original survey: (1) understorey (all herbaceous plants regardless of their height, and trees and shrubs <0.5 m high), (2) shrub layer (trees 0.5-6.0 m high, and all shrubs higher than 0.5 m), and (3) canopy layer (trees over 6 m high). However, due to the inaccurate nature of historical data on the shrub layer, we omitted this layer in our analysis. We relocated the plots using exact descriptions of their location given by Matuszkiewicz [43]. BF is divided into numbered square compartments (1066 × 1066 m), with borders maintained since the end of the 19th century. Descriptions presented the location of the research plots relative to the four cardinal directions, the corner of compartment (marked with concrete poles with the numbers of compartments), or the distance from the road or characteristic point (e.g., an old, well known oak). To maximize the accuracy of plot relocation and make it more objective, we used ArcGIS 10.3 software to determine the geographic coordinates of their central point, on the basis of their distance from the permanent field marks (usually corners of a forest compartment). To locate a plot in the field, we used the QField application for smartphones (ver. 1.1.0 Matterhorn), which is compatible with the ArcGIS program and in optimal conditions allows accuracy down to a few centimetres. Both series of observations-original (1949) and resurvey (2018)-were carried out in July and August of the respective years, so as to stay in the timeframe of a similar phenological period.

Data Analyses
Ellenberg's Indicator Values (cwmEIV) for light (EV_L), temperature (EV_T), moisture (EV_M), pH (EV_pH), nutrient availability (EV_N), and shade casting ability index (SCAI) (following Verheyen [46]) were calculated as the weighted means of the species composition (percentage cover) of the understorey and stand. The Shannon-Wiener Index was calculated by percentage cover. All statistical analyses were carried out in R 3.4.2 using Rstudio 1.1.383. The "ggplot2" [47] and "vegan" [48] packages were used to generate graphical results. The Mahalanobis distance was calculated using the mahalanobis function of R. Student's paired t-test was used for analysis differences between two series of observations for forest layer covers, community weighted means of Ellenberg's Indicator Values (cwmEIV), species richness (SR), Shannon-Wiener Index, and SCAI. Change in Mahalanobis distance was tested by Student's unpaired t-test. Principal component analysis (PCA) was prepared on the basis of species cover expressed by the Braun-Blanquet scale. PCA was selected to present data on the basis of the length of the PC1 axis gradient. Differences between frequencies (the number of plots where the species occurred) of understorey and stand species on plots were analysed with the use of the chi-squared test. To identify indicator species of the oak-lime-hornbeam forest types for both series of observations, we used the "indicspecies" [49] package and the "multipatt" function.

Results
There were 135 species of understorey vascular plants recorded at 76 plots in total-100 species in 1949 (the baseline survey) and 109 in 2018 (the resurvey). In contrast to the increase in total species richness, the mean species richness (SR) on all plots decreased by two species per plot (from 30.8 (±7.35) to 28.5 (±6.78); t = −2.1605, p = 0.033), but when it comes to particular type of forest community, the decrease in SR was significant only in QCStachyetosum (by 5.36 species; Table 1). The mean value of the Shannon-Wiener index on all plots decreased from 3.3 (±0.25) to 2.7 (±0.26) (t = −17.11, p < 0.001), and differences were similar among all types of studied forest. The frequency of species on all plots increased significantly for 15 species and decreased for 19 species (Table S1). Geophytes were among the species that showed a tendency to decline: Adoxa moschatellina L., Anemone nemorosa L., Asarum europaeum L., Hepatica nobilis L., and Dentaria bulbifera L. At the same time, graminoids were among those increasing in frequency: Brachypodium sylvaticum Huds., Calamagrostis arundinacea L., Dactylis glomerata L., Poa nemoralis L., and Carex remota L. (Table S1). This pattern was similar in all three types of studied forest community (Table S2).
The changes in species composition and frequencies resulted in a significant decrease in cwmEIV for temperature and an increase in cwmEIV for moisture and nutrients in all studied communities. The cwmEIV for pH decreased in QCStachyetosum and QCTypicum, whereas light cwmEIV decreased only in QCCaricetosum (Table 1).
Principal component analysis (PCA) revealed directional changes along the PC1 (Figure 2). Points representing both surveys were concentrated in clearly separated clouds on the PCA scatter plot. Simultaneously, distances between points representing resurveyed plots (year 2018) decreased significantly in comparison to their distribution in the original survey (1949), expressing their higher species composition similarity, both in total as well as inside each type of studied forest (mean Mahalanobis distance for 1949 was 0.259 and for 2018 was 0.012; t = 8.612, p < 0.001). The shift between the surveys was driven by pairs of the cwmEIV: pH and temperature vs. moisture and nutrients. Additionally, the cwmEIV for light separated QCCaricetosum from two other types of oak-linden-hornbeam forest along PC2 (Figure 2) in 1949, but not in 2018. The first two PCA axes explained 39% of variance (25% (PC1) and 14% (PC2)).
The indicator species analysis, based on the relationship between the species occurrence and abundance values from the set of sampled sites and the classification of the same sites into site groups, revealed a substantial decrease in the number of indicator species. Many plant species lost their diagnostic value-their number decreased from 25 to 7 in QCStachyetosum, from 9 to 2 in QCCaricetosum, and from 1 to 0 in QCTypicum. Only a few species gained the status of indicator species-three species in QCStachyetosum and one species in QCCaricetosum (Table 2). Table 2. List of indicator species for three types of oak-lime-hornbeam forest: Querceto-Carpinetum stachyetosum silvaticae, Querceto-Carpinetum caricetosum pilosae, and Querceto-Carpinetum typicum in two observation periods 1949 and 2018 (based on analysis of indicator species from the "indicspecies" R package). Stat-value-association value based on [49]; boldface-species that were indicators in 1949 and 2018.

Discussion
We revealed a strong homogenization of the studied types of mixed deciduous forests on two ecological levels. On the community level, their species composition became less heterogeneous inside each of them separately, and on the ecosystem level these three communities melted down into one quite homogenous set of species. This was well expressed on the PCA scatter (Figure 2), where the dispersion of the points representing the year 2018 was significantly lower than in the baseline survey. The homogenization among the studied communities was manifested by a complete overlapping in 2018 of the clouds of points representing the species composition of all communities, in contrast to 1949, when communities were clearly separated from each other on the scatter plot space and overlap only to a small extent.
There were several mechanisms leading to the homogenization of the studied oak-lime-hornbeam communities. One of them is change in the frequency of species that were indicator species in 1949 and lost their status or gained indicator species status in 2018. Some indicator species expanded and increased frequency in all three communities, which resulted in a decrease of their specificity. This was the case of Carex remota L. and Circea lutetiana L. Both were indicator species for QCStachyetosum in 1949, but due to their recent expansion and increased frequency in all studied communities (Table S2), they lost their status of indicator species (Table 2). Another example are species that were indicator species in the past, but reduced their frequency in the community for which they were specific, and this way lost their status as indicator species. An example of such a species was Carex digitata L., which in the past was the indicator species of QCCaricetosum, but its frequency decreased in this community in the period between the surveys, which resulted in the species being dropped from the list of indicator species in 2018. A special case was Galium odoratum Scop, which was the only indicator species for QCTypicum in 1949, where it occurred with the highest abundance; however, between the surveys, its frequency decreased to a level similar for all studied communities. In effect, from an analytical point of view, QCTypicum ceased to exist as a separate plant community due to the lack of indicator species. Four species that did not have a diagnostic value in 1949 became indicator species in 2018. This also took place via two paths. Lathyrus vernus Bernh. and Polygonatum multiflorum All., which were similarly abundant everywhere, reduced in frequency in two out of the three studied communities, becoming indicator species for the third community (QCStachyetosum and QCCaricetosum, respectively), where change was not significant. Another two species-Rumex obtusifolius L. and Euonymus europaea L.-became indicator species for QCStachyetosum by the significant increase of their frequency in this community, but not in the other two. The range of functional characteristics of these species is so wide that it is not possible to explain what the reasons were for the observed changes.
Another mechanism leading to homogenization were changes in the frequency of certain functional groups of species. In our case, it was caused by an increase in the frequency of graminoids in all types of studied communities: grasses (Brachypodium sylvaticum, Calamagrostis arundinacea, Dactylis glomerata, Deschampsia caespitosa, Milium effusum) and sedges (Carex remota, Carex sylvatica). An increase in the share of graminoids is often associated with two causes: increased grazing pressure or intensive eutrophication [50]. The impact of herbivores, in our opinion, was not the main reason for these changes due to their nature. With the high impact of herbivores, we would expect an increase in graminoid coverage but without major changes in attendance, because grazing by herbivores causes graminoids to compensate for their rapid growth [51,52]. The graminoids may also benefit from seed dispersal by large herbivores, as they account for 32% of plant species dispersed endozoochorically by the guild of ungulates in BF (Poaceae 19%, Cyperaceae 9%, and Juncaceae 4%). However, their share in the number of dispersed seeds is much lower [53], and the colonization success of endozoochorically dispersed seeds in the studied oak-lime-hornbeam forest communities is negligible [54], which again argues against their importance. Many publications suggest that the expansion of graminoids in forest ecosystems is caused by eutrophication, especially by increased nitrogen deposition [55][56][57]. Van der Wal et al. [50] pointed out that these two factors (herbivore pressure and nitrogen deposition) can propel each other on the basis of a positive feedback loop. The increase in fertility due to nitrogen deposition causes an increase in the share of grasses and sedges, which attract herbivores, and results in even greater fertilization of the soil by the depositing of faeces and the further increase in the share of graminoids [50]. The combined positive influence of these two factors on the graminoids was also possible in our study. However, the separation of their effects was not possible due to the fact that both factors in Białowieża Forest have shown increasing tendencies during the last half-century, and they provide feedback to each other [50]. The increasing densities of ungulates since the 1950s [58] are accompanied by an increase in nitrogen deposition since at least the 1980s, when direct measurement of pollutant deposition started [30]. These effects were also expressed in our study in the increase in the share of nitrogen-demanding species, which caused directional changes in the species composition of the studied forest communities. This observation is consistent with reports from all over Europe, where such directional changes in species composition are observed in forests [7,8,46]. Usually, eutrophication of habitats caused by increased nitrogen deposition is quoted as the main explanation for such changes [59].
Gilliam [60], in his work on the understorey plant response to an increase in nitrogen deposition, discussed many effects of eutrophication, including changes in the degree of herbivory, mycorrhizae, and the impact on plant diseases. He put forward the hypothesis that homogenization of forest communities is caused by an increase in nitrogen deposition in communities where nitrogen is a limiting factor. This causes a decrease in the spatial heterogeneity of nitrogen availability and drives homogenization of vegetation due to homogenization of soil fertility, which was also confirmed by Keith et al. [61]. In our study, homogenization of plant communities increased along with the increase in the original fertility of forest habitats, which was manifested by a significant increase in the share of species with a high Ellenberg's ecological indicator value for nutrients. The great importance of the increasing nitrogen deposition for explaining changes in the species composition of understorey plants and epiphytic lichens in Białowieża Forest was also reported by several other authors carrying out part of their research in Białowieża Forest [7,[62][63][64].
Another important process influencing the results of our study was the decrease in the frequency of early-spring species, especially geophytes Anemone nemorosa, Adoxa moschatellina, and Dentaria bulbifera. This change may be related to severe droughts during several years in the last decade, taking place from late spring and continuing till autumn in the studied area. Especially in April and May of the year of the resurvey (2018), very high temperatures and a very low amount of precipitation were recorded in the area of Białowieża Forest [65]. This brought serious consequences for geophytes, because their life cycle closes before the observation period (July/August) [40], and this year they were exposed to earlier senescence due to adverse conditions (drought and high temperatures). This probably happened with A. nemorosa, which in 1949 was one of the indicator species of QCCaricetosum, and thus its low detectability in the resurvey could have contributed to the homogenization of communities. However, this result could be biased by the ability of geophytes to survive unfavourable environmental conditions in the form of bulbs and rhizoms, among others [66], which may still be present in the community but were not recorded by us due to the shortening of their vegetation period. Shortening the vegetation period could be a consequence of severe environmental conditions (drought) or acceleration of phenology as it happens with flowering [67]. We were not able to verify this biasing effect because control resampling, planned for the following year (2019), was not possible due to even more severe weather conditions. The annual precipitation in 2019 was extremely low (459 mm), at the level of 73% of long-term mean annual precipitation (625 mm for the period 1985-2015 [31]) and extremely high temperatures (annual mean temperature 9.0 • C vs. 6.9 • C for the period 1950-2015 [31]).
Another factor that could also have had an impact on changes in the species composition of the studied communities involves changes in the stand, both in species frequency and average plot coverage. Over the past 70 years, many large-scale disturbances of the forest canopy took place in Białowieża Forest [34]. Dutch elm disease affected all native elm species (Ulmus glabra, Ulmus minor Mill., Ulmus laevis Pall.), ash dieback caused a significant reduction in the population of Fraxinus excelsior, and spruce bark beetle (Ips typographus L.) outbreaks significantly reduced the share of Norway spruce (P. abies), which resulted in a reduction in the frequency of these tree species on the studied plots. The declining species have been replaced mainly by small-leaved lime (T. cordata), a species with higher shade-casting ability than the replaced species [46]. As a result, the average shade-casting ability index of the forest canopy increased in each of the communities studied in contrast to the declined canopy cover. However, the reduction in the stand cover, oscillating around 10%, could be related to observer error, which, according to Morrison [68], is a common and acceptable value in this type of research. Similarly, the decrease in the cover of the understorey, although statistically significant (change by about 10%), cannot be explained by the changes in the environment, and is probably caused by the observer effect associated with its biased estimation. Usually, lower stand density, as well as an increase in nitrogen deposition and hence eutrophication of habitats, causes an increase in the cover of the understorey vegetation. In our case, the average coverage of the understorey decreased by approximately 10%, despite the parallel decrease in the coverage of shrub and canopy layers, which should result in better light access to the forest floor. The analysis of Ellenberg's ecological indicator values for light revealed the opposite-an increase in the shading of the forest floor, which had the highest impact on QCCaricetosum. The share of species with high Ellenberg's indicator value for light significantly decreased in this community (Table 1), which corresponds well with the increase in the shade-casting ability index of the canopy trees, which was the most open canopy in the past. This led to homogenization with the other studied types of communities ( Figure 2).
It is surprising that changes in species composition towards an increase in the share of cold-resistant and high moisture requirements species were contradictory with physical measurements of the climatic conditions. Both global [69] and local [31] climate characteristics indicate an increase in the mean annual temperature, as well as the more frequent occurrence of droughts and decreasing groundwater levels [31,32]. Thus, explanations for changes in species composition towards cooler and wetter conditions should be sought in the increase in SCAI. The increase in the mean shade-casting ability index of the stand is one of the factors increasing the ability of the forest to buffer temperatures under the canopy [70]. The lack of influence of climate warming on the species composition of epiphytic lichens, explained by the strong buffering effect of the tree canopy, was reported by Łubek et al. [62]. On the other hand, the increase in the mean shade-casting ability index clearly explains the decrease in the proportion of light-demanding species, such as Fragaria vesca L. and Plathantera bifolia Rich., which are usually associated with more open and relatively warmer forest communities [71].

Conclusions
The homogenization of forest communities is a poorly studied process, and the knowledge about its possible effects on biodiversity is limited. The number of publications on the homogenization of vegetation has increased with the increasing interest in global change effects [15]. According to the latest trends, most of the works related to the homogenization of plant communities concern the negative impact of alien species [20][21][22]. Only a small part of the homogenization works is associated with homogenization caused by the spread of native species [14,18,24].
The vegetation of the oak-lime-hornbeam forests of the Białowieża National Park has changed significantly during the 70 years that have passed since the baseline survey. The most visible change in the oak-lime-hornbeam forest communities was the homogenization of their species composition, and the melting of all three studied communities into one supercommunity. We consider the eutrophication of habitats and the expansion of native nitrophilic species as the reason for such changes. In addition, we conclude that changes in the stand could have contributed to changes in climatic conditions under the canopy of the stand, which results in changes in the species compositions of the understorey.
Supplementary Materials: The following are available online at http://www.mdpi.com/1999-4907/11/5/545/s1, Table S1: Changes in frequencies of species in all studied plots (number of occurrences, Nmax = 76). Table S2: Changes in frequencies (% of plots occupied) of plant species in three studied types of oak-lime-hornbeam forest in the Białowieża Forest, Eastern Poland. Table S3: Changes in frequencies of tree species in stand in all studied plots (number of occurrences, N = 76).

Conflicts of Interest:
The authors declare no conflicts of interest.