Published online 20 April 2005
Published in J Environ Qual 34:951-962 (2005)
DOI: 10.2134/jeq2004.0295
© 2005 American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
TECHNICAL REPORTS
Waste Management
Evaluation of Soils for Use as Liner Materials
A Soil Chemistry Approach
Tom M. DeSuttera,b and
Gary M. Pierzynskia,*
a Department of Agronomy, Throckmorton Hall, Kansas State University, Manhattan, KS 66506
b Now with the USDA-ARS National Soil Tilth Laboratory, Ames, IA 50011
* Corresponding author (gmp{at}ksu.edu)
Received for publication July 30, 2004.
 |
ABSTRACT
|
|---|
Movement of NH4+ below animal waste lagoons is generally a function of the whole-lagoon seepage rate, soil mineralogy, cations in the lagoon liquor, and selectivity for NH4+ on the soil-exchange sites. Binary exchange reactions (Ca2+K+, Ca2+NH4+, and K+NH4+) were conducted on two soils from the Great Plains and with combinations of these soils with bentonite or zeolite added. Binary exchanges were used to predict ternary exchanges
following the RothmundKornfeld approach and GainesThomas convention. Potassium and NH4+ were preferred over Ca2+, and K+ was preferred over NH4+ in all soils and soils with amendments. Generally, the addition of bentonite did not change cation selectivity over the native soils, whereas the addition of zeolite did. The RothmundKornfeld approach worked well for predicting equivalent fractions of cations on the exchanger phase when only ternary-solution phase compositions were known. Actual swine- and cattle-lagoon solution compositions and the RothmundKornfeld approach were used to project that native soils are predicted to retain 53 and 23%, respectively, of the downward-moving NH4+ on their exchange sites. Additions of bentonite or zeolite to soils under swine lagoons may only slightly improve the equivalent fraction of NH4+ on the exchange sites. Although additions of bentonite or zeolite may not help increase the NH4+ selectivity of a liner material, increases in the overall cation exchange capacity (CEC) of a soil will ultimately decrease the amount of soil needed to adsorb downward-moving NH4+.
Abbreviations: CEC, cation exchange capacity EFA, equivalent fraction of ammonium
G, Gibb's free energy of exchange RPD, relative percent difference S1, Soil 1 S1B, Soil 1 + bentonite S1Z, Soil 1 + zeolite S2, Soil 2 S2B, Soil 2 + bentonite S2Z, Soil 2 + zeolite
 |
INTRODUCTION
|
|---|
SOIL-LINED EARTHEN BASINS are often used to store and treat wastes from animal-feeding operations in the Great Plains. These lagoons range in size from 0.5 to 2.5 ha, with a depth of up to 6 m, and generally have compacted soil liners from 0.3 to 0.6 m thick (Ham and DeSutter, 2000). Liners may be constructed using soil present on-site or soil may be brought in from an external source. Bentonite may be added to the liner to reduce hydraulic conductivity. Historically, these soils have only been evaluated on their ability to restrict water movement and, although suggestions have been made that the CEC of the liner material should be included in the liner-evaluation procedure (Barrington and Madramootoo, 1989; Ham and DeSutter, 2000), chemical properties of the liner materials have gone largely unevaluated.
Recent reports have indicated that a substantial amount of NH4N can move into the soil underlying animal waste lagoons (DeSutter et al., unpublished data, 2005; Ham and DeSutter, 1999, 2000). Movement of NH4+ into the soil is dependent on the whole-lagoon seepage rate and the concentration of NH4+ in the liquor (Ham and DeSutter, 1999, 2000). Further movement of NH4+ into the underlying soil also would be a function of the mineralogy and subsequent CEC of the soil (Barrington and Madramootoo, 1989; Ham and DeSutter, 2000) and the effect of competing cations for exchange sites (Lance, 1972).
Movement of NH4+ into the soils underlying lagoons has been investigated by DeSutter et al. (unpublished data, 2005), Ham (2002), Ham and DeSutter (2000), Maule and Fonstad (1996), Miller et al. (1976), Parker et al. (1995), and Perschke and Wright (1998). In most instances, NH4+ concentrations were greatest near the lagoon floor and decreased with depth beneath the lagoon. Although most NH4+ concentrations decreased with depth, NH4+ concentrations were still above expected background levels at the lowermost sampling interval indicating that the true depth of NH4+ movement was not reached by sampling (DeSutter et al., unpublished data, 2005; Ham, 2002; Ham and DeSutter, 2000; Miller et al., 1976; Parker et al., 1995; Perschke and Wright, 1998).
Ham and DeSutter (1999) developed an export model to predict the depth of soil (Zeq), in m, needed to adsorb 100% of the downward moving NH4+ from animal waste lagoons:
 | [1] |
where [NH4+] is the ammonium charge concentration (cmolc m3) in the lagoon liquor; St is the total depth of seepage over the time period of interest (m);
b is the bulk density of the soil (kg m3); CEC is the cation exchange capacity of the soil (cmolc kg1); and EFA is the equivalent fraction of NH4+ on the soil's exchange sites. Although Ham and DeSutter (2000) used an EFA of 0.87 based on the work of Lance (1972), research suggests that the actual EFA in the upper 0.5 m of cattle and swine lagoons may actually be only as much as 0.22 and 0.44, respectively (DeSutter et al., unpublished data, 2005). The EFA of a soil is a function of the selectivity of the soil for NH4+ relative to that of other cations (i.e., Ca2+, Mg2+, K+, and Na+) in the waste for exchange sites (Ham, 2002).
The model was used to predict Zeq, in cm, and compare differing scenarios of CEC and EFA when using a constant concentration of NH4+ in the waste of 910 mg NH4N L1 (6500 cmolc m3), a bulk density of 1400 kg m3, a seepage rate of 1.1 mm d1, and a time of one year (Fig. 1)
. If a soil has an EFA of 0.1 and CEC of 10 cmolc kg1, about 180 cm of soil is needed to adsorb 100% of the downward-moving NH4+ per year, compared with 40 cm of soil if the EFA is 0.5. Also, a 75-cm reduction in Zeq is achieved when the EFA of a 20 cmolc kg3 soil is increased from 0.1 to 0.5. The model and Fig. 1 clearly show that increasing the EFA, CEC, or both can greatly reduce the amount of soil needed to adsorb downward-moving NH4+. Therefore, when evaluating soils for use as liners, choosing soils or soil amendments that improve protection of ground water quality by way of soil chemical properties should be considered. Thus, the objectives of this research were to (i) develop the methodology to evaluate two soils on their ability to selectively adsorb NH4+ by using binary (Ca2+K+, Ca2+NH4+, and K+NH4+) and ternary exchange (Ca2+K+NH4+) reactions; (ii) evaluate these soils with amendments of bentonite or zeolite; and (iii) use the ternary model and information from actual swine- and cattle-waste lagoons to predict the depth of soil needed to adsorb 100% of the downward-moving NH4+.

View larger version (22K):
[in this window]
[in a new window]
|
Fig. 1. Simulation of the depth of soil (Zeq; Eq. [1]) needed to adsorb 100% of the NH4+ moving into the underlying soil of a swine-waste lagoon under various soil conditions after 1 yr. Bd, bulk density of the soil; St, whole-lagoon seepage rate; CEC, cation exchange capacity.
|
|
 |
MATERIALS AND METHODS
|
|---|
Laboratory
Binary exchange relationships (Ca2+K+, Ca2+NH4+, and K+NH4+) were determined on two noncalcareous soils, hereafter referred to as exchangers, from the Great Plains region. Both soils were previously approved for use as compacted soil liners for animal-waste lagoon construction. Soil 1 (S1) was supplied by Kansas Livestock Association Environmental Services (Salina, KS), and was collected in Reno County, Kansas. Soil 2 (S2) was supplied by Agricultural Engineering Associates (Uniontown, KS), and was collected in Brown County, Nebraska. Neither soil had been exposed to waste materials before laboratory experiments. The same binary exchange studies were also performed on soilamendment mixtures, also referred to as exchangers. The two amendments were a Na-bentonite (Envirogel 200; Wyo-Ben, Billings, MT) and clinoptilolite zeolite (BR-40; ZEO, McKinney, TX). The bentonite used for this study was in powder form where 98% of the product was less than 0.18 mm and the zeolite used for this study had a diameter of less than 0.425 mm. The CEC values of the bentonite and zeolite materials were determined from Ca2+K+ exchange following the methods used to determine binary exchange reactions below and were 40 and 67 cmolc kg1, respectively. The determined CEC values for these bentonite and zeolite materials may differ from reported CEC values of other Na-bentonites and clinoptilolite zeolites based on the choice of indicator and replacement cations and also the length of time used for cation extraction (Kithome et al., 1998; Sposito et al., 1983). Individual amendments were mixed with each soil at a concentration of 60 g kg1, in accordance with bentonite mixture ratios used for sealing pond liners (USDA, 2002). Bentonite can be added during liner construction to decrease hydraulic conductivity and zeolite was selected because of its high values of CEC and its preference for NH4+ over Ca2+ and Na+ (Howery and Thomas, 1965; Ming and Mumpton, 1989; Vucinic and Tomanec, 1997).
The <2.0-mm fractions of S1 and S2 were analyzed for sand, silt, and clay following the pipet method outlined in Gee and Bauder (1986). The soils were pretreated in an 80°C water bath with 30% H2O2 and 1 M NaOAc to remove organic matter and carbonates before the determination of the clay concentration. The difference between 1000 g kg1 and the sum of the sand (>53 µm) and silt (253 µm) concentrations was determined to be the concentration of clay (<2 µm). Concentrations of sand, silt, and clay in the Soil 1 + bentonite (S1B), Soil 1 + zeolite (S1Z), Soil 2 + bentonite (S2B), and Soil 2 + zeolite (S2Z) soils were determined by using a modified hydrometer method of Gee and Bauder (1986). The clay fraction was determined by taking a hydrometer reading at a time determined by Stokes' law (about 8 h) after mixing. After the clay fraction was determined, the sand fraction was separated with a 53-µm sieve, and the mass of sand was determined by gravimetric methods after drying at 105°C for 24 h. Silt was determined to be the difference between 1000 g kg1 and the sum of the clay and sand fractions.
Mineralogies of the clay fraction of S1 and S2 were determined by using X-ray diffraction (XRD) techniques outlined in Ransom et al. (1988). The clay fractions from both S1 and S2 were composed of micasmectitekaolinite assemblages. Specific surface areas of the Mg-saturated clay fractions of S1 and S2 were determined according to a modified ethylene glycol monoethyl ether (EGME) method outlined in Carter et al. (1986) in which clays remained Mg-saturated instead of being Ca-saturated. Total C was determined by combustion (CN2000; LECO, St. Joseph, MI). Saturated hydraulic conductivities were determined once for each soil and soil plus amendment following ASTM Designation D 5084-90 (American Society for Testing and Materials, 1991).
Binary-exchange reactions were conducted on Ca- or K-saturated exchangers. One gram of the <1.0-mm size fraction of each exchanger was equilibrated with 4 L of a solution containing varying equivalences of Ca2+ and K+, Ca2+ and NH4+, and K+ and NH4+ (Tables 1 and 2) using a flow-through approach. Solutions were prepared using gravimetric methods from 2 M stock solutions of chloride salts and deionized water while maintaining a constant ionic strength of 0.05 mol L1. The flow-through system consisted of a 47-mm Teflon filter apparatus equipped with a 45-µm filter, peristaltic pump, and Mariotte bottle. For S1 and S2, eleven solutions were used for each binary-exchange reaction. Based on the information provided from the S1 and S2 reactions, S1B, S1Z, S2B, and S2Z were only equilibrated with five solutions. All equilibrations with exchangers were repeated twice and were conducted at 25°C. A 50-mL aliquot of each solution was removed before exposure with the exchanger and was saved for analysis. After equilibrations were achieved, each sample was washed with 100 mL of 95% ethanol. Cations were extracted from a 0.5-g subsample of each exchanger for 30 min using 10 mL of 0.5 M Na-acetate. Calcium and K+ concentrations in the solutions and the 0.5 M Na-acetate extract were determined with inductively coupled plasmaatomic emission spectroscopy (Accuris 141; Fisons Instruments, Beverly, MA). Ammonium N concentrations from solutions and exchanger extracts were determined by the salicylatehypochlorite method (Crooke and Simpson, 1971) and were analyzed with an autoanalyzer (Alpkem, Clackamas, OR). The CEC value of each exchanger (cmolc kg1) after equilibration with a solution was determined to be the sum of the charges of the two cations from the binary reactions. Individual CEC values of all the binary reactions from respective exchangers were summed and divided by the total number of binary reactions to derive the average CEC for each exchanger.
View this table:
[in this window]
[in a new window]
|
Table 1. Binary-exchange (CaK, CaNH4, and KNH4) solution and exchanger compositions for Soil 1 (S1 ) and Soil 2 (S2 ).
|
|
 |
THEORETICAL
|
|---|
Using the general binary exchange reaction, as outlined by Bond (1995) and Bond and Verburg (1997), of:
 | [2] |
where zi is the valence of cation Ii and (exch) and (aq) are the exchanger and solution phases, respectively, the Vanselow selectivity coefficients (KV) and equivalent fractions of cations on the exchanger (Ni) were computed for each binary exchange reaction as defined by Sparks (2003). Solution concentrations were used to calculate solution-phase ionic strengths and the Davies equation was used to calculate the single-ion activity coefficients used to estimate activities from concentrations (Sparks, 2003). Binary-exchange isotherms were constructed according to the approach outlined in Sparks (2003), in which the equivalent fraction on the exchanger phase was plotted as a function of the equivalent fraction in the solution phase. Each exchange isotherm is the average of the two replications. Nonpreference exchange isotherms for the divalentunivalent mixtures were calculated according to the approach outlined in Sposito (2000). Nonpreference isotherms for the univalentunivalent exchanges were 1:1 lines. Gibb's free energies (
Gex0) were calculated for each binary-exchange reaction from:
 | [3] |
and:
 | [4] |
where Kex is the equilibrium exchange constant; R is equal to 8.314 J mol1 K1; and T (K) is the reaction temperature (Sparks, 2003). Each reported Gibb's free-energy value is the average of two replications.
The RothmundKornfeld approach outlined in Bond (1995) and Bond and Verburg (1997) was followed, and ternary exchange between Ca2+K+NH4+ was predicted from GainesThomas binary-exchange reactions (Ca2+K+, Ca2+NH4+, and K+NH4+). The GainesThomas working form of the RothmundKornfeld equation is:
 | [5] |
where
i is the activity coefficient of cation Ii in the solution phase; C is the concentration of cation Ii in the solution phase in mol L1; and kij and nij are the RothmundKornfeld empirical exchange constants. Using the linear model:
 | [6] |
the RothmundKornfeld empirical exchange constants were determined by fitting the binary data to a plot of:
 | [7] |
where log kij and nij are the intercept and slope, respectively. The empirical exchange constants are then used to solve for the predicted ratios of cations on the exchanger phase NK/NNH4, NCa/N2K, and NCa/N2NH4 by:
 | [8a] |
 | [8b] |
 | [8c] |
where A is the solution phase activity (
iCi) in a ternary solution. Predicted equivalent fractions of each cation on the exchanger phase are then determined by describing the mass balance of the exchanger phase as:
 | [9] |
and solving the following quadratic equations:
 | [10a] |
 | [10b] |
 | [10c] |
where
= NK/NNH4, ß = NCa/N2K, and
= NCa/N2NH4. Through the use of Eq. [9] and [10], there are n + 1 equations for n unknowns, which makes it possible to calculate the three equivalent fractions of cations on the exchanger using only Eq. [9] and two of Eq. [10] (Bond and Verburg, 1997). Validity of the model was evaluated by equilibrating the exchanger phases with a ternary solution containing equal equivalences of Ca2+, K+, and NH4+. Also, each exchanger was equilibrated with solutions containing high equivalences (0.69) of each cation, Ca2+, K+, or NH4+, while keeping equivalent fractions of the remaining cations constant. All ternary solutions were prepared with a target ionic strength of 0.05 mol L1 and repeated twice. In addition, each exchanger was equilibrated with the previously described four solutions containing a background Na+ concentration of 0.013 mol L1. Additions of Na+ to binary-exchange reactions have been shown to not influence the exchangers' cationic preference (Sposito et al., 1983). This concentration of Na+ is similar to what one may expect in swine and cattle waste lagoons in Kansas (Ham, 2002). All equilibrium reactions were conducted twice, and the reaction temperature was held constant at 25°C. Finally, lagoon-liquor data from existing swine and cattle waste lagoons in Kansas were used to predict the equivalent fractions of Ca2+, K+, and NH4+ on the exchanger phase. These results were used in conjunction with Eq. [1] to predict the depth of soil needed (cm) to adsorb 100% of the downward-moving NH4+ after 20 yr of operation.
 |
RESULTS AND DISCUSSION
|
|---|
Characterization of Exchangers
The clay contents of S1 and S2 were 354 and 239 g kg1, respectively (Table 3). The addition of bentonite increased the concentration of clay in the S1 and S2 soils whereas the addition of zeolite decreased the concentration of clay. Addition of the amendments slightly decreased the total C in the S1, S1B, and S1Z exchangers, from 7.8 g kg1 to 7.2 and 7.0 g kg1, respectively (Table 3). The total C remained fairly constant at about 4.0 g kg1 in the S2, S2B, and S2Z exchangers.
View this table:
[in this window]
[in a new window]
|
Table 3. Select chemical and physical parameters of two Great Plains soils (S1 and S2) alone, and with additions of bentonite (B; addition of 60 g bentonite kg1) or zeolite (Z; addition of 60 g zeolite kg1).
|
|
The addition of the bentonite and zeolite amendments increased the CEC value of S1 from 17.9 cmolc kg1 to 19.4 and 21.9 cmolc kg1, respectively (Table 3). Although the addition of bentonite only increased the average CEC of S2 by 0.5 cmolc kg1, the addition of zeolite increased the CEC by 3.2 cmolc kg1 (Table 3). The saturated hydraulic conductivities of S1 and S2 were influenced by the addition of the amendments. The saturated hydraulic conductivity of S1 was 2.93 x 108 cm s1 (Table 3). The addition of bentonite and zeolite to S1 slightly increased and decreased the saturated hydraulic conductivities, respectively. As would be expected with the addition of bentonite to soil, the saturated hydraulic conductivity of S2B was 0.29 x 108 cm s1 lower than that of S2 (Table 3). Addition of the more granular zeolite material to S2 increased the saturated hydraulic conductivity by about 12 times, initially making S2Z undesirable for liner construction because it exceeded the critical value of 1 x 107 cm s1 (Table 3) (USDA, 1992). Although, Ham and DeSutter (1999) determined that laboratory-derived saturated hydraulic conductivities might be as much as five times greater than those calculated from whole-lagoon seepage test results. The specific surface areas of the clay fractions of S1 and S2 were 334 and 361 m2 g1, respectively (Table 3).
Binary Exchange
The repeatability of the cation concentrations on the exchangers was generally less than a 10% relative percent difference (RPD) from the means indicating that the flow-through method worked very well for equilibrating an exchanger phase with a solution phase. Thus, the average of the two replications was reported (Tables 1 and 2). For S1 and S2, complete exchange of the initial saturating cation was observed for all reactions when the saturating cation was not present in the equilibrating solutions (Table 1). Because of this behavior, equivalent fractions of 1.0 and 0.0 were added to the S1B, S1Z, S2B, and S2Z datasets to allow the
Gex0 to be calculated for all six exchangers to be comparable. The RPD from the means of the
G values, reported in Fig. 2 through 4
, indicate that the smallest RPD was from Ca2+K+ exchange reaction on the S2B exchanger (0.3%) and the largest RPD was from the K+NH4+ exchange reaction on the S1B exchanger (44.7%). In general, most RPDs for the
Gex0 values were less than 10%.
Binary-exchange isotherms for S1, S1B, and S1Z were all below the nonpreference isotherms, indicating exchanger preference for K+ over Ca2+ for the Ca2+K+ exchange reaction (Fig. 2a). Gibb's free-energy calculations supported this conclusion where
Gex0 values for S1, S1B, and S1Z exchangers were 4.66, 5.36, and 8.17 kJ mol1, respectively. The exchange isotherm and the
Gex0 value for the S1Z exchanger were lower and about 1.5 times more negative, respectively, than the exchange isotherms and
Gex0 values for the S1 and S1B exchangers, indicating a stronger preference for K+ over Ca2+ when zeolite was present. Exchange isotherms and
Gex0 values also indicate that K+ is preferred over Ca2+ on the S2, S2B, and S2Z exchangers (Fig. 2b). Again, S2Z had the greatest preference for K+ over Ca2+ compared with the S2 and S2B exchangers as observed by the exchange isotherm and a
Gex0 value of 8.99 kJ mol1.
Isotherms and
Gex0 values both indicate a preference for NH4+ over Ca2+ for the S1, S1B, and S1Z exchangers (Fig. 3a). The addition of zeolite to the S1 exchanger increased NH4+ preference over both the S1 and S1B exchangers (Fig. 3a). Exchanger S2Z had the greatest preference for NH4+ over Ca2+, as indicated by the exchange isotherm and a
Gex0 value of 8.50 kJ mol1 (Fig. 3b). The addition of bentonite to both the S1 and S2 exchangers had little affect on the preference of S1B and S2B exchangers for NH4+ over Ca2+ (Fig. 3a and 3b).
Exchange isotherms and
Gex0 values indicated a slightly greater preference for K+ over NH4+ for all six exchangers (Fig. 4). Addition of bentonite and zeolite to S1 numerically increased the
Gex0 values for K+ over NH4+, but statistical differences were not determined (Fig. 4a). Gibb's free energy values for the S2, S2B, and S2Z exchangers were very similar (Fig. 4b).
Ternary Exchange
Binary-exchange data were used to calculate the intercept and slope, log kij and nij, respectively, of the plot of Eq. [7] (Table 4). The r2 values of these plots indicate a good fit of the log-transformed data. Values of r2 ranged from 0.985 to 0.999 from all the plots. Thus, the determined values of kij and nij can be considered good estimates of the empirical exchange constants. Through the use of Eq. [8], [9], and [10], kij and nij were used to estimate Ni for each cation in a ternary system. These predicted Ni values were compared with the measured Ni for each cation in the ternary system. The relationship between the predicted and measured Ni on the S1 exchanger was very good, with an r2 of 0.966 and a root-mean square (RMS) of 0.038 (Fig. 5a) . The addition of bentonite did not hinder the predictability of Ni, but addition of zeolite decreased the r2 value to 0.881 and increased the RMS to 0.055. Adding a Na+ background to the equilibrium solutions did not decrease the model's ability to predict Ni in a ternary system (Fig. 5b). Both the r2 and RMS values of S1 and S1B slightly improved when Na was present in solution, although the r2 and RMS values from the S1Z exchanger slightly deteriorated with the presence of Na+ in the solution. No significant differences were observed between the exchangers with, and without, Na+ present in the solution phase (data not shown).
View this table:
[in this window]
[in a new window]
|
Table 4. RothmundKornfeld fitted parameters for the GainesThomas convention for two Great Plains soils (S1 and S2) alone, and with additions of bentonite (B; addition of 60 g bentonite kg1) or zeolite (Z; addition of 60 g zeolite kg1).
|
|

View larger version (30K):
[in this window]
[in a new window]
|
Fig. 5. (a) Equivalent fraction of cations on the Soil 1 (S1), Soil 1 + bentonite (S1B), and Soil 1 + zeolite (S1Z) exchanger phases predicted from the RothmundKornfeld equation and GainesThomas convention vs. laboratory-measured equivalent fractions on the S1, S1B, and S1Z exchangers. (b) Also with Na added to the ternary solutions at 0.013 mol L1. Root-mean square (RMS) = [( predicted N measured N)2/n]0.5.
|
|
The RothmundKornfeld approach also worked very well for predicting Ni on the S2, S2B, and S2Z exchangers, compared with the measured results. Values of r2 were 0.975 and 0.972 for the S2 and S2B exchangers, respectively, with RMS values less than 0.04 (Fig. 6a)
. Although the addition of zeolite decreased the r2 value to 0.859, the RothmundKornfeld method still provided a good estimate of Ni on the S2Z exchanger. The addition of 0.013 mol Na L1 to the solution phase slightly decreased the r2 value and increased the RMS value for the S2B exchanger (Fig. 6b). Although the S2 exchanger was relatively unchanged by addition of Na+, the r2 value slightly increased with the S2Z exchanger, to 0.918. Again, no significant differences were observed between the exchangers with, and without, Na+ present in the solution phase (data not shown).

View larger version (30K):
[in this window]
[in a new window]
|
Fig. 6. (a) Equivalent fraction of cations on the Soil 2 (S2), Soil 2 + bentonite (S2B), and Soil 2 + zeolite (S2Z) exchanger phases predicted from the RothmundKornfeld equation and GainesThomas convention vs. laboratory-measured equivalent fractions on the S2, S2B, and S2Z exchangers. (b) Also with Na added to the ternary solutions at 0.013 mol L1. Root-mean square (RMS) = [( predicted N measured N)2/n]0.5.
|
|
Application of the RothmundKornfeld Approach
The binary and ternary assessments of the six exchangers, in conjunction with Eq. [1], were applied to make predictions of the equivalent depth of soil needed to retain 100% of the NH4+ seeping from four animal waste lagoons in Kansas. This approach assumes piston displacement of the existing soil water with the incoming fluid where molecular diffusion and hydrodynamic dispersion are not considered. Swine 1 and Swine 2 were finishing and nursery operations, respectively, located in southwestern Kansas. Cattle 1 and Cattle 2 were runoff lagoons located in central and southwestern Kansas, respectively. The concentration of NH4+ in the swine lagoons was as much as six times greater than in the cattle lagoons, whereas Ca2+ concentrations were greater in the cattle lagoons (Table 5). Potassium concentrations were more consistent in the cattle lagoons than in the swine lagoons.
The information provided in Table 5 was inserted into the RothmundKornfeld model for each soil and soil-plus-amendment combination to predict the equivalent fraction of NH4+ on the exchangers (NNH4). The predicted NNH4 value for each exchanger was inserted into Eq. [1], along with the concentration of NH4+ in the waste, to estimate the equivalent depth of soil needed to retain 100% of the NH4+ that would be predicted to seep from each lagoon over 20 yr. The bulk density of the soil was assumed to be 1300 kg m3. Estimated exchanger NNH4 values from the Swine 1 lagoon were very similar, with values ranging from 0.51 to 0.60 (Table 6). From these exchangers, S1Z and S2Z had the highest NNH4 values, with 0.55 and 0.60, respectively. The range of exchanger NNH4 values from Swine 2, which had 2.57 and 2.86 cmolc L1 more NH4+ and K+, respectively, in the liquor than did Swine 1, were from 0.46 to 0.52. The addition of the amendments did not increase the NNH4 values above those of the native soils with Swine 2.
View this table:
[in this window]
[in a new window]
|
Table 6. Predicted equivalent fraction of NH4 (NNH4) on an exchanger in a ternary system according to the RothmundKornfeld approach, and the estimated depth of soil (Zeq) needed after 20 yr to retain 100% of the NH4N seeping from four animal waste lagoons.
|
|
The predicted NNH4 values for the cattle lagoons were about half of the swine-lagoon values (Table 6) due to the much lower NH4+ concentrations in the cattle-lagoon liquor than in the swine lagoon liquor (Table 5). Also, in conjunction with the NH4+ concentrations, overall NNH4 values from the Cattle 1 were less than values from Cattle 2. In Cattle 1 and Cattle 2, the addition of zeolite increased the NNH4 values for S1 and S2 by up to 0.05 and 0.1, respectively. But the addition of bentonite to S1 and S2 decreased or did not change the NNH4 values above values of the native soils. Overall, the addition of bentonite only increased NNH4 values 50% of the time whereas addition of zeolite increased the NNH4 values 75% of the time. However, the range of NNH4 values for S1 and S2 with and without amendments is relatively small, and indicates that the amendments used here did not greatly increase the exchangers' selectivity for NH4+.
Although the benefit of adding bentonite and zeolite to soils with respect to the predicted NNH4 values may be small, these slight increases (Table 6), coupled with increases in CEC (Table 3), will ultimately decrease the depth of soil needed to adsorb downward-moving NH4+. The value of Zeq required after 20 yr for the six exchangers ranged between 223 and 318 cm for Swine 1 and between 425 and 545 cm for Swine 2 (Table 6). For both Swine 1 and Swine 2, the benefit of adding bentonite to the native soils was very subtle, but the addition of zeolite to S1 and S2 decreased Z by 55 and 73 cm, respectively, in the Swine 1 lagoon. Although Zeq was as much as two times greater in Swine 2 than in Swine 1, the trend of results was very similar. Bentonite additions to S1 and S2 did very little in decreasing Zeq for Swine 2, but additions of zeolite reduced the required depth by up to 86 cm. Bentonite additions only decreased Zeq in the S2B exchangers over the S2 in both cattle lagoons by 10 cm or less (Table 6). Conversely, the S1B exchangers for both cattle lagoons slightly increased the predicted amount of soil, compared with S1. Again, zeolite additions to S1 and S2 greatly decreased Zeq for both lagoons by as much as 133 cm, compared with the Zeq of the native soils.
Overall, additions of bentonite to the native soils were not very effective at reducing Zeq in any of the lagoons. Considering all S1B and S2B combinations with lagoons, the average reduction in Zeq was only about 6 cm, although additions of zeolite to both soils decreased the Zeq from all lagoons by an average of about 78 cm when compared with that of the native soils. Thus, from a soil chemistry perspective, the benefits of adding bentonite to the soil to reduce Zeq may not be as important as the beneficial decrease of a soil's saturated hydraulic conductivity (Table 3). Conversely, zeolite additions to native soil greatly decreased the required Zeq under lagoons, which can mostly be attributed to increases in the CEC of the mixtures. Additions of zeolite, however, did tend to increase the saturated hydraulic conductivities of the soil and may not be appropriate for use with some soils.
 |
CONCLUSIONS
|
|---|
The CEC of native soils may be increased through the use of amendments such as bentonite and zeolite. Additions of zeolite at a rate of 60 g kg1 to two soils that were previously approved for use as compacted soil liners increased the CEC of the mixtures by as much as 4.0 cmolc kg1. Binary-exchange reactions between Ca2+K+, Ca2+NH4+, and K+NH4+ that were conducted on these soils, and soils plus amendments, indicated that the univalent cations, K+ and NH4+, were preferred over Ca2+ on the exchangers. Also, K+ was preferred over NH4+ in all instances. Mixtures of soil plus zeolite consistently had the greatest preference for K+ or NH4+ over Ca2+ and for K+ over NH4+, compared with preferences of the native soils and soils with bentonite added. But, preferences for K+ and NH4+ over Ca2+ and K+ over NH4+ in the native soils and soilbentonite mixtures were nearly identical, which suggests the presence of similar smectite clays in the soils and bentonite amendments.
The RothmundKornfeld approach, according to GainesThomas convention, yielded good predictions of the exchanger phase composition in a ternary system by use of simple binary-exchange reactions. Although these experiments only addressed Ca2+, K+, and NH4+, the addition of Na+ to the solution phase, essentially producing a quaternary system, did not hinder the ability of the model to effectively predict equivalent fractions of Ca2+, K+, or NH4+ on the exchanger phase. Research suggests that combining Mg2+ with Ca2+ (Feigenbaum et al., 1991) may also be a viable option for increasing the practicality of the model when investigating environmental issues.
The RothmundKornfeld approach, in conjunction with the knowledge of lagoon-liquor chemistries, may provide a useful tool to aid in the construction of soil liners for animal and municipal waste lagoons (Ham and DeSutter, 2000). This approach allows one to use soil from a proposed construction site and evaluate a soil's ability to selectively adsorb NH4+ on its exchange sites. Results from the ternary-prediction experiments indicate that soils from the Great Plains that are used for liner materials could retain about 53 and 23% of the NH4+ from swine and cattle lagoons on their exchange sites. Additions of bentonite or zeolite to these soils may not improve the equivalent fraction of NH4+ on the exchange sites, when used in conjunction with a swine-waste lagoon. When NH4+ concentrations are similar to those found in cattle lagoons, however, the addition of zeolite to soils could help increase the equivalent fraction of NH4+ on the exchange sites.
Another useful application of the RothmundKornfeld approach, when used in conjunction with Eq. [1], allows for the prediction of the equivalent depth of soil needed to adsorb 100% of the downward-moving NH4+ from waste lagoons. In general, swine lagoons needed a greater equivalent depth of soil, compared with cattle waste lagoons, to adsorb all the NH4+. The equivalent depth of soil needed to adsorb 100% of the NH4+ from a swine lagoon after 20 yr of operation was as much as 5 m, whereas the soil needed for a cattle lagoon was only as much as about 3 m.
Additional research that should be explored is how effectively the predicted Zeq values compare with actual depths of NH4+ movement below swine and cattle lagoons and how downward-moving organic C from animal waste lagoons affects an exchanger's cation selectivity preference and cation exchange capacity. Model validation through the investigation of decommissioned swine and cattle lagoons would allow researchers and regulators to better assess the potential costs for N cleanup and also better identify lagoon construction sites. Increases in CEC are expected with the addition of organic carbon (Parfitt et al., 1995), but some uncertainty in NH4+ selectivity still exists when organic C concentrations are increased (Chung and Zasoski, 1994). Another area of exploration is the evaluation of calcareous soils for use as soil liners. Dissolution of carbonates may increase the Ca2+ and/or Mg2+ content in calcareous soil and, thus, compete with the cations in the lagoon liquor for soil exchange sites and subsequently increase the amount of soil needed to adsorb downward-moving NH4+. Prediction of actual concentration profiles of NH4+ beneath the lagoon may be possible by incorporating exchange theory with transport theory (Mansell et al., 1988; Mansell et al., 1993; Bond, 1997). Investigating other amendments, such as vermiculite, that may have a greater ability to selectively adsorb NH4+ (Evangelou and Lumbanraja, 2002) may help reduce the downward movement of NH4+.
 |
NOTES
|
|---|
Contribution no. 04-385-J from the Kansas Agricultural Experiment Station, Manhattan, KS.
 |
REFERENCES
|
|---|
- American Society for Testing and Materials. 1991. ASTM Designation D 5084-90. Standard test method for measurement of hydraulic conductivity of saturated porous materials using a flexible wall permeameter. Annu. Book ASTM Stand. 4.08:10701077.
- Barrington, S.F., and C.A. Madramootoo. 1989. Investigating seal formation from manure infiltration into soils. Trans. ASAE 32:851856.
- Bond, W.J. 1995. On the RothmundKornfeld description of cation exchange. Soil Sci. Soc. Am. J. 59:436443.[Abstract/Free Full Text]
- Bond, W.J. 1997. Competitive exchange of K+, Na+, and Ca2+ during transport through soil. Aust. J. Soil Res. 35:739752.
- Bond, W.J., and K. Verburg. 1997. Comparison of methods for predicting ternary exchange from binary isotherms. Soil Sci. Soc. Am. J. 61:444454.[Abstract/Free Full Text]
- Carter, D.L., M.M. Mortland, and W.D. Kemper. 1986. Specific surface. p. 413423. In A. Klute (ed.) Methods of soil analysis. Part 1. 2nd ed. Agron. Monogr. 9. ASA and SSSA, Madison, WI.
- Chung, J.B., and R.J. Zasoski. 1994. Ammonium-potassium and ammonium-calcium exchange equilibria in bulk and rhizosphere soil. Soil Sci. Soc. Am. J. 58:13681375.[Abstract/Free Full Text]
- Crooke, W.M., and W.E. Simpson. 1971. Determination of ammonium in Kjeldahl digests of crops by an automated procedure. J. Sci. Food Agric. 22:910.
- Evangelou, V.P., and J. Lumbanraja. 2002. Ammonium-potassium-calcium exchange on vermiculite and hydroxyl-aluminum vermiculite. Soil Sci. Soc. Am. J. 66:445455.[Abstract/Free Full Text]
- Feigenbaum, S., A. Bar-Tel, R. Portnoy, and D.L. Sparks. 1991. Binary and ternary exchange of potassium on calcareous montmorillonitic soils. Soil Sci. Soc. Am. J. 55:4956.[Abstract/Free Full Text]
- Gee, G.W., and J.W. Bauder. 1986. Particle-size analysis. p. 383411. In A. Klute (ed.) Methods of soil analysis. Part 1. 2nd ed. Agron. Monogr. 9. ASA and SSSA, Madison, WI.
- Ham, J.M. 2002. Seepage losses from animal waste lagoons: A summary of a four-year investigation in Kansas. Trans. ASAE 45:983992.
- Ham, J.M., and T.M. DeSutter. 1999. Seepage losses and nitrogen export from swine-waste lagoons: A water balance study. J. Environ. Qual. 28:10901099.[Abstract/Free Full Text]
- Ham, J.M., and T.M. DeSutter. 2000. Toward site-specific design standards for animal-waste lagoons: Protecting ground water quality. J. Environ. Qual. 29:17211732.[Abstract/Free Full Text]
- Howery, D.G., and H.C. Thomas. 1965. Ion exchange on the mineral clinoptilolite. J. Phys. Chem. 69:531537.
- Kithome, M., J.W. Paul, L.M. Lavkulich, and A.A. Bomke. 1998. Kinetics of ammonium adsorption and desorption by the natural zeolite clinoptilolite. Soil Sci. Soc. Am. J. 62:622629.[Abstract/Free Full Text]
- Lance, J.C. 1972. Nitrogen removal by soil mechanisms. J. Water Pollut. Control Fed. 44:13521361.
- Mansell, R.S., S.A. Bloom, H.M. Selim, and R.D. Rhue. 1988. Simulated transport of multiple cations in soil using variable selectivity coefficients. Soil Sci. Soc. Am. J. 52:15331540.[Abstract/Free Full Text]
- Mansell, R.S., W.J. Bond, and S.A. Bloom. 1993. Simulating cation transport during water flow in soil: Two approaches. Soil Sci. Soc. Am. J. 57:39.
- Maule, C.P., and T.A. Fonstad. 1996. Ion movement immediately beneath earthen hog manure storages: How much, how deep, and how fast. ASAE Paper no. 962049. ASAE, St. Joseph, MI.
- Miller, M.H., J.B. Robinson, and D.W. Gallagher. 1976. Accumulation of nutrients in soil beneath hog manure lagoons. J. Environ. Qual. 5:279282.
- Ming, D.W., and F.A. Mumpton. 1989. Zeolites in soil. p. 873911. In J.B. Dixon and S.B. Weed (ed.) Minerals in soil environments. 2nd ed. SSSA Book Ser. 1. ASA and SSSA, Madison, WI.
- Parfitt, R.L., D.J. Giltrap, and J.S. Whitton. 1995. Contribution of organic matter and clay minerals to the cation exchange capacity of soils. Commun. Soil Sci. Plant Anal. 26:13431355.
- Parker, D.B., J.A. Nienaber, D.E. Eisenhauer, and D.D. Schulte. 1995. Unsaturated seepage from a feedlot runoff storage pond. ASAE Paper no. 951765. ASAE, St. Joseph, MI.
- Perschke, S.P., and P.E. Wright. 1998. Evaluation of soil under an earthen manure storage pond. ASAE Paper no. 982145. ASAE, St. Joseph, MI.
- Ransom, M.D., J.M. Bigham, N.E. Smeck, and W.F. Jaynes. 1988. Transitional vermiculite-smectite phases in aqualfs of southwestern Ohio. Soil Sci. Soc. Am. J. 52:873880.[Abstract/Free Full Text]
- Sparks, D.L. 2003. Environmental soil chemistry. 2nd ed. Academic Press, Boston.
- Sposito, G. 2000. Ion exchange phenomena. p. B-241 to B-264. In M.E. Sumner (ed.) Handbook of soil science. CRC Press, Boca Raton, FL.
- Sposito, G., C. Jouany, K.M. Holtzclaw, and C.S. LeVesque. 1983. Calcium-magnesium exchange on Wyoming bentonite in the presence of adsorbed sodium. Soil Sci. Soc. Am. J. 47:10811085.[Abstract/Free Full Text]
- USDA. 1992. Agricultural waste management field handbook. Natural Resources Conservation Service.
- USDA. 2002. Pond sealing or lining bentonite treatment. NRCS Conservation Practice Standard 521C-1. USDA, Washington, DC.
- Vucinic, D., and R. Tomanec. 1997. Ca2+
2NH4+ ion exchange reaction on clinoptilolite. p. 199206. In Kirov et al (ed.) Proc. Sofia Zeolite Meeting '95, Sofia, Bulgaria. 1825 June 1995. Pensoft, Moscow.
Related articles in JEQ:
- This Issue in Journal of Environmental Quality
JEQ 2005 34: ix.
[Full Text]
This article has been cited by other articles:

|
 |

|
 |
 
T. M. DeSutter, G. M. Pierzynski, and L. R. Baker
Flow-Through and Batch Methods for Determining Calcium-Magnesium and Magnesium-Calcium Selectivity
Soil Sci. Soc. Am. J.,
February 27, 2006;
70(2):
550 - 554.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
T. M. DeSutter, G. M. Pierzynski, and J. M. Ham
Movement of Lagoon-Liquor Constituents below Four Animal-Waste Lagoons
J. Environ. Qual.,
June 7, 2005;
34(4):
1234 - 1242.
[Abstract]
[Full Text]
[PDF]
|
 |
|