This is to certify that the dissertation entitled PARTITIONING 0F ZINC BETWEEN SOLID AND LIQUID IN SYNTHETIC %OIL MATERIALS presen ed y RICHARD CLEVELAND ZIELKE has been accepted towards fulfillinent of the requirements for Ph . D . SOIL CHEMISTRY degree in '/ Major professor Date 11/30/86 MS U is an Affirmative Action/Equal Opportunity Institution 0-12771 MSU LIBRARIES —:__. RETURNING MATERIALS: PIace in book drop to remove this checkout from your record. FINES wiII be charged if book is returned after the date stamped beIow. PARTITIONING OF ZINC BETWEEN SOLID AND LIQUID IN SYNTHETIC AND SOIL MATERIALS By Richard Cleveland Zielke A DISSERTATION Submitted to Michigan State University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Department of Crop and Soil Sciences 1986 “r or» L. a» s» 1 ABSTRACT PARTITIONING OF ZINC BETWEEN SOLID AND LIQUID IN SYNTHETIC AND SOIL MATERIALS By Richard Cleveland Zielke Metal activities in soil solutions and natural waters are primarily a function of partitioning ‘between colloidal and aqueous phases. For Zn2+ , solution activity is likely controlled by adsorption- desorption rather than precipitation-dissolution reactions. Frequently adsorption is irreversible and the resulting solution activity is much lower than analytical detection limits. A batch equilibrium technique based on isotopic dilution and mass balance was developed to study adsorption and desorption of trace Zn2+ by Ca smectite (montmorillonite # 25, Upton, Wyoming). The method was verified analytically and hence applicable to the study of desorption in clay minerals. Desorption hysteresis commonly occurred when <1% of the CEC was occupied by Zn. Three isotherms (Linear, Freundlich and Langmuir) were tested to determine their utility in characterizing Zn desorption. The Freundlich equation proved most useful and enabled adsorption to be viewed over 2+ several decades of Zn activity. Comparison of model ion exchangers (Ca. Dowex BOW-X8, Ca. Chelex-lOO) wit}: Ca. smectite and. Charity" clay Richard Cleveland Zielke (aeric. 'haplaquept) allowed. insight. to mechanisms operating in adsorption, desorption and hysteresis. For smectite, reversibility was attributed to an ion exchange mechanism while hysteresis was attributed to strong bonding at broken edge OH. Hysteresis was ubiquitous in Charity clay and was attributed to broken edge OH, A1 and Fe oxides, free carbonates and organic matter. The pH dependency of hysteresis at broken edge OH was investigated. Only at pH 3.5 was adsorption by Ca smectite completely reversible. As the pH increased, adsorption and hysteresis increased, however reversibility at high surface occupancy of Zn verified that smectites can mask precipitation in alkaline solution. The effect of initial exchange ion on hysteresis was also investigated using homoionic Na, K, Rb and Cs smectites. Adsorption followed the lyotropic series, however adsorption and desorption by Rb and Cs smectites was sterically inhibited at interlayer exchange sites because Rb+ and Cs+ have a demonstrated ability to reduce swelling of smectite minerals thus reducing the rate of Zn diffusion into and out of the interlayer positions. Hysteresis in this case was merely an artifact of nonequilibrium. DEDICATION TO MY MOTHER, BROTHERS, THEIR WIVES AND CHILDREN ii ACKNOWLEDGEMENTS Deepest appreciation is extended to my nmjor professor Dr. Boyd Ellis for his invaluable assistance in preparing this manuscript and the friendship provided throughout the course of study. Also gratefully acknowledged are Drs. M.M. Mortland, S.A. Boyd and T.J. Pinnavaia who served on the graduate committee. The friendship, expertise and unselfish dedication of Calvin Bricker was also invaluable throughout the course of my studies and hence warrants special thanks. Appreciation is also extended to the People of the State of Michigan who through taxation support the Michigan Agriculture Experiment Station which provided financial support for this project. iii TABLE OF CONTENTS LIST OF TABLES LIST OF FIGURES. INTRODUCTION . CHAPTER I. II. III. IV. VI. LITERATURE REVIEW . ADSORPTION AND DESORPTION OF ZINC BY SMECTITE: COMPARISON OF RADIOTRACER AND ANALYTICAL METHODS . Introduction . Materials and Methods. Results and Discussion . Literature Cited . DESORPTION OF ZINC FROM RESIN, SMECTITE AND CHARITY CLAY . Introduction . Materials and Methods. Results and Discussion . Literature Cited . EFFECT OF pH ON ADSORPTION AND DESORPTION OF ZINC BY SMECTITE. . . . . . . Introduction . Materials and Methods. Results and Discussion . Literature Cited . EFFECT OF INITIAL EXCHANGEABLE ION ON ADSORPTION AND DESORPTION OF ZINC BY SMECTITE . Introduction . Materials and Methods. Results and Discussion . Literature Cited . SUMMARY AND CONCLUSIONS. iv . vi . vii l6 l6 18 2O 29 31 31 33 36 61 62 62 64 66 79 8O 8O 81 83 97 98 APPENDIX. . . . . . . . . . . . . . . . . . . . 102 LIST OF REFERENCES. . . . . . . . . . . . . . . 111 LIST OF TABLES TABLE CHAPTER II 1. Comparison of Freundlich and Linear isotherm data at the lowest Zn level in the tracer study. CHAPTER III 1. Chemical properties and mass of adsorbent used in batch adsorption experiments . 2. Summary+ of linear regression equations for Zn adsorption by Ca chelex, Ca Dowex, Ca smectite and charity clay . 3. Values of D 1 for adsorption of Zn2+ on Ca Dowex, Cg Chelex, Ca smectite and charity clay . 4. Equilibrium pH and Zn activity for Ca Chelex adsorption isotherms . . . . . . . . CHAPTER IV 1. Initial Zn2+ levels and equilibrium pH's for Zn adsorption by Ca smectite . . . . . CHAPTER V 1. Initial Zn2+ concentrations and equilibrium pH's for the adsorption of Zn by homoionic smectites. . . 2. Effect of interlayer cation on basal d(001) spacings of Wyoming montmorillonite immersed in water . vi PAGE 27 34 37 . 43 . 49 . 66 . 83 . 9O LIST OF FIGURES FIGURE PAGE CHAPTER II 1. Comparative adsorption of Zn by Ca smectite using analytical and radiotracer methods . . . 21 2. Desorption of Zn from Ca smectite as measured by analytical methods . . . . . . . . 22 3. Desorption of Zn from Ca smectite as measured by radiotracer methods. . . . . . . . 23 4. An illustration of the Linear equation data variability for adsorption and desorption of Zn as measured by radiotracer methods. . . . . 26 CHAPTER III 1. Adsorption and desorption of Zn by Ca Dowex ion exchange resin fit to the Linear equation. 39 2. Adsorption and desorption of Zn from Ca Dowex ion exchange resin fit to the Freundlich equation . . . . . . . . . . . . . . . . . . . 4O 3. Adsorption and desorption of Zn by Dowex ion exchange resin fit to the Langmuir equation. . 41 4. Desorption of trace amounts of Zn by Ca Dowex ion exchange resin as measured by the Linear equation . . . . . . . . . . . . . . . . . . . 42 5. Adsorption and desorption of Zn by Ca Chelex ion exchange resin fit to the Linear equation. 45 6. Desorption of trace amounts of Zn by Ca Chelex ion exchange resin as measured by the Linear equation . . . . . . . . . . . . . . . . . . . 46 7. Adsorption and desorption of Zn from Ca Chelex ion exchange resin fit to the Freundlich equation. . . . . . . . . . . . . . . . . . . . 47 vii Adsorption and desorption of Zn by Ca Chelex ion exchange resin fit to the Langmuir Adsorption and desorption of Zn by Ca smectite at pH 6.6 fit to the Linear equation. Adsorption and desorption of Zn by Ca smectite at pH 6.6 fit to the Freundlich equation. Adsorption and desorption of Zn by Ca smectite at pH 6.6 fit to the Langmuir equation. Desorption of trace amounts of Zn by Ca smectite pH 6.5 as measured by the Linear Adsorption of Zn by Ca smectite desorption fit to the desorption fit to the desorption fit to the desorption pH 3.5. desorption pH 4.5. desorption pH 5.9. desorption pH 6.5. desorption pH 7.0 desorption pH 7.7 of Zn by a Charity Linear equation . of Zn by a Charity Freundlich equation . of Zn by a Charity Langmuir equation. of Zn by of Zn by of Zn by of Zn by of Zn of Zn by as a function Effect of time on the adsorption of Zn by pH 4.5 Ca smectite at three initial Zn levels . equation 9. 10. 11. 12. equation 13. Adsorption and clay at pH 7.4 14. Adsorption and clay at pH 7.4 15. Adsorption and clay at pH 7.4 CHAPTER IV 1. Adsorption and Ca Smectite at 2. Adsorption and Ca smectite at 3. Adsorption and Ca smectite at 4. Adsorption and Ca smectite at 5. Adsorption and Ca smectite at 6. Adsorption and Ca smectite at 7. of time and pH . 8. 9. Effect of time on adsorption and desorption of Zn by Ca smectite at pH 4.5. viii . 48 . 51 52 53 54 56 57 58 67 69 70 71 72 73 75 76 . 77 CHAPTER V Adsorption Adsorption Adsorption Adsorption Adsorption smectites as a function of time . Effect of time on smectite at three Effect of time on smectite at three Effect of time on smectite at three Effect of time on smectite at three the adsorption of initial Zn levels . the adsorption of initial Zn levels . the adsorption of initial Zn levels . the adsorption of initial Zn Levels . ix and desorption of Zn by Na smectite. and desorption of Zn by K smectite and desorption of Zn by Rb smectite. and desorption of Zn by Cs smectite. of Zn by alkali metal saturated Zn by Na Zn by K Zn by Rb Zn by Cs 84 . 86 87 88 . 89 . 91 . 92 . 93 . 94 INTRODUCTION Zinc retention by soils and soil materials has been the subject of extensive research due to the fact that Zn is both an essential micronutrient which has been deficient in many soils and a potential environmental hazard” Several soil components have been isolated and shown to effect removal of Zn from soil solution. These include Fe, A1 and Mn oxides, CaCO3, humic materials and clay minerals. In order to elucidate the mechanism. of adsorption, studies using synthetic and natural components ensued; however this approach has been criticized as being environmentally unrealistic owing to the complex interaction of soil colloids. Ironically, it was the comparison of single components and whole soils that revealed these interactions. Recently, there has been a resurgence of adsorption studies utilizing whole soils and/or soils which are treated to selectively remove individual fractions associated with Zn retention. In light of the heterogenous nature of soil, it is not surprising to find conflicting evidence for adsorption and precipitation. The interrelationship ‘between surface area, structure and reactivity of clays with metal emphasizes the importance of pure clay mineral studies. The role of organic matter in adsorption, chelation and mobilization of metals has also been demonstrated. The observation that metal retention by some soils closely resembles the behavior exhibited by amorphous Fe and Al oxides, and the subsequent role of surface hydroxyl groups in "specific adsorption" has furthered the debate over the mechanism of Zn 1 retention in soils. Indeed, progress in this area of research can be likened to Jenny's (1961) soil acidity merry-go-round in which each turn produces a different though not unfamiliar view of what was termed the soil chemistry landscape. There can be little doubt that no single component controls Zn retention and release in all environmental systems. Certainly, the often reported failure of solubility principles to characterize Zn2+ reactivity in soil is reason enough to pursue an adsorption approach. Reports that Zn is mainly associated with the mineral fraction in many soils and sediments, combined with the recent interest in clay minerals as sorbents for metallic and radioactive waste is the basis for selecting smectite clay as the adsorbent phase in the present study. Although adsorption of Zn by smectite minerals has been studied and numerous fractionation and extraction schemes have been proposed, few, if any contemporary studies have focused on true desorption over a wide range of surface compositions. Many such studies have failed or neglected to control pH or ionic strength, or have attempted desorption using unrealistic concentrations of acids, competing ions or chelating agents. The objective of this study was to characterize Zn desorption in response to pH and surface composition variables. More specifically, the objectives were to (i) verify the applicability of the principles of isotopic dilution to adsorption of trace amounts of Zn2+ by smectite, (ii) determine the applicability of selected adsorption isotherms in characterizing desorption phenomena for some model exchangers, (iii) determine the effect of pH on adsorption and desorption of Zn2+ by Ca- smectite, and (iv) determine the effect of initial exchange ion on Zn2+ adsorption and desorption by smectite. The results of these basic studies will hopefully have application to waste disposal engineering as well as conventional agriculture. LITERATURE REVIEW Mechanistically, removal of Zn2+ from solution by soils and clays affords extreme ambiguity. In absence of microrganisms, living plant roots and soluble organic ligands, two general processes are possible; precipitation and adsorption. For the .purpose of this discussion, adsorption is the general process whereby solute adheres to the surface of a soil particle. Chemisorption. is adsorption resulting from a specific chemical reaction, hence the term "specific adsorption". Ion exchange is the adsorption process by which charged solutes exchange with ions on a soil particle. Precipitation is a related process in that it is generally preceded by adsorption; the difference being that adsorption is a two dimensional process, while precipitation is a three dimensional process (Corey, 1981). The reactivity of solutes with soils and clays is generally characterized by an adsorption isotherm. Because of their empirical nature, adsorption isotherms fail to distinguish between the totality of processes resulting in the attainment of equilibrium, and hence cannot be used to distinguish the sorption mechanism (Veith and Sposito, 1977). A detailed description of various adsorption isotherms and their use is given by Travis and Etnier (1981). Of these, three have found extensive use in describing the equilibrium between reactive solutes and soil components. These are the Linear, Freundlich and Langmuir isotherms. The Linear adsorption isotherm. is an empirical equation that relates the quantity of solute adsorbed (X) per unit mass of adsorber (m) to an equilibrium solution concentration or activity (C). Thus, X/m - kC where k is the distribution coefficient. Few adsorption processes in soil obey this equation, but it has been used in conjunction with transport equations to describe radionuclide transport from geologic repositories (Burkholder, 1976) and 2+ in soil columns (Lai et al., 1978). adsorption of Na+ and .Mg Frequently the observed isotherm is nonlinear (Frost and Griffin, 1977; Maquire et al., 1981; Brummer et al., 1983; Harter, 1983), depicting a rapid rise in adsorption followed by a bending over or tailing off to a constant value. One consequence of the asymptotic behavior at high equilibrium solute levels is that it provides first-hand observation of the adsorption capacity. Also the isotherm allows direct viewing of the relationship between solute and solid. McLaren et a1. (1983) used this fact to characterize Cu2+ desorption hysteresis from humic acid, Fe and Mn oxides and montmorillonite. Several researchers have described Zn adsorption using the nonlinear equation proposed by Freundlich (1926), X/m - KCl/n where n and K are empirical constants. The equation may be linearized by taking the log of both sides: log(X/m) - l/n log(C) + log(K) Thus a plot of log(X/m) vs log(C) yields a straight line if the data conform to the Freundlich equation. The equation has one advantage in that it can be used to view adsorption over several orders of magnitude; however, like the Linear isotherm it fails to predict the adsorption capacity. Furthermore, log-log plots tend to minimize data variability and the two constants allow for easy curve fitting. The Freundlich isotherm has been used to characterize desorption hysteresis. Barney (1984) used the ratio of the Freundlich exponents Ns/Nd (where Ns and Nd are the measured exponents for sorption and desorption, respectively) as a measure of the magnitude of hysteresis in radionuclide desorption from Columbia River basalt materials. Elrashidi and O'Connor (1982) used the Freundlich equation to describe Zn sorption and desorption from nine soils of varying physical and chemical properties. Desorption of Zn from all soils was extremely hysteretic even after five desorption cycles using .01 N CaClZ. The Langmuir equation (1918), unlike the others, is derived from the kinetic theory of gases. It's adaptation to solid-liquid adsorption phenomena in soils was proposed by Boyd et a1. (1947). Of the many forms derived, two have found frequent usage: C/X/m - l/Kb + C/b X/m - KbC/(1+KC) where K is a measure of the bond strength, b is the adsorption maximum and C is the ion activity. In applying the Langmuir equation monolayer adsorption and uniform bonding energy are assumed. Based on the heterogeny of soils and clays, seldom are these assumptions met. NOnlinearity has been attributed to sites of differing bonding energy, hence two surface isotherms have been proposed (Shuman, 1975, 1977). This approach has been the subject of recent controversy (Sposito, 1982; Posner and Bowden, 1980), and it was suggested. that independent evidence for sites of different bonding energy be provided before attempting to use the two surface equation. Curvilinear isotherms have also been attributed to the a lack of consideration of desorbed ions in exchange processes (Griffin and Au, 1977; and Harter and Baker, 1977). Since extreme variability can result from the division of the measured quantities (C) and X/m, no attempts to use the Langmuir equation to describe desorption have been reported. Any discussion of Zn2+ adsorption by smectites should logically be preceded by a brief description of their chemical nature in relation to their properties. Smectites are layer lattice silicates formed through weathering of igneous rocks such as granite and basalt. Their behavior and abundance in soils and sediments make smectites an important class of minerals from both an agricultural and environmental standpoint. Among the most noted properties affecting ion adsorption are the tremendous surface area (8 x 105 m2 kg'l) and high CEC (80-150 cmol(+) kg'l). Smectites are comprised of a central alumina octahedra sandwiched between two silica tetrahedra. The tips of the two tetrahedra are oriented such that apical O are shared with O of the alumina octahedra. The result is a continuous two dimensional network extending in the a and b dimension with individual 2:1 units stacked in the c dimension (Grim, 1968). Depending upon the layer charge, 0 atoms of adjacent structural 'units form.'bonds which. may 'be cleaved to allow lattice expansion. The magnitude of the surface charge and occurrence of exchangeable cations between structural units, in accordance with the principal of electrical neutrality, is the result of localized electron density due to substitution of A13+ for Si4+ in the tetrahedral and Mg2+ and Fe2+ for A13+ in the octahedral layers. This type of surface charge is termed permanent CEC. Fbr the dioctahedral mineral montmorillonite #25 used in this study, the CEC (90 cmol(+) kg'l) is mostly attributed to octahedral substitution while the swelling properties are attributed to the 'moderate surface charge. and. hydrated. nature of exchangeable cations. In addition to lattice substitution, cation exchange sites are thought to arise from the existence of broken bonds in the a-b dimension. Such broken 0 bonds and edge OH's are considered to undergo ionization and exchange under alkaline conditions and hence are termed pH dependent sites (Borchardt, 1977). Helfferich (1962) generalized ion exchange on soil colloids as being an electrostatic phenomenon which is diffusion controlled, stoichiometric and reversible. There are several reports of non-reversible exchange of cations in pure mineral studies involving both neutral salts and transition metals. The inability to attain the same position of equilibrium when approached from either direction has been termed hysteresis or fixation. Vanselow (1932) observed hysteresis in Ca2+-NH4+ exchange on bentonite. Wiegner (1935) attributed irreversibility to the existence of sites having different geometries and. bonding energies. Based on electrostatic principles, broken edge OH's, octahedral substitution and tetrahedral substitution should provide for a range of sites with different bonding energy. Some researchers have attributed hysteresis to the imperfect dispersion of flocculated clays and resulting inaccessibility of exchange sites (Van Bladel and Laudelout, 1967; Tabikh et al., 1960). Marshall (1964) attributed hysteresis to the possibility of attaining minimum free energy at more than one substrate composition. More recently, Maes and Cremers (1975) considered hysteresis in Zn2+ exchange with smectite to be composition and pH dependent. Adsorption was completely reversible at low surface occupancy and pH, but was neither stoichiometric nor reversible at high surface occupancy. The reaction could be made reversible by decreasing the pH. Due to the apparent role of OH', hysteresis was attributed to specific bonding. The validity of this work must be questioned since specific bonding , if it occurs, is considered a high affinity phenomenon and hence is expected at low surface occupancy. The nonreversible sorption of Zn by bentonite was reported by Elgabaly (1950). He proposed that Zn not removed by 1 M neutral NH4(OAC) became fixed in alumina octahedra resulting in a permanent loss of CEC. DeMumbrum and Jackson (1956) furthered this concept by proposing specific adsorption. of Zn. by montmorillonite. Based on perturbations of OH in IR studies, Zn retained in excess of the CEC was thought to be associated with octahedral OH groups. Tiller and Hodgson 2+ (1962) proposed that Zn specifically adsorbed by nmmtmorillonite in 10 the presence of excess electrolyte (thereby‘ eliminating 'nonspecific electrostatic sorption) was of two forms. One species exchangeable with Cu2+, N12+, Mn2+ and Fe2+ was considered to be associated with broken edge OH's, while the smaller nonexchangeable fraction was the result of lattice penetration. Specific adsorption (of Cu2+) on smectite has been challenged by McBride and Mortland (1974). Using IR, XRD, ESR and thermal methods, Cu2+ was shown to enter the hexagonal holes of the silicate structure upon dehydration. Reduction in the CEC was attributed to further penetration into the unoccupied hole in the octahedral layer. Upon Cu2+ rehydration, in hexagonal positions was free to move, but octahedral Cu2+ remained fixed, resulting in permanent loss of layer charge. Since the temperature used by DeMumbrum and Jackson (1956) was not sufficient for total dehydration of Cu2+, lattice penetration could not have occurred. It was suggested that excess Cu retention was the result of incomplete removal of' Cu(OAC)2 ‘While the IR. shift was attributed to formation of K smectite during the preparation of a KBr pellet. As a result, specific adsorption via this mechanism at temperatures < 150 C seemed unlikely. The reaction of silicic acid with montmorillonite is also thought to result in the formation of specific adsorption sites for transition metals (Hingston et al.,l972), although this interaction is considered to be a surface adsorbed layer rather than lattice penetration (Tiller, 1967). 2+ Several factors have been shown to affect adsorption of Zn and other metals by smectite. Nelson and Melsted (1955) observed that the 11 2+ Ca2+ relative adsorption of different cations in clays was H+ > Zn > >Mg2+ > KI, and suggested that retention of applied Zn may be affected by ions on the sorption complex. It is likely that the H+ was really Al3+ in their experiments due to clay lattice breakdown. Garcia- Miragaya and Page (1977) determined the effect of initial exchangeable ion on sorption of Cd by montmorillonite. Sorption was related to ion size and polarizability factors and decreased along the series Na+ > K+ > Ca2+ >A13+. Increasing ionic strength has also been shown to reduce metal sorption through complexation with inorganic ligands. In addition, solution composition was shown to affect adsorption since some ligands (eg. 8042') have a larger tendency to form ion pairs in solution (Garcia-Miragaya and Page, 1976; Mattigod et al., 1979). Mattigod and Sposito (1977) compiled. experimental association constants for complexation of Zn2+ with several inorganic ligands. Such data allow prediction of solution species available for adsorption or transport based on the ionic strength and solution composition. Since the mechanism of adsorption cannot be determined from adsorption isotherms, indirect methods have been used. The fact that Zn2+ adsorption was complete in 15 minutes suggested a ion exchange mechanism (Bourg and Filby, 1974). In a later study this hypothesis was 2+ in solution was in kinetic confirmed by demonstrating that Zn equilibrium with Zn adsorbed on the surface of montmorillonite (Bourg and Filby, 1976). Bingham et a1. (1964) also suggested an ion exchange mechanism since Zn retained by H-montmorillonite did not exceed the CEC if the solution pH was low enough to prevent Zn(0H)2 formation. 12 Singhal and Kumar (1977) used the exchange isotherm approach of Gaines and Thomas (1953) to calculate thermodynamic parameters for Zn2+ exchange on Mg-bentonite. Zinc exchange was found to be endothermic, however the use of thermodynamic approaches has been questioned since exchange is not always stoichiometric or reversible for metals which form hydroxy salts (Maes and Cremers, 1975). Farrah and Pickering (1978) used selective extraction procedures to characterize nonexchangeable Zn on montmorillonite, illite and kaolinite. Solvents were classified as acid, protonated, ligand and competing species, and their efficiency at removing Zn2+ previously sorbed at pH 5 and 7 was evaluated. While EDTA could completely remove sorbed Zn at both pH's, competing species exhibited a pH dependence. At pH 5 where exchange mechanisms have been proposed, sorbed Zn could be removed completely with 0.1 M NaNO3 or 0.05 M CaC12. At more dilute concentrations extraction was incomplete, but the fact that some reagents resulted in total recovery is evidence of minimum irreversibility.‘ The pH dependence of extraction suggests that clays are polyfunctional, containing mostly traditional exchange sites with some sites exhibiting affinities for protons and metal ions. The pH dependency of Zn2+ reactions in clays has received considerable attention. While clays appear to act as nucleation centers for Zn-hydroxy species (Farrah and Pickering,1976), there is no known physical evidence for a precipitated solid phase. Much of the evidence is based on solubility principles wherein thermodynamic data is either unavailable or inaccurate due to analytical limitations. 13 In pure clay studies marked uptake of Zn2+ increased when the pH exceeded the formation constant of Zn(OH)2 (Farrah and Pickering, 1977. In other studies adsorption in excess of the CEC has been attributed to the formation of ZnOH+ in bulk solution (Hodgson et al., 1964), but based on calculations of the Zn2+/ZnOH+ ratio, this explanation assumes preference for the 'hydrolyzed. species 'which seems unlikely from an electrostatic viewpoint (Bingham et al., 1964). Frost and Griffin (1977) suggested that Zn2+ adsorption from landfill leachate occurred by ion exchange at pH's < 6, while precipitation as Zn(OH)2 resulted at higher pH's. Farrah and Pickering (1976) observed precipitation of Zn(OH)2 in pure solution at pH's > 6.5 followed by dissolution and formation of the zincate anion at pH's > 10.5. In the presence of montmorillonite the pH required for precipitation was reduced by 0.5 units. Furthermore, since deprotonation and dissociation did not occur at high pH, chemisorption, not simple precipitation was implied. It should be noted that the conclusions in both these studies were based solely on comparisons of isotherms obtained in the presence and absence of clays in solution. Up to this point the discussion has focused on reactions of Zn2+ involving pure clay minerals. Even in such homogeneous materials there is a general lack of distinction between adsorption and precipitation processes. Furthermore the occurrence of hysteresis, specific bonding and the debate over the role of surface OH's in Zn retention has not been fully resolved. The behavior of Zn2+ in the presence of whole soils and soil clays, though not completely different, is unique enough to present a brief review for contrast and comparison. More extensive 14 reviews are given by Lindsay (1972), Nicholas and Egan (1975) and Shuman (1980). Zinc retention in soils has been attributed to a variety of components besides clay minerals. Some studies have reported significant correlations with the CEC (Shuman, 1975), while others have not (Harter, 1983). Many researchers have emphasized the importance of organic matter (Tyler and McBride, 1982; Kuo and Baker, 1980). Retention of Zn2+ by acid subsoils indicates that Al, Fe and Mn oxides may control Zn adsorption (Kalbasi et al., 1979; Shuman, 1976, 1977; Kinniburg and Jackson, 1982; McBride and Blasiac, 1979). Adsorption of Zn on Ca and Mg—carbonates has also been reported (Jurniak and Bauer, 1956). Tiller et a1. (1984) fractionated several soils and identified different magnitudes of influence for Zn adsorption by organic matter, A1 and Fe oxides and clays. All these materials are presumed to sorb Zn2+ through different mechanisms (Cavallaro and McBride, 1984). Regardless of the mechanism of adsorption, soils exhibit a definite pH effect (Kurdi and Donner, 1983; Bar-Yosef, 1979; Harter, 1983). There is general agreement among researchers that Zn retention in acid soils occurs on a variety of soil components and that Zn2+ in the soil solution is controlled by adsorption-desorption reactions. Conversely, there is sufficient evidence for Zn precipitation in alkaline soils. Among the various Zn2+ controlling ’solid phases suggested are Zn-carbonate (Udo et al., 1970), Zn-phosphate (Kalbasi et al., 1978), Zn-silicate (Tiller and Pickering, 1974), ZnS (Kittrick, 1976) and Zn(OH)2 (Saeed and Fox, 1977). Lindsay (1979) reported that 15 2+ many of these compounds were too soluble to control Zn activity, and proposed the existence of the unknown compound "soil Zn". While precipitation of Zn undoubtedly occurs in some soils, distinction between the adsorption-precipitation boundary is difficult. Recently Brummer et a1. (1983) investigated this problem, and concluded that Zn2+ in equilibrium with whole soils below pH 7 was controlled exclusively by adsorption-desorption processes. In neutral to alkaline soils precipitation is likely only if sufficient Zn2+ is added to ensure complete saturation of the adsorption complex. Such observations emphasize the importance of clays, A1 and Fe oxides, carbonates, humus and soluble ligands in masking precipitation reactions in soils, and further imputes a need to understand adsorption and particularly desorption phenomenon in soils and clays. CHAPTER II ADSORPTION AND DESORPTION OF ZINC BY SMECTITE: COMPARISON OF RADIOTRACER AND ANALYTICAL METHODS Introduction Much of the evidence for removal of Zn2+ and other heavy metals from solutions by soils and clay is expressed in the form of an equilibrium adsorption isotherm. Failure to attain the same position of equilibrium when approached from the opposite direction is termed hysteresis, and there is extensive chemical evidence for hysteresis when heavy metals are adsorbed by soils (Kuo and Mikkelsen, 1979; Elrashidi and O'Connor, 1982) and pure clay minerals (Tabikh et al., 1960; Van Bladel and Laudelout, 1967; McLaren et al., 1973: Maes and Cremers, 1975). Because of the large surface area and high CEC, smectite clays are considered to contribute to Zn retention in many soils. Among the mechanisms proposed for Zn retention by smectites are ion exchange (Bourg and Filby, 1976; Bingham et al., 1964), specific adsorption (Brummer et al., 1983), and precipitation (Frost and Griffin, 1977). Generally, cation exchange is suggested if the reaction is rapid, stoichiometric and most importantly, reversible (Helfferich, 1962). Hysteresis is usually associated with non-reversibility in the low end 16 17 of the adsorption isotherm where precipitation seems unlikely (Tiller and Hodgson, 1960; Tiller et al., 1984). As a result such hysteresis has been termed "specific adsorption". Demumbrum and Jackson (1956) proposed that specific adsorption of 2+ 2+ Cu and Zn was the result of lattice penetration, but McBride and Mortland (1974) have demonstrated using IR, ESR and thermal methods that temperatures exceeding 150 C were required for this process. While these physical methods 'have proven useful in studying the surface chemistry of clays, they are not sensitive enough or not applicable to the study of Zn2+ at low surface occupancy where hysteresis is often observed. Similarly, the ability of clays to effect almost complete removal of metal ions from solution poses some unique analytical problems to characterization of hysteresis using chemical methods. In the present study 652n was used as a tracer to enable characterization of adsorption and desorption of Zn2+ by smectite in the trace region. In light of reports of irreversible sorption, precipitation and specific adsorption, there is no guarantee that isotopic equilibria will occur at all substrate levels. The objectives of this study were to determine if the principles of isotopic dilution are valid over a wide range of surface compositions, and to develop techniques to further study adsorption and desorption of trace amounts of Zn2+ by Ca smectite. To accomplish these objectives, batch adsorption experiments using 65Zn as a tracer were compared to results 2+ obtained over an analytical range in which Zn could be determined by standard AA methods. 18 Materials and Methods Preparation of Ca,smectite A given weight of smectite (Montmorillonite # 25, Upton, Wyoming) was dispersed in distilled, deionized H20 by stirring for 24 hr. The colloidal portion was obtained by repeated sedimentation and decantation of the < 2uM fraction. Ca smectite was prepared by equilibrating this fraction with 1 M Ca012 (pH 6.0) for a period of 24 hr. At the end of 5 equilibrations the Ca clay was dialyzed against distilled, deionized H20 until a negative test for 01' was obtained using AgNO3. The resulting clay was suspended in 0.01 M CaC12 and stored until ready for use. The CEC of the clay was determined to be 90 i 2 cmol (+) kg'1 using the NHAOAC-NaOAC method of Jackson (1958), and the pH was measured to be 6.5 i 0.2 using a standard reference electrode. XRD was conducted but revealed no other detectable minerals. Adsorption isotherms Batch adsorption experiments were initiated by pipetting 25 mL of the rapidly stirred clay suspension into tared 50 mL Oak Ridge centrifuge tubes. The clays were centrifuged at 3000 RPM for 5 min and the mass of clay used ( 3.3 i 0.1 x 10'4 kg) and solvent entrained ( 3.6 i 0.1 x 10'3 kg) were determined gravimetrically on 6 randomly selected tubes by drying at 105 C for 48 hr. In the nonradioactive study each of the 12 treatments were duplicated and received 21.4 mL 0.01 M CaC12 containing from 13.2 uM to 12.8 mM Zn2+ as ZnC12 + 1 mL supporting electrolyte. Three replications 2+ were run in the tracer study, and 16 initial Zn concentrations ranged 19 from 9.3 nM to 12.8 mM. In addition, each tube received 1 mL lmCi 65Zn L'l' The initial Zn2+ concentrations were calculated on a 26 mL basis. All treatments were dispersed on a Vortex mixer and equilibrated for 1 wk at 27 i 1 C on a wrist action shaker. The samples were separated by centrifugation as previously described, and 4ml was removed for analysis. For the analytical study Zn2+ in the equilibrium solution was determined by standard AA methods using a Varian AA6 with an air- acetylene flame. Radiometric data were obtained. Zn2+ > Ca2+ > Mgz+ > K+. Garcia-Miragaya and Page (1977) reported that Cd2+ sorption by different homoionic montmorillonites followed the order Na > K > Ca > A1. Shukla et a1. (1980) found that Zn adsorption. by ‘whole soils which. were saturated with different cations followed the order H(A1) < Ca < Mg < K ,Na. Adsorption of alkali metal cation on smectite followed the order Li < Na < K < Rb < Cs. The free energy of adsorption was related to ion size parameters and polarizability factors (Gast, 1969). Although the cation composition of soils varies widely, there can be little doubt: the initial exchange ion and. associated adsorption 8O 81 energy will affect the position of equilibria in heavy metal adsorption- desorption studies and in natural systems. Since different cations affect the basal spacings of smectite minerals, these effects may be steric as well as electrostatic (MacEwan and Wilson, 1980). The objective in the present study was to determine the effect of initial exchange ion on adsorption and desorption of Zn by a smectite clay. Homoionic alkali metal smectites were prepared in an attempt to systematically vary the hydrated ion size and adsorption energy of the initial exchange ion. In this fashion, the periodic nature of homoionic 2+ clays and their reactivity toward Zn could be investigated. Materials and Methods Preparation of adsorbents Sodium and K saturated smectites (montmorillonite # 25, Upton, Wyoming) were prepared by repeated saturation with the appropriate 1 M MCl salt at pH 5.0. Rubidium and Cs smectites were prepared by repeated saturation with a 0.1 M_MC1 salt. In all cases, smectites were prepared from the < 2 um fraction. which. was recycled from previous Ca-Zn adsorption. experiments. The saturated. alkali metal smectites were dialyzed against distilled H20 until a negative test for C1" was obtained using Ag(NO3)2. The resulting clay was diluted with the appropriate 0.03 MCl salt solution and stored until ready for use. 82 Adsorption isotherms Batch adsorption experiments were conducted in duplicate for each of the alkali metal saturated clays by pipetting 3 mL of the rapidly stirred clay suspensions into 5 mL RIA polypropylene test tubes. Each tube received 0.2 mL of an appropriate Zn solution (Table 1) and 0.2 mL of a 652n solution containing lmCi 65Zn L'l. The mass of each clay used was determined gravimetrically. The tubes were then sealed, dispersed on a vortex mixer and equilibrated on a wrist action shaker at 26 i 2 C for 2 wks. At the end of this period the samples were separated by centrifugation at 3000 rpm for 5 min, and a 1 mL aliquot was removed for radioassqy. Eight incremental desorption cycles were accomplished by replacing the 1 mL sample with the appropriate 0.03 M MCl solution and re-equilibrating for 24 hrs. Radiometric data were obtained on a Packard Tri—Carb spectrometer equipped with a gamma counter by comparing sample radioactivity with that obtained for known 65Zn standards of the smae geometry. Equilibrium Zn concentrations were calculated from principles of isotopic dilution and mass balance. the amount of Zn adsorbed was calculated from the initial concentration. The equilibrium pH was determined on the 1 mL aliquot using a combination electrode. The rate of Zn adsorption at the lowest Zn2+ level was investigated by equilibrating duplicate samples for periods of 1, 2, a, 8 and 16 hr. Adsorbed Zn was determined as described above. The effect of longer equilibration times on adsorption was investigated at the 3 2+ lowest Zn levels by equilibrating duplicate samples for 24, 72, 168 and 672 hrs. 83 Table 1. Initial Zn2+ concentrations and equilibrium pH's for the adsorption of Zn by homoionic smectites Initial Zn2+ Equilibrium le Na K Rb Cs uM L’1 .0097 5.49 5.66 5.73 5.45 .147 5.48 5.69 5.74 5.43 1.45 5.50 5.62 5.73 5.53 146 5.49 5.65 5.72 5.40 1450 l 5.13 5.21 5.36 5.18 i Means of duplicate samples Results and Discussion In this study an attempt was made to maintain the pH at a constant value so that adsorption-desorption effects could be attributed to the exchangeable ion rather than to differences in the equilibrium pH. This objective was attained as the pH in subsequent desorption cycles did not change by mnore than i 0.1 unit, however the equilibrium pH's at the highest Zn level were lower, possibly due to hydrolysis in the bulk solution. Adsorption data for the different homoionic clays are presented as Freundlich isotherms. For Na smectite, adsorption obeyed the Freundlich equation over 6 orders of magnitude and appears to be completely reversible over the entire range studied (Figure 1). Conformity to the Freundlich equation and. complete reversibility ‘was observed for Zn 84 .3386 mzsnfiuofiafinomannafiflfiomfl A «36$ 33 :5 3 M N —. O _.l NI .MI .vl it p......._._.nl 1 YIN-I- .IPI. no IF 1m 5:938 4 . 5:90.30 0 b — _ p P P — . — h r b p p n (SD/W0) peqmspo uz 60' 85 adsorption on K smectite as well (Figure 2). Adsorption of Zn by Rb smectite obeyed the Freundlich equation to a degree, but more importantly, desorption hysteresis was observed at log (X/m) ==() and below (Figure 3). Am: higher surface occupancies of' Zn, adsorption appeared to be completely reversible as is expected in normal exchange reactions. For Cs smectite, desorption hysteresis was extreme throughout the entire adsorption range (Figure 4). Furthermore, the adsorption data did not conform to the Freundlich equation. For' Na, the desorption. points are tightly grouped. around the adsorption points, indicating that Na does not displace adsorbed Zn very effectively. As we proceed along the series Na, K, Rb, and Cs, there is a trend for the desorption points to spread out over a larger concentration range. This is an indication of the increasing ability of these cations to displace adsorbed Zn, and is in line with the free energy of adsorption of these cation on Wyoming montmorillonite as measured ‘by Cast (1969). These observations are also in general agreement with data presented by Nelson and Melsted (1955), Garcia- Miragaya. and Page (1977) and Shukla et al. (1980); (i.e that the smaller, more electropositive ion is adsorbed to a greater extent). That desorption hysteresis of Zn occurs on the smectites saturated with the smaller, more competative Rb and Cs ions seems to be a contradiction to the commonly observed lyotropic series. However, these observations can be explained from a kinetic viewpoint. At the lowest Zn2+ level, adsorption as a function of time by the different clays follows the order Na > K > Rb > Cs (Figure 5). This order was observed at higher initial Zn2+ levels as well. Zn adsorption 86 .Sflumfixfiéuofinfioflcgfig .mflsma 323 TE 00. m N P o T. NI m... T. - n — h P . — p — p b P L » nil INI I_.II no IF I. IN 5:333 4 . 339030 0 — p b n — - — p — L b bi _ . n (o/wn) pquospo uz 60' 87 .3386 fififimofiflnuomoogfiag ASE TE 3 .n 650E 5:833 4 coautomno o b t _ . _ . _ . _ F _ r (g/wn) pquospo uz bog 88 .3385 8356839032483 :\23 E as n N F 0 Pl NI _ h _ .4 993m 5:983 4 3:90qu o b r _ b _ L _ . — . b . _ . (g/wn) pquospo uz 601 89 .gumoagmmmmoug 33.584334482432000 833804 .0 303m 0?: 0E: 8:03:53 mm 4W mm mm mm mm Lm or m o . 0::00Em 02 o .. 0000 o r 0::00E0 v. a 3:025 gm 4 [3:095 00 o . m - . Q r _ _ F _ _ h _ _ .VNOO.O (Q/wn) pquospo UZ 90 by Na (Figure 6) and K smectite (Figure 7) was > 99% complete in less than. 24 hrs. This rapid equilibrium. is characteristic of an ion exchange reaction (Helfferich, 1962). Adsorption equilibria of Zn for Rb (Figure 8) and Cs smectites (Figure 9) was not attained even after 672 hrs. The different adsorption rates of Zn can be attributed to the hydrated nature of the individual alkali metal ions and their effect on the basal spacings of smectites (Wyoming montmorillonite) in H20 (Table 2). Sodium and K spacings were large and indeterminate, and while data were not available for Rb, the value is expected to be slightly higher than that reported for Cs. Table 2. Effect of interlayer cation on basal d(001) spacings of Wyoming montmorillonite immersed in water. Interlayer cation d(001) spacing1 10‘3m Na >> K >> Rb -- Cs 12.5 Mg 19.5 Ca 19.0 Sr 19.0 Ba 18.8 l MacEwan and Wilson (1980) p. 203 91 .3363335833085 23500888383008.4040 0?: 0E: Estezzoo R 4m 5. m: m 0 G b — 1 F} 0 HT 0 <2: mm o .. <2: 5:0 <2: 9... 4 T41 4 4 _ _ _ _ _ .0 2502 00.0 Iood IN—d [mfo I4N.O 10nd 0nd (g/wn) pquospo uz 92 3335 33 8058 338.5 00350083883830008000 03:30 mfcv 0E: Ezt£=aoo hm 4m 5 m__. m... N: m m m o _ _ ml .9 mi 00.0 U W! m .. -80 <2: 50 o . - <2: 5: .0 -Eo <2: 04.. 4 I -20 .. [a 4 1&0 4 .. -000 _ _ _ _ _ _ _ . $0.0 (Q/wn) peqmspo uz 93 .maosmacugugooufiuoofluoflfl 00344008380838330008000 0?: 0E: €39.53 R 4m 8 m: 9 m. 0 0 0 _ _ $— _ 6_ _ ml Lml u 1 <2: :0 o I <2: 9: 0 <2: 04.: 4 0 0.93m rood 19.0 1&6 I4N.o Iond 0nd (Q/wn) pquospo uz 94 83400083883830.0880 03:3,: 0?: 0E: 2:25:53 :m 4: E m: m: N. 0 0 m. o O . _ Ia— _ a — _QIIIIII.—IIQI _ 00.0 ml [ml u U n - .000 - -S0 <2: :0 o . -4. <2: :10 I: 0 <2: 0.: 4 - L400 - 000 r _ _ _ _ r _ L 0nd (g/wn) DGQJOSPD uz 95 Interestingly, the values reported for alkaline earth ions are relatively constant. The value for Zn is expected to be very similar since its crystal radii and aqueous chemistry is similar to that for Mg (Cotton and Wilkinson, 1980). This would account for the fact that equilibrium.‘was attained. quite rapidly for Ca-Zn exchange (Zielke, Chapter IV). Conversely, desorption hysteresis was shown to be concentration and pH dependent in these earlier studies, and was attributed to strong bonding, possibly at broken edge OH's. In the present study, similar hysteresis would be expected on Na (and K) smectites since these ions could not compete for specific adsorption sites in other studies (Harsh and Doner, 1984). One possible explaination for this behavior was that Zn remained fixed at broken edges when Ca-Zn clays were recycled for use in this experiment. Regardless, migration of Zn into interlamelar exchange sites on Rb and Cs smectites appear to 'be sterically inhibited. by the limited swelling as indicated by the nonattainment of equilibria. Similarly, once adsorbed at interlayer positions, migration into the bulk solution would be equally inhibited. As a result, the occurrence of hysteresis on Rb and Cs smectites appears to be a consequence of non equilibrium rather than strong irreversible bonding at specific adsorption sites. While this was a testable hypothesis, similar desorption kinetic experiments at the lowest Zn level failed to confirm this under these experimental conditions because the rate of «J urption proceeded faster than the rate of desorption even after 8 wks. In summary, Zn adsorption by homoionic alkali metal smectites followed the order Na > K > Rb > Cs. Data obtained in kinetic 96 experiments revealed that both adsorption and desorption of Zn was sterically inhibited by the limited swelling of Rb and Cs smectites. Desorption. hysteresis was therefore attributed. to ‘nonattainment of equilibria rather than to specific bonding as was reported in earlier Ca-Zn adsorption-desorption experiments. Failure to consider the initial exchanger composition. and its effect on the attainment of equilibrium, and the ability of some homoionic smectites (eg Cs) to produce nonlinear isotherms ‘may result in erroneous conclusions as regards desorption hysteresis or adsorption at sites with different bonding energies. 97 Literature Cited Cotton, F.A. and G. Wilkinsen. 1980. Advanced inorganic chemistry. (4th ed). Interscience Publishers. John Wiley & Sons. New York Egozy, Y. 1980. Adsorption of cadimum and cobalt on montmorillonite as a function of solution composition. Clays and Clay Miner. 28:311-318. Elrashidi, M.A., and G.A. O'Connor. 1982. Influence of solution composition on sorption of zinc by soils. Soil Sci. Soc. Am. J. 46:1153-1158. Garcia-Miragaya, J . , and A. L. Page . 1977 . Influence of exchangeable cation on the adsorption of trace amounts of cadmium by montmorillonite. Soil Sci. Soc. Am. J. 41:718-721. Gast, R.G. 1969. Standard free energies of exchange for alkali metal cations on Wyoming bentonite. Soil Sci.l Soc. Am. Proc. 33:37- 41. Helfferich, F. 1962. Ion exchange. McGraw-Hill. New York. Harsh, J.G. and H.E. Doner. 1984. Specific adsorption of copper on an hydroxy-aluminum-montmorillonite complex. Soil Sci. Soc. Am. J. 48:1034-1039. MacEwan, D.M.C. and M.J. Wilson. 1980. Interlayer and intercalation complexes of clay minerals. p. 197-248. In G.W. Brindley and G. Brown (ed.). Crystal structures of clay minerals and other X-ray identification. Mineralogical Society Monograph No. 5. Miner. Soc., London. Nelson, J.L., and S.W. Melsted. 1955. The chemistry of added zinc to soils and clays. Soil Sci. Soc. Am. J. 19:323-326. Shukla, V.C., S.B. Mittal, and R.K. Gupta. 1980. Zinc adsorption in some soils as affected by exchangeable cations. Soil Sci. 129:366-370. CHAPTER VI SUMMARY AND CONCLUSIONS A batch equilibrium technique based on principles of isotopic dilution and mass balance was developed to study adsorption and desorption of trace amounts of Zn by Ca smectite (montmorillonite # 25, Upton, Wyoming) The technique was tested by comparison with results obtained using conventional atomic adsorption for Zn analysis. Based on these comparisons it was concluded that 65Zn in solution was in kinetic equilibrium with 65Zn adsorbed on the surface of Ca smectite, and hence the technique was applicable to the study of desorption hysteresis in clay minerals. Desorption hysteresis occurred when <1% of the CEC was occupied by Zn but at higher occupancies Zn adsorption was completely reversible. Three isotherms (Linear, Freundlich and Langmuir) were tested to determine their utility in characterizing desorption. Data were obtained for Zn adsorption-desorption with Ca smectite, Charity clay (aeric haplaquept), a strong acid cation exchange resin (Ca Dowex 50W- X8), and a weak acid cation exchange-chelating resin (Ca Chelex). None of the equations could quantitatively describe Zn adsorption or desorption by and from Ca smectite or Charity clay since nonlinear isotherms were observed, but the Linear and Freundlich equations were better than the Langmuir. The Freundlich equation was a convenient 98 99 2+ method to view adsorption over several decades of Zn activity, even though log-log plots tend to minimize data variability. Adsorption. of Zn. by Ca. Dowex. was completely reversible. and quantitatively described by the Freundlich equation. A two surface Freundlich isotherm was obtained for Chelex, and hysteresis was observed in the high affinity region. The "two surface" behavior was attributed to reduced selectivity of the resin for Zn in acid solutions. Hysteresis was attributed to strong bonding by the chelating resin. Comparison of the behavior of the model ion exchangers with the smectite and Charity clay Freundlich isotherms allowed insight to the mechanisms operating in adsorption, desorption and hysteresis. For smectite the reaction was completely reversible at high surface occupancies, but hysteresis was observed. at low surface occupancy. Reversibility was attributed to an ion exchange mechanism. For Charity clay, extreme hysteresis occurred at all Zn levels. Hysteresis on both materials resembled the behavior of the Chelex resin in the high affinity region. For smectite, hysteresis was attributed to strong bonding at broken edge OH's, while Zn retention on the Charity clay could be due to strong bonding on broken edges of clay minerals, Al and Fe oxides, free carbonates or organic matter. There was a strong interaction between pH and hysteresis in Zn adsorption-desorption on Ca smectite. At pH 3.5, a linear Freundlich isotherm was observed and adsorption was reduced probably because of competition with Al3+ that resulted from acidic dissolution of the clay mineral. However, the reaction was completely reversible at all Zn levels. In general, as the pH increased, adsorption and hysteresis 100 increased. At the highest initial Zn2+ level, adsorption was relatively constant at all pH's and was completely reversible. Fhrthermore, the ability' of smectite clay to mask 'precipitation. was observed since adsorption of large amounts of Zn was reversible even in alkaline solution. Thus hysteresis was shown to be pH dependent while adsorption in the higher end of the isotherm could be attributed to ion exchange on permanently charged sites. In combination these factors contributed to nonlinear Freundlich isotherms. The pH dependency of hysteresis must be attributed to strong bonding at broken edge OH sites since these sites are considered to be the only source of pH dependent charge in the system. The effect of the initial exchange ion on Zn adsorption was also investigated using homoionic Na, K, Rb and Cs smectites. Hysteresis was not observed on Na an K smectites as was expected from the occurrence of hysteresis in previous Ca-Zn studies. This anomaly was attributed to incomplete removal of Zn from hysteretic sites because the smectites were 'prepared from clays recycled from previous Ca-Zn experiments. Regardless, the extent of Zn adsorption followed the lyotropic series Na > K > Rb > Cs. The ability of these ions to compete with Zn2+ for permanently charged exchange sites also followed this order and was observed as a spreading out of the desorption points along the lyotropic series. Extreme hysteresis was observed for Rb and Cs smectites, and while this appeared to be contradictory to their highly competitive nature, kinetic studies revealed hysteresis was an artifact of nonequilibrium. Adsorption and desorption was sterically inhibited at interlayer 101 exchange sites because Rb and Cs ions have a demonstrated ability to reduce swelling in smectite minerals thus reducing the rate of diffusion of Zn into and out of interlayer positions. APPENDIX APPENDIX NUCALC 10 REM NUCALC BAS CALCULATES PPM ZN IN SOLUTION FROM RADIACTIVE CPM AND PRINCIPLES OF ISOTOPIC DILUTION VIA A MASS BALANCE APPROACH AND PRODUCES A HARD COPY OF THE MODEL PARAMETERS AND INPUT AND OUTPUT FILES 20 REM NUCALC.BAS ALSO CALCULATES VALUES FOR CONSTRUCTION OF GENERAL,FREUNDLICH AND LANGMUIR ISOTHERMS 30 REM NUCALC.BAS WAS WRITEN BY' RICHARD CLEVELAND ZIELKE GRADUATE ASSISTANT IN SOIL CHEMISTRY MICHIGAN STATE UNIVERSITY CROP AND SOIL SCIENCES DEPARTMENT JULY 23 1986 40 LPRINT CHR$(27) "N"CHR$(6) 50 LPRINT CHR$(27) "M"; 60 LPRINT CHR$(27) "P" CHR$(15); 7o WIDTH "LPTl:",150 80 LPRINT CHR$(27) "1"CHR$(10); 9O LPRINT CHR$(27) "Q"CHR$(140); 100 DEFDBL Z,G,F,L 110 OPTION BASE 1 120 INPUT "WHAT IS THE NUMBER OF TREATMENTS?",N 130 INPUT "WHAT IS THE NUMBER OF ADSORPTION/DESORPTION CYCLES?",D 140 DIM CSTD(N,D),CSAM(N,D),ZI(N,D),ZF(N,D),ZREM(N,D) 150 DIM ZADS(N,D),ZINI(N,D),ZSOL(N,D),CORR(N,D) 160 DIM GX(N,D),GY(N,D),FX(N,D),FY(N,D),LX(N,D),LY(N,D) 170 REM CSTD-STANDARD CPM CORRECTED FOR COUNTING GEOMETRY CSAM-SAMPLE CPM ZI-INITIAL PPM ZN ZF-FINAL PPM ZN ZREM-MG ZN REMAINING IN SYSTEM AFTER EQUIL ZADS-MG ZN ADSORBED ZINI-MG ZN IN SYSTEM BEFORE EQUIL ZSOL-MG ZN IN SOLUTION AT EQUIL 180 REM CORR-SAMPLE COUNT ADJUSTED FOR COUNTING GEOMETRY 19o REM GX-ZN ACTIVITY (uM/L) GY-X/M (uM/G) FX-LOG [ZN] FY-LOG X/M LX=ZN ACTIVITY LY-C/X/M (uM/L/uM/G) 200 INPUT "WHAT IS MESSAGE?",A$ 210 INPUT "WHAT IS THE SAMPLE VOLUME IN MLS?",VS 22o INPUT "WHAT IS THE ANALYTE VOLUME IN MLS?",VA 230 INPUT "WHAT IS IONIC STRENGTH IN M/L?",U 240 INPUT "WHAT IS THE SAMPLE WEIGHT IN GRAMS? ",w 250 K—10‘(-.5085*4*U‘.5/(1+1.969*U‘.5)) :REM EXTENDED DE BYE HUCKEL EQN TO CALCULATE ZN ACTIVITY COEFFICIENT 260 INPUT "WHAT IS THE AVERAGE BACKGROUND COUNT?",CBKGD 27o INPUT "WHAT Is THE AVERAGE STANDARD COUNT?",CSTD 280 STD-CSTD-CBKGD 2 90 LPRI NT n *************************PROGRAM I S NUCALC.BAS************************* ":LPRINT "IDENTIFYING MESSAGE: ”;:LPRINT AS 300 LPRINT " *****INPUT PARAMETERS FOR SPECIFIED FILENAME***** " 102 103 310 LPRINT "SAMPLE VOLUME IN MLS-";VS 320 LPRINT "ANALYTE VOLUME IN MLS-";VA 330 LPRINT "IONIC STRENGTH-":U 340 LPRINT "ZN ACTIVITY COEFFICIENT- ";K 350 LPRINT "SAMPLE WEIGHT IN GRAMS-";W 360 LPRINT "BACKGROUND COUNT IN CPM-";CBKGD 370 LPRINT "UNCORRECTED STANDARD COUNT IN CPM-";CSTD 380 CLS:FILES "A:*.DAT" 390 INPUT "WHAT IS FILENAME CONTAINING CPM DATA (OMIT EXTENSION)?",F$ 400 S$-F$+".DAT" 410 LPRINT ”INPUT FILENAME-";S$:LPRINT 420 OPEN "I",#1,S$ 430 FOR I-l TO NzFOR J-l TO D:INPUT#1,CSAM(I,J):NEXTzNEXT 440 LPRINT S$;:LPRINT "(UNCORRECTED FOR BACKGROUND CPM)" 450 FOR I-l TO NzFOR. J-l TO DzLPRINT CSAM(I,J);" ";:NEXT:LPRINT," ":NEXTzLPRINT 460 CLOSE #1 470 T$-F$+".COR" 480 FOR I-l TO N:FOR J-l TO D:CSAM(I,J)-CSAM(I,J)-CBKGD:NEXT:NEXT 490 OPEN "O",#2,T$ 500 FOR I-l TO NzFOR J-l TO D:PRINT#2,CSAM(I,J);" ";:NEXT:PRINT#2," ":NEXT 510 LPRINT T$;:LPRINT "(CORRECTED FOR BACKGROUND CPM)" 520 FOR I-l TO NzFOR J-l TO DzLPRINT CSAM(I,J);" ":tNEXTzLPRINT," ":NEXT 530 LPRINT: 540 CLOSE #2 550 OPEN T$ FOR INPUT AS #3 560 FOR I-l TO N:FOR J-l TO D:INPUT#3,CSAM(I,J):NEXTzNEXT 570 LPRINT "INITIAL PPM ZN FINAL PPM ZN INITIAL MG ZN MG ZN REMAININ MG ZN ADSORBED" 530 '*************************START DATA CALCULATIONS ************************** 590 FOR I-1 TO N 600 FOR J-l TO D 610 IF J-l THEN CSTD(I,J)-STD*VS/VA 620 IF J>1 THEN CSTD(I,J)-CSTD(I,J-1)-CSAM(I,J-1) 630 CORR(I,J)-CSAM(I,J)*VS/VA 640 IF J-l THEN INPUT "WHAT IS THE INITIAL PPM ZN?",ZI(I,J) 650 IF J>1 THEN ZI(I,J)-ZF(I,J-1)*(VS-VA)/VS 660 IF J-l THEN ZINI(I,J)-ZI(I,J)*VS/1000 670 IF J>l THEN ZINI(I,J)-ZREM(I,J-l) 680 ZSOL(I,J)-CORR(I,J)*ZINI(I,J)/CSTD(I,J) 690 ZF(I,J)-ZSOL(I,J)*lOOO/VS 700 ZREM(I,J)-ZINI(I,J)-ZF(I,J)*VA/1000 710 ZADS(I,J)-ZINI(I,J)-ZSOL(I,J) 720 GX(I,J)—ZF(I,J)*K/.06537 730 IF J-l THEN GY(I,J)-(ZI(I,J)-ZF(I,J))*VS*1000/(65370!*W) 740 IF J>l THEN GY(I,J)-GY(I,J-1)+(ZI(I,J)-ZF(I,J))*VS*1000/(65370!*W) 750 FY(I,J)-LOG(GY(I,J))/2.30256 760 FX(I,J)-LOG(ZF(I,J)/.06537)/2.30256 77o LX(I,J)-GX(I,J) 780 LY(I,J)-LX(I,J)/GY(I,J) 104 790 LPRINT USING" M.M ###.######### ###.######### M.M M.M" ;ZI(I,J) ,ZF(I,J) ,ZINI(I,J) ,ZREM(I,J) ,ZADS(I,J) 800 NEXT J 810 LPRINT 820 NEXT I 830 '**************************END DATA CALCULATIONS*************************** 840 CLOSE #3 850 KILL T$ 860 U$-F$+".OUT" 870 OPEN "O",#3,U$ 880 FOR. I-l TO 1N:FOR. J-l TO D:PRINT#3,USING "###.####### ";ZI(I,J);ZF(I,J);:NEXT:PRINT#3,:NEXT 890 ‘LPRINT ‘U$;:LPRINT "(ADJACENT COLUMNS ARE INITIAL .AND FINAL ZN REMAINING IN SOLUTION AT EQUILIBRIUM. THE # OF COLUMNS/2 -# OF DESORPTION CYCLES. ADJACENT ROWS ARE REPS.)” 900 Ks-"INI PPM ZN":L$-"FIN PPM ZN" 910 LPRINT USING "\ \";K$,L$ 920 FOR I-l TO N:FOR J-l TO DzLPRINT USING "###.####### ';ZI(I,J);ZF(I,J);:NEXTzLPRINTzNEXT 930 CLOSE #3 940 GOSUB 1490 950 REM SUBROUTINE TO PRINT AND FILE GENERAL,FREUNDLICH AND LANGMUIR ISOTHERMS 960 H$-"G"+F$+".SDF" 970 I$-"F"+F$+".SDF" 980 J$-"L"+F$+".SDF" 990 OPEN "O",#1,H$ 1000 OPEN "O",#2,I$ 1010 OPEN "O",#3,J$ 1020 FOR I-I TO N:FOR. J-l TO D:PRINT#1,USING "####.####### ";GX(I,J);GY(I,J);:NEXT:PRINT#1,:NEXT 1030 LPRINTzLPRINT H$;:LPRINT " GENERAL ADSORPTION ISOTHERM DATA" 1040 GOSUB 1250 1050 M$-"X(uM ZN/L)":N$-"Y(uM ZN/G)" 1060 LPRINT USING "\ \";M$,N$ 1070 FOR I-l TO NzFOR. J-l TO DzLPRINT USING "####.####### ";GX(I,J);GY(I,J);:NEXT:LPRINT:NEXT:LPRINT 1080 FOR I-l TO NzFOR J-l TO D:PRINT#2,USING "##.##### ";FX(I,J);FY(I,J);:NEXT:PRINT#2,:NEXT 1090 LPRINT I$;:LPRINT " FREUNDLICH ADSORPTION ISOTHERM DATA" 1100 GOSUB 1330 1110 OS-"x LOG[ZN]":P$-"Y LOG X/M" 1120 LPRINT USING "\ \";O$,P$ 1130 FOR I-1 TO N:FOR J-I T0 D LPRINT USING "##.##### ”;FX(I,J);FY(I,J);:NEXT:LPRINT:NEXT:LPRINT 1140 FOR I-1 TO NzFOR. J-l TO D:PRINT#3,USING "####.####### ";LX(I,J);LY(I,J);:NEXT:PRINT#3,:NEXT 1150 LPRINT J$; :LPRINT " LANGMUIR ADSORPTION ISOTHERM DATA" 1160 GOSUB 1410 1170 Q$-"X(uM ZN/L)":R$-"Y(uM ZN/L/uM/G)" 1180 LPRINT USING "\ \";Q$,R$ 1190 105 FOR I-l TO NzFOR J-l TO D:LPRINT USING "####.####### ";LX(I,J);LY(I,J);:NEXTzLPRINTzNEXTzLPRINT 1200 1210 1220 1230 1240 1250 1260 1270 CLOSE #1:CLOSE #2: CLOSE #3 LPRINT CHR$(27) "@" WIDTH "LPTl:",8O END REM *****SUBROUTINE TO PRINT REGRESSION EQUATIONS LPRINT "FOR THE GENERAL ISOTHERM (ZN) VS X/M " LPRINT "B-";B9;TAB(25),"A-";A3;TAB(4S);"R92-";R9 YS-B9*GX(1,1)+A3 1280 Y6-B9*GX(N,1)+A3 1290 LPRINT "X5-";GX(1,1);TAB(40),"Y5-";Y5 1300 LPRINT "X6-";GX(N,1);TAB(40),”Y6-";Y6 1310 LPRINT "Y-";B9;"X+";A3 1320 RETURN 1330 LPRINT ”FOR THE FREUNDLICH EQUATION LOG [ZN] VS LOG X/M" 1340 LPRINT "B-";B3;TAB(25),"A-";A1;TAB(45);"R“2-";R3 1350 Y1-B3*FX(1,1)+A1 1360 Y2-B3*FX(N,1)+A1 1370 LPRINT "X1-";FX(1,1);TAB(40),"Yl-";Y1 1380 LPRINT ”X2-";FX(N,1);TAB(40),"Y2-";Y2 1390 LPRINT "Y—";B3;"X+";A1 1400 RETURN 1410 LPRINT "FOR THE LANGMUIR ISOTHERM (ZN) VS C/X/M" 1420 LPRINT "B-";B12;TAB(25),"A-";A4;TAB(45),"R‘2-";R12 1430 Y7-B12*LX(1,1)+A4 1440 Y8-812*LX(N,1)+A4 1450 LPRINT "X7-";LX(1,1);TAB(40),"Y7-";Y7 1460 LPRINT ”X8—";LX(N,1);TAB(40),"Y8-";Y8 1470 LPRINT "Y-";B12;"X+";A4;TAB(40),"APPARENT ADSORPTION MAXIMUM= ";B14;” (uM Zn/G)" 1480 RETURN 1490 REM SUBROUTINE TO CALCULATE REGRESSION EQUATIONS 1500 LOXY-O 1510 PQXY-O 1520 PXXY-O 1530 XSUMLO-O 1540 YSUMLO-O 1550 XSUMPQ-O 1560 YSUMPQ-O 1570 XSUMPX-O 1580 YSUMPX-O 1590 XSSLO-O 1600 YSSLO-O 1610 XSSPQ-O 1620 YSSPQ-O 1630 XSSPX-O 1640 YSSPx-O 1650 FOR I-1 TO N 1660 REM *****GENERAL 1670 PXXY-PXXY+GX(I,1)*GY(I,1) 1680 1690 1700 XSUMPX-XSUMPX+GX(I,1) YSUMPX-YSUMPX+GY(I,1) XSSPx-XSSPX+GX(I,1)‘2 1710 1720 1730 1740 1750 1760 1770 1780 1790 1800 1810 1820 1830 1840 1850 1860 1870 1880 1890 1900 1910 1920 1930 1940 1950 1960 1970 1980 1990 2000 2010 2020 2030 2040 2050 2060 2070 2080 2090 2100 10 106 YSSPx-YSSPx+GY(I,1)‘2 REM *****FREUNDLICH LOXY-LOXY+FY(I,1)*FX(I,1) XSUMLO-XSUMLO+FX(I,1) YSUMLO-YSUMLO+FY(I,1) XSSLO-XSSLO+FX(I,1)‘2 YSSLO-YSSLO+FY(I,1)“2 REM *****LANGMUIR PQXY-PQXY+LX(I,1)*LY(I,1) XSUMPQ-XSUMPQ+LX(I,1) YSUMPQ-YSUMPQ+LY(I,1) XSSPQ-XSSPQ+LX(I,1)‘2 YSSPQ-YSSPQ+LY(I,1)“2 NEXT I REM *****TO CALCULATE GENERAL REGRESSION B7-PXXY-XSUMPX*YSUMPX/N BB-XSSPx-XSUMPx‘z/N B9-B7/BS A3-YSUMPX/N-B9*XSUMPX/N R7-B7“2 R8-B8*(YSSPx-YSUMPX“2/N) R9-R7/R8 REM *****T0 CALCULATE FREUNDLICH REGRESSION B1-LOXY-XSUMLO*YSUMLO/N 32-XSSLO-XSUMLO‘2/N B3-B1/B2 A1-YSUMLO/N-B3*XSUMLO/N R1-B1“2 R2-Bz*(YSSLO-YSUML0“2/N) R3-R1/R2 REM *****TO CALCULATE LANGMUIR REGRESSION B10-PQXY-XSUMPQ*YSUMPQ/N Bll-XSSPQ-XSUMPQAZ/N B12-BlO/B11 A4-YSUMPQ/N-B12*XSUMPQ/N R10-B10‘2 R11-BII*(YSSPQ-YSUMPq‘z/N) R12-R10/R11 B14-B11/B10 RETURN ZNCALC REM ZNCALC.BAS CALCULATES ‘VALUES MEASURED BY AA SPECTROSCOPY 20 REM ZNCALC.BAS WAS WRITEN BY' RICHARD CLEVELAND ZIELKE GRADUATE ASSISTANT IN SOIL CHEMISTRY MICHIGAN STATE UNIVERSITY CROP AND SOIL SCIENCES DEPARTMENT JULY 23 1986 30 LPRINT CHR$(27) "N"CHR$(6) 40 LPRINT CHR$(27) "M"; 50 LPRINT CHR$(27) "P" CHR$(15); CONSTRUCTION GENERAL,FREUNDLICH AND LANGMUIR ISOTHERMS USING ZN(I,J) AND ZF(I,J) AS 107 60 WIDTH 'LPTl:",150 7o LPRINT CHR$(27) "1"CHR$(10); 80 LPRINT CHR$(27) "Q"CHR$(140); 90 DEFDBL z,G,F,L 100 OPTION BASE 1 110 INPUT "WHAT Is THE NUMBER OF TREATMENTS?",N 120 INPUT "WHAT IS THE NUMBER OF ADSORPTION/DESORPTION CYCLES?",D 130 DIM ZI(N,D),ZF(N,D),ZREM(N,D) 140 DIM ZADS(N,D),ZINI(N,D),ZSOL(N,D),ZDES(N,D) 150 DIM GX(N,D),GY(N,D),FX(N,D),FY(N,D),LX(N,D),LY(N,D) 16o ZI-INITIAL PPM ZN ZF-FINAL PPM ZN ZREM-MG ZN REMAINING IN SYSTEM AFTER EQUIL ZADS-MG ZN ADSORBED ZINI-MG ZN IN SYSTEM BEFORE EQUIL ZSOLPMG ZN IN SOLUTION AT EQUIL ZDEs-MG ZN DESORBED 17o REM GX-ZN ACTIVITY (uM/L) GY-X/M (uM/G) FX-LOG [ZN] FY=LOG X/M LX=ZN ACTIVITY LY-C/X/M (uM/L/uM/G) 180 INPUT "WHAT IS MESSAGE?",A$ 19o INPUT "WHAT IS THE SAMPLE VOLUME IN MLS?",VS 200 INPUT "WHAT IS THE ANALYTE VOLUME IN MLS?",VA 210 INPUT "WHAT IS IONIC STRENGTH IN M/L?",U 220 INPUT "WHAT IS THE SAMPLE WEIGHT IN GRAMS? ",W 230 K:10“(-.5085*4*U“.5/(1+1.969*U“.5)) :REM EXTENDED DE BYE HUCKEL EQN TO CALCULATE ZN ACTIVITY COEFFICIENT 240 LPRINT n *************************PROGRAM IS ZNCALC.BAS************************* ":LPRINT "IDENTIFYING MESSAGE: ":zLPRINT A$ 250 LPRINT " *****INPUT PARAMETERS FOR SPECIFIED FILENAME*****” 260 LPRINT "SAMPLE VOLUME IN MLS-";VS 270 LPRINT "ANALYTE VOLUME IN MLS-";VA 280 LPRINT "IONIC STRENGTH-";U 290 LPRINT "ZN ACTIVITY COEFFICIENT- ":K 300 LPRINT "SAMPLE WEIGHT IN GRAMS-";W 310 INPUT "WHAT IS FILENAME CONTAINING ZF(I,J) DATA (OMIT EXTENSION)?",F$ 320 S$-F$+".DAT" 330 LPRINT "INPUT FILENAME-";S$:LPRINT 340 OPEN "I",#1,S$ 350 FOR I-1 TO NzFOR J-l TO D:INPUT#1,ZF(I,J):NEXTzNEXT 360 LPRINT "INITIAL PPM ZN FINAL PPM ZN INITIAL MG ZN MG ZN REMAININ MG ZN ADSORBED MG ZN DESORBED" 37o '*************************START DATA CALCULATI ONS ************************** 380 FOR I-l TO N 390 FOR J-l TO D 400 IF J-l THEN INPUT ”WHAT IS THE INITIAL PPM ZN?",ZI(I,J) 410 IF J>l THEN ZI(I,J)-ZF(I,J-1)*(VS-VA)/VS 420 IF J-l THEN ZINI(I,J)-ZI(I,J)*VS/1000 430 IF J>l THEN ZINI(I,J)-ZREM(I,J-1) 440 ZSOL(I,J)-ZF(I,J)*VS/1000 450 ZREM(I,J)-ZINI(I,J)-ZF(I,J)*VA/1000 460 ZADS(I,J)-ZINI(I,J)-ZSOL(I,J) 470 IF J-l THEN ZDES(I,J)-O ELSE ZDES(I,J)-ZADS(I,J-I)-ZADS(I,J) 480 GX(I,J)-ZF(I,J)*K/.06537 490 IF J-l THEN GY(I,J)-(ZI(I,J)-ZF(I,J))*VS*1000/(65370!*W) 108 500 IF J>1 THEN GY(I,J)-GY(I,J-1)+(ZI(I,J)-ZF(I,J))*VS*1000/(6S370!*W) 510 FY(I,J)-LOG(GY(I,J))/2 30256 520 FX(I,J)-LOG(ZF(I,J)/.06537)/2.30256 530 LX(I,J)-GX(I,J) 540 LY(I,J)-LX(I,J)/GY(I,J) 550 LPRINT USING" ###.#########= ###.######### :###.######### ###.######### ###.######### M.M";ZHLJ),ZF(I,J),ZINI(I,J),ZREM(I,J),ZADS(I,J),ZDES(I,J) 560 NEXT J 570 LPRINT 580 NEXT I 590 '**************************END DATA CALCULATIONS*************************** 600 CLOSE #1 610 U$-F$+".OUT" 620 OPEN "O",#3,U$ 630 FOR I-l TO NzFOR J-l TO D:PRINT#3,USING "##4##” ";ZI(I,J);ZF(I,J);:NEXT:PRINT#3,:NEXT 640 LPRINT U$;:LPRINT "(ADJACENT COLUMNS ARE INITIAL .AND FINAL ZN REMAINING IN SOLUTION AT EQUILIBRIUM. THE # OF COLUMNS/2 -# OF DESORPTION CYCLES. ADJACENT Rows ARE REPS.)" 650 K$-"INI PPM ZN':L$-"FIN PPM ZN" 660 LPRINT USING "\ \";K$,L$ 670 FOR I-I TO N:FOR J-l TO D:LPRINT USING "###.####### ";ZI(I,J);ZF(I,J);:NEXTzLPRINTzNEXT 680 CLOSE #3 690 GOSUB 1220 700 REM SUBROUTINE TO PRINT AND FILE GENERAL,FREUNDLICH AND LANGMUIR ISOTHERMS 710 H$-"G"+F$+".SDF" 720 I$-"F"+F$+".SDF" 730 J$-"L"+F$+".SDF" 740 OPEN "O",#1,H$ 750 OPEN "O",#2,I$ 760 OPEN "O",#3,J$ 770 FOR I—1 TO NzFOR J-l TO D:PRINT#1,USING "####.####### ":GX(I,J);GY(I,J);:NEXT:PRINT#1,:NEXT 780 LPRINT:LPRINT H$;:LPRINT " GENERAL ADSORPTION ISOTHERM DATA" 790 GOSUB 980 800 M$-"X(uM ZN/L)":N$-"Y(uM ZN/G)" 810 LPRINT USING "\ \";M$,N$ 820 FOR I-1 TO NzFOR J-l TO D:LPRINT USING "####.####### ";GX(I,J);GY(I,J);:NEXT:LPRINT:NEXT:LPRINT 830 FOR I-1 TO N:FOR J-l TO D:PRINT#2,USING "##.##### ";FX(I,J);FY(I,J);:NEXT:PRINT#2,:NEXT 840 LPRINT Is; LPRINT " FREUNDLICH ADSORPTION ISOTHERM DATA" 850 GOSUB 1060 860 O$-"X LOG[ZN]":P$-"Y LOG X/M" 870 LPRINT USING "\ \";O$,P$ 880 FOR I-1 TO NzFOR J-I TO D:LPRINT USING "##.##### ";FX(I,J);FY(I,J);:NEXT:LPRINT:NEXT:LPRINT 890 FOR I-l TO NzFOR J-1 TO D:PRINT#3,USING "####.####### ";LX(I,J);LY(I,J);:NEXT:PRINT#3,:NEXT 900 LPRINT J$;:LPRINT " LANGMUIR ADSORPTION ISOTHERM DATA" 109 910 GOSUB 1140 920 Q$-"X(uM ZN/L)":R$-"Y(uM ZN/L/uM/G)" 930 LPRINT USING "\ I-l 940 FOR \";Q$.R$ TO N:FOR J-l TO D:LPRINT USING "####.####### ":LX(I,J);LY(I,J);:NEXTzLPRINTzNEXTzLPRINT 950 CLOSE #1:CLOSE #2: CLOSE #3 960 END 970 REM *****SUBROUTINE TO PRINT REGRESSION EQUATIONS 980 LPRINT "FOR THE GENERAL ISOTHERM (ZN) VS X/M " 990 LPRINT "B-";B9;TAB(25),"A-";A3;TAB(45);"R“2-";R9 Y5-B9*GX(1,1)+A3 Y6-B9*GX(N,1)+A3 1000 1010 1020 1030 1040 1050 1060 1070 1080 1090 1100 1110 1120 1130 1140 1150 1160 1170 1180 1190 1200 LPRINT LPRINT LPRINT RETURN LPRINT LPRINT ”X5-";GX(1,1);TAB(40),"Y5-";Y5 "X6-”;GX(N,1);TAB(40),"Y6-";Y6 "Y-";B9;"X+";A3 "FOR THE FREUNDLICH EQUATION LOG [ZN] VS LOG X/M" "B-";B3;TAB(25),"A-";A1;TAB(45);"R“2-";R3 Y1-BB*FX(1,1)+A1 Y2-B3*FX(N,1)+A1 LPRINT LPRINT LPRINT RETURN LPRINT LPRINT "X1—";FX(1,1);TAB(40),”Y1—";Y1 "X2-";FX(N,1);TAB(40),"Y2-";Y2 ”Y5";B3;"X+";A1 "FOR THE LANGMUIR ISOTHERM (ZN) VS C/X/M" "B-";312;TAB(25),"A-";A4;TAB(45),"RA2-";R12 Y7-BlZ*LX(l,1)+A4 Y8-BlZ*LX(N,l)+A4 LPRINT "X7-";LX(1,1);TAB(40),"Y7-";Y7 LPRINT "X8-";LX(N,1);TAB(40),"Y8-";Y8 LPRINT "Y-" ;812 ; "X+" ;A4 ; TAB(40) , "APPARENT ADSORPTION MAXIMUM= ";B14;" (uM Zn/G)" 1210 1220 1230 1240 1250 1260 1270 1280 1290 1300 1310 RETURN REM SUBROUTINE TO CALCULATE REGRESSION EQUATIONS LOXY-O PQXY-O PXXY-O XSUMLO-O YSUMLO-o XSUMPQ-O YSUMPQ-O XSUMPX-o YSUMPX-O 1320 1330 1340 1350 1360 1370 1380 1390 1400 1410 1420 XSSLO-O YSSLO-O XSSPQ-O YSSPQ-O XSSPX-o YSSPX-O FOR I-1 TO N REM *****GENERAL PXXY-PXXY+GX(I,1)*GY(I,1) XSUMPX-XSUMPX+GX(I,1) YSUMPX-YSUMPX+GY(I,1) 1430 1440 1450 1460 1470 1480 1490 1500 1510 1520 1530 1540 1550 1560 1570 1580 1590 1600 1610 1620 1630 1640 1650 1660 1670 1680 1690 1700 1710 1720 1730 1740 1750 1760 1770 1780 1790 1800 1810 1820 1830 110 XSSPX-XSSPX+GX(I,1)‘2 YSSPX-YSSPX+GY(I,1)‘2 REM *****FREUNDLICH LOXY-LOXY+FY(I,1)*FX(I,1) XSUMLO-XSUMLO+FX(I,1) YSUMLO-YSUMLO+FY(I,1) XSSLO-XSSLO+FX(I,1)“2 YSSLO—YSSLO+FY(I,1)“2 REM *****LANGMUIR PQXY-PQXY+LX(I,1)*LY(I,1) XSUMPQ-XSUMPQ+LX(I,1) YSUMPQ-YSUMPQ+LY(I,1) XSSPQ-XSSPQ+LX(I,1)“2 YSSPQ-YSSPQ+LY(I,1)‘2 NEXT I REM *****TO CALCULATE GENERAL REGRESSION B7-PXXY-XSUMPX*YSUMPX/N BB—XSSPX-XSUMPX‘z/N B9-B7/BB A3-YSUMPX/N-B9*XSUMPX/N R7-B7“2 R8-BB*(YSSPX-YSUMPX“2/N) R9-R7/R8 REM *****TO CALCULATE FREUNDLICH REGRESSION B1-LOXY-XSUMLO*YSUMLO/N Bz-XSSLO-XSUMLo‘z/N B3-B1/B2 A1-YSUMLO/N-B3*XSUMLO/N R1-BI‘2 R2-B2*(YSSLO-YSUML0“2/N) R3-R1/R2 REM *****TO CALCULATE LANGMUIR REGRESSION BlO-PQXY-XSUMPQ*YSUMPQ/N B11-XSSPQ-XSUMPQ‘2/N B12—B10/B11 A4-YSUMPQ/N-B12*XSUMPQ/N R10-BIO‘2 R11-B11*(YSSPQ-YSUMPQ“2/N) R12-R10/R11 B14-Bll/BlO RETURN LI ST OF REFERENCES LIST OF REFERENCES Bar-Yosef, B. 1979. pH-dependent zinc adsorption by soils. Soil Sci. Soc. Am. J. 43:1095-1099. Barney, G.S., 1984. Geochemical behavior of disposed radioactive waste. ACS Symposium Series No. 246:3-23. American Chemical Society. Washington D.C. Bingham, F.T., A.L. Page, and J.R. Sims. 1964. Retention of Cu and Zn by H-montmorillonite. Soil Sci. Soc. Am. J. 28:351-354. Borchardt, G.A., 1977. Montmorillonite and other smectite minerals. p. 293-330. J.B. Dixon, and S.B. Weed (ed.) E Minerals in soil environments. Soil Science Society of America. Madison, WI. Bourg, A.C.M., and R.H. Filby. 1976. Isotopic exchange of Zn65 with stable Zn adsorbed on reference clay minerals. Geochim. Cosmochim. Acta. 40:1573-1574. Bourg, A.C.M. , and R.H. Filby. 1974 Adsorption isotherms for the uptake of zinc by clay minerals in a freshwater medium. In Proceedings of the International Conference on Persistent Chem. Aquat. Ecosystems, Vol II. National Research Council of Canada. *Ottowa, Ont. Boyd, G.E. J. Shubert, and A.W Adamson. 1947. The exchange adsorption of ions from aqueous solutions by organic zeolites I. Ion exchange equilibria. J. Am. Chem. Soc. 69:2818-2829. Brummer, G., K.G. Tiller, V. Herms, and P.M. Clayton. 1983. Adsorption- desorption and/or precipitation-dissolution processes of zinc in soils. Geoderma. 31:337-354. Burkholder, H.C. 1976. Methods and data for predicting nuclide migration in geologic media. pp. 658—666. In Proc. Int. Symp. on Manage. of Wastes from LWR Pnel Cycle. Denver, Colo. 11-16 July 1976. CONF- 760701. Natl. Tech. Information Serv., U.S. Dept. of Commerce, Springfield, Va. Cavallaro, N., and M.B. McBride. 1984. Zinc and copper sorption and fixation by an acid soil clay: effect of selective dissolutions. Soil Sci. Soc. Am. J. 48:1050-1054. Corey, R.B. 1981. Adsorption vs precipitation. p.161-182. 13 M.A. Anderson and A.J. Rubin (ed). Adsorption of inorganics at solid- liquid interfaces. Ann Arbor Science. Ann Arbor, Mich. 111 112 Cotton, F.A. and G. Wilkinson. 1980. Advanced inorganic chemistry. (4th ed) Interscience Publishers. John Wiley & Sons. New York. DeMumbrum, L.E. , and M.L. Jackson. 1956. Infrared adsorption evidence on exchange reaction mechanism of copper and zinc with layer silicate clays and peats. Soil Sci. Soc. Am. J. 20:334-337. Egozy, Y. 1980. Adsorption of cadmium and cobalt on montmorillonite as a function of solution composition. Clays and Clay Miner. 28:311- 318. Elgabaly, M.M. 1950. Mechanism of zinc fixation by colloidal clays and related minerals. Soil Sci. 69:167-174. Elrashidi, M.A., and G.A. O'Connor. 1982. Influence of solution composition on sorption of zinc by soils. Soil Sci. Soc. Am. J. 46:1153-1158. Farrah, H., and W.F. Pickering. 1976. The sorption of zinc species by clay minerals. Aust. J. Chem. 29:1649-1656. Farrah, H., and W.F. Pickering. 1977. Influence of clay-solute interactions on aqueous heavy metal ion levels. Water, Air, Soil Pollut. 8:189-197. Farrah, H., and W.F. Pickering. 1978. Extraction of heavy metal ion sorbed on clays. Water, Air, Soil Pollut. 9:491-498. Freundlidh, H. 1926. Colloid and capillary chemistry. Translated from the third German edition by H. Stafford Hatfield. Methven & Co. Ltd. London. Frost, R.R., and R.A. Griffin. 1977. Effect of pH on adsorption of copper, zinc and cadmium from landfill leachate by clay minerals. J. Environ. Sci. Health, Part A. 12:139-156. Fritz, J.S., and G.H Schenk. 1979. Quantitative analytical chemistry. Allyn an Bacon, Boston. Gaines, G.L., and H.C. Thomas. 1953. Adsorption studies on clay minerals, II: A formulation of the thermodynamics of exchange adsorption. J. Chem. Phys. 21:714-718. Gamble, D.S., and M. Schnitzer. 1973. Trace metals and metal organic interactions in natural waters. Ann Arbor Science, Ann Arbor Garcia-Miragaya, J . , and A.L. Page . 1977 . Influence of exchangeable cation on the adsorption of trace amounts of cadmium by montmorillonite. Soil Sci. Soc. Am. J. 41:718-721. Cast, R.G. 1969. Standard free energies of exchange for alkali metal cations on Wyoming bentonite. Soil Sci. Soc. Am. Proc. 33:37-41. Helfferich, F. 1962. Ion exchange. McGraw-Hill. New York. 113 Garcia-Miragaya, J., and A.L. Page. 1976. Influence of ionic strength and inorganic complex formation on the sorption of trace amounts of Cd by montmorillonite. Soil Sci. Soc. Am. J. 40:658-663. Garcia-Miragaya, J . , and A.L. Page . 1977 . Influence of exchangeable cation on the adsorption of trace amounts of cadmium by montmorillonite. Soil Sci. Soc. Am. J. 41:718-721. Griffin, R.A., and A.K. Au. 1977. Lead adsorption by montmorillonite using a competitive Langmuir equation. Soil Sci. Soc. Am. J. 41:880-882. Grim, R.B. 1968. Clay mineralogy. 2nd ed. McGraw-Hill, New York. Harsh, J.B, and H.E. Doner. 1984. Specific adsorption of copper on an hydroxy-aluminum-montmorillonite complex. Soil Sci. Soc. Am. J. 48:1034-1039. Harsh, J.B. and H.E. Doner. 1984. Specific adsorption of copper on an hydroxy-aluminum-montmorillonite complex. Soil Sci. Soc. Am. J. 48:1034-1039. Harter, R.D. 1983. Effect of soil pH on adsorption of lead, copper, zinc and nickel. Soil Sci. Soc. Am. J. 47:47-51. Harter, R.D., and D.E. Baker. 1977. Application and misapplications of the Langmuir equation to soil adsorption phenomena. Soil Sci. Soc. Am. J. 41:1077-1080. Helfferich, F. 1962. Ion exchange, McGraw-Hill. New York. Hendrickson, L.L., and R.B. Corey. 1983. A chelating resin method for characterizing soluble metal complexes. Soil Sci. Soc. Am. J. 47:467-474. Hendrickson, L.L., M.A. Turner, and R.B. Corey. 1982. Use of Chelex-100 to maintain constant metal activity and its application to characterization of metal complexation. Anal. Chem. 54:1633-1637. Helfferich, F. 1962. Ion exchange, McGraw-Hill. New York. Hingston, F.J., A.M. Posner, and J.P Quirk. 1972. Anion adsorption by goethite and gibbsite. I: The role of the proton in determining adsorption envelopes. J. Soil Sci. 23:177-192. Hodgson, J.F., K.G. Tiller, and M. Fellows. 1964. The role of hydrolysis in the reaction of heavy metals with soil forming materials. Soil Sci. Soc. Am. J. 28:42-46. Jackson, M.L., 1958. Soil chemical analysis. Prentice-Hall. New Jersy. Jenny, H. 1961. Reflections on the soil acidity merry-go-round. Soil Sci. Soc. Am. J. 25:428-432. ‘ .- 1"! 114 Jurinak, J.J., and N. Bauer. 1956. Thermodynamics of zinc adsorption in calcite, dolomite and *manesite-type minerals. Soil Sci. Soc. Am. J. 20:466-471 Kalbasi, M., G.J. Racz and L.A. Bowen-Rudgers. 1978. Mechanism of zinc adsorption by iron and aluminum oxides. Soil Sci. 125:146-150. Keren, R. 1980. Effects of hydroxy aluminum precipitation on the exchange capacity reduction and aggregate size distribution of montmorillonite hydroxy-aluminum complexes. Soil Sci. Soc. Am J. 44:1209-1212. Kinniburgh, D.G., and M.L. Jackson. 1982. Concentration and pH dependence of calcium and zinc adsorption by iron hydrous oxide gel. Soil Sci. Soc. Am. J. 46:56-61. Kittrick, J.A. 1976. Controls of Zn2+ in the soil solution by *sphalerite. Soil Sci. Soc. Am. J. 40:314-317. Kuo, S., and A.S. Baker. 1980. Sorption of copper, zinc and cadmium by some acid soils. Soil Sci. Soc. Am. J. 44:969-974. Kuo, S., and D.S. Mikkelsen. 1979. Zinc adsorption by two alkaline soils. Soil Sci. 128:274-279. Kurdi, F., and Doner. 1983. Zinc and copper sorption and interaction in soils. Soil Sci. Soc. Am. J. 47:873-876. Lai, S., J.J. Jurinak, and R.J. Wagenet. 1978. Multicomponent cation adsorption during convective-dispersive flow through soils: experimental study. Soil Sci. Soc. Am. J. 42:240-243. Lamm, C.G., E.H. Hansen, and J.A. Jorgensen. 1963. Isotopic exchange in soil fertility studies. Soil Sci. 94:16-2Maes, A., and A. Cremers. 1975. Cation-exchange hysteresis in montmorillonite: a pH dependent effect. Soil Sci. 119:198-202. Langmuir, I. 1918. The adsorption of gases on plane surfaces of glass, mica and platinum. J. Am. Chem. Soc. 40:1361-1403. Lindsay, W.L. 1979. Chemical equilibria in soils. Wiley, New York. Lindsay, W.L. 1972. Inorganic phase equilibria of micronutrients in soils. p.4l-57. In J.J. Morvedt, P.M. Giordano, and W.L. Lindsay (ed.) Micronutrients in agriculture. Soil Science Society of America. Madison, WI. Maes, A., and A. Cremers. 1975. Cation-exchange hysteresis in montmorillonite: a pH dependent effect. Soil Sci. 119:198-202. Maquire, M., J. Slavek, I. Vimpany, F.R. Higginson, and W.F. Pickering. 1981. Influence of pH on copper and zinc uptake by soil clays. Aust. J. Soils Res. 19:217-229. 115 MacEwan, D.M.C. and M.J. Wilson. 1980. Interlayer and intercalation complexes of clay minerals. p. 197-248. In G.W. Brindley and G. Brown (ed.). Crystal structures of clay minerals and other X-ray identification. Mineralogical Society Monograph No. 5. Miner. Soc., London. Marshall, G.E. 1964. The physical chemistry and mineralogy of soils. Wiley, New York. Mattigod, S.V., A.S. Gibali, and..A.L. Page. 1979. Effect of ionic strength and ion pair formation on the adsorption of nickel by kaolinite. Mattigod, S.V., and G. Sposito. 1977. Estimated association constants for some complexes of trace metals with inorganic ligands. Soil Sci. Soc. Am. J. 41:1092-1097. McBride, M.B., and J.J. Blasiak. 1979. Zinc and copper solubility as a function of pH in an acid soil. Soil Sci. Soc. Am. J. 43:866-870. McBride, M.B., and M.M Mortland. 1974. Copper (11) interactions with montmorillonite: evidence from physical methods. Soil Sci. Soc. Am. J. 38:408-414. McLaren, R.G., and J.G. Williams, and R.S. Swift. 1983. Some observations on the desorption and distribution behavior of copper with soil components. J. Soil Sci. 134:325-331. Nelson, J.L., and S.W. Melsted. 1955. The chemistry of added zinc to soils and clays. Soil Sci. Soc. Am. J. 19:323-326. Nicholas, D.J.D, and A.R. Egan. 1975. Trace elements in soil-plant- animal systems. Academic Press. New York. Pionke, H.B. and RAB. Corey. 1967. Relations between acidic aluminum and soil pH, clay and organic matter. Soil Sci. Soc. Am. Proc. 31:749-752. Posner, A.M, and J.W Bowden. 1980. Adsorption isotherms: should they be split? J. Soil Sci. 31:1-10. Saeed, M., and R.F. Fox. 1977. Relations between suspension pH and zinc solubility in acid and calcareous soils. Soil Sci. 124:199-203. Shukla, V.C., S.B. Mittal, and R.K. Gupta. 1980. Zinc adsorption in some soils as affected by exchangeable cations. Soil Sci. 129:366-370. Shuman, L.M. 1980. Zinc in soils. p. In J. Nriagu (ed) Zinc in the environment. Part II. Wiley, New York. Shuman, L.M. 1977. Adsorption of Zn by Fe and Al hydrous oxides as influenced by aging and pH. Soil Sci. Soc. Am. J. 41:703-706. 116 Shuman, L.M. 1976. Zinc adsorption isotherms for soil clays with and without iron oxides removed. Soil Sci. Soc. Am. J. 40:349-352. Shuman, L.M. 1975. The effect of soil properties on zinc adsorption by soils. Soil Sci. Soc. Am. J. 39:454-458. Singhal, J.P., and D. Kumar. 1977. Thermodynamics of exchange of zinc with Mg-bentonite and Mg-illite. Geoderma. 17:253-258. Sposito, G. 1982 On the use of the Langmuir equation in the interpretation of "adsorption" phenomenon :II. The "two surface" Langmuir equation. Soil Sci. Soc. Am. J. 46:1147-1152. Tabikh, A.A., I. Barshad, and R. Overstreet. 1960. Cation exchange hysteresis in clay minerals. Soil Sci. 104:219-226. Tiller, K.G. 1967. Silicic acid and the reaction of zinc with clays. Nature. 214:852. Tiller, K.G., and J.F. Hodgson. 1962. The specific adsorption of cobalt and zinc by layer silicates. p393-403. In Proceedings of the 9th (1960) National Conference on Clays and Clay Minerals. Pergamon Press. New York. Tiller, K.G., an J.G. Pickering. 1974. The synthesis of zinc silicates at 200 C and atmospheric pressure. Clays Clay Miner. 22:409-416. Tiller, K.G., J. Gerth, and G. Brummer. 1984. The relative affinities of Cd, Ni, and Zn for different soil clay fractions and goethite. Geoderma. 34:17-35. - Tiller, K.G., J.L. Honeysett, and M.P.C. de Vries. 1972. Soil zinc and its uptake by plants. I. Isotopic exchange equilibria and the application of tracer techniques. Aust. J. Soil Res. 10:151-164 Travis, C.C., and E.L. Etnier. 1981. A survey of sorption relationships for reactive solutes in soil. J. Environ. Qual. 10:8-17. Turner, M.A., L.L. Hendrickson, and R.B. Corey. 1984. Use of chelating resins in metal adsorption studies. Soil Sci. Soc. Am. J. 48:763- 769. Tyler, L.D., and M.B. McBride. 1982. Mobility and extractability of cadmium, copper, nickel and zinc in organic and mineral soil columns. Soil Sci. 134:198-205. Udo, E.T., H.L. Bohn, and T.C. Tucker. 1970. Zinc adsorption by calcareous soils. Soil Sci. Soc. Am. J. 34:405-407. Vanselow, A.P. 1932. Equilibria of base exchange reactions of bentonites, permutites, soil colloids and zeolites. Soil Sci. 33:95-113. 117 Veith, J.A., and G.Sposito. 1977. On the use of the Langmuir equation in the interpretation of fadsorption" phenomenon. Soil Sci. Soc. Am. J. 41:697-702. Van Bladel, R., and H. Laudelout. 1967. Apparent irreversibility of ion- exchange reactions in clay suspensions. Soil Sci 104:134-137. Weigner, G. 1935. Ionen-umtausCh und Struktur. Trans. 3rd Int. Congr. Soil Sci. Oxford.