ИЗМЕНЕНИЯ СОДЕРЖАНИЯ ТЯЖЕЛЫХ МЕТАЛЛОВ В ПОЧВАХ КОСТОМУКШСКОГО РЕГИОНА ПОД ВЛИЯНИЕМ ТЕХНОГЕНЕЗА

ИЗМЕНЕНИЯ СОДЕРЖАНИЯ ТЯЖЕЛЫХ МЕТАЛЛОВ В ПОЧВАХ КОСТОМУКШСКОГО РЕГИОНА ПОД ВЛИЯНИЕМ ТЕХНОГЕНЕЗА

Научная статья

Арестова И.Ю.^1, *, Опекунова М.Г.², Елсукова Е.Ю.³, Кукушкин С.Ю.⁴

^{1, 2, 3, 4} Санкт-Петербургский государственный университет, Санкт-Петербург, Россия

* Корреспондирующий автор (iarestova[at]mail.ru)

Аннотация

В статье представлены результаты мониторинговых исследований почв на участках Костомукшского государственного природного биосферного заповедника, города Костомукша и Костомукшского горно-обогатительного комбината. Отмечено низкое валовое содержание всех изученных металлов в естественных природных комплексах. Выявлены особенности вертикального распределения концентраций тяжелых металлов по почвенным горизонтам. Установлено резкое увеличение концентраций большинства изученных элементов во всех почвенных горизонтах в городской зоне и на территории комбината. Для ряда металлов отмечается изменение в характере распределения по горизонтам в сравнении с фоновыми почвами.

Ключевые слова: экологический мониторинг, загрязнение почв, тяжелые металлы, заповедник, Карелия.

CHANGES IN CONTENT OF HEAVY METALS IN SOILS OF KOSTOMUKSH REGION UNDER INFLUENCE OF TECHNOGENESIS

Research article

Arestova I.Yu.^1, *, Opekunova M.G.², Yelsukova E.Yu.³, Kukushkin S.Yu.⁴

^{1, 2, 3, 4} St. Petersburg State University, St. Petersburg, Russia

* Corresponding author (iarestova[at]mail.ru)

Abstract

The paper presents the results of monitoring studies of soils in the area of the Kostomuksha State Natural Biosphere Reserve, the city of Kostomuksha and the Kostomuksha mining and processing plant. The low gross content of all the studied metals in natural complexes was found. The authors revealed features of the vertical distribution of heavy metal concentrations over soil horizons. A sharp increase in the concentrations of most of the elements studied was found in all soil horizons in the urban zone and on the territory of the plant. For a number of metals, there is a change in the nature of the distribution over the horizons in comparison with background soils.

Keywords: environmental monitoring, soil pollution, heavy metals, nature reserve, Karelia. 

Введение

Территории, расположенные в зоне влияния крупных месторождений, давно являются объектом изучения содержания тяжелых металлов (ТМ) в сопряженных компонентах окружающей среды. Чаще всего в роли таких компонентов вступают почвы и растения, так как они отражают как краткосрочные, так и долгосрочные изменения, происходящие в природно-территориальных комплексах (ПТК). При этом в подобных исследованиях наибольший интерес представляет сравнение биогеохимических показателей участков, испытывающих постоянное воздействие, и территорий с минимальным уровнем нарушений, характеристики которых принимаются за фоновые.

Костомукшский горно-обогатительный комбинат (КГОК), расположенный на северо-западе Карелии, является ведущим комбинатом по добыче и переработке железной руды в России. Горные работы ведутся с 1982 года на базе Костомукшского месторождения, которое включает железистые кварциты, кристаллические сланцы различного состава и геллефлинты. В марте 2007 года начата добыча и переработка железной руды Корпангского месторождения. В 2008 году компания получила лицензию на разведку и обработку запасов руд Южно-Корпангского участка недр. В 2017 году добыча горной массы на комбинате составила более 51 млн м³, добыча руды – 36714000 т. [1]

На Костомукшском месторождении выделены четыре разновидности руд: щелочноамфиболовые, биотитовые, грюнеритовые и двуамфиболовые. Руды Корпангского месторождения имеют амфибол-магнетитовый состав. Содержание (Fe₃O4+SiO₂) в магнетитовых кварцитах составляет 87-92%. Оксиды Al₂O₃, MgO, CaO, K₂O+Na₂O содержатся в количестве n (0,1-1)%. Незначительную примесь составляют P2O5 (n 0,1%) и SrO (n 0,01%). Содержание S возрастает от богатых (0,03%) к бедным (0,9%) рудам и достигает 7% в силикатных кварцитах. Из всех элементов-примесей в рудах только содержания Mn (n 0,01%) и Ge (n 0,0001%) достигают кларковых значений. Руды содержат низкие концентрации экологически опасных элементов (в %): Zn, Ti – n 0,01; Cr, V, Pb, Ni, Y – n 0,001; Ga, Sc – n (0,001-0,0001); Co, Cu, Ge – n 0,0001; Ag – n (0,00001-0,000001). Содержание Au в рудах – n (1-0,001) г/т. [2] По имеющимся данным, на начало 2000-х годов аэротехногенные выбросы SO₂ предприятия составляли приблизительно 30 тыс. т/год, пыли 7-8 тыс. т/год. [3]

Жилые массивы города Костомукша расположены в десяти километрах к юго-западу от территории КГОКа. При строительстве жилых районов в городе были сохранены участки коренных лесных массивов, поэтому площадь зеленой зоны в городе более 30%. Площадь города составляет 15,5 км², население города на начало 2018 года - 29,4 тыс. человек.

В 30 км к юго-западу от города находится Костомукшский государственный природный биосферный заповедник общей площадью 49,2 тыс. га. Заповедник был создан в 1983 г, основной целью создания являлось сохранение северотаежных лесов европейской части России. До 40-х годов ХХ в. на данной территории располагались небольшие финские деревни и хутора, в настоящее время прямое антропогенное воздействие на охраняемые ПТК практически отсутствует. Кроме того, взаиморасположение заповедника, города и КГОКа таково, что все выбросы рассеиваются в противоположную от заповедника сторону. [4]

Исследования, связанные с воздействием выбросов комбината на различные компоненты окружающей среды, ведутся практически с начала его активной деятельности. С начала 1980-х годов в районе проводятся работы различными институтами Академии наук СССР (сейчас – РАН), научными коллективами С.-Петербургского Горного института, Петрозаводского университета и др. [5], [6], [7] Многолетние наблюдения состояния почв и растительности, включающие анализ биогеохимических характеристик, проводятся кафедрой геоэкологии и природопользования СПбГУ на участках, распложенных на территории Костомукшского государственного природного биосферного заповедника, в городе Костомукша и на территории Костомукшского ГОКа.

Основная часть проводящихся в районе исследований посвящена оценке содержания различных загрязняющих веществ в тех или иных компонентах окружающей среды в ПТК определенного уровня нарушенности – отдельно рассматриваются заповедные экосистемы, отдельно городские территории и территория ГОКа. Однако такая избирательность в выборе объектов исследования не позволяет оценить происходящие в среде изменения в аккумуляции загрязнений. Поэтому для выявления особенностей накопления ТМ почвами Костомукшского региона нами рассматривались ПТК всех уровней нарушенности.

Материалы и методы

Работы проводились в июне 2018 г. на площадках, относящихся к ПТК с разным уровнем нарушенности. В качестве фоновых были выбраны площадки на территории Костомукшского заповедника, расположенные недалеко от станции интегрированного мониторинга. В городе площадки выбирались в зеленых зонах и дворовых териториях, представляющих собой фрагменты исходных природных комплексов. На территории ГОКа площадки устанавливались в фрагментарных лесных сообществах вблизи источников техногенного воздействия. Всего было установлено 14 площадок в заповеднике, 14 площадок в городе и 6 площадок на территории ГОКа.

Пробоотбор почв осуществлялся в соответствии с ГОСТ 17.4.3.01-83. Размер пробных площадок составлял 20х25 м. На каждой площадке методом конверта отбирались пробы почв из горизонтов О (0-10 cм) и ВF (15-25 см). Названия горизонтов почв даны по Классификации почв России 2004 г. [9] На антропогенно нарушенных участках при отсутствии четко выраженного почвенного профиля пробы отбирались с глубины 0-5 см и 5-15 см. Для анализа вертикальной и латеральной миграции ТМ в почвах на нескольких площадках были заложены почвенные разрезы, пробы из которых отбирались по каждому горизонту. Всего было отобрано 66 проб почв из поверхностных горизонтов и 34 пробы почв из почвенных разрезов.

Определение подвижных форм ТМ в почвах проводилось с использованием почвенной вытяжки аммонийно-ацетатным буфером (pH=4,8) на атомно-абсорбционном спектрометре NOVAA-315. Анализ валового содержания тяжелых металлов в почвах проводился методом масс-спектрометрии с индуктивно-связанной плазмой (ИСП-МС) на приборе «ELAN-6100 DRC» с полным кислотным разложением проб по ПНД Ф 16.1:2.3:3.11-98. При обработке результатов измерений использовалась компьютерная программа обработки данных «TOTALQUANT», что обеспечивало пределы обнаружения элементов на уровне 0,0002-0,01 г/т. Анализ выполнялся в Центральной лаборатории ВСЕГЕИ им. А.П. Карпинского.

Для характеристики естественных особенностей содержания ТМ в почвах района рассчитывались коэффициент концентрации каждого элемента в изученных почвенных горизонтах (K_k=C_x:K_кларк) и коэффициент аккумуляции (K_a = C_x : ПДК (ОДК)_x), где С_х – концентрация элемента, Kкларк – кларк элемента в почве, ПДК – предельно допустимая концентрация элемента в почве, ОДК – ориентировочно допустимая концентрация элемента в почве.

Результаты и обсуждения

Определение валового содержания ТМ в почвах на территории заповедника (табл.1) показало, что в большинстве исследованных проб из поверхностных горизонтов концентрация Co ниже чувствительности метода.

 

Таблица 1 – Валовое содержание тяжелых металлов в почвах Костомукшского заповедника, мг/кг [10]

	Mn	Zn	Cu	Cr	Ni	Fe	Ba
Гор. O	263	91	6,56	10,81	7,45	4201	173
Гор. BF	211	54	5,53	38,36	15,63	20326	347
ПДК, (ОДК)	1500	(55)	(33)	–	(20)	–	–
Фон	282	–	–	47,30	–	17505	200

 

Изучение фоновых почвенных разрезов на подзолах иллювиально-железистых выявило накопление Со в нижних горизонтах BF и ВС. Cd в изученных образцах практически нигде не обнаружен. Исключение составляет единственная проба из поверхностного горизонта, где его содержание составляет 2,58 мг/кг и превышает установленную ОДК (0,5 мг/кг). Это может быть связано с тем, что Cd является очень подвижным элементом и способен переходить в больших количествах из твердой фазы в жидкую. Анализ вертикального распределения концентраций остальных изученных ТМ показал, что Cu распределяется по всему почвенному профилю равномерно, Mn, Pb и Zn в большинстве изученных почв тяготеют к верхним горизонтам, а Fe, Ba, Ni и Cr активно накапливаются в иллювиальных горизонтах.

Вместе с тем необходимо отметить естественный низкий региональный уровень содержания большинства изученных металлов в почвах по сравнению с кларками, что подтверждается рассчитанными коэффициентами концентрации ТМ для почвенных разрезов. Исключение составляют только Zn и Ba, концентрации которых в отдельных пробах из горизонта BF достигают кларка (табл. 2).

Таблица 2 – Значения K_k тяжелых металлов в почвах заповедника

	Mn	Zn	Cu	Cr	Ni	Fe	Ba
Гор. O	0,2-0,4	0,6-0,8	0,2-0,4	0,01-0,03	0,1	0,05	0,3
Гор. BF	0,4-0,5	0,5-1	0,2-0,4	0,2-0,3	0,3-0,4	0,4-0,6	0,9-1

 

Большинство изученных почв содержит ТМ в концентрациях ниже предельно допустимых. Это справедливо для всех изученных горизонтов в ненарушенных почвах. Исключения отмечены только для Zn и Ni, концентрации которых в отдельных пробах близки или даже незначительно превышают нормативы.

Анализ валового содержания ТМ в почвах городской зоны и территории комбината показал резкое увеличение концентраций большинства изученных элементов как в поверхностных горизонтах, так и в иллювиальных (табл. 3). В частности, на всех площадках и в черте города, и в зоне ГОКа активно накапливается Со. Единственный металл, концентрации которого в почвах заповедника выше, это Zn. Несмотря на это, его концентрация превышает ОДК во всех исследованных почвах.

Для ряда металлов отмечается изменение в характере распределения по горизонтам в сравнении с фоновыми почвами. Так Cu в городских почвах активно накапливается в поверхностных горизонтах, а в почвах ГОКа верхние горизонты обогащены Fe и Ni. Эти изменение связаны, скорее всего с прямым поступлением металлов с пылевыми выбросами от комбината и городского транспорта. Одновременно стоит отметить, что в большинстве почв комбината Pb активно накапливается именно в иллювиальном горизонте.

 

Таблица 3 – Валовое содержание тяжелых металлов в почвах г.Костомукша и КГОКа, мг/кг

	Mn	Zn	Cu	Cr	Ni	Fe	Ba	Co	Cd	Pb
	Почвы г.Костомукша
Гор. O	400	88	13,7	25,9	15,3	9240	279	3,99	0,3	26,3
Гор. BF	371	31,6	6,3	76,6	20,6	17045	461	9,83	0,05	13,3
	Почвы КГОКа
Гор. O	800	73,7	14,6	50,9	39,4	55679	250	12	0,16	8,71
Гор. BF	444	59,1	14,4	58,5	28	36855	396	9,77	0,07	12,5
ПДК, (ОДК)	1500	(55)	(33)	-	(20)	-	-	-	(0,5)	(32)

 

Анализ пространственного распределения загрязнений ТМ по территории города и ГОКа выявил несколько наиболее загрязненных участков, где содержания большинства изученных элементов резко возрастают. На территории ГОКа такой является площадка вблизи железнодорожной ветки, среди исследованных городских площадок наибольшее загрязнение поверхностных горизонтов почв отмечено во дворах жилых домов в центре и в юго-восточной части города.

Кроме того, в почвах этих площадок отмечены изменения в распределении ТМ по горизонтам в сравнении с фоновыми территориями. Так на площадках, расположенных в городе, такие металлы как Cr, Co, Ni, Zn, Cd, Pb сосредоточены в верхних горизонтах почв, что можно объяснить поступлением металлов из атмосферы от мобильных и стационарных источников. На территории ГОКа изменение накопления ТМ имеет более сложный характер, часть элементов закрепляется в верхних горизонтах (V, Cr, Cu, Fe), а другие, наоборот, активно накапливаются в нижних (Zn, Cd, Pb, Mn). Это связано с загрязнением поверхностных горизонтов в результате постоянного поступления металлов из атмосферы с газопылевыми выбросами от карьера, комбината и железнодорожного транспорта, а также с нарушением почвенного профиля при разработке и обустройстве территории.

Для анализа вертикальной и латеральной миграции ТМ в зоне антропогенного воздействия были проанализированы пробы из почвенных разрезов, расположенных в 5 км к юго-западу от территории КГОКа (табл. 4).

 

Таблица 4 – Валовое содержание ТМ в почвенных разрезах возле КГОКа, мг/кг

Горизонт	Cu	Zn	Ni	Sc	Cd	Cr	Pb	Co	Mn	Fe	V	Ba
Почвенный разрез 34-1. Вершина моренного холма.
O	10,3	43,2	35,7	2,42	0,35	34	14,5	4,71	310	37100	59,2	147
H	17,9	33,5	5,34	5,01	0,02	23,6	8,25	2,26	178	3885	30	480
BF	13,2	66,3	22,4	11,4	0,06	68,5	11,7	17,2	496	20405	99,5	486
BC	8,6	49,1	38,9	14,2	0,06	115	8,9	16,4	550	14945	94,5	451
Почвенный разрез 34-2. Нижняя часть склона моренного холма С-СЗ экспозиции
O	7,8	57,9	22,1	1,18	0,22	17,2	14	2,3	604	23100	36,8	119
H	<1	19,1	3,9	6,12	0,04	18,1	18	1,5	201	3255	32,4	652
BF	11,2	34,5	30,3	11,7	0,04	88,8	8,9	13,3	473	14350	81,4	496
Почвенный разрез 34-2. Подножье моренного холма
O	47,0	147	15,9	4	0,21	25,1	24,3	6,8	155	9940	32,3	242
H	2,3	19	14,5	8	0,02	50,5	10,4	5,9	294	6475	54,1	530
BF	5,8	27,5	27,6	11,6	0,03	91	7,6	11,9	434	12810	77,2	406
ОДК (ПДК)
	33	55	20	-	0,5	-	(32)	-	(1500)	-	150	-

Выявлено отличие в распределении ряда элементов по горизонтам от фоновых территорий. В частности, в поверхностных горизонтах данных разрезов идет накопление Fe, Pb, Ni и Cd, во всех горизонтах увеличивается содержание Co, в иллювиальных горизонтах идет активное накопление Mn. Накопление Zn в почвах имеет неустойчивый характер и зависит от положения в рельефе. Для остальных металлов характерно быстрое вымывание из поверхностных горизонтов, миграция и накопление именно в нижележащих иллювиальных горизонтах.

Исследования содержания подвижных форм ТМ в фоновых почвах на территории Костомукшского заповедника позволили выделить две группы металлов (табл. 5).

 

Таблица 5 – Содержание подвижных форм тяжелых металлов в почвах, мг/кг [11]

	Co	Mn	Zn	Cu	Cd	Pb	Cr	Ni	Fe
Гор. O	0,29	185	16,61	2,78	0,12	2,98	0,07	0,75	11
Гор. BF	0,37	15	2,86	2,88	0,04	0,72	0,51	0,26	212
ПДК	5	60	23	3	0,5	6	6	4	–

 

Первая группа, объединяющая Mn, Zn, Cd, Pb и Ni, характеризуется преимущественным накоплением в поверхностных горизонтах почв и резким снижением концентраций в иллювиальных горизонтах. Вторая группа включает Fe, Cо и Cr, чьи концентрации, наоборот, возрастают в иллювиальных горизонтах. Стоит отметить, что несмотря на различия в аккумуляции содержание подвижных форм большинства ТМ в почвах заповедника существенно ниже ПДК. Исключение составляет только Mn, концентрации которого в поверхностных горизонтах превышают ПДК в 3 раза Степень подвижности ТМ зависит от формы рельефа, а также от почвенного горизонта. В поверхностных горизонтах исследованных почв подвижность элементов выше, чем в иллювиальных. Это может быть обусловлено влиянием атмосферных осадков, а также возможным аэротехногенным переносом.

Заключение

Проведенные исследования подтвердили низкий природный региональный уровень валового содержания большинства изученных металлов в ненарушенных почвах по сравнению с кларками. Концентрации этих ТМ ниже предельно допустимых во всех почвенных горизонтах в фоновых природных комплексах. Исключения отмечены только для Zn и Ni, концентрации которых в отдельных пробах близки или даже незначительно превышают нормативы. На территориях, подверженных антропогенному воздействию происходит не только резкое увеличение концентраций большинства изученных ТМ, но и изменения в их аккумуляции различными горизонтами почв.

По распределению подвижных форм тяжелых металлов в почвенном профиле выделены две группы элементов. Металлы первой группы накапливаются в поверхностных горизонтах почв, тогда как металлы, относящиеся ко второй группе, тяготеют к более глубоким горизонтам.

Финансирование Исследования проведены при поддержке гранта РФФИ 18-05-00217 «Биогеохимические индикаторы техногенной трансформации потоков тяжелых металлов в ландшафтах».	Funding The studies were supported by the RFBR grant 18-05-00217 “Biogeochemical indicators of technogenic transformation of heavy metal flows in landscapes”.
Конфликт интересов Не указан	Conflict of Interest None declared

Список литературы / References

1. Официальная статистическая информация ПАО «Северсталь» [Электронныйресурс]  – URL.: http://karelskyokatysh.severstal.com/rus/customers/products/index.phtml (дата обращения: 23.11.2019).
2. Галахина Н.Е. Оценка воздействия предприятий горнодобывающей железорудной промышленности на водную среду с учетом природно-техногенных факторов формирования вод : дисс. … канд. хим. наук : 03.02.08 : защищена 18.04.18 : утв. 10.08.18 / Галахина Наталия Евгеньевна. – Петрозаводск: Ин-т водных проблем Севера КНЦ РАН, 2016. – 160 с.
3. Пантелеева Я.Г. Геохимические изменения окружающей среды в зоне влияния горнопромышленного комплекса ОАО "Карельский окатыш" : дисс. … канд. геол.-мин. наук : 25.00.09 : защищена 16.10.2009 : утв. 23.04.2010 / Пантелеева Яна Геннадьевна. – СПб: СПбГорный Университет, 2009. – 242 с.
4. Арестова И.Ю. Мониторинговые исследования фоновых территорий Костомукшского заповедника / Арестова И.Ю., Опекунова М.Г., Елсукова Е.Ю., Кукушкин С.Ю. // «Роль научно-исследовательской работы в управлении и развитии ООПТ» - Материалы Всероссийской научно-практической конференции, посвященной 50-летию со дня образования Байкальского государственного природного биосферного заповедника (пос. Танхой, 14–15 октября 2019 г.). – Иркутск: Изд-во Института географии им. В.Б. Сочавы СО РАН, 2019 г. – с. 28-32.
5. Германова Н.И. Динамика биологической активности лесных почв в зоне действия Костомукшского ГОКа / Германова Н.И., Медведева М.В. // Биоэкологические аспекты мониторинга. Петрозаводск: КНЦ РАН, Институт леса, 2001. — с. 192-201.
6. Федорец Н.Г. Воздействие эмиссий Костомукшского горно-обогатительного комбината на лесные подстилки сосняков в северотаежной подзоне Карелии / Федорец Н.Г., Солодовников А.Н. // Труды Карельского научного центра РАН, № 6. – с. 143–152.
7. Виноградова А.А. Антропогенная нагрузка на экосистемы Костомукшского природного заповедника: Атмосферный канал. / Виноградова А.А., Иванова Ю.А. // М., Изд-во Физматлит, 2013. – 84 с.
8. Новиков С.Г. Содержание тяжелых металлов в почвах города Костомукши / Новиков С.Г. // Современные проблемы науки и образования (электронный научный журнал), 2016, № 5.
9. Шишов Л.Л. Классификация и диагностика почв России / Л.Л. Шишов, В.Д. Тонконогов, И.И. Лебедева, М.И. Герасимова. – Смоленск: Ойкумена, 2004. – 342 с.
10. Елсукова Е.Ю. Биогеохимический мониторинг особо охраняемых природных территорий на примере Государственного природного заповедника "Костомукшский"/ Елсукова Е.Ю., Опекунова М.Г., Арестова И.Ю., Кукушкин С.Ю, Карандашева Н.И., Решетняк В.А. // Актуальные вопросы биогеографии: Материалы международной конференции (Санкт-Петербург, Россия, 9-12 октября 2018г.). – Санкт-Петербургский государственный университет. СПб, 2018.- с. 130-132.
11. Elsukova E. Ecological monitoring of special protected areas by the example of Kostomuksha Nature Reserve / Elsukova E., Opekunova M., Arestova I., Kukushkin S. // International Multidisciplinary Scientific GeoConference Surveying Geology and Mining Ecology Management, SGEM 17, Water Resources. Forest Marine and Ocean Ecosystems. – 2017. P. 805-812.

Список литературы на английском языке / References in English

1. Oficial'naya statisticheskaya informaciya PAO «Severstal'»  [Official statistical information of PAO Severstal] [Electronic resouce]  / –URL.: http://karelskyokatysh.severstal.com/rus/customers/products/index.phtml (accessed: 23.11.2019) [in Russian]
2. Galakhina N.E. Ocenka vozdejstviya predpriyatij gornodobyvayushchej zhelezorudnoj promyshlennosti na vodnuyu sredu s uchetom prirodno-tekhnogennyh faktorov formirovaniya [Assessment of the impact of mining enterprises in the iron ore industry on the aquatic environment, taking into account natural and technogenic factors of water formation] : PhD Thesis in Chemistry: 03.02.08 : defended 18.04.18 : approved 10.08.18 / Galakhina Natalia. – Petrozavodsk, Northern Water Problems Institute Karelian Research Centre RAS, 2016. - 160 p. [in Russian]
3. Panteleeva Ya.G. Geohimicheskie izmeneniya okruzhayushchej sredy v zone vliyaniya gornopromyshlennogo kompleksa OAO "Karel'skij okatysh" [Geochemical changes in the environment in the zone of influence of the mining complex of OJSC Karelsky Okatysh] : PhD Thesis in in Geology and Mineralogy : 25.00.09 : defended 16.10.2009 : approved 23.04.2010 / Panteleeva - St. Petersburg, Saint-Petersburg Mining University, 2009. – 242 p. [in Russian]
4. Arestova I.Yu. Monitoringovye issledovaniya fonovyh territorij Kostomukshskogo zapovednika [Monitoring research of the background territories of Kostomukshsky reserve] / Arestova I.Yu., Opekunova M.G., Elsukova E.Yu., Kukushkin S.Yu. // “The role of scientific research in the management and development of protected areas” - Materials of the All-Russian Scientific and Practical Conference dedicated to the 50th anniversary of the Baikal State Nature Biosphere Reserve (Tankhoi village, October 14–15, 2019). / Irkutsk: Publishing House of the Institute of Geography. V.B. Sochava SB RAS, 2019 - p. 28-32. [in Russian]
5. Germanova N.I. Dinamika biologicheskoj aktivnosti lesnyh pochv v zone dejstviya Kostomukshskogo GOKa [The dynamics of the biological activity of forest soils in the zone of operation of the Kostomuksha GOK] / Germanova N.I., Medvedeva M.V. / Bioecological aspects of monitoring. Petrozavodsk: KSC RAS, Forest Institute, 2001. - P. 192-201. [in Russian]
6. Fedorets N.G. Vozdejstvie emissij Kostomukshskogo gorno-obogatitel'nogo kombinata na lesnye podstilki sosnyakov v severotaezhnoj podzone Karelii [The impact of emissions of the Kostomuksha mining and processing plant on forest litter of pine forests in the northern taiga subzone of Karelia] / Fedorets N.G., Solodovnikov A.N. / Transactions of the Karelian Scientific Center of the Russian Academy of Sciences, No. 6. 2013. - P. 143–152. [in Russian]
7. Vinogradova A.A. Antropogennaya nagruzka na ekosistemy Kostomukshskogo prirodnogo zapovednika: Atmosfernyj kanal [Anthropogenic load on ecosystems of Kostomuksha natural reserve: Atmospheric canal] / Vinogradova A.A., Ivanova Yu.A. /M., Fizmatlit Publishing House, 2013. - 84 p. [in Russian]
8. Novikov S.G. Soderzhanie tyazhelyh metallov v pochvah goroda Kostomukshi [The content of heavy metals in soils of the city of Kostomuksha] / Novikov S.G. // Modern problems of science and education (electronic scientific journal), 2016, No. 5. [in Russian]
9. Shishov L.L. Klassifikaciya i diagnostika pochv Rossii [Classification and diagnostics of Russian soils] / L. Shishov, V. Tonkonogov, I. Lebedeva, M. Gerasimova. – Smolensk, 2004. – 342 p. [in Russian]
10. Elsukova E.Yu. Biogeohimicheskij monitoring osobo ohranyaemyh prirodnyh territorij na primere Gosudarstvennogo prirodnogo zapovednika "Kostomukshskij" [Biogeochemical monitoring of specially protected natural areas on the example of the Kostomuksha State Nature Reserve] / Elsukova E.Yu., Opekunova MG, Arestova I.Yu., Kukushkin S.Yu., Karandasheva N.I., Reshetnyak V.A. // Actual issues of biogeography: Materials of an international conference (St. Petersburg, Russia, October 9-12, 2018). - St. Petersburg State University. SPb, 2018.- p. 130-132. [in Russian]
11. Elsukova E. Ecological monitoring of special protected areas by the example of Kostomuksha Nature Reserve / Elsukova E., Opekunova M., Arestova I., Kukushkin S. // International Multidisciplinary Scientific GeoConference Surveying Geology and Mining Ecology Management, SGEM 17, Water Resources. Forest Marine and Ocean Ecosystems. – 2017. P. 805-812.