Quantification and isotherm modelling of competitive phosphate and silicate adsorption onto micro-sized granular ferric hydroxide

Adsorption onto ferric hydroxide is a known method to reach very low residual phosphate concentrations. Silicate is omnipresent in surface and industrial waters and reduces the adsorption capacity of ferric hydroxides. The present article focusses on the influences of silicate concentration and contact time on the adsorption of phosphate to a micro-sized iron hydroxide adsorbent (μGFH) and fits adsorption data to multi-component adsorption isotherms. In Berlin drinking water (DOC of approx. 4 mg L−1) at pH 7.0, loadings of 24 mg g−1 P (with 3 mg L−1 initial PO43−–P) and 17 mg L−1 Si (with 9 mg L−1 initial Si) were reached. In deionized water, phosphate shows a high percentage of reversible bonds to μGFH while silicate adsorption is not reversible probably due to polymerization. Depending on the initial silicate concentration, phosphate loadings are reduced by 27, 33 and 47% (for equilibrium concentrations of 1.5 mg L−1) for 9, 14 and 22 mg L−1 Si respectively. Out of eight tested multi-component adsorption models, the Extended Freundlich Model Isotherm (EFMI) describes the simultaneous adsorption of phosphate and silicate best. Thus, providing the means to predict and control phosphate removal. Longer contact times of the adsorbent with silicate prior to addition of phosphate reduce phosphate adsorption significantly. Compared to 7 days of contact with silicate (c0 = 10 mg L−1) prior to phosphate (c0 = 3 mg L−1) addition, 28 and 56 days reduce the μGFH capacity for phosphate by 21 and 43%, respectively.

Quantification and isotherm modelling of competitive phosphate and silicate adsorption onto micro-sized granular ferric hydroxide †

Introduction
Phosphate loads in up to 20% of European surface waters have to be drastically reduced to reach the good ecological status as demanded by the European water framework directive. 1 Critical phosphate concentrations for slowly owing or dammed waters are set as 0.1 mg L À1 total phosphorus (TP). 2 In Germany, the majority of surface waters exceed this value, especially shallow and polymictic lakes. Depending on the lake type, even lower values of 0.02-0.06 mg L À1 have been suggested in order to limit algae growth and to reach a good status. 3 Conventional techniques for P removal include precipitation with di-or trivalent metal ions and biological removal. Low residual phosphate concentrations are achieved with high doses of occulants. Besides high costs, the advanced treatment leads to an increased input of salts, increased sludge volumes and the need for large sedimentation tanks. A more simple and space saving alternative is the use of adsorption as it is a suitable method to reach low P concentrations. 4 Many studies on the adsorption of phosphate onto iron oxides have been conducted. [5][6][7][8] Adsorption of phosphate as a polyprotic anion takes place over a wide pH range and has its maximum at values below the pH PZC of the adsorbent.
Adsorbents are not specic for one target substance but accumulate a mixture of competing water constituents. Competitive adsorption of phosphate onto iron hydroxides was widely studied. 6,[8][9][10][11] While sulfate and chloride were shown to have no effects 12 and bivalent cations (e.g. Ca 2+ , Mg 2+ ) can have benecial surface complexation effects on phosphate, bicarbonate and silicate proved to have adverse impacts. 13 Silicate is released during wreathing of Si-containing minerals and occurs in natural waters in concentrations of 3 to 30 mg L À1 , mainly as Si(OH) 4 . 14 The sorption of silicate onto iron hydroxides receives ongoing attention. 9,[13][14][15][16][17][18][19] The adsorption is pH-dependent with increasing adsorption until the silicate species changes from H 4 SiO 0 4 to H 3 SiO 4 À at pK a 9.82. 13 The adsorption is accompanied by a net release of protons, which results in a decrease of the iso-electrical point of the adsorbent surface. At the pH PZC of the adsorbent the building of two different surface complexes is proposed by Hiemstra. 13 A binuclear bidentate complex with exchange of two ligands resulting in^(FeO) 2 Si(OH) 2 and a mononuclear monodentate complex^(FeOHFeO)Si(OH) 3 . At low concentrations and low pH values silicate is mostly present as monomers, whereas the amount of oligomeric and polymeric silicates increases steeply at circum-neutral pH values. 13 Polymerization was observed at molar ratios of Si/Fe below 0.1. 14,15,20 In competitive adsorption experiments phosphate was less affected by silicate than oppositely. 13 While at low pH-values the effect of phosphate on silicate was very large, it was less for higher pH-values. For the competitive adsorption of silicate and arsenate no inuence was found for pH values below 9, but silicate was able to replace arsenate when adsorbing onto ferrihydrite above pH 9. 20 However, polymerization of silicate was prevented when arsenate was added rst. 14 They conclude that polymerization is a surface controlled process. Christl et al. 18 suggest that a long pre-equilibration of the adsorbent with silicate might alter the surface conguration of the adsorbent, favouring the formation of oligomers and polymers and thus reducing reactivity. No binary adsorption isotherms to t this data have been developed yet to our knowledge.
Given the high relevance of phosphate for surface water quality, the adverse impacts of the ubiquitously present adsorptive competitor silicate requires further elucidation. The current study highlights the effects of variable silicate concentrations, preloading times, and iron-complexation on phosphate adsorption onto a micro-sized iron hydroxide (mGFH, characterized in a previous study 21 ). In addition, various bisolute models are tested for describing the observed competitive adsorption. Further, we also tested the potential reversibility of silicate and phosphate adsorption under practically relevant conditions. The provided results improve the understanding of competitive phosphate-silicate adsorption on iron hydroxide surfaces, thus providing better means for process control and safety.

Adsorbent
Micro-sized ferric hydroxide was obtained from GEH Wasserchemie (Osnabrück, Germany). Granular ferric hydroxide (GFH) is produced by precipitation of Fe(OH) 3 from a FeCl 3 solution and subsequent conditioning of the resulting sludge to obtain stable granules of up to 3 mm grain size. 22 In the current study, we used only particles <0.3 mm, potentially increasing adsorption kinetics and overall loadings. 23 Conventional GFH consists to 50-70% of akaganeite (b-FeOOH), ferrihydrite, and other iron oxides. 24 Saha et al. 25 reported a point of zero charge (pH PZC ) of 7.5 for GFH. The supplied mGFH material had a water content of approx. 53% and a specic surface area of approx. 300 m 2 g À1 dry weight. Sieve analyses showed a contribution of 30, 25, and 45% weight of the fractions < 63 mm, 63-120 mm and 120-300 mm, respectively, and particle counting with a particle analyser (PAMAS SVSS) revealed that 99% of the particles were smaller than 63 mm. 21 The mGFH was wet-sieved to the desired grain fractions and then air-dried at room temperature. The residual water content varied between 10 and 20%. Previous experiments showed no alteration of the material due to drying. A detailed characterization of the adsorbent was provided by Hilbrandt et al. 21

Experimental set-up
Batch experiments were conducted at room temperature (20 AE 2 C) to obtain adsorption isotherms. All chemicals and reagents used were of analytical grade or higher. NaCl was added to deionized water (DI) to set an ionic strength of 10 mmol L À1 . Initial phosphate and silicate concentrations (added as KH 2 PO 4 and Na 2 SiO 3 $5H 2 O) were adjusted to 4 mg L À1 P and 10 mg L À1 Si respectively (unless indicated otherwise). 2 mmol L À1 MES (2-(N-morpholino)ethanesulfonic acid, pK a ¼ 6.1), BES (N,N-bis(2hydroxyethyl)-2-aminoethanesulfonic acid, pK a ¼ 7.1) or TAPS (N-tris(hydroxymethyl)methyl-3-aminopropanesulfonic acid, pK a ¼ 8.4) were used to buffer the pH value, adjusted with NaOH. All buffer substances were purchased from Sigma Aldrich (Germany). No signicant adsorption of the buffer substances onto mGFH was found in preliminary tests and thus no competition with the target substances was expected. The isotherms were obtained by varying the adsorbent doses between 40 and 500 mg L À1 . Dry mGFH was added to 100 mL or 150 mL batches of test solutions containing the target substances and shaken at room temperature at 220 rpm. Aer 72 h, the adsorbent was separated using membrane ltration (0.45 mm, cellulose nitrate), as kinetic experiments showed no further adsorption aer that time period. 26 For desorption experiments the suspensions were centrifuged aer 72 h contact time (as described above, without adsorbent removal by ltration), following a method introduced by Aschermann et al. 27 95 mL of the solution were removed aer centrifugation and the same volume of desorption solution was added. The desorption solution corresponded to the adsorption solution without the target substances. The desorption batches were shaken for another 72 h prior to adsorbent removal by ltration and subsequent analyses.
For the inuence of silicate in simultaneous adsorption, mGFH was added to the above described test solution containing 3 mg L À1 PO 4 -P and varying concentrations of silicate (6, 10 and 18 mg L À1 ). In addition to competitive phosphate-silicate tests in DI-water, drinking water tests were conducted. The concentrations of 6, 10, 18 mg L À1 Si were added to drinking water which itself had a Si concentration of 6 mg L À1 , resulting in nal Si concentrations of 12, 16, 24 mg L À1 . Sequential adsorption experiments were carried out as described for the single solute tests but with the addition of the second substance from a highly concentrated stock solution aer a contact time of 72 h. The equilibrium concentrations of both substances were measured aer the cumulative adsorption time of 144 h. Further experiments were carried out with 7, 28 and 56 d contact time between silicate and mGFH with subsequent addition of phosphate and additional 7 d of contact time in the binary solution.

Analytical methods
mGFH doses in the batch experiments were gravimetrically controlled by membrane ltration through pre-washed 0.45 mm cellulose nitrate lters, which were dried for 24 h at 105 C and weighed before and aer ltration. Orthophosphate was quantied via ow injection analysis according to ISO. 28 Silicate concentrations were determined photometrically (Lambda 12, PerkinElmer) with test kits (Spectroquant, Merck) using 10 mm glass cuvettes. The uncertainty of the measurement is shown with a 95% condence interval expressing the summed up errors of the experimental and analytical processes.  (1)) with the equilibrium loading q i (mg g À1 ), the Freundlich coefficient K F,i (L n (g mg nÀ1 ) À1 ), the equilibrium concentration c i (mg L À1 ) and the Freundlich exponent n i (À), of Langmuir (eqn (2)) with the maximum loading q m,i (mg g À1 ) and the Langmuir coefficient

Multi-component isotherms.
For the description of adsorption in multi-solute systems, multi-component isotherm equations were developed by extending the above described single component isotherms. 29 The following equations were used in this work: (a) Non-modied Langmuir multi-component isotherm (NLMI) with the Langmuir parameters q m,i and K L,i derived from the single-component isotherms 30 (b) Modied Langmuir multi-component isotherm (MLMI) with the Langmuir parameters q m,i and K L,i derived from the single-component isotherms and an additional interaction term h i , 31 number of tting parameters: 2 (n i , n j ) (c) Extended Langmuir multi-component isotherm (ELMI) with a total maximum loading of both components q max (mg g À1 ) and Langmuir coefficients K i and K j in the binary solution, 32 number of tting parameters: 3 (q max , (d) Extended Freundlich multi-component isotherm (EFMI) with the Freundlich parameters K F and n from the singlecomponent isotherm and the binary coefficients x i , y i and z i , 31,33 number of tting parameters: 6 (x 1 , y 1 , z 1 , x 2 , y 2 , z 2 ) (e) Sheindorf-Rebuhn-Sheintuch equation (SRS) with the Freundlich coefficient from the single-solute isotherm K F and the interaction factor a ij , 34 number of tting parameters: 2 (a ij , a ji ) (g) Modied Redlich-Peterson multi-component isotherm (MRPMI) with the single-component isotherm parameters and an additional interaction factor h i , 35 number of tting parameters: 2 (n i , n j ) (h) Ideal adsorbed solution theory (IAST) with the mole fraction z i , the spreading pressure 4, and the Freundlich parameters K F,i and n i and the total loading q T , 29 number of tting parameters: 2 (4, q T ) 2.4.3. Parameter tting and determination. All models described above use isotherm parameters derived from single-solute ttings. The parameter tting was thus applied to the additional parameters for binary systems. The parameters were tted by minimizing Marquardt's percent standard deviation (MPSD). 36 The equation incorporates the degree of freedom into the geometric mean error distribution and was used by Srivastava et al. 31 for tting of binary adsorption isotherms of cadmium and nickel onto bagasse y ash. Through the incorporation of the degree of freedom a better comparison between models with different numbers of tting parameters can be assured. n m is the number of experimental data points and n p the number of parameters in the isotherm equation.
MPSD ¼ 100 ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 1 n m À n p X n i¼1 q e;i;exp À q e;i;cal q e;i;exp To avoid nding results related to a local minimum of the minimization criteria a matrix of starting values containing at least 500 different combinations of the parameters was used. Thus, it was made sure to nd the global minimum and the best tting parameters.

Effect of pH
The single-solute adsorption isotherms for phosphate and silicate are dependent on the pH value and the water composition as shown in Fig. 1. Fitting parameters were obtained through linearization of the single-solute models and are listed in Table  1. The mGFH surface sites with a pH PZC of 7.5 react as acid or base dependent on the pH. 25 At pH values below the pH PZC the mGFH surface is positively charged and therefore the adsorption of anions is electrostatically favored. In the pH range of 6.0 to 8.0 phosphate is present as HPO 4 À and H 2 PO 4 2À (pK a ¼ 7.2) and thus adsorption capacity is higher for lower pH values. As the pH increases, more surface groups of the adsorbent are uncharged and negatively charged. Phosphate becomes twofold negatively charged. Thus the adsorption capacity decreases for pH values above the pH PZC . Several studies revealed that phosphate forms strong inner-sphere complexes such as bidentate, binuclear or monodentate surface complexes with ferric (hydr)oxide. 6,13,37 Silicate (pK a ¼ 9.8) loadings increase with increasing pH, as shown in Fig. 1. Also, the percentage of polymeric (mostly dimeric) silicate rises with increasing pH. Davis et al. 38 suggested a direct sorption mechanism of monomeric and dimeric silicate to iron oxide surfaces, which leaves the surface with a net positive charge up to loadings of 1 mol mol À1 of adsorbent. The resulting negative surface charge acts electrostatically on the remaining free surface groups by increasing the fraction of protonated groups. 13 This explains the high loadings of 28 mg g À1 at initial concentrations of 10 mg L À1 silicate.

Effect of water composition
P loadings of approx. 17 mg g À1 were reached at an equilibrium concentration of 3 mg L À1 P at pH 7.0 with an ionic strength of 10 mmol L À1 . The adsorption in drinking water shows higher capacities with approx. 23 mg g À1 P. Several studies have shown that calcium, which occurs in concentrations up to 100 mg L À1 in Berlin drinking water, enhances the adsorption of phosphate 7,10 due to a suggested multi-layer coverage or the formation of ternary surface complexes. Furthermore, the possibility of surface precipitation of calcium phosphate was discussed. 10,39 Contrasting effects can be seen in the adsorption of silicate. In drinking water mGFH shows much lower silicate loadings than in DI water (17 mg g À1 compared to 27 mg g À1 at c e ¼ 9 mg L À1 ). The addition of Ca 2+ to DI water with NaCl led to a rise in adsorption of approx. 15%. The adsorption of calcium leads to a more positive surface charge of the adsorbent, resulting in higher adsorption capacities of silicate. Thus, the reduction in DW is linked to the presence of other competing anions. Bicarbonate was identied as main competitor for chromate adsorption onto mGFH 40 and might also adversely impact silicate adsorption.

Reversibility of adsorption
The desorption isotherms of phosphate and silicate in DI + NaCl water are shown in Fig. 2, together with Freundlich isotherm ttings. In the case of phosphate (Fig. 2a), adsorption is partly reversible as indicated by partially overlapping ad-and desorption isotherms. An irreversible portion of adsorbed phosphate can also be observed, shown by the desorption isotherm being located above the adsorption isotherm around liquid-phase concentrations of >0.1 mg L À1 . Similar results with 16% of irreversible bonds and slow re-diffusion out of micropores were described by Cornell and Schwertmann. 37 Atkinson et al. 41 explained the irreversibility by the formation of inert binuclear surface complexes. The adsorption isotherm of silicate is only partially shown in Fig. 2b. Independent of the loading reached during the adsorption phase, silicate desorbs down to an equilibrium concentration of approx. 2 mg L À1 in the desorption solution. Thus, high loadings remain in the adsorbent. High surface concentrations and contact times of silicate lead to the formation of oligomers and polymers and irreversible bonds to the hydroxide surfaces. This process is favoured by the neutral pH range. 13 Padungthon and Sen-Gupta 19 regenerated a silicate loaded xed-bed adsorber and state that desorption of silicate is more difficult due to polymerization. The mole ratios applied in the present study (Si/Fe ¼ 0.3 to 0.6) exceeded by far the ratio above which polymerization occurs (Si/Fe > 0.1). 20 Thus, the irreversible formation of Si-O-Si bonds is likely.

Multi-component adsorption modelling
Silicate competes with phosphate for adsorption sites and reduces P-loadings signicantly (Fig. 3). Phosphate loadings are reduced by 27, 33 and 47% (c eq ¼ 1.5 mg L À1 ) for 9, 14 and 22 mg L À1 Si respectively. The difference in loadings between the single-solute and the multi-component systems increases with decreasing adsorbent doses (and increasing residual Pconcentrations). Sorption competition is strongly inuenced by kinetics and the sorption kinetics of phosphate are considerably faster than the kinetics of silicate. 40 The fraction of oligomeric silicates on the surface of the adsorbent steadily increases with time and loading. 18 Especially at high loadings, silicate strongly competes with phosphate for limited adsorption sites. Silicate adsorption is also strongly affected by the presence of phosphate (Fig. 3). Through the addition of 5 mg L À1 P the slope of the resulting isotherm changes from positive to negative and loadings are strongly reduced.
The adsorption data of silicate and phosphate in simultaneous adsorption was tted using the different multi-component isotherm models (eqn (4)- (13)). Comparisons between experimentally obtained and calculated loadings are shown in the parity plots (Fig. 4). The closer the data points are to the bisecting line (45 line), the better the respective model ts the experimental data. An exemplary application of the described models for the binary adsorption of phosphate and silicate with initial concentrations of 5 and 15 mg L À1 respectively is shown in Fig. S2 † as 2D isotherms. The associated MPSD values for all tested models are given in Table 2 and the isotherms are shown in Fig. S3 and S4. † Fig. 1 Adsorption isotherms dependent on the pH (upper row) and on the water matrix (lower row) for phosphate and for silicate; errors indicate the 95% confidence interval; Freundlich isotherms as continuous and Langmuir as dotted lines fitted, with constants for pH 7.0 and DI provided in Table 1. Redlich-Peterson isotherm fittings are found in Fig. S1. † Table 1 Isotherm parameters for the removal of phosphate and silicate by mGFH (pH 7.0, DI). A complete list can be found in Table S1 Langmuir Freundlich Redlich-Peterson The non-modied Langmuir model (NMLM) creates a poor t with a MPSD value of 80.6 as it presumes the adsorbates to adsorb independently without interaction. The introduction of the interaction term h results in a better t of the MLMI (MPSD ¼ 61.0) with a deviance of phosphate loadings below 30%, but overestimated silicate loadings. Similar ttings were achieved using the extended Langmuir model (ELMI). The MLMI uses only two tting parameters additionally to the singlecomponent Langmuir parameters whereas three parameters are used in the ELMI. As both models show similar ttings to the experimental data but the MLMI uses fewer tting parameters it was chosen as the best Langmuir-type model.
However, the best t with a MPSD value of 16.6 was obtained using the extended Freundlich equation assuming multilayer coverage (EFMI). This is expected as mGFH has heterogeneous surface sites and single solute isotherms were also well represented using the Freundlich model. Also, the EFMI used six tting parameters additionally to the single-component Freundlich parameters. Thus, the interactive effects between phosphate and silicate are taken into account by the modication of the Freundlich equation. The resulting binary isotherms for phosphate and silicate are shown in Fig. 5. Phosphate loadings show a sharp increase for low phosphate or silicate concentrations, but only low additional adverse inuences of silicate when its concentrations exceed 8 mg L À1 . According to these results, even low silicate concentrations have strong negative effects on phosphate adsorption. Thus, an impact of silicate can be expected for most practical situations and water matrices. Silicate loadings are greatly reduced by the presence of phosphate, even at the lowest tested concentrations of 0.5 mg L À1 . For P-concentrations of 1-5 mg L À1 , the reduction of the loadings with silicate are almost constant.
The SRS model as another Freundlich derived model, but with only two tting parameters shows a very good t for silicate  This journal is © The Royal Society of Chemistry 2019 especially for higher loadings but an overestimation of up to 30% for phosphate loadings.
Both multi-component Redlich-Peterson isotherms do not reproduce the experimental values satisfactorily with errors exceeding 30% considerably. The modication with the interaction term decreases the MPSD value from 107 to 55.7.
The IAST as the only model based on thermodynamic considerations overestimates phosphate loadings by more than 30%. The poor t (MPSD of 55.1) is explained by the heterogeneity of adsorption sites of the iron hydroxide and the interaction of the adsorbates. Both phenomena are not considered in the model calculations.
Thus, the adsorption of phosphate in competition with silicate onto mGFH can be predicted using multi-component isotherms. Freundlich derived models were found to describe the system best with MPSD values of 38.4 and 16.6 for the SRS and the EFMI. The number of tting parameters used has to be taken into consideration.
In systems with dened matrices, resulting loadings of both adsorbents can be calculated and predicted with the described models. For drinking water (loadings shown in Fig. 1) the found correlations are not sufficient, as a multitude of water constituents affect phosphate adsorption. Thus, for complex waters, like drinking water or waste water main competitors have to be iden-tied and modelled in order to predict the removal of certain ions.

Inuence of contact time
The contact time of silicate with the adsorbent prior to the addition of phosphate inuences the loadings (Fig. 6). For 150 mg L À1 mGFH P-loadings are reduced by 21 and 43% for 28 and 56 d compared to loadings at 7 d. No signicant inuence (<5%) was observed at high adsorbent doses of 600 mg L À1 . At the same time, a signicant increase (18%) of silicate loadings in the single and binary system can be observed when increasing the contact time from 7 to 28 d. No further increase was measured for 56 d. This contradicts results from Christl et al. 18 who found a steady increase of silicate loadings on hematite over a time period of 210 days. Using ATR-FTIR spectroscopy, they measured a shi from monomeric to oligomeric and polymeric surface complexes over adsorption time. The competitive behaviour of chromate and silicate was studied by Zachara et al. 42 Chromate adsorption was greatly reduced by the presence of silicate and the reduction increased with  With mGFH as adsorbent and initial silicate concentrations of 10 mg L À1 , phosphate access to adsorption sites is not completely blocked, as can be inferred from only partially reduced phosphate loadings in Fig. 5. However, phosphate adsorption was reduced to 2.2 mg g À1 at low adsorbent doses and long contact times. As there is no difference in silicate adsorption between a contact time of 28 and 56 d, but phosphate adsorption is further reduced, the formation of a coating layer preventing access to internal pores and hindering adsorption is likely.  The order of addition does not lead to signicantly differing loadings for phosphate or silicate at pH values of 7.0 and 8.0 at contact times of 7 d (Fig. S5 †). At a pH of 6.0 however, the adsorbent surface is positively charged and thus favours anion adsorption. If phosphate (c eq ¼ 2 mg L À1 ) is added at a pH of 8.0 prior to silicate, resulting P-loadings are 26% higher than for later addition. Adsorption modelling using the double layer model (MINEQL5.0 with hydrous ferric hydroxide as adsorbent surface and c(Fe) ¼ 76 mg L À1 , c(P) ¼ 2 mg L À1 , c(Si) ¼ 10 mg L À1 , ionic strength ¼ 10 mmol L À1 ) predicts Fe-HPO 4 À (approx. 60%) and Fe-PO 4 2À (approx. 35%) as main resulting surface species (Fig. S6 †). Only at pH values above 8.5 Fe-PO 4

2À
predominates. Silicate adsorbs exclusively as Si(OH) 4 up to pH 8.0. When silicate adsorbs rst, it shis the pH PZC of iron hydroxides to lower values 14 and thus impedes phosphate adsorption. Also, through the addition of silicate, previously adsorbed phosphate is replaced.

Conclusions
Competing effects of silicate strongly inuence phosphate adsorption onto mGFH. Phosphate loadings on mGFH were reduced by 14, 23 and 41% in the presence of 6, 10 and 18 mg L À1 silicate, respectively. If the number of adsorption sites is limited, phosphate outcompetes silicate while the inuence of silicate rises with increasing numbers of available adsorption sites. Thus, monitoring of silicate concentrations and adjustments of adsorbent dose in phosphate elimination processes are necessary to ensure high P-removals. In contrast to phosphate, silicate bonds are mostly irreversible due to the proposed formation of Si-O-Si bonds on the surface of ferric hydroxides. Phosphate adsorption is favoured compared to silicate at lower pH values below the pH PZC . Binary isotherms could be tted using the extended Freundlich equation assuming multi-layer coverage. Thus, resulting loadings of both adsorbates in competition can be predicted using the tted model parameters. The models provide means to control phosphate elimination in water works applications. Real water samples contain a multicity of water constituents and additional effects are expected. However, the gained knowledge helps to understand the processes taking place in real water matrices.
Phosphate adsorption highly depends on the contact time of the iron hydroxide surface in silicate solution. An increase in contact time prior to phosphate addition from 7 to 28 and 56 d led to a decrease of 21 and 43% in P-loadings, respectively. Simultaneously, only minor silicate replacement by phosphate was observed, further underlining partially irreversible silicate adsorption. Thus, the formation of a coating layer resulting from silicate polymerization is likely to hinder access to adsorption sites reducing successful phosphate removal in applications with high silicate concentrations.

Conflicts of interest
The authors declare no conict of interest.