Balancing the weighting of various components of phosphorus loss in models is a critical but often overlooked step in accurate estimation of risk of P loss under field conditions. This study compared the P loss coefficients used to predict dissolved P losses from desorption from accumulated P in the soil, and those incidental to applications of P as fertilizer or manure, with extraction coefficients determined from actual P losses reported in literature for sites in Canada, with the addition of some sites with similar soils and climate from some northern states. The extraction coefficients for dissolved P measured in runoff water were 6.5× greater in year-round edge-of-field (EoF) measurements than in runoff boxes, indicating that models using P extraction coefficients derived from runoff box experiments will be underestimating the magnitude of losses from P accumulation in soil. Differences among the measurement methods (runoff box, rainfall simulator, or EoF) were not evident for incidental losses from applied P, but current models appear to overpredict the losses of applied P. Good agreement between measured and predicted dissolved P (DP) concentrations using the equations in the Annual Phosphorous Loss Estimator model were achieved by applying coefficients of 0.275 to the fertilizer equations and 0.219 to the manure equations, implying that 72.5% of fertilizer P and 78% of manure P are not available for runoff. This study underlines the importance of considering the relative weights of the various components of P loss as new models are developed and validated.
Introduction
Losses of both particulate and dissolved phosphorus (DP) from agricultural land have been identified as important drivers for surface water quality impairment by encouraging harmful (Steffen et al. 2014) and nuisance algae growth (Howell and Dove 2017) and hypoxic zones (Scavia et al. 2014). Recently, attention has focused on the disproportionate role of DP in algae blooms because of its high bioavailability and because the proportion of DP to total P exported from farmland appears to be increasing, particularly in the Lake Erie basin (Joosse and Baker 2011; Baker et al. 2014). This shift in the proportion of P fractions has been attributed to changes in agricultural management including a conversion to conservation tillage combined with surface broadcast of P fertilizers (Jarvie et al. 2017; Smith et al. 2017). DP losses from agricultural land are generally attributed to one of three sources: (i) desorption of a portion of the P that has accumulated on soil particles (from previous applications of mineral or organic fertilizers or native in the soil) into runoff water (desorbed P, DeP); (ii) incidental losses of DP from recent application of P fertilizer or manure while these remain on the surface of the soil (applied P, ApP) (Reid et al. 2018); or (iii) losses of P from the above-ground portion of overwintering vegetation (Liu et al. 2019). The third source has been identified as being important for prairie environments, but less so for parts of Canada with milder winters (Plach et al. 2019). This paper focuses on the P desorption from the soil and P loss incidental to fertilizer or manure application, as being important for the entire country. The third source (P losses from overwintering vegetation) has recently been thoroughly reviewed in Liu et al. (2019) and will not be discussed further in this paper.
Models are of critical importance in understanding the drivers for DP losses and for predicting the impacts of management practices on reducing DP exports (Bolster et al. 2012; Sharpley et al. 2012). There is a continuum of spatial and temporal scales for models of water and solute movement, but they can be classed very broadly into the bench-top or pedon scale mechanistic models, edge-of-field (EoF) models, and watershed-scale models. An example of a very detailed mechanistic model is Hydrus (Šimůnek et al. 2012), which requires many parameters at a high level of precision. It would be impractical to scale up to a field or regional scale and therefore, is out of scope for this study. At the opposite end of the scale are watershed models, such as SWAT (Radcliffe et al. 2009; Woodbury et al. 2014), SPARROW (Benoy et al. 2016), or CANWET (Singh et al. 2012). A common element in all of these models is the estimation of DP concentration and (or) loading in runoff water. Pferdmenges et al. (2020) have compiled a comprehensive review of the P source routines used in runoff models and point out that a large number of watershed models use the same conceptual model, initially developed for the EPIC model, for DP release into runoff and pore water. These quasi-mechanistic P routines assume partitioning between organic and mineral forms of P in classes of varying degrees of stability, with transfers between these pools controlled by zero- or first-order kinetics. Unfortunately, as Das et al. (2019) and Neumann et al. (2021) point out, the EPIC model (and, by extension, all of the models based on the EPIC P routine) does not accurately characterize P cycling across a wide range of soils, and better models of P release from soils are required. For a mechanistic model such as EPIC, this would require much greater detail regarding soil properties and initial soil P fractions, and at a much finer spatial scale, than is widely available. The P routines embedded in SWAT (Radcliffe et al. 2009) have also been shown to provide poor estimates of P losses from applications of concentrated P sources such as manure or fertilizer, in large part because these additions are assumed to immediately enter the soil pool with no allowance for any incidental losses from the applied materials (Collick et al. 2016; Menzies Pluer et al. 2019). The result is that most watershed models are much less sensitive to the impact of field-scale management decisions on dissolved P losses than on total P losses (that tend to be dominated by soil erosion and particulate P). This has been partially, but not completely, addressed in the Conservation Effects Assessment Program in the United States (CEAP 2011) by linking the Agricultural Policy Environmental Extender (APEX) Model (Williams and Izaurralde 2005; Gassman et al. 2009) to SWAT. This linkage allows better characterization of fertilizer and manure applications within the HRUs used in the SWAT model but does not correct the lack of accounting for incidental P losses within SWAT.
The third group of models focuses on P losses at the EoF and explicitly accounts for variability in soil characteristics and in soil, fertilizer, and manure management. The most common of these are the P Indexes (Sharpley et al. 2012; Sharpley et al. 2017), which have evolved over time from conceptual, additive models (Lemunyon and Gilbert 1993) to mixed empirical–mechanistic models with a component structure such as the N.C. P Loss Assessment Tool (N.C. PLAT Committee 2005). These indexes do not attempt to predict P losses from runoff events, but rather to indicate the relative risk of P loss from fields based on inherent soil and climate conditions combined with soil and nutrient management practices by farmers. Also included in the EoF group of models is the semiquantitative Annual Phosphorous Loss Estimator (APLE) (Vadas 2017), which has been proposed to be robust enough for use in validating other P loss models (Fiorellino et al. 2017).
Rather than a process-based model of DP release, the EoF models use empirical relationships between soil test P values (either an agronomic test or an environmental indicator such as water extractable P (WEP)) and DP concentration in runoff. Additional coefficients are used to predict the incidental losses of DP from applied fertilizer or manure. Discussion of the derivation of these coefficients follows later in the paper. This approach has been used in the national Indicator of Risk of Water Contamination by Phosphorus (IROWC-P) to predict the spatial and temporal variation in the risk of P loss (Reid et al. 2018). The coefficients for DP loss from accumulated P in soil have been updated from the previous versions of IROWC-P (van Bochove et al. 2010), and factors for the incidental losses from applied P have been added.
A key part of model development is validation of the model outputs against measured data to assess the accuracy and precision of the model predictions (Bolster and Vadas 2013; Vadas et al. 2013). One challenge with this validation step is avoiding being misled by equifinality, where models with different weightings of coefficients can give very similar results, so simply accepting the model with the most precise fit may not necessarily give the best characterization of the underlying processes (Hollaway et al. 2018). This could occur, for example, where the validation dataset includes fields with elevated soil test values that are also receiving applications of manure; models that attribute different relative weightings to the soil test P and the applied P could give similar results, and the validation data would not be able to ascertain which model was correct. This takes on great importance if the models are going to be used in different geographies or are going to be used to predict the impact of management changes on P losses. The authors are not aware, however, of any studies that have attempted to validate the relative weighting of the components in predictions of P losses within a comprehensive model, to ensure that not only are the model totals accurate but that the values of the individual components can be used confidently to identify appropriate management responses. This study aims to overcome that deficiency.
Methods
IROWC-P is a national-scale indicator of the spatial and temporal variability in the risk of P loss from Canadian agricultural land, derived from publicly available data on land use, nutrient application, and on-farm practices. The spatial scale of the assessment is the Soil Landscape of Canada (SLC) polygon. It uses a component structure, where each source of P loss from agricultural land (desorption of DP from soil; DP losses incidental to application of fertilizer or manure; DP leached from overwintering vegetation; and bioavailable particulate P from water erosion) is multiplied by a transport modifier appropriate for that P source, and then the products of these calculations are summed to give an overall risk of P loss at the edge of the field. The transport modifiers are derived from the surface runoff and tile flow predictions in the hydrology module of the DNDC model, which accounts for soil type, crop cover, and local precipitation. Delivery modifiers based on topography and distance to streams then predict the risk of P from agricultural land reaching surface water. The details of this modelling approach are found in Reid et al. (2018), and a very similar structure is used in the provincial P Loss Assessment Tool for Ontario (PLATO) (Ontario Ministry of Agriculture, Food and Rural Affairs 2021).
The goals of the recent updates to the IROWC-P model were to incorporate recent advances in our understanding of P losses from agricultural land while maintaining parsimony in data inputs and model structure. To accomplish this, updated coefficients were derived from current literature to predict DeP (Wang et al. 2010) and ApP for fertilizer (Vadas et al. 2008) and manure (Vadas et al. 2004; Vadas et al. 2009). The DeP coefficient is multiplied by soil test P (STP, converted to WEP as described below) and runoff (surface, plus the proportion of subsurface runoff that reaches tile drains through preferential pathways) to estimate the mass of DP that could be exported under average weather conditions (Reid et al. 2018). The DeP coefficient from Wang et al. (2010) was chosen as it represented a broad range of Ontario soil types and STP values and produced values in a similar range to those found by Vadas et al. (2005). Both studies showed strong correlations between DeP and soil WEP that appeared to be independent of soil type. The DeP coefficient should be validated in future for soils across Canada, to confirm that it is independent of soil type, but it represents the best currently available data.
The ApP coefficients were used in the manner of Vadas (2017), where a P distribution factor (PDF) is multiplied by the runoff: precipitation ratio between and an estimate of the WEP applied to estimate the mass of applied P that could be exported. The theoretical basis for this is that the runoff concentration of DP desorbed from freshly applied fertilizer or manure will be more concentrated than any runoff that has only interacted with the soil, so infiltration into the soil will trap the P from the runoff and reduce potential losses. The PDF values are different for fertilizer and manure, reflecting the differing nature of these two materials (Vadas et al. 2004; Vadas et al. 2008). For all materials, these calculations apply only to what remains on the soil surface, as fertilizer or manure that has been placed below the soil surface is assumed to be protected from desorption to runoff (Daverede et al. 2004). Differences in manure characteristics between livestock species are accounted for by basing the calculations on the water-extractable portion of the P in the manure rather than the total P applied.
These predictions of incidental P losses were integrated into larger models to allow their use for assessment of water quality impacts from agricultural activities. The ApP coefficients were developed using data from studies across the United States, as cited in Vadas et al. (2004, 2008, and 2009), and have been incorporated into the APLE model (Vadas 2017), which has been recognized as a suitable model for assessing the performance of P Indexes in the United States (Bolster et al. 2012; Bolster et al. 2014). Significant effort was expended in assessing the precision of each these coefficients during their development, so they were simply added together along with the soil desorption components when the APLE model was created. There have also been assessments of the accuracy of the APLE model that combines the predictions using these coefficients with other components to predict total P losses (e.g., Bolster and Vadas 2013; Kleinman et al. 2017; Sharpley et al. 2017).
A literature review was conducted to gather existing data on P runoff across a range of conditions. The geographic scope of this review was primarily Canada, but studies from some of the northern American States were included where, in the authors’ opinion, the soil and climatic conditions were similar enough to the adjacent Canadian regions to be relevant to IROWC-P (See Table 1 for a full list). Criteria for study selection included data for DP concentrations in runoff from individual treatments (or included both P loading and mean runoff volume so DP concentrations could be calculated), STP, soil texture, and P application data (source, rate, timing, incorporation method, and timing), as well as the type of runoff study. Some papers were excluded because the data was presented in aggregate to derive P loss coefficients so individual values for DP losses could not be determined (e.g., Wang et al. 2010; Vadas et al. 2005) or because only tile or surface runoff was reported but not both (e.g., Zhang et al. 2015). Runoff DP concentration data was categorized according to four methods of sample collection:
year-round EoF monitoring,
rainfall simulator studies in-situ on each plot at various times of the year, usually following a P application (RS),
measurements from individual rainfall events, and
runoff box studies from samples collected from the field and repacked into boxes for rainfall simulator studies in the laboratory (RB).
Table 1.
Summary of plots used for assessment of the IROWC-P and PLATO models.
![cjss-2021-0049tab1.gif](ContentImages/Journals/cjss/102/2/CJSS-2021-0049/graphic/WebImages/cjss-2021-0049tab1.gif)
Only one study reported measurements of DP concentrations from individual rainfall events, so it was excluded from further statistical analysis for DeP. Studies that included P applications were only considered if there was a check treatment without any added P, so that the runoff P incidental to applied manure or fertilizer could be separated from the inherent P losses from the soil. A total of 38 references had adequate data for inclusion in the assessment (Table 1).
The structure of a component P Index, like that used in IROWC-P and PLATO, calculates the contribution from DeP and ApP separately and then sums them. For consistency with this structure, the contribution from DeP was deducted from the total DP runoff by subtracting the P loading in runoff from check plots from the P loading in the corresponding amended plot, as follows:
where DPTreated is the mass of DP from the plots amended with fertilizer or manure, and DPCheck is the mass of DP from the corresponding unamended plot. This avoided the double counting of the contribution of P losses from soil desorption that would otherwise occur.
WEP concentrations were calculated from the STP concentrations, and where appropriate, soil textural class and soil pH, according to the method used in IROWC-P (details of this method are included in the Supplementary Material (cjss-2021-0049suppla.docx)). The coefficients used to convert STP to WEP varied from province to province because the standard soil tests vary between provinces (Mehlich-3 in Quebec and Atlantic Canada, Olsen in Ontario and Manitoba, and Modified Kelowna Saskatchewan, Alberta, and British Columbia); there are also cases of different measured conversion factors between provinces, e.g., Manitoba soils showed higher proportion of WEP in Olsen P than Ontario soils (Table 2 and Supplementary Fig. S1 (cjss-2021-0049suppla.docx) 1 1). Standardizing the soil P concentrations as WEP allowed comparison of a single extraction coefficient to estimate DeP in runoff water across the country. Where WEP was not reported directly in the American studies, values for STP were converted to WEP using the methodology for the nearest province. Predicted values for runoff P concentrations were calculated from the WEP data, and the conversion factors from WEP to DeP were used in the IROWC-P model. Measured and predicted values for runoff P losses from soil desorption and applied P were compared graphically, and goodness-of-fit statistics (R2, RMSE) were calculated using the Analysis Tool Pack add-in in Excel (Microsoft 2016).
Table 2.
Summary statistics of regression between measured DP concentrations in runoff from different study types (EoF, year-round edge-of-field; RS, rainfall simulated field plots; RB, runoff boxes) and soil WEP concentrations (no applied P).
![cjss-2021-0049tab2.gif](ContentImages/Journals/cjss/102/2/CJSS-2021-0049/graphic/WebImages/cjss-2021-0049tab2.gif)
Results and discussion
P loss from desorbed P
The results from the three different water quality sample collection methods (EoF, RS, and RB), where sufficient data was available in the literature, were compared to determine if they were predicting P losses in a similar manner, in terms of both total losses and relative to STP or applied P. When assessing losses of DeP, the measured values of DP in runoff appeared to fall into different populations with higher mean concentrations observed in EoF monitoring than in RB, although the differences were not statistically significant (Table 2). The range of STP in each of the measurement categories also varied, so that when the DP concentrations were regressed against soil test (expressed as WEP), the extraction coefficients, indicated by the slopes of the lines, were significantly different (Table 2 and Fig. 1). Year-round EoF measurements were felt by the authors to be more relevant to environmental impacts of P losses from agricultural land, as they represent the spatial scale most likely to influence adjacent water bodies and the temporal scale, which represents the range of runoff events over the year, so this was the primary focus of comparisons between measured and modelled DeP. The extraction coefficients in the literature cited for all of the measurement systems were much higher than those reported by Wang et al. (2010) or Vadas et al. (2005), so to maintain consistency with the original dataset used in the IROWC-P model, the relative difference between extraction ratios was used rather than applying the extraction coefficients directly from the literature values.
Fig. 1.
Comparison of observed and modelled DP concentrations in annual surface runoff from check plots (no added P) to soil WEP values from different measurement systems. Statistical values for these correlations are in Table 2.
![cjss-2021-0049f1.tif](ContentImages/Journals/cjss/102/2/CJSS-2021-0049/graphic/WebImages/cjss-2021-0049f1.jpg)
The impact of the choice of extraction coefficients in predicted DP concentrations is significant, with much lower predictions from the coefficients derived from runoff box data than from year-round EoF studies (Fig. 2).
Fig. 2.
Comparison of measured to modelled DP predictions from soil WEP using coefficients from year-round edge-of-field or runoff box experiments. Regression statistics are found in Table 2.
![cjss-2021-0049f2.tif](ContentImages/Journals/cjss/102/2/CJSS-2021-0049/graphic/WebImages/cjss-2021-0049f2.jpg)
The initial estimates of P desorption for use in IROWC-P were based on extraction coefficients derived from runoff box experiments by Wang et al. (2010), which had shown similar relationships to STP as Vadas et al. (2005). The study by Wang et al. (2010) showed a linear relationship between DP concentration in runoff and STP, with a coefficient of 0.0183 mg·L−1 DP in runoff water per mg·kg−1 of WEP. Applying this coefficient to WEP values reported in the literature showed a strong linear correlation between STP and predicted DP concentrations (Fig. 2), but the model was significantly underestimating the DP losses from soil desorption at EoF. The slope of the correlation was lower than the measured DP concentrations from these plots by a factor of 6.5 (the slope of the regression of modelled to observed EoF plots was greater than the slope of the RB plots by this factor; Fig. 1 and Table 2). Multiplying the original coefficient by the observed difference between the RB and the EoF values (6.5) increased the coefficient to a value of 0.114 and brought the predicted and measured DP concentrations into nearly perfect alignment. When the revised coefficients are used, there is a very close correspondence between the measured and modelled DP values from DeP (Fig. 3), with a slope that does not differ significantly from one and an R2 value of 75%.
Fig. 3.
Comparison of modelled to observed DP concentrations in surface runoff for the edge-of-field plots with no added fertilizer or manure.
![cjss-2021-0049f3.tif](ContentImages/Journals/cjss/102/2/CJSS-2021-0049/graphic/WebImages/cjss-2021-0049f3.jpg)
The reasons for the discrepancy between year-round EoF measurements and measurements based on artificial rainfall are unclear. The studies for all three measurement systems included a broad range of soil types, covering similar ranges of soil texture and pH and with no obvious differences in past management where information was provided. The runoff box studies included soils with higher values for WEP than the EoF or RS, but there is no evidence that the slope of the correlation with runoff P concentrations changes at the higher WEP values. Kleinman et al. (2004) and Vadas et al. (2005) found only small differences between DP concentrations measured in runoff box and field plots using a rainfall simulator (with a trend toward higher concentrations from the rainfall simulator studies), but neither of these studies included year-round measurement of DP concentration in runoff from natural precipitation. In contrast, Doody et al. (2006) found that DP concentrations in runoff water increased as the length of flow path increased, which could partly explain the pattern of higher concentrations measured at the EoF compared with artificial runoff experiments with small areas. Little et al. (2007) also noted that P extraction coefficients for small Alberta watersheds were higher than those measured in runoff box experiments. Studies from Ontario (Lozier et al. 2017) and the Canadian prairie provinces (Tiessen et al. 2010) have suggested that longer contact time between runoff water and soil, as could occur in snowmelt, would increase soluble reactions between soil and water. The implication of this observation is that the direct application of runoff coefficients from runoff boxes or rainfall simulator studies to models of DP loss from agricultural fields can result in a significant underestimation of the actual risk of P loss.
P loss from applied manure and fertilizer
When losses incidental to the application of P fertilizer or manure were compared among the three sample collection methods, each appeared to fall into distinct populations for both the amount of applied P and the measured ApP losses, in order of EoF < RS < RB (Table 3). The differences among the measurement methods in ApP, however, were not due to the measurement method but rather to the rate of P applied. When the ratio of P loss from ApP to mass of P applied as fertilizer or manure was calculated, there were no significant differences among the three sample collection methods for either fertilizer or manure, indicating that the proportion of applied P that was desorbed into runoff water was the same for all three methods. Since application rate is part of the prediction of ApP, the data from the sample collection methods was pooled for assessment of the model.
Table 3.
Mean values for DP concentration in runoff incidental to P application (ApP), rates of applied P as fertilizer or manure, and the DP concentration as a proportion of applied P.
![cjss-2021-0049tab3.gif](ContentImages/Journals/cjss/102/2/CJSS-2021-0049/graphic/WebImages/cjss-2021-0049tab3.gif)
Applying the coefficients derived by Vadas et al. (2009) to predict the proportion of applied P from both manure and fertilizer that was lost in runoff resulted in values that were significantly higher than those that were measured by any of the sampling methods. It should be noted that the derivation of the PDF values did not appear to segregate the DeP from ApP, as has been done here, which would result in higher apparent P loss from ApP. This effect would be largest in cases where the P applications were on soils that were already high in WEP (e.g., from historic manure applications in excess of crop P requirements). The PDFs from the APLE model (Vadas et al. 2009) have unique values for fertilizer and manure to account for the differences in the P release characteristics of the materials. These equations predicted that close to 100% of the available WEP would be lost in runoff at high runoff: precipitation ratios, but across all of the studies reported in the literature, WEP losses were less than one-third of this value.
The relative amount of P loss for both fertilizer and manure as the runoff: precipitation ratio increased appeared to be well characterized by the PDF equations, as indicated by statistically significant correlations between measured and predicted DP concentrations in runoff. Using the equations directly from the APLE model, however, resulted in significant overprediction of DP concentrations from applied P in runoff water (Figs. 4 and 5), by a factor of 3.63× for fertilizer and 4.57× for manure (Note: the assumption from the APLE model that 60% of liquid manure infiltrates immediately was not used in this analysis, so that solid manure was not assigned a 2.5× greater risk of P loss than liquid manure; this discrepancy in risk was not apparent in the data we examined). The relative amount of applied P that was lost in runoff water was very well predicted by the APLE equations, but the predicted results only matched the measured when the equations were multiplied by the inverse of the overprediction, resulting in factors of 0.275 for fertilizer and 0.219 for manure, implying that 72.5% of fertilizer P and 78% of manure P are not available for runoff. (Figs. 4 and 5). Regression analysis shows that the slope of the measured to modelled correlation was not different from 1. The intercept for the fertilizer prediction was not significantly different from 0, while for manure there was a significant positive intercept. This suggests that the model is predicting some risk of P loss inherent in manure application, although this value is small relative to the total risk. There was also greater variability in the measured P losses incidental to manure application. This is not unexpected given the inherently greater variability in manure P stemming from uncertainty regarding the precise proportion of both DP and total P in manure (Kleinman et al. 2005b) and variability in the amount remaining on the soil surface following incorporation due to the bulky nature of manure. Relations between measured and predicted DP concentrations for both fertilizer and manure were highly significant (p<0.001) with moderately high R2 values and low RMSE (Figs. 4 and 5). Residuals were evenly distributed with average = 0 and no significant slope. The adjustment to the fertilizer availability for loss is consistent with the original model proposed by Vadas et al. (2007), which included a factor for P absorption by the soil between the time of fertilizer application and the first runoff event. In the literature investigated, it appears that the greatest fertilizer loss expected would be 27.5% of the total applied.
Fig. 4.
Modelled DP concentrations from fertilizer applications compared with measured DP concentrations. Open circles and dashed line are predicted P losses incidental to fertilizer application using the equations in the APLE model, which overpredicted the P losses by 3.63×. The closed circles and solid line show the result of multiplying the APLE model equations by a factor of 0.275. Statistics are for the linear fit of the revised model.
![cjss-2021-0049f4.tif](ContentImages/Journals/cjss/102/2/CJSS-2021-0049/graphic/WebImages/cjss-2021-0049f4.jpg)
Fig. 5.
Modelled DP concentrations from manure applications compared with measured DP concentrations. Open circles and dashed line are predicted P losses incidental to manure application using the equations in the APLE model, which overpredicted the P losses by 4.57×. The closed circles and solid line show the result of multiplying the APLE model equations by a factor of 0.219. Statistics are for the linear fit with the revised model.
![cjss-2021-0049f5.tif](ContentImages/Journals/cjss/102/2/CJSS-2021-0049/graphic/WebImages/cjss-2021-0049f5.jpg)
The APLE model includes a factor for infiltration of liquid manure into the soil, assuming that 60% of the applied liquid infiltrated and therefore protected the P from loss in runoff (Vadas et al. 2004). This had not been included in the original IROWC-P model, to treat the P from solid manure in a similar manner to liquid manure as a source of P runoff. If infiltration was the only mechanism reducing the loss of P from applied manure, we would expect to see relatively greater losses from solid manure applications than from liquid, but this difference was not observed in the literature results. It suggests there is another mechanism in place. Possibly the temporal pattern of runoff generation, which is accounted for in the SWCS Curve Number approach as the “initial abstraction” (SCS 1985; Rawat et al. 2010), is carrying a disproportionate amount of desorbable P into the soil during the early part of a rain event and therefore, limiting the amount available for runoff. There could also be different mechanisms operating with liquid and solid manure (infiltration of liquid manure and sequestration of P within solid manure either in organic forms or physically protected from diffusion out of clumps) that result in similar apparent availability for loss from the two forms. This is an area that requires further investigation to fully understand the underlying processes of manure and fertilizer P dissolution and transfer to runoff water.
In the APLE model, P runoff from fertilizer or manure is estimated through the use of a PDF based on the runoff:precipitation ratio to partition DP between runoff and infiltration, assuming that the P that infiltrates into the soil is then adsorbed to the soil (Vadas et al. 2009). To align the predictions of P loss with the measured P concentrations in the literature, we suggest that the equations for the PDF for incidental losses from applied P need to be multiplied by the upper limit of P available for loss: 0.275 in the case of fertilizer and 0.219 in the case of manure.
Impact of model revisions on relative weighting of P loss components
A set of synthetic values was constructed using various combinations of soil erosion (Bio Available Particulate P, BAPP), soil test values, and manure and fertilizer application rates to test the relative contributions of each component using the original and revised coefficients (Fig. 6). The relative importance of STP to the overall scores increased, while manure and fertilizer were reduced, providing a much more balanced weighting for these components. It would be beneficial to carry out a sensitivity analysis to determine the impact of changes to individual components to the final P loss assessments, as well as to compare the overall scores to a range of water quality measurements from across the country to test the overall performance of the model, as these steps have not yet been completed.
Future work
Despite the broad range of conditions represented by the literature reviewed, further work is needed to increase confidence in the weightings of components of P loss. In particular, there needs to be emphasis on year-round EoF monitoring under a range of soil types, climatic conditions, P sources, and application rates. Performing these trials with a broad enough range of P management practices will increase confidence that the P losses are being partitioned correctly among the different components.
Conclusions
This paper compares the results of P runoff experiments relevant to Canadian soil and climatic conditions with modelled P concentrations, focusing in particular on the relative weights of P losses through desorption from the soil and incidental to P applications. The measurements of DeP varied significantly depending on what sampling methods were used. Within each sampling method, there were strong and consistent correlations between STP and DeP, but the DP concentrations measured in runoff boxes were significantly lower than in year-round EoF measurements. Models using DeP coefficients derived from runoff box studies are likely to be underestimating the contribution of elevated STP to risk of P loss from fields over the course of a year. Runoff box and rainfall simulator studies will continue to be valuable for understanding the underlying P release processes and the impacts of soil test P, soil texture and pH, and other soil parameters on the relative losses of DP in runoff. The controlled conditions of these studies are not reflective of EoF situations, and so scaling up their results to field scale models should be approached with caution.
Applying the same type of analysis to ApP showed much less variation between sampling methods, once the rate of P application had been considered. The amount of P desorption, however, did not approach that which would be expected if all of the applied P (or applied WEP in the case of manure) on the soil surface desorbed into runoff water. Models using the assumption of 100% availability for loss of P applied as manure or fertilizer appear to be overestimating the contribution of applied P to the risk of P loss within a component P index.
Ensuring correct balance among the components of P loss in any P index or P loss model will improve the overall accuracy of risk assessments and provide better guidance to land managers regarding the optimum choice of mitigation techniques if the risk of P loss is high. Further, it will improve the capacity of models to correctly predict the impact of applying various beneficial management practices on P loss reduction in both the short and long term, which is critical for meeting P reduction targets.
Acknowledgements
This project was funded as part of the Sustainability Metrics initiative of Agriculture and Agri Food Canada. Thanks to co-op student Alexander Nauta, who summarized the data from the literature sources and performed the initial data analysis. Thanks as well to the scientists who shared their data on the relations between STP and WEP (Wayne Pettipiece, Alvin Anderson, Norma Sweetland, Bob Eilers, Steve and Maria Sheppard, Wole Akinremi, Craig Drury, Tiequan Zhang, Gordon Barnett, Nicolas Tremblay, Denise Desrosiers, Pierre Audesse, Ken Webb, Delmar Holmstrom, as well as the Laboratoire de pédologie et d’agriculture de précision (Sainte-Foy, QC, Canada) and Tiequan Zhang and Yutao Wang for their assistance in data analysis.