Environmental Quality and Ecological Risk Assessment of Heavy Metals in the Zhuhai Coast, China

The rapid urbanization and industrialization in the coastal Zhuhai have resulted in severe heavy metal contamination. Heavy metals are extensively used as indicators of coastal contamination, but the comparability among different geographic environments needs to be further explored and scrutinized. In this study, heavy metals (i.e., As, Zn, Cd, Cr, Hg, Cu, and Pb) in surface seawater (20 sites), sediments (10 sites), and organisms (13 sites, shellfish, algae, fish, and crustaceans) were analyzed to assess the distribution, enrichment, and ecological risk in the coastal regions of Zhuhai (2017 and 2018). The results showed that Cu, Zn, and Pb exceeded the standards of seawater, indicating the presence of anthropogenic pollution. Meanwhile, the ecological risk results showed that Cd and Cu indicated moderately high environmental risks. The quantity and quality of heavy metals in organisms showed different trends in the sea, which affected marine organisms in the northeastern area. The main factors contributing to the ecological risk of these two heavy metals are the emissions from chemical companies. Therefore, preserving the high biological diversity in this area is important to ensure the health and productivity of the seawater ecosystem in the coastal area. These results may play a key role in promoting decision makers to manage the sustainability of the Zhuhai coast.


INTRODUCTION
Heavy metals are one of the most serious pollutants in the environment and have attracted widespread attention worldwide due to their inherent toxicity, persistence, and bioaccumulation properties (Abaya et al., 2018;Huang et al., 2019;Kong et al., 2021). With industrial, agricultural, and urban development, heavy metals discharged into the coastal marine environment through surface runoff are increasing and enriching in marine sediments and aquatic organisms (Sruthi et al., 2018). When the concentration of heavy metals in the coastal ecosystem reaches a certain level, it will pose a great danger to the marine ecosystem and human health . The coastal ecosystem is a very important link in the ecological protection of coastal areas and plays an important role in the marine environment, social development, and population growth of coastal areas (Majithiya et al., 2018;Rao et al., 2021). Coastal ecosystems can also provide ecological buffers to the coastal areas and protect them from external intrusion and hazards (Sruthi et al., 2018). Changes in coastal zone ecosystems are partly due to sea-level rise, climate change, and partly human activities, such as changing land properties for coastal development, aquaculture, agricultural, and industrial activities Wang L. et al., 2019).
Ecological risk assessment is developed from risk assessment, and in recent decades, its evaluation content, scope, and research scale have been greatly expanded. In general, the assessment tools can be used to qualitatively evaluate the degree to which the seawater-or sediment-associated chemical status might adversely affect aquatic organisms and are designed to assist in the interpretation of water/sediment quality (Kaikkonen et al., 2018;Quyen et al., 2021). The research scale of ecological risk assessment has been extended from single species to ecosystems, and the research objects have been extended from terrestrial to marine ecosystems (Wang Q. et al., 2019). The species sensitivity distribution (SSD) method, on the basis of differences in the sensitivity of different species to contaminations, was employed to develop environmental criteria for the conservation. The SSD calculates the toxicological data of different species for the whole ecosystem so as to understand the hazards of pollutants at the ecosystem level (Xia et al., 2020). The SSD method has become one of the advantageous methods of risk assessment internationally in recent years due to its advantages of conciseness and clear ecological significance Xu et al., 2021). However, only a few studies on this method have been reported in China and mainly focused on freshwater organisms and assessed the ecological risk of various organisms such as DDT, lindane, polycyclic aromatic hydrocarbons, chlorophenols, and diisooctyl phthalate to freshwater organisms (Li et al., 2018a;Cui et al., 2021;Xu et al., 2022). Ecological risk assessment based on species sensitivity differences has obvious advantages for solving complex ecosystems (Sruthi et al., 2018), such as integrated ecological risk classification in coastal seas (Kumar et al., 2019;Luo et al., 2019).
The Pearl River is the second most important river system in South China, ranking among the seven major rivers in China in terms of annual runoff. The stability and health of the ecological environment of the Zhuhai coast, which is one of the most important sea outlets in the Zhuhai watershed, are related to the healthy development of this special economic zone. In previous studies, large-scale data on the heavy metal concentration and ecological risk of coastal seawater and sediments along the coast of Guangdong province have been collected, especially focusing on the economically developed Pearl River estuary (Birch, 2018;. Previous studies found that the characterization of the distribution of heavy metals in these coastal marine areas has certain regularity. From far to near the coast, the concentration of heavy metals is gradually increasing, and the resulting risk is also gradually increasing (Table S1). Usually, the coastal special economic zones are experiencing high economic development, which inevitably leads to the discharge of industrial pollution. Analyzed from the environmental-geochemical point of view, the accumulation of heavy metals is gradually weakened with the extension of migration paths, which is the main reason for the serious excess of heavy metals on the coast of this study area.
The ecological risk studies based on SSD have been carried out less for the sea area of the Zhuhai coast. Therefore, in this study, the seawater, sediment, biological species, and productivity were mainly sampled and investigated by selecting the sea area of Gaolan on the Pearl River coast as the study area.
In China, the marine areas with the most serious heavy metal pollution are the bays and estuaries where industrial and mining industries are concentrated, some of them are contaminated with Cd, Cu, Pb, and Hg, and some of them are seriously contaminated with Cd and As (Fang et al., 2018;Hu et al., 2021). Heavy metal pollution in the ocean is bound to bring a great threat to the local ecological diversity (Reddy et al., 2016;Satheeswaran et al., 2019;Pandey et al., 2021;Pandey et al., 2022). Therefore, it is urgent to study the ecological risk of heavy metal pollution in China's oceans. Hence, the target of this research is to assess the potential impact of heavy metals on seawater, sediment, and marine organisms and provide ecological risk evaluation and risk management of heavy metal pollution in China's marine environment. Thus the SSD curves of heavy metals, the 5% hazard concentration (HC5) of marine organisms at different trophic levels, and the potential impact ratio (PAF) were calculated to support the completion of the above goals.

Description of Study Area
Zhuhai (21°45′-22°40′N, 113°03′-113°35′E) is an economic functional area with a total area of 380 km 2 and is the main coastal hub port in South China (Figure 1). Zhuhai and its coastal sea areas have a typical subtropical monsoon climate with an obvious winter monsoon. From 1962 to 2010, the linear warming rate of the annual average temperature in Zhuhai was 0.14°C/10a, the average temperature was 22.61°C (Li and Chen, 2014), and the SD was 0.44. The depth of seawater is about −8 to 10 m, which provides excellent conditions for industrial and agricultural production. As a national-level economic and technological development zone in China, Zhuhai is the gateway to the world from the South China Sea, and the Zhuhai coast is one of the most dynamic and potential areas of Guangdong's marine economy, with the largest berths of liquid chemicals and bulk cargo in the Pearl River Delta. The industrial types in Zhuhai are mainly large equipment manufacturing, petrochemical industry, and fine chemical industry, and aquaculture is also an important contributor to economic development.

Sampling Collection and Analysis
The location of the study area and the layout of the survey stations are shown in Figure 1. A total of 20 water quality survey sites, 10 sediment survey sites, and 13 ecological survey sites (partially overlapping, totaling 20 sampling sites) were studied in this survey. Sample collection and survey studies were conducted in November 2017 (Autumn) and March 2018 (Spring). Samples were collected only once per season. On each of the sampling sites, four subsamples were taken by randomly collecting fresh surfaces around the center sampling site (5 m). Combined with the center sample, five samples are collected and mixed thoroughly into one composite sample for post-processing analysis. Following the collection, all samples were stored in a refrigerator and transported to the laboratory by a professional transporter. The People's Republic of China for Marine Survey Specifications (GB/ T 12763-2007) and Marine Monitoring Specifications (GB17378-2007) were chosen as the references for steps and tools in the sample collection process (Li et al., 2018a;Liu et al., 2018;. The concentration of heavy metals was determined by inductively coupled plasma mass spectrometry and inductively coupled plasma atomic emission spectrometry. To ensure the quality of analysis, the control methods and laboratory quality assurance were implemented, including the use of standard operating procedures, analysis of reagent blanks, calibration with standards, and analysis of replicates. The relative SD indicating the precision of analytical procedures ranged from 3.5% to 8.5%. The percent recoveries varied between 86% and 95% for all metals. The seawater depth, transparency, and color were observed on site. Surface and bottom samples were collected at 0.1-1 m from the seawater surface and 2 m from the bottom. Sediment samples were collected using a 0.05-m 2 grab mud collector to collect sediment samples at 0-1 cm from the surface layer of seafloor sediments and were stored in sealed packages. After stones and other debris were removed, samples were air-dried at room temperature (20°C to 23°C). Then the samples were ground with an agate mortar and passed through 160 mesh nylon screens. Four categories of organisms, including locally common and representative shellfish, algae, fish, and crustaceans, were selected among intertidal organisms and benthic organisms for the quality survey, and the samples were collected, preserved, and transported according to the relevant regulations of GB/T 12763.6-2007 for the biological survey (Kim et al., 2018;Wu et al., 2018). A trawling survey was carried out on the sampling site, the net was set at a distance of 2 nautical miles from the station, and the towing speed was controlled at 2-3 knots. A sufficient number (about 1.5 kg) of intact organisms samples were selected and placed in a clean polyethylene bag for sealing and preservation. This bag along with the sample label in another polyethylene bag was sealed and refrigerated. The collected water samples were filtered using a 0.45-mm or 60mm microporous membrane with HNO 3 /H 2 SO 4 to pH < 2 for cryogenic refrigeration, the collected sediment samples were sealed and stored away from light, and the samples of various organisms were sorted and stored in polyethylene bags for freezing. As and Hg were determined by atomic fluorescence, while Cu, Cd, Cr, Pb, and Zn were determined by flameless atomic absorption spectrophotometry Anusuya et al., 2021).

Data Analysis
Several statistical tools, mainly including Pearson's correlation analysis, non-metric multidimensional scaling analysis (NMDS), Adonis test, Procrustes analysis, and constrained ordination principal coordinate analysis (CPCOA), were used to investigate the correlations or differences between heavy metals and samples. Network analysis was employed to explore the cooccurrence of heavy metals and organisms in marine based on their strong and significant correlation matrixes (p-value <0.05; Spearman's r > 0.8). Network visualization was performed employing Gephi (v0.9.2). Statistical analyses were conducted in R (v3.4.2) by utilizing several R packages, including gunifrac, car, and vegan. Kriging interpolation method of ArcGIS (v10.8, ESRI, Redlands, CA, USA) was employed for the spatial distribution characteristics of heavy metals in the sediment.

Partitioning and Bioaccumulation
The distribution coefficient between sediment and water (K d ), biotasediment accumulation factor (BSAF), and bioaccumulation factor (BAF) were employed to explore the accumulation and adsorption of heavy metals between different phases. The specific formulas were calculated as follows (Ashayeri and Keshavarzi, 2019;Liu et al., 2022): where C sed , C sea , and C org are the concentration of target heavy metal in sediments, seawaters, and marine organisms, respectively. The value of K d shows the possibilities of heavy metals preferentially retrained through sediments, and a higher K d value indicates greater liquidity. Both BSAF and BAF are employed to describe the capability of organisms accumulating heavy metals from other mediums in the surrounding. BSAF < 1 or BAF > 100 demonstrates a low accumulation effect between the surrounding environment and organisms (Liu et al., 2022).

Pollution Assessment
The standard index method (single water quality index (WQI) method) was used to assess the degree of pollution of each heavy metal in seawater, and the characteristic pollution factors were identified based on the ranking of the calculated results (Zhao et al., 2019). In order to assess the regional water quality condition, the WQI method was used for comprehensive water quality assessment (Xiao et al., 2021;. The calculation formula was as follows: where A i is the standard index of i pollutants; C i is the measured concentration of i pollutants; C si is the evaluation standard of i pollutants. The evaluation is carried out according to the corresponding category standards in the Seawater Quality Standard (GB3097-1997). WQI is a comprehensive index of water quality, A i is a single factor of the standard index, and n is the number of items of all the water quality parameters involved in the assessment. WQI < 1 means that the water quality is excellent; 1 < WQI < 2, good; 2 < WQI < 3, normal; 3 < WQI < 5, poor; and WQI > 5, very poor .

Biological Quality Assessment
In this study, primary productivity and diversity index were used to characterize the biological quality in the study sea, where primary productivity characterizes the richness of the marine environment and diversity index characterizes the dominance of marine organism populations. The primary productivity was estimated by the chlorophyll a method according to the simplified formula proposed by Cadeé and Hegeman (1974). The diversity index was calculated using the Shannon-Wiener index formula, as follows: where P is the primary productivity, mg·C/m 2 ·d; C a is the surface chlorophyll a concentration, mg/m 3 ; Q is the assimilation coefficient, mg·C/(mgChl-a·h), based on the results of the South China Sea survey, which was taken here as 3.70 . H′ is the species diversity index, S is the total number of species in the sample, and P i is the number of individuals and the total number of species in the ith sample.
Pi is the ratio of the number of individuals in the ith to the total number. Uniformity J is used to characterize the integrity of the marine ecosystem, and Pielou's uniformity formula is J = H ′/log 2 S, where J represents uniformity, H′ characterizes the species diversity index, and S is the total number of species in the sample.

Species Sensitivity Distribution Assessment
The SSD method was used to evaluate the ecological effects of heavy metal pollutants in seawater, which is an advanced statistical extrapolation method with high confidence based on the differences in the sensitivity of different species to pollutants, and has the advantages of conciseness and clear ecological significance. The Burr type III distribution was chosen to fit the exposure concentration accumulation curve, and the Burr type III distribution covers a large range in both skewness and kurtosis, which makes it more flexible to fit the parameters (Li et al., 2018b;Liu et al., 2020). The parametric equation of the Burr type III function F(x) is as follows: In the equations, b, c, and k are the three parameters of Burr type III distribution; HC5 is the cumulative concentration corresponding to the proportion of hazardous species on the SSD curve when it reaches 5%; and the smaller HC5 means that it represents a more toxic heavy metal. In equation (9), q is the corresponding protection level, and PAF is the proportion of hazardous species corresponding to each measured concentration on the SSD curve. In equation (10), x is the measured concentration of pollutants. In equation (11), msPAF is the proportion of hazardous species produced by the combination of multiple pollutants, which can reflect the joint pollution of multiple pollutants in the water body.

RESULTS AND DISCUSSION
3.1 Characteristics of Heavy Metals in Seawater, Sediment, and Marine Organism

Contamination Characteristics of Heavy Metals in the Seawater
The water quality evaluation results of the sea area calculated by using the factor evaluation method are shown in On the other hand, the water quality in the lower layer of seawater (B layer) is significantly better than that in the upper layer (A layer), mainly since the upper layer is the main area of influence of farming, industrial discharges, and other activities.
The main source of inorganic nitrogen and phosphate, as the main pollutants, is the marine agricultural farming activities near the sea area, and the release of large amounts of inorganic nitrogen and phosphate in this process is the main reason for the exceedance of these two pollutants in the overall sea area. The analysis results of a heavy metal evaluation in the coastal seawater of Zhuhai are shown in Figures 2, 3. From Figures 2, 3, it can be seen that the concentration of Cu, Zn, Cd, Cr, Hg, and As in the seawater in Autumn 2017 are in line with the first-class standard of seawater quality, and only one site has Pb concentration exceeding the second class standard of seawater quality. The concentration of Cu, Pb, Cd, Cr, Hg, and As in the seawater in spring 2018 is in line with the first-class standard of seawater quality, and only one site has a Zn concentration exceeding the second class standard of seawater quality. Zn concentration at one site exceeded the second standard for seawater quality. At the same time, the concentration of heavy metals in the seawater of the Gaolan Sea has a significant seasonal trend, and the concentrations of As, Cd, and Cr have significant seasonality. The concentrations of As, Cd, and Cr in autumn are higher than those in spring, and this result may be due to the higher rainfall in spring (569.6 mm) than in autumn (288.9 mm) (Cui et al., 2021) and the seasonal variation of river flow in Guangdong Province. This result suggests that dilution leads to lower concentrations of heavy metals in seawater; the seasonal distribution of precipitation in Guangdong is affected by the recession of cold high pressure. The rise and fall of summer winds leading to higher rainfall in spring were also reported in the studies of Zhang et al. (2017) and Liu et al. (2018). As, Cd, Cr, Cu, Hg, and Zn were detected at the highest concentrations at stations closer to the coastline, with a decreasing trend in the direction toward the ocean, while the highest concentrations of Pb were found at 4-5 km south of the coastline. The comparison of the distribution of heavy metal concentration in the sea area with other similar areas is shown in Table S1, which shows that the surface heavy metal hazard in the sea area of coastal Zhuhai is not very serious and is at a relatively acceptable level, but the proliferation of heavy metals should still be strictly prevented. The water quality evaluation results of Zhuhai coastal waters are shown in Figure 4. It can be found that in 2017, the water quality in the southwestern sea is worse than that in the northeast, while the two are opposite in 2018. This may be due to the seasonal industrial production and agricultural fishing activities. In 2017, there were more factory discharges along the southwestern coast, which led to the accumulation of heavy metals in seawater .

Contamination Characteristics of Heavy Metals in Sediments
The results of the heavy metals in sediments of coastal Zhuhai are shown in Table 1. It can be seen from Table 1 that the main exceedance factors of the marine protected area within the study area are Cu, Cr, and Cd, and the exceedance rates are 80%, 80%, and 40%, respectively. In contrast, the main exceedance factor in the industrial town discharge area is Cd, and the exceedance rate is 50%. The spatial distribution characteristics of heavy metals in the sediment of the Zhuhai coast were simulated by using the Kriging spatial interpolation method. From Figures 2, 3, it can be found that As in this study area is concentrated in the recreational tourism area in the southeast direction; Zn, Pb, Hg, and Cd are mainly concentrated in the discharge area of industrial towns in the east; and Cu and Cr are concentrated in the discharge area of industrial towns in the southwest direction. The highly concentrated sediment heavy metal distribution indicates that the main distribution of heavy metals in the sediment of this sea area is still influenced by the discharge of industrial towns under the influence of the oceanic monsoon climate. Although there will be oceanic monsoon currents driving the migration of heavy metals with seawater, in general, sedimentation activities are still dominant (Han et al., 2020). This requires the environmental and ecological management department to focus on the impact of industrial towns on the marine sediment environment in the later planning and management.

Contamination Characteristics of Heavy Metals in Marine Organism
The quality standards of shellfish in marine organisms refer to Marine Biological Quality (GB18421-2001), and the quality evaluation standards of other fish, crustaceans, mollusks, and other marine organisms adopt the standards in the Concise Regulations for the Comprehensive Survey of National Coastal Zone and Sea-coat Resources. The single index method was used to calculate the index of marine biological quality, and the results are shown in Table 1. The heavy metal Cu has the greatest impact on marine organisms, especially invertebrates and shellfish, followed by the heavy metal Pb. According to Pandey's study, organisms in the ocean near the coast are mostly affected by external factors (Pandey et al., 2021). Cu, Pb, Zn, and Cd in shellfish in the fall of 2017 survey area exceeded the first category of biological quality standards, and all other factors met the first category of biological quality standards. In addition, the Cd in the shrimp mayflies exceeded the standards stipulated in the Concise Regulations for the Comprehensive Survey of National Coastal Zone and Sea-coat Resources and the Technical Regulations for the Second National Marine Pollution Baseline Survey (Second Branch), while the quality offish, shrimp, and other marine organisms in other stations met the Concise Regulations for the Comprehensive Survey of National Coastal Zone and Sea-coat Resources and the Technical Regulations for the Second National Marine Pollution Baseline Survey (Second Branch). The quality of marine organisms in the survey area is good.
In the spring of 2018, there were no shellfish samples in the survey area, and only zinc in the nudibranchs exceeded the standard stipulated in the Concise Regulations for the Comprehensive Survey of National Coastal Zone and Sea-coat Resources, while the quality of other fish, shrimps, and other marine organisms all met the relevant index standards in the Concise Regulations for the Comprehensive Survey of National Coastal Zone and Sea-coat BSAF, biota-sediment accumulation factor; BAF, bioaccumulation factor. BSAF < 1 or BAF > 100 demonstrates low accumulation effect between surrounding environment and organisms (Liu et al., 2022).

Resources and the Technical Regulations for the Second National
Baseline Survey of Marine Pollution (the second volume), indicating that the quality of marine organisms in the survey area was good.

Marine Organism Quality Assessment
The results of the marine ecological characteristics survey in the sea area of coastal Zhuhai are shown in Figure 5. In terms of chlorophyll a and primary productivity, the variation of chlorophyll a concentration in the surface water column of the sea area surveyed in Autumn 2017 ranged from 0.74 to 5.66 mg/m 3 , with a mean value of 1.58 mg/m 3 , and the variation of primary productivity ranged from 5.99 to 75.62 mg·C/(m 2 ·d), with a mean value of 23.55 mg·C/(m 2 ·d). The variation of chlorophyll a concentration in the surface waters of the sea area in Spring 2018 ranged from 1.42 to 11.42 mg/m 3 with a mean value of 4.70 mg/m 3 , and the variation of primary productivity ranged from 25.97 to 464.00 mg·C/(m 2 ·d) with a mean value of 201.89 mg·C/(m 2 ·d). Overall, the productivity index in the northeast is higher than in other regions, indicating that this region is the priority production fishing area in the study area. The variation in the number of phytoplankton species at each station in Autumn 2017 ranged from 4 to 22 species, with an average of 12 species (Figure 6). The Shannon-Wiener diversity index ranged from 0.009 to 1.486, with an average of 0.366, and the distribution of biomass among species at each station was very heterogeneous. In contrast, the results of Spring 2018 survey showed that the number of phytoplankton species varied from 15 to 35 species at each station, with an average of 23 species. The Shannon-Wiener diversity index ranged from 1.836 to 2.487, with an average of 2.175, and Pielou's evenness index ranged from 0.601 to 0.878, with an average of 0.701. The distribution of biomass among species was relatively uniform, and the planar distribution trends of diversity index and evenness index were basically consistent.
Compared with plankton, benthic organisms showed less variation. The range of species variation of macro-benthic organisms at each quantitative sampling station in the Fall 2017 survey area was from 4 to 14 species/station, with an average of 9 species/station. The variation of the diversity index (H′) ranged from 0.21 to 2.39, with a mean value of 1.56 (Table S2), and the benthic biodiversity index in the surveyed sea area was moderately low. The range of uniformity was from 0.15 to 1.00, with an average of 0.72, reflecting the uneven distribution of species, and the variation in the number of macro-benthic species at each quantitative sampling station in the spring of 2018 ranged from 2 to 14 species/s tation, with an average of 9 species/station. The variation of diversity index (H′) ranged from 0.451 to 2.199, with a mean value of 1.421 (Table S2). The benthic biodiversity index of the surveyed sea area was moderately low. The evenness ranged from 0.268 to 0.945, with a mean value of 0.691, reflecting a relatively even distribution of species among the stations.

Interrelation and Source Analysis of Heavy Metals
The results of heavy metal and other particle correlation analyses in this study area are shown in Figure 7. The correlation matrix table shows that salinity was negatively correlated with As and Cr concentrations, indicating that salinity may hinder seawater mixing and maintain higher concentrations of As in surface seawater. Pearson's correlation analysis further showed that the concentration of Cu, Pb, and Zn in coastal sediments of the high bar sea area was significantly and positively correlated with the Gross Domestic Product (GDP) and wastewater discharge of each province (p = 0.01), suggesting that these three metals are more influenced by seasonal changes and fishing periods. Pb, Zn, Cd, and As were all significantly and positively correlated with seasons (p = 0.01), suggesting that the concentration of these four heavy metals in sediments is more influenced by local seasonal changes. Cu concentration was significantly and positively correlated with the gross industrial product (p = 0.005), indicating that Cu concentration was mainly influenced by the discharge of local industrial wastewater. K d is the specific coefficient employed to quantify the partition and transport of heavy metals between different phases such as soluble and solid. It can be affected by both the characteristics of heavy metals and the nature of sediment and seawater. The results of heavy metal enrichment in coastal Zhuhai are shown in Table 2. In this study, the order of value K d in both 2017 and 2018 was Pb > Cr > Cu > As > Zn > Hg > Cd, and the value of Pb was much larger than the value of other heavy metals, which suggested that the heavy metal Pb was most easily transported and enriched from seawater to the nearby sediment environment. From the temporal analysis, the overall K d values of each heavy metal in 2017 were smaller than the K d   Zhou et al., 2004 values obtained in 2018, which indicated that the lower seawater temperature was not conducive to the enrichment of heavy metals in the sediment environment. In terms of orientation, the distribution of K d values of heavy metals was clearly characterized. The K d values of points near the southwest and southeast coasts are much higher than those in the sea area far from the coast, which indicated that the more submerged the seawater depth was, the easier the heavy metals were enriched. From Table 3, we can see that the pollution of heavy metals in the sediments along the coast of Zhuhai generally showed a decreasing trend; the concentration of Hg, As, and Cr remained unchanged on the whole; and the concentration of Cu, Pb, Cd, and Zn decreased year by year except for the reasons of data errors in individual years. In general, due to the increased national attention to marine heavy metal pollution in recent years, marine heavy metal concentration may decrease in the next few years.
The BSAF values of each heavy metal except Pb in 2017 and 2018 were less than 1, indicating that only the heavy metal Pb is easily enriched in organisms from seawater or sediment environment. Other heavy metals are not easily enriched in organisms due to their metal properties. This required local managers to pay special attention to the management of heavy metal Pb to prevent elevated biological health risks. Geographically, Pb was more likely to be enriched in organisms in the inner bay of the northeastern part of the Zhuhai-Gaolan Peninsula compared to other locations, indicating that the mobility of seawater also had an influence on the bioconcentration of Pb. The highly mobile offshore area was less prone to Pb enrichment and was suitable for mariculture production activities in this area. The above results were also verified by BAF.
From the perspective of heavy metal classification, the main source of As was aquaculture in the sea area of coastal Zhuhai. In the southern sea area, aquaculture was the main source of the local economy, the large area of aquaculture caused the disorderly placement of aquaculture feed, and the discharged seawater will then cause the accumulation of As in the seawater flowing through the sea area . Pb, Zn, Cu, and other heavy metals were mainly from the local fine processing industries such as precision electronic instruments, chips, electronic product parts processing, and other industries. Coastal Zhuhai is a national industrial park established by the Chinese government in Zhuhai, equipped with a large number of precision electronic processing industries, which inevitably affect the surrounding seawater environment in the process of continuous development and growth, thus causing the accumulation of heavy metals in the sediment.

Ecological Risk Assessment
The potential biohazard effects of heavy metals in this sea area, calculated according to the species sensitivity analysis function, are shown in Figure 8 and Table S3. It can be seen that the biohazard coefficient of heavy metals in the study area was significantly higher in Autumn 2017 than in Spring 2018, which may be due to the weakened seawater activity and elevated heavy metal concentration due to low temperatures in autumn. From the analysis of heavy metal species, it was found that Cd, Cu, and Zn were the main heavy metals causing biohazard, while As, Hg, and Pb did not produce heavy metal biohazard. Taken together, the potential biohazard coefficients and ratios (more than 6%) in the northeastern sea area are significantly higher than those in other sea areas. Based on the heavy metal characterization of ocean dynamics, this was caused by seawater currents and coastal topography (Jahan and Strezov, 2018;Khalid et al., 2021). This requires the local ecological and environmental protection department to pay close attention to the protection of marine ecology in the northeast. In addition, the spatial distribution of bio-sensitivity of heavy metals in the sea area of coastal Zhuhai has obvious directionality. From Figure 9, it was found that the bio-sensitivity survey points were all located in the southeast and southwest areas of coastal Zhuhai, and these two areas were also the main precision instrument processing areas in coastal Zhuhai. A large number of metals processing several times the heavy metals discharged from the production of processing industries were accumulated in the sediment, which was the main reason for the significantly higher risk of heavy metal bio-sensitivity in these areas than in other areas. The biological species and biological quality in the sea area were greatly affected by the heavy metal concentration, and the number of biological populations and biological quality that suffer from the risk crisis can also seriously affect the ecological and environmental safety of the sea area in this region (Jawad et al., 2019;Han et al., 2020). So there is a need to carry out wastewater disposal for precision instrument processing enterprises in the southeast and southwest, which can cause the bio-sensitivity crisis collection to reduce the bioecological risk in the region.

Priority Control and Recommendations
In general terms, few previous studies showed spatial differences in heavy metal contamination on the Zhuhai coast, but while comparing the data with seasonal responses in this study, Cu, Pb, and Zn demonstrated higher values in winter in seawater and sediment, especially near the coast, and this might indicate that organisms near the coast suffered pro-oxidant heavy metal enrichment pressures during the cold season. This explained why organisms, especially invertebrates' biochemical responses, would be higher in winter, which is the beginning of the reserve accumulation period on the Zhuhai coast. From the perspective of aquaculture production, the survey of heavy metals distribution on the Zhuhai coast can further improve the quality of marine products and effectively evade the pollution areas. From the perspective of environmental management, the data of HM  distribution and risk characteristics can provide accurate information for precise heavy metal pollution control and avoid wasting government resources for ecological protection. Heavy metals are widespread in aquatic environments on the Zhuhai coast, and understanding their distribution at the local scare is necessary to advise marine environmental protection policymakers. The key issue for heavy metal contamination on the Zhuhai coast is the scarcity of emission control, and it is linked to coastal areas of industrial production, which do not have any option other than polluted rivers as sources for sewage discharge. Industrial wastewater treatment needs massive expenditure and reduces local economic revenue due to the installation of costly treatment installations. For better safeguarding and management of the marine environment including organisms, it is paramount that regional heavy metal polluting enterprises implement the most stringent wastewater treatment and discharge control measures in water bodies. Consequently, provisions and specialized supervision departments and supervisors should be introduced, consisting of GIS online monitoring in order to tap into local contamination information, as well as to recognize pollution sources and information requirements.
Our research exhibited that the cumulative effects of different heavy metals in seawater and sediment were different, indicating that a single evaluation method was not useful to evaluate heavy metal contamination on the Zhuhai coast. The heavy metal distribution characteristics and ecological risk curve evaluation system established in this study can effectively assess the cumulative effects of heavy metals in seawater, sediment, and organisms and can provide scientific and effective information support for the next phase of ecological restoration work. In addition, advance planning is needed based on Zhuhai's special geographical environment and industrial layout. This requires the local government to plan in advance the special disposal institutions and control areas for heavy metal wastewater in order to centralize the management and purification of heavy metal wastewater within the industrial park and to achieve zero pollution discharge. Therefore, the governmental environmental protection department, enterprise regulatory agencies, and local resident representatives need to participate in the development of the regional heavy metal control plan.

CONCLUSIONS
In the present study, multi-evaluation methods have been employed to investigate the heavy metal pollution on the Zhuhai coast. The results showed that the study area was heavily polluted by heavy metals, especially on the north coast of Gaolan Island. Comparative studies revealed that the main source of heavy metals was the discharge of local industrial wastewater, especially in offshore waters. The SSD evaluation of organisms showed that the northern part of the Zhuhai coast was in medium risk conditions, with high pollution levels of Cd and Cu.
Overall risk distribution patterns of heavy metals and zonation of SSD obtained from this study could support policymakers to take coherent action based on pollution control measures and marine aquaculture. Our research suggested that chemical fractionation should be considered to provide a more accurate assessment of the ecological risks of heavy metals in marine ecosystems. In addition, marine resources in the northern part of the Zhuhai coast may be more vulnerable to anthropogenic pollution. Our research results provided important information on the ecological response of organisms experiencing heavy metal contamination on the Zhuhai coast. Future studies should focus on integrated prevention and control strategies to ensure better conservation of the coastal and estuarine ecosystem.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
JW and SR developed the methodology. JW, SR, and MW wrote the original manuscript. JW and JL revised the manuscript. RL is responsible for contacting the editorial department and