Environ. Sci. Technol. 2000, 34, 2215-2223
Modeling Cd and Zn Sorption to Hydrous Metal Oxides P. TRIVEDI AND L. AXE* Department of Civil and Environmental Engineering, New Jersey Institute of Technology, University Heights, Newark, New Jersey 07102
The mobility and bioavailability of Cd and Zn in soils and sediments are affected by contaminant distribution mechanisms. One important process is sorption to hydrous aluminum, iron, and manganese oxides, which are ubiquitous in soils and sediments as both discrete particles and coatings and exhibit a high affinity for these metals. Mechanistic models are required for accurately assessing risks to populations and in the long-term management of contaminated soils and sediments. This research demonstrates intraparticle diffusion is the rate-limiting step in the sorption of Cd and Zn to microporous oxides. Furthermore, as much as 90% of the total sorption sites on the oxides reside on the micropore walls. Because longterm experiments require a lengthy period of time, predictive methods would be useful for determining surface diffusivities. Theoretically, surface diffusivities can be predicted from site activation theory, which is based on the random walk model where atoms or molecules vibrate at localized sites along the surface. Once the vibrating ion has sufficient energy, it will jump to a neighboring site. For a given metal, the associated activation energy was observed to be equivalent for all three oxides; in an effort to predict this energy, a correlation is presented between the adsorption enthalpy and the adsorbate hydrated radius. For each oxide, the Polanyi constant (R) that relates adsorption enthalpy and activation energy was equivalent for the transition metals studied.
Introduction Heavy metals such as cadmium and zinc are released into aquatic and soil environments largely from various anthropogenic activities and pose a serious threat to plants, animals, and humans because of their toxicity and persistence (1-4). Extensive research (5-9) has demonstrated that the mobility and bioavailability of these metals are impacted by sorption to hydrous oxides of aluminum (HAO), iron (HFO), and manganese (HMO), which are ubiquitous in soils and sediments. To assess risks to the surrounding environment and manage remedial activities requires accurate models and well-defined transport parameters. Although limited, studies have shown that sorption of heavy metals to hydrous oxides is a two-step process (10-16): a rapid adsorption of metal ions to the external surface is followed by slow intraparticle diffusion along the oxide micropore walls. This intraparticle diffusion is the rate-limiting mechanism in the sorption process. Studies including the present one demonstrate that * To whom correspondence should be addressed. Phone: (973) 596-2477. Fax: (973) 596-5790. E-mail:
[email protected]. 10.1021/es991110c CCC: $19.00 Published on Web 04/25/2000
2000 American Chemical Society
FIGURE 1. Adsorption edges of Cd and Zn adsorption to hydrous metal oxides at 25 °C and as a function of ionic strength (µ). For all Cd studies, the initial concentration was 5 × 10-5 M [except for Cd adsorption to HFO, 2 × 10-7 M (15)] and for all Zn studies the initial concentration was 5 × 10-5 M. HAO and HFO concentrations were 1 g L-1, and HMO was 0.1 g L-1. the long-term experiments used to observe and model diffusion require lengthy periods of time and therefore methods to predict diffusivities are needed. Axe and coworkers (15, 16) have estimated surface diffusivities using site activation theory and assuming a sinusoidal potential field along the pore surface (17):
Ds ) λ[Ea/2m]1/2 exp[-Ea/(RT)]
(1)
where λ is the mean distance between sites, m is the molecular weight of the diffusing species, Ea is the activation energy required for a sorbed ion to jump to the neighboring site, and exp[-Ea/(RT)] is the Boltzmann factor. The adsorption enthalpy (∆H°) is related to Ea through the Polanyi relation Ea ) R (∆H°), where R is the Polanyi constant (17, 18). In this paper, experimental surface diffusivities are determined, and theoretical ones are estimated. With the objective of developing a predictive method, two hypotheses are tested. (1) For a specific metal ion, EA is equivalent for all hydrous oxides. (2) For a specific hydrous oxide, metals from the one group of the Periodic Table form similar sorption complexes with the hydrous oxide surface and hence have the same R. Several studies have revealed the adsorption affinity for metal ions to oxides follows the order of Pb > Cu > Zn > Cd > Co > Mn > Sr > Ca (5-8, 19-25). Gray (26) compared VOL. 34, NO. 11, 2000 / ENVIRONMENTAL SCIENCE & TECHNOLOGY
9
2215
FIGURE 2. Isotherms of Cd and Zn adsorption to hydrous metal oxides at 25 °C and different pH values. Solid lines represent the linear distribution model; Kd (L g-1) represents the distribution coefficient. HAO and HFO (15) concentrations were 1 g L-1, and HMO was 0.1 g L-1. the trends of cation adsorption to manganese dioxide with Hofmeister series and showed that the adsorption sequence for alkaline earth metals and certain heavy metals follows the reverse order of their hydrated radii. Pan and Liss (27) proposed that Zn, a much harder Lewis acid than Cd2+, binds more strongly to goethite. In efforts to predict adsorption, Morgan and Stumm (5) and others (24, 27-31) have correlated the adsorption equilibrium constant to the first hydrolysis constant of the metal. Wehrli et al. (32) developed a linear relation between the rate constant for water exchange and the intrinsic adsorption rate constants. For the most part, correlations developed in previous studies (33-38) were based on one temperature; hence, the relations cannot be used for other conditions. In a number of studies (5, 27-29), however, adsorption of metal ions to oxides has been modeled 2216
9
ENVIRONMENTAL SCIENCE & TECHNOLOGY / VOL. 34, NO. 11, 2000
as a function of temperature to obtain the adsorption enthalpies. From the van’t Hoff equation, adsorption enthalpy is related to the adsorption affinity or equilibrium constant (39). Adsorption enthalpy represents the bonding energy between the adsorbate ion and the sorbent; it is primarily a function of the intermolecular forces of attraction (40, 41) and when electrostatic ones dominate:
∆H° ) f(Z2/RH)
(2)
where Z is valence charge of a metal ion and RH is its hydrated radius. In this paper, in an effort to predict enthalpies and activation energies, two correlations are presented with one relating adsorption enthalpies and hydrated cationic radii
and the other relating enthalpies and activation energies. These correlations can be used to help predict sorption behavior.
Materials and Methods All sorption studies were conducted with freshly precipitated oxides. HAO and HMO were synthesized according to the methods detailed previously (16). HFO was precipitated as described by Dzombak and Morel (42). These porous amorphous oxides have been characterized in earlier works (11, 16). Solutions were prepared with reagent-grade chemicals and Millipore-Q water following standard methods (43). Stock solutions were tagged with 109Cd and 65Zn; their activities in the samples of suspension and filtrate were measured with a Beckmann LS6000SE liquid scintillation counter (16). Two types of sorption experiments were conducted: (i) short-term (i.e., 4 h contact time) and (ii) long-term ones that were designed to evaluate diffusivities. Turbulent hydraulic conditions (Re g 2.5 × 105 with respect to the reactor diameter) were maintained in all experiments to minimize the external mass-transfer resistance (44, 45). Procedures for short-term studies of adsorption edges and isotherms are conventional and described elsewhere (11, 15, 16). In earlier short-term studies (11, 15, 16), a contact time of 4 h demonstrated equilibrium (or pseudo-equilibrium) was reached between the bulk aqueous phase and that sorbed to the external surface. For each system, adsorption edges were studied at 25 °C and two ionic strengths. Isotherm studies of Cd and of Zn were conducted at 4, 14, and 25 °C and at multiple pHs to evaluate equilibrium constants, site densities, and adsorption enthalpies. In addition, for select systems, a Phillips X’Pert X-ray diffractometer (XRD) with Ni-filtered CuKR radiation was employed to study Zn sorption to HMO and Cd sorption to HAO and to determine whether in the sorption process there were any detectable variations in the XRD patterns. The long-term or constant boundary condition (CBC) experiments were used to study the slow sorption process of intraparticle surface diffusion. In these studies, the metal ion concentration in the bulk aqueous phase was maintained approximately constant by continuously monitoring its concentration and adding adsorbate as needed (11, 15, 16). Therefore, the adsorbate concentration on the external surface of the particle was maintained constant, Cs. Under the Results and Discussion, the analytical solution (46) to the mass balance given the initial and boundary conditions is integrated over the volume of the spherical aggregated particles. For these experimental studies, surface diffusivity is the only fitting parameter in the model. By using the three ubiquitous amorphous oxides and two metals from the same group in the Periodic Table, this research involved testing the hypothesis that, in surface diffusion and based on site activation theory, the activation energy required to jump to the neighboring site is a function of the metal type and that this type behaves similarly with a specific oxide.
Results and Discussion X-ray diffraction profiles of hydrous metal oxides revealed that these oxides remain amorphous when Cd or Zn was sorbed to their surfaces. There was no evidence of a surface precipitate or solid solution formation in metal sorption to the hydrous oxides. The diffraction studies showed that HAO and HMO aged 6 months with and without the adsorbate remained amorphous. Many other studies (11, 47-49) have found that sorption may inhibit oxide crystallization; however, for the duration of our studies, we did not observe crystallization with or without contaminant sorbed. Adsorption edges for Cd and Zn adsorption to hydrous metal oxides (Figure 1) illustrate the sigmoid curve charac-
FIGURE 3. Effect of temperature on Cd and Zn adsorption to hydrous oxides. Solid lines represent Van’t Hoff model, and the data points represent adsorption equilibrium constant (K) evaluated from the isotherm studies. Reported values are adsorption enthalpies (kcal mol-1). teristic of transition metals; these results are consistent with others (7, 13, 15, 19, 50-54). The affinity of Cd and Zn ions for the hydrous oxide surfaces follows the order of HMO > HFO > HAO; this result is in agreement with others (13, 19, 51, 54). Comparing these results with earlier ones of strontium (16, 22), the affinity of cations to hydrous metal oxides follows the trend of Zn > Cd > Sr and is consistent with the cation affinity order discussed earlier. For all systems, the amount of Cd2+ and Zn2+ sorbed decreased with increase in ionic strength, which suggests that metal ions retain their waters of hydration upon sorption to amorphous oxides (55). Many other studies have found similar results (15, 16, 29, 56). Conversely, many others observed no significant effects of ionic strengths on cation adsorption to metal oxides (28, 55, 57-59); however, most of these studies were conducted with crystalline oxides (28, 55, 59). Krapiel et al. (29) attributed the decrease in adsorption of Zn and Ni to montmorillonite with increase in ionic strength to an electrostatic cation exchange process. Misak et al. (13) demonstrated that the decrease in the activity coefficient with increase in ionic strength could explain the adsorption edge shift to the right. Applying activity coefficients to the present studies could not explain the decrease in sorption with increasing ionic strength; one probable reason is different electrolytes. Misak et al. (13) used NaCl when studying Co and Zn sorption to ferric and stannic oxides, where metal-complexation increases with increasing electrolyte concentration. However, NO3- does not form metal complexes, and its contribution is limited to the suppression of the diffuse layer with increase in ionic strength. To study the effect of concentration on adsorption, isotherms for cadmium and zinc adsorption to hydrous oxides (Figure 2) demonstrated a linear relationship (R2 g 0.96) between the sorbed and the bulk aqueous concentrations. VOL. 34, NO. 11, 2000 / ENVIRONMENTAL SCIENCE & TECHNOLOGY
9
2217
FIGURE 4. CBC studies with 1 g L-1 HAO at pH 7 and 25 °C: (a) experimental data of Cd sorption [Cd]bulk ) 1.5 × 10-8 M, (b) internal sorption modeled for Cd, (c) experimental data of Zn sorption [Zn]bulk ) 1.3 × 10-9 M, and (d) internal sorption modeled for Zn. Dashed lines determined from the propagation of errors (POE) method provide the error or standard deviation (SD).
TABLE 1. Cd and Zn Sorption Parameters at pH 7 parametersa
HAO-Cdb
HAO-Znb
HFO-Cdc,d
HFO-Znc
HMO-Cdc
HMO-Znc
Ct (mol of Sr/g) K ∆G° (kcal/mol) ∆S° (kcal/mol K)e ∆H° (kcal/mol) Ea (kcal/mol) R (Polanyi relation constant) theor. Ds (cm2/s) expt. Ds (cm2/s)
0.011 5011 -5.03 0.07 15.98 13.10 0.82 2.51 × 10-13 2.5 × 10-13
0.011 5364 -5.07 0.08 18.80 14.50 0.77 5.41 × 10-14 5.4 × 10-14
0.025 2600 -4.64 0.09 21.50 13.30 0.62 1.2 × 10-14 1.0 × 10-14
0.025 4000 -4.90 0.09 22.70 15.70 0.69 8.53 × 10-16 8.49 × 10-16
0.034 134 653 -6.97 0.10 23.10 13.35 0.58 1.19 × 10-13 1.14 × 10-13
0.034 234 529 -7.30 0.11 24.99 15.20 0.61 9.65 × 10-15 9.5 × 10-15
a All parameters are based on 25 °C except the enthalpy. b (8% error. c (10% error. is calculated from the thermodynamic relation (39), ∆G° ) ∆H° - T ∆S°.
This relationship suggests one average type of adsorption site present on the oxide surface. During these studies, no significant change in pH was observed, however, the oxide 2218
9
ENVIRONMENTAL SCIENCE & TECHNOLOGY / VOL. 34, NO. 11, 2000
d
Data for Cd sorption to HFO obtained from ref 15. e ∆S°
may act as a buffer. In a number of studies, the Langmuir model has been employed to describe metal sorption to oxides: Zn and Pb adsorption to goethite (33, 34); Zn
FIGURE 5. CBC studies and modeling results at pH 7 and 25 °C, where the externally sorbed contribution has been subtracted: (a) Cd internal sorption to 1 g L-1 HFO with [Cd]bulk ) 1.0 × 10-8 M (two studies conducted)(15), (b) Zn internal sorption to 1 g L-1 HFO with [Zn]bulk ) 1.5 × 10-8 M, (c) Cd internal sorption to 0.1 g L-1 HMO with [Cd]bulk ) 1.5 × 10-9 M, and (d) Zn internal sorption to 0.1 g L-1 HMO with [Zn]bulk ) 1.2 × 10-9 M. Dashed lines determined from the propagation of errors (POE) method provide the error or standard deviation (SD). adsorption to Al and Fe hydrous oxides (50); Cd adsorption to hydrous manganese oxides (60); Pb adsorption to hydrous ferric oxide (61); Cr(III) sorption to hydrous Fe oxides (62); Sr sorption to ferrihydrite (63); iron-cyanide adsorption onto γ-Al2O3 (64); Cu, Cd, and Pb adsorption on red mud (mixed Al, Fe, Si, and Ti oxides) (65); and Pb sorption to biogenic Mn oxides (66). Because the site densities (Ct in moles per gram) of HAO, HFO, and HMO are approximately equivalent to the total number of available sites (Cv in moles per gram) (15, 16), the Langmuir isotherm reduces to the linear distribution model. The distribution coefficient Kd (L g-1) is then simply the product of the equilibrium constant and site density (11, 15, 16, 44). Similar results were observed at 4 and 11 °C with R2 g 0.94. Metal adsorption to the hydrous oxides increased with increasing temperature, indicating an endothermic adsorption reaction (Figure 3). This phenomenon has been
observed by many others: Sr and Cd sorption to hydrous manganese oxide (36, 38); Sr and Cd sorption to HFO (15); Cd, Pb, and Zn sorption to goethite (33-35); and Sr sorption to HAO and HMO (16). Typically, most hydrogen bonds range between 2 and 12 kcal mol-1, covalent or chemical bonds are greater than 50 kcal mol-1, and physical types of bonds are much less than chemical ones (40, 67). Because adsorption enthalpies in the present studies do not exceed 25 kcal mol-1, adsorption is a physical type of reaction where the cations do not lose their waters of hydration upon sorption (40, 41). In the slow sorption process, where a CBC was maintained, the initial amount of metal sorbed corresponded to the isotherm results representing adsorption to the external surface (Figure 4, panels a and c). With time, the amount of metal ion sorbed to the oxide gradually increased due to intraparticle surface diffusion. From our isotherms, we found VOL. 34, NO. 11, 2000 / ENVIRONMENTAL SCIENCE & TECHNOLOGY
9
2219
FIGURE 6. Predicting equilibrium at pH 7 and 25 °C for Sr, Cd, and Zn sorption to hydrous metal oxides given their respective particle size distributions. no evidence of more than one average type of site. Assuming that there is no difference between the type of sites located externally and internally and given the particle geometry and boundary conditions, the solution to the partial differential equation from the mass balance is integrated over the volume of the spherical aggregated particles resulting in the mass sorbed per particle at a given time (11, 15, 16, 46):
M ) 4πCs
[
(
)]
-Dn2π2t 6 ∞ 1 1exp 3 π2n)1n2 R2
R3
∑
(3)
where 2220
9
ENVIRONMENTAL SCIENCE & TECHNOLOGY / VOL. 34, NO. 11, 2000
D)
Ds
[ ( )] 1+
FKi
(4)
where Ds is the surface diffusivity and fitting parameter in the model, is the oxide porosity, F is the bulk density, and Ki is the distribution coefficient representing the equilibrium constant times the internal site density. As mentioned above, the oxides have been characterized in previous work (11, 16). The amount sorbed to the internal surface of a single particle times the number of particles present for given radius (R) was summed over the entire particle size distribution to obtain the concentration sorbed internally. The particle size
FIGURE 7. (a) Correlation of adsorption enthalpy with the inverse of the hydrated radius [RH (71)], and (b) correlation between activation energy and adsorption enthalpy based on Polanyi relation. distributions used have been characterized for all three oxides in situ (11, 16). This total plus that sorbed to the external surface found from the short-term sorption studies (Kd) provides the theoretical concentration. By minimizing the variance, the only fitting parameter is surface diffusivity; modeling results are shown in Figures 4 and 5 (and Table 1) and include Cd sorption to HFO (15). Errors associated with the model from the propagation of errors (POE) statistical method (68) are also included in Figures 4 and 5. The POE accounted for the standard deviations (SD) in the number of particles for the 32 different sizes as well as the error in the distribution coefficient describing the mass adsorbed to the external surface. All data fall within 1 or 2 standard deviations of the model. Best-fit surface diffusivities ranged from 10-16 to 10-12 cm2 s-1. Barrow et al. (56) estimated diffusion parameters for Zn and Cd in goethite to be 6.5 × 10-19 cm2 s-1 and 3 × 10-20 cm2 s-1, respectively. Fuller et al. (10) evaluated effective
diffusivities for As(V) in HFO, which were in the order of 10-11 cm2 s-1. Papelis et al. (12) found Cd and Se(IV) diffusivities in porous aluminas ranged between 10-12 and 10-10 cm2 s-1. Axe and Anderson (11, 15) evaluated Cd and Sr surface diffusivities in HFO that ranged from 10-14 to 10-13 cm2 s-1. Misak et al. (13) reported self-diffusion coefficients of Co and Zn in hydrous Fe(III) and Sn(IV) oxides to be in the order of 10-11 cm2 s-1. In a more recent work, diffusivities for Sr in HAO and HMO ranged from 10-12 to 10-11 cm2 s-1 (16). Furthermore, in the research presented here along with previous work (15, 16), of the total adsorption sites, internal ones constitute as much as 50% in HAO, 40% in HFO, and 90% in HMO. These site densities are based on modeling the systems studied here and are in agreement with those observed when studying Sr sorption and assessing site densities as detailed by Trivedi and Axe (16). For each oxide and its associated particle size distribution, Figure 6 illustrates internal sorption or intraparticle surface diffusion (at pH 7) VOL. 34, NO. 11, 2000 / ENVIRONMENTAL SCIENCE & TECHNOLOGY
9
2221
as a function of the square root of time, where sorption reaches equilibrium in approximately 3-4 months for Sr and 2-5 years for Cd. Interestingly, Zn, which has the greatest affinity for these oxides, requires approximately 5-10 years to reach equilibrium. Thus, in aquatic environments where the oxides are present as discrete particles and coatings, amorphous aluminum, iron, and manganese oxides act as a sink for trace metal contaminants. Conversely, during desorption, these amorphous oxide particles act as longterm source for contaminants in ecosystems. Most transport models that employ distribution coefficients and retardation factors are inadequate for describing metal and radionuclide mobility in subsurface environments (69, 70). Because the internal sites account for as much as 90% of the total, accurate modeling requires inclusion of this contribution. However, long-term studies require a lengthy period of time, and therefore predictive methods would be useful for determining surface diffusivities. In earlier studies (15, 16), surface diffusivities decreased with an increase in pH, demonstrating that they are a function of site density, which is consistent with site activation theory (15-17). In sorption studies with strontium (16), Ea was approximately equivalent for all three oxides. In this research, the activation energy for a given metal was found again to be equivalent [within the error (68)] for all oxides as well. On the basis of the experimental surface diffusivity, the best-fit activation energy for each system was evaluated using eq 1. As a result, for a given oxide, R was approximately equivalent for the transition metals Cd and Zn (Table 1). In an effort to predict the activation energies, adsorption enthalpies were correlated with the inverse of the hydrated radius (RH) of the adsorbate (Figure 7a) (71). The linear correlation suggests that metal adsorption to hydrous oxides is dominated by Coulombic forces, as in physical interactions the hydrated sheath would be retained. The correlation presented potentially allows for predicting adsorption enthalpy (and hence the equilibrium constant) for any metal sorbed to hydrous Al, Fe, and Mn oxides. This correlation and recent spectroscopic studies of Sr sorption to HFO (72) and HMO (73) suggest that metals such as Sr, Cd, and Zn are electrostatically complexed with 1-2 oxygens on the oxide surface. In Figure 7b, the relationship between activation energy and enthalpy is presented for Cd and Zn adsorption to each oxide. Interestingly, as shown in the figure, the activation energy for a metal is equivalent for all three oxides. Furthermore, for each oxide, R is approximately equivalent (given the error) for Cd and Zn, which suggests that for a particular oxide, R is constant for metals of the same group from the Periodic Table. Therefore, the activation energy can be estimated from the ∆H° - RH-1 correlation. Theoretical surface diffusivities presented in Table 1 suggest that the slow sorption process can be modeled using site activation theory to predict diffusivities. These surface diffusivities reveal that the slow sorption process may require anywhere from a few months to several years to reach equilibrium. The contribution of this slow intraparticle sorption accounts for as much as 90% of the adsorption sites.
Acknowledgments This research was supported through start-up funds granted to L. Axe from New Jersey Institute of Technology.
Literature Cited (1) Hesterberg, D. Agric., Ecosys. Environ. 1998, 67, 121-123. (2) U.S. EPA. ATSDR Information Center, CERCLA Priority List of Hazardous Substances, SARA, 1997. (3) ATSDR, Division of Toxicology, Public Health Statement: Cadmium, 1989. (4) ATSDR, Division of Toxicology, Public Health Statement: Zinc, 1989. (5) Morgan, J. J.; Stumm, W. J. Colloid Sci. 1964, 19, 347-359. 2222
9
ENVIRONMENTAL SCIENCE & TECHNOLOGY / VOL. 34, NO. 11, 2000
(6) McKenzie, R. M. Aust. J. Soil. Res. 1980, 18, 61-73. (7) Dzombak, D. A.; Morel, F. M. M. Surface Complexation Modeling: Hydrous Ferric Oxide; John Wiley & Sons: New York, 1990. (8) Martinez, C. E.; McBride, M. B. Environ. Sci. Technol. 1999, 33, 745-750. (9) O’Day, P. A.; Carroll, S. A.; Waychunas, G. A. Environ. Sci. Technol. 1998, 32, 943-955. (10) Fuller, C. C.; Davis, J. A.; Waychunas, G. A. Geochim. Cosmochim. Acta 1993, 57, 2271-2282. (11) Axe, L.; Anderson, P. J. Colloid Interface Sci. 1995, 175, 157165. (12) Papelis, C.; Roberts, P. V.; Leckie, J. O. Environ. Sci. Technol. 1995, 29, 1099-1108. (13) Misak, N. Z.; Ghonemy, G. H.; Morcos, T. N. J. Colloid Interface Sci. 1996, 184, 31-43. (14) Strawn, D. G.; Scheidegger, A. M.; Sparks, D. L. Environ. Sci. Technol. 1998, 32, 2596-2601. (15) Axe, L.; Anderson, P. J. Colloid Interface Sci. 1997, 185, 436448. (16) Trivedi, P.; Axe, L. J. Colloid Interface Sci. 1999, 218, 554-563. (17) Ka¨rger, J.; Ruthven, D. M. Diffusion in Zeolites and Other Microporous Solids, Wiley-Interscience: New York, 1992. (18) Boudart, M. Kinetics of Chemical Processes, Prentice Hall: Englewood Cliffs, NJ, 1968. (19) Gadde, R. R.; Laitinen, H. A. Anal. Chem. 1974, 46 (13), 20222026. (20) Kinniburgh, D. G.; Jackson, M. L.; Syers, J. K. Soil Sci. Soc. Am. J. 1976, 40, 796-799. (21) Benjamin, M.; Leckie, J. O. J. Colloid Interface Sci. 1981, 79, 209-221. (22) Lion, L. W.; Altmann, R. S.; Leckie, J. O. Environ. Sci. Technol. 1982, 16, 660-666. (23) Okazaki, M.; Takamidoh, K.; Yamane, I. Soil Sci. Plant Nutr. 1986, 32 (4), 523-533. (24) Garnier, J.; Benyahya, L. Environ. Sci. Technol. 1999, 33, 13981407. (25) Mason, Y.; Ammann, A. A.; Ulrich, A.; Sigg. L. Environ. Sci. Technol. 1999, 33, 1588-1597. (26) Gray, M. J. Effluent Water Treat. J. 1981, 21, 201-203. (27) Pan, G.; Liss, P. S. J. Colloid Interface Sci. 1998, 201, 77-85. (28) Coughlin, B.; Stone, A. Environ. Sci. Technol. 1995, 25, 24452455. (29) Kraepiel, A. M. L.; Keller, K.; Morel, F. M. M. J. Colloid Interface Sci. 1999, 210, 43-54. (30) Lothenbach, B.; Furrer, G.; Schulin, R. Environ. Sci. Technol. 1997, 31, 1452-1462. (31) Tamura, H.; Furuichi, R. J. Colloid Interface Sci. 1997, 195, 241249. (32) Wehrli, B.; Ibric, S.; Stumm, W. Colloids Surf. 1990, 51, 77-88. (33) Johnson, B. Environ. Sci. Technol. 1990, 24, 112-118. (34) Rodda, D. P.; Johnson, B. B.; Wells, J. D. J. Colloid Interface Sci. 1996, 184, 365-377. (35) Angove, M. J.; Johnson, B. B.; Wells, J. D. J. Colloid Interface Sci. 1998, 204, 93-103. (36) Mishra, S. P.; Tiwary, D. Radiochim. Acta 1995, 69, 121-126. (37) Mishra, S. P.; Singh, V. K.; Tiwary, D. Radiochim. Acta 1997, 76, 97-101. (38) Mishra, S. P.; Tiwary, D. Radiochim. Acta 1998, 80, 213-217. (39) Smith, J. M.; Van Ness, H. C. Introduction to chemical engineering thermodynamics, 4th ed.; McGraw-Hill: New York, 1987. (40) de Boer, J. H. The Dynamical character of Adsorption, Clarendon Press: Oxford, U.K., 1968. (41) Israelachvilli, J. Intermolecular & Surface Forces, 2nd ed.; Academic Press: New York, 1997. (42) Dzombak, D. A.; Morel, F. M. M. J. Colloid Interface Sci. 1986, 112 (2), 588-598. (43) Eaton, A. D.; Clesceri, L. S.; Greenberg, A. E. Standard Methods for the Examination of Water and Wastewater; APHA, AWWA, and WEF: Washington DC, 1995. (44) Fogler, H. S. Elements of Chemical Reaction Engineering, 2nd ed.; Prentice Hall: Englewood Cliffs, NJ, 1992. (45) Axe, L.; Anderson, P. R. In Adsorption of metals by Geomedia: Variables, Mechanisms, and Model Applications; Jenne, E. A., Ed.; Academic Press: New York, 1998; Chapter 8, pp 193-209. (46) Crank, J. The Mathematics of Diffusion, 2nd ed.; Clarendon Press: Oxford, U.K., 1975. (47) Waychunas, G. A.; Rea, B. A.; Fuller, C. C.; David, J. A. Geochim. Cosmochim. Acta 1993, 57, 2251-2269. (48) Jenne, E. A. In Trace Inorganics in Water; Gould, R. F., Ed.; ACS Advances in Chemistry, 1968; Vol. 73, pp 337-387.
(49) Baltpurvins, K. A.; Burns, R. C.; Lawrance, G. A.; Stuart, A. D. Environ. Sci. Technol. 1996, 30 (3), 939-944. (50) Shuman, L. M. Soil. Sci. Soc. Am. J. 1977, 41, 703-706. (51) Anderson, P.; Benjamin, M. M. Environ. Sci. Technol. 1990, 24 (5), 692-698. (52) Meng, X.; Letterman, R. D. Environ. Sci. Technol. 1993, 27, 970975. (53) Crawford, R. J.; Harding, I. H.; Mainwaring, D. E. Langmuir 1990, 9, 3050-3056. (54) Tamura, H.; Katayama, N.; Furuichi, R. Environ. Sci. Technol. 1996, 30, 1198-1204. (55) Lu ¨ tzenkirchen, J. J. Colloid Interface Sci. 1997, 195, 149-155. (56) Barrow, N. J.; Gerth, J.; Bru ¨ mmer, G. W. J. Soil Sci. 1989, 40, 437-450. (57) Leckie, J. O.; Benjamin, M. M.; Hayes, K. F.; Kaufman, G.; Altmann, S. Adsorption/Coprecipitation of Trace Metals from Water with Iron Oxyhydroxide; Electric Power Research Institute, CA, 1980; EPRI RP-910-1. (58) Cowan, C. E.; Zachara, J. M.; Resch, C. T. Environ. Sci. Technol. 1991, 25, 437-446. (59) Christophi, C. A.; Axe, L. J. Environ. Eng. 2000, 126 (1), 66-74. (60) Posselt, H. S.; Weber, W. J., Jr. In Symposium on the Chemistry of Water Supply, Treatment, and Distribution; Rubin, A. J., Ed.; Ann Arbor Science: Ann Arbor, MI, 1974; pp 89-108. (61) Gadde, R. R.; Laitinen, H. A. Environ. Lett. 1976, 5 (4), 223-235. (62) Charlet, L.; Manceau, A. A. J. Colloid Interface Sci. 1992, 148 (2), 443-458.
(63) C ˇ ernik, M.; Borkovec, M.; Westall, J. C. Langmuir 1996, 12, 61276137. (64) Huang, C.; Cheng, W. P. J. Colloid Interface Sci. 1997, 188, 270274. (65) Apak, R.; Gu ¨ c¸ lu ¨ , K.; Turgut, M. H. J. Colloid Interface Sci. 1998, 203, 122-130. (66) Nelson, Y. M.; Lion, L. W.; Ghiorse, W. C.; Shuler, M. L. Appl. Environ. Microbiol. 1999, 65 (1), 175-180. (67) Prausnitz, J. M.; Lichtenthaler, R. N.; de Azevedo, E. G. Molecular Thermodynamics of Fluid-Phase Equilibria, 2nd ed.; Prentice Hall International Series: Englewood Cliffs, NJ, 1986. (68) Ku, H. H. J. Res. Natl. Bur. Stand., Sect. C 1966, 70C (4), 263273. (69) Mahara, Y. J Environ. Qual. 1993, 22, 722-730. (70) Ohtsuka, Y.; Takebe, S. J. Nucl. Sci. Tech. 1990, 27, 62. (71) Nightingale, E. R., Jr. J. Can. Phy. Chem. 1959, 63, 1381-1387. (72) Axe, L.; Bunker, G. B.; Anderson, P.; Tyson, T. J. Colloid Interface Sci. 1998, 199, 44-52. (73) Axe, L.; Tyson, T.; Trivedi, P.; Morrison, T. J. Colloid Interface Sci. 2000, 224, 408-416.
Received for review September 28, 1999. Revised manuscript received March 13, 2000. Accepted March 15, 2000. ES991110C
VOL. 34, NO. 11, 2000 / ENVIRONMENTAL SCIENCE & TECHNOLOGY
9
2223