Гидротермальные извержения вызывают особый интерес у исследователей во всем мире [Moyle et al., 1990; Browne, Lawless, 2001; White, 1955; Eichelberger et al., 2020; Morgan et al., 2023; Rey-Devesa et al., 2023]. Обусловлено это тем, что гидротермальные извержения могут быть тесно связаны с рудообразованием, приводить к разрушениям сооружений геотермальных станций и трубопроводов; присутствие ряда гидротермальных минералов в выбросах позволяет получить информацию о субповерхностных температурах и типах гидротерм. Также при освоении ресурсов термальных полей современных действующих гидротермальных систем необходимо учитывать риски гидротермальных взрывов и сопутствующих явлений типа обвалов и оползней. Известны случаи, когда гидротермальные извержения приводили к гибели людей, как, например, трагические события на озерах Ниос и Манун (1984 и 1986 гг.) в Республике Камерун [Ionina, Kubeev, 2001]. Таким образом, прежде чем приступить к разработке новых гидротермальных месторождений, важно учитывать возможность и оценивать вероятность в этом районе таких событий.

Термальные проявления Южных Курильских островов изучены достаточно подробно, однако в научной литературе, посвященной Курильским островам, описание гидротермальных взрывов встречается редко и в основном связано с извержением в кратерных областях вулканических построек [Menyailov et al., 1969; Korzhinsky et al., 2004].

В 2022 г. исследователям стало известно о термальном озере Утиная Баня на о. Итуруп [Degterev et al., 2022]. В 2023 г. вышла статья, в которой высказывалось предположение о происхождении озера в результате серии гидротермальных взрывов [Khubaeva et al., 2023]. Спустя год, после проведения дополнительных полевых работ, у авторов появилось достаточно материала, который позволяет более уверенно говорить о генезисе оз. Утиная Баня.

Цель работы – определение генезиса уникального маарового термального оз. Утиная Баня по комплексу данных, оценка физико-химических характеристик озера и рисков гидротермальных взрывов.

Остров Итуруп является самым крупным островом Большой Курильской гряды. Состав и строение позднекайнозойских отложений Курильской дуги указывает на существование здесь интенсивных и длительных процессов вулканизма [Piskunov, 2004].

Остров относится к Большекурильской структурно-фациальной зоне с широко развитыми вулканическими и вулканогенно-осадочными образованиями неогенового и четвертичного возраста, которые прорываются многочисленными экструзивными, интрузивными, субвулканическими телами различного состава [State Geological Map…, 2002]. Наиболее широкое распространение на о. Итуруп получили породы неогенового возраста [Gorshkov, 1967]. Неогеновые отложения слагают фундамент четвертичных вулканических построек и обнажаются на поверхности. Самой древней в неогеновом разрезе (более 1200 м) острова является рыбаковская свита (средний – верхний миоцен). Породы представлены переслаивающимися пачками вулканогенно-осадочных и вулканокластических пород. На рыбаковской свите залегает фрегатский вулканический комплекс, сложенный эффузивно-пирокластическими образованиями среднего и основного состава (370 м). На Южных Курильских островах локально распространена камуйская свита, которая выполняет в основном плиоценовые палеодепрессии. Фрегатский вулканический комплекс и камуйская свита согласно перекрываются пемзовыми песками и песчаниками головинской свиты. Нижняя граница неогена в Большекурильской структурно-формационной зоне не установлена [Gladenkov et al., 1998]. Четвертичными образованиями сложены вулканические аппараты, отложения морских террас, делювиально-пролювиальные плащи под склонами гор.

Остров Итуруп имеет мозаичное строение. На острове выделяется серия продольных и поперечных разломов и разрывных нарушений (рис. 1), определяющих расположение термальных вод.

В центральной части острова выделяется два крупных дуговых разлома [Asaulov et al., 1982], с которыми генетически связаны термальные проявления и зоны формирования термальных вод. Первый разлом протягивается от с. Горячие Ключи до с. Рейдово. Второй расположен на 20–30 км южнее. С этим разломом связаны активные вулканы Иван Грозный, Мачеха, Тебенькова, Баранский (рис. 1).

В 2020–2021 гг. нами были проведены исследования в районе оз. Утиная Баня, полученные результаты позволили сделать предположение о том, что озеро имеет взрывное происхождение и представляет собой маар, образовавшийся в результате серии гидротермальных извержений [Khubaeva et al., 2023]. Однако полученных данных было недостаточно для того, чтобы делать выводы о происхождении озера, поэтому в 2022–2023 гг. были выполнены дополнительные исследования:

1) построена 3D модель тектонической раздробленности блока горных пород для центральной части о. Итуруп;

2) определен минеральный состав пород в шурфе 1;

3) проведена температурная съемка дна озера;

4) сделаны химические анализы свободных газов, выделяющихся со дна оз. Утиная Баня;

5) при помощи портативного газового анализатора ГАНК-4 НПО «ПРИБОР» ГАНК» (г. Москва, Россия) проведена газовая съемка над поверхностью озера;

6) получены дополнительные данные о химическом составе воды;

7) на основе данных эхолотной съемки [Khubaeva et al., 2023] и геометрического нивелирования построена 3D модель озерной котловины и вала вокруг нее.

Для идентификации разрывных нарушений на земной поверхности применялся широко распространенный в геологии метод дешифрирования космических снимков и топографических карт. При помощи метода воксельного моделирования была построена 3D модель тектонической раздробленности для центральной части о. Итуруп.

Термосъемка дна озера проводилась при помощи термопары и GPS-навигатора. Замеры температуры выполнялись по профилям с шагом 5–7 м. Всего было пройдено пять профилей вдоль озерной котловины и два профиля в крест.

Дешифрирование космоснимков. Метод дешифрирования спутниковых снимков заключается в распознавании объектов по их признакам, определении характеристик и установлении взаимосвязей с другими объектами. Использовался метод дешифрирования космических снимков разной детальности, современные разрывные тектонические нарушения выделялись с использованием методики Н.П. Костенко [Kostenko, Bryanceva, 2004].

Линеаментный анализ. На основе дешифрирования топографических карт масштаба 1:100000 и 1:50000, а также космических снимков различной детальности были построены схемы линеаментных сетей. Пространственное распределение, плотность и количество пересечений линеаментов на единицу площади являются ключевыми факторами для выделения высокопроницаемых зон, служащих путями миграции вещества [Taskin, 2017].

Воксельное моделирование (трехмерное моделирование тектонической раздробленности). Основой для построения 3D модели тектонической раздробленности послужили карты линеаментных сетей, полученных в результате дешифрирования спутниковых снимков и топографических карт масштаба 1:100000 и 1:50000. Цифровая обработка проводилась с использованием технологии геоинформационных систем (ГИС). Пространственная привязка и векторизация результатов дешифрирования выполнялись программами ArcMap и ArcView. По векторизованным данным создана единая цифровая карта линеаментов центральной части о. Итуруп.

Раздробленность геологической среды количественно оценивалась по величине удельной длины разломов (УДЛ) [Nechaev, 2010]. УДЛ представляется в виде трехмерной матрицы, из которой формируется объемное (voxel) изображение раздробленности исследованного блока. Тектоническая раздробленность кубического блока определяется отношением объема всех трещин внутри блока к его объему. Объем трещин определяется как сумма произведений длины (l), ширины раскрытия (b) и глубины проникновения (h) отдельных трещин. Для блока кубической формы с ребром а, в котором одной из граней является дневная поверхность, коэффициент тектонической раздробленности будет равен Σl·b·h/a3. По фотоизображению невозможно определить величины b и h, но для конкретной локализованной территории ширину трещин b можно считать постоянной и в расчете УДЛ не учитывать. Также, исходя из эмпирически установленного правила, что трещиноватость одной грани кубика породы адекватно отражает степень трещиноватости всего образца [Nechaev, 2010], и основываясь на принципах подобия (физические явления и объекты, обладающие подобием геометрии и физических переменных, подобны, деформация макроскопического уровня любого объема земной коры может быть рассмотрена как сумма составляющих его деформаций многочисленных объемов), можно аппроксимировать это правило и для отдельных блоков земной коры. Тогда глубину проникновения разломов h можно принять равной ребру а. Таким образом, определяя степень раздробленности как значение удельной длины линеаментов, равной Σl/а2 верхней грани блока, и увеличивая последовательно размер элементарной ячейки (куба), можно проследить раздробленность на глубину. Вычисленные значения УДЛ относятся к центру куба на глубине а/2.

Построение модели проводилось в несколько этапов. Сначала на карту линеаментов накладывалась квадратная сетка с размером ячейки 1 км2, в пределах каждой клетки подсчитывалась суммарная длина линеаментов, нормировалась на площадь клетки, полученные значения относились к центрам кубов. Ребро каждого куба составляло 1 км. Так была составлена матрица УДЛ для глубины 0.5 км. За нулевую отметку принимался уровень моря. На следующем этапе расчетов размер ребра элементарного куба увеличивался на 200 м, и вся процедура повторялась. Полученные данные по всем глубинным уровням соединены в одну 3D матрицу (Х, У, Z, УДЛ), которая использована для объемной визуализации. Удельная длина линеаментов выражена в км–1.

Метод построения 3D модели тектоничекой раздробленности блока пород подробно описан в работах [Nechaev, 2010; Taskin, Sidorov, 2014, 2015; Taskin, 2017; Sidorov, Taskin, 2017; Khubaeva et al., 2020, 2022].

Отбор образцов воды, газа. В период полевых работ в 2022 и 2023 гг. были собраны образцы термальных вод озера и источников (рис. 2).

Неконденсирующиеся газы анализировали хроматографическим методом, на базе газового хроматографа «Кристалл-5000».

Образцы воды отбирали в чистые пластиковые бутылки объемом 1 л, а для анализа на микроэлементный состав пробы отбирали отдельно в пластиковые пробирки объемом 50 мл, подкисляли для консервации азотной кислотой квалификации «осч».

Анализ водных образцов. Анализ микрокомпонентного состава проводили на квадрупольном ICP-MS масс-спектрометре Agilent 7500 c, допустимая ошибка составляла 10 %, пределы обнаружения для различных элементов варьировались от 0.0001 до 241 ppb с медианным значением 0.003 ppb (мг/т, мкг/кг). Анализы были выполнены в ЦКП «Микроанализ» при Лимнологическом институте СО РАН г. Иркутск (исполнитель к.г.-м.н. О.Г. Степанова, лаборатория междисциплинарных эколого-экономических исследований и технологий).

Определение макрокомпонентного состава водных образцов было проведено в Аналитическом центре Института вулканологии и сейсмологии ДВО РАН. Концентрации кальция, магния, натрия, калия, лития определяли методом атомно-абсорбционной и пламенно-эмиссионной спектрометрии, с помощью атомно-абсорбционного спектрофотометра Thermo Electron SOLAAR, методики ПНДФ 14:1:2:4.137-98 и 14:1:2:4.138-98.

Алюминий определяли фотометрически, по методике ПНДФ 14:1:2:4.166-2000. Общее содержание железа также определяли фотометрическим методом, по ПНДФ 14:1:2:4.50-96, измерения проводились с помощью спектрофотометра UV-Mini-1240 Shimadzu.

Определение pH было выполнено потенциометрически с помощью pH-метра HANNA Instruments HI 2211, по методике ПНДФ 14:1:2:3:4.121-97. Содержание фторид-ионов определено потенциометрическим методом с неопределенностью 18 %, методика соответствует ПНДФ 14:1:2:4.270-2012.

Концентрация кремнекислоты определена фотометрически, по ПНДФ 14:1:2:4:.215-06. Аммоний определен фотометрически, по методике ПНДФ 14:1:2:3.1-95. Содержание хлоридов определяли аргентометрически, по методике ПНДФ 14:1:2:3.96-97. Концентрация сульфат-иона определена турбидиметрическим методом, по образованию взвеси осадка сульфата бария, по методике ПНДФ 14:1:2:159-2000. Погрешности методик составляли 10–15 % от полученного значения, погрешность pH-метра не более 0.2 единицы pH.

Определение минерального состава. Для исследования минерального состава отложений была использована рентгеновская дифрактометрия. Дифрактограммы записаны с помощью порошкового рентгеновского дифрактометра XRD 7000 MAX (Shimadzu) в диапазоне 6–65° 2θ с шагом 0.1° 2θ, скорость сканирования 4 град/мин. Образцы растирались в агатовой ступке, съемка велась при вращении образца. Исследование проведено на базе Аналитического центра ИВиС ДВО РАН.

Для центральной части о. Итуруп была построена 3D модель тектонической раздробленности с целью уточнения проницаемых зон верхней части земной коры. В разное время в этом районе исследователями [Pchelkin et al., 1986, 1993; Talovina et al., 2022] выделялась серия крупных разрывных нарушений (см. рис. 1) СВ простирания, которые также прослеживаются на 3D модели тектонической раздробленности в виде зон с высокой тектонической деструкцией (рис. 3). Высокая плотность трещин в вулканических районах часто обнаруживается над внедренными магматическими телами. На границах таких зон нередко существуют области перехода участков питания метеорных и глубинных вод и, как следствие, на поверхность выходят термальные проявления (рис. 3).

Озеро Утиная Баня расположено на границе зоны с высокой плотностью трещин (рис. 3), трассирующей крупное дуговое разрывное нарушение (см. рис. 1). На вертикальном разрезе 3D модели (линия А–Б) (рис. 3, 4), в районе дугового разрывного нарушения (см. рис. 1), а также под оз. Утиная Баня, на глубине до 2000 м, отмечаются области аномально высокой плотности трещин, которые сливаются между собой, образуя подобие подводящих каналов, что может свидетельствовать о наличии интрузивного тела на глубине [Taskin, Sidorov, 2014; Taskin, 2016; Khubaeva et al., 2020, 2022], выполняющего роль источника теплового питания термальных выходов на дне озера. Однако такое предположение требует проверки.

Также данные эхолотной съемки котловины оз. Утиная Баня [Degterev et al., 2022; Khubaeva et al., 2023] указывают на наличие двух воронок в восточной и западной части озера. По результатам проведения термосъемки дна озера в 2023 г. была построена карта распределения температур (рис. 5). Наиболее высокие температурные значения совпадают с зонами воронкообразных углублений на дне озера [Khubaeva et al., 2023].

На рис. 5 проявлены два наиболее прогретых участка озерной котловины. Первый участок совпадает с наиболее глубокой частью озера. По данным эхолотной съемки в этом месте зафиксировано V-образное углубление, а также прослеживается акустическая брешь от подводного газогидротермального выхода [Khubaeva et al., 2023]. На поверхности воды наблюдаются выходы газа. Второй участок расположен в ЮЗ части озерной котловины и сопровождается бурлением. Температура на этих участках в зимний период колеблется от 17.0 до 18.7 ℃ (см. рис. 2). Средняя температура воды в озере в зимний период составляет около 15 ℃.

По химическому составу и по физико-химическим характеристикам (рис. 6; табл. 1) воды оз. Утиная Баня похожи на источники, расположенные в районах активных вулканов о. Итуруп: Старозаводское сольфатарное поле (влк. Баранского), Нагорные источники (влк. Тебенькова) и источник Лопастной (влк. Иван Грозный) (рис. 6) [Markhinin, Stratula, 1977; Zharkov, 2014; Pchelkin et al., 1986].

Обобщенный ионный состав воды озера, выраженный через формулу Курлова:

вода хлоридно-сульфатная натрий-кальциевая, кислая, основной растворенной солью является сульфат кальция. Минерализация воды составляет около 410 мг/л, а молочно-голубой цвет обусловлен наличием в толще взвешенных частиц. Во время полевых работ авторами наблюдалась биогенная сера в виде длинных волосоподобных колоний микроорганизмов.

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

Отметим повышенные относительно фоновых концентрации бора, мышьяка, цинка, лития и фосфора, бария и стронция, эти элементы характерны для термальных источников современных гидротермальных систем [Taran, Kalacheva, 2023; Kalacheva et al., 2016].

Кислые сульфатные кальциевые термы неоднократно наблюдались на юге Камчатки и на Курильских островах, и приурочены они как раз к активным вулканам, к которым прилегают современные гидротермальные системы [Taran, Kalacheva, 2023; Kalacheva et al., 2016]. Отложения гипса, бассанита и ангидрита, которые являются сульфатами кальция разной степени гидратированности, обнаруживаются на термальных полях современных гидротермальных систем Камчатки [Kiryukhin et al., 2023] и свидетельствуют о наличии сульфатных кальциевых вод.

По совокупности макро- и микрокомпонентного состава, а именно по низкому pH, высокому содержанию сульфата кальция в водах, заметному присутствию стронция и бария, мышьяка и цинка, можно сделать заключение о сходстве воды оз. Утиная Баня и вод других гидротермальных систем [Abrosimova, 2013].

В связи с этим мы полагаем, что воды озера подпитываются глубинными гидротермальными флюидами, с которыми сопряжена опасность гидротермальных взрывов при изменении вулканической активности. Воды озера кислые, в том числе за счет окисления поступающей с глубины серы, агрессивны к вмещающим породам и способствуют выщелачиванию ряда элементов (табл. 1). Поскольку оз. Утиная Баня представляет интерес не только как научный объект, но и как часть уникальных туристических, бальнеологических, стратегических ресурсов Дальнего Востока, необходимо оценивать риски повторения гидротермальных взрывов в перспективе. Для этого рассмотрим геотермометрию озера по химическому составу.

Применение классического Na/K-геотермометра, рассчитанного для равновесия полевых шпатов с водой, несмотря на удобство, наталкивается на сложности, связанные с влиянием вторичных катионообменных минералов на состав вод. На классификационном треугольнике Гиггенбаха воды оз. Утиная Баня лежат в области «магниевого угла» и относительно равновесий воды с полевыми шпатами, слюдами и хлоритом, для которых рассчитаны граничные линии [Giggenbach, 1988], являются незрелыми. Характерно то, что в областях современного вулканизма воды могут быть равновесны с вторичными катионообменными минералами, а относительно равновесий с полевыми шпатами, слюдами и хлоритами имеют «незрелый» облик. Вследствие этого в нашем случае классический Na/K-геотермометр можно использовать лишь формально, оценка температуры формирования состава вод с его помощью дает значения в пределах 230–290 °C, а для прибрежного источника UB3-1/22 (см. рис. 2) – в пределах 263–316 °C.

Высокие температуры отвечают низким значениям соотношения Na/K, которые характерны для вод озера вообще. В свою очередь, низкие соотношения Na/K могут быть связаны как с высокой температурой при формировании состава вод, так и с вторичной минерализацией, например с формированием толщ гидротермальных глин, в состав которых входит слоистый силикат монтмориллонит. Для монтмориллонита, который является эффективным катионообменником, характерны как раз невысокие соотношения Na/K в контактном растворе, и в присутствии смектита применение формулы для классического полевошпатового геотермометра дает завышенные значения, поэтому соотношение Na/K озерных вод обусловлено или контактом с монтмориллонитом, или повышенными температурами в недрах, возможны и обе причины. Давление водяного пара при этих температурах лежит в диапазоне порядка 25–75 атм, что в перспективе чревато разрушением верхнего водоупора и гидротермальным взрывом.

На рис. 7 показаны дифрактограммы отложений, взятых в ходе проходки шурфа 1 (см. рис. 2; рис. 7), с наложенными эталонными данными для дешифровки. Верхний горизонт представлен обломками пород, в составе которых определены типичные породообразующие минералы, а именно полевые шпаты ряда альбит – анортит с составом, сдвинутым в сторону альбита, и кварц, что отвечает скорее кислым исходным породам. Нижележащие слои представлены слабо окристаллизованным ромбическим гетитом FeO(OH), магнетитом Fe3O4 и гидратированным гематитом Fe2O3(H2O)n, при этом доминирует гетит, а гематит и магнетит появляются в основании разреза. Сечение второго шурфа, расположенного у основания вала, представлено суглинками [Khubaeva et al., 2023].

О возможном превращении гетита в гематит и магнетит свидетельствует размытость рефлексов гематита (рис. 7), в то время как обычно гематит из исходных пород характеризуется узкими четкими рефлексами. Размытость рефлексов может объясняться малой областью когерентного рассеяния, невысокой степенью кристалличности гематита, поэтому весьма вероятно, что гетит претерпел преобразование при локальном нагреве и трансформировался в гематит. Отметим, что промежуточной фазой такого превращения является гидрогематит или протогематит, который по структуре близок гематиту [Gualtieri, Venturelli, 1999; Wolska, Schwertmann, 1989]. Стратиграфические данные Вайотапу [Hedenquist, Henley, 1985] показывают, что гидротермальные извержения развивались снизу вверх в результате локальных изменений давлений в пределах верхних 300 м гидротермальной системы. В других системах (Каверау, Новая Зеландия) [Nairn, Solia, 1980] определяются аналогичные глубины.

Для сравнения, в кальдере Узон, где расположены многочисленные термальные источники, есть оз. Банное, расположенное неподалеку от кордона и по всем признакам имеющее мааровое происхождение. Состав его вод:

соотношение Na/K низкое, около 3.6, а температура, рассчитанная по геотермометру для полевых шпатов, порядка 280 °C. Но оз. Банное в кальдере Узон обрамлено толщей глинистых отложений, которые могут искажать значения температуры, полученные с помощью классического геотермометра, за счет присутствия монтмориллонита завышать их. Явных толщ аргиллизированных пород, глинистых плащей вокруг оз. Утиная Баня не найдено, что не исключает присутствия монтмориллонита в почвах и некоторого его влияния на состав вод [Sergeeva et al., 2023].

Воронкообразный профиль дна озера [Khubaeva et al., 2023], перевернутый/нарушенный порядок отложений на бортах озера, когда осадочные озерные отложения оказались погребены под слоем породообразующих минералов (рис. 7), свидетельствуют о том, что гидротермальные взрывы в истории озера случались. Характерные для гидротермальных взрывов температуры, составляющие сотни градусов, после взрыва вызвали прогрев гетитовых отложений и частичное превращение гетита в магнетит и гематит. Отметим восстановленный характер флюидов, поскольку имел место не только прогрев, но и переход гематита в магнетит.

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

Свободно выделяющиеся со дна озера газы по составу близки к составу атмосферы, но сдвинуты в сторону азот-аргоновой грани, а не в сторону тех составов воздуха, которые характерны для его водных растворов при разных температурах. Для сравнения на график (рис. 8) нанесены данные из работы [Kalacheva et al., 2022], и видно, что озерные газы являются результатом смешения глубинной вулканогидротермальной компоненты и воздуха атмосферы, так как составы газов источников лежат между составом глубинных газов и составом воздуха. Состав свободно выделяющихся газов указывает скорее на глубинный источник флюидов и глубинное тепловое питание. Содержание CO2 в свободно выделяющихся газах оз. Утиная Баня составляет порядка 0.16 об. %, а на источниках мыса Конакова не менее 5.7 об. %.

Полевая съемка газов, выделяющихся со дна озера, показала довольно заметные содержания метана и зафиксировала следы присутствия диметилсульфоксида (ДМСО C2H6SO), диметилсульфида (ДМС C2H6S), сероуглерода (CS2). Концентрации метана лежат в пределах 0.01–0.03 об. %, что совпадает с данными по свободным газам мыса Конакова [Kalacheva et al., 2022]. Метан, ДМС могут быть как биогенного, так и глубинного происхождения, а ДМСО является продуктом окисления ДМС. Концентрации диметилсульфоксида и диметилсульфида лежали в пределах 0.15–0.30 и 0.90–3.30 мг/м3 соответственно, сероуглерода – 0.08–0.50 мг/м3. Мониторинг малых компонент газов оз. Утиная Баня представляет интерес, в частности, для выяснения сезонных колебаний концентраций.

Формирование осадочных и взрывных отложений. Верхний горизонт около озера состоит из плотных агрегатов гетита [Khubaeva et al., 2023]. Слабо окристаллизованный гетит образовался при осаждении соединений железа из вод озера, при подпитке термальными источниками. Гетит – это типичная форма, в которой кристаллизуется окисленное железо в условиях разбавленных растворов в слабокислых – близнейтральных средах [Krauskopf, 1963], а эти условия отвечают исследуемому озеру, так как растворимость гетита значительно ниже растворимости сульфатов и хлоридов железа. Присутствие гематита и магнетита в нижней части гетитовых отложений может свидетельствовать о локальных прогревах основания толщи отложений, так как дегидратация, то есть потеря воды, идет при нагревании, а гематит является именно продуктом дегидратации гидроокислов железа и гетита.

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

Похожий порядок отложений наблюдался на Мутновском геотермальном месторождении (рис. 9) после гидротермального взрыва, произошедшего весной 2024 г.

В результате взрыва образовалась интенсивно парящая воронка с размерами примерно 200×140 м, глубиной порядка 30 м. Взрыв произошел на месте локализации скважин, через которые под давлением разгружалась горячая пароводяная смесь. Ранее длительная разгрузка пароводяной смеси привела к формированию вокруг скважин отложений аморфного кремнезема, который образовывался из парогазовой смеси, и было точно известно его происхождение. После гидротермального взрыва практически весь аморфный кремнезем оказался погребенным под рыхлым, вынесенным взрывом материалом. Проводя аналогию, мы считаем, что породообразующие минералы, залегающие над осадочным гетитом в районе оз. Утиная Баня, вынесены гидротермальным взрывом.

На основе полученных нивелирных профилей и данных батиметрической съемки [Degterev et al., 2022; Khubaeva et al., 2023] была построена 3D модель оз. Утиная Баня (рис. 10), на которой отчетливо виден вал, окружающий озеро, а также углубление на дне озера. По данным батиметрической съемки дно озера представляет собой две воронки с выходом газа к поверхности [Khubaeva et al., 2023]. Все это может свидетельствовать в пользу взрывного происхождения озера.

Следует отметить, что питание озера осуществляется за счет грунтовых вод и атмосферных осадков. Вал, обрамляющий озерную котловину, ограничен сопками, представляющими собой небольшую котловину, открытую на СЗ (см. рис. 2). Из озера вытекает ручей с расходом воды в зимний период около 1.9–2.0 м3/с. Площадь водосбора озера составляет 3600 м3. Сезонно понижение между озерным валом и бортами окружающих озеро сопок заполняется атмосферными осадками, формирующими временный водоток, питающий озеро (см. рис. 2). Не исключено, что были периоды, когда уровень воды поднимался настолько высоко, что озерный вал оказывался под водой. Этим можно объяснить наличие в верхних слоях шурфа (см. рис. 7) темно-зеленых суглинков небольшой мощности (до 20 см), которые могут интерпретироваться как мелководное дно, этим же можно объяснить осаждение гетита.

Концептуальная модель гидротермальной системы оз. Утиная Баня. По линии А–Б (рис. 11) построен профиль, на основе которого была составлена схема гидротермально-магматической системы: р. Аргунь – крупное разрывное нарушение СВ простирания – оз. Утиная Баня – разрывное нарушение СВ простирания, восходящий поток гидротерм (рис. 11).

Температура. Выделенные на схеме глубинного строения оз. Утиная Баня изотермы частично основаны на данных замеров дна озерной котловины, максимальные температуры здесь составляют 19–25 ℃, а также на данных, полученных в результате применения Na/K геотермометров. Расчеты по геотермометрам указывают на температуру 230–290 ℃, глубина этих изотерм неизвестна, поэтому авторы использовали данные поисковой скважины 2к, расположенной в р-не пос. Рейдово, в 6 км от оз. Утиная Баня [Pchelkin et al., 1986].

По результатам термометрии геотермический градиент в скважине составляет 5.48 ℃ на 100 м. Температура в устье скважины в 1985 г. составляла 21 ℃, тогда на глубине около 1200 м температура должна составлять около 252 ℃.

Воронки гидротермального взрыва. По данным эхолотной съемки, на дне оз. Утиная Баня выделяется два воронкообразных углубления, сопряженных с зонами разгрузок на поверхность термальных вод.

Малоглубинные гидротермальные взрывы являются обычными на активных геотермальных системах. Такие извержения могут быть обусловлены изменениями локального теплового потока в тех местах, где миграция потока пара изменялась в связи с отложениями минералов, или там, где происходили изменения субповерхностных гидрогеологических условий [Hedenquist, Henley, 1985; White, 1955].

Сейсмические события, способствующие формированию новых и открытию уже залеченных трещин или образованию усадочных трещин, могут привести к внезапному снижению давления и резкой смене агрегатного состояния вещества в гидротермальной системе с жидкого на газообразное и, как следствие, к гидротермальному взрыву. На возможность гидротермального взрыва указывает также присутствие газов (H2, CH4, СО2) в гидротермальной системе оз. Утиная Баня, наиболее вероятно имеющих вулканическое происхождение.

Разрывные нарушения и магматические очаги. Разрывные нарушения и расположение магматических очагов выделены на основе тектонической схемы, выполненной В.И. Пчелкиным и др. в 1986 г. [Pchelkin et al., 1986], а также по данным 3D модели тектонической раздробленности блока пород для центральной части о. Итуруп. Глубина магматических очагов не установлена.

В районе центральной части о. Итуруп выделяется крупный дуговой разлом (см. рис. 1), который тянется от с. Горячие Ключи до с. Рейдово. С этим разломом генетически связаны зоны формирования термальных вод (Горячие Ключи, Многоозерные, Глухие, Жаркие Воды). Вдоль разлома наблюдается большое количество интрузивных тел и термальных источников (см. рис. 1). Разлом проходит в стороне от современных вулканических систем. Очевидно то, что оз. Утиная Баня географически попадает в зону влияния этого разлома.

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

Давление в гидротермальной системе. По данным классического Na/K геотермометра определяется высокая температура (230–290 ℃) формирования состава вод, свидетельствующая о том, что давление водяного пара при этих температурах лежит в диапазоне порядка 25–75 атм, что может привести к разрушению верхнего водоупора.

Состав газов, выделяющихся со дна озера. Озерные газы являются результатом смешения глубинной вулканогидротермальной компоненты и воздуха атмосферы (см. рис. 8). Состав свободно выделяющихся газов указывает на глубинный источник флюидов и глубинное тепловое питание.

Полученные в ходе исследований данные говорят о том, что термальное оз. Утиная Баня образовалось в результате серии гидротермальных взрывов, о чем свидетельствует состав воды и свободно выделяющихся газов; высокая степень тектонической раздробленности вокруг озера; высокие температуры на глубине (по данным Na/K геотермометров); наличие термальных источников на дне озера; вал, окружающий озерную котловину; перевернутый/нарушенный порядок отложений на бортах озера; наличие слоя брекчированной породы в разрезе шурфа 1.