Characterizing the Organohalogen Iceberg: Extractable, Multihalogen Mass Balance Determination in Municipal Wastewater Treatment Plant Sludge

The large number and diversity of organohalogen compounds (OHCs) occurring in the environment poses a grand challenge to analytical chemists. Since no single targeted method can identify and quantify all OHCs, the size of the OHC “iceberg” may be underestimated. We sought to address this problem in municipal wastewater treatment plant (WWTP) sludge by quantifying the unidentified fraction of the OHC iceberg using targeted analyses of major OHCs together with measurements of total and extractable (organo)halogen (TX and EOX, respectively; where X = F, Cl, or Br). In addition to extensive method validation via spike/recovery and combustion efficiency experiments, TX and/or EOX were determined in reference materials (BCR-461 and NIST SRMs 2585 and 2781) for the first time. Application of the method to WWTP sludge revealed that chlorinated paraffins (CPs) accounted for most (∼92%) of the EOCl, while brominated flame retardants and per- and polyfluoroalkyl substances (PFAS) accounted for only 54% of the EOBr and 2% of the EOF, respectively. Moreover, unidentified EOF in nonpolar CP extracts points to the existence of organofluorine(s) with physical–chemical properties unlike those of target PFAS. This study represents the first multihalogen mass balance in WWTP sludge and offers a novel approach to prioritization of sample extracts for follow-up investigation.


■ INTRODUCTION
Of the millions of chemicals used in society, organohalogen compounds (OHCs) are a class of particular concern. Approximately 30,000 individual OHCs are registered for production and use on the global market. 1 The strength of the carbon−halogen bond imparts considerable stability and contributes to the environmental persistence, bioaccumulation, and toxicity of OHCs. 2 In silico screening of chemical inventories from Europe and North America revealed 3421 chemicals with both persistent and bioaccumulative properties, of which ∼52% were halogenated. 3 Well-known classes of highly halogenated pollutants include per-and polyfluoroalkyl substances (PFAS), chlorinated paraffins (CPs), and halogenated flame retardants (HFRs) such as polybrominated diphenyl ethers (PBDEs) and organophosphate esters (OPEs). However, only a small fraction of OHCs are regulated by the global UN Stockholm Convention on Persistent Organic Pollutants; many others are considered trade secrets and are either undisclosed, ambiguously described, or infrequently monitored. 1 Persistent OHCs used in society and industry eventually end up in wastewater and urban runoff, entering wastewater treatment plants (WWTPs). WWTP samples are therefore useful for monitoring stable chemicals used in society. Additionally, since sewage sludge may be released back into the environment via agriculture, recultivation/land reclamation, and landfilling (combined totaling ∼70% of all sludge in Europe 4 ), it is of great importance to ensure that harmful substances are identified and prevented from entering the wastewater system so that sludge can be reused in a sustainable manner. Capturing the large number and diversity of OHCs in a single targeted method is currently not possible, leading to concerns that OHC contamination may be underestimated. For this reason, there is growing interest in so-called "organohalogen mass balance" experiments, which seek to quantify the fraction of unidentified organohalogen in samples. This is assessed through paired analysis of target OHCs and extractable organohalogen (EOX; where X = F, Cl, or Br).
The earliest organohalogen mass balance experiments date back to the mid 1990s and early 2000s. Instrumental neutron activation analysis (INAA) was used to determine EOCl, EOBr, and EOI across a wide range of environmental samples, including biota, 5−10 incinerator ashes, 11 air, 12 and pine needles. 13 These studies mostly targeted legacy pollutants, including polychlorinated biphenyls (PCBs), dichlorodiphenyl trichloroethanes (DDTs), hexachlorocyclohexanes (HCHs), polychlorinated/brominated dibenzo-p-dioxins (PCDDs/ PBDDs), polychlorinated/brominated dibenzofurans (PCDFs/PBDFs), and PBDEs. More recently, PCBs, pesticides, and PBDEs were determined together with EOCl and EOBr in high trophic level mammals. 14 In 2007, the first fluorine mass balance experiments were performed, pairing target PFAS analysis with combustion ion chromatography (CIC)-based EOF determination in water and human blood. 15,16 Later, other matrices were analyzed, including biological samples, 17−22 consumer products, 23,24 and sewage sludge. 25−27 Collectively, these studies point to significant gaps in the EOX mass balance across a diverse range of sample types. However, despite nearly three decades of research activity, there remains fairly little data on the contribution of contemporary OHCs (e.g., CPs and emerging HFRs) to the total OHC burden, and even fewer data on TX and EOX in reference materials, which are necessary for widespread adoption and standardization of halogen mass balance methods.
In the current study, we aimed to combine measurements of EOF, EOCl, and EOBr, with analysis of emerging OHCs to identify and compare the known and unknown fractions of the organohalogen "iceberg" in WWTP sludge. This approach enables prioritization of extracts/extraction procedures which generate the largest fractions of unidentified halogen for further investigation. To achieve this objective, we developed and validated a multihalogen CIC method for determination of EOX and total X (TX; where X= F, Cl, or Br). This method, together with quantitative target analyses of emerging fluorinated, chlorinated, and brominated organic contaminants (e.g., PFAS, CPs, and HFRs), was applied to form the first multihalogen mass balance data set in both municipal WWTP sludge and standard reference materials.

■ MATERIALS AND METHODS
Sample Collection and Preparation. Digested (both aerobic and anaerobic) and dewatered sewage sludge was collected on November 14, 2019, at Henriksdal WWTP, situated in Stockholm, Sweden. The major treatment steps are illustrated in Figure S1 and include mechanical (coarse screens, grit chamber, primary sedimentation, sand filter), chemical (precipitation with iron), and biological (aerobic bioreactor, anaerobic digestion) treatments. The collected sludge is the final product after treatment. The Henriksdal WWTP has two inlets consisting mainly of domestic and commercial wastewater plus urban runoff, serving ∼870,000 people. The plant has an average sewage treatment capacity of 273,000 m 3 /day, 28 and the average residence time of the sludge is 19 days. After sampling, the sludge samples were freeze-dried and homogenized. Water content, determined by noting the change in mass before and after freeze-drying, was 67%. Organic matter (64%) was determined by loss on ignition, in which the freeze-dried sludge sample was burned overnight at 105°C followed by 2 h at 550°C. Thereafter, organic carbon (39%) was estimated using a conversion factor of 0.6. 29 The samples were stored in the freezer at −18°C prior to analysis.
Overview and Justification of Analytical Approach. While we initially considered a single holistic procedure for combined analysis of organofluorine, organochlorine, and organobromine substances, such a method is not currently available and would be extremely challenging to develop given the wide range of OHC physical−chemical properties and the vulnerability of OHC measurements to coextractable matrix interferences (which require treatment using specific cleanup procedures). Consequently, we focused our EOX mass balance on extracts produced from separate methods for measurement of three classes of OHCs: EOF/PFAS, EOCl/CPs, and EOCl +Br/HFRs. This approach compares EOX and target OHCs produced from the same (unique) extract and thus any gap between "known" and "total" EOX can reasonably be attributed to unknown substances with similar physical− chemical properties to the OHCs targeted by each method. We note that EOX is, by definition, an underestimate of "total organohalogen", since some organohalogens are inevitably not extracted and/or are lost during cleanup. Targets were chosen due to their elevated concentrations in WWTP sludge, 27,30−34 and by extension, expected contribution to the EOX mass balance. Moreover, targeted analytical methods for PFAS, CPs, and HFRs are well established in our laboratories. Other legacy halogenated classes such as PCBs and DDTs were not included, since their concentrations in previously analyzed sewage sludge from the same WWTP, sampled several times over the same year (2019) were consistently >45 times lower than medium chain length CPs (∑ 7 PCBs at 21 ng/g dry weight [dw] compared to 972−3540 ng/g dw for MCCPs), 30 and were not expected to contribute significantly to the organochlorine mass balance in this sludge sample.
Characterization of sludge samples involved four separate analyses for each OHC class ( Figure 1A). First, TX concentrations, which represent the sum of all organic and inorganic halogen (IX), were determined by direct combustion (i.e., no sample preparation). Second, targeted OHC analysis was performed on replicate (n = 3) sludge samples which were spiked with internal standard (IS) and then extracted (performed on all but CPs; see details in section on CP Extraction Procedure) and cleaned up. Concentrations derived from these extracts (Extract 1 in Figure 1A) are corrected for procedural losses and are used to present target OHC profiles. In parallel, separate portions (n = 3) of the same sludge were extracted and cleaned up without IS, and the resulting extracts (Extract 2 in Figure 1A) were each split into two fractions: one fraction was spiked with IS postextraction and analyzed for target OHCs (nonrecovery corrected, representing the "known" fraction of EOX), while the second fraction was analyzed for EOX by CIC (no IS; Figure 1A). Target OHC and EOX concentrations in Extract 2 (both of which are nonrecovery corrected) are directly comparable and are therefore used to assess the halogen mass balance.
While the extraction and cleanup procedures employed here remove inorganic halogens (confirmed via spike/recovery experiments; see below), they may also selectively exclude some nonextractable organohalogens (NEOX; including OHCs that are nonextractable and/or removed during cleanup; Figure 1B). For this reason, EOX concentrations are expected to underestimate total organohalogen concentrations. Moreover, halogens that are not extractable by one method, may be observed in another. For example, an extremely nonpolar organofluorine that is not extractable using our PFAS extraction method (methanol) might be observable in CP extracts. For this reason, in addition to pairing EOX measurements with their associated targeted OHCs (i.e., PFAS/EOF, CPs/EOCl, and HFRs/EOCl+EOBr; see Table S1 for a full list of substances and abbreviations, and Figure S2 for schematic workflow), we also performed EOX measurements on several extracts where a particular halogen was unexpected (e.g., EOF measurements on CP extracts).
PFAS Extraction Procedure. The PFAS extraction procedure was based on previous studies and is described in detail in the SI. 26,35 Briefly, ∼0.5 g freeze-dried sludge was fortified with ISs and then extracted with methanol, followed by EnviCarb cleanup (Extract 1; Figure 1A). The same procedure was used for EOF mass balance analysis (Extract 2; Figure 1A), but ISs were not added prior to extraction. Extract 2 was split: one portion was fortified with ISs for target PFAS analysis while the remainder was stored for EOF determination by CIC. Instrumental analysis for target PFAS was performed by ultraperformance liquid chromatography−tandem mass spectrometry (UPLC-MS/MS) and included 11 perfluoroalkyl carboxylic acids (PFCAs), 4 perfluoroalkyl sulfonic acids (PFSAs), and perfluorooctane sulfonamide (FOSA). Instrumental details are presented in our previous study, 21 and retention times, precursor and selected ions, ISs, and limits of quantification (LOQs) are presented in Table S2. Accuracy and precision of targeted analysis was assessed via replicate (n = 3) spike/recovery experiments for nonrecovery corrected Extract 2, performed by adding 5 ng of individual PFAS into sludge and analyzing them together with real samples (Table  S3). Additionally, replicate (n = 3) portions of standard reference material (SRM) 2781 ("domestic sludge", NIST, Gaithersburg, MD), with reported values for several PFAS, were analyzed for method performance comparison (Table  S4). Additional QC and method validation procedures are described below.
CP Extraction Procedure. The CP extraction procedure was based on accelerated solvent extraction (ASE) using an ASE350 (Dionex, US) using a previously developed method (see SI for details). 36 A separate targeted analysis (i.e., Extract 1 in Figure 1A; with ISs added prior to extraction) was not included because prior data have shown IS recoveries close to 100% (i.e., non-IS corrected values can be expected to be the same as IS-corrected values). 37 Briefly, ∼0.5 g portions of freeze-dried sludge were loaded into ASE cells. Extraction was performed at 100°C and 1500 psi using n-hexane:acetone (1:1) as the extraction solvent. The resulting extract was evaporated to near dryness with nitrogen and underwent three cleanup steps for lipid and elemental sulfur removal. The final extract was concentrated to ∼1 mL, after which a portion (150 μL) was spiked with 10 ng 13 C-labeled C 10 Cl 6 for target CP analysis; the remainder of the extract was analyzed for EOCl content by CIC. Targeted analysis of very short-, short-, medium-, and long-chained CPs (vSCCPs, SCCPs, MCCPs, LCCPs) was carried out by UPLC−atmospheric pressure chemical ionization−high resolution mass spectrometry (APCI-HRMS) using a QExactive HF Orbitrap (Thermo Fisher Scientific) in full scan (m/z 250−2000) mode with a resolution of 120,000 full width at half-maximum (fwhm). A total of 342 CP homologues (expressed as C n Cl m , n = 6−33 m = 1−16) were measured. Instrumental details are described elsewhere. 38 Quantification was carried out based on a C n Cl mprofile reconstruction method by Bogdal et al. 39 Sixteen commercial mixtures were used for quantification. The chlorine weight was calculated based on the semiquantitative determination of the degree of chlorination obtained by the instrument. The deviation from the manufacturer's values ranged from 0.6% to 2.5% Cl. LOQs were calculated from the average blanks levels plus 10 times the standard deviations (Table S5). Accuracy and precision of CP determination were assessed via spike/recovery experiments using portion of diatomaceous earth (n = 3) fortified with 196, 1736, and 1068 ng of SCCPs (55.5% Cl), MCCPs (52.0% Cl), and LCCPs (49% Cl), respectively (Table S6). Replicates (n = 3) of NIST SRM 2585 "organic contaminants in house dust" (NIST, Gaithersburg, MD) were also analyzed and compared to previous studies (Table S7). Additional QC and method validation procedures are described below.
HFR Extraction Procedure. For target HFR analysis, 2 g of freeze-dried sludge was weighed into a tube, fortified with isotope-labeled ISs, and then extracted using a previously published method (details in SI). 40 Thereafter, fractionation of the extract (Extract 1, Figure 1A) was performed using a silica column based on a previously published method 41 with slight modifications (see SI for details). The analytes were eluted in three fractions that were subsequently cleaned up individually ( Figure S3), prior to targeted analysis and EOX. Fraction I contained PBDEs, decabromodiphenyl ethane (DBDPE), and some emerging brominated flame retardants (BFRs). Fraction II included 2-ethylhexyl 2,3,4,5-tetrabromobenzoate (EH-TBB), 1,2-bis(2,4,6-tribromophenoxy)ethane (BTBPE), and bis(2-ethylhexyl) tetrabromophthalate (BEH-TEBP). Fraction IIIa contained tetrabromobisphenol A (TBBPA), α-β-, and γ-1,2,5,6,9,10-hexabromocyclododecane (HBCDD), and Fraction IIIb contained chlorinated OPEs. Fractions I and II were analyzed on a gas chromatography mass spectrometer (GC-MS) equipped with an electron capture negative ionization (ECNI) source. Since highly brominated PBDEs are prone to thermal degradation, 42 the octaBDEs, nonaBDEs, and decaBDEs were analyzed using a short column (15 m), while the remaining BDEs were analyzed on a 30 m column. Fraction IIIa was run on a LC-MS/ESI. Fraction IIIb analysis was carried out on a GC-MS coupled to an electron ionization (GC-MS/EI), as described previously. 41 Bromine and chlorine mass balance analysis was performed using the aforementioned method, but extracts (Extract 2, Figure 1A) were only fortified Environmental Science & Technology pubs.acs.org/est Article with ISs immediately prior to instrumental analysis. Retention times and LOQs are available in the SI (Table S8). IS yields in blanks and sludge samples are presented in Table S9. For targets with IS recoveries <20% or >150% (i.e., TBBPA, HBCDDs, and TCEP), no quantification was performed. Accuracy and precision in the nonrecovery corrected Extract 2 were evaluated using replicate (n = 3) spike/recovery experiments performed at concentrations ∼5 times higher than the measured/expected concentration in sludge (Table  S10). Additional QC and method validation procedures are described below. TX and EOX Analyses. TX and EOX measurements were performed using a CIC (Thermo-Mitsubishi), controlled by NSX-2100 software. Neat samples (∼5 mg) or extracts (100− 300 μL) were placed on a ceramic boat (for extracts, the boats contained glass wool for better dispersion). The boats, including glass wool, were baked prior to use. The samples were combusted slowly in a combustion furnace (HF-201, Mitsubishi) at 1100°C under a flow of oxygen (300 mL/min) and argon mixed with water vapor (200 mL/min) for approximately 5 min. Combustion gases were absorbed in Milli-Q water during the entire length of the combustion process using a gas absorber unit (GA-210, Mitsubishi). A 200 μL aliquot of the absorption solution was subsequently injected onto an ion chromatograph (Dionex Integrion HPIC, Thermo Fisher Scientific) equipped with an anion exchange column (2 mm × 50 mm guard column (Dionex IonPac AS19-4 μm) and 2 mm × 250 mm analytical column (Dionex IonPac AS19-4 μm)) operated at 30°C. Chromatographic separation was achieved by running a gradient of aqueous hydroxide mobile phase ramping from 8 to 60 mM at a flow rate of 0.25 mL/min (Table S11), and halogens were detected using a conductivity detector. Retention times were approximately 7.5, 11.5, and 14.4 min for F − , Cl − , and Br − , respectively.
CIC data were processed using Chromeleon 7.2, and quantification was performed using a dilution of the inorganic standard in deionized water (Combined Seven Anion Standard I, #056933, Dionex). The 9-point calibration ranged from 0.05 to 20 ng/μL for F − , 0.075−30 ng/μL for Cl − , and 0.25−100 ng/μL for Br − and showed good linearity (R 2 > 0.99). The mean concentrations in the boat blanks were subtracted from the samples. The instrumental detection limit (DL), defined as the mean concentration plus three times the standard deviation in the boat blanks, was 0.11 ng/μL for F − and 0.34 ng/μL for Cl − . Since Br − was not detected in instrumental blanks, it was determined by the lowest detectable calibration point at 0.25 ng/μL.
The quantified target OHC levels were converted to halogen equivalents (i.e., mass of halogen per mass of sample [ng X/g]) according to eq 1, where C OHC is the concentration of OHC (ng/g), n X is the number of halogens in the molecule, A X is the atomic mass of the halogen (g/mol), and MW OHC is the molecular weight of the OHC (g/mol).
Quality Control, Method Validation, and Reference Material Characterization. Halogen mass balance experiments require careful quality control and method validation in order to ensure comparability between target OHCs, EOX, and TX data. While reference materials are available for validating target OHC measurements across a wide range of matrices, very few of these materials exist for TX or EOX determination. In the present study, we performed a series of method validation experiments in order to ensure high quality halogen mass balance data. Moreover, we provide here some of the first TX and EOX data on commercially available standard reference materials (SRMs), in order to facilitate interlaboratory comparison and standardization of halogen mass balance experiments. An overview of both QC and method validation experiments, as well as reference material characterization is provided here: • Accuracy/precision: Replicate spike/recovery experiments were carried out for most target OHCs (details in individual Materials and Methods section). In addition, NIST SRM 2585 (house dust), NIST SRM 2781 (domestic sludge), and BCR-461 (fluorine in clay) were analyzed for CPs, PFAS, and TF, respectively, and compared to reference values, including an interlaboratory comparison on EOF measurements. 26 Finally, a comparison of target data generated with and without fortification of ISs prior to extraction was included for most targets (except CPs).  (Table  S7). 43−45 For targeted analysis of chlorinated flame retardants, IS recoveries were acceptable (46%; RSD: 23%) for TDCPP, but low (<20%) for TCEP; consequently, only the former target was quantified. TCIPP blank levels (22.6 ng) were significant and were therefore subtracted from the sludge samples. Finally, for targeted analysis of BFRs, no contamination was detected in the procedural blanks. Spike/recovery results for BFRs generally ranged from 68% to 133% (RSD 2%−26%), indicating reasonable accuracy and precision in most cases (Table S10). Exceptions were BDE196, BDE206, TBP-AE, BTBPE, and BEH-TEBP, which showed lower recoveries (38%, 59%, 6%, 42%, 22%, respectively), and DBDPE, PBBz, PBT, TBP-DBPE, HBB, DBHCTD, and DDC-CO (syn), which showed higher recoveries (178%, 199%, 261%, 251%, 407%, 161%, respectively), likely due to matrix effects. Nevertheless, most of these over-and under-recovered targets  Samples where organohalogen mass balance was performed are highlighted, and percentages indicate the fraction of EOX that was characterized. Individual target analyte concentrations are presented in Table S16. a Extraction blank subtracted (−47 ng F/g). N/A = not analyzed. N/A* = not analyzed due to high instrumental blank levels.
are not likely to affect the EOBr mass balance, since they account for a small proportion of the EOBr.
Overall, these results point to reasonable method accuracy and precision across all target OHCs included in the present work, providing confidence in assessment of known EOX concentrations in sludge samples.
Quality Control of Halogen Mass Balance. Samples spiked with inorganic halogens displayed percent recoveries of <4% for NaF (PFAS extraction) and <0.1% for NaCl (CPs extraction), with measured EOX concentrations that were not significantly different from unspiked samples. Fluoride concentrations were also measured in the final sludge extract using a fluoride ISE, and were not statistically different from the procedural blank (0.017 ± 0.006 vs 0.010 ± 0.002 μg F/g, respectively), nor did they make a substantial contribution to EOF (<5.6%; Table S12). Nevertheless, caution is warranted when using this extraction method procedure, since highfluoride content samples could lead to elevated fluoride in the final extract (Table S12).
Combustion efficiencies among organic and inorganic halogen standards, as well as in the presence of matrix (BCR-461 clay) revealed consistent halogen-specific responses ( Figure 2; Table S13), with most combustion efficiencies ranging from 84 ± 16% (TClPP) to 110 ± 31% (BDE209), and good accuracy of 104 ± 4% found for BCR-461 (fluorine only). The exception was for BDE183 which displayed a slightly higher and more variable combustion efficiency (137 ± 40%) at low concentrations. These data are in reasonable agreement with a prior study investigating combustion efficiencies of 13 PFAS, which showed consistent combustion efficiency for most substances. 46 Overall, combustion efficiencies for inorganic and organic F, Cl, and Br were sufficiently accurate and precise for carrying out OHC mass balance experiments, although there is a need for EOX interlaboratory studies and reference materials to further standardize these methods for widespread use.
Characterization of Standard Reference Materials. Direct combustion of NIST domestic sludge SRM 2781 resulted in concentrations of 813, 2225, and 120 μg X/g dw for F − , Cl − and Br − , respectively, and in NIST house dust SRM 2585, TX concentrations were 315 ± 8, 8820 ± 94, and 471 ± 111 μg X/g dw, F − , Cl − , and Br − , respectively (Table 1). In addition, TCl and TBr measurements in BCR-461 clay (which is certified for total fluorine) produced concentrations of 69 ± 9 μg Cl/g dw and nondetectable levels of Br. To the best of our knowledge, this is the first time NIST SRMs have been characterized for total halogens. Since TX measurements do not require sample handling, these data enable interlaboratory comparison and standardization of instrumental analysis methods.
In addition to TX determination, organohalogen mass balance measurements were carried out on NIST reference materials for the first time. SRM 2781 (sludge) was extracted with the PFAS extraction method and contained 3590 ± 311 ng F/g dw of EOF, of which ∼6% (i.e., 226 ± 10 ng F/g) was accounted for by the 16 PFAS included in the present work (Table 1; Table S4). SRM 2585 (dust) extracted using the CP extraction procedure contained 417 ± 77 ng F/g dw, 17,000 ± 1600 ng Cl/g dw, and 2340 ± 433 ng Br/g dw. Essentially all (∼106%; 17,900 ± 300 ng Cl/g) of the EOCl was accounted for by CPs (Table 1; Table S7). We note that there are a number of other chlorinated chemicals with certified or reference values present in this SRM (e.g., PCBs, DDE and DDT, chlordane, heptachlor, and nonachlor), which when summed produce a chlorine equivalent of ∼1.8 μg Cl/g dw (i.e., ∼11% of EOCl in the CP extract) and therefore would contribute considerably less than CPs. Furthermore, these compounds might not be extracted with the CP extraction method employed here. Finally, it is germane to note that EOX measurements determined using extraction methods other than those employed in the present work may generate different results. Nevertheless, these measurements represent a first step toward standardizing EOX measurements and will aid in interlaboratory comparison and external quality control.
Halogen Mass Balance in Sewage Sludge. PFAS. Six out of 16 target PFAS were detected in the sludge (Figure 3, Table S16). PFOS was the most prevalent target (7.2 ± 0.6 ng/g dw [sum isomers], recovery-corrected), which is consistent with prior measurements from the same location in 2019 and 2020 (7.7 and 11 ng/g dw, respectively). 30,33 The sum fluorine equivalent concentration from target PFAS  Table S1.

Environmental Science & Technology
pubs.acs.org/est Article measured here (6.4 ± 1.0 ng F/g dw, nonrecovery corrected) accounted for a mere 2% of the 304 ± 116 ng F/g dw in sludge EOF (Figure 4), indicating that the vast majority of the organofluorine is unidentified. Previous reports of EOF concentrations in sludge from the same WWTP showed considerable variability, (i.e., 535−1270 ng F/g reported for samples from 2004 to 2015, 27 154 ng F/g from 2016, 31 and ∼650 ng F/g from 2017 32 ), albeit with slightly different procedures for extraction and quantification. For example, while all three studies used a similar methanolic extraction procedure, two included an alkaline digestion step, two used an ion-pair cleanup, 27,32 and one used a weak anion exchange cleanup. 31 Moreover, two were extraction-blank sub-tracted, 27,31 and one was combustion-blank subtracted. 32 This variability highlights the need for reference materials and a more standardized procedure. 26 The report by Karrman et al. 32 analyzed sludge from the same WWTP and included an expanded target list of 78 PFAS, which despite the inclusion of more PFAS, could only characterize <10% of EOF. The main contributors of EOF in sludge were polyfluoroalkyl phosphate esters (PAPs), perfluooctane sulfonamidoacetates (FOSAAs), fluorotelomer carboxylic acids (FTCAs), and legacy PFCAs and PFSAs. 32 Other possible contributors to EOF include fluorinated pharmaceuticals, pesticides, and polymers. Fluoropharmaceuticals and fluoropesticides typically only have a few fluorine  Table S16.

Environmental Science & Technology
pubs.acs.org/est Article atoms, but at high concentrations, these substances may contribute significantly to the fluorine mass balance. Roughly 20% of commercial pharmaceuticals contain one or more fluorine atoms. 47 Fluorine has been introduced in active ingredients for crop protection during recent decades, and organofluorine represents a growing family of commercial agrochemicals. 48 In fact, nearly 70% of all pesticides introduced into the global market from 2015 to 2020 contained fluorine. 49 Another possibility are fluorinated polymers, which are known to be used throughout consumer products, including cosmetic products, 50 and textiles. 51 Letcher et al. 52 and Fredriksson et al. 53 previously observed side-chain fluorinated polymers in sewage sludge at elevated concentrations, but only by using a substantially different extraction procedure (i.e., acetone/ hexane) than that employed here (MeOH). Thus, the contribution of fluorinated polymers to EOF determined here remains unclear.
CPs. MCCPs (1127 ± 203 ng/g dw) made up the vast majority (i.e., 62%), of CPs in sludge, followed by LCCPs (502 ± 114 ng/g dw; 28%) and SCCPs (181 ± 21 ng/g dw; 10%) ( Figure 3, Table S16). Very short-chain CPs (C 6−9 H x Cl y ) were also observed (5.5 ± 0.5 ng/g dw) but accounted for only 0.3% of the ∑CP concentration ( Figure 3). The concentrations of SCCPs and MCCPs reported here aligned well with a recent report in Swedish sludge, 30 and the congener profile from the present work was very similar to the profile of settled dust from Norway (Table S14). 38 Prior CP measurements in Swedish sewage sludge between 2002 and 2010 have reported higher concentrations of SCCPs, MCCPs, and LCCPs (median concentrations were 1100, 3800, and 31000 ng/g dw for SCCPs, MCCPs, and LCCPs, respectively), 54 possibly due to changes in production of specific chain lengths over the last 20 years, 55 and/or differences in instrumental and quantification techniques. 56,57 Conversion to chlorine equivalents revealed that ∑CP concentrations (951 ± 169 ng Cl/g dw) accounted for nearly all (i.e., 92%) of the measured EOCl CP-extraction concentration (1030 ± 417 ng Cl/g; Figure 4). While the unidentified EOCl CP-extraction was not statistically significant, we cannot rule out the occurrence of relevant organochlorine compounds which appear minor when considered on a chlorine equivalent basis, yet are much more significant when considering molar concentrations. Many pharmaceuticals, for example, although usually more polar and might not be present in this particular extract, contain chlorine atoms (even more so than fluorine and bromine 58 ) and may occur at elevated molar concentrations in sludge. Furthermore, drinking water facilities use a small amount of chloramine for disinfection and could therefore contribute with chlorinated disinfection byproducts (DBPs). 59 Chlorinated OPEs. Following blank subtraction, recoverycorrected TCIPP and TDCPP concentrations were 60.3 ± 8.7 and 4.2 ± 0.5 ng/g dw, respectively ( Figure 3, Table S16). The EOCl HFR-extraction for this fraction was 634 ± 208 ng Cl/g ( Figure 4; Table 1). Due to low recoveries, the ∑Cl Cl-OPEs of the nonrecovery-corrected fraction could not be calculated as concentrations from all three Cl-OPEs were below quantification limits. OPEs have been analyzed previously in Swedish WWTP sludge, and concentrations reported here are lower than measured before (see overview in Table S15).
Possible sources of unidentified EOBr are numerous. For example, although more polar, hydroxylated, and methoxylated PBDEs (OH-and MeO-PBDEs) could end up in the same extracts as parent BDEs. 60 Polybrominated dibenzo-p-dioxins and dibenzofurans (PBDD/Fs) can form after combustion of PBDEs. 61 Another possible source could be from polymeric flame retardants. 62 Brominated drugs, pesticides, rodenticides, 63 and their metabolites could also contribute. Bromhexine, for example, is used in over-the-counter cough syrups and pastils, and bromocriptine, is used in the treatment of Parkinson's. Both have been detected in influent, effluent, and sludge from Sweden. 64,65 Several brominated herbicides are registered in the Swedish Chemical Agency's pesticide register, e.g., bromacil and bromophos. 66 However, these compounds are more polar and are less likely to be coextracted using the BFR extraction procedure employed here.
Prioritizing Extracts for Nontarget Investigation. Further characterization of samples investigated here using nontargetbased analytical approaches should prioritize extracts with both high EOX concentrations and large fractions of unidentified EOX. In the present work, PFAS extracts (304 ± 116 ng F/g dw EOF; 98% unknown) and HFR Fraction I (149 ± 10 ng Br/g dw EOBr; 37% unknown) and Fraction II (39 ± 4 ng Br/ g dw EOBr; 81% unknown) extracts meet these criteria. Tools such as the total oxidizable precursors assay (TOP; for EOF), 67,68 as well as high resolution mass spectrometry-based nontarget and suspect screening (for EOF and/or EOBr) could provide further insight into the unidentified fractions. 69 With regards to chlorinated extracts, nearly all EOCl in the CP extract was accounted for by target CPs. For HFR Fractions II and IIIb, we expect that CPs also explain the vast majority of EOCl, despite that CPs were not measured in this sample, as they have been found to elute in similar fractions previously. 70 Thus, on the basis of halogen equivalents, the HFR extracts from Fraction IIIb could be considered lower priority (albeit with some caveats, discussed in the following section).
In addition to the aforementioned halogen mass balances, we performed EOX measurements on several extracts where a particular halogen was unexpected (Table 1). For example, sludge CP extracts contained 54.0 ± 22 ng F/g dw, which was surprising considering the extraction procedure is designed to favor neutral, lipophilic substances, suggesting the presence of nonconventional PFAS or other fluorinated substances. Similarly, CP extracts from dust SRM 2585 contained 417 ± 77 ng F/g dw. There is a paucity of data on the occurrence of lipophilic organofluorine, with the exception of a single report of high EOF in the blubber of a Greenland Killer whale, 20 and a more recent article predicting the environmental occurrence of mixed halogenated n-alkanes (including those containing fluorine). 71 Similar observations for EOBr (e.g., 2340 ± 433 ng Br/g dw observed in CP extracts of SRM 2585) are also of interest but perhaps less remarkable since a wide range of lipophilic/nonpolar and lipophobic/polar organobromines are known. Thus, nontarget investigations of organofluorines in  73 and biomagnification factors, 74 but is unique in that it identifies specific halogen class(es) and extraction technique(s), rather than a specific feature of a mass spectrum. However, despite showing considerable potential for widespread adoption and applications to a wider suite of matrices, there remains a scarcity of reference material necessary for standardization of halogen mass balance methods. While this work offers some of the first measurements of TX and EOX in publicly available reference materials, consensus values are needed across a much wider suite of matrices.
The halogen mass balance approach also has several limitations that researchers should be aware of prior to adopting it as a prioritization strategy. First and foremost, quantifying on a halogen-equivalent basis reduces the importance of substances with a low number of halogen atoms (or none at all), which may nevertheless exist at high concentrations and/or be important toxicologically. Halogenated pharmaceuticals, for example, are designed to be bioactive, yet typically only contain a few halogen atoms. In a similar vein, we cannot rule out that organohalogens occurring naturally in the terrestrial and marine environment (which may or may not be toxicologically relevant) 76 may contribute to the halogen mass balance.
A fundamental limitation of halogen mass balance experiments is that EOX is method specific; i.e., different concentrations of EOX will be obtained using different methods since some relevant OHCs will be removed during extraction. 17 For this reason, EOX should not be confused with "total" organic halogen, the latter of which can be conceptualized but is difficult to determine experimentally, since it requires assumptions about the method's ability to extract all OHCs. While measurement of inorganic halogen concentrations and subtracting these from TX has been employed for "total" organohalogen determination in consumer products, 75 this assumes that all inorganic halogens exist as halides (i.e., X − ), which is unlikely for sludge. In comparison, EOX assumes that other inorganic halogens (if present) behave similarly to X − and are removed during the extraction procedure. While it could be argued that a more holistic extraction method with fewer cleanup steps, or performing the halogen mass balance after each cleanup step, could capture a large fraction of the "total" organic halogen, these approaches would inevitably compromise the quality of targeted OHC data due to coextraction of matrix interferences. For these reasons, we argue that pairing EOX and targeted analyses on the same extract is the most appropriate and practical approach to performing a halogen mass balance.
Additional information for PFAS, CP, and HFR extraction procedures; schematic illustration of the WWTP; extraction workflow; HFR cleanup workflow; overview target PFAS, CPs, and HFRs; RTs, precursor, quantitative and qualitative ions, ISs, and LOQs for PFAS; spike/recovery results for PFAS, PFAS in NIST sludge material; LOQs for CPs, spike/recovery results for CPs; CPs in NIST dust material; RTs, precursor and quantitative ions, ISs and LOQs for HFRs; IS yields for HFRs; spike/recovery results for HFRs; eluent gradient for EOX analysis; results for NaF spike/recovery experiment; combustion efficiencies for CIC; CP homologue profile; OPE concentrations in literature; mass balance overview (PDF)