3.1 Concentration and profiles
Eighteen individual PFASs all were analyzed in park A and B (Fig. 2). All the short chain PFASs (C<8) and F-53B in park A and B all were detected with the rate of 100%. Generally, the contamination status of PFASs in park A was higher than park B. Due to the industrial production of two parks, the total concentration of PFASs (∑PFASs) and the main individual-PFASs concentration exceeds other regions of Taihu lake(Ma et al. 2018), West lake(Lu et al. 2018), Jiaozhou bay (Han et al. 2020)and Bohai sea(Chen et al. 2017) (Table S6). In addition, comparing with other drinking water sources of central eastern China, mainly detected substances (PFBA, PFHxA and PFOA) and ∑PFASs were lower than our studies area(Lu et al. 2018). Some previous studies revealed serious contamination degree in local area. The range of ∑PFASs(15 target substances) in 2018 was 87-18,000 ng·L-1, and PFBA was 16-12,000 ng·L-1, followed by PFBS (4-5,800 ng·L-1) and PFOA with (16-1,400 ng·L-1), respectively(Chen et al. 2018). Then, Bao et al (2020) detected groundwater as local irrigation source that ∑PFASs(10 target substances) has exceed 35 µg L-1, the mainly detected PFASs concentrations of PFBS (32.40 µg L-1) and PFOA (2.47 µg L-1) (Bao et al. 2020)(Table S2).
In park A, ∑PFASs varies from 72.49 ng·L-1 in Ga7 to 68,142.16 ng·L-1 in Ga10 and the mean concentration degree was 9,104.03 ng·L-1. Although park A has been forbade to carry out industrial production and lacks out obvious point pollution for a long time, PFASs were detected with the high concentration due to itself environmental accumulation and stability of persistent organic pollutant (POPs)(Naile et al. 2013, Nascimento et al. 2018, Ojemaye &Petrik 2019). Short chain PFASs concentration was obviously stronger than long chain PFASs. The highest mean concentration of PFBS was detected with 8,734.27 ng·L-1, followed by PFOA(1,563.13 ng·L-1), PFBA(400.41 ng·L-1), PFHxA(136.09 ng·L-1) and PFHxS(106.87 ng·L-1), respectively (Table 1). The larger concentration contribution of C4-C7 (16%-63%) and C8 (34%-83%) were stronger than C>8 PFASs (≤3%) (Fig 3a). Raining activity in surface and ground system play an important impact on the migration of PFASs(Ahrens et al. 2016, Xu et al. 2021). The ground water flow and surface water supply will cause the increasing of PFASs concentration in superficial aquifer while along with time variation, a mass of storm runoff will dilute high PFASs concentration(Li et al. 2020a). PFSAs with functional groups of -SO3H in surface and groundwater emerge a larger concentration ratio (51%-80%) more than PFCAs with -COOH (20%-49%). In especial, PFSAs in Ga2 with 80% concentration contribution whereas PFCAs only accounts for 20% (Fig. S1a).
In park B, the PFASs contamination status in groundwater of WS is far stronger than in DS. The mean concentration was calculated with 801.68 and 714.64 ng·L-1 in WS and DS, respectively. The mean concentration of surface water was analyzed with 136.39 ng·L-1 in WS whereas mean concentration has reached 1,257.60 ng·L-1 in DS. Similarly, PFBS, PFBA and PFOA were also the mainly detected substances in surface and ground water of WS or DS (Table 1). The serve contamination of short chain PFASs in recent years results from their vast using as long chain alternatives. As Bao et al study indicated that PFBA, PFPeA and PFHxA were used to be the alternatives of PFOA and PFBS is the main substitute of PFOS(Bao et al. 2020). This sampling site of Rbb1 locates on the agricultural area, detected with serious individual concentration and ∑PFASs due to long term reception place from agricultural activities, even pesticides and fertilizers using (Fig.1, Fig. 2c and Table S7). Surface water will accelerate the migration of pollutants, and especially the supply also carries PFASs into groundwater(Chiavola et al. 2020, Feng et al. 2020). The PFOS substitute of F-53B was detected with low contamination degree, whereas concentration in surface water was higher in the groundwater. What in mentioned that long chain PFASs (C≥10) that PFDS, PFUdA, PFDoA, PFTeDA, PFTrDA and PFODA were detected certain concentration in all water body of DS without appearing in WS (nd). C4-C7 and C8 of PFASs show higher concentration contribution than C9 of PFASs in the surface and ground water of WS and DS. There was a difference that concentration rate of C4-C7 PFASs in DS (range of 64%-98%) was far lower than in WS (range of 95%-98%) (Fig. 3b,c). Time difference was better reflected on short chain PFASs than long chain PFASs. The reason was that short chain PFASs with greater hydrophobicity and lower lipophilicity will tend to be affected by season variation while long chain PFASs were adsorbed by organic matter and the suspended material in sediment(Ahrens et al. 2016, Gagliano et al. 2020, Xu et al. 2021, Zheng et al. 2020). For function groups, the concentration contribution of PFSAs was far than PFCAs in park B (Fig. S1b). Chemical groups between -SO3H and -COOH make PFASs in water possess the environmental behavior abundantly(Gobelius et al. 2018).
In our study, PFASs concentration have closely negative relationship with chain length in park A and B (Fig. S2). PFASs with unique physicochemical structure tend to dissolve in water body and migrate with a long distance. In addition, the functional groups are also an important factor(Liu et al. 2019). Functional groups difference between PFSAs and PFCAs are limited to affect by the water load variation in different term (Fig S1). Although other non-point sources may contribute to PFASs in surface and groundwater, the park production is the major contamination source. In addition, industrial emission in different time also affects PFASs profiles by discharge of sewage because traditional equipment is hard to removal effectively.
3.2 Spatial distribution
The source of PFASs in the groundwater may be assumed to two pathways. One was that industrial waste water was directly discharged into groundwater system and another was the ground pollutant source by some environmental behavior of the surface runoff, surface water recharge and rainwater deposition, which will carry PFASs and its precursors into groundwater(Ahrens et al. 2010, Bao et al. 2020, Blaine et al. 2014, Cao et al. 2019, Chen et al. 2015, Miranda et al. 2021). The spatial distribution of the pollutants is not only related to its chemical properties and structures but also the environmental factor and geographical condition should be taken into consideration(Catherine et al. 2019, Chen et al. 2015, Choi et al. 2011, Gao et al. 2020). Especially in intensive industry region, river direction, groundwater flow direction and agricultural activities, all closely affect the distribution of PFASs. In park A, mainly detected substances of individual PFASs and ∑PFASs show an increasing tendency along with the groundwater flow direction (Fig. S3). This result indicates that the groundwater flow direction is the key element to the PFASs polluted distribution. In addition, the mainly detected individual-PFASs of PFBA, PFBS, PFHpA, PFHxA, PFHxS, PFPeA and PFOA is the same with the tendency of ∑PFASs (Fig. S3). Different industrial types are the critical reason of causing PFASs distribution but PFASs in water body was also influenced by the upstream waterbody(Lee et al. 2020). Although industrial process is the main source, the groundwater flow direction also possesses a significant function to PFASs distribution and migration in park A. Other studies evidenced that the PFASs not only comes from the point source but also attributes to water mass migration(Lee et al. 2020).
In park B, this distribution of PFASs concentration was limited to be reflected, and the dominating detected substances and ∑PFASs also lack out obviously increasing or decreasing concentration tendency along with the groundwater flow direction. The concentration of PFASs in the downstream water is slightly higher than the upstream along with the groundwater flow direction but this tendency is largely weak (Fig S4). Some sites show abnormal concentration due to the location distribution of the waste water treatment plant (WWTP) and fire station (FS). These places all influence the spatial distribution of PFASs. Current biotreatment in WWTP can remove parts of PFNA, PFHxA, and PFPeA but has limited ability to treat other PFASs(Knutsen et al. 2019, Xu et al. 2021). In addition, FS as important pollution point has been widely reported. Graetz et al. stated that the past firefighting water additives were found to contain the PFASs that could persist in the environment resulting in the potential adverse effects to the biota(Graetz et al. 2020). Importantly, a FS in the park B is the most polluted sources. Daily training and drills will widely release the PFASs precursors into the environment (Gobelius et al. 2017, Graetz et al. 2020). Different areas were polluted by the point and the non-point source of PFASs and finally these substances entering into the groundwater system, even migrate to other places(Li et al. 2020a, Li et al. 2020b). PFASs distribution was hard to reveal the increasing and decreasing tendency due to tremendously industrial emission with different states of the waste air and water in park B. The industrial process and the production possess a strong effect on the distribution of PFASs far than geographical factor.
3.3 Health risk assessment
Environmental pollutants will enter into and finally accumulate in the human body by direct drinking(Gobelius et al. 2018). In our study, the groundwater as main drinking source is a serve exposure to local residents. The short chain PFASs were mainly detected in two parks whereas part of PFASs lacks out related evaluation messages, still all individual PFASs is hard to assess the health risk. Hence, combining with the mainly detected substances and available toxicity data, PFBA(C4), PFHpA(C7), PFHxA(C6), PFHxS(C6), PFOA(C8), PFOS(C8), PFNA(C9), PFUdA(C11) and PFDoA(C12) were assessed, respectively.
3.3.1 Individual PFASs (C≤8)
Along with the increasing of the aged groups, the RQ value show a decreasing tendency in park A and park B-WS and DS. This reveals that PFASs risk have a more remarkable influence on the infants and the children than the teenagers and the adults. RQ100%/20% value has a positive relationship with PFASs concentration in the water body and the negative relationship with DWEL. In park A, the highest DWEL value appears in PFBS (17,820,547.09 - 3,270,348.84 L·day-1), followed by PFHpA and PFHxA (654,069.01 – 3,564,109.42 L·day-1), respectively. But these three substances lack out determining the highest health risk whereas PFOA possesses a largest RQ value, followed by PFOS (Table S8). RQ100%/20% value of all assessed substances in park A is higher than park B and especially part of the RQ100%/20% value in park A has exceeded 1(Table S8). The RQ100%/20% value of PFBA, PFHpA, PFHxA and PFHxS in case 1-4 were in the reasonable range and this demonstrates that the max and the mean concentration of four substances is limited to exceed the RQ reference line of 0.2 and 1. Short chain PFASs (C≤7) in groundwater is hard to constitute the health risk to local residents. But PFOA and PFOS in part of sampling sites exceeding RQ=0.2, even RQ=1, has given rise to harm human body in case 1-4. The RQ100%/20% value of PFOA in case 1 and case 3 were in the range of 0.26-1.44 and 1.32-7.20 for aged groups, and the range of 0.23-0.64 for birth - 6years in case 4 (Fig. 4a). What counts is that total PFASs contamination degree is hard to pose health risk in case 2 for all aged groups. Causally, the water flow is the main impact to take away PFASs for the part of sampling sites with high concentration and this will dilute PFASs concentration, even decrease the toxic effect(Hongkachok et al. 2018, Ikkere et al. 2018, Jin et al. 2015). The ADI value difference between PFOS (0.15 μg/ kg day−1) and PFOA (1.5 μg/ kg day−1) has exceeded 1 orders of magnitude. However, the concentration of PFOS is 1-2 magnitude lower than PFOA which causes the RQ100%/20% value of PFOS is hard to exceed PFOA. The RQ100%/20% value of PFOS in case 1-2 and 4 was difficult to exceed RQ=0.2 expect for case 3, which is similar to PFOA. For different periods in park B, the degree of RQ100%/20% value in WS is higher in DS due to the high concentration was detected in WS (Fig 4b,c and Table S8). Based on all terms of different aged groups except the case 3, only PFOA in WS exceeds the control line of RQ=0.2 but is not exceed RQ=1. The range of RQ20% value for PFOA in case 3 was 0.22-0.62 that the max value of 0.62 for birth-3 months and the mix value of 0.22 for 3-6 years (Fig 4b,c). Many studies have revealed that PFASs can harm human body but different types of PFASs have various effect mechanisms(Graetz et al. 2020, He et al. 2018, Li et al. 2019, Miranda et al. 2021). Local residents regard the polluted water as drinking source, having larger exposure opportunity of touching these toxic substances.
3.3.2 Individual PFASs (C>8)
Long chain PFASs(C>8) were detected with lower concentration than short chain PFASs, which means the low concentration was also related to lower RQ value. Three individual-PFASs of PFNA(C9), PFUdA (C11) and PFDoA(C12) were assessed. The result reveal that three long chain PFASs is hard to pose potential risk to the local residents. The range of RQ value is that the negligible (ng)- 0.002, ng-0.0001 and ng-0.006 for PFNA, PFUdA and PFDoA, respectively (Table S8). These max RQ100%/20% value from three long chain PFASs were lower than 0.2 while the magnitude difference ranges from 1 to 3. Then in our study, extreme conditions are considered to perform operations to initially assess the risk to the human health. The group of birth-3 months can be easily affected by short chain PFASs in environmental media. If three assessed substances will cause health risk basing on the line of RQ=0.2, the concentration of 5.37, 6.53 and 0.13 ng·L-1 for PFNA, PFUdA and PFDoA were calculated, respectively. There was a huge difference of 1-2 magnitudes between the predicted and actual concentration. Thus, the long chain PFASs for assessed substances were difficult to threaten the human body.
3.3.3 Mixed risk of PFASs
The mixed risk assessment was increasingly concerned due to the different individual PFASs existing in the environmental media. But the related assessing method is limited to develop resulting from the complex mechanism in current environment system. The most serious situation appears in the park A of RQmix (4≤C≤12) =8.07 for birth-3 months (Table S6). ∑PFASs (4≤C≤7) and ∑PFASs (9≤C≤12) were hard to reach the line of RQ=0.2, and RQmix of these substances were limited to exceed that RQ=1 (Table S7). Expect for PFOS and PFOA, all PFASs (4≤C≤7 or 9≤C≤12) were analyzed with lower RQ contribution. Especially, the risk of the individual PFASs with low RQ value without exceeding RQ=0.2, then the RQmix value was hard to possess the risk to all aged groups. In this process, PFOS and PFOA risk values determine the crucial factor to the size of RQmix. In other words, if the risk from PFOS and PFOA can be effectively impaired, the RQmix of all the individual-PFASs will deeply decrease. Mixed RQmix value mainly relies on the PFOA and PFOS contribution rate (80%-90%). This indicates that PFOA and PFOS should be taken into consideration and listed as priority control. Notably, RQmix value from the group of teenagers (16-21 years) appears a small spike relative to two adjacent groups both 11-16 years and >21 years. In comparison with the infant and child, teenager group is the third group which is influenced by PFASs due to the weight and drinking intake.