JEQ Grow Your Career With ASA
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Published online 31 August 2007
Published in J Environ Qual 36:1452-1460 (2007)
DOI: 10.2134/jeq2006.0207
© 2007 American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Watson, C. J.
Right arrow Articles by Matthews, D. I.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Watson, C. J.
Right arrow Articles by Matthews, D. I.
Agricola
Right arrow Articles by Watson, C. J.
Right arrow Articles by Matthews, D. I.
Related Collections
Right arrow Surface Water Quality
Right arrow Phosphorus
Right arrow Field-Scale Studies

ENVIRONMENTAL ISSUES

Increase in Phosphorus Losses from Grassland in Response to Olsen-P Accumulation

C. J. Watsona,*, R. V. Smitha and D. I. Matthewsb

a Agriculture, Food and Environmental Science Division
b Biometrics Branch, Agri-Food and Biosciences Institute, Newforge Lane, Belfast BT9 5PX, United Kingdom

* Corresponding author (Catherine.Watson{at}afbini.gov.uk).

Received for publication May 26, 2006.
ABSTRACT

The Olsen-P status of grazed grassland (Lolium perenne L.) swards in Northern Ireland was increased over a 5-yr period (March 2000 to February 2005) by applying different rates of P fertilizer (0, 10, 20, 40, or 80 kg P ha–1 yr–1) to assess the relationship between soil P status and P losses in land drainage water and overland flow. Plots (0.2 ha) were hydrologically isolated and artificially drained to v-notch weirs, with flow proportional monitoring of drainage water and overland flow. Annually, the collectors for overland flow intercepted between 11 and 35% of the surplus rainfall. Single flow events accounted for up to 52% of the annual dissolved reactive phosphorus (DRP) load. The Olsen-P status of the soil influenced DRP and total phosphorus (TP) concentrations in land drainage water and overland flow. Annual TP loss was highly variable and ranged from 0.19 to 1.55 kg P ha–1 yr–1 for the plot receiving no P fertilizer and from 0.35 to 2.94 kg P ha–1 yr–1 for the plot receiving 80 kg P ha–1 yr–1. Despite the Olsen-P status in the soils ranging from 22 to 99 mg P kg–1, after 5 yr of fertilizer P applications it was difficult to identify a clear Olsen-P concentration at which P losses increased. Any relationship was confounded by annual variability of hydrologic events and flows and by hydrologic differences between plots. Withholding P fertilizer for over 5 yr was not long enough to lower P losses or to have an adverse effect on herbage P concentrations.

Abbreviations: AFWM, annual flow-weighted mean • AE, actual evaporation • DRP, dissolved reactive phosphorus • DUP, dissolved unreactive phosphorus • PP, particulate phosphorus • TDP, total dissolved phosphorus • TP, total phosphorus

DIFFUSE losses of phosphorus from agricultural land to surface water have been increasing in recent years, leading to increased eutrophication of many rivers, lakes, and impounding reservoirs worldwide (Foy and Bailey-Watts, 1998; Sharpley and Rekolainen, 1997). Evidence suggests that the increased loss of agricultural P is due to farms operating an annual P surplus of approximately 15 to 20 kg P ha–1 yr–1, and, as a result, their soils have become enriched with P (Edwards and Withers, 1998). In Northern Ireland, Foy et al. (1995) estimated that the rate of accumulation in soils in the Lough Neagh catchment over the last 50 yr was 10 kg total P ha–1 yr–1. Monitoring the six major rivers in this predominantly grassland catchment showed that the rates at which the dissolved reactive phosphorus (DRP) loadings were increasing were in the range of 13.0 to 14.5 g P ha–1 yr–1 over a 17-yr period (Foy et al., 1995). This increase is much greater than could be accounted for by changes in rural human population densities or industrial sources and was sufficient to counterbalance the benefits of a P removal program operating at sewage treatment works. Similar increases in P loads were found over a 28-yr period in Scandinavia by Krug (1993) and over a 15-yr period in Finland by Rekolainen (1989). These increases were postulated to be due to increases in surface and subsurface losses of P resulting from P accumulation in soils originating from an imbalance between agricultural inputs and outputs. It is estimated that the agricultural contribution to total P losses to inland waters in Northern Ireland is 58% and that the annual economic costs of eutrophication are in the range of 36 to £61 million (Smith et al., 2005). The target of the European Union's Water Framework Directive (European Commission, 2000) is to restore all surface water to good ecologic and chemical status by 2015.

Evidence that increases in DRP concentrations in land drainage water were a direct response to soil P accumulation has been shown by Smith et al. (2003) in a temporal study on a 9-ha grassland catchment at Greenmount Agricultural and Horticultural College, County Antrim, Northern Ireland, that had no domestic or industrial sources of P and no farm buildings. They showed that for the period 1989 to 1997, annual median DRP concentrations in land drainage water increased by 2.4 µg P L–1 yr–1 in response to a mean P surplus of 23.4 kg P ha–1 yr–1.

The question arises as to what soil P status is optimum to minimize environmental losses of P to surface waters while maintaining grass yields. Several dose-response curves have been produced where agricultural P losses are plotted against their soil P status. A breakpoint for DRP losses to drainage water of around 60 mg Olsen-P kg–1 was shown by Heckrath et al. (1995) from soils of the Broadbalk continuous wheat experiment at Rothamsted. However, Hesketh and Brookes (2000) showed that the breakpoint measured after laboratory extraction of different soils with CaCl2 ranged from 10 to 119 mg Olsen-P kg–1. They were unable to find any causal relationships between differences in soil chemical and physical properties or management and the different breakpoints. Knowledge of the relationship between P losses to water and soil P status, together with quantifying the contribution of overland flow to total P loss, is fundamental to developing a sustainable P management policy for agriculture. The aim of the current study was to increase the Olsen-P status of grazed grassland swards over a 5-yr period by applying P fertilizer and to assess the relationship between soil P status and P losses in overland flow and land drainage water. Particular emphasis is given to the DRP fraction because of its high bioavailability (Watson and Foy, 2001).

Materials and Methods

Study Site
The grassland experimental site was located at the Agricultural Research Institute of Northern Ireland at Hillsborough, County Down on a slightly gleyed, sandy clay-loam soil (48% sand, 31% silt, 21% clay, and 12% organic matter) overlying Silurian shale. Five plots, each 143 m x 14 m (0.2 ha), were laid side by side on the northerly aspect of a sloping drumlin. The plots were hydrologically isolated and artificially drained to separate v-notch weirs, which monitored flows in the range of 0.01 to 7.0 L s–1 to an accuracy of ±4%. Details of the installation of the drains have been described elsewhere (Watson et al., 2000). To prevent cross-plot flows of surface water, residual soil-spoil was placed on both sides of the drain trench. Interplot fences were erected on these spoil lines. The drains were installed during the spring and summer of 1987, and the pasture was plowed and reseeded to perennial ryegrass in August 1987. From March 1989 to February 1999, the plots received different rates of N fertilizer (0–500 kg N ha–1 yr–1) as ammonium nitrate–calcium carbonate (Watson et al., 2000) and were continuously grazed by beef steers from April to October of each year. The same, single annual application of P (as single superphosphate) and K (as muriate of potash) was applied to all plots, according to soil analysis and standard recommendations for grazing management (Ministry of Agriculture, Fisheries, and Food, 1994).

Ground limestone (54% neutralizing value) was surface applied to all plots in April 1999 at a rate of 5 t ha–1 to increase soil pH from 5.8 to 6.2. All plots received 250 kg N ha–1, 19 kg K ha–1, and 18 kg S ha–1 during 1999 and were under the same grazing management as in previous years. The N fertilizer was applied in six equal applications from March to August, the K fertilizer was applied as a single dressing in June, and the S fertilizer was applied with the N source on three occasions. No P fertilizer was applied during 1999. Because all plots received the same fertilizer applications during 1999, this acted as a break year from the previous 10-yr study, using different N fertilizer rates. There was no evidence that previous N application rates had any effect on P losses to drainage water (Watson, unpublished). The experimental P treatments were imposed from March 2000 to February 2005. The plots received 0, 10, 20, 40, or 80 kg P ha–1 as triple superphosphate (46% P2O5) in six equal applications between March and September each year. All plots received the same inputs of N (250 kg N ha–1 yr–1 applied with the P fertilizer), K (20 kg K ha–1 yr–1), and S (18 kg S ha–1 yr–1) and were continuously grazed by beef steers (7–12 mo old) from April to October to maintain a constant sward height of 7 cm. Output was estimated from live-weight gain of animals (measured every 20 d) allocated to each plot randomly. For most of the grazing season, the same two cattle grazed each plot, but the stocking rate was varied depending on growing conditions. The average number of grazing days for the 5-yr study period was 1718 d ha–1 yr–1 (or 344 d yr–1 per plot), and the average animal live-weight gain was 1226 kg ha–1 yr–1.

Sampling and Analysis
Drainage flow through each weir was monitored continuously and logged at 15-min intervals using Warren Jones WJ 460 FM flow monitors (Warren Jones Engineering, Bicester, Oxon, UK) fitted with ultrasonic probes capable of measuring water depths in the weirs to ±1 mm. Warren Jones Spartan-2 Magnum refrigerated autosamplers collected a daily composite water sample from each plot, comprising four 100-mL samples taken at 6-h intervals. The daily water samples were bulked to give a weekly composite on which various P fractions were analyzed.

A collector for overland flow was installed in a shallow trench across the width of each plot at the lowest end and consisted of a 185-mm diameter solid PVC pipe cut in half with a 45-mm flange attached to the top edge of one side. This was pushed horizontally into the soil to an approximate depth of 50 mm. The collector was connected to a mini v-notch weir (152 mm) and portable 6700 series ISCO water sampler (ISCO, Inc., Lincoln, NE) with a bubbler module. Flow data were recorded using ISCO's Flowlink software and downloaded using a 581 Rapid Transfer Device. Each collector for overland flow was fitted with an easily removable gray PVC lid (200 mm wide), which, when closed, rested on top of the grass stubble on the collection side. This prevented contamination of the pipe with bird droppings or debris but allowed overland flow (including shallow interflow to 50 mm) to enter the pipe. When overland flow occurred, 24 x 200 mL water samples were taken at 20-min intervals, and a composite sample was produced for the event. Phosphorus loads were calculated from recorded flows and measured P concentrations for each event. The overland flow collectors were protected from trampling by cattle with an electrified fence placed across the width of the plot approximately 0.3 m up from the collector pipe. The vegetation immediately around the collectors and water samplers was cut on a regular basis to a height of approximately 50 mm and removed.

Dissolved reactive P (DRP) was determined on a 0.45-µm, membrane-filtered (Millipore, UK) water samples by the acidic molybdate-ascorbic acid method of Murphy and Riley (1962). Total dissolved P (TDP) and total P (TP) were determined on filtered and unfiltered samples, respectively, by digestion with potassium persulphate and sulfuric acid, followed by analysis of the digest as for DRP (Eisenreich et al., 1975). Dissolved unreactive phosphorus (DUP = TDP – DRP) and particulate phosphorus (PP = TP – TDP) fractions were calculated. Annual flow-weighted mean (AFWM) concentrations were calculated for each P species by Eq. [1]:

Formula 1[1]
where w is week number (1–52), fw is the weekly flow through the weir, and cw is the concentration of the determinant in the weekly composite sample collected from each weir. Each experimental year ran from 1 March to 28 or 29 February of the following year.

Daily rainfall and actual evaporation (AE) data were obtained from the Central Climate Unit, Meteorological office, Exeter, Devon, UK using the Meteorological Office Rainfall and Evaporation Calculation System (MORECS) for Hillsborough. Daily rainfall was manually recorded at the site using a tipping bucket collector and was not significantly different from the MORECS data. Daily rainfall and AE values were summed to give monthly and annual total values. The efficiency with which the drains intercepted flow at the start of the 10-yr experimental period varied considerably (range, 72–103%) (Watson et al., 2000). For this reason, annual loads of P lost to drainage water were determined using the theoretical drainage volume calculated from meteorologic data (i.e., AFWM concentration was multiplied by annual rainfall minus AE). It was assumed that the volume intercepted on each plot was a representative sample of soil drainage water.

For soil sampling, each plot was divided into three sections (top, middle, and bottom), running down the slope. On a weekly basis, from January 2000 to February 2005, 30 cores (15 mm diameter x 75 mm depth) were randomly taken from each section and bulked. Soil samples were air dried at 30°C, ground to pass a 2-mm sieve, and analyzed for bicarbonate-extractable inorganic P (Olsen-P) according to the method of Olsen et al. (1954) using a soil to Olsen reagent ratio of 1:20. A linear model (Eq. [2]) was fitted to the weekly Olsen-P data over the 5 yr of the study:

Formula 2[2]
where Y is the estimated P status in mg kg–1, Pfert is the rate of P applied as kg P ha–1 yr–1, and Date is the date of the observation in Julian format. To allow for random variation in the observed values, the annual increase was derived by calculating the estimated Olsen-P value at the start and the end of the fertilizer year and subtracting the two values.

Grass clippings (>5 cm) were taken at approximately 2- to 3-wk intervals, depending on growing conditions, from March/April to September/October each year. The grass samples were dried at 102°C for 18 h and milled to pass through a 1-mm mesh. Samples were analyzed for N by dry combustion using a Carlo Erba NA1500 elemental analyzer and for P, K, Ca, Mg, and S by inductively coupled plasma spectroscopy after digestion in a (4:1) mixture of nitric/perchloric acid (Ministry of Agriculture, Fisheries, and Food, 1986).

Annual P balances (inputs minus outputs) were calculated for each plot over the 5-yr study. Total P inputs were the sum of P applied in fertilizer and wet deposition, which averaged 0.31 kg P ha–1 yr–1. Total P outputs were calculated as the sum of P lost in drain-flow and overland flow and removed in animal live-weight gain. Phosphorus removed in animal product was calculated as 8 g kg–1 empty body weight, which was determined by dividing the animal live-weight gain by 1.09 (Agricultural Research Council, 1980).

Data Analysis
Logarithmic transformations were performed on the weekly concentrations of DRP, DUP, and PP in drainage water before analysis against weekly drainage flow because their distribution was found to be non-normal using the Kolmogorov-Smirnov test (Mood et al., 1974). The average monthly P concentrations in drainage water were analyzed by ANOVA to test for differences between summer months (May–September) and winter months (October–April), having adjusted for yearly variation. A linear contrast was estimated to test whether there was an increase in P concentrations with the rate of fertilizer P applied. The Seasonal Kendall test (Hirsch et al., 1982) was used to determine whether there was an overall trend in average monthly P concentrations in drainage water and flows for each plot, from March 2000 to February 2005, having adjusted for seasonal variation.

Linear Mixed Model of Phosphorus Concentrations in Overland Flow
The P concentrations in overland flow were logarithmically transformed to normalize the distribution. A linear mixed model (McCulloch and Searle, 2001) was used to determine if differences occurred between plots:

Formula 3[3]
where LPrp is the log of P concentration for the overland flow event r on plot p; ar is a random intercept for each overland flow event (identified by the date of occurrence), which is distributed as N(ar,{sigma}a2); xp is the estimated fixed effects for each plot; yc is the estimated fixed effects of the years; and {varepsilon} is random variation which is normally distributed N(0,{sigma}2).

A residual maximum likelihood (Patterson and Thompson, 1971) method was used to obtain estimates for the marginal means and SEs. Paired comparisons of the marginal means for the different plots were estimated, and the significance values were adjusted for multiple comparisons using the Bonferroni method (Bonferroni, 1928).

Results

Over the 5-yr study, the average annual rainfall (March–February) and actual evaporation were 895 mm and 520 mm, respectively. The hydrologically effective rainfall (surplus rainfall) was calculated by difference and ranged from 625.7 mm in 2002/2003 to 160.1 mm in 2003/2004 (Table 1). The Olsen-P status of the plots receiving 0, 10, 20, 40, or 80 kg P ha–1 yr–1 was 23.5, 19.3, 21.9, 22.6, and 20.7 mg kg–1 (equivalent to 16.0, 13.1, 14.9, 15.4, and 14.1 mg P L–1), respectively, in February 2000. In February 2005, this had increased to 35.3, 41.2, 55.9, and 98.5 mg kg–1 for the plots receiving 10, 20, 40, and 80 kg P ha–1 yr–1, respectively (Fig. 1 ). The change in weekly Olsen-P status in the plots receiving zero P fertilizer and 80 kg P ha–1 yr–1 is shown in Fig. 2 from March 2000 to February 2005.


View this table:
[in this window]
[in a new window]

 
Table 1. Monthly rainfall from March 2000 to February 2005.

 

Figure 1
View larger version (31K):
[in this window]
[in a new window]

 
Fig. 1. Olsen-P (mg kg–1) status in February each year for plots receiving 0, 10, 20, 40, or 80 kg P ha–1 yr–1.

 

Figure 2
View larger version (33K):
[in this window]
[in a new window]

 
Fig. 2. Change in weekly Olsen-P (mg kg–1) status in the plots receiving zero P and 80 kg P ha–1 yr–1.

 
The yearly increase in Olsen-P predicted using Eq. [2] is shown in Table 2 for each plot. There was a highly significant (P < 0.001) increase in annual Olsen-P status in all plots receiving P fertilizer, which increased with the rate applied. In the plot where no fertilizer P was applied, there was a small but significant (P < 0.001) annual decrease in Olsen-P. The variability in Olsen-P status did not change with time in the zero P plot. There was a significant increase in variability of Olsen-P with time in the plots receiving P fertilizer, particularly at the highest application rate of 80 kg P ha–1 yr–1 (P < 0.001) (Fig. 2).


View this table:
[in this window]
[in a new window]

 
Table 2. Predicted increase in Olsen-P per year using data from March 2000 to February 2005 (estimated from Eq. [2]).

 
Increasing the rate of P fertilizer applied significantly increased (P < 0.001) the P concentration in the herbage. The 5-yr averages for the plots receiving 0, 10, 20, 40, or 80 kg P ha–1 were 3.2, 3.4, 3.6, 3.8, and 4.1 g P kg–1, respectively. Herbage P concentrations were lower during the summer months, but there was no evidence of P deficiency when no P fertilizer was applied. Herbage analysis throughout the 5-yr study showed satisfactory concentrations of other major nutrients. The average N, K, Mg, and S concentrations over all plots were 31, 25, 1.8, and 2.7 g kg–1, respectively. The Diagnosis and Recommendation Integrated System methodology was used to confirm that there was no nutrient imbalance in the herbage from any of the plots (Bailey et al., 1997a, 1997b).

Land Drainage Water
Generally there was a negative correlation between weekly DRP, DUP, and PP concentrations and weekly drainage flow (results not shown), probably reflecting a dilution effect. Figure 3 shows the weekly DRP concentrations in drainage water, for each calendar month, for the plots receiving zero P and 80 kg P ha–1 yr–1 for 5 yr. An ANOVA showed that there was significant (P < 0.001) seasonal variation in DRP concentrations for all plots, with higher concentrations occurring between May and September compared with the rest of the year (Fig. 3). This increase may be related to lower flows during the May to September period compared with the winter months. The seasonal difference in DRP concentrations was significantly greater for the plot receiving 80 kg P ha–1 yr–1 compared with the other plots. A linear contrast showed that there was a significant increase (P < 0.05) in DRP concentrations with the rate of P applied. There was also significant (P < 0.001) seasonal variation in DUP and PP concentrations in the drainage water of all plots with higher concentrations occurring during the summer months (May to September) compared with the other months. However, in contrast to the DRP concentrations, there was no significant difference between plots.


Figure 3
View larger version (18K):
[in this window]
[in a new window]

 
Fig. 3. Weekly dissolved reactive phosphorus (DRP) concentrations in land drainage water for each calendar month for the plots receiving (a) zero P and (b) 80 kg P ha–1 yr–1 for 5 yr (March 2000 to February 2005). The line is drawn through the mean values.

 
The Seasonal Kendall test (Hirsch et al., 1982) indicated that there was a significant upward trend in DRP concentrations in drainage water from March 2000 to February 2005 from the plots receiving 40 and 80 kg P ha–1 yr–1 (Table 3). Significantly lower (P < 0.05) weekly drain-flow with time may have been partly responsible for the increase in DRP concentrations for the plot receiving 40 kg P ha–1 yr–1 (results not shown). However, this was the only plot showing a significant trend in weekly drain-flow with time. It is not known why the hydrology of this plot was different to the other plots. The increase in DRP concentrations with time for the plot receiving 80 kg P ha–1 yr–1 was not flow related. Figure 4 shows the mean weekly difference in DRP concentrations in land drainage water for the period March 2002 to February 2005 compared with the first 2 yr (March 2000 to February 2002) of the study. There was a significant (P < 0.001) positive shift in concentrations for the plots receiving 40 and 80 kg P ha–1 yr–1. However, there was no significant trend in TP, DUP, or PP concentrations with time for any of the plots (results not shown).


View this table:
[in this window]
[in a new window]

 
Table 3. Seasonal Kendall test results of average monthly dissolved reactive phosphorus concentrations in land drainage water from March 2000 to February 2005.

 

Figure 4
View larger version (35K):
[in this window]
[in a new window]

 
Fig. 4. Histograms showing the mean weekly differences in DRP concentrations for the period March 2002 to February 2005 compared with the period March 2000 to February 2002 in land drainage water from plots receiving 0, 10, 20, 40, or 80 kg P ha–1 yr–1.

 
The CV (the ratio of SD divided by the mean) in DRP concentrations increased from 0.23 for the zero P plot to 0.30 for the plot receiving 80 kg P ha–1 yr–1. The CV in TP concentrations increased from 0.32 to 0.39 for the zero P and 80 kg P ha–1 yr–1 plots, respectively. There was no apparent trend in the CV in DUP or PP, which averaged 0.27 and 0.35, respectively. In contrast to the Seasonal Kendall test, which showed a significant increase in DRP concentrations in land drainage water with time on the plots receiving high fertilizer P application rates, AFWM concentrations did not show a consistent relationship with P application rate and were highly variable from year to year (Table 4). For example, on the plot receiving zero P fertilizer, AFWM concentration of DRP and TP ranged from 19.6 to 56.1 µg P L–1 and from 50.2 to 111.2 µg P L–1, respectively, over the 5 yr of the study. It was thought that annual median concentrations may have shown a better relationship with soil P status because they are a superior measure of central tendency than mean values and are resistant to the effects of outliers (Smith et al., 2003). This was not the case.


View this table:
[in this window]
[in a new window]

 
Table 4. Five-year annual flow-weighted minimum, mean, and maximum dissolved reactive phosphorus (DRP) and total P (TP) concentrations in land drainage water.

 
Drainage flows were highly correlated with rainfall (P < 0.001). However, as reported previously (Watson et al., 2000), there were differences in drain efficiencies between plots, which varied from year to year. Therefore, annual loads of P lost to drainage water were estimated by multiplying the AFWM concentration by the hydrologically effective rainfall (rainfall minus actual evaporation summed daily over the year). The range of annual DRP and TP loads in drainage water are shown in Table 5 for each plot over the 5-yr period. The highest annual DRP and TP loads occurred in 2002/03 because this was the wettest year of the study. Similarly to the AFWM concentrations, there was no consistent relationship between annual P loads in land drainage water and fertilizer P applied. In 4 out of the 5 yr, annual DRP and TP loads were higher from the plot receiving 40 kg P ha–1 yr–1 than the plot receiving 80 kg P ha–1 yr–1.


View this table:
[in this window]
[in a new window]

 
Table 5. Five-year minimum, mean, and maximum load of dissolved reactive phosphorus (DRP) and total phosphorus (TP) in land drainage water.

 
Overland Flow
There was considerable variation in the number of overland flow events from year to year. March 2002 to February 2003 was the wettest year of the 5-yr study with 1132 mm rainfall. High rainfall occurred in May, June, July, October, and November of 2002 (Table 1). From March 2002 to February 2003, there were 64 overland flow events, giving a total collection volume (average for plots) of 439000 L per plot. The overland flow collectors intercepted 35% of the hydrologically effective rainfall. March 2003 to February 2004 was the driest year, with 678 mm rainfall and only 14 overland flow events, giving a total collection volume (average for plots) of 36600 L per plot. This was equivalent to 11% of the hydrologically effective rainfall. Five of the 14 overland flow events occurred in May 2003, which was a relatively wet month with 100 mm of rainfall (Table 1).

The CV in DRP and TP concentrations in overland flow showed a marked increase with increasing fertilizer P application (Table 6). The CV in DUP and PP concentrations were higher in overland flow than in land drainage and averaged 0.42 and 0.50, respectively. The annual flow-weighted mean DRP and TP concentrations in overland flow were considerably higher than in land drainage water and increased with the rate of P fertilizer applied (Table 7). Maximum DRP and TP concentrations were generally recorded when overland flow occurred soon after fertilizer P application. For example, a maximum DRP concentration of 21500 µg P L–1 was recorded in overland flow from the plot receiving 80 kg P ha–1 yr–1 on 4 May 2003. Fertilizer P had been applied 2 d earlier. This single event accounted for 52% of the annual DRP load in overland flow for this plot.


View this table:
[in this window]
[in a new window]

 
Table 6. Coefficient of variation in dissolved reactive phosphorus (DRP) and total phosphorus (TP) concentrations in overland flow.

 

View this table:
[in this window]
[in a new window]

 
Table 7. Five-year annual flow-weighted minimum, mean, and maximum dissolved reactive phosphorus (DRP) and total phosphorus (TP) concentrations in overland flow.

 
A linear mixed model approach was used to determine the significance of pairwise differences between the plots for the Log10 of the mean DRP concentrations in overland flow (Table 8). The mean DRP concentration from the plot receiving 80 kg P ha–1 yr–1 was significantly higher (P < 0.001) than all other plots. Generally, there was no significant difference between plots in DUP or PP concentrations in overland flow.


View this table:
[in this window]
[in a new window]

 
Table 8. Multiple comparisons table showing the significance of pairwise differences between the plots for the log10 of mean dissolved reactive phosphorus (DRP) concentrations in overland flow (averaged over all years).{dagger}

 
The contribution of the various P fractions to total P loss was calculated on a weekly basis and averaged over the 5-yr period for each plot. A linear mixed model approach was used to determine differences between plots. The percentages of total P as DRP and PP in land drainage water were 48.2 and 31.1%, respectively, averaged over all plots. Although there were small differences between plots, there was no consistent trend with P application rate. There was no significant difference in the percentage of DRP and PP to total P lost in land drainage water between a relatively wet year (2002/03) and a dry year (2003/04). However, in overland flow, the contribution of DRP to total P was significantly (P < 0.001) higher on the plot receiving 80 kg P ha–1 yr–1 compared with the plots receiving 0, 10, or 20 kg P ha–1 yr–1. The percentage of total P as DRP on the plots receiving 80 and 10 kg P ha–1 yr–1 averaged 61.3 and 45.3%, respectively, over the study period. The percentage of PP in overland flow was significantly (P < 0.01) higher in a relatively wet year (2002/03) compared with a dry year (2003/04), and overall plots averaged 30.2 and 21.2%, respectively. This probably reflected a greater removal of soil particles from the soil surface with high rainfall intensity.

Total Phosphorus Loads in Land Drainage and Overland Flow
The annual average TP load (kg P ha–1 yr–1) in overland flow and drainage water for each plot for the 4-yr period (March 2001 to February 2005) is shown in Fig. 5 . Year 2000 data could not be included in Fig. 5 because the volume of overland flow was not measured during 2000. Total P lost (land drainage plus overland flow) was highly variable from year to year. On the zero P plot, TP lost ranged from 0.19 kg P ha–1 in 2003/04 to 1.55 kg P ha–1 in 2002/03. On the plot receiving 80 kg P ha–1 yr–1, the range was from 0.35 kg P ha–1 in 2001/02 to 2.94 kg P ha–1 in 2002/03. Consistent trends in TP loss with soil P status were not obvious. However, when all years were averaged, TP losses were higher from the plots receiving the highest fertilizer P applications. There was a problem in recording the overland flows on the plot receiving 40 kg P ha–1 yr–1 during the summer of 2002 due to instrumentation failure. This resulted in anomalously low P losses in that year and is the reason why the average TP load in overland flow from this plot was lower than the plot receiving 20 kg P ha–1 yr–1 (Fig. 5).


Figure 5
View larger version (13K):
[in this window]
[in a new window]

 
Fig. 5. Average annual total phosphorus (TP) load (kg P ha–1) in overland flow and land drainage water (March 2001 to February 2005) for plots receiving 0, 10, 20, 40, or 80 kg P ha–1 yr–1. P loads in overland flow from the plot receiving 40 kg P ha–1 yr–1 were anomalously low due to instrumentation failure during the summer of 2002. The average weekly Olsen-P concentration (mg P kg–1), from March 2001 to February 2005, is shown above each plot.

 
The contribution of overland flow to total annual P loss was generally similar to that from land drainage for the plot receiving no P fertilizer. However, at the highest application rate (80 kg P ha–1 yr–1), total P loss in overland flow exceeded that observed in land drainage by a factor of nearly two. For example, for the plot receiving 80 kg P ha–1 yr–1, total P lost in overland flow for 2002/03 was 1.91 kg P ha–1, compared with 1.03 kg P ha–1 from land drainage.

Phosphorus balances (inputs of fertilizer P and wet deposition minus outputs in drain-flow, overland flow, and removal in animal live-weight gain) for the plots receiving 0, 10, 20, 40, and 80 kg P ha–1 yr–1 averaged –8.4, 1.2, 10.8, 29.5, and 71.4 kg P ha–1 yr–1, respectively, over the 5-yr study period. Regression analysis showed a good linear relationship (R2 0.999) between annual P balance and the observed annual increase in Olsen-P status in the 0- to 75-mm soil depth as shown in Eq. [4]:

Formula 4[4]
where {Delta}Py is the yearly change in Olsen-P (mg P kg–1) and Pb is the annual P balance.

Discussion

In the current study, the Olsen-P status of grazed grassland swards was increased over a 5-yr period by surface applying different rates of P fertilizer. Regression analysis predicted that every 10 kg P ha–1 yr–1 surplus in the P budget would increase Olsen-extractable P in the top 75 mm of soil by 2.5 ± 0.17 mg P kg–1 yr–1. Assuming a soil bulk density in the field of 1 g cm–3 (Cruickshank, 1997), this indicates that 18.9% of the P surplus in soil to a depth of 75 mm is Olsen extractable. This increase in Olsen-P influenced P concentrations in land drainage water and overland flow.

There was significant seasonal variation in all P fractions in land drainage water, with the highest concentrations occurring between May and September each year. This may be related to low flows occurring during these months compared with other times of the year. Because the seasonal effect in DRP concentrations was greatest on the plot receiving 80 kg P ha–1 yr–1, compared with the other plots, it suggested that part of the seasonal increase may be due to fertilizer P applications. On the plot receiving 20 kg P ha–1 yr–1, daily TP concentrations in land drainage water were measured for part of the study period. A maximum concentration of 2674 µg P L–1 occurred on 10 Sept. 2004, 3 d after fertilizer P application and heavy rainfall (18 mm d–1) even though the soil was at a moisture deficit of 24.8 mm. This provides evidence of preferential flow pathways facilitating the loss of fertilizer P from the soil surface to drainage water at 1 m deep. Preferential flow pathways have been reported in many soil types, especially clay soils that are susceptible to cracking (Simard et al., 2000; Turner and Haygarth, 2000; Stamm et al., 2002).

Total phosphorus lost in land drainage water ranged from 0.12 to 1.49 kg P ha–1 yr–1 and was agronomically small relative to the amount of P being applied (<3.7%). These small losses are important factors in the eutrophication process leading to AFWM concentrations of TP of up to 237 µg P L–1. These concentrations are well above 100 µg P L–1, which would be classified as hypertrophic if they occurred in a lake environment (Organisation for Economic Cooperation and Development, 1982), but below the guideline of 300 µg P L–1 set by the European Union for water abstraction for drinking water purposes (Newman, 1988).

Dissolved reactive phosphorus and TP concentrations in overland flow were considerably higher than in land drainage water and increased with P application rate. Concentrations in overland flow were highly dependent on the timing of the event relative to fertilizer application. A maximum DRP concentration of 21500 µg P L–1 was recorded from the plot receiving 80 kg P ha–1 yr–1 when a large flow event occurred 2 d after fertilizer P was applied. Single overland flow events can account for a large percentage (up to 52%) of the annual DRP load. Overland flow is a major loss pathway for P from grazed grassland swards, and its contribution to total P loss can be well above that from land drainage water. There is substantial evidence linking high levels of soil P to increased losses of P in overland flow from other studies (Sharpley, 1995; Pote et al., 1999; McDowell et al., 2002), but few have quantified the relative contribution of overland flow to total P loss. Management practices should focus on controlling dissolved P losses in overland flow by considering not only the rate and timing of fertilizer applications but also on strategies to improve soil structure and hence water infiltration.

In contrast to other studies (Heckrath et al., 1995; Hesketh and Brookes, 2000), it was difficult to identify a clear Olsen-P concentration above which TP losses increased, despite the Olsen-P status in the soils ranging from 22 to 99 mg P kg–1, after 5 yr of fertilizer P applications. Any relationship was confounded by annual variation in hydrologic events and flows as well as hydrologic differences between plots. This emphasizes the importance of medium- to long-term studies. Phosphorus concentrations in land drainage water were significantly increased with fertilizer P application rates of 40 kg P ha–1 yr–1 or above. However, P concentrations in overland flow were increased with P applications above 10 kg P ha–1 yr–1.

There was little change in the Olsen-P status of the plot receiving zero P fertilizer over the 5 yr of the study (decrease of 1.10 mg P kg–1 yr–1 or 0.75 mg P L–1 yr–1) and no evidence of a decrease in DRP or TP concentrations in overland flow or land drainage water with time. AFWM concentrations of TP in drainage water from the zero P plot ranged from 50 to 111 µg P L–1 and were well above the 35 µg P L–1 believed to trigger eutrophic effects in lakes (Organisation for Economic Cooperation and Development, 1982). AFWM concentrations of TP in overland flow from this plot ranged from 400 to 1021 µg P L–1. Withholding P fertilizer for 6 yr (5 yr of study plus the equilibration year of 1999/2000) from a grazed grassland sward at the optimum Olsen-P status (24–37 mg kg–1 or 16–25 mg P L–1) was not long enough to lower soil P losses to drainage water or overland flow or to have an adverse effect on herbage P concentrations.

The average annual P balance (inputs minus outputs) for the plot receiving zero P was a deficit of 8.4 kg P ha–1. A nutrient management scheme was introduced for 3 yr (1998–2000) in the Greenmount Agricultural and Horticultural College catchment, which reduced the P balance (inputs minus outputs) from a surplus of 23.5 kg P ha–1 yr–1 to a deficit of 15.0 ha–1 yr–1 (Smith et al., 2003). This led to a significant reduction in concentrations of DRP in land drainage in 2000 compared with 1997. The major difference between the two studies was the removal of P in cut silage (18 kg P ha–1 yr–1) in the Greenmount catchment. In the current grazed study, total P output (loss to water plus P removal in animal product) averaged 8.8 kg P ha–1 yr–1.

Having built up a range of Olsen-P concentrations in the soils (from 22 to 99 mg P kg–1) over a 5-yr period, an interesting question is how quickly the Olsen-P status of the soils declines when no fertilizer P is applied and how loss rates of P to water are affected. Recent soil samples from intensive grassland farms in Northern Ireland showed that 51% of silage fields were above the agronomic optimum for Olsen-P (>37 mg kg–1 or >25 mg P L–1) (Jordan, personal communication, 2006). For Northern Ireland to comply with a range of European Union Directives (European Community, 1991; European Commission, 2000) on the aquatic environment, it has declared a "total territory" approach to control the agricultural contribution to nutrient losses. Part of the proposed action plan includes withholding fertilizer P application on soils of high Olsen-P status. Curtailing fertilizer P applications at the current site should demonstrate how rapidly environmental improvements can be achieved under grazing management and provide information to support policy decisions. Results suggest that improving water quality with respect to P may be a slow process unless P is redistributed within the soil profile by, for example, ploughing for the reseeding of pastures.

ACKNOWLEDGMENTS

We thank the staff in the Agriculture, Food, and Environmental Science Division, Agri-Food and Biosciences Institute, particularly Phil Dinsmore's team in the freshwater chemistry laboratory for their excellent technical support with sample collection and analysis. We also thank the staff at the Agricultural Research Institute for Northern Ireland at Hillsborough for managing the experimental site and Raymond Stewart and Briege McCarney for the collection of soil and herbage samples.

NOTES

All rights reserved. No part of this periodical may be reproduced or transmitted in any form or by any means, electronic or mechanical, including photocopying, recording, or any information storage and retrieval system, without permission in writing from the publisher.

REFERENCES




This article has been cited by other articles:


Home page
J. Environ. Qual.Home page
C. van der Salm, W. J. Chardon, G. F. Koopmans, J. C. van Middelkoop, and P. A.I. Ehlert
Phytoextraction of Phosphorus-Enriched Grassland Soils
J. Environ. Qual., February 25, 2009; 38(2): 751 - 761.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Watson, C. J.
Right arrow Articles by Matthews, D. I.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Watson, C. J.
Right arrow Articles by Matthews, D. I.
Agricola
Right arrow Articles by Watson, C. J.
Right arrow Articles by Matthews, D. I.
Related Collections
Right arrow Surface Water Quality
Right arrow Phosphorus
Right arrow Field-Scale Studies


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
The SCI Journals Agronomy Journal Crop Science
Journal of Natural Resources
and Life Sciences Education
Vadose Zone Journal
Soil Science Society of America Journal Journal of Plant Registrations The Plant Genome