InTechOpen uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Chemistry » "Statistical Approaches With Emphasis on Design of Experiments Applied to Chemical Processes", book edited by Valter Silva, ISBN 978-953-51-3878-5, Print ISBN 978-953-51-3877-8, Published: March 7, 2018 under CC BY 3.0 license. © The Author(s).

Chapter 6

Evaluation of Factors Affecting Chemical Extraction of Co Ions from Contaminated Soil

By Ivana Smičiklas
DOI: 10.5772/68066

Article top


System responses obtained for different extracting agents under experimental conditions defined in Table 2: (a) extraction efficiency of Co ions from the contaminated soil and (b) final pH of extracts.
Figure 1. System responses obtained for different extracting agents under experimental conditions defined in Table 2: (a) extraction efficiency of Co ions from the contaminated soil and (b) final pH of extracts.
Contour plots for cobalt extraction efficiency using: (a) EDTA, (b) CA, and (c) HCl solutions, presenting the relationships between calculated responses and statistically significant factors, at the constant intermediate value of the third factor (A—reagent concentration, B—soil/solution ratio, C—extraction time).
Figure 2. Contour plots for cobalt extraction efficiency using: (a) EDTA, (b) CA, and (c) HCl solutions, presenting the relationships between calculated responses and statistically significant factors, at the constant intermediate value of the third factor (A—reagent concentration, B—soil/solution ratio, C—extraction time).
Contour plots for the final pH of extracts: (a) EDTA, (b) CA, and (c) HCl, presenting the relationships between calculated responses and statistically significant factors, at the constant intermediate value of the third factor (A—reagent concentration, B–soil/solution ratio, C—extraction time).
Figure 3. Contour plots for the final pH of extracts: (a) EDTA, (b) CA, and (c) HCl, presenting the relationships between calculated responses and statistically significant factors, at the constant intermediate value of the third factor (A—reagent concentration, B–soil/solution ratio, C—extraction time).
Overlaid contour plots for determination of process conditions which lead to achievement of target system responses using (a) Na2EDTA, (b) CA, and (c) HCl (A—reagent concentration, B—soil/solution ratio, C—extraction time).
Figure 4. Overlaid contour plots for determination of process conditions which lead to achievement of target system responses using (a) Na2EDTA, (b) CA, and (c) HCl (A—reagent concentration, B—soil/solution ratio, C—extraction time).

Evaluation of Factors Affecting Chemical Extraction of Co Ions from Contaminated Soil

Ivana Smičiklas
Show details


Excessive concentrations of cobalt (Co) ions in the soil cause quality degradation and pose a significant hazard to biota. One of the options for the permanent separation of the pollutant from soil matrix is extraction by chemical reagents. In this study, response‐surface methodology (RSM) was applied to evaluate the factors affecting Co extraction from contaminated calcareous soil. Solutions of disodium ethylenediaminetetraacetate (Na2EDTA), citric acid (CA), and HCl were considered as leaching media. Reagent concentration, soil to solution ratio, and extraction time were selected as process variable, while Co extraction efficiencies and final pH values of extracts were the measured responses. The effect of factor variation between three levels was studied using Box‐Behnken experimental design. By statistical analysis, the most influential factors were determined for each reagent, and the model equations were proposed for the prediction of system responses. Overlaid contour plots were used for the analysis of the effect of process conditions on both responses simultaneously. Given that each case of contamination is unique and requires extensive research before the remediation is implemented full‐scale, it was shown that experimental design methodology is a smart approach for the assessment and comparison between the treatments.

Keywords: cobalt ions, soil contamination, chemical extraction, Box‐Behnken design, ANOVA

1. Introduction

Cobalt (Co) is one of the essential heavy metals, which is vital at trace levels for proper functioning of the human metabolism [1]. However, at concentrations higher than optimal, essential metals become toxic. The poisoning by Co is commonly a result of drinking water contamination, high ambient air and soil concentrations, and consequent entrance and bioaccumulation of Co in the food chain. Therefore, development and application of remediation methods are necessary in order to mitigate the negative impact of Co on the ecosystem.

The average concentration of Co in the Earth’s crust is 10 mg/kg [2]. The levels of Co in the soil are influenced by the pedogenic factors, as well as by the characteristics of a parent rock material. For different climatic zones, the values of 0.05 and 300 mg/kg were reported as the minimum and maximum naturally occurring concentrations in the soil [3]. Lower amounts were found in the soils from the northern regions, which originated from the glacial deposits. The average concentrations of Co in the northern parts of Ukraine, Russia [2], and Sweden [4] are 3.5, 5.5, and 7.1 mg/kg, respectively. In contrast, naturally higher values can be found in the areas with arid or semiarid climate, like in Egypt [3].

The property which affects the manifestation of both beneficial and toxic effects is the solubility of Co in a particular environment. In respect to total concentration of Co in the soil, the fraction available to plants is more important. The accessibility to living organisms is the highest for free Co ions and water soluble complexes, while its metallic form and insoluble compounds commonly exhibit very low bioaccessibility. The available content of Co in different soil types, determined using 2.5% acetic acid as an extracting agent, was found to vary from 0.05 to 1 ppm [3]. On one hand, the lack of available Co leads to Co‐deficiency in living organisms, while on the other, Co and its compounds are highly toxic in excessive amounts causing serious cell and tissue damages (LD50 values for intake of Co salts for rats are in the range 150–500 mg/kg of body weight [5]).

The main sources of environmental pollution with Co are industrial activities. Utilization of Co as catalyst in the chemical industry, in the production of dyes and pigments, magnetic recording media, alloys, batteries, etc., makes this element of the strategic importance for military, industrial, and commercial applications [6]. In numerous studies, Co concentrations in soil and sediments were measured to define its levels and identify the main pollution sources. Soil pollution by Co in Shenzhen soil (China) was attributed to uncontrolled discharge of industrial wastewater from factories that produce or use chemical compounds or alloys containing Co [7]. Furthermore, ceramic industry was highlighted as the source of the sediment contamination with Co and other heavy metals in Jiangsu Province (China) [8]. The highest concentrations were observed in the surface soil layer, within the 20–22 cm depth. In north greater Cairo (Egypt), direct discharge of industrial wastewaters to irrigation water canals during 30 years provoked significant contamination of soil with Co (146 mg/kg) [9]. The activity of a smelter in Zambia was found to result in dust with elevated concentrations of Co and other heavy metals, causing the contamination of soil and plants [10]. The highest measured Co concentration was 606 mg/kg found in the 0–2 cm of soil depth [11]. Furthermore, Co has been found in increased concentrations in at least 426, out of 1636 most serious hazardous waste sites in the USA, identified by Environmental Protection Agency (US EPA) [12]. The contamination of soil with Co is therefore a global problem with a tendency of steady increase.

Contaminated areas need to be treated using the techniques based either on the pollutant separation from the soil matrix or on its stabilization [13]. The levels and properties of heavy metal pollutants, as well as the physicochemical properties of the soil, govern the distribution of heavy metals which is an important factor for the selection of effective remediation method. When soil is exposed to Co‐containing solution, several processes can occur at the surface of soil constituents: electrostatic attraction, ion‐exchange, complexation, co‐precipitation, precipitation, occlusion, diffusion, and migration. Various single‐step leaching procedures and sequential extractions are in use for the determination of metal distribution and mobility in the soil. Sequential extraction of the soil by so‐called Tessier protocol is frequently applied [14]. This procedure is based on the consecutive application of selective chemical reagents in five extraction steps (F1–F5), with the aim to extract the pollutant associated with different soil fractions: ion‐exchangeable (F1), acid soluble (F2), bound to Fe, Mn‐oxides (F3), complexed by organic matter (F4), and incorporated in the residual fraction (F5). As the strength of the bonds between the metal and the soil constituents increases along the extraction scheme, the metal mobility decreases in the same manner.

Clean‐up of the contaminated soil can be completed by the method of chemical extraction if the pollutant is distributed between relatively mobile fractions F1–F3 [15]. Extraction processes require mixing of contaminated soil with leaching solutions (solutions of acids, inorganic salts, chelating agents, surfactants, etc.) that cause a transfer of pollutant from the soil matrix into the liquid phase. If the chemical extraction is chosen among other alternative methods, the optimization of the method performance requires extensive research on the effects of a large number of variables. Some of the previous studies have addressed the impact of reagent type and concentration, reaction time, soil/solution ratio, pH, temperature, etc., on the efficiency of heavy metals extraction from the soil matrix [16, 17, 18, 19]; however, the experiments were conducted by varying one factor at the time and little attention has been paid onto evaluation and comparison of factor effects, their interactions and process optimization.

Taking into account numerous potentially important factors, experimental design methodology (DOE) can be a useful strategy for the analysis of different soil treatments. In that sense, the present chapter aimed to explore the applicability of DOE approach for the analysis of the factors affecting chemical extraction of Co from contaminated soil and for the prediction of system responses.

2. Design of experiment (DOE)

The analysis and optimization of virtually any process can be conducted by experimental design methodology (DOE). Compared to the classical approach which implies variation of one factor at the time, DOE relies on simultaneous variation of all factors and reveals the most influential ones, the significant interactions between the factors, and the optimal levels of the factors. This methodology aims to describe or explain the variation of information under conditions that are hypothesized to reflect the variation. Therefore, DOE is a planned approach to find out the relationships between process variables and process responses, in relatively small number of experimental trials. Depending on the process under consideration, as well as on the number and the type of information of interest, different types of experimental design are in use, such as full and fractional factorial designs, response‐surface designs, mixture designs, random block, Latine squares designs, etc. [20].

The response‐surface methodology (RSM) was developed in order to optimize chemical reactions used on an industrial scale [21]. The application of RSM is of relevance in situations where several input variables influence the performance or quality characteristics of the process, which are called the responses. The input variables, usually denoted as independent variables, are controlled. Variation in experimental conditions provokes changes in the process response. Dependence between process responses and two process variables is a smooth surface called the response surface. The purpose of the response surface is to locate the conditions which lead to the achievement of the minimum/maximum response. Furthermore, in cases when several responses are of interest, RSM can give the optimal conditions for their target values [22]. The response values (y) are related to the process variables (x1, x2…) according to the model:

y= f(x1, x2,)+ε

In Eq. (1), f(x1, x2…) is called the response surface and represents the mean response at each x, while ε represents other sources of variability (measurement error, background noise, the effect of other variables, etc.), and it is usually treated as a statistical error.

The function f(x1, x2…) can be given as the first‐order [Eq. (2)], or the second‐order model [Eq. (3)]. The first‐order model describes experimental conditions which provide the responses without peaks, and it is applicable for the description of local area responses without function extremes.

f(xi)= β0+βixi+ε

In Eq. (2), βi denotes the linear effect of the process variable xi. Given that it includes only the main effects of the two or more variables, the form of the first‐order model given by Eq. (2) is called the main effects model. In some processes, the interactions between studied variables may be significant; thus, the interaction term must be added to the first‐order model, which introduce curvature into the response function.

The advantage of RSM is that the lack of fit can be estimated and the adequacy of the model can be estimated. The lack of fit test provides the answer whether the proposed model adequately describes system response (H0—null hypothesis) or not (HA—alternative hypotheses). If the proposed model is suitable, there is no lack of fit. When the first‐order polynomial regression model is not adequate, the applicability of the second‐order model is analyzed:

f(xi)= β0+βixi+βijxixj+βiixi2+ε

The term βi,j represents the interaction effect between two factors xi and xj, and βi,i is the quadratic effect of factor xi. The meanings of the other terms in Eq. (3) are the same as aforementioned. The designs that correspond to the surface response methodology are Box‐Behnken design, central composite design, and optimal designs.

The results of a literature survey showed that different types of experimental design were used for the investigation of Co sorption by potential sorbent materials such as apatite, zeolite, bauxite residue, and dried activated sludge [23, 24, 25]; however, DOE approach has not been applied, to this point, for the study of Co sorption or leaching from the soil matrix. In fact, just a few examples of RSM application for the analysis of a soil remediation process can be found [26, 27]. In this study, a three‐level‐three‐factor Box‐Behnken experimental design was selected for the evaluation of the effect of factors which influence the Co extraction from contaminated soil. Box‐Behnken design is a three‐level factorial design which considers factors at their low (–1) and high (+1) values, as well as at their arithmetic mean (0) [28]. This is the type of quadratic design which does not contain the factorial or fractional factorial designs but considers the midpoints of the edges of the process space and its center. In general, it is an alternative to full factorial design at three levels.

3. Materials and methods

3.1. Soil contamination and clean‐up experiments

The soil was sampled from the site of Vinča Institute of Nuclear Sciences (Belgrade, Serbia), from the surface layer (0–20 cm). Prior to contamination with Co ions, soil was dried at room temperature, homogenized, grained, and sieved in order to separate the fraction of particles with the diameter <2 mm. The soil at this locality was characterized as weak alkaline (pHH2O = 8, pHKCl = 7), with CaCO3 content of 5.4%, total organic content (TOC) of 2.1%, and cation exchange capacity (CEC) of 13 meq/100 g [29]. Mineralogical analysis revealed quartz, kyanite, and muscovite as main crystal phases, whereas mica, albite, kaolinite, and calcium, magnesium–carbonate were present in lower quantities [29].

In the first step, the soil was artificially contaminated by Co ions. The solution containing 0.0012 mol/L of Co(NO3)2(Merck, p.a.) was mixed with 100 g of dried soil at the solid/liquid ratio 1:20. After 24 h of equilibration, the liquid phase was separated from the soil by filtration and the contaminated sample was dried in the air atmosphere. The concentration of Co ions sorbed by the soil was calculated as the difference between its initial and the residual concentration.

Soil clean‐up was conducted using the method of chemical extraction. Based on the high extraction efficiencies observed in the previous study [18], the solutions of disodium ethylenediaminetetraacetate (Na2EDTA), citric acid (CA), and HCl were selected as extracting agents. The contaminated soil was mixed with the reagent solutions in the centrifuge tubes which were placed on the rotary shaker and agitated at the constant speed (10 rpm) at ambient temperature (20 ± 2°C).

The effects of the variation of reagent concentration, solid/liquid ratio, and contact time, among three levels, were investigated. Experimental variables and their levels are given in Table 1. The matrix with the experimental conditions for simultaneous variation of process parameters (Table 2) was generated using Minitab Release Software 13.1.

Independent factorSymbolLevel 1 coded value (−1)Level 2 coded value (0)Level 3 coded value (−1)
Reagent concentration (mol/L)A0.00050.050250.1
Soil/solution ratioB51525
Contact time (h)C13.56

Table 1.

Experimental factors and their levels used in the Box‐Behnken design.

Experimental runIndependent factor

Table 2.

Combinations of factors and their levels according to Box‐Behnken design matrix.

Extraction of Co ions was performed using each of the selected reagents, according to the conditions given in Table 2. After specified contact times, the suspensions were centrifuged for 10 min at 10,000 rpm and the extracted concentrations of Co ions were measured in clear supernatants. The chemical extraction experiments were conducted in duplicate.

In the sorption and extraction experiments, Co concentrations were determined by Perkin Elmer 3100 atomic absorption spectrometer (AAS), using a wavelength of 252.1 nm. Standards for instrument calibration were prepared by diluting certified Perkin Elmer standard (1000 mg/L), and the calibration was repeated after every 10‐sample measurement. The detection limit was 0.05 mg/L, whereas the deviations among five replicate measurements for each sample were lower than 3%.

The percentages of the extracted Co ions were used as the main response function of the process. In addition, pH values of the leaching solutions after the completion of the process (denoted as final pH) were measured using WTW InoLab pH‐meter, and also considered as system response.

3.2. Statistical analysis

Mean values of measured parameters, obtained from duplicate extraction experiments, were used as system responses for data interpretation and statistical analysis. Experimental point designation, analysis of variance (ANOVA), fitting of regression polynomial models, and graphical presentations (ternary plots) were performed using the statistical software MINITAB Release 13.2. The statistical analysis was done at the confidence level α = 95%.

4. Results and discussion

4.1. System responses at different experimental conditions

After the contamination step, the total amount of Co ions sorbed by the investigated soil was found to be 1390 mg/kg. High sorption affinities of soils and soil components toward Co are well documented [24, 30, 31, 32]. The sorbed amount is in agreement with the results of a previous study in which the soil from the same location was contaminated with Co under the range of experimental conditions [18].

The variation of chemical reagents and other experimental conditions according to the Box‐Behnken design resulted in variation of system responses (Figure 1).


Figure 1.

System responses obtained for different extracting agents under experimental conditions defined in Table 2: (a) extraction efficiency of Co ions from the contaminated soil and (b) final pH of extracts.

The mean values of Co extraction efficiency were found to fluctuate in the range 10–72% for Na2EDTA, 1–66% for CA, and 0–71% for HCl as extracting agent (Figure 1a). The results signify that, at certain levels of the considered factors, amounts of extracted Co can reach high values using tested reagents. It is evident that the change of experimental conditions differently affects the performance of tested reagents and that effect of HCl was the lowest when compared to other reagents under the same sets of experimental conditions.

The final pH values were also largely dependent on the type of reagent and process variables (Figure 1b). The initial pH values of extracting solutions at their lower, middle, and higher selected concentrations were: 4.4, 4.6, and 4.9 for Na2EDTA, 2.0, 2.2, and 3.5 for CA, and 1.0, 1.3, and 3.3 for HCl. After the reaction of contaminated soil with Na2EDTA solutions, pH values were in the range 4.4–8.8, while using solutions of CA and HCl, pH values varied between 2.8–7.8 and 1.6–7.6, respectively. The observed increase in pH, after interaction with the soil, can be attributed to the buffering capacity of the soil which mainly originates from its carbonate content [18].

The knowledge of the nature and the strength of the bonds established between added Co ions and soil constituents are important for interpretation and understanding the effect of different chemical reagents. By means of the sequential extraction analysis, the portions of Co bonded in different fractions of the investigated soil were previously analyzed [18]. Considering various contamination levels and aging times (1 h to 30 days), Co ions were found to be associated with Fe, Mn‐oxide fraction (F3), carbonate/acid soluble (F2), and ion‐exchangeable (F1) fraction (F1 + F2 + F3 > 92%). For the comparison, the majority of naturally occurring Co in the soil from the Vinča locality was found in F3 and residual (F5) phase [29], indicating that sorbed cobalt ions have significantly higher mobility. These results signify that chemical extraction can be considered for soil clean‐up [33].

Fe‐ and Mn‐oxides and hydroxides have been generally recognized as the major substrates for Co ions in the soil [34]. High affinities of such minerals toward Co ions were also revealed in the studies of Co sorption by Fe‐ and Mn‐oxides [35, 36, 37], as well as by geochemical modeling of Co speciation in different soil samples [38]. The amounts of Co found in F2 phase of the sequential extraction correspond to the fraction associated with carbonates, as well as the fraction chemisorbed on the surface of soil components. The crystal lattice of calcium‐carbonate can accommodate high quantities of Co and other divalent metals, implying that calcite is an important sorbent for metals in calcareous soils [39]. In addition, carbonate phase is of a crucial importance for high soil pH, which results in higher stability of sorbed metal ions [40]. Finally, the ion‐exchange between Co and other divalent cations contained in the soil minerals also contributes to Co retention in the soil [41].

Taking into account the effects achieved by investigated reagents (Figure 1a), it can be concluded that both complexing agents and the mineral acid can release Co ions sorbed within different soil fractions. The extraction potential of HCl solutions is governed by the changes of pH. As the pH becomes more acidic, the dissolution of soluble soil constituents occurs, as well as the exchange reactions between H+ and metal ions, and the protonation of active surface sites which become more positively charged [31, 42, 43]. Depending on the pH, these reactions cause the liberation of metal ions bonded by ion‐exchange mechanism, specifically sorbed, and associated with carbonates and oxide minerals. Solubility of metal ions in the soil is also enhanced in the presence of chelating anions. These organic ions form multiple bonds with the metal cation, essentially in the form of a ring, which exhibits high stability in aqueous media. The chelating agents have a potential to extract the metal ions bonded to all non‐residual fractions [44].

4.2. Evaluation of the effects of factors by Box‐Behnken design

By using actual values of process variables, the response variable (y) was fitted to a second‐order polynomial model as follows:


where y is system response, β0 is the intercept, βi is the linear effect, and βii is the quadratic effect of the process variable xi, while βi,j is the effect of the interaction between two independent process variables xi and xj.

Applicability of different models was tested by ANOVA. In order to determine which mathematical model adequately fits the obtained experimental results, full quadratic model was initially applied [Eq. (4)]. Statistical calculations were assessed on the basis of the F‐ and p‐values. F‐values are obtained from Fisher’s test and represent the ratio of the mean square due to regression and mean square due to residual error, whereas p‐value is defined as the smallest level of significance which leads to the rejection of null hypothesis. Therefore, higher F and lower p‐values (p < 0.05) represent statistically significant terms. By omitting the insignificant terms (p > 0.05) and considering the values of determination coefficients (R2) together with the results of lack of fit, the most appropriate models were found.

For Na2EDTA as Co leaching reagent, the most suitable model included only linear terms, following the equation:


According to the proposed Eq. (5), Co desorption efficiency was dependent on the linear effect of all selected factors. However, by observing the values of coefficients next to the respective factors, it can be concluded that the highest effect on the process response was achieved by the change of initial reagent concentration, followed by the effect of applied reagent volume, and finally the effect of contact time. The positive signs of the calculated coefficients indicate that leached amounts of Co increase with the increased levels of factors.

The coefficient of determination (R2) was calculated according to the equation:


where SSres denotes the sum of squares of the residuals, while SStot is the total sum of squares. SStot was calculated as the sum of squared differences between the different system responses values (yi) and its mean. Furthermore, the SSres is the sum of the squared differences between each system variable and its value predicted by the model. Knowing that the R2 values increase with the number of terms in the model even if their effect is insignificant, an incorrect model could be proposed if R2 magnitude is solely considered. Therefore, the value of the adjusted determination coefficient (R2adj) was included. In contrast to their effect on R2 values, insignificant terms in the model provoke the decrease in R2adj values.

The summary of ANOVA, F and p‐values obtained for significant model terms and the results of the lack of fit are given in Table 3. Using the experimental results obtained for Na2EDTA extraction efficiency, calculated R2 and R2adj values were 99.02 and 98.76%, respectively. Practically, the predicted model adequately describes even 99.02% of experimental results. Giving that F‐value was <1 and p‐value was >0.05, the proposed model can be regarded as adequate.

Extraction efficiency
Lack of fit4.480.1965.370.1603.130.310
Final pH of extracts
Lack of fit2.510.1516.370.065.140.18

Table 3.

Analysis of variance.

The following model was developed for CA as a leaching reagent:

Co(CA)=54.527+27.779 A+5.304 B+3.990 C23.238A2

The extraction of Co ions with CA solutions was found to be dependent on all investigated factors with the statistical significance. Changes in the level of reagent initial concentration have caused the highest impact on the leaching efficiency. Furthermore, the squared term in the derived equation (A2) indicates the existence of curvature in the response surface. The obtained Eq. (7) fitted 97.88% of the experimental results, with the calculated R2adj value of 97.03%.

The extraction of Co using HCl can be described by the model which includes the following linear and interaction terms:

Co(HCl)=28.31+24.15 A+19.20 B+12.39AB

For the proposed model, determined R2 was 94.69%, while R2adj was 93.25%. The analysis of factors affecting the efficiency of applied inorganic acid indicates that the changes of contact time did not significantly influence on the amount of leached cobalt. On the other hand, changes in the concentration and volume of the reagent were significant, as well as the interaction of these two parameters.

Furthermore, the models for the description of final pH values were developed. Considering the data obtained using Na2EDTA, the following equation was found suitable, with calculated values of R2 = 96.48% and R2adj = 95.52%:


The final pH values of the CA extracts can be described (R2 = 96.33%, R2adj = 95.33%) by Eq. (10):


Eqs. (9) and (10) indicate that variation of initial concentration and volume of both investigated complexing agents had an effect on final pH. The negative signs of the corresponding coefficients denote that increased concentration levels and increased levels of applied volume influenced the decrease in pH.

The complex mathematical model was obtained for describing the solution pH values after soil leaching with HCl:


The adequacy of Eq. (11) was confirmed by high R2 and R2adj values (98.32% and 97.06%, respectively). The model indicates that the variation of all studied factors, as well as AB and BC interactions, significantly affected the pH of resulting extracts. The HCl concentration and volume exhibited negative correlation with pH, while the contact time and pH were positively correlated.

Inside the range of experimental conditions covered by the design, the efficiency of Co extraction can be predicted using the mathematical models proposed by Eqs. (5), (7), and (8). The graphical interpretation of these equations, that help visualize the shape of the response surface, is given in Figure 2 in form of contour plots. The graphs express the relationships between calculated responses and statistically significant variables, at the constant intermediate value of the third parameter. The lines on the graphs represent constant values of the system response.


Figure 2.

Contour plots for cobalt extraction efficiency using: (a) EDTA, (b) CA, and (c) HCl solutions, presenting the relationships between calculated responses and statistically significant factors, at the constant intermediate value of the third factor (A—reagent concentration, B—soil/solution ratio, C—extraction time).

The parallel lines on the contour plots in Figure 2a refer to the linear effect of process variables on Co desorption efficiency by Na2EDTA, with the highest yield of extraction achieved at their maximum levels. The curvature in system responses occurred in case when CA was used as reagent, due to the quadratic term in the model equation, whereas in case of HCl due to the significant interaction effect between the reagent concentration and applied volume.

Similar predictions were made for final pH values as response functions, using Eqs. (9)(11), and the contour plots are given for statistically significant parameters (Figure 3).


Figure 3.

Contour plots for the final pH of extracts: (a) EDTA, (b) CA, and (c) HCl, presenting the relationships between calculated responses and statistically significant factors, at the constant intermediate value of the third factor (A—reagent concentration, B–soil/solution ratio, C—extraction time).

The contour graphs, as well as the ANOVA results, signify that increase of the initial reagent concentration provoked decrease in the final pH. Furthermore, pH declined with increased volume of the reagents. In contrast, longer reaction times contribute only to higher pH of HCl solution.

If multiple system responses are considered in the same experiment, their simultaneous analysis could be essential to optimize the process. Construction of overlaid contour plots enables determination of the ranges of process variable which lead to the achievement of the target effect. These plots are constructed by overlaying the contour plots of each considered response. The responses defined in the presented study, final pH values of the filtrate and extraction efficiency, are related: the lower the pH―the higher extraction yield (Figures 2 and 3). However, high acidity of extracting solutions may cause the degradation of soil vital properties due to significant degradation of both mineral and organic phase [43], making the soil inappropriate for on‐site disposal and re‐vegetation after completion of the treatment. Furthermore, the obtained filtrates represent liquid waste which requires further management (neutralization, metal recovery, reagent recovery, etc.). Thus, both the wastewater and the processed soil would need to be neutralized, increasing the complexity and cost of the treatment.

Therefore, the performance of different reagents was compared using near neutral final pH (5–7) and Co extraction efficiency >50%, as target values. Under these conditions, overlaid contour plots were constructed for Na2EDTA and CA (Figure 4a and b). On the other hand, desired response values could not be reached simultaneously using HCL, and the target value of Co extraction was decreased to >40% (Figure 4c).


Figure 4.

Overlaid contour plots for determination of process conditions which lead to achievement of target system responses using (a) Na2EDTA, (b) CA, and (c) HCl (A—reagent concentration, B—soil/solution ratio, C—extraction time).

The white areas in the constructed plots present the ranges of factors A and B that provide given responses, at different contact times (C). Using Na2EDTA, target response values can be obtained for contact times of 1 h, as well as 3.5 h, by selecting reagent concentration and volume in the designated ranges. Designated area is wider for 3.5 h, which practically means that if the duration of extraction is increased the desired responses can be gained at lower reagent concentrations. For the same extraction times, using CA instead of Na2EDTA (Figure 4b), target values of Co desorption efficiency and solution pH can be obtained using initial concentrations and volumes of the reagent in much lower ranges. The effect of HCl was limited since it is governed by pH decrease (Figure 4c). Approximately 40% of sorbed Co ions can be leached under conditions that simultaneously assure target pH of the extract, at higher contact times (6 and 3.5 h).

Such results may be indicative for further evaluation of reagents. Consideration of the advantages and disadvantages related to the application of a chemical reagent includes the environmental impact and overall costs. Although HCl is often used for chemical leaching at full‐scale, its effect is limited in calcareous soils. Consumption of high amounts of acid reagent would be necessary to provide efficient Co separation, which would in turn results in dissolution and degradation of soil matrix, creation of acidic waste water and soil residue. On the other hand, by selection of factor levels, complexing agents provide more efficient decontamination, even at near neutral pH conditions. Leaching of soil by chelating agents is potentially detrimental to its quality if a part of the chelating agent remains in the soil. As a natural compound, CA undergoes biodegradation process more easily compared to Na2EDTA [45]. Furthermore, CA is less expensive chelating reagent [46], which supports the selection of CA among the tested reagents.

5. Conclusion

The Box‐Behnken design was applied for screening and analyzing the factors affecting soil chemical extraction process using the sample of calcareous soil artificially contaminated with Co ions. Using Na2EDTA, CA, and HCl as reagents, the effect of variation of three factors (concentration of the chemical reagent, soil/liquid ratio, and the contact time between the phases) between three levels was analyzed in respect to extraction efficiency and final pH values. The adequacy of different mathematical models, with inclusion of linear or quadratic terms, was tested for the description of experimental results. The extraction efficiency was highly dependent of the applied chemical reagent. Analysis of variance of the chosen responses revealed that Co separation was predominantly affected by the variation of the reagents concentrations. Effect of applied reagents volume had smaller statistical significance, while the contact time played an important role in the performance of complexing agents. Mathematical models were developed to describe the effect of each independent parameter and their interactions on system response. Predicted values of Co recovery, obtained using model equations, were in good agreement with the experimental data. Treatments by different reagents were compared using overlaid contour plots, taking into consideration extraction efficiencies in near neutral range of final pH values. The results demonstrate limited effect of HCl in cancerous soil, while chelating agents can exhibit high efficiencies by selection of factor levels. CA stands out as the most suitable agent, due to its performance, price, and biodegradability. The applicability of RSM for fast assessment and comparison of soil treatments was verified, highlighting the significance of DOE in practical applications.


This work was supported by the Ministry of Education, Science and Technological Development of the Republic of Serbia (Project III43009).


1 - Yamada K. Cobalt: Its role in health and disease. In: Sigel A, Sigel H, Sigel RKO, editors. Interrelations Between Essential Metal Ions and Human Diseases (Metal Ions in Life Sciences). Vol. 13. Springer, Netherlands; 2013, pp. 295-320. DOI: 10.1007/987‐94‐007‐7500‐8_9
2 - Kabata‐Pendias A. Trace Elements in Soils and Plants. 4th ed. Boca Raton: CRC Press; 2010. p. 548. DOI: 10.1017/S0014479711000743
3 - Aubert H, Pinta M, editors. Development in Soil Science (Trace Elements in Soil). Vol. 7, Amsterdam: Elsevier; 1977. p. 394
4 - Eriksson J. Concentrations of 61 trace elements in sewage sludge, farmyard manure, mineral fertiliser, precipitation and in oil and crops. Report. Swedish Environmental Protection Agency, Stockholm, Sweden; 2001. Available from:‐6246‐8.pdf. [Accessed: December 20, 2016]
5 - Speijers GJA, Krajnc EI, Berkvens JM, Van Logten MJ. Acute oral toxicity of inorganic cobalt compounds in rats. Food and Chemical Toxicology. 1982;20:311-314. DOI:10.1016/S0278‐6915(82)80298‐6
6 - Donaldson JD, Beyersmann D. Cobalt and cobalt compounds. In: Ullmann’s Encyclopedia of Industrial Chemistry, Weinheim: Wiley‐VCH Verlag GmbH & Co. KGaA; 2005. DOI: 10.1002/14356007.a07_281.pub2
7 - Xu S, Lin C, Qiu P, Song Y, Yang W, Xu G, Feng X, Yang Q, Yang X, Niu A. Tungsten and cobalt‐dominated heavy metal contamination of mangrove sediments in Shenzhen, China. Marine Pollution Bulletin. 2015;100:562-566. DOI: 10.1016/j.marpolbul.2015.08.031
8 - Liao QL, Liu C, Wu HY, Jin Y, Hua M, Zhu BW, Chen K, Huang L. Association of soil cadmium contamination with ceramic industry: A case study in a Chinese town. Science of the Total Environment. 2015;514:26-32. DOI:10.1016/j.scitotenv.2015.01.084
9 - Lotfy SM, Mostafa AZ. Phytoremediation of contaminated soil with cobalt and chromium. Journal of Geochemical Exploration. 2014;144:367-373. DOI: 10.1016/j.gexplo.2013.07.003
10 - Kříbek B, Majer V, Knésl I, Nyambe I, Mihaljevič M, Ettler V, Sracek O. Concentrations of arsenic, copper, cobalt, lead and zinc in cassava (Manihot esculenta Crantz) growing on uncontaminated and contaminated soils of the Zambian Copperbelt. Journal of African Earth Sciences. 2014;99:713-723. DOI: 10.1016/j.jafrearsci.2014.02.009
11 - Ettler V, Mihaljevič M, Kříbek B, Majer V, Šebek O. Tracing the spatial distribution and mobility of metal/metalloid contaminants in Oxisols in the vicinity of the Nkana copper smelter, Copperbelt province, Zambia. Geoderma. 2011;164:73-84. DOI: 10.1016/j.geoderma.2011.05.014
12 - Agency for Toxic Substances and Disease Registry (ATSDR), “Toxicological Profile for Cobalt,” Atlanta, GA: US; 2004, Available from: [Accessed: 2016-12-15]
13 - Yao Z, Li J, Xie H, Yu C. Review on remediation technologies of soil contaminated by heavy metals. Procedia Environmental Sciences. 2012;16:722-729. DOI: 10.1016/j.proenv.2012.10.099
14 - Tessier A, Campbell PGC, Bisson M. Sequential extraction procedure for the speciation of particulate trace metals. Analytical Chemistry. 1979;51:844-851. DOI: 10.1021/ac50043a017
15 - Mulligan CN, Yong RN, Gibbs BF. Remediation technologies for metal‐contaminated soils and groundwater: An evaluation. Engineering Geology. 2001;60:193-207. DOI: 10.1016/S0013‐7952(00)00101‐0
16 - Zou Z, Qiu R, Zhang W, Dong H, Zhao Z, Zhang T, Wei X, Cai X. The study of operating variables in soil washing with EDTA. Environmental Pollution. 2009;157:229-236. DOI: 10.1016/j.envpol.2008.07.009
17 - Yoo JC, Lee C, Lee JS, Baek K. Simultaneous application of chemical oxidation and extraction processes is effective at remediating soil Co‐contaminated with petroleum and heavy metals. Journal of Environmental Management. 2017;186:314-319. DOI:10.1016/j.jenvman.2016.03.016
18 - Smičiklas I, Dimović S, Jović M, Milenković A, Šljivić‐Ivanović M. Evaluation study of cobalt(II) and strontium(II) sorption–desorption behavior for selection of soil remediation technology. International Journal of Environmental Science and Technology. 2015;12:3853-3862. DOI: 10.1007/s13762‐015‐0817‐y
19 - Naghipour D, Gharibi H, Taghavi K, Jaafari J. Influence of EDTA and NTA on heavy metal extraction from sandy‐loam contaminated soils. Journal of Environmental Chemical Engineering. 2016;4:3512-3518. DOI: 10.1016/j.jece.2016.07.034
20 - Lazić Ž, Design of Experiments in Chemical Engineering. Weinheim: Wiley‐VCH Verlag GmbH&Co.; 2004. p. 610. ISBN 3‐527‐31142‐4
21 - Box GEP, Wilson KB. On the experimental attainment of optimum conditions. Journal of the Royal Statistical Society. 1951;13:1-45. DOI: 10.1007/978‐1‐4612‐4380‐9_23
22 - Dean AM, Voss D. Design and analysis of experiments. Springer Texts in Statistics. Springer, New York; 1999; 768 p, DOI:10.1007/b97673
23 - Smičiklas I, Sljivić‐Ivanović M. Evaluation of factors influencing Co2+ removal by calcinated bone sorbent using experimental design methodology. Journal of Environmental Science and Health, Part A. 2012;47:896-908. DOI: 10.1080/10934529.2012.665006
24 - Šljivić‐Ivanović M, Smičiklas I, Dimović S, Jović M, Dojčinović B. Study of simultaneous radionuclide sorption by mixture design methodology. Industrial & Engineering Chemistry Research. 2015;54:11212-11221. DOI: 10.1021/acs.iecr.5b03448
25 - Frišták V, Remenárová L, Lesný J. Response surface methodology as optimization tool in study of competitive effect of Ca2+ and Mg2+ ions in sorption process of Co2+ by dried activated sludge. Journal of Microbiology, Biotechnology and Food Sciences. 2012;1:1235-1249
26 - MartínezÁlvarez LM, Lo Balbo A, Mac Cormack WP, Ruberto LAM. Bioremediation of a petroleum hydrocarbon‐contaminated Antarctic soil: Optimization of a biostimulation strategy using response‐surface methodology (RSM). Cold Regions Science and Technology. 2015;119:61-67. DOI:10.1016/j.coldregions.2015.07.005
27 - Li Y‐L, Fang Z‐X, You J. Application of Box‐Behnken experimental design to optimize the extraction of insecticidal Cry1Ac from soil. Journal of Agricultural and Food Chemistry. 2013;61:1464-1470. DOI: 10.1021/jf304970g
28 - Box GEP, Behnken DW. Some new three level designs for the study of quantitative variables. Technometrics. 1960;2:455-475. DOI: 10.1080/00401706.1960.10489912
29 - Dimović S, Smičiklas I, Šljivić‐Ivanović M, Dojčinović B. Speciation of 90Sr and other metal cations in artificially contaminated soils: The influence of bone sorbent addition. Journal of Soils and Sediments. 2013;13:383-393. DOI: 10.1007/s11368‐012‐0633‐7
30 - Smičiklas I, Dimović S, Plećaš I. Removal of Cs1+, Sr2+ and Co2+ from aqueous solutions by adsorption on natural clinoptilolite. Applied Clay Science. 2007;35:139-144. DOI: 10.1016/j.clay.2006.08.004
31 - Brown G, Parks G. Sorption of trace elements onto mineral surfaces: Modern perspectives from spectroscopic studies, and comments on sorption in the marine environment. International Geology Review. 2001;43:963-1073. DOI: 10.1080/00206810109465060
32 - Xu N, Hochella MF, Brown GE, Parks GA. Co(II) sorption at the calcite‐water interface: I. X‐ray photoelectron spectroscopic study. Geochimica et Cosmochimica Acta. 1996;60:2801-2815. DOI: 10.1016/0016‐7037(96)00133‐0
33 - USEPA. Selecting remediation techniques for contaminated sediment. Report EPA‐823‐B93‐001. Washington, DC, Cincinnati; 1993
34 - Bradl HB. Adsorption of heavy metal ions on soils and soils constituents. Journal of Colloid and Interface Science. 2004;277:1-18. DOI: 10.1016/j.jcis.2004.04.005
35 - Milenković A, Smičiklas I, Bundaleski N, Teodoro OMND, Veljović Ð, Vukelić N. The role of different minerals from red mud assemblage in Co(II) sorption mechanism. Colloids and Surfaces A. 2016;508:8-20. DOI: 10.1016/j.colsurfa.2016.08.011
36 - Al Abdullah J, Al Lafi AG, Al Masri W, Amin Y, Alnama T. Adsorption of cesium, cobalt, and lead onto a synthetic nano manganese oxide: Behavior and mechanism. Water, Air, & Soil Pollution. 2016;227:241. DOI: 10.1007/s11270‐016‐2938‐4
37 - Mukherjee J, Ramkumar J, Chandramouleeswaran S, Shukla R, Tyagi AK. Sorption characteristics of nano manganese oxide: Efficient sorbent for removal of metal ions from aqueous streams. Journal of Radioanalytical and Nuclear Chemistry. 2013;297:49-57. DOI: 10.1007/s10967‐012‐2393‐7
38 - Pourret O, Lange B, Houben D, Colinet G, Shutcha M, Faucon M‐P. Modeling of cobalt and copper speciation in metalliferous soils from Katanga (Democratic Republic of Congo). Journal of Geochemical Exploration. 2015;149:87-96. DOI: 10.1016/j.gexplo.2014.11.011
39 - Zachara JM, Cowan CE, Resch CT. Sorption of divalent metals on calcite. Geochimica et Cosmochimica Acta. 1991;55:1549-1562. DOI: 10.1016/0016‐7037(91)90127‐Q
40 - Rieuwerts JS, Thornton I, Farago ME, Ashmore MR. Factors influencing metal bioavailability in soils: Preliminary investigations for the development of a critical loads approach for metals. Chemical Speciation and Bioavailability. 1998;10:61-75. DOI: 10.3184/095422998782775835
41 - Harter RD. Competitive sorption of cobalt, copper, and nickel ions by a calcium‐saturated soil. Soil Science Society of America Journal. 1992;56:444-449. DOI: 10.2136/sssaj1992.03615995005600020017x
42 - Kuo S, Lai MS, Lin CW. Influence of solution acidity and CaCl2 concentration on the removal of heavy metals from metal‐contaminated rice soils. Environmental Pollution. 2006;144:918-925. DOI: 10.1016/j.envpol.2006.02.001
43 - Dermont G, Bergeron M, Mercier G, Richer‐Laflèche M. Soil washing for metal removal: A review of physical/chemical technologies and field application. Journal of Hazardous Materials. 2008;152:1-31. DOI: 10.1016/j.jhazmat.2007.10.043
44 - Wuana RA, Okieimen FE, Imborvungu JA. Removal of heavy metals from a contaminated soil using organic chelating acids. International Journal of Environmental Science and Technology. 2010;7:485-496. DOI: 10.1007/BF03326158
45 - Wen J, Stacey SP, McLaughlin MJ, Kirby JK. Biodegradation of rhamnolipid, EDTA and citric acid in cadmium and zinc contaminated soils. Soil Biology and Biochemistry. 2009;41:2214-2221. DOI: 10.1016/j.soilbio.2009.08.006
46 - Kim GN, Jung YH, Lee JJ, Moon JK, Jung CH. An analysis of a flushing effect on the electrokinetic‐flushing removal of cobalt and cesium from a soil around decommissioning site. Journal of Hazardous Materials. 2008;152:1-31. DOI: 10.1016/j.jhazmat.2007.10.043