Россия
Россия
УДК 51-7 Применение математических методов исследования в других областях знания
В статье представлены результаты моделирования переноса 210Pb вниз по торфяному керну для естественного болотного массива Иласский. В керне послойно исследованы активности радионуклидов 210Pb и 137Cs. Из-за неэкспоненциального снижения активности 210Pb возможны значительные ошибки в определении возраста и скоростей накопления торфа при применении классических моделей датирования. В статье приведены данные 210Pb-датирования, полученные по модели IP-CRS. Анализ данных показал, что вертикальное распределение 137Cs и 210Pb связано с параметрами торфа: зольностью и насыпной плотностью. Датирование керна по 210Pb показало для глубины 17—19 см возраст 1963-й год, что согласуется с данными 137Cs. Датировку данного горизонта косвенно подтверждают данные зольности и насыпной плотности торфа. Перечисленное подтверждает правильность и адекватность выбранной модели датирования естественного торфяника. Пик содержания антропогенных радионуклидов в торфяных отложениях соответствует 1963-му году — дате подписания Договора о частичном запрещении ядерных испытаний, который является ориентиром для геохронологических исследований. Приведённые в исследовании результаты линейной скорости накопления, скорости накопления массы торфяной залежи и величина атмосферного потока 210Pb хорошо согласуются с данными по торфяникам Северной Европы и Европейской Cубарктики России.
торфяной керн, скорость торфонакопления, 210Pb, 137Cs, миграция радионуклидов, моделирование IP-CRS, Иласский болотный массив
Омбротрофные торфяники являются уникальными экологическими архивами, изучение которых дает ценную информацию для понимания прошлых экологических и климатических событий, что особенно важно в связи с растущим интересом к глобальному изменению климата (Ndiaya et al., 2022; Sun et al., 2020). Дополнительно торфяники содержат данные о накоплении широкого спектра атмосферных загрязнителей, таких как тяжелые металлы и органические соединения (Magiera et al., 2021). Отметим, что с началом атомной эры на планете торфяники также стали концентрировать в себе техногенные радионуклиды и изотопы, выброшенные в атмосферу в результате испытаний ядерного оружия и радиационных аварий (Fialkiewicz-Koziel et al., 2014). Изучение последних в торфяниках является важнейшим инструментом для понимания хронологии торфа, развития и продуктивности торфяного месторождения во времени, изучения изменений климата и т. д. (Zhao et al., 2021). Основным радионуклидом атмосферных осадков, используемым в качестве изотопного маркера при изучении торфяников, является 137Cs. Однако из-за относительно высокой биодоступности 137Cs постоянно перемещается к растущим побегам сфагнового мха, где и накапливается (Mroz et al., 2017). В связи с этим восстановление точной хронологии торфа с использованием только 137Cs может быть затруднено. В качестве независимого параллельного геохронометра используется естественный радионуклид 210Pb (Fialkiewicz-Koziel et al., 2020). Предполагается, что скорость его поступления на поверхность торфяника относительно постоянна, а миграционная способность ограничена в силу химических свойств. В связи с этим в торфяной залежи, как правило, наблюдается экспоненциальное снижение активности 210Pb вниз по профилю, что даёт достаточно точную и непрерывную хронологическую последовательность торфа. В последние годы в некоторых исследованиях отмечается отклонение от типичного экспоненциального снижения, что приводит к ошибкам в оценке возраста и скорости накопления торфа (Cwanek et al., 2021). Возможной причиной увеличения миграции 210Pb является изменение комплекса физико-химических параметров торфяной залежи, которые зависят от конкретных природно-климатических условий местности и, в свою очередь, подлежат оценке. Несмотря на то, что торфяно-болотные экосистемы занимают огромные площади субарктического региона европейской части России, исследования радионуклидов в торфяниках этих территорий ограничены, поэтому, чтобы устранить недостаток знаний в области радиоэкологии торфяных болот для данного региона, нами проведено моделирование нисходящего переноса 210Pb на примере естественного торфяника — болотного массива Иласский.
Методы исследования
Торфяной керн ИСНО-1 (64°18¢55.3²N 40°41¢15.6²E) был отобран 12 августа 2020 года в пределах Иласского болотного массива, который является типичным представителем верховых болот северной тайги. Массив расположен в ~20 км к югу от Архангельска. Мощность торфяной залежи в месте отбора проб составляет 2.2 м. Залежь состоит из сфагновых мхов с примесью пушицы с увеличением её доли в нижнем слое керна. Степень деградации торфа увеличивается с глубиной и составляет 5—20 %. Колонка до глубины 49 см была извлечена с помощью ПВХ-трубы и далее разделена в лаборатории на слои по 2 см (исключение — слой 0—3 см) и подготовлена для дальнейших аналитических процедур. Определение физико-химических параметров торфа проводили в лаборатории болотных экосистем, радиохимическую подготовку и измерение проб — в лаборатории экологической радиологии ФИЦ комплексного изучения Арктики им. Н. П. Лавёрова УрО РАН.
Определение зольности торфа
Навеску торфа массой 2—3 г взвешивали на аналитических весах (𝑚𝑠𝑎𝑚𝑝𝑙𝑒) DA-224C (BEL ENGINEERING SRL, Италия) с пределом допускаемой погрешности ± 0.5 мг. Навеску помещали в кварцевый тигель, предварительно доведенный до постоянной массы прокаливанием в муфельной печи при 900 °С. Далее последовательно прокаливали тигель с пробой торфа при 525 и 900 °С до постоянной массы, которую фиксировали на аналитических весах. Расчет потерь при прокаливании (%) проводили по формуле согласно ГОСТ-1988.
Определение насыпной плотности торфа
Насыпную плотность торфа оценивали гравиметрическим методом для воздушно-сухих образцов с фракционным составом 0.5—2 мм. Определение проводили путем измерения объёма образцов с помощью мерного цилиндра с последующей фиксацией их массы на аналитических весах DA-224C в 5 параллельных повторностях для каждой пробы (ГОСТ-1988).
Радиохимическое определение радионуклидов
Определение радионуклида 137Cs проводили с применением низкофонового гамма-спектрометра производства CANBERRA Packard (США) с коаксиальным полупроводниковым детектором GX2018 на основе кристалла Ge(Li) и программным обеспечением Genie-2000. Удельную активность 137Cs определили по линии гамма-излучения 661.66 с квантовым выходом 89.90 %, принадлежащей его дочернему радионуклиду — 137mBa. Минимально измеряемые активности 137Cs при экспозиции t = 18 000 с для геометрии «плоский сосуд» объёмом 0.1 литр составили 0.1 Бк (Рекомендация МВИ, 2010).
Для определения 210Pb образец торфа обрабатывали согласно методике (Бахур и др., 2013). Полученный счётный образец с радионуклидами 210Po и 210Pb измеряют через 10 часов после его приготовления для распада дочерних альфа- и бета- излучателей (212, 214, 216, 218Po; 210Bi и 210Bi) на альфа-бета-радиометре «Абелия».
Результаты и обсуждения
Процедура датирования
При 210Pb-датировании нами использовалась модель IP-CRS (Initial Penetration — Constant Rate of Supply), для расчёта которой мы следовали процедуре, подробно изложенной у авторов (Olid, 2016). Эта двухэтапная модель описывает распределение атмосферного 210Pb (210Pbxs) в торфе, учитывая как поступление 210Pb в накапливающуюся торфяную матрицу, так и первоначальную миграцию 210Pb через самые верхние слои торфа вниз. Чтобы описать распределение атмосферного 210Pb (210Pbxs) в профиле торфа ИСНО, мы разделили торфяную колонку на n слоев, представленных расстоянием z (см) сверху вниз. Каждый слой простирается от глубины zi−1 до глубины zi (табл. 1).
В модели IP-CRS, которая является модификацией модели CRS, распределение объёмной активности 210Pbxs (Бк/см3) описывается термином смешивания, который включает опадание подстилки, биотурбацию и криотурбацию внутри акротельма. Предполагается постоянный коэффициент диффузии и непрерывное перемешивание в керне. Кроме того, учитываются вариации объёмной активности 210Pbxs в зависимости от скорости аккреции и её мигрирование вниз и перераспределение на глубине. Авторы аргументируют рассмотрение только начальной подвижности 210Pbxs из-за переноса жидкости тем, что Pb образует прочные комплексы с органическим веществом и переносится таким же образом, как и нерастворимые частицы после сорбции органическим веществом почвы. Наконец, учитываются также потери 210Pbxs при радиоактивном распаде (l) (Olid et al., 2016).
Ключевым моментом в модели IP-CRS является то, что количество 210Pbxs, которое было перемещено в объёме торфа, может быть оценено. Далее уже возможно применить модель CRS к скорректированному распределению 210Pbxs. Подробнее метод датирования и его валидность описаны у авторов (Olid et al., 2016).
Для расчёта хронологии IP-CRS мы выполнили следующие этапы:
Этап A. Ограничение глубины вымываемого отдела (zk ) торфяного разреза
Вертикальную миграцию 210Pbxs в профиле ИСНО количественно оценили, разделив его на две зоны. Зона вымывания от поверхности до глубины zk составила 20 см, здесь 210Pbxs мигрировал вниз, что моделируется предположением разных скоростей проникновения (ri) для каждого слоя. Толщина верхней вымываемой зоны варьируется и состоит из живых мхов и самых верхних слоев акротельма, который характеризуется высокой пористостью и гидравлической проводимостью, что способствует смыву вниз атмосферных отложений. Зона вмывания (накопления) — совокупный отдел от zk до ∞ (для керна ИСНО составила от 20 до 33 см), где слои торфа получали дополнительный запас 210Pbxs из самых верхних слоёв. Здесь слои акротельма имеют более высокую плотность и меньшую пористость. Уменьшение гидравлической проводимости, вероятно, приводит к образованию слоя, в котором накапливается вынесенный 210Pbxs. Нижняя граница зоны соответствует самой большой глубине с обнаруживаемой активностью 210Pb (Olid et al., 2016).
Этап В. Решение системы интегродифференциальных уравнений
Решили систему n интегродифференциальных уравнений (уравнения (1) и (2)) по разрезу торфа, выполнили соблюдение ограничивающих условий (уравнения (3) и (4)) и ограничений непрерывности между слоями (уравнения (5) и (6)) для получения скоростей проникновения (ri) и дробных коэффициентов усиления (fi). Здесь модель была реализована в программе Matlab. Вывод формул и расчёты подробнее указаны в дополнительных материалах к статье авторов (Olid et al., 2016).
Решения уравнений:
От 0 до zk:
(1)
где i =1, 2,…, k.
От zk до ∞:
(2)
где i = k + 1, k + 2, …., n.
Константы Ai и Bi определяются итеративно путем решения ограничивающих условий:
(3)
(4)
где F — атмосферное поступление 210Pb (Бк/м2·год) (F = l·I, где I — общий запас 210Pbxs (Бк/м2)), а ограничения непрерывности между слоями следующие:
(5)
(6)
Этап С. Определение запаса 210Pbxs , удаленного из каждого слоя зоны вымывания
После решения системы n оценивали запас 210Pbxs (Бк/м2), мигрировавшего из каждого слоя верхнего отдела согласно уравнению:
(7)
где Ti — время образования этого слоя.
Этап D. Определение общего запаса 10Pbxs
Количественно определили общий запас 210Pbxs, перенесенного вниз, суммируя 210Pbxs, полученный на предыдущем этапе.
Этап Е. Определение дополнительного 210Pbxs для каждого слоя
Количественно оценили запас 210Pbxs, включенного в качестве дополнительного источника в каждый слой, умножив общее удаленное количество 210Pbxs на соответствующее дробное усиление fi .
Этап F. Определение исправленных объёмных активностей 210Pbxs
1. Добавили объёмные активности 210Pbxs, полученные на этапе C, в соответствующие слои промываемой зоны; 2. Вычли объёмные активности 210Pbxs, полученные на этапе E, из соответствующих слоёв более глубокой зоны.
Этап G. Применение модели к исправленному профилю
Применили CRS-модель к скорректированному профилю объёмной активности 210Pbxs и дополнительно использовали метод Монте-Карло, подробно описанный в работах коллектива авторов (Sanchez-Cabeza et al., 2014; Яковлев и др., 2020).
Физико-химические параметры торфа
Величина показателя зольности для профиля ИСНО находится в интервале 1.10—5.77 %, что позволяет отнести его к малозольному типу (рис. 1). В вертикальном распределении в горизонте 20 см отмечается максимум содержания зольных компонентов, который может быть обусловлен как вымыванием минеральной части из вышележащих горизонтов, так и антропогенной запыленностью атмосферы в период формирования торфяной залежи. Ниже глубины 21 см наблюдается тенденция немонотонного снижения показателей с глубиной, связанная с атмосферным типом питания болота.
Согласно полученным результатам, насыпная плотность постепенно увеличивается с глубиной залегания и находится в интервале 0.075 ÷ 0.109 г/см3, что характерно для верховых торфяников в ненарушенном состоянии. Естественное влагосодержание находится в интервале 15 ÷ 24 г/г и снижается с глубиной, что отчасти является следствием постепенного уплотнения пористой структуры торфа.
Вертикальное распределение радионуклидов в торфяном профиле
Активность 137Cs в торфяном профиле ИСНО-1 варьирует в пределах от 3.2 до 45.6 Бк/кг. Согласно графику вертикального распределения (рис. 2), отмечается два выраженных пика активности 137Cs: первый пик с активностью 45.5 Бк/кг находится на глубине 19—21 см, второй пик с активностью 45.6 Бк/кг — на глубине 3—5 см. Отметим также высокие показатели 137Cs, которые составили 43.4 Бк/кг на глубине 17—19 см. Можно предположить миграцию радионуклида вниз по керну. В целом наибольшее распределение 137Cs в верхней части керна является типичным для верховых болот (Mroz et al., 2017).
Максимум активности 137Cs здесь связан отчасти с химическим сродством между Cs и K, которые совместно активно транспортируются по торфяному разрезу корнями растений. Другой причиной высокой подвижности 137Cs в верховых торфяниках указывается отсутствие подходящих минеральных частиц для адсорбции 137Cs (Fialkiewicz-Koziel et al., 2014). Однако мы наблюдаем помимо пика активности 137Cs в верхней части профиля сопоставимую активность на глубине от 17 до 21 см (45.6 Бк/кг в верхней части профиля, 43.4 Бк/кг на глубине 17—19 см и 45.5 Бк/кг на глубине 19—21 см соответственно). Схожая картина изменения активности цезия и содержания золы по глубине (рис. 1, 2), вероятно, связана с механизмом их поступления в торф с атмосферными осадками. Как можно видеть на рис. 2, максимальная глубина усваивания цезия растительностью составляет 11 см для исследуемого профиля, поскольку ниже, на глубине 19—21 см, мы наблюдаем пик активности, превышающий по значению максимальные верхние показания и обусловленный, вероятно, глобальными выпадениями.
Активность 210Pb в изученном торфяном профиле варьирует от 21.9 до 330.6 Бк/кг. Максимум активности 210Pb приходится на горизонт 3–5 см (рис. 2). Ниже горизонта 35—37 см активность 210Pb перестает изменяться и составляет во всех нижележащих горизонтах ~ 26 Бк/кг, т.к. здесь отсутствует избыточный атмосферный свинец. Наблюдаемая активность 210Pb в этих горизонтах является поддерживаемой за счет распада 226Ra. Вымывание 210Pb из верхних слоёв (0—3 см), вероятно, связано с атмосферными осадками (Olid et al., 2016).
Миграция 210Pb вниз по торфяному профилю происходит в коллоидной форме в виде комплексов с гумусовыми кислотами (фульво- и гуминовыми кислотами) и в ионной форме. Миграция 137Cs вниз по разрезу происходит только в ионной форме, а связывание 137Cs практически полностью обеспечивается ионно-обменным механизмом (Schleich et al., 2000).
Хронология и скорость накопления торфа
Результаты моделирования и датирования торфяного профиля ИСНО-1 представлены на рис. 3. и в табл. 1. Данные показали, что горизонт 17—19 см соответствует возрасту 1963-го года. Выделяемый пик активности 137Cs попадает на горизонт 19—21 см, это также совпадает с пиком показателя зольности. Как правило, пики антропогенных радионуклидов в торфяных отложениях указывают на конкретныe масштабныe события в истории радиоактивных выпадений, такие как Договор о частичном запрещении ядерных испытаний 1963 года и авария на Чернобыльской АЭС, обеспечивая ориентир для геохронологических исследований (Mroz et al., 2017; Putyrskaya et al., 2019). Однако 137Cs в связи с его высокой биодоступностью и, как следствие, миграцией в торфяных отложениях не всегда указывает на радиационные события в профиле торфа (Fialkiewicz-Koziel et al., 2014).
В нашем случае физико-химические условия торфяного профиля способствовали вымыванию до глубины 19—21 см части 137Cs, поступившего на поверхность торфяника в результате глобальных выпадений 1963 года. Полученные данные подтверждают правильность и адекватность выбранной нами модели датирования торфяного керна.
Результаты расчёта линейной скорости накопления s и скорости накопления массы r в изученной торфяной залежи представлены на рис. 4. Значения s варьировались от (0.073 ± 0.05) до (0.635 ± 0.05) см/год и составляли в среднем (0.28 ± 0.08) см/год, что несколько ниже ранее полученных данных для Европейской Субарктики России (Yakovlev et al., 2021) и ниже значений, указанных в работе Cwanek et al. (2021), где последние составили (0.38 ± 0.07) см/год. Скорость накопления массы торфа r находилась в диапазоне от (0.0019 ± 0.001) до (0.0199 ± 0.001) г/см2·год. Среднее значение r составило (0.0117 ± 0.001) г/см2·год.
В вертикальном распределении линейной скорости накопления s отмечается равномерное уменьшение с глубиной с выраженным пиком на глубине 19—21 см, что датируется 1964 годом и, возможно, объясняется изменением климатических условий произрастания торфа.
Картина распределения скорости накопления массы r имеет сложный характер с несколькими пиками в верхней части профиля и закономерным уменьшением с глубиной, поскольку там происходит уплотнение торфяной залежи.
Основываясь на датировании 210Pb, нами был оценен атмосферный поток 210Pb (58.98 ± 2.94) Бк/м2·год, который хорошо согласуется с данными по торфяникам Северной Европы. Атмосферный поток 210Pb для юга Финляндии составил 50—80 Бк/м2·год (Paatero et al., 1998), в то же время значение 210Pb для центральной Швеции составило 78 Бк/м2·год (Vassma et al., 2017). Дополнительно отметим, что полученные нами ранее данные находились в пределах 80—91 Бк/м2·год для Европейской Субарктики России (Yakovlev et al., 2021).
Заключение
Исследуемая часть профиля торфяной залежи сложена преимущественно сфагновым видом торфа. Физико-химические параметры изученной торфяной залежи находятся в интервалах, типичных для ненарушенных верховых торфяников Европейской Субарктики России. Установлено, что в разрезе выделяется горизонт на глубине 20 см, характеризующийся существенными колебаниям параметров зольности и насыпной плотности. Формирование этого горизонта может быть обусловлено как вымыванием минеральной части из вышележащих горизонтов, так и антропогенной запыленностью атмосферы в период формирования торфяной залежи.
Распределение 137Cs по глубине обнаруживает два выделяющихся максимума активности — в верхней части профиля и на глубине 19—21 см, обусловленных особенностями вертикальной миграции и динамикой поступления цезия с атмосферными выпадениями.
Для распределения 210Pb выявлено неэкспоненциальное снижение активности с глубиной. Наблюдаемое отклонение от типичного снижения 210Pb может приводить к серьёзным ошибкам в оценке возраста и скорости накопления торфа при датировании торфа классическими моделями, поэтому была использована разновидность моделирования CRS — IP-CRS, которая учитывает нисходящий перенос 210Pb в керне. Моделирование проведено на примере керна из естественного торфяника Иласского болотного массива Европейской Субарктики России.
Результаты датирования торфяного профиля моделью IP-CRS показали, что горизонт 17—19 см соответствует возрасту 1963-го года. Ориентиром для геохронологических исследований выступило подписание Договора о частичном запрещении ядерных испытаний в 1963 году. В то же время отметим, что пик активности 137Cs и максимум показателя зольности попадает на горизонт 19—21 см. В данном случае имеет место миграция 137Cs вниз по керну, т. к. его особенность — это высокая биодоступность в торфяных отложениях, поэтому он не всегда может указывать на радиационные события. Тем не менее результаты подтверждают правильность и адекватность выбранной нами модели датирования торфяного керна.
Результаты исследования позволяют расширить представления о миграционных процессах радионуклидов атмосферных выпадений в торфяниках Европейской Субарктики, что может послужить базисом для поддержания устойчивой экологической обстановки при освоении и использовании ресурсного потенциала торфяно-болотных экосистем.
1. Бахур А. Е., Мануилова Л. И., Зуев Д. М., Овсянникова Т. М., Трухина Т. П. Методика измерений удельной активности полония-210 (210Po) и свинца-210 (210Pb) в пробах почв, грунтов, донных отложений, горных пород и строительных материалов на их основе альфа-бета-радиометрическим методом с радиохимической подготовкой: Методика: ФР.1.40.2013.15381. М.: ВИМС, 2013. 17 с.
2. ГОСТ 27784-88. Почвы. Методы определения зольности торфяных и оторфованных горизонтов почв. М.: Изд-во стандартов, 1988. 7 с.
3. ГОСТ 28268-89. Почвы. Методы определения влажности, максимальной гигроскопической влажности и влажности устойчивого завядания растений. М.: Стандартинформ, 1989. 8 с.
4. Рекомендация МВИ. МH 3421-2010. Методика выполнения измерений объемной и удельной активности гамма-излучающих радионуклидов на гамма-спектрометрах с полупроводниковыми детекторами. Минск: Изд-во Белорусского гос. ин-та метрологии, БелГИМ, 2010.
5. Яковлев Е. Ю., Очеретенко А. А., Спиров Р. К., Мищенко Е. В., Жуковская Е. В. Применение метода моделирования по типу Монте-Карло для оценки погрешностей 210Pb-датирования торфяных отложений Европейской Субарктики России (на примере Архангельской области) // Вестник геонаук. 2020. № 7. С. 18-27. DOIhttps://doi.org/10.19110/geov.2020.7.3
6. Cwanek A., Łokas E., Mitchell E. A. D., Mazei Y., Gaca P., Milton J. A. Temporal variability of Pu signatures in a 210Pb-dated Sphagnum peat profile from the Northern Ural, Russian Federation // Chemosphere, 2021. No. 281. P. 130962. DOIhttps://doi.org/10.1016/j.chemosphere.2021.130962
7. Fialkiewicz-Koziel B., Kolaczek P., Piotrowska N., Michczynski A., Lokas E., Wachniew P., Woszczyk M., Sensula B. High-resolution age-depth model of a peat bog in Poland as an important basis for paleoenvironmental studies // Radiocarbon, 2014. No. 56, P. 109-125. DOIhttps://doi.org/10.2458/56.16467
8. Fialkiewicz-Koziel B., Lokas E., Galka M., Kolaczek P., De Vleeschouwer F., Le Roux G., Smieja-Krol B. Influence of transboundary transport of trace elements on mountain peat geochemistry (Sudetes, Central Europe) // Quaternary Science Reviews, 2020. No. 230, P. 106162. DOI:https://doi.org/10.1016/j.quascirev.2020.106162
9. Magiera T., Szuszkiewicz M. M., Michczyński A., Chróst L., Szuszkiewicz M. Peat bogs as archives of local ore mining and smelting activities over the centuries: A case study of Miasteczko Śląskie (Upper Silesia, Poland) // CATENA, 2021. No. 198, P. 105063. DOI:https://doi.org/10.1016/j.catena.2020.105063
10. Mroz T., Lokas E., Kocurek J., Gasiorek M. Atmospheric fallout radionuclides in peatland from Southern Poland // Journal of Environmental Radioactivity, 2017. No. 175-176, P. 25-33. DOI:https://doi.org/10.1016/j.jenvrad.2017.04.012
11. Ndiaye A., Bentaleb I., Favier C., Fourel F., Sebag D., Fall M., Giresse P., Diouf B. Reconstruction of the holocene climate and environmental changes of Niayes peat bog in northern coast of Senegal (NW Africa) based on stable isotopes and charcoals analysis. Quaternary Science Reviews, 2022. No. 289, P. 107609. DOI:https://doi.org/10.1016/j.quascirev.2022.107609
12. Olid C., Diego D., Garcia-Orellana J., Cortizas A. M., Klaminder J. Modeling the downward transport of 210Pb in Peatlands: Initial Penetration-Constant Rate of Supply (IP-CRS) model // Science of The Total Environment, 2016. No. 541. P. 1222-1231. DOI:https://doi.org/10.1016/j.scitotenv.2015.09.131
13. Paatero J., Jaakkola T., Kulmala S. Lichen (sp. Cladonia) as a deposition indicator for transuranium elements investigated with the Chernobyl fallout // Journal of Environmental Radioactivity, 1998. No. 38. P. 223-247. DOI:https://doi.org/10.1016/S0265-931X(97)00024-6
14. Putyrskaya V., Klemt E., Röllin S., Corcho-Alvarado J., Sahli H. Dating of recent sediments from Lago Maggiore and Lago di Lugano (Switzerland/Italy) using 137Cs and 210Pb // Journal of Environmental Radioactivity, 2020. No. 212. P. 106135. DOI:https://doi.org/10.1016/j.jenvrad.2019.106135
15. Sanchez-Cabeza J. A., Ruiz-Fernandez A. N., Ontiveros-Cuadras J. F., Perez Bernal L. H., Olid C. Monte Carlo uncertainty calculation of 210Pb chronologies and accumulation rates of sediments and peat bogs // Quaternary Geochronology, 2014. No. 23. P. 80-93. DOI:https://doi.org/10.1016/j.quageo.2014.06.002
16. Schleich B. N., Degering D., Unterricker S. Natural and artificial radionuclides in forest and bog soils: Tracers for migration processes and soil development // Radiochimica Acta, 2000. No. 88. P. 803-808. DOI:https://doi.org/10.1524/ract.2000.88.9-11.803
17. Sun A., Yang Y., Wu H., Ran M. Climate change on the northeastern Tibetan Plateau during the past~600 years inferred from peat pollen records // Review of Palaeobotany and Palynology, 2020. No. 276, P. 104194. DOI:https://doi.org/10.1016/j.revpalbo.2020.104194
18. Vaasma T., Karu H., Kiisk M., Pensa M., Isakar K., Realo E., Alliksaar T., Traczyk H. A. Pb-210 and fly ash particles in ombrotrophic peat bogs as indicators of industrial emissions // Journal of Environmental Radioactivity, 2017. No. 174. P. 78-86. DOIhttps://doi.org/10.1016/j.jenvrad.2016.07.027
19. Yakovlev E., Spirov R., Druzhinin S., Ocheretenko A., Druzhinina A., Mishchenko E., Zhukovskaya E. Atmospheric fallout of radionuclides in peat bogs in the Western Segment of the Russian Arctic // Environmental Science and Pollution Research, 2021, without a number. P. 1-11. DOI:https://doi.org/10.1007/s11356-020-12224-7
20. Zhao X., Hou X., Zhang D., Yang Y., Huang Z., Liu Q. Records of iodine isotopes (129I, 127I) in the Barkol peat bog from northwest China and their sources, transport and preservation // Chemosphere, 2021. No. 279, P. 130531. DOI:https://doi.org/10.1016/j.chemosphere.2021.130531
