ОБНАРУЖЕНИЕ МИКРОБНЫХ КОМПОНЕНТОВ В ВЫПАДАЮЩИХ ОСАДКАХ ЭНЦЕЛАДА

АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ, 2020, том 75, № 2, с. 188-199

УДК 523.46-87*2

ОБНАРУЖЕНИЕ МИКРОБНЫХ КОМПОНЕНТОВ В ВЫПАДАЮЩИХ

ОСАДКАХ ЭНЦЕЛАДА

(© 2020 Д. Котларж1*, У. Зеленкевич2, Н. Е. Залевска1, К. А Кубяк1

Поступила в редакцию 3 июня 2019 года; после доработки 19 января 2020 года; принята к публикации 19 января 2020 года

Благодаря наличию жидкой воды и доступной энергии Энцелад является привлекательной целью для поиска признаков жизни. В недавнем исследовании было доказано, что в выбросах гейзеров вблизи южного полюса Энцелада присутствуют водяной пар, водяной лед, углекислый газ, метан и молекулярный водород. Возможное сходство между физическими и химическими условиями на дне океана Энцелада и карбонатным минеральным раствором интенсивных подводных источников в «трубах» гидротермального поля Потерянного Города позволяет создать математическую модель процесса продвижения микробов вверх сквозь ледяную оболочку. В настоящей работе мы представляем первые результаты кинетических симуляций по методу «частица-в-клетке» для вычисления расстояний микробов в океане глубиной 10 км. Мы получили результаты для распределения микробного компонента — аналога Ме1И.апо$агс1па1е$$р. — характеризуемого массой 6.6 пг и диаметром 2.0 мкм в струйных выбросах вещества Энцелада. Мы предположили, что плотность теплового потока составляет 0.1 Вт м-2, толщина ледяной коры — 5 км, а концентрация клеток около дна океана — 105 клеток/мл. Мы подтвердили предположение исследовательской команды Порко о концентрации клеток около поверхности океана примерно в 104 клеток/мл и многообразии вертикальной плотности струй. Мы нашли, что оптимальная высота для детектирования микробного компонента составляет менее 1.0 км и что локальные измерения, сделанные ранее масс-спектрометром Кассини и предложенные для орбитальной миссии Энцелада, возможно, неэффективны на высоте 50 км.

Энцелад является привлекательной средой во внешней Солнечной Системе, где жизнь может появиться и процветать жизнь, благодаря жидкой воде и доступности энергии (McKay et al. 2008; Lunine et al. 2015). В последние годы были предложены, проанализированы и обсуждены концепции миссии для поиска жизни (Cable et al. 2017a;b).

Особенный интерес представляет то, что такой маленький астрономический объект, диаметром менее 500 км, термически активен. Источник этой активности до сих пор неизвестен. Рассматривается приливное нагревание от Сатурна; оно является следствием орбитального резонанса Эн-целада и еще одного спутника Сатурна, Дионы. В результате этих взаимодействий ледяная поверхность Энцелада трескается и на ней формируются проломы, которые, как в случае трещин под названием «Тигровые полосы», начинают проявлять

*E-mail: Jan.Kotlarz@ilot.edu.pl

активность при контакте с вакуумом. Теплые внутренние слои Энцелада под действием приливных процессов и радиоактивного распада прогревают внешний слой льда, что приводит к образованию «Тигровых полос», а также выделению водяного пара с кристаллами льда. В настоящий момент по причине активности Энцелада проявление крио-вулканизма только на южном полюсе этого спутника остается без объяснения.

Под поверхностью ледяной коры находится океан глубиной 10—50 км, являющийся источником вещества, которое выбрасывается в космос через «Тигровые полосы». Активность струй увеличивается в апоцентре орбиты Энцелада. Более точные наблюдения Кассини показали, что эти структуры удивительно схожи с океаническими хребтами на Земле. Помимо характерной «ступенчатой» структуры, наблюдалось нечто вроде зоны распространения и трансформации дефектов. Можно предположить, что «растрескивание» Энцелада связано с существованием гидротермальных расщелин, как и в случае земных разломов.

Недавние исследования доказали, что выбрасываемое в струях вещество формирует нечто вроде заслона. Похоже, что выбросы непрерывны и создают вокруг Энцелада плотный туман из мелкой ледяной пыли, которая поставляет материал для одного из колец Сатурна (кольцо E).

«Тигровые полосы» образованы 5 щелями, расположенными на расстоянии 35 км друг от друга, а толщина ледяной коры под ними составляет всего 5 км (Yin et al. 2016). Ученые предполагают, что соленый океан Энцелада (-Na, -Cl, -CO3) является щелочной средой с pH от 11 до 12 (Glein et al. 2015). Это позволяет нам сделать важные выводы о геохимических процессах внутри Энцелада. Такой высокий pH может быть следствием серпенти-низации оливина. В результате серпентинизации, помимо прочих, образуются такие вещества как H2 и CH4 (Sleep et al. 2004), которые считаются подходящей средой для возникновения жизни.

Глубина океана Энцелада по оценкам составляет 10—50 км (less et al. 2014; Davila et al. 2018). Любые микроорганизмы, живущие на его дне, должны будут преодолеть это расстояние, прежде чем начнут подниматься через линейные провалы «Тигровых полос» в южной полярной области этого спутника Сатурна. Основные параметры струй, такие как средний темп выработки льда, соотношение твердого вещества и пара, химический состав, градиент температуры, зависят от процессов подледяной корой. Зависимость других, например, среднего геотермического потока в море района южного полюса Энцелада, pH океана, содержание соли, толщина ледяной коры и океана, от градиентов их давления и температуры все еще обсуждается.

Возможное сходство физических и химических условий на дне океана Энцелада с карбонатным минеральным раствором интенсивных подводных источников в «колоннах» гидротермального поля Потерянного Города (LCHF), находящегося на 15 км западнее Срединно-Атлантического Хребта на вершине массива Атлантис на глубине примерно 750 м, дает шанс найти микробную жизнь и создать математическую модель процессов развития микробов и их подъема к поверхности через ледяную кору. Давила et al. (2018) рекомендовали включить разработку моделей формирования струй и выбросов в план будущей Астробиологической Стратегии NASA (Wright et al. 2019). Среда обитания подповерхностного океана Энцелада является объектом изучения Авиационного института Варшавы (Польша). В нашей работе представлены первые результаты кинетического («частица-в-клетке») моделирования параметров океана («мелководная» версия, около 10.0 км толщиной) и пути микробов — описанных как особый вид частиц — к поверхности океана.

Недавние исследования микробных семейств ЬСНР показали, что пористые внутренние стены колонны с наиболее интенсивным подводным источником населены простейшей биопленкой, в которой на 80-90% доминирует единственный филотип порядка Ме1Напо8агста1е8— МеШапо8агс1па1е8 Потерянного Города (Вга2е11;оп е! а1. 2010). Эта биопленка содержит 2.0 х 106 и 3.1 х 108 клеток/мл в областях, промываемых гидротермальными жидкостями (БИгепк 2017). Многочисленные Ме1Напо8агста1е8 Потерянного Города адаптировались к температуре примерно 90°С, рН 9—11 и жидкости источников с растворенным водородом (Н2) и метаном (СН4) (Вга2еИоп е! а1. 2010). Н2 и СН4 являются результатом серпентинизации. Количество СН4 довольно постоянно, а доля Н2 меняется, вероятно, по причине биологического окисления водорода (Pгoskuгowski е! а1. 2006). Физиология МеШапо8агс1па1в8 Потерянного Города неизвестна, однако анализ ряда клонов микробных семейств биопленок показал наличие генов, соответствующих процессам производства и окисления метана независимо от обилия Н2 (БЬгепк 2017), так что эти биоценозы могут выступать как в роли потребителя Н2, так и быть одновременно потребителями СН4 (анаэробное окисление метана) и его источниками (Вга2еИоп е! а1. 2006). Кроме того, анализ нуклеотидной последовательности гена 16Б гРЫА выявил штаммы, принадлежащие ТНегшососса1в8 и некультивированной Сгепагскаео(а, обнаруженных в жидкостях источников. В самых старых колоннах (с более низкими температурами и меньшим сбросом) в порах преобладала единая группа анаэробных метанотрофных архебактерий (Вга2еИоп е! а1. 2006). Редкие последовательности были обычно более изобильны в молодых колоннах, чем в более старых. Небольшое число видов также преобладает в бактериальных сообществах, в основном это метилотрофические виды (БЬгепк 2017; Вга2еИоп е! а1. 2010). Разнообразие и плотность микробного сообщества зависит от местоположения колонн (их внутренних и внешних стенок), где минералогия и химия жидкостей может сильно отличаться от внутренних частей колонн. Пористые внутренние стенки колонны поддерживают простейшие сообщества биопленки, а архе-бактерии с высокой степенью схожести с серо-и метано-окислителями гораздо более обильны в колоннах с источниками жидкости температурой вплоть до 90°С и высоким рН. Бактерии всегда были более разнообразны, чем простейшие виды во всех исследуемых выборках вне зависимости от места сбора образцов и возраста колонн (Вга2еИоп е! а1. 2010).

Использованные нами в качестве микробного компонента модели Methanosarcinales имеют склонность в ходе деления формировать иррегулярные скопления клеток (Balch et al. 1979) разных размеров, но обычно с толстой стенкой клетки (около 0.5 мкм). В кинетических моделях мы предполагаем, что эти клетки двух диаметров (размеров): в диапазонах 0.05—1.0 мкм и 1.5—2.5 мкм, а концентрация клеток составляет около 105 клеток/мл в обоих случаях.

Наблюдались также временные вариации состава газа струи и потоков вещества от отдельных источников на Энцеладе (Saur et al. 2008). Основная компонента струи — вода в двух состояниях: твердом и парообразном. Предыдущими миссиями были также обнаружены и другие органически и биологически важные составляющие. Кассини во время своего самого близкого прохода около поверхности Энцелада в июле 2005 года собрал данные с помощью масс-спектрометра ионов и нейтральных частиц. Наилучшая оценка состава дает: 91 ± 3% H2O, 3.2 ± 0.6% CO2, 4.0 ± 1.0% N2/CO и 1.6 ± 0.4% CH4. H2 был также обнаружен в паре струи (Waite J. H. et al. 2007). Органические компоненты — это те, которые содержат ковалент-но связанные комбинации углерода с одним или несколькими из пяти других основных элементов жизни: водорода, азота, кислорода, фосфора и серы, но следует отметить, что размер и сложность молекул не являются достаточными условиями для того, чтобы причислить их к формам жизни, так как некоторые очень большие молекулы на основе углерода являются абиотическими (Kahana et al. 2019).

Многократные звездные и солнечные измерения колонковой плотности пара, сделанные по ультрафиолетовым данным и данным масс-спектрометра, предполагают множественные газовые джеты с эмиссией нескольких видов: 1) высокоскоростной газ (т.е. термическое расширение газа сквозь узкие каналы, ускорение под действием давления) и 2) низкоскоростной термический (как сублимация твердого вещества или жидкости на поверхности) (Teolis и Perry, 2017). Согласно этим данным полный выход H2O в струе составляет от 180 кг с-1 до 250 кг с-1. Поданным масс-спектрометра ионов и нейтральных частиц (INMS), поток H2O, наоборот, сильнее варьируется в пространстве: от 200 кг с-1 до 1000 кг с-1 (Hedman et al. 2018). Все эти данные показывают, что и структура самой струи, и структура геологической формации, откуда она берет начало, могут меняться. Обнаружение натриево-солевых зерен льда привело некоторых исследователей к выводу, что струи могут содержать в том числе и вещество океанского дна (Hsu et al. 2015).

В последние годы многие исследователи заинтересовались возможностью существования зон обитания микроорганизмов в среде с этими разнообразными условиями (Bedrossian et al. 2017; Steel et al. 2017; Porco et al. 2017).

Мы можем выделить четыре основных типа теоретически обитаемых областей около южного полюса Энцелада:

• океан в зоне вулканической активности и около нее;

• поверхность океана внутри щелеподобных формирований;

• ледяная поверхность около струйных источников куда более тяжелые частицы попадают с внутренней части подповерхностного океана и, наконец,

• сами гейзеры.

В настоящей работе мы приняли во внимание четыре процесса, связанных с теоретически обитаемыми зонами: 1) формирование суспензии микроорганизмов и свойства клеток около дна океана, 2) процесс транспортировки клеток к поверхности океана, 3) транспортировка клеток в самой струе и 4) процесс выпадения микробных частиц.

Согласно Porco et al. (2017), мы предположили наличие микробных экосистем в ультрамафических, щелочных, гидротермальных источниках на дне океана с концентрацией клеток, аналогичной LCHF, около 105 клеток/мл (Brazelton et al. 2006). Мы применили кинетическое моделирование «частица-в-клетке» для оценки концентрации и параметров клеток на поверхности океана, затем использовали кинетическую симуляцию «частица-в-клетке» для моделирования формы и плотности струи и, наконец, вычислили распределение плотности клеток в струе и в процессе выпадения. Мы применили спектральные данные четырех типов микробных штаммов для оценки возможной разности отражательной способности у разных струй, с микробной компонентой и без нее.

Для выявления спектральных признаков микробов на основе модели с параметрами Methanosarcinales мы использовали для измерений разные бактериальные клетки. Получение спектральных данных для отобранных концентраций микроорганизмов в диапазоне видимого света и в ближнем инфракрасном диапазоне, согласно исследованию Porco et al. (2017), было выполнено с помощью спектрометра Evolution 220 и мультиспектральной камеры QUERCUS 6, разработанной и созданной

в институте авиации в Варшаве (Zalewska et al. 2017). Мы использовали четыре чистых штамма, выделенных из пласта биопленки на факультете биохимии микробов в институте биохимии и биофизики Польской академии наук в Варшаве. Все четыре штамма отличались по форме и цвету. Для спектральных измерений была приготовлена суспензия каждого штамма клеток в дистиллированной воде в концентрации 105 клеток/мл. Все измерения отражательной способности проводились при температуре 21°C с помощью а) спектрометра Evolution 220 в диапазоне 500—1100 нм с частотой дискретизации 1 нм (рис. 3, внизу слева) и б) мультиспектральной камеры QUERCUS 6 с шестью оптическими фильтрами в полосах: 532 нм, 570 нм, 660 нм, 740 нм, 810 нм и 880 нм (рис. 3, внизу справа). Изображения с использованием камеры QUERCUS 6 были получены с черным фоном (коэффициент отражения меньше 5%).

Мы разработали кинетическую модель «частица-в-клетке» шириной 10 км, глубиной 10 км и высотой 15 км с дном океана на z = -10 км, нижним краем ледяного щита на z = 0 км и поверхностью ледяной коры на z = +5 км. В центре ледяной коры находилась щель шириной 2 км и длиной 10 км, представляющая собой один из провалов «Тигровых полос» (см. рис. 1).

Согласно Porco et al. (2017), мы предполагаем поток тепла на дне океана схожим с LCHF: Ф ~ ~ 0.1 Втм-2, содержание соли около 1.0%, концентрацию микробов в гидротермальных каналах Энцелада около 105 клеток/мл. Необходимо отметить, что параметры кинетической модели соответствует нашим знаниям на текущий момент, однако известны они с большой неопределенностью.

Мы предполагаем, что на воду и на частицы микробной компоненты действуют четыре основных силы: гравитация Энцелада, подъемная сила, сопротивление воды и сила градиента давления, описываемые следующими параметрами: масса m, вектор скорости v, температура T, объем V и плотность р. Далее мы заменили непрерывные поля их дискретными аналогами, определенными на сетках клеток (в этой модели использовались 100 х 100 х х 150 кубических клеток с длиной стороны 100 м ).

Мы ввели конечное пространственное разрешение г ^ ¿Аг, разрешение момента р ^ jAp и конечное временное разрешение t ^ kAt, с целыми значениями ¿, j и k. Все дифференциальные уравнения потенциальных полей были заменены простыми разностями.

Мы предполагаем радиус Энцелада rsurface = = 250 км, толщину ледяного щита 5 км и глубину океана 10 км. Если внутренняя плотность

Энцелада постоянна и расстояние между центром Энцелада и дном океана гьо«от = 235 км, ускорение силы тяжести на дне океана на 6% меньше, чем на поверхности: д8иггасе = 0.1108 мс2 и дъо«от = 0.1042 мс-2 соответственно. Гравитационный потенциал был вычислен для каждой сетки согласно стандартному уравнению: и (ж, у, г) = — ^ М

—-, где G — гравитационная постог + гЪойот

янная, М — масса Энцелада, а г — высота над нижней частью ледяного щита. Гравитационный потенциал в нашей модели не зависит от ж и у. Мы предполагаем основную температуру океана постоянной (около 0°С). Подъемная сила зависит от плотности частицы и объема смещенной океанской воды, согласно стандартной формуле В =

— ррайю1е ~^1вр1асеё.

Сопротивление океана было введено согласно уравнению Гагена—Пуазейля АР = 8^Ь^/(пг4) где АР — разность давлений между соседними сетками клеток, ^ — динамическая вязкость воды, Ь — размер клетки, г — радиус частицы, V — объемный расход жидкости, вычисленный перемножением поперечного сечения частицы, ее скорости и временного шага итерации: V = А£пг2 |. Локальное давление для каждой клетки в сетке вычислялось с помощью интеграла:

Р = кт р (г) [1т2] д (г) ^ где д — локальное ускорение тяжести. Так как в нашей модели давление зависит только от силы сопротивления по оси г, оно также параллельно оси г.

Модифицированная симуляция была применена к модельной плотности струй, форме, составу, пространственному многообразию в случае теплового расширения. Для оценки скоростей и температур частиц струи мы использовали концентрацию частиц, массу и температуру поверхности океана. Масса частиц моделировалась рандомно согласно нормальному распределению со средней массой 5 х 10-6 мкг. Это значение вычислялось в предположении радиуса частицы водяного льда в струях 1 мкм, согласно фотометрическому поведению струи, наблюдавшемуся научной подсистемой получения изображений Кассини в 2005 году (т.е. на основе изображений от 27-го ноября) (Рогсо е! а1. 2006; Рогсо е! а1. 2017).

Мы предполагали единую силу, действующую на частицы струи, включая частицы микробной компоненты: гравитацию Энцелада. Здесь мы применили 100 х 100 х 250 клеток размерами 600 м (ж и у) и 400 м (г).

Ice shield (5 km thick) Ocean (10 km thick version) Ocean water particles Biological component particles

Hydrothermal field, ocean bottom

Рис. 1. Кинетическая модель «частица-в-клетке» океана Энцелада глубиной 10 км.

Скорости частиц были вычислены в предположении того, что тепловая энергия частицы океана равна кинетической энергии частицы в струе:

Мы предполагали температуру источника струи около 180 K согласно тепловым измерениям Кас-сини (Dong et al. 2014) и максимальную доступную кинетическую энергию, равную разности между тепловыми энергиями согласно температурам: 1) частиц микробной компоненты на поверхности океана (277.64 K) и 2) источника струи (180 K). Если эту энергию конвертировать в кинетическую, в предположении средней массы частицы 6.60 пг, то средняя скорость будет равна 904.5 мс-1 (85% скорости водяных частиц струи).

Мы рассматриваем четыре случая в нашем кинетическом моделировании:

Случаи 6F и 6S можно использовать как представление частиц Methanosarcinales. 1F и 1S были подобраны в качестве сравнения со случаем клеток гораздо меньшего размера.

Мы предположили, что вещество в струе неоднородно и средняя длина свободного пробега фотона варьируется в разных частях струи, но средняя длина свободного пробега может быть определена локально. Мы также использовали предположение о тонкости модельной атмосферы (оптическая толща т мала и е-т « 1 — т и, следовательно, интенсивность отражения может быть определена из уравнения (КиЫак М. А. 1994):

I» (т) = IV (0) — т1и (0)

в предположении о коэффициенте излучательной способности, равном 0. Спектральная отражательная способность была оценена на следующих высотах: 1 км, 5 км, 10 км, 50 км и 80 км и сравнена с изображениями струй, полученных Кассини в полосе 635 нм (рис. 2).

Для каждой высоты мы вычислили максимальную отражающую способность ЯМАХ и стандартное отклонение положения частиц струи ОХДалее мы предположили, что интенсивность отражающей способности имеет вид

R (х, alt.) = RMAXs\\teDX"^ Используя данные модели струи 6S, мы оценили потенциальную интенсивность отражающей способности струи с частицами Methanosarcinales согласно уравнению

R&(x&alt) = (1 + k (x)) RMAXalte

где к вычисляется просто как доля суммарной массы частиц МеШапозагста1вз и суммарной массы частиц воды в модели на конкретной высоте и х.

Рис. 2. Изображение струйного выброса Энцелада, полученное 27 ноября 2005 года с помощью MAC на борту Кассини в полосе 635 нм и использованное для сравнения с результатами моделирования. Отмечены средние высоты микробной компоненты для четырех кинетических моделей (6S—модель с вариациями начальной скорости Methanosarcinales).

Мультиспектральные данные, полученные со спектрометром и мультиспектральной камерой, представлены на рис. 3.

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

Основные параметры частиц микробной компоненты на дне, в середине и около поверхности океана представлены в таблице 2.

Зависимости между вертикальной скоростью частиц, температурой и массой около океанского дна представлены на рис. 4. Зависимости концентрации клеток от глубины представлены на рис. 5.

Среднее положение микробных и водяных частиц в моделировании струи (включая Ме(Напо8агста1е8—6Ъ) представлено в таблице 3 и на рис. 6.

Отражающая способность струи и параметры, измеренные по изображению Кассини (рис. 2), оцененные с помощью модели струи и дополнительных факторов отражения МеШапо8агста1е8 для высот 1—80 км, представлены в таблице 4.

в возможные метаболические процессы в струях Энцелада, еще не доказывает существования микробной формы жизни. Источниками многих из этих соединений могут быть абиотические процессы, возможно, происходящие на дне океана, одинаковые для Энцелада и Гидротермального поля Потерянного города. Под гидротермальными колоннами ЬСНР существует глубокий магматический резервуар, который нагревает вытекающую воду, обогащенную Н2 и СН4 в результате происходящей на перидотитах серпентинизации. Такая реакция и такой глубоко расположенный источник тепла генерируют относительно низкие температуры (40—90°С) и высокий рН растворов (9—11), что способствует выпадению СаСОз и созданию карбонатных колонн в несколько десятков метров.

Серпентинная реакция — это превращение оливина в серпентин через множественные химические реакции. Оливинный кремний в процессе реакции с гидротермальной водой становится водным кремнием с водородом (уравнение 1) (НоеЫег, 2005). Растворенный кремний формирует серпентин после реакции с магниевым оливином (уравнение 2) (ВппШеу, 1961)

Обнаружение Н2 в струе Энцелада масс-спектрометром ^а^е Л. Н. е! а1. 2017) было интерпретировано как возможный источник энергии микроорганизмов, в частности, метаногенных архебактерий (Рогсо е! а1. 2017). Однако важно отметить, что обнаружение 1) органических соединений и 2) химических соединений, вовлеченных

Несмотря на низкое содержание СО2, абиотический метан может производиться с помощью

# Чего #Bac1 #Bac2 #Bac3 #Bac4

■ f 1 t 1 S 1 I 1 J 1

#Bac4 0.019 0.014 0.027 0.045 0.070 0.070

Рис. 3. Вверху: сравнение отражающей способности бактериальных штаммов в оптических полосах: 660 нм (красный), 720 нм («красный край») и 810 нм (близкий ИК), по данным мультиспектральной камеры QUERCUS 6, созданной в Институте Авиации (проект EC Life 11 ENV/PL/000459). Для наглядности отсчеты DN (0—255) домножены на фактор 28. Внизу: Спектральная отражающая способность четырех штаммов в концентрациях согласно Porco et al. (2017): по измерениям на спектометре Evolution220 (слева) и мультиспектральной камере QUERCUS 6 (справа).

Таблица 1. Начальные параметры частиц микробной компоненты

Параметр Начальное значение Источник / Подтверждение

Масса 6.2 ±2.0 х Ю-12 г НПс1епЬгапс1 е1: а1. (2011)

Физический размер 2.0 ± 0.5 мкм НПс1епЬгапс1 е1: а1. (2011)

Температура 80° С ВгагеНоп е1: а! (2010)

Местонахождение 0 < х, у < 10 км г < 10 м

Скорость 0.1 ±0.05 мг1 Минимальное значение в симуляции

для начальной скорости 0.1 м с-1

водорода, сгенерированного в реакции серпентини-зации, во время дегазации магматических растворов (см. уравнение 3) (Sleep, Meibom, Fridriksson, Coleman, Bird, 2004).

С одной стороны, органические соединения вроде метана могут иметь микробное начало, с другой, они могут быть результатом абиотических процессов, например, серпентинизации. Следует заметить, однако, что эти процессы могут поддерживать возникновение среды, схожей с «Первичным бульоном» (КаИапа е! а1. 2019).

Концентрация клеток в жидкостях каналов LCHF была оценена как 109 клеток/мл (Steel et al. 2017) и, менее оптимистично, как 105 клеток/мл (Jakosky и Shock 1998) (Porco et al. 2017). Согласно этой меньшей концентрации, аналогичной LCHF, мы нашли, что концентрация микроорганизмов у основания ледяной коры может быть оценена в диапазоне от 103 до 104 клеток/мл. Мы также установили, что разнообразие в концентрации клеток наблюдается вплоть до нескольких километров над дном океана (3— 4 км). На этом расстоянии концентрация клеток медленно меняется как примерно 200 клеток/мл на каждый километр. Мы можем принять, что в случае

ОБНАРУЖЕНИЕ МИКРОБНЫХ КОМПОНЕНТОВ В ВЫПАДАЮЩИХ ОСАДКАХ ЭНЦЕЛАДА195 Таблица 2. Усредненные параметры частиц микробной компоненты в различных областях океана Энцелада

Параметр Океан Энцелада

Дно (z < -9.0 км) Середина (—5.5 < z < —4.5 км) Поверхность (z > -1.0 км)

Масса, пг 6.42 ± 1.77 6.57 ± 1.74 6.60 ± 1.74

Температура, °С 31.46 ±25.83 8.96 ±3.47 4.49 ± 1.50

Концентрация клеток, кл./мл 4.37 х 104 9.93 х 103 8.93 х 103

гТ 0.8

Vertical velocity, m s-1

Mass, pg

Рис. 4. Скорость, температура и масса частиц биологической компоненты на дне океана (г < —9 км) (вверху) и около поверхности (г > —1.0 км) (внизу).

моделей с более глубоким океаном концентрация клеток у основания ледяной коры будет схожей (т.е. 3 х 103 для модели океана глубиной 40 км). Эта величина параметра должна быть подтверждена в будущих работах аналогичными кинетическими симуляциями. Глубина океана Энцелада по оценкам составляет не менее 30 км (Thomas et al. 2016) и любой микроорганизм со дна океана должен будет подняться на это расстояние. Потенциальный микроорганизм, используемый в качестве модели частицы микробной компоненты, должен иметь возможность размножаться по пути подъема, характеризуемом различными параметрами океана, такими как давление и температура воды. Используемые в нашей модели метаногены могут выживать в океанской воде от 4°C на поверхности океана вплоть до 80°C около гидротермальных каналов. Резюмируя: кинетические модели типа «частица-в-клетке» дают нам возможность моделировать путь микроорганизмов в океане и должны

использоваться в качестве альтернативы жидкостным гидродинамическим моделям. Наш конечный результат — 104 клеток/мл — подтвердил оценку исследовательской команды Порко для потенциальной концентрации клеток на поверхности (Porco et al. 2017). Остаются открытыми вопросы о свойствах микробной компоненты, в частности: 1) форма, плотность и масса клеток, 2) форма микробной компоненты (одиночные клетки, фрагменты биопленки). Свойства Methanosarcinales, используемые в нашем исследовании, могут быть не единственно приемлемыми.

Результаты, полученные для модели струи, указывают на то, что проектируемые траектории научных пролетов орбитальной миссии Энцелада, согласно разработке концепции (Spencer 2010), с высотами около 50 км, могут быть неподходящими для сбора микробных клеток. Средняя высота Methanosarcinales в моделях струй: 0.50 км для модели с низкой начальной скоростью и 0.74 км для

Рис. 5. Концентрация клеток в океане Энцелада.

Таблица 3. Среднее положение частиц микробной и водной компонент в струе

Случай Микробиальный компонент Водный компонент

(ж), км (у), км (z), км (ж), км (у), км (г), км

модели с максимальной начальной скоростью, — фактически почти в 100 раз ниже, чем проектируемая высота миссии. Предложенная низкая высота в 10 км и орбитальная скорость около 170 мс-1 (Porco et al. 2017) подходит только для клеток очень малых размеров (модели 1S: (г) > ~20 км и 1F (г) ~ 30 км). Наш результат также согласуется с ранее высказанной гипотезой о том, что должно наблюдаться варьирование размеров частиц в струях в зависимости от высоты (Hedman et al. 2009) и единственная возможность наблюдать микробную активность в виде Methanosarcinales может состоять в удаленных сенсорных исследованиях источников струй и их непосредственного окружения, но этот результат должен быть подтвержден будущими симуляциями с более усовершенствованной моделью теплового расширения.

Предполагая потенциальную силу источника света около 1/100 солнечного света на земной орбите, параметры чувствительности мультиспек-тральной камеры, согласующиеся с характеристиками среднеугловой камеры MAC на космическом аппарате Кассини и потенциальную разность отражающих способностей льда/воды и клеток около 30% согласно измерениям в нашем примере, следует ожидать изменений от нескольких (на

высоте около 80 км) до 1000 фотонов в секунду (на высоте около 5 км) на каждые 400 квадратных метров струи. Выявление этих отличий по данным спектрального анализа изображений источников струй также представляет серьезную задачу. Даже для разности отражающей способности клеток в видимом свете (менее 0.05) и ближнем инфракрасном диапазоне (более 0.10), для максимального вклада отражающей способности в диапазоне от 10-7 до 10-8 все стандартные мультиспектральные техники, т.е. показатели разности, используемые для раннего обнаружения цветения воды (КиЫак е!; а1. 2016), могут быть недостаточны.

Многообещающей областью обнаружения являются окрестности источника струи, для которых в наших моделях 6Б и 6Р характерен феномен выпадения микробов. Среднее расстояние выпадающих микробов менее 1 км от источника струи и одновременно среднее расстояние падения частиц льда в диапазоне от 6 до 7 км, что может быть хорошим основанием сфокусироваться на данной области в удаленных сенсорных исследованиях. Из-за этой разницы расстояний выпадения концентрация клеток может быть выше, чем в предыдущих исследованиях (Кешр! е! а1. 2010). Приблизительно 68— 93% частиц струй Энцелада выпадают обратно на поверхность (Рогсо е! а1. 2017) и наблюдения этого

Рис. 6. Частицы воды (черный) и органического вещества (зеленый) в струе (вверху) и в наземных отложениях в результате процесса выпадения микробов (внизу) в случае маленьких (менее 1 пг) и больших (6.6 пг) органических частиц с максимально возможной начальной скоростью (904.5 м с-1, Р) и со стандартным равномерным распределением (0.0-904.5 м с-1, Б).

Таблица 4. Максимальная отражательная способность струи и ее горизонтальное стандартное отклонение, вычисленные по изображениям камеры MAC аппарата Кассини в полосе 635 нм и в рамках кинетической модели без микробной компоненты и с ее включением, на пяти высотах; АД —разность между отражательной способностью воды и микробной компоненты. В двух крайних правых столбцах приведены примеры расчетов для штамма #Bac1 в двух оптических полосах: 660 нм (красный) и 880 нм (ближний ИК).

По данным Модель струи без Дополнительная

Высота, Cassini MAC микробной компоненты микробная компонента

км RMAX DX, RMAX DX, RMAX/AR DX, RMAX ДМАХ

км км км (660 нм, #Вас1) т нм, #Вас1)

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

Концентрация микробной компоненты у источников струй Энцелада зависит от массы и плотности клеток. Модель океана глубиной в 10 км с гидротермальным полем, аналогичным LCHF в Атлантическом океана, содержанием соли около 1 %, концентрацией Methanosarcinales примерно 105 клеток/мл около гидротермальных каналов подтвердила потенциальное существование концентрации микробов около основания ледяной коры в диапазоне от 103 до 104 клеток/мл, предсказанной в работе Porco et al. (2017). Кинетическая модель возникновения клеток в океане позволяет нам исследовать зависимость концентрации от физических параметров клеток, таких как масса и плотность. Используя модель водяной струи «частица-в-клетке», мы нашли сильную зависимость средней высоты клеток Methanosarcinales в струе и их массы. Для массы 6.6 пг у нас получилось, что оптимальная область, где могут находиться микробы, не превышает 1 км по высоте и расположена близко к источникам струй (менее 1 км). Этот результат показал, что методы поиска биологической активности, предложенные для орбитальной станции Энцелада (масс-спектрометр и высота орбиты около 50 км), могут быть недостаточными для этого случая. Кроме того, при потенциальном использовании мультиспектральной камеры для этой цели следует сфокусироваться на области около источника струи, где должна быть более высокая концентрация частиц Methanosarcinales шириной 2 мкм и массой 6.6 пг.

Спектральные измерения четырех образцов штаммов бактерий показали разницу между отражающей способностью клеток в видимом и ближнем инфракрасном диапазонах. Однако она слишком мала (приблизительно 0.05 при концентрации около 107 клеток/мл) для стандартного удаленного метода сенсорной регистрации микробов. По нашим оценкам ожидается всего несколько фотонов отраженного света на каждые 400 м2 на высоте 80 км и несколько тысяч на высоте 5 км (разность отражающей способности приблизительно 10-7 —10-8 согласно параметрам накопления среднеугловой камеры Кассини MAC). Это значение может быть серьезной задачей при анализе мультиспектральных изображений, но отличие среднего расстояния, которое для падающих микробов составляет менее 1 км от источника

струи в сравнении с 6—7 км для выпадающих частиц льда и воды размером в 1 мкм, является весомой причиной, чтобы сфокусироваться в будущих удаленных сенсорных исследованиях на этой конкретной области.

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

Работа выполнена при финансовой поддержке института авиации. Мы благодарим наших коллег из отдела удаленного сенсорного обнаружения за проведенные консультации, которые сильно помогли в нашем исследовании. Мы благодарим профессоров Romana Ratkiewicz и Wojciech Konior за помощь в моделировании по методу «частица-в-клетке».

КОНФЛИКТ ИНТЕРЕСОВ

Авторы заявляют об отсутствии конфликта интересов.

СПИСОК ЛИТЕРАТУРЫ

Rev. 61 (1), 81 (1979).

Astrobiology 17 (9), 913 (2017).

Nat. Academy Sci. 107 (4), 1612(2010).

J. A. Baross, Appl. Environmental Microbiology 72 (9), 6257 (2006).

No. 2042, id. 4124 (2017a).

No. 1964, id. 2577 (2017b).

in the Subsurface Ocean of Enceladus. White paper submitted to the Committee on an Astrobiology Science Strategy for the Search for Life in the Universe (2018).

Microbial Component Detection in Enceladus Snowing Phenomenon

J. Kotlarz, U. Zielenkiewicz, N. E. Zalewska, and K. A. Kubiak

Enceladus is an attractive place to look for signs of life thanks to liquid water and the availability of energy. Recent research has proven that the ejected material of Enceladus south pole consists of water vapor, water ice, carbon dioxide, methane and molecular hydrogen. Possible similarities of physical and chemical conditions between Enceladus ocean bottom and the carbonate mineral matrix of actively venting chimneys of the Lost City Hydrothermal Field give an opportunity to create a mathematical model of microbial ascent process through the ice shell. In this study we present first results of particle-in-cell kinetic simulations of microbial distance through 10 km deep ocean. We have obtained results for microbial component— Methanosarcinalessp. analogue—characterized by 6.6 pg mass and 2.0 ^m diameter distribution in Enceladus plumes. We have assumed 0.1 Wm-2 heating process, 5 km ice shell and cells concentration near ocean bottom 105 cells/ml. We have confirmed assumption of Porco research team about cells concentration near ocean surface about 104 cells/ml and vertical density diversity in plumes. We have found that the optimal altitude for microbial component detection is less than 1.0 km and that in-situ measurements done previously by Cassini mass spectrometer and proposed for Enceladus Orbiter mission 50 km altitude would be ineffective.