The Effects of Soil Moisture on Harvesting Operations in Populus spp. Plantations: Speciﬁc Focus on Costs, Energy Balance and GHG Emissions

: Background: Poplar tree plantations for wood production are part of a worldwide growing trend, especially in moist soil sites. Harvesting operations in moist sites such as poplar plantations require more study for detailed and increased knowledge on environmental and economic aspects and issues. Methods: In this study, the effects of soil moisture content (dry vs. moist) on productivity, cost, and emissions of greenhouse gases (GHG) caused by operations of different harvesting systems (chainsaw-skidder and harvester-forwarder) were evaluated in three poplar plantations (two in Italy and one in Iran). Results: The productivity (m 3 h − 1 ) of both systems in the dry sites were signiﬁcantly higher (20% to 30%) than those in the moist sites. Production costs ( € m − 3 ) and GHG emissions (g m − 3 ) of both systems in the dry sites were also signiﬁcantly lower than those in the moist sites. The productivity of the harvester-forwarder system was about four times higher, and its production cost was 25% to 30% lower than that of the chainsaw-skidder system, but the calculated GHG emissions by harvester-forwarder system was 50–60% higher than by the chainsaw-skidder system. Conclusions: Logging operations are to be avoided where there are conditions of high soil moisture content (>20%). The result will be higher cost-effectiveness and a reduction in the emission of pollutants.


Introduction
Poplar planting has occurred around the world for a very long time. The plantations in Iran and Italy provide an important source of wood supply. At present, in both countries there are over 100,000 ha of these monospecific plantations (50,000 ha in Iran and 66,000 ha in Italy) and they mainly consist of Populus deltoides and P. euramericana [1][2][3]. Although poplar plantations cannot be currently considered among the main sources of wood in both countries, their importance is rapidly increasing [4,5]. Poplar wood shows interesting features, such as uniform mechanical properties and a high percentage of juvenile wood. These make it possible to obtain several products from plantations, i.e., building and veneering material, paper pulp and wood chips for bioenergy [2,6,7]. Moreover, the poplars in both the Italian and Iranian conditions can reach a growth rate of approximately 10-30 m 3 ha −1 year −1 , which is substantially higher than local tree species [8][9][10].
This study was carried out in three different geographical areas, two in Italy and one in Iran, investigating three different harvesting operations in poplar plantations.
The Iranian poplar plantation (IR) is located in the coastal area of the Caspian Sea in the Guilan province in northern Iran. The total area of the plantation is 60 ha, situated in flat terrain at an altitude range from 0 m to 20 m a.s.l. The average annual rainfall is from 1260 to 1340 mm and most of the precipitation occurs between the months of September and December. The average annual temperature is +15 • C, with the minimum during the winter at a few degrees below 0 • C, and the maximum at +25 • during the summer. The soil type is clay loam with poor drainage. This plantation was divided into two areas of 30 ha each. The harvesting operation in 30 ha for the dry site was performed in the first half of September 2019 before the rainfall. Harvesting operations were carried out in the second half of September 2019 after rainfall in the 30 ha moist site. The soil moisture content was 14.4% in the dry site and 34.6% in the moist site. Trees were felled by chainsaw, and whole trees were extracted to roadside landings by wheeled skidder Timberjack 450C. Finally, processing operations were motor-manually performed at the landing site.
The first Italian poplar plantation (IT1) was located in the Lazio region in central Italy. The total area of the plantation is 20 ha and it is situated in flat terrain and with an altitude range from 90 m to 110 m a.s.l. The average annual rainfall is from 830 to 900 mm and most of the precipitation occurs from October to December. The average annual temperature is +14.9 • C with a minimum during the winter at 2.5 • C, and a maximum at +30.7 • C during the summer. The soil type is clay loam with a low level of organic matter, nitrogen and phosphorus and with poor drainage. This plantation was divided into two areas of 10 ha each. The harvesting operations in 10 ha on the dry site were performed in the second half of June 2018. Harvesting operations were carried out in the first half of April 2018 after rainfall in the 10 ha moist site. The soil moisture content was 12.1% in the dry site and 36.8% in the moist site. As in the Iranian site, trees were motor-manually felled by chainsaw, and whole tree extraction was carried out by a wheeled skidder Timberjack 450C. In this case too, motor-manual processing with a chainsaw was carried out at the landing site.
The second Italian poplar plantation (IT2) was located in the Veneto region in North Italy. The total area of the plantation is 20 ha and situated in flat terrain with an altitude range from 10 m to 30 m a.s.l. The average annual rainfall is from 730 to 850 mm and most of the precipitation occurs from September to November. The average annual temperature is +13.6 • C with a minimum during winter of 0.1 • C, and a maximum of +29.2 • C during the summer. The soil type is clay loam texture with a low level of organic matter and with poor drainage. This plantation is divided into two areas of 10 ha each. The harvesting operations in the 10 ha dry site were performed in the first half of July 2018. Harvesting operations were carried out in the second half of September 2018 after rainfall in the 10 ha moist site. The soil moisture content was 15.0% in the dry site and 35.7% in the moist site. Trees were mechanically felled and processed by a harvester and cut to length through an extraction system by a forwarder. Technical characteristics of the machinery used in the various harvesting sites are given in Table 1. Average dendrometric characteristics of the three plantations are shown in Table 2. A preliminary analysis for dendrometric characteristics of the three different stands was done by one-way ANOVA to check for differences among the average values of the three plantations. There were no significant differences of dendrometric characteristics between the Italian sites (IT1 and IT2). However, density, basal area and standing volume of trees in the Iranian site were higher than in the Italian sites.

Data Collection and Analysis
Dendrometric data were obtained through systematic plot sampling. Grid dimension was 150 m × 150 m, the area of each circular plot was 1256 m 2 (20 m radius), and in total 20 plots in each area were established. Diameter at breast height (dbh) and height of tree species were measured by caliper and clinometer, respectively, in each plot. The volume of winched logs was calculated by Huber's formula (V = A m × L), where V is log volume (m 3 ), A m is the middle point cross-sectional area of log (m 2 ), and L is the length of log (m).
Soil samples were collected with a steel ring (inside diameter 5 cm, length 10 cm) and immediately put in hermetic plastic bags and labeled. The wet weight of all samples was measured before transfer to the laboratory (on the same sampling day). In the laboratory, soil samples were dried in an oven at 105 • C for 24 h until reaching a constant mass to determine the soil moisture content. A time-motion study was carried out to evaluate working productivity. Each working cycle was stop watched individually, separating productive time from delay time [33]. Calculated delay factor represents the quotient of delay time over net cycle time. Productivity was evaluated both on delay-free time and on actual total time, inclusive of all delays. Inclusion of delays was not capped on the basis of a maximum event duration. Scheduled Machine Hours (SMH) include all the time the machine is scheduled to work, whereas Productive Machine Hours (PMH) represent the time during which the machine actually performs work, excluding the time lost to both mechanical and non-mechanical delays.
The working group had 10 to 15 years of work experience with the machines and they were able to service and repair them.
The working cycles are reported, for the three areas, in Tables 3 and 4. Continuous time was recorded to the nearest second with a chronometer. The cycle times of the machines were divided into time elements (process steps) that were considered characteristic of this work. starts when the chainsaw operator moves from the last felled tree to the next to be felled and ends when the team cleans the tree stump before the felling starts when the harvester wheels start moving from one standing point and ends when they stop at the next standing point Felling (F) starts when the chainsaw operator turns on the chainsaw and performs the cut and ends with the fall of the tree starts when the harvester head grips the stem and ends when the tree falls onto the ground Processing (P) starts when the chainsaw operator cuts the first branch and ends when he finishes the cross cutting of the tree starts when the tree stem starts moving through the harvester head and ends when the harvester wheels start moving Table 4. Description of extraction cycle elements of harvesting in the three yards.

Time Elements IR & IT1 by Skidder (2 operators) IT2 by Forwarder (1 Operators)
Travel unloaded (TUL) begins when the skidder leaves the roadside landing area and ends when the skidder arrives at a suitable position (nearest distance from the logs) on the skid trail begins when the forwarder leaves the roadside landing area and ends when the forwarder arrives at the first suitable position (nearest distance from the first logs) Bunching-Loading (B) begins when the skidder driver releases the cable and ends when the winching phase is finished begins when the forwarder driver loads the first log and ends when the forwarder is fully loaded Details of the harvested trees and volume in each treatment are given in Table 5. The system boundaries for the study area were set to those of the harvesting operations, from the felling to the landing site. The Functional Unit (FU) for the analyses was the cubic meter of round wood (m 3 ); in Appendix A the data referring to another Functional Unit (FU) are shown (1 t of fresh mass, following the data shown in Table 2). This is important in order to compare these results more readily with other studies.
Operational costs were estimated according to the Miyata method [34] as previously explained in Spinelli et al. [35]. Economic evaluation of the different machines was carried out taking into consideration different periods of use. The skidder, harvester and forwarder depreciated by 1200 SMH per year [35,36] in a depreciation period of 10 years [37]. The chainsaw for the felling depreciated by 800 SMH per year in a depreciation period of 2 years [36]. Labor cost was set at € 15 SMH −1 inclusive of indirect salary costs [38]. Lubricant consumption was calculated as reported by Picchio et al. [39].
Costs for insurance, repair and service were obtained by literature analysis [35], while the fuel and lubricant prices were taken by a market survey (second semester 2019) conducted upon three company products. The calculated operational cost, as reported in similar studies [35], was increased by 10% to account for overhead costs [40].
Focusing instead on the environmental aspects, an energy consumption analysis was performed, applying the Gross Energy Requirement (GER) method [41]. Indirect input (MJ kg −1 ) of harvesting machinery was evaluated taking into consideration the average energy value of the raw materials. This is related to several parameters, i.e., quantitative presence (%), total mass of the machine (kg), overall service life of the machine (h m −3 ) and use of the machine during harvesting. Energy consumption related to human manpower was evaluated according to what was reported in previous works [42][43][44], through the application of a standard value equal to 0.030 MJ min −1 worker −1 .
To calculate the energy balance, the energy value of poplar wood was determined as Higher Heating Value (HHV) (CEN/TS 14918), on 30 random samples, through Parr calorimeter, model 6200 [45].
Regarding pollutant emissions during logging operations, emissions related to fuel were evaluated as the sum of the emissions during combustion (Efc) and the emissions produced within the production and logistic process (Efp). For Efc assessment of fuel energy content, the emission factor of the engine and the thermal efficiency of the combustion were taken into consideration, as reported in Klvac et al. [46] and Athanassiadis [47].
Dealing with Efp assessment, fuel energy content and emission factors were obtained by [46], but HC emission factor was taken from [47].

Statistical Analysis
The first step in statistical analysis was checking for normality using the Kolmogorov-Smirnov test and for homogeneity of variance using the Levene test. Averages of dendrometric characteristics, skid trail network, and average extraction cycle time elements in each area between the two site conditions (moist and dry soil) were compared by independent t test. A regression analysis of time study data was used to check the model's capability of predicting productivity as a function of statistically significant independent variables such as distance and load size. If the data were not normally distributed, a non-parametric Spearman's rank coefficient was applied to analyze the correlation between the variables.
A major focus was placed on extraction operations investigating the relationships between time elements and dendrometric characteristics of extracted timber, and between time elements and bunching-extraction distance. This investigation was performed by nonlinear regression analysis, performed by SPSS 19.0 software (New York, NY, United States).

Results
Results of the t-test showed no statistically significant differences regarding both extraction distance and mean volume per working cycle between moist and dry sites in all the three plantations. Bunching distance was statistically lower in the IT1 dry site than in the moist one, while in IR and IT2 no statistically significant difference was found for this parameter (Table 6). Table 6. Extraction trail and corridors' average features for the three areas in the two soil moisture conditions (mean ± SD), From the t-test for independent samples applied, statistically significant differences (p < 0.05) between the average values are highlighted (underlined text) (data reported to FU of 1 t of fresh mass showed in Table A2). Data regarding working time analysis of felling and processing operations are given in Table 7 and Figure 1. In all three plantations felling operation time (motor-manual in IR and IT1 and mechanized in IT2) did not show any statistically significant differences between dry and moist sites. There were some differences found in processing and moving time. In every study area the most time-consuming operation was processing.

IR
The soil moisture in IT2 did not affect working productivity, with no statistically significant difference among different soil moisture conditions for every phase, and consequently also for the overall working time. In IR1 and IT1, instead, both TET (Total Effective Time) and TGT (Total Gross Time) were significantly higher in moist conditions than in dry ones, with higher DT (Delay Time) also in the moist soil. However, as previously reported, such differences are related not to the felling operations, but only to moving (both IR and IT1) and processing (only IR).
Working time analysis data in bunching and extraction are reported in Table 8 and Figure 2. The only phase which did not show an effect of soil moisture on productivity was the landing operation (LO), while all the other phases were influenced by the moisture content of the soil. This led to a significant difference in overall working times (both TET and TGT), related to the different site conditions, in all the three yards. In particular, higher soil moisture negatively affected working productivity. Table 7. Description statistics of felling-processing operations for the harvesting sites studied referring to a single tree. From the t-test for independent samples applied, statistically significant differences (p < 0.05) between the average values are highlighted (underlined text). Statistical comparisons are between columns.    The analysis of the various factors influencing extraction time (Table 9) revealed that the parameter with the highest impact on time consumption was bunching-extraction distance, with R 2 values ranging from 0.6 to 0.8 for both dry and moist sites in all the three yards. Increased bunching-extraction distance obviously led to increased bunching-extraction time; however, it is interesting to notice (Figure 3) how this effect is less evident in forwarding operations (IT2) than it is in winching operations (IR and IT1). Focusing on the overall harvesting system productivity (Figure 4), moist soil showed negative effects in all three plantations. In detail, SMH in dry soil conditions was 5.671 m 3 h −1 in IR, 6.403 m 3 h −1 in IT1 and 23.761 m 3 h −1 in IT2; while, respectively, were 12.53%, 18.68% and 16.27% lower in the moist soil. Moreover, the higher moisture content of soil also resulted in a higher percentage difference between PMH and SMH, i.e., 2.73% vs. 3 Focusing on harvesting costs, the results of the economic evaluation carried out within the present study are given in Tables 10 and 11.

Sites
The details of hourly costs reported in Table 10 show how the harvesting machinery applied in IT2 (harvester and forwarder) presents higher hourly costs, mostly related to the higher purchase price. However, the higher productivity of this fully mechanized harvesting system allowed IT2 to have a lower cost per m 3 of timber produced (Table 11).
Regarding the influence of soil moisture conditions on harvesting costs, it is evident that the negative influence on working performance correlated to higher moisture also led to higher harvesting costs. In detail, this was about 16%, 26% and 16% higher in the moist site than in the dry one for IR, IT1 and IT2, respectively.  Figure A1).  Table 11. Harvesting costs for one cubic meter of wood and percentage of costs at two main operations (felling-processing and bunching-extraction to landing) in the studied sites (data reported to FU of 1 t of fresh mass showed in Table A3). Regarding environmental aspects, the results of the analysis of energy efficiency are given in Table 12. The highest energy input was reported for IT2, due to the complete mechanization of the overall harvesting operations. The effects of moisture on environmental performance can be observed in all of the three plantations where higher soil moisture led to lower energy efficiency, more exactly, 97.7% vs. 97.0% in IR; 98.1% vs. 96.9% in IT1 and 97.0% vs. 96.6% in IT2. In IR and IT1 a major part of the energy input is related to bunching-extraction operations in both moist and dry soil conditions. Instead, in IT2, felling operations via harvester were the reason for the highest portion of energy input in this yard, in both soil conditions ( Figure 5). Interestingly, moist soil led to an increased portion of energy input related to bunching-extraction in all three yards (77.2% vs. 81.6% in IR; 74.6% vs. 84.0% in IT1 and 25.8% vs. 32.7% in IT2), showing how this operation was most influenced by soil moisture when regarding environmental issues.

Description
As shown in Table 13 and Figures 6 and 7, soil moisture also showed negative effects regarding pollutant emissions in the three different yards, with increasing emissions for all the investigated parameters in moist soil conditions. What is more, mechanized felling via harvester (IT2) led to higher emissions in comparison to motor-manual felling (IR and IT1).   6. Percentage distribution of total PM10 emission in the harvesting sites studied, data shown for single operation (data reported to FU of 1 t of fresh mass showed in Figure A2).

Figure 7.
Percentage distribution of GHG emission in the harvesting sites studied, data shown for single operation and reported in CO 2 equivalent (data reported to FU of 1 t of fresh mass showed in Table A3).

Comparison of Harvesting Systems Performance
Several studies on work productivity and cost analysis in poplar plantations are available in the literature, even if most of these deal with Short Rotation Coppice (SRC) plants for bioenergy production [48]. Indeed, few studies have focused on productivity analysis in poplar plantations for timber production. In a poplar plantation located in Serbia, Danilovic et al. [49] reported a productivity for mechanized felling-processing via harvester of 30.3 to 34.7 m 3 h −1 , depending on working method and stem dimension. In the same year, Spinelli et al. [50] carried out an extensive productivity and cost analysis in a 25 year old poplar plantation in Italy, reporting an average work productivity (SMH) for motor-manual felling and processing via chainsaw of 6.3 m 3 h −1 and a value of 21.1 m 3 h −1 for the same operation performed via harvester.
The productivity values found in the present study are higher than reported in the above cited studies, both for motor-manual and mechanized felling processing. This difference is more pronounced in comparison to Spinelli et al. [50]. Such a gap concerning mechanized felling-processing can be partially related to the lower average dbh of the stems in the previous study (around 30 cm vs. 40.4 cm). However, the major difference between the present studies and the literature, which can explain the higher productivity found in the present analysis, is the lower percentage of delay times. In particular, delay percentage ranges from 1.7% in IR to 2.9% in IT2, while the average delay for motor-manual felling-processing in Spinelli et al. [50] was 29.6%, decreasing to 13.0% in mechanized operations, while in Danilovic et al. [49] delay accounted for 28.5 of the working time.
Concerning bunching-extraction operations, no study on winching and forwarding in poplar plantations for timber production were found in the literature. However, there are several studies on bunching and extraction via TimberJack 450 cable skidder, which analyzed work productivity in different forest stands. Lotfalian et al. [51] found a winching productivity of 20.2 m 3 h −1 in beech high stand thinning with an average extraction distance of 289 m, while Mousavi [52] reported a bunching-extraction productivity of about 11 m 3 h −1 in beech selection cutting with an average skidding distance of 439 m. Nikooy et al. [18] reported instead a lower value of 5.2 m 3 h −1 productivity of Timberjack 450C in timber extraction of path cutting in a pine plantation. Such lower productivity is probably related to the lower dimension of trees considering that skidding distance was comparable to those in both IR and IT1 [18]. Therefore, as a general trend, bunchingextraction productivity in the present study showed higher values than previous works reported in the literature for the same machinery in different stands and silvicultural interventions. This difference is related to both the type of forest intervention (clear cut) and to the flat terrain of IR and IT1, which facilitated logging operations.
Regarding forwarding operations, also in this case it is possible to make a comparison regarding work productivity only between different stands, considering the lack of studies on poplar plantations for timber production. The forwarder is the most commonly applied machinery in CTL (Cut to Length) harvesting operations, and it has been widely applied in artificial plantations, mostly of softwood species [53]. This machinery can reach a very high working productivity [54,55], even if a proper planning of the intervention is needed to reduce the impact which can occur considering the average dimension of a forwarder [13]. Comparing the findings of the present study with other forestry interventions in artificial plantations, it is evident that the substantial average dimension of stems, the short extraction distance and the flat terrain features in IT2 led to higher work productivity. In detail, Puttock et al. [56] reported a SMH productivity of 11.2 m 3 h −1 in a poplar-dominated mixed-wood stand in Southern Ontario during thinning interventions; Eriksson and Lindroos [57] showed PMH productivity for forwarding in clear cutting in pine and spruce stands of 21.4 m 3 h −1 . Another study performed in Romania reported a SMH productivity of 15.35 m 3 h −1 in a clear cut of spruce stand, with an average slope of 10% and an average extraction distance of 479 m [58]. In another study recently carried out in Poland, Magagnotti et al. [11] reported a forwarding productivity of 24.4 m 3 h −1 SMH for a poplar plantation.
Focusing on harvesting costs, it is possible to notice how felling and processing operations accounted for the major part of these in all the yards in both soil moisture conditions, as reported by previous literature [38,59]. Felling and processing costs, with values ranging from 4.22 € m −3 (IT2 dry) to 5.77 € m −3 (IR moist), were in line with the literature findings for several Italian poplar plantations for high value timber production, notwithstanding higher work productivity. Spinelli et al. [50] reported a unit cost of about € 5 m −3 for both motor-manual and mechanized felling-processing. This can be explained by the higher purchase costs of the machinery applied in the studied plantations. Skidding and forwarding costs are also in line with the literature, for similar harvesting systems but in different kinds of stands. In the present study the lowest cost for extraction was shown by IT2 dry at € 2.46 m −3 , while the highest cost was related to winching in IT1 moist (€ 5.28 m −3 ). Regarding winching operations through cable skidder, Jourgholami and Majnounian [33] reported 6.15 € m −3 as the cost of Timberjack 450C in timber extraction on a pine plantation, while the findings of Lotfalian et al. [51] showed 5.15 USD m −3 . Focusing on forwarding, extraction costs were assessed by Cabral et al. [60] at about 1.95 € m −3 , while about 7.5 € m −3 were reported by Kaleja et al. [61], and about 9.2 € m −3 by Magagnotti et al. [11].
Concerning environmental impact, a comparison was carried out between the findings of the present study and other similar studies, regarding both high and medium level of mechanization. In both cases, energy inputs and energy balance in the investigated poplar plantations were substantially higher [37,62].
System efficiency values were high in all three yards in both the soil conditions (ranging from 96.6% to 98.1%), and thus were in line with previous literature findings in other forest interventions [16,62,63].
GHG emissions, mostly regarding CO 2 , were lower than in previous literature findings, which reported a range between 3 and 33 kg CO 2 eq, while the values of the present work ranged from 1.6 to 2.8 kg CO 2 eq [5,16,24,46,64].

Influence of Soil Moisture on Harvesting Performance
There is a considerable amount of literature regarding working productivity evaluation, but a major part of the studies focused on the influence on working performance of parameters such as terrain features, working distance, age, species composition, labor skills, etc. However, not much attention has been directed towards soil moisture conditions and productivity. Studies were, however, conducted on soil impact related to logging activities [65][66][67][68].
In all of the three yards, higher soil moisture led to lower work productivity, thus to higher harvesting costs, in accordance to what was reported in the few studies dealing with this topic [69,70]. High soil moisture negatively affected both motor-manual and mechanized operations, resulting in higher working times, except for felling (both with chainsaw and harvester) and processing (only with harvester). Longer working time in higher soil moisture conditions for mechanized operations (bunching -extraction) is related to the lower driving speed which the machinery was able to achieve with moist soil. The high level of moisture of the terrain caused reduced tire grip and the operators had to reduce the working speed for safety reasons. Interestingly, this did not happen for felling-processing operations by harvester, which were not affected by soil moisture regarding working time. Soil moisture also negatively affected motor-manual operations, specifically, moving (IR and IT1) and processing (only IR). This is equally related to worker safety issues, with operators that had to be more cautious during the logging activities, due to the fact that high moisture in the soil made the trail slippery, and therefore prone to accidents.
High soil moisture showed negative effects on the environmental performance of logging activities in the studied poplar plantations. There were two reasons for these negative effects: the longer working time, which required more time in which motors were running, and the higher torque needed to move the machinery in moist soil, considering the lower grip and higher attrition. Thus, there were higher emissions of pollutants.

Conclusions
Although poplar plantations are important sources of timber in both Iran and Italy, few studies have focused on work productivity evaluation under different aspects related to sustainability in such kinds of stand. Moreover, these plantations are often located in plain or floodplain lands, therefore harvesting operations can occur in soil conditions with a high moisture content.
The aims of this study were: (i) to evaluate different harvesting systems in poplar plantations and (ii) to evaluate the influence of soil moisture on economic and environmental performance of logging operations.
In order to assess different harvesting systems in poplar plantations, from what was analyzed it is possible to state that a fully mechanized harvesting system (harvesting forwarder) is the most productive and economically sustainable with respect to semimechanical harvesting/processing and skidding extraction. However, in terms of energy balance and emissions, it is possible to state exactly the opposite, that the best harvesting system was semi-mechanical harvesting/processing and skidding extraction. These are aspects to be carefully considered during operations planning, but they must also be analyzed in terms of greater efficiency of the mechanization used, seeking to bring mechanical technologies that are increasingly efficient also in environmental terms to the forestry sector.
In order to assess the influence of soil moisture on economic and environmental performance of the logging operations, the findings revealed that high moisture content led to lower work productivity in all of the three investigated plantations, with detrimental effects on harvesting costs, which were found to be higher in moist soil conditions in all three yards. Moreover, environmental features related to pollutant emissions were higher in moist soil conditions, as a consequence of the longer time and the major torque required for the machinery to perform the logging activities.
It can be concluded from these findings that it is advisable to avoid logging operations in conditions of high soil moisture (>20%), to decrease the impact on the soil, to create higher cost-effectiveness, and to reduce the emissions of pollutants.

Data Availability Statement:
The data presented in this study could be available on request from the corresponding author.

Appendix A
The data assessed in the Appendix A are referred to another Functional Unit (FU) respect to how reported in the main text. The FU in this case was 1 t of fresh mass (following the data showed in Table 2). This was important in order to give more possibility to compare these results with other studies.  Figure A1. Average yard productivity (bars) and possible increase of performance (lines) from SMH to PMH for the six harvested sites. Figure A2. Percentage distribution of total PM10 emission in the harvesting sites studied, data shown for single operation, referring to fresh mass. Figure A3. Percentage distribution of GHG emission in the harvesting sites studied, data shown for single operation and reported in CO 2 equivalent, data referring to fresh mass.