## Decision making in forest road planning considering both skidding and road costs: a case study in the Hyrcanian Forest in Iran

iForest - Biogeosciences and Forestry, Volume 6, Issue 2, Pages 59-64 (2013)
doi: https://doi.org/10.3832/ifor0640-006

Technical Reports

Hyrcanian forest is the only forest designated for commercial timber production in Iran. Ground-based skidding is the most common extraction system used in the forest, though large parts of the forest are still inaccessible due to the low road density (1-2 m ha-1). To facilitate timber harvesting in the forest, it has been proposed to increase the road density up to 20 m ha-1. The aim of this study was to incorporate the estimated skidding costs through a time study into an existing transportation planning tool, NETWORK 2000, to help decision making on forest road building. We applied this method to the Gorazbon district in the Hyrcanian forest, where 244 links were built for the analysis. The best solution found by NETWORK 2000 indicates that all the proposed roads should be built to minimize the total skidding and road building costs.

# Introduction

Forest road building and timber harvesting are important parts of forest management activities. The Hyrcanian forest is the only commercial forest in Iran with an area of 1.9 millions ha. Ground-based skidding system is the most common mechanized timber harvesting system that has been used in Iran. The system is the process of extracting felled trees from the cutting site to a landing or roadside where the trees are processed into logs and loaded onto trucks for further transport to the sawmill. Many studies have been carried out on productivity and cost of timber harvesting, as well as influencing factors that affect the performance of each harvesting machinery ([21], [20], [15], [2], [5], [19], [18], [30], [9], [31], [4], [25], [28]). In north of Iran, Sobhani & Stuart ([27]) evaluated the harvesting system with a Clark 66BDS wheeled skidder, and several researchers have studied the productivity and cost of skidding operations ([10], [23]).

The road density in the forest is currently very low, ranging between 1 and 2 m ha-1. To improve accessibility and facilitate timber harvesting in the forest, it has been proposed to increase the road density up to 20 m ha-1 ([12]). Therefore, an extensive road planning and construction will take place in the forest to facilitate management activities and timber harvesting that can produce required wood for the country.

Roads are an essential structure that provides access to the forest for a variety of activities such as logging operation, recreation and firefighting. The profitability of timber harvesting is largely influenced by road network structure ([22], [29], [6]), and thus it is important that roads are properly planned to ensure the cost-effective transport of forest products, as well as to provide safe and comfort vehicle operations ([3]). However, planning a proper road network to access timber harvest sites is a difficult and time-consuming task. Road networks are currently designed manually using contour and thematic maps ([16], [26]), but the manual approach might result in an inefficient network because only limited number of alternatives can be examined ([8], [1]). Therefore, the development of a methodology for the identification of current and future access needs, as well as the selection of the optimal road network from an economic point of viewpoint, are still a fundamental area of research ([22], [29], [6]). For the Hyrcanian forest, the need for tools and techniques to assist in the development of road network systems in the forest, particularly in the context of accessing timber harvest sites, it has been well identified ([12]).

The aim of this study is to estimate skidding costs from the results of a time study and incorporate the estimated road and skidding costs into an existing transportation planning tool, NETWORK 2000 ([7]), to select a cost effective road network that minimizes the total skidding and road building costs. We applied this method to the Gorazbon district in the Hyrcanian forest, where 244 links were built for the analysis. Detailed methods and application results are presented below.

# Material and methods

## Study Area

The study was carried out in the Gorazbon district in the Hyrcanian forest, which is a part of the Educational and Experimental Forest of University of Tehran (Fig. 1). It is located on northern slopes of Elbourz Mountains, about 7 km from Nowshahr (36° 33’ N, 50°33’ E). The district covers about 1000 hectares and its elevation ranges from 550 to 1380 m a.s.l. Ground slopes range between 0 to 115 percentages, and the average slope in the whole district is about 15 percentage. This area was dominated by natural forests containing native mixed deciduous tree species such as Fagus orientalis Lipsky, Carpinus betulus L. and Acer velutinum Boiss. A mixed uneven-aged management system with single and group selection harvesting has been used in this district, resulting in an average growing stock of 338 m3 ha-1 in a multistory forest structure. Selected trees for harvesting are felled by chainsaws, limbed and topped at the stump with a crew consisting of two men (chainsawer and coworker). Felled trees are bucked with chainsaws into saw logs (with diameters > 30 cm) and pulp logs (< 30 cm). Logs with a 5 to 15 meter length are extracted by wheeled skidders with winch to roadside landings and the fuel woods from delimibing parts are extracted by mules. In steep terrain where skidders are hindered, logs are processed to sawn-lumber and then hauled by mules. Sawn-lumbers (flitched lumbers) in the Hyrcanian forest region are 2.2-2.8 m long, 30-34 cm wide and 10-15 cm thick ([14]).

Fig. 1 - Study area located in the Hyrca- nian Forest.

The total length of extant roads in the Hyrcanian district is 12.45 km. Hence, to meet the target road density in this district, 6.92 km of new roads have been proposed for construction in the next decade. Current forest road density is 19.4 m ha-1. Forest roads are categorized as permanently main forest roads used for trucking, with an average width of 6.5 m and a longitudinal slope of 3 to 8% ([13]). Three compartments (301, 302 and 303), located on steep slopes and designated as protected areas, were excluded from our analysis.

## Data collection

A time study was conducted to estimate skidding costs and production rate of a Timberjack 450C skidder during June and July of 2011. This machine is normally an articulated, four-wheel-drive vehicle weighing 10.3 ton (55 % on the front and 45 % on the rear axle) with engine power of 177 hp (132 kW) and engine model 6BTA5.9. It is equipped with a blade for light pushing of obstacles and stacking of logs. The skidder was fitted with size 24.5-32 tires inflated to 220 kPa on both front and rear axles. Times and operational variables were measured using a stopwatch and recorded on paper ([9]). The working cycle for each operation consisted of certain elemental functions and factors. The times for each function and the value of each factor were recorded in the field. The time elements measured in the skidding cycle time include travel empty, releasing winch, choker setting, winching, travel loaded, unhooking and decking. Technical, personal and operational delays were recorded during the time study up to 15 minutes. In addition to measuring skidder working cycle time, independent variables expected to affect skidder productivity were documented. The variables measured during the time study include skidding distances (meter), skid trail slope (%), number of logs per turn, and load volume. A total of 150 working cycles were collected using the Timberjack 450C.

Timber harvest volume data and road costs were obtained from the forest management plan of the district. A total of 40 249 m3 of timber is planned to be harvested in the next 10 years (Tab. 1). For road costs, we estimated an annual life cycle cost (USD year-1) based on the total construction cost (Tab. 2), an interest rate of 5%, and a life cycle time period of 50 years. The interest on investment and the depreciation cost during the life cycle time period were calculated as USD 844 and USD 2.111 per year per km, respectively. The annual life cycle cost was thus estimated at USD 2.955 per km. Also, the truck transportation cost per cubic meter as the variable cost in the study area was assumed to be USD 0.31 per km.

Tab. 1 - Estimated harvest volume in the study area. (1): Refer to Fig. 2 for map of the harvest compartments. Some compartments are split into LN1 and LN2 depending on skidding direction. (2): Year 0 indicates the current year.

Harvest Comp.(1) Volume
(m3)
Year(2) Harvest Comp. Volume
(m3)
Year Harvest Comp. Volume
(m3)
Year
309LN1 1519.6 0 313LN2 758.8 2 323LN1 2201.7 6
307LN1 1349.7 0 317LN2 1090.8 3 327LN1 1201.2 6
307LN2 490.8 0 318LN2 990 3 322LN1 999 6
309LN2 681.2 0 318LN1 810 3 324LN1 2699.9 7
310LN1 844.2 1 317LN1 969.6 3 316LN1 1401.3 7
311LN1 632.7 1 320LN1 2298.3 4 319LN1 2499.9 8
312LN1 1200.0 1 314LN1 728 4 305LN1 445.5 8
310LN2 656.6 1 304LN1 650 4 305LN2 405 8
311LN2 266.4 1 304LN2 600 4 315LN1 1000 9
312LN2 600.0 1 314LN2 318.5 4 321LN1 2200 9
313LN1 1788.6 2 325LN1 699.2 5 306LN1 583.8 9
308LN1 1275.0 2 326LN1 2701.9 5 306LN2 166.8 9
308LN2 525.0 2 - - - - - -

Tab. 2 - Elements of road construction costs.

elements
Cost (USD
per km)
Planning 1 030.9
Earthwork (cut and fill) 9 278.4
Culverts and ditch 5 154.6
Maintenance and repair (5% of roading costs) 2 010.3
Total 42 216.5

## Using NETWORK 2000

Forest road network planning typically involves decision making on road locations that provides access to pre-determined harvesting sites (i.e., landings) in a way that minimizes overall road building and timber harvesting costs. Matthews ([17]) presented a theoretical model of the optimum road network that minimizes both skidding and road costs, but the model applications have been limited due to the fact that forest roads are not always parallel ([8]). Road network planning is often done manually, but some mathematical optimization techniques, such as mixed integer programming and heuristics, have been also used to solve road network location problems ([11]).

Fig. 2 - Harvest compartments and road links used in the network analysis.

# Results

## Skidding cycle time

Fig. 3 - A scatter plot showing the relationship between skidding distance and cycle time.

A regression analysis with the stepwise method between independent variables (skidding distance, trail slope, number of logs, and load volume) was performed on the time study data collected for Timberjack 450C skidder. The stepwise analysis revealed that skidding distance significantly affects the skidding cycle time (Fig. 3); thus a skidding cycle time regression model was developed using skidding distance as predicting variable (eqn. 1):

$$T = 6.43235 + 0.02686 \cdot D$$

where T is the skidding cycle time without delay (min) and D is the skidding distance (m). The multiple correlation coefficient was 0.79 and the variance accounted for by the model was 63% of the total variance. Overall goodness-of-fit was tested using a F-test (p<0.0001) and the significance of individual parameters by a t-tests (Tab. 3, Tab. 4).

Tab. 3 - ANOVA table for the regression model developed for the Timberjack 450C skidder.

Factor Sum of
Squares
df Mean
Square
F R Prob
Regression 3956.34 1 3956.34 258.5 0.79 <0.0001
Residual 2295.63 150 15.3 - - -
Total 6251.97 151 - - - -

Tab. 4 - Results of the t-test for the predicting variable.

Factor Unstandardized
Coefficients
Standardized
Coefficients
t Prob
Coeff. Std. Error Beta
Skidding distance 0.027 0.002 0.79 16.07 <0.0001
Constant 6.432 0.78 - 8.2 <0.0001

Travel loaded is the most time consuming phase of the skidding operation. Similar to travel unloaded, travel loaded is strongly related to skidding distances. The results showed that the skidding distance was an important variable to construct the time consumption model of skidding cycle time.

## Machine productivity

An hourly production rate (PR, m3 h-1) of the Timberjack 450C skidder was estimated using the production volume and time data collected during the time study as follows (eqn. 2):

$$PR = \frac {TP} {TT}$$

where TP is the total extracted volume (m3), and TT represents the total machine hours.

Fig. 4 - A scatter plot showing the relationship between skidding distance and hourly production rates (delay free time) of Timberjack 450C skidder.

The gross productivity (the average productivity with delay time) and the net productivity (the average productivity without delay time) of the Timberjack 450C for different distances were 8.3 m3 h-1 and 11.1 m3 ha-1, respectively. It was found that the production rate of the skidder decreases as skidding distance increases, but the decrease rate slows down as skidding distance exceeds 300 m (Fig. 4). A logarithmic regression model was chosen to describe the relationship between skidding distance and hourly production rates (delay free time - eqn. 3):

$$Y = -10.04 \,\ln(X) + 72.943$$

where Y is the net productivity (m3 h-1), and X is the skidding distance (m). Overall goodness-of-fit was tested by F-test (p<0.0001) and the significance of individual parameters by t-tests (Tab. 5, Tab. 6).

Tab. 5 - ANOVA table for the regression model developed for the production without delay time.

Factor Sum of
Squares
df Mean
Square
F R Prob
Regression 5163.76 1 5163.76 69.91 0.54 <0.0001
Residual 12311.61 150 82.08 - - -
Total 17475.37 151 - - - -

Tab. 6 - Results of t-tests for the predicting variable.

Factor Unstandardized
Coefficients
Standardized
Coefficients
t Prob
Coeff. Std. Error Beta
Skidding distance -10.04 1.266 -0.544 -7.932 <0.0001
Constant 72.943 7.534 - 9.681 <0.0001

## Estimating skidding costs

Estimates of hourly machine costs (USD ha-1) of the skidding operation were computed using the machine rate method ([21]). This machine cost was then divided by the estimated production rate (m3 ha-1) at a given skidding distance to estimate skidding cost per unit of volume (USD m-3). A purchase price of USD 123.711 was used in the machine cost estimating model, and the annual interest rate was set at 19%. A machine life of 10 years was assumed. Insurance and tax rate and utilization rate were set at 5% and 71.4%, respectively. Productive Machine Hour (PMH) and Scheduled Machine Hour (SMH) for the skidder were considered to be 900 hours and 1260 hours, respectively, so that the utilization of the chain saw is 75 % (eqn. 4). Cost calculations were based on the assumption that the machine operator works the whole year except the rainy season when the logging area is not accessible. The skidding operation was considered for a total 180 days per years (eqn. 4):

$$U = \frac{SMH}{PMH} = 71.4$$

Labor cost for skidding operation depends on the number of persons involved in each phase, the salary of each worker and the duration of their contracts. In Iran, almost all workers in the companies are paid monthly. Hourly cost derives from monthly salary divided by annual production hours. The skidding operation using the Timberjack 450 C typically involves a crew of three people (eqn. 5):

$$LC = \frac{\text{Total Labor Cost}} {\text{Productive Hour}} = \frac{USD\;12764} {900} = 14.18$$

Tab. 7 - Machine cost calculation results for Timberjack 450C.

Cost elements Total Cost
(USD)
Cost per hour
(USD)
Fixed costs - 12.37
Depreciation 11 134.02 15.13
Interest 13 617.53 2.75
Tax and Insurance 2 475.15 -
Subtotal (Fixed) 27 226.7 30.25
Operating costs - -
Maintenance and Repair 10 017 11.13
Fuel and Lubricant 5 013 5.57
Tires 3 294 3.66
Cable 2 772 3.08
Subtotal (Operating) 21 105 23.44
Hourly labor cost - 14.18
Total Hourly Machine Rate
( Fixed + Operating + Labor)
- 67.87

Tab. 7 summarizes machine cost estimate for Timberjack 450C skidder. The hourly machine cost was estimated as USD 67.87. Unit cost of production in different work phases was calculated by dividing the system cost by the average productivity per hour (eqn. 6):

$$UC = \frac{SC}{AP}$$

where UC is the unit cost (USD m-3), SC is the system cost (USD hour-1), and AP is the average productivity (m3 hour-1). Hence, the skidding cost per cubic meter of timber volume including the delay time was 8.2 USD m-3, while the skidding cost without delay was estimated as 6.1 USD m-3. Approximately 25% of total operating hours were classified as delay times during the time study, which results in an average net productivity of 75%. Fig. 5 show the distribution of elemental times per skidding with the Timberjack.

Fig. 5 - Percentage of time distribution of skidding elements.

## NETWORK 2000 output

The best solution obtained by the use of the software NETWORK 2000 indicates that 3.63 km of two new roads (one is from R3_2 to R3-6 and the other is from R1AS_2 to R1A_3 in Fig. 2) should be built to minimize both skidding and road costs. The total variable costs including skidding and haul costs is estimated as USD 340 631 for harvesting 40 250 m3 of timber in 37 compartments, while the estimated total road cycle time cost for next 10 years is USD 101 321. The average harvesting and road costs per m3 are USD 8.46 and USD 2.52, respectively.

# Discussion

In conclusion, we presented a methodology to support decision-making in timber harvest planning using skidding cost estimates and NETWORK 2000 software tool. We were able to identify cost effective landing and new road locations for the 37 harvest compartments in the Hyrcanian forest analyzed in this study. As suggested by Chung et al. ([8]), the heuristic network algorithm implemented in the NETWORK 2000 software does not guarantee optimal solutions, though trade-off analysis between road and skidding costs may be carried out, obtaining the overall least cost solution in road network and harvesting design. To this regard, NETWORK 2000 can be a valuable tool for forest managers who want to analyze road alternatives and select cost effective road and harvesting systems for timber management.

# References

(1)
Abdi E, Majnounian B, Darvishsefat AA, Mashayekhi Z, Sessions J (2009). A GIS-MCE based model for forest road planning. Journal of Forest Science 55 (4): 171-176.
Online | Gscholar
(2)
Abeli WS (1993). Comparing productivity and costs of three upgrading machines. International Journal of Forest Engineering 5 (1): 33-39.
CrossRef | Gscholar
(3)
Abeli WS, Shemwetta TK, Ole Meiludie REL, Kachwele M (2000). Road alignment and gradient issues in the maintenance of logging roads in Tanzania. International Journal of Forest Engineering 11 (2): 151-159.
Online | Gscholar
(4)
Akay AE, Eradas O, Sessions J (2004). Determining productivity of mechanized harvesting machines. Journal of Applied Science 4 (1): 100-105.
CrossRef | Gscholar
(5)
Bjorheden R, Thompson MA (1995). An international nomenclature for forest work study. In: Proceeding of the “XX IUFRO World Congress”. Tampere (Finland) 6-12 August 1995, pp.16.
Gscholar
(6)
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.
CrossRef | Gscholar
(7)
Chung W, Sessions J (2003). NETWORK 2000: a program for optimizing large fixed and variable cost transportation problems. In: Proceeding of the “24 Council of Forest Engineering Conference” (Wang J, Wolford M, McNeel J eds). West Virginia University, Morgantown, WV, USA.
Gscholar
(8)
Chung W, Stuckelberger J, Aruga K, Cundy TW (2008). Forest road network design using a trade-off analysis between skidding and road construction. Canadian Journal of Forest Research 38: 439-448.
CrossRef | Gscholar
(9)
Egan AF, Baumgas J (2003). Ground skidding and harvested stand attributes in Appalachain hardwood stands in West Virginia. Forest Product Journal 53 (9): 59-65.
Gscholar
(10)
Feghhi J (1989). Assessment of two mechanized forest harvesting systems. M.Sc. Thesis, Faculty of Natural Resources, University of Tehran, Iran, pp. 125. [in Persian]
Gscholar
(11)
Ghafarrian R, Sobhani H (2007). Optimization of an existing forest road network using NETWORK 2000. Croatian Journal of Forest Engineering 28 (2): 185-193.
Online | Gscholar
(12)
Iranian Plan and Budget Organization (2000). Forest road project guidelines no. 131 (2 edn). Office of the Deputy for Technical Affairs. Bureau of Technical Affairs and Standards, Tehran, Iran, pp. 106. [in Persian]
Gscholar
(13)
Jourgholami M (2010). Harvesting plan of Gorazbon district. Faculty of Natural Resources, University of Tehran, Iran, pp. 240. [in Persian]
Gscholar
(14)
Jourgholami M (2012). Small-scale timber harvesting: mule logging in Hyrcanian Forest. Small-scale Forestry 11(2): 255-262.
CrossRef | Gscholar
(15)
Lanford BL, Sobhani H, Stokes BJ (1990). Tree length loading production rates for southern pine. Forest Product Journal 33 (10): 120-131.
Gscholar
(16)
Liu K, Sessions J (1993). Preliminary planning of road systems using digital terrain models. International Journal of Forest Engineering 4: 27-32.
CrossRef | Gscholar
(17)
Matthews DM (1942). Cost control in the logging industry. McGraw-Hill, New York, USA, pp. 374.
Gscholar
(18)
McDonald T, Rummer B (2002). Automating time study of feller-buncher. In: Proceedings of the “33 Annual Meeting of Council of Afforests Engineering (COFE)”. Corvallis, OR, USA, pp. 17.
Gscholar
(19)
McDonald T (1999). Time study of harvesting equipment using GPS-derived positional data. In: Proceedings of the Conference “Forestry Engineering for Tomorrow”. Edinburgh University, Edinburgh, Scotland, pp. 12.
Gscholar
(20)
Meng CH (1984). A model for predicting logging machine productivity. Canadian Journal of Forest Research 14: 191-194.
CrossRef | Gscholar
(21)
Miyata ES (1980). Determining fixed and operational costs of logging equipment. General Technical Report NC-55, USDA Forest Service, St. Paul, MN, USA, pp. 16.
Gscholar
(22)
Murray AT (1998). Route planning for harvest site access. Canadian Journal of Forest Research 28 (7): 1084-1087.
CrossRef | Gscholar
(23)
Naghdi R (2004). Study of optimum road density in tree length and cut to length systems. Ph.D. Thesis, Faculty of Natural Resources, University of Tarbiat Modarres, Iran, pp. 225. [in Persian]
Gscholar
(24)
Pentek T, Picman D, Potocnik I, Dvorscak P, Nevecerel H (2005). Analysis of an existing forest road network. Croatian Journal of Forest Engineering 26 (1): 39-50.
Online | Gscholar
(25)
Picchio R, Spina R, Maesano M, Carbone F, Lo Monaco A, Marchi E (2011). Stumpage value in the short wood system for the conversion into high forest of a oak coppice. Forestry Studies in China 13 (4): 252-262.
CrossRef | Gscholar
(26)
Rogers L (2005). Automating contour-based route projection for preliminary forest road designs using GIS. M.Sc. Thesis, University of Washington, DC, USA, pp. 125.
Gscholar
(27)
Sobhani H, Stuart WB (1991). Harvesting systems evaluation in Caspian forests. International Journal of Forest Engineering 2 (2): 21-24.
CrossRef | Gscholar
(28)
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
(29)
Spinelli R, Nati C, Magagnotti N (2007). Recovering logging residues: experiences from the Italian eastern Alps. Croatian Journal of Forest Engineering 28: 1-9.
Online | Gscholar
(30)
Wang J, Haarlaa R (2002). Production analysis of an excavator-based harvester: a case study in finnish forest operation. Forest Product Journal 53 (2): 85-90.
Online | Gscholar
(31)
Wang J (2003). A computer-based time study system for timber harvesting operations. Forest Product Journal 53 (3): 47-53.
Online | Gscholar

#### Authors’ Affiliation

(1)
M Jourgholami
E Abdi
Department of Forestry and Forest Economics, Faculty of Natural Resources, University of Tehran, Karaj (Iran)
(2)
W Chung
Department of Forest Management, College of Forestry and Conservation, University of Montana, Missoula (USA)

#### Corresponding author

M Jourgholami
mjgholami@ut.ac.ir

#### Citation

Jourgholami M, Abdi E, Chung W (2013). Decision making in forest road planning considering both skidding and road costs: a case study in the Hyrcanian Forest in Iran. iForest 6: 59-64. - doi: 10.3832/ifor0640-006

#### Paper history

Accepted: Nov 06, 2012

First online: Jan 21, 2013
Publication Date: Apr 02, 2013
Publication Time: 2.53 months

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

#### Breakdown by View Type

(Waiting for server response...)

#### Article Usage

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

Breakdown by View Type
HTML Page Views: 13675
Abstract Page Views: 382

Web Metrics
Days since publication: 2495
Overall contacts: 18231
Avg. contacts per week: 51.15

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

(No citations were found up to date. Please come back later)

#### 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