|
|
||||||||
Department of Soil Management and Soil Care, Ghent Univ., Coupure Links 653, B-9000 Ghent, Belgium
* Corresponding authors (donald.gabriels{at}UGent.be).
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: Ek, kinetic energy of the rain on the soil surface ER, enrichment ratio OC, organic carbon qs, unit sediment discharge
| INTRODUCTION |
|---|
|
|
|---|
Total soil loss by erosion is the result of different processes such as deposition, splash, and interrill, rill, and gully erosion. These processes can cause a preferential transport of OC. Ghadiri and Rose (1991a) found that the outer layers of soil aggregates contained more N, organic matter, and pesticides than the inner core. Therefore, "stripping" of the aggregates by raindrop impact caused an enrichment of these soil components in the eroded material. Most research has shown that the enrichment of OC in suspended sediment is related to a selective transport of clay particles, indicating that the enrichment mechanism is mainly controlled by the transport capacity of overland flow (e.g., Sharpley, 1985; Weigand et al., 1998). Enrichment of OC can also be caused by deposition: large particles and aggregates settle out more quickly, resulting in a physical enrichment of smaller particles and less dense aggregates in runoff (Davis et al., 1983; Lu et al., 1988; Beuselinck, 1999).
The selectivity of the erosion processes is often expressed as an enrichment ratio, ER (Massey and Jackson, 1952), calculated as the ratio of the concentration of a soil constituent in the eroded material to its concentration in the original soil. The enrichment ratio of OC, EROC, has been related to OC content in the topsoil (Young and Onstad, 1976; Avnimelech and McHenry, 1984) and to soil loss (Sharpley, 1985; Ghadiri and Rose, 1991b). In general, EROC values in suspended sediment have been found to be >1. For a sandy loam soil, Cogle et al. (2002) reported EROC values between 1.3 and 4.1, with larger variations in EROC at low sediment yields. Jacinthe et al. (2002) conducted laboratory rainfall simulations on small soil pans with undisturbed samples from fields under different tillage techniques. For no-till, moldboard plow, and chisel till, they reported EROC values of 1.4, 1.7, and 2.3 at sediment delivery rates of 17.4, 34.1, and 48.9 g m–2 h–1, respectively. This indicates an increase in the EROC with increasing soil loss rate, which is rather unusual and may be a result of the difference in tillage technique; however, the results of Cogle et al. (2002) indicated that the kind of tillage technique affected the sediment yields but did not directly influence the EROC. This was also observed by Basic et al. (2002) and Kisic et al. (2002) for enrichment of clay and OC, respectively. On the other hand, Fullen (1991) and Fullen et al. (1996) found lower OC, clay, and silt concentrations in eroded material compared with the loamy sand topsoil and attributed it to eluviation of clay and organic matter into the topsoil during runoff events. This in-washing process resulted in a "filtration pavement" immediately beneath the surface crust.
The enrichment in OC is also affected by the storm intensity. The results of Jacinthe et al. (2004) suggest a more selective detachment and sorting of labile OC fractions during low-intensity storms than high-intensity storms. They stated that the high soil losses associated with intense rainstorms can determine the overall impact of erosional events on losses of labile OC. On the other hand, field plot experiments conducted by Quinton et al. (2001) showed that the enrichment of sediment during small events should not be neglected in estimating the loss of soil constituents by erosion.
The cited research results indicate that, in general, suspended sediment is enriched in OC. To what extent different erosion processes account for this enrichment mechanism, however, is still unclear. This knowledge is required for modeling OC losses. Polyakov and Lal (2004a) concluded in their review that only a few soil organic matter models incorporate an erosion component in their routine. These erosion components are based on empirical models like the Universal Soil Loss Equation, and are mostly unable to represent erosion-induced soil organic matter loss and deposition as spatially and temporally distributed phenomena. Erosion and deposition processes are dependent on spatially variable factors. Therefore, assessment of total soil loss in a lumped way does not take into account the spatial variation in enrichment processes. Moreover, many models use constant values for ER or (usually logarithmic) equations that relate ER to average soil loss during a rainfall event or a whole year. In this way, the effect of temporal variations in erosion intensity on the ER is not taken into account.
The purpose of this study was to examine the enrichment of OC caused by rill and interrill erosion processes simulated in the laboratory. Flume experiments and rainfall simulations were conducted on small soil pans to find a relationship between OC losses and erosion. Based on these experiments, also the influence of rainfall impact and interrill and rill erosion on EROC was investigated. In this study, we only considered suspended material in runoff. Therefore, the term sediment refers to suspended sediment and not deposited material.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Interrill Erosion Experiments
The soil was air dried and ground (<8 mm) before use in the experiments. The sieved soil was spread out equally in soil pans (0.55 m long, 0.20 m wide, and 0.035 m deep) and slightly compacted with a wooden plate to obtain a bulk density of 1.33 Mg m–3, equaling the conditions of the field plots at Maarkedal. Each soil pan had a perforated bottom, allowing free drainage of percolating water. The soil pans were set at slopes of 10 and 30%.
To examine the influence of rainfall kinetic energy on the ER, a plastic cover was applied: small plastic sheets (20 by 20 mm) were placed on wooden toothpicks, which were put into the soil according to an alternating grid based pattern, leaving 10-mm space between the plastic cover and the soil surface. The use of plastic sheets with known dimensions allowed an accurate measurement of the covered area and better resembled the natural cover of germinating plants compared to mesh screens. The following surface covers were applied: 0, 0.25, 0.40, 0.50, 0.75, and 0.9 m2 m–2. The kinetic energy of the rain on the soil surface, Ek (J m–2 s–1), was calculated by multiplying the kinetic energy of the rain (J m–2 s–1) by the fraction of uncovered soil surface (m2 m–2).
The experiments were done using two different kinds of rainfall simulators to have a wider range of rainfall intensities and rainfall kinetic energies. In the first series of experiments, a rainfall simulator with capillary tubes was used. The water level above the capillary tubes determined the rainfall intensity. Two different water levels were used, resulting in intensities of 25 ± 5 and 40 ± 5 mm h–1, respectively. The diameter of the raindrops was 4.92 ± 0.12 mm and their fall height was 2.75 m. The nomograph constructed by Laws (1941) was used to estimate the terminal velocity of the drops. The calculation of the rainfall kinetic energy was based on the normal velocity of raindrop impact. The rainfall intensities of 25 and 40 mm h–1 provided rainfall kinetic energies of 0.15 and 0.23 J m–2 s–1, respectively. More details about the rainfall simulator were given by Gabriels et al. (1973). Rainfall intensities were kept constant during the experiments, which lasted for 100 min. The first runoff sample was taken after 20 min; afterward, sampling was done every 10 min. Runoff volume was measured and analyzed for sediment content. During the second half of the experiment, a steady runoff and sediment discharge was obtained. Rainfall simulations at 40 mm h–1 indicated that the ER of the suspended sediment reached a steady state when sediment discharge became constant (Schiettecatte, 2006). Therefore, only the runoff sampled at the end of the experiment was used to analyze the OC content in the sediment.
To provide higher rainfall intensities (and consequently higher runoff discharges), the rainfall simulator in the ICE wind tunnel (International Centre for Eremology, Ghent University, Ghent, Belgium) was used. A detailed description of the wind tunnel was given by Gabriels et al. (1997) and Cornelis et al. (2004a). Instead of using capillary tubes, the rainfall was simulated by pressurized nozzles (type TeeJet TG SS 14w), fixed at a height of 1.8 m. A nozzle pressure of 1.5 x 105 Pa was applied. The kinetic energy was determined indirectly using the splash-cup technique (Ellison, 1947; Salles and Poesen, 2000; Salles et al., 2000; Cornelis et al., 2004b,c). Because the median diameter of the raindrops produced by the nozzles was only 1 mm, the kinetic energy of the raindrops was much lower than the simulator with the capillary tubes. The rainfall simulation experiments lasted for 30 min, with a rainfall intensity of 70 mm h–1 and a rainfall kinetic energy of 0.035 J m–2 s–1 for conditions in which one nozzle was used. When two nozzles were used (at 1-m distance from each other), an intensity of 120 mm h–1 with a rainfall kinetic energy of 0.079 J m–2 s–1 was applied for 20 min. Every 2 min, the runoff volume was measured and analyzed for sediment content. The runoff sampled at the end of the experiment was used to determine the OC content in the sediment.
The combination of two slopes, six different covers, and four different intensities resulted in a total of 48 experiments.
Rill Erosion Experiments
Rill Experiment 1: Addition of Water to a Preformed Erodible Gully
A soil pan (0.96 m long, 0.56 m wide, and 0.20 m deep) was filled with air-dried soil aggregates (<8 mm) to obtain a bulk density of 1.33 Mg m–3. An artificial V-shaped rill was made with slopes of 30% at both sides. The soil pan was placed at slopes of 10, 20, and 30%. Sediment-free tap water was added at a constant discharge of 0.01 and 0.017 L s–1 at the upper edge of the rill. Runoff volume was measured and analyzed for sediment content every 60 s for the low discharge and every 30 s for the high discharge. At the end of the experiment, one more runoff sample was taken to determine OC content. Every experiment lasted for 5 or 10 min depending on the high or low runoff discharge. At the end of the experiments, the rill width was measured for calculation of the sediment discharge per unit width.
Rill Experiment 2: Addition of Water and Sediment to a Preformed Erodible Gully
A similar setup (same slopes and soil pan configuration) was also used with dry soil aggregates (<2 mm) added to the water inflow, resulting in an average input sediment concentration of 163 g L–1. Sampling and measurements were done following the same procedure as described above.
Rill Experiment 3: Rainfall Simulation on a Preformed Erodible Gully
An analogous experiment was done using the same soil pan configuration, but neither water nor sediment was added to the rill. Instead, rainfall simulations were done in the ICE wind tunnel, using intensities of 100 mm h–1. Kinetic energy of the rainfall was measured using the splash-cup technique. The soil pan was set at slopes of 10, 20, and 30%. The rainfall simulations lasted for 30 min and runoff was sampled every 30 s. At the end of the experiment, one more runoff sample was taken to determine OC content. At the end of the experiments, the rill width was measured for calculation of sediment discharge per unit width.
Analysis of Samples and Data
The sediment concentration of the runoff samples was determined by drying (at 105°C), using tared recipients. The OC content was determined using the method of Walkley and Black (1934). Because this method yields, on average, 75% of the total OC, the results of the analysis were divided by 0.75 to obtain the total OC content (De Leenheer and Van Hove, 1958). The fully dispersed particle-size analysis was obtained using the sieve pipette method (De Leenheer, 1966). The minimally dispersed particle-size distribution was obtained in a similar way without chemically dispersing the samples.
Aggregate stability was determined by wet sieving according to the procedure of De Leenheer and De Boodt (1959), using a fixed aggregate-size distribution with a mean weight diameter (MWD) of 4.45 mm. Briefly, the soil aggregates were prewetted by raindrops and incubated for 24 h. After wet sieving for 5 min, the MWD of the resulting aggregate size distribution was determined. The stability index was calculated as the reciprocal value of the difference between the original MWD (4.45 mm) and the MWD of the wet-sieved soil sample. After wet sieving and drying (at 105°C), the OC content of the different aggregate size fractions (0.3–0.5, 0.5–1.0, 1.0–2.0, 2.0–2.8, 2.8–4.8, and 4.8–8.0 mm) was determined.
Unless stated otherwise, the ER of a soil constituent was calculated by dividing the content of the constituent in the suspended sediment by its content in the original soil material. Analysis of variance and regression analysis were done using the statistical software SPSS Version 11.0.1 (SPSS, 2001).
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
![]() | [1] |
In conducting laboratory rainfall simulations on small soil pans (0.3 by 0.45 m) but with a rainfall intensity of 30 mm h–1, Jacinthe et al. (2002) found, for undisturbed samples from moldboard-plowed plots, EROC values close to 1.7 for qs rates of 3.7 and 4.9 mg s–1 m–1. Although these qs values are outside the range of our observations, the EROC values measured by Jacinthe et al. (2002) seem to be slightly lower than expected based on our observations.
Several (usually logarithmic) relationships have been given to predict ER values from soil loss data (e.g., Menzel, 1980). These are difficult to compare with Eq. [1] and the results in Fig. 1, however, as soil loss is usually expressed in kilograms per hectare. Rescaling soil losses from experimental plot areas to hectares may produce different results because of differences in plot width. Moreover, soil loss expressed in kilograms per hectare does not indicate the intensity of the erosion process because no time component is included. Low-intensity rainfall for a long period may cause as much erosion as a short, intense rainfall event, but enrichment processes will be different in the two cases. Weigand et al. (1998) found that for the same soil loss of 60 kg ha–1, the ER values of 137Cs varied between 0.51 and 2.46. Therefore, the applicability of an ER in models that assess soil loss on a long-term (annual) basis can be questioned, especially if the runoff events for a longer period show large variations in erosion intensities.
The reason for using unit sediment discharge, qs (kg s–1 m–1), which is the sediment transport per unit width and unit time, rather than soil loss (kg ha–1), in Eq. [1] is that it is a better expression for "erosion intensity" and can also be measured easily during experiments. Furthermore, qs is commonly used in sediment transport equations, relating runoff to erosion (e.g., Nearing et al., 1997). Equation [1] enables linking nutrient losses to runoff and sediment transport, which is of particular interest in the case of event-based modeling.
Despite the difficulties in comparing data due to the different expressions of soil loss, we verified our results with a relationship given by Sharpley (1985), who related the EROC to soil loss (kg ha–1) for different soil types. For a Durant soil, which has similar textural properties as the soil in our study, he found
![]() | [2] |
![]() | [3] |
To further validate our findings, we used data published by Polyakov and Lal (2004b). They conducted rainfall simulations on small soil pans (0.3 m wide and 1 m long) filled with sieved (<8-mm) silt loam soil. Four soil pans were positioned in a cascade system in which runoff from the upper soil pan flowed into the lower pan. A subsample of the runoff from each soil pan was taken at regular time intervals and analyzed for sediment and OC content. The subsequent pans were put under different slopes (8, 7, 1, and 0.5%), resembling a concave hillslope. This setup resulted in erosion on the 8 and 7% slopes, of which the sediment was partly deposited on the 1% slope. Because deposition is a different enrichment process than interrill erosion, we only used the data from the 8 and 7% slopes to validate our findings. As shown in Fig. 2 , the observations of Polyakov and Lal (2004b) agree well with ours: the best-fit power function for the data in Fig. 2 closely resembles Eq. [1]. Polyakov and Lal (2004b) further found that the EROC decreased with rainfall duration. They attributed the higher EROC values during the initial stage of the rainfall experiment to a flushing of loose organic particles and fine soil fractions. Based on Fig. 2, we may conclude that the higher EROC values are a result of the lower qs values. Due to the lower transport capacity of the runoff, a selective transport of small particles rich in OC occurs.
|
Ghadiri and Rose (1991a) found that "stripping" of soil aggregates was the main mechanism in nutrient enrichment of eroded particles in the case of a well-aggregated soil. This was also confirmed by Wan and El-Swaify (1998) in the case of a well-aggregated Oxisol. The soil used in our experiments has a stability index of 0.36, indicating a low aggregate stability. This is also confirmed by the particle size data (Table 1) showing a similar distribution for minimally dispersed and fully dispersed samples, except for the clay fraction, which was mainly aggregated to 2- to 10-µm sizes in the case of the samples that were minimally dispersed. The use of air-dried soil in the experiments may have enhanced the rapid breakdown of the aggregates, as dry aggregates are more susceptible to slaking and breakdown by rainfall impact than prewetted aggregates. Therefore "stripping" is considered to be an unimportant enrichment mechanism for this soil. If "stripping" would have been important, the EROC should decrease with increasing Ek because of a more complete breakdown of the aggregates. This trend could not be observed in our results (Fig. 1), however, and Ek was found to be an insignificant variable in the stepwise regression analysis. Analysis of the OC content in the different aggregate fractions obtained after wet sieving did not show a significant increase in OC (at the 0.05 level) with decreasing aggregate size (Table 2 ). This indicates a rather homogeneous distribution of OC within the aggregates. This was also observed by Schiettecatte (2006) for the OC content and textural composition of different aggregate classes, which had been separated by dry sieving the same soil used in our study. Therefore, we can conclude that for the soil under study, detachment by raindrop impact is an unselective process. This agrees with the results of Gabriels and Moldenhauer (1978), who found that for different soil types (fine sand, silt loam, silty clay loam, and silty clay), the particle size distribution of detached material by splash erosion was not significantly different from that of the original soil.
|
The results given in Fig. 1 further indicate that detachment by overland flow in rills is not influenced by raindrop impact. This is in agreement with the results of the interrill erosion experiments.
| CONCLUSIONS |
|---|
|
|
|---|
| NOTES |
|---|
|
|
|---|
Received for publication June 4, 2007.
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Crop Science | |||
| Journal of Natural Resources and Life Sciences Education |
Vadose Zone Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||