Metal pollution in the Pearl River Estuary and implications for estuary management: The influence of hydrological connectivity associated with estuarine mixing

Heavy Metal Accumulation in Oysters from an Aquaculture Area in the Luoyangjiang River Estuary

Oysters are a group of economically important bivalves in China, with estuaries serving as one of their primary cultivation areas. However, heavy metal pollution in these estuarine environments poses a potential threat to aquaculture by leading to the accumulation of heavy metals in farmed oysters, which could impact their safety and marketability. This study was conducted in the aquaculture area of the Luoyangjiang River estuary, where eight sampling sites were selected. Water, sediment, and oysters categorized by shell length were collected from each site. The concentrations of heavy metals (Ag, As, Cd, Cr, Cu, Ni, Pb, and Zn) were determined in both the environmental samples and oyster tissues. Additionally, multiplex species-specific PCR was used to identify oyster species. The results showed significant variations in dissolved-phase and suspended particulate matter (SPM) metal concentrations across different sampling sites, while sediment metal concentrations were more consistent but similar to those in SPM. The large oysters were comprised of 50% Magallana angulata and 50% Magallana gigas, while small oysters were identified as Magallana sikamea. The Cd, Cu, Pb, and Zn levels in both size groups of oysters exceeded data from previous studies, indicating contamination in the estuary. The observed differences in heavy metal concentrations between large and small oysters primarily reflect species-specific variability in metal accumulation, which may also be influenced by factors such as growth and exposure duration. Furthermore, the lack of significant correlation between metal concentrations in environmental media and oysters suggests that oysters may be exposed to multiple sources of metal contamination.

Marine coastal chemistry related to inland inputs in San Jorge gulf and the adjacent north coast

Coastal regions are sectors where human activities impact the marine ecosystem, and if necessary control measures are not taken, they can generate negative consequences for health and ecosystem services. Within the framework of the Pampa Azul initiative and under the One Health paradigm, the interconnection between the terrestrial and marine environments of the San Jorge Gulf and the adjacent north coast has been studied. In November, 2022, a campaign was carried out aboard the R/V "Mar Argentino" at thirty-four stations near the coast. There, for the first time, simultaneously with in-situ measurements of physical variables, macronutrients (NO3-, PO4-3, Si(OH)4 and NH4+), particulate silica (BSi and LSi), trace metals in the particulate material (Cd, Cu, Cr, Fe and Pb) and the phytoplankton community were analyzed. The results showed a high nutrient dynamic, with a significant influence of natural stratification and anthropogenic condition due to the discharge of effluents off the cities of Comodoro Rivadavia and Caleta Olivia. Under natural conditions, NO3- and Si(OH)4 limited the surface primary production by 47 % and 41 %, respectively. Additionally, due to the anthropogenic contribution, NH4+ concentration reached 3 μM, increasing the proliferation of nanophytoplankton, among other consequences. As a result of nutrient dynamics, the uptake of Si(OH)4, the growth rate of diatoms and their production of BSi were decoupled. Furthermore, a significant correlation between LSi and Fe in particulate matter was evidenced, opening new lines of research that relate dust storms to primary productivity in this marine environment. The measured concentrations of trace metals do not appear to be a biological risk; however, contamination by Cd (37.6 μg g- 1 d.w.) and Cu (214.97 μg g- 1 d.w.) off Camarones poses a significant concern that must be addressed in the immediate future.

Insights into pollution characteristics and human health risks of plasticizer phthalate esters in shellfish species

The ubiquitous application of phthalate esters (PAEs) as plasticizers contributes to high levels of marine pollution, yet the contamination patterns of PAEs in various shellfish species remain unknown. The objective of this research is to provide the first information on the pollution characteristics of 16 PAEs in different shellfish species from the Pearl River Delta (PRD), South China, and associated health risks. Among the 16 analyzed PAEs, 13 were identified in the shellfish, with total PAE concentrations ranging from 23.07 to 3794.08 ng/g dw (mean = 514.35 ng/g dw). The PAE pollution levels in the five shellfish species were as follows: Ostreidae (mean = 1064.12 ng/g dw) > Mytilus edulis (mean = 509.88 ng/g dw) > Babylonia areolate (mean = 458.14 ng/g dw) > Mactra chinensis (mean = 378.90 ng/g dw) > Haliotis diversicolor (mean = 335.28 ng/g dw). Dimethyl phthalate (DMP, mean = 69.85 ng/g dw), diisobutyl phthalate (DIBP, mean = 41.39 ng/g dw), dibutyl phthalate (DBP, mean = 130.91 ng/g dw), and di(2-ethylhexyl) phthalate (DEHP, mean = 226.23 ng/g dw) were the most abundant congeners. Notably, DEHP constituted the most predominant fraction (43.98 %) of the 13 PAEs detected in all shellfish from the PRD. Principal component analysis indicated that industrial and domestic emissions served as main sources for the PAE pollution in shellfish from the PRD. It was estimated that the daily intake of PAEs via shellfish consumption among adults and children ranged from 0.004 to 1.27 μg/kgbw/day, without obvious non-cancer risks (< 0.034), but the cancer risks raised some alarm (2.0 × 10-9-1.4 × 10-5). These findings highlight the necessity of focusing on marine environmental pollutants and emphasize the importance of ongoing monitoring of PAE contamination in seafood.

Low salinity influences the dose-dependent transcriptomic responses of oysters to cadmium

Species in estuaries tend to undergo both cadmium (Cd) and low salinity stress. However, how low salinity affects the Cd toxicity has not been fully understood. Investigating the impacts of low salinity on the dose-response relationships between Cd and biological endpoints has potential to enhance our understanding of the combined effects of low salinity and Cd. In this work, changes in the transcriptomes of Pacific oysters were analyzed following exposure to Cd (5, 20, 80 μg/L Cd2+) under normal (31.4 psu) and low (15.7 psu) salinity conditions, and then the dose-response relationship between Cd and transcriptome was characterized in a high-throughput manner. The benchmark dose (BMD) of gene expression, as a point of departure (POD), was also calculated based on the fitted dose-response model. We found that low salinity treatment significantly influenced the dose-response relationships between Cd and transcripts in oysters indicated by altered dose-response curves. In details, a total of 219 DEGs were commonly fitted to best models under both normal and low salinity conditions. Nearly three quarters of dose-response curves varied with salinity condition. Some monotonic dose-response curves in normal salinity condition even were replaced by nonmonotonic curves in low salinity condition. Low salinity treatment decreased the PODs of differentially expressed genes induced by Cd, suggesting that gene differential expression was more prone to being triggered by Cd in low salinity condition. The changed sensitivity to Cd in low salinity condition should be taken into consideration when using oyster as an indicator to assess the ecological risk of Cd pollution in estuaries.

A microcosm evaluation of metal cycling in an urbanized contaminated estuary varying with oxic-hypoxic-anoxic-reoxic transition: Behavior, fluxes, and mechanism

Water hypoxia and metal pollution are commonly co-existed in urbanized estuaries. This study focuses on the effect of an extended dissolved oxygen (DO) full-life dynamics (86 days) on metal behavior across the sediment-water interface through laboratory microcosms from two typical zones in Pearl River Estuary. Combining our time-series results of concentrations and fluxes, it showed that Co, Ni, and Zn consistently presented a release-precipitation-release trajectory with an oxic-hypoxic-anoxic-reoxic transition, characterized with highly variable behavior in the hypoxic-anoxic hotmoments. In parallel, changing DO dynamics significantly activated a repartitioning process of Co, Ni, and Zn among several species and elevated their risk in sediments, promoting the formation of more labile species in the 0-10 mm hotspots, where metals sensitively responded. Over DO transition, metal cycling was tightly co-related with Fe, Mn, and S elements. It was found that Mn was dominated in low oxygen-hypoxic period, but switched to S and Fe in anoxic stage, limiting sustained metal liberation to overlying water. Enlarging this experiment to practice, released Zn fluxes from sediments in hypoxic summer could contribute about ∼2.0% to their stocks in water column, while increase to 20% (1 m bottom water) in highly-stratified zones. This study has certain significance in understanding the long-term metal behavior and fate in estuarine regions, even lakes and reservoirs.

Endocrine-Disrupting Chemical Exposure Induces Adverse Effects on the Population Dynamics of the Indo-Pacific Humpback Dolphin

Cetaceans play a pivotal role in maintaining the ecological equilibrium of ocean ecosystems. However, their populations are under global threat from environmental contaminants. Various high levels of endocrine-disrupting chemicals (EDCs) have been detected in cetaceans from the South China Sea, such as the Indo-Pacific humpback dolphins in the Pearl River Estuary (PRE), suggesting potential health risks, while the impacts of endocrine disruptors on the dolphin population remain unclear. This study aims to synthesize the population dynamics of the humpback dolphins in the PRE and their profiles of EDC contaminants from 2005 to 2019, investigating the potential role of EDCs in the population dynamics of humpback dolphins. Our comprehensive analysis indicates a sustained decline in the PRE humpback dolphin population, posing a significant risk of extinction. Variations in sex hormones induced by EDC exposure could potentially impact birth rates, further contributing to the population decline. Anthropogenic activities consistently emerge as the most significant stressor, ranking highest in importance. Conventional EDCs demonstrate more pronounced impacts on the population compared to emerging compounds. Among the conventional pollutants, DDTs take precedence, followed by zinc and chromium. The most impactful emerging EDCs are identified as alkylphenols. Notably, as the profile of EDCs changes, the significance of conventional pollutants may give way to emerging EDCs, presenting a continued challenge to the viability of the humpback dolphin population.

Historical construction, quantitative source identification and risk assessment of heavy metals contamination in sediments from the Pearl River Estuary, South China

Historical reconstruction of heavy metals (HMs) contamination in sediments is a key for understanding the effects of anthropogenic stresses on water bodies and predicting the variation trends of environmental state. In this work, eighteen sediment cores from the Pearl River Estuary (PRE) were collected to determine concentrations and geochemical fractions of HMs. Then, their potential sources and the relative contributions during different time periods were quantitatively identified by integrating lead-210 (210Pb) radioisotope dating technique into positive matrix factorisation (PMF) method. Pollution levels and potential ecological risks (PERs) caused by HMs were accurately assessed by enrichment factors (EF) based on establishment of their geochemical baselines (GCBs) and multiparameter evaluation index (MPE). HMs concentrations generally showed a particle size- and organic matter-dependent distribution pattern. During the period of 1958-1978, HMs concentrations remained at low levels with agricultural activities and natural processes being identified as the predominant sources and averagely contributing >60%. Since the reform and opening-up in 1978, industrial and traffic factors become the primary anthropogenic sources of HMs (such as Cu, Zn, Cd, Pb, Cr, and Ni), averagely increasing from 22.1% to 28.1% and from 11.6% to 23.4%, respectively. Conversely, the contributions of agricultural and natural factors decreased from 37.0% to 28.5% and from 29.3% to 20.0%, respectively. Subsequently, implementation of environmental preservation policies was mainly responsible for the declining trend of HMs after 2010. Little enrichment of sediment Cu, Zn, Pb, Cr and Ni with EFs (0.15-1.43) was found in the PRE, whereas EFs of Cd (1.16-2.70) demonstrated a slight to moderate enrichment. MPE indices of Cu (50.7-252), Pb (52.0-147), Zn (35.5-130), Ni (19.6-71.5), Cr (14.2-68.8) and Cd (0-9.90) highlighted their potential ecological hazards due to their non-residual fractions and anthropogenic sources.

Assessment of water quality of Netravathi-Gurupur estuary, India through chemometric approach for fisheries sustainability

The present study aimed to assess the water quality dynamics (physicochemical properties, nutrient and chlorophyll-a) in the Netravathi-Gurupur estuary, India across the space and seasons and to simplify the complex water quality dataset through a chemometric approach. The results indicated that pH, EC, temperature, alkalinity, salinity, TDS, hardness, chloride and intense mixing of suspended solids, driven by the semidiurnal tides, are the major factors affecting water quality. Spatial heterogeneity and monsoon have profound impact on nutrient distribution revealing the following values (mg l-1): phosphate-P (0.015-0.105), nitrate-N (0.016-0.094), nitrite-N (0.001-0.012), and silicate (1.83-14.50). The estuary was evaluated for suitability for brackish water fisheries. The results indicated fair water quality during pre- and post-monsoon but marginal quality in monsoon, primarily due to dilution associated with reduced salinity. The outcome of this study can be suitably utilized for the sustainable development of estuaries and their feasibility for brackish water fisheries.

Unravelling the life-history patterns and habitat preferences of the Japanese eel (Anguilla japonica) in the Pearl River, China

Eels have fascinated biologists for centuries due to their amazing long-distance migrations between freshwater habitats and very distant ocean spawning areas. The migratory life histories of the Japanese eel, Anguilla japonica, in the waters of south China are not very clear despite its ecological importance, and the need for fishery regulation and management. In this study, strontium (Sr) and calcium (Ca) microchemical profiles of the otoliths of silver eels were measured by X-ray electron probe microanalysis based on data collected from different habitats (including freshwater and brackish habitats), in the large subtropical Pearl River. The corresponding habitat preference characteristics were further analysed using redundancy analysis (RDA). A total of 195 Japanese eels were collected over 6 years. The collected individuals ranged from 180 to 771 mm in total length and from 8 to 612 g in body weight. Two-dimensional pictures of the Sr:Ca concentrations in otoliths revealed that the A. japonica in the Pearl River are almost entirely river eels, spending the majority of their lives in fresh water without exposure to salt water, while the catadromous migration time has delayed about 1 month in the Pearl River estuary in the past 20 years. RDA analysis further indicated that juveniles and adults preferred water with high salinity and high tide levels. Youth preferred habitats with high river fractals. Our findings contribute to a growing body of evidence showing that the eels are extremely scarce currently and conservation measures against them are imminent, including the protection of brackish and freshwater areas where they live in south China.

Metagenomics-Based Microbial Ecological Community Threshold and Indicators of Anthropogenic Disturbances in Estuarine Sediments

Assessing the impacts of cumulative anthropogenic disturbances on estuarine ecosystem health is challenging. Using spatially distributed sediments from the Pearl River Estuary (PRE) in southern China, which are significantly influenced by anthropogenic activities, we demonstrated that metagenomics-based surveillance of benthic microbial communities is a robust approach to assess anthropogenic impacts on estuarine benthic ecosystems. Correlational and threshold analyses between microbial compositions and environmental conditions indicated that anthropogenic disturbances in the PRE sediments drove the taxonomic and functional variations in the benthic microbial communities. An ecological community threshold of anthropogenic disturbances was identified, which delineated the PRE sediments into two groups (H and L) with distinct taxa and functional traits. Group H, located nearshore and subjected to a higher level of anthropogenic disturbances, was enriched with pollutant degraders, putative human pathogens, fecal pollution indicators, and functional traits related to stress tolerance. In contrast, Group L, located offshore and subjected to a lower level of anthropogenic disturbances, was enriched with halotolerant and oligotrophic taxa and functional traits related to growth and resource acquisition. The machine learning random forest model identified a number of taxonomic and functional indicators that could differentiate PRE sediments between Groups H and L. The identified ecological community threshold and microbial indicators highlight the utility of metagenomics-based microbial surveillance in assessing the adverse impacts of anthropogenic disturbances in estuarine sediments, which can assist environmental management to better protect ecosystem health.

Assessment of metal bioavailability in sediments and bioaccumulation in edible bivalves, and phyto-remediation potential of mangrove plants in the tropical (Kali) estuary, India

The bioavailability of metals (Fe, Mn, Zn, Cu, Co and Ni) in sediment cores (K-1, K-2, K-3 and K-4) and bioaccumulation in edible bivalves were studied to determine the toxicity of metals in the Kali Estuary, India. Enrichment Factor (EF) construed anthropogenic sources of Zn, Co and Ni, while Geo-accumulation Index (Igeo) revealed pollution of Zn and Ni based on total metal analysis. The Pollution Load Index (PLI >1) supported anthropogenic origin of metals in estuary. Metal speciation study indicated bioavailability of metals in sediments. The bioavailable Mn and Co equalled/exceeded the Apparent Effect Threshold (AET) limit (cores K-1, K-3 and K-4) and indicated toxicity to estuarine biota. The metals in Metetrix casta (Fe, Mn, Zn, Cu, and Ni), Saccostrea cucullata (Fe, Mn, Zn, Cu, Co and Ni) and Villorita cyprinoides (Fe, Mn, Zn, Cu, Co and Ni) exceeded the permissible bioaccumulation limit. Thus, revealed metal toxicity to bivalves and labelled them un-safe for human consumption. Translocation Factor (TF > 1) indicated the use of Kandelia candel in phyto-remediation of Fe, Zn, Cu, Co and Ni at station K-3, and Sonnaretia caseoloris in phyto-remediation of Fe, Zn and Ni at station K-4.

Identification of key features relating to the coexistence mechanisms of trace elements and per- and polyfluoroalkyl substances (PFASs) in marine mammals

Organic and inorganic substances coexist in the livers of marine mammals and may correlate with one another; however, their coexistence mechanisms and relevant key features remain largely unknown. In this study, temporal variations (2011-2021) in the concentrations of nine trace elements and 19 per- and polyfluoroalkyl substances (PFASs) in the livers of Indo-Pacific humpback dolphins (Sousa chinensis) and finless porpoises (Neophocaena phocaenoides) were investigated. Interannual Cd in dolphins increased significantly whereas Pb concentrations decreased over the past decade (p < 0.05). Interannual levels of seven and four PFASs in dolphins and porpoises decreased significantly with time (p < 0.05). By further extending the timescale to 1993-2021, the sensitivity of trace elements to annual change further increased, whereas the sensitivity of PFASs remained relatively stable. Cu levels, similar to the majority of PFASs, were negatively correlated with the body length of the studied cetaceans, which led to positive correlations of Cu with six long-chain perfluoroalkyl carboxylic acids, perfluorodecane sulfonic acid, and perfluoroethylcyclohexane sulfonic acid. The concentrations of trace elements in the cetacean liver were closely correlated with cetacean sex, species, and body length, whereas PFAS concentration was responsive to time-related features such as stranded season and year. By further employing a machine learning method, we demonstrated that body length and a time-related factor (year) played a crucial role in predicting the concentrations of certain trace elements and PFASs, respectively, particularly Cu and perfluoroheptanoic acid.

Characterization of metal resistance genes carried by waterborne free-living and particle-attached bacteria in the Pearl River Estuary

Toxic metals can substantially change the bacterial community and functions thereof in aquatic environments. Herein, metal resistance genes (MRGs) are the core genetic foundation for microbial responses to the threats of toxic metals. In this study, waterborne bacteria collected from the Pearl River Estuary (PRE) were separated into the free-living bacteria (FLB) and particle-attached bacteria (PAB), and analyzed using metagenomic approaches. MRGs were ubiquitous in the PRE water and mainly related to Cu, Cr, Zn, Cd and Hg. The levels of PAB MRGs in the PRE water ranged from 8.11 × 10⁹ to 9.93 × 10¹² copies/kg, which were significantly higher than those of the FLB (p < 0.01). It could be attributed to a large bacterial population attached on the suspended particulate matters (SPMs), which was evidenced by a significant correlation between the PAB MRGs and 16S rRNA gene levels in the PRE water (p < 0.05). Moreover, the total levels of PAB MRGs were also significantly correlated with those of FLB MRGs in the PRE water. The spatial pattern of MRGs of both FLB and PAB exhibited a declining trend from the low reaches of the PR to the PRE and on to the coastal areas, which was closely related to metal pollution degree. MRGs likely carried by plasmids were also enriched on the SPMs with a range from to 3.85 × 10⁸ to 3.08 × 10¹² copies/kg. MRG profiles and taxonomic composition of the predicted MRG hosts were significantly different between the FLB and PAB in the PRE water. Our results suggested that FLB and PAB could behave differential response to heavy metals in the aquatic environments from the perspective of MRGs.

Impacts of periodic saltwater inundation on heavy metals in soils from the Pearl River Delta, China

The alluvial plain delta is vulnerable to complex land-sea interactions as a result of rising sea levels and increasing storm surges due to global environmental change. Topsoil samples (0-20 cm) from the Pearl River Delta (PRD) were subjected to periodic artificial saltwater inundation treatments with varying salinities (0, 35‰, 40‰, 50‰) for 50 days to explore the impacts of saltwater inundation on heavy metals (Cd, Pb, Zn) in soils. The inundation treatments reached dynamic equilibrium in approximately 20 days, and heavy metals were promoted to release into leachate. The extraction rate of heavy metals was highest with artificial saltwater at 40‰ salinity, which was generally attributed to pH variation, increasing ionic strength and reductive dissolution of Fe-Mn oxyhydroxides. However, once the salinity reached 50‰, a higher SO2- 4concentration could depress the release of heavy metals by providing more negative adsorption sites. Cd in soils was most likely to leach, followed by Zn, but Pb showed higher retention. After saltwater inundation, the bioavailability of heavy metals decreased in the order Cd > Zn > Pb. Redundancy analysis (RDA) results demonstrated that Cd and Zn are more susceptible to soluble salt ions in soils than Pb. The retention of Pb could be attributed to its larger ionic radius and reduced hydrated radius as well as the stable species under the pH conditions of the treatments. This study suggests that the migration of heavy metals could reduce the water quality and increase the ecological risk of the interaction zone between land and sea.

Risk assessment and spatio-temporal distribution of dissolved trace metals in Swarna, Sharavati and Kali estuaries, South-West Coast of India

Trace metals act as a limiting nutrient and prerequisite for primary productivity in marine environments. The distribution of metals in dissolved phase along the salinity gradients of Swarna, Sharavati and Kali estuaries in southwestern India, during post and pre-monsoon seasons, were studied. We have investigated the behaviour of trace metals in the estuarine environment and their extent of impact on human health and ecosystem. The study revealed, non-conservative behaviour of dissolved Mn, Fe, Ni, Cd and Co in the estuaries. Whereas Cu behaved non-conservatively in post-monsoon and conservatively in pre-monsoon seasons. Risk assessment studies revealed that higher chronic daily intake (CDI) in humans, through dermal pathway, in Swarna and Sharavati estuaries during post-monsoon, whereas it was during pre-monsoon season in the Kali estuary. Hazard Index values for the studied metals in adults and children are below risk thresholds, though children are more prone to health risk through the dermal pathway.

Spatiotemporal variation and ecological risk assessment of sediment heavy metals in two hydrologically connected lakes

Excessive accumulation of heavy metals in global lake sediments poses a serious threat to lake water quality and ecosystem security. However, there is still a knowledge gap in comparison of heavy metal variation and pollution in hydrologically connected lakes. In this study, concentrations of As, Cd, Cr, Cu, Hg, Pb, and Zn in sediments of two hydrologically connected lakes, Xingkai Lake and Xiaoxingkai Lake, were determined during the hydrologically connected periods (May and September) and disconnected period (January and July) in 2021. We found the range of As was 2.58∼14.35 mg/kg, Cd was 0.050∼0.21 mg/kg, Cr was 28.58∼262.3 mg/kg, Cu was 3.12∼28.05 mg/kg, Hg was 0.0030∼0.14 mg/kg, Pb was 10.87∼58.86 mg/kg, and Zn was 18.21∼90.73 mg/kg. Heavy metal concentrations were lower than grade I level in Chinese soil quality standards with significant spatial and temporal differences in the basin. Overall, most of the sampling sites in Xingkai Lake and Xiaoxingkai Lake were at the uncontaminated level and moderate ecological risk during the sampling period. Two lakes showed different heavy metal compositions, accompanied by higher contamination level and higher potential ecological risk in the small lake than those in the large lake based on analysis of the geo-accumulation index and potential ecological risk index. Besides, the contamination level and potential ecological risk in May and September were higher than those in January and July, mainly due to human activities and hydrological connectivity. The ecological risks were moderate for Cd and Hg, and low for As, Cr, Cu, Pb, and Zn. Correlation and PCA analyses showed that Cd mainly originated from anthropogenic sources, while other metals mainly came from natural sources. These findings elucidate the effects of agriculture and hydrological connectivity on heavy metals in sediments, and provide scientific basis for the reasonable management of lake ecosystem.

Source-specific ecological risk assessment and quantitative source apportionment of heavy metals in surface sediments of Pearl River Estuary, China

In this study, surface sediments of the Pearl River Estuary were collected from 29 stations and investigated the spatial distribution, pollution level, quantitative source apportionment, and source-specific ecological risk of 10 heavy metals. The mean concentrations followed the order of Mn > Zn > Cr > Cu > Ni > Pb > As > Co > Cd > Hg. In terms of spatial distribution, it showed that the heavy metals were enriched in the inner Pearl River Estuary with 'extremely high' level of Hg, whereas, Cd and Zn posed 'moderate to high' contamination potential. We apportioned four main sources using positive matrix factorization model, in which natural geogenic and industrial manufacturing sources accounted for 36.84% and 27.11% of the total, respectively. However, the source-specific risk assessment suggested that mixed anthropogenic sources were the main contributors, and ecological risks were strongly affected by anthropogenic imports from the surrounding cities.