Biological Risk Assessment of Heavy Metals in Sediments and Health Risk Assessment in Marine Organisms from Daya Bay, China

: The concentrations of heavy metals in sediments and marine organisms in Daya Bay were investigated, and the Monte Carlo method was used to analyze the uncertainty of the results of geo-accumulation characteristics and ecological and health risks. The mean concentrations of metal elements in sediments were in the following order: Zn > Cr > Cu > As > Cd > Hg, while those in marine organisms were Zn > Cu > As > Cr ≈ Cd > Hg. The geo-accumulation index (I geo ) indicated that the primary pollutant was Hg, with 5.46% moderately polluted, and 39.52% for unpolluted to moderately polluted. Potential ecological risks (RI) were between low and high risks, and the contributions of Hg, Cd, and As to ecological risks were 50.85%, 33.92%, and 11.47%, respectively. The total hazard coefﬁcients (THQ) were less than 1, but on the basis of total carcinogenic risks (TCR), the probability of children and adults exceeded the unacceptable risk threshold of 22.27% and 11.19%, respectively. Sensitivity analysis results showed that the concentrations of carcinogenic elements contributed to risk in the order of As > Cd > Cr. Therefore, in order to effectively control heavy metals contamination in Daya Bay, it is necessary to strengthen the management of Hg, Cd, and As emissions.


Introduction
The impact of heavy metal pollution on the sustainability of biogeochemical cycles and ecological risks has received increasing global attention [1]. Heavy metals have special eco-toxicological effects, such as toxicity, accumulation, and refractory degradation [2]. When heavy metals enter the aquatic environment, they eventually accumulate in sediments, so sediments are usually used as a main medium for evaluating heavy metal pollution [3]. However, when the environmental conditions change, these are released from the sediments, and then used as a source to return to the overlying water, threatening the aquatic environments and marine organisms [2,4,5]. Marine organisms are rich in proteins, low saturated fatty contents and Omega-3 fatty acids, and other nutrients needed by the human body [6]. However, due to the enrichment and amplification of heavy metals in marine organisms, seafood that originally provided the human body with abundant protein may become a carrier of concentrated poisons, thereby endangering human health [7][8][9][10][11]. Therefore, identifying the concentration distribution of heavy metals in sediments and marine organisms and assessing the ecological and health risks under heavy metal exposure are of great significance to the protection of biological habitats and human safety.

Study Area and Samples Collection
The study area is located in the southwestern region of Daya Bay in Guangdong Province, South China. In January and March 2020, 27 surface sediment samples and 20 species of marine organism samples were collected, respectively ( Figure 1). The sediment samples were collected by precleaned stainless steel grab sampler, and two sediment samples were collected at each site to prepare a mixed sample. Marine organisms were collected by horizontal trawl and intertidal surveys, and selected adult individuals as samples. The details of organism samples are shown in Table S1. All samples collected from Daya Bay were placed in polyethylene sealed bags and stored at −4 • C until processing and analysis.

Chemical Analysis and Quality Assurance
The sediment samples were freeze-dried to a constant weight, ground with an agate mortar and pestle, passed through a 63 µm sieve, and stored at −20 • C in brown glass bottles. Heavy metal microwave digestion was conducted in accordance with the EPA Method 3050B using an Ethos Plus Microwave Labstation (Milestone Inc., Italy).
The marine organism samples were first thawed naturally at room temperature, then rinsed with deionized water for 3-5 times. The edible tissues were removed from each organism with stainless steel scissors and a scalpel, tissues form organisms of the same species were mixed together thoroughly, crushed with a stainless steel mixer, freeze-dried to constant weight, and ground to give a homogeneous powder, which was stored in a brown glass bottle at −20 • C. An aliquot of each dried and homogenized sample was digested using a microwave digestion system. The metal extraction procedure is described in detail elsewhere [1].
The quality control and assurance of sediment samples were carried out by using the standard recovery method. The recoveries of the metals were between 83% and 111%. The Chinese National Standard material (offshore marine sediment, GBW 07314) were analyzed together.

Chemical Analysis and Quality Assurance
The sediment samples were freeze-dried to a constant weight, ground with an agate mortar and pestle, passed through a 63 μm sieve, and stored at −20 °C in brown glass bottles. Heavy metal microwave digestion was conducted in accordance with the EPA Method 3050B using an Ethos Plus Microwave Labstation (Milestone Inc., Italy).
The marine organism samples were first thawed naturally at room temperature, then rinsed with deionized water for 3-5 times. The edible tissues were removed from each organism with stainless steel scissors and a scalpel, tissues form organisms of the same species were mixed together thoroughly, crushed with a stainless steel mixer, freeze-dried to constant weight, and ground to give a homogeneous powder, which was stored in a brown glass bottle at −20 °C. An aliquot of each dried and homogenized sample was digested using a microwave digestion system. The metal extraction procedure is described in detail elsewhere [1].
The quality control and assurance of sediment samples were carried out by using the standard recovery method. The recoveries of the metals were between 83% and 111%. The Chinese National Standard material (offshore marine sediment, GBW 07314) were analyzed together.
Reagent blanks were used to correct the analyses results of marine organism samples. The concentration of all seven metals in the standard muscle tissue reference material 2976 (for trace elements and methylmercury) supplied by the United States National Institute of Standards and Technology (Gaithersburg, MD, USA) were determined, as well as on samples for quality control. The recoveries of the metals were between 80% and 120%. The Cu, Zn, Cd, and Cr concentrations in samples were measured by an atomic absorption spectrophotometer (AAS-Z2000; Hitachi, Tokyo, Japan) and the As and Hg concentrations Reagent blanks were used to correct the analyses results of marine organism samples. The concentration of all seven metals in the standard muscle tissue reference material 2976 (for trace elements and methylmercury) supplied by the United States National Institute of Standards and Technology (Gaithersburg, MD, USA) were determined, as well as on samples for quality control. The recoveries of the metals were between 80% and 120%. The Cu, Zn, Cd, and Cr concentrations in samples were measured by an atomic absorption spectrophotometer (AAS-Z2000; Hitachi, Tokyo, Japan) and the As and Hg concentrations were measured by an atomic fluorescence spectrophotometer (AFS-9700; Haiguang Instruments, Beijing, China).

Assessment of Sediment Contamination
The index I geo has been widely used to assess heavy metal contamination in sediment and soil [19,20]. The I geo values were calculated by using the following equation: where C n represents the measured concentration of heavy metals and B n represents the background reference concentration. Due to lack of data on background values for the study area, element abundances of the upper continental crust were used as the background reference [19]. For this, 1.5 was the background matrix correction factor due to lithogenic effects. The geo-accumulation index can be divided into 7 categories, as shown in Table S2.
The methodology developed by Hakanson was used to calculate the potential ecological risk (RI) of Daya Bay [21].
where C i and B i represent the determined concentration and reference standard value of heavy metals, respectively. C i f , the pollution index of an individual element; T i r , the toxicity coefficient of heavy metals, for Cu, Zn, Cr, Cd, As, and Hg were 5, 1, 2, 30, 10, and 40, respectively [22]; and E i r is the ecological risk coefficient of an individual element.

Health Risk Models
Marine organisms accumulate heavy metals through respiration and ingestion, and heavy metals may endanger human health through the food chains. This study used health risk assessment index to evaluate the edible health risks of Daya Bay organisms. The calculation method was as follows: where C m is the concentration of heavy metals in marine organisms (mg/kg) and IR is the ingestion rate of seafood (g/person/day). The daily consumption of seafood in Guangdong Province is 51 g for an adult, and the ingestion rate for children is estimated to be 26.5% of that of adults [23]. BW is the average body weight (adults: 70 kg; children: 15 kg) [24], EF and ED indicate the exposure frequency (365 day/year) and exposure duration (adults: 70 years; children: 6 years), respectively. AT is the total exposure time (adults: 25,550 days; children: 2190 days). AR m is the intestinal absorption rate of heavy metals, the values of AR m for Cu, Zn, Cd, Cr, Hg and As were 45%, 65%, 75%, 55%, 64%, and 75%, respectively [25,26]. Rf D is the chronic reference dose, the chronic reference doses of Cu, Zn, Cd, Cr, Hg and As were 0.04, 0.3, 0.001, 0.003, 0.0001, and 0.0003 mg/kg/day, respectively; SF is the carcinogenic slope factor, based on values of 0.38, 0.5, and 1.5 mg/kg/day for Cd, Cr, and As, respectively [10]. Moreover, EDI indicates the established daily intake, HQ is the noncarcinogenic hazard quotient of heavy metals to human body (THQ indicates the total hazard quotient), and CR is the carcinogenic risk quotient (TCR indicates the total carcinogenic risk). In addition, the As concentration in this study was the total content, but previous studies on the health risks of As reported that only inorganic As posed threats to humans, and the inorganic As was assumed to be 10% of the total As [27].

Monte Carlo Simulation
When the uncertainty arises owing to lack of knowledge of the metal concentrations and variability of toxic response among individuals, the probability approximate solution obtained by Monte Carlo simulation can enhance the understanding of the environmental behavior of contaminants and show the uncertainties [28,29]. To capture variations of risk analysis in Guangdong costal, Monte Carlo simulation was used to handle the uncertainty of heavy metal concentrations in surface sediments and marine organism samples collected in Daya Bay. Simultaneously, sensitivity analysis was carried out to estimate the contribution of the exposure variables of the total risks to discern critical exposure factors. The main steps of Monte Carlo simulation are as following: (1) determine the random variables selected by the evaluation model; (2) define the distribution models of the random variables; and (3) sample randomly from the above distribution of the variables and output the simulation results. The concentrations of heavy metals in sediments and organism samples were transformed by natural logarithm to satisfy normality, indicating that the concentration data obey lognormal distribution (Table S3). Although, IR and BW, exposure time parameters (e.g., EF, ED, and AT) and toxicity criteria of heavy metals (e.g., Rf D and SF) could also have been changed, the current study only considered the impact of exposed heavy metal concentrations on risks, thus other exposure parameters were single-point inputs.

Statistical Analysis
SPSS (version 24.00) was used to transform the metal concentration data into natural logarithm and Q-Q plot was used to evaluate the normality of the dataset of each variable. Monte Carlo simulation was performed using Oracle Crystal Ball (version 11.1.4323.0) loaded in Microsoft Excel to run 10,000 times to obtain stable risk outputs and sensitivity index. ArcGIS 10.2 and OriginPro 2017 were used for graphics delineation.

Heavy Metal Concentrations in Sediments
The concentrations of Cu, Zn, Cd, Cr, Hg, and As in sediments of Daya Bay were in the ranges of 5.3-44.5, 37.1-109.0, 0.06-0.22, 13.6-72.5, 0.014-0.171, and 4.8-9.7 mg/kg, respectively ( Table 2). The mean concentrations of these elements followed a decreasing order: Zn > Cr > Cu > As > Cd > Hg. From the perspective of coefficient of variation (CV), Hg and Cu reached the CV level of 49% and 47%, respectively, and the coefficients of variation of other elements were between 14% and 31%. The results of descriptive statistics indicated that the metal concentrations in sediments of Daya Bay were quite different, and there was a large uncertainty in the spatial distribution of the metal concentration levels. Compared with the limited values of Chinese Marine Sediment Quality Standard (GB18668-2002), most of the concentration levels were lower than the corresponding limits. The mean concentrations of Zn, Cd, Hg, and As were slightly higher than their background values, indicating the impact of human activities. Threshold effect level (TEL) and probable effect level (PEL) could be used to evaluate the environmental risk caused by heavy metals in sediments. When concentrations of metals in sediments are lower than TEL, adverse effects upon sediment dwelling fauna will be infrequent, while when the PEL is exceeded, adverse harmful effects on sediment-dwelling organisms are expected to occur frequently [30]. The concentrations of heavy metal in the sediments of Daya Bay were lower than TEL or between TEL and PEL.
Comparing the results of this study with the studies of the Daya Bay in different periods (Table 1), it could be found that all concentrations of metals in sediments after 2008 were higher than those in 1988 [5]. From 2008 to 2020, the metals concentrations showed a law of first decline, then increase, and then decrease [5,13,19]. The trend of regular fluctuations reflected the impact of urban development and the government's corresponding environmental management policies on heavy metal contamination in sediments of Daya Bay. However, the results of the one-way analysis of variance showed that there was no significant difference in the concentrations of heavy metals from 1988 to 2020 (p > 0.05), which seems to indicate that the contribution of nuclear power plants sources to the contents of heavy metals in the sediment is limited. Compared with other bays along the coast of China, the levels of heavy metal concentration in Daya Bay were slightly higher than Beibu Gulf [20], but lower than Liaodong Bay, Bohai Bay, Jiaozhou Bay, Zhoushan Bay, and Quanzhou Bay [31][32][33][34][35].   Table 3 lists the overall heavy metal concentrations of 20 organism species collected in Daya Bay. Heavy metal concentrations (mg/kg, dry weight) ranged from 0.7 to 102.0 for Cu, 12.5-274.0 for Zn, 0.011-1.690 for Cd, 0.10-1.99 for Cr, 0.4-4.5 for As, and 0.015-0.147 for Hg. The relative ranking order of the mean concentrations of overall marine organism samples were: Zn > Cu > As > Cr ≈ Cd > Hg. The concentrations of mollusks in this study were similar to the results of previous study by Yuan et al. [10]. Comparing the mean concentrations of different marine organism groups indicated that mollusks had the highest Cr and As concentrations; crustaceans had the highest concentrations of Cu, Zn, and Cd; and fish had the highest concentrations of Hg. The variance in the concentration of heavy metals in different groups of organisms is related to their feeding patterns and nutritional levels [8,26]. Additionally, this study did not determine the sexes and reproductive stages of marine organisms, and the accumulation of metals in gonad tissues as well as in eggs at different embryonic development stages will have a slight impact on the statistical heavy metal contents of marine organisms [37,38]. In terms of the mean concentrations, the heavy metals concentrations in the marine organisms from Daya Bay basically met the safety standards of the Chinese National Standard of Maximum Levels of Contaminants in Food specifies (GB 2762-2017) and the safety standards of the FAO guidelines [39]. When comparing this study with the concentrations of heavy metals in marine organisms in other regions of the world, it was found that the concentrations of heavy metals in marine organisms varied widely in terms of regions and species (Table 4). Overall, the concentrations of Cd, Cr, and Hg in the marine organisms from Daya Bay were comparable to those in other regions, and the concentrations of As were much lower than those in Xiangshan Bay, Sanmen Bay, and Bangladesh coast. In addition, the concentrations of Cu and Zn in mollusks were also lower than those in other regions except the South China Sea. The geographical variation in heavy metal concentrations in marine organisms was commonly related to the local pollution status [40]. The high concentrations of metal elements in seawater and sediments contributed to the bioaccumulation of marine organisms [32].     For the contamination characteristics of heavy metals in the sediments of Daya Bay, the statistical results of random sampling based on the elements' concentrations distribution through Monte Carlo simulation and the class distribution of I geo values are shown in Table  S4 and Figure 3. The average values of I geo were arranged in descending order: Hg ≈ As > Cd > Zn > Cu > Cr, and the average values of I geo of all elements were less than zero, which is similar to the evaluation results based on the geo-accumulation index method. However, the class distribution of the I geo values revealed that there was relative variance in the spatial distribution of heavy metal contamination. The I geo value of Hg was the highest, with a moderately polluted score of 5.46% and unpolluted to moderately polluted score of 39.52%, followed by Cd, As, and Zn, with unpolluted to moderately polluted accounting for 35.71%, 34.78%, and 12.76%, respectively. Hg has no biological function, its most important sources are fuel combustion and coal combustion [12]. Previous studies reported that Hg and As in Daya Bay were affected by the same anthropogenic input [19], and Cd originated from the discharge of domestic sewage [13]. Therefore, heavy metals contamination in Daya Bay was closely related to nearby industrial and agricultural activities. Figure 4 illustrated the potential ecological risks of the heavy metals in sediments of Daya Bay. In general, the RI values of sampling sites varied significantly, ranging from 58.76 to 202.48. According to the category of RI (Table S5), five sampling sites reached a moderate risk and 22 sampling sites reached a low risk. In addition, Figure 4 shows that the E i r values of Hg had the highest contribution to potential ecological risks, followed by Cd. For the potential ecological risk characteristics, the cumulative probability of RI of heavy metals in sediments of Daya Bay based on Monte Carlo simulation were 80.94% for low risk, 18.96% for moderate risk, and 0.1% for high risk, respectively. High I geo values also corresponded to high potential ecological risk indices. Hg, Cd, and As contributed 50.85%, 33.92%, and 11.47% of the potential ecological risk, respectively (Figure 5b). It has been found in many studies that heavy metals negatively affect microbial community diversity [46]. Indirect ecological risks posed by changes in microbial communities are difficult to predict [47]. However, the concentrations of heavy metals, such as Hg, measured in this study were relatively low, and did not exceed the primary standard of the China Marine Sediment Quality standards. Even at low concentrations, the comprehensive evaluation results of I geo and RI showed that Hg was the primary heavy metal pollutant in sediments, followed by Cd and As. These results were similar to the study of Ding et al. [31]. From this point of view, Hg emissions in the Daya Bay area should be effectively controlled, and it is necessary to conduct ecological risk analysis and evaluation on the uncertainty of heavy metal emissions in the Daya Bay area.

Health Risk Assessment
The EDI of the heavy metals determined as described above are shown in Table 5. For both adults and children, EDI values were below the acceptable daily intake (ADI) and HQ values were generally in decreasing order as follows: Hg > As > Cu > Cd > Zn > Cr, and the CR values for different metals generally followed as As > Cr > Cd. Monte Carlo simulation was used to evaluate the total noncarcinogenic risks (THQ) and carcinogenic risks (TCR) (Figure 6a,b). Potential risks were indicated by the probabilities corresponding to THQ values higher than 1 (THQ > 1). It was found that there was no potential noncarcinogenic risk for adults and children. But for total carcinogenic risks, the probability of adult and children exceeded the unacceptable risk value (1 × 10 −4 ) at 11.19% and 22.27%, respectively. Therefore, the continuous or excessive intake of these marine organisms could cause chronic carcinogenic effects. The health risk results were comparable to the results reported by Yuan et al. [10]. It should be noted that the THQ and TCR used for risk assessment were based on the concentration addition of different metals, and their calculation is a noninteractive process. Humans are exposed to more than one pollutant and suffer combined or interactive effects [48]. However, THQ and TCR cannot directly reflect the true risk because they do not define any dose-response relationship [49]. Furthermore, the loss of elements caused by cooking and human intestinal absorption rate was not considered here. .05 × 10 −2 1 Average wet weight concentrations calculated after conversion from dry weight. 2 Acceptable daily intake [6]. Figure 7 presents the sensitive analysis results of the metals concentrations of THQ and TCR estimates for humans. The highest contribution to THQ was Hg (28.45%), followed by Cd (22.75%), Cu (22.67%), and As (17.84%). Noncarcinogenic risks of Hg are also expected, because Hg exposure to human is believed to be significantly correlated with the consumption of fish among all food stuffs [8,50]. The contribution of As, Cd, and Cr to the TCR value was relatively similar, with the highest being As (37.16%), followed by Cd (32.47%) and Cr (30.37%), respectively. The carcinogenic health risk results of children and adults simulated by Monte Carlo warn people to pay more attention to the contents

Health Risk Assessment
The EDI of the heavy metals determined as described above are shown in Table 5. For both adults and children, EDI values were below the acceptable daily intake (ADI) and HQ values were generally in decreasing order as follows: Hg > As > Cu > Cd > Zn > Cr, and the CR values for different metals generally followed as As > Cr > Cd. Monte Carlo simulation was used to evaluate the total noncarcinogenic risks (THQ) and carcinogenic risks (TCR) (Figure 6a,b). Potential risks were indicated by the probabilities corresponding to THQ values higher than 1 (THQ > 1). It was found that there was no potential noncarcinogenic risk for adults and children. But for total carcinogenic risks, the probability of adult and children exceeded the unacceptable risk value (1 × 10 −4 ) at 11.19% and 22.27%, respectively. Therefore, the continuous or excessive intake of these marine organisms could cause chronic carcinogenic effects. The health risk results were comparable to the results reported by Yuan et al. [10]. It should be noted that the THQ and TCR used for risk assessment were based on the concentration addition of different metals, and their calculation is a noninteractive process. Humans are exposed to more than one pollutant and suffer combined or interactive effects [48]. However, THQ and TCR cannot directly reflect the true risk because they do not define any dose-response relationship [49]. Furthermore, the loss of elements caused by cooking and human intestinal absorption rate was not considered here. of carcinogenic elements in marine organisms. There were reports that the enrichment and amplification of heavy metals in organisms will harm human health through the food chain. For example, As can cause diseases of the cardiovascular system, central nervous system, and hematopoietic system [24]. Cd can lead to high blood pressure, high risk of fracture [51], and liver dysfunction [9]. Cr is involved in lipid metabolism and insulin function, and can cause damage to DNA and tissue structure [52]. Additionally, Table 2 presented that the concentrations of Cr, As, and Cd in mollusks were higher. Crustaceans have the highest concentrations of Cu, Zn, and Cd, and fish have the highest concentrations of Hg. Therefore, mollusks were the main species that caused carcinogenic risk. Since the food habits could not be easily changed, and seafood was also the main source of protein for residents in coastal areas, on the basis of the same overall consumption of seafood, the carcinogenic risk can be maintained below a safe level by reducing the intake of mollusks and supplementing the intake of crustaceans. Moreover, the intake of crustaceans also contributed to the supplement of Zn in the human body. Cu and Zn were trace elements necessary for life activities. The moderate intake of Zn from marine organisms may therefore promote health [53].    Figure 7 presents the sensitive analysis results of the metals concentrations of THQ and TCR estimates for humans. The highest contribution to THQ was Hg (28.45%), followed by Cd (22.75%), Cu (22.67%), and As (17.84%). Noncarcinogenic risks of Hg are also expected, because Hg exposure to human is believed to be significantly correlated with the consumption of fish among all food stuffs [8,50]. The contribution of As, Cd, and Cr to the TCR value was relatively similar, with the highest being As (37.16%), followed by Cd (32.47%) and Cr (30.37%), respectively. The carcinogenic health risk results of children and adults simulated by Monte Carlo warn people to pay more attention to the contents of carcinogenic elements in marine organisms. There were reports that the enrichment and amplification of heavy metals in organisms will harm human health through the food chain. For example, As can cause diseases of the cardiovascular system, central nervous system, and hematopoietic system [24]. Cd can lead to high blood pressure, high risk of fracture [51], and liver dysfunction [9]. Cr is involved in lipid metabolism and insulin function, and can cause damage to DNA and tissue structure [52]. Additionally, Table 1 presented that the concentrations of Cr, As, and Cd in mollusks were higher. Crustaceans have the highest concentrations of Cu, Zn, and Cd, and fish have the highest concentrations of Hg. Therefore, mollusks were the main species that caused carcinogenic risk. Since the food habits could not be easily changed, and seafood was also the main source of protein for residents in coastal areas, on the basis of the same overall consumption of seafood, the carcinogenic risk can be maintained below a safe level by reducing the intake of mollusks and supplementing the intake of crustaceans. Moreover, the intake of crustaceans also contributed to the supplement of Zn in the human body. Cu and Zn were trace elements necessary for life activities. The moderate intake of Zn from marine organisms may therefore promote health [53]. of carcinogenic elements in marine organisms. There were reports that the enrichmen and amplification of heavy metals in organisms will harm human health through the foo chain. For example, As can cause diseases of the cardiovascular system, central nervou system, and hematopoietic system [24]. Cd can lead to high blood pressure, high risk o fracture [51], and liver dysfunction [9]. Cr is involved in lipid metabolism and insuli function, and can cause damage to DNA and tissue structure [52]. Additionally, Table  presented that the concentrations of Cr, As, and Cd in mollusks were higher. Crustacean have the highest concentrations of Cu, Zn, and Cd, and fish have the highest concentra tions of Hg. Therefore, mollusks were the main species that caused carcinogenic risk. Sinc the food habits could not be easily changed, and seafood was also the main source of pro tein for residents in coastal areas, on the basis of the same overall consumption of seafood the carcinogenic risk can be maintained below a safe level by reducing the intake of mo lusks and supplementing the intake of crustaceans. Moreover, the intake of crustacean also contributed to the supplement of Zn in the human body. Cu and Zn were trace ele ments necessary for life activities. The moderate intake of Zn from marine organisms ma therefore promote health [53].

Limitations
For heavy metals in sediments, only certain sampling points were used to estimate the overall risk of the area, and the results were usually uncertain, especially when the study area is large and there are many pollution sources. A Monte Carlo method, based on mathematical statistics and probability theory, was used to estimate the distribution of metal concentrations for the purpose of evaluating uncertainty in environmental risk models caused by random sampling [54]. However, the traditional Hakanson risk assessment method does not incorporate the ecosystem (species and their ecological functions) into the RI as a whole [22]. Therefore, the evaluation results cannot determine the ecological components that may be affected by heavy metals, and specify the species that need to be clearly protected. In addition, due to the spatial heterogeneity of soil/sediment, although the evaluation results can reflect the uncertainty of the contamination degree in the study area to a certain extent, it cannot combine geographic information to reveal the temporal and spatial distribution of contamination in the entire region. Moreover, the background information, such as the influence of pH, Eh, and SPM in sediments on the toxic effects of heavy metals, should also be taken into consideration.
It should be pointed out that although this study collected representative economic marine organisms in Daya Bay, the evaluation results do not represent all marine species. In addition, heavy metals can be enriched through the food chains and endanger human health. Therefore, the bioavailability and chemical speciation of heavy metals are also key factors affecting risk assessment [55]. The measurement in this study was the total concentration of marine organisms. Although the concentration and bioavailability of heavy metals have been treated with reference to a large number of previous studies, the results of the risk assessment were still uncertain. Moreover, risk models and key exposure parameters (C m , BW, IR, Rf D, and SF) are also two main sources that affect the evaluation results [28]. The Rf D and SF based upon the Integrated Risk Information System and International Agency for Research on Cancer might not yield a precise picture of current toxic responses of the Chinese population [29]. On the other hand, more precise health risks should consider the biological half-life of heavy metals, as well as the differences in the interaction and metabolic processes of heavy metals in the human body.

Conclusions
This study combined Monte Carlo simulation with risk assessment methods, and converted the traditional single point calculations into the probability distribution of risk description, providing a more comprehensive and objective assessment result of the heavy metal contamination in Daya Bay. The contents of heavy metals in the sediments of Daya Bay were lower than that in other regions of China. However, from the evaluation results of I geo and RI, Daya Bay sediment was found to be most seriously contaminated with Hg, which reached moderate pollution and accounted for 50.85% of the potential ecological risk, followed by Cd and As, which were unpolluted to moderately polluted, accounting for 33.92% and 11.47% of the potential ecological risk, respectively. Although the consumption of marine organisms in Daya Bay would not cause noncarcinogenic risks to the human body, the probability that the carcinogenic risk of adults and children exceeded the threshold (1 × 10 −4 ) was 11.19% and 22.27%, respectively. Sensitivity analysis pointed out that the contribution of carcinogenic elements to risk can be arranged in descending order as As > Cd > Cr. In summary, efforts to control and remediate extensive heavy metal contamination in the Daya Bay should focus on Hg, As, and Cd. Taking into account the uncertainty of the concentration information and the parameters in the evaluation models, the risk assessment of the coastal marine environment and human health should focus on the analysis of the bioavailability and toxic interaction of heavy metals in the future.
Supplementary Materials: The following are available online at https://www.mdpi.com/2077-131 2/9/1/17/s1, Table S1: List and information of marine organisms sampled in the Daya Bay., Table S2: Seven classes of the geo-accumulation index., Table S3: Uncertain concentration (mg/kg) of heavy metals in sediments and marine organisms in Daya Bay., Table S4: Statistical results of I geo and RI based on Monte Carlo simulation., Table S5: Category of potential ecological risks. Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Data Availability Statement:
The data presented in this study are available on request from the corresponding authors.