Skip to main content

Factorial design study of total petroleum contaminated soil treatment using land farming technique


Land farming technique was used to treat hydrocarbon contaminated soil collected from a crude oil spill sites in Edo State, Nigeria. Calibrated standard auger was used to collect soil samples from the site at depth below 30 cm. The samples were characterized and classified. Cow dung and NPK fertilizer were added as additives to complement the nutriments of the soil samples before total petroleum hydrocarbon (TPH) quantification and remediation procedures. Factorial design was applied to vary the input parameters such as pH, mass of substrate, moisture content and turning times of land farming so to ascertain the optimal conditions for the procedure. The result revealed that the in-situ TPH value was 5000 mg kg− 1 on the average and after 90 d of treatment, TPH reduced to 646 mg kg− 1. The turning rate, pH, moisture content and mass of substrate hade 83, 4.36, 0.48 and 0.046% contribution, respectively, for the degradation process using land farming treatment. Numerical optimization techniques applied in the optimum point for land farming input parameters to achieve predicted maximum removal of 99% were evaluated as pH, mass of substrate, moisture content and turning rate to be 6.01, 1 kg, 10% and 5 times in a week, respectively. TPH removed at this optimum point was 98% reducing from 5000 to 636 mg kg− 1. The high coefficient of determination (r2 = 0.9865) as observed in the closeness of predicted and experimental values reflects the reliability of the model and hence, land farming practice with close attention on turning rate as revealed by this study, is recommended for TPH contaminated soil remediation.


Advancement in technology, continuous urban sprawling and improved standard of living have over the years, caused a corresponding increase on energy demand, which is largely used in powering automobile and other related machines and appliances. Energy from coal, fossil fuel and some renewable sources like solar and biomass have been widely used with fossil fuel being the most utilized among them [1, 2]. Fuel is one of the major products of processed crude oil, rich in hydrocarbon content and is largely sort after for effective running of human daily activities. These activities have in one way or the other hampered the chain procedures of crude oil - drilling, refining, treatment, transportation and utilization which on the long run result in spillage, thus distorting several ecosystems and rendering most lands useless [3]. In Nigeria (especially the Niger Delta region – comprising of nine states), there has been oil spills resulting in soil contamination due to poor operation and management practices [4, 5]. It is reported that about 13 Mt of hydrocarbons are spilled which is caused largely by pipeline vandalism, destructive crude oil theft, operational spills and engineering failure (such as pipeline rupture), and uncivilized refining conditions [6,7,8,9,10,11]. The severity of damage done to these soils by hydrocarbon spill is a function of diverse factors such as partition coefficient of the soil, permeability, absorption properties and chemical constituents of the hydrocarbon. Another source through which spillage occurs is through natural seeping in locations where hydrocarbon is found in sub-surface deposit to accidental discharge of crude oil onto ground surface and several other points of pollution, but irrespective of this source, once hydrocarbon spills into the soil, it alters both its physical and chemical properties [12,13,14], thus becoming harmful to plants, microorganism in the soils and humans.

Effective cleaning-up oil-contaminated soils by adopting some available technologies, is a viable option of remediation process and this is done to degrade hydrocarbon present in the soil. Hydrocarbon degradation is a process that involves the gradual weathering and removal of petroleum constituents especially the non-volatile compounds from the contaminated location by using physical, chemical and biological methods for remediation of contaminated soils [15,16,17]. For instance, bioremediation which involves the utilization of effective microbes for hydrocarbon degradation has increasingly gained researchers interest in recent decades. The most frequently isolated and utilized hydrocarbon degrading microbes are genus Pseudomonas which degrade complex chains of hydrocarbon into smaller and less toxic compound. Also, fungi in the genera of Fusarium, Rhizopus and Penicillium have gained acceptance in treating hydrocarbon contaminated soil since Exxon Valdez spillage in 1980 [17, 18]. Land farming has been acknowledged as an effective and low-cost technology for abstraction of total petroleum hydrocarbons (TPHs) from soil [18,19,20]. It is reckoned to use less energy and it is not harmful to the environment, with reduced residue disposal problems [21]. Land farming treatment is the application of calculated organic and inorganic substrates on contaminated soil in order to mineralize the toxic substances in the soil [22, 23]. Land farming is a concept that entails nutriments addition and replication of microbes, geared towards increasing the number and growth of microorganisms in order to accelerate bioremediation rate [20, 24, 25]. As microbes require sufficient major element like carbon, hydrogen, oxygen, nitrogen, and phosphorous for the development of macromolecules, fertilizer addition provides the bacterial with vital elements to thrive and reproduce. In some cases, sawdusts, animal dungs, and straws may supply bacterial with carbon sources [26]. Land farming techniques has been practiced in some regions of the world to bioremediate crude oil contamination in soils to minimize the health risk on human and the environment at large [27, 28]. It has been used successfully to remove petroleum hydrocarbons at large scale [23, 29], and because of its simplicity in implementation, Niger Delta has also employed it. Unfortunately, with the handful of its application within Nigeria, there is still dearth of information on the efficient practice of land farming treatment for crude oil contaminated soils for effective remediation.

The effectiveness of land farming can be enhanced when environmental circumstances allow the growth of microbes, and this depends on some certain environmental parameters such as pH, moisture content, nutrient availability, among others [23, 30]. Factorial design (FD) is normally used in screening variables (both dependent and independent) and also in optimizing response surfaces. The latter is frequently used for experimental designs involving experimental procedures [31]. FD has been employed in some oil biodegradation studies of constituent’s optimization that may induce the microbial debasement phenomenon hereby contributing to the progress of oil spill bioremediation process. Bhattacharya and Biswas [32] investigated the effect of various nutrients added to waste engine oil biodegradation of Bushnell-Haas medium using Ochrobactrum pseudintermedium bacterium. The data permit the development of an empirical model (p <  0.00672) through the application of a full FD for experimental work thus, describing the connection between dependent and independent variables. Jasmine and Mukherji [33] also assessed the treatment of refinery oily sludge using 2n full FD via bioaugmentation and biostimulation processes. FD was also applied in the bioremediation of artificially contaminated soil with weathered bonny light crude oil (WBLCO) using biostimulation and bioaugmentation processes. A statistically significant (p <  0.0001) second-order regression model with a coefficient of determination (R = 0.9996) was ultimately obtained for removal of WBLCO. Numerical optimization process was also carried out based on desirability function to optimize the bioremediation process [34]. Further researches are ongoing to develop and improve on FD methods for minimizing the experiment number and the interactions of their input variables/parameters. This has been achieved by utilizing design experiment procedure to generate information on direct effects, interactive pair effects and effects due to curvilinear variables. Some ample studies have been done on the application of FD in bioremediation of soil contamination using bioaugmentation and biostimulation techniques as presented above. From the resources available and accessible and to the utmost best of our knowledge, there are limited or no information on the optimization of land farming procedure using FD study which plays a major role in the adequate treatment of hydrocarbon contaminated soils. In this study, FD was applied to vary the input parameters such as pH, mass of substrate and moisture content in order to optimize them for best hydrocarbon removal.

Materials and methods

Site location

The site selected for this project is an oil field located in Ologbo community, Ikpoba Okha Local Government Area of Edo state in Southern Nigeria. Edo state is bounded to the right by Ondo State and to the lower left by Delta (Fig. 1). Ologbo as a major community is one of the oil producing area with multiple petroleum production facilities in Niger-Delta area of Nigeria. The community houses a gas plant operated by the Nigerian National Petroleum Development Corporation (a subsidiary of Nigerian National Petroleum Corporation) and some other petroleum facilities. It lies between longitude 05038′36.44″ E to 0504′26.56″ E and latitude 06004′28.17″ N to 06004′33.79″ N. It is about 32 km away from the south-western part of Benin-City and over 30 km from Nigeria National Petroleum Development Corporation access road, which is off Benin-Sapele highway. Within this location, crude oil spillage is frequent resulting from vandalism and sabotaging of oil pipes and equipment by militants and oil pilferers, thus leaving the land degraded and contaminated. Figures 1 and 2 give the location map of the study area and one of the contaminated spots in the study area, respectively.

Fig. 1
figure 1

Map of project location

Fig. 2
figure 2

Contaminated area in the study location

Preliminary investigation and TPHs quantification procedure

As a vital step towards a successful remediation process, reconnaissance survey was carried out on the contaminated site in order to minimize challenges during sample collection. A calibrated standard auger was used to collect samples at surficial depth not exceeding 30 cm, the samples were sun dried and homogenized (using mortar and pestle) before sieving through a 4 mm sieve. The homogenized samples were store in polythene bags at room temperature to prevent moisturizing. The soil samples were characterized so as to determine its physical, chemical and microbes’ constituents using British Standards BS 5930 (Table 1). The constituents of the soil are seen to fall below the recommended nutrients required for effective biodegradation process. Therefore, NPK fertilizer in ratio 20:10:10 and cow dung was added as additives to complement these nutriments for the remediation procedure. These fertilizers (organic and inorganic) used, have high nitrogen content which makes them suitable for remediation operations. Their compositions are shown in Table 2. Fresh samples from the contaminated site were taken to the laboratory for residual TPH quantification in accordance with USEPA [35] and ASTM [36]. TPHs were extracted from the samples by drying and passing them through a 4 mm sieve aperture size. The samples were placed in 40 mL centrifuge bottle with 25 mL of chloroform added. The samples were tightly closed and kept well in a sonicator bath for 60 min. During the process of extraction, iced deionized water was continuously added to maintain a temperature below 40 °C. On completion of extraction, samples were subjected again to centrifugal force for 11 min at 3000 rpm. The resultant extract was then placed in an Erlenmeyers flask where it was dried to achieve a specific weight. Bathing was done at 65 °C to evaporate volatile chloroform and the extract shows an average contamination concentration of 5000 mg kg− 1. This is equivalent to intervention level according to USEPA, hence the need to remediate the contaminated soil.

Table 1 Physical and chemical properties of the soil used
Table 2 Properties of NPK fertilizer and cow dung used

Experimental design and procedure

In the initiation of the treatment, 100 kg of sieve samples were placed in twenty buckets and labeled based on the treatment to be accommodated in the setup in accordance with USEPA [35] and ASTM [36]. The choice of input variables, range of variables and duration of the experiment as stipulated in the USEPA procedure were adopted for this study. Four major input variables were selected namely pH, moisture content, mass of substrate and turning rate and were varied in each of the buckets. Substrate used was cow dung and NPK fertilizer, with its application ranged from 0.6 to 1 kg. In every application, mass of substrate constitutes 50% of cow dung and 50% of NPK fertilizer in any experimental run to make up the total mass of substrate required. The pH and moisture content of each experimental run was adjusted to reflect the value to be used for that particular setup. The pH was adjusted using slaked lime and measured using pH meter while the moisture content was a percentage of the weight of each experimental setup. The batch to batch variation was controlled using the range of input variables presented in Tables 3 and 4. The sorption of hydrocarbon from the soil was carried out using laboratory examination in order to feasibly select factors controlling the biosorption process in land farming treatment. To be able to select the input variables with the highest significant contributions to the remediation process and determine their optimum values, FD of experiment was used for screening. The range and levels of the input variables used in designing the experiment is presented in Table 3. Runs 17–20 were used as control for the study and the treatment was carried out for 90 d after which the samples were taken from each bucket for residual TPH determination. According to USEPA [35] and ASTM [36], FD study of this nature with experimental setup of 2n + 1 < 100, should have four middle values (control) with the same input variables, hence runs 17–20 were designated as control while all the input variable had the same range of values as shown in Table 4. Petroleum degrading bacteria was enumerated through Mineral Salt Agar culture following the procedures of Sepahi [37].

Table 3 Ranges and levels of selected input variables
Table 4 Response of TPH on factorial design of experiment used for variable screening

Statistical analysis

The data obtained from the experimental procedures were statistically analyzed using Excel (Microsoft office product version 16), Design-expert and STATISTICA software. The suitability of the FD to screen the variables was carried out by computing the standard error, correlation matrix of regression coefficient and model leverages. Analysis of variance (ANOVA) and goodness of fit were also computed to validate the model significance. The major effects of the four-treatment variable as well as the interactions were interpreted jointly. In every 22 factorial designs, the F-tests is enough to reveal the interrelation in combined treatment procedures. It also tells the relationship between all the variables concentrations in the treatment parameters. The result reveals the main variable with the largest effect in the four combined parameters by comparing the means. The F-test procedures employed are shown in Eqs. (1), (2) and (3) respectively.

$$ {F}_{presentation}=\frac{SS_{presentation}/{df}_{presentation}}{SS_{error}/{df}_{error}} $$
$$ {F}_{difficulty}=\frac{SS_{difficulty}/{df}_{difficulty}}{SS_{error}/{df}_{error}} $$
$$ {F}_{interac\mathrm{t} ion}=\frac{SS_{interac tion}/{df}_{interac tion}}{SS_{error}/{df}_{error}} $$

where Fpresentation = main effect due to presentation, Fdifficulty = main effect due to difficulty, and Finteraction = main effects due to interaction.

Results and discussion

TPH biodegradation

The degradation for standard run 17–20 which serve as control for the study was similar having the same combined variables. TPH concentration degraded from 5000 mg kg− 1 to between 722 and 862 mg kg− 1 in 90 days duration while maintaining 3 d wk.− 1 turning rate. Standard runs 3, 4, 7, 8, 11, 12, 15, and 16; with moisture content of 50% the set-ups had hydrocarbon content floating on the surface. This made the admix semifluid in nature and easy to turn using hand trowel. When turning rate is effective and properly practice, the hydrocarbon contaminants becoming exposed to degrading agents and are therefore either degraded or mineralized [16]. This is attributed to the over 80% TPH reduction recorded in standard run 16 after 90 d treatment.

Substrate addition also enhanced TPH degradation as it serves as energizer for the microbes. This served mainly as catalyst in microbial reproduction processes and consequent consumption of the TPH contaminants. Although higher concentration of substrate does not guarantee high TPH degradation, but when suitably combined with other parameters such as turning rate, high pH value and average moisture content; then a better degradation result can be obtained [27]. As fertilizer application on crude oil contaminated site was well systemic and well calculated, TPH degradation result was less than 40%, this was mainly due to the low moisture content and high acidic values of the treated samples. Nwilo and Badejo [11] had similar results from their study in which NPK fertilizer was used in the treatment of soil collected from a spill site. TPH degradation was faster in samples with lower moisture content than samples with higher moisture content. The pH of the contaminated soil samples before treatment ranged from 2 to 5. An increase in pH values was observed as the treatment progress into day 15 to 70, the pH value ranged from 5.7 to 7.1 (neutral). The addition of fertilizer to hydrocarbon polluted soil samples had a catalyst effect on the treatment and the pH value increase from acidic range of 2 to neutral range of 6.8. The substrate applied caused an increase in the total nitrogenous content of the soil but as the treatment days increased from day 50–60, the nitrogen content decreased gradually. This could be linked to the soil bacterial consuming the nitrogen for the hydrocarbon degradation, thus reducing the available nitrogen as treatment time increases [9]. This process utilizes biochemical reduction and it is initiated by denitrifying bacterial in the soil [11]. In all the factorial setups for the hydrocarbon contamination treatment, there was significant TPH degradation and the bacterial population in all the setups increased exponentially. The petroleum degrading bacteria increased from 1.8E+ 01 to 3.6E+ 08 cfu g− 1 during the treatment period. This increase confirms the loss of nitrogen which usually accompany degradation procedures [11, 17]. This increment in petroleum degrading bacteria is in tandem with the findings of Oluwatuyi et al. [12] and Okonofua et al. [13]. FD analysis of results was then employed to determine the variable with the most significant contribution in the TPH degradation procedure.

FD of experiment

The response of TPHs on FD of experiment, used for variable screening hydrocarbon contaminant concentration of 5000 mg kg− 1 within a period of 90 d is presented in Table 4. The minimum value of TPH is given as 450 mg kg− 1 while the maximum value is 1393 mg kg− 1. The calculated mean value is 921 mg kg− 1 with the standard deviation of 302 mg kg− 1. In assessing the worthiness of FD in screening the input variables based on their fundamental and important contributions, model standard error analysis was used based on Montgomery [38]. Presented in Table 5 are the computed standard errors for the chosen response. From the result, a low model standard error of 0.25 was achieved for both the individual and combine terms and effects. According to Jasmine and Mukherji [33], standard errors must be akin within a coefficient and the lower the value is, the better. Similarly, the error values were lower than the model basic standard deviation (SD) of 1.0 suggesting that the FD was perfect for the screening process. To demonstrate for multicollinearity, the variance inflation factor (VIF) of the analysis was obtained all through as 1.0 representing a superb outcome as a perfect VIF should give 1.0. VIFs closer to 10 or greater are usually cause for concern, and this signifies that coefficients are basely calculated due to multicollinearity [39].

Table 5 Result of computed standard errors

Furthermore, the Ri-squared values also gives zero which perfectly match an ideal Ri-square as high Ri-squared especially values above 1.0 shows that design terms are correlated utimately resulting to poor models. Table 6 presents the correlation matrix of the regression coefficient. It can be seen that, off diagonal matrix, the lower values obtained point out the fact that the model is well fitted and it is strengthened enough to pilot the design space thus adequately optimizing the chosen response variable. Also, the model leveages were computed in order to better understand the influencial effect of individual design points on the model’s predicted value. According to Meloun and Militky [40], leverage point indicates the extent of influence of an individual design point on the model’s predicted values and it usually varies from 0 to 1. A leverage of 1 indicates that, the predicted value at a specific case will perfectly equal the observed value of the experiment, making the residual to be zero. The addition of leverage values in all cases equals the number of coefficients fit by the model, and the ultimate leverage an experiment can have is determined by 1/m, with m being the number of rounds the experiment was repeated. Leverage of 0.675 calculated in the factorial point indicates that, there is closeness between the predicted values and the experimental values. Hence, less or low residual value approves the sufficiency of the model.

Table 6 Correlation matrix of regression coefficients

Strength assessment of factorial model

To assess the strength of the factorial model towards an effective screening and optimization of the input variables, based on their significant contributions, one-way ANOVA was done for the response variable (Table 7). This was used to examine if the model is significant or not and to also measure the important contributions of individual variable. From the analysis in Table 7, the Model F-value of 56 connotes that the model is significant owing to the fact that there is only 0.01% probability that a “Model F-Value” with high value could occur due to noise. When the values of “Prob > F” are < 0.05, it indicate that the model terms are significant while values > 0.1 indicate the model terms are not significant [41].

Table 7 Analysis of variance table for input variable screening

Therefore, the terms A, D, AD, BC and CD are all significant model terms. Also, 22 from the “Curvature F-value” means that there exist significant curvature in the design space. This is mostly estimated by the difference between the average of the factorial points and that of the center points, and there is just 0.15% chance that a “Curvature F-value” with high value could occur as noise. Furthermore, 0.60 from the “Lack of Fit F-value” connotes that, it is not significant when compared with the pure error but on the other hand, there is a 71% probability that a “Lack of Fit F-value” could occur due to noise. In Table 8, the goodness of fit statistics was used to formalize the sufficiency of the factorial model regarding its potential to screen the input variables based on their significant contribution. From the statistical analysis, the “Predicted R-Squared” value of 0.9188 is in logical agreement with the “Adj R-Squared” value of 0.9684. According to Singh et al. [42], obtaining an adequate precision shows an adequate signal to noise ratio > 4 as desirable. Thus, the computed ratio of 20 as shown in Table 8 connotes an adequate signal. This model outcome therefore shows that it can be used to pilot the design space and properly screen the input variables while also determine their optimum value.

Table 8 Goodness of fit statistics for validating model significance

Input parameters and generated equation

The significant contributions of each input variables were determined using pareto chart. Pareto chart is a graphical presentation of input variables in order of their ranking. Statistical tool was used to generate Pareto’s chart (Fig. 3) for the selected input variables. The result shows that the variables contributed to the hydrocarbon degradation in varying proportion with turning rate, pH, moisture content and mass of substrate all contributing 83, 4.36, 0.48 and 0.046% respectively. Furthermore, the most fitting equation which depicts both the combine interactions and individual effects of the significant input variables (pH, moisture content, mass of substrate and turning rate) against the mesured response (TPH) is provided based on the coded variables and the actual factors which are shown in Eqs. (4) and (5). Either of these two equations can be used in the estimation of the predicted TPH values which is shown in column 3 of Table 9. The predicted TPH values are then compared with the measured values to obtain the residual and the cook’s distance shown in columns 4 and 9 in Table 9. In FD study, only terms without coefficients (zero coefficient) are left out in TPH evaluation using either coded or actual factors, hence the inclusion of AB and BD.

$$ \mathrm{TPH}=950.88-70.61\mathrm{A}-23.51\mathrm{B}+7.26\mathrm{C}-307.70\mathrm{D}-3.06\mathrm{AB}+30.24\mathrm{AC}-36.78\mathrm{AD}+52.63\mathrm{BC}+3.01\mathrm{B}\mathrm{D}-56.25\mathrm{CD} $$
$$ \mathrm{TPH}=1930.89-65.90\mathrm{A}-11.32\mathrm{B}-541.29\mathrm{C}+29.96\mathrm{D}-0.08\mathrm{AB}+75.59\mathrm{AC}-9.20\mathrm{A}\mathrm{D}+13.16\mathrm{BC}+0.08\mathrm{BD}-140.62\mathrm{CD} $$
Fig. 3
figure 3

Pareto chart showing the ranking of the selected input variables

Table 9 Diagnostics case statistics report of observed and predicted TPH

The symptomatic case statistics show the observed values of the response covariant (TPH) against the predicted values as shown in Table 9. This symptomatic case statistics vividly present a clear and deep understanding into the model strength and the adequacy of the FD model.

Model validation

To further evaluate the accuracy of the prediction and established the appropriateness of FD of experiment, the observed and predicted values of TPH were presnted via a reliability plot as shown in Fig. 4. The r2 = 0.9865 which represents the coefficient of determination was utilized in affirming the eligibility of the FD in reducing the TPH. An adequate statistical analysis output must first be used to check the satisfactoriness level of any model before its acceptance. Thus, to examine the statistical properties of the FD model, the normal probability plot of studentized residual shown in Fig. 5 was used to evaluate the regularity of the calculated residuals. The plot of residuals which represent the standard deviation of actual values based on the predicted values was adopted to ascertain if the residuals (observed–predicted) follow a normal distribution pattern. It was depicted that, the computed residuals are normally and approximately distributed which indicates the degree of satisfaction of the developed model. Furthermore, in the analysis, to determine the availability of a possible outlier, cook’s distance plot was generated (Fig. 6). This cook’s distance is to measure the degree at which the regression can change if the outlier is excluded from the analysis. A particular point having a high distance value relative to the other points can possibly be an outlier and should therefore be investigated [43]. From Fig. 6, the plot shows an upper bound and lower bound of 1.00 and 0.00. Therefore, experimental values below the lower bound (0.00) or above the upper bounds (1.00) are termed as outliers which must be adequately investigated. Fortunately, the data of this analysis are free of possible outliers thus showing forth the adequacy of the experimental data. A 3D surface response plot was also provided to study the effects of combine input variables on the response (Fig. 7). It can be seen that the plot depicts the connection between the input variables (pH and turning rate) and the response variable (TPH) and also provides a comprehensible concept of the factorial model. In addition, the colour of the surface gets darker towards the turning rate which connotes that a higher turning rate leads to a reduction in TPH. This observation is in tandem with the work of Agarry and Ogunleye [34].

Fig. 4
figure 4

Reliability plot of observed versus predicted TPH

Fig. 5
figure 5

Normal probability plot of studentized residuals for TPH

Fig. 6
figure 6

Generated cook’s distance for TPH

Fig. 7
figure 7

Effect of pH and turning rate on TPH

Numerical optimizaton

The numerical optimization was finally done to be sure of the desirability of the absolute model. Design expert was adopted in the numerical optimization phase in order to minimize the TPH level and determine the optimum pH, moisture content, mass of substrate and the turning rate. The numerical optimization interphase presents the objective function (Figure S1 of Supplemental Information) with production of twenty (20) optimal solutions (Table 10). From the analysis, turning rate of 5 times a week, with pH of 6.01, moisture content of 10% and substrate mass of 1 kg will result in a minimum TPH value of 636 with a reliability value of 98.6%. The ramp solution showing the graphical representation of the best solution (Figure S2) while the desirability chart depicting the veracity with which the model can predict the values of the chosen input variables and the similar response is presented in Fig. 8. From the outcome on the chart, it can be inferred that the developed and optimized model using FD and numerical optimization method respectively, predicted the TPH by an accuracy level of 97.83%.

Table 10 Optimal solutions of numerical optimization model
Fig. 8
figure 8

Prediction accuracy of numerical optimzation


This research has studied the remediation of total petroleum hydrocarbon using an environmental friendly method in order to create a clean environment. Factorial design was applied in varying the input parameters (pH, mass of substrate, moisture content and turning) of land farming treatment in order to ascertain the optimal conditions for the procedure. The significant contributions of each input variables which are pH, moisture content, mass of substrate and turning rate associated in the land farming treatment process revealed that, turning rate with 83% was the highest contribution while pH, moisture content and mass of substrate had 4.36, 0.48 and 0.046% contributions, respectively. The numerical optimization done to be sure of the desirability of the absolute model revealed that with initial contamination concentration of 5000 mg kg− 1; turning rate of 5 times weekly, pH of 6.01, moisture content of 10% and substrate mass of 1 kg will achieve a minimum TPH value of 636 mg kg− 1 with 98.6% reliability thus validating the factorial experimental design established for this study.

Availability of data and materials

Not applicable.


  1. IEA. 2012 Key World Energy Statistics. Paris: International Energy Agency; 2012.

    Google Scholar 

  2. Ajiboye AV, Lasisi KH, Babatola JO. Evaluation of the effect of sodium hydroxide solution on biogas yield of anaerobic digestion of poultry waste and the digestate. Int J Energ Water Res. 2018;2:23-31.

    Article  Google Scholar 

  3. Gharaee A, Khosravi-Nikou MR, Anvaripour B. Hydrocarbon contaminated soil remediation: a comparison between Fenton, sono-Fenton, photo-Fenton and sono-photo-Fenton processes. J Ind Eng Chem. 2019;79:181–93.

    Article  Google Scholar 

  4. Zabbey N, Uyi H. Community responses of intertidal soft-bottom macrozoobenthos to oil pollution in a tropical mangrove ecosystem, Niger Delta, Nigeria. Mar Pollut Bull 2014;82:167–74.

    Article  Google Scholar 

  5. Sam K, Coulon F, Prpich G. Management of petroleum hydrocarbon contaminated sites in Nigeria: current challenges and future direction. Land Use Policy. 2017;64:133–44.

    Article  Google Scholar 

  6. Ambituuni A, Amezaga J, Emeseh E. Analysis of safety and environmental regulations for downstream petroleum industry operations in Nigeria: problems and prospects. Environ Dev. 2014;9:43–60.

    Article  Google Scholar 

  7. Konne BR. Inadequate monitoring and enforcement in the Nigerian oil industry: the case of shell and Ogoniland. Cornell Int Law J. 2014;47:181–204.

    Google Scholar 

  8. Ite AE, Ibok UJ, Ite MU, Petters SW. Petroleum exploration and production: past and present environmental issues in the Nigeria’s Niger Delta. Ame J Environ Prot. 2013;1:78–90.

    Article  Google Scholar 

  9. Linden O, Palsson J. Oil contamination in Ogoniland, Niger Delta. Ambio. 2013;42:685–701.

    Google Scholar 

  10. UNEP. Environmental assessment of Ogoniland. Nairobi: United Nations Environment Programme; 2011.

    Google Scholar 

  11. Nwilo PC, Badejo OT. Impacts and management of oil spill pollution along the Nigerian coastal areas. In: Administering marine spaces: international issues. Copenhagen: International Federation of Surveyors; 2006. p. 119–33.

    Google Scholar 

  12. Oluwatuyi OE, Ajibade FO, Ajibade TF, Adelodun B, Olowoselu AS, Adewumi JR, et al. Total concentration, contamination status and distribution of elements in a Nigerian State dumpsites soil. Environ Sustain Indic. 2020;5:100021.

    Article  Google Scholar 

  13. Okonofua ES, Babatola JO, Ojuri OO. Evaluation of total petroleum hydrocarbon content in crude oil exploration field. J Niger Assoc Math Phys. 2018;47:229–34.

    Google Scholar 

  14. Lin QX, Mendelssohn IA. Potential of restoration and phytoremediation with Juncus roemerianus for diesel-contaminated coastal wetlands. Ecol Eng. 2009;35:85–91.

    Article  Google Scholar 

  15. Gomez F, Sartaj M. Field scale ex-situ bioremediation of petroleum contaminated soil under cold climate conditions. Int Biodeter Biodegr. 2013;85:375–82.

    Article  Google Scholar 

  16. Smith E, Thavamani P, Ramadass K, Naidu R, Srivastava P, Megharaj M. Remediation trials for hydrocarbon-contaminated soils in arid environments: evaluation of bioslurry and biopiling techniques. Int Biodeter Biodegr. 2015;101:56–65.

    Article  Google Scholar 

  17. Okonofua ES, Babatola JO, Ojuri OO, Lasisi KH. Determination of suitable TPH remediation approach via MANOVA and inferential statistics assessment. Remediation. 2020;30:75–87.

    Article  Google Scholar 

  18. Wang SY, Kuo YC, Hong A, Chang YM, Kao CM. Bioremediation of diesel and lubricant oil-contaminated soils using enhanced landfarming system. Chemosphere. 2016;164:558–67.

    Article  Google Scholar 

  19. Petavy F, Ruban V, Conil P, Viau JY, Auriol JC. Two treatment methods for stormwater sediments–pilot plant and landfarming–and reuse of the treated sediments in civil engineering. Environ Technol. 2009;30:825–30.

    Article  Google Scholar 

  20. Giasi CI, Morelli A. A landfarming application technique used as environmental remediation for coal oil pollution. J Environ Sci Heal A. 2003;38:1557–68.

    Article  Google Scholar 

  21. Hejazi RF, Husain T, Khan FI. Landfarming operation of oily sludge in and region–human health risk assessment. J Hazard Mater. 2003;99:287–302.

    Article  Google Scholar 

  22. Khan FI, Husain T, Hejazi R. An overview and analysis of site remediation technologies. J Environ Manage. 2004;71:95–122.

    Article  Google Scholar 

  23. Brown LD, Cologgi DL, Gee KF, Ulrich AC. Bioremediation of oil spills on land. In: Fingas M, editor. Oil spill science and technology. 2nd ed. Cambridge: Gulf Professional Publishing; 2017. p. 699–729.

    Chapter  Google Scholar 

  24. Besalatpour A, Hajabbasi MA, Khoshgoftarmanesh AH, Dorostkar V. Land farming process effects on biochemical properties of petroleum-contaminated soils. Soil Sediment Contam. 2011;20:234–48.

    Article  Google Scholar 

  25. Widada J, Nojiri H, Omori T. Recent developments in molecular techniques for identification and monitoring of xenobiotic-degrading bacteria and their catabolic genes in bioremediation. Appl Microbiol Biot. 2002;60:45–59.

    Article  Google Scholar 

  26. Chien HY, Kao CM, Jou CJ, Yang PY, Huang CC. Application of enhanced bioremediation to clean up diesel-oil contaminated soils: laboratory microcosm study. J Biotechnol. 2008;136:678–87.

    Google Scholar 

  27. Brown DM, Okoro S, van Gils J, van Spanning R, Bonte M, Hutchings T, et al. Comparison of landfarming amendments to improve bioremediation of petroleum hydrocarbons in Niger Delta soils. Sci Total Environ. 2017;596-97:284–92.

    Article  Google Scholar 

  28. Maila MP, Cloete TE. Bioremediation of petroleum hydrocarbons through landfarming: are simplicity and cost-effectiveness the only advantages? Rev Environ Sci Biotechnol. 2004;3:349–60.

    Article  Google Scholar 

  29. Heitzer A, Sayler GS. Monitoring the efficacy of bioremediation. Trends Biotechnol. 1993;11:334–43.

    Article  Google Scholar 

  30. Liu SH, Zeng GM, Niu QY, Liu Y, Zhou L, Jiang LH, et al. Bioremediation mechanisms of combined pollution of PAHs and heavy metals by bacteria and fungi: a mini review. Bioresour Technol. 2017;224:25–33.

    Article  Google Scholar 

  31. Adeniyi AA, Afolabi JA. Determination of total petroleum hydrocarbons and heavy metals in soils within the vicinity of facilities handling refined petroleum products in Lagos metropolis. Environ Int. 2002;28:79–82.

    Article  Google Scholar 

  32. Bhattacharya M, Biswas D. Enhancement of waste engine oil biodegradation by optimization of media using factorial design study. Indian J Biotechnol. 2014;13:293–300.

    Google Scholar 

  33. Jasmine J, Mukherji S. Evaluation of bioaugmentation and biostimulation effects on the treatment of refinery oily sludge using 2n full factorial design. Environ Sci-Proc Imp. 2014;16:1889–96.

    Google Scholar 

  34. Agarry SE, Ogunleye OO. Factorial designs application to study enhanced bioremediation of soil artificially contaminated with weathered bonny light crude oil through biostimulation and bioaugmentation strategy. J Environ Prot. 2012;3:748–59.

    Article  Google Scholar 

  35. USEPA. SW-846 Test Method 3540C: Soxhlet Extraction. Washington, DC: US Environmental Protection Agency; 1996.

    Google Scholar 

  36. ASTM. Standard Practice for Extraction of Solid Waste Samples for Chemical Analysis Using Soxhlet Extraction. West Conshohocken: American Society for Testing and Materials; 2005.

    Google Scholar 

  37. Sepahi AA, Golpasha ID, Emami M, Nakhoda AM. Isolation and characterization of crude oil degrading bacillus spp. J Environ Health Sci. 2008;5:149–54.

    Google Scholar 

  38. Montgomery DC. Design and analysis of experiments. 8th ed. Hoboken: Wiley & Sons, Inc; 2008.

    Google Scholar 

  39. Forthofer RN, Lee ES, Hernandez M. Biostatistics: a guide to design, analysis and discovery. 2nd ed. Burlington: Academic Press; 2007.

    MATH  Google Scholar 

  40. Meloun, M, Militky J. Statistical data analysis: a practical guide. New Delhi: Woodhead Publishing; 2011.

    Book  MATH  Google Scholar 

  41. Tavares APA, Coelho MAZ, Agapito MSM, Coutinho JAP, Xavier AMRB. Optimization and modeling of laccase production by Trametes versicolor in a bioreactor using statistical experimental design. Appl Biochem Biotech. 2006;134:233–48.

    Article  Google Scholar 

  42. Singh G, Ahuja N, Sharma P, Capalash N. Response surface methodology for the optimized production of an alkalophilic laccase from γ-proteobacterium JB. Bioresources. 2009;4:544–53.

    Google Scholar 

  43. Cook RD. Detection of influenctial observations in linear regression. Technometrics. 1977;19:15–18.

    MathSciNet  Google Scholar 

Download references


The authors express gratitude to the Federal Government of Nigeria for providing funding for this research which was undertaken at the Federal University of Technology Akure through the Tertiary Education Fund (TETfund) with Grant No REG/SSA/P.13735/75.


This research was funded by the Tertiary Education Fund (an organization under the auspices of the Federal Government of Nigeria) with Grant No REG/SSA/P.13735/75.

Author information

Authors and Affiliations



OES conceived and investigated the experiment, and was also a major contributor in writing the manuscript. LKH analyzed data and interpreted the result, and was also a major contributor in writing of the manuscript. AE supervised the research and also edit the manuscript. All authors read and approved the final manuscript.”

Corresponding author

Correspondence to Kayode H. Lasisi.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Interphase of numerical optimization model for minimizing total petroleum hydrocarbon. Figure S2. Ramp solution of numerical optimization.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Okonofua, E.S., Lasisi, K.H. & Atikpo, E. Factorial design study of total petroleum contaminated soil treatment using land farming technique. Sustain Environ Res 31, 5 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: