Acessibilidade / Reportar erro

Water erosion in surface soil conditions: runoff velocity, concentration and D50 index of sediments in runoff

ABSTRACT

Water erosion and contamination of water resources are influenced by concentration and diameter of sediments in runoff. This study aimed to quantify runoff velocity and concentration and the D50 index of sediments in runoff under different soil surface managements, in the following treatments: i) cropped systems: no-tilled soil covered by ryegrass (Lolium multiflorum Lam.) residue, with high soil cover and minimal roughness (HCR); no tilled soil covered by vetch (Vicia sativa L.) residue, with high soil cover and minimal roughness (HCV); chiseled soil after ryegrass crop removing the above-ground residues and keeping only the root system, with high roughness (HRR); chiseled soil after vetch crop removing the above-ground residues and keeping only the root system, with high roughness (HRV); ii) bare and chiseled soil, with high roughness (BHR). The research was conducted on a Humic Dystrupept under simulated rainfall. The design was completely randomized and each treatment was replicated twice. Eight rainfall events of controlled intensity (65 mm h−1) were applied to each treatment for 90 minutes. The D50 index, runoff velocity and sediment concentration were influenced by crop and soil management. Runoff velocity was more intensely reduced by cover crop residues than by surface roughness. Regardless of surface condition, the D50 index and concentration of sediment in runoff were lower under ryegrass than vetch crop. Runoff velocity and the D50 index were exponentially and inversely correlated with soil cover by residues and with surface roughness, while the D50 index was positively and exponentially correlated with runoff velocity.

hydric; erosion; soil management; conservation tillage

Introduction

Understanding the effect of water flow velocity, concentration and diameter of sediments in runoff is necessary to predict water erosion, mainly, the runoff capacity to transport chemicals that can contaminate water resources (Kuhn et al., 2012Kuhn, N.J.; Armstrong, E.K.; Ling, A.C.; Connolly, K.L.; Heckrath, G. 2012. Interrill erosion of carbon and phosphorus from conventionally and organically farmed Devon silt soils. Catena 91: 94–103.). Soil surface characteristics with greater influence on water erosion, related to runoff velocity and concentration and diameter of eroded sediments, are soil cover by crop residues and soil surface roughness (Cogo et al., 1984Cogo, N.P.; Moldenhauer, W.C.; Foster, G.R. 1984. Soil loss reductions from conservation tillage practices. Soil Science Society America Journal 48: 368-373.; Engel et al., 2009Engel, F.L.; Bertol, I.; Ritter, S.R.; Paz González, A.; Paz-Ferreiro, J.; Vidal Vázquez, E. 2009. Soil erosion under simulated rainfall in relation to phenological stages of soybeans and tillage methods in Lages, SC, Brazil. Soil & Tillage Research 103: 216-221.). These characteristics are usually found in conservative soil tillage (Bertol et al., 2003Bertol, I.; Mello, E.L.; Guadagnin, J.C.; Vedana Zaparolli, A.L.; Carrafa, M.R. 2003. Nutrient losses by water erosion. Scientia Agricola 60: 581-586.), represented by no-till and minimum tillage.

Crop residues can dissipate kinetic energy from raindrops and, in part, from runoff by slowing it down, hence, preventing soil disaggregation and sediment transport (Engel et al., 2009Engel, F.L.; Bertol, I.; Ritter, S.R.; Paz González, A.; Paz-Ferreiro, J.; Vidal Vázquez, E. 2009. Soil erosion under simulated rainfall in relation to phenological stages of soybeans and tillage methods in Lages, SC, Brazil. Soil & Tillage Research 103: 216-221.). Moreover, residues filter runoff, retaining coarse sediments and reducing sediment load. Roughness increases surface tortuosity and, therefore, decreases runoff velocity and its ability to detach the soil and transport sediments (Rodríguez-Caballero et al., 2012Rodríguez-Caballero, E.; Cantón, Y. Chamizo, S.; Afana, A.; Solé-Benet, A. 2012. Effects of biological soil crusts on surface roughness and implications for runoff and erosion. Geomorphology 145-146: 81–89.), however, this effect is ephemeral.

Grasses can be more efficient than leguminous crops to control water erosion. Grass biomass is usually composed of long and narrow leaves and stems and, for the same mass, it presents more efficient and persistent soil cover than broad-leaf plants do (Gilmour et al., 1998Gilmour, J.T.; Mauromoustakos, A.; Gale, P.M.; Norman, R.J. 1998. Kinetics of crop residue decomposition: variability among crops and years. Soil Science Society of American Journal 62: 750-755.). In addition, grass roots can improve soil aggregate stability and increase soil resistance to water erosion (Martinez-Trinidad et al., 2012Martinez-Trinidad, S.; Cotler, H.; Cruz-Cárdenas, G. 2012. The aggregates stability indicator <ASItest> to evaluate soil spatiotemporal change in a tropical dry ecosystem. Journal of Soil Science and Plant Nutrition 12: 363-377.).

Off-site erosion damage, especially related to surface water contamination, depends on the quantity and size of eroded sediments as well as their adsorption capacity of chemical species (Kuhn et al., 2012Kuhn, N.J.; Armstrong, E.K.; Ling, A.C.; Connolly, K.L.; Heckrath, G. 2012. Interrill erosion of carbon and phosphorus from conventionally and organically farmed Devon silt soils. Catena 91: 94–103.). Therefore, the higher the concentration of fine-size sediments in the flow, the higher the risk of environmental contamination (Barbosa et al., 2010Barbosa, F.T.; Bertol, I.; Luciano, R.V.; Paz-Ferreiro, J. 2010. Sediment size and organic carbon content in runoff and soil under two crops and two seed row directions. Revista Brasileira de Ciência do Solo 34: 1701-1710 (in Portuguese, with abstract in English).; Kuhn et al., 2012Kuhn, N.J.; Armstrong, E.K.; Ling, A.C.; Connolly, K.L.; Heckrath, G. 2012. Interrill erosion of carbon and phosphorus from conventionally and organically farmed Devon silt soils. Catena 91: 94–103.).

This study aimed to quantify and correlate runoff velocity, sediment concentration and the D50 index of sediments in runoff under different soil surface conditions. Therefore, a simulated rainfall study was performed in plots with soil covered by crop residues of ryegrass or vetch, soil chiseled after ryegrass or vetch crops where the above-ground residues were removed and, also, a bare and chiseled soil.

Materials and Methods

Description of the study site

The research was conducted in an experimental field with a mean slope of 0.134 m m−1 (moderately sloping surface), located in southern Brazil at coordinates 27º47’ S and 50o18’ W. The climate is Cfb type according to the Köppen classification with average annual rainfall of 1,535 mm (Schick et al., 2014Schick, J.; Bertol, I.; Cogo, N.P.; Paz González, A. 2014. Erosivity of rainfall in Lages, Santa Catarina, Brazil. Revista Brasileira de Ciência do Solo 38: 1890-1905 (in Portuguese, with abstract in English).). The altitude of the experimental site is 908 m and the soil is Humic Dystrupept.

Prior to this experiment, the study site was submitted to simulated rainfall with soybean (Glycine max L.) and maize (Zea mays L.) crops sown in contour and downhill in no tillage, except for the site where the bare soil treatment was installed, which was already uncropped. At the end of the last crop cycle, above-ground residues were removed and then ryegrass (Lolium multiflorum Lam.) and vetch (Vicia sativa L.) were broadcast sown and seeds were incorporated into the soil with light disking in the transverse direction of the slope. Ryegrass was sown at 40 kg ha−1, whereas vetch was sown at 60 kg ha−1.

The experimental plots were 11-m long in slope direction and 3.5-m wide. The lateral and upper plot boundaries were delimited by galvanized sheets with 0.2 m of height and were inserted 0.1 m into the soil. At the lower boundary, a runoff collector was installed, connected to a plastic tube with 75-mm diameter and 6-m long. The tube canalized runoff flow to a trench where measurements and runoff sampling were performed. The average plot slope was 0.134 m m−1, ranging from 0.124 m m−1 to 0.145 m m−1.

Experimental design

The experimental design was completely randomized. Five treatments were studied with two replications totaling ten experimental plots. The treatments with two replications were: i) no tilled soil with: a) ryegrass crop keeping residues on surface, with high soil cover (100 % of the soil covered at the beginning of the experiment) and minimum soil roughness (HCR); b) vetch crop keeping residues on surface, with high soil cover (90 % of the soil covered at the beginning of the experiment) and minimum soil roughness (HCV); ii) tilled soil with one chiseling transversely to the slope, resulting high roughness and low cover with: a) ryegrass crop, removing the above-ground residues and keeping the root system, resulting high roughness (20.5 mm) (HRR); b) vetch crop, removing the above-ground residues and keeping the root system resulting high roughness (17.4 mm) (HRV); iii) bare and no tilled soil, prepared for two years and chiseled transversely to the slope, with high roughness (14.6 mm) (BHR). The minimum roughness existing in the HCR and HCV treatments resulted from light harrowing performed for incorporating seeds into crop sowing. In the HRR and HRV treatments, the above-ground crop residues were removed after mowing the plants at full bloom stage. Afterward, the soil was chiseled on 13 December 2012, using a two-bar mounted device with 13 tines, spaced 0.25 m, operating at 0.15-m deep to produce high soil roughness.

Measurements

Eight simulated rainfall tests were applied to each treatment using a rotating-boom simulator, driven by hydraulic thrust, covering simultaneously two plots spaced 3.5 m apart. The rains were applied for 90 min, with constant planned intensity of 65 mm h−1, resulting in erosivity of 1,313 MJ mm ha−1 h−1. The first test was performed on 17 Dec 2011 and the others respectively on 01 Oct 2012, 07 Feb 2012, 10 Mar 2012, 11 May 2012, 18 Aug 2012, 02 Nov 2012 and 18 Dec 2012, for a year between the first and the eighth tests to reduce soil cover and surface roughness. Rainfall intensity was checked with 20 rain gauges previously installed on soil surface over the coverage area of the rainfall simulator that covered an area of approximately 300 m2. Between the rainfall tests, the plots were uncovered and subjected to the action of natural rain. Weeds were controlled by application of herbicides whenever necessary.

Surface roughness was determined immediately before soil tillage and prior to each simulated rainfall in the HRR, HRV and BHR treatments. The determinations were always carried out in the same place in the plot, marked by wooden stakes also used to level and support the micro relief meter. Surface roughness was obtained through a photographic record of a 20-bar set, where the bar heights in the photographs, representing the micro relief, were obtained by image interpretation, as described by Bertol et al. (2010)Bertol, I.; Zoldan Junior, W.A.; González, A.P.; Barbosa, F.T.; Werner, R.S. 2010. Sediment transport in runoff on rugous soil surface submitted to simulated rainfall Scientia Agricola 67: 591-597.. The random roughness index (RR) was calculated using the method proposed by Kamphorst et al. (2000)Kamphorst, E.C.; Jetten, V.; Guérif, J.; Pitkanen, J.; Iversen, N.V.; Douglas J.T.; Gonzáles, A.P. 2000. Predicting depressional storage from soil roughness. Soil Science Society of America Journal 64: 1749-1758. with no transformation of the data into logarithms or elimination of the extreme values.

Soil cover (SC) by crop residues was determined before each simulated rainfall in the HCR and HCV treatments using the marked rope method (Sloneker and Moldenhauer, 1977Sloneker, L.L.; Moldenhauer, W.C. 1977. Measuring the amounts of crop residue remaining after tillage. Journal of Soil Water and Conservation 32: 231-236.) with a 5-m long rope, containing marks at every 0.05 m, where soil coverage was obtained by counting for the total points with residues immediately below the marked points, using two replications per plot. The residue dry mass was also determined immediately before each simulated rainfall, in the HCR and HCV, collecting the residues from an area of 0.36 m2 in a single position randomly selected in the plot.

During the simulated rain, the runoff instantaneous rate was quantified by the methodology described by Cogo et al. (1984)Cogo, N.P.; Moldenhauer, W.C.; Foster, G.R. 1984. Soil loss reductions from conservation tillage practices. Soil Science Society America Journal 48: 368-373.. After starting, runoff was manually measured every five minutes with a graduated cylinder and a chronometer and a runoff sample was collected in a 0.75-L plastic pot to determine the concentration of sediments in runoff (CSR). To determine CSR, we made an average of five readings of runoff instantaneous rate with steady flood, because the erosion rate at this time was not influenced by previous discharge variation.

Surface runoff velocity (RV) was determined 70 min after the rainfall event began, when the measured runoff attained constant flow. The velocity was measured in the six central meters of each plot, marked with stakes (6 m are located 2.5 m above and 2.5 m below the extremities). At the upper extremity, methylene blue dye (2 %) was applied and time required for it to run this distance was recorded, as described by Engel et al. (2009)Engel, F.L.; Bertol, I.; Ritter, S.R.; Paz González, A.; Paz-Ferreiro, J.; Vidal Vázquez, E. 2009. Soil erosion under simulated rainfall in relation to phenological stages of soybeans and tillage methods in Lages, SC, Brazil. Soil & Tillage Research 103: 216-221..

Runoff samples to quantify the D50 index were collected after 80 min of each simulated rainfall event. The D50 index was quantified using a set of sieves with opening sizes of 4.75; 2; 1; and 0.25 mm. In that order, the sieves were placed on a bucket with 2 L capacity, placing the whole set under the flow until the bucket was completely filled. The content of the bucket was taken to the laboratory, where it was filtered through the sieves openings of 0.125; 0.053 and 0.038 mm, also considering the content that passed through the last sieve. All samples were transferred to an incubator at 50 °C for 72 h and then weighed. Thus, the following sizes of sediments were obtained: >4.75; 4.75-2; 2-1; 1-0.25; 0.25-0.125; 0.125-0.053; 0.053-0.038 and <0.038 mm. The D50 index was calculated using the procedure adopted by Gilley et al. (1987)Gilley, J.E.; Finkner, S.C.; Varvel, G.E. 1987. Size distribution of sediment as affected by surface residue and slope length. Transactions of the American Society of Agricultural Engineers 30: 1419-1424.. This index consists of a numerical value in which 50 % of the sediment mass in the runoff has a size greater than this value and the other 50 % of the sediment mass has a smaller size.

Statistical analysis

The results were analyzed using the statistical package SAS version 9.2 for Windows®, using a mixed model procedure. The Akaike information criterion (Wolfinger, 1993Wolfinger, R.D. 1993. Covariance structure selection in general mixed models. Communications in Statistics Simulation and Computation 22: 1079-1106.) was used to choose the variance and covariance matrices and detect the main effects such as causes of variance and interaction between them. The rain test effect (time) was considered fixed and the response variables (sediment concentration, runoff velocity and the D50 index) were the random effects. The means were compared by DMS of Fischer at 5 % significance level.

Results and Discussion

The average value of random roughness (RR) before the application of soil tillage was 4.1 mm, while after chiseling, RR increased to 17.5 mm (Table 1), agreeing with Bramorski et al. (2012)Bramorski, I.; Maria, I.C.; Silua. R.L.; Crestana, S. 2012. Relations between soil surface roughness, tortuosity, tillage treatments, rainfall intensity and soil and water losses from a Red Yellow Latosol. Revista Brasileira de Ciência do Solo 36: 1291-1297. who also worked with chiseling under natural rainfall. The lowest RR (9 mm) occurred in BHR, while in HRR and HRV treatments, the values were based on the average of the rainfall tests, 12.3 and 12.1 mm, respectively. This difference is explained by the absence of soil cultivation in BHR. The bare soil had lower organic carbon content and was less resistant to soil detachment, agreeing with Martinez-Trinidad et al. (2012)Martinez-Trinidad, S.; Cotler, H.; Cruz-Cárdenas, G. 2012. The aggregates stability indicator <ASItest> to evaluate soil spatiotemporal change in a tropical dry ecosystem. Journal of Soil Science and Plant Nutrition 12: 363-377..

Table 1
− Random soil surface roughness (RR) in treatments before tillage (BT) in each simulated rainfall test (T), soil cover with crop residues (SC) and dry mass of crop residues (DM) in a Humic Dystrupept.

The highest RR after chiseling was found in the HRR treatment, with 20.5 mm, followed by HRV and BHR treatments, with 17.4 mm and 14.6 mm, respectively (Table 1). The greatest RR observed in HRR after chiseling is due to a higher aggregation provided by ryegrass roots compared to vetch. According to Roisin (2007)Roisin, C.J.C. 2007. A multifractal approach for assessing the structural state of tilled soils. Soil Science Society of America Journal 71: 15-25., grass species increase water stability of soil aggregates compared to leguminous crops. The lower RR in BHR can be explained by the soil less resistance to breakdown (Martinez-Trinidad et al., 2012Martinez-Trinidad, S.; Cotler, H.; Cruz-Cárdenas, G. 2012. The aggregates stability indicator <ASItest> to evaluate soil spatiotemporal change in a tropical dry ecosystem. Journal of Soil Science and Plant Nutrition 12: 363-377.) and can be related to mechanical soil tillage and raindrop impact, increasing aggregate disruption in comparison with cropped treatments, as discussed by Engel et al. (2009)Engel, F.L.; Bertol, I.; Ritter, S.R.; Paz González, A.; Paz-Ferreiro, J.; Vidal Vázquez, E. 2009. Soil erosion under simulated rainfall in relation to phenological stages of soybeans and tillage methods in Lages, SC, Brazil. Soil & Tillage Research 103: 216-221..

RR decreased from the first to the eighth rainfall test in the HRR and HRV treatments, while in BHR, after reducing, RR increased slightly from the seventh test onward due to the formation of rill erosion on soil surface.

Residue biomass in the HCR and HCV treatments, determined immediately before the first test, was respectively of 4.7 and 4.4 Mg ha−1 (Table 1). After the first simulated rain test and before the second rain test, both treatments showed small reductions of residue biomass, probably due to the limited contact of this biomass with soil surface. After the second rain test, greater decay occurred in HCV (42 %) compared to HCR (13 %) due to the lower C/N ratio of vetch compared to ryegrass (Gilmour et al., 1998Gilmour, J.T.; Mauromoustakos, A.; Gale, P.M.; Norman, R.J. 1998. Kinetics of crop residue decomposition: variability among crops and years. Soil Science Society of American Journal 62: 750-755.). The greater decay observed in HCR occurred between the seventh and the eighth rain tests with a mass loss of 34 % (percentage of total decay).

On average, HCR covered 10 % more soil than HCV from the first to the fourth simulated rain tests (Table 1) ranging from 20 to 25 % between the fifth and the eighth. This difference is associated with a higher DM amount in HCR (Table 1), with its slower decomposition rate and with higher soil cover produced by ryegrass above-ground residues in relation to vetch for the same plant mass (Gilmour et al., 1998Gilmour, J.T.; Mauromoustakos, A.; Gale, P.M.; Norman, R.J. 1998. Kinetics of crop residue decomposition: variability among crops and years. Soil Science Society of American Journal 62: 750-755.).

The lowest concentration of sediments in runoff (CSR) (0.47 g L−1) was found in the HCR treatment (Table 2), followed by the HCV and HRR treatments (5.20 and 6.66 g L−1, respectively), while the highest values were found in the HRV and BHR treatments, 14.55 and 15.23 g L−1, respectively. These crop residues on the surface were effective to reduce CSR compared to bare soil, in agreement with the results found by Gilley et al. (2008)Gilley, J.E.; Sabatka W.F.; Eghball, B.; Marx D.B. 2008. Effect of Overland Flow Rate on Nutrient Transport from Land Application Areas. University of Nebraska, Lincoln, NE, USA., who reported that higher soil cover decreases soil losses. CSR was similar in the HCV and HRR treatments and protection against soil erosion caused by the ryegrass root system in HRR was equivalent to vetch residue cover in HCV. Furthermore, soil surface microrelief produced by chiseling transversely to the slope concentrates erosion in the channels, increasing sediment transport in these treatments, also reported by Zheng et al. (2004).

Table 2
− Sediment concentration in runoff (CSR), runoff velocity (RV) and the D50 index (D50) in treatments and tests of simulated rain (T) in a Humic Dystrupept.

Comparing the treatments with crop residues on soil surface, CSR was 91 % lower in HCR than in HCV (Table 2), a similar difference was reported by Zheng et al. (2004) using wheat (Triticum aestivum L.) residue that was more effective to reduce soil erosion compared maize residue. Among the chiseled treatments, HRR differed from HRV and BHR in terms of CSR. This is related to the greater soil aggregation in soil cropped with grass rather than legumes, also reported by Roisin (2007)Roisin, C.J.C. 2007. A multifractal approach for assessing the structural state of tilled soils. Soil Science Society of America Journal 71: 15-25. and Martinez-Trinidad et al. (2012)Martinez-Trinidad, S.; Cotler, H.; Cruz-Cárdenas, G. 2012. The aggregates stability indicator <ASItest> to evaluate soil spatiotemporal change in a tropical dry ecosystem. Journal of Soil Science and Plant Nutrition 12: 363-377.. Better soil structure increases soil resistance to detachment by raindrop impact and surface flow. The HRV and BHR treatments presented similar CSR due to low influence of vetch on soil aggregation, similar to that of bare soil.

On average, the lowest runoff velocity (RV) was found (Table 2) in the HCR treatment (0.063 m s−1) and the highest in BHR (0.190 m s−1). The low RV in HCR is explained by the formation of a physical barrier by crop residues, increasing flow tortuosity, also observed by Vidal Vázquez et al. (2008)Vidal Váquez, E.; Garcia Moreno, R.; Miranda, J.G.V; Díaz, M.C.; Saá Requejo, A.; Paz-Ferreiro, J.; Tarquis, A.M. 2008. Assessing soil surface roughness decay during simulated rainfall by multifractal analysis. Nonlinear Process Geophysical 15: 457-468.. Thus, in HCR, residue decomposition was slower and biomass cover for the same amount of residue was greater than in HCV, increasing soil protection against erosion.

The effect of soil surface conditions on RV was described by an exponential model (Figure 1). This RV reduction is important since it has great influence on the final flow energy, thus, on its capacity to detach and transport coarse soil particles and aggregates (Lambrechts et al., 2013Lambrechts, T.; François, S.; Lutts, S.; Muñoz-Carpena, R.; Bielders, C.L. 2013. Impact of plant growth and morphology and of sediment concentration on sediment retention efficiency of vegetative filter strips: flume experiments and VFSMOD modeling. Journal of Hydrology 511: 800–810.). In bare soil (BHR), the exponential model was weakly adjusted to the data and this treatment had the lowest efficiency to reduce RV due to its low soil aggregate resistance, resulting in fast flattening and sealing of soil surface, which accelerate and increase runoff flow energy, as reported for this soil condition by Bertol et al. (2010)Bertol, I.; Zoldan Junior, W.A.; González, A.P.; Barbosa, F.T.; Werner, R.S. 2010. Sediment transport in runoff on rugous soil surface submitted to simulated rainfall Scientia Agricola 67: 591-597..

Figure 1
− Relation (a) between runoff velocity (RV) and soil cover (SC) in no-tilled, cropped and covered soil treatments, with ryegrass (HCR) and vetch (HCV) residues and (b) between RV and random soil surface roughness (RR) in chiseled soil, cropped with ryegrass, removing above-ground residues and keeping the root system (HRR); chiseled soil, cropped with vetch, removing above-ground residues and keeping the root system (HRV); and bare and chiseled soil (BHR). These adjustments were based on the average of the rainfall tests.

Crop residues on soil surface in the HCR and HCV treatments were more effective to reduce RV than roughness caused by the action of chisel plow in HRR, HRV and BHR (Table 2, Figure 1). The decrease was on average 48 % in the treatments with residues compared to the treatments with chiseling. This is attributed to the more effective physical barrier to control runoff provided by residues compared to the action of roughness, increasing flow tortuosity, also reported by Vidal Vázquez et al. (2008)Vidal Váquez, E.; Garcia Moreno, R.; Miranda, J.G.V; Díaz, M.C.; Saá Requejo, A.; Paz-Ferreiro, J.; Tarquis, A.M. 2008. Assessing soil surface roughness decay during simulated rainfall by multifractal analysis. Nonlinear Process Geophysical 15: 457-468.. Among the residues, ryegrass in the HCR treatment was more effective than vetch in HCV, decreasing RV by 27 % (Figure 1A). This is supported by Lopes et al. (1987)Lopes, P.R.C.; Cassol, E.A.; Cogo, N.P. 1987. Influence of crop residue mulch on water flow velocity and size distribution of eroded soil sediments. Revista Brasileira de Ciência do Solo 11: 193-197 (in Portuguese with abstract in English). and Gilmour et al. (1998)Gilmour, J.T.; Mauromoustakos, A.; Gale, P.M.; Norman, R.J. 1998. Kinetics of crop residue decomposition: variability among crops and years. Soil Science Society of American Journal 62: 750-755. that reported changes in soil cover due to different crop species, considering the same amount of residue, which was confirmed by the percentage of soil cover in those treatments (Table 1) with similar residue mass.

Comparing the treatments in which the soil was chiseled, the lowest RV values were observed in the cropped soil (Table 2). In the HRR and HRV treatments, RV values were, on average, 37 % lower than in BHR. This result is justified by greater RR kept over time, in the cropped treatments, compared to BHR (Figure 1B). Greater roughness generated higher tortuosity over time in cropped and chiseled soils, slowing runoff flow, as reported by Engel et al. (2009)Engel, F.L.; Bertol, I.; Ritter, S.R.; Paz González, A.; Paz-Ferreiro, J.; Vidal Vázquez, E. 2009. Soil erosion under simulated rainfall in relation to phenological stages of soybeans and tillage methods in Lages, SC, Brazil. Soil & Tillage Research 103: 216-221..

The lowest D50 index of sediments (0.18 mm) occurred in the HCR treatment, considering the average of the rainfall tests (Table 2), confirming the protective effect of ryegrass residue to dissipate the kinetic energy of raindrops and reduce soil detachment. The fact that the soil was not tilled in this treatment also supports this result. Crop residue cover reduces runoff and increases sediment deposition and it also filters coarse sediments, as reported by Gilmour et al. (1998)Gilmour, J.T.; Mauromoustakos, A.; Gale, P.M.; Norman, R.J. 1998. Kinetics of crop residue decomposition: variability among crops and years. Soil Science Society of American Journal 62: 750-755.. However, fine-size sediments transported in the runoff originating from conservational soil tillage are generally richer in nutrients and organic C than coarser sediments are, as reported by Foster et al. (1985)Foster, G.R.; Young, R.A.; Neibling, W.H. 1985. Sediment composition for nonpoint source pollution analyses. Transactions of the American Society of Agricultural Engineers 28: 133-139., as coarse sediments have less absorption capacity of chemicals, according to Kuhn et al. (2012)Kuhn, N.J.; Armstrong, E.K.; Ling, A.C.; Connolly, K.L.; Heckrath, G. 2012. Interrill erosion of carbon and phosphorus from conventionally and organically farmed Devon silt soils. Catena 91: 94–103.. Bertol et al. (2003)Bertol, I.; Mello, E.L.; Guadagnin, J.C.; Vedana Zaparolli, A.L.; Carrafa, M.R. 2003. Nutrient losses by water erosion. Scientia Agricola 60: 581-586. found greater P concentration and losses in no-tillage in desiccated field. Zhang et al. (2011)Zhang, G.H.; Liu, G.B.; Wang, G.L.; Wang, Y.X. 2011. Effects of vegetation cover and rainfall intensity on sediment-associated nitrogen and phosphorus losses and particle size composition on the Loess Plateau. Journal of Soil and Water Conservation 66: 192-200. found a positive relation between N and P adsorption with clay content in eroded sediments.

Considering the average of the rainfall tests, the greater value of the D50 index of sediments (1.56 mm) occurred in the BHR treatment (Table 2). This value is nine-fold greater than that of HCR due to the action of soil cover by ryegrass residue and can be related to the lower aggregation and protection against erosion in the case of bare soil. These conditions increase soil detachment and sediment transport in BHR, as reported by De Baets and Poesen (2010)De Baets, S.; Poesen, J. 2010. Empirical models for predicting the erosion-reducing effect of plant roots during concentrated flow. Geomorphology 118: 425–432. and Martinez-Trinidad et al. (2012)Martinez-Trinidad, S.; Cotler, H.; Cruz-Cárdenas, G. 2012. The aggregates stability indicator <ASItest> to evaluate soil spatiotemporal change in a tropical dry ecosystem. Journal of Soil Science and Plant Nutrition 12: 363-377., who found influence of soil cultivation and aggregate stability on soil resistance to erosion. Thus, bare soil allows surface sealing and concentrates runoff in channels, resulting in greater flow velocity in BHR compared to the covered soil (Table 2).

The D50 index found in the HCV treatment was similar to that for HRR (Table 2), disagreeing with results reported by Lopes et al. (1987)Lopes, P.R.C.; Cassol, E.A.; Cogo, N.P. 1987. Influence of crop residue mulch on water flow velocity and size distribution of eroded soil sediments. Revista Brasileira de Ciência do Solo 11: 193-197 (in Portuguese with abstract in English)., who found that residues filter coarse sediments in the flow in relation to the tilled soil. This result can be related to the greater soil aggregation caused by ryegrass roots in HRR (Martinez-Trinidad et al., 2012Martinez-Trinidad, S.; Cotler, H.; Cruz-Cárdenas, G. 2012. The aggregates stability indicator <ASItest> to evaluate soil spatiotemporal change in a tropical dry ecosystem. Journal of Soil Science and Plant Nutrition 12: 363-377.) and by the lower efficiency of vetch residue to filter coarse sediments in HCV, as reported by Barbosa et al. (2010)Barbosa, F.T.; Bertol, I.; Luciano, R.V.; Paz-Ferreiro, J. 2010. Sediment size and organic carbon content in runoff and soil under two crops and two seed row directions. Revista Brasileira de Ciência do Solo 34: 1701-1710 (in Portuguese, with abstract in English)..

Regardless of the management practice, ryegrass generated a smaller D50 index of sediments than vetch did (Table 2). Ryegrass decreased D50 by 51 % in relation to vetch, when covered treatments (HCR and HCV) and chiseled (HRR and HRV) were grouped in terms of crop type. When considered separately, the HCR and HRR treatments decreased D50 by 78 % and 31 % in relation to HCV and HRV, respectively. Lopes et al. (1987)Lopes, P.R.C.; Cassol, E.A.; Cogo, N.P. 1987. Influence of crop residue mulch on water flow velocity and size distribution of eroded soil sediments. Revista Brasileira de Ciência do Solo 11: 193-197 (in Portuguese with abstract in English). found no difference in the D50 index comparing treatments with wheat, maize and soybean under different management practices. However, in a study with oats (Avena strigosa S.) and vetch cover crops, Barbosa et al. (2010)Barbosa, F.T.; Bertol, I.; Luciano, R.V.; Paz-Ferreiro, J. 2010. Sediment size and organic carbon content in runoff and soil under two crops and two seed row directions. Revista Brasileira de Ciência do Solo 34: 1701-1710 (in Portuguese, with abstract in English). found lower D50 values in oats compared to vetch crop, regardless of the sowing direction. These authors attributed the small D50 under oats crop to its positive effect on soil aggregation due to its high density of fine roots, agreeing with De Baets and Poesen (2010)De Baets, S.; Poesen, J. 2010. Empirical models for predicting the erosion-reducing effect of plant roots during concentrated flow. Geomorphology 118: 425–432., who reported that increased density of fine roots increased soil aggregation and its resistance to detachment.

The D50 index showed an inverse exponential relationship with soil cover and surface roughness (Figure 2), agreeing with Bertol et al. (2010)Bertol, I.; Zoldan Junior, W.A.; González, A.P.; Barbosa, F.T.; Werner, R.S. 2010. Sediment transport in runoff on rugous soil surface submitted to simulated rainfall Scientia Agricola 67: 591-597. and Panuska et al. (2008)Panuska, J.C.; Karthikeyan, K.G.; Miller, P.S. 2008. Impact of surface roughness and crusting on particle size distribution of edge-of-field sediments. Geoderma 145: 315-324.. For the treatments with surface residues, HCR and HCV, the coefficients of determination between the D50 index with soil coverage were 0.73 and 0.47, respectively (Figure 2A). For the chiseled treatments, HRR, HRV and BHR, the coefficients of determination between the D50 index with soil roughness were 0.72, 0.77 and 0.32, respectively (Figure 2B). The low coefficient of determination in BHR is related to the increase in soil surface roughness, observed in the seventh rainfall test, as previously discussed.

Figure 2
− Relation (a) between the D50 index and soil cover (SC) in no-tilled, cropped and covered soil treatments, with ryegrass (HCR) and vetch (HCV) residues and (b) between the D50 and random soil surface roughness (RR) in chiseled soil, cropped with ryegrass, removing above-ground residues and keeping the root system (HRR); chiseled soil, cropped with vetch, removing above-ground residues and keeping the root system (HRV); and bare and chiseled soil (BHR). These adjustments were based on the average of the rainfall tests.

Regardless of the treatment, RV increase was exponentially correlated with the increase in the D50 index (Figure 3), showing that RV increased sediment transport capacity, also verified by Bertol et al. (2010)Bertol, I.; Zoldan Junior, W.A.; González, A.P.; Barbosa, F.T.; Werner, R.S. 2010. Sediment transport in runoff on rugous soil surface submitted to simulated rainfall Scientia Agricola 67: 591-597.. The adjustment degree of these variables was higher in the HRR and HCR treatments in which the D50 index was explained by RV 79 % and 88 % respectively, compared to HCV and HRV, where the relations were 50 % and 25 %, respectively. In BHR, the adjustment degree of these variables was also high, indicating that 74 % of the D50 index depended on RV.

Figure 3
− Relation between the D50 index and runoff velocity (RV) in no-tilled, cropped and covered soil treatments, with ryegrass (HCR) and vetch (HCV) residues; in chiseled soil, cropped with ryegrass, removing above-ground residues and keeping the root system (HRR); chiseled soil, cropped with vetch, removing above-ground residues and keeping the root system (HRV); and bare and chiseled soil (BHR). These adjustments were based on the average of the rainfall tests.

Conclusions

Soil cover with crop residues is more effective than surface roughness to reduce runoff velocity, while chiseled soil after removing the above-ground residues and keeping only the root system of crops is more effective than the bare soil chiseled to reduce runoff velocity.

No-tilled and cropped soil with vetch has concentration and D50 index of sediments transported by runoff similar to those of chiseled soil cropped with ryegrass, removing above-ground residue and keeping the root system. Chiseled and cropped soil with vetch removing above-ground residue and keeping the root system has a sediment concentration similar to that of bare and chiseled soil.

Runoff velocity and the D50 index of transported sediments are exponentially and inversely correlated with soil cover with crop residues and with soil surface roughness, while the D50 index of sediments is exponentially and positively correlated with runoff velocity.

Acknowledgements

We would like to thank the Foundation for Research and Innovation of the State of Santa Catarina (FAPESC) and the Brazilian National Council for Scientific and Technological Development (CNPq) for the financial support.

References

  • Barbosa, F.T.; Bertol, I.; Luciano, R.V.; Paz-Ferreiro, J. 2010. Sediment size and organic carbon content in runoff and soil under two crops and two seed row directions. Revista Brasileira de Ciência do Solo 34: 1701-1710 (in Portuguese, with abstract in English).
  • Bertol, I.; Mello, E.L.; Guadagnin, J.C.; Vedana Zaparolli, A.L.; Carrafa, M.R. 2003. Nutrient losses by water erosion. Scientia Agricola 60: 581-586.
  • Bertol, I.; Zoldan Junior, W.A.; González, A.P.; Barbosa, F.T.; Werner, R.S. 2010. Sediment transport in runoff on rugous soil surface submitted to simulated rainfall Scientia Agricola 67: 591-597.
  • Bramorski, I.; Maria, I.C.; Silua. R.L.; Crestana, S. 2012. Relations between soil surface roughness, tortuosity, tillage treatments, rainfall intensity and soil and water losses from a Red Yellow Latosol. Revista Brasileira de Ciência do Solo 36: 1291-1297.
  • Cogo, N.P.; Moldenhauer, W.C.; Foster, G.R. 1984. Soil loss reductions from conservation tillage practices. Soil Science Society America Journal 48: 368-373.
  • De Baets, S.; Poesen, J. 2010. Empirical models for predicting the erosion-reducing effect of plant roots during concentrated flow. Geomorphology 118: 425–432.
  • Engel, F.L.; Bertol, I.; Ritter, S.R.; Paz González, A.; Paz-Ferreiro, J.; Vidal Vázquez, E. 2009. Soil erosion under simulated rainfall in relation to phenological stages of soybeans and tillage methods in Lages, SC, Brazil. Soil & Tillage Research 103: 216-221.
  • Foster, G.R.; Young, R.A.; Neibling, W.H. 1985. Sediment composition for nonpoint source pollution analyses. Transactions of the American Society of Agricultural Engineers 28: 133-139.
  • Gilley, J.E.; Finkner, S.C.; Varvel, G.E. 1987. Size distribution of sediment as affected by surface residue and slope length. Transactions of the American Society of Agricultural Engineers 30: 1419-1424.
  • Gilley, J.E.; Sabatka W.F.; Eghball, B.; Marx D.B. 2008. Effect of Overland Flow Rate on Nutrient Transport from Land Application Areas. University of Nebraska, Lincoln, NE, USA.
  • Gilmour, J.T.; Mauromoustakos, A.; Gale, P.M.; Norman, R.J. 1998. Kinetics of crop residue decomposition: variability among crops and years. Soil Science Society of American Journal 62: 750-755.
  • Kamphorst, E.C.; Jetten, V.; Guérif, J.; Pitkanen, J.; Iversen, N.V.; Douglas J.T.; Gonzáles, A.P. 2000. Predicting depressional storage from soil roughness. Soil Science Society of America Journal 64: 1749-1758.
  • Kuhn, N.J.; Armstrong, E.K.; Ling, A.C.; Connolly, K.L.; Heckrath, G. 2012. Interrill erosion of carbon and phosphorus from conventionally and organically farmed Devon silt soils. Catena 91: 94–103.
  • Lambrechts, T.; François, S.; Lutts, S.; Muñoz-Carpena, R.; Bielders, C.L. 2013. Impact of plant growth and morphology and of sediment concentration on sediment retention efficiency of vegetative filter strips: flume experiments and VFSMOD modeling. Journal of Hydrology 511: 800–810.
  • Lopes, P.R.C.; Cassol, E.A.; Cogo, N.P. 1987. Influence of crop residue mulch on water flow velocity and size distribution of eroded soil sediments. Revista Brasileira de Ciência do Solo 11: 193-197 (in Portuguese with abstract in English).
  • Martinez-Trinidad, S.; Cotler, H.; Cruz-Cárdenas, G. 2012. The aggregates stability indicator <ASItest> to evaluate soil spatiotemporal change in a tropical dry ecosystem. Journal of Soil Science and Plant Nutrition 12: 363-377.
  • Panuska, J.C.; Karthikeyan, K.G.; Miller, P.S. 2008. Impact of surface roughness and crusting on particle size distribution of edge-of-field sediments. Geoderma 145: 315-324.
  • Rodríguez-Caballero, E.; Cantón, Y. Chamizo, S.; Afana, A.; Solé-Benet, A. 2012. Effects of biological soil crusts on surface roughness and implications for runoff and erosion. Geomorphology 145-146: 81–89.
  • Roisin, C.J.C. 2007. A multifractal approach for assessing the structural state of tilled soils. Soil Science Society of America Journal 71: 15-25.
  • Schick, J.; Bertol, I.; Cogo, N.P.; Paz González, A. 2014. Erosivity of rainfall in Lages, Santa Catarina, Brazil. Revista Brasileira de Ciência do Solo 38: 1890-1905 (in Portuguese, with abstract in English).
  • Sloneker, L.L.; Moldenhauer, W.C. 1977. Measuring the amounts of crop residue remaining after tillage. Journal of Soil Water and Conservation 32: 231-236.
  • Vidal Váquez, E.; Garcia Moreno, R.; Miranda, J.G.V; Díaz, M.C.; Saá Requejo, A.; Paz-Ferreiro, J.; Tarquis, A.M. 2008. Assessing soil surface roughness decay during simulated rainfall by multifractal analysis. Nonlinear Process Geophysical 15: 457-468.
  • Wolfinger, R.D. 1993. Covariance structure selection in general mixed models. Communications in Statistics Simulation and Computation 22: 1079-1106.
  • Zhang, G.H.; Liu, G.B.; Wang, G.L.; Wang, Y.X. 2011. Effects of vegetation cover and rainfall intensity on sediment-associated nitrogen and phosphorus losses and particle size composition on the Loess Plateau. Journal of Soil and Water Conservation 66: 192-200.

Edited by

Edited by: Silvia del Carmen Imhoff

Publication Dates

  • Publication in this collection
    May-Jun 2016

History

  • Received
    11 Mar 2015
  • Accepted
    17 Sept 2015
Escola Superior de Agricultura "Luiz de Queiroz" USP/ESALQ - Scientia Agricola, Av. Pádua Dias, 11, 13418-900 Piracicaba SP Brazil, Phone: +55 19 3429-4401 / 3429-4486 - Piracicaba - SP - Brazil
E-mail: scientia@usp.br