Ecological risk assessment for organophosphate esters in the surface water from the Bohai Sea of China using multimodal species sensitivity distributions

Organophosphate esters (OPEs) as the foremost substitutes of brominated flame retardants have been ubiquitously found in the aquatic environment around the world. However, the information on the community-level risks induced by OPEs to the marine ecosystem remains scarce. This study adopted ten commonly used species sensitivity distribution (SSD) parametric statistical approaches coupled with the acute-to-chronic transformation for the toxicity data to fit the sensitivity distributions of different species to four major OPE congeners including triethyl phosphate (TEP), tri-n-butyl phosphate (TnBP), tri(2-chloroethyl) phosphate (TCEP), and tris(1-chloro-2-propyl) phosphate (TCPP) in the surface water of the Bohai Sea. All SSD models except Exponential for TnBP, TCEP, and TCPP fitted well the chronic toxicity data for the four OPE congeners. Discrepancies appeared among the best fitting models for different congeners, which also happened to the fitting results from the multiple SSD models for each congener. Based on the best fitting models, the hazard concentrations corresponding to the cumulative probability of 5% were 3.58 mg/L, 0.116 mg/L, 1.30 mg/L, and 1.44 mg/L for TEP, TnBP, TCEP, and TCPP, respectively. The risks induced by the four OPE congeners to the Bohai Sea ecosystem were negligible during the monitoring period because of both the risk quotients and the hazard indexes far <0.1. This study drew a clear picture of the joint ecological risks of TEP, TnBP, TCEP, and TCPP to the Bohai Sea environment. The application of multimodal SSD statistical methods will benefit the accurate derivation of water quality criteria and the community-level ecological risk assessment for pollutants.