Next Article in Journal
Pesticide Importation in Sierra Leone, 2010–2021: Implications for Food Production and Antimicrobial Resistance
Next Article in Special Issue
Application of Phosphate Materials as Constructed Wetland Fillers for Efficient Removal of Heavy Metals from Wastewater
Previous Article in Journal
Leaching of Polycyclic Aromatic Hydrocarbons from the Coal Tar in Sewage Wastewater, Acidic and Alkaline Mine Drainage
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Removal of Pb2+ from Aqueous Solution by Using Navel Orange Peel Biochar Supported Graphene Oxide: Characteristics, Response Surface Methodology, and Mechanism

1
School of Civil and Surveying & Mapping Engineering, Jiangxi University of Science and Technology, Ganzhou 341000, China
2
Jiangxi Provincial Key Laboratory of Environmental Geotechnology and Engineering Disaster Control, Ganzhou 341000, China
3
School of Live Sciences, Jinggangshan University, Ji’an 343009, China
4
School of Resources and Environmental Engineering, Jiangxi University of Science and Technology, Ganzhou 341000, China
*
Authors to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(8), 4790; https://doi.org/10.3390/ijerph19084790
Submission received: 25 January 2022 / Revised: 21 March 2022 / Accepted: 24 March 2022 / Published: 15 April 2022

Abstract

:
The value-added utilization of waste resources to synthesize functional materials is important to achieve the environmentally sustainable development. In this paper, the biochar supported graphene oxide (BGO) materials were prepared by using navel orange peel and natural graphite. The optimal adsorption parameters were analyzed by response surface methodology under the conditions of solution pH, adsorbent dosage, and rotating speed. The adsorption isotherm and kinetic model fitting experiments were carried out according to the optimal adsorption parameters, and the mechanism of BGO adsorption of Pb2+ was explained using Scanning Electron Microscope (SEM-EDS), X-ray Photoelectron Spectroscopy (XPS), X-ray Diffraction (XRD), and Fourier Transform Infrared Spectroscopy (FTIR). Compared with virgin biochar, the adsorption capacity of Pb2+ on biochar supported graphene oxide was significantly increased. The results of response surface methodology optimization design showed that the order of influence on adsorption of Pb2+ was solution pH > adsorbent dosage > rotating speed. The optimal conditions were as follows: solution pH was 4.97, rotating speed was 172.97 rpm, and adsorbent dosage was 0.086 g. In the adsorption–desorption experiment, the desorption efficiency ranged from 54.3 to 63.3%. The process of Pb2+ adsorption by BGO is spontaneous and endothermic, mainly through electrostatic interaction and surface complexation. It is a heterogeneous adsorption process with heterogeneous surface, including surface adsorption, external liquid film diffusion, and intra-particle diffusion.

1. Introduction

Excessive amounts of heavy metals released by the geochemical cycle [1,2] and human production activities (domestic sewage [3], industrial production [4], mining [5,6]) pose a major threat to the ecological environment system and may have a negative impact on the health of humans, plants, and microorganisms [7,8]. At present, there is no effective degradation pathway for heavy metals in the water environment, and the most common method is centralized recovery treatment after adsorption and desorption [9,10,11,12]. Due to the requirements of social environment and sustainable development, exploring high efficiency and environmentally-friendly adsorbents has become the current research hotspot [13,14,15]. The biochar extracted from various solid wastes was widely used in various environmental treatment [16,17,18] and has a strong application prospect to achieve the goal of “treating waste with waste”.
Biochar can be understood as a continuum of pyrolysis products with very complex chemical composition and characteristics [19]. It contains a variety of negatively charged functional groups [20,21], which can fix heavy metal ions on the surface of biochar. However, there are some problems in the virgin biochar materials, such as limited specific surface area, low porosity, and few adsorption sites, which limit its wide application [22,23,24,25]. Therefore, studies enhance the adsorption capacity of biochar by modifying and loading functional materials [26,27]. Graphene oxide (GO) is the main derivative of graphene. It has attracted wide attention because of its large specific surface area and abundant oxygen-containing functional groups [28,29]. Huang et al. (2017) [30] found that the pore volume, specific surface area, and adsorption performance of BGO were significantly increased. Chen et al. (2013) [31] found that graphene oxide chitosan composite material has a high treatment capacity for heavy metal pollutants.
There are abundant navel orange resources in South Jiangxi, which are deeply loved by people. However, the navel orange peels have low economic value and are often discarded, which releases a large amount of CO2 and causes environmental pollution [32]. The navel orange peel contains abundant plant fibers and functional groups, which is a good raw material for preparing biochar. Therefore, this study selected lead, a common and typical heavy metal pollutant in the water environment, as the adsorption object, and navel orange peel and flake graphite were used as the raw materials to prepare BGO composites. The optimal adsorption parameters under solution pH, rotating speed, and dosage were analyzed by response surface method. Thermodynamic and kinetic models were used to describe the adsorption process of BGO composites. Combined with Scanning Electron Microscope (SEM-EDS), X-ray Photoelectron Spectroscopy (XPS), X-ray Diffraction (XRD), and Fourier Transform Infrared Spectroscopy (FTIR), the adsorption mechanism of BGO for Pb2+ in aqueous solution was described. The results will provide some basic experimental data and material information for the effective treatment of heavy metals in sewage and provide new ideas for the resource reuse of navel orange peel.

2. Materials and Methods

2.1. Preparation of BGO Composites

Navel oranges were purchased from Ganzhou fruit wholesale market in Jiangxi Province. The peels were used as the experimental raw material, the graphene oxide was prepared by the modified hummer method [33], and the ratio of flake graphite: KNO3:KMnO4:H2SO4 was 1:1.2:6:46 [34]. The mass ratio of biomass to graphene oxide was determined as 20:1 according to the results of preliminary experiments. About of 73 g of biomass and 200 mL of graphene oxide suspension were put into ultrasonic cleaner for 2 h to obtain a uniformly dispersed and stable brown mixture. It was then freeze-dried for 48 h, and then placed in a tube furnace for pyrolysis. We adopted a slow heating and oxygen-limited pyrolysis conditions, heating to target temperature (300 °C, 500 °C, 700 °C) at a rate of 10 °C/min, and then the pyrolysis was carried out for 2 h. It used a N2 flow rate of 20 mL/min. After the pyrolysis, the samples were passed through a 100-mesh sieve and put into a sample bag for dry storage. According to the different pyrolysis temperature, the samples of virgin biochar were marked as B300, B500, and B700, the graphene oxide were marked as GO300, GO500, and GO700, and biochar supported graphene oxide were marked as BGO300, BGO500, and BGO700. The adsorbent information were shown shown in Table S1.

2.2. Experimental Methods

2.2.1. Single-Factor Adsorption Experiment

The single-factor adsorption experiment was carried out in the water bath constant temperature shaker. The influence of factors such as solution pH, rotating speed, and adsorbent dosage on the adsorption effect were investigated. BGO composite was put into a 100 mL polyethylene centrifuge tube, and 50 mL of Pb2+ solution (100 mg/L) was added. The experimental operating temperature was 25 ± 2 °C, and the shaking speed was set to 150 rpm. In order to achieve adsorption equilibrium, the reaction time was set to 24 h. The pH value was set adjusted to 5.0 ± 0.1 with 0.1 mol/L HNO3 and NaOH solution. This experiment adopted the control variable method—when a single factor changes, other reaction conditions remain unchanged. The design scheme was shown in Table 1. The supernatants were collected through 0.45-μm filter membrane. Finally, the concentrations of Pb2+ in each sample were tested after acidification with HNO3.

2.2.2. Response Surface Experiment

The factor of solution pH, rotating speed, and adsorbent dosage that have a significant impact on the adsorption performance were selected. The experimental design was optimized by Box–Behnken response surface [35]. The central composite design of the adsorption experiment was carried out at three levels of low, medium, and high, respectively. There were 17 groups of experiments, and the central point experiment was repeated for 5 groups. The experiment was repeated three times in each group, and the average value was taken as the corresponding response value.

2.2.3. Adsorption Isotherm and Kinetic Fitting Experiment

The single-molecular-layer Langmuir model, multi-molecular-layer Freundlich model, and diffusion interface Temkin model were used to perform fitting analysis on the adsorption of Pb2+ by BGO700. The Langmuir model can better represent the experimental results when the adsorption on solid surface is quite uniform, and the adsorption is limited to single molecular layer. The Freundlich model and Temkin model can describe the heterogeneous adsorption process of heterogeneous surfaces. The sample was placed into a 100 mL centrifuge tube, then 50 mL of Pb2+ solution with different initial concentration (50, 100, 150, 200, 250, 300, 350, 400, 500 mg/L) was added into each of them. The mixtures were shaken at 25, 35, and 45 °C for 24 h in the water bath constant temperature shaker. The pH of the solution and rotating speed are adjusted to 5 and 173 rpm, respectively, according to the optimal experiment conditions coupled with the response surface. The supernatants were collected through 0.45-μm filter membrane. Finally, the concentrations of Pb2+ in each sample were tested after acidification with HNO3.
The adsorption kinetics of BGO700 and B700 were fitted with quasi-first-order kinetics, quasi-second-order kinetics, and intra-particle diffusion models. The solution of Pb2+ with concentration of 100 mg/L was placed in 250 mL conical flask, then shaken in the water bath constant temperature shaker at 25 °C. The adsorbent dosage, solution pH, and rotating speed were adjusted to 0.086 g, 5 and 173 r/min. The water sample was taken out at 1, 5, 10, 30, 60, 90, 120, 240, 360, 600, 900, 1260, 1440, 2160, and 2880 min, respectively. The supernatants were collected through 0.45-μm filter membrane. Finally, the concentrations of Pb2+ in each sample was tested after acidification with HNO3.

2.2.4. Batch Adsorption–Desorption Experiments

Adsorption experiment was first tested to get saturated adsorbent. An amount of 0.086 g BGO700 adsorbent was added to 100 mL of 100 mg/L Pb2+ solution in conical flask. The mixed solution was stirred in a water-bath constant temperature shaker for 24 h and then the pollutant concentration was measured. The saturated adsorbent was collected by vacuum suction filtration. In the desorption experiment, 0.5 mol/L HCl solution was selected as the desorbent [36,37]. The saturated adsorbent was put into 0.5 mol/L HCl solution, the mixed solution was stirred for 12 h, and then the pollutant concentration was measured. The adsorbent of desorption experiment was rinsed with deionized water for several times until the pH of the filtrate was close to neutral. Next, the samples were dried in two different conditions, into the vacuum oven and under an ultraviolet light. The temperature of the vacuum oven was set to 60 °C for 4 h. The ultraviolet drying adopted 256 nm, intensity of 80 W, and height of 10 cm to illuminate for 4 h. The adsorption and desorption experiments were repeated for five times in the same way. There are four groups in this experiment, which are marked as UV/BGO700, UV/B700, Bake/BGO700, Bake/B700, respectively.

2.3. Analytical Instruments and Methods

The pH value of solution was tested in the potential method by using a pH meter (Sartorius, PB-10, Göttingen, Germany). The CHNS content of material was analyzed by an element analyzer (EA, vario El cube, Hanau, Germany), and the O element content was calculated by difference method. The microscopic morphology and surface composition of the sample were characterized by a field emission scanning electron microscope (FEI, MLA650F, Hillsboro, OR, USA) with an energy depressive spectroscopy (Bruker, Karlsruhe, Germany). The surface area and pore volume of the samples were determined by Brenner Teller (ASAP, 2020HD88, Norcross, GA, USA). The surface functional groups of the samples were characterized by Fourier Transform Infrared Spectroscopy (NICOLET5700, Waltham, MA, USA) with the wave number range of 4000~400 cm−1. The changes of crystal structure and functional groups before and after adsorption were analyzed by X-ray diffraction (DX-2700, Dandong, Liaoning, China) and X-ray photoelectron spectroscopy (Thermo Scientific K-Alpha+, Waltham, MA, USA). The concentrations of Pb2+ in solution were determined by using ICP-OES (AvioTM200, Fremont, CA, USA).

2.4. Quality Control

In this study, to ensure the quality of the sample analysis, before measuring the concentrations of Pb in solution, a set of Pb standard solution of different concentrations (0, 1.0, 2.0, 5.0, 10, 20, 50 mg/L) was used to correct the ICP-OES and obtain a calibration curve (R2 > 0.999). All samples were tested three times, and the error of the parallel determination results was within 10%.

2.5. Data Processing and Analysis

The thermodynamic parameters of adsorption mainly include Gibbs free energy (ΔG), enthalpy change (ΔH), and entropy change (ΔS).
The thermodynamic parameters of adsorption can be calculated by the following equation,
∆G = −R·T·ln(KL)
∆G = ∆H − T∆S
From Equations (1) and (2), the following Equation (3) can be obtained:
ln(KL)= ∆S/R − ∆H/(R·T)
where R stands for general gas constant, 8.314 J/(mol·K); T stands for thermodynamic temperature, K; KL stands for solid–liquid partition coefficient, the ratio of the concentration of Pb2+ at adsorption equilibrium (Qe) by BGO700 to the concentration of Pb2+ in the solution (Ce), is linearly fitted by taking Qe as the abscissa and ln(Qe/Ce) as the ordinate, and the intercept is ln(KL).

3. Results and Discussion

3.1. Characterization of Biochar Composites

3.1.1. The Analysis of Physical Properties

The physical properties and specific surface area analysis of biochar composites is shown in Table 2. The pH and C content are positively correlated with the pyrolysis temperature, while the H and O contents were decreased with the increase of pyrolysis temperature. The oxygen content was increased after loading graphene oxide, which indicates that the surface oxygen containing functional groups were increased [38]. The atomic ratios H/C, O/C, and (O + N)/C of the composites were decreased with the increase of the pyrolysis temperature, which indicates that the higher degree of carbonization and more complete π-conjugated aromatic structure with high pyrolysis temperature [39]. The specific surface area was increased, which indicates that the pyrolysis temperature is the key factor affecting the specific surface area.

3.1.2. Result of SEM-EDS

The surface micromorphology and component of B700, GO700, and BGO700 is shown in Figure 1. It can be found from Figure 1(a1,a2) that biochar is porous and has abundant specific surface area, and contains some mineral elements of Ca, Mg. It can be found from Figure 1(b1) that GO is a transparent flocculent structure with smooth surface, curly edge, and fold, and has good flexibility. The surface components mainly contain C and O elements, and the content of O elements is relatively high. The GO supported on biochar was shown Figure 1(c1), and there are smooth and folded flakes on the surface. At the same time, the percentage of O element content increased (Table 2), which indicates that the GO was successfully supported on the biochar.

3.1.3. Result of FT-IR

The FT-IR analysis of B700, GO700, and BGO700 is shown in Figure 2. It can be seen from the figure that the peak at 3434 cm−1 is the characteristic peak of −OH group [40]. The peak of biochar at 2922.8 cm−1 is caused by the stretching vibration of aliphatic hydrocarbons or cycloalkanes −CH3 and −CH2. After being compounded with GO, the composite material also showed a peak at 2357.49 cm−1, where it is unique to GO, and indicates that GO was successfully loaded on biochar. The wavelength around 1620.7 cm−1 is the peak of −COOH group [41]. Moreover, before 1700 cm−1 the composite material showed alkyl, aromatic, and some oxygen-containing groups. These functional groups can provide adsorption sites and enhance the adsorption performance of the composite material.

3.2. Single-Factor Adsorption Experiment of Material

3.2.1. Adsorption Capacity of the Material

The adsorption capacity of Pb2+ by different pyrolysis temperatures was shown in Figure 3. It can be seen from the figure that the pyrolysis temperature has a greater influence on the adsorption capacity of the material. The adsorption capacity of BGO700 was increased by 134.9% compared to BGO300, and B700 was increased by 187.1% compared to B300. This is due to the influence of pyrolysis temperature, which is consistent with the results of other studies [42,43]. In addition, orange peel biochar showed the worst adsorption capacity (B300: 21.7 mg/g; B500: 43.43 mg/g; B700: 62.3 mg/g), which was significantly improved by about 97%, 78%, and 62% when biochar was supported by the GO (BGO300: 42.82 mg/g; BGO500: 77.46 mg/g; BGO700: 100.61 mg/g). This indicates that GO has strong adsorption performance and can effectively improve the ability of the navel orange peel biochar for heavy metal ions, so that biomass resources can be fully utilized.

3.2.2. The Effect of Initial pH on Adsorption Capability

The pH of the solution is an important factor affecting the adsorption of Pb2+, and the electrostatic interaction with Pb2+ is affected by changing the surface charge of BGO. The effect of different pH on the adsorption of Pb2+ by BGO was shown in Figure 4a, which indicated the pH was 6 as the best condition for the adsorption of Pb2+. With the decrease of pH, the adsorption capacity of Pb2+ decreased, showing the lowest value at pH was 2. This is due to the high concentration of H+ in the solution under low pH conditions, which will compete with Pb2+. In addition, the radius of hydrated H+ ion is smaller and easily adsorbed than the ion Pb2+, thus competing for the adsorption site of BGO and inhibiting the adsorption of Pb2+ [44]. With the increase of solution pH, H+ ion decreased and the negative charge on the BGO surface increased, together with the electrostatic attraction between Pb2+ and BGO charged sites.
The morphological distribution of Pb2+ in aqueous solution under different pH conditions is analyzed by visual MINTEQ software. The results are shown in Figure 4b. It can be seen from Figure 4b that when the pH value of the solution is less than 6, Pb exists in ionic state, and with the increase of pH, ionic Pb hydrolyzes to form a variety of hydroxyl complexes. Therefore, when the pH value in the solution is high, Pb2+ is prone to hydrolyze and form Pb(OH)2 precipitation, which is beneficial to the removal of Pb2+ from pollutants. At the same time, the formation of hydroxyl complex will also reduce the average charge of Pb2+, resulting in a significant decrease in the secondary solvation energy, which is more beneficial to the adsorption of Pb2+ on the BGO surface through electrostatic attraction [45].

3.2.3. The Effect of Rotating Speed on Adsorption Capability

The influence of different rotating speed conditions on the adsorption of Pb2+ by BGO is shown in Figure 4c. It can be seen from Figure 4c that the efficiency of Pb2+ adsorption was slightly reduced under low and high-speed conditions. This may be that the diffusion rate of Pb2+ increased with the increase of rotational speed, which leads to an increase in the effective collision probability of the adsorbent BGO and Pb2+ [46]. When the rotating speed is low, the diffusion rate of Pb2+ is low, and the effective collision probability between adsorbent BGO and Pb2+ is low. When the rotating speed was increased to 200 rpm, additional hydrodynamic shear stress will be generated on the surface of BGO, which will destroy the bond bridge between BGO particles and Pb2+ in the gap, increasing the possibility of desorption of BGO adsorbing Pb2+, and reducing the adsorption efficiency of BGO for Pb2+.

3.2.4. The Effect of Adsorbent Dosage on Adsorption Capability

The dosage of adsorbent can determine the adsorption capacity of BGO for Pb2+ at a certain concentration, thereby affecting the adsorption efficiency of BGO for Pb2+ in water. The adsorption capacity of Pb2+ with different dosages of BGO is shown in Figure 4d. The results showed that with increasing of BGO dosage, the adsorption capacity of BGO first increases and then decreases. The reason for the increase may be that the dosage of BGO adsorbent increases, and it is relatively easier to be adsorbed by active sites in the collision process, and the number of adsorption sites is also more. When the content of adsorbent is low, Pb2+ is easy to be separated from the surface of the adsorbent under hydraulic shear. When the content of adsorbent is excessive, the adsorbents are easy to overlap with each other, but the effective active sites cannot be fully utilized.

3.3. Simulation and Optimization of Pb2+ Adsorption Process by Response Surface Methodology

Based on the results of single-factor experiments and environmental factors in engineering practice [47,48], the Box–Behnken method was used to explore the influence of solution pH, rotating speed, and adsorbent dosage on the adsorption capacity of Pb2+ by BGO700 in the coupling system. The Pb2+ adsorption capacity of BGO700 was used as the response value to explain the adsorption effect. The experimental design and results were shown in Tables S2 and S3.
In order to analyze the interaction between the factors and the response value in the experiment, the three-dimensional response surfaces were made according to the regression model equation. We compare the effect of rotating speed and adsorbent dosage against solution pH in Figure 5. The interaction effect between solution pH and rotation speed on Pb2+ adsorption capacity of BGO700 under the central value of adsorbent dosage is shown in Figure 5a. The solution pH has a significant effect on the adsorption of BGO700, while the speed has a small effect. The interaction effect between the solution pH and adsorbent dosage on Pb2+ adsorption capacity of BGO700 under the central value of rotation speed is shown in Figure 5b. The results showed that the dosage of adsorbent has a certain effect on the adsorption capacity of BGO700, but the effect was not significant. The interaction effect between rotation speed and adsorbent dosage on Pb2+ adsorption capacity of BGO700 under the central value of solution pH is shown in Figure 5c. The flat curvature shows that the effect of rotating speed and adsorbent dosage are negligible in comparison to the solution pH, but the influence of adsorbent dosage was slightly higher than the rotation speed. Therefore, the order of the degree of influence on the adsorption of Pb2+ was solution pH > adsorbent dosage > rotation speed. The analysis with the quadratic multiple regression model showed that the optimal conditions of obtaining the maximum Pb2+ adsorption capacity (125.165 mg/g) were: solution pH was 4.97, rotation speed was 172.97 r/min, adsorbent dosage was 0.086 g.
Considering the feasibility of the actual operation, the solution pH, rotation speed, and adsorbent dosage were set to 5, 173 r/min, and 0.086 g, respectively. According to the optimal conditions, three groups of parallel experiments were carried out on the adsorption of Pb2+. The experimental results showed that the adsorption capacity of Pb2+ was 118.32 mg/g, and the deviation from the predicted response value was only 5.47%, which indicates that the error between the actual experimental value and predicted value of model was small. At the same time, combined with the relationship between the actual value and the predicted value of the model, Figure 5d further illustrates that the model was more accurate in predicting the influence factors, and can better optimize the conditions of experiment.

3.4. Adsorption Thermodynamics and Kinetics after Optimization

The adsorption isotherm is a curve that describes the relationship between the dynamic equilibrium process and parameters of the adsorbent in the solution, which is of great significance for evaluating the practicality of adsorption. At present, isothermal models are commonly used in solid–liquid phase adsorption, including Langmuir, Freundlich, and Temkin [18]. The adsorption isotherms of Pb2+ by BGO700 at different temperatures are shown in Figure 6a. When the initial concentration of Pb2+ is the same, the adsorption capacity of Pb2+ by BGO700 is higher at the highest temperature. This is due to the enhanced thermal and diffusion movement of the molecules, which speeds up the movement rate of adsorbent active site and adsorbate.
The fitting results of Langmuir, Freundlich, and Temkin models are shown in Figure 6b–d, respectively. The relevant fitting parameters calculated by the three models are shown in Table 3. It can be seen from the table that the relative coefficient R2 of Langmuir model was less than 0.9, and the theoretical maximum adsorption capacity of Pb2+ fitted by Langmuir model was lower than the actual experimental value. The results showed that there is an interaction between Pb2+ in the solution, the adsorption process was not single molecule adsorption on uniform surface, and BGO700 adsorption process was not suitable for Langmuir model fitting. The Freundlich model and Temkin model can describe the heterogeneous adsorption process of heterogeneous surfaces. The value of 1/n in Freundlich model was between 0.079 and 0.108, which indicates that the adsorption process is easy to take place and multi-layer adsorption can occur. Therefore, the Freundlich model is more suitable for describing the adsorption process.
The adsorption thermodynamic equation and parameters were used to describe the state characteristics of the adsorption process. According to the experimental results of the isotherm adsorption of Pb2+ by BGO700, the linear equation of ln(KL)−1/T was fitted with 1/T as the abscissa and ln(KL) as the ordinate. The parameter fitting results are shown in Figure 7, and the correlation coefficient R2 = 0.989, which has a better linear relationship.
The calculation results of adsorption thermodynamic parameters are shown in Table 4. The calculated values of ∆S > 0, which indicates that the process of adsorption by BGO700 is an entropy increasing process, which may be due to the degree of confusion was increased in the reaction system after the addition of composite materials. The calculated values of ∆G at different temperatures (25 °C, 35 °C, and 45 °C) are negative, ranging from −3.934 to −5.094 kJ/mol, indicating that the adsorption process is spontaneous and endothermic. The absolute value of ∆G was increased with the increase of temperature, indicating that increasing the temperature is beneficial to the progress of adsorption.

3.5. Optimized Adsorption Kinetics Fitting

The variation of Pb2+ adsorption capacity with oscillation time is shown in Figure 8a, from which it can be seen that the adsorption process is divided into three stages: fast adsorption, slow adsorption, and adsorption equilibrium stage. At the initial stage of the reaction, there are many adsorption sites and surface functional groups on the surface of the adsorbent. With the progress of the adsorption reaction, the effective adsorption sites were gradually decreased. Meanwhile, the concentration of Pb2+ in the solution gradually decreases, and the mass transfer rate slows down. At slow adsorption stage, the adsorption rate depends on the speed at which Pb2+ enters the internal site from the outside of the adsorbent [49]. At the last stage, the solid–liquid two–phase reaches dynamic adsorption equilibrium.
The fitting curves of quasi first–order kinetics, quasi second–order kinetics, and intra-particle diffusion model are shown in Figure 8b–d, and the fitting parameters are shown in Table 5. The quasi second–order kinetics equation has a better fitting effect on the Pb2+ adsorption, and the fitting degree reached R2 ≥ 0.998. At the same time, the theoretical equilibrium adsorption capacity fitted by quasi second order–kinetic equation was 128.53 mg/g and 81.9 mg/g, which are very similar to the experimental equilibrium adsorption capacity, while the adsorption capacity was fitted by the quasi first–order kinetic equation which was quite different from the experimental value. This indicates that the adsorption process includes surface adsorption, external liquid film diffusion, and intra-particle diffusion. In addition, the intra–particle diffusion model was used to analyze the diffusion rate. As shown in Figure 8d, the adsorption process of Pb2+ in solution can be divided into two parts. The first stage is the migration of Pb2+ from solution to sample surface, which is controlled by molecular diffusion and membrane diffusion. The second stage may be caused by diffusion within the particles. The result indicates that membrane diffusion and intra–particle diffusion may occur simultaneously, which are the main rate control steps in the whole adsorption process.

3.6. Adsorption/Desorption Recycling Experiment

The adsorption/desorption cycle experiment can not only understand the utilization efficiency of the sample but can also help to further understand the adsorption mechanism. The adsorption/desorption capacity of sample for each number of experiments is shown in Figure 9. The adsorption capacity of UV/BGO700, bake/BGO700, UV/B700, and bake/B700 gradually decreased and finally reached a stable trend, which was 57.4 mg/g, 50.6 mg/g, 26.6 mg/g, and 28.5 mg/g, respectively. The reason for the low adsorption efficiency may be that the Pb2+, through the action of hydrogen bonds, π–π bonds, and oxygen–containing functional groups, was not completely desorbed by HCl. In addition, the adsorption capacity of the sample irradiated by ultraviolet light is slightly higher than that of the vacuum-dried sample, which is due to the introduction of some oxygen-containing functional groups on the surface of the sample by ultraviolet light irradiation [50]. According to Figure 9b, the desorption efficiency ranged from 54.3 to 63.3%.

3.7. Adsorption Mechanism of Pb2+ on Biochar-Supported Graphene Oxide

3.7.1. SEM-EDS Analysis after Adsorption

The structural characteristics (specific surface area, micropores) and surface properties (oxygen-containing functional groups) affect the adsorption process. At the same time, various interactions such as surface precipitation, electrostatic interaction, and π–π interaction are also conducive to the adsorption of Pb2+. The microscopic morphology and surface composition analysis after adsorption are shown in Figure 10. There are many small spots on the surface and pores of BGO700. According to surface composition analysis, these small spots are adsorbed Pb2+. This is because there are abundant oxygen-containing functional groups with negative charge on the surface of BGO700, which can be complex with Pb2+. At the same time, the more negative surface charges, the stronger the electrostatic interaction between Pb2+ and BGO700.

3.7.2. XPS Analysis before and after Adsorption

The XPS spectra of total elements and peaks before and after adsorption of Pb2+ by BGO700 is shown in Figure 11. It can be seen from Figure 11 that the Pb peaks were added in the total spectrum after adsorption, which are pb4f5/2 and pb4f7/2, respectively. The corresponding binding energies are 144.23 EV and 139.4 EV detected by XPS. According to the chemical state database of XPS, Pb exists in the form of Pb2+ on the surface of BGO700. The peaks of C1s, O1s, and N1s did not shift significantly before and after adsorption, which indicates that the chemical valence states of C, O, N elements and outer electron density did not change basically before and after adsorption of Pb2+ by BGO700.

3.7.3. FTIR and XRD Analysis before and after Adsorption

The analysis of FTIR and XRD before and after the adsorption of Pb2+ by BGO700 is shown in Figure 12. The oxygen–containing functional groups of BGO700 before the wavelength of 1700 cm−1 have changed after adsorption of Pb2+, in which the wavelength of C-O group has shifted from 1395.5 cm−1 to the left by 26.9 cm−1, and the absorption peak of C=O group has shifted from 1126.9 cm−1 to 1162.3 cm−1. The strength of the stretching vibration was significantly weakened, which indicates that the oxygen–containing functional groups participate in the adsorption process, and the surface complexation with Pb2+ occurred on the surface of BGO700.
The XRD diffraction pattern of the sample before and after adsorption is shown in Figure 12b. The broad diffraction peak appears near 2θ = 23.06°, which corresponds to the graphitic carbon structure. The broadened diffraction peaks indicates that the carbon structure is amorphous [51]. This indicates that the graphite crystal structure of the biomass was destroyed during the pyrolysis process, so the biochar composite material has an amorphous structure. In addition, a relatively weak diffraction peak appears at 2θ = 43.18°, which is a characteristic peak of the carbon structure. It is formed by the cracking of the graphene structure during the pyrolysis process [52]. The strong and narrow diffraction peaks were observed in the vicinity of 2θ = 24.32°, 26.82°, 33.66°, 39.86°, 42.36°, which indicates that the adsorption of Pb2+ exists as a crystal structure, and that biochar composites could adsorb Pb2+ by precipitation.

4. Conclusions

Graphene oxide significantly improves the ability of biochar to adsorb heavy metal ions. The pH of the solution has the greatest impact on the adsorption effect of biochar composite material. When the pH of the solution was less than 4, it inhibited the adsorption of Pb2+ by BGO. The rotational speed and the initial concentration of Pb2+ mainly affect the adsorption performance by affecting the contact rate between Pb2+ and the adsorbent surface. The response surface method optimization design experiment results showed that the order of influence on the adsorption of Pb2+ was solution pH > adsorbent dosage > rotating speed. The optimum conditions for the maximum Pb2+ adsorption capacity (125.17 mg/g) were as follows: solution pH was 4.97, rotation speed was 172.97 rpm, adsorbent dosage was 0.086 g. In the adsorption–desorption experiment, the desorption efficiency ranged from 54.3 to 63.3%. The process of Pb2+ adsorption by BGO was spontaneous and endothermic, mainly through electrostatic interaction and surface complexation. It was a heterogeneous adsorption process with heterogeneous surface, including surface adsorption, external liquid film diffusion, and intra-particle diffusion.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijerph19084790/s1, Table S1. The adsorbent information; Table S2. Experimental design and results; Table S3. Significance test and variance analysis for the quadratic regression model.

Author Contributions

All authors made contributions to this article. Z.L. and S.Y. conceived and designed the experiments; L.Z. and S.Y. performed the sampling and experiments; S.T. and Y.L. analyzed the data; S.Y. wrote the original draft; Z.L. and J.Z. revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the National Natural Science Foundation of China (No. 51464014; No. 52160019).

Institutional Review Board Statement

The studies not involving humans or animals.

Informed Consent Statement

The studies not involving humans.

Data Availability Statement

The study did not report any data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, Y.D.; Ouyang, W.; Wang, A.H.; Liu, L.H.; Lin, C.Y.; He, M.C. Synergetic Loss of Heavy Metal and Phosphorus: Evidence from Geochemical Fraction and Estuary Sedimentation. J. Hazard. Mater. 2021, 416, 125710. [Google Scholar] [CrossRef]
  2. Jia, Z.Z.; Li, S.Y.; Liu, Q.X.; Jiang, F.; Hu, J.T. Distribution and Partitioning of Heavy Metals in Water and Sediments of a Typical Estuary (Modaomen, South China): The Effect of Water Density Stratification Associated with Salinity. Environ. Pollut. 2021, 287, 117277. [Google Scholar] [CrossRef] [PubMed]
  3. Bezzina, J.P.; Ruder, L.D.; Dawson, R.; Ogden, M.D. Ion Exchange Removal of Cu(II), Fe(II), Pb(II) and Zn(II) from Acid Extracted Sewage Sludge—Resin Screening in Weak Acid Media. Water Res. 2019, 158, 257–267. [Google Scholar] [CrossRef] [PubMed]
  4. Su, Q.; Dai, H.C.; Chen, H.; Lin, Y.; Xie, Y.; Karthikeyan, R. General Equilibrium Analysis of the Cobenefits and Trade-Offs of Carbon Mitigation on Local Industrial Water Use and Pollutants Discharge in China. Environ. Sci. Technol. 2019, 53, 1715–1724. [Google Scholar] [CrossRef] [PubMed]
  5. Pateda, S.M.; Sakakibara, M.; Sera, K. Element Rich Area Associated with Human Health Disorders: A Geomedical Science Approach to Potentially Toxic Elements Contamination. Int. J. Environ. Res. Public Health 2021, 18, 12202. [Google Scholar] [CrossRef]
  6. Liu, Z.W.; Lu, C.B.; Yang, S.; Zeng, J.F.; Yin, S.Y. Release Characteristics of Manganese in Soil under Ion-absorbed Rare Earth Mining Conditions. Soil Sediment Contam. 2020, 29, 703–720. [Google Scholar] [CrossRef]
  7. Sharma, S.; Nagpal, A.K.; Kaur, I. Appraisal of Heavy Metal Contents in Groundwater and Associated Health Hazards Posed to Human Population of Ropar Wetland, Punjab, India and its Environs. Chemosphere 2019, 227, 179–190. [Google Scholar] [CrossRef]
  8. Chen, Y.G.; He, X.L.S.; Huang, J.H.; Luo, R.; Ge, H.Z.; Wołowicz, A.; Wawrzkiewicz, M.; Gładysz-Płaska, A.; Li, B.; Yu, Q.X.; et al. Impacts of Heavy Metals and Medicinal Crops on Ecological Systems, Environmental Pollution, Cultivation, and Production Processes in China. Ecotoxicol. Environ. Saf. 2021, 219, 112336. [Google Scholar] [CrossRef]
  9. Sarma, G.K.; Sen Gupta, S.; Bhattacharyya, K.G. Nanomaterials As Versatile Adsorbents for Heavy Metal Ions in Water: A review. Environ. Sci. Pollut. Res. 2019, 26, 6245–6278. [Google Scholar] [CrossRef] [PubMed]
  10. Politaeva, N.A.; Smyatskaya, Y.A.; Tatarintseva, E.A. Using Adsorption Material Based on the Residual Biomass of Chlorella Sorokiniana Microalgae for Wastewater Purification to Remove Heavy Metal Ions. Chem. Petrol. Eng. 2020, 55, 907–912. [Google Scholar] [CrossRef]
  11. Qiao, L.Z.; Li, S.S.; Li, Y.L.; Liu, Y.; Du, K.F. Fabrication of Superporous Cellulose Beads Via Enhanced Inner Cross-linked Linkages for High Efficient Adsorption of Heavy Metal Ions. J. Clean. Prod. 2020, 253, 120017. [Google Scholar] [CrossRef]
  12. Huggins, T.M.; Haeger, A.; Biffinger, J.C.; Ren, Z.J. Granular Biochar Compared with Activated Carbon for Wastewater Treatment and Resource Recovery. Water Res. 2016, 94, 225–232. [Google Scholar] [CrossRef] [Green Version]
  13. Meng, Y.; Xiao, L.; Muslim, A.; Hojiahmat, M. Improving the Adsorption of Poly(o-phenylenediamine) to Heavy Metal Ions in Aqueous Solution Through its Composite with Carbon Dots. J. Polym. Res. 2021, 28, 404. [Google Scholar] [CrossRef]
  14. Laysandra, L.; Ondang, I.J.; Ju, Y.H.; Putro, J.N.; Santoso, S.P.; Soetarejo, F.E.; Ismadji, S. An environment-friendly composite as an adsorbent for removal Cu (II) ions. Environ. Sci. Pollut. Res. 2019, 26, 22979–22989. [Google Scholar] [CrossRef] [PubMed]
  15. Pal, D.B.; Singh, A.; Jha, J.M.; Srivastava, N.; Hashem, A.; Alakeel, M.A.; Allah, E.F.A.; Gupta, V.K. Low-cost Biochar Adsorbents Prepared from Date and Delonix Regia Seeds for Heavy Metal Sorption. Bioresour. Technol. 2021, 339, 125606. [Google Scholar] [CrossRef] [PubMed]
  16. Odinga, E.S.; Waigi, M.G.; Gudda, F.O.; Wang, J.; Yang, B.; Hu, X.J.; Li, S.Y.; Gao, Y.Z. Occurrence, Formation, Environmental Fate and Risks of Environmentally Persistent Free Radicals in Biochars. Environ. Int. 2019, 134, 105172. [Google Scholar] [CrossRef] [PubMed]
  17. Shi, D.; Xie, C.; Wang, J.; Xiong, L. Changes in the Structures and Directions of Heavy Metal-Contaminated Soil Remediation Research from 1999 to 2020: A Bibliometric & Scientometric Study. Int. J. Environ. Res. Public Health 2021, 18, 7358. [Google Scholar] [CrossRef]
  18. Alam, M.S.; Gorman-Lewis, D.; Chen, N.; Flynn, S.L.; Yong, S.O.; Konhauser, K.O.; Alessi, D.S. Thermodynamic Analysis of Nickel(II) and Zinc(II) Adsorption to Biochar. Environ. Sci. Technol. 2018, 52, 6246–6255. [Google Scholar] [CrossRef]
  19. Wei, S.Y.; Zhu, M.B.; Fan, X.J.; Song, J.Z.; Peng, P.A.; Li, K.M.; Jia, W.L.; Song, H.Y. Influence of Pyrolysis Temperature and Feedstock on Carbon Fractions of Biochar Produced from Pyrolysis of Rice Straw, Pine Wood, Pig Manure and Sewage Sludge. Chemosphere 2019, 218, 624–631. [Google Scholar] [CrossRef] [PubMed]
  20. Yuan, H.R.; Lu, T.; Huang, H.Y.; Zhao, D.D.; Kobayashi, N.; Chen, Y. Influence of Pyrolysis Temperature on Physical and Chemical Properties of Biochar Made from Sewage Sludge. J. Anal. Appl. Pyrolysis 2015, 112, 284–289. [Google Scholar] [CrossRef]
  21. Panwar, N.L.; Pawar, A. Influence of Activation Conditions On the Physicochemical Properties of Activated Biochar: A Review. Biomass Conv. Bioref. 2020, 12, 925–947. [Google Scholar] [CrossRef]
  22. Wang, L.; Wang, Y.J.; Ma, F.; Tankpa, V.; Bai, S.S.; Guo, X.M.; Wang, X. Mechanisms and Reutilization of Modified Biochar Used for Removal of Heavy Metals from Wastewater: A review. Sci. Total Environ. 2019, 668, 1298–1309. [Google Scholar] [CrossRef] [PubMed]
  23. Intan, N.N.; Pfaendtner, J. The Composition of Oxygen Functionals Groups at the Surface of Carbon-Based Graphitic Anode. Phys. Chem. 2021. [Google Scholar] [CrossRef]
  24. Intan, N.N.; Pfaendtner, J. Effect of Fluoroethylene Carbonate Additives on the Initial Formation of the Solid Electrolyte Interphase on an Oxygen-Functionalized Graphitic Anode in Lithium-Ion Batteries. ACS Appl. Mater. Interfaces 2021, 13, 8169–8180. [Google Scholar] [CrossRef]
  25. Rustam, S.; Intan, N.N.; Pfaendtner, J. Effect of Graphitic Anode Surface Functionalization on the Structure and Dynamics of Electrolytes at the Interface. J. Chem. Phys. 2021, 155, 134702. [Google Scholar] [CrossRef] [PubMed]
  26. Jung, K.W.; Lee, S.Y.; Choi, J.W.; Lee, Y.J. A Facile One-pot Hydrothermal Synthesis of Sydroxyapatite/Biochar Nanocomposites: Adsorption Behavior and Mechanisms for the Removal of Copper(II) from Aqueous Media. Chem. Eng. J. 2019, 369, 529–541. [Google Scholar] [CrossRef]
  27. Inyang, M.I.; Gao, B.; Yao, Y.; Xue, Y.W.; Zimmerman, A.; Mosa, A.; Pullammanappallil, P.; Yong, S.O.; Cao, X.D. A Review of Biochar as A Low-cost Adsorbent for Aqueous Heavy Metal Removal. Crit. Rev. Environ. Sci. Technol. 2016, 46, 406–433. [Google Scholar] [CrossRef]
  28. Zhang, Y.; Cao, B.; Zhao, L.L.; Sun, L.L.; Gao, Y.; Li, J.J.; Yang, F. Biochar-supported Reduced Graphene Oxide Composite for Adsorption and Coadsorption of Atrazine and Lead Ions. Appl. Surf. Sci. 2017, 427, 147–155. [Google Scholar] [CrossRef]
  29. Zhao, G.X.; Li, J.X.; Ren, X.M.; Chen, C.L.; Wang, X.K. Few-Layered Graphene Oxide Nanosheets As Superior Sorbents for Heavy Metal Ion Pollution Management. Environ. Sci. Technol. 2011, 45, 10454–10462. [Google Scholar] [CrossRef] [PubMed]
  30. Huang, D.L.; Wang, X.; Zhang, C.Z.; Zeng, G.M.; Peng, Z.W.; Zhou, J.C.; Wang, M.; Wang, R.Z.; Hu, Z.X.; Qin, X. Sorptive Removal of Ionizable Antibiotic Sulfamethazine from Aqueous Solution by Graphene Oxide-coated Biochar Nanocomposites: Influencing Factors and Mechanism. Chemosphere 2017, 186, 414–421. [Google Scholar] [CrossRef]
  31. Chen, Y.Q.; Chen, L.B.; Bai, H.; Li, L. Graphene Oxide–chitosan Composite Hydrogels as Broad-spectrum Adsorbents for Water Purification. J. Mater. Chem. A 2013, 1, 1992–2001. [Google Scholar] [CrossRef]
  32. Huang, X.F.; Deng, Q.L.; Wang, X.Z.; Deng, H.Q.; Zhang, T.H.; Liao, H.W.; Jiang, J.L. High-efficient Removal of U(VI) from Aqueous Solution by Self-assembly Pomelo Peel/palygorskite Composite. Environ. Sci. Pollut. Res. 2021, 28, 17290–17305. [Google Scholar] [CrossRef]
  33. Hummers, W.S.; Offeman, R.E. Preparation of Graphitic Oxide. J. Am. Chem. Soc. 1958, 80, 1334–1339. [Google Scholar] [CrossRef]
  34. Zhou, X.; Liu, Z.P. A Scalable, Solution-Phase Processing Route to Graphene Oxide and Graphene Ultralarge Sheets. Chem. Comm. 2010, 46, 2611–2613. [Google Scholar] [CrossRef]
  35. Anyika, C.; Asri, N.A.M.; Majid, Z.A.; Jaafar, J.; Yahya, A. Batch Sorption–desorption of As(III) from Waste Water by Magnetic Palm Kernel Shell Activated Carbon Using Optimized Box–Behnken Design. Appl. Water Sci. 2017, 7, 4573–4591. [Google Scholar] [CrossRef]
  36. Mohammadnia, E.; Hadavifar, M.; Veisi, H. Kinetics and Thermodynamics of Mercury Adsorption onto Thiolated Graphene Oxide Nanoparticles. Polyhedron 2019, 173, 114139. [Google Scholar] [CrossRef]
  37. Nigam, M.; Rajoriya, S.; Singh, S.R.; Kumar, P. Adsorption of Cr (VI) ion from tannery wastewater on tea waste: Kinetics, equilibrium and thermodynamics studies. J. Environ. Chem. Eng. 2019, 7, 103188. [Google Scholar] [CrossRef]
  38. Peng, W.J.; Li, H.Q.; Liu, Y.Y.; Song, S.X. A Review on Heavy Metal Ions Adsorption from Water by Graphene Oxide and Its Composites. J. Mol. Liq. 2017, 230, 496–504. [Google Scholar] [CrossRef]
  39. Keiluweit, M.; Nico, P.S.; Johnson, M.G.; Kleber, M. Dynamic Molecular Structure of Plant Biomass-derived Black Carbon (Biochar). Environ. Sci. Technol. 2010, 44, 1247–1253. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Xue, Y.J.; Wang, C.; Hu, Z.H.; Zhou, Y.; Xiao, Y.; Wang, T. Pyrolysis of Sewage Sludge by Electromagnetic Induction: Biochar Properties and Application in Adsorption Removal of Pb(II), Cd(II) from Aqueous Solution. Waste Manag. 2019, 89, 48–56. [Google Scholar] [CrossRef]
  41. Wu, Z.Y.; Chen, X.X.; Yuan, B.L.; Fu, M.L. A Facile Foaming-polymerization Strategy to Prepare 3D MnO2 Modified Biochar-based Porous Hydrogels for Efficient Removal of Cd(II) and Pb(II). Chemosphere 2020, 239, 124745. [Google Scholar] [CrossRef]
  42. Wu, C.X.; Li, Y.G.; Chen, M.J.; Luo, X.; Chen, Y.W.; Belzile, N.; Huang, S. Adsorption of Cadmium on Degraded Soils Amended with Maize-Stalk-Derived Biochar. Int. J. Environ. Res. Public Health 2018, 15, 2331. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Singh, E.; Kumara, A.; Mishr, R.; You, S.; Singh, L.; Kumar, S.; Kumar, S. Pyrolysis of Waste Biomass and Plastics for Production of Biochar and Its Use for Removal of Heavy Metals from Aqueous Solution. Bioresour. Technol. 2021, 320, 124278. [Google Scholar] [CrossRef] [PubMed]
  44. Zhou, S.Y.; Xue, A.L.; Zhao, Y.J.; Wang, Q.W.; Yan, C.; Li, M.Y.; Xing, W.H. Competitive Adsorption of Hg2+, Pb2+ and Co2+ Ions on Polyacrylamide/attapulgite. Desalination 2010, 270, 269–274. [Google Scholar] [CrossRef]
  45. Deng, X.Y.; Wang, L.X.; Xiu, Q.H.; Wang, Y.; Han, H.; Dai, D.M.; Xu, Y.J.; Gao, H.T.; Liu, X. Adsorption Performance and Physicochemical Mechanism of MnO2-polyethylenimine-tannic Acid Composites for the Removal of Cu(II) and Cr(VI) from Aqueous Solution. Front. Chem. Sci. Eng. 2021, 15, 538–551. [Google Scholar] [CrossRef]
  46. Wang, Z.; Yao, M.J.; Wang, X.R.; Li, S.P.; Liu, Y.; Yang, G.H. Influence of Reaction Media on Synthesis of Dialdehyde Cellulose/GO Composites and Their Adsorption Performances on Heavy Metals. Carbohydr. Polym. 2020, 232, 115781. [Google Scholar] [CrossRef] [PubMed]
  47. Yetilmezsoy, K.; Demirel, S.; Vanderbei, R.J. Response Surface Modeling of Pb(II) Removal from Aqueous Solution by Pistacia Vera, L.: Box–Behnken Experimental Design. J. Hazard. Mater. 2009, 171, 551–562. [Google Scholar] [CrossRef] [PubMed]
  48. Sivamani, S.; Prasad, B.S.N.; Nithya, K.; Hosseini-Bandegharaei, A. Back-propagation Neural Network: Box–Behnken Design Modelling for Optimization of Copper Adsorption on Orange Zest Biochar. Int. J. Environ. Sci. Technol. 2021, 1–16. [Google Scholar] [CrossRef]
  49. Barati, A.; Moghadam, E.A.; Miri, T.; Asgari, M. Rapid Removal of Heavy Metal Cations by Novel Nanocomposite Hydrogels Based on Wheat Bran and Clinoptilolite: Kinetics, Thermodynamics, and Isotherms. Water Air Soil Pollut. 2014, 225, 2096. [Google Scholar] [CrossRef]
  50. Wojnárovits, L.; Földváry, C.M.; Takács, E. Radiation-induced Grafting of Cellulose for Adsorption of Hazardous Water Pollutants: A Review. Radiat. Phys. Chem. 2010, 79, 848–862. [Google Scholar] [CrossRef]
  51. Nikkhah, A.A.; Zilouei, H.; Asadinezhad, A.; Keshavarz, A. Removal of Oil from Water Using Polyurethane Foam Modified with Nanoclay. Chem. Eng. J. 2015, 262, 278–285. [Google Scholar] [CrossRef]
  52. Misnon, I.I.; Zain, N.K.M.; Aziz, R.A.; Vidyadharan, B.; Jose, R. Electrochemical Properties of Carbon from Oil Palm Kernel Shell for High Performance Supercapacitors. Electrochim. Acta 2015, 174, 78–86. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a1) surface micromorphology of B700, (a2) component spectra of B700; (b1) surface micromorphology of GO700, (b2) component spectra of GO700; (c1) surface micromorphology of BGO700, (c2) component spectra of BGO700. Yellow square in the images indicates the zone where the EDS spectra was performed.
Figure 1. (a1) surface micromorphology of B700, (a2) component spectra of B700; (b1) surface micromorphology of GO700, (b2) component spectra of GO700; (c1) surface micromorphology of BGO700, (c2) component spectra of BGO700. Yellow square in the images indicates the zone where the EDS spectra was performed.
Ijerph 19 04790 g001
Figure 2. FTIR spectra of the obtained materials: B = biochar, GO = graphene oxide, BGO = composite biochar-graphene oxide.
Figure 2. FTIR spectra of the obtained materials: B = biochar, GO = graphene oxide, BGO = composite biochar-graphene oxide.
Ijerph 19 04790 g002
Figure 3. Adsorption capacity of biochar (B), composite biochar-graphene oxide (BGO), and graphene oxide (GO) obtained at different pyrolyzed temperatures.
Figure 3. Adsorption capacity of biochar (B), composite biochar-graphene oxide (BGO), and graphene oxide (GO) obtained at different pyrolyzed temperatures.
Ijerph 19 04790 g003
Figure 4. (a) The effect of solution pH; (b) The morphological distribution of Pb2+ in aqueous solution; (c) The effect of rotating speed; (d) The effect of adsorbent dosage.
Figure 4. (a) The effect of solution pH; (b) The morphological distribution of Pb2+ in aqueous solution; (c) The effect of rotating speed; (d) The effect of adsorbent dosage.
Ijerph 19 04790 g004
Figure 5. The three–dimensional response surfaces of interaction effect; (a) Coupling of rotational speed and pH; (b) Coupling of dosage and pH; (c) Coupling of dosage and rotational speed; (d) Actual value and predicted value.
Figure 5. The three–dimensional response surfaces of interaction effect; (a) Coupling of rotational speed and pH; (b) Coupling of dosage and pH; (c) Coupling of dosage and rotational speed; (d) Actual value and predicted value.
Ijerph 19 04790 g005
Figure 6. (a) The adsorption isotherms of Pb2+ by BGO700; (b) The adsorption isotherm model of Langmuir; (c) The adsorption isotherm model of Freundlich; (d) The adsorption isotherm model of Temkin.
Figure 6. (a) The adsorption isotherms of Pb2+ by BGO700; (b) The adsorption isotherm model of Langmuir; (c) The adsorption isotherm model of Freundlich; (d) The adsorption isotherm model of Temkin.
Ijerph 19 04790 g006
Figure 7. The thermodynamic parameter fitting.
Figure 7. The thermodynamic parameter fitting.
Ijerph 19 04790 g007
Figure 8. The kinetic fitting curve; (a) The variation of Pb2+ adsorption capacity with oscillation time; (b) The quasi first–order kinetics; (c) The quasi second–order kinetics; (d) The intra–particle diffusion model.
Figure 8. The kinetic fitting curve; (a) The variation of Pb2+ adsorption capacity with oscillation time; (b) The quasi first–order kinetics; (c) The quasi second–order kinetics; (d) The intra–particle diffusion model.
Ijerph 19 04790 g008
Figure 9. (a) The adsorption capacity of cycle experiments; (b) The desorption capacity of cycle experiments.
Figure 9. (a) The adsorption capacity of cycle experiments; (b) The desorption capacity of cycle experiments.
Ijerph 19 04790 g009
Figure 10. The SEM-EDS images after adsorption.
Figure 10. The SEM-EDS images after adsorption.
Ijerph 19 04790 g010
Figure 11. The XPS analysis spectrum before and after adsorption: (a) total element, (b) C1s peak split, (c) O1s peak split, and (d) N1s peak split.
Figure 11. The XPS analysis spectrum before and after adsorption: (a) total element, (b) C1s peak split, (c) O1s peak split, and (d) N1s peak split.
Ijerph 19 04790 g011
Figure 12. (a) The FTIR of biochar composites before and after adsorption; (b) The XRD of biochar composites before and after adsorption.
Figure 12. (a) The FTIR of biochar composites before and after adsorption; (b) The XRD of biochar composites before and after adsorption.
Ijerph 19 04790 g012
Table 1. Design scheme of single-factor adsorption experiment.
Table 1. Design scheme of single-factor adsorption experiment.
NumberFactorsVariation Range
ASolution pH2, 3, 4, 5, 6
BRotating speed (rpm)100, 125, 150, 175, 200
CAdsorbent dosage (g)0.005, 0.01, 0.02, 0.05, 0.1
Table 2. The physical properties and element analysis of samples.
Table 2. The physical properties and element analysis of samples.
SamplespHElement Content (%)Atomic RatioSSA
(m2/g)
PV
(cm3/g)
PD
(nm)
Yield
(%)
Ash
(%)
CHNSO aH/CO/C(O + N)/C
B3008.8859.613.421.850.7623.660.0570.400.431.820.01221.75343.810.70
B5009.6867.272.261.800.6217.750.0610.240.275.390.01511.77333.510.30
B7009.9874.091.771.010.3812.400.0240.170.1818.500.01661.85126.610.35
BGO3009.4356.633.282.020.4924.390.0580.430.478.670.01551.78945.313.19
BGO5009.6459.422.281.750.6222.450.0380.380.4123.860.01731.99129.413.48
BGO70010.168.231.371.380.7414.940.0200.220.2438.230.01892.11427.411.34
“a” means that ω (O) was calculated by difference method; ω (O) = 100% − ω (C) − ω (H) − ω (N) − ω (S) − ω (ash). “SSA” means that specific surface area; “PV” means that pore volume; “PD” means that pore diameter.
Table 3. The fitting parameters of isothermal adsorption model.
Table 3. The fitting parameters of isothermal adsorption model.
ModelParameterTemperature (°C)
25 °C35 °C45 °C
q e , exp   =   130.21 q e , exp   =   135.88 q e , exp   =   148.25
Langmuir
q e = K L q max C e 1 + K L C e
q max (mg/g)127.66118.04123.28
K L (L/mg)0.1460.2830.281
R 2 0.8380.8480.744
RSS 65.1915.6762.82
Freundlich
q e = K F C e 1 / n
1 / n 0.1080.0790.091
K F (L/mg)66.1776.2775.91
R 2 0.9700.9450.894
RSS 12.155.72626.07
Temkin
q e = B + AlnC e
A 57.4972.2071.12
B 11.488.2689.709
R 2 0.9590.9350.875
RSS 16.476.68930.63
q e , exp is the actual experimental value (mg/g); q max is the theoretical maximum adsorption (mg/g); K L , 1 / n , K F , A, B are adsorption isotherm model equation parameters; R2 is the correlation coefficient; RSS is the sum of squares of residuals.
Table 4. The thermodynamic parameters of Pb2+ adsorption by BGO700.
Table 4. The thermodynamic parameters of Pb2+ adsorption by BGO700.
T(K)ΔG (KJ/mol)ΔH (KJ/mol)ΔS (KJ/(mol·K)
298−3.93413.350.058
308−4.514
318−5.094
Table 5. The fitting parameters of adsorption kinetics model.
Table 5. The fitting parameters of adsorption kinetics model.
Kinetic ModelQuasi First Order KineticsQuasi Second Order Kinetics
ln ( q e q t )   =   l n q e k 1 t t q t   =   1 k 2 q e 2 + 1 q e t
Parameter q e , exp q e , cal k 1 R 2 RSS q e , cal k 2 ( × 10 5 ) R 2 RSS
BGO700124.4248.210.001270.9740.028128.536.730.9990.026
B70079.137.8260.001440.9490.03681.99.540.9980.205
Kinetic ModelIntra-Particle Diffusion Model
q t   =   K P t 0.5 + C
Parameter q e , exp K P 1 C 1 R 1 2 RSS 1 K P 2 C 2 R 2 2 RSS 2
BGO700124.422.93659.520.98612.680.70489.420.90234.44
B70079.11.45339.630.9831.8180.50853.850.9288.036
q e , exp is the actual experimental value (mg/g); q e , cal is the theoretical value (mg/g); k 1 ,   k 2 ,   K P 1 , 2 , C 1 , 2 are kinetic equation fitting parameters; R2 is the correlation coefficient; RSS is the sum of squares of residuals.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, Z.; Yang, S.; Zhang, L.; Zeng, J.; Tian, S.; Lin, Y. The Removal of Pb2+ from Aqueous Solution by Using Navel Orange Peel Biochar Supported Graphene Oxide: Characteristics, Response Surface Methodology, and Mechanism. Int. J. Environ. Res. Public Health 2022, 19, 4790. https://doi.org/10.3390/ijerph19084790

AMA Style

Liu Z, Yang S, Zhang L, Zeng J, Tian S, Lin Y. The Removal of Pb2+ from Aqueous Solution by Using Navel Orange Peel Biochar Supported Graphene Oxide: Characteristics, Response Surface Methodology, and Mechanism. International Journal of Environmental Research and Public Health. 2022; 19(8):4790. https://doi.org/10.3390/ijerph19084790

Chicago/Turabian Style

Liu, Zuwen, Shi Yang, Linan Zhang, Jinfeng Zeng, Shuai Tian, and Yuan Lin. 2022. "The Removal of Pb2+ from Aqueous Solution by Using Navel Orange Peel Biochar Supported Graphene Oxide: Characteristics, Response Surface Methodology, and Mechanism" International Journal of Environmental Research and Public Health 19, no. 8: 4790. https://doi.org/10.3390/ijerph19084790

APA Style

Liu, Z., Yang, S., Zhang, L., Zeng, J., Tian, S., & Lin, Y. (2022). The Removal of Pb2+ from Aqueous Solution by Using Navel Orange Peel Biochar Supported Graphene Oxide: Characteristics, Response Surface Methodology, and Mechanism. International Journal of Environmental Research and Public Health, 19(8), 4790. https://doi.org/10.3390/ijerph19084790

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop