| Home | E-Submission | Sitemap | Contact Us |  
Environ Eng Res > Volume 21(2); 2016 > Article
Vasudevan, Ajithkumar, Singh, and Natarajan: Mass transfer kinetics using two-site interface model for removal of Cr(VI) from aqueous solution with cassava peel and rubber tree bark as adsorbents


Present study investigates the potential of cassava peel and rubber tree bark for the removal of Cr (VI) from aqueous solution. Removal efficiency of more than 99% was obtained during the kinetic adsorption experiments with dosage of 3.5 g/L for cassava peel and 8 g/L for rubber tree bark. By comparing popular isotherm models and kinetic models for evaluating the kinetics of mass transfer, it was observed that Redlich-Peterson model and Langmuir model fitted well (R2 > 0.99) resulting in maximum adsorption capacity as 79.37 mg/g and 43.86 mg/g for cassava peel and rubber tree bark respectively. Validation of pseudo-second order model and Elovich model indicated the possibility of chemisorption being the rate limiting step. The multi-linearity in the diffusion model was further addressed using multi-sites models (two-site series interface (TSSI) and two-site parallel interface (TSPI) models). Considering the influence of interface properties on the kinetic nature of sorption, TSSI model resulted in low mass transfer rate (5% for cassava peel and 10% for rubber tree bark) compared to TSPI model. The study highlights the employability of two-site sorption model for simultaneous representation of different stages of kinetic sorption for finding the rate-limiting process, compared to the separate equilibrium and kinetic modeling attempts.

1. Introduction

Due to limited water resources in the coming era, the role of water reuse, especially treated wastewater in agriculture is highly promising. The status of irrigated and rain-fed agriculture in the world is directly alarming us to the need for better irrigation scheduling to select optimum decision [13]. As far the material recovery and reuse strategies are concerned, the industrial effluents provide a better choice, in addition to the required environmental compliance. Cadmium, zinc, copper, nickel, lead, mercury, chromium are some of the metals often detected in industrial wastewater originating from mining, paint manufacturing, printing, smelting, pesticides, tanneries, petroleum, battery manufacturing, etc. [4, 5]. However, it is also important to note that toxic compounds like hexavalent chromium possess carcinogenic properties and can lead to corrosive effects in the intestinal tract of living beings [6, 7]. Because of their hazardous nature, they have to be removed before being discharged into the lakes, streams or any other water bodies.
Since the traditional methods for heavy metals removal such as chemical precipitation, solvent extraction, ion-exchange, reverse osmosis, etc. being very costly, industrialists and researchers are always interested in low-cost, eco-friendly and sustainable alternatives [810]. Considering the water resources constraints, reuse of treated wastewater for irrigation is demanding the continuous efforts to maintain good chemical quality in rivers and groundwater. In this regard, it is interesting to explore the viability of crop residues for the treatment of wastewater, which can, in turn, be used to irrigate the same crops. Among the physico-chemical treatment processes available, adsorption has been found to be effective in removing metals from aqueous solutions. Research on sustainable technology for the removal of Cr (VI) from wastewater has resulted in extensive evaluation of adsorption properties of natural materials such as modified wool [11, 12], orange peel [13], different types of sawdust [1421], exhausted coffee and tea [22], different types of moss [23, 24], peanut hull carbon [25, 26], palm [27], rice husk [28], durian shell waste [29], landfill clay [30], guar gum [31], Spirulina sp. [32], quaternary ammonium salt [33] and cork waste [34].
In this study, the efficacy of cassava peel and rubber bark powder for the removal of Cr (VI) metal ions has been analyzed which are largely grown in southern India as plantation crops. Considering the irrigation water requirement, these crops are highly tolerant to salinity, heavy metals and other water quality issues [3537]. Cassava waste consists of ligands such as hydroxyl, sulphur, cyano and amino groups which could bind heavy metal ions. Considering the usage of cassava peel as an adsorbent, recently a few studies have been reported about the biosorptive capacity of cassava peel [3843]. However, related to the adsorption studies using rubber tree bark, we could not find any literature except the study of chemically activated rubber wood sawdust for the removal of copper from industrial wastewater [4446].
While exploring novel adsorbents for developing an ideal adsorption system, it is essential to establish the most appropriate adsorption equilibrium correlation [47, 48], which is indispensable for reliable prediction of adsorption parameters. In this perspective, an understanding of the adsorption isotherm models can give insights into the pollutant-adsorbent interactions for optimizing the mass transfer mechanism. Linear regression method is the most widely approach to determine the most fitted isotherm model, which helps in evaluating the model parameters based on their goodness fit to the experimental data. However, use of kinetic adsorption models has become crucial for better understanding the mass transfer phenomena at the solid-liquid interphase which can simplify the estimation of rate-limiting step [4951].
The study of chemical kinetics even in homogeneous systems is quite complex and arduous. The investigation of rates of chemical reaction and molecular processes are important to determine the extent of equilibrium mass transfer as well as to explore the possibilities of mass transfer limitations [52, 53]. Since these kinetic models often describe the reaction at the solid-solution interface, varying surface phenomena has led to the development of multiple site models. However, as natural sorbents frequently exhibit functionally non-uniform surfaces, heterogeneous reactive site rate models have been suggested over simplified rate expressions [5456]. In several experiments conducting adsorption of heavy metals onto adsorbents, multi-site adsorption models were proposed to represent the temporal dependency of mass transfer kinetics [49]. It is based on the classification of the sorption sites into instantaneous and time-dependent fractions representing the effect of heterogeneity on mass transfer. The pseudo-second order kinetic equation has always been associated with the model of two-site-occupancy adsorption.
In this study, we demonstrate the applicability of different multiple-site models to explain temporal variation in the removal efficiency of the heavy metals with different bio-sorbents. The focus of this paper is to understand the kinetics of mass removal of Cr (VI) from aqueous solution using cassava peel and rubber tree bark. As far as the literature survey is concerned, cassava peel has been introduced for the adsorption studies only very recently [3843] but the usage of the same for the removal of Cr (VI) metal ions from the aqueous solution is yet to be understood. It is also observed from literature that the usage of rubber bark tree for the removal of heavy metals from wastewater is very limited. A comparison of the results with the classical isotherm models is also intended to reveal the importance of addressing the heterogeneity in mass transfer due to surface properties.

2. Materials and Methods

In order to evaluate the performance of the selected adsorbents in the aqueous solution of heavy metals, it is necessary to establish the optimum conditions of various process parameters through batch experiments. Realistic estimation of these variables is crucial in explaining the mass transfer rate and mechanism during biosorption.

2.1. Adsorbent Preparation

The adsorbent materials (cassava peel and rubber tree bark) were collected in their raw form in small strips (4–5 cm length) from the field and were subjected to a series of washing, dehydration and dehumidification in order to overcome the problems during wet grinding to make into powder form. Due to the surface modification resulted from temperature alterations, they were ground properly in this way and the materials are segregated over a series of 600, 425, 300 and 150 micron Indian Standard sieves. Based on the physical analysis, cassava peel possesses comparatively larger apparent density (7.8%) and surface area (24.2%) than rubber tree bark although they were ground to the same average particle size (Table 1). Similarly, stock solution of adsorbate (Cr (VI)) was prepared by taking 141.4 mg of potassium dichromate (Merck India) and dissolved up to 1000 ml mark of the volumetric flask to give a stock solution of concentration 50 mg/L. This stock solution was diluted to the desired concentration for different experiments.

2.2. Initial Batch Experiments

The critical environmental parameters which influence the adsorption phenomenon, viz. pH, adsorbent dose, time of contact and effect of concentration of Cr (VI) ions were investigated by conducting a series of batch experiments by varying one parameter at a time. These experiments were conducted using temperature-controlled orbital shaker (Trishul Equipment, Mumbai) at 170 rpm at a temperature of 30 ± 10°C by adding varying quantity of adsorbents into 25 mL of synthetic metal sample taken in 250 ml conical flask. The samples were filtered through Whatman A40 filter paper and the filtrates were analyzed for residual Cr (VI) by using Spectrophotometer (Genesys 20, Thermo Spectronic, India) at 540 nm wavelength after complexing with diphenylcarbazide.
The pH was optimized after studying for a range between 2 and 11 (Orion 420A+ pH meter, Thermo Electron Corporation, India) for maximum dissolved mass removal from the aqueous solution having an initial concentration of Cr (VI) as 50 mg/L, dose of cassava peel as 3.5 g/L and dose of rubber tree bark as 8 g/L. In the next step, the effect of adsorbent dose at the selected pH was investigated by varying biomass concentration in the range of 1 to 20 g/L for each adsorbent in a similar way. The effect of contact time was determined by drawing aqueous samples after 5, 10, 15, 20, 30, 45, 60, 90, 120, 150, 180, and 240 min which are filtered and analyzed for the residual metal concentration. Further, the impact of initial metal ion concentration on the adsorption efficiency was investigated by varying the initial Cr (VI) concentration from 60 to 350 mg/L for keeping the other conditions at optimum.

2.3. Batch Isotherm Experiments

Batch isotherm experiments were carried out in 250 m conical flask kept at 30°C in a horizontal shaker by varying the initial metal concentrations from 60 to 350 mg/L under optimal contact time (240 min) and optimum adsorbent dose (3.5 g/L of cassava peel and 8.0 g/L of rubber tree bark). By measuring the concentrations at the initial time and at equilibrium, maximum adsorption capacity as well as removal efficiency were calculated. Further experiments were conducted by changing the pH of the solution before addition of biomass by using 1N HCl, 0.1N HCl, 1N NaOH and 0.1N NaOH in order to check the suitability of various isotherm models under study for adsorption of Cr (VI) ions by both the adsorbents.

2.4. Batch Kinetic Experiments

Batch kinetic experiments were performed at an initial Cr (VI) concentration of 50 mg/L by varying contact time viz. 2, 5, 10, 15, 30, 45, 60, 90, 120, 180, 240 and 300 minutes. Samples were withdrawn at different time intervals and were analyzed for residual metal ions concentration. Further experiments were conducted for different initial Cr (VI) metal concentrations viz. 11.5, 27.63, 49.23 and 97.8 mg/L under optimum adsorbent dosage. The effect of variation of adsorbent dosage for different initial metal concentration on the kinetics of adsorption was also analyzed in this study. By conducting triplicate experiments, the percentage mean error during experiments was found to be about ± 10% from the average value.

3. Results and Discussion

3.1. Evaluation of Optimum Conditions for Maximum Mass Removal

Earlier studies indicated solution pH as an important parameter affecting biosorption of heavy metals [51, 54]. It was observed that maximum removal efficiency of Cr (VI) was found to be 99.88% for rubber bark and 98.12% for cassava peel at pH of 2 (Fig. 1). This behavior can be primarily attributed to the presence of surface-charged functional groups of the biomass which in turn influenced by the pH of the solution. Higher Cr (VI) adsorption at lower pH (~2) could be explained in terms of the electrostatic attraction between the negatively charged chromium species and positively charged functional groups on the surface of biomass cell wall. But at higher pH values, the number of functional groups carrying a net negative charge is more, which tend to repel the anions [54]. However, the reduction in removal (30%) is observed at pH values beyond 2, which is concerned with instantaneous overcharging of surfaces and subsequent charge reversals [5456]. The results also signify the need for acidification as a suitable pre-treatment for enhancing the mass removal at higher pH range.
The maximum biosorption of Cr (VI) ions was attained at adsorbent dose of 3.5 g/L and 8 g/L of the biomass for cassava peel and rubber tree bark respectively. The synchronous increase in removal efficiency is quite obvious because of the increase in adsorption site at a higher adsorbent dosage. Further, when the dosage is increased beyond optimum dose, there is little or no further adsorption of metal ions as a consequence of partial aggregation of biomass, resulting in decrease in effective surface area available for biosorption. However, it is important to maintain maximum removal at optimum dosage, thereby minimizing the costs involved and making it a favorable and cheap alternative.
It was observed that the removal efficiency of Cr (VI) ions onto both the adsorbents increased with increase in contact time and attained a maximum removal at 240 min (Fig. 2). Beyond optimum contact time, there was significant reduction in removal efficiency due to the possible desorption [51]. The effect of initial metal ion concentration was investigated by varying the initial concentration in the range of 60 to 350 mg/L for Cr (VI). It was observed that at lower concentrations, adsorption efficiency varied significantly and the optimum dose was found slightly higher than the expected doses. The results are in agreement with the similar studies in literature [22, 26].

3.2. Evaluation of Equilibrium Models

Evaluation of inter-phase mass transfer rate can be simplified by the assumption of local equilibrium (LEA) which represents idealized sorption with instantaneous mass transfer [57]. The differential uptake of heavy metals onto the surface of the biosorbent can be effectively characterized by separately evaluating the kinetic and equilibrium mass transfer. As seen from the physiographic properties of the selected biosorbents, it is hard to assume which isotherm they would follow since the activity of the sorption sites can invariably change during the mass transfer. Hence, a comparison was performed here with five equilibrium models such as Langmuir, Freundlich, Temkin, Dubinin-Radeshkevich and generalized isotherm to assess the most suitable equilibrium model for Cr (VI) adsorption onto cassava peel and rubber tree bark.

3.2.1. Langmuir isotherm model

Langmuir equation is based on the assumption of a structurally homogeneous adsorbent where all sorption sites are identical and energetically equivalent [58] which can be expresses as:
Where qe is solid phase sorbate concentration at equilibrium (mg/g), Ce is aqueous phase sorbate concentration at equilibrium (mg/L), Qm is Langmuir isotherm constant (mg/g), and b is Langmuir isotherm constant (L/mg). In adapting the Langmuir isotherm to solution chemistry, the idea of an upper limit to surface adsorption has been verified. The linearity of isotherm and high coefficient of determination (R2) values (0.986 for cassava peel and 0.997 for rubber tree bark) indicate that adsorption follows Langmuir model over the selected concentration ranges (Table 2). However, it was observed that the isotherm models were best fitted at pH 2 for Cr (VI) for both the adsorbents. The separation factor (RL= 1/(1+b Co)) was found to be varying between 0.21 and 0.58 for cassava peel, and between 0.22 and 0.62 for rubber tree bark, indicating the suitability of Langmuir isotherm for the range of initial Cr (VI) concentrations. This was further confirmed by the close fit between experimental and model values of Qm (Fig. 3).

3.2.2. Freundlich isotherm model

Freundlich isotherm model describes reversible adsorption and is not restricted to the formation of a monolayer unlike Langmuir isotherm model. The empirical equation can be expressed as:
Where KF is Freundlich constant (mg/g) and 1/n is the heterogeneity factor. Based on the experimental results, cassava peel showed a better adsorption (R2 = 0.989) compared to rubber tree bark (R2 = 0.907) (Table 2). This can be inferred from the difference in the physiographic properties of the adsorbents such as particle density and surface area (Table 1) and is confirmed by previous studies [59, 60].

3.2.3. Redlich-Peterson (R-P) isotherm model

A hybrid mechanism of adsorption was proposed by Redlich-Peterson (R-P) isotherm model by incorporating both the Langmuir and Freundlich equations, which can be expressed as:
Where KR is R-P isotherm constant (L/g), αR is R-P isotherm constant (l/mg) and β is the exponent which lies between 1 and 0. A general trial-and-error procedure was used to determine the best coefficient of determination (R2) as well as the three isotherm constants (KR, αR and β) from the linear plot of ln(Ce) against ln[KR(Ce/qe) – 1] by varying the value of KR, to yield maximum optimized value of R2 (Table 2, Fig. 3). Based on the linear regression analysis, R-P isotherm was most suitable for both cassava peel and rubber tree bark, followed by Langmuir and Freundlich isotherm. We hypothesis that, this can be either due to homogeneous adsorption independent of temporal variation in surface properties or existence of multi-sites for enhancing the mass transfer.

3.2.4. Temkin isotherm model

Temkin isotherm assumes that the fall in the heat of sorption to be linear rather than logarithmic, as implied in the Freundlich equation [51]. It can be expressed as:
Where T is the absolute temperature in Kelvin and R is the universal gas constant (8.314 J mol−1K−1). The constant bT is related to the heat of adsorption and AT is the equilibrium binding constant (Lmol−1) corresponding to the maximum binding energy. The applicability of the Temkin isotherm was analyzed using a linear plot between qe and ln Ce at constant temperature and the model constants (bT and AT) were evaluated (Table 2). The regression analysis indicated a lesser linear fit compared to the previous isotherm models for both adsorbents. This gives insight to the relative influence of chemical kinetics over the thermodynamic effect where the heat exchange is minimal compared to the inter-phase mass transfer [48].

3.2.5. Dubinin-Radeshkevich (D-R) isotherm model

The D-R model has been generally used to express adsorption based on pore filling mechanism with a Gaussian energy distribution on heterogeneous surface, and is written as:
qe=Qmexp (-Kɛ2)
Where K (mol2 kJ−2) is a constant which relates to the adsorption energy; and the model parameter ɛ can be calculated as:
ɛ=RTln (1+1Ce)
The values of Qm obtained from D-R model were found to be less compared to the Langmuir model prediction (27.2% and 19.6% for cassava peel and rubber tree bark respectively) (Table 2). It also resulted in poor linear fit for both adsorbents compared to other isotherms under consideration. This can be attributed to the specific nature of the adsorbents showing high surface uniformity where pore filling is not favored at high metal concentrations, thereby violating the assumptions of D-R model. The results also confirmed that the physical and chemical processes of biosorption cannot be distinguished easily with these model parameters.

3.2.6. Generalized isotherm model

Linear form of the generalized isotherm can be written as [51]:
log(Qmqe-1)=log Kg-Nblog Ce
Where KG is the saturation constant (mg/L); and Nb is the cooperative binding constant. The results showed that the model is fitting well for rubber tree bark (R2 = 0.985), while poorly fitting for cassava peel (R2 = 0.630) for the adsorption of Cr (VI) (Table 2). This can be due to the high Qm values for cassava peel as obtained from the Langmuir model showing the sensitivity of the model towards saturation capacity which is purely sorbent-specific.
From the values of coefficient of determination based on linear regression analysis (Table 2), it can be inferred that R-P model fitted well for adsorption of Cr (VI) ions onto both the adsorbents. However, maximum adsorption capacity obtained from this model for Cr (VI) adsorption was found to be higher for cassava peel (79.37 mg/g) followed by rubber tree bark (Qm = 43.86 mg/g) which were in close agreement to the experimental values (within 10% error margin), confirming the suitability of the Langmuir model for both cassava peel and rubber tree bark (Fig. 3). This is in close agreement with the published results on similar studies [38, 44].

3.3. Testing of Kinetic Adsorption Models

In order to investigate the mechanism of sorption such as mass transport and chemical reaction processes and to evaluate potential rate-controlling steps, six kinetic models (First Order kinetic model, Second Order kinetic model, Pseudo-First Order model, Pseudo-Second Order kinetic model, Intra-Particle Diffusion model and Elovich kinetic model) were tested in this study. The experiments were performed by varying the contact time, by keeping the adsorbent dose, initial concentration and temperature constant and the details of the linear regression fit were presented in Table 3.

3.3.1. First order (FO) kinetic model

From the point of view of fitting with linear regression models, a linear plot of t versus ln(1-(q/qe)) showed a lesser degree of fit (<0.95) for both adsorbents from the experimental data (Fig. 4). However, the adsorption of Cr (VI) ions exhibited saturation kinetics after the first hour as the quantity of bound metal ions reached saturation for both the adsorbents.
The variation in the initial Cr (VI) had resulted in significant increase in slope and intercept values for cassava peel compared to the rubber tree bark (Fig. 5). It showed that when the Co of Cr (VI) increased by 140.3%, 78.2% and 98.7%, the intercept of FO model for cassava peel also increased by 59%, 19.5% and 20.7% respectively, whereas the increase with rubber tree bark was quite marginal (8.8%, 12.9% and 14.3%). This can be attributed to the limit of saturation of sorption sites at higher dosage as observed in Fig. 4. Similar results were also found in literature regarding the impact of higher dosage on adsorption efficiency [42, 61].

3.3.2. Pseudo first order (PFO) kinetic model

As an alternative to the previous model, the PFO model was applied to the experimental results by plotting log (qeqt) versus t and the model fitness was compared in terms of R2 value and the prediction of rate constant, k1. It was observed that PFO model gave relatively a better fit to the experimental results on cassava peel than those for rubber tree bark (Table 3 and Fig. 6). Yet, it cannot be suggested as a representative model due to the validity over a limited period of time and large variation in equilibrium adsorption capacity (qe).

3.3.3. Second order (SO) kinetic model

It has been commonly observed that the sorption rate is very rapid at the beginning of the process and becomes slower as equilibrium is approached. In order to understand the significance of chemisorption, the experimental results were fitted with SO kinetic model. The experimental data was found to be best fitted for removal of Cr (VI) by cassava peel (R2 = 0.99), but to a lesser extent for rubber tree bark (R2 = 0.97) (Table 3, Fig. 7(a)).

3.3.4. Pseudo second order (PSO) kinetic model

Results from the linear plot between t/qt and t indicated a very close fitting (>0.999) of experimental data for PSO kinetic model (Table 3, Fig. 7(b)). In addition, the calculated qe values also agree well with the experimental qe values. Hence, the adsorption of Cr (VI) onto cassava peel and rubber tree bark can be described well by this model, indicating chemisorption to be the rate limiting step. This is in accordance with similar studies on adsorption with heavy metals [23, 24].

3.3.5. Intra-particle diffusion (IPD) model

As observed so far, the chemical potential of the sorbent-metal interactions were found to have an upper hand in controlling the mass transfer compared to physical mechanisms like pore filling or intra-film diffusion. On the basis of the change in availability of metal ion concentration over a period of time, the significance of intra-particle diffusion was analyzed by estimating the model constants ki from the slope of the plots of qt versus t1/2. The values of intercept represent the boundary layer thickness, which in effect account for the increase in adsorption capacity. Even though regression analysis resulted in lowest fit (< 0.85), the intercept was observed to be increasing with increasing initial metal concentration (Table 3) which can be interpreted as the rate-limiting step for the adsorption of Cr (VI) unto both the adsorbents.
As seen from Fig. 8, the trend of the experimental data revealed existence of multi-linearity where two or more steps may be involved in adsorption [62, 63]. The initial portion represents instantaneous adsorption stage, followed by a gradual adsorption stage where intra-particle diffusion seems to be rate-controlling. The third portion is the final equilibrium stage where intra-particle diffusion starts to slow down due to the low adsorbate concentrations in the solution. This indicates that transport of metal ions from the solution (through the sorbent-liquid interface) to the available surface of sorbent as well as to the micro-pores may be responsible for the high uptake of metal ions. Similar findings were also reported by Namasivayam and Ranganathan [9] for the uptake of Cadmium on a waste product from fertilizer industry. The deviation of the curves from the origin also indicated that intra-particle transport is not the only rate-limiting step in adsorption of Cr (VI) ions onto selected adsorbents in this study.

3.3.6. Elovich kinetic model

The kinetics of chemisorption has been adequately represented in the literature by the second-order equation and the Elovich Equation [64]. It was observed that Elovich kinetic model can adequately represent the adsorption of Cr (VI) onto both the adsorbents based on the regression analysis (Table 3, Fig. 9). On the basis of the excellent fitting of both the PSO and Elovich model, it can be proposed that the predominant sorption mechanism may be chemisorption.

3.4. Comparison with Two-Site Sorption Model

The two-site sorption conceptual model is based on assumption that a fraction (f) of the mass exchange sites of the sorbents favours instantaneous sorption (S1 = fKC) and another fraction (1 – f) is considered to be time-dependent (S2 = (1 – f)KC). Two different approaches were considered in this study depending on the mode of contact of the sorption sites with the liquid concentration namely, two-site sorption model with a serial interface (TSSI) and two-site sorption model with a parallel interface (TSPI) [24]. By considering the reversible sorption kinetics between the two types of sites, TSSI model can be expressed as:
Where K1 and K2 are the forward and backward sorption rate coefficients and S1 and S2 are the associated sorbed mass respectively. Similarly, TSPI model can be used to estimate equilibrium and non-equilibrium sorption parameters as given by:
Where α and K are the sorption rate constant (h−1) and distribution coefficient (L/kg−1), respectively. By applying these kinetic models to the batch system with the assumption that time-dependent fraction obey FO kinetic model, Lee et al. [24] developed analytical solutions for TSSI and TSPI models. We have further modified them to express in the linear form as:
ln(C*R-1)=-(Krβ)t+ln(1β-1)   (for TSSI model)
ln(C*-(1R))=-(αβ)t_ln(1-β)   (for TSSI model)
Where C*=Ce/C0; retardation factor, R=1+ Π MK/V; is the mass of sorbent (g); Kr is the sorption rate constant for TSSI model (h−1); V is the volume of the solution (mL); α is the sorption rate constant (h−1) for TSPI model and β is a model parameter corresponding to equilibrium fraction.
The possibility of two-site sorption models was evaluated based on the linear regression analysis of Eq. (10) and Eq. (11) for both adsorbents. It was observed that both sorbents exhibited very good fit for both two-site models even though TSPI model fitted more accurately for sorption unto cassava peel whereas TSSI model was more suitable for representing sorption unto rubber tree bark (Fig. 10). Considering the specific kinetic nature of sorption in terms of model parameters, TSPI model was found to be invariant for both the sorbents with high mass transfer rate (0.89 h−1), whereas TSSI model resulted in low mass transfer rate (5% for cassava peel and 10% for rubber tree bark) compared to TSPI model (Table 4). However, based on the surface properties of the adsorbents (size and shape of the particles), it is suggested that TSPI model can hold well under different conditions.
By comparing the model results with separate equilibrium and kinetic modeling attempts, it can be inferred the sequential process of mass transfer can be more effectively described by two-site model. During the initial stage, the exterior surface reached the saturation level due to the direct attachment to the surface; then the metal ions begin to attach to the adsorbent via the pores within the particles and were adsorbed onto the interior surfaces of the particles. Until this stage, isotherm model holds well, while the latter part of kinetic sorption can act as the rate-limiting step as explained by the chemi-sorption. With the decrease in metal ion concentration in the solution, the diffusion rate became constantly lower, and consequently, the adsorption process end up in equilibrium conditions.

4. Conclusions

Present study discusses about the kinetics of Cr (VI) removal from the aqueous solution by using two naturally available agro-waste materials (cassava peel and rubber tree bark). Apart from the estimation of removal efficiency based on batch (equilibrium and kinetic) adsorption experiments this study focuses on the suitability of various modeling approaches to evaluate the kinetics of mass transfer. A comparison of six isotherm models and six kinetic models were presented based on the linear regression fitting. The study also highlights the possibility of employing multi-site sorption model for simultaneous representation of different stages of kinetic sorption as well as identifying the rate-limiting process. The salient results of the present study are briefed here.
With the increase in the cassava peel adsorbent dose from 0.25 to 3.5 g/L the percentage removal of Cr (VI) ion increased from 26.33 to 99.95%. The percentage of adsorption increased for rubber tree bark from 1.83 to 100% as the dose increased from 0.25 to 8 g/L due to the availability of more unsaturated adsorption sites. When the dosage is increased beyond optimum dose, there is little or no further adsorption of metal ions as a consequence of partial aggregation of biomass resulting in decrease in effective surface area available for the biosorption.
Based on the comparison of isotherm models, Redlich-Peterson model and Langmuir model fitted well (R2 > 0.99) for both sorbents indicating that an upper limit existed for the adsorption over the external surface. The maximum adsorption capacity for cassava peel and rubber tree bark were found to be 79.37 mg/g and 43.86 mg/g respectively. Comparison with the Temkin isotherm model showed that the influence of chemical kinetics dominated the thermodynamic effect where the thermodynamic effect is minimal compared to the inter-phase mass transfer.
While comparing the kinetic sorption models, pseudo-second order model showed a very good fit as well as closeness to the experimental values of sorption capacity, indicating chemisorption to be the rate limiting step. This was further confirmed by the fitness to Elovich model. Unlike similar studies reporting intra-particle diffusion as the rate-limiting step, the diffusion model resulted in multi-linearity indicating the possibility of two or more steps in adsorption.
Two different multi-sites models namely, two-site series interfaces (TSSI) and two-site parallel interface (TSPI) were considered in this study depending on the type of interaction of sorption sites with the liquid concentration. It was observed that both sorbents exhibited very high fit for both two-site models. Considering the specific kinetic nature of sorption in terms of model parameters, TSPI model was found to be invariant for both the sorbents with high mass transfer rate (0.89 h−1), whereas TSSI model resulted in low mass transfer rate (lower by 5% for cassava peel and 10% for rubber tree bark) compared to TSPI model.
By comparing the model results with the conventional explicit equilibrium and kinetic modeling attempts, it can be inferred that the sequential process of mass transfer can be more effectively described by two-site model. It is also inferred that the adsorbents selected in this study (cassava peel and rubber tree bark) are efficient for the removal of Cr (VI) as a suitable alternative to other costly adsorbents.


1. Valipour M. A comprehensive study on irrigation management in Asia and Oceania. Arch Acker Pfl Boden. 2015;61:1247–1271.

2. Valipour M, Ziatabar AM, Raeini-Sarjaz M, et al. Agricultural water management in the world during past half century. Arch Acker Pfl Boden. 2015;61:657–678.

3. Valipour M. What is the tendency to cultivate plants for designing cropping intensity in irrigated area. Adv Water Sci Tech. 2015;2:1–12.

4. Williams CJ, Aderhold D, Edyvean GJ. Comparison between biosorbents for the removal of metal ions from aqueous solutions. Water Res. 1998;32:216–224.

5. Kadirvelu K, Thamaraiselvi K, Namasivayam C. Removal of heavy metal from industrial wastewaters by adsorption onto activated carbon prepared from an agricultural solid waste. Bioresour Technol. 2001;76:63–65.
crossref pmid

6. Hueper WC, Payne WW. Experimental studies in metal carcinogenesis. Arch Environ Health. 1962;5:445–462.
crossref pmid

7. Shanmugavalli R, Madhavakrishnan S, Kadirvelu K, Rasappan K, Mohanraj R, Pattabhi S. Adsorption studies on removal of Cr (VI) from aqueous solution using silk cotton hull carbon. J Ind Pollu Contr. 2007;23:65–72.

8. Viraraghavan T, Dronamrajum . Removal of copper, nickel and-zinc from wastewater by adsorption using peat. J Environ Sci HealthPart A. 1993;28:1261.

9. Namasivayam C, Ranganathan K. Removal of Pb (II), Cd (II) and Ni (II) and mixture of metal ions by adsorption onto waste Fe (III)/Cr (III) hydroxide and fixed bed studies. Environ Technol. 1995;16:851–860.

10. Ngah WS, Hanafiah MAKM. Removal of heavy metal ions from wastewater by chemically modified plant wastes as adsorbents: A review. Bioresour Technol. 2008;99:3935–3948.
crossref pmid

11. Netzer A, Hughes DE. Adsorption of Cr, Pb and Co by activated carbon. Water Res. 1984;18:927–933.

12. Reed BE, Arunachalam S. Use of granular activated carbon columns for lead removal. J Environ Eng. 1994;120:416–436.

13. Masri MS, Friedman M. Effect of chemical modification of woolon metal ion binding. J Appl Polym Sci. 1974;18:2367–2377.

14. Dikshit VP. Removal of chromium (VI) by adsorption using sawdust. Natl Acad Sci Lett. 1989;12:419–421.

15. Bryant PS, Petersen JN, Lee JM, Brouns TM. Sorption of heavy metals by untreated red fir sawdust. Appl Biochem Biotechnol. 1992;34–35:777–788.

16. Zarraa MA. A study on the removal of chromium (VI) from waste solutions by adsorption on to sawdust in stirred vessels. Adsorption Sci Technol. 1995;12:129–138.

17. Selvi K, Pattabhi S, Kadirvelu K. Removal of Cr (VI) from aqueous solution by adsorption onto activated carbon. Bioresour Technol. 2001;80:87–89.
crossref pmid

18. Baral SS, Das SN, Rath P. Hexavalent chromium removal from aqueous solution by adsorption on treated sawdust. Biochem Eng J. 2006;31:216–222.

19. Argun ME, Dursun S, Ozdemir C, Karatas M. Heavy metal adsorption by modified oak sawdust: thermodynamics and kinetics. J Hazard Mater B. 2007;141:77–85.

20. Bansal M, Singh D, Garg VK, Rose P. Mechanisms of Cr (VI) removal from synthetic wastewater by low cost adsorbents. J Environ Res Dev. 2008;3:228–243.

21. Vinodhini V, Das N. Mechanisms of Cr (VI) biosorption by Neem sawdust. Amer Euras J Scientif Res. 2009;4:324–329.

22. Orhan Y, Buyukgungor H. The removal of heavy metals by using agricultural wastes. Water Sci Technol. 1993;28:247–255.

23. Sharma DC, Forster CF. Removal of hexavalent chromium using sphagnum moss peat. Water Res. 1993;27:1201–1208.

24. Lee SM, Laldawngliana C, Tiwari D. Iron oxide nano-particles-immobilized-sand material in the treatment of Cu (II), Cd (II) and Pb (II) contaminated waste waters. Chem Eng J. 2012;195:103–111.

25. Periasamy K, Namasivayam C. Removal of nickel (II) from aqueous solution and nickel plating industry wastewater using an agricultural waste: Peanut hulls. Waste Manag. 1995;5:63–68.

26. Dubey SP, Gopal K. Adsorption of chromium (VI) on low cost adsorbents derived from agricultural waste material: a comparative study. J Hazard Mater. 2007;145:465–470.
crossref pmid

27. Iqbal M, Saeed A, Akhtar N. Petiolar felt-seath of palm: A new biosorbent for the removal of heavy metals from contaminated water. Bioresour Technol. 2002;81:151–153.
crossref pmid

28. Shraim AM. Rice is a potential dietary source of not only arsenic but also other toxic elements like lead and chromium. Arab J Chem. 2014;

29. Kurniawan A, Sisnandy VOA, Trilestari K, Sunarso J, Indraswati N, Ismadji S. Performance of durian shell waste as high capacity biosorbent for Cr (VI) removal from synthetic wastewater. Ecol Eng. 2011;37:940–947.

30. Ghorbel-Abid I, Trabelsi-Ayadi M. Competitive adsorption of heavy metals on local landfill clay. Arab J Chem. 2011;8:25–31.

31. Khan TA, Nazir M, Ali I, Kumar A. Removal of Chromium (VI) from aqueous solution using guar gum–nano zinc oxide biocomposite adsorbent. Arab J Chem. 2013;

32. Rezaei H. Biosorption of chromium by using Spirulina sp. Arab J Chem. 2013;

33. Yoganand KS, Umapathy MJ. Green methodology for the recovery of Cr (VI) from tannery effluent using newly synthesized quaternary ammonium salt. Arab J Chem. 2013;

34. Sfaksi ZN, Abdelwahab Azzouz A. Removal of Cr (VI) from water by cork waste. Arab J Chem. 2014;7:37–42.

35. Valipour M. Future of agricultural water management in Africa. Arch Acker Pfl Boden. 2015;61:907–927.

36. Valipour M. Land use policy and agricultural water management of the previous half of century in Africa. Ap Water Sci. 2015;5:367–395.

37. Valipour M. Pressure on renewable water resources by irrigation to 2060. Acta Adv Agri Sci. 2014;2:32–42.

38. Seepe L. The use of cassava waste in the removal of cobalt, chromium and vanadium metal ions from synthetic effluents. MS Thesis. Johannesburg: University of Witwatersrand; 2015.

39. Alinnor IJ, Nwachukwu MA. Adsorption of phenol on surface–– Modified cassava peel from its aqueous solution. Int J Environ Sci Manage Eng Res. 2012;1:68–76.

40. Simate GS, Ndlovu S. The removal of heavy metals in a packed bed column using immobilized cassava peel waste biomass. J Ind Eng Chem. 2015;21:635–643.

41. Isiuku BO, Horsfall M, Spiff AI. Removal of methyl red from aqueous solution by NaOH-activated cassava peels carbon in a fixed-bed column. Res J App Sci. 2014;9:238–243.

42. Gin WA, Jimoh A, Abdulkareem AS, Giwa A. Kinetics and Isotherm Studies of Heavy Metal Removals from Electroplating Wastewater Using Cassava Peel Activated Carbon. Int J Eng. 2014;3:

43. Owamah HI. Biosorptive removal of Pb (II) and Cu (II) from wastewater using activated carbon from cassava peels. J Mater Cycl Waste Manage. 2014;16:347–358.

44. Hanafiah MA, Ngah KM, Ibrahim WSW, Zakaria SC, Ilias H. Kinetics and thermodynamic study of lead adsorption onto rubber (Heveabrasiliensis) leaf powder. J Appl Sci. 2006;6:2762–2767.

45. Kalavathy MH, Karthikeyan T, Rajgopal S, Miranda LR. Kinetics and isotherm studies of Cu(II) adsorption onto H3PO4-activated rubber wood sawdust. J Colloid Interf Sci. 2005;292:354–362.

46. Masae M, Sikong L, Kongsong P, Phoempoon P, Rawangwong S, Sririkun W. Application of rubber wood ash for removal of Nickel and Copper from aqueous solution. Environ Nat Resour. 2013;11:17–27.

47. Abia AA, Didi OB, Asuquo ED. Modeling of Cd2+ sorption kinetics from aqueous solutions onto some thiolated agricultural waste adsorbents. J Appl Sci. 2006;6:2549–2556.

48. Foo K, Hameed BH. Insights into the modeling of adsorption isotherm systems. Chem Engg J. 2010;156:2–10.

49. Rudzinski W, Płazinski W. Theoretical description of the kinetics of solute adsorption at heterogeneous solid/solution interfaces on the possibility of distinguishing between the diffusional and the surface reaction kinetics models. Appl Surf Sci. 2007;253:5827–5840.

50. Chiban M, Carja G, Lehutu G, Sinan F. Equilibrium and thermodynamic studies for the removal of As (V) ions from aqueous solution using dried plants as adsorbents. Arab J Chem. 2011;

51. Gupta S, Babu BV. Removal of toxic metal Cr (VI) from industrial wastewater using sawdust as adsorbent: equilibrium and, kinetics and regeneration studies. J Chem Engg. 2008;141:1–20.

52. Natarajan N, Suresh Kumar G. Numerical modelling and spatial moment analysis of solute transport with Langmuir sorption in a fracture matrix-coupled system. ISH J Hydrau Eng. 2015;21:28–41.

53. Vasudevan M, Suresh Kumar G, Nambi IM. Numerical studies on kinetics of sorption and dissolution and their interactions for estimating mass removal of toluene from entrapped soil pores. Arab J Geosci. 2014;8:1–16.

54. Mohanty K, Jha M, Meikap BC, Biswas MN. Removal of chromium (VI) from dilute aqueous solutions by activated carbon developed from Terminaliaarjuna nuts activated with zinc chloride. Chem Eng Sci. 2005;60:3049–3059.

55. Singh RP, Zahra F, Savio W, Prasad SC. Axial dispersion and mass transfer controlled simulation study of chromium (VI) adsorption onto tree leaves and activated carbon. J Environ Eng. 2009;135:1071–1083.

56. Zhou M, Liu Y, Zeng G, Li X, Xu W, Fan T. Kinetic and equilibrium studies of Cr (VI) biosorption by dead Bacillus licheniformis biomass. World J Microbio Biotech. 2007;23:43–48.

57. Seagren EA, Rittmann BE, Valocchi AJ. A critical evaluation of the local-equilibrium assumption in modeling NAPL-pool dissolution. J Contam Hydrol. 1999;39:109–135.

58. Wong YC, Szeto YS, Cheung WH, McKay G. Effect of temperature, particle size and percentage deacetylation on the adsorption of acid dyes on chitosan. J Adsor. 2008;14:11–20.

59. Qiu H, Lv L, Pan BC, Zhang QJ, Zhang WM, Zhang QX. Critical review in adsorption kinetic models. J Zhejiang Univ Sci. 2009;10:716–724.
crossref pdf

60. Lee S, Kim DJ, Choi JW. Comparison of first-order sorption kinetics using concept of two-site sorption model. Environ Eng Sci. 2012;29:1002–1007.
crossref pmid pmc

61. Machida M, Kikuchi Y, Aikawa M, Tatsumoto H. Kinetics of adsorption and desorption of Pb (II) in aqueous solution on activated carbon by two-site adsorption model. Coll Surf A: Physic Eng Asp. 2004;240:179–186.

62. Annadurai G, Juang RS, Lee DJ. Use of cellulose-based wastes for adsorption of dyes from aqueous solutions. J Hazard Mater. 2002;92:263–274.
crossref pmid

63. Wu FC, Tseng RL, Juang RS. Adsorption of dyes and phenols from water on the activated carbons prepared from corncob wastes. J Environ Technol. 2001;22:205–213.

64. Ungarish M, Aharoni C. Kinetics of chemisorption: deducing kinetic laws from experimental data. J Chem Soc Faraday Trans. 1981;77:975–985.

Fig. 1
Effect of pH on biosorption of Cr (VI) onto Cassava peel and Rubber tree bark.
Fig. 2
Determination of effective time and optimum dose for an initial Cr (VI) concentration of (a) 25 mg/L for cassava peel, (b) 25 mg/L for rubber tree bark, (c) 50 mg/L for cassava peel, and (d) 50 mg/L for rubber tree bark.
Fig. 3
Comparison of computed qe from Langmuir isotherm model and R-P isotherm model with experimental qe for Cr (VI) removal by both adsorbents.
Fig. 4
(a) Kinetics of the adsorption and (b) First Order kinetic plot of Cr (VI) onto cassava peel and rubber tree bark
Fig. 5
First order kinetic plot for different initial concentrations onto (a) cassava peel, (b) rubber tree bark.
Fig. 6
Pseudo-first order kinetic plot for Cr (VI) adsorption onto both adsorbents.
Fig. 7
(a) Second Order kinetic plot and (b) Pseudo-Second Order kinetic plot for Cr (VI) adsorption onto cassava peel and rubber tree bark.
Fig. 8
Intra-particle diffusion plot for Cr (VI) adsorption onto both adsorbents.
Fig. 9
Elovich kinetic model plot for Cr (VI) adsorption onto both adsorbents.
Fig. 10
Comparison of multi-site sorption models for adsorption of Cr (VI) onto (a) cassava peel and (b) rubber tree bark. TSSI refers to Two-Site Series Interface and TSPI refers to Two-Site Parallel Interface.
Table 1
Physical Properties of Selected Adsorbents
Properties Cassava peel Rubber tree bark
Apparent density, kg/m3 619.01 574.12
Bed porosity 33.90 36.30
Particle density, kg/m3 937.03 902.36
Ash content, g/g 33.10 37.56
Surface area, cm2/g 897.28 722.54
Average particle diameter, μm 755.38 775.99
Table 2
Isotherm Constants and Regression Equation for Various Adsorption Isotherms for Adsorption of Cr (VI) on Selected Adsorbents
Sl. No Isotherm Adsorbent Linear equation R2 Model constants
1 Langmuir Cassava peel Ce/qe = 0.0123 Ce + 0.0535 0.986 Qm = 79.37, b =0.327
Rubber tree bark Ce/qe = 0.023 Ce + 0.0761 0.997 Qm = 43.86, b =0.294
2 Freundlich Cassava peel log qe = 0.237 log Ce +1.4653 0.989 Kf = 27.55, n =10.03
Rubber tree bark log qe = 0.198 log Ce +1.2371 0.907 Kf = 17.26, n =5.05
3 Redlich-Peterson Cassava peel ln((KRCe/qe)-1) = 1.06 ln Ce – 1.902 0.993 KR =50.0, aR =0.15, β =1.06
Rubber tree bark ln((KRCe/qe)-1) = 1.1568 ln Ce – 0.9254 0.996 KR =75.0, aR =0.40, β =1.16
4 Temkin Cassava peel qe =10.846 ln Ce + 29.104 0.947 bT =232.265, AT =14.635
Rubber tree bark qe =5.677 ln Ce + 16.202 0.942 bT =443.745, AT =17.357
5 Dubinin-Radeshke vich Cassava peel ln qe = -1×10−07ɛ2 + 4.056 0.681 Qm =57.754, K= 1 × 10−07
Rubber tree bark ln qe = -4×10−07ɛ2 + 3.618 0.756 Qm =35.252, K= 4 × 10−07
6 Generalized Cassava peel log((Qm/qe)-1) = −1.135 log Ce + 0.392 0.630 Nb =1.135, Kb = 2.465
Rubber tree bark log((Qm/qe)-1) = −0.96 log Ce + 0.4599 0.985 Nb =0.96, Kb = 2.883
Table 3
Kinetic Constants and Regression Equations for Kinetic Models for Cr (VI) Adsorption onto Cassava Peel and Rubber Tree Bark
Sl No. Kinetic model Adsorbent Linear equation R2 value Model constants
1 First order Cassava peel ln(1 – (q/qe) = −0.021t – 0.7993 0.959 K= 0.02
Rubber tree bark ln(1 – (q/qe) = −0.0212t – 0.7448 0.928 K= 0.02

2 Second order Cassava peel 1/(qe – qt) = 0.0118t + 0.089 0.993 K2 = 0.0118, qe=9.51
Rubber tree bark 1/(qe – qt) = 0.0346t + 0.0027 0.972 K2 = 0.0518, qe=2.13

3 Pseudo-first order Cassava peel log(qe – qt) = −0.0087t + 0.7418 0.978 k1=0.02, qe=5.51
Rubber tree bark log(qe – qt) = −0.014t + 0.6152 0.920 k1 = 0.032, qe=4.12

4 Pseudo-second order Cassava peel t/qt = 0.0707t + 0.4175 0.999 K2i =0.012, qe = 14.29
Rubber tree bark t/qt = 0.01594t + 1.0098 0.999 K2i =0.025, qe=6.29

5 Intra-particle diffusion Cassava peel qt =0.4176t1/2 + 8.0866 0.819 kid = 0.417
Rubber tree bark qt =0.2013t1/2 + 3.3723 0.803 kid = 0.2013

6 Elovich Cassava peel qt =1.5014lnt + 5.9346 0.956 β =0.67, α = 1121
Rubber tree bark qt =0.7368lnt + 2.288 0.973 β = 1.36, α =3.98
Table 4
Kinetic Constants and Regression Equations for two-Site Models (TSSI and TSPI) for Cr (VI) Adsorption onto Cassava Peel and Rubber Tree Bark
Two-site model Adsorbent Linear equation R2 value Model constants
TSSI Cassava peel ln(C*R–1) =−0.0002t + 0.253 0.968 β = 0.437; Kr = 0.053 (h−1)
Rubber tree bark ln(C*R–1) =−0.004t + 0.384 0.995 β = 0.406; Kr = 0.098 (h−1)
TSPI Cassava peel ln(C*–1/R) =−0.02t – 1.362 0.992 β = 0.744; α=0.892 (h−1)
Rubber tree bark ln(C*–1/R) =−0.02t – 1.368 0.988 β = 0.745; α=0.894 (h−1)
PDF Links  PDF Links
PubReader  PubReader
Full text via DOI  Full text via DOI
Download Citation  Download Citation
Web of Science
Editorial Office
464 Cheongpa-ro, #726, Jung-gu, Seoul 04510, Republic of Korea
FAX : +82-2-383-9654   E-mail : eer@kosenv.or.kr

Copyright© Korean Society of Environmental Engineers.        Developed in M2PI
About |  Browse Articles |  Current Issue |  For Authors and Reviewers