Next Article in Journal
Long-Term Soil Fertility and Site Productivity in Stem-Only and Whole-Tree Harvested Stands in Boreal Forest of Quebec (Canada)
Previous Article in Journal
Stakeholders’ Perception of the Impact of the Declaration of New Protected Areas on the Development of the Regions Concerned, Case Study: Czech Republic
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Short-Term Impacts of Harvesting Intensity on the Upper Soil Layers in High Karst Dinaric Fir-Beech Forests

1
Chair for Forest and Urban Greenery Ecology, Faculty of Forestry, University of Sarajevo, Zagrebačka 20, 71000 Sarajevo, Bosnia and Herzegovina
2
Slovenian Forestry Institute, Večna pot 2, 1000 Ljubljana, Slovenia
3
Department of Forestry and Renewable Forest Resources, Biotechnical Faculty, University of Ljubljana, Večna pot 83, 1000 Ljubljana, Slovenia
4
Faculty of Agriculture and Food Science, University of Sarajevo, Zmaja od Bosne 8, 71000 Sarajevo, Bosnia and Herzegovina
*
Author to whom correspondence should be addressed.
Forests 2021, 12(5), 581; https://doi.org/10.3390/f12050581
Submission received: 25 March 2021 / Revised: 26 April 2021 / Accepted: 29 April 2021 / Published: 6 May 2021
(This article belongs to the Section Forest Ecology and Management)

Abstract

:
The present study addresses the short-term effects of different harvest intensities under close-to-nature selective management on the upper soil layers in Slovenian and Bosnian Dinaric karst fir-beech forests. The different harvest intensities coincided with the single-tree and irregular shelterwood management, common in the region. The effect of harvesting intensity on the upper soil layers (Ol, Of, Ol and 0–10 cm mineral soil) was investigated by a repeated measurements experiment in Slovenia on 27 research plots in close-to nature managed forests. The properties of the upper layers (concentration of SOC and TN, C/N ratio, weights, BD and SOC stocks) were analyzed twice, before (2011) and after (2014) treatment of 50% and 100% harvest intensity in relation to the total standing growing stock of trees. As a control, we used no-treatment <20% harvesting intensity plots. To extend this experiment, we added three comparable plots from the Bosnian site: one in an old-growth forest with 0% harvest intensity and two in the managed forest with <20% harvest intensity. The results of the assessment of mean differences indicated a significant influence of harvesting intensity on the decrease in SOC, TN concentrations, weights and SOC stocks in the organic layers and the increase in BD and SOC stocks in the 0–10 cm mineral soil. The highest relative decreases in Ol, Of and Oh SOC stocks occurred in 50% (−10 and −38%) and 100% (−16 and −49%) harvest intensities. Negligible relative differences in both organic and 0–10 cm mineral layers were found for the <20% harvest intensity in the region. The change in forest light conditions resulting from differences in canopy openness as a function of applied harvest intensity explained the significant difference in the properties of the upper soil layers. The impact of the short-term losses in SOC stocks, in terms of overall soil productivity, may depend on the regeneration dynamics and melioration methods.

1. Introduction

The well-preserved condition of high karst Dinaric fir-beech forests is closely related to their low management intensity in the past [1]. Their uneven-aged structure resembles that of old-growth reserves [2], with elevation being the key factor controlling the microclimate in these temperate mountain forest stands [3]. Liberal selection of felling regimes in uneven-aged fir-beech forests, also known as close-to-nature silviculture [4], employs relatively low-intensity and small-scale felling regimes to mimic natural forest composition, structures and natural disturbances at the lower end of the stand-level disturbance intensity gradient [1,5], which is practiced in some European countries.
The result of such a forest management approach can be seen in high Dinaric karst fir-beech forests, which have an uneven-stand structure, multiple canopy layers, horizontal heterogeneity and a considerable amount of deadwood [6,7,8,9]. Compared to other geological substrates, karst ecosystems are considered vulnerable due to the rock–soil system [10] and forests growing on a calcareous bedrock, sensitive to changes in microclimatic conditions and extreme weather events [11,12]. However, the effects of different levels of stand changes in this region have not been adequately studied.
Soil organic matter properties are particularly important for the productivity of shallow mountain soils on limestone and dolomite in the Dinaric region. Soil properties, especially total soil organic carbon (SOC) and total nitrogen (TN), are used as elementary predictors of forest stand stability [13]. Losses in SOC stocks resulting from abiotic and biotic changes influence the decline in soil function. Although soil organic matter has some inherent resilience to environmental change through stabilization mechanisms [14,15], the upper layers, including the humus layer and upper mineral soil, respond to a change in stand conditions. Short-term or long-term reductions in the SOC stock depend on the degree of site disturbance, which varies between harvest intensities in the humus layer [16] and miner soil layers [17]. Despite the importance of knowing the current status and changes of SOC stocks depending on the applied forest management practices, little is known in this regard in the Dinaric karst region [18].
Selective harvesting and practices that promote continuous canopy cover are usually considered to maintain and increase SOC content. Modeled data show that the selection system causes an increase in SOC stocks [19]. The comparison between single-tree-selection management and conventional age-class management in Norway spruce stands in Austria showed an 11% increase in SOC stocks in the upper mineral soil [20]. Applied shelterwood harvesting systems in the Lenga forest of Chilean Patagonia reduced SOC at short-term levels [21]. Single-tree selective and irregular shelterwood harvests should increase the rate of decomposition at the short-term level without causing long-term SOC declines. There is an important influence by increasing the mineralization rate on forest stand recovery and increasing productivity. Because there is no clear evidence on how different harvesting practices under selective management affect soil organic matter properties, this is investigated in this study.
For the Dinaric region, it was reported that the SOC stock in 0.1-hectare canopy gaps decreases, compared to the old-growth reserves [22]. The creation of small canopy gaps can lead to a significant change in soil conditions [23]. Moreover, SOC declination has been reported in mountain forest soils of the Bavarian limestone Alps, where the single-tree selection system was applied [24,25]. We hypothesize that high-intensity harvesting, equivalent to 50 and 100% removal, even if applied to a small forest area, would significantly affect SOC stocks in high karst Dinaric forests. We hypothesize that a single-tree selection vs. an irregular shelterwood system, corresponding to different levels of harvest intensity, would cause different levels of site change, reflecting SOC storage. Therefore, we believe it is critical to evaluate the effect of different harvest intensities on the upper soil layers that are most susceptible to change. Understanding the mechanisms of post-harvest carbon loss is far from complete and requires deeper understanding [26], especially under the threat of extreme weather events and global climate change.
Studies of how soils on limestone and dolomite are affected by harvesting are rare, and more data and studies are needed for a better understanding of the SOC dynamic in these areas. Because soils on limestone and dolomite are regularly shallow and rocky, they have a lower water storage capacity [27,28,29]. Stable soil conditions are maintained under a continuous canopy cover, which may be an indicator of changes in soils. Thus, the forest canopy determines the light and thermal microsite conditions [30,31] and hence the SOC dynamics. It influences the soil organic layer through the amount of litterfall [32] and the mineral SOC content through root debris [33]. Studies show that higher concentrations of soil organic matter and mineral N coincide with tree canopy density [34]. Canopy openness, as a function of harvest, intensity can therefore be considered as an indicator of environmental change at the site, which is crucial for understanding the nature of such a relationship in assessing optimal openness for SOC conservation.
In the present study, it was hypothesized that the quantity and quality of SOM would be similar in all plots studied (1) and that the intensity of harvest would be reflected in light conditions (canopy openness, site light factors) and upper soil layer properties (pH values, SOC, TN concentrations, C/N ratio, SOC stock). The objectives of the present study are (1) to determine the quantity and quality of soil properties in the upper soil layers; (2) to evaluate the effect of harvest intensity on the upper soil layers, especially on SOC stocks and (3) to characterize the relationship between canopy openness, light conditions and upper soil layer properties, i.e., pH values, SOC, TN concentrations, C/N ratio and SOC stocks.

2. Materials and Methods

2.1. Site Description

The study sites were located in the high karst Dinaric fir-beech forests of Bosnia and Herzegovina (BA) and Slovenia (SL), starting from the southernmost site on the Bjelašnica Mountain (1420 m a.s.l.), through Kočevski Rog (871 m a.s.l.), Snežnik (870 m a.s.l.) to Trnovski Gozd (814 m a.s.l.) (Figure 1). All sites had similar ecological factors.
On each of these sites, the predominant tree species were European beech (Fagus sylvatica L.), silver fir (Abies alba Mill.) and Norway spruce (Picea abies (L.) H. Karst), belonging to the Omphalodo-Fagetum association. Each site is characterized by high plant diversity [35,36]. A comparable forest structure was attributed to each site, which is related to the applied close-to-nature forest management. A canopy cover of 95%, necessary for natural regeneration, is the typical characteristic of forest canopy conditions (Table 1). Furthermore, forest sites with growing stocks between 270 and 442 m3 ha−1 and increments between 6.2 and 9.4 m3 ha−1 yr−1 were described (BA: Forest Management Base for the Igmansko forest area for the period 2004 to 2014; SL: [37]). Climate and soil data indicate similar site conditions. Mean annual temperatures (MATs) varied between 7.6 and 11.3 °C and mean annual precipitation (MAP) between 1192 and 1619 mm for the period 1985–2015 (meteorological data were obtained from the website Royal Netherlands Meteorological Institute “Climate Explorer” http://climexp.knmi.nl, accessed on 2015). The sites were attributed to calcareous and dolomitic bedrock and the typical soil type complex of Leptosol, Rendzic Leptosol, Cambisol and Luvisol in BA [38] and in SL [37]. The soils found at the sites were predominantly shallow to moderately deep, rocky, rich in soil organic matter and clay content.
Harvesting intensity, which, in practice, depends on the single-tree or irregular shelterwood harvesting method, influenced the formation of small canopy gaps typical for close-to-nature forest management. The close-to-nature silviculture differs from plantation forestry and clear-cutting [4]; it is characterized by site-appropriate species composition, avoidance of clear-cutting, focus on stability, reliance on natural processes and focus on the development of individual trees and mixed and uneven-aged, structurally diverse forests [39]. Typical openness in these forests does not exceed 5%. This practice is defined in the Forestry Laws BA and SL, which prohibit larger areas of intensive logging in both countries (Forest Law, Official Gazette of the Federation of Bosnia and Herzegovina, No. 20/02 and Forest Act, Official Gazette of the Republic of Slovenia, No. 30/93).

2.2. Study Approach

The effect of harvest intensity was studied in a repeated-measures experiment. At SL, 27 research plots were studied (2011 and 2014) before and after the harvest treatments; eighteen with 50% and 100% (SL50% and SL100%) harvest intensity, representing harvested tree volume compared to total growing tree volume; and nine control plots (SL < 20%) with no treatment in 2011, representing a regularly managed forest with the harvest intensity up to 20% of the total growing tree volume. To gain a more complete insight into the effect of harvest intensity in a paired-plot experiment using an identical sampling methodology, we compared a 0% plot (unmanaged forest since 1960) and two <20% plots (managed forest in 2012). The control plot in BA 0% represented an old-growth forest without treatment. Surveys and sampling at BA were conducted in 2015. The area of each plot at the observed sites was 0.50 and 0.45 hectares, adjusted to the proposed minimum experimental area for fir-beech forests [40].

2.3. Field Work

The soil sampling design followed a number of samples suggested by the ICP Forests Forest Soil Co-ordinating Centre and the Expert Panel on Soil and Soil Solution to detect differences in treatment effects [41]. Organic (litter layer (Ol), fragmented/fermented layer (Of), humified layer (Oh)) and top mineral soil (top 0 to 10 cm (M10)) were collected at designated sampling points in each plot. The number of exact sampling points was 8 (0%) and 16 (<20%) in BA and 45 (per each < 20%, 50%, 100%) in SL (Figure 2). The distance between sampling points was at least 10 m. Three soil samples were collected at each point to create a composite sample. The Ol, Of, Oh and M10 soil layers were sampled separately. Each separate organic layer was excavated below a 25 × 25 cm frame; each mineral soil layer to a depth of 10 cm was sampled with a stainless-steel auger (diameter = 6.7 cm). For the organic layer, the determined thickness was delineated on four sides of the square frame. Due to the significant influence of the type on soil organic carbon (SOC) stocks [42], only Rendzic Leptosols and Cambisol were considered for sampling. Above-ground stoniness was estimated for each plot.
The canopy openness was analyzed using 40 digital color hemispherical photographs (HPs) at BA (at 0% 8 points and at 20% 24 points) and SL (at 100% 8 points) sites. Canopy condition is considered an indicator of change in the fir-beech forest canopy and was used to characterize the relationship between canopy openness and soil organic matter. The HPs were taken in the high growing season using a DSLR Canon EOS Rebel T3i digital camera with a calibrated fish-eye lens and analyzed using WinScanopy software (2010 pro-d). Exposure was fitted to the above-canopy conditions on each plot prior to shooting, as described by [43], with no significant noteworthy blooming effects [44]. Canopy openness (%), direct site factor (DSF%), indirect site factor (ISF%) and total site factor (TSF%) were presented using Surfer® (Golden Software, LLC, Golden, CO, USA).

2.4. Laboratory Analysis

Organic and mineral soil samples were air-dried immediately after removal from the field. Both were weighed and then ground after the removal of stones and roots, and an aliquot of the sample was analyzed for concentrations of SOC and total nitrogen (TN), by dry combustion using a CNS analyzer (Elemental Analyser LECO CNS 2000, St. Joseph, MI, USA). The volume of stones in the soil was evaluated. The ratio of C to N was calculated from their concentrations. The pH values of the soil were determined in 0.01 M CaCl2 solution. From the weight (kg m−2), the SOC stock of Ol, Of and Oh layers (Formula (1)) and from the bulk density of the fine soil (kg dm−3) the stock M10 SOC (Formula (2)) were calculated.
In addition, a possible carbonate content was determined for samples that reacted with hydrochloric acid or exceeded a pH value in CaCl2 of 6.5 on a calcimeter (Scheibler instrument, Eijkelkamp Soil & Water, Giesbeek, The Netherlands). Most samples did not contain carbonates, as soils on limestone and dolostone are rarely calcareous. Concentrations of SOC were calculated from the difference between total organic carbon and inorganic carbon concentrations. Samples with a pH less than 6.5 were assumed to have total carbon consisting entirely of organic carbon.
  • Ol/Of/Oh-layer:
SOC-stockOl/Of/Oh = SOC × weight × 10.
SOC-stockOl/Of/Oh = stock of soil organic carbon in the Ol/Of/Oh layer (g/m2/100 t/ha).
SOC = concentration of organic carbon in the organic layer (%).
Weight = dry weight of the organic layer (kg/m2).
  • M10-layer:
SOC-stock M10 = SOC × BD × d × CFst*10.
SOC-stock M10 = soil organic carbon stock in mineral soil (kg/m2 × 10 = t/ha).
d = depth class/horizon thickness (m).
SOC = concentration of organic carbon (%).
BD = bulk density (kg/dm3).
CFst = correction factor for stoniness, 100 − (% stones)/100.

2.5. Statistical Analysis

Statistical analyses included descriptive–statistical parameters and univariate statistical methods to compare soil properties between plots. Differences between treatments at the SL site for each variable (SOC, TN, C/N ratio, weight (organic layer) and SOC stocks for the observed layers (repeated measurements in 2011 and 2014)) were analyzed using a paired t-test (p < 0.05). For the BA site only, we used the t-test (p < 0.05) to compare soil properties between control and 20% harvest plots. Before performing univariate statistical analysis, the Shapiro–Wilk test was performed to check the normality of the sampling distributions. Homogeneity of variance was determined using Levene’s test. Variables that did not conform to the normal distribution were sqrt-transformed before further testing.
Relative differences in SOC, TN, C/N, weights or BD, SOC stocks were calculated between initial and after post-treatment values. A one-way test ANOVA was performed to determine significant differences (p < 0.05) between the effects of <20% BA, <20% SL, 50% SL, and 100% SL harvest intensities.
Redundancy analysis (RDA) aimed to present a general pattern of variation and relationships among plots with different canopy openness conditions and focused on determining the most distinguishing parameters among 0%, 20% and 100% harvest intensity. Response variables were log-transformed. We used the Vegan package developed in R Studio Version 1.3.1093 (RStudio Team, 2020). Pearson correlation was applied to examine the relationships between principal components and light parameters.

3. Results

3.1. Quantities and Quality of Soil Organic Matter

Based on the quantity and quality of the soil organic matter, observed by SOC and TN concentration, C/N ratio, weight/BD and SOC stocks, the BA and SL Dinaric forest sites showed high similarity of upper soil layers (Table S1). The organic layers Ol, Of and Oh layers showed high similarity, while the mineral layer M10 proved to be more variable. Moreover, the variability of soil organic matter properties increased with increasing soil depth in the Dinaric karst region.
However, despite the similarities, the SOC stocks in the upper soil layers varied in a wide range between 0.36 and 2.58 t ha−1 in Ol, 0.73 and 4.07 t ha−1 in Of, 2.54 and 7.43 t ha−1 in Oh and 25.2 and 41.9 t ha−1 in M10 (Figure 3, Table S1). These variations, starting from an old-growth forest BA0% to the close-to-nature managed forest plots BA < 20%, SL < 20%, SL < 50% and SL < 100%, affect the ability to detect the effect of harvesting.

3.2. The Effect of Increased Harvesting Intensity on the Upper Soil Layers

The results of the analysis of mean differences show the significant influence of harvesting intensity on the values SOC, TN concentrations, C/N ratio, weights or BD and C stocks (Figure 3). The significant level of differences between the initially measured variables and after treatment (50% and 100%) indicates a clear effect of increased harvesting intensity on the decrease in SOC, TN, weights and C stocks. The effect was not observed in the deeper part of the organic layer Oh or in a mineral M10 layer. This indicates that there is no short-term response in mineral soil to an increase in harvest intensity. On the contrary, no significant differences were found in the observed variables at BA after comparing the control plots with 0% (old-growth forest) and 20% harvest intensity. This means that the effect of 20% harvest intensity has negligible influence on both organic and top mineral soil layers.
The relative differences in the upper soil layer indicate the effect of harvest intensity on soil organic matter quality observed from TN concentrations and C/N ratios (Figure 4). The decreases were not only observed in the plots with high harvest intensity but also in the plots without treatment, which points to annual variation in the observed characteristics in the upper soil layer. A significant decrease is shown with respect to the individual (Ol, Of and Oh) and the sum of organic layers (∑Ol, Of, Oh) for harvest intensities of 50 and 100%.
Decreased weights and SOC stocks in the organic layer differed significantly (p < 0.01) between harvest intensity treatments in SL. Relative differences in SOC stocks ranged from −10 to −38% for 50% and from −16 to −49% for 100% harvest intensity. Increases in BD and SOC stocks of the upper mineral layer also differed significantly (p < 0.01) between harvest intensities.
As shown in Figure 4, there was a negligible difference in SOC stocks in the total organic and M10 layers at BA20%.

3.3. The Relationships between Forest Canopy and Light Conditions with the Upper Soil Layer

An increase in the canopy openness reflects light conditions and indicates thermodynamic changes in the forest stand. For this reason, we investigated the relationship between topsoil properties and canopy indices based on measurements from three BA (0% and 20%) and two SL (100%) plots.
Harvest intensity influenced different values of canopy openness (%), ranging from 6.87 to 8.40% (0%BA), 5.30 to 15.8% (20%BA) and 5.02 to 43.4% (100%SL, Kočevski Rog). These variations were also consistent with measured light parameters (ISF%, DSF% and TSF%). As shown in Table S2 in the order BA0%, BA20%, SL100%, the increases in calculated light values were followed by decreases in the upper soil layers (SOC, TN and SOCSTOCK). The quality of the organic soil matter, as indicated by a greater C/N ratio, was also influenced by the canopy openness and the light conditions.
Based on the RDA results of the upper soil layers the first two axes explained 15.2% (Ol) 15.5% (Of), 15.9% (Oh) and 16.1% (M10). The restricted proportion of variation explained by the first two axes was 99.8% (Ol), 98.0% (Of), 95.6% (Oh), and 97.7% (M10).
The properties of the upper soil layer pH values, SOC, TN, C/N ratio and SOC stocks scores are shown in Table 2. In the graphs shown in Figure 5, the closer the points in the RDA graphs, the higher the degree of similarity between them. For example, the topsoil layer in the 100% plot appeared to be different from the 0 and 20% harvest intensity plots, which appeared to be more similar to each other. The closer the distribution of the samples to the explanatory vectors in the graph, the more similar the sample is to openness (%), DSF(%), ISF(%) and TSF(%).
The arrow length of each soil property indicated its contribution to the order axes. High correlations between RDA1 and PC1 with explanatory vectors were found for both the Of and Oh layers (Table 3).

4. Discussion

4.1. Upper Soil Layer Properties in the Dinaric Karst Region

The estimated SOC stocks (Table S1) were in good agreement with previously determined values (Ol + Og + Oh: 2.8 and 5.3 t ha−1; M10: 25.8 to 52.5 t ha−1) for similar site conditions at Snežnik in Slovenia [22]. The SOC stock for the organic layer fits perfectly with the estimates of other relevant studies [45] but also with reports from calcareous alpine forest stands [24,46]. The almost equal composition of SOC and TN in the organic layers could be explained by similar stand conditions and species compositions. The data support the statement that shallow calcareous mountain soils contain the highest concentrations of SOC [18]. Minor variations in soil organic horizons were most likely influenced by variations in a mixture of herb species, microrelief and elevation in fir-beech forests in the region.

4.2. Harvest Intensity Impact

Our results confirm that a harvest intensity of up to 20%, corresponding to the single-tree harvesting method, has no effect on SOC stocks and changes in soil organic matter composition SOC, TN or the C/N ratio in the upper soil layers. In contrast, 50 and 100% harvest intensities affect short-term losses in SOC stocks in the organic layer and alter organic matter composition. The effect was not consistent with upper mineral soil.
Differences indicating losses in SOC stocks are explained by reductions in above-ground litter input and increased mineralization due to changes in abiotic factors in harvested plots after canopy removal. Higher solar radiation, more heat energy adapted from the soil and temporarily higher moisture might have caused more favorable microbial conditions for decomposition [47,48]. The decomposition process is expected to be most pronounced in labile organic horizons [49], which is confirmed by our results. Our estimation of the effect is supported by the reported soil CO2 efflux, which increased after harvest for both intensities on the same study plots [50]. We assume that soil conditions are not altered enough to prevent enzymatic activity in soil due to the protective function of forest surrounding the gap. The long-term impact is difficult to predict based on these findings. In the short term, the increased rate of organic matter turnover is a favorable outcome of harvesting that supports faster forest regeneration. Some reports from the region suggest that lower forest density, which can be associated with openness, means higher forest ecosystem stability [9], but this does not consider the gap formation. Such short-term losses in the amount of soil organic matter and accelerated litter decay have been reported in beech-dominated limestone stands of the Austrian Alps [51]. The effect of enhanced decomposition can last up to 4–5 years [52], so a longer period should be considered for repeated measurements.
Low-intensity harvesting up to 20%, commonly applied as selective, single-tree harvests, had an insignificant effect on soil organic matter according to our results from the Bosnian site on Mount Bjelašnica, which is consistent with the general understanding that no reduction [53,54,55,56] or even an increase [57] occurs after the application of low-intensity harvests. Moreover, no losses in the SOC stocks or even an increase were observed in the upper ten centimeters of the mineral soil. A meta-analysis evaluating the response of SOC stocks to timber harvesting at more than 20 forest sites showed an increase in SOC stocks in the mineral soil [17]. The applied forest practices in the Dinaric region are evaluated as SOC protective. This conclusion is supported by the practice of preserving branches and woody debris in forests [58] and protective measures considering soil disturbance [59].
The relationship between canopy openness and soil properties (Figure 5) supports the accelerated decomposition explanation. This openness may regulate soil organic matter redistribution and decomposition. Canopy openings of alpine fir forests of the Eastern Tibetan Plateau have been reported to increase litter decomposition, with the highest rates in the middle of gaps [60]. We hypothesize that one consequence, based on the character of the relationships, is the increased leaching of carbon and nitrogen into the mineral soil. This process may imply a shift to deeper horizons and stabilization of organic matter with a mineral component. Moreover, SOC can be sequestered in deep soil layers, among others [61]. We hypothesize that it is critical to find an optimal size of an area to harvest in order to preserve SOC and improve soil productivity. Impact detection is challenging due to the heterogeneity and high dispersion of mean values of SOC, which usually vary between 42% and 70% [62]. Given the importance of gap formation to soil productivity, further studies are needed to optimize canopy openness and better understand long-term impacts and possible shifts from SOC to more stable forms in mineral soils of the Dinaric Karst region. Gaps represent favorable sites for natural regeneration [31], where increased canopy openness could lead to strong changes in SOM quantity and quality.
There was a redistribution of SOC between the organic and mineral layers of the soil, indicating that the entire soil profile should be analyzed to better assess the effects of management [63]. In addition, losses of SOC from the organic horizon were accompanied by gains in the M10 layer, which can be attributed to the combined effect of the increased decomposition and transport of SOC deeper into the soil. This raises the question of the possible effect of forest thinning on the translocation of SOC from labile to more stable forms and mineral soils.

5. Conclusions

We combined soil data from a fir-beech forest in the Bosnian and Slovenian Dinaric region to investigate the influence of harvest intensity on soil organic matter properties, especially soil organic C stocks. A harvest intensity of <20% of the growing stock did not lead to a decrease in SOC stock, which is the most predictable result and indicates no change in soil and site productivity. Harvest intensities of 50 and 100% on 0.45 ha induced a SOC stock decrease in the organic layer (on average 24% to 32%) and increases in the 0–10 mineral layer, explained by changes in soil bulk density. The impact of short-term losses on soil productivity may depend on regeneration dynamics and restoration methods.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/f12050581/s1, Table S1: Soil organic matter quantity and quality in Bosnian (BA20%) and Slovenian sites (SL0%, SL50% and SL100%). Means and standard deviations of the soil organic carbon (SOC) and total N (TN) concentrations, C/N ratio, weights, bulk densities (BD) and SOC stocks are shown for the organic and 0–10 cm mineral (M10) layer.; Table S2: Means (bold) and standard deviations for openness (%), light condition parameters—direct site factor (DSF%), indirect site factor (ISF%), total site factor (TSF%) and soil organic matter properties pH values, soil organic carbon (SOC in %), total nitrogen (TN in %), C/N ratio, SOC stocks (t ha−1).

Author Contributions

Conceptualization, E.H., M.K. and P.S.; formal analysis, M.Č.; investigation, E.H., A.M. and M.F.; methodology, M.Č., M.F., D.Ž. and M.K.; project administration, P.S.; resources, D.Ž., H.Č. and P.S.; software, M.Č.; supervision, H.Č. and P.S.; validation, M.Č. and A.M.; visualization, E.H.; writing—original draft, E.H.; writing—review and editing, M.Č., A.M., M.F., D.Ž., H.Č. and P.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by LIFE09ENV/IT/000078 ManFor C.BD coordinated by the National Research Council through the Institute of Agro-environmental and Forest Biology (CNR-IBAF); COST Action FP1305 BioLink: 1305, COST-STSM-ECOST-STSM-FP1305-011016-080693, COST Action FP1206 EUMIX-FOR: FP1206 COST-STSM-FP1206-26763, LIFEGENMON (LIFE13 ENV/SI/000148), EUFORINNO (RegPot Nr. 315982), and Hojka Kraigher, Ines Strauss and Tanja Mrak. The preparation of the paper was partly supported by the Programme Groups P4-0107 and P4-0085, financed by the Slovenian Research Agency; COST Action TN1401 Training School for Successful Scientific Writing and especially to Professor Ceulemans Reinhart from the University of Antwerp Department of Biology.

Data Availability Statement

All data can be found through contacting author E.H. via e-mail ([email protected]) or Slovenian Forestry Institute in Ljubljana via e-mail (A.M. ([email protected])).

Acknowledgments

Many thanks to Hojka Kraigher and colleagues from Slovenian Forestry Institute. We also thank Alma Hajrudinović-Bogunić, Mirel Subašić, Milica Kašanin-Grubin and Senad Hukić for their constructive comments and discussions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Diaci, J.; Roženbergar, D.; Anić, I.; Mikac, S.; Saniga, M.; Kucbel, S.; Višnjić, Ć.; Ballian, D. Structural dynamics and synchronous silver fir decline in mixed old-growth mountain forests in Eastern and Southeastern Europe. Forestry 2011, 84, 479–491. [Google Scholar] [CrossRef] [Green Version]
  2. Čater, M.; Kobler, A. Light response of Fagus sylvatica L. and Abies alba Mill. in different categories of forest edge-Vertical abundance in two silvicultural systems. For. Ecol. Manag. 2017, 391, 417–426. [Google Scholar] [CrossRef]
  3. Körner, C.; Hiltbrunner, E. The 90 ways to describe plant temperature, Perspect. Plant. Ecol. 2018, 30, 16–21. [Google Scholar] [CrossRef]
  4. Schütz, J.-P.; Saniga, M.; Diaci, J.; Vrska, T. Comparing close-to-nature silviculture with processes in pristine forests: Lessons from Central Europe. Ann. For. Sci. 2016, 73, 911–921. [Google Scholar] [CrossRef] [Green Version]
  5. Schütz, J.P. Silvicultural tools to develop irregular and diverse forest structures. Forestry 2002, 75, 327–337. [Google Scholar] [CrossRef] [Green Version]
  6. Čavlović, J.; Božič, M.; Bončina, A. Stand structure of an uneven-aged fir-beech forest with an irregular diameter structure: Modelling the development of the Belevine forest, Croatia. Eur. J. For. Res. 2006, 125, 325–333. [Google Scholar] [CrossRef]
  7. Banković, S.; Medarević, M.; Pantić, D.; Petrović, N. The National Forest Inventory of the Republic of Serbia: The growing Stock of the Republic of Serbia; Ministry of Agriculture, Forestry and Water Management of the Republic of Serbia, Forest Directorate: Belgrade, Serbia, 2009; p. 238. Available online: http://www.upravazasume.gov.rs (accessed on 10 October 2019).
  8. Bončina, A. History, current status and future prospects of uneven-aged forest management in the Dinaric region: An overview. Forestry 2011, 84, 467–478. [Google Scholar] [CrossRef]
  9. Čater, M.; Levanič, T. Beech and silver fir’s response along the Balkan’s latitudinal gradient. Sci. Rep. 2019, 9, 16269. [Google Scholar] [CrossRef] [Green Version]
  10. Wang, S.; Liu, Q.; Zhang, D. Karst rocky desertification in southwestern China: Geomorphology, landuse, impact and rehabilitation. Land Degrad. Dev. 2004, 15, 115–121. [Google Scholar] [CrossRef]
  11. Seidl, R.; Rammer, W.; Lexer, M.J. Adaptation options to reduce climate change vulnerability of sustainable forest management in the Austrian Alps. Can. J. For. Res. 2011, 41, 694–706. [Google Scholar] [CrossRef] [Green Version]
  12. IPCC. Climate Change 2014: Synthesis Report. In Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Core Writing Team, Pachauri, R.K., Meyer, L.A., Eds.; IPCC: Geneva, Switzerland, 2014; p. 151. Available online: www.ipcc.ch (accessed on 3 September 2019).
  13. Samec, P.; Caha, J.; Zapletal, M.; Tuček, P.; Cudlín, P.; Kučera, M. Discrimination between acute and chronic decline of Central European forests using map algebra of the growth condition and forest biomass fuzzy sets: A case study. Sci. Total Environ. 2017, 599–600, 899–909. [Google Scholar] [CrossRef]
  14. Giardina, C.P.; Ryan, M.G. Evidence that decomposition rates of organic carbon in mineral soil do not vary with temperature. Nature 2000, 404, 858–861. [Google Scholar] [CrossRef]
  15. Schmidt, M.W.I.; Torn, M.S.; Abiven, S.; Dittmar, T.; Guggenberger, G.; Janssens, I.A.; Kleber, M.; Kögel-Knabner, I.; Lehmann, J.; Manning, D.A.C.; et al. Persistence of soil organic matter as an ecosystem property. Nature 2011, 478, 49–56. [Google Scholar] [CrossRef] [Green Version]
  16. Novak, J.; Slodicak, M.; Dusek, D. Thinning effects on forest productivity and site characteristics in stands of Pinus sylvestris in the Czech Republic. For. Syst. 2011, 20, 464–474. [Google Scholar]
  17. Johnson, D.W.; Curtis, P.S. Effects of forest management on soil C and N storage: Meta-analysis. For. Ecol. Manag. 2001, 140, 227–238. [Google Scholar] [CrossRef]
  18. Baritz, R.; Seufert, G.; Montanarella, L.; Ranst, E.V. Carbon concentrations and stocks in forest soils of Europe. For. Ecol. Manag. 2010, 260, 262–277. [Google Scholar] [CrossRef]
  19. Taylor, A.R.; Wang, J.R.; Kurz, W.A. Effects of harvesting intensity on carbon stocks in eastern Canadian red spruce (Picea rubens) forests: An exploratory analysis using the CBM-CFS3 simulation model. For. Ecol. Manag. 2008, 263, 31–38. [Google Scholar] [CrossRef]
  20. Pötzelsberger, E.; Hasenauer, H. Soil change after 50 years of converting Norway spruce dominated age class forests into single tree selection forests. For. Ecol. Manag. 2015, 338, 176–182. [Google Scholar] [CrossRef]
  21. Klein, D.; Fuentes, J.P.; Schmidt, A.; Schmidt, H.; Schulte, A. Soil organic C as affected by silvicultural and exploitative interventions in Nothofagus pumilio forests of the Chilean Patagonia. For. Ecol. Manag. 2008, 255, 3549–3555. [Google Scholar] [CrossRef]
  22. Urbančič, M.; Kobal, M.; Vilhar, U.; Simončič, P. Zaloge organske snovi v izbranih sestojih na bukovih rastiščih. In Trajnostna Raba Lesa v Kontekstu Sonaravnega Gospodarjenja z Gozdovi; Humar, M., Kraigher, H., Eds.; Silva Slovenica; Gozdarski Inštitut Slovenije: Ljubljana, Slovenia, 2009; pp. 19–29. [Google Scholar]
  23. Ritter, E.; Vesterdal, L. Gap formation in Danish beech (Fagus sylvatica) forests of low management intensity: Soil moisture and nitrate in soil solution. Eur. J. For. Res. 2006, 125, 139–150. [Google Scholar] [CrossRef]
  24. Christophel, D.; Spengler, S.; Schmidt, B.; Ewald, J.; Prietzel, J. Customary selective harvesting has considerably decreased organic carbon and nitrogen stocks in forest soils of the Bavarian Limestone Alps. For. Ecol. Manag. 2014, 305, 167–176. [Google Scholar] [CrossRef]
  25. Christophel, D.; Höllerl, S.; Prietzel, J.; Steffens, M. Long-term development of soil organic carbon and nitrogen stocks after shelterwood- and clear-cutting in a mountain forest in the Bavarian Limestone Alps. Eur. J. For. Res. 2015, 134, 623–640. [Google Scholar] [CrossRef]
  26. Mayer, M.; Prescott, C.E.; Abaker, W.E.A.; Augusto, L.; Cécillon, L.; Ferreira, G.W.D.; James, J.; Jandl, R.; Katzensteiner, K.; Laclau, J.P.; et al. Tamm Review: Influence of forest management activities on soil organic carbon stocks: A knowledge synthesis. For. Ecol. Manag. 2020, 466, 118–127. [Google Scholar] [CrossRef]
  27. Ewald, J. Ecological background of crown condition, growth and nutritional status of Picea abies (L.) Karst. in the Bavarian Alps. Eur. J. For. Res. 2005, 124, 9–18. [Google Scholar] [CrossRef]
  28. Kobal, M.; Grčman, H.; Zupan, M.; Levanič, T.; Simončič, P.; Kadunc, A.; Hladnik, D. Influence of soil properties on silver fir (Abies alba Mill.) growth in the Dinaric Mountains. For. Ecol. Manag. 2015, 337, 77–87. [Google Scholar] [CrossRef] [Green Version]
  29. Reger, B.; Gottlein, A.; Katzensteiner, K.; Ewald, J. Assessing the Sensitivity of Mountain Forests to Site Degradation in the Northern Limestone Alps, Europe. Mt. Res. Dev. 2015, 35, 139–151. [Google Scholar] [CrossRef]
  30. Canham, C.D.; Finzi, A.C.; Pacala, S.W.; Burbank, D.H. Causes and consequences of resource heterogeneity in forests: Interspecific variation in light transmission by canopy trees. Can. J. For. Res. 1994, 24, 337–349. [Google Scholar] [CrossRef]
  31. Čater, M.; Diaci, J.; Roženbergar, D. Gap size and position influence variable response of Fagus sylvatica L. and Abies alba Mill. For. Ecol. Manag. 2014, 25, 128–135. [Google Scholar] [CrossRef]
  32. Hirabuki, Y. Heterogeneous dispersal of tree litterfall corresponding with patchy canopy structure in a temperate mixed forest. Vegetatio 1991, 94, 69–79. [Google Scholar] [CrossRef]
  33. Tingey, D.T.; Phillips, D.L.; Johnson, M.G.; Rygiewicz, P.T.; Beedlow, P.A.; Hogsett, W.E. Estimates of Douglas-fir fine root production and mortality from minirhizotrons. For. Ecol. Manag. 2005, 204, 359–370. [Google Scholar] [CrossRef]
  34. Gallardo, A. Effect of tree canopy on the spatial distribution of soil nutrients in a Mediterranean Dehesa. Pedobiologia 2003, 47, 117–125. [Google Scholar] [CrossRef]
  35. Sebastia, M.T.; Casals, P.; Vojniković, S.; Bogunić, F.; Beus, V. Plant diversity and soil properties in pristine and managed stands from Bosnian mixed forests. Forestry 2005, 78, 297–303. [Google Scholar] [CrossRef]
  36. Kutnar, L.; Eler, K.; Marinšek, A. Effects of different silvicultural measures on plant diversity-the case of the Illyrian Fagus sylvatica habitat type (Natura 2000). iForest 2015, 9, 318–324. [Google Scholar] [CrossRef] [Green Version]
  37. De Cinti, B.; Bombi, P.; Ferretti, F.; Cantiani, P.; Di Salvatore, U.; Simončič, P.; Kutnar, L.; Čater, M.; Garfì, V.; Mason, F.; et al. From the experience of LIFE+ ManFor C.BD to the manual of best practices in sustainable forest management. Ital. J. Agron. 2016, 11, 1–175. [Google Scholar] [CrossRef] [Green Version]
  38. Topalović, M.; Travar, J. Tumač Pedološke Karte i Tipološke Karte Gospodarskih Jedinica “Igman”, “Zujevina” i “Mehina Luka” (Interpreter of the Pedological Map and Typological Map of Economic Units “Igman”, “Zujevina” and “Mehina Luka”); Institut za Šumarstvo i Drvnu Industriju, Zavod za Šumarstvo i Lovstvo (Institute for Forestry and Wood Industry, Department of Forestry and Hunting): Beograd, Serbia, 1978. (In Serbian) [Google Scholar]
  39. Čater, M.; Diaci, J. Scientific support for close-to-nature forestry. In Forests and forestry in Slovenia, Studia Forestalia Slovenica, 1st ed.; Čater, M., Železnik, P., Eds.; Slovenian Forestry Institute, The Silva Slovenica Publishing Centre: Ljubljana, Slovenia, 2020. [Google Scholar]
  40. Kotar, M. Zgradba, Rast in Donos Gozda na Ekoloških in Fizioloških Osnovah; Zveza Gozdarskih Društev Slovenije in Zavod za Gozdove Slovenije: Ljubljana, Slovenia, 2005; p. 500. (In Slovenian) [Google Scholar]
  41. Cools, N.; De Vos, B. Part X: Sampling and Analysis of Soil. In Manual on Methods and Criteria for Harmonized Sampling, Assessment, Monitoring and Analysis of the Effects of Air Pollution on Forests; UNECE ICP Forests Programme Coordinating Centre, Ed.; Thünen Institute of Forest Ecosystems: Eberswalde, Germany, 2016; p. 29. Available online: http://www.icp-forests.org/manual.htm (accessed on 20 October 2020).
  42. Grüneberg, E.; Schöning, I.; Kalko, E.K.V.; Weisser, W.W. Regional organic carbon stock variability: A comparison between depth increments and soil horizons. Geoderma 2010, 155, 426–433. [Google Scholar] [CrossRef]
  43. Zhang, Y.; Chen, J.M.; Miller, J.R. Determining digital hemispherical photograph exposure for leaf area index estimation. Agric. For. Meteorol. 2005, 133, 166–181. [Google Scholar] [CrossRef]
  44. Leblanc, S.G.; Chen, J.M.; Fernandes, R.; Deering, D.W.; Conley, A. Methodology comparison for canopy structure parameters extraction from digital hemispherical photography in boreal forests. Agric. For. Meteorol. 2005, 129, 187–207. [Google Scholar] [CrossRef] [Green Version]
  45. Haber, J. Beziehungen Zwischen dem Humusgehalt, Humusvorrat Bayerischer Bergwaldböden und Steuernden Standorts- und Bestandesfaktoren. Doctoral Dissertation, University of Bayreuth, Bayreuth, Germany, 1985. [in Christophel, D.; Spengler, S.; Schmidt, B.; Ewald, J.; Prietzel, J. Customary selective harvesting has considerably decreased organic carbon and nitrogen stocks in forest soils of the Bavarian Limestone Alps. For. Ecol. Manag. 2014, 305, 167–176, doi:10.1016/j.foreco.2013.05.054]. [Google Scholar]
  46. Prietzel, J.; Christophel, D. Organic carbon stocks in forest soils of the German Alps. Geoderma 2014, 221–222, 28–39. [Google Scholar] [CrossRef]
  47. Mayer, M.; Sandén, H.; Rewald, B.; Godbold, D.L.; Katzensteiner, K. Increase in heterotrophic soil respiration by temperature drives decline in soil organic carbon stocks after forest windthrow in a mountainous ecosystem. Funct. Ecol. 2017, 31, 1163–1172. [Google Scholar] [CrossRef]
  48. Morehouse, K.; Johns, T.; Kaye, J.; Kaye, M. Carbon and nitrogen cycling immediately following bark beetle outbreaks in southwestern ponderosa pine forests. For. Ecol. Manag. 2008, 255, 2698–2708. [Google Scholar] [CrossRef]
  49. Misson, L.; Tang, J.; Xu, M.; Mckay, M.; Goldstein, A. Influences of recovery from clear-cut, climate variability, and thinning on the carbon balance of a young ponderosa pine plantation. Agric. For. Meteorol. 2005, 130, 207–222. [Google Scholar] [CrossRef]
  50. Čater, M.; Darenova, E.; Simončič, P. Harvesting intensity and tree species affect soil respiration in uneven-aged Dinaric forest stands. For. Ecol. Manag. 2021, 480, 118638. [Google Scholar] [CrossRef]
  51. Katzensteiner, K. Effects of harvesting on nutrient leaching in a Norway spruce (Picea abies Karst.) ecosystem on a Lithic Leptosol in the Northern Limestone Alps. Plant Soil 2003, 250, 59–73. [Google Scholar] [CrossRef]
  52. Prescott, C.E.; Hope, G.D.; Blevins, L.L. Effect of gap size on litter decomposition and soil nitrate concentrations in a high-elevation spruce-fir forest. Can. J. For. Res. 2003, 3, 2210–2220. [Google Scholar] [CrossRef]
  53. Mund, M.; Schulze, E.D. Impacts of forest management on the carbon budget of European beech (Fagus sylvatica) forests. Allg. For. Jagdztg. 2006, 177, 47–63. [Google Scholar]
  54. Nave, L.E.; Vance, E.D.; Swanston, C.W.; Curtis, P.S. Harvest impacts on soil carbon storage in temperate forests. For. Ecol. Manag. 2010, 259, 857–866. [Google Scholar] [CrossRef]
  55. Nilsen, P.; Strand, L.T. Thinning intensity effects on carbon and nitrogen stores and fluxes in a Norway spruce stand (Picea abies (L.) Karst.) stand after 33 years. For. Ecol. Manag. 2008, 256, 201–208. [Google Scholar] [CrossRef]
  56. Clarke, N.; Gundersen, P.; Jönsson-Belyazid, U.; Kjønaas, O.J.; Persson, T.; Sigurdsson, B.D.; Stupak, I.; Verstedal, L. Influence of different tree-harvesting intensities on forest soil carbon stocks in boreal and northern temperate forest ecosystems. For. Ecol. Manag. 2015, 351, 9–19. [Google Scholar] [CrossRef]
  57. Kishchuk, B.E.; Morris, D.M.; Lorente, M.; Keddy, T.; Sidders, D.; Quideau, S.; Maynard, D. Disturbance intensity and dominant cover type influence rate of boreal soil carbon change: A Canadian multi-regional analysis. For. Ecol. Manag. 2016, 381, 48–62. [Google Scholar] [CrossRef]
  58. Jandl, R.; Lindner, M.; Vesterdal, L.; Bauwens, B.; Baritz, R.; Hagedron, F.; Johnson, D.W.; Minkkinen, K.; Byrne, K.A. How strongly can forest management influence soil carbon sequestration? Geoderma 2007, 137, 253–268. [Google Scholar] [CrossRef]
  59. Achat, D.L.; Deleuze, C.; Landmann, G.; Pousse, N.; Ranger, J.; Augusto, L. Quantifying consequences of removing harvesting residues on forest soils and tree growth—A meta-analysis. For. Ecol. Manag. 2015, 348, 124–141. [Google Scholar] [CrossRef]
  60. Wu, Q.; Wu, F.; Yang, W.; Zhao, Y.; He, W.; Tan, B. Foliar litter nitrogen dynamics as affected by forest gap in the alpine forest of Eastern Tibet Plateau. PLoS ONE 2014, 9, e97112. [Google Scholar] [CrossRef]
  61. Hoffmann, T.; Schlummer, M.; Notebaert, B.; Verstraeten, G.; Korup, O. Carbon burial in soil sediments from Holocene agricultural erosion, Central Europe. Glob. Biogeochem. Cycles 2013, 27, 828–835. [Google Scholar] [CrossRef]
  62. Wilding, L.P.; Drees, L.R.; Nordt, L.C. Spatial variability: Enhancing the mean estimate of organic and inorganic carbon in a sampling unit. In Assessment Methods for Soil Carbon; Lal, R., Kimble, J.M., Follett, R.F., Stewart, B.A., Eds.; Lewis Publishers: Boca Raton, FL, USA, 2001; pp. 69–86. [Google Scholar]
  63. Yanai, R.D.; Currie, W.S.; Goodale, C.L. Soil carbon dynamics after forest harvest: An Ecosystem Paradigm Reconsidered. Ecosystems 2003, 6, 197–212. [Google Scholar] [CrossRef]
Figure 1. Location of the sites in Bosnia and Herzegovina (BA, Bjelašnica Mt.—B) and Slovenia (SL, Kočevski Rog—K, Snežnik—S and Trnovski gozd—T).
Figure 1. Location of the sites in Bosnia and Herzegovina (BA, Bjelašnica Mt.—B) and Slovenia (SL, Kočevski Rog—K, Snežnik—S and Trnovski gozd—T).
Forests 12 00581 g001
Figure 2. Scheme of the experimental plots (BA0% 1; BA < 20% 2; SL < 20% 9, SL50% 9, SL100% 9) showing the total number of sampled points. The study approach combines paired-plot comparison (Bosnia—BA) and repeated measurements approach (Slovenia—SL). White color represents the forest canopy and black the canopy openings. Ol (litter), Of (fragmented/fermented), Oh (humified) and M10 (top mineral 0 to 10 cm) are sampled upper soil layers and N represents the number of sampling points per plot.
Figure 2. Scheme of the experimental plots (BA0% 1; BA < 20% 2; SL < 20% 9, SL50% 9, SL100% 9) showing the total number of sampled points. The study approach combines paired-plot comparison (Bosnia—BA) and repeated measurements approach (Slovenia—SL). White color represents the forest canopy and black the canopy openings. Ol (litter), Of (fragmented/fermented), Oh (humified) and M10 (top mineral 0 to 10 cm) are sampled upper soil layers and N represents the number of sampling points per plot.
Forests 12 00581 g002
Figure 3. Means and standard deviations of soil organic carbon (SOC; %), total nitrogen (TN; %), C/N ratio, weights or bulk densities (kg m−2, kg dm−3) and carbon stocks (SOC stock; t ha−1) distributed in the organic (Ol, Of, Oh) and mineral (M10) layers at the BA (Bosnian) and SL (Slovenian) sites. The green and violet dashed areas represent connected external mean values between sites. The labels 0, <20, 50 and 100% represent harvest intensity. * defines significant (p < 0.05) difference between initially measured values and after the treatment.
Figure 3. Means and standard deviations of soil organic carbon (SOC; %), total nitrogen (TN; %), C/N ratio, weights or bulk densities (kg m−2, kg dm−3) and carbon stocks (SOC stock; t ha−1) distributed in the organic (Ol, Of, Oh) and mineral (M10) layers at the BA (Bosnian) and SL (Slovenian) sites. The green and violet dashed areas represent connected external mean values between sites. The labels 0, <20, 50 and 100% represent harvest intensity. * defines significant (p < 0.05) difference between initially measured values and after the treatment.
Forests 12 00581 g003
Figure 4. Mean relative differences and standard errors of soil organic carbon (SOC), total nitrogen (TN), C/N ratio, weight or bulk density (BD) and SOC stocks in Ol (litter), Of (fragmented/fermented) and Oh (humified) layers calculated between control plots versus harvested plots (harvest intensity <20% in Bosnian site: <20% BA, harvest intensity <20% in Slovenian site with no treatment: <20% SL, harvest intensity 50% and 100% in Slovenian site: 50% SL and 100% SL). The green and violet dashed areas connect external mean values.
Figure 4. Mean relative differences and standard errors of soil organic carbon (SOC), total nitrogen (TN), C/N ratio, weight or bulk density (BD) and SOC stocks in Ol (litter), Of (fragmented/fermented) and Oh (humified) layers calculated between control plots versus harvested plots (harvest intensity <20% in Bosnian site: <20% BA, harvest intensity <20% in Slovenian site with no treatment: <20% SL, harvest intensity 50% and 100% in Slovenian site: 50% SL and 100% SL). The green and violet dashed areas connect external mean values.
Forests 12 00581 g004
Figure 5. Redundancy analysis (RDA) for upper soil layers (organic litter (Ol), fragmented/fermented (Of), humified Oh and mineral 0–10 cm (M10) layers). Response variables are soil organic carbon (SOC), total nitrogen (TN), C/N ratio, SOC STOCK and pH values and explanatory variables are openness (%), direct site factor (DSF, %), indirect site factor (ISF, %) and total site factor (TSF, %).
Figure 5. Redundancy analysis (RDA) for upper soil layers (organic litter (Ol), fragmented/fermented (Of), humified Oh and mineral 0–10 cm (M10) layers). Response variables are soil organic carbon (SOC), total nitrogen (TN), C/N ratio, SOC STOCK and pH values and explanatory variables are openness (%), direct site factor (DSF, %), indirect site factor (ISF, %) and total site factor (TSF, %).
Forests 12 00581 g005
Table 1. Dinaric karst fir-beech forest site characteristics.
Table 1. Dinaric karst fir-beech forest site characteristics.
Experimental SiteBA 1SL(K) 2SL(S) 2SL(T) 2
Latitude (°N)43.73845.66845.67245.989
Longitude (°E)18.25415.03314.4613.759
Altitude (m a.s.l.)1420871870814
MAT (°C)7.69.08.411.3
MAP (mm)1192146515731619
Growing stock (m3 ha−1)270352442292
Increment (m3 ha−1 yr−1)6.29.48.36.2
Surface rockiness (%)19 ± 1038 ± 1331 ± 1817 ± 11
Soil pHCaCl25.2 ± 0.85.3 ± 0.65.3 ± 0.74.6 ± 0.4
1 BA: Bosnia and Herzegovina Mt Bjelašnica—B; 2 SL: Slovenia sites Kočevski Rog—K; Snežnik—S; Trnovski gozd—T; a.s.l. = above sea level; MAT = mean annual temperature; MAP = mean annual precipitation. Surface rockiness and pHCaCl2 were determined on site.
Table 2. Sample scores for upper soil layer and variables pH values, soil organic carbon (SOC), total nitrogen (TN), C/N ratio and SOC stocks.
Table 2. Sample scores for upper soil layer and variables pH values, soil organic carbon (SOC), total nitrogen (TN), C/N ratio and SOC stocks.
LayerScorepH ValueSOCTNC/N RatioSOC Stock
OlRDA10.0030.0250.171−0.2880.543
RDA2−0.0030.0050.087−0.119−0.091
PC10.0050.008−0.0220.0471.416
PC2−0.004−0.010−0.3160.572−0.024
OfRDA1−0.0570.0620.096−0.1040.649
RDA20.0120.1140.183−0.202−0.069
PC10.016−0.060−0.021−0.0281.505
PC2−0.0910.2060.406−0.4920.006
OhRDA1−0.0280.3130.265−0.1170.444
RDA20.0050.1660.211−0.213−0.299
PC1−0.0280.3130.265−0.1170.444
PC20.0050.1660.211−0.213−0.299
M10RDA1−0.079−0.675−0.3220.268−0.153
RDA2−0.0810.172−0.0730.319−0.007
PC10.0981.1030.3780.176−0.015
PC20.028−0.013−0.1140.228−0.974
Note: Labels define organic litter (Ol), fragmented/fermented (Of), humified Oh and mineral 0–10 cm (M10) soil layers. The RDA1 and RDA2 are first and second redundancy analysis axes. The PC1 and PC2 are first and second principle components.
Table 3. Correlation coefficient (R2) values between vector values, first axe of redundancy analysis (RDA1) and first principle component (PC1) for the analyzed upper soil layers.
Table 3. Correlation coefficient (R2) values between vector values, first axe of redundancy analysis (RDA1) and first principle component (PC1) for the analyzed upper soil layers.
Layer Openness (%)DSF (%)ISF (%)TSF (%)
OlRDA1−0.477−0.477−0.477−0.477
PC10.8790.8790.8790.879
R20.410.570.240.55
OfRDA1−0.418−0.418−0.418−0.418
PC10.9080.9080.9080.908
R20.150.290.020.25
OhRDA1−0.470−0.470−0.470−0.470
PC1−0.883−0.883−0.883−0.883
R20.200.350.080.32
M10RDA10.6250.6250.6250.625
PC10.7800.7800.7800.780
R20.610.680.440.70
Note: Labels define organic litter (Ol), fragmented/fermented (Of), humified Oh and mineral 0–10 cm (M10) soil layers. Bold and underlined values represent significant (p < 0.05) correlation coefficients.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hukić, E.; Čater, M.; Marinšek, A.; Ferlan, M.; Kobal, M.; Žlindra, D.; Čustović, H.; Simončič, P. Short-Term Impacts of Harvesting Intensity on the Upper Soil Layers in High Karst Dinaric Fir-Beech Forests. Forests 2021, 12, 581. https://doi.org/10.3390/f12050581

AMA Style

Hukić E, Čater M, Marinšek A, Ferlan M, Kobal M, Žlindra D, Čustović H, Simončič P. Short-Term Impacts of Harvesting Intensity on the Upper Soil Layers in High Karst Dinaric Fir-Beech Forests. Forests. 2021; 12(5):581. https://doi.org/10.3390/f12050581

Chicago/Turabian Style

Hukić, Emira, Matjaž Čater, Aleksander Marinšek, Mitja Ferlan, Milan Kobal, Daniel Žlindra, Hamid Čustović, and Primož Simončič. 2021. "Short-Term Impacts of Harvesting Intensity on the Upper Soil Layers in High Karst Dinaric Fir-Beech Forests" Forests 12, no. 5: 581. https://doi.org/10.3390/f12050581

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