Activated carbon is often used in the drinking water advanced treatment process, which has good antibiotic removal capacity. However, the presence of nanoplastics (NPs) as carriers may increase the risk of antibiotic leakage in the carbon filtration column. We designed experiments with the polystyrene nanoplastics (PSNPs) concentration, norfloxacin (NOR) concentration, flow rate, and ionic strength as four orthogonal factors to investigate the effects of each factor on NOR removal by carbon filtration columns. The influence mechanism of PSNPs was inferred by combining with NOR transport curves and characterization analysis, and a prediction model of NOR removal efficiency was established through back-propagation (BP) network. The results showed that the increase of both PSNPs concentration and flow rate decreased the NOR removal efficiency. There was an optimal value of NOR concentration to maximize the NOR removal efficiency, while with increasing ionic strength, the NOR removal efficiency decreased, then increased, and finally decreased again in an inverted ‘N’ pattern. Furthermore, PSNPs can affect NOR removal efficiency via carrier function and aggregation on the activated carbon surface. On the other hand, the relative errors of the predicted and experimental values for two evaluated samples were 3.37 and 6.62%, respectively, indicating a good prediction effect.

  • Orthogonal experiments were designed to investigate the NOR removal efficiency by a carbon filtration column under different conditions.

  • The effect of each orthogonal factor on the NOR removal efficiency by a carbon filtration column was discussed and the role played by PSNPs was inferred.

  • A prediction model for the removal effectiveness of NOR by the carbon filtration column in the presence of PSNPs was developed.

Plastic products are widely used in various aspects of human life and social production, due to their low production cost and excellent properties, such as shopping bags, agricultural films, and medical components (Zhang et al. 2022). According to statistics, global plastic production has increased dramatically from 1.7 million tons in the 1950s to 368 million tons in 2019 (Yee et al. 2021). The huge plastic production has led to more than 13 million tons of plastic waste being discharged into rivers and oceans every year, and this quantity will continue to rise (Enfrin et al. 2019). Due to the difficulty in the degradability of materials, these plastic wastes usually accumulate in water bodies for decades. The larger plastic blocks are subjected to various physical, chemical, and biological effects, gradually breaking down into microplastics (MPs) with particle sizes less than 5 mm and then into nanoplastics (NPs) with particle sizes less than 1 μm (Peller et al. 2022). In comparison to MPs, NPs have smaller particle sizes and larger specific surface areas, which gives them higher surface energy and facilitates the adsorption of heavy metal ions and organic pollutants (Feng et al. 2022). Wang et al. (2019) investigated the adsorption capacity of polystyrene (PS) from 170 μm to 50 nm for phenanthrene and nitrobenzene. The results showed that the adsorption coefficients of both phenanthrene and nitrobenzene increased with decreasing PS particle size and reached the maximum at 235 nm. It indicated that NPs have stronger adsorption capacity than MPs and are more likely to be potential transport carriers of organic pollutants in water bodies. Pivokonsky et al. (2018) found considerable levels of particles smaller than 1 μm in size in raw water samples, with abundances ranging from 111 to 2,181 L−1 particles. Although such small particles could not be analyzed for proper material composition, it can be expected that an unknown percentage of them could be of plastic origin. In addition, once entering the organism, NPs may release contaminants that lead to the Trojan horse effect with increased toxicity (Binelli et al. 2017). It indicates that the adsorption capacity of NPs for contaminants and the potential mobility of NPs can influence the transport behavior of contaminants in the aquatic environment, which may pose a threat to drinking water quality safety.

Antibiotics are metabolites produced by bacteria, molds, and other microorganisms or analogs synthesized by artificial means (Inyinbor et al. 2018). They are widely used in human and veterinary medicine because of their excellent antibacterial effects (Kovalakova et al. 2020). While the use of antibiotics has brought many benefits to humans, excessive use has also caused a number of environmental problems. The most important of these problems is the generation of antibiotic-resistant genes (ARGs) and antibiotic-resistant bacteria (ARB), which negatively affect the effectiveness of antibiotics in treating human and animal pathogens. According to the chemical structures, antibiotics are mainly classified into β-lactams, macrolides, quinolones, sulfonamides, and tetracyclines (Kümmerer 2009). Among them, quinolone antibiotics are commonly used to treat infectious diseases such as urinary tract infections, respiratory tract infections, and gastrointestinal and abdominal infections caused by Escherichia coli or other gram-negative pathogens (Naeem et al. 2016). In China, the scale of production and consumption of quinolone antibiotics is huge, taking norfloxacin (NOR) as an example, the use of NOR in the whole country was up to 5,440 tons in 2013 (Zhang et al. 2015). On the other hand, since antibiotics are not completely absorbed by humans and animals, for example, 70% of ofloxacin and 30% of NOR are excreted unchanged in urine (Zhang et al. 2013), the residual antibiotics may be discharged through the municipal sewerage network into municipal wastewater treatment plants and eventually into the natural water bodies. However, the conventional drinking water treatment processes are also not effective in removing antibiotics and have been reported to remove less than 30% of NOR (He et al. 2015). It would lead to the increasing resistance of bacteria to quinolone antibiotics and the increased abundance of ARGs and ARB has been demonstrated in water distribution pipeline systems (Xu et al. 2016), which makes humans and animals more vulnerable to these microorganisms. Therefore, it is particularly important to control the content of quinolone antibiotics in drinking water treatment.

Filtration is one of the core processes in drinking water treatment, and the quality of the filtrated water is affected by various aspects such as filtration methods, pre-filtration water quality, and filtration media. For example, Liang et al. (2021) used ultrafiltration and two-stage reverse osmosis-integrated membrane filtration technology to treat pig wastewater, and the results showed that integrated membrane filtration could reduce more than 99.0% of conventional pollutants and more than 99.79% of ARGs, but it was not suitable for municipal drinking water treatment due to its high cost. Recent research also reported that sand filtration was superior to other conventional drinking water processes in terms of absolute abundance of ARGs removal (Tan et al. 2019), but its removal of antibiotics was not satisfactory, and furthermore, the filtered water did not directly meet the requirements of drinking water. Granular activated carbon (GAC) is an adsorbent material and filter material with a well-developed pore structure, high surface activity, and good mechanical and chemical stability, which is commonly used in the advanced treatment process of drinking water treatment to adsorb and remove fine particulate matter and antibiotics from water. The removal of quinolones by batch adsorption using GAC alone can achieve removal rates of more than 90% (Liu et al. 2017). However, when the water coexists with NPs, it has a strong adsorption effect on antibiotics due to its large specific surface area, which may increase the risk of antibiotic migration and diffusion in drinking water, and moreover, the static adsorption effect does not represent a treatment effect of the carbon filtration. The removal efficiency of carbon filter columns for quinolones in the presence of NPs was hardly reported, and the effect of NPs on the removal of quinolones by carbon filtration columns needs further study.

In this study, polystyrene nanoplastics (PSNPs) and NOR were selected as experimental materials and the orthogonal design method was used to design the experiment. The effects of four orthogonal factors, namely, PSNPs concentration, NOR concentration, filtration rate, and ionic strength, on NOR removal by the carbon filtration column were investigated. The possible influence mechanism was explored with characterization and transport curves. In addition, an orthogonal design and back-propagation (BP) network model were used to establish a prediction model for the removal effect of carbon filtration column on NOR under the influence of PSNPs. The model can predict the removal effect of carbon porous media on NOR under different conditions, and this model can greatly reduce the experimental volumes. The results contribute to further understand the adsorption characteristics of PSNPs on NOR on the surface of activated carbon and provide new ideas to predict the removal efficiency of quinolone antibiotics in the porous media of activated carbon when NPs and quinolones are coexisting. It also provides a theoretical basis for environmental and risk management of NPs and quinolone antibiotics.

Chemicals

PSNP microspheres were purchased from Big Goose Technology Co., LTD. (Tianjin, China). Using a transmission electron microscopy (TEM), the diameter of the PSNP microspheres was determined to be 98 ± 9 nm, with a density of 1.05 g/cm3 at 20 °C, and a specific surface area of 67 m2/g. Stock suspensions of 10 g/L were prepared from the original solution, and diluted with ultrapure water. PS stock suspension was sonicated for 3–5 min with a sonication bath before each extraction. The stock solution was stored in a dark place at a constant temperature of 4 °C and used for further experiments.

GAC of diameter 1.0–2.0 mm and pore volume of 0.357 cm3/g was purchased from Henan Huanyu Activated Carbon Co., LTD. (Henan, China). Then, the GAC was cleaned thrice by ultrasound for 5 min, and then rinsed with pure water for 3–5 times until the supernatant was clear.

NOR (98%) was purchased from Shanghai Aladdin Biochemical Technology Co., Ltd. (Shanghai China). NOR was dissolved in 4% sodium hydroxide solution, and then the NOR storage solution of 10 g/L was prepared with ultrapure water and stored in dark at 4 °C.

Characterization of GAC and detection of NOR

Scanning electron microscopy (SEM, FEI-F50) was used to observe the surface morphology of GAC before and after the experiment to help explain the mechanism of NOR removal in the carbon filtration column. The pore characteristics of the GAC were obtained by recording the adsorption isotherms of N2 on ASAP 2460 (Micrometrics, USA) at the boiling temperature of liquid nitrogen (77 K). The Brunauer–Emmett–Teller (BET) method was used to calculate the surface area and pore volume of GAC samples. The NOR concentration was measured by UV-5200 spectrophotometer produced by Shanghai Yuanxi Instrument Co., Ltd. (Shanghai, China). Standard curves are shown in Supplementary Material, Figure S1.

Orthogonal experimental design

Orthogonal experiments are able to reflect the influence of various factors on the removal of NOR from carbon filtration columns more comprehensively, and can greatly reduce the number of experiments that would otherwise be required. In this study, PSNPs concentration, NOR concentration, flow rate, and ionic strength were selected as orthogonal design factors, and four levels were determined for each factor. The orthogonal factor levels are shown in Table 1. Standard orthogonal tables L16(54) were used to design the experiments according to the determined four factors and four levels. Also, the concentration of NaCl was used uniformly to control the ionic strength in each experiment.

Table 1

Orthogonal experimental factors and levels

LevelsOrthogonal factors
NPNPs (mg/L)NOR (mg/L)Flow rate (mL/min)Ionic strength (mmol/L)
10 10 
10 15 10 15 
LevelsOrthogonal factors
NPNPs (mg/L)NOR (mg/L)Flow rate (mL/min)Ionic strength (mmol/L)
10 10 
10 15 10 15 

Adsorption and filtration experiments

Isothermal adsorption experiment of NOR by PSNPs

The mixture solutions containing 3 mg of PSNPs and different NOR concentrations (2, 4, 6, 8, 10, 12, 14, and 16 mg/L) were prepared in 50-mL sample bottles in sequence. The sample bottles were placed in a thermostatic shaker (150 rpm, 25 °C) and shaken for 16 h against light before sampling. Then, the sampled solutions were filtered through a pinhole filter membrane (0.22 μm) in time, and the remaining concentrations of NOR were determined by visible-UV spectrophotometer. Blank samples without PSNPs and two parallel samples were set up for each group of experiments.

NOR transport experiments under the influence of PSNPs

The experimental setup consists of magnetic agitator, peristaltic pump, carbon filtration column, and collector. Each part is connected by silicone tubes, and the carbon filtration column used in the experiment was made of Plexiglas acrylic column with an inner diameter of 15 mm and a height of 80 mm. The experimental setup sketch and more details of the setup are described in Supplemental Material, Figure S2. Before each experiment, pure water was pumped through the carbon filtration column using a peristaltic pump at a certain flow rate for 6 h to rinse impurities out of the column to reduce errors in NOR measurement. Also, prior to the transport experiments, 10 PV of background salt solution (5 mmol/L NaCl solution) was used to stabilize the experimental conditions, and the mixture of PSNPs (5 and 10 mg/L) and NOR (5 and 10 mg/L) was first equilibrated by dark adsorption in a constant temperature water bath shaker for 16 h. The experimental conditions are shown in Table 2. The transport experiments were performed by injecting 5 PV of prepared NOR solution or mixed solution of PSMPs and NOR, followed by 5 PV of a background solution containing the same ionic strength. The concentration of NOR in the effluent was measured for every 0.5 PV. A magnetic stirrer was used to continuously stir the mixed solution while the experiment was in progress to ensure that the PSNPs were uniformly dispersed in the solution. All solutions were injected into the column at 3 mL/min with empty bed contact time (EBCT) of 4.7 min and a hydraulic loading rate of 1.70 mL/min. Since this experiment is a dynamic process and the removal efficiency of GAC on NOR is high, the influence of PSNPs on the transport of NOR cannot be visualized by the outflow ratio of NOR, so the transport curve of NOR in the activated carbon filter column was plotted using the discharge concentration of NOR as the vertical coordinate and PV as the horizontal coordinate.

Table 2

Summary of experimental conditions for NOR transport

ColumnMaterialsBackground solutionColumn lengthFlow
(mm)(mL/min)
5 mg/L NOR 5 mM NaCl 80 
10 mg/L NOR 5 mM NaCl 80 
5 mg/L PSNPs 5 mg/L NOR 5 mM NaCl 80 
5 mg/L PSNPs 10 mg/L NOR 5 mM NaCl 80 
10 mg/L PSNPs 5 mg/L NOR 5 mM NaCl 80 
10 mg/L PSNPs 10 mg/L NOR 5 mM NaCl 80 
ColumnMaterialsBackground solutionColumn lengthFlow
(mm)(mL/min)
5 mg/L NOR 5 mM NaCl 80 
10 mg/L NOR 5 mM NaCl 80 
5 mg/L PSNPs 5 mg/L NOR 5 mM NaCl 80 
5 mg/L PSNPs 10 mg/L NOR 5 mM NaCl 80 
10 mg/L PSNPs 5 mg/L NOR 5 mM NaCl 80 
10 mg/L PSNPs 10 mg/L NOR 5 mM NaCl 80 

NOR removal experiments by a carbon filtration column

The NOR removal experiments continued using the setup shown in Supplemental Material, Figure S2. When the flow rate was 1, 3, 5, and 10 mL/min, the EBCT of the carbon filter column was 14.1, 4.7, 2.8, and 1.4 min and the hydraulic load was 0.56, 1.70, 2.83, and 5.66 cm/min at flow rates of 1, 3, 5, and 10 mL/min, respectively. Then, as in the preparation steps of the previous NOR transport experiment, the carbon filter column was rinsed with a certain flow rate of pure water for 6 h to reduce the measurement error of NOR by removing impurities. After that, 10 PV of the corresponding salt solution was used to stabilize the experimental conditions. The mixture solution of PSNPs and NOR was prepared in advance according to Table 3 and dark adsorbed in a shaker for 16 h before being taken out and added to the corresponding salt solution. A magnetic stirrer was used to stir the mixed solution to prevent the aggregation of PSNPs and the mixed solution was injected into the carbon filtration column with a peristaltic pump at the flow rate of each group of experiments. The absorbance of the effluent solution was measured every 10 min until the activated carbon adsorption reaches relative saturation (i.e., little change in the absorbance of the effluent water). The NOR removal rate at the relative saturation of the carbon filter column was calculated based on the change of NOR concentration, and the calculation formula was as the following:
(1)
where is the removal rate of NOR at relative saturation of the carbon filter column; is the initial NOR concentration, mg/L; is the concentration of NOR in the effluent water when the carbon filter column is relatively saturated, mg/L. All experimental conditions and results are summarized in Table 3. The groups 1–16 are the orthogonal experiments required to build the prediction model, and the groups 17–18 are the experiments to test the prediction performance of the model.
Table 3

Summary of experimental conditions for NOR removal

Experimental groupsNPNPsNORFlow rateIonic strengthNOR removal rate (nr)
(mg/L)(mg/L)(mL/min)(mmol/L)(%)
78.12 
79.55 
10 10 63.92 
15 10 15 41.43 
10 72.10 
15 87.06 
10 10 37.24 
15 44.05 
15 37.46 
10 10 10 75.41 
11 10 58.98 
12 15 45.68 
13 10 10 8.25 
14 10 42.41 
15 10 10 15 52.83 
16 10 15 10 70.58 
17 10 85.23 
18 56.22 
Experimental groupsNPNPsNORFlow rateIonic strengthNOR removal rate (nr)
(mg/L)(mg/L)(mL/min)(mmol/L)(%)
78.12 
79.55 
10 10 63.92 
15 10 15 41.43 
10 72.10 
15 87.06 
10 10 37.24 
15 44.05 
15 37.46 
10 10 10 75.41 
11 10 58.98 
12 15 45.68 
13 10 10 8.25 
14 10 42.41 
15 10 10 15 52.83 
16 10 15 10 70.58 
17 10 85.23 
18 56.22 

Construction of a BP neural network

The PSNPs concentration, NOR concentration, flow rate, and ionic strength were selected as input variables, and the NOR removal rate at the relative saturation of the carbon filter column was selected as the output variable. Experimental groups 1–16 in Table 3 were used as training samples to train the neural network, and groups 17–18 were used as prediction samples to test the reliability of the prediction model.

Firstly, the samples are normalized using the mapminmax function, and a three-layer neural network model is constructed. The transfer functions of the implicit layer and the output layer are Tan-sig and Purelin functions, respectively. Then, the learning objective function, learning rate, and momentum term coefficients are set to be 1 × 10−8, 0.1, and 0.95, respectively, and the number of iterations is set at 1,000. On the other hand, the number of input nodes and the number of output nodes of the network topology are determined by the number of input and output variables, respectively, so that the topology is determined to contain four input nodes and one output node. In addition, there is no explicit method to determine the number of nodes in the hidden layer, and the range of the number of nodes in the hidden layer is generally determined according to the following empirical formula:
(2)
where J is the number of nodes in the hidden layer; n is the number of nodes in the input layer; m is the number of nodes in the output layer; a ranges from 1 to 10. The number of hidden layers could be determined from 4 to 12 by the empirical formula, and then the number of nodes in the hidden layer was finally determined according to the trial-and-error method. The results are shown in Table 4, the training error was minimized when the number of nodes in the hidden layer was 5. Therefore, the neural network topology was finally determined to be 4 × 5 × 1, and its structure schematic is shown in Figure 1. The BP prediction model is established by MATLAB software, and the related code is shown in Supplementary Material, Figure S3.
Table 4

The relation between the hidden layer node number and training error

The number of nodes in the hidden layer456789101112
Training error 0.272 0.128 0.166 0.395 0.185 0.184 0.248 0.246 0.368 
The number of nodes in the hidden layer456789101112
Training error 0.272 0.128 0.166 0.395 0.185 0.184 0.248 0.246 0.368 
Figure 1

The sketch of prediction model network topology.

Figure 1

The sketch of prediction model network topology.

Close modal

Adsorption isotherms

The Langmuir and Freundlich models were used to fit the experimental data for the adsorption isotherms of NOR by PSNPs, and the results and related parameters are shown in Figure 2 and Table 5. At fixed PSNPs concentration, the equilibrium adsorption of NOR by PSNPs increased with increasing NOR concentration, but the adsorption rate gradually became slower. It was because the adsorption sites on the surface of PSNPs were not fully utilized at low NOR concentrations, and when the NOR concentration increased, the adsorption sites on the surface of PSNPs gradually reached saturation, resulting in the adsorption rate decreasing gradually. It can also be seen from Figure 2 and Table 5 that the Langmuir model can better fit the adsorption process of PSNPs to NOR with the correlation coefficient (R2) of 0.972, while the R2 fitted by the Freundlich isothermal model is 0.83. It indicated that the adsorption of PSNPs to NOR is a single molecular layer adsorption on a relatively homogeneous surface. The results of the Langmuir model calculation showed that the maximum adsorption of NOR by PSNPs was 22.9 mg/g. In addition, the adsorption isotherms of NOR by PSNPs were more in line with the nonlinear fit, which indicated that the adsorption process of NOR by PSNPs was dominated by several mechanisms, such as microporous filling, hydrogen bonding, π–π interaction, and electrostatic interaction (Pei et al. 2004).
Table 5

Isothermal model fitting parameters for the adsorption of NOR by PSNPs

Langmuir
Freundlich
KL (L/mg)Qm (mg/g)R2KF ((mg/g)/(mg/L)1/n)nR2
0.585 22.9 0.972 11.46 1.20 0.83 
Langmuir
Freundlich
KL (L/mg)Qm (mg/g)R2KF ((mg/g)/(mg/L)1/n)nR2
0.585 22.9 0.972 11.46 1.20 0.83 
Figure 2

Adsorption isotherms of NOR by PSNPs.

Figure 2

Adsorption isotherms of NOR by PSNPs.

Close modal

Effect of PSNPs on the transport of NOR in carbon filter columns

As seen in Figure 3, the concentrations of NOR in the effluent of carbon filtration columns were higher when PSNPs were transported along with NOR than that of NOR alone, regardless of the NOR concentration of 5 or 10 mg/L. It indicated that the presence of PSNPs facilitated the transport of NOR in the carbon filtration column. In Figure 3(a), compared with NOR transport alone, the NOR concentrations in the effluent of carbon filtration columns increased by 48 and 72.4 μg/L when NOR was co-transported with 5 mg/L of PSNPs and 10 mg/L of PSNPs, respectively. Similarly, compared to 10 mg/L NOR solution alone, the NOR concentration in the carbon column effluent increased by 44.6 and 76.2 μg/L when NOR was co-transported with 5 mg/L of PSNPs and 10 mg/L of PSNPs, respectively, as shown in Figure 3(b). It indicated that when the NOR concentration was certain, the greater the concentration of PSNPs in the mixed solution, the greater the effect on NOR transport. In contrast to PSNPs, GAC has a higher adsorption capacity for NOR. When NOR was transported along with PSNPs, PSNPs occupied the available adsorption sites on the activated carbon, which is the reason for the higher NOR concentration in the effluent when PSNPs were transported along with NOR than NOR alone. Moreover, at the same NOR concentration, the more PSNPs were co-transported, the more adsorption sites were occupied by PSNPs on the activated carbon, resulting in higher NOR concentration in the effluent.
Figure 3

Effect of PSNPs on the transport of NOR in carbon filter columns: (a) 5 mg/L NOR and (b) 10 mg/L NOR.

Figure 3

Effect of PSNPs on the transport of NOR in carbon filter columns: (a) 5 mg/L NOR and (b) 10 mg/L NOR.

Close modal

Analysis of variance and removal mechanism

Analysis of variance was performed with SPSS software based on the experimental results in Table 3. The results are shown in Table 6, and it can be seen that the degree of influence of each factor on the retention of NOR by the carbon filtration column was flow rate, NOR concentration, ionic strength, and PSNPs concentration, in that order. In addition, flow rate was a significant influence on NOR removal by carbon filtration column (p ≤ 0.05), and there was no significant difference in the effect of different levels of other orthogonal factors on NOR removal by the carbon filtration column (i.e., p > 0.05).

Table 6

Analysis of variance

Source of differenceDiscrepancy sum of squares (SS)Degree freedomMean square (MS)F-statisticSignificance (p)
PSNPs concentration 1,093.796 364.599 3.717 0.152 
NOR concentration 1,255.039 588.680 2.687 0.127 
Flow rate 1,522.039 507.541 2.317 0.049a 
Ionic strength 1,293.987 431.329 1.969 0.131 
Source of differenceDiscrepancy sum of squares (SS)Degree freedomMean square (MS)F-statisticSignificance (p)
PSNPs concentration 1,093.796 364.599 3.717 0.152 
NOR concentration 1,255.039 588.680 2.687 0.127 
Flow rate 1,522.039 507.541 2.317 0.049a 
Ionic strength 1,293.987 431.329 1.969 0.131 

aWhen the F-statistic corresponds to a concomitant probability of less than 0.05 (p ≤ 0.05), it means that there is a 95% probability that this source of variation has a significant effect on NOR removal by carbon filtration column.

The mean value of NOR removal rate of each orthogonal factor at each level was used as an indicator to represent the degree of influence of the factor on NOR removal by the carbon filtration column at that level and the value of is calculated as the following:
(3)
where represents the mean value of NOR removal rate of a certain orthogonal factor at level i; represents the sum of removal rates of a certain orthogonal factor at level i in the orthogonal experiment ( = 1, 2, 3, 4). The values of each orthogonal factor at different levels were calculated according to the formula, and the results are shown in Figure 4.
Figure 4

The effect of each orthogonal factor on NOR removal by a carbon filtration column.

Figure 4

The effect of each orthogonal factor on NOR removal by a carbon filtration column.

Close modal

As shown in Figure 4, the lower the PSNPs concentration, the more favorable the removal of NOR by the carbon filtration column (i.e., k-value decreased as the PSNPs concentration increases from level 1 to level 4). On the one hand, when the carbon filter column is relatively saturated, some of the adsorption sites that should be occupied by NOR are occupied or obstructed by PSNPs, resulting in an increase in NOR effluent concentration. Also because of the concentration gradient, the higher the concentration of PSNPs, the more adsorption sites are occupied or obstructed, which is why the concentration of NOR leakage increases with PSNPs concentration when the carbon filter column is relatively saturated. On the other hand, a part of NOR might be adsorbed on the PSNPs and removed out of the carbon filtration column along with the water flow. This was confirmed in NOR transport experiments as well. The value of k gradually increased with the increase in the NOR concentration to the peak at level 2, and subsequently the value of k decreased with the further increase in the concentration. It may be due to the fact that in the carbon filtration column, the weights of the filled activated carbon were almost equal, and the numbers of their active sites were also close to each other. When the NOR concentration was low, most of the NOR could be adsorbed on the active sites of the activated carbon. But when the NOR concentration exceeded a certain limit, the NOR generated by the concentration exceeding the limit could not be retained in the carbon filtration column because the active sites were saturated, resulting in a decrease in the removal rate of NOR with the increase in the NOR concentration.

As shown in Figure 4, the value of k decreased as the flow rate increased. It may be at a low flow rate, the proportion of horizontal diffusion of the mixed solution was higher than that at a high flow rate, which increased the chance of NOR contacting with the activated carbon media, so more NOR was adsorbed on the activated carbon media. While at a high flow rate, the shear force was enhanced and the mixed solution passed through the carbon filtration column in a shorter time, which made it easier for NOR to penetrate the carbon filtration column. As shown in Figure 4, the value of k decreased, then increased, and again decreased with the increase of ionic strength, forming an inverted ‘N’ trend. According to previous studies, the penetration of PSNPs in the carbon filtration column decreased with the increase in the ionic strength (Ji et al. 2023), and the presence of Na+ would form a competitive adsorption relationship with NOR. Therefore, when the ionic strength increased from level 1 to level 2, the value of k decreased slightly, probably because the amount of NOR penetration due to competitive adsorption of Na+ was greater than the amount of NOR carried by PSNPs and trapped in the carbon filtration column. Also, when the ionic strength increased from level 2 to level 3, the PSNPs were further trapped on the surface of the activated carbon and the number of NORs carried by PSNPs and trapped was dominant, thus showing an increase in k-value. However, when the ionic strength was further increased from level 3 to level 4, the effect of the increase in the ionic strength on the retention of PSNPs by activated carbon was weakened, and the competitive adsorption of Na+ with NOR played a dominant role, which eventually led to the decrease in the k-value. The exact relationship of action needs to be determined by further experiments.

Characterization analysis

The performance of activated carbon for NOR adsorption was also influenced by its particle size, specific surface area, and pore structure (AlHazmi et al. 2022). According to the pore size, pore structure can be classified as microporous (<2 nm), mesoporous (2–50 nm), and macroporous (>50 nm) (Hassan et al. 2020). Among them, the effectiveness of carbon-based materials for antibiotic adsorption is mainly influenced by the number of micropores and mesopores (Nakagawa et al. 2004). When the size of NOR molecules and the pore size of activated carbon were close to each other, the relative pore walls would have a ‘pore filling effect’, which would increase the adsorption of NOR (Ismadji & Bhatia 2001). In order to investigate the pore structure of activated carbon, N2 adsorption–desorption experiments were performed at 77 K on GAC samples. As shown in Figure 5(a), the isotherms of N2 adsorption and desorption of GAC were type I isotherms, which indicated that GAC adsorbs mainly through narrow pores, and its adsorption capacity was mainly controlled by the pore volume (Ji et al. 2023). The BET equation was also used to calculate the specific surface area and pore volume of GAC. The results showed that the specific surface area (SBET) of GAC was estimated to be about 759 m2/g, the total pore volume was calculated to be 0.357 cm3/g, and the average pore size was measured to be 3.09 nm, indicating that GAC had a large specific surface area and a well-developed pore structure. The pore size distribution of GAC was also shown in Figure 5(b), and it could be found that the pore structure of GAC was overwhelmingly in the form of micropores and mesopores, which was favorable for NOR adsorption.
Figure 5

(a) N2 adsorption–desorption isotherms for GAC at 77 K and (b) distribution of GAC pore size.

Figure 5

(a) N2 adsorption–desorption isotherms for GAC at 77 K and (b) distribution of GAC pore size.

Close modal
In order to further investigate the reasons for the decrease in NOR removal when PSNPs and NOR passed through the relatively saturated carbon filtration column synergistically. The activated carbon samples before and after the dynamic adsorption filtration experiment were analyzed by SEM. As shown in Figure 6(a), the SEM image showed the surface morphology of the activated carbon before the experiment. It can be seen from the image that the activated carbon surface had a large number of pores with different sizes and shapes, and these complex pore structures provided additional adsorption areas for PSNPs and NOR. However, due to the pore size, PSNPs can only be adsorbed in the macroporous structure of activated carbon, while NOR can be adsorbed in macroporous, mesoporous, and microporous structures. Meanwhile, the SEM image of the activated carbon after the experiment in Figure 6(b) shows that most of the PSNPs were adsorbed on the surface of the activated carbon, and a small portion formed aggregates that were trapped by the pores and were harder to be released from them. Therefore, it was speculated that the presence of PSNPs might hinder the binding of NOR to more micropore and mesoporous structures inside the activated carbon, which resulted in a decrease of NOR removal. On the other hand, the increase in ionic strength promoted the aggregation between the PSNPs particles to combine with the activated carbon pores to form new microporous or mesoporous structures as shown in Figure 6(b), which promoted the removal of NOR in the activated carbon. Also, the increased Na+ concentration would form a competitive adsorption relationship with NOR. Therefore, it was speculated that as the ionic strength increased, it showed an increase in k-value when the aggregation of PSNPs dominated the NOR removal, while a decrease in k-value when the competitive adsorption of Na+ dominated the NOR removal.
Figure 6

(a) Surface morphology of activated carbon before the dynamic adsorption filtration experiment and (b) surface morphology of activated carbon after the dynamic adsorption filtration experiment.

Figure 6

(a) Surface morphology of activated carbon before the dynamic adsorption filtration experiment and (b) surface morphology of activated carbon after the dynamic adsorption filtration experiment.

Close modal

BP prediction model

After the BP neural network model was established, MATLAB software was used for to program and the train training function was called to train the network model, and the training results are shown in Figure 7. From Figure 7(a), it can be seen that the calculated values of the model for the training set samples basically match with their true values, and the correlation coefficient is 0.9661. The model fit of the overall training samples of the neural network (training set, test set, and validation set) is shown in Figure 7(b) with a correlation coefficient of 0.9643, which indicated that the predicted values of the training samples were closer to the real ones and the model fit is good. Figure 7(c) shows the training process curve of the BP neural network, and the model predicted optimally after four iterations with a training error of 0.0046. In addition, it can also be seen from Figure 7(d) that for the sample points used for training, the errors between the model predictions and the experimental values were small, which were almost in the range of ±0.1.
Figure 7

(a) The model fitting degree of neural network training set; (b) the model fitting degree of neural network all training samples; (c) the training process curve of BP neural network; and (d) the error comparison between the predicted value and expected value of training samples.

Figure 7

(a) The model fitting degree of neural network training set; (b) the model fitting degree of neural network all training samples; (c) the training process curve of BP neural network; and (d) the error comparison between the predicted value and expected value of training samples.

Close modal

In order to further test the reliability of the prediction model, the experimental conditions of groups 17 and 18 were set by random combination for each orthogonal factor and each level to test the prediction performance of the model. The sim function was called to make model predictions for the 17 and 18 groups of experiments and to inverse normalize the results. Table 7 shows the comparison of experimental and predicted values for groups 17 and 18. The predicted values of the model in groups 17 and 18 were not significantly different from the experimental values, and the relative errors were 3.32 and 6.22%, respectively, which indicated that the network model had a good prediction effect.

Table 7

Comparison of sample prediction results with experimental values

Experimental groupsLevels
Sample experimental values (%)Model predicted value (%)Relative error (%)
PSNPsNORFlow rateIonic strength
17 85.23 82.36 3.37 
18 56.22 52.44 6.62 
Experimental groupsLevels
Sample experimental values (%)Model predicted value (%)Relative error (%)
PSNPsNORFlow rateIonic strength
17 85.23 82.36 3.37 
18 56.22 52.44 6.62 

Activated carbon filtration process is one of the advanced common drinking water treatment processes used to remove antibiotics. Due to the environmental risks of quinolone antibiotics and the transport properties of NPs as carriers, it is very important to investigate the removal efficacy of quinolone antibiotics under the influence of NPs. In this study, the isothermal adsorption lines of NOR by PSNPs and the effect of PSNPs on the transport of NOR in the carbon filter column were investigated. The PSNPs concentration, NOR concentration, flow rate, and ionic strength were selected as four orthogonal factors to design orthogonal experiments. Analysis of variance using SPSS software was used to investigate the effect of each orthogonal factor on NOR removal from carbon filtration columns and to speculate on the possible mechanism. The transport curves of PSNPs synergistic NOR were combined with characterization to verify the possible mechanism of PSNPs in NOR removal from carbon filtration columns. In addition, a prediction model with the four orthogonal factors as input variables and the removal efficiency of NOR at the relative saturation of the carbon filter column as output variables was also developed by the BP neural network, and the model was evaluated by two sets of prediction samples.

The results showed that the adsorption curves of PSNPs to NOR were more consistent with the Langmuir model. The increase of PSNPs concentration might hinder the binding of NOR to the active site of activated carbon, and the PSNPs might also carry a portion of NOR through the filtration column, resulting in a decrease in the removal rate of NOR. The degree of influence on NOR removal from the carbon filtration column was in the order of flow rate, NOR concentration, ionic strength, and PSNPs concentration, and the flow rate was a significant influence factor. The effects of other orthogonal factors were speculated by the analysis of variance. High flow rates lead to faster passage of the mixed solution through the carbon filtration column compared to low flow rates, resulting in a lower probability of NOR contact with the active site. Also, there was an optimum concentration for the NOR concentration on the removal of NOR by a carbon filtration column, lower or higher than this concentration, the removal rate of NOR decreased. However, the increase of ionic strength would promote the aggregation of PSNPs with each other and the formation of new microporous or mesoporous structures with activated carbon, which would promote the removal rate of NOR; on the other hand, the increase of ionic strength would also strengthen the competitive adsorption of Na+ with NOR, which would decrease the removal rate of NOR. These two aspects work together to make the NOR removal rate show an inverse ‘N’ effect, i.e., the NOR removal efficiency of the carbon filtration column showed a trend of decreasing, then increasing, and then decreasing with the increase of ionic strength. Moreover, the overall model fit was 0.9643, and the relative errors between the model predicted and sample experimental values were 3.37 and 6.62%, respectively, for the two sets of predicted samples, which indicated that the model prediction was good. This also provided a theoretical calculation basis and reference for the removal of NPs and quinolone antibiotic NOR in drinking water treatment plants using the activated carbon filtration process.

H.J., Z.L., X.X., W.J., S.W., B.W., and X.X. contributed to the study conception and design. Material preparation, data collection, and analysis were performed by H.J. and W.J. The first draft of the manuscript was written by H.J. and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

This work is supported by the Key Laboratory of Poyang Lake Environment and Resource Utilization, Ministry of Education, Nanchang University (No. 2022Y04).

All relevant data are included in the paper or its Supplementary Information.

The authors declare there is no conflict.

AlHazmi
G. A. A.
,
AbouMelha
K. S.
,
El-Desouky
M. G.
&
El-Bindary
A. A.
2022
Effective adsorption of doxorubicin hydrochloride on zirconium metal-organic framework: equilibrium, kinetic and thermodynamic studies
.
Journal of Molecular Structure
1258
,
132679
.
https://doi.org/10.1016/j.molstruc.2022.132679
.
Binelli
A.
,
Del Giacco
L.
,
Santo
N.
,
Bini
L.
,
Magni
S.
,
Parolini
M.
,
Madaschi
L.
,
Ghilardi
A.
,
Maggioni
D.
,
Ascagni
M.
,
Armini
A.
,
Prosperi
L.
,
Landi
C.
,
La Porta
C.
&
Della Torre
C.
2017
Carbon nanopowder acts as a Trojan-horse for benzo(α)pyrene in Danio rerio embryos
.
Nanotoxicology
11
(
3
),
371
381
.
https://doi.org/10.1080/17435390.2017.1306130
.
Enfrin
M.
,
Dumée
L. F.
&
Lee
J.
2019
Nano/microplastics in water and wastewater treatment processes – origin, impact and potential solutions
.
Water Research
161
,
621
638
.
https://doi.org/10.1016/j.watres.2019.06.049
.
Feng
H.
,
Liu
Y.
,
Xu
Y.
,
Li
S.
,
Liu
X.
,
Dai
Y.
,
Zhao
J.
&
Yue
T.
2022
Benzo[a]pyrene and heavy metal ion adsorption on nanoplastics regulated by humic acid: cooperation/competition mechanisms revealed by molecular dynamics simulations
.
Journal of Hazardous Materials
424
,
127431
.
https://doi.org/10.1016/j.jhazmat.2021.127431
.
Hassan
N.
,
Shahat
A.
,
El-Didamony
A.
,
El-Desouky
M. G.
&
El-Bindary
A. A.
2020
Mesoporous iron oxide nano spheres for capturing organic dyes from water sources
.
Journal of Molecular Structure
1217
,
128361
.
https://doi.org/10.1016/j.molstruc.2020.128361
.
He
K.
,
Soares
A. D.
,
Adejumo
H.
,
McDiarmid
M.
,
Squibb
K.
&
Blaney
L.
2015
Detection of a wide variety of human and veterinary fluoroquinolone antibiotics in municipal wastewater and wastewater-impacted surface water
.
Journal of Pharmaceutical and Biomedical Analysis
106
,
136
143
.
https://doi.org/10.1016/j.jpba.2014.11.020
.
Inyinbor
A. A.
,
Bello
O. S.
,
Fadiji
A. E.
&
Inyinbor
H. E.
2018
Threats from antibiotics: a serious environmental concern
.
Journal of Environmental Chemical Engineering
6
(
1
),
784
793
.
https://doi.org/10.1016/j.jece.2017.12.056
.
Ismadji
S.
&
Bhatia
S. K.
2001
A modified pore-filling isotherm for liquid-phase adsorption in activated carbon
.
Langmuir
17
(
5
),
1488
1498
.
https://doi.org/10.1021/la0009339
.
Ji
H.
,
Liu
Z.
&
Jiang
W.
2023
Transport behavior of nanoplastics in activated carbon column
.
Environmental Science and Pollution Research.
30
,
26256
26269
.
https://doi.org/10.1007/s11356-022-24056-8
.
Kovalakova
P.
,
Cizmas
L.
,
McDonald
T. J.
,
Marsalek
B.
,
Feng
M.
&
Sharma
V. K.
2020
Occurrence and toxicity of antibiotics in the aquatic environment: a review
.
Chemosphere
251
,
126351
.
https://doi.org/10.1016/j.chemosphere.2020.126351
.
Kümmerer
K.
2009
Antibiotics in the aquatic environment – A review – part I
.
Chemosphere
75
(
4
),
417
434
.
https://doi.org/10.1016/j.chemosphere.2022.135943
.
Liang
C.
,
Wei
D.
,
Zhang
S.
,
Ren
Q.
,
Shi
J.
&
Liu
L.
2021
Removal of antibiotic resistance genes from swine wastewater by membrane filtration treatment
.
Ecotoxicology and Environmental Safety
210
,
111885
.
https://doi.org/10.1016/j.ecoenv.2020.111885
.
Liu
P.
,
Wang
Q.
,
Zheng
C.
&
He
C.
2017
Sorption of sulfadiazine, norfloxacin, metronidazole, and tetracycline by granular activated carbon: kinetics, mechanisms, and isotherms
.
Water, Air, & Soil Pollution
228
(
4
),
129
.
https://doi.org/10.1007/s11270-017-3320-x
.
Naeem
A.
,
Badshah
S. L.
,
Muska
M.
,
Ahmad
N.
&
Khan
K.
2016
The current case of quinolones: synthetic approaches and antibacterial activity
.
Molecules
21
(
4
).
https://doi.org/10.3390/molecules21040268.
Nakagawa
K.
,
Namba
A.
,
Mukai
S. R.
,
Tamon
H.
,
Ariyadejwanich
P.
&
Tanthapanichakoon
W.
2004
Adsorption of phenol and reactive dye from aqueous solution on activated carbons derived from solid wastes
.
Water Research
38
(
7
),
1791
1798
.
https://doi.org/10.1016/j.watres.2004.01.002
.
Pei
Z.
,
Yang
S.
,
Li
L.
,
Li
C.
&
Zhang
S.
2004
Effects of copper and aluminum on the adsorption of sulfathiazole and tylosin on peat and soil
.
Environmental Pollution
184
,
579
585
.
https://doi.org/10.1016/j.envpol.2013.09.038
.
Peller
J. R.
,
Mezyk
S. P.
,
Shidler
S.
,
Castleman
J.
,
Kaiser
S.
,
Faulkner
R. F.
,
Pilgrim
C. D.
,
Wilson
A.
,
Martens
S.
&
Horne
G. P.
2022
Facile nanoplastics formation from macro and microplastics in aqueous media
.
Environmental Pollution
313
,
120171
.
https://doi.org/10.1016/j.envpol.2022.120171
.
Pivokonsky
M.
,
Cermakova
L.
,
Novotna
K.
,
Peer
P.
,
Cajthaml
T.
&
Janda
V.
2018
Occurrence of microplastics in raw and treated drinking water
.
Science of The Total Environment
643
,
1644
1651
.
https://doi.org/10.1016/j.scitotenv.2018.08.102
.
Tan
Q.
,
Li
W.
,
Zhang
J.
,
Zhou
W.
,
Chen
J.
,
Li
Y.
&
Ma
J.
2019
Presence, dissemination and removal of antibiotic resistant bacteria and antibiotic resistance genes in urban drinking water system: a review
.
Frontiers of Environmental Science & Engineering
13
(
3
),
36
.
https://doi.org/10.1007/s11783-019-1120-9
.
Wang
J.
,
Liu
X.
,
Liu
G.
,
Zhang
Z.
,
Wu
H.
,
Cui
B.
,
Bai
J.
&
Zhang
W.
2019
Size effect of polystyrene microplastics on sorption of phenanthrene and nitrobenzene
.
Ecotoxicology and Environmental Safety
173
,
331
338
.
https://doi.org/10.1016/j.ecoenv.2019.02.037
.
Xu
L.
,
Ouyang
W.
,
Qian
Y.
,
Su
C.
,
Su
J.
&
Chen
H.
2016
High-throughput profiling of antibiotic resistance genes in drinking water treatment plants and distribution systems
.
Environmental Pollution
213
,
119
126
.
https://doi.org/10.1016/j.envpol.2016.02.013
.
Yee
M. S.
,
Hii
L.-W.
,
Looi
C. K.
,
Lim
W.-M.
,
Wong
S.-F.
,
Kok
Y.-Y.
,
Tan
B.-K.
,
Wong
C.-Y.
&
Leong
C.-O.
2021
Impact of microplastics and nanoplastics on human health
.
Nanomaterials
11
(
2
).
https://doi.org/10.3390/nano11020496.
Zhang
H.
,
Liu
P.
,
Feng
Y.
&
Yang
F.
2013
Fate of antibiotics during wastewater treatment and antibiotic distribution in the effluent-receiving waters of the Yellow Sea, northern China
.
Marine Pollution Bulletin
73
(
1
),
282
290
.
https://doi.org/10.1016/j.marpolbul.2013.05.007
.
Zhang
Q.-Q.
,
Ying
G.-G.
,
Pan
C.-G.
,
Liu
Y.-S.
&
Zhao
J.-L.
2015
Comprehensive evaluation of antibiotics emission and fate in the river basins of China: source analysis, multimedia modeling, and linkage to bacterial resistance
.
Environmental Science & Technology
49
(
11
),
6772
6782
.
https://doi.org/10.1021/acs.est.5b00729
.
Zhang
Y.
,
Li
Y.
,
Su
F.
,
Peng
L.
&
Liu
D.
2022
The life cycle of micro-nano plastics in domestic sewage
.
Science of The Total Environment
802
,
149658
.
https://doi.org/10.1016/j.scitotenv.2021.149658
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).

Supplementary data