Next Article in Journal
Identification of Low-Light-Resistant Germplasm and Related Loci of Soybean
Next Article in Special Issue
Interactive Effects of Tillage Systems and Nitrogen Fertilizer Rates on the Performance of Mustard-Boro-aman Rice Cropping Systems under Conservation Agriculture Practices
Previous Article in Journal
Carbon Budget of Paddy Fields after Implementing Water-Saving Irrigation in Northeast China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Soil Penetration Resistance Influenced by Eucalypt Straw Management under Mechanized Harvesting

by
Judyson de M. Oliveira
1,2,*,
Cássio A. Tormena
2,
Gérson R. dos Santos
3,
Lincoln Zotarelli
1,
Raphael B. A. Fernandes
4 and
Teógenes S. de Oliveira
4
1
Horticulture Science Department, University of Florida/IFAS, 2550 Hull Rd, Gainesville, FL 32611, USA
2
Agronomy Department, State University of Maringá, Av. Colombo, n 5790, Maringá 87020-900, PR, Brazil
3
Department of Statistics, Federal University of Viçosa, Av. Peter Henry Rolfs, s/n, Viçosa 36570-900, MG, Brazil
4
Soil Science Department, Federal University of Viçosa, Av. Peter Henry Rolfs, s/n, Viçosa 36570-900, MG, Brazil
*
Author to whom correspondence should be addressed.
Agronomy 2022, 12(6), 1482; https://doi.org/10.3390/agronomy12061482
Submission received: 6 May 2022 / Revised: 10 June 2022 / Accepted: 18 June 2022 / Published: 20 June 2022

Abstract

:
This study aimed to evaluate the impacts of mechanized harvesting and soil tillage on soil penetration resistance (PR), influenced by the eucalypt straw management under sandy clay Oxisol in Southern Brazil. The study was conducted in a eucalyptus production area under Oxisol in Paraná State, Brazil. The treatments consisted of two harvesting systems: harvester + forwarder (HF) and feller + skidder (FS) both applied in areas under coppicing and stand renewal eucalypt cultivation systems. For stand renewal areas, eucalypt straw was managed on the soil surface at levels of 100, 50, and 0% before soil tillage. PR and soil moisture measurements were made in points distributed in regular grid for all treatments. This grid also was used to evaluate the geospatial behavior of PR in the stand renewal areas. During the measuring of PR, the averages (± confidence interval) of soil moisture up to 0–0.60 m depth were 0.20 ± 0.01 and 0.24 ± 0.01 in coppicing and stand renewal areas, respectively. In areas under coppicing, the PR mean ± confidence interval at 0–0.05 m layer in HF (1.28 ± 0.24 MPa) was lower than in FS treatment (2.11 ± 0.44 MPa). However, the PR values were similar between treatments in stand renewal areas, regardless of the forest straw level on the soil surface. For both harvesting systems, there was a lack of spatial dependence of PR up to 0.40 m soil depth, indicating some physical and mechanical homogenization induced by the soil tillage in the layer. Eucalypt straw contributed to mitigating effects of harvest traffic on PR level in coppicing forest systems. However, different levels of eucalypt straw managed before soil tillage did not influence PR levels in stand renewal forest systems.

1. Introduction

Planted forests occupy 9.93 million hectares, approximately 1.2% of the Brazilian territory, and are responsible for 91% of the wood for the industry. Eucalyptus plantation has stood out for its yield average of 36.8 m3 ha−1 year−1 in 2021 [1], with a harvesting cycle lasting less than seven years, which has intensified the mechanization use in soils under forest production [2]. However, the practices used for harvesting and transporting operations frequently do not consider the potential risk of soil compaction [3] and its distribution in depth [4]. Monitoring the impact of the mechanized operations on soil physical quality could reduce the risk of its compaction and enhance planted forest production and environmental protection.
In planted forest production, soil physical quality degradation can be intensified by the mechanical cutting of trees and dragging and lifting logs, requiring machines with large masses and sizes [5]. In addition, the heterogeneous pattern of the machine traffic across the field increases the physical disturbance in the soil surface, which can also affect deeper soil layers [4,6]. Thus, these changes in soil physical properties increased bulk density and resistance to root penetration and decreased macroporosity, water infiltration, and airflow, resulting in non-uniform stands and lower yields of forest species [7,8].
Mainly in stand renewal areas, soil tillage is carried out to mitigate compaction and create a suitable physical environment for root development [9,10], succeeding the subsequent management practices. In areas with high forest straw production, it is necessary to remove part or total of the forest straw to facilitate the displacement of machinery equipment and, consequently, increase their operational performance [11]. Although the deposition of forest straw on the soil surface contributes to attenuating the impacts of machinery traffic, cumulative effects of removal practices before soil tillage are not well known yet [2,12,13].
Many soil physical properties such as those mentioned above have been used to characterize and quantify the effects of forest management practices. Between them, soil penetration resistance (PR) is widely used as an indicator of soil compaction due to the relatively easy measurement in the field and being correlated with dendrometric variables (e.g., wood volume, diameter at breast height, total height, etc.). Therefore, there is a critical PR value, above which roots growth is strongly limited, and can be used as an indicator of soil physical quality under forest management systems. PR > 3 MPa have been reported as restrictive for the growth of roots of Pinus sp [14] and may reach 4.5 MPa [15] for eucalyptus. PR values between 2 and 3 MPa have been considered limiting for other forest species (e.g., Alnus cordata, Betula pendula, and Larix kaempferi) [16]. Such variations occurred because of the differences between species and edaphic conditions [17,18]. However, it is known that root growth rate decreases exponentially with the increment in PR, and root development can be affected, depending on the magnitude and intensity with which plants are subjected to the effects of the resistance [19].
Although the RP measurement is relatively straightforward, the soil moisture variation across the areas can negatively affect its results due to an inverse relationship between these variables [20,21]. However, PR measurement should be performed under homogenous soil moisture conditions, preferentially close to the soil field capacity [21,22]. Additionally, the geospatial evaluation of PR can allow the identification of compacted sites for specific management operations [23,24]. Such detailed PR information can improve soil physical quality evaluation in mechanized planted forest areas [24,25,26].
In this study, we evaluate the impact of mechanized operations (harvesting and soil tillage) on soil compaction measured by PR in eucalypt plantation harvested using two different methods (harvester + forwarder and feller + skidder), both applied to coppicing and stand renewal with eucalypt straw managed on the soil surface. We hypothesized that eucalypt straw permanence on the soil surface can reduce the effect of the harvesting traffic and tillage on soil compaction measured by PR. The objective was to evaluate the impacts of the mechanized harvesting and soil tillage on PR, influenced by the eucalypt straw management under sandy clay Oxisol in Southern Brazil.

2. Materials and Methods

2.1. Experimental Area and Treatments

The experimental area under eucalyptus plantation was located at Monte Alegre Farm (Block BRM-J9A), belonging to the Klabin Company, in the Telêmaco Borba county, Paraná, Southern Brazil (24°18′15′’ S, 50°36′39′’ W, ~770 m above sea level) (Figure 1a). According to Köppen’s classification, the climate of the region is Cfa and Cfb, a subtropical climate characterized by hot summers, infrequent frost, and a trend of rainfall concentration in the summer months. The average annual rainfall is between 1400 and 1600 mm, respectively [27]. The dominant soil is classified as an Oxisol, according to [28]. The relief is classified as gently undulating. The soil chemical and particle-size distribution characterization for the experimental area are shown in Table 1.
The experimental area has had a history of use with commercial stands of forests for approximately 30 years, with eucalyptus plantation with a cycle of 20 years and modified to eight-year cycles from September 2008. The last forest stand had been harvested in November 2016 and was composed of hybrids that originated from the cross between Eucalyptus grandis x Eucalyptus urophylla, clone 12/52. The spacing was 2.50 m between plants and 3.50 m between rows, totaling 1142 trees ha−1. This stand had the following indicators: average canopy height of 19.50 m, average diameter at breast height of 0.16 m, and average annual increment at seven years of 65 m3 ha−1 year−1. Then, on March 16, 2017, a new planting was carried out using eucalyptus seedlings, clone I-144, at a spacing of 1.82 m between plants and 3.30 m between rows, totaling 1666 trees ha−1.
The experiment was set up to represent two eucalypt harvesting systems: harvester + forwarder (HF) and feller + skidder (FS) applied to coppicing and stand renewal cultivation systems with straw management on the soil surface (Figure 1a). In the coppicing area, the soil tillage was not used and the level of eucalypt straw kept on the soil surface was the same as after harvesting. In the HF, there was an accumulation of eucalypt straw (bark, branches, leaves, tips, etc.), while for FS there was no accumulation of eucalypt straw on the soil surface. For stand renewal areas, the level of eucalypt straw was managed before the soil tillage considering differences in relation to the amount of straw, establishing the maintenance of approximately 100, 50, and 0% of forest straw on the soil surface (Figure 1b). For HF, 100% represents the eucalypt straw produced at harvest that was kept along with the remaining forest litter; 50% represents partial removal of straw produced at harvest along with the remaining forest litter; 0% refers to the total removal of eucalypt straw but keeping the remaining forest litter. In the FS, 100% refers to the eucalypt straw produced outside the plot (tips), which were returned to the area and spread on the remaining forest litter; 50% corresponded to the permanence of the remaining forest litter, and 0% represented the total removal of the forest litter (Figure 1b). After eucalypt straw management (removal or addition), the soil tillage was carried out by subsoiling down to 0.50 m depth using a bulldozer equipped with a subsoiler.
Thus, the treatments established consisted of two harvesting systems: harvester + forwarder and feller + skidder, both applied in coppicing and stand renewal areas with eucalypt straw level (100, 50, and 0%) managed before soil tillage, as described in Figure 1b.

2.2. Soil Sampling

The sampling for soil penetration resistance (PR) measurements was performed between July 15 and 30, 2017, approximately 240 days after harvest and 120 days after soil tillage. Measurements of PR were taken by adopting a systematic sampling scheme with readings in quadruplicates of 22 sampling points defined by the vertices of isosceles triangles (15 m of height and 30 m of base) demarcated in the experimental area (Figure 2). Each triangle had an area of 225 m2, and the sampling encompassed planting rows and between rows, alternately. In the areas of stand renewal, a grid of 66 points formed by the union of the 22 points of each one of the three straw levels (100, 50, and 0%), in both harvesting systems (HF and FS), was used to study the geospatial behavior of PR (Figure 1a).
Penetration resistance measurements were taken up to 0.60 m depth using a digital electronic penetrometer (penetroLOG/PLG1020) with a constant speed of approximately 30 mm s−1. In 50% of the PR sampling points, alternately selected, disturbed samples were taken for the determination of soil moisture, chemical, and particle-size distribution characterization. The samples were collected from the center of the 0–0.20, 0.20–0.40, and 0.40–0.60 m soil depth layers, using a Dutch auger. The gravimetric soil water content was determined by drying disturbed samples in an oven at 105 °C for 48 h, according to [29]. The sampling points were georeferenced with a Leica GS14 GNSS receiver, set in the metric coordinate system (UTM), Zone 22s, and Datum Sirgas2000.

2.3. Statistical Analyses

All statistical analyzes were performed using R software [30]. The confidence intervals (CI) to means were calculated for the variables with a confidence level of 90% [31]. The coefficient of variation (CV%) was calculated and classified as: low (CV < 15%), medium (15 ≥ CV ≤ 50%), and high (CV > 50%) [32]. Data distribution normality was also tested using the Shapiro-Wilk test (p > 0.05). The correlation between PR and soil moisture was verified by the Student’s t-test at a 95% confidence level using PR averages every 0.20 m in the profile.
The PR average of each 0.20 m soil depth layer was also used for the geospatial evaluation in the stand renewal areas. The data distribution asymmetry was calculated and classified according to the coefficient of asymmetry (As): asymmetric (|As| < 0.15), moderate (0.15 < |As| < 1.0), and strong (|As| > 1.0). As for the shape of the distribution curve (tailedness), it was classified according to the coefficient of kurtosis (Ks): mesokurtic (Ks = 0), platykurtic (Ks< 0), and leptokurtic (Ks > 0) [33]. Geostatistical analysis was performed using the R functions available in the gstat package. The semivariances were calculated according to the equation [31]: ɤ h = 1 / 2 N h i = 1 ) N h Z x i + h Z x i 2 ; where ɤ(h) is the semivariance for the interval class h; N(h) is the number of experimental pairs of observation; and Z(xi) and Z (xi + h) are the values of attributes measured at the positions xi and xi + h, separated by a vector h (distance between samples). From the semivariances analysis, the following parameters were obtained: range (Ao), nugget effect (Co), and sill (C), which were used to evaluate the degree of spatial dependence of the variable, by calculating the index of spatial dependence (SD), according to the expression: SD = [Co/(C + Co)] * 100, being classified as strong (SD ≤ 25%), moderate (25% < SD ≤ 75%), and weak (SD > 75%) [34].

3. Results and Discussion

3.1. Soil Moisture Conditions for Penetration Resistance Sampling

In coppicing areas, during the measure of PR the soil moisture averages (± confidence interval) were 0.21 ± 0.02 and 0.20 ± 0.01 kg kg−1 in the soil profile (0–0.60 m depth layer), for harvester + forwarder (FS) and feller + skidder (FS) treatments, respectively; being statistically similar (Table 2). The averages of soil moisture (0–0.60 m depth layer) were also similar between the harvest treatments in the stand renewal areas, being equal to 0.23 ± 0.01 and 0.25 ± 0.01 kg kg−1, for HF and FS treatments, respectively (Table 2). The CVs were classified as medium and low (15.92 and 9.60%) in the areas of coppicing and stand renewal, respectively (Table 2). Similar moisture conditions between treatments were due to the period of rainfalls before the sampling in the areas. Precipitation data is available in the supplementary material (Figure S1). The correlations between PR and soil moisture were not significant (p > 0.05), regardless of the cultivation system (coppicing and stand renewal), harvesting method (HF and FS), and eucalypt straw level (100, 50, 0%) (Table 2), indicating that the differences in PR are due to the effects of treatments.

3.2. Soil Penetration Resistance in Coppicing Areas

Under coppicing areas, the harvester + forwarder (HF) treatment showed a PR average significantly lower (1.28 ± 0.24 MPa) than the feller + skidder treatment (FS) (2.11 ± 0.44 MPa) at 0–0.05 m depth layer (Figure 3). The lower PR value in HF was associated with the physical effect of eucalypt straw kept on the soil surface, which attenuates the soil compaction. For the underlying layers, regardless of the treatments, the confidence interval (CI) limits of PR were superposed, suggesting no significant differences between them. These results suggest that only the superficial layer (0–0.05 m) was affected by the eucalypt straw management as also found by other studies [12,35] reporting positive effects of forest straw kept on the soil surface, promoting reduced soil compaction. According to [35], plant straw act as a protective layer that reduces the contact pressure of machine wheels on the soil surface, dissipating the applied load and reducing the risks of soil compaction.
Figure 3 shows also that PR averages in both harvest systems were below 3 MPa, which was established as critical for eucalyptus roots development based on [16]. Less intense machinery traffic for harvesting in areas of coppicing possibly contributed to maintaining PR < 3 MPa. In forest soils, other studies have suggested that significant increments in PR only occur after several events of heavy traffic [36,37]. Besides that, the Oxisols are characterized by their stable structure and well drainage, which result in a higher deformation resistant [35,38]. These structure characteristics may have contributed to lower soil compaction during mechanized harvesting operation in the investigated Oxisol. Well drained soils have higher load support capacity, which reduces the risks of compaction [39].

3.3. Soil Penetration Resistance in Stand Renewal Areas

In areas under stand renewal and harvested by harvester + forwarder (HF) method with three eucalypt straw levels (100, 50, and 0%) kept on the soil surface, the averages of PR at 0–0.60 m were 1.09 ±0.18, 1.21 ± 0.17, and 1.25 ± 0.18 MPa for areas with straw levels of 100%, 50%, and 0%, respectively; in areas harvested with feller + skidder (FS) method, RP averages were 1.14 ±0.16, 1.16 ± 0.15, and 1.24 ± 0.18 MPa for areas with straw levels of 100%, 50%, and 0%, respectively. There were no significant differences in PR between soil layers (Figure 4) in both HF and FS, regardless of the amount of eucalypt straw on the soil surface, suggesting that maintenance or removal of eucalypt straw prior to soil tillage does not significantly affect PR values.
Despite the lack of differences between treatments (Figure 4), it should be considered that in areas where the eucalypt straw was completely removed, the subsequent renewals will be carried out on an increasingly smaller amount of straw on the soil surface, especially in areas harvested by feller + skidder method in which the trees are processed outside the plots. For this harvesting system, the practice of returning the eucalypt straw after harvest (represented here by the treatment FS with 50% of a straw level) leads to higher operational costs as well as a reduction of raw material as a source of energy for the industry. As already mentioned, the presence of harvest straw deposited on the soil surface can attenuate the effects of machinery traffic on PR, so that under stand renewal with subsequent cycles, much more attention should be paid to monitoring soil compaction where the eucalypt straw was removed before soil tillage.
The decreased values of PR, as well as its relative homogeneity between HF and FS treatments, are consequences of the current soil tillage, which was carried out shortly after the operations of harvest and straw management. It is well known that, in forest soils, soil tillage is usually carried out up to 0.50 m depth [9] and can reach up to 1 m depth in situations where there are impeding layers for root growth [40]. The soil physics conditions promoted by tillage can increase forest stand uniformity in renewal areas by decreases in root resistance to penetration [9] and improvements in water storage capacity [41]. Besides that, there is evidence that the effect of soil tillage on soil physical properties may persist for decades in forest soils [42].
The analysis of the PR spatial distribution in the harvesting systems HF and FS were done without considering the individual treatments of eucalypt straw level. The PR data showed moderate asymmetric (0.15 < |As| < 1.0) with right and platykurtic distribution curves (Ks < 0), which reflects frequencies close to the mean value (Table 3). According to the parameters obtained in the geostatistical analysis of PR, the layers 0–0.20 and 0.20–0.40 m showed a pure nugget effect (C = C + Co) for both harvesting systems (HF and FS) (Table 3). In contrast, in the 0.40–0.60 m soil depth layer, there was a weak (SD = 90.50%) and moderate (SD = 72.77%) spatial dependence of PR for HF and FS, respectively. This lack of spatial dependence of PR up to 0.40 m soil depth did not allow the fitting of any semivariance model, as well as the PR mapping by the kriging method. According to [43], the existence of spatial dependence is a prerequisite for applying the kriging method. For the 0.40–0.60 m layers, a linear model was fitted with a coefficient of determination (R2) of 0.14 and 0.69 for HF and FS, respectively, which indicated the weakness of the semivariance model to estimate the PR spatial variability in the studied area. Besides that, the large spatial dependence ranges (≥139.96 m) found for both linear models (HF and FS) do not have any practical application to justify the PR mapping in this soil depth layer.
Furthermore, these results indicated that the soil tillage reduced the spatial dependence of PR by increasing the soil’s physical and mechanical homogenization, especially up to 0.40 m soil depth. The occurrence of some spatial dependence only in the layer less influenced by the management (0.40–0.60 m) corroborates the observations made. According to [34], the weak spatial dependence is related to the external factors (e.g., soil tillage and agricultural machinery traffic), while the well-defined spatial dependence is more related to intrinsic factors associated with the natural soil structure. In other similar studies [24,26,44], the authors found spatial dependence of PR in areas under eucalyptus production. However, in these studies, the PR readings were carried out in eucalyptus stands before harvesting and soil tillage. Even under these conditions, [26] also observed a decrease in the PR spatial dependence from the bottom to the soil’s top layer. In the present study, the geospatial assessment allowed inferences about the spatial distribution of the PR measured after soil tillage.

4. Conclusions

We found in this investigation that the eucalypt straw kept on the soil surface contributed to mitigating the effects of machinery traffic on the penetration resistance of the 0–0.05 m layer depth during harvesting operations applied in areas under coppicing cultivation system. However, the harvesting operations did not increase the PR up to the critical (3 > MPa) limit to eucalypt roots development regardless of the harvest method (HF and FS) applied.
In stand renewal areas, the PR measured was not influenced by the different levels of eucalypt straw managed on the soil surface before soil tillage, regardless of the harvesting method used. The geospatial analysis of PR indicated that the soil tillage caused physical and mechanical homogenization of the soil, leading to uniform PR conditions for eucalyptus roots development. Additionally, future investigations may explore other physical indicators also related to plant development to monitor the short and long term changes in physical quality promoted by the forest straw management on the soil surface.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/agronomy12061482/s1, Figure S1: Precipitation and temperature between 1 June and 31 July 31 2017, measured by the nearest (38 km) whether station from the experimental areas. The weather station belongs to the National Institute of Meteorology (INMET), located in Ventania, Parana, Brazil, lat: −24.28027777; long: −50.21027777; level: 1093.41 m. Data are available on the INMET website: https://bdmep.inmet.gov.br/ (accessed 27 May 2022); Table S1: Particle-size distribution at the 0-0.20, 0.20-0.40, and 0.40-0.60 m layers in eucalypt cultivation areas under coppicing and stand renewal with different straw levels (100, 50 and 0%) both harvested by harvester + forwarder and feller + skidder systems.

Author Contributions

Conceptualization, J.d.M.O., T.S.d.O. and R.B.A.F.; methodology, J.d.M.O., T.S.d.O. and G.R.d.S.; software, J.d.M.O.; validation, J.d.M.O.; formal analysis, J.d.M.O. and C.A.T.; investigation, J.d.M.O. and C.A.T.; resources, J.d.M.O., T.S.d.O. and L.Z.; data curation, J.d.M.O.; writing—original draft preparation, J.d.M.O.; writing—review and editing, C.A.T. and L.Z.; visualization, J.d.M.O.; supervision, T.S.d.O. and C.A.T.; project administration, T.S.d.O.; funding acquisition, T.S.d.O. All authors have read and agreed to the published version of the manuscript.

Funding

This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brazil (CAPES)—Finance Code 001.

Acknowledgments

We would like to acknowledge the National Council of Technological and Scientific Development—Brazil (CNPq) and the Klabin S.A. company for the financial support. We also thank the Klabin S.A. company for the great interest and concession of their areas, allowing this work could be carried out, and to its employees James Stahl and his team for their assistance with field operation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. IBÁ. Relatorio Anual Iba 2021; IBRE Instituto Brasileiro de Economia, Ed.; Indústria Brasileira de Árvores: São Paulo, Brazil, 2021. [Google Scholar]
  2. Reichert, J.M.; Morales, C.A.S.; de Bastos, F.; Sampietro, J.A.; Cavalli, J.P.; de Araújo, E.F.; Srinivasan, R. Tillage Recommendation for Commercial Forest Production: Should Tillage Be Based on Soil Penetrability, Bulk Density or More Complex, Integrative Properties? Geoderma Reg. 2021, 25, e00381. [Google Scholar] [CrossRef]
  3. Cambi, M.; Certini, G.; Neri, F.; Marchi, E. Forest Ecology and Management The Impact of Heavy Traffic on Forest Soils. For. Ecol. Manag. 2015, 338, 124–138. [Google Scholar] [CrossRef]
  4. Nazari, M.; Eteghadipour, M.; Zarebanadkouki, M.; Ghorbani, M.; Dippold, M.A.; Bilyera, N.; Zamanian, K. Impacts of Logging-Associated Compaction on Forest Soils: A Meta-Analysis. Front. For. Glob. Chang. 2021, 4, 780074. [Google Scholar] [CrossRef]
  5. Hartmann, M.; Niklaus, P.A.; Zimmermann, S.; Schmutz, S.; Kremer, J.; Abarenkov, K.; Lüscher, P.; Widmer, F.; Frey, B. Resistance and Resilience of the Forest Soil Microbiome to Logging-Associated Compaction. ISME J. 2014, 8, 226–244. [Google Scholar] [CrossRef] [PubMed]
  6. Horn, R.; Vossbrink, J.; Becker, S. Modern Forestry Vehicles and Their Impacts on Soil Physical Properties. Soil Tillage Res. 2004, 79, 207–219. [Google Scholar] [CrossRef]
  7. Cambi, M.; Hoshika, Y.; Mariotti, B.; Paoletti, E.; Picchio, R.; Venanzi, R.; Marchi, E. Compaction by a Forest Machine Affects Soil Quality and Quercus Robur L. Seedling Performance in an Experimental Field. For. Ecol. Manag. 2017, 384, 406–414. [Google Scholar] [CrossRef]
  8. Solgi, A.; Najafi, A. The Impacts of Ground-Based Logging Equipment on Forest Soil. J. For. Sci. 2014, 60, 28–34. [Google Scholar] [CrossRef] [Green Version]
  9. Gonçalves, S.B.; Lopes, E.D.S.; Fiedler, N.C.; Cavalieri, K.M.V.; Stahl, J.; Drinko, C.H.F. Efeito Da Profundidade de Trabalho Na Qualidade Da Operação de Subsolagem Para Implantação Florestal. Rev. Arvoré 2016, 40, 29–37. [Google Scholar] [CrossRef] [Green Version]
  10. Curzon, M.T.; D’Amato, A.W.; Palik, B.J. Harvest Residue Removal and Soil Compaction Impact Forest Productivity and Recovery: Potential Implications for Bioenergy Harvests. For. Ecol. Manag. 2014, 329, 99–107. [Google Scholar] [CrossRef] [Green Version]
  11. Bortolas, E.P.; Stahl, J.; Minatel, R. Preparo de Solo Em Áreas de Reforma de Pinus e Eucalyptus Com Idade Superior a 20 Anos. SÉRIE TÉCNICA IPEF 2013, 17, 93–98. [Google Scholar]
  12. Lopes, E.D.S.; de Oliveira, D.; Rodrigues, C.K.; Drinko, C.H. Compactação de Um Solo Submetido Ao Tráfego Do Harvester e Do Forwarder Na Colheita de Madeira. Floresta Ambient. 2015, 22, 223–230. [Google Scholar] [CrossRef] [Green Version]
  13. Mason, W.L.; McKay, H.M.; Weatherall, A.; Connolly, T.; Harrison, A.J. The Effects of Whole-Tree Harvesting on Three Sites in Upland Britain on the Growth of Sitka Spruce over Ten Years. Forestry 2012, 85, 111–123. [Google Scholar] [CrossRef]
  14. Sands, R.; Greacen, E.L.; Gerard, C.J. Compaction of Sandy Soils in Radiata Pine Forest I, a Penetrometer Study. Aust. J. Soil Res. 1979, 17, 101–113. [Google Scholar] [CrossRef]
  15. Misra, R.K.; Gibbons, A.K. Growth and Morphology of Eucalypt Seedling-Roots, in Relation to Soil Strength Arising from Compaction. Plant Soil 1996, 182, 11. [Google Scholar] [CrossRef]
  16. Sinnett, D.; Morgan, G.; Williams, M.; Hutchings, T.R. Soil Penetration Resistance and Tree Root Development. Soil Use Manag. 2008, 24, 273–280. [Google Scholar] [CrossRef]
  17. Furtado, B.F.; Morris, L.A.; Markewitz, D. Loblolly Pine (Pinus taeda L.) Seedling Growth Response to Site Preparation Tillage on Upland Sites. Soil Sci. Soc. Am. J. 2016, 80, 472–489. [Google Scholar] [CrossRef] [Green Version]
  18. Mohieddinne, H.; Brasseur, B.; Spicher, F.; Gallet-Moron, E.; Buridant, J.; Kobaissi, A.; Horen, H. Physical Recovery of Forest Soil after Compaction by Heavy Machines, Revealed by Penetration Resistance over Multiple Decades. For. Ecol. Manag. 2019, 449, 117472. [Google Scholar] [CrossRef]
  19. Strock, C.F.; Rangarajan, H.; Black, C.K.; Schäfer, E.D.; Lynch, J.P. Theoretical Evidence That Root Penetration Ability Interacts with Soil Compaction Regimes to Affect Nitrate Capture. Ann. Bot. 2022, 129, 315–330. [Google Scholar] [CrossRef]
  20. Amirinejad, A.A.; Kamble, K.; Aggarwal, P.; Chakraborty, D.; Pradhan, S.; Mittal, R.B. Geoderma Assessment and Mapping of Spatial Variation of Soil Physical Health in a Farm. Geoderma 2011, 160, 292–303. [Google Scholar] [CrossRef]
  21. Vaz, C.M.P.; Manieri, J.M.; de Maria, I.C.; Tuller, M. Modeling and Correction of Soil Penetration Resistance for Varying Soil Water Content. Geoderma 2011, 166, 92–101. [Google Scholar] [CrossRef]
  22. Goutal, N.; Boivin, P.; Ranger, J. Assessment of the Natural Recovery Rate of Soil Specific Volume Following Forest Soil Compaction. Soil Sci. Soc. Am. J. 2012, 76, 1426–1435. [Google Scholar] [CrossRef]
  23. Barik, K.; Aksakal, E.L.; Islam, K.R.; Sari, S.; Angin, I. Spatial Variability in Soil Compaction Properties Associated with Field Traffic Operations. Catena 2014, 120, 122–133. [Google Scholar] [CrossRef]
  24. Filho, G.R.; de Passos e Carvalho, M.; Montanari, R.; da Silva, J.M.; Siqueira, G.M.; Zambianco, E.C. Variabilidade Espacial de Propriedades Dendrométricas Do Eucalipto e de Atributos Físicos de Um Latossolo Vermelho. Bragantia 2011, 70, 439–446. [Google Scholar] [CrossRef] [Green Version]
  25. Gonçalves, S.B.; Lopes, E.S.; Fiedler, N.C.; Cavalieri-polizeli, K.M.V.; Stahl, J. Resistência Do Solo a Penetração Em Diferentes Profundidades de Subsolagem. Nativa 2017, 5, 224–229. [Google Scholar] [CrossRef]
  26. Lima, C.G.D.R.; Panosso, A.R.; Costa, N.R.; De Carvalho, M.B.; Júnior, N.G.; da Silva, P.R.T. Sample Arrangements and Spatial Variability Characterization of Dendometrics Parameters of Eucalyptus Camaldulensis and Physical Soil Attributes. Rev. Ceres 2021, 68, 597–608. [Google Scholar] [CrossRef]
  27. Nitsche, P.R.; Caramori, P.H.; Ricce, W.D.S.; Pinto, L.F.D. Atlas Climático Do Estado Do Paraná. Available online: www.iapar.br (accessed on 5 May 2022).
  28. Soil Survey Staff Natural Resources Conservation Service, United States Department of Agriculture. Available online: http://websoilsurvey.nrcs.usda.gov/app/ (accessed on 31 August 2019).
  29. EMBRAPA. Manual de Métodos de Analise de Solo, 3rd ed.; Teixeira, P.C., Fontana, G.K.D.A., Teixeira, W.G., Eds.; Empresa Brasileira de Pesquisa Agropacuária, Centro Nacional de Pesquisa de Solos: Brasília, Brazil, 2017; ISBN 9788570357717. [Google Scholar]
  30. R Core Team. R: A Language and Environment for Statistical Computing; R Core Team: Vienna, Austria, 2019. [Google Scholar]
  31. Isaaks, H.E.; Srivastava, R.M. An Introduction to Applied Geostatistics; Oxford University Press: New York, NY, USA, 1989. [Google Scholar]
  32. Warrick, A.W. Spatial Variability. In Environmental Soil Physics; Hillel, D., Ed.; Academic Press: New York, NY, USA, 1998; pp. 655–675. [Google Scholar]
  33. Coelho, D.S.; Cortez, J.W.; Olszevski, N. Variabilidade Espacial Da Resisttênncia Mecânica à Penetração Em Vertissolo Cultivado Com Manga No Perímetro Irrigado de Mandacaru, Juazeiro, Bahia, Brasil. Rev. Bras. Cienc. Do Solo 2012, 36, 755–763. [Google Scholar] [CrossRef] [Green Version]
  34. Cambardella, C.A.; Moorman, T.B.; Novak, J.M.; Parkin, T.B.; Karlen, D.L.; Turco, R.F.; Konopka, A.E. Field-Scale Variability of Soil Properties in Central Iowa Soils. Soil Sci. Soc. Am. J. 1994, 58, 1501–1511. [Google Scholar] [CrossRef]
  35. Holthusen, D.; Brandt, A.A.; Reichert, J.M.; Horn, R.; Fleige, H.; Zink, A. Soil Functions and in Situ Stress Distribution in Subtropical Soils as Affected by Land Use, Vehicle Type, Tire Inflation Pressure and Plant Residue Removal. Soil Tillage Res. 2018, 184, 78–92. [Google Scholar] [CrossRef]
  36. Ampoorter, E.; Goris, R.; Cornelis, W.M.; Verheyen, K. Impact of Mechanized Logging on Compaction Status of Sandy Forest Soils. For. Ecol. Manag. 2007, 241, 162–174. [Google Scholar] [CrossRef]
  37. Williamson, J.R.; Neilsen, W.A. The Influence of Forest Site on Rate and Extent of Soil Compaction and Profile Disturbance of Skid Trails during Ground-Based Harvesting. Can. J. For. Res. 2000, 30, 1196–1205. [Google Scholar] [CrossRef]
  38. da Veiga, M.; Horn, R.; Reinert, D.J.; Reichert, J.M. Soil Compressibility and Penetrability of an Oxisol from Southern Brazil, as Affected by Long-Term Tillage Systems. Soil Tillage Res. 2007, 92, 104–113. [Google Scholar] [CrossRef]
  39. Richart, A.; Tavares Filho, J.; Brito, O.R.; Llanillo, R.F.; Ferreira, R. Compactação Do Solo: Causas e Efeitos. Semin. Ciências Agrárias 2005, 26, 321–344. [Google Scholar] [CrossRef]
  40. Fabiola, N.; Giarola, B.; da Silva, A.P.; Imhoff, S.; Dexter, A.R. Contribution of Natural Soil Compaction on Hardsetting Behavior. Geoderma 2003, 113, 95–108. [Google Scholar] [CrossRef]
  41. Stape, J.L.; Binkley, D.; Ryan, M.G.; Fonseca, S.; Loos, R.A.; Takahashi, E.N.; Silva, C.R.; Silva, S.R.; Hakamada, R.E.; Ferreira, J.M. de A.; et al. The Brazil Eucalyptus Potential Productivity Project: Influence of Water, Nutrients and Stand Uniformity on Wood Production. For. Ecol. Manag. 2010, 259, 1684–1694. [Google Scholar] [CrossRef]
  42. Dias Junior, M.D.S.; Leite, F.P.; Lasmar Junior, E.; Araújo Junior, F. Traffic Effects on the Soil Preconsolidation Pressure Due to Eucalyptus Harvest Operations. Sci. Agric. 2005, 62, 248–255. [Google Scholar] [CrossRef]
  43. Vieira, S.R. Geoestatística Em Estudos de Variabilidade Espacial Do Solo. In Tópicos em Ciência Do Solo; Novais, R.F., Alvarez, V.V.H., Schaefer, G.R., Eds.; Sociedade Brasileira de Ciência do Solo: Viçosa, Brazil, 2000; Volume 1, pp. 1–54. [Google Scholar]
  44. Lima, C.G.D.R.; Carvalho, M.D.P.E.C.; Narimatsu, K.C.P.; Da Silva, M.G.; Queiroz, H.A. De Atributos Físico-Químicos de Um Latossolo Do Cerrado Brasileiro e Sua Relação Com Caracteristicas Dendrométricas Do Eucalipto. Rev. Bras. Cienc. Do Solo 2010, 34, 163–173. [Google Scholar] [CrossRef]
Figure 1. (a) Location of the harvester + forwarder (HF) and feller + skidder (FS) treatments in coppicing areas and stand renewed areas with different eucalypt straw levels (100, 50 and 0%). Yellow dots represent the place of soil penetration resistance measurements. (b) Flowchart illustrating the scheme of treatments. The solid and dashed black arrows indicate that the eucalypt straw was removed or added before soil tillage in renewal areas, respectively.
Figure 1. (a) Location of the harvester + forwarder (HF) and feller + skidder (FS) treatments in coppicing areas and stand renewed areas with different eucalypt straw levels (100, 50 and 0%). Yellow dots represent the place of soil penetration resistance measurements. (b) Flowchart illustrating the scheme of treatments. The solid and dashed black arrows indicate that the eucalypt straw was removed or added before soil tillage in renewal areas, respectively.
Agronomy 12 01482 g001
Figure 2. Sampling scheme for penetration resistance (PR) and soil moisture in coppicing and stand renewed with different straw levels 100, 50 and 0% areas, both preceded by eucalyptus harvest with Harvester + Forwarder (HF) and Feller + Skidder (FS). Sampling points are defined by the vertices of isosceles triangles (15 m of height and 30 m of base).
Figure 2. Sampling scheme for penetration resistance (PR) and soil moisture in coppicing and stand renewed with different straw levels 100, 50 and 0% areas, both preceded by eucalyptus harvest with Harvester + Forwarder (HF) and Feller + Skidder (FS). Sampling points are defined by the vertices of isosceles triangles (15 m of height and 30 m of base).
Agronomy 12 01482 g002
Figure 3. Soil penetration resistance (PR) up to 0.60 m depth in eucalypt cultivation areas under coppicing harvested with harvester + forwarder (HF) and feller + skidder (FS) systems. Bars refer to the confidence interval (90%) of the average (n = 22).
Figure 3. Soil penetration resistance (PR) up to 0.60 m depth in eucalypt cultivation areas under coppicing harvested with harvester + forwarder (HF) and feller + skidder (FS) systems. Bars refer to the confidence interval (90%) of the average (n = 22).
Agronomy 12 01482 g003
Figure 4. Soil penetration resistance (PR) up to 0.60 m depth in eucalypt cultivation areas under stand renewal with different straw levels (100, 50, and 0%) and harvested with (A) harvester + forwarder (HF) and (B) feller + skidder (FS) systems. Bars refer to the confidence interval (90%) of the average (n = 22).
Figure 4. Soil penetration resistance (PR) up to 0.60 m depth in eucalypt cultivation areas under stand renewal with different straw levels (100, 50, and 0%) and harvested with (A) harvester + forwarder (HF) and (B) feller + skidder (FS) systems. Bars refer to the confidence interval (90%) of the average (n = 22).
Agronomy 12 01482 g004
Table 1. Soil chemical characterization and particle-size distribution at the 0–0.20, 0.20–0.40 and 0.40–0.60 m layers, in eucalypt cultivation areas under coppicing and stand renewal with different straw levels (100, 50 and 0%), both harvested by harvester + forwarder (HF) and feller + skidder (FS) systems.
Table 1. Soil chemical characterization and particle-size distribution at the 0–0.20, 0.20–0.40 and 0.40–0.60 m layers, in eucalypt cultivation areas under coppicing and stand renewal with different straw levels (100, 50 and 0%), both harvested by harvester + forwarder (HF) and feller + skidder (FS) systems.
Chemical Attributes ςParticle Size Distribution
LayerspH (H2O) PK+Ca2+Mg2+Al3+H + AlSOM §ClaySiltFine SandCoarse Sand
(m) mg dm−3cmolc dm−3dag kg−1kg kg−1
Coppicing areas
harvester + forwarder
0.00–0.204.42.223.00.310.392.0610.203.630.320.110.300.26
0.20–0.404.21.38.00.000.121.676.301.750.320.110.290.25
0.40–0.604.21.14.00.020.121.476.001.480.300.120.310.27
feller + skidder
0.00–0.204.22.027.00.190.142.659.504.030.330.090.290.27
0.20–0.404.21.38.00.030.072.067.902.550.270.140.310.27
0.40–0.604.11.36.00.010.071.968.902.280.230.190.330.28
Stand renewal areas δ
harvester + forwarder
0.00–0.205.02.147.02.490.540.598.903.760.280.100.360.26
0.20–0.404.81.823.02.440.450.497.303.490.300.090.340.26
0.40–0.604.51.415.00.520.161.186.502.020.310.110.320.26
feller + skidder
0.00–0.204.31.615.00.230.131.769.103.090.350.060.270.32
0.20–0.404.21.412.00.180.091.768.602.420.370.070.270.29
0.40–0.604.11.06.00.020.061.577.902.280.390.070.260.28
ς Methods described in [29]. § SOM: Soil organic matter. δ The soil chemical and particle-size distribution characterization in stand renewal areas was done without considering the individual treatments of eucalypt straw level. Additionally, the complete soil particle-size distribution data considering the individual treatments of eucalypt straw level is available in the supplemental material (Table S1).
Table 2. Soil moisture averages (θm) in the 0–0.20, 0.20–0.40 and 0.40–0.60 m layers, in eucalypt cultivation areas under coppicing and stand renewal with different straw levels (100, 50 and 0%), both harvested by harvester + forwarder (HF) and feller + skidder (FS) systems.
Table 2. Soil moisture averages (θm) in the 0–0.20, 0.20–0.40 and 0.40–0.60 m layers, in eucalypt cultivation areas under coppicing and stand renewal with different straw levels (100, 50 and 0%), both harvested by harvester + forwarder (HF) and feller + skidder (FS) systems.
LayersAverage ςCV §Normality Test ¥PR = f(θm) δ
(m)(kg kg−1)(%)Wp-ValueRp-Value
Coppicing areas
harvester + forwarder
0.00–0.200.21 ± 0.0322.880.950.62ns−0.0070.98ns
0.20–0.400.21 ± 0.0215.940.850.05ns
0.40–0.600.22 ± 0.0216.50.850.05ns
feller + skidder
0.00–0.200.20 ± 0.0215.240.960.84ns−0.070.85ns
0.20–0.400.20 ± 0.0112.840.950.65ns
0.40–0.600.21 ± 0.0112.110.950.69ns
Stand renewal areas
harvester + forwarder with straw level 100%
0.00–0.200.20 ± 0.0113.540.880.12ns−0.0930.80ns
0.20–0.400.20 ± 0.0112.670.970.90ns
0.40–0.600.20 ± 0.018.140.950.59ns
harvester + forwarder with straw level 50%
0.00–0.200.21 ± 0.019.70.960.72ns0.5970.54ns
0.20–0.400.22 ± 0.019.390.90.15ns
0.40–0.600.21 ± 0.0112.860.950.69ns
harvester + forwarder with straw level 0%
0.00–0.200.26 ± 0.0213.050.890.15ns0.230.52ns
0.20–0.400.26 ± 0.0111.190.920.32ns
0.40–0.600.26 ± 0.019.750.890.15ns
feller + skidder with straw level 100%
0.00–0.200.23 ± 0.018.710.940.48ns0.540.09ns
0.20–0.400.23 ± 0.0110.080.940.48ns
0.40–0.600.25 ± 0.017.810.930.37ns
feller + skidder with straw level 50%
0.00–0.200.25 ± 0.017.070.880.09ns−0.0610.85ns
0.20–0.400.26 ± 0.0110.490.940.53ns
0.40–0.600.27 ± 0.017.10.940.56ns
feller + skidder with straw level 0%
0.00–0.200.25 ± 0.016.630.940.53ns−0.070.36ns
0.20–0.400.26 ± 0.0110.530.850.05ns
0.40–0.600.27 ± 0.014.010.910.23ns
ς Average ± confidence interval (90% of confidence); § Coefficient of variation (CV); ¥ Normality test: Shapiro-Wilk test for normal distribution, where: (ns) not significant; δ Indexes of the correlation between soil moisture and penetration resistance (PR) in the 0–0.60 m soil depth layer: Coefficient of linear correlation (R) and Student’s t-test, where: (ns) not significant.
Table 3. Geostatistical parameters for soil penetration resistance (PR) at the 0–0.20, 0.20–0.40 and 0.40–0.60 m soil depth layers, in eucalypt cultivation areas under stand renewal and harvested by systems harvester + forwarder (HF) and feller + skidder (FS).
Table 3. Geostatistical parameters for soil penetration resistance (PR) at the 0–0.20, 0.20–0.40 and 0.40–0.60 m soil depth layers, in eucalypt cultivation areas under stand renewal and harvested by systems harvester + forwarder (HF) and feller + skidder (FS).
ParameterHarvester + Forwarder Feller + Skidder
Layer (m)0.00–0.200.20–0.400.40–0.600–0.200.20–0.400.40–0.60
Average (MPa)0.991.271.321.051.241.25
Median0.600.460.300.580.370.23
CV (%) §60.5135.9722.6654.9729.4518.58
Asymmetry0.50−0.010.300.29−0.200.25
Kurtosis−0.81−0.65−0.51−1.30−0.57−0.12
Sill0.000.000.010.000.000.02
Nugget0.340.200.090.350.140.04
Sill + Nugget 0.000.000.090.000.000.06
Range (m)0.000.00163.220.000.00139.96
SD (%) ς0.000.0090.500.000.0072.77
Modelnonolinearnonolinear
R2 ¥nono0.14nono0.69
§ CV: coefficient of variation; ς SD: index of spatial dependence; ¥ R2: coefficient of determination.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Oliveira, J.d.M.; Tormena, C.A.; dos Santos, G.R.; Zotarelli, L.; Fernandes, R.B.A.; de Oliveira, T.S. Soil Penetration Resistance Influenced by Eucalypt Straw Management under Mechanized Harvesting. Agronomy 2022, 12, 1482. https://doi.org/10.3390/agronomy12061482

AMA Style

Oliveira JdM, Tormena CA, dos Santos GR, Zotarelli L, Fernandes RBA, de Oliveira TS. Soil Penetration Resistance Influenced by Eucalypt Straw Management under Mechanized Harvesting. Agronomy. 2022; 12(6):1482. https://doi.org/10.3390/agronomy12061482

Chicago/Turabian Style

Oliveira, Judyson de M., Cássio A. Tormena, Gérson R. dos Santos, Lincoln Zotarelli, Raphael B. A. Fernandes, and Teógenes S. de Oliveira. 2022. "Soil Penetration Resistance Influenced by Eucalypt Straw Management under Mechanized Harvesting" Agronomy 12, no. 6: 1482. https://doi.org/10.3390/agronomy12061482

APA Style

Oliveira, J. d. M., Tormena, C. A., dos Santos, G. R., Zotarelli, L., Fernandes, R. B. A., & de Oliveira, T. S. (2022). Soil Penetration Resistance Influenced by Eucalypt Straw Management under Mechanized Harvesting. Agronomy, 12(6), 1482. https://doi.org/10.3390/agronomy12061482

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop