Risk Assessment of a Coastal Ecosystem from SW Spain Exposed to CO 2 Enrichment Conditions

: The Weight-of-Evidence (WOE) approach uses multiple lines of evidence to analyze the adverse effects associated with CO 2 enrichment in two stations from the Gulf of C á diz (Spain) with different contamination degrees. Sediment contamination and metal (loid) mobility, toxicity, ecological integrity, and bioaccumulation from the samples exposed to different acidiﬁcation scenarios (pH gradient from 8.0 to 6.0) were used in the WOE. The experiments were conducted under laboratory conditions using a CO 2 -bubbling system. Different integration approaches such as multivariate analyses were used to evaluate the results. The results indicated that the adverse biological effects under pH 6.5 were related to the mobility of dissolved elements (As, Fe, Cu, Ni, and Zn). Furthermore, the pH reduction was correlated to the increase of bioaccumulation of As, Cr, Cu, Fe, and Ni in the tissues of mussels at pH 7.0. The noncontaminated sediment showed environmental degradation related to the acidiﬁcation at pH values of 7.0; whereas the sediment moderately contaminated showed both environmental risks, caused by acidiﬁcation and the presence and the increase of the bioavailability of contaminants. The WOE approach supposes an effective tool to identify and distinguish the causes of adverse effects related to the enrichment of CO 2 in marine environments.


Introduction
Understanding the impacts of climate change and global warming is one of the greatest challenges faced by the scientific community, as are promoting mitigation and adaptation strategies [1,2]. The release of carbon dioxide (CO 2 ) from human activities is considered a significant contributor to observed changes in global climate patterns and the acidification of the oceans [3]. In addition, ocean acidification associated with global warming is projected to amplify the adverse effects of warming (further at 2 • C), provoking impacts on the growth, development, survivorship, calcification, and richness/abundance of species [4]. Therefore, in the last years, many studies analyzed the potential effects of acidification and CO 2 enrichment in marine environments [5][6][7][8][9] among natural processes [10]. They demonstrated a relationship between adverse effects on marine organisms promoted by acidification in seawater and the increase in the bioavailable forms of toxic elements. However, not many of them were designed or executed from an integrative point of view integrating multiple lines of evidence (LOEs): contamination, toxicity, and bioavailability.
The integrative methods based on the weight-of-evidence (WOE) were applied to assess the environmental quality of aquatic ecosystems in the last decades [11][12][13][14][15]. This approach uses a different set of data obtained within the application of multiple lines of evidence (LOEs) that usually implies (1) sediment physical characteristics (grain-size distribution); (2) chemical factors (contaminants in specific phases); (3) ecotoxicologically based on the use of acute and chronic toxicity tests and biomarker responses in key organisms; (4) ecological factors through the analysis of the structure and function of benthic communities; and (5) the bioaccumulation assessment of contaminants in living organisms.
The risk assessment based on the WOE approach to assessing the pollution derived from CO 2 enrichment in aquatic environments was one of the adopted methods. It integrates the different responses, whilst it reduces the uncertainty associated with the complexity [10].
A recent study used the WOE to assess the environmental risks related to the enrichment of CO 2 in the coastal zones in Brazil [9]. Results showed that the elemental bioavailability increased from pH 7.0, being responsible for the adverse effects caused in marine organisms. In this sense, the current study uses the WOE integrating LOEs to understand the potential adverse effects associated with the enrichment of CO 2 in marine ecosystems, and to identify the source of contamination responsible for the pollution.
The main goal of the current study was to determine the environmental risk of CO 2 enrichment in marine ecosystems using the WOE approach. Therefore, different LOEs from two sampling sites with different contamination degrees were employed under the influence of acidification (pH gradient from 8.0 to 6.0).

Study Area
Two sediment samples were collected from the Bay of Cádiz ( Figure 1). The selected areas were previously studied and were classified according to the metal contamination levels [5,16,17]. The Rio San Pedro (36 • 31 52.90" N; 6 • 12 48.43" W) sampling site is classified as low metal concentration (control site), while the Trocadero (36 • 31 15.5" N; 6 • 12 25.1" W) is considered a sampling site with a moderate metal concentration (contaminated site) [16,18]. tribution); (2) chemical factors (contaminants in specific phases); (3) ecotoxicologically based on the use of acute and chronic toxicity tests and biomarker responses in key organisms; (4) ecological factors through the analysis of the structure and function of benthic communities; and (5) the bioaccumulation assessment of contaminants in living organisms.
The risk assessment based on the WOE approach to assessing the pollution derived from CO2 enrichment in aquatic environments was one of the adopted methods. It integrates the different responses, whilst it reduces the uncertainty associated with the complexity [10].
A recent study used the WOE to assess the environmental risks related to the enrichment of CO2 in the coastal zones in Brazil [9]. Results showed that the elemental bioavailability increased from pH 7.0, being responsible for the adverse effects caused in marine organisms. In this sense, the current study uses the WOE integrating LOEs to understand the potential adverse effects associated with the enrichment of CO2 in marine ecosystems, and to identify the source of contamination responsible for the pollution.
The main goal of the current study was to determine the environmental risk of CO2 enrichment in marine ecosystems using the WOE approach. Therefore, different LOEs from two sampling sites with different contamination degrees were employed under the influence of acidification (pH gradient from 8.0 to 6.0).
Water used in toxicity tests was collected from the surface (1 m depth) of the Rio San Pedro (RSP), during high tide, transported to the laboratory using 50 L gallons, and placed in a 400 L tank (pH 7.9 ± 1; Salinity 34 psu, dissolved oxygen >81%, and temperature 18 ± 1 °C). As previously mentioned, the RSP has been used as a reference point in other works since no contaminants are present.  Water used in toxicity tests was collected from the surface (1 m depth) of the Rio San Pedro (RSP), during high tide, transported to the laboratory using 50 L gallons, and placed in a 400 L tank (pH 7.9 ± 1; Salinity 34 psu, dissolved oxygen >81%, and temperature 18 ± 1 • C). As previously mentioned, the RSP has been used as a reference point in other works since no contaminants are present. Figure 2 shows, in detail, the schematic design used for each experiment included in the WOE: (i) design to assess the benthic integrity; (ii) design used to elutriate assays employed as samples prior to acidification; (iii), and (iv) acute toxicity tests with amphipods Appl. Sci. 2023, 13, 5805 3 of 15 and mussels design, respectively. In addition, it was adapted (iv) to be used for the bioaccumulation and NRRT assays. Appl. Sci. 2023, 13, Figure 2 shows in detail the schematic design used for each experiment included in the WOE: (i) design to assess the benthic integrity; (ii) design used to elutriate assays employed samples prior acidified; (iii) and (iv) acute toxicity tests with amphipods and mussels design, respectively. Also, it was adapted (iv) to be used for the bioaccumulation and NRRT assays.

Figure 2.
Schematic design of the CO2 injection system used in this study for the different LOEs: (i) the benthic integrity; (ii) the elutriate assays for chronic toxicity; (iii,iv) acute toxicity tests with amphipods and mussels, respectively, and (iv) it was also applied for bioaccumulation and NRRT assays.
The laboratory-scale injection systems employed ( Figure 2) for ecotoxicological and ecological integrity assessment were described by [16] and chemical analyses, acute and chronic toxicity tests, ecological integrity and metal bioaccumulation are further described in [18,19]. The CO2 injection system proposed mimics acidification by CO2 enrichment in marine ecosystems.

Chemical Analyses
Sediment and water samples were collected from each station (RSP and TRO, Figure  1) at the beginning and the end of the toxicity and mesocosms tests. The metal and metalloid concentrations (Al, As, Cr, Cu, Fe, Ni, and Zn) in the sediment were analysed in the total fraction. The dissolved metal concentration in overlying and elutriate seawater were determined in filtered (0.45 mm) and acidified (pH < 2 with ultrapure HNO3) subsamples for metal(loid) analyses.

Toxicity Tests
Different biological responses were selected: (i) acute toxicity tests: amphipods and mussels mortalities; (ii) chronic toxicity tests: embryo-larval development success of seaurchin and mussels; (iii) metal bioaccumulation in soft body of mussels; (iv) macro-   Figure 2. Schematic design of the CO 2 injection system used in this study for the different LOEs: (i) the benthic integrity; (ii) the elutriate assays for chronic toxicity; (iii,iv) acute toxicity tests with amphipods and mussels, respectively, and (iv) it was also applied for bioaccumulation and NRRT assays.
The laboratory-scale injection systems employed ( Figure 2) for ecotoxicological and ecological integrity assessments were described by [16] and chemical analyses, acute and chronic toxicity tests, ecological integrity, and metal bioaccumulation are further described in [18,19]. The CO 2 injection system proposed mimic acidification by CO 2 enrichment in marine ecosystems.

Chemical Analyses
Sediment and water samples were collected from each station (RSP and TRO, Figure 1) at the beginning and the end of the toxicity and mesocosms tests. The metal and metalloid concentrations (Al, As, Cr, Cu, Fe, Ni, and Zn) in the sediment were analysed in the total fraction. The dissolved metal concentration in overlying and elutriate seawater were determined in filtered (0.45 mm) and acidified (pH < 2 with ultrapure HNO 3 ) subsamples for metal(loid) analyses.

Toxicity Tests
Different biological responses were selected: (i) acute toxicity tests: amphipod and mussel mortalities; (ii) chronic toxicity tests: embryo-larval development success of sea urchins and mussels; (iii) metal bioaccumulation in the soft bodies of mussels; (iv) macrobenthic community alteration and; (v) biomarker response of cellular stress: neutral red retention time (NRRT) in the hemolymph of mussels.
The pH values of the treatments ranged between 8.0 and 6.0 (8.0, 7.5, 7.0, 6.5±, and 6.0) and were conducted in replicates. Table 1 summarises the description of the different toxicity tests performed in the study using different species and conditions and the ecological integrity assessment conditions.

Acute-Toxicity Tests
Individuals of amphipods Ampelisca brevicornis were collected from a noncontaminated area at the Bay of Cádiz. These are used as environmental quality indicators and are recommended by the specific protocol for dredged-sediment characterization [20]. Individuals of mussels Mytilus galloprovincialis were obtained from an aquaculture farm (Cetária, Cádiz, Spain). Both organisms were acclimated for a week before testing. Parameters were adjusted according to the optimal conditions established for each species employed in this study and kept until the end of the experiments ( Table 1). The experiments were adapted from [21,22].

Chronic Toxicity Using Elutriates
Specimens of sea urchin Paracentrotus lividus were collected from a rocky intertidal platform from a protected area in the Gulf of Cádiz (Southwestern Spain); whereas the specimens of mussels M. galloprovincialis were obtained from an aquaculture farm (Cetária, Cádiz, Spain).
Sea urchin gametes were obtained from direct extraction from the gonads. On the other hand, mussels were induced to spawn by thermal stimulation according to the protocol recommended by [23,24] and exposed to acidified sediment elutriate following [22].

Metal Bioaccumulation
M. galloprovincialis mussels were exposed to acidified sediments for 14 days. Assays were performed using whole sediment with a sediment-to-water ratio of 1:2 v/v as reported by [25]. The concentrations of metals Zn, Cu, Al, Fe, Mn, Ni, Cr, and the metalloid As in the soft tissue were determined using an inductively coupled plasma mass spectrometer (ICP-MS) (Thermo Elemental Series-X, Santa Clara, CA, USA).

Ecological Integrity
After the acclimation time, mesocosms were maintained for 21 days, natural photoperiod and constant aeration and temperature. Then, the macrobenthic community structure (total abundance, richness, and diversity) was determined using Margalef's and Shannon's methods, respectively.

Neutral Red Retention Time Assay
The NRRT assay was based on the principle that only lysosomes in healthy cells take up and retain the neutral red vital stain [26]. Lysosomal membrane damage caused by xenobiotics can decrease the NRRTs [27]. The experiments were performed following the method described by [28]. Figure 3 represents the four LOEs integration used to characterize and distinguish the impacts related to the contamination levels and acidification degrees.

Statistical Analysis (WOE Approach)
The obtained results from the application of each LOE were linked for each component assessed using three different methods (WOE): (a) tabular matrix based on statistical significance (p < 0.05) compared to the reference case (RSP sediment, pH 8.0); (b) pie charts, showing statistical differences (p < 0.05 significantly different and positive effect; 0.01 ≥ p > 0.05 significantly different but moderate effect; and p > 0.01 not significantly different) among the samples sites and pH values compared to the reference case (RSP sediment, pH 8.0); and (c) linking the results using a multivariate analysis using the factor analysis principal component as extraction method [29].
among the samples sites and pH values compared to the reference case (RSP sediment, pH 8.0); and (c) linking the results using a multivariate analysis using the factor analysis principal component as extraction method [29].
Significant differences found between the reference site and the pH treatments were determined by an ANOVA and Dunnet's tests (TOXSTAT software, version 3.5) to develop the pie charts and tabular matrix. The indexes of contamination, toxicity, bioaccumulation, and ecological integrity were calculated as outlined by [11], considering the ratio-to-reference (RTR) values. Results from each LOE and each pH treatment were referred to the reference case and then Significant differences found between the reference site and the pH treatments were determined by an ANOVA and Dunnet's tests (TOXSTAT software, version 3.5) to develop the pie charts and tabular matrix.
The indexes of contamination, toxicity, bioaccumulation, and ecological integrity were calculated as outlined by [11], considering the ratio-to-reference (RTR) values. Results from each LOE and each pH treatment were referred to the reference case and then derived as the minimum value. The RTRs were averaged for contamination in sediments and mobility, acute and chronic toxicity, bioaccumulation, and ecological integrity.
The contamination index was obtained regarding the metal concentration in sediments (RTR sediment ), and concentration in seawater and elutriates (RTR mobility ). The elements As, Cr, Cu, Ni, and Zn were considered to determine the indexes. The acute toxicity index was derived from the average of RTR calculated for amphipod and mussel mortalities and NRRT assays. The values of NRRT were transformed as their inverse to be treated with equal importance. The chronic toxicity index was estimated from the average RTR for larval development inhibition results of mussels and sea urchins. Regarding the index of macrobenthic community structure (ecological integrity), the values were also inverted to be treated with equal importance, and the average of RTR was calculated based on the results of total abundance, diversity, and richness. The bioaccumulation index was derived from the average RTR estimated for the metals bioaccumulated in the tissues of mussels.
Furthermore, to develop an integrated and more precise interpretation of the risks associated with CO 2 enrichment in the marine environment, a multivariate analysis (principalcomponent analysis, PCA-with-varimax rotation) was carried out. Therefore, the correlation between increasing pollution, pH reduction, and biological effects measured in the laboratory was further explained. Factor scores for each sampling station were calculated following the protocols outlined by [11,30] and using the STATISTICA ® software package (version 13.2). The principal factors were extracted and the eigenvalues above 1.0 were taken into account. Variables having loadings ≥0.40 to a particular factor were considered associated with the respective factor, following [31].

Results and Discussion
The results obtained for each of the LOEs used in this study are summarized in Table 2. The concentrations of most of the studied metal(loid)s in the TRO sediment sample were higher than in RSP sediment samples. Concentrations of Al, Cr, and Fe showed little variation among the pH values for both sediment samples (RSP and TRO). The Cu concentration in sediment decreased with pH for both sampling sites, whereas an increase in Zn concentration was observed in the TRO sediment at pH 6.0. However, greater variability in the concentration of metals was observed in the samples of overlying water and elutriates when compared with the sediment samples. Overall, the metal concentrations were higher in the elutriate treatment than in the overlying seawater. The concentrations of the studied elements were higher in sediment from TRO elutriates, while RSP sediments exhibited higher concentrations for most of the metals than for the overlying seawater samples. Table 2. Summarized results for the different lines of evidence (LOEs) used in this study: a) chemical concentrations in sediment (sed, in µg g −1 dry weight), overlying water (sw, in mg L −1 ), and elutriate (elu, in mg L −1 ), and the bioaccumulation of elements in the soft tissue of mussels (bio, in µg g −1 dry weight); b) acute and chronic toxicity tests and neutral red assays (mortality of amphipods (Mort_A), and mortality of mussels (Mort_M); inhibition of larval development of mussels (Inh_M) and sea urchins (Inh_SU); neutral red retention time (NRRT) for 24 h, 96 h, and 10 d of exposure are expressed in minutes); and d) ecological integrity (total abundance, richness, and diversity) exposed to the pH treatments (8.0, 7.5, 7.0, 6.5, and 6.0) using the RSP and TRO sample sites. Regarding the toxicity tests, the mortality of amphipods and mussels were 14% and 20%, respectively, for the control (not CO 2 injection, RSP sediment). Significant mortality (p < 0.05) of amphipods was found at pH 7.0, and for mussels at pH 6.0 when exposed to RSP sediment. However, exposure to TRO sediment control displayed a mortality of 37% for amphipods and 35% for mussels. The absence of acidification treatment with TRO sediments implies that contamination strongly influenced the obtained results [18].
However, for the embryo-larval development test, a percentage of inhibition of less than 20% was observed for both species (sea urchin and mussel) and both sediment and elutriates samples (RSP and TRO) for the control treatment (not CO 2 injection). The mussellarval development inhibition differed significantly between the control and pH 6.5 and 6.0 for RSP elutriates; though it was statistically significant (p < 0.05) for pH 7.5 and below for TRO samples. On the other hand, the sea urchin larval-development inhibition presented significant differences (p < 0.05) at pH 7.0 for both sediment samples.
In general, a different bioaccumulation trend for each metal(loid) was shown by M. galloprovincialis with acidification. Organisms exposed to RSP samples accumulated greater concentrations of elements than those exposed to TRO samples for all the studied elements, except for Zn and Cr. Results are further detailed in [32].
Results using the ecological integrity showed a decrease in the total abundance, richness, and diversity of benthic fauna for pH 6.0 when compared to the control for both sediment samples tested (Table 2), in agreement with [33].

Tabular Matrix
Statistically significant differences (p < 0.05) between the reference treatment (RSP sediment, pH 8.0-no CO 2 added) and the other treatments and stations are presented as showing a positive response in contamination, toxicity, integrity, or bioaccumulation (Table 3). In general, the metal concentration in sediment was not significantly affected by acidification for RSP sediment. Nevertheless, significant differences (p < 0.05) were found between the metal concentration in the TRO and the control. These results were expected since the metal contamination in the TRO station was higher than in the RSP station. These results agreed with those reported by [16], which characterized the RSP as a low metalconcentration station and, the TRO as a station with moderate metal concentration. On the other hand, the results obtained in the mobility index showed significant differences (p < 0.05) between the reference case and pH 6.5 and 6.0 for the RSP and starting at pH 7.5 for the TRO. These findings indicated that acidification could increase the mobility of metals from the sediment samples to the water column, which might cause toxic effects on marine organisms. Previous studies reported the pH as a variable that mostly might influence the availability of elements from sediments [17]. Furthermore, it can be suggested that the significant difference found in the toxicity indexes at pH 7.0, using the RSP sediment, could also be related to the increase in metal mobility. Nevertheless, a significant toxic effect at pH 7.5 was observed, and, according with previous results, showed significant effects on the larval development inhibition of mussels and amphipods, and clam mortalities associated with acidification [7,22,34,35].
The results of the integrity of the benthic community showed significant (p < 0.05) differences for all treatments tested using the TRO sediment. These results were expected since the TRO sediment sample presented higher contamination than the RSP sample and the lowest value of total abundance of macrobenthic found in the study. On the other hand, the results of the bioaccumulation of metals in the tissues of mussels showed a significant (p < 0.05) increase in the concentration of metals at pH 6.5 for both sediment samples.

Sector Charts
The obtained results of the present study were also expressed in sector charts that used different levels of the statistical p-value (Figure 4) to assess the influence of ocean acidification promoted by CO 2 enrichment.  Significant differences were represented with different colours and weights for the LOEs. The RSP site was selected as the reference site at control pH (not CO2, pH 8.0 ± 0.1). Then, each case (pH treatment in the sampling site) was compared to this reference case and the significant difference respect to the reference case is represented as black segments when p < 0.05, whereas when 0.01 ≥ p > 0.05 they were represented by grey segments, and  Black colour represents a significant difference (p < 0.05) against the control, grey colour represents differences using different p-values (0.01 ≥ p > 0.05), and white colour represents no significant difference concerning the reference case. Significant differences were represented with different colours and weights for the LOEs. The RSP site was selected as the reference site at the control pH (not CO 2 , pH 8.0 ± 0.1). Then, each case (pH treatment in the sampling site) was compared to this reference case and the significant difference, with respect to the reference case, is represented as black segments when p < 0.05, whereas when 0.01 ≥ p > 0.05, they were represented by grey segments, and no significant difference (p > 0.01) was represented by white segments (Figure 4).
This approach allows clear visualization of the influence of CO 2 enrichment on marine organisms concerning the pH tested and the contamination levels in each site used. According to [36], pie charts represent the strength of the statistical tool and the simplicity of interpretation, though without quantification. Results showed that both sediment contamination and contamination caused by acidification were associated with the biological responses in this study. Significant results (black area) were observed at pH 7.0 when using noncontaminated sediment. On the other hand, the results obtained using the TRO sediments showed significant effects even in the control treatment (pH 8.0, not CO 2 added) for macrobenthic and contamination by metals in sediment. These results demonstrated that the sediment sampling sites presented the initial differences in their adverse effects associated with the contamination levels being not related to the acidification. Grey areas (moderate effects) indicate the beginning of the adverse effects and are useful for classifying cases in which additional attention should be taken, especially for potential scenarios of acidification.

Multivariate Analysis
The principal-component analysis (PCA) identified the correlations between sediment contamination, biological responses, and the pH decrease in seawater caused by CO 2 enrichment in the ocean (Table 4). A total of 30 variables were used, including dissolved metals in overlying seawater (As, Cr, Cu, Fe, Ni, and Zn) and elutriate treatment (Cu, Fe, and Zn), toxicity to marine organisms (mortalities, embryo-larval development inhibition, and NRRT time after 24 h, 96 h, and 10 d of exposure), sediment characteristics (concentrations of organic matter, fines, and metals), concentrations of metals bioaccumulated in the whole body of mussels (As, Cr, Cu, Fe, Ni, and Zn) and benthic integrity (total abundance, diversity, and richness). The two main factors explained almost 75% of the original variance in the original data set. The criterion selected to interpret a variable associated with a particular factor was a loading of 0.4 or higher. Each factor is described according to the dominant group of variables.
The first factor (F1, 48.41% of variance) linked all dissolved metals (except Fe) in seawater and elutriates, the metal concentration in mussels (except Zn), mortalities of amphipods and mussels, embryo-larval development inhibition of sea urchins and mussels and the increasing proton (H + ) concentration, in the other words, the CO 2 acidification. Furthermore, the F1 also linked negative values to the results of NRRT at 24 h, 96 h, and 10 d, and the richness of the macrobenthic community. The second factor (24.93% of variance) linked all metals in sediment, dissolved Ni in the seawater and Cu in the elutriate, As and Zn bioaccumulated in the mussels, mortality of mussels, sea urchin larval development inhibition, and the values of total abundance and diversity of the macrobenthic community. Moreover, negative values were linked in factor two for Cr and Fe in seawater.
In general, the multivariate analysis showed that both causes can influence the adverse biological effect measured: acidification by enrichment of CO 2 and the metal concentration analyzed in the sediments. −0.74 The ecological integrity of the macrobenthic community related to the total abundance and diversity showed to be associated with the sediment contamination of the sampling sites, though not with the pH reduction. However, previous studies showed significant changes in the structure and reductions in diversity, abundance, and biomass of the benthic communities exposed to acidification scenarios [8,37,38]. These results confirm the obtained results (Table 3 and Figure 4) that showed a significant (p < 0.05) difference in the TRO benthic community at pH 8.0 when compared to the reference case. Therefore, the pH level probably was not correlated with these results, and the macrobenthic community could be affected by the CO 2 enrichment in the ocean only at highly-acidified conditions.
It is known that the forms in which metals are bound in the sediments and the sediment porosity stability and grain size may influence their bioavailability and toxicity [39][40][41][42][43]. Factor one suggests that the majority of metals can increase their mobility in acidified conditions, and might cause adverse effects on marine organisms. These results were also observed in the tabular matrix and pie charts (Table 3 and Figure 4, respectively) that showed a significant difference (p < 0.05) in the mobility of metals starting at pH 6.5 for the RSP sediment and at pH 7.5 for the TRO sediment.
Furthermore, acute and chronic toxicity were linked in factor one with the metal's mobility and the increase in the concentrations of H + (enrichment of CO 2 ). These results were indicated in the tabular matrix and pie charts as significantly differenced starting at pH 7.5 and 7.0 for TRO and RSP sediment, respectively, and corroborate with the results found by [9], which showed significant effects starting at pH 7.0 in no contaminated sediment sample. On the other hand, the multivariate analysis identified which metals were correlated with their concentration increase and the toxicity caused to the organisms (As, Ni, and Zn in seawater and Cu, Fe, and Zn in elutriate sediments). These results corroborate previous studies which showed adverse effects on marine organisms associated with increases in the metal concentrations and the pH reduction [10,16,22,34,35] reported with the mortality of amphipod Hyale youngi with seawater acidification due to an increase in Ni and As concentrations. In addition, an increase in the mobility of dissolved Zn promoted mortality in bivalves (adults and juveniles) [22,25] and amphipods [16].
Furthermore, it was reported that there was a correlation between the increase in the mobility of Cu, Zn, and As with the mortality of Hediste diversicolor polychaete [44]. Furthermore, it has been observed that there is variability in the feeding behaviour of marine amphipods (O. tuberculata) in response to ocean acidification [45].
Regarding the results of bioaccumulation of metals, it was correlated that factor one increases in the concentrations of As, Cr, Cu, Fe, and Ni in the tissues of mussels, increases in the mobility of the metals mentioned above, and increases in the concentrations of H + protons. These results were described in Table 3 and Figure 4, presenting significant differences from pH 6.5 for both sediment sampling sites. These findings are similar to [44], a study that showed a correlation between the increase of Al, Fe, Mn, Cu, Zn, and As accumulated in the tissues of H. diversicolor with the pH reduction. Furthermore, the bioaccumulation of As and Zn was shown in factor two which was correlated with the sediment contamination. These results are expected since mussels exposed to trace elements in seawater can accumulate high amounts in their tissue [6]. In addition, a previous study reported that the bioaccumulation in Tellina deltoidalis showed that the bioavailability of metals was likely to be a result of the continued cycling between the pore water and the surface sediments due to the bioturbation process [42].
According to [45], the bioaccumulation of metals may be influenced by some factors, such as metal specificity, environmental influences, exposure routes, and species-specific characterization. Results of metal bioaccumulation in mussels may be associated with different parameters such as CO 2 enrichment in seawater, sediment contamination, metals behaviour, etc.
To confirm the factors' descriptions and to establish the relationship between components associated with contamination, biological effects and acidification at each studied site, it was proposed as a representation of estimated factor scores per case ( Figure 5).
The factor scores have demonstrated that the contamination levels of sediment influence the obtained results. Overall, the RSP sediment that is considered a noncontaminated sediment site did not correlate with the contamination levels and biological responses (factor one) for the higher pH values. Moreover, the results of factor scores confirm those previously reported ( Table 2) with lower levels of metal contamination in sediments (factor two).
Furthermore, as expected, the results of factor two were correlated even at the control pH values using the TRO which also confirms the results obtained using the tabular matrix and pie charts. Furthermore, it was noted that factor two decreases with the concentration of protons increase (pH reduction) while the opposite results were demonstrated for factor one, which increases its value as the pH values decrease.
Furthermore, as expected, the results of factor two were correlated even at the control pH values using the TRO which also confirms the results obtained using the tabular matrix and pie charts. Furthermore, it was noted that factor two decreases with the concentration of protons increase (pH reduction) while the opposite results were demonstrated for factor one, which increases its value as the pH values decrease. Figure 5. Representation of the estimated factor scores for each sediment sampling site (RSP and TRO) and treatment (control, 7.5, 7.0., 6.5, and 6.0).
In this sense, the pH value of 6.5 was correlated to provoke adverse effects on marine organisms that could also be associated with the increase in the mobility of metals. Moreover, it is important to perform a complete risk assessment in areas affected by CO2 enrichment considering different LOEs associated with biological responses and sediment characterization for a proper interpretation of the impacts caused by the acidification in a marine environment.

Conclusions
The present study applied for the first time a complete and specific design of a weight-of-evidence approach (WOE) to address the potential impacts of CO2 enrichment in the marine environment using different lines of evidence. It has been demonstrated from the different lines of evidence (LOEs) employed in this study using a WOE that the acidification by enrichment of CO2 provokes adverse effects on the marine organisms, starting at pH 7.0, which considerably enhances at pH 6.5.
Metal contamination in sediments seems to play a very important role in the environmental risk assessment under the enrichment of CO2 scenarios. The results showed that dissolved elements such as As, Fe, Cu, Ni, and Zn increased their concentrations and contributed to causing toxic effects on marine organisms under acidified conditions. The WOE distinguished between contamination sources that are causing biological adverse effects in marine organisms. It differs between the extension of the effects produced by the acidification by enrichment of CO2, the extension of the effects associated with the levels of contaminants (mainly metals), and the combination of both causes of adverse effects, acidification, and metal contamination at the same time.
Author Contributions: M.C.P. and I.R. contributed to all aspects of the article from the conceptualization to the funding acquisition including methodology software, validation, formal analysis, investigation, resources, data curation, and writing and editing the paper from the start to the final version. E.B. contributes in most of the aspects including conceptualization, investigation, resources, data curation, writing-original draft preparation, writing-review, and editing and visualization. A.C. contributed to most of the aspects including methodology, validation, formal analysis, writing-review and editing, visualization, supervision, and funding acquisition. All authors have read and agreed to the published version of the manuscript. In this sense, the pH value of 6.5 was correlated to provoke adverse effects on marine organisms that could also be associated with the increase in the mobility of metals. Moreover, it is important to perform a complete risk assessment in areas affected by CO 2 enrichment considering different LOEs associated with biological responses and sediment characterization for a proper interpretation of the impacts caused by the acidification in a marine environment.

Conclusions
The present study applied for the first time a complete and specific design of a weightof-evidence approach (WOE) to address the potential impacts of CO 2 enrichment in the marine environment using different lines of evidence. It has been demonstrated from the different lines of evidence (LOEs) employed in this study using a WOE that the acidification by enrichment of CO 2 provokes adverse effects on the marine organisms, starting at pH 7.0, which considerably enhances at pH 6.5.
Metal contamination in sediments seems to play a very important role in the environmental risk assessment under the enrichment of CO 2 scenarios. The results showed that dissolved elements such as As, Fe, Cu, Ni, and Zn increased their concentrations and contributed to causing toxic effects on marine organisms under acidified conditions. The WOE distinguished between contamination sources that are causing biological adverse effects in marine organisms. It differs between the extension of the effects produced by the acidification by enrichment of CO 2 , the extension of the effects associated with the levels of contaminants (mainly metals), and the combination of both causes of adverse effects, acidification, and metal contamination at the same time.