The Impact of Changes in Species Richness and Species Replacement on Patterns of Taxonomic Homogenization in the Carpathian Forest Ecosystems

: The Carpathians are Eastern Europe’s largest contiguous forest ecosystem and a hotspot of biodiversity. However, not much is known about changes in species diversity in these forests. Here we investigate the long-time changes in the diversity of plant communities in the Eastern Carpathian beech forests. We also inquire about the similarity of plant communities between managed and protected forests. We conduct our analyses using not only the broad-sense measure of changes in β diversity ( β sor ) but also the measures that reveal whether the changes result from spatial turnover ( β sim ) or nestedness ( β nes ). Our study demonstrated that the mean species richness did not signiﬁcantly change over time in both forest types. However, we found a signiﬁcant decrease of species diversity in protected forests and the same trend, but on a smaller scale, in managed forests. In both forest types the decrease of species diversity was mainly caused by spatial turnover, while nestedness-related changes were relatively small. However, the direction of changes in β nes differed in managed and unmanaged forests in such a way that it reduced the decrease of species diversity in managed forests and ampliﬁed the decrease of species diversity in unmanaged forests. We discuss our ﬁndings in terms of biotic homogenization.


Introduction
The changes in species composition and diversity are widespread in nature [1]. In most ecosystems, depending on population dynamics and the current set of multiple disturbances, species diversity increases or decreases. The average species composition rarely reaches a steady-state equilibrium [2]. Changes in species diversity and the direction of such changes (increase/decrease) are natural processes, and as such cannot be treated as positive or negative ones. For example, natural succession in forest plant communities often involves a decrease in species diversity during the establishment phase of colonization and an increase in species diversity during steady-state phases [3].
However, beside natural disturbances, which are important forces in maintaining species diversity, forest ecosystems are also subjected to human impacts, often resulting in biotic homogenization [4]. Human-induced environmental changes may decrease habitat suitability for a high number of species, at the same time increasing the suitability for a limited number of biotas better adapted to anthropogenic altered habitats. The process of continuing replacement of native specialists by generalists or nonindigenous, locally expanding species is rapidly increasing the similarity of species composition [5][6][7]. As the stability of ecosystems in the face of disturbances may be positively related to the diversity in these ecosystems [8], biotic homogenization is perceived as a global problem.
The beech forest is one of the major forest types in Europe [9,10], and thus, understanding ecological processes in this woodland is very important. Despite the fact that the Carpathians are Eastern Europe's largest contiguous forest ecosystem and a hotspot of biodiversity [11,12], not much is known about long-term changes in the diversity of plant communities in these beech forests. Studies have focused on species diversity changes observed in the beech forests of Western and Northern Europe (reviewed in [10]). The conclusions of these studies cannot be generalized and applied to Eastern Carpathian beech forests for several reasons. First, forests in Western and Northern Europe have been subjected to higher environment pollution than the Eastern Carpathian woodlands [13,14]. They have also been more intensively timber-harvested, which considerably affected the forest structure and biodiversity [15,16]. Also, the forests in Western and Northern Europe have been exploited for a considerably longer time (some of them since the Roman Empire) than the Eastern Carpathian forests [17]. This resulted in both severe clearance and forest plantation of the former and a lower level of human-induced disturbance in the latter woodlands. There are not many primeval beech forests left in Western Europe and habitat fragmentation is a common issue [14]. In contrast, the Eastern Carpathian ancient beech forests are a precious ecosystem which constitutes the East Carpathian Biosphere Reserve (1992, ECBR). They are home to many endangered species, which are threatened with extinction or are no longer present in Western Europe. Many of them are iconic and/or the International Union for Conservation of Nature (IUCN) Red Listed animals and plants [18]. Thus, understanding the long-term changes and the differences in the diversity of plant communities between managed and protected forests will help to maintain their uniqueness.
Because forests are dynamic mixtures of habitats over time and across landscapes, the spatio-temporal analyses of changes in species diversity provide a useful diagnostic tool. We use such an approach in the present paper. The availability of past datasets enabled us to analyze and describe long-time changes in the understory vegetation of two beech forests in the Carpathian Mountains in Poland using re-inventory data based on semi-permanent plots. Moreover, in one of the forests, a part of the plots have received no forest management after the first inventory as it came under protection in 1973 and now it constitutes the Bieszczady National Park. This created an opportunity to inquire whether long-term protection resulted in different patterns of changes in species diversity than those observed in managed forests. The aim of our study was to investigate the changes in the diversity of plant communities in the Eastern Carpathian beech forests over the periods of 30 and 50 years. We also inquired about the similarity of plant communities between managed and unmanaged (protected) forests.
According to the intermediate disturbance hypothesis, which argues that species diversity is maximum at intermediate disturbance levels [19], one may expect that species diversity will decrease with decreasing disturbances. Thus, in case of the protected area, which recovers from human-induced disturbances, a decrease in species diversity may occur as a natural process. Similar trends may be observed in managed forests due to the gradual introduction of more sustainable management during the last decades. We addressed the above issues by testing the following hypotheses: (i) diversity of plant communities decreased over a long time-period in the Eastern Carpathian beech forests and (ii) patterns of changes in diversity of plant communities were different in managed and unmanaged forests. We also demonstrated that the partition of data, taking into account those two distinct methods of forest use, revealed patterns of changes in diversity that were vague on the pooled analysis.

The Site of Study and Data Collection
In this study we used two datasets collected from the Polish part of the Eastern Carpathians. The first dataset was obtained from sampling plots located in the Western Bieszczady Mountains (WBM, 1560 km 2 , maximum altitude 1346 m above sea level (a.s.l.), 49˝03 1 N-49˝21 1 N, 22˝05 1 E-22˝47 1 E). The second dataset was collected on the adjacent, more northern part of the Eastern Carpathians and covered the majority of the area of the Sanocko-Turczanskie Mountains (S-TM, approx. 460 km 2 , maximum altitude of 672 m a.s.l., a part of the Natura 2000 protection program, Słonne Mountains PLH180013, PLB180003, 49˝25 1 N-49˝38 1 N, 22˝10 1 E-22˝30 1 E). The investigated montane region is built of flysch, which consists mainly of sandstone. The main types of soils are brown soils: Eutric Cambisols (brown soils leached and gleyed) and Dystric Cambisols (acid brown soils) [20,21]. The mean annual temperature and precipitation are 5.2˝C and 1100 mm in WBM and 7.3˝C and 814 mm in S-TM areas [22,23]. The area is characterized by a high forest cover (ca. 70%). The montane zone is mainly covered by the fertile Carpathian beech forest (dominated by Fagus sylvatica L.) with a fir (Abies alba Mill.) and sycamore (Acer pseudoplatanus L.) admixture (phytosociological association Dentario glandulosae-Fagetum (D.g.-F)) [20,24].
There are four phytosociological units within the analyzed fertile Carpathian beech forest, assigned along the gradient of nutrient availability, moisture, and soil pH: D.g.-F. festucetosum drymeiae on the driest, most acidic, and nutrient-poor soils; D.g.-F. typicum and D.g.-F. lunarietosum on intermediate soil conditions; and D.g.-F. allietosum ursinii on the moist and most fertile sites with a slightly acidic or slightly alkaline pH. They represent the whole range of habitat conditions present in the montane forest zone in this part of the Carpathians [20,25].
Over the past 200 years forest management practices in the study area underwent significant changes. Due to an intense and unsustainable exploitation combined with an increase in human population, the forest covered only 30%-40% of the area in the 19th and at the beginning of the 20th centuries [26]. Due to forced resettlements after World War II (WWII), the density of the local population drastically dropped. In those years, forest stands often had a uniform age structure and an average age of 40-60 years [25]. After the socio-economic changes of the 1990s, an emphasis was placed on sustainable forest management. This was mainly achieved by substituting the uniform shelterwood system (the main management system since WWII) with the irregular shelterwood system. Additionally, in the mid-20th century, the common practice of traditional non-timber forest use by the locals, such as firewood and litter raking or wood pasture, disappeared [27,28]. A part of the forests in the WBM area came under protection in 1973 and now it constitutes the Bieszczady National Park. The WBM and S-TM areas are a part of the Eastern Carpathian Biosphere Reserve (1992, ECBR). Based on the data obtained from the State Forestry administration and the authorities of the Bieszczady National Park, the average age of the majority of forest stands in the S-TM and WBM sampling areas exceeded 100 years in the 2000s.
We used two datasets obtained by the repeated sampling of semi-permanent plots. The first dataset (62 plots, 50 years between samplings) characterized the beech forest in the WBM area Phytosociological databases collected since the early 20th century are an important source of ecological information as they enable us to study long-term spatio-temporal changes in ecosystems. The weak point of such databases is that different researchers may determine the plots somewhat subjectively, which may bias the results and does not show the real differentiation of vegetation [29,30]. Additionally, the analyses may be burdened with an error resulting from the imprecise locating of past inventories. However, our datasets depicted full variability of investigated beech forest vegetation and we put an effort into the exact repetition of sampling. Thus, we assume observer's bias, which seems unavoidable, may play some role but does not account for a significant error.
For both study areas, sampling plots were located based on the data from archive phytosociological documents [20,25] and the sampling was repeated. In the case of the S-TM area, the plots were located with the use of a 1:45,000 scale unpublished map, drawn by Dzwonko during his study; for the WBM area, the plots were located with the use of coarse coordinates and a small-scale map published by Zarzycki [25]. To make the newly obtained datasets comparable with the old ones, the same sampling method (Braun-Blanquet approach, [31]) was used and the conditions such as the size of the sampling plot (in most cases between 200 and 500 m 2 in WBM and 200 or 400 m 2 in S-TM), topography (altitude above sea level, exposure, and slope), and time of the year were matched with those described in the original reports by Dzwonko [20] and Zarzycki [25]. It was difficult to locate the same plots precisely; however, our plots were not more than 50 m away from the plots of Dzwonko [20] and, at last, in the same forest vegetation patches investigated by Zarzycki [25]. Thus, we assume that differences between the two samples in the species composition reflected real changes that took place in the studied forests.
According to the species-area theory, the number of species increases along with the sampled area. Consequently, differences in species composition between the two plots of different size may result from differences in each plot area. However, we found no significant correlation between species richness and plot size in the S-TM dataset (mean number of species = 25.14 on plots 200 m 2 and 21.65 on plots 400 m 2 ) and even a weak negative correlation in the WBM dataset (r s = 0.28, p = 0.27, mean number of species = 28.65 on plots 200-300 m 2 and 27.08 on plots 400-500 m 2 ). It demonstrates that larger plots do not consist of higher numbers of species. Thus, we made the assumption that the potential methodological bias resulting from using plots of different size did not affect the results of our study.
The sampling plots from the 1950s and 1970s were located in the area under active forest management. Due to the conversion of the best-preserved forests of WBM into a protected area in 1973, 32 of 62 plots from the 1950s were included in the protected forests, from here on termed 'unmanaged' forests. As a consequence the dataset from the WBM area was divided into managed (30 plots) and unmanaged forests (32 plots). Our analyses concerned species composition (presence-absence data) of the herbaceous layer (including ferns and flowering plants and excluding tree seedling and bryophytes), which is the most disturbance-sensitive and diverse plant community in forests [32].

Data Analysis
In order to quantify biotic homogenization, changes in the taxonomic similarity were analyzed with pairwise dissimilarity indices between the two sampling periods for each dataset (WBM and S-TM forests). We used three β diversity indices: (1) Sorensen's dissimilarity index (β sor ), which measures dissimilarity in species composition between sites due to both spatial turnover and differences in species richness [33], (2) Simpson's dissimilarity index (β sim ), which expresses spatial turnover without the influence of richness gradients, (3) nestedness-resultant dissimilarity index (β nes ), which expresses compositional differences in richness caused by nestedness.
Using β sim and β nes allowed us to demonstrate not only general changes in the diversity of plant communities but also more subtle patterns of these changes [33]. A difference between a spatial turnover and nestedness lies in the fact that turnover implies the replacement of some species by others, while nestedness occurs when species of one site are the subset of the species of another (richer) site. Thus, nestedness mostly implies a species extinction in the species-poor site without the replacement of the lost species by other ones [33]. As both processes (species replacement and species extinction) may lead to the decrease of the diversity of plant communities and the increase of homogenization, it seems reasonable to analyze them separately.
Sorensen's dissimilarity index [34] was calculated according to the equation: For this and further equations, a denotes the number of species common for the two sampling plots, b indicates the number of species only occurring in the first sampling plot of the pair, and c denotes the number of species occurring only in the second sampling plot of the pair. Simpson's dissimilarity index was calculated using the equation [35,36]: Nestedness-resultant dissimilarity index was calculated using the equation [33]: max pb, cq´min pb, cq 2a`min pb, cq`max pb, cqˆa a`min pb, cq We analyzed the data using two different approaches: (i) among-plot dissimilarity between the two surveys and (ii) within-plot compositional shifts over time, following a method adapted from Baeten et al. [37].
In the first approach, our analyses of biotic homogenization were based on the most frequently used comparison of mean values of pairwise dissimilarity indices between temporal censuses. This method is largely equivalent to the measure of β diversity proposed by Anderson et al. [38,39]. To do it we calculated average interplot dissimilarities (for each plot the mean of the pairwise dissimilarities against all other plots in terms of β sor , β sim , and β nes ). The changes in plot taxonomic similarity were tested by comparing the mean values of dissimilarities (β sor , β sim , and β nes ) of each plot between past and present inventory using a paired t-test with a bootstrap procedure (n = 10,000).
In the second approach we used the NMDS (nonmetric multidimensional scaling) on the full data matrix (i.e., the data of the past and present inventory) for each forest to demonstrate changes in all β indices on the plot level. For each studied forest (S-TM, WBM, WBM managed, and WBM unmanaged) we created three dissimilarity matrices based on the three dissimilarity measures (β sor , β sim , and β nes ). These matrices were used for the ordination of plots. Next, vectors indicating mean shifts between sampling periods were created. This was achieved by running the NMDS analysis until a stable result was obtained. This stable result was then used as the final result of the analysis. Based on the analysis of the stress value against the number of dimensions (scree test) and the clarity of the final configuration of the NMDS, we used two axes in the case of the β sor , β sim matrix and one axis in the case of the β nes [40]. The significance of directional shifts between the pairs of sampling plots from two sampling periods was assessed using a permutational multivariate analysis of variance (PERMANOVA) [41].

Results
Altogether 180 plant species of the herbaceous layer were found in the study area, including past (1950s and 1970s) and present (2000s) inventories. Of these species, 113 occurred both in the past and at present. Forty-seven species recorded in the past were lost and 20 new species were gained in the 2000s. The highest loss in the number of species occurred in the WBM forest ( Table 1). The largest difference between the number of species lost and gained occurred in unmanaged WBM forests (29 species) and the smallest in the S-TM forest (three species). Approximately 16% of the total number of species disappeared and 19% was gained (Table 1). Mean sampling plot species richness in the investigated area was 25 species, both in the past (1950s-1970s) and at present (2000s). The highest number of species was recorded in sampling plots of managed WBM forests (on average 28 species in the 1950s and 27 species in the 2000s). We did not find any significant changes in the mean species richness in forests in the S-TM and WBM areas ( Figure 1). However, the decrease in species richness in unmanaged WBM forests was near-significant ( Figure 1).  Table 1). The largest difference between the number of species lost and gained occurred in unmanaged WBM forests (29 species) and the smallest in the S-TM forest (three species). Approximately 16% of the total number of species disappeared and 19% was gained (Table 1). Mean sampling plot species richness in the investigated area was 25 species, both in the past (1950s-1970s) and at present (2000s). The highest number of species was recorded in sampling plots of managed WBM forests (on average 28 species in the 1950s and 27 species in the 2000s). We did not find any significant changes in the mean species richness in forests in the S-TM and WBM areas ( Figure 1). However, the decrease in species richness in unmanaged WBM forests was near-significant ( Figure 1). A comparison of the past and present species richness in the S-TM and WBM forests suggests that species richness decreased in species-rich plots and increased in species-poor plots in WBM unmanaged forests (Figure 2a). The species richness gradient decreased in unmanaged WBM forests and slightly increased in both managed forests (Figure 2). A comparison of the past and present species richness in the S-TM and WBM forests suggests that species richness decreased in species-rich plots and increased in species-poor plots in WBM unmanaged forests (Figure 2a). The species richness gradient decreased in unmanaged WBM forests and slightly increased in both managed forests (Figure 2).  According to Sorensen's dissimilarity index (βsor) the highest species diversity decrease occurred in unmanaged WBM forests, while in managed forests it was either small (S-TM forests) or insignificant (managed WBM forests) ( Table 2). Most of the changes in the broad-sense species diversity (βsor) were explained by βsim, while βnes accounted for relatively small changes in all studied forests. Our results showed that βsim decreased significantly in all cases but the direction of changes in βnes differed between forest types. In both managed forests βnes increased and in the unmanaged forest it decreased ( Table 2). Due to the additive effect of decreasing βsim and βnes in the unmanaged WBM forest, the overall species diversity decrease (βsor) was the highest in this type of the forest. The combined analysis of all WBM According to Sorensen's dissimilarity index (β sor ) the highest species diversity decrease occurred in unmanaged WBM forests, while in managed forests it was either small (S-TM forests) or insignificant (managed WBM forests) ( Table 2). Most of the changes in the broad-sense species diversity (β sor ) were explained by β sim , while β nes accounted for relatively small changes in all studied forests. Our results showed that β sim decreased significantly in all cases but the direction of changes in β nes differed between forest types. In both managed forests β nes increased and in the unmanaged forest it decreased (Table 2). Due to the additive effect of decreasing β sim and β nes in the unmanaged WBM forest, the overall species diversity decrease (β sor ) was the highest in this type of the forest. The combined analysis of all WBM data (from managed and unmanaged forests) blurred the differences in changes in β nes (Table 2). This sheds some light on the importance of analyzing differently managed forests separately. Table 2. The mean values of species diversity for the first and second inventory and the changes in species composition dissimilarity in the investigated part of the Sanocko-Turczanskie Mountains and the Western Bieszczady Mountains between the 1950s/1970s and 2000s in managed and unmanaged forests. S-TM-managed forest: n = 66; WBM-whole forest: n = 62; WBM-managed forest: n = 30; WBM-unmanaged forest: n = 32 (paired t-test with bootstrap procedure (n = 10,000)). The NMDS analysis of changes in the species composition on the sampling plot level, based on β sor and β sim , indicated marked plot-level compositional shifts between the past (1950s and 1970s) and present (2000s) plant communities (Table 3, Supplementary Figures S1 and S2). Changes in dissimilarity based on β sim and β nes of the pairs of sampling plots from the two sampling periods differed in managed and unmanaged forests. The compositional shifts on the sampling plot levels based on β sim were lower in unmanaged forests than in managed ones. On the other hand, compositional shifts on the sampling plot levels based on β nes were statistically significant only in the unmanaged forest.

Discussion
The results of our study revealed that the mean species richness of the herb layer did not significantly change over time both in managed and unmanaged Eastern Carpathian beech forests. However, we observed a near-significant decrease of species richness in the latter woodlands, which indicates that the ecological processes in both forest types may slightly differ.
Despite no significant changes in species richness, our analyses of the herb layer vegetation structure revealed the significant decrease of species diversity in the unmanaged forest. The same trend, but on a smaller scale, was observed in managed forests. This finding confirms our hypothesis assuming that the diversity of plant communities decreased over a long time-period in the Eastern Carpathian beech forests. Whether or not the decrease in species diversity was a result of decreased human-induced disturbances is disputable because some other factors, for example forest aging, may have similarly affected species diversity. However, our study demonstrated that the highest decrease in species diversity occurred in the least human-disturbed protected forests. Thus, it seems plausible that the pronounced decrease of species diversity in Bieszczady National Park is a natural process of recovery from the past human-induced disturbances. This finding is consistent with results of surveys of German submontane beech forests at 10-year intervals [46].
Our study also inquired about the patterns of changes in the diversity of plant communities in managed and unmanaged forests. To date, analyses on changes in plant diversity have indicated that homogenization is the most frequent process occurring in the European forests [47][48][49]. In the majority of the studies addressing this issue, classical, 'broad-sense' Sorensen's beta diversity or Jaccard's dissimilarity indices were used. However, as demonstrated by Baeten et al. [37], if the more specific components of beta diversity, such as beta turnover and beta nestedness, are not discriminated, the actual patterns of changes in the species composition may be concealed. The results of our study revealed that both in managed and unmanaged forests the decrease of species diversity was mainly caused by non-random replacement of some species by others (spatial turnover), while nestedness-related changes were relatively small. However, the direction of changes in β nes differed in managed and unmanaged forests in such a way that it balanced β sim in the former but acted additively with β sim in the latter. As a result, nestedness-related changes in species diversity reduced the decrease of species diversity in managed forests and amplified the decrease of species diversity in unmanaged forests. Thus, our hypothesis assuming different patterns of changes in the diversity of plant communities in managed and unmanaged forests was confirmed by the results of our study.
The analysis of compositional shifts on the sampling plot level indicated that the main cause of homogenization in managed forests was the decrease in spatial turnover. In contrast, the changes in beta nestedness were not uniform-in some plots species richness increased over time, but in other ones it decreased. As a result of the above-mentioned processes, managed forests underwent homogenization and differentiation of the plant community simultaneously. This is in accordance with the results of an earlier study of the same forests [50]. In contrast, homogenization in unmanaged forests was caused not only by the decrease of turnover but also by the nestedness-related decrease in species richness. As a consequence, the species composition of the present sampling plots in unmanaged forests was a subset of the species composition of these plots sampled in the past. A similar pattern of homogenization was reported in a forest reserve in the UK [37].
There are various environmental drivers of change in plant community composition. In the Polish part of the Eastern Carpathians, these drivers were most likely related to the decrease in the intensity of forest use and forest aging [50][51][52], which resulted in a lower level of disturbances as well as an increase in the density of the canopy and the volume of the forest stand [28,53,54]. These factors most likely contributed towards higher habitat stability, which led to the decrease in species diversity.
In managed aging beech forests the accumulating litter hinders plant development and, combined with wood extraction by forestry, causes soil acidification. This especially affects the nutrient-poor beech forest communities [55,56], which naturally consist of a small number of species [20,25]. The above process may have caused a decrease in species richness in species-poor areas due to the extinction of the species demanding a higher soil pH and fertility (e.g., Corydalis cava Schweigg. and Körte, Paris quadrifolia L., Pulmonaria obscura Dumort.) [50,57]. At the same time species richness may have increased in species-rich areas, for example, due to the higher amount of nutrients available in more diverse stands. As presented by Durak and Holeksa [50], such conditions may favor the expansion of species with a wide soil pH range, and a strong ability to spread laterally and induce vegetative regeneration or regeneration by seeds/spores dispersed on long distans (e.g., Dryoperis dilatata (Hoffm.) A. Gray, Rubus idaeus L.). These processes were likely to trigger an increase in the species richness differences between plots, which in turn may have resulted in a decrease in nestedness at the community level.
The cessation of non-timber and forestry practices in the protected forests caused a significant decrease of disturbances, allowing the recovery of forest stands. The resulting increase in canopy closure contributed to a reduction in light availability at the forest floor [57]. It is also likely that the amount of coarse woody debris (considered as significant carbon and nutrient pools) increased, which may have improved the fertility of habitats [58][59][60][61]. According to the succession theory, the increase in the density of the canopy layer and the fertility of the forest habitats is concomitant with an increase in competition. Consequently, light-demanding species are likely to be replaced by the shade-tolerant ones (e.g., Epilobium montanum L. and Urtica dioica L. vs. Galeobdolon luteum Huds. and Lysimachoa nemorum L. [57]). Our findings also suggest that, due to the recovery of forest habitats, the plant species richness increased in the species-poor areas and it decreased in the species-rich forests. As a result, the diversity decreased.
According to Peet and Christensen [3], the diversity and species richness should increase in steady-state phases of forest stand development. Thus, the decline in species diversity and richness suggests that, despite 30 years of strict protection, Bieszczady National Park still represents a transitional phase of succession. In conclusion, our results indicated that plant diversity changes in time appeared to be related to forestry practices. In managed forests, plant diversity changes led to homogenization resulting from the turnover. However, it is also important to note that the process of homogenization occurred in parallel to the differentiation. In unmanaged protected forests, homogenization which occurred during forest habitat recovery was driven by both turnover and nestedness. Therefore, to fully disentangle and assess the impact of turnover and nestedness on vegetation diversity, future analyses should include the characteristics of the species-facilitating β turnover and nestedness changes in the managed as well as the protected forests.