|
|
||||||||
a State Key Laboratory of Soil Erosion and Dryland Farming on Loess Plateau, Institute of Soil and Water Conservation, CAS and MWR, Northwestern Sci-Tech University of Agriculture and Forestry, 26 Xinong Road, Yangling, Shaanxi, 712100 China
b USDA-ARS National Soil Erosion Research Laboratory, 275 South Russell Street, Purdue University, West Lafayette, IN 47907-2077
* Corresponding author (flzh{at}ms.iswc.ac.cn)
Received for publication December 22, 2002.
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: C, control fertilizer treatment (no fertilizer input) L, low fertilizer treatment (40 kg P ha1, 100 kg N ha1) H, high fertilizer treatment (80 kg P ha1, 200 kg N ha1) FD, free drainage hydraulic gradient Sa, saturation hydraulic gradient Sp, artesian seepage without rain hydraulic gradient Sp + R, artesian seepage with rain hydraulic gradient
| INTRODUCTION |
|---|
|
|
|---|
Due to its high solubility, NO3N tends to be transported in drainage and subsurface flow. In southern Georgia, Hubbard and Sheridan (1983) reported that 20% of the applied N over a 10-yr period was lost in surface runoff and subsurface flow, and 99% of this loss occurred in subsurface flow. Other studies by Hubbard et al. (1991) and Lowrance (1992) at the same watershed and elsewhere by Alberts and Spomer (1985) with different crops showed the same trend with significantly greater NO3N movement in subsurface flow than in surface runoff.
Because P is strongly bound in soil and much less mobile than N, many efforts have been focused on relating different soil test P values to P in surface runoff (Daniel et al., 1993; Pote et al., 1996, 1999; Sharpley et al., 1996; Cox and Hendricks, 2000). Despite a general positive correlation between soil test P to runoff P, the slope of the trend line varied with methods of P extraction, soil type, organic matter, and soil management. Recent studies had been initiated to couple watershed hydrology and chemical transport using the variable source area (VSA) concept (Gburek and Sharpley, 1998; Gburek et al., 2000). Gburek and Sharpley (1998) studied P loss in east-central Pennsylvania and showed that zones of runoff production, and consequently, the areas that ultimately controlled most P transport, were the near-steam saturated areas of the watershed. Further analysis showed that most in-stream P came from soils within 60 m of the stream, rather than from the entire area of the watershed (Gburek et al., 2000).
In this study, we hypothesized that because most agricultural chemicals are in the surface layer, the near-surface hydraulic gradients could have a large influence on chemical transport. This study was initiated based on recent laboratory findings of increased sediment delivery and rilling from soil subject to saturation and artesian seepage or exfiltrating through flow (Huang and Laflen, 1996; Gabbard et al., 1998; Huang et al., 1999; Zheng et al., 2000). Artesian seepage is a common occurrence during the wet season such as early spring after a wet winter, when the shallow ground water emerges at the soil surface as return flow in the middle or lower portions of the hillslope (Fig. 1) . Soils with clayey subsoil or tight plow pan are prone to produce the return flow or seep (Whipkey and Kirkby, 1978; Dunne, 1978).
|
The objective of this study was to quantify the effect of near-surface hydraulic gradients on NO3N and PO4P transport in surface runoff water. A laboratory rainfall simulation study was conducted using a silt loam soil under four different near-surface hydraulic conditions: free drainage (FD), saturation (Sa), artesian seepage without rain (Sp), and artesian seepage with rain (Sp + R).
| MATERIALS AND METHODS |
|---|
|
|
|---|
Experimental Setup
The study was conducted on soil pans that were 45 cm long, 32 cm wide, and 35 cm deep. Each soil pan had six drainage holes at the bottom. A water supply system was designed to supply water to the soil pan from the bottom to control the near-surface hydraulic gradient (Fig. 2)
. A saturation condition was created when the supply water level was set at the soil surface and an artesian seepage condition was created when the supply water level was set higher than the soil surface, forcing water to flow out of the soil. In this study, the artesian seepage condition was created with the supply water level set 20 cm above the soil surface. For the free drainage treatment, the water supply tubes were not connected to let the soil pan drain freely under gravity.
|
Experimental treatments in this study included three rates of fertilizer input: control (no fertilizer input), low (40 kg P ha1, 100 kg N ha1), and high (80 kg P ha1, 200 kg N ha1), hereafter designated as C, L, and H, respectively. These fertilizer treatments were subjected to four near-surface hydraulic gradient treatments: free drainage (FD), saturation (Sa), artesian seepage without rain (Sp), and artesian seepage with rain (Sp + R). The detailed experimental treatments appear in Table 1. For each fertilizer level, three replicates were made for rainfall simulation for the FD and Sa treatments, and three additional soil pans were prepared and used to collect the before-run soil profile samples three days after prewetting under free drainage; four replicates were made for the Sp and Sp + R treatments, and four additional soil pans were prepared and used to collect the before-run soil profile samples three days after prewetting under free drainage.
|
After soil pan preparation, a prewetting rain (15 mm of rain at 9.5 mm h1 intensity) was applied at 0% slope. The 9.5 mm h1 rain intensity was selected to avoid the initiation of runoff during the prewetting event. To further avoid raindrop-induced splash and surface sealing, two layers of wire mesh were placed over the soil surface. This initial rain was applied to allow fertilizer leaching into the soil profile, create a uniform surface soil moisture condition before the experiment, and reduce surface variability from preparation. After the prewetting rain, the soil pan was covered with a plastic sheet and allowed to equilibrate under free-drained conditions for three days.
Rainfall Experiments
Three days after the prewetting rain, soil pans were set to 5% slope and subjected to the experimental hydraulic gradients. A simulated rainstorm of 50 mm h1 for 90 min was applied to the FD, Sa, and Sp + R treatments. For the Sa and Sp + R treatments, the rainfall was applied after the surface was either saturated or when seepage flow started. During rainfall simulation runs, runoff samples were collected at 10-min intervals. For the Sp treatment, due to the small quantity of seepage flow, runoff samples were collected at 15-min intervals for 90 min immediately after the seepage flow started. During each run, the rainfall amount was measured at least twice with a hyetometer on the right and left sides of the soil pan.
Immediately after each run, runoff samples were centrifuged (4000 rpm, 15 min) and filtered through 2.5-µm paper (Whatman [Maidstone, UK] no. 5). The filtered solution was stored in a refrigerator and analyzed 24 h later. Sediment samples at the bottom of the centrifuge bottle were washed into cans, and oven-dried at 55°C. This dry sediment was weighted to calculate sediment delivery.
Collection of Soil Samples from the Soil Pans
Soil samples taken before and after each run were analyzed for moisture and nutrient contents. The before-run soil profile samples were taken three days after prewetting under the free drainage. The after-run samples were taken one day after the runs were completed. Samples were taken at five depth intervals: 0 to 2, 2 to 5, 5 to 10, 10 to 15, and 15 to 20 cm. To ensure a proper representation of the profile distribution, soil moisture samples were collected from five different locations and nutrient samples from 12 different locations in each soil pan. Soil moisture samples from a given depth increment in a soil pan were combined in an aluminum can and oven dried at 105°C. Soil samples for determining soil nutrient concentration from a given depth increment in a soil pan were combined, air-dried, sieved (2 mm), and stored for chemical analysis.
Chemical Analysis
Concentrations of NO3N in surface runoff water were determined by the cadmium reduction method, and PO4P in surface runoff water was measured by the ascorbic acid method (Sparks, 1996). The P in soil was extracted using the BrayKurtz P-1 method and determined using the ascorbic acid method. Soil NO3N was extracted with KCl solution and measured by the cadmium reduction method (Sparks, 1996).
Calculation of Dissolved Nitrogen and Phosphorus Loss in Surface Runoff Water
For each run, NO3N and PO4P loss in surface runoff water for the event, Ln (mg), was calculated by the following formula:
![]() |
For each run, the average NO3N and PO4P concentration in surface runoff water, Cn (mg L1), was calculated as:
![]() |
For each hydraulic gradient treatment under the same fertilizer treatment, the average NO3N and PO4P concentration in surface runoff water, Cm (mg L1), was calculated as:
![]() |
Determination of Statistical Significance
Within each fertilizer treatment, the LSD test was used to determine whether differences in runoff, sediment, concentrations, and losses of NO3N and PO4P among the four hydraulic gradients were statistically significant at the 95% confidence level. Similarly, within each hydraulic gradient treatment, the LSD test was performed to determine whether differences in the concentrations and losses of NO3N and PO4P among the three fertilizer treatments were statistically significant.
| RESULTS AND DISCUSSIONS |
|---|
|
|
|---|
Rainfall, Runoff, and Sediment Delivery
Rainfall ranged from 67.8 to 80.8 mm (Table 2). Variability of 10% from the target intensity is expected for this type of rainfall simulator. Runoff from the Sa treatment was on average 10 mm greater than the applied rainfall. This could have been due to the release of stored soil water by pressure wave advance (Torres, 2002). This physical process may influence the chemical transport process.
|
Concentrations of Nitrogen and Phosphorus in Surface Runoff Water
The average concentrations of NO3N and PO4P in surface runoff water from the FD treatments were less than 0.3 mg L1 regardless of the fertilizer treatments (Table 3). Except for PO4P concentrations under the control treatment, when the near-surface hydraulic gradients shifted from FD to Sa, then to Sp, concentrations of NO3N and PO4P increased greatly. For the low and high fertilizer treatments, NO3N concentrations in surface runoff water from the Sp and Sp + R treatments were statistically greater than those from the Sa and FD treatments. Under the low fertilizer treatment, NO3N concentrations in runoff water from LSp and LSp + R averaged 6618 and 897 times greater than those from LFD, and 241 and 33 times greater than those from LSa, respectively. Under the high fertilizer treatment, NO3N concentrations in surface runoff water from HSp and HSp + R averaged 11665 and 1509 times greater than those from HFD, and 127 and 16 times greater than those from HSa, respectively. The NO3N concentrations from Sp and Sp + R were also significantly different. Due to rain water dilution, NO3N concentration from the Sp + R treatment decreased by 87% compared with the Sp treatment under the low and high fertilizer treatments.
|
Fertilizer application rate greatly influenced NO3N and PO4P concentrations in surface runoff water. Except for the FD treatment, concentrations of NO3N and PO4P in surface runoff water from Sa, Sp, and Sp + R under the high fertilizer treatment were statistically greater than those under the low fertilizer treatment (Table 4). Doubling the fertilizer rate increased NO3N and PO4P concentrations in surface runoff water by 13 and 145%, respectively, under the FD treatment. Under Sa, Sp, and Sp + R treatments, NO3N concentrations with the high fertilizer treatment increased 275, 98, and 89% from the low fertilizer treatment. Compared with the low fertilizer treatment, PO4P concentration in surface runoff under the high fertilizer treatment was 126, 415, and 181% greater under the Sa, Sp, and Sp + R treatments, respectively.
|
Average PO4P loss from CSp + R was 4.7 times greater than loss from CFD. Under the low fertilizer treatment, average PO4P losses from Sp + R and FD were statistically different. The PO4P loss from LSp + R was 9.6 times greater than that from LFD. However, PO4P losses from LSp and LFD were not statistical different. Under the high fertilizer treatment, average PO4P losses from HSa, HSp, and HSp + R were statistically higher than those from HFD. Average PO4P losses from HSa, HSp, and HSp + R were 5.6, 4.0, and 10.3 times greater, respectively, compared with HFD. The PO4P loss from Sp + R was significantly higher than that from Sp under both low and high fertilizer treatments, with a magnitude of 8.0 and 2.9 times greater, respectively (Table 3).
Under the FD treatment, NO3N loss was approximately 0.01% of the applied N. Under the Sp and Sp + R treatments, NO3N loss accounted for 11 to 16% of the applied N (Table 3). Seepage flow alone (9.013.0 mm) comprised only 9 to 13.5% of the total runoff from the Sp + R treatment, but it appeared to produce most of the NO3N loss in surface runoff water. These results show that seepage flow could make an important contribution to total chemical transport.
Temporal Trends of Nitrogen and Phosphorus Concentrations
Under the FD treatment, the NO3N concentrations in runoff were greatest in the initial runoff and decreased gradually during the run (Fig. 3)
. With seepage flow, a gradual increasing trend of NO3N in runoff water was observed as the run progressed. These responses were somewhat expected, because under the free drainage treatment, the downward movement of rain water leached NO3N deeper into the soil profile, while under the seepage treatments, NO3N in the profile was brought to the surface by the seepage flow. The temporal trend under the saturation treatment was between those from the drainage and seepage treatments. Without the fertilizer application, the control condition showed little change in NO3N loss as time progressed for all hydraulic gradients tested. The temporal variation of the PO4P concentration during the 90-min run was less pronounced than the NO3N variation (Fig. 4)
.
|
|
Distribution of Soil Nitrogen and BrayKurtz Phosphorus in the Soil Profile
Concentrations of soil NO3N and BrayKurtz P-1 in the soil profile before and after the run provided additional information for explaining different patterns of NO3N and PO4P losses in runoff water (Table 5). Three days after the 15-mm pre-rain, the highest NO3N concentration was found at the 5- to 10-cm layer under the low and high fertilizer treatments, and at the 10- to 15-cm layer under the control treatment.
|
The P movement in the soil profiles under all hydraulic gradients was less pronounced than NO3N. Soil BrayKurtz P-1 content in the 0- to 2-cm layer under all hydraulic gradients was very high, which could be attributed to fertilizer being applied in the top 2-cm soil layer and little vertical movement of the applied P with rain water. Soil BrayKurtz P-1 in the top 2-cm layer after the 90-min run significantly decreased under the LSa, LSp, LSp + R, HSa, HSp, and HSp + R treatments. Under the low fertilizer treatment, soil BrayKurtz P-1 content in the top 2-cm layer from LSa, LSp, and LSp + R decreased by 35, 53, and 50%, respectively, compared with the soil BrayKurtz P-1 content in the top 2-cm layer before the run. A similar degree of decrease of the soil BrayKurtz P-1 content in the surface layer after the 90-min run was observed for HSa, HSp, and HSp + R under the high fertilizer treatment. This decline appears to be a result of greater PO4P loss rates in runoff from these treatments except the LSp treatment.
| CONCLUSIONS |
|---|
|
|
|---|
The results of this study also challenge current rainfall simulation methods used in the field to quantify soil erodibility and chemical transport. Field studies are often performed when the fields are dry enough to make the runs, hence under a drainage gradient. Results from prior studies on near-surface hydraulic gradient effects on soil erosion and the current study on NO3N and PO4P transport show that the saturation and artesian seepage conditions can cause greater soil loss and chemical transport than the drainage condition. An examination of soil loss and chemical transport data at various soil moisture conditions allows us to identify critical conditions that would cause water quality problems at field scale.
Although extensive efforts have been made to determine factors contributing to P and N losses, the specific effects of near-surface hydraulic gradient have not been previously quantified. Incorporating the basic understanding of hydraulic factors may contribute to more effective control measures that can minimize the chemical loading to surface runoff at both hillslope and watershed scales.
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
Related articles in JEQ:
This article has been cited by other articles:
![]() |
P. J. A. Kleinman, M. S. Srinivasan, C. J. Dell, J. P. Schmidt, A. N. Sharpley, and R. B. Bryant Role of Rainfall Intensity and Hydrology in Nutrient Transport via Surface Runoff. J. Environ. Qual., July 1, 2006; 35(4): 1248 - 1259. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. W. Cox, J. Varcoe, D. J. Chittleborough, and J. van Leeuwen Using Gypsum to Reduce Phosphorus in Runoff from Subcatchments in South Australia J. Environ. Qual., November 7, 2005; 34(6): 2118 - 2128. [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 |
Soil Science Society of America Journal |