Watershed-level comparison of predictability and sensitivity of two phosphorus models

Modeling of phosphorus loss from field to watershed: A review

Phosphorus (P) losses from nonpoint sources into surface water resources through surface runoff and tile drainage play a significant role in eutrophication. Accordingly, the number of studies involving the modeling of agricultural P losses, the uncertainties of such models, and the best management practices (BMPs) supported by the modeling of hypothetical P loss reduction scenarios has increased significantly around the world. Many improvements have been made to these models: separate manure P pools, variable source areas allowing the determination of critical source areas of P loss, analyses of modeling uncertainties, and understanding of legacy P. However, several elements are still missing or have yet to be sufficiently addressed: the incorporation of preferential flow into models, the modification of P sorption-desorption processes considering recent research data (e.g., pedotransfer functions for labile, active, or stable P, along with P sorption coefficients), BMP parameterization, and scale-up issues, as well as stakeholder-scientist and experimentalist-modeler interactions. The accuracy of P loss modeling can be improved by (a) incorporating dynamic P sorption-desorption processes and new P subroutines for direct P loss from manure, fertilizer, and dung, (b) modeling preferential flow, connectivity between field and adjacent water bodies, and P in-stream processes, (c) including an assessment of model uncertainty, (d) integrating field and watershed models for BMP calibration and scaling field results up to larger areas, and (e) building a holistic interaction between stakeholders, experimentalists, and modelers.

Modeling the Impacts of Manure on Phosphorus Loss in Surface Runoff and Subsurface Drainage

Simulation of phosphorus (P) transfer from manured agricultural lands to water bodies via surface runoff and subsurface drainage is potentially of great help in evaluating the risks and effects of eutrophication under a range of best management practice scenarios. However, it remains a challenge since few models are capable of providing a reasonably accurate prediction of P losses under manure treatment. The Environmental Policy Integrated Climate (EPIC) model was applied to simulate the impacts on dissolved reactive P (DRP) losses through surface runoff and subsurface drainage from a solid cattle manure-amended corn ( L.)-soybean [ (L.) Merr.] rotation on a clay loam soil (Vertisol) located in the Lake Erie region. Simulations of DRP loss in surface runoff and tile drainage were satisfactory; however, EPIC did not consider DRP loss directly from manure, weakening its accuracy in the prediction of DRP loss in surface runoff. Having previously drawn on EPIC-predicted surface runoff to initiate SurPhos (Surface Phosphorus and Runoff Model) predictions of DRP losses strictly in surface runoff, no comparison had been made of differences in manure application impacts on EPIC- or SurPhos-predicted DRP losses-accordingly, this was assessed. The SurPhos improved the estimation of DRP loss in surface runoff (Nash-Sutcliffe coefficient, 0.53), especially when large rain events occurred immediately after or within 6 wk of manure application. Generally, EPIC can capture the impacts of manure application on DRP loss in surface runoff and subsurface drainage; however, coupling of the EPIC and SurPhos models increased the accuracy of simulation of runoff DRP losses.

Modeling phosphorus losses from soils amended with cattle manures and chemical fertilizers

While applied manure/fertilizer is an important source of P loss in surface runoff, few models simulate the direct transfer of phosphorus (P) from soil-surface-applied manure/fertilizer to surface runoff. The SurPhos model was tested with 2008-2010 growing season daily surface runoff data from clay loam experimental plots subject to different manure/fertilizer applications. Model performance was evaluated on the basis of the coefficient of determination (R²), Nash-Sutcliffe efficiency (NSE), percent bias (PBIAS), and the ratio of the root mean square error to the standard deviation of observed values (RSR). The model offered an acceptable performance in simulating soil labile P dynamics (R² = 0.75, NSE = 0.55, PBIAS = 10.43%, and RSR = 0.67) and dissolved reactive P (DRP) loss in surface runoff (R² ≥ 0.74 and NSE ≥ 0.69) for both solid and liquid cattle manure, as well as inorganic fertilizer. Simulated direct P loss in surface runoff from solid and liquid cattle manure accounted for 39% and 40% of total growing season DRP losses in surface runoff. To compensate for the unavailability of daily surface runoff observations under snow melt condition, the whole four years' (2008-2011) daily surface runoff predicted by EPIC (Environmental Policy Integrated Climate) was used as SurPhos input. The accuracy of simulated DRP loss in surface runoff under the different manure/fertilizer treatments was acceptable (R² ≥ 0.55 and NSE ≥ 0.50). For the solid cattle manure treatment, of all annual DRP losses, 19% were derived directly from the manure. Beyond offering a reliable prediction of manure/fertilizer P loss in surface runoff, SurPhos quantified different sources of DRP loss and dynamic labile P in soil, allowing a better critical assessment of different P management measures' effectiveness in mitigating DRP losses.

Comparing an Annual and a Daily Time-Step Model for Predicting Field-Scale Phosphorus Loss

A wide range of mathematical models are available for predicting phosphorus (P) losses from agricultural fields, ranging from simple, empirically based annual time-step models to more complex, process-based daily time-step models. In this study, we compare field-scale P-loss predictions between the Annual P Loss Estimator (APLE), an empirically based annual time-step model, and the Texas Best Management Practice Evaluation Tool (TBET), a process-based daily time-step model based on the Soil and Water Assessment Tool. We first compared predictions of field-scale P loss from both models using field and land management data collected from 11 research sites throughout the southern United States. We then compared predictions of P loss from both models with measured P-loss data from these sites. We observed a strong and statistically significant ( < 0.001) correlation in both dissolved (ρ = 0.92) and particulate (ρ = 0.87) P loss between the two models; however, APLE predicted, on average, 44% greater dissolved P loss, whereas TBET predicted, on average, 105% greater particulate P loss for the conditions simulated in our study. When we compared model predictions with measured P-loss data, neither model consistently outperformed the other, indicating that more complex models do not necessarily produce better predictions of field-scale P loss. Our results also highlight limitations with both models and the need for continued efforts to improve their accuracy.

Evaluation of the APEX Model to Simulate Runoff Quality from Agricultural Fields in the Southern Region of the United States

The Agricultural Policy Environmental eXtender (APEX) model has been widely applied to assess phosphorus (P) loss in runoff water and has been proposed as a model to support practical decisions regarding agricultural P management, as well as a model to evaluate tools such as the P Index. The aim of this study is to evaluate the performance of APEX to simulate P losses from agricultural systems to determine its potential use for refinement or replacement of the P Index in the southern region of the United States. Uncalibrated and calibrated APEX model predictions were compared against measured water quality data from row crop fields in North Carolina and Mississippi and pasture fields in Arkansas and Georgia. Calibrated models satisfactorily predicted event-based surface runoff volumes at all sites (Nash-Sutcliffe efficiency [NSE] > 0.47, |percent bias [PBIAS]| < 34) except Arkansas (NSE < 0.11, |PBIAS| < 50) but did not satisfactory simulate sediment, dissolved P, or total P losses in runoff water. The APEX model tended to underestimate dissolved and total P losses from fields where manure was surface applied. The model also overestimated sediments and total P loads during irrigation events. We conclude that the capability of APEX to predict sediment and P losses is limited, and consequently so is the potential for using APEX to make P management recommendations to improve P Indices in the southern United States.

Improved Simulation of Edaphic and Manure Phosphorus Loss in SWAT

Watershed models such as the Soil Water Assessment Tool (SWAT) and the Agricultural Policy Environmental EXtender (APEX) are widely used to assess the fate and transport of agricultural nutrient management practices on soluble and particulate phosphorus (P) loss in runoff. Soil P-cycling routines used in SWAT2012 revision 586, however, do not simulate the short-term effects of applying a concentrated source of soluble P, such as manure, to the soil surface where it is most vulnerable to runoff. We added a new set of soil P routines to SWAT2012 revision 586 to simulate surface-applied manure at field and subwatershed scales within Mahantango Creek watershed in south-central Pennsylvania. We corroborated the new P routines and standard P routines in two versions of SWAT (conventional SWAT, and a topographically driven variation called TopoSWAT) for a total of four modeling "treatments". All modeling treatments included 5 yr of measured data under field-specific, historical management information. Short-term "wash off" processes resulting from precipitation immediately following surface application of manures were captured with the new P routine whereas the standard routines resulted in losses regardless of manure application. The new routines improved sensitivity to key factors in nutrient management (i.e., timing, rate, method, and form of P application). Only the new P routines indicated decreases in soluble P losses for dairy manure applications at 1, 5, and 10 d before a storm event. The new P routines also resulted in more variable P losses when applying manure versus commercial fertilizer and represented increases in total P losses, as compared with standard P routines, with rate increases in dairy manure application (56,000 to 84,000 L ha). The new P routines exhibited greater than 50% variation among proportions of organic, particulate, and soluble P corresponding to spreading method. In contrast, proportions of P forms under the standard P routines varied less than 20%. Results suggest similar revisions to other agroecosystem watershed models would be appropriate.