Spatiotemporal Distribution Characteristics and Influencing Factors of Dissolved Potentially Toxic Elements along Guangdong Coastal Water, South China

: In this study, five potentially toxic elements (PTEs) (Hg, Pb, As, Zn, and Cu) and five physicochemical factors (pH, dissolved oxygen, salinity, chlorophyll a, and chemical oxygen demand) relating to surface seawater were measured along the coast of Guangdong Province (GD), China, during three seasons in 2019. Comprehensive analyses were conducted to explore the distribution characteristics, pollution degrees, and influencing factors of PTE. The results showed that the concentrations of PTE varied spatially and seasonally due to these elements’ complex sources and behaviors. Pb was the most abundant toxic element, reaching considerable contamination levels. Overall, the average pollution degrees of Cu, As, and Zn were higher around the east coast of GD, while Hg and Pb levels were higher in the west coast region of GD. The correlation analysis showed that the variation of different physicochemical parameters had different degrees of influence on PTE transport and transformation. This study can help environmental managers gain deeper insight into the influence of complex factors on PTE and improve the efficiency of pollution control in this significant subtropical coastal area.


Introduction
Coastal ecosystems are among the most important ecosystems on earth, providing a wide range of ecological services and being of high economic value [1][2][3].However, potentially toxic elements (PTEs) have accumulated rapidly in marine environments in recent decades due to intensive human activities conducted in coastal areas and upstream of estuaries [4][5][6].PTE pose great potential threats to aquatic ecosystems because of their wide distribution, persistence, non-degradability, and toxicity [7][8][9].Increasing pollution loads of PTE not only threaten the habitats of aquatic organisms but also affect human health and the structure and function of the entire ecosystem through accumulation and amplification in the food chain [10][11][12].In estuaries and coasts, the interaction of tidal currents and runoff results in the exchange of freshwater from rivers with seawater of higher salinity.Dramatic changes in various physical elements (water flow, velocity, suspended particulate matter, etc.) and chemical elements (salinity, macronutrient ions, dissolved oxygen (DO), dissolved organic carbon, pH, etc.) in estuarine waters, as well as various complex biogeochemical interactions in these waters, affect the transport, migration, flocculation, and settling of floating particulate matter and its biological processes, resulting in the transport and transformation of PTE in and into various forms [13,14].These processes affect the biogeochemical behavior of PTE distribution, transport, and accumulation in estuarine and coastal environments and have a direct impact on the self-purification of water bodies, the with a mean annual temperature of 22 °C and a mean annual precipitation of 1800 mm [40].Rain and heat coincide in the same season, with precipitation mainly concentrated in April to September.Along the coastline from south to north, there are 14 cities in GD Province, namely ZhanJiang (ZJ), MaoMing (MM), YangJiang (YJ), JiangMen (JM), ZhuHai (ZH), ZhongShan (ZS), GuangZhou (GZ), DongGuan (DG), ShenZhen (SZ), Hui-Zhou (HZ), ShanWei (SW), JieYang (JY), ShanTou (ST), and ChaoZhou (CZ) City, as well as two Special Administrative Regions: Hong Kong and Macao.These 14 cities have distinct development patterns and account for more than 80% of GD's GDP.
There are two major cross-provincial rivers along the coast that flow into the South China Sea.One is the Pearl River, a complex water system formed by the convergence of the Xijiang, Beijiang, Dongjiang and the rivers of the Pearl River Delta.It is the third-largest river in China and the largest water system flowing into the South China Sea.It has a total length of 2320 km, a basin area of 453,690 km 2 , and an annual runoff of more than 330 × 10 7 m 3 [42].The other is the Han River, with a main stream length of 470 km, a total catchment area of 30,112 km 2 , and a total annual runoff of 24.5 × 10 7 m 3 .

Sampling and Analysis
Three campaigns were carried out in April (spring), August (summer), and October (autumn) of 2019 along the coast of GD Province.Surface water samples were collected from 242 stations.A spatial distribution map of the sampling points is provided in Figure 1.In this work, the study area was mainly divided into three parts: the west zone, central zone, and east zone.The seawater samples were collected on board from 0.5 m below the water surface at each site using 5 L Niskin bottles and immediately filtered through a 0.45 µm cellulose acetate membrane.Additionally, the water samples were acidified to pH < 2 using nitric acid for Cu, Zn, and Pb and sulfuric acid for As and Hg.The filters and water samples were stored at −20 °C.The salinity, pH, and DO were measured onboard with respective meters.
Anodic stripping voltammetry (ASV) (797 VA, Metrohm, Herisau, Switzerland) was performed for Cu, Pb, and Zn analysis.ASV is the most commonly used technique known The seawater samples were collected on board from 0.5 m below the water surface at each site using 5 L Niskin bottles and immediately filtered through a 0.45 µm cellulose acetate membrane.Additionally, the water samples were acidified to pH < 2 using nitric acid for Cu, Zn, and Pb and sulfuric acid for As and Hg.The filters and water samples were stored at −20 • C. The salinity, pH, and DO were measured onboard with respective meters.
Anodic stripping voltammetry (ASV) (797 VA, Metrohm, Herisau, Switzerland) was performed for Cu, Pb, and Zn analysis.ASV is the most commonly used technique known as the stripping electroanalytical method, which involves reconcentrating an analyte at the electrode surface in order to lower the detection limit for that analyte.The characteristic peak voltages used for Cu, Pb, and Zn were −0.30V, −0.52 V, and −1.1 V, respectively.Atomic fluorescence spectrometry (AFS) (AFS-9562, Beijing Haiguang Instrument, Beijing, China) was performed to assess As and Hg.The samples were converted to gaseous atoms, and the element of interest was excited to a higher electronic energy level via a light source.
Following excitation, the atoms were deactivated by the emission of a photon.The measured fluorescence stems from this emission process.The AFS negative high voltage and lamp current were 265 V and 28 mA for Hg detection and 275 V and 60 mA for As detection.To ensure the precision and accuracy of the analysis's results, blanks, duplicate samples, and sample spikes were analyzed.The seawater-certified reference materials GBW(E)080040 (Second Institute of Oceanography, MNR) for Cu, Pb, and Zn; GBW(E)080042 (Second Institute of Oceanography, MNR) for Hg; and GBW(E)080117 (National Institute of Metrology, Beijing, China) for As were used for analysis validation.The relative standard deviation of duplicate samples was <11%, and the recovery rates of all the PTE varied from 88% to 105%.Chl a and COD were determined using spectrophotometry and alkalescent permanganate titration methods.The determination limits were 0.007 µg/L, 0.3 µg/L, 0.5 µg/L, 0.6 µg/L, 1.2 µg/L, and 0.11 µg/L for Hg, Pb, As, Cu, Zn, and Chl a, respectively.All the processes involved in pretreatment, preservation, and analysis were conducted according to the Code of Practice for Marine Monitoring Technology (GB 17378-2007) [43].

Pollution Assessment Methods
The Nemerow index has been widely used to evaluate environmental quality [44,45].In this method, first, the signal factor pollution index (P i ) is calculated; it is the ratio between the measured concentration of each PTE and its background concentration and indicates the individual pollution degrees of different PTEs.Then, the average and maximum values of P i were obtained to calculate the comprehensive water quality index (CWQI), which was used to reflect the total pollution degrees of 5 PTE at each site.
where c i and c s are the measured concentrations and corresponding background values of each PTE, respectively; P i is the signal-factor pollution index, which represents the exceeding multiple; CWQI is the comprehensive water quality index; and P i and P i max represent the average and maximum P i values of the 5 PTE, respectively.The grade-one seawater quality standard (GB 3097-1997) (Table S1) was used as the background value of the PTE [46].Thus, the c s values of Hg, Pb, As, Cu, and Zn are 0.05, 1, 20, 5, and 20 µg/L, respectively.Coastal water pollution could be divided into 4 classes according to P i : low pollution (P i < 1), moderate pollution (1 < P i < 3), considerable pollution (3 < P i < 6), and very high pollution (P i > 6).Based on the CWQI value, the water pollution situation could be divided into four levels: uncontaminated water (CWQI ≤ 1), minor contaminated water (1 < CWQI ≤ 2), moderately contaminated water (2 < CWQI ≤ 3), and heavily contaminated water (CWQI > 3).

Data Analysis
The spatial distribution of PTE in coastal water is directly related to the intensity of human activities along the coast and influenced by the geological and physical characteristics of seawater.The Kruskal-Wallis test, a non-parametric method, was used to explore the statistical differences among seasons and among areas.Spearman correlation analysis, principal component analysis (PCA), and two hierarchical cluster analyses were used to investigate the possible sources and explore the correlations among PTE and influencing factors.The normality and equal variance of the data were tested before further processing.The statistical analysis was conducted with the "SciPy" module in Python (Version: 3.9.12).Box plot was conducted with "seaborn" module in Python.

Spatial and Seasonal Distributions of Physicochemical Factors
Figure 2 exhibits the seasonal and regional variations in physicochemical factors.Figure S1 illustrates the spatial distributions of physicochemical factors over three seasons.Table S2 lists the p-value set for the Kruskal-Wallis test.The statistical data of physicochemical factors are shown in Table S3.COD is a comprehensive index reflecting the organic pollution status of water bodies, and it is more related to human activities than other physicochemical factors.The central area had the highest COD mean value (1.55 mg/L) compared with that of the west area (0.72 mg/L) and the east area (1.02 mg/L) in different seasons.In addition, the sites with higher COD values were mainly concentrated in the inner Pearl River estuary and Daya Bay, which are part of the central area (Figure S1).This phenomenon may be due to the heavy organic pollutants discharged from the large upstream watershed.However, this pollution had a limited influence on the offshore area due to the dilution and diffusion effects.In the central zone, COD corresponded to the following ranking: spring > summer > autumn; this result may be related to the continuous thunderstorms in the spring and drought in the autumn in 2019 [49].Furthermore, it was found that the west and east zones had relatively high COD levels in the spring and summer, respectively.This may be highly affected by the seasonal patterns of the coastal ocean current near GD.The currents flow northeastward in the summer and southwestward in the spring and winter [50,51].To an extent, the flow direction of a coastal current can influence the spread of the fresh water sourced from a river and the extent of pollution [34].The highest mean Chl a level was 6.93 µg/L, which appeared in the central area during the summer.This may be due to the abundant nutrients from the upper Pearl River and higher temperatures, which promote phytoplankton growth in the summer.The variation range between different seasons in According to the results of the Kruskal-Wallis test, the differences among the three areas and among the three seasons were statistically significant for each physicochemical factor.The highest seasonal mean salinity was found in the west area (29.67), followed by the east area (27.49) and then the central zone (22.38).This difference was statistically significant.The lowest salinity occurred in the central area in the summer (20.78), and this may be due to the fact that the Pearl River has the greatest amount of runoff in the summer.The salinity of estuary waters in the central zone increased gradually along the river in the longitudinal direction, and the corresponding variance range is relatively high (from <0.1 to >34).Thus, salinity partly reflects the distance of the sampling point from the estuary in the central zone to some extent [47].DO is an important parameter for evaluating seawater quality, which is influenced by temperature, biological and non-biological oxygen consumption [48], etc.In this study, DO showed a slight variation among the different seasons and areas.Several sampling sites in the narrow inner Pearl River Estuary had lower DO values (lower than 5), especially in the summer; this phenomenon can be attributed to the greater amount of dissolved oxygen being consumed by nutrient degradation in the Pearl River Estuary.Furthermore, the average dissolved oxygen content was lower than that in the other two seasons because of the high temperatures in the summer.Similarly, the central area had lower values and a higher variation in pH, with the lowest value occurring in the autumn (7.8).Higher mean pH levels were found in the east area, with the highest value occurring in the spring (8.15).The pH values of most sites range between 8 and 8.5.
COD is a comprehensive index reflecting the organic pollution status of water bodies, and it is more related to human activities than other physicochemical factors.The central area had the highest COD mean value (1.55 mg/L) compared with that of the west area (0.72 mg/L) and the east area (1.02 mg/L) in different seasons.In addition, the sites with higher COD values were mainly concentrated in the inner Pearl River estuary and Daya Bay, which are part of the central area (Figure S1).This phenomenon may be due to the heavy organic pollutants discharged from the large upstream watershed.However, this pollution had a limited influence on the offshore area due to the dilution and diffusion effects.In the central zone, COD corresponded to the following ranking: spring > summer > autumn; this result may be related to the continuous thunderstorms in the spring and drought in the autumn in 2019 [49].Furthermore, it was found that the west and east zones had relatively high COD levels in the spring and summer, respectively.This may be highly affected by the seasonal patterns of the coastal ocean current near GD.The currents flow northeastward in the summer and southwestward in the spring and winter [50,51].To an extent, the flow direction of a coastal current can influence the spread of the fresh water sourced from a river and the extent of pollution [34].The highest mean Chl a level was 6.93 µg/L, which appeared in the central area during the summer.This may be due to the abundant nutrients from the upper Pearl River and higher temperatures, which promote phytoplankton growth in the summer.The variation range between different seasons in the central zone (from 0.2 to 66.6 µg/L) was more apparent than that in the other two areas.Additionally, the east area had the lowest Chl a (1.7 µg/L) compared to that of the central (5.36 µg/L) and west areas (3.2 µg/L).

Spatial and Seasonal Distributions of PTE
Figure 2 exhibits the seasonal and regional variations in dissolved PTE.Table S2 lists the p-value set for the Kruskal-Wallis test, and Figure S2 illustrates the spatial distributions of PTE in the three seasons.The seasonal and regional statistical data on PTE are listed in Table S3.
According to the results of the Kruskal-Wallis test, the differences between the three areas and the three seasons were statistically significant for each PET.The spatial distribution of Hg exhibited a similar pattern throughout the seasons, with higher Hg values in the Leizhou peninsula, Daya Bay, and the east coast.The west zone had the highest mean value of Hg (0.024 µg/L), followed by the east zone (0.019 µg/L) and central zone (0.011 µg/L).In the spring and summer, the west zone had the highest Pb concentrations (1.33 µg/L and 1.08 µg/L), especially the area near YangJiang city and Leizhou peninsula, respectively.However, there was no clear pattern for Pb in the autumn.Compared with Hg and Pb, Cu had the highest values in the east zone (1.39 µg/L), followed by the central (0.96 µg/L) and west zones (0.61 µg/L).Relatively high Cu levels were present in the inner Pearl River Estuary in the summer, suggesting that the higher runoff of the Pearl River in the summer brought more Cu to the sea.Although the differences in the Zn mean value in each region were not large, they had obvious spatial distribution characteristics in the three seasons studied.Zn concentrations in the central region were low in the spring, but these values were higher in the summer and winter, and some extreme values were detected in the center zone during the summer.As levels were more evenly distributed spatially, with slightly higher concentrations in the east than in other regions in the autumn.The highest As concentration, 2.46 µg/L, was found in the east in the summer.In general, the average concentrations of Hg, Pb, Cu, Zn, and As in the whole study area were 0.018, 0.80, 0.99, 4.54, and 1.67 µg/L, corresponding to the following order: Zn > As > Cu > Pb > Hg.
Compared with the seawater quality standard (GB3097-1997) [46], nearly all the PTE measurement results did not exceed the Grade-II water quality standard, and most of them met the Grade-I water quality standard.For Hg, Pb, Cu, and Zn, about 94%, 77%, 99%, and 98% of all sampling sites, respectively, had concentrations lower than the Grade-I limit.Only one Pb measurement result was greater than the Grade-II limit.For As, the concentrations from all the sampling sites were lower than 20 µg/L, which is the Grade-I limit.
Figure 3 shows the P i values of five PTEs, and detailed data are shown in Table S4.This heavy pollution was caused by Pb (max P i :5.1), which reached considerable contamination levels.Although Zn was the most abundant element, its pollution degree reached moderate contamination levels like those of Hg and Cu.Hg had more sampling results classified as moderate contamination than Zn and Cu.In addition, since the P i values for As were significantly lower than 1, As falls into the low-contamination range.Furthermore, the average P i values of Cu, As, and Zn in the east area were generally higher, while those for Hg and Pb were higher in the west area; this may be due to the disparity of the industrial distribution characteristics of GD.The west area is rich in mineral resources, and the petrochemical and metal products industries are the major industries in this region.However, the most widely distributed industries in the eastern part of GD are the textile and electrical and electronic industries [29,52].The average and seasonal spatial distributions of CWQI are shown in Figure S3.It can be observed that the majority of sites belonged to uncontaminated water (CWQI ≤ 1).Minor contaminated water was detected in the vicinity of Leizhou Peninsula during the spring and summer seasons (1 < CWQI ≤ 2).Moderately contaminated water was found in the west area during spring (2 < CWQI ≤ 3).And there was no instance of heavily contaminated water.Overall, the level of comprehensive PTE pollution in the coastal area of western GD was significantly higher than that in the other areas.
Furthermore, the current PTE data were compared with previously measured data and recent sampling results obtained from adjacent and global sea water bodies (Table 1).It is notable that the levels of Cu, Pb, Zn, and As were lower than the values observed 13 years ago, especially regarding Zn, and this phenomenon is probably related to the continuous strengthening of environmental protection in recent years [41].Similarly, higher Cu, Zn, and As levels were detected in east GD, indicating that eastern GD has been more polluted by PTE than the west.In addition, the sorted order of the metals showed an analogous pattern: Zn > As > Cu > Pb > Hg.Compared with recent studies on GD's coastal areas, the concentrations of PTE in this study were in the range of the values in the seawater along the GD coast.In Honghai Bay, the concentrations of Zn and Hg in 2017 were approximately twice the mean values observed in the eastern area of this study [53].This could be attributed to the proximity of sampling sites to the coast in 2017, indicating a significant influence of land-based pollution on Zn and Hg levels.However, the concentration of As in this study was higher than that observed in 2017.In Daya Bay, Hg in 2020 was significantly higher than the mean value of the central area of this study [39].This could be attributed to the fact that sampling sites were located within the bay in 2020 and there was an evident positive correlation between Hg and oil [39].Considering the local economic composition, the petrochemical sewage, ship transportation and oil terminal contamination are likely to be the main sources of Hg in Daya Bay.In the Pearl River Estuary, previous studies have reported significantly higher levels of Zn compared to our findings [40].The fluctuation in PTE results may be influenced by the sampling time, area, and statistical methods [39,40,53].Compared with the reports from North China, the average levels of Cu, Pb, and Zn in GD were lower than those in Bohai Bay, the Shandong Peninsula, and the Yellow River Estuary [54,55].The Luanhekou Estuary has lower Cu and Pb levels and higher Zn levels [54].The lowest Cu, Pb, and Zn levels were found in the Yangtze River Estuary [56].Compared with Beibu Gulf to the west of GD, the concentrations of Cu, Hg, and Zn were lower in GD, while the levels of As were higher [57].Overall, the concentration levels of Cu, Pb, and Zn in the GD were at a low level compared to the results from other coastal waters studies conducted in China.Hg and As were present in moderate levels according to the listed concentrations.In addition, the average concentrations of PTE were lower than those in Tampamachoco in the Gulf of Mexico, the Gulf of Suez in Egypt, and the Meghna Estuary in Bangladesh, except with regard to Cu in Egypt [58][59][60].Palk Bay in India had the lowest levels of Cu, Pb, and Zn [61].comprehensive PTE pollution in the coastal area of western GD was significantly higher than that in the other areas.Furthermore, the current PTE data were compared with previously measured data and recent sampling results obtained from adjacent and global sea water bodies (Table 1).It is notable that the levels of Cu, Pb, Zn, and As were lower than the values observed 13 years ago, especially regarding Zn, and this phenomenon is probably related to the continuous strengthening of environmental protection in recent years [41].Similarly, higher Cu, Zn, and As levels were detected in east GD, indicating that eastern GD has been more polluted by PTE than the west.In addition, the sorted order of the metals showed an analogous pattern: Zn > As > Cu > Pb > Hg.Compared with recent studies on GD's coastal areas, the concentrations of PTE in this study were in the range of the values in the seawater along the GD coast.In Honghai Bay, the concentrations of Zn and Hg in 2017 were approximately twice the mean values observed in the eastern area of this study [53].This could be attributed to the proximity of sampling sites to the coast in 2017, indicating a significant influence of land-based pollution on Zn and Hg levels.However, the concentration of As in this study was higher than that observed in 2017.In Daya Bay, Hg in 2020 was significantly higher than the mean value of the central area of this study [39].This could be attributed to the fact that sampling sites were located within the bay in 2020 and there was an evident positive correlation between Hg and oil [39].Considering the local economic composition, the petrochemical sewage, ship transportation and oil terminal contamination are likely to be the main sources of Hg in Daya Bay.In the Pearl River Estuary, previous studies have reported significantly higher levels of Zn compared to our findings [40].The fluctuation in PTE results may be influenced by the sampling time, area, and statistical methods [39,40,53].Compared with the reports from North China, the

Statistical Analysis Results of Physicochemical Factors and PTE
Figure 4 presents the correlation coefficient and significance results between physicochemical factors and PTE concentrations.The correlations among physicochemical factors were similar among different seasons and areas.Salinity exhibited significant negative correlations with COD and Chl a in different seasons and areas but positive correlations with pH except in the west area.COD was positively correlated with Chl a. DO had weakly positive correlations with Chl a and pH and no significant correlations with salinity and COD.
Salinity, as a conservative factor, could reflect the behavior of PTE in the estuarine dilution-mixing process.Salinity was positively correlated with Hg, except in the autumn and in the west zone, and with Pb in the summer and autumn in the west zone, indicating that Hg and Pb were more influenced by the adsorption of suspended particles and atmospheric deposition.Leaded gasoline has been considered the most significant source of Pb pollution [62].With the abolition of leaded gasoline usage, the major sources of Pb in the environment shifted to coal combustion and non-ferrous metal metallurgy.High amounts of salt water can inhibit the removal of Hg from seawater, enhancing the desorption effect of suspended particles, and resulting in the transformation of Hg from suspended particles to a dissolved state [63].In addition, salinity had negative correlations with Cu and As in the central and east areas, indicating that the sources of Cu and As in this region were terrigenous.
COD is highly influenced by terrestrial pollution and discharged into the sea through rivers, and it can be used to indicate the contribution of upstream pollution.There were no significant correlations between COD and Hg, and the correlations between COD and Pb were negative in the east and in the summer, suggesting that terrestrial pollution might not be the major source of Hg and Pb.On the contrary, COD was positively correlated with Cu in the central and east areas and with As in the central zone.
Changes in the pH of seawater will affect the adsorption and desorption of PTE in suspended particulate matter.In general, the amount of PTE adsorbed by suspended particles in seawater increases with pH; conversely, hydrolytic desorption plays a dominant role, and the concentration of PTE in seawater increases [64].PH was weakly negatively correlated with Cu except in the autumn, with As in the central and eastern regions in the summer, and with Zn in the west and in the autumn.The pH effects on Hg and Pb were temporally and spatially different.
The concentration of DO in water can affect this medium's redox potential, thus affecting the migration and transformation of PTE in water.In general, PTE migrates more easily under aerobic conditions than anaerobic conditions [65].DO has positive correlations with Hg and Pb in the summer and with As in the east area, while it has varied correlations with Cu.Zn has no correlations with Do.In general, the effect of redox on PTE was limited.
Chl a is an indicator of phytoplankton biomass in the ocean.Phytoplankton have different reactions to the physicochemical properties of different PTEs.Zn and Cu are essential elements for life and are involved in biophysiological metabolic processes [24].As, as a trace element of the same main group as phosphorus, is absorbed by organisms as a phosphate analog in sea areas with low phosphorus content [66].Furthermore, Pb and Hg are highly biotoxic.Chl a levels were weakly negatively correlated with Hg in the central area in the summer, with Pb in the east area in the summer and autumn, and with Zn in the summer, indicating the biological effects on the distribution of these PTE.On the contrary, Chl a showed positive correlations with Cu and As in the east area.Table 2 lists the results of the PCA for physicochemical factors and PTE.The Kaiser-Meyer-Olkin (KMO) test was used to ensure the applicability of PCA analysis.The KMO statistics value was 0.668 and the significance probabilities of Bartlett's spherical test were less than 0.05, verifying the adequacy of PCA analysis.The analysis yielded three principal components (PCs) that accounted for a total of 65.63% variance.PC1, PC2, and PC3 explained 35.95%, 16.98%, and 12.70% of the total variance, respectively.PC1 exhibited positive loadings for salinity, pH, Do, and Hg and negative loadings for COD and Chl a, indicating the influence of freshwater discharge.Cu, As, and Zn mainly contributed to PC2, and Pb mainly contributed to PC3.The results of clustering analysis are illustrated in Figure 5 as a dendrogram.The ten items could be categorized into three clusters: (1) COD and Chl a, (2) As, Cu, and Zn, and (3) Pb, Hg, DO, salinity, and pH.The cluster results were in line with Spearman correlation and PCA analysis to some extent.As, Cu, and Zn may display similar origins.Hg and Pb showed different source patterns compared to As, Cu, and Zn.Hg exhibited a stronger association with terrestrial sources as it clustered together with DO, salinity, and pH.The origin characteristics of Pb differed from those of other elements.Luo found that atmospheric deposition is the main source of Pb in the Pearl River Estuary because of its volatile and ubiquitous character in the atmosphere.The fluxes of As and Zn from submarine groundwater discharge accounted for more than 54% of the total flux [40].In general, the specific distribution of PTE in different seasons is affected by complex factors.Variations in physicochemical factors tend to influence the transportation, transformation, etc., of PTE [67].However, the correlations between them were generally weak.The external inputs, including terrestrial sources from river discharge, atmospheric deposition from industries and fuel burning, and sedimentary sources, play a more important role.

Conclusions
In this study, the concentrations of five physicochemical factors and five PTE in the coastal water of GD were collected and analyzed in three seasons in 2019.Physicochemical factors showed obvious distribution characteristics that were highly affected by the Pearl River water.Salinity, DO, and pH levels were lower in the central area, while COD and Chl a levels were higher in that area.DO was lower in the summer, whereas COD and Chl a were higher in the summer.In general, the PTE concentrations showed the following decreasing ranking: Zn > As > Cu > Pb > Hg.However, Pb was the most abundant PTE, reaching considerable contamination levels.The concentrations of all measured PTE, except for Pb, were within Grade-II standards.Concentrations of As, Cu, Zn, Hg, and Pb were below the Grade-I limit at 100%, 99%, 98%, 94%, and 77% of all sampling points, respectively.The average pollution indices of Cu, As, and Zn in the east area were generally higher, while those for Hg and Pb were higher in the west area.The majority of sites belonged to uncontaminated water (CWQI ≤ 1).Minor and moderate contaminated water were detected in the west area.
Overall, two hierarchical clusters and PCA analysis revealed the potential source characteristics of PTE: Hg exhibited an association with terrestrial sources as it was grouped in the same cluster and PC along with DO, salinity, and pH.Pb appeared to be more influenced by the adsorption of suspended particles and atmospheric deposition.As, Cu, and Zn may display similar origins.It should be noted that the sources may vary across different locations and seasons.The correlation analysis showed that changes in physicochemical factors may influence the transportation, transformation, etc., of PTE.Although the level of PTE pollution in this region is slightly lower than it was before, the sources of PTE pollution are complex, and there is still a need to strengthen the control and management of PTE pollution.

Figure 1 .
Figure 1.Location of study area and sampling sites.

Figure 1 .
Figure 1.Location of study area and sampling sites.

J 16 Figure 2 .
Figure 2. Box plots of concentrations of physicochemical factors and PTE (The box extends from the first quartile to the third quartile of the data, with a line at the median.The whiskers extend from the box to the farthest data point lying within 1.5× the inter-quartile range from the box.The values outside of the whiskers can be treated as outliers).

Figure 2 .
Figure 2. Box plots of concentrations of physicochemical factors and PTE (The box extends from the first quartile to the third quartile of the data, with a line at the median.The whiskers extend from the box to the farthest data point lying within 1.5× the inter-quartile range from the box.The values outside of the whiskers can be treated as outliers).

Figure 3 .
Figure 3. Box plots of PTE pollution indices.

Figure 3 .
Figure 3. Box plots of PTE pollution indices.

Figure 5 .
Figure 5. Dual hierarchical cluster results of physicochemical parameters and PTE.Figure 5. Dual hierarchical cluster results of physicochemical parameters and PTE.

Figure 5 .
Figure 5. Dual hierarchical cluster results of physicochemical parameters and PTE.Figure 5. Dual hierarchical cluster results of physicochemical parameters and PTE.

Table 1 .
Concentrations of dissolved heavy metals s in Guangdong coastal surface seawater and comparison with other regions (µg/L).

Table 2 .
Principal component analysis results of physicochemical factors and PTE.