Бесплатный автореферат и диссертация по наукам о земле на тему

Солнечный блик как "инструмент" исследования океана из космоса

ВАК РФ 25.00.28, Океанология

Автореферат диссертации по теме "Солнечный блик как "инструмент" исследования океана из космоса"

Министерство образования и науки Российской Федерации

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ РОССИЙСКИЙ ГОСУДАРСТВЕННЫЙ ГИДРОМЕТЕОРОЛОГИЧЕСКИЙ УНИВЕРСИТЕТ

На правах рукописи УДК 551.465:621.396.969:551.46.09

Мясоедов Александр Германович

СОЛНЕЧНЫЙ БЛИК КАК «ИНСТРУМЕНТ» ИССЛЕДОВАНИЯ ОКЕАНА ИЗ

КОСМОСА

Специальность 25.00.28 - Океанология

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата физико-математических наук

2 7 МАР 2014

рггмы

Санкт-Петербург, 2014

005546412

005546412

Министерство образования и науки Российской Федерации

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ РОССИЙСКИЙ ГОСУДАРСТВЕННЫЙ ГИДРОМЕТЕОРОЛОГИЧЕСКИЙ УНИВЕРСИТЕТ

На правах рукописи удк 551.465:621.396.969:551.46.09

Мясоедов Александр Германович

солнечный блик как «инструмент» исследования океана из космоса

Специальность 25.00.28 - Океанология

автореферат

диссертации на соискание ученой степени кандидата физико-математических наук

сн-сх

рггмы

Санкт-Петербург, 2014

Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования "Российский государственный гидрометеорологический университет" (РГГМУ)

Научный руководитель:

доктор физико-математических Кудрявцев Владимир Николаевич наук

Официальные оппоненты:

доктор физико-математических Ермаков Станислав Александрович

доктор физико-математических Мелентьев Владимир Владимирович

наук, профессор

Ведущая организация:

Институт физики атмосферы им. A.M. Обухова Российской академии наук (ИФА РАН)

Защита состоится «24» Апреля 2014 г. в 15 час. 30 мин. на заседании Диссертационного совета Д-212.197.02 при Российском Государственном Гидрометеорологическом Университете по адресу: 195196, Санкт-Петербург, Малоохтинский пр., 98 Электронная почта: mag@rshu.ru

С диссертацией можно ознакомиться в библиотеке Российского Государственного Гидрометеорологического Университета или на сайте: http://www.rshu.ru/universitv/dissertations/

Автореферат разослан «_»_2014 г.

Ученый секретарь

наук

диссертационного совета, кандидат географических наук

Воробьев Владимир Николаевич

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ Актуальность темы диссертации

Современную научно-исследовательскую деятельность в области океанологии и метеорологии, которые плохо обеспечены контактными данными, уже невозможно представить без использования данных спутниковых наблюдений. Существующие методы обработки и анализа спутниковых измерений обеспечивают получение данных об «огромном» наборе параметров, характеризующих систему «атмосфера — земля — океан», которые в настоящее время широко используются в задачах мониторинга и прогноза окружающей среды. В настоящее время ряд спутниковых геофизических продуктов находится в открытом доступе. В приложении к океанографическим исследованиям, спутниковые данные могут быть получены, например, из центра данных Physical Oceanography Distributed Active Archive Center (PO.DAAC) -американский центр HACA (http://podaac.ipl.nasa.gov/'). Centre ERS d'Archivage et de Traitement - (CERSAT) - французский центр данных института IFREMER (http://cersat.ifrcmer.fr/'): информационный портал спутниковых данных РГГМУ (SATIN - http://satin.rshu.ru/).

В то же время, перспективы развития исследований Земли из Космоса неизбежно требуют создания новых подходов и методов обработки, анализа и использования спутниковой информации. Оптические методы исследования Земли являются наиболее развитыми и широко используемыми в оперативной практике. В настоящее время на орбите Земли находится большое количество сканеров, работающих в оптическом диапазоне (например, сканеры MODIS на спутниках Terra и Aqua, радиометры AVHRR на серии спутников NOAA). Одно из основных применений данных оптических сканеров, - изучение «цвета» Океана (содержание фитопланктона и минеральной взвеси, биогеохимические характеристики), а также температуры его поверхности.

При изучении оптических характеристик Океана, солнечная радиация, отраженная от морской поверхности, является шумом по отношению к радиации, рассеянной в верхнем слое Океана. В областях солнечного блика отражённая радиация составляет значительную часть регистрируемого излучения, что исключает возможность применения алгоритмов восстановления «цвета» Океана. Существование солнечного блика приводит к тому, что огромная часть спутниковых сканерных данных (до 30%) не может быть использована

в классических океанографических приложениях. Области, где восстановление параметров цвета Океана по спутниковым данным невозможно, маскируются для конечного пользователя и, таким образом, «выбрасываются в мусорный ящик». Пример маскировки областей поверхности океана, «засвеченных» солнечным бликом в данных сканера МОБК, приведен на Рис. 1.

Рис. 1. «Глобальная» маска блика от 21 Марта 2004 года для МСЮТБ/Адиа.

Повышение информативности данных измерений спутниковыми сканерами определяет актуальность данного исследования. Основная идея состоит в том, что отраженная солнечная радиация несёт информацию о характеристиках «шероховатости» поверхности Океана. В этом случае данные оптических сканеров могут быть использованы для исследования статистических характеристик ветрового волнения и их вариаций, вызванных различными океаническими процессами. Актуальность данного исследования определяется так же необходимостью разработки нового метода, позволяющего использовать отбрасываемые ранее данные оптических сканеров для исследования проявления различных динамических процессов на поверхности Океана. Предлагаемый подход, совместно с существующими радиолокационными (РЛ) методами наблюдения поверхности Океана, открывает новые возможности для мониторинга океанических явлений из Космоса по их поверхностным проявлениям.

Цели работы

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

Задачи исследования

• разработать метод восстановления пространственных вариаций среднеквадратичного наклона (СКН) морской поверхности по полю яркости солнечного блика;

• применить разработанный метод для анализа данных спутниковых оптических сканеров МОБК и МЕШЯ;

• исследовать поверхностные проявления биологических и нефтяных сликов в солнечном блике и в поле СКН морской поверхности, а также исследовать подобие и отличия аномалий «шероховатости» морской поверхности в сликах, измеряемых оптическими и радиолокационными методами;

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

• исследовать связь аномалий характеристик «шероховатости» морской поверхности с параметрами мезомасштабных течений на основе синергетического анализа оптических и радиолокационных изображений;

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

Научная новизна работы состоит в следующем:

• разработан новый метод восстановления пространственных вариаций СКН морской поверхности по изображениям солнечного блика, регистрируемого спутниковыми оптическими сканерами;

• исследованы проявления нефтяных сликов в солнечном блике и в поле СКН морской поверхности. Показано что контрасты СКН в нефтяных сликах систематически ниже

контрастов СКН в сликах биологического происхождения, эффективный коэффициент упругости тонкой нефтяной плёнки может быть задан как Е= 1 ЪмН/м\

• продемонстрировано, что наблюдения Океана в солнечном блике могут являться эффективным инструментом исследования внутренних волн (ВВ). Поверхностные проявления ВВ видны через модуляцию среднеквадратичного наклона морской поверхности. Усиление среднеквадратичного наклона (СКН) происходит в зонах конвергенции течения ВВ, а его подавление - в зонах дивергенции;

• на основе совместного анализа оптических и радиолокационных изображений установлено, что мезомасштабные течения проявляются на морской поверхности в виде пространственных вариаций СКН и обрушений волн;

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

Практическая и научная значимость

Полученные научные результаты реализованы в виде алгоритмов и элементов программного обеспечения и использованы для обработки данных радиолокаторов с синтезированной апертурой (РСА) и оптических изображений, а также восстановления статистических параметров поверхности Океана.

Предложенные алгоритмы и методики были апробированы и внедрены в Международном центре по окружающей среде и дистанционному зондированию им. Нансена (NIERSC), а также в Лаборатории Спутниковой Океанографии (ЛСО, на англ. SOLab)

РГГМУ, в виде элементов спутникового информационного портала SATIN (от англ. S/i7ellite Data Search and Manage /TVformation Portal), для поиска, получения, отображения, распространения и хранения данных дистанционного зондирования (http://satin.rshu.ru/), а также как элемент разрабатываемой синергетической платформы SYNTool (http://syntool.solab.rshu.ru/) ЛСО РГГМУ.

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

Положения, выносимые на защиту

• разработанный метод диагностики пространственных аномалий «шероховатости» поверхности Океана по спутниковым изображениям солнечного блика позволяет работать с различными оптическими спектрометрами благодаря использованию передаточной функции, которая напрямую зависит от наблюдаемых градиентов яркости солнечного блика, без априорного задания плотности распределения уклонов;

• контрасты СКН в нефтяных сликах систематически ниже контрастов СКН в сликах биологического происхождения;

• для одного и того же слика, сформированного тонкой нефтяной плёнкой, контрасты УЭПР примерно в 1.6 раза сильнее контрастов СКН;

• поверхностные проявления ВВ и мезомасштабных течений отчётливо проявляются в модуляциях уклонов морской поверхности в результате усиления среднеквадратичного наклона (СКН) в зонах

конвергенции течения, и его подавления в зонах дивергенции;

• аномалии характеристик ветрового волнения (СКН, обрушения) связаны с зонами дивергенции течений и пространственно привязаны к областям сильных градиентов завихренности полей квази-геострофических течений.

Апробация работы и публикации

Результаты работы докладывались на различных отечественных и международных конференциях и семинарах, в частности на 7-ой всероссийской открытой конференции «Современные проблемы Дистанционного зондирования Земли из Космоса» (Москва, Россия, Ноябрь 2009); 3rd SeaSAR workshop "Advances in SAR Oceanography from ENVISAT, ERS and ESA third party missions" (Frascati, Italy, 25-29 January 2010); ESA Living Planet Symposium (Bergen, Norway, 28 June -2 July 2010); 39th COSPAR Scientific Assembly (Mysore, India, 14-22 July 2012); 6th International Workshop on Science and Applications of SAR Polarimetry and Polarimetric Interferometry, POLinSAR 2013 (28 January - lFebruary 2013, Frascati (Rome), Italy); Asia-Pacific Conference on Synthetic Aperture Radar "Overcoming the Hardships: Responding to Disasters with SAR", (Tsukuba, Japan, 23-27 September 2013).

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

Автор работы принимал участие на всех этапах исследования от постановки задачи до анализа результатов, разрабатывал компьютерные программы, реализующие предложенные в работе методы и алгоритмы, производил обработку спутниковых данных. Структура и объём диссертации

Диссертационная работа состоит из введения, трёх глав, заключения, библиографии, включающей 78 наименований, из них 73 на иностранных языках. Общий объём работы — 118 машинописных страниц, включая 44 рисунка.

СОДЕРЖАНИЕ РАБОТЫ

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

В первой главе описывается метод восстановления пространственных вариаций СКН морской поверхности по солнечному блику, регистрируемому оптическими сканерами из космоса. Разработанный метод применяется к анализу данных спутниковых оптических спектрометров MODIS и MERIS. Описываются разработанный алгоритм и программное обеспечение для восстановления СКН.

В разделе 1.1 даны общие представления о физике рассеяния и распространения видимого излучения в морской среде.

В разделе 1.2 даётся обзор оптических методов исследования Океана из Космоса. В качестве примеров приводятся характеристики двух спектрометров MODIS и MERIS, данные которых использовались в работе. Известно, что отражённый от морской поверхности солнечный свет составляет основной вклад восходящей радиации и создаёт значительные трудности для разработчиков алгоритмов восстановления цвета Океана. Однако, в солнечном блике содержится ценная информация о статистических характеристиках шероховатости морской поверхности, её среднеквадратичном наклоне (СКН), асимметрии и эксцессе, как было впервые показано в работе Кокса и Манка [Сох, Münk, 1954], а также применительно к спутниковым данным - Бреоном и Хенриотом [Вгёоп, Henriot, 2006].

В разделе 1.3 описывается метод восстановления пространственных вариаций среднеквадратичного наклона (СКН) морской поверхности по солнечному блику, регистрируемому оптическими сканерами из космоса. Вариации СКН связаны с поверхностными проявлениями различных процессов, происходящих в верхнем слое Океана, например, искусственными и биологическими сликами, внутренними волнами, границами мезо-масштабных течений и фронтальными разделами. С этой точки зрения, предложенный метод может рассматриваться как определенный шаг в направлении развития методов диагностики Океана из Космоса.

Приводятся основные соотношения модели изображения

морской поверхности в области солнечного блика и связи вариаций среднеквадратичного наклона (СКН) s2 с вариациями яркости морской поверхности В:

I -7-4. <»

В s2 где Т - это передаточная функция:

Г = 1 + + (2)

2 рдт])

% = ZX/ s и т] = Z / 5, -- "безразмерная" плотность

вероятностей наклонов морской поверхности, нормированная на

среднеквадратичный наклон морской поверхности s2, Zx и Z -

наклоны морской поверхности, удовлетворяющие условиям зеркального отражения солнечного излучения в приемную апертуру прибора (см. соотношения (1-4) в работе [Сох, Münk, 1954]).

На спутниковых изображениях солнечного блика, получаемых различными приборами, в силу технических особенностей этих приборов, поля яркости могут иметь градиенты как в одном, таки в двух (ортогональных) направлениях. При наличии двумерного поля яркости солнечного блика (как, например, в данных сканера MODIS) передаточная функция Т может быть определена по осредненным градиентам яркости, которые берутся непосредственно из изображения солнечного блика. Градиенты безразмерной плотности вероятности наклонов р находятся из градиентов крупномасштабной яркости солнечного блика следующим образом:

т2

Рд% вд% 1+z;+z

7] dp 7] dB 4 Z

У 2

У_

(3)

рдт] В дт] l + Z2x+Z2

У

Градиент яркости в пространстве ( ^' ^ ) можно выразить через градиенты поля в двух ортогональных направлениях ( ' и п ) как

ldB=z У, 1п(Д cos ev) • V„Zy - V„ In В ■ V,Zy

' Д . (4)

Bdrj y A

где лапласиан Д = V,ZX • VnZy - VnZx • V,Z . В применении к данным сканера MODIS У ¡В соответствовало градиенту в направлении сканирования, а У „В - в перпендикулярном направлении

(вдоль полета спутника).

В ряде случаев поля яркости морской поверхности имеют характер одномерного поля (как например в данных сканера MERIS) и для определения Т необходимо задание модели плотности распределения вероятности морских уклонов. Для Гауссова распределения вероятности наклонов передаточная функция принимает следующий вид:

Т а2%2 + а,2¿77-ахг]2, (5) где коэффициенты а,, а2 и ап зависят от направления ветра и анизотропии наклонов поверхности; для изотропного распределения наклонов: а1=а2= \ и а12 —0.

В разделе 1.4 рассматриваются основные технические характеристики приборов MODIS и MERIS, а также особенности формирования спутниковых оптических изображений, полученных этими сканерами. Показаны их отличия и сходства в области солнечного блика. Обосновано, почему в работе используются исходные спутниковые данные (уровня LIB) в красном канале.

В разделе 1.5 описывается процедура и программа восстановления СКН по полям яркости спутниковых оптических изображений MODIS и MERIS. Приводится краткое описание языков программирования, использованных для реализации описанных алгоритмов. Описывается пошаговая процедура применения разработанного метода, включая этапы чтения, обработки и восстановления СКН на примере спутниковых снимков MODIS и MERIS района течения Гольфстрим у берегов Флориды. Дается информация об использовании вспомогательных данных для обработки изображений, включая файлы с геолокацией и

геометрией съёмки, а также вспомогательные данные о приводном ветре и маске Земли.

Во второй главе метод, описанный в первой главе, применяется для исследования морской поверхности, покрытой нефтяными плёнками. Приводится совместный анализ полученных результатов с данными радиолокаторов с синтезированием апертуры (РСА), и раскрываются преимущества синергетического подхода в исследовании поверхностных сликов.

В разделе 2.1 разработанный метод восстановления контрастов СКН применен к анализу поверхностных проявлений разливов нефти природного происхождения (так называемые «грифоны») в Мексиканском заливе. Показано, что вариации яркости блика, связанные с присутствием нефтяных сликов, бывают как положительными, так и отрицательными. Дается описание явления инверсии контрастов в солнечном блике, и приводятся причины его возникновения. Представлена зависимость контрастов СКН от скорости ветра по серии выбранных нефтяных сликов. Делается вывод, что контрасты СКН в нефтяных сликах систематически ниже контрастов СКН в биологических сликах.

В разделе 2.2 проанализирован катастрофический разлив нефти в результате взрыва на нефтяной платформе Дипвотер Хорайзон (англ. Deepwater Horizon) в Мексиканском заливе. Предложенный подход применён к изображениям солнечного блика, полученным приборами MODIS (MODIS/Terra, Май, 24, 2010, 16:45 GMT) и MERIS (MERIS/Envisat, Май 24, 2010, 16:17 GMT). На Рис.2 приведены

7 is2

контрасты СКН 0, восстановленные по изображениям яркости в солнечном блике.

Рис. 2. Аномалии СКН ^ ^ 5 , восстановленные по данным прибора МЕШБ (вверху) и МСЮ18 (внизу). Красными стрелками обозначены области инверсии контрастов, где восстановление контрастов СКН невозможно. Яркие участки контрастов СКН, обведённые жёлтым контуром, не относятся к особенностям шероховатости морской поверхности, а, скорее всего, отражают оптические свойства (цвет) толстой нефтяной плёнки, толщина которой значительно большей длины волны красного света.

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

восстановленные значения аномалий СКН противоречат ожидаемому эффекту подавления коротковолнового волнения в сликах, а контрасты яркости поверхности определяются оптическими свойства нефти (её цветом).

В разделе 2.3 рассмотрен случай синхронной съёмки нефтяного разлива 26 Апреля 2010г. приборами MERIS (15:58 GMT) и ENVISAT ASAR (15:56 GMT), представленный на Рис. 3.

Синхронные РСА и оптические изображения дают возможность проанализировать подобие и различия между проявлениями одних и тех же нефтяных сликов. Показано, что оптические и PJI-контрасты одного и того же слика, сформированного тонкой нефтяной плёнкой, хорошо коррелируют друг с другом. При этом радиолокационные контрасты примерно в 1.6 раза сильнее контрастов СКН. Для интерпретации наблюдаемых контрастов удельной эффективной площади рассеяния (УЭПР ) и СКН в сликах использовалась модель спектров коротких волн, предложенная в работах [Kudryavtsev и др., 2005; Yurovskaya и др., 2013]. Получено, что модельные контрасты согласуются с наблюдениями в том случае, если эффективный коэффициент упругости равен Е=\5мН/м.

Показано, что при совместном анализе контрастов яркости морской поверхности в солнечном блике и РСА контрастов в сликах, можно сделать качественное заключение о толщине пленки. Желтым контуром на Рис. 3 очерчена область слика, в которой толщина нефтяной пленки существенно превышает длину волны красного цвета. В этом случае, контрасты УЭПР отрицательны (в силу подавления коротких волн нефтяной плёнкой), а восстановленные значения контрастов СКН имеют «нереалистичные» положительные значения, обусловленные цветом нефти.

29.5° N

29.0° N

J

28.5° N

Рис. 3. Увеличенные фрагменты изображений ASAR (слева) и MERIS (справа). Изображения содержат нефтяной разлив и представлены в виде контрастов УЭПР (линейные единицы) и контрастов СКН. Толщина нефтяной плёнки в области внутри жёлтого контура значительно больше длины волны

красного света.

В третьей главе рассматриваются примеры исследования суб- и мезомасштабной динамики Океана по оптическим и радиолокационным изображениям. При взаимодействии волн и течений СКН морской поверхности может изменяться, следовательно появляется возможность идентификации течений по изображениям солнечного блика. Используется синергетический подход для исследования поверхностных проявлений мезомасштабных течений по оптическим (включая ИК-канал) и РСА изображениям Океана. Установлено, что аномалии «шероховатости» поверхности Океана, полученные по изображениям солнечного блика хорошо соотносятся с аномалиями на РСА изображениях. Поля аномалий «шероховатости» поверхности океана пространственно коррелируют с зонами дивергенции течений, расположенных в областях сильных градиентов температуры поверхности Океана (ТПО). Проводится анализ и интерпретация данных наблюдений на основе модельных представлений.

В разделе 3.1 разработанный алгоритм применяется к спутниковым оптическим изображениям внутренних волн (ВВ), — как простейшем типе течений. В качестве примера выбран район западно-экваториальной Атлантики, напротив устья реки Амазонки, - область регулярного возникновения интенсивных ВВ, формируемых

полусуточными приливными волнами (см. например, [Ivanov, Pelinovsky, Talipova, 1993]). Изображение MODIS/Aqua этого района, полученное 26 апреля 2009, 16:20, показано на Рис. 4 а. Для иллюстрации работы метода, на этом рисунке приведены: (а) -исходное изображение, (б) - вариации яркости блика, (в) -передаточная функция (см. ур-ние (2)) и (г) восстановленное поле контрастов СКН по уравнению (1). В исходном изображении и в поле СКН отчетливо наблюдаются следующие друг за другом ВВ.

Профиль вариаций СКН вдоль сечения А-В (отмеченного на Рис. 4 в) представлен на Рис. 5 вверху. Проявление ВВ в вариациях СКН имеет «биполярную» форму, с положительными и отрицательными аномалиями СКН. Для интерпретации наблюдений, предположим, что увеличение/уменьшение СКН морской поверхности связано с конвергенцией/дивергенцией поверхностных течений, вызванных ВВ,

т.е. Ks =s jsi ос ди/дх. В этом выражении коэффициент

пропорциональности является «подгоночной» константой, которая определяется сравнением наблюдаемого и модельного СКН. Тогда,

поверхностная скорость определяется как и(х) = сн£ {К3 — (К^сЬс

Определенная таким образом поверхностная скорость течений задавалась в качестве входного параметра для модельных расчетов поверхностных проявлений ВВ по модели формирования PJI-изображения RIM (от англ. Radar Imaging Model) [Kudryavtsev и др., 2005]. Далее, константа си определялась так, чтобы отклонение

модельных расчетов от наблюдений было минимальным. Результаты расчетов представлены на Рис. 5. Так как «профиль» модельных контрастов согласуется с измеренными вариациями СКН, можно заключить, что наблюдаемые модуляции СКН в действительности определяются конвергенцией и дивергенцией течений, индуцируемой ВВ на поверхности. На Рис. 5 так же представлены смещения термоклина, соответствующие «восстановленному» полю скорости, которые согласуются с измерениями ВВ в этом районе [Кудрявцев и др., 2008].

В разделе 3.2 проводится синергетический анализ изображений MODIS и ASAR на примере района течения мыса Игольный), характеризующийся интенсивной мезомасштабной динамикой, отчетливо видимой в изменчивости поля ТПО на Рис. 6 б. Исследования проводятся на фоне сильно меняющегося поля ветра, восстановленного по РСА, Рис. 6 а.

46* № " 46° VI/ 44° VI/

(в) (г)

Рис. 4. (а) Фрагмент исходного изображения МОБВМ^иа (26 апреля 2009, 16:20) района устья реки Амазонки с проявлениями ВВ, представлено в красном канале, (б) Контрасты яркости В / В . (в) Передаточная функция, (г)

2 / 2

Контрасты СКН 5 / £ . Белые области на изображениях - маска облаков. Линия А-В обозначает положение сечения, показанного на Рис. 5.

04

0-3

"й 0 2 у o.i

0

50 ICO 159 200 250 300 350

О 50 100 150 200 250 300 350

к. km

Рис. 5. Профиль контрастов СКН и ЮМ-моделирование сечения А-В, проходящего через ВВ на Рис. 4 в. (вверху) Профиль контрастов СКН (сплошная линия) вдоль сечения А-В. Пунктирная линия показывает результаты модельных расчетов, (внизу) Восстановленные смещения термоклина во внутренних волнах (жирная линия) и скорость течения на поверхности, вызванное

ВВ (тонкая линия).

Рис. 6. Поле ветра и поле поверхностной температуры Океана в районе течения мыса Игольный, юг Африки, (а) Поле ветра, полученное по изображению ASAR WS, с использованием алгоритма CMOD4. Направление ветра взято из модели NCEP. (б) поле ТПО, полученное по данным MODIS. Белые области -маска облаков. Юг Африканского континента выделен коричневым цветом.

На Рис. 7 представлено восстановленное поле контрастов СКН, и контрасты УЭПР РСА. На обоих полях видны проявления фронтального раздела, меандров и вихрей. Типичные значения амплитуды контрастов СКН, наблюдаемые на изображении, имеют величины около 20-30% с масштабами проявления от 1 до 10км. Сопоставляя поля контрастов СКН и УЭПР РСА, можно видеть пространственную «корреляцию» проявления одних и тех же мезомасштабных явлений, что является «неожиданным», поскольку физика формирования оптических и радиолокационных изображений весьма различна.

(а) (б)

Рис. 7. (а) Поле контрастов СКН, полученное из поля вариаций яркости, с применением разработанного алгоритма и поле контрастов УЭПР РСА (б) в районе течения мыса Игольный, юг Африки.

Для выяснения закономерности формирования поверхностных проявлений динамики Океана, необходимо иметь информацию о мезомасштабных особенностях полей течений. Для восстановления мезомасштабной изменчивости поля скорости геострофического течения (ГТ) используется модель, предложенная [Твеш-Ро^апе! и др.,

2008]. В рамках этой модели функция тока ГТ у/ (к, г) и поле ТПО Тз(к) в пространстве Фурье связаны следующим соотношением:

¿(k,Z) = gaJs(;k)exp(n0kZ), (6) fnbk

где n = N / f отношение Прандтля для частот Брента-Вяйсяля N0 и Nb, определяющих, соответственно, мезо- и субмезомасштабные свойства потока. Поле скорости FT определяется через функцию тока у/, как U - (-ikw,ikxi//), или в физическом

пространстве как U = (~ду/ / дх2,ду/ /<Эх,). На Рис. 8 приведено поле завихренности ГТ, восстановленное по полю ТПО, показанному на Рис. 6 б.

Однако, поле ГТ является бездивергентным. Как показано в [Kudryavtsev и др., 2005], влияние такого типа течений на интегральные параметры спектра ветровых волн (СКН, обрушения ветровых волн) является слабым, т.е. проявления ГТ на поверхности океана должны быть «невидимы». Поэтому, наблюдаемые на Рис. 7 поверхностные проявления океанических явлений, скорее всего, связаны с другими факторами.

Взаимодействие ветрового дрейфа с полем скорости ГТ может приводит к генерации достаточно сильной вторичной агеострофической циркуляции [Garrett, Loder, 1981]. Следуя модели [Klein, Hua, 1990] дивергенция агеострофической компоненты течения, генерируемого при взаимодействии Экмановского течения с ГТ, имеет вид:

v-«=-/"'»;rn,( 7)

ОХр

где Q. = dU0 /cbtj —дС/, /дх2 = Ai// - завихренность КГТ, / -параметр Кориолиса, ueß - компоненты Экмановской скорости течения.

Рис. 8. Завихренность поля КГТ и дивергенция поверхностного течения, восстановленные по наблюдаемому полю ТПО. (а) Завихренность поверхностного квазигеострофического течения (КГТ), полученное по полю ТПО, показанному на Рис. 6 б. (б) Поле поверхностной дивергенции, V • и , связанное с вторичной агеострофической циркуляцией, генерируемой в результате взаимодействия Экмановского течения с КГТ. Поле дивергенции инвертировано (показано поле —V ■ и ), поэтому яркие области на изображении соответствуют зонам конвергенции, а тёмные - зонам дивергенции.

Дивергенция поверхностных течений, рассчитанная по полю завихренности (Рис. 8) и Экмановскому дрейфу, соответствующему полю ветра на Рис. 6 а приведено на Рис. 8 б. Видно, что зоны дивергенции/конвергенции привязаны к зонам градиентов завихренности КГТ, и соответственно - к зонам резких перепадов ТПО.

Сопоставляя поля контрастов СКН и УЭПР на Рис. 7 с полями дивергенции поверхностных течений на Рис. 8 б (см. также Рис. 9), можно прийти к выводу, что усиление/ослабление поверхностной шероховатости происходит в зонах конвергенции/дивергенции поверхностных течений. Фрагмент поля дивергенции и контрастов РСА, представленный на Рис. 9, наиболее отчетливо демонстрирует этот факт. Заметим, что основным механизмом формирования поверхностных проявлений течений на РСА изображениях является модуляция обрушений ветровых волн. Поэтому, контрасты РСА изображенные на Рис. 9 а, отражают связь пространственных вариаций обрушений волн с полем дивергенции мезомасштабных течений.

Рис. 9. (а) Фрагмент изображения контрастов УЭПР РСА, полученное по изображению ASAR WS (18 Ноября 2007, 7:24 GMT), и (b) соответствующий фрагмент поля дивергенции поверхностного течения по данным MODIS (18 Ноября, 2007, 12:05 GMT). Яркие области на рисунке (Ь) соответствуют зонам конвергенции, а тёмные - дивергенции течения.

Хорошее соответствие между аномалиями шероховатости поверхности Океана и градиентами ТПО может быть рассмотрено в качестве «экспериментального подтверждения» того факта, что влияние дивергенции поверхностного течения на ветровые волны является основным механизмом формирования поверхностных проявлений мезомасштабных течений.

В разделе 3.3 дается интерпретация данных наблюдений на основе расчета по модели RIM [Johannessen и др., 2005; Kudryavtsev и др., 2005]. Результаты моделирования проявления поля поверхностного течения в полях СКН, обрушения ветровых волн и УЭПР приводятся на Рис. 10. Пространственные распределения этих контрастов подобны полю конвергенции поверхностного течения. Модельные расчеты экспериментально подтверждают обнаруженный факт, что доминирующим параметром, определяющим появление поверхностных проявлений течений, является дивергенция вектора скорости течения.

(В) (г)

Рис. 10 - Фрагменты восстановленного по ТПО поля дивергенции поверхностного течения и модельных значений: контрастов СКН, обрушений волн, и контрастов УЭПР. (а) Дивергенция поля поверхностного течения, полученная по полю ТПО (Рис. 6 б). Яркие области на рисунке (а) соответствуют конвергенции течения, а тёмные - дивергенции. На другие рисунках представлены модельные расчеты проявления вихревого образования в полях контрастов СКН (б), обрушений ветровых волн (в), и контрастов УЭПР (г).

В заключении сформулированы основные результаты работы:

• Предложен алгоритм восстановления пространственных вариаций среднеквадратичного наклона (СКН) морской поверхности по спутниковым изображениям солнечного блика.

• Связь вариаций яркости в солнечном блике с вариациями СКН осуществляется через передаточную функцию, которая может быть определена непосредственно по усредненным 2-мерным градиентам поля яркости солнечного блика (например, для спутникового сканера МСЮК) или на основе априорного задания модели плотности распределения вероятности наклонов морской поверхности (как, например, для спутникового сканера

МЕМБ).

Метод восстановления контрастов среднеквадратичного наклона (СКН) применён к анализу проявления нефтяных сликов естественного и техногенного происхождения по изображениям солнечного блика.

Установлено, что контрасты СКН в нефтяных сликах систематически ниже контрастов СКН в сликах биологического происхождения. Этот результат объясняется различием упругостей нефтяных плёнок и плёнок биологического происхождения. Показано, что эффективный коэффициент упругости для тонкой нефтяной плёнки может быть задан как Е=\ ЪмН/м. Показано, что УЭПР- и СКН-контрасты одного и того же слика, сформированного тонкой нефтяной плёнкой, хорошо коррелируют. При этом контрасты УЭПР поверхности примерно в 1.6 раза сильнее контрастов СКН поверхности в сликах.

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

Предложен синергетический подход для идентификации, восстановления и анализа параметров поверхностных проявлений мезо-масштабных океанических течений по оптическим и радиолокационным изображениям, получаемым из космоса.

В рамках предложенного подхода, поля геострофических течений (ГТ) и вторичных агеострофических течений, с которыми связаны зоны конвергенции и дивергенции течений, могут быть восстановлены по спутниковым полям ТПО полям РСА-ветра.

Установлено, что поверхностные проявления мезомасштабных течений в виде аномалий СКН и обратного рассеяния радиоволн «привязаны» к зонам

конвергенции и дивергенцию поверхностного течения.

• Полученные научные результаты реализованы в виде алгоритмов и элементов программного обеспечения для обработки РСА и оптических изображений и восстановления статистических параметров поверхности Океана. А также как элемент разрабатываемой синергетической платформы SYNTool

(http://syntool.solab.rshu.ru/) Лаборатории спутниковой океанографии (ЛСО) РГГМУ.

Список работ, опубликованных по теме диссертации

Издания из списка ВАК

Мясоедов А.Г., Кудрявцев В.Н.. Оценка контрастов поверхностных проявлений океанических явлений по изображениям солнечного блика. // Учёные записки Российского государственного гидрометеорологического университета No. 16. Научно-теоретический журнал. - СПб.: изд. РГГМУ, 2010. - С. 94-114.

2. Kozlov I.E., Kudryavtsev V.N., Johannessen J.A., Chapron В., Dailidiene I., Myasoedov A.G. ASAR imaging for coastal upwelling in the Baltic Sea. //Advances in Space Research. 2012. № 50. C. 1125-1137.

3. Kudryavtsev V., Myasoedov A., Chapron В., Johannessen J.A., Collard F. Imaging mcsoscale upper ocean dynamics using synthetic aperture radar and optical data. // Journal of Geophysical Research. 2012. №117. С. C04029.

4. Kudryavtsev V., Myasoedov A., Chapron В., Johannessen J.A., Collard F. Joint sun-glitter and radar imagery of surface slicks. // Remote Sensing of Environment. 2012. № 120. C. 123-132.

5. Kudryavtsev V.N., Chapron В., Myasoedov A.G., Collard F., Johannessen J.A. On Dual Co-Polarized SAR Measurements of the Ocean Surface. // IEEE Geoscience and Remote Sensing Letters. 2013. T. 10. № 4. C. 761-765.

6. Зимин A.B., Родионов A.A., Шапрон Б., Романенков Д.А., Здоровеннов Р.Э., Козлов И.Е., Мясоедов А.Г., Шевчук О.И. Работы с научно-исследовательского судна «Эколог» по проекту «Мегагрант» в белом море, выполненные в июле-августе 2012 года. // Учёные записки Российского государственного гидрометеорологического университета. 2013.

Патенты

Мясоедов А.Г., Кудрявцев В.Н. Свидетельство о государственной регистрации программы для ЭВМ №2011610860: GLAMOROS: Оценка контрастов поверхностных проявлений океанических явлений по изображениям солнечного блика. Заявка № 2010617455. Дата поступления 22 ноября 2010г. Зарегистрировано в Реестре программ для ЭВМ 20 января 2011г.

Мясоедов А.Г., Кудрявцев В.Н., Шапрон Б.Ж.А., Коллард Ф.Ж.К. Свидетельство о государственной регистрации программы для ЭВМ №2012660424 DUPOL: «Оценка характеристик поверхности Океана по двух-поляризационным радиолокационным изображениям из Космоса» Заявка № 2012617706. Дата поступления 9 октября 2012г. Зарегистрировано в Реестре программ для ЭВМ 19 ноября 2012г.

Мясоедов А.Г., Козлов И.Е. Свидетельство о государственной регистрации программы для ЭВМ №2012660637 INTERWAVE: «Определение характеристик океанских внутренних волн по их проявлениям в спутниковых радиолокационных снимках морской поверхности» Заявка № 2012618623. Дата поступления 10 октября 2012г. Зарегистрировано в Реестре программ для ЭВМ 26 ноября 2012г.

Зимин A.B., Козлов И.Е., Мясоедов А.Г., Мохнаткин Ф.Ю. Свидетельство о государственной регистрации базы данных №2012621188 «Проявления внутренних волн по данным контактных и спутниковых наблюдений в Белом море в 2010 году (ВВ БМ 2010)» Заявка № 2012621065. Дата поступления 9 октября 2012г. Зарегистрировано в Реестре баз данных 19 ноября 2012г.

Литература других авторов, использованная в автореферате

1. Bréon F.M., HenriotN. Spaceborne observations of ocean glint reflectance and modeling of wave slope distributions // J. Geophys. Res. 2006. Т. 111. № C06005. C. 1-10.

2. Cox C., Münk W. Measurement of the Roughness of the Sea Surface from Photographs of the Sun's Glitter // J. Opt. Soc. Am. 1954. T. 44. № п. C. 838-850.

3. Garrett C.J.R., Loder J.W. Dynamical aspects of shallow sea fronts // Philos. Trans. R. Soc. London Ser. A Math. Phys. Sci. 1981. T. 302. № 1472. C. 563-581.

4. Isern-Fontanet J. и др. Three-dimensional reconstruction of oceanic mesoscale currents from surface information // J. Geophys. Res. 2008. T. 113. № C09005. C. 17.

5. Ivanov V.A., Pelinovsky E.N., Talipova T.G. The Long-Time Prediction of Intense Internal Wave Heights in the Tropical Region of the Atlantic // J. Phys. Oceanogr. 1993. T. 23. № 9. C. 2136-2142.

6. Johannessen J.A. и др. On radar imaging of current features: 2. Mesoscale eddy and current front detection // J. Geophys. Res. 2005. T. 110. №C07017. C. 1-14.

7. Klein P., Hua B.L. The mesoscale variability of the sea surface temperature: An analytical and numerical model // J. Mar. Res. 1990. T. 48. № 4. C. 35.

8. Kudryavtsev V.N. и др. On radar imaging of current features: 1. Model and comparison with observations // J. Geophys. Res. 2005. T. 110. №C07016. C. 1-27.

9. Yurovskaya M. V. и др. Directional short wind wave spectra derived from the sea surface photography // J. Geophys. Res. Ocean. 2013. T. 118. №9. C. 4380^394.

10. Кудрявцев B.H. и др. Оценка контрастов спектра ветровых волн в сликах, вызванных биогенными и нефтяными пленками. Нижний Новгород: ИПФ РАН, 2008. Вып. Препринт №.

Мясоедов Александр Германович

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата физико-математических наук

СОЛНЕЧНЫЙ БЛИК КАК «ИНСТРУМЕНТ» ИССЛЕДОВАНИЯ ОКЕАНА ИЗ КОСМОСА

ЛР№ 020309 от 30.12.96.

Подписано в печать 07.03.14. Формат 60x90 1/16. Гарнитура Times. Печать цифровая. Усл. печ. л. 1,75. Тираж 100 экз. Зак. № 269/1. РГГМУ, 195196, Санкт-Петербург, Малоохтинский пр. 98. Отпечатано в ЦОП РГТМУ

Текст научной работыДиссертация по наукам о земле, кандидата физико-математических наук, Мясоедов, Александр Германович, Санкт-Петербург

Министерство образования и науки Российской Федерации

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ

МЯСОЕДОВ АЛЕКСАНДР ГЕРМАНОВИЧ

СОЛНЕЧНЫЙ БЛИК КАК «ИНСТРУМЕНТ» ИССЛЕДОВАНИЯ

ОКЕАНА ИЗ КОСМОСА

Специальность 25.00.28 — океанология

Диссертация на соискание учёной степени кандидата физико-математических наук

Научный руководитель: доктор физико-математических наук

Кудрявцев В. Н.

РГГМН

Санкт-Петербург, 2014

Содержание

Введение................................................................................4

1 Солнечный блик как "инструмент" исследования Океана из Космоса ... 11

1.1 Общие представления............................................................11

1.2 Оптические исследования Океана из Космоса ................................15

1.2.1 Солнечный блик как источник информации о поверхностных явлениях ......................................................................17

1.3 Модель изображения морской поверхности в области солнечного блика . 19

1.3.1 Основные соотношения..................................................19

1.3.2 Связь вариаций яркости с вариациями СКН..........................21

1.3.3 Восстановления СКН по полям яркости................................23

1.3.4 Пример восстановления СКН модельного блика......................26

1.4 Применение метода к анализу данных МОВ1Б и МЕШБ....................28

1.5 Программа восстановления СКН по полям яркости..........................31

1.5.1 Среда и языки программирования......................................31

1.5.2 Данные М001Б и МЕШБ, используемые в примере..................31

1.5.3 Процедура восстановления СКН по спутниковым оптическим изображениям ............................................................32

1.6 Выводы по главе..................................................................45

2 Оптические и радиолокационные исследования нефтяных сликов.....47

2.1 Нефтяные плёнки природного происхождения ................................49

2.2 Катастрофические нефтяные загрязнения......................................51

2.3 Контрасты СКН и УЭПР нефтяных сликов....................................61

2.4 Выводы по главе..................................................................67

3 Исследование суб- и мезомасштабной динамики Океана по оптическим

и радиолокационным изображениям......................69

3.1 Внутренние волны................................................................70

3.2 Мезомасштабные течения........................................................75

3.2.1 Наблюдения................................................................75

3.2.2 Процедура реконструкции квазигеострофической и агеострофиче-ской циркуляции по ТПО................................................78

3.2.3 Особенности мезомасштабных течений, восстановленные по ТПО,

и их связь с аномалиями РСА сигнала и СКН........................81

3.3 Интерпретация данных наблюдений на основе модельных представлений . 86 3.3.1 Результаты интерпретации данных......................................86

3.4 Выводы по главе..................................................................91

Заключение ......................................93

Список рисунков...................................96

Список таблиц.....................................99

Литература.......................................100

А Основные соотношения модели RIM.......................111

А.1 Приближённое решение задачи о трансформации волн...........114

А.2 Трансформация СКН и обрушений волн....................117

Введение

Современную научно-исследовательскую деятельность в области океанологии и метеорологии, которые плохо обеспечены контактными данными, уже невозможно представить без использования данных спутниковых наблюдений. Существующие методы обработки и анализа спутниковых измерений обеспечивают получение данных об «огромном» наборе параметров, характеризующих систему «атмосфера - земля -океан», которые в настоящее время широко используются в задачах мониторинга и прогноза окружающей среды. В настоящее время ряд спутниковых геофизических продуктов находится в открытом доступе. В приложении к океанографическим исследованиям, спутниковые данные могут быть получены, например, из центра данных Physical Oceanography Distributed Active Archive Center (PO.DAAC) - американский центр HACA (http://podaac.jpl.nasa.gov/), Centre ERS d'Archivage et de Traitement - (CERSAT) - французский центр данных института IFREMER (http: //cersat.ifremer.fr/); информационный портал спутниковых данных РГГМУ (SATIN -http://satin.rshu.ru/).

В то же время, перспективы развития исследований Земли из Космоса неизбежно требуют создания новых подходов и методов обработки, анализа и использования спутниковой информации. Оптические методы исследования Земли являются наиболее развитыми и широко используемыми в оперативной практике. В настоящее время на орбите Земли находится большое количество сканеров, работающих в оптическом диапазоне (например, сканеры MODIS на спутниках Terra и Aqua, радиометры AVHRR на серии спутников NOAA). Одно из основных применений данных оптических сканеров, - изучение «цвета» Океана (содержание фитопланктона и минеральной взвеси, биогеохимические характеристики), а также температуры его поверхности.

При изучение оптических характеристик Океана, солнечная радиация, отраженная от морской поверхности, является шумом по отношению к радиации рассеянной в верхнем слое Океана. В областях солнечного блика отражённая радиация составляет

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

Повышение информативности данных измерений спутниковыми сканерами определяет актуальность данного исследования. Основная идея состоит в том, что отраженная солнечная радиация несёт информацию о характеристиках «шероховатости» поверхности Океана. В этом случае данные оптических сканеров могут быть использованы для исследования статистических характеристик ветрового волнения и их вариаций, вызванных различными океаническими процессами. Актуальность данного исследования определяется так же необходимостью разработки нового метода, позволяющего использовать отбрасываемые ранее данные оптических сканеров для исследования проявления различных динамических процессов на поверхности Океана. Предлагаемый подход, совместно с существующими радиолокационными (РЛ) методами наблюдения поверхности Океана, открывает новые возможности для мониторинга океанических явлений из Космоса по их поверхностным проявлениям.

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

Для достижения поставленной цели в работе решаются следующие задачи:

1. разработать метод восстановления пространственных вариаций среднеквадратичного наклона (СКН) морской поверхности по полю яркости солнечного блика;

2. применить разработанный метод для анализа данных спутниковых оптических сканеров М0015 и МЕШБ;

3. исследовать поверхностные проявления биологических и нефтяных сликов в солнечном блике и в поле СКН морской поверхности, а также исследовать подобие и отличия аномалий «шероховатости» морской поверхности в сликах, измеряемых оптическими и радиолокационными методами;

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

5. исследовать связь аномалий характеристик «шероховатости» морской поверхности с параметрами мезомасштабных течений на основе синергетического анализа оптических и радиолокационных изображений;

6. разработать специализированное программно-математическое обеспечение, сопровождающее разработанные методы.

Научная новизна работы состоит в следующем:

1. разработан новый метод восстановления пространственных вариаций СКН морской поверхности по изображениям солнечного блика, регистрируемого спутниковыми оптическими сканерами;

2. исследованы проявления нефтяных сликов в солнечном блике и в поле СКН морской поверхности. Показано что контрасты СКН в нефтяных сликах систематически ниже контрастов СКН в сликах биологического происхождения, эффективный коэффициент упругости тонкой нефтяной плёнки может быть задан как Е=15мН/м;

3. продемонстрировано, что наблюдения Океана в солнечном блике могут являться эффективным инструментом исследования ВВ. Поверхностные проявления ВВ видны через модуляцию среднеквадратичного наклона морской поверхности. Усиление среднеквадратичного наклона (СКН) происходит в зонах конвергенции течения ВВ, а его подавление - в зонах дивергенции;

4. на основе совместного анализа оптических и радиолокационных изображений установлено, что мезомасштабные течения проявляются на морской поверхности в виде пространственных вариаций СКН и обрушений волн;

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

Практическая и научная значимость работы.

Полученные научные результаты реализованы в виде алгоритмов и элементов программного обеспечения и использованы для обработки данных радиолокаторов с синтезированной апертурой (РСА) и оптических изображений, а также восстановления статистических параметров поверхности Океана. Предложенные алгоритмы и методики были апробированы и внедрены в Международном центре по окружающей среде и дистанционному зондированию им. Нансена (NIERSC), а также в Лаборатории Спутниковой Океанографии (JICO, на англ. SOLab) РГГМУ, в виде элементов спутникового информационного портала SATIN (от англ. SATellite Data Search and Manage INformation Portal), для поиска, получения, отображения, распространения и хранения данных дистанционного зондирования (http://satin.rshu.ru/), а также как элемент разрабатываемой синергетической платформы SYNTool (http: //syntool.solab.rshu.ru/) JICO РГГМУ. В результате применения разработанных методов и алгоритмов, получена возможность использовать данные о яркости поверхности Океана внутри солнечного блика для исследования океанографических явлений по их поверхностным проявлениям, что, в свою очередь, позволило значительно расширить область применимости оптических сканеров. Показано что, применение си-нергетического подхода, основанного на совместном использовании РСА и оптических данных, позволяет лучше понять механизмы проявления океанических явлений на поверхности и выработать предложения по комбинации датчиков и спектральных каналов для повышения эффективности спутникового мониторинга морской среды.

Положениями, выносимыми на защиту, являются:

1. разработанный метод диагностики пространственных аномалий «шероховатости» поверхности Океана по спутниковым изображениям солнечного блика позволяет работать с различными оптическими спектрометрами благодаря использованию передаточной функции, которая напрямую зависит от наблюдаемых градиентов яркости солнечного блика, без априорного задания плотности распределения уклонов;

2. контрасты СКН в нефтяных сликах систематически ниже контрастов СКН в сли-ках биологического происхождения;

3. для одного и того же слика, сформированного тонкой нефтяной плёнкой, контрасты УЭПР примерно в 1.6 раза сильнее контрастов СКН;

4. поверхностные проявления ВВ и мезомасштабных течений отчётливо проявляются в модуляциях уклонов морской поверхности в результате усиления среднеквадратичного наклона (СКН) в зонах конвергенции течения, и его подавления в зонах дивергенции;

5. аномалии характеристик ветрового волнения (СКН, обрушения) связаны с зонами дивергенции течений и пространственно привязаны к областям сильных градиентов завихренности полей квази-геострофических течений.

Апробация работы и публикации.

Результаты работы докладывались на различных отечественных и международных конференциях и семинарах, в частности на 7-ой всероссийской открытой конференции «Современные проблемы Дистанционного зондирования Земли из Космоса» (Москва, Россия, Ноябрь 2009); 3rd SeaSAR workshop "Advances in SAR Oceanography from ENVISAT, ERS and ESA third party missions" (Frascati, Italy, 25-29 January 2010); ESA Living Planet Symposium (Bergen, Norway, 28 June - 2 July 2010); 39th COSPAR Scientific Assembly (Mysore, India, 14-22 July 2012); 6th International Workshop on Science and Applications of SAR Polarimetry and Polarimetric Interferometry, POLinSAR 2013 (28 January - 1 February 2013, Frascati (Rome), Italy); Asia-Pacific Conference on Synthetic Aperture Radar "Overcoming the Hardships: Responding to Disasters with SAR (Tsukuba, Japan, 23-27 September 2013).

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

Личный вклад автора.

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

Структура и объём диссертации.

Диссертационная работа состоит из введения, трёх глав, заключения, библиографии, включающей 78 наименований, из них 73 на иностранных языках. Общий объём работы - 118 машинописных страниц, включая 44 рисунка.

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

В первой главе описывается метод восстановления пространственных вариаций среднеквадратичного наклона (СКН) морской поверхности по солнечному блику, регистрируемому оптическими сканерами из космоса. Разработанный метод применяется к анализу данных спутниковых оптических спектрометров М001Б и МЕШБ. Описываются разработанный алгоритм и программное обеспечение для восстановления СКН.

В разделе 1.1 описаны общие представления о физике рассеяния и распространения видимого излучения в морской среде.

В разделе 1.2 даётся обзор оптических методов исследования Океана из Космоса. В качестве примеров приводятся характеристики двух спектрометров МОВ15 и МЕШБ, данные которых использовались в работе.

В разделе 1.3 описывается метод восстановления пространственных вариаций среднеквадратичного наклона (СКН) морской поверхности по солнечному блику, регистрируемому оптическими сканерами из космоса.

В разделе 1.4 рассматриваются основные технические характеристики приборов М001Б и МЕШБ, а также особенности формирования спутниковых оптических изображений, полученных этими сканерами.

В разделе 1.5 описывается процедура и программа восстановления СКН по полям яркости спутниковых оптических изображений М0015 и МЕШБ. Приводится краткое описание языков программирования, использованных для реализации описанных алгоритмов.

Во второй главе метод, описанный в первой главе, применяется для исследования морской поверхности, покрытой нефтяными плёнками. Приводится совместный анализ полученных результатов с данными радиолокаторов с синтезированием апертуры (РСА), и раскрываются преимущества синергетического подхода в исследовании поверхностных сликов.

В разделе 2.1 разработанный метод восстановления контрастов СКН применен к анализу поверхностных проявлений разливов нефти природного происхождения (так называемые «грифоны») в Мексиканском заливе.

В разделах 2.2 и 2.3 проведён совместный анализ РЛ и оптических изображений катастрофического разлива нефти в результате взрыва на нефтяной платформе Дипво-тер Хорайзон (англ. Deepwater Horizon) в Мексиканском заливе 26 Апреля 2010г.

В третьей главе рассматриваются примеры исследования суб- и мезомасштабной динамики Океана по оптическим и радиолокационным изображениям. При взаимодействии волн и течений СКН морской поверхности может изменяться, следовательно появляется возможность идентификации течений по изображениям солнечного блика. Используется синергетический подход для исследования поверхностных проявлений мезомасштабных течений по оптическим (включая ИК-канал) и РСА изображениям Океана. Установлено, что аномалии «шероховатости» поверхности Океана, полученные по изображениям солнечного блика хорошо соотносятся с аномалиями на РСА изображениях. Поля аномалий «шероховатости» поверхности океана пространственно коррелируют с зонами дивергенции течений, расположенных в областях сильных градиентов температуры поверхности Океана (ТПО). Проводится анализ и интерпретация данных наблюдений на основе модельных представлений.

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

В разделе 3.2 проводится синергетический анализ изображений MODIS и ASAR на примере района течения мыса Игольный), характеризующийся интенсивной