Blog

Analysis of influencing factors of phenanthrene adsorption by different soils in Guanzhong basin based on response surface method | Scientific Reports

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Scientific Reports volume  12, Article number: 20906 (2022 ) Cite this article Falcon Tube

Analysis of influencing factors of phenanthrene adsorption by different soils in Guanzhong basin based on response surface method | Scientific Reports

Adsorption desorption is an important behavior affecting the migration of phenanthrene in soil. In this study, three typical soils of loess, silts and silty sand in Guanzhong Basin, Shaanxi Province, China were used as adsorbents. Batch equilibrium experiments were carried out to study the adsorption desorption kinetics and isotherm of phenanthrene in different soils. Response surface method (RSM) was used to study the effects of temperature, pH, phenanthrene concentration and organic matter content on soil adsorption of phenanthrene. The results showed that after adsorption, the outline of soil particles became more blurred and the degree of cementation increased. The kinetic adsorption of phenanthrene by soil conforms to the quasi second-order kinetic model, and the adsorption desorption isotherm is nonlinear and conforms to the Freundlich model. Due to the difference of soil properties, the adsorption amount of phenanthrene by soil is loess > silty sand > silts. The thermodynamic results show that the adsorption of phenanthrene by soil is spontaneous and endothermic, and the desorption is spontaneous and exothermic. Through RSM, the interaction between phenanthrene concentration and soil organic matter in Loess and silts is significant, and the interaction between temperature and soil organic matter in silty sand is significant. Among the four factors affecting the adsorption rate of loess, silts and silty sand, soil organic matter is the most significant. The theoretical optimum adsorption rates of loess, silts and silty sand are 98.89%, 96.59% and 93.37% respectively.

Polycyclic aromatic hydrocarbons (PAHs) are a kind of toxic and harmful organic pollutants widely distributed in the natural environment where the soil bears more than 90% PAHs in the environment1,2,3,4,5. A series of transport and transformation behaviors will occur when PAHs enters the soil, including adsorption, desorption, biodegradation, photodegradation et al. Adsorption and desorption play an important role in influencing PAHs in soil environments6. PAHs entering the soil environment can be easily absorbed by soils. At the same time the absorbed PAHs can also be desorbed with the solvent. Absorbed PAHs in soils will enter into the food chain through the plants metabolism, posing a serious threat to the ecological environment and human health7,8. Understanding the transport and fate of PAHs in soil environment is paramount to determinations of the PAHs risks to human health and the environment. It is also to help providing strategies for soil contamination.

The adsorption and desorption mechanisms of PAHs in soils were mainly focused on the kinetics and isotherms theory9. The adsorption and desorption kinetics of PAHs was considered to be completed in two processes. One was the rapid adsorption process, which took several minutes to several hours to reach equilibrium10,13. The other was the slow adsorption process, which took months to several years14. Some research showed that adsorption kinetics of phenanthrene on soils conformed to the first-order kinetic model. Others are more consistent with the pseudo-second-order model15. In addition, the adsorption isotherm of PAHs on soils could be described quantitatively by fitting experimental data with different adsorption isotherms. The most commonly used isotherms were Freundlich16,17, Langmuir18, Temkin19 and Brunauer–Emmett–Teller (BET) models20.

The adsorption and desorption of PAHs in soils can be influenced by different factors such as PAHs properties, soil type, mineral composition, pH, soil porous structure, organic matter et al. The adsorption of PAHs on particulate was affected by the molecular weight and hydrophobicity of PAHs. PAHs with a greater octanol–water partition coefficient (Kow) tended to be highly affected by partitioning into organic matter, resulting in more difficult desorption from soil9. It is recommended that adsorption and desorption behavior of PAHs in soils were mainly based on mineral component and organic matter21,22. The desorption of phenanthrene decreased with the organic content increasing23. This was because organic matter in soil enhanced the sorption and affinity of PAHs due to partitioning into the soil organic phase rather than clay minerals surfaces24. Soil organic matter could be divided into dissolve phase and pore filling phase. The former was linear distribution model and the latter followed the adsorption isotherm model. Some studies showed that pH affected the adsorption of PAHs on soils and the adsorption capacity would be decreased with the increased pH25. In addition, temperature was another important factor in the adsorption and desorption behavior of PAHs. Zhu studied the adsorption characteristics of naphthalene and phenanthrene in soil at three different temperatures (25, 35, 45 °C). It was found that the adsorption rate of naphthalene and phenanthrene were negatively correlated with temperature. Ren26 suggested that temperature would affect the solubility of PAHs in water, thereby affect the adsorption strength of PAHs. Shi27 showed that the adsorption of naphthalene in soil was an endothermic reaction, which was driven by entropy increased. But in most cases, the adsorption of the PAHs was an exothermic reaction.

Guanzhong basin is located in Shaanxi province, China. It is the center of Shaanxi province and the main industrial and agricultural area. With the economic development and population growth, the municipal wastewater and industrial waste discharge were greatly increased, leading to the organic pollution risk of soil. However, little is known about the mechanisms in PAHs adsorption and desorption on soils in Guanzhong basin. In addition, the soils in Guanzhong basin have different physical and chemical properties28. The characteristics of soils could influence the adsorption and desorption behaviors of PAHs. The adsorption and desorption capacity of the soil can directly affect the downward migration of pollutants from the surface layer. If the pollutants are retained in the soil, thereby affecting the mobility of the pollutants, the process of continuing the infiltration of the pollutants is blocked, and it is not easy to infiltrate. It can reduce the pollution of the soil and groundwater. The effect of soil on pollutants is different due to different soil adsorption capacity. However, comparative studies on the adsorption and fixation ability of different soils in the same area are rare, and the main focus is on the clay layer and other rocks. Comparative studies of sexual block strength are relatively rare. In this paper, three typical soils in Guanzhong basin, named loess, silts and silty sand, were chosen, while the adsorption and desorption characteristics of phenanthrene on these soils were studied, including kinetics, isotherm and thermodynamics of adsorption and desorption. The kinetics and equilibrium adsorption model were used to describe the adsorption and desorption characteristics. In addition, the BBD experimental design response surface method was used to study the influence of phenanthrene concentration, soil organic matter content, temperature, pH and their interaction on soil adsorption of phenanthrene, and to determine the optimal adsorption rate and the influencing factors under this condition. To study the environmental behavior of PAHs in soil and provide the basis for pollution control, and analyze their adsorption mechanism through Fourier transform infrared spectroscopy and scanning electron microscopy.

Experimental reagents: phenanthrene (analytical pure, purity > 97%), purchased from Aladdin Industrial Company, methanol (chromatographically pure) from TEDA, calcium chloride, hydrochloric acid, sodium hydroxide and hydrogen peroxide were all analytically pure.

Experimental equipment: constant temperature water bath oscillator SHA-B, desktop centrifuge TGL-15B, high performance liquid chromatograph (Agilent 1260), constant temperature magnetic stirrer DF-101S, electronic balance LE204EI02, box-type resistance furnace SX2-4-10A.

The soil used in the experiment was selected from a vertical section of the Weihe River, and three different types of soils, namely loess plateau loess, first-order terrace silty sand and floodplain silts, were collected as the research object. The three typical soils were collected at a depth of 0–20 cm, and the sundries were removed. After air-drying and grinding, they were sieved with a 60-mesh sieve. The soil samples after sieving were stored in a sealed sample bag for future use. The basic properties of the soils are summarized in Table 1.

Preparation of soil samples with different organic matter content: Weigh 10.00 g of each of the three soil samples into a conical flask, and burn them at 0 °C, 200 °C and 400 °C respectively. The organic matter content was determined according to GB/T642-1999 using potassium dichromate method. The specific operations are as follows:

Weigh 3 g soil samples respectively and put them into a dry conical flask, slowly drop 10 mL of potassium dichrochrome standard solution and shake well, add 1 mL of concentrated sulfuric acid, and put a small funnel into the test tube mouth to collect condensed liquid. Put the conical flask into the oven at about 185 °C, control the temperature of 170 ~ 180 °C to boil for five minutes, and take out a little cooler. Rinse the funnel and the inner wall of the conical flask with pure water, control the pure water within 10 mL, and then add 3 drops of Aphaline indicator. The ferrous sulfate standard solution was titrated until the solution changed from yellow to brownish red after green as the end point. The amount of ferrous sulfate standard solution was recorded and estimated to 0.05 mL. Then, the soil sample burned at 500 °C for 5 h was used to replace the sample and the above steps were repeated for blank test.

The formula for calculating organic matter is as follows:

where, Om for organic matter content (%); c(Fe2+) for ferrous sulfate standard solution concentration (mol/L); V(Fe2+) is a blank titration ferrous sulfate dosage (mL); V(Fe2+) as the sample (ferrous sulfate dosage (mL); 0.003 is the molar mass kg/mol of 1/4 ferrous sulfate standard solution concentration; 1.724 is the factor of converting organic carbon to organic matter. The organic matter content under the conditions is shown in Table 2.

Adsorption experiment was carried out using a batch procedure. An appropriate soil was transferred to centrifuge tube. 10 mg/L of Phenanthrene electrolyte (in background solution of 0.01 mol/L CaCl2 to maintain a constant ionic strength) was added to the tube. All tubes were immediately sealed and then mechanically shaken for 24 h in a thermostatic oscillation incubator (200 r/min) at 25 °C. In thermodynamic experiments, the temperatures were adjusted to 15 and 35 °C. Samples were taken at different times and then the suspensions were centrifuged at 8000 r/min for 3 min. A 2 mL supernatant was filtered through the membrane with a pore size of 0.22 μm and was then analyzed by HPLC. Each processing was set for three repetitions. A blank (no soil) was prepared for each time point. Phenanthrene loss during adsorption and filtration was negligible. To determine the contact time required to achieve sorption equilibrium, the suspensions were sampled at 0, 5, 10, 60, 120, 300, 480 and 1080 min29,30. The amount of Phenanthrene adsorbed was calculated by Eq. (2):

where, qe (mg/g) is the amount of phenanthrene absorbed; C0 (mg/L) and Ce (mg/L) are the initial and equilibrium aqueous concentrations, respectively; V(mL) is the solution volume; and m(g) is the mass of soil in the centrifuge tubes.

After adsorption experiments, desorption experiments were conducted by replacing supernatant with phenanthrene-free background solution after drying. The tubes were agitated (200 r/min) for 24 h and then centrifuged at 8000 r/min for 3 min. The supernatants were filtered through 0.22 μm syringe filters and analyzed by HPLC31.

Phenanthrene concentrations in the liquid phase were measured using a HPLC (Agilent 1260) equipped with a C18 analytical column (150 × 4.6 mm), and detected by UV spectrophotometry. The operating conditions were: mobile phase 90/10(V/V) of methanol/water, flow rate 1 mL/min, column temperature 35 °C, and detector wavelength 254 nm. The retention time for Phenanthrene was 8.1 min.

In order to better use the response surface method to optimize the experimental conditions and study the adsorption behavior of different soil samples to phenanthrene, 4 variables of phenanthrene concentration (A), pH (B), temperature (C) and soil organic matter content (D) were selected in this experiment, respectively, single factor experiments were carried out to obtain the subsequent response surface adsorption equilibrium time and the parameter range of the experimental design. According to the previous experimental results, the adsorption equilibrium state can be reached after 18 h of soil adsorption of phenanthrene. Therefore, 18 h was selected as the adsorption equilibrium time in the subsequent experiments. Taking soil organic matter content of 0 (coded value), initial phenanthrene concentration of 20 mg/L, temperature of 30 °C and pH = 6 as the center points of the response surface experimental design. At the same time, according to the Box-Behnken (BBD) experimental principle, the 4-factor 3-level response surface method was used to explore the main effect and interaction of the above-mentioned 4 factors on soil adsorption of phenanthrene. The coding value of the experimental design is coded by Eq. (3):

where xi is the coded value of the independent variable, XI is the actual value of the independent variable, and X0 is the actual value of the central point of the independent variable, ΔXI is the change value.

Through the data of various mathematical models (linear, two factor interaction, quadratic and cubic) and their variance analysis, it is shown that this method is most suitable for fitting with a quadratic polynomial model. The phenanthrene adsorption rate is selected as the response value, and its predicted value and independent variable are expressed by Eq. (4):

where R is the predicted response value of phenanthrene adsorption rate; β0 is the intercept coefficient; βi is a linear parameter; βii and βij is interaction and second-order coefficient respectively; xi and xj are coded independent variables.

The coded and experimental values of the experimental design are shown in Table 2.In order to make the data clear, the coded values are used for organic matter.

Pseudo-second-order kinetics (Eq. (5)) and the Elovich model (Eq. (6)) were used to fit (stimulate) the adsorption kinetics. The equations can be expressed as:

where, qt (mg/g) is the amount of phenanthrene adsorbed at time t(min) and under equilibrium conditions, respectively; and k2 (g/(mg⋅ min)) is the pseudo-second-order model rate constant for adsorption.

The Freundlich and Langmuir models were used to fit the adsorption and desorption isotherms. The two models were described dynamic equilibrium of adsorption, as the adsorption reaction proceeds. The adsorption site on the adsorbent surface was occupied constantly, making the adsorption rate gradually slow down. As the desorption rate was gradually accelerated in the process of adsorption, eventually the two rates were equal and the adsorption reaches the dynamic equilibrium. The Langmuir equation (Eq. (7)) is:

where, q(mg/g) is the amount of phenanthrene adsorbed by a mass of soil; and KL(L/mg) represents the Langmuir constant related to the bonding force of sorption.

Freundlich is an empirical equation, which is widely used in cases where adsorbent surface is not uniform and interaction between adsorbents exists. The Freundlich equation (Eq. (8)) is:

where, KF((mg/g)(mg/L) − n) is the adsorption equilibrium constant representing the adsorption capacity; \(n\) is a constant indicative of adsorption intensity.

The experiments were carried out at 293 K, 303 K, and 313 K. The thermodynamic parameters for the adsorption processes were obtained by the following equations (Eqs. (9), (10), (11)) 32:

where, T(K)is the absolute temperature; \(R\) is the gas constant (8.314 J/(mol⋅ K)); ΔS0 (kJ/(mol⋅ K)) is the standard entropy change; ΔG0 (kJ/mol) is the standard Gibbs free energy;and ΔH0 (kJ/mol) is the standard enthalpy change; ρw is the water density (mg⋅ L−1). ΔH0 and ΔS0 were obtained from the slopes and intercepts of a linear regression of lnK vs. \(1/T\) .

SEM images were shown in Fig. 1. It showed that the contour of three soils were fairly clear before adsorption. But it became fuzzier and the degree of cementation was increased when phenanthrene was adsorbed on the soils. According to the surface morphology, the silty sand (A) had furrows on the surface before adsorption compared with the fairly smooth without any furrows after adsorption (B). The silts (C) were flaky and the lamellar accumulation decreased (D). The loess (E) had a smooth surface with some flaky and rod like structure, after adsorption (F), the surface of loess increased in clay-like structure.

SEM micrographs of the three soil samples. (A) Silty sand; (B) Adsorbing 5 h of Silty sand; (C) Silts; (D) Adsorbing 5 h of Silts; (E) Loess; (F) Adsorbing 5 h of Loess.

Adsorption kinetics is one of the most important characteristics governing solute uptake rate and represents adsorption efficiency33. The sorption and desorption kinetics of phenanthrene in three soils were shown in Fig. 2. The results showed that the adsorption processes among all soils were similar. The kinetics of phenanthrene in soils was completed in two steps: a “fast” adsorption and a “slow” adsorption. The adsorption amount increased during 0-18h. It was a rapid reaction from 0 to 200 minutes. From 200 to 600 minutes, the adsorption amount increased slightly into balance. This phenomenon was due to the adsorption of phenanthrene occurred on the surface of soil organic matter. With the increase of time, soil surface adsorption sites were gradually saturated, causing the decrease of adsorption rate until reaching the equilibrium. Phenanthrene was a hydrophobic substance. It was easy to reach the soil surface and adhere to the grain surface. The results were consistent with the study of had also found that the balance time was approximately 18h and the adsorption amount increased with the adsorption reaction time34. Under the same conditions, loess had the highest adsorption capacity, which was mainly due to the highest organic content 18. The maximum phenanthrene sorption capacities ranked as follows: loess > silty sand > silts. As shown in Fig. 2, phenanthrene desorption in soils was relatively quick and its desorption equilibrium time was 3h. To reach an adequate desorption balance while remaining consistent with the adsorption reaction time, the balance time of the adsorption–desorption experiment was set at 18h. Generally, PAHs below 4 cycles could reach the adsorption equilibrium for about 16~24h.

(a)Adsorption equilibration curves of phenanthrene sorption in soils. (b) Desorption equilibration curves of phenanthrene sorption in soils.

Pseudo-second-order and Elovich models were used to study the phenanthrene adsorption mechanism (Table 3). Phenanthrene sorption kinetics were satisfactorily described by a pseudo-second-order model with coefficients of determination (R2) ranging from 0.99875 to 0.99847, compared with R2 values of 0.26508–0.73901 for the Elovich model. This well-fitting pseudo-second-order model indicated that the rate-limiting step was chemical adsorption, including electronic forces through sharing or exchange of electrons35,36. Moreover, it suggested that sorption was governed by the availability of sorption sites on the soil surfaces instead of by the phenanthrene concentration in solution.

The isotherm was used for quantitative analysis of phenanthrene transport from liquid to solid phase and for understanding the nature of interactions between phenanthrene and the soil matrix. The sorption and desorption isotherms of phenanthrene in soils were shown in Fig. 3. The data showed that phenanthrene adsorption and desorption capacities of three soils varied markedly due to their different physicochemical properties. With the increase of phenanthrene concentration, the adsorbed amount increased. At the same temperature, the adsorption capacity of silty sand was minimum while loess was maximum. This is mainly related to the soil physicochemical properties. At the same initial concentration, the temperature increase from 20 °C to 40 °C showed that the adsorption and desorption capacity decreased with temperature increase. On the one hand, the rise of temperature can increase the phenanthrene solubility in the liquid phase. On the other hand, it could reduce various forces between the soil surface and phenanthrene37.

(a)20 °C adsorption isotherms for phenanthrene in soils. (b)30 °C adsorption isotherms for phenanthrene in soils. (c)40 °C adsorption isotherms for phenanthrene in soils. (d) 20 °C desorption isotherms for phenanthrene in soils. (e) 30 °C desorption isotherms for phenanthrene in soils. (f) 40 °C desorption isotherms for phenanthrene in soils.

The Freundlich isotherm was used mainly for adsorption surfaces with nonuniform energy distribution, and the Langmuir isotherm was used for monolayer adsorption on perfectly smooth and homogeneous surfaces38. The experimental data were fitted with the Langmuir and Freundlich adsorption models, and the isotherm parameters logKF, 1/n, KL, qmax and the coefficient of determination (R2) of phenanthrene in soils were listed in Table 4.

As shown in Table 4, according to the coefficients of determination (R2), all soils were better fitted with the Freundlich model, which assumes that phenanthrene sorption and desorption occurs on a heterogeneous surface with the possibility of sorption being multi-layered39. This phenomenon has also been observed in humic acid and nanometer clay mineral40. It showed that the soil adsorption of organic matter was not only surface adsorption but also the process of soil organic matter distribution41,42,43 reached the equilibrium isotherm fitted well with the Freundlich equation when studying the adsorption behavior of aromatic compounds by solids.

To clarify the adsorption mechanisms, the thermodynamic parameters mentioned earlier were calculated and presented in Table 5. Generally, the value of Gibbs free energy changeΔG0 indicated the spontaneity of a chemical reaction. Therefore, it could evaluate whether sorption was relate to spontaneous interaction44. Negative values of ΔG0 indicated that the feasibility and spontaneous nature. The research was under the temperature range about 293–313 K. For adsorption process, all soils ΔG0 was < 0 (Fig. 4). It indicated that the processes of adsorption were spontaneous reactions. All adsorption processes were ΔH > 0 and desorption ΔH < 0. It had been reported that adsorption of organic pollutants in a solid–liquid interface was generally caused by a combination of various adsorption forces45.

(a) Effect of temperature on phenanthrene adsorption free energy in the three tested soils. (b) Effect of temperature on phenanthrene desorption free energy in the three tested soils.

In this study, the adsorption rate response value was selected to present the Box-Behnken experimental design and experimental results. The specific results are shown in Table 6.

Experiments were carried out according to the design table, and three kinds of experimental data were obtained respectively. Design-expert software is used to optimize the adsorption rate experiment and treatment. According to the results of the experimental design, constant terms, linear terms (A, B, C, D), interaction terms (AB, AC, AD, BC, BD, CD) and square terms (A2, B2, C2, D2) on the adsorption rate. The corresponding equations for the following second-order polynomials can be derived from the experimental data obtained:

In the formula: A is the concentration of phenanthrene, mg/L; B is pH; C is temperature, °C; D is organic matter, g/kg. A second-order polynomial indicates that the effects of the 4 experimental factors on the response values are interactive, rather than a simple linear relationship (Supplementary Information).

The response surface model is evaluated by variance analysis and significance test to test whether the model can be used to optimize the experimental conditions. The statistical significance of the model equation is determined by the F value, and the significance of each regression coefficient is determined by the P value. The analysis of variance of the model and the fitting results of the quadratic regression equation are shown in Tables 7, 8, 9.

According to Tables 7, 8, 9, the adaptability of the three soil models was very significant (F > 1, P < 0.05). Among them, the organic matter P value of loess is less than 0.0001, indicating that it has a very significant effect on the adsorption rate of loess, and the P value of AD is less than 0.05, which indicates that the interaction between phenanthrene concentration and organic matter has a more significant effect on the adsorption rate, and its complex correlation coefficient of determination R2 = 0.8102, indicating that the fit of the response model is good, and the experimental error is within the acceptable range. The fitting degree of the silt prediction model is R2 = 0.9681, indicating that the model fits well with the experimental results and the experimental accuracy is high; the correction coefficient of determination R2Adj = 0.9363, indicating that about 93.63% of the response value changes can be explained by this model; at the same time, by From the F value, it can be obtained that the order of the influence of various factors on the silt adsorption rate is organic matter > temperature > phenanthrene concentration > pH. In the interaction, the phenanthrene concentration and organic matter have a significant effect on the silt adsorption rate. The coefficient of determination of the silt complex correlation is R2 = 0.9464, indicating that the response model has a good fit, and the experimental error is within the acceptable range. Adjusting the complex correlation coefficient R2 = 0.8982 indicates that the regression relationship can explain 89.82% of the change in the dependent variable. Therefore, this The model can be used to analyze and predict the effect of different factors on the adsorption rate of phenanthrene.

In response surface optimization, the three-dimensional response surface graph reflects the influence of the interaction of the other two variables on the response value, and the slope of the response surface reflects the significance of the interaction of the two variables on the response value. The more significant the interaction effect is on the response value, when the slope is gentle, the effect is not significant. If the contour map is elliptical, it indicates that the interaction between the two variables is significant, and if the contour map is circular, it is not significant46. In addition, the slope and density of the contour line also reflect the influence of the variable on the response value. The steeper the contour line and the greater the density, the greater the influence of the variable on the response value47.

(1) Loess Fig. 5 is a three-dimensional response surface diagram of the interaction between initial phenanthrene concentration and pH to phenanthrene adsorption on loess. It can be seen from the figure that the slope of the response surface graph is steep, and the contour line is an approximate circle, indicating that the interaction between phenanthrene concentration and pH is not significant for the response value. With the increase of pH, the adsorption rate of phenanthrene on loess showed a slow decline at first to the lowest point at 6, and then gradually increased. When the soil pH was close to 6, with the increase of the initial phenanthrene concentration, the adsorption rate of loess also showed a trend of first decreasing and then increasing. According to the F value, F = 0.337, P = 0.5532 > 0.05, it can be concluded that soil pH and initial phenanthrene concentration of the solution have no significant interaction on the adsorption rate of loess.

Figure 6 shows the effects of initial phenanthrene concentration and organic matter on phenanthrene adsorption on Loess under the condition that pH value and temperature are at the central point. It can be seen from the figure that the initial phenanthrene concentration and soil organic matter contour are steep, indicating that their interaction is significant. The range of phenanthrene adsorption rate is 70 ~ 95, and the change of surface is steep. From the Loess error analysis, it can be seen that if f value is 6.05 and P value is 0.0275 < 0.05, the interaction between initial phenanthrene concentration and soil organic matter on Loess adsorption rate is more significant. In the 3D response surface diagram, the influence of soil organic matter on the adsorption rate of loess is more significant, and the adsorption rate of loess also increases with the increase of organic matter content, that is, when the coding value of loess soil organic matter is 1 and the concentration of phenanthrene is 30 mg/L, it is the best reaction condition of loess.

It can be seen from Table 8 that after optimization by the response surface method, the best influencing conditions for the adsorption of phenanthrene in loess are pH = 10, initial phenanthrene concentration of 30 mg/L, temperature of 20.14 °C, soil organic matter content code value of 1, theoretical energy The best adsorption rate achieved was 98.89%48.

(2) Slits Fig. 7 is a three-dimensional response surface diagram of initial phenanthrene concentration and pH to phenanthrene adsorption by silt under central conditions of temperature and organic matter.It can be seen from Fig. 3 that the adsorption rate of silts to phenanthrene decreases first and then increases with the increase of pH value, indicating that the increase of pH is conducive to the collision between phenanthrene molecules and soil, and then the adsorption of phenanthrene to the soil increases. Adsorption rate; with the increase of phenanthrene concentration, the adsorption rate changed little, which may be because the concentration of phenanthrene concentration in the experimental area was higher, which exceeded the adsorption saturation of silts to phenanthrene, so the adsorption rate remained basically unchanged.

Figure 8 effects of initial phenanthrene concentration and organic matter on silts adsorption rate. It can be seen from Fig. 8 that the adsorption rate of silts to phenanthrene increases with the increase of organic matter content, because the increase of organic matter content can provide more adsorption sites to adsorb phenanthrene, thus making the adsorption rate increase; With the increase of concentration, the change of adsorption rate is small, which may be because the concentration of phenanthrene concentration in the experimental area is higher, which exceeds the adsorption saturation of silts to phenanthrene, so that the adsorption rate remains basically unchanged. According to the F value, F = 6.29, P = 0.0251 < 0.05, that is to say, it can be concluded that the initial phenanthrene concentration and organic matter interact significantly on the silts adsorption rate.

Figure 9 shows the effects of temperature and pH on the adsorption rate of silts under the initial phenanthrene concentration and soil organic matter in the central condition. It can be seen from Fig. 9 that the trend of the contour map of temperature and pH is relatively slow, which means that the interaction effect of temperature and pH on the adsorption rate of silts is low; the adsorption rate of silts to phenanthrene decreases with the increase of temperature. With the increase of pH, the adsorption of silts first decreased and then increased. When pH = 6, the adsorption rate of silts was the lowest. It can be seen from Table 8 that the F value of temperature and pH is 2.94, and the P value is 0.1087, indicating that the interaction between temperature and pH is not significant.

It can be seen from Table 10 that after the optimization of the response surface method, the optimal conditions for the adsorption of phenanthrene by the silts are pH = 10, the initial phenanthrene concentration is 30 mg/L, the temperature is 21.15 °C, the soil organic matter content code value is 1, and the theoretical value is 1. The best adsorption rate that can be achieved is 96.59%.

(3) Silty sand When the temperature is 30 °C and the organic matter is 0, the effect of initial phenanthrene concentration and solution pH on the adsorption rate of silty sand is shown in Fig. 10. The trend of silty sand adsorption rate and loess adsorption rate is roughly similar. With the increase of pH value, when the pH value of silty sand is close to 6, the adsorption rate of silty sand is the lowest. With the increase of initial phenanthrene concentration, the adsorption rate of silty sand alsodecreased first and then increased. When pH = 6 and phenanthrene concentration was 20 mg/L,the adsorption rate of silty sand to phenanthrene was the lowest.

Effects of initial phenanthrene concentration and pH on adsorption rate of loess.

Effects of initial phenanthrene concentration and organic matter on the adsorption rate of loess.

Effects of initial phenanthrene concentration and pH on silts adsorption rate.

Effects of soil organic matter and phenanthrene concentration on silts adsorption rate.

Effects of temperature and pH on the adsorption rate of silts.

Effects of initial phenanthrene concentration and pH on silty sand adsorption rate.

Figure 11 shows the effects of temperature and soil organic matter on the adsorption rate of phenanthrene under the condition of pH value and initial phenanthrene concentration at the center point. It can be seen from the figure that the temperature and soil organic matter contour is steep, indicating that the interaction is more significant. The range of adsorption rate of phenanthrene by silty sand is 70 ~ 85, and the surface change is steep. It can be seen from the analysis of silty sand variance that the F value is 3.46 and the P value is 0.0838. It can be seen from Fig. 11 that soil organic matter has a significant effect on the adsorption rate of silty sand, and the adsorption rate of silty sand increases with the increase of organic matter content. That is, when the organic matter coding value of silty sand soil is 1 and the temperature is 20 °C, it is the best reaction condition for silty sand.

Effects of soil organic matter and temperature on silty sand adsorption rate.

Figure 12 shows the effects of initial phenanthrene concentration and temperature on the adsorption rate of phenanthrene under the condition of pH value and organic matter at the central point. It can be seen from the figure that the interaction between initial phenanthrene concentration and temperature is not significant, the range of phenanthrene adsorption rate is 76 ~ 84, and the surface change is gentle. When the temperature is 20 °C, the initial concentration of phenanthrene is 30 mg/L, and the adsorption rate of phenanthrene decreases with the increase of temperature.

Effects of initial phenanthrene concentration and temperature on silty sand adsorption rate.

It can be seen from Table 10 that after optimization by the response surface method, the optimum conditions for the adsorption of phenanthrene in silty sand are pH = 2, the initial phenanthrene concentration is 29.96 mg/L, the temperature is 40 °C, and the code value of soil organic matter content is 1. The optimal adsorption rate that can be achieved theoretically is 93.37%.

In order to verify the reliability of the response model, three groups of experiments were designed for verification. Loess experimental conditions (1) Using the quadratic regression linear equation obtained by the Design Expert8.0 program, the optimal experimental conditions were obtained by simulation: pH = 10, initial phenanthrene concentration of 30 mg/L, temperature of 20.14 °C, soil organic matter content code The value is 1; (2) pH = 2, the initial phenanthrene concentration is 30 mg/L, the temperature is 20 °C, and the soil organic matter content code value is 1; (3) pH = 2, the initial phenanthrene concentration is 10 mg/L, and the temperature is At 20 °C, the coded value of soil organic matter content is − 1. Silts experimental conditions (1) Using the quadratic regression linear equation obtained by the Design Expert8.0 program, the optimal experimental conditions were obtained by simulation: pH = 10, the initial phenanthrene concentration was 30 mg/L, the temperature was 21.15 °C, and the soil organic matter content was The code value is 1; (2) pH = 2, the initial phenanthrene concentration is 30 mg/L, the temperature is 30.12 °C, and the code value of soil organic matter content is 1; (3) pH = 2, the initial phenanthrene concentration is 23.9 mg/L, The temperature is 40 °C, and the coded value of soil organic matter content is 1. Silty sand experimental conditions (1) Using the quadratic regression linear equation obtained by the Design Expert8.0 program, the optimal experimental conditions were obtained by simulation: pH = 2, the initial phenanthrene concentration was 29.96 mg/L, the temperature was 40 °C, and the soil organic matter The content code value is 1; (2) pH = 2, the initial phenanthrene concentration is 29.96 mg/L, the temperature is 40 °C, and the soil organic matter content code value is 0; (3) pH = 2, the initial phenanthrene concentration is 29.96 mg/L, the temperature is 40 °C, the soil organic matter content code value is − 1. In order to investigate the practicability and accuracy of the optimization results, three parallel tests were carried out to verify the adsorption of phenanthrene by soil under the best conditions.

The actual measured average values are as follows. The test results are shown in Table 10.

The model verification results are shown in Table 10. The relative errors between the experimental results and the predicted values of the adsorption rate of phenanthrene in loess, silts and silty sand are all less than 5%, which indicates that the model is suitable and effective. The adsorption process of PAHs has certain guiding significance.

Phenanthrene adsorption and desorption on three soils in Guanzhong Basin were studied using a batch technique. The adsorption and desorption kinetics data were better modeled using pseudo-second-order kinetic equations, indicating that the rate-limiting step was chemical adsorption and the sorption was governed by availability of sorption sites on soil surfaces. Phenanthrene adsorption and desorption isotherms in soils were well correlated with Freundlich models. Moreover, the calculated ΔH and ΔG values indicated that the adsorption process was endothermic reaction and desorption process was exothermic reaction. The influencing factors of adsorption experiment were optimized by response surface method. The results showed that the best influencing conditions for the adsorption of phenanthrene on loess were pH = 10, the initial phenanthrene concentration was 30 mg/L, the temperature was 20.14 °C, the coding value of soil organic matter content was 1, and the best adsorption rate that could be achieved in theory was 98.89%. The optimum conditions for phenanthrene adsorption by silts are pH = 10, initial phenanthrene concentration of 30 mg/L, temperature of 21.15 °C, code value of soil organic matter content of 1, and the best adsorption rate that can be achieved in theory is 96.59%. The best influencing conditions for phenanthrene adsorption by silty sand are pH = 2, the initial phenanthrene concentration is 29.96 mg/L, the temperature is 40 °C, the coding value of soil organic matter content is 1, and the best adsorption rate that can be achieved in theory is 93.37%.

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Lee, C.-C., Chen, C. L. S., Wang, Z.-X. & Tien, C.-J. Polycyclic aromatic hydrocarbons in 30 river ecosystems, Taiwan: Sources, and ecological and human health risks. Sci. Total Environ. 795, 148867–148867 (2021).

Khadhar, S., Higashi, T., Hamdi, H., Matsuyama, S. & Charef, A. Distribution of 16 EPA-priority polycyclic aromatic hydrocarbons (PAHs) in sludges collected from nine Tunisian wastewater treatment plants. J. Hazard. Mater. 183(1–3), 98–102 (2010).

Abdel-Shafy, H. I. & Mansour, M. S. M. A review on polycyclic aromatic hydrocarbons: Source, environmental impact, effect on human health and remediation. Egypt. J. Pet. 25(1), 107–123 (2016).

Wang, Y. et al. Polycyclic aromatic hydrocarbons in the atmosphere and soils of Dalian, China: Source, urban-rural gradient and air-soil exchange. Chemosphere 244(C), 125518 (2020).

Yang, L. et al. Polycyclic aromatic hydrocarbons and nitro-polycyclic aromatic hydrocarbons in five East Asian cities: Seasonal characteristics, health risks and yearly variations. Environ. Pollut. 287, 117360 (2021).

Mo, L. M., Wang, Q. H. & Bi, E. P. Effects of endogenous and exogenous dissolved organic matter on sorption behaviors of bisphenol A onto soils. J. Environ. Manag. 287, 112312 (2021).

Nguyen, T. N. T., Park, M.-K., Son, J.-M. & Choi, S.-D. Spatial distribution and temporal variation of polycyclic aromatic hydrocarbons in runoff and surface water. Sci. Total Environ. 793, 148339 (2021).

Łyszczarz, S., Lasota, J., Szuszkiewicz, M. M. & Błońska, E. Soil texture as a key driver of polycyclic aromatic hydrocarbons (PAHs) distribution in forest topsoils. Sci. Rep. 11(1), 14708–14708 (2021).

Article  PubMed  PubMed Central  Google Scholar 

Saeedi, M., Li, L. Y. & Grace, J. R. Desorption and mobility mechanisms of co-existing polycyclic aromatic hydrocarbons and heavy metals in clays and clay minerals. J. Environ. Manag. 214, 204–214 (2018).

Yin, M., Sun, J., He, H., Liu, J. & Chen, D. Uranium re-adsorption on uranium mill tailings and environmental implications. J. Hazard. Mater. 416(6), 126153 (2021).

Ouassel, S., Chegrouche, S., Nibou, D., Melikchi, R. & Mellah, A. Application of response surface methodology for u(vi) adsorption using hydroxyapatite prepared from eggshells waste materiel: Study of influencing factors, and mechanism. Water Sci. Technol. 83(5), 1198–1216 (2021).

Ouassel , S. , Chegrouche , S. , Nibou , D. , Aknoun , A. & Gisi , SD Adsorption of uranium(vi) onto natural Algerian phosphate: A study of influencing factors and mechanism.Arabic.J. Sci.Eng.46(7), 6645–6661 (2021).

Barkat, M., Nibou, D., Amokrane, S., Chegrouche, S. & Mellah, A. Uranium (vi) adsorption on synthesized 4a and p1 zeolites: Equilibrium, kinetic and thermodynamic studies. C. R. Chim. 18(3), 261–269 (2015).

Schneckenburger, T. & Thiele-Bruhn, S. Sorption of PAHs and PAH derivatives in peat soil is affected by prehydration status: The role of SOM and sorbate properties. J. Soils. Sediments 20, 1–12 (2020).

Wu, Xiao & Shao, Y. X. Study of kinetics mechanism of PAHs photodegradation in solution. Procedia Earth Planet. Sci. 17, 348–351 (2017).

Chandarana, H., Suganya, S. & Madhava, A. K. Surface functionalized Casuarina equisetifolia pine powder for the removal of hetero-polyaromatic dye: Characteristics and adsorption. Int. J. Environ. Anal. Chem. https://doi.org/10.1080/03067319.2020.1798418 (2020).

Osorio-López, C. et al. As (V) adsorption on forest and vineyard soils and pyritic material with or without mussel shell: Kinetics and fractionation. J. Taiwan Inst. Chem. Eng. 45(3), 1007–1014 (2014).

Lončarski, M. et al. Sorption behavior of polycyclic aromatic hydrocarbons on biodegradable polylactic acid and various nondegradable microplastics: Model fitting and mechanism analysis. Sci. Total Environ. 785, 147289 (2021).

Dehghan, P., Abbasi, M., Mofarahi, M. & Azari, A. Adsorption of synthetic and real Kinetic Hydrate Inhibitors (KHI) wastewaters on activated carbon: Adsorption kinetics, isotherms and optimized conditions. Sep. Purif. Technol. 56(13), 2266–2277 (2021).

Wang, P. et al. Adsorption behavior and mechanism of Cd (II) by modified coal-based humin. Environ. Technol. Innov. 23, 101699 (2021).

Vitale, C. M., Terzaghi, E., Zati, D. & Di Guardo, A. How good are the predictions of mobility of aged polychlorinated biphenyls (PCBs) in soil? Insights from a soil column experiment. Sci. Total Environ. 645, 865–875 (2018).

Nguyen, X. P. et al. Effects of pore water chemical composition on the hydro-mechanical behavior of natural stiff clays. Eng. Geol. 166, 52–64 (2013).

Zhang, W. H., Zhuang, L. W., Yuan, Y., Tong, L. Z. & Tsang, D. C. W. Enhancement of phenanthrene adsorption on a clayey soil and clay minerals by coexisting lead or cadmium. Chemosphere 83(3), 302–310 (2010).

Wu, M. J., Liu, H. Y. & Yang, C. P. Effects of pretreatment methods of wheat straw on adsorption of Cd(II) from waterlogged paddy soil. Int. J. Environ. Res. Public Health 16(2), 205–205 (2019).

Article  PubMed  PubMed Central  Google Scholar 

Tegegne, B., Chimuka, L. K., Chandravanshi, B. S. & Zewge, F. Molecularly imprinted polymer for adsorption of venlafaxine, albendazole, ciprofloxacin and norfloxacin in aqueous environment. Sep. Purif. Technol. 56(13), 2217–2231 (2021).

Ren, J. M. Thermodynamics and kinetics of the adsorption for phenanthrene onto CTMA-bentonite in aqueous solutions. Adv. Mater. Res. 1266(233–235), 1089–1092 (2011).

Shi, H. & Sun, Y. P. The adsorption behavior of naphthalene on soil and the influence of temperature were studied. Chin. J. Soil. Sci. 41, 308–313 (2010).

Wang, W. K. et al. Response of the groundwater system in the Guanzhong Basin (central China) to climate change and human activities. Hydrogeol. J. 26(5), 1–13 (2018).

He, H. X. et al. Surface characteristics of a recyclable and efficient adsorption material. Mater. Lett. 256(11), 126658 (2019).

Liang, X. R. et al. Study of the sorption/desorption behavior of chlortetracycline on sediments in the upper reaches of the Yellow river. Chem. Eng. J. 428, 131958 (2022).

Peng, H., Jiang, Anna & Wang, H. Adsorption and desorption characteristics of phosphorus on sediments in Panzhihua section of Jinsha river China. IOP Conf. Ser. 651(4), 042059 (2021).

Denys, A. S. R., Jaqueline, M. M., Guilherme, L. D., Tito, R. S. C. & Luiz, A. A. P. Preparation, characterization and dye adsorption/reuse of Chitosan-Vanadate films. J. Polym. Environ. 26(7), 2917–2924 (2018).

Yue, L. et al. Adsorption–desorption behavior of atrazine on agricultural soils in China. J. Environ. Sci. 57(07), 180–189 (2017).

Wang, P., Wu, T. H. & Zhang, Y. In situ investigation the photolysis of the PAHs adsorbed on mangrove leaf surfaces by synchronous solid surface fluorimetry. PLoS ONE 9(1), e84296 (2014).

Article  PubMed  PubMed Central  Google Scholar 

Qu, Y. Y., Li, H., Li, A., Ma, F. & Zhou, J. T. Identification and characterization of Leucobacter sp. N-4 for Ni (II) biosorption by response surface methodology. J. Hazard. Mater. 190(1), 869–875 (2011).

Shanmugaprakash, M., Sivakumar, V., Manimaran, M. & Aravind, J. Batch and dynamics modeling of the biosorption of Cr(VI) from aqueous solutions by solid biomass waste from the biodiesel production. Environ. Prog. Sustain. Energy 33(2), 342–352 (2014).

Wu, W. H., Huang, Y., Lin, D. H. & Yang, K. Sorption mechanism of naphthalene by diesel soot: Insight from displacement with phenanthrene/p-nitrophenol. J. Environ. Sci. 106(08), 136–146 (2021).

Asadullah, Kaewsichan, L. & Tohdee, K. Adsorption of hexavalent chromium onto alkali-modified biochar derived from Lepironia articulata: A kinetic, equilibrium and thermodynamic study. Water Environ. Res. 91(11), 1433–1446 (2019).

Pirsalami, S., Bagherpour, S., Ebrahim, B. M. & Riazi, M. Adsorption efficiency of glycyrrhiza glabra root toward heavy metal ions: Experimental and molecular dynamics simulation study on removing copper ions from wastewater. Sep. Purif. Technol. 275, 119215 (2021).

Li, W. B. et al. Composite modification mechanism of cationic modifier to amphoteric modified kaolin and its effects on surface characteristics. Int. J. Environ. Sci. Technol. 13(11), 2639–2648 (2016).

Qian, J. et al. Adsorption of perfluorooctane sulfonate on soils: Effects of soil characteristics and phosphate competition. Chemosphere 168, 1383–1388 (2017).

Wang, T. F. et al. Effect mechanism of litter extract from Alternanthera philoxeroides on the selective absorption of heavy metal ions by amphoteric purple soil. J. Environ. Manag. 321, 115970–115982 (2022).

He, H. X. et al. Adsorption of phenanthrene on the magnetized biochar organic-coatedby amphoteric modifier. Desalin. Water Treat. 148(3), 195–201 (2019).

Wang, W. et al. Adsorption of emerging sodium p-perfluorous nonenoxybenzene sulfonate (OBS) onto soils: Kinetics, isotherms and mechanisms. Pedosphere 31(4), 596–605 (2021).

Wu, Z. S. et al. Competitive adsorption of naphthalene and phenanthrene on walnut shell based activated carbon and the verification via theoretical calculation. RSC Adv. 10(18), 10703–10714 (2020).

Article  PubMed  PubMed Central  Google Scholar 

Zhou, R. J., Zhang, M., Li, J. Y. & Zhao, W. Optimization of preparation conditions for biochar derived from water hyacinth by using response surface methodology (RSM) and its application in Pb2+ removal. J. Environ. Chem. Eng. 8(5), 104198 (2020).

Ojeda-López, R. et al. Tailoring synthesis conditions of carbon microfibers to enhance the microporosity, CO2 and CH4 adsorption by using the response surface methodology. Microporous Mesoporous Mater. 305, 110333 (2020).

Nazari, F., Ghoreishi, S. M. & Khoobi, A. Bio-based Fe3O4/chitosan nanocomposite sensor for response surface methodology and sensitive determination of gallic acid. Int. J. Biol. Macromol. 160, 456–469 (2020).

This study was supported by the National Natural Science Foundation of China (Program No.41202176, 42177076), Key Research and Development Program of Shaanxi (Program No.2021SF-445, 2021ZDLSF05-09) and Technology Innovation Leading Program of Shaanxi (Program No.2022QFY06-04). We thank professor Zhao for revising our manuscript.

College of Geology and Environment, Xi’an University of Science and Technology, Xi’an, 710054, China

Hua Tian, ​​Qing Zhang, Zu-feng Xie, Fang Pu & Qian-ji Wang

Shaanxi Provincial Key Laboratory of Green Coal Development and Geological Guarantee, Xi’an, 710054, China

Zhongsheng Environmental Technology Development Co, Ltd, Xi’an, 710054, China

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

H.T.: Conceptualization, Methodology, Investigation, Resources, Writing—review & editing, Project administration, Funding acquisition. Q.Z.: Conceptualization, Methodology, Investigation, Formal analysis, Writing—original draft. X.T.: Methodology, Investigation. Z.-f.X.: Data curation. F.P.: Methodology. Q.-j.W.: Analysis, Validation.

The authors declare no competing interests.

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

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

Tian, H., Zhang, Q., Tian, X. et al. Analysis of influencing factors of phenanthrene adsorption by different soils in Guanzhong basin based on response surface method. Sci Rep 12, 20906 (2022). https://doi.org/10.1038/s41598-022-25293-0

DOI: https://doi.org/10.1038/s41598-022-25293-0

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Scientific Reports (Sci Rep) ISSN 2045-2322 (online)

Analysis of influencing factors of phenanthrene adsorption by different soils in Guanzhong basin based on response surface method | Scientific Reports

Conical Centrifuge Tubes Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.