We considered hydrologic and chemical factors controlling P export from a 39.5-ha mixed land use watershed in east-central Pennsylvania, focusing our evaluation on watershed vulnerability to P loss. The spatial variations of P source factors, soil P, and P inputs from fertilizer and manure were evaluated. Distribution of Mehlich-3 soil P on a 30-m grid over the watershed showed that soil P varied with land use. Soils in wooded areas had low Mehlich-3 P (<30 mg kg" 1); grazed pasture had Mehlich-3 P values between 100 and 200 mg kg '; and cropped fields receiving manure and fertilizer applications were mostly >200 mg kg" 1. Phosphorus sources and transport controls on P loss were evaluated by examining in-stream P concentrations during storm hydrographs. Phosphorus concentrations decreased 50% downstream from headwaters to watershed outlet, and were more closely related to near-stream (within 60 m) distribution of high-P soils than to that of the whole watershed. This suggests that near-stream surface runoff and soil P are controlling P export from the watershed. Based on these findings, we modified the Phosphorus Index (P1), a useroriented tool developed by the NRCS-USDA to identify critical source areas controlling P export from agricultural watersheds. The modification separately evaluates P source and transport factors, and incorporates the hydrologic return period to describe contributing areas. The modified P1 was applied to the watershed to illustrate interactions between P source and transport processes controlling P export, and approaches for managing P loss.

Phosphorus (P) loss from agricultural watersheds is generally greater in storm rather than base flow. Although fundamental to P-based risk assessment tools, few studies have quantified the effect of storm size on P loss. Thus, the loss of P as a function of flow type (base and storm flow) and size was quantified for a mixed-land use watershed (FD-36; 39.5 ha) from 1997 to 2006. Storm size was ranked by return period (<1, 1-3, 3-5, 5-10, and >10 yr), where increasing return period represents storms with greater peak and total flow. From 1997 to 2006, storm flow accounted for 32% of watershed discharge yet contributed 65% of dissolved reactive P (DP) (107 g ha(-1) yr(-1)) and 80% of total P (TP) exported (515 g ha(-1) yr(-1)). Of 248 storm flows during this period, 93% had a return period of <1 yr, contributing most of the 10-yr flow (6507 m(3) ha(-1); 63%) and export of DP (574 g ha(-1); 54%) and TP (2423 g ha(-1); 47%). Two 10-yr storms contributed 23% of P exported between 1997 and 2006. A significant increase in storm flow DP concentration with storm size (0.09-0.16 mg L(-1)) suggests that P release from soil and/or area of the watershed producing runoff increase with storm size. Thus, implementation of P-based Best Management Practice needs to consider what level of risk management is acceptable.

Many source and transport factors control P loss from agricultural landscapes; however, little information is available on how these factors are linked at a watershed scale. Thus, we investigated mechanisms controlling P release from soil and stream sediments in relation to storm and baseflow P concentrations at four flumes and in the channel of an agricultural watershed. Baseflow dissolved reactive phosphorus (DRP) concentrations were greater at the watershed outflow (Flume 1; 0.042 mg L(-1)) than uppermost flume (Flume 4; 0.028 mg L(-1)). Conversely, DRP concentrations were greater at Flume 4 (0.304 mg L(-1)) than Flume 1 (0.128 mg L(-1)) during stormflow. Similar trends in total phosphorus (TP) concentration were also observed. During stormflow, stream P concentrations are controlled by overland flow-generated erosion from areas of the watershed coincident with high soil P. In-channel decreases in P concentration during stormflow were attributed to sediment deposition, resorption of P, and dilution. The increase in baseflow P concentrations downstream was controlled by channel sediments. Phosphorus sorption maximum of Flume 4 sediment (532 mg kg(-1)) was greater than at the outlet Flume 1 (227 mg kg(-1)). Indeed, the decrease in P desorption between Flumes 1 and 4 sediment (0.046 to 0.025 mg L(-1)) was similar to the difference in baseflow DRP between Flumes 1 and 4 (0.042 to 0.028 mg L(-1)). This study shows that erosion, soil P concentration, and channel sediment P sorption properties influence streamflow DRP and TP. A better understanding of the spatial and temporal distribution of these processes and their connectivity over the landscape will aid targeting remedial practices.

No abstract

No abstract

This study by the U.S. Geological Survey (USGS), in cooperation with the Agricultural Research Service (ARS), U.S. Department of Agriculture, compared multiple methods for estimating groundwater recharge and base flow (as a proxy for recharge) at sites in east-central Pennsylvania underlain by fractured bedrock and representative of a humid-continental climate. This study was one of several within the USGS GroundWater Resources Program designed to provide an improved understanding of methods for estimating recharge in the eastern United States. Recharge was estimated on a monthly and annual basis using four methods-(1) unsaturated-zone drainage collected in gravity lysimeters, (2) daily water balance, (3) water-table fluctuations in wells, and (4) equations of Rorabaugh. Base flow was estimated by streamflow-hydrograph separation using the computer programs PART and HYSEP. Estimates of recharge and base flow were compared for an 8-year period (1994-2001) coinciding with operation of the gravity lysimeters at an experimental recharge site (Masser Recharge Site) and a longer 34-year period (1968-2001), for which climate and streamflow data were available on a 2.8-square-mile watershed (WE-38 watershed). Estimates of mean-annual recharge at the Masser Recharge Site and WE-38 watershed for 1994-2001 ranged from 9.9 to 14.0 inches (24 to 33 percent of precipitation). Recharge, in inches, from the various methods was: unsaturated-zone drainage, 12.2; daily water balance, 12.3; Rorabaugh equations with PULSE, 10.2, or RORA, 14.0; and water-table fluctuations, 9.9. Mean-annual base flow from streamflowhydrograph separation ranged from 9.0 to 11.6 inches (21-28 percent of precipitation). Base flow, in inches, from the various methods was: PART, 10.7; HYSEP Local Minimum, 9.0; HYSEP Sliding Interval, 11.5; and HYSEP Fixed Interval, 11.6. Estimating recharge from multiple methods is useful, but the inherent differences of the methods must be considered when comparing results. For example, although unsaturatedzone drainage from the gravity lysimeters provided the most direct measure of potential recharge, it does not incorporate spatial variability that is contained in watershed-wide estimates of net recharge from the Rorabaugh equations or base flow from streamflow-hydrograph separation. This study showed that water-level fluctuations, in particular, should be used with caution to estimate recharge in low-storage fractured-rock aquifers because of the variability of water-level response among wells and sensitivity of recharge to small errors in estimating specific yield. To bracket the largest range of plausible recharge, results from this study indicate that recharge derived from RORA should be compared with base flow from the Local-Minimum version of HYSEP.

No abstract

Abstract:A ®eld study site was installed in east-central Pennsylvania to examine processes controlling groundwater recharge. It was instrumented to monitor climatic inputs, soil water dynamics and groundwater response. Characterization of the layered fractured bedrock underlying the site by rock coring, seismic surveys and interval packer testing showed consistencies between layer depths, fracture frequencies, seismic velocities and hydraulic conductivities. Monthly summaries of rainfall and percolate over two years showed that percolate rates were generally high and closely related to precipitation during the dormant season. During the growing season, however, the relationship became erratic with large variabilities occurring between individual lysimeter measurements. Eight dormant season rainfall events were examined in detail. Smaller events produced similar responses from 1 m deep percolate lysimeters. Approximately 10±15 mm of rain was required to initiate percolate, with the time delay in response dependent on how long it took this depth to accumulate; 5 to 6 mm of the rain was retained in storage, with the remainder becoming percolate. Larger rains, from 30±110 mm, caused correspondingly larger depths of percolate and larger water table responses, but generally similar patterns of site response. Groundwater at the site was typically about 6 m below the land surface during the dormant season. It responded 1±2 hours after the onset of percolate, and reached its maximum elevation anywhere from 4 to 16 hours after that, even though percolate was still occurring. Based on causative depth of recharge and amount of water level rise in wells, the speci®c yield of the aquifer was found to be of the order of 0 . 01. This value is characteristic of fracture geometry rather than matrix properties of the bedrock.