Evaluation of Leaf Litter Mulching and Incorporation on Skid Trails for the Recovery of Soil Physico-Chemical and Biological Properties of Mixed Broadleaved Forests

Engineering applications can be used to mitigate the adverse effects of soil compaction and amend compacted soils. Previous literature has highlighted the beneficial effects of interventions such as litter mulching and incorporation on skid trails. However, little is known about the effectiveness of these alternatives in restoring forest soil quality after forest logging. The objective of this study was to properly elucidate the effects of the above mentioned soil protection methods, litter incorporation before skidding (LI) and litter mulching after skidding (LM), on the recovery of compacted soil’s physico-chemical and biological properties on skid trails over a 2-year period in the Hyrcanian forests of Iran to identify the best option for restoration intervention. The litter used in both methods consisted of dried leaves of the hornbeam and maple tree in three intensities of 3, 6, and 9 Mg ha−1. The results showed that the application of both methods (LI and LM) significantly improved the soil properties when compared to the untreated skid trail. Results showed that the recovery values of soil properties in the LI treatments were significantly higher than those of the LM. The recovery values of soil properties by 6 and 9 Mg ha−1 were significantly higher than those of 3 Mg ha−1, while the differences were not significant between 6 and 9 Mg ha−1. Our findings showed that soil properties were partially recovered (70–80%) over a 2-year period from treatment, compared to untreated, but the full recovery of soil properties required more time to return to the pre-harvest value. Overall, the results of this study demonstrated that the application of soil protection methods accelerates the process of recovering soil properties much faster than natural soil recovery, which can take more than 20 years in these forests.


Introduction
The increase in mechanized forest harvesting and the expansion of wood production areas raises questions about the environmental impact of machine traffic on forest soil and their subsequent impact on stand productivity [1][2][3]. Given the long-term loss of stand fertility due to the construction of skid trails, it is essential to evaluate the factors which could play a key role in achieving sustainable forest growth [4]. Scientific research has already highlighted that the passage of skidding machines can reduce the yield of forest stands as a consequence of soil compaction [5]. stability, as well as enhance infiltration rate, increase water storage capacity, and increase the levels of organic matter content [3,38]. Mulch is any organic and inorganic materials (i.e., litter, straw, leaves, plastic film, and gravel) spread on the surface to protect the soil, as opposed to materials that are incorporated into the soil profile [39]. Jourgholami et al. [3] concluded that 1.31 kgm −2 of leaf litter was an optimal level for spreading on skid trails to modify and enhance soil quality. Likewise, Jourgholami et al. [23] reported that the application of sawdust mulch resulted in a partial recovery of soil bulk density, total porosity, penetration resistance, and rut depth 6 years after logging interventions. Furthermore, Fernández et al. [14] found that the application of mulch on the surface of compaction-induced soil on skid trails resulted in an improvement in soil structure.
Similarly, Jourgholami et al. [17] concluded that the recovery of soil physical and chemical properties (except for soil C/N ratio) was significantly higher three years after logging in compacted soil treated with litter and sawdust as compared to untreated soil.
Furthermore, litter incorporation into soil showed interesting features in restoring the characteristics of compacted soil after logging operations [40]. However, the application rate and incorporation depth are the main drivers for effectiveness of organic mulch to improve soil properties [41].
Surface applied plant residues decreased soil temperature and increased soil moisture content more than the incorporation of plant residues into soil [42]. However, mulching the organic materials on the soil surface increased the possible leaching of nitrate compared to the incorporation of plant residues [43]. In an urban area, Cogger et al. [44] found that the soil bulk density was significantly higher when compost was used on the surface than the incorporated compost treatment over a 6-year period after compost application.
Similarly, Kranz et al. [41] concluded that compost incorporation significantly decreased soil bulk density. Previous studies highlighted that the decomposition rates of the residues incorporated into soil are faster than the ones of surface applied organic residues [42].
Notwithstanding the attention which scientific literature has put on both the assessment of soil impacts related to forest operations and the related possible techniques of restoration, there is still a lack of knowledge regarding the effectiveness of two known restoration methodologies (litter mulching and litter incorporation before harvesting operations) when applied to logging in mixed broadleaved stands. Comparative studies between these two different approaches are indeed scarce in the current literature. Hence, the objectives of this study were to: (1) elucidate the effects of mulching and incorporation of litter on the recovery of compacted soil properties in the skid trails over a 2-year period, compared to the untreated skid trail (UNT) and undisturbed (control) area (UND), and (2) assess the efficacy of three litter application rates of 3, 6, and 9 Mg ha −1 on the recovery process of physico-chemical, biological, and microbial soil properties.

Site Description
The study area is located within compartment n. 315 in the Gorazbon district of the Kheyrud forest research station (51 • 37 01"-51 • 37 17" E, 36 • 33 34"-36 • 33 15" N) in the Hyrcanian forests of northern Iran (Figure 1a). The altitude of the study area ranges from 1150 to 1250 m asl with a southern aspect. The study area has a humid climate with a total of 1260 mm annual rainfall. The maximum and minimum monthly rainfall occurred in October and July, respectively. The mean annual temperature is 12.8 • C and the hottest and coldest months are in July and January, respectively. According to the United States Department of Agriculture (USDA) soil taxonomy, the soil type in the study area is deep brown (Alfisols), derived from limestone and dolomite limestone, belonging to the upper Jurassic and lower Cretaceous periods. The study area is dominated by natural unevenaged stands of beech (Fagus orientalis Lipsky), with hornbeam (Carpinus betulus L.) and velvet maple (Acer velutinum Boiss). The average stem density and average growing stock of stands in the study were 227.7 tree ha −1 and 462.7 m 3 ha −1 , respectively. The average Land 2021, 10, 625 4 of 15 tree dbh (diameter at breast height) and height were 48 cm and 28 m, respectively. The single tree selection system was the main silvicultural treatment in the study area. In March 2017, a total of 98 trees with 676 m 3 were felled and processed by chainsaw, and extracted by a wheeled skidder of TAF E655 from stump to landings in July 2017. The main technical characteristics of a TAF E655 skidder are: articulated, four-wheel-drive vehicle weighing 6.8 metric tonnes, engine power of 65 hp (48 kW), a tire size of 18.4-26 inches inflated to 659 kPa on both front and rear axles, and an overall width of 2.85 m. The slope gradient of skid trails was in the range 8-15%, the average load volume was 3.4 m 3 , and operating trails had a width of 3.5 m.
Land 2021, 10, x FOR PEER REVIEW aged stands of beech (Fagus orientalis Lipsky), with hornbeam (Carpinus betulus L velvet maple (Acer velutinum Boiss). The average stem density and average growing of stands in the study were 227.7 tree ha −1 and 462.7 m 3 ha −1 , respectively. The averag dbh (diameter at breast height) and height were 48 cm and 28 m, respectively. The tree selection system was the main silvicultural treatment in the study area. In March a total of 98 trees with 676 m 3 were felled and processed by chainsaw, and extracte wheeled skidder of TAF E655 from stump to landings in July 2017. The main tec characteristics of a TAF E655 skidder are: articulated, four-wheel-drive vehicle wei 6.8 metric tonnes, engine power of 65 hp (48 kW), a tire size of 18.4-26 inches infla 659 kPa on both front and rear axles, and an overall width of 2.85 m. The slope gr of skid trails was in the range 8-15%, the average load volume was 3.4 m 3 , and ope trails had a width of 3.5 m.

Experimental Design
In total, 20 pure stands of hornbeam and maple in the study area were selecte 10 sample plots of 2 × 2 m were installed to collect the leaf litter in each stand. Th decomposed leaf litters were collected by hand in the fall, transported to the labor air-dried [45], re-weighed in each sample plot, and stored over winter before appli the following year. The dried leaf litter of the hornbeam and maple were manuall

Experimental Design
In total, 20 pure stands of hornbeam and maple in the study area were selected and 10 sample plots of 2 × 2 m were installed to collect the leaf litter in each stand. The undecomposed leaf litters were collected by hand in the fall, transported to the laboratory, air-dried [45], re-weighed in each sample plot, and stored over winter before application the following year. The dried leaf litter of the hornbeam and maple were manually and evenly distributed on the skid trail segments with three application rates of 3, 6, and 9 Mg ha −1 before and after skidding operations. In order to compare the effects of litter mulching on the surface soil (after skidding operations) and litter incorporation into the soil (before the skidding operations) on soil properties, three segments with dimensions of 4 × 5 m were randomly selected on the skid trails along the slope gradient, ranging from 10 to 14%, on the skid trails exposed to the equivalent of 10 machine passes in 2017. For litter incorporation into soil before the machine entrance over the skid trails, mixed litter of hornbeam and maple (combined as a weight ratio of 1:1) were evenly scattered over the selected segments with three application amounts of 3 (LI3), 6 (LI6), and 9 (LI9) Mg ha −1 for appropriate mixing of the litter with soil before the skidder traffic of 10 machine passes ( Figure 1). For litter mulch treatments on soil surface, litter mulches were manually distributed on the surface soil at the determined segments with three applications in the amounts of 3 (LM3), 6 (LM6), and 9 (LM9) Mg ha −1 immediately after 10 machine passes. Wire mesh was installed on the surface of the skid trail segments to capture autumn leaf fall.
In 2019, two years after incorporation on the treatment areas (LI and LM), 5 soil samples were collected in each treated segment on the mineral soil at the depth of 0-10 cm. In addition, 5 soil samples were collected in the untreated skid trails (UNT) and in the undisturbed (control) area 20 m away from the skid trails (UND). A total of 120 soil samples (i.e., 3 skid trails segments × 2 type of litter applying × 3 application rate × 5 soil samples + 3 skid trails segments × 2 references (UNT and UND) × 5 soil samples) were collected and then analyzed.

Data Collection and Laboratory Analysis
A thin-walled steel cylinder with a length of 100 mm long and diameter of 56 mm was used to take soil sample cores from the top mineral soil (depth of 10 cm) in each sample point. Then, soil cores were weighed, stored in bags, labelled, transported to the lab, and oven-dried to constant mass at 105 • C for 24 h to measure the soil bulk density and water content [46]. Soil samples were collected from a 20 × 20 × 10 cm volume, air-dried, and passed through a 2 mm sieve to measure the other soil physico-chemical properties. The soil particle size distribution was determined by the hydrometer method [47]. The water desorption method [48] was employed to determine the macroporosity. The analog handheld soil penetrometer (Eijkelkamp 06.01.SA penetrometer with a 60 • cone and a 1 m maximum measuring depth) was used to determine the soil penetration resistance (PR). The wet sieving procedure was applied [49] to determine the aggregate stability. The ASTM D854-00 2000 standard was used to determine the soil particle density, and then Equation 1 was used to calculate total porosity as follows [46]: where TP is the apparent total porosity (%), Ms is mass of soil (g), VC is the volume of the intact soil cores (246.30 cm 3 ), and 2.65 (g cm −3 ) is the particle density. The sieved soil was analyzed to determine soil chemical properties. Soil pH was determined by the Orion Ionalyzer (Model 901) pH meter in a soil: water ratio of 1:2.5. The Walkley-Black technique was employed to determine the soil organic C [50]. To determine soil total N [51], the Kjeldahl method was used. Soil C and N storage, at a depth of 0-10 cm, was determined by Equation (2): where the SO (C or N) s indicates the organic C or N storage in the soil (Mg ha −1 ); C or N is the organic C or N content (g kg −1 ); BD is the bulk density (g cm −3 ); e is the thickness of the layers (cm), and 0.1 is a depth conversion factor. The Olsen method, with a spectrophotometer, was employed to measure the available phosphorous (P). The available potassium (K), calcium (Ca), and magnesium (Mg) (by ammonium acetate extraction at pH 9) were measured by using an atomic absorption spectrophotometer [52]. Manual sampling was employed at the soil surface of 25 × 25 cm with a depth of 0-10 cm to determine the number and density of earthworms [52]. The earthworms were counted in the sampling area, transported to the laboratory, washed, euthanized by placement into hot water, oven-dried at 60 • C for 24 h to reach the constant mass, and then reweighed to determine the earthworm dry mass [53]. To measure fine root biomass, fine roots (<2 mm diameter) were separated from the soil in each sample area (20 × 20 × 10 cm) and dried at 70 • C to reach a constant mass [54]. Soil microbial respiration was measured by measuring the CO 2 evolved in a 3-day incubation experiment at 25 • C [51].

Statistical Analysis
A complete block design was randomly employed to show the effects of litter (mixed hornbeam and maple) mulching on the surface soil and incorporation into the soil with three application amounts (3, 6, and 9 Mg ha −1 ) on the skid trail to determine soil physicochemical, biological, and microbial properties as compared to the untreated skid trail (UNT) and the undisturbed (control) treatments (UND) over a 2-year period after treating and skidding operations. The Kolmogorov-Smirnov test was used to test the normality of data distribution. The Levene test was used to test the homogeneity of variance. To compare physico-chemical, biological, and microbial properties among treatments (UND, litter mulching on surface soil (i.e., LM3, LM6, and LM9), litter incorporation into soil (i.e., LI3, LI6, and LI9), and UNT), a one-way analysis of variance (ANOVA) was employed. The Tukey's test was applied as a post hoc for the treatment group means at p ≤0.05, after the ANOVA detected significant differences among treatments. The Pearson correlation was employed to test the relationship between physico-chemical, biological, and microbial properties. All statistical analyses were performed using the SPSS (release 17.0; New York, NY, USA) statistical package.

Soil Physico-Chemical Properties
The results of the ANOVA showed that both methods (litter incorporation during skidding and litter mulching after skidding) had a significant effect on all physical (Table 1), chemical (Table 2), and biological (Table 3) properties of the soil. The highest soil bulk density, penetration resistance, and amount of sand was found in the UNT, followed by LM3, LM6, and LM9, whereas the lowest soil bulk density, penetration resistance, and amount of sand was detected in the UND. Total porosity, macroporosity, and aggregate stability were significantly highest in the UND, followed by LI9 and LI6, whereas the lowest total porosity, macroporosity, and aggregate stability was measured at the UNT. The UND had the highest soil moisture, followed by LM9 and LM6, whereas the lowest soil moisture was found in the UNT. Silt under LM3 was lower as compared to other treatments, although the value was higher than the silt measured at the UNT. The highest clay content was observed on the LM9 followed by LI9, whereas the lowest clay was detected on UNT ( Table 1).
The highest soil pH was found under UNT, followed by the LM3 ≈ LM6 treatments, whereas the lowest pH was observed at the UND. The highest values of soil organic C, soil C/N ratio, and C storage were detected in UND followed by LM3, whilst the lowest values of soil C, soil C/N ratio, and C storage were noted under the UNT. Significantly, the highest amount of soil N, N storage, and available nutrients (P, K, Ca, and Mg) were found under the LI9, followed by LI6, whereas these values were still less than those of the values under the UND over a 2-year period after soil compaction. Two years after applying litter mulch and incorporation, the recovery values of soil N, N storage, and available nutrients (P, K, Ca, and Mg) were significantly lower in LM3, as compared to the UNT treatment ( Table 2).
Results showed that soil physico-chemical properties were partially recovered due to the application of both litter mulch on the surface of the soil and incorporation into the soil, as compared to the UNT, although soil physical properties were still lower than the UND. Regardless of the litter mulch and incorporation application rate, the recovery of soil physicochemical properties was higher in the incorporation treatments than the mulch treatments. For both litter mulch and incorporation, the recovery level of soil physico-chemical properties increased as the application of litter amount increased (Tables 1 and 2).

Soil Biological and Microbial Properties
Soil biological and microbial properties significantly differed (p < 0.001) among different litter mulching and incorporation, UND, and UNT treatments ( Table 3). The greatest values of earthworm density and dry mass, fine root biomass, and soil microbial respiration were found under the UND, followed by LI9 and LI6, whereas the significantly least amount of these values was measured under the UNT treatment (Table 3).
Due to the litter mulching on the soil surface and litter incorporation into the soil, the recovery values of earthworm density and dry mass, and fine root biomass as well as soil microbial respiration, were significantly higher in the incorporation treatments (i.e., LI9, LI6, and LI3) than the litter mulching treatments (i.e., LM9, LM6, and LM3). However, the values of earthworm density and dry mass, fine root biomass, and soil microbial respiration did not return to the pre-harvest levels as observed under UND over a 2-year period, but these values were higher than the values of the UNT (Table 3).  Soil bulk density significantly and positively correlated with soil penetration resistance, sand, silt, and pH, and negatively correlated with total porosity, macroporosity, soil moisture, aggregate stability, clay, soil C, soil N, C and N storage, available nutrients (P, K, Ca, and Mg), earthworm density and biomass, fine root biomass, and soil microbial respiration (Table 4).

Soil Physical Properties
Machinery induced soil compaction without any mitigation-restoration technique could lead to changes in the nutrient cycle by changing plant diversity, plant uptake, moisture, temperature fluctuations on the soil surface, decreased soil microbial activity, and water infiltration rate [39,55]. Soil compaction on the skid trails resulted in an alteration in the natural processes of soil biochemical parameters, which led to the loss of soluble nutrients in the downstream water. Nitrate, as an inorganic form of N, is a major nutrient in forest soils that naturally contributes to the nitrogen cycle in the environment, which is easily removed by soil leaching [56]. In line with the current study, several studies showed that disturbance and compaction of soil, and rutting following heavy machinery traffic during logging operations, drastically led to the destruction of the litter layer, which causes the soil to become a harsh environment for the forest, resulting in a reduction in the nutrient cycle [14,16].
The results of the current study revealed that the higher level recovery of soil bulk density, total porosity, and macroporosity were detected in the litter incorporation treatments (i.e., LI3, LI6, and LI9), which was more than the litter mulching treatments (i.e., LM3, LM6, and LM9), as compared to the UNT. The high recovery of soil bulk density in the litter incorporation treatments can be attributed to the faster decomposition rate of litter that was easily mixed with disturbed soil and subsequent release in soil N. For both litter mulching and incorporation, the recovery level of soil bulk density increased as the application amount of the litter increased from 3 to 9 Mg ha −1 , as compared to the UNT.
Results of the current study revealed that the recovery level of soil penetration resistance was higher in the litter incorporation treatments than those of the litter mulching treatments, as compared to the UNT. Previous studies concluded that litter mulch application significantly accelerated the recovery of soil penetration resistance on the skid trails, which could be attributed to an increase in soil biological activities on the organic substrates and soil-root interactions [3,17,23], resulting in a restoration of soil structure and improvement of the soil's physical properties [7,14]. By application of litter mulch at a rate of 18 Mg ha −1 , Jourgholami et al. [23] found that the highest recovery of soil bulk density, total porosity, macroporosity, soil penetration resistance, and aggregate stability was observed in the beech-hornbeam-maple litter treatment, as compared to the untreated trails. However, these values showed significant differences with the UND treatment, even over a 5-year period.
The recovery level of soil aggregate stability was higher in the litter incorporation than the observed level in the litter mulching due to the high fine root biomass in the INC treatments. According to the results of the current study, litter incorporation treatment with a different application amount appears to restore soil physical properties more than the litter mulching treatment. Litter incorporation into soil increases the litter decomposition rate, which augments the soil moisture content and water retention, and improves the nutrients viability, leading to the restoration of the soil's physical properties [41,42].
In addition, Jourgholami et al. [57] pointed out that the full recovery of the soil physical properties (i.e., bulk density, total porosity, and penetration resistance) did not occur with hardwood sawdust mulch (beech and hornbeam) over a 6-year period, when compared to the control, due to the high amount of 36.5 Mg ha −1 , which resulted in an increase in soil acidity from high lignin compounds. Furthermore, by applying four litter mulch amounts (i.e., levels of 4.2, 8.1, 13.1, and 16.9 Mg ha −1 ), Jourgholami et al. [3] found that the highest recovery rate for soil bulk density, total porosity, macroporosity, penetration resistance, and aggregate stability were observed in the litter level of 16.9 Mg ha −1 followed by a litter amount of 13.1 Mg ha −1 .
Results of the current study revealed that clay particles were higher in the LI9 and LM9 than in those of the treatment with a low application amount, as was also reported by Jourgholami et al. [3]. Similar results were found by Jourgholami et al. [17], since the increase in the rate of litter mulching and incorporation from 3 to 9 Mg ha −1 resulted in the formation of a protective layer; this layer intercepted the raindrop erosive power and reduced clay particle detachment.
Incorporation of mulch with soil before the movement of harvesting vehicles on the surface of the skid trails led to the formation of a buffer layer between soil and machine wheels, which was able to absorb and reduce the pressure of the machine on the soil surface.

Soil Chemical, Biological and Microbial Properties
Results showed that soil pH values were higher in the litter mulching treatments than in the litter incorporation treatments. This aspect can be attributed to the high organic matter content on the surface of the soil due to the high-quality litter. High-quality litter increases the decomposition rate and enhances the cation exchange capacity, which favors soil biological activities as reported by Diao et al. [35].
Results of the current study showed that the highest recovery level for soil organic C and C/N ratio was observed in the litter incorporation treatment followed by the litter mulching treatment. According to previous studies, the highest level of soil organic C and soil C/N ratio can be attributed to the slower rate of mineralization for litter mulching treatment and the increase in soil pH in the LM3, LM6, and LM9 treatments [24,25]. The decomposition processes of carbon and nitrogen compounds following litter mulching and incorporation, as well as subsequent mixing with mineral particles, led to improvement in soil C and N. Soil organic C mainly contributed to enhance the structure of soil and sustained the food substrates, which led to an increase in microbial and biological activities.
This study showed that the recovery values for N and N storage were significantly higher in the litter incorporation treatments (i.e., LI9, LI6, and LI3) than the litter mulching treatments (i.e., LM9, LM6, and LM3), as compared to the UNT. The high-quality litter of hornbeam-maple, which was incorporated into the soil, resulted in an increase in the decomposition rate, improved soil pH, and, ultimately, accelerated N release. Results showed that the application of the high-quality litter of tree species, including hornbeam and maple litters, as mulch on the soil surface or incorporation into the soil increased the N delivery into soil layers, which ultimately regulated nutrient cycling [17]. Moreover, litter mulch C/N ratio was an important factor that regulated the decomposition rate and subsequent nutrient release into the soil layers [21]. Instead, applying litter with recalcitrant compounds regulated the nitrogen, lignin, and the C/N and lignin/nitrogen ratios in litter [58], which not only increased soil C/N ratio, but also led to the accumulation of the organic matter content on the soil surface [36,58].
The results showed that litter mulching on the soil surface was a determining factor in the improvement of the C storage (3.7-4.36 Mg ha −1 ), more so than the litter incorporation treatments (2.89-3.33 Mg ha −1 ). Hence, application of litter mulch and incorporation in the compacted soil can significantly increase the carbon sequestration, which plays a key role in accelerating the restoration of disturbed soil after heavy machinery traffic.
The results of the present study elucidated that the highest recovery level of available nutrients (P, K, Mg, and Ca) was observed in the litter incorporation treatments, including LI9, LI6, and LI3, followed by LM9, LM6, and LM3, when compared to the UNT. However, these values, for both litter incorporation and mulching, were still lower than the value of the undisturbed treatment (UND). Falling leaves on the surface of the soil, especially highquality litter, appeared to have a crucial role in regulating the food resources, including organic matter content and soil C, and as an important key factor for microbial activities [59]. In accordance with the results of the current study, previous studies showed that alteration in litter and topsoil properties, highly related to native and non-native tree species, were mainly attributed to the quality of litter [26,35]. After application of litter mulch ranging from 4.2-16.9 Mg ha −1 , Jourgholami et al. [3] determined that soil N, N sequestration, and available nutrients (i.e., P, K, Ca, and Mg) were at the highest level in litter mulch applied in the amount of 16.9 Mg ha −1 , and at the lowest level in litter treatment in the amount of 4.2 Mg ha −1 .
The highest recovery level of available soil nutrients (i.e., P, K, Ca, and Mg) in litter incorporation into the soil can be attributed to the labile components of hornbeam-maple litter, the low content of lignin, and the high content of N leading to the greater lignin/nitrogen ratio, which decomposed faster than the litter mulching on the soil surface. Similar observations were made by Maggard et al. [60], Aponte et al. [22], Langenbruch et al. [24], and Diao et al. [35].
In the current study, the recovery values of the earthworm density and dry mass were higher in the litter incorporation treatments (i.e., LI9, LI6, and LI3) than those of the litter mulch treatments (i.e., LM9, LM6, and LM3) which can be attributed to the lower C and C/N ratio, higher N, and available nutrients [3]. However, Bottinelli et al. [10] reported that the recovery of earthworm density and biomass did not occur over a 4-year period due to the harsh environment in the compaction-induced soil after logging operations. Both litter mulching on the soil surface and incorporation into the soil led to the stimulation of soil micro-organisms (i.e., algae, mosses, fungi, and bacteria) and macro-organisms like earthworms [61]. Induced compaction is thought to degrade the soil quality and may accelerate the formation of a hostile environment for populations of soil organisms, with an undesirable impact on ecosystem engineers such as earthworms [10]. Instead, adding litter mulch enhances the food substrates, which encourage earthworm movement into the soil, leading to the enhancement of soil structure and the nutrient cycle [61].
Apart from the effects of organic mulches on improving soil quality conditions, another important aspect is their role in soil and water conservation. Previous studies have revealed that skid trails were considered an important non-point source to produce runoff and soil erosion [3,26,29,30,62]. Mulch application can protect the surface of bare mineral soil and play a crucial role in soil and water conservation [38]. Previous studies indicated that application of mulch increased the soil cover, absorbed the raindrops hitting the soil surface, and enhanced the ground roughness of the soil surface, which resulted in an increase of water infiltration into the soil layer, reduced surface runoff, decreased soil detachment, and reduced sediment and nutrient losses [17,23,26,57]. Mulch application can led to altered hydrological connectivity at the headwater catchment [62]. In the current study, splash erosion was the main erosion process, rather than other erosion processes such as rill and interrill, due to the shorter length of plot as reported by Bagarello and Ferro [63]. In the Hyrcanian forests, Jourgholami et al. [26] concluded that the application of leaf litter with rate of 4.2-16.9 Mg ha −1 resulted in a decrease in runoff from 49% to 79% and sediment yield from 76% to 93%. By reviewing the studies in agricultural areas, Li [38] concluded that a mulch application rate of 6-8 Mg ha −1 can be effective at mitigating runoff and soil loss. However, López-Vicente et al. [62] demonstrated that check-dam installation in the skid trails areas can play an important role in suppressing sediment delivery into water bodies.
Results of this study demonstrated that the recovery levels of fine root biomass were higher in the litter incorporation treatments than in those of the litter mulching treatments, as compared to the UNT. Similar to our findings, Jourgholami et al. [39] reported that increasing the application amount of litter resulted in an increase of the recovery values of fine root biomass. Likewise, Jourgholami et al. [23] concluded that litter mulching on the soil surface led to a restoration of the values of fine root biomass. The results showed that induced compaction significantly reduced the amount of organic matter and depth of litter layer, which negatively affected the soil microbial biomass. The findings of this study revealed that the recovery values of soil microbial respiration were higher in the litter incorporation into soil treatment than in the litter mulching on the soil surface. The faster recovery values observed in the litter incorporation treatments can be related to the high-quality litter, with lower levels of carbon and lignin, which resulted in a faster decomposition rate and higher nutrient release.
The period of 1-2 years after machinery traffic on the skid trail is a very crucial period in which to apply best management practices (BMPs), such as mulching, to mitigate surface runoff and soil erosion. However, organic mulches, such as leaf litter, were decomposed one or two years after the fall. Hence, the effect of organic mulch on soil quality can appear one to two years after mulching, although the long-term effects of soil compaction may take a few years to several decades to dissipate.

Conclusions
In the present study, the effects of litter (a mixture of hornbeam and maple) mulching on the surface soil and incorporation into the soil, with three application amounts (3, 6, and 9 Mg ha −1 ), on the skid trails to restore the soil physico-chemical, biological, and microbial properties were tested, and compared to the untreated skid trail (UNT) and the undisturbed (control) treatments (UND) over a 2-year period after treating and skidding operations. Our findings demonstrated that the recovery values of physico-chemical, biological, and microbial soil properties were significantly higher in the litter incorporation treatments (i.e., LI9, LI6, and LI3) than those of the litter mulching treatments on the soil surface (i.e., LM9, LM6, and LM3). Results of the current study elucidated that physico-chemical, biological, and microbial soil properties were partially recovered over a 2-year period from treatments, as compared to the untreated skid trail (UNT), but the full recovery of soil properties required more time before it returned to the pre-harvest value or values of the undisturbed treatments (UND). Regardless of litter mulching or incorporation, as litter amount increased from 3 to 9 Mg ha −1 , the recovery level for physico-chemical, biological, and microbial soil properties increased compared to the UNT, but these values were still different from the values of the undisturbed treatments (UND) even 2-years after soil compaction and litter applications. According to the current results, litter incorporation into the soil was more effective for soil restoration than litter mulching. Some BMPs that can be applied to mitigate the adverse impacts and regulate soil and water conservation are: • Application of the harvesting residues and organic mulch before the forestry vehicle traffic or incorporation instead of mulching after harvesting operations; • Litter incorporation into the soil with the application of an amount ranging between 6-9 Mg ha −1 as an effective method to restore the compaction-induced soil to the same conditions before the logging operations on the skid trails; • Mulching and incorporation of organic materials on the skid trails to mitigate surface runoff and rill and interrill erosion.