vol. 6, pp. 255-261
Copyright © 2013 by the Italian Society of Silviculture and Forest Ecology
doi: 10.3832/ifor0773-006

Research Articles

S. Grigolato, M. Pellegrini, R. Cavalli

# Introduction

Forest road network investigations provide a primary method for evaluating forest accessibility in terms of on-road transportation and forest operations. Forest roads play a dual role in providing both access to the forest and transportation mobility and, together with extraction routes, normally form a relatively even network in the forest. A forest road network is therefore the result of the combination between the main forest roads and the timber extraction routes ([16]).

Forest road network, which includes both forest and secondary public roads, has received much attention in recent years due to its function and effect on forested rural landscapes and the related environment ([22], [15]) and its primary role in the efficiency of logging operations and forest transportation ([16]). As a consequence, the importance of the multi-functional role of the forest road network has been widely recognized, being the interface between forested lands and the roads. It provides access not only for timber extraction and management but also for water protection and wildlife habitat improvement, fire control and many recreational activities ([1]).

On the other hand, the presence of roads has many undesirable effects on both biotic and abiotic components of the forest. Several reviews have been published in the literature about the negative effects of forest roads ([8], [11]). As an example, there has been considerable interest in the analysis of the forest road network as a source of accelerated soil erosion and the production of sediment ([5], [14]), also using GIS models ([34]).

The maintenance and management of forest road networks are also an essential element to be investigated ([26]), and have recently been analyzed by Coulter et al. ([9]) with the intent of protecting road infrastructure and limiting negative environmental effects by considering non-economic benefits. Furthermore, an assessment of the impact of vehicle traffic moving through a designated area via a road network involves specific problems of forest management.

Traffic increases sediment generation rates through detachment, abrasion and crushing of the road surface material ([33]). The assessment and subsequent prevention of soil erosion from forest roads requires an understanding of how road design and maintenance affect sediment production ([20]). Empirical evidence providing an insight into how sediment yield is affected by traffic is reported by Bilby et al. ([5]) who highlighted that the amount of sediment produced is also related to traffic rate. Researches by Reid ([30]), Reid & Dunne ([31]), Grayson et al. ([13]), Luce & Black ([21]), Ramos-Scharron & MacDonald ([28]), Sheridan & Noske ([32]), Ziegler et al. ([36]) were specifically aimed at determining the effects of traffic on road erosion, confirming that increasing erosion rates are related to increased traffic.

Despite the fact that the production of sediment and the maintenance of the forest road network are strictly connected to the traffic loads, few papers consider the analysis and quantification of the traffic load generated by forest operations ([27]). Several examples are reported in the literature on traffic load analysis based on the planned wood extraction volume ([19]) or on the potential yield located within a buffer area around the forest road segment ([24]).

In this study a 21-year period of logging operations are investigated with the aim of verifying the spatial and temporal distribution of the traffic load generated by wood transportation in a mountainous area of the Alps. The first goal was to examine and describe the logging operations in terms of extraction distance and terrain gradient. A novel GIS approach has therefore been developed and applied to characterize the spatial and temporal distribution of the traffic load. The distribution of the traffic was then analyzed according to the characteristics of the forest road network.

The investigation also aims to compare the traffic loads generated by the effective logging operations (actual harvest) vs. the potential traffic load on the forest road network by considering the average wood mass extracted yearly as calculated from the prescribed yield (maximum allowable harvest).

# Materials and methods

## Study area

The study area is located in Chiese valley in north-eastern Italy. The managed forest area covers 25 658 ha, partitioned into 27 Forest Management Plans (FMP). The productive forest area covers 14 905 ha with a prescribed annual yield of 17 000 m3. The productive forest area is managed using a close-to-nature silvicultural approach, ranging from single-tree selection cut to irregular group cut. Cutting operations include the harvest of large, economically mature trees, the tending of the intermediate size classes and the removal of trees without prospect of further growth. In the largest groups of adult trees, plant selection follows the shelterwood method with a long regeneration period.

According to the official forest road register and the indications of local stakeholders, the road network within the study area has not been substantially modified in the last 25 years in terms of expansion but merely in terms of construction standards.

## Data collection of the forest road network

The existing public (primary and secondary) roads and forest roads were extracted from the respective road network databases, merged together and then verified and updated by GPS surveys. Landing and vehicle constraints, such as road gradient and width were also surveyed as points of interest using GPS. As for the forest road network, the vehicle constraints were categorized according to the most common Italian classification currently in use ([17]) as reported in Cavalli & Grigolato ([7]). As indicated in the Provincial Law LP 11/07 ([2]), the forest roads were identified according to their use: “exclusive use for forest activities” and “not exclusive use for forest activities”. The remaining roads with open access and used for forest operations and forest area access were simply classified as secondary public roads according to the National Road Code.

## Data collection and spatial localization of the logging operations

Data on logging operations were collected from the local Forest Service offices. For each logging operation between 1987 and 2007, information was collected on cutting and merchantable wood volume, harvested area and the identification code of the stands. The collected data included 2231 harvesting operations.

The 27 FMP were unified in a single geographic FMP database, ad the logging operations grouped into a single geographic database by the identification codes of each stand and harvesting site. This resulted in all 2231 logging operations being spatially located.

In order to characterize the logging operations by location and topography, the skidding distance and the terrain gradient of the harvesting area were determined by GIS analysis. The analysis was based on a high resolution (up to 1.0 m) Digital Elevation Model (DEM) with an altitude accuracy ranging from 30 to 60 cm. The DEM was derived from LiDAR (Light Detection and Ranging) aerial survey data acquired at an altitude of 1800 m a.s.l. in 2006 and 2007 by an OPTECH ALTM 3100 sensor mounted on an aircraft with a point density of 0.5 per square meter ([3]).

According to Nakazawa et al. ([23]), Ranta ([29]) and Pentek et al. ([25]), the center of each harvesting site was supposed to correspond to the center of mass of the polygons identifying the harvesting area. A first geometrical extraction distance was calculated by computing the distance from the center of mass to the closest point of the forest road network. The x and y coordinates of the closest point of the forest road network were also recorded in order to identify the most likely landing sites ([18]). However it has to be pointed out that the method only estimates the positions of the landing sites, therefore there could be slight differences between the estimated and reallocation of landings.

High resolution DEM was used to extract the difference in elevation between the hypothetical landing sites and the centers of mass of the polygons. The mean extraction distance was then recalculated by means of the Pythagorean theorem into a more accurate distance that considered the change in elevation.

The mean terrain gradient of each harvesting site was calculated by zonal statistical analysis on the terrain gradient value. The percentage terrain gradient was previously calculated by a spatial analysis of the elevation change of the surface in the horizontal (dz/dx) and vertical (dz/dy) directions based on Horn’s method ([6]) with a resolution of 1 m of the DEM. The final database containing logging operations for each stand was thus spatially located.

## Determination of the traffic load

The traffic analysis consisted of a GIS analysis using raster data. First, the entire road network (including the forest and secondary public roads) was divided into 72 sub-networks by identifying their end points at the junction with the nearest primary road. The traffic load analysis concerned the traffic generated by trucks without trailer as this is the most used transport system in the area and because the use of tractor and trailer involves a minimal part of the road network with minimum construction standards and an episodic use. Maximum gradient was assumed as 18%, carriageway width 3 m and curve minimum radius 8 m. The field surveys identified 33 critical key constraints, mostly steep gradients, limiting the transit of trucks.

The 2231 centers of mass of the polygons identifying each harvesting area were then snapped to the closest forest road or secondary public road (sometimes used as support for the forest operations). The snapped positions were assumed to be the most likely landing sites for the logging operations ([18]). The snapped centers of mass were subsequently converted into raster data with a resolution of 1 m. A set of 21 rasters was thus generated from the logging operation database. The 21 rasters represented the extracted mass at landings for each of the 21 analyzed years. The extracted mass equivalent of the merchantable wood volume was calculated by considering green wood with a density of 860 kg m-3 for Norway spruce, 920 kg m-3 for fir, 900 kg m-3 for larch, 880 kg m-3 for Scots pine and 1005 kg m-3 for broadleaves ([12]).

The forest road network was converted into raster data and re-classified according to the minimum accumulated route distance from the end points to the farthest points, also considering any interruptions caused by transportation barriers. A flow accumulation analysis was thus computed by calculating the accumulated flow of mass from the furthest point to the end point as the accumulated weight of all cells flowing into each down slope cell. The output raster consequently identified the wood mass accumulated along the forest road network according to the direction previously identified by the route distance analysis. This methodology was adapted from the well-known hydrological method used for the calculation of water flow accumulation to categorize the stream network in a watershed ([35]). The determination of the flow accumulation of the extracted wood mass was therefore based on the direction that was identified by the distance instead of the down slope as used by the hydrological method. The calculation was repeated for each of the 21 years in order to characterize each road segment by the flowed wood mass.

The number of truck round trips for the transportation of the extracted mass from each logging operation was then calculated according to the following formula proposed by Nevečerel et al. ([24] - eqn. 1):

$$n = \frac{M_{land}} {m_{max} - m_{empty}}$$

where n is number of truck round trips, mland is the extracted wood mass (t) at landing, mmax is the gross vehicle weight (t), mempty is the mass (t) of the empty truck.

The number of truck round trips was first calculated for each year and then as the sum of the 21 years. The mass of empty trucks and the maximum mass allowed was considered according to the results of interviews with nine logging companies, which regularly operate within the study area; it was found that the most common truck for roundwood transportation is a single three-axle truck with double wheels at the rear end axles, with a gross vehicle weight of 26 t and a maximum payload of 13.5 t. The traffic load of each cell of the forest road network raster was therefore calculated considering the mass of the truck according to (eqn. 2):

$$Tl=n \cdot (m_{max}+m_{empty})$$

where Tl is the traffic load, n is the number of truck round trips, mmax is the gross vehicle weight (t), mempty is the mass (t) of the empty truck.

The analyses on the spatial localization of logging operations and traffic load were conducted by ASDGIS Desktop 9.3 with ASDIOGP extension (ESRI Inc., Redland CA) for GIS analysis.

## Statistical analysis

As for logging operations, the statistical analysis was aimed to verify if the means were:

• equal for different extraction distances (m) by considering 6 intervals of distance: <100, 101-200, 201-300, 301-400, 401-500, >500;
• equal at different terrain gradient (%) by considering 5 intervals of the gradient: <20, 21-40, 41-60, 61-80, >80.

For traffic load, the statistical analysis aimed to verify if the means were:

• equal in each road network segment during the 21 analyzed years;
• equal in relation to the use of the roads: “exclusive use for forest activities”, “not exclusive use for forest activities” and “secondary public roads”;
• equal in relation to two types of road surface: gravel and asphalt.

Furthermore, in order to evaluate the average yearly effective traffic loads generated by forest operations in the 21-year period vs. the yearly potential traffic load, the traffic load analysis also considered the average yearly extracted wood mass as calculated from the expected yield. The traffic load generated by the recorded logging operations of each forest road segment was therefore compared to the potential traffic load that would be generated by the logging operations for the theoretical prescribed yield as reported by the FMP. Hence, the traffic load procedure was applied also considering the theoretical extracted wood mass to the landing.

For each of the above hypothesis, goodness-of-fit to normal distributions was verified by the Kolmogorov-Smirnov test, and the equality of variances among different groups by the Levene test. One-Way ANOVA was applied to test the hypothesis that the groups described above were not significantly different. In order to determine which groups differ from each another, the post-hoc test considered the unbalanced samples’ size. Tukey’s honestly significant difference test (HSD) for unbalanced sample size was thus applied.

Statistical analyses were performed using the SPSS 18.0 statistical package (SPSS Inc., Chicago, IL) with α = 0.05.

# Results

## Logging operations between 1987 and 2007

The collected information covered 2231 harvesting sites in 21 years for a total cutting volume of 444 381 m3, corresponding to 406 531 m3 of merchantable wood (with bark). Norway spruce represented 61.9% of the merchantable wood volume, 14.3% was fir, 12.7% larch, 4.5% Scots pine and 6.6% broadleaves (mainly beech).

Over the 21 years considered, the average extracted volume per logging operation was 182.27 ± 105.45 m3 (SD). A remarkable portion of logging operations (69.2%) had an extracted volume of up to 300 m3 (Fig. 1). There is also evidence suggesting that a relevant part (57.3%) of the wood volume extracted every year was obtained in logging operations with a merchantable wood volume higher than 300 m3 (Fig. 2).

Fig. 1 - Yearly distribution of logging operations according to the extracted wood volume.
Fig. 2 - Yearly distribution of total extracted volume according to the extracted wood volume per site.

Regarding the extraction distance, the ratio between the average geometrical distance and the average slope extraction distance was 1.61.

Average slope extraction distance per logging operation was 154.15 ± 112.76 m2 over the 21 years considered. As depicted in Fig. 3, the majority of logging operations (89%) over the 21 years had an extraction distance shorter than 300 m. As a consequence, the percentage of extracted merchantable wood volume within 300 m (Fig. 4) was considerable (72% of the total extracted volume).

Fig. 3 - Yearly distribution of logging operations according to the extraction distance.
Fig. 4 - Yearly distribution of extracted volume according to the extraction distance.

Unclear results on the homogeneity of variances among the six extraction distance classes considered were obtained by applying the Levene test. The application of the more robust Brown-Forsythe test confirmed the homoscedasticity in the extracted wood volume among the above groups (p-value = 0.091). One-way ANOVA did not revealed significant differences in the extracted wood volumes among the six groups considered (p-value = 0.610 - Tab. 1).

Tab. 1 - Descriptive statistics of the extracted wood volume at different extraction distance intervals. Different letters in the means values indicate significant differences at the 5% level.

In the time span considered, logging operations were carried out often on steep slopes (average: 62.03 ± 12.51 %). A large percentage of the 2231 logging operations were located on terrain with a gradient higher than 40% (Fig. 5). The merchantable wood extracted from slopes higher than 40% was about 95% of the total extracted volume (Fig. 6).

Fig. 5 - Yearly distribution of logging operations according to the terrain gradient.
Fig. 6 - Yearly distribution of extracted wood volume according to the terrain gradient.

As for the merchantable wood extracted in the mentioned slope classes, no departures from normal distributions within groups (p-value=0.063) or from equality of variances among groups (p-value = 0.085) were found. Significant differences among slope class means were found by applying One-Way ANOVA (p-value = 0.024). The HSD test was thus applied in order to compare class means. The results demonstrated that extracted wood volumes on slopes up to 20% and beyond 80% were significantly lower than those on slopes within the range 20-80% (Tab. 2).

Tab. 2 - Descriptive statistics of the extracted wood volume at terrain gradient intervals. Different letters in the means values indicate significant differences at the 5% level.

## Traffic load between 1987 and 2007

Forest roads currently represent 53.5% (364.8 km) of the total road network (681.9 km). Some 38.9% of the forest road network is classified as “exclusive use for forest activities” and most of that (84.5%) is suitable for the transit of trucks without trailer.

Fig. 7 - Yearly distribution of traffic loads along the forest road network.

As for the means of each traffic load in the different road network segment, no significant differences from normal distribution were found (p-value = 0.072). As the null hypothesis of equal group variances was rejected (p-value = 0.014), the variance of the median was verified by the non-parametric Kruskal-Wallis test, which also rejected the null hypothesis that the medians were equal. This supports the hypothesis that traffic loads generated by forest operations could take different yearly intensities and distributions (data not shown), probably because of the different locations and intensities of logging operations.

Considering the total sum of traffic loads for each road network segment over the 21 years analyzed, a comparison of traffic loads was made according to the use (exclusive forest use, not exclusive forest use and secondary public road), surface (asphalt and gravel) and gradient (<6%, 6-12% and >12%) of the road segment (Tab. 3).

Tab. 3 - Means of the total traffic loads per road segment trafficked in the period 1987-2007 in relation to the use, surface and gradient of the road. Different letters in the means values indicate significant differences at the 5% level.

As for the means of traffic loads in relation to road use, no significant departures from normal distribution (p-value = 0.152) and from homoscedasticity (p-value = 0.062) were found. Significant differences among groups were found after One-Way ANOVA (p-value = 0.008). The HSD test revealed that the traffic load on secondary public roads was 260.1% higher than on the forest roads used exclusively for forest operations, and 37.4% higher than on forest roads not exclusively used for forest operations (Tab. 3).

With regard to the traffic load on asphalt and gravel roads, equal variances among groups were confirmed by Levene’s test (p-value = 0.102). The One-Way ANOVA rejected the null hypothesis that the means were equal (p-value = 0.000). The traffic load along asphalt roads was 133.8% higher than on gravel roads.

As for the traffic load on roads grouped according to their average slope (three classes: <6%, 6-12% and >12%), homoscedasticity was confirmed by Levene’s test (p-value = 0.851). In this case, the application of one-way ANOVA confirmed the null hypothesis that means were equal (p-value = 0.897 - Tab. 3).

The effective traffic loads recorded over the period 1987-2007 due by logging operations (and corresponding to the effective harvested yield) was thus compared with the potential traffic loads that would be generated in the same period by harvesting the prescribed yield indicated by the FMP. As depicted in Fig. 8, the effective traffic load along the road network was mostly lower than the potential traffic load, except in 36 cases.

Fig. 8 - Effective traffic loads (actual harvest) vs. potential traffic loads (maximum allowable harvest) by the prescribed yield in the period 1987-2007.

# Discussion

In this study a period of 21 years of logging operations was investigated with the aim of verifying the spatial and temporal distribution of the traffic load generated by wood transportation in a Alpine area in Italy.

The study highlighted that at least 60% of the logging operations had an extracted volume of up to 300 m3. The extraction distance of most of the logging operations was within 500 m from the road network. This situation is consistent with the results reported by Baldini et al. ([4]) indicating that most of the forest area (78.6%) in central-northern Italy is located within a distance of 500 m from the road network. Moreover, most of the logging operations, and consequently most of the extracted wood volume, are located on terrains with a slope gradient higher than 40% and lower than 80%.

The logging operation results are quite consistent with a similar analysis carried out in Japan ([23]), reporting that thinning operations tend to ignore remote forest areas (more than 300 m apart from the forest road) and declivity steeper than 90%.

With regard to the traffic load analysis, previous works have documented the function and effects of the forest road network, pointing out the relationships between the traffic and sediment productions ([30], [31], [20], [21], [33], [32]). However, few studies have considered the assessment of the traffic load due by forest operations. These were examples of traffic load analysis based on either the planned wood extraction volume ([19]) or the potential yield located within a buffer area around the forest road segment ([24]).

In this study, a 21-years period of logging operations was analyzed and mapped in order to investigate the long-term effect of the traffic load. To this purpose, the GIS analysis has taken into consideration the segments of the secondary public and forest road network suitable for roundwood transportation by trucks. Our results showed that most of the road segments with at least one transported load of roundwood had traffic load intensities generally less than 500 tons per year, and only 11.7% of the analyzed road network had a traffic load higher than 3000 tons per year.

The highest traffic loads were recorded on secondary public roads closer to the main public road network, where the end point of the traffic flow analysis was set. This confirmed the results obtained by Nevečerel et al. ([24]) in Croatia, although in that area a consistent part of the road network had a traffic load higher than 30 000 tons in 10 years.

The results of the study also showed that the highest traffic loads were registered on the road segments with asphalt surface, while there was no significant difference between road segments with different slopes.

In addition, the distribution of traffic loads over time and features of the road network showed that the traffic load due by logging operations may vary in intensity and distribution over the years, probably due to the different locations and intensities of the logging operations.

This investigation also revealed that the forest road network considered has an actual traffic load lower than its potential. This is due to the difference between the effective and prescribed harvest yield, the former being only 45 % of the latter (as reported by the FMP). Underexploitation of prescribed yields is quite common in north-eastern Italy, as reported in Emer et al. ([10]). On the other hand, the effective traffic load recorded was rarely higher than the potential traffic loads. By analyzing the notes reported in the harvesting schedules describing the logging operations considered in this study, most situations concerned the harvesting of a large amount of wind- or snow-thrown timber, thus explaining the anomaly in the effective traffic load on local road networks.

# Conclusions

In summary, the GIS methodology developed in this study for the traffic load determination using time data of logging operations could be useful for further analyses to identify solutions for reducing the impact of the traffic load due by logging operations on the forest environment, such as sediment production, and to identify the management priorities for the forest road network. The temporal analysis of traffic loads on the forest road network can also provide a prediction of road maintenance and eventual upgrading costs.

# Acknowledgments

The Authors express their appreciation to Dott. Andrea Bagattini for his cooperation in the collection of data. A special thanks go to Prof. Peter Ackerman from the University of Stellenbosch (South Africa) for his helpful suggestions.

# References

(1)
Akay AE, Sessions J (2004). Roading and transport operations. In: “Encyclopedia of Forest Sciences” (Burley J, Evans J, Youngquist J eds). Elsevier Academic Press, Amsterdam, The Netherlands, pp. 259-269.
(2)
Autonomous Province of Trento (2007a). LP 11/07. Provvedimenti per il potenziamento delle aree forestali e delle loro risorse. Provincia Autonoma di Trento, Trento, Italy, pp. 1-79.
(3)
Autonomous Province of Trento (2007b). Portale cartografico trentino. Provincia Autonoma di Trento, Trento, Italy.
(4)
Baldini S, Cavalli R, Piegai F, Spinelli R, Di Fulvio F, Fabiano F, Grigolato S, Laudati G, Magagnotti N, Nati C, Picchio R (2009). Prospettive di evoluzione nel settore delle utilizzazioni forestali e dell’approvvigionamento del legname. In: Terzo congresso nazionale di “Selvicoltura per il miglioramento e la conservazione dei boschi italiani”. Italian Academy of Forest Sciences, Florence 2: 717-728. [in Italian]
(5)
Bilby RE, Sullivan K, Duncan SH (1989). The generation and fate of road-surface sediment in forested watersheds in southwestern Washington. Forest Science 35: 453-468.
(6)
Burrough PA, McDonell RA (1998). Principles of geographical information systems. Oxford University Press, New York, USA.
(7)
Cavalli R, Grigolato S (2010). Influence of characteristics and extension of a forest road network on the supply cost of forest woodchips. Journal of Forest Research 15: 202-209.
(8)
Coffin AW (2007). From roadkill to road ecology: a review of the ecological effects of roads. Journal of Transport Geography 15: 396-406.
(9)
Coulter ED, Sessions J, Wing MG (2006). Scheduling forest road maintenance using the analytic hierarchy process and heuristics. Silva Fennica 40:143-160.
(10)
Emer B, Grigolato S, Lubello D, Cavalli R (2011). Comparison of biomass feedstock supply and demand in Northeast Italy. Biomass and Bioenergy 35: 3309-3317.
(11)
Forman RTT, Alexander LE (1998). Roads and their major ecological effects. Annual Review of Ecology, Evolution and Systematics 29: 207-232.
(12)
Giordano G (1981). Il legno: caratteristiche e lavorazioni industriali. UTET, Turin, Italy. [in Italian]
(13)
Grayson RB, Haydon SR, Jayasuriya MDA, Finlayson BL (1993). Water quality in mountain ash forests - separating the impacts of roads from those of logging operations. Journal of Hydrology 150: 459-480.
(14)
Gucinski H, Furniss M, Ziemer R, Brookes M (2001). Forest Roads: a synthesis of scientific information. Report no. PNW-GTR-509, Pacific Northwest Research Station, USDA Forest Service, Portland, OR, USA.
(15)
Gumus S, Acar HH, Toksoy D (2008). Functional forest road network planning by consideration of environmental impact assessment for wood harvesting. Environmental Monitoring and Assessment 142: 109-116.
(16)
Heinimann HR (1998). A computer model to differentiate skidder and cable yarder based road network concepts on steep slopes. Journal of Forest Research 3: 1-9.
(17)
Hippoliti G (1976). Sulla determinazione delle caratteristiche della rete viabile. L’Italia Forestale e Montana 31: 241-255. [in Italian]
(18)
Kinoshita T, Inoue K, Iwao K, Kagemoto H, Yamagata Y (2009). A spatial evaluation of forest biomass usage using GIS. Applied Energy 86: 1-8.
(19)
Krč J (2006). A model for evaluating forest road load by forest operations. In: “Precision forestry in plantations, semi-natural and natural forests” (Ackerman PA, Langin, DW, Antonides MC eds). International Precision Forestry Symposium, Stellenbosch University, Stellenbosch, South Africa, pp. 437-450.
(20)
Luce CH, Black TA (1999a). Sediment production from forest roads in western Oregon. Water Resources Research 35: 2561-2570.
(21)
Luce CH, Black TA (1999b). Changes in erosion from gravel surfaced forest roads through time. In: Proceedings of the Symposium “International Mountain Logging and 10 Pacific Northwest Skyline”. International Union of Forestry Research Organizations and Oregon State University, Corvallis, OR, USA, pp. 204-218.
(22)
Lugo AE, Gucinski H (2000). Function, effects, and management of forest roads. Forest Ecology and Management 133: 249-262.
(23)
Nakazawa M, Matsumoto T, Yamada Y (2004). Analysis of current forest operations in the Okumikawa forestry area by location and topography. Journal of Forest Research 9: 187-193.
(24)
Nevečerel H, Pentek T, Pičman D, Stakić I (2007). Traffic load of forest roads as a criterion for their categorization-GIS analysis. Croatian Journal of Forest Engineering 28 (1): 27-37.
(25)
Pentek T, Dragutin P, Dvorščak P, Nevečerel H (2005). Analysis of an existing forest road network. Croatian Journal of Forest Engineering 26 (1): 39-50.
(26)
Potočnik I, Yoshioka T, Miyamoto Y, Igarashi H, Sakai H (2005a). Maintenance of forest road network by natural forest management in Tokyo University Forest in Hokkaido. Croatian Journal of Forest Engineering 26 (2):71-78.
(27)
Potočnik I, Pentek T, Pičman D (2005b). Impact of traffic characteristics on forest roads due to forest management. Croatian Journal of Forest Engineering 26 (1): 51-57.
(28)
Ramos-Scharron CE, MacDonald LH (2005). Measurement and prediction of sediment production from unpaved roads, St. John, US Virgin Islands. Earth Surface Processes and Landforms 30: 1283-1304.
(29)
Ranta T (2005). Logging residues from regeneration fellings for biofuel production - A GIS-based availability analysis in Finland. Biomass and Bioenergy 28: 171-182.
(30)
Reid LM (1981). Sediment production from gravel-surfaced roads. Clearwater basin, Washington. Publication no. FRI-UW-8108, University of Washington Fisheries Research Institute, Seattle, WA, USA.
(31)
Reid LM, Dunne T (1984). Sediment production from forest road surfaces. Water Resources Research 20: 1753-1761.
(32)
Sheridan GJ, Noske PJ (2007). A quantitative study of sediment delivery and stream pollution from different forest road types. Hydrological Processes 21 (3): 387-398.
(33)
Sheridan GJ, Noske PJ, Whipp RK, Wijesinghe N (2006). The effect of truck and road water content on sediment delivery from unpaved roads. Hydrological Processes 20: 1683-1699.
(34)
Skaugset A, Surfleet C, Meadows M, Amann J (2011). Evaluation of erosion prediction models for forest roads. Transportation Research Record 2203: 3-12.
(35)
Tarboton DG, Bras RL, Rodriguez-Iturbe I (1991). On the extraction of channel networks from digital elevation data. Hydrological Processes 5: 81-100.
(36)
Ziegler AD, Sutherland RA, Giambelluca TW (2001). Interstorm surface preparation and sediment detachment by vehicle traffic on unpaved mountain roads. Earth Surface Processes and Landforms 26: 235-250.

### Cited by

 Grigolato S, Pellegrini M, Cavalli R (2013). Temporal analysis of the traffic loads on forest road networks iForest - Biogeosciences and Forestry 6: 255-261. - doi: 10.3832/ifor0773-006 Close
 Original Size Reduce the image size Enlarge the image < > > < First Previous Next Last
Close