
The paper is a logical continuation of the study of the chemical composition of natural waters in various regions of the Murmansk region and their impact on human health. The relevance of such studies is driven by the elevated incidence of diseases in the Murmansk region. Monitoring of groundwater in the Apatity-Kirovsky district has been carried out using inductively coupled plasma mass spectrometry (ICP-MS). The paper presents the results of a hygienic assessment of groundwater from the Khibiny Massif based on the concentrations of macro- and microelements. The aim of the study is to evaluate the chemical composition of groundwater in relation to the criteria for physiologically adequate drinking waters, selected from the underground water intake and after transportation through the water supply networks. It has been established that the drinking water supplied to the population of the Apatity-Kirovsky district, both in its natural state and at the point of delivery to consumers, does not have a balanced, physiologically complete composition in terms of the content of major macro- and microelements (magnesium, calcium, bicarbonates) and the value of total mineralization, and therefore cannot be recognized as physiologically complete. Calcium and magnesium concentrations are about three orders of magnitude lower than the recommended values for waters of different categories and by the criterion of completeness. Due to the low calcium concentration, there is a need to revise the regional maximum permissible concentrations (MPC) for strontium. The waters contain rare and rare-earth elements (Zr, La, Ce, Th, etc.) for which MPCs have not been developed, and the effects of these elements on the human organism are still not sufficiently studied.
The integration of modernized emanation surveys and seismic exploration is an effective tool for diagnosing the geodynamic activity of rock masses. In 2023, experimental and methodological integrated geophysical surveys were conducted at the Severopeschanskaya mine site (Sverdlovsk region, Krasnoturyinsk). The unique properties of radon and thoron determine their use as natural indicators of recent geodynamics. An improved emanation survey technique allows for the rapid (without lengthy monitoring) identification of the geodynamic component of radon field formation, zones of tensile and compressive deformation in rock masses, and the delineation and ranking of active deformation zones by hazard. Seismic exploration studies provide information on the distribution of elastic properties of rocks in the section and the location of structural heterogeneities. The results of comprehensive experimental and methodological geophysical studies demonstrate that radon anomalies correlate with deep structural heterogeneities identified by seismic exploration; the obtained data are confirmed by prior geological and geodetic information. The proposed integrated approach eliminates the limitations of radon measurements and significantly improves the information content of geodynamic diagnostics performed to ensure the safety of subsoil use facilities.
Polycyclic aromatic hydrocarbons have marker properties and serve as indicators of pollution sources associated with natural or man-made processes. During the extraction of PAHs from model mixtures using various solvents, it has been found that the use of n-hexane ensures high efficiency of compound extraction from surface water samples taken in areas with different anthropogenic loads (extract weight of 400–1900 µg/L). Using high-performance liquid chromatography, individual PAHs have been identified and their total content in the samples has been determined. Water from the area affected by the oil spill contains 0.25–0.37 µg/L PAHs, 78 % of which are low-molecular compounds (naphthalene, fluorene, phenanthrene). In waters from the area where oil production facilities are located, the total content is determined in the range of 0.40–1.10 µg/L; the main contribution (77 %) is made by high-molecular compounds [fluoranthene, pyrene, chrysene and benzo(a)pyrene]. Surface waters from coal mining areas contain PAHs in amounts of 0.24–0.88 µg/l; the presence of dibenz(a,h)anthracene (0.01–0.05 µg/l) has been noted in them. When calculating indicator ratios, it has been determined that the sources of PAH entry into surface waters from oil spill and coal mining areas are mainly natural, petrogenic and pyrogenic in nature. In areas with a load from the activities of oil producing enterprises, the entry of PAHs is caused by both petrogenic (fluid influx from strata) and anthropogenic (spill or combustion of oil products) processes.
The stability of an earth dam may be compromised by water erosion and seepage, which primarily affect the distribution of specific electrical resistivity within the dam's soil. This paper presents the results of complex geophysical electrical surveys of the dam body conducted to monitor potential water seepage through its clay screen. The study presents investigations using two electrical survey methods: the self-potential method and electrical tomography. The self-potential method is a passive electrical survey method that determines the potentials arising in the rock mass due to the transport of charges by moving water flows. Analysis of the self-potential field method data has identified a zone of high potential values – the length of this section exceeds 400 m, and the potential values exceed 60 mV. Electrical tomography allows the structure of the rock mass to be determined based on the distribution of the specific electrical resistivity of the rocks. The results of the electrical tomography have indicated that the dam body and the bedrock at its foundation are saturated with water almost everywhere. Processing and interpretation of the measurement data have made it possible to identify potential problem areas for water seepage.
The paper considers the possible correlation of the rocks of the Chervurt Formation in the Kolmozero-Voron'ya zone with its constituent supracrustal complexes and the formations of the adjacent Keivy domain, as well as the role that metasedimentary rocks of the upper part of this zone's section played in the formation history of the Keivy median massif. The study has employed the authors' original method for assessing statistically significant patterns in petrogeochemical characteristics, which has yielded good results in solving such problems, as previously demonstrated by the authors in a number of publications devoted to the correlation of supracrustal complexes of the Bolshiye Keivy and Malyye Keivy structures. The presented data allow us to conclude that the Chervurt rocks of the Kolmozero-Voron'ya zone are most similar to the metasedimentary complexes of the Keivy domain, and within them, to the rocks of the Upper Chervurt Subformation. This is also indicated by other clearly manifested patterns previously noted by the authors for the Upper Chervurt Subformation of the Keivy structure. It is also shown that the rocks of the Chervurt formation of the Kolmozero-Voron'ya belt are compositionally most distant from other complexes constituting this zone, primarily the Lyavozero complex, with which some researchers have attempted to correlate them. Taking into account these data, as well as information on the age of the rocks composing the Keivy domain and the Kolmozero-Voron'ya greenstone belt, and the authors' previously obtained data on the possible formation settings of certain Keivy formations, a model for the final stage of the formation of the Keivy median massif has been proposed.
Large-scale environmental contamination caused by the Baikal Pulp and Paper Plant and its industrial waste determines the relevance of studying the chemical composition of abiotic and biotic components in the impact area of the anthropogenic source. The aim of the research is to investigate the chemical composition and assess the level of soil contamination in the Solzan landfill and adjacent areas, including the coastal area of Lake Baikal. The study has determined the concentrations of Hg, As, Pb, Cu, Zn, Co, Ni, Cr, Mn, Sr, Ba, Rb, and TOC (total organic carbon) in soils near and far from the lignin sludge storage cell in the Solzan landfill. The obtained results show an extremely uneven distribution of elements across the study area. Abnormal element concentrations exceeding Clarke values, regional background concentrations, and soil sanitary and hygienic assessment criteria have been revealed near the lignin sludge storage cell. The uneven distribution of elements, high organic matter content, and wide range of pH values in the soils indicate an anthropogenic transformation of their chemical composition. The concentration of chemical elements in the soil decreases with distance from the impact area. Widespread soil contamination with zinc has been identified, with concentrations exceeding regional background levels and sanitary and hygienic standards both near and far from the source of contamination. According to the total pollution index, anthropogenically transformed soils in the study area mainly have an acceptable level of contamination. Nevertheless, moderately hazardous and hazardous levels of soil contamination have been detected both near and far from the Solzan landfill.
The paper presents a comprehensive study of the formation conditions, composition, and genesis of surface oxalate minerals discovered on technogenically disturbed soils in the Udachnaya kimberlite pipe mining area in Western Yakutia. The main objects include two soil profiles (R-29 and R-29/1) located on a gentle slope between the slurry pipeline and the tailings storage facility of the Udachninsky Mining and Processing Plant, directly above the deep disposal site for highly mineralized drainage brines. In profile P-29 located in the micro-depression an extensive continuous white coating was observed, whereas no similar formations were found in the adjacent micro-elevation (P-29/1). A complex of methods (X-ray diffraction, gas chromatography-mass spectrometry, atomic absorption spectrometry) reveals that the coating consisted of well-crystallized calcium oxalate dihydrate – weddellite (CaC2O4·2H2O). GC-MS analysis has confirmed biogenic nature of the formation: free oxalic acid has not been detected in the soil, although a significant amount of bound oxalates (2.7 mg/g dry weight) was found in the lichen biomass from the P-29 site. It has been concluded that coating formation is a protective response of lichen biota to technogenic impact due to a slurry pipeline failure. The major natural conditions of the process include the location in low microrelief parts, increased humidity due to seasonal permafrost thawing, a slightly alkaline soil environment, and the presence of carbonates in the parent rock. The trigger for abnormally intense oxalate formation is a technogenic pressure associated with emergency situations in the processing plant environment. The formation of visible macroscopic weddellite accumulations accompanied by lichen suppression is proposed to be utilized as a diagnostic indicator of ecological imbalance in permafrost zones subjected to intensive industrial impact.
Recycling and safe disposal of non-ferrous metallurgy wastes is an important and pressing environmental issue. The flotation tailings of cast copper smelting slags are a finely dispersed (less than 0.05 mm), mechanically activated material containing an increased concentration of copper, zinc and other chalcophilic elements. Bacterial leaching of the tailings has been carried out using silicate bacteria isolated from various sources (mineral aggregates of Sokoliny Kamen mountain, stale flotation tailings of copper smelting slags, mosses and lichens) and grown on Zak's medium at a temperature of 25±2 °C for a month. Inductively coupled plasma mass spectrometry (NexION 300S) is used to analyze the aqueous solution. It has been established that the solubility of the tailings components is influenced by the intensity of mixing, time factor, and type of strain. The maximum increase in the concentration of elements in the solution has been observed for iron on the fifth day, and for phosphorus, potassium, calcium, copper and lead on the nineteenth day. The migration activity of zinc and cadmium increases over time, which is associated with the disintegration of glass by silicate bacteria. The maximum mobility of silicon, zinc, cadmium, and phosphorus in aqueous solutions has been observed in the presence of silicate bacteria isolated from lichens. Iron and lead are more mobile in the presence of bacteria isolated from mosses. Sodium, potassium, calcium, and magnesium are more intensively leached in the presence of silicate bacteria isolated from stale tailings. Additionally, a comparative analysis of the solubility of flotation tailings in distilled water, in Zak's nutrient medium with and without silicate bacteria has been carried out. The obtained results are of interest for the development of methods for the disposal and safe burial of mineral tailings from non-ferrous metallurgy.
Orogenic granitoid massifs contain a variety of mafic inclusions (autoliths, xenoliths, mingling dikes). These include predominantly diorite-series inclusions, which visually contrast with the host tonalite-granite associations but have similar mineral and chemical compositions. The origin of these inclusions is controversial, and even the use of a range of modern precision methods does not always resolve the issue. The paper provides a petrographic, mineralogical, and geochemical review of orogenic granitoids from three massifs in two fold belts of the world: Adamello (Southern Alps), Akhun (Southern Urals), and Shabry (Middle Urals). The massifs are composed of tonalites and granites containing uniform angular and rounded diorite inclusions. The results obtained have shown that all rock varieties are composed of the same set of minerals of fairly similar composition (Ca-Na plagioclase, quartz, potassium feldspar, biotite, and amphibole), differing only in their ratios. Variations in mineral composition, particularly biotite, amphibole, and plagioclase indicate the significant influence of a fluid phase released during crystallization of a water-saturated melt. The use of amphibole-plagioclase thermobarometers has shown that the trapped diorites and host tonalite-granites crystallized under similar conditions (usually in the range of 700–800 °C and a pressure of 2–4 kbar). Diorite inclusions are autoliths, i. e., fragments of either crystallized rocks (angular inclusions) or trapped portions of partially crystallized melt (rounded inclusions) from the lower part of a single intrusive chamber. Patterns of their distribution and chemical composition in matrix tonalite-granites and diorite inclusions indicate their genetic unity. Patterns of change in the chemical composition of rocks of all massifs indicate a uniform course of their evolution. Binary diagrams show that similar geochemical patterns of rocks from all massifs indicate the uniformity of their evolutionary course.
This paper examines the results of magnetotelluric soundings in the eastern Fennoscandian Shield, in the Murmansk region and the Republic of Karelia, conducted as part of the FENICS experiment series in 2014, 2019, and 2024. During this series of deep electromagnetic sounding experiments, a magnetotelluric data set has been obtained using modern measuring stations. Observation points have been presented in the form of three profiles. Based on polar diagrams, as well as the Swift and Barr parameters, natural field data in the AMT-MT range corresponding to the two-dimensional approximation have been selected. The study has revealed the general nature of the apparent resistivity and impedance phase curves. A two-dimensional inversion of the magnetotelluric soundings along three profiles has been performed using the MARE2DEM program. At that two zones of low specific resistance have been identified, which are consistent with previous magnetotelluric studies along the SVEKA profile in the area of the Ladoga electrical conductivity anomaly, as well as with the general nature of regional geological structures.
This study examines the application of a suite of geophysical methods for assessing the physical and mechanical properties of rock masses in an open-pit chrysotile asbestos mining operation. The investigation employs three complementary techniques: seismic profiling (determination of compressional and shear wave velocities), electrical resistivity tomography (mapping of bulk electrical resistivity), and induced polarization surveying. The key innovation is a developed algorithm for synergistic integration of these methods, wherein results from each technique serve as prior information for interpreting data from other methods, achieving a 15–20 % increase in reliability compared to single-method approaches. Based on measured seismic velocities, Poisson's ratio and dynamic Young's modulus have been calculated. Electrical resistivity tomography cross-sections are used to verify seismic boundaries and refine density models, while induced polarization data enable identification of zones correlating with weakened segments of the rock mass. The resulting geophysical cross-sections visualize the spatial distribution of electrical resistivity, chargeability, and elastic modulus throughout the rock mass. Results from integrated interpretation are presented as two-dimensional sections and correlate with borehole geological data. The methodology has been tested at the Zhitikara chrysotile asbestos deposit. The obtained geophysical data are transmitted to the operations department for incorporation into the technological parameters of drilling and blasting operations. The methodology is recommended for implementation in the mining production cycle.
In modern petrological and geochemical literature, considerable attention is paid to the issue of the transformation of mantle rocks during their interaction with fluids and melts of various origins. An additional approach to building a picture of the mantle metasomatism process is to conduct modeling, both experimentally in high-pressure facilities and numerically, using specialized software packages for modeling phase equilibria in systems with a given gross composition, one of which is Perple_X. In order to simulate the reactions of phlogopitization and amphibolization during mantle metasomatism, an experimental study of the interaction of a model garnet-harzburgite with a water-carbon dioxide fluid H2O – CO2 – (K, Na)Cl at a pressure of 5 GPa and a temperature of 1 100 °C has been conducted. It has been found that the introduction of alkaline fluid components into the system provokes modal metasomatism: clinopyroxene and spinel are formed due to the decomposition of primary orthopyroxene and garnet, phlogopite appears in potassium systems, and amphibole appears in sodium systems. At the same time, the content of Al2O3 in the clinopyroxene increases, and a consistent change in the parameters Ca/(Ca+Mg+Fe) and Cr/(Cr+Al) in the garnet has been observed. The obtained results are compared with thermodynamic modeling (Perple_X), which reproduces the observed phase associations and clarifies the influence of H2O and K2O activities on the stability of minerals.
There is currently a lack of fundamental research on the impact of invasive aquatic species on the ecosystems of the Baltic Sea basin. The aim of the study is to review and summarize the existing knowledge on invasive aquatic species in the Baltic Sea basin. The present study has shown that at least 76 aquatic (or semi-aquatic) invasive species have been recorded on the eastern coast of the Baltic Sea. The microalga Prorocentrum minimum is so far the only invasive planktonic microalga known in the aquatic environment of the Baltic Sea basin. Researchers associate it with the accumulation of toxic substances, but increased toxicity of Prorocentrum has not been confirmed. Green algae Cladophora glomerata are one of the introduced macroalgae species found in the coastal areas of the Baltic. The list of invasive aquatic, semi-aquatic and coastal vascular plants includes 16 species. 30 invasive invertebrates have been recorded in marine coastal and inland freshwater environments of the Baltic Sea basin. Most of them are benthic and nektobenthic species, the rest are zooplankton (3 species) and parasites (1 species). A total of 28 invasive fish species belonging to eight families have been registered in the Baltic Sea basin, 27 of which were intentionally introduced for breeding. About half of the species (14–15) originate from different Asian regions (mainly from the Russian Far East and Siberia), five species from the Ponto-Caspian region, and four species from North America. The available information on invasive species requires the initiation of systematic and well-coordinated studies of their vital functions to provide the necessary scientific recommendations (e.g., vulnerability assessments) for the study of aquatic invaders. The present study should serve as a catalyst for further serious studies of invasive aquatic organisms in the Baltic Sea basin.
Knowledge of the temperature coefficient density distribution function for muons in the atmosphere is required to correctly account for the temperature effect in data from cosmic-ray muon telescopes. Temperature coefficient density distributions were found using continuous muon intensity observations obtained by muon telescopes at the surface and underground at depths of 7, 20, and 40 m.w.e. for the period from 2016 to 2018, and temperature data at 17 isobaric atmospheric levels obtained from upper-air soundings over the same period. Multivariate data analysis for the temperature effect has been conducted using principal component analysis and latent structure projection (LSP). Temperature coefficient density distributions are determined using one and two principal components, followed by a comparison of the results. With two principal components, the information content increases by ~4 % to 97.4 % (the majority of the original sample's variations are attributed to the first principal component, which contains 93.5 % of the information). However, the likelihood of noise contribution increases, while noise contribution is minimal with a single LC. The study has demonstrated that, in the presence of high noise levels (with low initial data quality), PLS analysis should be performed using only one principal component to increase the reliability of the final result.
As part of geophysical studies on the structure of the upper sedimentary sequence and the topography of the Arctic seas' shelves, a geophysical survey has been conducted on a section of the seabed in the coastal zone of the Rybachy Peninsula (Kola Peninsula). The analysis of seismic-acoustic data has revealed three seismic complexes, separated by unconformity surfaces. These complexes are tentatively identified using methods of sequence stratigraphy and reflect events of Quaternary glaciation and interglacial periods. Bathymetric data have allowed for the identification of multiple traces of ice scouring, which are of both modern and relict origin. The internal structure of the sedimentary sequence is not always typical of glacial deposits but is distinguished by considerable diversity and significant lateral variability. Furthermore, numerous signs of degassing in the form of pockmarks are evident. However, the absence of gas flares and stable migration pathways suggests that these pockmarks are currently inactive or of an impulsive nature.
The study analyzes the selenium content in the soils and plants of six municipal districts of the Chelyabinsk region: Chesmensky, Troitskiy, Plastovskiy, Uvelskiy, Yemanzhelinskiy and Etkulskiy. The samples were taken during the ripening period of agricultural crops (August – September 2023–2024) and included cereals, buckwheat and oilseeds. The soils are represented by ordinary and leached chernozems. The results have shown that the average selenium content varies from 0.1204 mg/kg in the Plastovsky district to 0.2679 mg/kg in the Yemanzhelinsky one. 42 % of samples are in the optimum zone, while 27 and 31 % of samples are in the zones of marginal insufficiency and selenium deficiency, respectively. Leached chernozems have higher selenium content (0.2009 mg/kg) compared to ordinary chernozems (0.1344 mg/kg). The analysis of correlations has revealed a negative relationship between the selenium content and the pH level of the soil, especially in the Chesmensky and Uvelsky districts. In other areas, this relationship is less pronounced, indicating the influence of other factors such as soil composition and anthropogenic influences. The study has shown that wheat and buckwheat accumulate selenium more efficiently than sunflower. The coefficient of biological absorption of selenium varies depending on the type of soil and the type of plant. The dependence of selenium content on the content of heavy metals (lead, zinc) in soils has been established. The influence of heavy metals on the selenium content varies depending on the area, which is related to the peculiarities of soil conditions and anthropogenic load. The data obtained emphasize the need for further research to develop differentiated strategies for improving the selenium status of agricultural land and increasing crop yields.
Rare Earth Elements (REE) and yttrium phosphates monazite-(Ce) and xenotime-(Y), fluorine carbonates synchysite-(Ce) and parisite-(Ce), silicate allanite-(Ce) have been studied in the metamorphic rocks in the area of the Pellapahk Cu-Mo porphyry occurrence in the north-western part of the Kolmozero-Voronya greenstone belt. Monazite-(Ce) and xenotime-(Y) form regular dissemination of fine grains in interstitions in rock-forming silicates and quartz in high-alumina gneisses and schists of the Porosozero series. In the altered granite-porphyry with Cu-Mo mineralization, REE minerals form zonal grains with monazite or xenotime in the core, surrounded by apatite-(F), then by fluorine carbonates synchysite-(Ce) and parisite-(Ce), and by silicates allanite-(Ce) and epidote in the outer parts. Besides, monazite-(Ce) and xenotime-(Y) have been found in inclusions in pyrite metacrystals. In the high-alumina schists, monazite and xenotime are of metamorphic genesis, they are formed due to re-distribution of REE and Y from isomorphic impurities in silicates to proper minerals. In the altered granite-porphyry, monazite and xenotime are supposed to be of magmatic genesis, but altered later during hydrothermal processes and partly or fully substituted by apatite, REE carbonates and silicates. Two schemes of heterovalent isomorphic substitutions are defined in monazites. Combination of "cheralite" and "huttonite" types of substitution (3LREE3+ + P 5+ ↔ Ca2+ + 2Th4+ + Si4+) prevails in monazites from the high-alumina schists, and "huttonite" type (REE3+ + P5+ ↔ Th2+ + Si2+) is more typical for monazites in altered granite-porphyry. Monazites of metamorphic genesis in high-alumina schists differ from primary magmatic in lower content of total REE and higher thorium impurity.
The paper considers one of the acute problems of the mining industry – the mechanization of mining processes during the development of flat thin ore deposits. The use of mobile complexes of self-propelled technological equipment (SPTE) using traditional stoping excavation schemes is limited by the requirement to ensure the value of extracted reserves close to their balance value. As a technical solution to this problem, selective excavation of ore mass and rock can be used. To increase the efficiency of this mining technique, selective excavation by double-cut stoping with blasting of the soil waste rock in the second cut has been considered. The given technical solution of the mining technology by double-cut stoping is characterized as an independent direction for improving the mining methods of flat thin ore deposits. For the conditions of the Karnasurt mine, the cost of performing the processes of blasting and ore mass haulage in each of the two considered options for stoping excavation has been calculated. At the same time, the option using of the SPTE complex has been considered with the width of the first cut of 3.5 and 7.0 m, and the second cut in both cases – 3.5 m. An analysis of the calculations performed has shown that doubling the width of the first cut makes it possible to reduce the cost difference between the blasting and haulage processes in the compared variants by 5 and 6 %.
The results of the assessment of dust deposited on the territory of Apatity are based on data obtained during previously performed numerical experiments to predict atmospheric pollution with fine dust (taking into account the deposition rate) from a man-made enrichment tailings storage facility. The models are developed in the COMSOL program for the ANOF-2 – Apatity tailings storage area. In the course of the study, an algorithm for calculating deposited dust on a horizontal surface has been implemented and a description of modern approaches to assessing secondary dust generation on urban highways has been presented. The interval and total dust masses accumulated during 1 hour on the territory of Apatity are calculated taking into account variations in the dusting area of 2–10 ha (with a random selection of dusty areas on the surface of the tailings dump) and wind speeds of 5–23 m/s. For an analytical description of the average total mass of dust accumulated in 1 hour under the conditions of the indicated ranges of model parameters, a power function of wind speed and a linear dependence on the dusting area are valid. The surface density of РM10 dust accumulated on the motor roads of the town of Apatity for 1 hour during dusting at the enrichment tailings storage facility has been estimated under conditions of variation in wind speed and dusting area. The above estimates of the hourly accumulation of dust on the roads of Apatity are necessary for environmental monitoring and the development of measures to reduce the impact of dust on public health.
Assessment of the natural waters' state requires an integrated approach including the determination of a wide range of physico-chemical parameters, on the basis of which specialized indexes and classifications are developed. Traditional analysis methods do not always take into account the effects of a number of pollutants, such as polyaromatic hydrocarbons (PAHs), which can have a significant impact on ecosystem processes and the quality of the aquatic environment. The paper presents the results of studies of small lakes in the Tomsk district (Tomsk region) affected by anthropogenic activity, which are freshwater reservoirs with electrical conductivity values of no more than 0.34 ms/cm, pH values up to 7.2. The content of inorganic macro components (cations and anions) does not exceed the maximum permissible concentrations (MPC), which indicates a satisfactory chemical composition of the water. The content of photopigments in bottom sediments has been analyzed in detail. A wide range of chlorophyll a concentrations from 0.35 to 2.3 µg/g has been recorded for all the studied samples, which may be due to the heterogeneity of the spatial distribution of phytoplankton and the trophic status of reservoirs. The individual and total content of 13 priority PAHs in surface waters has been studied in detail. Their total concentrations range from 0.37 to 0.54 µg/l, which indicates the presence of a moderate level of contamination; naphthalene, fluorene, phenanthrene and benz[a]anthracene predominate in the composition, which have the ability to bioaccumulate. Based on the data obtained, the water pollution index (WPI) has been calculated showing that all the lakes studied belong to the category of moderately polluted, approaching the border with polluted reservoirs. Thus, the study has made it possible to assess the impact of unintended anthropogenic impact on the chemical composition of surface waters of small reservoirs in the Tomsk region. Special attention has been paid to the role of polyaromatic hydrocarbons in the formation of the ecological state of water bodies and their contribution to changes in the values of the water pollution index. The results obtained demonstrate the importance of such studies for monitoring and protecting aquatic ecosystems.