Effect of Soil Diversity on Forest Plant Species Abundance: A Case Study from Central-European Highlands

: Plant distribution is most closely associated with the abiotic environment. The abiotic environment affects plant species’ abundancy unevenly. The asymmetry is further deviated by human interventions. Contrarily, soil properties preserve environmental inﬂuences from the anthropogenic perturbations. The study examined the supra-regional similarities of soil effects on plant species’ abundance in temperate forests to determine: (i) spatial relationships between soil property and forest-plant diversity among geographical regions; (ii) whether the spatial dependencies among compared forest-diversity components are inﬂuenced by natural forest representation. The spatial dependence was assessed using geographically weighted regression (GWR) of soil properties and plant species abundance from forest stands among 91 biogeographical regions in the Czech Republic (Central Europe). Regional soil properties and plant species abundance were acquired from 7550 national forest inventory plots positioned in a 4 × 4 km grid. The effect of natural forests was assessed using linear regression between the sums of squared GWR residues and protected forest distribution in the regions. Total diversity of forest plants is signiﬁcantly dependent on soil-group representation. The soil-group effect is more signiﬁcant than that of bedrock bodies, most of all in biogeographical regions with protected forest representation >50%. Effects of soil chemical properties were not affected by protected forest distribution. Spatial dependency analysis separated biogeographical regions of optimal forest plant diversity from those where inadequate forest-ecosystem diversity should be increased alongside soil diversity. more of saturation


Introduction
Distribution of plants on Earth's surface is most closely associated with abiotic environment diversity (geodiversity) and most significantly conditions the diversity of other living forms. The relationship between terrestrial plants and geodiversity is not uniform but rather spatially differentiated along climatic transitions, soil properties and the intensity of human interventions [1]. While climate defines supra-regional differences in the relationships between geodiversity and plants, soil properties define local differences in ecosystem relationships [2]. Soil properties preserve geodiversity effects on plants, even though human interventions have disrupted ecosystem relationships. Nevertheless, positive soil diversity effects on the relationship tightness among general ecosystem-diversity components affected by human activities at a supra-regional level are not taken into account [3].
The total ecosystem diversity includes geo-, bio-and functional diversity. While biodiversity is formed by species abundance of life forms, functional diversity is the interconnection complexity of metabolic transfers among all ecosystem components [4]. Biodiversity varies along variously favourable site conditions. In ecosystems that rely availability [21]. Supporting nutrient availability increases both biodiversity and ecosystem resilience [18]. Plants support nutrient accessibility by promoting their rhizosphere functions, which increases biodiversity by creating microbial niches, but also increases soil component diversity of geodiversity (pedodiversity).
In response, however, uneven local effects on pedodiversity development differentiate plant dependence on the abiotic environment proportionally to the extent of bedrock transformations during the Quaternary glaciation [22].
Uneven plant dependence on geodiversity is conditioned by geographical location and the extent of anthropogenic loading. Human activities affect all ecosystem components. Although local biodiversity is significantly declining due to human activities, species richness of biomes has declined little due to invasive-species introduction and transitioncommunity formation in a fragmented landscape [5]. Biodiversity is globally used for evaluation of ecosystem health status since it controls key ecological processes [23][24][25]. However, replacing diverse transitions with ecosystems simplified by exploitation, crop cultivation or pollution has reduced the original diversity as a result of species extinction [26,27]. The effects of biodiversity loss on individual ecosystem characteristics differ together with varying nutrient availability [16]. The effect of nutrient availability on ecosystem diversity within the cultural landscape was an initial prerequisite of the study aim to assess the dependence of semi-natural plant-community diversity on soil diversity. Soil fertility differentiates the ways in which interactions between plants and bedrock develop [4,9,28]. On the other hand, the similarities between the environment and plantcommunity interactions suggest the importance of soil influences for biodiversity increase on a supra-regional scale [29]. Diversity was estimated in temperate forests, representing the most complex terrestrial ecosystems directly affected by Quaternary environmental changes [30]. Cultural landscape forests are communities in which the relationship between soil and plant diversity develops naturally, unless it is deviated by non-native tree species' afforestation or environmental change [31]. Such variations in the soil and plant diversity relationship caused by forest management indicate the impact of human activity on ecosystem integrity [32]. The dependence of forest diversity components was assessed in biogeographical regions, whose similarities suggest a supra-regional influence of local relationships [33].
Bioregions are basic territorial units used for individual division of biotic communities. The Czech Republic is divided into 91 bioregions appertained to 2 intervening provinces and 4 subprovinces [38]. The Central European province includes 71 bioregions in the Hercynian subprovince, 4 bioregions in the Polonian subprovince and 11 bioregions in the West Carpathian subprovince, while the Pannonian province includes 5 bioregions in a single North Pannonian subprovince ( Figure 1). The Hercynian subprovince covers 84.7% of the Czech Republic territory and includes 83.3% of the SPTs. In contrast, the Polonian subprovince stretches over the smallest part of the Czech Republic with the total of 2.2% area cover and only 0.8% of protected areas. The West Carpathian subprovince covers almost 9% of the Czech Republic and includes 14.9% of protected territories. The Pannonian province includes the remaining 4.1% of the Czech territory and 1% of the SPTs [39]. Within the Czech Republic, the Hercynian subprovince forest cover amounts to almost 35.9%, the West Carpathian to more than 37.8%, the Polonian to 12.4%, and the North Pannonian comprises a slightly smaller share of 12.3% [34].
the Hercynian subprovince, 4 bioregions in the Polonian subprovince and 11 bioregions in the West Carpathian subprovince, while the Pannonian province includes 5 bioregions in a single North Pannonian subprovince ( Figure 1). The Hercynian subprovince covers 84.7% of the Czech Republic territory and includes 83.3% of the SPTs. In contrast, the Polonian subprovince stretches over the smallest part of the Czech Republic with the total of 2.2% area cover and only 0.8% of protected areas. The West Carpathian subprovince covers almost 9% of the Czech Republic and includes 14.9% of protected territories. The Pannonian province includes the remaining 4.1% of the Czech territory and 1% of the SPTs [39]. Within the Czech Republic, the Hercynian subprovince forest cover amounts to almost 35.9%, the West Carpathian to more than 37.8%, the Polonian to 12.4%, and the North Pannonian comprises a slightly smaller share of 12.3% [34].

Data
Evaluated data consist of forest-division and forest-diversity component classifications. Forests were divided by overlapping biogeographical regions and SPTs. Forestcover distribution was acquired from the geodatabase of regional forest development plans administered by the Forest Management Institute Brandýs nad Labem [40]. Bioregion distribution, including geodiversity components of relief and bedrock-type subpolygons, was adopted from the Unified Environmental Information System (UEIS) administered by the Ministry of the Environment of the Czech Republic [41]. Bedrock-geomorphological bodies (BGMBs), comprehensively capturing the altitudinal variability of the climate and the soil fertility, were derived by overlapping the types of relief and bedrock [42].
Czech SPTs are categorized into reserves and large-area units. While reserves predominantly include natural ecosystem remnants, large protected landscape areas and na-

Data
Evaluated data consist of forest-division and forest-diversity component classifications. Forests were divided by overlapping biogeographical regions and SPTs. Forest-cover distribution was acquired from the geodatabase of regional forest development plans administered by the Forest Management Institute Brandýs nad Labem [40]. Bioregion distribution, including geodiversity components of relief and bedrock-type subpolygons, was adopted from the Unified Environmental Information System (UEIS) administered by the Ministry of the Environment of the Czech Republic [41]. Bedrock-geomorphological bodies (BGMBs), comprehensively capturing the altitudinal variability of the climate and the soil fertility, were derived by overlapping the types of relief and bedrock [42].
Czech SPTs are categorized into reserves and large-area units. While reserves predominantly include natural ecosystem remnants, large protected landscape areas and national parks integrate natural communities with seminatural for renaturalization purposes [43]. Polygons of specially protected territories were acquired from the UEIS geodatabase [44].
Unlike the forest-division systems, forest-diversity components were acquired from the National Forest Inventory (NFI) plot database surveyed in 2011-2015. Forest diversity was acquired from 7550 plots in a 4 × 4 km grid [45]. The grid size was set to capture site effect on temperate vegetation on a regional scale [46]. Inventory plot sites were simplified to BGMB and soil environment. The soil environment was characterized by group and chemical properties of the dominantly developed horizon. Soil groups were evaluated according to the WRB-ISSS-ISRIC [47] and soil chemical properties according to the EMEP-LRTAP criteria [48]. Soil properties were determined as physico-chemical and budget-based. Physico-chemical properties were characterized by pH, determined by acidometry and by sorption properties: cation exchange capacity (CEC) and base saturation (BS) determined by ion-selective extraction. Substance budget determination addressed C org , N tot , Al 2 O 3 , CaO, MgO a P 2 O 5 , which are together with physical and physico-chemical soil properties most significantly differentiated in the Czech Republic [49].
Forest-plant communities were characterized by species abundance of the moss, herb, shrub and tree storeys [2]. The bryophyte composition was simplified to indicatively important species [50]. Vascular plants were defined using a system by Kubát et al. [51], with the exception of the tree-species characterization, which was simplified according to the Czech Forest Management Information Standard [40]. The NFI database includes records on 20 bryophyte indicator species, 815 herbaceous species, 139 species of shrubs and 62 woody species capable of reaching the tree layer.

Statistical Evaluation
The evaluation was divided into an assessment of soil effects and a comparison of management impact on forest plant communities. Soil effects were determined using global linear and geographically weighted regression (GWR) between soil diversity and forest plant species diversity (phytodiversity) at p < 0.05. Soil diversity (pedodiversity) was divided into soil-group diversity and diversity of chemical properties. Influence of pedodiversity components was compared with the geodiversity determined from the BGMB representation [42]. Diversity was estimated using the Shannon-Wiener (H ) index, the calculation of which was adapted to the representation of polygons or to evaluated component numbers from inventory plots. The diversity of BGMBs as well as of the plant species' abundance in the particular forest stand storeys was determined using the equation: where A n is the area of forest ecosystem component type in the bioregion and A is the area sum of the selected ecosystem component types in the bioregion [29]. Plant abundance was measured through H -indices from percentual coverage of individual species in particular stand storey. Total forest phytodiversity (H tfpd ) was derived as the sum of plant species abundances in all stand storeys: where H m is the abundance of indicator bryophytes, H h is the grass-herbaceous storey species abundance, H s is the shrub storey species abundance and H t is the tree species abundance. In contrast, the diversity of soil-group representations or intervals of chemical property values of dominantly developed soil horizons was determined by adjusting: where g k is the interval value count of selected soil property in the bioregion and g is the number of intervals of the whole selected soil property classification range [49]. Similar to H tfpd , soil-chemical diversity (H sc ) was also determined by summing the interval diversities of the individual chemical properties (x): Subsequently, the GWR was approximated either as a simple model for phytodiversity dependence on the BGMBs or soil groups, or as a multiple model with a linear combination of chemical property diversities: where Y m is a dependent variable in the bioregion m, q is the number of independent variables, k is the independent variable order, H mk is the k-squared of independently variable diversity value of soil properties in the bioregion m, β m0 is the intercept parameter in the bioregion m, β mk is the local k-squared independent variable coefficient and m is a random error in the bioregion m [52]. Geographically weighted models were assessed using the Akaike Information Criterion (AIC) and determination indices (R 2 ). The R 2 < 0.5 indices demonstrated lowly significant regressions, while more significant R 2 ≥ 0.5 were divided into five equally wide intervals [53]. Assessment of the forest management impact on the GWR size was based on managed and seminatural ecosystem differentiation. The GWR evaluation was performed using a linear regression of spatially weighted Residual Sum of Squares (RSS) dependence on the specially protected forest territory representation in bioregions. RSS showed the actual detected phytodiversity deviation from the model corresponding with geodiversity. The plant RSS dependence with protected territories confirmed the geodiversity influence on phytodiversity differing between natural and managed forests.

Results
Representation of specially protected forests (SPF) in the Czech Republic is directly proportional to the total forest cover of biogeographical subprovinces. Less than 30% of Czech forests are located in SPTs. The highest proportion of SPFs is found in the West Carpathian subprovince (45.8%). The lowest proportions of SPFs are found the Polonian (7.9%) and North Pannonian subprovinces (10.7%), while the largest Hercynian subprovince includes only 27.9% of SPFs.
The SPF representation has increased the correlations among ecosystem diversity components. The effects of managed forests on correlations were manifest by uneven differences between means and standard deviations (SD) with respect to the overall range of ecosystem properties. The distribution of values did not significantly affect the differences between the mean values and the SD. Area sizes of bedrock-geomorphological bodies in bioregions, soil group totals, soil cation exchange capacity, basic saturation, and the C org , N tot , CaO, MgO a P 2 O 5 contents showed smaller averages than SD. On the contrary, the distribution of soil sorption and the herbaceous to shrub-storey covers of the inventory plots were significantly sloping; however, their averages and SD were balanced even at a wide range of values ( Table 1).
The diversity of vascular forest plants was higher than the diversity of soil environment properties. Out of the total forest-plant diversity, the highest was among the herbal storey and the lowest among the mosses. Despite disrupted normality of the H -indices, the shrub-storey species diversity was on average higher than that of the tree storey. The BGMB diversity and the soil-group representation resembled each other by almost identical H -index ranges. Both groups of soil bodies showed higher diversity than individual soil chemical properties. Soil sorption was the most diversified out of all chemical properties. BS diversity was higher than CEC despite a similar range of H -indices. The diversity of pH was lower than that of soil sorption; however, it was higher than that of mineral-content diversity. Moreover, the diversity of soil organic matter composition was higher than the diversity of mineral contents. Nutrient contents were more diversified than Al 2 O 3 . The diversity of mineral contents decreased in the order of MgO > P 2 O 5 > CaO > Al 2 O 3 .
The total forest phytodiversity was linearly conditioned the most by the BGMB diversity. Although the BGMB diversity in bioregions and the soil-group representation in inventory plots were similar, the diversity of soil chemical properties corresponded more significantly with the soil-group diversity. Plant diversity was linearly less significantly dependent on the soil group representation, but more significantly dependent on soil chemical diversity. In descending order, the total forest plant diversity depended on pH > BS > CEC > CaO > P 2 O 5 . The dependence rate of total plant diversity or species diversity of individual forest ecosystem storeys on soil chemical properties differed due to the unevenness of the cover. The diversity of individual forest plant storeys depended most on the soil CaO content and less on pH, BS, CEC and N tot . However, the species diversity of bryophytes depended most on the BGMBs. The forest herb and shrub species diversity reflected the BGMB and soil chemical property diversity to varying levels. While herbs were more dependent on the BGMB representation, shrubs were more dependent on pedochemical diversity ( Table 2). Geographically weighted regression confirmed global linear dependencies between forest ecosystem components with similar value variability. Spatial dependencies between the representation of polygons within bioregions and the recurrence of properties from inventory plots were differentiated. Due to the differences between bioregions, the total forest-plantdiversity dependence on the geodiversity components was closer than the dependence among particular soil-environment components. However, unlike the relationships found by global linear regression, spatial-plant diversity depended more on the soil group diversity than on the BGMBs. Similarly, the dependence of plant diversity with soil chemical properties was significantly differentiated due to the bioregion dissimilarities (Table 3).
Geographical distribution of the diversity of forest-soil chemical properties simultaneously affects geodiversity and certain parts of phytodiversity. Soil-chemical diversity is most significantly related to the diversity of the forest moss and tree storeys. Seemingly similarly significant spatial dependencies are formed by soil-chemical properties with BGMBs and trees, and by soil groups and bryophytes. Soil-chemical properties influence the soil-group diversity the most. Similarly, the most important relationship between the soil and vegetation occurs between the soil-chemical properties and the bryophyte diversity. However, even similarly significant dependencies differed mostly by parameter opposite proportionalities. Soil-chemical diversity did not significantly determine the vascular-plant diversity distribution into the shrub storey. However, in contrast to the undergrowth storeys, the distribution of forest-tree diversity is significantly dependent on soil-chemical properties. Dependences of the BGMB and the tree-species diversity distribution differed in the parameters of soil pH, sorption properties and mineral contents. On the contrary, the soil group and the bryophyte-diversity dependences were similar due to the soil pH, CEC, BS, C org , CaO and MgO parameters, and differed in N tot , Al 2 O 3 and P 2 O 5 . Only the soil C org content affected the components of ecosystem diversity, within similarly significant regressions, by the same proportion (Table 4). Table 3. Geographically weighted regression among forest-geodiversity components and to total phytodiversity (TPD). BGMB-bedrock-geomorphological body; FSG-forest-soil group; FSC-forest-soil chemistry; A-parameter; A 0 -interception; -random error; w-weight; AQ-adaptive quartile; AIC-Akaike information criterion; R 2determination index.  Table 4. Geographically weighted regression of soil chemical properties diversity with geodiversity components and forest phytodiversity (for detail explanatory of abbreviations see Table 3). Pedodiversity affected forest diversity more than overall geodiversity. Natural forests retain a higher influence of soil groups. In contrast, the influence of the soil-chemicaldiversity components is not increased in the SPTs. Natural environment changes the influence of geodiversity on phytodiversity only to a slight extent. Natural-forest represen-Forests 2021, 12, 534 9 of 15 tation separated the extent and proportion of the geo-and pedodiversity effects. While significant influence of pedodiversity on forest-plant diversity occurs mainly in bioregions with a higher proportion of protected territories, the influence of geodiversity on the diversity of forest plants is slightly higher in bioregions with a smaller proportion of protected territories (Table 5).
Smaller effect of soil-chemical diversity on forest plants has been interchanged by the soil-group representations, although the effect of a linear-chemical-diversity combination on soil-group diversity is more significant than the effect of BGMBs. Geographical dependence between the chemical properties and the forest-soil group diversity attained R 2 in the range of 0.63-0.70. The influence of soil-group diversity on forest-plant diversity attained the R 2 interval of 0.05-0.77. In the Czech Republic, soil-diversity-component dependencies are concentrated in two separate clusters of very high correlation in the west and medium-high correlation north of the Hercynian, North Pannonian and Carpathian subprovince boundary (Figure 2). Similarly, moderate correlations are concentrated in two clusters in the centre of the Hercynian subprovince and in the northeast of the Czech Republic. Very high correlation of pedodiversity includes 33.8% of forests, high correlation 9.1%, medium-high correlation less than 14.2%, most-extensive medium correlation almost 35.4% and moderate correlation only 7.5%. At the same time, the influence of pedodiversity on forest plant diversity is divided into a larger number of clusters separated by decreasing correlations. Very high correlations between pedo-and phytodiversity are found in the west, south and northeast of the Czech Republic. Clusters of significant correlations are segregated by low correlation comprising 43.9% of forests. Very high correlation includes less than 10.7% of forests, high correlation 9.9%, medium-high correlation 14.1%, medium correlation 8.6% and moderate correlation 12.8% (Figure 3). Significant spatial correlation agreement appeared in 11 bioregions in the Hercynian or West Carpathian subprovince. The area of identically evaluated GWRs is forested by better-than-average 40.7% and comprises a 51.9% share of protected forests. found in the west, south and northeast of the Czech Republic. Clusters of significant correlations are segregated by low correlation comprising 43.9% of forests. Very high correlation includes less than 10.7% of forests, high correlation 9.9%, medium-high correlation 14.1%, medium correlation 8.6% and moderate correlation 12.8% (Figure 3). Significant spatial correlation agreement appeared in 11 bioregions in the Hercynian or West Carpathian subprovince. The area of identically evaluated GWRs is forested by better-than-average 40.7% and comprises a 51.9% share of protected forests.

Discussion
Diversity of the soil groups significantly affects the diversity of forest-plant species Figure 2. Geographically weighted regression between forest-soil group and chemical-properties diversity in the Czech Republic.  Plant indication consists of a narrow range of environmental properties in plant species or community adaptation [54]. A comparison of global and localized statistical models pointed to differentiated relationships of forest plants with individual components of geodiversity. Whereas globally forest-plant diversity appeared to be downwardly dependent on the BGMBs, soil-chemical properties, and soil-group representation, the most significant spatially weighted dependence was due to the soil groups and the least significant was due to the soil-chemical properties. Moreover, globally significant influences of the soil physicochemical properties and the CaO and P2O5 contents did not affect forest vegetation in individual bioregions equally. Diversity of the soil-chemical properties most affected the moss and the forest-tree diversity. Nevertheless, their dependencies were affected by species abundancy, but also by the distribution of identified species. Small numbers of species compensated for the spatial differences between bioregions, where the observed dependencies were incorrectly highlighted. On the other hand, the identified dependencies were optimal due to the preserved normal distribution of the species abundance.
Global and localized model differences follow the fact that vegetation in nutrientpoor areas is more strongly differentiated along rock-type transitions, while vegetation in nutrient-rich areas is more strongly differentiated with altitude [2,55,56]. In the Czech Republic, the variability of vegetation dependence on rock or relief type was reflected in the boundary between the Hercynian to Polonian subprovinces and in the North Pannonian

Discussion
Diversity of the soil groups significantly affects the diversity of forest-plant species in temperate conditions. The influence of the soil groups' diversity is more significant than the influence of the bedrock-geomorphological bodies or the soil-chemical properties. Significant impact of pedodiversity on forest plants included more than half of the evaluated Central European forests. The impact of soil diversity is higher in specially protected forests and lower in managed forests. The impact of forest management on the pedo-and phytodiversity association foreshadow the effects on landscape biodiversity. Similarly, variability in the relationships among forest-diversity components affects the ways of evaluating properties of soil-environment indication, ecosystem functions and integrity [44].
Plant indication consists of a narrow range of environmental properties in plant species or community adaptation [54]. A comparison of global and localized statistical models pointed to differentiated relationships of forest plants with individual components of geodiversity. Whereas globally forest-plant diversity appeared to be downwardly dependent on the BGMBs, soil-chemical properties, and soil-group representation, the most significant spatially weighted dependence was due to the soil groups and the least significant was due to the soil-chemical properties. Moreover, globally significant influences of the soil physicochemical properties and the CaO and P 2 O 5 contents did not affect forest vegetation in individual bioregions equally. Diversity of the soil-chemical properties most affected the moss and the forest-tree diversity. Nevertheless, their dependencies were affected by species abundancy, but also by the distribution of identified species. Small numbers of species compensated for the spatial differences between bioregions, where the observed dependencies were incorrectly highlighted. On the other hand, the identified dependencies were optimal due to the preserved normal distribution of the species abundance. Global and localized model differences follow the fact that vegetation in nutrientpoor areas is more strongly differentiated along rock-type transitions, while vegetation in nutrient-rich areas is more strongly differentiated with altitude [2,55,56]. In the Czech Republic, the variability of vegetation dependence on rock or relief type was reflected in the boundary between the Hercynian to Polonian subprovinces and in the North Pannonian to the West Carpathian subprovince regions. The Hercynian to Polonian subprovince area consists of old reliefs on nutrient-poorer bedrock with smaller differences in altitude. The area of the North Pannonian to West Carpathian subprovince is part of the younger alpine orogen with bigger altitude differences [38]. The Hercynian subprovince is characterized mainly by low correlations of forest plants, while the Pannonian-Carpathian subprovince is a transect of significant correlation transitions of forest ecosystem components.
The geographically weighted relationship of forest plants with the environment was closer than the relationship of soil-chemical properties with the soil-group representation. The relationship between forest plants and soil groups was broader than the relationships among soil properties. It reached both low and at the same time the most significant correlations out of the evaluated forest diversity properties. High GWRs of the plant and soil diversity comprised 34.7% of forests in the above-average forested regions. The high GWR concentrations were conditioned not only by the higher numbers of inventory plots enhancing the spatial correlation, but also by the higher proportion of protected territories [39,43,44]. Natural forests significantly maintain the soil-group relationships with plants but reduce the relief effects. Forest management disrupts natural relationships of soils with plants up to a point where they are replaced by more general geodiversity effects. However, the effects of soil-chemical properties on plants were not disrupted by the management. The disruption of relationships with soil groups forced plants to differentiate interactions with the soil-chemical properties along the BGMB boundaries. The effects of BGMBs coincide with the indicator plant species distribution both along the biogeographical subregion boundaries and along the transitions of soil-chemical properties.
The diversity of forest-soil chemical properties was most influenced by base saturation, organic matter content, MgO and P 2 O 5 . The soil carbon content formed equal relations with all the compared components of forest diversity. It evenly affected the plant diversity as well as mirrored the effects of the environment. Such equal proportionality in soil-carbon relationships is due to the plant and the atmosphere connection effects on the nutrient cycles [57]. At the same time, carbon increases forest diversity by supporting other soil properties along with all ecosystem functions [58]. Soil carbon improves physico-chemical sorption and nutrient availability. Its sorption efficiency and nutrient substrate formation increase ecosystem bioproduction and provide habitat for more biodiversity-enhancing organisms [59]. High GWR of pedodiversity components comprised more than 57.1% of the forest stands in the Czech Republic. The high soil GWR occurrences were not affected by large SPTs, which supports the assumption that soil-chemical diversity is unaffected by forestry management.
The soil-group representation has affected the diversity of forest-plant communities most out of the geodiversity components as a result of natural site protection. Yet again, the undisturbed effects of soil-chemical properties by forest management suggest that plants indicate the quality of individual environmental properties equally reliably in differently structured ecosystems. The correlation closeness was higher in bioregions with a higher share of protected forests. The connection between protected forest cover and total forest cover in the Czech Republic is proportionally related to total SPT representation and biogeographical subprovince size [38]. The range of high correlations, on the one hand, adhered to the bioregion forest cover size, and on the other hand was defined along different GWRs between the pedo-and phytodiversity components. Due to the variable GWRs among the forest, diversity components were high correlations isolated into several distant clusters. The degree of dependence between the ecosystem components indicates the effectiveness of ecosystem functions in climate-change adaptation [43]. Similarly, the most influential soil-chemical properties included compounds determinable of forest-health deterioration [59].
The natural dependence closeness between the ecosystem diversity components is the result of ecological integrity [44]. Ecological integrity is the ability of an ecosystem to support a living community appropriate to site conditions on a regional scale [60]. The ecosystem integrity is formed by the species composition, structure and functional connections between species in the community. Promoting integrity consists of processes that increase the ecosystem diversity, in which new food opportunities arise. Ecosystem diversity is increased mainly by [4]. However, in addition to the disturbance recurrences, the degree of integrity is also indicated by the individual component sizes of ecosystem diversity [14]. The spatial dependence of forest-plant diversity on soil diversity indicates the integrity of the community with the bedrock. The correlations of the Central European forest diversity suggest distribution of their response to the variability of growth conditions along a gradient of soils with different fertility [49]. Soils of bedrock-derived fertility are less diversified and display lower resistance to changes in other environmental properties. Soils minerally unrelated to the bedrock are more diversified, though simplified species composition of managed forests reduced the adaptability of the ecosystem to the climate change. Forests with a composition corresponding with the soil diversity, on the other hand, appear to be adaptable.

Conclusions
Soil-group diversity affects the diversity of forest plants in temperate conditions more than the bedrock or soil-chemical properties. Plant diversity and soil-group dependence is significant in areas with an above-average forest cover and a share of protected territories >50%. Nevertheless, plant and soil-chemical property dependence was not affected by forest management. Spatial dependencies pointed to supra-regional overlaps of local relations between forest plants and soils, differentiating environment indication and ecosystem integrity. Environment indication was divided between nutrient-poor soils with low plant dependencies and more fertile soils with higher plant dependencies. Forest plants most distinctly reflect the diversification of soil base saturation and carbon, MgO and P 2 O 5 contents. Spatial-dependency outstretches indicated ecological integrity, which can be consolidated by increasing forest plant diversity in agreement with geodiversity.