Received 16 November 2015; accepted 4 January 2016; published 7 January 2016
In 1996, over half of surveyed waters in the U.S. were nutrient impaired  . Nearly 20 years later, nutrient enrichment remains a major impairment to designated uses of fresh and coastal waters of the U.S.   , which have led to major initiatives to reduce nutrient losses from major watersheds, such as the Chesapeake Bay Watershed  and Mississippi River Basin  . As nonpoint sources are far more difficult to quantify and control, the contribution of excess nutrients to surface waters from agriculture, in particular intensive livestock and crop production, has received increased attention  .
The U.S. Environmental Protection Agency (USEPA) recently focuses on the potential for dust and feathers exhausted fans from poultry broiler houses to be a nonpoint source of nutrients to nearby surface waters  . With the backdrop of ongoing litigation between Oklahoma and poultry integrators headquartered in northwest Arkansas   , Region 6 USEPA (i.e., Arkansas, Louisiana, New Mexico, Oklahoma, and Texas) sampled soils and standing water around poultry production facilities in the Illinois River Watershed in northwest Arkansas and northeast Oklahoma in 2010 through 2013. The results of these inspections are used to encourage independent producers to implement management strategies to reduce nutrient runoff from their facilities  .
Producing 1.0 billion birds per year, Arkansas is third to Georgia and Alabama (1.3 and 1.1 billion birds∙yr−1, respectively) in broiler production  . The majority of this production is located in northwest Arkansas and northeast Oklahoma. To accommodate a large production density, broiler houses are constructed on compacted clay pads with a minimum of two adjacent houses per farm. Between each house is a drainage swale for rainfall runoff to be transported away from the houses. Grass is established in this area for erosion control and to facilitate water infiltration and reduce surface runoff volume and energy.
One of the few studies investigating dust emissions from poultry house operation is conducted by Burns et al.  in Kentucky, who report that a single poultry house can emit in excess of 1700 kg of dust annually. The study utilizes sensors on the fans to measure size and concentration of dust particles prior to being discharged from the house. Based on visual examinations of production facilities on aerial photos, dust deposition can extend 3 to 9 m from the fan outlet. Fans are operated for different purposes and at different times during a flock production cycle, which greatly affects the amount of dust exhausted by each fan. There are two modes of fan operation, desired in-house temperature and ventilation to provide required air quality needed for optimum bird health and weight gain.
While dust deposited adjacent to fans has been observed at several poultry facilities in northwest Arkansas by USEPA inspections, no studies since Burns et al.  have investigated the potential for fan dust to actually enrich runoff from adjacent to poultry houses with P and N. Given the water quality concerns in the Illinois River Watershed   and ongoing litigation affecting poultry producers and integrators   , there is an urgent need for more information on the extent and magnitude of the potential of ventilation fan dust to contribute nutrients to runoff from production facilities.
The objectives of this study were to quantify P and N in runoff from sites immediately adjacent to poultry house fans in order to evaluate the potential for P and N in dust emitted during bird growth to be transported off site during surface runoff. Using simulated rainfall-runoff adjacent to house ventilation fans, the collected data were intended to represent the greatest potential nutrient concentrations from the dust from a single flock of birds. The dust deposition pattern allowed for installation of simulated rainfall plots centered over the thickest dust deposits outside a poultry house. Rainfall simulations were conducted after each flock was removed for five consecutive flocks to measure nutrient runoff for different seasons over one year from both sidewall and tunnel fans. It was hypothesized that simulated rainfall-runoff P and N concentrations differed seasonally and between fan types.
2. Materials and Methods
2.1. Site Description
This study was conducted at the University of Arkansas’ Savoy broiler production facility (36˚7'43.41''N, 94˚18'22.84''W) in western Washington County, northwest Arkansas (Figure 1). This facility is a small-scale, commercial production facility located in the Ozark Highlands (i.e., Major Land Resource Area 116A)  and consists of a set of four broiler houses. The natural soils immediately surrounding this facility, which is located in a summit landscape position at an elevation of approximately 384 m, are mapped as Nixa very gravelly silt loam (loamy-skeletal, siliceous, active, mesic Glossic Fragiudults)  . The Nixa series soils are very slowly permeable soils that formed in cherty limestone  , thus natural rainfall is prone to surface runoff from the area surrounding the broiler houses. Runoff from the area flows to nearby upland headwater tributaries that flow into the Illinois River, which is approximately 3.3 km directly west of the study site.
Figure 1. Study area location.
2.2. Experimental Design
Rainfall-runoff simulations were conducted adjacent to house #4 at the facility (Figure 2) to evaluate the effects of fan type (i.e., sidewall ventilation fans versus tunnel ventilation fans) and fan operation time (i.e., season, warm months versus cool months) on runoff quantity and runoff water quality; thus the experiment was treated as a block design. Rainfall-runoff simulations were conducted adjacent to only one house to minimize variability in nutrient concentrations in the dust on the soil surface adjacent to each fan. Though the measured data were specific to this particular single house and this production facility, the results were assumed to be representative of other broiler houses with both sidewall and tunnel ventilation fans throughout the Ozark Highlands region.
2.3. Ventilation Fan Configuration and Operation
At the study site, sidewall fans 1 (S1) and 4 (S4) operate on a timed schedule at the beginning of each flock during cool weather (Figure 2(a)). Sidewall fans S1, S3, and S4 operate on a timed schedule at the beginning of each flock during warm weather. If the in-house temperature is 1.1˚C above the desired level, all sidewall fans, including sidewall fan S2, automatically operate to maintain house temperature and humidity. Sidewall-fan use dominates during cool weather. When the in-house temperature rises to 1.6˚C, 2.2˚C, 3.1˚C, or 3.3˚C above the desired level, tunnel ventilation fans T1, T2, T3, and T4, respectively, automatically start operating (Figure 2(b)). Prior to T4 operating, airflow is directed across the house. Once T4 begins operating, airflow is pulled down the house through cooling cells at the opposite end of the house.
Tunnel fan use is greater during warm weather and each tunnel fan (34,000 m3 min−1) exhausts air at an appreciably greater rate than a sidewall fan (20,400 m3∙min−1), resulting in more dust being exhausted from the house via tunnel than sidewall ventilation fans. Based on the amount of electricity used by each fan to maintain optimum in-house humidity and temperate, the volume of air exhausted by each fan was determined (Figure 3).
Figure 2. Location of simulated rainfall-runoff plots adjacent to sidewall ((a): S1, S2, S3, and S4) and tunnel ventilation fans ((b): T1, T2, T3, and T4). (a) Location of sidewall ventilation fan plots; (b) Location of tunnel ventilation fan plots.
Figure 3. Volume of air exhausted by sidewall and tunnel fans during each flock of birds in 2012 and 2013.
2.4. Surface Runoff Plot Installation
Identically sized plots (3 m2) were installed outside broiler house #4 adjacent to each fan (Figure 4). As the grassed areas adjacent to the fans are exposed to natural rainfall, structures were built to cover the plot areas to protect deposited dust from natural rainfall (Figure 4(a)). For each fan, plots were positioned within the high-volume dust trail and in areas with consistent slopes (~3%) to obtain uniform runoff (Figure 4(b)). Dust volumes and deposition uniformity were expected to differ at least slightly among plots due to variability in fan operating time. The runoff plots were covered to avoid disturbance by natural rainfall between flocks (i.e., when there were birds in the house). Roof structures were removed prior to each rainfall simulation study, and replaced immediately following rainfall. The only disturbance to the grass at the site was cutting grass to a 5-cm height following each rainfall simulation.
Rainfall simulation plots 1.5-m wide by 2-m long with 3% slope bounded by 0.2-cm thick and 15-cm wide metal frames, were installed in the center of the dust plume adjacent to each fan (Figure 4(b)). Runoff collection gutters at the downslope edge of the plots diverted runoff to a below-ground container, from which the runoff was pumped into a 190-L barrel using a peristaltic pump. The frames were designed to fit under the rainfall simulator so the areas enclosed would receive equal amounts of rainfall with a coefficient of uniformity of >90% and followed National Phosphorus Research Project and U.S. Department of Agriculture-Natural Resources Conservation Service protocols  -  .
2.5. Surface Runoff Plot Installation
Rainfall simulation studies were conducted following bird removal for five consecutive flocks over a one-year period. Simulations on sidewall fan plots were conducted between April 2012 and January 2013 and on tunnel fan plots between April 2013 and January 2014. The portable rainfall frame, which was centered and leveled over the plot, was constructed of hollow metal tubing, upon which a single fixed nozzle (TeeJetTM 1/2HH-SS-
Figure 4. Simulated rainfall-runoff plot borders (a) and protective shelter (b).
5OWSQ), designed for a flow rate of 210 mL∙s−1, was mounted central to the top of the frame and 3 m above the surface of the plot  . Windscreens were attached to three sides to restrict wind impact on rainfall uniformity. Water was supplied by a 25-mm diameter pipe from the farm’s well. Samples of simulated rainfall were collected prior to each rainfall simulation for chemical analysis.
A regulator maintained water pressure at 28.3 kPa (4.1 PSI) to produce raindrops of similar size to natural rainfall   . Rainfall distribution and volume were measured prior to each simulation study using five collection pans located in each corner and the center of the plot, equidistant from each other. Pressure was maintained at 28.3 kPa (4.1 PSI) and rainfall was collected for 5 minutes prior to measurement.
Rainfall simulations were conducted at an intensity of 7 cm∙hr−1, which is equivalent to a 2-year, 24-hour storm, and runoff was collected for a total of 30 minutes according to the National Phosphorus Research Project  . Rainfall simulations were conducted on the same day for all plots of a particular fan type.
2.6. Control Data
The chemical composition of runoff from pastures adjacent to the farm, which had not received fertilizer nor had been grazed in the prior 15 years, were used as unamended control observations for comparison  .
2.7. Chemical Analyses
Total runoff collected during the rainfall simulations was weighed to determine volume and subsampled for analyses. A 40-mL subsample was filtered (0.45 µm) immediately after collection in the field and stored at 4˚C until analyzed. The concentration of dissolved P (DP) was determined on filtered runoff by the colorimetric molybdenum-blue method of Murphy and Riley  and ammonia-N (NH3-N) and nitrate-N (NO3-N) were measured on a Lachat Instruments QuikChem 8500 flow injection analysis system (Milwaukee, WI). A 125-mL subsample was acidified with 12 drops of concentrated sulfuric acid for sample preservation. Total phosphorus (TP) and total nitrogen (TN) were determined on unfiltered runoff samples following persulfate/autoclave digestion  . Total solids (TS) content was determined gravimetrically on a 125-mL sample of unfiltered runoff after over-drying (105˚C) for 12 hours. Rainwater was sampled, prepared, and analyzed identically.
2.8. Flow-Weighted Concentrations
Flow-weighted concentrations were determined by calculating the load (i.e., event concentration multiplied by total volume from 30 minutes of runoff) for each fan from a rainfall simulation runoff event. Loads for all of the plots for the runoff event were summed and flow for all of plots and runoff events were summed. Finally, the sum of loads was divided by the sum of runoff to obtain flow-weighted concentrations for the event for sidewall and tunnel fan plots. The flow-weighted concentrations represent the worst-case scenario if a poultry house was located adjacent to a waterway.
2.9. Statistical Analyses
Separate one-factor analyses of variance (ANOVA), were conducted to evaluate the effects of season and fan location/type (i.e., sidewall and tunnel fans) on runoff water quality (i.e., DP, TP, NH3-N, NO3-N, TN, and TS concentrations). Since sidewall and tunnel fans did not operate at the same times of the year, the interaction effects of season and fan type could not be evaluated. All relationships are reported as statistically significant at the P < 0.05 level, unless noted otherwise. Wolfram Mathematica (online beta version 2014, Wolfram, Champaign, Illinois) was used for all statistical analyses.
3. Results and Discussion
3.1. Ventilation Fan Use
The volume of air exhausted from house #4 by sidewall fans S1, S2, S3, and S4 and tunnel fans T1, T2, T3, and T4 were determined from in house electricity use records over 2012 and 2013 (Figure 3). A greater volume of air was exhausted by the four sidewall (52 and 62%) than then the four tunnel (48 and 38%) fans in the mid-winter flocks of 2012 (24th February to 11th April) and 2013 (5th February to 18th March), respectively (Figure 3). In contrast, during the warmer summer months (June through August), tunnel fans exhausted a greater volume of air from house #4 (89 and 84%) than from the sidewall fans (11 and 16%) during 2012 and 2013, respectively (Figure 3).
On an annual basis, the four sidewall fans exhausted 35% and tunnel fans 65% of the total amount of air exhausted in 2012. However, in 2013 sidewall fans exhausted a greater proportion of air than in 2012, exhausting 43% of the annual volume of air with tunnel fans exhausting 57% (Figure 3). The greater tunnel fan use relative to sidewall fans in 2012 than in 2013, reflects the appreciably warmer June, July and August of 2012 (average high of 34˚C) than 2013 (average high of 29˚C from National Oceanic and Atmospheric Administration http://www.srh.noaa.gov/). The difference in air volume between fan type and over time of the year when birds are present in the house, influences the amount of dust exhausted and deposited adjacent to house ventilation fans.
3.2. Effect of Season
The concentration of P and N forms in runoff varied with time of the year that rainfall was applied for sidewall (Figure 5 and Figure 6) and tunnel fans (Figure 7 and Figure 8). For sidewall fans, P and N concentrations were lower (P < 0.01) in runoff from fans S2, S3, and S4 for the August rainfall simulation (9.0 and 56 mg∙L−1 for TP and TN, respectively) than in runoff for the other four rainfall simulation dates (2.5 and 16.4 mg∙L−1). This seasonal difference reflects the fact that sidewall fans were used less than tunnel fans in warmer, more humid months, when tunnel fans operate longer to more effectively achieve desired in-house temperature and humidity.
Figure 5. Average concentration of dissolved and total P in runoff from plots adjacent to sidewall ventilation fans as a function of season with error bars given.
This shift in fan use with season is reflected in greater P and N concentrations in runoff from the tunnel fans in May and August (8.9 and 86.1 mg∙L−1, respectively) compared with runoff from rainfall simulations in October, January, and April (2.9 and 31 mg∙L−1, respectively; Figure 7 and Figure 8). Clearly, the potential for dust emitted from fans can enrich the P and N concentration of runoff occurring immediately adjacent to the fans. This potential is directly related to fan use, which changes during the year, with seasonal fluctuations in in-house temperature and humidity triggering fan use.
3.3. Effect of Fan Location
The concentration of DP and TP in runoff from sidewall fan S1 tended to be greater than from S2, S3, and S4 for each simulations (Figure 5). Unlike P, N did not vary consistently among fans (Figure 6). Sidewall fan S1 operates longer than the other sidewall fans, irrespective of season, with only fan S1 operating for the first week of a new flock, when birds occupy a quarter of the house. Thus, the amount of dust emitted from fan S1 is expected to be greater than from the other three fans, which related directly to the nutrient runoff trends observed.
3.4. Flow-Weighted Concentrations
Concentrations of P and N were converted to a flow-weighted basis within rainfall simulation events by calculating loads at each fan, summing the loads, then dividing the total by the total flow for all runoff plots (Figure 9 and Figure 10). The average for runoff concentrations from sidewall and tunnel fans allowed a comparison
Figure 6. Average concentration of ammonia-N, nitrate-N, and total N in runoff from plots adjacent to sidewall ventilation fans as a function of season with error bars given.
among fans and with nutrient runoff from nearby pastures that had not been fertilized or grazed in 15 years and pasture that had received poultry litter (Table 1). Fluctuations in average flow-weighted concentrations of P and N forms in runoff from plots adjacent to sidewall and tunnel fans were consistent with changes in fan use (Figure 9 and Figure 10).
The flow-weighted concentration of DP in runoff from plots adjacent to sidewall (4.9 mg∙L−1) and tunnel fans (2.7 mg∙L−1) was within the range of values measured in runoff from area pastures that received varying amounts of poultry litter (0.7 to 33 mg∙L−1; Table 1). Total P concentration in runoff from sidewall (8.1 mg∙L−1) and tunnel fans (4.7 mg∙L−1) tended to be at the lower limits of the range measured in fertilized pasture runoff (3.8 to 37 mg∙L−1; Table 1). As expected, runoff from plots impacted by fan dust were appreciably greater than concentrations measured in runoff from pasture that had not been fertilized or grazed for at least 15 years (0.07 and 0.10 mg∙L−1; Table 1).
Figure 7. Average concentration of dissolved and total P in runoff from plots adjacent to tunnel ventilation fans as a function of season with error bars given.
Table 1. Phosphorus runoff from plots treated with poultry litter and from unfertilized ungrazed pasture.
Figure 8. Average concentration of ammonia-N, nitrate-N, and total N in runoff from plots adjacent to tunnel ventilation fans as a function of season with error bars given.
4. Summary and Conclusions
The calculation of P and N loss in runoff from 3-m2 plots using simulated rainfall (7 cm∙hr−1 for 30 min) does not represent loss on a field basis under natural rainfall; thus loads are not calculated  . Use of the plot and simulated rainfall in this research does, however, allow a comparison of nutrient runoff adjacent to ventilation fans of poultry broiler houses under a constant rain fall event. In this research, fan type and time of year are evaluated for their independent effects on simulated rainfall-runoff water quality.
Tunnel fan use is greater during warm weather and each tunnel fan exhausts at a rate that is 2.5 times greater than a sidewall fan, resulting in more dust being exhausted from the house during August, the peak period of tunnel fan use. The current research is conducted to obtain baseline data to evaluate the potential for P and N enrichment of runoff occurring from around ventilation fans of poultry broiler houses. Additional testing will be required to determine if there are statistically significant differences in runoff concentrations among fan location and type during the same growing season.
Phosphorus concentrations in runoff from fan plots are also compared to data from area pastures to which poultry litter has been applied. The concentration of DP and TP from pastures receiving litter ranges from 0.7 to
Figure 9. Flow-weighted concentrations of P and N from runoff plots adjacent to sidewall ventilation fans with error bars given.
Figure 10. Flow-weighted concentrations of P and N from runoff plots adjacent to tunnel ventilation fans with error bars given.
33 mg∙L−1 and 3.8 to 37 mg∙L−1, respectively. In contrast, concentrations from fan plots ranges from 0.2 to 12 mg∙L−1 for DP and 1.0 to 26 mg∙L−1 for TP. Based on this comparison, runoff concentrations from the fan plots are similar to that observed in runoff from fields that have received poultry litter.
This study provides baseline data for the potential for P and N runoff adjacent to poultry house ventilation fans. Although the subject facility is not discharging into a waterway, it is important to recognize that runoff concentrations for each fan plot, regardless of the season, are consistently above current regulated limits for wastewater treatment plant discharge into the Illinois River (1.0 mg∙TP∙L−1). If a facility is proximate to a waterway, poultry house fan dust runoff concentrations at the nutrient source can be elevated and pose a threat to the quality of receiving waters.
This research was supported by a grant from the USGS 104B program funding.
 United States Environmental Protection Agency (1996) Environmental Indicators of Water Quality in the United States. EPA 841-R-96-002. U.S. EPA, Office of Water (4503F), U.S. Govt. Printing Office, Washington DC.
 Dale, V.H., Kling, C.L., Meyer, J.L., Sanders, J., Stallworth, H., Armitage, T., Wangsness, D., Bianchi, T.S., Blumberg, A., Boynton, W., Conley, D.J., Crumpton, W., David, M.B., Gilbert, D., Howarth, R.W., Lowrance, R., Mankin, K.R., Opaluch, J., Paerl, H.W., Reckhow, K., Sharpley, A.N., Simpson, T.W., Snyder, C. and Wright, D. (2010) Hypoxia in the Northern Gulf of Mexico. Springer Series on Environmental Management. Springer Science, New York.
 Schindler, D.W., Hecky, R.E., Findlay, D.L., Stainton, M.P., Parker, B.R., Paterson, M.J., Beaty, K.G., Lyng, M. and Kasian, S.E.M. (2008) Eutrophication of Lakes Cannot Be Controlled by Reducing Nitrogen Input: Results of a 37-Year Whole-Ecosystem Experiment. Proceedings of the National Academy of Sciences of the United States of America, 106, 11254-11258.
 Dubrovsky, N.M., Burow, K.R., Clark, G.M., Gronberg, J.M., Hamilton, P.A., Hitt, K.J., Mueller, D.K. and Munn, M.D. (2010) The Quality of Our Nation’s Waters—Nutrients in the Nation’s Streams and Groundwater, 1992-2004: United States Geological Survey Circular 1350. 174 p.
 National Research Council (2008) Nutrient Control Actions for Improving Water Quality in the Mississippi River Basin and Northern Gulf of Mexico. Committee on the Mississippi River and the Clean Water Act: Scientific, Modeling and Technical Aspects of Nutrient Pollutant Load Allocation and Implementation. National Research Council, Washington, District of Columbia. 75 p.
 United States Environmental Protection Agency (2010) Guidance for Federal Land Management in the Chesapeake Bay Watershed. Chapter 2: Agriculture. EPA 841-R-10-002. U.S. EPA, Nonpoint Source Pollution, Office of Wetlands, Oceans, and Watersheds. Washington DC. 247 p.
 United States Environmental Protection Agency (2011) Findings of Violation and Order of Compliance, in the Matter of Mrs. Lois Alt, Old Fields, West Virginia. Administrative Order CWA-03-2012-0023DN. U.S. EPA Region 3, Philadelphia. 12p.
 State of Oklahoma vs. Tyson Foods Inc. (2009) The Court Should Exclude Nutrient Management Plans from Watersheds Other than the Illinois River Watershed.
 Burns, R., Li, H., Moody, L., Xin, H., Gates, R., Overhults, D. and Earnest, J. (2008) Quantification of Particulate Emissions from Broiler Houses in the Southeastern United States. Proceedings of the 2008 Livestock Environment VIII Conference, American Society of Agricultural and Biological Engineers, Publication 701P0408.
 Herron, S., Sharpley, A.N., Watkins, S. and Daniels, M. (2012) Poultry Litter Management in the Illinois River Watershed of Arkansas and Oklahoma. Cooperative Extension Service, Division of Agriculture, University of Arkansas. Fact Sheet FSA 9535, 4 p.
 Sharpley, A.N., Richards, P., Herron, S. and Baker, D. (2012) Comparison between Litigated and Voluntary Nutrient Management Strategies. Journal of Soil & Water Conservation, 67, 442-450.
 Soil Survey Staff, Natural Resources Conservation Service, United States Department of Agriculture (2015) Official Soil Series Descriptions.
 Sharpley, A.N. and Kleinman, P.J.A. (2003) Effect of Rainfall Simulator and Plot Scale on Overland Flow and Phosphorus Transport. Journal of Environmental Quality, 32, 2172-2179.
 Southern Extension-Research Activity-17 (2015) National Phosphorus Research Project for Simulated Rainfall—Surface Runoff Studies. Southern Extension-Research Activity-Information Exchange Group 17, Virginia Tech. University, Blacksburg, VA.
 Spaeth, K., Weltz, M., Pierson, F., Blackburn, W.H., Fox, D., Merz, D., Shaver, P., Brackley, G., DeGarmo, H., Whited, M., Flanagan, M., Hood, K. and Franks, C. (1995) Small Plot Rainfall Simulation: Background and Procedures. Technical Note Rangeland Hydrology. US Department of Agriculture, Natural Resources Conservation Service, Washington DC.
 Humphry, J.B., Daniel, T.C., Edwards, D.R. and Sharpley, A.N. (2002) A Portable Rainfall Simulator for Plot-Scale Runoff Studies. Applied Engineering in Agriculture, 18, 199-204.
 Daniel, T.C., Stead, J., DeLaune, P., Sharpley, A.N. and Verkler, T. (2008) Edge-of-Field Water Quality Monitoring from Various Management Practices in the Ozark Highlands. Arkansas Natural Resources Commission, 319 Project 05-1300. Final Report, 189 p.
 Murphy, J. and Riley, J.R. (1962) A Modified Single Solution Method for the Determination of Phosphate in Natural Waters. Analytica Chimica Acta, 27, 31-36.
 Patton, C.J. and Kryskalla, J.R. (2003) Methods of Analysis by the US Geological Survey National Water Quality Laboratory: Evaluation of Alkaline Persulfate Digestion as an Alternative to Kjeldahl Digestion for the Determination of Total and Dissolved Nitrogen and Phosphorus in Water. US Geological Survey, Water Resources Investigations Report 03-4174. US Geological Survey, Branch of Information Services, Federal Center, Denver, 33 p.
 DeLaune, P.B., Moore Jr., P.A., Carmen, D., Sharpley, A.N., Haggard, B.E. and Daniel, T.C. (2004) Factors Affecting Phosphorus Loss from Poultry Litter: Development of Phosphorus Source Component. Journal of Environmental Quality, 33, 2183-2191.
 Edwards, D.R. and Daniel, T.C. (1993) Effects of Poultry Litter Application Rate and Rainfall Intensity on Quality of Runoff from Fescue Grass Plots. Journal of Environmental Quality, 22, 361-365.
 Harmel, R.D., Torbert, H.A., Haggard, B.E., Haney, R. and Dozier, M. (2004) Water Quality Impacts of Converting to a Poultry Litter Fertilization Strategy. Journal of Environmental Quality, 33, 2229-2242.
 Kleinman, P.J.A., Sharpley, A.N., Bryant, R., Czymmek, K. and Ketterings, Q. (2001) National Phosphorus Runoff Project: New York. USDA-ARS, Pasture Systems and Watershed Management Research Unit, University Park, PA and Department of Crop and Soil Sciences, Cornell University, Ithaca, NY, 7 p.
 Sistani, K.R., Bolster, C.H., Way, T.R., Torbert, H.A., Pote, D.H. and Watts, D.B. (2010) Influence of Poultry Litter Application Methods on the Longevity of Nutrient and E. coli in Runoff from Fescue Pasture. Water, Air, & Soil Pollution, 206, 3-12.
 Smith, D.R., Moore Jr., P.A., Miller, D.M., Haggard, B.E. and Daniel, T.C. (2004) Decreasing Phosphorus Runoff Losses from Land-Applied Poultry Litter with Dietary Modifications and Alum Addition. Journal of Environmental Quality, 33, 2210-2216.