Next Article in Journal
Behavioral Impairment in Aquatic Organisms Exposed to Neurotoxic Pollutants
Next Article in Special Issue
Biofilm and Rivers: The Natural Association to Reduce Metals in Waters
Previous Article in Journal
COVID-19 Lockdowns—Effect on Concentration of Pharmaceuticals and Illicit Drugs in Two Major Croatian Rivers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of the Anthropogenic Impact and Distribution of Potentially Toxic and Rare Earth Elements in Lake Sediments from North-Eastern Romania

by
Laurentiu Valentin Soroaga
1,2,3,†,
Cornelia Amarandei
1,2,3,†,
Alina Giorgiana Negru
2,3,
Romeo Iulian Olariu
1,2,* and
Cecilia Arsene
1,2,*
1
Faculty of Chemistry, “Alexandru Ioan Cuza” University of Iasi, 11 Carol I, 700506 Iasi, Romania
2
ICI-CERNESIM, “Alexandru Ioan Cuza” University of Iasi, 11 Carol I, 700506 Iasi, Romania
3
RECENT-AIR, “Alexandru Ioan Cuza” University of Iasi, 11 Carol I, 700506 Iasi, Romania
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 15 April 2022 / Revised: 5 May 2022 / Accepted: 7 May 2022 / Published: 10 May 2022
(This article belongs to the Special Issue Metal Oxidative Stress in Polluted Inland Water)

Abstract

:
Chemical analysis was performed on sediment samples collected in two sampling sessions (July and October) from Podu Iloaiei Dam Lake, one of the most important water resources used for aquaculture in north-eastern Romania. The concentration of 15 trace elements (TEs), 8 refractory elements (REs), and 15 rare earth elements (REEs)—determined using inductively coupled plasma mass spectrometry—showed variability largely dependent of the sampling points and collection time. Manganese was the most abundant TE, V and Cr were the most abundant REs, while Ce was one of the most abundant REEs. The cerium negative anomaly and Gd positive anomaly were observed in the Chondrite-normalized distributions. In October, the Ce anomaly showed significant negative correlation with Mn, emphasizing the water body oxidation potential. The identified positive Gd anomaly was most likely associated with the use of Gd-chelating agents in magnetic resonance imaging in Iasi, the largest medical hub in north-eastern Romania. Principal component analysis extracted three factors explaining 96.0% of the observed variance, i.e., rock weathering, leaching from soil surface, contributions from urban stormwater and atmospheric deposition (50.9%), pedological contributions (23.7%), and mixed anthropogenic sources (e.g., traffic, waste discharge, agricultural activities; 21.4%). The evaluation of pollution indices highlighted low and moderate degrees of contamination for most of the elements and a considerable degree of contamination for Cd. Assigned Cd sources included fertilizers and pesticides used in the near agricultural areas or the high traffic road located near the lake. Since contamination of aquatic ecosystems with harmful elements is a human health concern, further monitoring of specific vectors in the food chain of the investigated dam lake will be of the utmost importance.

1. Introduction

The nature of element sources in sediments, combined with their different mobility, persistence, and toxicity, are critical factors in the assessment of their environmental impact [1]. Determining the concentration and distribution of elements provides an accessible method to evaluate anthropogenic impact on lake sediments and may represent an extent of biological stress caused by contaminants [2].
In urbanized areas, several heavy metals at trace levels are assigned as major pollutants. The direct discharge of industrial effluents, domestic wastewater into the environment, and use of agricultural chemicals has led to elevated heavy metal concentrations in urban soils and surface water bodies [3,4,5]. In lake sediments, the accumulation of heavy metals may negatively impact the status of the biota. Heavy metal accumulation in fish tissue through the food chain may pose a health risk to humans [6,7].
The extensive use of relatively new technologies is a major source of rare earth elements (REEs) in the environment, and the emission of REEs into natural waters is increasing [8]. The most important input of REEs in the environment is comprised by the fertilizers used in agriculture [9,10], intensive use in the industry, specific medical applications [9,11], and dry deposition of aerosol particles [8]. The REEs are characterized by very similar physicochemical behaviors [12,13], are very good tracers of chemical [14,15,16] and geochemical [9,16,17] processes, and may represent a new class of micropollutants that are not fully understood and could impact the environment.
Although, in aqueous environments, the REEs chemistry involves almost exclusively the trivalent oxidation state, stable compounds formed by Ce4+ and Eu2+ are of high interest [13,18], as they may fractionate from the trivalent state REEs, according to the water redox potential [12,13,19]. In well-oxygenated aqueous environments, bacterial activity mediated oxidation of Ce3+ to Ce4+ may occur [20], with the tetravalent state associated mainly with the water insoluble solid phase [19]. In this case, some sort of Ce depletion is observed in the aqueous phase, leading to a different distribution than the normal crust abundances of REEs, a process known as negative Ce anomaly [12,18]. Both the positive and negative Ce anomaly (Ce/Ce*) can be used for the evaluation of the redox conditions of aquatic environments [12,13,21]. Gadolinium is another REE susceptible to anomalies, and its positive anomaly (Gd/Gd*) is most likely associated with anthropogenic sources related to the use of Gd-chelating agents in magnetic resonance imaging (MRI) [11,22,23].
The spatial distribution characteristics of metals in sediment may indirectly give an indication of the pollution sources localization. Multivariate statistics techniques (e.g., principal component analysis, correlation analysis) are powerful tools for pollution source assessments [4,24]. These techniques were applied for the investigation of pollutants in soil, sediment [25,26,27,28], and natural waters [29,30].
Sediment contamination across a wide range of environments is reported in several studies [1,2,3,31]. Contamination with various elements is often assessed using approaches such as contamination factor (CF), pollution load index (PLI), geoaccumulation index (Igeo), and enrichment factor (EF) [32,33,34,35]. While PLI standardizes the contribution from the elements and is used to indicate bioavailability [33], EF allows to identify natural variations of a given element in order to detect even a very small anthropogenic impact and to find out the extent of the element contaminants to be evaluated in toxicological studies [36,37]. In the estimation of the EF, geochemical elements present in large quantities (e.g., Mn) [1,34] are usually used as reference elements to determine anthropogenic enrichment. Moreover, freshwater sediment quality guidelines (SQGs) are useful in assessing the sediment quality and the potential of toxic elements to cause negative effects on aquatic organisms [38,39]. Consensus-based SQGs, namely the Threshold Effect Concentration (TEC) and Probable Effect Concentration (PEC) [1,3,39], can be used for assessing element effects on the environment. No adverse effects of the contaminant on organisms are expected to occur at concentration levels below TEC. Potential toxic effects may or may not occur for values greater than or equal to TEC but lower than PEC. At concentrations above PEC, probable toxic effects of contaminants would be very likely to occur [39,40].
To our knowledge, limited information on lake sediment quality in north-eastern Romania is available. Strungaru et al. [6] studied the potential contamination with Cu, Cd, Pb, Cr, and Ni of water and sediment in the Stanca-Costesti freshwater reservoir from north-eastern Romania. Although the information is valuable, the number of investigated contaminants is scarce and the potential sources were not evaluated. Sandu et al. [41] focused on the REEs and some trace element geology and geochemistry evaluation in river sediments in the Bistricioara hydrographic basin (Eastern Carpathian Mountains, Romania). In this case, the investigated remote location is less likely to be influenced by anthropogenic urban contributions.
The present study aimed both to evaluate the distribution of potentially toxic heavy metals and rare earth elements and to assess the anthropogenic impact on lake sediments from Podu Iloaiei Dam Lake. The investigated lake is located on the Bahluet River from the Bahlui river basin and is considered as one of the most anthropized hydrological basins of the country in north-eastern Romania. The evaluation of possible contamination of the lake water and sediment is of the utmost importance since it is an important water resource for fish farming in the region. The objectives of this study were: (1) to investigate the spatial distribution characteristics of 15 trace elements (TEs), 8 refractory elements (REs), and 15 rare earth elements (REEs) in Podu Iloaiei Dam Lake sediments in north-eastern Romania; (2) to evaluate element sources in sediments and driving environmental processes using multivariate statistical methods; (3) to assess the pollution and to identify the harmful potential of these elements in the investigated area. The results obtained in this study may provide a reliable scientific database for the estimation of the local ecological status and implementation of pollution remediation strategies.

2. Materials and Methods

2.1. Site Description and Sample Collection

The Podu Iloaiei Dam Lake is located on the Bahluet River from the Bahlui hydrographic basin, north-eastern Romania. More geographical details of the hydrographic area are elsewhere reported [42]. The surroundings of the lake and the sampling points are presented in Figure 1. In general, the water level of Podu Iloaiei Dam Lake (water surface of 251 ha) fluctuates considerably throughout the year, with lower values during the summer (June, July, and August) and higher values during major rainfall events in the autumn (September, October, and November), with an average annual rainfall of 589 mm. In 2017 (the year of sampling), the highest average precipitation values were recorded in May (44 mm) and June (54 mm). The lowest level of precipitation was recorded in October, with an average of 27 mm. The Podu Iloaiei Dam Lake is mainly used for fish farming and is also an important source of water for fish ponds downstream of the dam [43]. The water is drained for the wintertime.
In the Prut–Barlad hydrographic area, including the Bahlui hydrographic basin, agricultural (71%), forest (14%), and residential (12%) areas are the dominant forms of land use. According to the estimate of the National Institute of Statistics, only a small area is occupied by surface waters (3%) [44].
A total number of 28 sediment samples were collected from Podu Iloaiei Dam Lake within s two session measurement, on 16 July 2017 and 5 October 2017, and the distribution of the sediment sampling sites is presented in Figure 1. Data on 88 water samples collected during the same sampling sessions are presented elsewhere [42].

2.2. Sample Preparation of the Sediment

The solid-state dissolution procedure and inductive coupled plasma-mass spectrometry (ICP-MS) analyses were applied for the metal elemental analysis in the sediment samples. Sodium peroxide (Na2O2) was used as the oxidizing agent, as it is one of the strongest available [45,46]. A ratio of 1 to 6.5 for sample to Na2O2 reagent proved to be suitable for complete sample dissolution.
The sediment samples were dried slowly and, subsequently, were finely grinded using an agate mortar and pestle. Around 0.1 g of sediment sample was weighed, and 50 µL of internal standard solution (Indium—200 µg g−1) was added to each sample. The samples were further mixed with the oxidizing agent, followed by transfer to the porcelain crucibles lined with aluminum foil. Thermal processing was performed at 460 °C. After cooling, the sintering mix was dissolved in HNO3 3% solution prepared in deionized water and stored in closed polypropylene vials until analysis.

2.3. Chemical Analysis

The concentrations of trace elements (TEs: Li, Be, Mn, Co, Zn, Ga, As, Rb, Zr, Cd, Ba, Tl, Pb, Bi, and U), refractory elements (REs: V, Cr, Nb, Mo, Hf, Ta, W, and Ir), and rare earth elements (REEs, divided in Light REEs—LREEs: La, Ce, Pr, Nd, Sm, and Eu; and Heavy REEs—HREEs: Y, Gd, Tb, Dy, Ho, Er, Tm, Yb, and Lu) in the sediment samples were determined by the ICP-MS technique using the ICP-MS 7700X system (Agilent Technologies, Santa Clara, CA, USA).
Analyte quantification was performed by an internal standard calibration procedure. Certipur® Indium ICP standard from Merck KGaA (Darmstadt, Germany) was used as an internal standard. Certipur® ICP multi-element standard solution XXI for MS (Merck KGaA, Darmstadt, Germany) was used for TE quantification, Specpure® refractory metals plasma standard solution (Alfa Aesar—Thermo Fisher GmbH, Erlenbachweg, Kandel, Germany) was used for RE quantification, and TraceCERT® rare earth element mix for ICP (Sigma-Aldrich—Merck KGaA, Darmstadt, Germany) was used for REE quantification. Nitric acid 65% Suprapur® (Merck KGaA, Darmstadt, Germany) was used for HNO3 3% preparation. Ultrapure water obtained from an Advantage A10 water purification system (Milli-Q—Merck-Millipore, Darmstadt, Germany) was used for standards and sample dilution.
Special emphasis was placed on carefully controlling the entire array of preparatory and analysis steps in order to prevent the possible contributions brought by sample loss or contamination. For quality control assurance, blanks and control samples were frequently analyzed.

2.4. Statistical Analysis

2.4.1. Statistical Approaches for Source Assessment

The geostatistical approach was used in the present work to comprehensively describe the variation of elements concentration in the sediment of Podu Iloaiei Dam Lake. For the quantified elements, analysis of the geospatial distribution was used in order to discriminate between potential contributions of natural or anthropogenic sources. Google Earth Pro and ArcMap (included in the ArcGIS 10.2 package) programs were used for geospatial distribution analyses. For data processing inverse distance weighting (IDW), interpolation was used.
Pearson correlation analysis and principal component analysis (PCA) were used to identify the possible sources of elements in sediment samples from Podu Iloaiei Dam Lake. Varimax normalized PCA was used to determine the association of the different variables by reducing the dimensionality of the dataset. The concentrations of the quantified elements from sediment samples were used as variables for the PCA. The validity of the PCA was checked with the Kaiser-Meyer-Olkin (KMO) statistic and Bartlett’s test using SPSS (version 22.0). The PCA was performed to determine possible chemical associations that could explain TE, RE, and REE behaviors and the relationships between them.
The Kolmogorov-Smirnov test, carried out in OriginLabPro 2022, was performed to check for possible differences between the sampling sessions. Differences were considered of significance at the level of p < 0.05.
The correlation matrix, including Pearson correlation coefficient values, pie charts, pollution indices heat map representations, and the boxplots of the potentially harmful elements, were graphically represented in OriginLabPro 2022.

2.4.2. REEs Normalization and Anomalies

The normalization to Chondrite [17] and the upper continental crust (UCC) [47] was performed for the evaluation of the REE distribution in the investigated sediment samples and for a better comparison of the two sampling sessions. Although North American Shell Composites (NASC), Post-Archean Australian Shale (PAAS), and European Shell (ES) normalizations were also performed, in terms of distributions, the differences compared to UCC were almost insignificant and therefore were of no interest in the present study.
The Ce and Gd anomalies were estimated as identified in the literature with Equation (1) for Ce/Ce* [12,21,48,49] and Equation (2) for Gd/Gd* [22,48,50]:
Ce / Ce *   =   2   Ce N / ( L a N + P r N )
Gd / Gd *   =   3   Gd N / ( S m N + 2   T b N )
where N represents the normalized concentrations to the UCC values [47].
The bell-shape index (BSI) was evaluated using the formula presented in Equation (3) [12], with the main aim to determine the enrichment of the environment in medium REEs. The lREE, mREE, and hREE factors were determined using Equations (4)–(6), respectively. The concentration of REEs was normalized to UCC, and Ce was excluded from the lREE factor calculation due to anomalous behavior.
B S I   =   2   m R E E s / ( l R E E s + h R E E s )
l R E E s   =   ( L a N + P r N + N d N ) / 3
m R E E s   =   ( S m N + E u N + G d N + T b N + D y N ) / 5
h R E E s   =   ( H o N + E r N + T m N + Y b N + L u N ) / 5

2.4.3. Sediment Contamination Assessment

The contamination factor (CF), pollution load index (PLI), geoaccumulation index (Igeo), and enrichment factor (EF) were estimated in the present work in order to assess the Podu Iloaiei Dam Lake sediment contamination with TEs and REs. The equation and the criteria for assessing the pollution indices are presented in Table 1. In the estimation of the EF, Mn was used as a reference element to determine the anthropogenic element enrichment.
Because background element concentrations in the Podu Iloaiei Dam Lake sediment are not available, the Average Shale Values (ASVs) from [51] were used in the present work to calculate the pollution indices.
The potential of the contaminants to cause adverse effects on aquatic organisms was determined by comparing the concentrations of Cr, Zn, As, Pb, and Cd with the Sediment Quality Guideline (SQG) values for freshwater sediments.

3. Results and Discussion

3.1. Elements Concentration in Sediment Samples

Table 2 shows the mean, minimum, and maximum concentration (μg g−1) values for the TEs, REs, and REEs determined in the sediment samples collected in July and October 2017 from Podu Iloaiei Dam Lake. REE concentration values in the water and the water-soluble fraction from sediments (all below the limits of quantification) are presented in Table S1 in the Supplementary Materials (SM).
In the sediment samples, the average of ∑TE concentrations was higher in July (1755 µg g−1) compared to October (1540 µg g−1). For the ∑RE concentrations, the difference was almost insignificant, with 276 µg g−1 in July versus 262 µg g−1 in October. While variation of the emission sources between the sampling sessions could be responsible for TE abundances, more stable sources could explain the RE behavior. For the ∑REEs, an opposite trend could be observed with a higher average concentration in October. However, more detailed aspects regarding the determined concentrations of REEs and potential sources will be discussed in Section 3.3.
Regarding the TE concentrations in both sampling sessions, the highest levels were found for Mn, followed by Ba, Zn, and Rb, an observation which may indicate a natural accumulation process [55]. Zinc concentration average values (132 µg g−1 in July and 124 µg g−1 in October) were consistent with the values reported by Sandu et al. [41] for stream sediments (Bistricioara, Romania). The mean concentration of Cr was 103 µg g−1 in July and 94.6 µg g−1 in October (Table 2), which are higher than the values reported in lake sediments from Bafa Lake, Turkey (75.90 ± 13.22 µg g−1 in summer season and 80.97 ± 29.74 µg g−1 in autumn season) [3].
The present study reports significantly higher concentrations for Cd (0.73–1.98 μg g−1 for July and 1.35–2.06 μg g−1 for October) compared with sediment from Lake Balaton, Hungary (0.1–0.7 μg g−1) reported by Nguyen et al. [56]. In the Nguyen et al. [56] study, lower concentrations for Cr (5.7–66 μg g−1) and Mn (160–760 μg g−1), similar concentrations for Co (1.7–17 μg g−1) and Zn (13–150 μg g−1), and higher concentrations for Pb (2.4–160 μg g−1) were presented. By comparison with a remote location, Bory Tucholskie National Park, Poland [49], used for the geochemical background values estimation, the REE concentrations in Podu Iloaiei Dam Lake sediment were higher and may indicate anthropogenic contribution in our case.
The Kolmogorov-Smirnov test indicated different distributions for Bi, Ce, Hf, Ir, Ta, U, Zr, Ba, Be, and Mn between the two sampling sessions. The element concentration range was considerably wider in July compared to October, and this might suggest a mixed source contribution to the element content in July (Table 2).
The relative distributions of quantified elements in July and October are presented in Figure 2. TEs, REs, and REEs contributed with 81.0%, 12.7%, and 6.3% in July (Figure 2a) and with 78.7%, 13.4%, and 8.0% in October (Figure 2c). In the ΣREEs, the LREEs and HREEs contributed with 71.8% and 28.2% in July (Figure 2b), respectively, 75.2% and 24.8% in October (Figure 2d). A decrease in TE contribution was observed from July to October, together with a slight increase for REEs. This is related to the increase in Ce contribution to the REEs from 24.4% in July to 31.3% in October. The individual contribution of other elements to ΣREEs was relatively constant.
The relative contributions of TEs for both sampling sessions varied in the following order: Mn > Ba > Zn > Rb > Zr > Li > Ga > Pb > As > Co > U > Be > Cd > Tl > Bi. As presented in Figure 2a,c, for Mn contribution, a slight decrease was observed from July to October.
For the RE group, the relative contributions varied in the following order: V > Cr > Nb > Ir > Hf > Mo > W > Ta, with V presenting a relatively constant contribution. For REEs, particularly Ce, La and Nd brought high contributions to the LREEs, while Y represented about half of the HREEs in both sampling sessions.

3.2. Spatial Distribution of Elements

IDW interpolation was performed on the element concentrations as a mapping method to obtain information about their spatial distributions. The selected TEs (Mn, Co, Zn, As, Cd, Ba, Tl, and Pb), REs (V, Cr, and Mo), Ce, LREEs* (LREEs without Ce), and HREEs spatial distributions in the sediment of Podu Iloaiei Dam Lake for July and October are presented in Figure 3. The distributions of other TEs (Li, Be, Ga, Rb, Zr, Bi, and U) and REs (Nb, Hf, Ta, W, and Ir) are presented in Figure S1 from the SM.
As observed from Figure 3, for Co, Zn, Cd, Tl, Pb, V, Cr, and Mo, lower concentration levels were identified in the south of the lake (S-4, S-7, and S-13 in Figure 1), close to the meadow and forest. Higher concentrations of these elements are observed in the middle and the north lake areas (S-3, S-5, and S-6 in Figure 1), with more susceptibility to be influenced by contributions of the residential area, railway, and transport road.
The profile of the spatial distribution of Mn, Ba, V, Cr, Mo, HREEs, and LREEs may indicate a displacement of the sediment from the input water area (in the west side) especially in July.
Similar patterns were observed in July for the distributions of As, Cd, and Pb. In October, As seemed to have a more pronounced contribution from the river water source, while Cd and Pb seemed to be associated with contributions from the road and railway traffic. Cadmium, Pb, and other elements’ distribution in soil samples collected from nearby roads were also highlighted in recent studies [57,58]. High concentrations of these elements were observed in the feeding area and near the dam. This may indicate the existence of a punctual and persistent source of pollution related to food supplied for fish in both sampling sessions.
Wastewater discharge and traffic might be important contributors for Co and Zn in both sampling sessions. Additionally, the feeding area presented higher concentrations of these elements.
Moreover, it is worth mentioning that V, Cr, and Mo presented higher concentrations in the feeding zone in both sampling sessions, with the higher concentration area being broader in October compared to July.
A complementary spatial distribution of Mn and Ce was observed. This indicated that the oxidative dissolution of Mn from the sediment particles can generate Ce enrichment from Ce4+ [59,60]. The cerium concentration showed the highest variability along the lake and between the two sampling periods, indicating different transport and deposition processes of the sediments in addition to redox conditions of the lake.
Element concentration variability was higher in July compared to October, with the coefficients of variation from 10 to 33% in July (19% in average) and from 5 to 26% in October (10% in average). This variability indicated that the urbanization, accumulation, and distribution of the potential toxic elements in lake sediments may lead to an increased spatial heterogeneity. For almost all REEs, the coefficient of variation values were higher in July (10 to 18%) compared to October (5 to 8%). These aspects indicated a homogeneous spatial distribution for these elements, except for Ce (with 33% in July and 26% in October). Differences in the input water quality, the weathering and leaching from soil surface, the local biochemical processes in the area with low water levels, and direct anthropogenic pressure could be important factors controlling the variability in the chemical composition of sediment samples from July. In October, the increase of the minimum concentration led to a reduced variability.
For measured REEs, the variations observed between sampling sessions may be related to the input water source, which may change the pH and salinity of the lake, seasonal changes in water redox conditions, and also due to an enhancement of the adsorption–desorption reactions and REEs complexation processes [15].

3.3. REEs Normalization and Anomalies

The distribution of Chondrite- and UCC-normalized REE concentrations for all sediment samples is presented in Figure 4. Negative anomalies can be observed for Ce and Eu elements, while the positive anomaly is identified for Gd in the Chondrite-normalized distribution. The UCC-normalized REE concentrations show non-flat patters, indicating that some fractionation may be present. The similar profiles of the REE normalized concentrations over the investigated area indicated common sources in sediments [16]. However, the presence of the Eu negative anomaly is relatively common in the upper crust [49,50]. The numerical evaluation of the Ce anomaly was performed using, in Equation (1), the UCC-normalized concentrations of its neighbors, La and Pr. The approach could not be applied for Gd because its left neighbor, Eu, is itself susceptible to anomalous behavior. Therefore, for the Gd anomaly evaluation, the Sm normalized concentration was used in Equation (2).
The negative Ce anomaly (Ce/Ce* < 1) was determined for all investigated samples, with values ranging from 0.33 to 0.78 in July and from 0.37 to 0.92 in October. The individual values of the Ce anomaly are included in Table S2 from the SM. Data presented in Figure 4 show that few samples from the October sampling session (i.e., S-2, S-5, S-10, and S-12 in Figure 1) have an almost insignificant Ce anomaly. It is worthy to mention that S-10 and S-12 samples are associated with nutrient rich feeding areas with higher susceptibility to present more intense bacterial activity (that enhance Ce3+ oxidation to Ce4+) on the surface of the water in October than in July. Under these circumstances, the sedimentation of Ce4+ due to scavenging by Fe-Mn oxide particles [15] may induce a reduced anomaly in the sediment samples, as observed in the present study. No correlation between Ce/Ce* and Mn content was observed for the July samples (Pearson correlation coefficient of −0.0005), while in October the correlation coefficient was as high as −0.7577, reflecting an actually significant negative correlation of the Ce anomaly and the Mn content. The correlation between the Ce anomaly and Mn content is in agreement with other reports [9,12,20,61,62].
Considering that the Mn concentration in sediment samples was fairly similar for both sampling sessions, it is assumed that, in October, the variation of the Ce anomaly was mainly due to the more pronounced oxidation potential of the water body. Only then, the formation of Ce4+ is possible, followed by Mn-containing particle scavenging, which finally leads to a diminished Ce anomaly in the sediment due to sedimentation.
A positive Gd anomaly was obtained for both July and October sediment samples. The values obtained for the Gd anomaly varied from 1.15 to 1.24 in July and from 1.18 to 1.28 in October. The individual values are presented in Table S2 from the SM. Although a slight increase from July to October can be observed, the geochemistry and redox condition of the lake were not controlling factors of Gd anomaly. Gadolinium-chelating agents are known to have applications in MRI [11,22,23]. For the investigated dam lake water, micropollution with Gd was associated with the largest medical hub in north-eastern Romania (Iasi, 30 km away from the investigated area, in average with about 10 MRI/1000 inhabitants especially during the last years). Similar results for water resources situated even at about 50 km away from medical hubs or highly populated areas are presented also in other studies [22,48]. Since the use of MRI in medical imaging is continuously growing (see Figure S2 from Supplementary Materials), the Gd micropollution in water resources situated nearby is expected to increase. Additionally, the discharge of treated wastewater into the environment can play a significant role for the positive Gd anomaly, as Gd-chelating agents are not successfully being removed in water treatment plants [22]. Moreover, Gd can have potential toxic effects due to its similar ionic radius, similar electronegativity, and bigger binding affinity to that of Ca2+ [15]. It is already suggested that Gd can block the Ca ion channels in cells [48,63] or may alter the catalytic effect of Ca2+ in enzymes [63,64].
As presented in Figure 4, for the UCC-normalized REE concentration distribution, medium mass REE concentrations higher than the values of the UCC may indicate an enrichment process. The Bell-Shaped index (BSI) was used to numerically evaluate the element enrichment in the investigated dam lake sediment. The individual values of the BSI are presented in Table S2, from the Supplementary Materials, together with the values of the Ce and Gd anomalies. The BSI values ranging from 1.21 to 1.28 in July (average of 1.25) and from 1.19 to 1.28 in October (average of 1.25) suggested stability in terms of the sampling points and sampling sessions. The enrichment of mREEs (Sm, Eu, Gd, Tb, and Dy) associated with phosphate minerals [65] was most probably related to geological processes [16,66,67] but may also indicate anthropogenic sources, such as phosphogypsum pollution [23] or the use of P-fertilizers. Although the anthropogenic contribution cannot be accurately estimated, since about 71% of the surrounding area of the Podu Iloaiei Dam Lake is used for agriculture [44], P-fertilizers can be suggested as a possible important contamination source.

3.4. Correlation Analysis and Principal Component Analysis

The correlations between the investigated elements (Pearson correlation coefficients) are presented in Figure 5, with major differences observed between the sampling sessions. Significant correlations (p < 0.05) were observed for As with Mn (0.60) and Pb (0.79) in July. In October, the correlations extended to Co (0.66), Zn (0.59), and some REEs (La:−0.58, Sm:−0.56, Gd:−0.54, Dy:−0.60, Ho:−0.61, Yb:−0.56, and Lu:−0.57). Cadmium was significantly correlated (p < 0.05) with Co (r = 0.86, 0.89), Zn (r = 0.86, 0.87), Bi (r = 0.79, 0.73), and Pb (r = 0.61, 0.71) in both sampling sessions (Figure 5). Significant correlations at p < 0.05 were observed between Mn and Ga (r = 0.69, 0.86), V (r = 0.63, 0.72), Cr (r = 0.76, 0.59), and Mo (r = 0.87, 0.83). Negative significant correlations were observed between Mn and Zr (r = −0.68), Hf (r = −0.68), and Ce (r = −0.69) in October. These observations suggested that the reductive dissolution of Mn from the sediment can generate Ce, Hf, and Zr enrichment because these elements are redox condition sensitive [59,60].
Significant strong correlations at p < 0.01 were observed between Zr and Hf, with a correlation coefficient of 0.99 in both sampling sessions. Such behavior is expected since Zr and Hf elements are very similar in their size and chemical properties due to the lanthanide contraction [68].
High values of Pearson coefficients were obtained between REEs (r > 0.81, p < 0.05), except for Ce in July. For October, a different behavior was observed, with a clear separation between LREEs and HREEs. Figure 5 shows significant correlations between LREEs (r > 0.77, p < 0.05), except for Ce, and between HREEs (r > 0.81, p < 0.05), except for Gd and Tb (Figure 5). These last two elements are more likely correlated with LREEs (r > 0.68, p < 0.05). In both sampling sessions, Nb and Ta elements were highly correlated with most of the REEs.
Furthermore, the PCA was performed on the database in order to identify the most important sources for TEs, REs, and REEs in the sediment of Podu Iloaiei Dam Lake. The results from the PCA applied to the database are presented in Figure 6 and Tables S3 and S4 from the Supplementary Materials. For optimal identification of the statistical solutions, the PCA was performed on the data set obtained by merging the results from both sampling sessions. The elements were divided into two groups, the first one including TEs, REs, Ce, LREEs*, and HREEs, while the second group was represented by all REEs. The measure of sampling adequacy by KMO statistics provided the values of 0.62 for the first group and 0.85 for the last one, while the significance of Bartlett’s test of sphericity was < 0.001 in both cases. Both tests led to the conclusion that the database size was suitable to be evaluated by the PCA. The number of factors with eigenvalues greater than 1.0 changed from three for the first group of elements to two for the second one, with 85.6% and, respectively, 93.5% of the variance explained. As presented in Table 3, the factor loadings were classified into strong and moderate, with the loading values of >0.75 and 0.75–0.50, respectively.
Figure 6 presents the loading of PC1 versus PC2 factors from the PCA applied to the investigated lake. Factor 1 in the first group of elements (TEs, REs, Ce, LREEs*, HREEs) had strong loadings of Rb, Cr, Li, Ba, Ga, Be, V, and W and moderate loadings of HREEs, Co, Mn, Mo, Nb, Tl, and Zn. This factor can be attributed to a mixed source of rock weathering and leaching from the soil surface. Work performed in the research group allowed for suggesting a potential additional contribution to this factor from atmospheric deposition (dry or wet). Stormwater residuals from the high traffic road located in the proximity of the lake can be another source [69,70]. In Factor 1, strong loading was observed for W (0.81) with moderate loadings for Mo (0.64) and Zn (0.59). This agrees with other studies that indicated W, Mo, and Zn to be tracers of anthropogenic pollution from the urban road dust [69,70].
Factor 2 was comprised of Ta, LREEs*, Ce, Hf, Zr, Bi, U, Nb, and Tl, corresponding to pedological characteristics of the investigated area. This may also highlight the Ce, Hf, and Zr elements sensitivity to the redox status of the sediment [71]. Factor 3 had strong loadings of Pb, As, and Cd with moderate loadings of Zn, Co, and Mo. This factor can be attributed to mixed sources of wastewater discharge, leaching from agricultural lands, or waste accumulation in unauthorized areas surrounding the Bahluet River.
The PCA analysis applied to the REE database discriminated two factors. Factor 1 had strong loadings of HREEs (Er, Tm, Y, Yb, Lu, Ho, and Dy), except for Gd and Tb, which presented strong loadings in Factor 2 among the LREEs (Ce, Pr, Nd, Sm, and La). A slight distancing of Ce from the Factor 2 elements can be observed in Figure 6b, indicating additional specific influences.

3.5. Contamination Assessment and Potentially Harmful Elements

Figure 7a presents the contamination degree assessment of the individual contaminant, by CF, at the interest sampling points. For a high number of elements (Li, Be, Co, Rb, Zr, Ba, Tl, U, and Mo) the contamination factor values were smaller than 1, indicating a low degree of contamination in the sediment samples. Vanadium, Cr, Nb, and Hf showed moderate contamination, with CF values close to 1. Low or moderate degrees of contamination (CF < 3) were observed for all elements, except for Cd, in all sampling points. CFs for Cd showed considerable contamination (CF = 3 ÷ 6) at almost all the sampling points. The S-4, S-7, S-13, and S-14 sampling points were exceptions, with CF values between considerable and moderate contamination (CF = 1 ÷ 3). The variation in CFs of Cd, a very toxic element in the aquatic ecosystem, which tends to bioaccumulate in sediment [72], reflected different source contributions between the sampling points. Cadmium contamination in Podu Iloaiei Dam Lake may be caused by fuel used for fishing boats or phosphate fertilizers and pesticides used in agricultural areas [3]. The railway line (606 magistral) and high road traffic (E58 road), located near the lake, might be important sources too.
The CFs of Zn, Ga, As, and Pb showed moderate contamination for all sampling points, with the highest values obtained for S-8, S-9, S-10, S-11 and S-12, corresponding to the feeding area. Natural weathering of minerals is also a known source for Zn [73]. The use of Pb-As insecticides and arsenic-based herbicides [74] on agricultural land in neighboring areas could contribute to the Pb and As concentrations in the sediment from Podu Iloaiei Dam Lake.
Figure 7b presents the results from the PLI analysis, with values ranging from 0.72 to 1.04 for TEs, and from 0.82 to 1.20 for REs. Several sampling points (S-1, S-9, and S-11 in Figure 1) from the Podu Iloaiei Dam Lake had PLI values for TEs very close to 1, indicating a slightly progressive deterioration of the environment in these locations. With the exception of the S-13 and S-14 samples, the PLI values for REs were above 1 in all locations, mainly due to sediment contamination with V, Cr (from fuel combustion, [75,76]), Nb, and Ta (from e-wastes, [77]). The sampling points associated with the feeding area (S-09, S-10, S-11, and S-12) seemed to have some of the highest PLI values for REs, highlighting a supplementary pollution source. This may be related to boat fuel, which is used more frequently in the feeding area. Domestic wastewater from settlements with no sewage infrastructure near S-1, S-3, and S-5 could be the main source of pollutants in this part of the lake.
The PLI values obtained for the S-1 and S-2 sampling points (Figure 1) located near a large area covered by plants highlighted the potential role of aquatic plants in absorbing and accumulating heavy metals [3].
The results obtained with the geoaccumulation index (Igeo) are shown in Figure 7c. The Igeo values indicate that the sediment was uncontaminated (Igeo ≤ 0) for most of the investigated TEs and Res. None to moderate contamination (Igeo = 0 ÷ 1) with Zn, Ga, As, and Ta was observed. For Cd, element moderate contamination (Igeo = 1 ÷ 2) was found at 10 sampling locations, confirming the contributions from anthropogenic pollution sources. Only at the S-4, S-7, S-13, and S-14 sampling locations, the Igeo values for Cd were below 1. The highest Igeo values for Zn and Ga were found at S-3. Sediments at the S-9, S-10, and S-11 locations were hinted as uncontaminated to moderately contaminated with Ga, As, and Ta. In specific agricultural activities related to crops of sunflower, wheat, corn, and rapeseed, fruit-bearing orchards in the area, the use of agrochemicals (e.g., herbicides, insecticides) generally used at the European level can be responsible for As pollution, as suggested by Defarge et al. [74]. Another possible source for As in north-eastern Romania may be represented by atmospheric depositions. Galon-Negru et al. [78] reported the largest concentration of As in summer and autumn, associated with a bursting increase in the coal combustion process. Moreover, the study showed that the water-soluble As from the fine particulate matter in the Iasi area had the largest contribution to inhalation and carcinogenicity risks. Possible contamination with Ga and Ta may be attributed to improper management of electronics in the residential area, with these elements being assigned as contaminants from e-wastes [79].
The influence of anthropogenic sources on the elements in the sediment samples was also assessed by determining the EF values. The calculated element EFs are shown in Figure 7d. For all sampling points, EF values less than 1.5 were obtained for the investigated elements, indicating their most probable natural origin, with few exceptions being observed for Zn, Ga, As, and Ta. Cadmium showed the highest EF values, with moderate anthropogenic modification (3 < EF < 5) for the S-3, S-8, S-9, and S-11 sampling points and minor anthropogenic modification (1.5 < EF < 3) for the other ten sampling points from Podu Iloaiei Dam Lake. The obtained results confirmed the anthropogenic source of Cd, with emissions associated mainly with fuel used for boats, railways, and road traffic.
The potentially harmful element (Cr, Zn, As, Cd, and Pb) concentration in the sediment samples of Podu Iloaiei Dam Lake, along with threshold effect concentration (TEC) and probable effect concentration (PEC) for each element, are shown in Figure 8.
Among the investigated elements, Cd and Pb showed concentrations that were mainly below TEC in all samples. The obtained results indicated that these elements are very unlikely to pose a threat to the organisms living in the bottom of the aquatic system [3] in the investigated location. For Cr and As, the concentrations were higher than the TECs for all sampling points, both in the July and October sampling sessions. The occurrence of different anthropogenic sources or processes during the two sampling sessions could be responsible for the variability in Zn concentrations, which were higher than the TEC.
The concentrations of Zn, As, Cd, and Pb did not exceed the PEC values for any of the investigated samples. Values higher than the PEC thresholds were identified for Cr in several samples (five sampling points in July and one in October), and such behavior may hint that Cr is one of the most important elements that could lead to toxicity in the investigated area. The sources of Cr pollution in Podu Iloaiei Dam Lake may be related to domestic wastewater from villages without centralized sewerage and fuel combustion associated with road traffic, both sources being located in the northern part of the lake.

4. Conclusions

High variabilities of TE and RE concentrations in the Podu Iloaiei Dam Lake sediments were identified. For both TEs and REs, indication of their main geogenic origin was obtained. The concentrations of Pb, As, Cd, and Zn were mainly related to anthropogenic sources.
The REEs had low variability, except for Ce. Normalized distributions of REEs and the observation of an enrichment of mREEs (Sm, Eu, Gd, Tb, and Dy), evaluated by BSI, indicated shared sources such as phosphates that may have both natural and anthropogenic origins. The statistical evaluation of redox-sensitive elements, such as Ce and Mn, and the Ce anomaly highlighted a higher oxidizing potential of the aqueous environment in October compared to July. A positive Gd anomaly was observed for the Chondrite-normalized concentration in both the July and October sampling sessions, which indicated micropollution and may be related to the use of Gd-chelating agents in MRI, a medical investigation technique with a growing use in the Iasi medical hub.
The assessment of the pollution indices emphasized that the sediment samples were highly enriched in Cd and contaminated with elements such as Zn, Ga, As, Pb, V, Cr, Nb, and Ta, related to anthropogenic sources. Sediment quality guidelines indicated that Cr posed a considerable threat to the aquatic biota.
The present study brought high contribution to understanding the specific sources of REEs, TEs, and REs in the aquatic system used for intensive fish farming, from north-eastern Romania. The use of agrochemicals (fertilizers, herbicides, and insecticides), road and railway traffic, and domestic wastewater discharge were found to be the major contributors to the potential harmful element contents in the Podu Iloaiei Dam Lake sediment.
A close monitorization of the Bahluet River, which is the main water input on the investigated Podu Iloaiei Lake, and even extending the investigation on the Bahlui hydrographic basin may be necessary in order to fully understand the evolution of metal pollutants and micropollutants.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/toxics10050242/s1, Figure S1. Spatial distribution of Be, Bi, Ga, Hf, Ir, Li, Nb, Rb, Ta, W, U, and Zr in sediment of Podu Iloaiei Dam Lake for July and October sampling sessions; Figure S2. Number of MRI exams performed per 1000 inhabitants in some European countries. Data compiled from OECD (2022), “Magnetic resonance imaging (MRI) exams” (indicator), https://0-doi-org.brum.beds.ac.uk/10.1787/1d89353f-en (accessed on 29 March 2022) and EUROSTAT, Medical technologies—examinations by medical imaging techniques (CT, MRI and PET), http://appsso.eurostat.ec.europa.eu/nui/submitViewTableAction.do (accessed on 02 May 2022); Table S1. Limit of quantification (LoQ) for inductively coupled plasma mass spectrometry determined REEs in the water and sediment samples; Table S2. The numerical evaluation of Ce and Gd anomalies and the BSI for the investigated sediment samples; Table S3. Principal component analysis data for TE, RE, Ce, ΣLREE*, and ΣHREE quantified in sediment samples of Podu Iloaiei Dam Lake. Bold values present loadings > 0.75; Table S4. Principal component analysis data for REE quantified in sediment samples of Podu Iloaiei Dam Lake. Bold values present loadings > 0.75.

Author Contributions

Conceptualization, C.A. (Cecilia Arsene) and R.I.O.; methodology, C.A. (Cornelia Amarandei), L.V.S., A.G.N., C.A. (Cecilia Arsene) and R.I.O.; validation, C.A. (Cornelia Amarandei), L.V.S., A.G.N., C.A. (Cecilia Arsene) and R.I.O.; formal analysis, C.A. (Cornelia Amarandei), L.V.S. and A.G.N.; investigation, C.A. (Cornelia Amarandei) and L.V.S.; resources, C.A. (Cecilia Arsene) and R.I.O.; data curation, C.A. (Cornelia Amarandei), L.V.S., A.G.N., C.A. (Cecilia Arsene) and R.I.O.; writing—original draft preparation, C.A. (Cornelia Amarandei), L.V.S. and A.G.N.; writing—review and editing, C.A. (Cornelia Amarandei), L.V.S., A.G.N., C.A. (Cecilia Arsene) and R.I.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors express their gratitude for the Administration of the Podu Iloaiei Dam Lake facilitating the access in the area for sampling purposes. Acknowledgment is given by R.I.O. and C.A. (Cecilia Arsene) to infrastructure support from the Operational Program Competitiveness 2014–2020, Axis 1, under POC/448/1/1 Research infrastructure projects for public R&D institutions/Sections F 2018, through the Research Center with Integrated Techniques for Atmospheric Aerosol Investigation in Romania (RECENT AIR) project, under grant agreement MySMIS no. 127324. Authors are thankful to Romanian Ministry of Research, Innovation and Digitization, within Program 1—Development of the national RD system, Subprogram 1.2—Institutional Performance—RDI excellence funding projects, Contract no.11PFE/30.12.2021, for financial support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Paul, V.; Sankar, M.S.; Vattikuti, S.; Dash, P.; Arslan, Z. Pollution assessment and land use land cover influence on trace metal distribution in sediments from five aquatic systems in southern USA. Chemosphere 2021, 263, 128243. [Google Scholar] [CrossRef] [PubMed]
  2. Birch, G.F.; Olmos, M.A. Sediment-bound heavy metals as indicators of human influence and biological risk in coastal water bodies. ICES J. Mar. Sci. 2008, 65, 1407–1413. [Google Scholar] [CrossRef]
  3. Algul, F.; Beyhan, M. Concentrations and sources of heavy metals in shallow sediments in Lake Bafa, Turkey. Sci. Rep. 2020, 10, 11782. [Google Scholar] [CrossRef] [PubMed]
  4. Gao, L.; Wang, Z.; Shan, J.; Chen, J.; Tang, C.; Yi, M.; Zhao, X. Distribution characteristics and sources of trace metals in sediment cores from a trans-boundary watercourse: An example from the Shima River, Pearl River Delta. Ecotoxicol. Environ. Saf. 2016, 134, 186–195. [Google Scholar] [CrossRef]
  5. Wu, S.; Peng, S.; Zhang, X.; Wu, D.; Luo, W.; Zhang, T.; Zhou, S.; Yang, G.; Wan, H.; Wu, L. Levels and health risk assessments of heavy metals in urban soils in Dongguan, China. J. Geochem. Explor. 2015, 148, 71–78. [Google Scholar] [CrossRef]
  6. Strungaru, S.A.; Nicoara, M.; Teodosiu, C.; Baltag, E.; Ciobanu, C.; Plavan, G. Patterns of toxic metals bioaccumulation in a cross-border freshwater reservoir. Chemosphere 2018, 207, 192–202. [Google Scholar] [CrossRef]
  7. Yi, Y.; Yang, Z.; Zhang, S. Ecological risk assessment of heavy metals in sediment and human health risk assessment of heavy metals in fishes in the middle and lower reaches of the Yangtze River Basin. Environ. Pollut. 2011, 159, 2575–2585. [Google Scholar] [CrossRef]
  8. Sojka, M.; Choinski, A.; Ptak, M.; Siepak, M. The variability of lake water chemistry in the Bory Tucholskie National Park (Northern Poland). Water 2020, 12, 394. [Google Scholar] [CrossRef] [Green Version]
  9. Liu, H.; Guo, H.; Pourret, O.; Wang, Z.; Sun, Z.; Zhang, W.; Liu, M. Distribution of rare earth elements in sediments of the North China Plain: A probe of sedimentation process. Appl. Geochem. 2021, 134, 105089. [Google Scholar] [CrossRef]
  10. Hu, Z.; Richter, H.; Sparovek, G.; Schnug, E. Physiological and biochemical effects of rare earth elements on plants and their agricultural significance: A Review. J. Plant Nutr. 2004, 27, 183–220. [Google Scholar] [CrossRef]
  11. Hatje, V.; Bruland, K.W.; Flegal, A.R. Increases in anthropogenic Gadolinium anomalies and rare earth element concentrations in San Francisco Bay over a 20 Year record. Environ. Sci. Technol. 2016, 50, 4159–4168. [Google Scholar] [CrossRef] [PubMed]
  12. Bodin, S.; Godet, A.; Westermann, S.; Föllmi, K.B. Secular change in northwestern Tethyan water-mass oxygenation during the late Hauterivian-early Aptian. Earth Planet. Sci. Lett. 2013, 374, 121–131. [Google Scholar] [CrossRef]
  13. Nozaki, Y. Rare earth elements and their isotopes in the Ocean. In Encyclopedia of Ocean Sciences, 1st ed.; Steele, J.H., Ed.; Academic Press: Cambridge, MA, USA, 2001; pp. 2354–2366. [Google Scholar] [CrossRef]
  14. Elderfield, H.; Greaves, M.J. The rare earth elements in seawater. Nature 1982, 296, 214–219. [Google Scholar] [CrossRef]
  15. Laveuf, C.; Cornu, S. A Review on the potentiality of rare earth elements to trace pedogenetic processes. Geoderma 2009, 154, 1–12. [Google Scholar] [CrossRef]
  16. Zhang, X.; Zhang, F.; Chen, X.; Zhang, W.; Deng, H. REEs fractionation and sedimentary implication in surface sediments from eastern South China Sea. J. Rare Earths 2012, 30, 614–620. [Google Scholar] [CrossRef]
  17. Mclennan, S.M. Rare earth elements in sedimentary rocks: Influence of provenance and sedimentary processes. In Geochemistry and Mineralogy of Rare Earth Elements; Lipin, B.R., McKay, G.A., Eds.; De Gruyter: Berlin, Germany, 1989; pp. 169–200. [Google Scholar] [CrossRef]
  18. Greenwood, N.N.; Earnshaw, A. Chemistry of the Elements, 2nd ed.; Butterworth-Heinemann: Oxford, UK, 1997; pp. 1227–1249. [Google Scholar] [CrossRef]
  19. Elderfield, H.; Whitfield, M.; Burton, J.D.; Bacon, M.P.; Liss, P.S. The oceanic chemistry of the rare-earth elements. Phil. Trans. R. Soc. Lond. A 1988, 325, 105–126. [Google Scholar] [CrossRef]
  20. German, C.R.; Holliday, B.P.; Elderfield, H. Redox cycling of rare earth elements in the suboxic zone of the Black Sea. Geochim. Cosmochim. Acta 1991, 55, 3553–3558. [Google Scholar] [CrossRef]
  21. Tostevin, R.; Shields, G.A.; Tarbuck, G.M.; He, T.; Clarkson, M.O.; Wood, R.A. Effective use of Cerium anomalies as a redox proxy in carbonate-dominated marine settings. Chem. Geol. 2016, 438, 146–162. [Google Scholar] [CrossRef] [Green Version]
  22. Ebrahimi, P.; Barbieri, M. Gadolinium as an emerging microcontaminant in water resources: Threats and opportunities. Geosciences 2019, 9, 93. [Google Scholar] [CrossRef] [Green Version]
  23. Wang, L.; Han, X.; Ding, S.; Liang, T.; Zhang, Y.; Xiao, J.; Dong, L.; Zhang, H. Combining multiple methods for provenance discrimination based on rare earth element geochemistry in lake sediment. Sci. Total Environ. 2019, 672, 264–274. [Google Scholar] [CrossRef]
  24. Gu, X.; Lin, C.; Wang, B.; Wang, J.; Ouyang, W. A comprehensive assessment of anthropogenic impacts, contamination, and ecological risks of toxic elements in sediments of urban rivers: A case study in Qingdao, East China. Environ. Adv. 2022, 7, 100143. [Google Scholar] [CrossRef]
  25. Xia, F.; Hu, B.; Zhu, Y.; Ji, W.; Chen, S.; Xu, D.; Shi, Z. Improved mapping of potentially toxic elements in soil via integration of multiple data sources and various geostatistical methods. Remote Sens. 2020, 12, 3775. [Google Scholar] [CrossRef]
  26. Agyeman, P.C.; Ahado, S.K.; Boruvka, L.; Biney, J.K.M.; Sarkodie, V.Y.O.; Kebonye, N.M.; Kingsley, J. Trend analysis of global usage of digital soil mapping models in the prediction of potentially toxic elements in soil/sediments: A bibliometric review. Environ. Geochem. Health 2021, 43, 1715–1739. [Google Scholar] [CrossRef] [PubMed]
  27. Alexakis, D.; Kokmotos, I.; Gamvroula, D.; Varelidis, G. Wildfire effects on soil quality: Application on a suburban area of West Attica (Greece). Geosci. J. 2021, 25, 243–253. [Google Scholar] [CrossRef]
  28. Doufexi, M.; Gamvroula, D.E.; Alexakis, D.E. Elements’ content in stream sediment and wildfire ash of suburban areas in West Attica (Greece). Water 2022, 14, 310. [Google Scholar] [CrossRef]
  29. Marzola, E.L.; Rosa, L.M.T.; Toppa, R.H.; Martines, M.R.; Pitombo, L.M.; Cavagis, A.D.M.; do Carmo, J.B.; Botero, W.G.; de Oliveira, L.C. Spatial modeling applied to environmental monitoring: Identifying sources of potentially toxic metals in aquatic system. Environ. Geochem. Health 2019, 41, 1339–1350. [Google Scholar] [CrossRef]
  30. Oleiwi, A.S.; Al-Dabbas, M. Assessment of contamination along the Tigris River from Tharthar-Tigris canal to Azizziyah, middle of Iraq. Water 2022, 14, 1194. [Google Scholar] [CrossRef]
  31. 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]
  32. 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. Meeresun. 1980, 33, 566–575. [Google Scholar] [CrossRef] [Green Version]
  33. Angulo, E. The Tomlinson Pollution Load Index applied to heavy metal, “Mussel-Watch” data: A useful index to assess coastal pollution. Sci. Total Environ. 1996, 187, 19–56. [Google Scholar] [CrossRef]
  34. Loska, K.; Cebula, J.; Pelczar, J.; Wiechula, D.; Kwapulinski, J. Use of enrichment, and contamination factors together with geoaccumulation indexes to evaluate the content of Cd, Cu, and Ni in the Rybnik water reservoir in Poland. Water. Air. Soil Pollut. 1997, 93, 347–365. [Google Scholar] [CrossRef]
  35. Kowalska, J.B.; Mazurek, R.; Gasiorek, M.; Zaleski, T. Pollution indices as useful tools for the comprehensive evaluation of the degree of soil contamination—A review. Environ. Geochem. Health 2018, 40, 2395–2420. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Buat-Menard, P.; Chesselet, R. variable influence of the atmospheric flux on the trace metal chemistry of oceanic suspended matter. Earth Planet. Sci. Lett. 1979, 42, 399–411. [Google Scholar] [CrossRef]
  37. Barbieri, M. The importance of enrichment factor (EF) and geoaccumulation index (Igeo) to evaluate the soil contamination. J. Geol. Geophys. 2016, 5, 1000237. [Google Scholar] [CrossRef]
  38. Long, E.R.; MacDonald, D.D. Recommended uses of empirically derived, sediment quality guidelines for marine and estuarine ecosystems. Hum. Ecol. Risk Assess. 1998, 4, 1019–1039. [Google Scholar] [CrossRef]
  39. MacDonald, D.D.; Ingersoll, C.G.; Berger, T.A. Development and evaluation of consensus-based sediment quality guidelines for freshwater ecosystems. Arch. Environ. Contam. Toxicol. 2000, 39, 20–31. [Google Scholar] [CrossRef]
  40. 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]
  41. Sandu, M.C.; Soroaga, L.V.; Balaban, S.I.; Chelariu, C.; Chiscan, O.; Iancu, G.O.; Arsene, C.; Olariu, R.I. Trace elements distribution in stream sediments of an abandoned U mining site in the Eastern Carpathians, Romania, with particular focus on REEs. Geochemistry 2021, 81, 125761. [Google Scholar] [CrossRef]
  42. Amarandei, C.; Negru, A.G.; Soroaga, L.V.; Cucu-Man, S.M.; Olariu, R.I.; Arsene, C. Assessment of surface water quality in the Podu Iloaiei Dam Lake (North-Eastern Romania): Potential implications for aquaculture activities in the area. Water 2021, 13, 2395. [Google Scholar] [CrossRef]
  43. WATER PIN Limited Liability Company. Securing the Podu Iloaiei Reservoir, on the Bahluet River, Iasi County, Synthesis Memoir. 2011. Available online: https://www.scribd.com/document/375298910/Memoriu-de-Sinteza-Ac-Podu-Iloaiei (accessed on 9 March 2022).
  44. Climate Change—Department for Waters Forests and Fisheries National Administration “Romanian Waters” Management-Bucharest, National Institute of Hydrology and Water, Prut-Barlad River Basin Management Plan, 2016–2021, in Accordance with Article No. 2015. Available online: https://vdocuments.mx/plan-management-prut-barlad-31032015.html (accessed on 14 March 2022).
  45. Lu, Y.; Li, G.; Liu, W.; Yuan, H.; Xiao, D. The application of microwave digestion in decomposing some refractory ore samples with solid fusion agent. Talanta 2018, 186, 538–544. [Google Scholar] [CrossRef]
  46. Potts, P.J. Inductively coupled plasma—mass spectrometry. In A Handbook of Silicate Rock Analysis; Potts, P.J., Ed.; Springer: Boston, MA, USA, 1992; pp. 575–586. [Google Scholar] [CrossRef]
  47. Taylor, S.R.; McLennan, S.M. The geochemical the continental evolution crust. Rev. Geophys. 1995, 33, 241–265. [Google Scholar] [CrossRef]
  48. Pedreira, R.M.A.; Pahnke, K.; Böning, P.; Hatje, V. tracking hospital effluent-derived gadolinium in Atlantic coastal waters off Brazil. Water Res. 2018, 145, 62–72. [Google Scholar] [CrossRef] [PubMed]
  49. Sojka, M.; Choinski, A.; Ptak, M.; Siepak, M. Causes of variations of trace and rare earth elements concentration in lakes bottom sediments in the Bory Tucholskie National Park, Poland. Sci. Rep. 2021, 11, 244. [Google Scholar] [CrossRef]
  50. Song, H.; Shin, W.-J.; Ryu, J.-S.; Shin, H.S.; Chung, H.; Lee, K.-S. Anthropogenic rare earth elements and their spatial distributions in the Han River, South Korea. Chemosphere 2017, 172, 155–165. [Google Scholar] [CrossRef] [PubMed]
  51. Turekian, K.K.; Wedepohl, H. Distribution of the elements in some major units of the Earth’s Crust. Geol. Soc. Am. Bull. 1961, 72, 175–192. [Google Scholar] [CrossRef]
  52. Hakanson, L. An ecological risk index for aquatic pollution control. A sedimentological approach. Water Res. 1980, 14, 975–1001. [Google Scholar] [CrossRef]
  53. Harikrishnan, N.; Ravisankar, R.; Gandhi, M.S.; Kanagasabapathy, K.; Prasad, M.V.R.; Satapathy, K. Heavy metal assessment in sediments of east coast of Tamil Nadu using energy dispersive X-ray fluorescence spectroscopy. Radiat. Prot. Environ. 2017, 40, 21–26. [Google Scholar] [CrossRef]
  54. Abata, E.O.; Aiyesanmi, F.; Adebayo, A.O.; Ajayi, O.O.; Kazuhiko, T.; Hiroshi, S. Contamination and ecological risk assessment of heavy metal in the sediment of Ala River, South—West Nigeria: An index—Analysis approach. Int. J. Sci. Eng. Res. 2016, 7, 1392–1408. [Google Scholar]
  55. Saygı, Y.; Yigit, S.A. Heavy metals in Yeniçaǧa Lake and its potential sources: Soil, water, sediment, and plankton. Environ. Monit. Assess. 2012, 184, 1379–1389. [Google Scholar] [CrossRef]
  56. Nguyen, H.L.; Leermakers, M.; Osan, J.; Torok, S.; Baeyens, W. Heavy metals in Lake Balaton: Water column, suspended matter, sediment and biota. Sci. Total Environ. 2005, 340, 213–230. [Google Scholar] [CrossRef]
  57. Alexakis, D.E. Multielement contamination of land in the margin of highways. Land 2021, 10, 230. [Google Scholar] [CrossRef]
  58. Hong, N.; Guan, Y.; Yang, B.; Zhong, J.; Zhu, P.; Ok, Y.S.; Hou, D.; Tsang, D.C.W.; Guan, Y.; Liu, A. Quantitative source tracking of heavy metals contained in urban road deposited sediments. J. Hazard. Mater. 2020, 393, 122362. [Google Scholar] [CrossRef] [PubMed]
  59. Cabral, A.R.; Zeh, A.; Vianna, N.C.; Ackerman, L.; Pašava, J.; Lehmann, B.; Chrastný, V. Molybdenum-isotope signals and cerium anomalies in Palaeoproterozoic manganese ore survive high-grade metamorphism. Sci. Rep. 2019, 9, 4570. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Kechiched, R.; Laouar, R.; Bruguier, O.; Kocsis, L.; Salmi-Laouar, S.; Bosch, D.; Ameur-Zaimeche, O.; Foufou, A.; Larit, H. Comprehensive REE + Y and sensitive redox trace elements of Algerian Phosphorites (Tébessa, eastern Algeria): A geochemical study and depositional environments tracking. J. Geochem. Explor. 2020, 208, 106396. [Google Scholar] [CrossRef]
  61. Goldstein, S.J.; Jacobsen, S.B. Rare Earth Elements in River Waters. Earth Planet. Sci. Lett. 1988, 89, 35–47. [Google Scholar] [CrossRef]
  62. Kakuwa, Y.; Matsumoto, R. Cerium negative anomaly just before the Permian and Triassic boundary event—The upward expansion of anoxia in the water column. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2006, 229, 335–344. [Google Scholar] [CrossRef]
  63. Sherry, A.D.; Caravan, P.; Lenkinski, R.E. Primer on gadolinium chemistry. J. Magn. Reson. Imaging 2009, 30, 1240–1248. [Google Scholar] [CrossRef] [Green Version]
  64. Darnall, D.W.; Birnbaum, E.R. Lanthanide ions activate α-amylase. Biochemistry 1973, 12, 3489–3491. [Google Scholar] [CrossRef]
  65. Hannigan, R.E.; Sholkovitz, E.R. The development of middle rare earth element enrichments in freshwaters: Weathering of phosphate minerals. Chem. Geol. 2001, 175, 495–508. [Google Scholar] [CrossRef]
  66. Lecuyer, C.; Reynard, B.; Grandjean, P. Rare earth element evolution of Phanerozoic seawater recorded in biogenic apatites. Chem. Geol. 2004, 204, 63–102. [Google Scholar] [CrossRef]
  67. Sholkovitz, E.R. The aquatic chemistry of rare earth elements in rivers and estuaries. Aquat. Geochem. 1995, 1, 1–34. [Google Scholar] [CrossRef]
  68. Orians, K.J.; Merrin, C.L. Refractory Metals. In Encyclopedia of Ocean Sciences, 1st ed.; Steele, H., Thorpe, S., Turekian, K.K., Eds.; Academic Press: Cambridge, MA, USA, 2001; Volume 4, pp. 2387–2399. [Google Scholar] [CrossRef]
  69. Nazzal, Y.; Rosen, M.A.; Al-Rawabdeh, A.M. Assessment of metal pollution in urban road dusts from selected highways of the Greater Toronto Area in Canada. Environ. Monit. Assess. 2013, 185, 1847–1858. [Google Scholar] [CrossRef] [PubMed]
  70. Kasimov, N.S.; Vlasov, D.V.; Kosheleva, N.E. Urban climate enrichment of road dust particles and adjacent environments with metals and metalloids in Eastern Moscow. Urban Clim. 2020, 32, 100638. [Google Scholar] [CrossRef]
  71. Censi, P.; Raso, M.; Yechieli, Y.; Ginat, H.; Saiano, F.; Zuddas, P.; Brusca, L.; D’Alessandro, W.; Inguaggiato, C. Geochemistry of Zr, Hf, and REE in a wide spectrum of Eh and water composition: The case of Dead Sea Fault system (Israel). Geochem. Geophys. Geosyst. 2017, 18, 844–857. [Google Scholar] [CrossRef]
  72. World Health Organization, Nordic Council of Ministers. Cadmium Review. Report No. 1, Issue No. 4, Geneva, 2003. Available online: https://cdn.who.int/media/docs/default-source/chemical-safety/cadmium/nmr-cadmium.pdf?sfvrsn=36032d8_2&download=true (accessed on 23 March 2022).
  73. Robson, T.C.; Braungardt, C.B.; Rieuwerts, J.; Worsfold, P. Cadmium contamination of agricultural soils and crops resulting from sphalerite weathering. Environ. Pollut. 2014, 184, 283–289. [Google Scholar] [CrossRef] [Green Version]
  74. Defarge, N.; Spiroux de Vendomois, J.; Seralini, G.E. Toxicity of formulants and heavy metals in glyphosate-based herbicides and other pesticides. Toxicol. Rep. 2018, 5, 156–163. [Google Scholar] [CrossRef]
  75. Hope, B.K. A dynamic model for the global cycling of anthropogenic vanadium. Glob. Biogeochem. Cycles 2008, 22, 1–16. [Google Scholar] [CrossRef] [Green Version]
  76. Fiedler, S.; Siebe, C.; Herre, A.; Roth, B.; Cram, S.; Stahr, K. Contribution of oil industry activities to environmental loads of heavy metals in the Tabasco Lowlands, Mexico. Water Air Soil Pollut. 2009, 197, 35–47. [Google Scholar] [CrossRef]
  77. Bu-Olayan, A.H.; Thomas, B.V. Bourgeoning impact of the technology critical elements in the marine environment. Environ. Pollut. 2020, 265, 115064. [Google Scholar] [CrossRef]
  78. Galon-Negru, A.G.; Olariu, R.I.; Arsene, C. Size-resolved measurements of PM 2.5 water-soluble elements in Iasi, North-Eastern Romania: Seasonality, source apportionment and potential implications for human health. Sci. Total Environ. 2019, 695, 133839. [Google Scholar] [CrossRef]
  79. Filella, M.; Rodushkin, I. A concise guide for the determination of less-studied technology-critical elements (Nb, Ta, Ga, In, Ge, Te) by inductively coupled plasma mass spectrometry in environmental samples. Spectrochim. Acta-Part B At. Spectrosc. 2018, 141, 80–84. [Google Scholar] [CrossRef]
Figure 1. Map of the Podu Iloaiei Dam Lake and the distribution of sediment sampling sites (the Image © 2022 CNES/Airbus from Google Earth Pro was incorporated in the present figure).
Figure 1. Map of the Podu Iloaiei Dam Lake and the distribution of sediment sampling sites (the Image © 2022 CNES/Airbus from Google Earth Pro was incorporated in the present figure).
Toxics 10 00242 g001
Figure 2. Relative contribution averages for July ((a)—total quantified elements, (b)—rare earth elements) and October ((c)—total quantified elements, (d)—rare earth elements) in sediment samples collected from Podu Iloaiei Dam Lake.
Figure 2. Relative contribution averages for July ((a)—total quantified elements, (b)—rare earth elements) and October ((c)—total quantified elements, (d)—rare earth elements) in sediment samples collected from Podu Iloaiei Dam Lake.
Toxics 10 00242 g002
Figure 3. Spatial distribution of selected TEs (As, Ba, Cd, Co, Mn, Pb, Tl, and Zn), REs (Cr, Mo, and V) and Ce, ΣLREE*, and ΣHREE in the sediment of Podu Iloaiei Dam Lake for the July and October sampling sessions.
Figure 3. Spatial distribution of selected TEs (As, Ba, Cd, Co, Mn, Pb, Tl, and Zn), REs (Cr, Mo, and V) and Ce, ΣLREE*, and ΣHREE in the sediment of Podu Iloaiei Dam Lake for the July and October sampling sessions.
Toxics 10 00242 g003aToxics 10 00242 g003b
Figure 4. The Chondrite-(up) and UCC-normalized (down) REE concentration distribution in the investigated samples for both sampling sessions.
Figure 4. The Chondrite-(up) and UCC-normalized (down) REE concentration distribution in the investigated samples for both sampling sessions.
Toxics 10 00242 g004
Figure 5. Pearson correlation coefficient matrix for the relationships between the element concentrations in the Podu Iloaiei Dam Lake sediment in July (left) and October (right). The color and the diameter of the circles are related to the correlation coefficient value (according to the color scale on the right) and the p-value (inversely proportional), respectively. The number of white stars is related to the p-value: * (p < 0.05), ** (p < 0.01), and *** (p < 0.001).
Figure 5. Pearson correlation coefficient matrix for the relationships between the element concentrations in the Podu Iloaiei Dam Lake sediment in July (left) and October (right). The color and the diameter of the circles are related to the correlation coefficient value (according to the color scale on the right) and the p-value (inversely proportional), respectively. The number of white stars is related to the p-value: * (p < 0.05), ** (p < 0.01), and *** (p < 0.001).
Toxics 10 00242 g005
Figure 6. Loading of PC1 versus PC2 factors from TE, RE, Ce, ΣLREE*, and ΣHREE (a) and REE (b) quantified in the sediment samples of Podu Iloaiei Dam Lake (loading values are presented in Tables S3 and S4 in the Supplementary Materials).
Figure 6. Loading of PC1 versus PC2 factors from TE, RE, Ce, ΣLREE*, and ΣHREE (a) and REE (b) quantified in the sediment samples of Podu Iloaiei Dam Lake (loading values are presented in Tables S3 and S4 in the Supplementary Materials).
Toxics 10 00242 g006
Figure 7. Contamination factors (CFs), (a) pollution load indices (PLIs), (b) geoaccumulation indices (Igeos) (c), and enrichment factors (EFs) (d) calculated for the sediment sampling points from Podu Iloaiei Lake.
Figure 7. Contamination factors (CFs), (a) pollution load indices (PLIs), (b) geoaccumulation indices (Igeos) (c), and enrichment factors (EFs) (d) calculated for the sediment sampling points from Podu Iloaiei Lake.
Toxics 10 00242 g007
Figure 8. Box plots of the potentially harmful element (Cr, Zn, As, Cd, and Pb) concentration in the Podu Iloaiei Lake sediment, along with the threshold effect concentration (TEC) and probable effect concentration (PEC)) (thin black line—median; white square—mean; box—25–75% percentiles; length of the whiskers plot—10% and 90% of observed concentrations).
Figure 8. Box plots of the potentially harmful element (Cr, Zn, As, Cd, and Pb) concentration in the Podu Iloaiei Lake sediment, along with the threshold effect concentration (TEC) and probable effect concentration (PEC)) (thin black line—median; white square—mean; box—25–75% percentiles; length of the whiskers plot—10% and 90% of observed concentrations).
Toxics 10 00242 g008
Table 1. Sediment pollution indices and evaluation criteria for potential contamination assessment.
Table 1. Sediment pollution indices and evaluation criteria for potential contamination assessment.
Pollution IndexEquationEvaluation Criteria
Contamination factor (CF) [34,52,53] C F = C n B n <1, low degree of contamination
1 ÷ 3, moderate degree of contamination
3 ÷ 6, considerable degree of contamination
>6, very high degree of contamination
Pollution load index (PLI) [32,33] PLI = ( C F 1 × C F 2 × × C F n ) 1 n 0, no pollution
0 ÷ 1, only baseline pollutants
>1, progressive deterioration
Geoaccumulation index (Igeo) [1,37,54] I g e o = l o g 2 C n 1.5   ×   B n ≤0, no contamination
0 ÷ 1, none to moderate contamination
1 ÷ 2, moderate contamination
2 ÷ 3, moderate to heavy contamination
3 ÷ 4, heavy contamination
4 ÷ 5, heavy to extreme contamination
>5, extreme contamination
Enrichmentfactor (EF) [1,2,31,34,36] E F = C n / C M n B n / B M n 0 ÷ 1.5, natural processes
1.5 ÷ 3, minor anthropogenic modification
3 ÷ 5, moderate anthropogenic modification
5 ÷ 10, severe anthropogenic modification
>10, very severe anthropogenic modification
Note: Cn = measured concentration of element n; Bn = background concentration of element n; CMn = measured concentration of Mn; BMn = background concentration of Mn.
Table 2. Descriptive statistics for the TE, RE, and REE concentrations in sediment (μg g−1, n = 14/session) samples collected in July and October 2017 from Podu Iloaiei Dam Lake, Romania. Data are presented as mean (minimum–maximum).
Table 2. Descriptive statistics for the TE, RE, and REE concentrations in sediment (μg g−1, n = 14/session) samples collected in July and October 2017 from Podu Iloaiei Dam Lake, Romania. Data are presented as mean (minimum–maximum).
ElementJulyOctober
As19.8 (15.0–25.5)20.7 (15.2–25.0)
Ba358 (299–413)319 (280–381)
Be2.25 (1.23–2.97)1.96 (1.60–2.25)
Bi0.27 (0.12–0.35)0.33 (0.25–0.43)
Cd1.62 (0.73–1.98)1.81 (1.35–2.06)
Co12.4 (6.11–15.6)11.1 (9.40–13.3)
Ga37.1 (19.8–47.3)32.0 (25.6–39.6)
Li45.6 (20.4–62.9)39.9 (31.2–47.8)
Mn922 (694–1196)757 (587–1012)
Pb22.1 (13.3–27.3)23.2 (19.9–28.5)
Rb125 (59.2–171)107 (94.4–130)
Tl0.65 (0.35–0.83)0.70 (0.61–0.75)
U2.87 (1.58–3.49)3.20 (2.75–3.68)
Zn132 (62.3–194)125 (100–174)
Zr74.0 (25.5–115)97.8 (65.3–127)
ΣTE1755 (1241–2205)1540 (1311–1830)
Cr103 (45.3–143)94.6 (80.2–111)
Hf2.27 (0.90–3.33)3.02 (2.05–3.84)
Ir3.30 (1.68–4.78)4.38(3.40–6.03)
Mo1.83 (1.26–2.41)1.62 (1.12–2.30)
Nb13.2 (9.15–15.1)13.5 (11.1–14.7)
Ta1.19 (0.81–1.38)1.38 (1.11–1.49)
V149 (57.3–213)142 (117–168)
W1.79 (1.12–2.35)1.72 (1.39–1.95)
ΣRE276 (118–383)262 (231–306)
Ce33.5 (15.5–52.0)48.8 (23.5–63.5)
Eu1.20 (0.80–1.38)1.23 (1.09–1.33)
La28.5 (16.2–32.9)29.6 (24.1–32.7)
Nd23.8 (16.0–26.7)25.1 (22.3–27.6)
Pr6.48 (4.10–7.33)7.00 (6.03–7.65)
Sm5.10 (3.34–5.76)5.41 (4.83–5.78)
ΣLREE98.5 (57.5–124)117 (85.3–138)
Dy4.49 (3.32–5.05)4.66 (4.08–4.99)
Er2.59 (1.99–2.90)2.66 (2.29–2.99)
Gd5.33 (3.64–5.90)5.75 (4.91–6.34)
Ho0.89 (0.66–0.99)0.93 (0.80–1.04)
Lu0.36 (0.27–0.41)0.38 (0.33–0.45)
Tb0.76 (0.54–0.83)0.81 (0.68–0.84)
Tm0.37 (0.29–0.43)0.39 (0.36–0.45)
Y21.4 (16.5–24.5)20.5 (17.8–22.5)
Yb2.36 (1.73–2.74)2.50 (2.18–2.90)
ΣHREE38.6 (28.9–43.0)38.6 (33.5–42.2)
ΣREE137 (88.9–163)155 (118–178)
Table 3. Distribution of the elements in the factors extracted according to PCA for loadings > 0.50, loadings > 0.75 in bold, and the total variance (%) explained by each factor.
Table 3. Distribution of the elements in the factors extracted according to PCA for loadings > 0.50, loadings > 0.75 in bold, and the total variance (%) explained by each factor.
Group 1Group 2
FactorElementsFactorElements
Factor 1
(50.9%)
Rb, Cr, Li, Ba, Ga, Be, V, W; ΣHREEs, Co, Mn, Mo, Nb, Tl, Zn.Factor 1
(81.9%)
Er, Tm, Y, Yb, Lu, Ho, Dy; Tb, Eu, Sm, La, Gd, Nd.
Factor 2
(23.7%)
Ta, ΣLREEs *, Ce, Hf, Zr, Bi, U; Nb, Tl.Factor 2
(11.6%)
Ce, Pr, Gd, Nd, Sm, La, Tb; Eu, Dy.
Factor 3
(21.4%)
Pb, As, Cd; Zn, Co, Mo.
* Sum of LREEs without Ce.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Soroaga, L.V.; Amarandei, C.; Negru, A.G.; Olariu, R.I.; Arsene, C. Assessment of the Anthropogenic Impact and Distribution of Potentially Toxic and Rare Earth Elements in Lake Sediments from North-Eastern Romania. Toxics 2022, 10, 242. https://0-doi-org.brum.beds.ac.uk/10.3390/toxics10050242

AMA Style

Soroaga LV, Amarandei C, Negru AG, Olariu RI, Arsene C. Assessment of the Anthropogenic Impact and Distribution of Potentially Toxic and Rare Earth Elements in Lake Sediments from North-Eastern Romania. Toxics. 2022; 10(5):242. https://0-doi-org.brum.beds.ac.uk/10.3390/toxics10050242

Chicago/Turabian Style

Soroaga, Laurentiu Valentin, Cornelia Amarandei, Alina Giorgiana Negru, Romeo Iulian Olariu, and Cecilia Arsene. 2022. "Assessment of the Anthropogenic Impact and Distribution of Potentially Toxic and Rare Earth Elements in Lake Sediments from North-Eastern Romania" Toxics 10, no. 5: 242. https://0-doi-org.brum.beds.ac.uk/10.3390/toxics10050242

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