Simulating endosulfan transport in runoff from cotton fields in Australia with the GLEAMS model

Endosulfan (6,7,8,9,10,10-hexachloro-1,5,5a,6,9,9a-hexahydro-6,9methano-2,4,3-benzodioxathiepin 3-oxide), a pesticide that is highly toxic to aquatic organisms, is widely used in the cotton (Gossypium hirsutum L.) industry in Australia and is a risk to the downstream riverine environment. We used the GLEAMS model to evaluate the effectiveness of a range of management scenarios aimed at minimizing endosulfan transport in runoff at the field scale. The field management scenarios simulated were (i) Conventional, bare soil at the beginning of the cotton season and seven irrigations per season; (ii) Improved Irrigation, irrigation amounts reduced and frequency increased to reduce runoff from excess irrigation; (iii) Dryland, no irrigation; (iv) Stubble Retained, increased soil cover created by retaining residue from the previous crop or a specially planted winter cover crop; and (v) Reduced Sprays, a fewer number of sprays. Stubble Retained was the most effective scenario for minimizing endosulfan transport because infiltration was increased and erosion reduced, and the stubble intercepted and neutralized a proportion of the applied endosulfan. Reducing excess irrigation reduced annual export rates by 80 to 90%, but transport in larger storm events was still high. Reducing the number of pesticide applications only reduced transport when three or fewer sprays were applied. We conclude that endosulfan transport from cotton farms can be minimized with a combination of field management practices that reduce excess irrigation and concentration of pesticide on the soil at any point in time; however, discharges, probably with endosulfan concentrations exceeding guideline values, will still occur in storm events.     

Phosphorus budgets and riverine phosphorus export in northwestern Ohio watersheds

Phosphorus (P) budgets for large watersheds are often used to predict trends in riverine P export. To test such predictions, we calculated annual P budgets for 1975-1995 for soils of the Maumee and Sandusky watersheds of northwestern Ohio and compared them with riverine P export from these watersheds. Phosphorus inputs to the soils include fertilizers, manure, rainfall, and sludge while outputs include crop removal and nonpoint-source export via rivers. Annual P inputs decreased due to reductions in fertilizer and manure inputs. Annual outputs increased due to increasing crop yields. Net P accumulation decreased from peak values of 13.4 and 9.5 kg P ha(-1) yr(-1) to 3.7 and 2.6 kg P ha(-1) yr(-1) for the Maumee and Sandusky watersheds, respectively. Thus, P budget analysis suggests that riverine P export should have increased throughout the study period, with smaller increases during more recent years. However, detailed water quality studies show that riverine export of total phosphorus (TP) has decreased by 25 to 40% and soluble reactive phosphorus (SRP) by 60 to 89%, both due primarily to decreases from nonpoint sources. We suggest that these decreases are associated with farmers' adoption of practices that minimize transport of recently applied P fertilizer and of sediments via surface runoff, coupled with changes in winter weather conditions. In comparison with most Midwestern watersheds, rivers draining these watersheds have high unit area yields of TP, low unit area yields of SRP, and high ratios of nonpoint source- to point source-derived P.     

Fate and toxicity of endosulfan in Namoi River water and bottom sediment

Endosulfan (6,7,8,9,10,10,-hexachloro-1,5,5a,6,9,9a-hexahydro-6,9-methano-2,4,3-benzodioxathiepine-3-oxide) sorption (standardized to 1% total organic carbon and dry weight) was significantly (P < 0.05) more concentrated on the large (>63 microm) particle fraction compared with smaller size fractions (<5 microm and 5-24 microm) of bottom sediments from the Namoi River, Australia. Following completion of the particle size fractionation (6 to 12 wk) and a sediment toxicity assessment (2 wk), the sediments showed large decreases in concentrations of alpha-endosulfan that coincided with an increase in endosulfan sulfate concentrations and minimal changes in beta-endosulfan concentrations. In the Namoi River, similar patterns were observed in the composition of total endosulfan in monthly measurements of bottom sediments and in passive samplers placed in the water column following runoff from cotton (Gossypium hirsutum L.) fields. The toxicity of endosulfan sulfate in river water indicated by the nymphs of the epibenthic mayfly Jappa kutera, was more persistent than the alpha- and beta-endosulfan parent isomers due to its longer half-life. This suggests that endosulfan sulfate would contribute most to previously observed changes in population densities of aquatic biota. Measured concentrations of total endosulfan in river water of up to 4 microg L(-1) following storm runoff, exceed the range of the 96-h median lethal concentration (LC50) values in river water for both alpha-endosulfan (LC50 = 0.7 microg L(-1); 95% confidence interval [CI] = 0.5 to 1.1) and endosulfan sulfate (LC50 = 1.2 microg L(-1); 95% CI = 0.4 to 3.3). In contrast, the 10-d LC50 value for total endosulfan in the sediment toxicity test (LC50 = 162 microg kg(-1); 95% CI = 120 to 218 microg kg(-1)) was more than threefold higher than the highest measured concentration of total endosulfan in field samples of bottom sediment (48 microg kg(-1)). This suggests that pulse exposures of endosulfan in the water column following storm runoff may be more acutely toxic to riverine biota than in contaminated bottom sediment.     

Sources of nitrate yields in the Mississippi River Basin

Riverine nitrate N in the Mississippi River leads to hypoxia in the Gulf of Mexico. Several recent modeling studies estimated major N inputs and suggested source areas that could be targeted for conservation programs. We conducted a similar analysis with more recent and extensive data that demonstrates the importance of hydrology in controlling the percentage of net N inputs (NNI) exported by rivers. The average fraction of annual riverine nitrate N export/NNI ranged from 0.05 for the lower Mississippi subbasin to 0.3 for the upper Mississippi River basin and as high as 1.4 (4.2 in a wet year) for the Embarras River watershed, a mostly tile-drained basin. Intensive corn (Zea mays L.) and soybean [Glycine max (L.) Merr.] watersheds on Mollisols had low NNI values and when combined with riverine N losses suggest a net depletion of soil organic N. We used county-level data to develop a nonlinear model ofN inputs and landscape factors that were related to winter-spring riverine nitrate yields for 153 watersheds within the basin. We found that river runoff times fertilizer N input was the major predictive term, explaining 76% of the variation in the model. Fertilizer inputs were highly correlated with fraction of land area in row crops. Tile drainage explained 17% of the spatial variation in winter-spring nitrate yield, whereas human consumption of N (i.e., sewage effluent) accounted for 7%. Net N inputs were not a good predictor of riverine nitrate N yields, nor were other N balances. We used this model to predict the expected nitrate N yield from each county in the Mississippi River basin; the greatest nitrate N yields corresponded to the highly productive, tile-drained cornbelt from southwest Minnesota across Iowa, Illinois, Indiana, and Ohio. This analysis can be used to guide decisions about where efforts to reduce nitrate N losses can be most effectively targeted to improve local water quality and reduce export to the Gulf of Mexico.     

Off-site movement of endosulfan from irrigated cotton in New South Wales

The fate and transport of endosulfan (6,7,8,9,10,10-hexachloro-1,5, 5a,6,9,9a-hexahydro-6,9-methano-2,4,3-benzodioxathiepin 3-oxide) applied to cotton (Gossypium hirsutum L.) fields were studied throughout three consecutive years on two selected locations in New South Wales (Australia). Rates of dissipation from foliage and soil, volatilization from the field, and transport of residues in irrigation and/or storm runoff waters were measured in order to estimate a total field balance. Dissipation of endosulfan from both foliage and soil is best explained by a two-phase process rather than by a first-order decay. Half-lives of total endosulfan toxic residues (alpha- and beta-endosulfan and the sulfate product) in the first phase were 1.6 d in foliage and 7.1 d in soil, and could be explained by the rapid volatilization of the parent isomers in the first 5 d (up to 70% of endosulfan volatilizes). In the second phase, half-lives were 9.5 d in foliage and 82 d in soil, mostly due to the persistence of the sulfate product. Concentration of endosulfan residues in runoff water varied from 45 to 2.5 microg L(-1) depending on the residue levels present on field soil at the time of the irrigation or storm events. These in turn are related to the total amounts applied, the cotton canopy cover at application, and the time since last spraying. Most of the endosulfan in runoff was found in the water phase (80%), suggesting it was bound to colloidal matter. Total endosulfan residues in runoff for a whole season accounted for no more than 2% of the pesticide applied on-field.     

Phosphorus concentrations, loads, and sources within the Illinois River drainage area, northwest Arkansas, 1997-2008

In the Ozark Highlands and across the United States, effluent phosphorus (P) sources often have a profound impact on water column concentrations and riverine transport. This study evaluated (i) annual P loads at the Illinois River at Arkansas Highway 59 from calendar year 1997 through 2008, (ii) the relative contribution of effluent P sources to annual riverine P transport, (iii) longitudinal gradients in water column P concentrations downstream from several wastewater treatment plant effluent discharges, and (iv) changes in monthly P loads over the last decade. This study showed that annual P loads have ranged from 64,000 kg to over 426,000 kg and that P transport was positively correlated to hydrology (i.e., the amount of water delivered downstream). The relative contribution of P inputs from municipal facilities has decreased from 40% of the annual P load at the Illinois River at Arkansas Highway 59 to < 15% in recent years. Elevated P concentrations during base flow conditions were traced 45 river km upstream to one municipal effluent discharge, but all effluent discharges influenced P concentrations in the receiving streams. Most important, flow-adjusted monthly P loads showed two distinct trends over time. Flow-adjusted loads significantly increased from 1997 through 2002 and significantly decreased from 2002 through 2008. The concentrations and transport of P within the Illinois River drainage area are significantly decreasing from all the watershed management changes that have occurred, and monitoring should continue to determine if this decrease continues at the same rate over the next several years.     

Ecological Linkages Between Headwaters and Downstream Ecosystems: Transport of Organic Matter,Invertebrates, and Wood Down Headwater Channels1

Abstract: Headwater streams make up a large proportion of the total length and watershed area of fluvial networks, and are partially characterized by the large volume of organic matter (large wood, detritus, and dissolved organic matter) and invertebrate inputs from the riparian forest, relative to stream size. Much of those inputs are exported to downstream reaches through time where they potentially subsidize river communities. The relative rates, timing, and conversion processes that carry inputs from small streams to downstream reaches are reasonably well quantified. For example, larger particles are converted to smaller particles, which are more easily exported. Also, dissolved organic matter and surface biofilms are converted to larger particles which can be more easily intercepted by consumers. However, the quality of these materials as it affects biological activity downstream is not well known, nor is the extent to which timing permits biological use of those particles. These ecological unknowns need to be resolved. Further, land uses may disrupt and diminish material transport to downstream reaches by removing sources (e.g., forest harvest), by affecting transport and decomposition processes (e.g., flow regulation, irrigation, changes in biotic communities), and by altering mechanisms of storage within headwaters (e.g., channelization). We present conceptual models of energy and nutrient fluxes that outline small stream processes and pathways important to downstream communities, and we identify informational gaps that, if filled, could significantly advance the understanding of linkages between headwater streams and larger rivers. The models, based on empirical evidence and best professional judgment, suggest that navigable waters are significantly influenced by headwater streams through hydrological and ecological connectivities, and land use can dramatically influence these natural connectivities, impacting downstream riverine ecosystems.     

The impact of agricultural runoff on the quality of two streams in vegetable farm areas in Ghana

A study of two small streams at Akumadan and Tono, Ghana, was undertaken during the rain and dry season periods between February 2005 and January 2006 to investigate the impact of vegetable field runoff on their quality. In each stream we compared the concentration of current-use pesticides in one site immediately upstream of a vegetable field with a second site immediately downstream. Only trace concentrations of endosulfan and chlorpyrifos were detected at both sites in both streams in the dry season. In the wet season, rain-induced runoff transported pesticides into downstream stretches of the streams. Average peak levels in the streams themselves were 0.07 microg L(-1) endosulfan, 0.02 microg L(-1) chlorpyrifos (the Akumadan stream); 0.04 microg L(-1) endosulfan, 0.02 microg L(-1) chlorpyrifos (the Tono stream). Respective average pesticide levels associated with streambed sediment were 1.34 and 0.32 microg kg(-1) (the Akumadan stream), and 0.92 and 0.84 microg kg(-1) (the Tono stream). Further investigations are needed to establish the potential endosulfan and chlorpyrifos effects on aquatic invertebrate and fish in these streams. Meanwhile measures should be undertaken to reduce the input of these chemicals via runoff.     

Comparison of field-scale herbicide runoff and volatilization losses: an eight-year field investigation

An 8-yr study was conducted to better understand factors influencing year-to-year variability in field-scale herbicide volatilization and surface runoff losses. The 21-ha research site is located at the USDA-ARS Beltsville Agricultural Research Center in Beltsville, MD. Site location, herbicide formulations, and agricultural management practices remained unchanged throughout the duration of the study. Metolachlor [2-chloro--(2-ethyl-6-methylphenyl)--(2-methoxy-1-methylethyl) acetamide] and atrazine [6-chloro--ethyl--(1-methylethyl)-1,3,5-triazine-2,4-diamine] were coapplied as a surface broadcast spray. Herbicide runoff was monitored from a month before application through harvest. A flux gradient technique was used to compute volatilization fluxes for the first 5 d after application using herbicide concentration profiles and turbulent fluxes of heat and water vapor as determined from eddy covariance measurements. Results demonstrated that volatilization losses for these two herbicides were significantly greater than runoff losses ( < 0.007), even though both have relatively low vapor pressures. The largest annual runoff loss for metolachlor never exceeded 2.5%, whereas atrazine runoff never exceeded 3% of that applied. On the other hand, herbicide cumulative volatilization losses after 5 d ranged from about 5 to 63% of that applied for metolachlor and about 2 to 12% of that applied for atrazine. Additionally, daytime herbicide volatilization losses were significantly greater than nighttime vapor losses ( < 0.05). This research confirmed that vapor losses for some commonly used herbicides frequently exceeds runoff losses and herbicide vapor losses on the same site and with the same management practices can vary significantly year to year depending on local environmental conditions.     

Spray drift of pesticides arising from aerial application in cotton

This paper presents results from field studies carried out during the 1993-1998 Australian cotton (Gossypium hirsutum L.) seasons to monitor off-target droplet movement of endosulfan (6,7,8,9,10,10-hexachloro-1,5,5a,6,9,9a-hexahydro-6,9-methano-2,4,3-benzodioxathiepin 3-oxide) insecticide applied to a commercial cotton crop. Averaged over a wide range of conditions, off-target deposition 500 m downwind of the field boundary was approximately 2% of the field-applied rate with oil-based applications and 1% with water-based applications. Mean airborne drift values recorded 100 m downwind of a single flight line were a third as much with water-based application compared with oil-based application. Calculations using a Gaussian diffusion model and the U.S. Spray Drift Task Force AgDRIFT model produced downwind drift profiles that compared favorably with experimental data. Both models and data indicate that by adopting large droplet placement (LDP) application methods and incorporating crop buffer distances, spray drift can be effectively managed.     

A spatial analysis of phosphorus in the Mississippi river basin

Phosphorus (P) in rivers in the Mississippi River basin (MRB) contributes to hypoxia in the Gulf of Mexico and impairs local water quality. We analyzed the spatial pattern of P in the MRB to determine the counties with the greatest January to June P riverine yields and the most critical factors related to this P loss. Using a database of P inputs and landscape characteristics from 1997 through 2006 for each county in the MRB, we created regression models relating riverine total P (TP), dissolved reactive P (DRP), and particulate P (PP) yields for watersheds within the MRB to these factors. Riverine yields of P were estimated from the average concentration of each form of P during January to June for the 10-yr period, multiplied by the average daily flow, and then summed for the 6-mo period. The fraction of land planted in crops, human consumption of P, and precipitation were found to best predict TP yields with a spatial error regression model ( = 0.48, = 101). Dissolved reactive P yields were predicted by fertilizer P inputs, human consumption of P, and precipitation in a multiple regression model ( = 0.42, = 73), whereas PP yields were explained by crop fraction, human consumption of P, and soil bulk density in a spatial error regression model ( = 0.49, = 61). Overall, the Upper Midwest's Cornbelt region and lower Mississippi basin had the counties with the greatest P yields. These results help to point out specific areas where agricultural conservation practices that reduce losses to streams and rivers and point source P removal might limit the intensity or spatial occurrence of Gulf of Mexico hypoxia and improve local water quality.     

Pesticide transport pathways from a sloped Litchi orchard to an adjacent tropical stream as identified by hydrograph separation

This study was performed to identify the transport pathways of pesticides from a sloped litchi ( Sonn.) orchard to a nearby stream based on a three-component hydrograph separation (baseflow, interflow, surface runoff). Dissolved silica and electrical conductivity were chosen as representative tracers. During the study period (30 d), 0.4 and 0.01% of the applied mass of atrazine and chlorpyrifos, respectively, were detected in the stream after 151 mm of rainfall. Baseflow (80-96%) was the dominant hydrological flow component, followed by interflow (3-18%) and surface runoff (1-7%). Despite its small contribution to total discharge, surface runoff was the dominant atrazine transport pathway during the first days after application because pesticide concentrations in the surface runoff flow component declined quickly within several days. Preferential transport with interflow became the dominant pathway of atrazine. Because chlorpyrifos was detected in the stream water only twice, it was not included in the hydrograph separation. A feature of the surface runoff pathway was the coincidence of pesticide and discharge peaks. In contrast, peak concentrations of pesticides transported by interflow occurred during the hydrograph recession phases. Stormflow generation and pesticide transport depended on antecedent rainfall. The combination of high-resolution pesticide concentration measurements with a three-component hydrograph separation has been shown to be a suitable method to identify pesticide transport pathways under tropical conditions.     

Antibiotic Transport via Runoff and Soil Loss

Research has verified the occurrence of veterinary antibiotics in manure, agricultural fields, and surface water bodies, yet little research has evaluated antibiotic runoff from agricultural fields. The objective of this study was to evaluate the potential for agricultural runoff to contribute antibiotics to surface water bodies in a worst-case scenario. Our hypothesis was that there would be significant differences in antibiotic concentrations, partitioning of losses between runoff and sediment, and pseudo-partitioning coefficients (ratio of sediment concentration to runoff concentration) among antibiotics. An antibiotic solution including tetracycline (TC), chlortetracycline (CTC), sulfathiazole (STZ), sulfamethazine (SMZ), erythromycin (ERY), tylosin (TYL), and monensin (MNS) was sprayed on the soil surface 1 h before rainfall simulation (average intensity = 60 mm h(-1) for 1 h). Runoff samples were collected continuously and analyzed for aqueous and sediment antibiotic concentrations. MNS had the highest concentration in runoff, resulting in the highest absolute loss, although the amount of loss associated with sediment transport was <10%. ERY had the highest concentrations in sediment and had a relative loss associated with sediment >50%. TYL also had >50% relative loss associated with sediment, and its pseudo-partitioning coefficient (P-PC) was very high. The tetracyclines (TC and CTC) had very low aqueous concentrations and had the lowest absolute losses. If agricultural runoff is proven to result in development of resistance genes or toxicity to aquatic organisms, then erosion control practices could be used to reduce TC, ERY, and TYL losses leaving agricultural fields. Other methods will be needed to reduce transport of other antibiotics.     

THE IMPACT OF RUNOFF GENERATION MECHANISMS ON THE LOCATION OF CRITICAL SOURCE AREAS1

ABSTRACT: Identifying phosphorus (P) source areas and transport pathways is a key step in decreasing P loading to natural water systems. This study compared the effects of two modeled runoff generation processes ‐ saturation excess and infiltration excess ‐ on total phosphorus (TP) and soluble reactive phosphorus (SRP) concentrations in 10 catchment streams of a Catskill mountain watershed in southeastern New York. The spatial distribution of runoff from forested land and agricultural land was generated for both runoff processes; results of both distributions were consistent with Soil Conservation Service‐Curve Number (SCS‐CN) theory. These spatial runoff distributions were then used to simulate stream concentrations of TP and SRP through a simple equation derived from an observed relation between P concentration and land use; empirical results indicate that TP and SRP concentrations increased with increasing percentage of agricultural land. Simulated TP and SRP stream concentrations predicted for the 10 catchments were strongly affected by the assumed runoff mechanism. The modeled TP and SRP concentrations produced by saturation excess distribution averaged 31 percent higher and 42 percent higher, respectively, than those produced by the infiltration excess distribution. Misrepresenting the primary runoff mechanism could not only produce erroneous concentrations, it could fail to correctly locate critical source areas for implementation of best management practices. Thus, identification of the primary runoff mechanism is critical in selection of appropriate models in the mitigation of nonpoint source pollution. Correct representation of runoff processes is also critical in the future development of biogeochemical transport models, especially those that address nutrient fluxes.     

Phosphorus species and fractionation--why sewage derived phosphorus is a problem

Phosphorus (P) inputs to sewage treatment works (STW) come from a variety of sources and filtration of treated wastewater prior to discharge into receiving waters is a common practice. This means P in treated wastewaters may be present in forms that are potentially more bioavailable and mobile. We conducted a 2-year study to determine P species up and downstream of two STW outfalls into two tributaries of the River Thames. Downstream of the outfalls, P concentrations in both rivers were frequently greater by an order of magnitude for all species of P. A high proportion of total P (TP) in the downstream waters was determined as dissolved, which was largely comprised of soluble reactive P (SRP) - considered as the most bioavailable P species. Furthermore no significant difference in SRP was found in receiving waters passed through 0.45 and 0.10 μm filters. This means that P from STWs occurs in <0.1 μm fraction size, which will not readily settle to the channel bed and is more easily assimilated by biota. This distinguishes STW inputs from agricultural runoff where a high proportion of P occurs as particulate P which is both less bioavailable and more likely to settle to the channel bed. This implies that STWs derived P is likely to have a greater adverse impact on the receiving river than agricultural runoff.     

Seasonal variation of herbicide concentrations in prairie farm dugouts

Prairie farm dugouts are frequently constructed for use as potable water sources. Consequently, cumulative pesticide inputs via atmospheric deposition and surface runoff may constitute a risk to human health. Since, relative to other pesticides, herbicides are used in greatest amount on the Canadian prairies, herbicide concentrations were intensively monitored in three dugouts over three growing seasons. Herbicides were detected in the water of all three dugouts each growing season which may reflect cumulative inputs from atmospheric and surface processes over the lifetimes of the dugouts, which varied from 11 to 22 yr. Detections, which were not continuous, tended to be seasonal in nature. During the 3-yr study, detections were most frequent during the spring application period and late fall following dugout turnover. Between these periods, herbicide concentrations generally decreased to below detection limits. The reappearance of herbicides in the dugout water during fall turnover and in concentrations generally greater than those present during the spring application period suggest that, under appropriate environmental conditions, the bottom sediments may act as a source of herbicides to the water column. In general, herbicide inputs due to deposition of application drift did not result in detectable concentrations of herbicides in the dugouts. In the only year that winter samples were monitored, herbicides were also detected during ice cover. On the basis of monthly sampling over each growing season, median concentrations of 9 of the 10 herbicides monitored were less than 0.05 microg L(-1). The exception, 2,4-D, which has been used extensively on the Canadian prairies for more than 50 yr and in greatest amounts, was the most frequently detected herbicide. In no case did herbicide concentrations exceed Canadian drinking water guidelines; however, on occasion maximum herbicide concentrations did exceed aquatic life and irrigation water guidelines.     

MODELING THE IMPACT OF MULTICOMPONENT VOCs ON GROUND WATER USING THE STEFAN‐MAXWELL EQUATION1

ABSTRACT: Computer programs that model the fate and transport of organic contaminants through porous media typically use Fick's first law to calculate vapor phase diffusion. Fick's first law, however, is limited to the case of a single, dilute species diffusing into a stagnant, high concentration, bulk vapor phase. When dealing with more than one diffusing species and at higher concentrations, the multicomponent coupling effects on vapor phase diffusion and advection of the various constituents become significant. VLEACH, a one‐dimensional finite difference model developed for the U.S. Environmental Protection Agency (USEPA), is typical of the models using Fick's first law to model vapor‐phase diffusion. The VLEACH model was modified to accommodate up to 10 components and to calculate the binary diffusion coefficients for each of the components based on molecular weight, molecular volume, temperature and pressure, and to address the coupling effects on multiple component vapor phase diffusion and its impact on ground water. The resulting model was renamed MC‐CHEMSOIL. At low vapor phase concentrations, MC‐CHEMSOIL predicts identical ground water impacts (dissolved phase loading) to those from VLEACH 2.2a. At higher vapor phase concentrations, however, the relative difference between the models exceeded 20 percent.     

A model for phosphorus transformation and runoff loss for surface-applied manures

Agricultural P transport in runoff is an environmental concern. An important source of P runoff is surface-applied, unincorporated manures, but computer models used to assess P transport do not adequately simulate P release and transport from surface manures. We developed a model to address this limitation. The model operates on a daily basis and simulates manure application to the soil surface, letting 60% of manure P infiltrate into soil if manure slurry with less than 15% solids is applied. The model divides manure P into four pools, water-extractable inorganic and organic P, and stable inorganic and organic P. The model simulates manure dry matter decomposition, and manure stable P transformation to water-extractable P. Manure dry matter and P are assimilated into soil to simulate bioturbation. Water-extractable P is leached from manure when it rains, and a portion of leached P can be transferred to surface runoff. Eighty percent of manure P leached into soil by rain remains in the top 2 cm, while 20% leaches deeper. This 2-cm soil layer contributes P to runoff via desorption. We used data from field studies in Texas, Pennsylvania, Georgia, and Arkansas to build and validate the model. Validation results show the model accurately predicted cumulative P loads in runoff, reflecting successful simulation of the dynamics of manure dry matter, manure and soil P pools, and storm-event runoff P concentrations. Predicted runoff P concentrations were significantly related to (r2=0.57) but slightly less than measured concentrations. Our model thus represents an important modification for field or watershed scale models that assess P loss from manured soils.