2. Materials and Methods
The study area is located in the west of the Kola Peninsula in the Monchegorsk area of the Murmansk region. Five lakes located in the vicinity of buildings and structures that belong to MS KGMK and Monche Bay of Lake Imandra were studied (
Figure 1). In addition, one sample of groundwater was obtained from a self-flowing well near of Kumuzhye Lake (number 7 in
Figure 1).
All the industrial and transport infrastructures in the area have been intensively developed since 1938, and the great amounts of sulfur dioxide and HMs emitted into the atmosphere have resulted in degradation of the soil and vegetation cover (from lichen depression to the complete depletion of soils to technogenic barren formation in the lakes’ landscape). The transformation of the Lake Kumuzhie lowlands as a result of anthropogenic impact was insignificant, and mainly the north-taiga vegetation—typical for the area—was preserved there. The investigated lakes represent a single system and have a water catchment area of 90.54 km
2. Connected to its northern part through pipes, the southern part of the Lake Nyudjavr basin, which is separated from the main water basin by a dam, serves as a technological basin to settle saline wastewater transported from the smelter and domestic waste transported from MS KGMK (
Figure 1). The wastewaters belong to the chloride-sulfate class, and to the sodium group. In 2017, suspended solids in the sewages amounted up to 503 mg/dm
3 and contained extremely high concentrations of sulfate ions (up to 2720 mg/dm
3), sodium ions (up to 2085 mg/dm
3), chloride ions (up to 1246 mg/dm
3) and borate ions (up to 63 mg/dm
3). Concentrations of acid-soluble forms of HMs in wastewaters are also high: Cu up to 20.5 mg/g, Ni up to 156 mg/g, and Co up to 0.96 mg/g.
A channel connects Lake Pyslyschimjavr to the northern part of Lake Nyudjavr, through which waste is transported around the southern part of Lake Nyudjavr. The Nyuduai river basin discharges into Monche Bay of Lake Imandra. Municipal sewage from the Monchegorsk town pours into the northern part of Monche Bay.
We carried out water sampling in 2016 and 2017. There were two water samples taken from the surface and the bottom water layers from each lake and Monche Bay annually. Water samples were taken in polyethylene bottles from Nalgen®, the material of which does not have sorbing properties. Previously, the bottles were thoroughly washed in the laboratory. When sampling water, the bottles were rinsed twice with lake water, then placed in dark containers and refrigerated (~+4 °C) in a short time transported to the laboratory. Filtration of water samples was carried out in the laboratory during discharge at the Millipore phase separation unit from high-density polypropylene through glass and polycarbonate membrane filters of Millipore HVLPO 4700, Schleicher and Schuell ME 25/21 ST, Whatman GF/A with a pore size of 0.45 μm. The following parameters were determined in unfiltered samples: pH, conductivity, alkalinity, NH
4+, NO
3−, PO
43+, metals (after acidification with concentrated nitric acid). In the filtered samples, Si, P, color, Cl
−, SO
42+ were determined. Water chemistry was determined at the resource sharing center of the Institute of North Industrial Ecology Problems (INEP) Kola SC RAS using the same methods described in
Guidelines for Methods… [
13] and
Standard Method… [
14]. For the quality control of the measurements of pH, alkalinity, concentrations of chlorides, sulfates, alkali, and alkaline-earth metals, the specialized ALPEFORM software suite was used, including an assessment of ion balance, measured and calculated electrical conductivity. The quality of the laboratory analytical measurements was evidenced by annual international verification [
15].
A hydro-chemical analysis was carried out by Analytical Laboratory INIEP KSC following the techniques described in
Guidelines for Methods… [
13] and
Standard Method… [
14]. We compared the results of the analysis to the averaged hydro-chemical indices of the lakes in the north-taiga zone of the Murmansk region, taking these indices conditionally as the background ones [
16,
17,
18,
19]. We analysed the chemical composition of groundwater from a self-discharged well located near Lake Kumuzhie (
Figure 1). We sampled and analysed phytoplankton in the summer periods of 2016 and 2017 following the Russian state’s standards [
20] and the standard techniques recommended in
Guidance on Methods… [
21] and
Guidance on Hydrobiological Monitoring… [
22], according to the scheme defined in the INEP [
23]. We sampled phytoplankton with a typical Ruttner water sampler with a 2.2-litre capacity. Plankton was concentrated with a plankton net (29 μm) and the final sample was placed into a 50 ml plastic flask and fixed with a Lugol solution layer. The phytoplankton samples were concentrated in the laboratory through sedimentation (settling). We calculated the phytoplankton biomass using the counting-volume technique, using the counting-volume technique by calculating the volume of an individual cell (or dense colonies) of each species [
24,
25,
26]. We calculated the abundance and determined the taxonomic identification of the algae in a 0.1 mL Nageotte chamber through a light microscope, Motic BA300, with an immersion lens. We merged the taxa names with the international algae database [
27]. The algae communities’ floristic similarity was assessed using the Sörensen–Czekanowski coefficient [
28,
29] and the analysis was carried out by the program module Graphs [
30]. We assessed the species diversity using the Shannon–Weaver index and the Simpson’s reverse (polydominance) index [
31,
32]. To assess water quality, we calculated the saprobic index (S) using the Pantle–Buck method modified by Sladechek [
33,
34], in accordance with the Russian state’s standards [
20]. The ecological characteristics of the discovered taxa were taken from Barinova et al. [
35]. To determine the concentration of chlorophyll “a” in the plankton, we put the water samples (600 mL) through a membrane filter with a pore diameter of 0.47 μm with the help of a Millipore injector and a filtering nozzle. To avoid photosynthetic pigments changing in the water samples in transport, we filtered the water near the water body. We extracted chlorophyll with pure acetone solution (90%, PA) and measured the optical density of the extracts using the Hitachi UV-VIS 181 spectrophotometer. We calculated the photosynthetic pigment concentrations using standard methods, generally accepted in international and domestic practice, following the scheme employed at INIEP KSC RAS [
23,
36,
37].
We sampled and analysed the zoobenthos according to the Russian state’s standards, using standard techniques [
17,
19]. The bottom sediment samples were taken from deep zones with an Ekman–Birge grab, the catchment area being 290 cm
2. For the littoral zone, we took quality samples with a hydrobiological netscraper. Quality and quantity zoobenthos sampling was carried out in the late summer in lakes Kumuzhie, Sopchjavr, Pyslyschimjavr, Travyanoe, Nyudjavr (August 2016 and 2017) and Monche Bay of Lake Imandra (August 2013). Invertebrates were identified by the key to zooplankton and zoobenthos of the fresh waters of European Russia [
38] and the key to freshwater invertebrates in Russia and adjacent territories [
39].
We sampled and analysed the ichthyology material according to Russian state standards applied by the research institutions during ichthyologic researches [
40,
41,
42,
43,
44]. The trout organs were analysed for HM accumulation, and the samples were frozen to be further analysed. To achieve the constatnt weight, the organ samples were dried in a laboratory drying oven at 90 °C. Then, the organic matrix was removed in a concentrated nitric acid solution in a microwave decomposition system (Multiwave 3000, Anton Paar, Austria), followed by further filtration. We determined HM concentration in the fish using atomic adsorption spectroscopy with atomisation in flame or in graphite pipes using a PerkinElmer 5000 in graphite oven. The HM concentration was calculated in µg/g of dry tissue weight. We compared the solution samples used to identify the metals with the certified standard sample of Fluka Chemie GmbH, Switzerland, and controlled quality using DORM-2 standard samples.
3. Results and Discussion
3.1. Hydro-Chemical Indices
All five lakes exceeded the regional background level for all the basic pollutants—HMs Cu, Ni, Pb, and Cd—and for nutrients—nitrogen and phosphorus compounds (
Table 1). The highest pollutant concentrations were recorded in Lake Nyudjavr, which received wastewater in which the chloride, sulfate, copper, and nickel concentrations were two orders of magnitude higher than that of the conditional background [
7]. The hydro-chemical parameters of the other water bodies were mainly the result of aero-technogenic pollution. It was suggested [
9,
16] that the lake ecosystems develop under the influence of two main factors: biogenic and toxic loads. It seemed that the presence of available biogenic elements could increase water organisms’ resistance to the impact of toxicants. Groundwater was characterized by the hydro-chemical parameters that were most like the background ones. The nickel and sulfate concentrations exceeded the background ones (
Table 1).
The inflow of comparatively less polluted groundwater to Lake Kumuzhie was likely to have a positive effect on the water body’s ecosystem. The water chemical composition in Monche Bay of Lake Imandra was affected by the polluted discharge from the Nyuduai River and by sewage works located near the town of Monchegorsk. The nutrients and HM concentrations in Monche Bay were found to be lower than those in the lakes being studied, which can be explained by the high water exchange, the connection with the main water basin of Lake Imandra, and the greater water volumes than those of the five lakes in question.
3.2. Phytoplankton
In total, there were 89 algae and Cyanobacteria taxa found in all the water bodies studied. They were below the rank of the following systematic groups: Cyanobacteria (8), Dinophyceae (5), Cryptophyta (1), Bacillariophyta (43), Chlorophyta (15), Charophyta (12), Chrysophyceae (3), and Euglenophyta (1). All the lakes differed greatly in taxonomic composition, dominating taxa, and species richness (
Figure 2).
In Lake Kumuzhie, the typical representatives were Cyanobacteria Pseudanabaena limnetica (Lemm.) Kom.; dinoflagellates Peridinium bipes Stein; diatoms Aulacoseira alpigena (Grun.) Kramm., Tabellaria flocculosa; green algae Chlorella mucosa Korsh.; Charophyta Spondylosium secedens (De Bary) Arch.; and Spirogyra sp. Green filamentous algae Ulothrix sp. were found in the phytoplankton of Lake Sopchjavr, representatives of Chlamydomonas and dinoflagellates Peridinium bipes. Dominating in Lake Pyslyschimjavr were diniflagellate Peridinium bipes. Diatoms Aulacoseira alpigena (Grun.) Kramm. were also found.
Highly dominating in Lake Travyanoe was
Pseudanabaena limnetica (Lemm.) Kom. Great taxonomic richness was characteristic of Lake Nyudjavr, where diatoms dominated. They included
Fragilaria capucina subsp.
rumpens (Kütz.) Lange-Bert.,
F. capucina var. gracilis (Oestr.) Huste.,
Surirellabre bissonii Kramm., and Lange-Bert., green algae
Stigeoclonium longipilum Kütz., and Cyanobacteria
Anabaena sp. The composition of phytoplankton from Monche Bay of Lake Imandra corresponded to that typical of the eutrophicated zones of Lake Imandra [
45]—
Asterionella formosa Hass.,
A. islandica (Müll.) Simons.,
Fragilaria crotonensis Kitt.—and differed greatly in taxonomic richness.
According to floristic analysis, there was slight similarity in the quality characteristics (the taxonomic composition) and community structure (taking each species’ richness into account) of the phytoplankton flora sampled from the lakes studied. The phytoplankton communities of Lake Nyudjavr and Monche Bay, where diatom algae dominated, as well as those of Lakes Sopchjavr and Pyslyschimjavr, in which dinoflagellates dominated, were found to be the closest in taxonomic composition. The most significant difference in the taxonomic composition was in the algal flora of Lake Travyanoe, where Cyanobacteria dominated (
Figure 2d,e).
As a whole, the phytoplankton species’ composition and the community structure in all the water bodies differed from those in the lakes located in the background areas of Murmansk region: the Cyanobacteria and dinoflagellate levels, indicators of toxic load, were relatively high. The phytoplankton in small lakes included filamentous algae characteristic of fouling. Some features of natural water bodies were provided by algal flora in Monche Bay of Lake Imandra: the diatom level was high, and one could find the Chrysophyceae Dinobryon divergens Imh. and Synura sp. due to the large size of the water body and intensive water-exchange processes.
The state of the phytoplankton algae communities in Monche Bay and of the sensitive-to-pollution Charophyta algae of Lake Kumuzhie was found to be the closest to the natural state. The algal flora in Lake Travyanoe was found to be the most transformed, with utter domination of Cyanobacteria in both amount and biomass. In the slightly saline, alkaline Lake Nyudjavr, under both toxic load and anthropogenic eutrophication, we observed the formation of specific alkaliphilic algal flora in the relatively stable-to-pollution diatom species Surirella brebissonii Kramm. and Lange-Bert, Fragilaria capucina var. gracilis (Oestr.) Hust., and F. capucina subsp. rumpens (Kütz.) Lange-Bert. The presence of toxic substances (primarily HMs) suppressing photosynthetic activity seemed to prevent abundant phytoplankton mass development. The presence of nutrients was likely to increase algae resistivity to toxic load.
The average amount of phytoplankton biomass was small, ranging from 0.13 to 3.43 g/m
3; its abundance, however, varied greatly—from 88 to 56979 ind./L (
Figure 2). The greatest amount of biomass was recorded in Monche Bay. In other water bodies, the quantitative indices of phytoplankton corresponded to those of oligotrophic water. The highest amount of phytoplankton was recorded in Lake Travyanoe, in which one could observe mass development of the Cyanobacteria
Pseudanabaena limnetica. The amount of chlorophyll “a” in the phytoplankton ranged from 0.51 to 6.98 mg/m
3. The highest photosynthesis activity was observed in the algocoenoses in Monche Bay of Lake Imandra, and the lowest activity was observed in Lake Kumuzhie (
Table 2).
We calculated the degree of saprobity (S) using the phytoplankton indices, and all the water bodies were classified from xenosaprobic (x) to beta-oligosaprobic (β-o). We also classified the water quality, which ranged from I (very clear) to III (moderately polluted) within the limnosaprobic category. The lowest S-index (0.36) was characteristic of Lake Sopchjavr, and the highest index S (1.60) was characteristic of Lake Travyanoe. It was clear that the intensive impact produced by non-ferrous plants limited the supplication of the S-index to determine water quality in water bodies for two reasons: firstly, the amount of saprobionts was not enough to be calculated due to the extremely small number of taxa; secondly, the presence of an intensive toxic load suppressed algae development, even with sufficient nutrients (nitrogen and phosphorus compounds). In the present study, the objectively lowered values of the S-index under the biogenic load could be taken as indicators of toxic water pollution.
The presence of enough biogenic elements could be one of the factors affecting plankton organism resistivity to toxic load. Barinova S.S. [
46] showed that nitrate nitrogen is depleted by photoautotrophes in all cases when any impact—for instance, toxic pollution—does not prevent photosynthesis. In most water bodies studied, the NO
3 concentration exceeded the background indices (
Table 1), which indicates that the degree of photosynthesis activity suppression was affected by the toxic impact of HMs. At the same time, the critical concentrations of toxicants hindered the eutrophication processes, supporting the phytoplankton indices’ indication that the trophic status of the lakes was low [
8].
3.3. Zoobenthos
The HM concentration in the water and bottom sediments was the limiting factor for benthos invertebrates. For instance, in the toxic environment of the small lakes of the Pechenga area, the species diversity and amount and biomass of zoobenthos sharply decreased. Of all the systematic groups, only nematods, chironomids, hemiptera, caddisflies, and alderflies (Sialidae) were found in water in which Ni and Cu concentrations were high [
47]. The typical trends were also found in small lakes near the Monchegorsk KGMK site. The macrozoobenthos sampled in the five water bodies was characterized by depleted taxonomic composition and extremely low quantity indices. No bottom organisms were found in Lake Pyslyschimjavr, located in the vicinity of the Monchegorsk KGMK site. The basis of zoobenthos of the remaining water bodies and streams located in this area was formed by chironomids Psectrocladius, Procladius, Cricotopus, and Orthocladius, usually found in water bodies polluted with HMs.
Larvas Psectrocladius dominated the benthos composition of Lake Sopchjavr (95% of the total zoobenthos amount and biomass), and Polypedilum (Pentapedilum) and Procladius dominated in Lake Travyanoe (97%). The amount of bottom invertebrates in the water bodies did not exceed 40 ind./m2. The biomass was equal to 0.2 g/m2. Single chironomids Procladius (Holotanypus) choreus gr, widely found in freshwater bodies of the subarctic zone, were found in the bottom sediments of the deep zone of Lake Nyudjavr.
The benthos fauna diversity was much greater and developed throughout the stony littoral zone of Lake Nyudjavr. Found there were chironomids Psectrocladius, Glyptotendipes, Procladius (Holotanypus) choreus gr., Cricotopus, and Orthocladius. Larvas Psectrocladius dominated in the composition of the communities (66% of the total amount and 54% of the total zoobenthos biomass), and the representatives of Cricotopus sub-dominated (15% and 12%, respectively). Also found there were water bugs: boat bugs (Sigara sp.) and water striders (Gerris sp.). The average amount of the littoral zoobenthos in Lake Nyudjavr was equal to 812 ind./m2. That of the biomass was 1.1 g/m2. The bottom biocenosis in Moncha Bay, Lake Imandra, degraded greatly. In the study periods, no benthos invertebrates were found in the littoral zone of the Nyuduai River’s mouth. In the zoobenthos communities in deep zones, we only found chironomids Chironomus—typical of Moncha Bay.
In the profundal zone of Lake Kumuzhie, the basis of zoobenthos was formed by chironomids Procladius (Holotanypus) choreus gr., Chironomus sp., Polypedilum (Pentapedilum) sp., (>80% of the total macrozoobenthos amount and biomass), and amphipoda Gammarus lacustris Sars, 1863 (18%). In the littoral zone, the dominating complex was complemented with caddis Polycentropus flavomaculatus Pictet, 1834, Molanna sp., and the Diptera larvae (Rhagionidae and Dolihopodidae). In the littoral zone, the quantity indices were low: The zoobenthos amount accounted for 90 ind./m2, and the biomass accounted for 1.8 g/m2; in the profundal zone, these accounted for 154 ind./m2 and 1.2 g/m2, respectively. The macrozoobenthos composition of Lake Kumuzhie was relatively similar to that of the water bodies located in unpolluted areas in the north-taiga zone of the Murmansk region.
Most distinguished zoobenthos systematic groups were rather resistant to HM water pollution and were widely found in both “background” and polluted water bodies and streams in the Murmansk region, excluding amphipod
G. lacustris, which was found only under very low Ni (7 µg/L) and Cu (<5 µg/L) concentrations in water [
47]. This situation was highly likely to occur due to running water in the lake and due to the presence of several permanent wells located in the water basin and on the banks of the water body, whose water, which runs to the lake, helps dilute polluted water and reduce toxic load. The trophic status of the five water bodies was assessed as oligotrophic, which, to a significant extent, was due to the impact exerted on the water bodies by aero-technogenic discharges produced by the smelter. When there was a general slowing down of the bioproduction processes in the toxic environment, these discharges induced water oligotrophication.
3.4. Ichthyofauna
Fish fauna were found only in Lake Nyudjavr (European vendace
Coregonus albula L., 1758) and in Lake Kumuzhie (trout
Salmo trutta L., 1758; common minnow
Phoxinus phoxinis L., 1758; and burbot
Lota lota L., 1758). The trout found in Lake Kumuzhie was the most valuable representative of the fish community of the five water bodies. The amount of fish caught and the size and weight indices were extremely low. With a mass of 93–207 g (160 g on average) and a length of 19.7–26.7 cm (23.8 cm on average), the fish ages varied from 3+ to 5+ (
Table 3).
The trout found in Lake Kumuzhie were likely to represent the local population spread in this lake and in an adventive lake located upstream. It is necessary to point out that in several water bodies located in areas of industrial pollution [
42,
48], no pathological transformations of inner and outer organs—otherwise typical in the fish in Lake Kumuzhie—were found.
Only isolated vendace were found in Lake Nyudjavr, which could be indicative of unfavourable conditions for their inhabitance. The fish ranged from 1.8 to 5.9 g (4 g on average) in mass and from 6 to 9.2 cm (7.9 cm on average) in length (
Table 3). Comparative analysis showed that the growth rate of vendace in Lake Nyudjavr was lower compared to the similar index of fish from the central part of Lake Imandra. The size-and-weight indices of vendace compared to those from various areas of Lake Imandra were also lower (
Figure 3).
The absence of other kinds of fish in Lake Nyudjavr was found to be due to intensive water body transformation and the long-term impact caused by MP KGMK.
The levels of HM accumulation in fish organs and tissues can serve as an index of the anthropogenic load placed in the five water bodies. We carried out a comparative analysis to determine the HM concentration in the organs of trout found in Lake Kumuzhie (
Table 4).
It was found that the HM and Al accumulated in muscular tissue to a lesser extent in comparison to other analysed organs. The highest concentrations of copper were found in the liver; those of nickel, cadmium, lead, and mercury in the kidneys; and those of zinc and manganese in the gills and bones (
Table 3). The absolute values of copper accumulation, which was a priority pollutant, were high (824.3 µg/g of dry weight). Thus, although Lake Kumuzhie the water body and its catchment area are located rather far from MS KGMK and are protected by the Monchegorskoye tundra massif from the smoke emissions produced by MS KGMK (~3 km), the water body ecosystem was intensively polluted. The comparative analysis of HMs found in the organs and tissues of trout caught in Lake Kumuzhie and different zones of Lake Imandra showed that the concentrations of Cu and Ni, as priority pollutants, were higher in fish near the smelter. The trout from Lake Kumuzhie had the highest copper concentration in their liver (
Figure 4). The nickel and lead in their kidneys had similar regularity (
Figure 4).
It was found that the concentrations of nickel and lead (in the kidneys) and copper (in the liver) exceeded those in fish of this species in Lake Imandra severalfold. At the same time, zinc accumulation in the fish organisms had an opposite tendency (
Figure 4). The highest zinc concentrations were found in fish caught in the southern part of Lake Imandra and the central part of Lake Imandra. We suggested that the zinc was an element resulting from global pollution and its load on the water bodies was not produced by MS KGMK.
4. Conclusions
The ecosystems of the five water bodies developed in conditions with levels of priority pollutants that exceeded those of the conditional background severalfold (HMs and sulphates). Lake Nyudjavr was also exposed to anthropogenic eutrophication. At present, the hydro-chemical conditions are transformed most in Lake Nyudjavr and least in Lake Kumuzhie. The phytoplankton is characterized by depleted species composition and number indices. The taxonomic structure of the phytoplankton differs sharply from that typical of the arctic lakes, with the proportion of the Cyanobacteria and dinoflagellates in the community composition being high. Small lakes located in the vicinity of the sources of pollution differed in their extremely depleted taxonomic composition (less than 10 taxa below the rank of genera). Therefore, in some cases, the abundance could be significant: as a result of the mass development of Pseudanabaena limnetica in Lake Travyanoe, the amount of phytoplankton was as high as 60 thousand ind./L. We found that the size of the lake, combined with intensive water exchange and the presence of excessive biogenic elements (phosphorus and nitrogen compounds), increased algae resistance to toxic pollution. Due to low species-saprobiont numbers, the use of the S-index, traditionally used to qualify water, was limited for the lakes being studied. The S-index used to qualify water affected by smelters should take the toxic load into account, which significantly understates the calculated data. The S-index may function as an indirect indicator of the toxic effect on algal cenosis. Nutrient excess seems to be one of the factors affecting phytoplankton resistivity to toxic load. In this case, the extreme toxicant concentrations slow down the process of eutrophication, which demonstrates low trophic status in the lakes.
The macrozoobenthos communities found in the five water bodies were developed under conditions of toxification, resulting in depleted species composition and depleted biomass in the bottom fauna.
The state of the fish found in Lake Nyudjavr was assessed and showed that, despite the significant transformation of the lake ecosystem, European vendace was found. This could be due to their penetration to Lake Nyudjavr from Lake Imandra via the Nyuduai River. The openness of Lake Nyudjavr and its basin system, its connection with Lake Imandra—the greatest water body in the region—and the sufficient amount of plankton organisms make it possible to reduce high anthropogenic load and promote European vendace penetration into Lake Nyudjavr. Being sufficiently mobile, large animals, fish are able, on the one hand, to avoid highly toxic habitats, and, on the other hand, adapt greatly to them.