Soils Developed from Dolomitic Shale in the Yichang Area, China and Adsorption Characteristics for Phenol ()
1. Introduction
Phenol, also known as carbolic acid, is a common contaminant in the chemical industrial wastewater. As a weakly ionizable aromatic organic compound, phenol is notorious for its acute toxicity on aquatic organisms [1], potential cancerogen, mutagenicity and teratogenicity on human beings [2] [3], as well as refractory biodegradation with respect to wastewater treatments [4]. The discharge of phenolic effluents into water bodies without a prior satisfactory treatment would pose significant public health risk as well as serious threat to the aquatic ecosystems. In China, the Ministry of Ecology and Environment (MEE) has set a discharge standard of less than 0.3 ppm of phenol in the effluents from municipal wastewater treatment plants [5]. As a result, effective treatment of phenolic effluents is imperative before reuse or discharge [3]. There are a number of water purification techniques for treating phenolic effluents, such as adsorption [6] [7] [8], oxidation [9] [10] [11], extraction [12] [13], membrane separation [14] [15], photocatalytic degradation [16] [17], biological treatment [3] [18], and electrodialysis [19]. Among all the mentioned techniques, adsorption is one of the simplest, effective and economical methods for the removal of low concentrations of organic pollutants from large volumes of potable water, process effluents, wastewater, and aqueous solutions [20] [21].
Over the last century, activated carbon (AC) is consistently proved to be effective for the removal of phenol and its derivatives from industrial effluents [22] [23]. However, the disadvantage associated with AC lies in, the need for a costly regeneration system, the generation of carbons fines due to the brittle nature of carbons as well as its high initial manufacture cost [6] [7]. Furthermore, the production of AC generally involves high energy consumption processes and CO2 emissions, which make it difficult to meet the needs of low economic development and led many investigators to search low-cost adsorbents to remove phenols from wastewater.
As natural raw materials, minerals usually have thin or layered structure, and the space between layers helps them to act as carries to store gas [24] [25] or absorb organic [26] [27]. Shale is a fine-grained sedimentary rock and also a carrier of shale gas. Most previous studies on shale were conducted with a focus on methane adsorption and/or desorption behavior on shale matrix [24] [25]; however, there are few reports on the removal of phenolic contaminants by soils originated from natural shales. In this paper, soils developed on dolomitic shale (Yichang area, China) were used without any pretreatment to adsorption of phenol from simulated phenolic wastewaters. Effects of operating temperature, adsorbent dosage, contact time and pH on the phenol removal were experimentally investigated. Laboratory batch kinetics and isotherm studies as well as adsorption thermodynamics were conducted to evaluate the adsorption capacity of the soils.
2. Experimental
2.1. Chemicals and Soil Samples
Phenol, 4-aminoantipyrine, potassium ferricyanide, hydrochloric acid, sodium hydroxide and other inorganic chemicals were supplied by Aladdin Industrial Corporation (Shanghai, China), as analytical grade reagents. Shale outcrops of the Ediacaran Doushantuo Formation are distributed across a large part of South China [28]. The studied soils were developed on shale in the Yichang area, South China (as shown in Figure 1). All samples were air-dried, crushed with a mechanical rock crusher, ground in a ceramic mill, and then sieved to the desired particle size of ~180 μm. Generally, the studied soils were kept in a desiccator over silica gel desiccant until use in characterization and adsorption experiments.
Figure 1. Geological map of the Yichang area, South China. 1—Quaternary; 2—Palaeocene-Neocene; 3—Cretaceous; 4—Jurassic; 5—Triassic; 6—Devonian-Permian; 7—Silurian; 8—Cambian-Ordovician; 9—Nanhua-Sinian; 10—Neoproterozoic granite; 11—Syncline axis; 12—Normal fault; 13—Unidentified fault; 14—Boundary of the structural unit; 15—Well name; 16—Area name; 17—Sampling site of shale. F1—Tongcheng river fault; F2—Wudu river fault; F3—Tianyangping fault; F4—Xiannüfeng fault [29].
2.2. Analytical Method
The powder X-ray diffraction (XRD) data were collected at room temperature with a PANalytical X’Pert PRO X-ray diffractometer (PANalytical Co., Almelo, Netherlands). The 2θ Bragg angles were scanned over a range of 5˚ - 85˚ at a scanning speed of 6.25˚/min. To obtain a quantitative measurement of the mineral components of a given sample, the maximum intensity of each identified mineral was measured and compared to the standard intensity obtained from a pure mineral sample. The FT-IR spectra were obtained in a diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) mode using the micro sampling cup of a Spectra-Tech diffuse reflectance accessory against a KBr background on a Nicolet 750 FT-IR spectrometer (Thermo Nicolet Corporation, Waltham, MA, DTGS detector; Nichrome source; KBr beam splitter). Morphology of the shale soils was observed by a field emission-scanning electron microscope (FE-SEM, Model MIRA3-XMU, TESCAN Czech Republic) operated at an acceleration voltage of 15 kV. The N2 adsorption-desorption isotherms can provide information on the total pore volume, specific surface area and pore size distribution of the shale soils. Samples were degassed at 140˚C under a pressure of 0.334 MPa for 2 h in nitrogen atmosphere before testing.
2.3. Batch Adsorption Experiments
In batch adsorption experiments, accurately weighed shale soils (1.0 - 3.0 g) and 50 mL of the aqueous solution containing different initial concentrations of phenol (10 - 300 mg/L) were put into a glass-stoppered flask. The flask was darkly brown-colored to prevent photooxidation. The pH of the solution was adjusted with HCl or NaOH solution by using a pH meter (Model PB-10, Shanghai Cany Precision Instrument Co., Ltd., China). The flask was subsequently capped with the glass stopper and shaken in a temperature-controlled water bath shaker (Model SHA-B, Tianjin Saidlis Experimental Analytical Instrument Factory, China) for 2 h at the temperature range of 20˚C - 60˚C. At preset contact times, the concentrations of phenol in the solution were analyzed for residual concentration of phenol at λmax = 500 nm using a 722N UV-vis spectrophotometer (Shanghai AuCy Scientific Instrument Co., Ltd, China) by the 4-aminoantipyrene method [30]. The amount adsorbed was calculated from the concentration of the phenol solution before and after the adsorption experiments using a Beer’s law plot to interpolate concentrations. Phenol uptake at equilibrium, qe (mg/g), was calculated by the equation, qe = (C0 − Ce) V/W. The C0 and Ce (mg/L) are the concentrations of phenol at initial and at equilibrium, respectively. V is the volume of the solution (L) and W is the mass of dry adsorbent used (g). The calibration curves between absorbance and the concentration of the phenol solution were established. The calibration plot of absorbance versus concentration of phenol showed a linear variation. The standard curve was given as y = 0.0919x − 0.0002, R2 = 0.992, where y is the absorbance and x is the phenol concentration. All experiments were run in triplicates under identical conditions to ensure reproducibility, and sometimes repeated again and the mean values have been reported. Standard deviations were found to be within ±3%.
3. Results and Discussions
3.1. Characterization of the Shale Soils
The FT-IR spectra of the shale soil were presented in Figure 2. The broad and strong band in the range of 3150 - 3557 cm−1 was assigned to the -OH stretching vibration in hydroxyl groups of constitution water and interlayer water of clay minerals. The peaks at 1420 cm−1, 872 cm−1 and 720 cm−1 were denoted as the stretching vibration of C-O in the carbonate anion (CO32−) of dolomite [31]. The strong band in the range of 1000 - 1120 cm−1 was related to the Si-O vibration in quartz and/or feldspar [32]. The peaks at 456 - 1000 cm−1 were assigned to the vibration of Si-O, Al-O-H and Si-O-Al [32]. The results of FTIR are further proved by XRD, as shown in Figure 3. The samples are mainly composed of mineral-like materials, and the quantitative mineralogical analyses demonstrated that the typical shale soil was mainly composed of dolomite (41.12 wt%), quartz (28.78 wt%), kaolinite (23.03 wt%), analcime (2.30 wt%) and plagioclase (4.76 wt%). This shale sample collected from outcrops of the Ediacaran Doushantuo Formation, south China is classified as dolomitic shale [29]. The SEM micrograph in Figure 4 was used to assess the morphological structure of the shale soil. It was demonstrated that the studied shale soil exhibited some visible fractures and layered structures with different diameters (Figure 4), possibly owing to the disintegration of the shale matrix. The well-developed structure of the shale soil was expected to potentially benefit absorption of phenols when it was applied to phenolic wastewater treatment.
The BET specific surface areas of the shale soil were determined as 19.287 m2/g. The low-temperature nitrogen adsorption-desorption isotherms of the shale soil were shown in Figure 5, taking on an invert “S” shape. The low-temperature nitrogen adsorption isotherms of the shale soil (Figure 5) exhibited a BET type
Figure 2. FT-IR spectra of the shale soil.
Figure 3. X-ray diffraction patterns of the shale soil.
Figure 4. The SEM micrograph of the shale soil.
Figure 5. Nitrogen adsorption-desorption isotherms of the shale soil.
IV [33]. Within the lower range of 0 - 0.05 for the relative pressure p/p0, the adsorption isotherms of the shale soil protruded slightly upward, indicating monolayer adsorption occurring on the surface of the shale soil (Figure 5). When the relative pressure gradually increased to 0.05 - 0.6, this change in slope of adsorption isotherm branch may be considered indicative of transition from monomolecular to multimolecular layer. When p/p0 > 0.8, the adsorption branch ascended sharply, but did not reach a stable state even when p/p0 ≈ 1 (Figure 5). This implied that nitrogen vapor confined in the pore of sample particles condensed at a pressure lower than its saturation pressure, which was called “capillary condensation”. This phenomenon indicates the presence of a certain amount of mesopores and macropores in the studied shale soils [34] [35] [36]. As evidenced in Figure 5, desorption branches did not completely coincide with the corresponding adsorption branches, especially at the relatively high-pressure section (p/p0 > 0.6). The adsorption and desorption isotherms formed a hysteresis loop. This hysteresis loop was classified as type H3 according to IUPAC classification [33], implying micropores developed in the shale soil [37]. Pore size distributions (PSDs) were calculated using the BJH (Barrett-Joyner-Halenda) method and shown in Figure 6. It was found that the PSDs of the shale soil was mainly limited within the range of 1.78 - 5.54 nm and peaked at 2.74 nm (Figure 6). An average pore size of 10.99 nm was obtained by averaging BJH values from the isotherms.
Figure 6. Pore size distribution of the shale soil.
3.2. The Orthogonal Test of Phenol Removal by Natural Shale Soils
The orthogonal test (Table 1) with four factors and five levels was designed to analyze the main influencing factors on the phenol removal by the shale soil. The selected factors included adsorbent dose, contact time, operating temperature and pH. It was found that the operating temperature and adsorbent dose with range values of 31.16 and 12.98, respectively, were the key factors affecting the phenol adsorption on the shale soil (Table 2). pH did not significantly affect the phenol removing within the scope of 3.0 - 11.0. Furthermore, the contact time of 60 min was sufficient to assure progress to the adsorption of phenol by the shale soil, and the extra time would be not needed. These contact time and pH showed subordinate factors on the phenol removing, which can be proved from their low range values of 7.31 and 3.96, respectively (Table 2). Therefore, the order of influencing factors on the phenol removal can be obtained as follows: operating temperature > adsorbent dosage > contact time > pH.
Table 1. The orthogonal test scheme of 4 factors and 5 levels.
Table 2. Range analysis of the orthogonal test.
Note: K: sum of each level; Rang = Kmax − Kmin; removal efficiencies are average values.
3.3. Adsorption Kinetics
The adsorption kinetic characteristics are important for providing valuable insights into the adsorption mechanisms [38]. The experimental kinetic data of adsorption of phenols onto the shale soil were examined using the pseudo-first-order and pseudo-second-order equations [39]. The pseudo-first-order kinetic equation is expressed by
(1)
The pseudo-second-order kinetic equation is expressed by
(2)
where qe and qt are the amount of phenol adsorbed per unit mass of the adsorbent (mg/g) at equilibrium and time t respectively. k1 (min−1), and k2 (g(mg·min)) are the pseudo-first-order and pseudo-second-order rate constant, respectively. The fitting results of the adsorption kinetics were shown in Figure 7 and Table 3.
Figure 7. Kinetic plots of phenol adsorption on shale soil: (a) Pseudo-first-order model; (b) Pseudo-second-order model.
It was found that the correlation coefficient (R2) values of the pseudo-second-order were higher than that of the pseudo-first-order kinetics model (Table 3). This demonstrated that the pseudo-second-order model fitted well with the adsorption experimental data and can be used for the entire adsorption process (Figure 7). Our kinetic results for the adsorption of phenol by the shale soil are consistent with some previous studies using natural soil [40] and natural zeolite [41]. The overall rate of adsorption may involve surface diffusion, pore diffusion and adsorption on the interior sites of the sorbent [42]. The kinetic process that is fitted by the pseudo-second-order kinetic model has chemisorption as the rate-determining step and involves valency forces by sharing or through the exchange of electrons between the shale soils and phenols as covalent forces [43]. The actual adsorption sites of the soils may be associated with the organic matter within the shale soil, in which case phenol would be attracted by their similar aromatic character. Generally, the adsorption capacity of the natural shale soil for phenol increases with the increasing time and becomes constant after the equilibrium time is reached. Evidently, a large number of active sites are available on the surface of the shale soil in the initial stages, resulting in its high phenol adsorption capacity and the occurrence of rapid adsorption. Thereafter, the uptake of phenols becomes slower before a dynamic equilibrium is attained. This may be due to the gradual diminishing of the surface sites and the increasing steric repulsion among the phenol molecules on the soil surface and the incoming ones from the solution.
Table 3. Parameters of adsorption kinetics.
3.4. Adsorption Isotherms
Several isotherm models have been reported [44] to describe experimental data of adsorption equilibriums including Langmuir, Freundlich, Redlich-Peterson, Sip isotherms, etc. The Langmuir isotherm is the most widely used two-parameter equation, and assumes the following hypotheses: uniform adsorption energies along the homogeneous adsorbent surface, equal solute affinity in all the adsorption sites, no interaction among adsorbed molecules, single adsorption mechanism, and formation of a monolayer on the free surface [22] [30] [45] [46] [47]. The Freundlich isotherm is used for modelling the adsorption on heterogeneous surfaces. This experimental model can be applied to nonideal sorption on heterogeneous surfaces as well as multi-layer sorption [41] [48] [49] [50]. The Redlich-Peterson isotherm model is an empirical isotherm incorporating three parameters. It combines elements from both the Langmuir and Freundlich equations, and the mechanism of adsorption is a hybrid and does not follow ideal monolayer [40]. The Sip model equation follows Freundlich isotherm at lower solute concentration and obeys Langmuir isotherm at higher solute concentration [40]. A statistical analysis of the adsorption isotherms reveals that Langmuir and Freundlich are the most examined models in fitting the experimental adsorption data [8]. Accordingly, in this study the sorption data of phenols were correlated with these two isotherms (Equations (3) and (4)).
Langmuir equation:
(3)
where Ce is the equilibrium concentration of phenol in solutions (mg/L), qe is the amount of phenol adsorbed per unit mass of the adsorbent at equilibrium state (mg/g), qmax is the maximum adsorption capacity (mg/g). KL is the Langmuir equilibrium constant related to the affinity of binding sites and energy of adsorption (L/(mg·min)). KF is the Freundlich constant (L/(mg·min)), indicating the adsorption capacity of the adsorbent and n (dimensionless) is a constant related to the intensity of adsorption.
Freundlich equation:
(4)
The fitting results of the adsorption isotherms were shown in Figure 8 and Table 4. The R2 values of the Freundlich model were generally higher than those
Figure 8. Adsorption isotherms of phenol on the shale soil: (a) Langmuir model; (b) Freundlich model.
Table 4. Parameters of adsorption isotherm models at various temperatures.
of the Langmuir model for the adsorption of phenol onto the natural shale soil (Table 4). Evidently, the Freundlich model can be applied to describe experimental data of adsorption equilibriums during phenol removal by the shale soil. This is in accordance with the heterogeneity of the shale soil structures characterized by FE-SEM (Figure 4). Moreover, the increasing values of qe and n with the increase of temperatures in Table 4 indicate that phenol removal by the shale soil is favored at elevated temperatures under the constrained batch experiments [51]. Nayak and Singh [30] presented that higher temperatures promote the rupture of hydrogen bonds between the phenol molecule and water, causing more phenol adsorbed by the adsorbent rather remaining in the aqueous solution.
3.5. Adsorption Thermodynamics
The feasibility of adsorption of phenols on the shale soil surface can be evaluated by exploring the thermodynamics parameters. Energetic changes during adsorption in terms of considering amount of heat being absorbed or released, dispersal of energy and amount of energy associated with the adsorption are quantified by standard enthalpy (Ho), standard entropy (So) and Gibbs free energy (Go) respectively. These parameters are essential for indicating the spontaneity and/or possibility of adsorption mechanism and defining the exothermic and endothermic adsorption processes [8].
According to Van’t Hoff Equation (5), the Gibbs free energy change of adsorption ΔG˚ (kJ/mol) for the adsorption of phenol onto the shale soil was calculated using Equation (5) and shown in Table 5.
(5)
where R is the universal gas constant (8.314 J/(mol·K)), T is the absolute temperature and KF is the Freundlich constant. Temperature T was plotted against ΔG˚ based on the following basic thermodynamic equation:
(6)
The values of the enthalpy (ΔH˚, kJ/mol) and entropy (ΔS˚, J/(mol·K)) could be obtained from the slope and intercept of the line in Figure 9 and then listed in Table 5.
Values of ΔG˚ are significantly smaller than 50 KJ/mol (Table 5). This indicates that the adsorption of phenol onto the studied shale soil surface is thermodynamically possible. Moreover, the magnitude of ΔGo decreased with an increase in the temperature denoting that the adsorption is more favorable at higher temperature range of 25˚C - 45˚C. The value of ΔHo is positive, implying the endothermic nature of the adsorption of phenol onto the shale soils. One possible explanation of endothermicity of the enthalpy of adsorption is due to the reason that hydrogen bonds were formed between phenol and water molecules [52]. Phenols would not be adsorbed by the shale soil until the hydrogen bonds were destroyed. This rupture of hydrogen bonds between phenol and water molecules obviously required energy [53]. Therefore, it is reasonably assumed that this energy of hydrogen-bond breaking exceeded the exothermicity of the phenol attaching to the surface of the shale soil. The positive values of ΔSo indicate increase in randomness at the interface of the solid soil/liquid water during the adsorption process and show the good affinity of the soils for the adsorbed phenols. It should be noted that a chemisorbed layer could be followed by the physically adsorbed layer which may not be excluded in this study.
Table 5. Thermodynamic parameters for adsorption.
3.6. Comparison with Other Studies Using Natural Materials
Various chemical, thermal and/or biological pretreatments had been utilized to improve the adsorption capacity of the low-cost materials for phenol removal. The modification methods always resulted in the increase of manufacturing cost as well as difficulties in regeneration/disposal of the adsorbents [40]. Notably, some previous researchers investigated the adsorption of phenol by naturally occurring adsorbents under the constrained experimental conditions. Typical experimental results were listed in Table 6. Isotherm models are selected to report based on their fitting performance with the experimental data. As shown in Table 6, the Langmuir and Freundlich models are the most frequently employed to describe the relationship between the amount of phenol adsorbed by naturally occurring adsorbents and its equilibrium concentration in solutions under different experimental conditions. The present study revealed that the Freundlich adsorption isotherm model could better describe the adsorption of phenol on the shale soil surface than the Langmuir adsorption isotherm model. Our finding is in agreement with some previous studies using natural bentonite [49] [54], clinoptylolite [55], siltstone [48], raw lignite [50], kaolinite [47] [56], spent oil shale [57], soil [56], peat [58], natural zeolite [41] and zeolitic tuff [41]. In comparison with studies on adsorption isotherms, adsorption kinetics for phenol removal by natural materials were far less documented. It can be summarized from Table 6 that adsorption mechanism of phenol is best described by pseudo-second-order [40] [41], pseudo-first-order [22] [50] and Elovich equations [59]. Considering different experimental conditions (e.g., temperature, contact
Table 6. Comparison of phenol adsorption on various naturally occurring adsorbents.
time, pH, concentration ranges of phenol) utilized in adsorption experiments, it is not ready to compare adsorption performance of various natural materials. However, the studied shale soils showed a relatively higher equilibrium capacity (of an order of magnitude of 10 mg/g, Table 6) than those of other natural materials including clinoptylolite [55], siltstone [48], spent oil shale [57], granite residual soil [56], clays [47] [54] [56] [58] and coals [58] [60].
It should be pointed out that direct utilization of the studied shale soils without surface modification does not seem viable for the effective treatment of high concentration phenolic effluents from chemical and allied industries (200 - 2000 mg/L) [22]. Nevertheless, the studied soils, which are developed on the natural shale widespread in South China, are promising candidates for advanced removal of trace phenol from mass production of potable and mineral waters, considering the limit for phenols in potable and mineral waters is 0.5 μg/L (0.5 ppb) recommended by the European Union [61].
4. Conclusion
In the present study, laboratory experiments were performed to evaluate the use of shale soils as an adsorbent for phenols in aqueous solution. It was found that the factors affecting the adsorption of phenol by the shale soil are temperature > dosage > adsorption time > pH. The studied shale soils exhibited small fragments and well-developed pore structures with different diameters. The pore size distributions of the samples were mainly limited within the range of 1.78 - 5.54 nm and peaked at 2.74 nm. Thermodynamic analysis showed that the adsorption of phenol by the shale soil is the endothermic process, and the higher temperature in the range of 25˚C - 45˚C favored the adsorption process. The entire adsorption kinetic process follows the pseudo-second-order model, which has chemisorption as the rate-determining step. Investigation of adsorption thermodynamics indicates that Freundlich isotherm model gives a better description for the isotherm curve than the Langmuir model, consistent with the heterogeneity of shale soils in physical-chemical properties. Results of our study indicate potential utilization of natural shale soils as a promising low-cost absorbent for trace phenol removal from mass production of potable and mineral waters. Utilization of naturally occurring adsorbents undoubtedly benefits the maintenance of a healthy surrounding environment as well as the sustainable development of economy.
Acknowledgements
We gratefully acknowledge financial support by the National Natural Science Foundation of China (No. 41472095), the Open Fund of Key Laboratory of Exploration Technologies for Oil and Gas Resources (Yangtze University), Ministry of Education (No. K2018-05) and the Open Fund of Shandong Provincial Key Laboratory of Depositional Mineralization & Sedimentary Minerals, Shandong University of Science and Technology (Grant No. DMSM2018041).