Features of carbon forms distribution in peat deposits of an oligotrophic bog in the Kondinsky Lakes Nature Park

Полный текст

Используемые сокращения

РОУ – растворенный органический углерод

РОВ – растворенное органическое вещество

УБВ – уровень болотных вод

ГВ – коэффициент соотношения гумусовых веществ

ГК – гуминовые кислоты

ФК – фульвокислоты

МВ – коэффициент среднего молекулярного веса

УВ – увлажнение

БЗ – активное богатство и засоление почв

ВВЕДЕНИЕ

Торфяники являются крупнейшим резервуаром долгосрочного накопления атмосферного углерода. Благодаря наличию торфа и болотной растительности болота связывают в несколько раз больше углерода на единицу площади, чем другие наземные экосистемы, включая леса [Minayeva, Sirin, 2012]. Баланс углерода природных экосистем определяется соотношением процессов поглощения углерода растениями и его выделения [Golovatskaya et al., 2022]. Поток углерода из болот, особенно в северных широтах, значительно влияет на климат [Lee, Romero, 2023]. Различные колебания климата, например повышение температуры, могут увеличивать поток углерода из болот, что ведет к положительной обратной связи – дополнительному нагреву атмосферы [Lee, Romero, 2023]. Таким образом, составление прогноза изменения климата планеты и построение единой модели, учитывающей содержание парниковых газов, является сложной задачей.

В процессе роста болотной растительности, а также под воздействием редуцентов углерод частично возвращается в атмосферу в виде углекислого газа и метана [Leonova, 2023]. Оставшаяся часть углерода депонируется на длительный срок в торфяной залежи. Для торфяников северного полушария средняя скорость накопления торфа оценивается в 0.067±0.018 см год⁻¹ (в диапазоне от 0.013 до 0.332 см год⁻¹), а скорость накопления углерода составляет 38.56±12.21 г м⁻² год⁻¹ (в диапазоне от 28.46 до 57.91 г м⁻² год⁻¹) [Turunen et al., 2001].

В России наибольшая площадь болотных экосистем (592 440 км²) сосредоточена на территории Западной Сибири [Sheng et al., 2004], где накопленное количество углерода оценивается в 473-621 Гт [Yu et al., 2010]. Наиболее древние слои торфа имеют возраст 9-12 тысяч лет, а средняя его мощность составляет 3-4 м [Zemcov et al., 1998; Macdonald et al., 2006].

Динамика накопления углерода не остается постоянной на протяжении всего существования болот. На данный процесс оказывают влияние многие факторы, главным из которых является смена растительности и изменение гидрологического режима болот [Tsyganov et al., 2021; Zarov et al., 2023]. Таким образом, комплексное изучение структуры болот позволяет оценить скорость накопления углерода и воссоздать картину палеоклимата (по составу растительных остатков) и палеогидрологических условий (по остаткам раковинных амеб) [Tsyganov et al., 2021; Halaś et al., 2023; Marcisz et al., 2020]. Данные, полученные при исследовании болот, позволят понять взаимообратные связи между биотой и климатом и прогнозировать реакцию экосистем на климатические изменения в будущем. Стратиграфический анализ торфяных залежей позволяет выявить региональные аспекты генезиса и развития болотных массивов.

В процессе разложения органического вещества часть элементов из подстилки и торфа переходит в растворимую форму, что приводит к перераспределению растворенной органики по торфяной толще. Фракции, проходящие через фильтр с размером пор 0.45 мкм, известны как растворенный органический углерод (РОУ) [Kolka et al., 2008]. Типы растительных сообществ влияют на процессы формирования торфа и определяют состав и концентрацию РОУ [Armstrong et al., 2012]. Также концентрация РОУ зависит от среднегодовой температуры воздуха [Freeman, 2001], концентрации углекислого газа в атмосфере и уровня болотных вод (УБВ) [Kolka et al., 2008]. Кроме того, существуют локальные факторы, контролирующие вынос РОУ с торфяников: размер водосборной площади, глубина залегания грунтовых вод и тип почвы [Worral et al., 2008].

Изучение концентрации РОУ важно для определения углеродного баланса торфяников [Waddington, Roulet, 1997]. Ожидаемое глобальное потепление может способствовать увеличению концентрации РОУ в воде, стекающей с заболоченных водосборов [Freeman, 2001]. Это в свою очередь приведет к увеличению концентрации РОУ в водах более крупных водотоков, где его часть минерализуется и возвращается в атмосферу в виде парниковых газов [Frey, Smith, 2005]. При этом часть растворенной органики с высокой долей ароматических соединений в составе более устойчива к разложению и способна мигрировать на значительные расстояния. Ароматичность является важным показателем, характеризующим состав растворенного органического вещества (РОВ), повышение значений которого отражает увеличение скорости переноса и распределения загрязняющих веществ в водные системы [William et al., 2016]. Для оценки ароматичности РОВ используется показатель SUVA₂₅₄ (specific ultraviolet absorption) – соотношение нормализованного значения поглощения при 254 нм к содержанию РОУ [Weishaar et al., 2003; Pittaway, 2014]. Следовательно, чем выше содержание ароматических соединений в составе органических веществ, тем выше поглощение в данном диапазоне. В состав РОВ входят гумусовые вещества, которые представляют собой смесь фракций с различной молекулярной массой и химическим составом; они обладают высокой миграционной способностью и участвуют в углеродном цикле биосферы [Bolin et al., 1979; Inisheva et al., 2013]. Известно, что гумусовые вещества традиционно делят на гуминовые кислоты (ГК) и фульвокислоты (ФК). ГК являются гидрофобной фракцией с более высокой молекулярной массой и размером от 2 до 50 нанометров, которая подавляет подвижность и способствует его накоплению. ФК определяют как гидрофильную низкомолекулярную фракцию с размером около 2 нанометров. Данное соединение характеризуется высокой подвижностью [Bai et al., 2018]. Доля ФК в болотных водах на порядок превышает содержание ГК, что связано с высоким вкладом в их структуру карбоксильных групп и фенольных оксигрупп [Thurman, 1985; Bolin et al., 1979]. В ходе роста торфяной залежи в глубинных слоях происходит потеря неустойчивых ФК и накопление ГК, которые со временем зрелости торфа становятся более концентрированными и устойчивыми [Savicheva et al., 2003; Moore, 1987].

Средний молекулярный вес РОВ является ключевым фактором, влияющим на физические, химические и биологические характеристики. Он определяет фотореакционную способность, оптические свойства и биодоступность РОВ, а также сорбционную способность тяжелых металлов и гидрофобные связывающие свойства [Chin, 1997]. Данный показатель изменяется в зависимости от качества и количества молекул, а также от фракций органических и неорганических частей и функциональных групп, образующих эти соединения [Uyguner, Bekbolet, 2005, Lou, Xie, 2006; Olefeldt, Roulet, 2012].

Значения pH имеют значительное влияние на высвобождение РОУ во время стока, поскольку он определяет сорбционную способность и растворимость органического углерода [Worrall F et al., 2008].

На территории природного парка «Кондинские озера» им. Л.Ф. Сташкевича проведены единичные работы, описывающие скорость прироста мхов [Kosykh et al., 2017], мониторинг дереворазрушающих микосообществ [Stavishenko, 2008], анализ геохимических особенностей торфяных залежей [Golovackaya, Smirnova 2021], флористические исследования [Bespalova, 2022; Korotkikh et al., 2022]. При этом отсутствуют исследования, направленные на оценку содержания углерода в торфяной залежи, концентрации РОУ и спектральные характеристики. Несмотря на то, что разный ботанический состав и экологические условия торфообразования являются основной причиной широкой изменчивости химического состава органического вещества торфов, работы в данном направлении весьма ограничены.

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

Целью данного исследования была оценка влияния ботанического состава и физико-химических параметров торфяной залежи на количество и состав РОУ.

Для достижения данной цели были поставлены следующие задачи:

1. Определить физико-химические показатели торфа;
2. Изучить стратиграфию торфяной залежи;
3. Определить содержание органического углерода в торфяной залежи;
4. Определить концентрацию и состав РОУ в торфяной залежи.

ОБЪЕКТЫ И МЕТОДЫ ИССЛЕДОВАНИИ

Исследования проводились в Советском районе Ханты-Мансийского автономного округа – Югры на территории природного парка «Кондинские озера» им. Л.Ф. Сташкевича. Территория расположена на левом берегу реки Конды и относится к Кондинской физико-географической провинции в среднетаежной подзоне тайги [Voronov et al., 1971].

Климат природного парка «Кондинские озера» характеризуется коротким, сравнительно теплым летом и суровой продолжительной зимой. Продолжительность зимнего периода составляет около 25 недель. Самый холодный месяц – январь со средней температурой -19°С. Самый теплый месяц – июль, средняя температура которого достигает 17.1°С. Природный парк находится на территории с достаточным увлажнением. Среднегодовое количество осадков составляет 500-550 мм, основное количество которых выпадает в теплый период года – с апреля по октябрь (82%) с максимумом в июле – августе. Минимальные значения выпадающих осадков приходятся на зимние месяцы [Levasheva et al., 2010].

В растительном покрове преобладают сосновые зеленомошные, лишайниковые и елово-кедровые зеленомошные леса. Среди болотных сообществ выделены следующие типы: сосново-кустарничково-сфагновые сообщества (рямы); кустарничково-сфагновые сообщества с чередующимися грядами и мочажинами; осоково-сфагновые сообщества на транзитных водотоках и в комплексных мезотрофных болотах; вахтово-сфагновые сообщества в участках выклинивания грунтовых вод и топей; осоково-вахтово-моховые сообщества в увлажненной части топей [Korotkikh et al., 2018].

Торфяная залежь имеет глубину 3-4 м с минеральным дном, сложенным супесями, песками и суглинками [Syr'eva, 2017; Lapshina et al., 2020]. Основными растениями-торфообразователями в современном растительном покрове являются Sph. divinum, Sph. balticum, Sph. angustifolium, E. vaginatum, Scheuchzeria palustris, C. limosa, C. rostrata и древесные остатки.

Отбор торфа

Торфяные керны были отобраны в весенний период (март 2023 года) с помощью пробоотборника типа «Русский бур» (производства Eijkelkamp, длина 0.5 м, ширина 5 см) на всю глубину торфяной залежи. В мочажинных сообществах было отобрано три торфяных керна: Cond-1 (60.852855° с.ш 63.502405° в.д., общая глубина 3.4 м), Cond-2 (60.854300° с.ш. 63.505760° в.д., общая глубина 2.9 м) и Cond-3 (60.855770° с.ш. 63.509037° в.д., общая глубина 3.4 м). Каждый образец торфяного керна длиной 50 см был помещен в С-образные кассеты из пенополиэтилена и обтянут стретч-пленкой, промаркирован и отправлен на хранение при температуре -20°С в лабораторию (Югорский государственный университет) [Batrshina et al., 2023]. Верхние 50 см не были отобраны для анализа в связи с трудностью отбора проб.

Пробоподготовка

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

Влажность торфа

Сегмент торфа отправляли в сушильный шкаф VO400cool (Memmert, Германия) для сушки при температуре 105℃ в течение суток, после чего высушенный торф взвешивался. Влажность торфа была определена по разнице массы влажного и сухого образца.

Элементный анализ (CNH)

Высушенный торф был перемолот до однородного состояния на измельчителе NIMA NM-8300 (Китай). Определение общего содержания органического углерода производилось на элементном анализаторе EA-3000 (EuroVector, Италия) [Batrshina et al., 2023]. Торфа верховых болот Западной Сибири не содержат карбонатов, поэтому общее содержание углерода равно его органической доле [Zarov, 2013]. В ходе анализа образец окислялся на катализаторе (чистая медь и оксид хрома) при температуре 1020°С в токе гелия и избытке чистого кислорода. Разделение смеси происходило хроматографическим методом с определением компонентов на детекторе по теплопроводности (TCD). Калибровка прибора осуществлялась с использованием стандартов Atropine (C = 70.56%, N = 4.84%, H = 8.01%, O = 16.59%) и Alfalfa (С = 42.35%, H = 5.95%, N = 2.91%). Всего было проанализировано 97 торфяных проб.

Ботанический анализ и степень разложения

Образец торфа промывали проточной водой через сито с ячейкой 0.25 мм. Для определения остатков растений использовался бинокулярный микроскоп (увеличение 10-40×; Zeiss Axiostar, Йена, Германия). Методика выполнялась в соответствии с протоколами [Mauquoy, Van Geell, 2013] с применением банка данных ключевых образцов, то есть коллекции остатков растений, которые были найдены в регионе и использовались для определения ботанического состава торфа. Наличие каждого типа растительных остатков было выражено в процентах от общей массы, а типы торфа были определены на основе доминирования видов растений по Матухину [Matukhin et al., 2000]. Степень разложения торфа выражалась как доля гумифицированных веществ к общему количеству неразложившихся растительных остатков. Всего было проанализировано 97 торфяных проб.

Расчет запаса органического углерода в слое торфа

Запас органического углерода в каждом 10-сантиметровом слое торфяной залежи рассчитывался как произведение доли органического углерода и объемного веса торфа, умноженное на мощность торфяного слоя (10 см) на площадь 1 м². Послойное суммирование запасов углерода каждого 10-сантиметрового слоя торфяной залежи позволило определить запас углерода на 1 м² [Leonova, 2023; Batrshina et al., 2023].

Показатели увлажнения (УВ) и активного богатства и засоления почвы (БЗ)  по шкалам Л.Г. Раменского

Для каждого образца были рассчитаны показатели увлажнения (УВ) – оценки водообеспеченности растений в различных природных зонах и местообитаниях и активного богатства и засоления почвы (БЗ) – наличие питательных элементов для растений в доступной форме. Для расчета показателей увлажнения был использован метод взвешенного усреднения согласно экологическим шкалам Л.Г. Раменского. Расчет производился в интегрированной ботанической информационной системе IBIS 7.2 [Zverev, 2007]. Данная методика была применена к ботаническим описаниям торфа в количестве 97 образцов.

Определение РОУ

Существуют частичные рекомендации относительно методов отбора проб и анализа, что ограничивает получение стандартизированных результатов [Peacock M. et. al., 2015]. Для извлечения РОУ из торфяного сегмента навеска торфа (~10-20 г) помещалась в коническую колбу объемом 100 мл, в которую добавлялось 50 мл деионизированной воды. Содержимое колб помещалось в перемешивающее устройство (ЭКРОС 6410) при частоте 200 об/мин в течение 30 минут. Затем содержимое колб переливалось на сито (диаметр ячеи 0.25 мм), установленное на керамической чаше, и отжималось керамической лопаткой [Shanyova, 2023]. Полученный экстракт отбирался в шприц и пропускался через мембранный фильтр (0.45 мкр) в стеклянные виалы (2 мл) [Leenheer, Croue, 2003; Peacock, 2015].

Измерения проводились на элементном анализаторе Flash-2000 (Thermo Scientific, США). Принцип определения содержания органического углерода основан на сжигании образца на катализаторе (чистая медь и оксид меди) при высокой температуре в присутствии избытка кислорода в токе гелия. Газообразные продукты горения проходят через систему восстановительных реакторов и поглотителей, разделяются на хроматографической колонке и детектируются на катарометре. Калибровка осуществлялась по стандарту Urea (C = 22.56%, N = 7.84%, H = 15.01%, O = 18.59%). Всего было проанализировано 97 водных экстрактов.

Определение спектральных характеристик РОВ

Отфильтрованный экстракт отбирался автоматической лабораторной пипеткой (Thermo Scientific до 10 мл) объемом 1 мл и разбавлялся 2 мл деионизированной воды в стеклянной пробирке [Freeman C. et al., 2004]. Подготовленный образец помещался в рабочую кварцевую кювету с толщиной пропускания 1 см, а аналогичная кювета сравнения заполнялась деионизированной водой [Rakhova, 2023]. Измерения проводились на приборе UV/Visible Spectrophotometer T8DCS (PERSEE, Китай) при длинах волн: 250, 254, 365, 400 и 600 нм [Peuravuori, Pihlaja, 1997; Uyguner, Bekbolet, 2005; Olefeldt, Roulet, 2012]. В программе UV-Win строилась базовая линия по отношению к деионизированной воде, определяющая нулевое светопоглощение. Значения показателей образцов воды были получены как отклонение значений от базовой линии. В качестве выходного сигнала был использован показатель оптической плотности.

Для изучения показателя SUVA₂₅₄ используется соотношение между поглощением при 254 нм к содержанию РОУ. Показатель увеличивается прямо пропорционально доле ароматической части в составе органического соединения [Weishaar et al., 2003]. Коэффициент соотношения гумусовых веществ (ГВ) в РОВ определялся как соотношение значений поглощения при длинах волн 400 и 600 нм (E4 : E6). Состав гумусовых веществ принципиально разделяется на две фракции: фульвокислоты и гуминовые кислоты. Коэффициент среднего молекулярного веса (МВ) РОВ определялся по соотношению значений длин светопоглощения между 250 и 365 нм (E2 : E3). Данный показатель имеет прямую зависимость от показателя содержания ГВ [Thurman, 1985; Lou, Xie, 2006; Olefeldt, Roulet, 2012]. Всего было проанализировано 97 водных экстрактов.

Определение значений pH

Показатель pH торфа измерялся потенциометрическим методом (HANNA Instruments, Edge, США) в суспензии образца торфа. К высушенной в количестве 1 г навеске торфяного образца (в случае минерального образца навеска составляла 10 г) добавляли 25 мл дистиллированной воды. Полученная смесь взбалтывалась в шейкере в течение 1 часа и после отстаивалась в течение 1 часа. Измерения проводили при температуре 20±2°С. Электроды погружались в верхнюю часть суспензии. Показания рН-метра фиксировались после стабилизации. Всего было проанализировано 97 водных экстрактов.

Визуализация данных

Обработка и визуализация данных была проведена на языке программирования R [R Core Team, 2020] с применением пакета ‘ggplot2’ [Wickham, 2016] для визуализации данных, ‘analogue’ [Simpson, Oksanen, 2020] и ‘rioja’ [Juggins, 2017] для построения стратиграфического профиля торфяной залежи. Стратиграфические зоны торфяных кернов были выделены с применением кластерного анализа [Grimm, 1987] и модели ‘broken-stick’ [Bennett, 1996]. Для определения уровня корреляции между переменными был использован тест Спирмана с уровнем значимости p<0.01. Корреляция считалась значимой при значении коэффициента Спирмана больше r>0.7.

РЕЗУЛЬТАТЫ

На основе полученных данных были составлены графики и выявлены зависимости показателей между собой (Рис. 1-3).

 

Рис. 1. Ботанический состав и физико-химические показатели торфяного керна «Cond-1»

(С_орг – органический углерод (%), Запас – запас углерода (кг/м²); УВ – увлажнение (ступень увлажнения); БЗ – активное богатство и засоление почв (ступень богатства); РОУ – растворенный органический углерод (мг/л); SUVA – показатель ароматичности; МВ – коэффициент среднего молекулярного веса; ГВ – коэффициент соотношения гумусовых веществ; pH – водородный показатель; Влажность (%), Степень разложения (%).

Fig. 1. Macrofossil and physicochemical parameters of the peat core "Cond-1"

(С_орг – organic carbon (%), Запас – carbon reserve (kg/m²); УВ – moisture (moisture level); БЗ – active soil richness and salinization (richness level); РОУ – dissolved organic carbon (mg/l); SUVA – aromaticity index; МВ – average molecular weight coefficient; ГВ – humic substances ratio coefficient; pH – hydrogen index; Влажность (%) – humidity (%), Степень разложения (%) – decomposition degree (%).

 

Рис. 2. Ботанический состав и физико-химические показатели торфяного керна «Cond-2»

(С_орг – органический углерод (%), Запас – запас углерода (кг/м²); УВ – увлажнение (ступень увлажнения); БЗ – активное богатство и засоление почв (ступень богатства); РОУ – растворенный органический углерод (мг/л); SUVA – показатель ароматичности; МВ – коэффициент среднего молекулярного веса; ГВ – коэффициент соотношения гумусовых веществ; pH – водородный показатель; Влажность (%), Степень разложения (%).

Fig. 2. Macrofossil and physicochemical parameters of the peat core "Cond-2"

(С_орг – organic carbon (%), Запас – carbon reserve (kg/m²); УВ – moisture (moisture level); БЗ – active soil richness and salinization (richness level); РОУ – dissolved organic carbon (mg/l); SUVA – aromaticity index; МВ – average molecular weight coefficient; ГВ – humic substances ratio coefficient; pH – hydrogen index; Влажность (%) – humidity (%), Степень разложения (%) – decomposition degree (%).

 

Рис. 3. Ботанический состав и физико-химические показатели торфяного керна «Cond-3»

(С_орг – органический углерод (%), Запас – запас углерода (кг/м²); УВ – увлажнение (ступень увлажнения); БЗ – активное богатство и засоление почв (ступень богатства); РОУ – растворенный органический углерод (мг/л); SUVA – показатель ароматичности; МВ – коэффициент среднего молекулярного веса; ГВ – коэффициент соотношения гумусовых веществ; pH – водородный показатель; Влажность (%), Степень разложения (%).

Fig. 3. Macrofossil and physicochemical parameters of the peat core "Cond-3"

(С_орг – organic carbon (%), Запас – carbon reserve (kg/m²); УВ – moisture (moisture level); БЗ – active soil richness and salinization (richness level); РОУ – dissolved organic carbon (mg/l); SUVA – aromaticity index; МВ – average molecular weight coefficient; ГВ – humic substances ratio coefficient; pH – hydrogen index; Влажность (%) – humidity (%), Степень разложения (%) – decomposition degree (%).

 

Для первого керна было определено 5 зон (Рис. 1): зона 1 характеризовалась эвтрофной стадией с древесно-травяным типом торфа, где доминирующими видами являлись вахта, осоки (C. lasiocarpa и C. rostrata), хвощ, а также остатки древесных пород, представленных сосной, кедром и березой. В зоне 2 наблюдалась смена стадии – от эвтрофной к мезотрофной: в ботаническом составе преобладали древесные виды, преимущественно сосна, кедр и береза. В зоне 3 фиксировалась смена растительных сообществ: доминирующую роль приобретали кустарнички, осока (C. rostrata), пушица и шейхцерия. В зоне 4 наблюдался переход к олиготрофной стадии, характеризующийся снижением доли осоковых и появлением сфагнового торфа, включающего в состав Sph. balticum, Sph. angustifolium и Sph. divinum. Зона 5 характеризовалась резким переходом растительного покрова к современному состоянию – мочажинному сообществу с преобладанием Sph. balticum.

Степень разложения торфа имела неизменно низкие значения для поверхности, до глубины 120 см не превышала 10%, а затем увеличивалась по всей глубине, достигая пика в 60% на 290 см, и снижалась до 45-50% до минерального ложа.

В нижнем слое стратиграфического профиля были обнаружены остатки хвоща. Выше по разрезу, на глубине 300-220 см, присутствовала примесь шейхцерии и вахты, отмечалось постепенное нарастание содержания коры и древесины сосны, ели и берёзы. С 300 до 60 см возрастало содержание пушицы, при этом на глубине 150-130 см была обнаружена прослойка вересковых кустарничков. С глубины 220-180 отмечалось увеличение доли сфагновых мхов – Sph. angustifolium и Sph. divinum. На глубине 170-60 см в ботаническом составе доминировали Sph. balticum.

Статистические характеристики физико-химических показателей торфа для керна «Cond-1» представлены в таблице 1.

 

Таблица 1. Статистические характеристики по торфяному керну «Cond-1»: минимальное (Min), максимальное (Max), среднее значения (X) и стандартное отклонение (STD).

Table 1. Statistical characteristics of the peat core “Cond-1”: minimum (Min), maximum (Max), average (X) and standard deviation (STD).

Показатель	Min	Мax	X	STD
Концентрация органического углерода (%)	19.4	60.7	51.5	8.4
Запас углерода (кг/м²)	0.9	11	6.4	3
УВ (ст.)	82.8	95.4	91.1	2.7
БЗ (ст.)	1.6	6.5	3.5	1.7
РОУ (мг/л)	84.5	314.7	181	57.8
SUVA254	0.2	0.9	0.6	0.2
МВ	2.9	5.1	4.2	0.5
ГВ	2.2	7.3	3.8	1.3
pH	3.2	4.2	3.7	0.3

 

Содержание органического углерода в торфе изменялось от 19.4% до 60.7%, при этом минимальное значение было обнаружено на глубине 320-330 см. Ниже по горизонту происходило увеличение содержания углерода, и к придонному слою концентрация составляла 38.9%. Максимальное значение (60.7%) наблюдалось в середине торфяной колонки (240-250 см). Среднее значение по всей глубине составило 51.5±8.4%.

Запас органического углерода в десятисантиметровом сегменте керна возрастал ближе к середине. Минимальный запас С был обнаружен на глубине 60 см – 0.9 кг/м²; максимальный запас – на глубине 210-220 см и составлял 11 кг/м²; среднее значение в десятисантиметровом сегменте составило 6.4±3 кг/м², а общий запас содержания углерода по всей глубине (до 400 см) составил 225 кг/м².

Для показателя УВ характерно неравномерное распределение значений по глубине – максимальное увлажнение наблюдалось на 240-250 см – 95.4 ступени (ст.), что характерно для болотной среды обитания. Минимальное значение 82.8 ст. (сыро-луговая среда обитания) было обнаружено на глубине 290-300 см. Среднее значение показателя УВ по всей торфяной колонке составило 91.1±2.7 ст.

Показатель БЗ имел тенденцию повышения значений к придонному слою. Минимальное значение 1.6 ст. (олиготрофная среда обитания) было обнаружено на глубине 50-60 см, максимальное значение в 6.5 ст. (мезотрофная среда обитания) наблюдалось на 360-370 см. Среднее значение показателя для всего керна составило 3.5±1.7 ст.

Содержание РОУ волнообразно изменялось по всему торфяному керну. Наибольшая амплитуда значений была обнаружена в пределах глубин 100-200 см. Максимальные значения обнаружены на глубине 190-200 см и составляли 315 мг/л. Минимальные значения наблюдались в поверхностных и придонных слоях (на глубине 60-120 см и 300-400 см), где концентрация изменялась в пределах 110÷200 мг/л.

Показатель SUVA₂₅₄ с глубиной увеличивался, где максимальный пик (0.9) зафиксирован на глубинах 160-180 см. Минимальное значение 0.2 наблюдалось на глубине 150-160 см. На глубинах 280-400 см было отмечено повышение значений показателя с 0.5 до 0.8±0.1.

Распределение значений МВ органических соединений имело резкие скачки от минимальных до максимальных значений (2.9÷4.9) на глубинах 200-300 см. Средние значения МВ на всей глубине торфяной залежи составили 4.2±0.5.

Значения ГВ в торфах были распределены неравномерно. Максимальные значения обнаружены на глубинах 270-340 см и изменялись в диапазоне 3.9÷4.2, а минимальные значения (3.2÷3.8) характерны для глубин 100-230 см.

Показатели pH и влажности имели взаимно обратную тенденцию: значение pH с глубиной увеличивалось, достигая максимального значения 4.2 на глубине 340 см и снижаясь до минимальных значений 3.8 в придонном слое; при этом значение влажности имело тенденцию к уменьшению с 95% на поверхности до 90% в придонном слое с пиковым снижением до 70% на 340 см.

Для второго керна было выделено 7 зон (Рис. 2): зона 1 характеризовалась эвтрофной стадией с древесно-травяным типом торфа, в котором преобладали осоки (C. lasiocarpa и C. rostrata), пушица, а также остатки древесных растений – ива, береза и сосна. В зоне 2 наблюдался переход к мезоолиготрофной стадии – в ботаническом составе осоки сменились на пушицу и кустарнички с примесью сфагновых мхов (Sph. divinum и Sph. balticum). Зона 3 характеризовалась переходом к олиготрофному типу питания и снижению доли пушицы, увеличению доли кустарничков и древесной растительности с заменой Sph. divinum на Sph. balticum. В зоне 4 были обнаружены пики пушицы и сфагнового мха (Sph. fuscum). Зоны 5 и 7 характеризовались преобладанием Sph. balticum в составе растительности с примесью Sph. divinum в зоне 6.

Степень разложения имела низкие значения для верхних 110 см, достигая 15%. С глубиной наблюдалось повышение до уровня в 30% на глубине 250 см. Максимальное значение зафиксировано в придонном слое – 70%.

Для торфяного керна «Cond-2» на глубине 320 см была обнаружена минеральная примесь. С глубины 320 см до 280 см в составе торфа увеличивалась доля осок (C. rostrata и C. lasiocarpa). Доля коры и древесины сосны доминировала в промежутке с 330 до 170 см, при этом остатки березы и ивы появлялись на глубине 330 см и пропадали на 240 см. Единично наблюдалась шейхцерия на глубинах 260-250 см, 200 см, 160-140 см и 90 см. С 200 см до 60 см шло преобладание Sph. balticum с небольшой примесью Sph. divinum на глубинах 290 см, 240 см, 200-180 см и 80-60 см и единичная встречаемость Sph. fuscum на 280 и 130 см. С глубины 330 см до глубины 60 см увеличивалось содержание пушицы и кустарничков.

Статистические характеристики физико-химических показателей торфа для керна «Cond-2» представлены в таблице 2.

 

Таблица 2. Статистические характеристики по торфяному керну «Cond-2»: минимальное (Min), максимальное (Max), среднее значения (X) и стандартное отклонение (STD).

Table 2. Statistical characteristics of the peat core “Cond-2”: minimum (Min), maximum (Max), average (X) and standard deviation (STD).

Показатель	Мin	Мах	X	STD
Концентрация органического углерода (%)	0.6	59.7	46.6	18.3
Запас углерода (кг/м²)	0.04	11.7	7	3.6
УВ (ст.)	86	94.3	91.1	1.9
БЗ (ст.)	1.7	4.6	2.6	0.9
РОУ (мг/л)	114.3	454.3	265.9	76.1
SUVA254	0.1	0.8	0.5	0.2
МВ	4.2	5.8	4.9	0.4
ГВ	3.6	32	8	6.2
pH	3.3	3.7	3.4	0.1

 

Концентрация органического углерода в торфе увеличивалась с глубиной, принимая минимальные значения (46.4%) на границе с минеральным слоем на глубине 300-310 см, а максимальные значения 59.7% – на глубине 290 см. Среднее значение концентрации составило 46.6±10.2%.

Кривая запаса органического углерода по глубине имела волнообразную форму. Зафиксировано увеличение показателя с глубины 60 см до 310 см, после чего значение уменьшалось. Минимальное значение (0.04 кг/м²) наблюдалось в придонном слое на глубине 330 см. Наибольшее значение наблюдалось на глубинах 120-150 см с показателем 11.7 кг/м². Среднее значение в десятисантиметровом сегменте составило 7.0±2.8 кг/м². Запас органического углерода по всей глубине керна (до 330 см) составил 209 кг/м².

УВ характеризовался высоким значением в верхних слоях – на глубине 90 см составил 94.3 ст. (болотная среда обитания); минимальное значение (86 ст.) обнаружено на глубине 270 см (сыро-луговая среда обитания). Среднее значение показателя УВ по всей глубине составило 91.1±1.9 ст.

Для БЗ было характерно увеличение значений с глубиной залегания торфяной залежи. Минимальные показатели наблюдались для поверхностных слоев 60-170 см с минимумом в 1.7 ст. (олиготрофная среда обитания) на глубине 110 см. Максимальные значения 4.6 ст. (мезоолиготрофная среда обитания) были обнаружены на глубинах 290-310 см. Среднее значение для всего керна составило 2.6±1.9 ст.

Содержание РОУ увеличивалось с глубиной залегания торфа. Максимальные значения концентраций (266÷454 мг/л) были обнаружены на глубинах 200-320 см. Минимальные значения концентрации наблюдались в поверхностных слоях торфяной залежи на глубине 60-100 см и составили 114÷243 мг/л.

Показатель SUVA₂₅₄ имел максимальные значения 0.38÷0.82 на глубинах 100-200 см. Снижение значений приходилось на глубины 210-350 см, и минимальные значения 0.2±0.1 были обнаружены на глубинах 320-330 см.

Значения МВ органических соединений имели максимальный пик 5.81 на глубине 210 см, а минимальные значения 4.2±0.3 – на глубинах 140 и 200 см. Средние значения МВ по всей глубине торфяной залежи составили 5.0±0.4.

Значения ГВ в торфе распределялись относительно равномерно с глубиной. Максимальный пик 32.0 наблюдался на глубине 210 см, а минимальные значения составили 3.6÷4.0.

Значение pH увеличивалось с глубиной с минимальных значений 3.25-3.40 у поверхности до максимальных значений 3.5±0.2 на глубинах 250-300 см. Показатель влажности имел тенденцию на снижение с глубиной залегания с 93% до 82%, принимая максимальные значения 91-92% на глубинах 60-100 см, а минимальные 82% – на глубине 310 см.

Для третьего керна выделено 6 зон (Рис. 3). В зоне 1 преобладала древесно-травяная группа торфа, характеризующаяся значительным содержанием пушицы, осок (C. lasiocarpa и C.rostrata), кустарничков, а также остатков древесных растений, представленных ивой и березой. В зоне 2 происходила смена стадий – эвтрофная сменилась на мезоолиготрофную, к перечисленным растительным остаткам добавились Sph. fuscum и остатки сосны. Зоны 3 и 4 отличались наличием древесно-травяной группы торфа (пушица, кустарнички, береза и сосна) с преобладанием в составе сфагновых остатков, включающих Sph. balticum, Sph. majus и Sph. divinum. В зонах 5 и 6 наблюдалось незначительное снижение доли пушицы и кустарничков при одновременном увеличении сфагновых мхов.

Степень разложения не превышала 15% до глубины 160 см с максимальным значением 80% в минеральном слое.

 На третьей точке минеральный слой начинался с глубины 400 см. На глубине 350 см увеличивалась доля коры и древесины ивы с примесью пушицы, осоки (C. lasiocarpa, C. rostrata) и Sph. fuscum. На глубинах 290-140 см присутствовали вересковые кустарнички, а на глубине 200 см в ботаническом составе появились кора и древесина березы и сосны. В диапазоне глубин 220-60 см в ботаническом составе преобладали Sph. balticum с примесью Sph. majus, Sph. divinum, кустарничками и пушицей.

Статистические характеристики физико-химических показателей торфа для керна «Cond-3» представлены в таблице 3.

 

Таблица 3. Статистические характеристики по торфяному керну «Cond-3»: минимальное (Min), максимальное (Max), среднее значения (X) и стандартное отклонение (STD).

Table 3. Statistical characteristics of the peat core “Cond-3”: minimum (Min), maximum (Max), average (X) and standard deviation (STD).

Показатель	Min	Мax	X	STD
Концентрация органического углерода (%)	45.3	56.2	52.3	2.7
Запас углерода (кг/м²)	2.3	8.2	5.2	1.3
УВ (ст.)	88.4	95.2	92.2	1.7
БЗ (ст.)	1.8	3.3	2.4	0.4
РОУ (мг/л)	142.7	365.6	276.9	69.9
SUVA254	0.4	0.6	0.5	0.1
МВ	3.6	6.8	5.5	0.5
ГВ	2.1	27	9.6	6.5
pH	3.4	4	3.7	0.1

 

Минимальные значения содержания органического углерода составили 45.3% на глубине 400 см; максимальные значения 56.2% были обнаружены на глубине 230 см. Среднее значение содержания органического углерода по глубине составило 52.3±2.7%.

Максимальный запас углерода (8.2 кг/м²) установлен на глубине 170 см. Минимальное значение (2.3 кг/м²) наблюдалось на глубине 100 см. Среднее значение запаса углерода в десятисантиметровом сегменте керна составило 5.2±1.3 кг/м². Общий запас углерода по всему торфяному керну составил 182 кг/м².

Значение УВ уменьшалось с глубиной. Среднее значение УВ составило 92.2±1.7 ст.; минимальное значение (88.3 ст.) наблюдалось на глубине 320 см (сыро-луговая среда обитания). Максимальное значение 95.2 ст. (болотная среда обитания) обнаружено на глубине 90 см.

Увеличение БЗ происходило с глубиной залегания торфяного керна. Среднее значение для всего керна составило 2.3±0.4 ст. Минимальное значение 1.7 ст. (олиготрофная среда обитания) наблюдалось на глубине 70 см; максимальное значение 3.2 ст. (олиготрофная среда обитания) – на глубине 370 см.

Максимальные концентрации РОУ были зафиксированы в пределах 190-390 см с пиком в 365.6 мг/л на глубине 310 см. Низкие концентрации характерны для верхних 180 см, где значения изменялись в пределах 142.8÷274.8 мг/л.

Показатель SUVA₂₅₄ имел максимальные значения в диапазоне 0.6÷0.7 с пиком на глубине 280 см. Каждые 100 см отмечались скачки значений. Минимальное значение 0.4 находилось на глубине 220-230 см.

Значения МВ увеличивались в диапазоне 4.68÷6.8 на глубинах 200-270 см с максимальным пиком на 230 см и минимальным значением на глубине 90 см. Среднее значение коэффициента по всей глубине торфяной залежи составило 5.6±0.5.

Значения ГВ на глубине 110 см имели высокое значение 27.0; увеличение показателя до 12.0÷12.8 было обнаружено на глубинах 140 см и 200-210 см соответственно. С глубины 210 см и до минерального дна значения коэффициента выходили за пределы обнаружения.

Значение показателя pH уменьшалось от 3.73 на поверхности до 3.44 на глубине 270 см с последующим ростом до 3.95 у минерального дна.

Максимальное значение показателя влажности торфа 93.8% было обнаружено на глубине 70 см с последующим снижением до минимального значения (88.8%) на глубине 360 см.

ДИСКУССИЯ

На особенности физико-химических свойств торфяных почв верхового болота «Кондинских озер» существенное влияние оказывает ботанический состав торфа. Эвтрофная стадия в развитии болотного массива характеризуется обилием остатков осок (C. lasiocarpa и C. rostrata), вахты и древесины (кора и древесина сосны, ели, березы и ивы). На исследуемом участке данный тип торфа формирует древесно-травяной слой мощностью до 1 метра. Эвтрофный горизонт на всех изученных участках входит в одну стратиграфическую зону. Мезотрофная стадия также занимает около 1 метра и подразделена на две стратиграфические зоны в Cond-1, а в Cond-2-3 – на одну стратиграфическую зону. Эти стадии сменяются олиготрофной стадией, для которой характерно преобладание сфагновых мхов, формирующих комплексный торф (Sph. balticum, Sph. majus, Sph. divinum и Sph. fuscum). Определена высокая положительная корреляция содержания вахты с неопределенными травяными остатками (r = 0.61, p<0.01).

По всем кернам сфагновые мхи в верхнем слое характеризуются меньшей степенью разложения, что отражается на свойствах органического вещества: в верхних слоях концентрация органического углерода и его содержание (в слое 10 см) имели средние значения 47.9% и 3.6 кг/м². Постепенное повышение значений содержания органического углерода 55.1% и его запаса (в слое 10 см) 7.5 кг/м² наблюдалось к середине торфяной залежи (с 160 до 310 см). Увеличение концентрации органического углерода напрямую связано с повышением степени разложения, а также с содержанием пушицы и кустарничков в составе торфа. Наличие в придонном слое минеральной примеси отображается в малом содержании органического углерода – 38.8%; запас органического углерода в слое 10 см составил 3.8 кг/м². Среднее значение концентрации по глубине всех трех кернов составило 50±9.8%, запас органического углерода в слое 10 см составил 5 кг/м².

Таким образом, торфы территории природного парка отличаются большими запасами органического углерода по сравнению с торфами верхового болота территории Мухрино, на что непосредственное влияние оказывает состав торфа. Так, наиболее распространенным видом торфа для территории парка «Кондинские озера» является древесно-травяная группа (около 60% от всех исследованных образцов). Для аналогичной подзоны тайги на верховом болоте Мухрино наиболее преобладающей группой торфа является сфагновый (22.5% от всех исследованных образцов) [Zarov et al., 2023].

Показатели УВ и БЗ имеют обратно пропорциональную зависимость. От верхнего горизонта до средних глубин торфяной залежи идет увеличение УВ до 91.9 ст. с дальнейшим уменьшением значения до 89 ст. Показатель БЗ увеличивается с глубиной до 4.4 ст. при среднем значении в 2.8 ст., что связано с изменением трофности растений-торфообразователей и доступностью минерального питания на дне. Определена высокая положительная корреляция между SUVA₂₅₄ и БЗ (коэффициент корреляции 0.5, p<0.01).

Для исследуемого участка на территории парка «Кондинские озера» концентрация РОУ увеличивается с глубиной, похожая тенденция прослеживается для верхового болота Мухрино [Shanyova, 2023]. Количественный состав РОУ зависит от содержания органического вещества в торфе и от степени разложения торфа [Kalbitz et al., 2002]. Основными факторами, которые определяют концентрацию РОУ, являются ботанический состав и тип торфа [Dalva M. et al., 1991]. Высокие концентрации РОУ совпадали с понижением показателя УВ и увеличением содержания пушицы, кустарничков и шейхцерии в стратиграфическом профиле. Снижение концентраций РОУ происходит с преобладанием остатков Sph. balticum, Sph. divinum в торфе. Среднее содержание РОУ для исследуемой торфяной залежи имеет значения в пределах 241.27 мг/л. Выявлена корреляция между содержанием пушицы и РОУ с коэффициентом r = 0.6, p<0.01.

Взаимосвязь между концентрациями РОУ в речных водах и запасами органического углерода в почве была исследована в различных пространственных масштабах системы реки Ди (Шотландия). В данном исследовании рассматривается возможность прогнозирования потока РОУ при условии наличия детализированных данных о почве [Aitkenhead et al., 1999]. Определено, что не только биологические процессы контролируют концентрацию РОУ, но и значения pH имеют значительное влияние на скорость высвобождения РОУ из торфа по причине влияния на сорбцию и растворимость органического углерода. Например, в болотах Великобритании низкие концентрации РОУ взаимосвязаны с подкислением торфяных вод до значения pH 3.5÷4.1 [Clark et. al., 2005].

По всем торфяным кернам значения pH увеличивались с глубиной залегания торфяной залежи. Увеличение значения pH предположительно связано с преобладанием кустарничков и древесных остатков в составе торфа. Определена корреляция pH c содержанием сосны (коэффициент r = 0.6, p<0.01) и вахты (коэффициент r = 0.6, p<0.01). В данном исследовании показана прямая зависимость между pH и степенью разложения торфяной залежи. На значение показателя pH влияет смешивание дождевой воды с торфяными водами. Это приводит к нарушению баланса растворимости органического углерода и, как следствие, к увеличению или уменьшению содержания РОУ в воде [Scott et al., 1998].

На болотах Швеции различия в химическом составе и оптических свойствах торфов определяет преобладание осок – концентрация РОУ и показатель SUVA₂₅₄ становятся ниже, чем на болотах с преобладанием сфагнового мха [Olefeldt, Roulet, 2012]. В исследовании на Скандинавских островах в водно-болотных угодьях уровень SUVA₂₅₄ увеличивался менее чем на 1 единицу за 3 метра. На глубине от 50-90 см SUVA₂₅₄ составлял 2.5, а затем увеличивался до 4.5 у поверхности, на глубине 6-12 см. При этом содержание РОУ и показатель SUVA₂₅₄ зависят от ботанического состава, в частности, наличие древесных остатков снижает концентрацию РОУ и имеет стабильное значение для показателя SUVA₂₅₄ (~3.7) [Agren et al., 2008].

Полученные нами данные по показателю SUVA₂₅₄ оказались ниже в 4 раза и имели тенденцию на увеличение к середине торфяной залежи, при этом на аналогичном верховом болоте Мухрино показатель SUVA₂₅₄ у поверхностных вод оказался значительно выше, изменяясь в пределах 2.0÷3.8 [Rakhova, 2024]. Обнаружено резкое снижение показателей SUVA₂₅₄ и РОУ с преобладанием в торфяной залежи Sph. fuscum и Sph. balticum. Придонные слои торфа имеют наименьшие значения SUVA₂₅₄ в связи с присутствием минеральной примеси. Кривая распределения показателя SUVA₂₅₄ по глубине характеризуется наличием пикового значения через каждые ~100 см. Определена высокая положительная корреляция SUVA₂₅₄ c содержанием вахты (показатель r = 0.6, p<0.01).

Соотношение различных фракций гуминовых веществ в торфяной залежи «Кондинских озер» распределено неравномерно, что связано с количеством органических веществ: лигнины, углеводы, липиды, полисахариды, а также содержание карбоксильных и фенольных оксигрупп [Thurman, 1985; Bolin et al., 1979]. Низкое значение показателя отражает высокую молекулярную массу гуминовых веществ в РОВ и степень гумификации [Traversa et al., 2014]. Гуминовые кислоты состоят из высокогумифицированных органических соединений, которые устойчивы к разложению и, следовательно, имеют тенденцию к накоплению [Machado et al., 2020].

Исследование, выполненное на северо-востоке Германии, подчеркивает значимую роль микроорганизмов в разложении и гумификации растительных остатков. Таким образом, возрастающая степень гумификации органического вещества связана с возрастающей микробной трансформацией водорастворимых органических веществ углерода [Kalbitz et al., 2002]. В нашем исследовании средний молекулярный вес имеет средние значения по глубине 4.2÷5.6, что, вероятно, взаимосвязано с содержанием Sph. balticum и Sph. majus. Выявлена корреляция среднего молекулярного веса с гуминовыми веществами (коэффициент r = 0.7, p<0.01).

Таким образом, обнаруженная динамика изучаемых показателей обусловлена специфическими условиями исследуемой территории, что требует проведения ряда дополнительных исследований.

ВЫВОДЫ

В данной работе были описаны три стратиграфических профиля, определен показатель УВ, значение БЗ, распределение концентрации РОУ, определены спектральные характеристики РОВ и содержание органического углерода в торфе на ненарушенных участках олиготрофного болота территории парка «Кондинские озера». Первая стадия развития болота характеризовалась преобладанием эвтрофной растительности, отлагающей древесно-травяной (шейхцерия, вахта и осоки) вид торфа, с долей гипновых мхов, встречающихся нерегулярно в озерах и местах выхода грунтовых вод. Обеднение питающих вод способствует переходу болот к мезотрофному этапу со сменой древесно-травяных видов торфа на травяные (пушица). Переход в олиготрофную стадию обусловлен минеральным питанием лишь за счёт атмосферных осадков с малым содержанием минеральных веществ, что приводит к развитию сфагновых видов мхов, отлагающих одноименный вид торфа. В целом болото на территории природного парка «Кондинские озера» характеризуется как равномерно развивающееся и захватывающее все стадии формирования болот.

 На основе выполненного исследования установлено, что содержание органического углерода и РОУ в торфах не всегда увеличивается с глубиной залегания. Различия в характере РОУ между торфяными колонками были обусловлены ботаническим составом торфа. Установлено, что для показателя SUVA₂₅₄ характерна прямая зависимость от содержания древесных остатков и кустарничков в торфе. Также высокие концентрации РОУ совпадают с понижением показателя УВ и увеличением содержания пушицы, кустарничков и шейхцерии в стратиграфическом профиле.

Отмечается корреляция соотношения гуминовых веществ и среднего молекулярного веса. Значимыми факторами являются состав ГВ – более сложные структуры могут характеризоваться более высоким средним молекулярным весом и устойчивость фракции – гуминовые вещества, образующиеся в процессе разложения органических материалов, которые менее подвержены деструкции и фотодеградации.

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

БЛАГОДАРНОСТИ

Работа выполнена при поддержке гранта Правительства Тюменской области в соответствии с программой Западно-Сибирского межрегионального научно-образовательного центра мирового уровня в рамках национального проекта «Наука».

Выражаем искреннюю благодарность за помощь в проведении исследований доктору биологических наук Елене Лапшиной и старшему научному сотруднику Евгению Зарову, а также лаборантам Артему Кулику и Леониду Литвинову.

Об авторах

V. S. Shanyova

ФГБОУ ВО «Югорский государственный университет»

Автор, ответственный за переписку.
Email: SHANYOVA.VIKA@mail.ru
Россия, г. Ханты-Мансийск

V. R. Batrshina

ФГБОУ ВО «Югорский государственный университет»

Email: lera.batrshina@icloud.com
Россия, г. Ханты-Мансийск

S. E. Rakhova

ФГБОУ ВО «Югорский государственный университет»

Email: sonya.rahova00@mail.ru
Россия, г. Ханты-Мансийск

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

Дополнительные файлы