Effect of Heavy Machine Trafﬁc on Soil CO 2 Concentration and Efﬂux in a Pinus koraiensis Thinning Stand

: Mechanized timber harvesting is cost efﬁcient and highly productive. However, mechanized harvesting operations are often associated with several environmental problems, including soil compaction and disturbance. Soil compaction impedes air circulation between the soil and atmosphere, which in turn results in increased concentrations of CO 2 within soil pores. In this study, we investigated the impact of forest machine trafﬁc on soil conditions to determine soil CO 2 efﬂux ( F c ), and soil CO 2 concentrations ( S c ). Field measurements included soil bulk density (BD), soil temperature (ST), soil water content (SWC), F c , and S c over a 3-year period at a speciﬁc thinning operation site (Hwacheon-gun) in the Gangwon Province of Korea. To assess the soil impacts associated with machine trafﬁc, we established four machine-treatment plots (MT) characterized by different geographical and trafﬁc conditions. The results revealed that BD, S c , and SWC increased signiﬁcantly on the disturbed track areas ( p < 0.05). Furthermore, reduced soil F c values were measured on the soil-compacted (machine disturbed) tracks in comparison with undisturbed (control) areas. Vari-ations in BD, SWC, and S c were signiﬁcantly different among the four MT plots. Additionally, in comparison with undisturbed areas, lower F c and higher S c values were obtained in compacted areas with high soil temperatures.


Introduction
Mechanized equipment is necessary for more efficient harvesting operations [1]. However, heavy machine traffic during harvesting operations can cause disturbance to soil ecosystems, which results in long-term impacts on soil production [2]. Soil compaction is one of the major negative impacts associated with forest machinery traffic [3,4]. Heavy machine traffic may cause more severe soil compaction in forests than in croplands, since it involves machine traffic as well as the cutting, pulling, pushing, lifting, and transport of timber during the harvesting operation [2]. Froehlich et al. [5] reported that soil compaction can persist for around 15 years in temperate forests that have no frosting of soil water in the winter season. DeArmond et al. [6] found persistent compaction after a period of 30 years following logging operations. Additionally, it has been reported that soil biological disturbance can be recovered within 20 years in sandy neutral soil following compaction [7]. Severe soil disturbances and impacts on soil physical properties, such as soil compaction, rutting, and soil displacement can occur during mechanized harvesting operations [8][9][10].

Site Description
To study the effects of soil compaction caused by the traffic of forestry machines on various soil properties, a 40-year-old Pinus koraiensis Siebold and Zucc. (Korean pine) stand located in the experimental forest of National Institute of Forest Science (38 • 00 58.7 N, 127 • 48 31.2 E; 370 m elevation) in the Gangwon Province, Korea was chosen as the study site. Terrain slope in the study area ranged between 21% and 57%, while the soil type corresponded to the Mui Series (coarse loamy, mixed, Typic Humudepts) based on U.S. Soil Taxonomy. Soils in the area consisted of dark brown sandy loam or loamy sand A horizons. The organic layer above the surface soil had a mean thickness of 4 cm and a mean density of 1.8 kg/m 2 . The dominant forest species in the study area was Pinus koraiensis with a stocking density of 448.6 trees per ha and a timber volume of 279.9 m 3 per ha. The 30-year mean annual precipitation was 1358.7 mm, with average minimum and maximum temperatures of -18.5 • C and 35.2 • C, respectively.

Experimental Design
The route and traffic intensity (number of machine passes) were collected using the GPS of the forestry machine at the thinning operation site. Each pass comprised a round trip of the machine. The thinning operation applied the whole tree method and was conducted in May 2018. Two forestry machines were used during the thinning operation: a harvester and a skidder. The harvester had a total weight of 6030 kg, with the machine (DX55 A-MT) and the grapple saw (GMT035) having a weight of 5800 kg and 230 kg, respectively. The contact area of each harvester wheel was 15,580 cm 2 (205(L) × 38(W) × 2EA), with a contact pressure of 0.38 kgf/cm 2 . The skidder weighed 5760 kg. Each wheel of the machine had a contact pressure of 0.36 kgf/cm 2 , and the contact area was the same as the harvester.
The soil-compacted areas (disturbed tracks) were visually identified by the tire marks over an average length of 25 m in each MT plot. An undisturbed soil area (control) of the same length and within 5 m of the disturbed track was selected in each plot to reduce the effect of environmental conditions. Based on the results from previous studies that have reported the impact of traffic intensity and terrain slope on compaction, four machinetreatment (MT) plots were selected: MT1-with one round-trip performed by the harvester on a 54% terrain slope, MT2-with one round-trip performed by the harvester on a 34% terrain slope, MT3-with one round-trip performed by the harvester and skidder on a 22% terrain slope of 22%, and MT4 with three round-trips performed by the harvester and two round-trips performed by the skidder on a 33% terrain slope (Table 1, Figure 1).  However, it was difficult to select a MT plot with a specific traffic intensity and terrain slope in the thinning stand since the out-row thinning method is limited in most of Korean forest areas that have a rugged terrain. Therefore, four MT plots were selected to determine the effect of low and high machine traffic intensity on soil compaction.
Finally, five measurement points along the disturbed track and control areas were selected in each MT plot; at each point, measurements of ST, SWC, F c , and S c were performed.

Soil Physicochemical Properties
Soil samples used for the analysis of initial soil pH, OMC, and texture in the study site were collected during May 2018 (prior to the mechanized timber harvesting operation) ( Table 2). The collected soil samples (>1 kg) were transported to the laboratory, and subsequently analyzed at the Korea Forestry Promotion Institute in Seoul, Korea. Soil pH ranged from 4.7 to 5.1 at a soil depth of 0-30 cm ( Table 2). The original OMC in each study plot is provided in Table 1. The highest mean OMC was found in MT2 followed by that in MT3 (Table 2). Soil was classified as mostly sandy loam and loamy sand through an analysis of the soil physical properties (Table 1). Similarly, OMC was measured to determine the effect of compaction caused by the machine traffic. Soil samples for the analysis of OMC were collected during October 2020. In addition, soil BD was measured to assess the effects of compaction. Soil samples collected for the analysis of BD were obtained from five measuring points at locations other than those used to measure ST, SWC, F c , and S c . These soil samples were collected from the disturbed and undisturbed (control) tracks on three occasions at soil depths of 0-10, 10-20, and 20-30 cm in September 2018, 2019, and 2020. The soil samples were dried for 48 h in a laboratory oven at 105 • C and in containers of 100 cm 3 . BD (representing soil dry bulk density) was calculated using Equation (1): where BD is dry bulk density, BD w is wet bulk density, and M w is the mass of water (g cm −3 ).

Soil Temperature, Water Content, and CO 2 Efflux and Concentration
In each MT plot, measurements of ST, SWC, F c , and S c were taken at five selected measurement points from the disturbed tracks and undisturbed (control) areas. These variables were measured once a month from September 2018 to September 2020 (a total of 18 measurement days). This excluded measurements in the winter season (December, Jan-uary, and February), since in this period the soil is generally frozen; also, no measurements were performed in March 2020, owing to COVID-19-related travel restrictions. ST and SWC were measured to determine the response of F c and S c resulting from soil compaction. ST was measured using a digital thermometer TP3001, and SWC was measured with probe sensors (TDR 300, FieldScout) in the form of volumetric water content (%). At each of the measurement points, ST and SWC were measured five times.
F c and S c were measured to assess gas diffusivity and soil CO 2 productivity resulting from soil compaction caused by machine traffic. F c was measured at the same time as ST and SWC (total of 18 days) at each measurement point. All measurements were carried out during daytime (11:00 to 16:00). The PVC chambers, 13 cm in diameter and 16 cm in length, were embedded 1.4 cm into the soil surface. After closing the chamber, a carbon dioxide probe (GMP343, Vaisala CARBOCAP ® , Helsinki, Finland) was inserted, and F c was measured for 5 minutes. F c was calculated from the change rate in CO 2 concentration (ppm) using the closed dynamic chamber method [39] as per Equation (2): where V is the total chamber volume of the system, A is the area covered by the chamber, m w is the molecular weight, and m v is the CO 2 volume. S c was measured on the same day as the other variables. In July 2018, a total of 40 vapor tips (mesh covered for air suction) were installed at a soil depth of 10-15 cm, and silicon tubes (3-mm internal diameter) connected to the vapor tip were pulled out of the soil surface. To avoid potential disturbance of these pre-treatment procedures on the subsequent measurement of S c , data collection commenced only after two months of completing the pre-treatment in September 2018. The pump (GM70, Vaisala CARBOCAP ® , Helsinki, Finland) was connected to the other side of the tube to suck air from the soil, and S c in soil pores was measured using a carbon dioxide probe (GMP222, Vaisala CARBOCAP ® , Helsinki, Finland). The vapor tips, made of metal and meshed for sucking air in the soil, were installed at a soil depth of 10 cm and connected to a silicon tube. The other side of the tube was connected to a pump (GM70, Vaisala CARBOCAP ® , Helsinki, Finland) for sucking soil air, which was allowed to pass through the carbon dioxide probe (GMP222, Vaisala CARBOCAP ® , Helsinki, Finland). At each of the measurement points, S c was measured for 5 min after F c measurements had been completed to avoid disturbances while sucking the soil air.

Data Analysis
All measurements were conducted at five points on each machine disturbed track and control area. Monthly average data for F c , S c , ST, and SWC, as well as yearly average data for BD were calculated as the mean of the five replicated data points on the machine disturbed track and control area to avoid pseudo-replication.
To assess the effect of soil compaction caused by machine traffic, we calculated the relative change between the soil compaction values on the disturbed track and control area using Equation (3): The Anderson-Darling normality test and Bartlett test were used to evaluate the normality and homogeneity of variance, respectively. To assess the effects of machine traffic-related factors (i.e., traffic intensity and terrain slope) and/or year on BD, F c , S c , ST, SWC, analyses of variance (ANOVA) with repeated measures were performed using a generalized linear model (GLM). A Post-Hoc test was conducted with Tukey's pair-wise comparisons. All analyses of variance and normality were conducted in SPSS v. 25, and significance was set at p ≤ 0.05.
To assess the effects of ST on F c and S c , a first-order exponential function Q 10 was fitted to Equation (4): where y represents the measured F c or S c , β 0 and β 1 are the fitted parameters. Q 10 values were calculated using Equation (5): F c and S c were measured on the same day during daylight hours (11:00 to 16:00), while ST and SWC were measured at the same time as F c and S c . Therefore, a normalization of soil CO 2 values to daily mean ST was not required. However, F c and S c were normalized to a soil temperature of 10 • C during the measurement period to assess the effect of SWC on soil CO 2 and to avoid changes in soil CO 2 associated with the seasonal variation in soil temperature.
In addition, a polynomial expression was fitted as per Equation (6) using F 10 and S 10 values to assess the effect of SWC on F c and S c .
Multiple nonlinear regression was performed using a SigmaPlot demo (Systat Software, Inc., San Jose, CA, USA) for Windows (Microsoft, Redmond, WA, USA).
The covariance structures of regression Equation (6) was chosen by minimizing the corrected Akaike information criteria (AICC) [26].

Soil Physical Properties following Machine Traffic
Average BD following compaction was significantly higher (average increase of 14%) at p = 0.001 on the disturbed track with respect to the control area (Tables 3 and 4). This change appears to be the effect of compaction due to machine traffic. Soil porosity decreases in soils with increased BD, which results in higher SWC values [12,40]. In comparison with the control plot, we detected an increase in BD on disturbed tracks, although these differences were not statistically significant in any MT plot at different soil depths (Table 3). In addition, although dissimilar changes in BD following compaction were found among the four MT plots, these differences were not statistically significant at p = 0.133 (Table 3). The largest increase in BD (23.0%) on disturbed tracks was found in MT4 (23.0%), followed by MT1 (12.3%). Table 3. Statistical significance probability values for the effect of variables on soil bulk density (BD), soil CO 2 efflux (F c ), soil CO 2 concentration (S c ), soil temperature (ST), and soil water content (SWC) from ANOVA using general linear models (GLM).

Variable
Pr   In comparison to the control area, BD on disturbed tracks increased more in areas closer to the surface soil across all MT plots, although these differences were not statistically significant at p ≤ 0.05 (Table 4). In a previous study, BD was reported to increase by 28% (0-10 cm soil depth), 12% (10-20 cm soil depth), and 34% (20-30 cm soil depth) as a consequence of forwarder traffic (1 and 5 round-trips) on silt loam in a mixed forest (consisting of 60% Pinus koraiensis) in Korea [41]. Goutal et al. [1] observed an increase in BD of 6% (30-40 cm soil depth) and 27% (0-10 cm soil depth) following the traffic of a forwarder, while Jankovský et al. [10] reported an increase in BD ranging from 35 to 38% following machine traffic. These studies also confirmed that increased BD on the soil surface can lead to poor gas diffusion between soil and atmosphere [42].
Excepting in MT2, OMC was higher on disturbed tracks than in control areas, although the differences were not statistically significant at p ≤ 0.05 (Table 4). In addition, there was no statistically significant difference (p ≤ 0.05) in OMC between MT1 and MT4 following compaction (Table 4).
Following compaction, SWC increased significantly (average increase of 29.8%) on the disturbed tracks with respect to the control areas in all MT plots (p ≤ 0.001) ( Tables 3 and 5). In addition, SWC increased in all seasons but there was no significant difference among the seasonal variation of SWC on the disturbed tracks (Table 5).
An increase in SWC led to a slight decrease in ST on the disturbed tracks in comparison to that in the control areas, but this change was not statistically significant (p = 0.053). In addition, there was no seasonal difference of changes in ST following compaction (Table 5). Table 5. Seasonal average ± standard error of soil temperature (ST), soil water content (SWC) obtained on disturbed tracks (T) and control areas (C) for the four machine-treatment (MT) plots following machine traffic during the thinning operation in May 2018. The data were collected every month from September 2018 to September 2020 except during the winter months (December, January, February).

Soil CO 2 Efflux and CO 2 Concentration following Machine Traffic
The increase in BD and SWC impedes gas diffusivity following machine traffic, resulting in a statistically significant (p = 0.025) average reduction of 23.7% in F c (Table 3). In addition, the correlation was less obvious and more complex on the disturbed tracks ( Figure 2). Additionally, S c on disturbed tracks increased significantly at p = 0.002 (average increase of 10.6%) with respect to the control areas; this was a consequence of the impact on gas diffusivity following mechanized timber harvesting (Figure 2). The effect in F c was statistically significant (p ≤ 0.001) for the interaction of soil compaction and measurement period (Table 3). This may be because F c had a substantial decline on disturbed tracks from June to September in comparison to other periods ( Figure 3a). However, there was no significant difference (p = 0.514) among the four MT plots for the interaction effect of soil compaction and measurement period (Table 3, Figure 3a).
S c was high from July to September both on disturbed tracks and control areas (Figure 3b). In addition, the increase in S c on disturbed tracks was higher from July to September in comparison to the other periods; this effect was associated with the largest decrease in F c that occurred on disturbed tracks in the same periods (Figure 3a). However, the interaction effect between soil compaction and measurement period on S c was not significant (Table 3, Figure 3b).

Effects of Machine Traffic Conditions on Soil Compaction
The effects of compaction caused by machine traffic on S c and SWC depended significantly on machine traffic conditions (traffic intensity and terrain slope) at p = 0.024 and p = 0.012, respectively (Table 3, Figure 4). The change rate in SWC and S c differed significantly among the MT plots ( Figure 4). The highest increase in BD was observed in MT4 (23.0%), followed by MT1 (12.3%). SWC also showed the largest increase in MT4 (46.4%), followed by a 22.9% increase in MT1, 26.6% in MT2, and 33.8% in MT3. Although there were no significant differences in F c among the MT plots, the largest change rate in F c occurred on the disturbed track of MT1. The largest increase in S c was observed in MT1 ( Figure 4). While MT4, had the largest increase in BD and SWC, it also had the lowest F c reduction rate and the second largest increase in S c (Figure 4). These results can be explained to a greater extent by the change rate in BD and S c .

Yearly Difference Associated with the Impact of Machine Traffic on Soil Characteristics
Significant effects on BD's change rate were observed for the interaction between compaction (C), machine traffic condition (MT), and measurement period (P) (p = 0.004) (Table 3, Figure 5a). Dissimilar results were obtained for BD's change rate on disturbed tracks in comparison to control areas; thus, the change rate increased in MT1 and MT2 throughout the year and decreased in MT3 and MT4 during the same period (Figure 5a). MT4 had the largest average BD increase on disturbed tracks with respect to control areas; in 2018 this increase reached 40.4% after compaction (Figure 5a). Since 2018, the growth rate in BD declined significantly in MT4 (p ≤ 0.05) (Figure 5a). In MT3, the growth rate in BD decreased over the years (Figure 5a).  The interaction effect for the variables compaction (C), machine traffic condition (MT) and measurement period (P) was not determined except for BD; however, the change rate in S c was significantly different during the three years that follow the thinning in each MT plot (Table 3, Figure 5b). In MT4, the growth rate in S c increased while BD decreased throughout the year (Figure 5b). Thus, the effect of compaction on S c was not simply caused by the effect of pore compression, but also by the interruption of CO 2 gas exchange between soil and atmosphere; this latter effect was a consequence of the diminished air circulation that resulted from the interruption in pore continuity [12,13]. On the other hand, BD increased by 40.4% and S c was lower on the disturbed track with respect to the control area in MT4 during 2018.

Effects of Compaction on ST and Soil CO 2
Unfortunately, soil CO 2 measurements in our study could not be performed in winter since the soil was frozen during this season. However, the Q 10 values showed that F c and S c decreased and increased following compaction in the summer season, respectively (Table 6). Furthermore, the exponential function between ST and F c or S c had a distinctive positive relationship, revealing that soil disturbance occurred mainly when the soil temperature was high after compaction.

Changes in Gas Diffusion following Soil Compaction
In general, as CO 2 production increases, so does CO 2 efflux from the soil surface [18]. Therefore, F c and S c probably have an obvious positive correlation in machine-disturbed soils. In this study, F c correlated well with S c in the control area. As shown in Figure 2, gas diffusion was inhibited in all MT plots following compaction. The results are similar to those reported by Goutal et al. [1] and reveal impediments of gas diffusion following compaction.
Increased BD results in reduced macropores, and the disconnection of pores could impede gas diffusion in compacted soils [18,34] which might lead to an increased concentration of CO 2 in soil pores. F c values were lower on disturbed tracks than in control areas following compaction [3,43,44]. Greater values of F c in control areas may be the result of CO 2 production in the soil; this effect is more predominant active under aerobic conditions than anaerobic conditions [45]. The changes in soil CO 2 obtained in our study are in agreement with those from previous studies, which have reported higher S c values in disturbed than in undisturbed areas following compaction [1,17,34,46]. Therefore, significant increases in S c in this study are a consequence of the impact on gas diffusivity following mechanized timber harvesting (Figure 2).

Impact of Machine Traffic on Soil Compaction
The largest increase in BD was observed on the disturbed track of MT4, followed by that of MT1 (Figure 4), which reveals that machine traffic on steep terrain and high machine traffic intensity can result in severe soil compaction. Agherkakli et al. [28] studied the change in rut depth and BD after 1, 5, and 9 passes of forestry machines on terrain slopes of <20% and >20%. They reported significantly higher BD and rut depth values when terrain slope exceeded 20%. They also showed that traffic intensity resulted in increased BD values. Najafi et al. [27] also reported that machine traffic on slopes greater than 20% led to significant soil disturbance.
During 2018, F c in MT4 was lower on the disturbed track than in the control area, while root respiration and microbial activities were severely hampered by soil compaction.
In a previous mechanized harvesting study, a reduction in root depth of more than 20% was reported on disturbed tracks with respect to undisturbed control areas [47]. As a result of S c , the effects of compaction were not retained in MT2 and MT3. However, in MT1 and MT4, compaction was significantly higher on the disturbed track than in the control area, and S c continued to increase from 2018 to 2020. Therefore, even when the number of machine passes was just one, the effect of compaction was quite evident in areas of steep terrain; the same effect was observed when the number of trips was larger than five.
Our study did not include cases where the number of roundtrips of machine was three or more. However, the results can be used to estimate values in scenarios where the number of roundtrips is equal to five and in cases when the number of passes is just one and terrain slope is high.

The Soil Respiration Variance in Soil Compacted Area
Increased S c in soil pores can hinder future root growth and soil microbial activities [40]. Conlin and Driessche [46] concluded that since the change in S c following compaction is not significantly related to whether organic matter is removed or not, the origin of most respiratory activities might be associated with autotrophic rather than heterotrophic activity. However, Hanson et al. [48] reported that, on average, root respiration accounts for about 45% of F c in forested ecosystems, with values that range between 20% and 80%. In addition, Striegl and Wickland [3] reported a value of 35% for a mature Saskatchewan jack pine forest. In this study, no measure was taken to distinguish between root respiration and microbial activity, which made it difficult to determine if the origin of soil respiration accounted for the majority. However, it can be assumed that the effect of compaction on roots was greater than that on microorganisms because the OMC was not significantly difference between disturbed tracks and control area after compaction.

The Seasonal Soil CO 2 Impacted by Soil Compaction
Buyanovsky and Wagner [49] reported that the exponential function explained the increase in S c as soil temperature went up from 10 to 20 • C. Epron et al. [18] concluded that a reduction in F c represents a degradation in soil CO 2 production, while an increase in S c means a slow CO 2 diffusivity. Although root growth or microbial activities were not measured in our study, a Q 10 reduction in F c measured on disturbed tracks following compaction may have decreased especially in the summer. The restricted gas exchange between soil and atmosphere following compaction increase S c and decrease the oxygen concentrations in the soil [46]. In addition, Allman et al. [40] reported that a critical S c is believed to fluctuate around 0.6%; other studies have indicated that an increase in S c over this limit might impede the root growth of seedlings [19,20]. The main growing season in our study site was summer when both air temperature and precipitations are high. Regardless of when compaction occurs, our results suggest that soil compaction can lead to reduced soil productivity owing to reduced root growth and microbial activities.

Soil Recovery after Soil Compaction
The process of root penetration and annual freeze-thaw cycles may allow the recovery of macroporosity on the topsoil with the corresponding decrease in BD [50,51]. However, Goutal et al. [1] and Epron et al. [18] reported that there was no evidence of recovery either immediately or within five years after compaction. Jakobsen [52] observed that the effect of soil compaction persisted for 30 years following machine harvesting in a temperate forest without freeze-thaw cycles. In addition, Labelle and Jaeger [50] observed an increase in BD following compaction, which persisted for 5 years in soils exposed to swelling and shrinking; this was due to the high variation in the precipitation rates throughout the year and the occurrence of the freeze-thaw cycles. In our study, BD recovery following compaction was not observed (Figure 5a). In addition, growth rate of S c in MT1 and MT4 increased year by year after compaction. This suggests that machine traffic on steep slopes or high-intensity machine traffic during mechanized harvesting should be avoided.

Conclusions
This study aimed to quantify the compaction of forest soil caused by mechanized timber harvesting. After machine traffic, soil bulk density on disturbed tracks increased significantly with respect to control areas. This resulted in a reduction in soil porosity and pore connectivity. Consequently, soil water content also increased after soil physical degradation due to the reduction in infiltration. In addition, an increased soil bulk density and water content led to a reduction in gas diffusivity. This observation was supported by a decrease in soil CO 2 efflux and a significant increase in soil CO 2 concentration after machine traffic. In addition, temperature coefficient (Q 10 ) values were lower for soil CO 2 efflux and higher for soil CO 2 concentration across all-disturbed tracks in comparison to the control areas, regardless of the machine traffic conditions. Soil recovery from compaction was not observed within three years after the mechanized thinning operation. We conclude from these results, that machine traffic during mechanized timber harvesting can lead to a reduction of soil productivity associated with root growth and soil microbial activity.
Soil compaction caused by machine traffic was significantly affected by traffic intensity (number of passes) and terrain slope; this was confirmed by the change rate in soil bulk density, water content, and CO 2 concentration on disturbed tracks in comparison with control areas among the MT plots. In our study, disturbed tracks with moderate terrain slope and high machine traffic (5 round-trips) showed a significant increase in soil bulk density, water content, and CO 2 concentration; the same was observed for disturbed tracks on steep terrain and low machine traffic (1 roundtrip). On the contrary, disturbed tracks with a moderate slope and machine traffic (1 or 2 roundtrips) did not show any significant increase in soil CO 2 concentration. Based on our research, we can conclude that soil compaction is mainly caused by a high machine traffic intensity and operation of the machine on steep terrain during mechanized harvesting. To reduce soil compaction in small scale areas of steep terrain, it is highly recommended to plan and design the skid trails properly considering environmental and economic factors.
In our study, field data was collected from a forest stand that had been subjected to mechanized harvesting in order to assess the effects of soil compaction following machine traffic. Nevertheless, the study has two main limitations. The first is that the various conditions associated with machine traffic were not repeatedly investigated to determine the degree of soil compaction in relation to traffic intensity, terrain slope, and type of forestry machine. The second is that the 3-year measurement period was not long enough to observe soil recovery from the disturbance caused by the machine operations. Future research should consider long-term studies that include plots of similar conditions (i.e., traffic intensity and terrain slope) to the ones of the present study. These studies should also include measurements and analyses regarding the response of soil microbes and roots to compaction.

Conflicts of Interest:
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.