Modern wood forwarding operations require heavy machines to move across sensitive terrain and soil conditions. In this regard, improper timing by season and weather can lead to substantive soil compression and rutting across and along flow channels, on ridge tops, through depressions with moist to wet soils. As soils become compacted and/or displaced, studies have shown that:
1) Soil porosity is reduced affecting the oxygen levels and gas distributions within the soil ( Czyz , 2004
2) Increasing soil bulk density affects plant root distribution and decrease hydraulic conductivity thereby reducing plant available water while enhancing runoff and soil erosion (Horn et al., 1995; Startsev & McNabb, 2000; Jamshidi et al., 2008) ;
4) Soil compaction has lasting effects on soil, with natural soil-structure recovery varying from a few years to decades ( Prose, 1985 ; Brady & Weil, 2008 ; Labelle & Jaeger, 2011 ; Ezzati et al., 2012) .
This article focuses on monitoring and examining machine-to-ground clearances in relation to number of passes, wood forwarding machine speed, and dry to wet ground conditions involving 54 harvest blocks in New Brunswick, Canada. For this purpose, two forwarders (John Deere 1110E and 1510E) and one grapple skidder (Tigercat 635D) were equipped with ultrasonic distance sensors, and GPS data loggers. The wood-forwarding operations involved retrieving logs from clear cuts, shelter-wood cuts, and commercially thinned forest plantations. The thinning operations involved laying out brushmats to reduce soil compression impacts (Labelle & Jaeger, 2012) . The dry to wet variations in ground conditions were related to changing seasons and weather, and to topography across each harvest block as revealed by way of the metric cartographic depth-to-water index ( Murphy et al., 2011 ; White et al., 2012) .
2.1. Site Description
This harvest-block study was spread across four ecoregions in New Brunswick, focusing on a North-western area near Saint-Quentin, and amid-western area centralized near Juniper and Dorn Ridge, as described below (Figure 1). For this study, the areas were split into three groups―Northwestern uplands (NWU), Midwestern uplands (MWU), and lowlands (LL).
The NWU study area (elevation range 230 to 450 m)is located on the Chaleur uplands in northwest New Brunswick, north of Grand Falls and encompasses the both the highlands and the northern part of the Southern Uplands ecoregion. The forest cover includes tolerant hardwoods, mixed woods and forest plantations (spruce species). Dominant tree species refer to sugar maple (Acer saccharum Marsh.), balsam fir (Abies balsamea L.), yellow birch (Betula alleghaniensis Britt.), red spruce (Picea rubens Sarg.), black spruce (Picea mariana Mill.) and white spruce (Picea glauca (Moench) Voss). Mean annual air temperature from
1990-2016 amounted to 3.6˚C, with mean annual January and July air temperatures at −5.3˚C and 12.5˚C. Mean annual precipitation was 1140 mm, with 310 mm as snow (Department of Environment and Climate Change Canada, 2016) . BB is underlain by Ordovician-Silurian carbonates and Ordovician-Silurian-Devonian deep-water clastics. The terrain varies from rolling to hummocky, interspersed by steeply incised valleys. Surficial deposits vary from residuals to stony ablation and loamy lodgement tills, glaciofluvial deposits (moraines, kames, eskers), and alluvium.
The MWU study range (elevations 250 to 580 m) is located in mid-western New Brunswick, near Juniper and encompasses the lower part of the Southern Uplands ecoregion. The forest cover consists of tolerant hardwoods with sugar maple and white birch (Betula papyrifera Marsh), mixed with balsam fir and black spruce, and interspersed by black and white spruce forest plantations. Mean annual air temperature from 1990-2016 amounted to 5.2˚C, with mean annual January and July air temperatures at −3.5˚C and 13.9˚C. Mean annual precipitation was1180 mm, with 280 mm as snow (Department of Environment and Climate Change Canada, 2016) .
The LL study range (elevations 120 to 290 m) is located around 50 km northwest of Fredericton and covers both the Continental Lowlands and Eastern Lowlands ecoregions. This area is predominantly covered by tolerant hardwoods consisting of yellow birch, beech (Fagus grandifolia Ehrh.), and sugar maple, mixed with balsam fir, Eastern white cedar (Thuja occidentalis L.) and black spruce. The mean annual air temperature from 1990-2016 amounted to 5.5˚C, with mean annual January and July air temperatures at −2.8˚C and 13.8˚C. Mean annual precipitation amounted to 1100 mm, with 250 mm as snow (Department of Environment and Climate Change Canada, 2016) . Both MWU and LL are located on the Miramichi Caledonia highlands at the northeast stretch of the Appalachian Mountain range. Surficial deposits on rolling to moderate terrain vary from bouldery loamy lodgement till to moraines, eskers, kames and sandy glaciofluvial outwash plains. Bedrock formations on LL mainly refer to Cambrian-Ordovician-Silurian deep-water clastics, and on MWU mainly Carboniferous to Silurian felsic to mafic extrusions.
2.2. Machinery and Sensor Installations
In NWU, two John Deere (JD) forwarders, i.e., model JD 1110E equipped with tire chains (referred to as JD11 below), and model JD 1510E with front and back tire tracks (referred to as JD15 below) were used for tracking wood-forwarding machine clearances (Table 1, Figure 2). In MWU and LL, a Tigercat model TC 635Dgrapple skidder (referred to as TC below) was used, with chained tires in front and tracked tires in the rear. Machine clearances were tracked using custom-built data loggers (c/o FP Innovations Ltd., Montreal, Canada, Figure 3) to record time-stamped GPS locations and machine-to-ground clearances during back and forth machine travel within up to 54 harvesting blocks from February 2012 through November 2014, all at time-stamped 10-second intervals. The data loggers were installed within the operator cabs, with polycarbonate GPS antennas firmly attached to the cab roofs. Clearance tracking involved two ultrasonic sensors, each vertically placed in an open steel pipe welded or bolted to the least vulnerable position on the outside chassis of each machine, with sensor wires safely guided to the data loggers. On the two forwarders, the ultrasonic sensors were mounted on the right and left of the chassis (145 cm for JD11and 140 cm for JD15 from the ground up to sensor). On the TC, the sensors were mounted on the front and back of the chassis (170 cm and 152 cm from the ground to sensor, respectively). The sensors signalled machine clearances to
Table 1. Machine specifications.
ground, stumps, protruding rocks and brush piles directly below the sensors, adjacent to the tracks (Figure 4).
Operations involved JD11 and JD15 wood forwarding from commercially thinned forest plantations in northwest NB, mostly centred on NWU, and TC wood forwarding from shelter-wood and clear cutting operations in central west NB, mostly centered on MWU and LL. The commercial thinning operations involved laying out brushmats (Labelle & Jaeger, 2012) . The shelter-wood and clear cutting operations proceeded on bare ground covered by organic forest floor accumulations about 5 to 10 cm thick. JD15 operations were clearance
Figure 2. Left; wood-forwarding equipment: JD 1110E (top), JD 1510E (middle), and TC 635D (bottom). Right; wood harvesting-forwarding operations: commercial thinning (white spruce plantation, top), clearcutting (middle), hardwood shelter-wood cutting (bottom).
Figure 3. Left: Ultrasonic protective mount for the ultrasonic distance sensors. Right: datalogger box connected to GPS and clearance sensors.
Figure 4. Various sensor distance measurements and error potential from the ultrasonic distance sensor. Left to right: rutting due to soil compaction, protruding rocks, stumps, branches and slash, and soil displacement.
monitored year-round, including operations on snow accumulations. The JD11 and TC operations were monitored spring to fall.
2.3. Data Production
The logged data generated over 4.1 million data points across 54 harvesting blocks (Table 2). Data due to sensor malfunction, fixed distance-to-ground obstructions, machine idling, or traveling on paved and unpaved roads were removed. The remaining data were entered to determine point-for-point machine elevation (m), direction, slope (m/m), and speed (m/sec) along each track. Data were catalogued for each forwarding (loaded) and returning (unloaded) pass towards the loading zones. The number of passes per same forwarding and returning track was also determined. All data were processed through ArcGIS, which also included determining track densities and number of passes per same track using point buffering and overlapping tools (Buja, 2012) . The sensor-to-ground distance data (S) were normalized relative to the sensor-to-ground distance on solid ground (S0), referred to below as normalized sensor clearances (S/S0), or normalized clearances for short. These numbers can be converted into actual machine clearances (M) by setting M = S − (S0 − M0), where M0 is the machine-to-solid-ground (chassis) clearance (Table 1; Figure 4). Hence, at zero machine clearance (M = 0), S/S0 = 1 − M0/S0. For tracks with S < S0 − M0, machines bottoms would have sunk below the signal reflecting surface, as would be the case when the machines grind into wet soil or move across deep snow accumulations.
2.4. DTW Delineation
Each data point was placed into its geospatial elevation context using digital elevation models (DEMs) and associated depth-to-water maps (DTW, Figure 5; Murphy et al., 2009 ). The DTW map delineates the extent of the least elevation rise next to the nearest water bodies such as streams, rivers, and lakes (Murphy et al., 2009; White et al., 2012) . The DTW map, when derived for the end-of-summer water level for water bodies emulates soil drainage across the mapped areas from very poor (DTW < 10 cm), to poor (10 < DTW < 25 cm), imperfect (25 < DTW < 50 cm), moderate (50 < DTW < 100 cm), well (1 < DTW < 20 m) and excessive (DTW > 20 m) (Murphy et al., 2009) . The end-of-summer DTW condition was emulated using 4 ha of upslope flow accumulation area for
Table 2. Block description by machine, forest, operation, operation date, and ocular DTW % area coverage.
SW: Softwood, HW: Hardwood, MX: Mixedwood, Unk: Unknown, BF: Balsam Fir, WS: White Spruce, BS: Black Spruce, NS: Norway Spruce, RS: Red Spruce, WB: White Birch, BI: Birch, BE: Beech, SM: Sugar Maple, MA: Maple, PO: Poplar, CC: Clear cut, CT: Commercial Thinning; SHW: Shelterwood cut.
permanent stream flow initiation. For visualizing the varying DTW extent away from local flow channels by weather and season, ephemeral stream flow was initiated using 1 and 0.25 ha of upslope flow accumulation areas.
2.5. Hydrological Modelling
For contextual evaluation purposes, daily temporal variations in upland soil moisture, snowpack depth and frost depth were simulated for the NWU and for the combined MWU and LL areas using weather records for daily precipitation (rain, snow) and air temperature (daily means) from September 2011 to end of 2014 from Edmundston and Fredericton weather stations (Department of Environment and Climate Change Canada, 2016) (Appendix A3-1). Soil properties were set as follows: NWU = soil depth 1.5 m, loam to sandy loam, organic matter content 5% - 1%, coarse fragment content 20% - 30%; forest floor depth 10 cm; MWU = soil depth 1.2 m, sandy loam, organic matter content 1% - 10%, coarse fragment 25%, forest floor depth 5 cm. LL = soil depth 1.0 m, loamy sand, organic matter content 1% - 5%, coarse fragment 30%, forest floor depth 10 cm.
2.6. Data Processing
For geospatial visualization purposes, the individual data points per harvest block were mapped on top of the hill-shaded bare-earth DEM, with DTW drainage classes overlaid. The normalized clearance data were compiled into a single spreadsheet to enable the machine clearance analysis in relation to: 1) month of harvesting, 2) weather-induced soil wetness, 3) machine speed, 4) machine load, 5) back and forth track direction, 6) number of passes, 7) machine type, 8) harvesting type (clear cutting, commercial thinning, shelter-wood cutting), 9) slope, and 10) DTW. This was done for each harvest block to allow for detailed per block analyses in terms of 1) histogram and frequency distribution assessment, and 2) to determine how the clearance data were affected by pass number, machine speed, and by DTW classes.
2.7. Histogram and Frequency Distribution Assessment
The histograms for the normalized clearance data were clustered about 1, with values > 1 trailing off sharply, while values < 1 trailed off slowly towards 0, following a probability distribution function (pdf) given by (Devore, 1999) :
where , for which cmax is the maximum normalized clearance value (range of 0.2 - 0.6), Γ(k) is the gamma function, k is its shape parameter, and θ and f are scale parameters for x and c, respectively. For the pdf non-linear least-squares fitting process, cmax, k, θ and f were used as adjustable parameters, with f kept in common across all blocks. Using Equation 1 implies the following properties for : mean = k θ, variance = k , skewness = 2/k0.5, mode for k ≥ 1 = (k − 1)θ (Devore, 1999) . Hence, increasing k and θ implies widening the distribution function.
2.8. Multiple Regression Analysis
The normalized clearance data with ≤1 were summarized by way of a pivot table using the following class specifications:
1) speed (m/sec);
2) number of passes(10, 20, 30, etc.);
3) DTW (1 = 0 - 0.5 m, 2 = 0.5 - 1 m, 3 = 1 - 2 m, 4= 2 - 4 m, 5 = 4 - 8 m, 6 = 8 - 16 m, 7 = 16 - 32 m, 8 ≥ 32 m);
4) forwarding (1) versus returning (0).
The data so assembled were examined using multivariate regression analysis, using the ≤1 normalized clearance data as dependent variables, and number of passes, speed, and DTW classes as independent variables, by harvest block. Other topographic derivatives were analyzed (elevation, slope, aspect), but were found to be insignificant in comparison to the machine parameters.
3. Results and Discussion
3.1. Wood-Forwarding Track Patterns
Figures 6-8 show the wood forwarding and end-of-summer drainage (DTW)
Figure 6. Wood-forwarding track patterns within JD11 harvest blocks, showing clearance sensor points overlaid on the DEM-derived DTW patterns, based on 4 ha upslope flow-initiation areas.
patterns per harvest block, by machine type. As mapped, drainage varied in area by block from poor to imperfect (DTW < 0.5 m), and moderate (0.5 < DTW <1 m) to well drained (DTW > 1 m). Within the DTW < 0.5 m zone, soils were generally wet to moist. Within the 0.5 < DTW < 2 m zone, soil wetness tended to be transitional from wet to dry depending on extent and weather-dependent upslope water seepage. Some of the main multiple-pass TC wood-forwarding tracks straddled across these transitions, and were rutted extensively.
The average number of passes, as well as speed, and normalized clearance are listed in Table 3. Typically, machine clearance increased towards 1 with increasing pass numbers along the same track regardless of loading. However the initial unloaded passes decreased the clearances the most (Figure 9). Conversely, machine speeds were significantly faster when empty than when loaded (p = <0.001). By machine type, wood forwarding speed was lower for JD11 and JD15 than for TC, while clearances increased with increasing machine speed, and especially so for JD11 (Figure 9).
Figure 7. Wood-forwarding track patterns within JD15 harvest blocks, showing clearance sensor points overlaid on the DEM-derived DTW patterns, based on 4 ha upslope flow-initiation areas.
3.2. Normalized Clearance Distribution Patterns
The histograms of the standardized clearances per block followed a left-skewed pattern, as shown in Appendix 1 (Figure A1-1, Figure A1-2, and Figure A1-3). The corresponding Equation (1) generated probability distributions were overlaid on these histograms, with corresponding best-fitted k, cmax and θ values listed in Table 4, by harvest block. Examining these values revealed that log10k correlates with θ such that
In addition, the combination of k and θ as in log10kθ correlates with cmax and machine type, i.e.
(JD11, JD15 and TC coded 1 when present, otherwise 0).
Figure 8. Wood-forwarding track patterns within TC harvest blocks, showing clearance sensor points overlaid on the DEM-derived DTW patterns, based on 4 ha upslope flow-initiation areas.
In general, increasing kθ values reflect a broadening of the clearance frequency distributions, and Equation (3) implies that this broadening increases with increasing cmax, with further increases observed for the JD11 operations. This can be attributed to differences in machine operation, as follows:
1) JD11 was mainly used for commercial thinning, which involved generating brushmats from tree delimbing and topping by single-grip harvesters along
Figure 9. Mean machine specific normalized clearance versus number of passes (first 10 and all passes), empty versus full loads, and machine speed versus empty or full loads.
Table 3. Variability in normalized clearance, number of passes, speed by machine and operation type.
Table 4. Best-fitted Equation (1) and Equation (5) regression results for the clearance frequency distributions and normalized clearances ≤1, by machine and block; speed, passes, and DTW coefficients (×10−2); f = 0.122 ± 0.002.
1Blocks 1, 5, 6, 22, and 27 had missing sensor data and were therefore removed from the analysis.
wood forwarding trails; the broadening of the clearance distributions would be due to the ultrasonic signals bouncing off machine-induced brush-mat sagging and lifting.
2) By machine type, the kθ product varied as follows: JD11 = 3.9 ± 0.4, JD15 = 2.8 ± 0.1; TC = 2.5 ± 0.1; i.e., JD11 significantly higher than JD15 and TC (p-value < 0.001), and JD15 is significantly higher than TC (p-value = 0.14).
3) Using the <10% tail of the standardized clearance marker produced the following sequence: JD11 = 0.67 ± 0.02; JD15 = 0.80 ± 0.02, TC = 0.834 ± 0.01 (p-value < 0.05). Hence, forwarding wood on brushmats produced longer standardized clearance trails towards zero-clearance than forwarding on bare ground. By Equation 4, the 10th percentile clearances per block are directly relatable to the best-fitted Gamma distribution parameters as follows:
4) Trails for the shelterwood and clear-cut operations without brushmats had standardized clearance peaks at or near 1. Clearances greater than 1 were due to machine movements over stumps, rocks and uneven ground. The cmax ranges were similar by machine type, as follows: 1.04 < JD11 < 1.22; 1.14 < JD15 < 1.30; 1.13 < TC < 1.29.
3.3. Box Plots
The normalized data with clearances < 1 are presented in Appendix 2 by way of box plots showing 10th, 25th, 50th, 75th and 90thnormalized clearance percentiles and associated outliers per block and machine type by number of passes, speed and DTW classes (Figure A2-1, Figure A2-2, and Figure A2-3). These plots show that machine clearances varied from block to block in relation by machine operations. For TC, these variations were in part attributable to the variable June to November cumulative precipitation pattern from 2012 to 2014 (Figure 10), with lower clearances less prominent in the fall of 2012 following dry summer conditions, but more prominent in summer blocks where DTW < 1 m (Blocks 32, 35, 36; Figure A2-3).
Figure 10. Weekly and cumulative May to December rainfall for 2012-2014 (top), also showing 10th, 25th, 50th, 75th and 90th percentile box plots for the normalized TC sensor clearances ≤ 1.
3.4. Normalized Clearance ≤ 1
Regressing the normalized clearances ≤ 1 values against number of passes, speed and DTW classes by way of
produced the best-fitted interception and regression coefficient results also listed in Table 4, with R2 values ranging from near 0 to 0.67. These results are summarized as follows:
1) Normalized clearances increased with increasing machine speed, with only one reduction registered on a wet TC trail (block # 32).
2) Increasing the number of wood-forwarding passes per track had mostly positive to no clearance effects, with blocks # 18 and #25 being exceptions. For the TC wood-skidding operations, the influence of pass numbers per track on clearances varied from mainly positive during dry conditions to negative during wet conditions.
3) The influence of DTW on TC clearances also varied from positive to negative. The positive trend occurred when the ground was dry due to increased soil resistance to compaction. The negative trend occurred when operating on wet ground during and after soil-saturating rain events, due to decreased soil resistance to compaction and tire slippage.
Analyzing the normalized clearance intercepts in Table 4 in terms of their associated speed, number of passes and DTW regression coefficients generated the following multiple regression result:
This equation implies that the normalized clearance intercept:
1) is equal to 0.910 ± 0.04, or 75% to 78% actual machine-to-ground clearances on average, when machine speed, number of passes and DTW have no clearance influence; in terms of actual depth, this number amounts to about 12 cm; a considerable part of this would be due to the compression of the organic forest floor accumulations;
2) decreases below 0.910 on dry to moist and wet soils with increasing number of passes and/or machine speed, this would be due to repeating soil compaction along the same track and increased shearing stress with increased tire rotation;
3) decreases below 0.910 as the DTW-projected influence on soil resistance to compaction becomes stronger from very poor (where DTW is near 0) to well and excessive (where DTW is >1 m); i.e., a positive DTW influence causes the regression intercept to be lowest at DTW = 0; a negative DTW influence would do the opposite and may occur where DTW > 0 due to increased slope-induced wheel slippage.
Equation (6) suggests that the lowest normalized intercept for class-averaged normal clearances ≤ 1 would amount to 0.910 − 1.23 × (kpass = 0.19) − 0.041 × (kspeed = 2.34) − 0.052 × (kDTW = 3.65) = 0.39. The lowest average normalized clearance intercept per number of pass, speed and DTW classes registered in Table 4 was 0.60, for block # 30 (JD15, Nov. 2012). This is approximately equal to chassis clearance, and this occurred on a well-drained soil (DTW > 2 m) following a rain event. The lowest non-averaged normalized signal-to-ground surface distance was 0.43, i.e., equivalent 60.1 cm, thereby amounting to 19.5 cm below the reflecting ground surface. This occurred on block # 15 in January 2013 when the ground was covered with snow.
The correlation matrix in Table 5 pertaining to the four Equation (6) variables indicates that the influence of number of passes and machine speed on the clearance intercepts were positively correlated with one another, while both were negatively correlated to DTW clearance influence. Hence, DTW-influenced softening of the ground tends to decrease machine speed and reduces repeat traffic along the same track.
3.5. Block-Specific Examples
An example of low to extensive rutting is shown in Figure 11 for TC block # 41 (Sept.-Oct. 2014), where soil conditions varied from dry to moist to wet due to season- and weather-induced variations in local DTW patterns. Moderate soil compression occurred along trails where soils remained dry at DTW > 1 m, with compression decreasing towards zero with increasing number of passes. On the landing site where DTW was <0.5 m, soil ruts were deeper and number of passes per track increased up to 270. For the trail through transitional and somewhat sloped DTW zones, normalized clearances decreased with increasing number of passes per track. In general, the extent of soil compaction was highest at or near the plastic soil moisture limits, while soil displacement was highest at and above the liquid soil moisture limits.
The data for machine clearance, number of passes per track and speed are illustrated in Figure 12 for harvest blocks 10, 32, 33, 43. The mapped dots so shown demonstrate the following effects:
1) Normalised clearances along the central wood forwarding road in block 10 approached 1 at high track numbers and speed.
2) Rutting as indicated by persistently low clearances occurred along the most
Table 5. Correlation matrix (with p-values) for Equation (6) variables and associated block entries in Table 4.
n = 49.
Figure 11. (a): TC block #41 showing weather-affected DEM-derived DTW patterns based on assuming 4 ha (top-left), 1 ha (top-right) and 0.25 ha (bottom-left) of upslope flow accumulation area for channel-flow initiation. The high-resolution surface image (bottom-right) shows ruts with low machine clearances along high multiple-pass connector trails; (b): block 41 normalized clearance versus increasing number of passes along dry and wet forwarding trails and on landing site; (c): dry trail (no rutting), wet trail (deep rutting), and landing site (braided rutting) images. Normalized clearance in terms of chassis- to sensor-to solid-ground distances (front): M0/S0 = 0.63.
frequently used TC skidding passes as shown for blocks 43, 32, and 33; elsewhere, TC clearances tend to be near 1. In contrast, the JD11 clearances along brush-matted wood forwarding tracks were lower due to brush-mat flexing.
3) TC machine track speeds varied more than JD11 track speeds, with TC being faster on higher DTW ground. JD11 speeds were consistently slower along the commercial wood forwarding tracks.
(a) (b) (c) (d)
Figure 12. Visual correlation between DTW and normalized clearance (a), number of passes per track (b); and machine speed (c); for blocks 47, 32, 34, and 10. Ortho imagery of blocks (d) showcasing harvest type and soil disturbance (Imagery Source: Esri, DigitalGlobe, GeoEye, Earthstar Geographics, CNES/Airbus DS, USDA, USGS, AEX, Getmapping, Aerogrid, IGN, IGP, swisstopo, and the GIS User Community).
3.6. Further Observations
3.6.1. Clearances by Machine Type and Number of Passes
Some of the normalized clearance differences would be caused by differences in machine weight, loads, tire pressure, number of tires, and whether the tires were chained or tracked (Table 1). For example, TC had a higher front footprint pressure about 2 to 3 times higher than the JD11 and JD15 forwarders. Hence―according to the Table 1 entries―soil compression should have been deeper with TC than with JD11 and JD15 operations. This was indeed the case along repeated TC forwarding tracks on moist to wet ground (Figure 12). Under dry conditions, however, normalized clearance patterns peaked around one, except for the lower JD11 brush-mat clearances. Compared to JD15, JD11-exerted footprint pressure was about double (Table 1). Consequently, and consistent with the experimental machine-load and brush-mat observations by Labelle and Jaeger (2012) , JD11brush-mat clearances were not only lower than for JD15 brush-mat clearances, but JD11 clearance were also significantly lower during forwarding than returning (Figure 9). In terms of increasing number of passes, brushmat clearances flattened through repeated crushing and consolidation.
3.6.2. Machine Speed
Machine speeds are affected by soil-tire interactions, including soil compaction and soil displacement (Shmulevich et al., 1998) . In this regard, Liu et al. (2009) reported increasing soil displacement with increased tire rotation on dry ground. On soft ground, however, slow traffic increased soil compression and soil displacement ( Grahn , 1991
3.6.3. Bare-Ground Operations
When operating on bare ground during clear-cutting and shelter-wood harvesting, off-road clearances can also be expected to decrease with increasing footprint pressure and increasing number of passes due to increasing soil compression ( Jamshidi et al., 2008
3.6.4. Season and Weather Details
With varying soil wetness by weather and across seasons, there is a general correspondence between TC clearances and ground conditions (Figure 10). For TC, lower clearances and deeper soil rutting were incurred within or across the harvest blocks on account of: 1) low DTW locations, 2) snow accumulations, 3) snow melt and prolonged rain events, and 4) low evapotranspiration ( Raven et al., 1999 ; Ács et al., 2011 ; Jones & Arp, 2017) before leaf-out (May-June) and after leaf-fall (October-November). In contrast, JD11 and JD15 clearances were less relatable to the varying ground conditions due to low footprint pressure (Table 1) and brush-mat operations.
An example of weather-induced soil rutting on dry ground (4 < DTW < 32 m) is shown in the form of normalized clearances < 0.8 in Figure 13. This occurred during TC operations in block #48 two days after a 110 mm rain event in July 2014 (Hurricane Arthur). Similarly, Block #41 operations immediately led to extensive soil rutting after a 30 mm rain event on wet October soils led to extensive soil rutting (Figure 11). To avoid such occurrences, wood forwarding operations were generally deferred to occur on dry ground during spring, summer and fall,
Figure 13. Normalized clearances across block 2, showing low TC clearances on DTW > 4 m ground following a 110 mm per day storm event in July 2014, with two close-ups (bottom panels; A, B red boxes in top panel).
and on frozen ground during winter. During the JD11 and JD15 winter operations (block #1, 12 - 15, 17 - 19, and 28 - 30), the normalized clearances were mostly due to ultrasonic reflections from snow surfaces. These clearances decreased with increasing number of passes due to snow compaction and track widening. This effect was most notable at NWU where the snowpack was almost twice as deep than at the southern locations (Figure A3-1). While snow compaction has no effect on frozen soil rutting, it has been shown to cause soils along the tracks to freeze deeper and longer, thereby delaying soil thawing along the established tracks ( Grady, 1982 ; Garcia et al., 2015) .
4. Conclusions and Concluding Remarks
The assessment of the machine clearance data provided the following insights:
1) The normalized clearances were affected by brushmat versus bareground operations, with the former producing lower but still broader clearance distributions than the latter.
2) All of the block-based clearance histograms followed asymmetric Gamma frequency distributions.
3) Clearances increased systematically with increasing passes towards 1 on dry ground due to successive soil compression and track widening. On wet ground, clearances would decrease because of successive soil re-displacements.
4) The TC centralizing wood-forwarding pattern following clear-cutting and shelter-wood operations included long tracks with number of passes exceeding 100 per track, which would lead to rutting along wet and wet-to-dry transitional ground conditions.
5) Number of passes, machine speed and low to high DTW classes all affected the normalized clearance results in terms of block and weather-specific conditions, with least effects registered along brushmat and dry-ground tracks, all based on fairly low machine footprint pressures in the order JD15 < JD11 < TC.
6) Clearances along the same bare-ground tracks do not necessarily reflect rut depths due to track broadening and using machine-fixed positions for machine-to-ground distance monitoring. On brushmats, the clearance data reflect brushmat re-conditioning due to repeated track traffic.
7) Actual machine clearance can be obstructed due to ultrasonic sensor blockage by way of debris, and by the presence of snowpack accumulations.
8) The extent of soil rutting versus machine clearance needs to be evaluated separately by way of, e.g., high-resolution block surveys following completion of block operations (Salmivaara et al., 2018) , or through machine-mounted LiDAR-based ground scanning (Giannetti et al., 2017) .
9) GPS tracking of tire rotation in connection with machine speed would assist in determining actual to potential soil displacement in terms of tire slippage and tread design.
This research was supported by J. D. Irving Limited (JDI) as a part of the NSERC-sponsored Cooperative Research Development (CRD) Project on Forest Soil Trafficability. Special thanks go Greg Adams at JDI for facilitating this project, and to FP Innovations personnel for the development and installations of the ultrasonic sensor data loggers. Also many thanks to the Shane Furze, Doug Hiltz, and John Paul Arp for help with data retrieval and GIS processing.
Normalized clearance histograms perharvest block by machine type, with best-fitted Gamma distribution functions overlaid.
Figure A1-1. Histograms for the JD11 normalized clearance data, by block. Normalized clearance in terms of chassis- to sensor-to solid-ground distances: M0/S0 = 0.57.
Figure A1-2. Histograms for the JD15 normalized clearance data, by block. Normalized clearance in terms of chassis- to sensor-to solid-ground distances: M0/S0 = 0.57.
Figure A1-3. Histograms for the TC normalized clearance data, by block. Normalized clearance in terms of chassis- to sensor- to solid-ground distances (front): M0/S0 = 0.63.
10th, 25th, 50th, 75th and 90th percentiles for the ≤ 1 normalized clearance data per block and machinetype by pass number, speed and DTW classes.
Figure A2-1. Normalized JD11 clearance ≤1box plots by pass number, speed class, and DTW classes per block. Blocks #1, #5 and #6 generated no clearance data due to sensor malfunction. Normalized clearance in terms of chassis- to sensor-to solid-ground distances: M0/S0 = 0.57.
Figure A2-2. Normalized JD15 clearance ≤1 box plots by pass number, speed class, and DTW classes per block. Blocks #32 generated no clearance data due to sensor malfunction. Normalized clearance in terms of chassis- to sensor-to solid-ground distances: M0/S0 = 0.55.
Figure A2-3. Normalized TC clearance ≤ 1 percentile box plots by pass number, speed class, and DTW classes per block. Normalized clearance in terms of chassis- to sensor-to solid-ground distances: M0/S0 = 0.63.
Year-round hydrological weather and ground conditions 2011-2014.
Figure A3-1. Daily variations in air temperature and precipitation, with modelled stream discharge, volumetric soil moisture content including field capacities for the wider NWU (left) and MWU + LL (right) areas, and snowpack depth, frost depth, and volumetric soil moisture content for the wider NWU (left) and MWU + LL (right) areas.
 Balland, V., Pollacco, J. A. P., & Arp, P. A. (2008). Modeling Soil Hydraulic Properties for a Wide Range of Soil Conditions. Ecological Modelling, 219, 300-316.
 Berisso, F. E., Schjønning, P., Keller, T., Lamandé, M., Etana, A., de Jonge, L. W., Iversen, B. V., Arvidsson, J., & Forkman, J. (2012). Persistent Effects of Subsoil Compaction on Pore size Distribution and Gas Transport in a Loamy Soil. Soil & Tillage Research, 122, 42-51. https://doi.org/10.1016/j.still.2012.02.005
 Ezzati, S., Najafi, A., Rab, M. A., & Zenner, E. K. (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, 521-538. https://doi.org/10.14214/sf.908
 Garcia, M., Özdogan, M., & Townsend, P. A. (2015). Impacts of Forest Harvest on Cold Season Land Surface Conditions and Land-Atmosphere Interactions in Northern Great Lakes States. Journal of Advances in Modeling Earth Systems, 6, 923-937.
 Giannetti, F., Chirici, G., Travaglini, D., Bottalico, F., Marchi, E., & Cambi, M. (2017). Assessment of Soil Disturbance Caused by Forest Operations by Means of Portable Laser Scanner and Soil Physical Parameters. Soil Science Society of America Journal, 81, 1-23.
 Grahn, M. (1991). Prediction of Sinkage and Rolling Resistance for Off-the-Road Vehicles Considering Penetration Velocity. Journal of Terramechanics, 28, 339-347.
 Horn, R., DomżaÖ, H., SÖowińska-Jurkiewicz, A., & van Ouwerkerk, C. (1995). Soil Compaction Processes and Their Effects on the Structure of Arable Soils and the Environment. Soil & Tillage Research, 35, 23-36.
 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.
 Jones, M.-F., & Arp, P. A. (2017). Relating the Cone Penetration and Rutting Resistance of Soils to Variations in Soil Properties and Daily Moisture Variations over Time. Open Journal of Soil Science, 7, 149-171. https://doi.org/10.4236/ojss.2017.77012
 Labelle, E. R., & Jaeger, D. (2011). Soil Compaction Caused by Cut-to-Length Forest Operations and Possible Short-Term Natural Rehabilitation of Soil Density. Soil Science Society of America Journal, 75, 2314-2329. https://doi.org/10.2136/sssaj2011.0109
 Liu, K., Ayers, P. D., Howard, H., & Anderson, A. (2009). Influence of Turning Radius on Wheeled Military Vehicle Induced Rut Formation. Journal of Terramechanics, 46, 49-55. https://doi.org/10.1016/j.jterra.2009.02.004
 Murphy, P. N. C., Ogilvie, J., & Arp, P. A. (2009). Topographic Modelling of Soil Moisture Conditions: A Comparison and Verification of Two Models. European Journal of Soil Science, 60, 94-109. https://doi.org/10.1111/j.1365-2389.2008.01094.x
 Murphy, P. N. C., Ogilvie, J., Meng, F.-R., White, B., Bhatti, J. S., & Arp, P. A. (2011). Modelling and Mapping Topographic Variations in Forest Soils at High Resolution: A Case Study. Ecological Modelling, 222, 2314-2332.
 Salmivaara, A., Miettinen, M., Finér, L., Launiainen, S., Korpunen, H., Tuominen, S., Heikkonen, J., Nevalainen, P., Sirén, M., Ala-Ilomäki, J., & Uusitalo, J. (2018). Wheel Rut Measurements by Forest Machine-Mounted LiDAR Sensors—Accuracy and Potential for Operational Applications? International Journal of Forest Engineering, 29, 1-12.
 Taghavifar, H., & Mardani, A. (2014). Effect of Velocity, Wheel load and Multipass on Soil Compaction. Journal of the Saudi Society of Agricultural Sciences, 13, 57-66.
 Vega-Nieva, D. J. D., Murphy, P. N. C., Castonguay, M., Ogilvie, J., & Arp, P. A. (2009). A Modular Terrain Model for Daily Variations in Machine-Specific Forest Soil Trafficability. Canadian Journal of Soil Science, 89, 93-109. https://doi.org/10.4141/CJSS06033
 White, B., Ogilvie, J., Campbell, D. M. H., Hiltz, D., Gauthier, B., Chisholm, H. K. H., Wen, H. K., Murphy, P. N. C., & Arp, P. A. (2012). Using the Cartographic Depth-to-Water Index to Locate Small Streams and Associated Wet Areas across Landscapes. Canadian Water Resources Journal, 37, 333-347. https://doi.org/10.4296/cwrj2011-909