- Download PDF-
SEEFOR 9 (1): early view
Article ID: 166
Original scientific paper
Soil Measurements in the Context of Planning Harvesting Operations and Variable Climatic Conditions
Andreja Đuka1*, Tomislav Poršinsky1, Tibor Pentek1, Zdravko Pandur1, David Janeš1, Ivica Papa1
(1) University of Zagreb, Faculty of Forestry, Department of Forest Engineering, Svetošimunska 25, HR-10000 Zagreb, Croatia
* Correspondence: e-mail: email@example.com
Citation: ĐUKA A, PORŠINSKY T, PENTEK T, PANDUR Z, JANEŠ D, PAPA I 2018 Soil Measurements in the Context of Planning Harvesting Operations and Variable Climatic Conditions. South-east Eur for 9 (1): early view. DOI: https://doi.org/10.15177/seefor.18-04
Received: 6 Oct 2017; Revised: 3 Apr 2018; 16 Apr 2018; Accepted: 26 Apr 2018; Published online: 10 May 2018
Background and Purpose: Terrain trafficability and vehicle mobility during timber extraction are highly influenced by terrain characteristics, as well as by soil bearing capacity. Insufficient planning, too heavy vehicles, excessive traffic, poor soil knowledge, together with bad weather conditions, cause severe disturbance to forest soil. Damage to the forest ecosystem arises due to felling and timber extraction operations, regardless of the technical means used in this process. Traffic intensity plays an essential role in soil compaction because deformations can increase with the number of passes, which may lead to excessive soil disturbance.
Materials and Methods: The research was conducted in Zalesina, a hilly and pre-mountainous part of Gorski kotar, Croatia, in two management units (MUs) comprising of 569.64 ha of selective beech and fir forests. Physical and mechanical properties of forest soil, essential for vehicle mobility and terrain trafficability during timber extraction are based on measuring the following soil parameters: 1) soil penetration resistance, 2) soil shear strength, and 3) current soil moisture. Measurements were made during one calendar year, on five sampling plots (three in MU "Belevine" and two in MU "Kupjački vrh") chosen on the prevailing soil types in both MUs.
Results: The highest values of cone index and shear strength were recorded in July which was the warmest and driest month with only 7 rainy days. Current soil moisture was the lowest in July at φ=15.26% vol, while the lowest values of cone index were from October to January with precipitation at its climax from September to January with cumulative 1232 mm of rain and 846 cm of snow. Lower values of penetration resistance were in MU "Kupjački vrh" at 0.96 MPa, which is related to overall terrain structure of high karst and surface roughness, making the sampling of plots difficult. By comparing cone index values (CI15) and the mean values between 5 and 25 cm of depth (CI5-25), the variation coefficient indicates a lower variability of CI5-25 values on four out of five sampling plots. However, t-test showed no statistically significant difference between these soil parameters.
Conclusions: Data regarding soil bearing capacity in a map form as an additional layer to other maps of stand characteristics would ensure better opportunities for planners or operators to reduce and avoid rutting and soil compaction. Measured data on penetration resistance, soil shear strength and current soil moisture have not only given in a simple and a fast manner the in situ stand conditions, but they have also shown their dependence on climatic conditions.
Keywords: digital penetrometer, cone index, soil shear strength, precipitation, soil moisture content
The planning of harvesting operations usually starts with a thematically prepared maps and a pre-planned marking of stand borders and/or marking of trees. A planner besides marking cutting area borders, environmentally delicate areas and information on stand inventory needs to evaluate the best driving routes for vehicles on either temporary forest roads or trails , or on the existing primary and secondary forest traffic infrastructure using Geographic Information System (GIS). One of the most significant applications of LiDAR-derived terrain models has been in facilitating the mapping of areas of anticipated high moisture and potentially high susceptibility to soil damage by vehicles . Terrain trafficability and vehicle mobility during timber extraction are highly influenced by terrain characteristics, as well as by soil bearing capacity. Insufficient planning, too heavy vehicles, excessive traffic, poor soil knowledge, together with bad weather conditions, cause severe disturbance to forest soil. Damage to the forest ecosystem arises due to felling and timber extraction operations, regardless of the technical means used in this process  and traffic intensity plays an essential role in soil compaction because deformations can increase with the number of passes, which may lead to excessive soil disturbance [4, 5]. Therefore, one of the most important issues in ground-based timber extraction in forestry is to minimise the soil damage caused by heavy forestry machines .
Soil bearing capacity is the ability of the ground to withstand external forces (wheel or track movement), and is determined by compaction, rutting and surface disturbance under external load. In forestry, soil bearing capacity is defined as the maximum permissible tire pressure of the vehicle without damaging the soil , depending on the type and texture of the soil, humus layer and skeletal particles (permanent parameters of soil), and the only variable parameter – soil moisture . Soil bearing capacity is expressed in kPa (MPa, N·cm-2), but there is no prescribed measurement method, due to the difficulty of accurate and straightforward procedures caused by the natural condition of forest soil, rapid changes of its moisture and the influence of the vegetation roots . In terra-mechanics and wheel – soil interaction, Mohr-Coulomb's failure criterion is the most often applied one, where soil strength is quantified as the maximum permitted strain in the soil before its deformation (cracking, crushing or flowing).
Poršinsky  states that the share of sand and clay particles in the soil composition affects the boundary of the shear strength of soil. During vehicle movement on deformable ground, the transfer of force is relevant during load and deformation of the soil in the tangential direction, provided there is a friction between the ground and the tire. As tires of forest vehicles have mandatory treads, tangential stress is ensured in case of tread penetration into the ground. In this case, the most important feature of soil deformation is the strain because force can be achieved due to the shear strength of the soil. In general, it can be concluded that soil with the highest shear strength is over-consolidated clay or well-compacted sand, and soil with the least shear strength are poorly graded sands [10-14].
In forestry, soil bearing capacity is usually determined by the penetration of the cone into the ground and is defined as the ratio of force required to press the standardised cone, as well as varying ground resistance to penetration depending on its depth. The ground penetration curve contains data on the estimation of soil strength depending on the depth of the cone penetration, caused by the horizon condition of certain soil types. In practical application, a standardised value of the cone penetration measurement  is used in a 15 cm ground depth called the cone index (CI). The same standard recommends the use of a cone index when assessing terrain trafficability and the vehicles’ mobility as a parameter used in modelling vehicle traction performance. However, an altered definition of the cone index can be found in literature, where the cone index is the average value of the penetration resistance measurement in the range from 0 to 15 cm of soil depth . Horvat  recommends repeated measurements on forest soils due to their inhomogeneous structure, and considers penetrometer measurements as suitable for fast comparative mensuration, mainly when performed with the same penetrometer.
The first forest ground bearing capacity classification based on cone index (CI), shear strength and deformation modulus (E value), as each of these can be measured at once in the field using portable equipment, was made by the EcoWood project, which gave special attention to environmentally efficient timber harvesting operations on sensitive soils . Poršinsky  and Poršinsky et al. , based on measurements of multiple parameters of forest soil bearing capacity in Croatia, expressed doubts about the range of the values representing high ground bearing capacity (H-GBC) at >500 kPa for cone index value; >60 kPa of shear strength. In accordance with the above-mentioned suspicion, Poršinsky et al.  state that Bygdén's (2012) guideline for ground bearing capacity in forestry is more suitable for the conditions of Croatian forestry, where CI values range between 700 and 1000 kPa, which is categorised as medium hard soil (may soften after rain), and CI above 1000 kPa, which represents hard soil with a good bearing capacity.
Saarilahti  recommends that, during soil defrosting, constraints in the production system should be identified due to increased moisture and reduced soil bearing capacity, so major soil damage is reduced. Shoop et al.  state that the snow cover reduces soil damage, but also the systems’ productivity. As the basic features of the snow cover, the authors state its thickness, temperature and density. The vehicle's mobility stops at a snow thickness of 30% of the wheel diameter, while chains on wheels enlarge mobility up to 50% of wheel diameter .
Access to soil data on the European level is possible via European Soil Data Centre (ESDAC). ESDAC gives topsoil physical properties for Europe (based on LUCAS topsoil data) for ca 20,000 samples across EU. Some layers for soil properties have been created based on the data from the European Soil Database in combination with data from the Harmonized World Soil Database (HWSD) and Soil-Terrain Database (SOTER) in 10×10 km resolution. The available layers include total available water content, depth available to roots, clay content, silt content, sandy content, organic carbon, bulk density, and coarse fragments . On the national level, an existing Pedological Base Map of Croatia in 1:300,000 scale  is still in use.
Talbot et al.  state that the development of remote and proximal sensing technology and techniques will provide a previously inconceivable amount of data. Authors continue that the machine-mounted sensors that unceasingly collect vast amounts of data will especially give the forest operation researcher a large and continually increasing basis from which to extract useful information.
These data can, with the application of sensible analytical approaches, provide significant opportunities for decision support, as well as for operation monitoring and evaluation. Besides the advantage of the availability and use of high-resolution data, prominent future development of GIS applications in forest operations could be an entirely integrated 4D spatial dataset and Virtual Reality (VR) systems .
This article deals with measuring soil bearing capacity through measurements of cone index, shear strength and soil moisture content to plan harvesting operations in the right period and gain more knowledge on soil properties specific for operational forestry purposes.
MATERIALS AND METHODS
The research was conducted in Zalesina, a hilly and pre-mountainous part of Gorski kotar, Croatia, in two management units comprising of 569.64 ha of selective beech and fir forests. The average growing stock is 446 m3·ha-1, with an annual increment of 6.83 m3·ha-1 and a total allowable cut in the period from 2010 until 2019 at 30,565 m3.
According to Köppen climate classification, the research area belongs to the Cfsbx climatic type .
The essential features of the Cfsbx climatic type are cool summers with the mean temperature of the hottest month below 22oC, while the average monthly cool summer temperature ranges from -3oC to +18oC. In the research area chilled air masses are retained due to its geographic position in a wet and closed plateau, leading to a low average annual temperature of just 6.8oC.
Zalesina is abundant with precipitation, with annual average values up to 2074 mm in 136 days of precipitation and with two yearly minimums, one at the end of winter – in March with average precipitation of 140 mm, and the other one in summer – in July with average precipitation of 143 mm . Climate elements affect harvesting operations, which may be best reflected by the data on the highest daily precipitation. Thus, in one day up to 154 mm of rain may fall (data from September 1951) in the form of strong showers, which is the same amount of precipitation in June. Therefore, soil protection from erosion and torrents is necessary. A significant share of annual precipitation is significant for late autumn and winter, which consists of 98 days with snow, while the snowfall lasts on average 188 days (from the end of October to the end of June). There is also a significant proportion of snowfall in annual precipitation, up to 26%, so that in spring the melting of snow significantly soaks the forest soil.
In the two adjacent management units (MUs) "Belevine" and "Kupjački vrh", both with natural fir and beech forests, managed in selective cuts, variations in terrain conditions are evident. The relief of MU "Belevine" is mildly undulating and fan-like striated with many beds of water streams with their steep sides affecting timber extraction routes. MU does not contain stoniness, rockiness or similar ground obstacles. Terrain slope on 55.73% of the MU is between 0–20% . Terrain conditions and harvesting system in use resulted in a skid trail network for timber extraction by a skidder. In the MU "Kupjački vrh", karst phenomena and cliffs are common, with a centrally located peak at 997 m and with hillsides descending in all directions. More than 50% of the MU has terrain slopes higher than 34%, and 65.87% of MU is categorised as having uneven to very rough terrain regarding stoniness and rockiness . This all resulted in a built skid road network without which skidders’ mobility in this MU would be unattainable .
Physical and mechanical properties of forest soil essential for vehicle mobility and terrain trafficability during timber extraction are based on measuring the following soil parameters: 1) soil penetration resistance, 2) soil shear strength, and 3) current soil moisture. Measurements were made during one calendar year, on five sampling plots (three in MU "Belevine" and two in MU "Kupjački vrh") chosen on the prevailing soil types in both MUs (Figure 1):
- BEL 1 – dystric colluvial brown soil
- BEL 2 – dystric brow soil (typically brown podzolic)
- BEL 3 – deep dystric brown soil
- KUP 1 – shallow rendzina
- KUP 2 – shallow and medium deep rendzina
A digital penetrometer Eijkelkamp Penetrologger was used to measure penetration with a cone surface range of 2 cm2 and the angle of cone at 30° according to the standard ASAE S.313.3 at a penetration speed of 2 s·cm-1. The measurements were performed on a monthly basis in 15 repetitions per plot at a depth of maximal 80 cm. Soil bearing capacity was defined by cone index as a derivative of soil penetration resistance at 15 cm depth .
For measuring soil shear strength, Eijkelkamp Field inspection vane tester was used with 20×40 mm wing dimension, a measuring range of 0 to 130 kP, and a reading accuracy of 2 kPa. The measurements were performed on a monthly basis in 15 repetitions per plot at a depth of 15 cm.
The current soil moisture was measured by a Fieldscount TDR 300 soil moisture meter at a depth of 15 cm. Measurements were made in repetitions of 15 on all five surface plots, with each of the 15 measurements being the average of 10 measurements so that 150 current moisture data were measured on each plot.
From the State Meteorological and Hydrological Service data on average daily temperatures (°C); average daily wind power (m·s-1); daily height of total snow (cm); daily precipitation (mm); and daily insolation (hours) were taken for the analysis and better understanding of the soil bearing capacity dynamics during the year.
According to the national hydro-meteorological data, there were no discrepancies in temperature for autumn and winter in the measuring period, but the amount of precipitation varied. Spring was characterised as warm to very warm (higher values by 0.8-1.5°C) and precipitation was normal (101%) to dry (67%). Summer was extremely warm (temperatures higher by 1.7°C) with normal amount of precipitation at 101%.
Soil Moisture Content
By a repeated ANOVA analysis made in software Statistica 13.1, a statistically significant difference between sampling plots and the period of measurements were determined (Table 1).
Tukey Post Hoc test (Table 2) showed statistically significant difference between some sampling plots with MS=26.747 and df=70.
Figure 1 shows difference in values regarding the time of measurements (month) and location (sampling plot). The highest average monthly moisture values were recorded in October, and the lowest values in July. There were no field measurements in December due to high snow which prevented approach to the sampling plots.
Figure 3 comprises the climatic data (precipitation amount and snow) together with soil moisture data during twelve months of measurements. The peak of soil moisture content was from September to November on average φ=39.10%, which correlates to the high amount of precipitation above 272 mm.
Precipitation showed more considerable influence on soil moisture content than snow, whose climax was in January and which remained present in the forest stand till April.
By repeated ANOVA analysis of soil penetration resistance, a statistically significant difference between sampling plots and the period of measurements was determined (Table 3).
Tukey Post Hoc test showed statistically significant differences between all five sampling plots at MS=0.20571 and df=2170.
Figure 4 represents average values of soil penetration resistance per sampling plot in a given period of measurements. The peak of measurements was in July at 3 MPa, and the lowest value was measured at KUP 2 sampling plot in October – 0.62 MPa. The lowest average measurements for all five sampling plots were recorded in January at 0.86 KPa. According to Bygdén’s ground bearing capacity classification, all sampling plots (except KUP 2 in October) regardless the month, can be classified as medium hard to hard soils.
Average values of soil penetration resistance per plot and the period of measurements but depending of soil depth are given in Figure 5. By increasing soil depth, penetration resistance increases (Figure 5), and correlation coefficients show strong relationship (0.82-0.99) between these two variables. The highest average values were recorded in July and August from 1.62 MPa to 3 MPa from the surface up to 30 cm depth. The lowest average values at the soil surface were in February and March, 0.35 MPa and 0.30 MPa respectively. At 15 cm and 30 cm depth minimal values were recorded in October and January, 0.91 MPa and 1.14 MPa respectively. The period of measurements and soil depth show lower influence on soil penetration resistance for sampling plots in MU "Kupjački vrh". On KUP 1 sampling plot the difference in soil penetration resistance between 10 cm and 30 cm depth ranges up to 0.45 Ma, while on KUP 2 sampling plot the difference is smaller and ranges up to 0.25 MPa. On both sampling plots, the highest jump in penetration resistance is at first 10 cm of soil depth. On sampling plots in MU "Belevine" in a month of measurements there has been higher influence on penetration resistance ranging up to 1.0 MPa.
By comparing the values of cone index at a depth of 15 cm (CI15) and mean values between 5 and 25 cm of depth (CI5–25), the coefficient of variation (CV) indicates a lower variability of CI5–25 values at four sampling plots from 14.56% to 22.84%. At KUP 2 sampling plot variation coefficients of both cone indexes were the same.
T-test, however, showed no statistically significant difference between the mean values of the cone index at a depth of 15 cm and 5 to 25 cm.
By repeated ANOVA analysis of shear strength a statistically significant difference between sites and the period of measurements was recorded (Table 4). Tukey Post Hoc test showed statistically significant difference between some sampling plots (Figure 6) at MS=234.07 and df=70, while there were no differences between sampling plots BEL 1 and BEL 2, nor between BEL 3 and KUP 1.
The highest values of soil shear strength at 15 cm depth were recorded in July at τ=130 kPa, while the lowest values were recorded in May at τmax=61 kPa. If monthly mean values are compared to all five sampling plots, data deviation is visible, except in the summer months (July and August) when variation is at its minimum. The lowest mean values of the soil shear strength were measured on sampling plot KUP 1 (τ=62 kPa), and the absolute minimum was measured in October (τ=32 kPa).
When comparing cone index and shear strength to soil moisture (Figure 7) on all sampling plots during one-year period and in 15 cm depth, a consistency is evident. The years’ peak of cone index (CI=2109 kPa) and shear strength values (τ=129 kPa) in July was followed by a soil moisture minimum (φ=15.26%), after which soil moisture content started to rise steadily, only to reach its high point in October at φ=41.21%. In this period both the cone index and shear strength were descending and in October reached 891 kPa and 87 kPa, respectively. After October soil moisture content plunged steeply only to be φ=18.05% in April. Cone index and shear strength fluctuated in this period from 997±120 kPa and 100±7 kPa, respectively.
The correlation between soil moisture, precipitation and temperature data has been discussed widely [30, 31] and projects such as Global Soil Wetness Project (GSWP) reproduced the seasonal cycle of soil moisture . To retrieve soil moisture, the effects of vegetation, surface roughness, and heterogeneous land cover must be taken into account , and in the territory of the Republic of Croatia, these data specific to forestry operations are still missing. There have been several studies on soil penetration resistance having water content or moisture content (MC) as an input variable [34, 35]. The highest values of cone index and shear strength were recorded in July which was the warmest (above average) and driest (within average) month with only 7 rainy days. Current soil moisture was the lowest in July at φ=15.26% vol. Climate characteristics in Croatia have shown a trend of decreasing precipitation and increasing temperatures during all seasons in the 20th century . The lowest values of cone index were recorded from October to January with precipitation at its climax from September to January with cumulative 1232 mm of rain and 846 cm of snow. Even though soil moisture is a versatile parameter, data recorded by previous research and models combining cone index and soil moisture by Freitag  and Hinze’s  are non-linear models which all give dependence between these parameters. Data analysis of penetration resistance showed statistically significant differences between five sampling plots, with their average at 1.12 MPa, which according to Bygdén  are hard soils with good bearing capacity. Lower values of penetration resistance were recorded in MU "Kupjački vrh" at 0.96 MPa, which is related to the overall terrain structure of high karst and surface roughness  in which the placement of sampling plots was difficult. Two chosen areas of KUP 1 and KUP 2 sampling plots were based on a search of deep soil exceptions so that measurements with digital penetrometer, vane tester and soil moisture meter could be achieved.
By comparing cone index values (CI15) and the mean values between 5 and 25 cm of depth (CI5-25), the variation coefficient indicates a lower variability of CI5-25 values on four out of five sampling plots which correlate with a research by Pandur . However, t-test showed no statistically significant difference between these soil parameters.
Even though soil penetration resistance was measured up to 80 cm of depth, the majority of changes in penetration resistance were in top 10 to 15 cm of soil after which curves stabilised and showed little to no difference depending on soil depth. Measurements up to 30 cm of soil depth are sufficient for operational forestry purposes [40-42] even though Jansson and Johansson  found that traffic-increased bulk density of a silt loam podzol was down to 40–50 cm of soil depth.
Vane tester in a simple way measures the unconfined shear strength of soil – soil cohesion , and gives the maximal strength that must be used to provoke a free movement of the soil. Formerly established dependence between penetration resistance and soil shear strength introduced by Meek  has shown its influence during field measurements (Figure 7). The lowest values of soil shear strength were recorded from April to June on average as τ=66 kPa, which according to Ward et al.  places these soils to high ground capacity category (H-GBC). Due to expressed doubts regarding the applicability of the mentioned categorisation in Croatian forestry, at least 80 kPa should be considered as medium hard soils, and values above 100 kPa as hard soils.
Understanding terrain factors is critical for planning forest operations. Various stand conditions affect the type of machines, harvesting system and the activity itself regarding eco-efficiency, productivity, and revenue. Data regarding soil bearing capacity in a map form as an additional layer to other maps of stand characteristics would ensure better opportunities for planners or operators to reduce and avoid rutting and soil compaction. The measured data on penetration resistance, soil shear strength and current soil moisture have not only given in a simple and fast manner the in situ stand conditions, but have also shown their dependence on climatic conditions. In the absence of these data on the national level, particularly intended for forestry operations in which top 300 mm of soil are crucial, future research challenges should be addressed regarding: 1) the development of unified and straightforward methods for describing forest soil bearing capacity, 2) linking forest soil parameters to long-term monitoring of climatic elements, 3) connecting forest soil-bearing capacity to eco-indicators (occurrence of hygrophytes in forest phytocenoses), and 4) quantifying forest soil bearing capacity classification suitable for a wide range of conditions on the national level. Gaining knowledge on versatile soil bearing capacity throughout the calendar year would support route planning for primary and secondary traffic infrastructure network, as well as defining sensitive parts of logging areas and landing sites, and alarming machine operators to equip vehicles with wider tyres, chains or bogie tracks in specific conditions.
The study was carried out within the framework of "Optimisation of harvesting systems and forest traffic infrastructure on strategical and tactical level of planning" financed by the Croatian Ministry of Agriculture and non-market forest values.
- BYGDÉN G 2012 GIS for Operative Support. In: Okia CA (ed) Global Perspectives on Sustainable Forest Management. IntechOpen Limited, London, UK, pp 217-222. DOI: https://doi.org/10.5772/33230
- TALBOT B, PIERZCHAŁA M, ASTRUP R 2017 Applications of remote and proximal sensing for improved precision in forest operations. Croat J For Eng 38 (2): 327-336
- CUDZIK A, BRENNENSTHUL M, BIAŁCZYK W, CZARNECKI J 2017 Damage to Soil and Residual Trees Caused by Different Logging Systems Applied to Late Thinning. Croat J For Eng 38 (1): 83-95
- OZTURK T 2016 The effects on soil compaction of timber skidding by tractor on skid road in plantation forest in northern Turkey. Sumar list 140 (9-10): 485-490
- SOLGI A, NAGHDI R, TSIORAS PA, ILSTEDT U, SALEHI A, NIKOOY M 2017 Combined Effects of Skidding Direction, Skid Trail Slope and Traffic Frequency on Soil Disturbance in North Mountainous Forest of Iran. Croat J For Eng 38(1): 97-106.
- CAMBI M, GRIGOLATO S, NERI F, PICCHIO R, MARCHI E 2016 Effects of forwarder operation on soil physical characteristics: a case study in the Italian alps. Croat J For Eng 37 (2): 233-239
- SAARILAHTI M 2002 Modelling of the wheel and soil. Soil interaction model. Development of a protocol for ecoefficient wood harvesting on sensitive sites (Ecowood), Appendix report no 8, 37 p
- PORŠINSKY T, STANKIĆ I, BOSNER A 2011 Ecoefficient timber forwarding based on nominal ground pressure analysis. Croat J For Eng 31 (1): 345-356
- LÖFFLER HJ 1984 Terrain classification for forestry. Report TIM/EFC/WP.1/R.51, 24 August 1984, EU Timber Committee and FAO-ILO, 55 p
- PORŠINSKY T 2005 Efficiency and environmental evaluation of Timberjack 1710 forwarder on roundwood extraction from Croatian lowland forests. PhD thesis, University of Zagreb, Faculty of Forestry, Zagreb, Croatia, 170 p
- BALADI GY 1987 Terrain evaluation for off-road mobility. J Terramechanics 24 (2): 127-140. DOI: https://doi.org/10.1016/0022-4898(87)90003-6
- HORVAT D 1996 Estimate of forest vehicle tractive performance based on cone penetrometer soil characteristic. Protection and utilisation of the Croatian forest resources from 1846 to 1996, Croatian Forestry Society, Zagreb, Croatia, pp 253-264
- SAARILAHTI M 2002 Soil interaction model. Development of a protocol for ecoefficient wood harvesting on sensitive sites (Ecowood), Project deriverable D2, pp 1-87
- KVASNIČKA P, KOVAČEVIĆ ZELIĆ B, DOMITROVIĆ D, VEINOVIĆ Ž, STRGAR H 2007 Soil mechanics 1. Faculty of Mining, Geology and Petroleum Engineering, Carnet, Zagreb, Croatia, 216 p
- ASAE 1999 Procedures for using and reporting data obtained with the soil cone penetrometer (ASAE EP 542 FEB99). ASAE Standards 2000: 986-989.
- SAARILAHTI M 2002 Soil interaction model. Development of a protocol for ecoefficient wood harvesting on sensitive sites (Ecowood), 87 p
- WARD SM, OWENDE PMO 2003 Development of a protocol for eco-efficient wood harvesting on sensitive sites. In: Proceedings of the 2nd International Scientific Conference "Forest and Wood-Processing Technology vs. Environment – Fortechenvi Brno 2003", Brno, Czech Republic, 26-30 May 2003. Mendel University of Agriculture and Forestry Brno & IUFRO WG 3.11.00, Brno, Czech Republic, 473-482
- PORŠINSKY T, SRAKA M, STANKIĆ I 2006 Comparison of Two Approaches to Soil Strength Classifications. Croat J For Eng 27 (1): 17-26
- SAARILAHTI M 2002 Dynamic terrain classification – Modelling of the seasonal variation of the trafficability on forest sites. Soil interaction model. Development of a protocol for ecoefficient wood harvesting on sensitive sites (Ecowood), Appendix report no. 1, 22 p
- SHOOP S, YOUNG B, ALGER R, DAVIS J 1994 Effect of test method in winter traction measurements. J Terramechanics 31 (3): 156-162. DOI: https://doi.org/10.1016/0022-4898(94)90012-4
- NUTTAL CJ, MCGOWAN RP 1962 Scale models of vehicles in soils and snows. In: Proceedings of 1st ISTVS Conference, Torino, Italy, pp 656-667
- PANAGOS P, VAN LIEDEKERKE M, JONES A, MONTANARELLA L 2012 European Soil Data Centre: Response to European policy support and public data requirements. Land Use Policy 29 (2): 329-338. DOI: https://doi.org/10.1016/j.landusepol.2011.07.003
- BOGUNOVIĆ M, VIDAČEK Ž, RACZ Z, HUSNJAK S, SRAKA M, 1997 The practical aspects of soil suitability map of Croatia. Agronomski glasnik: Glasilo Hrvatskog agronomskog društva 59 (5-6): 363-399
- GRIGOLATO S, MOLOGNI O, CAVALLI R 2017 GIS applications in forest operations and road network planning: an overview over the last two decades. Croat J For Eng 38 (2): 175-186
- SELETKOVIĆ Z, KATUŠIN Z 1992 Croatian Climate. In: Rauš Đ (ed) Forests in Croatia, Faculty of Forestry and Croatian forests Ltd., Zagreb, Croatia, pp 13-18
- ANON 2009 Management program of forests with a specified use, forest training and research center Zalesina, MU Belevine 2010 – 2019. University of Zagreb, Faculty of Forestry, Zagreb, Croatia, 178 p
- ĐUKA A, PORŠINSKY T, VUSIĆ D 2015 DTM models to enhance planning of timber harvesting. Bulletin of the Faculty of Forestry - University of Belgrade 2015: 35-44. DOI: https://doi.org/10.2298/GSF15S1035D
- ĐUKA A, VUSIĆ D, PORŠINSKY T 2017 Terrain roughness evaluation for timber extraction by cable skidder. Gozdarski vestnik 75 (1): 21-35.
- ĐUKA A, GRIGOLATO S, PAPA I, PENTEK T, PORŠINSKY T 2017 Assessment of timber extraction distance and skid road network in steep karst terrain. iForest 10 (6): 886-894. DOI: https://doi.org/10.3832/ifor2471-010
- HENDERSON-SELLERS A 1996 Soil moisture: A critical focus for global change studies. Global Planet Change 13 (1-4): 3-9. DOI: https://doi.org/10.1016/0921-8181(95)00034-8
- DAI A, TRENBERTH KE, QIAN T 2004 A global dataset of Palmer Drought Severity Index for 1870-2002: Relationship with soil moisture and effects of surface warming. J Hydrometeorol 5 (6): 1117-1130. DOI: https://doi.org/10.1175/JHM-386.1
- ENTIN JK, ROBOCK, A, VINNIKOV KY, ZABELIN V, LIU S, NAMKHAI A, ADYASUREN T 1999 Evaluation of global soil wetness project soil moisture simulations. J Meteorol Soc Jpn 77 (1B): 183-198. DOI: https://doi.org/10.2151/jmsj1965.77.1B_183
- WAGNER W, SCIPAL K, PATHE C, GERTEN D, LUCHT W, RUDOLF B 2003 Evaluation of the agreement between the first global remotely sensed soil moisture data with model and precipitation data. J Geophys Res-Atmos 108 (D19): 1-17. DOI: https://doi.org/10.1029/2003JD003663
- AYERS PD, PERUMPRAL JV 1982 Moisture and density effect on cone index. Transactions of the ASAE 1982 (5): 1169-1972. DOI: https://doi.org/10.13031/2013.33691
- SULLIVAN PM 1999 Impact of weather data accuracy on NRMMII vehicle performance predictions. In: Proceedings of the 13th International Conference of the ISVTS, Munich, Germany, 14-17 September 1999, pp 55-62
- ZANINOVIĆ K, GAJIĆ-ČAPKA M 2000 Changes in components of the water balance in the Croatian lowlands. Theor Appl Climatol 65 (1-2): 111-117. DOI: https://doi.org/10.1007/s007040050008
- FREITAG DR 1985 Soil dynamics as related to traction and transport systems. In: Proceedings of International Conference on Soil Dynamics, Auburn, Alabama, vol. 4: 605-629
- HINZE D 1990 The influence of seasonal moisture changes on soil strength. In: Proceedings of the 10th International Conference of the ISVTS, Kobe, Japan, 20-24 August 1990, vol. 1: 107-116
- PANDUR Z 2013 Application of commercial monitoring system for the purpose of study forwarder extracting features. PhD thesis, University of Zagreb, Faculty of Forestry, Zagreb, Croatia, 312 p
- WILLIAMSON JR, NEILSEN WA 2000 The influence of forest site on rate and extent of soil compaction and profile disturbance of skid trails during ground-based harvesting. Can J For Res 30 (8): 1196-1205. DOI: https://doi.org/10.1139/x00-041
- LANDSBERG JD, MILLER RE, ANDERSON HW, TEPP JS 2003 Bulk density and soil resistance to penetration as affected by commercial thinning in north-eastern Washington. USDA For. Serv. Res. Pap. PNW-RP-551.
- HAN S-K, HAN H-S, PAGE-DUMROESE DS, JOHNSON LR 2009 Soil compaction associated with cut-to-length and whole-tree harvesting of a coniferous forest. Can J For Res 39 (5): 976-989. DOI: https://doi.org/10.1139/X09-027
- JANSSON K, JOHANSSON J 1998 Soil changes after traffic with a tracked and a wheeled forest machine: a case study on a silt loam in Sweden. Forestry 71 (1): 57-66. DOI: https://doi.org/10.1093/forestry/71.1.57
- MEEK P 1996 Effects of skidder traffic on two types of forest soils. Forest Engineering Research Institute of Canada, Pointe-Claire, QC. Technical Report No. TR 117, 12 p
© 2018 by the Croatian Forest Research Institute. This is an Open Access paper distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0).