Pharmaceutical pollution of the world's rivers.

Environmental exposure to active pharmaceutical ingredients (APIs) can have negative effects on the health of ecosystems and humans. While numerous studies have monitored APIs in rivers, these employ different analytical methods, measure different APIs, and have ignored many of the countries of the world. This makes it difficult to quantify the scale of the problem from a global perspective. Furthermore, comparison of the existing data, generated for different studies/regions/continents, is challenging due to the vast differences between the analytical methodologies employed. Here, we present a global-scale study of API pollution in 258 of the world's rivers, representing the environmental influence of 471.4 million people across 137 geographic regions. Samples were obtained from 1,052 locations in 104 countries (representing all continents and 36 countries not previously studied for API contamination) and analyzed for 61 APIs. Highest cumulative API concentrations were observed in sub-Saharan Africa, south Asia, and South America. The most contaminated sites were in low- to middle-income countries and were associated with areas with poor wastewater and waste management infrastructure and pharmaceutical manufacturing. The most frequently detected APIs were carbamazepine, metformin, and caffeine (a compound also arising from lifestyle use), which were detected at over half of the sites monitored. Concentrations of at least one API at 25.7% of the sampling sites were greater than concentrations considered safe for aquatic organisms, or which are of concern in terms of selection for antimicrobial resistance. Therefore, pharmaceutical pollution poses a global threat to environmental and human health, as well as to delivery of the United Nations Sustainable Development Goals.

In Situ Insight into the Availability and Desorption Kinetics of Per- and Polyfluoroalkyl Substances in Soils with Diffusive Gradients in Thin Films.

The physicochemical exchange dynamics between the solid and solution phases of per- and polyfluoroalkyl substances (PFAS) in soils needs to be better understood. This study employed an in situ tool, diffusive gradients in thin films (DGT), to understand the distribution and exchange kinetics of five typical PFAS in four soils. Results show a nonlinear relationship between the PFAS masses in DGT and time, implying that PFAS were partially supplied by the solid phase in all of the soils. A dynamic model DGT-induced fluxes in soils/sediments (DIFS) was used to interpret the results and derive the distribution coefficients for the labile fraction (Kdl), response time (tc), and adsorption/desorption rates (k1 and k-1). The larger labile pool size (indicated by Kdl) for the longer chain PFAS implies their higher potential availability. The shorter chain PFAS tend to have a larger tc and relatively smaller k-1, implying that the release of these PFAS in soils might be kinetically limited but not for more hydrophobic compounds, such as perfluorooctanesulfonic acid (PFOS), although soil properties might play an important role. Kdl ultimately controls the PFAS availability in soils, while the PFAS release from soils might be kinetically constrained (which may also hold for biota uptake), particularly for more hydrophilic PFAS.

In Situ Understanding of the Effect of Manure on the Availability of Sulfonamide Antibiotics in Soils Using DGT.

In this study, the effects of manure on the availability of sulfonamide antibiotics (SAs) in soils were explored in situ by the Diffusive gradients in thin films (DGT) technique. Five antibiotics, including sulfadiazine (SDZ), sulfamethoxazole (SMX), sulfamethazine (SMZ), sulfachloropyridazine (SCP), and sulfadimethoxine (SDM), were selected as target compounds. Results showed that the manure application to soil could reduce the antibiotic availability indicated by DGT. DGT measurement (CDGT) showed good correlations with the soil solution concentrations (Cd). Manure application can suppress the fluxes of SAs from the soil to the soil solution. Using the DGT-induced soil/sediment flux model (DIFS), the labile pool size (Kdl), the rate constants (k1, k-1) of adsorption and desorption and response time (Tc) of SAs in soils were obtained. The addition of manure increased extractable fraction, labile pool size (Kdl) and k1 but decreased k-1. Together with the nonlinear relationship between DGT fluxes and the reciprocal of diffusive layer thickness (Δg), these findings suggested that the release of SAs from soil particles into the soil solution is thermodynamically and kinetically limited, and the manure application could enhance this limitation. This study offers insight into antibiotic availability in soils caused by manure application.

Passive Sampling: A Greener Technique for the 'Dual Carbon' Goal While Implementing the Action Plan for Controlling Emerging Pollutants.

In response to global climate change, China recently committed to achieving peak carbon emissions by 2030 and carbon neutrality by 2060. Carbon emission reduction should be considered in every sector of society including environmental monitoring. As an energy-saving technique in environmental monitoring, passive sampling has many advantages, such as in-situ sampling and a reduction in time/labour requirements. This perspective illustrates the "greenness" of passive sampling techniques, by comparison with traditional sampling methods, and its potential contribution to the 'dual carbon' goal. At the same time passive sampling can provide key support for the action plan for controlling emerging pollutants in China.

Diffusive gradients in thin films (DGT) probe for effectively sampling of per- and polyfluoroalkyl substances in waters and sediments.

The passive sampling technique, diffusive gradients in thin films (DGT) has attracted increasing interests as an in-situ sampler for organic contaminants including per- and polyfluoroalkyl substances (PFAS). However, its effectiveness has been questioned because of the small effective sampling area (3.1 cm2). In this study, we developed a DGT probe for rapid sampling of eight PFAS in waters and applied it to a water-sediment system. It has a much larger sampling area (27 cm2) and as a result lower method quantification limits (0.15 - 0.21 ng/L for one-day deployment and 0.02 - 0.03 ng/L for one-week deployment) and much higher (by > 10 factors) sampling rate (100 mL/day) compared to the standard DGT (piston configuration). The sampler could linearly accumulate PFAS from wastewater, was sensitive enough even for a 24 hr deployment with performance comparable to grab sampling (500 mL). The DGT probe provided homogeneous sampling performance along the large exposure area. The use of the probe to investigate distributions of dissolved PFAS around the sediment-water interface was demonstrated. This work, for the first time, demonstrated that the DGT probe is a promising monitoring tool for trace levels of PFAS and a research tool for studying their distribution, migration, and fate in aquatic environments including the sediment-water interface.

Molecular Dynamics and Light Absorption Properties of Atmospheric Dissolved Organic Matter.

The light-absorbing organic aerosol referred to as brown carbon (BrC) affects the global radiative balance. The linkages between its molecular composition and light absorption properties and how environmental factors influence BrC composition are not well understood. In this study, atmospheric dissolved organic matter (ADOM) in 55 aerosol samples from Guangzhou was characterized using Fourier transform ion cyclotron resonance mass spectrometry and light absorption measurements. The abundant components in ADOM were aliphatics and peptide-likes (in structure), or nitrogen- and sulfur-containing compounds (in elemental composition). The light absorption properties of ADOM were positively correlated with the levels of unsaturated and aromatic structures. Particularly, 17 nitrogen-containing species, which are identified by a random forest, characterized the variation of BrC absorption well. Aggregated boosted tree model and nonmetric multidimensional scaling analysis show that the BrC composition was largely driven by meteorological conditions and anthropogenic activities, among which biomass burning (BB) and OH radical were the two important factors. BrC compounds often accumulate with elevated BB emissions and related secondary processes, whereas the photolysis/photooxidation of BrC usually occurs under high solar radiance/•OH concentration. This study first illuminated how environmental factors influence BrC at the molecular level and provided clues for the molecular-level research of BrC in the future.

Tissue Distribution of Several Series of Cationic Surfactants in Rainbow Trout (Oncorhynchus mykiss) Following Exposure via Water.

Bioaccumulation assessment is important for cationic surfactants in light of their use in a wide variety of consumer products and industrial processes. Because they sorb strongly to natural surfaces and to cell membranes, their bioaccumulation behavior is expected to differ from other classes of chemicals. Divided over two mixtures, we exposed rainbow trout to water containing 10 alkyl amines and 2 quaternary alkylammonium surfactants for 7 days, analyzed different fish tissues for surfactant residues, and calculated the tissues' contribution to fish body burden. Mucus, skin, gills, liver, and muscle each contributed at least 10% of body burden for the majority of the test chemicals. This indicates that both sorption to external surfaces and systemic uptake contribute to bioaccumulation. In contrast to the analogue alkylamine bases, the permanently charged quaternary ammonium compounds accumulated mostly in the gills and were nearly absent in internal tissues, indicating that systemic uptake of the charged form of cationic surfactants is very slow. Muscle-blood distribution coefficients were close to 1 for all alkyl amines, whereas liver-blood distribution coefficients ranged from 13 to 90, suggesting that the dominant considerations for sorption in liver are different from those in blood and muscle. The significant fraction of body burden on external surfaces can have consequences for bioaccumulation assessment.

Deriving in Vivo Bioconcentration Factors of a Mixture of Fragrance Ingredients Using a Single Dietary Exposure and Internal Benchmarking.

Chemicals in mixtures that are hydrophobic with Log KOW > 4 are potentially bioaccumulative. Here, we evaluate an abbreviated and benchmarked in vivo BCF measurement methodology by exposing rainbow trout to a mixture of eight test chemicals found in fragrance substances and three benchmark chemicals (musk xylene (MX), hexachlorobenzene (HCB) and PCB52) via a single contaminated feeding event followed by a 28-day depuration period. Concentrations of HCB and PCB52 in fish did not decline significantly (their apparent depuration rate constants, kT, were close to zero), whereas kT for MX was 0.022 d-1. The test chemicals were eliminated much more rapidly than the benchmark chemicals ( kT > 0.117 d-1). The bioconcentration factors (BCFA) for the test chemicals were in the range of 273 L kg-1 (8-cyclohexadecen-1-one (globanone)) to 1183 L kg-1 (α-pinene); the benchmarked BCFs (BCFG) calculated relative to HCB ranged from 238 L kg-1 (globanone) to 1147 L kg-1 (α-pinene). BCFG were not significantly different from BCFA but had smaller standard errors. BCFs derived here agreed well with values previously measured using the OECD 305 test protocol. We conclude that it will be feasible to derive BCFs of chemicals in mixtures using a single dietary exposure and chemical benchmarking.

DGT Passive Sampling for Quantitative in Situ Measurements of Compounds from Household and Personal Care Products in Waters.

Widespread use of organic chemicals in household and personal care products (HPCPs) and their discharge into aquatic systems means reliable, robust techniques to monitor environmental concentrations are needed. The passive sampling approach of diffusive gradients in thin-films (DGT) is developed here and demonstrated to provide in situ quantitative and time-weighted average (TWA) measurement of these chemicals in waters. The novel technique is developed for HPCPs, including preservatives, antioxidants and disinfectants, by evaluating the performance of different binding agents. Ultrasonic extraction of binding gels in acetonitrile gave good and consistent recoveries for all test chemicals. Uptake by DGT with HLB (hydrophilic-lipophilic-balanced) as the binding agent was relatively independent of pH (3.5-9.5), ionic strength (0.001-0.1 M) and dissolved organic matter (0-20 mg L-1), making it suitable for applications across a wide range of environments. Deployment time and diffusion layer thickness dependence experiments confirmed DGT accumulated chemicals masses are consistent with theoretical predictions. The technique was further tested and applied in the influent and effluent of a wastewater treatment plant. Results were compared with conventional grab-sampling and 24-h-composited samples from autosamplers. DGT provided TWA concentrations over up to 18 days deployment, with minimal effects from biofouling or the diffusive boundary layer. The field application demonstrated advantages of the DGT technique: it gives in situ analyte preconcentration in a simple matrix, with more quantitative measurement of the HPCP analytes.

Desorption kinetics of sulfonamide and trimethoprim antibiotics in soils assessed with diffusive gradients in thin-films.

Although sorption/desorption of antibiotics in soils affects their mobility and availability, with consequences for risks to the surrounding environment, the dynamics of these processes are not well-known. In this study, diffusive gradients in thin-films devices suitable for measuring polar organic compounds (o-DGT) were deployed in two soils for a range of times (5 h to 20 d) to measure the distribution and rates of exchange between solid phase and solution of three sulphonamides (SAs; sulfamethoxazole, SMX; sulfamethazine, SMZ; and sulfadimethoxine, SDM) and trimethoprim (TMP). o-DGT continuously removes antibiotics to a XAD gel layer after passage through a well-defined diffusion layer and therefore perturbs their concentration in the adjacent soil solution. This induces a remobilization flux from the solid phase, which is related to the concentration of antibiotics in the soil solution, their diffusional supply, and the exchange kinetics between dissolved and sorbed antibiotics. A dynamic model of solute interactions called DIFS (DGT induced fluxes in soils) was used to derive distribution coefficients for labile antibiotics (Kdl) and the rate constant for supply of antibiotics from solid phase to solution, expressed as a response time (Tc). Larger labile solid phase pools were observed for TMP than for SAs. The soils could resupply TMP so rapidly that in one soil, where Tc = 2 min, supply was controlled by diffusion. Response times for SAs were generally longer (>27 min), particularly for SDM (>3 h), implying that the supply of SAs to o-DGT samplers was limited by the desorption release rate. A wider implication of this study is that similar solid phase release kinetics may control the uptake of antibiotics by biota.

POPs identification using simple low-code machine learning.

Effectively identifying persistent organic pollutants (POPs) with extensive organic chemical datasets poses a formidable challenge but is of utmost importance. Leveraging machine learning techniques can enhance this process, but previous models often demanded advanced programming skills and high-end computing resources. In this study, we harnessed the simplicity of PyCaret, a Python-based package, to construct machine-learning models for POP screening based on 2D molecular descriptors. We compared the performance of these models against a deep convolutional neural network (DCNN) model. Utilising minimal Python code, we generated several models that exhibited superior or comparable performance to the DCNN. The most outstanding performer, the Light Gradient Boosting Machine (LGBM), achieved an accuracy of 96.20 %, an AUC of 97.70 %, and an F1 score of 82.58 %. This model outshone the DCNN model. Furthermore, it excelled in identifying POPs within the REACH PBT and compiled industrial chemical lists. Our findings highlight the accessibility and simplicity of PyCaret, requiring only a few lines of code, rendering it suitable for non-computing professionals in environmental sciences. The ability of low code machine learning tools (e.g. PyCaret) to facilitate model comparison and interpretation holds promise, encouraging prompt assessment and management of chemical substances.

Bisphenol A sorption on commercial polyvinyl chloride microplastics: Effects of UV-aging, biofilm colonization and additives on plastic behaviour in the environment.

Chemical additives are important components in commercial microplastics and their leaching behaviour has been widely studied. However, little is known about the potential effect of additives on the adsorption/desorption behaviour of pollutants on microplastics and their subsequent role as vectors for pollutant transport in the environment. In this study, two types of commercial polyvinyl chloride (PVC1 and PVC2) microplastics were aged by UV irradiation and biotic modification via biofilm colonization to investigate the adsorption and desorption behaviour of bisphenol A (BPA). Surface cracks and new functional groups (e.g., O-H) were found on PVC1 after UV irradiation, which increased available adsorption sites and enhanced H‒bonding interaction, resulting in an adsorption capacity increase from 1.28 µg/L to 1.85 µg/L. However, the adsorption and desorption capacity not showed significant changes for PVC2, which might be related to the few characteristic changes after UV aging with the protection of light stabilizers and antioxidants. The adsorption capacity ranged from 1.28 µg/L to 2.06 µg/L for PVC1 and PVC2 microplastics, and increased to 1.62 µg/L-2.95 µg/L after colonization by biofilms. The increased adsorption ability might be related to the N-H functional group, amide groups generated by microorganisms enhancing the affinity for BPA. The opposite effect was observed for desorption. Plasticizers can be metabolized during biofilm formation processes and might play an important role in microorganism colonization. In addition, antioxidants and UV stabilizers might also indirectly influence the colonization of microorganisms' on microplastics by controlling the degree to which PVC microplastics age under UV. The amount of biomass loading on the microplastics would further alter the adsorption/desorption behaviour of contaminants. This study provides important new insights into the evaluation of the fate of plastic particles in natural environments.

Unveiling per- and polyfluoroalkyl substance contamination in Chinese paper products and assessing their exposure risk.

The contamination characteristics, migration patterns and health risks of per- and polyfluoroalkyl substances (PFAS) were investigated in 66 Chinese paper products by using target and non-target screening methods. Among 57 target PFASs, 5 and 6 PFASs were found in the hygiene paper products (

Short-term mass loads of per- and polyfluoroalkyl substances in a wastewater treatment plant from South China.

Wastewater treatment plants (WWTPs) are one of the most important sources and sinks for per- and polyfluoroalkyl substances (PFAS). However, limited studies have evaluated short-term temporal variability of PFAS in WWTPs, particularly for their intra-day variations. For this purpose, a time-composite sampling campaign was carried out at a WWTP influent from South China for 1 week. Five out of ten PFAS were found in the influent, i.e., perfluoroheptanoic acid (PFHpA), perfluorooctanoic acid (PFOA), perfluorononanoic acid (PFNA), perfluorobutane sulfonic acid (PFBS), and perfluorooctanesulfonic acid (PFOS). PFOA was the most domain PFAS whereas PFOS was detected occasionally, which might be associated with the prohibition of PFOS use in China. For the first time, we observed significant intra-day fluctuations in mass fluxes for PFOS. Different from a morning peak of pharmaceuticals reported previously, PFOS mass loads fluctuated sharply at noon and night on the weekdays. Furthermore, the mass fluxes of PFOA on the weekend were significantly elevated. For the other PFAS detected, no significant diurnal variations in mass loads were identified. Correlation analysis indicated that domestic activities (e.g., home cleaning) are likely to be the major source of these perfluorocarboxylic acids especially PFOA. In addition, flow fluxes had little effects on these PFAS mass load. These results can aid in future sampling campaigns and optimizing removal strategies for PFAS in wastewater.

Predicting the new psychoactive substance activity of antitussives and evaluating their ecotoxicity to fish.

The misuse of antitussives preparations is a continuing problem in the world, and imply that they might have potential new psychoactive substances (NPS) activity. However, few study focus on their ecological toxicity towards fish. In the present study, the machine learning (ML) methods gcForest and random forest (RF) were employed to predict NPS activity in 30 antitussives. The potential toxic target, mode of action (MOA), acute toxicity and chronic toxicity to fish were further investigated. The results showed that both gcForest and RF achieved optimal performance when utilizing combined features of molecular fingerprint (MF) and molecular descriptor (MD), with area under the curve (AUC) = 0.99, accuracy >0.94 and f1 score > 0.94, and were applied to screen the NPS activity in antitussives. A total of 15 antitussives exhibited potential NPS activity, including frequently-used substances like codeine and dextromethorphan. The binding affinity of these antitussives with zebrafish dopamine transporter (zDAT) was high, and even surpassing that of some traditional narcotics and NPS. Some antitussives formed hydrogen bonds or salt bridges with aspartate (Asp) 95, tyrosine (Tyr) 171 of zDAT. For the ecotoxicity, the MOA of these 15 antitussives in fish was predicted as narcosis. The prenoxdiazin, pholcodine, codeine, dextromethorphan and dextrorphan exhibited very toxic/toxic to fish. It was necessary to pay close attention to the ecotoxicity of these antitussives. In this study, the integration of ML, molecular docking and ECOSAR approaches are powerful tools for understanding the toxicity profiles and ecological hazards posed by new pollutants.

Screening androgen receptor agonists of fish species using machine learning and molecular model in NORMAN water-relevant list.

Androgen receptor (AR) agonists have strong endocrine disrupting effects in fish. Most studies mainly investigate AR binding capacity using human AR in vitro. However, there is still few methods to rapidly predict AR agonists in aquatic organisms. This study aimed to screen AR agonists of fish species using machine learning and molecular models in water-relevant list from NORMAN, a network of reference laboratories for monitoring contaminants of emerging concern in the environment. In this study, machine learning approaches (e.g., Deep Forest (DF)), Random Forests and artificial neural networks) were applied to predict AR agonists. Zebrafish, fathead minnow, mosquitofish, medaka fish and grass carp are all important aquatic model organisms widely used to evaluate the toxicity of new pollutants, and the molecular models of ARs from these five fish species were constructed to further screen AR agonists using AlphaFold2. The DF method showed the best performances with 0.99 accuracy, 0.97 sensitivity and 1 precision. The Asn705, Gln711, Arg752, and Thr877 residues in human AR and the corresponding sites in ARs from the five fish species were responsible for agonist binding. Overall, 245 substances were predicted as suspect AR agonists in the five fish species, including, certain glucocorticoids, cholesterol metabolites, and cardiovascular drugs in the NORMAN list. Using machine learning and molecular modeling hybrid methods rapidly and accurately screened AR agonists in fish species, and helping evaluate their ecological risk in fish populations.

Evidence and recommendations to support the use of a novel passive water sampler to quantify antibiotics in wastewaters.

A novel passive water sampler (diffusive gradients in thin-films for organics, o-DGT) was previously developed and successfully tested in the laboratory, but has not yet been validated in the field. Here, o-DGT samplers were deployed in the influent and effluent of a typical UK wastewater treatment plant (WWTP); the influent was also sampled with a conventional automatic sampler (Auto) and by grab (Grab) sampling. All the samples were analyzed by LC-MS/MS for 40 target antibiotics (including 16 sulfonamides (SAs), 12 fluoroquinolones, 6 macrolides, 2 ionophores, 2 diaminopyimidines, 1 aminocoumarin, and 1 lincosamide). The diffusion coefficients (D) of these antibiotics in o-DGT, measured in the laboratory, ranged from 0.58 × 10(-06) to 6.24 × 10(-06) cm(2) s(-1). The derived surface area normalized sampling rates (RS/A, 0.54-5.74 mL d(-1) cm(-2)) were comparable with those for another passive sampler called POCIS. Fourteen antibiotics were detected in the actively sampled water samples, with 10 of the 14 detected in o-DGT devices deployed for more than 7 days. Most of the antibiotics detected in o-DGT, except sulfapyridine, were continually accumulated by o-DGT for ~10 days. Deployment for 7 days is recommended to integrate ambient concentrations over time, without risks of reaching capacity and significant biofouling. Diffusive boundary layer (DBL) thickness had less effect on the o-DGT measurement than reported for other passive samplers. The comparison between o-DGT and Auto and Grab samplings showed that o-DGT was more efficient in terms of cost, time, and labor. This study demonstrates for the first time in a real environment that o-DGT is an effective tool for the routine monitoring of antibiotics in wastewaters and provides a powerful approach to studying their occurrence, fate, and behavior in the environment.

Sorption and desorption of bisphenols on commercial plastics and the effect of UV aging.

Plastics gradually degrade in the natural environment from the effect of irradiation, which can change the surface properties of plastics and affect the migration behaviour of pollutants. Up to now, studies on the sorption/desorption behaviour of organic pollutants on aged plastics are still limited. In this study, several types of commercial plastics (polyurethane (PU), polyamide (PA), polyvinyl chloride (PVC), expanded polystyrene (EPS)) were selected to investigate the sorption and release behaviour for four kinds of bisphenols (bisphenol-F, A, B, AP). The results from Raman spectroscopy and scanning electron microscopy (SEM) analysis showed evidence of oxidization and surface cracks of plastics after irradiation. The sorption behaviour for both fresh and aged plastics were dominated by hydrophobicity. In addition, the electrostatic force, H-bonding interaction, and π-π interaction were also the important factors impacting the sorption process. The desorption kinetics behaviour indicates that desorption becomes faster after aging. Hydrophobicity is also an important factor that affects desorption behaviour. This study showed that sorption capacity for most fresh and aged plastics was enhanced by the impact of salinity and dissolved organic matter (DOM). Increased temperature could increase the desorption of bisphenols on both fresh and aged plastics, which illustrated that warm environments would promote more pollutants be released from plastics to water bodies.

Assessing the bioavailability of antibiotics in soil with the diffusive gradients in thin films (DGT).

The diffusive gradients in thin films (DGT) technique is an excellent method for investigating the dynamic processes of antibiotics in soils. However, whether it is applicable in antibiotic bioavailability assessment is yet to be disclosed. This study employed DGT to determine the antibiotic bioavailability in soil, and compared the results with plant uptake, soil solutions, and solvent extraction methods. DGT exhibited predictive capability for plant taking in antibiotics proved by the significant linear relationship between the DGT based concentration (CDGT) and antibiotic concentration in roots and shoots. Although the performance of soil solution was acceptable based on linear relationship analysis, its stability was weaker than DGT. The results based on plant uptake and DGT indicated the bioavailable antibiotic contents in different soils were inconsistent because of the distinct mobility and resupply of sulphonamides and trimethoprim in different soils, as represented by Kd and Rds, which were affected by soil properties. Plant species played an important role in antibiotic uptake and translocation. Antibiotic uptake by plants depends on antibiotic, plant and soil. These results confirmed the capability of DGT in determining antibiotic bioavailability for the first time. This work provided a simple and powerful tool for environmental risk evaluation of antibiotics in soils.

Occurrence, fate and ecological risks of 90 typical emerging contaminants in full-scale textile wastewater treatment plants from a large industrial park in Guangxi, Southwest China.

Recent industrial relocation in China causes lots of environment concerns including risks of emerging contaminants (ECs). Herein, the occurrence, fate, removal and ecological risks of 34 per- and polyfluoroalkyl substances (PFAS), 17 endocrine disrupting chemicals (EDCs), 16 phthalate esters (PAEs), and 23 polycyclic aromatic hydrocarbons (PAHs) were investigated in two textile WWTPs (conventional and Fenton-modified) from a large textile industrial park in Southwest China. Totally 50 ECs were identified and the levels followed the order of PAEs > EDCs > PFAS ≈ PAHs. The EDCs predominated in textile washing and rinsing wastewater whereas the PAEs did in desizing wastewater. Biphasic correlations of log Kd and log P, molecular weight, and numbers of rings (r2 = 0.63-0.66, p < 0.01) were observed for PAHs, suggesting that hydrophobicity might not facilitate adsorption of super-hydrophobic PAHs onto activated sludge. 63-69% of detected ECs were effectively removed by two textile WWTPs with removal efficiencies ≥ 80%, which were much higher than previous reports. Fenton processing enhanced the removal efficiencies for long-chain PFAS rather than short-chain PFAS. The PAEs and EDCs posed a medium-to-high risk to aquatic organisms and were screened as the priority ECs. To date, such a comprehensive investigation for ECs has not been previously conducted in textile WWTPs and this study provides basic information about regional chemical emission inventory of ECs.