An Assessment of Mass Flows, Removal and Environmental Emissions of Bisphenols in a Sequencing Batch Reactor Wastewater Treatment Plant

This study analyzed 16 bisphenols (BPs) in wastewater and sludge samples collected from different stages at a municipal wastewater treatment plant based on sequencing batch reactor technology. It also describes developing an analytical method for determining BPs in the solid phase of activated sludge based on solid-phase extraction and gas chromatography-mass spectrometry. Obtained concentrations are converted into mass flows, and the biodegradation of BPs and adsorption to primary and secondary sludge are determined. Ten of the sixteen BPs were present in the influent with concentrations up to 434 ng L−1 (BPS). Only five BPs with concentrations up to 79 ng L−1 (BPA) were determined in the plant effluent, accounting for 8 % of the total BPs determined in the influent. Eleven per cent of the total BPs were adsorbed on primary and secondary sludge. Overall, BPs biodegradation efficiency was 81%. The highest daily emissions via effluent release (1.48 g day−1) and sludge disposal (4.63 g day−1) were for BPA, while total emissions reached 2 g day−1 via effluent and 6 g day−1 via sludge disposal. The data show that the concentrations of BPs in sludge are not negligible, and their environmental emissions should be monitored and further studied.


Introduction
Bisphenols (BPs) are a group of synthetic organic compounds used in the production of epoxy resins and polycarbonate. These polymers are used in various applications, ranging from food contact materials, compact discs, construction materials, thermal paper, dental composites, medical equipment, water pipes, toys and sports equipment [1]. The most widely-known and produced BP is Bisphenol A (4,4 -(propane-2,2-diyl)diphenol, BPA), a known endocrine-disrupting compound. Concerns over its safety have resulted in restrictions on its use in many countries and its gradual replacement by other BPs in specific applications [2]. However, since all BPs share a common structure, i.e., two hydroxyphenyl groups (Table S1), they also have the potential to be endocrine-disrupting, as is the case for BPS, BPF and BPAF [3].
Bisphenols can enter the environment via different routes, although wastewater treatment plant (WWTP) effluents are considered one of the primary sources [3]. However, their fate during sewage treatment depends on the treatment technology and their physicochemical properties [4]. The end products of wastewater treatment are wastewater effluent and sewage sludge. Treated effluent is commonly discharged to surface water (river or sea) but is sometimes used to irrigate crops in areas that suffer a water deficit [5]. In addition, 2 of 15 depending on the WWTP, primary and secondary sludge can be further treated; for example, by aerobic or anaerobic digestion, which are used to stabilize the sludge. In Europe, sewage sludge is either spread on agricultural soils (50%), incinerated (28%), landfilled (18%), or disposed of through other methods (4%) [6].
Although BPs have been investigated in wastewater, a literature review reveals that only a limited number of studies address the presence of BPs in both wastewater (WW) and sludge (Table 1). For instance, only Huang et al. (2021) [7], Sun et al. (2017) [8] and Xue and Kannan (2019) [9] have investigated the fate of BPs during different stages of WW treatment. Huang et al. (2020) [10] report the presence of BPA, BPAF, BPE, BPF, BPS, BPB, BPZ, BPAP, BPP, BPBP, BPC, BPG, BPPH and BPTMC of up to 4537 ng L −1 (total) in influent, up to 569 ng L −1 in effluent and up to 878 ng g −1 in sludge. The most common BPs in influent, effluent and sludge at the highest concentrations were BPA, BPS and BPF. Except for Xue and Kannan (2019) [9], where the removal of BPA, BPF and BPS in two WWTPs was ≤29%, most authors report >87% removal of BPs from WW [10]. According to Hu et al. (2019) [11], the main removal mechanisms are biodegradation and adsorption.
The removal of BPs from sludge during anaerobic digestion remains poorly researched, but limited data suggest low (35%) or negligible removal of BPA, BPS and BPAF [12][13][14]. Interestingly, Abril et al. (2020) [15] found a five-fold increase in BPA concentrations during anaerobic digestion. Overall, the data reveal that BPs are generally more biodegradable under aerobic than anaerobic conditions [11], except for BPAF, which contains strong C-F bonds [4]. Despite these studies, gaps in the knowledge remain, such as the behavior of BPs at different stages of WW treatment and during different sludge treatments. Moreover, among all BPs, only the fate of BPA has been investigated in detail at different points of the technological processes of WW and sludge treatment [12]. In addition, no group has investigated the fate of BPs in WWTP utilizing sequencing batch reactors (SBR), nor the fate of BPs (except BPA) in primary and secondary sludge separately or in anaerobically stabilized sludge.
To address this gap, we studied the fate of BPs in a WWTP utilizing SBR technology and anaerobic sludge digestion. This work involved (1) developing a method for determining 16 BPs in sludge by gas chromatography-mass spectrometry (GC-MS), (2) analyzing their concentrations at different stages of WWTP, (3) determining the adsorption of BPs onto the primary and secondary sludge, (4) calculating their removal from WW and during the anaerobic sludge digestion and (5) evaluating the emissions of BPs into the environment via effluent release and sludge disposal.

Sample Collection
All samples were collected at the Domžale-Kamnik WWTP, Slovenia (149,000 population equivalents), which accepts municipal and industrial WW. The WWTP includes primary treatment (rakes, grease and sand traps, primary settler), secondary treatment (four SBRs) and sludge treatment using anaerobic sludge digestion. The plant accepts external WW and sludge from cesspits and small WWTPs, which enter directly to the primary settler or anaerobic digesters, respectively, and treats hazardous and non-hazardous liquid waste (which are first pre-treated by electrocoagulation) and then enter the mechanical stage. Before disposal, the anaerobically stabilized sludge is centrifuged. The centrate is then treated in a deammonification plant and returned to the mechanical stage. The yearly measurements of the volumetric flows show that external WW from cesspits and small WWTPs, effluents from electrocoagulation and the centrate from mechanical dewatering of anaerobically stabilized sludge represent only 0.25% of the total WWTP influent. For this reason, they were not sampled.
The sampling scheme for BPs monitoring is shown in Figure 1. Twenty-four-hour composite samples (V = 2 L) of WWTP influent (WWTP inf ), primary settler influent (PSE inf ), primary settler effluent (PSE eff ) and WWTP effluent (WWTP eff ), considering the WWTP hydraulic retention time (24-h), were sampled using automatic samplers in March 2021. Moreover, on the same day, 6-h composite samples (V = 2 L) of primary (PS) and secondary sludge (SS) and a grab sample (V = 2 L) of anaerobically stabilized sludge (AS) were collected manually. Since the sludge retention time of AS in the anaerobic digester is 30 times longer than that of PS and SS, only a grab sample was taken for the preliminary studies. In total, seven samples were collected. The WWTP operational parameters recorded during the sampling period are presented in Table S7. All of the aqueous samples were stored at −20 • C prior to analysis. Sludge samples were transferred to 50 mL centrifuge tubes and centrifuged at 10,000 RCF for 15 min. The aqueous phase was separated and stored at −20 • C. The remaining sludge was then lyophilized for 72 h at 75 mbar and −55 • C (CHRIST Gamma 1-16 LSCplus). All samples were stored at −20 • C.  [17]. All samples were defrosted overnight and then filtered sequentially through glass fiber (MN GF-2, Machery-Nagel, Dueren, Germany) and cellulose nitrate filters (0.45 µm, Sartorius, Göttingen, Germany) to prevent clogging of the solid-phase extraction (SPE) cartridges. In the case of the aqueous sludge phase, the samples were filtered using MN GF-4 filters. The filtered samples (250 mL) were then spiked with 25 µL of the internal standard mixture 13 C 12 -BPS, 13 C 12 -BPF, 13 C 12 -BPB, BPA-d 16 (c = 1 µg mL −1 ), and 500 µL of concentrated HCl were added. The samples were loaded onto Oasis Prime HLB cartridges (60 mg, 3 mL; Waters, Massachusetts, USA) at a flow rate of 3 mL min −1 . After loading, the sorbents were washed with 3 mL of 10% MeOH in water and dried under vacuum (−1.33 kPa) for 45 min. The elution step was performed using 5% FA in EtAc (3 × 0.6 mL). The solvent was evaporated under nitrogen at 40 • C.

Extraction of BPs from the Aqueous Phase
Aqueous samples were prepared following the method of Kovačič et al. (2019) [17 All samples were defrosted overnight and then filtered sequentially through glass fib (MN GF-2, Machery-Nagel, Dueren, Germany) and cellulose nitrate filters (0.45 μm, Sa torius, Göttingen, Germany) to prevent clogging of the solid-phase extraction (SPE) ca tridges. In the case of the aqueous sludge phase, the samples were filtered using MN G 4 filters. The filtered samples (250 mL) were then spiked with 25 μL of the internal stan ard mixture 13 C12-BPS, 13 C12-BPF, 13 C12-BPB, BPA-d16 (c = 1 μg mL −1 ), and 500 μL of conce trated HCl were added. The samples were loaded onto Oasis Prime HLB cartridges (6 mg, 3 mL; Waters, Massachusetts, USA) at a flow rate of 3 mL min −1 . After loading, th sorbents were washed with 3 mL of 10% MeOH in water and dried under vacuum (−1.3 kPa) for 45 min. The elution step was performed using 5% FA in EtAc (3 × 0.6 mL). Th solvent was evaporated under nitrogen at 40 °C.

Extraction of BPs from the Solid Phase
Optimization of SPE (Solid Phase) The method was optimized for extraction solvent, centrifugation and number of e tractions. Prior to clean-up with SPE, two additional clean-up steps using the QuEChER method (Bond Elut, Fruits and Veg, EN, Agilent Technologies, Santa Clara, CA, USA) an filtration (CHROMAFIL ® Xtra PTFE-45/25, Macherey-Nagel, Düren, Germany) were i vestigated. Full details are in the Supplementary Material S2.
Two types of SPE cartridges were assessed for acidified (conc. HCl, pH 2) and no acidified samples. These included (1) Oasis Prime HLB cartridges (60 mg, 3 mL; Water Massachusetts, USA), which are based on water-wettable, hydrophilic-lipophilic-ba anced divinylbenzene-N-vinylpyrrolidone copolymer and (2) Affinimip ® SPE cartridg designed for extracting bisphenols (100 mg, 6 mL; Affinisep, Petit Couronne, France which are based on molecularly imprinted polymers. In the case of Oasis Prime HLB ca tridges, the tested parameters (washing step and elution solvents) are shown as a sch matic in Figure S7.
For the Affinimip ® SPE Bisphenols cartridges, two procedures were tested accordin to the manufacturer's instructions. In the first procedure ( Figure S9a), the solvent w evaporated under N2 at 40 °C, reconstituted in 1 mL of MeOH/MeCN (1:1, v/v), mixe and 9 mL of water were added. Cartridges were placed on a vacuum manifold (Agile Technologies) and conditioned with 10 mL of 2% FA in MeOH (v/v), 4 mL of MeCN an

Extraction of BPs from the Solid Phase Optimization of SPE (Solid Phase)
The method was optimized for extraction solvent, centrifugation and number of extractions. Prior to clean-up with SPE, two additional clean-up steps using the QuEChERS method (Bond Elut, Fruits and Veg, EN, Agilent Technologies, Santa Clara, CA, USA) and filtration (CHROMAFIL ® Xtra PTFE-45/25, Macherey-Nagel, Düren, Germany) were investigated. Full details are in the Supplementary Material S2.
Two types of SPE cartridges were assessed for acidified (conc. HCl, pH 2) and nonacidified samples. These included (1) Oasis Prime HLB cartridges (60 mg, 3 mL; Waters, Massachusetts, USA), which are based on water-wettable, hydrophilic-lipophilic-balanced divinylbenzene-N-vinylpyrrolidone copolymer and (2) Affinimip ® SPE cartridges designed for extracting bisphenols (100 mg, 6 mL; Affinisep, Petit Couronne, France), which are based on molecularly imprinted polymers. In the case of Oasis Prime HLB cartridges, the tested parameters (washing step and elution solvents) are shown as a schematic in Figure S7.
For the Affinimip ® SPE Bisphenols cartridges, two procedures were tested according to the manufacturer's instructions. In the first procedure ( Figure S9a Solid phase sludge samples were prepared as follows: a known amount of lyophilized sludge (0.2 g) was transferred to a 10 mL centrifuge tube and spiked with 100 µL of the internal standard mixture 13 C 12 -BPF, 13 C 12 -BPS, 13 C 12 -BPB (c = 50 ng mL −1 ) and 160 µL of BPA-d 16 (c = 100 ng mL −1 ). An aliquot of MeOH (1 mL) was added, and the contents were homogenized (vortex) for 1 min. The samples were then shaken (oscillating shaker) at 100 rpm overnight, after which 5 mL of MeOH/MeCN (1:1, v/v) were added. The samples were then ultrasonicated (VWR MODEL 250D, 180 W) for 15 min, centrifuged at 12,000 RCF for 15 min, and the supernatant was removed. The extraction process was repeated twice, the extracts were combined, diluted with H 2 O (90 mL) and acidified (pH 2) using concentrated HCl (50 µL).
For the extraction step, Oasis Prime HLB cartridges were chosen since they outperformed AFFINIMIP ® SPE in terms of recovery (by 4%), repeatability (by 7%), time (2-fold) and cost (3-fold). The samples were loaded onto the cartridges at a flow rate of 3 mL min −1 , washed with 3 mL of 10% MeOH in H 2 O and dried under vacuum (−1.33 kPa) for 30 min. The elution step was performed using 5% FA in EtAc (3 × 0.6 mL). The extract was then dried under N 2 at 40 • C.

Instrumental Analysis
The method followed the protocol of Kovačič et al. (2019) [17]. Briefly, before analysis, all samples were derivatized by adding 50 µL of MSTFA and 50 µL of pyridine to the dry residue, followed by heating at 80 • C for 1 h. The samples were then analyzed using a 7890B series gas chromatograph (GC) coupled to a 5977A single quadrupole mass selective detector (Agilent, USA). Chromatographic separation was achieved using a DB-5 MS capillary column (30 m × 0.25 mm × 0.25 µm; Agilent, USA) with helium as the carrier gas (1 mL min −1 ). Samples (1 µL) were injected in splitless mode at 250 • C. The GC oven temperature program was as follows: an initial temperature of 120 • C was ramped at 20 • C min −1 to 250 • C and held for 6 min, then ramped to 300 • C at 10 • C min −1 and held for 3 min. The total runtime was 22 min. The mass spectrometer was operated in electron impact (EI) mode at 70 eV. The target compounds were identified and quantified using selected ion monitoring (SIM mode). The monitored SIM ions and retention times (RTs) for the derivatized BPs and internal standards are presented in Table S2. The data were processed using MassHunter software (Agilent Technologies).

Method Validation
Method performance for sludge analysis was assessed regarding recovery, linearity, accuracy, the limit of detection (LOD) and quantification (LOQ), sensitivity, precision (method and instrumental repeatability) and matrix effect. Recovery at low concentrations ranged from 62% (BPBP) to 107% (BPA) and from 81% (BPA) to 131% (BPAP) at high concentrations. The LOD ranged from 0.01 ng g −1 (BPPH) to 8.43 ng g −1 (BPA), and the LOQ from 0.03 ng g −1 (BPPH) to 28.10 ng g −1 (BPA). Full details are provided in the Supplementary Material (S4). Method validation of the aqueous phases was performed previously by Kovačič et al. (2019) [17]. In this case, recovery at low concentrations ranged from 41% (BPAF) to 115% (BPAP) and from 70% (BPPH) to 106% (BPAF) at high concentrations. The LOD ranged from 0.10 ng g −1 (BPC) to 5.22 ng g −1 (BPA), and LOQ ranged from 0.34 ng g −1 (BPC) to 17.39 ng g −1 (BPA). Again, full details are provided in Supplementary Material (S5). All glassware was cleaned according to standard laboratory practice for trace analysis, and procedural blanks were prepared for each experimental setup. Contamination was assessed by evaluating the ratio of the areas between the quantifier ions of the targeted BPs and the internal standards in the procedural blanks, which revealed the presence of BPA and BPS. All sample results were blank corrected, and solvent blanks (EtAc) were analyzed after every tenth sample to evaluate potential carryover.

Calculations of Mass Flows and BPs Removal
The mass flows of BPs in WWTP inf , PSE inf , PSE eff and WWTP eff were calculated using Equation (1): . .
.  Table S6), the most abundant in the WWTP influents/effluents were BPA and BPS and, to a lesser extent, 22BPF, 24BPF, 44BPF and BPE. The highest concentrations were determined for BPS in WWTP inf ≤434 ng L −1 , PSE inf ≤578 ng L −1 , PSE eff ≤591 ng L −1 and BPA in WWTP eff ≤79 ng L −1 . Below the LOQ (Table S6)  Concentrations of BPs in WW influent and effluent were much higher than repo by Česen et al. (2018) [2], who analyzed WW from the same WWTP. The obtained concentrations in the influent and effluent are seven times higher compared with li ture values (Table S9), while the obtained BPA concentrations are usually lower tha ported. In our study, the levels of BPS in WWTPinf, PSEinf and PSEeff are higher than of BPA. However, this is usually not the case, although Caban and Stepnowski (2020) observed higher levels of BPS than BPA [21]. There are several possible origins of BP WW, but a likely source is from recycled paper used to make toilet paper [22]. How in the case of the Domžale-Kamnik WWTP, several possible industrial sources have identified, including wastewater from a paint and lacquer manufacturer, a pharmac cal factory and a textile cleaning company [2].

BPs in the Aqueous Phase of Sludge
The most abundant BPs in the aqueous sludge phase (Figure 3, Table S6) were BPS, 44BPF and 24BPF. The highest concentrations of total BPs were in the anaerobi stabilized sludge, followed by the primary sludge, and the lowest total amount of BPs in the secondary sludge in both aqueous and solid phases. One reason is that only a s grab sample of anaerobically stabilized sludge was collected, and since its retention is 30 days, the concentrations cannot be directly related to those in the primary and ondary sludge. Adsorption and desorption of the compounds in the anaerobic diges deconjugation of conjugated compounds and possibly the addition of external sludge goes directly to anaerobic stabilization could also be factors. Moreover, anaerobic d tion of sludge can increase the concentration of trace organic compounds in the anae cally stabilized sludge due to solids reduction [13]. The total concentration of BPs in the effluent from the first part of the mechanical stage (rakes, grease and sand trap) is higher than the BPs levels in WWTP influent (Figure 2). This increase is possible since three potential sources of BPs are introduced to influent prior to the primary settler: (1) external WW (cesspits) delivered by trucks, (2) effluent from a deammonification plant treating centrate from mechanical dewatering of the anaerobically stabilized sludge and (3) pre-treated wastewater deriving from the treatment of hazardous and non-hazardous liquid wastes from the electrocoagulation plant. Moreover, the total amount of BPs is also higher in the effluent from the primary settler. Given that the inflow and outflow should be approximately equivalent and the sample was a 24-h composite sample, a likely explanation for the observed increase is a result of desorption processes from the solid particles, given that anaerobic conditions could occur in the settled sludge at the bottom of the primary settler [19]. Gu et al. (2021) [20] observed a similar trend of increasing concentration through the mechanical stage for BPA.
Concentrations of BPs in WW influent and effluent were much higher than reported byČesen et al. (2018) [2], who analyzed WW from the same WWTP. The obtained BPS concentrations in the influent and effluent are seven times higher compared with literature values (Table S9), while the obtained BPA concentrations are usually lower than reported. In our study, the levels of BPS in WWTP inf , PSE inf and PSE eff are higher than that of BPA. However, this is usually not the case, although Caban and Stepnowski (2020) also observed higher levels of BPS than BPA [21]. There are several possible origins of BPS in WW, but a likely source is from recycled paper used to make toilet paper [22]. However, in the case of the Domžale-Kamnik WWTP, several possible industrial sources have been identified, including wastewater from a paint and lacquer manufacturer, a pharmaceutical factory and a textile cleaning company [2].

BPs in the Aqueous Phase of Sludge
The most abundant BPs in the aqueous sludge phase (Figure 3, Table S6) were BPA, BPS, 44BPF and 24BPF. The highest concentrations of total BPs were in the anaerobically stabilized sludge, followed by the primary sludge, and the lowest total amount of BPs was in the secondary sludge in both aqueous and solid phases. One reason is that only a single grab sample of anaerobically stabilized sludge was collected, and since its retention time is 30 days, the concentrations cannot be directly related to those in the primary and secondary sludge. Adsorption and desorption of the compounds in the anaerobic digesters, deconjugation of conjugated compounds and possibly the addition of external sludge that goes directly to anaerobic stabilization could also be factors. Moreover, anaerobic digestion of sludge can increase the concentration of trace organic compounds in the anaerobically stabilized sludge due to solids reduction [13].

BPs in the Solid Phase of Sludge
The most abundant BP in the solid sludge phase was BPA and BPS, 22BPF, 24B and 44BPF, to a lesser extent (Figure 4, Table S6). In other research studies, the leve BPs are mainly reported in the solid phase of sludge. In the present study, the total c centrations of BPs in the solid phase of primary and secondary sludge are compar with previous studies (Table 1). Also, the total concentrations of BPs in primary and ondary sludge are lower than those measured in the domestic WWTPs in Korea by Le al. (2015) [23]. The authors also suggest that the paper and textile industries are source BPA and BPS, while concentrations of BPF are associated with domestic activity. Gu et al. (2015) [24] reported BPA levels of 1300 ng g −1 and 520 ng g −1 in stabilized sludge f thermophilic and mesophilic anaerobic digestors, respectively. These values are hig than the 365 ng g −1 reported by Abril et al. (2020) [15] in anaerobically digested and de drated sludge.

BPs in the Solid Phase of Sludge
The most abundant BP in the solid sludge phase was BPA and BPS, 22BPF, 24BPF, and 44BPF, to a lesser extent (Figure 4, Table S6). In other research studies, the levels of BPs are mainly reported in the solid phase of sludge. In the present study, the total concentrations of BPs in the solid phase of primary and secondary sludge are comparable with previous studies (Table 1). Also, the total concentrations of BPs in primary and secondary sludge are lower than those measured in the domestic WWTPs in Korea by Lee et al. (2015) [23]. The authors also suggest that the paper and textile industries are sources of BPA and BPS, while concentrations of BPF are associated with domestic activity. Guerra et al. (2015) [24] reported BPA levels of 1300 ng g −1 and 520 ng g −1 in stabilized sludge from thermophilic and mesophilic anaerobic digestors, respectively. These values are higher than the 365 ng g −1 reported by Abril et al. (2020) [15] in anaerobically digested and dehydrated sludge.
BPs are mainly reported in the solid phase of sludge. In the present study, the total concentrations of BPs in the solid phase of primary and secondary sludge are comparable with previous studies (Table 1). Also, the total concentrations of BPs in primary and secondary sludge are lower than those measured in the domestic WWTPs in Korea by Lee et al. (2015) [23]. The authors also suggest that the paper and textile industries are sources of BPA and BPS, while concentrations of BPF are associated with domestic activity. Guerra et al. (2015) [24] reported BPA levels of 1300 ng g −1 and 520 ng g −1 in stabilized sludge from thermophilic and mesophilic anaerobic digestors, respectively. These values are higher than the 365 ng g −1 reported by Abril et al. (2020) [15] in anaerobically digested and dehydrated sludge.

Mass Flows
Mass flows of BPs in different WW and sludge were calculated from their concentrations (Chapter 3.1), volumetric flow, TSS and TSS% (Table S7). Mass flows were not calculated for the following compounds: (a) BPFL, since its concentrations were below LOQ in all cases; (b) BPB, BPC, BPZ and BPPH, which were detected only in the sludge; and (c) BPC2, since its concentration was <LOQ in the PSEinf and the PSEeff. For other compounds

Mass Flows
Mass flows of BPs in different WW and sludge were calculated from their concentrations (Chapter 3.1), volumetric flow, TSS and TSS % (Table S7). Mass flows were not calculated for the following compounds: (a) BPFL, since its concentrations were below LOQ in all cases; (b) BPB, BPC, BPZ and BPPH, which were detected only in the sludge; and (c) BPC2, since its concentration was <LOQ in the PSE inf and the PSE eff . For other compounds <LOQ, the concentrations were substituted with a value equal to LOQ/2 [18]. Therefore, mass flows were calculated for BPS, 22BPF, 24BPF, 44BPF, BPE, BPA, BPAF, BPAP, BPBP and BPP ( Figure 5). The data show that the highest bisphenols in the WWTP inf were below 8.15 g day −1 (BPS), ≤10.85 g day −1 (BPS) in the PSE inf , ≤11.09 g day −1 (BPS) in the PSE eff , ≤1.48 g day −1 (BPA) in the WWTP eff , ≤1.75 g day −1 (BPA) in primary sludge, ≤0.17 g day −1 (BPA) in the secondary sludge and ≤4.63 g day −1 (BPA) in the anaerobically stabilized sludge. Among all the BPs, the highest mass flows were observed for BPS (≤11.09 g day −1 ) in PSE eff and BPA (≤9.07 g day −1 ) in PSE eff , which confirms their high usage [4].
When plotted (Figure 5), the BPs mass flow data indicate a similar trend to BPs concentrations, i.e., increase during the mechanical treatment stage and reaching a maximum in the PSE eff . The reason is that the volumetric flows of WW were generalized as one constant volumetric flow. For sludge, the highest mass flows of BPs were present in the anaerobically stabilized sludge, lower in the primary sludge and lowest in the secondary sludge. Possible explanations are single grab sampling of anaerobically stabilized sludge, adsorption and desorption of the compounds in the anaerobic digesters, deconjugation of conjugated compounds, the addition of external sludge and the increase due to solids reduction.
The mass flows of BPA in primary (1.75 g day −1 ), secondary (0.160 g day −1 ) and anaerobically stabilized sludge (4.63 g day −1 ) were high compared to the other BPs. Due to high BPS mass flows in WW, higher BPS mass flows were expected in sludge than were detected. Possible explanations can be a difference in logK ow , since BPA has higher logK ow (3.64) and is, therefore, more likely to adsorb to sludge than BPS (logK ow = 1.65). The observation is in contrast to Karthikraj and Kannan (2017) [16], who suggest that BPS has a higher affinity for particulate matter/sludge than BPA, but confirms the results of Xue and Kannan (2019) [9], who observed that BPA has a higher affinity towards particulate matter/sludge than BPS and BPF.
detected. Possible explanations can be a difference in logKow, since BPA has higher lo (3.64) and is, therefore, more likely to adsorb to sludge than BPS (logKow = 1.65). The servation is in contrast to Karthikraj and Kannan (2017) [16], who suggest that BPS h higher affinity for particulate matter/sludge than BPA, but confirms the results of Xue Kannan (2019) [9], who observed that BPA has a higher affinity towards particulate ter/sludge than BPS and BPF. The composition of BPs in selected wastewater flows and sludge ( Figure 6) sh that the most abundant BPs in WWTPinf, PSEinf and PSEeff are BPS and BPA, whi WWTPeff, BPA is more abundant than BPS. In PS, SS and AS, BPA is more abundant BPS. In the case of SS and AS, 44BPF is more abundant than BPS, which agrees with et al., 2019 [25]. The total mass flow distribution between the aqueous and solid phas primary, secondary and anaerobically stabilized sludge (Table S8) reveals that most are in the solid phase (85-95%), while only 5-15% are in the aqueous phase of sludge The composition of BPs in selected wastewater flows and sludge ( Figure 6) shows that the most abundant BPs in WWTP inf , PSE inf and PSE eff are BPS and BPA, while in WWTP eff, BPA is more abundant than BPS. In PS, SS and AS, BPA is more abundant than BPS. In the case of SS and AS, 44BPF is more abundant than BPS, which agrees with Zhu et al., 2019 [25]. The total mass flow distribution between the aqueous and solid phases in primary, secondary and anaerobically stabilized sludge (Table S8)

Distribution and Removal of BPs from Wastewater
The removal of BPs from wastewater can result from biodegradation, adsorption, chemical degradation or photodegradation, although the latter is unlikely since the sequencing batch reactors are covered at the WWTP Domžale-Kamnik. Biodegradation of BPs (Figure 7) ranged from 45% (22BPF) to 96% (BPE); while between 1% (BPBP) and 45% (22BPF) were adsorbed to the primary sludge, less than 5% (BPAP) adsorbed to secondary sludge, and 2% (24BPF, 44BPF, BPE) to 31% (BPP) remained in the WWTPeff. In total, 10% of the BPs were adsorbed to primary sludge, 1% to secondary sludge, 8% remained in the WWTPeff, and 81% were biodegraded. Total removal, the sum of biodegradation and adsorption, was 92%, which is consistent with 87-94% removals obtained in China and India (Tables 1 and S9) and >96% removal (BPA not included) obtained in Slovenia in a previous study [2]. A comparison of the removals of all BPs with the literature data can be found in

Distribution and Removal of BPs from Wastewater
The removal of BPs from wastewater can result from biodegradation, adsorption, chemical degradation or photodegradation, although the latter is unlikely since the sequencing batch reactors are covered at the WWTP Domžale-Kamnik. Biodegradation of BPs (Figure 7) ranged from 45% (22BPF) to 96% (BPE); while between 1% (BPBP) and 45% (22BPF) were adsorbed to the primary sludge, less than 5% (BPAP) adsorbed to secondary sludge, and 2% (24BPF, 44BPF, BPE) to 31% (BPP) remained in the WWTP eff . In total, 10% of the BPs were adsorbed to primary sludge, 1% to secondary sludge, 8% remained in the WWTP eff , and 81% were biodegraded. Total removal, the sum of biodegradation and adsorption, was 92%, which is consistent with 87-94% removals obtained in China and India (Table 1 and Table S9) and >96% removal (BPA not included) obtained in Slovenia in a previous study [2]. A comparison of the removals of all BPs with the literature data can be found in Table S9.

Removal of BPs in Anaerobic Digesters
Since sludge retention time in the anaerobic digesters (30 days) is much longer t a single treatment cycle of WW (24 h), it is impossible to compare the results of AS m flow with the PS and SS mass flows directly. Moreover, only a single grab sample of anaerobically stabilized sludge was taken. Still, measuring the removal of BPs in ana bic digesters can provide useful information. For instance, the removal of BPAF among the highest of all BPs, which was expected since halogenated compounds known to be more biodegradable under anaerobic conditions [26]. The removal of B 22BPF, BPE, BPAF, BPAP and BPP (Table 2) ranged from 11% (BPE) to 82% (BPAP), w an increase in mass flows after anaerobic digestion was observed for 24BPF, 44BPF, and BPBP. The latter is likely due to a single grab sampling of anaerobically stabil sludge, adsorption and desorption of the compounds in the anaerobic digesters, deco gation, external sludge addition or, most likely, volatile solids reduction [13]. In o studies, Samaras et al. (2013) [12] reported low removal of BPA (35%), while Phan e (2018) [14] reported negligible removal of BPA; in addition to BPA, Choi et al. (2021) reported negligible removal of BPS and BPAF.
Interestingly, Abril et al. (2020) [15] reported a five-fold increase in the levels of during anaerobic digestion, i.e., from 45 ng g −1 in primary sludge and 100 ng g −1 in ondary sludge to 245 ng g −1 in anaerobically stabilized sludge. This finding agrees w 2.4-fold increase in mass flows of BPA observed in this study, although in our case le of BPA were much higher. The results also suggest that BPAP and 22BPF are remo more under anaerobic than under aerobic conditions. The reason could be due to the c pounds' structure or unknown anaerobic biodegradation mechanisms.  The lowest biodegradation was observed for 22BPF, BPAF, BPAP and BPP. BPAF is a halogenated bisphenol, and its low biodegradation is consistent with compounds with strong C-F bonds being poorly biologically degradable under aerobic conditions [8,17]. BPAP and BPP had a similar biodegradation rate to BPAF, and all had low mass flows. However, the way biodegradation is calculated, considering that LOQ/2 was used when the measured concentration was below LOQ (Chapter 2.6), can affect their calculated distribution, so further investigation is needed to confirm their biodegradation rates. A possible reason for the low biodegradation of 22BPF compared to its other two isomers, 24BPF and 44BPF, could be a consequence of the different positions of the two hydroxyl groups. According to Noszczyńska and Piotrowska-Seget (2018) [1], the degradation pathway of 44BPF described for bacteria Sphingobium yanoikuyae FM-2 starts with the formation of bis(4-hydroxyphenyl)methanol after hydroxylation of the 44BPF bridging carbon. Since the two hydroxyl groups of 22BPF are in the second position on the aromatic rings, they are likely to hinder the hydroxylation of the bridging carbon, which may explain why 22BPF is less biodegraded than its isomers.

Removal of BPs in Anaerobic Digesters
Since sludge retention time in the anaerobic digesters (30 days) is much longer than a single treatment cycle of WW (24 h), it is impossible to compare the results of AS mass flow with the PS and SS mass flows directly. Moreover, only a single grab sample of the anaerobically stabilized sludge was taken. Still, measuring the removal of BPs in anaerobic digesters can provide useful information. For instance, the removal of BPAF was among the highest of all BPs, which was expected since halogenated compounds are known to be more biodegradable under anaerobic conditions [26]. The removal of BPS, 22BPF, BPE, BPAF, BPAP and BPP (Table 2) ranged from 11% (BPE) to 82% (BPAP), while an increase in mass flows after anaerobic digestion was observed for 24BPF, 44BPF, BPA and BPBP. The latter is likely due to a single grab sampling of anaerobically stabilized sludge, adsorption and desorption of the compounds in the anaerobic digesters, deconjugation, external sludge addition or, most likely, volatile solids reduction [13]. In other studies, Samaras et al. (2013) [12] reported low removal of BPA (35%), while Phan et al. (2018) [14] reported negligible removal of BPA; in addition to BPA, Choi et al. (2021) [13] reported negligible removal of BPS and BPAF.  [15] reported a five-fold increase in the levels of BPA during anaerobic digestion, i.e., from 45 ng g −1 in primary sludge and 100 ng g −1 in secondary sludge to 245 ng g −1 in anaerobically stabilized sludge. This finding agrees with a 2.4-fold increase in mass flows of BPA observed in this study, although in our case levels of BPA were much higher. The results also suggest that BPAP and 22BPF are removed more under anaerobic than under aerobic conditions. The reason could be due to the compounds' structure or unknown anaerobic biodegradation mechanisms.

Conclusions
This study investigated the occurrence and mass flows of 16 BPs in WWTP based on SBR technology. According to the research goals, this study can be summarized as follows. A method for determining 16 BPs in sludge based on SPE and GC-MS with good average repeatability, recovery and low LOQ values was successfully developed, revealing that BPA and BPS were the most abundant compounds measured in WWTP influent, effluent and sludge and represent 80-90% of all the amounts of BPs in streams. The data also show that the total concentration of BPs increases during the mechanical stage of treatment, with the highest values in the primary settler effluent. In sludge, the highest concentrations of BPs were found in the anaerobically stabilized sludge, followed by primary and secondary sludge. The study also found that significant amounts of BPs remain in the WWTP effluent (8%) and the primary (10%) and secondary sludge (1%), with the majority being biodegraded (81%) in the sequencing batch reactors. Overall, removal was 92%, with the highest daily emissions from the WWTP being 1.48 g day −1 and 4.63 g day −1 for BPA via effluent and anaerobically stabilized sludge, and the sum of the mass flows of all BPs in the WWTP eff was 2.13 g day −1 , and in the anaerobically stabilized sludge it was 6.03 g day −1 . Finally, given the potential toxicity of BPs, these results could prove useful when assessing risk regarding the emissions into the environment and the reuse of wastewater and sludge in agriculture.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/molecules27238634/s1. Table S1: Names, chemical structures, and physicochemical properties of the studied compounds; Figure S1: Recoveries obtained using five different extraction solvents for 16 BPs; Figure S2: Recoveries obtained using different centrifugation parameters at each repetition of extraction for 16 BPs; Figure S3: Recoveries for 16 BPs obtained using solutions containing different MeOH contents in water for sorbent washing; Figure S4: Recoveries for 16 BPs obtained when the sample was nonacidified and acidified before loading onto the cartridges; Figure S5: Recoveries for 16 BPs obtained using five different elution solvents; Figure S6: Recoveries obtained without an additional clean-up step, cleaned with filtration and QuEChERS; Figure S7: Schematic of test parameters using Oasis Prime HLB cartridges; Figure S8: Comparison of the recoveries obtained using Oasis HLB Prime and Affinimip ® SPE Bisphenols cartridges following two different procedures and their variations; Figure S9: Schematic of Affinimip ® SPE Bisphenol extraction procedures; Table S2: A list of bisphenols and internal standards, monitored ions [m/z] (quantifier ion in bold and two qualifier ions) of derivatised compounds, and retention times (RT); Table S3: Validation parameters (method repeatability, instrumental repeatability, accuracy, k and R2 values, SPE recovery, LOD and LOQ) of an analytical method for determination of 16 BPs in the solid phase of sludge; Table S4: Matrix effect; Table S5: Validation parameters (k and R2 values, recovery, LOD, LOQ) of an analytical method for the determination of 16 BPs in the aqueous phase of sludge and WW; Table S6: Concentrations of BPs in WW and sludge at different points and in different phases; Table S7: Sampling time, volumetric flows, total suspended solids, and other basic parameters in wastewater and sludge flows; Table S8: Total BPs mass flows and their distribution in different types of sludge in aqueous and solid phase; Table S9: Comparison of mean concentrations of BPs in influent, effluent, sludge and their removal with literature [7,[9][10][11]16,17,27].