Critically raw materials as potential emerging environmental contaminants, their distribution patterns, risks and behaviour in floodplain soils contaminated by heavy metals
HomeHome > Blog > Critically raw materials as potential emerging environmental contaminants, their distribution patterns, risks and behaviour in floodplain soils contaminated by heavy metals

Critically raw materials as potential emerging environmental contaminants, their distribution patterns, risks and behaviour in floodplain soils contaminated by heavy metals

Jun 20, 2023

Scientific Reports volume 13, Article number: 9597 (2023) Cite this article

686 Accesses

Metrics details

The expanding demand for new critical raw materials can lead to their increased release to the environment in the form of emerging environmental contaminants (EECs). However, there has never been a comprehensive study that takes into account the total EEC content, the content of various EEC fractions, their behaviour in floodplain soils, and potential ecological and human health risks. The occurrence, fractions, and influencing factors of the seven EECs (Li, Be, Sr, Ba, V, B, Se) originating from historical mining in floodplain soils of various ecosystems (arable lands, grasslands, riparian zones, contaminated sites) were investigated. Based on the evaluation of the overall levels of EECs (potentially toxic elements) in comparison to the soil guideline values set by European legislation for Be, Ba, V, B, and Se, it was found that only Be did not exceed the recommended limits. Among the elements analyzed, Li had the highest average contamination factor (CF) of 5.8, followed by Ba with 1.5 and B with 1.4. Particularly concerning was the discovery of a potential serious health risk associated with Li exposure for children, as indicated by hazard quotients ranging from 0.128 to 1.478. With the exception of Be and Se, the partitioning of the EECs into the different fractions revealed that the EECs are primarily bound with the residual fraction. Be (13.8%) had the highest percentage of exchangeable fraction as the most bioavailable in the first soil layer, followed by Sr (10.9%), Se (10.2%), Ba (10.0%), and B (2.9%). The most frequently observed correlations were between EEC fractions and pH/KCl, followed by soil organic carbon and manganese hydrous oxides. Variance analyses confirmed the impact of different ecosystems on EEC total content and fractions.

Globally, the demand for new critical raw materials (CRMs) is growing very quickly. The 2020 EU list contains 30 materials including lithium, beryllium, strontium, vanadium, baryte (barium sulphate) and borates (naturally occurring minerals containing boron)1. Selenium was identified as a critical material in the 2011 assessment2. This growing global importance and demand for CRMs can lead to undesirable increased in their release to the environment. Therefore, many of these CRMs have also attracted attention from the scientific community as emerging environmental contaminants (EECs). Environmental contamination of these elements can be related to their extensive use, e.g. Li in mobile phones and mood-stabilizing drugs3; Be in electronics, telecommunication, nuclear power and industry4; barium compounds are used as fillers or additives in industrial products; vanadium is used as an additive in steel and titanium alloys and as a catalyst for chemicals; borates are important ingredients in a variety of household and commercial products. To date, research has not succeeded to provide a clear picture of how increased production of CRMs as EECs adversely affects all compartments of the environment3. There is a wide variety of inorganic and organic EECs5.

The release of industrial inorganic EECs into the environment can also be related to mining activities that are mainly responsible for the release of traditional environmental contaminants (TECs). EECs that were not recognized before or attracted less attention do not exert acute toxicity but exert their effects in more hidden ways6. Additionally, increased use and disposal of CRMs can contribute to increased exposure of EECs. Drinking water and plants often provide the primary pathway of exposure for many EECs to the food chain7. Both depend on the behaviour of EECs in soil which can be entirely different from TECs6. There is a lack of information on the distribution patterns and behaviour of many EECs in the soil. The areas suitable for this research are those in which past mining activities have left negative environmental footprints related to soil contamination by heavy metals. Heavy metal accessories can simply be EECs with different chemical properties. In this regard, special attention should be paid to floodplain soils. In floodplains, contamination can occur from both geogenic and anthropogenic sources. Additionally, drainage systems and water flows can contribute to the transport and spread of contamination. Furthermore, contamination can negatively affect even distant rare natural and artificial alluvial ecosystems with a possible negative impact on human health.

The purpose of this research is to look into the distribution patterns and behavior of selected EECs (Li, Be, Sr, Ba, Ti, V, B, and Se) in the Štiavnica River (SR) floodplain soils that have been contaminated by heavy metals from historical mining. In this study, we aimed at: (1) to specify the total content and the different fractions of EECs in soils; (2) to calculate soil contamination indices; (3) to assess human health risk; (4) to evaluate the distribution patterns of EECs; (5) to find the relationships between soil characteristics and the EEC fractions; (6) to assess the impact of different factors on the behaviour of EECs in soil. Our study is novel because selected EECs have not yet been studied in such an intensive and thorough way in floodplain soils, taking into account the total content, the various fractions, their relationships in floodplain soils, and potential ecological and human health risks.

Sampling sites were selected in the floodplain along the SR in central Slovakia. The riverbanks consist of fluvial deposits (sand and silt). The SR is 55 km long, with a 443 km2 large catchment area. The Štiavnické vrchy Mountains, Podunajská pahorkatina Hills, and Krupinská planina Plain are the main physiographic divisions that characterize this drainage basin. In the upstream part of the basin in the Štiavnické vrchy Mounatins of volcanogenic origin, the relief is dominated by hills and valleys with forests and partly by permanent grasslands used as meadows. The Štiavnické vrchy Mountains are a type of neovulcanic mountain range. They consist of pyroclastic rocks, primarily andesite tuff agglomerates, andesites and rhyolites. The region hosts rich metal, silver and gold mineralizations that have been the basis for a long lasting mining tradition. In the eighteenth century, the region became the biggest mining centre in the Habsburg Monarchy8. In the downstream part of the basin in the Podunajská pahorkatina Hills, the relief is dominated by a flat surface. The bedrock consists of clay- and sand-stones on slopes, sand and gravel in areas of floodplains. The land is used as arable land.

We collected a total of 42 soil samples from two layers (0–10 and 20–30 cm) from two contaminated sites (CS, classified as Ekranic Technosols), 1 reference site (RS, classified as Haplic Cambisol), and 12 alluvial sites (AS, classified as Haplic Fluvisols). The 12 AS were located along the continuous elevation gradient of the SR, from the highest elevation (659 m a.s.l.) down to an elevation of 123 m a.s.l. (Fig. 1, Table 1).

Isometric map of the sampling points and photos from selected study sites (photos R. Kanianska). The map was generated using 3D Map Generator, ver. 1.5, https://www.3d-map-generator.com/3d-map-generator-terrain/.

The sampling points were located on the river bank within the riparian zones (RZ) in the Štiavnické vrchy Mountains, and on arable land (AL) in the Podunajská pahorkatina Hills. In the case of a larger floodplain, we also collected samples distant from the water flow at the end of the floodplain, that were used as permanent grasslands (GL) in the Štiavnické vrchy Mountains, and as arable land in the Podunajská pahorkatina Hills (at 5 of 12 alluvial sites). Outside the river floodplain, we selected the two contaminated sites (CS), which are poorly developed on metallurgical slag and shallow soils. As CS sites, we selected the area of the Maximilian shaft that is registered in the Register of environmental burdens of the Slovak republic as probable environmental burden (PEB), and the hard rock mine tailing Lintich that is registered as a confirmed environmental burden (CEB), where hazardous substances caused by mining activities can pose a risk to human health or to the environment, soil and groundwater. The Maximilian shaft area was located approximately 1 km and the Lintich mine tailing approximately 100 m from the Štiavnica River.

The RS site represents a site with a ‘normal’ level of contaminants and arising from a combination of natural and diffuse pollution contributions. The RS site was used to determine the background values. The RS site was located about 400 m away from the Štiavnica River outside the floodplain on a slope and is used as a meadow. The soil type was classified as Haplic Cambisol.

Soil samples were collected in the study area in September 2018. Five subsamples were taken with a stainless-steel blade, and then mixed thoroughly to get a composite soil sample. The collected soil samples were stored in a polyethylene bag, transported to the laboratory where they were air-dried at ambient temperature, crushed, and sieved through a 2-mm mesh.

The volumetric soil moisture (SM in %) and salinity (SAL in mS/m) were measured at the same time and places of soil samples collection with WET Sensor (type HH2, Delta-T Devices Ltd., Cambridge, UK), by taking the average of five measurements. Similarly, the redox potential (Eh) was measured by the Oxygen Diffusion meter (Ejkelkamp).

The soil chemical parameters were measured in November and December 2018. We determined soil chemical properties such as pH in H2O (with a ratio of 1:2.5) and 1 M KCl solution (with a ratio of 1:2.5), total soil organic carbon (SOC) using the oxidimetric method according to Tyurin that is similar to the Walkley–Black oxidation method. The Tyurin titrimetric method is a wet combustion method, where soil organic matter is oxidized by 0.2 M K2Cr2O7 with H2SO4. After oxidation, excess dichromate is determined by titration with a Mohr´s salt solution. The content of humic and fulvic acids was determined by the short fractionation method according to Kononova and Belchikova using extraction with 0.1 M Na4P2O7·10H2O. Oxalate-extractable Al, Fe, and Mn as an indicator of hydrous oxides, were determined by the acid oxalate extraction method as described by Pansu and Gautheyrou with a (NH4)2-oxalate solution adjusted to pH 3, shake for 4 h in the dark, and centrifuge for 10 min at 10,000g. Particle size analysis was performed by pipette method using soil particle sedimentation and particle size fractions (sand, silt, clay) were classified according to the United States Department of Agriculture (USDA) system.

We investigated the total content and content in four chemical forms of seven EECs from various periodic table groups (Li, an alkali metal; Be, Sr, Ba, alkaline earth metals; V, a transition metal; B, metalloid; Se, a nonmetal). Total dissolution of soil samples was performed using a microwave unit (Anton-Paar, Austria). The dissolution method used was an EPA 3015A9: 9 mL of HNO3, 3 mL of HCl and 1 mL of HF were added to a 0.250-g dry soil sample.

The Community Bureau of Reference Sequential Extraction Procedure (BCR SEP) as described by Rauret et al.10 was used. The exchangeable phase (I) was leached with HOAc, the reducible fraction (II) bound to Fe and Mn oxides was leached with NH2OH-HCl, the oxidisable fraction (III) bound to organic matter and sulphides was leached with H2O2 in NH4OAc, and the residual fraction (IV) bound to silicate minerals was calculated as the difference between the total content of EECs and the sum of the three fractions. The procedure followed was inductive coupled plasma optical emission spectrometry (ICP-OES).

We evaluated soil contamination by EECs using the contamination factor (CF) according to Varol11 and the pollution load index (PLI) according to Tomlinson et al.12. The CF was calculated according to the formula:

where Ci is the measured content of EEC in soil and Cb is the background value of the EEC content of the reference site. According to Varol11, CF < 1 means low contamination, 1–3 is moderate contamination, 3–6 is considerable contamination, and CF > 6 means very high contamination.

The PLI was calculated according to the formula:

where n is the number of elements. The PLI > 1 indicates pollution.

Potential human health risk was also evaluated. In the study, ingestion of soils contaminated with EECs was considered as the main pathways for risk assessment. We calculated the average daily dose (ADD, mg element kg−1 bodyweight day−1) for children, adult male and adult females:

where C is total soil element content; IR is soil ingestion rate (children: 200; adults: 100 mg dust day−1); EF is exposure frequency (children: 350; adults: 250 days year−1); ED is exposure duration (children: 6 years; adults: 25 years); BW is bodyweight (children: 15 kg; adults males: 68 kg; adults females: 58 kg); AT is averaging time (children: six 365-day years = 2190 days; adults: 9125 days). Values were obtained as per by Rinklebe et al.13, as originally reported in USEPA14.

Then hazard quotient (HQ) was calculated to evaluate non-carcinogenic human health risk:

where RfD (in mg element kg−1 bodyweight day−1) is the EEC oral reference dose: for Be = 0.002, Ba = 0.07, V = 0.007, Se = 0.00514, Sr = 0.615, B = 0.216. The RfD for lithium is not available on the Integrated Risk Information System of USEPA. However, the USEPA has developed a provisional peer reviewed toxicity value for lithium with provisional p-RfD of 0.002. The provisional subchronic and chronic RfD for lithium was derived from the lowest observed adverse effect level (LOAEL) from the therapeutically used lithium of 2.1 mg/kg-day for adverse effects in several organs and systems. Dividing the LOAEL of 2.1 mg/kg-day by an uncertainty factor of 1000 yields a subchronic and chronic pRfD of 0.002 mg/kg-day17. We think that this value can be used for health risk assessment of soil because exposure through environmental compartments should be minimal to zero.

Then we calculated hazard index (HI) as the sum of HQ values of all EECs. Values of HQ and HI of > 1 indicate high probability of the occurrence of adverse health effects.

The study analyzed the data using the Pearson correlations, principal component analysis (PCA), one-way analysis of variance (ANOVA), and hierarchical cluster analysis (HCA). The SPSS Statistics 28 software was used.

In Slovakia, it is allowed to collect soil samples for research purposes.

The basic soil properties are reported in Table 2.

The total content of EECs in the two soil layers of the SR floodplain together with soil guideline values (SGV) as a range in European legislation18 is presented in Table 3. At least one lower SGV limit set up for 5 selected EECs (Be, Ba, V, B, Se) was exceeded in the case of all selected EECs with the exception of Be. Some exceedances were recorded in the CS, AS and RS localities, in both layers, and in all ecosystem types (CS, AL, GL and RZ). For Ba and Se, the upper SGV limit was even exceeded, but only in one locality. The contamination risk in the CS and AS localities was confirmed by the CF and PLI along the whole SR floodplain for 5 EECs (Li, Sr, Ba, V, B) (Table 4). The CF and PLI were not calculated for 2 EECs (Be, Se) due to the appearance of zero values. We measured the lowest content of these two elements in soil with no SGV exceedance in the case of Be and only a few exceedances in the case of Se. In the CS localities, the average CF showed moderate contamination by Li. At AS localities, the average CF showed considerable contamination by Li, moderate contamination by Ba and B. The highest CF was quantified for Li (15.6 in 13AS-AL2), Ba (2.6 in 15AS-AL2) and B (2.3 in 8AS-RZ). In many AS localities, contamination by EECs was similar to or even higher than in CS localities. The PLI indicated that all AS sites were contaminated by EECs (PLI > 1) except one AS locality (3AS-RZ).

The potential health risk was assessed by considering oral topsoil ingestion for 7 EECs (Table 5). We discovered possible serious potential health risks only in children, a vulnerable population. The two adult groups are not under the potential health risk by EECs. In the two adult groups, the mean value of HQ for individual EECs ranged from 0.000 to 0.043 for males, and from 0.000 to 0.050 for females. The resultant mean HI values for female persons were 0.040 at CS localities and 0.080 at AS localities. The mean HI values for male persons were 0.034 at CS localities and 0.068 at AS localities. The mean value of HQ for individual EECs ranged from 0.002 to 0.545 for children. The resultant mean HI value for children was 0.435 at CS localities and 0.868 at AS localities. In children, the mean HQ (0.545) at AS localities was found high for Li. Lithium was followed by V (0.238), and Ba (0.070); all the rest of the EECs were much lower and ranged from 0.002 (Sr, B) to 0.005 (Be, Se). The risk for children indicated by HQ for Li was significant at 5 of 15 AS localities where HQ was found high (ranging from 1.085 to 1.478). Two of these 5 localities are used as grasslands, two as arable land and one is a riparian zone.

The content of seven EECs leached in the sequential extraction procedure in the CS, AS, and RS localities and in both depths is summarized in Fig. 2a–e. We recorded different distribution patterns among all four fractions between different elements and related chemical groups of EECs. Between the two different layers among the AS and RS sites, more similar patterns occurred at the AS sites.

Emerging environmental contaminants distribution in 4 fractions in the depth of 0–10 and 20–30 cm in contaminated, alluvial, and reference sites. I—exchangeable fraction, II—reducible fraction, III—oxidisable fraction, IV—residual fraction.

V and B showed the highest percentage in the residual fraction (IV) in all CS, AS, and RS sties, and in both depths. This fraction may be considered as biologically unavailable. On the other hand, among all the elements observed, only one element, Sr, showed the highest percentage in the first extractable fraction, which was considered the most mobile in the RS site at the second depth. The third (oxidisable) fraction occurred in the highest percentage in the cases of Li (at CS and in the first depth at RS), Ba (in the first layer at AS and in the second layer at RS), and Se (in the first and second depths at RS). Be appeared in three different fractions (II, III, IV) in the highest percentage at different sites, followed by Li, Sr, Ba, and Se in two different fractions in the highest percentage at different sites. V and B occurred in only one (IV) fraction of the highest percentage at all sites.

Alkali earth metals (Be, Sr, and Ba) reflected the most equable distribution between fractions together with B and Se in some cases. The distribution of the alkali earth metals in the different fractions mostly associated with the first two most mobile fractions (I and II) indicated that the mobility declined on the order of Ba, Sr, and Be, at CS sites; Ba, Be, and Sr at AS sites; Be, Sr, and Ba at the RS site in the first depth and Sr, Be, and Ba in the second depth, respectively.

At AS localities in the first depth, the average content of Ba represented 10.0% in exchangeable, 22.5% in reducible, and 9.6% in oxidisable fraction; the content of V represented 0.2% in exchangeable, 2.7% in reducible, and 8.8% in oxidisable fraction; the content of B represented 2.9% in exchangeable, 1.6% in reducible, and 7.8% in oxidisable fraction; and the content of Se represented the same.

The correlation coefficient matrix for the EEC fractions with the chemical and physical parameters of the soil is shown in Table 6. The most numerous significant correlations with chemical properties occurred in the case of B (17), while the least occurred in the case of Li (6). The most numerous significant correlations with physical properties occurred in the case of Sr (12) while no correlations occurred in the case of Se (0). Among chemical properties, we most often observed correlations with pH/KCl (18 out of 35), followed by SOC (14 out of 35), Mno (13 out of 35), pH/H2O (10 out of 35), HA/FA (10 out of 35), Feo (5 out of 35), Eh (4 out of 35) and Alo (3 out of 35). Only in one case of boron, there was correlations between all four fractions, including the total content and the oxalate-extractable Mno. Among physical properties, we most typically observed correlations with the silt and sand content (9 out of 35), followed by the clay content and salinity (7 out of 35).

The PCA showed that the first two factors can explain more than 43% of the total variance for all EECs, from 43.2% in the case of Be to 50.9 in the case of B (Fig. 3). Salinity and pH/H2O were positively and Eh negatively correlated with PC2 for Li and Ba. pH/KCl was clearly positively correlated with PC1 for Be. pH was positively correlated with PC2 and silt was negatively correlated with PC1 for Sr. Soil moisture was positively correlated with PC2 and pH negatively correlated with PC1 for V. pH was clearly strongly positively correlated with PC1 and Eh with PC2 for B and Se.

PCA relationships between EECs and soil properties. I—exchangeable fraction, II—reducible fraction, III—oxidisable fraction, IV—residual fraction, T—total content.

To evaluate the significance of the difference between mean values of each EEC fraction, one-way ANOVA was used. The independent variables were two different ecosystem types (riparian zone and agricultural ecosystem represented by arable land and grassland) and soil depths (0–10 and 20–30 cm), while the dependent variables were the EEC fractions. The result of ANOVA test is presented in Table 7. The statistical analysis showed significant differences between the ecosystem types and all Ba fractions, including the total content. Significant differences appeared between the different ecosystem types for Li and Sr in fractions III, IV and total content, for B and Se in fractions II and III, and for V in fraction I. The ANOVA showed significant differences between the soil depth and the EEC fractions only in the case of Li (I).

The HCA (Fig. 4) has highlighted two large clusters (A and B). Cluster A contains sites that are primarily used as agricultural land (GL and AL). Two clusters of lower hierarchical order are recognised in B (B1 and B2). Cluster B1 includes predominantly RZ sites (16 out of 19). Cluster B2 includes five sites, two of which are contaminated sites.

Dendrogram for the hierarchical clustering of sites (Ward’s method).

Based on the available SGVs, we discovered four EECs (Ba, V, B, Se) contaminating the entire SR floodplain. The SGV value exceeded was not only reported in the case of Be. The SGVs were not available for Li and Sr. According to the calculated contamination indices for five EECs (Li, Sr, Ba, V, and B), Li is the most serious EEC of the SR floodplain, followed by Ba and B. Based on the calculated contamination indices for 5 EECs (Li, Sr, Ba, V, B), mainly Li and fewer Ba, B are the most serious EECs of the SR floodplain. Based on the content of selected EECs in soil and calculated contamination and pollution indices, many AS localities are more contaminated by selected EECs than the confirmed environmental burden (4CS-CEB). In addition, some of them are used intensively for agricultural production. Although the contents of the selected EECs exceeded some SGVs, the contamination indices indicate mainly moderate contamination. Therefore, they are probably not of interest for further industrial use in the sense of soil mining. The demand for new raw materials is growing very quickly, and other raw or natural materials are already being used today for the extraction of CRM. In recent times, many CRMs are seldom mined as the main ore. They are usually coproducts of the principal ore. Therefore, historical mining regions may be of such interest19. Many tailings did not meet the requirements of further industrial production and were only piled up as waste20. In addition, tailing deposits have become potential contamination source21. Thus, they impacted the surrounding environmental compartments. The proof of this is the SR floodplain contaminated by heavy metals as reported by Kanianska et al.22 and other pollutants. The HMs together with EECs released during mining processes were transported and spread along the floodplain. Rivers and alluvial soils are affected by discharges from mines and tailings23. In addition, polluted aquifers act as long-term pollution source to the surrounding gaining rivers and soils, even after the mining activities were stopped24. Despite the considerable scale, there are many such territories with limited data availability25. However, they can pose ecological and health risks with respect to the cumulative effect of possible additional contamination from existing and potential new pollution sources. In our study, Li exhibited the highest mean contamination factor with a value of 5.8, followed by Ba with 1.5 and B with 1.4. We also found possible potential health risk for children in the case of Li at 5 of 15 alluvial localities with HQ higher than 1.

Different contaminants with different chemical properties can contribute to these risks and interact with each other under the simultaneous effect of the soil environment and its specific properties. Therefore, research on the most serious EEC behaviour in soil in relation to soil properties is one of the key factors to prevent or decrease ecological and health risks. The BCR SEP method allowed us to study the mobility and availability of EECs. This method is accepted by many researchers dealing with toxic metals in mining waste (e.g., Fernandez-Ondono et al.26). Although there are authors who believe that the application of this procedure can lead to ambiguous and misleading interpretation27. This method is still being increasingly applied in EEC analyses (e.g. Chen et al.28, Mikavica et al.29, Xiong et al.30).

Lithium does not belong to the essential nutrients group. But there is evidence that low Li levels have beneficial effects on living organisms31. Lithium possesses unique essential and neurophysiological characteristics in biota and humans32. However, at high doses, it can be toxic to humans33.

The Li content in the soil is variable from place to place. Upper layers usually contain less Li than the underlying layers34, as markedly the case at the RS locality. In the AS localities, the total Li content was mostly higher in the first depth, but at 7 out of 19 localities it was the opposite. The average Li content in European agricultural soils is 11.4 mg/kg, with a predominantly low concentration (median 6.4 mg/kg) in northern Europe and higher values (median 15 mg/kg) in southern Europe35. We recorded a total Li content of 14.78 mg/kg in the RS locality and 29.29 mg/kg as a mean value in the CS localities. A significantly higher total Li content was recorded in AS localities used in agriculture with a mean value of 152.52 mg/kg compared to localities in the riparian zone with a mean value of 36.46 mg/kg. Such content can pose an environmental risk. In the CS localities, the average CF (2.0) showed moderate contamination by Li. At AS localities, the average CF (5.8) showed considerable contamination by Li.

We also discovered serious potential health risks for children, a vulnerable population. The mean HQ (0.545) at AS localities was found to be high for Li. The risk for children indicated by HQ for Li was significant at 5 of 15 AS localities where HQ was found high (ranging from 1.085 to 1.478). Lithium toxicity results in abnormalities and dysfunctions in several metabolic pathways, causing serious problems for human health3.

The toxicity of EECs is also related to the distribution of their fractions. Li was heterogeneously distributed in fractions among different ecosystem types. Exchangeable and reducible fractions showed up in all land use types and in both layers with the lowest Li content, accounting for ˂ 1% Li, which is in line with the findings of other authors36. The residual fraction demonstrated the highest Li content at the AS localities in both depths and at the RS locality in the deeper layer. At the CS and in the first depth at the RS locality, Li demonstrated the highest content in the oxidisable fraction. This was confirmed by correlation analyses, where we found positive correlations between the reducible and oxidizable Li fractions and the SOC, similarly as with pH/KCl. Motowicka-Terelek et al.37 also found correlations between soil Li content and SOC, but there are also authors who did not36,38. From the physical properties, the fine granulometric fraction (silt) showed a positive correlation with the residual fraction and the total Li content. Kabata-Pendias and Mukherjee39 refer to the texture, respectively, the fine clay granulometric fraction as the most important factors controlling the Li status of soil. They consider SOC and pH to be factors of lesser importance.

Barium is not an essential nutrient for humans. But it probably fulfils biochemical functions important for animals and plants40. The content of barium in food does not usually pose a health risk. Its long-term impact on human health is questionable41. Barium toxicity is related to the solubility of the compound42.

The general threshold and guideline values for Ba have not been developed. Only a few countries in the world have adopted them. The Ba limit for food safety has not yet been established43. Furthermore, geographic variation is little known44. Anthropogenic releases are primarily associated with industrial activities, including mining, and pesticide application in agriculture45. Barium compounds are quickly dissolved and can be widely spread, as was also the case with the SR.

In AS localities, the total Ba content was mostly higher in the first depth, but at 6 out of 19 localities, it was the opposite. The median content of Ba in European agricultural soil is 62 mg/kg, with predominantly low concentrations (median 43 mg/kg) in Northern Europe and higher values (median 74 mg/kg) in Southern Europe18. We recorded a total Ba content of 255.34 mg/kg in the RS locality and 66.91 mg/kg as a mean value in the CS localities. A higher total Ba content was recorded in AS localities with a mean value of 384.33 mg/kg (median value of 394.47 mg/kg), ranging from 108.34 to 659.82 mg/kg in the first soil layer. González-Valoys et al.46 measured in the mine soils in Panama an average total Ba content of 200.6, and in the surrounding soils 349.7 mg/kg. Norini et al.43 measured in the mine soil of a settling pond of the Ag–Pb ore mine in France 1701 mg/kg of total Ba and 27.2 mg/kg of water-soluble and cation exchange Ba fraction. At the AS localities, we measured the mean content of Ba in the I exchangeable fraction of 38.58 in the upper layer, ranging from 10.27 to 79.77 mg/kg.

The oxidisable fraction demonstrated the lowest Ba content in all types of land use and in both layers, accounting for 4.14–10.81% Ba. The Ba is thought to be incapable of forming organic ligands or complexes with fulvic and humic acids47. This was also confirmed by correlation analysis, where we did not find correlations between the Ba fractions and the SOC. We recorded only negative correlations between the oxidisable Ba fraction and the HA/FA ratio. In the AS and RS localities, the residual fraction showed the highest content, accounting for 51.7 to 85.44% Ba. Ba was heterogeneously distributed in extractable and reducible fractions between different types of land use and soil layers. The reducible and extractable fractions showed the highest Ba content in CS localities. The extractable fraction demonstrated on average 10.0% Ba at AS localities in the first layer and 11.6% in the second layer. Cappuyns44 found low Ba mobility in alluvial soils of the Schelter and Demer rivers in Flanders (˂ 3% of the total content). Ba adsorbs into Fe hydroxides48, and we also found a positive correlation between the reducible Ba fraction and Feo. The pH of the soil can influence barium behaviour and is more mobile under acidic conditions. The pH showed negative correlations with the total Ba content, oxidisable fractions, and residual fractions. In contrast to Frohne et al.49, we did not find correlations with Eh. Similarly, we did not find correlations with clay particles, as was found by Atun and Bascetin50. The reducible and oxidisable Ba fractions were positively correlated with the content of the sand particles. Only the IV Ba fraction was positively correlated with the silt content and negatively correlated with the content of the sand particles.

Boron at low concentrations is probably an essential element for plants51 humans and animals, while it is toxic at high concentrations52. But recently, B was assumed to be rather toxic to plant cells at low levels53. Boron has an extremely narrow range between deficiency and toxicity54.

In AS localities, the total B content was mostly higher in the first depth, but at 6 out of 19 localities it was the opposite. The median content of B in European agricultural soil is 2.42 mg/kg, with a predominantly low concentration (median 1.7 mg/kg) in northern Europe and higher values (median 3 mg/kg) in southern Europe18. Near a large boron mine in Turkey, Alan and Kara55 measured the total B content in soil from 36.5 to 820.5 mg/kg. We recorded a total B content of 27.49 mg/kg at the RS locality and 25.56 mg/kg as the mean value at the CS localities. In AS localities, the mean value was 37.23 mg/kg (median value 32.46 mg/kg), ranging from 20.27 to 64.57 mg/kg in the first soil layer.

Although the total B content exceeds the soil guideline value in most AS localities, sequential extraction showed the highest percentage of residual fraction (IV) in all CS, AS and RS sties, and at both depths. On the other hand, the percentage of mobile and other fractions was low (2.9% in exchangeable, 1.6% in reducible and 7.8% in oxidisable fraction). Similarly, Alan and Kara55 found the highest B content in the residual fraction. Therefore, most of the boron that occurs in the soil is unavailable to plants56.

We found the most correlations between the B fractions and the physicochemical properties. However, the Eh was negatively correlated with extractable and reducible fractions of B. The pH, SOC, and Mno oxides were positively correlated with the extractable, reducible, and oxidisable fractions of B. The extractable B fraction was also positively correlated with HA/FA. Hou et al.57 also found positive correlations between readily soluble B and pH and between the oxidisable fraction and SOC. From the physical properties, the extractable fraction correlated positively with the clay content of the soil in the study area.

We detected some level of EEC contamination, particularly by Li, Ba, and B in the floodplain soils of the SR. This contamination could originate simultaneously with HM contamination, mainly from historical mining. The downstream transport could be a source of EECs as accessories to HMs for distant soils.

Sequential extraction technique was used to assess the environmental risk of EECs. The partitioning of the EECs showed that the EECs are mainly associated with the most stable (biologically unavailable) residual fraction, except for Be and Se. The risk of these two elements is reduced by their low (even zero) total contents in the soil. The exchangeable fraction as the most bioavailable occurred in the first soil layer at the AS localities in the highest percentage for Be (13.8%) followed by Sr (10.9%), Se (10.2%), Ba (10.0%) and B (2.9%). The Li and V showed the lowest percentage in the fraction (less than 0.3%). However, the potentially available fractions represented by the II and III fractions represented a higher percentage in the case of Ba (22.5% in reducible and 9.6% in the oxidizable fraction) and Li (25.5% in the oxidizable fraction).

The soil parameters were selectively correlated with the EEC fractions, and we most often observed correlations with pH/KCl, followed by SOC and manganese hydroxides.

The highest mean value of the contamination factor was calculated for Li (5.8), followed by Ba (1.5) and B (1.4). We also found a possible potential health risk for children in the case of Li (HQ ranged from 0.128 to 1.478). The results underscore the significance of recognizing the interconnections between environmental and human health. Such integrating framework of environmental risk assessment posed by EECs comprising chemical and ecotoxicological lines of evidence can improve decision-making processes and help to manage contaminated sites.

The datasets used for the current study are available from the corresponding author on request.

European Commission (EC). Critical raw materials resilience: Charting a path towards greater security and sustainability. COM(202)474 final. https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:52020DC0474&from=EN (2020a).

European Commission (EC). Study on the EU’s list of Critical Raw Materials. Final report. https://doi.org/10.2873/11619 (2020b).

Bolan, N. et al. From mine to mind and mobiles—Lithium contamination and its risk management. Environ. Pollut. 290, 118067. https://doi.org/10.1016/j.envpol.2021.118067 (2021).

Article CAS PubMed Google Scholar

Islam, M. R. et al. Beryllium in contaminated soils: Implication of beryllium bioaccessibility by different exposure pathways. J. Hazard. Mater. 421, 126757. https://doi.org/10.1016/j.jhazmat.2021.126757 (2022).

Article CAS PubMed Google Scholar

Ramírez-Malule, H., Quinones-Murillo, D. & Manotas-Duque, D. Emerging contaminants as global environmental hazards. A bibliometric analysis. Emerg. Contam. 6, 179–193. https://doi.org/10.1016/j.emcon.2020.05.001 (2020).

Article Google Scholar

Richardson, S. D. & Kimura, S. Z. Emerging environmental contaminants: Challenges facing our next generation and potential engineering solutions. Environ. Technol. Innov. 8, 40–56. https://doi.org/10.1016/j.eti.2017.04.002 (2017).

Article Google Scholar

Franzaring, J., Schlosser, S., Damsohn, W. & Fangmeier, A. Regional differences in plant levels and investigations on the phytotoxicity of lithium. Environ. Pollut. 216, 585–865. https://doi.org/10.1016/j.envpol.2016.06.059 (2016).

Article CAS Google Scholar

Liščák, P. & Antalík, M. Information system of important geosites in the Slovak republic. Slovak. Geol. Mag. 18(1), 5–68 (2018).

Google Scholar

USEPA. Method 3051A: Microwave Assisted Acid Digestion of Sediments, Sludges, Soils, and Oils. (Camgridge University Press, 2007). https://www.epa.gov/sites/default/files/2015-12/documents/3051a.pdf.

Rauret, G. et al. Improvement of the BCR three step sequential extraction procedure prior to certification of new sediment and soil reference materials. J. Environ. Monit. 1, 57–61. https://doi.org/10.1039/A807854H (1999).

Article CAS PubMed Google Scholar

Varol, M. Assessment of heavy metal contamination in sediments of the Tigris river (Turkey) using pollution indices and multivariate statistical techniques. J. Hazard. Mater. 195, 355–364. https://doi.org/10.1016/j.jhazmat.2011.08.051 (2011).

Article CAS PubMed Google Scholar

Tomlinson, D., Wilson, J., Harris, C. & Jeffrey, D. Problems in the assessment of heavy-metal levels in estuaries and the formation of a pollution index. Helgoländer Meeresuntersuchungen 33(1–4), 566–575. https://doi.org/10.1007/BF02414780 (1980).

Article Google Scholar

Rinklebe, J., Antoniadis, V., Shaheen, S. M. & Rosche, O. Health risk assessment of potentially toxic elements in soils along the Central Elbe River Germany. Environ. Int. 126, 76–88. https://doi.org/10.1016/j.envint.2019.02.011 (2019).

Article CAS PubMed Google Scholar

USEPA. Supplemental Guidance for Developing Soil Screening Levels for Superfund Stes. Document OSWER 9355. 4–24, December 2002. Office of Emergency and remedial Response, Washington, DC. https://nepis.epa.gov/Exe/ZyPDF.cgi/91003IJK.PDF?Dockey=91003IJK.PDF (2002a).

USEPA. Strontium; CASRN 7440-24-6. Integrated Risk Information system (IRIS). U.S. Environmental Protection Agency Washington, DC. https://iris.epa.gov/static/pdfs/0550_summary.pdf (2002b).

USEPA. Boron and Compounds; CASRN 7440-42-8. Integrated Risk Information system (IRIS). U.S. Environmental Protection Agency Washington, DC. https://cfpub.epa.gov/ncea/iris/iris_documents/documents/subst/0410_summary.pdf (2004).

USEPA. Provisional Peer reviewed Toxicity values for Lithium. CASRN 7439-93-2. December 2008. U.S. Environmental Protection Agency Washington, DC. https://cfpub.epa.gov/ncea/pprtv/documents/Lithium.pdf (2008).

Reimann, C. et al. GEMAS: Establishing geochemical background and threshold for 53 chemical elements in Eruopean agricultural soil. Appl. Geochem. 88, 302–318. https://doi.org/10.1016/j.apgeochem.2017.01.021 (2018).

Article ADS CAS Google Scholar

Négrel, P. et al. GEMAS: Geochemical background and mineral potential of emerging tech-critical elements in Europe revealed from low-sampling density geochemical mapping. Appl. Geochem. 111, 104425. https://doi.org/10.1016/j.apgeochem.2019.104425 (2019).

Article ADS CAS Google Scholar

Zhang, Y. et al. Extraction of lithium and aluminium from bauxite mine tailings by mixed acid treatment without roasting. J. Hazard. Mater. 404, 124044. https://doi.org/10.1016/j.jhazmat.2020.124044 (2021).

Article CAS PubMed Google Scholar

Qin, W., Han, D., Song, X. & Liu, S. Sources and migration of heavy metals in a karst water system under the threats of an abandoned Pb-Zn mine, Southwest China. Environ. Pollut. 277, 116774. https://doi.org/10.1016/j.envpol.2021.116774 (2021).

Article CAS PubMed Google Scholar

Kanianska, R., Varga, J., Benková, N., Kizeková, M. & Jančová, Ľ. Floodplain soils contamination assessment using the sequential extraction method of heavy metals from past mining activities. Sci. Rep. 12, 2927. https://doi.org/10.1038/s41598-022-06929-7 (2022).

Article ADS CAS PubMed PubMed Central Google Scholar

Davidson, Ch. M. et al. The long-term environmental behaviour of strontium and barium released from former mine workings in the granites of the Sunart region of Scotland, UK. Chemosphere 58, 793–798. https://doi.org/10.1016/j.chemosphere.2004.08.101 (2005).

Article ADS CAS PubMed Google Scholar

Wang, Y., Dong, R., Zhou, Y. & Luo, X. Characteristics of groundwater discharge to river and related heavy metal transportation in a mountain mining area of Dabaoshan, Southern China. Sci. Total Environ. 679, 346–358. https://doi.org/10.1016/j.scitotenv.2019.04.273 (2019).

Article ADS CAS PubMed Google Scholar

Jones, A., Rogerson, M., Greenway, G., Potter, H. A. B. & Mayes, W. M. Mine water geochemistry and metal flux in a major historic Pb-Zn-F orefield, the Yorkshire Pennines, UK. Environ. Sci. Pollut. Res. 20, 7570–7581. https://doi.org/10.1007/s11356-013-1513-4 (2013).

Article CAS Google Scholar

Fernandez-Ondono, E. et al. Use of BCR sequential extraction procedures for soils and plant metal transfer predictions in contaminated mine tailings in Sardinia. J. Geochem. Explor. 172, 133–141. https://doi.org/10.1016/j.gexplo.2016.09.013 (2017).

Article CAS Google Scholar

Caraballo, A. M. et al. Uncertainty in the measurement of toxic metals mobility in mining/mineral wastes by standardized BCR SEP. J. Hazard. Mater. 360, 587–593. https://doi.org/10.1016/j.jhazmat.2018.08.046 (2018).

Article CAS PubMed Google Scholar

Chen, L., Liu, J., Hu, W., Gao, J. & Yang, J. Vanadium in soil-plant system: Source, fate, toxicity, and bioremediation. J. Hazard. Mater. 405, 124200. https://doi.org/10.1016/j.jhazmat.2020.124200 (2021).

Article CAS PubMed Google Scholar

Mikavica, I. et al. Concentration and mobility of trace elements (Li, Ba, Sr, Ag, Hg, B) and macronutrients (Ca, Mg, K) in soil-orchid system on different bedrock types. Environ. Sci. Pollut. Res. 30, 979–995. https://doi.org/10.1007/s11356-022-22110-z (2022).

Article CAS Google Scholar

Xiong, Y., Qi, H., Hu, R., Xiao, Y. & Wei, L. Lithium isotope behavior under extreme tropical weathering: A case study of basalts from the Hainan Island, South China. Appl. Geochem. 140, 105295. https://doi.org/10.1016/j.apgeochem.2022.105295 (2022).

Article ADS CAS Google Scholar

Jacobsson, E. et al. Towards a unified understanding of lithium action in basic biology and its significance for applied biology. J. Membr. Biol. 250, 587–604. https://doi.org/10.1007/s00232-017-9998-2 (2017).

Article CAS Google Scholar

Kot, F. S. On the rubidium and lithium content and availability in the sub-arid south-eastern Mediterranean: Potential health implications. Environ. Geochem. Health 40, 1841–1851. https://doi.org/10.1007/s10653-018-0134-8 (2018).

Article CAS PubMed Google Scholar

Aral, H. & Vecchio-Sadus, A. Toxicity of lithium to humans and the environment—A literature review. Ecotoxicol. Environ. Saf. 70, 349–356. https://doi.org/10.1016/j.ecoenv.2008.02.026 (2008).

Article CAS PubMed Google Scholar

Robinson, B. H., Yalamanchjali, R. & Reiser, R., Dicinson, N.M.,. Lithium as an emerging environmental contaminant: Mobility in the soil-plant system. Chemosphere 197, 1–6. https://doi.org/10.1016/j.chemosphere.2018.01.012 (2018).

Article ADS CAS PubMed Google Scholar

Négrel, P., Reimann, C., Ladenberger, A., Birke, M. Distribution of lithium in agricultural and grazing land soils at European continental scale (GEMAS project). Geophys. Res. Abst. 19, EGU2017-15340. https://meetingorganizer.copernicus.org/EGU2017/EGU2017-15340.pdf (2017).

Li, W., Liu, X. M. & Chadwick, O. A. Lithium isotope behaviour in Hawaiian regoliths: Soil–atmosphere–biosphere exchanges. Geochim. Cosmochim. Acta 285, 175–192. https://doi.org/10.1016/j.gca.2020.07.012 (2020).

Article ADS CAS Google Scholar

Motowicka-Terelak, T., Terelak, H., Pietruch, C. Content of lithium in soils of some regions of Poland. Zeszyty Problemowe Postepow Nauk Rolniczych 448b, 217–224. (1997) (in Polish).

Kashin, V. K. Lithium in sols and plants of western Transbaikalia. Eurasian Soil Sci. 52, 359–369. https://doi.org/10.1134/S1064229319040094 (2019).

Article ADS CAS Google Scholar

Kabata-Pendias, A., Mukherjee, A.B. Trace elements of group 1. In Trace Elements from Soil to Human. 87–104 (Springer, 2007). https://doi.org/10.1007/978-3-540-32714-1_8 (2007).

Peana, M. et al. Environmental barium: Potential exposure and health-hazards. Arch. Toxicol. 95, 2605–2612. https://doi.org/10.1007/s00204-021-03049-5 (2021).

Moffett, D. et al. Toxicological profile for barium and barium compounds. Agency for Toxic Substances and Disease Registry, US. https://www.atsdr.cdc.gov/toxprofiles/tp24.pdf (2007).

Gad, S.C. Barium. In Encyclopedia of Toxicology. 368–370. https://doi.org/10.1016/B978-0-12-386454-3.00819-8 (2014).

Norini, M. P. et al. Mobility of Pb, Zn, Ba, As and Cd toward soil pore water and plants (willow and ryegrass) from a mine soil amended with biochar. J. Environ. Manag. 232, 117–130. https://doi.org/10.1016/j.jenvman.2018.11.021 (2019).

Article CAS Google Scholar

Cappuyns, V. Barium (Ba) leaching from soils and certificied reference materials. Appl. Geochem. 88, 68–84. https://doi.org/10.1016/j.apgeochem.2017.05.002 (2018).

Article ADS CAS Google Scholar

Ezzeldin, H. A., Sallam, A. A. A., Helal, T. Y. & Fouad, H. A. Effect of some materials on Seasamia cretica infesting some maize and sorghum varieties. Arch. Phytopathol. Plant Prot. 42(3), 277–290. https://doi.org/10.1080/03235400601037180 (2009).

Article CAS Google Scholar

González-Valoys, A. C. et al. Ecological and health risk assessments of an abandoned gold mine (Remance, Panama): Complex scenarios need a combination of indices. Int. J. Environ. Res. Public Health 18, 9369. https://doi.org/10.3390/ijerph18179369 (2021).

Article PubMed PubMed Central Google Scholar

WHO. Environmental health criteria 107: Barium. International programme on chemical safety, Geneva. https://apps.who.int/iris/bitstream/handle/10665/39400/9241571071-eng.pdf?sequence=1 (1990).

Carbonell, A. A., Pulido, R., De Laune, R. D. & Patrick, W. H. Soluble barium in barite and phosphogypsum amended Mississippi river alluvial sediment. J. Environ. Qual. 28, 316–321. https://doi.org/10.2134/jeq1999.00472425002800010039x (1999).

Article CAS Google Scholar

Frohne, T., Diaz-Bone, R. A., Laing, G. D. & Rinklebe, J. Impact of systematic change of redox potential on the leaching of Ba, Cr, Sr, and V from a riverine soil into water. J. Soils Sediments 15, 622–633. https://doi.org/10.1007/s11368-014-1036-8 (2015).

Article CAS Google Scholar

Atun, G. & Bascetin, E. Adsorption of barium on kaolinite, illite and montmorillonite at various ionic strengths. Radiochim. Acta 91(4), 223–228. https://doi.org/10.1524/ract.91.4.223.19964 (2003).

Article CAS Google Scholar

Pereira, G. L. et al. Boron: More than an essential element for land plants?. Front. Plant Sci. 11, 610307. https://doi.org/10.3389/fpls.2020.610307 (2021).

Article PubMed PubMed Central Google Scholar

Uluisik, I., Karakaya, H. C. & Koc, A. The importance of boron in biological systems. J. Trace Elem. Med Biol. 45, 156–162. https://doi.org/10.1016/j.jtemb.2017.10.008 (2018).

Article CAS PubMed Google Scholar

Landi, M., Margaritopoulou, T., Papadakis, I. E. & Aranti, F. Boron toxicity in higher plants: An update. Planta 250, 1011–1032. https://doi.org/10.1007/s00425-019-03220-4 (2019).

Article CAS PubMed Google Scholar

Brdar-Jokanovič, M. Boron toxicity and deficiency in agricultural plants. Int. J. Mol. Sci. 21, 1424. https://doi.org/10.3390/ijms21041424 (2020).

Article CAS PubMed PubMed Central Google Scholar

Alan, M. & Kara, D. Comparison of a new sequential extraction method and the BCR sequential extraction method for mobility assessment of elements around boron mines in Turkey. Talanta 197, 189–198. https://doi.org/10.1016/j.talanta.2018.10.030 (2019).

Article CAS Google Scholar

García-Sánchez, F. et al. Multiple stresses occuring with boron toxicity and deficiency in plants. J. Hazard. Mater. 397, 122713. https://doi.org/10.1016/j.jhazmat.2020.122713 (2020).

Article CAS PubMed Google Scholar

Hou, J., Evans, L. J. & Speirs, G. A. Boron fractionation in soils. Commun. Soil Sci. Plant Anal. 25(9–10), 1841–1853. https://doi.org/10.1080/00103629409369157 (2008).

Article Google Scholar

Download references

This work was supported by the Slovak Scientific Grant Agency—VEGA under contract no. VEGA 1/0184/21. The authors acknowledge the effort of the laboratory engineer Marie Králová assisting with the geochemistry analyses (Czech University of Life Sciences Prague).

Faculty of Natural Sciences, Matej Bel University Banská Bystrica, Tajovského 40, 974 01, Banská Bystrica, Slovakia

Radoslava Kanianska, Marek Drimal & Jozef Varga

Faculty of Environmental Sciences, Czech University of Life Sciences Prague, Kamýcká 129, 165 00, Praha-Suchdol, Czech Republic

Michael Komárek

Faculty of Agrobiology, Food and Natrual Resources, Czech University of Life Sciences Prague, Kamýcká 129, 165 00, Praha-Suchdol, Czech Republic

Samuel Kudjo Ahado

Department of Applied and Landscape Ecology, Faculty of AgriSciences, Mendel University in Brno, Zemědelská 1, 613 00, Brno, Czech Republic

Milada Šťastná

National Agricultural and Food Centre, Research Institute of Plant Production,Grassland and Mountain Agriculture Institute, Mládežnícka 36, 974 21, Banská Bystrica, Slovakia

Miriam Kizeková & Ľubica Jančová

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

R.K. conceptualization, methodology, formal analysis, data analysis, investigation, writing original draft; M.D. methodology, data analysis; J.V. formal analysis, investigation; M.Ko. conceptualization, methodology, validation; S.K.A. methodology, data analysis, M.S. conceptualization, methodology, validation; M.Ki. methodology, investigation; L.J. formal analysis, investigation. All authors reviewed the manuscript.

Correspondence to Radoslava Kanianska.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Kanianska, R., Drimal, M., Varga, J. et al. Critically raw materials as potential emerging environmental contaminants, their distribution patterns, risks and behaviour in floodplain soils contaminated by heavy metals. Sci Rep 13, 9597 (2023). https://doi.org/10.1038/s41598-023-36477-7

Download citation

Received: 08 February 2023

Accepted: 04 June 2023

Published: 13 June 2023

DOI: https://doi.org/10.1038/s41598-023-36477-7

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.