1 Test Тест

2 Восстановление высотных профилей оптических характеристик марсианского аэрозоля по лимбовым измерениям спектрометра OMEGA миссии MARS-EXPRESS Майоров Б. С. (ИКИ РАН) Научный руководитель: зав. лаб., к. ф.-м. н. Засова Л. В. (ИКИ РАН). Содокладчики: Васильев А. В. (НИИФ СПбГУ), Bibring J.-P. ( L'Institut d'Astrophysique Spatiale, CNRS-Universite de Paris ) апреля 2008 г. Россия, г. Москва, Институт космических исследований РАН. V конференция молодых учёных «Фундаментальные и прикладные космические исследования» Секция II. Атмосфера, гидросфера и литосфера Земли и планет

3 Описание миссии и прибора OMEGA. Европейская миссия Mars-express (запуск 2 июня 2003 г.; космодром Байконур) Орбитальный космический аппарат: с 26 декабря 2003 г. искусственный спутник Марса. Картирующий спектрометр видимого и инфракрасного спектральных диапазонов для исследования минерального состава поверхности Марса. Российское участие: сканирующее устройство

4 Описание прибора OMEGA. Некоторые характеристики спектрометра: Спектральный диапазон: VNIR канал SWIR канал Спектральный диапазон: 0.37 ÷ 1.05 мкм 0.93 ÷ 2.73 и 2.55 ÷ 5.1 мкм Спектрал. дискретизация: мкм мкм мкм Простран. дискретизация: 0.4 мрад 1.2 мрад (мгновенное FOV) на лимбе (с 300 км): 0.6 км 2 км

5 Лимбовые измерения OMEGA. Примеры лимбовых измерений: - Орбита Первое наблюдение на орбите 0291 Координаты лимба: долгота: 13° E широта: - 44° N Пиксель # 2 сканера Пиксели сканера

6 Модель планеты, приближения. - Сферическая форма твёрдого тела планеты с радиусом R > 0; - Сферически-симметричные свойства атмосферы и поверхности. C – центр планеты Сферическая симметрия:

7 Геометрия орбитальных наблюдений. Базируется на точке наблюдения D (точка детектора или наблюдения). Детектор не может быть внутри твёрдого тела планеты ( ); - единичный вектор внешней нормали в точке D.

8 Геометрия орбитальных наблюдений. - единичный вектор направления на Солнце. (Солнце бесконечно далеко от области наблюдений: нет параллакса).

9 Геометрия орбитальных наблюдений. - единичный вектор направления визирования. 4 независимых параметра полностью определяют геометрию наблюдений.

10 Общепринятая орбитальная геометрия. OMEGA и др. орбитальные приборы используют геометрию (i, e, φ, s, h t ): Преобразования координат:

11 Особенности орбитальной геометрии. Общепринятое описание не чувствительно к: - направлению FOV; - моде наблюдений (поверхность или лимб) без использования h t. Критерий для моды наблюдений (только для тупого угла z v ; IFOV=0): - наблюдение поверхности, - разрыв интенсивности и атмосферного пропускания, - Лимбовые наблюдения

12 Лимбовые измерения OMEGA.

13 - Феноменологический подход (излучение – поток энергии) Спектральная (по, монохроматическая) интенсивность (в ДПСК Cxyz) : Общая теория и приближения - Нет поляризации. - Линейная теория (для процессов ослабления и генерации излучения). - Нет перераспределения энергии по длинам волн. Скалярное уравнение переноса излучения (в дифференциальной форме): n – показатель преломления; - коэффициент экстинкции; - коэффициент излучения. - Единичный вектор направления; t – время. - Стационарность поля излучения: -, нет рефракции. Скалярное стационарное уравнение переноса излучения в инвариантной форме:

14 Другие приближения - Падающее солнечное излучение представляет собой пучок параллельно распространяющихся фотонов. - Элементарный объём среды изотропен для падающего излучения, а индикатриса рассеяния аксиально-симметрична.

15 Радиационный код Радиационный код SCATRD [Васильев, 2006] для вычисления интенсивности рассеянного солнечного излучения в сферической атмосфере. Адаптация к орбитальной моде спектрометрических наблюдений: код SCATRD-OFOS [Майоров и Васильев, 2006] Интернет-сайт: Особенности кода: - оптические параметры атмосферы непрерывные кусочно-линейные функции высоты(неоднородные слои); - таблично-заданные (произвольные) индикатрисы рассеяния; - изотропно отражающая поверхность; - возможность вычисления приближении однократного и двукратного рассеяния по точным формулам (для многократного рассеяния метод Монте-Карло); - тепловое излучение не учитывается; - детальные настройки расчётов; - подробное описание; - работа над совершенствованием кода продолжается.

16 Цель исследования. Надирные измерения: существенный вклад поверхности в спектр (свойства изменяются в пространстве, если нет пылевой бури). Лимбовые измерения: для определения свойств аэрозоля (атмосферы). Получить из этих данных оптические свойства марсианской аэрозольной среды (разработка алгоритма)

17 Анализ лимбовых данных OMEGA. Орбита 291, первое наблюдение. - λ = μm узла по высоте, задающих свойства («аэрозольной») атмосферы: от поверхности (h=0 км) до верхней границы (h = 53 км). - Параметризация аэрозоля: [Ockert-Bell M. E. et a., 1997] Henye-Greenstein индикатриса - - Радиус Марса R = 3395 км ([ Allen, 1973 ]). - Ошибка расчётов методом Монте-Карло 1 %. - FOV не учитывалось. [Mayorov et al., 2006]

18 Анализ лимбовых данных OMEGA. First rough estimation: calculation for exponentially distributed aerosol:

19 Результаты анализа данных OMEGA. Retrieving vertical distribution of aerosol: analogically to "onion peeling" technique.

20 Выбор аэрозольных каналов. Были проанализированы данные спектрометра OMEGA в диапазоне длин волн от 0.4 до 3.4 мкм, в котором были выбраны участки свободные от полос поглощения газов (CO 2, H 2 O, CO и др.) с учётом разрешения прибора. Выбраны аэрозольные каналы для анализа данных спектрометра OMEGA.

21 Список аэрозольных каналов. Выбранные аэрозольные каналы для анализа данных спектрометра OMEGA Условный номерДлина волны, мкмКанал OMEGA VIS VIS VIS VIS VIS VIS IR C IR C IR C IR C IR C IR C

22 Лимбовые измерения OMEGA.

23 Параметризация свойств аэрозоля Аэрозольная оптическая модель может быть записана как функция Аэрозоль: однородные сферические частицы => оптические свойства полностью определяются микрофизическими параметрами (комплексным показателем преломления вещества фракций; концентрацией частиц отдельных мод фракций, их функциями распределения по размерам). I – общее число мод функции распределения аэрозольных частиц по размерам (с учётом всех фракций); n i (h) – счётные концентрации частиц каждой моды, зависящие от высоты h; m i (λ) – комплексный показатель преломления вещества каждой фракции; p i,,k (h) - параметры функции распределения частиц по размерам для каждой моды k= 1,…, K(i) параметров; x( ) – известный алгоритм расчёта оптических параметров ансамбля частиц для единичной концентрации. Мало-параметрическая модель (минимизируем количество параметров): - в составе аэрозолей доминирует одно вещество (m i (λ) - фиксированная функция); - число параметров функции распределения частиц по радиусам = 2; - число мод распределения = 1.

24 Параметризация свойств аэрозоля Параметризация спектральной зависимости характеристик аэрозолей Марса. Параметра(-функции) всего три: профили концентрации частиц n(h) и профили параметров функции распределения p 1 (h), p 2 (h). При этом зависимость от концентрации линейная. Это позволяет перейти непосредственно к параметризации спектрального хода оптических аэрозольных характеристик – функции подлежащие определению параметры модели уже не зависят от длины волны. Задав на определённой высоте концентрацию аэрозолей и параметры функции распределения, находим требуемые оптические характеристики X. Учитывая возможности современных компьютеров, для указанной аппроксимации предлагается предварительная табуляция функции на заданных сетках параметров. Задача параметризации спектрального хода оптических аэрозольных характеристик решена не в традиционном стиле аналитической аппроксимации, а в современном, основанным на возможностях сегодняшней вычислительной техники.

25 Компьютерные коды CompAM, реализующие параметрическую аэрозольную модель. [Майоров и Васильев, 2008]. Интернет-страница: Особенности кода: - сначала генерируется база данных, хранящая таблицу, а затем используется и собственно функция параметризации, выдающая значения оптических характеристик x по конкретно заданным значениям λ, p 1, p 2. - табуляция характеристик аэрозоля с заданной точностью (неравномерная трёхмерная таблица). - извлечение данных для произвольных значений осуществляется трёхмерной линейной интерполяцией по таблице. -реализация параметризации выполнена в достаточно общем виде: предусмотрено задание произвольных аэрозольных веществ и нескольких функций распределения частиц по размерам. -детальные настройки расчётов; -подробное описание; -работа над совершенствованием кода продолжается. Код параметризации свойств аэрозоля

26 Для уменьшения числа параметров мы зафиксировали при параметризации значение [Korablev и др., 2005], т. е. рассматривали функцию распределения с двумя параметрами: r m - модальный радиус (в мкм), α - безразмерный (отвечающий за ширину распределения). Параметризация свойств аэрозоля Марса Для рассматриваемой нами задачи в качестве аэрозольного вещества для атмосферы Марса использовались данные по комплексному показателю преломления из полуэмпирической модели [Ockert-Bell и др., 1997]. Распределение частиц по радиусам – модифицированное гамма-распределение

27 Анализ лимбовых данных OMEGA. Экспоненциальное распределение аэрозоля (n 0 = 2 см -3, H 0 = 10 км, r m = 0.05 мкм, α = 1.5)

28 Результаты анализа данных OMEGA. Экспоненциальное распределение аэрозоля (n 0 = 2 см -3, H 0 = 10 км, r m = 0.05 мкм, α = 1.5) Фиксированное значение параметра ширины распределения частиц по размерам α = 1.5)

29 Выводы. - Радиационный код SCATRD успешно адаптирован к орбитальным спектрометрическим измерениям (код SCATRD-OFOS). - Предложена параметризация спектральной зависимости оптических характеристик аэрозолей с малым числом параметров (код CompAM) (для задач дистанционного зондирования атмосферы Марса). - На основе созданных авторами вышеуказанных кодов разработан алгоритм восстановления высотных профилей микрофизических параметров аэрозоля из орбитальных (лимбовых) спектрометрических измерений (для спектрометра OMEGA миссии Mars-Express).

30 Дальнейшая работа. - Радиационный код SCATRD и его применение к орбитальным спектрометрическим измерениям (код SCATRD-OFOS). - Код параметризации спектральной зависимости оптических характеристик аэрозоля CompAM (прежде всего, для задач дистанционного зондирования атмосферы Марса). -Восстановление высотных профилей микрофизических параметров аэрозоля из орбитальных (лимбовых) спектрометрических измерений (для спектрометра OMEGA миссии Mars-Express) в полностью автоматическом режиме с применением производных и МНК. -Учёт FOV.

31 Спасибо за внимание! Вопросы?