Strip Clear-Cutting Application and Logging Typologies for Renaturalization of Pine Afforestation — A Case Study

Renaturalization treatment in black pine afforestation is an important topic that should be considered. There is a need to favor the evolution of artificial pine forests toward natural forest systems. Overall, this study focused on pine forests, and suggests one typology of clear-cutting (dismantling cutting) on strips, which is associated with different extraction management techniques. Some ecological and environmental aspects associated with renaturalization treatments that have been applied by different mechanizations in black pine afforestation have been highlighted, as well as how renaturalization and the active ecological management of these stands could affect soil and vegetation. The main objectives of this research were to: (1) analyze the impact of silvicultural treatment and logging activities on forest soil, and (2) assess tree regeneration and floristic biodiversity in an ecological management system, in terms of both quantity and quality characteristics. These analyses were planned to obtain an overview of the environmental impact related to a multifunctional approach to the forest management of black pine afforestation. Essentially, the answers to the main research questions are: (1) less invasive extraction systems seem to use a cable yarder and forest winch; (2) a clear soil recovery trend with good capabilities is visible, in particular for the two extraction systems by cable; however, over a three-year period, only a partial but substantial recovery has been shown; (3) in general, silvicultural treatment showed qualitative and quantitative improvement in terms of tree regeneration; in particular, the extraction systems by forest winch and cable yarder showed better results; (4) silvicultural treatment seems not to have led to improvement at the level of the herbaceous and shrubby layers; however, clear differences are shown among the different harvesting systems. Significant recovery after use of the cable yarder was observed.


Introduction
Among the different coniferous species, black pine (Pinus nigra Arnold) has been successfully used in Italy since the beginning of last century for afforestation, as well as the recolonization of pastures Forests 2018, 9, 366 2 of 24 and degraded environments due to its favorable frugal characteristics [1,2].Black pine is appropriate for afforestation in shallow soils with severe climatic conditions such as prolonged summer drought, due to its root system extension [1].In 2005, Pinus nigra forests covered an area of approximately 240,000 ha (about 3% of the forested area in Italy) [3].This afforestation provides soil protection and the restarting of dynamic and evolutionary processes of vegetation.
Tonon et al. [4] highlighted that black pine afforestation modified microclimate and soil conditions, making them suitable for more discerning species.Nevertheless, today many of these afforested areas are in a poor physical and biological state with no dynamic processes.This is due to the combination of several factors (i.e., the inadequate treatment and biotic and environmental adversity).Thus, these areas require urgent silvicultural action in the form of renaturalization [5,6] to ensure the permanence of land cover, with the gradual replacement of pine with late successional tree species that are typical of more mature stages of evolution (native broadleaves).Furthermore, the woody material that was required as a source of renewable energy opens up interesting prospects for the use of biomass obtainable by the renaturalization of these types of forests.
The application of silvicultural treatments has important effects on all of the ecological aspects [7][8][9].Soil characteristics and tree functioning may be affected by harvesting due to the modified input of light (quality and quantity), temperature, and water availability.Moreover, it is important to highlight that logging activities, with machine and load traffic, may cause soil compaction, soil horizon mixing, and topsoil removal [10].Compaction reduces soil porosity and the connectivity of pores, thus increasing soil density and shear strength [9,11,12].Thus, soil degradation could reduce tree growth [13], and carbon dioxide efflux from the soil may change significantly [14].After a renaturalization treatment, the restocking of the stand is left to two main natural regeneration processes: through seedlings (gamic) and sprouts (agamic).Due to the high sprouting potential of some broadleaves growing under the canopy cover of the pine forest, the dynamics of recruiting new individuals to existing stands has largely been overlooked [15].This aspect will be of particular interest for these new forests that are composed of a mosaic of artificial and natural stands.
Biodiversity conservation has long been a goal of European conservation policy [16,17].However, despite more than 25% of European land having been given some level of protection for conservation, biodiversity continues to decline [18].One factor contributing to this decline may be unsuitable management practices in protected areas, or, more specifically, an insufficient application of evidence-based conservation recommendations [18,19]; for example, a lack of diversity of management systems with a total abandonment of traditional systems, without considering the possibility of improving them.Biodiversity monitoring is essential in order to support management decisions in maintaining multiple forest ecosystem functions [20].A better understanding of the roles of the components of biological diversity that support the provision of multiple forest ecosystem services is necessary [21,22].
Logging systems may differ depending on silvicultural management and the final products.The technical and economic utilization of forests depend on various factors concerning terrain conditions, transportation networks, and harvesting technologies, as well as silvicultural treatment and forest operation systems [23,24].Although in recent times there have been significant innovations in forest operations [8,25] in terms of both technology and methodology, logging activities in Italy are still developed by traditional methods, i.e., motor-manual felling with chainsaws and the use of mules and/or agricultural tractors for extraction [25][26][27].
Actually, studies focusing on the effects of silvicultural treatment and utilization on soil and plant biodiversity are rare, or are overly guided by extreme currents of thought.Without some clear effects and recovery times, it is very difficult to respect the concept of sustainability.Also for these reasons, companies are seldom required to take into account the impact of their operations on the land and on forest sustainability, or, rather, to consider the real application of sustainable forest management as suggested by Forest Europe.For these reasons as well as the fundamental purposes of renaturalization, it is important to know and assess the possible impacts on the soil and biodiversity of silvicultural treatments and logging operations, considering that one of the most important problems of the forest sector is to minimize ground damage caused by forest operations [28] and biodiversity loss.In general, the effects of harvesting include changes in vegetation, nutrient availability, and soil microclimate and structure, as well as litter quantity and quality [28].In particular, forest operations, such as forwarding and skidding, have a high potential for soil compaction [29][30][31].However, adequately managed forest ecosystems are claimed to be highly resilient in the long term [32].
The modern silvicultural approaches lead to aggregated retention, which is the practice of leaving unlogged patches within logged areas, and is a form of retention forestry that has gained increasing use globally as an alternative to clear-cut logging [33][34][35].In contrast to clear-cutting, retention forestry maintains habitats for species affiliated with a closed forest-thus mitigating the negative effects of timber harvest-while also providing a habitat for early seral species [36,37].These approaches are characterized by a set of fundamental principles, including an avoidance of clear-cutting, an emphasis on structural diversity and small-scale variability, the deployment of mixed species with natural regeneration, and an avoidance of intensive site-preparation methods [38].However, although these modern approaches find real justification in the management of "natural" forests, they have some complications in the renaturalization practices or management of forest monospecific plantations.In these cases, which are similar to problems with new silvicultural approaches [38], a suite of ecological, economic, logistical, informational, cultural, and historical constraints currently hinder the wider adoption of alternative silvicultural practices.Individual contexts display their own unique combinations and the relative significance of these constraints, and accordingly, targeted efforts, such as regulations and incentives, may help overcome specific challenges.In a broader context, the possibility of clear-cutting treatments on adequate surfaces and with proper methodologies might provide an additional flexibility and facilitate renaturalization within a broader set of ecological circumstances.
In order to improve silvicultural management and logging methods, better knowledge of the long-term impact of forest operations is needed [26,39].Reducing the negative effects of felling and extraction is one of the main goals of sustainable forest management [40][41][42].
The main objectives of this research were to (1) analyze the impact of silvicultural treatment and logging activities on forest soil, and (2) assess tree regeneration and floristic biodiversity in an ecological management system, in terms of both quantity and quality characteristics.
To analyze the impact on soil and short-term recovery in particular, in addition to physical and chemical analyses of soil (i.e., pH, organic matter, bulk density, penetration and shear resistance) [30], an innovative arthropod-based Soil Biological Quality index, QBS-ar [43,44], was used.This is a valuable tool in ecosystem restoration programs for monitoring the development of soil functions and biodiversity, and preventing the negative effects of soil compaction that result from logging activities.
Moreover, other ecological and environmental aspects associated with renaturalization treatments were evaluated in black pine afforestation in relation to tree regeneration characteristics and vegetation dynamics.In particular, the most important points in investigating the renaturalization and active ecological management of these Pinus nigra forests were as follows:

•
Three mainly different logging methodologies and their immediate impact on some soil features • Possible recovery time of the impact on soil features over a range of three years

•
The effect of silvicultural treatment and logging activities on natural tree regeneration • First analysis of silvicultural treatment and logging activities on plant structure and functioning These analyses were planned so as to obtain an overview of the environmental impacts related to a multifunctional approach to the forest management of black pine afforestation.

Study Site
The studied areas are located near Passo delle Capannelle Municipality of Pizzoli (AQ 42 • 26 49 N, 13 • 20 14 E) (Figure 1).Abruzzo is a region in Italy that has a fairly developed forest sector and a huge forest surface (about 41% of the region) accounting for about 4% of the entire Italian forest surface.In this region, there are about 19,158 ha of coniferous plantations, and of these, black pine afforestation covers about 13,004 ha [44].The pine afforestation that was studied covers about 27 ha along the middle mountain slope.The two blocks that were chosen for the experiment are located on the southeastern slope in the altitudinal range of 1200-1300 m a.s.l.(above sea level), with an average slope of about 50%.On the slope are three different calcareous formations dating back to the Triassic-Jurassic: lower limestone calcareous, ammonite green, limestone, Posidonia marl, and Corniola.In some portions, there is generally stratified and cemented Pleistocene stratum debris consisting of gravels and breccias, predominantly limestone.

Study Site
The studied areas are located near Passo delle Capannelle Municipality of Pizzoli (AQ 42°26′49 N, 13°20′14 E) (Figure 1).Abruzzo is a region in Italy that has a fairly developed forest sector and a huge forest surface (about 41% of the region) accounting for about 4% of the entire Italian forest surface.In this region, there are about 19,158 ha of coniferous plantations, and of these, black pine afforestation covers about 13,004 ha [44].The pine afforestation that was studied covers about 27 ha along the middle mountain slope.The two blocks that were chosen for the experiment are located on the southeastern slope in the altitudinal range of 1200-1300 m a.s.l.(above sea level), with an average slope of about 50%.On the slope are three different calcareous formations dating back to the Triassic-Jurassic: lower limestone calcareous, ammonite green, limestone, Posidonia marl, and Corniola.In some portions, there is generally stratified and cemented Pleistocene stratum debris consisting of gravels and breccias, predominantly limestone.A general climatic analysis of the area was conducted on the basis of data available from the meteorological stations of Montereale and Assergi, relating to the years 1950-2000.Rainfall has a bimodal regime, with two annual maximums, in April and November.The driest month is usually July, with about 40 mm of precipitation, and only six rainy days.The winter cold has remarkable peaks, with monthly averages of minimum air temperatures lower than 0 °C for four months; average daily temperatures below 10 °C are exhibited for six months.From the study of the Mitrakos diagrams, it appears that the maximum cold stress was reached only in the month of January, whereas aridity stress was in July.The classification in Rivas-Martínez [45] was used for the bioclimatic description, which defines the considered area as a continental temperate bioclimate, and more precisely as within the attenuated subcontinent subtype.The thermo type was found to be supratemperate lower and ombrotype humid.Additional details regarding temperature and precipitation values throughout the years of the study period are reported in Table 1.A general climatic analysis of the area was conducted on the basis of data available from the meteorological stations of Montereale and Assergi, relating to the years 1950-2000.Rainfall has a bimodal regime, with two annual maximums, in April and November.The driest month is usually July, with about 40 mm of precipitation, and only six rainy days.The winter cold has remarkable peaks, with monthly averages of minimum air temperatures lower than 0 • C for four months; average daily temperatures below 10 • C are exhibited for six months.From the study of the Mitrakos diagrams, it appears that the maximum cold stress was reached only in the month of January, whereas aridity stress was in July.The classification in Rivas-Martínez [45] was used for the bioclimatic description, which defines the considered area as a continental temperate bioclimate, and more precisely as within the attenuated subcontinent subtype.The thermo type was found to be supratemperate lower and ombrotype humid.Additional details regarding temperature and precipitation values throughout the years of the study period are reported in Table 1.The black pine (Pinus nigra Arnold subsp.nigra var.italica Villetta Barrea) plantation studied was about 60 years old.The soil preparation was done in steps (called "gradoni", or small terraces) 2 m to 4 m apart, with numerous stone artefacts.The planting was carried out with bare root black pine transplants at a distance of 1 m in the step.It was a pure stand, with poor social differentiation and a high slenderness ratio, which is a clear sign of a lack of thinning.The degree of coverage was high (90-100%) and homogeneous.The distance between the trees was similar to the original planting scheme.The shrubby undergrowth was very poor, while dense brachypodium carpets and mosses were intertwined with needles and cones.An analysis conducted on the wooden bores of some trees showed that the annual growth had undergone a significant slowdown in the last eight years (fractions of mm), accentuating the slowdown that occurred 16 years ago, and making it difficult to hypothesize a reaction to any thinning.The limit of 60-65 years represents a threshold beyond which the reaction capacity of the stand to the thinning is greatly reduced [1].In the considered afforestation, no significant meteoric damage had occurred, there were no obvious signs of fungal and insect attacks, dead wood snags were substantially absent, and logs were not consistent.

Renaturalization Treatment and Logging Methods
Among the possible renaturalization techniques [46], taking into account the stand age, the need to combine harvesting cost-effectiveness, that it was solid ground despite being on a high slope, and that the black pine forest was sheltered from northern winds, clear-cutting (dismantling cutting) on strips was chosen.This consisted of a series of strips 15 m wide and 100 m long, alternating with uncut strips (to eliminate 50% of the surface of the pine plantation).The width was proportional to the tree height, while the length was proportional to the logging systems and planning.The choice of the two experimental blocks was made through a design-based approach, which is a statistical approach that establishes the methods of choice and use of the sites, allowing possible pseudoreplication problems to be overcome [47].The two experimental blocks (replicates) were delimited upstream of a forest road, and they differed in altitude location, forest growth conditions, and soil texture.The first experimental block, which was at an altitude of 1200 m a.s.l.(east-southeast), consisted of 12 strips that were 100 m long (according to the lines of maximum slope) and 15 m wide.This block was surrounded on all sides, excluding the track, with a protection buffer that was a minimum of 20 m wide.The second block, with similar characteristics, was realized slightly lower, at an altitude of 1100 m a.s.l.(southeast).The logging treatments within each block were randomly assigned, while the silvicultural treatments were systematically assigned (one uncut strip and one clear-cut strip).Within each block, the silvicultural treatment was applied on six alternating strips (one harvested and one untouched) to have six silvicultural repetitions and six control areas.In the harvested strips, three logging treatments Forests 2018, 9, 366 6 of 24 were applied, with two repetitions for every replicate (block).The different extraction methods were characterized by different machines and techniques (detailed in Table 2): (a) animal skidding downhill by TPR (heavy rapid skidding) horse, (b) winching downhill by a forestry-fitted farm tractor using snatch block, and (c) yarding downhill by a cable yarder with automatic carriage.The two blocks (replicates) differed in altitude, soil texture, and dendrometric characteristics.

Analytical Methods
A dendrometric analysis was performed by measuring the diameters at breast height (DBH) of all of the trees, and the plant height of 10% of the trees, which allowed constructing the hypsometric curve and then determining the growing stock by applying the model tree method.
In order to determine the soil particle size distribution for the studied areas, six soil samples in each block were randomly sampled from the top 30 cm of the mineral soil, which was considered to be a crucial indicator of vulnerability to soil compaction [48].Rock fragments (particles >2 mm in diameter) were removed from the air-dried samples for particle size distribution.Afterwards, three sand fractions-2.00-0.50mm (coarse), 0.50-0.25 mm (medium), and 0.25-0.05mm (fine)-were separated by wet sieving.Finally, silt and clay were determined using a hydrometer [31].
Three linear transects, consisting of 1 m × 30 m rectangles for every harvested strip, were tracked to estimate the tree regeneration community composition and perform quantitative analysis.A systematic sampling method was applied as shown in Figure 2.These post-operation analyses were conducted using research methods based on internationally shared protocols, elaborated, and adapted to this context of study as proposed in [8,9,25,26,49].

Analytical Methods
A dendrometric analysis was performed by measuring the diameters at breast height (DBH) of all of the trees, and the plant height of 10% of the trees, which allowed constructing the hypsometric curve and then determining the growing stock by applying the model tree method.
In order to determine the soil particle size distribution for the studied areas, six soil samples in each block were randomly sampled from the top 30 cm of the mineral soil, which was considered to be a crucial indicator of vulnerability to soil compaction [48].Rock fragments (particles >2 mm in diameter) were removed from the air-dried samples for particle size distribution.Afterwards, three sand fractions-2.00-0.50mm (coarse), 0.50-0.25 mm (medium), and 0.25-0.05mm (fine)-were separated by wet sieving.Finally, silt and clay were determined using a hydrometer [31].
Three linear transects, consisting of 1 m × 30 m rectangles for every harvested strip, were tracked to estimate the tree regeneration community composition and perform quantitative analysis.A systematic sampling method was applied as shown in Figure 2.These post-operation analyses were conducted using research methods based on internationally shared protocols, elaborated, and adapted to this context of study as proposed in [8,9,25,26,49].The vegetation study was set up according to the phytosociological method [44].The surveys were performed in summer, when most of the species were flowering.The analysis was carried out in the entire slope of the forest to investigate the potential vegetation and identify the relationships The vegetation study was set up according to the phytosociological method [44].The surveys were performed in summer, when most of the species were flowering.The analysis was carried out Forests 2018, 9, 366 7 of 24 in the entire slope of the forest to investigate the potential vegetation and identify the relationships between plant associations and topography that were established between the different vegetation series.Vascular species were recorded with the percentage of coverage-abundance for each area.The phytosociological surveys carried out in the reforestation area allowed us to characterize the current state of the underbrush.The analysis of the groups of species related to the phytosociological syntaxa allowed us to interpret the current phase of reforestation dynamics and its location in the vegetation series that characterizes the area in which the pinewood was located.The floristic study took place through sample collections before cutting.The samples were determined using the main national [50][51][52][53][54] and European flora [54,55].
The tree composition of natural regeneration was analyzed by the Species Importance Value (SIV) index, which was calculated for each species, as reported in the literature [48,56] (Table 3).The regeneration species diversity index was computed using the Shannon-Wiener information function (Table 3).The natural regeneration and floristic biodiversity of tree composition were also assessed using species richness.The Shannon index and the evenness index were used for floristic biodiversity, without considering the tree layer before cutting.Biodiversity indices were compared by year and type of intervention.Species richness was defined as the total number of species sampled.The Shannon index [57] is based on information theory and the degree of difficulty in accurately predicting the next species sampled.This diversity index is sensitive to changes in rare species, has good discriminant ability, and is well represented in the literature [58].Biodiversity indices were compared by year and type of intervention.The structural evenness index was calculated as reported in [59] (Table 2).This index varies between 0 and 1, where a value of 1 symbolizes an exact uniform distribution.The Shannon index is a model that measures species diversity and the degree of homogeneity in species abundance.One of its applications is to correctly estimate the anthropogenic impact on ecological systems.
The impact on soil due to silvicultural treatment and forest operations was assessed by a systematic sampling method selecting three sample plots (SPs) for each harvested strip, one every 30 m along the strip axis, starting at 15 m from the lower border.Each SP consisted of a circular area 12 m in diameter.In each SP, three tests for every parameter were done, which measured: bulk density (BD), pH, organic matter content (OM), inorganic carbon content (IC), penetration resistance (PR), and shear resistance (or strength) (SR), and QBS-ar index as detailed in [43][44][45][46][47].For the two blocks, two different control areas were selected at the edge of the blocks, which were at least 30 m from the edges.In each area, three SPs were randomly selected, and three tests for all of the parameters were done on each of them.
Natural regeneration analyses were done in 2016 and 2017, which were one and two years after harvesting, respectively.Soil impact analyses were done in 2015, immediately after the logging operations, and in 2016 and 2017, one and two years after harvesting, respectively.All of the analyses were replicated in the same months.

Statistics
Statistical analyses were carried out using Statistica 7.1 (2007) software (StatSoft Inc., Tulsa, OK, USA).As a first step, data distribution was plotted and checked for normality (Lilliefors) and homogeneity of variance (Levene test).All of the data points then underwent a t-test, an ANOVA, or a MANOVA to test the effects of the different treatments.In order to determine the relation between QBS-ar, BD, PR, and SR, a non-parametric correlation analysis (Spearman correlation matrix) was applied.The data, which were not normally distributed and had insufficient homogeneity of variance, were statistically processed using the non-parametric ANOVA Kruskal-Wallis test.Principal component analysis (PCA) was applied to investigate any linear correlations between the expressed soil conditions of the main six characteristics studied, and also between the tree regeneration situations of the main four characteristics studied.To minimize the scaling effect due to the different measurement units, the data corresponding to each independent variable were standardized using Box-Cox lambda.

Silvicultural Treatment and Dendrometric Analysis
The silvicultural treatment applied was clear-cutting (dismantling cutting) on strips, with harvesting of 50% of the surface of the pine plantation.The strip dimensions were proportioned to the tree height and length to the logging system and planning.This treatment aimed to achieve the renaturalization of these forests within a sustainable management system.
The dendrometric data collected before the harvesting of the two blocks showed average values, which were statistically similar for trees age, DBH tree height, basal area, tree density, and aboveground biomass stock (Table 4).

Analysis of Stand Regeneration
No statistical difference was observed between the two blocks, and the nine species were indifferently present (p > 0.05).The three logging treatments and the control showed different percentages (ANOVA p < 0.05) of species distribution (Figure 3).The tree species most represented were Sorbus aria L., Fraxinus ornus L., Quercus pubescens Willd., and Quercus cerris L. According to the Tukey test, the strips extracted by cable yarder were statistically grouped with the control data, while the other two treatments were statistically different.Pinus nigra regeneration was largely present in Forests 2018, 9, 366 9 of 24 the strips extracted by horse (9%) and winch (7%).There was an exotic tree present in these strips (Robinia pseudoacacia).
Forests 2018, 9, x FOR PEER REVIEW 9 of 24 present in the strips extracted by horse (9%) and winch (7%).There was an exotic tree present in these strips (Robinia pseudoacacia).One year after harvesting, the highest natural regeneration density was found for the strips extracted by winch and cable yarder: 23,580 trees/ha and 21,150 trees/ha, respectively (Figure 4).The strips extracted by horse showed the lowest value, 6250 trees/ha, and the control had stable values ranging between 10,260-11,150 trees/ha.Two years after harvesting, the highest regeneration density was found for the strips extracted by winch: 20,650 trees/ha, which was slightly lower than in 2016.The strips extracted by horse and cable yarder showed similar values ranging between 18,850-19,170 trees/ha; the first treatment had a considerable increase, while the second one had a slight decrease compared with 2016.One year after harvesting, the highest natural regeneration density was found for the strips extracted by winch and cable yarder: 23,580 trees/ha and 21,150 trees/ha, respectively (Figure 4).The strips extracted by horse showed the lowest value, 6250 trees/ha, and the control had stable values ranging between 10,260-11,150 trees/ha.Two years after harvesting, the highest regeneration density was found for the strips extracted by winch: 20,650 trees/ha, which was slightly lower than in 2016.The strips extracted by horse and cable yarder showed similar values ranging between 18,850-19,170 trees/ha; the first treatment had a considerable increase, while the second one had a slight decrease compared with 2016.present in the strips extracted by horse (9%) and winch (7%).There was an exotic tree present in these strips (Robinia pseudoacacia).One year after harvesting, the highest natural regeneration density was found for the strips extracted by winch and cable yarder: 23,580 trees/ha and 21,150 trees/ha, respectively (Figure 4).The strips extracted by horse showed the lowest value, 6250 trees/ha, and the control had stable values ranging between 10,260-11,150 trees/ha.Two years after harvesting, the highest regeneration density was found for the strips extracted by winch: 20,650 trees/ha, which was slightly lower than in 2016.The strips extracted by horse and cable yarder showed similar values ranging between 18,850-19,170 trees/ha; the first treatment had a considerable increase, while the second one had a slight decrease compared with 2016.The SIV detailed analysis for 2017 (Figure 5) shows that the most important tree regeneration was for the species Q. pubescens and Sorbus aria.The worst SIV was shown by Q. cerris and Fraxinus ornus.Fraxinus ornus only showed an SIV peak in the strips extracted by horse.The analysis of the total SIV (Table 5) showed that the strips harvested by horse and winch increased compared to the control, while for the strips harvested by cable yarder, the trend was positive, but the value was lower than the control.
Tree species diversity, as tested by the Shannon-Wiener and the evenness indices, showed different situations.First, there was a general slightly positive trend from 2016 to 2017 for the control and the strips harvested by winch and cable yarder, while the index was stable for the strips harvested by horse.The highest diversity was found for the strips harvested by cable yarder and winch (in particular for the evenness index, amounting to 0.88 ± 0.07 and 0.85 ± 0.09, respectively, as shown in Table 5), while the lowest diversity was found in the control (Shannon index ranged from 1.68 to 1.70 and the evenness index ranged from 0.73 to 0.77, Table 5).The SIV detailed analysis for 2017 (Figure 5) shows that the most important tree regeneration was for the species Q. pubescens and Sorbus aria.The worst SIV was shown by Q. cerris and Fraxinus ornus.Fraxinus ornus only showed an SIV peak in the strips extracted by horse.The analysis of the total SIV (Table 5) showed that the strips harvested by horse and winch increased compared to the control, while for the strips harvested by cable yarder, the trend was positive, but the value was lower than the control.
Tree species diversity, as tested by the Shannon-Wiener and the evenness indices, showed different situations.First, there was a general slightly positive trend from 2016 to 2017 for the control and the strips harvested by winch and cable yarder, while the index was stable for the strips harvested by horse.The highest diversity was found for the strips harvested by cable yarder and winch (in particular for the evenness index, amounting to 0.88 ± 0.07 and 0.85 ± 0.09, respectively, as shown in Table 5), while the lowest diversity was found in the control (Shannon index ranged from 1.68 to 1.70 and the evenness index ranged from 0.73 to 0.77, Table 5).

Analysis of Floristic Biodiversity
The most consistent group of species, by both qualitative and quantitative analysis, taking into account coverage, was arid secondary grazed meadows (Figure 6).The pool of diagnostic species of the Apennine Phleo ambigui-Bromion erecti alliance was highlighted, among which the graminaceous Brachypodium rupestre was relevant, forming a fairly uniform mat.It was accompanied by Cytisus spinescens, Asperula purpurea subsp.purpurea, and Galium lucidum subsp.lucidum.The group of plants belonging to the superior categories of the suborder Artemisio albae-Bromenalia, order Phleo ambigui-Brometalia erecti, and class Festuco-Brometea was very abundant.Among these, the species forming the structure of secondary grasslands in the Apennines was Bromus erectus.Even if this vegetation was attributable to the Phleo ambigui-Bromion erecti alliance, the difference was the presence in the pinewood of a community dominated by Brachypodium rupestre, which was absent in the pasture.This vegetation settled and became dominant in the stations with a more advanced soil type and a higher water content than that in the Apennines.Therefore, it was less xeric, and had a more closed structure compared to the xeric meadows, highlighting the effect of the tree cover.

Analysis of Floristic Biodiversity
The most consistent group of species, by both qualitative and quantitative analysis, taking into account coverage, was arid secondary grazed meadows (Figure 6).The pool of diagnostic species of the Apennine Phleo ambigui-Bromion erecti alliance was highlighted, among which the graminaceous Brachypodium rupestre was relevant, forming a fairly uniform mat.It was accompanied by Cytisus spinescens, Asperula purpurea subsp.purpurea, and Galium lucidum subsp.lucidum.The group of plants belonging to the superior categories of the suborder Artemisio albae-Bromenalia, order Phleo ambigui-Brometalia erecti, and class Festuco-Brometea was very abundant.Among these, the species forming the structure of secondary grasslands in the Apennines was Bromus erectus.Even if this vegetation was attributable to the Phleo ambigui-Bromion erecti alliance, the difference was the presence in the pinewood of a community dominated by Brachypodium rupestre, which was absent in the pasture.This vegetation settled and became dominant in the stations with a more advanced soil type and a higher water content than that in the Apennines.Therefore, it was less xeric, and had a more closed structure compared to the xeric meadows, highlighting the effect of the tree cover.The results show a constant increase in species richness after cutting (Table 6).This increase does not correspond to an increase in the values of the biodiversity indices, which showed a sharp decline after cutting.However, these values reached higher values after two years before pre-cutting.
A decrease in the values of floristic biodiversity was more evident in the strips harvested by horse, as shown in Figure 7. On the contrary, in the strips where the collection of material was carried out by cable yarder, there was no decrease in Shannon index values (Table 6), but even after the second year, the greatest increases occurred for both indices.The results show a constant increase in species richness after cutting (Table 6).This increase does not correspond to an increase in the values of the biodiversity indices, which showed a sharp decline after cutting.However, these values reached higher values after two years before pre-cutting.
A decrease in the values of floristic biodiversity was more evident in the strips harvested by horse, as shown in Figure 7. On the contrary, in the strips where the collection of material was carried out by cable yarder, there was no decrease in Shannon index values (Table 6), but even after the second year, the greatest increases occurred for both indices.

Physical and Chemical Analyses of Soil
Soil texture was characterized by a high content of rocky skeleton for both blocks (48.6 ± 17.3%), while sand, silt, and clay content were statistically different between the two blocks (p < 0.05).In block 1, sand was 59.2 ± 1.4%, silt was 33.7 ± 1.3%, and the clay was 7.0 ± 1.1%.Thus, the soil can be assigned to the sand loam (SaLo) soil material class.In block 2, the sand was 36.6 ± 3.7%, the silt was 59.4 ± 3.6%, and the clay was 4.0 ± 1.3%.Thus, the soil can be assigned to the silt loam (SiLo) soil material class.Relative to the granulometric values and the soil water method (K.Saxton), the soil field capacity (CC) was calculated at 24% and 27%, respectively.
For each year, sampling was done once in spring, once in summer, and once in autumn.During the sampling periods (2015, 2016, and 2017) soil moisture showed significant differences (p < 0.05) between spring and autumn (average moisture 31 ± 5%) and summer (average moisture 22 ± 3%).
Soil BD showed significant differences only among treatments and years (Table 7).In particular, in the strips harvested by horse, BD was higher than the control and the strips harvested by winch and cable yarder.Nevertheless, in 2017, a recovery trend was highlighted.In the strips harvested by winch, the BD was similar to that of strips harvested by cable yarder, but in 2017, no clear recovery trend was shown.In the strips harvested by cable yarder, BD was higher than that of the control, but lower than that of the strips harvested by horse and winch, although in 2017, a complete recovery was shown.

Physical and Chemical Analyses of Soil
Soil texture was characterized by a high content of rocky skeleton for both blocks (48.6 ± 17.3%), while sand, silt, and clay content were statistically different between the two blocks (p < 0.05).In block 1, sand was 59.2 ± 1.4%, silt was 33.7 ± 1.3%, and the clay was 7.0 ± 1.1%.Thus, the soil can be assigned to the sand loam (SaLo) soil material class.In block 2, the sand was 36.6 ± 3.7%, the silt was 59.4 ± 3.6%, and the clay was 4.0 ± 1.3%.Thus, the soil can be assigned to the silt loam (SiLo) soil material class.Relative to the granulometric values and the soil water method (K.Saxton), the soil field capacity (CC) was calculated at 24% and 27%, respectively.
For each year, sampling was done once in spring, once in summer, and once in autumn.During the sampling periods (2015, 2016, and 2017) soil moisture showed significant differences (p < 0.05) between spring and autumn (average moisture 31 ± 5%) and summer (average moisture 22 ± 3%).
Soil BD showed significant differences only among treatments and years (Table 7).In particular, in the strips harvested by horse, BD was higher than the control and the strips harvested by winch and cable yarder.Nevertheless, in 2017, a recovery trend was highlighted.In the strips harvested by winch, the BD was similar to that of strips harvested by cable yarder, but in 2017, no clear recovery trend was shown.In the strips harvested by cable yarder, BD was higher than that of the control, but lower than that of the strips harvested by horse and winch, although in 2017, a complete recovery was shown.Soil PR showed statistically significant differences only among treatments and years (Table 8).In particular, in the strips harvested by horse, PR was higher than that of the control as well as the strips harvested by winch and cable yarder.Nevertheless, in 2017, no clear recovery trend was shown.In the strips harvested by winch, the PR was similar to that of the control from 2015.In the strips harvested by cable yarder, in 2015, the PR was higher than that of the control and strips harvested by winch, but lower than that of the strips harvested by horse, although in 2017, a complete recovery was shown.
Table 8. Results of the ANOVA and Tukey test for soil penetration resistance during the three years after felling (df 3, 144; average ± SD), difference tested between the soil of the harvested area, and control.The soil SR showed statistically significant differences only among treatments and years (Table 9).In particular, in the strips harvested by horse, in 2015 and 2016, SR was higher than that of the control as well as the strips harvested by winch and cable yarder.Nevertheless, in 2017, a complete recovery was shown.In the strips harvested by winch and cable yarder, in 2015, SR was higher than that of the control, although already in 2016, a complete recovery was shown.

Treatments
Table 9. Results of the ANOVA and Tukey test for soil shear resistance during the three years after felling (df 3, 144; average ± SD), difference tested between soil of harvested area, and control.The pH is a very important parameter for the correct functioning of soil, and its variations influence various soil parameters and processes [60].The pH values showed statistically significant differences between blocks, treatments, and years.However, this parameter did not seem to be clearly influenced by silvicultural treatment or logging operations (Table 10).The pH variations were recorded only for 2017.

Treatments
Soil organic matter content (OM) showed significant differences only among treatments and years (Table 11).In particular, in the strips harvested by horse, in 2015, OM was lower than that of the control as well as the strips harvested by winch and cable yarder, and the trends in 2016 and 2017 were negative.In 2015, in the strips harvested by winch and cable yarder, the OM was lower than the control, but for the strips harvested by winch, the trends in 2016 and 2017 were negative.There was a complete recovery for the strips harvested by cable yarder in 2017.Soil inorganic carbon content (IC) analysis was done only in 2017.The results showed significant differences only among the treatments and control (Table 12).In particular, in the strips harvested by horse, the IC was lower than that of strips harvested by cable yarder, but it was higher than that of the strips harvested by winch, and similar to the control.Different letters show significant differences among values in a column (Tukey test).

Soil Biodiversity Analysis
The QBS-ar index (Table 13) showed significant differences only among treatments and years.In particular, in the strips harvested by horse and by winch, in 2015, QBS-ar was higher than that of the strips harvested by cable yarder, and lower than that of the control.In 2016, a strong decrease was recorded for the strips harvested by horse, and a slightly decrease was recorded for those harvested by winch compared to the control.Nevertheless, both showed complete recovery in 2017.For the strips harvested by cable yarder, in 2015 and 2016, QBS-ar was unchanged and lower than the control, but in 2017, its value appreciably increased.Microarthropod density (Table 14) showed statistically significant differences among treatments and years.In particular, in the harvested strips, the density values were always lower than the control, but the trends were positive from 2015 to 2017.In 2017, density that was more similar to the control was shown in the strips harvested by winch.More time is necessary for the recovery of the strips harvested by horse and by cable yarder.

Discussion
The management of pine plantations has strong and variable effects on plant species occurrence and diversity due to plantation and treatment operations, and the alteration of ecological processes caused by changes in the landscape and stand structure [61,62].Plantations contribute to biodiversity conservation in various ways.In Poorbabaei [63], a high similarity in species composition between plantation and the adjacent natural forest, which is the main source of seed in plantations, was considered.Neighboring plantation and natural forest has resulted in the dispersion of hardwood tree seeds within the plantation.
The results highlight that the good density and richness of tree species in this pine plantation indicate the high potential reached by the stand for biodiversity restoration.Similar results are shown in studies where plantations of native and/or exotic tree species increased biodiversity by promoting woody understory regeneration [64,65].The silvicultural treatment applied in this research showed positive effects on density, richness, and biodiversity of tree species in three years.Logging operations have important effects on regeneration and tree biodiversity, which are the main topics for restoration and renaturalization purposes.
From the results concerning stand regeneration, different taxonomic compositions of the tree forest community among the logging treatments are shown, in particular in the percentage of distribution.A clear simplification with respect to the others is shown in the strips harvested by horse, with the presence of about 45% Fraxinus ornus.The general finding shows that ground-based logging systems allowed for the presence of Robinia pseudoacacia and only marginal Pinus nigra regeneration.In general, felling produced an abundance of light-demanding species due to the increase of solar radiation [48].
In terms of regeneration density, silvicultural treatment showed a positive effect with greater consistency in the strips harvested by winch and cable yarder.When compared to the control, they had increases of 85% and 72%, respectively.The strips harvested by horse showed a positive trend, with an increase of about 69% compared to the control.
In this study, attention was paid to other important aspects, such as the richness and diversity of tree species.In particular, silvicultural treatment more positively influenced the richness (an evenness index increase ranging from 5.2% to 14.3%) and marginally influenced diversity.The strips harvested by winch and by cable yarder had higher richness values than the control and the strips harvested by horse.
Floristic biodiversity, which consists of the numerically and structurally consistent presence of species referable to herbaceous vegetation, and is also present in other pine woods of the Aterno Valley [55], suggested a slow evolution or in some cases no evolution in the undergrowth of artificial pine forests.However, one sign of an ongoing evolutionary process can be found in the analysis of species packages of the shrub layer with the presence of entities referable to mantle vegetation, which are shrub-like formations that set up at the edge of the forest or colonize fields and abandoned pastures.These formations had particular importance in the analysis of the dynamic relationships between the various types of Apennine vegetation.This is confirmed by the significant presence of Cytisophyllum sessilifolium, Rosa canina, Amelanchier ovalis subsp.ovalis, Juniperus communis, Juniperus oxycedrus subsp.deltoides, Viburnum lantana, and Prunus spinosa.
In accordance with results of other similar interventions [66,67] in the years following cutting, we observe an increase in the floristic richness linked to the opening of the tree layer.Increased species richness following canopy disturbance is largely attributed to early seral, shade-intolerant herbs and shrubs invading sites to take advantage of increased light conditions [68,69].
It is important to note that the data presented so far concern a limited period of time that does not allow for more structured statistical analyses.More time is needed to evaluate whether the cutting effect on biodiversity will last long, and observe whether non-forest species are able to regenerate when the crown cover tends to close [70,71].The future effects on biodiversity of each treatment will also require more time to be assessed.
The findings show that soil BD, PR, and SR were influenced by both silvicultural treatment and logging operations.In particular, silvicultural treatment had an impact mainly in the first year post-harvesting, and after two and three years post-harvesting.Recovery was shown, as also found in other studies [44,48].Among the strips harvested by winch and by cable yarder, only a few differences were found, and the values after three years were similar to the control.The strips harvested by horse showed an important impact in the first year, and the recovery was much slower.Similar values and situations have been observed in other studies, where logs were skidded and vehicles or animals were moved across the forest soil [9,31].
From the data analysis concerning pH, OM, and IC content, it was possible to note different behaviors among the variables.The pH values, as found in other studies [44,48,72], did not show any statistical relation to treatments, periods, or blocks.OM content was affected by both silvicultural treatment and logging operations.In particular, silvicultural treatment had an impact mainly in the first year post-harvesting only for the strips harvested by cable yarder, and recovery was reached after two years post-harvesting.For the strips harvested by horse and by winch (ground-based extraction), a negative effect was highlighted in the first year post-harvesting.The trends in the following years were negative, which was mainly attributable to the partial remixing of the superficial soil layers.These negative effects in the first year were due to the silvicultural treatment and the logging operations, while in the following years, the ground-based extraction systems negatively affected OM content, due to partial litter removal [73].IC content seemed to be affected mainly by the logging operations, and ground-based systems in particular.In the short term, the ground-based extraction systems that cause the mixing of the litter with the soil superficial layers (horse and mainly winch) seem to have a negative effect on the IC content.The loss of tree cover, on the contrary, where no litter mixing (cable yarder) is observed, seems instead to increase the IC content in the short term.
The QBS-ar values show a significant positive correlation with soil physical parameters according to the literature [44,74].The observed variation is explained by the different degrees of soil compaction and the abundance of litter associated with sudden stand removal.This mainly affected the strips harvested by cable yarder.In addition to the QBS-ar index, population density was evaluated during sampling in terms of individuals per dm 2 .As can be observed from the data gathered, microarthropod density was inferior in all of the areas involved in the impact caused by the extraction activities.There was a significant difference between the logging operations, for both the QBS-ar index and population density.The strips harvested by winch had higher values than those harvested by horse and cable yarder.In particular, the values recorded in the strips extracted by cable yarder were lower with respect to the control and the other extraction systems.This is mainly due to the concomitant effect of removing the tree cover and maintaining the pine litter in its original condition.In this study, soil compaction and uncovering led to the rarefaction of specialized groups such as Protura, Diplura, and Pauropoda, and between the logging operations, no particular difference is shown.However, regarding microarthropod density, this index did not efficiently describe the conditions in terms of soil biodiversity.The values that it assumes should be interpreted in terms of trend over time.
Two principal component analyses (PCAs) were carried out to investigate any linear correlations between the expressed soil conditions of the main six studied traits, and between the tree regeneration situations of the main four studied traits.Data corresponding to each independent variable were standardized using Box-Cox lambda to minimize the scaling effect due to the different measurement units.For the soil conditions, the principal components, PC1 and PC2, explained 68% and 22% of the total variance, respectively.The PC1 and PC2 scores for the three logging operations and the control are shown in Figure 8.Each logging operation condition is distinct from the other two and from the control based on the score plot, as depicted in Figure 8.In general, a positive trend is shown from 2015 to 2017.In particular, strips harvested by winch and cable yarder showed higher recovery, and for some parameters, the recovery was complete.In detail, the strips harvested by cable yarder had a higher recovery in their physical soil parameters, which implies that these parameters were closely associated with their PC1 scores, according to Marchi [75].Similarly, the ground-based extraction systems (horse and winch) had a higher recovery of the biological soil parameters that were closely associated with their PC2 scores.
For the tree regeneration situations, the principal components, PC1 and PC2, explained 66% and 30% of the total variance, respectively.The PC1 and PC2 scores for the three different logging operations and the control are shown in Figure 9.The tree regeneration situation was similar among the strips harvested by cable yarder and winch.The situation for the strips harvested by horse was similar to the control in 2016, and showed an implementation in 2017 that approached the values of the other two logging operations based on the score plot (Figure 9).In general, a positive trend can be associated with stand harvesting.In detail, the strips harvested by cable yarder and winch had a higher score for tree regeneration density and in the evenness index.These parameters are closely associated with the PC2 axis.SIV values are always associated with the PC2 axis, but with an inverse relation.The Shannon index is closely associated with PC1 scores, and it only shows the best situation for the strips harvested by winch.
Properly managed and programmed clear-cutting within silvicultural systems, such as the strip clear-cutting described here, can positively influence sustainability.This work indicates that less invasive extraction systems (i.e., cable yarder, forest winch) also inflict the least damage on forest soils.For the tree regeneration situations, the principal components, PC1 and PC2, explained 66% and 30% of the total variance, respectively.The PC1 and PC2 scores for the three different logging operations and the control are shown in Figure 9.The tree regeneration situation was similar among the strips harvested by cable yarder and winch.The situation for the strips harvested by horse was similar to the control in 2016, and showed an implementation in 2017 that approached the values of the other two logging operations based on the score plot (Figure 9).In general, a positive trend can be associated with stand harvesting.In detail, the strips harvested by cable yarder and winch had a higher score for tree regeneration density and in the evenness index.These parameters are closely associated with the PC2 axis.SIV values are always associated with the PC2 axis, but with an inverse relation.The Shannon index is closely associated with PC1 scores, and it only shows the best situation for the strips harvested by winch.
Properly managed and programmed clear-cutting within silvicultural systems, such as the strip clear-cutting described here, can positively influence sustainability.This work indicates that less invasive extraction systems (i.e., cable yarder, forest winch) also inflict the least damage on forest soils.

Conclusions
Renaturalization treatments in black pine afforestation is an important topic that should be considered.There is a need to favor the evolution of artificial pine forests toward natural forest systems.The original main purpose of these plantations was to maintain and improve soil characteristics to encourage more complex forest systems.This shift from even-aged pine monocultures to a mixed forest type through natural regeneration is the final goal of these plantations.In Italy in recent years, this substitution has been planned by means of thinning, whereas clear-cutting has been used less.However, the limit of 60-65 years represents the threshold beyond which the reaction capacity of the stand to thinning is greatly reduced, especially in low-fertility areas.
On the whole, this study focused on pine forests and tried to suggest one typology of clear-cutting (dismantling cutting) on strips associated with different extraction management techniques.Some ecological and environmental aspects associated with renaturalization treatments, including techniques applied by different mechanizations in black pine afforestation, have been highlighted as well as how the renaturalization and the active ecological management of these stands could affect soil and vegetation.An outline of answers to the main research questions is as follows:

•
In terms of immediate impact on soil features, the less invasive extraction systems seem to be cable yarder and forest winch.

•
A clear soil recovery trend with good capability was visible, in particular for the two extraction systems by cable.However, over a three-year period, only a partial but substantial recovery was shown.The main physical soil characteristics reached quick recovery in particular for the strips harvested by cable yarder and forest winch; the organic matter content highlighted a heavy impact, but only for the strips harvested by horse and forest winch.The biological soil characteristics, in particular arthropod communities, still showed a clear impact, which was more highlighted in the strips harvested by cable yarder • Silvicultural treatment, particularly logging activities, affected natural tree regeneration in a different way.In general, silvicultural treatment showed qualitative and quantitative improvement in terms of tree regeneration.In particular, the extraction systems by forest winch and cable yarder showed better situations • Silvicultural treatment seems not to have led to improvements at the level of the herbaceous and shrubby layers.However, clear differences are shown among the different harvesting systems.
In particular, the strips harvested by cable yarder showed clear recovery, in terms of positive compositions and quantities, whereas the situation was slightly worse for the strips harvested by forest winch, and definitely worse for the strips harvested by horse.
It is important to highlight that the data concern a limited period of time that does not allow for more structured statistical analyses.More time is needed to evaluate whether cutting has more effects.This first step of the research was planned in order to obtain an overview in terms of the environmental impacts related to multifunctional approaches to the forest management of black pine afforestation.

Figure 1 .
Figure 1.Geographical information and detailed scheme of the treatment area.

Figure 1 .
Figure 1.Geographical information and detailed scheme of the treatment area.

Figure 2 .
Figure 2. Scheme of the systematic sampling method applied in every harvested strip.The black points are only the starting points of every linear transect (grey lines).

Figure 2 .
Figure 2. Scheme of the systematic sampling method applied in every harvested strip.The black points are only the starting points of every linear transect (grey lines).

Figure 3 .
Figure 3. Tree species percentage distribution, average values grouped for logging treatments and control.

Figure 4 .
Figure 4. Tree regeneration density, average values with standard deviation for the different treatments, one year after harvesting (2016) and two years after harvesting (2017).ANOVA results shown statistically significant differences, p < 0.05 and the Tukey test applied found five groups, showed by different letters.

Figure 3 .
Figure 3. Tree species percentage distribution, average values grouped for logging treatments and control.

Figure 3 .
Figure 3. Tree species percentage distribution, average values grouped for logging treatments and control.

Figure 4 .
Figure 4. Tree regeneration density, average values with standard deviation for the different treatments, one year after harvesting (2016) and two years after harvesting (2017).ANOVA results shown statistically significant differences, p < 0.05 and the Tukey test applied found five groups, showed by different letters.

Figure 4 .
Figure 4. Tree regeneration density, average values with standard deviation for the different treatments, one year after harvesting (2016) and two years after harvesting (2017).ANOVA results shown statistically significant differences, p < 0.05 and the Tukey test applied found five groups, showed by different letters.

Figure 5 .
Figure 5. SIV detailed analysis between the treatments and the control for the year 2017.Figure 5. SIV detailed analysis between the treatments and the control for the year 2017.

Figure 5 .
Figure 5. SIV detailed analysis between the treatments and the control for the year 2017.Figure 5. SIV detailed analysis between the treatments and the control for the year 2017.

Figure 6 .
Figure 6.Phytosociological surveys carried out in the afforested area and analysis of the groups of species related to the phytosociological syntaxa.

Figure 6 .
Figure 6.Phytosociological surveys carried out in the afforested area and analysis of the groups of species related to the phytosociological syntaxa.

Figure 7 .
Figure 7. Floristic evenness index trend related to the different harvesting systems.

Figure 7 .
Figure 7. Floristic evenness index trend related to the different harvesting systems.

Figure 8 .
Figure 8. Score plot of the components PC1 and PC2 of the principal component analysis referred to soil characteristics.C (orange): strips harvested by cable yarder; W (pink): strips harvested by winch; H (violet): strips harvested by horse; Control (green): stand not harvested.The numbers 16 and 17 represent the years 2016 and 2017, respectively.

Figure 8 .
Figure 8. Score plot of the components PC1 and PC2 of the principal component analysis referred to soil characteristics.C (orange): strips harvested by cable yarder; W (pink): strips harvested by winch; H (violet): strips harvested by horse; Control (green): stand not harvested.The numbers 16 and 17 represent the years 2016 and 2017, respectively.

Figure 9 .
Figure 9. Score plot of the components PC1 and PC2 of the principal component analysis referred to tree regeneration situations.C (orange): strips harvested by cable yarder; W (pink): strips harvested by winch; H (violet): strips harvested by horse; Control (green): stand not harvested.The numbers 16 and 17 represent the years 2016 and 2017 respectively.

Figure 9 .
Figure 9. Score plot of the components PC1 and PC2 of the principal component analysis referred to tree regeneration situations.C (orange): strips harvested by cable yarder; W (pink): strips harvested by winch; H (violet): strips harvested by horse; Control (green): stand not harvested.The numbers 16 and 17 represent the years 2016 and 2017 respectively.

Table 1 .
Detailed information regarding temperature and precipitation values throughout the years of the study period.

Table 2 .
Different extraction methods applied, characterized by different machines and techniques.
H: animal skidding; W: winching downhill by a forestry-fitted farm tractor using snatch block; C: yarding downhill by a cable yarder with automatic carriage.
Forests 2018, 9, x FOR PEER REVIEW 6 of 24 every replicate (block).The different extraction methods were characterized by different machines and techniques (detailed in Table 2): (a) animal skidding downhill by TPR (heavy rapid skidding) horse, (b) winching downhill by a forestry-fitted farm tractor using snatch block, and (c) yarding downhill by a cable yarder with automatic carriage.The two blocks (replicates) differed in altitude, soil texture, and dendrometric characteristics.

Table 2 .
Different extraction methods applied, characterized by different machines and techniques.
H: animal skidding; W: winching downhill by a forestry-fitted farm tractor using snatch block; C: yarding downhill by a cable yarder with automatic carriage.

Table 3 .
Detailed formulas and descriptions of the Species Importance Value (SIV) index, Regeneration species diversity index and structural evenness index.
(Number of plots containing a species × 100)/sum of frequencies of all species Relative dominance (RDo) (Sum of the height of a species × 100)/sum of total height of all species Regeneration species diversity index (Shannon diversity index) H = −Σ(ni/n)•Ln (ni/n) ni = denotes the SIV of a species and n = denotes the sum of the total SIV of all the species Structural evenness index J = H /Hmax H is the number derived from the Shannon diversity index and Hmax is the maximum value of H
cDifferent letters show significant differences among values in a column (Tukey test).
cDifferent letters show significant differences among values in a column (Tukey test).

Table 6 .
Floristic species diversity indexes for the studied areas.

Table 6 .
Floristic species diversity indexes for the studied areas.

Table 7 .
Results of the ANOVA and Tukey test for soil bulk density (BD) during the three years after felling (df 3, 144; average ± SD), difference tested between the soil of the harvested area, and control.
Different letters show significant differences among values in a column (Tukey test).

Table 7 .
Results of the ANOVA and Tukey test for soil bulk density (BD) during the three years after felling (df 3, 144; average ± SD), difference tested between the soil of the harvested area, and control.
Different letters show significant differences among values in a column (Tukey test).

Table 10 .
Results of the ANOVA and Tukey test for pH (df 3, 72; average ± SD), difference tested between the soil of the harvested area, and control.

Table 11 .
Results of the ANOVA and Tukey test for organic matter (df 3, 72; average ± SD), difference tested between the soil of the harvested area, and control.
Different letters show significant differences among values in a column (Tukey test).

Table 12 .
Results of the ANOVA and Tukey test for inorganic carbon (df 3, 24; average ± SD), difference tested between soil of harvested area and control.

Table 13 .
Results of the Kruskal-Wallis and Tukey test for Soil Biological Quality (QBS-ar) index data (df 3, 144; median), difference tested between disturbed, undisturbed, and control soil.

Table 14 .
Results of the Kruskal-Wallis and Tukey test for soil microarthropod density data (df 3, 144; median), difference tested between disturbed, undisturbed, and control soil.