CC BY
83
я
Физическая география и биогеография, география почв и геохимия ландшафтов
УДК 551.8, 550.4
СКОРОСТИ ОСАДКОНАКОПЛЕНИЯ И ВЫДЕЛЕНИЕ ГРАНИЦ ГОДОВЫХ СЛОЕВ ГОДИЧНО-СТРАТИФИЦИРОВАННЫХ ОТЛОЖЕНИЙ ПРИЛЕДНИКОВОГО ОЗЕРА КУЧЕРЛИНСКОЕ (АЛТАЙ) ПО ДАННЫМ СКАНИРУЮЩЕГО МИКРО-РЕНТГЕНОФЛУОРЕСЦЕНТНОГО АНАЛИЗА НА СИНХРОТРОННОМ
ИЗЛУЧЕНИИ
Б01: 10.24411/1728-323Х-2018-13070
Ф. А. Дарьин, м. н. с. Института географии РАН, Москва; аспирант Института ядерной физики им. Г. И. Будкера СО РАН, [email protected]. Новосибирск, О. Н. Соломина, академик РАН, директор Института географии РАН, Москва, А. М. Грачев, м. н. с. Института географии РАН, Москва,
А. В. Дарьин, с. н. с. Института геологии и минералогии им. В. С. Соболева СО РАН, Новосибирск,
Я. В. Ракшун, ученый секретарь Института ядерной физики им. Г. И. Будкера СО РАН, Новосибирск,
Д. С. Сороколетов, м.н.с. Института ядерной физики им. Г. И. Будкера СО РАН, Новосибирск
Рассмотрены методические вопросы проведения микроаналитических исследований донных осадков приледниковых озер, содержащих ежегодно ламинированные слои («ледниковые глины»). Из образцов донных осадков высокогорного ледникового озера Кучерлинское (Алтай) с визуально различимыми годовыми слоями были изготовлены твердые препараты, пропитанные эпоксидной смолой. Процедура подготовки образцов позволяла сохранить исходную структуру и состав осадка. Были получены данные об изменениях породообразующих и микроэлементов в годовом цикле осадкообразования. Сканирующий рентгенофлуоресцентный микроанализ образцов проводился на экспериментальной установке в ЦКП «Сибирский центр синхротронного и терагерцового излучения» (ИЯФ СО РАН, Новосибирск). Профили сканирования с шагом 100 мкм, содержащие данные о распределении более чем 20 элементов, были совмещены с фотографиями поверхности исследуемых образцов. Получен исходный экспериментальный материал для оценки скорости и выявления литолого-геохимических особенностей процесса осадконакопления, а также для выявления количественных характеристик связи состава и структуры осадков с данными инструментальных метеонаблюдений.
Methodological issues of conducting micro-analytical studies of the bottom sediments of the proglacial lakes, containing annually laminated layers ("glacial clays") are considered. Solid samples impregnated with epoxy resin were made from the samples of bottom sediments of high mountain periglacial lake Kucherla (the Altai) with visually discernible annual layers. The procedure for preparing the samples allowed us to preserve the original structure and composition of the sediment. The data were obtained on the changes of rock forming elements and microelements in the annual cycle of sedimentation. Scanning X-ray fluorescence microanalysis of the samples was carried out at the experimental facility at the Collective Use Center "Siberian Center for Synchrotron and Terahertz Radiation" (INP SB RAS, Novosibirsk). The scanning profiles with the step of 100 цш, containing data on the distribution of more than 20 elements, were combined with the photographs of the surface of the investigated samples. The initial experimental data for estimating the sedimentation rate and for uncovering the lithologic-geochemical features of the sedimentation process have been obtained. This material is also suitable for determining the quantitative characteristics of the relationship between the composition and structure of the sediment and the data of instrumental meteorological observations.
Ключевые слова: палеореконструкции, осадки, приледниковые озера, рентгено-флуоресцентный анализ, синхротронное излучение, осадконакопление, элементный состав.
Keywords: paleoreconstructions, sediments, periglacial lakes, X-ray fluorescence analysis, synchrotron radiation, sedimentation, elemental composition.
Важной проблемой в современных науках о Земле является реконструкция условий прошлого, как климатических, так и геохимических, а также биотических. Для этой цели используются различные природные архивы накапливающиеся непрерывным образом на больших временных промежутках. Одним из таких архивов являются озерные осадки.
Формирование донных осадков (как в прибрежных морях, так и в континентальных озерах) происходит в условиях меняющейся внешней среды. Изменения температуры и количества атмосферных осадков влияют на вещественный состав и структуру донных отложений. Изменения химического состава донных отложений можно проследить, исследуя керн осадков от поверхности (граница вода — осадок) вглубь отложений. Двигаясь от современности в прошлое, мы получаем информацию об изменениях природной среды с детальностью, зависящей от дискретности опробования. При традиционных м етодах исследования проводится нарезка керна на отдельные образцы с шагом 5—10 мм. Используя сканирующие методы, можно уменьшить шаг опробования до величин
-Ф-
-Ф-
10—100 мкм, что позволяет строить временные ряды литолого-геохимических данных с годовым (а в ряде случаев и с сезонным) временным разрешением.
Объектом исследования в данной работе являлось высокогорное озеро Кучерлинское (Алтай) с ежегодно слоистыми (варвными) донными отложениями. Термин «варвы» (уагуе8) впервые использовал шведский геолог [1] для описания слоистых приледниковых озерных осадков. В российской литературе используется термин «ленточные глины» — отложения приледниковых озер, состоящие из чередующихся слоев тонкозернистого песка, и глины — продуктов осаждения ледниковой мути. Варвы распространены в Скандинавии, Северо-Западной части России, Польше, Германии, США и Канаде. В литературе описания ламинированных осадков Азиатской части России встречаются редко. Рассмотрение распределения годично-слоистых отложений на карте мира, представленное в обзорной публикации [2], скорее отражает степень изученности территории, чем их реальное наличие. Необходимо отметить, что годично-слоистые осадочные отложения рассматриваются мировым палеоклиматическим сообществом как один из важнейших палеоархи-вов. В данном случае исследователи располагают точной хронологической привязкой слоев в сочетании с исключительно высоким временным разрешением.
Слоистость в осадке обусловлена неравномерным привносом обломочного материала в разные сезоны года (песчанистые — летом, глинистые — зимой). Варвы образуются благодаря сезонному поступлению взвешенных частиц, в случае при-ледниковых озер источниками являются водные потоки с ледников. Варвы состоят из крупных частиц (песка, ила) и мелкозернистой глины, осаждающейся в зимний период.
Использование синхротронного излучения для рентгенофлуоресцентного анализа (РФА-СИ) донных осадков имеет ряд преимуществ, по сравнению с другими методами. Сканирование в пучке СИ позволяет получить чрезвычайно подробную летопись по распределению породообразующих и микроэлементов в осадках. Первые предложения по использованию СИ в качестве чрезвычайно яркого источника рентгеновского излучения для физических исследований появились в 1970-х годах, тогда как до этого данное излучение рассматривалось исключительно как паразитное. Институт ядерной физики им. Г. И. Буд-кера СО РАН (ИЯФ СО РАН) по праву входит в число мировых лидеров в этой области: первые эксперименты с СИ здесь были проведены еще в 1973 году [3]. Работы по сканирующему микро-
анализу геологических (в том числе осадочных) образцов были начаты в конце 1980-х годов [4, 5] и продолжены в серии работ по исследованию озерных донных отложений [6—10].
В отличие от обычного лабораторного метода РФА с рентгеновской трубкой (ХИБ), метод РФА-СИ позволяет снизить предел обнаружения в 10—100 раз и определять в несколько раз больше элементов. Аналитические преимущества данного подхода относительно других методов рассмотрены в ряде работ [11—13]. Уникальные возможности РФА СИ обусловлены непрерывным спектром и малой угловой расходимостью СИ, его естественной поляризованностью и чрезвычайно высокой яркостью. При использовании нескольких оптимально подобранных режимов измерения можно определять практически всю таблицу Менделеева в сканирующем режиме с пространственным разрешением 0,1—10 мм. В настоящее время метод специально адаптирован для палео-климатических задач, разработаны и аттестованы методики подготовки и исследования кернов донных осадков [14, 15]. В последние годы развитие получил вариант метода с использованием фокусирующей рентгеновской оптики — поликапиллярных линз [16, 17]. При этом пространственное разрешение анализа достигает 10 мкм, что позволяет проводить детальные исследования отдельных годовых слоев в донных отложениях. На рис. 1 представлены измерительные возможности ЦКП «Сибирский центр синхротронного и тера-герцового излучения» (ИЯФ СО РАН, Новосибирск) по набору элементов в стандартных методиках.
Основной целью настоящего исследования являлось получение надежных аналитических данных о вариациях элементного состава донных осадков в годовом цикле осадконакопления. Эти данные необходимы для проведения дальнейших исследований, предполагающих выделение границ отдельных годовых слоев по геохимическим маркерам, подсчет слоев, определение их мощности.
Объект исследования. Озеро Кучерлинское расположено у подножия северного склона Ка-тунского хребта в верховьях реки Кучерла (рис. 2). Озеро находится на высоте 1790 м над уровнем моря. Длина озера 5 км, ширина — 1 км, максимальная глубина достигает 45 м. На рис. 2 представлено место расположения озера и точки отбора кернов на глубинах 45 и 30 м.
Озеро находится в труднодоступной местности и подвержено минимальному антропогенному воздействию, что делает его перспективным объектом для палеоклиматических реконструкций.
Отбор проб и подготовка образцов. Керны осадков озера Кучерлинское были получены с ис-
"Ф
-Ф-
пользованием ящика-пробоотборника (Ьох-согег) и ударных трубок в 2009 и 2018 гг. (табл. 1). Использование ящика-пробоотборника дало возможность получения сохранными верхних слоев керна вплоть до года отбора.
Керны после вскрытия в лабораторных условиях разрезались вдоль оси отбора на две половины. Одна продольная половина керна ис-
пользовалась для приготовления твердых блоков, зафиксированных в эпоксидной смоле. Размер каждого блока составлял 170 х 30 х 5 мм. Первый блок покрывает 0—170 мм от верха керна, второй — следующие 17 см с перекрытием и т. д. Блоки отбирались при помощи специальных алюминиевых форм, затем высушивались (лиофиль-ная сушка) и пропитывались смолой. Из блоков
1 н Породообразующие элементы 1-10% * | Породообраз. элементы 0,1-10% * Микроэлементы (до 10"б%) * Микроэлементы (до 10'4%) ** 2 Не
3 и 4 Ве * определяются по К-сериям ** определяются по ¿-сериям 5 В 6 с 7 N 8 О 9 Р 10 №
11 N3 12 Мд 13 А1 14 в! 15 Р 16 Б 17 С1 18 Аг
19 К 20 Са 21 Эс 22 Т\ 23 V 24 Сг 25 Мп 26 Ре 27 Со 28 № 29 Си 30 гп 31 ва 32 ве 33 Ав 34 Эе 35 Вг 36 Кг
37 ЯЬ 38 Эг 39 У 40 гг 41 № 42 Мо 43 Тс 44 1*11 45 46 Рс1 47 Ад 4« са 49 |П 50 Эп 51 эь 52 Те 53 54 Хе
56 Сэ 56 Ва 71 1-11 72 Ж 73 Та 74 \Л/ 75 Ре 76 Об 77 1г 78 Р1 79 Аи 80 нд 81 Т1 82 РЬ 83 В1 84 Ро 85 А1 86
87 Рг 88 103 1.Г 104 105 оь
-Ф
Рис. 1. Измерительные возможности РФА-СИ
-Ф-
Рис. 2. Объект исследования — озеро Кучерлинское
-Ф-
Рис. Схема станции «Элементный анализ»
затем выпиливались тонкие пластины с размерами 170 х 5 х 2 мм, которые анализировали в пучке СИ.
Схема эксперимента, оборудование. Эксперименты выполнялись в Институте ядерной физики им. Г. И. Будкера СО РАН в Сибирском центре синхротронного и терагерцового излучения. Использовалась станция «Элементный анализ», подробно описанная в работе Дарьин, Ракшун (2013). Упрощенно оборудование состоит из системы коллимации и фокусировки пучка, детектора и сканирующего устройства, которое протягивает керн с заданным шагом. Источником СИ служит накопитель ВЭПП-3. Как правило, определяется не менее 20 элементов, при этом в большинстве случаев предел обнаружения составляет от 1 до 5 г/1000 кг. Схема станции представлена на рис. 3. Для управления станцией используется специализированное программное обеспечение, разработанное в ИЯФ СО РАН. Программа обеспечивает подачу образца с заданным шагом сканирования и архивирование РФА спектров для последующей количественной обработки (см. более подробное методическое описание в работе Дарьин и др. 2013).
Были исследованы верхние 120 мм керна Ку-черла-2009 и 170 мм керна Кучерла-2018-2. Параметры эксперимента представлены в таблице 2. Определялись элементы: K, Ca, Ti, V, Cr, Mn, Fe, Ni, Zn, Ga, As, Br, Rb, Sr, Y, Zr, Nb, U.
Обработка результатов. Аналитические сигналы, полученные в результате измерений, представляют собой площади пиков аналитов, выделенные из спектрограммы при помощи программы обработки сложных РФА спектров (AXIL), предназначенной для энергодисперсионного спектрометрического анализа.
Программа обеспечивает пользователя информацией о площадях «чистых» пиков интен-сивностей (с учетом фона и частичного или полного наложения пиков) различных исследуемых элементов. Для учета нестабильности первичного возбуждающего излучения применяется аналитический сигнал, получаемый нормировкой площади пика аналита на площадь пика, рассеянного на образце излучения. Расчет содержаний каждого элемента в пробе производится м етодом внешнего стандарта с использованием государственного стандартного образца состава БИЛ-1 (Байкальский ил).
После совмещения профиля сканирования с фотографией исходного (влажного) керна было
Таблица 1
Керны озера Кучерлинское
Год Глубина места Длина керна, Маркировка
отбора отбора, м мм керна
2009 45 130 Кучерла-2009
2018 30 1120 Кучерла-2018-1
2018 45 865 Кучерла-2018-2
Таблица 2
Параметры сканирования исследуемых образцов кернов
Керн Энергия возбуждения, кэВ Размеры пятна возбуждающего излучения (мкм), высота — ширина Шаг сканирования (мкм) по высоте керна
Кучерла-2009 19,5 50-100 100
Кучерла-2009 23 500-2000 500
Кучерла-2018-2 19,5 50-100 100
Кучерла-2018-2 23 500-2000 500
-Ф-
-Ф-
проведено (рис. 4). Учитывая, что регистрация излучения определяемых микроэлементов происходит не только с поверхности, но и из глубины образца, что приводит к усреднению сигнала и «размыванию» границ слоев, совпадение экстремумов в распределении геохимических индикаторов с высокой точностью совпадает с визуально выделяемыми слоями. Таким образом, проведенный подсчет слоев дает возможность построения временной шкалы для современных осадков озера Кучерлинское на исследуемом 50-летним интервале 2008—1958 гг.
Результаты и обсуждение. Доказательством того, что выделенные слои являются действительно годовыми, служит совпадение построенной временной модели для керна Кучерла 2009 с данными по распределению изотопов 137С8 на этом интервале керна (рис. 4). Максимум активности в распределении 137С8 маркирует временной интервал 1961—1963 гг., что связано с глобальными выпадениями после атмосферных ис-
пытаний термоядерного оружия в конце 1961 г. Также небольшое повышение активности 137Cs в интервале 60—70 мм может являться следом Чернобыльской катастрофы 1986 г. Таким образом, полученная информация дает возможность определить координаты границ годовых слоев, начиная от верха керна, т. е. построить временную модель с годовым разрешением и рассчитать толщину каждого годового слоя.
Предполагая, что основным веществом донных отложений в озере Кучерлинское является терригенный материал, поставляемый паводковыми водами, было проведено сравнение мощностей годовых слоев с инструментальными данными ближайшей (рис. 2) метеостанции Кара-Тюрек, полученными из базы данных KNMI Climate Explorer (https://climexp.knmi.nl). Сравнение проводилось для временного интервала 1959—2006 гг. (рис. 5). Наблюдается хорошее визуальное совпадение общего характера двух графиков, что подтверждается значимой корреля-
g § S 14
$ к
о g
о о к
н о
1,2 1 0,8 0,6 0,4 0,2 0
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 105 110 115 120 125
Глубина керна, мм
Й 160
о
120 80
о
я 40
К
% 0
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 105 110 115 120 125
Глубина керна, мм
Рис. 4. Сопоставление экстремумов в нормированных распределениях КЬ/Бг отношения и содержания стронция
с границами визуально выделяемых слоев
-Ф-
-Ф-
§ 1 £ 0,9 | 0,8
Толщина слоя Атмосферные выпадения
1955 1960 1965 1970 1975
1980 1985 Годы, н. э.
1990 1995 2000 2005 2010
Рис. 5. Сравнение мощностей годовых слоев с инструментальными данными ближайшей метеостанции Кара-Тюрек
за период 1959—2006 гг.
цией между ними (коэффициент корреляции для 47 точек составлял величину +0,44).
Учитывая возможную погрешность в абсолютной датировке годовых слоев, был проведен расчет коэффициента корреляции для сдвига временной шкалы донных осадков на величину 2 года (рис. 6). Резкое падение коэффициента корреляции при сдвиге временной шкалы говорит о правильной датировке годовых слоев. Попытка найти связь мощности годовых слоев с температурой региона не принесла успеха — значимая корреляция не выявлена.
Временной интервал, перекрываемый керном Кучерла-2018-2, по полученным данным составляет 488 лет (2018—1530 гг.) при средней скорости осадконакопления 2,5 мм/год. Следующий этап работ предполагает проведение сравнения данных инструментальных метеонаблюдений не только с мощностью годовых слоев, но и с составом осадков, накопленных в течение годового цикла. Для этого необходимо провести математическую обработку данных микро-РФА сканирования с учетом границ выделяемых годовых слоев и пересчетом концентраций элементов в точках измерений в величины, количественно характеризующие процесс накопления осадочного материала в течение годового цикла. Подготовленные временные ряды литолого-геохимических дан-
-10 12 Сдвиг временной оси, года
Рис. 6. Расчет коэффициента корреляции для сдвига временной шкалы донных осадков на величину 2 года
ных, как правило, имеют более статистически значимые связи с инструментальными метеонаблюдениями.
Выводы. На препаратах кернов донных осадков приледникового озера Кучерлинское (Алтай) выполнены серии измерений (сканирующий микро-РФА) на энергиях возбуждения 19,5 и 23 кэВ с различным шагом сканирования с одновременным определением более 20 химических элементов. Проведена обработка спектров и построены профили распределения породообразующих и микроэлементов. Совмещение фотографических изображений поверхности образцов с профилями сканирования позволило получить оценку скоростей современного осадконакопле-ния в исследуемых озерах путем подсчета годовых слоев, выделяемых визуально и по геохимическим индикаторам. На основе полученных данных построены временные модели глубина керна — возраст осадка для последних десятилетий. Подготовлен эмпирический материал, позволяющий строить модели годового цикла современного осадкообразования в исследованном ледниковом озере и вести поиск аналитических зависимостей состава донных осадков от погодно-климатичес-ких факторов с годовым (сезонным) временным разрешением. Представляется перспективным проведение дальнейших исследований полученных кернов для построения реконструкций изменения гидрологических режимов озера за последние 450—500 лет.
Высокая детальность и количественная оценка реконструируемых параметров климата потенциально значительно расширяет наши представления о климатических изменениях последних тысячелетий и дает возможность обнаружения природных циклов и оценки их периодичности.
Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (проекты № 17-35-50115 и № 18-55-53016).
-Ф-
-Ф-
Библиографический список
1. De Geer, G. Zeitschrift fur allegemeine Geologie, 1912 v. 3, p. 457—471.
2. Ojala, A. E. K., Francus, P., Zolitschka, B., Besonen, M., Lamoureux, S. F., 2012. Quat. Sci. Rev. 43, 45—60.
3. Золотарев К. В., П. А. Пиминов, А. Д. Николенко, В. А. Трунова, Н. В. Полосьмак, К. Э. Купер, 2015 // Наука из первых рук. — 2015. — Т. 2(62). С. 10—18.
4. Baryshev V. B., Gavrilov N. G., Daryin A. V., Zolotarev K. V. et а1. Scanning x-ray fluorescent microanalysis of rock samples. Rev. Scientific Instruments. 1989. V. 60. № 7. Pt. II. P. 2456—2457.
5. Daryin A. V., Baryshev V. B., Zolotarev K. V. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment. 1991. V. 308. (1, 2). Р. 318—320.
6. Грачев М. А. и др. Геология и геофизика. 1997. Т. 38. (5). С. 957—980.
7. Phedorin M. A., Zolotarev K. V., Bobrov V. A. Nuclear Instruments and Methods in Physics Research. Section A: Accelerators, Spectrometers, Detectors and Associated Equipment. 1998. Т. 405. (2—3). С. 560—568.
8. Zolotarev K. V., Goldberg E. L., Kondratyev V. I., et al. Nuclear Instruments and Methods in Physics Research. Section A: Accelerators, Spectrometers, Detectors and Associated Equipment. 2001. Т. 470. (1—2). С. 376—379.
9. Goldberg E. L., Grachev M. A., Phedorin M. A. et al. Nuclear Instruments and Methods in Physics Research. Section A: Accelerators, Spectrometers, Detectors and Associated Equipment. 2001. Т. 470. (1—2). С. 388—395.
10. Daryin A. V., Kalugin I. A., Maksimova N. V. et al. // Nucl. Instruments and Meth. Phys. Res. 2005. V. A543. P. 255—258.
11. Phedorin, M. A., Bobrov, V. A., Chebykin, E. P., Goldberg, E. L., Melgunov, M. S., Filippova, S. V., Zolotarev, K. V. // Geostandards Newsletter: The Journal of Geostandards and Geoanalysis 2000a, 24, 205—216.
12. Phedorin, M. A., Bobrov, V. A., Goldberg, E. L., Navez, J., Zolotaryov, K. V., Grachev, M. A. // Nuclear Instruments and Methods in Physics Research A: Accelerators, Spectrometers, Detectors and Associated Equipment. — 2000. — Vol. 448. — P. 394—399.
13. Трунова В. А. Дисс. докт. хим. наук. — 2017. ИГУ: Иркутск.
14. Дарьин А. В., Ракшун Я. В. // Научный вестник НГТУ, 2013. № 2 (51). С. 119.
15. Дарьин А. В., Калугин И. А., Ракшун Я. В. // Изв. РАН. Сер. физ. 2013. Т. 77. № 2. С. 210.
16. Сороколетов Д. С., Ракшун Я. В., Дарьин Ф. А. Автометрия. 2015. Т. 51. № 3. С. 94—103.
17. Darin F. A., Kalugin I. A., Darin A. V., Rakshun Ya. V. Acta Geologica Sinica (English Edition). 2014. Т. 88. № S1. С. 5—6.
THE RATES OF SEDIMENTATION AND THE DETERMINATION OF THE BOUNDARIES OF ANNUAL LAYERS OF THE ANNUALLY-STRATIFIED DEPOSITS OF PROGLACIAL LAKE KUCHERLA (THE ALTAY) ACCORDING TO THE DATA FROM SCANNING MICRO-X-RAY FLUORESCENCE ANALYSIS BASED ON SYNCHROTRON RADIATION
F. A. Dar'in, Junior Researcher, Institute of Geography of the RAS, Moscow; Ph. D. Student, G. I. Budker Institute of Nuclear Physics, SB RAS, [email protected], Novosibirsk,
0. N. Solomina, Director, Institute of Geography of the RAS, Moscow,
A. M. Grachev, Senior Researcher, Institute of Geography of the RAS, Moscow,
A. V. Darin, Senior Researcher, V. S. Sobolev Institute of Geology and Mineralogy, Novosibirsk,
1. V. Rakshun, Scientific Secretary, G. I. Budker Institute of Nuclear Physics, SB RAS, Novosibirsk, D. S. Sorokoletov, Junior Researcher, G. I. Budker Institute of Nuclear Physics, SB RAS, Novosibirsk
References
1. De Geer, G. Zeitschrift fur allegemeine Geologie, 1912 v. 3, p. 457—471.
2. Ojala, A. E. K., Francus, P., Zolitschka, B., Besonen, M., Lamoureux, S. F., 2012. Quat. Sci. Rev. 43, 45—60.
3. Zolotarev K. V., P. A. Piminov, A. D. Nikolenko, V. A. Trunova, N. V. Polos'mak, K. E. Kuper Nauka iz pervikh ruk, 2015, P. 10—18 [in Russian]
4. Baryshev V. B., Gavrilov N. G., Daryin A. V., Zolotarev K. V. et al. Scanning x-ray fluorescent microanalysis of rock samples. Rev. Scientific Instruments. 1989. Vol. 60. No 7. Part II. P. 2456—2457.
5. Daryin A. V., Baryshev V. B., Zolotarev K. V. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment. 1991. Vol. 308. (1, 2). P. 318—320.
6. Grachev M. A. et al. Geologiya i geophizika. 1997. Vol. 38. (5). p. 957—980 [in Russian]
7. Phedorin M. A., Zolotarev K. V., Bobrov V. A. Nuclear Instruments and Methods in Physics Research. Section A: Accelerators, Spectrometers, Detectors and Associated Equipment. 1998. Vol. 405. (2—3). P. 560—568.
8. Zolotarev K. V., Goldberg E. L., Kondratyev V. I., et al. Nuclear Instruments and Methods in Physics Research. Section A: Accelerators, Spectrometers, Detectors and Associated Equipment. 2001. Vol. 470. (1—2). P. 376—379.
9. Goldberg E. L., Grachev M. A., Phedorin M. A. et al. Nuclear Instruments and Methods in Physics Research. Section A: Accelerators, Spectrometers, Detectors and Associated Equipment. 2001. Vol. 470. (1—2). P. 388—395.
10. Daryin A. V., Kalugin I. A., Maksimova N. V. et al. Nucl. Instruments and Meth. Phys. Res. 2005. V. A543. P. 255—258.
11. Phedorin, M. A., Bobrov, V. A., Chebykin, E. P., Goldberg, E. L., Melgunov, M. S., Filippova, S. V., Zolotarev, K. V. Geostandards Newsletter: The Journal of Geostandards and Geoanalysis 2000a, 24, 205—216.
12. Phedorin, M. A., Bobrov, V. A., Goldberg, E. L., Navez, J., Zolotaryov, K. V., Grachev, M. A. // Nuclear Instruments and Methods in Physics Research A: Accelerators, Spectrometers, Detectors and Associated Equipment. — 2000. — Vol. 448. — P. 394—399.
13. Trunova V. A. Thesis for Dr. Habil. 2017 [in Russian]
14. Darin A. V., Rakshun I. V. Nauchniy vestnik NGTU, 2013. No. 2 (51). P. 119 [in Russian]
15. Darin A. V., Kalugin I. A., Rakshun I. V. Isvestiya RAN. Seriya fizicheskaya, 2013. Vol. 77. No. 2. p. 210 [in Russian]
16. Sorokoletov D. S., Rakshun I. V., Darin A. V. Avtometriya, 2015. Vol. 51. No. 3. P. 94—103 [in Russian]
17. Darin F. A., Kalugin I. A., Darin A. V., Rakshun Ya. V. Acta Geobgica Sinica (English Edition). 2014. Vol. 88. No. S1. P. 5—6.
