|
|
||||||||
a Alberta Agriculture, Food and Rural Development, Irrigation Branch, 100, 54011st Avenue South, Lethbridge, AB, Canada T1J 4V6
b Agriculture and Agri-Food Canada, Lethbridge Research Centre, 54031st Avenue South, Lethbridge, AB, Canada T1J 4B1
* Corresponding author (joanne.little{at}gov.ab.ca)
Received for publication February 14, 2005.
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: CU, cultivator CUC, cultivator control CUM, cultivator manure DD, double disk DDC, double disk control DDM, double disk manure DRP, dissolved reactive phosphorus FWMC, flow-weighted mean concentration MP, moldboard plow MPC, moldboard plow control MPM, moldboard plow manure NI, no incorporation NIC, no-incorporation control NIM, no-incorporation manure STP, soil test phosphorus TN, total nitrogen TP, total phosphorus TSS, total suspended solids
| INTRODUCTION |
|---|
|
|
|---|
Rainfall and overland runoff generally interact with a very thin layer of surface soil. The depth of interaction is influenced by rainfall intensity, runoff energy, and soil slope (Ingram and Woolhiser, 1980). Sharpley (1985) reported that the effective depth of interaction between surface soil and runoff increased from 1.3 to 37.4 mm with an increase in rainfall intensity from 50 to 160 mm h1 and an increase in soil slope from 2 to 20%. Consequently, differences in soil erosion and runoff water quality are strongly related to the depth of tillage and the degree of mixing with surface soil.
Surface applications of fertilizer or manure without incorporation are extremely vulnerable to losses into surface waters, particularly when runoff occurs shortly after application (Hansen et al., 2002; Tabbara, 2003). The potential for nutrient losses in overland flow differs by manure source (Kleinman et al., 2002), and increases with the rate of manure application (Kleinman and Sharpley, 2003), decreases with time after application (Edwards and Daniel, 1994; Eghball et al., 2002), and decreases with successive rainfall events (Sharpley, 1997; Kleinman and Sharpley, 2003).
Plowing or cultivation of soils to incorporate fertilizers or manure reduces the risk of direct transmission of nutrients to surface water; however, incorporation of crop residues and soil amendments increases the potential for soil erosion. Tillage systems have been developed to maintain crop residues near the soil surface, thereby reducing soil erosion and runoff and limiting sediment and nutrient losses from agricultural soils (Hansen et al., 2002). Accumulation of nutrients from fertilizers and crop residues near the soil surface, however, can result in significant concentrations of dissolved nutrients in overland flow from these systems (Seta et al., 1993).
Previous research on the effects of different tillage methods and manure on phosphorus losses in surface runoff has been conducted with dairy, beef, or swine manure on corn (Zea mays L.) or cornsoybean [Glycine max (L.) Merr.] cropland (Mueller et al., 1984b; Ginting et al., 1998b; Bundy et al., 2001; Daverede et al., 2004). Tillage and manure application effects on nitrogen and phosphorus losses in overland flow or leaching have also been investigated with beef manure in grain sorghum [Sorghum bicolor (L.) Moench]winter wheat (Triticum aestivum L.) cropping systems (Eghball and Gilley, 1999) and in corn production systems (Zhao et al., 2001). Solid beef manure from cattle feedlots in Alberta is generally applied to silage barley stubble at nitrogen-based rates and is usually incorporated with a double disk or heavy-duty cultivator.
The objective of this study was to compare tillage methods for incorporation of beef cattle manure through evaluation of nutrient and sediment losses in surface runoff and nutrient concentrations in subsurface leachate during field rainfall simulations conducted within a week of manure application in a silage barley cropping system.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Treatments consisted of two rates of solid beef cattle manure (0 and 60 Mg ha1, wet weight) and four tillage methods (surface application without incorporation, and incorporation with a single pass of a double-disk, a heavy-duty cultivator, or a moldboard plow). Individual treatments were designated as cultivator control (CUC), cultivator manure (CUM), double disk control (DDC), double disk manure (DDM), moldboard plow control (MPC), moldboard plow manure (MPM), no-incorporation control (NIC), and no-incorporation manure (NIM). The depth of tillage was approximately 10 to 15 cm for the double-disk and cultivator, and about 25 to 30 cm for the moldboard plow. Tillage was completed (parallel to the slope) on the day of manure application. Plots were harrowed and then seeded (both perpendicular to the slope) to barley (cv. Duke) in early May each year. Crops were harvested at the silage stage in late July. Irrigation water was applied using a wheel-move sprinkler system to meet the water use requirements of the crop. Manure was applied by hand to one replicate at a time following harvest each year. A garden rake was then used to redistribute the manure uniformly within the main plot.
Beef cattle manure was obtained from the feedlot at the Lethbridge Research Station. Five manure samples were collected from the manure source before the application on each replicate. Available nitrogen (NO3N and NH4N) in the manure was determined by extracting field-moist samples (10 g manure and 200 mL of 2 M KCl) and NH4N was determined by automated salicylate (Rhine et al., 1998), while NO3N was measured by hydrazine reduction (Kempers and Luft, 1988). The remainder of each sample was air-dried, ground (<2 mm), and extracted for chemical analyses. Available phosphorus (ortho-P) extracts were prepared by shaking 1 g of manure and 25 mL of Kelowna extract (Van Lierop, 1988). Finely ground samples (<150 µm) were prepared for total C, N, and P analyses. Total P was determined using a wet-oxidation procedure (Parkinson and Allen, 1975) and ortho-P was measured on an autoanalyzer by the ascorbic acid reduction method (Murphy and Riley, 1962). Total carbon (TC) and total nitrogen (TN) were measured by the Dumas automated combustion method (McGill and Fiqueiredo, 1993) on a carbonnitrogensulfur analyzer (Carla Erba, Milan, Italy).
Before annual rainfall simulations, soil samples were collected from the 0- to 2.5-, 2.5- to 7.5-, and 7.5- to 15-cm depth intervals using an excavation method. Subsamples were taken and dried to determine antecedent soil moisture content. Surface roughness was also measured using the chain method (Saleh, 1993). Deep soil cores (00.15, 0.150.30, 0.300.60, 0.600.90, 0.901.20, and 1.201.50 m) were taken before the simulations in the fall of 2000. Soil samples were air-dried, ground (<2 mm), and analyzed. Nitrate N and NH4N were determined on extracts from 1 g soil and 25 mL of 2 M KCl. Soil test phosphorus (STP) was determined by measuring the ortho-P content of a Kelowna extract (Van Lierop, 1988). Total C, N, and P were determined on finely ground soil samples according to the same methods as the manure samples.
Rainfall Simulations
Simulated rainfall was applied to a 1- x 1-m area within the plots using a portable Guelph rainfall simulator (Tossell et al., 1987). A stainless-steel frame (0.2 m deep) was installed at a depth of about 0.10 m on the upslope end and both sides of the test area, a triangular metal apron was installed at the lower end, and a hole was excavated at the apex of the apron to allow collection of runoff in 1-L containers. A Plexiglas cover was placed over the metal apron to prevent rainfall directly onto the apron. Rainfall was applied using a Fulljet 3/8GG 20W nozzle (Spraying Systems Co., Wheaton, IL) from a height of 0.8 m and at an intensity of approximately 100 mm h1 when operated at a pressure of 96.5 kPa (Miller, 2003). This rainfall intensity for a 30-min duration has a 1-in-50-yr return period in the Lethbridge area. Deionized water was used for all simulations and tests were usually performed within 5 d of manure application. Runoff generated from the plots was collected for 30 min following initiation of continuous runoff.
Runoff samples were collected in 1-min intervals. Volume was measured in the field for each sample and total runoff depth (mm) was determined on the basis of the plot area. Ten subsamples from each run were collected for analysis at 1, 2, 3, 5, 7, 9, 12, 15, 20, and 30 min after commencement of continuous runoff. Water samples were analyzed for dissolved reactive phosphorus (DRP), TP, NO3N, NH4N, and total nitrogen (TN). Samples analyzed for DRP were centrifuged (10000 rpm for 10 min) and filtered through 0.45-µm membrane filters within a few minutes following collection. These DRP samples were stored at 4°C until analysis, which occurred within 4 d of collection. Total phosphorus, TN, NO3N, and NH4N samples were acidified and then frozen. Ammonium N, NO3N, and ortho-P (DRP) were determined by the same methods as the soil extracts. Total N and total P were determined using the persulfate digestion method (Methods 4500-Norg D and 4500-P B, respectively; American Public Health Association, 1995). A composite sample for the 30-min runoff interval was analyzed for total suspended solids (TSS) by filtering a 100-mL aliquot of runoff water through a 0.45-µm membrane filter and then drying and weighing the filter. Total mass loads were estimated by multiplying subsample concentrations by their respective volumes, and loads were interpolated for intervening times by integrating the area under the curve. Flow-weighted mean concentration (FWMC) was calculated by dividing the total mass load by the total runoff volume.
In 2001 and 2002, four 25-mm-diameter zero-tension lysimeters (Simmons and Baker, 1993) were installed in the plots at a depth of approximately 60 cm to determine the concentration of constituents in subsurface flow. Subsurface water samples were analyzed for the same chemical constituents and by the same methods as the runoff samples. Subsurface measurements were restricted to concentrations of constituents, since subsurface flow volumes were not measured in this study and mass losses could not be estimated.
Data Analysis
Runoff Analysis
Variables were log transformed to log10 (x + 1) values, as necessary, to satisfy conditions of normality. Proc Mixed (SAS Institute, 2000) was used to analyze runoff data using tillage treatment, manure application, and year as fixed effects, and replicate as a random effect. A repeated year statement using a compound symmetry model was employed for FWMCs and mass loads. Each year of the study was analyzed separately when the tillage by year or manure by year interaction was significant for a given variable. Mean comparisons were performed using Tukey's studentized range test at P < 0.05 (SAS Institute, 2000). The Univariate procedure in SAS (SAS Institute, 2000) was used to evaluate residuals for normality and to check for obvious outliers.
Soil Analysis
Soil physical and chemical data obtained for site characterization were analyzed using an analysis of variance procedure (Proc GLM) to determine the initial homogeneity of the site (SAS Institute, 2000). Treatment effects were subsequently evaluated using the same mixed-model analysis as the runoff data, with mean comparisons performed using Tukey's studentized range test at P < 0.05 (SAS Institute, 2000). Soil variables were transformed to log10 (x + 1) values, as necessary, before analysis to satisfy conditions of normality. The Univariate procedure in SAS (SAS Institute, 2000) was used to evaluate residuals for normality and to check for obvious outliers. The regression procedure (SAS Institute, 2000) was used to examine relationships between soil nutrient levels and nutrient FWMCs in runoff.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
|
Total P and TN content of the upper 2.5 cm were reasonably uniform for all of the control treatments (Fig. 1d and 1e). Soil test phosphorus in the upper 2.5 cm ranged from 35.3 to 78.0 mg kg1 and mineral nitrogen (NH4N and NO3N) content was less than 30 mg kg1 for the control treatments (data not shown). Nutrient content of the manure treatments was highest in the upper 2.5 cm and generally decreased with depth, except for the MPM treatment that had higher nutrient content with depth (data below 2.5 cm not shown). Manure treatments had STP levels that were 9 to 42% of TP and mineral nitrogen content ranged from 1 to 3% of TN. The NIM treatment generally had the highest nutrient content in the upper 2.5 cm, except for NO3N (data not shown), and the MPM treatment generally had the lowest nutrient content at the surface. The CUM and DDM treatments had similar nutrient content within the upper 15 cm (data below 2.5 cm not shown).
Surface Runoff
The rates of runoff increased gradually throughout the 30-min runoff interval and were highest for the NI treatments and lowest for the MP treatments, regardless of manure application (Fig. 2) . Runoff was greatest from NI treatments and was reduced by 20% with a double disk (DD), by 34% using a cultivator, and by 56% with a moldboard plow. Mean runoff depths for the different tillage treatments during the 30-min rainfall simulations ranged from 10.6 to 24.1 mm (Fig. 3a)
. This represents from 21 to 48% of the irrigation water applied during this 30-min interval. The lowest runoff depths were observed on the MP treatments and the highest overall runoff depths were obtained from the NI treatments. Manure treatments did not have a significant effect on runoff depth. Runoff depths were quite variable from year to year and among the tillage and manure treatments.
|
|
Direct comparisons with hydrological data from previous runoff studies are complicated by differences in site hydrology, variation in natural precipitation or in rainfall simulation methodology, variability in tillage methods, and differences in manure characteristics and the timing of incorporation. These factors are particularly important in determining runoff rates and depths, but they can also affect other water quality variables.
Several other studies have reported significantly lower runoff volumes from conventional tillage compared with no-till or reduced-till systems (Mueller et al., 1984a; Daverede et al., 2003). Differences in the volume of runoff among various tillage treatments have generally been attributed to differences in crop residue (Daverede et al., 2003), manure application and incorporation (Mueller et al., 1984a; Bundy et al., 2001), and higher infiltration rates and surface roughness in some cultivated soils (Blevins et al., 1990). In contrast, Seta et al. (1993) reported lower runoff volumes and mean runoff rates from no-till treatments than from chisel-plow or moldboard-plow treatments. The higher infiltration rate in their study for the no-till treatment was attributed to less surface sealing and more undisturbed macropores.
Differences in runoff volumes tend to be most apparent in spring, but lessen with time due to the development of the crop canopy and additional cultivation (Hansen et al., 2000b) or to surface crusting (Mueller et al., 1984a). Since our plots were treated after harvest, results are likely comparable to spring conditions.
Manure tends to increase infiltration; therefore, application of manure before tillage generally resulted in lower runoff volumes compared to the control (Mueller et al., 1984a; Ginting et al., 1998a), as was observed in our study.
Total Suspended Solids
Total suspended solids (TSS) concentrations were greatest from MP treatments and least from the NI treatments (Table 3). However, TSS loads were not significantly different among the tillage treatments (Table 3). Manure treatments significantly reduced both TSS loads and concentrations compared to the control treatments. The TSS concentrations and loads were significantly lower in 2000 than in subsequent years (Table 3).
|
Phosphorus
Dissolved Reactive Phosphorus
Tillage methods were evaluated separately for each manure treatment because the tillage by manure interaction was significant (Table 3, Fig. 4)
. For both the control and manure treatments, the DRP FWMC was highest for the NI treatments, but there were no significant differences observed among the control tillage treatments. The DRP FWMC for the MPM treatment was significantly lower than the other three manure treatments, and was similar to the control treatments. Mass loads of DRP showed a similar pattern to the DRP FWMC results, except that a significant year effect was observed, with DRP loads significantly greater in 2002 than in 2001 (Table 3).
|
In general, incorporation of crop residue and manure reduced losses of DRP compared to NI treatments; however, there were no significant differences among the most common tillage methods, CU and DD, and the NI treatments. Most previous studies have compared losses between no-incorporation (unincorporated manure) or reduced tillage (usually ridge tillage) and conventional tillage (chisel or moldboard plow), with mixed results. Most of these studies reported higher DRP losses from no-incorporation or reduced-tillage methods (Mueller et al., 1984b; Ginting et al., 1998b; Eghball and Gilley, 1999; Bundy et al., 2001; Zhao et al., 2001; Daverede et al., 2004); nevertheless, Hansen et al. (2000a) reported the opposite trend. Higher DRP losses under no-till or reduced-tillage systems have been attributed to leaching from crop residues, retention of snow, and the buildup of phosphorus near the soil surface in these treatments. We observed the same trends among the NIM treatments and the other incorporation methods; however, given that average TP and STP levels in the uppermost soils were very similar among the DD and CU treatments (Fig. 1d, Fig. 5) , it is not surprising that DRP losses were comparable among the DD and CU treatments.
|
Total Phosphorus
A significant interaction was observed between tillage and manure treatments so control and manure treatments were analyzed separately for TP FWMC (Table 3, Fig. 4c). The TP FWMC was greatest from the NIM treatment, followed by the CUM, DDM, and MPM treatments. The TP FWMC was significantly greater in 2002 than in 2000 and 2001 (Table 3).
The NI treatments had the greatest TP loads followed by the DD and CU treatments (Table 3). The TP losses for the MP treatments were significantly lower than the other three tillage treatments. Manure treatments had significantly higher TP losses than the control treatments. The TP loads were significantly greater in 2002 than in 2000 and 2001. An increase in 2002 may be attributed to a combination of the accumulation of TP near the soil surface with time in most of the manure treatments (Fig. 1d) and the wetter conditions in 2002.
Bundy et al. (2001) reported that no-till and unincorporated manure applications generally reduced TP loads in runoff compared to unamended soils or incorporated manure due to lower sediment losses. In contrast, numerous studies report higher TP losses from conventional tillage following manure application, due to higher sediment losses than from no-till or reduced-tillage systems (Ginting et al., 1998b; Eghball and Gilley, 1999; Hansen et al., 2000a). In our study, TSS concentrations were greatest from the MP treatments; nevertheless, the MPM treatment had lower TP losses than other treatments due to the lower runoff volumes. Similarly, CUM and DDM treatments had lower TP losses than the NIM treatment, despite higher sediment concentrations. Although TP losses in runoff were lower in the MPM treatment, accumulation of phosphorus within the root zone (data not shown) and increased sediment losses make this incorporation method unsuitable for most soils.
Most previous studies on tillage methods have been conducted during a single year (Eghball and Gilley, 1999; Bundy et al., 2001). Our results suggest that there may be considerable year-to-year variability in phosphorus losses from different incorporation methods. These differences are likely due to a combination of factors, including climatic factors, variability in manure phosphorus concentration and application rates, and the inherent variability in generating runoff using a rainfall simulator (Mueller et al., 1984b).
Relationships between STP and the log (TP FWMC + 1) of both manure and control treatments were not significant in this study. Total P concentrations in runoff from some soils are directly related to sediment concentrations in runoff (Cox and Hendricks, 2000; Aase et al., 2001; Andraski and Bundy, 2003; Kleinman et al., 2004). In our study, a weak but significant linear relationship was only detected between TSS and log(TP FWMC + 1) for the control treatments (Fig. 6) .
|
|
For total NH4N loads, both tillage and manure had significant effects (Table 4). Losses of NH4N were greatest from the NI treatments and both the NI and DD treatment losses were significantly greater than the MP treatments due to higher runoff volumes. Ammonium N losses were also significantly greater from manure treatments than from the control treatments, which can be attributed to higher concentrations of NH4N from manure treatments.
Eghball and Gilley (1999) also reported significantly higher NH4N losses from no-till treatments than from disked treatments during a field rainfall simulation experiment, with generally higher NH4N losses from compost, manure, and fertilizer treatments than from control treatments. Conversely, in a 3-yr monitoring study in two watersheds in Maryland, Angle et al. (1984) observed significantly greater NH4N losses from a conventional-tilled watershed than from a no-till watershed.
Nitrate N
The NO3N FWMCs were significantly affected by tillage and manure (Table 4). Nitrate N FWMCs were highest from the DD treatments and were significantly lower from the MP treatments. Manure treatments had NO3N FWMCs that were significantly greater than the control treatments. None of the NO3N FWMCs exceeded the drinking water guideline of 10 mg L1 in the surface runoff.
As in our study, Eghball and Gilley (1999) reported greater NO3N concentrations from disked plots compared to no-till plots, which they attributed to high concentrations of NO3N in the surface soil and the disturbance from tillage. Eghball and Gilley (1999) also reported much higher values for NO3N concentrations in runoff from manured plots (average: 22.426.6 mg L1); however, the source water used in their study averaged 21 to 23 mg NO3N L1. Nitrate N FWMCs in our study were similar to those reported in a laboratory rainfall simulation of well-mixed, manured soils from Alberta (Wright et al., 2003).
The DD treatments also had the highest mass losses of NO3N, followed by the NI, CU, and MP treatments (Table 4). The NO3N losses from the MP treatments were significantly lower than the other treatments. Nitrate N losses were significantly greater from the manure treatments than from the control treatments.
Zhao et al. (2001) also reported lower NO3N losses in surface runoff from MPM treatments compared with other tillage methods (ridge-tillage); however, losses of NO3N in surface runoff accounted for only a small portion of total NO3N losses. No differences in NO3N losses were detected between disked and no-till treatments in an initial dry run, but greater losses were observed from disked plots than no-till plots under saturated conditions (Eghball and Gilley, 1999).
Total Nitrogen
Total N FWMCs were greatest from the DD treatments, followed by the CU, NI, and MP treatments; however, differences were only significant for the DD and MP treatments (Table 4). Total N FWMCs in runoff were significantly greater in 2000 and 2001 than in 2002.
Total N losses were significantly lower from the MP treatments than from the CU, DD, and NI treatments (Table 4) and were greatest from the DD treatments. Manure treatments had significantly greater TN losses than the control treatments in 2002 (data not shown). Eghball and Gilley (1999) reported no significant differences in TN losses between disked and no-till treatments from an initial dry run; however, losses were significantly higher from disked plots in a second run on saturated plots.
Phosphorus and Nitrogen Leaching
Nearly all simulation tests produced subsurface leachate volumes that exceeded the capacity of the lysimeters, except for a few occasions on NIC plots where runoff occurred very rapidly and no sample was collected. Tillage had no significant effect on subsurface DRP or TP concentrations at a depth of 60 cm (Table 5). The DRP concentrations were 21 to 30% of the TP concentrations for the tillage treatments. Manure treatments had significantly higher subsurface DRP and TP concentrations than the control treatments (Table 5). The mean concentration of TP as DRP was about 6% for the control treatments and 37% for the manure treatments. Subsurface TP concentrations were significantly higher in 2002 than in 2001 and tended to be greater from NI plots. This finding suggests that leaching of P may increase with time after repeated manure applications and that undisturbed soil macropores may allow greater leaching of phosphorus, especially particulate forms (Beauchemin et al., 1998). Soil moisture conditions were also much wetter in 2002 than in 2001 due to above normal precipitation in 2002.
|
Regardless of manure application, NH4N and NO3N concentrations were greatest in subsurface runoff from MP treatments. Manure treatments also had significantly greater subsurface NH4N and NO3N concentrations than control treatments. The differences were only significant for NO3N concentrations in 2002, a wet year.
The CU, DD, and MP treatments had significantly higher subsurface TN concentrations in 2001 than the NI treatments (Fig. 7) . In 2002, the MP treatments had significantly greater TN concentrations than the other tillage treatments. Manure treatments had significantly higher subsurface TN concentrations than the control treatments (Table 5).
|
| CONCLUSIONS |
|---|
|
|
|---|
Tillage had no significant effect on subsurface DRP or TP concentrations, but manure treatments had significantly higher subsurface DRP and TP concentrations than control treatments. Tillage effects on subsurface NH4N, NO3N, and TN concentrations were not consistent among tillage treatments and between years, but MP treatments had significantly higher NO3N and TN concentrations than the other tillage treatments in 2002. Manure treatments generally contained significantly higher subsurface NH4N, NO3N, and TN concentrations than the control treatments, except for subsurface NO3N concentration in 2001.
Cultivated soils are susceptible to sediment and nutrient losses following land application of livestock manure. Based on our study, partial incorporation of manure with a cultivator or double disk was the most favorable option for minimization of combined surface and subsurface nutrient losses during rainfall simulations immediately after annual applications of beef cattle manure in a silage barley cropping system.
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
D. M. Butler, D. H. Franklin, M. L. Cabrera, A. S. Tasistro, K. Xia, and L. T. West Evaluating Aeration Techniques for Decreasing Phosphorus Export from Grasslands Receiving Manure J. Environ. Qual., May 1, 2008; 37(3): 1279 - 1287. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Franklin, C. Truman, T. Potter, D. Bosch, T. Strickland, and C. Bednarz Nitrogen and Phosphorus Runoff Losses from Variable and Constant Intensity Rainfall Simulations on Loamy Sand under Conventional and Strip Tillage Systems J. Environ. Qual., May 7, 2007; 36(3): 846 - 854. [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] |
||||
![]() |
J. J. Miller, E. C. S. Olson, D. S. Chanasyk, B. W. Beasley, F. J. Larney, and B. M. Olson Phosphorus and nitrogen in rainfall simulation runoff after fresh and composted beef cattle manure application. J. Environ. Qual., July 1, 2006; 35(4): 1279 - 1290. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||