SciELO - Scientific Electronic Library Online

vol.34 número1El depósito de psefítas "El Peladero", en los acantilados del litoral marino de la provincia de Río Negro: Argentina índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados



  • No hay articulos citadosCitado por SciELO

Links relacionados

  • En proceso de indezaciónCitado por Google
  • No hay articulos similaresSimilares en SciELO
  • En proceso de indezaciónSimilares en Google



versión On-line ISSN 1852-7744

Geoacta v.34 n.1 Ciudad Autónoma de Buenos Aires ene./jun. 2009


Fluoride retardation from quartz sand-packed column tests

Eduardo Usunoff, Pablo Weinzettel, and Sebastián Dietrich

CIC-UNCPBA, Rep. de Italia 780, Azul
CIC-UNCPBA, Rep. de Italia 780, Azul,
Becario CONICET, Rep. de Italia 780, Azul,


Inasmuch as both low and high concentrations of F- in groundwater have different detrimental effects on human health (increased dental caries, and mottled enamel or even severe structural bone deformations, respectively), many efforts have focused on the movement of such anion in aqueous systems. It is so because water drinking is the main intake of F- by humans. This paper presents the results of seven dynamic experiments in which solutions of varying [F-], pH, and flow velocities circulated through columns packed with clean quartz sand. The breakthrough data were analyzed by means of a computer code adapted to the estimation of equilibrium and non-equilibrium solute transport parameters from miscible displacement experiments in a steady-state, uniform flow field using a pulse-type or continuous source. It was found that larger retardation factors (R) are associated with low pH, low [F-] in the feed solutions, and larger flow velocities. Such results appear to be related to the form of the F species, the rather weak bond between the adsorbate and the quartz sand, and stronger anion repulsion at low pore velocities. The estimated values for R vary between 1,22 and 1,50, whereas the distribution coefficients were in the range of 0,1 to 0,05 L kg-1. It should be said that the breakthrough curves display hysteresis, leading to a desorption behavior that proceeds faster than the adsorption process.

Keywords: Column experiments; F- adsorption/desorption; Retardation


Puesto que tanto valores bajos como altos de F- en las aguas subterráneas tienen efectos adversos en la salud de los humanos (mayor cantidad de caries, y manchas del esmalte dental e incluso deformaciones óseas, respectivamente), han sido muchos los esfuerzos destinados a dilucidar el movimiento de especies fluoradas en ambientes acuosos. Ello es así porque el consumo de agua por parte de seres humanos representa la mayor fuente de ingreso de F-. Este trabajo presenta los resultados de siete experimentos dinámicos en los que soluciones de diferentes [F-], pH y velocidades de flujo percolaron a través de columnas empacadas con arena cuarzosa. Los datos de la curva de salida fueron analizados con un programa que permite la estimación de parámetros del transporte de solutos en condiciones de equilibrio o no equilibrio a partir de experimentos de columnas, en un campo de flujo homogéneo y en estado estacionario, creado por un impulso acotado o continuo del trazador de interés. Se halló que los mayores valores del factor de retardo (R) se asociaron a bajos valores de pH, bajas concentraciones de F- de las soluciones de inicio, y mayores velocidades de flujo. Los resultados se relacionaron con la forma de las especies de F, el débil enlace entre la especie adsorbida y la superficie de los granos de la arena cuarzosa, y la mayor repulsión aniónica a bajas velocidades porales. Los valores estimados de R variaron entre 1,22 y 1,50, en tanto que los coeficientes de distribución estuvieron en el rango de 0,1 a 0,05 L kg-1. Debe indicarse que las curvas de salida muestran cierto efecto de histéresis, que hace que el fenómeno de desorción se produzca más rápidamente que el de adsorción.

Palabras clave: Experimentos de columna; Adsorción/desorción de F-; Retardo


It has been known for a long time that the fluorine (F) content in water supplies and diets has a profound effect in the development of human bones and teeth. The intake of F is mainly determined by the concentration of fluoride ions (F-), the most conspicuous fluorine-carrying dissolved species, in the local water supply. One important fact is that both low and high concentrations of F- may affect human health. Indeed, F- acts as an inhibitor of the bacterial enzymes which, in turn, produce the acid that initiate the enamel caries (Rahamatulla, 1977). Hodge and Smith (1965) show that the absence of F- in drinking water or its presence in low concentrations (below 0,5 mg L-1) increases remarkably the number of dental caries in children. Gómez Artola et al. (1983) found that the low F- content in public drinking water explains the epidemic dental caries in Madrid, Spain. Fabiani et al. (1995) demonstrated that water consumption with an average [F-] of 1,54 mg L-1 reduces notably the risk of osseous fractures, with a considerable saving of public money that can be invested in other sanitary treatments, and that offers a valid alternative to the conventional treatment of the osteoporosis disease. On the other hand, the ingestion of excess F- (1-4 mg L-1) is responsible for the disease known as fluorosis (mottled enamel), characterized by a brownish-yellow discoloration of teeth (Hodge and Smith, 1965). Waters containing high F- (above 4-5 mg L-1), if consumed for a long time, can produce severe structural deformations of human bones (Krishnamachari, 1977).
Assuming that most of water supplies rely on groundwater resources, and given the limitations pointed out above, the purpose of this paper is to present and discuss dynamic laboratory experiments data where solutions of varying [F-], pH, and v (mean pore velocity) percolated columns packed with clean quartz sand. The final goals were to observe the adsorption-desorption behavior, and to estimate the retardation factor, R (i.e., a dimensionless quantity that indicates whether the solute travels with the same velocity of the water).

MATERIALS AND METHODS Quartz sand preparation

To eliminate impurities, the quartz sand samples, commercially known as bond sand, were boiled in 10% HCl for ten minutes and then washed with excess distilled water until the pH of the leachate equalled that of the distilled water. The samples were then oven-dried at 105 ºC and kept in close containers.

[F-] and pH measurements

The F- concentrations were measured with a Fisher Accumet pH/mV (model 620) meter and a specific electrode (Orion, model 94-09). Before measuring a given set of samples, fresh solutions of known [F-] were prepared to obtain the calibration curve, and a recheck of such a curve was done between samples measurements, following the guidelines given by Orion Research Inc. (1973) and by the U.S. Environmental Protection Agency (1979). Both samples and standards were mixed in proper amounts with TISAB, a commercially available solution designed to partially buffer the pH between 5 and 6, and to provide a high total ion strength background that swaps out variations in total ion strength between samples (Orion Research Inc., 1973). For the pH measurements, the pH meter mentioned above was used with a standard glass-combination electrode, which was calibrated against commercial buffer solutions of known pH. All [F-] and pH readings were taken at 23 ± 1 ºC.

Column preparation

A plexi-glass column was packed with the quartz sand above mentioned. Because of the many precautions involved in the packing tasks and the tests needed before launching the experiments, the reader is referred to Usunoff (1993) for further details. Suffice to say that the following data were measured/obtained:

Column length: 30,48 cm
Column internal diameter: 3,15 cm
Mass of quartz sand in column: 369,5 g
Θ (porosity): 39% (estimated at packing time)

Dimensionless variables
Ci/C0 (relative concentration): ratio between the outlet concentration at time i and the input concentration of tracer (C0)
T (pore volume or PV) = v t / L, where v is the mean pore velocity (L T-1), t is the time for C/C0 = 0,5, and L is the column length
Experimental mechanics
The flow, which took place from bottom to top, was driven by a peristaltic pump connected through a two-way valve to two containers with fluids at the same pH (one carrying the F- solution, and the other distilled water).
Typically, a run proceeded as follows: distilled water was circulated until trie effluent pH stabilized, time at which the two-way valve was switched and the circulation of trie tracer solution (F) started. Samples of the effluent were taken until their [F] was about the same as the input [F]. Then, the double-position valve was switched to allow the entrance of the pH-adjusted distilled water. Sampling continued until the effluent [F] was negligible. Such experiments yielded data to construct the BTC (breakthrough curve). Breakthrough curve data handling
The adsorption-desorption data were the input of the computer program CFITIM (van Genuchten, 1981), which was later included as one of the options of the more robust set of programs STANMOD (Simunek et al., 1999). CFITIM is a program for estimation of equilibrium and non-equilibrium solute transport parameters from miscible displacement experiments in a steady-state, uniform flow field using a pulse-type or continuous source. Such estimations use a non-linear, least-squares curve fitting procedure.
If the relationship between the solute adsorbed and its concentration in solution can be described by a linear isotherm, the one-dimensional convective-dispersive equation of solute transport in saturated porous media can be written as:


where R is the retardation factor (dimensionless), C is the solute concentration (ML3), t is the time (T), D is the longitudinal dispersión coefficient (L2T_1), x is the space (L), and v is the pore velocity (LT1). The retardation factor R is given by (Freeze and Cherry, 1979):


where pb is the bulk mass density (ML3), and Kd is the distribution or partition coefficient (L3M_1). One of the advantages of using CFITIM is that it estimates R, which can be plugged directly into equation 1. Moreover, given that valúes of pb and 0 are available, Kd can be calculated. Notice that if R = 1, the solute travels with the same velocity as the carrier fluid, whereas if R > 1, the solute interacts with the solid, and its movement is said to be retarded with respect to the bulk flow.


Seven runs were carried out, numbered 1 to 7, in which the [F] and the pH of the input solutions were varied, as well as the flow velocity. Such valúes are shown in Table 1.

Table 1. Column experiments: measured and estimated data.

As for the breakthrough data (measured and fitted), their graphical display is shown in Figures 1 through 7.

Figure 1. Breakthrough curve for column experiment 1 (dots: data; solid line: fitted by CFITIM).

Figure 2. Breakthrough curve for column experiment 2 (dots: data; solid line: fitted by CFITIM).

Figure 3. Breakthrough curve for column experiment 3 (dots: data; solid line: fitted by CFITIM).

Figure 4. Breakthrough curve for column experiment 4 (dots: data; solid line: fitted by CFITIM).

Figure 5. Breakthrough curve for column experiment 5 (dots: data; solid line: fitted by CFITIM).

Figure 6. Breakthrough curve for column experiment 6 (dots: data; solid line: fitted by CFITIM).

Figure 7. Breakthrough curve for column experiment 7 (dots: data; solid line: fitted by CFITIM).

Runs 1, 2, and 3 were designed to assess the relationship between F- adsorption and pH. Runs 4, 5, and 6 were to test the degree of F- adsorption as the [F-] in the feed solutions was varied. Finally, runs 1, 4, and 7 were intended to reveal the influence of the mean pore velocity in the adsorption phenomenon.


Comments on the retardation factor values

As presented in the Material and Methods section of this paper, the estimation of the retardation factor R made use of the fitting capabilities of the program CFITIM (van Genuchten, 1981) to BTC data. Although it is a well-known program developed by a notorious scientist, the reader may like to have an independent check on how trustable those values are. Bovin (2008) indicated that when the BTC data are plotted in the conventional way of C/C0 (relative concentration) vs. dimensionless time T (pore volume), the T value on the x-axis is equal to R for C/C0 = 0,5. By looking at Figures 1 to 7 one may conclude that the so-estimated values for R resemble quite well the fitted values shown in Table 1, which reinforces the suitability of Bovin´s rule of thumb.

[F-] and flow velocity (almost) constant, pH varied

Taking the retardation factor R as the diagnostic variable, R is seen to increase as the pH of the input solutions decreases. Notice, also, the slight increment of the effluent pH. That is to say, low pH values would enhance the F adsorption on the quartz sand particles. Such observations may be related to at least two factors:

(1)  The zero point charge of quartz, for which Stumm and Morgan (1981) have given a value of 2. Indeed, as the pH increases above 2, the surface of quartz particles becomes negatively charged, with a corresponding decrease in anion sorption.
(2)  The fluorine complexation as the pH is lowered. At pHs below neutrality and in the absence of cations such as Al3+ and Fe3+, the F species able to be present are HF0, HF-2, and F- (Hem, 1968; Roberson and Barnes, 1978). The concentrations of HF0 and HF-2 (at 25 °C) can be calculated from (Broene and De Vries, 1947): [HF0] = [H+] [F-] 103,17, and [HF-2] = [H+] [F-]2 103,76
Hence, unless the pH is extremely low and the total [F-] is very high, the amount of HF-2 in solution can be neglected. Of importance, however, are HF0 and F-, whose percentages for runs 1, 2, and 3 were as follows:

Run 1: almost all added F is in the F- form
Run 2: [HF0] ≈ 93% of total F, and [F-] ≈ 7% of total F
Run 3: [HF0] ≈ 1,5% of total F, and [F-] ≈ 98,5% of total F

As for the possible species that F can form with silicon, Roberson and Barnes (1978) concluded that SiF-6 is the most important one. However, the solubility of quartz at 25 °C is so low (6 mg L-1) (Morey et al., 1962; Marion et al., 1976; Roberson and Barnes, 1978) that, even for the column run at the lowest pH (run 2), the [SiF-6] barely exceeds 10-8 M, and it can safely be disregarded. Based on (1) or (2) above, or on both, the conclusion is that the F- adsorption depends on the type of F species present, which, in turn, are linked to the prevailing pH in the system, with low pH values (below 4,5 or 5) promoting the formation and adsorption of HF0. It may be pertinent to recall that Hingston (1981) suggested the consideration of the role of both anion (F-) and conjugated acid (HF) in adsorption studies.

pH and flow velocity (almost) constant, [F-] varied

The observed trend indicates that adsorption occurred, with R values becoming smaller as the input [F-] was increased (Table 1). Such results are rather puzzling because, in general terms, the adsorption of inorganic species onto mineral surfaces has been found to increase with the concentration of adsorbate in solution. Inasmuch as [HF0] should have been minimal at the pH of the three input solutions, the conclusion is that the F- ion did not interact much with the mineral surface, or that the adsorption process was overcome by the highly dynamic nature of the experiment. It may be interesting to add that Fluhler et al. (1982) reported that an earlier breakthrough of F- at the top of columns as the [F-] of the feed solutions was increased could be interpreted as the non-linear nature of the adsorption isotherm.

pH and [F-] (almost) constant, flow velocity varied

Although care was taken to select velocities so that the flow would be laminar, the results indicate that larger retardation factors are associated with faster flow velocities (Table 1). Again, this finding is contrary to the expected outcome. It can be thought, at least theoretically, that an enlarged time of contact between the adsorbate and the adsorbent (i.e., lower velocities) would result in better chances for the species in solution to attach to the mineral surface. Such an explanation, however, ignores the status of the surface charge of the adsorbent when adsorption takes place. As stated earlier, the quartz particles are negatively charged above a pH of 2, which would certainly inhibit anion adsorption. If so, and in the light of the results already discussed, such repulsion seems to be more effective at lower pore velocities.

Changes in pH and the adsorption-desorption behavior

Two common characteristics were found in the column experiments, namely, an increase in the pH of the effluent solutions as compared to the pH of the input solutions (Table 1), and the existence of hysteresis in the otherwise reversible adsorption process. Allowance was not made to account for the buffering effect of quartz, although it is assumed to be negligible at the pHs used in the experiments. Therefore, the observed systematic increase in pH can be ascribed to the release of OH- from the exchange sites on the sand particles.
The recession limb of the BTC indicates that, for all cases, desorption proceeded at a faster pace than adsorption. Such an effect may represent the lack of equilibrium during washing (Hingston et al., 1974) and, at the very least, gives a clear indication of the weak bond between the adsorbed species and the sand particles.

The distribution or partition coefficient Kd

Inasmuch as the distribution coefficient Kd is embodied in the definition of the retardation factor R (equation 2), the fitted values for the latter can be used to calculate the magnitude of the former. If the assumption is made that the mass density of the quartz sand is 2,65 g cm-3, then the bulk mass density ρb that corresponds to a porosity θ of 39% is around 1,62 cm-3, and the distribution coefficient can be calculated from the expression for R. Having done so, the figures shown in Table 2 were obtained for the various runs.

Table 2. Estimated distribution coefficients for the column experiments.

The underlying assumption in estimating Kd is that the values for both R and θ are correct and, hence, not subject to uncertainties. However, that may not be strictly so. If, instead a porosity of 39%, a range of porosities between 35 and 45% is considered (i.e., a change of 29%) with other variables remaining constant, a 52% change in the values for Kd is calculated. On the other hand, when θ is kept constant and R is varied between 1,1 and 1,6 (i.e., a 45% change), the corresponding values for Kd change by 400%. Such a simple exercise demonstrate quite explicitly the important effect that uncertainties in the measured values of θ and the fitted values for R can have on the estimation of the distribution coefficient. Extreme care, therefore, should be taken by those applying values for Kd reported in the literature, particularly when those values are to be used in environmental-impact studies.


Solutions of variable pH, concentration of F-, and mean flow velocity were forced to circulate through columns packed with quartz sand, which provided data for estimating the retardation factor R (fitted by means of a non-linear, least-squares curve fitting procedure). The outlet measurements indicate that a weak interaction between the solute and the quartz sand takes place, and that the early breakthrough of outflowing F- (hysteresis) could point out to a lack of equilibrium during the displacement of the adsorbate species.
It would appear that lower pH input solutions enhance the adsorption of F species, most probably in the HF0 form. Surprisingly enough, increased amounts of F- in the input solutions did not result in increased adsorption, which indicates a hardly noticeable interaction between the adsorbent and the adsorbate and the non-linear nature of the adsorption phenomenon. Increased flow velocities resulted in larger retardation factors, suggesting that longer contact times between the adsorbent and the adsorbate promote anion repulsion.
The estimated values for R vary between 1,22 and 1,50, whereas the distribution coefficients were in the range of 0,1 to 0,05 L kg-1.


1. Bovin, T. 2008. Environmental Hydrogeology. Lecture Notes. Downloaded on March 18, 2008 from:        [ Links ]

2. Broene, H. and T. De Vries. 1947. "The thermodynamics of aqueous hydrofluoric acid". Journal of the American Chemical Society 69, pp. 1644-1646.        [ Links ]

3. Fabiani, L., M. Vitali, M. Parafati, and S. Rodolico. 1995. "Fluoride in waters as a protective factor for bone fractures: Preliminary data of an epidemiological study in Italy". IAHS Publication 223, pp. 135-140.        [ Links ]

4. Fluhler, H., J. Polomski, and P. Blaser. 1982. "Retention and movement of fluoride in soils". Journal of Environmental Quality 11, N° 3, pp. 461-468.         [ Links ]

5. Freeze, R. and J. Cherry. 1979. Groundwater. Prentice-Hall, 604 p.        [ Links ]

6. Gómez Artola, C., J. Borregón Martínez, M. Llamas Madurga y J. Sánchez Sáez. 1983. "Análisis preliminar del contenido de flúor de las aguas de Madrid y de su incidencia en las caries dentales de los niños". Actas del III Simposio Nacional de Hidrogeología, pp. 157-166.        [ Links ]

7. Hem, J. 1968. "Graphical methods for studies of aqueous aluminum hydroxide, fluoride, and sulfate complexes". Geol. Survey Water-Supply Paper 1827-B, 33 p.        [ Links ]

8. Hingston, F. 1981. "A review on anion adsorption". In Anderson, M. and A. Rubin (Editors). 1981. Adsorption of Inorganics at Solid-Liquid Interfaces. Ann Arbor Science, pp. 51-89.        [ Links ]

9. Hingston, F., A. Posner, and J. Quirk. 1974. "Anion adsorption by goethite and gibbsite. II. Desorption of anions from hydrous oxide surfaces". Journal of Soil Science 25, N° 1, pp. 16-26.        [ Links ]

10. Hodge, H. and F. Smith. 1965. Fluorine Chemistry, Volumen IV". Academic Press, 789 p.        [ Links ]

11. Krishnamachari, K. 1977. "Some new aspects of fluorosis in South India. A new approach to prevent the problem". Proceedings of the 1974 Symposium on Fluorosis, Indian Academy of Geosciences, pp. 519-529.        [ Links ]

12. Marion, G., D. Hendricks, G. Dutt, and W. Fuller. 1976. "Aluminum and silica solubility in soils". Soil Science 121, N° 2, pp. 76-84.        [ Links ]

13. Morey, G., R. Fournier, and J. Rowe. 1962. "The solubility of quartz in the temperature interval from 25 °C to 300 °C". Geochimica et Cosmochimica Acta 26, pp. 1029-1043.        [ Links ]

14. Orion Research, Inc. 1973. Instruction manual. Fluoride electrodes models 96-09 and 94-09. Orion Res., Inc., 28        [ Links ]

15. Rahamatulla, Md. 1977. "Dental fluorosis as a dental health problem". Proceedings of the 1974 Symposium on Fluorosis, Indian Academy of Geosciences, pp. 487-495.        [ Links ]

16. Roberson, C. and R. Barnes. 1978. "Stability of fluoride complex with silica and its distribution in natural water systems". Chemical Geology 21, pp. 239-256.        [ Links ]

17. Simunek, J., M. Th. van Genuchten, M. Sejna, N. Toride, and F. J. Leij, 1999. "The STANMOD computer software for evaluating solute transport in porous media using analytical solutions of convection-dispersion equation". Versions 1.0 and 2.0, International Ground Water Modeling Center, Colorado School of Mines, Golden, Colorado, 32p.        [ Links ]

18. Stumm, W. and J. J. Morgan. 1981. Aquatic Chemistry. An Introduction Emphasizing Chemical Equilibria in Natural Waters. John Wiley and Sons, New York, New York.        [ Links ]

19. U.S. Environmental Protection Agency. 1979. Methods for chemical analysis of water and wastes, 1978. Env. Mon. and Support Lab., Office Res. and Develop., 441 p.        [ Links ]

20. Usunoff, E. 1993. "Aspectos prácticos de experimentos con columnas de laboratorio". CICPBA, Situación Ambiental de la Provincia de Buenos Aires, III(24):21-33.        [ Links ]

21. van Genuchten, M. Th. 1981. "Non-equilibrium transport parameters from miscible displacement experiments". U.S. Salinity Lab. Research Report 119, 88 p.        [ Links ]