Published in J. Environ. Qual. 33:107-113 (2004).
© ASA, CSSA, SSSA
677 S. Segoe Rd., Madison, WI 53711 USA
TECHNICAL REPORTS
Ecological Risk Assessment
An Environmental Threshold for Degree of Phosphorus Saturation in Sandy Soils
V. D. Nair*,
K. M. Portier,
D. A. Graetz and
M. L. Walker
Statistics Department, 522 McCarty Hall, University of Florida, Institute of Food and Agricultural Sciences, Gainesville, FL 32611-0510
* Corresponding author (vdna{at}ifas.ufl.edu).
Received for publication March 14, 2003.
 |
ABSTRACT
|
|---|
There is critical need for a practical indicator to assess the potential for phosphorus (P) movement from a given site to surface waters, either via surface runoff or subsurface drainage. The degree of phosphorus saturation (DPS), which relates a measure of P already adsorbed by a soil to its P adsorption capacity, could be a good indicator of that soil's P release capability. Our primary objective was to find a suitable analytical protocol for determining DPS and to examine the possibility of defining a threshold DPS value for Florida's sandy soils. Four farmer-owned dairy sprayfields were selected within the Suwannee River basin and soil profiles were randomly obtained from each site, as well as from adjacent unimpacted sites. The soil samples were divided either by horizon or depth, and DPS was determined for each soil sample using ammonium-oxalate (DPSOx), Mehlich-1 (DPSM1), and Mehlich-3 (DPSM3) extracts. All methods of DPS calculations were linearly related to one another (r2 > 0.94). Relationships between water-soluble P and DPS indicate that the respective change points are: DPSOx = 20%, DPSM1 = 20%, and DPSM3 = 16%. These relationships include samples from Ap, E, and Bt horizons, and various combinations thereof, suggesting that DPS values can be used as predictors of P loss from a soil irrespective of the depth of the soil within a profile. Taking into consideration the change points, confidence intervals, agronomic soil test values, and DPS values from other studies, we suggest replacing Mehlich-1 P values in the Florida P Index with the three DPS categories (DPSM1 = <30, 3060, and >60%) to assign different P loss ratings in the P Index.
Abbreviations: DPS, degree of phosphorus saturation DPSM1, degree of phosphorus saturation calculated as [M1-P/0.5(M1-Fe + M1-Al)] x 100 DPSM3, degree of phosphorus saturation calculated as [M3-P/0.5(M3-Fe + M3-Al)] x 100 DPSOx, degree of phosphorus saturation calculated as [Ox-P/0.5(Ox-Fe + Ox-Al)] x 100 M1-Al, M1-Fe, and M1-P, Mehlich 1extractable aluminum, iron, and phosphorus, respectively M3-Al, M3-Fe, and M3-P, Mehlich 3extractable aluminum, iron, and phosphorus, respectively Ox-Al, Ox-Fe, and Ox-P, oxalate-extractable aluminum, iron, and phosphorus, respectively STP, soil test phosphorus WSP, water-soluble phosphorus
 |
INTRODUCTION
|
|---|
SANDY SOILS in the Suwannee River basin of northern Florida have little ability to adsorb P, and yet many dairies in the Suwannee River basin routinely apply P-rich lagoon effluent onto permanent sprayfields for waste disposal and nutrient recycling. Increased P loading to these sites may lead to P loss through runoff and subsurface drainage, contributing to surface water quality degradation. Improved P management of these fields requires the evaluation of soil P concentrations relative to the soil's ability to adsorb P.
Recent studies suggest that the DPS, which relates ammonium oxalateextractable P to the sum of oxalate-extractable Fe and Al (DPSOx), is a good indicator of a soil's potential to release P (Hooda et al., 2000). This concept was first introduced in the Netherlands, where it has been shown that P concentrations in the soil solution can exceed a critical concentration well before the soil is completely saturated with P (Breeuwsma and Silva, 1992). This critical concentration is determined by local conditions and generally reflects local surface water criteria for P. The Netherlands has established a water quality goal for ground water of 0.15 mg total P L1, and their studies have shown that leaching of P could occur especially from manure-contaminated soils (Breeuwsma et al., 1995). In the Netherlands, soils with DPSOx of >25% were identified as contributing to ground water pollution with P (Breeuwsma et al., 1995). They calculated DPSOx as [(Ox-P)/0.5(Ox-Fe + Ox-Al)] x 100 where P, Fe, and Al were measured in an oxalate extract.
Oxalate extraction is not frequently performed in soil test laboratories in Florida (Nair and Graetz, 2002) or in other parts of the USA (Sims et al., 2002) due to practical difficulties in the measurement of parameters in the DPS calculations. More common soil tests include Mehlich-1 and Mehlich-3 extractions. The use of these routine agronomic soil tests to calculate DPS would simplify the measurement of DPS, and provide a more accessible analytical tool for P management.
The objectives of this study were to (i) calculate the DPS for manure-impacted and unimpacted sandy soils using ammonium-oxalate (DPSOx), Mehlich-1 (DPSM1), and Mehlich-3 (DPSM3) extractions; (ii) determine the relationship between DPSOx and DPSM1 and DPSM3; (iii) evaluate the relationship between water-soluble phosphorus (WSP) (assessed using either deionized water or 0.01 M CaCl2) and each method of DPS calculation; and (iv) examine the possibility of defining a threshold DPS value for Florida's sandy soils.
 |
MATERIALS AND METHODS
|
|---|
Study Site
The Suwannee River basin was selected for this study. Many of the dairies found in the middle Suwannee River basin, proximal to the Suwannee and Santa Fe Rivers, are situated atop a geomorphic zone classified as the Chiefland Limestone Plain. The upper surface of the aquifer system is relatively close to the surface and layers above it are thin and unconfined (Andrews, 1992). The surficial aquifer system is largely recharged by rainfall that percolates downward through the loose surficial clastic sediments. Water naturally discharges from the aquifer through evaporation, transpiration, spring flow, and downward seepage into the underlying Floridan aquifer system. Phosphorus, nitrogen, and other surface-applied material with leaching potential can move vertically through the soil profile and then both vertically and laterally in the surficial aquifer system.
The middle Suwannee River basin is approximately 25 km long by 25 km wide at the widest point, covering an area of about 1.3 million hectares. The dominant soils of the basin are Entisols, such as Penney (thermic, uncoated Typic Quartzipsamments), Kershaw (thermic, uncoated Typic Quartzipsamments), Ortega (thermic, uncoated Typic Quartzipsamments), or Ridgewood (thermic, uncoated Aquic Quartzipsamments); Ultisols such as Blanton (loamy, siliceous, thermic Grossarenic Paleudults); or Alfisols such as Otela (loamy, siliceous, thermic Grossarenic Paleudalfs) (Soil Survey Staff, 1999). Cropping systems within the sprayfields include rotations of corn (Zea mays L.)perennial peanut (Arachis hypogaea L.)rye (Secale cereale L.), bermudagrass [Cynodon dactylon (L.) Pers.]rye, perennial peanutrye, cornbermudagrassrye, cornsorghum [Sorghum bicolor (L.) Moench]rye, as well as sole crops such as bermudagrass or ryegrass (Lolium perenne L.).
Soil Sampling
Soil samples from each horizon within a 0- to 2-m profile were collected from manure-impacted sites on four farmer-owned dairy sprayfields and from adjacent unimpacted sites in the middle Suwannee River basin using a 5-cm-diameter auger. A total of 57 soil profiles were collected from the manure-impacted sites, and 12 profiles were collected from adjacent unimpacted sites of similar soil types. The soil profile locations (a minimum of three each) were selected to represent different vegetation types and management practices (such as irrigation levels) within the manure-impacted sites. The soil samples were collected either by depth or by horizon, depending on the nature of the soil profile. If the depth of any horizon was >25 cm, then that horizon was subdivided and two or more samples were obtained from the horizon. All samples were air-dried and passed through a 2-mm sieve.
For Dairy 1, the depths sampled were: 1 = 0 to 36 cm, 2 = 36 to 51 cm, 3 = 51 to 71 cm, 4 = 71 to 97 cm, and 5 = 97 to 122 cm. Dairies 2, 3, and 4 were sampled by horizon, with each horizon being associated with a different depth increment. For these three dairies, the sampling depth was at least 2 m and, whenever possible, sampling included part of the underlying Bt horizon. Some typical soil profiles were Ap, E1, E2, E3, E4, and E5 (no Bt within the surface 2-m depth); or Ap, E1, E2, E3, E/Bt, and Bt or various combinations of these horizons up to 2 m in depth.
Soil Characterization
Soil pH was determined using a 1:2 soil and water suspension. Oxalate-extractable Al (Ox-Al), Fe (Ox-Fe), and P (Ox-P) were determined by extraction with 0.1 M oxalic acid + 0.175 M ammonium oxalate (pH = 3.0) (McKeague and Day, 1966). The suspension was equilibrated for 4 h in the dark with continuous shaking, centrifuged, filtered through a 0.45-µm filter, and analyzed for Al, Fe, and P. Mehlich 1, or double acidextractable (0.0125 M H2SO4 + 0.05 M HCl) P (M1-P), Fe (M1-Fe), and Al (M1-Al) were obtained using a 1:4 soil to double acid ratio (Mehlich, 1953). Mehlich-3 extractions for determination of P (M3-P), Fe (M3-Fe), and Al (M3-Al) were performed as proposed by Mehlich (1984). All metals and P in the Mehlich-1 and oxalate solutions were determined using inductively coupled argon plasma spectroscopy (Thermo Jarrel Ash ICAP 61E; Thermo Elemental, Franklin, MA).
Water-soluble P was determined by extracting each soil sample with water at a 1:10 soil to water ratio for 1 h, and determining P on the filtrate collected after passing through a 0.45-µm filter. The CaCl2extractable P was also measured using the suggested method for animal manure (Self-Davis et al., 2000), using a 1:10 soil to 0.01 M CaCl2 solution. Total P was determined by ashing 1.0 g of soil for 2 h at 823 K, and then solubilizing with 6 M HCl (Anderson, 1976). Water-soluble P and total P concentrations were determined by an autoanalyzer (USEPA, 1983; Method 365-1) by the Murphy and Riley (1962) procedure. Total C and N contents of the air-dried samples were determined by an automated combustion procedure using a CNS Analyzer (Carlo Erba, Milan, Italy).
Calculation of the Degree of Phosphorus Saturation
The following methods of calculation were adopted:
where
is an empirical factor that compares different soils with respect to P saturation. The value of
for the current studies was taken as 0.50 (Beauchemin and Simard, 1999; Breeuwsma and Silva, 1992; Koopmans et al., 2003; Schoumans, 2000; Sims et al., 2002). The
value was close to the value of 0.55 for Spodosols in Florida (Nair and Graetz, 2002).
Statistical Analysis
Mean concentrations were computed for each variable by dairy and impact status. Concentrations for all horizons below the surface horizon were averaged for Dairies 2 through 4 and unimpacted sites to produce a subsurface concentration. Comparisons between surface and subsurface average concentrations were made using general linear models for impacted and unimpacted sites separately due to large differences in residual variances. Comparisons between impacted and unimpacted average concentrations were used to simply indicate trends and hence two-sample t tests assuming unequal variances were used.
The relationship between DPS and WSP was modeled as a segmented line (Eq. [1]), with parameters estimated using nonlinear least squares. The change point (d0) in the fitted segmented-line model was directly estimated. To ensure that the two line segments joined at the change point, the slope of the left-hand line is estimated as a function of the change point and other model parameters (Eq. [2]). Standard errors were estimated from the Fisher information matrix and confidence intervals are constructed using these standard errors and an appropriate t distribution critical value. Computations were performed in SAS (SAS Institute, 2001) using a NLIN procedure.
 | [1] |
 | [2] |
 |
RESULTS AND DISCUSSION
|
|---|
Soil Characterization
Texture analysis of selected soil samples (n = 37) representative of all horizons gave mean values of 96% sand, 2% silt, and 2% clay (data not shown). Some chemical properties of the soils used for this study are presented in Table 1. For Dairy 1, the various chemical properties for the subsurface soils were averaged by depth since the depths for all soil profiles were identical. For the other three dairies, mean values of chemical properties (Table 1) for all horizons below Ap were considered together although the samples were analyzed separately. However, extraction values from individual soil samples were used in the calculation of threshold DPS values. Our intention was not to evaluate depth distribution of DPS in these soil profiles, but to determine whether DPS in subsurface samples (irrespective of horizon type) could be related to WSP as well.
View this table:
[in this window]
[in a new window]
|
Table 1. Mean values for selected chemical characteristics of the soils from four dairy sprayfields and two unimpacted sites.
|
|
The sum of Ox-Al and Ox-Fe provides an indicator of the sorption capacity of an acid soil (Breeuwsma and Silva, 1992; Nair et al., 1998). The surface soils in this study had Ox-Al + Ox-Fe values ranging from 15 to 40 mmol kg1 compared with the range of values for Dutch agricultural soils of 50 to 110 mmol kg1 (Schoumans and Groenendijk, 2000). For 465 soils in Delaware, Sims et al. (2002) reported mean Ox-Al + Ox-Fe values of 52.5 mmol kg1 with a median value of 37.6 mmol kg1. Thus, the surface sandy soils of the Suwannee River basin would probably have poorer adsorbing capacities than the sandy soils of the Netherlands and comparable or poorer adsorbing capacity than the soils of the Mid-Atlantic coastal plain. Florida soils were formed from coastal plain (mainly marine) pre-weathered sediments dominated by quartz sand and low in sources of Fe and Al (Brown et al., 1990).
The pH values of manure-impacted soils were invariably higher than for unimpacted soils, with high Ca and/or Mg concentrations being typical of dairy manureimpacted soils (Nair et al., 1995). Calcium and Mg concentrations were significantly higher in the subsurface horizons compared with the surface horizons for impacted soils (P < 0.01 using all data after accounting for site differences in overall mean concentrations) suggesting manure constituent movement through the soil profile. Total C concentrations in the soil samples were variable, with a tendency toward higher values in the manure-impacted compared with the unimpacted soils at the surface (P = 0.04) but less so at the subsurface (P = 0.055) (Table 1). Total N concentrations were below detection limits at several of the sites (Table 1).
Mean values for all P parameters (WSP, M1-P, M3-P, and Ox-P; Table 2) indicate higher concentrations in the surface horizons compared with the subsurface horizons (P < 0.01 for all data adjusted for overall site average levels). The concentrations for a given soil vary as M3-P > Ox-P > M1-P > WSP. Concentrations of WSP, M1-P, M3-P, and Ox-P for Dairy 1 all decrease with depth (P < 0.001, trend confirmed using Waller LSD separation). Mehlich-1 P concentration, which is Florida's soil test phosphorus (STP), shows a concentration of 30 mg kg1 at a 51- to 71-cm depth, indicating that the soil has P concentrations at this depth that are above the agronomic critical level (Kidder et al., 2002).
View this table:
[in this window]
[in a new window]
|
Table 2. Mean values for water-soluble phosphorus (WSP), Mehlich 1extractable phosphorus (M1-P), Mehlich 3extractable phosphorus (M3-P), and oxalate-extractable phosphorus (Ox-P) for the soils from four dairy sprayfields and two unimpacted sites.
|
|
Analytical Protocol for Degree of Phosphorus Saturation Calculations
Relationships among the various methods of DPS calculations were linear with r2 values of >0.94 (Table 3). These relationships suggest that Mehlich-1 and Mehlich-3 solutions can be used for calculation of DPS in states where they are used routinely for agronomic soil tests. The DPSOx involves measurements of P, Fe, and Al in an oxalate solution. In Florida, and in many parts of the USA, institutional and private laboratories may not always have the facilities to measure these parameters. The oxalate extraction method requires extraction in the dark and measurements of the elements in the solution within a week (Schoumans, 2000), making it difficult to perform the analyses on a routine basis. On the other hand, Mehlich 1 is the current soil test P in Florida and conducting a couple of additional analyses in the extract such as Fe and Al is not a major problem. Further, the time involved in Mehlich-1 extraction is short (5 min) compared with the 4-h extraction time for oxalate extractions. Mehlich 3 is widely used as a soil test P in the Mid-Atlantic states (Sims et al., 2002), so we included calculations of DPS from Mehlich-3 analyses. Knowing the relationships between the three methods of DPS calculations would allow DPS data conversion of comparable soils (e.g., sandy Entisols, Ultisols, or Alfisols) from DPS determinations conducted in most analytical laboratories (soil test solutions) to information available from research laboratories (oxalate solutions).
View this table:
[in this window]
[in a new window]
|
Table 3. Relationships among the various methods of degree of phosphorus saturation (DPS) calculations (DPSM1 and DPSOx, DPSM3 and DPSOx, DPSM1 and DPSM3).
|
|
Relationship between Water-Soluble Phosphorus (Deionized Water or 0.01 M CaCl2) and Degree of Phosphorus Saturation
Relationships between water-soluble P (deionized water) and DPS, calculated as DPSOx (Fig. 1)
, DPSM1 (Fig. 2)
, and DPSM3 (Fig. 3)
, each gave a "change point." A change point DPS value may be defined as that value above which there is a rapid increase in WSP and therefore a likelihood of a negative impact of P in the soil on water quality. For DPSOx, the change point is at 20% (95% confidence limits: 1724%); for DPSM1, the change point is at 20% (95% confidence limits: 1029%); and for DPSM3, the change point is at 16% (95% confidence limits: 1121%). Parameter estimates including standard errors and R2 values for the fitted nonlinear relationship models are given in Table 4. Soil samples from all soil horizons, including the upper part of the Bt horizon, were included in the calculations to determine the relationships. Separating the surface and subsurface soils to determine the change point gave almost identical change points by the three methods of DPS determinations. Thus, a measure of DPS can be used to predict the potential for P release from surface soils or from soils at any depth within a profile. All surface soils in the manure-impacted sprayfields have WSP concentrations greater than the change points calculated using any one of the three methods of DPS calculations. Several of the subsurface soils are also above the change points (Fig. 1, 2, and 3), suggesting P movement throughout the soil profile. Mean DPS values for unimpacted soils for all horizons were DPSOx = 18%, DPSM1 = 18%, and DPSM3 = 16%, all at or below the respective change points.

View larger version (21K):
[in this window]
[in a new window]
|
Fig. 1. Relationship between the concentration of water-soluble phosphorus (WSP) and the degree of P saturation calculated using an oxalate extraction (DPSOx) for manure-impacted surface and subsurface soils from the middle Suwannee River basin.
|
|

View larger version (22K):
[in this window]
[in a new window]
|
Fig. 2. Relationship between the concentration of water-soluble phosphorus (WSP) and the degree of P saturation calculated using a Mehlich-1 extraction (DPSM1) for manure-impacted surface and subsurface soils from the middle Suwannee River basin.
|
|

View larger version (23K):
[in this window]
[in a new window]
|
Fig. 3. Relationship between the concentration of water-soluble phosphorus (WSP) and the degree of P saturation calculated using a Mehlich-3 extraction (DPSM3) for manure-impacted surface and subsurface soils from the middle Suwannee River basin.
|
|
View this table:
[in this window]
[in a new window]
|
Table 4. Parameter estimates (with standard errors in parentheses) and R2 values for the fitted nonlinear relationship models.
|
|
The DPSOx is closely related to P concentrations in leachate waters (Leinweber et al., 1999; Maguire and Sims, 2002), suggesting that DPSOx can be a suitable tool for predicting subsurface P losses. Nair and Graetz (2002) showed that DPSM1 can be used as an indicator of soluble P for both surface A horizons and subsurface Bh (spodic) horizons of Spodosols from the Lake Okeechobee basin. Soils with DPSOx of >25% contributed to ground water pollution by P in the Netherlands (Breeuwsma et al., 1995). The 25% value corresponds to the Netherlands water quality goal of 0.15 mg total P L1 for ground water P concentrations. No similar comparisons are available for our studies. Values for DPSOx of >30% in topsoils have been shown to be a threat to water quality degradation in Mid-Atlantic U.S. soils (Paulter and Sims, 2000), and also to be associated with P losses in runoff (Pote et al., 1996).
Maguire and Sims (2002) defined the Mehlich-3 P saturation ratio (M3-PSR) as the ratio between Mehlich-3 P and the sum of Mehlich 3extractable Fe and Al. Conversion of M3-PSR in their studies (range of 0.100.15) to DPSM3 (expressed as a percentage) requires introducing a factor of 200 into their equation. The resulting DPSM3 is 20 to 30%. The range slightly exceeds our DPSM3 change point of 16% (95% confidence limits: 1121%). The range of DPSM3 calculated for the soils used by Maguire and Sims (2002) corresponds to 30 to 45% DPSM1 using the conversion equation in Table 3. This information was also used in our selection of threshold ranges for Florida soils.
We also calculated change points using 0.01 M CaCl2 solution (Table 5) because CaCl2 extraction has been shown to be a useful P leaching indicator (McDowell and Sharpley, 2001). Change points were detected (Table 5), but they were higher than those values obtained using deionized water. The range among DPS values calculated using the different methods was greater using the CaCl2 extraction (2638%) compared with the water extraction (1620%). Maguire and Sims (2002) and McDowell and Sharpley (2001) reported that CaCl2P concentrations were generally less than WSP concentrations. Mean CaCl2P concentrations in our studies were only one-third of the mean WSP concentrations.
View this table:
[in this window]
[in a new window]
|
Table 5. Change points and their confidence intervals for DPSOx, DPSM1, and DPSM3 using 0.01 M CaCl2P instead of water-soluble P.
|
|
The Agronomic Soil Test Factor for Setting Environmental Phosphorus Limits
Mehlich-1 P is the agronomic soil test P currently used in Florida, and it is appropriate to relate the soil test to environmental parameters such as the DPS. In Florida, Mehlich-1 P values above 30 mg kg1 are considered high from an agronomic standpoint and a value above 60 mg kg1 is considered very high (Kidder et al., 2002). The University of Delaware has rated soils with Mehlich-1 P values of >50 mg P kg1 as excessive (Paulter and Sims, 2000). We therefore examined corresponding DPSOx values for these Mehlich-1 P concentrations. A Mehlich-1 P concentration of 30 mg P kg1 corresponds to a DPSOx value of 22%, whereas a 60 mg P kg1 value corresponds to a DPSOx value of 28% (Fig. 4)
. These values agree well with the DPSOx value of 25% corresponding to 50 mg P kg1 (Paulter and Sims, 2000). Paulter and Sims (2000) used an
value of 0.68 in their equation to calculate DPS whereas we used a value of 0.50. Given the empirical nature of
, we included the factor in our calculations, primarily for comparison with DPSOx values that have been used in recent literature. The
value for Spodosols of the Lake Okeechobee basin in Florida is 0.55 (Nair and Graetz, 2002), which is close to the 0.5 value used in the current studies.

View larger version (24K):
[in this window]
[in a new window]
|
Fig. 4. Relationship between Mehlich-1 P and the degree of P saturation calculated using an oxalate extraction (DPSOx) for manure-impacted surface and subsurface soils from the middle Suwannee River basin. Horizontal dotted lines indicate agronomic high (30 mg P kg1) and very high (60 mg P kg1) soil test P concentrations in Florida (Kidder et al., 2002).
|
|
A Threshold Degree of Phosphorus Saturation Value for Florida's Sandy Soils
Based on the change points, confidence intervals, and agronomic soil test values measured herein, we recommend a threshold DPSM1 of 30% for Florida sands. Values for DPSM1 of 31 to 60% warrant caution with regard to further addition of P to a land-use system, and DPSM1 values of >60% suggest soils as contributors to water quality impairment. While these ranges are arbitrary, they are useful as indicators of P loss potential from agricultural systems. The suggested ranges may be subjected to changes based on added field and laboratory information. The recommended threshold value of 30% is above agronomic fertility requirements (Fig. 4), and further P additions to the soils are unnecessary for plant growth. This threshold value of 30% is also comparable with the Netherlands-recommended DPS value of 25% and the calculated DPSM1 range between 30 and 45% for the soils Maguire and Sims (2002) used in their studies. Dairy soils from the Suwannee River basin, introduced as a check of our data, fitted the WSPDPS relationship (Fig. 3) (Nair et al., 2002b). Also, the relationship held for soils of the Okeechobee basin of Florida, in spite of the extremely poor retention capacity of these soils. Therefore, the threshold DPS ranges recommended based on the best professional judgment of the authors were applicable to the Okeechobee soils as well (Nair et al., 2002a).
Maguire and Sims (2002) concluded that an environmental soil limit (such as DPSM1) should be used as an initial indicator of potential water quality problems, with a more detailed scheme, such as the P Index, being used to access the risk of P contamination at any given site. The Florida P Index (http://efotg.nrcs.usda.gov/popmenu3FS.aspx?Fips=12001&MenuName=menuFL.zip; select "Section IV" then "B. Tools" then "Florida Phosphorus Index"; verified 26 June 2003) was developed as a field-based index to assess site conditions and potential P loss vulnerability. The index includes consideration of transport factors such as soil erosion, soil runoff class, leaching potential, and distance from a water course, along with management factors such as STP, P application method, and source and rate of P application. Agronomists often consider STP as an inappropriate factor for evaluating environmental P losses as STP was originally calibrated for agronomic purposes (Sharpley et al., 1999). DeLaune et al. (2002) recently showed that STP is not the most reliable indicator of P in runoff when animal manure (poultry manure) was surface applied to a plot. The relationship between WSP and M1-P was linear in our study, WSP = 0.612(M1-P) + 0.3976; r2 = 0.9072; n = 405; p < 0.001. Also, WSP was related to Mehlich-3 P as WSP = 0.038(M3-P) 0.6606; r2 = 0.8720; n = 405; p < 0.001 (data not shown). No change point could be identified in the WSP relationship with either Mehlich-1 or Mehlich-3 P.
We recommend replacing the STP factor in the Florida P Index with DPSM1 for the fertility index value. The three ranges for DPSM1 (<30, 3060, and >60%) would then be assigned different P loss ratings. At present we are considering only surface soil DPS values to replace surface Mehlich-1 P concentrations in the Florida P Index.
The current Florida P Index attempts to incorporate leaching potential based primarily on visual observation of the Bt horizon within a soil profile (Nair and Graetz, 2002). This study shows that DPS can be related to WSP for all soil samples throughout a soil profile, including samples of the Bt horizon. However, it may not be practical for the field evaluator to determine DPS at regular intervals throughout a soil profile during evaluation of the P Index. We are currently evaluating the possibility of incorporating subsurface DPS values into the P Index using simple field tests that could be related to DPS throughout a soil profile.
 |
CONCLUSIONS
|
|---|
Strong correlations exist between DPSOx and DPSM1, DPSM3 and DPSM1, and DPSM3 and DPSOx, indicating that the three methods are equally appropriate for DPS calculations. For the sandy soils of Florida, as well as for parts of the USA where Mehlich-1 P is used routinely as the STP, DPSM1 could be a convenient indicator of P loss from an agricultural system. Where Mehlich-3 P is the routine STP, then DPSM3 may be the appropriate indicator. Relationships between WSP and DPS for Florida soils support change points of DPSOx = 20%, DPSM1 = 20%, and DPSM3 = 16%. The relationships include soils from all horizons (Ap, E, Bt, and various combinations thereof), indicating that DPS values can be related to P loss from a soil irrespective of the depth of the soil within a profile. Various factors employed for the calculation of DPS, including both confidence intervals and agronomic factors, suggest that threshold DPS values should be used with caution. However, a comprehensive scheme for potential P loss vulnerability, such as a P Index that includes DPS as a factor, appears appropriate.
 |
ACKNOWLEDGMENTS
|
|---|
We dedicate this manuscript to the memory of Dr. E.C. (Tito) French, 19451999. Sincere thanks are due to Dr. W.G. Harris and Dr. G.A. O'Connor for their comments and suggestions during the various stages of this work. We would also like to thank Dr. K.R. Woodard and Mr. Tony Sweat for their assistance, particularly with soil sampling, and Ms. Dawn Lucas for assistance with laboratory analysis.
 |
NOTES
|
|---|
This research was supported by the Florida Agricultural Experiment Station and a grant from the Florida Department of Environmental Protection, and approved for publication as Journal Series no. R-09384.
 |
REFERENCES
|
|---|
- Anderson, J.M. 1976. An ignition method for determination of total phosphorus in lake sediments. Water Res. 10:329331.
- Andrews, W.J. 1992. Reconnaissance of water quality at nine dairy farms in North Florida, 19901991. Water Resour. Investigations Rep. 92-4058. U.S. Geol. Survey, Reston, VA.
- Beauchemin, S., and R.R. Simard. 1999. Soil phosphorus saturation degree: Review of some indices and their suitability for phosphorus management in Canada. Can. J. Soil Sci. 79:615625.
- Breeuwsma, A., J.G.A. Reijerink, and O.F. Schoumans. 1995. Impact of manure on accumulation and leaching of phosphate in areas of intensive livestock farming. p. 239251. In K. Steele (ed.) Animal waste and the land-water interface. Lewis Publ.-CRC, New York.
- Breeuwsma, A., and S. Silva. 1992. Phosphorus fertilisation and environmental effects in the Netherlands and the Po Region (Italy). Rep. 57. Agric. Res. Dep., The Winand Staring Centre for Integrated Land, Soil, and Water Res., Wageningen, the Netherlands.
- Brown, R.B., E.L. Stone, and V.W. Carlise. 1990. Soils. p. 3969. In R.L. Myers and J.J. Ewel (ed.) Ecosystems of Florida. Univ. of Central Florida Press, Orlando.
- DeLaune, P.B., P.A. Moore, Jr., D.K. Carman, A.N. Sharpley, B.E. Haggard, and T.C. Daniel. 2002. Using the phosphorus index to reduce phosphorus runoff from soils in the Eucha/Spavinaw watershed. Project #445. U.S. Poultry and Egg Assoc., Tucker, GA.
- Hooda, P.S., A.R. Rendell, A.C. Edwards, P.J.A. Withers, M.N. Aitken, and V.W. Truesdale. 2000. Relating soil phosphorus indices to potential phosphorus release to water. J. Environ. Qual. 29:11661171.[Abstract/Free Full Text]
- Kidder, G., C.G. Chamblis, and R. Mylavarapu. 2002. UF/IFAS standard fertilization recommendations for agronomic crops. SL129. Soil & Water Sci., Coop. Ext. Serv., IFAS, Gainesville, FL.
- Koopmans, G.F., W.J. Chardon, J. Dolfing, O. Oenema, P. van der Meer, and W.H. van Riemsdijk. 2003. Wet chemical and phosphorus-31 nuclear magnetic resonance analysis of phosphorus speciation in a sandy soil receiving long-term fertilizer or animal manure applications. J. Environ. Qual. 32:287295.[Abstract/Free Full Text]
- Leinweber, P., R. Meissner, K.U. Eckhardt, and J. Seeger. 1999. Management effects on forms of phosphorus in soil and leaching losses. Eur. J. Soil Sci. 50:413424.
- Maguire, R.O., and J.T. Sims. 2002. Soil testing to predict phosphorus leaching. J. Environ. Qual. 31:16011609.[Abstract/Free Full Text]
- McDowell, R.W., and A.N. Sharpley. 2001. Approximating phosphorus release from soils to surface runoff and subsurface drainage. J. Environ. Qual. 30:508520.[Abstract/Free Full Text]
- McKeague, J.A., and J.H. Day. 1966. Dithionate and oxalate-extractable Fe and Al as aids in differentiating various classes of soils. Can. J. Soil Sci. 46:1322.
- Mehlich, A. 1953. Determination of P, Ca, Mg, K, Na, and NH4. Soil Testing Div. Publ. 1-53. North Carolina Dep. of Agric., Raleigh.
- Mehlich, A. 1984. Mehlich 3 soil test extractant: A modification of Mehlich 2 extractant. Commun. Soil Sci. Plant Anal. 15:14091416.
- Murphy, J., and J.P. Riley. 1962. A modified single solution method for the determination of phosphate in natural waters. Anal. Chim. Acta 27:3136.
- Nair, V.D., and D.A. Graetz. 2002. Phosphorus saturation in Spodosols impacted by manure. J. Environ. Qual. 31:12791285.[Abstract/Free Full Text]
- Nair, V.D., D.A. Graetz, and K.M. Portier. 1995. Forms of phosphorus in soil profiles from dairies of south Florida. Soil Sci. Soc. Am. J. 59:12441249.[Abstract/Free Full Text]
- Nair, V.D., D.A. Graetz, and K.R. Reddy. 1998. Dairy manure influences on phosphorus retention capacity of Spodosols. J. Environ. Qual. 27:522527.[Abstract/Free Full Text]
- Nair, V.D., D.A. Graetz, and M.L. Walker. 2002a. Degree of phosphorus saturation (DPS): An indicator for predicting water quality degradation. Paper 222. In Trans. 17th World Congr. of Soil Sci., Bangkok, Thailand [CD ROM]. 1421 Aug. 2002.
- Nair, V.D., D.A. Graetz, M.L. Walker, and K.M. Portier. 2002b. Threshold degree of phosphorus saturation for Florida's sandy soils. In 2002 Annual meeting abstracts [CD ROM]. ASA, CSSA, and SSSA, Madison, WI.
- Paulter, M.C., and J.T. Sims. 2000. Relationships between soil test phosphorus, soluble phosphorus, and phosphorus saturation in Delaware soils. Soil Sci. Soc. Am. J. 64:765773.[Abstract/Free Full Text]
- Pote, D.H., T.C. Daniel, A.N. Sharpley, P.A. Moore, Jr., D.R. Edwards, and D.J. Nichols. 1996. Relating extractable soil phosphorus to phosphorus losses in runoff. Soil Sci. Soc. Am. J. 60:855859.[Abstract/Free Full Text]
- SAS Institute. 2001. SAS Version 8.1. SAS Inst., Cary, NC.
- Schoumans, O.F. 2000. Determination of the degree of phosphorus saturation in non-calcareous soils. p. 3134. In G.M. Pierzynski (ed.) Methods of phosphorus analysis for soils, sediments, residuals, and waters. Southern Coop. Ser. Bull. 396. North Carolina State Univ., Raleigh.
- Schoumans, O.F., and P. Groenendijk. 2000. Modeling soil phosphorus levels and phosphorus leaching from agricultural land in the Netherlands. J. Environ. Qual. 29:111116.
- Self-Davis, M.L., P.A. Moore, Jr., and B.C. Joern. 2000. Determination of water- and/or dilute salt-extractable phosphorus. p. 2426. In G.M. Pierzynski (ed.) Methods of phosphorus analysis for soils, sediments, residuals, and waters. Southern Coop. Ser. Bull. 396. North Carolina State Univ., Raleigh.
- Sharpley, A.N., T. Daniel, T. Sims, J. Lemunyon, R. Stevens, and R. Parry. 1999. Agricultural phosphorus and eutrophication. ARS-149. USDA Agric. Res. Serv., Washington, DC.
- Sims, J.T., R.O. Maguire, A.B. Leytem, K.L. Gartley, and M.C. Paulter. 2002. Evaluation of Mehlich 3 as an agri-environmental soil phosphorus test for the Mid-Atlantic United States of America. Soil Sci. Soc. Am. J. 66:20162032.[Abstract/Free Full Text]
- Soil Survey Staff. 1999. Soil taxonomy: A basic system of soil classification for making and interpreting soil surveys. 2nd ed. USDA-NRCS Agric. Handb. 436. U.S. Gov. Print. Office, Washington, DC.
- USEPA. 1983. Methods of chemical analysis of water and wastes. Environ. Monitoring and Support Lab., Cincinnati, OH.
This article has been cited by other articles:

|
 |

|
 |
 
A. K. Obour, M. L. Silveira, M.B. Adjei, J. M. Vendramini, and J. E. Rechcigl
Cattle Manure Application Strategies Effects on Bahiagrass Yield, Nutritive Value, and Phosphorus Recovery
Agron. J.,
September 1, 2009;
101(5):
1099 - 1107.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
Y. C. Newman, S. Agyin-Birikorang, M. B. Adjei, J. M. Scholberg, M. L. Silveira, J. M. B. Vendramini, J. E. Rechcigl, and L. E. Sollenberger
Nitrogen Fertilization Effect on Phosphorus Remediation Potential of Three Perennial Warm-Season Forages
Agron. J.,
August 31, 2009;
101(5):
1243 - 1248.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
A. Mukherjee, V. D. Nair, M. W. Clark, and K. R. Reddy
Development of Indices to Predict Phosphorus Release from Wetland Soils
J. Environ. Qual.,
March 25, 2009;
38(3):
878 - 886.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
S. Agyin-Birikorang, O. O. Oladeji, G. A. O'Connor, T. A. Obreza, and J. C. Capece
Efficacy of Drinking-Water Treatment Residual in Controlling Off-Site Phosphorus Losses: A Field Study in Florida
J. Environ. Qual.,
March 25, 2009;
38(3):
1076 - 1085.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
Z. Dou, C. F. Ramberg, J. D. Toth, Y. Wang, A. N. Sharpley, S. E. Boyd, C. R. Chen, D. Williams, and Z. H. Xu
Phosphorus Speciation and Sorption-Desorption Characteristics in Heavily Manured Soils
Soil Sci. Soc. Am. J.,
January 21, 2009;
73(1):
93 - 101.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
J. Siemens, K. Ilg, H. Pagel, and M. Kaupenjohann
Is Colloid-Facilitated Phosphorus Leaching Triggered by Phosphorus Accumulation in Sandy Soils?
J. Environ. Qual.,
October 23, 2008;
37(6):
2100 - 2107.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
L. R.F. Alleoni, S. R. Brinton, and G. A. O'Connor
Runoff and Leachate Losses of Phosphorus in a Sandy Spodosol Amended with Biosolids
J. Environ. Qual.,
January 4, 2008;
37(1):
259 - 265.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
O. O. Oladeji, G. A. O'Connor, J. B. Sartain, and V. D. Nair
Controlled Application Rate of Water Treatment Residual for Agronomic and Environmental Benefits
J. Environ. Qual.,
October 16, 2007;
36(6):
1715 - 1724.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
M. Chrysostome, V. D. Nair, W. G. Harris, and R. D. Rhue
Laboratory Validation of Soil Phosphorus Storage Capacity Predictions for Use in Risk Assessment
Soil Sci. Soc. Am. J.,
August 9, 2007;
71(5):
1564 - 1569.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
J. L. Little, S. C. Nolan, J. P. Casson, and B. M. Olson
Relationships between Soil and Runoff Phosphorus in Small Alberta Watersheds
J. Environ. Qual.,
July 17, 2007;
36(5):
1289 - 1300.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
K. R. Woodard, L. E. Sollenberger, L. A. Sweat, D. A. Graetz, V. D. Nair, S. J. Rymph, L. Walker, and Y. Joo
Phosphorus and Other Soil Components in a Dairy Effluent Sprayfield within the Central Florida Ridge
J. Environ. Qual.,
May 25, 2007;
36(4):
1042 - 1049.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
A. L. Shober and J. T. Sims
Integrating Phosphorus Source and Soil Properties into Risk Assessments for Phosphorus Loss
Soil Sci. Soc. Am. J.,
March 12, 2007;
71(2):
551 - 560.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
J. P. Casson, D. R. Bennett, S. C. Nolan, B. M. Olson, and G. R. Ontkean
Degree of Phosphorus Saturation Thresholds in Manure-Amended Soils of Alberta
J. Environ. Qual.,
October 27, 2006;
35(6):
2212 - 2221.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
C. R. Bond, R. O. Maguire, and J. L. Havlin
Change in Soluble Phosphorus in Soils following Fertilization is Dependent on Initial Mehlich-3 Phosphorus
J. Environ. Qual.,
August 9, 2006;
35(5):
1818 - 1824.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
B. L. Allen and A. P. Mallarino
Relationships between Extractable Soil Phosphorus and Phosphorus Saturation after Long-Term Fertilizer or Manure Application
Soil Sci. Soc. Am. J.,
February 2, 2006;
70(2):
454 - 463.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
D. V. Ige, O. O. Akinremi, and D. N. Flaten
Environmental Index for Estimating the Risk of Phosphorus Loss in Calcareous Soils of Manitoba
J. Environ. Qual.,
October 12, 2005;
34(6):
1944 - 1951.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
N. O. Nelson, J. E. Parsons, and R. L. Mikkelsen
Field-Scale Evaluation of Phosphorus Leaching in Acid Sandy Soils Receiving Swine Waste
J. Environ. Qual.,
October 12, 2005;
34(6):
2024 - 2035.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
M. Mamo, S. C. Gupta, C. J. Rosen, and U. B. Singh
Phosphorus Leaching at Cold Temperatures as Affected by Wastewater Application and Soil Phosphorus Levels
J. Environ. Qual.,
June 7, 2005;
34(4):
1243 - 1250.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
K. Ilg, J. Siemens, and M. Kaupenjohann
Colloidal and Dissolved Phosphorus in Sandy Soils as Affected by Phosphorus Saturation
J. Environ. Qual.,
April 20, 2005;
34(3):
926 - 935.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
P. A. Vadas, P. J. A. Kleinman, A. N. Sharpley, and B. L. Turner
Relating Soil Phosphorus to Dissolved Phosphorus in Runoff: A Single Extraction Coefficient for Water Quality Modeling
J. Environ. Qual.,
March 1, 2005;
34(2):
572 - 580.
[Abstract]
[Full Text]
[PDF]
|
 |
|