Occurrence, Distribution, and Risk of Organophosphate Flame Retardants in Sediments from Jiulong River Estuary and Adjacent Western Taiwan Strait, China

Abstract

Organophosphate ester flame retardants (OPFRs) are widely prevalent in the environment and are of significant concern because of their potential toxicity to human health and wildlife. In this study, the concentration, frequency, spatial distribution, potential sources, and ecological risks of OPFRs in sediments from the Jiulong River estuary and the adjacent western Taiwan Strait were investigated. Concentrations of four of the five studied OPFRs were between <LOD and 36.6 ng/g. The distribution of all OPFRs, except 2-Ethylhexyl diphenyl phosphate (EHDPP), remained highly consistent with hydrological (salinity) trends. Furthermore, a significantly positive correlation between EHDPP and total concentrations suggested that it may be the dominant contaminant at both sites. Principal element analysis indicated multiple sources of OPFRs, which were categorized as emissions from road runoff and surface traffic, effects of atmospheric deposition and hydrologic conditions, and a combination of industrial and population effects. Ecological risk indicates that tris (chloroethyl) phosphate (TCEP) and triphosphate ester (2,3-dibromopropyl) (TDBPP) have almost no risk, tris (clorisopropyl) phosphate (TCPP) generally has low risk, while EHDPP has moderate risk with the highest value of 0.487 in the sediments from both sites. Meanwhile, TCPP and TCEP exhibit lower theoretical health risks but are still not negligible. Overall, this work provides data to support global pollutant studies and facilitate the implementation of pollutant control strategies.

1. Introduction

Organophosphate ester flame retardants (OPFRs) are widely used as alternatives to brominated flame retardants in various commercial products, including furniture, textiles, electronics, vehicles, and petroleum, and for plasticizing, defoaming, and improving fire safety [1,2]. From 1992 to 2015, the global consumption of OPFR increased from 100,000 t to 680,000 t [1,3]. In 2007, the total OPFR consumption was approximately 85,000 t in Europe and 70,000 t in China, with an annual growth rate of 15% [4,5]. Moreover, OPFRs were detected in different environmental matrices, such as indoor dust, air, water, soil, and sediments [6,7,8,9,10]. Their long-distance mobility caused detectable OPFR pollution in Canada’s Arctic region and the central Arctic Ocean region [11,12]. Additionally, OPFR residues were detected not only in marine and freshwater biota but also in humans [13,14,15,16]. Due to the increasing demand for OPFRs by industries, their ubiquitous distribution in the environment has raised concerns about their potential toxicity and risk to human health and wildlife [17].

Previous studies showed that exposure to OPFR leads to adverse effects on animals and humans, including acute, reproductive, and developmental toxicity, neurotoxicity, organ toxicity, genotoxicity, mutagenicity, and endocrine disruption [18]. For example, tris (1,3-dichloroisopropyl) phosphate (TDCPP) affects liver cells and neurons, while tris (clorisopropyl) phosphate (TCPP) affects fertility [19]. In addition, OPFRs (especially water-soluble OPFR compounds) have high acute toxicity and endocrine-disrupting effects on aquatic organisms, such as algae, daphnia, and fish [20,21]. The U.S. Environmental Protection Agency (EPA) estimated that 80% of human exposure to OPFRs comes from indoor dust [15]. A population epidemiological study found that human exposure to TCEP via indoor dust was positively associated with papillary thyroid cancer [odds ratios = 2.42 (95% CI: 1.10, 5.33), p = 0.03] [17]. In the past decade, almost all produced OPFRs were detected in marine and freshwater animals. In Canada, multiple OPFRs were detected in lake trout (Salvelinus namaycush) and walleye (Sander vitreus) collected in different lakes, with detected concentrations ranging from <0.07 to 9.8 ng/g (WW) [13]. OPFRs enter the marine environment through human activities and accumulate in the ocean via rivers, lakes, and atmospheric deposition [20,22,23]. Therefore, it is important to determine the sources of OFPRs and their mode of distribution, migration, and environmental transformation to protect marine aquatic ecosystems.

Most studies on OPFRs tend to focus on dust, air, and surface water. A limited number of studies have focused on sediments, which are an important component of the marine environment and are closely related to aquatic fauna and ecological cycles [12,17,24]. In the Asia-Pacific region, OPFRs in sediments were primarily reported in the Yellow Sea, Bohai Sea, Laizhou Bay, Pearl River Delta, and offshore waters of Xiamen [17,23,25]. The Jiulong River estuary and adjacent western part of the Taiwan Strait have developed industries, high urbanization rates, and significant concentrations of halogenated flame retardants [17]. Although OPFRs are gradually being phased out, their hazards cannot be ignored [17,26,27]. In this study, we sampled inland estuaries and land–sea junctions of the aforementioned regions to study the content, distribution, pollution characteristics, and migration of OPFRs in sediments and analyzed their sources and potential risks.

2. Materials and Methods

2.1. Chemicals and Reagents

OPFR standards, including tris (chloroethyl) phosphate (TCEP), triphosphate ester (2,3-dibromopropyl) (TDBPP), tris (1,3-dichloroisopropyl) phosphate (TDCPP), 2-Ethylhexyl diphenyl phosphate (EHDPP), tris (clorisopropyl) phosphate (TCPP), and d15-triphenyl phosphate (d15-TPP) were purchased from Cambridge Isotope Laboratories Inc. (Andover, MA, USA). Detailed structures and physicochemical properties of the OPFRs are listed in Table S1. Acetonitrile, methanol, and formic acid (HPLC grade) were purchased from Anpel Laboratory Technologies Inc. (Shanghai, China). Ultra-pure water was produced using a Milli-Q water purification system (Millipore, Bedford, MA, USA).

2.2. Sampling Collection

The study area included the Jiulong River Estuary (JRE) and the adjacent western part of the Taiwan Strait (WTPS), as shown in Figure 1. The surface sediment samples were collected using a Van Veen grab sampler (35 cm × 25 cm × 30 cm) at a depth of 0−5 cm. Samples from JRE (n = 16) were collected in January and February of 2019, and samples from the WTPS (n = 12) were collected in November and December of 2018. The sampling locations were identified accurately using a GPS locator. Furthermore, the collected sediment samples were stored in solvent-cleaned aluminum foil, sealed in polyethylene bags, placed on ice and transported to the laboratory, and frozen at −20 °C until further analysis. The physicochemical properties of the sediment overlying water were determined using YSI-EXO2 (Palo Alto, CA, USA), as shown in Table S2. The Jiulong River Basin is located in the southeast of Fujian Province, where the economy is relatively developed. The Jiulong River Basin has a total population of 5.43 million and a cultivated area of 1500 km². The downstream Zhangzhou Plain covers an area of 567 km². The Jiulong River flows through the dense industrial parks, agricultural areas, and aquaculture areas. Most of the production wastewater is discharged into the Jiulong River Basin, which causes great harm to the water quality of the Jiulong River, resulting in the increasingly serious pollution problem. As the receiving water body of the Jiulong River, the water quality of the Taiwan Strait is bound to be affected to some extent.

2.3. Sample Preparation and Instrumental Analysis

Samples were prepared following the method described by Li et al. [28]. Briefly, sediment samples were freeze-dried and sieved; approximately 0.5 g dry weight (dw) of each sample was spiked with an internal deuterated standard (d15-TPP). Then, an ultrasound-assisted extraction was performed thrice with 20 mL of an acetonitrile–water mixture (6:4). After centrifuging this mixture at 4000 rpm, ~40 mL of water was added to ~60 mL of the supernatant. Furthermore, the vial was filtered with a disposable 0.22 μm hydrophilic polytetrafluoroethylene (PTFE) syringe filter (Anpel Laboratory Technologies Inc. Shanghai, China). Solid-phase extraction (SPE) cartridges (Oasis HLB 200 mg, Waters, Milford, MA, USA) were conditioned with 4 mL of acetonitrile and 4 mL of Milli-Q water. After loading the sample at a rate of 0.5 mL/min, the cartridge was rinsed with 3 mL of Milli-Q water and vacuumed for 30 min until it was dry. Target compounds were eluted twice with 8 mL of acetonitrile (ACN). The combined eluate was passed under a gentle nitrogen flow to near dryness, reconstituted with 1 mL ACN, and filtered through a 0.22 μm PTFE syringe filter.

The analyses of target OFPRs were performed using an Agilent 1260 liquid chromatograph coupled to an Agilent 6460 triple quadrupole mass spectrometer (Palo Alto, CA, USA). The LC system was equipped with an Agilent EC-C18 (3.0 mm × 100 mm, 2.7 μm particle size) that was maintained at 35 °C. The mobile phases consisted of water (A) and ACN (B) at a flow rate of 0.5 mL/min. The following gradient was employed: 85% B ramped to 90% B in 2.50 min (linear), followed by a linear increase to 95% B in 1.75 min (held for 4.20 min), and then changed to 5% B for 4.08 min. Moreover, a 10 μL aliquot of the sample was injected into the LC system for analysis. The detailed LC-MS/MS analytical parameters for the determination of the target OPFRs are shown in Table S3.

2.4. Quality Assurance and Quality Control

This study applied a quality control procedure for data analysis [29]. D15-TPP was added to the samples as a surrogate to assess method recovery. Recovery of the surrogate was between 67% and 119%. The limit of detection (LOD) for each OPFR was defined as signal-to-noise (S/N) ratios of 3. Instrumental calibration standards were analyzed after every five samples to monitor instrument stability, and regression coefficients (R²) of the calibration curves were greater than 0.99.

2.5. Date Analysis

Using the sediment no-effect concentration (PNEC) as a water quality criterion to assess the ecological risk of OPFRs in the Jiulong River Estuary and Taiwan Strait, the predictions of five OPFRs were calculated using Equation (1) in combination with data from ECHA (European Chemicals Agency) [30].

(1)$RQ=\frac{MEC}{PNEC}$

where RQ is the ecological risk, and MEC is the actual concentration.

To assess the human exposure dose, OPFRs in sediments were converted to concentrations in the aqueous phase by multi-media imputation simulations, which were then converted to doses based on the USEPA exposure equation (USEPA, 2017). The equation is expressed as follows

(2)$CW=CS\times \rho$

(3)$EX{P}_D=\frac{CW\times IR}{BW}$

where CS is the concentration in sediment (ng/g), ρ is the conversion factor between the aqueous phase and sediment based on a provided, CW is the concentration in water (ng/L), EXP_D is the exposure dose (ng/kg/bw/day), IR is the uptake rate (L/d), and BW is the body weight (kg). The RQ values for health risks were determined using the following equation.

(4)$RQ=\frac{EX{P}_D}{RfD}$

where RfD is the USEPA oral reference dose.

Spearman correlation analysis was performed using IBM SPSS statistical software (IBM, Armonk, NY, USA) for the pollutant-related parameters at both sites. Descriptive statistics (i.e., minimum, maximum, mean, and standard deviation) were used to analyze the relationship between residual concentrations.

3. Results and Discussions

3.1. Concentration of OPFRs in Sediments

As shown in Table 1 and Table 2, the concentration ranges of all OFPRs ranged between <LOD and 36.63 ng/g; furthermore, their concentration level was 10² ng/g, which was high according to China’s standards [31]. Previous studies estimated that seven OPFRs in the sediments of Taihu Lake in China had concentration ranges of 3.4–14 ng/g, and five OPFRs in the southwest coast of Taiwan had concentration ranges of 1.0–13 ng/g; these values were lower than the values obtained in the present study [32,33]. Meanwhile, the concentration range in the Yangtze River was similar to the present study, ranging from 3.37 to 29.65 ng/g [34]. However, studies in the Pearl River Delta and in Spain found that OPFR concentrations could reach 10³ orders of magnitude, with maximum concentrations of 470 ng/g and 824 ng/g, respectively [17,33,35]. Of the five selected OPFRs, only TDCPP was not detected, probably due to local policies prohibiting the use of this compound.

Three of the four detected compounds, namely TCEP, TCPP, and EHDPP, had concentrations greater than 50%, indicating that they were more common in the sediments of the sampling sites. The detection rates of TCEP and TCPP in sediments obtained from the Yangtze River and Taihu Lake in China were close to 100% [7,34], and the detection rates of these two substances from the United States and Korea also reached 100%, 100% and 60%, 80% [36,37], respectively, indicating that their results were consistent with the present study. EHDPP in sediments is relatively less studied. Previous studies estimated that their detection rate in Korean sediments was 50%, while their detection rate in Guangxi, China, was only 16.7% [33,36]. Therefore, in this study, EHDPP may have a specific source of contamination. TDBPP studies tend to focus on water media because of the high detection rate of the substance in the water column, which was only 35.7% for sediments in this study [38].

The concentration distribution of the detected OFPRs was consistent with the detection rate. Detection rates of TCEP, TCPP, and EHDPP were one order of magnitude higher than those of TDBPP. The average concentration of EHDPP was 36.6 ng/g, which was significantly higher than its concentration in the sediment of Taihu Lake (ND−0.94 ng/g) [7]. EHDPP was reported to be more prominent in processed foods; thus, the surrounding population density and food industry could be important influencing factors [40]. Concentration ranges of TCEP and TCPP were similar and were consistent with previous reports on OPFRs [41,42].

3.2. Spatial Distribution of OPFRs in Sediments of the Jiulong Estuary and Western Taiwan Strait

As seen in Figure 2, TDBPP, TCEP, and TCPP showed similar distribution characteristics in both regions, with their concentrations ranging from 0 to 20 ng/g.

The most obvious difference was observed in EHDPP, which were two orders of magnitude higher in the Taiwan Strait region than in the Jiulong region. Electronic waste is an important source of pollutants, and the Beijiang River basin upstream of the Jiulong River traverses the largest electronics factory waste recycling area in Guangzhou [33,43]. As a result, pollutant concentrations tend to be high in the Jiulong River basin. In addition, the neighboring city of Xiamen, which is highly developed, densely populated, and the industrial area, also plays a role. High concentrations of certain OPFRs were found in P₂, P₃, and P₄. TCEP concentrations of 14.3 and 19.7 ng/g were found in P₂ and P₃, respectively, and TCPP exceeded 10¹ orders of magnitude in P₃ and P₄. TCPP and TCEP are often used as flame retardants in industries that produce automotive, rubber, polyurethane foam (PUF), and textile coatings [17,44]. Pollutants from these industries may be discharged into water bodies and migrate into sediments [28].

TDBPP and TCPP were dominant in A₁-A₃, TCEP and TDBPP in A₄-P₄, and EHDPP and TCEP in Q₁-X₃. TCEP is being gradually replaced by TCIPP in many countries due to its carcinogenic and viscous nature [28]. The high level of TCEP in the sampling sites indicates its widespread use and low removal efficiency of the nearby wastewater treatment.

Several studies showed differences in OPFR concentrations in estuarine and coastal areas [17,28,45]. The concentrations generally show a decreasing trend from west to east owing to a strong dilution effect. When the pollutants reach the river–sea interface, they become disturbed by the mixing of freshwater and seawater and are deposited with the sediments at the boundary. The data for TOC, TN, and ammonia nitrogen in Table S2 also illustrate this phenomenon well. As the geographical location changes from west to east, the three data show a decreasing trend, especially for TOC, which decreases more significantly by one order of magnitude. The decrease in the concentrations of TDBPP, TCPP, and TCEP from the estuary to the coast was not obvious, indicating that river runoff was not their only source. However, the concentration of EHDPP showed a sharp increase in X₁-X₃ (marine area), which was much higher than the other points in the upper Taiwan Strait. This may be due to the input of coastal currents from winter to spring in Fujian and Zhejiang [28]. According to the relevant literature, industrial wastes discharged from coastal areas are a potentially important source of OPFRs in the ocean [46,47]. In addition, atmospheric deposition caused by the East Asian monsoon has an important influence on the spatial distribution of organic pollutants [28].

The spatial distribution of OPFRs in the Jiulong area may also be influenced by hydrological conditions. Watersheds can be roughly divided into three zones based on salinity gradients and local conditions [43]. The estuarine zone (RR, salinity < 5‰) is located at the end of the river and is mainly influenced by river runoff. The river–sea zone (RMR, 5‰ ≤ salinity ≤ 25‰) is influenced by both rivers and oceans according to the tidal cycle of the Taiwan Strait. The marine zone (MR, salinity ≥ 25‰), located at the end of the ocean, is heavily influenced by seawater. Based on geographical location and watershed division, sites A₁-A₃ are located in the RR, sites A₄-A₁₄ and B₁-B₂ in the MR, and sites X₁-X₃ are in the RMR. From Figure 3b, it can be seen that the distribution of the three OPFRs, except EHDPP, is consistent with salinity; that is, the distribution of OPFRs showed a decreasing trend with increasing salinity, which was similar to the result of previous studies [28]. It was reported that the large-sized particles in the RR zone were remobilized by deposition, while the small-sized particles were transported to the RMR and MR zones [28]. The small particles allowed OPFRs to occur mainly in water and suspensions, thus decreasing their concentrations in the sediments [7]. This indicated that the distribution of OPFRs in the sediments of the Jiulong River and Taiwan Strait was not only controlled by river runoff but also influenced by hydrodynamic, oceanic, and anthropogenic factors, which was consistent with the results of other studies [17,25,28].

3.3. Source Analysis of OPFRs

The correlations between the four OPFRs and total concentrations were analyzed using SPSS statistical software.

Table 3 show the correlation coefficients and significance levels. It was observed that the correlations among the OPFRs were weak. However, there was a good correlation between TCEP and TCPP, which indicated significant variability in their origin. This conclusion was consistent with the discussion in Section 3.2. Except for TDBPP, all three OPFRs were positively correlated with the total mean concentration. TCEP and TCPP were moderately correlated (0.3–0.5), whereas EHDPP had the most pronounced correlation of 0.769. Usually, a correlation coefficient between 0.5 and 0.8 is considered strong. Thus, EHDPP, TCEP, and TCPP were the highest contributors of OPFRs in the Jiulong River Estuary and Taiwan Strait.

To further analyze the possible contaminant sources, principal component analysis (PCA) was performed on the parameters associated with the OPFRs. As seen in Figure 4, the three main factors, PC₁, PC₂, and PC₃, had values of 34.0, 29.1, and 18.8%, respectively, accounting for 81.9% of the total variance. TCPP and TCEP had high loadings (>0.6) in PC₁. Furthermore, they were reported to be highly enriched in the drinking water and surface waters of China, and their sources were closely related to emissions from vehicles and marine traffic [46,48,49]. Therefore, PC₁ can be considered to be a combination of emissions from road runoff, vehicles, and marine traffic.

PC₂ accounted for 29.1% of the total variance, with TDBPP (0.73907) having the highest value. TDBPP is often used as a flame retardant in cellulose, triacetate, and polyester fabrics [28]. It is commonly found in fabrics [17]. In addition, TDBPP solid waste from textile processing plants is an important source of pollution. Therefore, PC₂ can be defined as a combination of emissions from industry and population.

The loads of EHDPP and TCPP in PC₃ were more pronounced. Given that they are commonly present in the air, atmospheric particles, and dust [34], PC₃ was a combination of the effect of atmospheric deposition and hydrology.

3.4. Risks of OPFRs

In this study, ecological risks were assessed based on risk quotients (RQs) calculated by the PNEC. The PNEC values provided by ECHA for TCPP, TCEP, TDBPP, and EHDPP are 292 ng/g, 200 ng/g, 92.9 ng/g, and 81.4 ng/g, respectively [50]. The risk assessment criteria were as follows: low risk (0.01 < RQ < 0.1), medium risk (0.1 ≤ RQ < 1), and high risk (RQ ≥ 1).

Figure 5a show that all detected OPFRs have low or no risk at the sampling points. TCEP is almost risk-free at most sampling sites and is a low risk only in zone Q₄, with a value of 0.0199. TDBPP is similar, reaching 0.01 only in zone A₂. TCPP has a more balanced risk distribution, ranging from 0 to 0.0565, and generally has a low risk. In contrast, the ecological risk of EHDPP was more severe, with three areas falling into the medium risk zone. The highest risk values for EHDPP were found in the X zone, with 0.457, 0.180, and 0.487, respectively, which were far from the risk values of the other three substances. This was consistent with the results of the study on Taihu Lake [7]. Moreover, the risk of pollutant coexistence in the basin should not be underestimated. The high risk of certain substances increases the risk value of the total concentration. Thus, there is a need for long-term detection and control when the overall risk level in the two zones is high.

As shown in Table S4, the high water solubility of OPFRs in this study suggests that most of these compounds tend to be distributed in the aqueous phase. Therefore, the higher ecological risk of TCEP, TCPP, and EHDPP in sediments suggests that these compounds are also at high risk in water. In addition, the studied watersheds are also used for drinking water and fisheries, and OPFRs in sediments can enter humans indirectly through drinking water and organisms, which has important implications for human exposure and human health and therefore requires further assessment of human health risks [30]. At the same time, the sensitivity analysis of the organism is an important parameter for subsequent studies when conducting risk assessments.

The theoretical health risks of OPFR in adult males and females in this study were analyzed using the relevant formulae (Figure 5b). Since there were no other RfD reference values for OPFR, only the health risks for TCPP and TCEP were calculated. The RfD values for TCPP and TCEP were 10ng/g/day and 7ng/g/day, respectively, with conversion factors ρ of 0.0412 and 0.544 [51,52]. The remaining parameters are shown in Table S5. As can be seen in Figure 5b, most of the sampling sites had almost no risk. Only the TCPP located in parts A₇-P₄ had low risk with the highest value of 0.0439. In addition, the risk values were generally higher for males relative to females. In conclusion, the OPFR in both the Jiulong River Basin and the Taiwan Strait have low health risks but are still not negligible.

4. Conclusions

In this study, OPFRs were detected and analyzed in the sediments obtained from the Jiulong River Basin and the Taiwan Strait; their concentrations, frequencies, spatial distributions, and ecological risk profiles were systematically described. Four of the five OPFRs were detected in the samples. Furthermore, TCPP, TCEP, and EHDPP were the major contaminants, given their high frequency of use. The spatial distribution of OPFRs was mainly influenced by river runoff, monsoon, and hydrology. Correlation and PCA showed that the sources of OPFRs were multifaceted, as follows: (1) emissions from road runoff and surface traffic, (2) influence of atmospheric deposition and hydrological conditions, and (3) a combined effect of industry and population. The ecological and health risks indicate that the overall risk values of the OPFR are not negligible, despite the low individual risk values. This study can facilitate the development and implementation of future pollution control strategies for the relevant sectors in two selected regions and provide data support for global pollutant studies.

Author Contributions

Conceptualization, L.C. and C.P.; methodology, L.C.; software, F.Z.; validation, L.C. and S.W.; formal analysis, C.P.; investigation, L.C.; resources, J.M.; data curation, Y.S. and S.W.; writing—original draft preparation, F.Z. and C.P.; writing—review and editing, F.Z., C.P. and J.M.; visualization, J.D.; supervision, M.Y.; project administration, J.M.; funding acquisition, M.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was sponsored by the National Natural Science Foundation of China (Grant No.31971523) and the Natural Science Foundation of Fujian province, China (Grant No.2018J05066).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated and analyzed during this study are included in this published article and its supplementary information files.

Conflicts of Interest

The authors declare the nonexistence of any conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Tables

Figure 1. Sampling sites in the estuary of the Jiulong River and the western part of the Taiwan Strait.

Figure 2. Distribution of OPFRs concentrations at sampling sites at (a) the Jiulong River estuary and (b) the western Taiwan Strait.

Figure 3. (a) Percentage distribution of OPFRs in the Jiulong River estuary and Taiwan Strait sampling sites; (b) Variation in OPFRs concentrations in different regions.

Figure 4. PCA master element analysis diagram. (JL: Jiulong River Estuary, TW: Taiwan Strait).

Figure 5. (a) Ecological risk assessment of all sampling sites; (b) Health risk assessment of all sampling sites.

Table 1

Concentrations and frequencies of the five OPFRs in the sediments of Jiulong River Estuary and adjacent Western Taiwan Strait.

OPFRs	Min	Max	Range	Med	Average	SD	DF%
TDBPP	0.00	0.926	<LOD-0.926	0.0971	0.289	0.343	35.7
TCEP	0.00	19.7	<LOD-19.7	0.771	3.34	4.92	75.0
TCPP	0.00	16.5	<LOD-16.5	1.61	4.48	5.55	60.7
EHDPP	0.00	36.6	<LOD-36.6	0.448	5.57	11.7	67.9
TDCPP	NA	NA	<LOD	NA	NA	NA	0.00
∑OPFRs	0.118	39.9	-	4.61	9.12	11.2	-

Table 2

Global concentrations of OPFRs (ng/g).

Region	TDBPP	TCEP	TCPP	EHDPP	Year	Reference
Yangtze River	-	3.13–4.08	3.37–29.7	-	2018	[34]
Taihu Lake	ND	1–3.17	0–2.19	ND-0.94	2018	[7]
Qinzhou Bay	-	0–3.076	-	0–94.4	2021	[39]
Pearl River Delta	-	ND-58	0.91–185	-	2016	[35]
United States	-	0.168–5.6	0.146–36.8	-	2016	[37]
Australia	-	0–160	33–170	-	2021	[17]
Korea	-	0–60	0–216	0–50	2018	[36]

Table 3

Correlation analysis between the components as well as total concentration.

	TDBPP	TCEP	TCPP	EHDPP	∑OPFRs
TDBPP	1.000	−0.177	−0.054	−0.134	−0.193
TCEP	−0.177	1.000	0.218	−0.100	0.405 *
TCPP	−0.054	0.218	1.000	−0.194	0.348
EHDPP	−0.134	−0.100	−0.194	1.000	0.769 **
∑OPFRs	−0.193	0.405 *	0.348	0.769 **	1.000

Note. “**”, Significant correlation at p  <  0.01; “*”, significant correlation at p  <  0.05.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijerph19042449/s1, Table S1: Structures and physicochemical properties of the target organophosphate flame retardants, Table S2: Sediment sampling information, Table S3: LC-MS/MS analytical parameters for determination of the target organophosphate flame retardants, Table S4: Acute and chronic HC5 of five OPFRs, Table S5: The parameters for health risk assessment.

References

1. Wang, Q.; Lam, J.C.-W.; Man, Y.-C.; Lai, N.L.-S.; Kwok, K.Y.; Guo, Y.Y.; Lam, P.K.-S.; Zhou, B. Bioconcentration, metabolism and neurotoxicity of the organophorous flame retardant 1,3-dichloro 2-propyl phosphate (TDCPP) to zebrafish. Aquat. Toxicol.; 2015; 158, pp. 108-115. [DOI: https://dx.doi.org/10.1016/j.aquatox.2014.11.001] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25461749]

2. Marklund, A.; Andersson, B.; Haglund, P. Screening of organophosphorus compounds and their distribution in various indoor environments. Chemosphere; 2003; 53, pp. 1137-1146. [DOI: https://dx.doi.org/10.1016/S0045-6535(03)00666-0]

3. Medical Advisory, S. Energy delivery systems for treatment of benign prostatic hyperplasia: An evidence-based analysis. Ont. Health Technol. Assess. Ser.; 2006; 6, pp. 1-121.

4. Ou, Y. Developments of organic phosphorus flame retardant industry in China. Chem. Ind. Eng. Prog.; 2011; 30, pp. 210-215.

5. Wang, X.; Liu, J.; Yin, Y. The Pollution Status and Research Progress on Organophosphate Ester Flame Retardants. Prog. Chem.; 2010; 22, pp. 1983-1992.

6. Chen, M.; Gan, Z.; Qu, B.; Chen, S.; Dai, Y.; Bao, X. Temporal and seasonal variation and ecological risk evaluation of flame retardants in seawater and sediments from Bohai Bay near Tianjin, China during 2014 to 2017. Mar. Pollut. Bull.; 2019; 146, pp. 874-883. [DOI: https://dx.doi.org/10.1016/j.marpolbul.2019.07.049]

7. Chen, M.; Liu, Y.; Guo, R.; Xu, H.; Song, N.; Han, Z.; Chen, N.; Zhang, S.; Chen, J. Spatiotemporal distribution and risk assessment of organophosphate esters in sediment from Taihu Lake, China. Environ. Sci. Pollut. Res.; 2018; 25, pp. 13787-13795. [DOI: https://dx.doi.org/10.1007/s11356-018-1434-3]

8. Cristale, J.; Ramos, D.D.; Dantas, R.F.; Machulek Junior, A.; Lacorte, S.; Sans, C.; Esplugas, S. Can activated sludge treatments and advanced oxidation processes remove organophosphorus flame retardants?. Environ. Res.; 2016; 144, pp. 11-18. [DOI: https://dx.doi.org/10.1016/j.envres.2015.10.008]

9. Xu, L.; Hu, Q.; Liu, J.; Liu, S.; Liu, C.; Deng, Q.; Zeng, X.; Yu, Z. Occurrence of organophosphate esters and their diesters degradation products in industrial wastewater treatment plants in China: Implication for the usage and potential degradation during production processing. Environ. Pollut.; 2019; 250, pp. 559-566. [DOI: https://dx.doi.org/10.1016/j.envpol.2019.04.058]

10. Zhu, F.; Pan, J.; Zou, Q.; Wu, M.; Wang, H.; Xu, G. Electron beam irradiation of typical sulfonamide antibiotics in the aquatic environment: Kinetics, removal mechanisms, degradation products and toxicity assessment. Chemosphere; 2021; 274, [DOI: https://dx.doi.org/10.1016/j.chemosphere.2021.129713]

11. Sun, Y.; Gong, X.; Lin, W.; Liu, Y.; Wang, Y.; Wu, M.; Kannan, K.; Ma, J. Metabolites of organophosphate ester flame retardants in urine from Shanghai, China. Environ. Res.; 2018; 164, pp. 507-515. [DOI: https://dx.doi.org/10.1016/j.envres.2018.03.031] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29604578]

12. Wu, Y.; Wang, X.; Ya, M.; Li, Y.; Hong, H. Distributions of organochlorine compounds in sediments from Jiulong River Estuary and adjacent Western Taiwan Strait: Implications of transport, sources and inventories. Environ. Pollut.; 2016; 219, pp. 519-527. [DOI: https://dx.doi.org/10.1016/j.envpol.2016.05.081] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27307268]

13. Hou, L.; Jiang, J.; Gan, Z.; Dai, Y.; Yang, P.; Yan, Y.; Ding, S.; Su, S.; Bao, X. Spatial Distribution of Organophosphorus and Brominated Flame Retardants in Surface Water, Sediment, Groundwater, and Wild Fish in Chengdu, China. Arch. Environ. Contam. Toxicol.; 2019; 77, pp. 279-290. [DOI: https://dx.doi.org/10.1007/s00244-019-00624-x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30963189]

14. Krzeminski, P.; Schwermer, C.; Wennberg, A.; Langford, K.; Vogelsang, C. Occurrence of UV filters, fragrances and organophosphate flame retardants in municipal WWTP effluents and their removal during membrane post-treatment. J. Hazard. Mater.; 2017; 323, pp. 166-176. [DOI: https://dx.doi.org/10.1016/j.jhazmat.2016.08.001] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27566858]

15. Zhang, Z.; Shao, H.; Wu, M.; Zhang, J.; Li, D.; Li, J.; Wang, H.; Shi, W.; Xu, G. Occurrence, Distribution, and Potential Sources of Organophosphate Esters in Urban and Rural Surface Water in Shanghai, China. Arch. Environ. Contam. Toxicol.; 2019; 77, pp. 115-126. [DOI: https://dx.doi.org/10.1007/s00244-019-00633-w]

16. Pan, C.; Zhu, F.; Wu, M.; Jiang, L.; Zhao, X.; Yang, M. Degradation and toxicity of the antidepressant fluoxetine in an aqueous system by UV irradiation. Chemosphere; 2022; 287, 132434. [DOI: https://dx.doi.org/10.1016/j.chemosphere.2021.132434]

17. Chen, M.-H.; Ma, W.-L. A review on the occurrence of organophosphate flame retardants in the aquatic environment in China and implications for risk assessment. Sci. Total Environ.; 2021; 783, 147064. [DOI: https://dx.doi.org/10.1016/j.scitotenv.2021.147064]

18. Bai, X.-y.; Lu, S.-y.; Xie, L.; Zhang, B.; Song, S.-m.; He, Y.; Ouyang, J.-p.; Zhang, T. A pilot study of metabolites of organophosphorus flame retardants in paired maternal urine and amniotic fluid samples: Potential exposure risks of tributyl phosphate to pregnant women. Environ. Sci.-Processes Impacts; 2019; 21, pp. 124-132. [DOI: https://dx.doi.org/10.1039/C8EM00389K]

19. Feng, L.; Ouyang, F.; Liu, L.; Wang, X.; Wang, X.; Li, Y.-J.; Murtha, A.; Shen, H.; Zhang, J.; Zhang, J.J. Levels of Urinary Metabolites of Organophosphate Flame Retardants, TDCIPP, and TPHP, in Pregnant Women in Shanghai. J. Environ. Public Health; 2016; 2016, 9416054. [DOI: https://dx.doi.org/10.1155/2016/9416054]

20. Dasguta, S.; Vliet, S.M.F.; Cheng, V.; Mitchell, C.A.; Kirkwood, J.; Vollaro, A.; Hur, M.; Mehdizadeh, C.; Volz, D.C. Complex Interplay Among Nuclear Receptor Ligands, Cytosine Methylation, and the Metabolome in Driving Tris(1,3-dichloro-2-propyl)phosphate-Induced Epiboly Defects in Zebrafish. Environ. Sci. Technol.; 2019; 53, pp. 10497-10505. [DOI: https://dx.doi.org/10.1021/acs.est.9b04127]

21. Giraudo, M.; Douville, M.; Houde, M. Chronic toxicity evaluation of the flame retardant tris (2-butoxyethyl) phosphate (TBOEP) using Daphnia magna transcriptomic response. Chemosphere; 2015; 132, pp. 159-165. [DOI: https://dx.doi.org/10.1016/j.chemosphere.2015.03.028] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25855008]

22. McGee, S.P.; Cooper, E.M.; Stapleton, H.M.; Volz, D.C. Early Zebrafish Embryogenesis Is Susceptible to Developmental TDCPP Exposure. Environ. Health Perspect.; 2012; 120, pp. 1585-1591. [DOI: https://dx.doi.org/10.1289/ehp.1205316] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23017583]

23. Hwang, J.H.; Van, S.P.; Choi, B.-J.; Chang, Y.S.; Kim, Y.H. The physical processes in the Yellow Sea. Ocean Coast. Manag.; 2014; 102, pp. 449-457. [DOI: https://dx.doi.org/10.1016/j.ocecoaman.2014.03.026]

24. Pantelaki, I.; Voutsa, D. Organophosphate flame retardants (OPFRs): A review on analytical methods and occurrence in wastewater and aquatic environment. Sci. Total Environ.; 2019; 649, pp. 247-263. [DOI: https://dx.doi.org/10.1016/j.scitotenv.2018.08.286]

25. Vasseghian, Y.; Alimohamadi, M.; Khataee, A.; Dragoi, E.-N. A global systematic review on the concentration of organophosphate esters in water resources: Meta-analysis, and probabilistic risk assessment. Sci. Total Environ.; 2022; 807, 150876. [DOI: https://dx.doi.org/10.1016/j.scitotenv.2021.150876]

26. Lai, N.L.S.; Kwok, K.Y.; Wang, X.-h.; Yamashita, N.; Liu, G.; Leung, K.M.Y.; Lam, P.K.S.; Lam, J.C.W. Assessment of organophosphorus flame retardants and plasticizers in aquatic environments of China (Pearl River Delta, South China Sea, Yellow River Estuary) and Japan (Tokyo Bay). J. Hazard. Mater.; 2019; 371, pp. 288-294. [DOI: https://dx.doi.org/10.1016/j.jhazmat.2019.03.029]

27. Qi, C.; Yu, G.; Zhong, M.; Peng, G.; Huang, J.; Wang, B. Organophosphate flame retardants in leachates from six municipal landfills across China. Chemosphere; 2019; 218, pp. 836-844. [DOI: https://dx.doi.org/10.1016/j.chemosphere.2018.11.150]

28. Li, R.; Liang, J.; Duan, H.; Gong, Z. Spatial distribution and seasonal variation of phthalate esters in the Jiulong River estuary, Southeast China. Mar. Pollut. Bull.; 2017; 122, pp. 38-46. [DOI: https://dx.doi.org/10.1016/j.marpolbul.2017.05.062]

29. Zhu, F.; Wang, S.; Liu, Y.; Wu, M.; Wang, H.; Xu, G. Antibiotics in the surface water of Shanghai, China: Screening, distribution, and indicator selecting. Environ. Sci. Pollut. Res.; 2021; 28, pp. 9836-9848. [DOI: https://dx.doi.org/10.1007/s11356-020-10967-x]

30. Sun, R.; Wu, M.; Tang, L.; Li, J.; Qian, Z.; Han, T.; Xu, G. Perfluorinated compounds in surface waters of Shanghai, China: Source analysis and risk assessment. Ecotoxicol. Environ. Saf.; 2018; 149, pp. 88-95. [DOI: https://dx.doi.org/10.1016/j.ecoenv.2017.11.012]

31. Ren, G.; Chu, X.; Zhang, J.; Zheng, K.; Zhou, X.; Zeng, X.; Yu, Z. Organophosphate esters in the water, sediments, surface soils, and tree bark surrounding a manufacturing plant in north China. Environ. Pollut.; 2019; 246, pp. 374-380. [DOI: https://dx.doi.org/10.1016/j.envpol.2018.12.020] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30577005]

32. Hu, Y.-X.; Sun, Y.-X.; Li, X.; Xu, W.-H.; Zhang, Y.; Luo, X.-J.; Dai, S.-H.; Xu, X.-R.; Mai, B.-X. Organophosphorus flame retardants in. mangrove sediments from the Pearl River Estuary, South China. Chemosphere; 2017; 181, pp. 433-439. [DOI: https://dx.doi.org/10.1016/j.chemosphere.2017.04.117] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28458218]

33. Hu, M.; Li, J.; Zhang, B.; Cui, Q.; Wei, S.; Yu, H. Regional distribution of halogenated organophosphate flame retardants in seawater samples from three coastal cities in China. Mar. Pollut. Bull.; 2014; 86, pp. 569-574. [DOI: https://dx.doi.org/10.1016/j.marpolbul.2014.06.009] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24974164]

34. Zha, D.; Li, Y.; Yang, C.; Yao, C. Assessment of organophosphate flame retardants in surface water and sediment from a freshwater environment (Yangtze River, China). Environ. Monit. Assess.; 2018; 190, 222. [DOI: https://dx.doi.org/10.1007/s10661-018-6587-5]

35. Tan, X.-X.; Luo, X.-J.; Zheng, X.-B.; Li, Z.-R.; Sun, R.-X.; Mai, B.-X. Distribution of organophosphorus flame retardants in sediments from the Pearl River Delta in South China. Sci. Total Environ.; 2016; 544, pp. 77-84. [DOI: https://dx.doi.org/10.1016/j.scitotenv.2015.11.089]

36. Lee, S.; Cho, H.-J.; Choi, W.; Moon, H.-B. Organophosphate flame retardants (OPFRs) in water and sediment: Occurrence, distribution, and hotspots of contamination of Lake Shihwa, Korea. Mar. Pollut. Bull.; 2018; 130, pp. 105-112. [DOI: https://dx.doi.org/10.1016/j.marpolbul.2018.03.009]

37. Maruya, K.A.; Dodder, N.G.; Sengupta, A.; Smith, D.J.; Lyons, J.M.; Heil, A.T.; Drewes, J.E. Multimedia screening of contaminants of emerging concern (CECS) in coastal urban watersheds in southern California (USA). Environ. Toxicol. Chem.; 2016; 35, pp. 1986-1994. [DOI: https://dx.doi.org/10.1002/etc.3348]

38. Hu, L.; Tao, Y.; Luo, D.; Feng, J.; Wang, L.; Yu, M.; Li, Y.; Covaci, A.; Mei, S. Simultaneous biomonitoring of 15 organophosphate flame retardants metabolites in urine samples by solvent induced phase transition extraction coupled with ultra-performance liquid chromatography-tandem mass spectrometry. Chemosphere; 2019; 233, pp. 724-732. [DOI: https://dx.doi.org/10.1016/j.chemosphere.2019.05.242]

39. Zhang, L.; Lu, L.; Zhu, W.; Yang, B.; Lu, D.; Dan, S.F.; Zhang, S. Organophosphorus flame retardants (OPFRs) in the seawater and sediments of the Qinzhou Bay, Northern Beibu Gulf: Occurrence, distribution, and ecological risks. Mar. Pollut. Bull.; 2021; 168, 112368. [DOI: https://dx.doi.org/10.1016/j.marpolbul.2021.112368]

40. Luo, Q.; Shan, Y.; Muhammad, A.; Wang, S.; Sun, L.; Wang, H. Levels, distribution, and sources of organophosphate flame retardants and plasticizers in urban soils of Shenyang, China. Environ. Sci. Pollut. Res.; 2018; 25, pp. 31752-31761. [DOI: https://dx.doi.org/10.1007/s11356-018-3156-y]

41. Cao, D.; Guo, J.; Wang, Y.; Li, Z.; Liang, K.; Corcoran, M.B.; Hosseini, S.; Bonina, S.M.C.; Rockne, K.J.; Sturchio, N.C. et al. Organophosphate Esters in Sediment of the Great Lakes. Environ. Sci. Technol.; 2017; 51, pp. 1441-1449. [DOI: https://dx.doi.org/10.1021/acs.est.6b05484] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28084727]

42. Giulivo, M.; Capri, E.; Kalogianni, E.; Milacic, R.; Majone, B.; Ferrari, F.; Eljarrat, E.; Barcelo, D. Occurrence of halogenated and organophosphate flame retardants in sediment and fish samples from three European river basins. Sci. Total Environ.; 2017; 586, pp. 782-791. [DOI: https://dx.doi.org/10.1016/j.scitotenv.2017.02.056] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28215802]

43. Lin, B.H.; Lei, H.Y.; Guan, B.C.; Wang, M.G.; Liu, H.R. Elemental Characteristics and Geochemical Significance of Surface Sediment in Jiulong River Estuary. J. Xiamen Univ. (Nat. Sci.); 2009; 48, pp. 450-455.

44. Han, Q.F.; Zhao, S.; Zhang, X.R.; Wang, X.L.; Song, C.; Wang, S.G. Distribution, combined pollution and risk assessment of antibiotics in typical marine aquaculture farms surrounding the Yellow Sea, North China. Environ. Int.; 2020; 138, [DOI: https://dx.doi.org/10.1016/j.envint.2020.105551] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32155507]

45. Sun, J.; Huang, J.; Zhang, A.; Liu, W.; Cheng, W. Occurrence of phthalate esters in sediments in Qiantang River, China and inference with urbanization and river flow regime. J. Hazard. Mater.; 2013; 248, pp. 142-149. [DOI: https://dx.doi.org/10.1016/j.jhazmat.2012.12.057]

46. Kim, U.-J.; Oh, J.K.; Kannan, K. Occurrence, Removal, and Environmental Emission of Organophosphate Flame Retardants/Plasticizers in a Wastewater Treatment Plant in New York State. Environ. Sci. Technol.; 2017; 51, pp. 7872-7880. [DOI: https://dx.doi.org/10.1021/acs.est.7b02035]

47. Li, H.; La Guardia, M.J.; Liu, H.; Hale, R.C.; Mainor, T.M.; Harvey, E.; Sheng, G.; Fu, J.; Peng, P.a. Brominated and organophosphate flame retardants along a sediment transect encompassing the Guiyu, China e-waste recycling zone. Sci. Total Environ.; 2019; 646, pp. 58-67. [DOI: https://dx.doi.org/10.1016/j.scitotenv.2018.07.276]

48. Pang, L.; Yang, H.; Wang, Y.; Luo, X.; Liu, S.; Xiao, J. Organophosphate flame retardants in total suspended particulates from an urban area of zhengzhou, China: Temporal variations, potential affecting factors, and health risk assessment. Ecotoxicol. Environ. Saf.; 2019; 176, pp. 204-210. [DOI: https://dx.doi.org/10.1016/j.ecoenv.2019.03.092]

49. Ma, J.; Zhu, H.; Kannan, K. Organophosphorus Flame Retardants and Plasticizers in Breast Milk from the United States. Environ. Sci. Technol. Lett.; 2019; 6, pp. 525-531. [DOI: https://dx.doi.org/10.1021/acs.estlett.9b00394]

50. Li, J.; He, J.; Li, Y.; Liu, Y.; Li, W.; Wu, N.; Zhang, L.; Zhang, Y.; Niu, Z. Assessing the threats of organophosphate esters (flame retardants and plasticizers) to drinking water safety based on USEPA oral reference dose (RfD) and oral cancer slope factor (SFO). Water Res.; 2019; 154, pp. 84-93. [DOI: https://dx.doi.org/10.1016/j.watres.2019.01.035]

51. Yin, H.; Li, L.; Li, S.; Deng, X.; Liu, Q. Multimedia fate simulation of typical OPEs in Chengdu. Acta Sci. Circumstantiae; 2017; 37, pp. 3489-3497.

52. Mcgoldrick, D.J.; Letcher, R.J.; Barresi, E.; Keir, M.J.; Small, J.; Clark, M.G.; Sverko, E.; Backus, S.M. Organophosphate flame retardants and organosiloxanes in predatory freshwater fish from locations across Canada. Environ. Pollut.; 2014; 193, pp. 254-261. [DOI: https://dx.doi.org/10.1016/j.envpol.2014.06.024] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25063913]