|
|
||||||||
a INRA-Science du sol, B.P. 20619, 45166 Olivet Cedex, France
b USDA-ARS, National Soil Erosion Research Lab., 275 S. Russell St., Purdue Univ., West Lafayette, IN 47907-2077
* Corresponding author (Frederic.Darboux{at}orleans.inra.fr)
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: WEPP, Water Erosion Prediction Project
| INTRODUCTION |
|---|
|
|
|---|
This kind of roughness effect has been incorporated into erosion assessment tools such as Universal Soil Loss Equation (USLE) and its revised version (Revised USLE or RUSLE) (Renard et al., 1997). This commonly accepted roughness scenario compounds the runoff production into sediment production, that is, the reduced erosion is caused by a reduced water runoff. Hence it does not differentiate between the roughness effect on water runoff and the roughness effect on sediment production.
Even in a process-based model, such as the Water Erosion Prediction Project (WEPP), where the process of runoff production supposedly has been isolated from the sediment production, an increased surface roughness also results in an overall reduction in sediment delivery. In WEPP, an increased surface roughness causes a decrease in interrill sediment delivery and an increase in critical shear resistance in the rills (Flanagan and Nearing, 1995).
Despite the dominance of research results and predictive models showing that an increased roughness decreases erosion, there is evidence pointing the other direction. Burwell et al. (1968) and Burwell and Larson (1969) showed that after runoff had initiated, a rougher surface might not have the distinctly higher infiltration as a smooth surface as shown before runoff. The laboratory study of Helming et al. (1998) showed that while runoff was marginally affected, rough surfaces did show a greater soil loss than smooth surfaces because flow concentration may cause a localized increase in erosion. On the other hand, surface depressions that trap sediment and surface mounds that increase flow meandering (or resistance) may lead to a reduced sediment delivery. Therefore, the net roughness effect on sediment delivery depends on the balance between these opposing processes, and erosion can either increase of decrease as soil roughness is increased.
The effect of surface roughness on runoff was often associated with surface storage capacity, that is, volume of water puddles (Mitchell and Jones, 1976; Moore and Larson, 1979; Onstad et al., 1984; Moran and Vézina, 1993; Hansen et al., 1999; Kamphorst et al., 2000). Conceptually, the rainfall-runoff process can be divided into three stages and surface roughness may affect each of them. Stage 1 is mainly for surface wetting and depression filling and ends when runoff starts at the point of observation. Time to runoff is usually used to characterize this stage. Stage 2 is mainly associated with the rising portion of the hydrograph as the runoff contributing area expands. At Stage 3, runoff reaches a plateau or an apparent steady state when the entire surface is contributing runoff. Prior research indicates the importance of quantifying soil erosion at apparent steady-state runoff when the full detachment and transport potentials have been reached (Huang, 1998; Huang et al., 1999; Zheng et al., 2000). However, most prior roughness studies considered a fixed amount of rain and compared the runoff and sediment productions from rough vs. smooth surfaces. Such an approach lumps the effects of at least the first two runoff stages and does not allow a comparison of the true roughness effect under full runoff, or Stage 3. Since sediment production is closely linked to runoff production, factors affecting runoff generation during a rainfall event, such as initial wetting, depressional storage filling and infiltration, need to be isolated first before a meaningful comparison in sediment production can be made.
The roughness effect on erosion can be further compounded by surface and subsurface factors affecting soil erosion, because erosion process itself also causes a change in surface morphology or microtopography. Recently, the near-surface hydraulic gradient, that is, drainage and seepage, has been shown to significantly affect erosion (Bryan and Rockwell, 1998; Huang and Laflen, 1996; Huang, 1998; Owoputi and Stolte; 2001). How the roughness effect interacts with the hydrologic condition in runoff and sediment production has not yet been quantified.
Soil surface roughness is usually partitioned into oriented roughness and random roughness (Römkens and Wang, 1986). In previous studies, changes in random roughness were mainly from changes in aggregate-size distribution. The objective of this study was to compare runoff and sediment productions from two types of roughness, that is, a smooth surface and a surface with mound-and-depression pattern, under different near-surface hydraulic gradient and surface flow conditions. Runoff samples were collected and analyzed from side-by-side, smooth vs. rough plots under simulated rainfall with and without run-on. Runs made with rain only simulate the conditions prevailing at the upper boundary of a hillslope where there is no upstream contributing area. By adding run-on, we aim to reproduce the flow conditions encountered along a hillslope. The present approach studies separately the effect of roughness on the duration of the first stage (or time to runoff initiation or runoff delay) and on the water and particle fluxes at apparent steady state, that is, during Stage 3. An analysis of the roughness effect under a range of surface and hydrologic conditions provides a better understanding on how soil roughness actually affects the runoff and particle production.
| MATERIALS AND METHODS |
|---|
|
|
|---|
The experimental setup consisted of two soil boxes up and down slope to each other that could be either run independently or connected together. The upslope feeder box was used to vary the inflow to the downslope study box (Fig. 1). Each box had separate rainfall simulators mounted above, thus, enabling us to rain simultaneously on both boxes with different rainfall intensities.
|
A baffle plate was set up 1 m from the outlet of the study box to reduce the edge effects due to excessive seepage at the lower part of the box. Metal plates inserted 10 cm into the soil surface divided the study box along its length into two separate plots each 0.6 m wide and 5 m long. This arrangement allowed us to prepare and make rain events on a pair of contrasting smooth and rough surfaces simultaneously. The metal divide extended 5 cm above the soil.
The surface roughness was measured with an instantaneous-profile laser scanner (Darboux and Huang, 2003) on the lower 3.9-m portion of the study box with a horizontal resolution of 1.5 mm and a vertical resolution of 0.5 mm. Surface storage capacity was computed from the laser scanner data using the algorithm developed by Planchon and Darboux (2001).
Experimental Procedure
Box Preparation
The boxes were initially filled with air-dried surface soil. Before the series of experiment was started, seepage and drainage conditions were alternated to stabilize the soil structure.
Before each experiment, box preparation started with air drying of the soil surface using a fan. After the soil surface appeared dry, approximately 5 cm of the surface soil was turned using a hand trowel to help additional drying. Aggregates bigger than 5 cm were manually broken down using the hand tool. During box preparation, new soil was added to compensate for soil loss from the previous experiment to keep a similar amount of soil in both feeder and study boxes.
The prepared surface resembled a fine seedbed with no aggregates larger than 1 cm. The whole soil surface was smoothed down to obtain an even surface. The metal plates were then inserted to divide the study box in two equal areas. On one side of the study box, the surface was kept smooth. On the other side, depressions were molded by hand (Fig. 2). Depressions had a circular shape with 10 to 12 cm in diameter, a depth around 2 cm, and a density of approximately 40 depressions per squared meter.
|
Experiment
A total of eight experiments were conducted (Table 1). The slope of the study box was kept constant at 5% and the feeder box was set to 10%. Experiments were conducted with either seepage with water level at the watering trough maintained 5 cm above the soil surface or under free drainage condition. Experiments with seepage and drainage conditions were alternated to avoid a long-term evolution of soil properties during the experiment series. The sides with smooth surface and surface with depressions were also alternated to avoid a systematic bias due to potential differences in lateral conditions.
|
A rain event consisted of a sequence of operations performed while rainfall was being applied to each of the two sides of the study box. Soil roughness was measured by a laser scanner before each rain event. Rain was first applied on the right side of the study box, with the left side protected from the rain with corrugated metal sheets (Fig. 3). Initially, the feeder box was disconnected from the study box (no runonStage A in Fig. 3). Runoff samples from the study box were collected in 1-L bottles at 1-min time step. Depending on the flow rate, samples were collected for the full 1 min or up to three-fourth full if the bottle would overflow (in this case, collection duration was recorded). After the runoff flux from the study box reached an apparent steady state, defined as three successive bottles with weight differences of <50 g, eight samples were taken simultaneously at the outlets of both boxes. Then, the outlet of the feeder box was connected to the upslope end of the study box (with runonStage B in Fig. 3). After the apparent steady state was reached at the outlet of the study box, another eight samples were collected. Then the boxes were disconnected (no runonStage C in Fig. 3) and, after few minutes, four samples were collected again at each outlet to check the similarity of flow rates before and after the connection.
|
Runoff samples were weighted at collection time. After the end of an experiment, 3 to 5 mL of saturated alum [AlK(SO4)2] was added to each 1-L sample bottle to flocculate the solid fraction. The next day, clear supernatant was poured off and bottles were oven dried at 105°C. The dry bottles were weighted and water and sediment masses were calculated by subtracting the bottle tare weight. Runoff and particle fluxes were calculated from these data and adjusted for the duration of sample collection.
Statistical Analysis
We used the R Statistical Software (R Development Core Team, 2004) to analyze the results. To estimate the statistical significance of differences in time to runoff initiation between smooth surfaces and surfaces with depressions, paired t tests were computed for each of the subsurface conditions (drainage and seepage). The null hypothesis H0 was the equality of the differences to zero. The alternative hypothesis was that the time to runoff was larger for the surfaces with depressions than for the smooth surfaces. Similar statistical procedures were also used for differences between apparent steady state fluxes at Stages A and C. We considered a Type-I statistical risk of 5%.
To identify the conditions where roughness had a significant effect on runoff at the outlet of the study box, an analysis of variance was conducted. The three variables were water flux, particle flux, and particle concentration. The full linear model was:
![]() |
To better identify the conditions where the roughness had a significant effect, submodels were run on subdatasets. First, the dataset was split according to the subsurface condition and then each of these subsets was split according to the upstream flow condition (with runon or without runon). We computed Type-III sums of squares (also known as Yates' weighted squares-of-means) and considered a Type-I statistical risk of 5%. Postulates of linear modeling (residuals with no bias, equality of residual variances, independence of residuals and Gaussian distribution of residuals) were checked graphically.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
Storage Capacity
Initial storage capacity was clearly a function of surface condition, and for all surfaces it decreased with added rainwater (Fig. 4). Differences between initially smooth surfaces and initial surfaces with depressions continued up to the second rain, at least. The depressions appeared to be more persistent for drainage than for seepage condition. The faster decline in storage capacity under seepage condition could be caused by both greater soil erodibility and runoff rates (Huang and Laflen, 1996; Huang, 1998).
|
|
To separate the relative effects of storage and infiltration under drainage conditions, a time to depression filling was computed from roughness data and compared with the actual time to runoff initiation. This hypothetical time to depression filling was estimated by dividing the storage capacity by the area of the surface and the rainfall intensity. This assumes that there is no infiltration and that runoff starts only after all depressions are filled. The percentages of time to runoff initiation explained by the time to depression filling are displayed in Table 3. They were statistically significant for the first two rain events (Table 3). For the initially smooth surfaces, water storage by surface depressions explains in average only 8% of the runoff delay for the first rain and 2% for the second rainfall. In the cases with initial depressions, surface storage explains, in average, 50% of the time to runoff initiation for the first rainfall and 70% for the second rainfall.
|
Although a single rainfall intensity was used in this research, it must be noted that the proportion of delay directly explained by depressional storage capacity is also a function of rainfall intensity. For lower rainfall intensities (but identical infiltrability and storage capacity), it is expected that more water would infiltrate before depressions get filled. In this case, calculations would show smaller direct effect of depressions on time to runoff initiation. On the other hand, indirect effects of depressions would be larger because depressions trap water and prevent runoff, hence allowing more water to infiltrate. The relative balance between the direct and indirect effects of depressions on runoff initiation remains to be specified.
Fluxes and Concentration at Apparent Steady State
For each rain event and surface condition, two types of apparent steady state were reached at the outlet of the study box. The first type of apparent steady state was reached when the boxes were not connected together (Fig. 3). There was no run-on inflow to the study box. Runoff water was due to the applied rain and to seepage flow when it was applied. For a given surface and a given rain event, this no-runon apparent steady state was reached once before the two boxes are connected together (Stage A) and once after they have been disconnected (Stage C). By comparing the fluxes for these two disconnected stages, we could assess the quality of the no-runon apparent steady state. Statistical analysis showed both water and particle fluxes were not statistically different before and after connection. Consequently, flux data from Stages A to C were averaged and analyzed as a single no-runon steady state data set (Tables 4, 5, and 6).
|
|
|
The experiment was designed to analyze the effect of soil surface roughness on various flow conditions after an apparent steady state had been reached. Performances of the linear models for nested datasets are outlined in Table 7. All of these models are significant at the 5% probability level, but not all the factors are significant at that probability level.
|
|
At apparent steady state, roughness had a significant effect on water flux except for the condition "drainage with runon" (Table 8). The surface with initial depressions yielded larger water flux than the smooth side (Table 4). Nevertheless, it should be pointed out that the difference in fluxes between the two roughness conditions was mostly in the range of 10% (Table 4). The effect of roughness on particle flux and concentration was mostly nonsignificant. If such effect exists, its amplitude is probably very low and could not be characterized with the current dataset.
Two or three successive rainfalls were applied on each surface. For a given surface, the main difference between the initial conditions of each rain event was soil surface properties. So, the factor "rain event" could be used to assess the evolution of the roughness effect with successive rainfalls. Among successive rainfalls, the upstream input of sediment was varied, possibly modifying the sediment output. The factor "rain event" is significant on water flux, particle flux, and particle concentration only when the whole dataset is considered, but no conclusion can be drawn about the effect of the successive rain events because values of "rain event" coefficients lack statistical significance. Overall, it appears that successive rains did not change much the runoff characteristics and that effect of initial roughness was larger that the effect of successive rains.
Overall, the analysis points out initial depressions had a continuing effect by increasing the water flux at apparent steady state. In the meantime, the storage capacities of the surfaces with initial depressions decreased sharply with the successive rain events. This leads to the rejection of any hypothesis connecting the higher runoff fluxes with the storage capacity of the surfaces. Therefore, the decrease in infiltration should be related to the initial presence of depressions and not to their persistence.
Results of the present study are not in total agreement with previous experiments reported in the literature. Prior studies showed that an increased roughness either decreased water runoff (Johnson et al., 1979; Cogo et al., 1984) or had no significant effect on it (Burwell et al., 1968; Burwell and Larson, 1969; Helming et al., 1998). With regard to particle transfer, it was found an increased roughness either decreased soil loss (Johnson et al., 1979; Cogo et al., 1984) or increased it (Helming et al., 1998).
The reason for these different results probably lies in differences in soil properties, roughness characteristics, and experiment setups. In the current experiment, aggregates had small diameters and their size was kept somewhat identical for all roughness conditions. The roughness of the surface was dominated by the macro-scale depressions that our experiment specifically targeted. The other studies varied the aggregate size (either by tillage or by sieving), also varying the size of the pores open to the surface. The larger the aggregates were, the larger the pores, probably explaining the decrease in runoff rate with the increase in aggregate size. Such phenomenon was not likely to occur in our experiment.
| CONCLUSIONS |
|---|
|
|
|---|
Considering the steady-state water flux is uncommon in field conditions, and even at this stage, the effect of surface depressions is limited to a slight increase in water flux, we conclude that the only assured soil and water conservation benefit from soil surface roughness is the delay in runoff initiation before the entire field is fully saturated and contributing to runoff.
In this study, the relationship between roughness, infiltration and runoff initiation was assessed for a single rainfall intensity and proved to be significant. To achieve a better implementation of this relationship in water erosion models, the present results will need to be extended for a range of rainfall intensities.
In future experiments, it may be important to partition random roughness into subcomponents such as aggregate size and mound-and-depression pattern to better explain roughness effect on runoff and erosion. Comparison of the present results reveals the complexity of the interaction between roughness, overland flow and erosion. At this point, our knowledge is still insufficient to offer a mechanism for the roughness effects on overland flow and sediment detachment and transport.
| ACKNOWLEDGMENTS |
|---|
Received for publication December 2, 2003.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
L. R. Ahuja, L. Ma, and D. J. Timlin Trans-Disciplinary Soil Physics Research Critical to Synthesis and Modeling of Agricultural Systems Soil Sci. Soc. Am. J., February 2, 2006; 70(2): 311 - 326. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Crop Science | |||
| Vadose Zone Journal | Journal of Plant Registrations | ||||
| Journal of Natural Resources and Life Sciences Education |
Journal of Environmental Quality |
||||