|
|
||||||||
a USDA-ARS, U.S. Salinity Lab., 450 W. Big Springs Rd., Riverside, CA 92507 USA
sgoldberg{at}ussl.ars.usda.gov
| ABSTRACT |
|---|
|
|
|---|
-4, B surface complexes are postulated, consistent with experimental spectroscopic results. Boron surface complexation constants for Al and Fe oxides and kaolinites are not statistically significantly different from each other. Boron surface complexation constants for kaolinites are statistically significantly different from those for 2:1 clays and soils. Boron surface complexation constants for 2:1 clays and soils are not statistically significantly different from each other, reflecting the dominance of 2:1 clay minerals in B adsorption reactions in arid-zone soils. Average sets of B surface complexation constants provided adequate descriptions of B adsorption behavior on all adsorbents studied, indicating some predictive capability. The constant capacitance model was able to predict B adsorption behavior on additional arid-zone soils using the average set of B surface complexation constants.
Abbreviations: ATR-FTIR, Attenuated Total Reflectance Fourier Transform Infrared SA, surface area suspension SD, suspension density
| INTRODUCTION |
|---|
|
|
|---|
Various modeling approaches have been used to describe B adsorption on soils and soil minerals. Historically, adsorption isotherm equations such as the Langmuir equation and the Freundlich equation provided descriptions of B adsorption on soil materials (e.g., Elrashidi and O'Connor, 1982; Goldberg and Forster, 1991). Adsorption isotherms are empirical models whose parameters are only valid for the particular conditions under which the experiment was conducted. More recently, various surface complexation models have been used to describe B adsorption on soils and soil minerals (Goldberg and Glaubig, 1985, 1986a, 1986b; Bloesch et al., 1987; Singh and Mattigod, 1992; Toner and Sparks, 1995). Surface complexation models are chemical models that define surface species, chemical reactions, mass balances, and charge balance. Surface complexation models contain molecular features that can be given thermodynamic significance (Sposito, 1983).
The constant capacitance model was developed to describe ion adsorption at the oxidesolution interface (Stumm et al., 1980) and has been used to describe B adsorption on oxides (Goldberg and Glaubig, 1985), clay minerals (Goldberg and Glaubig, 1986b), and soils (Goldberg and Glaubig, 1986a). Conventionally, surface complexation reactions for anion adsorption in the constant capacitance model are written in terms of undissociated acids (Stumm et al., 1980). The number of anion surface complexes considered is equal to the number of dissociations undergone by the acid. Because H3BO3 is a weak monobasic acid
, Goldberg and Glaubig (1985, 1986a, 1986b) considered one B surface complex in their model descriptions of B adsorption. For Al and Fe oxides and kaolinite minerals, good fits of the model to the data were obtained by optimizing only the B surface complexation constant (Goldberg and Glaubig, 1985). To adequately describe B adsorption on the 2:1 clay minerals montmorillonite and illite and soil samples, it was necessary to optimize values for the protonationdissociation constants along with that of the B surface complexation constant. Even so, chemically unrealistic values of the protonationdissociation constants were obtained in some cases, reducing the chemical significance of the model application to a curve-fitting procedure (Goldberg and Glaubig, 1986a, 1986b).
Recent work by Su and Suarez (1995) has provided direct spectroscopic evidence for the presence of B in both trigonal, B(OH)3, and tetrahedral, B
-4, geometries on the surface of amorphous Fe oxide. The authors analyzed mineral pastes using Attenuated Total Reflectance Fourier Transform Infrared (ATR-FTIR) spectroscopy. Because water was present in the samples, the results of Su and Suarez (1995) are directly applicable to natural soil systems.
The objectives of the present study are: (i) to reevaluate the ability of the constant capacitance model to describe B adsorption on oxides, clay minerals, and soils including microscopic experimental information, and (ii) to test the ability of the constant capacitance model to predict B adsorption on additional soils. The B surface complexes used in the model application will be constrained to those observed in the spectroscopic study of Su and Suarez (1995). Inclusion of experimental spectroscopic observations into surface complexation models maximizes the chemical significance of model parameters. Identical surface complexation reactions will be used for all adsorbing materials. The ability of the constant capacitance model to predict B adsorption on additional soils using the fitted parameters obtained with the improved modeling approach will also be investigated. Predictive ability is an important criterion for incorporation of adsorption models into chemical speciation and transport models.
| Materials and methods |
|---|
|
|
|---|
Boron adsorption isotherms were carried out on four additional soils using the procedure described by Goldberg and Glaubig (1985). Surface samples of the Porterville (fine, montmorillonitic, thermic Typic Chromoxerert), Wasco (coarse-loamy, mixed, nonacid, thermic Typic Torriorthent), Wyo (fine-loamy, mixed, thermic Mollic Haploxeralf), and Hesperia (coarse-loamy, mixed, nonacid, thermic Xeric Torriorthent) soil series consisted of the <2-mm fraction. The background electrolyte consisted of 0.1 M NaCl. Physical and chemical characteristics for these soils are provided in Goldberg et al. (1988). The ability of the constant capacitance model to predict B adsorption on these soils using average B surface complexation constants was evaluated. A detailed explanation of the theory and assumptions of the constant capacitance model is given in Goldberg (1992).
In the present application of the constant capacitance model to B adsorption, the following surface complexation reactions are included:
![]() | (1) |
![]() | (2) |
![]() | (3) |
![]() | (4) |
-4, and/or a positive surface site, XOH+2, could be written. Significant in this modeling approach is that the formation of both trigonal and tetrahedral B surface species is considered, consistent with the experimental spectroscopic results of Su and Suarez (1995).
Intrinsic equilibrium-constant expressions for the surface complexation reactions are:
![]() | (5) |
![]() | (6) |
![]() | (7) |
![]() | (8) |
Mass balance for the reactive surface functional group is:
![]() | (9) |
![]() | (10) |
has units of molc L-1. The computer program FITEQL 3.2 (Herbelin and Westall, 1996) was used to fit B surface complexation constants to the experimental data. FITEQL 3.2 uses a nonlinear least squares optimization technique to optimize equilibrium constants to experimental data. The FITEQL program contains the constant capacitance model of adsorption and can also be used as a chemical speciation model to evaluate predictions using previously determined equilibrium constants.
Input values for constant capacitance model parameters were the same as for the original applications (Goldberg and Glaubig, 1985 for oxides; Goldberg and Glaubig, 1986b for clays; Goldberg and Glaubig, 1986a for soils). The capacitance density was fixed at
. The protonation constant, log K+, and the dissociation constant, log K-, were set at: log
, log
for Fe oxides and log
, log
for Al oxides and clays, and log
, log
for soils. The total number of reactive hydroxyl groups, [XOH]T, was obtained experimentally from maximum B adsorption as was done in previous studies (Goldberg and Glaubig, 1985, 1986a, 1986b).
| Results and discussion |
|---|
|
|
|---|
|
|
|
|
|
|
|
![]() | (11) |
-Al2O3 (Alon, Cabot Corp., Boston, MA; Aluminium Oxid C, Degussa, Teterboro, NJ), and 9.9 for pseudoboehmite. These model-derived pKa values compare favorably with the surface pKa value of 8 to 9 for H3BO3 on amorphous Fe oxide postulated by Suarez et al. (1998) on the basis of analysis of the ATR-FTIR spectra obtained by Su and Suarez (1995). The model-derived pKa values are slightly higher than the pKa value for H3BO3 in aqueous solution (9.2), suggesting a somewhat greater proportion of B(OH)3 on the surface than predicted from its speciation in solution.
The fits obtained with the constant capacitance model (Fig. 1 and 2) are very similar in quality to those obtained by Goldberg and Glaubig (1985)(see Fig. 2 and 1, respectively) optimizing only the B surface complex, XOB(OH)2. Since the model predicts XOB(OH)2 to be the dominant B surface species until pH 9.7, addition of the surface species XOB
-3 provides only a slight improvement in fit. Nonetheless, the present modeling approach is considered superior since the B surface species postulated are consistent with experimental spectroscopic observations. Because values of the protonation and dissociation constants were also used as adjustable parameters in the earlier study, the fit of the model to B adsorption on amorphous Al oxide presented by Goldberg and Glaubig (1985)(Fig. 1d) describes the data more closely than the present model fit (Fig. 2d).
Figure 3 presents the fit of the constant capacitance model to B adsorption on three kaolinite clays. The model describes the data very well on all three materials. As expected when increasing the number of adjustable parameters, the fits in Fig. 3, using both XOB(OH)2 and XOB
-3, are improved over those presented by Goldberg and Glaubig (1986b)(Fig. 1) using XOB(OH)2. Values of the B surface complexation constants and average values of log KB and log KB- for kaolinites are presented in Table 1. Average values of the B surface complexation constants for kaolinites are not statistically significantly different from those for Al and Fe oxides at the 95% level of confidence. Dashed lines in Fig. 3 indicate the model predictions of B adsorption obtained using an average set of B surface complexation constants for kaolinites. The model predictions describe the B adsorption data equally well as the model fits, indicating that the average set of surface complexation constants can be used to describe B adsorption behavior on other kaolinite minerals. Surface pKa values for reaction [11] are: 9.8 for well crystallized kaolinite, 10.1 for poorly crystallized kaolinite, and 9.7 for Georgia kaolinite. These values are comparable in magnitude to those for Al and Fe oxides, reflecting the similarity in values of the B surface complexation constants of all of these materials.
The ability of the constant capacitance model to describe B adsorption on the 2:1 clay minerals, montmorillonite and illite, is presented in Fig. 4. Overall, the model was able to describe B adsorption on 2:1 minerals, fitting a smooth curve through the data. The shape of the model fit is parabolic in nature, leading to some underprediction of the data near the pointed adsorption maxima. For SWy-1 montmorillonite it was necessary to optimize the protonationdissociation constants along with the B surface complexation constant to obtain a fit to the data. Values of the B surface complexation constants and average values of log KB and log KB- for 2:1 clays are provided in Table 1. Average values of the B surface complexation constants for 2:1 clays are statistically significantly different from those for kaolinites at the 95% level of confidence. Model predictions using the average set of B surface complexation constants are indicated in Fig. 4 by dashed lines. The model predictions represent the data less closely than the model fits but provide an adequate description. These results suggest some predictive utility of average sets of B surface complexation constants for describing B adsorption behavior on additional 2:1 clay minerals.
The model fits presented in Fig. 4 are approximately similar in quality to those obtained by Goldberg and Glaubig (1986b)(Fig. 2 and 3). It is significant to note that the model fits in the present study were obtained by optimizing constants for two B surface species XOB(OH)2 and XOB
-3-, while the model fits of Goldberg and Glaubig (1986b) were obtained by optimizing three constants, XOB(OH)2 for the B surface species and values for the protonation constant log K+, and the dissociation constant log K-. In the optimization process, Goldberg and Glaubig (1986b) obtained unrealistic values of log K+ and log K- for one montmorillonite and one illite. The present modeling approach is considered superior because: (i) one less adjustable parameter is optimized; (ii) identical values of log K+ and log K- are used for all 2:1 clays; (iii) unrealistic values of log K+ and log K- are avoided. Surface pKa values for adsorbed B are: 7.5 for SAz-1 montmorillonite, 7.4 for STx-1 montmorillonite, 7.0 for SWy-1 montmorillonite, 7.1 for IMt-1 illite, 8.3 for Fithian illite, and 8.5 for Morris illite. These pKa values are significantly lower than those for Al and Fe oxides and kaolinite clays, suggesting a much greater role for the surface species XOB
-3 in B adsorption on 2:1 clays.
The ability of the constant capacitance model to describe B adsorption as a function of solution pH on 14 arid-zone soil samples is indicated in Fig. 5. The model provides a quantitative description of B adsorption behavior on all soils except for Altamont soil and subsoil, Arlington subsoil, and Imperial subsoil. For these samples, the model describes the data qualitatively. For Arlington subsoil, and Imperial soil and subsoil it was necessary to optimize the protonationdissociation constants along with the B surface complexation constant to obtain a fit to the data; for these soils log KB(int) did not converge. Table 1 provides values of the B surface complexation constants and average values of log KB and log KB- for soils. Average values of the B surface complexation constants for soils are not statistically significantly different from those of 2:1 clays at the 95% level of confidence. This result is consistent with the fact that clay minerals are important B-adsorbing sites in arid-zone soils (Elrashidi and O'Connor, 1982).
The model fits obtained in Fig. 5 are comparable with those obtained for these soils by Goldberg and Glaubig (1986a). As for the 2:1 clay minerals, it is significant to note that the present model fits were obtained by optimizing two surface constants (KB and KB-), while the model fits of Goldberg and Glaubig (1986a) were obtained by optimizing three surface constants (KB, K+, and K-). In carrying out their optimization, Goldberg and Glaubig (1986a) obtained unrealistic values of log K+ for one soil and of log K- for five soils. As for 2:1 clays, the modeling approach of the present study is considered superior because one less parameter is adjusted and identical log K+ and log K- values are used for all soils, thus avoiding unrealistic values of these parameters. Surface pKa values for adsorbed B range from 6.8 for Pachappa subsoil to 9.0 for Altamont soil. The distribution of B surface species for these two soils are presented in Fig. 6. For the Pachappa subsoil the XOB
- 3species dominates across most of the pH range, while for the Altamont soil XOB(OH)2 is dominant across most of the pH range.
Dashed lines in Fig. 5 indicate model predictions of B adsorption obtained using an average set of B complexation constants for soils. The model predictions represent the data slightly less closely than the model fits but still provide a reasonable description. The ability of the average set of B surface complexation constants to predict B adsorption behavior was tested on four additional arid-zone soils and is presented in Fig. 7 . The model was able to predict B adsorption behavior on all four soils, fitting a smooth parabolic curve through the data. These results indicate the utility of the average set of B surface complexation constants for predicting B adsorption on other arid-zone soils. The capability of these constants to describe B adsorption on soils of higher oxide and organic matter content remains to be evaluated.
|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
S. Goldberg, L. J. Criscenti, D. R. Turner, J. A. Davis, and K. J. Cantrell Adsorption Desorption Processes in Subsurface Reactive Transport Modeling Vadose Zone J., August 1, 2007; 6(3): 407 - 435. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Goldberg, D. L. Corwin, P. J. Shouse, and D. L. Suarez Prediction of Boron Adsorption by Field Samples of Diverse Textures Soil Sci. Soc. Am. J., August 4, 2005; 69(5): 1379 - 1388. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Goldberg Modeling Boron Adsorption Isotherms and Envelopes Using the Constant Capacitance Model Vadose Zone J., May 1, 2004; 3(2): 676 - 680. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Goldberg, D. L. Suarez, N. T. Basta, and S. M. Lesch Predicting Boron Adsorption Isotherms by Midwestern Soils using the Constant Capacitance Model Soil Sci. Soc. Am. J., May 1, 2004; 68(3): 795 - 801. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Matsi and V. Z. Keramidas Alkaline Fly Ash Effects on Boron Sorption and Desorption in Soils Soil Sci. Soc. Am. J., July 1, 2001; 65(4): 1101 - 1108. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Goldberg, S. M. Lesch, and D. L. Suarez Predicting Boron Adsorption by Soils Using Soil Chemical Parameters in the Constant Capacitance Model Soil Sci. Soc. Am. J., July 1, 2000; 64(4): 1356 - 1363. [Abstract] [Full Text] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Crop Science | |||
| Vadose Zone Journal | Journal of Plant Registrations | ||||
| Journal of Natural Resources and Life Sciences Education |
Journal of Environmental Quality |
||||