Оценка величины зенитной задержки радионавигационных сигналов в юго-восточной части Охотского моря по данным метеоизмерений вблизи навигационного приемника

УДК 537.876.23

Ф.Н. Захаров, С.Г. Госенченко, М.В. Крутиков

Оценка величины зенитной задержки радионавигационных сигналов в юго-восточной части Охотского моря по данным метеоизмерений вблизи навигационного приемника

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

Ключевые слова: псевдодальность, задержка навигационного сигнала, тропосфера, индекс преломления, среднеквадратическая ошибка.

В статье рассматривается задача расчёта поправки за счёт влияния тропосферы при оценке псевдодальности в навигационных системах. При работе навигационного приёмника в автономном режиме одним из основных факторов, влияющих на точность определения координат, является задержка навигационного сигнала в тропосфере [1], которая возникает вследствие отличия скорости радиоволны в тропосфере от скорости света.

В случае, когда угол возвышения навигационного космического аппарата (НКА) над горизонтом составляет 90 град, а высота приёмника радионавигационного сигнала составляет hs, для расчета тропосферной путевой задержки Tj используется вертикальный профиль индекса преломления N(h) [2]:

H

tt = J 10-6 N

hs

где H - высота над уровнем моря в метрах, на которой значение индекса преломления можно считать равным нулю. Задержка, вычисленная по формуле (1), называется зенитной. При меньших углах возвышения задержка увеличивается за счёт удлинения пути в тропосфере и определяется с помощью так называемой функции отображения [2, 3].

Известно [2, 3], что в практике использования навигационной аппаратуры потребителя (НАП) зенитная задержка сигнала оценивается по приземным значениям метеоэлементов (давление, влажность, температура) в месте расположения НАП. Алгоритмы оценивания зенитной задержки основаны на знании моделей вертикального профиля метеоэлементов [4] или профиля N(h) [5-7]. Уточнение расчетов с их использованием ведет к усложнению моделей, прежде всего к увеличению их числа в соответствии с особенностями местных условий. Например, в численной модели распределения температуры и влажности для Северного полушария [8] для приемлемого описания случайных вариаций вертикальных профилей этих метеоэлементов применено 20 локальных статистических моделей зимой и 17 - летом [9].

Целью настоящего сообщения является оценка величины и флуктуаций зенитных задержек в нижнем слое атмосферы 0...100 км в юго-западной части Охотского моря, полученных расчетным путем по вертикальным профилям индекса преломления Nj(h) (j - номер профиля в объеме экспериментальных данных), и задержек, рассчитанных по моделям N(h) [6].

Экспериментальные данные. В качестве исходных данных для расчетов экспериментальных профилей Nj(h) были взяты массивы результатов аэрологического зондирования трех расположенных выше 45° с.ш. аэрологических станций (АС) района Охотского моря - АС «Поронайск», АС «Уруп», АС «Южно-Курильск» - по 20 зондирований в каждом пункте, полученных как в летнее время (июнь-август), так и зимнее время (январь-февраль) в разные сроки зондирования в 1973, 1974 и 1975 гг. Высоты расположения АС hs над уровнем моря составляли: для АС «Уруп» - 70 м, для АС «Поронайск» - 2 м, для АС «Южно-Курильск» - 40 м.

Для расчёта задержек во всей толще тропосферы тТ. и сопоставления результатов по данным различных АС был проведен пересчет измерений ^-ф,), выполненных на высотах расположения АС Н,, к значению на уровне моря N0.

В связи с отсутствием данных зондирования до высот 100 км оценки зенитных задержек в настоящей работе определяются суммированием задержек в слое 0...30 км, полученным по профилям стандартного аэрологического зондирования, и задержки в слое 30... 100 км, вычисленной в соответствии с Эталонными стандартами атмосферы [10].

Оценки задержек для слоя атмосферы 0-30 км. Экспериментально измеренный профиль Щ(Н) представляет собой ряд значений индекса преломления, вычисленных по формуле Смита-Вейнт-рауба [11] с использованием метеовеличин, измеренных аэрологическим зондом на высотах к{.

N. )=НА

3 Т

РОн) +4810 е

т

где Т(Нг) - измеренное значение температуры на высоте Н. в градусах Кельвина; Р(Нг) - давление в мбар; е(Н) - влажность в мбар.

Расчёт задержки сигнала по экспериментальному профилю N.(Н) осуществлялся по формуле (1) при линейно-ломаной аппроксимации промежуточных значений индекса преломления. Кроме того, расчеты хТз- с использованием формулы (1) проведены также для различных математических моделей вертикального профиля индекса преломления при условии измерений метеовеличин вблизи НАП. В качестве математического описания профиля можно рассмотреть модель профиля N(0) в виде ряда гармоник [12]. Однако в этой модели нет привязки к наземным измерениям N(0), а используются только многолетние усреднённые данные. Мы ограничимся в расчетах двумя наиболее распространёнными моделями вертикального профиля N(0) из работы [6].

Экспоненциальная модель. Основная идея этой модели заключается в предположении о том, что в среднем вертикальный (высотный) профиль индекса преломления над земной поверхностью для всей толщи атмосферы описывается одной общей экспоненциальной зависимостью [6].

N

где Н - высота над уровнем моря.

Модель имеет два параметра. N0 и Ь = 1п(^/ЭД), где N0, N - значения индекса преломления на уровне моря и на высоте 1 км соответственно. Распределение среднемесячного параметра N по разным районам земного шара для разных сезонов года представлено в виде мировых карт [13]. Необходимые для расчетов таблицы соотношений между N0, N1, и Ь приведены в [6].

Модель модифицированного эффективного радиуса Земли. Эта более сложная модель индекса преломления предполагает следующее [6].

1. В пределах высот до 1 км над поверхностью Земли соблюдается концепция эффективного радиуса Земли, при которой существует линейная высотная зависимость ^Н) до 1 км.

N

где Д^ = -N1) = 7,32ехр(0,005577#5) - разность значений индекса преломления у поверхности

Земли и на высоте 1 км, определенная усреднением данных ежедневных наблюдений за 6-8 лет.

2. Выше 1 км индекс преломления изменяется по экспоненциальному закону до некоторого значения N9 на высоте 9 км над уровнем моря

N

где Ь = [(8 - % )1п(Щ/Щ)]-1.

3. Выше 9 км ^профиль описывается единой экспонентой

N

В соответствии с рекомендациями [10] при расположении НАП на широтах более 45° с.ш. значение N9 составляет летом 103,2 №ед., а зимой 99,8 №ед.

Для удобства будем называть такую модель трехэлементной [6].

Расчет с использованием экспоненциальной модели. Для определения параметра Ь были использованы карты из [14] с изолиниями параметров N0 и N - N1. В исследуемом районе для февраля и августа использовались значения параметров (в №ед.) N0 = 310, (N0 - N1) = 40 и N0 = 340, (N0 - N1) = 50 соответственно, в результате значения Ь составили 0,1382 км 1 для зимы и 0,1591 км 1 для лета.

Мгновенное значение N0j каждого профиля должно быть найдено из условия принадлежности этого значения модельному вертикальному профилю N(h) по формуле N0j = N(hs)exp(bhs).

Расчет зенитной задержки от уровня моря до высоты H = 30 км по экспоненциальной модели проводился с использованием параметров N0j и b и соотношения (1) по формуле:

6 3° 10-6 No j tt = 10-6 N0 j f exp

0 b

Расчет с использованием трехэлементной модели. Для применения трехэлементной модели (4)-(6) необходимо определить приповерхностное значение N0j по значению Ns на некоторой высоте, из уравнения:

Ns = N0 j + hs AN0 j,

где AN0. = —7,32hs exp(0,005577N0.). Вычисления AN0j, с учетом линейной связи AN0j и N0j в интервале значений от 280 N-ед. до 350 N-ед., выполнялись по формуле

AN0j = hs (0,2375N0j -32,108).

Отсюда

N0j = Ns -32,108hs .

J 1-0,2375hs

Таким образом, расчет задержки ttj с помощью трехэлементной модели сводится, во-первых, к приведению индекса преломления Ns к приповерхностному N0j по формуле (8); во-вторых, к последующему интегрированию на трех участках толщи атмосферы от высоты уровня моря до высоты H = 30 км с использованием соотношений (4)-(5) по формулам:

( 1 9 30 Л

tTj ttj = 10-6I f Nj

V 0 1 9 )

где Ij = 0,5((V0 j + Nj, IIj = 8(N1j - N9 ))ln(N1j- /N9), III = N9 {1-exp [-0,1424(H - 9)]}/0,1424.

Результаты расчета задержек в слое 0-30 км. Для экспериментальных профилей и рассматриваемых моделей по полученным рядам значений ttj рассчитывались среднее значение Тт и СКО. Результаты расчетов использовались для построения гистограмм (полигонов) ttj. Кроме того, оценивались отклонения Attj как разница между экспериментальными и модельными значениями ttj.

Статистические характеристики зенитной задержки ttj (среднее %т и СКО флуктуаций стТг) и отклонений Attj (среднее Атт и СКО флуктуаций аАТт ) приведен^1 в табл. 1 и 2 соответственно.

Гистограммы ttj представлены на рис. 1 и 2 для лета и зимы соответственно, а гистограммы Attj для экспоненциальной и трехэлементных моделей - на рис. 3 и 4 для лета и зимы соответственно.

Таблица 1 показывает, что в слое 0-30 км экспериментальные данные о ttj проявляют сезонную зависимость средней задержки, имеющей летом величину 246,4 см, зимой 231,7 см, и СКО флуктуаций задержки летом 5,6 см, а зимой 2,1 см. В соответствии с табл. 2 лучшей моделью вертикального профиля индекса преломления по критерию минимума полной ошибки стошибки является трехэлементная модель с ошибкой стошибки 5,1 см и 1,9 см летом и зимой соответственно.

Таблица 1

Статистические характеристики зенитных задержек TTj в слое 0-30 км_

№ п/п Форма представления вертикального профиля ЩИ) Среднее значение Атт , см СКО стТт, см

лето зима лето зима

1 Экспериментальная оценка по данным аэрологического зондирования и линейно-ломаной аппроксимации профиля ЩИ) 246,4 231,7 5,6 2,1

2 Экспоненциальная модель 208,6 220,8 4,8 3,2

3 Трехэлементная модель 244,5 232,1 2,7 1,4

Таблица 2

Статистические характеристики отклонений Дту модельных оценок тт от экспериментальных в слое 0-30 км

№ п/п Вертикальный профиль Среднее значение Дтт, см СКО флуктуаций стДхт , см Погрешность оценок модели, см ¡2 2 с°шибки =уДтт +°ДхТ

лето зима лето зима лето зима

1 Экспоненциальная модель -39,8 -11,9 6,6 3,2 40,3 12,3

2 Трехэлементная модель -1,9 0,5 4,7 1,8 5,1 1,9

25 * 20

рс

Эксперимент а льные д анные Экспоне ндаа лвна л модель Трехэлементная мод ель

.Л

Г+ г:.' Нмк А.

г  £- \ V V  • А * • г -

200

210

250

220 230 240

Задержка,™

Рис. 1. Вероятность задержек xTJ в слое тропосферы 0-30 км, лето

260

35

10 -

25 -

i-

с 20 -

и 15 -

а 10 -

5 -

0 -

♦ Эксперимент альные д анные Я Экспоненциальная модель —д— Трехэлементная мод ель

1V

п

/* ! д /

А V ^ { к

* 1

210

220

230

240

Задержка, см

Рис. 2. Вероятность задержек сигнала т^. в слое тропосферы 0-30 км, зима

Гистограммы на рис. 1 и 2 поясняют указанные оценки. Так, экспериментальные значения задержки в летний сезон принимают практически равновероятно значения от 238 до 258 см, в то время как модельные оценки задержки по экспоненциальной модели имеют существенное смещение в область значений от 202 до 222 см с явным максимумом вероятности в области 205-207 см. Зимой смещение расчетов для экспоненциальной модели уменьшается, а значения задержек распределены в области 215-228 см. В то же время при трехэлементной модели вертикального профиля зимой и летом модельные оценки вписываются в пределы существования экспериментальных значений, отличаясь более узким распределением в сравнении с экспериментальными данными.

25

Я

I 15

5 0

-60

Я Экспоненциальная модель —*— Трехэлементная модель t

/1

р   \ / i

я щ I 1 1 k J Ь

-.Г % Ы4 Й 7 а

-50

40 -30 -20 -10 Отклонение оценки, сы

10

Рис. 3. Вероятность отклонений Дтп модельных и экспериментальных задержек в слое 0-30 км, лето

30 25 Й 20

в 15

рс

10 5 0

—Ш— Экспоненциальная модель д Трехэлементная модель

Л/\ / Л \

1 * V

1 \

/ / --   -У  1 МИ -4

-20

-15

-10 -5 О

Отклонение оценки, см

10

Рис. 4. Вероятность отклонений Дтт. модельных и экспериментальных задержек в слое 0-30 км, зима

Распределения отклонений Дт. в летнем и зимнем сезонах (рис. 3 и 4) указывают на наибольшую близость расчетов по трехэлементной модели к экспериментальным данным.

Оценки задержек в слое атмосферы 30-100 км. С использованием приведенных в [10] значений метеовеличин и их градиентов были рассчитаны средние профили N(0) для каждого сезона года. Затем задержки в слое атмосферы 30-100 км оценивались по формуле (1), в которой нижний предел интегрирования равнялся 30 км, а верхний - 100 км. В результате расчетов оказалось, что средняя задержка в этом слое атмосферы, по данным использованных для расчетов АС в широтах выше 45°с.ш., составила 3,6 летом и 3,0 см зимой. В связи с тем, что указанная рекомендация не определяет флуктуации задержки, то они были оценены эмпирически величиной 10% от величины средней задержки и составили 0,36 и 0,30 см для лета и зимы соответственно. Заметим, что 95% средней задержки определяется нижним участком атмосферы до высоты 50 км, что подтверждается работой [15].

Результаты оценки зенитной задержки в слое 0-100 км для различных сезонов года приведены в табл. 3. Таким образом, в юго-восточной части Охотского моря среднее значение задержки летом может составлять 250 см, а зимой 234,7 см при СКО флуктуаций 5,61 и 2,12 см соответственно.

Таблица 3

Статистические характеристики зенитной задержки тт в слое 0-100 км, см

№ п/п Форма представление вертикального профиля N(0) Среднее значение тт СКО стт тт

лето зима лето зима

1 Экспериментальная оценка по данным аэрологического зондирования и линейно-кусочной аппроксимации профиля N(0) в слое 0-30 км 246,40 231,70 5,60 2,10

2 Расчетная оценка по данным модели эталонной атмосферы в слое 30-100 км 3,60 3,00 0,36 0,30

3 Ожидаемая задержка в слое 0-100 км 250,00 234,70 5,61 2,12

Погрешности экспериментальной оценки зенитной задержки складываются из аппаратурных ошибок измерения метеовеличин и наличия невыявленных выбросов в реальном профиле N(0) вследствие слоистого строения нижнего слоя тропосферы.

Аппаратурная ошибка оценки зенитной задержки ст б в значительной мере зависит от качества аэрологических зондов и принятой методики измерений и при наличии предстартовой калибровки по наземным приборам может быть порядка стN прибор = 0,5 №ед. [16]. В работе [17] показано,

что систематические составляющие погрешности измерений по совокупности зондов имеют нулевое значение. Поэтому аппаратурную погрешность оценки зенитной задержки тт можно оценить из

30

соотношения стприбор = 10-6 ^ СТN прибор& , из К°торого следует что стприбор = 1,5 см.

0

Для оценки влияния крупномасштабных неоднородностей вертикального профиля N(0) на измерения тт. воспользуемся результатами работы [18]. В этой работе приводятся параметры неодно-

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

июнь, перечисленные в табл. 4 с указанием среднего числа неоднородностей q в каждом из пяти почти одинаковых по толщине слоев нижней тропосферы до высот 5 км и среднеквадратического уклонения ст^ (дополнительного отклонения относительно среднего градиента индекса преломления внутри неоднородностей). Там же представлены результаты расчетов среднеквадратического отклонения индекса преломления в слое стм = , величины фазовых флуктуаций длины (задержки) трассы стф в каждом километровом слое нижнего слоя тропосферы и суммарные флуктуации в слое тропосферы от 0 до 5 км ст^_5.

Таблица 4

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

нижняя тропосфера до 5 км

№ п/п Параметр Значение

1 Нижняя и верхняя границы слоя тропосферы над уровнем моря, км 0,2-1 1-2 2-3 3-4 4-5

2 Толщина слоя г, м 800 1000 1000 1000 1000

3 СКО градиента уклонения стg, №ед./м 0,58 0,19 0,06 0,009 0,008

4 Среднее число неоднородностей q 5,35 4,6 2,9 2,0 1,9

5 Средняя толщина неоднородности 1е, м 149,5 217,4 344,6 500,0 526,3

6 СКО флуктуаций индекса преломления в слое ст N =ст 21е, ^ед. 86,7 41,3 20,7 4,5 4,2

7 СКО стф фазовых флуктуаций в слое высот, см 2,4 1,5 0,9 0,3 0,2

8 СКО флуктуаций длины трассы (задержки) в слое вьгсот 0-5 км СТо_5 = ^/ст0,2_1 + ст12 + ст2_ 3 + ст32_ 4 + ст4_5 3,00 см

СКО флуктуаций зенитной задержки радиоволн в неоднородностях рассчитывались по формуле

[19]

ynmN 10-12 4

аф = 2V2 '

где r - толщина слоя тропосферы с крупными неоднородностями.

Например, для слоя высот 0,2-1 км флуктуации длины (задержки) трассы за счет неоднородностей составят

2 yi^ral 10-12/е 4П-800• 86,72 • 10"12 • 149,5 п ппп.„ 2 0 ,

<0,2-1 =-272- = 100--- = 0,000563 м2 или Сф^^ = 2,4 см.

В работе [18] указывается, что при толщине крупномасштабных неоднородностей 200 м и менее отечественная система зондирования не измеряет их параметров, влияние таких слоев в конечный результат расчетов TTj не вносится. Однако поскольку эти слои фактически влияют на распространение радиоволн, то возникает дополнительная ошибка в оценках TTj. Как видно из табл. 4, не

регистрируемые системой аэрологического зондирования слои располагаются на высотах до 2 км. Таким образом, дополнительные флуктуации в оценке задержки TTj за счёт неучтённых крупномасштабных слоёв составляют примерно акруп = 2,83 см.

Влияние мелкомасштабных неоднородностей тропосферы на задержку определим по данным [19], где размер неоднородности (внешний масштаб турбулентных неоднородностей) принят равным 1 км, величина СКО флуктуаций индекса преломления aN »(10_3...10~2)N или 0,3-3 N-ед. Выполнив расчеты флуктуаций по формуле (10) для aN » SN-ед., получим, что дополнительный вклад мелкомасштабных неоднородностей во флуктуации тропосферной задержки в слое 0-5 км составляет амелк = 0,5 см и менее, этот вклад учитывается в ходе аэрологического зондирования, если длина

пути зонда внутри неоднородности более 200 м. В противном случае дополнительный вклад мелкомасштабных неоднородностей в флуктуации задержки ст не превышает стмелк = 0,1 см.

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

I 2 2 I 2 2"

Стнеод = "\/сткруп + Стмелк = V2,83 + 0,1 - 2,83 см и менее, а зимой порядка Стнеод - Стмелк = 0,1 см.

Анализ результатов расчетов. Учёт дополнительных флуктуаций за счет тропосферных слоистых и мелкомасштабных неоднородностей и устранения аппаратурной погрешности позволяет дать оценки СКО флуктуаций задержки, более точно характеризующие флуктуации в тропосферном канале распространения радиоволн. Эти значения стТтррв составляют 6,1 см летом и 1,5 см зимой и

могут быть достигнуты при использовании более точной и менее инерционной аппаратуры зондирования.

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

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

I 2 2

стхтмод =устДтт +СТттРРВ .

Результаты расчётов по формуле (11) без учёта постоянного смещения Дтт представлены в табл. 5.

Таблица 5

Результирующие значения СКО оценок задержек, см

№ п/п Вертикальный профиль СКО оценки задержки Сезон

лето зима

1 Экспериментальный Экспериментальное СКО ст , см 5,6 2,1

СКО флуктуаций задержки в тропосферном канале ¡2 2 2 СТтт РРВ = у СТтт + СТнеод — СТприбор 6,1 1,5

2 Экспоненциальный СКО оценки задержки относительно эксперимента СТДтт 6,6 3,2

СКО оценки задержки по модели 1 2 2 СТттмод ~ уСТДтт + СТттРРВ 9,0 3,5

3 Трехэлементный СКО оценки задержки относительно эксперимента стДт 4,7 1,8

СКО °ценки задержки по модели ст^мод =^стДХт + ст2Хтррв 7,7 2,3

Выводы. Представленные выше результаты позволяют сделать следующие выводы.

1. Оценка тропосферной задержки по данным аэрологического зондирования в слое атмосферы 0-30 км в юго-восточной части Охотского моря показывает, что средняя зенитная задержка в слое атмосферы 0-30 км составляет 246 см летом и 232 см зимой, а в слое 0-100 км увеличивается до 250 см летом и 235 см зимой. СКО флуктуаций тропосферной задержки составляют 5,6 см летом и 2,1 м зимой, а в слое 0-100 км могут достигать значений 5,61 и 2,12 см соответственно.

2. СКО фактических флуктуаций тропосферной задержки с учетом того, что аэрологическая аппаратура не регистрирует слои менее 200 м, в слое 0-30 км должно быть увеличено на величину 2,83 см летом и 0,1 см зимой и уменьшено на 1,5 см за счет учета аппаратурных флуктуаций. В итоге СКО составит 6,1 см летом и 1,5 см зимой и может быть получено на безынерционной и точной аппаратуре.

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

экспериментальных данных составила 4,7 см летом и 1,8 см зимой. Среднее смещение при этом для трехэлементной модели составило —11,9 см летом и 0,5 см зимой.

4. Ошибка оценки тропосферной задержки без учёта постоянного смещения при расчетах с помощью моделей с привязкой к приземному значению индекса преломления для экспоненциальной модели равняется 9,0 см летом и 3,5 см зимой, для трехэлементной модели - 7,7 см летом и 2,3 см зимой. Приведённые значения СКО, по всей видимости, являются пределом точности оценки зенитной тропосферной задержки по описываемым моделям высотного профиля индекса преломления с привязкой к приземному значению, которое достигается при исключении аппаратурной погрешности.

Литература

1. Захаров Ф.Н. Сравнение точности оценки времени задержки навигационных сигналов при использовании различных моделей высотного профиля индекса преломления тропосферы / Ф.Н. Захаров, М.В. Крутиков // Доклады Томского государственного университета систем управления и радиоэлектроники. - 2014. - № 2(32). - С. 7-12.

2. Антонович К.М. Использование спутниковых навигационных систем в геодезии: в 2 т. -Т. 1. - М.: ФГУП «Картгеоцентр», 2005. - 334 с.

3. Перов А.И. ГЛОНАСС. Принципы построения и функционирования / А.И. Перов, В.Н. Хари-сов. - М.: Радиотехника, 2010. - 800 с.

4. Саастамойнен Ю. Тропосферная и стратосферная поправки радиослежения ИСЗ / Ю. Саа-стамойнен // Использование искусственных спутников в геодезии / под ред. С. Хенриксена, А. Ман-чини, Б. Човица. - М.: Мир, 1975. - С. 349-356.

5. Hopfield H.S. Two-quartic tropospheric refractivity profile for correcting satellite data / H.S. Hop-field // J. of Geophysical Research. - 1969. - Vol. 74, № 18. - P. 4487-4499.

6. Справочник по радиолокации / под ред. М. Сколника; пер.с англ. под общ. ред. К.Н. Трофимова: в 4 т. - Т. 1. Основы радиолокации / под ред. Я.С Ицхоки. - М.: Сов. радио, 1976. - 456 с.

7. Колосов М.А. Рефракция электромагнитных волн в атмосферах Земли, Венеры и Марса / М.А. Колосов, А.В. Шабельников. - М.: Сов. радио, 1976. - 221 с.

8. Комаров В.С. Справочник. Статистические характеристики полей температуры и влажности в атмосфере Северного полушария. - Ч. 5. - М.: Гидрометеоиздат, 1981. - 87 с.

9. Красюк Н.П. Влияние тропосферы и подстилающей поверхности на работу РЛС / Н.П. Кра-сюк, В.Л. Коблов, В.Н. Красюк. - М.: Радио и связь, 1988. - 216 с.

10. Рекомендация МСЭ-R P.835-4. Эталонные стандарты атмосферы [Электронный ресурс]. -Режим доступа: www.itu.int/dms_pubrec/itu-r/rec/p/R-REC-P.835-4-200503-S! IPDF-R.pdf, свободный (дата обращения: 14.02.2014).

11. Бин Б.Р. Радиометеорология / Б.Р. Бин, Е.Дж. Даттон. - Л.: Гидрометиздат, 1971. - 363 с.

12. Новиков А.В. Модель пространственно-временного поля индекса преломления приземного слоя атмосферы / А.В. Новиков, Ю.П. Акулиничев // Доклады Томского государственного университета систем управления и радиоэлектроники. - 2010. - № 1(21). - С. 36-44.

13. Рекомендация МСЭ-R P.453-10. Индекс рефракции радиоволн: его формула и данные о рефракции [Электронный ресурс]. - Режим доступа: http://www.itu.int/dms_pubrec/itu-r/rec/p/R-REC-P.453-10-201202-I!!PDF-R.pdf, свободный (дата обращения: 14.06.2014).

14. Радиоклиматический тропосферный атлас Тихого океана / под ред. Г.С. Шарыгина. - Томск: Томский гос. ун-т систем управления и радиоэлектроники, 2000. - 171 с.

15. Першин Д.Ю. Сравнительный анализ моделей тропосферной задержки в задаче определения местоположения высокой точности в спутниковых навигационных системах ГЛОНАСС/GPS // Вестник Новосибирского государственного университета. Серия: Информационные технологии. -2009. - Т. 7. № 1. - С. 84-91.

16. Павлов Н.Ф. Аэрология, радиометеорология и техника безопасности. - Л.: Гидрометеоиздат, 1980. - 432 с.

17. Фридзон М.Б. Методология радиозондирования атмосферы и достоверность измерений вертикальных профилей температуры и влажности до высот 35-40 км: дис. ... д-ра техн. наук. - М., 2004. - 323 с.

18. Казаков Л.Я. Неоднородности коэффициента преломления воздуха в тропосфере / Л.Я. Казаков, А.И. Ломакин. - М.: Наука, 1976. - 164 с.

19. Рытов С.М. Введение в статистическую радиофизику. - Ч. II. Случайные поля / С.М. Рытов, Ю.А. Кравцов, В.И. Татарский. - М., 1978. - 464 с.

Захаров Фёдор Николаевич

Аспирант каф. радиотехнических систем (РТС) ТУСУРа

Тел.: (382-2) 41-38-89

Эл. почта: fzakharov89@gmail.com

Госенченко Сергей Григорьевич

Канд. физ.-мат. наук, ст. науч. сотр. НИИ РТС ТУСУРа

Тел.: (382-2) 41-38-89

Эл. почта: rwplab@sibmail.com

Крутиков Михаил Владимирович

Зав. лаб. распространения радиоволн НИИ РТС ТУСУРа

Тел.: (382-2) 41-39-69

Эл. почта: rwplab@sibmail.com

Zakharov F.N., Gosenchenko S.G. Krutikov M.V

Estimation of zenith tropospheric delay time of navigation signals according to meteorological measurements near the navigation receiver at the South-Eastern area of Okhotsk Sea

In this paper corrections to the measured values of the tropospheric delay of a navigational signal are calculated using the refractive index height profiles obtained as a result of upper-air sounding of the atmosphere in the southeastern part of the Sea of Okhotsk, as well as two model profiles. It is shown that the three-element model is more consistent with the experimental profiles than the exponential one. Values of the tropospheric delay RMS for these two models are given.

Keywords: navigation signal delay, troposphere, refraction coefficient, mean-square error.