Институт радиотехники и электроники им. В.А. Котельникова РАН
Фрязинский филиал
Исследования параметров распространения и рассеяния сантиметровых радиоволн на радиозатменных трассах спутник-спутник применительно к задаче идентификации гидрометеорных образований
Е.В. Сухонин
Аннотация
Выполнены оценки ослабления А и деполяризации XPD радиоизлучения на волне 2 см при его прохождении по радиозатменной трассе спутник-спутник с гидрометеорами при заходе одного из спутников в область радиотени Земли с целью изучения возможности идентификации дождей. Результаты расчетов ослабления радиоволн в дожде сравнивались с экспериментальными данными, полученными на радиозатменных трассах аппарат МИР-геостационарный спутник. Показано, что в качестве критерия идентификации можно использовать характерный вид зависимости ослабления А от высоты лучевой линии в тропосфере Н. Провал на кривой A(H) может служить признаком существования либо обложного дождя с зоной таяния осадков, либо слоистых протяженных облаков. Различать эти гидрометеорные образования можно по поведению сигнала при H>0 km. Получено соотношение, связывающее между собой А и XPD в дождях на волне 2 см, позволяющее идентифицировать их более надежно. Установлено, что поляризационный метод дает возможность обнаруживать в атмосфере грозовые разряды.
Ключевые слова: радиоволны, радиозатменные трассы спутник-спутник, мониторинг параметров атмосферы, идентификация дождей, ослабление.
В связи с появлением систем глобального позиционирования и низкоорбитальных спутников начинают развиваться работы по мониторингу метеопараметров атмосферы с помощью радиозатменного метода [1,2]. Хотя этот метод уступает по точности наземным методам, он является глобальным.
Из анализа [3,4] опубликованных работ и приведенной в них библиографии следует, что радиозатменный метод может применяться для определения и мониторинга содержания в атмосфере водяного пара и водности облаков по поглощению сантиметровых радиоволн на трассах спутник-спутник при заходе одного из спутников в область радиотени Земли. Что касается дождей, то в ряде случаев их идентификация и определение их интенсивности по ослаблению волн затруднительны, поскольку ослабление сантиметровых волн в дождях и мощных облаках часто бывает одинаковым. Поэтому требуются соответствующие теоретически обоснованные обработка и интерпретация экспериментальных данных по ослаблению в дождях радиоволн на спутниковых радиозатменных трассах, а также правильный выбор моделей дождя и методов анализа расчетов, чтобы можно было разделить вклады дождя и облаков без дождя в ослабление радиоволн.
Целью работы являлось изучение возможности идентификации дождей и получения их отличительного признака при сравнении с облаками на основе исследований теоретических расчётных зависимостей ослабления А и деполяризации XPD сантиметровых волн в дождях разного типа на затменных трассах спутник-спутник от расстояния Н между земной поверхностью и трассой распространения радиоволн. Предварительные оценки показывают, что эти зависимости могут иметь "характерный" вид для дождей разного типа и по этой причине могут использоваться для их идентификации. В [5] было показано, как можно идентифицировать облака и оценивать поглощение в них радиоволн, однако вопрос об идентификации дождей остался открытым.
В работе анализировались результаты расчетов зависимостей А(Н) и XPD(Н) в дождях на волне 2 см. При расчетах ослабления А капли дождя рассматривались как сферические, а при расчетах деполяризации XPD как сфероиды, наклоненные по отношению к линии горизонта [6,7]. При расчетах использовались три модели структуры дождя (М1, М2 и М3) [8,9], при этом только одни водные дождевые капли в чистом виде существуют в самом нижнем слое атмосферы толщиной 4 км, верхняя граница которого соответствует уровню замерзания воды 0°С.
В случае всех трёх моделей над нижним слоем, состоящим из одних дождевых капель, находится слой облаков толщиной 2 км с содержанием воды 500 г/м2. В модели М2, в отличие от М1, слой облаков включает также дождевые ледяные капли. Модели М1 и М2 соответствуют обложным дождям. Модель М3, которая более всего соответствует сильным ливням конвективного происхождения ("глубокое дождевое облако с гидрометеорами в смешанной фазе" [9] ), включает в себя ещё два более высоких слоя - один с ледяными кристаллами, замёрзшими каплями дождя и жидкими каплями дождя, и другой, самый высокий слой, только с ледяными кристаллами. Кроме того, слой облаков в модели М3 включает также жидкие и замёрзшие капли дождя
В случае моделей М1 и М2 обложные дожди представляют собой однородные образования с интенсивностью дождя R = 0.5…1 мм/ч, ограниченные сверху уровнем нулевой изотермы атмосферы и занимающие в пространстве десятки и даже сотни километров. Ещё одна особенность таких дождей (модель М2) - это существование в них зоны таяния осадков, находящейся ниже уровня нулевой изотермы атмосферы и достигающей по размеру нескольких сотен метров. Зона таяния осадков представляет собой мокрый снег и снег с дождём, ослабление в которых в два-три раза больше, чем в дожде с такой же интенсивностью [7].
Результаты расчета представляют собой имитацию радиопросвечивания дождей на волне 2 см при изменении минимальной высоты лучевой линии над земной поверхностью Н от НZ до Н=0 км. Результаты подробного анализа зависимости А(Н) для всех трех моделей дождя М1, М2 и М3 приведены в [10].
Судя по А(Н) на радиозатменных трассах в [5], характерный вид экспериментальной кривой А(Н) свидетельствует о наличии на радиозатменной трассе ливневого дождя, описываемого моделью дождя М3. Следовательно, такая характерная особенность А(Н) может служить одним из критериев идентификации ливневых дождей.
Зависимость A(H) для обложных дождей (модель М2) с провалом кривой A(H) при H ? 4 км обусловлена зоной таяния осадков. В модели М1 зависимость A(H) полностью совпадает с A(H) для модели М2 при H ? 3.5 км, при этом провал в кривой A(H) при H > 3.5 км отсутствует, т.к. в этом случае зона таяния осадков не существует. Расчёты зависимости A(H) проводились в предположении, что длина трассы распространения радиоволн в обложном дожде равна 150 км, дождь является однородным с R = 1 мм/ч, толщина зоны таяния осадков 500 м, а коэффициент ослабления в ней СМ радиоволн втрое превышает коэффициент ослабления в дожде.
Зависимость A(H) с провалом можно рассматривать как характерную особенность ослабления в слабых обложных дождях. Однако использовать её в качестве единственного критерия для идентификации обложных дождей типа М2, строго говоря, нельзя, поскольку аналогичные зависимости A(H) с провалом могут существовать и в случае слоистых протяжённых облаков [5].
С целью повышения надежности идентификации дождей на радиозатменных трассах следует одновременно с А измерять деполяризацию XPD, учитывая, что величины XPD могут заметно превышать величины ослабления А, если преобладают гидрометеоры в виде ледяных кристаллов. Анализ экспериментальных и теоретических результатов [11] показал, что между А и XPD существует функциональная связь, не зависящая от интенсивности дождя и длины трассы распространения радиоволн.
С учетом того, что в рассматриваемом нами случае f=15 ГГц (л=2 см), угол поляризации волны 45о, угол возвышения трассы 0о, то в соответствии с рекомендациями МККР [12]
XPD = 35.3 - А (1)
В случае другой частоты и поляризации величина первого члена в правой части (1) изменится. Равенства типа (1) можно использовать для идентификации дождя на радиозатменных трассах. При одновременных измерениях ослабления А и деполяризации XPD волн возможны три случая: (а) левая и правая части в (1) равны друг другу, (б) левая часть больше правой и (в) правая часть больше левой. Первый случай означает, что на затменной трассе идет дождь, в котором деполяризация радиоволн обусловлена только каплями дождя. Во втором и третьем случаях на величину XPD в дожде влияет слой ледяных кристаллов, вызывая, в частности, дополнительный рост деполяризации при малой величине ослабления А.
По резкому изменению XPD во время грозовых разрядов можно также определять области атмосферы с грозовыми разрядами. Проведенные исследования показали, что измеряемые ослабление и деполяризация сантиметровых радиоволн на радиозатменных трассах в гидрометеорных образованиях позволяют идентифицировать дожди разных типов, облака и ледяные кристаллы.
На рис.1 представлена геометрия задачи. Земная поверхность изображена в виде полуокружности с радиусом а , в точках S1 и S2 расположены спутники с передающей и приёмной аппаратурой, OD = OC = a , DB = h , CA = H - высота лучевой линии S1ABS2 над земной поверхностью. Поскольку согласно [5] угол рефракции радиоволн в тропосфере мал, то лучевая линия S1ABS2 практически представляет собой прямую линию.
Рис.1. Геометрия задачи.
Длина пути радиоволн при их распространении на трассе AB в общем случае зависит от h и определяется приближённо выражением
Эффективная длина в тропосфере трассы спутник-спутник ~500 км. Диаметр пучка радиоволн с длиной 2 см, формирующего луч, в соответствии с формулой Френеля составляет ~200 м. Размеры зоны ливневых и обложных дождей обычно не превышают соответственно 20…30 и 100…150 км. Поэтому решаемая в данной работе прямая задача расчёта ослабления радиоволн на затменной трассе может рассматриваться как плоско-слоистая, а не сферически-слоистая.
При распространении на радиозатменной трассе спутник-спутник радиоволны испытывают поглощение парами воды, кислородом и облаками, ослабляются в осадках и вследствие эффектов рефракции. В [5] был выполнен анализ всех факторов, влияющих на поглощение волн на указанных трассах, за исключением дождей.
Обычно используются три модели, характеризующие форму капель дождя различной степени сложности [6]. В самой простой модели предполагается, что капли дождя являются сферическими, поэтому применима теория Ми, описывающая рассеяние радиоволн на сферических частицах.
Во второй модели капли являются сфероидальными и образуют анизотропную среду. В отличие от первой модели вторая модель позволяет учитывать зависимость ослаблений рассеяния от направления распространения радиоволн.
Третья модель учитывает в некоторой степени динамику капель дождя и связанную с ней деполяризацию радиоволн. При падении капли сплющиваются, причём вследствие существования в земной атмосфере высотного градиента скорости ветра сфероидальные капли дождя наклонены по отношению к горизонтали [7]. Поэтому наблюдается деполяризация сигнала в дожде как результат дифференциального ослабления и дифференциального сдвига фазы между компонентами сигнала с векторами поляризации параллельными большой и малой осям сфероида. Хотя согласно [13], вытянутые и сплющенные сфероиды могут иметь место с равной вероятностью и даже возможны колебания капель, при теоретическом исследовании деполяризации волн в дожде предполагается, что преобладают сплющенные сфероиды [6]. Первая и вторая модели являются частными случаями третьей модели.
При расчётах чаще всего используется распределение капель дождя по размерам Лоуса-Парсонса [14]. Это позволяет более точно сравнивать результаты расчётов разных авторов.
На рис.2 и 3 приведены результаты расчётов зависимости ослабления A радиоволн с длиной волны 2 см в дождях на радиозатменных трассах спутник-спутник от H соответственно для случаев дождей конвективного происхождения (модель М3) и обложных дождей (модели М1 и М2). Приводимые на обоих рисунках пунктирные кривые 1 заимствованы из работы [5]; они характеризуют молекулярное поглощение в водяном паре в чистой стандартной атмосфере на трассах без дождей и облаков. Все величины ослабления A, лежащие ниже кривой 1, связаны с существованием в атмосфере гидрометеорных образований.
На рис.2 кривые 2, 3, 4 и 6 соответствуют интенсивностям ливневых дождей R = 5, 25, 50 и 100 мм/ч. Представленная в виде кружков на рис.2 экспериментальная кривая 5 A(H) была получена из работы [5] на основе анализа результатов двух измерений A(H) на радиозатменной трассе аппарат Мир-геостационарный спутник. При этом учитывалось, что в [5] из измеряемых величин A(H) вычитались величины ослабления в кислороде и вследствие рефракции. (В [15] показано, в принципе рефракционное ослабление можно измерять, используя сигналы системы глобального позиционирования GPS). Это позволило нам провести сравнение экспериментальной кривой 5 с теоретическими кривыми на рис.2, при расчёте которых данные факторы ослабления не учитывались. Кривые 2, 3, 4 и 6 вместе с кривой 1, описывающей молекулярное поглощение в парах воды, характеризуют дополнительный вклад в ослабление в водяном паре, вносимое ливневыми дождями. Начиная с H і 5 км, этот вклад стремится к нулю, т.к. в верхней части дождей преобладают ледяные кристаллы и капли, которые очень слабо поглощают cантиметровые радиоволны [16]. Судя по рис.2 и анализу атмосферного поглощения на радиозатменных трассах в [5], характерный вид экспериментальной кривой A(H) свидетельствует о наличии на радиозатменной трассе ливневого дождя, следовательно, такая характерная особенность может служить одним из критериев идентификации ливневых дождей. Результаты расчётов A(H) на рис.2 в принципе позволяют оценить интенсивность дождя R по измеряемой зависимости A(H). (Например, приводимая на рис.2 экспериментальная зависимость A(H) даёт R = 50…100 мм/ч). Однако для повышения надёжности такой оценки требуются дополнительные критерии идентификации дождей.
деполяризация радиоизлучение спутник гидрометеор