## Influence of slope on physical soil disturbance due to farm tractor forwarding in a Hyrcanian forest of northern Iran

iForest - Biogeosciences and Forestry, Volume 7, Issue 5, Pages 342-348 (2014)
doi: https://doi.org/10.3832/ifor1141-007

Research Articles

Tractor-based systems are the most common type of small-scale timber harvesting equipment in the Hyrcanian forest, where farm tractors equipped with a 2-wheel trailer typically perform forwarding operations for pulpwood and fuelwood. The objective of this study was to evaluate the effects of different frequencies of machine passes on soil bulk density (BD), penetration resistance (PR) and total porosity (TP) in different slope gradient and forwarding direction (i.e., uphill vs. downhill). The study, performed on a brown soil with a silt loam to loamy texture, was designed as a factorial experiment with three slope gradient/forwarding direction classes (0-10% downhill forwarding slope, 10-20% downhill forwarding slope and 0-10% uphill forwarding slope) and four traffic intensity classes (0 [undisturbed control area], 1-4, 5-10 and >10 machine passes). Responses were assessed at four soil profile depth classes (5, 10, 15 and 20 cm). BD and PR significantly increased and TP significantly decreased with increasing machine traffic, with magnitudes of change that were not consistently greater in the upper 5 cm compared to the 20 cm soil profile depth class. Whereas the majority of changes in BD and TP occurred after fewer than 5 passes, large increases in PR occurred even after 10 passes. Changes of BD, PR and TP in response to machine traffic differed significantly among slope gradient/forwarding direction classes. Increases in BD and PR and decreases in TP were largest on the 0-10% uphill forwarding slope, followed by the 10-20% downhill and 0-10% downhill forwarding slope. We conclude that the application of farm tractors in forwarding operations should be limited to gentle slopes and uphill forwarding with these machines should be avoided.

# Introduction

The entry of forestry machineries into forest has resulted in some cases in levels of soil disturbance that may have adverse effects on forest ecosystems and subsequently on long-term tree growth ([27], [25], [21], [29]). Numerous studies have documented changes in physical soil properties following harvesting machine traffic, such as increased soil compaction, soil displacement, rutting and soil profile disturbances ([1], [44], [32], [34], [12], [40], [13], [20], [42], [46], [4], [3], [31]). For instance, soil compaction, in which soil pores (especially macropores) are destroyed and surface aggregates are broken down ([4]), depend on the soil texture and moisture at the time of harvest ([26], [6]), but also on the number of machines passing on a specific point. Generally, soil compaction occurs in the first ten passes of a vehicle, but most of the compaction occurs in the first three trips; subsequent passes generally have little additional effect ([16], [12], [13], [4]). So it would be appropriate to concentrate machine traffic on carefully designed skid trails in such a way to minimize the affected area and leave most trees undisturbed ([45], [4]).

It is unclear, however, to what extent the slope gradient (i.e., longitudinal and transversal slope) of the skid trail and the direction of machine traffic (i.e., uphill vs. downhill) could modify the impacts of vehicle traffic on the soil. Whereas some studies have failed to document an effect of slope on bulk density following harvesting traffic (e.g., [38], [22]), others have shown that soil disturbances increase with slope gradient ([14]), possibly because machines slip more and remain in a given place for a longer time on steep terrain, causing soil pulling and dragging ([16]). Thus, steepness may have an even stronger effect on soil disturbance than that of the soil moisture ([24]), but possibly less so than that of the traffic intensity ([14]).

Soil disturbance is also strongly affected by the type of harvesting equipment used. Small-scale harvesting technology, such as farm tractors or better small skidder or other specific machines (such as ATV), is thought to result in lower levels of residual stand and soil damage because the equipment is typically smaller and lighter than regular skidders ([41], [36], [39]). A recent study dealing with soil compaction has found that the highest relative increase in bulk density was associated with the use of bulldozers, the lowest with the use of a forwarder, and intermediate levels with the use of a farm tractor and a rubber-tired cable skidder ([9]). However, farm tractors have been designed primarily for use in an agricultural rather than a forest environment and must be modified for use as wheeled skidders, such as through the addition of a logging winch, a roll-over canopy, a blade, belly pan protection, valve stem protection, wheel chains and wheel weights ([10], [33], [37], [11], [19], [35]).

In the Hyrcanian forest of northern Iran, where rough topography and use of the sheltercut system favor small-scale harvesting operations, farm tractors often pull an additional two-wheeled trailer to forward the wood raised completely off the ground. The felling and on-site processing of felled trees with large diameters are typically done by chainsaw ([23]). Mechanization is largely limited to the transport function, whereby logs of 5-15 meter length are extracted with rubber-tired skidders to roadside landings and farm tractors and mules are the most common timber extraction systems for pulpwood and fuelwood. The advantages of using relatively light farm tractors with large-dimension tires in the rear axle is the large contact area of the tire with the soil and the reduced ground pressure. Although the addition of a trailer increases the payload per trip, thereby reducing the number of equipment passes with its associated soil disturbances, the trailer also reduces tractor maneuverability, potentially creating deep ruts in steep terrain and in uphill forwarding.

We hypothesized that the extent of the soil degradation from farm tractor forwarding operations along the upper 20 cm of the soil profile would increase with traffic frequency, but that these increases would vary with slope gradient and direction of forwarding (uphill or downhill). To test this hypothesis, we quantified soil bulk density, penetration resistance and total porosity in skid trails on different slope gradients with varying forwarding direction and traffic intensities.

# Materials and methods

## Study sites

The research was carried out in Compartments 311 and 319 of the Gorazbon district in the Kheyrud Forest in the Hyrcanian forest region of northern Iran. The studied forests range in elevation from 1050-1270 m a.s.l. with a prevailing southern aspect. Average rainfall ranges from 1150 to 1260 mm year-1, with peaks in summer and autumn. Average daily temperatures range from a few degrees below 0 °C in December, January, and February to +25 °C during the summer. Soils are classified as well-drained forest brown soils (Alfisols) with a soil texture ranging from silt loam to loamy. The relief of the study area is characterized by carbonate rocks (limestone and dolomite) and their numerous karst forms: crests, carks in limestone, under-ground flows and sink hole structure. The study area is dominated by natural forests with native mixed deciduous tree species, including Fagus orientalis Lipsky, Carpinus betulus L., Alnus subcordata C.A.Mey, and Quercus castaneifolia C.A. Mey. The prevalent silvicultural treatment is a combination of the group selection and single tree selection methods, resulting in uneven-aged stands. Based on field data collected in October 2012, the average growing stock in compartments 311 and 319 was 289.1 m3 ha-1 and 555.2 m3 ha-1, respectively. A total of 240 trees (9 trees ha-1) with a total volume of 940.2 m3 (33.8 m3 ha-1) were marked in compartment no. 311 and a total of 210 trees (5 trees ha-1) with a total volume of 1177.3 m3 (27.5 m3 ha-1) were marked in compartment no. 319.

## Experimental design and data collection

A farm tractor equipped with a 2-wheel trailer (2-driving wheel) was used for forwarding pulp- and fuelwood. Technical specifications of the farm tractor were: UTB/ Universal 650 Engine UTB, diesel, 4-cylinder; Displacement: 4.8 l; Fuel capacity: 98.0 l; Mass: 3999 kg; Engine Power: 46.2 kW; Length: 414 cm; Width: 199 cm; Height: 260 cm; 4WD; and trailer Mass: 1700 kg. At the time of forwarding, weather conditions had been very dry and warm, and these conditions remained constant during the extraction operation. The soil moisture was 20-65% and the organic matter content was 11%. The slope of the skid trail ranged from 0 to 20% and extraction distances to the landing area near the forest roadside were about 220 m. A Suunto inclinometer was used for measuring slope in the skid trails. Calculation of slope was weighted as follows ([28] - eqn. 1):

$$S = \frac{\sum_{i=1}^{n}{d_{i}s_{i}}} { \sum_{i=1}^{n} {d_{i}} }$$

where S is the average slope along the skid trail (percent), d is the distance between two points in the sample (m), s is the slope between two points in the sample (percent), i is the sample number, and n is the number of samples.

The average forwarded volume in each trip was 2280 kg. Sampling transects were placed on different slope gradient classes along the designated skid trail (Fig. 1a). To ensure that traffic was restricted to the same tracks, painted sticks were implanted in the center of the skid trail as a guide for machinery operators. Before the forwarding operation with the farm tractor began, sampling plots were established in different skid trail segments on three slope gradients (i.e., 0-10% downhill, 10-20% downhill, and 0-10% uphill forwarding) that would be exposed to different levels of machine traffic, i.e., 0 (undisturbed), 1-4, 5-10 and >10 machine passes. A machine pass implies a drive both back and forth (i.e., with and without a wood load) on the selected trail segment. Each combination of slope by traffic frequency was replicated three times. Prior to and after the forwarding operations, soil response measurements (i.e., bulk density, resistance to penetration and total porosity) were obtained in each replicate at soil depths of 5, 10, 15 and 20 cm in the wheel rut as well as in control sample points adjacent to the skid trail (Fig. 1a). Organic horizons were removed from the soil surface prior to bulk density and penetration resistance measurements, so that measurements were referenced to the mineral soil surface. Soil sample cores were taken from the top mineral soil using a thin walled steel cylinder, 40 mm long and 56 mm in diameter, driven into the soil by a hammer-driven device. After extracting the steel cylinder, soil cores were trimmed flush with the cylinder end and extruded into a plastic bag for transport to the laboratory. Samples were weighed on the day they were collected and again after oven-drying at 105 °C for 24 h to determine water content and bulk density. Soil penetration resistance (PR) was measured using a hand-held soil penetrometer (Fig. 1b). The study yielded 144 PR readings (samples) before and after treatment. When the penetrometer had hit a buried rock or tree roots, the readings were ignored.

Fig. 1 - Experimental design of the field trail for the BD and PR sampling (a) and measurement of PR using a hand-held penetrometer (b).

Soil porosity (TP) was determined by the following equation (eqn. 2):

$$TP= \frac{(PD-BD)}{PD} \cdot 100$$

where PD is the particle density measured by a pycnometer on the same soil samples used to determine the bulk density (BD). A pycnometer is a device used to determine the density of a liquid. Such instrument is usually made of glass, with a close-fitting ground glass stopper with a capillary tube through it, so that air bubbles may escape from the apparatus. This device enables a liquid’s density to be measured accurately by reference to an appropriate working fluid, such as water or mercury, using an analytical balance. Particle density (PD) was calculated as follows (eqn. 3):

$$PD = \frac{dw \cdot Ws}{Ws - (Wsw - Ww)}$$

where dw is the density of water (g cm-3) at the temperature observed, Ws is the weight of soil sample (oven dry), Wsw is the weight of pycnometer, soil and water, and Ww is the weight of pycnometer and water.

## Statistical analysis

The experimental design was a factorial arrangement of treatments conducted in a completely randomized design. We applied general linear modeling (GLM) to relate bulk density, penetration resistance and total porosity to the slope gradient (three classes), machine passes (four classes), soil depth (four classes) and their interactions. Since no departure from normal distribution of the data was observed after Kolmogorov-Smirnov test (α = 0.05), standard parametric analyses were carried out. Homogeneity of variance among treatments was verified by the Levene’s test (α = 0.01). Post-hoc comparison of group means was performed using the Duncan’s multiple range test with a 99% confidence level. Treatment effects were considered significant when P<0.01.

All statistical analyses were carried using the SPSS (release 15.0) statistical package.

# Results

## Bulk density

Slope gradient (SG), traffic intensity (T) and soil profile depth (D - depth readings were referenced to the mineral soil surface), as well as the interaction effects of SG×T, SG×D, SG×D and SG×T×D significantly affected BD (Tab. 1 and Tab. 2). Average BD varied significantly among slope gradients and ranged from 1.20 g cm-3 on the 0-10% downhill forwarding slope, 1.35 g cm-3 on the 10-20% downhill forwarding slope and 1.46 g cm-3 on the 0-10% uphill forwarding slope. Average BD increased significantly with increasing frequency of machine passes and reached the highest level of 1.45 g cm-3 after 10 passes. Average BD also increased significantly with increasing soil profile depth. However, the major increases following forwarding occurred at the soil surface (0-5 cm), particularly on the 0-10% uphill forwarding slope and the 10-20% downhill forwarding slope.

Tab. 1 - Means ± standard deviation of different physical soil properties on skid trails after forwarding with a farm tractor. (BD): bulk density; (TP): total porosity; (PR): penetration resistance. Different letters indicate statistically significant differences within columns for each slope class separately (P<0.01), based on analysis of variance (GLM).

Slope
(%)
Traffic
(# of passes)
Depth
(cm)
BD
(g cm-3)
PR
(kPa)
TP
(%)
0-10 %
(down)
0 0-5 1.02 ± 0.01d 1367.7 ± 35.7b 55.3 ± 1.1a
5-10 1.09 ± 0.02c 1366.0 ± 24.2b 53.6 ± 4a
10-15 1.18 ± 0.01b 1382.7 ± 16.2ab 51.1 ± 2.7a
15-20 1.23 ± 0.01a 1429.3 ± 19a 45.3 ± 1.7b
<5 0-5 1.09 ± 0.01d 1487.3 ± 2.9a 53.5 ± 2.5a
5-10 1.17 ± 0.01c 1470.3 ± 44.7a 49.4 ± 1.6ab
10-15 1.21 ± 0.02b 1496.7 ± 86.6a 43.7 ± 6.5b
15-20 1.28 ± 0.01a 1562.3 ± 48.1a 44.8 ± 3.3b
5-10 0-5 1.07 ± 0.01d 1614.0 ± 45.6a 55.0 ± 1.6a
5-10 1.21 ± 0.02c 1553.7 ± 18.3ab 45.6 ± 0.7b
10-15 1.25 ± 0.01b 1541.3 ± 43b 44.4 ± 1.7b
15-20 1.33 ±0.01a 1498.0 ± 18.2b 40.0 ± 2.9c
>10 0-5 1.15 ± 0.02c 1595.7 ± 27.4b 48.8 ± 1.7a
5-10 1.25 ± 0.01b 1643.3 ± 20.8ab 47.1 ± 1.9a
10-15 1.34 ± 0.03a 1690.0 ± 45.8a 43.3 ± 2.2a
15-20 1.37 ± 0.04a 1695.3 ± 32.3a 44.3 ± 6.4a
10-20 %
(down)
0 0-5 1.03 ± 0.01d 1343.3 ± 32.1c 56.2 ± 2.4a
5-10 1.09 ± 0.01c 1366.7 ± 30.6c 53.8 ± 3.7ab
10-15 1.19 ± 0.01b 1433.3 ± 28.9b 50.7 ± 2.3bc
15-20 1.24 ± 0.01a 1486.7 ± 5.8a 47.2 ± 2.4c
<5 0-5 1.23 ± 0.02d 1672.0 ± 41.1b 47.6 ± 2.4a
5-10 1.37 ± 0.01c 1676.7 ± 11.6b 41.2 ± 1.7b
10-15 1.43 ± 0.02b 1717.7 ± 42.4ab 37.5 ± 1.8b
15-20 1.47 ± 0.03a 1778.7 ± 51.4a 38.2 ± 4.6b
5-10 0-5 1.29 ± 0.02d 1763.3 ± 36.1a 45.8 ± 6.4a
5-10 1.41 ± 0.02c 1776.7 ± 38.3a 43.3 ± 1.1a
10-15 1.45 ± 0.01b 1752.0 ± 18.4a 40.4 ± 1.3ab
15-20 1.52 ± 0.01a 1807.7 ± 51.9a 33.6 ± 3.9b
>10 0-5 1.36 ± 0.03c 1895.3 ± 49.5a 41.2 ± 3.0a
5-10 1.45 ± 0.01b 1943.7 ± 37.6a 39.1 ± 5.4a
10-15 1.51 ± 0.02a 1916.7 ± 30.6a 34.0 ± 7.8a
15-20 1.52 ± 0.01a 1953.3 ± 15.3a 35.5 ± 3.0a
0-10 %
(up)
0 0-5 1.04 ± 0.02d 1376.7 ± 20.8b 57.8 ± 0.9a
5-10 1.10 ± 0.01c 1356.7 ± 32.1ab 50.8 ± 1.0b
10-15 1.17 ± 0.02b 1380.0 ± 17.3ab 48.7 ± 1.6bc
15-20 1.24 ± 0.01a 1423.3 ± 25.2a 47.1 ± 3.0c
<5 0-5 1.41 ± 0.02a 1577.0 ± 28.6c 38.1 ± 5.5a
5-10 1.50 ± 0.02b 1674.0 ± 50.5b 30.8 ± 4.4ab
10-15 1.57 ± 0.01a 1727.3 ± 34.1b 35.8 ± 2.5ab
15-20 1.58 ± 0.02a 1826.7 ± 23.1a 29.9 ± 2.7b
5-10 0-5 1.49 ± 0.02c 1953.3 ± 15.3c 34.0 ± 1.2a
5-10 1.58 ± 0.01b 2067.3 ± 50.0b 30.1 ± 2.2a
10-15 1.61 ± 0.03ab 2078.3 ± 44.1b 26.6 ± 2.7ab
15-20 1.64 ± 0.02a 2166.7 ± 54.6a 30.7 ± 1.8b
>10 0-5 1.55 ± 0.02b 2197.7 ± 81.0a 34.1 ± 1.8a
5-10 1.63 ± 0.02a 2152.7 ± 10.5a 34.8 ± 4.7a
10-15 1.63 ± 0.02a 2196.7 ± 70.1a 30.8 ± 3.9a
15-20 1.65 ± 0.02a 2193.0 ± 37.4a 30.2 ± 4.5a

Tab. 2 - Analysis of variance (ANOVA) for the effects of slope gradient (SG), traffic intensity (T), soil profile depth (D) classes and their interactions on BD after forwarding with a farm tractor.

Source Sum of
Squares
df Mean
Square
F-value P-value
SG 1.61 2 0.81 3224.34 < 0.001
T 2.11 3 0.7 2811.7 < 0.001
D 0.79 3 0.26 1049.24 < 0.001
SG × T 0.51 6 0.08 339.78 < 0.001
SG × D 0.02 6 0 12.31 < 0.001
T × D 0.02 9 0 8.2 < 0.001
SG × T × D 0.02 18 0 3.92 < 0.001

The effects of significant interactions of slope gradient with traffic intensity and soil profile depth on BD are shown in Fig. 2. Average BD increased with increasing traffic intensity, regardless of slope gradient, but throughout the profile BD increase was particularly strong in plots with fewer than 5 passes on the 0-10% uphill and 10-20% downhill forwarding slope compared with control plots. The highest average BD values were observed on the 0-10% uphill forwarding slope with >10 machine passes (Tab. 1), which also experienced the greatest amount of relative change (Fig. 3). Further, average BD increased with increasing soil profile depth in each slope gradient class (Tab. 1), but relative changes were greatest in the upper 10 cm than of the soil profile (Fig. 2). The highest average BD values were observed on the 0-10% uphill forwarding slope at a soil depth of 15-20 cm (Tab. 1).

Fig. 2 - Change in BD (%) after forwarding by slope gradient class, traffic intensity class and soil profile depth class. Changes are relative to measurements from the same profile depth class in untrafficked control plots.

Fig. 3 - Change in PR (%) after forwarding by slope gradient class, traffic intensity class and soil profile depth class. Changes are relative to measurements from the same profile depth class in untrafficked control plots.

## Penetration resistance

Slope gradient (SG), traffic intensity (T) and soil profile depth (D), as well as the interaction effects of SG×T, SG×D, SG×D and SG×T×D significantly affected PR (Tab. 3). Average PR by slope gradient was 1525 kPa on the 0-10% downhill, 1705 kPa on the 10-20% downhill and 1834 kPa on the 0-10% uphill forwarding slope. Average PR increased significantly with increasing the frequency of machine passes and reached the highest level of 1922 kPa after 10 passes. Average PR also increased significantly with increasing soil profile depth.

Tab. 3 - Analysis of variance (ANOVA) for the effects of slope gradient (SG), traffic intensity (T), soil profile depth (D) classes and their interactions on PR after forwarding with a farm tractor.

Source Sum of
Squares
df Mean
Square
F-value P-value
SG 2321852 2 1160926 860.29 < 0.001
T 5643747 3 1881249 1394.08 < 0.001
D 134033.1 3 44677.71 33.11 < 0.001
SG × T 1286932 6 214488.7 158.95 < 0.001
SG × D 29403.57 6 4900.59 3.63 < 0.001
T × D 38066.81 9 4229.65 3.13 < 0.001
SG × T × D 110708.2 18 6150.45 4.56 < 0.001

Average PR increased significantly with increasing traffic intensity, regardless of slope gradient, with particularly strong PR increases on the 10-20% downhill and the 0-10% uphill forwarding slope compared to control plots; the increase in PR on the trafficked 0-10% downhill forwarding slope was lower at all traffic intensities and soil depths (Fig. 3). Compared to undisturbed control areas, which had PR values of 1386 kPa on the 0-10% downhill, 1408 kPa on the 10-20% downhill and 1384 kPa on the 0-10% uphill forwarding slope, average PR values after >10 machine passes increased to 1656 kPa on the 0-10% downhill, 1927 kPa on the 10-20% downhill and 2185 kPa on the 0-10% uphill forwarding slope. Soil penetration resistance continually increased with traffic frequency in each soil depth and was still increasing at >10 passes, particularly on the 0-10% uphill and 10-20% downhill forwarding slope (Fig. 3). Although the highest average PR values were found on the 0-10% uphill forwarding slope at a soil depth of 10-20 cm, PR did not increase significantly with soil depth in any of the slope gradients.

## Total porosity

Slope gradient (SG), traffic intensity (T) and soil profile depth (D) as well as the interaction of SG×T significantly affected TP (Tab. 4). Average TP by slope gradient was 47.8% on the 0-10% downhill, 42.8% on the 10-20% downhill and 36.9% on the 0-10% uphill forwarding slope (Fig. 4). Average TP decreased significantly with increasing the traffic intensity and reached the lowest level of 38.6% after 10 passes. However, differences in TP were small among trafficked plots. TP also decreased significantly with increasing soil profile depth (Tab. 1).

Tab. 4 - Analysis of variance (ANOVA) for the effects of slope gradient (SG), traffic intensity (T), soil profile depth (D) classes and their interactions on TP after forwarding with a farm tractor.

Source Sum of
Squares
df Mean
Square
F-value P-value
SG 2873.63 2 1436.81 126.25 < 0.001
T 3955.4 3 1318.47 115.85 < 0.001
D 1443.86 3 481.29 42.29 < 0.001
SG × T 1013.74 6 168.96 14.85 < 0.001
SG × D 61.63 6 10.27 0.9 0.50
T × D 128.83 9 14.31 1.26 0.27
SG × T × D 281.26 18 15.63 1.37 0.16

Fig. 4 - Change in TP (%) after forwarding by slope gradient class, traffic intensity class and soil profile depth class. Changes are relative to measurements from the same profile depth class in untrafficked control plots.

Average TP was highest in control plots and decreased with increasing traffic intensity, regardless of slope gradient (Tab. 1). Particularly strong TP decreases were observed in plots with fewer than 5 passes on the 0-10% uphill and 10-20% downhill forwarding slope as compared to control plots. Generally, only small additional decreases in TP were observed with additional passes (Fig. 4). The lowest average TP values were observed on the 0-10% uphill forwarding slope after the first 5 passes (29.9%), where TP was 17.5% lower compared to the controls (Tab. 1).

# Discussion

In this study, the investigated physical soil properties exhibited consistent and largely predictable responses to different tractor traffic intensities throughout the upper 20 cm of the soil profile. The significant impact of traffic intensity on BD, PR and TP on all slope gradients is consistent with findings from previous studies that have documented a marked increase of bulk density, soil strength and loss of macropores with increasing traffic intensity ([38], [16], [12], [40], [4], [31]). The large initial increase in compaction with fewer than 5 machine passes and the modest increase thereafter reflects the typical logarithmic relationship between the degree of soil compaction and the traffic intensity, where the greatest impacts occur during the first three passes and traffic intensities beyond 10 passes contribute very little to further compaction ([18], [16], [12], [13], [20], [22], [7], [4], [31]). One explanation for this logarithmic relationship is that the bearing capacity of the soil (maximum load without soil structure failure) grows with increasing BD, protecting the surface layer against further compaction with increasing traffic ([7], [4], [31]). The increase in PR beyond 10 passes observed in this study is in agreement with findings that PR does not taper off after a few passes but continues to increase with traffic frequencies well beyond the traffic intensity seen in this study ([46]). PR is a more sensitive measure of logging traffic intensity than bulk density, because particle rearrangement after soil disturbance can increase PR but not bulk density ([15]) and greater increases in PR than bulk density are commonly found for the same logging intensity ([2], [8]). The maximum PR increase of up to 60% observed in this study is similar to the 50% increase reported by Picchio et al. ([31]), but maximum PR values of 2.3 MPa did not reach levels where root growth of trees becomes restricted (e.g., 2.5-3 MPa; [18], [7]).

The reduction in TP of up to 45% in this study is similar to that reported by previous analyses, and is likely due to a loss or substantial decrease (50-60%) of macropores in favor of smaller pores after the first few machine passes ([4], [31]). In the soil compaction process, soil pores are typically compressed or destroyed, surface aggregates are broken down, and macropores are transformed into meso- (0.2<d<50 µm) and micropores (d<0.2 µm). As substantial destruction of macropores can already occur after the first few passes, this would explain why additional machine traffic only minimally decreased TP beyond the initial reduction.

Even though the tendency of increased BD and PR and decreased TP with increased traffic intensity was observed in all slope classes, the observed changes were highest when the forwarding direction was uphill (even on moderate slopes and after just a few passes) or on steeper slopes when the forwarding direction was downhill. Differences in levels of compaction between uphill and downhill forwarding may be explained by an uneven load distribution between the downhill and uphill tires of the skidder ([22]) or increased wheel slippage and vibration in uphill compared to downhill forwarding. Because stresses applied at the soil surface are always transmitted three-dimensionally, they cause not only soil compaction but also shear effects ([20]) and exposure of mineral subsoil, which naturally has a higher BD and PR than the surface layer. Further, slope gradient has been found to be an important factor for the potential level of compaction ([38]), presumably because machines can slip continuously and remain longer in a given spot in steeper terrain, thus puddling and dragging the soil ([16]).

An increase in BD and PR and a decrease in TP with soil depth, as observed in control plots on all slope classes in this study, is typical for soil profiles. In untrafficked forest soils, BD and PR are typically low and TP is high in the surface horizons due to large amounts of organic matter; as the organic matter content decreases with soil depth, BD and PR increase and TP decreases ([17], [13]). Further, biological activity of roots and animals can reduce BD and PR and increase TP in the upper soil, while fine soil texture, gravel content and soil structure may increase BD and PR and lower TP at lower soil depths ([18], [4], [31]). The pore volume of a medium- to fine-textured soil consists mainly of meso- and micropores that easily hold water against gravitational forces. The smaller the pore diameter, the more the adhesion of water molecules to soil colloids exceeds the force of gravitational capillarity. So in a saturated state, all pores are filled with water that cannot be compressed. Soil moisture content determines the proportion between soil compaction and plastic deformation (in the form of rutting) after the application of machine forces ([3]). Tractor traffic with its associated wheel or track slip can directly affect the soil structure and alter physical soil properties down to deeper depths. In this study, increased soil compaction occurred throughout the soil profile with associated increases in BD and PR and reduction of TP, which has also been observed previously ([18], [38], [16], [25], [30], [5], [13]). However, relative to values from different soil depths in the control plots, changes in BD, PR and TP following machine traffic were often (but not always) highest in the upper 5 cm layer and lowest in the 20 cm layer of the soil profile. Compaction effects are generally most distinct in the upper soil layer, since the exerted pressure is maximal at the soil surface and declines with increasing depth as the total pressure is spread out over a growing area ([18], [4], [31]). As soils are compacted and lose large pores after a few passes, compaction leads to higher soil strength in the upper soil layers, which prevents much further compaction of the surface layer with additional passes and shifts forces to deeper soil layers that are then increasingly compacted ([4], [31]). Overall, damage to the soil (including deeper soil layers) was modest in this study, which might be bacause the intensity of machine traffic was generally at the low end compared to other studies (e.g., [46]) and the very dry soil conditions on this site at the time of the forwarding operation were better able to resist to the mechanical soil compaction and rutting (e.g., [43], [7]).

# Conclusion

This study clearly demonstrates that changes in physical soil properties in response to farm tractor traffic are significantly affected by slope gradient and forwarding direction (i.e., uphill vs. downhill forwarding). Uphill forwarding, even on gentle slopes below 10% inclination, consistently caused more soil damage throughout the upper 20 cm of the soil profile than downhill forwarding on even steeper slopes. The pattern of the largest negative impact for uphill forwarding was the same for all three soil response variables, bulk density, resistance to penetration and total porosity, and at all traffic intensities along the skid trail. As a consequence, we recommend that skid trails have to be designed to avoid uphill forwarding by this farm tractor-2-wheel-trailer system. If possible, the application of this forwarding system should be limited to gentle slopes and alternative means of forwarding should be sought on steeper slopes. Given that impacts would have been even more severe had soils not been dry during the forwarding operation, logging operations should be restricted to periods when soil conditions are dry in order to minimize compaction and rutting. Future designs of skid trails should not only take into account the gradient of the slope of a skid trail, but also the direction of forwarding used by the chosen harvesting system.

# Acknowledgments

The authors acknowledge the financial support by the University of Tehran for this research under grant number 28514/1/3.

# References

(1)
Adams PW, Froehlich HA (1984). Compaction of forest soils. Extension Publication PNW 217, Pacific Northwest Research Station, USDA Forest Service, Portland, OR, USA, pp. 13.
Gscholar
(2)
Alban DH, Host GE, Elioff JD, Shadis D (1994). Soil and vegetation response to soil compaction and forest floor removal after aspen harvesting. Research Paper RP-NC-RP-315, North Central Forest Experiment Station, USDA Forest Service, St. Paul, MN, USA, pp. 8.
Online | Gscholar
(3)
Ampoorter E, de Schrijver A, Van Nevel L, Hermy M, Verheyen K (2012). Impact of mechanized harvesting on compaction of sandy and clayey forest soils: results of a meta-analysis. Annals of Forest Science 69: 533-542.
CrossRef | Gscholar
(4)
Ampoorter E, Van Nevel L, De Vos B, Hermy M, Verheyen K (2010). Assessing the effects of initial soil characteristics, machine mass and traffic intensity on forest soil compaction. Forest Ecology and Management 260: 1664-1676.
CrossRef | Gscholar
(5)
Ares A, Terry TA, Miller RE, Anderson HW, Flaming BL (2005). Ground-based forest harvesting effects on soil physical properties and douglas-fir growth. Soil Science Society of America Journal 69: 1822-1832.
CrossRef | Gscholar
(6)
Bock MD, Van Rees KCJ (2002). Mechanical site preparation impacts on soil properties and vegetation communities in the Northwest Territories. Canadian Journal of Forest Research 32: 1381-1392.
CrossRef | Gscholar
(7)
Bolding MC, Kellogg LD, Davis CT (2009). Soil compaction and visual disturbance following an integrated mechanical forest fuel reduction operation in Southwest Oregon. International Journal of Forest Engineering 20 (2): 47-56.
Gscholar
(8)
Brais S (2001). Persistence of soil compaction and effects on seedling growth in northwestern Quebec. Soil Science Society of America Journal 65: 1263-1271.
CrossRef | Gscholar
(9)
Bustos O, Egan A (2011). A comparison of soil compaction associated with four ground-based harvesting systems. Northern Journal of Applied Forestry 28 (4): 194-198.
Online | Gscholar
(10)
Conway S (1984). Logging practice, principles of timber harvesting systems. Miller Freeman Publications, San Francisco, USA, pp. 465.
Gscholar
(11)
Dykstra DP, Heinrich R (1996). FAO model code of forest harvesting practice. FAO, Rome, Italy, pp. 85.
Gscholar
(12)
Eliasson L (2005). Effects of forwarder tire pressure on rut formation and soil compaction. Silva Fennica 39: 549-557.
Online | Gscholar
(13)
Eliasson L, Wasterlund I (2007). Effects of slash reinforcement of strip roads on rutting and soil compaction on a moist fine-grained soil. Forest Ecology and Management 252: 118-123.
CrossRef | Gscholar
(14)
Ezzati S, Najafi A, Rab MA, Zenner EK (2012). Recovery of soil bulk density, porosity and rutting from ground skidding over a 20-year period after timber harvesting in Iran. Silva Fennica 46 (4): 521-538.
Online | Gscholar
(15)
Froehlich HA, McNabb DH (1984). Minimizing soil compaction in Pacific Northwest forests. In: “Forest soils and treatment impacts” (Stone EL ed). Proceedings of the “6th North American forest soils conference”. Soil Science Society of America, Madison, WI, USA, pp. 159-192.
Gscholar
(16)
Gayoso J, Iroume A (1991). Compaction and soil disturbances from logging in Southern Chile. Annals of Forest Science 48: 63-71.
CrossRef | Gscholar
(17)
Gent JA, Morris LA (1986). Soil compaction from harvesting and site preparation in the upper Gulf Coastal Plain. Soil Science Society of America Journal 50: 443-446.
CrossRef | Gscholar
(18)
Greacen EL, Sands R (1980). Compaction of forest soil, a review. Australian Journal of Soil Research 18: 163-189.
CrossRef | Gscholar
(19)
Greulich FR, Hanley DP, McNeel JF, Baumgartner D (1999). A primer for timber harvesting. EB1316, WSU Cooperative Extension, Pullman, Washington, USA, pp. 33.
Gscholar
(20)
Horn R, Vossbrink J, Peth S, Becker S (2007). Impact of modern forest vehicles on soil physical properties. Forest Ecology and Management 248: 56-63.
CrossRef | Gscholar
(21)
Ilstedta U, Malmerb A, Nordgrenb A, Liauc P (2004). Soil rehabilitation following tractor logging: early results on amendments and tilling in a second rotation Acacia mangium plantation in Sabah, Malaysia. Forest Ecology and Management 194: 215-222.
CrossRef | Gscholar
(22)
Jamshidi R, Jaeger D, Raafatnia N, Tabari M (2008). Influence of two ground-based skidding systems on soil compaction under different slope and gradient conditions. International Journal of Forest Engineering 19: 9-16.
Online | Gscholar
(23)
Jourgholami M (2012). Small-scale timber harvesting; mule logging in Hyrcanian Forest. Small-scale Forestry 11 (2): 255-262.
CrossRef | Gscholar
(24)
Krag R, Higgingbotham K, Rothwell R (1986). Logging and soil disturbance in southeast British Columbia. Canadian Journal of Forest Research 16: 1345-1354.
CrossRef | Gscholar
(25)
Landsberg JD, Miller RE, Anderson HW, Tepp JS (2003). Bulk density and soil resistance to penetration as affected by commercial thinning in northeastern Washington. Research Paper PNW-RP-551, Pacific Northwest Research Station, USDA Forest Service, Portland, OR, USA, pp. 35.
Online | Gscholar
(26)
McNabb DH, Startsev AD, Nguyen H (2001). Soil wetness and traffic level effects on bulk density and air-filled porosity of compacted boreal forest soils. Soil Science Society of America Journal 65: 1238-1247.
CrossRef | Gscholar
(27)
Miller D, Anderson H (2002). Soil compaction: concerns, claims, and evidence. In: Proceedings of the Meeting “Small Diameter Timber: Resource Management, Manufacturing, and Markets”. Spokane (Washington, USA) 25-27 Feb 2002. Washington State University Cooperative Extension, Friday Harbor, WA, USA, pp. 97-106.
Gscholar
(28)
Mousavi R (2009). Comparison of productivity, cost and environmental impacts of two harvesting methods in northern Iran: short-log vs. long-log. PhD thesis, Faculty of Forest Sciences, University of Joensuu, Finland, pp. 93.
Gscholar
(29)
Murphy G, Firth JG, Skinner MF (2004). Long-term impacts of forest harvesting related soil disturbance on log product yields and economic potential in a New Zealand forest. Silva Fennica 38 (3): 279-289.
Online | Gscholar
(30)
Nugent C, Kanali C, Owende PMO, Nieuwenhuis M, Ward S (2003). Characteristic site disturbance due to harvesting and extraction machinery traffic on sensitive forest sites with peat soils. Forest Ecology and Management 180: 85-98.
CrossRef | Gscholar
(31)
Picchio R, Neri F, Petrini E, Verani S, Marchi E, Certini G (2012). Machinery-induced soil compaction in thinning two pine stands in central Italy. Forest Ecology and Management 285: 38-43.
CrossRef | Gscholar
(32)
Rab MA (1996). Soil physical and hydrological properties following logging and slash burning in the Eucalyptus regnans forest of southeastern Australia. Forest Ecology and Management 84: 159-176.
CrossRef | Gscholar
(33)
Rodriguez EO (1986). Wood extraction with oxen and agricultural tractors. FAO Forestry Paper No. 49, Rome, Italy, pp. 120.
Online | Gscholar
(34)
Rohand K, Kalb AA, Herbauts J, Verbrugge JC (2004). Changes in some mechanical properties of a loamy soil under the influence of mechanized forest exploitation in a beech forest of central Belgium. Journal of Terramechanics 40: 235-253.
CrossRef | Gscholar
(35)
Russell F, Mortimer D (2005). A review of small-scale harvesting systems in use worldwide and their potential application in Irish forestry. COFORD, Dublin, Ireland, pp. 48.
Online | Gscholar
(36)
Savelli S, Cavalli R, Baldini S, Picchio R (2010). Small scale mechanization of thinning in artificial coniferous plantation. Croatian Journal of Forest Engineering 31 (1): 11-21.
Online | Gscholar
(37)
Shaffer RM (1992). Farm tractor logging for woodlot owners. College of Natural Resources, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA, pp. 5.
Gscholar
(38)
Sidle RC, Drlica DM (1981). Soil compaction from logging with a low-ground pressure skidder in the Oregon Coast Ranges. Soil Science Society of America Journal 45: 1219-1224.
CrossRef | Gscholar
(39)
Spinelli R, Magagnotti N (2012). Wood extraction with farm tractor and sulky: estimating productivity, cost and energy consumption. Small-scale Forestry 11 (1): 73-85.
CrossRef | Gscholar
(40)
Susnjar M, Horvat D, Seselj J (2006). Soil compaction in timber skidding in winter conditions. Croatian Journal of Forest Engineering 27: 3-15.
Online | Gscholar
(41)
Updegraff K, Blinn CR (2000). Applications of small-scale forest harvesting equipment in the United States and Canada. Staff Paper Series no. 143, College of Natural Resources and Minnesota Agricultural Experiment Station, University of Minnesota, St. Paul, MN, USA, pp. 51.
Online | Gscholar
(42)
Wang J, LeDoux CB, Edwards P (2007). Changes in soil bulk density resulting from construction and conventional cable skidding using preplanned skid trails. Northern Journal of Applied Forestry 24: 5-8.
Gscholar
(43)
Williamson JR, Neilsen WA (2000). The influence of forest site on rate and extent of soil compaction and profile disturbance of skid trails during ground-based harvesting. Canadian Journal of Forest Research 30: 1196-1205.
CrossRef | Gscholar
(44)
Wronski EB, Murphy G (1994). Responses of forest crops to soil compaction. In: “Soil Compaction in Crop Production” (Soane BD, van Ouwerkerk C eds). Elsevier, Amsterdam, The Netherlands, pp. 317-342.
Gscholar
(45)
Zenner EK, Berger AL (2008). Influence of skidder traffic and canopy removal intensities on the ground flora in a clearcut-reserves northern hardwood stand in Minnesota, USA. Forest Ecology and Management 256: 1785-1794.
CrossRef | Gscholar
(46)
Zenner EK, Fauskee JT, Berger AL, Puettmann KJ (2007). Impacts of skidding traffic intensity on soil disturbance, soil recovery, and aspen regeneration in north central Minnesota. Northern Journal of Applied Forestry 24: 177-183.
Online | Gscholar

#### Authors’ Affiliation

(1)
Shoukuh Soltanpour
Department of Forestry and Forest Economics, Faculty of Natural Resources, University of Tehran, Zob-e-Ahan Street, 315854314 Karaj (Iran)
(2)
Eric K Zenner
College of Agricultural Sciences, Pennstate University, 305 Forest Resources Building, University Park, 16802 PA (United States of America)

#### Corresponding author

mjgholami@ut.ac.ir

#### Citation

Jourgholami M, Soltanpour S, Etehadi Abari M, Zenner EK (2014). Influence of slope on physical soil disturbance due to farm tractor forwarding in a Hyrcanian forest of northern Iran. iForest 7: 342-348. - doi: 10.3832/ifor1141-007

#### Paper history

Accepted: Nov 25, 2013

First online: Apr 17, 2014
Publication Date: Oct 01, 2014
Publication Time: 4.77 months

© SISEF - The Italian Society of Silviculture and Forest Ecology 2014

#### Breakdown by View Type

(Waiting for server response...)

#### Article Usage

Total Article Views: 13639
(from publication date up to now)

Breakdown by View Type
HTML Page Views: 9787
Abstract Page Views: 643

Web Metrics
Days since publication: 2185
Overall contacts: 13639
Avg. contacts per week: 43.69

Article citations are based on data periodically collected from the Clarivate Web of Science web site
(last update: Jan 2020)

Total number of cites (since 2014): 14
Average cites per year: 2.00

#### iForest Database Search

Search By Author

Search By Keyword

Citing Articles

Search By Author

Search By Keywords

#### PubMed Search

Search By Author

Search By Keyword