The Transformation of Per-and Polyfluoroalkyl Substances in the Aquatic Environment of a Fluorochemical Industrial Park

: Jiangsu High-Tech Fluorochemical Industrial Park in Changshu City, Jiangsu Province, is the largest fluorochemical industrial park in Asia. The occurrence of per-and polyfluoroalkyl substances (PFASs) in surface water and widespread local plants was investigated in Jiangsu High-Tech Fluorochemical Industrial Park. Thirty-two target PFASs were detected in dissolved-phase, particle-phase and plant samples. The concentrations of total PFASs ranged from 1650 to 8250 ng/L in the dissolved-phase samples, 132 to 6810 ng/g dw in the particle-phase samples and 25.8 to 9460 ng/g dw in different plant tissues. Perfluorohexanoic acid (PFHxA), Perfluorooctanoate acid (PFOA) and 6:2 fluorotelomer carboxylic acid (6:2 FTCA) were predominant PFASs and contributed 80 − 91% to Σ PFAS in water samples. A total of 67 emerging PFAS were identified in all samples using nontargeted analysis. Typha orientalis showed better accumulation ability, with an average Σ PFAS concentration of 3450 ng/g dw and the highest root concentration factor (RCF) of 171. Typha orientalis , Eichhornia crassipes and Alternanthera sessilis have potential for use in PFAS phytostabilization.


Introduction
Per-and polyfluoroalkyl substances (PFASs) are a class of man-made persistent organic pollutants (POPs) with the general chemical formula C n F 2n+1 -R, where at least one hydrogen atom on a carbon atom in the carbon chain is replaced by a fluorine atom [1].Due to the different number of carbon atoms and functional groups attached, the number of PFASs is more than 16,000 and is still increasing [2,3].Previous research has demonstrated that PFASs exhibit environmental persistence and bioaccumulation, posing a significant threat to human health, including hepatotoxicity, immunotoxicity, and potential carcinogenicity [4][5][6][7][8][9][10].Recognizing these concerns, 3M, the world's largest PFAS producer, initiated the voluntary phase-out of perfluorooctane sulfonate (PFOS) and its products in 2000 [11].From 2009 to 2019, with the restriction of PFOS, perfluorooctanoic acid (PFOA) and perfluorohexanesulfonic acid (PFHxS) under the Stockholm Convention on Persistent Organic Pollutants, short-carbon-chain PFASs and emerging PFASs emerged to meet the huge market demand [12][13][14][15].Emerging PFAS compounds contain one or more ether bonds, -CH 2 groups, or phosphate (phosphonate) groups and are more easily biotransformed or chemically transformed, posing a potential ecological risk [16].Suspect and nontarget screening methods are used to identify emerging PFASs due to the lack of reference standards for emerging PFASs, and there have been many studies using these methods to discover hundreds of emerging PFAS congeners in environmental and biological samples worldwide [17][18][19][20][21][22].
Water plays a pivotal role as a primary carrier for PFASs because of their water solubility, facilitating the transport of PFASs across diverse environmental media and biological populations.This phenomenon substantially contributes to the widespread distribution of PFASs.Surface water surrounding Zhangjiang High-Technology Park in Shanghai exhibited ΣPFAS concentrations ranging from 115 to 600 ng/L [23].Aquatic environments surrounding fluorochemical industrial parks are often characterized by high levels of PFAS contamination due to the direct discharge of fluorochemicals.The concentrations of total perfluoroalkyl acids (PFAAs) ranged from 4570 to 11,900 ng/L in the lakes surrounding a fluorochemical production base in Wuhan, China [24].Aquatic plants are crucial vectors for the transfer of contaminants from the aquatic environment to the food chain.In hydroponic experiments, aquatic plants displayed potential to remove PFASs from the aquatic environment [25].Greger et al. treated leakage water from a landfill with plants, and the results showed that the removal capacity of aquatic plants (42.1%) was higher than that of terrestrial plants (7.7%) [26,27].The aquatic plant removal rate of PFOA and PFOS was more than 82% after treatment for 15 days in the study by Chen et al. [28].Studies have reported that aquatic plants are effective species for PFAS phytoremediation given their high bioaccumulation (RCF > 1) and translocation capacity (transport factor (TF) > 1) [29,30].Alternanthera sessilis was considered a good phytoremediator of C4-C7 PFCAs in surface water, with mean RCFs and TFs greater than one [31].To date, research on PFAS uptake and distribution in plants has focused on terrestrial plants and agricultural species.Notably, conclusions obtained from soil environments may vary when extrapolated to aquatic ecosystems, due to the differences in environmental factors such as temperature, pH and dissolved oxygen [32].
Jiangsu High-Tech Industrial Park is located north of Changshu City, Jiangsu Province, China, adjacent to the Yangtze River.The park has focused on the development of the fluorine material industry since 1999.It is the largest fluorine material production, import and export base in Asia, and it is a fluorine material industrial park with the longest fluorine industrial chain and the highest agglomeration degree in the world [33].Some scholars have conducted research at Jiangsu High-Tech Fluorochemical Park [33][34][35][36][37][38].The ΣPFAS concentrations were 276 and 120 ng/g dw in tree leaf and bark samples in Jiangsu High-Tech Fluorochemical Industrial Park [35].In vegetable and fruit samples collected from downstream regions, the maximum values were 11.5 and 10.5 ng/g ww [34].However, few studies have investigated the distribution and transport of PFASs within aquatic plants growing in contaminated sites.
In this study, an investigation was conducted at a fluorochemical industrial park in Changshu, China.In total, 32 target PFASs were quantified and 67 suspect and nontarget PFASs were identified.The objectives of this study were to (1) investigate the concentrations and distribution of PFASs in the surface water and aquatic plants in the fluorochemical industrial park; (2) identify emerging PFASs by suspect and nontarget screening; and (3) explore the distribution and transport of PFASs within aquatic plants.

Sampling Campaign
All samples were collected from Jiangsu High-Tech Fluorochemical Industrial Park in October 2022.The sampling sites and main fluorochemical factories are provided in Figure 1.Surface water samples were collected at sites CS1-CS6 (Table S1).Plant samples were widespread local species and collected at sites CS1 (Alternanthera sessilis), CS4 (Eichhornia crassipes) and CS6 (Typha orientalis).
Surface water samples were collected in 1L Polypropylene (PP) bottles.Plant samples were stored in PP plastic bags.Upon returning to the laboratory, the plant samples were washed with tap water and Milli-Q water and separated into roots, stems, and leaves.For further extraction, the plant tissues were freeze-dried and milled into powder.All the samples were frozen at −20 °C before pre-treatment.

PFAS Extraction
Surface water samples were filtered through 0.7 µm glass fiber filters (47 mm, Whatman, Maidstone, UK, pre-combusted at 450 °C for 6 h) to separate the particle phase from water.The dissolved phase was spiked with 2 ng of IS compounds and extracted with Oasis WAX cartridges (30 µm, 6 mL, 150 mg, Waters Corporation, Milford, MA, USA).The cartridges were pre-conditioned with 0.2% ammonia methanol solution, 4 mL of methanol, and 4 mL of Milli-Q water.Then, 10 mL of 0.2% ammonia methanol solution was used for eluting.The eluent was concentrated to 1 mL under a gentle stream of nitrogen gas (>99.999%) and transferred to a Supelclean Envi-Carb cartridge (1 mL, 100 mg, Merck, Darmstadt, Germany) for further purification.Finally, 200 µL of eluent was stored in a sampler vial at −20 °C before analysis.After freeze-drying, the particle-phase samples were extracted with 10 mL of 0.2% ammonia methanol solution by sonication for 20 min and then centrifuged (3000 rpm) for 5 min.The procedure was repeated twice, and the combined extracts were processed in the same way as the dissolved-phase samples.Surface water samples were collected in 1L Polypropylene (PP) bottles.Plant samples were stored in PP plastic bags.Upon returning to the laboratory, the plant samples were washed with tap water and Milli-Q water and separated into roots, stems, and leaves.For further extraction, the plant tissues were freeze-dried and milled into powder.All the samples were frozen at −20 • C before pre-treatment.

PFAS Extraction
Surface water samples were filtered through 0.7 µm glass fiber filters (47 mm, Whatman, Maidstone, UK, pre-combusted at 450 • C for 6 h) to separate the particle phase from water.The dissolved phase was spiked with 2 ng of IS compounds and extracted with Oasis WAX cartridges (30 µm, 6 mL, 150 mg, Waters Corporation, Milford, MA, USA).The cartridges were pre-conditioned with 0.2% ammonia methanol solution, 4 mL of methanol, and 4 mL of Milli-Q water.Then, 10 mL of 0.2% ammonia methanol solution was used for eluting.The eluent was concentrated to 1 mL under a gentle stream of nitrogen gas (>99.999%) and transferred to a Supelclean Envi-Carb cartridge (1 mL, 100 mg, Merck, Darmstadt, Germany) for further purification.Finally, 200 µL of eluent was stored in a sampler vial at −20 • C before analysis.After freeze-drying, the particle-phase samples were extracted with 10 mL of 0.2% ammonia methanol solution by sonication for 20 min and then centrifuged (3000 rpm) for 5 min.The procedure was repeated twice, and the combined extracts were processed in the same way as the dissolved-phase samples.
Plant samples were extracted by the method described by Zhang et al. with a few minor modifications [39].In total, 10 mL of NaOH (0.4 M) was added to a 50 mL PP tube containing the freeze-dried plant tissue samples.The mixture was spiked with 2 ng of IS compounds and refrigerated at 4 • C overnight.Then, 2 mL of TBAHS (0.5 M) and 4 mL of Na 2 CO 3 /NaHCO 3 buffer (0.25 M, pH10) were added into the mixture for extraction.After vortexing, 10 mL of MTBE was added and the mixture was shaken vigorously for 20 min.The organic and aqueous layers were separated by centrifugation (3000 rpm, 10 min), and the MTBE layer was transferred to another PP tube.The aqueous phase was further extracted twice with MTBE.The combined organic phase was evaporated under nitrogen gas and then reconstituted in 1 mL of methanol.After diluting with 9 mL of water, the sample was extracted with Oasis WAX cartridges and eluted by 4 mL of methanol and 4 mL of 0.2% ammonia methanol solution.The eluent was reduced to 1 mL under a gentle stream of nitrogen gas and then filtered through 0.22 µm nylon organic phase filters (13 mm, CNW, Berlin, Germany) before being transferred into a sampler vial.

Target Analysis
Thirty-two PFASs were measured in this study; details of the mass parameters for the compounds are given in Table S2.The quantification of PFASs was performed using ultra-high-performance liquid chromatography coupled with triple-quadrupole tandem mass spectrometry (UPLC-MS/MS, SHIMAZU, Kyoto, Japan) under multiple reaction monitoring (MRM) modes with negative electrospray ionization (ESI).A ZORBAX Eclipse Plus C18 analytical column (2.1 × 100 mm, 1.8 µm, Agilent, Santa Clara, CA, USA) was used at 50 • C. The mobile-phase constituents were Milli-Q water (A) and methanol (B) at a flow rate of 0.3 mL/min.The gradient was as follows: t = 0.1 min, 25% B; t = 5 min, 75% B; t = 10 min, 100% B. The injection volume was 1 µL.
Samples were analyzed by intelligent data acquisition mode in AcquireX (Thermo Fisher Scientific, Waltham, MA, USA), which was based on the data-dependent acquisition (DDA) mode.AcquireX was used to acquire abundant MS2 information, and the data processing workflow was run by Compound Discoverer 3.3 (CD 3.3, Thermo Fisher Scientific, Waltham, MA, USA).For suspect screening, the raw data were pre-filtered following suspect lists, which combined the US EPA CompTox Chemistry Dashboard [40], the Norman Suspect List Exchange [41], and the Nist Suspect List of Possible Per-and Polyfluoroalkyl Substances (PFASs) [42].When the suspect lists were unavailable, emerging PFAS compounds were discovered by extracting peaks with feature fragment ions [43][44][45].For further confirmation in nontarget screening, the peak shape, peak response intensity and rationality of retention time were checked to comprehensively judge whether the feature was retained [36].The confidence level (CL) of 2-4 was determined according to the study of Schymanski et al. [46].Semi-quantitative analysis was conducted to determine the estimated concentrations of nontarget PFASs using target PFASs with similar structures [36].More details can be found in Text S2.
2.6.Quality Assurance and Quality Control (QA/QC) Teflon and polytetrafluoroethylene (PTFE) materials were avoided during the study to minimize the background contamination.The correlation coefficient of the calibration curve consisted of a concentration gradient (0.5, 1, 2, 5, 10, 20, 50, 100 and 200 ng/mL) of PFAS standards and was at least 0.999 for all target analytes.Procedure blanks were prepared using Milli-Q water and standard silica sand; they were routinely analyzed to check for contamination during PFAS extraction.The method spike recovery experiment was verified by adding a target PFAS to blank matrices including Milli-Q water and standard silica sand, and the matrix spike recovery experiment was verified by adding a target PFAS to plants with low background (Table S3).The method spike recoveries of all PFASs were in the range of 60-120%.Matrix effects were deemed not obvious, with the matrix effect recovery rates within the acceptable range of 70% to 130%.The limit of detection (LOD) and the limit of quantification (LOQ) were defined as the concentrations giving signal-to-noise ratios of 3 and 10, respectively (Table S3).The average spike recoveries of the IS compounds in the instrumental analysis are provided in Table S4.

Statistical Analysis
The partition coefficient, K d , of PFASs between particles and dissolved phases was calculated by Equation ( 1): C particle phase : the concentrations of PFASs in the particle phase (ng/g dw); C dissolved phase : the concentration of PFASs in the dissolved phase [24,47].
The root concentration factor (RCF), stem transfer factor (TF s ) and leaf transfer factor (TF l ) were calculated by Equations ( 2)-( 4): C root : the concentrations of PFASs in the roots of plants (ng/g dw); C stem : the concentrations of PFASs in the stems of plants (ng/g dw); C leaf : the concentrations of PFASs in the leaves of plants (ng/g dw) [48,49].

Concentration and Distribution of PFASs in Surface Water
The concentrations of PFASs in the dissolved phase and particle phase are summarized in Table 1.In the dissolved phase, 29 PFASs were detected among the 32 target PFASs analyzed in this study.The ΣPFAS ranged from 1650 to 8250 ng/L, with an average value of 3990 ng/L.This was higher compared to a 2013 study in Jiangsu High-Tech Fluorochemical Industrial Park, which reported a mean concentration of 2800 ng/L [33].This indicated the influence of fluorochemical plants on the surrounding surface water over time.In the particle phase, the ΣPFAS ranged from 113 to 6800 ng/g dw, with a mean concentration of 1140 ng/g dw.The particle levels (g dw) in the particle phase are given in Table S5.All Perfluorocarboxylic acids (PFCAs) in this study were detected at a frequency of 100%, except perfluoropentanoic acid (PFPeA) (33%) and perfluorotridecanoic acid (PFTriA) (67%).Meanwhile, in the particle phase, C4-C14 PFCAs were detected with detection frequencies (DFs) of 100%, followed by PFBA (83%) and PFPeA (67%).The mean (median) concentrations of perfluorononanoic acid (PFNA), perfluorodecanoic acid (PFDA), perfluoroundecanoic acid (PFUnDA), perfluorododecanoic acid (PFDoDA), PFTriA and perfluorotetradecanoic acid (PFTreA) with longer chain lengths (C > 8) were in the range of <0.290-42.6 ng/L (<0.290-5.26ng/L), which were much lower than those of PFCAs with short chain lengths (C ≤ 6).The mean (median) levels of perfluorinated sulfonic acids (PFSAs) were in the range of <0.228-9.15ng/L (<0.228-10.86 ng/L) and relatively lower than those of PFCAs.For PFSAs in the particle phase, PFPeS and PFOS were 100% detected, with average concentrations of 0.0900 ng/g dw and 2.04 ng/g dw.Long-chain PFTreA and PFDS, not detected in the dissolved phase, were found in the particle phase at concentrations ranging from 0.240 to 2.25 ng/g dw and 0.440 to 1.12 ng/g dw, respectively.This suggested a stronger affinity of long-chain PFCAs and PFSAs to the particle phase [50].For emerging PFAS, 6:2 FTCA, 6:2 fluorotelomer sulfonic acid (6:2 FTSA), 6:2 chlorinated polyfluorinated ether sulfonic acid (6:2 Cl-PFESA), hexafluoropropylene oxide dimer acid (HFPO-DA) and trimer acid (HFPO-TA) were frequently detected with DFs of 100%, and the mean (median) value of them ranged from 0.29 to 1840 ng/L (0.300-1060 ng/L).As substitutes for PFOA, HFPO-DA and HFPO-TA were present in total concentrations ranging from 16.4 ng/L to 601 ng/L (with an average of 179 ng/L) and 76.8 ng/L to 731 ng/L (with an average of 240 ng/L), respectively.The PFAS levels in this study were classified as moderate compared to other regions.The average concentrations of ΣPFAS in fluorinated industrial parks in Zibo, China and Fuxin, China were 354,000 ng/L and 3050 ng/L [36].As for the river water collected upstream and downstream from a fluorochemical production plant in the Netherlands, the average concentrations of ΣPFAS occurred in a range from 38.0 to 876 ng/L [51].Figure 2a,b show that PFHxA, PFOA and 6:2 FTCA were the predominant PFASs in the dissolved phase and particle phase in Jiangsu High-Tech Fluorochemical Industrial Park, with both contributing 84% to ΣPFAS.This was consistent with the research by Song et al. [36].In the dissolved phase, PFHxA accounted for the highest proportion (45%) at site CS6.PFOA accounted for the highest proportion at site CS3 (27%) and CS6 (27%), which was the most widely used PFAS during the emulsion polymerization of PTFE [52].This may be because PFHxA has been used by Arkema Factory to produce Fluorine Kautschuk material (FKM) since 2008 [53].Chemours (part of the DuPont fluorochemicals division) and the fluoropolymer plant of Arkema were reported to use PFOA to manufacture PTFE [52].It was shown that 6:2 FTCA constituted 18-59% and 22-53% of the total PFASs in two phases, which may have originated from industrial wastewater discharge during its production process, and the perfluorocarbon C6 alkyl structure products produced by the DuPont Branch through the telomerization reaction may be important precursors of 6:2 FTCA [36,54].Emerging PFASs, particularly 6:2 FTSA, HFPO-DA and HFPO-TA, dominated the dissolved phase and particle phase with total concentrations of 2520 ng/L and 1630 ng/g dw.These three compounds can be used as processing aids in the manufacture of fluorinated polymers, including polyvinylidenefluoride (PVDF) and PTFE, which are the main products of Arkema, Chemours and Sanaifu [53,55].The concentrations of ΣPFAS were comparatively higher at sampling sites CS3, CS4 and CS5 compared to other sampling points.The highest concentration was detected at CS4, which represented the confluence of two rivers.The mixing of two rivers may lead to a release of PFASs in the sediment, thus resulting in higher concentrations at the confluence [56].The concentrations of ΣPFAS at CS1 and CS2 were 1790 ng/L and 1910 ng/L; this may have been due to the well-developed textile and apparel industries and the proximity of some sampling points to the fluorochemical factory [57].
The occurrence of nontarget PFASs in different media is shown in Figure 3 based on semi-quantitative calculations.A total of 67 emerging PFASs were identified in all samples using nontargeted analysis (Tables S7 and S8).There were nine classes in the CF 2 Kendrick mass defect plot (Figure S1), where same line implied homologues with CF 2 units.Twenty-seven emerging PFASs were identified as Level 2, including ultra-short-chain PFCAs, unsaturated perfluorocarbonic acids (UPFCAs), perfluoroether carbonic acids (PFE-CAs), hydrogenated perfluorocarbonic acids (H-PFCAs), chlorinated perfluorocarbonic acids (Cl-PFCAs), hydrogenated perfluorosulfonic acids (H-PFSAs), ketone perfluorosulfonic acids (Ke PFSAs) and perfluoroether sulfonic acids (PFESAs), totaling eight classes.Figures S2-S17 provide structures and MS/MS spectra evidence of each class of PFAS.The estimated concentrations of nontarget PFASs (Level 2) in surface water (dissolved phase and particle phase) were in the ranges of 61.3-649 ng/L and 3160-62,600 ng/g dw (Figure 3).The total estimated concentrations of nontarget PFASs exceeded those of target PFASs in 15 out of 22 samples, indicating the severe emerging PFAS pollution.In the dissolved phase, UPFCAs, H-PFCAs and Ether PFCAs were the main pollutants, accounting for 28%, 28% and 23% of the total nontarget PFAS concentration, while H-PFSAs, UPFCAs and H-PFCAs dominated in the particle phase with proportions of 53%, 31% and 10%.H-PFCAs could be an alternative polymerization agent of PVDF production or an impurity in the oligomerization process using ammonium salt of C 9 HF 17 O 4 as the main active ingredient [36].In this study, a small amount of C 9 HF 17 O 4 was identified (2-3% of the total nontarget PFAS concentration), so H-PFCAs may be mainly released from the production of PVDF.
Water 2024, 16, x FOR PEER REVIEW 8 of 16 In this study, a small amount of C9HF17O4 was identified (2-3% of the total nontarget PFAS concentration), so H-PFCAs may be mainly released from the production of PVDF.Kd values varied greatly among different sampling points, which could have been due to the dynamic aquatic ecosystem, in which it is difficult to achieve a steady-state equilibrium in the environment [47,58].Therefore, the mean value of log Kd was used to In this study, a small amount of C9HF17O4 was identified (2-3% of the total nontarget PFAS concentration), so H-PFCAs may be mainly released from the production of PVDF.Kd values varied greatly among different sampling points, which could have been due to the dynamic aquatic ecosystem, in which it is difficult to achieve a steady-state equilibrium in the environment [47,58].Therefore, the mean value of log Kd was used to

Partitioning of PFASs in Dissolved Phase and Particle Phase
K d values varied greatly among different sampling points, which could have been due to the dynamic aquatic ecosystem, in which it is difficult to achieve a steady-state equilibrium in the environment [47,58].Therefore, the mean value of log K d was used to represent the partition coefficient in this study.As in shown in Figure 4, the average log K d of PFCAs ranged from 2.04 to 3.41, exhibiting a positive linear correlation (r = 0.7117, p < 0.05) with increasing carbon chain length (C4-C14) from PFBA to PFTriA.Meanwhile, the average log K d of PFSAs (2.00-3.34)showed a positive linear correlation (r = 0.8333, p < 0.05) with increasing carbon chain length (C5-C8) from PFPeS to PFOS.This observed trend aligns with previous research findings [47,59].The particle phase may influence the distribution of PFASs in aquatic environments.Long-chain PFCAs and PFSAs were more likely to interact with the particle phase, promoting their long-distance transport and deposition in sediments via water movement.PFECHS had a relatively high average log K d value (3.42), suggesting that the hexagonal chemical ring structure of PFECHS enhances its association with the particle phase [47].

Absorption of PFASs in Aquatic Plants
Figure 2c,f show the concentrations and distributions of PFASs in aquatic plants.A total of 32 PFASs were detected in aquatic plant samples from Jiangsu High-Tech Fluorochemical Industrial Park, of which, 10 PFASs had DFs of 100%, including PFBA, PFHxA, PFHpA, PFOA, PFNA, PFUnDA, PFDA, PFUnDA, PFOS, 8:2 FTCA, 6:2 FTSA and HFPO-DA.PFOA and PFHxA were the main substances in plant samples, with mean concentrations in three species of 415 ng/g dw and 159 ng/g dw, indicating that plants may be affected by river pollution.Short-chain PFAAs (C4-C7) were more abundant than longchain PFAAs (C > 8) in aquatic plants.On average, short-chain PFAAs were found 24 times more than long-chain PFAAs.This aligned with previous research indicating that there was a higher uptake rate of PFASs with shorter-chain carbons in plants [34,60].The main emerging PFASs were 8:2 FTCA and HFPO-DA, with average concentrations of 23.6 ng/g
RCF was used to estimate a plant's ability to accumulate PFASs from water.RCF > 1 was considered to indicate effective accumulation [29,64].As is shown in Figure 5 and Table S7, RCFs were in the range of 0.00460-1680, 0.000234-60.4 and 0.00665-61.8 in Typha orientalis, Eichhornia crassipes and Alternanthera sessilis, with average RCFs of 5.54, 0.0314 and 0.0352.Among the three plants, Typha orientalis had a relatively strong root bioaccumulation ability and was most effective in taking up PFHpA, PFDoDA, PFTriA, PFPeS, PFHxS, PFHpS, PFOS, 6:2 FTSA and PFECHS with higher RCFs.PFCAs followed a "U"-shaped trend of RCFs decreasing as chain length increased from PFBA to PFOA, and increasing from PFOA to PFTriA, except slightly increased RCFs of PFHpA in Typha orientalis and PFPeA in Eichhornia crassipes.In contrast, RCFs increased with increasing chain lengths of PFSAs (Figure 5).Similar phenomena have been observed in multiple studies [65][66][67][68].According to the previous indoor exposure experiments, the root accumulation of shortchain PFCAs (C < 7) and PFSAs (C < 6) may be due to their small molecular size and strong hydrophilicity, making them easy to be absorbed by plants' roots along with water.For long-chain PFCAs and PFSAs, this may be due to their strong hydrophobicity, making them easy to be adsorbed on the root surface [68][69][70].But in this study of an actual aquatic environment, there were more complex mechanisms; the growth environment and plant species may have been possible influencing factors.For emerging PFASs, 6:2 Cl-PFESA had a higher RCF than PFOS in Eichhornia crassipes, which may have been due to the hydrophobic group Cl reducing the hydrophilicity of the PFASs, thereby weakening their migration ability [71].Based on previous studies, nontarget PFASs may be by-products of industrial production processes and degradation products formed by the biological transformation of PFAS precursors in plants [36,62,63].RCF was used to estimate a plant's ability to accumulate PFASs from water.RCF > 1 was considered to indicate effective accumulation [29,64].As is shown in Figure 5 and Table S7, RCFs were in the range of 0.00460-1680, 0.000234-60.4 and 0.00665-61.8 in Typha orientalis, Eichhornia crassipes and Alternanthera sessilis, with average RCFs of 5.54, 0.0314 and 0.0352.Among the three plants, Typha orientalis had a relatively strong root bioaccumulation ability and was most effective in taking up PFHpA, PFDoDA, PFTriA, PFPeS, PFHxS, PFHpS, PFOS, 6:2 FTSA and PFECHS with higher RCFs.PFCAs followed a "U"shaped trend of RCFs decreasing as chain length increased from PFBA to PFOA, and increasing from PFOA to PFTriA, except slightly increased RCFs of PFHpA in Typha orientalis and PFPeA in Eichhornia crassipes.In contrast, RCFs increased with increasing chain lengths of PFSAs (Figure 5).Similar phenomena have been observed in multiple studies [65][66][67][68].According to the previous indoor exposure experiments, the root accumulation of short-chain PFCAs (C < 7) and PFSAs (C < 6) may be due to their small molecular size and strong hydrophilicity, making them easy to be absorbed by plants' roots along with water.For long-chain PFCAs and PFSAs, this may be due to their strong hydrophobicity, making them easy to be adsorbed on the root surface [68][69][70].But in this study of an actual aquatic environment, there were more complex mechanisms; the growth environment and plant species may have been possible influencing factors.For emerging PFASs, 6:2 Cl-PFESA had a higher RCF than PFOS in Eichhornia crassipes, which may have been due to the hydrophobic group Cl reducing the hydrophilicity of the PFASs, thereby weakening their migration ability [71].

Transport of PFASs in Aquatic Plants
The transport pathways of PFASs in aquatic plants were compound-specific, and there were differences in concentrations among different tissues even within the same

Transport of PFASs in Aquatic Plants
The transport pathways of PFASs in aquatic plants were compound-specific, and there were differences in concentrations among different tissues even within the same species [72].Among the four plant tissues (roots, stems, leaves and inflorescences) of Typha orientalis grown at CS6, the highest total concentration of the target PFAS was observed in the roots (9460 ng/g dw), followed by stems (4170 ng/g dw), leaves (112 ng/g dw) and inflorescences (51.2 ng/g dw).The fast growth rate, large size and significant below-surface components made Typha orientalis more likely to absorb pollutants from the environment [73].It was noted that the highest concentration of PFOS was detected in the roots of Typha angustifolia (8760 ng/g d.w.), accounting for 93% of ΣPFAS.This may have been due to the fact that the root cell structure of Poaceae Barnhart changed with growth, increasing the contact area between their roots and pollutants in the aquatic environment [74].Moreover, previous research had reported that PFOS accumulated largely in the roots (49-96%) of wetland plants [75].In contrast, for Eichhornia crassipes collected from CS4, leaves exhibited the highest total target PFAS concentration (1800 ng/g dw), followed by stems (624 ng/g dw) and roots (259 ng/g dw).It is noteworthy that the concentrations of PFSAs in the roots of Eichhornia crassipes were 5 times higher than those in the leaves, demonstrating a stronger affinity for binding to root cells.This observation aligned with the results of hydroponic experiments conducted by Phung et al. [76].The trend of leaf accumulation was also observed in Alternanthera sessilis samples from CS1, where leaves (140 ng/g dw) accumulated the most target PFASs, followed by roots (48.7 ng/g dw) and stems (25.8 ng/g dw).Across the three plant species examined, stems exhibited a relatively weaker capacity for PFAS accumulation, contributing 1-23% to ΣPFAS.This may have been because when PFASs in the root system migrate through the xylem of the stem via transpiration, stems thus exhibit a relatively weaker capacity to retain these PFASs as primary transport organs compared to other plant tissues [77].In the leaves of Eichhornia crassipes and Alternanthera sessilis, short-chain PFCAs accounted for 58% and 63% of the total PFCAs, indicating that short-chain PFCAs with stronger hydrophilicity were more likely to migrate upwards in floating plants through transpiration.
In this study, TFs were used to measure the ability of plants to absorb and transport pollutants [78].When TF > 1, this indicates that PFASs have the ability to migrate and transport between different environmental media [49].On average, TF s followed the following order: Typha orientalis (31.6) > Eichhornia crassipes (2.64) > Alternanthera sessilis (0.800).Meanwhile, TF l showed a different trend: Eichhornia crassipes (7.83) > Alternanthera sessilis (5.02) > Typha orientalis (1.65).Among PFCAs (C < 7) in floating plants (Eichhornia crassipes and Alternanthera sessilis), the average TF s and TF l were higher than those in PFCAs (C ≥ 8), indicating that these two plants had a stronger ability to transport shortchain PFCAs.For PFSAs, only PFBS and PFPeS had TF s values greater than 1.Two possible explanations for the declining TF s with increasing chain length were that (1) largermolecular-weight PFASs require more energy for transport, leading to reduced mobility; (2) longer carbon chain lengths may hinder their ability to cross the Kjeldahl band, a barrier within the plant [68,79].As an alternative to PFOA, both TF s and TF l of HFPO-DA in Typha orientalis were higher than PFOA (∆TF = 0.860-7.16),and it was the opposite situation in another two species.Similar patterns of TF were also observed in a previous study, where the TF value of HFPO-DA was significantly higher than that of PFOA across the four plants studied, with ∆TF = 0.9-3.2(Phyllanthus urinaria, Justicia procumbens, Eleusine indica, and Aster indicus) [68].This trend may indicate that ether bonds could enhance the migration and transformation of pollutants, and plant species also have an impact on the migration of pollutants in their bodies [78].Intriguingly, the RCF values for PFHxA exhibited minimal variation across Typha orientalis, Eichhornia crassipes and Alternanthera sessilis (0.0300, 0.0200 and 0.0100), despite TF s values differing (139, 36.5, and 0.760).Conversely, the RCF values for PFOS showed a vast difference of 10 3 orders of magnitude between the three plants (1680, 0.840 and 0.180), yet TF s values remained relatively similar (0.0100, 0.710 and 1.10).
These observations highlighted how both plant species and the inherent properties of PFAS compounds themselves influenced plants' abilities to transport these chemicals.
Plant species with both RCFs and TFs greater than one have the potential to be used for phytoremediation [80,81].Although no PFAS phytoremediators were found in this study, PFAS-tolerant species with high BCF (BCF > 1) and low TF (TF < 1) values can be used for the phytostabilization of contaminated sites [82,83].In this study, such plants included Typha orientalis for C6-C8 PFSAs, 6:2 FTSA and PFECHS; Eichhornia crassipes for 8:2 FTSA, PFECHS, FHxSA and FOSA; and Alternanthera sessilis for 8:2 FTCA and FHxSA (Figure 5).More plant species should be studied to find PFAS phytoremediators in actual aquatic environments.Also, more influencing factors corresponding to the translocation and bioaccumulation of PFASs within plant compartments merit further investigation.

Conclusions
In the Jiangsu High-Tech Fluorochemical Industrial Park, PFHxA, PFOA and 6:2 FTCA were the primary contaminants in surface water.The spatial distribution of ΣPFAS in dissolved-and particle-phase samples suggests that fluorochemical plants may be potential sources.Long-chain PFCAs and PFSAs have a higher affinity for interacting with the particle phase.Emerging PFAS contamination within fluorochemical industrial parks has warranted serious attention.Typha orientalis, Eichhornia crassipes and Alternanthera sessilis have the potential for use in PFAS phytostabilization.The uptake and transport of PFASs in plants exhibit species and compound specificity.The root bioaccumulation ability of Typha orientalis was relatively strong.Short-chain PFCAs were more likely to migrate upwards in floating plants (Eichhornia crassipes and Alternanthera sessilis) through transpiration.
Author Contributions: Methodology, data analysis, manuscript drafting and formal analysis, J.H.; writing-review, editing and funding acquisition, Z.Z.; software, J.L.; validation, S.L.All authors have read and agreed to the published version of the manuscript.
Funding: This research was funded by the Natural Science Foundation of Shanghai, 22ZR1427600 and National Natural Science Foundation of China, NSFC41977310.

Figure 1 .
Figure 1.Sampling sites around the fluorochemical industrial park.

Figure 1 .
Figure 1.Sampling sites around the fluorochemical industrial park.

Figure 2 .
Figure 2. Percentages of target PFASs in dissolved phase (a), particle phase (b) and plants (c) in the fluorochemical industrial park.Concentrations of target PFASs in dissolved phase (d), particle phase (e) and plants (f) in the fluorochemical industrial park.

Figure 3 .
Figure 3. Percentages of nontarget PFASs in the dissolved phase (a), particle phase (b) and plants (c) in the fluorochemical industrial park.Concentrations of nontarget PFASs in the dissolved phase (d), particle phase (e) and plants (f) in the fluorochemical industrial park.

Figure 2 .
Figure 2. Percentages of target PFASs in dissolved phase (a), particle phase (b) and plants (c) in the fluorochemical industrial park.Concentrations of target PFASs in dissolved phase (d), particle phase (e) and plants (f) in the fluorochemical industrial park.

Figure 2 .
Figure 2. Percentages of target PFASs in dissolved phase (a), particle phase (b) and plants (c) in the fluorochemical industrial park.Concentrations of target PFASs in dissolved phase (d), particle phase (e) and plants (f) in the fluorochemical industrial park.

Figure 3 .
Figure 3. Percentages of nontarget PFASs in the dissolved phase (a), particle phase (b) and plants (c) in the fluorochemical industrial park.Concentrations of nontarget PFASs in the dissolved phase (d), particle phase (e) and plants (f) in the fluorochemical industrial park.

Figure 3 .
Figure 3. Percentages of nontarget PFASs in the dissolved phase (a), particle phase (b) and plants (c) in the fluorochemical industrial park.Concentrations of nontarget PFASs in the dissolved phase (d), particle phase (e) and plants (f) in the fluorochemical industrial park.

Water 2024 ,
16,  x FOR PEER REVIEW 9 of 16 represent the partition coefficient in this study.As in shown in Figure4, the average log Kd of PFCAs ranged from 2.04 to 3.41, exhibiting a positive linear correlation (r = 0.7117, p < 0.05) with increasing carbon chain length (C4-C14) from PFBA to PFTriA.Meanwhile, the average log Kd of PFSAs (2.00-3.34)showed a positive linear correlation (r = 0.8333, p < 0.05) with increasing carbon chain length (C5-C8) from PFPeS to PFOS.This observed trend aligns with previous research findings[47,59].The particle phase may influence the distribution of PFASs in aquatic environments.Long-chain PFCAs and PFSAs were more likely to interact with the particle phase, promoting their long-distance transport and deposition in sediments via water movement.PFECHS had a relatively high average log Kd value(3.42),suggesting that the hexagonal chemical ring structure of PFECHS enhances its association with the particle phase[47].

Figure 4 .
Figure 4. (a) Log Kd values of PFASs in the water dissolved phase-particle phase system in the fluorochemical industrial park.The linear relationship between the carbon chain length of PFCAs (b) and PFSAs (c) and log Kd values in the water dissolved phase-particle phase system.Red line represented that there were only one log Kd value.

Figure 4 .
Figure 4. (a) Log K d values of PFASs in the water dissolved phase-particle phase system in the fluorochemical industrial park.The linear relationship between the carbon chain length of PFCAs (b) and PFSAs (c) and log K d values in the water dissolved phase-particle phase system.Red line represented that there were only one log K d value.

Table 1 .
Concentrations (ng/L and ng/g dw) of PFASs in water (dissolved-phase and particle-phase) samples.