Soil Microbial Functional Diversity under the Single-Season Inﬂuence of Traditional Forest Management in a Sessile Oak Forest of Central Europe

: This one-year study focuses on the responses of a soil environment to the implementation of traditional forest management practices in oak–hornbeam stands with the following treatments: cut (C), cut + litter raking (CR), cut + grazing (CG), cut + litter raking + grazing (CRG) and control (Ctrl). The cut was conducted in 2018 through extremely heavy thinning. In autumn of 2017 and 2018, we sampled the soils, focusing on microbial functional diversity (FD) assessments using BIOLOG Ecoplate TM . After one season, the FD was the highest in the Ctrl stand and the lowest in the CRG stand. Furthermore, we detected signiﬁcant seasonal differences in soil reaction, nitrate nitrogen content, phosphatase activity and microbial biomass among the treatments. In particular, the Ctrl stand was deﬁned via FD indices and biochemical and biological soil properties that contrasted mainly with those of the CRG stand deﬁned by the content of mineral nitrogen forms. The soil properties did not differ substantially in the remaining treatments. Of the 31 carbon sources deﬁning FD, 6 were treatment-speciﬁc (putrescine, L-arginine, L-serine, L-threonine, D-cellobiose and glycogen), while the remaining carbon sources mainly displayed either uniform high or low activity across the treatments.


Introduction
Forests are characterized by the presence of long-lived tree species. In connection with these tree species, forest soils typically display arrangements of forest floor horizons (litter). Most organic matter, and thus the biological activity in the soil, is concentrated in the upper soil horizons [1]. Therefore, the upper horizons represent an integral part of the forest soil body [2].
Intensive forest management significantly affects the soil environment. Traditional forest management measures, including coppicing, grazing and litter raking, represent a variety of anthropogenic disturbances [3] that, according to some studies, cause long-term impact on the soil level [4][5][6]. These methods are controversially discussed and frequently abandoned due to the resulting soil degradation and gradual reductions in the production capacities of habitats. According to certain studies, these methods inhibit the succession of forest ecosystems, and thus, in a way, are beneficial in terms of the presence of rare taxa [7]. Some abandoned traditional management types [8] are currently applied again [9]. The purpose is not to acquire natural resources, but to promote species diversity. Some traditional management measures are prohibited in the Czech Republic [10], but their impact on both the soil environment and understory vegetation has been intensively studied [11].
In Central Europe, coppicing was traditional management for centuries [20][21][22]. In the Czech Republic, it was systematically eliminated during the second half of the 20th century [23,24]. Grazing and litter raking are prohibited due to their significant degradation effects on the soil environment. The negative effects of litter raking on nutrient uptake capacity, acid-neutralizing capacity and the C/N ratio of the soil were already known [2] during the first half of the 19th century. Some local contemporary works attempted to assess the soil microbiota response to changing external conditions. Nonetheless, effects of litter raking and grazing on biological activity and the functional diversity of the soil microbial community have not yet been clarified sufficiently.
The rate and turnover intensity of organic matter conversion in forest ecosystems involve participation of soil microorganisms and depend on the composition of the microbial community [25,26]. Therefore, the function of soil microbiota is substantial for the carbon cycle and the cycles of other elements.
The structure and diversity of soil microbial communities reflect climate conditions, weather patterns (temperature and moisture regimes), trophic level, composition and diversity of plant communities [27], anthropogenic influences and the habitat disturbance regime. Higher microbial functional diversity can be expected in species-diverse communities or in communities with equilibrium balance at the food web level [28]. Thus, functionally more diverse ecosystems are more resistant to stress and disturbances. This resistance can be enhanced by mixing species [29][30][31]. In addition, the adaptive capacity of an ecosystem can be increased by alternative forest management measures [32]. On the contrary, disturbances can (temporarily) promote specific functional groups by an increased source of carbon [33]. The exact nature of the interaction between the soil environment-soil microbial functional system and the plant management approach remains unclear.
The soil microbial functional diversity can be expressed as the diversity index of carbon utilization patterns from known carbon sources using the BIOLOG EcoPlate TM system. The BIOLOG EcoPlate TM system expresses the impact of the external environment on microbial functional diversity under various types of environments, management measures and tree species compositions [25,34].
The aim of our study was to assess the influence of coppicing, grazing and litter raking and their combinations on a wide range of soil properties, especially the microbial functional diversity over one season. We formulated the following hypotheses: The effects of particular treatments on soil properties would not differ significantly. The treatments would affect the soil microbial functional diversity and lead to promotion of certain microbial functional groups (as revealed also in Tomao et al.'s study of fungal communities [35]). This would thus lead to overall decrease of diversity indices, whereas in the intact forest stands (control) the diversity index values would be higher.

Study Area
The study area is situated at the Training Forest Enterprise Masaryk Forest Křtiny of the Bílovice forest district, in the southeastern part of the Czech Republic ( According to the Czech forest ecosystem classification [36], the forest sites are mostly classified as nutrient-rich beech-oak (Fageto-Qurcetum eutrophicum), rarely transitioning to nutrient-rich oak-beech (Querceto-Fagetum eutrophicum) on the north-oriented slopes and, to loamy beech-oak (Fageto-Quercetum illimerosum trophicum) on the rare loess loam soils, and to water-deficient beech-oak (Fageto-Quercetum subxerothermicum) forest site complexes on the very rare sunny slopes.

Experimental Design
In 2017, we established six 60 m × 40 m experimental blocks (blocks 1-6; Figure 1). In January 2018, extremely heavy thinning (C-cut) was performed on each block with a buffer zone of 20 m (for the elimination of edge effects) with an intensity of 88%, and the remaining 12% of trees were preserved as standards (trees preserved on the plot after thinning) [39]. Each experimental block was divided into two 30 m × 40 m plots (plots A and B, 12 plots in total). Blocks 2, 3 and 5 were assigned as cut + grazed on the 2B, 3B and 5B plots (CG) and as cut + grazed + raked on the 2A, 3A and 5A plots (CRG). Blocks 1, 4 and 6 were assigned as cut on the 1B, 4B and 6B plots (C) and as cut + raked on the 1A, 4A and 6A plots (CR). The assignment of these treatment to the blocks was performed randomly. We also established three 30 m × 40 m control plots (representing undisturbed forests) assigned as Ctrl.
The litter (L, F and part of the H horizons) was raked in the CR and CRG plots in April 2018. The total weight of the dry mass of the removed litter was 30.5 t ha −1 (sd 8.74). From June to September 2018, blocks 2, 3 and 5 were grazed by Šumavka breed sheep under the conditions of one 14-day rotation (4-pc herd) and two 10-day rotations (6-pc herd).

Field and Laboratory Work
The soil taxonomic units were defined (1) using deep soil probes (10 in total) to sample the soil down to the soil-forming substrate and distinguish the diagnostic horizons (41 samples in total) and (2) using a probe rod (26 in total) to the verify soil taxonomic units according to the World Reference Base [37].
Each of the 15 plots was divided into four quadrants for soil sampling. From each quadrant, one mixed soil sample was collected (four samples per plot, 60 samples in total) from the 25 cm × 25 cm mineral A horizon subsurface area (depth of 5-10 cm). In this sampling survey, we avoided a 15 m × 15 m area in each plot center, which was reserved for geobotanical inventory (not included in this work). We also marked the sampling quadrants in each plot and for each sample to allow comparisons between 2017 and 2018. The samples were taken in autumn 2017 (prior to thinning) and autumn 2018 (one season after thinning). The functional diversity was determined in 2018, and the other soil properties were analyzed in 2017 and 2018.
The activity of soil proteases was determined by hydrolysis with casein as a substrate, and the amount of l-tyrosine produced was measured according to the methodology by Rejšek et al. [47]. The ammonium and nitrate nitrogen (N-NH 4 + and N-NO 3 − ) contents in the soil extracts were determined using 1 M KCl extract according to the method described by Kučera et al. [48]. The urease activity of the soil was determined by the Kandeler and Gerber method [49] after incubation of the urea samples. The acid phosphomonoesterase activity was determined using the Rejšek method [50]. The determination of the amount of microbial carbon (C mic ) in the soil biomass was performed by the fumigation-extraction method according to Zbíral et al. [43].
The microbial functional diversity was assessed using the BIOLOG EcoPlate TM [51,52] technique with 31 organic substrates (each with 3 repetitions per plate) used as carbon sources for microbial community cultivation. The carbon compounds were coupled to capture electrons by colorless tetrazolium salts, forming reduced purple formazans that could be readily measured spectrophotometrically via absorbance. The degree of substrate oxidation was monitored through the reduction of a tetrazolium dye. The degree of staining depended on the amount and activity of the microbial cells. The microplate inoculation was based on a fresh soil sample (10 g) shaken for 20 min (at 200 rounds per min) with 90 mL 0.40% NaCl solution and centrifuged for 8 min at 4000 rpm. The supernatant was diluted 1:1000 and 150 µL of this extract was incubated in BIOLOG EcoPlate TM at 23 • C for 384 h. The absorbance was recorded every 12 h during the first two days and every 24 h during the remaining incubation period using a Tecan Infinte M200Pro Nano+ microplate reader at 590 nm against a control well (pure water) as a blank.

Data Analysis
The data were processed using Statistica CZ 12 software (with graphical outputs of post hoc Tukey multiple comparison test) [53] and R software for statistical computing (R version 3.6.2 and RStudio version 1.3.1056) (remaining statistical analysis) [54]. The functional diversities obtained based on the microplates were assessed using the (1) Shannon-Wiener diversity index (S-W index), calculated as follows [55,56]: where p i is the mean absorbance of the three repetitions of the substrate per microplate. The average well color development index (AWCD) (2) was calculated as follows: where A i is the absorbance within each well and A 0 is the absorbance of the control well.
For the data processing, we used the whole data set of n = 60 (12 values for each treatment and 4 values for each plot).
We used paired sample t-tests to compare the soil parameter values obtained in 2017 and 2018. The pairs were determined using the zone sampled within each sampling plot in 2017 or 2018. Parametric ANOVA followed by Tukey s HSD test (with the treatments set as categorical variables) was employed to compare (1) the soil properties in each year separately and (2) the diversity indices (S-W index, AWCD) at each absorbance reading time of ecoplate incubation. For the diversity indices, we finally used the 384 h incubation period with the most significant differences in the diversity index values. A two-tailed correlation analysis was performed for n = 60 samples (Pearson critical value was 0.250) using the corrplot function as a graphical output. Principal component analysis (PCA) was performed with the vegan package version 2.5-6 [57] after data standardization using the scale function. For the PCA, the treatments were applied for categorical variables and for continuous variables: (1) the S-W index and AWCD; (2) the soil properties determined in 2018 (Nt_2018, etc.); and (3) the differences in soil properties between 2018 and 2017, set as the 'delta' values (Nt_delta, etc.). The number of variables was reduced in the instance of a significant correlation, as assessed from the ordination plot. In 2018, and again in the 2018-2017 delta soil data, we followed an equal list of soil properties (AWCD, S-W index, N t , DegHum, BS, AS, K-rat.law, N-NH 4 + , N-NO 3 and C mic ). The individual treatment factor levels are shown on the ordination plot of the first two axes using the ordispider function. The functions ggplot and qplot ( ggplot2 package version 3.3.0 [58]) were applied as graphic instruments to show the extent to which each substrate contributed to the total index. The pattern of utilization of the 31 carbon substrates was scaled to 5 categories in relation to the AWCD index plate (1: <1st quartile; 2: <median; 3: AWCD; 4: <3rd quartile; 5: >3rd quartile). All the tests were performed at the significance level of alpha = 0.05.

Differences between 2017 and 2018
Before the thinning, the soils were typical with slightly acidic soil reactions, a slight predominance of C-FA, and a contribution of approximately 25% of humus substances to the total organic matter content ( Table 1). The nutrient contents corresponded to the mesotrophic habitat as well as the development of humus forms, with a slight predominance of bivalent base cations over potassium and optimal base saturation in mineral A horizon above 50%.
The results of the paired t-test detected differences between 2017 and 2018 in the monitored parameters within the treatments. In all plots, increases in N t , Ca and N-NO 3 and decreases in H + and acid phosphomonoesterase activity between 2017 and 2018 could be observed. Statistically significant increases were observed in the soil response associated with the CRG and Ctrl plots, an increase in the humus fraction content was observed in connection with the CRG plots, and an increase in N-NH 4 + within the C between 2017 and 2018. The greatest increase in the N-NO 3 concentration occurred in the context of the CRG (nearly a tenfold increase), and the N-NH 4 + concentration increased almost comparably (at the statistically significant threshold) in conjunction with the CG and CRG plots (cf. Ctrl without any significance). Table 1. The results of the paired t-test comparison of the soil parameter values in relation to the individual treatments in 2018 and 2017. n-s > α; 0.05-0.1; * 0.01-0.05; ** 0.001-0.01; *** < 0.001.

Specifications of the Treatments
ANOVA showed significant differences among the treatments for the delta values of five soil parameters between 2017 and 2018 ( Table 2). The greatest increases in the soil reaction values occurred within the CRG and Ctrl plots, while in connection with the CR plots, a very slight decrease was measured. The largest increase in the nitrate nitrogen concentration was detected in the CRG plots, and the smallest increase was observed in the Ctrl plots. The smallest decreases in the acid phosphomonoesterase activity with respect to the baseline were shown by the CRG and Ctrl plots, while the greatest decrease was observed in the C plots. The greatest increase in the Cmic content was recorded in the CRG plots, and the greatest decrease was recorded in the C plots. The CRG evinced the highest annual change rate in all instances where the ANOVA null hypothesis was rejected, while in Ctrl plots the response did not show such an unequivocal trend.
The treatments differed in the spatial heterogeneity of the blocks (Table 3), which was systematically higher in pH and lower in Cmic in 2017 than in 2018. The variability increased in the concentration of nitrate nitrogen, especially within plots and after their aggregation to blocks, except for the CG plots. The highest increase in variability was observed in the CRG plots, especially in the 3A and 2A plots. The variability in phosphatase activity decreased, especially in plots within the CR, CRG and Ctrl sites. Both diversity indices were the most homogenous in the Ctrl plots in 2018 (Table 3), in which the variability was the lowest and most balanced. In contrast, the highest AWCD and S-W variabilities were in the CRG plots.

Diversity Indices and Multiple Relations
The diversity index development over the whole incubation time is shown in Figure S1. The AWCD index showed the highest values in the Ctrl sites during the whole incubation time and the lowest values since the 168th hour of reading in the CRG. At 384 h of incubation, both the AWCD and S-W indices showed similar results ( Table 2). Similar and symmetrical distributions of these values were observed especially in the AWCD values, with the lowest variability and the highest median recorded in the Ctrl plots. The S-W index was characterized by a slightly right-skewed distribution under the CR treatment. The highest soil spatial heterogeneity was in the CRG plots ( Figure 2).   PC2 was positively correlated with enzymatic activity (represented by Cmic), with the functional diversity indices and with K-ratio law. PC2 was negatively correlated with the mineral nitrogen contents. The treatment levels were clustered around the center of the factorial plane except for the Ctrl plots, which corresponded with the positive PC2 values.   The measurement of substrate absorbance participation on the AWCD index is shown graphically in Figure 5 and in a numerical form in the Supplementary Material (Table S1). The substrate consumption was found to be heterogeneous mainly for compound groups of amino acids (L-arginine, L-serine, L-threonine) and then for putrescine, D-cellobiose and glycogen, and partially also for 4-hydroxybenzoic acid and D-galactonic acid-gammalactone. As demonstrated in Figure 6 and Figure S2, the weakest reliance of diversity index on activity of particular substrate was in Ctrl, which contradicted specifically with the CRG plots as well as with all the management practices. This situation was analogous in the dependency of the S-W index on the overall activity of the compound groups ( Figure S3).

Discussion
The microbial communities of forest soils are affected by several factors [59][60][61]. In our study, we measured early responses of soil parameters and soil microbial functional diversity to traditional forest management practices. The outcomes of this study confirm our hypotheses that in many soil parameters the effects of different management practices are ambiguous after one season. This is due to the different robustness or 'sensitivity' levels of soil properties. At the same time, management changes result in reductions in microbial functional diversity. As the outputs of this short-term study were not directly recognizable, we explain the results, including several issues, in the sections that follow.

Soil under Influence of Traditional Management
The traditional management of forests in the European region, both due to the intensity and duration of management measures, has substantially influenced fertility and overall soil chemistry in a long-term outlook.
According to numerous authors [62][63][64][65][66], coppicing leads to root-to-shoot ratio imbalances, the thinning of organic soil layers due to accelerated mineralization and, consequently, nutrient losses from the forest floor. Coppicing is considered to be the cause of soil degradation due to acidification and soil genetic processes leading to accelerated soil profile differentiation, as reported in historical studies starting in the 1950s [23,67]. On the other hand, Hölscher et al. [68] showed increases in both pH and the concentrations of exchangeable cations in a 20-year-old coppice forest on acidic soils in the organo-mineral horizons, corresponding with our results, albeit in a shorter time interval. Coppicing does not limit the accessibility of nutrients in soils, nor does it impoverish the soil environment [68].
A study from Mediterranean oak stands in Italy [69] also showed increased pH after three years of coppicing, as opposed to the organic carbon content. This increase in pH was explained by decreased root exudation associated with increasing H + concentration in soil solution, and by release of nutrients from organic bonds in the forest floor. The work also reported the response of soil microbiota to coppicing in terms of their functional diversity. The microbiota responded to active management by increasing their amount of microbial biomass and their respiratory and enzymatic activity by up to 47.8%. In our research, increase in Cmic was observed in three of the five treatments, including the nonintervention treatment, but these increases were not statistically significant.
Similar conclusions were presented in a short-term study from hardwood forests near the town of Saranac, New York, that were converted to pastures and silvopastures [18]. During the two years following these conversions, the percentage of nitrogen and phosphorus increased compared to the woodlot and pretreatment levels, while the observed decrease in pH was not statistically significant. Similarly, other nutrient contents did not manifest significant trends that were dependent on the treatment, even when compared with the pretreatment values.
Litter raking may be equivalent in effect to acidification by atmospheric deposition in the 20th century [16]. In the study conducted in the Slavkov Forest and Ore Mountains (Czech Republic) with Norway spruce (Picea abies (L.) H. Karst.), one litter raking depleted the soil exchangeable base cation pool by up to 31%, and the overall nutrient loss exceeded the annual inputs by atmospheric deposition and weathering. Some authors mention the possibility of compensating for soil depletion due to litter raking through atmospheric nitrogen deposition [2,7]. Strong pH decreases and base cation leaching following litter removal were observed in the leaf-rich Luzulo-Fagetum habitat of the Schwarzwald forest [70], which we could not confirm via single-season data. However, these acidification and related degradation processes were also confirmed by Dzwonko and Gawroński [2], who investigated the effect of litter raking on soils in a 16-year-old study. They revealed the expected substantial nutrient impoverishment of the soil (mainly P, Mg and Ca together with the cation exchange capacity) in the H and A horizons and the lessivage layer.
The impacts of litter raking not only include nutritional impoverishment but also include changes in organic carbon inputs, soil moisture levels, surface temperatures, soil surface fauna and understory vegetation. In our study, litter raking combined with coppicing had an insignificant effect on the diversity index after one season. The lowest microbial functional diversity was found merely in the plots in which all three types of management were combined.

Factors Influencing Soil Response Measure
In our study, increased concentrations of phosphorus and bivalent bases were detected after one year, with the highest increase observed in the CRG plots.
Despite the general assumption that traditional management significantly affects the soil environment, in our treatise, management cannot be clearly attributed to this influence for the following reasons: (1) The soil environment is characterized not only by temporal heterogeneity but also by spatial heterogeneity within blocks and plots. Another reason (2) is the characteristics of short-term data, as the changes in soil did not yet display significant differences over the short time. The third reason (3) involves microclimatic conditions, i.e., the extremely low precipitation totals and soil moisture levels recorded in 2017 and 2018.

Soil Spatial Diversity
Variability in the data (Table 3) was mostly higher in the instance of blocks (where n = 12) than among plots (where n = 4). This result has also been confirmed by several authors [18,[71][72][73]. The spatial heterogeneity of blocks was higher in 2017 than in 2018 (Table 3), especially in the CR plot. We attribute this result to the intensive disturbance conducted in this plot, which led to the homogenization of soil under equal intensities of external factors. We attribute the prominent increase in N-NO 3 variability in the CRG block to the reduced root exudation [74] and to the additional input of N from feces. Compared with the remaining soil properties, the S-W index demonstrates a more homogenous soil environment (expressed by the carbon source consumption) in the Ctrl plot; this plot showed very similar values when comparing the variability among blocks and plots.

Temporal Aspect of Changes in Soil Environment
This study monitored the immediate response of the soil environment. Management changes had a systemic effect on some soil properties after one season, while other properties showed no apparent trends.
Responses in microbial biomass and metabolic activity were monitored by Chen et al. [61] following a management change after more than five years, and significant differences were measured using a paired comparison design. The results showed a decrease of microbial biomass due to increased nutrient uptake by trees. This trend clearly did not occur with Cmic in our study; however, it can be expected in the subsequent hypotheses. We noted distinctive decreases in phosphatase, urease and protease activities and in microbial biomass in C (42, 27, 29 and 27%, respectively) during one season, despite the other treatments in which the control confirmed no disturbance and reflected common stand conditions and soil seasonal dynamics.
As also reported by Bürgi and Gimmi [12], no rapid soil chemistry changes can be expected after litter raking. At the same time, we lack long-term data on the response of the soil environment (a fortiori, soil microbiota) to traditional management practices, as studies of traditional management have thus far given insufficient attention to the plant-soil relationship from the landscape perspective [14,22]. Soil changes over a 40-year horizon were addressed in four different soil taxonomical units by Hédl and Rejšek [75] at the soil texture, soil reaction, carbonates and organic carbon content levels. In addition to the incompatibility of laboratory methods, which represent a problem in pedological long-term studies [76,77], the authors concluded that the temporal changes were insignificant, and neither the leaching of carbonates nor decreases in soil reaction joined with acidification or soil environment impoverishment occurred due to nutrient depletions. Alterations in soil properties on a physicochemical basis can be expected with an uncertain trend due to the buffering capacity [78].

Microclimate Limitations
Changes in the climate and seasonal dynamics can influence the common behaviors of soil microbiota in the context of disturbances and forestry adaptation strategies. In our study, the soil moisture fell below 10% during the 2018 growing season [71]. The nature of the causes of low soil moisture varied between 2017 (when the low soil moisture originated from below-average precipitation and was escalated by the water suction function of the forest stand involved) and 2018 (when the water deficit came mainly from low precipitation totals). Thus, the lack of water as a key soil medium did not mediate sufficient cation exchanges, as did the expected eradication of nutrient leaching in the mineralization process.
The relationship between functional diversity and ecological stability, underlined by environmental conditions (e.g., the moisture regime), was specified in Hallett's [79] study focused on grazed ecosystems. Functional diversity was found to rapidly decline under wet conditions of grazed systems, while it was maintained (not elevated) under drier conditions (cf., [80]). In addition, grazed systems with greater functional diversity have been found to be more environmentally stable across varied moisture conditions. The responses of microbial communities are markedly associated with environmental conditions but less with management (as in our research paper) or plant species compositions, according to Howes study [81]. Nevertheless, he worked with the soil-limiting systems of Michigan (USA) dunes, i.e., a more preliminary system development. He revealed significant relations between functional diversity and the organic matter content, temperature, distance from a forest and surface detritus thickness.
The environmental limitation of microbial activity was also discussed in a study by Rosíková et al. [82]. In that work, the amelioration effect (effect of forest management practice) expressed as CO 2 efflux from the soil surface was obscured by environmental limitations such as drought and low temperatures. Hence, the results showed substantial seasonal dynamics that, however, could not be attributed to the treatments of the amelioration experiment.

Substrate Consumption Level, Diversity Indices and Metabolic Profile
The approaches used to assess soil functional diversity differ depending on the use of different methods. The abilities and limits of the EcoPlate TM method have been discussed in several studies [83,84]. In general, EcoPlate TM metabolic profiles provide a better understanding of in situ ecological and ecosystem functions that are based on metabolic activities. The BIOLOG EcoPlate TM system expresses the metabolism of carbon sources as a key factor associated with the development of microbial biomass [52,85,86]. However, fungi are only included in the metabolic profile to a limited extent (they experience slow growth due to the incubation duration). In addition, we obtained information only about culturable microorganisms, while during incubation, changes may have occurred in the structure of the microbial community. This is also why we applied the determination of the enzymatic activity. The EcoPlate TM system better reflects the presence of low-molecularweight organic compounds in the soil that are bound to the exudation of the physiologically active root system. Enzymatic activity captures the ability of the soil to hydrolyze complex polymers into simple compounds, capturing the state of the soil following disturbances associated with extremely heavy thinning (88% of the stock volume in our case) or other management changes (in our case, litter raking and grazing).
In a three-year-long study by Pignataro et al. [86] conducted in central Italy (Quercus cerris spp.), a decrease in the total organic carbon and an increase in pH were observed. The AWCD index was higher in aged coppice stands than in three-year-old coppice stands. In contrast, the activity of microbial enzymes increased in the coppiced soil, suggesting that different analytical approaches target diverse ecological processes depending on the quality of the available substrates. In contrast, the results of our study suggest higher diversity indices in the nonintervention treatment, demonstrating nonextremeness and the balance of food and metabolic soil relations at the microbiological level. As we stated in our hypotheses, the disturbances led to reductions in the diversity indices and thus reductions in specific functional groups of soil microorganisms, but the enzymatic activity itself did not yet show any clear response.
Regarding the partial substrate consumption uniformity, the results of which contradict those of the Pignataro study [87], the majority of substrates revealed rather similar reactions when the treatments were compared. The exceptions include the substrates mentioned above, which were consumed differently among the treatments. Some carbon sources were detected as displaying either uniform high activity or low activity, for which the uniform presence or absence of the substrate consumers was detected, respectively ( Figure 6; see also Figure S2 for the carbon sources and Figure S3 for the compound groups).

Multivariate Soil Responses and Disturbance Regime
The relations between disturbance events and soil microbial community responses were studied by Gömöryová et al. [88] on windthrow plots in the Tatra Mountains of Slovakia. They revealed the highest microbial activity and biomass in reference and fire plots with burnt organic matter in comparison with extracted and nonextracted plots. However, the trend of microbial response development corresponded mainly with time after the windthrow season × treatment interactions. The results were attributed to 'unstressed' microbial communities in the undisturbed reference plot and to a reduction in the surface organic layer and enrichment of organo-mineral horizons by released nutrients in the burnt plot. Compared with our results, the microbial diversity did not show any significant differences despite such extremely different habitats.
The complexity of soil systems cannot be denied, as shown in Figure 4 via the correlation matrices. The correlation matrices demonstrate the acidifying effect of the increased nitrate nitrogen concentration. The correlation was significant before treatments (Figure 4a) but not after the disturbances (Figure 4b). The correlation of the S-W index with the initial soil state was slightly statistically significantly ( Figure 4a) positive with the K-ratio law and, subsequently, negative with the Ca/K and Mg/K ratios at similar correlation coefficient values (0.28; −0.31 and −0.31, respectively). In all three cases, these results indicate a positive response in the diversity index in relation to the amount of potassium in the soil. The S-W index was also (insignificantly) positively correlated with the C-HA/FA ratio and the degree of humification (with coefficients of 0.20 and 0.22, respectively). Simultaneously, the AWCD index was significantly positively correlated with the calcium content.
Likewise, in the study of Vuong et al. [89], strong correlations between soil diversity indices and enzyme activities were not detected. In 2017, the enzyme activity was positively correlated with the K, C org and N t contents, with CEC and especially with pH/H 2 O, which also confirmed the findings in the studies of Błońska et al. [90] (urease activity), Adamczyk et al. [91] (protease, phosphatase activities) and Klose and Tabatabai [92] (phosphatase activity). In contrast, negative correlations between the soil reaction and enzyme activity were described by Tan et al. [93], and specifically with phosphatase activity by Šarapatka and Kršková [94]. Therefore, as reported in other studies [86], the chemical properties and the characteristics of the soil sorption complex, which indicate the characteristics of soil chemistry, together determine the microbial environment and the state and composition of the soil microbial communities. Regarding Figure 4b, the relationships substantially depend on the measure of disturbance, leading to decreases or even to changes in the signs of the correlations.
Likewise, in the work of Pignataro [87], treatments influenced the functional diversity by causing deviations from the balanced microbial community composition with a high diversity index. Numerous studies have reported a direct link between biological diversity and ecosystem stability [52,95,96]. This relation can be perceived as a feedback mechanism. Therefore, disturbances may result in reduced ecosystem stability and hence in reduced soil microbial functional diversity.
In our study, six months after the experiment was established, the individual treatments could still be perceived simply as thinned areas. A gradual increase in the importance of the treatments for the spatial differentiation of the soil environment in relation to specific management practices can be expected. The changes were linked to the disturbances and to the consequent mineralization of organic matter in the forest floor or the enrichment of nitrogen components by nitrification or grazing [97] (the changes in Cmic were statistically insignificant in our study). We also expect the following changes in the composition and metabolic profile of the microbial community. Zhang et al. [97] noted these changes after three years in semiarid temperate steppe conditions, which characterize a water-limiting environment in terms of biological activity. In our study, this trend was already demonstrated after a short time response (Figure 3), in which the concentration values of the mineral forms of nitrogen and the diversity index values were negatively correlated.
Although some differences could already be seen after such a short period of time, it was not possible to predict clear or unchanging trends in the dynamics of soil properties with certainty [64]. However, silvopastures are perceived by many authors [18,19] as beneficial and preserving in terms of agroforestry systems being an adaptive strategy under global climate change.

Conclusions
This study focused on the responses of soils to extremely heavy thinning (88% of stock volume) and subsequent litter raking and grazing. An alteration in microbial functional diversity occurred after one observation season. Microbial functional diversity was reduced compared to the parental forest stand and was differentiated according to the treatments in some parameters. In addition to the microbial functional diversity, the interventions also affected other soil characteristics, such as soil reaction (significant increase in stand in which a combination of thinning, litter raking and grazing was applied), nitrate nitrogen content (similar reactions with more distinct differentiations related to the nonintervention status of the stand), acid phosphomonoesterase activity (decreased activity was observed, especially in the thinned stands) and the microbial carbon content (the most significant increases compared to the nonintervention state were observed at the plots where thinning, litter raking and grazing were combined).
A better comprehension of soil metabolic processes contributes to the consideration of appropriate adaptation strategies during environmental changes. In subsequent studies, more unequivocal and increasing differentiations in relation to individual traditional management types can be expected. Further studies may also be focused on the closely related water regime of the forest soil in a particular habitat, which also closely affects soil microbiota activities.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.