Next Article in Journal
Evidence of Constrained Divergence and Conservatism in Climatic Niches of the Temperate Maples (Acer L.)
Previous Article in Journal
Caffeine Interactions with Wood Polymers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Department of Geology and Soil Science, Faculty of Forestry and Wood Technology, Mendel University in Brno, Zemědělská 3, CZ-613 00 Brno, Czech Republic
2
Global Change Research Institute CAS, Belidla 986/4a, CZ-603 00 Brno, Czech Republic
3
Forest Management Institute Brandýs nad Labem, Nábřežní 1326, CZ-250 01 Brandýs nad Labem, Czech Republic
*
Author to whom correspondence should be addressed.
Forests 2021, 12(5), 534; https://doi.org/10.3390/f12050534
Submission received: 30 March 2021 / Revised: 23 April 2021 / Accepted: 24 April 2021 / Published: 26 April 2021
(This article belongs to the Section Forest Ecology and Management)

Abstract

:
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 influences 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 influenced 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 significantly dependent on soil-group representation. The soil-group effect is more significant 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.

1. 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 on solar radiation, biodiversity derives from the growth of plants. Environmental effects on plants systematize biodiversity through the diversification of food opportunities and habitats for animals. Terrestrial environment effects are limited by the continent size within which biodiversity is climate-controlled [5]. Biodiversity evolves through community isolation, landscape penetrability and habitat diversity. Continent biodiversity is directly proportional to the number of endemics resulting from a sufficiently long isolation period. Larger continents, therefore, provide greater environmental diversity with more unique sites of greater total ecosystem [6].
While the climate drives generalized trends in the progression of biodiversity on the greatest scale, in smaller areas many local effects apply. Altitudinal zonation divides the species composition of plant communities more significantly than horizontal climatic transitions [7]. Site-specific effects can completely reverse the trend of biodiversity progression [8]. The relationship between environmental diversity and biodiversity is either direct or indirect, depending on whether the progression of the ecosystem is predominantly influenced by climate or local characteristics [9]. The relationship proportionality between the environment and biodiversity is differentiated between areas affected by ice ages and not affected by Quaternary glaciation [1]. In areas without the direct effect of glaciation, low geodiversity is accompanied by high biodiversity and vice versa. In areas directly affected by environmental changes during glaciation, low geodiversity is accompanied by low biodiversity, and high geodiversity is accompanied by high biodiversity [10].
Different impacts of Quaternary climate changes amidst tropical-to-polar latitudes have differentiated the development coherence between plant-community diversity and geodiversity into declining proportionality between species abundance and environmental diversity [11]. The greatest plant diversity remained concentrated across humid tropics and subtropics, where it constitutes a response to the most significant transitions of environmental properties [5]. The mildest transitions between environmental properties and plant communities are found in the subpolar latitudes, where plant communities are also the most homogeneous in the largest areas [12]. The polar regions are not distinguished by any plant geoelements, yet the northern temperate zone border includes up to 80% of relevant geoelements. In contrast, the southern temperate boundaries include only 15.5% of plant genera due to dominant species overlap from different climatic regions. Tropical regions include up to 80% of genera found on the continent [13].
Biodiversity regulates productivity and protects the ecosystem from damage during climate change [14]. Productivity is modified by functionally substitutable organisms and by the selection of species in competition. Substitutability and selection increase community competition that supports ecosystem resilience to environmental change [15]. Biodiversity protects the ecosystem proportionally to the number of species that complementarily stabilize nutrient accessing while increasing the productivity [16]. Ecosystem resilience consists in the length, asynchronicity and individual species response form to environmental fluctuations. Effective protection is manifest by delaying productivity fluctuation effects [17]. Ecosystem resilience increase in spite of changes in the environment is achieved through asynchronous intraspecific response, differential species response rates and decreased competition. The intra- and interspecies response processes take place as a temporal substitutability of species, while the decrease in competition is controlled by substitutability in functional groups of organisms [18]. Due to these protective processes, changes in species composition have greater impacts on ecosystems than changes in the abiotic environment [19].
Although biodiversity shapes ecosystems through functional groups of organisms, its impact is not uniform due to the persisting dependence on the supra-regional to local properties of the abiotic environment. Such feedback between biodiversity and the abiotic environment complicate distinction among the causes of variability at ecosystem characteristics. The abiotic environment and repeated disturbances affect most ecosystem characteristics more than plant species’ diversity [20]. Plant diversity is not as effective in increasing ecosystem diversity as dominant species that independently support nutrient 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 transition-community 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 plant-community 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].

2. Material and Methods

2.1. Investigated Area

The dependence of forest-plant species’ diversity on soil properties was evaluated among Czech biogeographical regions (bioregions). The Czech Republic is a Central European country in the transient range of temperate climate (78,866 km2, 115‒1602 m n.m., 2.0–9.4 °C; 250–1470 mm) [34]. Potential forest cover of the Czech Republic exceeds 98%, and natural forest-less land occurs in alpine, rocky or peatland areas. Current forest cover is 34%. Natural representation of woody plants comprises European beech (Fagus sylvatica) (40%), oaks (Quercus sp.) (18%), Silver fir (Abies alba) (16%), Norway spruce (Picea abies) (15%) and Scots pine (Pinus sylvestris) (3%) [35], while actual composition is formed by 44% of Norway spruce, almost 10% of Scots pine, more than 10% of European beech and by 8% of oaks [36]. Approximately 16.7% of the Czech Republic was designated as Specially Protected Territory (SPT), out of which almost 60.5% is under forest cover [37].
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].

2.2. 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 Corg, Ntot, Al2O3, CaO, MgO a P2O5, 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.

2.3. 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:
H = ( A n A ) · log 2 ( A n A ) ,
where An 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:
Htfpd = Hm + Hh + Hs + Ht
where Hm is the abundance of indicator bryophytes, Hh is the grass-herbaceous storey species abundance, Hs is the shrub storey species abundance and Ht 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:
H f = ( g k g ) · log 2 ( g k g ) ,
where gk 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 Htfpd, soil–chemical diversity (Hsc) was also determined by summing the interval diversities of the individual chemical properties (x):
H sc = H 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:
Y m = β m 0 + k = 1 q β m k · H m k + ϵ m
where Ym is a dependent variable in the bioregion m, q is the number of independent variables, k is the independent variable order, Hmk 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 (R2). The R2 < 0.5 indices demonstrated lowly significant regressions, while more significant R2 ≥ 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.

3. 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 Corg, Ntot, CaO, MgO a P2O5 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 Al2O3. The diversity of mineral contents decreased in the order of MgO > P2O5 > CaO > Al2O3.
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 > P2O5. 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 Ntot. 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-plant-diversity 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, Corg, CaO and MgO parameters, and differed in Ntot, Al2O3 and P2O5. Only the soil Corg content affected the components of ecosystem diversity, within similarly significant regressions, by the same proportion (Table 4).
Pedodiversity affected forest diversity more than overall geodiversity. Natural forests retain a higher influence of soil groups. In contrast, the influence of the soil-chemical-diversity components is not increased in the SPTs. Natural environment changes the influence of geodiversity on phytodiversity only to a slight extent. Natural-forest representation 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 R2 in the range of 0.63–0.70. The influence of soil-group diversity on forest-plant diversity attained the R2 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.

4. 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 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 nutrient-poor 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 P2O5. 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.

5. 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 P2O5 contents. Spatial-dependency outstretches indicated ecological integrity, which can be consolidated by increasing forest plant diversity in agreement with geodiversity.

Author Contributions

Conceptualization, P.S., P.C. and M.K.; methodology, M.K.; formal analysis, P.S.; resources, P.C.; data curation, M.K.; writing—original draft preparation, P.S.; writing—review and editing, J.V.; visualization, P.S.; supervision, P.C.; project administration, P.C.; funding acquisition, P.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded on the basis of Regulation 247/2009 by Government of the Czech Republic.

Acknowledgments

The authors gratefully acknowledge the support by the Project LM2018123 CzeCOS of the Ministry of Education, Youth and Sports of the Czech Republic.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Thom, D.; Rammer, W.; Dirnböck, T.; Müller, J.; Kobler, J.; Katzensteiner, K.; Helm, N.; Seidl, R. The impacts of climate change and disturbance on spatio-temporal trajectories of biodiversity in a temperate forest landscape. J. Appl. Ecol. 2017, 54, 28–38. [Google Scholar] [CrossRef]
  2. Vávrová, E.; Cudlín, O.; Vavříček, D.; Cudlín, P. Ground vegetation dynamics in mountain spruce (Picea abies/L./Karsten) forests recovering after air pollution stress impact. Plant Ecol. 2009, 205, 305–321. [Google Scholar] [CrossRef]
  3. Pimm, S.L. Biodiversity: Climate Change or Habitat Loss—Which Will Kill More Species? Curr. Biol. 2008, 18, R117–R119. [Google Scholar] [CrossRef] [Green Version]
  4. De La Riva, E.G.; Lloret, F.; Pérez-Ramos, I.M.; Marañón, T.; Saura-Mas, S.; Díaz-Delgado, R.; Villar, R. The importance of functional diversity in the stability of Mediterranean shrubland communities after the impact of extreme climatic events. J. Plant Ecol. 2017, 10, 281–293. [Google Scholar] [CrossRef] [Green Version]
  5. Kier, G.; Mutke, J.; Dinerstein, E.; Ricketts, T.H.; Küper, W.; Kreft, H.; Barthlott, W. Global patterns of plant diversity and floristic knowledge. J. Biogeogr. 2005, 32, 1107–1116. [Google Scholar] [CrossRef]
  6. Whittaker, R.J.; Willis, K.J.; Field, R. Scale and species richness: Towards a general hierarchical theory of species diversity. J. Biogeogr. 2001, 28, 453–470. [Google Scholar] [CrossRef] [Green Version]
  7. Cox, C.B. The biogeographic regions reconsidered. J. Biogeogr. 2001, 28, 511–523. [Google Scholar]
  8. Carnaval, A.C.; Moritz, C. Historical climate modelling predicts patterns of current biodiversity in the Brazilian Atlantic forest. J. Biogeogr. 2008, 35, 1187–1201. [Google Scholar] [CrossRef]
  9. Zeng, C.; Yang, L.; Zhu, A.-X.; Rossiter, D.G.; Liu, J.; Liu, J.; Qin, C.; Wang, D. Mapping soil organic matter concentration at different scales using a mixed geographically weighted regression method. Geoderma 2016, 281, 69–82. [Google Scholar] [CrossRef]
  10. Barthlott, W.; Biedinger, N.; Braun, G.; Feig, F.; Kier, G.; Mutke, J. Terminological and methodological aspects of the mapping and analysis of global biodiversity. Acta Bot. Fenn. 1999, 162, 103–110. [Google Scholar]
  11. Weigelt, P.; Steinbauer, M.J.; Cabral, J.S.; Kreft, H. Late Quaternary climate change shapes island biodiversity. Nature 2016, 532, 99–102. [Google Scholar] [CrossRef]
  12. Olson, D.M.; Dinerstein, E.; Wikramanayake, E.D.; Burgess, N.D.; Powell, G.V.N.; Underwood, E.C.; D’Amico, J.A.; Itoua, I.; Strand, H.E.; Morrison, J.C.; et al. Terrestrial Ecoregions of the World: A New Map of Life on Earth. BioScience 2001, 51, 933–938. [Google Scholar] [CrossRef]
  13. Qian, H.; Song, J.-S.; Krestov, P.; Guo, Q.; Wu, Z.; Shen, X.; Guo, X. Large-scale phytogeographical patterns in East Asia in relation to latitudinal and climatic gradients. J. Biogeogr. 2003, 30, 129–141. [Google Scholar] [CrossRef] [Green Version]
  14. Hisano, M.; Searle, E.B.; Chen, H.Y.H. Biodiversity as a solution to mitigate climate change impacts on the functioning of forest ecosystems. Biol. Rev. 2018, 93, 439–456. [Google Scholar] [CrossRef]
  15. Loreau, M. Linking biodiversity and ecosystems: Towards a unifying ecological theory. Philos. Trans. R. Soc. B 2010, 365, 49–60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Hooper, D.U.; Chapin, F.S., III; Ewel, J.J.; Hector, A.; Inchausti, P.; Lavorel, S.; Lawton, J.H.; Lodge, D.M.; Loreau, M.; Naeem, S.; et al. Effects of biodiversity on ecosystem functioning: A concensus of current knowledge. Ecol. Monogr. 2005, 75, 3–35. [Google Scholar] [CrossRef]
  17. Yachi, S.; Loreau, M. Biodiversity and ecosystem productivity in a fluctuating environment: The insurance hypothesis. Proc. Natl. Acad. Sci. USA 1999, 96, 1463–1468. [Google Scholar] [CrossRef] [Green Version]
  18. Loreau, M.; de Mazancourt, C. Biodiversity and ecosystem stability: A synthesis of underlying mechanisms. Ecol. Lett. 2013, 16, 106–115. [Google Scholar] [CrossRef]
  19. Chapin, F.S., III; Zavaleta, E.S.; Eviner, V.T.; Naylor, R.L.; Vitousek, P.M.; Reynolds, H.L.; Hooper, D.U.; Lavorel, S.; Sala, O.E.; Hobbie, S.E.; et al. Consequences of changing biodiversity. Nature 2000, 405, 234–242. [Google Scholar] [CrossRef]
  20. Loreau, M. Ecosystem development explained by competition within and between material cycles. Proc. R. Soc. Lond. B 1998, 265, 33–38. [Google Scholar] [CrossRef]
  21. Augusto, L.; Ranger, J.; Binkley, D.; Rothe, A. Impact of several common tree species of European temperate forests on soil fertility. Ann. For. Sci. 2002, 59, 233–253. [Google Scholar] [CrossRef] [Green Version]
  22. Dempster, M.; Dunlop, P.; Scheib, A.; Cooper, M. Principal component analysis of the geochemistry of soil developed on till in Northern Ireland. J. Maps 2013, 9, 373–389. [Google Scholar] [CrossRef] [Green Version]
  23. Solomou, A.D.; Sfougaris, A. Predicting Woody Plant Diversity as Key Component of Ecosystems: A Case Study in Central Greece. Int. J. Agric. Environ. Inf. Syst. 2019, 10, 1–20. [Google Scholar] [CrossRef] [Green Version]
  24. Ortiz, A.M.D.; Outhwaite, C.L.; Dalin, C.; Newbold, T. A review of the interactions between biodiversity, agriculture, climate change, and international trade: Research and policy priorities. One Earth 2021, 4, 88–101. [Google Scholar] [CrossRef]
  25. Proutsos, N.; Solomou, A.; Karetsos, G.; Tsagari, K.; Mantakas, G.; Kaoukis, K.; Bourletsikas, A.; Lyrintzis, G. The Ecological Status of Juniperus foetidissima Forest Stands in the Mt. Oiti-Natura 2000 Site in Greece. Sustainability 2021, 13, 3544. [Google Scholar] [CrossRef]
  26. Giam, X. Global biodiversity loss from tropical deforestation. Proc. Natl. Acad. Sci. USA 2017, 114, 5775–5777. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Plexida, S.; Solomou, A.; Poirazidis, K.; Sfougaris, A. Factors affecting biodiversity in agrosylvopastoral ecosystems with in the Mediterranean Basin: A systematic review. J. Arid Environ. 2018, 151, 125–133. [Google Scholar] [CrossRef]
  28. Ibañez, J.J.; De-Alba, S.; Bermúdez, F.F.; García-Álvarez, A. Pedodiversity concepts and tools. Catena 1995, 24, 215–232. [Google Scholar] [CrossRef]
  29. Guo, Y.; Gong, P.; Amundson, R. Pedodiversity in the United States of America. Geoderma 1995, 117, 99–115. [Google Scholar] [CrossRef] [Green Version]
  30. Bohn, F.J.; Huth, A. The importance of forest structure to biodiversity–productivity relationships. R. Soc. Open Sci. 2017, 4, 160521. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Lukac, M. Soil biodiversity and environmental change in European forests. Cent. Eur. For. J. 2017, 63, 59–65. [Google Scholar] [CrossRef] [Green Version]
  32. Schröter, D.; Cramer, W.; Leemans, R.; Prentice, I.C.; Araújo, M.B.; Arnell, N.W.; Bondeau, A.; Bugmann, H.; Carter, T.R.; Gracia, C.A.; et al. Ecosystem Service Supply and Vulnerability to Global Change in Europe. Science 2005, 310, 1333–1337. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Fotheringham, A.S.; Brunsdon, C.; Charlton, M. Geographically Weighted Regression: The Analysis of Spatially Varying Relationships; John Wiley & Sons: Chichester, UK, 2002. [Google Scholar]
  34. Chytrý, M. Vegetation of the Czech Republic: Diversity, ecology, history and dynamics. Preslia 2012, 84, 427–504. [Google Scholar]
  35. Neuhäuslová, Z.; Blažková, D.; Grulich, V.; Husová, M.; Chytrý, M.; Jeník, J.; Jirásek, J.; Kolbek, J.; Kropáč, Z.; Ložek, V.; et al. Map of Potential Natural Vegetation of the Czech Republic; Academia: Prague, Czech Republic, 1998. [Google Scholar]
  36. Zouhar, V. Dřevinná skladba. In Oblastní Plány Rozvoje Lesů 2; ÚHÚL Brandýs nad Labem, Česká lesnická společnost, Ministerstvo zemědělství České Republiky: Praha, Czech Republic, 2018; pp. 20–22. [Google Scholar]
  37. Céza, V.; Čermáková, E.; Kochová, T.; Mertl, J.; Pokorný, J.; Přech, J.; Rollerová, M.; Vlčková, V. Zpráva o Životním Prostředí České Republiky 2018; CENIA: Praha, Czech Republic, 2019. [Google Scholar]
  38. Culek, M. Biogeographical provinces, subprovinces and bioregions of the Czech Republic. J. Landsc. Ecol. 2013, 6, 5–16. [Google Scholar] [CrossRef] [Green Version]
  39. Divíšek, J.; Chytrý, M.; Grulich, V.; Poláková, L. Landscape classification of the Czech Republic based on the distribution of natural habitats. Preslia 2014, 86, 209–231. [Google Scholar]
  40. Bošeľa, M.; Redmond, J.; Kučera, M.; Marim, G.; Adolt, R.; Gschwantner, T.; Petráš, R.; Korhonen, K.; Kuliešis, A.; Kulbokas, G.; et al. Stem quality assessment in European National Forest Inventories: An opportunity for harmonised reporting? Ann. For. Sci. 2016, 73, 635–648. [Google Scholar] [CrossRef] [Green Version]
  41. Vlčková, V.; Buček, A.; Machar, I.; Kiliánová, H. The Application of Geobiocoenological Landscape Typology in the Modelling of Climate Change Implications. J. Landsc. Ecol. 2015, 8, 69–81. [Google Scholar] [CrossRef] [Green Version]
  42. Kubalíková, L. Geomorphological heritage and Geoconservation in the Czech Republic. In Landscapes and Landforms of the Czech Republic; Pánek, T., Hradecký, J., Eds.; Springer International Publishing: Berlin/Heidelberg, Germany, 2016; pp. 387–398. [Google Scholar]
  43. Simon, J. Management strategies on territories with special status of protection in the Czech Republic. J. For. Sci. 2004, 50, 510–513. [Google Scholar] [CrossRef] [Green Version]
  44. Cudlín, P.; Pechanec, V.; Štěrbová, L.; Cudlín, O.; Purkyt, J. Integrated approach to the mitigation of biodiversity lost in Central Europe. In Ecological Integrity and Land Use. Sovereignty, Governance, Displacements and Land Grabs; Westra, L., Bosselmann, K., Zabrano, V., Eds.; Nova Science Publishers: New York, NY, USA, 2019; pp. 75–86. [Google Scholar]
  45. Adolt, R.; Ene, L.; Fejfar, J.; Kohn, I.; Morneau, F.; Pesty, B.; Reidel, T.; Lanz, A. Scalable pan-European model-assisted biomass estimation. In A Century of National Forest Inventories—Informing Past, Present and Future Decisions; Tomppo, E., McRoberts, R., Fernandez, C.A., Alberdi, I., Breidenbach, J., Eds.; Norsk Institutt for Bioøkonomi Ås: Ås, Norway, 2019; Available online: https://nibio.pameldingssystem.no/auto/1/NFI-100-abstract/220_abstract-oslo-2019_v2.pdf (accessed on 22 May 2019).
  46. Affleck, D.L.R.; Gregoire, T.G.; Valentine, H.T. Edge effects in line intersect sampling with segmented transects. J. Agric. Biol. Environ. Stat. 2005, 10, 460–477. [Google Scholar] [CrossRef]
  47. Schad, P.; van Huysteen, C.; Michéli, E.; Vargas, R. (Eds.) World Reference Base for Soil Resources 2014. International Soil Classification System for Naming Soils and Creating Legends for Soil Maps; World Soil Resources Reports No. 106; FAO: Rome, Italy, 2014. [Google Scholar]
  48. Gauger, T.; Anshelm, F.; Schuster, H.; Erisman, J.W.; Vermeulen, A.T.; Draaijers, G.P.J.; Bleeker, A.; Nagel, H.D. Mapping of Ecosystem Specific Long-Term Trends in Deposition Loads and Concentrations of Air Pollutants in Germany and Their Comparison with Critical Loads and Critical Levels; Umweltbundesamt: Berlin, Germany, 2002. [Google Scholar]
  49. Samec, P.; Zapletal, M.; Horáček, M. Forest transformation urgency for topsoil diversity optimization during environmental change. J. Landsc. Ecol. 2020, 13, 82–106. [Google Scholar]
  50. Frego, K.A. Bryophytes as potential indicators of forest integrity. For. Ecol. Manag. 2007, 242, 65–75. [Google Scholar] [CrossRef]
  51. Kubát, K.; Hrouda, L.; Chrtek, J., Jr.; Kaplan, Z.; Kirchner, J.; Štěpánek, J. (Eds.) Klíč ke Květeně České Republiky; Academia: Praha, Czech Republic, 2002. [Google Scholar]
  52. Gollini, I.; Lu, B.; Charlton, M.; Brundson, C.; Harris, P. GWmodel: An R Package for Exploring Spatial Heterogeneity using Geographically Weighted Models. J. Stat. Softw. 2015, 63, 1–52. [Google Scholar] [CrossRef] [Green Version]
  53. Finstad, K. Response Interpolation and Scale Sensitivity: Evidence against 5-Point Scales. J. Usability Stud. 2010, 5, 104–110. [Google Scholar]
  54. Pauli, H.; Gottfried, M.; Dullinger, S.; Abdaladze, O.; Akhalkatsi, M.; Alonso, J.L.B.; Coldea, G.; Dick, J.; Erchbamer, B.; Calzado, R.F.; et al. Recent Plant Diversity Changes on Europe’s Mountain Summits. Science 2012, 336, 353–355. [Google Scholar] [CrossRef] [Green Version]
  55. Gégout, J.-C.; Križová, E. Comparison of indicator values of forest understorey plant species in Western Carpathians (Slovakia) and Vosges Mountains (France). For. Ecol. Manag. 2003, 182, 1–11. [Google Scholar]
  56. Sporbert, M.; Keil, P.; Seidler, G.; Bruelheide, H.; Jandt, U.; Aćić, S.; Biurrun, I.; Campos, J.A.; Čarni, A.; Chytrý, M.; et al. Testing macroecological abundance patterns: The relationship between local abundance and range size, range position and climatic suitability among European vascular plants. J. Biogeogr. 2020, 47, 2210–2222. [Google Scholar] [CrossRef]
  57. Lal, R. Forest soils and carbon sequestration. For. Ecol. Manag. 2005, 220, 242–258. [Google Scholar] [CrossRef]
  58. Cools, N.; Vesterdal, L.; De Vos, B.; Vanguelova, E.; Hansen, K. Tree species is the major factor explaining C:N ratios in European forest soils. For. Ecol. Manag. 2014, 311, 3–16. [Google Scholar] [CrossRef]
  59. Bakker, M.R.; Brunner, I.; Ashwood, F.; Bjarnadottir, B.; Bolger, T.; Børja, I.; Carnol, M.; Cudlín, P.; Dalsgaard, L.; Erktan, A.; et al. Belowground Biodiversity Relates Positively to Ecosystem Services of European Forests. Front. For. Glob. Chang. 2019, 2, 6. [Google Scholar] [CrossRef] [Green Version]
  60. Wurtzebach, Z.; Schultz, C. Measuring Ecological Integrity: History, Practical Applications, and Research Opportunities. BioScience 2016, 66, 446–457. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Distribution of forest reserves, protected landscape areas (PLA) and national parks (NP) in the Czech Republic.
Figure 1. Distribution of forest reserves, protected landscape areas (PLA) and national parks (NP) in the Czech Republic.
Forests 12 00534 g001
Figure 2. Geographically weighted regression between forest-soil group and chemical-properties diversity in the Czech Republic.
Figure 2. Geographically weighted regression between forest-soil group and chemical-properties diversity in the Czech Republic.
Forests 12 00534 g002
Figure 3. Geographically weighted regression between forest-soil group and plant diversity in the Czech Republic.
Figure 3. Geographically weighted regression between forest-soil group and plant diversity in the Czech Republic.
Forests 12 00534 g003
Table 1. Diversity intervals of compared forest ecosystem characteristics in biogeographical regions of the Czech Republic. BGMB—bedrock-geomorphological body; FSG—forest soil group; FSC—forest soil chemistry; CEC—cation exchange capacity; BS—base saturation; H′—Shannon–Wiener index; SD—standard deviation; E—elevation; A—skewness (bold significantly disrupted data normality at p < 0.05).
Table 1. Diversity intervals of compared forest ecosystem characteristics in biogeographical regions of the Czech Republic. BGMB—bedrock-geomorphological body; FSG—forest soil group; FSC—forest soil chemistry; CEC—cation exchange capacity; BS—base saturation; H′—Shannon–Wiener index; SD—standard deviation; E—elevation; A—skewness (bold significantly disrupted data normality at p < 0.05).
DiversityAverage ± SDMin–MaxH′ ± SDHminHmaxEA
BGMB (km2)22.13 ± 43.490.00–506.912.56 ± 0.880.01–4.170.34−0.72
FSG (n)25 ± 431–4092.50 ± 0.760.00–4.410.41−0.26
FSC (1)--10.45 ± 2.400.00–15.312.81−0.80
pH [−logc(H+)]4.84 ± 0.864.11–7.571.07 ± 0.460.00–1.92−0.12−0.44
CEC (cmol+/kg)12.08 ± 12.272.70–421.101.38 ± 0.420.00–2.122.04−1.15
BS (%)33.71 ± 33.737.40–99.951.66 ± 0.470.00–2.263.41−1.69
Corg (%)2.02 ± 3.920.29–9.911.29 ± 0.530.00–2.36−0.60−0.24
Nt (%)0.13 ± 0.160.03–0.321.41 ± 0.440.00–2.360.44−0.45
Al2O3 (g/kg)29.71 ± 17.988.35–81.070.53 ± 0.570.00–2.120.251.04
CaO (g/kg)3.08 ± 17.890.13–87.700.80 ± 0.570.00–2.00−0.890.28
MgO (g/kg)5.71 ± 7.520.92–17.801.18 ± 0.610.00–2.28−0.58−0.29
P2O5 (g/kg)0.75 ± 0.820.21–3.681.12 ± 0.400.00–2.250.93−0.28
Tree (%)44.17 ± 10.4524.49–69.192.65 ± 0.711.14–4.06−0.55−0.22
Shrub (%)18.43 ± 9.901.17–57.903.32 ± 0.630.59–4.312.90−1.13
Herb (%)17.93 ± 7.414.69–37.614.54 ± 0.761.44–5.882.79−1.29
Moss (%)21.04 ± 15.241.75–67.001.70 ± 0.820.00–3.03−0.29−0.74
Plant (%)76.05 ± 28.3547.75–100.0012.20 ± 1.854.85–15.791.73−0.93
Table 2. Linear correlations between compared Czech forest ecosystem diversity components (bold significant at p < 0.05). BGMB—bedrock-geomorphological body; FSG—forest soil group; FSC—forest soil chemistry; CEC—cation exchange capacity; BS—base saturation.
Table 2. Linear correlations between compared Czech forest ecosystem diversity components (bold significant at p < 0.05). BGMB—bedrock-geomorphological body; FSG—forest soil group; FSC—forest soil chemistry; CEC—cation exchange capacity; BS—base saturation.
ComponentDiversityBGMBFSGFSCpHCECBSCorgNtAl2O3CaOMgOP2O5
AbioticBGMB-0.290.18−0.020.360.400.060.13−0.24−0.050.050.40
FSG--0.620.190.490.430.440.430.36−0.030.560.14
BioticTree−0.18−0.160.150.380.08−0.12−0.040.060.080.51−0.13−0.17
Shrub0.370.190.440.510.340.380.040.220.030.410.090.19
Herb0.480.090.270.330.270.310.080.22−0.180.30−0.090.19
Moss0.590.310.15−0.080.210.430.050.03−0.02−0.340.250.32
Plant0.520.180.380.420.350.400.050.20−0.040.310.050.22
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; A0—interception; ϵ—random error; w—weight; AQ—adaptive quartile; AIC—Akaike information criterion; R2—determination index.
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; A0—interception; ϵ—random error; w—weight; AQ—adaptive quartile; AIC—Akaike information criterion; R2—determination index.
XYAA0ϵwAQAICR2
BGMBFSG0.22 ± 0.061.95 ± 0.190.02 ± 0.6951.49 ± 3.950.362280.11 ± 0.04
FSC0.55 ± 0.789.13 ± 2.22−0.06 ± 2.0214.67 ± 2.460.084670.28 ± 0.05
TPD1.14 ± 0.739.38 ± 2.000.00 ± 1.209.23 ± 2.090.053630.59 ± 0.19
Soil groupFSC2.04 ± 0.055.18 ± 0.10−0.01 ± 1.8772.38 ± 1.790.634400.39 ± 0.01
TPD0.36 ± 0.8611.10 ± 2.490.09 ± 1.146.31 ± 1.840.033650.62 ± 0.18
Soil chemistryTPD0.34 ± 0.308.54 ± 3.700.02 ± 1.198.82 ± 1.970.043640.59 ± 0.20
Table 4. Geographically weighted regression of soil chemical properties diversity with geodiversity components and forest phytodiversity (for detail explanatory of abbreviations see Table 3).
Table 4. Geographically weighted regression of soil chemical properties diversity with geodiversity components and forest phytodiversity (for detail explanatory of abbreviations see Table 3).
FormulaAtributeBGMBFSGTreeShrubHerbMossPlant
ParameterpH−0.38 ± 0.21−0.23 ± 0.070.40 ± 0.280.40 ± 0.120.14 ± 0.19−0.19 ± 0.290.74 ± 0.37
CEC0.74 ± 0.110.37 ± 0.10−0.04 ± 0.270.17 ± 0.060.19 ± 0.100.38 ± 0.480.65 ± 0.11
BS0.47 ± 0.070.43 ± 0.03−0.21 ± 0.180.22 ± 0.040.33 ± 0.120.57 ± 0.380.91 ± 0.25
Corg−0.22 ± 0.270.47 ± 0.10−0.09 ± 0.13−0.32 ± 0.02−0.20 ± 0.140.08 ± 0.45−0.54 ± 0.51
Ntot0.17 ± 0.280.16 ± 0.110.04 ± 0.270.29 ± 0.090.26 ± 0.22−0.15 ± 0.410.42 ± 0.70
Al2O3−0.05 ± 0.070.11 ± 0.080.18 ± 0.21−0.14 ± 0.05−0.33 ± 0.16−0.18 ± 0.28−0.47 ± 0.25
CaO−0.20 ± 0.04−0.05 ± 0.030.40 ± 0.080.15 ± 0.010.17 ± 0.03−0.54 ± 0.390.24 ± 0.19
MgO0.11 ± 0.150.42 ± 0.09−0.19 ± 0.130.06 ± 0.09−0.08 ± 0.240.16 ± 0.16−0.03 ± 0.40
P2O50.32 ± 0.08−0.21 ± 0.12−0.24 ± 0.070.06 ± 0.070.05 ± 0.170.44 ± 0.280.28 ± 0.36
01.20 ± 0.360.41 ± 0.112.78 ± 0.252.15 ± 0.103.58 ± 0.090.41 ± 0.368.96 ± 0.28
ϵ0.00 ± 0.550.01 ± 0.420.00 ± 0.45−0.01 ± 0.47−0.02 ± 0.62−0.01 ± 0.47−0.04 ± 1.38
Characteristicsw53.24 ± 3.7447.38 ± 4.1738.74 ± 4.2362.91 ± 3.1153.11 ± 3.7737.66 ± 4.1355.57 ± 3.53
AQ0.380.320.240.480.370.240.40
AIC194137156158220165389
R20.58 ± 0.050.67 ± 0.020.57 ± 0.060.46 ± 0.060.42 ± 0.100.66 ± 0.080.45 ± 0.06
Table 5. Linear regression of specially protected forest (SPF) representation effect on residual range from relationship between geo- and phytodiversity.
Table 5. Linear regression of specially protected forest (SPF) representation effect on residual range from relationship between geo- and phytodiversity.
GeodiversityR2FpSPFParameter
BGMB0.010.530.470.000310.444
Soil group0.011.330.25−0.000650.591
Soil chemistry0.000.020.89−0.000050.502
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Samec, P.; Volánek, J.; Kučera, M.; Cudlín, P. Effect of Soil Diversity on Forest Plant Species Abundance: A Case Study from Central-European Highlands. Forests 2021, 12, 534. https://doi.org/10.3390/f12050534

AMA Style

Samec P, Volánek J, Kučera M, Cudlín P. Effect of Soil Diversity on Forest Plant Species Abundance: A Case Study from Central-European Highlands. Forests. 2021; 12(5):534. https://doi.org/10.3390/f12050534

Chicago/Turabian Style

Samec, Pavel, Jiří Volánek, Miloš Kučera, and Pavel Cudlín. 2021. "Effect of Soil Diversity on Forest Plant Species Abundance: A Case Study from Central-European Highlands" Forests 12, no. 5: 534. https://doi.org/10.3390/f12050534

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop