Assessment of the heavy metal contamination in the surficial sediments of Lake Kalimanci (Macedonia): a preliminary study Ocena onesnaženosti recentnega sedimenta iz Kameniškega jezera (Makedonija) s težkimi kovinami - preliminarni rezultati Nastja Rogan Šmuc1' *, Petra Vrhovnik1, Tadej Dolenec1, Todor Serafimovski2, Goran Tasev2 & Matej Dolenec1 1Faculty of Natural Sciences and Engineering, University of Ljubljana, Department of Geology, Aškerčeva cesta 12, SI-1000 Ljubljana, Slovenia 2Faculty of Mining, Geology and Polytechnics,University "Goce Delcev - Stip", Goce Delcev 89, 2000 Stip, Macedonia Corresponding author. E-mail: nastja.rogan@guest.arnes.si Received: October 23, 2009 Accepted: November 23, 2009 Abstract: A major disaster happened in the eastern Macedonia (Europe), when part of the Sasa Mine tailings dam collapsed and caused an intensive flow of tailings material which ended up in Lake Kalimanci. The aim of this study was to assess the heavy metal contamination (Ag, As, Cd, Cu, Mo, Pb, Sb and Zn) of the surficial sediments from Lake Kalimanci. The pollution status of the sediments was evaluated by application of two environmental indices: the enrichment factor and a geoaccumulation index. The enrichment factor values (Cd 735; Pb 386; Zn 151; Ag 146; As 48; Cu 21; Sb 10 and Mo 2) revealed extremely elevated heavy metal concentrations in the investigated samples. The geoaccumulation index indicated that surfi-cial sediments of Lake Kalimanci can be treated as moderately to strongly polluted with Pb and Cd, moderately polluted with Zn and Ag, unpolluted to moderately polluted with Cu and As and unpolluted with Mo and Sb. Izvleček: Na območju vzhodne Makedonije se je v deponijskem predelu rudnika Sasa zgodila huda ekološka nesreča. Del jalovinskega nasipa se je zrušil in tako povzročil silovit tok flotacijsko-jalovinskega materiala, ki je zgrmel v dolino reke Kamenice, vse do jezera Ka- menica. Namen študije je bil oceniti onesnaženje površinskega sedimenta iz jezera Kamenica s težkimi kovinami (antimon, arzen, baker, cink, kadmij, molibden, srebro in svinec). Stopnjo onesnaženja smo določili z uporabo dveh okoljskih indeksov: faktorja obogatitve in geoakumulacijskega indeksa. Izračunane vrednosti faktorja obogatitve v raziskovanih vzorcih so bile zelo visoke (Cd 735; Pb 386; Zn 151; Ag 146; As 48; Cu 21; Sb 10 in Mo 2). Rezultati indeksa geoakumulacije so razkrili podobne vrednosti in označili površinske sedimente iz jezera Kamenica kot srednje do visoko onesnažene z Pb in Cd, srednje onesnažene z Zn in Ag, neonesnažene do srednje onesnažene s Cu in As ter neonesnažene z Mo in Sb. Key words: heavy metals, surficial sediments, enrichment factor, geoa-ccumulation index, Lake Kalimanci Ključne besede: težke kovine, površinski sedimenti, faktor obogatitve, indeks geoakumulacije, jezero Kalimanci Introduction Tailings dams are a special type of dam built to store mill and tailings material from mining activities. Currently, thousands of tailings dams worldwide contain billions of tonnes of waste material from mineral processing activity at mine sites. They are supposed to last forever, but since 1970 there have been 35 major failures reported around the world (Diehl, 2001; Goransson et al, 2001; Grimalt et al, 1999; UNEP/ OCHA, 2000). A major environmental disaster also happened in the eastern Republic of Macedonia. In the afternoon of 30 August 2003, part of the Sasa Mine tailings dam collapsed and caused an intensive flow of tailing materials through the Kamenica River valley. The estimated height of the tailings flow was around ten metres and the length of the flow was approximately 12 km. Some 70 000-100 000 m3 of tailings material was discharged and spread through the Kamenica River valley, down to the city of Kamenica and into Lake Kalimanci. The damaging tailings flow comprised a large amount of heavy metals (Dolenec et al, in press) and seriously affected the surrounding environment, especially Lake Kalimanci. When heavy metals enter an aquatic environment, they are redistributed throughout the water column, deposited or accumulated in sediments and consumed by biota (Long et al, 1996; Fichet et al, 1998; Linnik & Zubenko, 2000; Singh et al, 2005; Ghrefat & Yusuf, 2006 and Khaled et al, 2006). A fundamental characteristic of heavy metals is their lack of biodegradability and lake sediments usually operate as pollutant storage tanks that reflect long-term impacts (Maltby, 1992 and Schmitt et al, 2003). Therefore it is important to assess the extent of heavy metal accumulation in lake sediment and to determine its environmental threat. Our objective in the preliminary study was to define the total heavy metal contamination (Ag, As, Cd, Cu, Mo, Pb, Sb and Zn) of the surficial sediments from Lake Kalimanci and to evaluate enrichment factor (EF) and geoaccumulation index (I ) values for the heavy metals present in the sediment samples. Materials and methods Study area Kalimanci artificial lake is located in eastern Macedonia, in the 2000 m high Osogovo Mountains, near the city Makedonska Kamenica and the Sasa mine (Figure 1). The surface of the lake is 42 km2, the maximum depth is 85 m and it encompasses around 120 • 106 m3 of water. It is supplied by two tributaries: the Bregalnica and the Kamenica River, which flows directly from the Sasa mine and drains a large amount of mining effluents and tailings material which have discharged into Lake Kalimanci. Figure 1. Geographical map of the investigated area (merilo 1 : 625 000) The sediment samples from Lake Ka-limanci are dominated by: quartz, pla-gioclases, K-feldspars, muscovite, il-lite and clinochlorite. Subordinate minerals are hornblende, gypsum, bassanit, calcite, dolomite, smithsonite, pyrite, marcasite, haematite, goethite and dia-spor (Dolenec et al, unpublished). The Sasa-Toranica ore district is situated 10 km N of Lake Kamenica in the Os-ogovo Mountains and occupies an area of about 200 km (Figure 1). It is established as one of the largest ore districts within the Besna Kobila Osogovo Tas-sos metallogenetic zone. The important Pb and Zn ore bodies are usually found in quartz-muscovite-graphitic schists and also in greenschists and marbles. The ore bodies are always accompanied by variable amounts of Cu, Au, Ag, Mo and Sb. The Sasa mine has been in production for over 45 years, yielding 90 000 t of high quality Pb-Zn concentrate annually and numerous tons of tailings material (Serafimovski et al, 2005). The tailings material from the Sasa mine is made up of quartz, pyrite, ga-lenite, gypsum, hornblende, actinolite, albite, anortite, biotite and orthoclase. The geochemical analysis of tailings material indicated an average range of Ag 0.004 pg/g, As 69.2 pg/g, Cd 84.3 pg/g, Cu 279.3 pg/g, Mo 2.9 pg/g, Pb 5595.2 pg/g, Sb 4.2 pg/g and Zn 6970.2 pg/g (Dolenec et al, unpublished). Sampling and analysis Seventeen surficial sediment samples from Lake Kalimanci were taken in September 2007, three years after the tailings dam accident happened. The surficial sediments were composed of silt. The chosen sampling locations were formed into eight profiles (Figure 2), covering the area around the River Kamenica tributary in the northern site of the lake which was mostly affected by the tailings flow. The samples were collected with the plastic corers (tube 20 cm long with a 10 cm internal diameter), tightly packed into plastic bags and stored in the laboratory at 4 °C. The sediment samples were dried at 50 °C for 48 h. They were then sieved through a 0.315 mm polyethylene sieve to remove plant debris and homogenised by a mechanical agate grinder to a fine powder for subsequent analysis. The mineralogy of the sediment samples was determined at the Department of Geology, Ljubljana (Slovenia) by X-ray powder diffractometry with a Philips PW 3710 diffractometer and CuKa radiation. The diffraction patterns were identified with the data from Powder Diffraction File (1977) - JPDS system. The geochemical analysis of the following heavy metals, Ag, As, Cd, Cu, Mo, Pb, Sb and Zn, was obtained in a certified commercial Canadian labo- alysed by ICP Mass Spectrometry. The ratory (Acme Analytical Laboratories, analytical precision and the accuracy Ltd.) and 0.5 g of samples were lea- were better than ±5 % for the analysed ched in hot (95 °C) Aqua Regia and an- elements. Figure 2. Sampling locations in the Lake Kalimanci Statistical analysis Basic statistical analysis of data was performed by use of the original statistical software program Statistica 6. Pearson multiple correlation analysis was also applied to all results. Results and discussion Heavy metals (Ag, As, Cd, Cu, Mo, Pb, Sb and Zn) were examined for their abundance in the surficial sediments from Lake Kalimanci. When the accident occured, the damaging tailings flow mostly affected the northern part of the Kalimanci Lake. Consequently, the maximum concentrations of all heavy metals were determined in the northern part of the basin (Table 1, Figure 2). The maximum level of the Zn (20 900 pg/g), Pb (16 300 pg/g), Cu (1 162 pg/g), Cd (136 pg/g), Ag (17.3 pg/g), Mo (4.6 pg/g), and Sb (3.6 pg/g) was recorded in location II-3 and the highest amount of As was found in location II-1 (Table 1, Figure 2). Table 1. Concentrations c/(pg/g) of the heavy metals in the investigated samples Element c/ (pg/g) Ag As Cd Cu Mo Pb Sb Zn Location I-2 2 54.8 48.1 238 2.2 2721 1.4 5913 I-4 10.3 68.8 74 627.3 2.9 9357 2.4 10700 II-1 3.6 86.7 77.2 415 3.1 4461 2.2 9596 II-3 17.3 77.2 136 1162 4.6 16300 3.6 20900 II-5 15.2 70.1 111.6 928.7 4.2 13800 3.2 17600 III-3 11.1 73.2 89.6 723 3.7 10900 2.7 14000 III-6 10.4 57.6 86.8 692.5 3.3 9800 2.1 14000 IV-1 IV-4 5.1 66.1 59 413.1 3 6695 1.7 8105 4.2 61.6 47.6 341.4 2.3 5343 1.6 6734 V-1 3.6 61.6 40 303.2 2.4 4447 1.3 5677 V-7 8.1 66 81.1 595.7 3.7 9472 2 12600 VI-7 5.6 128.2 54.2 315.4 3.6 7885 3.3 7181 VI-11 7.1 58.6 77.6 545.7 2.5 7557 1.6 11600 VII-1 3.6 77.8 61.3 372.5 1.9 4575 1.1 10400 VII-12 4.5 66.4 53.5 398.1 2.7 5144 1.2 9326 VIII-1 3.8 69.8 47.9 311.9 2.2 5091 1.4 7224 VIII-8 4.1 61.5 46.8 327.6 2.3 4893 1.4 7056 Basic statistical information are given in the Table 2. Correlation analysis A Pearson correlation matrix was used to assess element associations and metal origins (Table 3). The correlation matrix presents a significantly positive correlation between Ag, Cd, Cu, Mo, Pb, Sb and Zn, demonstrating possible co-contamination from similar sources: in our case, that would be tailings material from the Sasa mine. Conversely, no correlations were noted between As and the other elements (except Sb), suggesting that As contamination might be from a different source or it had a natural origin but was not from the same pollution source or had a different depositional nature. Element Valid N Mean Minimum Maximum Std. Dev. Ag 17 7.04 2.000 17.30 4.414 As 17 70.94 54.800 128.20 16.859 Cd 17 70.14 40.000 136.00 25.773 Cu 17 512.42 238.000 1162.00 250.845 Mo 17 2.98 1.900 4.60 0.775 Pb 17 7555.35 2721.000 16300.00 3664.172 Sb 17 2.01 1.100 3.60 0.783 Zn 17 10506.59 5677.000 20900.00 4248.268 Table 3. Pearson correlation matrix showing inter-elemental relationship (n = 17) Ag As Cd Cu Mo Pb Sb Zn Ag 1.00 As 0.03 1.00 Cd 0.93* 0.06 1.00 Cu 0.98* -0.03 0.97* 1.00 Mo 0.85* 0.34 0.84* 0.83* 1.00 Pb 0.99* 0.14 0.91* 0.96* 0.90* 1.00 Sb 0.80* 0.55* 0.76* 0.74* 0.91* 0.84* 1.00 Zn 0.94* -0.00 0.98* 0.98* 0.80* 0.92* 0.69* 1.00 *p < 0.05 Table 2. Descriptive basic statistics of the heavy metals in the surficial sediments of Lake Kalimanci Enrichment factor The enrichment factor (EF) is a useful indicator reflecting the condition of environmental contamination. It was employed to evaluate possible anthropogenic input of metals to observed sediment (Brady, 1984), calculated according to the equation: EF = (M/Al) ,/ (M/Al) v y sample v A (1) where M , and M „ are the levels sample crust of the investigated metals (Ag, As, Cd, Cu, Mo, Pb, Sb and Zn) in the sediment samples and in the uncontaminated crust material, respectively, and Al sample and Al are the values of the Al in crust the sediment samples and in the un-contaminated crust material. Metal to aluminium ratios are widely adopted, presumably because the concentration of Al in weathering products and their parent materials are generally comparable. Al is also the normalising element assumed not to be enriched owing to local contamination. Baseline values for M were adopted from Taylor and crust McLennan (1985). EF-values lower than and around 1.0 indicate that the element in the sediment originates predominantly from the crustal material and/or weathering processes (Zhang & Liu, 2002), whereas EF-values much greater than 1.0 display the anthropogenic origin of the element (Szefer et al, 1996). According to Chen et al (2007), EF < 3 indicates minor enrichment (anthropogenic impact), EF = 3-5 moderate enrichment, EF = 5-10 moderately severe enrichment, EF = 10-25 severe enrichment, EF = 25-50 very severe enrichment, and EF > 50 extremely severe enrichment. The calculation of enrichment factors showed that all studied heavy metals were enriched in surficial sediments of Lake Kalimanci. Cd had the highest average EF-value (735) among the investigated heavy metals, which represents extremely severe enrichment. EF-values of Pb, Zn and Ag (average value 387, 151 and 146, respectively) also signified extremely severe enrichment. Calculated EF-values for As, Cu and Sb determined very severe enrichment (average value 48) with As and severe enrichment with Cu and Sb (average values 21 and 10, respectively). Mo exhibited the lowest EF-values among the heavy metals studied (average value 2) and displayed moderate enrichment (Table 4). Geoaccumulation index The geoaccumulation index (/ ) was also used to assess heavy metal pollution in surficial sediments of Lake Kalimanci. It is expressed by the following equation (Müller, 1969): I = log? (C /1.5 B) geo c? v n n' (2) geo ®2 ^ n where C is the measured concentration Table 4. Enrichment factors (EF) values (average) in surface sediments of Lake Kali-manci Element Ag As Cd Cu Mo Pb Sb Zn EF values (average) 146 48 735 21 2 387 10 151 Table 5. Igeo values (average) in surface sediments of Lake Kalimanci Element Ag As Cd Cu Mo Pb Sb Zn I values geo (average) 1.8 0.6 2.2 0.9 -0.1 2.4 -0.05 1.9 of the heavy metal (n) in the sediments, Bn is the geochemical background value in the average shale of element n, and 1.5 is the background matrix correction factor owing to lithogenic effects (Loska et al, 1997; Ghrefat & Yusuf, 2006; GoNzALES-MAciAS et al, 2006; Chen et al, 2007). According to Loska et al (1997) and GoNzALES-MAciAS et al (2006), the ge- oaccumulation index (I ) can be catv geo' egorised in a scale ranging from 1 to 6: I < 0 unpolluted, I < 1 unpolluted geo r 5 geo r to moderately polluted, Igeo < 2 moderately polluted, Igeo < 3 moderately to strongly polluted, Igeo < 4 strongly polluted, Igeo < 5 strongly to very strongly polluted, and Igeo > 5 very strongly polluted. (1.9) and Ag (1.8), unpolluted to moderately polluted with Cu (0.9) and As (0.6) and unpolluted with Mo (-0.1) and Sb (-0.05) (Table 5). The heavy metals in sediments are derived from two sources: natural and anthropogenic. In Lake Kalimanci, however, extremely elevated concentrations of Cd, Pb, Zn, Ag, As, Cu, Sb, Cu, Sb, and Mo in the surficial sediments clearly represent anthropogenic impact originating from the mining activities (mining effluents and tailings material) and the tailings dam accident in the Sasa mine. This anthropogenic influence was also evaluated and confirmed by the calculation of enrichment factors and the index of geoaccumulation. The results (average values of the Igeo) revealed that surficial sediments of Lake Kalimanci were moderately to strongly polluted with Pb (2.4) and Cd (2.2), moderately polluted with Zn Conclusion In the present study we applied two ecological indices, enrichment factor and index of geoaccumulation, to as- sess and evaluate the distribution of heavy metals (Ag, As, Cd, Cu, Mo, Pb, Sb and Zn) in surficial sediments from Lake Kalimanci (Macedonia). The results indicate that all heavy metals were enriched in the sediment samples studied. According to the calculated values of I the surficial sediments were geo strongly polluted with Pb and Cd and unpolluted with Mo and Sb. This pollution impact clearly originates from the mining activities (acid mine drainage) of the Sasa mine and from the Sasa tailings dam accident. Acknowledgments The research was financially supported by the Slovenian Research Agency (ARRS), contract number 1000-05310229. References Chen, C. W., Kao, C. M., Chen, C. F., & Dong, C. D. (2007): Distribution and accumulation of heavy metals in the sediments of Kaohsiung Harbor, Taiwan. Chemosphere; 66(8), 1431-1440. Diehl, P. (2001): World Information Service on Energy (WISE) Uranium project [online]. Available from: http://www.antenna.nl/wise/urani-um/. Fichet, D., Radenac, G., & Miramand, P. (1998): Experimental studies of impacts of harbour sediments resuspension to marine invertebrates larvae: bioavailability of Cd, Cu, Pb and Zn and toxicity. Marine Pollution Bulletin; 36, 509-518. Ghrefat, H., & Yusuf, N. (2006): Assessing Mn, Fe, Cu, Zn, and Cd pollution in bottom sediments of Wadi Al-Arab Dam, Jordan. Chemosphere; 65, 2114-2121. Günzáles-Macías, C., Schifter, I., Llu-CH-CoTA, D. B., MÉNDEZ-RüDRÍGU-ez, L., & Hernández-Vázquez, S. (2006): Distribution, enrichment and accumulation of heavy metals in coastal sediments of Salina Cruz Bay, Mexico. Environmental Monitoring and Assessment; 118, 211-230. Goransson, T., Benckert, A., Lindvall, M., & Ritzén, R. (2001): Dam failure at the Aitik mine: investigations, conclusions and measures taken, in: Proceedings of Securing the Future: International Conference on Mining and the Environment, Sweden. Grimalt, J. O., Ferrer, M., & Macpher-SON, E. (1999): The mine tailing accident in Aznalcollar. The science of the Total Environment; 242, 3-11. Khaled, A., El Nemr, A., & El Sikaily, A. (2006): An assessment of heavy-metal contamination in surface sediments of the Suez Gulf using geoaccumulation indexes and statistical analysis. Chemistry and Ecology; 22(3), 239-252. Linnik, p. M., & Zubenko, I. B. (2000): Role of bottom sediments in the secondary pollution of aquatic environments by heavy-metal compounds. Lakes Reservoirs Research Management; 5, 11-21. Long, E. R., Robertson, A., Wolfe, D. A., Hameedi, J., & Sloane, G. M. (1996): Estimates of the spatial extent of sediment toxicity in major US estuaries. Environmental Science & Technology; 30, 35853592. Loska, K., Cebula, J., Pelczar, J., Wi- ECHULA, D., & Kwapulinski, J. (1997): Use of enrichment, and contamination factors together with geoaccumulation indexes to evaluate the content of Cd, Cu, and Ni in the Bybnik water reservoir in Poland. Water, Air and Soil Pollution; 93, 347-365. Maltby, E. (1992): Soil and wetland functions. In Gerakis, P. A. (ed.), Conservation and Management of Greek Wetlands. The IUCN Wetlands Programme: 9-60. Müller, G. (1969): Index of geoaccumulation in sediments of the Rhine River. Geojournal; 2, 108-118. Schmitt, D., Saravia, F., Frimmel, F. M., & Schuessler, W. (2003). NOM-facilitated transport of metal ions in aquifers: importance of complex-dissociation kinetics and colloid formation. Water Resources; 37:3541-3550. Serafimovski, T., & Aleksandrov, M. (1995): Lead and zinc deposits and occurences in the Republic of Macedonia. Macedonia, Skopje: Special edition of RGF. Singh, K. P., Mohan, D., Singh, V. K. & Malik, A. (2005): Studies on distribution and fractionation of heavy metals in Gomti river sediments - a tributary of the Ganges. Journal of Hydrology; 312, 14-27. Szefer, P., Glasby, G.P., Szefer, K., Pemp-kowiak, J., & Kaliszan, R. (1996): Heavy-metal pollution in superficial sediments from the southern Baltic Sea of Poland. Journal of Environmental Science and Health; 31A, 2723-2754. Taylor, S. R., & McLennan, S. M. (1995): The geochemical evaluation of the continental crust. Reviews of Geophysics; 33,241-265. UNEP/OCHA (2000): Cyanide spill at Baia Mare, Romania, UNEP/ OCHA Assessment Mission, UNEP/Office for the Coordination of Humanitarian Affairs, OCHA. Zhang, J., & Liu, c. L. (2002): Riverine composition and estuarine geochemistry of particulate metals in China—weathering features, anthropogenic impact and chemical fluxes. Estuarine and Coastal Shelf Science; 54, 1051-1070.