Kinetics and statistical analysis of the bio-stimulating effects of goat litter in crude oil biodegradation process

The kinetics and statistical analysis of crude oil (CO) degradation in CO-contaminated soil (COCS) using goat litter (GL) were investigated. The data obtained from the CO degradation process was fitted to the first- and second-order kinetic models. The effects of process parameters such as temperature and the initial CO concentrations on the CO degradation process were also investigated. The one-way ANOVA and Turkey’s post-hoc analysis were also used to study the statistical significance of the process parameters on the CO degradation process. The microbial count showed that the GL contained a total viable count (TVC), coliform, and mold counts of 2.6 × 107 CFUg−1, 2.6 × 107 CFUg−1, and 6.9 × 103 CFUg-1, respectively. The error and linear regression analysis between experimental and model-predicted values revealed that the first-order kinetic model gave a better explanation of the CO degradation process. The rejection of the null hypothesis was evident from one-way ANOVA and Turkey’s post-hoc analysis as the P values at a temperature of 30 °C and initial CO concentrations of 70 gL−1 and 90 gL−1 were less than the significant level of 0.05. Notable organic nutrients in the GL which were beneficial in the COCS treatment process as indicated by the Fourier-transform infrared spectrophotometer (FTIR) analysis were phosphorous and nitrogen. It may be concluded that GL could be used as an effective organic treatment for COCS at CO initial concentrations of 70 and 90 gL-1 and a temperature of 30 °C.


Background
The activities of humans in the exploration, refining, transportation, storage, and usage of crude oil (CO) are some of the many ways by which CO is steadily being introduced into the environment. CO in its refined state could be used in the generation of energy and production of valuable industrial raw materials. In recent years, the spillage of CO and its refined products into the environment has been a source of environmental concern as the available arable soil is steadily being contaminated [1,2].
There have been lots of researches and innovations in the area of environmental remediation due to the increasing demands from public and government regulations on the need to scale down CO spill within the environments. A number of conventional technologies, which have been used for the remediation of hydrocarbon contaminated environments, include soil flushing, solvent extraction, electro-kinetic remediation, photo-catalytic degradation, etc. [3,4]. These conventional technologies are associated with high capital intensive, ecological unacceptable, and might not viable for developing nations like Nigeria.
However, microbial degradation of hydrocarbon has begun to gain wider acceptance as an efficient alternative for hydrocarbon degradation due to its simplicity and environmentally acceptable [5,6]. The different catabolic pathways developed by microorganisms enable them to efficiently utilize the insoluble hydrocarbon. One of the pathways developed by microorganism includes the high affinity to hydrophobic surfaces, which enable them to directly absorb and utilize the hydrocarbon contaminant [7,8]. Also, during microbial metabolisms, enzymes such as dioxygenases, catalases, and peroxidases are released [9]. These enzymes could oxidize and subsequently transform the hydrocarbon contaminant into less toxic compounds [10].
Bio-stimulation, an aspect of bioremediation entails the addition of organic materials rich in the microbial population, micro-and macronutrients in the contaminated soil [11]. The advantage of bio-stimulation is that the bioremediation process will be carried out by both the microbes in the added organic materials and the contaminated soil [2]. Furthermore, the added organic materials during bio-stimulation could improve the fertility of the contaminated soil [1,12,13].
In this present study, goat litter (GL) was used as biostimulating organic material to study the kinetics and statistical analysis of CO degradation. The application of GL as organic material has been known to improve the physicochemical properties of the soil, beneficial to soil microorganism, and enhance the overall fertility of the soil. Awodun et al. [14] reported that GL was better than other livestock manure as its moisture content is much lower, thus leading to a slow dissipation of its nutrients into the soil. Similarly, the neutral pH of the GL was also found to be conducive for microbial growth and proliferation [14]. However, GL is a major source of nitrogen, phosphorous, potassium, calcium, magnesium, and organic matter, which are necessary for crop production [15].
Accordingly, the objectives of this study were to investigate the kinetics and statistical analysis of CO degradation in the crude oil-contaminated soil (COCS). The effects of process parameters such as the initial CO concentrations and temperature were also investigated. The prevalent functional groups in the GL and COCS were also investigated using Fourier transform infrared (FTIR). Furthermore, one-way analysis of variance (ANOVA) and Turkey's post-hoc analysis were used to evaluate the statistically significant (P < 0.05) effects of the process parameters on the CO degradation process.

Samples collection
The GL used in this study was collected from a farmhouse located at Ogwofia-Ozom Mgbagbu Owa in Enugu state Nigeria. The GL was air-dried for 2 weeks, ground and passed through 5-mm sieve for homogeneity before use. The CO was collected from Warri petrochemical refinery in Delta State Nigeria. An agricultural soil with no history of hydrocarbon contamination was collected from the soil science garden of the University of Nigeria Nsukka [16]. The samples were taken to the soil science laboratory of the University of Nigeria Nsukka for further analysis.

Preparation of CO stock solution
The CO stock solutions used in this experiment were prepared by weighing out 50, 70, and 90 g of CO using PCE analytical weighing balance (PCE-6000). Each of the CO mass samples was dissolved in 1.0 L of distilled water to give initial CO concentrations of 50 gL -1 , 70 gL -1 , and 90 gL −1 [16]. The soil was artificially contaminated by spiking the various prepared CO concentrations on 200 g of the clean soil sample. The crude oil-contaminated soil (COCS) samples were allowed to stay for 21 days to permit the volatilization and sorption of the CO into the soil matrix. After 21 days, the COCS was supplemented with 100 g of GL as bio-stimulating organic material and was designated as (COCS-GL).

CO degradation procedure
The prepared 200 g of COCS was placed in a plastic container with dimension 6.5 × 6.5 × 3.5 cm. Control samples were also prepared with the same quantity of soil and CO concentrations but with no GL [16]. Both the control samples and treated samples were prepared in triplicates, aerated, and mixed properly twice a week. This was to enhance the supply of oxygen, which was required to maintain microbial respiration during the CO degradation process [13,17]. The contaminated soils in the treatment containers were covered with polyethylene foil after mixing to prevent the influence of atmospheric weather conditions [17]. The performance of the process was monitored and samples were collected on the 7, 14,21,28,35,42,49, and 56 days of treatment for analysis and also to calculate the percentage CO degraded in these periods.

Determination of CO percentage degradation
The solvent extraction procedure of gas chromatography equipped with a splitless injector and flame ionization detector was used to estimate the residual CO concentrations [2,12]. The dried sub-samples were extracted for 2 h using acetone/dichloromethane (1:1 v/v) as the extracting solvent. The solvent was left to evaporate after the extraction process, and the residue (extract) was dissolved in 5 mL dichloromethane. The temperature range was maintained at 40°C for 2 min before increasing it to 320°C at a rate of 7°C/min. The rate was increased to 20°C/min until the temperature reached 400°C, after which, the temperature was constant for 10 min. The recovery of CO using this method was higher than 68%. The CO concentrations were determined after proper calibration of the method with standard CO samples at different concentrations. The percentage of CO degraded was calculated using the expression in Eq. 1: Where; the initial CO concentration (gL −1 ) is at time = 0, and the final CO concentration (gL −1 ) is at time = t

Microbial count
A total of 1 g of the samples (GL and COCS) was placed in a test tube containing 9 mL Ringer solution. The Ringer solution had the following compositions (g/mL); 6.5 NaCl, 0.25 CaCl 2 , and 0.2 NaHCO 3 . The mixtures were mixed thoroughly at 150 rpm on a rotary shaker. The mixture was further diluted by transferring 1 mL into three test tubes containing a 9 mL Ringer solution using a sterile pipette. Thereafter, 0.1 mL of the diluted solution was transferred from each dilution test tube into three sterile Petri dishes using a sterile pipette. For total viable count, (TVC) mold count, and coliform count, 15 mL of sterile nutrient agar, Sabouraud Dextrose Agar (SDA) and sterile violet red bile agar, respectively were poured into each of the Petri dishes and gently rotated on the table to mix. The nutrient agar, SDA, and sterile violet red bile agar in the Petri dishes were allowed to gel for 10 min before sealing and were incubated in an inverted position at 37°C for 24 h. The colonies formed were counted using a Quebec colony counter (Reichert dark field 13332500/13332600) and was expressed as the colony-forming unit per gram (CFUg¯1).

Effect of temperature
An incubation experiment was set up in a water bath (Stuart SBS40 Shaking water bath) with a temperature regulator and adjustable shaking speed. The water bath platform accommodated four Erlenmeyer flasks of 250 mL. Each of the four flasks contained 100 g of the COCS-GL. The samples were incubated at 10°C, 20°C, and 30°C. The adjustable shaking speed of the water bath was set at 150 rpm. A polycarbonate water bath cover was used to cover the water bath in order to prevent heat loss by evaporation in either of the stainless steel interior or from the incubated samples. The water bath was capable of measuring temperatures within the ranges of 10°C to 90°C (±0.25). The set temperature of the water bath was monitored through the easy to read light-emitting diode (LED) display screen for a contact time of 30 min. The experiments were conducted in triplicates and the statistical mean of the results was recorded.

Statistical analysis
The significant effects of the initial CO concentration and temperature on CO degradation were evaluated using one-way ANOVA and Turkey's post-hoc multiple comparisons. This was in order to estimate the variance within and between the group means. The one-way ANOVA and Turkey's post-hoc tests were processed using IBM SPSS (version 22). However, the Turkey's post-hoc was used to predict the honest significant difference that must have existed in distance between two groups mean. The statistically significant difference between the two group's means was ascertained at P < 0.05.

Error analysis of the first-and second-order kinetic models
The error functions in Table 1 were employed to evaluate the suitable kinetic model that best describes the experimental data on CO degradation. The minimization of the fractional error across the entire CO concentration range was achieved using the average relative error deviation (ARED) [18,19]. The root-mean-square error (RMSE) was used to evaluate the residuals errors between the experimental and model-predicted values [20]. The standard error of prediction (SEP) measures the accuracy of the model prediction in comparison with the experimental values ( Table 1). The magnitudes of errors between the predicted and experimental values were estimated using the mean absolute error (MAE). The correlation coefficient (R 2 ) was used to measure how close the data are to the fitted regression line [21]. In summary, the goodness of fit of either model (first-and second-order) was evaluated with higher R 2 values and lower values of ARED, RMSE, MAE, and SEP [22].

Microbial count
The total viable count (TVC), coliform, and mold counts were used to ascertain the available microbial population in the GL and COCS (Table 2). Balba et al. [23] emphasized that the soil analysis of the total heterotrophic microbial count would indicate the degree at which the soil could support contaminant degradation.

Factor sensitivity analysis 3.3.1 Effect of initial CO concentration on CO degradation
CO concentration influenced the microbial degradation of CO as observed in both COCS-GL and the control. The effect of the initial CO concentration on CO degradation in the contaminated soil was investigated at different CO concentrations (50, 70, and 90 gL −1 ). The Values represent ± standard deviation (SD) of triplicates measurement, TVC is the total viable count in CFUg −1 Fig. 3 a-c Effect of initial CO concentration on CO degradation (bars represents ± SD of triplicates measurements) Ani and Chukwuma Beni-Suef University Journal of Basic and Applied Sciences (2020) 9:29 Page 5 of 13 temperature and quantity of GL used were fixed at 30°C and 100 g, respectively. Figure 3a-c showed the levels of CO reductions in the CO concentrations during the remediation period. After 56 days of treatment, the CO concentrations of 50, 70, and 90 gL −1 for COCS-GL reduced to 26.1, 21.2, and 15.6 gL −1 , respectively. These reductions resulted in CO degradation rates of 49.4%, 69.6%, and 89.8%, respectively. However, in the control, the initial CO concentrations reduced to 44.4, 59.1, and 70.8 gL −1 for 50, 70, and 90 gL −1 , which corresponded to 11.2%, 15.5%, and 21.3% degradation rates, respectively (Fig. 3a-c).

Effects of temperature on CO degradation
It is well established that temperature is a critical environmental factor affecting microbial degradation [24]. Batch experiments were performed with COCS-GL and the control at three different temperatures (10°C, 20°C, and 30°C) to determine the effect of this factor on the microbial degradation of CO. The water bath shaking speed was maintained at 150 rpm for 30 min contact time and constant CO concentration. The TVC of soil microbial population at the incubation temperatures were observed to determine the optimum temperature with higher TVCs. In Fig. 4, the TVCs of 15.8 × 10 6 CFUg −1 and 6 × 10 6 CFUg −1 were observed for COCS-GL and control, respectively at 30°C. The CO degradation rates of (52.1%) and 25% were recorded at 30°C for COCS-GL and control, respectively (Fig. 5).

Statistical analysis 3.4.1 One-way ANOVA of temperature effect on CO degradation
The result from the one-way ANOVA for temperature effect on CO degradation was presented in Table 3. The one-way ANOVA was used to judge the statistically significant difference between and within the group means. The one-way ANOVA results in Table 3 indicate a statistical significance difference between and within the group means (P < 0.05), thus excluding the null hypothesis. The result further indicated that the mean range of temperatures (10°C, 20°C, and 30°C) on CO degradation was statistically significant at P < 0.05.
The statistical significance (P < 0.05) of all pairs in the group was further examined using the Turkey's post-hoc analysis as shown in Table 4. From the Turkey's posthoc analysis, asterisks (*) were used to denote the pairs that were statistically significant at P < 0.05. The Turkey's post-hoc analysis in Table 4 revealed that the mean difference (I-J) within the temperature (I) of 30°C were all significant (P < 0.05). Consequently, the pairs within 10°C and 20°C were not all significant (P > 0.05). The statistical significance of the mean temperature   difference at 30°C suggests that the microbial degradation of CO was enhanced at 30°C. Table 5 shows the one-way ANOVA result for the effect of initial CO concentration on CO degradation. Accordingly, the rejection of the null hypothesis was as a result of the statistical significance difference (P < 0.05) that exists between and within the group means. The differences between the CO concentration means were due to the variation in the initial CO concentrations. Further studies of the statistical group means were investigated using the Turkey's post-hoc analysis ( Table 6). These investigations aimed to identify the specific group pairs that were indeed statistically significant (P < 0.05). From the Turkey's post-hoc analysis in Table 6, asterisks (*) were used to identify statistically significant pairs with P < 0.05. The Turkey's post-hoc results indicated that the mean difference (I-J) within the CO concentrations (I) of 70 and 90 gL −1 were all statistically significant (P < 0.05). However, the Turkey's post-hoc multiple comparisons within the 50 gL −1 initial CO concentration were not all statistically significant (P > 0.05). The above result revealed a significant CO degradation at 70 and 90 gL −1 initial CO concentrations.

Process kinetics 3.5.1 First-and second-order CO degradation kinetics
In order to evaluate the CO degradation kinetic parameters for COCS-GL and the control, experimental data were fitted to the nonlinear and linear forms of first-and second-order kinetics models (Table 7). From Table 7, the initial and final CO concentration in the COCS at time t is represented by [Co] and [Ct], respectively. K 1 and K 2 are the first-and second-order CO degradation rate constants (day −1 ), respectively while t is the time in days [26]. The first order plot of ln[Ct] versus t yields a straight line with K 1 and ln[Co] as the slope and intercept, respectively [27]. On the other hand, the second-order plot of 1 ½Ct versus t yields a straight line with K 2 and 1 ½Co as the slope and intercept, respectively ( Table 7).
The nonlinear CO degradation rate constants corresponding to each initial CO concentration were obtained using the linear parameters as the initial guess [25]. The variation in experimental and theoretically obtained data as a function of initial CO concentration and time for the nonlinear regression is presented in Figs. 6a-c, 7a-c, 8a-c, and 9a-c. The kinetic parameters at different initial CO concentrations for first-and second-order kinetic models are presented in Tables 8 and 9.
From Tables 8 and 9, the first-and second-order rate constants obtained for CO degradation were higher in the COCS-GL. Within the range of CO concentrations (50 to 90 gL −1 ) used in this study, the CO degradation rate constants for the first-and second-order kinetic models (K 1 and K 2 ) were observed to be highest at 90 gL −1 while lower rate constants were obtained at 50 gL −1 (Tables 8  and 9). The biological half-lives, which were the time taken by the indigenous microbial population to degrade half of the initial CO concentrations, were also calculated using Eqs. (2) and (3). The first order biological half-life is dependent on the first-order CO degradation rate constant (K 1 ) according to Eq. (2). In contrast, the secondorder biological half-life is dependent on the second-order CO degradation rate constant (K 2 ) and the initial CO concentration as expressed in Eq. (3): The results obtained from the biological half-lives indicated that the microbial population took more time to degrade the CO at 50 gL −1 initial CO concentration. For example, in COCS-GL and for first-order kinetic model, the microbial population took 22 days to degrade half of the initial CO concentration of 50 gL −1 ; whereas, it took 17 and 11 days, degrade half of the 70 and 90 gL −1 initial CO concentrations, respectively.

Error analysis of the kinetic models
The predictive capabilities of both first-and secondorder kinetic models for CO degradation using the error functions in Table 1 were compared. The statistical error analysis employing the error functions such as RMSE, ARED, MAE, SEP, and R 2 for both the first-and second-order kinetic models have been calculated and presented in Table 10.
The results showed that the data obtained for the second-order kinetic model were higher than the firstorder kinetic model (Table 10). These observations revealed that lower values of the RMSE, ARED, MAE, and SEP were obtained for the first-order kinetic model. The results further indicated that the first-order kinetic model gave a better fit to the experimental data on CO degradation. Furthermore, linear regression analysis between the CO degradation data predicted by the first-and secondorder kinetic models with their corresponding experimental values was investigated ( Fig. 10a and b). The R 2 values were used to judge how close the data are to the fitted regression line. In Fig. 10a and b, the first-order model predictions lie much closer to the fitted regression line in comparison with the second-order kinetic model. The R 2 values for the first-order kinetic model were 0.977 and 0.971 for COCS-GL and control, respectively. On the other hand, the R 2 values for the second-order kinetic model were 0.923 and 0.887 for COCS-GL and control, respectively. These results revealed that the first-order kinetic model performed better. However, the R 2 values for COCS-GL were higher than the control for both the first-and second-order kinetic models.

Characterizations of the COCS and GL
The prevalent functional group of the COCS was investigated using the FTIR technique as presented in Linear second-order  [28] attributed the peak located at 3471.98 cm −1 to the first overtone vibration of -CH 2 symmetric deformation in total petroleum hydrocarbon (TPH). This band was not observed in the spectrum pattern of GL and could be attributed to the amount of spiked CO in the COCS [28]. Also, in the spectrum pattern of TPH, Stuart [29] attributed the peak at 1605.79 cm −1 to the aliphatic hydrocarbon compound of an alkene. However, in hydrocarboncontaminated soil, Mohd et al. [30] reported that hydroxyl and alkyne group with stretching vibrations of O-H and C≡C were observed at 1033.88 cm −1 and 728.16 cm −1 . The series of weak absorption peaks in the spectrum of COCS (Fig. 1) were related to TPH and were attributed to the vibrational frequencies of terminal methyl (-CH 3 ) and alkyl halides [28]. However, the FTIR peaks of the COCS indicated the presence of aliphatic and aromatic hydrocarbon compounds (Fig. 1).
The spectrum pattern of GL as indicated in Fig. 2 shows the functional groups of known organic nutrients as reported previously. These nutrients could be  contributing to the CO degradation process. Pavia et al. [31] reported that N-H stretching of amine could also be assigned at the absorption bands at 3770.96 and 3421.83 cm −1 . Also, the peaks at 2939.61 and 2396.63 cm −1 were indicative of the presence of phosphorous and ester P-H stretching [29]. However, Stuart [29] and Pavia et al. [31] noted that the bands at 2939.61 and 2396.63 cm −1 represent N-CH 2 stretching vibration of tertiary amine (Fig. 2). Likewise, Stuart [29] attributed the band at 1641.48 cm −1 to the vibration of the N-H primary amine group, which was indicative of nitrogen may be due to the carbonyl stretching vibration. Stevenson [32] and Stuart [29] assigned the absorption band at 1424.48 cm −1 , 1241.23 cm −1 to the C-O axial deformation of polysaccharides. From the observed FTIR peaks, it was indicative that phosphorous and nitrogen compounds were present in GL (Fig. 2).

Microbial count
Observations made from the microbial count results showed that the addition of GL in COCS could affect the general nutrient cycling and microbial biomass of the contaminated soil. This was as the microbial count (TVC, coliform, and mold) in GL was high (Table 2). Previously, Czurak-Dainard [33] and Alkatib et al. [34] reported that soils with histories of pig and goat manure application generally contained higher concentrations of microbial biomass carbon compared with soil with no manure. The low microbial counts in the control could be attributed to the microbial competition for the scarce nutrient in the COCS [35,36]. However, the low microbial population and organic matter content of CO-contaminated soils have been reported previously [37].

Effect of initial CO concentration
The CO degradation rate was observed to have increased at 90 gL −1 in both the COCS-GL and control. Accordingly, 89.8% and 21.3% CO degradation rates were obtained at 90 gL −1 for COCS-GL and control, respectively. The bioavailability of the CO at 90 gL −1 might be sustaining the required microbial action on CO. The high rate of CO degradation in the COCS-  GL indicated that the COCS treated with GL enhanced CO degradation. Similar results on the use of organic nutrients to enhance the bioremediation of contaminated environments have been previously reported [12,13,17]. Furthermore, the FTIR result for GL revealed that nitrogen and phosphorous were the dominant organic nutrients in GL [2]. These nutrients may have contributed to the CO degradation due to the increase in the microbial count (TVC, coliform, and mold). Also, previous studies reported that these nutrients were responsible for enhanced contaminant degradation by increasing the growth rate of the available microbial population [38].

Effects of temperature
The majority of the soil microorganisms were mesophilic with maximum growth temperature between 25°C and 45°C [24,39,40]. This was in accordance with the TVC (15.8 × 10 5 CFUg −1 ) observed in COCS-GL at 30°C, where it was obvious that the added organic material (GL) favored the TVC in COCS (Fig. 4). However, in the relationship between soil microbial population and count, Critter et al. [41] demonstrated that different organic amendments significantly affected the microbial quantity and count. The CO degradation examined at the respective temperatures (10°C, 20°C, and 30°C) showed that the maximum CO degradation (52.1%) was obtained at 30°C (Fig. 5). This could be attributed to the maximum TVC recorded at 30°C. Similarly, Muftah et al. [42] reported that at a temperature of about 30°C, the activity of P. putida and its ability to degrade phenol was optimized.

Process kinetics for CO degradation
The biological half-live for the first-order kinetic model indicated that 90 gL −1 (11 days) degrades faster than 50 gL −1 (22 days) (Tables 8 and 9). This observation revealed that the 90 gL −1 of CO could be satisfying the microbial demands for carbon while 50 gL −1 might not be enough for microbial action on CO. However, these observations were in conformance with the works of Zappi et al. [43] where the degradation of poly-aromatic hydrocarbon did not occur at low concentration even an additional carbon source was added. The CO degradation rate constant for the first-order kinetic model (K 1 ) showed a faster CO degradation at 90 gL −1 (0.069 day −1 ) in comparison with the rate constant obtained at 50 gL −1 (0.045 day −1 ) for both the COCS-GL and control. A similar observation was made in the second-order kinetics model (Tables 8 and 9). The values of the rate constants and half-lives obtained in this study were lower than those obtained previously in the bioremediation of total petroleum hydrocarbon (TPH) [11,44]. This could be due to different experimental processes, analysis, and organic material used in the remediation process. Coupled with the addition of GL in the COCS, the aeration and mixing process of the COCS treatment might also be enhancing CO degradation [45]. The aeration and mixing processes gave rise to the increased microbial activity and aerobic nature of CO degradation [45]. Also, Nilanjana and Preethy [46] emphasized that the most rapid and complete degradation of the majority of organic pollutants is brought about under aerobic conditions. The CO degradation in the control could be attributed to the stimulation of the indigenous microbial population through the aeration and mixing processes as GL was not added. Similarly, Mohsen et al. [45] reported that the degradation of hydrocarbon in the control samples could be explained by the activity of the microbial population due to the aeration, and mixing processes.

Error analysis of the kinetic models
The lower values of the error analysis for the first-order kinetic model showed that the first-order kinetic model  performed better in the modeling of the CO degradation. Previous works reported that the lower the values of the error functions, the better the model goodness of fit (Kitanovic et al. 2008). Also, in the study of hydrocarbon modeling and kinetics, Anders and Haller [47] reported that the simplicity of the first-order kinetic model permits its application in degradation studies as the degradation rate is directly dependent on substrate concentration. Similarly, Jorgensen [48] emphasized that the degradation of the hydrocarbon compound is governed by the first-order kinetics model where the degradation rate of the hydrocarbon compound is proportional to its concentration.

Conclusion
The bio-stimulation of crude oil-contaminated soil using GL has been achieved. The microbial characterizations showed that GL contained active microbial populations, which were found useful during the degradation of CO.
The FTIR results revealed that the inherent organic compounds in GL were phosphorous and nitrogen. The effects of process parameters such as temperature and initial CO concentration showed that the optimum conditions for the CO degradation process were found at 30°C and 90 gL −1 , respectively. The one-way ANOVA and Turkey's post-hoc analysis revealed that the CO initial concentrations of 70 and 90 gL −1 were statistically significant (P < 0.05). Higher CO degradation rate constants with corresponding lower biological half-lives were obtained using the first-order kinetic model. The CO modeling data according to the linear regression suggested that the first-order kinetic model gave a better fit. Therefore, based on the findings of this present study, it may be concluded that GL could be used as an effective organic treatment for COCS, while the firstorder kinetic model was suitable for modeling the CO degradation.