Role of Silviculture in the Formation of Norway Spruce Forests along the Southern Edge of Their Range in the Central Russian Plain

: East European forests dominated by Norway spruce ( Picea abies (L.) H. Karst.) in the broad-leaved–coniferous zone should be considered as secondary communities formed under the inﬂuence of centuries-long activities (logging, ﬁres and planting) and extended outside their natural range. The study raises an issue—how stable is the current state of Norway spruce forests in the center of the Russian plain and what is the e ﬀ ect of silviculture on the forest cover of the large urban agglomeration—the Moscow Region? Current study is based on multidisciplinary research and consequently concerns the age dynamics of spruce plantation forests, the species and typological diversity of mature spruce forests and spatial pattern of spruce forests along the south edge of their range. The composition and structure of Norway spruce plantations have been studied for various age classes and compared with mature natural spruce forests and pine plantations on the basis of ﬁeld data. Remote sensing data and modeling approach were applied to estimate the spatial structure of spruce forests. It is found that mature plantations (over 80 years) become similar to natural forests in terms of structure and composition. The relationship between the distribution of spruce formations and the climatic and geomorphological conditions are conﬁrmed. The proportion of spruce and spruce–aspen / birch communities follows the pattern of zones—the transition from the coniferous and broad-leaved forest zone to the broad-leaved forest zone. Despite the signiﬁcant anthropogenic impact and the high proportion of plantations in the composition of Norway spruce forests (about 60–80%), their ﬂoristic and typological diversities correspond to such properties of zonal broad-leaved–coniferous communities. Over-matured plantations can provide valuable habitats for the re-establishment of native typological diversity. This makes it possible to use silviculture stands as an accelerated alternative to the natural recovery of disturbed habitats. nemoral boreal moss species— Atrichum undulatum (3%), Cirriphyllum piliferum (2%), Eurhynchium angustirete (8%), Plagiochila porelloides (4%), Plagiomnium undulatum (6%), P. a ﬃ ne (3%), Pleurozium schreberi (5%), Sciuro-hypnum curtum (3%), S. starkei (2%), Rhytidiadelphus triquetrus (4%), Hylocomium splendens (5%).

Despite considerable research experience in the region, for instance, studying the typological diversity of forest communities and their dynamics [40][41][42][43][44][45], there is no up-to-date information on the structure and composition of forests in the Moscow Region. The datum on biodiversity development during silviculture activities are also insufficient. Therefore, the south edge of the spruce forests natural range in the broad-leaved-spruce forests zone is important to study in terms of its possible changes and prospects of sustainable growth. Considering the level of anthropogenic impact, the Moscow Region has feasible potential as a case study area. The objective of this paper is to assess the input of silviculture in the formation of Norway spruce forests along the south edge of their range within the central Russian Plain (case study the Moscow Region). The multidisciplinary approach included the following tasks. We compared the process of spruce plantation forests maturation with natural spruce forests and with pine plantations. We than studied the current composition of mature spruce forests. Finally, we estimated the spatial structure and the pattern of the edge of spruce forests range within the case study area. Based on the results we proposed a set of general forest management recommendations considering spruce plantation and natural forests.

Study Area, History of Formation and Current Landscape
The Moscow Region is located in the central part of the East European (Russian) Plain-35°10′-40°15′ East, 54°12′-56°55′ North-and covers an area of 4.43 million hectares (less the area of Moscow) (Figure 1a). After expanding its administrative borders in 2012, Moscow became the 6th largest city in the world in terms of area, while the pressure on the region increased considerably due to rapid transport and housing development. The forests of study area which constitute the so-called green belt, are extremely important to maintain the environment and perform social functions taking into account active urban planning in the region, [46][47][48]. Currently, the forests in the Moscow Region have certain distinctive features:  considerable recreational pressure related to a high population density;

Study Area, History of Formation and Current Landscape
The Moscow Region is located in the central part of the East European (Russian) Plain-35 • 10 -40 • 15 East, 54 • 12 -56 • 55 North-and covers an area of 4.43 million hectares (less the area of Moscow) (Figure 1a). After expanding its administrative borders in 2012, Moscow became the 6th largest city in the world in terms of area, while the pressure on the region increased considerably due to rapid transport and housing development. The forests of study area which constitute the so-called green belt, are extremely important to maintain the environment and perform social functions taking into account active urban planning in the region, [46][47][48]. Currently, the forests in the Moscow Region have certain distinctive features: • considerable recreational pressure related to a high population density; • limitation on commercial cutting; in accordance with the Forest Code of the Russian Federation, forests of study area are protected (the status of protected forests means that clear cutting is prohibited except in cases of forest damage by disease, insects and fire, when there is a need for clear cutting); • poor forest sanitation conditions and an unstable environmental situation related to a lack of mandatory forest management actions.
A brief comparison of the climate and landscape features of the physiographic provinces of the Moscow Region is given below. The climate and terrain were compared using the WorldClim [51] and SRTM spatial databases [52]. The comparison of soils and underlying deposits was mainly grounded in the literature [49].
The Moscow Region has a temperate continental climate according to the map of climate zones [53]. The average annual air temperature is 2.7-3.8 • C and the precipitation is 479-644 mm (Table A1). In general, there is a latitudinal gradient of temperatures and precipitation. The winter temperature diagram indicates a warming impact of Moscow [54,55]. The elevated provinces in the north and northwest are characterized by a slight decline in summer temperatures compared to the general trend-which is most likely related to the impact of vertical zonal sequence.
The relief of the study area is gently sloping, with heights ranging from 90 to 320 m a.s.l.; the average height is 174 m a.s.l., while the average slope is 2.06 • (0 • -30.9 • ). Provinces #1 and #5 are related to a depression in the crystalline basement and, consequently, to active fluvioglacial flows in the Quaternary period. They have the lowest elevation. The regional subdivide between the rivers Volga and Oka passes through provinces # 2, #3 and #7. These provinces, located at higher elevations, have more rugged relief and steeper slopes.
The gradient change and mosaic of sandy and sandy loam substrates is a crucial factor shaping the edaphic conditions of the Moscow Region. The elevated northern and western provinces of the study area (#2 and #3) are dominated by moraine deposits which consist of sandy middle and heavy loam with sand lenses. In lowland, fluvioglacial deposits (provinces #1, #4 and #5) sands with loam interlayers prevail. In the western provinces (#3 and #4), nutrient-rich lacustrine and fluvioglacial deposits occur. Loess carbonate loam mantles prevail in the southern part of the area. Subdominant and azonal types of deposits in the study area are silty loam mantles, alluvial in river valleys, swamp, deluvial and aeolian deposits. Sands and sandy loams are poor in nutrients, but they are well-drained and permeable for air. Moraine deposits are nutrient rich; however, they have low permeability and lack drainage. The combination of nutrient content, substrate drainage, groundwater level and other conditions has given rise to a variety of soils, i.e., sod-podzolic, podzolic, swamp-podzolic and grey forest soils.

Vegetation and Phytogeographical Zoning
The main part of the Moscow Region is located in the broad-leaved coniferous forest zone in accordance with the phytogeographical zoning layout. The borderwith the broad-leaved forest zone is located in the southern part of the region [56][57][58]. Zone boundaries, here, are associated with a change in geomorphological conditions determined by the location of the region in relation to the borders of the Moscow glaciation: sandy-loamy plateaus with poor drainage in the west are replaced by well-drained eroded interfluve areas to the east.
The Moscow Region has different types of forests: broad-leaved coniferous forests and spruce subnemoral and boreal forests, as well as broad-leaved forests, the range of which depends on the landscape structure of the territory [59]. The upland forests are mainly mixed spruce and spruce-aspen/birch forests. Broad-leaved forests (oak and lime) cover the elevated surfaces with richer soil habitats, parts of the hill slopes and river terraces. Boreal spruce forest types are not widespread. Pine forests on the uplands are usually plantations.
Silviculture in the Moscow Region started in the first half of the 19th century [60]. It was difficult to promote silviculture and the area of plantations prior to 1917 was small-14,600 ha (1.7% of the total forest area) [61]. Silviculture was really scaled up only in the second half of the 20th century, arising from the need to tackle a shortage of wood after the war. In the late 1940s, pine plantations were widespread, but starting from the 1960s, there was a transition to spruce plantation, to which primarily ungulates, as well as diseases and insects, cause less damage [62,63]. The spruce forests in the region have recently suffered severely from outbreak of European spruce bark beetle (Ips typographus) [40], which is the reason for planting the plantations in the area which was made free after cutting down the trees. Overall, 4-8 thousand hectares of new plantations are created every year [31].
According to the Federal Agency for Forestry, forests including spruce currently account for 25% of the total forest area in the Moscow Region [31]. In terms of age, young spruce forests (younger than 40 years) prevail 39%, whereas the forests aged between 40 and 60 years, maturing forests (60-80 years old) and mature and overmatured forests (older than 80 years) cover a similar area, i.e., 23%, 22% and 17%, respectively. Forests older than 80 years are no more considered as plantation forest, therefore, they are not included into plantation forest statistics. The most part of spruce forests presented by plantations the in the region. The opinion of some experts is consistent with this data [39].

Methods
We design the following field work to answer the objectives: (1) compare the process spruce plantation forests maturation with natural spruce forests and with pine plantations, (2) estimate the structure and composition spruce forests in region and (3) classification of field relevés for (4) compiling a subsequent training set for mapping. The field data were collected using a standard methodology [64], in forest sample plots of 400-625 m 2 . A full species composition and the vertical structure on the site were recorded. The species coverage is given in percent. Layers are designated as follows: A-the tree layer (generative and senile trees), depending on the community structure, this was subdivided into two sublayers (A1 and A2); B-undergrowth (virginile trees-B1) and the shrub layer (virginile and generative shrubs, 1-10 m high-B2); C-herb and low bushes (including immature trees and bushes lower than 1 m); and D-the moss layer. Latin names of vascular plants are given according to [65,66]. Latin names of bryophytes are given according to [67][68][69].
For comparative study of the formation of spruce plantation forests, spruce natural forests and pine plantation forests, we designed the following approach. The main patterns of changes of the composition and structure of spruce plantations for the territory of the southwestern sector of the Moscow Region were studied in the course of age dynamics. As a comparison with spruce plantation forests, data were given on the formation of pine plantation forests, starting from the stage of middle-aged plantations (over 40 years old) as well as with mature natural spruce communities (over 80 years old). We compared ground layer cover and tree layer cover for these three formations.
The data on plantations are based on forest mensuration and on field and satellite data (monodominance, the age of communities and signs of planting in rows). However, judging from our experience, it is quite difficult to figure out, based on the visual features of the structure and composition of mature communities, whether the spruce communities are natural forests or silviculture stands. The information about old plantations was not often available from the forest management authorities. Therefore, the typological composition of mature communities was analyzed regardless of their origin.
To assess the advantages of forming the coniferous plantations of different species, we compared the composition of spruce (73 relevés) and pine plantations (83 relevés) of different age classes [14]. Relevés were made in similar habitat conditions-on uplands with sandy loam soils at sites with an area of (up to) 100 m 2 . As there were no young pine trees in the researched area, we were not able to monitor the entire cycle of their development (Table 1). To establish the natural origin of mature spruce forests, we used the old forest taxation data.
For assessing the age structure of the silviculture stands, we divided the coniferous species into groups based on their age, according to forest mensuration [70]. To provide a more detailed analysis of the beginning of reforestation, the first age class was divided into two stages. Therefore, seven age stages of forest development were established (Table 1). Total 83 86 19 The floristic similarity between the community groups was determined using the Sörensen Index [71]: where c is the number of common species in relevés A and B, whereas a and b are the numbers of species in relevés A and B, respectively.
To study the composition and the structure of mature spruce communities, we selected relevés with spruce dominating in the community. In spruce forests, we selected relevés where spruce accounted for 75% and more and in spruce-aspen/birch forests, we selected relevés where the proportions of spruce and aspen/birch were almost equal. Spruce communities are mentioned in 372 relevés, while spruce-aspen/birch communities are included in 237 relevés (the total number is 609). To study the typological structure and composition of spruce communities, we classified field relevés into 2 formations and each formation into 4 association groups. Classification of relevés was performed on a dominant basis, using the matrix principle [72]. Formations of the dominant forest species were identified. By formation, we mean the unit of classification, which is defined by the dominant species of trees. Each formation is represented by communities with different combinations of common dominating species in the understorey. Identification of syntaxa at the level of association groups was performed based on the prevailing ecological and morphological groups of plants in subordinate layers. In the studied area, each spruce and spruce-aspen/birch forest formation is represented by four association groups: • dwarf shrubs-small herb-green moss (DShG); • small herb (Sh); • small herb-broad herb (ShBh); • broad herb (Bh).
To determine significant differences between the species composition of various classification units (formations and association groups) and between the group means of description coordinates in the NMDS ordination space, the Duncan criterion was used (or Duncan s multiple comparison post-test) [73].
To interpret the ecology of the groups, indirect ordination methods were applied-nonmetric multidimensional scaling (NMDS ordination) in the R software environment [74], using square root transformation, Wisconsin double standardization and the Bray-Curtis distance [75].
During ordination and interpretation of the axes for each description, its Ellenberg's ecological indicator values were calculated [76,77]. Indicators such as light (L), nutrients (N), soil reaction (R) and moisture (M) were calculated using Juice 7.0 [78]. Thus, factor values calculated for each relevé were used to position vectors in the ordination space and to assess the differences between the community The syntaxa of spruce and spruce-aspen/birch forests are described using groups of indicator, dominant and constant species. The indicator species for association groups were identified using the IndVal method [79]. The species with IndVal > 25% were considered as indicator ones.
The spruce forests spatial structure was estimated by mapping the association groups. The association groups were mapped using the supervised classification approach in two stages: preparation of a training sample [80][81][82][83] and classification. At the first stage, the relevés were filtered and decimated for further use as training sample. The second stage included supervised classification, where the supervised data are environmental variables including a set of multispectral reflectance and coefficients of the digital elevation model. The classification was performed using linear discriminant analysis (LDA) [84]. We used Landsat Level-2 Science Products including surface reflectance bands 1-7 and the spectral indices NDVI, EVI, MSAVI, SAVI, NDMI and NBR [85]. The weather of the Moscow Region is very cloudy except in September and March and it is covered by multiple Landsat scenes. Therefore, we performed mosaicking of images for 5 periods: for Landsat 8: . We used SRTM DEM and 12 topographic coefficients: elevation, slope, aspect, RMS error, insolation, Laplacian and convexity (cross-sectional, longitudinal, plan, profile and maximum and minimum convexity) [86]. The total number of environmental variables was 84. To reduce the impact of intercorrelation between environmental variables, we performed principal component analysis. After removing principal components with eigenvalue <1, we obtained 14 principal components for modeling. The recognition quality of LDA was equal to 47% (including other non-forest land cover classes). The recognition quality for the 8 association groups discussed in current paper varies from 2.9% to 89.5% (Table A2).
We demonstrate the example of a site index plot (with an area of 51,000 hectares) in the southwestern part of the Moscow Region, in accordance with forest mensuration data, silviculture stands were digitized and mapped.

Formation of Spruce Plantations at Early Stages (1-40 Years)
In the areas disturbed by logging and preparing soil for planting, the diversity of biotopes determines the floristic richness of communities at the early stages of their development. In young spruce plantations (up to 7 years old), this is manifested within the high cover of the herb-dwarf shrub layer (up to 95%) and high species saturation of the relevés (Figure 2a,b). This stage (1-7 years) differs most in species composition from both natural spruce forests (Cs = 54%) and older plantations. The initial stages of community development are characterized by a large proportion of weed and meadow-edge species, such as Chamaenerion angustifolium, Hypericum maculatum, Leucanthemum vulgare, Potentilla erecta, Calamagrostis epigeios, etc.
The processes of community development accompanied by canopy closing (Figure 2a) and shading of subcanopy space lead to the gradual elimination of light-demanding species of disturbed habitats. Forest species already present there as well as those introduced from nearby communities (Ajuga reptans, Dryopteris carthusiana, Fragaria vesca, Rubus saxatilis, etc.) start spreading. At the age of 21-40 years, young plantations and some aspen/birch trees reach 14-20 m in height and the maximum cover of the tree layer and undergrowth is observed, while the cover of the shrub and herb-dwarf shrub layers is at a minimum (Figure 2a,b). At this stage, plantations with a sparse herb layer often develop.
Thus, closing the tree canopy radically transforms the postcutting communities. The canopy begins to shape the structure and composition of the herb-dwarf shrub and moss-lichen layers. Canopy closure changes the climate conditions to less deep freezing of the soil, later thawing, later snow cover melting and less warming of the soil and lower air layers during the warm period, which eventually shortens the growing season and weakens the vegetation. The composition and structure of plantation communities becomes gradually similar to the natural forests of the region with age. The impact of spruce leads to age-dependent differences in the structure and floristic composition of the studied plantations. This primarily arises from the morphological and physiological features of the spruce stand. Spruce litter is characterized by a low pH, which leads to the formation of peaty forest floor and coarse humus with an acidic reaction. Together, these features of spruce forests activate podzol formation processes. Thus, closing the tree canopy radically transforms the postcutting communities. The canopy begins to shape the structure and composition of the herb-dwarf shrub and moss-lichen layers. Canopy closure changes the climate conditions to less deep freezing of the soil, later thawing, later snow cover melting and less warming of the soil and lower air layers during the warm period, which eventually shortens the growing season and weakens the vegetation. The composition and structure of plantation communities becomes gradually similar to the natural forests of the region with age. The impact of spruce leads to age-dependent differences in the structure and floristic composition of the studied plantations. This primarily arises from the morphological and physiological features of the spruce stand. Spruce litter is characterized by a low pH, which leads to the formation of peaty forest floor and coarse humus with an acidic reaction. Together, these features of spruce forests activate podzol formation processes.
The formation and closing of plantation stand reduce the availability of light and soil moisture, leading to restructuring of the plant community. Analysis of the herb-dwarf shrub layer cover showed that it moderately depends on the cover of layers A and B (r = 0.5). A significant downsizing of the herbs and shrubs cover in spruce plantations at the age of 21-40 years is due to the both canopy closure ( Figure 3) and competition for resources. It is worthy of note that the species saturation of such plantations is markedly reduced-on average, 24 species per 100 m 2 ( Figure 2) and the species composition changes significantly compared with earlier stages. The formation and closing of plantation stand reduce the availability of light and soil moisture, leading to restructuring of the plant community. Analysis of the herb-dwarf shrub layer cover showed that it moderately depends on the cover of layers A and B (r = 0.5). A significant downsizing of the herbs and shrubs cover in spruce plantations at the age of 21-40 years is due to the both canopy closure ( Figure 3) and competition for resources. It is worthy of note that the species saturation of such plantations is markedly reduced-on average, 24 species per 100 m 2 ( Figure 2) and the species composition changes significantly compared with earlier stages. Further changes in spruce plantations are more related to the community structure than to the floristic composition-the similarity coefficient Cs does not change, starting from the age of 21-40 years. At the age of 41-80 years, the differentiation of the tree layer into sublayers occurs. Small-

Comparative Analysis of the Formation of Middle-Aged and Mature Spruce and Pine Plantations (over 40 Years Old)
Further changes in spruce plantations are more related to the community structure than to the floristic composition-the similarity coefficient Cs does not change, starting from the age of 21-40 years. At the age of 41-80 years, the differentiation of the tree layer into sublayers occurs. Small-leaved species (aspen and birch) gradually disappear from the community (if they were not removed mechanically at previous stages during logging and thinning) and spruce starts to dominate in the stand (Figure 4).

Comparative Analysis of the Formation of Middle-Aged and Mature Spruce and Pine Plantations (over 40 Years Old)
Further changes in spruce plantations are more related to the community structure than to the floristic composition-the similarity coefficient Cs does not change, starting from the age of 21-40 years. At the age of 41-80 years, the differentiation of the tree layer into sublayers occurs. Smallleaved species (aspen and birch) gradually disappear from the community (if they were not removed mechanically at previous stages during logging and thinning) and spruce starts to dominate in the stand (Figure 4). When the processes of forest stand fall begin with the formation of windfall gaps, spruce plantations older than 80 years become closer to natural forests in terms of their floristic (Cs = 77%) and structural indicators. The gaps are sometimes occupied by lime, which usually accompanies spruce in forests of natural origin. The cover of shrub and herb layers increases in the gaps.
In mature communities (aged 80 and older), the height of the A1 subcanopy layer ranges between 27 and 33 m. In spruce plantations, it mainly consists of spruce and birch, while aspen occurs more rarely and sometimes, there are broad-leaved species such as oak and lime, but their proportion in the canopy is less than 1% (4). The proportion of the A2 subcanopy layer in spruce plantations is small (the canopy density is about 15%-20%) and it is not present in a third of all relevés. Most often, When the processes of forest stand fall begin with the formation of windfall gaps, spruce plantations older than 80 years become closer to natural forests in terms of their floristic (Cs = 77%) and structural indicators. The gaps are sometimes occupied by lime, which usually accompanies spruce in forests of natural origin. The cover of shrub and herb layers increases in the gaps.
In mature communities (aged 80 and older), the height of the A1 subcanopy layer ranges between 27 and 33 m. In spruce plantations, it mainly consists of spruce and birch, while aspen occurs more rarely and sometimes, there are broad-leaved species such as oak and lime, but their proportion in the canopy is less than 1% (4). The proportion of the A2 subcanopy layer in spruce plantations is small (the canopy density is about 15-20%) and it is not present in a third of all relevés. Most often, it includes spruce trees of silviculture origin, which have low growth rate and birch trees. Regeneration of broad-leaved species is active on some sites.
Comparison of the organization of spruce and pine plantations of more than 40 years old with natural spruce communities revealed special patterns in their composition formation. The canopy density assessment showed that this indicator varies considerably in different communities. Pine plantations have the largest values (Figure 3a), which are related to pine morphological traits. In silviculture pine stands, the average canopy density is 50%, whereas in natural communities and spruce plantations, the average canopy density is 65%. The communities with regeneration of spruce under the canopy of pine have the highest canopy density.
In general, pine and spruce plantations are able to maintain the canopy species composition of the broad-leaved coniferous forests in the region. In natural forests, the proportion of aspen is slightly larger than in spruce plantations. In pine plantations, aspen did not occur. The proportions of broad-leaved species (lime and oak) are equal in the second subcanopy layer of spruce plantations and natural communities, whereas the proportion of lime is slightly smaller than in natural forests and the proportion of oak is somewhat larger in pine plantations.
The understorey and undergrowth (B layer) were well developed in all the analyzed communities. The cover was almost equal in all those community types. Regeneration of spruce was more active in pine plantations (8%) and less active in natural spruce forests (4%). Species such as Corylus avellana, Lonicera xylosteum, Viburnum opulus and Sorbus aucuparia are typical of the layer. A comparative analysis of species in the understorey and undergrowth showed that Daphne mezereum, Populus tremula and Viburnum opulus occur more often in natural broad-leaved coniferous forests. Pine plantations are grown in lighter and moister sites with more heliophilous and hydrophilous species such as Padus avium and Rubus idaeus, whereas Frangula alnus and Populus tremula almost do not occur there. The floristic composition of the layer B in spruce plantations is more similar to natural spruce forests than in pine plantations. The broad-leaved species such as Quercus robur and Tilia cordata occur more often in spruce plantations and their average cover is thicker in natural forests. This indicates that the older spruce plantations are closer to the natural zonal spruce forest communities.
With changes in age, under the canopy of plantations, especially pine, there is active regeneration of tree species which grow in the natural community. At the same time, spruce plantations maintain the floristic composition of the herb-dwarf shrub layer, which is closer to the initial composition; however, its regeneration is faster than in pine plantations.
The species of the herb-dwarf shrub layer (C layer) in spruce plantations are closer to natural forests in the region than the species of the layer in pine plantations. Some species occur more often in spruce plantations (Actaea spicata, Carex digitata, Circaea alpina, Impatiens noli-tangere, Luzula pilosa, etc.), while some of them are found only here (Adoxa moschatellina, Coccyganthe flos-cuculi, Daphne mezereum, Sanicula europaea, etc.)-the species included in the Red Data Book of the Moscow Region [87]. In spruce plantations, the number of species is larger, but the herb-dwarf shrub layer is not as thick as in natural spruce forests. It is worthy of note that the floristic similarity index of spruce plantations and both natural forests in the region and pine plantations is also 77%.
In the herb-dwarf shrub layer of pine forests, the species such as Geum rivale and Stellaria nemorum, which grow in moist and damp sites and Athyrium filix-femina, Dryopteris carthusiana and Ranunculus cassubicus and some other forest species that occur on sites with sufficient insolation and moist fertile soils, occur more often than in other communities.
Analysis of the number of vascular plant species showed that their largest number is in spruce plantations (Figure 2b). The number of species in the sampling areas ranges between 26 and 50 for spruce plantations, 15-42 for pine plantations and 22-44 species for natural communities.
Therefore, the research reveals that in the Moscow Region, the number of vascular plant species and bryophytes in the communities of coniferous plantations aged 80 and older is as large as the number of species in natural spruce forests. Mature and overmatured plantations (older than 80) are close to natural forests in terms of their structure and flora. In the case of timely care of planting (regular thinning, etc.), this is also possible at the age of 60 or 70.
In general, changes in the plant cover in spruce plantations occur faster than in natural secondary succession, as small-leaved species (aspen and birch) do not dominate for a long time.

Composition and Structure of Spruce and Spruce-Aspen/Birch Forests of Moscow Region
Our study demonstrated that it is hard to distinguish between natural or artificial origin of spruce communities when using visual patterns of community structure and composition. The data on old plantation activity often lacks in forest management administrations. In this respect, we analyzed species and typological composition regardless of mature community origin. This approach is substantiated in Section 3.1.

Typological Structure
As a result of the ecological and phytocoenotic classification of spruce and spruce-aspen/birch forest communities, we identified 8 association groups, based on the composition of the tree and subordinate layers ( Table 2). Representative images of association groups are provided in Figure 5.
To properly analyze the composition of the selected syntaxa at the level of formations and association groups, differences were identified for vegetation combinations of different layers. In particular, the differences between the two formations in the tree layer composition are significant according to the Mann-Whitney U-test (p < 0.01). Analysis of the composition of the tree layer between the 8 association groups using ANOVA revealed rather weak differences: F (7595) = 3.37, p = 0.001. This could account for the similarity of the tree species composition and the cover of species in the communities of spruce forests (groups 1-4) and spruce-aspen/birch forests (groups [5][6][7][8]. When comparing the vegetation composition of the 4 community groups (1-DShG forests, 2-Sh forests, 3-ShBh forests and 4-Bh forests), grouped based on the subordinate layers-C and D, excluding the tree layer, significant differences by the Duncan criterion were observed for all pairs. The value of the F-statistic was also high: F (3605) = 239.41, p < 0.01.
Furthermore, in general, there were significant differences in the ecological coordinate positions of the relevés when analyzing the total species composition of the 8 association groups: F (7601) = 35.85, p < 0.01. Significant differences in pairwise comparison were not confirmed (Duncan test) for the pairs of syntaxa: 1-2, 1-5, 1-7, 2-7, 3-8 and 5-6. This indicates the similarity of their composition (Table A5).
Thus, based on the composition characteristics of the two formations, the composition of the tree layer will be considered within all 8 association groups. The similarity of the vegetation composition of the subordinate layers of spruce and spruce-aspen/birch forests allows us to analyze the characteristics of the shrub (B2), herb-dwarf shrub (C) and moss (D) layers within the four community groups.
The organization of the spruce and spruce-aspen/birch forests is considered below according to a unified layout. In general, changes in the plant cover in spruce plantations occur faster than in natural secondary succession, as small-leaved species (aspen and birch) do not dominate for a long time.

Composition and Structure of Spruce and Spruce-Aspen/Birch Forests of Moscow Region
Our study demonstrated that it is hard to distinguish between natural or artificial origin of spruce communities when using visual patterns of community structure and composition. The data on old plantation activity often lacks in forest management administrations. In this respect, we analyzed species and typological composition regardless of mature community origin. This approach is substantiated in Section 3.1.

Typological Structure
As a result of the ecological and phytocoenotic classification of spruce and spruce-aspen/birch forest communities, we identified 8 association groups, based on the composition of the tree and subordinate layers ( Table 2). Representative images of association groups are provided in Figure 5.  To properly analyze the composition of the selected syntaxa at the level of formations and association groups, differences were identified for vegetation combinations of different layers. In

Tree Layer
The tree layer composition of spruce and spruce-aspen/birch forests differs in the ratio of coniferous and deciduous trees (Figure 6a). Spruce-aspen/birch forests include communities where spruce and birch are present in equal proportions (with a small mixture of aspen).
of the subordinate layers of spruce and spruce-aspen/birch forests allows us to analyze the characteristics of the shrub (B2), herb-dwarf shrub (C) and moss (D) layers within the four community groups.
The organization of the spruce and spruce-aspen/birch forests is considered below according to a unified layout.

Tree Layer
The tree layer composition of spruce and spruce-aspen/birch forests differs in the ratio of coniferous and deciduous trees (Figure 6a). Spruce-aspen/birch forests include communities where spruce and birch are present in equal proportions (with a small mixture of aspen). The spruce stands are defined by an uneven age (on average A1-95 and A2-65 years) and consist, as a rule, of 2 layers with an average cover of 65% and 32%, respectively (Table A2). In A1, the ratio of spruce and aspen/birch species (mainly birch) is stable at 4:1 (Figure 6a). However, in A2, the situation is different-in the group of Bh spruce forests, the proportion of spruce plummets, while that of birch, on the contrary, soars (Figure 6b). This fact points to the silvicultural origin of the Bh spruce forests without spruce regeneration. In the future, they will be replaced by broad-leaved species. This is confirmed by the nemoral composition of the ground layer.
The proportion of broad-leaved tree species is negligible. There is a trend that, from group 1 to group 4 in the second layer of the stand (A2), the abundance of broad-leaved tree species varies from The spruce stands are defined by an uneven age (on average A1-95 and A2-65 years) and consist, as a rule, of 2 layers with an average cover of 65% and 32%, respectively (Table A2). In A1, the ratio of spruce and aspen/birch species (mainly birch) is stable at 4:1 (Figure 6a). However, in A2, the situation is different-in the group of Bh spruce forests, the proportion of spruce plummets, while that of birch, on the contrary, soars (Figure 6b). This fact points to the silvicultural origin of the Bh spruce forests without spruce regeneration. In the future, they will be replaced by broad-leaved species. This is confirmed by the nemoral composition of the ground layer.
The proportion of broad-leaved tree species is negligible. There is a trend that, from group 1 to group 4 in the second layer of the stand (A2), the abundance of broad-leaved tree species varies from 3% to 15%. Thus, in DShG spruce forests, only single trees of Tilia cordata (5%) and Quercus robur (3%) occur in A2, but in Bh spruce forests, lime (15%) and oak (11%) occur regularly and lime is actively renewing (B-12%).
The communities of spruce-aspen/birch formations are younger, A1-75 years old and A2-about 50 years old, which corresponds to the successional status of this group of forests. The total cover of the tree layer is 70%; the cover of A1 is 60% and the cover of A2 is about 40% (Table A2). The proportions of spruce and small-leaved species are approximately equal; although in the A1 layer, the cover of small-leaved species (mainly birch) is higher than that of spruce, in A2, on the contrary, spruce significantly prevails over birch in terms of abundance. Typically, the broad-leaved species cover in the Bh group of associations (#4) is higher, especially in A2 and the abundance of spruce and birch is lower compared to other association groups (Figure 6b).
The total species richness of the community groups (including vascular plants and bryophytes of all layers of the community) was relatively high (341 species). That richness included: 27 trees, 21 shrubs, 5 small shrubs, 227 herbaceous plants and 61 bryophytes. The maximum number of bryophyte species in the Bh community group (#4) was 61 (Table A4).
The communities of the four groups considering only the species of subordinate layers are divided based on the species composition in the ordination space ( Figure 8). The leading differentiation factors characterized according to the Ellenberg scales are light, soil nutrients, soil reaction and moisture. Changes in soil reaction, nutrients and light are related to the first variation axis, whereas moisture is connected with the second axis. The correlation with soil reaction (r 2 = 0.51) and nutrients (r 2 = 0.35) is the most considerable (Table A6).
While the species composition of the community groups differs significantly, the indicators of species diversity are quite similar (Table A3). The number of species per unit area (species richness) of the herb-dwarf shrub layer varies slightly across community classes, ranging from 29 (boreal species) to 26 species (nemoral species). The number of bryophyte species compared in the community groups does not differ significantly. The minimum species diversity (201 species) due to the low number of species in the herb layer was characteristic of the Bh community group (#4). The cover of the herb-dwarf shrub layer was maximum here (Figure 7) and significantly differed from other community groups F(3.605) = 5.278, p < 0.01. The maximum species richness (232 species) was observed in the ShBh group (#3), which naturally reflects the mixed floristic composition of these communities.
The total species richness of the community groups (including vascular plants and bryophytes of all layers of the community) was relatively high (341 species). That richness included: 27 trees, 21 shrubs, 5 small shrubs, 227 herbaceous plants and 61 bryophytes. The maximum number of bryophyte species in the Bh community group (#4) was 61 (Table A4).
The communities of the four groups considering only the species of subordinate layers are divided based on the species composition in the ordination space ( Figure 8). The leading differentiation factors characterized according to the Ellenberg scales are light, soil nutrients, soil reaction and moisture. Changes in soil reaction, nutrients and light are related to the first variation axis, whereas moisture is connected with the second axis. The correlation with soil reaction (r 2 = 0.51) and nutrients (r 2 = 0.35) is the most considerable (Table A6). Analysis of the four forest community groups (DShG, Sh, ShBh and Bh) showed the special features of their variety in terms of species and typological diversity, with the transition from boreal species composition to the nemoral composition. For instance, within the sequence from the DShG group (#1) to the Bh group (#2) in the ecological spectrum, the proportion of boreal species decreased by a factor of 3.5, whereas the proportion of the nemoral species group tripled.

Spatial Distribution of Forests
According to our estimates, the total forest cover of the studied area is 51%, the proportion of spruce forests is 18% of the total area of the forest territory, the proportion of spruce-aspen/birch forest is 19%, the proportion of pine forest is 14%, the proportion of broad-leaved forest is about 8% and the proportion of broad-leaved-spruce forest is 1%. The maximum proportion belongs to the secondary forest communities with dominance of birch and aspen-it is more than 30% [80]. The distribution of formations and association groups is related to natural conditions-primarily to the soil-geomorphological and climatic conditions (Figure 9). spruce forests is 18% of the total area of the forest territory, the proportion of spruce-aspen/birch forest is 19%, the proportion of pine forest is 14%, the proportion of broad-leaved forest is about 8% and the proportion of broad-leaved-spruce forest is 1%. The maximum proportion belongs to the secondary forest communities with dominance of birch and aspen-it is more than 30% [80]. The distribution of formations and association groups is related to natural conditions-primarily to the soil-geomorphological and climatic conditions (Figure 9). The extent of spruce forests is the maximum in the elevated northern and western provinces (#2 and #3), which correlates with the thickness of the loamy terminal moraine deposits. The extent of spruce-aspen/birch forests is significant in northern lowland province #1, which is characterized by The extent of spruce forests is the maximum in the elevated northern and western provinces (#2 and #3), which correlates with the thickness of the loamy terminal moraine deposits. The extent of spruce-aspen/birch forests is significant in northern lowland province #1, which is characterized by lighter and better drained sand and loam-sand fluvioglacial deposits. The extent of spruce-aspen/birch forests in central provinces #4 and #5 dominates over the extent of spruce forests due to the low thickness of the loamy moraine deposits. DShG and Sh forests are typical for the northern lowland province #1, where the loam-sand fluvioglacial deposits as well as wetland habitats prevail. The extent of ShBh forests is greater in the south and southwest, which is connected with the high nutrient content of the underlying sediments (combined with drainage). The distribution of the Bh group is not entirely consistent with the distribution of the ShBh group, which can be explained by the contrast of hydromorphic conditions. In province #2, the soils have low water permeability because of the moraine loams and in province #5, it is due to the impermeability of the Jurassic clay bedrock. This distribution of boreal and nemoral type communities in the latitude direction reflects the zonal patterns of vegetation in the study area. A major extent of nemoral spruce communities in provinces #2 and #5 is connected with the "warming" effect of the megalopolis (Figure 9a,b). Yet, this assumption requires further research.

Spatiotemporal Dynamics of Spruce Plantations: An Example of the South-Western Part of the Moscow Region
As an example, we considered the occurrence of forest plantations in the forests on the moraine and moraine-fluvioglacial plains of the Smolenskaya province (#3) with the area of 51,500 hectares. Forests cover 56% of that area. The forests are mainly represented by natural forests (76% of the forested area). Plantations account for 22% of the forested area and the proportion of felled areas (without creating plantations) is 2% (Figure 10). In other parts of the region, the proportion of plantations is higher and reaches to 32% [88,89].

Spatiotemporal Dynamics of Spruce Plantations: An Example of the South-Western Part of the Moscow Region
As an example, we considered the occurrence of forest plantations in the forests on the moraine and moraine-fluvioglacial plains of the Smolenskaya province (#3) with the area of 51,500 hectares. Forests cover 56% of that area. The forests are mainly represented by natural forests (76% of the forested area). Plantations account for 22% of the forested area and the proportion of felled areas (without creating plantations) is 2% (Figure 10). In other parts of the region, the proportion of plantations is higher and reaches to 32% [88,89]. Spruce is the most widespread species used in silviculture (80% of all plantations), whereas pine and heterogeneous plantations account for a mere 6.5% and 5%, respectively. Besides, there are oak and larch plantations. Forest plantations often occur under the canopy of various species and their combinations: birch, birch and aspen, oak, spruce, spruce and birch, grey alder, aspen and lime. There is one case when forest plantations were created to replace old plantations aged 100.
According to forest inventory data, the oldest plantation in the area under study is dated 1891. Overall, according to the forest inventory data, prior to 1941 (before WWII), 281,500 hectares of forest were planted. Spruce is the most widespread species used in silviculture (80% of all plantations), whereas pine and heterogeneous plantations account for a mere 6.5% and 5%, respectively. Besides, there are oak and larch plantations. Forest plantations often occur under the canopy of various species and their combinations: birch, birch and aspen, oak, spruce, spruce and birch, grey alder, aspen and lime. There is one case when forest plantations were created to replace old plantations aged 100.
According to forest inventory data, the oldest plantation in the area under study is dated 1891. Overall, according to the forest inventory data, prior to 1941 (before WWII), 281,500 hectares of forest were planted.
Spruce and spruce-aspen/birch forests account for 57% in that area, while the proportion of plantations is only 20% of those forest communities. However, for instance, the proportion of plantations is 27% among spruce Sh communities and among spruce-aspen/birch ShBh communities, it is 24%. It is worthy of note that even in aspen/birch forests the proportion of plantations is large; for instance, in birch Bh forests, one-third of the forests are silviculture stands and the forest inventory data show this. These are probably the plantations which did not survive without the necessary care and birch dominated them. As spruce and spruce-aspen/birch communities have a large proportion of stands aged 80 and older, the proportion of silviculture stands should be much bigger.
In general, the distribution of plantations in terms of age in this model area is typical of the Moscow Region. Stages II and III (aged 21-40 and 41-60) and stage I (aged below 20) coniferous plantations cover the largest areas (64% and 28%, respectively). Stage V and VI plantations (aged over 80) occur most rarely-8% (Figure 10), which may be related to the fact that silviculture stands are excluded from the category of plantations aged over 80 in official sources (according to the forest inventory data).

Discussion
The regeneration of spruce communities with aspen and birch within the study area can proceed by either of two scenarios. According to the first scenario, the spruce population is formed during the primary overgrowth of felling or arable land with aspen and birch trees in the first 20-30 years or during the late regeneration of spruce that gradually enters the main canopy [90,91]. Spruce-aspen/birch forests, which include mainly spruce (Picea abies) and small-leaved tree species such as birch (Betula pubescens s. l., B. pendula) and aspen (Populus tremula), are regarded by many researchers [92][93][94][95] as a short-term secondary stage of spruce forests formed at logging sites as a result of both spontaneous succession and the development of spruce plantations. Aspen and birch trees preceding or accompanying the development of the spruce population are characterized by a high intensity of the production process [96,97]. Our study highlights another scenario: the mixed composition is formed as a result of active birch regeneration in spruce plantations with poor quality of tending. The frequent absence of tending in the Moscow Region in the last two decades has determined a large share of spruce-aspen/birch forests, the mixed composition which results from active birch regeneration in poorly tended spruce plantations [91]. According to experienced foresters [39] and old official data, the majority of spruce forests are plantations and covers 60-80% in the region.
The dynamic processes in spruce plantations studied in the example of the study area are generally consistent with the well-known "model of stand dynamics" developed for forests in Western North America [98]. After the disturbance, the community goes through the stages of initialization, canopy closing, biomass accumulation, maturation and vertical and then horizontal differentiation, reaching the stage of maturity. According to our study, the canopy closing changes the course of succession. The understorey of mature spruce plantation communities becomes similar to natural forests which comply to [99,100]. Our data showed that mature and elderly spruce plantation forests (over 80 years old) are similar to natural forests in terms of their structural and floristic characteristics, which is confirmed by the data from other studies [14,101,102].
The species diversity of vascular plants and bryophytes in communities of spruce plantations older than 80 years is comparable to that in natural spruce forests. The total species richness of the spruce communities (including species of vascular plants and bryophytes at all levels of communities) was relatively high-341 species. Out of these species, 27 were trees, 21-shrubs, 5-dwarf shrubs, 227-herbaceous plants and 61-bryophytes. The maximum number of bryophyte species in the Bh community group was 61. This data complies with well-known studies [14,[103][104][105].
The different data on species diversity in plantation forests are determined by various factors during their exploitation. The management regime should also be taken into consideration. For the forest plantations of Japan (where their proportion is 40% of area), it is shown that for uneven-aged stands, thinning and line cutting in plantations could improve the species richness or diversity of plant species [106,107], while after repeated rotation, the species richness decline [108]. Additionally, there is a recreational factor that influences the species richness. For example, the species richness in suburban forests of the Czech Republic is high and mostly depends on the human influenced landscape structure (distance to settlements, roads, etc.) [109]. For suburban Fagus sylvatica stands, it is demonstrated that the species richness, also the cover and height, of the herb and moss layer is lower in regularly visited forests [110]. Innovative forest management practices compared with traditional thinning result in higher stand productivity along with better structural and functional variables [111]. The data presented in current study allow us to challenge the hypothesis of low species diversity in silviculture stands [17,112]. However, the abovementioned studies conducted in the UK and Canada are justified for fast-growing plantations of introduced species that are cut down already at the age of 40. The authors note that a significant part of the floristic diversity, including vascular plants, in natural forests is due to the parcel and gap structure of forests. The importance of gap dynamics is also confirmed for temperate secondary forest ecosystems in Northeast China [113]. The age of young plantations (20-40 years old) is not sufficient for the internal structure of communities to develop. According to Ulanova [101], in the first 25 years of spruce plantation development after logging, the species diversity decreases, while there is a significantly high negative relationship (r = −0.90) between the number of species at the site and the time after logging. Our data did not show such a significant correlation in the study of plantations with the age up to 20 years (r = −0.48). We assume that the relatively high heterogeneity of the habitat and the appearance of weed and meadow-edge species [101] is due to the lack of proper care of plantations. In addition, the impact of forestry, namely, weeding and thinning of plantations, reduces the canopy closure of communities, giving the light-demanding species an opportunity to remain in the communities.
For the spruce and spruce-aspen/birch plantation forests, the well-known regularities of boreal and nemoral forest distribution with respect to latitude are confirmed, which reflects the zonal features of vegetation in the study region and is consistent with the studies of the spatial distribution of other types: broad-leaved-coniferous [114] and pine forests [115]. The majority of the communities dominated by spruce on the uplands in the center of the Russian plain in the zone of coniferous and broad-leaved forests should be regarded as secondary communities formed as a result of centuries-long forest management (logging, fires and forest planting), often alternating with agricultural activities. According to our estimates, the proportion of spruce and spruce-aspen/birch communities in the forest-covered area of the territory under study reaches about 37%. This indicator varies from 30% to 0.2% from northwest to southeast within the boundaries of the physiographic provinces. This regularity follows the change of zones-the transition from the coniferous and broad-leaved forest zone to the broad-leaved forest one.
The community composition was determined to form mainly based on growth conditions. Thus, in the study area, the DShG spruce forests are found on slopes of moraine rocks with middle sod-podzolic-middle loamy and sandy soils, with a groundwater level of 2-4 m [45]. Sh spruce forests are usually widespread on gentle slopes of low hills on fluvioglacial plains. According to the same authors, spruce forests with blueberry are more widespread in habitats with better drainage, in upper parts of moraine rocks with middle sod-podzolic and light loamy soils and with groundwater levels exceeding 4 m. ShBh spruce forests occur on arching and leveling parts of upland plateaus with middle podzolic light loamy and well-drained soils and are located close to uplands and moraine hills, mainly on slopes.
Spruce predominates on the elevated loamy final moraine deposits, whereas spruce-aspen/birch forests prevail on loam-sandy glaciofluvial deposits. In particular, DShG and Sh communities are the most characteristic of the Verkhnevolzhskaya (#1) province as, on the one hand, loam-sandy glaciofluvial deposits dominate there and on the other hand, the extent of bogging is noticeable.
Dispersion analysis of the distribution of various spruce forest and terrain types showed that the spruce forest types, which we identified, are determined by terrain and are statistically significant [116]. For instance, DShG spruce forests occur in the lowest and flattest areas with slight slopes. Sh spruce forests are spread in the lowest and concave terrain shapes. ShBh spruce forests are found on elevated, slightly concave surfaces with gentle slopes, whereas nemoral Bh spruce forests are located on the most elevated and arching terrain types with steep slopes.
Evidence suggests that mixed plantations (pine and spruce) are more stable and 40% more productive compared with monodominant ones [117,118]. As severe climate conditions are becoming more and more frequent in the studied region, according to the stress gradient hypothesis (SGH) [119], mixed plantations can be a good solution to maintain the ecological potential of forest cover. In particular, the reestablishment of pine plantations is appropriate, especially on river terraces on sandy loam soils, while mixed plantations of pine and spruce are promising on upland loam soils. It is worthy of note that the conditions in coniferous plantations may be favorable for the protected lichen species [17]. Our data and other studies [120] confirm the conclusion that it is necessary to involve ecosystem-oriented forest management afforestation actions when developing an optimum model of forest management.

Conclusions
We studied the spruce (Picea abies (L.) H. Karst.) forests within the Moscow Region. The Moscow region is feasible case study area for spruce forests on the south boundary of their range in the center of the Russian plain. Most of spruce forests are silviculture forests. The study area is also currently under weak forest management (frequent absence of tending in plantation forests). Our study is based on multidisciplinary research project and concerns the biodiversity and pattern of spruce forests, as well as age dynamics of spruce plantation forests.
The succession of spruce communities under poor tending begins with active birch regeneration and forming of mixed composition. From the forest management point of view, this is undesirable, as it reduces the forestry value of the main coniferous species; at the worst, spruce plantations collapse. However, the early successional species (birch and aspen) improve the edaphic properties of the soil, enriching it with biogenic elements and accumulating litter biomass, which is important to conserve biodiversity and normalize environmental processes.
During the age dynamics under the canopy of plantations, especially the pine ones, active restoration of tree species characteristic of natural forests was observed. Spruce plantations, in turn, maintain the floristic composition of the herb-dwarf shrub layer, which is the closest to the original and its recovery is faster than in pine plantations. The processes of vegetation dynamics in spruce plantation forests are faster than natural secondary succession as the stage when aspen and birch prevail is absent. Mature spruce and spruce-aspen/birch communities (the previous successional variant) are similar in structure and composition to the conditionally natural types of boreal and nemoral spruce forests. Importantly, the successional dynamics of forest plantations at the initial stages is mainly determined by the tending measures.
Despite significant anthropogenic impact and a high proportion of silviculture stands, these communities match the composition of zonal spruce communities of broad-leaved-coniferous forests in terms of floristic richness and typological diversity.
Analysis of the physiographical conditions confirmed the relationship between the distribution of spruce formations and the climatic gradient, as well as with the changes in geomorphological conditions associated with the position of the area relative to the limit of the Moscow glaciation. The relationship between nemoral (ShBh) forests and the drainage gradient of the underlying rocks may not be the sole key factor. In particular, the larger proportion of nemoral spruce communities in the area surrounding the Moscow urban agglomeration may be connected with the "warming" effect of the megalopolis, leading to the nemoralization of community composition. Yet, this hypothesis requires further research.
Based on the above analysis, we propose general recommendations for the management of sylvicultural spruce forests. Spruce plantation forests, especially in areas with a strong anthropogenic impact, with free seed access, have a potential for restoring the floristic diversity of natural forests, performing recreational and ecological functions. The widespread phenomenon of intensive decay of mature spruce forests in recent years, as well as the vulnerability of spruce plantation forests to dendrophagous insects, justifies the need to evaluate their combination with pine plantation forests. Spruce forests, as well as those combined with pine plantation forests, are promising for cultivation in the moraine plain on fresh, well-drained soils in the Moscow Region. Over-matured plantations can provide valuable habitats for the re-establishment of native biodiversity when seed sources are available. This enables silviculture stands to be used as a better alternative to natural overgrowth of disturbed habitats. In this case, plantations can also serve as buffers or ecological corridors between fragments of disturbed territories.   Bold-significant differences.  Note. mean ± standard error of mean. *-to correctly compare the species richness of coenofloras in different community types and to make the sample size equal, the number of relevés for groups 3 and 4 was reduced by random selection (up to 70 relevés). Table A5. Results of a pairwise comparison of the total species composition of the eight association groups according to Duncan criterion.