Kinetic modelling of in situ treatment of petroleum hydrocarbon contaminated soil using bone char and NPK fertilizers

This study investigates the effectiveness of bone char (organic) and inorganic nitrogen, phosphorous, and potassium (NPK) fertilizers as stimulants in the degradability of petroleum hydrocarbon contaminants on soil. The physicochemical properties of the hydrocarbon sludge were used to assess the effectiveness of this process over an 8-week period using 0.5–3.5 kg mass of each fertilizer at different experiments. A first order kinetic model was used to estimate the rate of degradation of the total hydrocarbon content (THC) and total organic carbon (TOC) contaminants and the half-life of the remediation process. The microbial population within the period was also determined. The p-value (p < 0.05) indicated that these fertilizers were effective in degrading these contaminants on the soil, because of the significant difference between the treated and the control soil samples. A direct relationship was observed between with the mass and performance of the fertilizers. With 3.5 kg mass of the fertilizers, rate constants of 0.018 and 0.019 d− 1 were obtained for the removal of the THC and TOC contaminants, respectively, using the bone char fertilizer, whereas NPK fertilizer gave rates of 0.03 and 0.023 d− 1 respectively. The performance of the NPK fertilizer is attributable to its nitrogen and phosphorous content. The model adequately described the process and showed the effectiveness of both fertilizers in the remediation process.


Introduction
Crude oil exploration and exploitation has given rise to the problem of environmental pollution in many parts of the world. The large-scale production and oil exploration activities in the Niger Delta region of Nigeria, renowned for its huge deposit of crude oil and natural gas reserve, has given rise to oil spillage which has impacted negatively on the environment leading to air, water, and land pollution. The effects on terrestrial and aquatic lives are equally huge and devastating, and these consequences include surface layer poisoning, distortion of macro and micro ecology, as well as other environmental degradation and health impacts. The land desolation, soil and water body degradation associated with these spills lead to reduced agricultural production [1].
To this end, the need for remediation of these soils has become vital. Bioremediation has been shown to be an efficient and viable approach in the management of oil spilled soils. This technique involves the use of microorganisms such as bacteria, yeast, and fungi to degrade toxic and nontoxic substances. Hence, it involves the breakdown of these substances using living organisms [2,3]. Various techniques of this approach exist and these include: bio-stimulation, bioaugmentation, bioventing, land farming, phytoremediation, intrinsic bioremediation and controlled natural attenuation, and they could be applied in-situ or ex-situ depending on other factors of consideration [4][5][6][7][8].
This study which considered the bio-stimulation approach involved the enhancement of the hydrocarbon degradation process through microbial activities which were stimulated through the addition of nutrients. Bundy et al. [3] and Chorom et al. [9] reported that water, temperature, pH, oxygen content of the soil as well as availability of the contaminants to microorganisms could affect the degradation rate. The most widely reported application of this technique involves the use of nitrogen, phosphorous, and potassium (NPK) fertilizers. Increase in soil pH, mineral leaching, and other environmental effects are among the adverse effects reported with the use of inorganic fertilizers, in addition to the huge cost consideration which remains a significant factor of consideration with these materials. On the other hand, organic fertilizers have been reported to fare better in these areas [10]. This study investigates the effectiveness of bone char (organic) and NPK (inorganic) fertilizers as stimulants in the remediation of an oil contaminated site with a view to determining the kinetics of the process.

Study area
An experimental site located within the University of Port Harcourt (4.9069°N, 6.9170°E), Rivers State-Nigeria was used for this study after due authorization by the University research and ethics committee. The area is a tropical rainforest region, which has an average daily minimum temperature of 23°C and a mean temperature of 31.5°C. It has about seven (7) months of rainfall annually (April to October), an average annual rainfall of 2400 mL, and a mean relative humidity of 85% per month.

Bone char fertilizer
Bone char sourced from an abattoir within the study area was subjected to pyrolysis treatment at the Chemical Engineering laboratory, University of Port Harcourt, to obtain an organic fertilizer.

Crude oil sample
The crude oil sample used for this experiment was sourced from an oil spill site at Abua/Odual local government area, Rivers state-Nigeria using a plastic container. 200 L of the crude sample was collected from recovery tanks installed by the clean-up contractor. The sample was collected under the supervision and assessment of department of petroleum resources to ensure conformity and compliance with product and environment standards. Other chemicals used in this experiment were purchased from a retail outlet, Austino laboratory services ltd. Choba, Port Harcourt (Nigeria) and were of analytical grade.

Soil preparation and in-situ bioremediation procedure
The land portion, with no known history of crude oil contamination, was segmented into 7 cells of similar dimension demarcated from each other with a 10 cm wooden structure. Each cell, having a dimension of 1.5 × 1.5 m and a depth of 30 cm (contaminated soil depth), was polluted with 20 L of the crude oil sample to form a sludge which was allowed to stand for 2 wk. under continuous monitoring. At the end of this period, the vegetations on the study site were observed to have died. Thereafter, different concentrations of the bone char and NPK fertilizers were administered on these cells (Table 1), and the soil tilled to expose its surface area for treatment and to enhance oxygen availability to bacteria and fungi for effective chemical reactions. Tilling of the soil was done once in every 3 days whereas a daily irrigation of the site helped ensure that the optimal soil moisture condition of the site was maintained. Moreover, the experiment was conducted in an uncontrolled environment during the dry season. The effect of bioremediation on the contaminated soil was observed for a period of 8 wk. Previous studies by Daka et al. [11] and Ayotamuno et al. [12] employed a similar approach under different experimental conditions.
The analytical soil samples were collected at depths of 6, 18, and 27 cm from the cells at the experimental site using a spade. The soil samples were, thereafter, mixed into composite samples, labeled and sealed in small polyethylene bags for analysis.

Characterization of soil sample
Physicochemical properties such as pH, total hydrocarbon content (THC), organic carbon, extractable phosphorous and total nitrogen content were used as indices to evaluate the contamination levels of the soil. Samples collected from contaminated site were air dried and homogenized by passing them through a 2 mm mesh sieve. The physiochemical properties were determined as follows: 2.4.1 THC 10 g of the representative soil sample was shaken vigorously with 10 mL of n-hexane for 10 min. The extract was analyzed with 420 nm Bausch & Lamb spectronic-70 spectrophotometer (Bamko-surplus, Texas, United States). The concentration of the extract was obtained from the spectrophotometer and compared with a standard curve of another substance with known concentration. The THC was evaluated with reference to the American Society for Testing and Materials D 9071 procedure and multiplied by the appropriate dilution factor in line with the study of Ikiogha et al. [13].

Total nitrogen (TN)
2 g of the soil sample was air died, weighed and added into a 150 mL beaker. 10 mL each of the following were added alongside: selenium, copper sulphate, sodium sulphate, and sulphuric acid (H 2 SO 4 ). These reagents were of analytical grade. The mixture was agitated by gentle swirling and digested in a Foss Tecator digester (Kjedahl™ Model 2300, Pittsburgh, USA) for 15 min until an ash coloured solution was formed. The setup was heated for 1 h and cooled for about same time. The digest was emptied into a 250 mL capacity beaker and made up with distilled water. An auto-analyzer was used to determine the ammonium nitrogen of the aliquots. The percentage nitrogen content was thereafter evaluated with due consideration of the dilution factors.

pH
This was done using a corning pH meter model 7 (Corning Scientific, NY, USA) by dipping the electrode into a 1:25 sample and water suspension which had been stirred and allowed to equilibrate for about 1 h.

Phosphorous content
The phosphorus content was determined using the Bray and Kurtz method [14]. 2.85 g of sample was weighed into a clean and dry sample bottle. 20 mL of Bray No. 1 extractant was added (0.025 N HCl + 0.03 N NH 4 F) and shaken for 60 s. Thereafter, 10 mL filtrate was added into a 50 mL capacity volumetric flask using a pipette. The mixture was diluted (made up) with about 20 mL of distilled water, after which 4 mL ascorbic acid solution was added. This setup was allowed to rest for 30 min, during which a change in colour was observed, before being analyzed with a spectronic-20 spectrophotometer operating at 660 nm wavelength. This analysis was carried out in accordance with the manufacturer's specification.

Total organic carbon (TOC)
This was evaluated using the wet combustion method by measuring 1 g of the finely ground soil sample into 3 different 500 mL beakers with 10 mL of potassium dichromate solution added alongside amidst continuous stirring using a magnetic stirrer. 20 mL H 2 SO 4 was added to the content of the beakers using a graduated cylinder and, thereafter, the beaker was rotated and allowed to stand for 10 min while distilled water was equally added for dilution, by making up the solution to the 200 mL mark of the beaker. 25 mL of 10.5 N ferrous ammonium sulphate and 0.4 N potassium permanganate were also added under strong light.

Soil texture
This was evaluated using the Bouyoucous hydrometer [15] in accordance with the method of Day [16].

Evaluation of microbial population
This process was carried out using the methods of Barrow and Feltham [17] and Garrity et al. [18]. Prior to the experimental process, all the media and diluents were first sterilized in an autoclave operating at a temperature condition of 120°C for 15 min. The glass wares were dried in a hot air oven for 3 min. The saline solution was prepared by weighing 0.85 sodium chloride solution which was thereafter transferred into the various 9 mL capacity test tubes used for the experiments for sterilization under a temperature condition of 120°C for 15 min. The solution was also used for the serial dilution of the soil samples. Cultivation and enumeration of the soil bacteria was done by adding 1 g of the air-dried soil sample into 1 mL of the sterile distilled water amidst vigorous agitation to obtain a homogeneous mixture. 1 mL of the 10 − 1 solution was added into a 9 mL test tube containing the diluents and serially diluted stepwise to obtain a 10 − 3 dilution. A spatula was used to transfer 0.1 mL aliquot of the solution into the sterilized nutrient agar plates in such a way that most of the surface area was exposed for reaction. The inoculated plates were incubated at a temperature of 37°C for 24 h, before the nutrient agar plates were examined for bacterial growth. The total viable heterotrophic bacteria were evaluated based on the number of the Colony Forming Units (CFU) using the colony counting technique to determine the number of cells capable of splitting [19]. The CFU was evaluated using the correlation: CFU mL − 1 = number of colonies × dilution factor/volume of culture plate.

Identification and isolation of soil bacteria
The bacteria culture was isolated by aseptically streaking colonies of different culture seen on the plate to a freshly prepared nutrient agar plate. The setup was incubated at 28°C and allowed to stand for 24 h to allow the bacteria grow. Based on the colony morphology and exhibition of zones of inhibition on the primary culture, isolates were selected. Isolate purification was done using dilute nutrient broth and agar media and stocked at − 80°C. Characterization of these isolates was done using Gram stain (methyl and violet test) to confirm the bacterial growth while the viable count was conducted manually. The isolates were observed under an inverted microscope at × 100 oil immersion following the method of Khan et al. [20] while the CFU was evaluated.

Isolation and identification of soil fungi
This was done using dilution plate method on potato dextrose agar. Isolates of soil fungi which were examined macroscopically and microscopically were identified by matching the characteristics of the fungus obtained from observations in accordance with the work of Umboh et al. [21].

Kinetics of the process
A description of the process was evaluated using a firstorder kinetic model in the form of Eq. (1). where: FE = fractional efficiency. C o = initial concentration. C t = concentration at time t. t = elapsed time.

Degradation
Degradation of the contaminants was evaluated from the initial and residual concentrations of the THC and TOC, respectively using Eq. (3).

Bio-stimulation efficiency
The bio-stimulation efficiency was used to estimate all the treatability options provided by the optimization rates of the bio-stimulants and was obtained using Eq. (4). This was done to ascertain the effect of the biotic and abiotic factors in the bioremediation process.
where: BE = bio-stimulation efficiency. Ct = percentage of crude oil removed from the treated soil.
Cu = percentage of crude oil removed from the untreated soil (control).

Statistical analysis
t-test was used to analyze the obtained THC data to evaluate the bioremediation process on the contaminated soil. The choice of the test instrument was based on the sample size (n < 30). The treated soil samples were compared with that of the control site and the following conditional statements were considered: i. If p > 0.05 accept null hypothesis (H o ) (no appreciable effect on the process) and reject alternative hypothesis (H i ) (significant effect on the process). ii. If p < 0.05 reject null hypothesis (H o ) (no appreciable effect on the process) and we accept alternative hypothesis (H i ) (significant effect on the process).
The Pearson correlation coefficient and t-stat values were also obtained. Microsoft excel and SPSS version 17 were used for this process.

Soil properties
The particle size analysis of the soil texture within the depth of interest (30 cm) prior to treatment showed it to comprise of the following: soil-75%, clay-15% and silt-10%. The soil texture, which is a classification of the soil based on its physical characteristics, was used to assess the soil sample to ascertain the air and water retention capability, rate of water movement in and out of the soil, as well as the ease of working with the soil sample. The texture triangle (Fig. 1) indicates the soil texture to comprise of both sand and loamy sand and is consistent with the work of Odokuma and Dickson [22] for the top soil of tropical rainforest. The physicochemical properties and microbial properties of the sandy loam soil ( Table 2) were adversely affected by the contamination of the petroleum hydrocarbon. Both the trace metal contents and the microbial population of the soil were greatly depleted. The pH conditions in this study were optimum and favourable for microbial growth and higher than the minimum growth pH from previous findings [23]. Slight changes in the pH could alter the ionization of aminoacid functional groups and interfere with the hydrogen bonding leading to their denaturation thereby rendering them ineffective for use in the biodegradation of the hydrocarbon [24]. This study corroborates previous findings on the effect of petroleum hydrocarbon contamination on soil moisture [25]. Increase in the electrical conductivity and cation exchange capacity (CEC) indicates the dissociation of ion charges into the soil due to crude oil contamination. These ions are attracted to the water molecules in the soil which carry polar charges. The ionic solutions become electrolytes capable of conducting electricity.

Fertilizer properties
The physiochemical properties of the bone char and NPK fertilizers (Table 3) are presented. The constituents of these fertilizers and their physical properties are critical to their performance in the degradation of the hydrocarbon contaminant. The effect of the structural-mechanical properties of fertilizer granules has been reported [26] and these are fundamental in the exposure of the reactant surface area necessary for these reactions to occur as well as for the interaction of the microbes with the substrate. Various authors have shown the effect of these fertilizers on growth characteristics and soil nutrients [27][28][29][30] and bioremediation [31][32][33][34][35]. This study corroborates these findings.

Effect of in-situ treatment of the contaminated soil
The effect of both grades of fertilizers (organic and inorganic) on the experimental site was evaluated through an analysis of the physicochemical properties of the soil samples (Table 2) and the THC (Fig. 2a) of the soil.
Both grades of fertilizer drastically reduced the THC of the contaminated soil over time. This indicates that bioremediation was effective in reducing the THC from the contaminated soil. There is a link between the mass of fertilizer and performance in the bioremediation process for both grades of fertilizers. 3.5 kg mass of both fertilizers had the highest reduction of the THC and TOC contaminants; efficiency for all treated samples on the THC was above 60% (Fig. 3a). Comparatively, the NPK grade gave a better performance than the bone char. The effectiveness of the bone char grade of fertilizer could be linked with its ion content leading to the degradation of the hydrocarbon contaminant [36]. A similar pattern was observed for the TOC (Fig. 2b and  3b) within this period which confirms the bioremediation effects of these grades of fertilizers.
The bacillus and pseudomonas species were predominant on the contaminated soil. Addition of these fertilizers resulted in a significant increment in the heterotrophic bacteria and heterotrophic fungi. Total Bacterial Count (TBC) was used to obtain a quantitative idea of the microorganisms present in the soil. TBC is indicative of the number of CFU of the soil sample and formed the basis to analyze the degradation rate of the crude oil in the polluted soil. Prolonged period of the bioremediation treatment with both grades of fertilizers favoured the production of the microbes (Fig. 4a). The NPK fertilizer grade gave higher bacterial count than the bone char whereas higher mass of both fertilizers favoured the TBC and consequently their performance in the remediation process. Chorom et al. [9] reported that excessive application of fertilizers makes survival conditions of these microbes critical. Hydrocarbon Utilizing Bacteria (HUB) degrades aliphatic alkanes through various metabolic activities. They utilize the linear hydrocarbon chains in petroleum as their main energy source under aerobic conditions. Additionally, when exposed with other sufficient limiting nutrients like nitrogen and phosphorous, they produce surfactant glycolipids which help reduce water tension and the uptake of hydrocarbon, e.g., Alcanivorax borkumensis. Performance of the both grades of fertilizer in the isolation of the HUB was shown to be time and dose-dependent (Fig. 4b). The higher performance of the NPK fertilizer in the remediation process could be attributed to its higher nitrogen and phosphorous content which makes it possible to isolate more of the desirable microbe population to produce the surfactant glycol lipids that reduce surface tension thereby making it possible to degrade the contaminant hydrocarbon. The biostimulation efficiency (Fig. 5) further corroborates this.
Similarly, hydrocarbon utilizing fungi (HUF) has proven to be very effective in the removal of a wide array of toxins or contaminants from the environment. This is achieved through the degradation of the hydrocarbon components to non-toxic substances. The growth rate of the HUF (Fig. 4c) followed a similar pattern as that of the HUB, and the NPK fertilizer favoured the production of the HUF as it consistently gave a higher count of the HUF than the bone char within the period of study. Performance of the both grades of fertilizer in the isolation of the HUF has been shown to be time and dosedependent. Kumar et al. [37] and Chikere et al. [38] reported the species of this as being capable of degrading petroleum hydrocarbons in contaminated soils: Penicillium sp., Candida sp., A. Niger, Rhizopus sp., Mucor sp., Talaromyces sp., Saccharomyces sp., Cladosporium sp., Fusarium sp., Rhodotorula sp. and Aspergillus flavus. The steady increase in the growth rate of these microorganisms with time entails a rapid reproduction of the crude oil degrading mechanisms which were responsible for the attenuation of the soil contaminants and this corroborates the works of Chikere et al. [31].

Statistical analysis
The results of the statistical analysis are shown ( Table 4). The R 2 , t-stat and p-values were used to study the effectiveness of the bioremediation treatment in the removal of the THC and TOC contents. The analysis was done in comparison with the control sample based on the pvalue for significance level (p < 0.05). Appreciable significant differences were observed for all the samples. Hence we reject the null hypothesis (H 0 ) and accept the alternative hypothesis (H i ). There was significant effect of the remediation treatment process on the contaminated soil which gave rise to an appreciable reduction of THC and TOC.

Kinetics of the process
The first order kinetic model, which proposes that rate of degradation of substrate is directly proportional to its concentration, was used to model the degradation level of THC and TOC in the contaminated soil. High k value implies high degradation rate (Table 5), while the fractional efficiency was plotted against time for all samples ( Fig. 6a-f).
Positive correlations were observed for the degradation of THC and TOC in the soil with respect to time. In both cases, 3.5 kg NPK treatment had the highest degradation rate constant (0.03 and 0.023 d − 1 ) and   consequently half-life of 23 and 30 d respectively. This is further proof of the effectiveness of this condition in the bioremediation process. The high values of the R 2 statistic across most of the treatment conditions indicate the sufficiency of the model in predicting these processes. It also shows that the models adequately describe the variability of the removal of the THC and TOC contaminants to the extent of their correlation. Physical inspection of the site indicates tremendous effect of the remediation process on the contaminated soil (Fig. 7). The general trend in the THC degradation profile was observed to follow similar trend reported by Ofoegbu et al. [39]. The kinetics of the remediation process also corroborated the physical evidence of the process (Fig.  7). The kinetic studies of the first order constant obtained in this study, using NPK as biostimulant for the remediation of the soil corroborates the range of 0.007-0.042 d − 1 earlier reported [39]. Results from this study also show the significant effect and contribution of bone char as a bio-stimulant in the remediation of the hydrocarbon contaminated site. These effects were comparable to those of the NPK grade which indicate their suitability as a viable alternative. More so, the first order kinetics is comparable to that obtained for the NPK grade. Although the higher performance of the NPK grade fertilizer has been attributed to its nitrogen and phosphorous content, leaching of the soil nutrients and loss of soil texture are some of the side effects associated with its use in addition to the risk of toxicity to aquatic life. With both grades of fertilizer, the bio-stimulation efficiency was above 70% for both THC and TOC in crude oil contaminated soil, whereas the degradation was above 60% which is a strong indication of the effectiveness of the process. The R 2 value obtained from the kinetic models indicates the reliability of the models for future predictions and the attenuation of the samples. Based on the indices considered in this study, it can be deduced that enrichment of the bone char fertilizer can make it a credible alternative for NPK (inorganic fertilizer).

Conclusions
This work has successfully shown the efficacy of the bone char (organic) and NPK (inorganic) fertilizers in the remediation of hydrocarbon contaminated soil. Vital information on the kinetics of the process, degradation, and biostimulation efficiency has been evaluated. Although the inorganic fertilizer improved the contaminated soil, this study has shown that this organic fertilizer can provide a viable alternative. The amount of fertilizer had a direct effect on soil improvement and microbial growth. This gave rise to an appreciable reduction of the THC and TOC. Enrichment of the organic fertilizer will boost its capacity for improving the recovery of oil spilled sites across this region.