Next Article in Journal
Development of Technology for Identification of Climate Patterns during Floods Using Global Climate Model Data with Convolutional Neural Networks
Next Article in Special Issue
Source Identification and Ecological Risk of Potentially Harmful Trace Elements in Lacustrine Sediments from the Middle and Lower Reaches of Huaihe River
Previous Article in Journal
Salinity Shapes the Microbial Communities in Surface Sediments of Salt Lakes on the Tibetan Plateau, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Paleolimnological Perspective on Arctic Mountain Lake Pollution

1
Institute of the North Industrial Ecology Problems of Kola Science Center of RAS, 184209 Apatity, Russia
2
Institute of Limnology of St. Petersburg Federal Research Center of RAS,196105 Saint Petersburg, Russia
*
Author to whom correspondence should be addressed.
Water 2022, 14(24), 4044; https://doi.org/10.3390/w14244044
Submission received: 11 November 2022 / Revised: 8 December 2022 / Accepted: 9 December 2022 / Published: 11 December 2022
(This article belongs to the Special Issue Geochemistry of Water and Sediment III)

Abstract

:
The chemical composition of sediments from the Arctic mountain Lake Bolshoy Vudjavr, situated in the western part of the Russian Arctic zone, was studied. The lake has been under intense anthropogenic load for more than 90 years since the development of the richest apatite–nepheline deposits in the world started. A 27 cm thick sediment core was sampled in the central part of the lake at the maximum depth of 37.4 m. The concentrations of more than 50 elements were analyzed by the mass spectral method, ICP-MS. The lake sedimentation rate established from the change in the content of the radioactive isotope 210Pb was 2.3 mm/yr. The effluent from apatite–nepheline production and atmospheric fallout enrich the sediments of Lake Bolshoy Vudjavr with alkali and alkaline earth metals, N, P, Mn, Fe, Al compounds, rare earth elements, and trace elements (Sb, Cu, Zn, Pb, Bi, Nb, Ta, Th). Analysis of the forms of elements in the lake sediments showed that the studied elements are mainly found in stable fractions—mineral, acid-soluble, and associated with organic matter. The pollution of the sediments of Lake Bolshoy Vudjavr was assessed by the integral index PLI (Pollution Load Index) and CF (contamination factor). The PLI value sharply increased after the “Apatite” Plant had been launched and a large amount of wastewater from the mines had been released into the lake. The highest PLI values were detected in the sediment layers accumulated during the period 1990s–2000s. Sb (18.2), P (10.3), Sr (7.8), and La (6.0) have the maximum CF values among all the studied elements.

1. Introduction

The negative impact of human activities, including the extraction of minerals, on the environment has been known since the Iron Age and the Roman Empire [1,2,3]. Evidence of anthropogenic metal pollution associated with mining and metallurgical activities dates back to at least 5000 years ago [4,5,6]. The pollution level rose during ancient Greek and Roman times [7] as a consequence of the increased mining, smelting, and trading activities [8,9]. Metal production technologies during this period were simple and highly polluting, releasing high concentrations of metals into the atmosphere, hydrosphere, sediments, and soils [10,11,12]. The 20th century saw the peaks of metal pollution in many European countries, due to heavy metal production after World War II, rebuilding and development of cities and industry, and maximum consumption of leaded gasoline around the 1970s [13].
During the 20th century, the Russian Arctic, and Murmansk region in particular, experienced rapid industrialization [14]. The mining industry has been and remains a key factor in the economic development of Murmansk region with numerous mines, processing plants and tailings that pollute soils, rivers, and lakes [15]. The development of the world’s richest apatite–nepheline deposits of the Khibiny alkaline massif [14] has led to the formation of runoff and emissions of the main rock-forming and trace elements and their subsequent deposition in the previously intact mountain environment [16,17,18,19]. Mining, crushing, ore dressing, and tailings discharges are the largest sources of environmental pollution [15,16,20,21,22]. The input of elements in elevated contents relative to natural background contents adversely affects the water quality in the lakes surrounding regional mining centers disrupt the biogeochemical cycles and production processes, which in turn decreases the species diversity of aquatic organisms [15,18,19,21,23,24]. Elevated levels of elements accumulated in the sediments of reservoirs are a sensitive indicator of human activity [15,16,18,22]. Due to the low ability of reservoirs to self-purify [15,18,19,25], the anthropogenicimpact irreversibly changes the natural environment in the conditions of the Arctic: the transformed chemical composition of water and sediments of the reservoirs alters the complexes and groups of hydrobionts living in them [15,18,21,25]. Arctic lakes can provide a long-term perspective on environmental change, including long-range atmospheric transport [25,26,27] and pollutant deposition trends derived from sediment core studies [15,16,17,28]. Studying elemental deposition in the Khibiny lakes is critically important in the light of current pollution and climate change trends reported from other European mountain ecosystems [29]. However, to date, only a limited number of studies have been devoted to the change in the elements accumulation in recent centuries in the lakes of the Russian Arctic and the Khibiny Mountains in particular (for example, [30,31]).
Several studies around the world investigated trace element contamination of water, sediments, soil, and plants as a result of mining activities [15,16,17,19,22,25,26,27,32,33,34]. These studies show that the levels of element contamination around mines depend on the geochemical characteristics, history, and pace of mining of mineral deposits [35,36,37,38,39,40,41]. High pollution levels of terrestrial and aquatic ecosystems persist for many years after the cessation of mining activities, especially if no measures are taken to restore the environment [42,43,44,45,46,47,48].
Trace elements, such as Pb, Sr, rare earth elements (REE), preserved in lacustrine sediments allow for a historical reconstruction of element accumulation in the Khibiny lakes [15,17,18,19,49]. This method offers a retrospective approach in the absence of long-term monitoring data [26,27,31,50]. Although, the Khibiny have a century-long history of pre-industrial [18,19], industrial [14],and modern anthropogenic activity [15,17], the lake sediments in the Khibiny were not studied with the determination of the sedimentation rate and a large list of analyzed elements.
This study of sediments in Lake Bolshoy Vudjavr focuses on the industrial era. In order to reconstruct the history of pollution and assess the intensity of anthropogenic impact on the largest Khibiny lake, the input of a large list of elements into the lake located in close proximity to the mines and receiving their effluents was analyzed. The studied sediments accumulate major rock-forming and trace elements associated with mining activities. Further studies of element accumulation in mountain lakes are required to determine the extent of environmental transformation in the Khibiny massif.

2. Materials and Methods

2.1. Research Area

The Khibiny Mountains are the largest mountain range in Murmansk region. They are located in the center of the region and their area is 1300 km² (Figure 1) [15,51,52]. The tops of the mountains are plateau-like, the slopes are steep with separate snowfields. The highest point of the Khibiny is Mount Judychvumchorr (1200.6 m above sea level). The Kukisvumchorr and Chasnachorr plateaus are located in the center of the Khibiny. Cities of Apatity and Kirovsk are at the foot of the Khibiny.
The climatic conditions in the Khibiny are harsh arctic. The outer slopes of the mountains experience a significant softening effect of the climate of the surrounding plains, and the microclimate of the central part of the massif is much more severe. The mountains are covered with snow from October to June. The average annual temperature in the Khibiny is −2.5 °C. The average temperatures in January-February in the valley areas are –13 °C, in July—no higher than +13 °C. The air temperature drops with altitude decreasing by about 0.5 °C for every 100 m of altitude. The polar night lasts from 40 to 42 days. Cyclones and sudden changes in atmospheric pressure occur frequently. Hurricane winds can blow at speeds up to 50 m/s in the open spaces of the peaks. Summer is short, in the mountains it lasts 60–80 days. The period with an average daily temperature above 10 °C in the foothills lasts about 70 days. Summer also receives the maximum amount of precipitation. The polar day lasts 50 days. Precipitation ranges from 600–700 mm in the valleys to 1600 mm on the mountain plateaus of the Khibiny. Precipitation is distributed almost evenly during the year, a little more in summer, a little less in winter [15].

2.2. Geological Structure

The Khibiny alkaline massif is a large intrusive body of complex shape and composition. According to U-Pb, Rb-Sr, and Sm-Nd dating, the main types of rocks were formed from 380 to 360 million years ago (the Devonian period) [53]. A characteristic feature of the Khibiny massif is a ring (in plan) structure, which has a number of analogies among some other alkaline massifs [15,51,52]. The rock complexes making up the massif have the shape of arcs folded into each other and opening to the east due to the magma intrusion along the alternating ring and cone faults. The Khibiny alkaline massif consists of nephelinesyenite rocks, the rock-forming minerals of which are potassium–sodium feldspars, aegirine NaFe3+(Si2O6) and nepheline (Na,K)AlSiO4 [51]. As a result of the weathering rocks, nepheline is destroyed rather than feldspars, and alkali metal ions Na+ and K+ enter the lake in increased concentrations.
About 500 minerals have been found on the territory of the Khibiny massif, more than 100 of them have been discovered here, 110 are found nowhere else [52]. Many minerals have practical value. Apatite, nepheline, titanite, molybdenite and rinkite are or have been mined. Astrophyllite, aegirine, and eudialyte are used as ornamental stones. The world’s largest deposits of apatite–nepheline ores are located on the territory of the Khibiny massif. The main minerals mined in the Khibiny are apatite, nepheline, sphene, aegirine, feldspar, titanomagnetite. Lovchorrite was previously mined. The massif contains the largest reserves of zirconium raw materials (zircon, eudialyte) and its accompanying hafnium (zircon), which are not currently mined. At the same time, a significant amount of this raw material is currently stored in the tailings of the apatite–nepheline factory.
According to the long-term plan, already in 1944, four mines were to operate in the “Apatite” Plant system: two apatite mines at Kukisvumchorr and Yuksporr, a sphene mine at Yuksporr, and a hibinite mine at Aikuaivenchorr. The ore they mined was to be processed by eight processing plants: two apatite, two nepheline, one sphene, one lovchorrite, one aegirine, and one rare earth. It was supposed to develop the production of phosphorus, to establish cement production [14]. However, the plan was not implemented due to the Second World War.
The following mines are currently operating: in the southeastern part of the Khibiny massif (catchment of the largest lake in the Murmansk region—Imandra)—Kirovsky (Kukisvumchorr and Yuksporr deposits) and Rasvumchorrsky (Apatite Circus and Rasvumchorr Plateau deposits) (Figure 1),in the southwestern part of the Khibiny massif (catchment of the deepest lake in the Murmansk region—Umbozero)—Vostochny (Koashva and Nyorkpakhk deposits) and the recently discovered Oleniy Ruchey (Koashva deposit) (outside of Figure 1). Mining is carried out both underground and open pit. The number of open pit mines is declining and soon the deposits will be developed only by the underground method.
In response to the plans of the North-Western Phosphorus Company, following the opening of the Oleniy Ruchey Mining and Processing Plant, to begin developing a new mine in the center of the Khibiny on the shores of Lake Goltsovoe (to date, forests have been cut down at the site of the future construction), a public environmental movement was formed advocating assigning the status of national park to the Khibiny and prohibiting their further development [54]. In 2018, the Khibiny alkaline massif was assigned the status of the Khibiny National Park (Decree of the Government of the Russian Federation of 8 February 2018 No. 130 “On the establishment of the Khibiny National Park”. Date of access: 28 March 2018).

2.3. Morphometric and Hydrological Characteristics of the Lake

Bolshoy Vudjavr is a dam-type lake, the largest inland water body of the Khibiny (Figure 1). The lake is oval in shape and occupies a depression among the surrounding mountains in the northern taiga zone [15]. The soil cover is sparse, and mountain–tundra vegetation prevails in the catchment area of the lake. Mountains are absent only in the southern part of the lake. Here, the depression is closed by a ridge of sediments of the terminal moraine, cut through by a narrow valley flowing in the southwest of the lake by the Belaya River, which inflows into Lake Imandra. The lake is fed by surface and underground runoff. The absolute elevation of the lake is 312 m, the area of the water surface is 3.9 km2, and the volume of water is 0.0579 km3. The amplitude of fluctuations in the water level in the lake does not exceed 1.1 m (from 312.2 to 311.3 m above sea level). The maximum level is observed in the flood, the minimum—at the end of winter. River Juksporrjok (discharge 0.24–12 m3/s) flows into Lake Bolshoy Vudjavr from the northeast, and River Vudjavrjok flows from the northwest (discharge 0.05–8 m3/s). In the distribution of depths, there is a depression with a depth of 37.6 m in the northeastern part of the lake. The total annual amount of atmospheric precipitation in the catchment area is in the range of 570–1060 mm, 820 mm on average. The temperature regime of Lake Bolshoy Vudjavr is determined by an ice-free period lasting 5–6 months. Studying the thermal regime of the lake showed that in summer there is a direct temperature stratification (from June to October), and in the winter—the reverse. Full spring and autumn homothermy usually sets in June and October, respectively. In winter, the temperature in the surface water layer is 0.4 °C, and is at a depth of 34 m –2 °C. In summer, the temperature of the surface layer on warm days can reach 15 °C.

2.4. History of Development and Anthropogenic Activity

The industrial development of the Khibiny began with the construction of the Murmansk railway during the First World War, which subsequently enabled exploring the geology and developing the mineral resources of the region. The geological surveys under the guidance of Academician A.E. Fersman in the early 1920s resulted in the discovery of the world’s largest apatite–nepheline deposits with huge reserves of phosphorus raw materials in the Khibiny mountain range [51]. In 1929, the Production Trust “Apatite” was organized to develop the mineral resources of the Khibiny. Since 1931, processing apatite ore began at the processing plant (ANOF-I) on the shore of Lake Bolshoy Vudjavr; since 1934 the second stage was launched at this enterprise. From that time on, sewage began to flow into Lake Bolshoy Vudjavr and further along River Bolshaya Belaya without treatment. Municipal sewage from City of Kirovsk was also discharged without treatment. The apatite industry steadily increased its capacity and by 1938 produced more than 1 million tons of apatite concentrate [14]. The waters of the Juksporrjok River and further Lake Bolshoy Vudjavr were enriched with nitrogen and phosphorus compounds and silicic acid [55]. There is a sedimentation basin for wastewater from the Kirovsky and Rasvumchorrsky mines in the northeastern part of the lake, which comes with the flows of Rivers Juksporrjok and Saamskaya, fenced off from the main lake water area by a dam (Figure 1). Population growth and inefficient treatment facilities have led to the pollution of the lake by municipal wastewater from City of Kirovsk, as well as the Polar Alpine Botanical Garden. The ichthyofauna of Lake Bolshoy Vudjavr, numbering 5 species of fish, had practically disappeared by 1931 due to the pollution of the Bolshaya Belaya River, along which fish ascended to spawn from Lake Imandra. At present, the lake is inhabited by Arctic char and nine-spined stickleback [23], as well as European smelt, which has been actively migrating from Lake Imandra in recent years [24].
The chemical composition of the lake sediments depends on the natural (geochemical, hydrological, climatic, hydrobiological) features of the watershed and the lake itself, as well as direct and aerotechnogenic pollution of mining and metallurgical enterprises (“Apatite” Joint-Stock Company and the Kola Mining and Metallurgical Company (KMMC), a subsidiary of “Norilsk Nickel” Joint-Stock Company), as well as atmospheric transboundary transfers of pollutants. The influx of sewage and dust emissions from mines and ANOF-I of “Apatite” JSC, containing mainly nepheline and apatite Ca5(PO4)3(F,OH,Cl), enrich the water and sediments of the lake with Na, K, Ca, P, Al, Sr, F, and other elements that are part of apatite–nepheline ore and products of their processing. The Khibiny deposits are represented mainly by fluorapatite, but there are also other varieties of it. Apatite also contains impurities of Mn, Fe, Th, rare earth elements (REE), calcium carbonate CaCO3 (carbonateapatite), and other impurities [51].
Atmospheric emissions from metallurgical plants can be transported over a distance of several tens or even hundreds of kilometers [56,57]. The Severonickel Copper andNickel Plant (Monchegorsk site of KMMC), located 45 km northwest of Lake Bolshoy Vudjavr, specializes in the production of electrolytic Ni and Cu, cobalt concentrate, sulfuric acid, and other products. Minerals such as pentlandite (Fe,Ni)9S8, chalcopyrite CuFeS2, cobaltite (Co,Ni)AsS, and nickeline NiAs are involved in the process of enrichment and metallurgical processing [58]), so all the above mentioned elements are present in high concentrations.

2.5. Research Methods, Sampling and Analysis of Sediments

2.5.1. Sediment Sampling

Sediment cores of 23 and 27 cm thickness were sampled on April 26, 2018 in the central part of Lake Bolshoy Vudjavr at a depth of 37.4 m (Figure 1). Sampling was carried out using the open gravity corer Limnos from the ice cover, allowing to extract stratified sediment core samples up to 60 cm long and separate these into 1 cm or thicker layers. After sampling, undisturbed sediment cores were transported to the laboratory, where they were divided into layers of 1 cm to determine the total content of elements and into layers of 5 cm (0–5, 5–10, 10–15 and 18–23 cm) to analyze the element fractions.Then, the samples were placed in a refrigerator, in which they were stored until dry at a temperature of about 4 °C according to the methodological guidelines [59,60,61,62,63]. For further study, the sediment samples were dried to an air-dry condition at room temperature and then to acompletely dry condition in an oven at a temperature of approximately 105–110 °C [64].

2.5.2. Granulometric Analysis of Sediments

The particle size distribution of sediments was determined using a multifunctional particle analyzer of the LS13 320 series (Beckman Coulter, Life Sciences, Indianapolis, IN, USA). This analyzer determines the particle size of sediments by laser diffractometry in water based on the physical principle of scattering of electromagnetic waves of different lengths. The technical features of the device allow the analysisof particles ranging in size from 0.04 µm to 2.0 mm in accordance with the ISO 13320-1 Standard.

2.5.3. Sediment Dating

When the collected sedimentsaccumulated in Lake Bolshoy Vudjavr, accordingly, the sedimentation rates were calculated from the change in the content of the 210Pb isotope with depth in the core (nonequilibrium 210Pb method) [65]. The model of constant initial concentration estimating only the average rate of sedimentation was used. The contents of 210Pb in individual layers of sedimentary core were determined from the alpha activity of its granddaughter isotope 210Po, which is in radioactive equilibrium with it and isolated from samples by the standard radiochemical method on nickel disks [66]. Sediment samples were dried at room temperature to a constant weight, followed by grinding in an agate mortar to a powdery state. The sample was placed in a solution of 15 mLof 7M HNO3 and kept for a day, and then was evaporated to dryness and treated with 15 mL of concentrated HF to completely dissolve the silicate fraction of the precipitate. The procedure was repeated thrice, each time evaporating the solution to dryness. At the next stage, the resulting residue was converted into the nitric acid form and evaporated again, after which the sample was dissolved in 100 mLof 0.1 M HCl with the addition of 2 mL of concentrated H2O2 with gentle heating until oxygen evolution ceased. By adding 100 mg of ascorbic acid, trivalent forms of iron were reduced and H2O2 residues were decomposed. Further, a pre-prepared nickel disk was placed into the resulting solution, fixed in a special fluoroplastic cassette, which ensures the release of 210Po on one side of the disk, and kept under stirring with a magnetic stirrer and 65 °C for 3 h. Further, the solution was poured off, and the removed disk was washed with distilled water and dried at 100 °C. The measuring targets prepared in this way were placed in radiometric equipment (an ORTEC Alpha-Duo spectrometer with an ion-implanted silicon detector) to obtain alpha spectra, which were used to measure the 210Po activities in each sample, and from them, the 210Pb concentrations were determined.

2.5.4. Determination of Loss on Ignition

Loss on ignition (LOI), which characterizes the content of organic matter in sediments (organic compounds are burnt in the temperature range from 100 to 500 °C), was determined by the gravimetric method after igniting an absolutely dry sample of each layer at a temperature of 550 °C in a muffle furnace to constant mass [60].

2.5.5. Chemical Analysis of Elements

The content of the main rock-forming elements (oxides of Al, Ca, Fe, K, Mg, Mn, Na, P, Ti, Si) in lake sediments was determined using an ARL ADVANT’X X-ray fluorescence spectrometer (Thermo Fisher Scientific). Before X-ray fluorescence analysis, the sediment samples were dried in an oven at a temperature of 105 °C to a constant weight.
To estimate the total concentrations of trace elements in lake sediments, open acid decomposition of sediment samples was carried out using HF, HNO3, and HCl. Analytical weights 0.1 g of samples were used for the analysis. As a blank sample, a mixture of decomposition acids was used that underwent the same sample preparation procedure as sediment samples and a standard (control) sample—the chemical composition of the bottom silt of Lake Baikal BIL-1—GSO 7126-94. The concentrations of trace elements in sediment samples were determined by mass spectrometry on an XSeries-2 ICP-MS instrument. This technique is described in detail in a previously published article [63].
All laboratory work was based in the analytical center of the Institute of Geology of the Karelian Scientific Center of the Russian Academy of Sciences (Petrozavodsk). Programs Statistica (version 10), Microsoft Office Excel 2016, ArcGIS for Desktop 10.4.1 were used for statistical data processing and illustrations.

2.5.6. Determining the Chemical Fractions of Elements

To determine the fractions of various elements, a number of sequential extraction methods for soils have been adapted for sediments [67,68,69,70]. The schemes differ in the number of isolated fractions, as well as in the reagents used. The method for sequential extraction of element forms was used in this work, including the identification of the following:
-
Water-soluble fractions (reagent H2O).
-
Available (exchangeable) fractions (reagent NH4CH3COO, pH 4.8). The elements can desorb in case of changing of pH or ionic composition of water.
-
Fractions bound to Fe and Mn hydroxides (reagents 0.04 M NH2OH+HCl in 25% CH3COOH). Iron and manganese oxides exist in sediments as nodules, concretions, cement between particles. They are excellent scavengers for metals. The fraction is thermodynamically unstable under anoxic conditions (low Eh).
-
Fractions bound to organic matter (reagents 0.02 M HNO3 + 30% H2O2 and 3.2 M NH4CH4COO in 20% HNO3). Trace metals may be associated to various forms of organic matter. The complexation ability of humic substances in natural ecosystems is studied previously [71]. Organic matter can be degraded under oxidizing conditions leading to a release of trace metals in water.
-
Acid-soluble (residual) fractions (reagent HNO3); the most stable (mineral) fraction. The metals are associated with residual compound and not available to living organism. Dissolution of the fraction is possible only under a very strong anthropogenic impact on the reservoir.
-
Mineral (silicate) fractions obtained by deducting the combined concentration of all the above fractions from the total concentrations. The metals bounded with the crystal structure of silicate minerals are not expected to be released in water over a reasonable time span under the conditions normally encountered in nature.
We decomposed the sediment samples to determine the total concentrations of heavy metals in the sediment in the core samples, and subdivided them into 5 cm layers by an acid breakdown using HF, HNO3, and HCl in an open system.

2.5.7. Intensity of Lake Pollution

L. Håkanson’s method [64] was chosen to assess the contamination of the sediments of Lake Bolshoy Vudjavr. To identify the intensity of contamination by substances, the values of the contamination factor (CFi) were determined:
CFi = Cmi/Cni,
where Cmi is the content of the element in the layer m (cm) of sediments, Cni is the pre-industrial value for this element, determined as the background value for the sediments of Lake Bolshoy Vudjavr in its deepest part (24–27 cm), deposited before the start of the work of the “Apatite” Plant, i.e., until 1929. This approach adhered to the following classification: CFi < 1—low contamination factor (indicating low contamination of sediments by the studied element); 1≤CFi < 3—moderate contamination factor; 3≤CFi < 6—significant contamination factor; CFi ≥ 6—high contamination factor.
The ecological interpretation of the geochemical data was based on calculating the integrated Pollution Load Index (PLI) [72], which takes into account the content of all the studied elements in the lake sediment core:
PLI = (CF1 * CF1* CFn)1/n,
where CFn is the contamination factor (the ratio of element concentration to the background value of this element). Depending on the value of the PLI, the pollution level can be assessed as low (PLI ≤ 1), moderate (1 < PLI < 2), high (2 < PLI < 3), and extremely high (3 < PLI).

3. Results and Discussion

3.1. Granulometric Composition of Sediments

Analysis of the granulometric composition of the sediments of Lake Bolshoy Vudjavr revealed that fractions smaller than 0.05 mm (88.4% of the total mass of the sample) prevail in the sediments and fractions in the range from 0.05 to 2 mm account for 11.6% of the total mass of the sample. Thus, the studied sediments can be attributed to the silty type of lacustrine sediments [73]. The particle size distribution was previously assessed by the pipette method based on the difference in the settling rates of particles of different sizes in water, the prevalence of smaller size fractions (clayey less than 0.005 mm) in the sediments of Lake Bolshoy Vudjavr was revealed and the deposits were classified as clayey [17]. This is probably due to different methods for analyzing the particle size distribution of sediments in these studies.

3.2. Dating of Sediments

As mentioned above, the distribution of the 210Pb isotope along the profile of the sampled sediment core was used to date the sediments of Lake Bolshoy Vudjavr. The decrease in the 210Pb content as a whole has a consistent character, reaching a basic value at a depth of about 22 cm, which almost does not change further with depth (Figure 2). This allows concluding that the excess of 210Pb (210Pbex. in Figure 2) in the upper seven studied layers of the lacustrine sediment core is not supported by radioactive equilibrium, and the concentrations determined for the three lower layers are considered to be maintained by the equilibrium in the 238U series. Therefore, their weighted average value (39.11 ± 0.59 Bq/kg) can be applied to calculate the concentration of 210Pbex. in the overlying layers. The data on the content of the 210Pbex. isotope in the layers from 0 to 20 cm are well approximated by the exponent (r = −0.852), which is obviously due to the decrease in its concentration according to the law of radioactive decay with a half-life of T1/2 = 22.2 years at a relatively constant sedimentation rate, which means that it allows applying the calculation model of constant initial concentration [65]. Calculations using this model showed that the sediment layer within the horizon of 2–20 cm accumulated over 75.23 ± 2.11 years, and the average sedimentation rate over the calculated period of time can be estimated as 2.26 ± 0.06 mm/yr. Notably, the obtained sedimentation rate is close to the sedimentation rates (2.5–3.0 mm/yr) in small lakes in the taiga, forest–tundra, and steppe landscapes of Siberia [74]. The similar sedimentation rate (1.5–3 mm/yr) was obtained in Lakes Imandra and Kuetsjärvi (northwestern Murmansk Region) affected by apatite–nepheline and copper–nickel production runoff [16,75]. At the same time, the calculated sedimentation rate in this work is higher than the average sedimentation rate (0.3-1.25 mm/yr) in the lakes of Northern Fennoscandia and Murmansk region, which are not directly polluted by industrial enterprises [26,62,76,77], and some lakes in the southern part of the Republic of Karelia [78].

3.3. The Main Rock-Forming Elements

The Khibiny massif is the world’s largest intrusion of nepheline syenites, melteigite–urtites, and apatite–nepheline rocks with an area of more than 1300 km2 [52]. The listed rocks, as well as the main part of other rocks that make up the Khibiny, are composed mainly of aluminosilicate minerals. Therefore, in the chemical composition of the studied sediments of Lake Bolshoy Vudjavr, oxides of Si (average content 44%), Al (16%), Ca (5.9%), and Fe (5.3%) prevail among the main components (Figure 3). In addition, according to the LOI, it can be concluded that the content of total organic matter in the sediments is 17.4%. According to the classification given in the work of D.A. Subetto [73], these deposits should be referred to siliceous siltstone sapropels in terms of particle size distribution, content of silicon and organic material. The organic material content in the studied sediments can be considered rather low for the lakes of the taiga zone of Russia, which can be explained by the geographical location of the reservoir—inside the mountain range with sparse soil-vegetation cover and significant filtration capacity of the clastic material of the upper rock layers, where the main material coming from the watershed into the lakes is represented by a terrigenous component, as well as by the influx of a large amount of suspended mineral particles with apatite–nepheline production effluents from the mines and the processing plant. The organic matter content in the sediments of lakes located in the plain (boggy or forest) regions in the north of Russia most often varies from 30 to 80% [79]. A similar content of organic material was recorded in the sediments of lakes in the Murmansk region [80].
The macroelements composing apatite Ca5(PO4)3(F,Cl,OH) and nepheline (Na,K)AlSiO4 are found in high concentrations in the surface sediment layers of Lake Bolshoy Vudjavr. On the one hand, the Si oxide content decreases along the core of the studied sediments of Lake Bolshoy Vudjavr from 51–52% in the lower layers to 38–39% towards the sediment surface (Figure 4, Table 1). On the other hand, there is a noticeable increase in the concentrations of P (by an order of magnitude) (Figure 5), Ca (by a factor of 3.5), K, Na, and Ti (the last three elements by more than a factor of 2) oxides, which, according to early works on this lake [18,19] is due to the extraction and processing of apatite–nepheline raw materials by the “Apatite” Plant since the 1930s. The organic material content (LOI) in the surface sediment layers decreases (Figure 4), which is associated with the suppression of biological production in the lake due to severe pollution and restructuring of hydrobionts, as well as the influx of a large amount of suspended mineral particles with wastewater from apatite–nepheline production [18,19].
A significant increase in P concentrations in the studied sediment core begins in the 20–21 cm layer, which corresponds to the beginning of the 1930s, and steadily increases from 0.3 to 3.4% in the upper layers of the lacustrine sediments (Figure 5). The noticeable increase in the P content can be traced in the upper part of the core after the collapse of the Soviet Union and the onset of the economic crisis after the 1990s, although the production at the “Apatite” Plant decreased by more than three times. Perhaps this is due to the decrease in the enterprise’s expenses for environmental protection during this period and the decrease in the efficiency of sedimentation of suspended matters in a fenced dam on the lake.
The content of the main rock-forming elements in the sediments of Lake Bolshoy Vudjavr exceeds their concentrations in the surface and background layers of the sediments of Lake Imandra (with the exception of Si, Fe, and Mn), which is associated with the influence of the runoff of the “Apatite” Plant (in the surface sediment layers) and the geochemical features of the territory catchment area (background sediment layers) in Lake Bolshoy Vudjavr (Table 1). Furthermore, the content of many main elements (Si, Al, Ca, K, Mg) in the surface part of the sediments is comparable to their content in the tailings of apatite–nepheline ore processing, and the P content is even higher than in the tailings (Table 1).
Based on studying the content of the main rock-forming elements, the thickness of lake sediments was divided into threelayers corresponding to the stages of formation of the chemical composition of sediments—the pre-industrial period (Figure 4,layer I), after the start of development of apatite–nepheline deposits and a constant increase in production (layer II), the current state of the lake (layer III). The vertical distribution of elements in the sediments of the lake will be further described according to these stages of development of the enterprise and sedimentation.

3.4. Trace Alkali and Alkaline Earth Metals

Trace alkali and alkaline earth metals were called metals, which, unlike the main rock-forming metals (Na, K, Mg, Ca), are found in the sediments of Lake Bolshoy Vudjavr in much smaller quantities, but in increased contents relative to other lakes in Murmansk region, both in surface and in the background layers (Table 2). The concentrations of trace alkali and alkaline earth metals in the sediment core of Lake Bolshoy Vudjavr have sharply increased since the extraction and processing of apatite–nepheline raw materials started in 1930s (Figure 6). The content of trace alkali and alkaline earth metals also increased after the collapse of the Soviet Union and the economic crisis of 1990, which is associated with an increase in production in the 2000s. The content of all alkali and alkaline earth metals (r = 0.87–0.97) was closely correlated, exception for Li, whose distribution differed from other alkali and alkaline earth metals. Nevertheless, even this lightest metal sharply increased its content in the early 1930s, when the development of apatite–nepheline ores started (Figure 6). Lithium isomorphically replaces potassium in widespread rock-forming minerals, and it is also found in some minerals of the Khibiny massif, for example, manganneptunite KNa2LiMn2+2Ti4+2[Si8O24]. The largest increase in the content in the surface layers relative to the background values was noted for Sr (by a factor of 8), while the concentrations of other trace alkali and alkaline earth metals increased not so significantly (from 1.5 to 2.7) (Figure 6, Table 2).
Feldspars (rock-forming minerals of the Khibiny rocks) in addition to potassium, sodium, and calcium also contain other alkali and alkaline earth metals. Feldspars are used as a raw material to extract rubidium. Potassium-barium feldspars (hyalophanes) are known. The rocks of apatite–nepheline deposits contain a large number of minerals made up of trace alkali and alkaline earth metals, for example, varieties of the mineral pyrochlore (Na,Ca)2Nb2O6(OH,F) widespread in the Khibiny—bariopyrochlore (Ba,Sr)2(Nb, Ti)2(O,OH)7, strontium pyrochlore Sr2Nb2O7, as well as gadolinite Ce2Fe2+Be2[SiO4]2O2, belovite Sr3Na(Ce,La)[PO4]3F, barite Ba[SO4], ankylite Sr(Ce,La)[CO3]2(OH)·H2O, strontianite Sr[CO3], burbankite (Na,Ca)3(Sr,Ba,Ca,REE)3[CO3]5, and many others [52].
The main ore mineral fluorapatite contains a large amount of strontium (SrO up to 10%) and other alkali and alkaline earth metals. There is even a variety of fluorapatite, strontioapatite Sr3Ca2[PO4]3F, in which the content of strontium oxide is about 50% [52]. The correlation coefficients of trace alkali and alkaline earth metals (except for Li) with phosphorus, the main element of the mineral apatite, are very high (r = 0.91–0.98, except for Li). Therefore, the main source of these metals in the sediments of Lake Bolshoy Vudjavr is probably mining and processing effluents of apatite raw material, although numerous trace alkali and alkaline earth metal minerals recovered from the mines during ore extraction also supply a certain amount of these trace metals to the sediments.

3.5. Rare Earth Elements (REE)

REE are a group of 17 elements including scandium, yttrium, and the lanthanides (from lanthanum to lutetium). The anthropogenic factors (industrial emissions and effluents, as well as dust from construction, deterioration of residential and road structures, weathering of soils, and Quaternary sediments in urban areas) play an important role in the accumulation of REE in the lake sediments in industrial areas and urban water bodies [83]. Environmental pollution with REE is associated with the activities of mining companies, industrial enterprises, as well as with the combustion of coal containing impurities of various metals, including REE, at energy enterprises [84].
An increase in the concentrations of all REE from two to six times is recorded in the upper sediment layers of Lake Bolshoy Vudjavr (the largest increase was noted for light REE, for example, for La—6.0), which results from the extraction and processing of apatite–nepheline ores (Figure 7). A similar distribution of REE was recorded in the sediments of Lake Imandra in the zone of wastewater inflow of the “Apatite” Plant [49,81,82], but sediment layer with elevated REE contents in Lake Imandra is two times smaller (10 cm), which indicates a lower sedimentation rate in it. There is an increase in the content of REE in the mid-1960s, and then a decrease in the concentration of REE, which is probably associated with the commissioning of the ANOF-I tailing dump and a decrease in the anthropogenic load on the lake. A large amount of REE is contained in fluorapatite (Y2O3 up to 3.3%, La2O3 up to 1.7%, Ce2O3 up to 4.1%, Nd2O3 up to 1.4%, etc.) [52]. There are also numerous minerals containing REE [52], for example, lovchorrite (Na(Ca,Na)2(Ca,Ce)4TiO2F2(Si2O7)2 (REE oxides TR2O3 up to 17%), in which calcium can be replaced by Sr (up to 3.6%), Ce and U, loparite (Ce,Na,Ca)2(Ti,Nb)2O6) (Ce2O3 up to 20%, La2O3 up to 13%, Nd2O3 up to 3.7%), gadolinite Ce2Fe2+Be2[SiO4]2O2 (Ce2O3 up to 57%, La2O3 up to 18%, Nd2O3 up to 5%), nordite Na3Sr(Ce,La)Zn[Si6O17] (Ce2O3 up to 20%, La2O3 up to 10%), rinkite Na2Ca4(Ce,Y)Ti[Si2O7]2(F,OH,O)4 (Ce2O3 up to 11%, La2O3 up to 5%, Nd2O3 up to 4%), and cerite (Ce,La,Ca)9(Fe3+,Mg)[SiO4] 6[SiO3(OH)](OH)3 (Ce2O3 up to 30%, La2O3 up to 24%, Nd2O3 up to 5.6%).
The distribution profiles of all REE are very similar (except for Sc, the content of which decreases from background to surface layers), and the vertical distribution of the light REE content in the sediments of Lake Bolshoy Vudjavr are shown in Figure 7 as an example. A similar distribution of all REE (with the exception of Sc) is indicated by the correlation coefficient between REE, which in most cases, is equal to 1.00. Scandium with all REE shows a correlation coefficient ranging from –0.91 to –0.95.
The content of REE in the surface sediment layers of Lake Bolshoy Vudjavr is the highest of all known values for the lakes of Murmansk region (Table 2) and the Republic of Karelia [83], which indicates the inflow of effluents with a high content of REE from JSC “Apatite” into the lake. The background REE contents in the sediments of lakes Bolshoy Vudjavr and Rybachy are comparable (Table 2), which indicates the similarity of the natural geochemical setting of these two lakes in Murmansk region. The highest accumulation of REE took place in the sediments of lakes in the tundra zone, while the lowest accumulation of REE was noted in the lakes of the forest zone (taiga) [84]. Lake Bolshoy Vudjavr is located in the mountain tundra zone, and Lake Rybachy is located in the flat tundra zone. In addition, REE-Th-U mineralization of pegmatoid granites within the Litsevsky uranium ore region, including the Rybachy and Sredny Peninsulas [85], is manifested in the elevated REE contents in the background sediments of Lake Rybachy.

3.6. Trace Elements That Reduce the Content throughout the Sediment Core

The Khibiny alkaline massif contains in its depths the largest reserves of zirconium raw materials (zircon Zr[SiO4], eudialyte Na12Ca6Fe2+3Zr3[Si3O9]2[Si9O24(OH)3]2, baddeleyite ZrO2 [52]) and its accompanying hafnium which are currently not processed. In total, 3.7 thousand tons of eudialyte concentrate were planned to be produced at the “Apatite” Plant during the five-year plan 1933–1937 [14]. The increase in the content in the distribution profile of Zr and Hf in the sediments of Lake Bolshoy Vudjavr was recorded in the 1930–1940s (Figure 8), which indicates the consequences of the extraction of zirconium raw materials at that time. The content of Zr and Hf decreases to the sediment surface since the extraction of zirconium raw materials was terminated. The content of Zr and Hf correlates very closely with each other (r = 0.97).
The Mo minerals molybdenite MoS2 molybdite MoO3, wulfenite PbMoO4, and ferrimolybdite Fe2(MoO4)2·8H2O are fairly widespread in the Khibiny mountain range [52]. The conditions to find molybdenite were described in detail by A.N. Labuntsov, who discovered the Takhtarvumchorr molybdenite deposit in 1927 [80]. A decision was made in 1933 to build a molybdenum mine and a molybdenum processing plant [14]. Molybdenite was mined in the early 1930s at the Takhtarvumchorr mine, and then the mining stopped due to “insignificance of reserves and poverty of ore”, as well as due to the discovery of more promising deposits of molybdenum ores in Kazakhstan and the Caucasus [86]. There is a significant increase in the content of Mo in the 1920s, which is most likely associated with the exploration and development of molybdenum ores from the Takhtarvumchorr deposit (Figure 9).
The content of Mo closely correlates with that of As (r = 0.86), and their distributions in the sediments of Lake Bolshoy Vudjavr completely coincide (Figure 9). The rocks of the Khibiny alkaline massif contain a number of minerals containing arsenic—arsenopyrite FeAsS, cobaltite CoAsS, lellingite (Fe,Co)As2, and safflorite (Co,Fe)As2 [52]. Probably, these minerals are also common in the rocks of molybdenum ores of the Takhtarvumchorr deposit.
The concentrations of Cr, Ga, and Tl in the sediments of Lake Bolshoy Vudjavr significantly decreased after the development of apatite–nepheline deposits started, which is associated with their lower content in the ore and tailings (Table 2), and, accordingly, in the effluents of mines and the concentration plant.

3.7. Trace Elements That Reduce the Content in the Upper Sediment Layers

An increase in the concentrations of almost all the studied elements is observed in the vertical distribution of elements in the sediment core of Lake Bolshoy Vudjavr, starting from the layer of 20-21 cm, which indicates the period of the beginning of mining and processing of apatite raw materials (1930s). At the same time, a large group of elements stands out, the content of which has been decreasing in recent decades for various reasons (Figure 10).
The 10–15-fold excess of Sb in the uppermost layers of the studied sediments is especially distinguished, which makes this element one of the most potentially dangerous for the urban lake ecosystem. In general, this is a common phenomenon for urban lakes in Murmansk region and Karelia. Thus, Sb was noted as a priority pollutant for small lakes in the cities of Petrozavodsk and Murmansk [87,88], for example, a 30-fold increase in the content of Sb in the surface sediment layer compared to the background content (Table 2). Antimony is wellknown to accumulate in different environments as a result of the long-range transport of pollutants [89]; however, in the case of Lake Bolshoy Vudjavr, there may be a strong influence of a local coal-fired thermal power plant located 10 km from the reservoir. The enterprise began operation in the late 1950s and is still operating, which can be seen from the Sb distribution in the sediment core of Lake Bolshoy Vudjavr—an increase in the concentration of this metal is recorded just in the middle of the 20th century (Figure 10). As Sb is one of the most carboniferous metals [90], the coals used at the local thermal power plant are likely to have an increased content of this element.
Contamination of the studied sediments with V and Ni is also directly related to the heating supply of Kirovsk,until a boiler house began operating in the city in 2013, using masut as fuel. The content of V and Ni in the upper sediment layer reached the background values (Figure 10). Masut ash is known to be enriched in these elements, which enter the atmosphere during fuel combustion and then are washed into soil, water, and reservoir sediments [91,92]. Currently, the boiler house is not functioning, and the heating supply of Kirovsk is provided by the coal-fired thermal power plant in Apatity. An increase in the content of V and Ni was recorded in the sediments of Murmansk, which is associated with the influence of masut boilers and thermal power plants of the city, which began to operate on this type of fuel in the 1960s (coal was previously used) [93]. A distribution similar to these metals in the sediments of Lake Bolshoy Vudjavr was noted for W (Figure 10), which closely correlates with V and Ni (r = 0.92). Abnormally high contents of W and Sb were found in the coal samples obtained from the deposits of the Northern Urals [94].
The content of U in the sediments of Lake Bolshoy Vudjavr sharply increased in the 1970s and 1980s, possibly due to the peaceful nuclear explosions of Dnepr-1 (4 September 1972) and Dnepr-2 (27 August 1984) at the Kuelporsky mine to increase the fragmentation of the ore body and reduce the cost of extracting the ore. It is quite probable that a certain amount of the nuclear explosion products reached the lake, although U is found in elevated concentrations in fossil fuels, including coal [94]. Limited amounts of U are found in the minerals of apatite–nepheline deposits, for example, rinkite, monazite, lovchorrite, and uranium pyrochlore (U,Ca,Ce)2(Nb,Ti,Ta)2O6(OH,F) [52]. Since the 1990s, there has been a constant decrease in the U content, and its level in the surface sediment layers is comparable to the background (Figure 10).
Aerotechnogenic transfer of substances also contributes to the pollution of Lake Bolshoy Vudjavr. These substances come from the plant processing of copper–nickel ores of Kola MMC, located in Monchegorsk, 45 km from the studied reservoir [95]. First of all, this applies to such elements as Cu, Ni, Co, and Zn. Moreover, the impact of emissions from Kola MMC on the environmental objects located hundreds of kilometers from the enterprise is wellknown [96]. The patterns of metal accumulation in the lake catchment area are different. For example, it was found that the input of Ni in the wastewater from the Severonickel Plant is eight times higher than the input from the catchment area (i.e., in the form of aerotechnogenic fallout), while the Cu input from the catchment area is 7.5 times higher than that from the wastewater [16]. This was reflected in the distribution of Cu and Ni in the sediment core of Lake Bolshoy Vudjavr, since the Cu concentration in the upper sediment layer exceeds the Ni concentration by a factor of seven, and Ni reduces the content after 1990, while Cu concentrations continue to increase during this period (Figure 10 and Figure 11). A similar distribution pattern was noted for Co and Zn—if the first element decreases its content in the upper sediment layers, the latter has been stable for the last 30 years (Figure 10 and Figure 11).
Metals such as chalcophile Pb, Cd, and Hg classified as global pollutants [97], actively accumulate in lake sediments as a result of long-range atmospheric transport, which is confirmed by numerous studies in Murmansk region and neighboring regions [98,99,100]. At the same time, Pb from vehicle exhaust gases could also actively enter the ecosystem of Lake Bolshoy Vudjavr, since the studied reservoir has an urban status. From the 1930s to the early 2000s, tetraethyl lead is known to have been added to gasoline to improve fuel quality [101]. The decrease in Pb concentrations in the upper sediment layers of Lake Bolshoy Vudjavr is most likely due to the ban on the production of leaded gasoline in Europe and in Russia in particular. A similar pattern of Pb distribution can also be noted in the sediments of lakes in other cities of Murmansk region and the Republic of Karelia [88,93]. Pb, along with Sb, is also usually ranked among the priority pollutants of small lakes in the urbanized environment of the north of Russia.
The vertical distribution of Hg in the sediments of Lake Bolshoy Vudjavr clearly shows a surface maximum of up to 1.1 μg/g at a depth of 13–16 cm (Figure 10), which is more than 20 times higher than the average background concentrations of this extremely toxic chalcophile element in the sediments of lakes in Murmansk region [21,98]. A similar distribution of Hg was found in the sediments of Belaya Guba of Lake Imandra [102]. A comparison with the vertical profiles of other pollutants leads to the conclusion that the increase in Hg contents in lakes Bolshoy Vudjavr and Imandra coincides with the onset of the wastewater inflow from “Apatite” JSC (P, Al, Ca, Na, Sr, K) and “Severonickel” (Ni, Cu, Co), but the maximum Hg concentrations were recorded earlier in time than those of the main polluting heavy metals (Ni and Cu). It can be assumed that the maximum influx of Hg occurred in the 1930s–1950s (Figure 10). During blasting operations at the first stages of development of apatite–nepheline deposits, fulminate mercury (mercury salt of fulmic acid Hg(CNO)2) was used in blasting caps. During the Second World War, on the basis of the “Apatite” Plant, there was a workshop for manufacturing incendiary phosphorus bombs, in which the fulminate mercury was also used as an igniter.
Another chalcophile element in the sediments of Lake Bolshoy Vudjavr, tin, shows a distribution pattern similar to Pb and Cd (the correlation coefficient of Sn with these metals is 0.83 and 0.88, respectively)—a rather sharp increase in the early 1930s, and a decrease after the 1990s (Figure 10). At the same time, high values of the correlation coefficient Sn were recorded with all alkali (except Li) and alkaline earth metals and REE (r>0.9). Sn possibly enters the lake sediments both as a part of mine and concentrating plant effluents, and as part of emissions from fossil fuel combustion and copper–nickel production. The latter conclusion is supported by a more than 50-fold increase in the Sn content in the sediments of Lake Rybachy, located 50 km from Pechenganickel plant, and a 15-fold increase in Sn in the sediments of Lake Severnoye compared to the background concentrations (Table 2) [27,83].

3.8. Trace Elements That Increasethe Content in the Upper Sediment Layers

Minerals containing niobium are widely distributed in the rocks of the Khibiny alkaline massif. Niobium is always accompanied by tantalum, and the content of both metals is very closely correlated (r = 0.98) in the sediments of Lake Bolshoy Vudjavr. The close chemical properties of Nb and Ta determine their joint presence in the same minerals and participation in common geological processes. A number of titanium-containing minerals: titanite (sphene) CaTi[SiO4]O, perovskite CaTiO3, ilmenite FeTiO3, ilmenorutil (Ti,Nb,Fe)3O6, rutile TiO2, and rinkite Na2Ca4(Ce,Y)Ti[Si2O7]2(F,OH,O)4, are minerals common in the Khibiny [52]. Niobium replaces Ti in titanium-containing minerals, as evidenced by the very close correlation of Ti with Nb and Ta (r = 0.97). The form of finding niobium can be different: scattered (in rock-forming and accessory minerals of igneous rocks) and mineral. In total, more than a hundred minerals containing niobium are known, and many of them are found in the Khibiny massif. Of these, the following minerals are of industrial importance in the Khibiny: pyrochlore (Na, Ca)2Nb2O6(OH, F) (Nb2O5 up to 62%), loparite (Ce,Na,Ca)2(Ti,Nb)2O6) (Nb2O5 up to 20%), lovchorrite (Na(Ca, Na)2(Ca, Ce)4TiO2F2(Si2O7)2 (Nb2O5 up to 2%, ThO2 up to 0.5–1%, and rare earth oxides TR2O3 up to 17%) [52]. The concentrations of Nb and Ta in the surface and background layers of the sediments of Lake Bolshoy Vudjavr are 1–2 orders of magnitude higher than the content of these elements in lakes outside the Khibiny (Table 2, Figure 12), which indicates the natural enrichment of this mountain massif with these metals.
Fifty thousand tons of sphene concentrate at the “Apatit”Plant were planned to be produced during the five-year plan of 1933–1937 [14]. Between 1934 and 1939, a lovchorrite mine was developed, and during this time, 19.4 thousand tons of lovchorrite ore were mined [86]. The content of Nb and Ta in the sediments of Lake Bolshoy Vudjavr, as well as Ti, increased sharply in the early 1930s since the beginning of the development of apatite–nepheline deposits and then only increased.
Thorium is present in small amounts (from fractions to a few %) in almost all minerals containing Nb and Ta (lovchorrite, rinkite, loparite, pyrochlore). Thorium is almost always present in the minerals of REE, which serve as one of the sources of its production (for example, in monazite Ce[PO4]). Therefore, Th closely correlates (r = 0.94) with Nb and Ta in the sediments of Lake Bolshoy Vudjavr, as well as with Ti (r = 0.87). All these four elements have similar distribution pattern (Figure 12).
Bismuth and tellurium show a very similar distribution pattern in the sediments of Lake Bolshoy Vudjavr, resembling that of other chalcophile elements (Figure 11), and their sources are probably emissions from fossil fuel combustion and copper–nickel production. The latter conclusion is supported by a more than 17-fold increase in the Bi content in the sediments of Lake Rybachy and a 7-fold increase in the sediments of Lake Severnoe compared to the background concentrations (Table 2) [27,83]. Since Bi and Te are not included in any of the described minerals of the Khibiny massif [52], these elements are most unlikely to have come into the lake with the effluents of the mines and the processing plant.

3.9. Trace Element Fractions in Sediments

Analysis of the forms of elements in the sediments of Lake Bolshoy Vudjavr showed that the studied elements are mainly found in stable compounds—mineral, acid-soluble, and organic. For example, Ni, Cu, Zn (Figure 13), Cd, W, V, Pb, Sb, and Sn are predominantly found in mineral phases throughout the entire depth of the studied sediment core, that is, in the most forms difficult to be involved in biogeochemical cycles. At the same time, the predominant phase for chalcophile elements Ni, Cu, Zn, Cd, Pb, W, and V, in addition to mineral forms, is the fraction associated with organic matter. Copper exhibits the most pronounced organophilic properties throughout the sediment core of Lake Bolshoy Vudjavr (Figure 13).
Among all the metals studied in this work, the largest relative amount of mobile and most potentially bioavailable (the sum of one and two fractions) forms over the entire depth of the core was found for Cd. In contrast to Cd, Sn is the most strongly bound metal, for which the absolute predominance of mineral forms was revealed throughout the core. For individual metals, Zn and Pb among them, the bond with oxides and hydroxides of Fe and Mn is essential throughout the entire depth of the investigated lake sediment core. In many respects, the established patterns are close to those obtained for the same elements in the sediments of Lake Komsomolskoye (Monchegorsk, Murmansk region) and in the sediments of small lakes in Petrozavodsk and Murmansk [71,77,103].

3.10. Anthropogenic Load on the Lake

The contamination of the sediments of Lake Bolshoy Vudjavr was assessed by the integral indicator PLI (Pollution Load Index) and the contamination factor (CF), where all the analyzed elements were taken into account, and the average concentrations of elements in the sediment layers of 24–27 cm, accumulated before the development of apatite–nepheline deposits started, were used as the background level. The value of the PLI increases from the lower layers to the upper ones (Figure 14). The contamination of the sediment layers of Lake Bolshoy Vudjavr, formed before the onset of active anthropogenic impact on the reservoir (20–27 cm), can be assessed as low. After the start of the activity of the “Apatite” Plant and the release of a large amount of runoff from the mines into the lake, a sharp increase in the PLI occurred and the sediment layers accumulated from the 1930s to the mid-1950s (14–20 cm),characterized by a moderate level of pollution. From the mid-1950s to the late 1980s, there was a gradual increase in the PLI values associated with the commissioning of new mines and an increase in the capacity for mining and processing apatite–nepheline ore, as well as a sharp increase in using automobile and fossil fuels (coal and masut), containing a large amount of heavy metals (Pb, Ni, V, W, see Figure 10). The highest PLI values (on the border between moderate and high pollution) were detected in the layers from 2 to 7 cm, which were accumulated from the early 1990s to the mid-2000s. The maximum concentrations of most of the elements present in the apatite–nepheline ore, including P, Ti, Nb, Ta, REE, alkali and alkaline earth metals, as well as heavy metals, are found precisely in these layers of the studied sediments (Figure 5, Figure 6, Figure 7, Figure 11 and Figure 12). Since the beginning of the 2000s (from 4 cm and above), the level of pollution has been gradually decreasing in the most recent sediment layers of Lake Bolshoy Vudjavr, which correlates well with a decrease in the content of individual elements, such as REE and heavy metals, in the surface sediment layers. Nevertheless, in general, the PLI values in the uppermost sediment layers remain on the border between moderate and high pollution (Figure 14).
Antimony (18.2), P (10.3), and Sr (7.8) have the maximum values of the contamination factor CF (high according to the classification [64]), as well as La (6.0) from the REE. The values of the contamination factor CF for the rest of the REE as well as Pb (4.6), Cu (4.0), U (3.5), and Bi (3.4) are also considerable.

4. Conclusions

Being the largest inland reservoir of the Khibiny alkaline mountain range, Lake Bolshoy Vudjavr has been experiencing intense anthropogenic pressure for more than 90 years, since the development of the richest apatite–nepheline deposits in the world started in 1929. A significant volume of raw materials containing ore minerals is currently stored in the tailings of apatite–nepheline factories. The development of apatite–nepheline deposits of the Khibiny alkaline massif led to the formation of runoff and emissions of the main rock-forming and trace elements into the environment and their subsequent sedimentation in the previously untouched mountain environment.
As a result of studying the sediments of Lake Bolshoy Vudjavr, elevated concentrations of alkali and alkaline earth metals, rare earth elements (REE), P, many trace elements, including heavy metals (Pb, Sn, Sb, Cu), associated with the extraction and processing of apatite–nepheline ores and the combustion of fossil fuels, were found in the surface sediment layers. Along with the macroelements extracted from apatite–nepheline ores at the processing plants at the present time, the increased content of Zr, Hf, and Mo was detected in the sediments of Lake Bolshoy Vudjavr in the 1930s, which is associated with the exploration and development of Zr raw materials and Mo ores at that time. The content of these elements decreased further to the surface, as the development of the raw materials was terminated. The content of Nb, Ta, Ti and Th increased sharply in the early 1930s due to the mining of sphene, loparite, and lovchorrite ores. Unlike the first three elements (Zr, Hf, and Mo), the content of the latter four in the lake sediments continues to increase to the surface layers, despite the cessation of mining of the loparite and lovchorrite ores.
Analysis of the forms of the studied elements in the sediments of Lake Bolshoy Vudjavr showed that they are mainly found in hard-to-reach (stable) fractions—mineral, acid-soluble, and organic, which indicates that they came into the lake sediments mainly as part of the runoff from apatite–nepheline mines and atmospheric fallout after fuel combustion. For example, Ni, Cu, Zn, Cd, W, V, Pb, Sb, and Sn are predominantly found in mineral phases throughout the entire studied sediment core, i.e., in the most difficult forms to be involved in biogeochemical cycles. At the same time, for Cd, Cu, Pb, Zn, Ni, W, and V, in addition to mineral forms, the predominant phase is the fraction associated with organic matter.Therefore, the latter elements can have a negative impact on hydrobionts, primarily on benthic organisms.
An assessment of the degree of pollution of the sediments of Lake Bolshoy Vudjavr using the integrated Pollution Load Index (PLI) revealed a high level of pollution. A sharp increase in the PLI value occurs first in the activity of the mining enterprise from the 1930s. This pollution reached its maximum in the late 1980s during the period of highest extraction of apatite–nepheline ores and has yet to decrease.
The diversity of the mineral composition of the apatite–nepheline deposits of the Khibiny and the mining activity predetermined the entry into Lake Bolshoy Vudjavr of a large list of elements that from the beginning of the deposits development to the present violate the natural geochemical cycles of elements and alter the chemical composition of the water column and lake sediments. As a result, the habitat of the flora and fauna of the mountain lake is changing, which, due to the natural environment, is vulnerable and slowly restoring natural resources.
Violation of the natural geochemical cycles of macroelements and rare elements, including alkaline, alkaline earth, and rare earth elements, causes serious changes in biological systems of various levels of their organization: from the organism to the ecosystem. The disintegration of rocks during technological processes increases the mobility of these elements, and increases their concentration in water and sediments. These elements are able to accumulate in the lake ecosystem, enter the food chain, and have a prolonged negative impact on living organisms, including the population health. Nevertheless, many elements accumulating in sediments (both macro- and microelements) are not given due attention in environmental practice. Many of these elements are not included in the lists of controlled elements; their behavior in natural environments and their impact on living organisms has not been studied. Therefore, a large number of uncontrolled elements can not only have a synergistic suppressive effect on lake hydrobionts, but also pose a serious threat to the health of the city’s population. Lake Bolshoy Vudjavr is a mining wastewater collector, and at the same time, the lake is an additional source of drinking water supply for Kirovsk. The problem of pollution of water supply sources in settlements is relevant for the population of the Murmansk region due to the presence of a large number of mining enterprises.
One of the solutions to the problem of intense pollution of Lake Bolshoy Vudjavr and other water bodies polluted by mining effluents can be the integrated development of the subsoil, which is one of the components of the theory of sustainable development. A large number of components that are valuable and necessary for the country’s economy can be extracted from the rocks of the apatite–nepheline deposits of the Khibiny mountain alkaline massif (raw materials for the extraction of Ti, Zr, Mo, Sr, F, Nb, Ta, Th, REE, etc., in addition to the currently extracted apatite and nepheline). The integrated use of mineral raw materials will reduce the load on the environment due to a decrease in the volume of industrial waste.

Author Contributions

Conceptualization, V.D., Z.S. and D.D.; methodology, Z.S., V.D. and A.G.; software, D.D. and Z.S.; validation, V.D., Z.S. and D.D.; investigation and resources, V.D., Z.S., D.D. and A.G.; data curation, V.D. and Z.S.; writing—original draft preparation, V.D. and Z.S.; writing-review and editing, V.D. and Z.S.; visualization, V.D., Z.S. and D.D., supervision, V.D.; project administration and funding acquisition, V.D. and Z.S. All authors have read and agreed to the published version of the manuscript.

Funding

The research is supported by the Russian Science Foundation (project No. 22-27-00131).

Acknowledgments

The authors sincerely thank the colleagues from Petrozavodsk and Apatity: P.M. Terentjev for his help in lake sediment sampling, as well as to our colleagues O.P. Korytnaya, A.S. Paramonov, S.V. Burdyukh, M.V. Ekhova, V.L. Utitsina, and A. Fedorov for conducting high-quality analytical studies.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Delile, H.; Blichert-Toft, J.; Goiran, J.-P.; Keay, S.; Albarède, F. Lead in ancient Rome’s city waters. Proceed. NAS USA 2014, 111, 6594–6599. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Renberg, I.; Wik-Persson, M.; Emteryd, O. Pre-industrial atmospheric lead contamination detected in Swedish lake sediments. Nature 1994, 368, 323–326. [Google Scholar] [CrossRef]
  3. Renberg, I.; Bindler, R.; Brännvall, M.L. Using the historical atmospheric lead-deposition record as a chronological marker in sediment deposits in Europe. Holocene 2001, 11, 511–516. [Google Scholar] [CrossRef]
  4. Carvalho, F.; Schulte, L. Reconstruction of mining activities in the Western Alps during the past 2500 years from natural archives. Sci. Total Environ. 2021, 750, 141208. [Google Scholar] [CrossRef]
  5. Killick, D.; Fenn, T. Archaeometallurgy: The study of preindustrial mining and metallurgy. Annu. Rev. Anthropol. 2012, 41, 559–575. [Google Scholar] [CrossRef]
  6. Martínez Cortizas, A.; López-Merino, L.; Bindler, R.; Mighall, T.; Kylander, M. Early atmospheric metal pollution provides evidence for Chalcolithic/Bronze Age mining and metallurgy in Southwestern Europe. Sci. Total Environ. 2016, 545–546, 398–406. [Google Scholar] [CrossRef] [Green Version]
  7. McConnell, J.R.; Wilson, A.I.; Stohl, A.; Arienzo, M.M.; Chellman, N.J.; Eckhardt, S.; Thompson, E.M.; Pollard, A.M.; Steffensen, J.P. Lead pollution recorded in Greenland ice indicates European emissions tracked plagues, wars, and imperial expansion during antiquity. Proc. Natl. Acad. Sci. USA 2018, 115, 5726–5731. [Google Scholar] [CrossRef] [Green Version]
  8. Durali-Mueller, S.; Brey, G.P.; Wigg-Wolf, D.; Lahaye, Y. Roman lead mining in Germany: Its origin and development through time deduced from lead isotope provenance studies. J. Archaeol. Sci. 2007, 34, 1555–1567. [Google Scholar] [CrossRef]
  9. Edmondson, J.C. Mining in the later Roman Empire and beyond: Continuity or disruption? J. Rom. Stud. 1989, 79, 84–102. [Google Scholar] [CrossRef]
  10. Martínez Cortizas, A.; López-Merino, L.; Bindler, R.; Mighall, T.; Kylander, M. Atmospheric Pb pollution in N Iberia during the late Iron Age/Roman times reconstructed using the high-resolution record of La Molina mire (Asturias, Spain). J. Paleolimnol. 2013, 50, 71–86. [Google Scholar] [CrossRef]
  11. Preunkert, S.; McConnell, J.R.; Hoffmann, H.; Legrand, M.; Wilson, A.I.; Eckhardt, S.; Stohl, A.; Chellman, N.J.; Arienzo, M.N.; Friedrich, R. Lead and antimony in basal ice from Col du Dome (French Alps) dated with radiocarbon: A record of pollution during antiquity. Geophys. Res. Lett. 2019, 46, 4953–4961. [Google Scholar] [CrossRef] [Green Version]
  12. Shotyk, W.; Krachler, M. Atmospheric deposition of silver and thallium since 12370 14C years BP recorded by a Swiss peat bog profile, and comparison with lead and cadmium. J. Environ. Monit. 2004, 6, 427–433. [Google Scholar] [CrossRef] [PubMed]
  13. Weiss, D.; Shotyk, W.; Appleby, P.G.; Cheburkin, A.K.; Kramers, J.D. Atmospheric Pb deposition since the Industrial Revolution recorded by five Swiss peat profiles: Enrichment factors, fluxes, isotopic composition, and sources. Environ. Sci. Technol. 1999, 33, 1340–1352. [Google Scholar] [CrossRef]
  14. Barabanov, A.V.; Kalinina, T.A.; Kiselev, A.A.; Krasnobaev, A.I. Giant in the Khibiny; Ore and Metals: Moscow, Russia, 1999; 287p. (In Russian) [Google Scholar]
  15. Kashulin, N.A.; Denisov, D.B.; Sandimirov, S.S.; Dauvalter, V.A.; Kashulina, T.G.; Malinovsky, D.N.; Vandysh, O.I.; Ilyashuk, B.P.; Kudryavtseva, L.P. Anthropogenic Changes in the Water Systems of the Khibiny Mountains (Murmansk Region); KSC RAS: Apatity, Russia, 2008; 282p. (In Russian) [Google Scholar]
  16. Dauvalter, V.; Moiseenko, T.; Kudryavtseva, L.; Sandimirov, S. Accumulation of heavy metals in Lake Imandra because of its pollution with industrial waste. Water Resour. 2000, 27, 279–287. [Google Scholar]
  17. Yuoguy, V.S.; Dauvalter, V.A. Characteristic of bottom sediments of Lakes Bolshoy and Malyi Vudjavr. In The Quaternary in All of Its Variety. Basic Issues, Results, and Major Trends of Further Research, Proceedings of the VII All-Russian Quaternary Conference, Apatity, Russia, September 12–17, 2011; Korsakova, O.P., Kolka, V.V., Eds.; RAS: Saint-Petersburg, Russia, 2011; Volume 2, pp. 333–336. (In Russian) [Google Scholar]
  18. Denisov, D.B.; Dauvalter, V.A.; Kashulin, N.A.; Kagan, L.Y. Long-term Changes of Diatom Assemblages in Subarctic Lakes under the Anthropogenic Influence (According to the Data of Diatom Analysis). Inland Water Biol. 2006, 1, 24–30. [Google Scholar]
  19. Denisov, D.B. Changes in the Hydrochemical Composition and Diatomic Flora of Bottom Sediments in the Zone of Influence of Metal Mining Production (Kola Peninsula). Water Resour. 2007, 34, 682–692. [Google Scholar] [CrossRef]
  20. Adriano, D. Trace Elements in the Terrestrial Environment; Springer-Verlag: New York, NY, USA, 2001; pp. 18–20. [Google Scholar]
  21. Kashulin, N.A.; Dauvalter, V.A.; Denisov, D.B.; Valkova, S.A.; Vandysh, O.I.; Terentjev, P.M.; Kashulin, A.N. Selected aspects of the current state of freshwater resources in the Murmansk Region, Russia. J. Environ. Sci. Health Part A 2017, 52, 921–929. [Google Scholar] [CrossRef]
  22. Schneider, L.; Mariani, M.; Saunders, K.M.; Maher, W.A.; Harrison, J.J.; Fletcher, M.-S.; Zawadzki, A.; Heijnis, H.; Haberle, S.G. How significant is atmospheric metal contamination from mining activity adjacent to the Tasmanian Wilderness World Heritage Area? A spatial analysis of metal concentrations using air trajectories models. Sci. Total Environ. 2019, 656, 250–260. [Google Scholar] [CrossRef]
  23. Denisov, D.B.; Kashulin, N.A.; Terentjev, P.M.; Valkova, S.A. Modern trends in the biota of freshwater ecosystems in the Murmansk Region. Proceed. Murm. St. Techn. Univer. 2009, 12, 525–538. (In Russian) [Google Scholar]
  24. Terentjev, P.M.; Kashulin, N.A.; Zubova, E.M. The role of the European smelt Osmeruseperlanus (Linnaeus) in the structure of the ichthyofauna of the lake basin Imandra (Murmansk region). Proceed. Zoolog. Inst. 2017, 321, 228–243. (In Russian) [Google Scholar]
  25. Moiseenko, T.I.; Gashkina, N.A.; Dinu, M.I.; Kremleva, T.A.; Khoroshavin, V.Y. Water Chemistry of Arctic Lakes under Airborne Contamination of Watersheds. Water 2020, 12, 1659. [Google Scholar] [CrossRef]
  26. Norton, S.A.; Appleby, P.G.; Dauvalter, V.; Traaen, T.S. Trace Metal Pollution in Eastern Finnmark, Norway and Kola Peninsula, Northeastern Russia as Evidences by Studies of Lake Sediment; NIVA-Report 41/1996; NIVA: Oslo, Norway, 1996; p. 18. [Google Scholar]
  27. Dauvalter, V.A.; Terentiev, P.M.; Denisov, D.B.; Udachin, V.N.; Filippova, K.A.; Borisov, A.P. Reconstruction of pollution of the territory of the Rybachy Peninsula of the Murmansk region with heavy metals. Proceed. Fersman Sci. Sess. Geolog. Inst. KSC RAS 2018, 15, 441–444. [Google Scholar]
  28. Köck, G.; Muir, D.; Yang, F.; Wang, X.; Talbot, C.; Gantner, N.; Moser, D. Bathymetry and Sediment Geochemistry of Lake Hazen (Quttinirpaaq National Park, Ellesmere Island, Nunavut). Arctic 2012, 65, 56–66. [Google Scholar] [CrossRef]
  29. Catalan, J.; Curtis, C.J.; Kernan, M. Remote European mountain lake ecosystems: Regionalisation and ecological status. Freshw. Biol. 2009, 54, 2419–2432. [Google Scholar] [CrossRef]
  30. Solovieva, N.; Jones, V.J.; Nazarova, L.; Brooks, S.J.; Birks, H.J.B.; Grytnes, J.-A.; Appleby, P.G.; Kauppila, T.; Kondratenok, B.; Renberg, I.; et al. Palaeolimnological evidence for recent climatic change in lakes from the northern Urals, arctic Russia. J. Paleolimn. 2005, 33, 463–482. [Google Scholar] [CrossRef]
  31. Solovieva, N.; Tarasov, P.E.; MacDonald, G. Quantitative reconstruction of Holocene climate from the Chuna Lake pollen record, Kola Peninsula, northwest Russia. Holocene 2005, 15, 165–172. [Google Scholar] [CrossRef]
  32. Beck, K.K.; Mariani, M.; Fletcher, M.-S.; Schneider, L.; Aquino-Lopez, M.A.; Gadd, P.S.; Heijnis, H.; Saunders, K.M.; Zawadzki, A. The impacts of intensive mining on terrestrial and aquatic ecosystems: A case of sediment pollution and calcium decline in cool temperate Tasmania, Australia. Environ. Pollut. 2020, 265, 114695. [Google Scholar] [CrossRef]
  33. Blake, J.M.; Brown, J.E.; Ferguson, C.L.; Bixby, R.J.; Delay, N.T. Sediment record of mining legacy and water quality from a drinking-water reservoir, Aztec, New Mexico, USA. Environ. Earth Sci. 2020, 79, 404. [Google Scholar] [CrossRef]
  34. Dauvalter, V.A.; Denisov, D.B.; Slukovskii, Z.I. Impact of Wastewaters from Apatite–Nepheline Production on the Biogeochemical Processes in an Arctic Mountain Lake. Geochem. Int. 2022, 60, 1014–1028. [Google Scholar] [CrossRef]
  35. Fuge, R.; Paveley, C.F.; Holdham, M.T. Heavy metal contamination in the Tanat Valley, North Wales. Environ. Geochem. Health 1989, 11, 127–135. [Google Scholar] [CrossRef]
  36. Horowitz, A.J.; Elrick, K.A.; Robbins, J.A.; Cook, R.B. A summary of the effects of mining and related activities on the sediment-trace element geochemistry of Lake Coeur d’Alene, Idaho, USA. J. Geochem. Explor. 1995, 52, 135–144. [Google Scholar] [CrossRef]
  37. Horowitz, A.J.; Elrick, K.A.; Robbins, J.A.; Cook, R.B. Effect of mining and related activities on the sediment trace element geochemistry of Lake Coeur d’Alene, Idaho, USA part II: Subsurface sediments. Hydrol. Process. 1995, 9, 35–54. [Google Scholar] [CrossRef]
  38. Jung, M. Environmental Contamination of Heavy Metals in Soils, Plants, Waters and Sediments in the Vicinity of Metalliferous Mine in Korea. Unpublished. Ph.D Thesis, University of London, London, UK, 1995; pp. 189–317. [Google Scholar]
  39. Jung, M.C. Heavy metal contamination of soils and waters in and around the Imcheon Au–Ag mine, Korea. Appl. Geochem. 2001, 16, 1369–1375. [Google Scholar] [CrossRef]
  40. Jung, M.C.; Thornton, I. Heavy metal contamination of soils and plants in the vicinity of a lead-zinc mine, Korea. Appl. Geochem. 1996, 11, 53–59. [Google Scholar] [CrossRef]
  41. Jung, M.C.; Thornton, I. Environmental contamination and seasonal variation of metals in soils, plants and waters in the paddy fields around a Pb–Zn mine in Korea. Sci. Total Environ. 1997, 198, 105–121. [Google Scholar] [CrossRef]
  42. Langman, J.B.; Ali, J.D.; Child, A.W.; Wilhelm, F.M.; Moberly, J.G. Sulfur Species, Bonding Environment, and Metal Mobilization in Mining-Impacted Lake Sediments: Column Experiments Replicating Seasonal Anoxia and Deposition of Algal Detritus. Minerals 2020, 10, 849. [Google Scholar] [CrossRef]
  43. Lee, C.G.; Chon, H.-T.; Jung, M.C. Heavy metal contamination in the vicinity of the Daduk Au–Ag–Pb–Zn mine in Korea. Appl. Geochem. 2001, 16, 1377–1386. [Google Scholar] [CrossRef]
  44. Merrington, G.; Alloway, B. The transfer and fate of Cd, Cu, Pb and Zn from two historic metalliferous mine sites in the UK. Appl. Geochem. 1994, 9, 677–687. [Google Scholar] [CrossRef]
  45. Pestana, M.H.D.; Formoso, M.L.; Teixeir, E.C. Heavy metals in stream sediments from copper and gold mining areas in southern Brazil. J. Geochem. Explor. 1997, 58, 133–143. [Google Scholar] [CrossRef]
  46. Thornton, I. Geochemical aspects of heavy metal pollution and agriculture in England and Wales. Ref. Book Minist. Agric. Fish. Food 1980, 326, 105–125. [Google Scholar]
  47. Wurtsbaugh, W.A.; Leavitt, P.R.; Moser, K.A. Effects of a century of mining and industrial production on metal contamination of a model saline ecosystem, Great Salt Lake, Utah. Environ. Pollut. 2020, 266, 115072. [Google Scholar] [CrossRef] [PubMed]
  48. Bouzekri, S.; El Fadili, H.; El Hachimi, M.L.; El Mahi, M.; El Lotfi, M. Assessment of trace metals contamination in sediment and surface water of quarry lakes from the abandoned Pb mine Zaida, High Moulouya-Morocco. Environ. Develop. Sustain. 2020, 22, 7013–7031. [Google Scholar] [CrossRef]
  49. Dauvalter, V.; Moiseenko, T.; Rodyushkin, I. Geochemistry of rare earth elements in Imandra Lake, Murmansk Area. Geochem. Int. 1999, 37, 325–331. [Google Scholar]
  50. Cooke, C.A.; Abbott, M.B. A paleolimnological perspective on industrial-era metal pollution in the central Andes, Peru. Sci. Tot. Environ. 2008, 393, 262–272. [Google Scholar] [CrossRef] [PubMed]
  51. Fersman, A.E. Our Apatite; Nauka: Moscow, Russia, 1968; 136p. (In Russian) [Google Scholar]
  52. Yakovenchuk, V.N.; Ivanyuk, G.Y.; Pakhomovsky, Y.A.; Men’shikov, Y.P. Khibiny; Laplandia Minerals: Apatity, Russia, 2005; p. 468. [Google Scholar]
  53. Arzamastsev, A.A.; Arzamastseva, L.V.; Shamatrina, A.M.; Travin, A.V.; Belyatsky, B.V.; Antonov, A.V.; Larionov, A.N.; Rodionov, N.V.; Sergeev, S.A. Duration of formation of magmatic system of polyphase paleozoic alkaline complexes of the Central Kola: U-Pb, Rb-Sr, Ar-Ar data. Dokl. Earth Sci. 2007, 413, 432–436. [Google Scholar] [CrossRef]
  54. Kola Ecological Center. To Prevent the Resumption of the Conflict in the Khibin. Available online: https://www.change.org/ (accessed on 1 December 2022).
  55. Kanygina, A.V. Hydrobiological and Hydrochemical Study of Bolshoy and Maly Vudjavr Lakes; Kola Research Base of Academy of Sciences of USSR: Apatity, Russia, 1939; 206p. (In Russian) [Google Scholar]
  56. Vinogradova, A.A.; Kotova, E.I.; Topchaya, V.Y. Atmospheric transport of heavy metals to the regions of the North of the European territory of Russia. Geogr. Natl. Resour. 2017, 1, 108–116. [Google Scholar] [CrossRef]
  57. Dauvalter, V.A. Geochemistry of Lakes in a Zone Impacted by an Arctic Iron-Producing Enterprise. Geochem. Int. 2020, 58, 933–946. [Google Scholar] [CrossRef]
  58. Gregurek, D.; Melcher, F.; Pavlov, V.A.; Reimann, C.; Stumpf, E.F. Mineralogy and mineral chemistry of snow filter residues in the vicinity of the nickel-copper processing industry, Kola Peninsula, NW Russia. Miner. Petrol. 1999, 65, 87–111. [Google Scholar] [CrossRef]
  59. Håkanson, L. Sediment sampling in different aquatic environments: Statistical aspects. Water Resour. Res. 1984, 20, 41–46. [Google Scholar] [CrossRef]
  60. Håkanson, L. Determination of characteristic values for physical and chemical lake sediment parameters. Water Resour. Res. 1981, 17, 1625–1640. [Google Scholar] [CrossRef]
  61. Taft, R.A.; Jones, C. Sediment Sampling Guide and Methodologies; Ohio Environmental Protection Agency: Columbus, OH, USA, 2001; p. 36. [Google Scholar]
  62. Dauvalter, V.A. Geoecology of Lake Sediments; Murmansk St. Techn. Univer.: Murmansk, Russia, 2012; 242p. (In Russian) [Google Scholar]
  63. Slukovskii, Z.I. Background concentrations of heavy metals and other chemical elements in the sediments of small lakes in the south of Karelia, Russia. Proceed. Murm. St. Techn. Univer. 2020, 23, 80–92. [Google Scholar] [CrossRef]
  64. Håkanson, L. An ecological risk index for aquatic pollution control—a sedimentological approach. Water Resear. 1980, 14, 975–1001. [Google Scholar] [CrossRef]
  65. Shukla, B.S.; Joshi, S.R. An evaluation of the CIC model of 210Pb dating of sediments. Environ. Geol. Water Sci. 1989, 14, 73–76. [Google Scholar] [CrossRef]
  66. Blanchard, R.L. Rapid Determination of Lead-210 and Polonium-210 in Environmental Samples by Deposition on Nickel. Anal. Chem. 1966, 38, 189–192. [Google Scholar] [CrossRef]
  67. Tessier, A.; Campbell, P.G.; Bisson, M. Sequential extraction procedure for the speciation of particulate trace metals. Anal. Chem. 1979, 51, 844–851. [Google Scholar] [CrossRef]
  68. Moore, F.; Nematollahi, M.J.; Keshavarzi, B. Heavy metals fractionation in surface sediments of Gowatr bay-Iran. Environ. Monit. Assess. 2015, 187, 4117. [Google Scholar] [CrossRef] [PubMed]
  69. McLaren, R.G.; Crawford, D.W. Studies on soil copper 1. the fractionation of copper in soils. J. Soil Sci. 1973, 24, 172–181. [Google Scholar] [CrossRef]
  70. Ure, A.M.; Quevauviller, P.; Muntau, H.; Griepink, B. Speciation of Heavy Metals in Soils and Sediments. An Account of the Improvement and Harmonization of Extraction Techniques Undertaken Under the Auspices of the BCR of the Commission of the European Communities. Int. J. Environ. Anal. Chem. 1993, 51, 135–151. [Google Scholar] [CrossRef]
  71. Guzeva, A.; Slukovskii, Z.; Dauvalter, V.; Denisov, D. Trace element fractions in sediments of urbanised lakes of the arctic zone of Russia. Environ. Monit. Assess. 2021, 193, 378. [Google Scholar] [CrossRef]
  72. Tomlinson, D.L.; Wilson, J.G.; Harris, C.R.; Jeffrey, D.W. Problems in the assessment of heavy-metal levels in estuaries and the formation of a pollution index. Helgoländer Meeresunters. 1980, 33, 566–575. [Google Scholar] [CrossRef] [Green Version]
  73. Subetto, D.A. Bottom Sediments of Lakes: Paleolimnological Reconstructions; Rus. St. Pedag. Univer.: Saint-Petersburg, Russia, 2009; 343p. (In Russian) [Google Scholar]
  74. Strakhovenko, V.D. Geochemistry of Sediments of Small Continental Lakes in Siberia. Dr. Sci. Thesis, Novosibirsk St. Univer., Novosibirsk, Russia, 2011; 24p. (In Russian). [Google Scholar]
  75. Dauvalter, V. Impact of mining and refining on the distribution and accumulation of nickel and other heavy metals in sediments of subarctic lake Kuetsjärvi, Murmansk region, Russia. J. Environ. Monit. 2003, 5, 210–215. [Google Scholar] [CrossRef] [PubMed]
  76. Rognerud, S.; Norton, S.A.; Dauvalter, V. Heavy Metal Pollution in Lake Sediments in the Border Areas between Russia and Norway; NIVA-Report 522/93; NIVA: Oslo, Norway, 1993; p. 18. [Google Scholar]
  77. Rognerud, S.; Hongve, D.; Fjeld, E.; Ottesen, R.T. Trace metal concentrations in lake and overbank sediments in southern Norway. Environ. Geol. 2000, 39, 723–732. [Google Scholar] [CrossRef]
  78. Slukovskii, Z.; Medvedev, M.; Siroezhko, E. Long-range transport of heavy metals as a factor of the formation of the geochemistry of sediments in the southwest of the Republic of Karelia, Russia. J. Elementol. 2020, 25, 125–137. [Google Scholar] [CrossRef]
  79. Demidov, I.N.; Shelekhova, T.S. Diatomites of Karelia (Features of Formation, Distribution, Prospects of Use); Karelian Research Centre of RAS: Petrozavodsk, Russia, 2006; 89p. (In Russian) [Google Scholar]
  80. Kashulin, N.A.; Sandimirov, S.S.; Dauvalter, V.A.; Kudryavtseva, L.P.; Terentyev, P.M.; Denisov, D.B.; Vandysh, O.I.; Korolova, I.M.; Valkova, S.A.; Kashulina, T.G. Annotated Ecological Catalogue of Lakes of the Murmansk Region: The Central and Southwest Areas of the Murmansk Region (Basins of the Barents Sea, the White Sea and the Bothnia Gulf of the Baltic Sea); KSC RAS: Apatity, Russia, 2013; 298p. (In Russian) [Google Scholar]
  81. Moiseenko, T.; Dauvalter, V.; Rodyushkin, I. Geochemical Migration and Covariation of Elements in the Imandra Lake, Barents Region; Research Report TULEA; Univer. Techn.: Luleå, Sweden, 1996; Volume 22, 108p. [Google Scholar]
  82. Moiseenko, T.I.; Dauvalter, V.A.; Rodyushkin, I.V. Geochemical Migration of Elements in Subarctic Freshwater Body (on the Example of the Imandra Lake); KSC RAS: Apatity, Russia, 1997; 127p. (In Russian) [Google Scholar]
  83. Slukovskii, Z.I.; Guzeva, A.V.; Dauvalter, V.A. Rare earth elements in surface lake sediments of Russian arctic: Natural and potential anthropogenic impact to their accumulation. Appl. Geochem. 2022, 142, 105325. [Google Scholar] [CrossRef]
  84. Balaram, V. Rare earth elements: A review of applications, occurrence, exploration, analysis, recycling, and environmental impact. Geosci. Front. 2019, 10, 1285–1303. [Google Scholar] [CrossRef]
  85. Kaulina, T.V.; Afanasieva, E.N.; Il’chenko, V.L.; Nerovich, L.I.; Lialina, L.M.; Nitkina, E.A.; Mokrushina, O.D. Litsa Uranium-Ore District; GEOS: Moscow, Russia, 2021; 136p. (In Russian) [Google Scholar]
  86. Krasotkin, I.S.; Voytekhovsky, Y.L.; Leskov, A.L.; Udobina, V.S. Abandoned molybdenite mine Takhtarvumchorr. Proceed. Fersman Sci. Sess. Geolog. Inst. KSC RAS 2005, 2, 10–14. [Google Scholar]
  87. Slukovskii, Z.I.; Dauvalter, V.A. Features of Pb, Sb, Cd accumulation in sediments of small lakes in the south of the Republic of Karelia. Proceed. Karelian SC RAS 2020, 4, 75–94. (In Russian) [Google Scholar] [CrossRef]
  88. Slukovskii, Z.; Dauvalter, V.; Guzeva, A.; Denisov, D.; Cherepanov, A.; Siroezhko, E. The Hydrochemistry and Recent Sediment Geochemistry of Small Lakes of Murmansk, Arctic Zone of Russia. Water 2020, 12, 1130. [Google Scholar] [CrossRef]
  89. Krachler, M.; Zheng, J.; Koerner, R.; Zdanowicz, C.; Fisher, D.; Shotyk, W. Increasing atmospheric antimony contamination in the northern hemisphere: Snow and ice evidence from Devon Island, Arctic Canada. J. Environ. Monit. 2006, 7, 1169–1176. [Google Scholar] [CrossRef]
  90. Yudovich, Y.E.; Ketris, M.P. Toxic Trace Elements in Fossil Coals; Ural Branch RAS: Yekaterinburg, Russia, 2005; 648p. (In Russian) [Google Scholar]
  91. Teng, Y.; Ni, S.; Zhang, C.; Wang, J.; Lin, X.; Huang, Y. Environmental geochemistry and ecological risk of vanadium pollution in Panzhihua mining and smelting area, Sichuan, China. Chin. J. Geochem. 2006, 25, 379–385. [Google Scholar] [CrossRef]
  92. Mejia, J.A.; Rodriguez, R.; Armienta, A. Aquifer Vulnerability Zoning, an Indicator of Atmospheric Pollutants Input? Vanadium in the Salamanca Aquifer, Mexico. Water Air Soil Pollut. 2007, 185, 95–100. [Google Scholar] [CrossRef]
  93. Slukovskii, Z.I. Accumulation level and fractions of heavy metals in sediments of small lakes of the urbanized area (Karelia). Vestn. St. Petersburg Univer. Earth Sci. 2020, 65, 171–192. (In Russian) [Google Scholar] [CrossRef]
  94. Shpirt, M.Y.; Rashevsky, V.V. Trace Elements of Fossil Fuels; Kuchkovo Pole: Moscow, Russia, 2010; 348p. (In Russian) [Google Scholar]
  95. Denisov, D.; Terentjev, P.; Valkova, S.; Kudryavtzeva, L. Small Lakes Ecosystems under the Impact of Non-Ferrous Metallurgy (Russia, Murmansk Region). Environments 2020, 7, 29. [Google Scholar] [CrossRef] [Green Version]
  96. Vinogradova, A.A.; Ivanova, Y.A. Pollution of central Karelia environment under long-range atmospheric transport of anthropogenic pollutants. Izv. RAS Geogr. Ser. 2013, 5, 98–108. (In Russian) [Google Scholar]
  97. Pacyna, J.M.; Pacyna, E.G. An assessment of global and regional emissions of trace metals to the atmosphere from anthropogenic sources worldwide. Environ. Rev. 2001, 9, 269–298. [Google Scholar] [CrossRef]
  98. Dauvalter, V.A. Chalcophile elements (Hg, Cd, Pb, and As) in bottom sediments of water bodies of the White Sea catchment area on the Kola Peninsula. Geochem. Int. 2006, 44, 205–208. [Google Scholar] [CrossRef]
  99. Dauvalter, V.A.; Kashulin, N.A. Biogeochemical Features of the Distribution of Chalcophilic Elements (Hg, Cd, Pb, As) in Water Bodies of the North of the European Part of Russia; Murmansk St. Techn. Univ.: Murmansk, Russia, 2015; 136p. (In Russian) [Google Scholar]
  100. Starodymova, D.P.; Shevchenko, V.P.; Kokryatskaya, N.M.; Aliev, R.A.; Bychkov, A.Y.; Zabelina, S.A.; Chupakov, A.V. Geochemistry of bottom sediments of a small lake (drainage area of Lake Onega, Arkhangelsk region). Succes. Mod. Natl. Sci. 2016, 9, 172–177. (In Russian) [Google Scholar]
  101. Thomas, V. The elimination of lead in gasoline. Ann. Rev. Energy Environ. 1995, 20, 301–324. [Google Scholar] [CrossRef]
  102. Dauvalter, V.A.; Kashulin, N.A. Mercury pollution of Lake Imandra sediments, Murmansk Region, Russia. Int. J. Environ. Res. 2018, 12, 939–953. [Google Scholar] [CrossRef]
  103. Slukovskii, Z.I.; Dauvalter, V.A.; Denisov, D.B.; Siroezhko, E.V.; Cherepanov, A.A. Geochemistry features of sediments of small urban arctic Lake Komsomolskoye, Murmansk region. IOP Confer. Series: Earth Environ. Sci. 2020, 467, 012004. [Google Scholar] [CrossRef]
Figure 1. Map of the study area.
Figure 1. Map of the study area.
Water 14 04044 g001
Figure 2. Contents of 210Pb and 210Pbex. (Bq/kg) in separate layers of the sediment core from Lake Bolshoy Vudjavr with indication when the individual layers were formed, calculated from the average sedimentation rate.
Figure 2. Contents of 210Pb and 210Pbex. (Bq/kg) in separate layers of the sediment core from Lake Bolshoy Vudjavr with indication when the individual layers were formed, calculated from the average sedimentation rate.
Water 14 04044 g002
Figure 3. The average content (%) of the main elements and organic matter (LOI) in the modern sediments of Lake Bolshoy Vudjavr.
Figure 3. The average content (%) of the main elements and organic matter (LOI) in the modern sediments of Lake Bolshoy Vudjavr.
Water 14 04044 g003
Figure 4. Vertical distribution of the content (%) of the main elements and organic matter (LOI) in the sediment core of Lake Bolshoy Vudjavr.
Figure 4. Vertical distribution of the content (%) of the main elements and organic matter (LOI) in the sediment core of Lake Bolshoy Vudjavr.
Water 14 04044 g004
Figure 5. Vertical distribution of phosphorus content (%) in the sediments of Lake Bolshoy Vudjavr.
Figure 5. Vertical distribution of phosphorus content (%) in the sediments of Lake Bolshoy Vudjavr.
Water 14 04044 g005
Figure 6. Vertical distribution of the content (µg/g) of alkali and alkaline earth metals in the sediments of Lake Bolshoy Vudjavr.
Figure 6. Vertical distribution of the content (µg/g) of alkali and alkaline earth metals in the sediments of Lake Bolshoy Vudjavr.
Water 14 04044 g006
Figure 7. Vertical distribution of the content (µg/g) of light rare earth elements in the sediments of Lake Bolshoy Vudjavr.
Figure 7. Vertical distribution of the content (µg/g) of light rare earth elements in the sediments of Lake Bolshoy Vudjavr.
Water 14 04044 g007
Figure 8. Vertical distribution of the content (µg/g) of Zr and Hf in the sediments of Lake Bolshoy Vudjavr.
Figure 8. Vertical distribution of the content (µg/g) of Zr and Hf in the sediments of Lake Bolshoy Vudjavr.
Water 14 04044 g008
Figure 9. Vertical distribution of the content (µg/g) of Mo and As in the sediments of Lake Bolshoy Vudjavr.
Figure 9. Vertical distribution of the content (µg/g) of Mo and As in the sediments of Lake Bolshoy Vudjavr.
Water 14 04044 g009
Figure 10. Vertical distribution of the content of elements (µg/g) that reduce the content in the upper sediment layers of Lake Bolshoy Vudjavr.
Figure 10. Vertical distribution of the content of elements (µg/g) that reduce the content in the upper sediment layers of Lake Bolshoy Vudjavr.
Water 14 04044 g010
Figure 11. Vertical distribution of the content (µg/g) of Cu, Zn, Te, and Bi in the sediments of Lake Bolshoy Vudjavr.
Figure 11. Vertical distribution of the content (µg/g) of Cu, Zn, Te, and Bi in the sediments of Lake Bolshoy Vudjavr.
Water 14 04044 g011
Figure 12. Vertical distribution of the content (µg/g) of Nb, Ta, Ti, and Th in the sediments of Lake Bolshoy Vudjavr.
Figure 12. Vertical distribution of the content (µg/g) of Nb, Ta, Ti, and Th in the sediments of Lake Bolshoy Vudjavr.
Water 14 04044 g012
Figure 13. Concentrations (µg/g) of various forms of Ni, Cu, and Zn in the sediments (layers 0–5, 5–10, 10–15, and 18–23 cm) of Lake Bolshoy Vudjavr, forms: I—water-soluble forms, II—exchangeable cations (mobile forms), III—forms associated with Fe and Mn hydroxides, IV—forms associated with organic matter, V—acid-soluble forms, and VI—mineral forms.
Figure 13. Concentrations (µg/g) of various forms of Ni, Cu, and Zn in the sediments (layers 0–5, 5–10, 10–15, and 18–23 cm) of Lake Bolshoy Vudjavr, forms: I—water-soluble forms, II—exchangeable cations (mobile forms), III—forms associated with Fe and Mn hydroxides, IV—forms associated with organic matter, V—acid-soluble forms, and VI—mineral forms.
Water 14 04044 g013
Figure 14. Distribution of Pollution Load Index (PLI) value in the sediment core of Lake Bolshoy Vudjavr.
Figure 14. Distribution of Pollution Load Index (PLI) value in the sediment core of Lake Bolshoy Vudjavr.
Water 14 04044 g014
Table 1. The content of organic matter (LOI) and macroelements (in %) in apatite ore, tailings of apatite production, and surface and background layers of sediments in lakes Bolshoy Vudjavr and Imandra according to the results of X-ray fluorescence analysis.
Table 1. The content of organic matter (LOI) and macroelements (in %) in apatite ore, tailings of apatite production, and surface and background layers of sediments in lakes Bolshoy Vudjavr and Imandra according to the results of X-ray fluorescence analysis.
ElementsApatite Ore1Tailings1Bolshoy VudjavrImandra1
0–1 cm26–27 cm0–1 cm18–19 cm
LOI14.0818.3629.3313.69
SiO226.639.837.6851.0969.4375.27
Al2O312.117.715.8716.359.519.98
Fe2O36.498.145.744.456.413.99
Na2O6.239.743.481.950.951.77
CaO249.729.682.772.982.31
K2O2.834.283.651.671.011.07
MgO1.271.521.411.820.801.72
MnO0.1880.2470.190.120.920.14
TiO23.014.071.330.620.460.41
Notes: 1 [81,82].
Table 2. The content of trace elements (in µg/g) in apatite ore, tailings of apatite production, and surface and background layers of sediments of lakes Bolshoy Vudjavr, Imandra, Rybachy, and Severnoye according to the results of mass spectrometry with inductively coupled plasma (ICP-MS).
Table 2. The content of trace elements (in µg/g) in apatite ore, tailings of apatite production, and surface and background layers of sediments of lakes Bolshoy Vudjavr, Imandra, Rybachy, and Severnoye according to the results of mass spectrometry with inductively coupled plasma (ICP-MS).
ElementApatite Ore 1Tailings1Bolshoy VudjavrImandra1Rybachy2Severnoye3
0–1 cm26–27 cm0–1 cm18–19 cm0–1 cm26–27 cm0–1 cm37–39 cm
Li4.416.4519.4618.3415.113.413.28.7012.635.07
Be3.216.395.483.671.840.751.301.490.640.49
P737591609333590321334881128
Sc2.520.769.1517.052.627.964.915.2212.189.76
Ti719101191093932173313621257479
V26130313711994.510963.153.552840
Cr32.833.747.3187.7994.911261.551.158.4934.77
Mn23440016249641932203174114
Co19.617.619.1515.4329.616.340.03.9820.0813.85
Ni13.49.3832381751591272417127
Cu79.110122555429963703811125
Zn4771245126132781485641664
Ga26.439.927.5764.4427.730.87.525.90
As2.800.863.786.1119.19.16.281.89
Rb689022010449.837.432.418.719.283.83
Sr139003710518665294118715316915753
Y16172.910430402117361011
Zr31758460664516312353452611
Nb168288191108368.55.43.73.01.2
Mo66.538102411.46.111.48.323.5
Cd0.360.0810.710.571.290.340.500.451.180.30
Sn<24<243.902.1716.3013.5032.10.617.110.80
Sb85704.820.290.500.1382.680.09
Cs6.842.551.980.910.830.25
Ba6031100125760668548619816925878
La1285280535891365729702620
Ce21304607791532368150983927
Pr1904479.215.424.813.66.615.65.14.9
Nd64015428354955325621919
Sm81.321.643.69.013.67.75.1711.633.233.31
Eu24.67.112.212.383.911.831.082.410.750.71
Gd98.126.834.87.014.107.253.848.922.732.78
Tb8.162.514.611.071.460.860.561.310.360.38
Dy33.111.422.175.698.684.112.986.781.741.92
Ho5.331.963.700.991.430.690.531.110.310.36
Er11.404.649.092.963.451.961.733.360.871.04
Tm1.220.5451.060.390.320.300.200.380.120.13
Yb6.042.795.552.442.571.761.342.620.690.88
Lu0.610.3090.610.300.270.240.200.350.100.12
Hf0.730.3289.3511.084.364.101.691.340.590.31
Ta0.610.1737.343.822.120.640.190.090.230.06
W<2434.54.002.8238.1031.40.780.338.701.00
Hg0.0510.0240.3080.0500.2820.034
Tl0.120.190.140.090.490.29
Pb2.662.9924.69.928.04.582.010.744.06.3
Bi0.350.100.980.0580.240.07
Th13.714.923.811.38.07.94.76.42.31.2
U2.763.2512.89.48.015.14419077204
Notes: 1 [81,82], 2 [27], 3 [83].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dauvalter, V.; Slukovskii, Z.; Denisov, D.; Guzeva, A. A Paleolimnological Perspective on Arctic Mountain Lake Pollution. Water 2022, 14, 4044. https://doi.org/10.3390/w14244044

AMA Style

Dauvalter V, Slukovskii Z, Denisov D, Guzeva A. A Paleolimnological Perspective on Arctic Mountain Lake Pollution. Water. 2022; 14(24):4044. https://doi.org/10.3390/w14244044

Chicago/Turabian Style

Dauvalter, Vladimir, Zakhar Slukovskii, Dmitry Denisov, and Alina Guzeva. 2022. "A Paleolimnological Perspective on Arctic Mountain Lake Pollution" Water 14, no. 24: 4044. https://doi.org/10.3390/w14244044

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop