Abstract
In the present work, an economical, more selective and a bio-compatible fluoride removal technique has been developed by using zirconium impregnated
Introduction
Fluoride prevents dental caries at lower concentrations but at concentrations above 1.5 ppm, it causes serious health problems to humans and animals through drinking water (Ayoob and Gupta, 2006; Maiti et al., 2011). Fluoride enters the groundwater due to slow dissolution of fluorine-containing rocks (Agarwal et al., 2003; Bhatnagar et al., 2011; Sujana et al., 2009). Industries like glass and ceramic, semiconductor, electroplating, coal-fired power stations, beryllium extraction plants, brick and iron works and aluminum smelters increase the fluoride contamination in groundwater (Bhatnagar et al., 2011).
Defluoridation of water using both physico-chemical and biological methods together have proven advantageous when compared to specific methods (Mekonen et al., 2001). Researchers have employed plain algae
In this work, we prepared an eco-friendly, novel biosorbent by impregnating zirconium (IV) onto algae, and investigated its fluoride adsorption behaviour. Besides, the effects of key operation parameters, such as initial fluoride concentration, fluoride solution pH, biosorbent dose, contact time and coexisting anions were optimized using a central composite design (CCD) based on response surface methodology (RSM). The RSM is an experimental modelling system that evaluates the relation between the experimentally controlled group of variables and an obtained response (Vasconcelos et al., 2000). In addition, the isotherms and kinetics of the adsorption process were studied, and the possible mechanism of fluoride removal by the adsorbent was also proposed based on the Fourier transform infrared (FT-IR) spectroscopy, scanning electron microscope–energy dispersive spectroscopy (SEM-EDS) and X-ray photoelectron spectroscopy (XPS) analysis results.
Experimental procedure
Materials
Analysis of fluoride ion concentration
The concentration of fluoride ion in the solution was determined by Mettler Toledo fluoride ion selective electrode (perfectION™ combined fluoride electrode make) and Mettler Toledo ion analyzer (SevenCompact pH/ion meter S220 make). A 5 mL of TISAB III (total ionic strength adjustment buffer, pH 5–5.5) was added to 50 mL fluoride solutions to prevent the interference of other ions with the fluoride measurements.
The instrument was calibrated using fluoride standards provided by Mettler Toledo to make working solutions of required concentration. These solutions were used to calibrate the instrument to measure fluoride concentrations.
Procurement and culturing of algal strains
Composition of Bold’s basal medium.
Composition of Bold’s basal medium (modified).
Composition of trace metal solution.
Cultivation of algae and preparation of the algal powder
The 300 mL of sterile medium in 1-L Erlenmeyer flasks were inoculated with algal cells and placed inside an illumination box at 24℃. The light intensity was maintained in the box by using a suitable number of fluorescent tube lights
Preparation of algal biosorbents by zirconium impregnation
In this study, the powdered algal species were doped with zirconium by adding 5% ZrOCl2.10H2O solution to the algal powder in the ratio of 3:1 (wt ratio). This ratio was fixed based on the optimum removal efficiency of adsorbent as indicated in online Figure SI-1. The mixture was equilibrated for 72 h at 25℃. The zirconium-doped algae were then filtered, washed with water to remove free zirconium ions, dried in an oven at 100℃ and subsequently used to investigate its defluoridation capacity (Rajan and Alagumuthu, 2013).
Characterization of zirconium-impregnated algal biomass
The chemical groups of the samples before and after fluoride biosorption were characterized using FT-IR spectroscopy (Spectrum Two™ from PerkinElmer, USA). The surface morphology of the samples was observed using scanning electron microscope (SEM, Ultra 55 FESEM model from Carl Zeiss). Elemental analysis of the samples was carried out using energy dispersive spectroscopic (EDS) detector from Oxford Instruments. Kratos Axis Ultra spectrometer (UK) with Al Kα anode radiation source was used for XPS studies. The binding energies in XPS studies were referenced to the C1s peak of the surface adventitious carbon at 284.8 eV. Mastersizer 2000 (Malvern Instruments, Inc., UK) was used to determine the particle size of the biosorbent. Brunauer–Emmett–Teller (BET) surface area was analysed using TriStar 3000 V6.05 A.
Zeta potential measurement
The zeta potential (ζ) of the surface was found as per a previously reported procedure (Zhang, 2005) for 10 g/L of biosorbent suspension with 100 ppm fluoride and without fluoride, in the pH range of 2–10 using a Zetasizer 2000 (Malvern Instruments, Inc.). The 0.01 mol/L NaNO3 was added as the background electrolyte and aged for 24 h at 25℃. The initial pH was adjusted using HNO3 or NaOH. Biosorbent suspensions with or without F− at the preferred pH were shaken at 25℃ and 180 r/min for 24 h. The final pH was analysed, and then the suspension was introduced into the electrophoretic cell for determining the zeta potential values. The pH at the point of zero charge was determined by interpolating the zeta potential data to the zero potential (Dou et al., 2012).
Batch biosorption experiments
The batch studies were conducted for the optimization of pH, initial fluoride concentration, biosorbent dosage and contact time. The experiments were performed at 30℃ in 250-mL polypropylene flasks. The initial pH of the samples was maintained using 0.1 M HCl or 0.1 M NaOH. A known amount of biosorbent was added to the feed solution, and the samples were shaken in a rotary shaker incubator for a known period at 190 r/min. The samples were then filtered using Whatman filter paper no. 42 and the clear solution was analysed to determine the residual fluoride concentration using the fluoride ion selective electrode.
The percentage biosorption was calculated using equation (1).
Experimental design and data analysis
The experimental factors and levels for the CCD.
The percentage fluoride removal was fitted into equation (2) and ANOVA was also performed to obtain a correlation between the independent variables and the response. The coefficient of determination,
Equilibrium studies
Adsorption isotherm experiments were carried out in 250-mL polypropylene flasks on 100 mL of solution containing fixed initial fluoride concentrations as 100 ppm. The adsorbent dose was varied from 1 to 10 g/L. The pH was adjusted and held at pH 2. These flasks were agitated at 190 r/min and maintained at 30℃ for 4 h. The solutions were filtered, and analysed for the residual fluoride concentrations.
Adsorption kinetics
A 10 g sample of the algal biosorbent in 250-mL polypropylene flasks containing 100 mL of fluoride solution (100 ppm) were mixed and placed in a temperature-controlled rotary shaker maintained at 190 r/min and sampled at different time intervals. The biosorbent was finally removed by filtration and the fluoride concentration was determined. The contact time required for complete fluoride biosorption was determined for each biosorbent.
Desorption studies
To increase the feasibility and economy of the biosorption method, the regeneration of the biosorbents were studied. The adsorption–desorption cycles of fluoride were repeated four times with 1% NaOH, distilled water and 0.1 N HCl. In batch desorption experiments, 10 g of fluoride-loaded biosorbent in 250 mL Erlenmeyer flasks was contacted in series with 100 mL of 1% NaOH, 80 mL of distilled water and 100 mL of 0.1 N HCl solutions consecutively, at room temperature (27℃ ± 2℃). Each mixture was agitated on orbital shaker at 160 r/min for 30 min. The biosorbent was removed and supernatant was analysed for fluoride concentration by ion selective electrode.
Results and Discussion
Characterization of the biosorbents
Particle size of the algal biosorbents.
BET surface area and specific gravity.
Figure 1(a) and (b) depict the SEM images of zirconium-doped SEM images of zirconium-doped EDS of EDS of 


The FTIR spectra of zirconium-doped FTIR spectra of zirconium-doped (a) 
The XPS wide spectra of the zirconium-doped XPS plot of (a) 
The new species responsible for the peak shifting was determined by using peak deconvolution method on the Zr3d spectra before and after adsorption. The method reported by Zhang et al. (1995) was used by employing a computer algorithm which varied the relative intensity (peak area) and binding energy position of two identical Zr3d peaks to generate a ‘best fit’ to the experimental data. The resolved spectra are given in Figure 6(a) to (d) and summarized in Table 7. During the fitting process, the peak area ratio remained identical for the two corresponding doublets, that is A(Zr3d3/2 doublet for Zr-F)/A(Zr3d3/2 doublet for Zr-O) = A(Zr3d5/2 doublet for Zr-F)/A(Zr3d5/2 doublet for Zr-O) (where A means peak area).
Deconvoluted XPS spectra of Zr3d and P2s for (a) The deconvolution of Zr3d and P2s spectra for the zirconium ion before and after fluoride adsorption. Binding energy (BE). The full width at half maximum (FWHM). Gaussian:Lorentzian.
The zeta potential of the zirconium-doped algal biomass in the absence (0 ppm F−) and the presence of fluoride (100 ppm F−) are depicted in Figure 7. The point of zero charge values obtained by linear interpolation of the experimental data is 4.5 and 6 for pure zirconium-doped Zeta potential values of the zirconium-doped (a) 
Equilibrium studies
The equilibrium concentrations in the solution (
The data so obtained were used to determine which model fits best (Tan et al., 2008). The four adsorption models employed in the present study were: the Freundlich, Langmuir, Dubinin–Radushkevich (D-R) and Temkin adsorption.
The Freundlich model is given by equation (4) and represents adsorption on a heterogeneous surface (Tan et al., 2008).
Adsorption isotherm plots of (a) Langmuir, (b) Freundlich, (c) D–R and (d) Temkin models for zirconium-doped Adsorption isotherm constants of Freundlich, Langmuir, D–R and Temkin models for zirconium-doped

The Langmuir adsorption isotherm model is given by equation (5) (Tan et al., 2008).
The adsorbent uptake capacity (
Based on the correlation coefficient values, the equilibrium data of fluoride ion biosorption were best described by the Langmuir equation for zirconium-doped
D-R isotherm was applied to understand the adsorption mechanism in terms of Gaussian energy distribution onto a heterogeneous surface. The linearized D-R isotherm equation is represented by equation (7) (Tan et al., 2008).
The Temkin isotherm was applied to take into account the adsorbent–adsorbate interactions. The model ignores the extremely low and significant value of concentrations and assumes that the heat of adsorption (a function of temperature) of all molecules in the layer would decrease linearly rather than logarithmically with coverage (Tan et al., 2008).
The model is represented by the equation (9).
Adsorption kinetics
The mechanism of adsorption can be explained using kinetic constants obtained from different kinetic models. To investigate the sorption mechanism of the algal biosorbents, the first-order, pseudo-second-order, Elovich and intra-particle diffusion kinetic models were selected.
The kinetics of fluoride adsorption was analysed by the first-order rate expression, which is used when the liquid and solid phases attain equilibrium conditions (Chhipa et al., 2013). The first order rate equation is given by equation (10).
The pseudo-second-order kinetic model considers that there is chemical bond formation between the ion and the adsorbent surface. The linear form of the model is given by equation (11).
Elovich equation considers the sorbent surface to be energetically heterogeneous and is primarily used for chemisorption process (Chhipa et al., 2013). The linearized equation is given by equation (12).
The batch studies involve vigorous shaking, which may cause bulk transport of fluoride ions into the pores of the biosorbent as well as biosorption at the outer surface of the sorbent. There is a possibility of either film diffusion or intra-particle diffusion to be controlling the adsorption rate. Hence, the flow of fluoride ions into the pores of the biosorbent was tested by the intra-particle diffusion or the Weber and Morris model. The intra-particle diffusion model by Weber and Morris is represented by equation (13) (Chhipa et al., 2013).
The pseudo-second-order kinetic model has shown highest regression coefficients for both zirconium-doped Adsorption kinetic models (a) first-order (b) pseudo-second-order (c) Elovich (d) intra-particle diffusion models for the sorption of fluoride by zirconium-doped Adsorption kinetic parameters for the adsorption of fluoride by zirconium-doped 
Statistical analysis
The central composite design matrix was used to determine the effect of variables: pH (2.0–8.0), initial fluoride concentration (10–100 mg/L), adsorbent dosage (1–10 g/L) and contact time (30–240 min) on the response (
Validation of response surface models and statistical analysis
Second-order polynomial equations drew the relationship between independent variables and response. The regression equation coefficients were evaluated and fitted to a second-order polynomial equation for fluoride removal using zirconium-doped algal strains
The empirical relationship between the response and process variables was obtained using one of the statistical testing tools called
Optimization of variables for the removal of fluoride ions
The effect of independent variables: pH, initial fluoride concentration, adsorbent dosage and contact time, on the biosorption of fluoride ions, were evaluated using the empirical equation obtained by RSM (equations (14) and (15)).
The pH was one of the significant factor (
Mechanism of fluoride removal
Till now, only a few published articles have been successfully explained the sorption mechanism. The algal cell walls generally contain either polysaccharides or a variety of glycoproteins or both. The main reason behind the removal of fluoride was considered to be the protonation of hydrogen ions and the presence of amine groups on the algal biomass. Maximum fluoride removal was at acidic pH due to the protonation of the biosorbent and the speciation of fluoride ions (Guibal et al., 2001). The important mechanisms namely electrostatic attraction, ion exchange and adsorption are considered to influence the fluoride removal. The selective separation of anions is further increased by surface modification of the algal biomass using tetravalent zirconium ions. The zirconium ions bonds strongly with negatively charged groups present on the algal cell wall and also have high affinity towards fluoride ions. The sorption mechanism is effective even at low fluoride concentrations. The FTIR analysis results confirm that the amine groups present on the biomass are the primary adsorption sites. The adsorption mechanism of the biomass was also investigated using XPS and SEM-EDS. The XPS and EDS spectra established the interaction between fluoride and zirconium ions on the zirconium-impregnated biomass. From the initial investigations of XPS, it appears that fluoride ions bind with zirconium and phosphorus groups present on the surface of algal biomass.
Desorption studies
Maximum desorption of fluoride ion (93.4%) was observed during first adsorption–desorption cycle (Figure 10). To understand the reusability of the biosorbents, adsorption–desorption cycle of fluoride was repeated four times using the same biosorbent. It was noticed that the adsorption capacities of the biosorbents reduced by 6.55% during the first regeneration and after that decreases gradually.
Adsorption/desorption cycles for fluoride removal by zirconium-doped 
Effect of co-anions on the fluoride adsorption
The presence of co-ions can either increase or decrease the adsorption of anions (Errais, 2011). The adsorption of ions may sometimes not be affected by the addition of ions (Baghriche et al., 2008). Hence, the effect of the presence of co-ions such as sulphate, nitrate, chloride, bicarbonate, and phosphate on the fluoride adsorption phenomena was studied by considering mixtures of 10 and 100 ppm concentrations of salts and 10 ppm initial fluoride concentration. The results of the study of the adsorption in the presence of the co-anions are shown in Figure 11. It was observed that the fluoride adsorption was least affected by chloride and was most affected by the presence of bicarbonates in the solution mixture. This can be interpreted by the fact that the bicarbonate ions compete more with the active sites present on the biomass than the other anions and hence the number of active sites for adsorption gets reduced. A similar result was observed by Aravind and Elango (2006) and Karthikeyan et al. (2004). The chances of binding of the co-ions to the biomass when compared to fluoride depends on several factors such as the chemistry of the ions, pH of the solution, the nature of the binding sites, the amount of binding sites, the diversity of chemical species, fluoride ion concentration and the selectivity of the biomass to bind certain species.
Effect of co-ions on the fluoride adsorption.
Comparison with other adsorbents
Comparative evaluation of various adsorbents for fluoride removal.
Conclusion
The CCD of RSM was used to optimize the various factors including pH, initial fluoride concentration, adsorbent dosage and contact time, which influenced the adsorption of fluoride onto zirconium-doped
Footnotes
Acknowledgements
The authors are thankful to the Principal and Management of Siddaganga Institute of Technology, Tumkur, Karnataka, India for their constant support and encouragement. The authors would also like to thank Dr Madhu Chennabasappa, Department of Physics, Siddaganga Institute of Technology, Tumkur for his help in interpreting the XPS spectra.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research work is partially supported by R&D grant from The Institution of Engineers (India).
