Next Article in Journal
Ballast Water Management Strategy to Reduce the Impact of Introductions by Utilizing an Empirical Risk Model
Previous Article in Journal
Appraisal of Environmental Health and Ecohydrology of Free-Flowing Aghanashini River, Karnataka, India
Previous Article in Special Issue
Environmental Aspects of Historical Ferromanganese Tailings in the Šibenik Bay, Croatia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Partitioning and Availability of Metals from Water Suspended Sediments: Potential Pollution Risk Assessment

1
Departamento de Recursos Naturales, Facultad de Zootecnia y Ecología, Universidad Autónoma de Chihuahua, Perif. Ramiro Almada km 1, Chihuahua 31453, Mexico
2
Centro de Investigación en Materiales Avanzados, CIMAV, Miguel de Cervantes 120, Chihuahua 31109, Mexico
3
Departamento de Física Aplicada II, Facultad de Arquitectura, Universidad de Sevilla, Av. Reina Mercedes, 41004 Sevilla, Spain
*
Author to whom correspondence should be addressed.
Water 2022, 14(6), 980; https://doi.org/10.3390/w14060980
Submission received: 7 February 2022 / Revised: 13 March 2022 / Accepted: 18 March 2022 / Published: 20 March 2022
(This article belongs to the Special Issue Geochemistry of Water and Sediment Ⅱ)

Abstract

:
The water management initiatives in freshwater systems focus on water availability to preserve this resource for human uses and the health of aquatic ecosystems. This work presents an assessment of the potential pollution risk caused by the metal availability in suspended sediments. The objective of this study was to determine the partitioning, association, and geochemical fractionation of metals in suspended sediments from a surface water body. Additionally, the environmental assessment for this reservoir was estimated using geoaccumulation, enrichment, and pollution indices of metals and the related potential risk by their elemental availability (RAC). Chemical, mineralogical, and morphological characterizations were obtained by inductively coupled plasma spectrometry, alpha spectroscopy, X-ray crystallography, and scanning electron microscopy. Clay, quartz, montmorillonite, and calcite were the main minerals of suspended sediments. Chemical fractionation was the parameter affecting the concentrations of metals in suspended sediments. The sediment composition is of natural origin; however, these finer particles can promote the scavenging of toxic metals. It contributes to obtaining moderate to high levels for enrichment/contamination indices. Although Ca, Mg, and U were the most accessible metals for aquatic biota, metals such as Sr, Mn, Li, Cu, and Ni in the exchangeable phase of suspended sediments are the potentially toxic elements in this aquatic ecosystem.

1. Introduction

There is an urgent need to balance social, economic, and environmental sustainability globally. Water scarcity is a severe problem due to dwindling surface and groundwater supplies caused by natural and human sources; these water changes adversely affect its availability [1]. In nature, water availability depends on climatic and geological factors [2], where their interaction has been called hydrogeologic environment [3].
Water quality conditions are essential for functioning in the ecosystem and meeting human commitments. Those systems naturally contain trace elements responsible for a fraction of the biota’s dose from the environment [4,5]. The indices such as enrichment factor, Geoaccumulation Index, contamination factor, and potential ecological risk index are used to determine the pollution risk given by metals [6,7,8,9,10,11]. However, most of those studies have been widely evaluated using the bottom sediments. Nevertheless, low-weight sediments remain in suspension, which can be mobilized through long distances, being a valuable transport source of pollutants [12,13,14].
The concentration of toxic elements is generated through lixiviation, chemical speciation, size of suspended sediments, and movement, resulting in their transport and distribution within a water body [15]. In this regard, processes such as diffusion and dispersion, transport, exchange in solid-solution phases, and sorption/desorption control the migration of the contaminants in water bodies [7,16,17,18]. Furthermore, environmental conditions influence the sorption of contaminants on suspended sediments, which can change to soluble states and be available to other organisms [19,20]. Chemically, element fractionation is controlled by chemical properties such as oxidative state, chemical affinity, chemical species competition, among others. Therefore, evaluating the availability of elements in suspended sediments should be a priority to establish the potential toxicity. The particle fractions can be [21,22,23,24]: bioavailable (exchangeable/soluble metals in water), reduced (metals linked to oxides freed when conditions change from oxic to anoxic), oxidized (metals linked to organic components and sulfides, freed under oxidant conditions), and residual (elements strongly connected to the crystalline minerals structures). Due to the chemical species, contents in gross sediments can be a poor indicator of metal releasing in aquatic ecosystems; it is critical to identify the bioavailability of metals in suspended sediments under variations of environmental conditions. The objective of the present study was to determine the partitioning, association, and geochemical fractionation of metals (including 238,234U and 232,230Th) in suspended sediments from surface water to assess the potential pollution risk for an aquatic environment.

2. Materials and Methods

2.1. Study Site

The present study was performed in the San Marcos reservoir, located between 28°55′ and 28°40′ latitude N and between 106°10′ and 106°25′ longitude W, in the municipality of Chihuahua. This reservoir is one of the oldest reservoirs in Chihuahua, Mexico, built-in 1905, before the Mexican revolution. This region is placed in a volcanic caldera, showing rhyolitic tuffs and some Upper Cenozoic volcanic rocks with intermediate to basic composition [25,26]. Likewise, basalts and a few calcareous rock outcrops can be found [25]. Due to its high topographic elevations, the San Marcos region serves as a source of water recharge for the Chihuahua–Sacramento aquifer. In this semiarid zone, short periods of rain are present with annual precipitation from 200 to 600 mm, where temperatures can reach up to 38 °C [27].
The San Marcos reservoir is fed with water from two rivers (Blanco and San Marcos rivers) and a spring, where San Marcos River is the primary source of recharge to this reservoir. This river comes in contact with two uranium deposits (Victorino and San Marcos outcrops) present in the area, contributing natural radioisotopes to this reservoir [28]. Several decades ago, some mining studies were carried out in these mineral deposits to benefit uranium extraction; however, it was concluded that there was no feasibility for uranium mining. The spring is located on the north-eastern slope of the reservoir and provides abnormal uranium and radon concentrations due to its groundwater origin [28,29]. This zone uses the surface and ground waters for agricultural and domestic activities and human consumption. The location of the study area is shown in Figure 1.

2.2. Sampling

Eighteen sampling sites were collected along the flux water of San Marcos reservoir in May 2016, before the rainy season; nine of those were taken at the surface level (40 cm below from the topwater) and the other nine at a deep level (1 m above from the reservoir bottom). Input water corresponds to sampling points 1 (San Marcos River), 3 (Blanco River), and 8 (spring). Every water sample was collected in plastic containers of 50 L; before sampling, these containers were acid washed in the lab and in situ pre-equilibrated with water from every sampling. According to Mexican regulation, 1 L of water was also collected to determine water quality parameters in every sampling point [30].

2.3. Sample Analysis

2.3.1. Water and Sediments Samples

Parameters such as pH, temperature (T), oxidation-reduction potential (ORP), and total dissolved solids (TDS) were recorded in situ. While pH, T, and TDS were obtained using multiparametric equipment (Hanna, Hl98130), the ORP was measured with an ORP tester (OAKTON, Waterproof ORPTestr Double junction). Total alkalinity, chlorides, and total hardness were determined by the methodology proposed in the Mexican regulations [31,32,33]. The analysis of sulfates and nitrates were performed through Hach procedures by 8039 (Cadmium reduction method) and 8051 (SulfaVer 4) methods and measured using a Hach spectrophotometer (Hach, DR 2010) [34], whereas Ca, Mg, Na, and K were determined using an atomic absorption spectrophotometer (GBC Avanta Sigma, Melbourne, Australia), following the methodology suggested in the regulation [35].
The 50 L water samples were exposed to successive filtrations to obtain the suspended sediments and dissolved phase. Although suspended sediments are composed of particles of different diameters, in this work, two sizes of particles were analyzed: coarse (>25 µm) and fine (25 µm > particles > 11 µm). Initially, water samples were passed through a cellulose filter with a pore size of 25 µm (coarse particles); subsequently, that water was filtered with a pore size filter of 11 µm to obtain the fine particles. For purposes of this work, the liquid phase obtained from the final filtration was also analyzed and called the dissolved phase.

2.3.2. Characterization of Suspended Sediments and Dissolved Phase

Scanning Electron Microscopy

The coarse and fine particles were characterized using SEM in a JEOL 6460LV, obtaining secondary electron images (SEI) and back-scattered electron images (BEI); with a maximum resolution of 3.5 nm. This device was coupled to a microprobe EDX and equipped with a Beryllium window ATW2 (resolution 137 EVAT 5.9 KeV) for semi-quantitative analysis using the Oxford INCA® software.

X-ray Diffraction

The suspended sediments were subjected to a mineral analysis through X-ray diffraction, using an XPERT-PRO PANalytical® (Radiation of Cu pipes with Kα radiation, 40 mA and 40 kV). The diffraction patterns were obtained at an angular interval of 2θ from 4 to 70°, with steps of 0.05° and 2 s exploration time. The pattern analysis was performed using the Data Collector® software.

Chemical Analysis

The coarse and fine particles underwent a chemical procedure to obtain the bioavailable, reduced, oxidized, and residual fractions. The sequential extraction procedure to obtain each fraction is described by [22,23,36,37]. Subsequently, the concentration of major and trace elements in coarse and fine particles (bulk particles), particle fractions, and dissolved phase was determined by ICP-OES, whereas U and Th isotopes were obtained using alpha spectrometry.

Element Characterization

Aliquots of bulk particles, particle fractions, and dissolved phase were digested using a microwave (ON 60 Hz Anton Paar, Graz, Austria). The element concentrations were obtained using an ICP-OES Optima 8300® (Perkin Elmer®, USA; software version 1.10); a calibration standard of 21 elements (Perkin Elmer, Waltham, MA, USA, LLC STD-AS QC 21 Elements) to calibrate this equipment. For the quality/accuracy control, sediment samples were analyzed by triplicate, and a certified reference material (IAEA-375) was measured; it was treated under the same sequential extraction procedure, fractionation of elements. The element recovery (%) was obtained by comparing the concentration of each element in bulk samples and the sum of concentrations in four fractions. The measurement precision was below 30% of RSD (relative standard deviations), while average recovery rates of element concentrations were between 80 and 110%. In this study, the elements reported were those in which most concentrations were above their respective detection limit.

Radiological Characterization

Radiological characterization was carried out only in suspended sediments at the bottom level. Thus, the contents of U and Th isotopes were determined in bulk particles, particle fractions, and dissolved phase. Internal standards of 232U (117 ± 0.6 mBq/g) and 229Th (104.04 ± 0.6 mBq/g) were added, through which the activities of the isotopes of interest were obtained. Each aliquot underwent a process of digestion, extraction, separation, and purification (deposition) of radioisotopes. The radioisotopes U and Th were obtained through alpha spectrometry with surface barrier detectors (Canberra ®Alpha, Atlanta, GA, USA), and the analysis of alpha spectrums was performed with the Genie 2000® software. The spike recoveries ranged from 84 to 96% and from 55 to 78% for U and Th, respectively.

2.4. Distribution Coefficient (kd)

The distribution of elements between suspended sediments and dissolved phase was obtained by the distribution coefficient, kd [38,39], Equation (1):
k d = Cis   ( mg kg ) Ciw   ( mg L ) ,
The kd is the element concentration ratio, Cis is the concentration of the i-element in the suspended sediments, and Ciw is the concentration of the i-element in the dissolved phase.

2.5. Statistical Analysis

Compositional data analysis (CoDA) methods are necessary due to the compositional nature of geochemical data [40]. Therefore, a Shapiro–Wilk test was employed to analyze the normal distribution of these data. As non-normal distribution was observed, results were log-transformed using the R software (www.rproject.org (accessed on 1 January 2022)) and analyzed by: (a) an ANOVA, to identify the interaction effects of sampling point and water depth in the physicochemical parameters, applying the procedure GLM using the package SAS® [41]; (b) an ANOVA, to identify the interaction effects of sampling point, particle size, fractionation of elements, and water depth in the elemental composition of suspended sediments, applying the procedure GLM using the package SAS®; (c) a distribution of elements using cluster heatmaps (isometric log-ratio transformation) using the R-package robCompositions (pivotCoord and gplots functions); (d) a correlation analysis between the implicit continuous variables, using the R-procedures ggplot2 and corrplot; and (e) a multivariate analysis, Factorial Analysis (FA), using the isometric log-ratio (ilr) transformation methods through the R-function “factanal” [42,43,44].

2.6. Potential Pollution Risk Assessment

Element concentrations from suspended sediments were used to identify their enrichment, pollution levels, and the potential pollution risk [6,7,45,46]. Due to element recoveries being greater than 80%, the analysis of indices was carried out using the sum of element concentrations obtained by fractions (F1–F4). Furthermore, to provide an improved prediction of the potential risk for aquatic biota, the element concentrations from the available sediment fractions (bioavailable, F1; reduced, F2; and oxidized, F3) were used to calculate the Risk Assessment Code (RAC).

2.6.1. Contamination Factor (Cf)

The Cf is a single index used to monitor the element status in sediment [7,11,16]. This index was obtained by Equation (2):
Cf = Ci/Ciearth crust
where Ci is the i-element concentration in sediment and Ciearth crust is the i-element concentration in the upper continental crust. According to [11], Cf is classified as Cf < 1 low contamination factor, 1 ≤ Cf < 3 moderate contamination factor, 3 ≤ Cf < 6 considerable contamination factor, and Cf > 6 very high contamination factor.

2.6.2. Enrichment Factor (Ef)

The Ef was obtained by the normalized ratio from sediment element concentrations to those in the average upper continental crust [47,48]. This normalization takes into consideration variations in sediment grain size and mineralogy. The Ef was calculated by Equation (3):
Ef = (Ci/Re)sediments/(Ci/Re)earth crust
Ci is the i-element concentration in sediment, and Re is the reference element. The reference element must be stable and associated with finer particles; elements such as Al, Fe, Mn, Rb, and Ti are widely used [16,45]. In this work, Fe was used as Re due to is one of the major constituents in clay material [11,49,50,51]. Ef is classified as Ef < 1 deficient/not enriched, 1 < Ef < 3 minor enrichment, 3 < Ef < 5 moderate enrichment, 5 < Ef < 10 moderate-severe enrichment, 10 < Ef < 25 severe enrichment, 25 < Ef < 50 very severe enrichment, and 50 < Ef extremely severe enrichment.

2.6.3. Geoaccumulation Index (Igeo)

The Igeo defines the metal concentration in fine fraction particles and detects the anthropogenic influence [45]. This index is expressed by Equation (4):
Igeo = log2 (Ci/1.5 × Bn)
where Ci is the i-element concentration in sediment, Bn is the geochemical background concentration of the element, and 1.5 is the constant to analyze natural fluctuations in the element contents. Its classification is uncontaminated Igeo < 0, uncontaminated-moderately contaminated 0 < Igeo < 1, moderately contaminated 1 < Igeo < 2, moderately heavily contaminated 2 < Igeo < 3, heavily contaminated 3 < Igeo < 4, heavily extremely contaminated 4 < Igeo < 5, and extremely contaminated 5 < Igeo.

2.6.4. Metal Availability Assessment

Metals in available fractions of sediments are the main responsible for pollution risk. The secondary to primary metal phases ratio (SPR) and the Risk Assessment Code (RAC) were evaluated to obtain the potential pollution risk by releasing metals from suspended sediments.
The SPR was calculated by the ratio of metal concentrations in secondary fractions (F1, F2, and F3) to those in primary fraction (F4) [52]. The SPR and RAC were obtained using the Equations (5) and (6), respectively:
SPR = CF1-3/CF4
where CF1-3i is the sum of i-element in F1 (bioavailable), F2 (reduced), and F3 (oxidized); and CF4i is the concentration of the i-element in residual fraction (F4).
RAC (%) = CF1-2/CTi
where CF1-2 is the sum of i-element concentrations in the most labile fractions (F1 and F2), and CTi is the total concentration of i-element obtained by the sum of their concentrations from all fractions. RAC is classified as no risk < 1, low risk 1–10, medium risk 11–30, high risk 31–50, and very high risk > 50 [6,36].

3. Results

3.1. Bulk Water Quality

The physicochemical parameters and chemical species obtained in bulk samples from the San Marcos reservoir are summarized in Table S1 (Supplementary Material). The pH was neutral at the surface and deep levels of this water body, as well as the ORP showed mainly sub-oxic conditions. In addition, alkalinity at surface level ranged from 14 to 17 mg CaCO3/L, whereas at high depth, 34 to 36 mg CaCO3/L was found. Likewise, the dissolved ions were found in the following abundance SO42− > Na+ > Cl > K+ > NO2 > Ca2+ > Mg2+.

3.2. Mineral Composition of the Suspended Particulate Matter

The mineral analysis of the fine and coarse particles showed that, in general, the main minerals contain Ca, Si, feldspars, and clays. The main mineral species found were Calcite [Ca(CO3)], Quartz [(SiO2)], Albite [(Al1.2Ca0.2Na0.8O8Si2.8)], and Montmorillonite [Ca (Al2Si4O11 (OH))].

3.3. Morphology and Elemental Composition

The results of the particle morphology are shown in Figure 2. The coarse and fine particles are amorphous with a rough and porous surface. Coarse particles show dimensions in a range from 30 to 150 µm. These particles showed differences on their surface: accumulation of small particles (Figure 2a), lightly smooth angles (Figure 2b), and sharp edges (Figure 2c). In addition, fine particles are composed of particle aggregates with dimensions from 13 to 227 µm (Figure 2d–f); it indicates a high surface charge for this particle size, attracting each other to form fine particle agglomerates. Furthermore, coarse particles are mainly composed of Si (35%), Na (23%), and Fe (14%), whereas fine particles are composed of Si and Ca with 41% and 20%, respectively.

3.4. Partitioning of Chemical Species

The concentrations of major and trace elements in suspended particles and dissolved phase are summarized in Table S2 (Supplementary material) at the surface and deep level.

3.4.1. Physicochemical Parameters

From one factor analysis, the concentrations of elements were significantly different according to depth (p < 0.05), but they were not for the sampling point. The alkalinity (CO3=, HCO3), sulfates (SO42−), calcium (Ca2+), and sodium (Na+) showed higher concentrations at the deeper water level, except for Na+. Likewise, the interaction effects depth-sampling point did not show significant differences.
Moreover, some distribution differences of physicochemical parameters related to levels were identified with the heatmaps of correlation coefficients. Figure 3 shows log-transformed parameters sorted according to cluster analysis [53]. The distribution of physicochemical parameters to both levels indicated two main subgroups by clustering (p < 0.05). At the surface level, the NO2, Mg2+, Ca2+, K+, Na+, and ORP form the first subgroup; this can be seen at the lower-left corner of Figure 3a. The highest positive correlations of those parameters were NO2 − Mg2+ and Ca2+ − K+ (R > 0.70), whereas Na+ − ORP showed a weak correlation (R = 0.5). The second subgroup comprises TDS, Alkalinity, pH, SO42−, Cl, T, and hardness (Figure 3a, upper right corner). In the deep level, the subgroup to the upper right corner in Figure 3b showed a high number of parameters that were positively correlated (pH, Alkalinity, Hardness, TDS, T, SO42−, ORP, Mg2+, Na+, and Cl), whereas NO2 and K+ remain clustered as another subgroup (R > 0.62).

3.4.2. Particle Composition

In Table 1 are shown the results of principal components and the kd. Most of the elements are concentrated in suspended sediments showing kd >> 1 Lkg−1. These metals showed kd in the following sequence Ti > Mn > Fe ≈ Ni ≈ V ≈ Cu ≈ Li ≈ Pb ≈ Zn > Mg ≈ Ca ≈ Sr. Statistically, major and trace elements explain >91% of the variance in PC1, showing a high correlation (R > 0.9, p < 0.05) for both particle sizes and at both depth levels. It confirms the natural origin of these elements in suspended sediments.
From ANOVA analysis, the element concentrations did not show significant differences (p < 0.05) according to depth and particle size, see Figure S1 in supplementary material; however, their concentrations were different by particle fractions. Heatmap, based on symmetric coordinates for element concentration by particle fractions, are presented in Figure 4. Uranium and thorium were not included in these analyzes.
The fractionation and co-occurrence of major elements indicate that Ca and Mg tend to be in the bioavailable fraction (F1), Fe and Mn in the reduced fraction (F3), and Ti in the residual fraction (F4), see Figure S1, supplementary material. Besides, trace elements such as Sr, Li, and Zn were mainly concentrated in F1, while the highest contents of Pb were in F2. Finally, Cu, Ni, and V tended to be in the residual fraction.

3.5. Isotopic Determination

U radioisotopes showed higher concentrations in fine particles than in coarse particles. In contrast, Th radioisotope concentrations were similar in both particle sizes (p < 0.05). The concentrations of 238,234U and 232,230Th are shown in Table 2. Likewise, Figure S3 in supplementary material shows the 238,234U and 232,230Th contents by chemical fractions. While the highest mean concentrations of 238U and 234U were found in the bioavailable fraction of fine particles, high concentrations of 230Th and 232Th were found in residual fractions of both particles.

3.6. Potential Pollution Risk

Figure 5 shows the results of the indices that evaluate the status of metals in suspended sediments. Elements such as Cu, Li, Mn, Pb, and Zn showed contamination indices values (Ef, Igeo, and Cf) << 1, mainly in fine particles. These metals reached a very high level of Cf, severe enrichment, and tended to show moderately heavily contaminated levels by Igeo values. Conversely, Ni, Sr, and V showed lower values for these indices.
On the other hand, the SPR ranged from 0.4 to 238 values, where Ca, Mn, and Pb showed the highest SPR-values, see Figure 6. Likewise, the highest RAC percent (>20%) were shown by Mg, Ca, Mn, Cu, Li, Sr, and Zn.

4. Discussion

The availability and quality of surface water are the current challenges to ensure its fair use. Suspended sediments are a comprehensive source of information to obtain the contribution of potentially toxic elements in aquatic ecosystems.
The physicochemical parameters regulate the chemical nature of metals that can participate in chemical reactions and be responsible for the accessibility of some of those chemical species to the biota. ORP showed lower values than those published [28,29], with values at a deeper level from 80 to 88 mV; suggesting a possible induction of reduction conditions at the bottom of this reservoir.
Suspended sediments are, according to grain size [54], in the interval of silts (0.004–0.063 mm) and clay (<0.004 mm). These suspended sediments are composed of calcite [Ca(CO3)], quartz [(SiO2)], albite [(Al1.2Ca0.2Na0.8O8Si2.8)], and montmorillonite [Ca(Al2Si4O11 (OH))]. Several authors have reported that clays (montmorillonite) are the main components of the suspended sediments in a river [14,55]. However, the composition of the sediment material in a reservoir can be attributed to the classification of the local rocks/soils, which under weathering process produce particles of different sizes, including some dissolved chemical species by leaching from these rocks. Previous studies of this reservoir found the same species in the sedimentary material at the bottom level [56]. Chemical partitioning help to understand the changes that metals suffer towards labile forms, allowing the metal availability for the environment. Some metals such as Cu, Ni, Pb, and Zn analyzed in these suspended sediments showed similar contents and distribution by partitioning to those published by other studies [6,52].
The adequate statistical treatment of the data to precisely describe the behavior of chemical species in environmental matrices is a challenge. The concentrations of elements in the geochemical assessment generally present a non-normal behavior. Thus, suspended sediments are a significant source of information about the contents and relationships of elements, which need to be adequately analyzed. According to PCA results, the elements are from natural origin, where depth and particle size do not influence metals’ contents in the sediment matrix (Table 1, Figure S1 supplementary material). Nevertheless, particle size may be a relevant factor of scavenging under ion conditions of elements commonly present in dissolved fractions [1,8,12,19,55,57].
As shown in Figure 4, the concentrations of each element are not modified independently and freely, but were affected by the chemical behavior of the other elements and the environmental characteristics. Due to the elements being concentrated in suspended sediments, their composition was intimately related to the particle fractions. Furthermore, the elements are clustered within every fraction (Figure 4), indicating the co-occurrence of metals. Thus, Ca (>66%), Mg (>53%), Sr (>62), and Li (>32%) were dominantly found in the bioavailable fraction (F1). Moreover, metals such as Ca, Li, Mg, Fe, Cu, and Zn are grouped in this fraction (Figure 4a), where Li, Ca, and Mg are highly correlated (R > 0.92, p value < 0.05). According to the BCR scheme, those metals are linked to carbonates and are easily freed to the aquatic environment [58,59].
Major metals in the earth’s crust tend to be linked to oxides released from the mineral matrix when conditions change from oxic to anoxic. The highest concentrations of Fe (>54%), Mn (>49%), Pb (>60%), and V (>44%) were in the reduced fraction (F2), most of them with a similar oxidation state (2+, 3+). In this fraction of particles, the elements were strongly correlated and distributed in three groups (Figure 4b): (a) Mn, Fe, and V (R > 0.95); (b) Sr, Ca, and Mg (R > 0.90); and (c) Li, Cu, Pb, and Zn (R > 0.90). According to cluster distributions, the two first groups are positively correlated and represent the major elements forming the mineral matrix of suspended sediments. The third group contains the trace elements and is distributed negatively to the first group (a) and weakly correlated to the second group (b); it probably indicates the adsorption of these trace elements onto the mineral matrix. The co-occurrence of Fe, Mn, and V in minerals is well known because of their similar geochemical process on earth’s crust [57]. Likewise, it has been published that Pb concentrations are mainly found in the reduced fraction of bottom sediments [6,22,36]. In this fraction, some trace elements may be bounded to Fe-Mn oxides that would be released under reductive conditions, which suggests that Li, Cu, Pb, and Zn in suspended sediments are being released from Fe/Mn oxides/hydroxides, as in bottom sediments [21,22,60].
No metals were predominant in the oxidized fraction (F3); however, elements such as Li, Mg, Cu, Zn, and Sr are co-occurring in this fraction. Previous studies have demonstrated a high association (>40%) between Cu and organic components [6,52,61], but in this study, copper contents in the oxidized fraction are in a low range (from 14 to 18%). Finally, Ti was concentrated (49–80%) in the residual fraction (F4), indicating its conservative behavior as part of the mineral matrix. Two dominant groups were observed (R > 0.8); the first is formed by Sr, Zn, Ca, Li, and Cu, while the second is grouped by Mg, Mn, Ti, Fe, and Ni. These results coincide with those determined by the XRD and SEM-EDX analyses, where the presence of Ca- and Fe- minerals are observed.
Nevertheless, it is noticeable that considerable contents of Fe (12–40%) and V (36–45%) are also present in the residual fraction, while Mn contents (30–43%) are in the bioavailable fraction. This behavior has been reported previously, where the residual fraction of sediments and soil shown high iron concentrations [6,21,22,62] and, even a high percent of Mn has been detected in the reducible fraction of sediments from city streets [62], this metal is also bounded to carbonates in soil particles [21]. In addition, Zn also tended to be in exchangeable (26–36%) and reducible (24–42%) fractions, as has been reported in other studies, [52] and [6], respectively. Finally, Ti concentrations were also ranged from 20 to 40% in oxidized fraction. Therefore, a significant correlation between bioavailable-oxidized and reducible-residual fractions was observed (p value < 0.05) due to the similarity of metal co-occurrence in those sediment phases; it is shown in Figure 7.
Here, the association and distribution of elements can be observed more clearly due to metal affinity. Thus, elements such as Ca, Cu, Mg, and Zn are grouped in the bioavailable-oxidized fraction (R > 0.85), which are elements classified as nutrients and tend to be present in organic material (Figure 7a). Likewise, the reduced-residual fraction shows three main clusters (R > 0.9; Figure 7b): (a) Ca-Sr, (b) Fe-Mn-Ti, and (c) trace elements. According to these groups, major elements are well defined in two groups, Ca-oxides/carbonates and Fe/Mn-oxyhydroxides, forming the mineral matrix of these sediments. Therefore, it may suggest that suspended sediments are particles produced by the physical and chemical erosion from both sedimentary rocks (limestone, dolomites, among others) and volcanic rocks (intermediate and basic composition) present in the study area [25,26].
On the other hand, lithology is a crucial parameter of high contents of isotopic uranium and thorium concentrations in water [15]. Thus, the radioisotopes 238U, 234U, 232Th, and 230Th were highly correlated (R > 0.8 and p < 0.05). In addition, isotopic uranium was found three-fold more concentrated in fine particles than in coarse particles (p < 0.05), whereas 230Th and 232Th did not show partitioning by size particle (p < 0.05).
The 234U/238U activity ratio is similar in both particles; this disequilibrium can be attributed mainly to the alpha recoil and the particle scavenging. Under oxidizing conditions, U is mainly in its mobile form (hexavalent) [15,37,63,64,65,66,67] and is maintained in solution due to the presence of CO32− and HCO3 [68]. Conversely, Th is usually found in its immovable form (+IV) under these conditions. At a deep level of this reservoir, sub-oxic conditions tend to affect the uranium and thorium concentrations in coarse particles (R > 0.7). Moreover, alkalinity also controls the mobility of isotopic uranium in fine particles (R > −0.7). Likewise, isotopic U and Th increased their concentrations in these suspended sediments by increasing most metal contents, specifically Cu, Fe, Ti, and Zn. Moreover, 230Th/232Th activity ratio was above the unit, indicating thorium from the 238U disintegration chain (alpha recoil).
The radioactive decay of 238U atom in a sediment particle produces alpha particles and daughter atoms, which can be ejected closer to the surface particle by the recoil energy. The alpha-recoil process promotes a radioactive daughter’s mobilization from its initial position. Moreover, uranium might migrate to oxidation states (+6) during the alpha decay, obtaining a high solubility. In suspended sediments, >43% of isotopic uranium was found in the bioavailable fraction, which can be attributed to the alpha-recoil of uranium. Conversely, isotopic thorium was mainly found in the residual fraction of suspended sediments, ranging from 70% to 84% for 232Th and from 49% to 70% for 230Th, showing their natural origin.
The 234U/238U and 230Th/238U activity ratios were in secular equilibrium within the residual fraction of coarse particles. However, this was different in fine particles where 234U/238U and 230Th/238U were 4 and 8, respectively; it corroborates the metal scavenging by finer particles. Furthermore, the 232Th/238U activity ratio ranges from 2 to 8 in the residual fraction, according to the abundance ratio in the earth’s continental crust [47], corroborating the conservative behavior immobile condition for thorium.

Potential Pollution Risk

The degree of metals mobilization in aquatic ecosystems is controlled by their chemical and physical conditions and environmental parameters.
While Ca and Mg can be essential nutrients for living organisms, other metals as Pb and Cu are chemically persistent and can be accumulated by biota, transferring these elements to the food chain [69].
Mobile elements show similar or lesser concentrations in the solid phase than in the solution phase, showing kd values of tens [14], whereas kd values of 100 and 1000 correspond to intermediate and highly immobile elements, respectively. From the total concentrations of elements in suspended sediments, mobile metals were Mg, Ca, Sr, and U; intermediate mobile metals were Fe, V, Ni, Pb, Cu, Li, and Zn, while Ti, Mn, and Th were highly immobile (Table 1).
Due to their chemical affinity, Ca, Mg, and Sr are mainly associated in available phases, Figure 6b. The co-occurrence of those elements is also observed within the sediment matrix or forming oxisalts (Figure 6b), giving them a high mobility status. These mobile divalent metals were correlated with hardness and alkalinity at the surface level, R = 0.9 and R = 0.6, respectively. The water temperature in this level ranged from 16 to 20 °C, leading to the dissociation of Ca/Mg/Sr carbonates from suspended sediments, which, by releasing carbonates, increases the alkalinity of water. Furthermore, when Ca Mg CO3 in water are in the presence of atmospheric CO2 might give a temporary hardness. Meanwhile, on a deeper level, Ca and Mg were highly correlated with chlorides (R = 0.85), providing a total hardness.
On the other hand, although the element mobility can substantially describe the metal behavior in sediments, it is crucial to understand the co-occurrence of these chemical species within solid phases. For instance, Ti and Mn were detected as immobile elements (kd > 1000) and are commonly associated with the mineral matrix; however, considerable amounts of them were obtained in exchangeable phases of suspended sediments; Mn (F2) and Ti (F3). Iron is a chemical reducing agent co-occurring with those elements (Figure 7), leading them to form oxides/hydroxides; thus, Mn and Ti can modify their availability depending on the presence and interrelationship with other elements in solid phases. Therefore, these transition metals (Figure 7b) should be established as elements with intermediate mobility.
Most metals showed concentrations >60% within available fractions of suspended sediments. The metal release due to its availability generates a potential environmental risk [6,36,69,70]. Although metal concentrations in coarse and fine particles are higher than those reported in bottom sediments from surface water bodies and sea [14,60,69,71], these particles showed similarities of element partitioning in sediment fractions [6,7,9,22,36,49,57,58,72].
Due to particle morphology and fractionation of metals, most of them provided contamination factors ranging from moderate to very high levels; even these metals show heavy contamination levels by geoaccumulation, reaching severe enrichment (Figure 5); Zn, Li, Cu, Pb, and Mn contents in fine particles reached extremely severe enrichment. Previously, Pb was found to be moderately enriched in the bottom sediments of this reservoir [56] with a low coefficient of variation on its concentration. Thus, sediment resuspension from the bottom affects Pb contents in suspended sediments, which shows moderate-severe enrichment levels.
Similarly, Igeo values ranged from moderately to heavily contaminated, where the main contributors were Cu, Mn, Li, Pb, and Zn. Commonly, contamination by trace elements is associated with industrial activities and urban development [52,73,74,75,76]. However, in some sediments from areas with high population density and significant industrial activity, pollution by heavy metals was not detected [69].
The Ef and Igeo use gross information of the average element concentrations from the upper earth’s crust and the metal normalization [11,45,49], but it does not consider the intrinsic processes of water bodies and the chemical behavior of elements within solid phases. In an aqueous medium, some new solid phases might be formed when rocks suffer chemical weathering, particularly from carbonate and silicate rocks [77,78,79]. Due to that, these new/secondary solid phases tend to have low weight and small diameters. The finer particles show a high surface area with binding sites that easily trap toxic metals [19]. Morphologically, the coarse particles have well-defined shapes, showing rounded or sharp surfaces (Figure 2a–c); conversely, fine particles showed agglomerates formed by the attraction of charges from their surfaces (Figure 2d–f).
Consequently, some metals can be adsorbed/desorbed from sediments, accordingly their oxidation state. Due to this, these particles can promote the scavenging of these trace metals from the dissolved phase. Therefore, natural processes can also modify the metal concentrations in suspended sediments, turning them into a contamination source for water bodies [80,81].
Metal toxicity is often a function of metal concentration and availability [6,10,58]. In this reservoir, metals are concentrated in the available fractions of suspended sediments above the unit, where Ca, Mn, and Pb are highly concentrated in these secondary phases (SPR); the metals in those phases can only be released under specific conditions, oxidation-reduction conditions in the water. Nevertheless, the metals in bioavailable fraction (F1) are the most accessible metals for biota; some of those are essential nutrients for organisms. Finally, the metals such as Sr (≥50%) > Mn=Li (≥30%) > Cu=Ni (≥20%) are being easily released, and consequently, might be the main contributor to the potential pollution risk for aquatic biota.

5. Conclusions

In this study, the partitioning, association, and geochemical fractionation of metals in suspended sediments were evaluated and used to establish the potential pollution risk of an aquatic ecosystem.
Suspended sediments showed a natural origin based on these results, produced by physical and chemical erosion from the surrounding rocks. However, some metals showed high concentrations linked to the exchangeable phases of suspended sediments.
The sediment phases were the main parameter affecting the partitioning and availability of metals, which are co-occurring due to their chemical affinity; nevertheless, the finer suspended sediments can act as scavengers of metals.
It was determined that metal concentration in sediments is highly dependent on the interrelation and co-occurrence of elements, inducing metals showing low solubility to modify their chemical conditions and migrate to phases that can be released into the surrounding environment. Although the distribution coefficient provides an overview of the mobility of elements, it is crucial to assess their availability within the solid phases. It will provide an accurate assessment of metal behavior and availability.
Pollution indices are a helpful tool for determining the potential pollution risk; however, these should be modified according to the availability of chemical species in the exchangeable phases of sediments. Thus, due to their concentrations and availability, metals such Sr, Mn, Li, Cu, and Ni are the main contributors to potential pollution risk for the aquatic environment under study.
Finally, the adequate statistical treatment of the geochemical data helps obtain a result closer to the occurrence and behavior of metals in the environmental evaluation. Likewise, chemical fractionation methodology improves the prediction of the potential pollution risk caused by metals(oids) in aquatic ecosystems.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w14060980/s1, Table S1: Physiochemical parameters and chemical species in bulk water of the reservoir San Marcos; Table S2: Mean contents of major (% weight) and trace (mg kg−1) elements and their range in suspended sediments and dissolved phase; Figure S1: Association of elements by particle size and depth: (a) surface coarse particles, (b) surface fine particles, (c) deep coarse particles, and (d) deep fine particles; Figure S2: The metal percent by particle fractions: (a) Coarse particles at surface level, (b) Fine particles at the surface, (c) Coarse particles at deep level, and (d) Fine particles at a deep level. F1 (bioavailable fraction); F2 reduced fraction; F3 oxidized fraction; F4 residual fraction; Figure S3: Isotopic U and Th percent by fractions of coarse and fine particles at a deep level. F1 (bioavailable fraction); F2 reduced fraction; F3 oxidized fraction; F4 residual fraction.

Author Contributions

M.C.-L.: conceptualization, methodology, validation, formal analysis. M.R.-V.: conceptualization, methodology, validation, formal analysis, resources, investigation, visualization, supervision, project administration, funding acquisition. A.M.-M.: methodology, formal analysis. Z.O.-C.: methodology, formal analysis. E.M.-C.: conceptualization, methodology, validation. I.V.: conceptualization, methodology. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by CONACyT, in the Proyect CB-2011-01-16697 (Grant number/296368).

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to acknowledge the technical support given by Manuel Reyes Cortés.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Islam, M.S.; Ahmed, M.K.; Raknuzzaman, M.; Habibullah-Al-Mamun, M.; Islam, M.K. Heavy metal pollution in surface water and sediment: A preliminary assessment of an urban river in a developing country. Ecol. Indic. 2015, 48, 282–291. [Google Scholar] [CrossRef]
  2. Toro, M.; Robles, S.; Avilés, J.; Nuño, C.; Vivas, S.; Bonada, N.; Prat, N.; Alba Tercedor, J.; Casas, J.; Guerrero, C.; et al. Calidad de Las Aguas de Los Ríos Mediterráneos del Proyecto GUADALMED; Características Físico-Químicas; Limnetica: Madrid, Spain, 2002; pp. 63–75. [Google Scholar]
  3. Sophocleous, M. Interactions between groundwater and surface water: The state of the science. Hydrogeol. J. 2002, 10, 52–67. [Google Scholar] [CrossRef]
  4. Lu, G.Y.; Wong, D.W. An adaptive inverse-distance weighting spatial interpolation technique. Comput. Geosci. 2008, 34, 1044–1055. [Google Scholar] [CrossRef]
  5. Ruiz-Fernández, A.C.; Páez-Osuna, F.; Urrutia-Fucugauchi, J.; Preda, M. 210Pb geochronology of sediment accumulation rates in Mexico City Metropolitan Zone as recorded at Espejo de los Lirios lake sediments. Catena 2005, 61, 31–48. [Google Scholar] [CrossRef]
  6. Ferrans, L.; Jani, Y.; Burlakovs, J.; Klavins, M.; Hogland, W. Chemical speciation of metals from marine sediments: Assessment of potential pollution risk while dredging, a case study in southern Sweden. Chemosphere 2021, 263, 128105. [Google Scholar] [CrossRef] [PubMed]
  7. Williams, J.A.; Antoine, J. Evaluation of the elemental pollution status of Jamaican surface sediments using enrichment factor, geoaccumulation index, ecological risk and potential ecological risk index. Mar. Pollut. Bull. 2020, 157, 111288. [Google Scholar] [CrossRef] [PubMed]
  8. Gloaguen, T.V.; Motta, P.N.S.D.; Couto, C.F. A grain-size correction for metal pollution indexes in river sediments. Int. J. Sediment Res. 2020, 36, 362–372. [Google Scholar] [CrossRef]
  9. Saravanan, P.; Krishnakumar, S.; Silva, J.D.; Pradhap, D.; Vidyasakar, A.; Radhakrishnan, K.; Godson, P.S.; Arumugam, K.; Magesh, N.S. Elemental concentration and potential ecological risk assessment of reef associated surface sediments of Appa Island, Gulf of Mannar Biosphere Reserve, Southeast coast of India. Mar. Pollut. Bull. 2018, 128, 398–407. [Google Scholar] [CrossRef] [PubMed]
  10. Tomlinson, D.L.; Wilson, J.G.; Harris, C.R.; Jeffrey, D.W. Problems in the assessment of heavy-metal levels in estuaries and the formation of a pollution index. Helgoländer Meeresunters. 1980, 33, 566–575. [Google Scholar] [CrossRef] [Green Version]
  11. Hakanson, L. An ecological risk index for aquatic pollution control. A sedimentological approach. Water Res. 1980, 14, 975–1001. [Google Scholar] [CrossRef]
  12. Rousseau, T.C.C.; Roddaz, M.; Moquet, J.-S.; Handt Delgado, H.; Calves, G.; Bayon, G. Controls on the geochemistry of suspended sediments from large tropical South American rivers (Amazon, Orinoco and Maroni). Chem. Geol. 2019, 522, 38–54. [Google Scholar] [CrossRef]
  13. Oliveira, M.L.S.; Saikia, B.K.; da Boit, K.; Pinto, D.; Tutikian, B.F.; Silva, L.F.O. River dynamics and nanopaticles formation: A comprehensive study on the nanoparticle geochemistry of suspended sediments in the Magdalena River, Caribbean Industrial Area. J. Clean. Prod. 2019, 213, 819–824. [Google Scholar] [CrossRef]
  14. Viers, J.; Dupré, B.; Gaillardet, J. Chemical composition of suspended sediments in World Rivers: New insights from a new database. Sci. Total Environ. 2009, 407, 853–868. [Google Scholar] [CrossRef] [PubMed]
  15. Chabaux, F.; Bourdon, B.; Riotte, J. Chapter 3 U-Series Geochemistry in Weathering Profiles, River Waters and Lakes. In Radioactivity in the Environment; Porcelli, D., Krishnaswami, S., Cochran, J.K., Eds.; Elsevier: Amsterdam, The Netherlands, 2008; Volume 13, pp. 49–104. [Google Scholar]
  16. Looi, L.J.; Aris, A.Z.; Yusoff, F.M.; Isa, N.M.; Haris, H. Application of enrichment factor, geoaccumulation index, and ecological risk index in assessing the elemental pollution status of surface sediments. Environ. Geochem. Health 2019, 41, 27–42. [Google Scholar] [CrossRef]
  17. IAEA. Handbook of Parameter Values for the Prediction of Radionuclide Transfer in Terrestrial and Freshwater Environments; IAEA: Vienna, Austria, 2010; p. 200. [Google Scholar]
  18. Ani, E.-C.; Wallis, S.; Kraslawski, A.; Agachi, P.S. Development, calibration and evaluation of two mathematical models for pollutant transport in a small river. Environ. Model. Softw. 2009, 24, 1139–1152. [Google Scholar] [CrossRef]
  19. Batley, G.E.; Simpson, S. Sediment Quality Guidelines. In Encyclopedia of Aquatic Ecotoxicology; Férard, J.-F., Blaise, C., Eds.; Springer: Dordrecht, The Netherlands, 2013; pp. 1015–1024. [Google Scholar]
  20. Warren, N.; Allan, I.J.; Carter, J.E.; House, W.A.; Parker, A. Pesticides and other micro-organic contaminants in freshwater sedimentary environments—A review. Appl. Geochem. 2003, 18, 159–194. [Google Scholar] [CrossRef]
  21. Alan, M.; Kara, D. Comparison of a new sequential extraction method and the BCR sequential extraction method for mobility assessment of elements around boron mines in Turkey. Talanta 2019, 194, 189–198. [Google Scholar] [CrossRef] [PubMed]
  22. Nemati, K.; Bakar, N.K.A.; Abas, M.R.; Sobhanzadeh, E. Speciation of heavy metals by modified BCR sequential extraction procedure in different depths of sediments from Sungai Buloh, Selangor, Malaysia. J. Hazard. Mater. 2011, 192, 402–410. [Google Scholar] [CrossRef] [PubMed]
  23. Pérez-López, R.; Álvarez-Valero, A.M.; Nieto, J.M. Changes in mobility of toxic elements during the production of phosphoric acid in the fertilizer industry of Huelva (SW Spain) and environmental impact of phosphogypsum wastes. J. Hazard. Mater. 2007, 148, 745–750. [Google Scholar] [CrossRef] [PubMed]
  24. Tessier, A.; Campbell, P.G.C.; Bisson, M. Sequential extraction procedure for the speciation of particulate trace metals. Anal. Chem. 1979, 51, 844–851. [Google Scholar] [CrossRef]
  25. Reyes-Cortés, M.; Reyes-Cortés, I.A.; Espino Valdez, S.; Rentería-Villalobos, M.; Burillo Montúfar, J.C.; Montero-Cabrera, M.E. Origen y distribución de la radiactividad natural en la zona norte de la cuenca de Chihuahua, México. Rev. Mex. Cienc. Geológicas 2012, 29, 659–675. [Google Scholar]
  26. Ferríz, H. Uranium mineralization in the San Marcos volcanic center Chihuahua, Mex. In Proceedings of the Uranium Deposits in Volcanic Rocks, El Paso, TX, USA, 2–5 April 1984; pp. 197–216. [Google Scholar]
  27. INEGI. Compendio de Información Geográfica Municipal; INEGI: Aguascalientes City, Mexico, 2017. [Google Scholar]
  28. Rentería-Villalobos, M.; Cortés, M.R.; Mantero, J.; Manjón, G.; García-Tenorio, R.; Herrera, E.; Montero-Cabrera, M.E. Uranium in the Surrounding of San Marcos-Sacramento River Environment (Chihuahua, Mexico). Sci. World J. 2012, 2012, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Burillo Montufar, J.C.; Reyes Cortés, M.; Reyes Cortés, I.A.; Espino Valdez, M.S.; Hinojosa de la Garza, O.R.; Nevárez Ronquillo, D.P.; Herrera Peraza, E.; Rentería Villalobos, M.; Montero Cabrera, M.E. Uranium-series isotopes transport in surface, vadose and ground waters at San Marcos uranium bearing basin, Chihuahua, Mexico. Appl. Geochem. 2012, 27, 1111–1122. [Google Scholar] [CrossRef]
  30. NOM-014-SSA1-1993. Sanitary Procedures for the sampling of water for human use and consumption in public and private water supply systems. (Procedimientos Sanitarios para el muestreo de agua para uso y consumo humano en sistemas de abastecimiento de agua públicos y privados). Official Mexican Standard. Enviromental Health. Mexico, H.a.A.D., Ed.; DOF, Diario Oficial de la Federación, Secretaria de Economía: Mexico City, Mexico, 1993. [Google Scholar]
  31. NMX-AA-072-SCFI-2001. In Water Analysis—Determination of Total Hardness in Natural, Wastewaters and Wastewaters Treated—Test Method; DOF, Diario Oficial de la Federación, Secretaría de Economía: Mexico City, Mexico, 2001; pp. 1–14.
  32. NMX-AA-073-SCFI-2001. In Water Analysis—Determination of Total Chlorine in Natural Water, Wastewaters and Wastewaters Treated—Test Method; DOF, Diario Oficial de la Federación, Secretaria de Economía: Mexico City, Mexico, 2001; pp. 1–13.
  33. NMX-AA-036-SCFI-2001. In Water Analysis—Determination of Acidity and Total Alkalinity in Natural, Drinking, Wastewaters and Treated Wastewaters; Secretaría de Comercio y Fomento Industrial, Secretaria de Economía: Mexico City, Mexico, 2001; pp. 1–22.
  34. Hach Company. Procedures Manual, DR/2010 Spectrophotometer: Photometric, Titration, and Microbiological Procedures, 3rd ed.; Hach Company: Loveland, CO, USA, 2000; p. 874. [Google Scholar]
  35. NMX-AA-051-SCFI-2001. In Water Analisis—Determination of Metals By Atomic Absorption in Natural, Drinking, Wastewaters and Wastewaters Treated—Test Method; DOF, Diario Oficial de la Federación, Secretaria de Economía: Mexico City, Mexico, 2001; pp. 1–47.
  36. Mondal, P.; Schintu, M.; Marras, B.; Bettoschi, A.; Marrucci, A.; Sarkar, S.K.; Chowdhury, R.; Jonathan, M.P.; Biswas, J.K. Geochemical fractionation and risk assessment of trace elements in sediments from tide-dominated Hooghly (Ganges) River Estuary, India. Chem. Geol. 2020, 532, 119373. [Google Scholar] [CrossRef]
  37. Pérez-Moreno, S.M.; Gázquez, M.J.; Pérez-López, R.; Bolivar, J.P. Validation of the BCR sequential extraction procedure for natural radionuclides. Chemosphere 2018, 198, 397–408. [Google Scholar] [CrossRef]
  38. Degueldre, C.; Kline, A. Study of thorium association and surface precipitation on colloids. Earth Planet. Sci. Lett. 2007, 264, 104–113. [Google Scholar] [CrossRef]
  39. Chabaux, F.; Riotte, J.; Dequincey, O. U-Th-Ra Fractionation during Weathering and River Transport. In Uranium Series Geochemistry: Reviews in Mineralogy and Geochemistry; Bourdon, B., Henderson, G.M., Lundstrom, C.C., Turner, S.P., Eds.; Mineralogical Society of America: Chantilly, VA, USA, 2003; pp. 1–19. [Google Scholar]
  40. Kynčlová, P.; Hron, K.; Filzmoser, P. Correlation Between Compositional Parts Based on Symmetric Balances. Math. Geosci. 2017, 49, 777–796. [Google Scholar] [CrossRef] [Green Version]
  41. SAS. Statistical Analysis System Users’ Guide; Statistical Analysis System Institute, Inc.: Cary, NC, USA, 2006. [Google Scholar]
  42. Boente, C.; Martín-Méndez, I.; Bel-Lán, A.; Gallego, J.R. A novel and synergistic geostatistical approach to identify sources and cores of Potentially Toxic Elements in soils: An application in the region of Cantabria (Northern Spain). J. Geochem. Explor. 2020, 208, 106397. [Google Scholar] [CrossRef]
  43. Tarvainen, T.; Sapon, S.; Jarva, J. Applying heatmaps in interpretation of geochemical baseline data on urban soils in Finland. J. Geochem. Explor. 2019, 205, 106345. [Google Scholar] [CrossRef]
  44. Reimann, C.; Filzmoser, P.; Garrett, R.G. Factor analysis applied to regional geochemical data: Problems and possibilities. Appl. Geochem. 2002, 17, 185–206. [Google Scholar] [CrossRef]
  45. Barbieri, M. The Importance of Enrichment Factor (EF) and Geoaccumulation Index (Igeo) to Evaluate the Soil Contamination. J. Geol. Geophys. 2016, 5, 237. [Google Scholar] [CrossRef]
  46. Sutherland, R.A.; Tolosa, C.A.; Tack, F.M.G.; Verloo, M.G. Characterization of Selected Element Concentrations and Enrichment Ratiosin Background and Anthropogenically Impacted Roadside Areas. Arch. Environ. Contam. Toxicol. 2000, 38, 428–438. [Google Scholar] [CrossRef] [PubMed]
  47. Rudnick, R.L. Earth’s Continental Crust. In Encyclopedia of Geochemistry: A Comprehensive Reference Source on the Chemistry of the Earth; White, W.M., Ed.; Springer International Publishing: Cham, Switzerland, 2018; pp. 392–418. [Google Scholar]
  48. Rudnick, R.L.; Gao, S. Composition of the Continental Crust. Treatise Geochem. 2003, 3, 659. [Google Scholar] [CrossRef]
  49. Birch, G.F. An assessment of aluminum and iron in normalisation and enrichment procedures for environmental assessment of marine sediment. Sci. Total Environ. 2020, 727, 138123. [Google Scholar] [CrossRef]
  50. Lin, Q.; Peng, X.; Liu, B.; Min, F.; Zhang, Y.; Zhou, Q.; Ma, J.; Wu, Z. Aluminum distribution heterogeneity and relationship with nitrogen, phosphorus and humic acid content in the eutrophic lake sediment. Environ. Pollut. 2019, 253, 516–524. [Google Scholar] [CrossRef]
  51. Mora, A.; Jumbo-Flores, D.; González-Merizalde, M.; Bermeo-Flores, S.A. Niveles de Metales Pesados En Sedimentos de La Cuenca Del Río Puyango, Ecuador. Rev. Int. Contam. Ambient. 2016, 32, 13. [Google Scholar] [CrossRef] [Green Version]
  52. Xia, F.; Zhang, C.; Qu, L.; Song, Q.; Ji, X.; Mei, K.; Dahlgren, R.A.; Zhang, M. A comprehensive analysis and source apportionment of metals in riverine sediments of a rural-urban watershed. J. Hazard. Mater. 2020, 381, 121230. [Google Scholar] [CrossRef] [PubMed]
  53. Reimann, C.; Filzmoser, P.; Hron, K.; Kynčlová, P.; Garrett, R.G. A new method for correlation analysis of compositional (environmental) data—A worked example. Sci. Total Environ. 2017, 607–608, 965–971. [Google Scholar] [CrossRef] [PubMed]
  54. Franz, G.; Leitão, P.; Pinto, L.; Jauch, E.; Fernandes, L.; Neves, R. Development and validation of a morphological model for multiple sediment classes. Int. J. Sediment Res. 2017, 32, 585–596. [Google Scholar] [CrossRef]
  55. Tian, S.; Li, Z.; Wang, Z.; Jiang, E.; Wang, W.; Sun, M. Mineral composition and particle size distribution of river sediment and loess in the middle and lower Yellow River. Int. J. Sediment Res. 2020, 36, 392–400. [Google Scholar] [CrossRef]
  56. Méndez-García, C.G. Reconstrucción Histórica de la Contaminación Por Arsénico, Plomo y Uranio en Los Sedimentos de Las Represas San Marcos y Luis L. León, en Chihuahua; Original Research; Centro de Investigación en Materiales Avanzados, CIMAV: Chihuahua, Mexico, 2014; p. 148. [Google Scholar]
  57. Telfeyan, K.; Breaux, A.; Kim, J.; Cable, J.E.; Kolker, A.S.; Grimm, D.A.; Johannesson, K.H. Arsenic, vanadium, iron, and manganese biogeochemistry in a deltaic wetland, southern Louisiana, USA. Mar. Chem. 2017, 192, 32–48. [Google Scholar] [CrossRef]
  58. Soliman, N.F.; El Zokm, G.M.; Okbah, M.A. Risk assessment and chemical fractionation of selected elements in surface sediments from Lake Qarun, Egypt using modified BCR technique. Chemosphere 2018, 191, 262–271. [Google Scholar] [CrossRef]
  59. Fernández-Ondoño, E.; Bacchetta, G.; Lallena, A.M.; Navarro, F.B.; Ortiz, I.; Jiménez, M.N. Use of BCR sequential extraction procedures for soils and plant metal transfer predictions in contaminated mine tailings in Sardinia. J. Geochem. Explor. 2017, 172, 133–141. [Google Scholar] [CrossRef]
  60. Yuan, C.G.; Shi, J.B.; He, B.; Liu, J.F.; Liang, L.N.; Jiang, G.B. Speciation of heavy metals in marine sediments from the East China Sea by ICP-MS with sequential extraction. Environ. Int. 2004, 30, 769–783. [Google Scholar] [CrossRef] [PubMed]
  61. Bruder-Hubscher, V.; Lagarde, F.; Leroy, M.J.F.; Coughanowr, C.; Enguehard, F. Application of a sequential extraction procedure to study the release of elements from municipal solid waste incineration bottom ash. Anal. Chim. Acta 2002, 451, 285–295. [Google Scholar] [CrossRef]
  62. Kartal, Ş.; Aydın, Z.; Tokalıoğlu, Ş. Fractionation of metals in street sediment samples by using the BCR sequential extraction procedure and multivariate statistical elucidation of the data. J. Hazard. Mater. 2006, 132, 80–89. [Google Scholar] [CrossRef]
  63. Riedel, T.; Kübeck, C. Uranium in groundwater–a synopsis based on a large hydrogeochemical data set. Water Res. 2018, 129, 29–38. [Google Scholar] [CrossRef] [PubMed]
  64. Priestley, S.C.; Payne, T.E.; Harrison, J.J.; Post, V.E.A.; Shand, P.; Love, A.J.; Wohling, D.L. Use of U-isotopes in exploring groundwater flow and inter-aquifer leakage in the south-western margin of the Great Artesian Basin and Arckaringa Basin, central Australia. Appl. Geochem. 2018, 98, 331–344. [Google Scholar] [CrossRef]
  65. Noli, F.; Kazakis, N.; Vargemezis, G.; Ioannidou, A. The uranium isotopes in the characterization of groundwater in the Thermi-Vasilika region, northern Greece. Isot. Env. Health Stud. 2016, 52, 405–413. [Google Scholar] [CrossRef] [PubMed]
  66. Yamamoto, M.; Tomita, J.; Sakaguchi, A.; Ohtsuka, Y.; Hoshi, M.; Apsalikov, K. Uranium isotopes in well water samples as drinking sources in some settlements around the Semipalatinsk Nuclear Test Site, Kazakhstan. J. Radioanal. Nucl. Chem. 2010, 284, 309–314. [Google Scholar] [CrossRef] [Green Version]
  67. Reyes, E.; Marques, L.S. Uranium series disequilibria in ground waters from a fractured bedrock aquifer (Morungaba Granitoids–Southern Brazil): Implications to the hydrochemical behavior of dissolved U and Ra. Appl. Radiat. Isot. 2008, 66, 1531–1542. [Google Scholar] [CrossRef] [PubMed]
  68. Osmond, J.K.; Ivanovich, M. Uranium-series mobilization and surface hydrology. In Uranium-Series Disequilibrium: Applications to Earth, Marine, and Environmental Science; Ivanovich, M., Harmon, R.S., Eds.; Clarendon Press: Oxford, UK, 1992; pp. 259–289. [Google Scholar]
  69. Fan, J.; Jian, X.; Shang, F.; Zhang, W.; Zhang, S.; Fu, H. Underestimated heavy metal pollution of the Minjiang River, SE China: Evidence from spatial and seasonal monitoring of suspended-load sediments. Sci. Total Environ. 2020, 760, 142586. [Google Scholar] [CrossRef] [PubMed]
  70. Yu, Z.; Liu, E.; Lin, Q.; Zhang, E.; Yang, F.; Wei, C.; Shen, J. Comprehensive assessment of heavy metal pollution and ecological risk in lake sediment by combining total concentration and chemical partitioning. Environ. Pollut. 2021, 269, 116212. [Google Scholar] [CrossRef] [PubMed]
  71. Burton, J.G.A. Sediment quality criteria in use around the world. Limnology 2002, 3, 65–76. [Google Scholar] [CrossRef]
  72. Jorfi, S.; Maleki, R.; Jaafarzadeh, N.; Ahmadi, M. Pollution load index for heavy metals in Mian-Ab plain soil, Khuzestan, Iran. Data Brief 2017, 15, 584–590. [Google Scholar] [CrossRef]
  73. Guan, Q.; Cai, A.; Wang, F.; Wang, L.; Wu, T.; Pan, B.; Song, N.; Li, F.; Lu, M. Heavy metals in the riverbed surface sediment of the Yellow River, China. Environ. Sci. Pollut. Res. 2016, 23, 24768–24780. [Google Scholar] [CrossRef]
  74. Aiman, U.; Mahmood, A.; Waheed, S.; Malik, R.N. Enrichment, geo-accumulation and risk surveillance of toxic metals for different environmental compartments from Mehmood Booti dumping site, Lahore city, Pakistan. Chemosphere 2016, 144, 2229–2237. [Google Scholar] [CrossRef] [PubMed]
  75. Zahra, A.; Hashmi, M.Z.; Malik, R.N.; Ahmed, Z. Enrichment and geo-accumulation of heavy metals and risk assessment of sediments of the Kurang Nallah—Feeding tributary of the Rawal Lake Reservoir, Pakistan. Sci. Total Environ. 2014, 470–471, 925–933. [Google Scholar] [CrossRef] [PubMed]
  76. Varol, M. Assessment of heavy metal contamination in sediments of the Tigris River (Turkey) using pollution indices and multivariate statistical techniques. J. Hazard. Mater. 2011, 195, 355–364. [Google Scholar] [CrossRef]
  77. Gaillardet, J.; Viers, J.; Dupré, B. 7.7—Trace Elements in River Waters. In Treatise on Geochemistry, 2nd ed.; Holland, H.D., Turekian, K.K., Eds.; Elsevier: Oxford, UK, 2014; pp. 195–235. [Google Scholar]
  78. Gaillardet, J.; Dupré, B.; Allègre, C.J. Geochemistry of large river suspended sediments: Silicate weathering or recycling tracer? Geochim. Cosmochim. Acta 1999, 63, 4037–4051. [Google Scholar] [CrossRef]
  79. Sun, H.; Han, J.; Li, D.; Zhang, S.; Lu, X. Chemical weathering inferred from riverine water chemistry in the lower Xijiang basin, South China. Sci. Total Environ. 2010, 408, 4749–4760. [Google Scholar] [CrossRef] [PubMed]
  80. Zhang, Y.; Li, H.; Yin, J.; Zhu, L. Risk assessment for sediment associated heavy metals using sediment quality guidelines modified by sediment properties. Environ. Pollut. 2021, 275, 115844. [Google Scholar] [CrossRef] [PubMed]
  81. Almeida, L.C.; da Silva Júnior, J.B.; dos Santos, I.F.; de Carvalho, V.S.; de Santana Santos, A.; Hadlich, G.M.; Ferreira, S.L.C. Assessment of toxicity of metals in river sediments for human supply: Distribution, evaluation of pollution and sources identification. Mar. Pollut. Bull. 2020, 158, 111423. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Satellite image of the study area. Distribution of the sampling points, from 1 to 9 along the flux water, in the San Marcos reservoir.
Figure 1. Satellite image of the study area. Distribution of the sampling points, from 1 to 9 along the flux water, in the San Marcos reservoir.
Water 14 00980 g001
Figure 2. Morphological analysis with secondary electrons of the suspended particulate matter. Coarse particle dimensions: (a) 130 × 150 µm, (b) 30 × 35 µm, and (c) 50 × 60 µm. Fine particles dimensions: (d) 160 × 227 µm, (e) 51 × 84 µm, and (f) 13 × 17 µm.
Figure 2. Morphological analysis with secondary electrons of the suspended particulate matter. Coarse particle dimensions: (a) 130 × 150 µm, (b) 30 × 35 µm, and (c) 50 × 60 µm. Fine particles dimensions: (d) 160 × 227 µm, (e) 51 × 84 µm, and (f) 13 × 17 µm.
Water 14 00980 g002aWater 14 00980 g002b
Figure 3. Heatmaps of physicochemical parameters for: (a) surface level and (b) deep level. The strength of color indicates the strength of correlation. Blue colors indicate positive correlations, light colors (yellow/orange) indicate negative correlations.
Figure 3. Heatmaps of physicochemical parameters for: (a) surface level and (b) deep level. The strength of color indicates the strength of correlation. Blue colors indicate positive correlations, light colors (yellow/orange) indicate negative correlations.
Water 14 00980 g003
Figure 4. Fractionation of metal: (a) bioavailable fraction (F1); (b) reduced fraction (F2); (c) oxidized fraction (F3); and (d) residual fraction (F4).
Figure 4. Fractionation of metal: (a) bioavailable fraction (F1); (b) reduced fraction (F2); (c) oxidized fraction (F3); and (d) residual fraction (F4).
Water 14 00980 g004aWater 14 00980 g004b
Figure 5. Results of Cf, Igeo, and Ef produced by metals in suspended sediments.
Figure 5. Results of Cf, Igeo, and Ef produced by metals in suspended sediments.
Water 14 00980 g005
Figure 6. Secondary–primary phases ratio (SPR) and RAC values.
Figure 6. Secondary–primary phases ratio (SPR) and RAC values.
Water 14 00980 g006
Figure 7. Co-occurrence of metals within integrated sediment fractions: (a) bioavailable-oxidized fraction and (b) reduced-residual fraction. Distinctive groups are highlighted by circles.
Figure 7. Co-occurrence of metals within integrated sediment fractions: (a) bioavailable-oxidized fraction and (b) reduced-residual fraction. Distinctive groups are highlighted by circles.
Water 14 00980 g007
Table 1. Principal component analysis and distribution coefficient (kd; L/kg) of the suspended sediment samples.
Table 1. Principal component analysis and distribution coefficient (kd; L/kg) of the suspended sediment samples.
VariableComponentskd
Coarse ParticlesFine ParticlesCoarse ParticlesFine Particles
PC1PC2PC1PC2
Surface level
Fe0.284−0.0060.3020.059.0 × 1027.4 × 102
Mg0.284−0.0340.301−0.0021.3 × 1012.0 × 101
Ca0.284−0.0160.3020.0389.9 × 1002.6 × 101
Mn0.2840.0110.2960.0981.1 × 1045.1 × 103
Pb0.2790.2530.302−0.0525.2 × 1011.7 × 102
V0.2750.0720.278−0.1091.7 × 1022.4 × 102
Ni0.207−0.937−0.0420.732.5 × 1022.3 × 102
Ti0.2830.0970.3020.045--
Sr0.284−0.0220.3010.091.0 × 1013.3 × 101
Cu0.2840.0620.3020.0597.1 × 1013.0 × 102
Li0.2810.1820.3020.0456.8 × 1013.1 × 102
Zn0.2840.0290.3010.0884.3 × 1011.5 × 102
U *------
Autocorrelation12.360.5310.891.73
Variance (%)95.14.183.813.3
Accumulated variance (%)95.199.283.897.1
Deep level
Fe0.2890.0860.2810.1412.6 × 1023.9 × 103
Mg0.2890.0760.281−0.1456.8 × 1009.4 × 101
Ca0.2870.1450.28105.7 × 1003.7 × 101
Mn0.2880.1220.2790.2487.7 × 1031.0 × 105
Pb0.279−0.1630.276−0.2681.4 × 1011.3 × 102
V0.135−0.9240.266−0.6555.1 × 1018.2 × 102
Ni0.283−0.0170.270.0491.5 × 1025.5 × 102
Ti0.2870.1140.2720.512
Sr0.2890.0530.281−0.1282.9 × 1002.0 × 101
Cu0.2890.0530.2790.2253.5 × 1011.3 × 102
Li0.2880.0530.2810.0792.5 × 1011.2 × 102
Zn0.2880.0960.2810.133.2 × 1011.3 × 102
U * 8.7 × 1002.3 × 101
Autocorrelation11.890.9112.610.25
Variance (%)91.57971.9
Accumulated variance (%)91.598.59798.9
* It was calculated as total uranium, sum of 238U and 234U concentrations.
Table 2. Partitioning of 238,234U and 232,230Th, mBq kg−1, and their activity ratios, in fractions of suspended sediments at high depth.
Table 2. Partitioning of 238,234U and 232,230Th, mBq kg−1, and their activity ratios, in fractions of suspended sediments at high depth.
CoarseCoarsesum *FineFinesum *Dissolved 1
238U45651251975.3
26–5716–9090–15278–4822–9.7
234U8711114232030
44–15930–171100–189175–47110–70
232Thnm a744979nm
18–10431–7443–133
230Thnm6416978nm
7–92101–24643–160
234U/238U2.11.71.22.16.3
1.2–3.01.5–1.90.8–1.51–3.51.5–10.5
230Th/232Th-0.83.61.0-
0.4–1.02.3–5.00.8–1.2
232Th/238U-1.10.40.6-
1.0–1.30.2–0.80.1–1.0
230Th/238U-0.91.50.6-
0.4–1.20.7–2.70.1–1.0
* The total concentration was obtained from the average concentration in every fraction and averaging the mean from each sampling point. 1 The concentration of elements is given in mBq L−1. a No measured.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cabral-Lares, M.; Rentería-Villalobos, M.; Mendieta-Mendoza, A.; Ortíz-Caballero, Z.; Montero-Cabrera, E.; Vioque, I. Partitioning and Availability of Metals from Water Suspended Sediments: Potential Pollution Risk Assessment. Water 2022, 14, 980. https://doi.org/10.3390/w14060980

AMA Style

Cabral-Lares M, Rentería-Villalobos M, Mendieta-Mendoza A, Ortíz-Caballero Z, Montero-Cabrera E, Vioque I. Partitioning and Availability of Metals from Water Suspended Sediments: Potential Pollution Risk Assessment. Water. 2022; 14(6):980. https://doi.org/10.3390/w14060980

Chicago/Turabian Style

Cabral-Lares, Magaly, Marusia Rentería-Villalobos, Aurora Mendieta-Mendoza, Ziury Ortíz-Caballero, Elena Montero-Cabrera, and Ignacio Vioque. 2022. "Partitioning and Availability of Metals from Water Suspended Sediments: Potential Pollution Risk Assessment" Water 14, no. 6: 980. https://doi.org/10.3390/w14060980

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop