AbstractThis study aimed to develop a robust predictive model for tetracycline (TC) adsorption onto biochar (BC) by employing machine learning techniques to investigate the underlying driving factors. Four machine learning algorithms, namely Random Forest (RF), Gradient Boosting Decision Tree (GBDT), eXtreme Gradient Boosting (XGBoost) and Artificial Neural Networks (ANN), were used to model the adsorption of TC on BC using the data from 295 adsorption experiments. The analysis revealed that the RF model had the highest predictive accuracy (R2 = 0.9625) compared to ANN (R2 = 0.9410), GBDT (R2 = 0.9152), and XGBoost (R2 = 0.9592) models. This study revealed that BC with a specific surface area (S (BET)) exceeding 380 cm3·g−1 and particle sizes ranging between 2.5 and 14.0 nm displayed the greatest efficiency in TC adsorption. The TC-to-BC ratio was identified as the most influential factor affecting adsorption efficiency, with a weight of 0.595. The concentration gradient between te dsorbate and adsorbent was demonstrated to be the principal driving force behind TC adsorption by BC. A predictive model was successfully developed to estimate the sorption performance of various types of BC for TC based on their properties, thereby facilitating the selection of appropriate BC for TC wastewater treatment.
IntroductionTetracycline (TC) is extensively employed as an antimicrobial agent and feed supplement in agriculture and animal husbandry1. However, researchers have recently paid significant attention to the issues of incomplete metabolism and TC emissions2,3. Due to its persistence as an organic pollutant, TC is frequently detected in surface water, groundwater, and drinking water. TC can induce endocrine disruption in target organisms and can also contribute to the dissemination of antibiotic resistance genes, thereby posing serious human health concerns and environmental hazards4,5. Given the inhibitory effect of TC on mcroorganisms, the removalof TC from water bodies using conventional biological water treatment methods proves challenging6.Currently, the principal methods for treating TC in wastewater include chemical oxidation, biological treatment, and physical removal2. Adsorption, on account of its inherent advantages, such as simplicity, low cost, and high efficiency, is viewed as an excellent technology for the treatment of TC. Among the various adsorbents, biochar (BC) has been extensively researched as an adsorbent for removing pollutants from wastewater due to its unique characteristics, such as a large specific surface area, uniform pore distribution, and high concentration of surface functional groups7.The uptake of tetracycline onto biochar mainly involves physical interactions such as van der Waals forces and hydrogen bonding, as well as chemical reactions including covalent and ionic bonding8. Therefore, the adsorption process primarily depends on the properties of biochar, adsorption onditions, and the ratio o adsorbate to absorbent. Several traditional kinetic and isothermal adsorption models have been extensively evaluated in previous studies9,10,11. Findings suggest that the possible adsorption mechanisms include π–π interactions, electrostatic interactions, and chemisorption. Although a typical controlled-variable experimental approach can determine the relationship between each influencing factor and the amount of sorption within the same framework, traditional batch sorption experiments are time-consuming and inefficient when selecting suitable biochar12. Therefore, there is an urgent need to develop practical tools for predicting adsorption efficiency, optimizing process parameters, and comprehending the adsorption mechanism.Machine learning (ML)-assisted modeling has been proposed as a potential approach to reduce the cost and time associated with laboratory contaminant removal processes. Previous research has utilized machine learning (ML) algorithms on slected carbon-based materias to adsorb tetracycline (TC)13,14, yet the accuracy of the models could be enhanced. Zhu et al.’s13 study employed carbon-based materials such as activated carbon and biochar, which have distinct compositions. Thus, developing prediction models for both materials represents a significant challenge due to the high variability; in addition, the study has a limited database, and the highest achieved R2 value was only 0.8944, highlighting the necessity to optimize the machine learning (ML) model15,16. This study will also use the results of the machine learning model to explore the driving factors for the adsorption of tetracycline on biochar. To evaluate the prediction effectiveness, generalized adsorption models must be utilized to predict TC adsorption on a single biochar, particularly integrated learning models. Ensemble learning is a typical ML algorithm that integrates the modeling outcomes of all models by building multiple models from the data17. The mot typical ensemble learning lgorithms used in assessing TC adsorption on biochar include random forest (RF), gradient boosting decision tree (GBDT), and eXtreme gradient boosting (XGBoost)18. In addition to ensemble learning, this study will also incorporate the most popular deep learning algorithms as a point of comparison.The integration of machine learning as an advanced algorithm within the field of environmental remediation employing biochar remains in its nascent stage, considering the widespread occurrence, substantial ecological risk, and unique properties of toxic compounds (TC) in the environment. This research was conducted with the aims of: (i) devising universal machine learning models to forecast the sorption capacity of TC on biochar (BC), contingent on BC attributes and sorption conditions; (ii) investigating the primary factors contributing to BC adsorption of TC; (iii) assessing the impact of various factors on the relative significance of BC sorption capacity and asertaining the combined effectof each factor on BC sorption capacity; and (iv) constructing an accessible web-based user interface for engineers. The machine learning-driven model devised in this investigation establishes a theoretical foundation for the pragmatic treatment of TC, delivering an all-encompassing comprehension of TC sorption on biochar relative to its features and sorption milieu.Materials and methodsCum biochar sorption capacity predictions layoutExperimental data for the adsorption of tetracycline by biochar were collected from ten papers, including 22 biochar species and 295 sets of experimental adsorption data9,19,20,21,22,23,24,25,26,27,28. Without author bias, the related articles were selected randomly and data were extracted from published papers using Plot Digitizer v3 (https://plotdigitizer.com/#download)29. Detailed data are provided in the supplementary materials (Tables S1, S2).To predict the sorption capacity of BC for TC, expressed as the equilibrium sorpton capacity Qe (mg·g−1), 12 crtical factors were considered and divided into three categories: (i) biochar properties: Brunauer–Emmett–Teller surface area [S (BET), m2·g−1], pH of the biochar in water (pH_H2O), total carbon in the biochar (C, w%), molar ratio of oxygen and nitrogen to carbon [(O + N)/C], molar ratio of oxygen to carbon (O/C), molar ratio of hydrogen to carbon (H/C), ash content (Ash, w%), pore volume (V, cm3·g−1), and biochar pore diameter (D, nm); (ii) adsorption conditions: adsorption temperature (T, °C) and solution pH (pH_sol); and (iii) initial concentration ratio of tetracycline to biochar (C0, mmol·g−1). In the ML section, data not provided in the published paper were replaced with K-Nearest Neighbor (KNN) Algorithm, while pH_H2O and ash were missing more often than the other values and was decided to remove them. The TC characteristics are listed in Table S3.The following equation was used to obtain C030:$$egin{array}{c}{ ext{C}}_{0} ext{=} , frac{{ extC}}_{ ext{TC}}/444.4}{{ ext{C}_{ ext{BC}}},end{array}$$
(1)
where ({ ext{C}}_{ ext{TC}}) (mg·L−1) is the initial concentration of TC and ({ ext{C}}_{ ext{BC}}) (g·L−1) is the initial concentration of BC.Pre-processing of dataThe linear correlation between any two randomly selected variables or between variables and target values was measured using the Pearson correlation coefficient using the following equation30,31:$$egin{array}{c}{r}_{xy} = frac{{sum }_{i=1}^{n}left({x}_{i}-overline{x }
ight){sum }_{i=1}^{n}left({y}_{i}-overline{y }
ight)}{sqrt{{sum }_{i=1}^{n}{left({x}_{i}-overline{x }
ight)}^{2}}sqrt{{sum }_{i=1}^{n}{left({y}_{i}-overline{y }
ight)}^{2}}}, end{array}$$
(2)
where (overline{x }) and (overline{y }) denote the mean of variable x or y.Construction of ML modelsEnsemble learning is a popular machine-learning algorithm that integrates multiple models (base estimators) to form an ensemble estimator that soves complex problems with specific rules17,32. Itegrated models act as an integrated platform to automatically manage the weaknesses and enhancements of individual models to achieve higher prediction accuracy. Three integration algorithms exist: bagging, boosting, and stacking.RF is a representative bagging integration algorithm and is the most commonly used algorithm for predicting poorly understood processes33. The outcome of an RF prediction is a combination of the predicted outcomes of each decision tree, so the critical step in Random Forest prediction is the formation of a decision tree and a forest. This principle is depicted in Fig. S1.The Gradient Boosted Decision Tree (GBDT) algorithm, an iterative decision tree algorithm, consists of multiple decision trees, and the conclusions of all the trees are summed to arrive at the final answer. As shown in Fig. S2, the GBDT algorithm uses the negative gradient value of the loss function of the base model in round i as an approximatin of the loss value of the base model in that roud34. The next step is to construct round i + 1 of base models based on this approximation to make the solution of the objective function more convenient34.XGBoost is an improved version of GBDT. It has an engineering goal of pushing the computational power of boosting trees to a limit to achieve fast computation and superior performance35. With many improvements over traditional gradient boosting algorithms, XGBoost can be performed faster than other comprehensive algorithms that use gradient boosting and is recognized as an advanced evaluator with ultrahigh performance in both classification and regression.In this study, Artificial Neural Network (ANN) was utilized due to its ability to simulate the connections and signal propagation between neurons, allowing the adjustment of weights using the backpropagation algorithm to learn patterns and relationships in the data. ANN consists of an input layer, hidden layers, and output layer, makng it suitable for various tasks such as predictio, classification, and pattern recognition.All machine-learning algorithm codes were obtained from the open-source Scikit-learn library. All datasets were divided into training and test data at a ratio of 70:30, with the random state set to 40. Tenfold cross-validation was used to select the best hyperparameters from the data. The test data were used to evaluate model performance. All data were normalized before training. All input and output parameters in this study are listed in Table S4.Modeling performance evaluationThe performance of the model was assessed using the coefficient of determination (R2) and the root mean squared error (RMSE)13,30,35.$$egin{array}{c}{R}^{2}=1- frac{{sum }_{i=1}^{N}{left({Y}_{i}^{exp}- {Y}_{i}^{pred}
ight)}^{2}}{{sum }_{i=1}^{N}{left({Y}_{i}^{exp}- {Y}_{ave}^{exp}
ight)}^{2}},end{array}$$
(3)
$$egin{array}{c}RMSE= sqrt{frac{1}{N}sum_{i=1}^{N}{left({Y}_{i}^{exp}-{Y}_{i}^{pred
ight)}^{2},}end{array}$$
(4)
where ({Y}_{i}^{exp}) and ({Y}_{i}^{pred}) are the experimental and predicted values, and ({Y}_{ave}^{exp}) is the average of the experimental values.Results and discussionStatistical results of biochar characteristicsThis study utilized a combination of box plots and normal distribution curves to illustrate the distribution patterns of continuous data (see Fig. 1). The composite plot comprises two sections—the left segment illustrates the box plot, whereas the right segment manifests the normal distribution curve of the data. The box plot depicts the median, signified by a white dot, the interquartile range, denoted by the box, and the whiskers, which represent the remaining data. Outliers are designated by circular points or alternative symbols. The probability density of data at each value is displayed on the right portion of the plot, with elevated values indicating a comparatively higher probability of data occurrence at that poit. Please refer to Table S5 for specific values.Figure 1Visualizatio of biochar properties through box plots and normal distribution curves.Full size imageThe majority of biochars exhibit alkalinity in water, a phenomenon predominantly correlated with the depletion of acidic functional groups and an augmentation in aromatic carbon at elevated temperatures, but also linked to the accumulation of alkaline ions such as Na+, Ca2+, K+, and Mg2+ in the biochar30. The mean pH of all biochar samples was 9.14. Nevertheless, a minority displayed weak acidity, which can be accounted for by the incomplete liberation of alkali salts from the biochar matrix at lower pyrolysis temperatures36,37. Given that pH exerts a significant impact on the adsorption of tetracycline onto biochar, the pH of the solution was adjusted using either an acid or a base in all experiments selected for this investigation. Consequently, subsequent analyses will disregard the pH of the biochar in an aqueous solution (pH_H2) and instead utilize the pH of the solution (pH_sol) during batch adorption studies.As depicted in Fig. 1, biochar exhibits a high carbon content, with an estimated average of approximately 60% and a maximum value of 92%. Existing research has demonstrated that the carbon content escalates in correlation with increasing pyrolysis temperatures within a specified range20,23,38. Consequently, the pyrolysis process concentrates carbon within the biomass feedstock39,40. The H/C, O/C, and (N + O)/C ratios serve as indicators of aromaticity, hydrophilicity, and polarity indices, respectively39,40. The H/C, O/C, and (N + O)/C ratios indicate the aromaticity, hydrophilicity, and polarity indices, respectively41. A lower H/C ratio in BC corresponds to higher aromaticity; a lower O/C ratio indicates reduced hydrophilicity; and an elevated (N + O)/C ratio signifies increased polarity42,43. The median O/C and H/C ratios were 0.14 and 0.2, respectively. According to the International Biochar Inititive (IBI) Standards, the H/C ratio for biochar should be less than 0.. Therefore, values exceeding 0.7 in the dataset can be eliminated to enhance the accuracy of biochar data analysis. Nevertheless, Table S2 contains only limited data.It has been posited that an increased S (BET) is advantageous for TC adsorption onto biochar, while D and V exhibit minimal direct influence, though an optimal range for each exists. The ash content of biochar samples exhibited a broad range, spanning from 1.50% to 55.27%, attributable to variations in feedstock type and pyrolysis conditions, which alter the physicochemical properties and spatial distribution of organic matter30. However, the role of ash in TC adsorption onto biochar remains a subject of debate.Statistical outcomes of data correlation analysisFigure 2 demonstrates that there exists a significant positive correlation between C0, S (BET), and V with respect to Qe. This relationship can be elucidated by the transfer equation: when the adsrption value remains constant, an increase in the amount of adsorbent lads to a higher adsorption capacity per unit mass of adsorbed material. Research conducted by Wang et al., Zhu et al., and Kim et al.20,23,30 substantiates the connection between S (BET) and Qe, suggesting that an elevated specific surface area permits more adsorbates to be adsorbed per unit mass of the adsorbent. Additionally, a higher number of pores per unit of adsorbent contributes to an increased adsorption capacity, as evidenced by the correlation between V and Qe.Figure 2Correlation coefficients and corresponding significant levels.Full size imageMoreover, in accordance with prior research findings30, S (BET) displays a positive correlation with carbon content, yet a negative correlation with ash content. This observation implies that elevated carbon content results from the removal of volatile substances, while a higher ash content can cause micropore filling, consequently reducing the surface area6,13. Furhermore, the inverse correlation between S (BET) and D lends support to revious studies suggesting that an increase in S (BET) corresponds to a decrease in D. The correlation coefficient between ash and carbon was − 0.91, and the correlation coefficients between (O + N)/C and O/C were 1. Thus, one variable from each pair must be excluded to prevent collinearity.Selection of machine models for tetracycline adsorption on biocharIn order to assess the potential of machine learning in predicting the adsorption of tetracycline on biochar, a tenfold cross-validation was executed to determine hyperparameters. Figure 3 depicts the two most significant parameter combinations, ‘n_estimators’ and ‘max_depth’ (‘hidden layer’ and ‘learning rate’) with the color on the surface model representing the model’s performance quality. The redder the color, the higher the model’s accuracy, as per the hyperparameter selection principle. The parameters selected in this study were determined to be equal to orgreater than the threshold value of two.Figure 3Schematic illustration ofhyperparameter selection (a RF, b GBDT, c XGBoost, d ANN).Full size imageAs shown in Table S4, the highest mean R2 value of 0.9625 was achieved by Random Forest (RF), followed by XGBoost at 0.9592, while Gradient Boosted Decision Trees (GBDT), had the lowest score of 0.9152, and the artificial neural network (ANN) had the second-to-last score of 0.9410. The RMSE values for RF, GBDT, ANN, and XGBoost were 18.02, 25.97, 22.42, and 20.99, respectively. As illustrated in Fig. 4, both the training and test sets exhibited R2 and RMSE values of 0.9703, 0.9625, and 14.91, and 18.02, respectively, for RF and XGBoost. The model did not appear to be overfitted. Although XGBoost is generally regarded as the superior model, in this study, RF outperformed XGBoost. One possible explanation is that the results were influenced by the features and the nature of the problem. Therefore, while XGBoost is among the better models, it my not always be the optimal choice, exemplifying the “no free lunch theore”.Figure 4Scatter plot of RF-model-predicted adsorption values and experimental data (a training data, b testing data).Full size imageAccording to Fig. 3, the predicted values for all data ranged from 0 to 350 mmol·g−1. The RMSE of the RF model (18.02 mmol·g−1) accounted for 0.0542 of the predicted range (0 to 350 mmol·g−1), indicating a high level of accuracy. Combining the accuracy of R2 and RMSE, the RF model was selected for subsequent analysis since it provided higher accuracy than the XGDT and XGBoost models. The best parameters for each model are listed in Table S4.Feature importance analysisThis section is based on the RF model. Feature importance analysis serves as a potent instrument for discerning the relevance of input features in target prediction. Employing machine learning to comprehend tetracycline (TC) adsorption on biochar (BC) can substantially curtail the time-consuming and costly experimenta design process by leveraging feature importance to select a limited numberof features for model training, thereby reducing time and cost while enhancing accuracy44,45. Although machine learning models are formidable tools for generating precise predictions, they frequently function as “black box” models, rendering the comprehension of their inner mechanisms and decision-making processes arduous. Nevertheless, feature importance analysis offers an efficacious approach for pinpointing the most crucial input variables in a model and comprehending their contributions to the output.A primary advantage of employing SHAP (SHapley Additive exPlanations) for feature importance analysis is its ability to furnish a visual representation of each feature’s contribution to the output prediction. Figure 5a exhibits a SHAP feature importance visualization, supplying an in-depth dissection of the weight and influence exerted by each input feature on the predicted outcome. This information can be harnssed to identify vital control parameters, optimize the experimental design rocess, and bolster model accuracy.Figure 5Relative importance of input variables on sorption capacity using SHAP (a SHAP value; b feature importance).Full size imageFigure 5b displays the specific values of feature importance derived from SHAP. The results suggest that C0 (0.695) is the most critical factor affecting Qe, signifying that BC’s adsorption capacity for TC predominantly depends on the TC-to-BC ratio. This phenomenon, referred to as the value transfer process, implies that the concentration gradient between the adsorbate and the adsorbent constitutes the principal driving force for TC adsorption by BC.In examining the biochar characteristics, surface area (S (BET)), (O + N)/C, and H/C ratios displayed notably significant effects (0.162, 0.036, and 0.032, respectively), suggesting that S (BET) is the most crucial factor influencing biochar’s adsorption properties. It can be deduced that the sites prvided by S (BET) also play a vital role as driving factors in the adsorption f target compounds (TC) by biochar (BC). Concerning adsorption conditions, temperature (T) and pH_sol contributed to 0.14 and 0.06 of the characteristic importance, respectively. The impact of each factor on Qe is discussed in “Analysis of partial dependence plots (PDP)” section.Analysis of partial dependence plots (PDP)Figure 6(A1) presents the single-factor PDP, which reveals a partial dependence of the initial concentration (C0) on the equilibrium adsorption capacity (Qe), demonstrating an initial increase in the adsorption rate followed by stabilization. This trend can be attributed to the gradual filling of adsorption sites on the biochar surface as the relative content of TC increases. Upon reaching C0 = 2 mmol·g−1, the adsorption sites on the biochar surface become fully occupied, resulting in maximum adsorption capacity. The saturation of adsorption sites limits any further increase in the removal rat. These findings highlight the critical role of biochar surface area and capacty in effective pollutant removal21.Figure 6The PDP of TC adsorption on significant variables (spikes on the x-axis represent data density).Full size imageBiochar’s adsorption efficacy is contingent upon its specific surface area (S (BET)). As depicted in Fig. 6(A2), a rapid rise in adsorption capacity is observed, followed by a gradual increase. The enhancement in the number of sorption sites with an increase in S (BET) leads to a higher sorption uptake capacity of biochar. Nonetheless, the adsorption efficacy of biochar is constrained by other factors, such as the initial concentration of the target compound (C0), when there is a certain increase in the number of adsorption sites20,21. S (BET) predominantly affects the driving force of chemisorption occurring between BC and TC. A larger S (BET) results in a more significant number of chemisorption sites, thus leading to higher adsorption efficiency. Furthemore, the indirect influence of biochar’s specific surface area on physical adsrption should not be overlooked. The adsorption of the target compound (TC) on biochar is the result of various driving forces. Based on Fig. 6(A2), it can be inferred that biochar with an S (BET) greater than 380 cm3·g−1 offers superior adsorption efficiency.The distribution of tetracycline (TC) species is influenced by the pH value. Table S3 presents the dissociation constants (pKas) of TC as 3.3, 7.7, and 9.7. Within the pH range of 2–3.3, TC+ is the predominant form of TC; at pH 3.3–7.7, TC0 is prevalent; at pH 7.7–9.7, TC− dominates, and at pH above 9.7, TC converts to TC2−. Figure 5(B1) illustrates that the adsorption of TC on biochar (BC) is most favorable when the solution pH (pH_sol) is approximately 5.5. Intriguingly, these findings contradict a report by Ref.21, which posited that biochar and tetracycline primarily undergo electrostatic adsorption. Despite the biochar’s negative charge, it assume a positive charge when the pH ranged between 3 and 3.3. Figure 6(A3) indicates hat, in addition to electrostatic adsorption, hydrogen bonding and π–π electron donor–acceptor interactions influence TC adsorption on biochar19,20,21.The adsorption efficacy of biochar is also partially dependent on the D value. Figure 6(B2) indicates that adsorption efficiency increases, stabilizes, and then significantly decreases. Biochar with a D value of 2.5–14.0 nm demonstrates a higher tendency to adsorb TC. Nonetheless, this change is insignificant when compared to other factors. This study affirms the findings of Zhang et al.21 that the adsorption performance of adsorbents is optimal when the molecular size of the adsorbent’s D is 1.7–3.0 times larger than that of the adsorbate. However, these outcomes should be interpreted with caution.Data for this investigation were gathered at temperatures between 15 and 40 °C (Fig. 6(C1)), with temperature fluctuations displaying minor, partially dependent vriations from 15 to 35 °C. The adsorption efficiency’s partial dependence on temprature increases substantially within a higher temperature range of 35–40 °C, which aligns with previous findings that the adsorption process is thermodynamically favorable at elevated temperatures6,24,25. This phenomenon may be attributed to the diffusion of TC and enhanced interfacial chemistry.The adsorption efficiency’s partial dependence on the chemical composition factors C, H/C, and (O + N)/C is relatively insignificant (Fig. 6(C2–D2)) and exhibits considerably less variation compared to changes in other factors. Consequently, this aspect will not be further explored in the present study.A two-factor analysis was also employed to investigate the impact on TC adsorption, a representative depiction of which is provided in Fig. 7. (i) At a fixed specific surface area (S (BET)), the partial dependence increases with initial concentration (C0) and tends to stabilize when C0 > 2 mmol·g−1. (ii) The effct of pH_sol was less pronounced when comparing the partial dependences of C0 and H_sol. (iii) At a fixed C0, the partial dependence tends to rise sharply when the temperature exceeds 35 °C.Figure 7Bivariate PDP of TC adsorption on any two significant input variables and the interaction between the two variables.Full size imageSignificance and drawbacks of this studyThe investigation of biochar preparation for tetracycline (TC) adsorption employing cost-effective biomass has recently emerged as a focal point of research, owing to its robust capacity to eliminate organic contaminants from aqueous media. Typically, controlled-variable experimental methodologies are employed to ascertain factors such as material characteristics and adsorption conditions. Nonetheless, conventional batch experiments are labor-intensive, expensive, and lack generalizability. In this study, Random Forest (RF) was demonstrated as a beneficial machine learning (ML) instrument for predicting the quantity of TC dsorbed by biochar, thereby showcasing its potential to directly forecast experimenal outcomes based on pre-established conditions. Moreover, discerning the most crucial factors impacting TC sorption and their influence on the process offers invaluable insights for selecting or devising TC removal techniques. Consequently, the requisite number of experiments can be considerably diminished, and the exploration of biochar (BC) applications for TC adsorption can be markedly expedited.This research revealed that the RF algorithm provides a reasonable prediction of TC adsorption quantities by biochar. However, the predictive performance of ML was hindered by data imbalance and scarcity, and several considerations must be addressed. The data utilized for model training solely predicted the sorption of TC by biochar, without accounting for the sorption of other antibiotics by the same material.ConclusionThis study successfully established a model for biochar adsorption of tetracycline (TC) uing a comprehensive learning-based approach. The principal findings are as follows:
(1)
The Random Forest (RF) algorithm proved to be an accurate and effective predictor of TC adsorption by biochar (BC), achieving a coefficient of determination (R2) value of 0.9625, slightly outperforming the GBDT, ANN, and XGBoost algorithms.
(2)
The ratio of tetracycline to biochar significantly influenced the sorption process, with a weight of 0.595.
(3)
The primary driving force for the adsorption of TC by BC is the concentration gradient between the adsorbate and the adsorbent.
(4)
Biochar with initial concentration (C0) greater than 2 mmol·g−1, specific surface area (S (BET)) exceeding 380 cm3·g−1, and adsorbent diameter (D) ranging from 2.5 to 14.0 nm exhibited the highest propensity for adsorbing TC.
The model developed in his study has significant implications for minimizing redundant experimentation and facilitating the selection of appropriate biochar. Moreover, it will guide the proper application of iochar in tetracycline wastewater treatment technologies.
Data availability
All data generated or analysed during this study are included in this published article and its Supplementary Information files.
AbbreviationsML:
Machine learning
RF:
Random forest
RMSE:
Root mean square error
PCC:
Pearson correlation coefficient
GBDT:
Gradient boosting decision tree
XGBoost:
XExtreme gradient boosting
D:
Particle size
S (BET) (m2·g−1):
Brunauer–Emmett–Teller surface area
V (cm3·g−1):
Total pore volume
C0 (mmol·g−1):
Initial concentration ratio of tetracycline to biochar
Qe (mg·g−1):
Equilibrium adsorption capacity of tetracyclie on biochar
TC:
Tetracycline
BC:
Biochar
PDP:
Partial dependence plots
C:
Total carbon in the biochar
pH_H2O:
pH of the biochar in water
pH_sol:
Solution pH
(O + N)/C:
Molar ratio of oxygen and nitrogen to carbon
H/C:
Molar ratio of hydrogento carbon
Ash:
Ash content
ANN:
Artificial neural networks
ReferencesGopal, G., Alex, S. A., Chandrasekaran, N. & Mukherjee, A. A review on tetracycline removal from aqueous systems by advanced treatment techniques. RSC Adv. 10, 27081–27095. https://doi.org/10.1039/d0ra04264a (2020).Article
ADS
CAS
PubMed
PubMed Central
Google Scholar
Phoon, B. L. et al. Conventional and emerging technologies for removal of antibiotics from wastewater. J. Hazard. Mater. 400, 122961. https://doi.org/10.1016/j.jhazmat.2020.122961 (2020).Article
PubMed
Google Scholar
eng, G., Liu, Y., Ma, X. & Fan, Y. Fabrication of magnetic multi-template molecularly imprinted polymer composite for the selective and efficient removal of tetracyclines from water. Front. Environ. Sci. Eng. 15, 1–12. https://doi.org/10.1007/s11783-021-1395-5 (2021).Article
CAS
Google Scholar
Bilal, M., Mehmood, S., Rasheed, T. & Iqbal, H. M. Antibiotics traces in te aquatic environment: Persistence and adverse environmental impact. Curr. Opin. Environ. Sci. Health 13, 68–74. https://doi.org/10.1016/j.coesh.2019.11.005 (2020).Article
Google Scholar
Zhang, X., Yan, S., Chen, J., Tyagi, R. & Li, J. 3-Physical, chemical, and biological impact (hazard) of hospital wastewater on environment: Presence of pharmaceuticals, pathogens, and antibiotic-resistance genes. Biotechnol. Bioeng. https://doi.org/10.1016/B978-0-12-819722-6.00003-1 (2020).Article
PubMed
PubMed Central
Google Scholar
Zhu, T. et al. Insights intothe fate and removal of antibiotics and antibiotic resistance genes using biological wastewater treatment technology. Sci. Total. Environ. 776, 145906. https://doi.org/10.1016/j.scitotenv.2021.145906 (2021).Article
ADS
CAS
Google Scholar
Akhil, D. et al. Production, characterization, activation and environmental applications of engineered biochar: A review. Environ. Chem. Lett. 19, 2261–2297. htts://doi.org/10.1007/s10311-020-01167-7 (2021).Article
CAS
Google Scholar
Thangaraj, B. & Solomon, P. R. Immobilization of lipases—A review. Part I: Enzyme immobilization. ChemBioEng Rev. 6, 157–166. https://doi.org/10.1002/cben.201900016 (2019).Article
CAS
Google Scholar
Chen, T. et al. Sorption of tetracycline on H3PO4 modified biochar derived from rice straw and swine manure. Bioresour. Technol. 267, 431–437. https://doi.org/10.1016/j.biortech.2018.07.074 (2018).Article
CAS
PubMed
Google Scholar
Jang, H. M. & Kan, E. Engineered biochar from agricultural waste for removal of tetracycline in water. Bioresour. Technol. 284, 437–447. https://doi.org/10.1016/j.biortech.2019.03.131 (2019).Article
CAS
PubMed
Google Scholar
Liu, C. et al. Response surface methodology for the optimization of the ultrasonic-assisted rhamnolipid treatment of oily sludge. Arab. J. Chem. 14, 102971. https://doi.org/10.1016/j.arabjc.2020.102971(2021).Article
CAS
Google Scholar
Li, X. et al. Characterization of biochars from woody agricultural wastes and sorption behavior comparison of cadmium and atrazine. Biochar 4, 1–12. https://doi.org/10.1007/s42773-022-00132-7 (2022).Article
CAS
Google Scholar
Zhu, X. et al. Machine learning for the selection of carbon-based materials for tetracycline and sulfamethoxazole adsorption. Chem. Eng. J. 406, 126782. https://doi.org/10.1016/j.cej.2020.126782 (2021).Article
CAS
Google Schoar
Taoufik, N. et al. The state of art on the prediction of efficiency and modeling of the processes of pollutants removal based on machine learning. Sci. Total. Environ. 807, 150554. https://doi.org/10.1016/j.scitotenv.2021.150554 (2022).Article
ADS
CAS
PubMed
Google Scholar
Lijian, L. et al. Machine learning predicting wastewater properties of the aqueous phase derived from hydrothermal treatment of biomass. Bioresour. Technol. 258, 127348. https:/doi.org/10.1016/j.biortech.2022.127348 (2022).Article
CAS
Google Scholar
Yang, Y. et al. Artificial intelligence-enabled detection and assessment of Parkinson’s disease using nocturnal breathing signals. Nat. Med. 28, 2207–2215. https://doi.org/10.1038/s41591-022-01932-x (2022).Article
CAS
PubMed
PubMed Central
Google Scholar
Nguyen, N. & Guo, Y. Comparisons of sequence labeling algorithms and extensions. J. ACM. https://doi.org/10.1145/1273496.1273582(2007).Article
Google Scholar
Kiangala, S. K. & Wang, Z. An effective adaptive customization framework for small manufacturing plants using extreme gradient boosting-XGBoost and random forest ensemble learning algorithms in an Industry 4.0 environment. Mach. Vis. Appl. 4, 100024. https://doi.org/10.1016/j.mlwa.2021.100024 (2021).Article
Google Scholar
Jang, H. M., Yoo, S., Choi, Y.-K., Park, S. & Kan, E. Adsorption isotherm, kinetic modeling and mechanism of tetracyclineon Pinus taeda-derived activated biochar. Bioresour. Technol. 259, 24–31. https://doi.org/10.1016/j.biortech.2018.03.013 (2018).Article
CAS
PubMed
Google Scholar
Wang, H. et al. Sorption of tetracycline on biochar derived from rice straw and swine manure. RSC. Adv. 8, 16260–16268. https://doi.org/10.1039/C8RA01454J (2018).Article
ADS
CAS
PubMed
PubMed Central
Google Scholar
Zhang, P., Li, Y., Cao, Y. & Han, L. Charcteristics of tetracycline adsorption by cow manure biochar prepared at different pyrolysis temperatures. Bioresour. Technol. 285, 121348. https://doi.org/10.1016/j.biortech.2019.121348 (2019).Article
CAS
PubMed
Google Scholar
Choi, Y.-K. et al. Adsorption behavior of tetracycline onto Spirulina sp. (microalgae)-derived biochars produced at different temperatures. Sci. Total. Environ. 710, 136282. https://doi.org/10.1016/j.scitotenv.2019.136282 (2020).Article
ADS
CAS
PubMed
Google Scholar
Kim, J.E. et al. Adsorptive removal of tetracycline from aqueous solution by maple leaf-derived biochar. Bioresour. Technol. 306, 123092. https://doi.org/10.1016/j.biortech.2020.123092 (2020).Article
CAS
PubMed
Google Scholar
Shen, Q. et al. Removal of tetracycline from an aqueous solution using manganese dioxide modified biochar derived from Chinese herbal medicine residues. Environ. Res. 183, 109195. https://doi.org/10.106/j.envres.2020.109195 (2020).Article
CAS
PubMed
Google Scholar
Shisuo, F. et al. Preparation of tea residue biochar and its removal characteristics of tetracycline in solution. Environ. Sci. Technol. 41, 1308–1318 (2020).
Google Scholar
Xu, D. et al. Application of biochar derived from pyrolysis of waste fiberboard on tetracycline adsorption in aqueous solution. Front. Chem. 7, 943. https://doi.org/10.3389/fchem.2019.0094 (2020).Article
ADS
PubMed
PubMed Central
Google Scholar
Chen, Y. et al. Preparation of Eucomia ulmoides lignin-based high-performance biochar containing sulfonic group: Synergistic pyrolysis mechanism and tetracycline hydrochloride adsorption. Bioresour. Technol. 329, 124856. https://doi.org/10.1016/j.biortech.2021.124856 (2021).Article
CAS
PubMed
Google Scholar
Zheng, Z. et al. Preparation of mesoporous batatas biochar via soft-template method for high efficiency removal oftetracycline. Sci. Total. Environ. 787, 147397. https://doi.org/10.1016/j.scitotenv.2021.147397 (2021).Article
ADS
CAS
PubMed
Google Scholar
Wilschut, R. A. et al. Combined effects of warming and drought on plant biomass depend on plant woodiness and community type: A meta-analysis. Proc. R. Soc. B 289, 20221178. https://doi.org/10.1098/rspb.2022.1178 (2022).Article
PubMed
PubMed Central
Google Scholar
Zhu, X., Wang, X. & Ok, Y. S. The application of machine learning methods for prediction of metal sorption onto biochars. J. Hazard. Mater. 378, 12077. https://doi.org/10.1016/j.jhazmat.2019.06.004 (2019).Article
CAS
PubMed
Google Scholar
Oh, S. et al. Effects of biochar addition on the fate of ciprofloxacin and its associated antibiotic tolerance in an activated sludge microbiome. Environ. Pollut. 306, 119407. https://doi.org/10.1016/j.envpol.2022.119407 (2022).Article
CAS
PubMed
Google cholar
Matloob, F. et al. Software defect prediction using ensemble learning: A systematic literature review. IEEE. Access 9, 98754–98771. https://doi.org/10.1109/ACCESS.2021.3095559 (2021).Article
Google Scholar
Speiser, J. L., Miller, M. E., Tooze, J. & Ip, E. A comparison of random forest variable selection methods for classification prediction modeling. Expert Syst. Appl. 134, 93–101. https://doi.org/10.1016/j.solener.2017.04.066 (2019).Article
PubMed
PubMed Central
Google Scholar
Persson, C., Bacher, P., Shiga, T. & Madsen, H. Multi-site solar power forecasting using gradien boosted regression trees. Sol. Energy 150, 423–436. https://doi.org/10.1016/j.solener.2017.04.066 (2017).Article
ADS
Google Scholar
Alabdrabalnabi, A., Gautam, R. & Sarathy, S. M. Machine learning to predict biochar and bio-oil yields from co-pyrolysis of biomass and plastics. Fuel 328, 125303. https://doi.org/10.1016/j.fuel.022.125303 (2022).Article
CAS
Google Scholar
Chen, X. et al. Adsorption of copper and zinc by biochars produced from pyrolysis of hardwood and corn straw in aqueous solution. Bioresour. Technol. 102, 8877–8884. https://doi.org/10.1016/j.biortech.2011.06.078 (2011).Article
CAS
PubMed
Google Scholar
Chen, H. Biogenic silica nanoparticles derived from rice husk biomass and their applications. Ceram. Int. 41, 275–281 (2013).CAS
Google Scholar
Yu, J. et al. Influence of temperature and particle size on structural characteristics of chars from Beechwood pyrolysis. J. Anal. Appl. Pyrol. 130, 127–134. https://doiorg/10.1016/j.jaap.2018.01.018 (2018).Article
CAS
Google Scholar
Windeatt, J. H. et al. Characteristics of biochars from crop residues: Potential for carbon sequestration and soil amendment. J. Environ. Manag. 146, 189–197. https://doi.org/10.1016/j.jenvman.2014.08.003 (2014).Article
CAS
Google Scholar
Fang, Y. et al. Concentrated solar thermochemical gasification of biomass: Principles, applications, and development. Renew. Sustain. Energy Rev. 150, 111484. https://doi.org/10.1016/j.rser.2021.111484 (2021).Article
CAS
Google Scholar
Eduah, J. O. et al. Nonlinear sorption of phosphorus onto plant biomass-derived biochars at different pyrolysis temperatures. Environ. Technol. Innov. 19, 100808. https://doi.org/10.1016/j.eti.2020.100808 (2020).Article
Google Scholar
Ahmad, M. et al. Trichloroethylene adsorption by pine needle biochars produced at various pyrolysis temperatures. Bioresour. Technol. 143, 615–622. https://doi.org/10.1016/jbiortech.2013.06.033 (2013).Article
CAS
PubMed
Google Scholar
Usman, A. R. et al. Biochar production from date palm waste: Charring temperature induced changes in composition and surface chemistry. J. Anal. Appl. Pyrol. 115, 392–400. https://doi.org/10.1016/j.jaa.2015.08.016 (2015).Article
CAS
Google Scholar
Hamza, M. A. et al. Gaussian process regression and machine learning methods for carbon-based material adsorption. Adsorp. Sci. Technol. 2022, 3901608. https://doi.org/10.1155/2022/3901608 (2022).Article
CAS
Google Scholar
Kang, L.-L. et al. Removal of pollutants from wastewater using coffee waste as adsorbent: A review. J. Water Process. Eng. 49, 103178. https://doi.org/10.1016/j.jwpe.2022.103178 (2022).Article
Google Scholar
Download referencesAcknowledgementsThe authors thank all our colleagues and students involved in this study for their unremitting efforts. They thank Editage for the providing language help. They also thnk anonymous reviewers for their comments and suggestions.FundingThis work was supported by the President’s Foundation of Tarim University (TDZKSS202152; TDZKSS202156; TDZKCQ202002), Bingtuan Science and Technology Program (2021DB019; 2022CB001-01). Natioal Natural Science Foundation of China (42275014).Author informationAuthors and AffiliationsKey Laboratory of Tarim Oasis Agriculture (Tarim University), Ministry of Education, Xinjiang, 843300, ChinaPengyan Zhang, Chong Liu, Dongqing Lao, Xiaoyan Qian, Xuefei Hu, Yongjun You & Fayong LiCollege of Water Resources and Architectural Engineering, Tarim University, Xinjiang, 843300, ChinaPengyan Zhang, Chong Liu, Xiaoyan Qian, Xuefei Hu, Yongjun You & Fayong LiCollege of Information Engineering, Tarim University, Xinjiang, 843300, ChinaDongqing LaoInstitution of Research and Development, Duy Tan University, Da Nang, 550000, VietnamXuan Cuong NguyenDepartment of Biotechnology and Medical Engineering, National Institute of TechnologyRourkela, Odisha, 769008, IndiaBalasubramanian ParamasivanDepartment of Physical Sciences, Industrial Chemistry Programme, Landmark University, Omu-Aran, Kwara State, NigeriaAdejumoke Abosede InyinborAuthorsPengyan ZhangView author publicationsYou can alo search for this author in
PubMed Google ScholarChong LiuView author publicationsYou can also search for this author in
PubMed Google ScholarDongqing LaoView author publicationsYou can also search for this author in
PubMed Google ScholarXuan Cuong NguyenView author publicationsYou can also search for this author in
PubMed Google ScholarBalasubramanian ParamasivanView author publicationsYou can also search for this author in
PubMed Google ScholarXiaoyan QianView author publicationsYou can also search for this author in
PubMed Google ScholarAdejumoke Abosede InyinborView author publicationsYou can also search for this author in
PubMed Google ScholarXuefei HuView author publicationsYou can also search for this uthor in
PubMed Google ScholarYongjun YouView author publicationsYou can also search for this author in
PubMed Google ScholarFayong LiView author publicationsYou can also search for this author in
PubMd Google ScholarContributionsP.Z.: experiments, data analysis, writing. C.L.: experiments, data analysis, writing. D.L.: data analysis, project funding. X.Q., X.H. and Y.Y.: writing, reviewing and editing. L.F.: project management, funding acquisition, writing, review and editing. X.C.N., B.P. and A.A.I.: review and guidance.Corresponding authorCorrespondence to
Dongqing Lao.Ethics declarations
Competing interests
The authors declare no competing interests.
Additional informationPublisher’s noteSpringer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Supplementary InformationSupplementary Information.Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 nternational 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 he 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/.
Reprints and PermissionsAbout this articleCite this articleZhang, P., Liu, C., Lao, D. et al. Unveiling the drives behind tetracycline adsorption capacity with biochar through machine learning.
Sci Rep 13, 11512 (2023). htts://doi.org/10.1038/s41598-023-38579-8Download citationReceived: 05 April 2023Accepted: 11 July 2023Published: 17 July 2023DOI: https://doi.org/10.1038/s41598-023-38579-8Share tis articleAnyone you share the following link with will be able to read this content:Get shareable linkSorry, a shareable link is not currently available for this article.Copy to clipboard
Provided by the Springer Nature SharedIt content-sharing initiative
CommentsBy 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.