Abstract
Adsorption process plays an important role in the removal of phosphorus from aqueous solutions. A laboratory experiment was conducted to investigate the adsorption characteristics of phosphorus onto natural zeolite and to find out the relative importance of some controllable treatments in phosphorus adsorption process using the Taguchi optimization methodology. Results showed that the adsorption of phosphorus in the presence of Fe3+ and Al3+ was higher than that in the absence of these two cations probably due to the adsorption of phosphorus-bearing anions on opposite charges of these cations. Also, increase in contact time tended phosphorus adsorption to be increased. The addition of base and acid treatments caused an increase and a decrease, respectively, on phosphorus adsorption. The order of effectiveness of treatments on the values of phosphorus adsorption was as follows: acid/base treatment >sorbent to sorbate ratio > modification with aluminium (Al)/iron (Fe) >contact time >phosphorus concentration. Phosphorus adsorption data well fitted to the Freundlich isotherm model. The pseudo-second order was the best model describing phosphorus adsorption kinetics. According to the results reported herein, it is assumed that the main mechanism controlling phosphorus adsorption onto natural zeolite is chemisorption.
Introduction
The discharge of excess phosphorus (P) into freshwater ecosystems can cause eutrophication which is a concern from environmental standpoint demonstrating the importance of P removal from wastewater via physical, chemical, and/or biological techniques (Lürling et al., 2014; Wang et al., 2016). Adsorption plays an important role in the removal of nutrients including P (Zhang et al., 2014). Due to simplicity, being cost-effective and environmentally friendly (Googerdchian et al., 2018), and relatively high efficient nature, this technique is one of the most prevalent physical techniques for the removal of a wide range of pollutants (Callery et al., 2016; Oguz, 2004).
Zeolites, naturally occurring hydrated alumino-silicates, are kinds of plentiful, high chemical and mechanical stability, high specific surface area, abundant reserve, low cost, and environmentally friendly materials which have been widely used as adsorbents for environmental cleanup (Šljivić et al., 2009; Wang and Peng, 2010; Zhan et al., 2017). Natural zeolites especially clinoptilolite types are known as an efficient ion exchanger and ion selective materials and substrates for fixing aqueous ionic species through their wide range of channels present in their crystalline structure and also their strong surface reactivity (Altin et al., 1998; Mahabadi et al., 2007). Guo et al. (2017) studied P adsorption by parent Al-modified eggshell and they observed that the value of P adsorption in Al-modified eggshell was significantly higher than parent eggshell. Similar results for Fe3+ and Na+ in P removal by ferric-modified zeolite reported by Gao et al. (2018). Presence of Ca2+ and Na+ enhanced P adsorption by zirconium-modified zeolite (Zhan et al., 2017).
Adsorption can be affected by different factors including solute concentration, sorbate to sorbent ratio, contact time, temperature, sorbent type and modification, ionic strength, and pH (Wang et al., 2016; Zhang et al., 2015; Zolfaghari et al., 2011). Adsorption characteristics of P can be investigated using linear and nonlinear kinetic and isotherms studies (Zhang et al., 2015). Commonly used Langmuir and Freundlich models are useful in describing adsorption isotherms and providing effective factors in predicting adsorption characteristics of nutrients (Cheng et al., 2014).
Wang et al. (2016) found that the Langmuir isotherm well predicts the patterns of P adsorption by fly ash, red mud, and ferric-alum water treatment residues. They reported that the pseudo-second order was the best model describing the kinetics of P adsorption onto different solid wastes.
To achieve maximum removal of pollutant, adsorption factors should be optimized. Therefore, using an experimental design based on statistical techniques of optimization has been found to be an effective method for analysis and evaluation of the best conditions of different investigations (Googerdchian et al., 2018). The Taguchi approach is a statistical experimental method designed based on orthogonal arrays which are much efficient in comparison with full factorial experiments since the number of experiments and hence time and cost of the experiment can be drastically reduced by using this method (Madan and Wasewar, 2017). It has been proved to be successful over the last 15 years for the improvement of product quality and process performance (Roy, 2002). Taguchi optimization methods have been used in aqueous solution to investigate removal of different pollutants such as removal of nickel (Ni) from waste water by calcined ouster shell powders (Yen and Li, 2015), adsorption of P onto Al/Ca-impregnated granular clay material (Yu et al., 2015), Cd adsorption by spent coffee (Yen and Lin, 2016), removal of Th (IV) from aqueous solution by deoiled karanja seed cake (Varala et al., 2016), removal of copper and Ni by growing Aspergillus sp. (Pundir et al., 2016), removal of benzeneacetic acid from aqueous solution using CaO2 nanoparticles (Madan and Wasewar, 2017) and lead and mercury by nanostructure and modified zinc-oxide (Zolfaghari et al., 2011).
Research on the P removal from aquatic systems using low-cost materials is of great interest for researchers. Limited research (if any) has been conducted on Taguchi optimization for P adsorption using modified zeolite. The main objectives of present study were to investigate the adsorption characteristics of P onto natural zeolite and to find out the relative importance of initial P concentration, sorbent to sorbate ratio, contact time, modification of sorbent and acid/base treatments on the retention of P onto natural zeolite using Taguchi optimization methodology.
Experimental analysis
Taguchi experimental method
According to the Taguchi method, experimental design is arranged based on the orthogonal arrays which consist of the controllable factors. Five controllable factors including sorbent to sorbate ratio (SR), P concentration (C), contact time (T), acid/base addition (AB) and sorbent modification using Al/Fe (AF) were considered in the present study (Table 1).
Controllable factors and associated levels.
The Taguchi experimental design (L27 orthogonal array) (Minitab, Inc., version 14, USA) was used to determine the optimal conditions for the adsorption of P from aqueous solution (Table 2). A total number of 35 = 243 experiments would require for a full factorial experiment. By using the Taguchi method, the number of experiments was reduced to 27 experiments. Therefore, only about 10% of the total number of samples was considered by the orthogonal array.
Design of experiments.
T: contact time; C: phosphorus concentrations; SR: sorbent to sorbate ratio; AB: acid/base addition; AF: sorbent modification using Al/Fe.
The signal to noise ratio (S:N) is needed for the evaluation of the experimental results and quality determination. Three types of S:N ratio analysis are embedded: (1) lower is better (LB), (2) nominal is best (NB), and (3) higher is better (HB). The S:N ratio of HB was considered since the target of this study was to maximize P removal efficiency (equation 1)
As the goal of this study was to remove P by zeolite, the quality characteristic selected was HB of P removal defined by equation (2).
ANOVA statistical method is also used to have a better understanding of how the observed results are reliable (Googerdchian et al., 2018). ANOVA performs to investigate the influence of each controllable factor on the removal of P and to determine percentage contribution of each factor (Madan and Wasewar, 2017). One-way ANOVA using Duncan’s multiple range test was employed to determine the significant differences (α = 0.05) in relevant parameters during modification. The result rechecked by ANOVA variance statistical approach using SPSS 17.0. The percentage contribution of each factor, RF, was calculated using equation (3).
Three levels of initial P concentrations (60, 75, and 90 mg L−1) as potassium dihydrogen phosphate in 0.01 M calcium chloride solution were used. Zeolite used in the experiment belongs to a natural zeolite mine located in Semnan province, Iran (Table 3). Three levels of SR including 1:100, 1:150, and 1:250 considered in the experiment. Three levels of acid/base treatments including 0 M hydrochloric acid, 0.5 M hydrochloric acid, or 0.5 M sodium hydroxide (1:10 zeolite mass) were also considered. Blank (nontreated zeolite) aluminum hydroxide (Merck, product number 1091) and iron metal (BDH, product number 28602) used as Al/Fe treatments for zeolite modification (1:5 zeolite mass) because of P-bearing anions have negative charges and could be easily adsorbed onto positively charged cations. Zeolite treated by Al/Fe or acid/base for 24 hours in an oven at 65°C. Acid/base and Al/Fe treatments add to the samples with abovementioned concentrations and shaken with an end over end shaker for periods of 5, 10 and 20 hours as T treatments. The selection of these periods was due to the fact that adsorption of P occurred mostly through the first shaking periods up to 20 hours (Figure 2). Similarly, Wang et al. (2016) observed that a considerable increase in P removal usually occurred during first shaking periods up to about 1 day. After each shaking period, samples centrifuge for 15 minutes and filtrated with Whatman no. 42 filter paper. The content of P was measured according to the Olsen method (1954) by spectrophotometer (Shimadzu AA-670G).
The values of zeolite characteristics.
The values of phosphorus removal and S:N ratio for different treatments.
Kinetics experiment
To investigate the kinetics of P adsorption, zeolite equilibrates with five concentrations of P (30, 45, 60, 75, and 90 mgL−1) for periods of 5, 10, 20, 30, and 40 hours. This experiment performs with three SR ratios of 1:100, 1:150, and 1:250. The samples were centrifuged immediately and then filtered to separate the liquid and solid phase. The amount of P in supernatant was measured by using spectrophotometer.
Kinetics equations used in the experiment were as follows (equations (4) to (7))
The goodness of fits was determined by comparing the values of coefficients of determination (R2), and standard errors of estimate (SE). Kinetic models with relatively high R2 and low SE values were considered as the best-fitted models.
Adsorption isotherms
To investigate the influence of different SR on P adsorption, 0.12, 0.2, and 0.3 g of zeolite equilibrated with 30 mL solution containing 30, 45, 60, 75, or 90 mg P L−1 as potassium dihydrogen phosphate in 0.01 M CaCl2 on an end over end shaker for 40 hours. The time 40 hours considered as equilibrium time in present study (Figure 1). Similarly, in previous studies, it has been observed that after 35 hours, the absorption of P is steady and without changes (Wang et al., 2016). Samples centrifuged and filtered and the amount of P in supernatants determined. The amount of P adsorption at equilibrium was calculated by equation (8)

Time effect in amount of phosphorus (P) adsorbed by zeolite in different P concentrations.
The values of n and kf, the Freundlich coefficients, calculated from the slope and intercept of the log qe versus log Ce in the linear form of Freundlich isotherm. It has been shown that the adsorption mechanism is chemical, physical, and linear for n being less than, greater than and equal to unity, respectively (Bansal et al., 2005).
Result and discussion
Taguchi analysis
Taguchi method used to find the ideal conditions for removing P from aqueous solution. Different S:N ratios for the evaluation of the significance of treatments were shown in Figure 2. As seen in Figure 2(a), the values of S:N ratios increased by increasing in contact times as 20 hours treatment was the maximum S:N ratio resulted in higher P removal by zeolite. A similar result of P adsorption in zeolite was reported by Gan et al. (2016). They observed that greatest P removal from aqueous solution by modified zeolite at 20–25 hours. According to Figure 2(b), the value of S:N ratio increased as initial P concentration decreased. Yu et al. (2015) and Zhang et al. (2014) observed similar results for P removal percentage in aqueous solutions. Decreasing initial Pb concentration by increasing the values of S:N ratio reported by Googerdchian et al. (2018).

Estimated individual factors effect on percentage of S:N by zeolite from aqueous solution.
The effect of SR on P adsorption is presented in Figure 2(c). As SR increased, S:N ratios decreased probably due to the reduction in adsorption sites. Similar results have been reported by Yen and Lin (2016). Maximum dosage for P adsorption in the bentonite-alum system was 0.1 g or SR (1:300) (Mahadevan et al., 2018) which was relatively close to the value of SR in the present study (1:250).
Increase in Cd and Hg adsorptions as a result of an increase in pH value has been reported by Zolfaghari et al. (2011) which is in agreement with results of a present study which demonstrated the significant role of NaOH in increasing P adsorption (Figure 1(d)). According to the results of Oliveira et al. (2015), zeolite was an efficient and promising adsorbent for P removal. They observed that maximum P removal occurred in higher pH. Al-Zboon (2017) showed that high P removal efficiency by activated carbon–silica nanoparticles composite, kaolin, and olive cake at high pH conditions. Gao et al. (2018) investigated phosphate removal by ferric-modified zeolite and reported modification treatment media such as Na salt and HCl increased phosphate sorption removal in order Na + Fe > Fe > HCl + Fe.
Effect of sorbent modification using Al/Fe treatments was presented in Figure 2(e). The sorbent modification of zeolite using Fe was an effective operation which considerably increased the value of P adsorption by sorbent from aqueous solution. It appears that Fe is also effective in the stabilization of cadmium (Saffari et al., 2015). Oliveira et al. (2015) observed high P removal ability in Fe-modified zeolite followed by Al-modified zeolite. Shanableh et al. (2016) studied P adsorption using Al/Fe modified bentonite. They observed that using Fe-modified bentonite increased both adsorption capacity and adsorption rate as compared to using Al-modified bentonite.
The highest significant factor influencing P adsorption onto zeolite was AB and the order of treatments importance was as follows AB > SR > AF > T > C. Yu et al. (2015) observed that the type of adsorbent was the most significant factor affected P removal. The results of ANOVA showed the same order for the influence of different treatments on P adsorption. Sum of square (SSF) and percentage contribution (PF) of each controllable factor are given in Table 5. The highest SSF and PF were related to AB treatment, followed by SR, AF, T, and C, respectively. The adsorption of P in the presence of Fe3+ and Al3+ was higher than that in the absence of these two cations, probably due to the fact that P-bearing anions have negative charges and could be easily adsorbed onto positively charged cations, resulting in increasing percentage of P removal from the solution (Zhou et al., 2011). Other researchers have reported similar results in their research on other adsorbents along with Fe and Al oxides (Liu et al., 2007; Oliveira et al., 2015; Pengthamkeerati et al., 2008; Xiong and Peng, 2008).
Response table of S:N ratio and contribution of each controllable factor.
DOF: degrees of freedom; SSF: sum of square of factor; PF: percentage contribution of factor.
Kinetics of P adsorption
The patterns of P adsorption at different periods evaluate by different models including pseudo-second-order, pseudo-first-order, elovich, and power function models (Table 6). The values of R2 and SE for different models are given in Table 6. The kinetics of P adsorption from aqueous solution by zeolite well fitted to the models in order of pseudo-second order, pseudo-first order, power function, and elovich (Table 6). The fit of experimental data to elovich equation was not satisfied due to high values of SE for this model (Table 6). The values of P adsorption decreased as the values of SR increased (Figure 3). The conformity of the experimental data to the pseudo-second-order model indicates that the P adsorption process is chemically controlled (Ofomaja et al., 2010). Similarly, Aljbour et al. (2017) observed that phosphate removal from aqueous solutions by natural Jordanian zeolitic tuff well conformed to pseudo-second-order model. As shaking periods increased, P sorption on zeolite surfaces increased which is in agreement with other studies (Wang et al., 2016; Zhang et al., 2016). Lin et al. (2018) reported P adsorption onto chitosan fitted to pseudo-second-order model indicating the occurrence of chemical adsorption. Pseudo-second-order parameters (qe and k) were calculated and shown in Table 7. The comparison of qe and initial concentrations of P indicated that qe increased as initial P concentration increased. The value of qe derived from pseudo-second order was higher with 0.2 g zeolite (1:150 SR) which is in accordance with the results of isotherm study. Zhang et al. (2014) reported relatively similar results and observed that the optimum dosage for P removal by zeolite was 0.25 (g L−1).
The values of R2 and SE for different kinetic models.
C: concentration; SE: standard error; SR: sorbent to sorbate ratio.

Phosphorus (P) adsorption at different sorbent to sorbate ratio using the nonlinear pseudo-second-order kinetic model, the P concentration in solution was 90 (mg L−1).
Mean R2 and kinetic parameters for pseudo-second-order model at different concentration (C) and sorbent to sorbate ratio (SR).
Adsorption isotherms
The P adsorption data were fitted to the Langmuir and Freundlich adsorption isotherm models (Table 8). Both the Freundlich and Langmuir models well described the patterns of P adsorption onto zeolite at different SRs as evidenced by their relatively high R2 and low SE values. Langmuir showed lower SE values and also lower R2 values as compared to the Freundlich. According to Kim et al. (2018) nf values between 1 and 10 showed that the adsorption was effective. They believed that as kf value increased the affinity of sorbate toward the sorbent increased. The value of nf derived from the Freundlich isotherm was in the following order: D2 (1:150) > D3 (1:200) > D1(1:100), indicating that 0.2 g zeolite had higher P adsorption efficiency (P adsorption per area unit). The largest kf were in D1 (1:100) zeolite which was attributed to the presence of the highest sorption site under this condition. Mahadevan et al. (2018) reported the values of nf ranged from 0.7 to 2.5 for P removal by the bentonite-alum system. Kim et al. (2018) reported the same results for P adsorption by fly ash. According to Novais et al. (2018), the values of nf ranged from 1 to 3 for P adsorption onto different types of biochar. A similar relationship was observed in Langmuir model between P adsorption capacity (bL) and zeolite dose. The zeolite dose of 0.2–0.3 (g ml−1) was the best dose for Cd adsorption as reported by Selim et al. (2017). The values of maximum adsorption capacity criterion (bL) ranged from about 40 to 50 mg g−1. Similar ranges reported by Abdelhay et al. (2017) for the adsorption of phosphate onto natural reed.
Langmuir and Freundlich isotherm parameters for phosphorus adsorption.
SE: standard error; SR: sorbent to sorbate ratio.
Conclusion
Results of Taguchi method shown that acid/base treatment was the most effective treatment influencing the values of P adsorption by zeolite from aqueous solution. The order of treatments importance was as acid-base addition >sorbent to sorbate ratio >sorbent modification using Al/Fe >contact time > phosphorus concentration. The optimization conditions for P removal in the present study are contact time (20 hours), P initial concentration (60 mg L−1), sorbent to sorbate ratio (1:250), and modification of zeolite with NaOH and Fe. The Freundlich model was better in describing the P adsorption on different sorbent to sorbate ratios as compared to the Langmuir. The kinetics of P adsorption from aqueous solution by zeolite well fitted to the pseudo-second-order model. Based on the results of kinetic studies, it appears that adsorption of P by zeolite is chemically controlled. Results reported herein demonstrated that the modification of zeolite with Fe along with base addition could be promising for efficient P removal from aqueous solutions.
Footnotes
Acknowledgements
Authors would like to appreciate Shiraz University for providing research facilities.
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) received no financial support for the research, authorship, and/or publication of this article.
