Влагосодержание тропосферы над пунктом наблюдения Усть-Баргузин

При использовании высокоточных спутниковых навигационных систем (GPS, ГЛОНАСС, GALILEO) обнаруживаются погрешности позиционирования, связанные с задержкой навигационных сигналов в ионосфере, стратосфере, тропопаузе и тропосфере. Задержка сигналов в нижней части атмосферы, которая является нейтральной, называется тропосферной задержкой. В случае зенитного направления сигнала задержка называется тропосферной зенитной задержкой (ZTD). Скорость распространения сигнала cзависит от коэффициента преломления среды nв нейтральной атмосфере: c= c0/n, где c0 - скорость света в вакууме.

Индекс рефракции в тропосфере (N) подчиняется закону аддитивности и записывается в виде суммы «сухой» и «влажной» компонент [Smith, Wein- traub, 1953]:

где п - коэффициент преломления; к = 77,6 - первая рефракционная константа, К/мбар (1 мбар = 102 Н/м2); к2 = 3,73-105 - вторая рефракционная константа, К2/мбар; Т - абсолютная температура, К; р - атмосферное давление, мбар; e- упругость водяного пара, мбар.

Здесь индекс рефракции для сухого воздуха, зависящий от изменений температуры и дав-ления воздуха, - индекс рефракции для водяного пара.

На рисунке 1 показан годовой ход приземных параметров N, NDRYиNWETдля пос. Усть-Баргузин за 2017 г. «Влажная» компонентаNWETвносит существенно меньший вклад в определение индекса рефракции (коэффициент преломления), особенно в зимний период. «Сухая» компонентаNDRYсоставляет значительную часть индекса рефракции.

Рис. 1. Значения N, NDRYи NWETпо месяцам 2017 г. для пос. Усть-Баргузин

Формула (1) определяет индекс рефракции Ы, если известны температура, давление и упругость водяного пара в пункте наблюдения. Из формулы для N следует, что 2ТЕ> также является суммой «сухой» ^НП) и «влажной» (ZWD) составляющих. Полная задержка пути сигнала от ОР8-спутника до приемной антенны равна разности геометрического расстояния и реального пути сигнала в атмосфере:

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

где к3 - высота расположения пункта наблюдения над уровнем моря, №ОКТ,

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

Для расчета параметров 2ИЕ> и ZWD по приземным метеоданным наиболее точной и обоснованной является модель Саастамоинена [Баа81:атотеп, 1972]:

Здесь (ц, ) 1 0, 00266 cos 2ц 0, 00028 SS

f h =- -h, где ц - географическая широта расположения приемника в градусах; hS- высота приемника над уровнем моря, км. Для пункта Усть-Баргузин знаменательf (ц,hS)в формулах (4) и (5) равен приблизительно 1,0005, и в расчетах ZHD и ZWD этим множителем можно пренебречь.

Суммарный водяной пар и осаждаемая вода в тропосфере

Текущее значение «сухой» компоненты ZHD может быть достаточно точно рассчитано с использованием значений приземного давления по формуле (4), учитывая приближение /(ф, к3) » 1. При обработке первичных измеренных ОРБ-данных учитываются поправки, связанные с неоднородностью нижней части атмосферы в виде тропосферной зенитной задержки ZTD [Каш-кин, 2013]. В этом случае для выявления текущего значения «влажной» составляющей ZWD достаточно вычесть значение «сухой» компоненты из измеренного значения ZTD: ZWD = ZTD -- ZHD . При наличии приземных метеорологических данных значение ZWD рассчитывается непосредственно по формуле (5). В качестве примера на рис. 2 показаны графики годовых изменений полной тропосферной зенитной задержки ZTD и ее компонент ZHD и ZWD, рассчитанных для пункта наблюдения Усть-Баргузин за 2017 г.

Рис. 2. Годовой ход рассчитанных по метеоданным значений 2ТБ, ХИБ и для пункта наблюдения Усть-Баргузин за 2017 г.

Выражение входящее в формулу (1), вставим в фор-мулу (3) для параметра ZWD:

В работах [Диагностика содержания атмосферного … , 2015; GPS meteorology: remote … , 1992; Geodesybyradiointerferometery … , 1985], посвященных тематике влагосодержания, описывается параметр средневзвешенной температуры по упругости водяного параTmнад пунктом наблюдения:

Для нормальных атмосферных условий, учитывая формулу (7) и соотношение , де сWV - плотность водяного пара, кг/м3, RW» 4, 6 - газовая постоянная для водяного пара, Дж·К-1кг-1, формулу (6) можно записать в виде

Здесь - суммарный водяной пар в килограммах в расчете на единицу площади. Осаждаемая вода определяется в виде колонки жидкой воды PW, м. Имеет место соотношениеPW = IWV/ с, где с - плотность жидкой воды, кг/м3.

Формула (8) дает возможность определять над пунктом наблюдения суммарный водяной пар (1ЖУ), а также потенциальный уровень осаждаемой воды в тропосфере (РЖ). В теплое время года количество суммарного водяного пара имеет максимальное значение в тропосфере. Измеряемые величины 2ТЕ> практически пропорциональны в летнее время уровню влагосодер-жания над пунктом наблюдения. Поэтому имеется возможность использования в задачах метеорологии значений тропосферной задержки, получаемых в результате обработки первичных данных постоянных ОРБ/ГЛОНАСС-измерений, для изучения вариативности влагосодержания нижней части атмосферы.

Вертикальные профили температуры воздуха

Наиболее громоздким при нахождении значений суммарного водяного пара 1ЖУ по данным «влажной» компоненты тропосферной задержки является нахождение значений средневзвешенной температуры по упругости водяного пара по приземным данным температуры воздуха. Для расчета по формуле (7) необходимо знание данных по температуре и упругости водяного пара в вертикальном столбе.

Помимо запусков метеозондов для определения вертикальных профилей температуры, давления и относительной влажности воздуха в настоящее время активно используются современные космические технологии с помощью инструментального комплекса ЛТОУБ на полярно-орбитальных спутниках КОЛЛ [Кашкин, 2013; Кашкин, Владимиров, Клыков, 2014]. При отсутствии облачности сканирующий спектрофотометр, входящий в комплекс ЛТОУБ, позволяет определять профили метеорологических данных до высоты 40 км над поверхностью Земли (http://ready.arl.noaa.gov).

На рисунке 3 а, б показаны сравнительные графики измеренных вертикальных профилей температуры, полученных для разных сезонных условий над пунктом наблюдения Усть-Баргузин по данным запусков радиозондов Аэрологической службой Росгидромета и по данным зондирования спутниками КОЛЛ. Наблюдается достаточно хорошее соответствие полученных результатов, которые могут быть использованы для целей определения рефракционных параметров тропосферы.

На рисунке 4 приведены линейная (Тт = 43 + 0,83- Т) и квадратичная (Тт = 431 - 2,05- Т + 0,0053-Т2) регрессии соотношений приземной температуры (Т) и средневзвешенной температуры (Тт) по упругости водяного пара. Регрессии выявлены по данным 2715 вертикальных зондирований атмосферы за период с 2013 по 2017 г. и рассчитаны по формуле (7) для пос. Усть- Баргузин. Показано, что средневзвешенная температура до высоты нижней части стратосферы может быть определена по данным приземной температуры воздуха.

Рис. 3. Сравнения высотных температурных профилей, полученных методами радиозондирования и спутникового зондирования КОДА 6 марта (а) и 6 июня 2018 г. (б)

Рис. 4. Соотношение параметров Т и Тт, выявленных по данным метеорологических зондирований над пос. Усть-Баргузин с 2013 по 2017 г.

Численные результаты