Опыт экспериментального изучения петромагнитных индикаторов археологических палеокострищ на примере лёссов Киргизии
Обсуждаются перспективы применения петромагнитного метода для характеристики степени воздействия антропогенных костров на вмещающий субстрат и реконструкции различных параметров кострищ для целей археологии. Впервые проведено экспериментальное моделирование кострищ на лёссовом субстрате с использованием четырех типов топлива и проанализировано изменение магнитных свойств отложений в результате трансформации железосодержащих минералов. Повышенные значения магнитных параметров по сравнению с необожженными лёссами характерны и для пеплов, и для термически измененного лёсса, что делает их надежными петромагнитными идентификаторами древних кострищ.
First experimental study of rock-magnetic indicators of archaeological paleo-hearths on the example of Kyrgyzstan loess.pdf Введение Свидетельства использования огня являются важным источником археологической информации об особенностях быта и адаптационных стратегиях древнего человека. Зачастую пепловые прослои, остающиеся на местах кострищ, обнаруживаются археологами по изменению цвета на фоне вмещающих отложений (изменение в сторону красноват^хх, черных или серых оттенков) и по нахождению остатков угля. Действительно, эти данные являются прямыми визуальными свидетельствами существования огня, хотя его происхождение может быть как естественным, вызванным природными пожарами, так и антропогенным. Однако сохранения угольных остатков может не произойти, а интенсивность изменения цвета может слабеть с течением времени или этого изменения может вовсе не быть [Morinaga et al., 1999]. Это приводит к тому, что в полевой практике археологических исследований для обнаружения древних кострищ все большее применение начинают находить магнитн^1е методы, такие как съемка аномалий магнитного поля и измерение магнитной восприимчивости in situ. Установлено, что древние кострища создают над собой положительные аномалии магнитного поля и характеризуются увеличением значений магнитной восприимчивости по сравнению с вмещающими породами [Gibson, 1986; Jrad et al., 2014]. Начиная с ранних работ Le Borgne [1955, 1960], широко известен факт, что нагрев изменяет магнитную структуру почв и увеличивает их магнитную восприимчивость. Это свойство вызвало интерес в междисциплинарных исследованиях и стало находить широкий отклик в геоархеологических работах [Tite, Mullins, 1971; Mullins, 1974]. В дальнейшем исследователи стали формулировать методологические подходы, позволяющие находить следы существования огня на археологических стоянках по магнитным и петромагнитным данным, даже если часть первоначальных признаков воздействия огня не сохранилась в результате изменения и переработки отложений [Barbetti, 1986]. Развитием этих работ стал поиск методов для распознания костров природного происхождения (природные пожары) и антропогенного генезиса (кострища, очаги, печи и другие виды контролируемых человеком костров) [Barbetti, 1986; Bellomo, 1993]. © Кулакова Е.П., Дедов И.Е., Мещерякова О. А., Курбанов Р.Н., 2021 D0I: 10.17223/25421379/18/9 Одним из направлений применения магнитных методов в решении археологических задач является геофизическое площадное картирование стоянок по аномалиям магнитного поля и значениям магнитной восприимчивости [Матасова и др., 2013, 2016]. Магнитная съемка стала важным методом в практике археологии, являясь недеструктивным методом с возможностью картирования больших площадей. Было установлено, что положительные аномалии магнитного поля над древними очагами связаны с повышенными концентрациями магнитных минералов [Bonhomme, Stanley, 1986; Gibson, 1986; Linford, Canti, 2001; Jrad et al., 2014]. Локальн^хе магнитн^хе аномалии и зоны с повышенными значениями магнитной восприимчивости in situ являются признаками возможных древних очагов и определяют наиболее перспективные участки для археологических раскопок [Смекалова и др., 2007]. Другим направлением является лабораторное пет-ромагнитное изучение образцов с целью характеристики отложений из кострищ и изменений, происходивших в ходе нагрева. Применение разнообразных петромагнитных параметров стало возможным благодаря развитию технической базы и увеличению чувствительности приборов. Данное направление тесно связано с экспериментальн^1м моделированием, проводимым для сравнения полученных результатов с природн^1ми объектами и дальнейшего решения обратных задач. Особый интерес представляют пепло-вые отложения, изучение которых проводится для определения типа топлива [Peters et al., 2001, 2002; Church et al., 2007] и, соответственно, адаптационн^хх стратегий древнего человека. В связи с тем, что прямые свидетельства существования костра (пепловые прослои или угольные остатки) могут быть эродированы или перемещены в результате деятельности человека или постседиментационных процессов, идентификация кострищ может быть основана на обнаружении термически измененных пород, служивших субстратом для древнего кострища. Многими авторами проводилось экспериментальное моделирование костров с дальнейшим изучением магнитных свойств термически измененных пород для установления степени воздействия на исследуемый субстрат (температура прогрева и минералогические изменения) [Morinaga et al., 1999; Linford, Canti, 2001; Jrad et al., 2014; Aldeias et al., 2016]. Стоит отметить, что эти параметры довольно вариативны и зависят от продолжительности горения, максимальной температуры костра, исходного минералогического состава пород, влажности и плотности субстрата, окислительно-восстановительных условий среды, повторного воспроизведения костров на том же месте. Таким образом, на данный момент разработана общая методика обнаружения древних кострищ по петромагнитным данным и обозначены основные изменения, которые могут происходить при воздействии высоких температур. Имеется опыт успешного применения экспериментального моделирования в решении задач при исследовании природных археологических объектов [Maki et al., 2006; Jrad et al., 2014]. Однако в преобладающем большинстве работы предшественников имеют лишь частное применение и могут быть использованы для археологических стоянок с аналогичными литологическими и палеокли-матическими условиями. Учитывая разнообразие типов отложений, широкий спектр климатических условий, характерных для того или иного археологического объекта, а также многообразие возможных источников топлива, метод экспериментального моделирования является индивидуальным для каждого археологического контекста и актуальн^хм. Экспериментальное моделирование является необходимым условием для определения особенностей изменения петромагнитных свойств отложений. Целью исследования является оценка перспектив применения петромагнитных методов для реконструкции характеристик древних кострищ, встречающихся в лёссово-почвенн^1х сериях различн^хх районов Центральной Азии. В своей работе мы представляем результаты экспериментального моделирования костров на лёссовом грунте с использованием различного типа топлива и анализ их магнитных свойств. В качестве исходного субстрата может выступать материал как из горизонтов лёсса, так и погребенных почв, выбор которого определяется конкретной геоархеологической задачей. В данном исследовании эксперимент проводился с позднеплейстоценовыми лёссами (голодностепский комплекс Ферганской долины) из разреза в районе стоянки Обишир (Баткенская область, Киргизия). Материал был отобран из горизонта лёссов между современной почвой и педокомплексом, отвечающим третьей изотопно-кислородной стадии, ввиду частого обнаружения в подобных осадках стоянок позднего палеолита и эпипалеолита [Додонов, 2002; Шнайдер, 2015]. В связи с необходимостью определения характеристик кострищ при изучении палеолитических памятников Центральной Азии, нами предпринята попытка разработки методики проведения подобных работ на основе анализа изменения петромагнитных свойств отложений в ходе контролируемого эксперимента. Для этого необходимо выявить причину увеличения магнитной восприимчивости, регистрируемого полевым каппаметром, и определить возможности использования этого параметра как средства быстрой идентификации древних кострищ в полевых условиях. Отдельной задачей является поиск путей для определения степени влияния типа топлива на магнитные свойства пепла и термически измененного субстрата. Методика исследований и используемый материал Экспериментальные костры. Нами б^хли проведен контролируемый эксперимент, в ходе которого изучены петромагнитные характеристики четырех костров с использованием различного вида топлива (рис. 1). В качестве топлива были выбраны материалы, доступные и используемые в географическом регионе исследования: кость (сухая), терескен (кустарник, распространенный на территории Центральной Азии), древесина (крупные ветки и поленья), кизяк (высушенн^1й навоз). Для первоначального розжига костра с использованием кости как основного вида топлива также использовались ветки терескена и деревьев. В качестве подложки для всех костров был использован лёсс. В связи с тем, что изучался материал из средней части разреза, а также для создания единых условий влажности и пористости, нами были отобраны лёссы с одной глубины, из которых мы подготовили искусственные насыпи. Каждое экспериментальное кострище представляло собой круглую насыпь лёсса диаметром 30-35 см и высотой 10-12 см с углублением на 3-5 см в центре. В ходе эксперимента наблюдались различия в долготе и интенсивности горения, а соответственно, и в достигаемых температурах прогрева нижележащих отложений в связи с различной калорийностью используемого топлива. Для дальнейшего петромаг-нитного изучения было отобрано по три образца из каждого костра: образец лёссовой подложки до начала эксперимента, образец обожженного лёсса из приповерхностного слоя (0-2 см) непосредственно под местом горения и образец пепла, образовавшегося в результате горения. Рис. 1. Экспериментальное моделирование костров a-c - фотографии различн^1х стадий эксперимента на примере костра с использованием терескена как топлива: a - лёссовый субстрат (подложка) до костра, b - горение топлива, c - кострище, образование пеплового слоя; d - общий вид на костры в течение эксперимента Fig. 1. Experimental fire modeling a-c - photographs of various stages of the experiment with using winterfat as a type of fuel: a - loess substrate before the fire, b - fuel combustion, c - hearth, formation of an ash layer; d - general view of the fires during the experiment Петромагнитные измерения. Петромагнитн^хе исследования являются быстрым, относительно дешевым и простым в измерении инструментом для идентификации изменений, происходящих в субстрате при его прогреве. В ходе исследований может быть получена информация о магнитной минералогии пород, концентрации и размере магнитных минералов [Evans, Heller, 2003]. Нами были проведены измерения следующих магнитных параметров. 1. Удельная магнитная восприимчивость (х) - величина, характеризующая связь намагниченности единицы массы вещества с магнитным полем в этом веществе. Ферромагнитные и парамагнитные минералы имеют положительные значения магнитной восприимчивости и вносят основной вклад в суммарную восприимчивость, в то время как диамагнитные минералы, обладая отрицательными значениями магнитной восприимчивости, занижают значения суммарной восприимчивости. 2. Температурная зависимость магнитной восприимчивости показывает изменение значений магнитной восприимчивости в ходе нагрева до заданной температуры и последующего охлаждения от нее (нагрев до 700 °С). Анализ термомагнитн^хх кривых позволяет определить основные минералы-носители намагниченности и помогает отследить минеральные преобразования в ходе нагрева и охлаждения. Измерения магнитной восприимчивости, ее температурной и частотной зависимостей осуществлялись на каппабридже MFK1-FA производства AGICO (Чехия). 3. Частотная зависимость магнитной восприимчивости (xfd) является величиной, отражающей различие значений магнитной восприимчивости при измерении ее на низкой и высокой частотах. Данный параметр чувствителен к присутствию очень мелких зерен магнетита (менее 0,03 мкм), находящихся в суперпарамагнитном состоянии (SP) [Bearing et al., 1996]. Частотная зависимость xFD в абсолютных величинах и в процентном выражении определяется следующим образом: xfd = xlf - Xhf, Xfd (%) = (xlf - xhf)/ xlf X 100, где xlf и xhf - значения магнитной восприимчивости на низкой и высокой частоте соответственно. При определении частотной зависимости на кап-пабридже MFK1-FA измерения магнитной восприимчивости осуществлялись на рабочих частотах LF = 976 Гц и HF = 15616 Гц. В связи с тем, что в истории петромагнитных исследований долгое время использовался иной инструмент, Bartington MS-2 Susceptibility Meter, с рабочими частотами LF = 465 Гц и HF = 4650 Гц, имеется необходимость в ведении поправки к получаемым нами данным, учитывающей различие в используемых частотах. Поэтому для дальнейшего сравнения с литературными данными мы будем использовать вместо хрв такой параметр, как xfb, предложенный в работе [Hrouda, 2011] и определяемый как In 10 (1) -Xfd , Xfb = In/HF -In/LF где fhf и flf - значения высокой и низкой частот, на которых производят измерения. Параметр xfb соотносит логарифмическую разницу между используемыми частотами к 10, т.е. отношению частот в Bar-tington. Очевидно, что для измерений, проведенных на Bartington, xfb = Xfd, а для иных инструментов параметр xfb будет отличаться согласно отношению используемых частот. 4. Идеальная (безгистерезисная) остаточная намагниченность (ARM) - остаточная намагниченность, созданная в убывающем от некоторого значения до нуля переменном поле (AF) в присутствии постоянного магнитного поля (DC). В наших экспериментах ARM создавалась при параметрах DC = 0,05 мТл и AF = ±130 мТл и измерялась на криогенном SQUID магнитометре производства фирмы 2G Enterprises (США). ARM чувствительна как к концентрации магнитн^1х минералов, так и к их размеру. 5. Остаточная намагниченность насыщения (SIRM) - величина остаточной намагниченности, образующейся в образце после мгновенного воздействия на него постоянного магнитного поля, равного или превосходящего поле магнитного насыщения образца. Является крайним случаем изотермической остаточной намагниченности (IRM), которая образуется подобным образом, но в магнитных полях меньших полей насыщения. SIRM создавалась на образцах, предварительно размагниченных переменным полем (AF = 130 мТл), при помощи установки импульсного намагничивания ASC IM-100 посредством воздействия коротким импульсом магнитного поля DC = 1 Тл. Созданная остаточная намагниченность насыщения измерялась на криогенном SQUID магнитометре. Величина SIRM может быть использована для оценки общей концентрации ферромагнитных минералов в породе. 6. S-ratio - параметр, определяющий относительный вклад магнитных минералов с разной коэрцитив-ностью в остаточную намагниченность насыщения. Данный параметр является эффективной мерой при оценке соотношения в образце количества «магнито-мягких» (например, магнетит и (или) маггемит) и «магнитожестких» (например, гематит и (или) гетит) минералов. S-ratio определяется следующим образом: прямое поле интенсивностью 1 Тл прикладывается к образцу, после чего измеряется его остаточная намагниченность (зачастую это SIRM), далее прикладывается обратное поле (в практике петромагнитн^хх исследований используют обратные поля, соответствующие -100 мТл и -300 мТл), и также измеряется остаточная намагниченность образца (IRM). Параметр рассчит^1вается следующим образом: S-ratiOlooмТл = - (IRM_looмТл/SIRM+loooмТл), S-ratioзooмТл = - (IRM_зooмТл/SIRM+loooмТл). Расчет S-ratio выполнен для обоих обратных полей. Измерения намагниченности произведены на вибромагнитометре (VSM) PMC MicroMag 3900. 7. Показатель магнитной жесткости HIRM отражает абсолютный вклад магнитожестких минералов в общую остаточную намагниченность образца. При его расчете используются те же значения остаточной намагниченности в прямом и обратных полях, что и при определении S-ratio: HIRM.looмТл = (SIRM+loooмТл + IRM_looмТл)/2, HIRM.зooмТл = (SIRM+loooмТл + IRM_зooмТл)/2. 8. Гистерезисные параметры включают в себя намагниченность насыщения (Ms), остаточную намагниченность насыщения (Mrs), коэрцитивную силу (Bc), получаемые непосредственно в процессе снятия петли гистерезиса, а также остаточную коэрцитивную силу (Bcr). Гистерезисные параметры помогают в определении магнитной минералогии и дают информацию о доменном состоянии (размере) зерен. Измерения производились на вибромагнитометре (VSM) PMC MicroMag 3900 и включали снятие петель гистерезиса, кривой нормального намагничивания и остаточной намагниченности в обратных полях. По данным гисте-резисных параметров была построена диаграмма Дэя-Данлопа, показывающая зависимость Mrs/Ms от Bcr/Bc [Day et al., 1977; Dunlop, 2002a, 2002b]. 9. Анализ спектров коэрцитивности образцов по методу «cumulative log-Gaussian analysis» (CLGA) кривой нормального намагничивания [Kruiver et al., 2001; Heslop et al., 2002] производился в веб-приложении MAX UnMIX [Maxbauer et al., 2016]. Все петромагнитные исследования проводились в лаборатории Главного геомагнитного поля и петро-магнетизма Института физики Земли РАН имени О.Ю. Шмидта. Для каждого исследуемого образца в связи с его естественной неоднородностью были сделаны дубли в количестве 3-4 штук. Основные значения петромагнитных параметров, таким образом, являются осреднением показателей дублей. Результаты и интерпретация 1. Экспериментальные костры. В нашем эксперименте наблюдались различные излучательные световые интенсивности костров, процессы горения (тление или открытый огонь), ощущаемый жар от костров и скорости накопления пеплового прослоя в зависимости от используемого типа топлива. Наибольшие излучение света и ощущаемый жар исходили от костров из древесины и терескена, отличающихся также стабильностью и легкостью поддержания. Сильнейший жар был и у костра с использованием сухих костей, однако его горение было нестабильным, а выгорание костей довольно быстрым. Костер с использованием кизяка характеризовался непродолжительным и малым в высоту пламенем, но очень длительной по продолжительности стадией тления с высоким тепловыделением. Костры из терескена и кизяка для долговременного поддержания требовали большой объем материала, что привело к быстрому формированию мощного пеплово-го прослоя. В связи с различной скоростью выгорания и разного первоначального объема материалов длительность горения костров была различной (табл. 1). В обобщающей работе [Aldeias, 2017] приводятся следующие величины максимальных температур, достигаемых в экспериментальных кострах, использующих различные типы топлива. В кострах на основе костей с небольшой примесью древесины в среднем достигаются температуры 605-825 °С вплоть до 900 °С даже при продолжительности горения порядка 1-2 ч. Костры из древесины при непродолжительном горении в первые часы дают средние температуры в диапазоне 465-760 °С, однако при подержании таких костров более суток могут быть достигнуты пиковые температуры 900-1 000 °С. Данн^1е по использованию высушенного навоза крупного рогатого скота как топлива свидетельствуют о пиков^хх температурах 630-800 °С, информация о продолжительности горения не приводится. Данные о кострах с использованием кустарникового типа топлива в литературе не приводятся, но ввиду особенностей строения стеблей кустарников по сравнению с деревьями мы предполагаем, что средние температуры не должны превосходить таковых для древесины и, вероятно, являются более низкими. 2. Магнитная минералогия. Значения основн^хх петромагнитных параметров приведены для образцов исходного лёссового субстрата и обожженных отложений в табл. 2, для образцов пепла - в табл. 3. Для всех образцов лёссов, подвергнутых температурному воздействию, и образцов пепла характерны более высокие значения удельной магнитной восприимчивости (х), частотной зависимости магнитной восприимчивости (xfd, xfb) и остаточн^1х намагниченностей SIRM и ARM по сравнению с этими показателями для лёссов. Степень увеличения магнитных свойств вариативна от образца к образцу. Исходный лёсс имеет средние значения SIRM - 5,0 мА м2/кг, ARM - 0,047 мА м2/кг. Величин^! SIRM и ARM возрастают в 1,1-4,0 и 1,13,6 раз в образцах обожженного лёсса и в 1,2-2,9 и 1,63,4 раз выше для образцов пепла (рис. 2). Параметр S-ratio оценивает пропорцию менее коэрцитивных минералов к более коэрцитивным, где порог разграничения минералов по жесткости определяется выбранным значением поля (100 или 300 мТл). Результаты показывают (табл. 2), что подверженные температурному воздействию образцы лёссов обладают более высокими значениями S-ratio100мТл, в то время как параметр S-ratioзooмТл имеет незначительное увеличение. Кроме того, намагниченность, связанная с магнитожесткими минералами (HIRM), в обратном поле 100 мТл для практически всех обожженных пород лёсса значимо уменьшается, а в обратном поле 300 мТл испытывает лишь небольшое понижение значений (за исключением образца Кс-2). Это свидетельствует, что наблюдаемое увеличение магнитной восприимчивости и остаточных намагниченностей обожженных пород связано одновременно и с повышением концентрации низко-коэрцитивн^1х минералов (с коэрцитивностью до 100 мТл) и с разрушением более магнитожестких (с преимущественной коэрцитивностью 100-300 мТл, но также и более) по сравнению с исходным субстратом. Для образцов из костров с терескеном и кизяком, где рост намагниченностей незначителен, происходит более интенсивное разрушение относительно магнито-жестких минералов, нежели новообразование магни-томягких. Таблица 1 Table 1 Параметр Тип топлива Кость (+кустарник) Древесина Кизяк Терескен Длительность горения Высота пламени Скорость в^1горания топлива Остатки 2 ч средняя (20-40 см) 4 кг/ч Карбонатизированные остатки костей, пепла мало 6 ч высокая (40-100 см) 2 кг/ч Фрагмент^! углей, чистый бел^1й пепел 4 ч низкая (0-20 см) 3 кг/ч Много серого неплотного пепла 1,5 ч высокая (40-100 см) 10 кг/ч Много плотного пепла с небольшими угольными фрагментами Таблица 2 Значения петромагнитных параметров для образцов исходного лёссового субстрата (1) и обожженных отложений непосредственно под местом разведения экспериментальных костров (2) Характеристика экспериментальных костров по данным наблюдения Characteristics of experimental fires according to observation data Values of rock-magnetic parameters for samples of the initial loess substrate (1) and burnt sediments directly under the place of the experimental fires (2) Значения петромагнитных параметров для образцов пепла Values of rock-magnetic parameters for ash samples Table 2 Образец Описание X (^ 10-6 м3/кг) %FB, % (X 10-8 м3/кг) SIRM (мА м2/кг) ARM (мА м2/кг) SratiOlooмТл SratioзooмТл HIRMlooмТл (мА м /кг) HIRMзooмТл (мА м2/кг) Исходный лёссовый Кс-1 субстрат (до разведения костра) Обожженный лёсс 0,41 5,2 2,6 5,0 0,045 0,51 0,92 1,09 0,17 Кс-2 (приповерхностн^1й слой под костром из кости) 1,55 7,4 13,9 11,5 0,142 0,87 0,97 0,88 0,20 Исходный лёссовый Др-1 субстрат (до разведения костра) Обожженный лёсс 0,41 5,1 2,5 4,8 0,044 0,52 0,93 1,20 0,19 Др-2 (приповерхностн^1й слой под костром из древесины) 1,72 5,8 12,0 19,4 0,162 0,89 0,99 1,19 0,15 Исходный лёссовый Тр-1 субстрат (до разведения костра) Обожженный лёсс 0,42 5,1 2,6 5,0 0,046 0,52 0,92 1,18 0,19 Тр-2 (приповерхностн^1й слой под костром из терескена) 0,54 5,8 3,8 5,4 0,057 0,71 0,93 0,70 0,17 Исходный лёссовый Кз-1 субстрат (до разведения костра) Обожженный лёсс 0,46 5,5 3,1 5,2 0,052 0,54 0,93 1,19 0,19 Кз-2 (приповерхностн^1й слой под костром из кизяка) 0,54 5,0 3,3 6,2 0,055 0,74 0,94 0,75 0,17 Таблица 3 Table 3 Образец X (X 10-6 м3/кг) %FB, % %FD (X 10-8 м3/кг) SIRM2 (мА м /кг) ARM (мА м2/кг) S-ratiOlooмТл S-ratioзooмТл HIRMlooмТл (мА м2/кг) HIRMзooмТл (мА м2/кг) Кс_пепел 1,20 7,3 10,5 9,0 0,088 0,86 0,97 0,61 0,14 Др_пепел 1,07 7,1 9,1 9,3 0,096 0,83 0,96 0,65 0,13 Тр_пепел 1,71 7,7 15,6 14,3 0,155 0,85 0,96 0,74 0,18 Кз пепел 0,69 6,5 5,4 6,1 0,070 0,84 0,97 0,81 0,16 Рис. 2. Величины петромагнитных параметров для образцов исходного лёсса, обожженного лёсса и пепла Магнитная восприимчивость (х) в^хражена в м3/кг, SIRM и ARM - мАм2/кг Fig. 2. Values of rock-magnetic parameters for samples of original loess, heated loess and ash Magnetic susceptibility (x) is expressed in m3/kg, SIRM and ARM are expressed in mAm2/kg Рис. 3. Диаграмма значений удельной магнитной восприимчивости (х) и частотной зависимости магнитной восприимчивости (xFD) для всех дублей образцов Fig. 3. Plot of frequency dependence of the magnetic susceptibility (xFD) against magnetic susceptibility (х) for all specimens 2.1. Магнитная восприимчивость. Результат^! экспериментального моделирования обнаруживают ожидаемое увеличение значений магнитной восприимчивости (МВ) как для образцов обожженного лёсса, так и для пепла. На рис. 3 приведено сравнение значения удельной МВ (х) и частотной зависимости МВ (хко). Изначальный лёссовый субстрат показывает среднее значение удельной МВ 0,43 x 10-6 м3/кг и диапазон значений частотной зависимости МВ 2,0-3,5 x -8 3 -8 3 тт 10 м /кг со средним значением 2,7 x 10 м /кг. Показатели выраженной в процентах частотной зависимости МВ (xfb) составляют 4,7-6,0 % со средним значением 5,2 %, что свидетельствует о значительном присутствии мельчайших магнитных частиц, находящихся в суперпарамагнитном состоянии, в первичном лёссе [Bearing et al., 1996]. Все пепловые образца! обнаруживают повышенные значения как по удельной, так и по частотной зависимости МВ по сравнению с лёссовым субстратом и имеют значения, варьирующие в диапазоне 0,61-1,87 x 10 6 м3/кг для удельной МВ, 5,0-18,6 X 10-8 м3/кг и 6,4-8,3 % - для абсолютной и процентной частотной зависимости МВ соответственно. Удельная МВ пеплов четко коррелирует с частотной зависимостью МВ (рис. 3). Это свидетельствует о том, что повышение МВ вызвано увеличением магнитной концентрации суперпарамагнитных частиц. Образцы обожженного лёсса тоже характеризуются повышением значений удельной МВ и абсолютной частотной зависимости МВ (за исключением образца обожженного лёсса из костра на кизяке), однако степень их корреляции ниже, чем у пеплов, что говорит о более разнообразном изменении магнитной минералогии. Значения удельной МВ варьируют в широком диапазоне 0,52-1,83 X 10-6 м3/кг и зависят, скорее всего, от степени прогрева лёсса. Повышение МВ связано с образованием в ходе нагрева новых магнитных фаз (магнетит/маггемит) из немагнитных или слабомагнитных (см. подраздел 2.2). 2.2. Температурная зависимость магнитной восприимчивости. Температурная зависимость магнитной восприимчивости является чувствительным параметром как к размеру зерен, так и к изменениям в магнитной минералогии, происходящим в ходе цикла нагрева и охлаждения. Для всех образцов, полученных в ходе эксперимента, в первом цикле кривая охлаждения идет значительно выше кривой нагрева, свидетельствуя об образовании новых сильных магнитных фаз в ходе прокаливания (рис. 4). Второй цикл нагрева, в свою очередь, не вызывает значительного увеличения магнитной восприимчивости. Поведение температурных кривых первого цикла исходного лёссового субстрата схоже с таковым для лёссов Центральной Азии [Zan et al., 2012; Song et al., 2018], Китая [Liu et al., 2005; Zan et al., 2017] и лёссов Юго-Восточной Украина! [Jelenska et al., 2010]. Из-за существенной разниц^! в вертикальном масштабе между кривыми нагрева и охлаждения кривая нагрева также показана отдельно на врезке (рис. 4). Для исходного лёссового субстрата поведение кривой нагрева описывается следующим образом: небольшой рост до температуры ~310 °С, выраженное падение значений МВ в интервале 310460 °С, увеличение МВ и пик в районе 560 °С, резкий спад значений в области температур 560-590 °С и слабое постепенное уменьшение МВ от 590 до 700 °С. Первоначальн^1й небольшой рост МВ до температуры ~310 °С может быть обусловлен несколькими причинами: постепенное деблокирование мелких однодоменных частиц и их переход в суперпарамагнитное состояние при повышении температуры [Liu et al., 2005]; дегидратация гидроксидов железа, таких как лепидокрокит и гетит, с последующим образованием маггемита и (или) гематита соответственно [Evans, Heller, 2003]. Мы предполагаем, что в изучаемых нами породах может реализо-вываться комбинация описанн^хх выше механизмов. Падение магнитной восприимчивости в интервале 310-460 °С интерпретируется как переход тонкозернистого температурно-нестабильного маггемита в более устойчивый, но характеризующийся низкими значениями МВ гематит. Данный спад на температурной кривой и его связь с маггемитом надежно установлены в лёссово-почвенных сериях Китая [Sun et al., 1995; Deng et al., 2001; Liu et al., 2005]. Дальнейший рост магнитной восприимчивости вплоть до пика при ~560 °С, отчетливо проявляющегося на кривой нагрева, связан с новообразованием магнетита, который также уверенно идентифицируется по резкому спаду МВ в диапазоне температур 560-590 °С. Постепенное уменьшение МВ до 700 °С указывает на гематит. Кривая первого охлаждения исходного лёсса стремительно растет от 600 °С до температур 400-450 °С, свидетельствуя об новообразовании значительного количества магнетита в ходе эксперимента. От температуры 400 °С кривая охлаждения имеет волнообразный вид и зачастую второй небольшой пик в районе 200 °С, имеющий отношение, скорее всего, к структурным переходам доменного состояния, нежели к магнитн^1м фазам. Для крив^хх охлаждения (а также кривых второго нагрева и охлаждения) всех образцов характерен наклон в сторону более низких величин МВ на начальных температурах, который свидетельствует о присутствии большого количества ультратонких новообразованных частиц, находящихся в пограничном между SP и SD состоянии. Все минеральные трансформации первого нагрева происходят в восстановительной обстановке, обусловленной содержанием гумуса в лёссе, что также было подтверждено мессбауэрской спектроскопией для лёссов Юго-Восточной Украины [Jelenska et al., 2010]. Восстановительные условия способствуют т- 3+ переходам Fe содержащих минералов в магнетит, что в результате сильно увеличивает конечные значения МВ по сравнению с исходными. Стоит отметить, что восстановительная обстановка при первом нагреве наблюдается для всех типов образцов (исходный лёсс, обожженный лёсс, пепел). Пепел содержит органические остатки, а для обожженного субстрата либо имеется смешение с более «свежим» лёссом из-за резкого падения температуры с глубиной, либо происходит проникновение мельчайших органических остатков из вышележащего пепла. Температурные кривые для образцов обожженного лёсса визуально практически идентичны кривым исходного лёсса для костров из кизяка и терескена, свидетельствуя о прогреве субстрата до температуры не более 400 °С в экспериментальном костре (вероятно, даже не более 250 °С и до первых минералогических преобразований). Образцы обожженного лёсса из-под костров с использованием древесины и кости имеют исходные более высокие значения МВ и меньшую конечную разницу в величине МВ после первого нагрева. По всей видимости, для этих образцов были достигнуты температуры свыше 400 °С и началось образование магнетита, но все же прогрев был недостаточным для полного восстановления гематита в магнетит. Рис. 4. Кривые температурной зависимости магнитной восприимчивости для всех типов образцов из экспериментальных костров с использованием различного типа топлива Сплошные линии соответствуют кривым нагрева (красные) и охлаждения (синие) первого цикла, пунктирные линии - второго цикла для того же образца Fig. 4. High-temperature variations of magnetic susceptibility for each types of samples from experimental fire with using different fuel types The solid lines correspond to the heating (red) and cooling (blue) curves of the first cycle, the dashed lines - the second cycle for the same sample Образцы пепла в большинстве своем имеют изначально более высокие значения МВ, чем исходн^хй лёсс. Для всех образцов кривая первого нагрева до ~500 °С имеет слегка «выпуклый» с плавн^хм ростом и уменьшением вид (рис. 4, врезки), свидетельствующий, по всей видимости, об отсутствии или незначительном присутствии нестабильных при умеренных температурах минералов (как гидроксиды и силикаты железа). Мы полагаем, что присутствие нестабильных гидроксидов вызвано контаминацией пепла с лёссом, произошедшей при отборе, поэтому кривые первого нагрева образцов наследуют в разной степени черты кривых для лёсса. Мы считаем, что образец пепла из костра с использованием терескена (см. рис. 4, врезка) в меньшей степени подвержен контаминации и поэтому позволяет утверждать, что основным исходным магнитным минералом для образцов пепла является тонкозернистый магнетит и, возможно, гематит, чье восстановление приводит к образованию нового магнетита и увеличению МВ. Поведение кривых второго нагрева и охлаждения свидетельствует о более простой минералогии, включающей превалирующее количество ультратонкого магнетита и для некоторых образцов небольшое количество гематита. Кривые второго нагрева характеризуются постепенным ростом МВ до 440-460 °С с последующим резким спадом, однако для некоторых образцов наблюдается небольшое уменьшение МВ в интервале 250-400 °С, которое еще сильнее подчеркивает пик. Для этих же образцов отмечается более высокая, чем для магнетита, температура Кюри ~600 °С. Вероятно, в ходе второго нагрева таких образцов, водного пара, образовавшегося при дегидратации минералов при первом нагреве, недостаточно для поддержания восстановительной обстановки, и условия меняются на окислительные из-за воздействия воздуха. Образованный в ходе первого цикла магнетит частично окисляется и замещается маггемитом новой генерации. 2.3. Гистерезисные параметры. Петли гистерезиса (нормированн^1е на массу образцов и скорректированные за диа- и парамагнитную составляющие) для образцов обожженного лёсса и пепла визуально практически полностью закрываются к 150 мТл и окончательно к 300-350 мТл (рис. 5). Это свидетельствует о преимущественном присутствии низкоко-эрцитивн^1х минералов (магнетит, маггемит) и небольшом количестве среднекоэрцитивн^хх (мелкозернистый гематит). Преобладание низкокоэрцитивных минералов справедливо и для исходного лёсса, однако полное закрытие петель происходит при более высоких значениях поля (около 550 мТл и даже выше). Это говорит о присутствии высококоэрцитивного магнитного минерала (вероятно, гетита), по всей видимости, преобразующегося в низкокоэрцитивные фазы при температурном воздействии. Все образцы попадают в поле псевдооднодоменных частиц на диаграмме Дэя-Данлопа (рис. 5), располагаясь параллельно кривой смешения одно- и многодоменных частиц (SD-MD mixing curve). Смещение вправо, в область более высоких значений Bcr/Bc, обусловлено присутствием суперпарамагнитных частиц во всех типах образцов, а также может быть связано с наличием средне-высококоэрцитивного мелкозернистого гематита. Образцы исходного лёсса (до нагрева) показывают более сильное отклонение вправо, связанное с одновременным присутствием в образцах значительного количества высококоэрцитивного минерала, гетита. Образцы пепла более компактно сгруппированы, чем образцы обожженного лёсса, широкое распределение последних, по-видимому, связано с различиями в минералогическом составе. 2.4. Компонентный анализ IRM. Одним из методов разделения магнитных минеральных ассоциаций на отдельные фазы является анализ кривой приобретения изотермической остаточной намагниченности (IRM), целью которого является выделение компонент по их магнитной коэрцитивности. В данном методе осуществляется математическое моделирование кривой первой производной зависимости IRM от приложенного поля с использованием ряда отдельных логнормальных функций плотности вероятности, каждая из которых отвечает отдельной магнитной фазе (метод «cumulative log-Gaussian analysis» (CLGA), описанный в работах [Kruiver et al., 2001; Heslop et al., 2002]). Мы произвели анализ для представительных образцов всех отобранных типов: для исходного, неизмененного нагревом лёссового субстрата (рис. 6, a) на трех дублях, для образцов обожженного лёсса из каждого экспериментального костра (рис. 6, b), для образцов пепла каждого типа топлива (рис. 6, c). Результаты компонентного анализа представлены в табл. 4. Полученные данные по разделению магнитных фаз выражены через три параметра: 1) относительный вклад каждой модельной компоненты, выраженный в процентном содержании; 2) B1/2, поле насыщения, в котором магнитная фаза приобретает половину своей остаточной намагниченности насыщения (SIRM); 3) DP - параметр дисперсии, отображающий стандартное отклонение логнормального распределения каждой магнитной фазы. Для образцов неизмененного лёсса сгагисгиче-ская модель лучшим образом описывается четырьмя кoмпoненгами, в то время как для образцов обожженного лёсса и пепла она преимущесгвеннo содержит три компоненты (рис. 6 и табл. 4). Первая кoмпoненга вносит наибольший вклад в остаточную намагниченнoсгь насыщения (SIRM) для образцов всех типов, составляя в среднем 6890 %, и имеет значения медианного поля насыщения B1/2 в диапазоне 27-50 мТл. Данные значения B1/2 характерны для магнетита и (или) маггемита, которые, по всей вероятности, являются носителями первой компоненты. Среди всех образцов данная ком-пoненга имеет более низкие значения B1/2 (2735 мТл) для образцов, подвергшихся воздействию сильных нагревов (как Кс-2 и Др-2), и образцов пепла. Меньшие поля насыщения характерны для зерен меньшей размерности, поэтому, по всей видимости, первая компонента в данных образцах преимущественно связана с ультратонким магнетитом. В свою очередь, образцы Кз-2 и Тр-2 (~43 мТл) имеют значения B1/2 близкие к первичному лёссу (47-50 мТл). Первая компонента этих образцов в большей степени связана с педогенным маггемитом, регистрируемым на термокривых МВ, хотя образцы Кз-2 и Тр-2 могут содержать в себе и небольшое количество новообразованного ультратонкого магнетита. Рис. 5. Диаграмма Дэя-Данлопа (Mrs/Ms от Bcr/Bc), отображающая доменную структуру магнитных минералов, и петли гистерезиса, характерные для каждого вида образцов [Day et al., 1977; Dunlop, 2002a, 2002b] SD - однодоменные зерна, PSD - псевдооднодоменные зерна, MD - многодоменные зерна, SP - суперпарамагнигн^хе зерна. Кривые смешения SD + MD и SD + SP нанесен^! по данн^хм [Dunlop, 2002a] Fig. 5. Day-Dunlop plot (Mrs/Ms versus Bcr/Bc) showing the domain structure of magnetic minerals and hysteresis loops for each type of samples [Day et al., 1977; Dunlop, 2002a, 2002b] SD - single-domain grains, PSD - pseudosingle-domain grains, MD - multidomain grains, SP - superparamagnetic grains. Mixing curves SD + MD and
Ключевые слова
петромагнетизм,
магнитная восприимчивость,
археология,
палеокострищаАвторы
Кулакова Екатерина Петровна | Институт физики Земли имени О.Ю. Шмидта РАН | научный сотрудник, Лаборатория главного геомагнитного поля и петромагнетизма | ek.kula@yandex.ru |
Дедов Игорь Евгеньевич | Институт археологии и этнографии СО РАН | лаборант | 11.dedov.com@gmail.com |
Мещерякова Ольга Андреевна | Институт физики Земли имени О.Ю. Шмидта РАН | ведущий инженер, Лаборатория главного геомагнитного поля и петромагнетизма | oliya@ifz.ru |
Курбанов Реджеп Нурмурадович | Московский государственный университет имени М.В. Ломоносова; Институт географии РАН | кандидат географических наук, старший научный сотрудник, научно-исследовательская лаборатория новейших отложений и палеогеографии плейстоцена; отдел палеогеографии четвертичного периода | roger.kurbanov@gmail.com |
Всего: 4
Ссылки
Zan J., Fang X., Yan M., Li B. New insights into the palaeoclimatic interpretation of the temperature dependence of the magnetic susceptibility and magnetization of Mid-Late Pleistocene loess/palaeosols in Central Asia and the Chinese Loess Plateau // Geophysical Journal International. 2017. V. 208 (2). P. 663-673. DOI: 10.1093/gji/ggw419
Zan J., Fang X., Nie J., Teng X., Yang S. Rock magnetism in loess from the middle Tian Shan: Implications for paleoenvironmental interpretations of magnetic properties of loess deposits in Central Asia // Geochemistry, Geophysics, Geosystems. 2012. V. 13. Q10Z50. DOI: 10.1029/2012GC004251
Vidic N.J., TenPas J.D., Verosub K.L., Singer M.J. Separation of pedogenic and lithogenic components of magnetic susceptibility in the Chinese loess/paleosol sequence as determined by the CBD procedure and a mixing analysis // Geophys. J. Int. 2000. V. 142. P. 551-562.
Tite M.S., Mullins C. Enhancement of the magnetic susceptibility of soils on archaeological sites // Archaeometry. 1971. V. 13 (2). P. 209-219. DOI: 10.1111/j.1475-4754.1971.tb00043.x
Sun W.W., Banerjee S.K., Hunt C.P. The role of maghemite in theenhancement of magnetic signal in the Chinese loess-paleosol sequence: An extensive rock magnetic study combined with citrate-bicarbonatedithionite treatment // Earth and Planetary Science Letters. 1995. V. 133. P. 493-505.
Song Y., Li Y., Wang Q., Dong H., Zhang Z., Orozbaev R. Effect of chemical pretreatments on magnetic susceptibility of loess from Central Asia and the Chinese Loess Plateau // Royal Society of Chemistry. 2018. V. 8. P. 11087-11094.
Peters C., Thompson R., Harrison A., Church M.J. Low temperature magnetic characterisation of fire ash residues // Physics and Chemistry of the Earth. 2002. V. 27(25-31). P. 1355-1361. DOI: 10.1016/S1474-7065(02)00133-X
Özdemir Ö., Banerjee S.K. High temperature stability of maghemite (γ-Fe2O3) // Geophysical Research Letters. 1984. V. 11 (3). P. 161-164. DOI: 10.1029/GL011i003p00161
Peters C., Church M.J., Mitchell C. Investigation of fire ash residues using mineral magnetism // Archaeological Prospection. 2001. V. 8 (4). P. 227-237. doi: 10.1002/arp.171
Mullins C.E. The magnetic properties of the soil and their application to archaeological prospecting // Archaeo-Physika. 1974. V. 5. P. 143-347.
Morinaga H., Inokuchi H., Yamashita H., Ono A., Inada T. Magnetic Detection of Heated Soils at Paleolithic Sites in Japan // Geoarchaeology - An International Journal. 1999. V. 14 (5). P. 377-399.
Maxbauer D.P., Feinberg J.M., Fox D.L. MAX UnMix: A web application for unmixing magnetic coercivity distributions // Computers and Geosciences. 2016. V. 95. P. 140-145. DOI: 10.1016/j.cageo.2016.07.009
Maki D., Homburg J.A., Brosowske S.D. Thermally activated mineralogical transformations in archaeological hearths: Inversion from maghemite γFe2O4 phase to haematite αFe2O4 form // Archaeological Prospection. 2006. V. 13 (3). P. 207-227. DOI:10.1002/arp.277
Liu Q., Deng C., Yu Y., Torrent J., Jackson M.J., Banerjee S.K., Zhu R. Temperature dependence of magnetic susceptibility in an argon environment: Implications for pedogenesis of Chinese loess/palaeosols // Geophysical Journal International. 2005. V. 161 (1). P. 102-112. DOI: 10.1111/j.1365-246X.2005.02564.x
Le Borgne E. Influence du feu sur les propriétés magnétiques du sol et sur celles du schiste et du granite // Annales de Geophysique. 1960. V. 16. P. 159-195.
Linford N.T., Canti M.G. Geophysical evidence for fires in antiquity: Preliminary results from an experimental study // Paper given at the EGS XXIV General Assembly in The Hague, April 1999. Archaeological Prospection. 2001. V. 8 (4). P. 211-225. DOI:10.1002/arp.170
Le Borgne E. Susceptibilitee magnetique anomale du sol superficiel // Annale Geophysics. 1955. V. 11. P. 399-419.
Kruiver P.P., Dekkers M.J., Heslop D. Quantification of magnetic coercivity components by the analysis of acquisition curves of isothermal remanent magnetisation // Earth and Planetary Science Letters. 2001. V. 189 (3-4). P. 269-276. DOI: 10.1016/S0012- 821X(01)00367-3
Jrad A., Quesnel Y., Rochette P., Jallouli C., Khatib S., Boukbida H., Demory F. Magnetic Investigations of Buried Palaeohearths Inside a Palaeolithic Cave (Lazaret, Nice, France) // Archaeological Prospection. 2014. V. 21 (October 2013). P. 87-101. DOI: 10.1002/arp.1469
Ji J.F., Balsam W., Chen J., Liu L.W. Rapid and quantitative measurement of hematite and goethite in the Chinese loess-paleosol sequence by diffuse reflectance spectroscopy // Clays Clay Miner. 2002. V. 50. P. 208-216.
Heslop D., Dekkers M.J., Kruiver P.P., van Oorschot I.H.M. Analysis of isothermal remanent magnetization acquisition curves using the expectation-maximization algorithm // Geophysical Journal International. 2002. V. 148 (1). P. 58-64. DOI: 10.1046/j.0956- 540x.2001.01558.x
Hrouda F. Models of frequency-dependent susceptibility of rocks and soils revisited and broadened // Geophysical Journal International. 2011. V. 187 (3). P. 1259-1269. DOI: 10.1111/j.1365-246X.2011.05227.x
Jeleńska M., Hasso-Agopsowicz A., Kopcewicz B. Thermally induced transformation of magnetic minerals in soil based on rock magnetic study and Mössbauer analysis // Physics of the Earth and Planetary Interiors. 2010. V. 179 (3-4). P. 164-177. doi: 10.1016/j.pepi.2009.11.004
Fine P., Singer M.J., La Ven R., Verosub K., Southard R.J. Role of pedogenesis in distribution of magnetic susceptibility in two California chronosequences // Geoderma. 1989. V. 44 (4). P. 287-306. DOI: 10.1016/0016-7061(89)90037-2
Gibson T.H. Magnetic prospection on prehistoric sites in Western Canada // Geophysics. 1986. V. 51. P. 553-560. DOI:10.1190/1.1442109
Evans M.E., Heller F. Environmental Magnetism - Principles and Applications of Enviromagnetics. San Diego : Academic Press, 2003. 293 p.
Dunlop D.J., Ozdemir O. Rock Magnetism: Fundamentals and Frontiers. Cambridge University Press, 1997.
Egli R. Analysis of the field dependence of remanent magnetization curves // Journal of Geophysical Research. 2003. V. 108 (B2). P. 281. doi: 10.1029/2002JB002023
Dunlop D.J. Theory and application of the Day plot (Mrs/Ms versus Hcr/Hc) 1. Theoretical curves and tests using titanomagnetite data // Journal of Geophysical Research. 2002a. V. 107 (B3). 2056. DOI: 10.1029/2001JB000486
Dunlop D.J. Theory and application of the Day plot (Mrs/Ms versus Hcr/Hc) 2. Application to data for rocks, sediments, and soils // Journal of Geophysical Research. 2002b. V. 107 (B3). 2057. DOI: 10.1029/2001jb000487
Ding Z.L., Ranov V., Yang S.L., Finaev A., Han J.M., Wang G.A. The loess record in southern Tajikistan and correlation with Chinese loess // Earth and Planetary Science Letters. 2002. V. 200. P. 387-400.
Deng C.L., Zhu R.X., Jackson M.J., Verosub K.L., Singer M.J. Variability of the temperature-dependent susceptibility of the Holocene eolian deposits in the Chinese Loess Plateau: A pedogenesis indicator // Physics and Chemistry of the Earth. Part A: Solid Earth and Geodesy. 2001. V. 26. P. 873-878.
Day R., Fuller M., Schmidt V.A. Hysteresis properties of titanomagnetites: grain-size and compositional dependence // Physics of the Earth and Planetary Interiors. 1977. V. 13. P. 260-267.
de Boer C.B., Dekkers M.J., van Hoof T.A.M. Rock-magnetic properties of TRM carrying baked and molten rocks straddling burnt coal seams // Physics of the Earth and Planetary Interiors. 2001. V. 126 (1-2). P. 93-108. DOI: 10.1016/S0031-9201(01)00246-1
Dearing J.A., Dann R.J.L., Hay K., Lees J.A., Loveland P.J., Maher B.A., O’Grady K. Frequency-dependent susceptibility measurements of environmental materials // Geophysical Journal International. 1996. V. 124 (1). P. 228-240. doi: 10.1111/j.1365- 246X.1996.tb06366.x
Church M.J., Peters C., Batt C.M. Sourcing fire ash on archaeological sites in the Western and Northern Isles of Scotland, using mineral magnetism // Geoarchaeology. 2007. V. 22 (7). P. 747-774. DOI: 10.1002/gea.20185
Bonhomme T., Stanley J. Magnetic mapping of prehistoric Aboriginal fireplaces at Bunda Lake, Belarabon Station, New South Wales // Aust. Arch. 1986. V. 21. P. 63-73.
Carrancho Á., Villalaín J.J. Different mechanisms of magnetisation recorded in experimental fires : Archaeomagnetic implications // Earth and Planetary Science Letters. 2011. V. 312 (1-2). P. 176-187. DOI: 10.1016/j.epsl.2011.10.006
Bellomo R.V. A methodological approach for identifying archaeological evidence of fire resulting from human activities // Journal of Archaeological Science. 1993. V. 20 (5). P. 525-553. DOI: 10.1006/jasc.1993.1033
Aldeias V., Dibble H.L., Sandgathe D., Goldberg P., Mcpherron S.J.P. How heat alters underlying deposits and implications for archaeological fire features : A controlled experiment // Journal of Archaeological Science. 2016. V. 67. P. 64-79. DOI: 10.1016/j.jas.2016.01.016
Barbetti M. Traces of fire in the archaeological record, before one million years ago? // Journal of Human Evolution. 1986. V. 15 (8). P. 771-781. DOI: 10.1016/S0047-2484(86)80009-4
Aldeias V. Experimental approaches to archaeological fire features and their behavioral relevance // Current Anthropology. 2017. V. 58. S191-S205. DOI: 10.1086/691210
Шнайдер С.В. Туткаульская линия развития в мезолите западной части Центральной Азии : дис. ... канд. ист. наук. Новосибирск, 2015. С. 290.
Смекалова Т.Н., Восс О., Мельников А.В. Магнитная разведка в археологии. 12 лет применения Оверхаузеровского градиентометра GSM-19WG. СПб. : СПбГУ, 2007. С. 74.
Матасова Г.Г., Казанский А.Ю., Позднякова О.А. Опыт использования петромагнитного метода для оценки перспектив применения магниторазведки на территории археологических памятников Барабинской лесостепи // Физика Земли. 2016. № 6. С. 86-99. DOI: 10.7868/s0002333716060077
Додонов А.Е. Четвертичный период Средней Азии: Стратиграфия, корреляция, палеогеография. М. : ГЕОС, 2002. С. 250.
Матасова Г.Г., Казанский А.Ю., Позднякова О.А., Молодин В.И., Мыльникова Л.Н., Нестерова М.С., Кобелева Л.С. Итоги и перспективы применения петромагнитного метода для исследования археологических памятников Барабинской лесостепи // Проблемы археологии, этнографии, антропологии Сибири и сопредельных территорий. 2013. № 19. С. 251-254.