<<

Article Performance Analysis of Log Extraction by a Small Operation in Steep Forests of South Korea

Eunjai Lee 1 , Sang-Kyun Han 2 and Sangjun Im 3,*

1 National Institute of Forest Science, Forest Technology and Management Research Center, 498 Gwangneungsumogwon-ro, Soheul-eup, Pocheon 11186, Korea 2 Department of and Landscape Architecture, Korea National College of Agriculture and Fisheries, 1515 Kongjwipatwji-ro, Deokjin-gu, Jeonju 54874, Korea 3 Department of Forest Sciences and Research Institute of Agriculture and Life Sciences, Seoul National University, 1 Gwanak-ro, Gwanak-gu, Seoul 08826, Korea * Correspondence: [email protected]; Tel.: +82-2-880-4768

 Received: 7 May 2019; Accepted: 11 July 2019; Published: 13 July 2019 

Abstract: In South Korea, logs for low-value products, such as pulpwood and fuelwood, are primarily extracted from harvest sites and transported to roadside or landing areas using small . Previous studies on log extraction, however, have focused on cable yarding operations with the goal of improving productivity on steep slopes and inaccessible sites, leaving small-shovel operations relatively unexamined. Therefore, the main objectives were to determine small-shovel extraction productivity and costs and to evaluate the impact of related variables on productivity. In addition, we developed a model to estimate productivity under various site conditions. The study took place in 30 case study areas; each area has trees with stems at a diameter at breast height ranging from 18 to 32 cm and a steep slope (greater than 15%). The areas ranged from 241 to 1129 trees per hectare, with conifer, deciduous, and mixed stands. Small-shovel drives ranged from 36 to 72 m per extraction cycle from stump to landing. The results indicated that the mean extraction productivity of small-shovel operations ranged between 2.44 to 9.85 m3 per scheduled machine hour (including all delays). At the forest level, the estimated average stump-to-forest road log production costs were US $4.37 to 17.66/m3. Small-shovel productivity was significantly correlated with stem size (diameter at breast height and tree volume) and total travelled distance (TTD). However, a Pearson’s correlation analysis indicated that stand density and slope did not have a significant effect on productivity. Our findings provide insights into how stem size and TTD influence small shovel performance and the predictive ability of productivity. Further, this information may be a valuable asset to forest planners and managers.

Keywords: shovel ; follow-up study; cut-to-length extraction; productivity; cost; scheduled machine hour

1. Introduction Logging operations (a.k.a., primary transportation), in which logs are transported from stumps to a designated roadside or central landing area using various extraction methods, are an important part of the timber harvesting process, but they can be extremely expensive and more time-consuming in practice than felling and processing [1–4]. Many studies have examined the performance of extraction methods on sites with differing characteristics in order to establish the logistics of extraction activities, such as ground-based extraction (skidding and forwarding) [5–8] and cable yarding [8–10]. It can be concluded from these studies that logging practices should be economically determined while maintaining a deep understanding of both the potential and limitations of the chosen method.

Forests 2019, 10, 585; doi:10.3390/f10070585 www.mdpi.com/journal/forests ForestsForests 2019,, 10,, 585x FOR PEER REVIEW 22 of of 14

Time-and-motion studies have been widely used to evaluate the performance of individual loggingTime-and-motion equipment as well studies as entire have harvesting been widely systems used [11, to evaluate12]. This type the performanceof study has been of individual essential forlogging predicting equipment machine as well productivity as entire harvestingand utilizatio systemsn rates [11 in, 12various]. This scenarios type of study under has similar been essentialworking conditionsfor predicting [13, machine14]. However, productivity this approach and utilization is limited rates in terms in various of data scenarios availability under due similarto the relatively working shortconditions period [13 of,14 data]. However, collection this and approach high costs is limited of field in work terms [14 of– data16]. Additionally availability due, a number to the relatively of past studiesshort period have ofassess dataed collection performance and high using costs the of follow field work-up method [14–16]., Additionally, which utilizes a numberhistoric ofoutput past records,studies have such assessed as productivity performance and usingcosts. theThis follow-up can provide method, more which accurate utilizes information historic output on long records,-term performancesuch as productivity than time and-and costs.-motion This canstudies provide [17–1 more9]. accurate information on long-term performance thanIn time-and-motion the Republic of studies Korea (a.k.a., [17–19 ].South Korea), forest land covers an area of 6.3 million hectares (64%In of thethe Republic total land of area), Korea and (a.k.a., approximately South Korea), 80% forest of all land forested covers areas an area are ofon 6.3 steep million terrain hectares with (64%slopes of greater the total than land 40% area), [20 and]. In approximately 2017, the stand 80% density of all forested was 154.1 areas m are3/h ona, withsteep conifers terrain with being slopes the 3 dominantgreater than species, 40% [20 namely]. In 2017, the theKorean stand red density pine was (Pinus 154.1 densiflora m /ha,), with Japanese conifers larch being (Larix the dominant kaempferi [Lamb.]species, namelyCarrière), the and Korean pitch red pine pine (Pinus (Pinus rigida densiflora Mill). ),The Japanese volume larch of timber (Larix harvested kaempferi [Lamb.] has considerably Carrière), increasedand pitch pinein recent (Pinus years: rigida fromMill). 1.3 The million volume m3 ofin timber2013 to harvested 2.2 million has m considerably3 in 2017 [20].increased in recent 3 3 years:Two from harvesting 1.3 million methods m in 2013 are to 2.2commonly million m usedin 2017 in South [20]. Korea: cut-to-length small-shovel productionTwoharvesting (CTL-S) methodsand tree- arelength commonly cable-yarder used in production South Korea: harvesting cut-to-length (TL- small-shovelC). There are production a number (CTL-S)of past studies and tree-length that investigate cable-yarder the TL production-C method to harvesting describe, (TL-C). understand, There and are improve a number upon of past the efficiencystudies that of investigate log production the TL-C and method associated to describe, operational understand, decisions. and Other improve studies upon focusthe effi onciency the productivityof log production and and operation associated efficiency operational of individual decisions. machines Other studies [21,22 focus], comparing on the productivity extraction performancesand operation among efficiency different of individual cable-yarder machines technologies [21,22], comparing [23,24], and extraction the effect performances of yarding direction among (uphilldifferent vs. cable-yarder downhill) technologies on productivity [23,24 and], and cost the [10, eff2ect5] ofin yarding order to direction support (uphill operational vs. downhill) decisions. on productivityHowever, the and use cost of the [10 ,TL25]-C in method order to for support log hauling operational activities decisions. remains However, limited thesince use this of themethod TL-C requiresmethod fornot log only hauling a high activities skill level, remains but also limited an inherently since this methodhigh level requires of investment not only [ a2 high6]. Thus, skill cable level, yardingbut also systems an inherently have rarely high level been of implemented investment [on26]. steep Thus, slopes cable and yarding remote systems areas havedue to rarely operation been costimplementeds. on steep slopes and remote areas due to operation costs. In South South Korea, Korea, CTL CTL-S-S is isthe the preferred preferred system, system, and and it has it replaced has replaced TL-C TL-Con steep on steepslopes. slopes. After Afterfelling felling and processing and processing trees with trees a , with a chainsaw, 2–4 m logs 2–4 are m logsextracted are extracted using a small using shovel, a small which shovel, is 3 whicha small is-sized a small-sized hydraulic hydraulicexcavator (5.0 excavator metric (5.0tons metricin weight tons with in weight a 0.2 m with3 bucket) a 0.2 with m abucket) log grapple with (i.e.,a log a grapple small shovel (i.e., a or small wood shovel grab; or Figure wood 1). grab; Generally, Figure1 ).small Generally,-shovel small-shovel(SS) extraction (SS) activities extraction use gravitationalactivities use energy gravitational to transport energy logs to from transport the stump logs fromto the the roadside stump/landing to the roadside area. For/landing example, area. on steepFor example, terrain, on gravity steep terrain, can be gravity useful in can assisting be useful with in assisting throwing, with sliding, throwing, and sliding, rolling and logs rolling downhill logs [3,downhill27]. Thus, [3,27 the]. Thus, small the shovel small- shovel-basedbased logging logging method method could could be beincreasingly increasingly app appliedlied in in Korea’s Korea’s harvesting operations, operations,but but the the productivity productivity level level and and costs costs associated associated with with extracting extracting 2–4 m 2– logs4 m using logs usingSS activities SS activities remains remain unclear.s unclear.

Figure 1. LogLog extraction extraction using using a a small small shovel shovel-mounted-mounted steel steel-track-track excavator with a log grapple (1.3 m mmaximumaximum jawjaw openingopening and and 0.1 0.1 m m closed closed jar jar gap); gap); the the boom boom length length of the of excavator the excavator is 5.1 mis (Photo5.1 m ( create:Photo create:E. Lee). E. Lee).

Therefore, in this study, the overall objective was to determine the performance of SSs in various types of forest. In particular, this study sought to: (1) determine the productivity (m3/SMH) and costs

Forests 2019, 10, 585 3 of 14

Therefore, in this study, the overall objective was to determine the performance of SSs in various Forests 2019, 10, x FOR PEER REVIEW 3 of 14 types of forest. In particular, this study sought to: (1) determine the productivity (m3/SMH) and costs 3 (US(US $$/m/m 3) of extracting logs through the follow-upfollow-up method,method, (2) establish the influentialinfluential variables in SS extraction productivity, productivity, and and (3) (3) develop develop regression regression models models to topredict predict SS productivity. SS productivity. Further, Further, the theresults results of this of this study study will will lead lead to to better better-informed-informed SS SS technology technology decisions decisions and and more more eefficientfficient production of timber products.products.

2. Materials and Methods In collaboration with a group of logging companies and contractors under the KKoreaorea Wood Products Association (KWPA), we conducted a follow-upfollow-up study o off SS extraction activity, which is becoming thethe prevailingprevailing extractionextraction technology. technology. The The study study focused focused on on the the production production of CTLof CTL clear-cut clear- harvestcut harvest units units (CHUs) (CHUs) in the in the Central Central Northeast Northeast region regi ofon Southof South Korea Korea (Gangwon-do, (Gangwon-do, Gyeonggi-do, Gyeonggi- Chungcheongbuk-do,do, Chungcheongbuk-do, and and Gyeongsangbuk-do; Gyeongsangbuk-do; Figure Figure2). The 2). The main main characteristics characteristics of the of CHUsthe CHUs are presentedare presented in Table in Table1. The 1. unitsThe units are located are located on relatively on relatively steep steep terrain terrain (ranging (ranging from 13from to 64%,13 to with64%, anwith average an avera slopege slope of 49%); of 49%); trees trees in these in these areas areas have have a DBH a DBH (diameter (diameter at breast at breast height) height) of up of toup 32 to cm 32 (minimumcm (minimum DBH, DBH, 18 cm). 18 cm). On average,On average, the unitsthe units have have 560 trees560 trees per hectare per hectare (TPH), (TPH), and the and range the range is from is 241from to 241 1129 to TPH.1129 TPH. The data The set data covered set covered the dominant the dominant forest standforest types:stand types: conifer, conifer, deciduous, deciduous, and mixed. and Themixed. total The traveled total traveled distance distance (TTD; a.k.a., (TTD; total a.k.a., driven total distance)driven distance) ranged fromranged 36 from to 72 36 m andto 72 average m and roadaverage density road was density 108 mwas/ha 108 (ranged m/ha from (ranged 32 to from 188 m32/ha). to 188 We m/ha) defined. We TTD define for small-shoveld TTD for small operations-shovel asoperation the distances as the the distance extractor the travels, extractor starting travels, when starting the small when shovel the small leaves shovel the forestleaves road the forest or landing road areaor landing and ending area and when ending it returns when to it thereturn landings to the (Figure landing3). The(Figure study 3). areasThe study were areas selected were to selected cover a wideto cover range a wide of forest range conditions; of forest conditions; this was done this to was achieve done an to enhanced achieve an understanding enhanced un ofderstanding small shovel of extractionsmall shovel performance extraction performance across South Korea.across South Korea.

Figure 2. Map of study sites in South Korea showing the harvesting units.

Forests 2019, 10, 585 4 of 14

Table 1. Summary of information collected for each unit: post-harvesting stand characteristics, average total travelled distance, and productivity.

Average Harvest Stand Average 3 Stand Average Average Total Productivity a Volume (m ) Unit No. Type DBH Density Slope Travelled (m3/day) 1 1 (cm) tree− ha− (Trees/ha) (%) Distance (m) 1 Deciduous 22 0.26 112 435 64 58 50.9 2 Deciduous 22 0.26 119 466 60 60 58.4 3 Conifer 20 0.20 155 756 42 64 59.8 4 Conifer 22 0.24 271 1129 52 68 58.6 5 Conifer 24 0.34 148 436 42 44 44.1 6 Conifer 18 0.14 155 1,107 64 42 26.3 7 Conifer 24 0.30 188 627 32 69 75.5 8 Deciduous 32 0.56 157 278 64 55 52.7 9 Conifer 22 0.25 190 760 32 61 60.3 10 Deciduous 18 0.14 105 755 64 43 19.5 11 Mix 22 0.25 113 450 42 45 40.2 12 Mix 18 0.15 96 650 13 43 47.0 13 Mix 22 0.23 97 418 32 58 66.8 14 Mix 24 0.30 97 320 32 55 55.8 15 Deciduous 24 0.28 68 241 52 56 57.4 16 Deciduous 24 0.29 138 476 42 45 33.7 17 Deciduous 20 0.19 83 442 64 46 35.0 18 Deciduous 22 0.23 69 296 60 46 45.1 19 Conifer 18 0.16 118 716 42 53 26.1 20 Deciduous 24 0.28 72 255 42 42 37.7 21 Mix 18 0.18 161 909 40 45 42.4 22 Mix 32 0.60 169 282 64 72 78.8 23 Mix 32 0.56 162 287 60 65 71.8 24 Mix 20 0.19 121 626 52 36 33.8 25 Deciduous 18 0.16 134 819 64 61 50.9 26 Deciduous 24 0.27 124 459 52 39 43.3 27 Deciduous 18 0.15 62 408 52 41 22.8 28 Mix 24 0.27 154 580 64 44 43.4 29 Mix 22 0.22 77 344 22 45 48.9 30 Conifer 26 0.36 201 558 64 55 60.7 a DBH: Diameter at breast height.

CTL-S clear-cutting operations in each unit were performed using a semi-mechanized system that employs a chainsaw for felling, delimbing, and bucking trees into 2–4 m logs, which are mostly used as pulpwood. (An alternative to the use of the SS, especially on steep slopes, is manipulating all logs at the stump, as described by Lee et al. [10,28].) The SS operation utilized the gravity extraction technique, which involves throwing, rolling, and pushing logs, to move the logs to the roadside or landing area. When SS travels up and down, the machine moves along the slope direction. For each unit, timber harvesting had a similar target: (1) support sustainable forest management to ensure future availability and (2) benefit domestic timber industries by increasing the self-sufficiency rate. Three years of historical data, from the period of 2015 to 2017, were manually collected from logging companies. The collected dataset includes a detailed description of harvest unit characteristics and the net production rate (m3/day, based on an 8-hour work day) of log extraction by SS operation (Table1). SS extraction productivity was determined by two processes: (1) the sorting of 2–4m logs when the SS commences travel into the felling site from the roadside/landing area and (2) the throwing and pushing of logs while the SS returns to the roadside/landing area. Forests 2019, 10, 585 5 of 14 Forests 2019, 10, x FOR PEER REVIEW 4 of 14

Figure 3. Vertical pattern for small shovel extraction activity cut-to-length to forest road/landing area Figure 3. Vertical pattern for small shovel extraction activity cut-to-length to forest road/landing area and for operation in steeply sloped forests in South Korea. and for operation in steeply sloped forests in South Korea. In this study, the follow-up data collection method conducted did not involve a time-and-motion measuringCTL-S clear device.-cutting This operations study relied in solely each onunit historical were performed data from extractionusing a semi operations-mechanized performed system that employsby forest contractors.a chainsaw for There felling, was nodelimbing, information and available bucking on trees delay into times, 2–4 m including logs, which mechanical, are mostly usedoperational, as pulpwood. or personal (An alternative delays. Therefore, to the use scheduled of the SS, machine especially hours on (SMH) steep were slopes, used is to manipulating evaluate the all logs productivityat the stump and, as cost described of SS operations. by Lee et al. [10,28].) The SS operation utilized the gravity extraction technique,The which costs forinvolves owning throwing, and operating rolling, SSs were and calculated pushing with logs, the to method move developedthe logs to by the Miyata roadside [29], or landingwhich area. is a When commonly SS travel accepteds up machine and down, rate calculation the machine technique. moves In along addition theto slope costs direction associated. withFor each machine ownership, operation, and labor, we included machine delays (mechanical, operational, and unit, timber harvesting had a similar target: (1) support sustainable forest management to ensure personal) and warm-up costs. This is necessary because the productivity data, which is divided into future availability and (2) benefit domestic timber industries by increasing the self-sufficiency rate. two categories: machine operation and idle time, is per SMH. The machine utilization rate, labor cost, andThree fuel years consumption of historical rate were data collected, from the from period the KWPA of 2015 (Table to 22017). The, were overhead, manually profit allowance,collected from loggingand transportation companies. costs The associated collected with data SSset were includes not obtained. a detailed description of harvest unit characteristicsThe SPSS and packagethe net (IBMproduction Co., New rate York, (m3/day, NY, USA,based v. on 22.0) an 8 was-hour used work for day) statistical of log analysis. extraction by SSPearson’s operation correlation (Table 1). test SS wasextraction conducted productivity to clarify how was the determined independent by variablestwo processes: (DBH, slope,(1) the and sorting of 2–TTD)4m logs affect when SS productivity. the SS commences Based on thetravel value into of the correlationfelling site coe fromfficients, the roadside a predictive/landing equation area of and (2) theproductivity throwing per and SMH pushing was developed of logs while using the ordinarySS returns least to squares the roadside regression/landing technique. area. Two-thirds ofIn thethis follow-up study, the data follow were-up randomly data collection selected for method model conducted development, did while not involve the remaining a time one-third-and-motion measuringof the data device was.applied This study for validation relied solely of the proposedon historical model. data A two-samplefrom extraction t-test was operations used to compare performed by forestthe predicted contractors. and observed There was values no and information to describe anyavailable statistical on di delayfferences times, for model including verification. mechanical, operational, or personal delays. Therefore, scheduled machine hours (SMH) were used to evaluate the productivity and cost of SS operations. The costs for owning and operating SSs were calculated with the method developed by Miyata [29], which is a commonly accepted machine rate calculation technique. In addition to costs associated with machine ownership, operation, and labor, we included machine delays (mechanical, operational, and personal) and warm-up costs. This is necessary because the productivity data, which is divided into two categories: machine operation and idle time, is per SMH. The machine utilization rate, labor cost, and fuel consumption rate were collected from the KWPA (Table 2). The overhead, profit allowance, and transportation costs associated with SS were not obtained.

Forests 2019, 10, x FOR PEER REVIEW 6 of 14

Table 2. Cost components and estimated hourly cost to own and operate a small shovel.

Cost Component Small Shovel Purchase price (US $) 54,000.00 Salvage value (%) 20 Economic life (year) 7 Scheduled machine hour per year 1400 Interests (%) 10 Insurance (%) 3 Forests 2019, 10, 585 Taxes (%) 2 6 of 14 Fuel consumption rate (liter/hour) 9.0 Fuel cost (US $/liter) 1.20 Table 2. Cost components and estimated hourly cost to own and operate a small shovel. Lubrication (% of fuel cost) 40 RepairCost Component and maintenance (% of depreciation)Small Shovel90 LaborPurchase (US price $/hour) (US $) 54,000.0017.00 FringeSalvage benefit value (%) (% of labor) 20 22 FixedEconomic cost life (US (year) $/hour) 7 8.21 OperatingScheduled machine cost (US hour $/hour) per year 140014.52 LaborInterests cost (%) (US $/hour) 1020.35 Insurance (%) 3 TotalTaxes (%)operation cost (US $/hour) 245.87 Fuel consumption rate (liter/hour) 9.0 The SPSS packageFuel cost(IBM (US Co., $/liter) New York, NY, USA, v. 22.0) was 1.20 used for statistical analysis. Pearson’s correlationLubrication test was conducted (% of fuel cost) to clarify how the independent 40 variables (DBH, slope, and TTD) affect SS productivity.Repair and Based maintenance on the value (% of depreciation) of the correlation coefficients, 90 a predictive equation Labor (US $/hour) 17.00 of productivity per SMHFringe was benefit developed (% of labor) using the ordinary least squares 22 regression technique. Two- thirds of the follow-Fixedup data cost were (US $ /randomlyhour) selected for model development, 8.21 while the remaining one-third of the dataOperating was applied cost for (US validation $/hour) of the proposed model. 14.52 A two-sample t-test was used to compare the predictedLabor costand (US observed $/hour) values and to describe any 20.35statistical differences for model verification. Total operation cost (US $/hour) 45.87

3.3. Results Results TheThe study study was was designed designed to to test test the the effi efficacycacy of of SS SS extraction extraction from from stumps stumps to to roadside roadside through through 30 30 didifferefferentnt case case studies. studies. Overall, Overall, the the SS SS was was capable capable of of a productivitya productivity rate rate of of 2.44 2.44 to to 9.85 9.85 m m3/SMH3/SMH at at a a costcost of USof US$4.37 $4.37 to 17.66 to /17.66m3 (Figure/m3 (Figure4). There 4). were There large were variations large in variations productivity in productivity and cost evaluations and cost acrossevaluations the harvest across units. the harvest units.

Figure 4. Box and whisker plots of productivity between the current study and previous studies: Kim and Park [30] and Lee et al. [31]; both studies evaluated the shovel logging productivity. The X indicates mean value.

Productivity may increase or decrease with variation in stem size, such as DBH and tree/log volume. Through Pearson’s correlation test, we found that DBH and tree volume had a considerable impact on productivity (p < 0.001). Productivity and stem size, including DBH (r = 0.6168) and tree volume (r = 0.6161), were directly related (Figures5 and6). Forests 2019, 10, x FOR PEER REVIEW 7 of 14 Forests 2019, 10, x FOR PEER REVIEW 7 of 14 Figure 4. Box and whisker plots of productivity between the current study and previous studies: Kim Figureand Park 4. Box [30 ]and and whisker Lee et plots al. [31 of] ;productivity both studies between evaluated the thecurrent shovel study logging and previous productivity. studies: The Kim X andindicates Park mean [30] andvalue. Lee et al. [31]; both studies evaluated the shovel logging productivity. The X indicates mean value. Productivity may increase or decrease with variation in stem size, such as DBH and tree/log volume.Productivity Through mayPearson’s increase correlation or decrease test, wewith found vari ationthat DBH in stem and size, tree volumesuch as hadDBH a considerableand tree/log volume.impact on Through productivity Pearson’s (p < correlation0.001). Productivity test, we found and stem that DBHsize, includingand tree volume DBH ( rhad = 0.6168) a considerable and tree Forestsimpactvolume2019 on (r, 10 =productivity ,0. 5856161), were (p directly < 0.001). related Productivity (Figures and 5 and stem 6). size, including DBH (r = 0.6168) and7 tree of 14 volume (r = 0.6161), were directly related (Figures 5 and 6).

Figure 5. Productivity, which was calculated using follow-up data from 30 study areas, and its Figurerelationship 5. Productivity,Productivity, with DBH (Diameter which which was was at calculated calculatedBreast Height). using using Productivity follow follow-up-up data data had from froma moderate 30 study correlation areas, and with its relationshipDBH. withwith DBH DBH (Diameter (Diameter at at Breast Breast Height). Height). Productivity Productivity had had a moderate a moderate correlation correlation with DBH.with DBH.

Figure 6. 6. Productivity,Productivity, which which was was calculated using followfollow-up-up data from 30 study areas, and its Figurerelationship 6. Productivity, withwith treetree volume.volume. which Productivity was calculated had using a moderate follow correlationcorrelation-up data from withwith 30treetree study volume.volume. areas, and its relationship with tree volume. Productivity had a moderate correlation with tree volume. In most harvest harvest units, units, a ahigher higher stand stand density density is significantly is significantly associated associated with with small small stem stem size size(p < (0.001;p < 0.001In DBH most; DBH ofharvest r of= r−0.5830= units,0.5830 anda higher and tree tree standvolume volume density of ofr = ris =− significantly0.5322).0.5322). However, However, associated stand stand with density density small had hadstem a a weaksize weak (p to < − − 0.001;moderate DBH negative of r = − correlationcorrelation0.5830 and withwith tree SSSSvolume productivityproductivity of r = (− (r0.5322).r= = −0.2214;0.2214; However, Figure 7 7stand).). AsAs aadensity result,result, SSSShad productivityproductivity a weak to − wasmoderate inversely negative correlatedcorrelated correlation withwith standstandwith SS density.density. productivity (r = −0.2214; Figure 7). As a result, SS productivity was inversely correlated with stand density.

Forests 2019, 10, 585 8 of 14

Forests 2019, 10, x FOR PEER REVIEW 8 of 14

Figure 7. Productivity,Productivity, which which was was calculated calculated using using follow follow-up-up data data from from 30 30 study areas, and its relationship with stand density. Productivity hadhad nono correlationcorrelation withwith standstand density.density.

TThishis study implied that an increased TTD would be positively and significantly significantly correlated with SS productivity (p < 0.001;0.001; rr == 0.8262;0.8262; Figure Figure 8). On the other hand, slope had no significant significant correlation with productivity ((rr == −0.1060;0.1060; Figure Figure 9).). Although our data is limited to divide intointo SSSS cyclecycle elementalelemental with productivity (r = −0.1060;− Figure 9). Although our data is limited to divide into SS cycle elemental time, we found that the time spent throwing and rolling logs may be considerably longer than the travel time time.. Further, Further, in in terms terms of of extraction extraction activity, activity, there was no common pattern across Korea for selecting SS operations in steeply-slopedsteeply-sloped forests.forests.

Figure 8. Productivity,Productivity, which which was was calculated calculated using using follow follow-up-up data data from from 30 30 study areas, and its relationship with total travelled distance. Productivity Productivity had had a a strong correlation with average total travelled distance.

Forests 2019, 10, 585 9 of 14 Forests 2019, 10, x FOR PEER REVIEW 9 of 14

Figure 9. Productivity,Productivity, which which was was calculated calculated using using follow follow-up-up data data from 30 study areas, and its relationship with slope. Productivity had no correlation with slope.

The SSSS productivity productivity regression regression equation equation was derivedwas derived from the from follow-up the follow data-up to predict data to extraction predict productivityextraction productivity (Table3). The (Table independent 3). The variables independent ranged variables from 0.14 ranged to 0.56 from m 3 for0.14 tree to volume0.56 m3 and for 36 tree to volume69 m for and TTD. 36 This to 69 model m for was TTD tested. This formodel assumptions was tested of for normality, assumptions independence, of normality, and independence, equal variance andto confirm equal variance the validity to confirm of the analysis.the validity The of TTD the analysis. was significant The TTD (p 0.05). Thus,value the(p > obtained 0.05). Thus, model the mayobtained be quite model accurate. may be quite accurate. Table 3. Productivity regression model for small-shovel operations extracting 2–4 m logs. Productivity Tableis in m 3 3./SMH Productivity (Scheduled regression Machine Hour).model forA paired small- t-testshovel was operations used for model extracting verification 2–4 m across logs. Productivityobserved data. is in m3/SMH (Scheduled Machine Hour). A paired t-test was used for model verification across observed data. Model Model t-Test Average Productivity Estimator SE t p-Value 2 adj. R p-ValueModel (tp-Test-Value) p- Model adj. Average= 2.1861 Productivity Estimator 0.9966SE t 3.0810 0.0396 0.7517 <0.01p- 0.2935(p- − − Value R2 Value Value) + 3.6873 Average tree volume (m3) 1.8400 2.2827 0.0581 = −2.1861 × 0.9966 −3.0810 0.0396 0.7517 <0.01 0.2935 + 0.1404 Average total travelled + 3.6873 ×× Average tree volume (m3) 0.02051.8400 2.2827 6.9908 0.0581<0.01 distance (m) + 0.1404 × Average total travelled 0.0205 6.9908 <0.01 distance (m) The SS productivity model showed that productivity improved as tree volume and travelled distanceThe increased;SS productivity the average model treeshowed volume that rangedproductivity from 0.15improved to 0.50 as m tree3, and volume the mean and TTDtravelled was distancebetween 40increased; and 70 m the (Figure average 10). tree The datavolume indicated ranged that from when 0.15 the to volume 0.50 m of3, and trees the extracted mean TTD increased was betweenfrom 0.15 40 m 3 andto 0.50 70 m m3 (Figure(and the 10 TTD). The ranged data between indicated 40 andthat 70 when m), productivity the volume increased of trees extracted by up to increased30%. This from result 0.15 may m be3 to explained 0.50 m3 (and by the the fact TTD that ranged efficiency between improved 40 and with 70 m), larger productivity tree volumes increased and an byincreased up to 30%. number This of result logs at may the stumpbe explained area. by the fact that efficiency improved with larger tree volumes and an increased number of logs at the stump area.

Forests 2019, 10, 585 10 of 14 Forests 2019, 10, x FOR PEER REVIEW 10 of 14

Figure 10. Productivity for the small-scalesmall-scale extraction operations in relation to treetree volumevolume (m(m33) and total travelled distance (m).(m).

4. Discussion Small-shovelSmall-shovel logging method has become widespread and has had increasing use across South Korea in steeply-slopedsteeply-sloped forests. With steep terrain extraction, the slope grade helps in hauling the logs downhill with the support of gravitygravity from the fellingfelling areaarea [[28,28,3232].]. We conducted a follow-upfollow-up study of the SS extraction method to evaluat evaluatee the productivity and cost under various types of forests and established the influential influential variables variables in in extraction extraction productivity. productivity. In addition, we developed a regression model model for for the the SS toSS estimate to estimate hourly hourly productivity productivity in SMH. in The SMH. results The show results that show the estimated that the 3 stump-to-forestestimated stump road-to-forest log productivityroad log productivity was between was between 2.44 and 2.44 9.85 and m 9.85/SMH m3/SMH at US at $17.66 US $17.66 and 4.37,and respectively.4.37, respectively. The stem The size stem and size TTD and statistically TTD statistically have an ha influenceve an influence on the productivity on the productivity of SS extraction of SS operation (p < 0.001), but the slope had no static correlation (r = 0.1060). Thus, the model showed extraction operation (p < 0.001), but the slope had no static correlation− (r = –0.1060). Thus, the model thatshowed productivity that productivity improved improved as tree volume as tree and volume TTD increased. and TTD increased. Further, we Further, examined we a examined sensitivity a analysissensitivity to analysisevaluate tothe evaluate impact of the the impact tree volume of the and tree TTD volume on the and model. TTD on The the result model. indicated The result that, 3 whenindicated the volume that, wh ofen trees the extracted volume of increased trees extracted from 0.15 increased to 0.50 m from, with 0.15 the to TTD 0.50 ranging m3, with between the TTD 40 andrang 70ing m, between productivity 40 and increased 70 m, productivity by up to 30%. increased by up to 30%. 3 In our study, mean SS SS extraction extraction productivity productivity was was 6.03 6.03 ± 1.901.90 m m3/SMH,/SMH, which which is isvery very similar similar to ± 3 tothe the findings findings by by Kim Kim and and Park Park [30 [30] and] and Lee Lee et et al. al. [31 [31],], whose whose results results were were 5.21 5.21 and 6.57 mm3/SMH,/SMH, respectively (Figure4 4).). TheseThese previousprevious studiesstudies likelylikely operatedoperated underunder similar similar conditions: conditions: anan averageaverage 3 DBH of 22 toto 2626 cm,cm, treetree volumevolume ofof 0.240.24 toto 0.270.27 mm3,, andand slopeslope ofof 3636 toto 44%.44%. Thus, the follow-upfollow-up evaluation isis deemeddeemed reasonablereasonable and and acceptable acceptable in in terms terms of of accurately accurately determining determining the the productivity productivity of theof the SS SS extraction extraction process. process. The productivityproductivity andand costs costs of of SSs SSs varied varied when when extracting extracting 2–4 2 m–4 logs m logs in steep in steep terrain terrain due todue a wide to a rangewide range of forest of forest types. types Numerous. Numerous studies studies have pointed have pointed out that out di thatfferences differences in productivity in productivity and cost and of extractioncost of extraction activities activities were due were to due locally to locally variable variable conditions, conditions, including including the volume the volume of extracted of extracted logs, stemlogs, stem size, extractionsize, extraction distance, distance, and workingand working conditions conditions [33–35 [3].3– Our35]. dataOur data significantly significantly showed showed that stemthat stem size andsize TTDand TTD were were important important key factorskey factor in SSs in productivity SS productivity prediction, prediction, while while the slopethe slope was was not statisticallynot statistically correlated. correlated. ShovelShovel logging logging productivity productivity is affected is affected by the by number the number of swings of swingsand road and spacing road [3 spacing6]. Although [36]. Althoughour data was our limited data was to determine limited to the determine effect of number the effect of ofswings number and ofroad swings spacing and, we road tested spacing, how wethe testedDBH, howslope, the and DBH, TTD slope, affect andSS productivity TTD affect SS using productivity Pearson’s using correlation. Pearson’s We correlation. found that Westem found size thathad a stem statistical size had impact a statistical on productivity. impact onNumerous productivity. studies Numerous have pointed studies out that have stem pointed size has out been that stema primary size hasfactor been in aextraction primary factorproductivity in extraction [11,34,37, productivity38], because [11 a,34 larger,37,38 log], because size increases a larger the log piece size volume, payload, and productivity. Therefore, DBH and tree volume may be two of the main factors influencing SS productivity. Further, tree volume, instead of DBH, was selected in the regression

Forests 2019, 10, 585 11 of 14 increases the piece volume, payload, and productivity. Therefore, DBH and tree volume may be two of the main factors influencing SS productivity. Further, tree volume, instead of DBH, was selected in the regression model because this variable was evaluated as more stable and applicable [17,39,40]. Both the Berendt et al. [41] and Han et al. [42] studies found that the most commonly used variable is tree volume to estimate the productivity. Another important variable that influences productivity was significantly related to extraction distance since the load travel time accounted for over one quarter of the variation in cycle time [43,44]. Thus, extraction productivity increases with a decrease in extraction distance [45–47]. However, our follow-up data provide only an average total travelled distance of SS across the harvest unit, while the extraction distance was not collected. We found that productivity and TTD were strongly correlated. These results are consistent with many published studies, such as Matthews et al., [48], Kumazawa et al., [49], and Berg et al., [50]. These studies posited that the driving unloading speed increases with an increase in driving distance. During extraction operation, operator’s stress on long distances could lead to increase driving speed. In addition, extraction productivity may be influenced by the number of logs [46,50]. This pattern could be explained by the fact that productivity increases with a decrease in the number of load stops and driven distances. The operation of SS in South Korea, unlike the excavator extraction method (a.k.a., shovel logging), uses gravity. During the operation, the throwing and rolling activity, which is done to transfer logs from the felling area down to the landing area, may be associated with a large number of logs and long TTD. Thus, the productivity of SS may depend on the number of throwing stops and TTDs. Slope is shown to be a primary factor in productivity because it affects the accessibility of the harvesting machine [18,43,46,51,52]. However, this was also observed in Berendt et al. [41], Ghaffariyan et al. [53], and Walsh and Strandgard [54]. These studies found that the slope did not significantly impact loaded travel time during extraction activities. This study showed that slope had no statistically significant influence on SS productivity, since the time spent throwing and rolling logs may be considerably longer than the travel time. Further, in terms of extraction activity, there was no common pattern across Korea for selecting SS operations in steeply-sloped forests.

5. Conclusions In conclusion, this study’s findings showed that SS efficiency in extracting 2–4 m logs from stumps to roadside/landing areas varied across harvest units in mountain forests. The mean productivity was 6.03 1.90 m3/SMH, with a minimum of 2.44 m3/SMH and maximum of 9.85 m3/SMH. According to ± the cost deduction, the corresponding extraction costs were estimated to be from US $4.37 to 17.66/m3. Productivity was significantly impacted by DBH, tree volume, and TTD, whereas stand density and slope were observed to be non-significant. Assessing the productivity of forest operations is a challenging task due to rough and unstructured environmental conditions. In this study, an SS productivity prediction model was derived using follow-up data taken from 30 forest sites across South Korea. Two-thirds of the total dataset were used to build the model, while the remaining one-third of the data was used for validation. This model is expected to be used by harvest planners, forest managers, and decision makers to improve SS management and extraction operations. However, these results were limited to applications across South Korea. Further study is needed using a broader and more updated range of data with a tree volume greater than 0.6 m3 to develop an applicable prediction model that can be used in other countries.

Author Contributions: Conceptualization, E.L., S.H. and S.I.; Data collection and curation, E.L.; Formal analysis, E.L.; Funding acquisition, S.H. and S.I.; Investigation, S.H. and S.I.; Methodology, E.L., S.H. and S.I.; Writing-original draft, E.L.; Writing-review & editing, S.H. and S.I. Funding: This research was funded by the Korea Forest Service (Korea Forestry Promotion Institute), grant number 2013069B10-1919-AA03. Forests 2019, 10, 585 12 of 14

Conflicts of Interest: The authors declare no conflict of interest.

References

1. Senturk, N.; Ozturk, T.; Demir, M. Productivity and costs in the course of timber transportation with the Koller K300 cable system in Turkey. Biosyst. Eng. 2007, 42, 2107–2113. [CrossRef] 2. Ero˘glu,H.; Ozkaya, M.S.; Acar, H.H.; Karaman, A.; Yolasi˘gmaz,H.A. An investigation on roundwood extraction of Fagus orientalis lipsky, Abies nordmanniana (Stew.) Spach. and Picea orientalis (L.) Link. by Urus Mill forest skyline on snow. Afr. J. Biotechnol. 2009, 8, 1082–1089. 3. Gumus, S.; Ascar, H.H. Evaluation of consecutive skylines yarding and gravity skidding systems in primary forest transportation on steep terrain. J. Environ. Biol. 2010, 31, 213–218. [PubMed] 4. Melemez, K.; Tunay, M.; Emir, T. A comparison of productivity in five small-scale harvesting systems. Small Scale For. 2014, 13, 35–45. [CrossRef] 5. Stanki´c,I.; Poršinsky, T.; Tomaši´c,Ž.; Tonkovi´c,I.; Frnti´c,M. Productivity models for operational planning of timber forwarding in Croatia. Croat. J. For. Eng. 2012, 33, 61–78. 6. Visser, R.; Stampfer, K. Expanding ground-based harvesting onto steep terrain: A review. Croat. J. For. Eng. 2015, 36, 321–331. 7. Kulak, D.; Sta´nczykiewicz,A.; Szewczyk, G. Productivity and time consumption of timber extraction with a grapple in selected pine stands. Croat. J. For. Eng. 2017, 38, 55–63. 8. Proto, A.R.; Skoupy, A.; Macri, G.; Zimbalatti, G. Time consumption and productivity of a medium size mobile tower yarder in downhill and uphill configurations: A case study in Czech Republic. J. Agric. Eng. 2016, 47, 216–221. [CrossRef] 9. Spinelli, R.; Magagnotti, N.; Lombardini, C. Performance, capability and costs of small-scale cable yarding technology. Small Scale For. 2010, 9, 123–135. [CrossRef] 10. Lee, E.; Im, S.; Han, S.-K. Productivity and cost of a small-scale cable yarder in an uphill and downhill area: A case study in South Korea. For. Sci. Technol. 2018, 14, 16–22. [CrossRef] 11. Wang, J.; Long, C.; McNeel, J. Production and cost analysis of a feller-buncher and grapple skidder in central Appalachian hardwood forests. For. Prod. J. 2004, 54, 159–167. 12. Standards for Time Studies for the South African Forest Industry. Available online: http://www. forestproductivity.co.za/?page_id=678 (accessed on 17 August 2018). 13. Kellogg, L.D.; Bettinger, P. Thinning productivity and cost for a mechanized cut-to-length system in the north-west Pacific Coast region of the USA. J. For. Eng. 1994, 5, 43–53. 14. She, J.; Chung, W.; Kim, D. Discrete-event simulation of ground-based timber harvesting operations. Forests 2018, 9, 683. [CrossRef] 15. Olsen, E.; Kellogg, L. Comparison of time-study techniques for evaluating logging production. Trans. ASAE 1983, 26, 1665–1668. [CrossRef] 16. Strandgard, M.; Mitchell, R. Automated time study of using GPS and a vibration sensor. Croat. J. For. Eng. 2015, 36, 175–184. 17. Sirén, M.; Aaltio, H. Productivity and costs of thinning harvesters and -forwarders. Int. J. For. Eng. 2003, 14, 39–48. [CrossRef] 18. Spinelli, R.; Hartsough, B.R.; Magagnotti, N. Productivity Standards for Harvesters and Processors in Italy. For. Prod. J. 2010, 60, 226–235. [CrossRef] 19. Eriksson, M.; Lindroos, O. Productivity of harvesters and forwarders in CTL operations in northern Sweden based on large follow-up datasets. Int. J. For. Eng. 2014, 25, 179–200. [CrossRef] 20. Korea Forest Service. Available online: https://www.forest.go.kr (accessed on 10 October 2018). 21. Cho, K.-H.; Cho, M.-J.; Han, H.-S.; Han, S.-K.; Cha, D.-S. Harvesting cost of tree-length thinning in a Larix leptolepis stands. J. Korean For. Soc. 2015, 104, 221–229. [CrossRef] 22. Cho, M.; Cho, K.; Choi, B.; Cha, D. Yarding productivity of tree-length harvesting using small cable-yarder in steep slope, South Korea. For. Sci. Technol. 2018, 14, 132–137. [CrossRef] 23. Han, W.S.; Han, H.-S.; Kim, N.-H.; Cha, D.S.; Cho, K.S.; Min, D.H.; Kwon, K.C. Comparison of harvesting productivity and cost of cable yarding systems. J. Korean For. Soc. 2014, 103, 87–97. [CrossRef] 24. Kim, M.-K.; Baek, S.-A.; Cho, K.-H.; Jung, D.-H. Productivity, cost, and optimal forest road network density of tree-length yarding operations with tower yarder. J. Korean For. Soc. 2017, 106, 300–309. Forests 2019, 10, 585 13 of 14

25. Han, W.S.; Han, H.-S.; Kim, Y.-S.; Shin, M.Y. The yarding productivity and cost cable yarding operation by yarder attached on tractor: Using the Chuncheon tower-yarder. J. Korean For. Soc. 2008, 97, 641–649. 26. Spinelli, R.; Visser, R.; Thees, O.; Sauter, H.; Krajnc, N.; Riond, C.; Magagnotti, N. contract rates in the Alps: The effect of regional variability and technical constraints. Croat. J. For. Eng. 2015, 36, 195–203. 27. Ero˘glu,H.; Acar, H.H.; Ozkaya, M.S.; Tilki, F. Using plastic chutes for extracting small logs and short pieces of wood from forests in Artvin, Turkey. Build. Environ. 2007, 42, 3461–3465. [CrossRef] 28. Lee, E.; Li, Q.; Eu, S.; Im, S. Assessing the impacts of traditional mechanized logging operations on soil physical and hydrological properties in the Republic of Korea. Heliyon 2019. under review. 29. Miyata, E.S. Determining Fixed and Operating Costs of Logging Equipment; General Technical Report Nr. NC-55; Northern Central Forest Experiment Station, US Forest Service: Saint Paul, MN, USA, 1980; p. 16. 30. Kim, M.-K.; Park, S.-J. An analysis of the operational time and productivity in whole-tree and cut-to-length logging operation system. J. Korean For. Soc. 2012, 101, 344–355. 31. Lee, C.G.; Choi, Y.S.; Ju, Y.M.; Lee, S.Y.; Nam, Y.S.; Cho, L.H.; Kim, J.M.; Oh, K.C.; Cho, M.J.; Lee, E.J.; et al. A study on the productivity and cost analysis of the timber and logging residue in CTL system of excavator yarding for using the woody resources. Korean Soc. New Renew. Energy 2016, 12, 51–58. [CrossRef] 32. Kaakkurivaara, N.; Kaakkurivaara, T. Productivity and cost analysis of three timber harvesting extraction methods on steep terrain in Thailand. Croat. J. For. Eng. 2018, 39, 213–221. 33. Sowa, J.M.; Szewczyk, G. Time consumption of skidding in mature stands performed by winches powered by farm tractor. Croat. J. For. Eng. 2013, 34, 255–265. 34. Hiesl, P.; Benjamin, J.G. Applicability of international harvesting equipment productivity studies in Maine, USA: A literature review. Forests 2015, 4, 898–921. [CrossRef] 35. Lindroos, O.; Cavalli, R. Cable yarding productivity models: A systematic review over the period 2000–2011. Int. J. For. Eng. 2016, 27, 79–94. [CrossRef] 36. Session, J.; Boston, K. Optimization of road spacing for log length shovel logging on gentle terrain. Int. J. For. Eng. 2006, 17, 67–75. [CrossRef] 37. Landford, B.L.; Stokes, B.J. Comparison of two thinning systems. Part 2. Productivity and costs. For. Prod. J. 1996, 46, 47–53. 38. Jiroušek, R.; Klvaˇc,R.; Skoupý, A. Productivity and costs of the mechanized cut-to-length wood harvesting system in clear-felling operations. J. For. Sci. 2007, 53, 476–482. [CrossRef] 39. Visser, R.; Spinelli, R. Determining the shape of the productivity function for mechanised felling and felling-processing. J. For. Res. 2012, 17, 397–402. [CrossRef] 40. Mederski, P.S.; Bembenek, M.; Karaszewski, Z.; Łacka, A.; Szczepa´nska-Álvarez, A.; Rosi´nska,M. Estimating and modelling harvester productivity in pine stands of different ages, densities and thinning intensities. Croat. J. For. Eng. 2016, 37, 27–36. 41. Berendt, F.; Fortin, M.; Suchomel, C.; Schweier, J. Productivity, costs, and selected environmental impacts of remote-controlled mini forestry crawlers. Forests 2018, 9, 591. [CrossRef] 42. Han, H.; Chung, W.; She, J.; Anderson, N.; Wells, L. Productivity and costs of two beetle-kill salvage harvesting methods in northern Colorado. Forests 2018, 9, 572. [CrossRef] 43. Ghaffariyan, M.R.; Stampfer, K.; Sessions, J. Forwarding productivity in Southern Austria. Croat. J. For. Eng. 2007, 28, 169–174. 44. Strandgard, M.; Mitchell, R.; Acuna, M. Time consumption and productivity of a operating on a slope in a cut-to-length harvest system in a Pinus radiata D. Don pine plantation. J. For. Sci. 2017, 63, 324–330. 45. Spinelli, R.; Cacot, E.; Mihelic, M.; Nestorovski, L.; Mederski, P.; Tolosana, E. Techniques and productivity of coppice harvesting operations in Europe: A meta-analysis of available data. Ann. For. Sci. 2016, 73, 1125–1139. [CrossRef] 46. Nurminen, T.; Korpunen, H.; Uusitalo, J. Time consumption analysis of the mechanized cut-to-length harvesting system. Silva Fenn. 2006, 40, 335–363. [CrossRef] 47. Danilovi´c,M.; Stojni´c,D.; Kari´c,S.; Su´cevi´c,M. Transport of technical roundwood by forwarder and tractor assembly from poplar plantation. Nova Mech. Šumarstva 2014, 35, 11–21. 48. Matthews, G.; Dorn, L.; Hoyes, T.W.; Davies, D.R.; Glendon, A.I.; Taylor, R.G. Driver stress and performance on a driving simulator. Hum. Factors 1998, 40, 136–149. [CrossRef][PubMed] Forests 2019, 10, 585 14 of 14

49. Kumazawa, Y.; Fujita, M.; Yamasaki, A.; Koyama, K.; Ichihara, K.; Oka, M. Consideration about safety and efficient log transportation by forwarders. J. Jpn. For. Eng. Soc. 2011, 26, 181–186. 50. Berg, S.; Ersson, B.T.; Manner, J. Distance driven and driving speed when forwarding during final felling in central Sweden. J. For. Sci. 2019, 6, 183–194. [CrossRef] 51. Adebayo, A.B.; Han, H.-S.; Johnson, L. Productivity and cost of cut-to-length and whole-tree harvesting in a mixed-conifer stand. For. Prod. J. 2007, 57, 59–69. 52. Norihiro, J.; Ackerman, P.; Spong, B.D.; Längin, D. Productivity model for cut-to-length harvester operation in South African Eucalyptus pulpwood plantations. Croat. J. For. Eng. 2018, 39, 1–13. 53. Ghaffariyan, M.R.; Sessions, J.; Brown, M. Machine productivity and residual harvesting residues associated with a cut-to-length harvest system in southern Tasmania Southern Forests. J. For. Sci. 2012, 74, 229–235. 54. Walsh, D.; Strandgard, M. Productivity and cost of harvesting a stemwood biomass product from integrated cut-to-length harvest opeartions in Australian Pinus radiata plantations. Biomass Bioenergy 2014, 66, 93–102. [CrossRef]

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).