Validation of the method for Nb, Ga, Ge, W, and Te determination in sediment
From ten TCE measured in this work, validation of the analytical method was earlier published for Li, Sc, Sb, Y, and La, which are more commonly measured in geological samples (Fiket et al. 2017). Here we present the results of analysis of Ga, Ge, Nb, Te and W in several different geological CRMs (soil, stream sediment and marine sediment), of which two are certified for all these elements and one for three of them (Table 2). For most of these elements, several isotopes and different instrumental resolutions are tested to define the best combination regarding accuracy and sensitivity. Namely, many of the isotopes of these elements are suffering from serious spectral interferences during ICP-MS measurement (Filella and Rodushkin 2018) and these authors provided a concise guide for the choice of isotopes and resolutions when using HR-ICPMS, and presented results of their measurements in CRM of river water. Following their recommendations, we made a choice of isotopes and resolution for analysis of Ga, Ge, Te, and Nb in geological CRMs (Table 2), while for W the most abundant isotope 184W, which should be free of major interferences, was used (Morisson et al. 2016). For Te analysis, in addition to 125Te and 126Te recommended by Fillela and Roduskhin (2018) and Yung et al (2013), we also tested isotope 128Te, which was recently used for the analysis of Te in freshwater sediments from various lakes in Canada (Wiklund et al. 2018). Recoveries for these elements in three different CRMs (Table 2) showed that for Ga both isotopes 69Ga and 71Ga can be used on MR, for Ge either 74Ge on MR or 72Ge on HR, while 93Nb can be measured both on LR and MR. The most demanding element is Te, because of its very low concentration in geological samples and isobaric interferences of Ba on the most abundant isotope 130Te. Yang et al. (2013) recommended analysis of 126Te in geological samples, as they found that 125Te suffer more from the interference of AgO+. Our results of CRMs analysis (Table 2) confirmed that 125Te is suffering from interferences (AgO+ and 85Rb40Ar+) on both LR and MR, especially in CRM for soil, where Te level is low (0.035 mg kg− 1). Isotope 128Te, which may suffer from the interference of 88Sr40Ar (Filella and Rodushkin 2018), gave too high recovery on LR for soil CRM, but could be used on MR for both CRMs. Analysis of different Te isotopes in our sediments confirmed that isotope 125Te gives higher concentrations (up to two times) than isotope 126Te. Isotope 128Te was suffering from the interferences on LR in the Plitvice sediments (giving 2–4 times higher concentrations than with 126Te) and on both LR and MR in sediments of the Lake Mir (giving 7 times higher concentrations than with 126Te), due to much higher level of Sr (2000–6000 mg kg− 1) in sediments from the Lake Mir than in the Plitvice Lakes (60–80 mg kg− 1). According to these results, we have chosen isotope 126Te on LR (which gives better LOD and RSD than 126Te on MR) for the determination of Te in our sediment samples. Detection limits for the determination of all analyzed TCE are given in Table 3. They were sufficiently low for the determination of low levels of these elements in the carbonate-rich sediments.
Distribution of TCE in the lake sediments
Concentrations (ranges, mean and median values) of measured TCE in sediments of the investigated lake systems are presented in Table 3 and vertical profiles of their concentration in sediment cores are displayed in Figs. 2, 3 and 4. Overall ranges of concentrations, which included all samples, are presented in Table 4, together with some literature data for the same elements in the marine and freshwater sediments. For comparison, we have used data for the coastal sediments (as we did not find data for marine lake sediments) from Japan (Otha et al. 2010) and the USA (Dolor et al. 2009), and freshwater lakes from Finland (Augustsson et al. 2010), China (Jin et al. 2006), Canada (Wiklund et al. 2018, 2020) and Macedonia (Vrhovnik et al. 2013). As we could not find more data for freshwater lake sediment for Ge, we used data on the long sediment core from the Baltic sea, which includes both older freshwater and more recent marine sediment (Ingri et al. 2014). Most of these areas are remote and not exposed to direct anthropogenic influence.
Table 4
Comparison of some literature data for TCE in marine and freshwater sediments with results from this work.
Element
|
Literature data, ranges and mean values (mg kg− 1)
|
This work (mg kg− 1)
|
Coastal sediment
|
Freshwater lake 1 (Baltic for Ge)
|
Freshwater lake 2
|
FOREGS (median)
|
All investigated lakes (ranges and median)
|
Li
|
5.2–68 (27)1
|
15–32 (23)3
|
26–49 (35)4
|
20.8 (Sed)
|
0.8–147 (18.4)
|
Nb
|
0.7–2.8 (4.7)1
|
2.1–4.3 (3.0)3
|
12–16 (13)4
|
13 (Sed)
|
0.19–17.4 (4.0)
|
Sc
|
1.6–67 (11)1
|
4.7–8.3 (6.8)3
|
9–13 (11)4
|
8.2 (Soil)
|
0.6–15.4 (3.2)
|
Ga
|
2.2–20.9 (12.3)1
|
5.4–10 (7.6)3
|
13–19 (16)4
|
12 (Sed)
|
0.19–19.8 (4.0)
|
Y
|
2.5–34 (14)1
|
12–27 (17)3
|
17–26 (22)4
|
25.7 (Sed)
|
0.42–41.7 (9.2)
|
La
|
2.9–31 (11)1
|
33–98 (50)3
|
36–45 (40)4
|
32.5 (Sed)
|
0.44–42.3 (10.6)
|
Sb
|
0.4–1.49 (1.2)2
(10)*
|
0.6–3.6 (1.7)5
|
0.06–1.2 (0.68)6
|
0.62 (Sed)
|
0.02–1.81 (0.45)
|
Ge
|
1.6–2.3 (2.0)2
(4.2)*
|
1.3–1.77
|
0.27–1.25 (0.65)6
|
1 (Soil)
|
0.01–1.5 (0.37)
|
W
|
1.11–1.89 (1.4)2
(3.7)*
|
1.9–2.4 (2.2)4
|
0.05–0.43 (0.12)6
|
1.24 (Sed)
|
0.07–2.6 (0.42)
|
Te
|
0.04–0.17 (0.1)2
(75)*
|
< 0.02–0.078
(10)*
|
< 0.02–0.420 (0.16)6
|
0.03 (Soil)
|
0.004-0.10 (0.031)
|
* – range includes locations more distant from pollution, in parentheses with *is maximum at polluted locations |
1 Otha et al. 2010; 2 Dolor et al. 2009; 3 Augustnsson et al. 2010; 4 Jin et al. 2006; 5 Vrhovnik et al. 2013; 6 Wiklund et al. 2020; 7 Ingri et al. 2014; 8 Wiklund et al. 2018. |
The ranges of concentrations for all TCE, when considering all samples from this study, are very wide, much wider than in the cited literature. The minimal concentrations of TCE are much lower (up to 10 times) than the minimum values in the literature. Median values for our sediments were also lower than the average (or median) values from the literature, which was used for comparison. They were also 2–3 times lower than FOREGS data for median values in the stream sediments in Europe (or soil for those elements for which data for sediments are not available), except for Y, Sb and Te, for which our median values were similar to those in FOREGS Atlas. Large ranges of TCE concentrations are a consequence of a great difference in sediments composition of the investigated lakes, from the almost entirely carbonate sediments in the Lake Kaluđerovac (Plitvice Lakes) to the clayey-rich material from the Visovac Lake (Table 1). For all TCE median and mean values are the highest in the Visovac Lake sediment (Table 3), comparing to the sediment from the Plitvice Lakes and the Lake Mir, which showed comparable average concentrations. The only exception is Y, which median concentration in Plitvice and Visovac sediments are similar.
Plitvice Lakes
Figure 2 displays vertical distribution of TCE concentrations in sediments of the Plitvice Lakes. For all TCEs concentrations decreased in order: P4 > P1 > K2 > Kal, in accordance with an increase of carbonate content and decrease of aluminum content in these sediment cores (Table 1). In sediments of the Lake Prošće (core P4) concentrations of all TCE were 20–40 times (depending on the element) higher than in sediments of the Lake Kaluđerovac. Variations in vertical profiles were similar to variations in Al profiles (Fig. S1) and were especially pronounced in the core P4. Core P4 is located close to the main tributary Matica which represents the most important sources of water and terrigenous material for the Plitvice Lakes system. Previous investigations have shown that surrounding source rocks (dolomites with small amounts of quartz and aluminosilicates) and soils are the main source of materials in sediments at locations P4, while the share of these minerals is low in sediments at location P1, where a significant amount of calcite was determined (Horvatinčić et al. 2017). Calcite was also the dominant mineral phase in the sediments where terrigenous input is small (the Lake Kozjak), or negligible (the Lake Kaluđerovac) (Horvatinčić et al. 2014, 2107). The morphological features of the calcite showed the structures typical for inorganically precipitated carbonates (Sondi et al. 2011), indicating that calcite in lakes is formed by active authigenic precipitation in water rather than by the erosion of surrounding carbonate rocks. Significant variations in the vertical profile of the core P4 relate to variations in the input of terrigenous material which depends on the flowrate of the tributaries. For example, the maximum at the 20–25 cm layers was a consequence of an extremely rainy period in spring 1981 (Horvatinčić et al. 2017).
The relationship between Al and metals can be used to distinguish between enrichment of geogenic and anthropogenic origin, as Al is of terrigenous origin and its concentration in lake sediments are not generally affected by anthropogenic influence (Roberts et al. 2019). Correlations of TCE with Al in all sediment cores are presented in Fig. 5. These correlations showed that sediments of the Lake Prošće have different geochemical composition than sediments from other Plitvice Lakes (Kozjak, Kaluđerovac), as well as the other investigated lakes, with proportionally lower content of Li, and much higher content of Y and partly La, in both P4 and P1 cores. This may be connected with weathering of surroundings soils or unique geological background of this lake, which was previously shown to be the main source of other trace elements, including some ecotoxic metals (Cd, Zn, Ni and Tl) to the Plitvice Lakes system (Dautović et al. 2014). However, numerous travertine barriers between the lakes, as well as efficient removal of the dissolved elements in the lentic parts of the system, mainly by co-precipitation with authigenic calcite and Mn oxides (Dautović et al. 2014), prevent the transport of elements to the lower lakes. For other TCE significant deviations from the general correlation with Al (Fig. 5) are not visible. In all four sediment cores elements Li, Nb, Sc, Ga, Y, and La were highly positively correlated with Al with correlation coefficient ( r ) being higher than 0.95, while this correlation was weaker for Sb, Ge, W and Te, with r ranging from 0.01–0.6 in cores P1, K2 and Kal, and 0.7–0.8 in the core P4. Weaker correlations between Al and elements Sb, Ge, W and Te may suggest a partly anthropogenic origin of these elements in the Plitvice Lakes sediments.
Lake Visovac
On Fig. 3 are shown vertical profiles of TCE concentrations in sediments of the Visovac Lake. For all TCE concentrations were similar in cores V1, V3, and V6 and the highest in the core V5. Variations in vertical profiles were similar to variations in Al profiles (Fig. S1) and were greatest in the core V5, where a pronounced maximum is visible in the central part of the core. Generally, concentrations of both Al and TCE are increasing with the sediment depth, indicating a continuous decrease of input of terrigenous material to the lake from the past to the present. Decreasing input of terrigenous material in the last decades, demonstrated also in other parts of the Krka River, was ascribed to advanced afforestation and changes in agricultural practices in the surrounding area (Cukrov et al. 2013). Vertical profiles of 137Cs, which was used for dating of V1, V3, and V5 sediment cores, indicated that V5 core was disturbed, and that the pronounced maximum in the central part should be a consequence of a massive input of surrounding soil to the lake in the past at this location (Mikac et al. 2017), so dating of this core was not possible. Higher concentrations of Al and all TCE in these layers are due to their higher concentrations in surrounding soils than in the lake sediment, where dilution by authigenic and biogenic carbonates decreases the metals content. The multi-element analysis of the same sediment cores (Mikac et al. 2017) showed partly different geochemical composition of material transported to the Visovac Lake by the rivers Krka (location V1) and Čikola (location V6), with the presence of dolomites in the drainage area of the Krka River. However, this difference was not visible in the relationship of measured TCEs with Al in cores V1 and V6 (Fig. 5). However, the concentrations of Nb and La were below TCE-Al line in the core V3, in the central part of the lake (Fig. 5), where the share of the fine authigenic calcite particles prevails, comparing to more coarse sediment, with lower share of carbonates, at the entrances to the lake (cores V1 and V6) (Mikac et al. 2017). This may indicate that dilution effect caused by carbonates has a greater influence on Nb and La compared to Al, which is more abundant in the finer fraction. In all four sediment cores very good positive correlations between Al and all TCE were obtained (r > 0.9), with little bit weaker correlations (r = 0.6–0.9) for Sb in the cores V1 and V3 and Te (r = 0.4) in the core V1.
Lake Mir
Figure 4 shows vertical profiles of TCE concentration in sediment core from the lake Mir. All measured TCE showed significant variations within the long sediment core, with RSD from 35 to 45 %. Variations of TCE along the core followed variations in the vertical profile of Al in the same core (Fig. S1). Previous investigations (Mlakar et al. 2015; Fiket et al. 2018) of metals distribution in the surface sediments of the lake Mir showed that the composition of the lake sediments is primarily influenced by natural processes: geological background dominated by carbonates and red soils, particle size distribution reflecting the bathymetry of the basin, and overall biological production (Mlakar et al. 2015). Comparison of REY (REE + Y) patterns in surrounding soils and lake sediment confirmed that the red soils were the main source of REY and other trace elements to the lake sediment (Fiket et al. 2018). The more recent investigation, showed that an additional important process is taking place in the lake Mir, the intensive authigenic precipitation of aragonite (I. Sondi, in the preparation according to the personal communication). Therefore, variations in the vertical profiles of TCE should be a consequence of the irregular input of terrigenous material into the lake (which is dependent on the amount of rainfall), in combination with seasonal authigenic precipitation of aragonite. Correlations with Al (Fig. 5) showed that Li, Nb, Sc, Ga, Y and La were highly positively correlated with Al (r > 0.95), correlation with Al was weaker for W (r = 0.55), Ge (r = 0.68) and Te (r = 0.55) and was not present for Sb (r = 0.09).
Evaluation of anthropogenic impact
Correlations of TCE with Al indicated that concentrations of some elements (Sb, W, Ge, Te) could be influenced by anthropogenic contribution. In order to evaluate the extent of possible anthropogenic impact, we calculated enrichment factors for TCE in the sediment cores for which background concentrations could be estimated, e.g., for which bottom sediment was deposited in the late 19th or early 20th century, before the major anthropogenic influence has started. Thus, enrichment factors were calculated for sediment cores P1, K2, and Kal from the Plitvice Lakes and core V3 from the Visovac Lake, in which the oldest sediments were deposited during that period (Table 1). EFs were not evaluated for two dated sediment cores, V1 and P4, in which, due to relatively high sedimentation rate, the bottom of the cores corresponded to the sediment deposited after 1950 (Table 1). However, we calculated EFs also for sediment core from the Lake Mir, in spite of the fact that it was not dated. Namely, as there is no direct input of material into the lake by streams, except occasional transport of surrounding soil with big rains, sedimentation rate should be low and not higher than 0.5 cm/year, making 90 cm long core at least 180 years old. Birch (2017) has defined the following EF thresholds: 1.0-1.5, 1.5-3.0, 3.0–5.0 and 5.0–10.0 to represent evidence of pristine conditions, minimal, moderate, and considerable enrichment, respectively. We used this scale to evaluate the degree of pollution in our sediment cores. In Table 5 are presented calculated EFs for all measures TCE in cores P1, K2, Kal, V3, and Mir. For elements Li, Nb, Sc, Ga, Y and La they ranged from 0.8 to 1.3 and confirmed that anthropogenic enrichment of these elements is not present in any of the investigated lakes. However, for other elements, maximal EFs were higher and EFs varied in the ranges 0.6–3.5 for Sb, 0.7–2.5 for Ge, 0.4-2.0 for W and 0.7-2.0 for Te, indicating the minimal or moderate anthropogenic influence on these TCE levels. Temporal trends of anthropogenic impact can be evaluated from EFs profiles of TCEs in the sediment cores (Fig. 6), where also EF profiles for Pb (from the same sediment cores, results will be published elsewhere) are presented for comparison.
Table 5
Ranges of enrichment factors (EF) in sediment cores (P1, K2, Kal, V3 and Mir) in which background concentrations could be determined.
Element
|
Range of EFs in sediment cores
|
Plitvice Lakes
|
Lake Visovac
|
Lake Mir
|
P1
|
K2
|
Kal
|
V3
|
MIR
|
Li
|
0.8-1.0
|
1.0-1.1
|
0.9–1.2
|
0.9-1.0
|
0.9–1.2
|
Nb
|
0.9–1.1
|
1.0-1.1
|
1.0-1.2
|
1.0-1.3
|
0.9–1.3
|
Sc
|
0.9–1.1
|
0.9–1.2
|
0.9-1.0
|
1.0–1.0
|
0.9–1.3
|
Ga
|
1.0-1.1
|
1.0–1.0
|
0.8–1.1
|
1.0–1.0
|
0.9–1.3
|
Y
|
1.0-1.3
|
0.8–1.1
|
0.9–1.2
|
1.0-1.2
|
0.9–1.2
|
La
|
1.0-1.1
|
0.9–1.1
|
0.9–1.2
|
1.0-1.2
|
0.9–1.3
|
Sb
|
0.9–1.9
|
0.7–2.5
|
0.9-2.0
|
0.8–1.5
|
0.6–3.5
|
Ge
|
0.8–1.5
|
0.7–1.4
|
0.8–1.5
|
0.9–1.2
|
0.7-2.0
|
W
|
0.9–1.5
|
0.7–1.4
|
0.8–1.4
|
0.9–1.3
|
0.4-2.0
|
Te
|
0.8–1.9
|
0.8–1.6
|
0.8–1.3
|
0.8–1.2
|
0.7-2.0
|
The most pronounced anthropogenic impact was found for Sb, for which EFs indicated, according to the scale of Birch (2017), minimal enrichment (EF = 1.5-3.0) in all 3 cores from the Plitvice Lakes (P1, K2, and Kal) and moderate enrichment (EF = 3.0–5.0) in the sediment core from the Mir Lake. The increase of anthropogenic impact on Sb level in the Plitvice Lakes started more intensively after the Second world war and reached maximum about 1990, after which it started to decrease. Vertical profiles of EFs for Sb followed closely profiles of EFs for Pb (Fig. 6), but EFs for Sb were much lower than those for Pb. Similar pollution trends of Pb and Sb were shown previously for the sediment of the lake Kozjak (Mikac et al. 2011) and were interpreted as a consequence of the use of leaded gasoline which started more intensively after 1950, peaked in the period 1980–1990 and was phased out after 2000, resulting in the emission of alkyl-lead compounds into the atmosphere and significant contamination of surface waters by lead and organolead compounds in Croatia (Mikac and Branica 1994a; Mikac and Branica 1994b). Smelting of lead sulfide ores, which are rarely monomineralic, also releases other elements, and lead pollution in atmospheric particles is regularly accompanied by other trace elements like Zn, Cu, Sb, Bi, Ag, Tl (Krachler et al. 2008). Due to its many commercial uses, especially those in automotive brakes, Sb was found to be one of the most enriched elements in the urban environment and one of the most heavily influenced by anthropogenic activity (Fillela et al. 2002; Sen et al. 2012; Fillela and He 2020; Wiklund et al. 2020). The highest enrichment for Sb was found in sediments of the Lake Mir, where an increase of EFs starts at a depth of about 30 cm, continues to the dept of 10–15 cm and then decreases toward the surface, closely following the profile for Pb (Fig. 6). To define when these pollutions started, sedimentation rate in the lake should be determined. The close relationship between Pb and Sb suggests that the sources of pollution for these elements are the same, and possibly related to the road traffic. However, more detailed interpretation of the anthropogenic impact on the Lake Mir needs further research, since this lake is located on the island separated from the urban areas. In spite of the evident anthropogenic enrichment of Sb, it should be emphasized that maximal Sb concentrations are low in comparison with existing sediment quality criteria for this element (Maccready et al. 2006) which define 2 mg kg− 1 for Sb as a limit under which adverse effects are rarely expected.
The next element which showed maximal EFs between 1.5-2.0 (minimal anthropogenic impact) in both Plitvice and Mir sediments was Te (Table 5, Fig. 6), with the maximal EFs recorded in similar layers (around 1990) as for Sb and Pb (Fig. 6). In a recent study of atmospheric tellurium contamination of numerous lakes in Canada (Wiklud et al. 2018) it was shown that the main sources of anthropogenic Te were metal (copper) smelters and coal mining/combustion facilities and that long-range atmospheric transport (more than 1000 km) and deposition of Te to remote lakes were significant. Taking into account this knowledge, we may presume that, in addition to coal combustion, the mining and smelting of copper in the town of Bor in Eastern Serbia (500 km from the Plitvice Lakes), which is possibly the most important mining and metallurgical center in southeast Europe (Dimitrijević et al. 2009; Urošević et al. 2018), could be a distant source of Te. Mining of the copper ores in Bor started in 1903, but smelting began only in 1935, and environmental protection measures were introduced first in 2015 (Urošević et al. 2018). As Te is rarely measured, there are no data on Te levels in this area, but severe pollution of air with S and As (Dimitrijević et al. 2009) and soil with Cu and Pb (Dimitrijević et al. 2015) was found.
Elements Ge and W did not show anthropogenic enrichment in the Plitvice Lakes (Table 5, Fig. 6), but demonstrated EFs between 1.5 and 2 in the upper 10 cm of the sediment of the lake Mir, similarly as for Te. As it was discussed previously, to interpret these results and define possible contamination sources of these elements further geochemical study of this lake is needed.
Any of TCE which is prone to anthropogenic influence (Sb, Ge, Te, W) did not show elevated EFs in sediments of the Lake Visovac (Table 5, Fig. 6). In sediments of this lake low anthropogenic input of Pb, Cr, and Cu by the Čikola River was demonstrated, as well as significant enrichment in deeper sediment layers with Cd and Zn in the past (Mikac et al. 2017), brought to the Visovac Lake by the Krka River from the zinc electroplating facilities located upstream in the area of the town of Knin (Cukrov et al. 2008; Filipović Marijić et al. 2018). Obviously, these contamination sources did not result in enrichment with elements Sb, Ge, Te and W.
Although for some TCE anthropogenic influence is demonstrated, it is very low and concentrations of all these elements are much lower than in similar studies from the literature, proving that the investigated lakes can be considered uncontaminated with these elements. This is in accordance with previous studies which showed that the investigated lakes are unpolluted regarding other trace elements like Cd, Pb, Cu, and Zn (Mikac et al. 2011; Cukrov et al. 2008; Mlakar et al. 2015).