This Article


Creative Commons License
Except where otherwise noted, this work is licensed under Creative Commons Attribution-NonCommercial 4.0 International License.

Application of Heterogeneous Catalytic Ozonation Process with Magnesium Oxide Nanoparticles for Toluene Degradation in Aqueous Environments

1 Health Promotion Research Center, Zahedan University of Medical Sciences, Zahedan, IR Iran
2 Analytical Research Laboratory, Department of Chemistry, University of Sistan and Baluchestan, Zahedan, IR Iran
*Corresponding author: Edris Bazrafshan, Health Promotion Research Center, Zahedan University of Medical Sciences, Zahedan, IR Iran, E-mail:
Health Scope. 5(4): e40439 , DOI: 10.17795/jhealthscope-40439
Article Type: Research Article; Received: Jul 1, 2016; Revised: Aug 25, 2016; Accepted: Sep 16, 2016; epub: Sep 27, 2016; collection: Nov 2016


Background: Volatile organic compounds (VOCs) removal like toluene is an important concern for health of humans and the environment.

Objectives: In this study, the efficiency of heterogeneous catalytic ozonation process (COP) with MgO nanoparticles for toluene degradation in aqueous solutions was investigated.

Materials and Methods: Response surface methodology (RSM) using central composite design (CCD) was applied to study the interactions between critical operating factors and to find their optimum levels. Four independent factors namely pH, catalyst dosage, ozonation time and concentration of toluene were transformed to coded values and consequently a second order quadratic model was prepared for responses prediction.

Results: The analysis of variance (ANOVA) of the toluene degradation by COP/MgO nanoparticles showed F-value of 19.53 for the quadratic model indicating that the prepared model is significant. The predicted result showed that maximum degradation of toluene (~99.99%) could be achieved at pH 12, MgO dose 0.5 g, ozonation time of 50 minutes and initial toluene concentration of 10 mg/l. The correlation coefficient (R2 = 0.9480) indicates a good correlation between the experimental and predicted values.

Conclusions: According to the results, COP is an efficient and rapid method for removing of toluene from aqueous solutions.

Keywords: Toluene degradation; Heterogeneous Catalytic Ozonation; MgO/O3

1. Introduction

Volatile organic compounds (VOCs) are normally considered as toxic chemicals that cause several health problems. Between the common VOCs, aromatic VOCs embody an important group of hazardous contaminants. Most of these pollutants are considered malodorous, poisonous, possibly mutagenic, and also confirmed carcinogenic to human (1, 2).

Monoaromatic compounds such as Benzene, Toluene, Ethylbenzene and Xylenes (BTEX compounds) were recognized as common VOCs and are very important (3). Furthermore, BTEX compounds constitute aromatic components in many petroleum products that are commonly found in groundwater sources. The presence of BTEX compounds in groundwater sources is usually caused by leaking underground storage tanks and pipelines, improper disposal of waste and leakage from landfill sites (4-6). Additionally, because of great solubility in aqueous environments comparative to other common hydrocarbons especially petroleum hydrocarbons, BTEX combinations can spread in water resources distant from the sources of contamination (7).

Eicosapentaenoic (EPA) has been categorized as BTEX compounds, priority contaminants that incline to cause significant threat to the environment and also health of human because of their poisonous and carcinogenic features. Even at minimum concentrations, BTEX compounds can cause kidney and liver damage (4, 5). As presented by the world health organization (WHO), the maximum allowable concentration of benzene, toluene, ethylbenzene and xylene in drinking water is 0.01, 0.7, 0.3 and 0.5 mg/L, respectively (8). In addition, the BTEX-contaminated water bodies particularly groundwater sources indicate a severe pollution issue due to great carcinogenetic affects (9, 10). Thus, the presence of these contaminants in water sources is objectionable, and removal of these compounds with reasonable approaches is necessary.

Among the BTEX components, toluene (or methylbenzene) is a typical VOC, which is harmful to humans (11). In particular, Toluene as a member of this group, is used as one of the most common solvents to dissolve paints, paint thinners, rubber and adhesives and causes numerous adverse effects on health of humans and the environment (12).

During the past few years, considerable efforts have been made on the development of Advanced Oxidation Processes (AOPs) for water and wastewater treatment, particularly for the treatment of refractory organic compounds and recalcitrant compounds (13). Based on previous studies, AOPs are quite an effective technique for pollutants decomposition because of the hydroxyl radical (•OH) generation as the predominant species (14). Advanced Oxidation Processes are an economic process, quite effective and easy to use (15).

Due to easy recycling (16), lower cost (17), greater efficiency of pollutants degradation (18), and ease of operation, catalytic ozonation process (COP) is considered as a novel AOP for water and wastewater treatment (19, 20). Indeed, COP can be regarded as a favorable AOP for purification of the environment. In this method, various catalysts decompose ozone molecules into hydroxyl radicals (•OH), which will react with organic materials in aqueous solutions in a non-selective manner (21, 22). On the other hand, COP promotes the hydroxyl radicals (HO•) generation from decomposition of ozone molecules, consequently overcoming the ozonation process limitations, including byproducts production and also inconsiderable efficiency for ozone consumption (23). It is well known that catalyst addition to the ozonation process increases the rate of oxidation, decreases the contact time and finally reduces the cost of treatment (24). A variety of metal oxides, including rare oxides, such as MnO2, CeO2, Ni2O3, TiO2, MnO2, Al2O3, CuO, V2O5, Co3O4, La2O3, Pr6O11, TiO2, Al2O3, ZnO, MgO, etc. have been selected and used as catalysts in COP (15, 25-27).

Various parameters such as pH of the solution, temperature, nature of pollutants, dose and type of catalyst, concentration of ozone molecules and also shape of the reactor can affect the performance of COP (27). Catalytic ozonation has been classified as homogenous and heterogeneous (24, 28, 29). Each of these methods by direct oxidation (reaction of organic materials with ozone molecules) or indirect oxidation (via hydroxyl radicals generation) can stabilize different organic compounds (16, 24, 30). Heterogeneous COP is a kind of innovative AOPs, which has the capability to improve the ozonation of recalcitrant contaminants in the presence of various catalysts (31, 32).

Furthermore, heterogeneous COP is an interesting method for treatment of different wastewaters because of its minimum cost, possible recovery of solid catalyst, and also minimum generation of secondary contamination. In fact, this process combines the molecules of ozone with the oxidative and adsorptive properties of metal oxide catalysts to attain mineralization of different organic pollutants at ambient temperature (33). Catalytic mechanisms are adopted for the COP, including acceleration of the production of hydroxyl radicals, complex catalysis, and the sorption of ozone molecules and organic matter on the surface of catalyst (15, 26).

Magnesium Oxide (MgO), as a heterogeneous catalyst, has very good potential for destruction of different organic contaminants (34, 35). It is an economical, nontoxic and also environment-friendly compound. MgO has already been widely used for the treatment of wastewater (36). For example, catalytic ozonation using MgO nanoparticles have been used efficiently for removal of dye, phenol and formaldehyde (19).

2. Objectives

In the current study, the degradation of toluene using heterogeneous catalytic ozonation with MgO nanoparticles was investigated in a batch reactor. The toluene degrading performance was checked by determining changes in the degradation efficiency when the reactor was operated under various conditions. In addition, for this study, to design the experiments, analysis of data and optimal condition determination, Central composite design (CCD) and response surface methodology (RSM) were used. Several important factors including pH, catalyst dose (MgO), time of ozonation and initial toluene concentration were converted to coded values (as variables) and consequently second order quadratic model was prepared for response prediction (toluene degradation).

3. Materials and Methods

3.1. Materials

Toluene (C7H8, 99.5%) and methanol (CH3-OH, 99.7%)) used in this work were procured from Merck. MgO nanoparticles were purchased from the Sigma-Aldrich company. The entire chemicals were used as received without extra pretreatments. For removal of any impurities, glassware were washed with chromic acid and rinsed with distilled water. Nanosized MgO was < 50 nm and a BET surface area analyzer was used for determination of surface area of MgO nanoparticles (120 m2/g).

3.2. Experimental Apparatus and Procedure

Catalytic ozonation was performed in a cylindrical glass column (Figure 1). The experiments were done at room temperature (22 ± 3°C) under batch mode, with 0.5 L of synthetic solution (with predefined concentration of toluene). Several critical parameters such as solution pH, MgO nanoparticles dosage, ozonation time and initial concentration of toluene were investigated.

Figure 1.
Schematic of the Used Catalytic Ozonation Reactor

In the present study, MgO nanoparticles were used as catalyst for oxidation of toluene in the ozonation process. Furthermore, efficiency of single ozonation and adsorption by MgO nanoparticles for toluene removal were determined. The ozone gas was generated from cylinder pure oxygen by the ozone generator (Pasargad model). Samples after predetermined reaction times were analyzed for residual toluene concentration. The degradation rate of toluene can be calculated via Equation 1, where Ct is toluene concentration after several regular intervals of time (t) and Ct = 0 is the initial concentration of toluene.

Equation 1.

3.3. Analytical Methods

In this study, toluene concentration as pollutant was analyzed using a headspace-gas chromatograph (Agilent 7890 A, Palo Alto, CA, USA) equipped with a flame ionization detector (FID), and an HP-5 capillary column (30-m length, 0.25-µm thickness, 0.32-mm I. D.). The pH of the solution was determined by a pH meter (WTW inoLab 7310).

3.4. Design of Experiments and Statistical Analysis

The RSM is a collection of mathematical and statistical methods designed for the modeling and analysis of problems wherein a response of interest is influenced by numerous affecting factors (37, 38). Determination of the optimal work settings of the system or determination of a region that satisfies the operational conditions are the key goal of RSM (38). The RSM has been used in various branches of science. For example, Gonen and Aksu (37), Anupam et al. (39) and Jo et al. (5), have efficiently used RSM for optimization of several factors for degradation of different pollutants.

In the present study, four factors and three levels of full factorial CCD based on RSM was used and in total 30 experiments were done (including eight axial points, 16 factorial points and six replicates at the center points) to optimize the chosen variables, including pH, MgO nanoparticles dosage, time of ozonation and initial concentration of toluene. Based on the initial analysis, the range and levels used in the experiments are designated and scheduled in Table 1. In addition, the experiment settings, which were designed by using central composite experimental design (CCD) for survey of the effect of several factors and optimization of process parameters, are illustrated in Table 2.

Table 1.
Experimental Range and Level of the Independent Variables

Table 2.
Experimental Design Matrix for Toluene Degradation

Each response of Y can be represented by a mathematical equation that correlates the response surface. The response (Y) can be expressed as a polynomial model based on the following quadratic equation (Equation 2):

Equation 2.

Where Y represents the predicted response (the toluene removal by the catalytic ozonation); β0 and βi are the constant coefficient and linear coefficient, βii and βij are the quadratic coefficients and the interaction coefficients, and Xi and Xj are the coded values of the independent process factors. This equation states the correlation between the independent variables in coded values and predicted response based on Tables 1 and 2. Furthermore, analysis of variance (ANOVA) was used for evaluation of the significance of quadratic regression model. Also, assessment of model terms was performed by the P value with 95% confidence level. Coefficient parameters were assessed by response surface regression analysis using the software Design-Expert (version 7) and Minitab (version 17); also applied to determine the residuals, 3-D surface and 2-D contour plots of the response models.

4. Results

4.1. Model Fitting and Statistical Analysis

The optimum conditions for degradation of toluene under catalytic ozonation process (COP) were determined by means of the CCD under RSM. In this way, the toluene removal by heterogeneous COP, adsorption and SOP could be expressed using the following equation:

Y1 (toluene removal by COP) = 68.30 +3.96 X1 + 5.11 X2 + 2.06 X3 - 12.11 X4 + 1.10 X1X2 - 1.30 X1X3 + 0.49 X1X4 + 1.84 X2X3 + 3.80 X2X4 +0.93 X3X4 - 2.76 X12 - 14.18 X22+6.53 X32 + 25.61 X42

Y2 (toluene removal by adsorption) = 58.49 + 2.62X1 + 3.54X2 + 1.79X3 - 9.00X4 + 2.20X1X2 - 1.86X1X3 + 0.49X1X4 + 0.86X2X3 + 2.89X2X4 + 0.023X3X4 + 9.53X12 - 14.64X22 + 5.25X32 + 11.83X42

Y3 (toluene removal by SOP) = 57.23+ 5.97X1 + 8.16X2 + 4.83X3 - 15.15X4 - 0.25X1X2 - 4.75X1X3 + 2.99X1X4 - 0.42X2X3 + 6.71X2X4 + 3.01X3X4 - 7.45X12-15.01X22 + 9.28X32 + 27.95X42

where, X1, X2, X3 and X4 are pH, MgO dose, ozonation time and initial concentration of toluene factors, respectively. The quality of the fitted model was expressed by R2. In the present study, RSM models were confirmed statistically for acceptability by using analysis of variance (ANOVA). The results are presented in Table 3.

Table 3.
Analysis of Variance for Response Surface Quadratic Model for Toluene Degradation by Catalytic Ozonation Process (COP/MgO)a

Figure 2 indicates a plot of predicted value versus experimental value of percentage removal from response surface design, where the coefficient of determination (R2) was 0.9480. As illustrated in Figure 2, the data points were well distributed near a straight line, which suggested a good correlation between the experimental and predicted values of the response and showed a good agreement with the model. Also, as presented in Figure 3 the quadratic model used in this work well gratified the ANOVA.

Figure 2.
Parity Plot for the Experimental and Predicted Values of Toluene Degradation (%) by Catalytic Ozonation Process

Figure 3.
Normal Probability Plot for the Residuals of Toluene Degradation by Heterogenic Catalytic Ozonation (COP) with MgO Nanoparticles

4.2. Interactive Effects of Different Parameters of COP/MgO on Toluene Degradation

In the current study, the three-dimensional response surface plots were used to determine the interaction among four variables for optimization of COP with MgO nanoparticles. The results are displayed in Figure 4A - F.

Figure 4.
Response Surface (3D) Plots for Optimization of Catalytic Ozonation Process with Catalyst Dose Nanoparticles

4.2.1. Effect of Solution pH

Figure 4A presents 3D plots showing the effects of solution pH (X1) and MgO dosage (X2) on the degradation of toluene using COP with MgO nanoparticles, while the ozonation time (X3) and toluene concentration (X4) values were fixed at their middle level (X3 = 27.5 min, X4 = 105 mg/L). As shown in Figure 4A, with the increase of pH and MgO dose, toluene degradation efficiency steadily increases. At initial pH of 12, with consumption of more than 0.3 g of MgO nanoparticles, > 70% degradation efficiency was achieved.

As it is clear from Figure 4B, with simultaneous increase of pH and time of ozonation, the degradation rate of toluene increased as well. Also, at pH values more than 12 and ozonation time more than 30 minutes, the increase of each of these two variables alone increased the removal efficiency, so that in alkaline conditions and ozonation time of 50 minutes, toluene degradation efficiency reached > 75%. Additionally, as shown in Figure 4C, maximum toluene degradation was achieved at high pH (= 12), and lowest toluene concentration (10 mg/L) under fixed conditions occurred at ozonation time of 27.5 minutes and MgO dose of 0.3 g. On the other hand, toluene degradation decreased with increase of toluene concentration.

4.2.2. Effect of MgO Dosage

As illustrated in Figure 4D, increased toluene degradation was observed with increased amount of MgO nanoparticles and ozonation time. At initial toluene concentration equal to 105 mg/L and pH 7.5, for reaction time of more than 50 minutes and MgO dose more than 0.3 g, toluene degradation efficiency reached about 75%.

Figure 4E indicates the simultaneous effect of MgO dose as a catalyst and concentration of toluene on the rate of toluene degradation. Under fixed condition for initial pH of the solution (7.5) and time of ozonation (27.5 minutes), increasing the toluene concentration had a negative effect on toluene degradation so that at concentrations greater than 105 mg/L of toluene, despite an increase in the amount of MgO nanoparticles, process efficiency dropped.

4.2.3. Effect of Ozonation Time

The 3D response surface plots according to two variables, ozonation time (minute) and toluene concentration (mg/L) is illustrated in Figure 4F, whereas the other two variables were kept at a central level (pH = 7.5 and MgO dose = 0.3 g). As indicated, the percentage of toluene removal increases with increase of ozonation time, so, the catalytic potential of MgO increased from 70.2% to 89.8% with an increase in the time of ozonation from 5 to 50 minutes, while, toluene concentration had a negative impact on toluene degradation.

5. Discussion

5.1. Model Fitting and Statistical Analysis

At present study, RSM models have been confirmed statistically for acceptability by using analysis of variance (ANOVA). Analysis of variance as a statistical method can be used as an exploratory tool to explain observations and testing hypotheses on the parameters of the model (40). The results are presented in Table 3. The F value was used for determination of statistical significance of the second-order regression model. As statted in Table 3, F value of 19.53 for the quadratic model indicating the significance of the model. In addition, a “Prob > F” lesser than 0.0001 for the second-order polynomial fitting shows that the model is statistically highly significant, and at 95% probability level, the model terms are significant. Values of “Prob > F” less than 0.05 show that the model terms are significant.

Actually, the larger F values and the smaller P values showed significance effects of the corresponding coefficients (41, 42). Findings of the present work revealed that the regressions for toluene degradation was statistically significant since P values were < 0.05 with larger F values. The table of ANOVA acquired from the response surface quadratic model for COP in degradation of toluene demonstrated that the linear model terms (pH, dose, time of ozonation and initial concentration of toluene) and interactive model term (Dose*Conc) were significant (P < 0.05) while other interactive terms were insignificant to the response.

Comparison of the experimental error of the replicate quantities can be performed by lack of fit (LOF). As stated in Table 3, the P value of LOF is more than 0.05 (= 0.1489), which confirm the significance of the used model. Additionally, “LOF F value” of 2.63 shows the LOF is non significant compared to the pure error. There is a 14.89% chance that a “LOF F value” this large can happen as a result of noise. Likewise, non-significant LOF is good. The LOF will be significant, if the model does not fit the data well.

The closer the R2 value to one, the better the model will be as it will give predicted values, which are near the actual values for the response. Furthermore, a rational agreement with adjusted R2 is needed (43). The value of predicted (model) multiple correlation coefficients (R2 = 0.7142) is in accordance with the value of adjusted R2 (= 0.8994). “Adeq Precision” measures the signal to noise ratio. A ratio more than four is appropriate. Signal to noise ratio for the present study was equal to 14.608, which represents adequate signal (43). Hence, optimization of the operational parameters can be performed by the quadratic model.

Figure 2 indicates a plot of predicted value versus experimental value of percentage removal from response surface design, where the coefficient of determination (R2) was 0.9480. As illustrated in Figure 2, the data points were well distributed near to a straight line, which suggested a good correlation between the experimental and predicted values of the response and showed good agreement with the model.

As presented in Figure 3 the quadratic model used in this work well gratified the ANOVA. The residual plots revealed a standard distribution, lending support to the assumption that pH, dose, time of ozonation, initial concentration of toluene and Dose*Conc. are significant terms, whereas, pH*Dose, pH*Time, pH*Conc., Dose*Time, and Time*Conc. are non significant terms. In fact, if the points on the normal probability plot are placed near the straight line, it suggests that the model fits the data well. As indicated by Fig. 3, almost all points lie near the straight line. Hence, there is an excellent correlation between developed regression model and experimental values and consequently the regression models can be applied to calculate the degradation of toluene in the initial experimental conditions.

5.2. Effect of Solution pH

The pH of the solution is a critical parameter in the catalytic ozonation process. In fact, the solution pH can affect the form in which organic compounds exist, the decomposition of the ozone molecules, and also the catalyst transformation (44). Also, it is well known that, ozone molecules can exclusively react with complexes of particular functional groups through selective reactions at acidic pH values. Under alkaline conditions, ozone molecules decompose yielding hydroxyl radicals, a highly oxidizing species that react with organic matters in a non-selective mode (15). On the other hand, it can be concluded that the alkaline conditions as well as the existance of a catalyst could stimulate hydroxyl radicals generation, which is satisfactory for organic matters decomposition (45, 46). Increase in the pH of the solution can accelerate the ozone mass transfer and also ozone decomposition rate through heterogeneous (catalytic) reaction. These reactions caused the formation of highly active radicals (HO• radical and others including OH, HO2 and HO3) and therefore can enhance the degradation rate of various pollutants (47, 48).

Based on Figure 4A, B and C it was found that, the best pH for ozone oxidation of toluene by COP is equal to 12. In addition, evaluating the potential of COP demonstrates the significant catalytic potential of MgO on ozonation process in toluene degradation, against SOP and adsorption processes, especially at higher pH (~12). As reported by Zhao et al. (49), an increase in the rate of nitrobenzene degradation in COP using manganese as a catalyst was observed with an increase in the solution pH from 3 to 11. Similar findings were reported by Moussavi and Mahmoudi (27) on degradation of a reactive red dye using COP/MgO and also Asgari et al. (50) on humic acid degradation by COP using MgO. Nevertheless, optimum pH for phenol removal by magnetic heterogeneous catalytic ozonation was found to be at neutral and weak alkaline (pH 7 - 8) (51). Freshour et al. (52) and Wu et al. (53) described that increase in available hydroxide ions and also acceleration of the ozone molecules decomposition at alkaline pH values can enhance the generation of hydroxyl and MgO-hydroxyl radicals (as oxidative radicals) and hence increase the degradation efficiency of pollutants.

The following mechanisms have been proposed to form hydroxyl and MgO-hydroxyl radicals (27):

O3 + (MgO - s) → (MgO - sO3 ) (3)

(MgO - sO3) → (MgO - sO•) + O2 (4)

(MgO - sO•) + 2H2O + O3 → (MgO - s˙OH) + 3˙OH + O2 (5)

(It must be considered that the suffix s in MgO-s indicates the Lewis acid site on the surface of MgO nanoparticles).

As stated in the above equation, radicals produced in the reaction, all lead to a progress in oxidation of toluene. In fact, the formed radicals in the solution under these conditions initiated toluene oxidation (as indirect) with greater rates than those obtained with SOP (15, 24). Similar mechanisms were reported by other researchers for organic compounds degradation (27, 54).

5.3. Effect of MgO Dosage

As shown, the effect of MgO dose was less than the ozonation time factor at minimum reaction times, so that at MgO dose of more than 0.3 g, and ozonation time of 11 minutes, 74.9% degradation efficiency was achieved. The increase of toluene degradation rate with an increase in catalyst dose can be attributed to the expansion of surface area of the catalyst and also the availability of more active sites for ozone molecules decomposition and subsequent reactions (18). In one previous study (21) on the catalytic ozonation of amoxicillin (with NH4Cl) and also phenol removal using magnetic heterogeneous catalytic ozonation by Dadban Shahamat et al. (51), similar results were reported. In addition, as presented by Rosal et al. (55), at fixed flow rate of ozone, greater increase of catalyst concentration can accelerate the ozone transfer from the gas media to the aqueous phase and also increase the decomposition rate of ozone molecules.

5.4. Effect of Ozonation Time

Reaction time, which is the time required to achieve the preferred objectives of purification, is one of the most critical parameters that affect the design and operation of an oxidation process. On the other hand, Figure 4F demonstrates that with an increased ozonation time and a decreased toluene concentration, degradation of toluene increases.

Similar conclusions were made (21) about amoxicillin antibiotic degradation and mineralization by the catalytic ozonation (with NH4Cl) and 4-chlorophenol removal using COP/MgO by Chen et al. (22). The results indicated that MgO nanoparticles extensively accelerated the rate of toluene degradation and are consequently a good catalyst to be applied for catalyzing the ozonation of model refractory contaminants. This catalytic potential of MgO in degradation of toluene could result in the presence of high density of basic functional groups that result in the decomposition of ozone molecules on the surface of MgO and also in the formation of active radicals (23).

5.5. Process Optimization

The optimum conditions for predicting maximum response (toluene degradation) have been achieved by using Minitab response optimizer as a global solution. Based on the software optimization step, the preferred objective for the operational conditions (pH of solution, MgO dosage, ozonation time and initial concentration of toluene) is chosen as “within” the range. The toluene degradation as response is defined as “maximum” to attain maximum efficiency. The program combined individual desirability to a single number and afterwards searched to improve this function according to the response goal. As illustrated in Figure 5, under optimized operational conditions, approximately 87.6% of toluene degradation rate was predicted by the model. The desirability function is found to be 0.6504 for the mentioned conditions as optimum. Furthermore, an extra experiment is then accomplished to approve these optimum findings, revealing agreement with the predicted response values.

Figure 5.
Optimum Conditions Derived by Response Surface Methodology for Toluene Degradation by Catalytic Ozonation Process /Catalyst Dose

5.6. Conclusion

In the present study, the efficiency of heterogenic catalytic ozonation with MgO nanoparticles on toluene degradation in synthetic solutions was investigated. A central composite design with response surface methodology was used to calculate the relationships between major critical factors, such as pH of solution, MgO dosage, time of ozonation and initial concentration of toluene to recognize the best conditions for the operation.

Based on the experimental data of this study, toluene degradation increased with increasing pH value, MgO dose and ozonation time. Nevertheless, performance of the process was decreased by increase of toluene concentration. The predicted result showed that maximum degradation of toluene (100%) could be achieved at pH 12, MgO dose 0.5 g, ozonation time 50 minutes and initial toluene concentration 10 mg/L. At this condition, experimental value was 99.99% for toluene degradation by COP. The correlation coefficient (R2 = 0.9480) indicates an excellent agreement between the predicted values and experimental ones. Additionally, the ANOVA of the quadratic model demonstrated that the used model was highly significant (P < 0.0001). At the optimum condition after process optimization by RSM, the efficiency for toluene degradation was found to be 87.6%. Finally, it was concluded that MgO nanoparticle is an efficient catalyst for toluene degradation and mineralization in the COP.


This article was derived from the PhD thesis of Ms. L. Mohammadi, PhD student of health promotion research center, Zahedan University of Medical Sciences. All authors wish to thank Dr. Hossein Kamani, Dr. Ferdos Kord Mostafapou and Ms. Shahnaz Sargazi for their support during analysis of the study.


Authors’ Contribution: The overall implementation of this study including design, experiments and data analysis, and manuscript preparation was the result of joint efforts by individuals, who are listed as coauthors of this paper.
Financial Disclosure: None declared.
Funding/Support: The authors are grateful to the Zahedan University of Medical Sciences for the financial support of this study (Project No. 7391).


  • 1. Zhao Z, Wang S, Yang Y, Li X, Li J, Li Z. Competitive adsorption and selectivity of benzene and water vapor on the microporous metal organic frameworks (HKUST-1). Chem Eng J. 2015;259:79-89. [DOI]
  • 2. Hinojosa-Reyes M, Arriaga S, Diaz-Torres LA, Rodríguez-González V. Gas-phase photocatalytic decomposition of ethylbenzene over perlite granules coated with indium doped TiO2. Chem Eng J. 2013;224:106-13. [DOI]
  • 3. Niri VH, Bragg L, Pawliszyn J. Fast analysis of volatile organic compounds and disinfection by-products in drinking water using solid-phase microextraction-gas chromatography/time-of-flight mass spectrometry. J Chromatogr A. 2008;1201(2):222-7. [DOI] [PubMed]
  • 4. Murata M, Tsujikawa M, Kawanishi S. Oxidative DNA damage by minor metabolites of toluene may lead to carcinogenesis and reproductive dysfunction. Biochem Biophys Res Commun. 1999;261(2):478-83. [DOI] [PubMed]
  • 5. Jo MS, Rene ER, Kim SH, Park HS. An analysis of synergistic and antagonistic behavior during BTEX removal in batch system using response surface methodology. J Hazard Mater. 2008;152(3):1276-84. [DOI] [PubMed]
  • 6. Barbee GC, Brown KW. Movement of xylene through unsaturated soils following simulated spills. Water, Air, and Soil Pollution. 1986;29(3):321-31. [DOI]
  • 7. Choi EJ, Jin HM, Lee SH, Math RK, Madsen EL, Jeon CO. Comparative genomic analysis and benzene, toluene, ethylbenzene, and o-, m-, and p-xylene (BTEX) degradation pathways of Pseudoxanthomonas spadix BD-a59. Appl Environ Microbiol. 2013;79(2):663-71. [DOI] [PubMed]
  • 8. WHO.. Guidelines for drinking-water quality World Health Organization, Distribution and Sales. Geneva, Switzerland: World Health Organization; 2004.
  • 9. Nadim F, Hoag GE, Liu S, Carley RJ, Zack P. Detection and remediation of soil and aquifer systems contaminated with petroleum products: an overview. J Petrol Sci Eng. 2000;26(1-4):169-78. [DOI]
  • 10. Seifi L, Torabian A, Kazemian H, Bidhendi G, Azimi AA, Nazmara S, et al. Adsorption of BTEX on Surfactant Modified Granulated Natural Zeolite Nanoparticles: Parameters Optimizing by Applying Taguchi Experimental Design Method. CLEAN - Soil, Air, Water. 2011;39(10):939-48. [DOI]
  • 11. Zhan S, Yang Y, Gao X, Yu H, Yang S, Zhu D, et al. Rapid degradation of toxic toluene using novel mesoporous SiO2 doped TiO2 nanofibers. Catalysis Today. 2014;225:10-7. [DOI]
  • 12. Solsona B, García T, Sanchis R, Soriano MD, Moreno M, Rodríguez-Castellón E, et al. Total oxidation of VOCs on mesoporous iron oxide catalysts: Soft chemistry route versus hard template method. Chem Eng J. 2016;290:273-81. [DOI]
  • 13. Deng Y, Zhao R. Advanced Oxidation Processes (AOPs) in Wastewater Treatment. Current Pollution Reports. 2015;1(3):167-76. [DOI]
  • 14. Kutschera K, Bornick H, Worch E. Photoinitiated oxidation of geosmin and 2-methylisoborneol by irradiation with 254 nm and 185 nm UV light. Water Res. 2009;43(8):2224-32. [DOI] [PubMed]
  • 15. Kasprzyk-Hordern B. Catalytic ozonation and methods of enhancing molecular ozone reactions in water treatment. Applied Catalysis B: Environmental. 2003;46(4):639-69. [DOI]
  • 16. Liotta LF, Gruttadauria M, Di Carlo G, Perrini G, Librando V. Heterogeneous catalytic degradation of phenolic substrates: catalysts activity. J Hazard Mater. 2009;162(2-3):588-606. [DOI] [PubMed]
  • 17. Qi F, Xu B, Chen Z, Feng L, Zhang L, Sun D. Catalytic ozonation of 2-isopropyl-3-methoxypyrazine in water by γ-AlOOH and γ-Al2O3: Comparison of removal efficiency and mechanism. Chem Eng J. 2013;219:527-36. [DOI]
  • 18. Zhao L, Ma J, Sun Z, Liu H. Influencing mechanism of temperature on the degradation of nitrobenzene in aqueous solution by ceramic honeycomb catalytic ozonation. J Hazard Mater. 2009;167(1-3):1119-25. [DOI] [PubMed]
  • 19. Moussavi G, Khavanin A, Alizadeh R. The investigation of catalytic ozonation and integrated catalytic ozonation/biological processes for the removal of phenol from saline wastewaters. J Hazard Mater. 2009;171(1-3):175-81. [DOI] [PubMed]
  • 20. Zeng YF, Liu ZL, Qin ZZ. Decolorization of molasses fermentation wastewater by SnO(2)-catalyzed ozonation. J Hazard Mater. 2009;162(2-3):682-7. [DOI] [PubMed]
  • 21. Shahidi D, Roy R, Azzouz A. Advances in catalytic oxidation of organic pollutants – Prospects for thorough mineralization by natural clay catalysts. Applied Catalysis B: Environmental. 2015;174-175:277-92. [DOI]
  • 22. Chen J, Tian S, Lu J, Xiong Y. Catalytic performance of MgO with different exposed crystal facets towards the ozonation of 4-chlorophenol. Applied Catalysis A: General. 2015;506:118-25. [DOI]
  • 23. Xing S, Lu X, Liu J, Zhu L, Ma Z, Wu Y. Catalytic ozonation of sulfosalicylic acid over manganese oxide supported on mesoporous ceria. Chemosphere. 2016;144:7-12. [DOI] [PubMed]
  • 24. Wu CH, Kuo CY, Chang CL. Decolorization of C.I. Reactive Red 2 by catalytic ozonation processes. J Hazard Mater. 2008;153(3):1052-8. [DOI] [PubMed]
  • 25. Wang J, Lou Y, Xu C, Song S, Liu W. Magnetic lanthanide oxide catalysts: An application and comparison in the heterogeneous catalytic ozonation of diethyl phthalate in aqueous solution. Separat PurifiTech. 2016;159:57-67. [DOI]
  • 26. Nawrocki J, Kasprzyk-Hordern B. The efficiency and mechanisms of catalytic ozonation. Applied Catalysis B: Environmental. 2010;99(1-2):27-42. [DOI]
  • 27. Moussavi G, Mahmoudi M. Degradation and biodegradability improvement of the reactive red 198 azo dye using catalytic ozonation with MgO nanocrystals. Chem Eng J. 2009;152(1):1-7. [DOI]
  • 28. Hammad Khan M, Jung JY. Ozonation catalyzed by homogeneous and heterogeneous catalysts for degradation of DEHP in aqueous phase. Chemosphere. 2008;72(4):690-6. [DOI] [PubMed]
  • 29. Dong Y, He K, Zhao B, Yin Y, Yin L, Zhang A. Catalytic ozonation of azo dye active brilliant red X-3B in water with natural mineral brucite. Catalysis Communications. 2007;8(11):1599-603. [DOI]
  • 30. He K, Dong YM, Li Z, Yin L, Zhang AM, Zheng YC. Catalytic ozonation of phenol in water with natural brucite and magnesia. J Hazard Mater. 2008;159(2-3):587-92. [DOI] [PubMed]
  • 31. Ma J, Sui M, Zhang T, Guan C. Effect of pH on MnOx/GAC catalyzed ozonation for degradation of nitrobenzene. Water Res. 2005;39(5):779-86. [DOI] [PubMed]
  • 32. Ikhlaq A, Brown DR, Kasprzyk-Hordern B. Catalytic ozonation for the removal of organic contaminants in water on alumina. Applied Catalysis B: Environmental. 2015;165:408-18. [DOI]
  • 33. Huang WJ, Fang GC, Wang CC. A nanometer-ZnO catalyst to enhance the ozonation of 2,4,6-trichlorophenol in water. Colloid Surf A. 2005;260:45-51.
  • 34. Gulková D, Šolcová O, Zdražil M. Preparation of MgO catalytic support in shaped mesoporous high surface area form. Microporous and Mesoporous Materials. 2004;76(1-3):137-49. [DOI]
  • 35. Moussavi G, khavanin A, Alizadeh R. The integration of ozonation catalyzed with MgO nanocrystals and the biodegradation for the removal of phenol from saline wastewater. Applied Catalysis B: Environmental. 2010;97(1-2):160-7. [DOI]
  • 36. Cadigan CA, Corpuz AR, Lin F, Caskey CM, Finch KBH, Wang X, et al. Nanoscale (111) faceted rock-salt metal oxides in catalysis. Catal Sci Technol. 2013;3(4):900-11. [DOI]
  • 37. Gonen F, Aksu Z. Use of response surface methodology (RSM) in the evaluation of growth and copper(II) bioaccumulation properties of Candida utilis in molasses medium. J Hazard Mater. 2008;154(1-3):731-8. [DOI] [PubMed]
  • 38. Srinivasan SV, Murthy DV. Statistical optimization for decolorization of textile dyes using Trametes versicolor. J Hazard Mater. 2009;165(1-3):909-14. [DOI] [PubMed]
  • 39. Anupam K, Dutta S, Bhattacharjee C, Datta S. Adsorptive removal of chromium (VI) from aqueous solution over powdered activated carbon: Optimisation through response surface methodology. Chem Eng J. 2011;173(1):135-43. [DOI]
  • 40. Segurola J, Allen NS, Edge M, Mahon AM. Design of eutectic photo initiator blends for UV/curable curable acrylated printing inks and coatings. Prog Org Coat. 1999;37:23-37.
  • 41. Demim S, Drouiche N, Aouabed A, Benayad T, Couderchet M, Semsari S. Study of heavy metal removal from heavy metal mixture using the CCD method. J Ind Eng Chem. 2014;20:512-20.
  • 42. Zheng Y, Wang A. . Removal of heavy metals using polyvinyl alcohol semi-IPN poly (acrylic acid)/tourmaline composite optimized with response surface methodology. Chem Eng J. 2010;162:186-93.
  • 43. Nordin MY, Venkatesh VC, Sharif S, Elting S, Abdullah A. Application of response surface methodology in describing the performance of coated carbide tools when turning AISI 104 steel. J Mater Process Technol. 2004;145:46-58. [DOI]
  • 44. Xiao H, Xu Y, Yu M, Zhang Q. Enhanced mineralization of 2,4-dichlorophenol by ozone in the presence of trace permanganate: effect of pH. Environ Technol. 2010;31(11):1295-300. [DOI] [PubMed]
  • 45. Yan H, Lu P, Pan Z, Wang X, Zhang Q, Li L. Ce/SBA-15 as a heterogeneous ozonation catalyst for efficient mineralization of dimethyl phthalate. J Molecular Catalysis A: Chem. 2013;377:57-64. [DOI]
  • 46. Yan H, Chen W, Liao G, Li X, Ma S, Li L. Activity assessment of direct synthesized Fe-SBA-15 for catalytic ozonation of oxalic acid. Separat Purifi Techno. 2016;159:1-6. [DOI]
  • 47. Valdes H, Murillo FA, Manoli JA, Zaror CA. Heterogeneous catalytic ozonation of benzothiazole aqueous solution promoted by volcanic sand. J Hazard Mater. 2008;153(3):1036-42. [DOI] [PubMed]
  • 48. von Gunten U. Ozonation of drinking water: Part I. Oxidation kinetics and product formation. Water Research. 2003;37(7):1443-67. [DOI]
  • 49. Zhao L, Ma J, Sun Z, Zhai X. Catalytic ozonation for the degradation of nitrobenzene in aqueous solution by ceramic honeycomb-supported manganese. Applied Catalysis B: Environ. 2008;83(3-4):256-64. [DOI]
  • 50. Asgari G, Faradmal J, Zolghadr Nasab H, Seidmohammadi A. Catalytic Potential of Nano-Magnesium Oxide on Degradation of Humic Acids From Aquatic Solutions. Avicenna J Environ Health Eng. 2014;1(1) [DOI]
  • 51. Shahamat YD, Farzadkia M, Nasseri S, Mahvi AH, Gholami M, Esrafili A. Magnetic heterogeneous catalytic ozonation: a new removal method for phenol in industrial wastewater. J Environ Health Sci Eng. 2014;12(1):1.
  • 52. Freshour AR, Mawhinney S, Bhattacharyya D. Two-phase ozonation of hazardous organics in single and multicomponent systems. Water Res. 1996;30(9):1949-58. [DOI]
  • 53. Wu CH, Kuo CY, Chang CL. Homogeneous catalytic ozonation of C.I. Reactive Red 2 by metallic ions in a bubble column reactor. J Hazard Mater. 2008;154(1-3):748-55. [DOI] [PubMed]
  • 54. Faria PCC, Órfão JJM, Pereira MFR. Activated carbon catalytic ozonation of oxamic and oxalic acids. Applied Catalysis B: Environmental. 2008;79(3):237-43. [DOI]
  • 55. Rosal R, Gonzalo MS, Rodríguez A, Perdigón-Melón J, García-Calvo E. Catalytic ozonation of atrazine and linuron on MnOx/Al2O3 and MnOx/SBA-15 in a fixed bed reactor. Chem Eng J. 2010;165(3):806-12. [DOI]

Table 1.

Experimental Range and Level of the Independent Variables

Independent Variables Coded Symbol Coded Levels And Actual Factor Levels
-1 0 +1
pH X1 3 7.5 12
MgO dose, mg/L X2 0.1 0.3 0.5
Ozonation time, min X3 5 27.5 50
Toluene concentration, mg/L X4 10 105 200

Table 2.

Experimental Design Matrix for Toluene Degradation

Run X1 X2 X3 X4 SOP Efficiency, % Ads. Efficiency, % COP Efficiency, %
Exp. Pred. Exp. Pred. Exp. Pred.
1 3 0.5 5 10 78.90 80.22 70.42 75.68 89.90 88.09
2 3 0.1 50 200 50.70 53.72 59.00 49.18 62.00 61.59
3 7.5 0.16 27.5 105 52.60 51.21 52.51 46.15 63.00 57.77
4 12 0.5 50 200 80.00 81.51 73.40 76.97 90.00 89.38
5 12 0.5 5 200 80.98 74.46 75.98 69.92 87.00 82.32
6 7.5 0.3 27.5 105 64.50 65.54 57.20 58.32 72.00 68.29
7 12 0.5 5 10 94.20 91.96 82.90 87.42 99.99 99.83
8 7.5 0.3 27.5 171.5 67.90 65.81 60.30 60.75 75.00 72.37
9 3 0.5 50 200 72.90 73.02 65.30 68.48 80.00 80.88
10 7.5 0.3 43.25 105 72.00 66.38 63.80 61.32 79.00 72.94
11 7.5 0.3 27.5 105 51.90 65.54 52.40 58.32 63.00 68.29
12 3 0.1 5 200 10.00 48.82 50.00 44.28 59.00 56.69
13 12 0.1 5 10 89.00 90.81 80.00 86.27 98.70 98.67
14 3 0.1 50 10 81.50 84.67 76.57 80.13 87.00 92.53
15 7.5 0.3 27.5 105 49.80 65.54 59.00 58.32 65.00 68.29
16 3 0.1 5 10 83.00 83.48 78.00 78.94 91.30 91.35
17 12 0.1 50 200 50.80 57.80 54.40 53.27 63.00 65.67
18 7.5 0.3 27.5 105 66.10 65.54 55.80 58.32 68.50 68.29
19 7.5 0.44 27.5 105 49.90 58.37 51.80 53.31 62.00 64.93
20 10.65 0.3 27.5 105 51.90 63.16 62.30 58.10 67.60 69.72
21 7.5 0.3 27.5 105 53.50 65.54 61.90 58.32 67.30 68.29
22 7.5 0.3 27.5 105 48.00 65.54 58.70 58.32 65.90 68.29
23 7.5 0.3 27.5 38.5 76.70 82.77 69.95 77.71 89.00 89.33
24 12 0.1 5 200 51.50 58.12 55.00 53.58 63.70 65.99
25 3 0.5 50 10 92.60 88.77 83.40 84.24 99.50 96.64
26 3 0.5 5 200 53.00 60.74 57.10 56.20 64.82 68.61
27 7.5 0.3 11.75 105 54.30 63.49 60.00 58.44 66.30 70.06
28 12 0.5 50 10 93.00 95.29 85.00 90.76 99.99 100.0
29 4.35 0.3 27.5 105 58.00 57.62 65.70 52.56 68.60 64.18
30 12 0.1 50 10 88.70 86.77 80.00 82.23 99.00 94.64

Table 3.

Analysis of Variance for Response Surface Quadratic Model for Toluene Degradation by Catalytic Ozonation Process (COP/MgO)a

Source SS DF MS RMS Prob > F, (P Value)
Model 5534.36 14 395.31 19.53 < 0.0001 Highly significant
X 1 265.63 1 265.63 13.12 0.0025 Significant
X 2 443.71 1 443.71 21.92 0.0003 Significant
X 3 72.02 1 72.02 3.56 0.0788 Highly significant
X 4 2490.93 1 2490.93 123.03 < 0.0001 Significant
X 1 *X 2 19.49 1 19.49 0.96 0.3420
X 1 *X 3 27.25 1 27.25 1.35 0.2641
X 1 *X 4 3.90 1 3.90 0.19 0.6670
X 2 *X 3 54.32 1 54.32 2.68 0.1222
X 2 *X 4 230.58 1 230.58 11.39 0.0042 Significant
X 3 *X 4 13.84 1 13.84 0.68 0.4213
Residual 303.69 15 20.25 - - -
Lack of fit 255.15 10 25.52 2.63 0.1489
Pure error 48.53 5 9.71 - - -
Cor Total 5838.04 29 - - - -
Abbreviations: SS, sum of squares; DF, degree of freedom; MS, mean square; RMS, ratio of mean square (F value).
a Adeq Precision = 14.608, R2 Adjusted = 0.8994, R2 = 0.9480, Predicted R2 = 0.7144.

Figure 1.

Schematic of the Used Catalytic Ozonation Reactor

Figure 2.

Parity Plot for the Experimental and Predicted Values of Toluene Degradation (%) by Catalytic Ozonation Process

Figure 3.

Normal Probability Plot for the Residuals of Toluene Degradation by Heterogenic Catalytic Ozonation (COP) with MgO Nanoparticles

Figure 4.

Response Surface (3D) Plots for Optimization of Catalytic Ozonation Process with Catalyst Dose Nanoparticles

Figure 5.

Optimum Conditions Derived by Response Surface Methodology for Toluene Degradation by Catalytic Ozonation Process /Catalyst Dose

Equation 1.

Equation 2.