Next Article in Journal
Study on the Head Loss of the Inlet Gradient Section of the Aqueduct
Next Article in Special Issue
Comparison of Phenol Adsorption Property and Mechanism onto Different Moroccan Clays
Previous Article in Journal
Water Distribution Network Partitioning Based on Complex Network Theory: The Udine Case Study
Previous Article in Special Issue
Applying Linear Forms of Pseudo-Second-Order Kinetic Model for Feasibly Identifying Errors in the Initial Periods of Time-Dependent Adsorption Datasets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Factorial Design Statistical Analysis and Optimization of the Adsorptive Removal of COD from Olive Mill Wastewater Using Sugarcane Bagasse as a Low-Cost Adsorbent

1
Laboratory of Marine Biotechnologies and Environment, Sciences Faculty, Chouaïb Doukkali University, El Jadida 24000, Morocco
2
Science Engineer Laboratory for Energy, National School of Applied Sciences, Chouaïb Doukkali University, El Jadida 24000, Morocco
3
Chemical & Biochemical Sciences, Green Process Engineering, CBS, Mohammed VI Polytechnic University, Ben Guerir 43150, Morocco
4
Multidisciplinary Research and 12 Innovation Laboratory, FP Khouribga, Sultan Moulay Slimane University of Beni Mellal, Beni-Mellal 23000, Morocco
5
EA 4297 TIMR Laboratory, Department of Chemical Engineering, Centre de Recherche de Royallieu, University of Technology of Compiegne—Sorbonne University, Rue du Dr Schweitzer, CS 60319, 60200 Compiegne, France
6
Dynamical Systems Laboratory, Sciences Faculty, Chouaïb Doukkali University, El Jadida 24000, Morocco
*
Author to whom correspondence should be addressed.
Water 2023, 15(8), 1630; https://doi.org/10.3390/w15081630
Submission received: 20 March 2023 / Revised: 18 April 2023 / Accepted: 19 April 2023 / Published: 21 April 2023
(This article belongs to the Special Issue Adsorption Technology for Water and Wastewater Treatments)

Abstract

:
This work highlights the elimination of chemical oxygen demand (COD) from olive mill wastewater using sugarcane bagasse. A 25−1 fractional factorial design of experiments was used to obtain the optimum conditions for each parameter that influence the adsorption process. The influence of the concentration of sugarcane bagasse, solution pH, reaction time, temperature, and agitation speed on the percent of COD removal were considered. The design experiment describes a highly significant second-order quadratic model that provided a high removal rate of 55.07% by employing optimized factors, i.e., a temperature of 60 °C, an adsorbent dose of 10 g/L, a pH of 12, a contact time of 1 h, and a stirring speed of 80 rpm. The experimental data acquired at optimal conditions were confirmed using several isotherms and kinetic models to assess the solute interaction behavior and kind of adsorption. The results indicated that the experimental data were properly fitted with the pseudo-first-order kinetic model, whereas the Langmuir model was the best model for explaining the adsorption equilibrium.

1. Introduction

Olive mill wastewater (OMW) is an important environmental problem due to the strong color, low pH, and high concentrations of chemical oxygen demand (COD), biochemical oxygen demand (BOD), and phenolic compounds [1,2]. The discharge of this effluent into soil or rivers without treatment leads to damages to the environment that are attributed to their potential risk to flora or depletion of clean water reservoirs [3]. Hence, the treatment of OMW is a great challenge in a problematic environment. There are numerous approaches that can be used to remove organic pollutants and reduce COD levels from OMW such as coagulation–flocculation, aerobic and anaerobic biodegradation, solar distillation, adsorption, and infiltration percolation [4,5,6,7,8,9,10]. Among these techniques, adsorption is a well-established technique that benefits from advantages such as a simple design, low-cost, eco-friendly, high ability, non-generation of secondary pollutants, and reusability [11]. Activated carbon is one of the most efficient adsorbents used to treat many kinds of effluents. Nevertheless, it is relatively expensive, and its regeneration is difficult, which has necessitated the search for alternative adsorbents [12].
Therefore, several industrial and agricultural wastes exhibit considerable potential for the removal of pollutants. Several agricultural by-products are investigated for COD removal from wastewater, including bamboo [13], date pit [14], raw bagasse [15], areca catechu fronds [16], and date palm waste [17]. Among the efficient agro-industrial residues is sugarcane bagasse, which has a high cellulose (45%), hemicellulose (28%), and lignin (18%) content [18]. These materials contain reactive functional groups such as carboxylic and hydroxyl groups with the potential to adsorb organic loads [19]. Several studies have been conducted to investigate the sugarcane bagasse as an adsorbent for different pollutants, including Congo red, Pb(II), phenol from water, and phenol from OMW [20,21,22,23,24].
Although several studies have been conducted using sugarcane bagasse to remove several organic dye contaminants, no investigation has been carried out on the mathematical modeling and statistical optimization of the COD removal from OMW using sugarcane bagasse. Therefore, the present work aimed at the statistical optimization and modeling of key experimental parameters using fractional factorial design and response surface methodology (RSM) to attain optimum COD removal. Moreover, it recognizes the effects of five parameters, such as adsorbent dosage, pH, agitation time, stirring speed, and temperature as well as their interactions on the adsorption process of COD by sugarcane bagasse. Additionally, the adsorption kinetics, equilibrium, and thermodynamics were also studied.

2. Materials and Methods

2.1. OMW Sample

An OMW solution was obtained from an olive oil mill in Marrakech City, Morocco, during the campaign 2018/2019. This factory operates with the modern two- and three-phase olive oil extraction technology. The OMW used is characterized by a COD of 347.8 g(O2)/L, a conductivity of 15.3 ms/cm, a phenolic content of 15.29 g/L and a pH of 4.27. The supplied effluent was immediately kept refrigerated at −4 °C to prevent any alteration in its physicochemical characteristics.

2.2. Sugarcane Bagasse Material

Sugarcane bagasse was collected from a local plant in the region of El Jadida, Morocco. The collected adsorbent was oven-dried for 24 h at 70 °C before being ground with a domestic electrical miller and sieved through a 50 µm screen. The obtained powder was repeatedly washed with distilled water to remove all impurities, and it was then oven dried at 100 °C for 48 h (Figure 1). The dried powder was used in experiments without any further treatment.
The sugarcane bagasse functional groups were characterized by FTIR using Perkin-Elmer 1720-x spectrometer. The sample was mixed with KBr at a mass ratio of 1:100 and finely powdered and pressed to pellets. The infrared spectra were recorded in the range of 4000–500 cm−1 with 2 cm−1 resolution. Therefore, X-ray fluorescence analysis was employed to investigate the chemical composition of the adsorbent using the “Axion” spectrometer.

2.3. Batch Adsorption Experiments

All reagents used in the preparation and the adsorption studies were of analytical grade. The elimination of COD from OMW using sugarcane bagasse was studied in batch mode using shut-top pyrex bottles comprising 100 mL of OMW and an appropriate mass of adsorbent, which were stirred during the desired time in an incubator (Stuart SI50). The remaining COD content was determined from the UV-Vis absorbance characteristic using a Jenway 6320D UV/Vis spectrophotometer. The wavelength of maximum absorption (λmax) was 620 nm for this measurement. The quantity adsorbed and the removal efficiency of COD were calculated by measuring the solution’s concentration before and after adsorption using the following equations:
q e = ( C o C e ) V m
% Rem   = ( C o C e ) m × 100
where qe is the amount of COD removed by the adsorbent (mg/g); Co and Ce are, respectively, the concentrations of COD before and after the batch adsorption study (mg/L); m is the mass of sugarcane bagasse (g); and V is the volume of the solution (L).
The nonlinear optimization method was used to fit equilibrium and kinetics data to the corresponding models. For this study, non-linear regression was applied using 8.0 software for fitting the curve. The best fit of the chosen non-linear models was determined by the use of four well-known error functions, namely the coefficient of determination (R2), adjusted determination coefficient (adj-R2), reduced chi-square (red-χ2), and Bayesian information criterion (BIC). A model with a higher adj-R2 value, lower red-χ2 and BIC value indicates a better fitting than the others [25].
The equations of all error functions used are expressed as follows:
R 2 = 1 ( q i , exp q i , model ) 2 ( q i , exp q i , exp mean ) 2
adj R 2 = 1 ( 1 R 2 ) × ( N 1 DOF )
red χ 2 = ( q i , exp q i , model ) 2 DOF
BIC = N × ln ( ( q i , exp q i , model ) 2 N ) + Pln ( N )
where qi,exp is the experimental adsorbed capacity value; qi,model is the modeled value; qi,exp–mean is an average of qi,exp values used for modelling; N is the number of experimental points; P is the number of model parameters; and DOF is the degrees of freedom.

2.4. Experimental Design

Batch experiments based on a 25−1 fractional factorial design were conducted randomly to study the influence of the experimental variables on the percentage of COD removal (% Rem). The studied factors which are focused on this research; sugarcane bagasse dose (X1), solution pH (X2), contact time (X3), stirring speed (X4) and temperature (X5). Table 1 shows the five experimental variables and their chosen levels. After performing batch experiments for the optimization process, the regression analysis was conducted to attain the study’s statistical parameters with 95% confidence intervals using the Minitab 18 statistical software. The variance of the regression equation (mathematical relation between dependent and independent variables) is the most important analysis by the ANOVA method to find the desired function of COD adsorption. RSM is used as a sequential process to show the relation between the studied independent factors and the response to determine the set of optimal experimental parameters.

3. Results and Discussion

3.1. Analysis of Factorial Design

25−1 Fractional factorial design was adopted using Minitab 18 to optimize the influence of the investigated parameters; sugarcane bagasse dose (X1), solution pH (X2), contact time (X3), stirring speed (X4) and temperature (X5) on the elimination of COD from OMW. This design yields in 16 experiments with all possible combinations of X1, X2, X3, X4 and X5. COD removal efficiency (Y) was measured for each of these experiments as shown in Table 2. The response obtained was correlated using the second-order polynomial model, expressed by Equation (7):
Y = b0 + b1X1 + b2X2 + b3X3 + b4X4 + b5X5 + b12X1X2 + b13X1X3 + b14X1X4 + b15X1X5 + b23X2X + b24X2X4
+b34 X3X4 + b25X2X5 + b35X3X5 + b45X4X5
where Y is the COD removal efficiency response, b0 is a constant, bi correspond to linear coefficient of Xi, and bij is the interaction coefficient.
The statistical calculations and regression analysis were conducted to fit the response function with the experimental data. The regression coefficient values obtained are given in the final regression equation, after putting the values of all coefficients, as follows (Equation (8)):
Y (%) = 37.85 − 2.183X1 + 0.39X2 − 0.562X3 + 0.05958X4 + 0.2182X5 − 0.01712X1X3 − 0.0529X1X3 + 0.102X1X4 − 0.0528X1X5 − 0.05223X2X3 − 0.01976X2X4 − 0.01153X2X5 − 0.003092X3X4 − 0.00436X3X5 − 0.000623X4X5
This equation expresses the COD removal efficiency as a function of investigated experimental factors and enables fixing experimental conditions for each targeted COD removal efficiency. From Table 2, the values predicted obtained by the relation (8) are compared with those of experimental results, which show a good agreement between the two sets of values. In Equation (8), the positive signs of the coefficients for X2 and X4 factors and the X1X4 interaction imply their positive effect on the response, while the negative signs of the coefficients for X1, X3, and X5 factors as well as the X1X2, X1X3, X1X5, X2X3, X2X5, X3X4, X3X5, and X4X5 interactions represent the negative effect on the response. Based on the equation, the most influential factor for response was the sugarcane bagasse dose with a coefficient value of 2.183. The sign (-) means that each one-point decrease will have an effect of 2.183 on the COD removal efficiency value.
To check the quality of the model fitting, an ANOVA analysis was performed using the F-value and p-value (Table 3). Significant effects of the model terms were identified based on the p-values less than 0.05. Since the F-value for the 95% confidence interval, one degree of freedom, and 16 factorial runs, was equal to 4.54, all terms of the model with a F-value greater than 4.54 are considered statistically significant. The COD removal model was determined to be greatly significant from Fisher’s test (F-value of 22.24) and a smaller p-value (<0.0001). Besides that, the significance of the model terms can be determined through the Pareto chart (Figure 2). The Student’s t-test was performed to evaluate the significance of the regression coefficients. For the 95% confidence interval and one degree of freedom, it was observed that the t-value was equal to 2.12. Student’s t-test values for model terms are displayed in horizontal columns. All term effects are significant if their absolute values exceed the vertical line (2.12). According to the obtained F value, p-value and Pareto chart, it seems that the main effects X1, X2, X3, X4 and X5 factors as well as the X2X3, X1X4, X3X5, X1X2, X3X4, X4X5 and X1X3 interactions are statistically significant. In this way, the COD adsorption by sugarcane bagasse could be expressed using the following equation (Equation (9)):
Y (%) = 21.5 − 0.2439 X1 + 1.494 X2 + 1.387 X3 − 0.0451 X4 + 0.4404 X5 − 0.01712 X1X2
− 0.00497 X1X3 + 0.001892 X1X4 − 0.09223 X2X3 − 0.001131 X3X4 − 0.01268 X3X5 − 0.000743 X4X5
To graphically verify the validity of the regression model obtained for COD removal, four validation indicator plots were created by plotting the differences between the predicted (model) and the observed (experimental) values (Figure 3).
The normal probability plot (Figure 3a) reveals that all points are reasonably near a straight line, suggesting that the predicted values of COD removal and the actual experimental data were in agreement, evidencing the normal distribution of the data and the validity of the regression model. The graphic plot of the residuals (Figure 3b) displays the predicted and observed values. The residual (vertical axis) is the difference between the experimental and the fitted values [26]. As shown in Figure 3b, the experimental points are reasonably aligned around zero, suggesting normal distribution. The histogram (Figure 3c) shows a random distribution of values with no noticeable shape or trend across all 16 runs conducted. The residuals versus the observation orders (Figure 3d) show that the residuals seem to be randomly scattered around zero.

3.2. Response Surface Analysis

The significant interaction effect between each two input parameters on the % removal was presented by the 3D response surface and the contour of the plots (Figure 4A–G). The interaction between the solution pH (X2) and contact time (X3) is shown in Figure 4A. The graph demonstrated that better COD removal is achieved at acidic pH and long contact time. Further, the % removal decreased from 39 to 30% with increasing pH. This can be explained by the impact of pH on the ionic configuration of the functional groups presented in the sorbent. Previous pH drift tests indicated that the pHpzc (point of zero charge) of sugarcane bagasse is equal to 5.0, which expresses that the surface of the bagasse is positively charged at a pH below 5 and negatively charged at a pH above 5 [21,27]. As a consequence, the electrostatic interaction between the adsorbent surface and the organic matter could be enhanced, which resulted in high adsorption of COD [28]. The interaction effect between the adsorbent dose (X1) and agitation speed (X4) shown in Figure 4B reveals that there will be a good COD removal equal to 44% when the adsorbent dose is 10 g/L and also at the stirring speed of 80 rpm. Nevertheless, an adsorbent dose greater than 10 g/L and an agitation speed above 80 rpm result in a substantial decrease in % removal. This phenomenon is due to the aggregation and glomeration of the sorbent particles and the reduction in the total surface area of the adsorbent [29].
Figure 4C describes the interaction between the temperature (X5) and reaction time (X3) on the % removal. The % removal increases with an increase in temperature, whereas it decreases with increasing contact time. There are more active surface sites when the temperature increases because the adsorbent swells more as well. This result suggests that the COD adsorption is an endothermic process [30]. As shown in Figure 4D, the interaction between the pH (X2) and adsorbent dose (X1) has a negative effect on % removal. An increase in either of these two parameters reduces the COD removal efficiency. The possible reason for such observations is that the changes in pH could lead to changes in the properties of the surface of adsorbent [31], and the agglomeration of the adsorbate particles occurs with an increase in the adsorbent dose. In Figure 4E, the % removal decreased with increasing contact time (X3) and stirring speed (X4). The % removal was very rapid in the beginning stages of contact time, reaching about 38%. This result was related to the availability of more active sites on the adsorbent and the fact that the gradient of concentration between the adsorbate molecules in the solution and the adsorbate molecules on the adsorbent is high, which improves COD diffusion to the adsorbent surface [32,33]. After 60 min, the removal decreases over time due to adsorption site saturation, and the adsorbate molecules may bind poorly to the active receptors on the adsorbent [34,35].
Furthermore, the obtained results represented in Figure 4F show that the optimum of % removal was found at a temperature of 60 °C and agitation speed of 80 rpm, giving 39% of the recovery. This tendency may be due to the breakdown of the expanding chain and flocs. However, the higher stirring speed can encourage the process of agglomeration [36]. Additionally, the % removal was decreased by increasing the adsorbent dosage (X1) and contact time (X3) (Figure 4G). The % removal is initially higher and more rapid, reaching up to more than 38% within the first 60 min and 10 g/L, which could be attributed to the greater availability of adsorption sites to bind COD and a greater adsorbent active sites/COD ratio [37]. Subsequently, it decreased slightly with the prolonging of the contact time and the increase in adsorbent dose, which could be caused by an aggregation of the adsorbent and the reduction of the surface area available to COD [20].

3.3. Optimization Process

The response optimization study used to determine the combination of optimal values of the factors for maximum COD removal is shown in Figure 5. The validation of the optimal conditions is regarded as the final step in the modeling approach to examine the precision and robustness of the investigated model. The optimization requests to determine the desired pH, contact time, temperature, stirring speed, and adsorbent concentration to attain significant desirability. The model optimization shows that the efficiency of COD removal increases with an increasing pH, contact time and temperature, while decreasing with an increasing adsorbent dose and stirring speed. The optimum conditions indicate that the experimental modeling yields a desirability close to 1 with 55.07% COD removal efficiency, in the subsequent conditions: 10 g/L doses of sugarcane bagasse, 1 h of contact time, 80 rpm of stirring speed, 60 °C of temperature, and pH 12 of the OMW solution.

3.4. Characterization of Sugarcane Bagasse

Sugarcane bagasse characterization is an important analysis for understanding the behavior or the mechanism of COD removal on its surface. X-ray fluorescence elemental analysis was used to identify and quantify the amount of the elements contained in sugarcane bagasse in order to ultimately determine its elemental composition. The mineral composition of sugarcane bagasse is presented in Table 4. The results demonstrate that sugarcane bagasse contains a high amount of silicon dioxide (SiO2), as high as 62.23%, and low amounts of alkaline oxide, Al2O3 and P2O5.
Therefore, the FTIR was utilized to investigate the functional groups of the sugarcane bagasse (Figure 6). The bands at 3330 and 2890 cm−1 indicate the presence of the O-H functional group and C–H stretching, respectively [20,38]. The peak at 1632 cm−1 is assigned to C=O vibrations in hemicellulose [39]. These groups are thought to play a very important role in the process of adsorption [40]. The bands that appeared at wave numbers between 1471 cm−1 and 1366 cm−1 are related to C=C–H indicating several bands in cellulose and xylose [41]. The bands appeared at 1241 cm−1 and 1029 cm−1 can be attributed to the CH=CH stretching of lignin [39] and C-O stretching in cellulose and hemicellulose [42], respectively. The bands appeared at 640 and 593 cm−1 in the FTIR spectra are attributed to the vibration of O-H groups out of the plane deformation [43].

3.5. Adsorption Isotherms

The specific relationship established between the COD remaining in solution (Ce) and the COD adsorbed by sugarcane bagasse (qe) at equilibrium was analyzed by the models of Langmuir and Freundlich. The investigated adsorption isotherms for the COD adsorption on sugarcane bagasse are presented in Figure 7. Table 5 shows the calculated parameters for each of the sorption isotherm models obtained from non-linear regression forms. The best fit of the experimental data was analyzed based on R2, adj-R2, red-χ2 and BIC. The obtained results displayed that the COD adsorption fitted well with the Langmuir isotherm model with the highest adj-R2 of 0.994 and the lowest red-χ2 and BIC values of 78.33 and 25.95, respectively. This result indicates that the adsorption process occurs on a homogeneous surface, and all adsorption sites are identical and energetically equivalent [44]. The maximum adsorption capacity (qmax) of COD onto sugarcane bagasse was determined to be 331.92 mg/g from the Langmuir model.
The adsorption properties of sugarcane bagasse are compared with other adsorbents used for COD sorption reported in the literature, which are given in Table 6. The table indicates that the qmax value obtained in the present study is higher than that of other materials from previous studies. The qmax value found in this study reveals a very good adsorption capacity of the sugarcane bagasse, which falls as a promising adsorbent.

3.6. Kinetic Studies

Kinetic modeling was undertaken to determine the rate of COD adsorption on the sugarcane bagasse and examine the controlling mechanisms of the adsorption process. The kinetic studies were carried out from 0 to 24 h, and the recorded data were studied with the pseudo-first-order and pseudo-second-order kinetic models. The rate constants values were valued from the non-linear plots and shown in Figure 8 and are also summarized in Table 7. The table indicated that the R2 value obtained from the pseudo-first order kinetic equation was found to be higher than that of pseudo-second order. In addition, the calculated qe value (337.45 mg/g) for the pseudo-first-order model is closer to the experimental value (326.29 mg/g) compared to the qe value (405.9 mg/g) calculated for the pseudo-second-order model. The pseudo-first-order model also presented the highest adj-R2 (0.967), lowest red-χ2 (467.41) and BIC (43.01) values. This result advises that the experimental data demonstrated the best fit to the pseudo-first order model, and its applicability also indicates that a physical process might control the sorption process.

4. Conclusions

In this work, the sugarcane bagasse prepared was utilized as an adsorbent for the elimination of COD from OMW by the adsorption. Sugarcane bagasse characterization was performed by FTIR and X-ray fluorescence. The 21−5 fractional design associated to response surface methodology was successfully applied to optimize the effects of the operating variables of COD removal from OMW using the prepared adsorbent. The optimal conditions were found to be pH 12, the adsorbent dose of 10 g/L, a stirring speed of 80 rpm, a contact time of 1 h, and a temperature of 60 °C with a percentage of removal of 55.07% and a desirability close to 1. The experimental data are well correlated by the Langmuir isotherm model with an R2 of 0.996, and the maximum sorption capacity of 331.92 mg/g under optimal conditions. The kinetic data of the COD removal process were properly fitted with the pseudo-first-order model instead of pseudo-second-order model.

Author Contributions

F.E. carried out the experiments and prepared the draft manuscript. S.E.: and W.B. interpreted and discussed the results, and calculated the statistical parameters. N.B. (Nadia Beniich) analyzed and checked the statistical results. N.B. (Noureddine Barka) verified the analytical methods. C.E.A. and E.L. aided in interpreting the results and worked on the manuscript. M.A. corrected and wrote the final version of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bayramoglu, M.; Kobya, M.; Can, O.T.; Sozbir, M. Operating cost analysis of electrocoagulation of textile dye wastewater. Sep. Purif. Technol. 2004, 37, 117–125. [Google Scholar] [CrossRef]
  2. Sayan, E. Optimization and modeling of decolorization and COD reduction of reactive dye solutions by ultrasound-assisted adsorption. Chem. Eng. J. 2006, 119, 175–181. [Google Scholar] [CrossRef]
  3. Rupani, P.F.; Singh, R.P.; Ibrahim, M.H.; Esa, N. Review of current palm oil mill effluent (POME) treatment methods: Vermicomposting as a sustainable practice. World Appl. Sci. J. 2010, 11, 70–81. [Google Scholar]
  4. Oladipo, A.A.; Adeleye, O.J.; Oladipo, A.S.; Aleshinloye, A.O. Bio-derived MgO nanopowders for BOD and COD reduction from tannery wastewater. J. Water Process Eng. 2017, 16, 142–148. [Google Scholar] [CrossRef]
  5. Elayadi, F.; El Adlouni, C.; Achak, M.; El Herradi, E.; El Krati, M.; Tahiri, S.; Naman, M.; Naman, F. Effects of raw and treated olive mill wastewater (OMW) by coagulation-flocculation, on the germination and the growth of three plant species (wheat, white beans, lettuce). Mor. J. Chem. 2019, 7, 111–122. [Google Scholar]
  6. Khatamian, M.; Divband, B.; Shahi, R. Ultrasound assisted co-precipitation synthesis of Fe3O4/bentonite nanocomposite: Performance for nitrate, BOD and COD water treatment. J. Water Process Eng. 2019, 3, 100870. [Google Scholar] [CrossRef]
  7. Aquilanti, L.; Taccari, M.; Bruglieri, D.; Osimani, A.; Clementi, F.; Comitini, F.; Ciani, M. Integrated biological approaches for olive mill wastewater treatment and agricultural exploitation. Int. Biodeterior. Biodegrad. 2014, 88, 162–168. [Google Scholar] [CrossRef]
  8. Mastoras, P.; Vakalis, S.; Fountoulakis, M.S.; Gatidou, G.; Katsianou, P.; Koulis, G.; Thomaidis, N.S.; Haralambopoulos, D.; Stasinakis, A.S. Evaluation of the performance of a pilot-scale solar still for olive mill wastewater treatment. J. Clean. Prod. 2022, 365, 132695. [Google Scholar] [CrossRef]
  9. Esteves, B.M.; Morales-Torres, S.; Madeira, L.M.; Maldonado-Hodar, F.J. Specific adsorbents for the treatment of OMW phenolic compounds by activation of bio-residues from the olive oil industry. J. Environ. Manag. 2022, 306, 114490. [Google Scholar] [CrossRef]
  10. Alaoui, S.B.; Lamy, E.; Achak, M. Assessment of the impact of diluted and pretreated olive mill wastewater on the treatment efficiency by infiltration-percolation using natural bio-adsorbents. Environ. Sci. Pollut. Res. 2022, 30, 16305–16320. [Google Scholar] [CrossRef]
  11. Jawad, A.H.; Abdulhameed, A.S.; Mastuli, M.S. Acid-factionalized biomass material for methylene blue dye removal: A comprehensive adsorption and mechanism study. J. Taibah Univ. Sci. 2020, 14, 305–313. [Google Scholar] [CrossRef]
  12. Mohan, D.; Pittman, C.U., Jr. Activated carbons and low cost adsorbents for remediation of tri- and hexavalent chromium from water. J. Hazard. Mater. 2006, 137, 762–811. [Google Scholar] [CrossRef] [PubMed]
  13. Hameed, B.H. Spent tea leaves: A new non-conventional and low-cost adsorbent for removal of basic dye from aqueous solutions. J. Hazard. Mater. 2009, 161, 753–759. [Google Scholar] [CrossRef]
  14. El-Naas, M.H.; Al-Zuhair, S.; Abu Alhaija, M. Reduction of COD in refinery wastewater through adsorption on date-pit activated carbon. J. Hazard. Mater. 2010, 173, 750–757. [Google Scholar] [CrossRef]
  15. Low, L.W.; Teng, T.T.; Alkarkhi, A.F.M.; Ahmad, A.; Morad, N. Optimization of the Adsorption Conditions for the Decolorization and COD Reduction of Methylene Blue Aqueous Solution using Low-Cost Adsorbent. Water Air Soil Pollut. 2011, 214, 185–195. [Google Scholar] [CrossRef]
  16. Ismail, M.N.; Aziz, H.A.; Ahmad, M.A.; Yusoff, N.A. Optimization of Areca cathecu Fronds as Adsorbent for Decolorization and COD Removal of Wastewater through the Adsorption Process. Sains Malaysiana. 2015, 11, 1609–1614. [Google Scholar]
  17. Nayl, A.E.A.; Elkhashab, R.A.; El Malah, T.; Sobhy, M.; Yakout Mohamed, A.; El-Khateeb Mahmoud, M.S.; Ali Hazim, M.A. Adsorption studies on the removal of COD and BOD from treated sewage using activated carbon prepared from date palm waste. Environ. Sci. Pollut. Res. 2017, 24, 22284–22293. [Google Scholar] [CrossRef]
  18. Mpatani, F.M.; Aryee, A.A.; Kani, A.N.; Wen, K.; Dovi, E.; Qu, L.; Li, Z.; Han, R. Removal of methylene blue from aqueous medium by citrate modified bagasse: Kinetic, Equilibrium and Thermodynamic study. Bioresour. Technol. Rep. 2020, 11, 100463. [Google Scholar] [CrossRef]
  19. Ge, M.; Du, M.; Zheng, L.; Wang, B.; Zhou, X.; Jia, Z.; Hu, G.; Jahangir Alam, S.M. A maleic anhydride grafted sugarcane bagasse adsorbent and its performance on the removal of methylene blue from related wastewater. Mater. Chem. Phys. 2017, 192, 147–155. [Google Scholar] [CrossRef]
  20. Surafel, M.B.; Venkatesa, P.S.; Tsegaye, S.T.; Abraham, G.A. Sugarcane bagasse based activated carbon preparation and its adsorption efficacy on removal of BOD and COD from textile effluents: RSM based modeling, optimization and kinetic aspects. Bioresour. Technol. Rep. 2021, 14, 100664. [Google Scholar] [CrossRef]
  21. Zhang, Z.; Moghaddam, L.; O’Hara, I.M.; Doherty, W.O.S. Congo Red adsorption by ball-milled sugarcane bagasse. Chem. Eng. J. 2011, 178, 122–128. [Google Scholar] [CrossRef]
  22. Akl, M.A.A.; Dawy, M.B.; Serage, A.A. Efficient removal of phenol from water samples using sugarcane bagasse based activated carbon. J. Anal Bioanal Tech. 2014, 5, 189. [Google Scholar] [CrossRef]
  23. Abdelhafez, A.; Jianhua, L. Removal of Pb (II) from aqueous solution by using biochars derived from sugar cane bagasse and orange peel. J. Taiwan Inst. Chem. Eng. 2016, 61, 367–375. [Google Scholar] [CrossRef]
  24. Elayadi, F.; Boumya, W.; Achak, M.; Chhiti, Y.; M’hamdi Alaoui, F.E.; Barka, N.; El Adlouni, C. Experimental and modelling studies of the removal of phenolic compounds from olive mill wastewater by adsorption on sugarcane bagasse. Environ. Chall. 2021, 4, 100184. [Google Scholar] [CrossRef]
  25. Tran, H.N. Applying Linear Forms of Pseudo-Second-Order Kinetic Model for Feasibly Identifying Errors in the Initial Periods of Time-Dependent Adsorption Datasets. Water 2023, 15, 1231. [Google Scholar] [CrossRef]
  26. Tawfik, A.; Saleh, S.; Omobayo, A.; Asif, M.; Dafalla, H. Response Surface Optimization and Statistical Analysis of Phenols Adsorption On Diethylenetriamine modified Activated Carbon. J. Clean. Prod. 2018, 182, 960–968. [Google Scholar]
  27. Moubarik, A.; Grimi, N. Valorization of olive stone and sugar cane bagasse by-products as biosorbents for the removal of cadmium from aqueous solution. Food Res. Int. 2015, 73, 169–175. [Google Scholar] [CrossRef]
  28. de Sales, P.F.; Magriotis, Z.M.; Rossi, M.A.L.S.; Resende, R.F.; Nunes, C.R. Optimization by Response Surface Methodology of the adsorption of Coomassie Blue dye on natural and acid-treated clays. J. Environ. Manag. 2013, 130, 417e428. [Google Scholar] [CrossRef]
  29. Leili, M.; Faradmal, J.; Kosravian, F.; Heydari, M. A Comparison study on the removal of phenol from aqueous solution using Organomodified Bentonite and commercial activated carbon. Avicenna J. Environ. Health Eng. 2015, 2, 2698. [Google Scholar] [CrossRef]
  30. Mondal, N.K.; Roy, S. Optimization study of adsorption parameters for removal of phenol on gastropod shell dust using response surface methodology. Clean Technol. Environ. Policy 2016, 18, 429–447. [Google Scholar] [CrossRef]
  31. Maleki, S.; Karimi-Jashni, A. Optimization of Ni(II) adsorption onto Cloisite Na+ clay using response surface2 methodology. Chemosphere 2019, 246, 125710. [Google Scholar] [CrossRef] [PubMed]
  32. Olusegun, S.J.; de Sousa Lima, L.F.; Santina Mohallem, N.D. Enhancement of adsorption capacity of clay through spray drying and surface modification process for wastewater treatment. Chem. Eng. J. 2018, 334, 1719–1728. [Google Scholar] [CrossRef]
  33. Hafshejani, L.D.; Hooshmand, A.; Naseri, A.A.; Mohammadi, A.S.; Abbasi, F.; Bhatnagar, A. Removal of nitrate from aqueous solution by modified sugarcanebagasse biochar. Ecol. Eng. 2016, 95, 101–111. [Google Scholar] [CrossRef]
  34. Yangui, A.; Abderrabba, M. Towards a high yield recovery of polyphenols from olive mill wastewater on activated carbon coated with milk proteins: Experimental design and antioxidant activity. J. Food Chem. 2018, 262, 102–109. [Google Scholar] [CrossRef] [PubMed]
  35. Mandal, A.; Bar, N.; Kumar Das, S. Phenol removal from wastewater using low-cost natural bioadsorbent neem (Azadirachta indica) leaves: Adsorption study and MLR modeling. Sustain. Chem. Pharm. 2020, 17, 100–308. [Google Scholar] [CrossRef]
  36. Geetha Devi, M.; Shinoon Al-Hashmi, Z.S.; Chandra Sekhar, G. Treatment of vegetable oil mill effluent using crab shell chitosan as adsorbent. Int. J. Environ. Sci. Technol. 2012, 9, 713–718. [Google Scholar] [CrossRef]
  37. Pourmortazavi, S.M.; Sahebi, H.; Zandavar, H.; Mirsadeghi, S. Fabrication of Fe3O4 nanoparticles coated by extracted shrimp peels chitosan as sustainable adsorbents for removal of chromium contaminates from wastewater: The design of experiment Compos B. J. Eng. 2019, 175, 107130. [Google Scholar] [CrossRef]
  38. Palin, D., Jr.; Rufato, K.B.; Linde, G.A.; Colauto, N.B.; Caetano, J.; Alberton, O.; Jesus, D.A.; Dragunski, D.C. Evaluation of Pb (II) biosorption utilizing sugarcane bagasse colonized by Basidiomycetes. Environ. Monit Assess. 2016, 188, 188–279. [Google Scholar] [CrossRef]
  39. Lei, S.; Dongmei, C.; Shungang, W.; Zebin, Y. Adsorption Studies of Dimetridazole and Metronidazole onto Biochar Derived from Sugarcane Bagasse: Kinetic, Equilibrium, and Mechanisms. J. Polym. Environ. 2017, 26, 765–777. [Google Scholar] [CrossRef]
  40. Srinivasan, A.; Viraraghavan, T. Decolorization of dye wastewaters by biosor-bents: A review. J. Environ. Manag. 2010, 91, 1915–1929. [Google Scholar] [CrossRef]
  41. Sharma, P.; Kaur, H. Sugarcane bagasse for the removal of erythrosin B and methylene blue from aqueous waste. Appl. Water Sci. 2011, 1, 135–145. [Google Scholar] [CrossRef]
  42. Kumar, R.; Mago, G.; Balan, V.; Wyman, C.E. Physical and chemical character-izations of corn stover and poplar solids resulting from leading pretreatment technologies. Bioresour. Technol. 2009, 100, 3948–3962. [Google Scholar] [CrossRef] [PubMed]
  43. Granados, D.A.; Ruiz, R.A.; Vega, L.Y.; Chejne, F. Study of reactivity reduction in sugarcane bagasse as consequence of a torrefaction process. Energy 2017, 139, 818–827. [Google Scholar] [CrossRef]
  44. Vimonses, V.; Lei, S.M.; Jin, B.; Chowd, C.W.K.; Saint, C. Kinetic study and equilibrium isotherm analysis of Congo Red adsorption by clay materials. J. Chem. Eng. 2009, 148, 354–364. [Google Scholar] [CrossRef]
  45. Ahmad, A.A.; Hameed, B.H. Reduction of COD and color of dyeing effluent from a cotton textile mill by adsorption onto bamboo-based activated carbon. J. Hazard. Mater. 2009, 172, 1538–1543. [Google Scholar] [CrossRef]
  46. Akter, S.; Sultana, F.; Kabir, M.R.; Brahma, P.P.; Tasneem, A.; Sarker, N.; Mst, M.; Akter, M.; Kabir, M.; Uddin, M.K. Comparative Study of COD Removal Efficacy from Pharmaceutical Wastewater by Areca Nut Husk Produced and Commercially Available Activated Carbons. Asian J. Environ. Ecol. 2021, 14, 47–56. [Google Scholar] [CrossRef]
Figure 1. Sugarcane bagasse material steps.
Figure 1. Sugarcane bagasse material steps.
Water 15 01630 g001
Figure 2. Pareto chart for standardized effects.
Figure 2. Pareto chart for standardized effects.
Water 15 01630 g002
Figure 3. Validation indicator plots of the model for the experiments (a) Normal probability plot, (b) Versus fits (c) Histogram and (d) Versus order.
Figure 3. Validation indicator plots of the model for the experiments (a) Normal probability plot, (b) Versus fits (c) Histogram and (d) Versus order.
Water 15 01630 g003
Figure 4. (AG): Response surface and contour plots for the % Rem.
Figure 4. (AG): Response surface and contour plots for the % Rem.
Water 15 01630 g004aWater 15 01630 g004b
Figure 5. Optimization of the model obtained by desirability function.
Figure 5. Optimization of the model obtained by desirability function.
Water 15 01630 g005
Figure 6. The FTIR spectra of sugarcane bagasse.
Figure 6. The FTIR spectra of sugarcane bagasse.
Water 15 01630 g006
Figure 7. Isotherm plots for COD removal by sugarcane bagasse.
Figure 7. Isotherm plots for COD removal by sugarcane bagasse.
Water 15 01630 g007
Figure 8. Kinetic plots for COD removal by bagasse.
Figure 8. Kinetic plots for COD removal by bagasse.
Water 15 01630 g008
Table 1. Process factors and their levels.
Table 1. Process factors and their levels.
FactorsVariableUnitLevels
Low (−)High (+)
X1Adsorbent dose g/L1060
X2pH-212
X3Contact timeh124
X4Stirring speedrpm80300
X5Temperature°C2560
Table 2. Matrix design with experimental and predicted COD removal values.
Table 2. Matrix design with experimental and predicted COD removal values.
X1X2X3X4X5Y (Experimental)Y (Predicted)
1−1−1−1−1+136.3937.46
21−1−1−1−141.2643.47
3−11−1−1−154.2751.69
411−1−1+144.6943.75
5−1−11−1−135.2435.54
61−11−1+142.4043.24
7−111−1+135.5330.66
8111−1−141.2340.85
9−1−1−11−140.7140.53
101−1−11+142.1243.90
11−11−11+133.7332.73
1211−11−139.8239.12
13−1−111+133.2329.86
141−111−131.1234.75
15−1111−139.2538.14
161111+142.1241.79
Table 3. Analysis of ANOVA for COD removal. The underline shows the p-values < 0.05.
Table 3. Analysis of ANOVA for COD removal. The underline shows the p-values < 0.05.
SourceDFAdj SSAdj MSF-Valuep-Value
Model163214.58200.91222.240.000
Linear5881.04176.20719.510.000
X11224.88224.87624.900.000
X21192.57192.56721.320.000
X31146.82146.82316.250.001
X41260.52260.52528.840.000
X5156.2556.2456.230.025
Interactions102327.91232.79125.770.000
X1 × X21146.54146.53816.220.001
X1 × X3165.3565.3507.230.017
X1 × X41865.84865.83995.860.000
X1 × X515.315.3140.590.455
X2 × X31900.00899.99899.640.000
X2 × X410.340.3430.040.848
X2 × X515.265.2600.580.457
X3 × X4165.5065.5027.250.017
X3 × X51208.31208.30623.060.000
X4 × X5165.4665.4647.250.017
Error15135.499.033
Total sum of squares313350.07
S squareR2R2 (ajust)R2 (prev)
6.253899.98%87.54%62.23%
Table 4. X-ray fluorescence analysis of sugarcane bagasse.
Table 4. X-ray fluorescence analysis of sugarcane bagasse.
CompositionConcentration (%)
SiO262.23
SO324.00
CaO11.50
MgO0.68
Na2O0.60
K2O0.51
Al2O30.28
P2O50.1
Fe2O30.10
Table 5. Parameters and error function data for sorption isotherm models obtained from non-linear regression forms.
Table 5. Parameters and error function data for sorption isotherm models obtained from non-linear regression forms.
ModelsEquation aParametersValues
Langmuir q e = q m a x K L C e ( 1 + K L C e ) qm (mg/g)331.92
KL0.019
R20.996
Adj-R20.994
Red-χ278.33
BIC25.59
Freundlich q e = K F C e 1 / n 1/n0.31
KF51.13
R20.991
Adj-R20.985
χ2208.11
BIC31.45
Note: a qe (mg/g): mass of adsorbed molecule per unit mass of sugarcane bagasse, Ce (mg/L): concentration of no-adsorbed molecules, qm and KL are constants of the Langmuir model, and KF and n are constants of Freundlich isotherm model.
Table 6. Comparison of the qmax of COD removal by various adsorbents.
Table 6. Comparison of the qmax of COD removal by various adsorbents.
Adsorbent Typesqmax (mg/g)Ref.
Sugarcane bagasse331.92 This study
Date-pit activated carbon252.81[14]
Activated carbon41.3[17]
Raw bagasse77.95[15]
Bamboo-based activated carbon24.39 [45]
Areca Nut Husk64.94[46]
Table 7. Parameters and error functions data for kinetic models studied obtained from non-linear regression forms.
Table 7. Parameters and error functions data for kinetic models studied obtained from non-linear regression forms.
ModelEquation bParametersValue
Pseudo-first-order q t = q e ( 1 e k 1 t ) qe.exp (mg/g)326.92
qecal (mg/g)337.45
k10.20
R120.978
Adj-R20.967
χ2467.41
BIC43.01
Pseudo-second-order q t = K 2 q e 2 t ( 1 + K 2 q e t ) qecal (mg/g)405.9
k2 0.0006
R220.976
Adj-R20.964
χ2503.52
BIC43.52
Note: b qe and qt (mg/g) are the amounts of COD removed at equilibrium and at time t, respectively; k1 and k2 are the adsorption rate constants.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Elayadi, F.; Achak, M.; Boumya, W.; Elamraoui, S.; Barka, N.; Lamy, E.; Beniich, N.; El Adlouni, C. Factorial Design Statistical Analysis and Optimization of the Adsorptive Removal of COD from Olive Mill Wastewater Using Sugarcane Bagasse as a Low-Cost Adsorbent. Water 2023, 15, 1630. https://doi.org/10.3390/w15081630

AMA Style

Elayadi F, Achak M, Boumya W, Elamraoui S, Barka N, Lamy E, Beniich N, El Adlouni C. Factorial Design Statistical Analysis and Optimization of the Adsorptive Removal of COD from Olive Mill Wastewater Using Sugarcane Bagasse as a Low-Cost Adsorbent. Water. 2023; 15(8):1630. https://doi.org/10.3390/w15081630

Chicago/Turabian Style

Elayadi, Fatima, Mounia Achak, Wafaa Boumya, Sabah Elamraoui, Noureddine Barka, Edvina Lamy, Nadia Beniich, and Chakib El Adlouni. 2023. "Factorial Design Statistical Analysis and Optimization of the Adsorptive Removal of COD from Olive Mill Wastewater Using Sugarcane Bagasse as a Low-Cost Adsorbent" Water 15, no. 8: 1630. https://doi.org/10.3390/w15081630

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop