|
|
||||||||
a USDA-ARS, National Soil Tilth Laboratory, Ames, IA 50011
b University of Georgia, Department of Crop and Soil Sciences, Athens, GA 30602
* Corresponding author (schroeder{at}nstl.gov)
Received for publication December 2, 2003.
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: DRP, dissolved reactive phosphorus R1, sufficient rainfall to produce 30 min of runoff immediately after litter application R2, no rainfall for 30 d after manure application, then sufficient rainfall to produce 30 min of runoff R3, small rainfall events every 7 d (5 min at 75 mm h1) for 30 d and then sufficient rainfall to produce 30 min of runoff STP, soil test phosphorus TP, total phosphorus WSP, water-soluble phosphorus
| INTRODUCTION |
|---|
|
|
|---|
The state of Georgia is one of the top broiler chicken (Gallus gallus domesticus) producing regions of the United States with 1.29 billion broilers raised in 2002 (National Agricultural Statistics Service, 2002). Besides broiler production, agriculture in central and northern Georgia is generally limited to beef cattle and hay production. Thus, broiler, cattle, and hay production are often integrated, and broiler litter (manure and bedding) serves as an organic fertilizer on pastures and hayfields. Traditionally, broiler litter is applied to meet forage N needs, but surface application greatly increases the risk of P loss in surface runoff because of low (2:1) N to P ratio. Forms of P in runoff include dissolved inorganic and organic P forms and particulate P associated with mineral or organic particles transported in runoff. Because erosion rates from hayfields and pastures are low, the dissolved fraction is usually the dominant form of P in runoff from hayfields and pastures (Sharpley et al., 1992).
A strong relationship exists between the rate of manure application and the concentration of total phosphorus (TP), dissolved reactive phosphorus (DRP), and particulate P in runoff (McLeod and Hegg, 1984; Mueller et al., 1984; Edwards and Daniel, 1994; Vervoort et al., 1998; Wood et al., 1999; Pote et al., 2001; Kleinman and Sharpley, 2003) where P concentrations in runoff typically increase with an increase in manure application rate. The majority of the published research shows that most of the P is lost in the first runoff event from fields where manure has been surface-applied (Edwards et al., 1994; Sharpley, 1997; Sauer et al., 1999). Several researchers have shown that after an initial spike, P concentration in runoff declines with time or number of rainfall events, often remaining greater than background concentration for an extended period (Heathman et al., 1995; Sharpley, 1997). In a multiple-year study, Pierson et al. (2001) reported increased P concentration in runoff from natural rainfall for up to 18 mo after an application of poultry litter to a small watershed in Georgia. In contrast, some small plot studies have shown P concentrations in runoff were similar to background concentrations after two artificial rainfall events (Edwards and Daniel, 1994; Sauer et al., 1999).
Because most of the P lost from manure applications is lost in the first runoff event, several researchers investigated the effect of time between manure application and a runoff-producing rainfall on P loss. Several small plot or lab studies have shown a decrease in P concentration with an increase in time to a runoff-producing rainfall event with liquid poultry manure (Westerman and Overcash, 1980) or incorporation of poultry litter into the soil (Sharpley, 1997). However, when dry poultry litter is surface-applied to pastures and hayfields, a surface layer of thatch is likely to prevent direct contact between the litter and the soil, reducing the possibility that P in the manure will be adsorbed immediately by the soil. Edwards et al. (1994) found that time to artificial rainfall producing runoff (414 d) had no effect on either the concentration or mass of ortho-P or TP in runoff. In a more recent study, Pierson et al. (2001) observed that P concentrations from 0.75 ha were less than 5 mg L1 DRP for intervals of up to 7 mo and greater than 5 mg L1 DRP when events occurred within 12 d of poultry litter application.
The objectives of this research were twofold. The first was to determine the effects of application rate and initial runoff timing on the long-term loss of P from poultry litter surface-applied to pastures and hayfields. The second was to develop equations that may be used to predict P loss from surface-applied poultry litter and evaluate the effectiveness of these equations with runoff data from a reference site.
| MATERIALS AND METHODS |
|---|
|
|
|---|
The three litter rates and rainfall scenarios were applied to 1- x 2-m plots in a 3 x 3 randomized complete block design with three replications. In addition to these treatments, three control plots were included (one for each block) to allow for correction of P loss not associated with litter application. Therefore, the experiment consisted of thirty plots (3 x 3 x 3 + 3 = 30). Following the implementation of the initial rainfall scenarios, all plots received sufficient rainfall to produce 30 min of runoff on a biweekly basis for 5 mo from May through October 2001. Due to several days of unseasonably warm weather in January 2002, we had the opportunity to conduct additional rainfall simulations on the 13 Mg ha1 plots and control plots.
The experimental site was a hayfield with a fairly uniform 8% slope at the University of Georgia Plant Sciences Farm near Athens, GA. The soil in the study area is a Cecil sandy loam (fine, kaolinitic, thermic Typic Kanhapludults) with average pH of 5.7 and Mehlich IIIextractable P of 22 mg kg1 in the A horizon. In late February 2001, Kentucky 31 tall fescue (Festuca arundinacea Schreb.) was planted on the site to supplement the relatively thick stand of fescue present. At that time, 14312 (NPK) starter fertilizer was applied at a rate of 448 kg ha1. Approximately 48.5 mm of natural rain fell on the experimental area between February and May.
From 3 Mar. to 1 May 2001, thirty 1- x 2-m plots were installed at the site in three blocks of 10 plots. The blocks were positioned so that the long axis of all plots was oriented down slope. Plot borders, consisting of 15-cm-tall, approximately 0.3-cm-thick sheet metal, were pressed into the ground to a depth of at least 7 cm to isolate runoff. Aluminum flumes were installed at the down-slope edge of each plot to divert surface runoff to a collection point. The litter used in this study was collected on 10 May 2001 from a broiler farm in northern Georgia. Poultry litter was applied by hand on 14, 15, and 16 May to plots in Blocks 1, 2, and 3, respectively.
During artificial rainfall simulations, collection of runoff from the field plots began after steady runoff commenced and continued for 30 min. Runoff was collected in toto, and a 500-mL subsample was taken and immediately placed on ice. Clear polymer covers were placed over each plot between simulations to prevent the plots from receiving natural rainfall. These covers consisted of polyethylene film stretched over wooden frames that were pitched to direct rainwater away from the plots. Wooden blocks, 10 cm tall, were used to support the covers above the plots. This resulted in a gap of at least 10 cm at each end of the plot. Paired thermocouples (in and out of the plot) installed in five plots showed that the soil temperature under the covers differed from ambient soil temperature by less than 0.5°C throughout the experiment. All plots were mowed to a height of approximately 10 cm every 2 wk for the duration of the experiment, and clippings were removed to prevent the loss of P to runoff from the decaying grass. Local well water was used as the water source (P concentration < 0.01 mg L1). After the last rainfall event, soil samples (composites of 10 random samples) were collected from the 0- to 10-cm depth within each plot.
Sample Analysis
Soil samples were collected from the surface 10 cm of each plot following the final runoff event. Samples were air-dried and ground to pass a 2-mm sieve. Soil pH was determined in a 1:2 soil to water mixture using a glass electrode. Water-soluble and Mehlich IIIextractable P was determined (Mehlich, 1984; Pote et al., 1996). Total P in unfiltered runoff samples was determined colorimetrically (Murphy and Riley, 1962) following micro-Kjeldahl digestion (Baker and Thompson, 1992). Immediately after collection, 125 mL of each runoff sample was filtered (0.45-µm pore diameter) to remove particulate matter and stored at 20°C until analyzed. The DRP concentration of filtered runoff water samples was also determined colorimetrically (Murphy and Riley, 1962).
Total P content of the poultry litter was determined colorimetrically following micro-Kjeldahl digestion (Baker and Thompson, 1992). Water-soluble phosphorus (WSP) in poultry litter was determined colorimetrically after shaking 20 g of manure in 4 L of deionized water for 4 h. Total P, water-soluble P, pH, and moisture content of the poultry litter were 23.98 g kg1, 6.2 g kg1, 8.46, and 54.2%, respectively. The pH of the litter was determined in a 1:5 litter to water mixture using a glass electrode. Litter moisture content was determined after oven-drying at 65°C for 48 h.
Reference Site
Data used to assess the effectiveness of P loss prediction equations presented as "observed" data were originally collected and reported by Pierson et al. (2001). In the Pierson et al. (2001) study, five fescuecommon Bermudagrass [Cynodon dactylon (L.) Pers.] fields (0.720.79 ha) were bordered by earthen berms and fitted with H-flumes and Isco (Lincoln, NE) refrigerated samplers. Soil series present at the Eatonton, GA, sites include Cecil, Altavista (fine-loamy, mixed, semiactive, thermic, Aquic Hapludults), Helena (fine, mixed, semiactive, thermic Aquic Hapludults), and Sedgefield (fine, mixed, active, thermic Aquultic Hapludalfs). Precipitation and runoff volume were recorded at 5-min intervals. During the two years studied, poultry litter was applied four times: 16 Mar. 1995 (102 kg P ha1), 30 Oct. 1995 (112 kg P ha1), 4 Mar. 1996 (174 kg P ha1), and 25 Sept. 1996 (103 kg P ha1). Litter samples were analyzed for TP and WSP by Kuykendall et al. (1999) by the same methods used in the present study. Runoff samples were filtered (0.45 µm) and analyzed for DRP by the molybdate blue method (Murphy and Riley, 1962). Thirty-nine runoff events from the reference site from 1 Jan. 1995 to 31 Dec. 1996 were used as observations for modeling purposes. The initial Mehlich-I soil test phosphorus (STP) level reported by Pierson et al. (2001) of 13 mg kg1 was converted to 25 mg kg1 Mehlich III (Mehlich, 1984; Shuman et al., 1988).
Because runoff P concentration can be strongly influenced by P associated with the soil, as well as P from surface-applied litter, any attempt to model P loss must include some estimate of the P contribution from the soil P pool. The contribution of P from the soil P pool, estimated by STP, was modeled with the following equation (Schroeder et al., 2004), which describes the relationship between STP (mg kg1) and P in runoff from similar soils:
![]() | [1] |
Statistical Analysis
Total and soluble P mass losses from each plot were calculated for each runoff event using P concentration and runoff volume. Overall total P loss and cumulative P losses by runoff event were calculated. Average TP and DRP losses from control plots were subtracted from treatment plot P loss so that only P loss associated with litter application would be analyzed. Statistical analyses were performed using the Statistical Analysis System (SAS Institute, 1994). Analysis of variance (ANOVA) techniques were used to determine treatment effects and to check for interaction across all 10 runoff events. The least significant difference method was used to separate treatment means. Nonlinear regression was used to develop predictive equations relating P loss in runoff from surface-applied poultry litter to P application rate, runoff depth, cumulative rainfall, days since manure application, antecedent soil water content, and temperature. Regression analysis was also employed to determine if STP levels were related to P application rate, runoff depth, cumulative rainfall, or pH.
To evaluate modeling performance, the following measures were used: correlation coefficient (r), a measure of the linear correlation between observed and simulated results; root mean square error (RMSE), an estimate of the inherent error in the simulation; and the relative RMSE (RRMSE = RMSE/observed mean x 100), a measure of error in relation to the mean. In addition to the above analysis, we also regressed observed against simulated results and analyzed the intercepts and slopes to determine if they were different from 0 and 1, respectively (SAS Institute, 1994).
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
Cumulative TP and DRP loss were greatest from the treatment (R1) with artificial rainfall producing runoff almost immediately after litter application (Table 1). The other treatments (R2 and R3) did not produce significantly different cumulative TP or DRP losses (p = 0.05). The fact that the most realistic initial rainfall treatment (R3) produced less cumulative TP and DRP loss (4.08 and 3.38 kg ha1, respectively) than R1 (7.40 and 5.62 kg ha1, respectively) suggests that under "real world" conditions (R3) P losses from surface-applied manure may be considerably less than the "worst case" scenario (R1). The effect of rainfall timing was most pronounced in the first runoff event where the R1 treatment showed both the greatest TP and DRP loss (4.41 and 3.22 kg ha1, respectively) and the highest percentage P loss (59.6 and 57.3%, respectively). The R3 treatment produced the smallest TP and DRP loss as well as the smallest percent TP and DRP losses in the first event. Over the remainder of the events there was little difference in P loss or percentage P loss among the rainfall treatments.
|
As expected, the 13 Mg ha1 litter treatments produced much higher cumulative TP and DRP loss than the 2 or 7 Mg ha1 treatments (Table 2). The greatest litter application rate produced the greatest TP and DRP loss for all individual runoff events (Table 2). The first runoff event showed the most dramatic differences among the three litter application rates with the 13, 7, and 2 Mg ha1 rates producing TP losses of 4.8, 2.4, and 1.3 kg ha1, and DRP losses of 3.6, 2.0, and 0.7 kg ha1, respectively. By the second runoff event, P losses from all three treatments decreased considerably to 1.4, 0.5, and 0.3 kg ha1 TP and 1.0, 0.4, and 0.3 kg ha1 DRP. For the second and later events DRP losses from the 7 and 2 Mg ha1 treatments were not significantly different.
|
13 Mg ha1). In fact, Pierson et al. (2001) observed DRP concentrations in excess of 1 mg L1 for more than 18 mo following 4 yr of poultry litter application. The fact that after 10 runoff events only 11.2, 6.2, and 5.7% of applied P was lost from the 2, 7, and 13 Mg ha1 litter treatments, respectively, combined with the relatively modest increases in STP (reported below), indicates that a significant portion of the applied P remained on or very near the soil surface. This residual surface P will continue to solubilize over time and may produce elevated P levels in runoff for a significant period (Pierson et al., 2001).
The results of this experiment both affirm and contradict previous research on P loss following poultry litter application. It is quite clear from our research that a 30-d delay between litter application and runoff significantly reduced both the initial TP and DRP loss and the overall mass of P lost when compared with runoff immediately following manure application. These results appear to agree with results published by Westerman and Overcash (1980) and Sharpley (1997). However, these two studies differ from the present study in that they either used liquid manure (Westerman and Overcash, 1980) or incorporated the manure into the soil (Sharpley, 1997). These two differences would probably amplify the delay effect because of the close contact between soil and manure. Contradictory findings were reported by Edwards et al. (1994), who concluded that delay intervals of 4, 7, and 14 d did not affect TP and DRP loss in runoff from poultry litter surface-applied to fescue plots. They reasoned that delay interval did not affect P loss because the grass cover limited contact between the litter and the soil (i.e., conditions were not optimal for soil adsorption of litter P). The differences between the Edwards et al. (1994) study and the present study are probably due to several factors, including the shorter delay interval (14 vs. 30 d), lower litter application rate (5.6 vs. 13 Mg ha1), and the fact that Edwards et al. (1994) did not include an immediate runoff treatment in their study. It is possible that the delay effect occurs in a short period due to rapid adsorption of P, so that without an immediate runoff treatment, this effect is not observed.
Soil Phosphorus Levels
After the final rainfall event, soil samples were collected from the upper 10 cm of each plot (Table 3). Soil pH ranged from 5.6 to 5.9, deionized waterextractable P ranged from 0.3 to 6.8 mg kg1, and Mehlich IIIextractable P ranged from 21.8 mg kg1 for the unamended controls to 65.9 mg kg1 for the 13 Mg ha1 treatment. Changes in Mehlich-III STP were related to the total phosphorus in the litter (LTP) applied (R2 = 0.57) by the following equation:
![]() | [2] |
|
Due to the relatively small amount of P extracted by deionized water (mean = 2.80 mg kg1) and the high degree of variability (standard deviation = 3.10 mg kg1) within treatments, no significant relationships were found between deionized waterextractable P and any of the independent variables used.
Curve Fitting
The contribution of surface-applied poultry litter to P loss from small plots was best described by a first-order decay equation:
![]() |
The coefficients related to the number of days since litter applications were similar for the R2 and R3 equations. However, the k values for the R1 equations significantly were larger than the k values of the other treatments. The larger k values for the R1 equations represent the effect of runoff occurring immediately after litter application. The fact that the k values for both the R2 and R3 equations are smaller than the k for the R1 equations indicates that, in the absence of runoff, as the number of days since application increases, P loss in runoff decreases. This effect may be due to immobilization of P in the litter over time as previously discussed. However, in soils with high STP levels, this effect may be less apparent because of reduced soil P sorption capacity.
In these equations the first term, P0A, represents the maximum P concentration in runoff when t is zero. The value of A in the equations for the R3 treatment is about half the value of A in the equations for the R1 treatments. This reflects the combined effect of the delay in runoff and the application of small rainfall events before runoff. Figure 1 shows the decay equation for TP and DRP developed for the R1 and R3 scenarios. The best fit was obtained for the R1 equations followed by the R3 equations (Table 4).
|
|
|
|
|
|
The above discussion implies that observed DRP concentration did not follow a predictable decrease over time following litter applications. Since the prediction equations used in this study were classic decay equations, they could not predict these increases in DRP concentration. The question follows, Why does the DRP concentration increase in the absence of added P, and why did the increased concentrations coincide with large volume runoff events? We propose two possible explanations that may be related.
First, variable source area (VSA) may play a role in this phenomenon. The concept of VSA is that for any given field there is a limited area that contributes runoff to a stream and that the size of this area changes over time (Dunne and Black, 1970). The size of the source area depends on the size of the storm, antecedent moisture, topography, and soil type. Based on the VSA concept, increases in DRP concentration could be due to the larger source area supplying runoff during larger storms. In effect, frequent small storms remove P from the same area, thereby depleting the soluble P pool. This results in a decrease in DRP concentration over time. However, when a large storm occurs, the VSA increases and areas that may have a large pool of soluble P contribute to runoff, increasing DRP concentration in runoff. Alternatively, the small, low DRP concentration events may be the result of overland flow from very small areas close to the collection flumes.
Second, microbial biomass turnover due to prolonged soil desiccation during periods without rain and rapid rewetting during rainfall events may contribute significantly to the size of the soluble P pool in the VSAs. In a study of C and N fluctuations, Van Gestel et al. (1993) observed that desiccation and rewetting contributed to C and N mineralization. Kieft et al. (1987), who studied microbial response to rapid increases in water potential, reported that 17 to 58% of soil biomass C was released on rapid wetting of dry soils. They concluded that a rapid water potential increase could be a potent catalyst for the turnover of soil C, as well as other nutrients such as N and P. More recently, Turner and Haygarth (2001) reported increases in water-extractable soil P of 185 to 1900% on drying and rewetting, which they attributed to microbial cell lysis. The P in cellular components, such as phospholipids and nucleic acids, that are released on cell lysis may pass through the typical 0.45-µm filter used to differentiate dissolved P from particulate P. Several studies have shown that organic P may be hydrolyzed by and react with molybdate (Ron Vaz et al., 1993; Tarapchak, 1993; Haygarth et al., 1997). In effect, the increases in DRP observed may be partly due to a flux of dissolved organic P released from lysed microbial cells.
In addition to the above-mentioned phenomena, some of the lack of model fit is probably related to the difficulties encountered when applying a model developed from small plot data to data from larger scales.
| CONCLUSIONS |
|---|
|
|
|---|
Nonlinear regression equations that may be useful in predicting P losses in runoff from surface-applied poultry litter were developed. These equations were able to explain 68 to 91% of the variability seen in P losses. Because the levels of P loss seemed to rebound after several months as reflected in the runoff from the January 2002 rainfall simulations, future research in this area should be directed at determining the effect of the interval between runoff events on P loss. The results indicate that the equations developed in this study were reasonably effective in predicting DRP loss from poultry litter that was surface-applied to the pastures and hayfields at the reference site and may be of use elsewhere under similar conditions. This prediction was most accurate for runoff events that occurred soon after litter application, which consequently are the events that produce the greatest DRP loss. The largest source of variation between observed and simulated DRP concentration was associated with instances where observed DRP concentration increased in the absence of additional P application. These counterintuitive increases in DRP concentration may be explained by a combination of processes including variable source area and microbial turnover. The results of this study also indicate that in the presence of surface-applied manure the use of equations that relate P in runoff to STP may not be appropriate. Additionally, further study of the dynamics of P cycling in the surface horizon and thatch layer of soils where poultry liter has been applied should be pursued.
| REFERENCES |
|---|
|
|
|---|
Related articles in JEQ:
This article has been cited by other articles:
![]() |
M. S. Srinivasan, P. J. A. Kleinman, A. N. Sharpley, T. Buob, and W. J. Gburek Hydrology of Small Field Plots Used to Study Phosphorus Runoff under Simulated Rainfall J. Environ. Qual., October 24, 2007; 36(6): 1833 - 1842. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. A. Volf, G. R. Ontkean, D. R. Bennett, D. S. Chanasyk, and J. J. Miller Phosphorus Losses in Simulated Rainfall Runoff from Manured Soils of Alberta J. Environ. Qual., April 5, 2007; 36(3): 730 - 741. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 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 | |||