Published in J. Environ. Qual. 32:2172-2179 (2003).
© ASA, CSSA, SSSA
677 S. Segoe Rd., Madison, WI 53711 USA
TECHNICAL REPORTS
Landscape and Watershed Processes
Effect of Rainfall Simulator and Plot Scale on Overland Flow and Phosphorus Transport
Andrew Sharpley* and
Peter Kleinman
USDA Agricultural Research Service, Pasture Systems and Watershed Management Research Unit, University Park, PA 16802-3702
* Corresponding author (andrew.sharpley{at}ars.usda.gov).
Received for publication January 6, 2003.
 |
ABSTRACT
|
|---|
Rainfall simulation experiments are widely used to study erosion and contaminant transport in overland flow. We investigated the use of two rainfall simulators designed to rain on 2-m-long (2-m2) and 10.7-m-long (32.6-m2) plots to estimate overland flow and phosphorus (P) transport in comparison with watershed-scale data. Simulated rainfall (75 mm h-1) generated more overland flow from 2-m-long (20 L m2) than from 10.7-m-long (10 L m2) plots established in grass, no-till corn (Zea mays L.), and recently tilled fields, because a relatively greater area of the smaller plots became saturated (>75% of area) during rainfall compared with large plots (<75% area). Although average concentrations of dissolved reactive phosphorus (DRP) in overland flow were greater from 2-m-long (0.50 mg L-1) than 10.7-m-long (0.35 mg L-1) plots, the relationship between DRP and Mehlich-3 soil P (as defined by regression slope) was similar for both plots and for published watershed data (0.0022 for grassed, 0.0036 for no-till, and 0.0112 for tilled sites). Conversely, sediment, particulate phosphorus (PP), and total phosphorus (TP) concentrations and selective transport of soil fines (<2 µm) were significantly lower from 2- than 10.7-m-long plots. However, slopes of the logarithmic regression between P enrichment ratio and sediment discharge were similar (0.2810.301) for 2- and 10.7-m-long plots, and published watershed data. While concentrations and loads of P change with plot scales, processes governing DRP and PP transport in overland flow are consistent, supporting the limited use of small plots and rainfall simulators to assess the relationship between soil P and overland flow P as a function of soil type and management.
Abbreviations:
, volumetric soil moisture DRP, dissolved reactive phosphorus NPRP, National Phosphorus Research Project PER, phosphorus enrichment ratio PP, particulate phosphorus TP, total phosphorus WEPP, Water Erosion Potential Predictor
 |
INTRODUCTION
|
|---|
PHOSPHORUS is an essential nutrient for crop and animal production that can also accelerate freshwater eutrophication (Carpenter et al., 1998; Sharpley, 2000). In the USA, eutrophication is one of the most widespread water quality impairments (USEPA, 1996), with agriculture a primary source of P in the surface waters of many watersheds (United States Geological Survey, 1999). Concern over eutrophication has prompted most states to develop recommendations for land application of P and watershed management based on the potential for P loss in overland flow (USDA and USEPA, 1999). As part of these recommendations, states are assigning soil P thresholds above which the potential for dissolved P loss in overland flow becomes unacceptable. These environmental thresholds are used in site assessment indices to identify agricultural fields that are most susceptible to soil P loss and indicate when current management should be reevaluated (Sharpley et al., 2003).
Despite the widespread implementation of P-based management strategies across the USA, limited field information exists on which to base environmental thresholds of soil P. Consequently, the National Phosphorus Research Project (NPRP) was launched to assess the effects of soil properties, particularly soil test P level, and management on P loss in overland flow (Sharpley et al., 2002b). The NPRP is composed of a consortium of federal and state agencies, as well as land grant universities, representing more than 20 states. To expedite data collection, promote comparability of results, and attempt to maintain field relevancy, studies are conducted on an in situ basis using paired, 2-m2 runoff plots with 2-m flow-path length. Overland flow from the NPRP plots is generated by a portable rainfall simulator based on the operating designs of Shelton et al. (1985) and Miller (1987) and described by Humphry et al. (2002).
In the early 1990s, the Water Erosion Prediction Project (WEPP) used a rotating-boom rainfall simulator (15.2-m-diameter boom) to update soil loss factors in the Universal Soil Loss Equation (USLE) and to improve erosion prediction technology (Simanton and Renard, 1992). Two 36.2-m2 plots with 10.7-m flow-path length are used with the WEPP simulator. The raindrop size, energy, distribution uniformity, and intensities of both NPRP and WEPP simulators have undergone extensive testing and comparison with the properties of natural rainfall (Humphry et al., 2002; Simanton and Renard, 1992). Both NPRP and WEPP simulators use the same nozzles and water pressure to produce rainfall of comparable characteristics.
To date, little research exists on the effect of plot scale on P transport in overland flow. There is, however, a growing body of research on scale-related trends in hydrology and erosion that is relevant to P transport. Wauchope and Burgoa (1995), reviewing pesticide-related research, observed that runoff volume per unit area decreased with increasing plot size, while the length of rainfall required to initiate overland flow increased. Bloschl and Sivapalan (1995), reviewing hydrologic literature, hypothesized that as slope length increases, overland flow processes change from infiltration excess processes to saturation excesses. Le Bissonnais et al. (1998) noted that the selective erosion of fine particles increased with plot size and length.
The main objective of this research was to compare the effects of two rainfall simulators (i.e., WEPP and NPRP) on overland flow patterns (e.g., time to initiation, volume, discharge rate, peak flow), P concentrations (dissolved and total P), and sediment discharge, in relation to plot size and watershed data. Both simulators rely upon the same rainfall generation system, nozzles, and intensity (75 mm h-1), although plot length (2 m) and area (2 m2) of the NPRP simulator is much less than that of the WEPP simulator (10.7 m and 32.7 m2, respectively). Assuming the WEPP rainfall simulator represents field-scale processes reasonably well, due to its large size and history of extensive testing, we attempted to determine if the more portable NPRP simulator can accurately represent overland flowrainfall interaction, soil Poverland flow extraction, P transport processes, and sediment discharge. An additional objective of this study was to establish the impact of scale and geometry of rainfall simulators on their use in examining P transport in overland flow.
 |
MATERIALS AND METHODS
|
|---|
Study Area and Management
Simulated rainfall studies were conducted within the mixed land-use watershed, FD-36, a 39.5-ha subwatershed of Mahantango Creek, a tributary to the Susquehanna River and ultimately the Chesapeake Bay (Fig. 1)
. Three sites in the FD-36 watershed were selected to cover a range in soils, P status, and tillage. The slope of all sites was between 4 and 5%. Sites included a Berks loam (loamy-skeletal, mixed, active, mesic Typic Dystrudept) in orchardgrass (Dactylis glomerata L.); a Berks loam that was tilled prior or seeding with orchardgrass; and a Watson clay loam (fine-loamy, mixed, active, mesic Typic Fragiudult) in no-till corn. The tilled site was chisel-plowed to a depth of about 25 cm and disked to prepare seedbed surface before planting orchardgrass in mid-April 2000. Rainfall simulations were conducted during months in which the majority of overland flow naturally occurs (April and May; Pionke et al., 2000). At the tilled site, simulations were conducted one week after plowing, in late April 2000. Simulations at the grassed Berks and no-till Watson sites (before corn planting) were conducted in early May 2001.

View larger version (17K):
[in this window]
[in a new window]
|
Fig. 1. Location of the study area within the Mahantango Creek watershed in relation to the Susquehanna River and Chesapeake Bay watershed.
|
|
At each site, duplicate adjacent overland flow plots were established for each simulator. Both NPRP and WEPP rainfall simulations were conducted at the same time and on the same day (Fig. 2)
.

View larger version (140K):
[in this window]
[in a new window]
|
Fig. 2. Simultaneous National Phosphorus Research Project (NPRP) and Water Erosion Potential Predictor (WEPP) rainfall simulations on the no-till Berks soil.
|
|
National Phosphorus Research Project Rainfall Simulator Protocol
Overland flow plots, each 1 by 2 m, with the long axis orientated down the slope, were delineated by metal borders installed 5 cm above and below ground level to isolate overland flow. Rainfall was applied with one TeeJet 2HH-SS50WSQ nozzle (Spraying Systems Co., Wheaton, IL) approximately 2.5 m above the soil to achieve terminal velocity (Sharpley et al., 2002b). The nozzle, associated plumbing, in-line filter, pressure gauge, and electrical wiring are mounted on a 3- by 3- by 3-m aluminum frame fitted with canvas walls to provide a windscreen. A coefficient of uniformity of >85% was obtained for rainfall over a 4-m2 footprint, which encompasses one pair of abutting plots. An average rainfall intensity of 75 mm h-1 was applied until 30 min of runoff was obtained. This rainfall intensity and duration has an approximate 10-yr return frequency in south-central Pennsylvania. Local ground water was used as the water source for the simulator and had a DRP concentration of <0.01 mg L-1, total phosphorus (TP) of 0.02 mg L-1, nitrate N of 3.1 mg L-1, pH of 5.7, and electrical conductivity of 0.02 S m-1. Electrical conductivity measurements were made using a digital conductivity meter. Before rainfall, volumetric soil moisture (
) was determined using a capacitance sensor at five locations within a plot (Theta Probe; Delta-T Devices, Cambridge, UK).
Two simulated rainfalls, each of 75 mm h-1 for 30 min of runoff, were applied on consecutive days coinciding with WEPP simulations. Overland flow was collected in metal gutters at the downslope edge of each plot and pumped to 200-L plastic containers. Total overland flow was measured by weighing the containers. A runoff sample was collected from each container after thorough agitation to resuspend and mix sediments.
Water Erosion Prediction Project Rainfall Simulator Protocol
The WEPP simulator is trailer-mounted, with 10 rotating booms (each 7.6 m long) radiating from a central center of rotation, which rotate at about 4 rpm (Simanton and Renard, 1992) (Fig. 2 and 3)
. The arms support 30 TeeJet 2HH-SS50WSQ nozzles positioned at various distances from the stem. The nozzles spray downward from an average height of 2.4 m, apply rainfall at an average intensity of 75 mm h-1, and produce drop-size distributions similar to natural rainfall. Simulator energies are about 77% of those of natural rainfall and the simulator produces intermittent rainfall impulses at the plot surface as the booms pass over the plot. The spatial distribution of rainfall over each plot has a coefficient of uniformity of >90%. Changes in rainfall intensities are produced by regulating the number of open nozzles: 15 nozzles are used to achieve an intensity of 75 mm h-1. Two plots, 3.05 by 10.7 m (32.6 m2), were covered by the simulator (Fig. 3). Before rainfall,
was determined by a capacitance sensor at 15 locations within each plot.
Two rainfall simulation runs were made on each plot pair on consecutive days. Rainfall application rate was measured with a recording rain gauge and rainfall distribution on each plot measured with six nonrecording rain gauges. Plot overland flow was measured by precalibrated 0.12-m HS flumes (4 L s-1 maximum capacity) equipped with water-level recorders, designed to measure small flows with a high degree of accuracy (<1% error; Brakensiek et al., 1979). Continuous hydrographs were produced using the flume's depthdischarge rating table. During a run, times of ponding (half of the plot surface had standing water), runoff initiation, sediment samples, and end of runoff were recorded on field notes for later comparison with recorder charts. Although overland flow samples were collected at 5-min intervals starting 2.5 min after flow initiation, data in this paper reflect composite samples (flow-weighted) of bulked flow to enable direct comparison with NPRP simulator results.
Chemical Analyses
After both WEPP and NPRP rainfall simulations, a minimum of 10 soil samples (05 cm) was collected within each plot, composited, air-dried, sieved (2 mm), and stored for physical and chemical analysis. For all runoff samples, a subsample was immediately filtered (0.45 µm) and stored at 277 K. Filtered samples were analyzed within 24 h of collection and unfiltered samples no more than 7 d after the completion of the rainfall simulation.
Soil particle size analysis was conducted by the hydrometer method after dispersion with sodium hexametaphosphate (Day, 1965). Organic C was determined by combustion using a LECO (St. Joseph, MI) C/N analyzer and pH using a glass electrode at a 1:2.5 soil to water ratio (w/w). Mehlich-3 soil P concentration was determined by 5 min of end-over-end extraction of 1 g soil with 10 mL mixture of 0.2 M CH3 COOH, 0.25 M NH4 NO3, 0.015 M NH4F, 0.013 M HNO3, and 0.001 M EDTA (Mehlich, 1984). Total soil P was determined following digestion with a semimicro Kjeldahl procedure (Bremner, 1996).
The concentration of DRP (subsequently referred to as DRP) in overland flow was determined on the 0.45-µm filtered sample. The concentrations of both total dissolved phosphorus (TDP) and TP were determined on filtered and unfiltered runoff samples, respectively, following digestion with a semimicro Kjeldahl procedure (Bremner, 1996). Particulate phosphorus (PP) was calculated as the difference between TP and TDP. Phosphorus in all soil extracts, filtrates, and neutralized digests was measured by the colorimetric method of Murphy and Riley (1962).
The suspended sediment concentration of each overland flow event was measured in duplicate as the difference in weight of 250-mL aliquots of unfiltered and filtered (0.45 µm) runoff samples after evaporation (378 K) to dryness.
Statistical Analyses
Comparisons of overland flow properties between simulators were conducted by Student's t test, while comparisons between soils were conducted by analysis of variance (ANOVA). Least squares regression was used to describe associations between individual variables, with all reported r2 values significant at the 0.05 probability level. Differences between regression coefficients or slopes of any two regressions (P < 0.05) were determined by testing the homogeneity of regression coefficients (Gomez and Gomez, 1984). All analyses were performed with SPSS Version 10.0 (SPSS, 1999).
 |
RESULTS AND DISCUSSION
|
|---|
Physical and chemical properties of the grassed and tilled Berks soils were generally similar across plots, while the no-till Watson soils had higher pH, clay, and organic C (Table 1). Mehlich-3 P and total P concentrations of all soils ranged widely, reflecting differences in manure application history; P concentrations in Berks soils were generally higher than in Watson soils.
View this table:
[in this window]
[in a new window]
|
Table 1. Mean and range (in parentheses) values of physical and chemical properties of the Berks and Watson soils.
|
|
Overland Flow Response
Response time between start of rainfall and initiation of overland flow was similar for the NPRP and WEPP simulators across all sites (Table 2). In addition, temporal trends in overland flow were consistent across soils for both simulators; it took appreciably longer for overland flow to occur and in lower volumes from the coarser-textured Berks soil than from the Watson soil. Differences in overland flow response time reflect the greater permeability of the Berks soil (15 to 150 mm h-1) than the Watson soil (1.5 to 15 mm h-1) (Eckenrode, 1985), as well as the influence of vegetative cover and tillage management on infiltration (Nash et al., 2002). For both simulators, overland flow occurred more rapidly after the start of rainfall on Day 2 than on Day 1 (Table 2), as soil moisture before rainfall was greater on Day 2 than Day 1. On the grassed Berks soil,
before rainfall increased from 28% on Day 1 to 39% on Day 2, while
before rainfall increased from 40 to 50% on the no-till Watson soil and from 23 to 35% on the tilled Berks soil. The standard error of
was consistently less than ±2% of mean values for each plot.
View this table:
[in this window]
[in a new window]
|
Table 2. Overland flow response to rainfall produced by the National Phosphorus Research Project (NPRP) (2-m2 plots) and Water Erosion Potential Predictor (WEPP) (32.6-m2 plots) simulators, both using an intensity of 75 mm h-1 for 30 min of overland flow. Data presented are averages for all plots on each soil and management type.
|
|
Total overland flow volume per unit area (L m-2) and peak flow rates per unit area (L m-2 min-1) were significantly greater under the NPRP simulator than under the WEPP simulator (Table 2). These differences suggest varying source areas of overland flow between the two plot scales. One possibility is that overland flow tends to be produced at the lower end of the plots, where subsurface lateral flow accumulates. For instance, given the low slope gradient, subsurface lateral flow may be intercepted by the lower plot boundary. As that water accumulates in the surface soil of the lower plot area, "infiltration" or "saturation excess" overland flow, as described by Nash et al. (2002), ensues. Under this process of overland flow generation, areas producing overland flow would be expected to expand upslope as rainfall continues.
A variety of observations support a variable source area hypothesis of overland flow generation. Both Berks and Watson soils in this study possess prominent plow pans that may impede vertical infiltration below the surface horizon and promote a component of lateral subsurface flow (Kleinman et al., 2003). Comparison of representative unit-width hydrographs (i.e., hydrographs normalized on a 1-m width basis) shows that flow from NPRP and WEPP plots increases at approximately the same rate (Fig. 4)
. In this case, flow from the NPRP plot plateaus with roughly 100% of rainfall converted to overland flow, while flow from the WEPP plot continues to rise for approximately 10 min, consistent with the expansion of the overland flowproducing zone beyond a 2-m flow path.

View larger version (15K):
[in this window]
[in a new window]
|
Fig. 4. Unit width hydrographs for National Phosphorus Research Project (NPRP) and Water Erosion Potential Predictor (WEPP) simulators for the first day of rainfall on the no-till Watson soil. Hydrographs are normalized on a 1-m-width basis for comparison. Bars represent standard errors of the mean values.
|
|
If overland flow occurred uniformly across the entire plot area, then different overland flow response times and hydrograph characteristics would be expected between the NPRP and WEPP simulators (Dingman, 1994). Assuming that overland flow is a function of saturated areas within the plots (i.e., saturation excess overland flow), unit area maximum flow rates indicate that 22 to 66% of the WEPP plot area was contributing overland flow, in comparison with 73 to 100% for the NPRP plots.
Spatial variability in surface soil infiltration properties may also contribute to the lower unit area overland flow from WEPP plots. Overland flow generated within the WEPP plots may thus have a greater opportunity to infiltrate in nonsaturated areas than in the fully saturated NPRP plots. Indeed, Srinivasan et al. (2002) observed that overland flow from a Berks soil infiltrated within several meters of its origin. The apparent plateau of flow from the WEPP plot at a rate representing less than 75% of incoming rainfall (Fig. 4) suggests that there is a greater potential for infiltration in the larger WEPP than NPRP plots.
Observed differences in overland flow volume between plot scales run contrary to the findings of Gascho et al. (1998), who found no significant difference between 6-m2 (3-m-long) and 622-m2 (43-m-long) plots under simulated rainfall. Soils studied by Gascho et al. (1998) were loamy sand, as opposed to silt loam and clay loam textures of the Berks and Watson soils, respectively. These loamy sands were also crust-forming soils, which may not promote near-surface lateral flow, a key process controlling variable source area hydrology. Results of that study may also reflect greater homogeneity of soil infiltration properties than in this study. For instance, the lower rainfall intensity (25 mm h-1) and longer rainfall periods (2 h) in Gascho et al. (1998) may have contributed to more uniform soil moisture distribution within plots.
Phosphorus Transport
Dissolved reactive P concentration of overland flow was significantly greater with the NPRP simulator than with the WEPP simulator (Table 3). For both simulators, DRP concentrations were lower on Day 2 than on Day 1, probably due to a temporary dilution of the pool of P that is released to overland flow in the interacting depth of surface soil. While the greater infiltration of rainfall within the WEPP plots may have resulted in greater vertical translocation of soluble P from the effective depth of soiloverland flow interaction, differences in DRP concentration between simulators cannot be attributed solely to dilution, as overland flow volume was greater with the NPRP than WEPP simulator (Table 2). Rather, significantly lower sediment discharge from the NPRP than WEPP plots (Table 3) probably resulted in a lower readsorption of P by particulates during overland flow (Sharpley et al., 1981; McDowell and Sharpley, 2002). Differences in sediment discharge between NPRP and WEPP plots are consistent with the shorter flow-path length of the NPRP (2 m) than WEPP (10.7 m) plots.
View this table:
[in this window]
[in a new window]
|
Table 3. Overland flow, P concentration, and sediment discharge from each soil and management type using the National Phosphorus Research Project (NPRP) (2-m flow-path) and Water Erosion Potential Predictor (WEPP) (10.7-m flow-path) rainfall simulator for an intensity of 75 mm h-1 for 30 min of overland flow.
|
|
The lower sediment discharge with the NPRP simulator is reflected in significantly lower PP concentrations in overland flow from the NRPR plots (Table 3). Again, it is likely that the longer flow-path length of the WEPP plots caused an increase in overland flow velocity, erosivity, and subsequent entrainment of particulates and associated P relative to the NPRP plots (Table 3). The selective erosion of fine particles also increased with flow-path length (Table 4). For instance, for the grassed Berks soil, sediment eroded from the WEPP plots was comprised of an average 26% clay-sized particles (<2 µm) during Day 1 overland flow, while sediment from the NPRP plots averaged only 19% clay (Table 4). This probably explains the wider PP to DRP ratio in Day 1 overland flow for the WEPP (4.7) than for the NPRP simulator (2.3) (Table 4). Similar enrichments of clay-sized particles and wider PP to DRP ratios were observed for the tilled Berks and no-till Watson soils and during Day 2 overland flow for all sites (Table 4). The difference in PP to DRP ratio results from the fact that P associated with clay-sized particles is less desorbable than P sorbed by sand-sized particles (McDowell et al., 2001).
View this table:
[in this window]
[in a new window]
|
Table 4. Percent clay-sized material (<2 µm) in eroded sediment and particulate phosphorus (PP) to dissolved reactive phosphorus (DRP) ratio averaged for overland flow from the Berks and Watson sites for the National Phosphorus Research Project (NPRP) and Water Erosion Potential Predictor (WEPP) simulators.
|
|
The effect of flow-path length on PP entrainment in overland flow and the erosional process was evaluated by calculating the phosphorus enrichment ratio (PER) as the P concentration of sediment discharged (mg PP kg-1 sediment) divided by that of source soil (mg TP kg soil-1). The PER was a function of sediment discharge for both NPRP and WEPP simulators (Fig. 5)
. However, PER at any given sediment discharge was significantly greater (P < 0.05) for WEPP than for NPRP plots for each of the soil treatments (grassed Berks and tilled Berks and no-till Watson). This results from the greater selective erosion of clay-sized particles in overland flow from WEPP plots. Although basic differences in WEPP (rotating boom and pulsed rainfall) and NPRP (stationary and continuous rainfall) simulators may have affected particle-size sorting, the data suggest this to be minimal. Even so, analysis of homogeneity of the PER and sediment discharge regression showed that slopes were statistically similar (P < 0.05) for NPRP and WEPP plots (Fig. 5). This suggests that although the selective removal of fines was greater at longer flow-path lengths, similar erosion processes governing PP transport in overland flow were operating at both plot scales.

View larger version (24K):
[in this window]
[in a new window]
|
Fig. 5. Relationship between P enrichment ratio and sediment discharge of overland flow generated by the National Phosphorus Research Project (NPRP) and Water Erosion Potential Predictor (WEPP) rainfall simulators for grassed and tilled Berks and no-till Watson soils.
|
|
At a much broader scale, Sharpley et al. (1985) calculated PER for up to 50 overland flow events from 2- to 6-ha watersheds in Oklahoma over two years. In that study, slopes of the Ln PERLn sediment discharge regression averaged 0.286 for six native grass and six conventionally tilled wheat (Triticum aestivum L.) watersheds. In this study, regression slopes were 0.281 and 0.301 for the NPRP and WEPP simulators, respectively.
Relationship between Overland Flow Phosphorus and Soil Phosphorus
The concentration of DRP in overland flow increased with Mehlich-3 P concentration of plot soils, although DRP concentrations were greater from the smaller NPRP than WEPP plots (Fig. 6)
. Analysis of homogeneity of regression coefficients for the relationship of DRP and Mehlich-3 P showed that absolute values of regression slopes for Berks (grassed and tilled) and Watson soils were statistically similar (P < 0.05) for NPRP and WEPP simulators (Fig. 6). For the grassed and tilled Berks soil, regression slopes using the NPRP simulator (0.0022 and 0.0105) were similar (P < 0.05) to those associated with the WEPP simulator (0.0021 and 0.0099). For the Watson soil, slopes were 0.0030 for NPRP and 0.0029 for WEPP simulators (Fig. 6).

View larger version (28K):
[in this window]
[in a new window]
|
Fig. 6. Relationship between the dissolved reactive P concentration of overland flow generated by the National Phosphorus Research Project (NPRP) and Water Erosion Potential Predictor (WEPP) rainfall simulators for grassed and tilled Berks and no-till Watson soils.
|
|
The similarity in DRPMehlich-3 P regression slopes between NPRP and WEPP simulators suggests that overland flow processes controlling soil P release and transport are independent of simulator type, flow-path length, and plot size. However, differences in DRP concentrations and overland flow volumes between simulators, as reflected by the y intercepts of the regression equations (Fig. 6), limits their use in quantifying P loss from agricultural landscapes due to the obvious effects of scale on P loss. This is consistent with NPRP objectives, which are not to quantify P losses as a function of field-scale agricultural management, but to determine factors controlling the relationship between overland flow P and soil P and then define this relationship for a wide range of soils (Sharpley et al., 2002b). Rainfall simulators and small plots cannot reproduce flow process occurring over a landscape or hillslope and as such must be limited to elucidating flowsoil P interdependencies.
Even so, some critical processes controlling overland flowsoil P interdependency appear to be analogous for overland flow generated from rainfall simulators and larger-scale watersheds (Sharpley et al., 2002a). For example, the slopes of the overland flowMehlich-3 soil P regression were 0.0027 for ten 2- to 6-ha grassed watersheds and 0.0135 for eight 2- to 3-ha conventionally tilled wheat watersheds in Oklahoma (Sharpley et al., 1991; Smith et al., 1991). In a recent review, Sharpley et al. (2002a) found that regression slopes were lower for 38 grasses (0.00150.0035, mean of 0.0022 and standard error of 0.0007) than for no-till (0.00230.0044, mean of 0.0036 and standard error of 0.0007) or cultivated watersheds (0.00900.0152, mean 0.0112 and standard error of 0.0007). These regression slopes and the relative differences between tillage practices are similar to values obtained in this study with NPRP and WEPP simulators for grassed (0.0022 and 0.0021), no-till (0.0030 and 0.0029), and tilled sites (0.0105 and 0.0099) (Fig. 6). The consistency of regression slopes derived from small-plot (<40 m2) to watershed (>2 ha) scales is evidence that similar soil P release and transport processes operate across these scales.
 |
CONCLUSIONS
|
|---|
Plot length influences hydrology, sediment discharge, and concentration of DRP, PP, and TP in overland flow. Differences in overland flow volume, peak flow, and hydrograph characteristics indicate that a proportionally smaller area of the WEPP plots contributed overland flow, and that infiltration for the WEPP plots was greater under equilibrium conditions than within the NPRP plots. Even so, trends in overland flow over time and between soils were generally consistent for the two simulators.
From the standpoint of water quality, DRP concentrations were consistently greater in overland flow from 2- than 10.7-m-long plots. In contrast, sediment discharge and PP concentrations were lower in overland flow from 2- than 10.7-m-long plots. This suggests that the greater enrichment of clay-sized particles in overland flow from longer plots controls not only the form but the concentration of P in overland flow.
No comparison was made of rainfall simulator and field-scale estimates of overland flow characteristics, P concentration, and sediment discharge, due to a lack of data for natural rainfall and broader scales (e.g., fields, hillslopes) that covered the range of soil P levels and tillage treatments used with the simulators. However, the similarity of overland flow Psoil P relationships derived from plot and watershed-based studies, supports the conclusion that processes controlling P release from soil and transport as DRP are similar for overland flow generated from rainfall simulators and natural rainfall. Comparable PER values for simulator and natural overland flow also indicate that analogous processes control PP detachment and entrainment. Furthermore, the same DRP and PP transport processes appear to be operating across plot scales; NPRP and WEPP simulators used 2- and 32.6-m2 plots while cited watershed studies ranged from 2 to 6 ha.
Because of the effects of flow-path length and landscape hydrology on overland flow and P transport, this study shows that inferences from simulated rainfall-runoff experiments on small plots should be limited to quantifying relationships between overland flow DRP and soil P and in assessing sediment discharged, PP, and P enrichment. This is the case for grassed and no-till situations where plant and residue cover protects the surface soil from flow interaction, as well as for tilled situations.
 |
ACKNOWLEDGMENTS
|
|---|
The contributions of Earl Jacoby (deceased, formerly at USDA-ARS, Klingerstown, PA), Todd Strohecker (USDA-ARS, Klingerstown, PA), and Terry Troutman (USDA-ARS, Klingerstown, PA) to rainfall simulations, M.S. Srinivasan to hydrologic evaluations, and Lou Saporito for statistical support and analysis are gratefully appreciated.
 |
NOTES
|
|---|
Mention of trade names does not imply endorsement by the USDA.
 |
REFERENCES
|
|---|
- Bloschl, G., and M. Sivapalan. 1995. Scale issues in hydrological modeling: A review. Hydrol. Processes 9:251290.
- Brakensiek, D.L., H.B. Osborn, and W.J. Rawls. 1979. Runoff. p. 75214. In Field manual for research in agricultural hydrology. USDA Agric. Handb. 224. U.S. Gov. Print. Office, Washington, DC.
- Bremner, J.M. 1996. NitrogenTotal. p. 10851121. In D.L. Sparks (ed.) Methods of soil analysis. Part 3. SSSA Book Ser. 5. SSSA, Madison, WI.
- Carpenter, S.R., N.F. Caraco, D.L. Correll, R.W. Howarth, A.N. Sharpley, and V.H. Smith. 1998. Nonpoint pollution of surface waters with phosphorus and nitrogen. Ecol. Appl. 8:559568.
- Day, P.R. 1965. Particle fractionation and particle size analysis. p. 545567. In C.A. Black, D.D. Evans, J.L. White, L.E. Ensminger, and F.E. Clark (ed.) Methods of soil analysis. Part 1. Agron. Monogr. 9. ASA, Madison, WI.
- Dingman, S.L. 1994. Physical hydrology. Prentice Hall Int., Englewood Cliffs, NJ.
- Eckenrode, J.J. 1985. Soil survey of Northumberland County, Pennsylvania. USDA Soil Conserv. Serv., U.S. Gov. Print. Office, Washington, DC.
- Gascho, G.J., R.D. Wauchope, J.G. Davis, C.C. Trueman, C.C. Dowler, J.E. Hook, H.R. Sumner, and A.W. Johnson. 1998. Nitrate-nitrogen, soluble, and bioavailable phosphorus runoff from simulated rainfall after fertilizer application. Soil Sci. Soc. Am. J. 62:17111718.[Abstract/Free Full Text]
- Gomez, K.A., and A.A Gomez. 1984. Statistical procedures for agricultural research. 2nd ed. John Wiley & Sons, New York.
- Humphry, J.B., T.C. Daniel, D.R. Edwards, and A.N. Sharpley. 2002. A portable rainfall simulator for plot-scale runoff studies. Appl. Eng. Agric. 18:199204.[ISI]
- Kleinman, P.J.A., B.A. Needelman, A.N. Sharpley, and R.W. McDowell. 2003. Using soil profile data to assess phosphorus leaching potential in manured soils. Soil Sci. Soc. Am. J. 67:215224.[Abstract/Free Full Text]
- Le Bissonnais, Y., H. Benkhadra, V. Chaplot, D. Fox, D. King, and J. Daroussin. 1998. Crusting, runoff and sheet erosion in silty loam soils at various scales and upscaling from m2 to small catchments. Soil Tillage Res. 46:6980.
- McDowell, R.W., and A.N. Sharpley. 2002. Effect of plot scale and upslope phosphorus source on phosphorus loss in overland flow. Soil Use Manage. 18:112119.
- McDowell, R.W., A.N. Sharpley, and G. Folmar. 2001. Phosphorus export from an agricultural watershed: Linking source and transport mechanisms. J. Environ. Qual. 30:15871595.[Abstract/Free Full Text]
- Mehlich, A. 1984. Mehlich 3 soil test extractant: A modification of Mehlich 2 extractant. Commun. Soil Sci. Plant Anal. 15:14091416.
- Miller, W.P. 1987. A solenoid-operated, variable intensity rainfall simulator. Soil Sci. Soc. Am. J. 51:832834.[Abstract/Free Full Text]
- Murphy, J., and J.P. Riley. 1962. A modified single solution method for the determination of phosphate in natural waters. Anal. Chim. Acta 27:3136.
- Nash, D., D. Halliwell, and J. Cox. 2002. Hydrological mobilization of pollutants at the field/slope scale. p. 225242. In P.M. Haygarth and S.C. Jarvis (ed.) Agriculture, hydrology and water quality. CABI Publ., Wallingford, UK.
- Pionke, H.B., W.J. Gburek, and A.N. Sharpley. 2000. Critical source area controls on water quality in an agricultural watershed located in the Chesapeake Basin. Ecol. Eng. 14:325335.
- Sharpley, A.N. (ed.) 2000. Agriculture and phosphorus management: The Chesapeake Bay. CRC Press, Boca Raton, FL.
- Sharpley, A.N., P.J.A. Kleinman, R.W. McDowell, M. Gitau, and R.B. Bryant. 2002a. Modeling phosphorus transport in agricultural watersheds: Processes and possibilities. J. Soil Water Conserv. 57:425439.
- Sharpley, A.N., P.J.A. Kleinman, R.J. Wright, T.C. Daniel, B. Joern, R. Parry, and T. Sobecki. 2002b. The National Phosphorus Project: Addressing the interface of agriculture and environmental phosphorus management in the USA. Section 2.3. Indicators for environmental performance. p. 95100. In J. Steenvoorden, F. Claessen, and J. Willems (ed.) Agricultural effects on ground and surface waters: Research at the edge of science and society. Int. Assoc. of Hydrol. Sci. Publ. 273. IAHS Press, Wallingford, UK.
- Sharpley, A.N., R.G. Menzel, S.J. Smith, E.D. Rhoades, and A.E. Olness. 1981. The sorption of soluble P by soil material during transport in runoff from cropped and grassed watersheds. J. Environ. Qual. 10:211215.[Abstract/Free Full Text]
- Sharpley, A.N., S.J. Smith, W.A. Berg, and J.R. Williams. 1985. Nutrient runoff losses as predicted by annual and monthly soil sampling. J. Environ. Qual. 14:354360.[Abstract/Free Full Text]
- Sharpley, A.N., S.J. Smith, J.R. Williams, O.R. Jones, and G.A. Coleman. 1991. Water quality impacts associated with sorghum culture in the Southern Plains. J. Environ. Qual. 20:239244.[Abstract/Free Full Text]
- Sharpley, A.N., J.L. Weld, D.B. Beegle, P.J.A. Kleinman, W.J. Gburek, P.A. Moore, Jr., and G. Mullins. 2003. Development of phosphorus indices for nutrient management planning strategies in the U.S. J. Soil Water Conserv. 58:137152.
- Shelton, C.H., R.D. von Bernuth, and S.P. Rajbhandari. 1985. A continuous-application rainfall simulator. Trans. ASAE 28:11151119.
- Simanton, J.R., and K.G. Renard. 1992. Upland erosion research on rangeland. p. 335375. In A.J. Parsons and A.D. Abrahams (ed.) Overland flowHydraulics and erosion mechanics. John Wiley & Sons, New York.
- Smith, S.J., A.N. Sharpley, J.W. Naney, W.A. Berg, and O.R. Jones. 1991. Water quality impacts associated with wheat culture in the Southern Plains. J. Environ. Qual. 20:244249.[Abstract/Free Full Text]
- SPSS. 1999. SPSS Version 10.0. SPSS, Chicago, IL.
- Srinivasan, M.S., W.J. Gburek, and J.M. Hamlett. 2002. Dynamics of stormflow generationA hillslope-scale field study in east-central Pennsylvania, USA. Hydrol. Processes 16:649665.
- United States Geological Survey. 1999. The quality of our nation's waters: Nutrients and pesticides. USGS Circ. 1225. USGS Info. Serv., Denver.
- USDA and USEPA. 1999. Unified national strategy for animal feeding operations [Online]. Available at http://cfpub.epa.gov/npdes/afo/ustrategy.cfm (verified 12 Aug. 2003). USEPA, Washington, DC.
- USEPA. 1996. Environmental indicators of water quality in the United States. EPA 841-R-96-002. USEPA Office of Water, U.S. Gov. Print. Office, Washington, DC.
- Wauchope, R.D., and B. Burgoa. 1995. Pesticide runoff studies: Towards a new protocol. p. 273282. In M.L. Leng, E.M.K. Leovey, and P.L. Zubkoff (ed.) Agrochemical environmental fate: State of the art. Lewis Publ., Boca Raton, FL.
Related articles in JEQ:
- This Issue in Journal of Environmental Quality
JEQ 2003 32: 1931-1938.
[Full Text]
This article has been cited by other articles:

|
 |

|
 |
 
D. J. Wagner, H. A. Elliott, R. C. Brandt, and D. Jaiswal
Managing Biosolids Runoff Phosphorus Using Buffer Strips Enhanced with Drinking Water Treatment Residuals
J. Environ. Qual.,
June 23, 2008;
37(4):
1567 - 1574.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
M. S. Srinivasan, P. J. A. Kleinman, A. N. Sharpley, T. Buob, and W. J. Gburek
Hydrology of Small Field Plots Used to Study Phosphorus Runoff under Simulated Rainfall
J. Environ. Qual.,
October 24, 2007;
36(6):
1833 - 1842.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
T. K. Udeigwe, J. J. Wang, and H. Zhang
Predicting Runoff of Suspended Solids and Particulate Phosphorus for Selected Louisiana Soils Using Simple Soil Tests
J. Environ. Qual.,
July 17, 2007;
36(5):
1310 - 1317.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
B. A. Ebel, K. Loague, W. E. Dietrich, D. R. Montgomery, R. Torres, S. P. Anderson, and T. W. Giambelluca
Near-surface hydrologic response for a steep, unchanneled catchment near Coos Bay, Oregon: 1. sprinkling experiments
Am J Sci,
April 1, 2007;
307(4):
678 - 708.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
D. M. Butler, D. H. Franklin, N. N. Ranells, M. H. Poore, and J. T. Green Jr.
Ground Cover Impacts on Sediment and Phosphorus Export from Manured Riparian Pasture
J. Environ. Qual.,
October 27, 2006;
35(6):
2178 - 2185.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
A. R. Guidry, F. V. Schindler, D. R. German, R. H. Gelderman, and J. R. Gerwing
Using Simulated Rainfall to Evaluate Field and Indoor Surface Runoff Phosphorus Relationships
J. Environ. Qual.,
October 27, 2006;
35(6):
2236 - 2243.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
C. J. Penn, G. L. Mullins, L. W. Zelazny, and A. N. Sharpley
Estimating Dissolved Phosphorus Concentrations in Runoff from Three Physiographic Regions of Virginia
Soil Sci. Soc. Am. J.,
September 20, 2006;
70(6):
1967 - 1974.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
H. Zhang, J. L. Schroder, R. L. Davis, J. J. Wang, M. E. Payton, W. E. Thomason, Y. Tang, and W. R. Raun
Phosphorus Loss in Runoff from Long-term Continuous Wheat Fertility Trials
Soil Sci. Soc. Am. J.,
December 2, 2005;
70(1):
163 - 171.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
P. J. A. Kleinman, A. N. Sharpley, T. L. Veith, R. O. Maguire, and P. A. Vadas
Evaluation of Phosphorus Transport in Surface Runoff from Packed Soil Boxes
J. Environ. Qual.,
July 1, 2004;
33(4):
1413 - 1423.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
C. J. Penn, G. L. Mullins, L. W. Zelazny, J. G. Warren, and J. M. McGrath
Surface Runoff Losses of Phosphorus from Virginia Soils Amended with Turkey Manure Using Phytase and High Available Phosphorus Corn Diets
J. Environ. Qual.,
July 1, 2004;
33(4):
1431 - 1439.
[Abstract]
[Full Text]
[PDF]
|
 |
|