CC BY
38
ЭКСПЕРИМЕНТАЛЬНЫЕ И ЭКСПЕДИЦИОННЫЕ ИССЛЕДОВАНИЯ
УДК 504.064.3:551.466.1
Влияние обрушений ветровых волн на обратное рассеяние миллиметровых радиоволн морской поверхностью
© 2015 Ю.Ю. Юровский1, И.А. Сергиевская2, С.А. Ермаков2,3,4, Б. Шапрон5, И.А. Капустин2,3,4, О.В. Шомина2
1 Морской гидрофизический институт РАН, Севастополь, Россия
E-mail: yyyurovsky@gmail. com
2 Институт прикладной физики РАН, Нижний Новгород, Россия
3 Нижегородский государственный университет им. Н.И. Лобачевского, Нижний Новгород, Россия
4 Волжский государственный университет водного транспорта, Нижний Новгород, Россия 5 Российский государственный гидрометеорологический университет, Санкт-Петербург, Россия
Поступила в редакцию 21.01.2015 г.
Представлены результаты натурных исследований обратного радиолокационного рассеяния на морской поверхности на длине радиоволны 8 мм. Радиолокационные наблюдения сопровождались синхронной видеорегистрацией облучаемой поверхности с целью определения параметров обрушений ветровых волн, влияющих на интенсивность обратного рассеяния, а также для измерения скорости поверхностного течения, вызывающего доплеровский сдвиг частоты сигнала. Показано, что вклад обрушений в удельную эффективную площадь рассеяния зависит от скорости ветра и поляризации излучения и не превышает 20 - 30%. Средние значения мгновенного доплеровского сдвига, не связанного с поверхностным течением, в 1,5 - 2 раза выше ожидаемого доплеровского сдвига от свободных брегговских волн. Анализ видеоданных показал, что полученное несоответствие не может быть объяснено влиянием только видимых обрушений. Возможными причинами наблюдаемого эффекта могут быть пространственная модуляция внутри пятна облучения, влияние связанных паразитных капиллярных волн, микрообрушения, не вызывающие видимых проявлений на морской поверхности.
Ключевые слова: обратное рассеяние, миллиметровые волны, доплеровский сдвиг, поверхностное течение, обрушения ветровых волн, паразитные волны.
Введение
Рассеяние микроволнового излучения на морской поверхности, как правило, описывается двумя физическими механизмами (см., например, обзор [1]): зеркальными отражениями при малых углах падения и резонансным (брегговским) рассеянием при умеренных углах падения (угол падения - угол между вертикалью и волновым вектором падающей радиоволны). Класси -ческой теорией, объединяющей оба механизма, стала двухмасштабная модель, в рамках которой морская поверхность рассматривается в виде резонансных волн, расположенных на более длинных волнах [2]. Такая модель адекватно воспроизводит экспериментальные данные, получаемые в различ-
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 4 2015
37
ных радиодиапазонах и условиях измерений [1, 3, 4]. Вместе с тем наблюдается ряд эффектов, появление которых невозможно описать обозначенными выше механизмами. В их числе можно отметить случаи превышения сигнала обратного рассеяния горизонтальной поляризации излучения и приема (ГГ-поляризация) над сигналом вертикальной поляризации (ВВ-поляризация) [5], а также различия в доплеровских сдвигах частоты принимаемого сигнала в зависимости от поляризации излучения. В сантиметровом диапазоне длин радиоволн указанные эффекты хорошо изучены и обычно ассоциируются с рассеянием на обрушивающихся ветровых волнах [5 - 7].
С другой стороны, развитие микроволновых технологий сделало возможным изготовление радиолокационных станций (РЛС) миллиметрового диапазона и их использование в задачах океанологии. К преимуществам такого рода систем можно отнести сравнительно низкую себестоимость, малые габариты и высокую разрешающую способность по сравнению с более длинноволновыми аналогами. Однако работы, посвященные исследованию особенностей рассеяния радиоволн миллиметрового диапазона, как правило, ограничиваются измерениями при малых углах падения и изучением только интенсивностей рассеяния [8 - 10], что связано прежде всего с традиционным использованием миллиметровых волн в альтиметрических системах. Построение же универсальной модели обратного рассеяния, необходимой для анализа как скаттерометрических данных, так и «.^-изображений, требует дальнейшего более детального экспериментального изучения свойств обратного рассеяния миллиметровых волн.
Цель данной работы - представить некоторые результаты натурных экспериментов по изучению обратного рассеяния на длине радиоволны 8 мм. Радиолокационные наблюдения сопровождались синхронной видеосъемкой морской поверхности с целью определения роли обрушений ветровых волн в обратном рассеянии. Особое внимание уделено исследованию доплеровских сдвигов частоты принимаемого сигнала, поскольку в используемом диапазоне брегговская длина волны лежит в капиллярной области спектра, для которой характерно существование миллиметровой паразитной ряби. Благодаря модуляции промежуточными дециметровыми волнами паразитные волны могут иметь фазовую скорость, отличную от фазовой скорости свободных ветровых волн, вследствие чего наблюдаемые доплеровские сдвиги могут быть завышены [11]. Для изучения этого эффекта в работе анализируются мгновенные доплеровские сдвиги, связанные только с кинематическими характеристиками рассеивателей и не зависящие от их вклада в удельную эффективную площадь рассеяния (УЭПР) морской поверхности.
Оборудование и эксперимент
Измерения проводились на стационарной океанографической платформе Морского гидрофизического института, расположенной в Голубом заливе на южной оконечности Крымского п-ова (пгт Кацивели). Платформа зафиксирована на морском дне на расстоянии ~500 м от берега. Глубина моря вблизи платформы составляет 25 - 30 м.
38
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 4 2015
Радиолокационные записи сделаны с помощью макета доплеровской РЛС непрерывного излучения, работающей на частоте 37,5 ГГц. Станция оснащена двумя идентичными рупорными антеннами для передачи и приема сигнала. Излучаемая волна имеет наклонную поляризацию, в приемном тракте принимаемая волна разделяется на вертикальную и горизонтальную компоненты, тем самым формируя сигналы горизонтальной и вертикальной поляризаций.
Для видеорегистрации облучаемого участка морской поверхности использована цифровая видеокамера. Частота записи кадров составляла 25 Гц при размере изображения 1440 х 1080 пикселей. Углы обзора в горизонтальной и вертикальной плоскостях были равны 47 и 28° соответственно. Видеокамера жестко фиксировалась на корпусе РЛС так, что положение пятна облучения в кадре было неизменным.
Скорость и направление ветра определялись с помощью стандартного чашечного анемометра, установленного на высоте 21 м.
Экспериментальный материал, используемый в данной работе, был получен в ходе нескольких экспедиций в 2011 - 2013 гг. Угол падения изменялся в пределах 0 - 70°, азимут относительно направления ветра менялся от 0° («на ветер») до 180° («под ветер»).
Допустим, что морская поверхность облучается РЛС с заданной двунаправленной диаграммой направленности О (рис. 1). Рассеивающая поверхность может быть представлена в виде совокупности независимых рассеива-телей, имеющих различные локальные скорости и УЭПР а0. Элементарные рассеиватели переносятся поверхностным течением (со скоростью Уйг) и орбитальными движениями волн, длина которых больше, чем характерный масштаб рассеивателя. Согласно [12], такие волны можно разделить на два типа: длинные и промежуточные. Длинные волны имеют длину, намного превышающую размер элемента разрешения РЛС, и их характеристики можно считать постоянными в пределах пятна облучения. Поэтому влияние таких волн на сигнал проявляется только в его временных вариациях. Масштаб промежуточных волн больше, чем размер рассеивателя, но меньше, чем размер пятна облучения. Поэтому они вызывают пространственную модуляцию локальных характеристик рассеивателей внутри пятна облучения и влияют на мгновенные значения сигнала.
Мгновенную и среднюю УЭПР морской поверхности можно представить в виде
где х - пространственная координата; ? - время; Т - время осреднения; - площадь пятна облучения.
Анализ данных
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 4 2015
39
Облучаемая область Рис. 1. Иллюстрация механизма формирования обратно рассеянного сигнала
Мгновенная доплеровская скорость (/) представляет собой лучевую скорость элементарных рассеивателей, осредненную по пятну облучения с весом, пропорциональным локальной УЭПР рассеивателя а0 (х, t) :
V) = ^гяпв +
I х, t )о-0( х, t )0( х)&х х, t )G( x)dx
где в - угол падения.
Особого внимания заслуживает определение среднего значения до-плеровской скорости. Можно рассматривать две величины: среднее от
мгновенной доплеровской скорости у0 и средневзвешенную с учетом мгновенной УЭПР доплеровскую скорость у0 :
— 1
IV ,
I VD(t )o■o(t ^ )dt
В оба параметра вносит вклад корреляция между скоростью элементарных рассеивателей х.С(х, ^ и их локальной УЭПР сг0(х, t) , которая обусловливается геометрической, гидродинамической и аэродинамической модуляцией и описывается в рамках модуляционных передаточных функций [13, 14]. Важное отличие между ^ и vD состоит в том, что последняя включает в себя вклад от волн всех масштабов, в то время как первая - только от промежуточных волн. Средневзвешенная оценка доплеровской скорости vD , по сути, является величиной, измеряемой реальными аэрокосмическими радиолокационными инструментами, элемент разрешения которых намного больше масштаба самых длинных ветровых волн, присутствующих на морской по-
40 МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 4 2015
^ =
верхности. Невзвешенная средняя доплеровская скорость vD может быть измерена специализированными системами, имеющими высокое пространственное разрешение (порядка масштаба рассеивателя) и расположенными, как правило, вблизи изучаемой поверхности. В этом случае доплеровская скорость отражает только среднюю скорость рассеивателей, независимо от их вклада в УЭПР, что облегчает и делает более наглядным анализ данных ввиду отсутствия необходимости рассматривать связь между скоростью рассеи-вателя и его УЭПР. Поэтому в данной работе мы ограничимся только анализом средней мгновенной доплеровской скорости vD .
Обработка радиолокационных данных.
Мощность записанного сигнала преобразовывалась в УЭПР морской поверхности путем калибровки РЛС с использованием специально изготовленных мишеней (металлических уголковых отражателей и сфер). Измерение двунаправленной диаграммы направленности G проводилось путем наблюдения сферической мишени на фоне неба под разными азимутами и углами склонения относительно геометрической оси РЛС, которые определялись с помощью видеокамеры.
Записанный сигнал каждой из поляризаций РЛС представляет собой синфазную и квадратурную компоненты, являющиеся действительной и мнимой частью комплексной амплитуды сигнала A(t). Мгновенный допле-
ровский спектр сигнала S(®D, t) есть функция времени и доплеровской частоты, запишем его в виде
S (oD, t ) =
j t+T/2
^ JA(t)exp(-/®Dt)dt
t-T/2
где ®D - круговая доплеровская частота. Положительные значения ®D соответствуют приближающимся к наблюдателю целям, отрицательные - удаляющимся. Круговая доплеровская частота связана с лучевой скоростью цели: vD = coD /2kr, где kr - волновое число радиоволны. При условии, что средняя вертикальная составляющая скорости рассеивателей равна нулю, лучевая скорость может быть использована для оценки средней горизонтальной скорости движения мишени VD = coD / 2kr sin в . Для удобства анализа мгновенные спектры S(aD, t) были преобразованы в S(VD, t), поэтому в дальнейшем вместо доплеровской частоты coD будем рассматривать горизонтальную до-плеровскую скорость VD.
Имеющиеся записи были разбиты на пятиминутные фрагменты, для каждого из которых рассчитывался мгновенный спектр S(VD, t). Мгновенные значения УЭПР a0(t) и доплеровской скорости VD(t) определялись путем расчета нулевого и первого, нормированного на нулевой, моментов мгновенного спектра S(VD,t) .
2
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 4 2015
41
Обработка видеоданных.
Оценка скорости поверхностного течения. Для определения собственной скорости рассеивателя из измеренного доплеровского сдвига должен быть вычтен сдвиг, обусловленный скоростью поверхностного течения, которая может быть принята постоянной как в пятне облучения, так и в пределах рассматриваемого пятиминутного фрагмента. Вектор скорости поверхностного течения определялся по имеющимся синхронным видеозаписям пятна облучения РЛС.
Идея метода заключается в отслеживании отдельных пузырей и пены, остающихся после обрушения ветровых волн. Сами обрушения имеют время жизни порядка долей секунды и, очевидно, некоторую собственную скорость, в то время как остаточная пена существует на поверхности в течение нескольких секунд и может рассматриваться в качестве пассивного трассера (лагранжева дрифтера). Как показала предварительная обработка, фильтрация событий с временем жизни более 4 с эффективно отделяет пену от коротко-живущих барашков и солнечных бликов.
Исходные видеозаписи были разбиты на одноминутные блоки (1500 кадров). Для «выравнивания» поля яркости из каждого кадра вычиталась матрица средней по блоку яркости изображения. Затем изображения бинаризирова-лись по порогу, величина которого определялась эмпирически так, чтобы все возможные области, которые нельзя отнести к невозмущенной морской поверхности, были обнаружены. Выделенные области объединялись в связанные трехмерные группы пикселей (две пространственные координаты и время). События, время жизни которых было меньше 4 с, отбрасывались, тем самым в рассмотрении оставались только свободно дрейфующие пузыри и фрагменты пены. Пиксельные координаты выделенных областей трансформировались в координаты на невозмущенной морской поверхности. Вектор скорости поверхностного течения, соответствующий данному блоку кадров, вычислялся как средний вектор скорости всех обнаруженных событий.
Выделение обрушений ветровых волн. Видеосъемка морской поверхности, ведущаяся синхронно с радиолокационными наблюдениями, позволяет определить роль обрушений ветровых волн в формировании обратно рассеянного сигнала. Как было показано выше, бинаризация видеоизображений позволяет обнаруживать различные особенности морской поверхности, а именно: обрушения ветровых волн, остаточную пену и солнечные блики. Однако, как показывает анализ видеозаписей и непосредственное визуальное наблюдение моря, блики и не образующие пену обрушения самого малого масштаба практически неразличимы. При этом, очевидно, блики в отличие от обрушений не могут заметно влиять на сигнал микроволновой РЛС. Для того чтобы избежать необходимости разделять эти два вида особенностей морской поверхности, при совместном анализе радиолокационных данных и видео-
42
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 4 2015
изображений рассматривались только кадры, на которых не наблюдалось каких-либо заметных проявлений обрушения гребней волн. При такой выборке данных можно считать, что рассеяние происходит на поверхности без обрушений. Тогда отличия параметров сигнала для этой выборки от тех же параметров для всей записи будут обусловлены обрушениями.
Исходные видеозаписи вплоть до бинаризации обрабатывались по алгоритму для выделения пассивных трассеров, описанному выше, однако фильтрация по времени жизни и размеру выделенной области не проводилась. После этого вычислялась мгновенная доля поверхности Q(t), покрытой обнаруженными областями:
Ц0(х, у)I(х, у, t)<Ыу
Q(t) = -:
|| 0(х, y)dxdy
где I (х, у, {) - бинарное изображение морской поверхности.
Средние значения Q сравнивались с известной эмпирической зависимостью из работы [15] для выявления неудачных записей, которые могли попасть в выборку вследствие неподходящего освещения морской поверхности или облачности.
Анализ видеоданных показал, что определенное малое количество выделяемых по яркости событий практически всегда присутствует на изображении, что мешает получению хорошей статистики для «поверхности без обрушений». Поэтому для получения достаточной выборки для Q(t) был выбран порог, составляющий 0,5%. По сути, эта процедура добавляет в выборку «поверхность без обрушений» некоторое количество событий, масштаб которых легко оценить. Площадь пятна облучения £ей- и его видеоизображения в зависимости от угла падения лежит в диапазоне 5 - 10 м2. Пороговое значение 0,5% для Q(t) оставляет в выборке области, площадь которых меньше, чем £ей- х 0,005 = 0,025 - 0,05 м2. Известно, что обрушения обладают свойством автомодельности и имеют универсальную форму эллипса с эксцентриситетом ~0,9 [16]. Таким образом, минимальная длина обнаруживаемых обрушений лежит в пределах 0,25 - 0,4 м.
Пример совместной обработки видеоизображений и радиолокационных данных приведен на рис. 2. Видно, что всплески в Q(t) случаются на определенной фазе длинной волны, что хорошо демонстрирует известный эффект модуляции обрушений длинными волнами [17]. Поэтому исключение из выборки моментов обрушений волн не даст возможность рассматривать определенную фазу длинной волны, на которой происходят обрушения. Это, в свою очередь, вызовет уменьшение или увеличение УЭПР и доплеровских сдвигов, рассчитанных по такой выборке. Для того чтобы этого избежать,
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 4 2015
43
вычисление средних по выборке параметров было выполнено по специальной методике.
Рис. 2. Пример временной реализации: фазы длинной волны - а; мгновенной доплеровской скорости - б; УЭПР - в; доли поверхности, покрытой обрушениями, - г. Угол падения 55°, азимут - вслед волне, скорость ветра 13 мс"1
Сдвиг доплеровской частоты, вариации которого отражают орбитальные движения длинных волн, был использован для оценки текущей фазы длинной волны Ф . Зависимости между параметрами сигнала и фазой Ф были осред-нены по 10-градусным интервалам, что в результате дало диаграммы распределения параметров сигнала по профилю длинной волны. Пример такого распределения показан на рис. 3.
Рис. 3. Пример распределения мгновенных УЭПР и доплеровской скорости по профилю длинной волны на ГГ-поляризации: а, в - вся запись (красные линии); б, г - «поверхность без обрушений» (синие линии). Горизонтальными линиями обозначены средние значения по всем фазам. Условия наблюдения те же, что на рис. 2
44 МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 4 2015
Описанная процедура осреднения была проведена для всех записей и выборок из них, соответствующих «поверхности без обрушений». В последнем случае среднее по диаграмме отличается от простого арифметического среднего, рассчитанного по выборке, поскольку фазы длинной волны, ассоциированные с обрушениями, вместо исключения из рассмотрения заменяются случаями, когда на данной фазе обрушений не было.
Результаты Вклад обрушений в УЭПР.
Сопоставление средних значений параметров сигнала со значениями, рассчитанными по выборкам «поверхность без обрушений», позволяет судить о вкладе обрушений в сигнал РЛС. Запишем УЭПР морской поверхности в виде
а0 =CTreg(l" Q) + awb Q ,
где а - УЭПР морской поверхности без обрушений; awb - УЭПР возмущенной обрушением поверхности [7].
Поскольку Q является малой величиной и 1 — Q«1, разность
Act = ст0 — а можно рассматривать как компоненту УЭПР, связанную с
рассеянием на обрушениях.
Относительный вклад обрушений в суммарную УЭПР, рассчитанный как Аа / а0, показан на рис. 4 как функция скорости ветра f/21. Данные разделены на три группы, соответствующие малым, умеренным и большим углам падения. Видно, что для скоростей ветра менее 10 мс"1 вклад обрушений в сигнал колеблется в пределах -5 ... +5% на ВВ-поляризации и ~5 - 10% - на ГГ-поляризации. С увеличением скорости ветра вклад возрастает и при [/21 = 15 мс-1 может составлять 20 - 30% на ГГ-поляризации, что приблизительно в два раза больше, чем на ВВ-поляризации.
Рис. 4. Вклад обрушений ветровых волн в УЭПР на ВВ-поляризации (а) и ГГ-поляризации (б) в зависимости от скорости ветра
МОРСКОЙ ГИДРОФИЗИЧЕСКИИ ЖУРНАЛ № 4 2015
45
Доплеровские скорости морской поверхности без обрушений.
Доплеровская скорость рассеивателя зависит от механизма рассеяния, который в первую очередь определяется углом падения в . При малых в рас-сеивателем является участок поверхности, ориентированный нормально к РЛС (зеркальная точка) и движущийся с фазовой скоростью несущей его волны. При умеренных в доминирующим считается резонансный (бреггов-ский) механизм рассеяния. Элементарным рассеивателем в этом случае выступает пакет, состоящий из нескольких периодов брегговских волн, волновое число которых кы удовлетворяет условию кы = 2кг sin в . Скорость такого рассеивателя равна фазовой скорости брегговских волн сьг = / кЬг + укЬт )1/2, где g - ускорение свободного падения, у - коэффициент поверхностного натяжения. Кроме распространения с собственной скоростью, рассеиватели переносятся орбитальными движениями длинных волн, нелинейность которых в среднем приводит к волновому (стоксову) дрейфу поверхности, входящему в измеряемую нами скорость поверхностного течения. Поэтому осредненная по времени мгновенная доплеровская скорость представляет собой сумму радиальной компоненты скорости поверхностного течения Vdr и фазовой скорости брегговских волн сЬг. Следовательно, интерес представляет разность AVD = (VD - Vdr) , которая позволяет судить о собственной скорости рассеивателей и следовательно - о механизме рассеяния.
Отношение измеренных значений AVD к ожидаемой брегговской скорости Сь показано на рис. 5 как функция угла падения при скоростях ветра 4 -
12 м с-1. Положительные величины соответствуют наветренным азимутам, отрицательные - подветренным (в интервале 30° от направления ветра). Приведенные данные хорошо демонстрируют вклад быстрых зеркальных точек при малых углах падения в < 40°. С дальнейшим ростом в отношение А V) / cbr убывает, хотя и не достигает единицы, как этого следовало бы ожидать из классической двухмасштабной модели рассеяния. Кроме того, допле-ровские сдвиги на ГГ-поляризации систематически выше, чем на ВВ-поляризации.
4г
10 20 30 40 50 60 70 80 10 20 30 40 50 60 70 80
Рис. 5. Отношение измеренных значений дкс / сЬг в зависимости от угла падения в для всех
случаев (синие символы) и для выборки «поверхность без обрушений» (красные символы): а -ВВ-поляризация, б - ГГ-поляризация
46
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 4 2015
Вклад обрушений в доплеровскую скорость может быть показан тем же способом, что и вклад в УЭПР (см. выше). Значения AV> /cbr, соответствующие выборке «поверхность без обрушений», изображены на рис. 5 символами «+». Количество таких точек меньше, чем общее число измерений, поскольку информация об обрушениях не всегда могла быть извлечена из видеоданных вследствие неподходящих условий съемки. Как видно на рис. 5, после исключения из выборки обрушений каких-либо существенных изменений в средней мгновенной доплеровской скорости не происходит.
Обсуждение и выводы
Анализ совместных радиолокационных наблюдений и видеозаписей морской поверхности показал, что вклад обрушений ветровых волн в обратное рассеяние 8-миллиметровых радиоволн составляет ~5 - 10% на ВВ-поляри-зации и 10 - 30% на ГГ-поляризации при скорости ветра более 10 мс-1 в диапазоне углов падения 15 - 70°.
Разработанная методика выделения пассивных трассеров на видеозаписях морской поверхности позволила определить скорость поверхностного течения в области облучения РЛС, что дало возможность проанализировать собственные скорости рассеивателей. При малых углах падения ( в < 40° ) средние мгновенные доплеровские скорости практически одинаковы на обеих поляризациях и значительно выше, чем скорость брегговских волн, что свидетельствует о действии зеркального механизма отражения. При умеренных углах падения ( в > 40° ) доплеровская скорость приближается к бреггов-ской, но при наветренных азимутах наблюдения может превышать ее в 1,5 -2 раза, причем это превышение на ГГ-поляризации более существенно, чем на ВВ-поляризации.
Анализ видеоданных показал, что обнаруженные отличия не удается объяснить влиянием только лишь регистрируемых на видеозаписи обрушений. Следует подчеркнуть, что методика обработки видеоданных, использованная в настоящей работе, позволила выделять обрушивающиеся гребни, длина которых больше 0,25 - 0,4 м. Согласно обзору [18], волны короче 25 см обрушиваются без образования пены и не могут быть зафиксированы средствами видеоконтроля. Тем не менее такие события, так называемые микрообрушения, формируют на морской поверхности определенную шероховатость, которая может оказывать влияние на обратное рассеяние радиоволн. Вне зависимости от механизма рассеяния микрообрушения могут иметь скорость выше, чем скорость свободных брегговских волн, в результате чего наблюдаемая доплеровская скорость будет также завышена.
С другой стороны, важным источником коротких волн является генерация паразитной ряби [7, 11]. Фазовые скорости паразитной капиллярной волны и ее несущей равны по определению, но это справедливо только в точках, МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 4 2015 47
где горизонтальная проекция орбитальной скорости несущей волны равна нулю. На других фазах паразитная рябь оказывается сжатой или растянутой орбитальными движениями несущих волн, что ведет к уширению спектра паразитных волн, имеющих одну и ту же фазовую скорость. Поэтому при фиксированном брегговском волновом числе скорость рассеивателей, связанных с паразитной рябью, может быть как больше, так и меньше фазовой скорости свободных брегговских волн. Однако несущие волны с более высокой фазовой скоростью имеют меньшее значение волнового числа и в силу характерной спадающей формы спектра ветровых волн обладают более высокой спектральной плотностью, а значит, имеют более весомый вклад в сигнал по сравнению с более медленными несущими волнами. Кроме того, одним из возможных механизмов увеличения доплеровской скорости может выступать пространственная модуляция брегговской ряби волнами, масштаб которых не превышает размер пятна облучения [19].
Подводя итог, можно сказать, что отличие наблюдаемой средней допле-ровской скорости от фазовой скорости брегговских волн нельзя объяснить только видимыми обрушениями ветровых волн. Завышенные значения до-плеровских скоростей могут быть вызваны влиянием микрообрушений и паразитной ряби, возникновение которой, по сути, является специфическим механизмом обрушения волн малого масштаба.
Авторы выражают благодарность профессору Российского государственного гидрометеорологического университета В.Н. Кудрявцеву за ценные замечания и комментарии.
Работа выполнена при поддержке Российского фонда фундаментальных исследований в рамках проектов 13-05-90912 «мол-ин-нр», 13-05-90429 «укр-ф-а», 14-45-01559 «р-юг-а», 14-05-00876, Правительства Российской Федерации в рамках гранта 11.G34.31.0078, Министерства образования и науки Российской Федерации в рамках ФЦП «Исследования и разработки по приоритетным направлениям развития научно-технологического комплекса России на 2014 - 2020 годы» (уникальный идентификатор проекта RFMEFI57714X0110) и частично поддержана грантом (соглашение от 27 августа 2013 г. № 02.В.49.21.0003 между МОН РФ и ННГУ).
СПИСОК ЛИТЕРАТУРЫ
1. Valenzuela G.R. Theories for the interaction of electromagnetic and ocean waves - a review // Bound. Layer Meteorol. - 1978. - 13. - P. 61 - 85.
2. Bass F.G., Fuks I.M., Kalmykov A.I. et al. Very high frequency radio wave scattering by a disturbed sea surface. Part I: Scattering from a slightly disturbed boundary. Part II: Scattering from an actual sea surface // IEEE Trans. Antenn. Propagat. - 1968. - 16, № 5. - P. 554 - 568.
48
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 4 2015
3. Plant W.J. A two-scale model of short wind-generated waves and scatterometry // J. Geophys. Res. (Oceans). - 1986. - 91, issue C9. - P. 10735 - 10749.
4. Donelan M.A., Pierson W. Radar scattering and equilibrium ranges in wind-generated waves with application to scatterometry // Ibid. - 1987. - 92, issue C5. - P. 4971 - 5029.
5. Kalmykov A.I., Pustovoytenko V. V. On polarization features of radio signals scattered from the sea surface at small grazing angles // J. Geophys. Res. (Oceans and Atmospheres). - 1976. -81. - P. 1960 - 1964.
6. Phillips O.M. Radar returns from the sea surface - Bragg scattering and breaking waves // J. Phys. Oceanogr. - 1988. - 18. - P. 1063 - 1074.
7. Kudryavtsev V., Hauser D., Caudal G. et al. A semiempirical model of the normalized radar cross-section of the sea surface. 1. Background model // J. Geophys. Res. (Oceans). - 2003. -108, issue C3. - P. FET2-1 - FET2-24.
8. Masuko H., Okamoto K., Shimada M. et al. Measurement of microwave backscattering signatures of the ocean surface using X band and Ka band airborne scatterometers // Ibid. - 1986. -91, issue C11. - P. 13065 - 13084.
9. Nekrasov A., Hoogeboom P. A Ka-Band Backscatter Model Function and an Algorithm for Measurement of the Wind Vector Over the Sea Surface // IEEE Geosci. Rem. Sens. Lett. -2005. - 2, issue 1. - P. 23 - 27.
10. Tanelli S., Durden S.L., Im E. Simultaneous Measurements of Ku- and Ka-Band Sea Surface Cross Sections by an Airborne Radar // Ibid. - 2006. - 3, issue 3. - P. 359 - 363.
11. Ermakov S.A., Kapustin I.A., Sergievskaya I.A. Tank study of radar backscattering from strongly nonlinear water waves // Bull. Russ. Acad. Sci. Phys. - 2010. - 74, № 12. -P. 1695 - 1698.
12. Hasselmann K., Raney R.K., Plant W.J. et al. Theory of synthetic aperture radar ocean imaging: A MARSEN view // J. Geophys. Res. (Oceans). - 1985. - 90, issue C3. - P. 4659 - 4686.
13. Plant W.J. The Modulation Transfer Function: Concept and Applications // Radar Scattering from Modulated Wind Waves / Edd. G.J. Komen, W.A. Oost. - Netherlands, Springer, 1989. - P. 155 - 172.
14. Kudryavtsev V., Hauser D., Caudal G. et al. A semiempirical model of the normalized radar cross section of the sea surface. 2. Radar modulation transfer function // J. Geophys. Res. (Oceans). - 2003. - 108, issue C3. - P. FET3-1 - FET3-16.
15. Monahan E.C., Woolf D.K. Comments on "Variations of Whitecap Coverage with Wind stress and Water Temperature" // J. Phys. Oceanogr. - 1989. - 19. - P. 706 - 709.
16. Миронов А.С., Дулов В.А. Статистические характеристики событий и диссипация энергии волн при обрушении // Экологическая безопасность природных и шельфовых зон и комплексное использование ресурсов шельфа. - Севастополь: МГИ НАН Украины, 2008. - Вып. 16. - С. 97 - 115.
17. Dulov V.A., Kudryavtsev V.N., Bolshakov A.N. A Field Study of White Caps Coverage and its Modulations by Energy Containing Waves // Gas Transfer at Water Surface, Geophysical Monograph Series / Edd. M.A. Donelan, W.M. Drennan, E.S. Saltzman et al. - AGU, 2002. -127. - P. 187 - 192.
18. Babanin A.V. Breaking of ocean surface waves // Acta Phys. Slov. - 2009. - 59, № 4. -P. 305 - 535.
19. Ermakov S.A., Sergievskaya I.A., Schegolkov Y.B. Laboratory study of strong modulation of radar signals due to long waves on the water surface covered with a surfactant film // Radi-ophys. Quant. Electron. - 2002. - 45, № 12. - P. 942 - 957.
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 4 2015
49
Influence of breaking wind waves on a millimeter-wave radar backscattering by the sea surface
Yu.Yu. Yurovsky I.A. Sergievskaya 2, S.A. Ermakov 2'3'4, B. Chapron 5, I.A. Kapustin 2'3'4, O.V. Shomina 2
1 Marine Hydrophysical Institute, Russian Academy of Sciences, Sevastopol, Russia e-mail: yyyurovsky@gmail com 2 Institute of Applied Physics, Russian Academy of Sciences, Nizhni Novgorod, Russia 3 N.I. Lobachevsky State University of Nizhni Novgorod, Russia 4 Volga State University of Water Transport, Nizhni Novgorod, Russia 5 Russian State Hydrometeorological University, Saint-Petersburg, Russia
Results of field studies of radar backscattering from the sea surface on 8 mm radio wave length are represented. Radar observations are accompanied by synchronous video registration of the illuminated area with the purpose to define parameters of breaking of wind waves influencing backscattering intensity and to measure velocity of surface current producing the signal frequency Doppler shift. It is shown that contribution of breaking waves to NRCS depends on the wind speed and polarization and does not exceed 20 - 30%. Average values of an instantaneous Doppler shift (not related to a surface current) are 1.5 - 2 times as much as the expected Doppler shift resulted from free Bragg waves. Analysis of video data shows that the obtained discrepancy cannot be attributed to influence of visible wave breakings only. Possible reasons of the observed effect can consist in the following: spatial modulation within the radar footprint, influence of bound parasitic capillary ripples and micro-scale breaking undetectable on the sea surface.
Keywords: radar backscattering, millimeter waves, Doppler shift, surface current, wind wave breaking, parasitic ripples.
50
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 4 2015
