SSSAJ Journal of Natural Resources and Life Sciences Education
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (6)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Needelman, B. A.
Right arrow Articles by Kleinman, P. J. A.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Needelman, B. A.
Right arrow Articles by Kleinman, P. J. A.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Needelman, B. A.
Right arrow Articles by Kleinman, P. J. A.
Related Collections
Right arrow Watershed and Landscape Processes
Right arrow Runoff
Right arrow Surface Hydrology
Right arrow Nutrient Management
Published in Soil Sci. Soc. Am. J. 68:914-923 (2004).
© 2004 Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA

DIVISION S-6—SOIL & WATER MANAGEMENT & CONSERVATION

Surface Runoff along Two Agricultural Hillslopes with Contrasting Soils

Brian A. Needelman*,a, William J. Gburekb, Gary W. Petersenc, Andrew N. Sharpleyb and Peter J. A. Kleinmanb

a Dep. of Natural Resource Sci. and Landscape Architecture, Univ. of Maryland, 1112 H.J. Patterson Hall, College Park, MD 20742
b Pasture Systems and Watershed Management Research Unit, USDA-ARS, Curtin Rd., University Park, PA 16802
c Dep. of Crop and Soil Sci., The Pennsylvania State Univ., 116 A.S.I. Building, University Park, PA 16802

* Corresponding author (bneed{at}umd.edu).


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
The targeting of critical surface runoff-producing zones should account for the influence of subsurface soil characteristics. In this study we assessed the runoff response of contrasting colluvial and residual soils. The study was conducted along two hillslopes within a 39.5-ha mixed land use watershed in Pennsylvania. Six sites (four colluvial, two residual) were monitored for runoff, hydraulic head, water table depth, and soil water content. A total of 111 rainfall events were monitored during the periods of July to December 2000, April to December 2001, and April to December 2002. Two high-intensity (5-min peak > 8 cm h–1) events had return periods of 2.5 and 4 yr. The colluvial soils are somewhat poorly and moderately well drained with fragipans and high clay content (37–44%) argillic horizons (fine, mixed, semiactive, mesic Aquic Fragiudalfs); the residual soils are well drained with moderate clay content (24%) argillic horizons (fine-loamy, mixed, semiactive, mesic Typic Hapludults). Across all events, overall runoff yields averaged 2.4% from the four colluvial sites and 0.01% from the two residual sites. The two colluvial sites with the greatest runoff production were located at the base of a primarily colluvial hillslope. The largest events at these sites occurred during periods of surface saturation (soil surface to a depth of at least 30 cm). These results suggest that nonwinter P management for these residual soils should focus on rare, large events. Nutrient management planning could be improved if runoff estimation methods were to better integrate information on subsurface and upslope soil hydrologic properties.

Abbreviations: VSA, variable source area


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
PHOSPHORUS IS TRANSPORTED primarily through surface runoff from agricultural lands and is a major cause of eutrophication in surface waters (Carpenter et al., 1998). In the long term, regional nutrient balances may be necessary for environmental and agricultural sustainability (Magdoff et al., 1997). In the short term, however, farmers require an economical manure utilization strategy. Landscape areas with a high probability of P loss via surface runoff, termed critical source areas, must be identified and targeted with specific P management strategies. This study was developed in response to the need to identify surface runoff source areas within agricultural watersheds to improve the targeting of pollutant loss prevention strategies.

Variable source area (VSA) hydrology has been proposed as an effective framework to describe the response of humid temperate watersheds to precipitation (Ward, 1984). Under VSA hydrology, a small percentage of the landscape is the source of stormflow; source areas contract and expand within and between events. Variable source areas are a function of topography, soils, geology, climate, and management. Within VSA hydrology, surface runoff generation is classified as either infiltration excess or saturation excess (Sklash, 1990). Infiltration excess, or Hortonian overland flow, occurs when rainfall intensities exceed the infiltration capacity of a soil. Variable infiltration capacities in the landscape cause partial areas of infiltration-excess surface runoff (Betson, 1964). Saturation excess occurs when the water table rises to saturate the soil profile, filling storage zones and minimizing infiltration capacity. Ideally, P management strategies should differ depending on the dominant surface runoff generation mechanism in a watershed. An infiltration-excess-based approach should target soils with a low infiltration capacity, while a saturation-excess-based strategy should target near-stream and other zones that are subject to surface saturation (Gburek et al., 1996) regardless of infiltration capacity at these sites.

The saturation-excess mechanism has been studied primarily under forest and grassland vegetation (Bonell, 1993). In most temperate forested landscapes, infiltration excess is a minor mechanism of runoff generation; rather, stormflow is dominated primarily by shallow subsurface flow, with contributions from saturation-excess surface runoff (Anderson and Burt, 1990). The lack of infiltration-excess surface runoff in forested landscapes is due to the high infiltration capacity of forest surface soil horizons, including O horizons. Few field studies have been conducted in agroecosystems to investigate the relative occurrence of infiltration-excess and saturation-excess runoff (Betson and Marius, 1969). Typically, agricultural soils, particularly seasonally bare soils, exhibit lower infiltration capacities and less macroporosity than do forest soils. These factors favor the infiltration-excess mechanism over the saturation-excess mechanism (Burch et al., 1987).

This field study was located within a subwatershed of the Mahantango Creek watershed. The Mahantango Creek watershed has been reported as a VSA watershed (Pionke et al., 1996). Variable source areas are thought to occur primarily in the near-stream zones in response to the close proximity of the water table to the land surface, which causes seep zones and high antecedent soil water content contents (Gburek and Sharpley, 1998). Stormflow from these VSAs is thought to be primarily saturation-excess surface runoff and rapid subsurface flow (Gburek et al., 1996). The VSA-based simulation modeling of a 26-ha mixed land use watershed located near the field site for the present study predicted that 98% of the runoff volume is produced from 14% of the area (Zollweg et al., 1995). These studies in the Mahantango Creek watershed have been conducted only under grassland vegetation. In cropped soils, cultivation, wheel traffic, and other management practices may decrease infiltration capacities so that infiltration-excess runoff may be of greater relative importance in runoff generation.

Subsurface hydrologic properties may be critical to the occurrence of surface saturation because of the central role of shallow water table dynamics. Subsurface properties that may influence the occurrence of surface saturation include the presence of a fragipan or other soil layers that may cause episaturation, the hydraulic conductivity of subsurface horizons, depth to bedrock, depth of topsoil, and the depth to a water table. The occurrence of the water table rise that causes surface saturation may be correlated to soil hydromorphology and natural drainage class (Engman and Rogowski, 1974). Runoff and drainage class have been indirectly linked through soil water content measurements in a study of similar soils as those in the present study (Henninger et al., 1976). A good correlation was observed between the extent of somewhat poorly drained soils and the maximum seasonal extent of the surface-saturated zone in a small catchment in northeastern Vermont (Dunne et al., 1975; Moore et al., 1976).

The objectives of this study were to assess the surface runoff response to rainfall at four colluvial and two residual soil sites and investigate the factors influencing runoff generation at these sites including rainfall event characteristics, surface saturation, management, and antecedent soil water content. For additional information on this study see Needelman (2002).


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Study Location
The study was conducted along two hillslopes in watershed FD-36 (Fig. 1) , a 39.5-ha subwatershed of the Mahantango Creek, which is a tributary to the Susquehanna River and ultimately the Chesapeake Bay (see Gburek and Sharpley, 1998, for additional details on this watershed). Watershed FD-36 is located within a valley portion of the folded and faulted Ridge and Valley Physiographic Province. The watershed has mixed land use: 50% soybean [Glycine max (L.) Merr.], wheat (Triticum aestivum L.), or corn (Zea mays L.), 20% pasture, and 30% woodland (Gburek and Sharpley, 1998). Information on management was obtained through an annual farmer survey with information tabulated on a field-by-field basis.



View larger version (56K):
[in this window]
[in a new window]
 
Fig. 1. Land use and flume locations at 39.5-ha watershed FD-36.

 
Watershed FD-36 lies within watershed WE-38 (7.3 km2), a primary research site of the USDA-ARS for more than 30 yr. Rainfall, streamflow, and meteorological data have been collected for WE-38 since 1967, and ground water data since 1972. Climate is temperate and humid, average rainfall is approximately 103 cm yr–1, and streamflow is about 45 cm yr–1 (Gburek and Sharpley, 1998). Research within the larger watershed WE-38 has shown that both the surface and subsurface flow systems are predominantly self contained at the scale of watershed FD-36 (Gburek and Folmar, 1999).

Instrumentation
Two hillslopes were instrumented and continuously monitored during nonwinter periods to compare the hydrologic response of colluvial and residual soils. The study could not be conducted during the winters because many of the sensors are susceptible to frost damage. Three sites were established on each hillslope at distances of about 15, 35, and 65 m from the stream channel (Fig. 2 and Table 1). Sites on the south side of the channel are labeled S1, S2, and S3, with S1 closest to the stream, and on the north side were labeled N1, N2, and N3, with N1 closest to the stream.



View larger version (36K):
[in this window]
[in a new window]
 
Fig. 2. Location of surface runoff flumes, topography, and soils (parent material and drainage class) at FD-36 watershed. Contour interval is 0.5 m.

 

View this table:
[in this window]
[in a new window]
 
Table 1. Selected site characteristics at monitoring sites.

 
Management at the sites was performed by local farmers following standard practices on surrounding fields. We encouraged the farmers to come as close as possible to the sensors with their equipment. Areas that were not accessible with the farmer's equipment were managed by researchers, including small-engine tillage, herbicide application, hand seeding, hand tillage, and weeding. Areas immediately around the sensors were kept bare to avoid disturbance.

Sites were monitored during the periods of July through December 2000, April through December 2001, and April through December 2002. At each site, a surface H-flume was installed to monitor surface runoff (Fig. 3) . A berm was constructed to give the flume a runoff collection width of 10 m. This width was chosen to be sufficiently large to aggregate small-scale runoff processes such as rill flow, yet small enough to allow for intense data collection in the upslope source area.



View larger version (17K):
[in this window]
[in a new window]
 
Fig. 3. Layout of theta probe, piezometers, flume, and paired wells and surface runoff indicators. Wells are on a 10- x 10-m grid.

 
Two piezometers and a soil water content probe were installed 2 m upslope from each flume (Fig. 3). The piezometers were installed at 30- and 60-cm depths and were monitored with a pressure transducer (Water Log Series, Model H-310, H2O FX, Logan, UT). Piezometer installation holes were creating by removing 6.25-cm-diam. soil cores, which were retained for description. Piezometers were constructed with 6.25-cm-diam. PVC tubes horizontally slotted across a depth of 5 cm, and were backfilled with soil and sealed with bentonite to the surface. Hydraulic head was referenced to the soil surface rather than to a common elevation datum because the data were interpreted at the measurement point rather than across the landscape. Surface soil water content was monitored with a theta probe (Type mL2x, Delta-T Devices, Cambridge, UK) installed next to the piezometers to a depth of 5 cm (Fig. 3). Dielectric voltage readings were converted to volumetric water content values by a soil-specific calibration. Note that while a single soil water content sensor is not sufficient to precisely estimate the soil water content of a site, it does provide a comparative measure to assess trends across time and between sites.

To prevent intersite effects, sites were staggered along the hillslopes. This was not possible between Sites S1 and S2 because of local topography (Fig. 2). A network of roughly 7- x 7-cm ditches were excavated to redistribute runoff water passing through the flume at S2 into the area upslope of flume at S1.

Eight paired 45-cm wells and surface runoff indicators were installed upslope from each flume (Fig. 3). Wells were constructed and installed in the same manner as were the piezometers, except tube diameters were 3.75 cm and slotting extended from 5 to 60 cm (rather than 5 cm of slotting at 30 or 60 cm). The first two well runoff indicator pairs were installed 2 m on both sides and 1 m upslope of the piezometers. The other six pairs were installed on a 10- x 10-m grid upslope from the flume. Wells were monitored with subsurface saturation sensors following Srinivasan et al. (2000). These sensors indicate whether the water table is at or above 1, 5, 15, 30, or 45 cm from the soil surface. Surface runoff indicators (Srinivasan et al., 2000) are miniature V-notch weirs (40-cm width) with an electrical conductivity sensor indicating the presence or absence of surface water. All data were collected on a 5-min interval with a Campbell CR-10 datalogger to capture rapid water table changes during events (Calmon and Day, 1999).

Data Analysis
This study was observational rather than inferential. It was not feasible to monitor a sufficient number of sites with the intensity of data collection performed at each site to address objectives with statistical hypothesis testing. Data were analyzed on a site-by-site basis across time, comparing runoff response under varying precipitation and site conditions. Conclusions on the comparative runoff generation between sites were based on broad, clearly discernible differences.

A topographically based surface-runoff-contributing area was calculated for each flume (Table 1). All flumes except S1 were located in a hillslope position with a plan curvature near zero. Therefore, contributing areas were calculated by extending the berm width to a road, which forms a surficial topographic boundary on each hillslope. Site S1 is located within a convergent position. The S1-contributing area was delineated by the upslope road and a topographic bench in the landscape.

Note that the topographic contributing area conflicts with the concept of a VSA. Nonetheless, the topographic contributing area is necessary for the calculation of runoff yields, an important comparative measure between sites. Runoff yields were calculated as the volume of runoff divided by the volume of rainfall falling within the topographically defined contributing area.

Soils of Study Site
An intensive soil survey was conducted along the two hillslopes through a cooperative agreement with the Pennsylvania Natural Resources Conservation Service (NRCS) (Needelman, 2002). Survey work included pit and soil core descriptions and laboratory characterization. The soils are formed in shale, siltstone, and sandstone residuum and colluvium. Two contrasting soil groups are found along the hillslopes, differing primarily in subsurface hydrologic characteristics including the presence of a fragipan, the clay content of the argillic horizons, and natural drainage class (Table 2). Ap horizons are generally 25 to 30 cm thick and are underlain directly by Bt horizons. The upper portions of both hillslopes have well-drained residual soils classified as fine-loamy, mixed, semiactive, mesic Typic Hapludults (Leck Kill series). The bottom of the northern hillslope has moderately well-drained colluvial soils classified as fine-loamy, mixed, active, mesic Oxyaquic Fragiudalfs (Hustontown series). The southern hillslope has somewhat poorly drained colluvial soils classified as fine-loamy, mixed, semiactive, mesic Aquic Fragiudalfs (Albrights series). On the southern hillslope, the colluvial soils extend about 90 m from the stream channel. There is a 15-m transition zone where the fragipan grades out and redoximorphic features decrease. Sites S1, S2, and S3 are colluvial soil sites (Fig. 2). On the north side of the channel, the bottom 15 m of the hillslope are colluvial soils, followed by a 15-m transition zone; the remaining hillslope is residual soils. Site N1 is a colluvial soil site while N2 and N3 are residual soil sites (Fig. 2).


View this table:
[in this window]
[in a new window]
 
Table 2. Selected soil characteristics at monitoring sites.

 
We have chosen to use parent material (colluvial vs. residual) to group the soils in this study because parent material is the dominant soil-forming factor influencing soil variability in this watershed. The colluvial soils are relatively young (periglacial). There is little evidence that the climate, organisms, or time factors vary within colluvial or residual units within the study area (e.g., there are no observable aspect effects). Topography does influence drainage class and related hydrologic parameters, but is included within the study because individual sites are defined by their landscape position. Colluviation most likely occurred during the most recent glaciation, about 18000 yr ago (Ciolkosz et al., 1989). Colluvial materials are thinly stratified, alternating between redder materials derived from red fine-grained sandstones and gray materials derived from gray mudstones. These materials are commonly intermixed. There are strongly contrasting rock fragment, sand, and clay percentages between the stratified colluvial layers. Fragipans tend to form within the loamy, channery materials while the argillic horizons are formed within silty, clayey materials. The high clay contents observed in the colluvial argillic horizons are only partly due to illuviation. Rather, these argillic horizons are formed in discrete colluvial layers (differentiated by rock fragment characteristics and matrix colors) with fine textures. The sedimentary rock fragments in these layers are fine-grained mudstones, which weather directly to clay-sized particles. These fine-textured layers are not present in the surficial residuum.

Detailed Site Characteristics
Cropping for the sites is presented in Table 1. Sites S2 and S3 were moldboard plowed and disc harrowed each spring before planting. The N1 site was chisel and disc tilled in the spring of 2000 and 2001 and was not tilled in 2002. The N2 site was disc harrowed in the fall of 1999. The N3 site was chisel and disc tilled in the spring of 2000 and the fall of 2001.

Site S1
Site S1 was converted to mowed grass in 1996 when a former runoff generation study was initiated (Fig. 2). In the former study, a network of saturation detectors were monitored throughout the grassed area. Results of the previous study remain unpublished. The conclusion of the previous study was that saturation-excess runoff occurred throughout the sampled area during major events, indicating that the extent of area monitored was not large enough to delineate the VSA. Fragipans at site S1 are typically shallow with minimum depths of 44 cm in some locations.

Site S2
Site S2 is in a toeslope position. Fragipan depth is generally deep in this area (80–150 cm), though depth is difficult to estimate because of broken argillic-fragipan boundaries and the presence of clay lenses within the fragipans in this area.

Site S3
Site S3 is located in the footslope position of a slope portion that grades from the shoulder with skeletal residual soils 40 m above the flume (upper extent of wells and runoff indicators) to lower rock fragment colluvial soils near the flume. While installing instrumentation, a French drain was discovered 2 m upslope from the flume running parallel to the berm. The bank of sensors located above the flume is just downslope of the French drain. The infrequent occurrence of a high water table during the study, despite the presence of redoximorphic features directly below the Ap horizon, may be because of this drain. This site therefore represents the artificial condition of a drained soil that is naturally somewhat poorly drained. Analysis of data from this site cannot be extrapolated to similar sites without French drains; nonetheless, this site proved worthwhile to elucidate the role of surface saturation in runoff generation. French drains are ubiquitous in this area and have not been mapped; it is possible that there are other drains in the study area. However, the sites with off-drained soils (S1, S2, and N1) had sustained high water tables during the study period, an indication that the effect of artificial drainage is minimal if present.

Site N1
Site N1 is located within a colluvial draw at the base of a primarily residual hillslope. Colluvium extends about 25 m from the stream at this location. A minor seep is located about 15 m directly upslope of the flume. Although this seep did not generate substantial overland flow during the study period, there is likely a channelized natural or artificial subsurface flow system that feeds the ground water in the site area.

Sites N2 and N3
Sites N2 and N3 have similar residual soils and are topographically similar: Both are located along an upland backslope that extends nearly to the watershed divide, though a road forms the upper boundary of the surface runoff contributing area. Site N2 was removed from the study following the 2001 monitoring period because of farmer participation and resource constraint limitations.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Precipitation Summary
A total of 111 rainfall events were monitored during the study period (July–December 2000, April–December 2001, and April–December 2002). The 2000 and 2001 periods were generally dry (24 and 31% lower than average precipitation) (30-yr average from Waltman et al., 1997). The 2002 period had a wet spring (15% greater than average), a dry summer (45% less than average), and a very wet fall (62% greater than average). Two brief events during the study period exceeded a 1-yr intensity-duration return period; both events had high peak intensities (>8.0 cm h–1). The event on 12 June 2001 was a 4-yr 20-min storm, and the event on 25 July 2001 was a 2.5-yr 30-min storm (Aron et al., 1986). The greatest rainfall depths recorded during the study period were 4.9 and 5.1 cm, both occurring during long (>24 h) events.

Surface Runoff
Aggregated runoff results at the six monitoring sites are given in Table 3 based on all 111 rainfall events, except that Site N2 was not in operation in 2002. Total runoff volumes were greatest at Sites S1 and S2 and minimal at the residual soil sites. Rainfall characteristics and runoff volumes for larger events (>2.0-cm depth) and selected smaller events (≤2.0 cm) are given in Table 4. The colluvial soils, Sites S1, S2, S3, and N1, produced dramatically greater runoff volumes than did the residual soils, Sites N2 and N3, for all major runoff events during the study period. For larger rainfall events, runoff volumes were generally greatest at Sites S1 and S2. The two sites located on residual soils, N2 and N3, produced little or no runoff during any event. Three events in the larger event group (1 Sept. 2000, 1 June 2001, 25 Nov. 2001, and 15 Sept. 2002) produced little or no runoff at all sites. There were several smaller events that produced runoff volumes equal to or greater than some of the large events. For example, at Site S2 the 22 June 2001 (2) event generated the largest volume of any smaller event, yet only had a rainfall depth of 0.94 cm.


View this table:
[in this window]
[in a new window]
 
Table 3. Hydrologic results for monitoring sites summarized from all events (111) during study period.

 

View this table:
[in this window]
[in a new window]
 
Table 4. Rainfall characteristics and runoff volumes for larger events (>2.0-cm depth) and selected smaller events (≤2.0 cm). Surface saturation status is shown in italic font (surface saturated) or normal font (non-surface saturated).

 
Overall runoff yields averaged 2.4% from the four colluvial sites and 0.01% from the two residual soil sites (Table 3); runoff yields were greater at colluvial sites for all runoff-generating events (Table 4). Runoff yields were greatest at Sites S1 and S2, with three larger and one smaller storms yielding >10% at both sites. At the residual soil sites, runoff yields were <2% for all events.

All sites except S1 had similar topographic contributing areas; the topographic contributing area for Site S1 is five to seven times larger than the area of the other sites. However, the actual contributing areas, as indicated by the runoff indicators, only included a small percentage of the full topographic contributing area. This follows the basic precepts of VSA hydrology, that contributing areas are dynamic during events and may not be directly related to the topographic contributing area.

Factors Influencing Surface Runoff Generation
Surface Saturation
Surface saturation was operationally defined as cases when at least two of the eight 45-cm wells recorded a water table within 5 cm of the soil surface. Sites S1, S2, and N1 were subject to surface saturation for many of the larger flow events during the study period (Table 4). At Site S1, there was a one-to-one correspondence between surface saturation and the occurrence of any surface runoff. Runoff was recorded at Site S2 under both surface-saturation and non-surface-saturation conditions; the largest runoff-producing events (>500 L) were all surface-saturating events at Site S2. There were four surface-saturating events that did not generate substantial runoff volumes (<300 L). The largest volumes at Site N1 were also restricted to surface-saturating events. Sites S3 and N2 each experienced surface saturation during only one event during the study period; Site N3 did not experience surface saturation at any point.

To quantify the effect of surface saturation, total runoff yields were calculated separately for surface-saturated and unsaturated events (Table 3). Sites S2 and N1 were the only sites with multiple surface-saturated and unsaturated runoff events. The ratio of total runoff yield generated under saturated conditions to total runoff yield generated under unsaturated conditions was 40 and 15 for Sites S2 and N1, respectively. These ratios are dramatic, yet they are somewhat biased because larger events tended to generate both surface saturation conditions and larger runoff yields. The effect of surface saturation can better be observed by comparing similarly sized events with differing saturation status. This analysis is presented in Fig. 4 , in which rainfall depths are plotted against peak rainfall intensity with symbol size proportional to runoff volume at Site S2. The effect of surface saturation can be assessed independent of rainfall depth by comparing the events with rainfall depths between 0.9 to 2.0 cm: the range of intermediate-sized events smaller than those that were exclusively saturation-inducing and larger than those that produced little runoff. There were 28 events in this size class, 10 of which caused surface saturation at Site S2. The non-surface-saturating events in this group had an average runoff yield of 0.1% with a maximum yield of 1.0%. In contrast, the surface saturation events in this size class had an average runoff yield of 8.5% with two events yielding over 20%.



View larger version (24K):
[in this window]
[in a new window]
 
Fig. 4. Rainfall depth and peak intensity with symbol size proportional to runoff volume at Site S2. Points with cross exhibited surface saturation, empty points did not. Vertical gray lines delimit 0.9- to 2.0-cm rainfall depth group discussed in text.

 
There was a close temporal correspondence during many events between surface saturation and runoff initiation. An example of the correspondence at Site S2 is displayed in Fig. 5 : the large rainfall depth, low intensity event on 24 Sept. 2001 that generated the largest runoff observed at this site during the study period. For this event, surface runoff did not begin at the flume until the well indicated surface saturation, despite relatively constant rainfall intensity.



View larger version (25K):
[in this window]
[in a new window]
 
Fig. 5. Temporal correspondence between runoff and surface saturation at Site S2 during 24 Sept. 2001 event.

 
A comparison between Sites S2 and S3 demonstrates the effect of the French drain on surface runoff generation. Upon installation of sensors at Site S3, a French drain was uncovered at about 40 to 75 cm running normal to the slope, located 1 m above the bank of piezometers and other sensors. The soil morphology at Sites S2 and S3 is very similar. Somewhat poorly drained soils extend 12 m upslope from the flume at S3. Surface saturation was only observed at this study during one event. We can use these two sites to distinguish the case when infiltration rate is limited by the hydraulic conductivity profile vs. when the rate is limited by surface saturation. When examining events that did not generate surface saturation at Site S2 (Table 4), Sites S2 and S3 had total runoff yields of 0.2 and 0.3%, respectively. These yields were low and none of these events caused substantial streamflow volumes. In fact, channel interception and near-stream ground water discharge can explain all of the stormflow for these events (McDonnell, 1990). For events that did cause surface saturation at Site S2, the total runoff yields for Site S2 and S3 were 8.0 and 1.8%, respectively. In general, these events were larger and more intense than those that did not result in surface saturation, so a greater mean yield was expected. The lower yields at Site S3 compared with Site S2 appear to be caused by the French drain, which prevents surface saturation.

The occurrence of surface saturation was strongly related to runoff generation at these sites. Surface saturation was observed almost exclusively at colluvial soil sites without subsurface drainage.

Management
Surface runoff from agricultural lands can be affected by surface cover, tillage, and other management practices. Management decisions were made by farmers in this on-farm study (Table 1). Therefore, site management was an uncontrolled variable in the experiment, which is a limitation of this study. Although the design of this study does not allow for a detailed analysis of management influences, broad trends were discernible.

There were three generalized cropping practices applied to the study sites: cultivated row crops, cultivated small grains and legumes, and managed grassland (Site S1). There was no discernible difference in runoff production between sites under cultivated row crops vs. those under cultivated small grains and legumes (Tables 1, 4). There were some minor management influences discernible within sites. For example, crop management may have caused the reduced runoff production at Site N3 during 2002 in comparison with 2000 and 2001. However, these differences were not substantial within the study (Table 4).

The greatest management variation was the grassland vegetation at Site S1. If this site had little runoff production, it would not have been clear whether the cause was management or other factors. However, this site was among the largest runoff-producing sites (along with S2). There was a one-to-one correspondence between surface runoff and the surface saturation at this site; the soil-dependent occurrence of surface saturation outweighed the runoff mitigating effect of the grassland vegetation.

The strongest runoff generation factor in this study, surface saturation, was independent of management practices and was related only to the presence of colluvial or residual soils and subsurface drainage at these sites (Table 4). For example, both Sites S2 and N3 were in small grains in 2002, yet three events with runoff yields >10% were recorded at Site S2 while there wasn't a single runoff-generating event observed at N3 (Table 4).

Antecedent Soil Water Content
Antecedent soil water content was associated with runoff production at these sites indirectly; the occurrence of surface saturation was dependent on antecedent soil water content for all except the largest events. The two large-volume, high-intensity events (>2.0 cm, >8.0 cm h–1 peak intensity) both caused surface saturation regardless of antecedent soil water content, though surface saturation was delayed. In all cases of surface saturation occurrence during smaller events (≤2.0 cm), there was a minimum antecedent surface soil water content of 0.35 m3 m–3 at the given site. The occurrence of a surface-saturated layer under these small rainfall depths was likely related to the presence of a capillary fringe before the event (Gillham, 1983; Novakowski and Gillham, 1988; Abdul and Gillham, 1989). The S1 and S2 sites exhibited the highest levels of surface soil water content during most of the study period. In general, surface soil water content levels are correlated to soil drainage class in this area (Henninger et al., 1976). These observations are consistent with positive correlations observed between soil water content levels and clay content, the presence of fragipan-perched water tables (Knuteson et al., 1989), and unsaturated flow of soil water to footslope and toeslope positions (Hewlett and Hibbert, 1963; Jackson, 1992).


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
The two contrasting soil groups found along the hillslopes (colluvial and residual) differ in subsurface morphological characteristics including the presence of a fragipan, the clay content of argillic horizons, and drainage class. The colluvial soil sites produced dramatically greater runoff than did the residual soil sites: overall runoff yields averaged 2.4% from the four colluvial sites and 0.01% from the two residual soil sites (Table 3), with greater runoff from the colluvial sites for all significant events (Table 4). The two residual soil sites (N2 and N3) generated minimal runoff for all events. These residual soils are the most common soil in the FD-36 watershed, comprising 76% of the area. Large, intense rainfall events may indeed generate runoff from the residual soils. For example, a 5-yr return period event that occurred just before the full establishment of the field equipment for this study generated runoff yields of 5 and 4% at Sites N2 and N3, respectively (though the yield at the colluvial sites was substantially greater: 23 and 18% at Sites S2 and S3, respectively).

We were not able to study winter conditions such as frozen soils or snowmelt events in this study because of equipment limitations. This is a major limitation of this study because hydrologic conditions during the winter (frozen soil, snowmelt events, high water table and soil water content) are favorable for large runoff events. For example, winter runoff from residual soils may have occurred during the study period. The results of this study suggest that nonwinter P management for the residual soils in this watershed should focus on rare, large events.

The methods common in current runoff indexing methods (such as curve number, soil hydrologic group, slope class) do not capture the important role of the subsurface soil properties and upslope soil characteristics that was observed in this study. For example, the method used in the Pennsylvania P index, a site assessment tool to rank agricultural fields by their potential for P export (Lemunyon and Gilbert, 1993), estimates that all of these sites are in the moderate class, except S1 which is low (Weld et al., 2002). This study was not designed to test the P index, and therefore we cannot make conclusions on the success of the P index at the monitoring locations. However, there are several implications of these results regarding indexing methods that warrant discussion. The flumes at Sites S1 and N1 are located within the same soil mapping unit and both frequently experienced surface saturation, but the S1 site generated substantially greater runoff yields than did the N1 site. This can be attributed to the contrasting soils upslope of these flumes (S1 primarily colluvial, N1 primarily residual). This comparison has important applied ramifications. For example, the Pennsylvania P index would estimate that Sites N1 and S1 have the same transport potential because they are the same distance to a stream and lie within the same soil mapping unit (therefore, identical Runoff Class and soil erodibility factor) (Gburek et al., 2000). However, the total runoff yield at Site S1 was nearly five times that of N1 (Table 3). For larger events, the contributing area for S1, composed of colluvial soils, was hydrologically active (as detected by distributed runoff indicators) while the N1 upslope residual soils were relatively inactive. To accurately estimate the runoff production of a site, upslope hydrologic characteristics should be integrated with other site-characteristic information.

While designing this study, we considered the rise of the fragipan-perched water table to the soil surface as the likely mechanism for surface saturation under the saturation-excess runoff generation mechanism. During the course of the study, we recognized and observed a second possibility, the perching of a water table within the argillic horizon. An argillic-perched transient water table may have developed as the advancing wetting front met sharply decreasing hydraulic conductivity rates at the Ap-Bt1 horizon boundary (p. 220–221 in Daniels and Hammer, 1992). A good example of this phenomenon occurred during the 12 June 2001 event at Site S1 (Table 4). The watershed was relatively dry before this event, stream baseflow was low, both piezometers were dry, and surface soil water content at the S1 site was 0.29 m3 m–3. Forty-eight hours after this event, the surface soil water content at this site was 0.43 m3 m–3. There is a monitored tile drain that outlets 5 m downslope of Flume 3. Before this event, this tile drain was not flowing (a rare occurrence during the study period). We deduced from this set of conditions that the subsurface soil horizons were not near field capacity before the event, and therefore a rapid rise of the water table caused by the wetting of a capillary fringe could not occur (Gillham, 1983). After only 2.0 cm of rainfall, surface saturation was recorded extending from 30-cm (piezometer response) to the soil surface (well data). Realistic estimates of porosity and antecedent profile water content are not consistent with the rise of a water table from the fragipan (60 cm) to the soil surface through a water-depleted profile with only 2.0 cm of rainfall. Rather, this free water must have been perched near the Ap-Bt1 boundary, which is normally 25 to 30 cm deep at these sites.

We may assume that there is a sharp hydraulic conductivity decrease at the Ap-Bt1 horizon boundary due to the strongly contrasting soil texture, bulk density, and structure at this boundary (Bonell, 1993). The degree of this drop is likely very strong in the colluvial soils and weaker in the residual soils because of the greater clay contents of the argillic horizons in the colluvial soils (37–44% in colluvial soils, 24% in residual soils). There is also likely a more general decrease in hydraulic conductivity with depth within the argillic horizon, because of both increasing clay content and increasing size of the subangular and angular blocky structure. A saturated layer similar to this argillic-perched water table was observed along the A-B horizon boundary in soils with thin A horizons in a field study on an agricultural watershed (Betson and Marius, 1969). It seems clear that surface saturation induced from the rise of a fragipan-perched water table is a case of saturation-excess runoff generation. However, the case of an argillic-perched water table inducing surface saturation has characteristics of both saturation-excess and infiltration-excess runoff generation.

All the colluvial soil sites were hydrologically active throughout the study period; the residual soil sites were hydrologically inactive. The lower portion of the colluvial footslope (Sites S1 and S2) produced the greatest runoff volumes and yields. The primary factor related to runoff production at the study sites was the occurrence of surface saturation, which was restricted to undrained colluvial soils. Antecedent soil water was also related to runoff production, both directly and because the occurrence of surface saturation was dependent on high antecedent soil water contents for smaller events. Agricultural management practices had minimal influence on runoff production at these sites in comparison with the strong influence of surface saturation. Nutrient management planning could be improved if surface runoff indexing methods, such as runoff class, were to better integrate information on the upslope and subsurface hydrologic properties (fragipans and high clay content argillic horizons) that differentiate the colluvial and residual soils in this study.


    ACKNOWLEDGMENTS
 
The authors would like to thank Jim Richards and Terry Troutman of the USDA-ARS Pasture Systems and Watershed Research Unit in Klingerstown, PA, for performing much of the design, installation, and data collection for the field study. We thank Jake Eckenrode of the Pennsylvania USDA-Nature Resources Conservation Service for conducting the soil survey. We also thank the three reviewers for their insightful comments and criticisms. Contribution from the U.S. Department of Agriculture, Agricultural Research Service, in cooperation with the Pennsylvania Agricultural Experiment Station, the Pennsylvania State University, University Park, PA, USA.

Received for publication March 3, 2003.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 




This article has been cited by other articles:


Home page
Soil Sci.Home page
H. E. Winzeler, P. R. Owens, B. C. Joern, J. J. Camberato, B. D. Lee, D. E. Anderson, and D. R. Smith
Potassium Fertility and Terrain Attributes in a Fragiudalf Drainage Catena
Soil Sci. Soc. Am. J., September 1, 2008; 72(5): 1311 - 1320.
[Abstract] [Full Text] [PDF]


Home page
J. Environ. Qual.Home page
Z. M. Easton, M. T. Walter, and T. S. Steenhuis
Combined Monitoring and Modeling Indicate the Most Effective Agricultural Best Management Practices
J. Environ. Qual., August 8, 2008; 37(5): 1798 - 1809.
[Abstract] [Full Text] [PDF]


Home page
J. Environ. Qual.Home page
A. N. Sharpley, P. J. A. Kleinman, A. L. Heathwaite, W. J. Gburek, J. L. Weld, and G. J. Folmar
Integrating Contributing Areas and Indexing Phosphorus Loss from Agricultural Watersheds
J. Environ. Qual., June 23, 2008; 37(4): 1488 - 1496.
[Abstract] [Full Text] [PDF]


Home page
J. Environ. Qual.Home page
A. N. Sharpley, P. J. A. Kleinman, A. L. Heathwaite, W. J. Gburek, G. J. Folmar, and J. P. Schmidt
Phosphorus Loss from an Agricultural Watershed as a Function of Storm Size
J. Environ. Qual., March 1, 2008; 37(2): 362 - 368.
[Abstract] [Full Text] [PDF]


Home page
J. Environ. Qual.Home page
R. W. McDowell
Phosphorus and sediment loss in a catchment with winter forage grazing of cropland by dairy cattle.
J. Environ. Qual., March 1, 2006; 35(2): 575 - 583.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (6)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Needelman, B. A.
Right arrow Articles by Kleinman, P. J. A.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Needelman, B. A.
Right arrow Articles by Kleinman, P. J. A.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Needelman, B. A.
Right arrow Articles by Kleinman, P. J. A.
Related Collections
Right arrow Watershed and Landscape Processes
Right arrow Runoff
Right arrow Surface Hydrology
Right arrow Nutrient Management


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